如何将Z - statistics的值作为来自glm对象的向量？例如,我有
fit <- glm(y ~ 0 + x,binomial)
如何Pr(>|z|)以与系数估计相同的方式访问列fit$coef？
Ben*_*ker 11
我相信
coef(summary(fit))[,"Pr(>|z|)"]
会得到你想要的.(summary.glm()返回一个具有coef()返回系数表的方法的对象.)(顺便说一句,如果存在访问器方法,最好使用它们而不是直接访问拟合模型的组件 - 例如coef(fit),优于fit$coef.)
从线性回归中拉出p值和r平方给出了类似的答案.
我建议methods(class="summary.glm")找到可用的访问器方法,但它实际上有点棘手,因为默认方法(在这种情况下coef.default())也可能是相关的...
PS如果你想要Z值,coef(summary(fit))[,"z value"]应该这样做(你的问题有点含糊不清:通常当人们说"Z统计"时他们的意思是想要测试统计的值,而不是p值)