--- Benjamin Allaire <[email protected]> wrote:
> I get this bizarre error when I try to predict the outcomes of my
> multinomial logit.
>
> It is an 18 cell outcome variable.
>
> It goes something like this:
>
> svy: mlogit mnlcat age agesq age3, baseoutcome(3)
>
> (output omitted)
>
> estimates store est1
>
> (other program stuff)
>
> use sample_pop, clear
>
> estimates restore est1
> foreach X of numlist 1/18 {
> predictnl p`X' = predict(outcome(`X')), se(p_se`X')
> }
>
> Maximum number of iterations exceeded.
> r(498);
The maximum iterations refer to iterations that take place while
computing a numerical derivative within -predictnl-. The easiest way to
solve this is to compute these probabilities and their standard errors
with -prgen- within the -spost- package: see -findit spost-, and:
Jun Xu and J. Scott Long (2005) "Confidence intervals for predicted
outcomes in regression models for categorical outcomes" The Stata
Journal, 5(4): 537--559.
http://www.stata-journal.com/article.html?article=st0094
Hope this helps,
Maarten
-----------------------------------------
Maarten L. Buis
Department of Social Research Methodology
Vrije Universiteit Amsterdam
Boelelaan 1081
1081 HV Amsterdam
The Netherlands
visiting address:
Buitenveldertselaan 3 (Metropolitan), room Z434
+31 20 5986715
http://home.fsw.vu.nl/m.buis/
-----------------------------------------
Send instant messages to your online friends http://uk.messenger.yahoo.com
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
