在数据框中打印循环的结果

Question

我编写了一个代码，根据因子a拆分数据帧data，并且对于因子的每个级别，都会返回因子b的方差分析表。

for (i in 1:length(levels(data$a))){
  levels<-levels(data$a)
  assign(paste("data_", levels[i], sep = ""), subset(data, a==levels[i]))
  print (levels[i])
  print(anova(lm(var~b, subset(data, a==levels[i]))))
}

结果正是我想要的，但是我希望将所有的anova表汇总起来并作为唯一列表或数据框返回。

任何人都可以帮忙吗？

Answer 1

显然，此代码可以解决问题：

result_anova<-data.frame()
for (i in 1:length(levels(data$a))){
  levels<-levels(data$a)
  assign(paste("data_", levels[i], sep = ""), subset(data, a==levels[i]))
  result<-as.data.frame(anova(lm(var~b, subset(data, a==levels[i]))))
  result_anova[i, 1]<-levels[i]
  result_anova[i, 2]<-result[1, 1 ]
  result_anova[i, 3]<-result[1, 2 ]
  result_anova[i, 4]<-result[1, 3 ]
  result_anova[i, 5]<-result[1, 4 ]
  result_anova[i, 6]<-result[1, 5 ]
  result_anova[i, 7]<-result[2, 1 ]
  result_anova[i, 8]<-result[2, 2 ]
  result_anova[i, 9]<-result[2, 3 ]
  result_anova[i, 10]<-result[2, 4 ]
  result_anova[i, 11]<-result[2, 5 ]
  colnames(result_anova_genos)<-c ( "genotype", "Df_fac", "Sum_Sq_fac",  "Mean_Sq_fac", "F_value_fac", "Pr(>F)_fac", "Df_res", "Sum_Sq_res",  "Mean_Sq_res", "F_value_res", "Pr(>F)_res")
}

请对此答案投票或让我知道是否可以改进此代码。