mixed - Stata
mixed - Stata
mixed - Stata
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>mixed</strong> — Multilevel <strong>mixed</strong>-effects linear regression 7<br />
organize them by model level. Residual-variance parameter estimates are also displayed in their<br />
original estimation metric.<br />
noheader suppresses the output header, either at estimation or upon replay.<br />
nogroup suppresses the display of group summary information (number of groups, average group<br />
size, minimum, and maximum) from the output header.<br />
nostderr prevents <strong>mixed</strong> from calculating standard errors for the estimated random-effects parameters,<br />
although standard errors are still provided for the fixed-effects parameters. Specifying this option<br />
will speed up computation times. nostderr is available only when residuals are modeled as<br />
independent with constant variance.<br />
nolrtest prevents <strong>mixed</strong> from fitting a reference linear regression model and using this model to<br />
calculate a likelihood-ratio test comparing the <strong>mixed</strong> model to ordinary regression. This option<br />
may also be specified on replay to suppress this test from the output.<br />
display options: noomitted, vsquish, noemptycells, baselevels, allbaselevels, nofvlabel,<br />
fvwrap(#), fvwrapon(style), cformat(% fmt), pformat(% fmt), sformat(% fmt), and<br />
nolstretch; see [R] estimation options.<br />
✄<br />
✄ <br />
EM options<br />
<br />
These options control the expectation-maximization (EM) iterations that take place before estimation<br />
switches to a gradient-based method. When residuals are modeled as independent with constant<br />
variance, EM will either converge to the solution or bring parameter estimates close to the solution.<br />
For other residual structures or for weighted estimation, EM is used to obtain starting values.<br />
✄<br />
emiterate(#) specifies the number of EM iterations to perform. The default is emiterate(20).<br />
emtolerance(#) specifies the convergence tolerance for the EM algorithm. The default is<br />
emtolerance(1e-10). EM iterations will be halted once the log (restricted) likelihood changes<br />
by a relative amount less than #. At that point, optimization switches to a gradient-based method,<br />
unless emonly is specified, in which case maximization stops.<br />
emonly specifies that the likelihood be maximized exclusively using EM. The advantage of specifying<br />
emonly is that EM iterations are typically much faster than those for gradient-based methods.<br />
The disadvantages are that EM iterations can be slow to converge (if at all) and that EM provides<br />
no facility for estimating standard errors for the random-effects parameters. emonly is available<br />
only with unweighted estimation and when residuals are modeled as independent with constant<br />
variance.<br />
emlog specifies that the EM iteration log be shown. The EM iteration log is, by default, not<br />
displayed unless the emonly option is specified.<br />
emdots specifies that the EM iterations be shown as dots. This option can be convenient because<br />
the EM algorithm may require many iterations to converge.<br />
✄<br />
Maximization<br />
<br />
maximize options: difficult, technique(algorithm spec), iterate(#), [ no ] log, trace,<br />
gradient, showstep, hessian, showtolerance, tolerance(#), ltolerance(#),<br />
nrtolerance(#), and nonrtolerance; see [R] maximize. Those that require special mention<br />
for <strong>mixed</strong> are listed below.<br />
For the technique() option, the default is technique(nr). The bhhh algorithm may not be<br />
specified.<br />
matsqrt (the default), during optimization, parameterizes variance components by using the matrix<br />
square roots of the variance–covariance matrices formed by these components at each model level.