Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down on April 23, and its replacement, **statalist.org** is already up and running.

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

From |
"Alba J. Collart-Dinarte" <ACollart@ag.tamu.edu> |

To |
"statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu> |

Subject |
st: MLE for Tobit I |

Date |
Mon, 24 Oct 2011 16:27:18 +0000 |

Dear Stata listers, I have a question about programming the likelihood function for Tobit models (type I, censored at zero) in STATA 10.1. This is the code I've written: program mytobit1 version 10.1 args lnf beta sigma quietly replace `lnf'= log(1-normal(`beta'/`sigma')) if $ML_y1==0 quietly replace `lnf'= log((1/`sigma')*normalden(($ML_y1-`beta')/`sigma')) if $ML_y1>0 end ml model lf mytobit1 (y1 = x) /sigma ml check ml maximize, nolog but it does not run. I believe that there is a problem in the way I define the likelihood function. I have looked into Gould, Pitblado and Poi (2010), but they mainly discuss probit, normal and weibull models. I have also tried: quietly replace `lnf'= cond($ML_y1==0,log(1-normal(`beta',`sigma')),log((1/`sigma')*normalden($ML_y1,`beta',`sigma'))) and using normalden($ML_y1,`beta',`sigma') - instead of - normalden(($ML_y1-`beta')/`sigma') but no success. Any suggestions on how can I fix the code would be greatly appreciated. Thank you very much in advance! Alba Collart * * 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/

**Follow-Ups**:**Re: st: MLE for Tobit I***From:*Matthew Baker <matthew.baker@hunter.cuny.edu>

- Prev by Date:
**Re: st: Fwd: Comparing marginal effects of two subsamples** - Next by Date:
**re:st: FE over two columns** - Previous by thread:
**st: FE over two columns** - Next by thread:
**Re: st: MLE for Tobit I** - Index(es):