青海田鼠鼠疫流行情況分析與預(yù)測(cè)

王芳+譚穎

利用青海田鼠體蚤檢測(cè)數(shù)據(jù)和青海田鼠鼠巢監(jiān)測(cè)數(shù)據(jù)，對(duì)鼠害的影響因子進(jìn)行因子分析化為較少因子后，利用因子得分作為自變量，分別以細(xì)菌性檢驗(yàn)陽(yáng)性率、牧犬血清陽(yáng)性率為因變量建立鼠巢監(jiān)測(cè)數(shù)據(jù)逐步回歸模型。以細(xì)菌性檢驗(yàn)陽(yáng)性率、田鼠血清陽(yáng)性率為因變量建立田鼠體蚤檢測(cè)數(shù)據(jù)多元回歸模型。

鼠疫流行因子分析逐步回歸模型多元回歸模型

青海田鼠疫源地是四川省1997年由原四川省衛(wèi)生防疫站在石渠縣首先發(fā)現(xiàn)的，2000年后被納入國(guó)家級(jí)鼠疫監(jiān)測(cè)點(diǎn)，2001年青海省稱多縣證實(shí)為青海田鼠鼠疫自然疫源地。了解和評(píng)估石渠縣的鼠疫流行狀況，對(duì)影響鼠疫流行的因素進(jìn)行科學(xué)地評(píng)估以及對(duì)預(yù)測(cè)未來(lái)鼠疫的變化態(tài)勢(shì)，對(duì)于四川省未來(lái)的鼠疫防治工作具有重要的意義。為此，根據(jù)實(shí)際需要，對(duì)石渠縣2001-2013年鼠疫數(shù)據(jù)進(jìn)行分析，建立數(shù)學(xué)模型。

1材料與方法

1.1材料

1.1.1數(shù)據(jù)來(lái)源

四川省石渠縣青海田鼠鼠疫自然疫源地鼠疫監(jiān)測(cè)數(shù)據(jù)。

1.2方法

1.2.1影響鼠疫流行因素篩選方法

利用因子分析法對(duì)青海田鼠鼠巢監(jiān)測(cè)數(shù)據(jù)指標(biāo)進(jìn)行篩選。

1.2.1數(shù)學(xué)建模方法

采用青海田鼠鼠巢檢測(cè)數(shù)據(jù)篩選出影響鼠疫流行的因子，進(jìn)行因子分析后將影響化為少量的因子后，用銀子得分代替影響因子作為自變量，分別用牧犬血清陽(yáng)性率、細(xì)菌性檢測(cè)陽(yáng)性率作為因變量建立逐步回歸方程；采用青海田鼠體蚤檢測(cè)數(shù)據(jù)篩選出影響鼠疫流行的因子，進(jìn)行因子分析后將影響化為少量的因子后，用因子得分代替影響因子作為自變量，分別用細(xì)菌性檢測(cè)陽(yáng)性率、田鼠血清檢驗(yàn)陽(yáng)性率作為因變量建立多元回歸方程。

2模型的建立

2.1鼠巢監(jiān)測(cè)數(shù)據(jù)分析

2.1.1數(shù)據(jù)標(biāo)準(zhǔn)化和指標(biāo)篩選

通過(guò)查閱相關(guān)文獻(xiàn)，結(jié)合實(shí)際應(yīng)用，選取巢數(shù)、蚤指數(shù)、細(xì)鉤蓋蚤數(shù)、直緣雙蚤指明亞種數(shù)、五側(cè)纖蚤鄰近亞種數(shù)作為影響因子。作Z標(biāo)準(zhǔn)化后進(jìn)行因子分析，提取3個(gè)因子，其特征根分別為：2.231、1.487和1.139，均大于1，且這3個(gè)因子的累積方差貢獻(xiàn)率為97.144%。因此選擇這3個(gè)因子。第一個(gè)因子中系數(shù)絕對(duì)值大的主要有直緣雙蚤指明亞種數(shù)、五側(cè)纖蚤鄰近亞種數(shù)，這兩個(gè)指標(biāo)主要反映鼠疫主要的傳播媒介；第二個(gè)因子中系數(shù)絕對(duì)值大的有：細(xì)鉤蓋蚤數(shù)，這個(gè)指標(biāo)是最重要的鼠疫傳播媒介；第三個(gè)因子中系數(shù)絕對(duì)大的主要有取巢數(shù)、蚤指數(shù)，這兩個(gè)指標(biāo)反映的是鼠巢監(jiān)測(cè)的總體信息。旋轉(zhuǎn)后的3個(gè)因子之間的相關(guān)系數(shù)為0，可以明顯的看出各因子之間互不相關(guān)，說(shuō)明提取的四個(gè)因子是合理的、科學(xué)的。

因子分析結(jié)束后，可以將原始觀測(cè)量轉(zhuǎn)換成這些因子的的得分，用這些因子得分代替原來(lái)的觀測(cè)量進(jìn)行進(jìn)行下一步的統(tǒng)計(jì)分析。

2.1.2牧犬血清陽(yáng)性率逐步回歸模型

用因子分析篩選出來(lái)的3個(gè)因子作為自變量，牧犬血清檢測(cè)陽(yáng)性率作為因變量，數(shù)據(jù)集中觀測(cè)值在這3個(gè)因子上的得分作為數(shù)據(jù)輸入，并采用F檢驗(yàn)，以確保每次引入新的變量之前回歸方程中只包含顯著的變量，建立逐步回歸分析模型。其R為0.702，R方為0.492，調(diào)整R方為0.446。模型的F值為10.668，P值為0.008。模型變量系數(shù)分別為：方程的常量為0.069，第一個(gè)因子的系數(shù)為0.045，常量的t值為5.161，第一個(gè)因子的t值為3.266.

逐步回歸分析法排除了第二個(gè)因子和第三個(gè)因子，說(shuō)明牧犬血清陽(yáng)性率與第一個(gè)因子有關(guān)。根據(jù)上文中對(duì)因子分析結(jié)果的解釋，第一個(gè)因子與直緣雙蚤指明亞種數(shù)、五側(cè)纖蚤鄰近亞種數(shù)這兩個(gè)指標(biāo)有關(guān)。因此，牧犬血清陽(yáng)性率主要與直緣雙蚤指明亞種數(shù)、五側(cè)纖蚤鄰近亞種數(shù)這兩個(gè)指標(biāo)相關(guān)。

2.1.3細(xì)菌性檢測(cè)陽(yáng)性率逐步回歸模型

用因子分析篩選出來(lái)的3個(gè)因子作為自變量，青海田鼠細(xì)菌性檢測(cè)陽(yáng)性率作為因變量，建立逐步回歸分析模型。其R為0.616，R方為0.380，調(diào)整R方為0.324。模型的F值為6.739，P值為0.025。模型變量系數(shù)分別為：方程的常量為0.007，第一個(gè)因子的系數(shù)為0.003，常量的t值為6.069，第一個(gè)因子的t值為2.596。

逐步回歸分析法排除了第二個(gè)因子和第三個(gè)因子，說(shuō)明牧犬血清陽(yáng)性率與第一個(gè)因子有關(guān)。根據(jù)上文中對(duì)因子分析結(jié)果的解釋，第一個(gè)因子與直緣雙蚤指明亞種數(shù)、五側(cè)纖蚤鄰近亞種數(shù)這兩個(gè)指標(biāo)有關(guān)。因此，青海田鼠細(xì)菌檢驗(yàn)陽(yáng)性率主要與直緣雙蚤指明亞種數(shù)、五側(cè)纖蚤鄰近亞種數(shù)這兩個(gè)指標(biāo)相關(guān)。

2.2青海田鼠體蚤數(shù)據(jù)分析

2.2.1數(shù)據(jù)指標(biāo)篩選

通過(guò)查閱相關(guān)文獻(xiàn)可知，影響鼠害流行的主要因素有：主要宿主和主要傳播媒介。青海田鼠鼠疫自然疫源地鼠疫的主要宿主為青海田鼠，主要傳播蚤類為細(xì)鉤蓋蚤和直緣雙蚤指明亞種。因此選取捕鼠數(shù)、鼠密度、染蚤率、蚤指數(shù)、直緣雙蚤指明亞種數(shù)、細(xì)鉤蓋蚤數(shù)為影響因子。作Z標(biāo)準(zhǔn)化后進(jìn)行因子分析，得到2個(gè)因子，其特征根分別為：2.741和2.594，均大于1，且這3個(gè)因子的累積方差貢獻(xiàn)率為88.92%。因此選擇這2個(gè)因子。

查看各成份在2個(gè)因子載荷矩陣可知：第一個(gè)因子中系數(shù)絕對(duì)值大的主要有鼠密度、染蚤率、蚤指數(shù)，第二個(gè)因子中系數(shù)絕對(duì)值大的有：捕鼠數(shù)、直緣雙蚤指明亞種數(shù)、細(xì)鉤蓋蚤數(shù)。旋轉(zhuǎn)后的2個(gè)因子之間的相關(guān)系數(shù)為0，可以明顯的看出各因子之間互不相關(guān)，說(shuō)明提取的四個(gè)因子是合理的、科學(xué)的。

2.2.2青海田鼠體蚤染蚤率和蚤指數(shù)相關(guān)性分析

青海田鼠體蚤染蚤率和蚤指數(shù)變化情況如圖1所示：

圖1青海田鼠體蚤染蚤率和蚤指數(shù)變化圖

用典型相關(guān)分析法分析這兩個(gè)變量的相關(guān)性，變量染蚤率和蚤指數(shù)的相關(guān)系數(shù)為0.799，說(shuō)明兩個(gè)變量顯著相關(guān)。

2.2.3細(xì)菌性檢驗(yàn)陽(yáng)性率多元回歸模型

用因子分析篩選出來(lái)的2個(gè)因子作為自變量，細(xì)菌性檢驗(yàn)陽(yáng)性率作為因變量，細(xì)菌學(xué)檢驗(yàn)陽(yáng)性率作為因變量，青海田鼠體蚤觀測(cè)值作為數(shù)據(jù)輸入，并采用F檢驗(yàn)，建立多元回歸分析模型。其R為0.588，R方為0.345，調(diào)整R方為0.214。模型的F值為2.638，P值為0.120 .模型變量系數(shù)分別為：方程的常量為5.632，第一個(gè)因子的系數(shù)為-0.387，t值為-1.514，第二個(gè)因子的系數(shù)為-0.442，t值為-1.728.

2.2.4血清學(xué)檢驗(yàn)陽(yáng)性率多元回歸模型

用因子分析篩選出來(lái)的2個(gè)因子作為自變量，血清學(xué)檢驗(yàn)陽(yáng)性率作為因變量，細(xì)菌學(xué)檢驗(yàn)陽(yáng)性率作為因變量，青海田鼠體蚤觀測(cè)值作為數(shù)據(jù)輸入，并采用F檢驗(yàn)，建立多元回歸分析模型。其R為0.748，R方為0.560，調(diào)整R方為0.472.模型的F值為6.359，P值為0.017。模型變量系數(shù)分別為：方程的常量為4.565，第一個(gè)因子的系數(shù)為-0.238，t值為-1.136，第二個(gè)因子的系數(shù)為0.709，t值為0.007。

3總結(jié)

本文分別對(duì)青海田鼠體蚤監(jiān)測(cè)數(shù)據(jù)和青海田鼠鼠巢檢測(cè)數(shù)據(jù)進(jìn)行多元統(tǒng)計(jì)分析，并得出相應(yīng)的回歸模型，對(duì)青海田鼠鼠疫疫情的監(jiān)測(cè)和預(yù)測(cè)有重要的知道作用。參考文獻(xiàn)：

[1]黃建忠，鄧祖軍.2000-2003年四川省石渠縣鼠疫自然疫源地初探[J].預(yù)防醫(yī)學(xué)情報(bào)雜志，2015，31（1）：26-32.

[2]汪立茂，宋曉玉，李光清等.2000-2004四川省青海田鼠鼠疫流行特點(diǎn)分析[J].預(yù)防醫(yī)學(xué)情報(bào)雜志，2006，22（2）：146-149.

[3]劉振才，周曉磊，張博宇等.動(dòng)物鼠疫預(yù)測(cè)模型及預(yù)警指標(biāo)的建立[J].中國(guó)地方病防治雜志，2015，30（1）：1-3.

[4]劉振才，海榮，李富忠等.青藏高原青海田鼠鼠疫自然疫源地的發(fā)現(xiàn)與研究[J].中國(guó)地方病防治雜志，2001，16（1）：321-327.

[5]祁騰，楊孔，汪立茂等.石渠縣2001-2013年鼠疫疫源地流行病學(xué)分析[J].中國(guó)人獸共患病學(xué)報(bào)，2015，31（5）：485-488.

[6]張立軍，任英華.多元統(tǒng)計(jì)分析實(shí)驗(yàn)[M].北京：中國(guó)統(tǒng)計(jì)出版社，2009.122-141.

[7]何曉群，劉文卿.應(yīng)用回歸分析[M].北京：中國(guó)人民大學(xué)出版社，2011.56-77.

[8]汪立茂，李光清，許光榮.青海田鼠鼠疫流行特點(diǎn)分析[J].中國(guó)地方病防治雜志，2004，19（3）：178-181.

[9]楊俊平.沙鼠鼠疫預(yù)測(cè)的數(shù)學(xué)模型[J].內(nèi)蒙古科技與經(jīng)濟(jì)，2001（1）：68.

[10]李仲來(lái)，周方孝等.達(dá)烏爾黃鼠鼠疫預(yù)測(cè)預(yù)報(bào)的數(shù)學(xué)模型（Ⅶ）[J].中國(guó)地方病防治雜志，200217（3）： 129-131.

基金項(xiàng)目：中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目，優(yōu)秀學(xué)生培養(yǎng)工程項(xiàng)目，項(xiàng)目編號(hào)：2016ZYXS16。