Skip to menu

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


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을 구하는 데에서 막여 문의 드립니다.


감사합니다.



No. Subject Author Date
1113 ROC 비교 관하여 질문드립니다!(delong's test) 나현둥 2022.08.09
1112 안녕하세요교수님! VIF 관련하여 질문드립니다. file kms9057 2022.08.08
1111 cox 에서 AIC모델 구하는 방법이 있나요? [1] 쥬쥬즈주 2022.08.03
1110 PSM 중 server disconnected [3] joonho 2022.08.01
1109 서로 다른 set 에서의 kappa 값 비교 papahong 2022.07.30
1108 PSM에 대하여 [1] 김습습 2022.07.29
1107 엑셀의 index 기능을 R 에서 구현하려고 합니다 file emilien 2022.07.28
1106 survival 분석 중 PPT 다운로드 오류 secret 곰곰11 2022.07.26
1105 결과 ppt 다운로드 중 오류발생문제 [2] 곰곰11 2022.07.25
1104 Control/experimental Group이 있을 때 유의미한 결과 도출 방법 [3] 민됴 2022.07.25
1103 Cox uni multivaraite [2] zsdfa 2022.07.14
1102 생존분석 문의드립니다 [1] 최홍재 2022.07.13
1101 인용 관련하여 [1] kms9057 2022.07.12
1100 메타 회귀 분석 시 ggplot 출력에 관해 문의드립니다. endostat 2022.07.09
1099 생존곡선에서 색깔 변화 문의드립니다. file 케로스 2022.07.09
1098 후향적 임상시험에서의 검체수 산정 방법 질문 라라보라 2022.07.05
1097 Kaplan Meier survival curve 질문입니다 [2] cetla 2022.07.03
1096 mytalbe( )패키지 사용중 문의드립니다. 하은 2022.07.02
1095 심평원 R studio 내 lubridate 설치 오류 [4] 꾸이뭉 2022.06.27
1094 안녕하세요 r 코드 중에 도저히 해결이 안나는 문제가 있어서 여기에 문의 드립니다 [1] file 슈퍼맨 2022.06.24