Hi Chris, Thanks for your quick reply. I replicate the example from http://www.bristol.ac.uk/cmm/media/r2mlwin/MCMCGuide13.R . data(alevchem, package = "R2MLwiN") alevchem$gcseav <- double2singlePrecision(alevchem$gcse_tot/alevchem$gcse_no - 6) model1 <- runMLwiN(a_point ~ 1 + gcseav + I(g...
Hi, I'm following the suggestion from this post https://www.cmm.bristol.ac.uk/forum/viewtopic.php?f=7&t=2010&p=4572&hilit=predict&sid=a398d4a33ea73179434b8589781fa976#p4572 to replicate the standardised residual vs predicted value within R. I am two models. Model 1 is a 2-level rando...
Hi, I am running a cross-classified cumulative logit model for my ordinal outcome. I used the following codes to store the residual: estoptions = list(resi.store = TRUE, resi.store.levs = c(2, 3, 4), xc = TRUE, EstM = 1 However, I have this error: Error: memory exhausted (limit reached?) In this cas...
Hi, I would like to calculate the population-averaged probability for a cross-classified ML model. I have read this post: https://www.cmm.bristol.ac.uk/forum/viewtopic.php?f=3&t=1558&p=3867&hilit=predicted+probability&sid=16e1ae467bd040ff6cbae53ea3c523fb#p3867 , 2. Depends on whether...
Hi, I'm using data (9.1wsz) from Module 9. The data is converted to dta and then imported to R. The categorical variable (occtype) is converted into dummy variables. Here is the syntax for replicating the random slope model: data9.1$occtype0 <- dummy.code(data9.1$occtype) RSmodel <- runMLwiN(logit(e...