메뉴 건너뛰기

웹에서 하는 R 통계

완전 초보가 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.
>
> pooled
Class: mipo    m = 5 
               estimate        ubar            b           t  dfcom       df          riv       lambda          fmi
ATRH         1.00980260 0.318820965 1.102597e-04 0.318953277 999999 958369.3 0.0004150029 0.0004148307 0.0004169167
log_protein  0.46618838 0.013448747 1.000273e-05 0.013460750 999999 833549.7 0.0008925200 0.0008917241 0.0008941214
BMI1        -0.07057015 0.003690066 5.613955e-06 0.003696803 999999 545840.1 0.0018256438 0.0018223169 0.0018259743
CCV_Hx       0.23557745 0.244975757 2.626380e-05 0.245007274 999999 995749.7 0.0001286518 0.0001286352 0.0001306435

pooled를 확인해 볼시 dfcom이 999999로 나타나서 문제가 있는 것으로 생각했습니다 완전 생초보가 교수님께 고견을 여쭙습니다.... 언제나 감사합니다.


번호 제목 글쓴이 날짜 조회 수
공지 묻고답하기 게시판을 운영합니다 [1] cardiomoon 2015.04.20 592
393 pubmed wordcloud new kaiross 2018.11.18 0
392 상관계수 문의 [1] seyi0918 2018.11.15 6
391 세 군 이상일 때의 KP curve의 해석에 대해 문의드립니다. [1] osis 2018.11.08 18
390 로지스틱 회귀분석 질문드립니다. [2] 윤이아빠 2018.11.06 43
» 교수님 질문드립니다 MICE에 관련한 질문입니다. [2] file vfrider 2018.11.05 11
388 multivariable logistic regression에서 sample size와 multivirate logistic regression에 대해 문의드립니다. [1] namuhanayeyo 2018.11.05 24
387 다중회귀분석 로그변환 후 결과값에 대한 문의 나다니엘 2018.11.05 19
386 logistic regression 으로 구한 predicted prababllity 그래프에서 합치는 방법 문의드립니다. [2] file jh100 2018.11.02 29
385 ROC 커브 관련 질문 [2] PedEndo 2018.11.02 21
384 변수 넣고 빼기 [1] Yoonida 2018.11.01 19
383 교수님..ggBar를 시행하면 계속 에러메시지가 뜹니다..? [4] file 스튜던트 2018.10.30 41
382 ggplot에서 y axis limit정하는 방법 여쭙니다. [2] fascinate23 2018.10.26 27
381 다변량 분석 결과 작성시 기술방법에 대하여 질문 드립니다. [2] drted 2018.10.26 30
380 결측치 메우는 방법 문의 [1] 김태오 2018.10.19 102
379 boxplot 질문드립니다. [2] file jjunglammy 2018.10.17 34
378 출간올유 funnel plot 에 대해서 [2] file 노을 2018.10.15 28
377 교수님! cox 회귀분석 큐빅 스플라인 그래프 작성 질문드립니다! [1] susan 2018.10.15 36
376 메타분석 출간오류에 대하여 문의드립니다. file 노을 2018.10.14 40
375 교수님, 질문 드립니다. kjg 2018.10.06 74
374 로지스틱 회귀분석 김밥친구 2018.10.04 83