试验设计两因子无交互作用方差分析的r语言代码

shuju=read.table("C:/Users/lenovo/Desktop/shuju.txt",header=T)

shuju  #读取数据集

y=c(shuju[ ,1],shuju[ ,2],shuju[ ,3]);y      #拉长数据，使之成为一列

g1=rep(c(1:4),3);g1     #分组的列

g2=rep(c(1:3),rep(4,3));g2     #另外一个分组的列

data=cbind(g1,g2,y)

data   #将数据造成了用于做实验分析的形式。也可以一开始就将数据写成这种格式

f1=factor(g1);mode(g1)

f2=factor(g2);mode(g2)  #将g1和g2的类型变为factor分类的

obj=lm(y~f1+f2)     #做回归

anova(obj)      #做回归分析中的方差分析

obj1=lm(y~g1+g2);anova(obj1)  #要是不将g1和g2变换的话，方差分析的自由度错误

#下面三行是另外一种在多元回归分析中使用的方差分析方法，先安装包，但是发现结果一样

install.packages("car")  #请先安装car包，下面的Anova函数在car包中

library(car)

Anova(obj)      #用Anova函数

#下面两行是有交互作用的方差分析

objnew=lm(y~f1+f2+f1*f2)
anova(objnew) 

附：数据原样

      B1    B2   B3 

A1 58.2 56.2 65.3

A2 49.1 54.1 51.6

A3 60.1 70.9 39.2

A4 75.8 58.2 48.7