资源描述:
《医学统计学(扬州大学)协方差分析》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、协方差分析一个协变量的协方差分析例:为研究三种饲料(A1(萨1),A2(g=2),A3(g=3))对猪催肥效果,用每种饲料喂养8头猪,实验用猪的初始体重(x)未控制。喂养一段时间后,观察小猪的增重(y)。所得资料如表2—1,试分析三种饲料对猪催肥效果是否相同。资料结构:(文件名covariancel.dta)Xyg158511383111651127611280116911148411790117972169021810021895221103222106219992189422289324913208332395325100327102330105
2、3321103对于不考虑初始体重影响而评价三种饲料的统计分析为单因素方茅协方差分析一个协变量的协方差分析例:为研究三种饲料(A1(萨1),A2(g=2),A3(g=3))对猪催肥效果,用每种饲料喂养8头猪,实验用猪的初始体重(x)未控制。喂养一段时间后,观察小猪的增重(y)。所得资料如表2—1,试分析三种饲料对猪催肥效果是否相同。资料结构:(文件名covariancel.dta)Xyg158511383111651127611280116911148411790117972169021810021895221103222106219992189422
3、2893249132083323953251003271023301053321103对于不考虑初始体重影响而评价三种饲料的统计分析为单因素方茅分析(One-wayANOVA),由于小猪的增重与初始体重有关,因此在分析三种饲料对增重的关系时,应该考虑校正初始体重对增重的影响。并假定初始体重与增重呈线性统计关系以及要求初始体重与饲料不构成交互作用。称校正变量(初始体重)为协变量,分组变量为因子变量。因此可用协方差分析上述统计问题,相应的角模型如下:Al(g=l)A2(g=2)A3(g=3)不校正初始体重校正初始体重“..+0+"p+yx用STATA命令
4、为:anovaygxg*x,class(g)Numberofobs=24R-squared=0.9297RootMSE=3.15855AdjR-squared=0.9102SourcePartialSSdfMSFProb>FModel2376.38195475.2763847.640.0000g24.4661579212.2330791.230.3168X830.4154071830.41540783.240.0000g*x48.0381359224.0190682.410.1184Residual179.576433189.97646848Tota
5、l2555.9583323111.128623由g*x项的卩值=0.1184>0.05,说明初始体重与饲料不构成交互作用。anovaygx,class(g)Numberofobs=24R-squared37353AdjR-squared=0.9109二0.8976RootMSE=3.SourcePartialSSdfMSFProb>FModel2328.343763776.1145886&200.0000g707.2187652353.60938231.070.0000X1010.7604311010.7604388.810.0000Residual
6、227.6145682011.3807284Total2555.9583323111.128623regressSource
7、SSdfMSNumberofobs=24+F(3,20)68.20Model
8、2328.343763776.114588Prob>F0.0000Residual227.6145682011.3807284R-squared0.9109AdjR-squared0.8976tTotal
9、2555.9583323111.128623RootMSE3.3735yCoef.Std.Err.tP>
10、t
11、[95%Conf.Interval]
12、_cons35.935186.5754715.470.00022.2189949.65137g112.793243.4089893.750.0015.68221419.90427217.335592.4091517.200.00012.3101922.360993(dropped)X2.401569・25483329.420.0001.8699962.93314235.93518,%=12.79324,a2=:17.33559,?=2.401569AlvsA3:(〃+e+/X)一(“+")=e,Ho:oti=0vsHi:cX]HO对应的P值为0.00
13、1<0.05,因此认为两组总体均数不同,由8的95%可信区间可认为A1的均数大于A3的均数,差别有统计意义。