> library(mice)
> library(lavaan)
> library(semTools)
>
> bemmiss <- read.csv("~/desktop/bem.csv",header=T)
>
>
>
>
> bemmiss$assert[bemmiss$assert == 1] <- NA
> bemmiss$softspok[bemmiss$softspok == 7] <- NA
> bemmiss$stand[bemmiss$stand >=6] <- NA
>
> bemmiss.model <- ' masculine =~ assert+forceful+dominant+stand
+ feminine =~ softspok+warm+truthful+tender+gentle '
>
> bemmiss_imp <- mice(bemmiss, m = 20, method="pmm")
iter imp variable
1 1 assert stand softspok
1 2 assert stand softspok
1 3 assert stand softspok
1 4 assert stand softspok
1 5 assert stand softspok
1 6 assert stand softspok
1 7 assert stand softspok
1 8 assert stand softspok
1 9 assert stand softspok
1 10 assert stand softspok
1 11 assert stand softspok
1 12 assert stand softspok
1 13 assert stand softspok
1 14 assert stand softspok
1 15 assert stand softspok
1 16 assert stand softspok
1 17 assert stand softspok
1 18 assert stand softspok
1 19 assert stand softspok
1 20 assert stand softspok
2 1 assert stand softspok
2 2 assert stand softspok
2 3 assert stand softspok
2 4 assert stand softspok
2 5 assert stand softspok
2 6 assert stand softspok
2 7 assert stand softspok
2 8 assert stand softspok
2 9 assert stand softspok
2 10 assert stand softspok
2 11 assert stand softspok
2 12 assert stand softspok
2 13 assert stand softspok
2 14 assert stand softspok
2 15 assert stand softspok
2 16 assert stand softspok
2 17 assert stand softspok
2 18 assert stand softspok
2 19 assert stand softspok
2 20 assert stand softspok
3 1 assert stand softspok
3 2 assert stand softspok
3 3 assert stand softspok
3 4 assert stand softspok
3 5 assert stand softspok
3 6 assert stand softspok
3 7 assert stand softspok
3 8 assert stand softspok
3 9 assert stand softspok
3 10 assert stand softspok
3 11 assert stand softspok
3 12 assert stand softspok
3 13 assert stand softspok
3 14 assert stand softspok
3 15 assert stand softspok
3 16 assert stand softspok
3 17 assert stand softspok
3 18 assert stand softspok
3 19 assert stand softspok
3 20 assert stand softspok
4 1 assert stand softspok
4 2 assert stand softspok
4 3 assert stand softspok
4 4 assert stand softspok
4 5 assert stand softspok
4 6 assert stand softspok
4 7 assert stand softspok
4 8 assert stand softspok
4 9 assert stand softspok
4 10 assert stand softspok
4 11 assert stand softspok
4 12 assert stand softspok
4 13 assert stand softspok
4 14 assert stand softspok
4 15 assert stand softspok
4 16 assert stand softspok
4 17 assert stand softspok
4 18 assert stand softspok
4 19 assert stand softspok
4 20 assert stand softspok
5 1 assert stand softspok
5 2 assert stand softspok
5 3 assert stand softspok
5 4 assert stand softspok
5 5 assert stand softspok
5 6 assert stand softspok
5 7 assert stand softspok
5 8 assert stand softspok
5 9 assert stand softspok
5 10 assert stand softspok
5 11 assert stand softspok
5 12 assert stand softspok
5 13 assert stand softspok
5 14 assert stand softspok
5 15 assert stand softspok
5 16 assert stand softspok
5 17 assert stand softspok
5 18 assert stand softspok
5 19 assert stand softspok
5 20 assert stand softspok
> mice.imp <- NULL
> for(i in 1:20) { mice.imp[[i]] <- complete(bemmiss_imp, action=i, inc=FALSE)
+ }
>
> # Run lavaan with previously imputed data using runMI