呂 游，柯資能，顧植山，王昌忠，孫 明

（1.中國科學(xué)技術(shù)大學(xué) 科技史與科技考古系，安徽 合肥 230026；2.江陰致和堂中醫(yī)藥研究所，江蘇 江陰 214400；3.山東理工大學(xué) 教務(wù)處，山東 淄博 255049）

0 引言

2014年以來，埃博拉病毒肆掠西非各國，有研究表明，天氣和氣候的變化與人類健康有密切的關(guān)系.［1］筆者檢查 Wolfram 氣象數(shù)據(jù)官方網(wǎng)站［2］提供疫區(qū)氣象資料發(fā)現(xiàn)，疫區(qū)中心幾內(nèi)亞自2011年深秋以來，日平均氣溫嚴(yán)重失常，從2012年初開始，高達(dá)95.7%的日平均氣溫低于相應(yīng)的日平均氣溫，同樣情況也出現(xiàn)在其他埃博拉疫區(qū)，如圖1～8所示，但程度輕微得多，而其他幾個(gè)隨機(jī)對(duì)照組（埃塞俄比亞、馬里、尼日爾、佛羅里達(dá)州、德克薩斯州）則沒有這個(gè)現(xiàn)象.在全球變暖這一背景下，幾內(nèi)亞的日平均氣溫長期低于同日的十年平均值，是否與埃博拉病毒流行相關(guān)？筆者試圖采用統(tǒng)計(jì)學(xué)的方法探討氣溫與埃博拉病例數(shù)的關(guān)系.

1 資料來源和方法

1.1 資料來源

2014／9／7 －2015／1／1西非六地（幾內(nèi)亞［科納克日里］、塞拉利昂［佛里敦機(jī)場(chǎng)］、尼日利亞［阿布賈］、埃塞俄比亞［默特哈拉］、馬里［巴馬科］、尼日爾［尼亞美機(jī)場(chǎng)］和美國佛羅里達(dá)州［波卡拉頓機(jī)場(chǎng)］、德克薩斯州［休斯頓機(jī)場(chǎng)］）的平均氣溫通過 Wolfram氣象數(shù)據(jù)官方網(wǎng)站獲得；2014／8／4－2014／10／5幾內(nèi)亞、尼日利亞和塞拉利昂埃博拉病例資料取自中國疾病預(yù)防中心發(fā)布的埃博拉疫情通報(bào).

1.2 數(shù)據(jù)處理

1.2.1 數(shù)據(jù)修復(fù)

為了充分利用氣象數(shù)據(jù)和保持?jǐn)?shù)據(jù)的連續(xù)性，本研究使用了對(duì)氣象數(shù)據(jù)求滑動(dòng)平均值的方法，即參考了某一時(shí)間點(diǎn)前15日的平均氣溫.對(duì)這十年的平均氣溫進(jìn)行滑動(dòng)平均處理，修復(fù)了數(shù)據(jù)缺失的問題.

1.2.2 積溫差值計(jì)算

利用2014／9／7－2005／1／1的平均氣溫，將相同日月的平均氣溫分別相加并求其平均值，得到逐日平均氣溫b1，b2，b3，…，bi，…，b365（1≤i≤365），結(jié)合表1埃博拉病例數(shù)的日期，計(jì)算t2＿k＝∑365i＝1bi＿k；將2014／9／7－2005／1／1的平均氣溫，以年為單位，進(jìn)行累 加， 得 到t1，t2，t3，…，t9，t10＿k， 其 中t10＿k為2014／8／4，2014／8／6，2014／8／11，2014／8／16，2014／8／18，2014／8／20，2014／8／26，2014／8／31，2014／9／5，2014／9／6，2014／9／7這12個(gè)時(shí)間點(diǎn)的累積平均溫度.具體公式如下所示：

表1 幾內(nèi)亞、尼日利亞和塞拉利昂埃博拉病例數(shù)（存量＋新發(fā)）Tab.1 Number of Ebola Cases in Guinea and Nigeria（original＋new）

1.2.3 日平均氣溫差值小于零的頻度計(jì)算

利用幾內(nèi)亞、尼日利亞、埃塞俄比亞、馬里、尼日爾、塞拉利昂、佛羅里達(dá)州、德克薩斯州7個(gè)地區(qū)從2011／1／1到2014／9／7日平均氣溫ai＿j分別減去同期的逐日平均氣溫bi＿g得到日平均氣溫差值T3＿g（g為8個(gè)地區(qū)），并計(jì)算小于零的天數(shù)頻度，得出幾內(nèi)亞的頻度為95.7%，尼日利亞的頻度為64.79%，塞拉利昂的頻度為82.12%，埃塞俄比亞的頻度為52.12%，馬里的頻度為51.98%，尼日爾的頻度為50.82%，佛羅里達(dá)州的頻度為44.92%，德克薩斯州的頻度為42.89%，并根據(jù)T3＿g繪制出差值圖（如圖1～8所示）.

圖1 幾內(nèi)亞日平均氣溫差值Fig.1Difference of Guinea daily temperature

圖2 尼日利亞日平均氣溫差值Fig.2Difference of Nigeria daily temperature

圖3 埃塞俄比亞日平均氣溫差值Fig.3Difference of Ethiopia daily temperature

圖4 馬里日平均氣溫差值Fig.4Difference of Mali daily temperature

圖5 尼日爾日平均氣溫差值Fig.5Difference of Niger daily temperature

圖6 佛羅里達(dá)州日平均氣溫差值Fig.6Difference of Florida daily temperature

圖7 德克薩斯州日平均氣溫差值Fig.7Difference of Texas daily temperature

圖8 塞拉利昂日平均氣溫差值Fig.8Difference of Sierra Leone daily temperature

1.3 Pearson相關(guān)性分析和一元線性回歸分析［4］

利用1.2.2中所求的積溫差值T3＿k和對(duì)應(yīng)時(shí)間點(diǎn)的疫區(qū)三個(gè)地方（幾內(nèi)亞、尼日利亞、塞拉利昂）埃博拉病例數(shù)，運(yùn)用SPSS 19.0作Pearson相關(guān)性分析和一元線性回歸分析.

2 結(jié)果

2.1 日平均氣溫差值小于零的天數(shù)和頻度結(jié)果

從圖1～8、頻度計(jì)算結(jié)果可以看出，幾內(nèi)亞的日平均氣溫從2011年深秋、塞拉利昂從2012年、尼日利亞從2013年開始長期持續(xù)低于該地區(qū)的逐日平均氣溫，對(duì)照組地區(qū)的日平均氣溫在逐日平均氣溫附近浮動(dòng).對(duì)比西非六個(gè)地方的埃博拉病例數(shù)看，截至2014年9月7日，幾內(nèi)亞最多時(shí)達(dá)到861人，尼日利亞一直保持在21人左右，塞拉利昂最多時(shí)達(dá)到了1412，其他地方暫未發(fā)現(xiàn)疫情.

2.2 Pearson相關(guān)性分析結(jié)果

對(duì)幾內(nèi)亞的積溫差值與發(fā)病人數(shù)作Pearson相關(guān)性分析，結(jié)果如表2所示.從表2可以看出，積溫差值與病例數(shù)之間存在顯著相關(guān)性，特別是幾內(nèi)亞p＝0.000，R＝0.974，塞拉利昂次之.

表2 積溫差值與病例數(shù)Pearson相關(guān)性分析結(jié)果Tab.2 Result of Pearson correlation between Difference of Accumulated Temperature and Cases

2.3 一元線性回歸分析結(jié)果

2.3.1 一元線性回歸方程的建立

利用2014／9／7－2005／1／1幾內(nèi)亞和尼日利亞的積溫差值和病例數(shù)，作兩者的一元線性回歸分析，推導(dǎo)出回歸方程，結(jié)果為：積溫差值和病例數(shù)存在一元線性關(guān)系，模型有統(tǒng)計(jì)學(xué)意義（p＜0.05或p＜0.01），所得到的回歸直線方程是有效的，如表3.

表3 積溫差值與發(fā)病人數(shù)回歸方程Tab.3 Regression Equation of Difference of Accumulated Temperature and Cases

2.3.2 一元線性回歸方程預(yù)測(cè)

以幾內(nèi)亞為例，根據(jù)1.2.2所推導(dǎo)出來的方程，結(jié)合中國疾病預(yù)防中心發(fā)布的最新埃博拉疫情通報(bào)（截至2014年10月5日）的埃博拉病例數(shù)，結(jié)合公式（1）、（2）算出幾內(nèi)亞對(duì)應(yīng)時(shí)間點(diǎn)的積溫差值，作一元線性回歸預(yù)測(cè)，結(jié)果圖9所示.

圖9 幾內(nèi)亞埃博拉病例數(shù)預(yù)測(cè)圖Fig.9Predictive cases of Guinea

從圖9可以看出，幾內(nèi)亞實(shí)際病例數(shù)和預(yù)測(cè)病例數(shù)的變化趨勢(shì)一致，因此該回歸方程可作為計(jì)算總積溫差值恢復(fù)到一定閥值，并預(yù)測(cè)埃博拉病例數(shù)的一種方法.

利用塞拉利昂和尼日利亞的一元線性回歸方程作預(yù)測(cè)分析，也得到同樣的結(jié)果，此處不作累述.

3 討論

從Pearson相關(guān)性分析可以看出，埃博拉的病例數(shù)和積溫差值之間具有極其顯著的相關(guān)性；從積溫差值小于零的頻度和一元線性回歸方程可以看出，積溫差值越大，即2012／1／1到2014／9／7的日平均氣溫低于逐日平均氣溫的頻度越低，病例數(shù)越少；積溫差值越小，即2012／1／1到2014／9／7的日平均氣溫低于逐日平均氣溫的頻度越高，病例數(shù)越多.

從埃博拉病例數(shù)預(yù)測(cè)圖（圖9）可以看出，埃博拉病例數(shù)隨著積溫差值的變化而變化，埃博拉的實(shí)際病例數(shù)與一元回歸預(yù)測(cè)方程的病例數(shù)高度同步，保持穩(wěn)定.日平均氣溫偏離逐日平均氣溫的頻度越小，大規(guī)模發(fā)生埃博拉疫情的可能越小，且當(dāng)總積溫差值恢復(fù)到一定閾值附近，埃博拉病例數(shù)將趨于穩(wěn)定，隨后逐步平息.由此可以推測(cè)日平均氣溫與逐日平均氣溫同期偏離不大、近年氣候比較穩(wěn)定的地區(qū)即便出現(xiàn)輸入性散發(fā)病例，也不容易發(fā)生大規(guī)模擴(kuò)散.

4 結(jié)語

通過對(duì)西非六地和美國兩地的積溫差值的研究，發(fā)現(xiàn)在全球變暖這一大背景下，氣溫變化不符合這一規(guī)律的地方較容易發(fā)生有一定規(guī)模的疫情，在疾病的預(yù)測(cè)預(yù)警中，應(yīng)該認(rèn)識(shí)并考慮到氣候變化對(duì)疫病發(fā)生影響的重要性.筆者只采用了積溫差值這一氣象因子進(jìn)行研究，其他氣象因子與埃博拉病例數(shù)的相關(guān)度有待進(jìn)一步研究.

［1］陸晨.腦卒中與氣象要素關(guān)系的探討［J］.氣象科技，2003，31（6）：397－399.

［2］Wolfram氣象數(shù)據(jù)官方網(wǎng)站.Wolfram Language ＆System Documentation Center［OL］.http：／／reference.wolfram.com／language／note／WeatherDataSourceInformation.html.

［3］邱飛，饒仕吉.年積溫演變狀態(tài)與氣候分析［J］.貴州氣象，2006，30（2）：25－26.

［4］朱星宇，陳勇強(qiáng).SPSS多元統(tǒng)計(jì)分析方法及應(yīng)用［M］.北京：清華大學(xué)出版社，2011：39－54，136－208.