食物中毒是指食用了被有毒有害物质污染的食品或者食用了含有毒有害物质的食品后出现的急性、亚急性疾病^[1]。近年来，尽管政府部门、食品生产加工企业、餐饮单位、食品安全研究人员及消费者给予了高度重视和关注，但每年均有重大食物中毒事件发生^[2]。食物中毒的发生严重危害人民群众的身心健康，影响人民群众的正常生活和工作，甚至有可能危及社会的稳定。因此，如何有效地防控食物中毒事件的发生，为人民群众提供良好的食品安全保障，已经成为全社会广泛关注的问题。

研究表明，利用数学统计方法对重大食物中毒事件发生的规律及其原因进行分析和探讨，对于控制和减少食物中毒事件的发生具有指导性意义^[2]。徐君飞和张居作^[3]对2001-2010年中国食源性疾病暴发监测与报告系统中资料进行了统计和分析，金连梅和李群^[4]分析了2004-2007年全国食物中毒事件，陈东军^[5]对2006-2012年全国食物中毒情况进行分析。这些研究虽然对食物中毒事件数据进行了统计并分别得出了食物中毒事件在季节、中毒因素、中毒场所中的病死率，但缺少更深入的数据统计分析。因此，本研究建立在前人研究的基础之上，首次对我国2001-2013年重大食物中毒事件数据进行较为系统的多元数学统计分析，以期深入了解我国重大食物中毒事件发生的特征，为有效地预防和控制食物中毒事件的发生提供科学的参考依据。

本研究数据来源于：中华人民共和国卫生部(http://www. moh. gov. cn)关于2001-2013年全国重大食物中毒事件数据情况的通报^{[6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]}。

季节分解法是时间序列分析预测法中的一种。其中，季节因子序列可用于描述事件的发生在1个年度内依据一定周期性变化。本文参考季节分解方法^{[22, 23, 24]}，运用EViews软件对2001-2013年重大食物中毒起数、死亡人数所有数据进行分析，探究季节因素对全国重大食物中毒事件发生的影响。

变异系数(CV)是标准差相对于平均数的百分数，是衡量资料变异程度的统计量^[25]。本文通过对季节因素、中毒因素和中毒场所变异系数的分析，判断各因素对重大食物中毒事件发生变化情况的影响。

本文运用t值检验法分别对各因素分类下的食物中毒数据进行显著性检验^[26]，分析各中毒因素、中毒场所对重大食物中毒事件发生影响的大小。

据卫生部重大食物中毒事件统计数据显示(图1)：2001-2013年期间，重大食物中毒上报起数最少年份为2002年(共128起)，最多年份为2006年(共596起)；死亡人数最少年份为2013年(共109人)，最多年份为2003年(共323人)。自2006年起，食物中毒的发生起数、死亡人数均呈下降趋势。

图1(Fig.1)

图1 2001-2013年重大食物中毒事件发生起数、死亡人数折线图 Fig.1 The number of occurrence and death of major food poisoning from 2001 to 2013

运用季节分解法对2001-2013年不同季节的重大食物中毒事件发生起数和死亡人数的数据进行作图分析。如图2所示，2001-2013年全国重大食物中毒事件在每年每季均有发生，且呈现先上升后下降的变化趋势，均为第三季度>第二季度>第四季度>第一季度。第三季度是食物中毒事件及死亡人数最多的季度，应给予重点预防和监控。

图2(Fig.2)

图2 2001-2013年每年全国食物中毒发生起数(a)、死亡人数(b)的季节分布 Fig.2 Number of occurrence (a) and death (b) of major food poisoning of each season from 2001 to 2013

为了进一步探究以上数据是否随季节变化而变化，本研究采用季节分解法分析了重大食物中毒事件发生起数和死亡人数的季节因子序列，该数据分析方法可用于描述一个年度内重大食物中毒事件发生的周期性变化规律。如图3所示，2001-2013年期间，每年全国重大食物中毒事件发生起数和死亡人数随季节变化均呈现明显的季节性周期特征。本研究所揭示的这种变化规律可为今后针对不同季节制定相应防控食物中毒事件发生的措施提供理论依据。

图3(Fig.3)

图3 2001-2013年重大食物中毒发生起数(a)、死亡人数(b)季节因子序列分析结果 Fig.3 The seasonal factor sequence of the number of occurrence(a)and death of major food poisoning(b) from 2001 to 2013

食物中毒事件的发生存在着许多不确定因素，如季节、中毒因素及中毒场所等。季节分解法虽然能够揭示全国重大食物中毒事件发生的季节周期性特征，但不能分析和探究不同因素对食物中毒事件发生起数和死亡人数的影响。因此，为了更好地了解在不同因素影响下重大食物中毒事件数量的变化情况，本文分别对各中毒因素、中毒场所与重大食物中毒事件之间的关系进行变异程度分析^[27]，结果如表1、表2所示。

表1(Tab. 1)

表1 食物中毒事件按中毒因素分类发生起数、死亡人数分析Tab. 1 The classified analysis of number of occurrence and death by poisoning factors 年份year 发生起数 number of occurrence 死亡人数 number of death 微生物性microorganism 化学性chemistry 有毒动植物plant and animal toxicants 其它others 微生物性microorganism 化学性chemistry 有毒动植物plant and animal toxicants 其它others 2001 49 90 33 13 20 81 35 10 2002 33 57 27 11 1 99 28 10 2003 97 146 81 55 31 180 94 18 2004 136 138 75 48 34 165 69 14 2005 51 84 65 56 10 106 76 43 2006 265 103 151 77 18 78 85 15 2007 174 89 189 54 5 74 167 12 2008 172 79 125 55 5 57 80 12 2009 118 55 81 17 20 66 93 2 2010 81 40 77 22 16 48 112 8 2011 78 30 53 28 14 57 51 15 2012 56 21 72 25 16 19 99 12 2013 49 19 61 23 1 26 79 3 百分比/%percentage 34.99 24.49 28.06 12.46 7.67 42.43 42.91 6.99 均值mean value 104.50 73.20 83.80 37.20 14.70 83.20 82.20 13.40 标准偏差standard deviation 67.16 41.12 45.81 21.05 10.39 47.61 35.35 10.00 CV/% 64.24 56.22 54.64 56.53 70.70 58.61 43.02 74.68

表1 食物中毒事件按中毒因素分类发生起数、死亡人数分析 Tab. 1 The classified analysis of number of occurrence and death by poisoning factors

表2(Tab. 2)

表2 食物中毒事件按中毒场所分类发生起数、中毒人数分析Tab. 2 The classified analysis of number of occurrence and death by poisoning places 年份year 发生起数 number of occurrence 死亡人数 number of death 集体食堂collective canteens 家庭home 饮食服务restaurant 其他others 集体食堂collective canteens 家庭home 饮食服务restaurant 其他others 2001 62 64 17 42 4 106 5 31 2002 55 54 10 9 6 77 43 12 2003 137 168 34 40 4 275 6 38 2004 133 188 21 55 15 232 0 35 2005 62 123 43 28 6 201 0 28 2006 237 181 86 92 5 162 3 26 2007 142 219 61 84 1 228 4 25 2008 162 147 64 58 4 132 2 16 2009 50 145 51 25 7 154 5 15 2010 37 106 27 50 0 145 3 36 2011 44 86 28 31 4 99 3 31 2012 42 96 12 24 3 128 8 7 2013 37 81 22 12 3 95 1 10 百分比/%percentage 30.90 42.69 12.26 14.16 2.49 81.72 3.33 12.45 均值mean value 92.31 127.54 36.62 42.31 4.77 156.46 6.38 23.85 标准偏差standard deviation 63.19 51.65 23.02 25.29 3.63 60.81 11.24 10.62 CV/% 68.46 40.50 62.88 59.78 76.16 38.87 176.11 44.54

表2 食物中毒事件按中毒场所分类发生起数、中毒人数分析 Tab. 2 The classified analysis of number of occurrence and death by poisoning places

目前，我国食物中毒的主要原因分为微生物性因素、化学性因素、有毒动植物因素等。由表1统计分析可知，化学性因素所引起食物中毒事件死亡人数均值最大，且变异系数小，说明在我国化学性因素引起的食物中毒事件对人民生命财产造成的危害相对较大、更为严重。而微生物性因素，中毒事件发生起数均值、变异系数相对较大，原因可能是微生物存在于食品之中具有难以控制的不确定性与变异性，且微生物生长的繁殖具有一定程度的季节周期性^{[28, 29]}。

由表2可知，家庭场所是食物中毒事件发生最多的场所，也是重大食物中毒事件死亡人数均值最多且变异度较小的场所。集体食堂中毒起数均值虽较高，但其造成的生命危害并不大，且变异度相对较大。因此，家庭区域是更应增强预警、防控的场所。

由2.2分析可知，第三季度是我国食物中毒事件发生的主要时间段。为了预防未来重大食物中毒事件发生，特对第三季度的食物中毒原因、中毒场所数据进行t值检验显著性分析(P<0.05)。

由表3中第三季度食物中毒起数的分析数据可知，通过分析微生物性、有毒动植物、化学性因素、未知因素之间的概率值关系，可得这4种因素对重大食物中毒事件发生起数的影响顺序依次为：微生物 > 化学性 ≈ 有毒动植物 > 未知。第三季度食物中毒死亡人数显著性分析结果可知(表3)：4种因素对重大食物中毒事件死亡人数的影响顺序依次为：化学性 ≈ 有毒动植物 > 微生物 ≈ 未知。以上数据表明，第三季度中，微生物是引起食物中毒的最主要因素，而死亡人数主要是由化学性食物中毒造成的，这一结果与李婷婷的研究^[30]保持一致，说明微生物和化学性食物中毒均为我们公众健康的重大威胁，其中微生物引起食物中毒危害之广，化学因素引起食物中毒危害之深，均应引起足够的重视，并应作为我国未来食物中毒因素的重点监测和防控对象。

表3(Tab. 3)

表3 食物中毒起数、死亡人数中毒因素显著性分析Tab. 3 Significance analysis of poisoning causes about number of occurrence and death 因 素factor 微生物性microorganism 化学性chemistry 有毒动植物plant and animal toxicants 其他others a b a b a b a b 微生物microorganism - - 0.164 <0.001** 0.368 <0.001** 0.004** 0.736 化学性chemistry 0.164 <0.001** - - 0.537 0.956 0.012* <0.001** 有毒动植物plant and animal toxicants 0.368 <0.001** 0.537 0.956 - - 0.003** <0.001** 未知others 0.004* 0.746 0.012* <0.001** 0.003** <0.001** - - 注：a. 相对应因素与其它各因素引起的食物中毒起数显著性分析；b. 相对应因素与其它各因素引起的食物中毒死亡人数显著性分析。 Note: a. significant analysis of the number of occurrence between the corresponding factor and other factors; b. significant analysis of the number of death between the corresponding factor and other factors.

表3 食物中毒起数、死亡人数中毒因素显著性分析 Tab. 3 Significance analysis of poisoning causes about number of occurrence and death

第三季度食物中毒起数中毒场所显著性分析如表4所示。通过分析集体食堂、家庭、饮食服务场所、其他场所间的概率值关系，可得这4种场所中重大食物中毒事件发生起数的显著顺序为：家庭≈集体食堂>饮食服务场所≈其他场所。集体食堂和家庭场所引起的食物中毒起数没有表现出显著性差异。食物中毒死亡人数发生场所显著性分析结果显示仅有集体食堂和饮食服务场所间的概率值(0.627)大于显著水平0.05，其他因素间P值均表现为极显著(P<0.001)。综合考虑概率值关系，4种场所的显著顺序为：家庭>其他场所>集体食堂≈饮食服务场所。家庭是食物中毒事件报告起数和死亡人数最多的场所。

表4(Tab. 4)

表4 食物中毒起数中毒场所显著性分析Tab. 4 Significance analysis of poisoning place about number of occurrence and death 场所place 集体食堂collective canteens 家庭home 饮食服务restaurant 其他others 集体食堂collective canteens - 0.133 0.009** 0.018* 家庭home 0.133 - <0.001** <0.001** 饮食服务场所restaurant 0.009** <0.001** - 0.554 其他others 0.018* <0.001** 0.554 -

表4 食物中毒起数中毒场所显著性分析 Tab. 4 Significance analysis of poisoning place about number of occurrence and death

食物中毒作为一种常见的突发公共卫生事件，存在着许多不确定因素，不同的季节、不同的致病因素、不同的就餐场所及不同的人群均存在一定的差异，但也并非无规律可循。本文通过对2001-2013年全国食物中毒的基本情况分析发现，近十三年全国食物中毒事件发生起数和中毒人数整体上均呈先下降后上升再下降的趋势，尤其是2006年以后全国重大食物中毒发生起数、死亡人数均呈现逐年下降趋势，表明我国食品安全状况已有明显改善，这可能是因为近年来我国颁布了一系列与食品安全相关的法律法规，如2003年《突发公共卫生事件应急条例》、2006年《中华人民共和国农产品质量安全法》、2009年《中华人民共和国食品安全法》等，加大了对生产不合格食品企业、执法不利的相关部门、机构等方面的处罚力度，这些在一定程度上起到了减少食物中毒事件死亡人数的作用^[2]。

本文分析显示2001-2013年全国重大食物中毒事件具有明显的季节性周期特征，每年食物中毒发生起数和死亡人数随季节变化均呈现先增大后减小的变化规律。本研究所揭示的这种变化规律可为今后针对不同季节制定相应防控食物中毒事件发生的措施提供理论依据。第三季度是食物中毒事件发生起数最多、造成危害程度最为严重、变异幅度相对其他季节较小的季度，可能是因气候改变而给食品安全带来的影响。意大利MIRAGLIA等^[28]认为，气候改变对食品安全影响不容乐观，尤其在发展中国家这种影响程度更大。美国TIRADO等^[29]认为气候改变的影响可以涉及整个食品链。第三季度还是各种有毒植物的收获期，缺乏辨别能力的采食者很容易中毒^[30]。

化学因素造成的死亡人数最高，是我国2001-2013年对人民生命财产造成危害较为严重的因素，造成这一现象的原因可能是目前我国仍存在不法分子以化学性物质取代食品添加剂应用于食品中，或者没有任何专业知识及没有接受过任何食品行业有关从业知识培训的人员直接从事食品一线生产等问题。现阶段，造成食物中毒的化学性因素已经引起人们的重视，并提出了一系列的法律法规对其进行控制^{[1, 30]}，而占食品安全事件起数百分比最大的微生物因素却缺乏行之有效的方案，因此我们认为微生物因素应是未来食品安全的防控重点。微生物所引起的食品安全事件具有波及范围广、致病力强等特点，往往导致极其严重的后果和不可估量的经济损失，而微生物自身存在于食品中的不确定性和变异性，也导致其在食品安全事件中难以检测或溯源，进一步加大了其监管的难度。因此，在我国的发展水平越来越接近发达国家的当今，政府及相关部门不仅需要继续保持对引起食物中毒的化学性因素进行监管和控制，而且还应给予微生物因素足够的重视，并配合越来越完善的检测技术和风险评估措施，从而遏制化学性和微生物因素所导致的食物中毒事件的发生，以保障消费者的食品安全和公众健康。

家庭是食物中毒事件报告起数和死亡人数最多的场所，而家庭食物中毒又主要发生在农村地区^[12]。导致这一现象可能是由于农村群众缺少食品安全知识，对自行采食的某些动植物缺乏相应的鉴别能力，同时，农村地区的医疗救治水平有限、交通不便，这更加剧了发生大规模食品安全事故的可能性^[13]。

基于以上分析讨论，本研究对我国未来食物中毒的监管和防控提出以下几点建议：

(1)应保证食品安全相关法律法规的有效实施，构建完善的食品安全风险评估体系，加强相关部门的执法力度，从而保障消费者的食品安全和公众健康；

(2)应区分不同季节进行相应的食品安全监测。其中，第三季度是我国食物中毒事件发生的主要时间段，应给予更高程度的关注；

(3)政府及相关部门不仅需要继续保持对引起食物中毒的化学性因素进行监管和控制，而且还应给予微生物因素足够的重视；

(4)针对农村城镇等食物中毒暴发相对严重的区域，应加强食品安全知识的普及，提高该区域的医疗措施和食物中毒的救治能力。