# NOT RUN {
data(gcsemv1, package = "R2MLwiN")
(mymodel <- runMLwiN(c(written, csework) ~ 1 + female + (1 | school) + (1 | student),
D = "Multivariate Normal", estoptions = list(EstM = 1), data = gcsemv1))
# }
Run the code above in your browser using DataLab