Search found 1344 matches
- Mon Oct 10, 2022 1:39 pm
- Forum: R2MLwiN user forum
- Topic: bug when setting starting parameters because of R update
- Replies: 3
- Views: 25681
Re: bug when setting starting parameters because of R update
Thank you very much for reporting this problem. It should be fixed in the development version of the package, available from https://github.com/r-forge/r2mlwin/tree/master/R2MLwiN . We are hoping to do a new CRAN release that contains this and some other fixes soon, however we still have to do a few...
- Tue Aug 16, 2022 9:11 am
- Forum: runmlwin user forum
- Topic: Endogeneity issues
- Replies: 3
- Views: 18006
Re: Endogeneity issues
I asked George about this and his reply was as follows: I guess you could apply ideas from fixed versus random effects to panel data which are motivated by endogeneity concerns. So you could enter region as fixed-effects dummy variables instead of as random effects. The effects of the firm covariate...
- Tue Jul 19, 2022 12:31 pm
- Forum: runmlwin user forum
- Topic: svyset compatibility with runmlwin
- Replies: 1
- Views: 8067
Re: svyset compatibility with runmlwin
I asked George about this and he suggested that the weighting options in MLwiN might do what you want. Unfortunately we don't provide much documentation for this, but you can see how to specify them by looking at the options in the runmlwin help file, i.e.
Code: Select all
help runmlwin
- Mon Jul 04, 2022 2:36 pm
- Forum: MLwiN user forum
- Topic: 2 responses using MLwiN3.01
- Replies: 1
- Views: 7645
Re: 2 responses using MLwiN3.01
The reason for this is that your have chosen Binomial responses. This response type has a fixed level-1 variance (level-2 in this case as everything is shifted up a level for multivariate responses), which therefore cannot be estimated. If you had fitted a univariate Binomial response you would find...
- Fri Jun 17, 2022 1:51 pm
- Forum: MLwiN user forum
- Topic: Calculating level-three residuals in multilevel model
- Replies: 2
- Views: 7959
Re: Calculating level-three residuals in multilevel model
If you are using MLwiN then the residuals returned from the Model>Residuals menu are already shrunken (see chapter 3 of the MLwiN User's Guide ). If you need an unshrunken version then the REFLATE command (see the MLwiN on-line help) is provided to calculate this, and more details of the formula use...
- Tue Jun 14, 2022 11:05 am
- Forum: R2MLwiN user forum
- Topic: informative priorParam in random slope
- Replies: 2
- Views: 9623
Re: informative priorParam in random slope
As you state the 'prior variance matrix is not positive definite' message usually occurs when the starting values for the covariance matrix, used as the prior for MCMC, is invalid. Normally this matrix is obtained from the IGLS estimates, however you can specify it manually via the startval option. ...
- Tue Apr 12, 2022 1:14 pm
- Forum: runmlwin user forum
- Topic: Level 2 unit problem
- Replies: 1
- Views: 7635
Re: Level 2 unit problem
Although you have checked the sorting, this does look like the sort of behaviour you would get if the data was not sorted by the identifiers defining the unit hierarchy. Can you confirm that the data is sorted by all of these ('n', 'wave_num' and 'wave')? You could also try recoding 'n' to go from 1...
- Thu Mar 10, 2022 9:04 pm
- Forum: R2MLwiN user forum
- Topic: how to extract credible intervals
- Replies: 1
- Views: 8371
Re: how to extract credible intervals
For MCMC models the parameter chains are stored in a chains slot of the returned object, in your case this would be modelstim@chains . You can use the functions provided in the coda package to manipulate these, for example: library(coda) summary(modelstim@chains) will list summary statistics for eac...
- Fri Feb 25, 2022 1:41 pm
- Forum: runmlwin user forum
- Topic: melogit verse runmlwin
- Replies: 4
- Views: 8196
Re: melogit verse runmlwin
I asked George if he could think of an explanation for this and he suggested running the model with PQL2 instead of MQL1, as MQL1 variances are biased downwards.
- Thu Feb 24, 2022 1:10 pm
- Forum: runmlwin user forum
- Topic: melogit verse runmlwin
- Replies: 4
- Views: 8196
Re: melogit verse runmlwin
You should be able to get this to work if you remove the predictors from the level-1 specification, i.e.: generate id = _n runmlwin depression cons age10, level2(mcsid: cons) level1(id:) discrete(distribution(binomial) link(logit) denominator(cons)) nopause You can find examples of this and other mo...