Skip to menu

Cox proportional analysis에서 subgroup분석에 대해 여쭙니다

김규 2017.09.30 10:04 Views : 12842

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


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


감사합니다.



I want to
List
No. Subject Author Date
1214 안녕하십니까 교수님 PSM 결과에 대해 궁금한 것이 있어 여쭤봅니다 [2] secret 김도훈 2023.06.29
1213 교수님 지난 번 문의 드렸던 PSM 재문의 드립니다. [3] secret 긍정적인마음가짐 2023.06.27
1212 안녕하세요 ^^ odds ratio 와 관련된 문의드리고자 합니다 ^^ [2] 성열 2023.06.26
1211 mytable() 함수에서 na처리와 pvalue 계산 문의드립니다. [1] hjk34 2023.06.23
1210 교수님 안녕하세요 ^^ web-R table 작성관련 문의드립니다. [2] file mkw7654 2023.06.21
1209 Methods for continuous variables과 P value corrected for multiple comparisons 방법 궁금합니다. [1] 호호셀 2023.06.19
1208 ordinal logistic regression analysis 문의 [1] 버블 2023.06.18
1207 e-mail 변경 문의 [1] secret kimdo 2023.06.16
1206 안녕하세요 number at risk 관련 문의 드립니다 [1] secret cebukim 2023.06.13
1205 Propensity matching score 에 대해 문의드립니다. [1] saga 2023.06.08
1204 콕스비례위험모델의 적용에 관련해서 질문 드립니다. [2] secret kng1032 2023.06.07
1203 기초통계 계산 [2] secret leesav 2023.06.07
1202 정회원 후원하였습니다. [2] secret toast 2023.06.04
1201 logistic regression analysis 가 에러가 뜹니다. [1] secret Ikarus86 2023.06.02
1200 PMS caliper 에 대해 여쭤봅니다. [2] secret pjhaaa 2023.05.30
1199 ggGam에 대하여 [2] leesav 2023.05.22
1198 R 에서 PFS 의 confidence interval 을 구하고 싶을 때 어떻게 하면 될까요? [1] ddeeny 2023.05.21
1197 메타분석 관련 문의드립니다 secret 화이팅 2023.05.12
1196 교수님, 새로운 앱 개발 문의드립니다! [4] swpapa 2023.05.08
1195 교수님, 안녕하세요? 전용 앱 관련 문의드립니다! [1] swpapa 2023.05.07
Write Tag