从“by”对象中提取表输出（如Summary和ANOVA）

Question

我有一个用by创建的回归模型。我知道我可以使用sapply为每个因素提取模型的特定部分，但如果我想要整个summary，anova等等，会怎样？

 model <- with(data, by(data, factor, function(data) lm(y ~ x, data=data)))

Answer 1

sapply会将summary.lm和anova.lm的结果强制转换为矩阵。我想你可能想要使用lapply，它在summary生成的列表中的每个元素上应用一个函数（此处为by），并返回一个列表。

models <- by(warpbreaks, warpbreaks$tension, function(x){
  lm(breaks ~ wool, data = x)
  })

lapply(models, summary)