st: maximum likelihood estimation similar to Poi 2002

[Alex Olssen <alex.olssen@gmail.com>]

Dear Statalist,

I am trying to manually estimate a system similar to the QUAIDS system
estimated in Poi's article in The Stata Journal, 2002, 2 (4).

A crucial difference is that I do not what to impose symmetry on my
coefficients like Poi did.

However my code doesn't seem to quite work.  I have fixed a couple of
problems but now when I try the command

ml search

I get the error

"b not found"

I suspected I had not put `' around a local somewhere but have checked
several times for this.

I was wondering if anybody would be interested enough to quickly look
through my code - altogether it is less than 80 lines long and this
includes spacing and comments for readability.
It makes a call to the program vec_sum.ado which was written by Poi
and can be accessed through the command

net sj 2-4
net install st0029



Finally, I am aware that in The Stata Journal, 2008, 8 (4), that Poi
provided another way to estimate QUAIDS using nlsur.

However I am particularly interested in the log-likelihood value
associated with my estimation and do not know of away to recover it
following nlsur.

Kind regards,

Alex Olssen
Motu Research
Wellington
New Zealand
*
*   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/