统计分析在动物疫病流行病学调查中的应用

流调、监测工作中常用的统计分析方法

一、统计分析的重要性

统计分析是流行病学调查灵魂所在。在动物疫病流行病学调查中，统计分析具有以下四个方面的重要意义。

1、统计分析是信息交流的需要

2、是调查方案设计的需要

3、是正确分析调查结果的需要

4、是深入分析调查结果的需要

二、统计分析在流行病学调查中的应用

1、统计数据的计算

2、分布情况的描述

3、确定样本量

4、从样本的测量值估计总体情况

5、分析各组数据之间的差异是否显著

6、分析某些变量是否相关

三，流行病学调查中常用的统计分析方法

1、样本量的确定

(1)至少发现1个阳性病例的样本量

公式n=ln(1-a)/ln(1-p×s)

ln为自然对数，a为置信区间，p为预测感染率，s为检测的灵敏度。

如预测某猪场猪瘟感染率为10%，RT-PCR检测的敏感性为90%，在置信区间为95%，确保至少能检出一例阳性动物需采集多少样本？

n=ln(1-0.95)/ln(1-0.1×0.9)=32

假如32份样品检测结果全部为阴性，表明该猪场猪瘟感染率低于10%。

(2)调查总体率的样本量计算

在调查某种疫病在某个群体中的流行率或免疫抗体合格率。

公式：n=Z×Z×(1-w×L)÷(w×L×ε×ε)

n为需要采集的样本量，Z是置信度（P）的标准正态分布参数（P=95%时，Z=1.96；P=99%时，Z=2.58），w为估计的流行率或合格率，L为检测的灵敏度，ε为最大的相对误差。

如某县猪O型口蹄疫免疫合格率估计为70%，间接血凝检测的灵敏度为90%，最大相对误差为10%，置信度P=95%，则检测该县猪O型口蹄疫免疫合格率，需要抽取多少样本？

n=1.96×1.96×(1-70%×90%)÷(0.7×90%×10%×10%)

n≈226

如估计的免疫合格率为90%，则需采集90头份，估计合格率为50%，则需采集470头份。

对于养殖场采样数量可参考下表：

置信度为95%，按估计的阳性率/合格率、群的大小计算的采样数量

群的大小流行百分率%

50 40 30 25 20 15 10 5 2 1 0.5 0.1

20 4 6 7 9 10 12 16 19 20 20 20 20

30 4 6 8 9 11 14 19 26 30 30 30 30

40 5 6 8 10 12 15 21 31 40 40 40 40

50 5 6 8 10 12 16 22 35 46 50 50 50

60 5 6 8 10 12 16 23 38 55 60 60 60

70 5 6 8 10 13 17 24 40 62 70 70 70

80 5 6 8 10 13 17 24 42 68 79 80 80

90 5 6 8 10 13 17 25 43 73 87 90 90

100 5 6 9 10 13 17 25 45 78 96 100 100 150 5 6 9 11 13 18 27 49 95 130 148 150 200 5 6 9 11 13 18 27 51 105 155 190 200 500 5 6 9 11 14 19 28 56 129 225 349 500 1000 5 6 9 11 14 19 29 57 138 258 450 950

2、统计数据的计算

发病率、死亡率、合格率、阳性率、平均值、离散度等指标

与去年同期比较（同比）、与上个月/季度比较（环比）、与国家要求比较（基线比），来描述疫情的现状和发展趋势。

Excell软件的应用：在Excell工具栏中点加载宏，在弹出的对话框中选中分析工具库和分析工具库-VBA函数，确认，以后使用t检验、相关系数分析、方差分析等在工具栏中选中数据分析，然后选择相应分析方法即可。

Excell中最常用的统计分析方法是自己输入相应的计算公式和选择插入函数。较为常用的函数有：

平均值－average-在插入函数中选定average，选定统计区域，确定；

标准差－stdev-同上；

离散度＝标准差÷平均值×100

样本数－count-同平均值的方法；

合格数－count if-在插入函数中选定count if，选定统计区域，设定合格参数，确认；

合格率＝合格数÷样本数×100。

免疫抗体合格率描述了群体对每一种疫病免疫状态的一个层面，不能反映其整体面貌。合格率、平均值、离散度相结合基本可以反映群体的免疫状况。

例1：A养鸡场新城疫检测效价分别为：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.9。

但是A场离散度为22.46，而B场为80.17，什么A场免疫状况良好，只是近期需要补免而已，而B场离散度较高，极有可能现在正在发生该病的疫情。

序号A鸡场B鸡场C鸡场D鸡场

1 4 8 5 11

2 4 8 5 12

3 4 8 5 4

4 4 8

5 4

5 4 8 5 4

6 4 8 6 4

7 4 8 6 5

8 4 8 6 5

9 3 3 3 0

10 3 3 3 0

平均值 3.80 7.00 4.90 4.90

标准差0.42 2.11 1.10 3.93

合格数8 8 8 8

离散度11.10 30.12 22.46 80.17

样本数10 10 10 10

合格率80.00 80.00 80.00 80.00

离散度低于40％，表示抗体水平较为一致，整齐度高；介于40％－60％之间，表示个体间有一定的差异，如此时平均值较高，则说明免疫状况良好，如平均值略高于保护临界值，应立即补免；介于60％－80％之间，表明个体间离散度较大，群体