Skip to menu

생존분석에서 missing value가 있는 경우 coxph 모델링은 어떻게 하는지요..

PipeDragoN 2017.01.18 15:35 Views : 412

* 교수님, 아래와 같이 model을 만들어 보았을때, 207개의 AchRAbtiter의 자료 중에 7개의 자료가 missing value가 있어서 variable이 0과1로 두개로 나누어지는 거 같습니다. 그래서 web-R의 결과와 똑같이 R console에서 구현이 안되는데요...

coxph는 na.rm=T와 같은 추가할 수 있는 항목이 없는지요?? Web-R의 결과를 R에서 구현해보고자 합니다. 

* web-R에서는 다름과 같이 표현되었습니다.

For stepwise backward elimination, data containing NA value were removed.

Total nrow: 207-> After na.omit : nrow: 200 

Total  7  row(s) were removed.


> cox2=coxph(Surv(MMfreeSurvival,MM==1)~AchRAbTiter_6.7_or_more,method="breslow")

> summary(cox2)

Call: coxph(formula = Surv(MMfreeSurvival, MM == 1) ~ AchRAbTiter_6.7_or_more, method = "breslow")

  n= 207, number of events= 134 

                            coef exp(coef) se(coef)     z Pr(>|z|)

AchRAbTiter_6.7_or_more0 0.66269   1.94001  0.52583 1.260    0.208

AchRAbTiter_6.7_or_more1 0.09358   1.09810  0.51150 0.183    0.855


                         exp(coef) exp(-coef) lower .95 upper .95

AchRAbTiter_6.7_or_more0     1.940     0.5155    0.6922     5.437

AchRAbTiter_6.7_or_more1     1.098     0.9107    0.4029     2.993


Concordance= 0.555  (se = 0.02 )

Rsquare= 0.039   (max possible= 0.998 )

Likelihood ratio test= 8.28  on 2 df,   p=0.01596

Wald test            = 9.01  on 2 df,   p=0.01105

Score (logrank) test = 9.26  on 2 df,   p=0.009772


I want to
List
No. Subject Author Date
1014 생존분석에서 multivariate analysis에 대해 여쭙니다. [3] joon 2016.12.23
1013 R통계와 그래프 책을 보며 따라하는데... [2] 쿠싸길스패츌러 2016.12.23
1012 cox ph 모형 자동화 [1] hawk1227 2017.01.05
1011 교수님 문의드립니다~ [3] neurokim 2017.01.06
1010 matchit function 의 caliper 관련 질문입니다. 자라선생 2017.01.11
1009 생존나무분석에서 예측과 실제의 비교와 관하여 질문드립니다. (죄송하게도 또 글을 올립니다.) 자라선생 2017.01.11
1008 기본적인 data 파일 전환 방법 문의 [1] sunny112 2017.01.16
1007 구조방정식 강의 [1] file cardiomoon 2017.01.17
1006 1-생존함수 plot 문의드립니다. [3] file PipeDragoN 2017.01.18
» 생존분석에서 missing value가 있는 경우 coxph 모델링은 어떻게 하는지요.. [2] PipeDragoN 2017.01.18
1004 1-생존함수에 관하여.. [1] PipeDragoN 2017.01.21
1003 교수님 문의드립니다~ [5] chamchi 2017.01.22
1002 생존분석에서의 interaction term 관련 여쭙습니다. Ridenc 2017.01.29
1001 로지스틱 회귀분석 질문 있습니다. [1] chamchi 2017.02.01
1000 Trend test [1] namuhanayeyo 2017.02.04
999 ANCOCA 질문드립니다. [1] BlcAng 2017.02.04
998 R notebook 에서 mytable 결과가 이상하게 나옵니다. [4] file Dilemma 2017.02.05
997 질문드립니다. [1] Dr.Han 2017.02.05
996 word cloud : pubmed error에 관하여 문의 드립니다. [3] file rhrroddl1 2017.02.06
995 web-r의 ggplot2를 사용해 보고 [2] 장미공원 2017.02.07
Write Tag