항상 빠른 답변 감사드립니다.
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 |
---|---|---|---|
1053 | 통계와 관련한 기초적인 질문 드립니다.. [2] | 감사합니다 | 2021.12.10 |
1052 | 메타분석 해석 문의 [2] | 루아로 | 2021.12.07 |
1051 | 초보라 질문드립니다 ㅜ [2] | 쵸쿄쵸쿄 | 2021.12.06 |
1050 | ㅠㅠㅠ 초보 질문입니다 [1] | 만고 | 2021.11.26 |
1049 | 라벨붙이기가 안됩니다 [1] | 찬미맘 | 2021.11.22 |
1048 | moonBook ANOVA 질문드립니다. [1] | 타쿠 | 2021.11.20 |
1047 | 3원 분할표? 에서 모델 만드는법 | 플랫슬립 | 2021.11.20 |
1046 | [회귀분석 초급 질문] [1] | 서빱스 | 2021.11.14 |
1045 | [ANOVA 분석 질문] [1] | 서빱스 | 2021.11.14 |
1044 | 대조군1과 치료군2,3,4 의 PSM 에 관해서 질문있습니다. [4] | yessul | 2021.11.11 |
1043 | 서버접속이 안됩니다. Reload [2] | silkywater | 2021.11.11 |
1042 | 메타분석 실행이 안됩니다. [1] | 애플K | 2021.11.08 |
1041 | pair plot 관련하여 질문드립니다. | ulysses11 | 2021.10.29 |
1040 | 메타분석 결과파일 다운로드 하려면 [1] | mono | 2021.10.23 |
1039 | 기본 baseline character table문의드립니다. [1] | ulysses11 | 2021.10.22 |
1038 | 이 error는 뭐가 문제인 걸까요 [1] | CSJU | 2021.10.19 |
1037 | logistic regression에서 p>0.05인 변수의 OR, p 값... | 브이콘65 | 2021.10.14 |
1036 | R 문제가 도저히 안풀려요.. [1] | 자하연미필 | 2021.10.09 |
1035 | 구조방정식에서 종속변수가 0,1로 코딩된 데이터 사용이 가능한가요? [1] | 순수엄마 | 2021.10.08 |
1034 | roc.test를 통해 AUC 비교 [2] | 나현둥 | 2021.10.07 |