library(simFrame)
bar<-rnorm(100,10,1) # data
bar2<-as.data.frame(bar) # 形式记得转化,否则报错
nc<-NAControl(NArate=0.2) #set missing rate
x<-setNA(bar2,nc)
填补categorical数据
library(missMDA) data(vnf) nb <- estim_ncpMCA(vnf,ncp.max=5) ## Time-consuming, nb = 4 res <- MIMCA(vnf, ncp=4,nboot=10)