data(wage1, package = "R2MLwiN")
(mymodel <- runMLwiN(logearn ~ 1 + age_40 + numjobs + (1 | company) + (1 | id),
estoptions = list(EstM = 1,
mm = list(list(mmvar = list("company", "company2", "company3", "company4"),
weights = list("weight1", "weight2", "weight3", "weight4")), NA)),
data = wage1))
Run the code above in your browser using DataLab