메뉴 건너뛰기

웹에서 하는 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 396
255 메타분석 결과 관련 질문드립니다. [2] seokzzn 2017.10.18 38
254 메타분석 질문 [2] 남일성 2017.10.06 70
» Cox proportional analysis에서 subgroup분석에 대해 여쭙니다 김규 2017.09.30 54
252 ROC 곡선 문의드립니다. [2] file kimgeunhhyo 2017.09.28 72
251 데이터 전처리 관련 문의드립니다. [7] 화모기 2017.09.28 84
250 Optimal Cut-Off Value에 대해서 두꺼비다 2017.09.25 51
249 연속변수를 사용한 생존 분석 [2] selteen 2017.09.25 56
248 Propensity score matching 에 관하여 질문 드립니다. [1] coma 2017.09.20 88
247 Sigmoid curve (회귀분석) 질문드립니다. [1] file chanheebox 2017.09.19 61
246 Error in f(x, na.rm = TRUE) : 사용되지 않은 인자 (na.rm = TRUE) [2] twogama2 2017.09.18 49
245 교수님, global chi square statistic 관련 문의 입니다. file 16njh 2017.09.17 51
244 venn diagram [1] dave 2017.09.16 40
243 로지스틱 회귀분석에 대해 질문 올립니다. [12] SOL 2017.09.15 142
242 Logistic regression 에 대하여 질문 올립니다. [2] coma 2017.09.14 58
241 밑에 댓글로 한 질문을 새로 드립니다. [2] limjk1114 2017.09.12 63
240 웹에서 cox proportional analysis 하는 방법 [2] limjk1114 2017.09.08 67
239 교수님 오프라인 강의들었던 수강생입니다. [1] 멋진남자 2017.09.05 72
238 Regression analysis와 relative importance에 대하여 문의드립니다. [1] file rokstone 2017.09.05 51
237 Upload파일의 javascript 에러 [1] DrDreamer 2017.08.31 53
236 교차비의 95% 신뢰구간에 관한 질문입니다. [3] 통계초보보보 2017.08.31 70