RE: st: ml porgramming

["emre ekinci" <eekinci2610@hotmail.com>]

sorry,I forgot the code in the previous thread.
here is the code:

capture program drop myll

program define myll
version 9


args todo b lnf

tempvar theta1 c1 c2 c3 c4 c5 c6 c7 c8 c9 a1

tempvar lnfj


mleval `theta1' = `b'

quietly gen double  `c1' = (0.628222008/(`theta1'))
quietly gen double  `c2' = (0.04836303/(`theta1'))
quietly gen double  `c3' = (0.009742825/(`theta1'))
quietly gen double  `c4' = (0.021884783/(`theta1'))
quietly gen double  `c5' = (0.242897686/(`theta1'))
quietly gen double  `c6' = (0.008396969/(`theta1'))
quietly gen double  `c7' = (0.008133649/(`theta1'))
quietly gen double `c8' = (0.016618391/(`theta1'))
quietly gen double  `c9' = (0.015740659/(`theta1'))


constraint define 1 `c1' + `c2' + `c3' = 0.688412327
constraint define 2 `c4' + `c5' + `c6' = 0.270318976
constraint define 3 `c7' + `c8' + `c9' = 0.041268697
constraint define 4 `c1' + `c4' + `c7' = 0.655623075
constraint define 5 `c2' + `c5' + `c8' = 0.311162907
constraint define 6 `c3' + `c6' + `c9' = 0.033214018


quietly gen double `lnfj' = ln(`c1') if ............
quietly replace  `lnfj' = ln(`c2') if ............
quietly replace  `lnfj' = ln(`c3') if .................
quietly replace  `lnfj' = ln(`c4') if .................
quietly replace  `lnfj' = ln(`c5') if ...............
quietly replace  `lnfj' = ln(`c6') if ..............
quietly replace  `lnfj' = ln(`c7') if...............
quietly replace  `lnfj' = ln(`c8') if ..............
quietly replace  `lnfj' = ln(`c9') if ..................

mlsum `lnf' = `lnfj'

end



> From: "emre ekinci" <eekinci2610@hotmail.com>
> Reply-To: statalist@hsphsun2.harvard.edu
> To: statalist@hsphsun2.harvard.edu
> Subject: st: ml porgramming
> Date: Fri, 18 May 2007 00:19:11 +0000
>
> dear all,
>
> I am trying to estimate a 3x3 transition matrix.each cell contains a known 
> probability as the numerator and a function g of which parameters are to be 
> estimated.the model is identified only if we impose some constraints.to 
> estimate the model,I wrote a code shown below;yet I cannot figure out how I 
> can incorporate these nonlinear constraints.the code below always gives 
> some errors-mainly error code 111.
> I really appreciate any advices.
> thank you in advance.
> best.
>
> _________________________________________________________________
> En etkili ve g�venilir PC Korumayi tercih edin, rahat edin! 
> http://www.msn.com.tr/security/
>
> *
> *   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/
_________________________________________________________________
Siz siz olun MSN'den hava durumunu �grenmeden evden �ikmayin! 
http://www.msn.com.tr/havadurumu/

*
*   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/