R語(yǔ)言進(jìn)行方差分析

1、R語(yǔ)言進(jìn)行方差分析方差分析在試驗(yàn)科學(xué)中有重要的地位，今天談?wù)勅绾斡肦做方差分析前提假設(shè)：獨(dú)立、正態(tài)、方差齊次（各水平間）例子：xbb得到結(jié)果如下（gl指定因子，5是水平，4是重復(fù)次數(shù)）x a1 12 13 14 15 26 27 28 29 310 311 312 313 414 415 416 417 518 519 520 5在進(jìn)行方差分析之前先對(duì)幾條假設(shè)進(jìn)行檢驗(yàn)，由于隨機(jī)抽取，假設(shè)總體滿足獨(dú)立、正態(tài)，考察方差齊次性（用bartlett檢驗(yàn)） (xa,data=b)Bartlett test of homogeneity of variancesdata: x by aBartletts

2、K-squared = , df = 4, p-value =方差齊次性符合下面進(jìn)行方差分析m1summary(m1)Df Sum Sq Mean Sq F value Pr(F)a 4 *Residuals 15Signif. codes: 0 * * * . 1從這個(gè)結(jié)果看出差別顯著接下來(lái)考察具體的差異（多重比較）， TukeyHSD(m1)Tukey multiple comparisons of means95% family-wise confidence levelFit: aov(formula = x a, data = b)$adiff lwr upr p adj2-13-1

3、4-15-13-24-25-24-35-35-4除了5、2和5、4間外，其他之間的差異是不顯著的。方差分析之二具有交互效應(yīng)的雙因素方差分析昨天談了單因素方差分析，今天談?wù)劸哂薪换バ?yīng)的雙因素分析，舉一個(gè)實(shí)例，例子來(lái)源于湯銀才所著R語(yǔ)言與統(tǒng)計(jì)分析。有一個(gè)關(guān)于檢驗(yàn)毒品強(qiáng)弱的試驗(yàn)，給48只白鼠注射3種毒藥（因素A），同時(shí)有4種治療方案（因素B），這樣的試驗(yàn)在每一種因素組合下都重復(fù)四次測(cè)試?yán)鲜蟮拇婊顣r(shí)間，試分析毒藥和治療方案以及它們的交互作用對(duì)老鼠存活時(shí)間有無(wú)顯著影響。數(shù)據(jù)如下time,toxcurebb結(jié)果如下time tox cure1 1 12 1 13 1 14 1 15 1 26 1 27

4、1 28 1 29 1 310 1 311 1 312 1 313 1 414 1 415 1 416 1 417 2 118 2 119 2 120 2 121 2 222 2 223 2 224 2 225 2 326 2 327 2 328 2 329 2 430 2 431 2 432 2 433 3 134 3 135 3 136 3 137 3 238 3 239 3 240 3 241 3 342 3 343 3 344 3 345 3 446 3 447 3 448 3 4方差分析m1summary(m1)Df Sum Sq Mean Sq F value Pr(F) tox 2 *cure 3 *tox:cure 6Residuals 36Signif. codes: 0 * * * . 1由此可以看出tox,cure之間的交互作用不顯著，在說(shuō)明tox,cure的作用是否顯著前，最好先考慮方差齊次性（因?yàn)閺南湫蛨D可以看出方差有很大差異）(timetox+cure,data=b)Bartlett test of homogeneity of variancesdata: time by tox by cureBartlet