obtaining level 2 residuals from multilevel ordinal model
Posted: Mon Jan 06, 2020 10:47 pm
Hi all,
I am trying to obtain the level 2 residuals from a null two-level ordinal model using runmlwin.
These models work fine:
However, when I add in the code for residuals as I have done with multilevel logistic models with runmlwin, I get the following error:
Is it possible to get the level 2 residuals with multilevel ordinal models using runmlwin? I looked at the sample code for runmlwin, and interestingly, this seemed to be the only model that didn't add the residuals(u) line in. Thanks in advance!
I am trying to obtain the level 2 residuals from a null two-level ordinal model using runmlwin.
These models work fine:
Code: Select all
* Fit model using IGLS MQL1
runmlwin hiv cons, ///
level2(group: (cons, contrast(1/2))) ///
level1(idcode:) ///
discrete(dist(multinomial) link(ologit) denom(cons) basecategory(0)) ///
nopause
* Fit model using MCMC
runmlwin hiv cons , ///
level2(group: (cons, contrast(1/2))) ///
level1(idcode:) ///
discrete(dist(multinomial) link(ologit) denom(cons) basecategory(0)) ///
mcmc(orthogonal) initsprevious ///
nopause
Code: Select all
* Fit model using IGLS MQL1
runmlwin hiv cons, ///
level2(group: (cons, contrast(1/2))) ///
level1(idcode:) ///
discrete(dist(multinomial) link(ologit) denom(cons) basecategory(0)) ///
nopause
* Fit model using MCMC
runmlwin hiv cons , ///
level2(group: (cons, residuals(u) contrast(1/2))) ///
level1(idcode:) ///
discrete(dist(multinomial) link(ologit) denom(cons) basecategory(0)) ///
mcmc(orthogonal) initsprevious ///
nopause
Cannot specify reset or residual options at level 1, as no random part predictor variables are specified
r(198);