双因素方差分析

谷小燕老师

#双因素分析部分数据如下

#有交互作用（*）
  mudle<-read.csv("double_ANOVA.csv",header = FALSE)
  #因素要是字符串型
  mudle[,1:2]<-lapply(mudle[,1:2],as.character)
  #因素数据统计
  table(mudle$V1,mudle$V2)
  #aov 双因素分析
  mudle_aov<-aov(mudle$V3~mudle$V1 * mudle$V2,data=mudle)
  summary(mudle_aov)

  #多重比较（Tukey氏固定差距检验法）
  TukeyHSD(mudle_aov)

＃无交互作用（＋）
　mudle_aov<-aov(mudle$V3~mudle$V1 + mudle$V2,data=mudle)
  summary(mudle_aov)
#多重比较（Tukey氏固定差距检验法）
  TukeyHSD(mudle_aov)

谷小燕老师

  #方差齐
  bartlett.test(V3~V1,data=mudle)
  bartlett.test(V3~V2,data=mudle)
  #正态
  aggregate(mudle$V3,by=list(mudle$V2),FUN=shapiro.test)
  aggregate(mudle$V3,by=list(mudle$V1),FUN=shapiro.test)