-------------------------------------------------------------------------------
       log:  /Applications/Stata/771ml1.smcl
  log type:  smcl
 opened on:  19 Jan 2004, 11:19:34


. . adotype fishy1_lf ~/ado/personal/fishy1_lf.ado:

program fishy1_lf version 8.2 args lnf theta quietly replace `lnf' = -`theta' + log(`theta')*$ML_y1 - lnfact($ML_y1) end

. . clear

. input y

y 1. 5 2. 0 3. 1 4. 1 5. 0 6. 3 7. 2 8. 3 9. 4 10. 1 11. end

. . ml model lf fishy1_lf (y=)

. . ml check

Test 1: Calling fishy1_lf to check if it computes log likelihood and does not alter coefficient vector... Passed.

Test 2: Calling fishy1_lf 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 fishy1_lf 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 fishy1_lf 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 + feasible: log likelihood = -31.105155 improving initial values +.....+... improve: log likelihood = -19.61942 rescaling entire vector +. rescale: log likelihood = -19.535825

restarting tests... ------------------------------------------------------------------------------

Test 1: Calling fishy1_lf to check if it computes log likelihood and does not alter coefficient vector... Passed.

Test 2: Calling fishy1_lf again to check if the same log likelihood value is returned... Passed.

Test 3: Calling fishy1_lf to check if 1st derivatives are computed... test not relevant for method lf.

Test 4: Calling fishy1_lf again to check if the same 1st derivatives are returned... test not relevant for method lf.

Test 5: Calling fishy1_lf to check if 2nd derivatives are computed... test not relevant for method lf.

Test 6: Calling fishy1_lf again to check if the same 2nd derivatives are returned... test not relevant for method lf.

------------------------------------------------------------------------------ Searching for alternate values for the coefficient vector to verify that fishy1_lf returns different results when fed a different coefficient vector:

Searching... initial: log likelihood = -<inf> (could not be evaluated) searching for feasible values +

feasible: log likelihood = -18.544773 improving initial values ......+... improve: log likelihood = -18.39317

continuing with tests... ------------------------------------------------------------------------------

Test 7: Calling fishy1_lf to check log likelihood at the new values... Passed.

Test 8: Calling fishy1_lf requesting 1st derivatives at the new values... test not relevant for method lf.

Test 9: Calling fishy1_lf requesting 2nd derivatives at the new values... test not relevant for method lf.

------------------------------------------------------------------------------ fishy1_lf HAS PASSED ALL TESTS ------------------------------------------------------------------------------

Test 10: Does fishy1_lf produce unanticipated output? This is a minor issue. Stata has been running fishy1_lf with all output suppressed. This time Stata will not suppress the output. If you see any unanticipated output, you need to place quietly in front of some of the commands in fishy1_lf.

-------------------------------------------------------------- begin execution ---------------------------------------------------------------- end execution

. . ml maximize

initial: log likelihood = -18.39317 rescale: log likelihood = -18.39317 Iteration 0: log likelihood = -18.39317 Iteration 1: log likelihood = -18.379288 Iteration 2: log likelihood = -18.379268 Iteration 3: log likelihood = -18.379268

Number of obs = 10 Wald chi2(0) = . Log likelihood = -18.379268 Prob > chi2 = .

------------------------------------------------------------------------------ y | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _cons | 2 .4472136 4.47 0.000 1.123477 2.876523 ------------------------------------------------------------------------------

. . test _cons = 1

( 1) [eq1]_cons = 1

chi2( 1) = 5.00 Prob > chi2 = 0.0253

. . log close log: /Applications/Stata/771ml1.smcl log type: smcl closed on: 19 Jan 2004, 11:19:39 -------------------------------------------------------------------------------