항상 빠른 답변 감사드립니다.
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을 구하는 데에서 막여 문의 드립니다.
감사합니다.
Comment 0
No. | Subject | Author | Date |
---|---|---|---|
1114 | 성향점수 맞추기 [3] | alse | 2016.05.25 |
1113 | 아래에 올렸었는데요,, 5Y overall survival 구하는 법이요. [2] | 난둥 | 2016.05.26 |
1112 | 로지스틱 회귀 분석에서.. [1] | 최00 | 2016.05.27 |
1111 | 기본표를 만드는데 file이 이상하게 읽히는 것 같습니다. | 매일이5월 | 2016.06.02 |
1110 | 파일 불러오기 에러 [2] | 매일이5월 | 2016.06.02 |
1109 | 파일 불러오기 오류 [2] | 매일이5월 | 2016.06.02 |
1108 | RPubs에 있는 ROC curve analysis by cardiomoon [2] | hawk1227 | 2016.06.02 |
1107 | R 그래프 옮기기 [2] | hawk1227 | 2016.06.03 |
1106 | 생존분석 그래프 x축 시간 [3] | alse | 2016.06.07 |
1105 | 자료업로드에러 [3] | dejavu730 | 2016.06.07 |
1104 | 학생강의 project1 | cardiomoon | 2016.06.07 |
1103 | 결측치 처리에 대한 질문입니다. [2] | 매일이5월 | 2016.06.07 |
1102 | IBM의 R 컨소시엄 가입 | Ben | 2016.06.08 |
1101 | 서버 접속이 안됩니다 [2] | kim | 2016.06.08 |
1100 | binary data 군집분석 [2] | 노가리 | 2016.06.08 |
1099 | 생존 분석시 문제 [3] | selteen | 2016.06.09 |
1098 | 생존 분석을 하는데 질문이 있어 게시판에 글을 올립니다. [1] | namuhanayeyo | 2016.06.11 |
1097 | 독립(설명)변수 비정규 분포시. logistic regression 문의입니다. [1] | 매일이5월 | 2016.06.16 |
1096 | Simple linear regression 에서 설명변수의 특성 [3] | 매일이5월 | 2016.06.16 |
1095 | Rstudio 에러 문제 [2] | 화모기 | 2016.06.16 |