www.cmm.bristol.ac.uk/forum

It looks as if the estimates replay command does not carry across the display options set when running the command, however you can re-specify them when displaying the results, for example: . use http://www.bristol.ac.uk/cmm/media/runmlwin/bang, clear . logit use age, or Iteration 0: log likelihood ...

The files that you provided ran successfully in both MLwiN versions 3.04 and 3.05 for me. You do not have to open a dataset prior to running the macro as all the values used are simulated. I set the current directory to the location of the macro files via Options>Directories, and opened the simu.txt...

I am not sure where the dashes in that syntax are coming from. Would it be possible to post the generated script, and if possible the inputs you gave to MLPowSim to create it?

If you are running a binary logit model then the variance is fixed at (π^2)/3 (≈3.29) rather than being an estimated parameter, and this is why you cannot add this to your model.

See the following FAQ for more information about variances: Partitioning variation across levels

I showed this to George, but he wasn't able to say for sure. His advice was that unless the models ran too slowly or couldn't be fitted with mixed you should use that. If you can't do that then he suggests initially running a subsample or simplified version of your model using both mixed and your ma...

This forum is primarily for questions related to the "R2MLwiN" R package. You might get an answer if you ask about related packages, but there are often more appropriate places for these questions.

If you are running the model with MCMC ( EstM=1 ) and do not specify starting values then it runs IGLS first and uses the result from this as starting values. Otherwise it uses the values that you provide. You can see this in lines 1496-1520 of https://github.com/r-forge/r2mlwin/blob/master/R2MLwiN/...

Bill gives some reasons why the variance can increase in his first comment in the following thread. It might be worth seeing whether that helps to explain what you are seeing too:

https://www.cmm.bristol.ac.uk/forum/vie ... f=1&t=3624

You can find an example of setting starting values for parameters in section 5.5 of the MCMC manual replication materials. If you do not specify these then the model will be run with (R)IGLS (up to the specified maximum number of iterations) to generate starting values before the MCMC run begins. Fo...

R2MLwiN does not currently provide an implementation of the VarCorr() function, however it does return these values within the output of coef() . The following shows the same model being fitted with glmer , brms , MLwiN(IGLS) and MLwiN(MCMC) : data(mmmec, package = "R2MLwiN") library(lme4...