Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
From | Richard Williams <richardwilliams.ndu@gmail.com> |
To | statalist@hsphsun2.harvard.edu, statalist@hsphsun2.harvard.edu |
Subject | Re: st: "Can Your Results be Replicated?" (Stata error?) |
Date | Thu, 19 Sep 2013 20:03:28 -0400 |
At 02:56 PM 9/19/2013, Jeff Pitblado, StataCorp LP wrote:
As for Richard's experience with perfect predictors, we believe you can expect one of the following outcomes when dealing with perfect predictors: 1. The model will fail to converge because of the lack of a concave Hessian. 2. The model "converged", but there are some extreme valued coefficients with just as extreme standard errors. Sometimes, but not always, the fitted intercept will be extreme too. For -mlogit- with the -rrr- option, these extreme values tend to translate to zero relative-risk ratios with missing standard errors.
Along those lines -- if you run the following code webuse repair, clear hetprob foreign b3.repair, het(b3.repair) oglm foreign b3.repair, het(b3.repair) link(p) loghetprob drags on for 16,000 iterations before giving an incorrect answer. oglm, on the other hand, gives an incorrect answer almost immediately. Which I guess means that oglm is better. ;-)
Anyway, you may want to add hetprob to the todo list.
--Jeff jpitblado@stata.com * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/
------------------------------------------- Richard Williams, Notre Dame Dept of Sociology OFFICE: (574)631-6668, (574)631-6463 HOME: (574)289-5227 EMAIL: Richard.A.Williams.5@ND.Edu WWW: http://www.nd.edu/~rwilliam * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/