www.cmm.bristol.ac.uk/forum

The first part is because I'm having to remote desktop into a computer that is not the computer I am posting from to use MLwiN. I did encounter the issue of the model typo and fixed it in the code on the remote desktop (but not the copy on this computer).

I think I probably breezed through the text ...

It's all the same length and still sequential.

In fact, if I run through the process manually i.e. estimate an IGLS model for starting values then switch to MCMC and set the classification to multiple membership the model estimates just fine.

Something has occurred to me in the interim - some of ...

And just in case it comes up, my membership/weight columns are sequential so that should not be an issue.

> names(mlwin_dat)
[1] "ID" "MM_s1" "MM_s2" "MM_s3" "S6Read"
[6] "S3Read" "wFL" "S3Read_gm" "S3Read_clus" "wcMean"
[11] "w1" "w2" "w3" "intercept"

Following the format of the vignette for multiple membership data, I ran this code:


mm <- list(list(
mmvar = list("MM_s1", "MM_s2", "MM_s3"),
weights = list("w1", "w2", "w3")),NA)

Funcet <- S6Read ~ 1 + S3Read + wFL + (1|MM_s1)

uncentlm <- runMLwiN(Formula = Funcent, data = mlwin_dat ...