小编Bet*_*eth的帖子

我有一个包含 40 个数据列（40 种不同的营养素）的数据表，还有用于绘图数字和因子的附加列。我想自动循环每个列名称并为每个列生成一个线性模型和摘要。数据列从第 10 列开始。

for(i in 10:ncol(df)) {       # for-loop over columns
  mod2<-aov(i~block+tillage*residue+Error(subblock),data=df)
  summary(mod2)
}

目前，这正在产生错误Error in model.frame.default(formula = i ~ subblock, data = df, drop.unused.levels = TRUE) : variable lengths differ (found for 'subblock') Variable lengths are concent，所以我想我循环不正确。

数据如下所示（开头有更多分类列），营养成分列从第 10 列开始。