www.cmm.bristol.ac.uk/forum

Hi, I have a similar problem before. However, the problem was solved after I upgraded the MLwiN to version 3.02 and installed R2MLwiN using http://www.bristol.ac.uk/cmm/software/r2mlwin/


Vivian

I get what you mean. Thanks!!

Do you mean in this case I need to manually add 4 columns of Intercept into the data frame?

Thanks,
Vivian

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(gcseav ...

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 random slope model ...

Thanks a lot.

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 ...

Thanks Billb, your answers clarify my questions!

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 you want ...

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 ...