河北省长爪沙鼠鼠疫疫源地属于内蒙古高原长爪沙鼠鼠疫自然疫源地的一部分，自2000年以来共发生过两次鼠间鼠疫，分别为2002年和2005年^[1]，之后进入了隐伏期，该疫源地主要宿主为长爪沙鼠，其在本地鼠疫发生与维持过程中起着关键作用^[2]，通常长爪沙鼠鼠密度升高，鼠疫发生的可能性也增加。长爪沙鼠是荒漠草原的代表性优势种动物，但也可分布于农业区，而且不冬眠，一年四季可见其活动。本研究利用2001－2013年河北省长爪沙鼠动物间鼠疫监测数据分析长爪沙鼠在不同年份、生境和季节的分布，并采用多重对应分析探讨它们之间的相互关系，为鼠疫防治提供科学依据。

1. 资料来源：来源于河北省鼠疫防治所2001－2013年河北省长爪沙鼠国家级和省级两个疫源地的动物间鼠疫监测数据。本研究的主要宿主是使用日夹法捕获的长爪沙鼠，分析指标为鼠密度，鼠密度＝捕鼠数（只）/样方面积（公顷）。国家《鼠疫控制及其考核原则与方法》要求长爪沙鼠密度控制标准为≤0.20只/公顷，据此本研究以年为观察单位定义鼠密度≥1.00只/公顷为高危年组、0.20～1.00只/公顷为预警年组、≤0.20只/公顷为标准年组。本研究主要分析的因素有监测生境，分为草地、耕地、荒地、山坡；监测月份为4、5、6、7和10－11月。

2. 统计学分析：用Excel 2010软件进行数据录入和整理，使用SPSS 19.0 软件进行统计学分析。采用非参数Kruskal-Wallis H秩和检验按年份组比较长爪沙鼠生境、月份分布的差异，并对这些变量进行多重对应分析，比较的指标用M、极大值（Max）和极小值（Min）表示，P＜0.05表示差异有统计学意义；多重对应分析中＜0的数据，将每个数据加一常数0.10进行分析^[3]。

1. 一般情况：2001－2013年河北省长爪沙鼠疫源地长爪沙鼠密度有2个高峰，分别为2002年和2004年，2005年之后长爪沙鼠密度直线下降。2001－2005年的鼠密度均＞1.00只/公顷属于高危年组，2005年后明显下降，但2008和2013年的鼠密度均＞0.20只/公顷，2006－2007年和2009－2012年鼠密度＜0.20只/公顷，见表 1。

表 1 2001－2013年长爪沙鼠捕鼠数和密度

表选项

2. 长爪沙鼠生境分布：分别比较高危年、预警年和标准年不同生境的鼠密度分布情况，结果显示，高危年长爪沙鼠密度在不同生境分布不同（H＝13.997，P＝0.003），耕地鼠密度M最高，但荒地的密度Max最高，为20.50只/公顷；两两比较，耕地和荒地（P＝0.012）、耕地和山坡鼠密度差异有统计学意义（P＝0.009）。预警年不同生境长爪沙鼠密度分布不同，两两比较差异无统计学意义（H＝7.899，P＝0.048），其中耕地的鼠密度M最高。标准年不同生境长爪沙鼠密度分布不同，差异有统计学意义（H＝12.017，P＝0.007），两两比较，耕地的鼠密度M高于荒地（P＝0.012），其他生境鼠密度差异无统计学意义，见表 2。

表 2 3组长爪沙鼠不同生境的密度（只/公顷）

表选项

3. 长爪沙鼠月份分布：Kruskal-Wallis H秩和检验显示，高危年鼠密度在不同月份分布不同，差异有统计学意义（H＝12.898，P＝0.012），鼠密度M＞1只/公顷的月份有3个，密度最高值为4月份的20.50只/公顷，6月份和10－11月份次之。预警年不同月份鼠密度分布差异无统计学意义。标准年不同月份鼠密度分布不同，差异有统计学意义（H＝10.599，P＝0.031），5个月份的鼠密度M均小于国家标准，见表 3。

表 3 3组长爪沙鼠不同月份的密度（只/公顷）

表选项

4. 长爪沙鼠密度多重对应分析：按照年份组、生境和月份3个变量对鼠密度进行多重对应分析，前两个变量的累积贡献率为84.1%（表 4）。对应分析图可见，鼠密度高危年相关因素有荒地和4、6月份；而预警年与耕地和7月份有关联；标准年与5月份有关联（图 1）。

表 4 不同年份组、生境、月份的长爪沙鼠密度多重对应分析（只/公顷）

表选项

图 1 不同年份组、月份和生境长爪沙鼠密度多重对应分析

图选项

本研究分析了2001－2013年河北省长爪沙鼠疫源地长爪沙鼠密度，结果显示2001－2005年连续5年长爪沙鼠密度均＞1.00只/公顷，与本地区2002年11月份和2005年4月份发生的长爪沙鼠鼠间鼠疫在时间上相吻合^[1]。提示鼠密度升高容易引起动物鼠疫发生流行，与已有研究结果一致^[4-5]。这5年中长爪沙鼠密度的两个高峰（2002年和2004年），推测鼠密度出现峰值之后的1年动物鼠疫很可能发生，即长爪沙鼠密度远高于上年时要警惕鼠疫的发生。鼠疫发生后通常会积极采取措施，控制传染源、切断传播途径、大量消灭长爪沙鼠，鼠密度降低。本研究结果也显示，2005－2013年长爪沙鼠密度均＜1.00只/公顷。但仍有2008年和2013年长爪沙鼠密度在0.20～1.00 只/公顷，大于《鼠疫控制及其考核原则与方法》中长爪沙鼠控制标准（0.20 只/公顷），这期间虽然没有动物鼠疫流行但仍需密切观察作为预警提示。

随着统计学的发展，越来越多的方法被用于研究鼠疫^[6-8]，本研究尝试使用多重对应分析来研究长爪沙鼠鼠密度年份、生境和季节间的相关关系。对应分析又叫相应分析，是由理查森和库德在1933年提出，之后经法国数学家Jean-Paul Benzecri和日本数学家Hayashi Chikio 改进而发展的一种多元统计方法^[9]。对应分析能综合考虑样品之间、变量之间及它们相互之间的关系，在同一因子平面上对变量和样品一起进行分类，从而揭示多研究变量和样品间的内在联系^[10]。简单对应分析的行列分类变量只有1个，当研究的分类变量＞2时采用多重对应分析。多重对应分析可以将多个分类变量的对应关系通过对应分析图显示出来，从而直观地展示不同分类变量间的相关关系^[11]。

本研究从多重对应分析图直观发现高危年长爪沙鼠鼠密度与荒地和4月份、6月份有较强关联。长爪沙鼠喜居于干旱沙质土壤地区，因而，常见于荒漠草原，但也分布于干草原和农业地区。在干草原的长爪沙鼠除非鼠疫大流行，一般情况下，仅有零星的分布，在撂荒地上往往能形成较高的密度。在农业地区主要栖居于田埂、水渠垄背和人工林边的荒地^[12]。本研究高危年生境单因素分析也发现每公顷鼠密度最高的生境为荒地。同样高危年月份单因素分析结果也显示长爪沙鼠鼠密度最高值4月份和6月份较高。长爪沙鼠繁殖高峰在2－5月份，这个时期的繁殖量可作为夏季宿主数量增加的预警信息。4月份天气变暖，幼鼠出洞，鼠密度增加。繁殖高峰过后，幼鼠增多，再加上长爪沙鼠喜食植物，6月份鼠密度也随之增加。由此可见对于长爪沙鼠除平均鼠密度外，密度极值也应引起重视。多重对应分析结果还显示预警年长爪沙鼠鼠密度相关因素为耕地和7月份，单因素分析的结果发现预警年不仅耕地的长爪沙鼠密度Max最高，而且M也最高，而且在高危年和标准年的M也均较高。在耕地生境中，由于受耕作措施的影响，耕地中地貌整齐划一，长爪沙鼠生活环境更稳定，加之耕地为长爪沙鼠提供了良好的觅食、储粮和隐蔽条件，在相对优越的环境条件下，生存竞争压力较轻，因此长爪沙鼠密度相对较高^[13-15]。虽然对应分析预警年与7月份出现了聚集，但预警组月份单因素的分析没有发现月份分布的差异，可能与样本量少有关，因为本研究的预警年只有2008和2013两个年份；另外，对应分析这种统计方法也容易受极端值的影响^[16]，本研究多重对应分析的列联表显示了耕地组的7月份鼠密度值最小。因此进一步研究预警年长爪沙鼠密度的分布规律需收集更多年份的数据。

综上所述，河北省长爪沙鼠鼠密度＞1.00只/公顷时与研究地区同期动物鼠疫发生情况一致，提示动物鼠疫发生流行的危险性增加。高危年和预警年的生境和月份关联的多重对应分析结果提示，今后监测工作应重点关注耕地和荒地及4和6月份长爪沙鼠密度的变化，防止动物鼠疫的发生和扩散，避免波及人间。