我正在使用coxme
来修复混合效应 Cox 模型。我在模型中有两个随机效应(截距 [基线风险] 和处理斜率)。运行模型后,我得到了固定系数和随机效应输出:
Fixed coefficients
coef exp(coef) se(coef) z p
treatment 0.420463710 1.5226675 0.28313885 1.49 0.14
Random effects
Group Variable Std Dev Variance Corr
Study Intercept 1.3979618 1.9542971 -0.6995296
treatment 0.4815433 0.2318839
我还使用ranef
来获取每项研究的随机效应系数。输出:
Intercept treatment
Study1 -1.591497 0.5479255
Study2 1.276003 -0.4392570
Study3 -1.051512 0.3621938
Study4 1.367005 -0.4708623
我想知道如何分别获得每项研究的治疗效果(系数)和 95% CI。我可以通过总结固定效应和研究特定随机效应的总体系数来获得点估计值(例如Study1 的0.420463710
+ 0.5479255
)?但是 95%CI 呢?
有什么方法可以获得研究特定的治疗效果和 95% CI(以及研究的权重),就像ipdforest
在 Stata 中一样?
非常感谢。