I am running several models in R2MLwiN. When I run the models with 1 MCMC chain, I am able to use the resi.store.levs = 2 argument as shown below:
Code: Select all
(mymodel3 <- runMLwiN(math ~ 1 + y + z + (1 | unique_s1) + (1 | student_id), estoptions = list(EstM = 1, mcmcMeth=
list(burnin=500, iterations=5000, thinning=1),
mm = list(list(mmvar = list("s1", "s2", "s3", "s4"),
weights = list("w1", "w2", "w3", "w4")), NA), resi.store = TRUE,
resi.store.levs = 2), data = dataset))
I need to be able to save the level 2 residuals and ideally I would like to run more than one chain. Is this possible? Any help would be greatly appreciated. Thank you!