精神分裂症是一种常见的严重精神障碍疾病，每年所造成疾病负担在所有疾病中排名前10位，2017年全球范围内约有2 000万精神分裂症患者^[1]，因病伤残导致的寿命损失年超过1 270万人年^[2]。2018年中国预计有840万的精神分裂症患者，其加权患病率约为0.6%^[3]。宁波市登记在册的精神分裂症患者超过2.6万，其造成的潜在减寿年数达3.02万人年^[4]。精神分裂症是一种基因和环境因素共同导致的疾病，基因因素和环境因素在精神分裂症的发病中起着重要作用^[5]。空间相关分析和时空扫描分析方法是研究疾病时间和空间分布规律的统计方法，在传染病和部分肿瘤的发生中得到广泛应用^[6-8]，但在精神分裂症中研究较少^[9]。本研究以镇（街道、乡）为单位，采用空间相关分析和时空扫描分析方法，分析2018-2022年宁波市精神分裂症的报告发病情况，为精神卫生资源合理配置和精神分裂症综合防治提供科学依据。

资料与方法

1. 资料来源：2018-2022年宁波市精神分裂症报告发病资料来源于宁波市精神卫生信息管理系统。内容包括患者性别、年龄、确诊日期、现居住地行政代码等。各镇（街道、乡）人口学资料来源于宁波市统计局，电子矢量地图来源于宁波市测绘局。本研究通过宁波市康宁医院伦理委员会审查（批准文号：NBKNYY-2023-LC-22）。

2. 相关指标及定义：①全局空间自相关分析：以全局Moran's I系数反映发病在整个区域的聚集程度，其取值范围为-1~1，当取值为正数时代表空间呈正相关，有相同的变化趋势；当取值为负数时代表空间呈负相关，有不同的变化趋势。取值绝对值越接近于1时代表空间聚集性越强，取值为0时代表空间上不存在空间聚集性，呈随机分布。全局空间自相关分析用Z检验进行显著性评价。本研究中Moran's I系数随机化置换均选择999次。②局部空间自相关性分析：主要反映相邻区域间的空间相关程度，常用局部Moran's I系数进行检验。局部空间自相关性可分为高-高、低-低、低-高、高-低4种聚集模式，高-高聚集模式区域即为热点区域。③时空扫描分析：从时间和空间两个维度探究疾病的分布情况。统计量为对数似然比（LLR）和RR，可以通过扫描窗口内外的实际病例数、理论病例数计算。当LLR差异有统计学意义时（P < 0.05），可认为该地区存在聚集。

3. 统计学分析：使用SPSS 22.0软件进行描述性分析和趋势χ²检验。定量资料呈偏态分布时用M（Q₁，Q₃）表示，计数资料用例数及构成比（%）描述。使用GeoDa 1.14.0软件进行空间自相关分析。使用SaTScan 9.4软件进行时空扫描分析，以离散泊松分布模型为基础进行分析，扫描时间范围为2018-2022年，最大空间聚集范围设定为总人口的30%。唐咸艳等^[10]研究发现20%~30%比例的扫描窗口为最适宜的空间扫描窗口，本研究发现且设定为30%时，LLR最大。使用ArcGIS 10.3软件实现分析结果可视化。双侧检验，检验水准α=0.05。

结果

1. 流行特征：2018-2022年宁波市共报告精神分裂症4 133例，年均报告发病率为9.76/10万。

（1）时间分布：宁波市精神分裂症年报告发病例数和年报告发病率均呈逐年下降趋势，年报告发病率差异有统计学意义（趋势χ²=29 008.00，P < 0.001）。见图 1。

（2）人群分布：2018-2022年宁波市精神分裂症报告发病在各个年龄段均有分布，报告病例最小年龄为7岁，最大为98岁，年龄M（Q₁，Q₃）为53（39，65）岁。男性报告1 745例，女性报告2 388例，男女性别比为0.73∶1。见表 1。

（3）地区分布：2018-2022年宁波市精神分裂症年平均报告发病率最高的5个镇（街道、乡）依次为海曙区江厦街道（29.83/10万）、鄞州区福明街道（29.19/10万）、海曙区龙观乡（28.80/10万）、鄞州区横溪镇（27.34/10万）、海曙区望春街道（24.42/10万）。见图 2。

2. 空间自相关分析：

（1）全局空间自相关分析：2018、2019和2020年宁波市精神分裂症报告发病率的全局Moran's I值分别为0.270、0.161和0.165（均P < 0.05）。2021和2022年精神分裂症报告发病率的Moran's I值差异无统计学意义（均P > 0.05）。见表 2。

（2）局部空间自相关分析：2018年精神分裂症报告发病率的高-高聚集区有12个镇（街道、乡），为宁海县大佳何镇、宁海县梅林街道、宁海县茶院乡、宁海县桥头胡街道、鄞州区云龙镇、鄞州区塘溪镇、鄞州区横溪镇、江北区庄桥街道、江北区洪塘街道、江北区慈城镇、海曙区高桥镇、海曙区望春街道；2019年高-高聚集区有6个镇（街道、乡），为海曙区章水镇、海曙区鄞江镇、海曙区鼓楼街道、鄞州区东柳街道、鄞州区邱隘镇、象山县墙头镇；2020年高-高聚集区有10个镇（街道、乡），为海曙区望春街道、海曙区古林镇、海曙区高桥镇、鄞州区云龙镇、鄞州区东钱湖镇、江北区庄桥街道、江北区甬江街道、奉化区裘村镇、奉化区莼湖镇、镇海区蛟川街道。见图 3。

3.时空聚集性分析：时空聚集区聚集时间为2018-2019年，集群中心为海曙区洞桥镇，覆盖范围半径为25.65 km，包含附近42个镇（街道、乡）。该聚集区域实际病例数为972例，期望病例数为456.64例，RR值为2.46，LLR值为256.89，P < 0.001。无其他次级时空聚集区。见图 4。

讨论

宁波市精神分裂症年报告发病例数和年报告发病率均呈逐年下降趋势，与全球精神分裂症发病趋势一致^[2]。本研究中精神分裂症年均报告发病率为9.76/10万，低于中国1990年的精神分裂症发病率（20.08/10万）和2017年发病率（19.81/10万），也低于全球1990年的精神分裂症发病率（14.98/10万）和2017年发病率（14.39/10万）^[2]。宁波市精神分裂症在各个年龄段均有发病，其中发病人数最多为51~60岁，与He等^[2]研究结果不一致，其发现全球范围内精神分裂症在男女性中均为20~29岁组发病人数最多。从性别分布来看，宁波地区男性精神分裂症报告发病数低于女性，该结果与全球数据相反^[2]。

全局空间自相关分析表明，地理位置不同可能和精神分裂症报告发病数相关，本研究结果显示，2018、2019和2020年精神分裂症报告发病分布存在空间聚集性特征，而2021和2022年则无空间聚集性。宁波市精神卫生信息管理系统2018年正式启用^[4]，专科医院医生直接通过系统进行个案登记，提高报卡的及时性和准确性^[11]，且信息系统会根据诊断结果直接提示医生报卡，减少漏报、错报的发生。2021年后无空间聚集性，可能因为新型冠状病毒感染疫情的发展导致精神疾病门诊量出现减少^[12-13]。

时空扫描分析同时考虑到精神分裂症的发病时间和发病地区两个因素，较单纯的时间或空间聚集性分析，更易发现疾病分布的规律^[14]。本研究发现宁波市精神分裂症报告发病率在时间和空间上均存在明显的聚集性，聚集时间为2018-2019年，集群中心为海曙区洞桥镇及周边的42个镇（街道、乡），可作为精神分裂症的重点防控区域。

本研究存在局限性。首先，精神分裂症发病数据未能对年龄进行标化；其次，本研究采用报告发病率进行分析，由于漏报病例的存在，发病率较实际情况偏低；最后，精神分裂症发病是由家族史、遗传因素、社会心理因素等共同作用^[15-16]，本研究中只收集到宁波市2018-2022年各地区的发病数据，未对其他影响因素进行分析，后续可以进一步分析标化社会经济和人口特征后的报告发病率。

综上所述，本研究采用空间相关分析和时空扫描分析的方法在镇（街道、乡）小尺度上探索精神分裂症报告发病率的时空聚集区域，可为合理分配精神卫生资源提供依据，也为后续进一步开展流行病学研究、探索精神分裂症发病与环境因素之间的关联奠定基础。

利益冲突  所有作者声明无利益冲突

作者贡献声明  李链：整理分析数据、论文撰写；洪航：研究设计/指导、论文审阅；周东升：行政/技术支持；方挺：行政/技术支持、指导；杨红英：论文审阅；边国林、许国章：数据获取、研究设计/指导、论文审阅