st: Return code 199 after redoing the example in stata reference

[jjc.li@utoronto.ca]

Hi,

Since my own codes don't work, I copyed the codes from Stata reference  
I-p, page471, example 5.

I use the same codes, same data set. I only modified two things.  
First, I changed -version 10- to -version 10.1-, which is the version  
I am using. Second, I rename lnexp to lnm in the data set, because the  
variable is named lnm in their program codes. I am quite sure the  
other codes are exactly the same as the one in the book.

However, when I typed the conmand in the last line, -nlsur- bla bla,  
there the error:

nlsuraids returned 199
verify that nlsuraidss is a function evaluator program
r(199);


That's so wireld! I list my code following:



program nlsuraids

        version 10.1

        syntax varlist(min=8 max=8) if, at(name)

        tokenize `varlist'
        args w1 w2 w3 lnp1 lnp2 lnp3 lnp4 lnm

        tempname a1 a2 a3 a4
        sca1ar `a1'=`at'[1,1]
        sca1ar `a2'=`at'[1,2]
        sca1ar `a3'=`at'[1,3]
        sca1ar `a4'= 1-`a1'-`a2'-`a3'

        tempname b1 b2 b3
        sca1ar `b1'=`at'[1,4]
        sca1ar `b2'=`at'[1,5]
        sca1ar `b3'=`at'[1,6]


        tempname g11 g12 g13 g14
        tempname g21 g22 g23 g24
        tempname g31 g32 g33 g34
        tempname g41 g42 g43 g44
        sca1ar `g11'=`at'[1,7]
        sca1ar `g12'=`at'[1,8]
        sca1ar `g13'=`at'[1,9]
        sca1ar `g14'=-`g11'-`g12'-`g13'

        sca1ar `g21'=`g12'
        sca1ar `g22'=`at'[1,10]
        sca1ar `g23'=`at'[1,11]
        sca1ar `g24'=-`g21'-`g22'-`g23'

        sca1ar `g31'=`g13'
        sca1ar `g32'=`g23'
        sca1ar `g33'=`at'[1,12]
        sca1ar `g34'=-`g31'-`g32'-`g33'

        sca1ar `gel'=`g14'
        sca1ar `g42'=`g24'
        sca1ar `g43'=`g34'
        sca1ar `g44'=-`g41'-`g42'-`g43'


        quietly {
                tempvar lnpindex
                gen double `lnpindex' =
5+`a1'*`lnp1'+`a2'*`lnp2'+`a3'*`lnp3'+`a4'*`lnp4'
                forva1ues i=1/4 {
                        forva1ues j = 1/4 {
                                replace  
`lnpindex'=`lnpindex'+0.5*`g`i'`j''*`lnp`i''*`lnp`j''
                        }
                }

                replace `w1' =  
`a1'+`g11'*`lnp1'+`g12*`lnp2'+`g13'*`lnp3'+`g14'*`lnp4'+`b1'*(`lnm'-`lnpindex')

                replace `w2' =  
`a2'+`g21'*`lnp1'+`g22*`lnp2'+`g23'*`lnp3'+`g24'*`lnp4'+`b2'*(`lnm'-`lnpindex')


                replace `w3' =  
`a3'+`g31'*`lnp1'+`g32*`lnp2'+`g33'*`lnp3'+`g34'*`lnp4'+`b3'*(`lnm'-`lnpindex')



        }

end

nlsur aids @ w1 w2 w3 lnp1 lnp2 lnp3 lnp4 lnm, parameters (a1 a2 a3 b1  
b2 b3 g11 g12 g13 g22 g32 g33) neq(3) ifgnls



Jingjing


引用 Kit Baum <baum@bc.edu>：

> <>
> According to help nlsur what you have typed below is not the proper
> format of the command.
>
>
>         nlsur func_prog @ depvar_1 depvar_2 ...  [varlist] [if] [in]
> [weight] , nequations(#)
>                        {parameters(namelist)|nparameters(#)} [options]
>
> so you should list the depvars for the several (4?) equations, then
> give a varlist of variables to be passed to the routine, then comma,
> nequations() option. In what you have written below you have a ) where
> you should not and are missing the comma. Also, either parameters() or
> nparameters() MUST be given.
>
> Look at the examples at the foot of the help file for the simpler
> example of a function evaluator program, and make sure your command
> follows the appropriate syntax.
>
> Kit Baum   |   Boston College Economics and DIW Berlin   |
> http://ideas.repec.org/e/pba1.html
> An Introduction to Stata Programming   |
> http://www.stata-press.com/books/isp.html
> An Introduction to Modern Econometrics Using Stata   |
> http://www.stata-press.com/books/imeus.html
>
>
>
> On Mar 14, 2009, at 02:33 ,Jingling wrote:
>
> > nlsur wellll @ lnc sl sk sm lnpl lnpk lnpm lnpe lnq halflnqlnq lnpllnq
> > lnpklnq lnpmlnq lnpelnq lnqt lnqd1 lnqd2 lnqd3 lnqd4 lnqd5 lnqd6 t
> > halftt lnplt lnpkt lnpmt lnpet d1 d2 d3 d4 d5 d6 td1 td2 td3 td4 td5
> > td6 lnpld1 lnpld2 lnpld3 lnpld4 lnpld5 lnpld6 lnpkd1 lnpkd2 lnpkd3
> > lnpkd4 lnpkd5 lnpkd6 lnpmd1 lnpmd2 lnpmd3 lnpmd4 lnpmd5 lnpmd6 lnped1
> > lnped2 lnped3 lnped4 lnped5 lnped6, parameters(al ak am aq aqq blq bkq
> > bmq bqt bq1 bq2 bq3 bq4 bq5 bq6 at att blt bkt bmt d1 d2 d3 d4 d5 d6
> > bt1 bt2 bt3 bt4 bt5 bt6 bl1 bl2 bl3 bl4 bl5 bl6 bk1 bk2 bk3 bk4 bk5
> > bk6  bm1 bm2 bm3 bm4 bm5 bm6 be1 be2 be3 be4 be5 dkl dkk dml dmk dmm)
> > neq(4) ifgnls
>
> *
> *   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/




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