内蒙古北部荒漠草原沙鼠鼠疫的预报模型
作者：米景川 王成国 吕卫东 刘俊 王兰芳 夏连续
单位：米景川、吕卫东、刘俊、王兰芳、夏连续 (010031 呼和浩特市 内蒙古流行病防治研究所)；王成国 (内蒙古放射卫生防护所)
关键词：
中华预防医学杂志990334 根据1971～1993年内蒙古北部荒漠草原区的天文、气象和动物鼠疫监测资料，应用计算机分析技术，对沙鼠鼠疫的预报模型作初步探讨。
一、材料与方法
1.资料来源：天文资料来源于美国空间环境服务中心。气象资料来源于达茂旗百灵庙草原气象站，动物鼠疫调查资料来源于该地1971～1993年的疫情报告和监测资料。因子包括：x1-太阳黑子数、x2-气温(℃)、x3-地表温度(℃)、x4-降水量(mm)、x5-相对湿度(%)、x6-日照时间(小时)、x7-日照百分率(%)、x8-水热系数、x9-温湿系数、x10-晴雨系数和x11-长爪沙鼠密度(只/公顷)。y-疫情动态值分级法：流行为1，未流行为0，预报结果判别：当预报值yt+1>0.5为流行，当yt+1<0.5为未流行。yt+1值越大，发生流行的可能性越大。
, http://www.100md.com
2. 方法：①多元回归分析和逐步回归分析。②AIC(akaike information crilerition)准则选元法。
二、结果与分析
1. 全因子回归模型：用1971～1992年因子资料与1972～1993年的动物鼠疫流行动态值进行多元回归分析，得到用t时段因子预报t+1时段(未来1年)动物鼠疫的11项因子回归模型。
yt+1(全)=-2.956+0.003 1x1+0.757 4x2+0.013 4x3+0.008 1x4-0.153 1x5-0.003 6x6+0.128 2x7-1.052 6x8+0.695 9x9+51.040 1x10+0.059 1x11；
, http://www.100md.com
R2=0.861 4,F=2.61,P<0.05；拟合符合率=95.46%。
2. 主因子回归模型：①逐步回归选元法：当f=1.6时，筛选得到x1，、x5，、x6，、x10，、 x11的5个主因子；②AIC选元法：选取AIC(Js)=0.75时,最优组合为x1，x3，x5，x6，x11的5因子组合列；经多元回归分析共同的主因子x1，x5，x6，x11，得到了预报未来1年动物鼠疫的主因子回归模型(表1)。
yt+1(主)=-0.055 8+0.002 6x1+0.073 3x5-0.001 1x6+0.040 8x11
, http://www.100md.com
R2=0.742 5,F=5.22,P<0.01；拟合符合率=90.91%。
表1 超前1年预报拟合表 年份
Yt
Yt+1(全)
Yt+1(主)
1971
1
-
-
1972
1
, http://www.100md.com
0.839 6
0.556 5
1973
0
0.027 8
0.360 2
1974
0
0.211 3
0.298 2
1975
0
-0.190 1
, 百拇医药
-0.312 7
1976
0
0.434 4
0.402 4
1977
1
0.470 5
0.421 3
1978
1
1.066 7
, 百拇医药
0.805 3
1979
0
0.018 9
0.370 9
1980
1
1.047 5
0.656 2
1981
1
0.951 2
0.906 3
, 百拇医药
1982
1
1.121 8
1.308 3
1983
1
0.633 5
0.998 0
1984
1
0.617 2
0.856 3
1985
, http://www.100md.com
1
0.998 7
0.814 9
1986
0
0.435 5
0.422 5
1987
0
0.051 2
0.165 6
1998
0
, 百拇医药
-0.193 1
-0.207 5
1989
0
0.270 1
0.348 1
1990
1
0.702 5
0.553 1
1991
1
1.107 8
, http://www.100md.com
0.877 3
1992
0
0.299 4
0.600 0*
1993
1
1.035 3
0.798 6
*为非符合年份
三、讨论
近年来，国内外已注意到太阳黑子、降雨量、水热系数等自然因素的变化可能对鼠疫流行有一定影响，但尚无一致性结论。我们用两种选元方法，对包括天文、气象和生物因子在内的11个因子进行筛选，得到与当年流行关联密切的5个主因子；与未来一年流行关联密切的4个主因子。其中太阳黑子数、相对湿度和长爪沙鼠密度3个主因子相同，其余为复合因子。建立的全因子回归模型和主因子回归模型，可以直接用上一年的因子资料来预测下一年动物鼠疫的流行动态，且从理论上分别达到了95.46%和90.91%的预报准确率。并于1994～1996年运用全因子回归模型对该疫源地的核心地区达茂旗和四子王旗6年次的动物间疫情进行了预报和调查验证，符合率为100%,初步达到了较好的预报效果。
收稿日期：1998-05-29 修回日期：1998-11-03), 百拇医药
单位：米景川、吕卫东、刘俊、王兰芳、夏连续 (010031 呼和浩特市 内蒙古流行病防治研究所)；王成国 (内蒙古放射卫生防护所)
关键词：
中华预防医学杂志990334 根据1971～1993年内蒙古北部荒漠草原区的天文、气象和动物鼠疫监测资料，应用计算机分析技术，对沙鼠鼠疫的预报模型作初步探讨。
一、材料与方法
1.资料来源：天文资料来源于美国空间环境服务中心。气象资料来源于达茂旗百灵庙草原气象站，动物鼠疫调查资料来源于该地1971～1993年的疫情报告和监测资料。因子包括：x1-太阳黑子数、x2-气温(℃)、x3-地表温度(℃)、x4-降水量(mm)、x5-相对湿度(%)、x6-日照时间(小时)、x7-日照百分率(%)、x8-水热系数、x9-温湿系数、x10-晴雨系数和x11-长爪沙鼠密度(只/公顷)。y-疫情动态值分级法：流行为1，未流行为0，预报结果判别：当预报值yt+1>0.5为流行，当yt+1<0.5为未流行。yt+1值越大，发生流行的可能性越大。
, http://www.100md.com
2. 方法：①多元回归分析和逐步回归分析。②AIC(akaike information crilerition)准则选元法。
二、结果与分析
1. 全因子回归模型：用1971～1992年因子资料与1972～1993年的动物鼠疫流行动态值进行多元回归分析，得到用t时段因子预报t+1时段(未来1年)动物鼠疫的11项因子回归模型。
yt+1(全)=-2.956+0.003 1x1+0.757 4x2+0.013 4x3+0.008 1x4-0.153 1x5-0.003 6x6+0.128 2x7-1.052 6x8+0.695 9x9+51.040 1x10+0.059 1x11；
, http://www.100md.com
R2=0.861 4,F=2.61,P<0.05；拟合符合率=95.46%。
2. 主因子回归模型：①逐步回归选元法：当f=1.6时，筛选得到x1，、x5，、x6，、x10，、 x11的5个主因子；②AIC选元法：选取AIC(Js)=0.75时,最优组合为x1，x3，x5，x6，x11的5因子组合列；经多元回归分析共同的主因子x1，x5，x6，x11，得到了预报未来1年动物鼠疫的主因子回归模型(表1)。
yt+1(主)=-0.055 8+0.002 6x1+0.073 3x5-0.001 1x6+0.040 8x11
, http://www.100md.com
R2=0.742 5,F=5.22,P<0.01；拟合符合率=90.91%。
表1 超前1年预报拟合表 年份
Yt
Yt+1(全)
Yt+1(主)
1971
1
-
-
1972
1
, http://www.100md.com
0.839 6
0.556 5
1973
0
0.027 8
0.360 2
1974
0
0.211 3
0.298 2
1975
0
-0.190 1
, 百拇医药
-0.312 7
1976
0
0.434 4
0.402 4
1977
1
0.470 5
0.421 3
1978
1
1.066 7
, 百拇医药
0.805 3
1979
0
0.018 9
0.370 9
1980
1
1.047 5
0.656 2
1981
1
0.951 2
0.906 3
, 百拇医药
1982
1
1.121 8
1.308 3
1983
1
0.633 5
0.998 0
1984
1
0.617 2
0.856 3
1985
, http://www.100md.com
1
0.998 7
0.814 9
1986
0
0.435 5
0.422 5
1987
0
0.051 2
0.165 6
1998
0
, 百拇医药
-0.193 1
-0.207 5
1989
0
0.270 1
0.348 1
1990
1
0.702 5
0.553 1
1991
1
1.107 8
, http://www.100md.com
0.877 3
1992
0
0.299 4
0.600 0*
1993
1
1.035 3
0.798 6
*为非符合年份
三、讨论
近年来，国内外已注意到太阳黑子、降雨量、水热系数等自然因素的变化可能对鼠疫流行有一定影响，但尚无一致性结论。我们用两种选元方法，对包括天文、气象和生物因子在内的11个因子进行筛选，得到与当年流行关联密切的5个主因子；与未来一年流行关联密切的4个主因子。其中太阳黑子数、相对湿度和长爪沙鼠密度3个主因子相同，其余为复合因子。建立的全因子回归模型和主因子回归模型，可以直接用上一年的因子资料来预测下一年动物鼠疫的流行动态，且从理论上分别达到了95.46%和90.91%的预报准确率。并于1994～1996年运用全因子回归模型对该疫源地的核心地区达茂旗和四子王旗6年次的动物间疫情进行了预报和调查验证，符合率为100%,初步达到了较好的预报效果。
收稿日期：1998-05-29 修回日期：1998-11-03), 百拇医药