eff=c(58.2,52.6,56.2,41.2,65.3,60.8,49.1,42.8,54.1,50.5,51.6,48.4,60.1,58.3,70.9,73.2,39.2,40.7,75.8,71.5,58.2,51.0,48.7,41.4)
tuijin=rep(c(1,1,2,2,3,3),4)
ranliao=rep(1:4,rep(6,4))
tuijin=as.factor(tuijin)
ranliao=as.factor(ranliao)
datas=data.frame(eff,tuijin,ranliao)
fit <- aov(eff~tuijin*ranliao,data=datas)
# <=> fit <- aov(eff~ tuijian + ranliao + tuijin:ranliao,data=datas)
summary(fit)
# 查看均值
model.tables(fit,'means')