메뉴 건너뛰기

웹에서 하는 R 통계

 항상 빠른 답변 감사드립니다.


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을 구하는 데에서 막여 문의 드립니다.


감사합니다.



번호 제목 글쓴이 날짜 조회 수
공지 묻고답하기 게시판을 운영합니다 [1] cardiomoon 2015.04.20 574
260 mytable 사용 시 결측치의 표시 [1] twogama2 2017.10.30 91
259 data오류에 대해 문의드립니다. [1] file FAdavid 2017.10.30 79
258 교수님, 메타분석에 대한 질문 몇가지 올립니다. SOL 2017.10.26 121
257 cut off value 구하기 KCRS_LeeJM 2017.10.26 334
256 ggplot 이용시 자꾸 오류가 납니다~ [1] 돌님 2017.10.23 221
255 메타분석 결과 관련 질문드립니다. [2] seokzzn 2017.10.18 108
254 메타분석 질문 [2] 남일성 2017.10.06 281
» Cox proportional analysis에서 subgroup분석에 대해 여쭙니다 김규 2017.09.30 290
252 ROC 곡선 문의드립니다. [2] file kimgeunhhyo 2017.09.28 324
251 데이터 전처리 관련 문의드립니다. [7] 화모기 2017.09.28 194
250 Optimal Cut-Off Value에 대해서 두꺼비다 2017.09.25 304
249 연속변수를 사용한 생존 분석 [2] selteen 2017.09.25 131
248 Propensity score matching 에 관하여 질문 드립니다. [1] coma 2017.09.20 286
247 Sigmoid curve (회귀분석) 질문드립니다. [1] file chanheebox 2017.09.19 128
246 Error in f(x, na.rm = TRUE) : 사용되지 않은 인자 (na.rm = TRUE) [2] twogama2 2017.09.18 189
245 교수님, global chi square statistic 관련 문의 입니다. file 16njh 2017.09.17 79
244 venn diagram [1] dave 2017.09.16 83
243 로지스틱 회귀분석에 대해 질문 올립니다. [12] SOL 2017.09.15 1653
242 Logistic regression 에 대하여 질문 올립니다. [2] coma 2017.09.14 301
241 밑에 댓글로 한 질문을 새로 드립니다. [2] limjk1114 2017.09.12 104