打印数据框中循环的结果

Question

我编写了一个代码，根据因子 a 拆分数据帧 data，对于因子 returns 的每个级别，因子 b 的方差分析 table。

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]))))
}

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

有人可以帮忙吗？

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")
}

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

打印数据框中循环的结果

Print results of a loop in a dataframe

loops

anova

lm