amices / mice

Multivariate Imputation by Chained Equations
https://amices.org/mice/
GNU General Public License v2.0
447 stars 108 forks source link

Pool for cox hazard model #29

Closed gajyakichi closed 7 years ago

gajyakichi commented 7 years ago

Thank you for great package, mice.

In cox proportional hazard model pooling, the result showed est was less than zero.

It might be confused with linear regression model.

How is the best method to pool with cox proportional hazard model?

library(mice)
library(survival)

#Making database
Tempmice<-mice(imputedmodel,
m=10,maxit=50,meth='pmm', seed=500)

#Cox proportional hazard model
CoxModel.1<- coxph(Tempmice,Surv(FU, Censor==1)~ Factor1+Factor2, method="breslow")

#Pool
fit<-with(data=Tempmice,CoxModel.1)
coxphimputed<-summary(pool(fit))
coxphimputed
stefvanbuuren commented 7 years ago

The sequence of statements seems OK to me. The Cox model allows for negative coefficients.