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

【问题标题】：Print results of a loop in a dataframe在数据框中打印循环的结果
【发布时间】：2020-03-22 04:14:32
【问题描述】：

我编写了一个代码，它根据因子 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]))))
}

结果正是我想要的，但我希望将所有方差分析表合并并作为唯一列表或数据框返回。

有人可以帮忙吗？

【问题讨论】：

标签： loops lm anova

【解决方案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")
}

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

【讨论】：