如何从 "aov" 模型中提取系数的标准误差
How to extract coefficients' standard error from an "aov" model
我做了一个aov模型，我只想提取每个系数的标准误差。
model <- aov(Molecule ~ Comorbidity + Age + BMI + Sex, data = mydata)
我可以看到估计的边际均值：
allEffects(model)
但我无法提取每个系数的标准误差。我在互联网上看到 se.coef() 函数，但它不起作用。 summary(model)$coefficients[, "Std. Error"] 也不行。
我已经阅读了 effects 和其他软件包，但没有找到我想要的。有什么想法吗？
对summary使用lm方法：
coef(summary.lm(model))
这将为所有可识别系数提供系数 table / 4 列矩阵（平均值、标准误差、t 值、p 值）。然后你可以提取标准错误的第二列。
aov returns 主要对象 class "aov" 但次要对象 class "lm"，因此 summary.aov 和 summary.lm 适用但给出不同的东西。当你简单地做 summary(model) 时，前者作为 S3 方法调度的结果被调用。
我做了一个aov模型，我只想提取每个系数的标准误差。
model <- aov(Molecule ~ Comorbidity + Age + BMI + Sex, data = mydata)
我可以看到估计的边际均值：
allEffects(model)
但我无法提取每个系数的标准误差。我在互联网上看到 se.coef() 函数，但它不起作用。 summary(model)$coefficients[, "Std. Error"] 也不行。
我已经阅读了 effects 和其他软件包，但没有找到我想要的。有什么想法吗？
对summary使用lm方法：
coef(summary.lm(model))
这将为所有可识别系数提供系数 table / 4 列矩阵（平均值、标准误差、t 值、p 值）。然后你可以提取标准错误的第二列。
aov returns 主要对象 class "aov" 但次要对象 class "lm"，因此 summary.aov 和 summary.lm 适用但给出不同的东西。当你简单地做 summary(model) 时，前者作为 S3 方法调度的结果被调用。