R语言回归模型中的OR/HR及95%CI的获得
2018-10-20 本文已影响668人
drlee_fc74
- 在R语言中我们可以通过glm及coxph计算logistic regression和cox regression。但是R语言并没有提供简单的输出其OD/HR及95%值得方法。需要进一步转换
- 通过计算或者其结果
library(survival)
data(pbc)
## Fit a Cox regression model
objCoxph <- coxph(formula = Surv(time, status == 2) ~ trt + age + albumin + ascites, data = pbc)
##得到HR值
exp(coef(model1)
##得到95%CI
exp(confint(model1)
##如果需要专门的一项的话则可以
exp(coef(objCoxph)["trt2"])
exp(confint(objCoxph, parm = "trt2"))
- 使用tableone包
## Load
library(tableone)
## Load Mayo Clinic Primary Biliary Cirrhosis Data
library(survival)
data(pbc)
## Check variables
head(pbc)
## Fit a Cox regression model
objCoxph <- coxph(formula = Surv(time, status == 2) ~ trt + age + albumin + ascites,
data = pbc)
## Show the simple table
ShowRegTable(objCoxph)
###chang digits
ShowRegTable(objCoxph, digits = 3)
