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

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

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

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

allEffects(model)

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

我已经阅读了effects和其他软件包，但是找不到我想要的东西。任何想法？

最佳答案

对lm使用summary方法：

coef(summary.lm(model))

这将为所有可识别系数提供4列的系数表/矩阵（平均值，标准误差，t值，p值）。然后，您可以提取第二列以获取标准误差。

aov返回主类“ aov”的对象，但返回辅助类“ lm”的对象，因此，summary.aov和summary.lm都适用，但给出了不同的内容。当您简单地执行summary(model)时，前者将作为S3方法分派的结果被调用。