Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down at the end of May, and its replacement, statalist.org is already up and running.


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

st: to hijack the stata ml?


From   Sun Yutao <yutao.sun.statalist@outlook.com>
To   "Statalist" <statalist@hsphsun2.harvard.edu>
Subject   st: to hijack the stata ml?
Date   Tue, 16 Oct 2012 02:58:32 +0200

Hi all,

I’m wondering if there is a way to “hijack” the stata ml command, sorry for using this word but I have a very big panel probit with fixed effects and lots of cross-sections to maximize (I know the incidental parameter problem and let’s don’t talk about that right now)

What I’m thinking is a partitioned mle which I either have to set constraints on some of the parameter in every iteration (e.g. the leapfrog) or to use a different approach of updating the parameter vector (e.g. Chamberlain 1980)

Both can be done, of course, but with problems: For the first method, the number of constraints are 1999 in stata se 12. I can put constraints in a matrix though, but it can be a huge matrix and it’s time consuming to operate it. For the second method, I can use mata to bypass the stata ml routine but then I think I will lose a lot features from the stata ml command, which I don’t really like.

I'm wondering if there is a way to "hijack" (sorry for this ugly word again :) ) the stata ml so it computes the derivatives and hessian, and updates the parameter vector in my own way but does everything else conventionally.

It will also be appreciated if someone can provide me a detailed document (really detailed) on the stata ml ado package, because I tried to read it but it's really complicated.



Best regards,
Sun Yutao




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


© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index