[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
KONSTANTARAS KONSTANTINOS <dinokon@otenet.gr> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re:Re: st: GLLAMM multinomial: tremendous instability |

Date |
Wed, 3 Dec 2008 18:28:39 +0200 (EET) |

Dear statalisters and personally Stas, Maarten and J. Verkuilen: Apologies for my late reply to your helpful suggestions -due to a strong (non PC) virus. The problem has been resolved by Sophia Rabe-Hesqueth a few days ago; it has been the result of reversing the order of the expand() option variables, namely patt and chosen. With the correct order, instability has stopped. As a matter of fact, the random effects correlation equal to 1 problem has not ceased to exist and I would be much obliged if you could post me any references you might have come accross on this. Much obliged, Dino Konstantaras .................................................................................................................................... On 11/24/08, Stas Kolenikov [skolenik@gmail.com] wrote: Dino, Give your -mlogit- results, too; did you feed them as starting values to -gllamm-? Keep in mind that introducing random effects changes the scaling of your coefficients; see http://www.citeulike.org/user/ctacmo/article/3057661. I'd fully agree with earlier reactions that your model does not appear to be identified well. Correlation of 1 are troublesome; in your particular case this may mean that there is no difference between two of the three alternatives -- is this making sense? If one of the categories is quite rare, or present in some weird patterns wrt units producing random effects, it may create those numeric problems, too. On 11/24/08, KONSTANTARAS KONSTANTINOS <dinokon@otenet.gr> wrote:

Dear statalisters,In my model runs (using Intercooled Stata 9.2) I experiencetremendous instability in gllamm results, higher than has beenelsewhere reported atstatalist: I run a simple multinomial logit model with random effects,using the results as initial matrix for exactly the same model, to getas a result an extremely different log likelihood and coefficientsetc. The resulting loglikelihood differs by large margins (3-9%). Ihave also tried to use exactly the same initial values from mlogit butnot only does the overall loglikelihood continues to be unstable butalso the same happens with the loglikelihood obtained without randomeffects –running gllamm with init option-. Nevertheless the modelinvariably converges, albeit to an entirely different estimate.

* * 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/

**Follow-Ups**:**RE: Re: st: GLLAMM multinomial: tremendous instability***From:*"Verkuilen, Jay" <JVerkuilen@gc.cuny.edu>

- Prev by Date:
**Re: st: New version of -rsource- on SSC** - Next by Date:
**st: running a backward stepwise multivariate analysis on my specific dataset** - Previous by thread:
**st: Agreement statistic, gold standard not binary** - Next by thread:
**RE: Re: st: GLLAMM multinomial: tremendous instability** - Index(es):

© Copyright 1996–2016 StataCorp LP | Terms of use | Privacy | Contact us | What's new | Site index |