Re: st: Error 198 in function evaluator program with nl

[Manuel Flores <mflores06@gmail.com>]

Thanks Nick for your reply.

I made the change Nick suggests, but I keep getting the same error message.

I really need to run the program. Thanks in advance for any other help.

Manuel Flores
mflores06@gmail.com


2013/11/11 Nick Cox <njcoxstata@gmail.com>:
> The lines
>
> scalar `delta' = `at'[1,17]
>
> gen double `wterm' = ln(exp(exp({`delta'})*`v16')-1) `if'
>
> look problematic to me. The problems start because `delta' is a
> tempname. You may want something more like
>
> gen double `wterm' = ln(exp(exp({delta=`delta'})*`v16')-1) `if'
>
> although my wild guess is that other fixes are needed.
> Nick
> njcoxstata@gmail.com
>
>
> On 11 November 2013 19:11, Manuel Flores <mflores06@gmail.com> wrote:
>> Dear statalist users,
>> I'm trying to estimate a LS nonlinear model with the nl command in
>> Stata MP/12.1. I have a gravity model of trade, with almost 300
>> variables (including country fixed effects), so I need to do it
>> running a FEP. As I'm having troubles with the program I started with
>> a shorter version of the model without the dummies, keeping just 17
>> independent variables and 18 parameters.
>> However, I keep getting a 198 error code, and I can't realize where
>> the bugs are.
>> Thanks in advance for any help!
>>
>> Here is the program:
>> ------------------------------
>> capture program drop nlgrav
>> set trace on
>> program nlgrav
>>       version 12
>>       syntax varlist (min=18 max=18) [aw fw iw] if, at(name)
>>       local lexp: word 1 of `varlist'
>>       local v1: word 2 of `varlist'
>>       local v2: word 3 of `varlist'
>>   local v3: word 4 of `varlist'
>>   local v4: word 5 of `varlist'
>>   local v5: word 6 of `varlist'
>>   local v6: word 7 of `varlist'
>>   local v7: word 8 of `varlist'
>>   local v8: word 9 of `varlist'
>>   local v9: word 10 of `varlist'
>>   local v10: word 11 of `varlist'
>>   local v11: word 12 of `varlist'
>>   local v12: word 13 of `varlist'
>>   local v13: word 14 of `varlist'
>>   local v14: word 15 of `varlist'
>>   local v15: word 16 of `varlist'
>>   local v16: word 17 of `varlist'
>>   local v17: word 18 of `varlist'
>>       // Retrieve parameters out of at matrix
>>       tempname b0 b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12 b13 b14 b15 delta bun
>>       scalar `b0' = `at'[1,1]
>>   scalar `b1' = `at'[1,2]
>>   scalar `b2' = `at'[1,3]
>>   scalar `b3' = `at'[1,4]
>>   scalar `b4' = `at'[1,5]
>>   scalar `b5' = `at'[1,6]
>>   scalar `b6' = `at'[1,7]
>>   scalar `b7' = `at'[1,8]
>>   scalar `b8' = `at'[1,9]
>>   scalar `b9' = `at'[1,10]
>>   scalar `b10' = `at'[1,11]
>>   scalar `b11' = `at'[1,12]
>>   scalar `b12' = `at'[1,13]
>>   scalar `b13' = `at'[1,14]
>>   scalar `b14' = `at'[1,15]
>>   scalar `b15' = `at'[1,16]
>>   scalar `delta' = `at'[1,17]
>>       scalar `bun' = `at'[1,18]
>>       // Some temporary variables
>>       tempvar linterm1 linterm2 linterm3 wterm heckterm
>>       gen double `linterm1' = (`b0' + `b1'*`v1' + `b2'*`v2' +
>> `b3'*`v3' + `b4'*`v4' + `b5'*`v5') `if'
>>   summarize `linterm1'
>>   gen double `linterm2' = (`b6'*`v6' + `b7'*`v7' + `b8'*`v8' +
>> `b9'*`v9' + `b10'*`v10') `if'
>>   gen double `linterm3' = (`b11'*`v11' + `b12'*`v12' + `b13'*`v13' +
>> `b14'*`v14' + `b15'*`v15') `if'
>>   gen double `wterm' = ln(exp(exp({`delta'})*`v16')-1) `if'
>>   gen double `heckterm' = (`bun'*`v17') `if'
>>   // Compute the first part of the function
>>   replace `lexp' = (`linterm1' + `linterm2' + `linterm3' + `wterm' +
>> `heckterm')  `if'
>> end
>>
>> -------------------------------
>> I'm creating a matrix with initial values:
>>
>> matrix A = J(1, 18, 0)
>> matrix A[1,1]=10
>> matrix A[1,2]=1
>> matrix A[1,17]=1
>> matrix A[1,18]=1
>>
>> And running the program with the following syntax:
>>
>> nl grav @ lexp_posit  ///
>>   ldistcap ///
>>   border ///
>>   island_both ///
>>   island_to_nonisl ///
>>   island_fr_nonisl ///
>>   landlock_both ///
>>   landl_to_nonldlk     ///
>>           landl_fr_nonldlk ///
>>   comleg ///
>>   comlang_off ///
>>   from_heg_to_col ///
>>   from_col_to_heg ///
>>   comcur ///
>>   rta ///
>>   rel_coinc ///
>>   z_bar_hat ///
>>   eta_hat2  ///
>>    , parameters(b0 b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12 b13 b14 b15
>> delta bun) at(A)
>>
>>
>> Best,
>> Manuel Flores
>> *
>> *   For searches and help try:
>> *   http://www.stata.com/help.cgi?search
>> *   http://www.stata.com/support/faqs/resources/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/faqs/resources/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/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/