Search found 1335 matches

by ChrisCharlton
Thu Feb 24, 2022 1:10 pm
Forum: runmlwin user forum
Topic: melogit verse runmlwin
Replies: 4
Views: 5254

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...
by ChrisCharlton
Mon Feb 21, 2022 4:34 pm
Forum: R2MLwiN user forum
Topic: Saving Level 2 Residuals When Running Multiple Chains
Replies: 2
Views: 3521

Re: Saving Level 2 Residuals When Running Multiple Chains

Thanks for letting us know about this, the error appears to be due to a bug. I have made a fix in the development version of the package which once the r-forge mirror ( https://github.com/r-forge/r2mlwin/tree/master/R2MLwiN ) has synchronised you should be able to install with the commands: library(...
by ChrisCharlton
Mon Feb 14, 2022 2:21 pm
Forum: runmlwin user forum
Topic: MCMC: calculating estimates, CrIs for combinations of terms/interactions
Replies: 7
Views: 5926

Re: MCMC: calculating estimates, CrIs for combinations of terms/interactions

The p-values reported here are one-tailed bayesian p-values, and essentially give the proportion of the chain that crosses zero (or one if the variable is exponentiated and it is aware of this). You can find the code for this in runmlwin_mcmcdiag.ado , and is as follows: local compval 0 if "`ef...
by ChrisCharlton
Fri Feb 11, 2022 1:55 pm
Forum: runmlwin user forum
Topic: MCMC: calculating estimates, CrIs for combinations of terms/interactions
Replies: 7
Views: 5926

Re: MCMC: calculating estimates, CrIs for combinations of terms/interactions

It should be fairly easy to modify mcmcsum to do this. You can see where the r-class is filled in on lines 535-552 or 410-427 of mcmcsum . To add Bayesian p-values you just need to decide whether you want the mean, median or mode version and add an extra line to return this, i.e. return scalar pvalu...
by ChrisCharlton
Wed Feb 09, 2022 3:48 pm
Forum: MLwiN user forum
Topic: mlnscript for Rocky Linux/OpenBLAS
Replies: 5
Views: 4444

Re: mlnscript for Rocky Linux/OpenBLAS

Thank you for the update, it's useful to know that this method works in this situation.
by ChrisCharlton
Mon Feb 07, 2022 10:43 am
Forum: runmlwin user forum
Topic: MCMC: calculating estimates, CrIs for combinations of terms/interactions
Replies: 7
Views: 5926

Re: MCMC: calculating estimates, CrIs for combinations of terms/interactions

I think that you should be able to do this with the normal variable manipulation commands such as generate and egen . Once you have the new variable representing the result of the calculation you should be able to get the credible intervals using the summarize command, or the mcmcsum provided with r...
by ChrisCharlton
Tue Feb 01, 2022 4:32 pm
Forum: MLwiN user forum
Topic: mlnscript for Rocky Linux/OpenBLAS
Replies: 5
Views: 4444

Re: mlnscript for Rocky Linux/OpenBLAS

You could try using LD_PRELOAD to switch to OpenBLAS as described in Is it possible to switch between BLAS libraries without recompiling program?, i.e.

Code: Select all

LD_PRELOAD=/usr/lib64/libopenblas.so.0 mlnscript --run script.mac
I have not tested how well this works however.
by ChrisCharlton
Tue Feb 01, 2022 1:06 pm
Forum: MLwiN user forum
Topic: mlnscript for Rocky Linux/OpenBLAS
Replies: 5
Views: 4444

Re: mlnscript for Rocky Linux/OpenBLAS

I have added a package version for Rocky Linux 8, which can be requested from https://www.cmm.bristol.ac.uk/clients/softwaredownload/ . Note however that this will still need the libblas.so.3 and liblapack.so.3 files to be present in order to run: ldd /usr/bin/mlnscript linux-vdso.so.1 (0x00007ffdc9...
by ChrisCharlton
Mon Jan 17, 2022 12:17 pm
Forum: runmlwin user forum
Topic: Calculating statistical significance of correlations in multivariate models
Replies: 2
Views: 4996

Re: Calculating statistical significance of correlations in multivariate models

I asked George about this and his response was that you can use the nlcom command for this, i.e.:

Code: Select all

nlcom [RP2]cov(cons_1\cons_2)/sqrt([RP2]var(cons_1)*[RP2]var(cons_2)
by ChrisCharlton
Fri Jan 07, 2022 12:15 pm
Forum: Realcom user forum
Topic: Start Analysis error in MLwiN after retrieving REALCOM impute
Replies: 4
Views: 11798

Re: Start Analysis error in MLwiN after retrieving REALCOM impute

When you click "Start Analysis" what should happen is that messages related to reading in the various data files will be displayed in the status bar at the bottom of the window. If these files are small or reading is complete before the status bar updates then you may not see these however...