如何从glm
对象获取Z-统计值作为矢量?
例如,我有
fit <- glm(y ~ 0 + x,binomial)
如何使用
Pr(>|z|)
获得系数估计值的方式与访问fit$coef
列相同? 最佳答案
我相信
coef(summary(fit))[,"Pr(>|z|)"]
会得到你想要的东西。 (
summary.glm()
返回一个对象,该对象具有一个返回系数表的coef()
方法。)(顺便说一句,如果存在访问器方法,则比直接访问拟合模型的组件要好用它们,例如coef(fit)
比fit$coef
更好。 )pull out p-values and r-squared from a linear regression给出了类似的答案。
我建议使用
methods(class="summary.glm")
来找到可用的访问器方法,但实际上要比这复杂一点,因为默认方法(在这种情况下为coef.default()
)也可能是相关的...PS 如果您想要Z值,
coef(summary(fit))[,"z value"]
应该这样做(您的问题有点模棱两可:通常,当人们说“Z statistic”时,他们的意思是想要测试统计量的值,而不是p值)关于r - 如何从glm对象获取Z值-统计信息?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/12341338/