Re: st: Statsby equivalent with -ml-

["E. Paul Wileyto" <epw@mail.med.upenn.edu>]

This seems to run when I do it:

. statsby _b _se, by(runnum) clear: ml model lf cure_lik (hazard: _t= 
drug ) (nocure:
> _d= drug ) (shape: _t0=), max difficult iter(50)
(running ml on estimation sample)

command: ml model lf cure_lik (hazard: _t= drug ) (nocure:_d= drug ) (shape:
_t0=), max difficult iter(50)
by: runnum

Statsby groups
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
........

What error are you getting?

Paul



Nuno wrote:

> Hi everyone,
>
> I have a maximum likelihood model that I want to apply to a set of
> difference indviduals and record the estimated coefficients. Normally I use
> the  -statsby- command to record the coefficients on a dta file, but it
> seems not to work with ml (or at lesat I was unable to do so). The setup is
> the following:
>
>
> program testml
> ...
> end
>
> ml model lf testml (mu1: Y=X1) (mu2: Y=X2) /sigma1 /sigma2 ml maximize
>
>
>
> My data is organized in the following way:
> id	Y	X1	X2
>
>
> and I need to run the testml program for each id and record the coefficients
> in a .dta file.
>
> Any ideas of how to solve this issue?
>
> Best,
>
> Nuno
>
> *
> *   For searches and help try:
> *   http://www.stata.com/support/faqs/res/findit.html
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>
>   
--
E. Paul Wileyto, Ph.D.
Assistant Professor of Biostatistics
Tobacco Use Research Center
School of Medicine, U. of Pennsylvania
3535 Market Street, Suite 4100
Philadelphia, PA  19104-3309

215-746-7147
Fax: 215-746-7140
epw@mail.med.upenn.edu
http://mail.med.upenn.edu/~epw/ 
*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/