r - R从cbind中的glm返回p值

Question

Statistics 和 R noob 想知道是否有办法将 glm 中的 p 值添加到由以下命令产生的输出的末尾：

exp(cbind(OR = coef(mod1), confint(mod1)))

也许是这样的：

summary(mod1)$coefficients[,4]

我意识到这在某种程度上是一个“化妆品”问题，但它仍然很方便。

谢谢

Answer 1

您可以保存 summary(mod1) 的结果，然后使用 访问系数表coefficients。

您可以编写一个函数来为您完成整个过程......

OR.summary <- function(x){
 # get the summary
  xs <- summary(x)
 # and the confidence intervals for the coefficients 
  ci = confint(x)
 # the table from the summary object
  coefTable <- coefficients(xs)
 # replace the Standard error / test statistic columns with the CI
  coefTable[,2:3] <- ci
 # rename appropriatly
  colnames(coefTable)[2:3] <- colnames(ci)
# exponentiate the appropriate columns
  coefTable[,1:3] <- exp(coefTable[,1:3])
# return the whole table....
  coefTable

}

一种更强大的方法是使用像rms....这样的包。