1.雙因素方差分析
這里我們使用內(nèi)置的的ToothGrowth數(shù)據(jù)集（60只豚鼠分為兩大組，分別喂養(yǎng)橙汁和維生素C钝荡，每種飼喂方法中抗壞血酸含量有三個水平，分別是0.5mg/day端辱，1mg/day，2mg/day个盆，每組有十只豚鼠颊亮；最終測試的牙齒長度是因變量）來進行演示门扇。如果只是看哪個效果更加明顯可以直接計算平均值悯嗓，如果要看差別的顯著水平則需要用方差分析。supp表示飼喂方法合武，dose表示飼喂的三種水平。
圖1 數(shù)據(jù)集截圖
attach(ToothGrowth)#用來加載數(shù)據(jù)红淡，用attach函數(shù)就不需要用變量賦值操作
table(supp,dose)#統(tǒng)計分組情況，結(jié)果是一個二聯(lián)表
aggregate(len, by=list(supp,dose), FUN=mean)#分組統(tǒng)計平均數(shù)
ToothGrowth$dose<-factor(ToothGrowth$dose)#將該列轉(zhuǎn)化為因子類型谅畅，supp列已經(jīng)是因子類型
fit<-aov(len~supp*dose,data=ToothGrowth)#進行雙因素方差分析
summary(fit)#總結(jié)結(jié)果
圖2 查看分組，均值以及方差分析結(jié)果
結(jié)果顯示芯咧，藥物種類和藥物劑量對牙齒長度數(shù)據(jù)均有明顯影響，P值小于0.05邪铲。
2.繪圖
可以用HH包中的interaction.plot函數(shù)進行繪圖
install.packages("HH")
library(HH)#載入包
interaction.plot(dose,supp,len,type="b",col=c("red","blue"),pch=c(16,18),main="Interaction between Dose and Supplement Type")#繪出雙因素分析圖
圖3 雙因素分析圖
