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
1194 PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] 긍정적인마음가짐 2023.05.06
1193 첫번째 행을 라벨로 사용 버튼 클릭시 에러 문의 [3] file 수담이 2023.05.04
1192 안녕하세요? web-R table 그리고 관련하여 문의드립니다. [4] mkw7654 2023.05.02
1191 메타분석 U3 그래프 관련 문의드립니다 secret 화이팅 2023.05.02
1190 교수님 안녕하세요~~ secret 메시12 2023.04.29
1189 연속형 변수인데 범주형으로 인지됩니다. [2] secret shararaku 2023.04.28
1188 통계마당 [2] secret ejo 2023.04.27
1187 moonBook 오류 관련 재문의드립니다. [1] secret holyhyun 2023.04.24
1186 moonBook 사용시 오류메시지로 문의드립니다. [1] holyhyun 2023.04.24
1185 에러 문의 [2] secret 화이팅 2023.04.19
1184 Cox regression 에 대한 문의입니다. [2] file saga 2023.04.18
1183 PSM 문의 [3] secret starwars 2023.04.13
1182 ggPS, ggPSMSummary function [6] 로렌쵸메디치 2023.04.12
1181 Survival 분석에 대해서 문의 드립니다. 정훈신 2023.04.05
1180 안녕하세요. 분석 툴에 대해 문의드립니다. [1] 스매직맨 2023.04.04
1179 드디어 publish 되었네요. [2] 정지윤 2023.04.04
1178 연구동향분석 시 논문 데이터 정제 관련 [5] vtlbtl 2023.03.31
1177 안녕하세요 교수님 에러 문의드립니다 [2] secret 소소 2023.03.26
1176 코딩값이 여러개인 변수에 대한 전체 p값을 알 수 있는 방법이 궁금합니다 쥬쥬즈주 2023.03.22
1175 PSM 매칭 버튼 보이지 않음. secret swpapa 2023.03.19
Write Tag