资源描述：
《多元统计分析方法 covariance 协方差分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、协方差分析一个协变量的协方差分析例：为研究三种饲料(A1(g=1)，A2(g=2)，A3(g=3))对猪催肥效果，用每种饲料喂养8头猪，实验用猪的初始体重(x)未控制。喂养一段时间后，观察小猪的增重(y)。所得资料如表2－1，试分析三种饲料对猪催肥效果是否相同。资料结构：(文件名covariance1.dta)xyg1585113831116511276112801169111484117901179721690218100218952211032221062199921894222893249132083323953251003271023301053321103对于不考虑初始体重影响而
2、评价三种饲料的统计分析为单因素方差分析(One-wayANOVA)，由于小猪的增重与初始体重有关，因此在分析三种饲料对增重的关系时，应该考虑校正初始体重对增重的影响。并假定初始体重与增重呈线性统计关系以及要求初始体重与饲料不构成交互作用。称校正变量(初始体重)为协变量，分组变量为因子变量。因此可用协方差分析上述统计问题，相应的角模型如下：A1(g=1)A2(g=2)A3(g=3)不校正初始体重校正初始体重用STATA命令为：anovaygxg*x,class(g)Numberofobs=24R-squared=0.9297RootMSE=3.15855AdjR-squared=0.9102
3、Source
4、PartialSSdfMSFProb>FModel
5、2376.38195475.2763847.640.0000
6、g
7、24.4661579212.2330791.230.3168x
8、830.4154071830.41540783.240.0000g*x
9、48.0381359224.0190682.410.1184
10、Residual
11、179.576433189.97646848Total
12、2555.9583323111.128623由g*x项的P值＝0.1184>0.05，说明初始体重与饲料不构成交互作用。anovaygx,class(g)Numberofobs=24R-squ
13、ared=0.9109RootMSE=3.37353AdjR-squared=0.8976Source
14、PartialSSdfMSFProb>FModel
15、2328.343763776.11458868.200.0000g
16、707.2187652353.60938231.070.0000x
17、1010.7604311010.7604388.810.0000Residual
18、227.6145682011.3807284Total
19、2555.9583323111.128623regressSource
20、SSdfMSNumberofobs=24-------------+-------------
21、-----------------F(3,20)=68.20Model
22、2328.343763776.114588Prob>F=0.0000Residual
23、227.6145682011.3807284R-squared=0.9109-------------+------------------------------AdjR-squared=0.8976Total
24、2555.9583323111.128623RootMSE=3.3735yCoef.Std.Err.tP>
25、t
26、[95%Conf.Interval]_cons35.935186.5754715.470.00022.21899
27、49.65137g112.793243.4089893.750.0015.68221419.90427217.335592.4091517.200.00012.3101922.360993(dropped)x2.401569.25483329.420.0001.8699962.933142，，，A1vsA3：，H0：a1=0vsH1:a1¹0对应的P值为0.001<0.05，因此认为两组总体均数不同，由a1的95％可信区间可认为A1的均数大于A3的均数，差别有统计意义。A2vsA3：，H0：a2=0vsH1:a2¹0对应的P值为0.001<0.05，因此认为两组总体均数不同，由a2的95％
28、可信区间可认为A2的均数大于A3的均数，差别有统计意义。A1vsA2：，H0：a1-a2=0vsH1:a1-a2¹0test_b[g[1]]-_b[g[2]]=0(1)g[1]-g[2]=0.0F(1,20)=4.70Prob>F=0.0424对应的P值为0.0424<0.05，因此认为两组总体均数不同，由于点估计为：<，P值小于0.05，因此可认为A2的均数大于A1的均数，差别有统计意义。结论：1)A2饲料喂养的小猪增
