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);