流調(diào)、監(jiān)測工作中常用的統(tǒng)計分析方法

一、統(tǒng)計分析的重要性

統(tǒng)計分析是流行病學調(diào)查靈魂所在。在動物疫病流行病學調(diào)查中,

統(tǒng)計分析具有以下四個方面的重要意義。

1、統(tǒng)計分析是信息交流的需要

2、是調(diào)查方案設計的需要

3、是正確分析調(diào)查結果的需要

4、是深入分析調(diào)查結果的需要

二、統(tǒng)計分析在流行病學調(diào)查中的應用

1、統(tǒng)計數(shù)據(jù)的計算

2、分布情況的描述

3、確定樣本量

4、從樣本的測量值估計總體情況

5、分析各組數(shù)據(jù)之間的差異是否顯著

6、分析某些變量是否相關

三，流行病學調(diào)查中常用的統(tǒng)計分析方法

7,料漳量確實定

(1)至少發(fā)現(xiàn)1個陽性病例的樣本量

公式n=ln(1-a)/In(1-pXs)

In為自然對數(shù)，a為置信區(qū)間，p為預測感染率，s為檢測的靈敏

度。

如預測某豬場豬瘟感染率為10%,RT-PCR檢測的敏感性為90%,在置

信區(qū)間為95%,確保至少能檢出一例陽性動物需采集多少樣本？

n=ln(l-0.95)/ln(l-0.1X0.9)=32

假設32份樣品檢測結果全部為陰性，說明該豬場豬瘟感染率低于

10%o

(2)調(diào)查總體率的樣本量計算

在調(diào)查某種疫病在某個群體中的流行率或免疫抗體合格率。

公式：n=ZXZX(1-WXL)4-(WXLXeX￡)

n為需要采集的樣本量，Z是置信度(P)的標準正態(tài)分布參數(shù)(P=95%

時，Z=1.96；P=99%時，Z=2.58),w為估計的流行率或合格率,L為檢測

的靈敏度，￡為最大的相對誤差。

如某縣豬0型口蹄疫免疫合格率估計為70%,間接血凝檢測的靈敏度

為9096,最大相對誤差為10%,置信度P=95%,那么檢測該縣豬0型口蹄

疫免疫合格率，需要抽取多少樣本？

n=L96X1.96X(l-70%X90%).(0.7X90%X10%X10%)

4226

如估計的免疫合格率為90%,那么需采集90頭份，估計合格率為50%,

那么需采集470頭份。

對于養(yǎng)殖場采樣數(shù)量可參考下表：

置信度為9596,按估計的陽性率/合格率、群的大小計算的采樣數(shù)量

流行百分率%

504030252015105210.50.1

2046791012161920202020

3046891114192630303030

40568101215213140404040

50568101216223546505050

60568101216233855606060

705G81013172440G2707070

80568101317244268798080

90568101317254373879090

10056910131725457896100100

150569111318274995130148150

2005691113182751105155190200

50056911M192856129225349500

10005691114192957138258450950

2、足計撤據(jù)的計算

發(fā)病率、死亡率、合格率、陽性率、平均值、離散度等指標與去年

同期比擬（同比）、與上個月/季度比擬（環(huán)比）、與國家要求比擬（基線

比），來描述疫情的現(xiàn)狀和開展趨勢。

Excell軟件的應用：在Excell工具欄中點加載宏，在彈出的對話框

中選中分析工具庫和分析工具庫-VBA函數(shù)，確認，以后使用t檢驗、相

關系數(shù)分析、方差分析等在工具欄中選中數(shù)據(jù)分析，然后選擇相應分析

方法即可。

Excell中最常用的統(tǒng)計分析方法是自己輸入相應的計算公式和選擇

插入函數(shù)。較為常用的函數(shù)有：

平均值一average-在插入函數(shù)中選定average,選定統(tǒng)計區(qū)域，確定；

標準差一stdev-同上；

離散度=標準差?平均值義100

樣本數(shù)一count-同平均值的方法；

合格數(shù)一countif-在插入函數(shù)中選定countif,選定統(tǒng)計區(qū)域，設

定合格參數(shù)，確認；

合格率=合格數(shù):樣本數(shù)X100。

免疫抗體合格率描述了群體對每一種疫病免疫狀態(tài)的一個層面，不

能反映其整體面貌。合格率、平均值、離散度相結合根本可以反映群體

的免疫狀況。

例1：A養(yǎng)雞場新城疫檢測效價分別為：4、4、4、4、4、4、4、4、

3、3,合格率為80%；

B場分別為：8、8、8、8、8、8、8、8、3、3,合格率也為80%；但

A場抗體效價平均值為3.8,而B場平均值為7.0,說明A場已低于保護臨

界線，應盡快進行補免，而B場僅個別雞只低于保護臨界線，需密切觀

察。

例2：A場新城疫免疫抗體分別為：5、5、5、5、5、6、6、6、3、3,

合格率為80%,平均值為4.9。

B場：11、12、4、4、4、、4、，5、5、0、0,合格率為80%,平均值

為4.9o

但是A場離散度為22.46,而B場為80.17,什么A場免疫狀況良好,

只是近期需要補免而已，而B場離散度較高，極有可能現(xiàn)在正在發(fā)生該

病的疫情。

序號A雞場B雞場C雞場D雞場

148511

248512

34854

44854

54854

64864

74865

84865

93330

103330

平均值3.807.004.904.90

標準差0.422.111.103.93

合格數(shù)8888

離散度11.1030.1222.4680.17

樣本數(shù)10101010

合格率80.0080.0080.0080.00

離散度低于40%,表示抗體水平較為一致，整齊度高；介于40%—

60%之間，表示個體間有一定的差異，如此時平均值較高，那么說明免

疫狀況良好，如平均值略高于保護臨界值，應立即補免；介于60%—80%

之間，說明個體間離散度較大，群體有可能已被感染；離散度大于80%,

那么表示個體間差異非常大，多數(shù)情況為群體已被感染。

合格率符合要求，平均值越高于保護臨界值、離散度越低，表示群

體免疫越好。

3.就針臺新演弧的刷作

應用Excell可以制作各種酶聯(lián)方法檢測出的0D值換算為抗體值或

阻斷率模板，以利檢測結果的計算。還可以制作檢測結果統(tǒng)計分析模板,

只要將檢測結果分組輸入，即可自動計算出合格率/陽性率、平均值、離

散度等相關結果和圖表。

久寺祈各俶盛據(jù)之間埼妻異是否支得

4.1、卡方檢驗

x2=[(ad_bc)2]Xn4-[(a+b)X(c+d)X(a+c)X(b+d)]

例1:甲乙2個養(yǎng)殖場，在某一時段的存欄數(shù)和死亡數(shù)見下表，試分

析2個場間的死亡率是否存在差異。

場名死亡數(shù)存活數(shù)存欄數(shù)死亡率

甲30(a)250(b)28010.7

乙22(c)398(d)4205.2

合計526487007.4

n=a+b+c+d

X2=7.32

從卡方分布界值表查得，a=0.05,v=l時，卡方值為3.84；a=0.01,

v=l時，卡方值為6.63,2個臨界值都小于本例計算的卡方值。所以，無

論在a=0.05水平上，還是在a=0.01水平二，甲乙2個養(yǎng)殖場的死亡率

都有顯著差異。

4.2、t檢驗

用于2組數(shù)據(jù)差異性檢驗。如同場、相同飼養(yǎng)條件下，60頭體重相

近的仔豬隨機分為2組，1組使用豬瘟細胞苗，1組使用豬瘟脾淋苗，免

疫后20天分別采集血清進行免疫效價測定，比擬哪種疫苗免疫的效價水

三高即可應用t檢驗。

編號12345678910

匆胞苗組5465645645

脾淋苗組6756754766

使用Excell表有2種進行t檢驗的方法。

A、上述數(shù)據(jù)輸入Excell表中，在菜單“插入”欄選擇“函數(shù)”，從

口選擇“TTEST”，將上述數(shù)據(jù)輸送到TTEST的界面中，尾值（Tails）設

置為2（雙側(cè)檢驗），上例檢驗類型［Type）設置為2（屬于等方差的2

人獨立樣本t檢驗）【檢驗類型（Type）設置為1〔實驗動物為1組）（屬

二配對資料的t檢驗，一般用于使用某項技術前與使用后的比照）】。

計算結果為0.0401,此值小于0.05,即P<0.05,說明豬瘟脾淋苗免

疫豬產(chǎn)生的免疫抗體與細胞苗產(chǎn)生的抗體有顯著差異。

如P<0.01,那么差異極顯著。

B、上述數(shù)據(jù)輸入Excell表中，在工具欄選中“數(shù)據(jù)分析”，選中“t

檢驗等方差假設”，將上述數(shù)據(jù)輸如即可。

P（T<=t）雙尾值為0.0401。

A法中：檢驗類型（Type）設置為1（屬于配對資料的t檢驗，一般

月于使用某項技術前與使用后的比照）的在B法中使用“t檢驗-平均值

的成對二樣本分析:

各組樣本數(shù)超過100時，也可使用Z檢驗，多組之間差異性分析可

以用t檢驗、Z檢驗，也可以用q檢驗。

5、今祈慕的愛/是否依關

相關系數(shù)：是用以反映變量之間相關關系密切程度的統(tǒng)計指標。相

關系數(shù)是按積差方法計算，同樣以兩變量與各自平均值的離差為根底，

通過兩個離差相乘來反映兩變量之間相關程度；著重研究線性的單相關

系數(shù)。

相關系數(shù)分為011正相關1和（TTi負相關）。一般來說相關系數(shù)

的絕對值大于0.7,表示高度相關，0.5-0.7表示中度相關，0.3-0.5為

弱相關，小于0.3表示不相關。

例1

2007年我們調(diào)查的豬“高熱病”發(fā)病率、死亡率以及豬場豬瘟免疫

合格率見下表。

序號豬瘟免疫合格率高熱病發(fā)病率高熱病死亡率

19500

296.772010

377.366050

在Excell表中錄入上述數(shù)據(jù)，選擇工具欄中的“數(shù)據(jù)分析）在對

話框中選擇“相關系數(shù)”，輸入相關數(shù)據(jù)，即得以下數(shù)據(jù)。

豬瘟免疫合格率高熱病發(fā)病率高熱病死亡率

豬瘟免疫合格率1

高熱病發(fā)病率-0.9147002091

高熱病死亡率-0.9630517140.9897433191

豬瘟免疫合格率與高熱病的發(fā)病率、死亡率的相關系數(shù)分別為

-0.9K-0.96,絕對值均大于0.7,說明豬瘟免疫合格率與高熱病的發(fā)病

玄和死亡率之間存在高度的負相關，也就是說養(yǎng)豬場豬瘟免疫的越好，

高熱病的發(fā)病率和死亡率就越低。

6、比值比和風險比

病例對照調(diào)查和暴露對照調(diào)查可以獲得重要的原始數(shù)據(jù)，依據(jù)這些

原始數(shù)據(jù),可以計算比值比（0R）、風險比（RR）。

OR值得意義

OR值被調(diào)查的因素與被調(diào)查的疾病的關聯(lián)程度

有較強的關聯(lián)

<0.4

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，且可能性較大

有中等程度的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，可能性一般

有較弱的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，但可能性較小

無關聯(lián)

有較弱的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，但可能性較小

有中等程度的關聯(lián)

1.5-3

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，可能性一般

有較強的關聯(lián)

>3

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，且可能性較大

RR值的意義

RR值被調(diào)查的因素與被調(diào)查的疾病的關聯(lián)程度

有很強的關聯(lián)

<0.1

被請查的因素可能防止被調(diào)查疾病的發(fā)生，且可能性很大

有較強的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，且可能性較大

有中等程度的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，可能性一般

有較弱的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，但可能性較小

無關聯(lián)

有較弱的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，但可能性較小

有中等程度的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，可能性一般

有較強的關聯(lián)

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，且可能性較大

有很強的關聯(lián)

>10

被調(diào)查的因素可能防止被調(diào)查疾病的發(fā)生，且可能性很大

A、OR值的計算

0R=(A+B)+(C+D)

暴露非暴露合計

病例ABG=A+B

非病例CDH=C+D

合計E=A+CF=B+DN=A+B+C+D

例1：調(diào)查奶?？谔阋甙l(fā)生情況時，設定一組為發(fā)病牛場，二組是沒

有發(fā)病的牛場，各組含有的牛場規(guī)模、飼養(yǎng)方式、衛(wèi)生狀況相似，都調(diào)

查了50個牛場。結果為：發(fā)病組接種0、亞I、A型疫苗的為36個，未

發(fā)病組接種上述疫苗的為