메뉴 건너뛰기

웹에서 하는 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 431
272 web r을 이용한 box plot 그리기 중 질문드립니다. [4] file 민세권세민 2017.12.09 27
271 Trend test에 대한 질문 하가지만 드립니다. [5] file namuhanayeyo 2017.11.30 22
270 한가지 문의드리고자 합니다. [2] namuhanayeyo 2017.11.30 8
269 비비교 연구에 대한 두 메타분석을 서로 비교할 수 있는지 궁금합니다. [3] SOL 2017.11.24 22
268 표만들기 결과 값이 이해가 안되어 조언 부탁드리고자 올립니다. [1] file 멋진홍이 2017.11.23 21
267 Table에서 출력되는 수치의 정확도 관련 [9] o121 2017.11.21 58
266 mytable() 함수에 대한 질문입니다. [3] 브리디온 2017.11.19 37
265 교수님, 다중공선성 질문드립니다. [2] selteen 2017.11.16 45
264 Data favorite preprocessing 문의드립니다. [2] 멋진홍이 2017.11.09 25
263 안녕하세요 교수님. 여쭙고 싶은 것이 있어 글 올립니다. [2] yahello 2017.11.08 44
262 ZCOR 값이 안나옵니다! [2] seokzzn 2017.11.02 35
261 R에서 변수변환을 통한 파생변수 생성 질문드립니다.. [1] file JJSWHY 2017.11.01 37
260 mytable 사용 시 결측치의 표시 [1] twogama2 2017.10.30 29
259 data오류에 대해 문의드립니다. [1] file FAdavid 2017.10.30 23
258 교수님, 메타분석에 대한 질문 몇가지 올립니다. SOL 2017.10.26 27
257 cut off value 구하기 KCRS_LeeJM 2017.10.26 38
256 ggplot 이용시 자꾸 오류가 납니다~ [1] 돌님 2017.10.23 49
255 메타분석 결과 관련 질문드립니다. [2] seokzzn 2017.10.18 72
254 메타분석 질문 [2] 남일성 2017.10.06 112
» Cox proportional analysis에서 subgroup분석에 대해 여쭙니다 김규 2017.09.30 134