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