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
1094 안녕하세요 r 코드 중에 도저히 해결이 안나는 문제가 있어서 여기에 문의 드립니다 [1] file 슈퍼맨 2022.06.24
1093 Box plot 그리프 설정에 관해 문의합니다. papahong 2022.06.20
1092 cluster data 의 risk difference 에 대해 여쭙고자 합니다. lagrima 2022.06.19
1091 Pair plot을 할 경우 [1] PipeDragoN 2022.06.13
1090 PSM에서 matching 버튼 위치 문의 통계어려워 2022.06.13
1089 PSM (4그룹) 관련 문의드립니다. emksh77 2022.06.08
1088 정회원후원하였습니다. [1] secret rosemary 2022.06.08
1087 PSM 문의 [2] secret sfumatoi 2022.06.07
1086 multipleROC 패키지에서 multivariable의 cutoff 구하기 file 찬쓰 2022.06.07
1085 survival analysis관련 문의 드립니다 [1] cebukim 2022.06.04
1084 autoReg 설치 관련 오류 [5] 닥터강 2022.05.24
1083 survival tree (RPA) 분석 관련 문의 드립니다. [2] 곰곰11 2022.05.23
1082 통계 질문입니다. [1] file 봄봄봄 2022.05.23
1081 데이터 불러오기 문제 [3] strat 2022.04.25
1080 r프로그래밍 문자는 factor로 읽을 수 없나요? [1] file 채채채추ㅐㅊ채야채 2022.04.21
1079 R 코드와 관련하여 문의 드려도 될까요..? [5] file 꾸이뭉 2022.04.19
1078 R 업데이트 후 mytable 함수 에러가 있습니다. [1] 카이토르 2022.04.18
1077 Survival curve 에서 cumulative incidence % 로 표현 가능할까요 [1] nyknicks 2022.04.16
1076 multiple logistic regression test에서 설명변수 선택에 대해 여쭙습니다. [2] file kms9057 2022.04.16
1075 PSM 문의 드립니다 [1] caesar 2022.04.07
Write Tag