Skip to menu

* 교수님, 아래와 같이 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