流行性乙型脑炎垂直分布数学模型.docVIP

流行性乙型脑炎垂直分布数学模型 　【关键词】 流行 　　摘要：为探讨乙型脑炎垂直分布的数学模型，采用对数曲线=a+blgx对开江县不同海拔高度与乙型脑炎年均(1974~2003年)发病率进行数学模拟，得曲线方程：=148.023776-47.920219lgx,R2=0.953133；经直线化方程假设检验F=122.02，Plt;0.01，按α=0.05水平拒绝H0，证明不同海拔高度与乙型脑炎年均发病率的曲线关系成立，理论估计值与实际观察值y分布一致。 　　关键词：乙型脑炎；垂直分布；对数曲线； 拟合优度 　　流行性乙型脑炎(简称乙脑)是一种分布广泛的自然疫源性疾病。迄今，描述乙脑时间、空间、人群分布的文献颇多［1~5］，而记述该病垂直分布的较少［6，7］，未见定量分析乙脑发病与海拔高度的理论模式。本研究选取地居东径107°40~108°05，北纬30°47~31°19，乙脑流行区的川东北低山丘陵地区的开江县。将该县不同海拔高度的乙脑的年均(1974~2003年)发病率资料作对数曲线=a+blgx拟合分析，以探求乙脑垂直分布的规律性。 　　1 资料来源 　　开江县1974~2003年乙脑疫情资料报告，并经个案调查核实，排除地址不明和重复病例，共计903例为本次研究对象。按各村海拔高度进行分组归纳统计，见表1。表1 开江县不同海拔高度与1974~2003年乙脑年均发病率（略） 　　21 求曲线方程［8］ 　　将表1的海拔高度组中值(m)为x，乙脑年均发病率为y,对x作对数(lg)变换，用最小二乘法求得对数曲线方程=148.023776-47.920219lgx。将x值分别代入方程得相应的值。结果与y值符合，见图1。 　　22 曲线拟合优度检验 　　221 直线化方程的假设检验［9］ 　　曲线拟合是利用直线化的原理完成的，证明直线化方程成立，则曲线化关系亦成立。H0:总体回归系数β=0；H1:β≠0；α=0.05。 SS总=Lyy=εy2-(εy)2n=459.4522SS回=L2yy Lxx=437.918888SS剩=SS总-SS回=21.533312 F=SS回/1 SS剩/(n-2)=122.020864令n1=1，n2=8-2=6，查方差齐性检验用F值表示［10］，得：F=122.02gt;F0.01=18.64，按α=0.05水平拒绝H0，可以认为x、y呈直线关系，则=148.023776-47.920219lgx对本例乙脑年均发病率(y)与海拔高度(x)的曲线关系能成立。 　　222 拟合优度用决定系数［11］ 　　剩余均方对Lyy的比例越小，曲线拟合则越优，用决定系数R2表示。代入算式得R2=1-1剩余 Lyy1-ε(y-)2 ε(y-)2=0.953133。R2愈接近1，表明曲线拟合愈优。 　　23 曲线理论估计值区间判定 　　先求出估计标准差(S)=ε(y-)2 n-2=21.533312 8-2=1.89，再用±S或±2S求出各海拔高度乙脑发病率的近似估计值区间，见表2。表2 曲线方程估计结果及区间判定（略）注：* y落在±S区间判“+”，落在±2S区间判“++”，落在gt;或lt;±2S区间判“-”。 　　　3 讨论 　　为了描绘乙脑垂直分布的理论模式，我们应用对数曲线=a+blgx对川东北低山丘陵地区乙脑流行区的开江县不同海拔高度脑年均(1974~2003年)发病率进行了数学模型拟合分析，经假设检验结果表明，该县不同海拔高度与乙脑年均发病率拟合得：=148023776-47.920219lgx的曲线关系成立，决定系数R2=0.953133，接近1，证明对数曲线拟合非常良好。同时，表2表示，875%的y值落在该曲线方程±S的区间范围，12.5%落在±2S的区间，没有超出±2S区间的情况，表明乙脑发病率与海拔高度间呈对数曲线关系。 　　对数曲线的数学性质，曲线递降速度由慢而快，逐渐加速，理论曲线与横座标可相交于零点。本模型外推，当x=1150时，=1.35/万，继续外推，发病率可降至0。从此曲线对该例的拟合及外推结果表明，定量描述了乙脑发病率有随海拔高度的升高而递降的自然规律。乙脑垂直分布的这种特征，可能与当地气温高低、蚊媒数量多少及人群免疫水平等相关条件有着密切的联系。一般而言，海拔越高，气温越低，越不利于蚊媒的繁殖孳生。1987年7月，我们选择海拔800m的山区(灵岩乡)、600m的丘陵(天师镇)和400m的平坝(任市镇)进行蚊密度调查，结果蚊密度山区为96只/人工小时，丘陵为180只/人工小时，平坝为204只/人工小时，与乙脑发病高低分布完全一致；另一方面，平坝地区人口密集，人群集中，流动儿童相对较多，给乙脑疫轩接种工作带来了一定的困难，形成了儿童乙脑免疫薄弱环节。这就确定了乙脑垂直分布呈对数曲线