r - 如何从 "aov"模型中提取系数的标准误差

Question

我做了一个 aov模型，我只想提取每个系数的标准误差。

model <- aov(Molecule ~  Comorbidity + Age + BMI + Sex, data = mydata)

我可以看到估计的边际平均值:

allEffects(model)

但我无法提取每个系数的标准误差。我在网上看到的功能 se.coef()但它不起作用。 summary(model)$coefficients[, "Std. Error"]也不起作用。

我已阅读 effects和其他包，但我没有找到我想要的。任何的想法？

Answer 1

使用 lm summary 的方法:

coef(summary.lm(model))

这将为所有可识别系数提供 4 列(平均值、标准误差、t 值、p 值)的系数表/矩阵。然后您可以提取标准误差的第二列。
aov返回主要类“aov”但次要类“lm”的对象，因此两者都是 summary.aov和 summary.lm申请但给出不同的东西。当你只是做 summary(model) ，前者作为 S3 方法调度的结果被调用。