Dear Gauri,
Your program (mylog) solves a probit model. Is this what
you want? I don't understand what loglinear means for you.
If the answer is yes, I assume that your dependent variable
is 0 or 1 and I suggest to use probit instead.
Otherwise, check -glm- command take a loof of options
family() and link() maybe gaussian and log is what you need
Or provide the model that you want to estimate.
Rodrigo.
----- Original Message -----
From: "Deepankar Basu" <basu.15@osu.edu>
To: <statalist@hsphsun2.harvard.edu>
Sent: Thursday, July 06, 2006 4:37 PM
Subject: Re: st: ML estimation, no feasible values-Gauri, 6th July
I think the only think that is missing is another empty bracket for
-theta1- in your -ml model- command. Try the following:
ml model lf mylog (outputpsp=corareapsp labourpsp manurepsp
fertiliserpsp tracpsp oxpsp irripsp) ()
Nothing else needs to be changed.
HTH,
Deepankar
On Thu, 2006-07-06 at 20:18 +0000, Gauri Khanna wrote:
> Dear Statalist,
>
> I am trying to run an ml maximisation to calculate loglikehood values for
> my
> model which has only one equation. It is a loglinear equation.
> unfortunately
> I run into problems as no feasible values are found, error message r(491).
> I
> even tried the ml search, repeat(10) and also with ml maximize,difficult
> and
> also with ml maximize, ltol(1e-4) but I got the same error message r(491)
> that I have pasted below. I only want to obtain the loglikelihood value
> estimated, nothing else!
>
> MY COMMANDS:
>
> capture program drop mylog
> program mylog
> version 8
> args lnf theta1
> quietly replace `lnf' = ln(norm(`theta1')) if $ML_y1==1
> quietly replace `lnf' = ln(norm(-`theta1')) if $ML_y1==0
> end
>
> ml model lf mylog (outputpsp=corareapsp labourpsp manurepsp fertiliserpsp
> tracpsp oxpsp irripsp)
> ml check
> ml search
> ml maximize
>
>
> OUTPUT GENERATED IN STATA
>
> . capture program drop mylog
>
> . program mylog
> 1. version 8
> 2. args lnf theta1
> 3. quietly replace `lnf' = ln(norm(`theta1')) if $ML_y1==1
> 4. quietly replace `lnf' = ln(norm(-`theta1')) if $ML_y1==0
> 5. end
>
> .
> . ml model lf mylog (outputpsp=corareapsp labourpsp manurepsp
> fertiliserpsp
> tracpsp oxpsp
> > irripsp)
>
> . ml check
>
> Test 1: Calling mylog to check if it computes log likelihood and
> does not alter coefficient vector...
> Passed.
>
> Test 2: Calling mylog again to check if the same log likelihood value
> is returned...
> Passed.
>
> ------------------------------------------------------------------------------
> The initial values are not feasible. This may be because the initial
> values
> have been chosen poorly or because there is an error in mylog and it
> always returns missing no matter what the parameter values.
>
> Stata is going to use ml search to find a feasible set of initial values.
> If mylog is broken, this will not work and you will have to press Break
> to make ml search stop.
>
> Searching...
> initial: log likelihood = -<inf> (could not be evaluated)
> searching for feasible values
> ...........................................................
> >
> ......................................................................
> >
> .......................................................................................
> >
> .......................................................................................
> >...........................................................................
>
> could not find feasible values
> r(491);
>
>
> Please HELP, what can I do with my model "mylog" to make it find feasible
> values?
>
> Thanking you in advance.
>
> Gauri
>
> _________________________________________________________________
> Don't just search. Find. Check out the new MSN Search!
> http://search.msn.click-url.com/go/onm00200636ave/direct/01/
>
> *
> * 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/
*
* 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/
*
* 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/