Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

RE: st: predictions from ML program


From   Maarten buis <[email protected]>
To   [email protected]
Subject   RE: st: predictions from ML program
Date   Wed, 11 Jun 2008 09:26:37 +0100 (BST)

--- marcel spijkerman <[email protected]> wrote:
> Predict is indeed useful but what I need is the contribution to the
> likelihood for each of the four terms in the likelihood (in:
> 
> ln(p1*L1 + p2*L2 + p3*L3 + (1-p1-p2-p3)*L4)
> 
> The model is a bivariate survival model with discrete heterogeneity
> distribution; and the terms in the likelihood are contributions for
> each person of the two duration processes. Can I obtain the separate
> contributions L1, L2, L3 and L4 by using predict or other
> postestimation commands?

These elements will be different equations in your likelihood, so you
can access them separately using -predict- with the -eq()- option. 

The Stata -ml- facilities are great, and make your live as a programmer
a lot easier, except when you are comming from another language. In
that case there are a lot of peculiarities you'll have to learn and
unlearn before you can really make the best use of these facilities. A
great place to start is: http://www.stata.com/bookstore/mle.html .

Hope this helps,
Maarten




-----------------------------------------
Maarten L. Buis
Department of Social Research Methodology
Vrije Universiteit Amsterdam
Boelelaan 1081
1081 HV Amsterdam
The Netherlands

visiting address:
Buitenveldertselaan 3 (Metropolitan), room Z434

+31 20 5986715

http://home.fsw.vu.nl/m.buis/
-----------------------------------------


      __________________________________________________________
Sent from Yahoo! Mail.
A Smarter Email http://uk.docs.yahoo.com/nowyoucan.html
*
*   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/



© Copyright 1996–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index