首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >R在cbind中从glm返回p值

R在cbind中从glm返回p值
EN

Stack Overflow用户
提问于 2014-03-05 01:40:08
回答 1查看 823关注 0票数 0

统计数据和R想知道是否有一种方法可以将p值从glm添加到由以下命令产生的输出的末尾:

代码语言:javascript
运行
复制
exp(cbind(OR = coef(mod1), confint(mod1)))

也许是这样的:

代码语言:javascript
运行
复制
summary(mod1)$coefficients[,4]

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

谢谢

EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2014-03-05 02:02:40

您可以保存摘要结果(Mod1),然后使用coefficients访问系数表。

你可以写一个函数为你完成整个过程.

代码语言:javascript
运行
复制
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这样的软件包.

票数 2
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/22186817

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档