如何从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/

10-12 13:59