Cox模型危险比置信区间

Question

我使用MICE估算数据并使用生存来运行Cox模型。

输出是：

(est - 1.96 X se, est + 1.96 x se). 

我理解SE与est有关，所以我知道我会做CI的估算：

(exp(est) -1.96 X se, exp(est) + 1.96 X se)

（这是lo 95和hi 95列）。我可以做什么具体的计算来获得HR的CI（exp（est））？

我可以这样做：

{{1}}获取人力资源的CI而非est？

更明显的是，如果我想使用此表发布结果（HR = X，95％CI = X，P Val = X），我该如何从上表中进行处理？

由于

Answer 1

要获得风险比的CI，您应该取幂 limit ，而不是基值，即

exp(c(est - 1.96*se, est + 1.96* se))

或

 exp(est+c(-1.96,1.96)*se)

（对于单个CI;如果est和se是向量，则后一种方法将不工作）

对于输出表，例如：

with(output,
   data.frame(HR=exp(est),lwr=exp(est-1.96*se),upr=exp(est+1.96*se),
              pval=`Pr(>|t|)`))

（我不是100％肯定你想要的。）注意这些是正常（不是基于学生t）的置信区间，但是如果df> 100，它只会产生很小的差异。