《基于运气学说的流感样病例预测》由会员分享，可在线阅读，更多相关《基于运气学说的流感样病例预测（7页珍藏版）》请在金锄头文库上搜索。
1、KKME-专业医学搜索引擎 http:/ 作者单位：(中国中医科学院中医药信息研究所,北京 100700)【摘要】 目的 尝试将运气学说应用于流感疾病的预测预报,建立流感预警模型。方法 以天津市 2005-2009 年的流感样病例数据为研究对象,通过分析气象因素(最高温度、平均温度、平均湿度等)以及推算运气学说中的客主加临情况,探索流感和气象因素及运气因子之间的关系。结果 流感的流行与最高温度、最大海平面压力、运气因子有一定的关联关系。结论 构建了一个基于运气学说的流感预警自回归移动平均模型,并且具有较好的流感样病例预测能力。【关键词】 流感样病例;运气学说;客主加临;气象因素;自回归移动平均
2、模型Abstract：Objective To apply Yunqi theory into the influenza epidemic forecast and set up a model of influenza early warning. Methods The influenza-like illness (ILI) data came from Tianjin city of the year 2005-2009 were adopted, meteorological factors (max temperature, average temperature, averag
3、e humidity, etc) were analyzed, and the situation of guest qi adding to host qi of Yunqi theory was calculated. The correlation between influenza and meteorological factors and Yunqi factors were explored. Results The max temperature, max sea level pressure and Yunqi factors had close relationship w
4、ith influenza epidemic. KKME-专业医学搜索引擎 http:/ An ARIMA influenza forecasting model is built based on Yunqi theory, which provides some reference for the dynamic surveillance of influenza.Key words：influenza-like illness;Yunqi theory;guest qi adding to host qi;meteorological factors;ARIMA model在“人与天地相
5、参应“这一整体观念指导下,运气学说着眼于天地自然的运动变化,特别是气候的变化对人以及发生于人身上的疾病的影响。运气学说以时间为变量,通过一系列既定的运算法则,可以推衍出未来某一时段的气候和疾病的发展变化趋势。2009 年上半年甲型 H1N1 流感流行,一些学者运用运气学说对流感的中医病机和证候特征进行了分析1,并对其流行趋势进行了预测。笔者尝试将运气学说应用到流感的预警中,重点讨论客主加临中的“相得“、“相克“与流感的关联关系,探索运气学说对流感疾病的预测预报,采用自回归移动平均模型(ARIMA)建立流感预警模型。1 资料与方法1.1 资料来源1.1.1 流感样病例数据流感样病例(ILI)按照
6、 WHO 和国家流感中心推荐的定义2：发热(腋下体温38 ),伴咳嗽或咽痛之一,缺乏其他实验室确定诊断依据。本研究采用流感样病例百分比(ILI%)作为判断流感的轻重程度的依据,目的是希望通过 ILI 的初期筛查,为流感的早期诊断提供帮助。ILI%资料来自中国国家流感中心公开发布的流感检测数据3,资KKME-专业医学搜索引擎 http:/ 2005 年 10 月 3 日-2009 年 9 月 20 日的数据(共 151 周)。1.1.2 气象数据气象资料来自中国气象科学数据共享服务网4,为同期天津市逐日地面平均气压、平均气温、最小相对湿度、降水量、平均风速等 14 个气象因素。为了与周流感 IL
7、I%对应,气象因素均以周为单位求平均值。1.1.3 运气学说数据处理运气推演主要参考程士德主编的内经讲义5,根据年干支推算公式,推算出 2005-2009 年相应的年干支。应用“天干纪运,地支纪气“理论,推演出每年相应的司天、在泉、主气和客气。司天和在泉见表 1,由于 5 年涉及到的主气和客气数据较多,在此不一一列举。表 1 2005-2009 年流感发病相应的司天、在泉之气表运气学说涉及的数据为文本数据,为了简化对运气因子的数据处理,本文仅讨论运气学说中的“客主加临“情况。客主加临的主要内容为：主气与客气相得;主气与客气不相得,逆;主气与客气不相得,顺。为此,我们对这 3 种“客主加临“情况
8、的处理是：把文本型数据转换成数值型,分别用 0、1、-1 表示。以 2006 年客主加临情况的数据处理为例,结果见表 2。表 2 2006 年客主加临情况的数据处理1.2 统计方法本研究初步选用的气象因素有 14 个,运气因子 1 个,先对这15 个因素与 ILI%进行相关分析,得到显著相关的几个因子。然后采KKME-专业医学搜索引擎 http:/ ARIMA 模型对 ILI%和这 15 个变量进行建模,对模型的残差序列做自相关函数(ACF)图和偏自相关函数(PACF)图进行模型的诊断检验。2 结果2.1 描述性统计分析2005 年 10 月 3 日-2009 年 9 月 20 日期间,天津市
9、最高的 ILI%为25,最低的 ILI%为 9.7。平均气温为 13.14 ,平均海平面压力为 1 018.79 hPa,平均风速为 21.3 km/h,日常温度曲线和海平面压力曲线成正态分布。气象因素变量与 ILI%直观的折线图见图 1(因最大海平面压力数值较高,故单独做图)。图 1 2005-2009 年天津市气象因素变量与 ILI%折线图2.2 相关分析Spearman 相关分析显示,ILI%与运气因子的相关系数最大,而与最高温度、最大海平面压力、最大风速这 3 个气象因素也相关,见表 3。其中 ILI%与最高温度成正相关,即最高温度相对升高时,流感样病例数相对增多。而 ILI%与最大海
10、平面压力、最大风速成负相关,即最大海平面压力相对升高、最大风速相对增大时,流感样病例数相对减少。表 3 ILI%与气象因素变量的 Spearman 相关分析2.3 时间序列分析2.3.1 时间序列模型的识别时间序列模型的建立一般可分为模型识别、模型拟合、控制和预测三大步骤。首先可以使用 ACF 和 PACF 图进行初步识别。KKME-专业医学搜索引擎 http:/ ACF 和 PACF 见图 2。图 2 模型残差 ACF 和 PACF 图ACF 和 PACF 图提示绝大部分相关函数都在可信限范围内,且 Box-Ljung 统计量无显著性,则可认为残差是一个“白噪声“序列,提示所估计的模型已包含
11、了原始时间序列的所有季节性及周期性变动,可接受所估计的模型。故构建 ARIMA(1,0,2)(1,1,2)12 的 ARIMA模型。结果显示气象因素中最高温度、最大海平面压力与 ILI%有关,运气因子与 ILI%也相关,结果见表 4。可以看出,当最高温度升高 1 、最大海平面压力降低 1 hPa,将分别导致 ILI%升高 0.46%、0.61%。表4 ARIMA(1,0,2)(1,1,2)12 模型参数表2.3.2 时间序列模型的预测分析用模型 ARIMA(1,0,2)(1,1,2)12 对 2005 年 10 月 3 日-2009年 9 月 20 日天津市 ILI 各月报告数进行回代预测(组
12、内回代),对 2009年 9 月 21 日-11 月 30 日的 ILI%进行组外回代,结果显示预测值与实际值较吻合,预测的动态趋势与实际情况基本一致,可见模型拟合的效果很好。拟合及预测的具体结果见图 3,预测值见表 5。图 3 2005 年10 月 3 日-2009 年 9 月 20 日天津市 ILI%拟合及预测结果表 5 2009年 9 月 21 日-11 月 30 日天津市 ILI%预测值本研究用 ARIMA 模型预测 2009 年 9 月 21 日-11 月 30 日的 ILI,从预测的数据对比来看,总的趋势是一致的,误差不大,均在置信区间内,可见预测效果良好。仔细观察表 5 数据,实
13、际数据要比预测数据稍大一些,说明天津市 2009 年下半年 ILI%增长的势头比以往要强劲。由于原时间序列只有 5 年的范围,对预测精度是个不利的因素,随KKME-专业医学搜索引擎 http:/ 讨论目前,美国疾病控制与预防中心(CDC)以每年 10 月的平均周病毒分离数作为基准值来判定流感流行季节和非流行季节6。而一些研究表明,ILI%比病毒分离结果更具有早期预警意义7,ILI%可以作为预测流感流行的指标8。因此,本研究采用 ILI%作为判断流感流行程度的依据。可通过 ILI 的初期筛查,为流感的早期诊断提供帮助。基于天津市 2005-2009 年的 ILI%数据,本研究构建了基于气象因素的
14、流感样病例 ARIMA 预测模型,时间序列预测法的自变量是气象因素和“客主加临“的运气因子,因变量是待预测的 ILI%,最后给出了基于(1,0,2)(1,1,2)12 模型的预测,得到了气象因素如最高温度、最大海平面压力等与 ILI%的相关关系,从而为增强流感监测水平及其动态变化规律提供了参考。为了提高应对流感大流行的能力,本研究将运气学说引入到流感样病例 ARIMA 预测模型中,重点解决主气和客气这 2 个运气因子对流感发病的影响。结果表明,客主加临运气因子与 ILI%有着明显的相关关系,研究结果有助于丰富和发展中医运气学说理论。今后考虑在模型中加入主运、客运、胜复等多个因素,逐步深化预警模
15、型,希望该模型能服务于公共卫生平台。【参考文献】1 顾植山.从中医五运六气理论谈对当前甲型 H1N1 流感KKME-专业医学搜索引擎 http:/ 全国流感/人禽流感监测实施方案(2005-2010)R.中国疾病预防控制中心,2005.3 中国国家流感中心EB/OL.2010-04-17.http:/ chn/down/?typied=20&page=9.4 国家气象信息中心EB/OL.2010-04-15.http:/ index.jsp.5 程士德.内经讲义M.上海：上海科学技术出版社,1984：191-203.6 Simonsen L, Fukuda K, Schonberger LB. The impact of influenza epidemics on hospitalizationsJ. JID,2001,181：831-837.7 张顺祥,谷利妞,程小雯,等.季节性流感监测中病毒分子变异分析的意义J.疾病控制杂志,2009,13(1)：23-28.8 Yang L, Wong CM, Lau EH, et
