if (FALSE) {
data(TitanicNA)
## First the number of components has to be chosen
## (for the reconstruction step)
## nb <- estim_ncpMCA(TitanicNA) ## Time-consuming, nb = 5
## Multiple Imputation
res.mi <- MIMCA(TitanicNA, ncp=5, verbose=TRUE)
## First completed data matrix
res.mi$res.MI[[1]]
## Analysis and pooling with mice
require(mice)
imp<-prelim(res.mi,TitanicNA)
fit <- with(data=imp,exp=glm(SURV~CLASS+AGE+SEX,family = "binomial"))
res.pool<-pool(fit)
summary(res.pool)
}
Run the code above in your browser using DataLab