資源描述：
《多元統(tǒng)計(jì)分析方法 covariance 協(xié)方差分析》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、協(xié)方差分析一個(gè)協(xié)變量的協(xié)方差分析例：為研究三種飼料(A1(g=1)，A2(g=2)，A3(g=3))對(duì)豬催肥效果，用每種飼料喂養(yǎng)8頭豬，實(shí)驗(yàn)用豬的初始體重(x)未控制。喂養(yǎng)一段時(shí)間后，觀察小豬的增重(y)。所得資料如表2－1，試分析三種飼料對(duì)豬催肥效果是否相同。資料結(jié)構(gòu)：(文件名covariance1.dta)xyg1585113831116511276112801169111484117901179721690218100218952211032221062199921894222893249132083323953251003271023301053321103對(duì)于不考慮初始體重影響而
2、評(píng)價(jià)三種飼料的統(tǒng)計(jì)分析為單因素方差分析(One-wayANOVA)，由于小豬的增重與初始體重有關(guān)，因此在分析三種飼料對(duì)增重的關(guān)系時(shí)，應(yīng)該考慮校正初始體重對(duì)增重的影響。并假定初始體重與增重呈線性統(tǒng)計(jì)關(guān)系以及要求初始體重與飼料不構(gòu)成交互作用。稱校正變量(初始體重)為協(xié)變量，分組變量為因子變量。因此可用協(xié)方差分析上述統(tǒng)計(jì)問(wèn)題，相應(yīng)的角模型如下：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項(xiàng)的P值＝0.1184>0.05，說(shuō)明初始體重與飼料不構(gòu)成交互作用。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:a110對(duì)應(yīng)的P值為0.001<0.05，因此認(rèn)為兩組總體均數(shù)不同，由a1的95％可信區(qū)間可認(rèn)為A1的均數(shù)大于A3的均數(shù)，差別有統(tǒng)計(jì)意義。A2vsA3：，H0：a2=0vsH1:a210對(duì)應(yīng)的P值為0.001<0.05，因此認(rèn)為兩組總體均數(shù)不同，由a2的95％
28、可信區(qū)間可認(rèn)為A2的均數(shù)大于A3的均數(shù)，差別有統(tǒng)計(jì)意義。A1vsA2：，H0：a1-a2=0vsH1:a1-a210test_b[g[1]]-_b[g[2]]=0(1)g[1]-g[2]=0.0F(1,20)=4.70Prob>F=0.0424對(duì)應(yīng)的P值為0.0424<0.05，因此認(rèn)為兩組總體均數(shù)不同，由于點(diǎn)估計(jì)為：<，P值小于0.05，因此可認(rèn)為A2的均數(shù)大于A1的均數(shù)，差別有統(tǒng)計(jì)意義。結(jié)論：1)A2飼料喂養(yǎng)的小豬增
