1.3 计算模型与分析方法 (1)在本研究中,鸡群数量为 n,鸡的平均质量为 m,每批鸡肉质量为 mt,则在每个鸡群选择屠宰鸡的数量为nk=Round[mt(/ nm)] 式中 k=1,2,…,n。 (2)鸡群疫病流行率为 p1,从感染鸡群随机选择鸡制成的鸡肉产品传带禽伤寒沙门氏菌的概率为 q0,则该批鸡肉携带疫病的概率为p=1-{1-p1 [1-(1-q0)nk]} (3)定义: X1=1 表示从感染鸡群选择某只鸡携带病菌,0表示未携带病菌; X2=1 表示该鸡在屠宰后携带病菌,0 表示未携带病菌; X3=1 表示分割后的鸡肉携带病菌,0 表示未携带病菌; X4=1 表示鸡肉在冷冻后携带病菌,0 表示未携带病菌; Y1=1 表示该鸡在宰前检验中被认为是病鸡,0表示不被认为是病鸡; Y1=1 表示该鸡在宰后检验中被认为是病鸡,0表示不被认为是病鸡。 根据全概率公式 q0=P(X4=1,Y=0) = ∑ P(X4=1,Y2=0|X2=i,Y2=0)P(X2=i,Y2=0) (3)式中:P(X4=1,Y2=0|X2=1,Y2=0)=p8(1-p10) (4) P(X4=1,Y2=0|X2=1,Y2=0)=p9(1-p10) (5)利用全概率公式P(X2=1,Y2=0)= ∑ P(X2=1,Y2=0|X1=i,Y1=0)P(X1=i,Y1=0) =(1-p6p7)p2(1-p3p4)+p5(1-p6p7)(1-p2)(1-p3pf)(6)P(X2=0,Y2=0)= ∑ P(X2=0,Y2=0|X1=i,Y1=0)P(X1=i,Y1=0) =P(X2=0,Y2=0|X1=0,Y1=0)P(X1=0,Y1=0) =(1-p5)(1-p6pf)(1-p2)(1-p3pf) (7)将式(4)(5)(6)带入式(3),并化简,即得q0=(1-p10){p8(1-p6 p7)[p2(1-p3 p4)+p5(1-p2)(1-p6 pf)]+ p9(1-p2)(1-p3 pf)p5(1-p6 pf)} (8) 利用上述数学模型及确定的风险参数值,采用固定参考值和蒙特卡洛方法,模拟分析进口鸡肉携带禽伤寒沙门氏菌的风险水平,以及不同风险因子对进口鸡肉携带禽伤寒沙门氏菌风险的影响。 2、结果 2.1 固定参考值分析 将表 1 中的参考值带入式(2)和式(8),得 q0=0.367 4,p=0.001 0,即在设定的参考值条件下,进口鸡肉携带禽伤寒沙门氏菌的概率为0.1%,也即每 1 000 批次进口产品中携带该病菌的批次数期望值为 1。 2.1.1 各风险因子对感染鸡群生产鸡肉带菌率的影响 通过固定其他参数值,分析表 1 中相关参数对感染鸡群生产的鸡肉产品传带禽伤寒沙门氏菌概率 [ 式(8)中 q0] 的影响。结果显示,在设定的风险参数取值范围内,q0 值随病菌感染率和鸡肉带菌率的上升而线性上升(图 2~3)。当鸡群感染率达100% 时,产品携带病菌的概率为 48.6%;当感染鸡鸡肉带菌率为 100% 时,产品携带病菌的概率达73.4%。其余风险参数对 q0 值影响很小。
|
上一篇:我国水产品质量安全的风险评估
下一篇:养禽场感染压力的风险评估