항상 빠른 답변 감사드립니다.

R에서 subgroup analysis를 통한 cox분석을 하려 하는데 p-value for interaction을 구하는 방법을 잘 몰라서 부끄럽지만 여쭙니다.

Statin이 심혈관계 합병증에 미치는 영향에 대하여 분석 중 DM에 대한 subgroup 을 다음과 같이 구하였습니다.

> COXMACE=coxph(Surv(MACE.day,MACCE==1)~Statin.use,data=ESP,subset=DM==0)

> summary(COXMACE)

Call:

coxph(formula = Surv(MACE.day, MACCE == 1) ~ Statin.use, data = ESP, 

    subset = DM == 0)

  n= 747, number of events= 73 

              coef exp(coef) se(coef)      z Pr(>|z|)    

Statin.use -0.8978    0.4075   0.2627 -3.417 0.000632 ***

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

           exp(coef) exp(-coef) lower .95 upper .95

Statin.use    0.4075      2.454    0.2435    0.6819

Concordance= 0.612  (se = 0.033 )

Rsquare= 0.017   (max possible= 0.676 )

Likelihood ratio test= 13.01  on 1 df,   p=0.0003097

Wald test            = 11.68  on 1 df,   p=0.0006325

Score (logrank) test = 12.48  on 1 df,   p=0.0004118

-------------------------------------------------------------------------------------------------------------------

> COXMACE=coxph(Surv(MACE.day,MACCE==1)~Statin.use,data=ESP,subset=DM==1)

> summary(COXMACE)

Call:

coxph(formula = Surv(MACE.day, MACCE == 1) ~ Statin.use, data = ESP, 

    subset = DM == 1)

  n= 333, number of events= 50 

              coef exp(coef) se(coef)      z Pr(>|z|)  

Statin.use -0.6009    0.5483   0.2888 -2.081   0.0375 *

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

           exp(coef) exp(-coef) lower .95 upper .95

Statin.use    0.5483      1.824    0.3113    0.9657

Concordance= 0.585  (se = 0.039 )

Rsquare= 0.013   (max possible= 0.786 )

Likelihood ratio test= 4.44  on 1 df,   p=0.03509

Wald test            = 4.33  on 1 df,   p=0.03746

Score (logrank) test = 4.46  on 1 df,   p=0.03467

이후에 DM이 있고 없고에 따른 statin효과의 차이를 보고 싶은데 p-value interaction을 구하는 데에서 막여 문의 드립니다.

감사합니다.