R语言双因素方差分析

Posted 2021-09-06 Data+Science+Insight

R语言双因素方差分析

条件：

各个样本是相互独立的随机；

各个样本来自正态总体；

具有方差齐性；

用途：

检验两个或多样本均数间的差异有无统计学意义；注：本均数的比较可以采用 t检验或 F检验，两个以上样本均数的比较只 能用 F检验。
回归方程的线性假设检验；
检验两个因素或者多个因素之间有无交互作用；

无法进行方差分析时：

进行变量换，以达到方差齐或正态的要求；
采用非参数法（秩和检验）；
使用近似 F检验；

# Listing 9.6 - Two way ANOVA
attach(ToothGrowth)
table(supp,dose)
aggregate(len, by=list(supp,dose), FUN=mean)
aggregate(len, by=list(supp,dose), FUN=sd)
dose <- factor(dose)
fit <- aov(len ~ supp*dose) # fit1<- aov(len~supp+dose+supp:dose)
summary(fit)

# plotting i