利用数学模型探讨计划免疫前后麻疹发病季节特征.pdfVIP

利用数学模型探讨计划免疫前后麻疹发病季节特征 王成科，袁卫华，吕勇，李革 四川省蓬溪县卫生防疫站，四川蓬溪（629100 ） E-mail ：Cheng_ke568@163.com 摘 要：本文应用圆形分布构成比法与余弦模型对某地开展计划免疫前后 12 年麻疹发病季 节特征进行分析。求得计划免疫前麻疹的发病高峰时点为4 月中旬，高峰时区为2 月上旬～ 6 月下旬；计划免疫后麻疹发病的高峰时点为3 月下旬，高峰时区为元月下旬～5 月中旬， 显示开展计划免疫后，麻疹发病高峰有所提前，表明控制措施也应提前，通过对求得模型的 拟合，并求得第二谐量三角多项式的决定系数均在 96%以上，表明用该模型拟合该资料是 可行的。 关键词：圆形分布 构成比 余弦模型 麻疹 发病季节特征 疾病的发生是以时、日、月、季、年为观察单位，具有周而复始性，没有真正零点的圆 形变量。圆形分布构成比法是在圆形分布的基础上，为了资料的保密，在此基础上采用构成 比来分析周期性资料。余弦模型分析是周期现象的简单模型，可用于分析角度、生物节律或 时间呈周期性圆变量[1~6] 。本文应用此两种方法对某地开展计划免疫前后（简称计免前后）， 麻疹发病季节变动进行拟合分析，效果良好，叙述如下。 1 材料与方法 本文资料来源于某地26 年疫情资料。将计免前后各年份麻疹发病按月合并，求月均发 病率，见表1，并取对数。为了计算方便，每月按30 天计，取每月15 日为时点。圆形分布 构成比法以各月发病率占全年12 个月发病率总和的百分率（R ）。求平均角和拟合余弦曲线 来分析有无季节变动或变动模式。 表 1 某地开展计免前后麻疹月平均发病率（10-5 ） 1 2 3 4 5 6 7 8 9 10 11 12 合计 计免前 2.3057 2.9082 3.3691 5.6754 5.2169 4.6533 2.6005 1.1328 0.7486 0.8428 1.1971 1.8933 32.544 计免后 2.1128 2.8393 4.1377 3.2708 1.8217 1.0467 0.6470 0.3560 0.2364 0.3325 0.8253 0.1363 17.790 1.1 求平均角 将各月份换算成角度α，求平均角α 2 2 1/2 X=ΣRC , Y=ΣRS , r=(X +Y ) li li 式中R 为各月发病率占全年12 个月总发病率的百分率，C =cosα,S =sinα,r 表示集中趋 li li 势的程度，查《圆分布r 值表》，若有显著性则各月发病率有集中趋势，此时平均角对应的 − 时间，即为疾病的高峰时点，高峰时区为x ±s 。 S （圆标准差）=122.9548(-lg)1/2 平均角ā 的求解有下列三条原则； -1 tg (Y/X)，当X0 时 ɑ= tg-1(Y/X)+180º，当X0 时 -1 tg (Y/X)+360º，当X0,Y0 时 例外情况 -1- 90º 如X=0 ，Y0 ɑ= 270º 如X=0 ，Y=0 不能确定，如X=0 ，Y=0 1.2 拟合余弦曲线 一条余弦曲线可用中值M、振幅A 、峰值位相