Я использую coxme, чтобы исправить модель Кокса со смешанным эффектом. У меня есть два случайных эффекта в модели (пересечение [базовая опасность] и уклон лечения). После запуска модели я получил фиксированные коэффициенты и выходные данные случайных эффектов:
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% ДИ для каждого исследования в отдельности. Можно ли получить точечную оценку путем суммирования общего коэффициента фиксированного эффекта и случайного эффекта для конкретного исследования (например, 0.420463710 + 0.5479255 для исследования 1)? Но как насчет 95% ДИ?
Есть ли способ получить эффект лечения для конкретного исследования и 95% ДИ (и вес исследования), как в ipdforest в Stata?
Большое спасибо.