항상 빠른 답변 감사드립니다.
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 |
---|---|---|---|
1234 | 안녕하세요 에러 관련 문의드립니다 [1] | 화이팅 | 2023.09.19 |
1233 | 선형회귀분석 시 보정변수 사용 관련 질문드립니다. [2] | 전한울 | 2023.09.18 |
1232 | 그래프 관련 문의드립니다. | nanak | 2023.09.14 |
1231 | 영수증 메일로 받을수 있을까요? [1] | 딸기공주 | 2023.09.13 |
1230 | mycsv 관련 문의드립니다. [1] | 도삼 | 2023.09.12 |
1229 | matched group 을 대상으로 stratified cox | saga | 2023.09.06 |
1228 | IPTW 시 SMD 계산 관련 여쭙습니다. | 장한나 | 2023.09.04 |
1227 | 회귀분석 관련 문의드립니다. [2] | 전한울 | 2023.09.04 |
1226 | PSM matching 실행시 disconnected 됩니다... [1] | 맹구모션 | 2023.08.25 |
1225 | 메타분석_출간오류분석 관련 질문드립니다. [2] | 블루레몬 | 2023.08.21 |
1224 | plot 문의 드립니다. | 1122 | 2023.08.20 |
1223 | semMediation 설치문의 [2] | 마술샘 | 2023.08.11 |
1222 | 로지스틱 회귀분석 문의드립니다. [2] | JKJ | 2023.08.09 |
1221 | 지출 문의 드립니다. [1] | 블루레몬 | 2023.08.08 |
1220 | autoReg 실행 중 에러 [1] | pkjsdf123 | 2023.08.08 |
1219 | 영수증 처리 [1] | britnepak | 2023.07.28 |
1218 | multipleROC 커브 글을 따라가다 궁금한 것이 생겨서 질문드립니다. [2] | rgppwls | 2023.07.20 |
1217 | disconnected from server 오류 [2] | JKJ | 2023.07.13 |
1216 | mytable에서 변수 순서와 ztable적용시 변수가 숫자로 변형 [1] | 둥니맨 | 2023.07.12 |
1215 | box plot 시 X 축 label이 다 나오지 않아 질문드립니다. [2] | 난둥 | 2023.07.03 |