완전 초보가 web-R로 통계 입문해서 많은 도움을 받고 있습니다. 감사합니다 교수님.
"의학논문작성을 위한 R 통계와 그래프"책을 보고 결측치 처리를 위해 MICE를 사용하여 처리를 연습하고 있습니다만
MICE함수에 cox regression 함수를 집어넣어서 처리하는 것이 가능한지요? 경고메시지가 뜨면서 처리값중 df가 제대로 처리되지 않는 것처럼 보였습니다.
첨부된 파일은 전체파일로 돌리면 너무 느려서 연습용으로 아주 일부만 발췌한 파일입니다만,
첨부파일중 SUBJECT_NO를 제외하고 나머지 변수를 투입해서 RE_EVENT를 사건, RE_YEAR를 생존기간으로 해서 cox regression을 할 계획이었습니다.
책을 참고하여 넣은 명령문은 다음과 같았습니다. 여기서는 들어가는 변수를 4개만 넣어서 연습해 보았었습니다.
data=read.csv("CKD_COX100.csv")
data
require(mice)
require(survival)
imp=mice(data,m=5,seed=1234)
fit1=with(imp,coxph(Surv(RE_YEAR,RE_EVENT==1)~ATRH+log_protein+BMI1+CCV_Hx))
pooled=pool(fit1)
summary(pooled)
pooled
경고메시지는 pool함수 사용시 발생했고 다음과 같았습니다.
> pooled=pool(fit1)
Warning messages:
1: Unknown or uninitialised column: 'df.residual'.
2: In pool.fitlist(getfit(object), dfcom = dfcom) : Large sample assumed.
| |
|
Comment 2
-
cardiomoon
2018.11.07 11:50
-
vfrider
2018.11.10 12:26
친절한 답변 감사드립니다 교수님
말씀해주신 레퍼런스 공부해보고 열심히 답을 찾아보겠습니다^^
이 부분은 제가 공부가 부족해서 지금 답변드리기는 어렵네요. 하지만 mice패키지 홈페이지에 가보니 다음 책이 reference로 달려있고 원본을 모두 읽을 수 있습니다. 책 이름을 클릭하셔요.
Van Buuren, S. (2018). Flexible Imputation of Missing Data. Second Edition.. Chapman & Hall/CRC. Boca Raton, FL.