www.cmm.bristol.ac.uk/forum

I have a variable with 5 levels. Level 5 is now automatically taken as the reference category (that is: I get intercept_1 thru intercept_4), but I want level 1 to be the reference category (i.e., have intercept_2 thru intercept_5 estimated). What code do I have to add to change that? If this informa...

That makes sense, thank you.

I am attempting to use ordered logistic regression using R2MLwiN. As a first step, I am looking at the R2MLwiN script for the MCMC manual example. In chapter 13.3, pupil is set as the level 1 identifier. In the R2MLwiN script, pupil is not included: ##IGLS (mymodel <- runMLwiN(logit(a_point, cons, 6...

Hi Bill,

If you could take a look that would be great. Just PB me with your email address.
I will try to fiddle with the parameter expansion a bit as well.
I have to note that I already use orth=1, as this seemed to help in some other models.

That does help, thank you.

Just to be clear, the nesting of classrooms within schools, and nesting of ID in both classT1 and classT2, will be automatically recognized?

For me the only solution was running the model via R2mlwin, because there I could tweak the starting values.

I have a dataset with the following levels: School (99 schools) Class * (see below) ID (4249 children) Time (3 timepoints, T1, T2, T3) For Class, there is cross-nesting going on: At T1, children were in their classT1 (N=250), At T2 and T3, children were in classT2 (N=232). To be precise, we measured...

As a first step in my analyses, I am running relatively simple logistic regression models using MCMC estimation. My simplest model has one outcome, N=4249 at level 1, N=232 at level 2, N=99 at level 3. Even for this model, Nhat keeps increasing and convergence does not seem to be able to be reached ...

Thanks! Unfortunately, as I am not used to programming with R code, relatively easy is still not feasible for me (I can adopt a manual to my own data, but that's about it...).
Could you give an example of how to extract, for instance, Nhat and quartiles for a specific parameter are derived?

Is it possible to save the output from sixway function (accuracy statistics and summary diagnostics) without the figures in an output? I have to run about a million analyses (slight exaggeration) and I would just like 1 output with all of the numbers. I already save the trajectories and sixway outpu...