Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: Could someone help me with my nl evaluation program? I update my program.


From   张志强 <[email protected]>
To   "[email protected]" <[email protected]>
Subject   st: Could someone help me with my nl evaluation program? I update my program.
Date   Thu, 30 Jan 2014 22:51:14 +0800

Dear all
I try to estimate a nonlinear lease square function in stata,It need
to estimate a compound summation terms like
lnwage_ij=ko+k1*ln(∑(income*exp(k2*d_ij) , d_ij is the distance between
individual observation i and observation j. My problems is how to deal
with the ∑ in my nl evaluation function. My naive code are as follows,
here I just assume only three region, there would be a lots of  dyadic
distance for different region in my sample.

 the code is as follows:
program nlwage1
version 12
syntax varlist(min=5 max=5) if, at(name)
local lnwage : word 1 of `varlist'
local income : word 2 of `varlist'
local dist1: word 3 of `varlist'
local dist2: word 4 of `varlist'
local dist3:word 5 of `varlist'
tempname k0 k1 k2
scalar `k0' = `at'[1, 1]
scalar `k1' = `at'[1, 2]
scalar `k2' = `at'[1, 3]

tempvar eterm

gen double `eterm'=`income'*exp(`k2'*`dist1')+`income'*exp(`k2'*`dist2')+`income'*exp(`k2'*`dist3')
   `if'

replace `lnwage' = `k0' +`k1'*ln(`eterm')  `if'
end

nl wage1 @ lnwage income prov1 prov2 prov3, parameters(k0 k1 k2)
initial(k0 0.05 k1 0.02 k2 0.03)


it return error

nlwage1 returned 199
verify that nlwage1 is a function evaluator program

What is wrong with my stata code?
I sincerely hope that someone can help me!

Yours

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


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index