항상 빠른 답변 감사드립니다.
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 |
---|---|---|---|
1193 | 첫번째 행을 라벨로 사용 버튼 클릭시 에러 문의 [3] | 수담이 | 2023.05.04 |
1192 | 안녕하세요? web-R table 그리고 관련하여 문의드립니다. [4] | mkw7654 | 2023.05.02 |
1191 | 메타분석 U3 그래프 관련 문의드립니다 | 화이팅 | 2023.05.02 |
1190 | 교수님 안녕하세요~~ | 메시12 | 2023.04.29 |
1189 | 연속형 변수인데 범주형으로 인지됩니다. [2] | shararaku | 2023.04.28 |
1188 | 통계마당 [2] | ejo | 2023.04.27 |
1187 | moonBook 오류 관련 재문의드립니다. [1] | holyhyun | 2023.04.24 |
1186 | moonBook 사용시 오류메시지로 문의드립니다. [1] | holyhyun | 2023.04.24 |
1185 | 에러 문의 [2] | 화이팅 | 2023.04.19 |
1184 | Cox regression 에 대한 문의입니다. [2] | saga | 2023.04.18 |
1183 | PSM 문의 [3] | 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] | 소소 | 2023.03.26 |
1176 | 코딩값이 여러개인 변수에 대한 전체 p값을 알 수 있는 방법이 궁금합니다 | 쥬쥬즈주 | 2023.03.22 |
1175 | PSM 매칭 버튼 보이지 않음. | swpapa | 2023.03.19 |
1174 | group-based trajectory modeling | swpapa | 2023.03.19 |