항상 빠른 답변 감사드립니다.
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 |
---|---|---|---|
1253 | R서버 설치 문의드립니다. [4] | 화모기 | 2015.05.13 |
1252 | 초초 초보자 문의 드립니다^^ [2] | 매직하마 | 2017.12.22 |
1251 | mytable() 함수에 대한 질문입니다. [3] | 브리디온 | 2017.11.19 |
1250 | 묻고답하기 게시판을 운영합니다 [1] | cardiomoon | 2015.04.20 |
1249 | ANOVA, Kruskal-Wallis 검정 및 다중비교 [3] | 소오름 | 2016.03.31 |
1248 | ios | bnmy6581 | 2018.12.28 |
1247 | 예제 데이터 열 제목 편집방법 [2] | 신돌석 | 2016.04.25 |
1246 | rstudio 사이트에 연결할 수 없음 문제 [1] | 지원 | 2019.09.10 |
1245 | 로지스틱 회귀분석에 대해 질문 올립니다. [12] | SOL | 2017.09.15 |
» | Cox proportional analysis에서 subgroup분석에 대해 여쭙니다 | 김규 | 2017.09.30 |
1243 | rJava 에러 문제 해결에 대해 문의드립니다 [1] | mdhjcho | 2017.06.29 |
1242 | 3D 형태 그래프에 대해 질문드립니다. | 백지백지 | 2019.04.26 |
1241 | 메타분석의 변수지정 [1] | 야스오 | 2021.05.24 |
1240 | Logistic regression [3] | namuhanayeyo | 2017.03.02 |
1239 | 생존 분석 시 Cox-regression vs. Log-rank [2] | begabung | 2017.05.22 |
1238 | 표본설계시 subset 명령어 문의드립니다. [1] | 네잎클로바 | 2020.05.09 |
1237 | Rstudio 사용법` [3] | 방랑자 | 2015.05.02 |
1236 | Survdiff로 log-rank test 시 P-value 소수점 [3] | begabung | 2019.02.15 |
1235 | length 오류질문 [1] | 촙춉 | 2020.05.26 |
1234 | 메타분석-출판편향 [1] | 날아라 | 2019.08.26 |