www.cmm.bristol.ac.uk/forum

got it to work, thank you so much for your help!

Actually, I can get credible intervals around a group-specific effect using this code:


bysort group (iteration): egen p_mean = mean(p)
bysort group (iteration): egen p_lo = pctile(p), p(2.5)
bysort group (iteration): egen p_hi = pctile(p), p(97.5)


where p is a variable for the group specific ...

thank you so much! so that makes sense, although it looks like I can only use mcmcsum (with the variables option) if I use all the chains; when I try to use it on a subset (e.g., on the residual values for one group) I get an output that looks like this:

. mcmcsum u if group==1, variables ...

To give you some more context: After I fit the models and pull residuals:


runmlwin outcome cons , ///
level2(group: cons) ///
level1(id:) ///
discrete(distribution(binomial) link(logit) denominator(denominator) pql2) ///
rigls maxiterations(100) ///
nopause

runmlwin outcome cons ...

Thank you for your reply! I have tried mcmcsum, and I get the following error:


Estimate results must come from runmlwin


Furthermore, when I load the .dta file where the parameter chains I stored, the only variables included are: iteration, deviance, FP1_cons, RP2_var_cons_, and OD_bcons_1; I ...

Hello! I am fitting a null two-level multilevel logistic model using the following example code:



runmlwin outcome cons , ///
level2(group: cons) ///
level1(id:) ///
discrete(distribution(binomial) link(logit) denominator(denominator) pql2) ///
rigls maxiterations(100) ///
nopause ...

Hello,

I am trying to fit a two-level multilevel logistic model in runmlwin using survey weights. The model runs fine using PQL2 estimation methods, however, I get the following error when trying to run the same model using MCMC methods: "Weights are only valid for univariate models estimated ...

yes, thank you so much!! it's always the comma placement that gets me

Hi all,

I am trying to obtain the level 2 residuals from a null two-level ordinal model using runmlwin.

These models work fine:


* Fit model using IGLS MQL1
runmlwin hiv cons, ///
level2(group: (cons, contrast(1/2))) ///
level1(idcode:) ///
discrete(dist(multinomial) link(ologit) denom(cons ...

Hi Chris,

Oh yes, I see now where I can attach a file. Thanks!
I created the script as you suggested, and got the following error:


./mlnscript: error while loading shared libraries: libmln.so: cannot open shared object file: No such file or directory


However, that file (libmln.so) is in ...