## Program to run Mplus for the simulation ## function to calculate the BC for mplus bc.ci<-function(res, est){ iB<-length(res) ## BC interval alpha<-c(.025, .975) z0<-qnorm(sum(res