Loops 在数据帧中打印循环的结果

我编写了一个代码，根据因子a拆分数据帧

数据

，对于因子的每个级别，返回因子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]))))
}

for（1中的i：长度（级别（数据$a）））{
级别显然，此代码实现了以下功能：

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

result\u anova显然，这个代码起到了关键作用：

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

result\u方差分析