![Statalist](/includes/images/head_statalist.jpg) |
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
Re: st: Statsby equivalent with -ml-
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
[email protected]
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/