猪的健康指数模型完整版.docx
- 文档编号:24669530
- 上传时间:2023-05-30
- 格式:DOCX
- 页数:16
- 大小:202.49KB
猪的健康指数模型完整版.docx
《猪的健康指数模型完整版.docx》由会员分享,可在线阅读,更多相关《猪的健康指数模型完整版.docx(16页珍藏版)》请在冰豆网上搜索。
猪的健康指数模型完整版
HENsystemofficeroom【HEN16H-HENS2AHENS8Q8-HENH1688】
猪的健康指数模型
2011年商丘师范学院数学建模模拟练习
承诺书
我们仔细阅读了商丘师范学院数学建模模拟练习的竞赛规则。
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与本队以外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们愿意承担由此引起的一切后果。
我们的参赛报名号为:
25
参赛组别(本科或专科):
本科
参赛队员(签名):
队员1:
陈玉松
队员2:
王小端
队员3:
郭文娜
2011年商丘师范学院建模模拟练习
编号专用页
参赛队伍的参赛号码:
(请各个参赛队提前填写好):
竞赛统一编号(由竞赛组委会送至评委团前编号):
竞赛评阅编号(由竞赛评委团评阅前进行编号):
2011年商丘师范学院数学建模模拟练习
题目利用多元回归求动物健康指数
一、摘要
本文主要应用数学建模中的多元回归模型,来拟合多个影响因素对一个变量的影响。
利用健康猪群及病态猪群血液成分中白细胞、红细胞、血小板、血红蛋白的数量与中毒指数、病毒指数、细菌指数、过敏指数、免疫指数,考察它们与健康指数的关联性。
用概率论与数理统计的方法找出健康指数和各个因素之间的近似线性关系,建立各个与健康指数的多元线性回归模型。
通过SPSS程序与MATLAB的统计工具箱中的regress命令,拟合出参数估计值,,并利用MATLAB对回归模型的残差及参数置信区间进行分析,找到原模型的一些不足,进而对模型进行改进,得出方程效果良好的结论,指出模型的应用价值。
在此基础上进一步提出检验猪健康状况的标准,同时给饲养者以相应的建议。
二、关键字
多元回归模型加权健康指数残差分析SPSS软件
三、问题的重述
当前猪场一般针对某一疫病进行抗体检测来观察疫苗免疫效果,检测病原或抗原确定是否存在病毒感染,但不能判断猪场猪群的健康情况,亦不能确定疫病可能发生的风险有多大,不能进行快速的诊断。
大量的健康动物机体携带病原,但并不发病。
实验室检测也充分证明,猪群带有猪繁殖与呼吸综合征病毒,即使表现为病毒血症,猪群亦表现为健康状态,这给临床判断带来极大困难。
血液是动物和人体内不可缺少的组成部分,对于平衡机体和环境之间的矛盾起着极其重要的缓冲作用,机体血液中含有各种各样与机体抵抗力相关的成分,血液之中的各种成分会随着机体健康状态的改变而改变。
所以,通过血液成分的变化来判断机体的健康状态就有了科学的依据。
四、符号说明
符号
含义
随机误差
影响健康指数的主要因素
健康指数
对健康指数影响的权重
9
变量因素个数
1
线性方程个数
五、基本假设
1、本文注重影响健康指数的主要因素——白细胞,红细胞,血小板,血红蛋白,中毒指数,病毒指数,细菌指数,过敏指数,免疫指数;大胆假设他们与健康指数的关系依次为线性关系;
2、忽略了次要的及相对微弱因素,例如贫血指数。
为综合评判后的健康指数,
,
…
为各主要因素,
,
…..
为对健康指数影响的权重;
3、各血液成分给予一定的加权,计算出一个综合数值,称为健康指数忽略外界温度、湿度、营养等环境因素的影响。
六、问题的分析
机体血液中含有各种各样与机体抵抗力相关的成分,血液之中的各种成分会随着机体健康状态的改变而改变。
机体健康状态的改变直接影响血液成分的变化,通过测量机体内部各血液成分的变化可以预测机体的健康指数,我们可以给予各血液成分一定的加权,利用多元线性回归模型计算出一个综合数值,即健康指数。
下面首先给出加权的定义及多元线性回归模型。
名词解释:
加权
统计学认为,在统计中计算平均数等指标时,对各个变量值具有权衡轻重作用的数值就称为权数变量大小对平均数起决定作用它的大小决定着平均数的大小。
权数大小对平均数起权衡轻重的作用,它的比重大小影响平均数的大小,使平均数趋于权数大的变量值。
一般说的平均数,就是把所有的数加起来,再除以这些数的总个数。
表示为:
(p1+p2+p3+…..+pn)/n;
但有的数据记录中有一些相同的数据,在计算的时候,那一个数有几个相同数,就把这个数乘上几,这个几,就叫权,加权,就是乘上几后再加。
平均数还是要除以总个数。
还是以上面的各个数为例:
它们每个数都有一些相同数,表示为:
k1,k2,k3…….kn;
加权平均的公式是:
(k1p1+k2p2+k3p3+……knpn)/(k1+k2+k3+…..kn)
对于本题,由于各个因素不同,则考虑健康指数的侧重比重亦不同,但是都是具体常量。
所以我们大胆假设,健康指数与各主要因素之间是线性关系。
多元线性回归:
有一个或一组非随机变量估计或预测某一随机变量的观测值时,所建立的数学模型称为回归模型,根据回归模型进行的统计分析,叫回归分析,如果这个模型是线性的,则成为线性回归分析。
设可控或不可控的因素为
;目标函数为
。
主要因素分段描述:
………………………………………….白细胞所占的加权数
…………………………………………..红细胞所占的加权数。
.................
…………………………………………….免疫指数所占的加权数
则健康指数为
其中
为常数系数。
建立线性模型
................
其中
为测量误差,相互独立。
令
,
,
,
上式称为线性回归方程的数学模型。
利用最小二乘法估计或极大似然估计,最小二乘解
,
,
…..
所满足的方程就是
可得系数
的值。
将得出的系数代入方程
由于实际意义,机体在被细菌或病毒感染时,血液内的细胞数量和形态会立即发生变化,各血液成分给予一定的加权,计算出一个综合数值,称为健康指数。
我们大胆假设,健康指数与各主要因素之间是线性关系。
由于各个因素不同,则考虑健康指数的侧重比重亦不同,但是都是具体常量。
七、模型的建立与求解
综合上面的分析,建立如下的数学模型:
直接通过MALAB统计工具箱得到如下结果:
参数
参数估计值
参数置信区间
a0
[]
a1
[]
a2
[]
a3
[]
a4
[]
a5
[]
a6
[]
a7
[]
a8
[]
a9
[]
F=R2=P=0
表1
应用SPSS软件得到如下结果:
相关性
健康指数
白细胞
红细胞
血小板
血红蛋白
中毒指数
病毒指数
细菌指数
过敏指数
免疫指数
Pearson相关性
健康指数
.311
.183
.965
白细胞
.267
.376
.114
.748
.464
红细胞
.311
.267
.239
.088
.219
血小板
.183
.376
.239
.191
.201
血红蛋白
.114
.071
.083
.479
.428
中毒指数
.071
.549
.117
病毒指数
.083
.549
细菌指数
.748
.088
.191
.479
.594
过敏指数
.464
.428
.117
.594
免疫指数
.965
.219
.201
表2
图1残差分析图
图2
从残差图1可以看出,除第十四、二十、三十二个数据外,其余数据的残差离零点均较近,且残差的置信区间均包含零点,这说明回归模型:
8能较好的符合原始数据,而第十四、二十、三十二个数据可视为异常点。
而表1显示,R2=指因变量健康指数的%可由模型确定,F=,P=0远小于也说明此模型从整体看是可用的。
且通过对已知数据的估计,影响健康指数的其他因素作用都包含在随机变量
中,由图2,随机误差
大致服从均值为零的正态分布。
八、模型的改进
通过对表1表2的观察分析可知参数的置信区间包含零点,说明模型存在一定的不可靠性,且红细胞和中毒指数、免疫指数和过敏指数存在很高的相关性,二者肯存在一定的交互作用,我们可以对模型进行一定的改进,模型记作:
直接通过MALAB统计工具箱得到如下结果:
参数
参数估计值
参数置信区间
a0
[,]
a1
[,]
a2
[,]
a3
[,]
a4
[,]
a5
[,]
a6
[,]
a7
[,]
a8
[,]
a9
[,]
a10
[,]
a11
[,]
F=R2=P=0
表3
图3
8、模型的分析
表3与表2相比,R2有所提高,通过图3可以看出,残差与原模型相比,与零点较接近,说明改进后的模型比原模型要精确。
而且,该模型得出了影响猪的健康指数的各种因素的权重,给如何判定猪的健康状况提供了一个可靠的依据,且如果实验正确,对猪群的健康检测不会有太大的误差。
九、模型的优缺点
优点
1、本文的基础为多元回归算法,由于有成熟的算法,使得计算过程十分流畅,算法效率很高。
2、忽略了很多次要的及相对微弱因素,例如贫血指数,尽量减少未知量,使计算过程更加简便。
缺点
1、忽略了其他因素之间的相关性,只考虑到相关性很高的红细胞与中毒指数,过敏指数与免疫指数;
2、该模型给定了一定的假设条件,存在一定的理想性,可能与实际有一定的偏差。
十、参考文献
[1]姜启源等,数学模型(第三版),高等教育出版社
[2]缪铨生等,概率与统计(第三版),华东师范大学出版社
[3]萧树铁等,数学实验(第二版),高等教育出版社
附录一:
1、利用MATLAB统计工具箱得到初步的回归方程
设回归方程为:
实现本功能的MATLAB程序为:
>>B=[;;15860021422780171326267344547011337712821749944648126782396312273116613123471330606579180;;000000000000;0000000000;000000000000;00000000;];
>>A=[;;;;000000000000000000;0000000000000000000;000000000000000000;000000000000000;];
>>C=[A,B]
C=
+003*
Columns1through9
00000000
000000000
00000000
0000000
Columns10through18
00000000
00000000
00000000
000000
Columns19through27
000
0000
0000000
00000
Columns28through36
00
00
00000
0000
Columns37through45
0000
00000
00
0
Columns46through50
00000
0
0
>>c=C';
>>a=ones(50,1);
>>X=[a,c];
>>Y=[];
>>y=Y';
>>alpha=;
[b,bint,r,rint,stats]=regress(y,X,alpha)
b=
2、对于寻求模型改进方向的方法通常也可采用残差分析法。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 健康 指数 模型 完整版
![提示](https://static.bdocx.com/images/bang_tan.gif)