如何从“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方法调度的结果。