是推个没有价格但是包括了些家庭属性的需要系统
capture program drop nlsurqaidss
program nlsurqaidss
version 10
syntax varlist(min=17  max=17) [if], at(name)
tokenize`varlist'
args w1 w2 w3 w4 w5 w6 w7 w8 w9 w10 w11 lnexp z1 z2 z3 z4 z5
tempname d1 d2 d3 d4 d5 d6 d7 d8 d9 d10 d11
scalar `d1'=`at'[1,1]
scalar `d2'=`at'[1,2]
scalar `d3'=`at'[1,3]
scalar `d4'=`at'[1,4]
scalar `d5'=`at'[1,5]
scalar `d6'=`at'[1,6]
scalar `d7'=`at'[1,7]
scalar `d8'=`at'[1,8]
scalar `d9'=`at'[1,9]
scalar `d10'=`at'[1,10]
scalar `d11'=`at'[1,11]
tempname a0 a01 a02 a03 a04 a05
scalar `a0'=`at'[1,12]
scalar `a01'=`at'[1,13]
scalar `a02'=`at'[1,14]
scalar `a03'=`at'[1,15]
scalar `a04'=`at'[1,16]
scalar `a05'=`at'[1,17]
tempname a11 a12 a13 a14 a15
scalar `a11'=`at'[1,18]
scalar `a12'=`at'[1,19]
scalar `a13'=`at'[1,20]
scalar `a14'=`at'[1,21]
scalar `a15'=`at'[1,22]
tempname a21 a22 a23 a24 a25
scalar `a21'=`at'[1,23]
scalar `a22'=`at'[1,24]
scalar `a23'=`at'[1,25]
scalar `a24'=`at'[1,26]
scalar `a25'=`at'[1,27]
tempname a31 a32 a33 a34 a35
scalar`a31'=`at'[1,28]
scalar`a32'=`at'[1,29]
scalar`a33'=`at'[1,30]
scalar`a34'=`at'[1,31]
scalar`a35'=`at'[1,32]
tempname a41 a42 a43 a44 a45
scalar`a41'=`at'[1,33]
scalar`a42'=`at'[1,34]
scalar`a43'=`at'[1,35]
scalar`a44'=`at'[1,36]
scalar`a45'=`at'[1,37]
tempname a51 a52 a53 a54 a55
scalar`a51'=`at'[1,38]
scalar`a52'=`at'[1,39]
scalar`a53'=`at'[1,40]
scalar`a54'=`at'[1,41]
scalar`a55'=`at'[1,42]
tempname a61 a62 a63 a64 a65
scalar`a61'=`at'[1,43]
scalar`a62'=`at'[1,44]
scalar`a63'=`at'[1,45]
scalar`a64'=`at'[1,46]
scalar`a65'=`at'[1,47]
tempname a71 a72 a73 a74 a75
scalar`a71'=`at'[1,48]
scalar`a72'=`at'[1,49]
scalar`a73'=`at'[1,50]
scalar`a74'=`at'[1,51]
scalar`a75'=`at'[1,52]
tempname a81 a82 a83 a84 a85
scalar`a81'=`at'[1,53]
scalar`a82'=`at'[1,54]
scalar`a83'=`at'[1,55]
scalar`a84'=`at'[1,56]
scalar`a85'=`at'[1,57]
tempname a91 a92 a93 a94 a95
scalar`a91'=`at'[1,58]
scalar`a92'=`at'[1,59]
scalar`a93'=`at'[1,60]
scalar`a94'=`at'[1,61]
scalar`a95'=`at'[1,62]
tempname a101 a102 a103 a104 a105
scalar`a101'=`at'[1,63]
scalar`a102'=`at'[1,64]
scalar`a103'=`at'[1,65]
scalar`a104'=`at'[1,66]
scalar`a105'=`at'[1,67]
tempname a111 a112 a113 a114 a115
scalar`a111'=`at'[1,68]
scalar`a112'=`at'[1,69]
scalar`a113'=`at'[1,70]
scalar`a114'=`at'[1,71]
scalar`a115'=`at'[1,72]
tempname b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11
scalar `b1'=`at'[1,73]
scalar `b2'=`at'[1,74]
scalar `b3'=`at'[1,75]
scalar `b4'=`at'[1,76]
scalar `b5'=`at'[1,77]
scalar `b6'=`at'[1,78]
scalar `b7'=`at'[1,79]
scalar `b8'=`at'[1,80]
scalar `b9'=`at'[1,81]
scalar `b10'=`at'[1,82]
scalar `b11'=`at'[1,83]
tempname c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11
scalar `c1'=`at'[1,84]
scalar `c2'=`at'[1,85]
scalar `c3'=`at'[1,86]
scalar `c4'=`at'[1,87]
scalar `c5'=`at'[1,88]
scalar `c6'=`at'[1,89]
scalar `c7'=`at'[1,90]
scalar `c8'=`at'[1,91]
scalar `c9'=`at'[1,92]
scalar `c10'=`at'[1,93]
scalar `c11'=`at'[1,94]
replace `w1'=`d1'+`a11'*`z1'+`a12'*`z2'+`a13'*`z3'+`a14'*`z4'+`a15'*`z5'+`b1'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c1'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w2'=`d2'+`a21'*`z1'+`a22'*`z2'+`a23'*`z3'+`a24'*`z4'+`a25'*`z5'+`b2'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c2'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w3'=`d3'+`a31'*`z1'+`a32'*`z2'+`a33'*`z3'+`a34'*`z4'+`a35'*`z5'+`b3'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c3'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w4'=`d4'+`a41'*`z1'+`a42'*`z2'+`a43'*`z3'+`a44'*`z4'+`a45'*`z5'+`b4'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c4'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w5'=`d5'+`a51'*`z1'+`a52'*`z2'+`a53'*`z3'+`a54'*`z4'+`a55'*`z5'+`b5'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c5'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w6'=`d6'+`a61'*`z1'+`a62'*`z2'+`a63'*`z3'+`a64'*`z4'+`a65'*`z5'+`b6'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c6'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w7'=`d7'+`a71'*`z1'+`a72'*`z2'+`a73'*`z3'+`a74'*`z4'+`a75'*`z5'+`b7'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c7'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w8'=`d8'+`a81'*`z1'+`a82'*`z2'+`a83'*`z3'+`a84'*`z4'+`a85'*`z5'+`b8'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c8'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w9'=`d9'+`a91'*`z1'+`a92'*`z2'+`a93'*`z3'+`a94'*`z4'+`a95'*`z5'+`b9'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c9'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w10'=`d10'+`a101'*`z1'+`a102'*`z2'+`a103'*`z3'+`a104'*`z4'+`a105'*`z5'+`b10'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c10'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
replace `w11'=`d11'+`a111'*`z1'+`a112'*`z2'+`a113'*`z3'+`a114'*`z4'+`a115'*`z5'+`b11'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')+`c11'*(`lnexp'-`a0'-`a01'*`z1'-`a02'*`z2'-`a03'*`z3'-`a04'*`z4'-`a05'*`z5')^2
end

结果stata出来这个err

. nlsur qaidss @  w1  w2  w3  w4  w5  w6  w7  w8  w9  w10  w11  z1  z2  z3  z4
> z5  lnexp , parameters(d1 d2 d3 d4 d5 d6 d7 d8 d9 d10 d11 a0 a01 a02 a03 a04
> a05 a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a
> 44 a45 a51 a52 a53 a54 a55 a61 a62 a63 a64 a65 a71 a72 a73 a74 a75 a81 a82 a8
> 3 a84 a85 a91 a92 a93 a94 a95 a101 a102 a103 a104 a105 a111 a112 a113 a114 a1
> 15 b1b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11) nequ
> ations(11) ifgnls
nlsurqaidss returned 199
verify that nlsurqaidss is a function evaluator program
r(199);

求指导，不知道哪里写错了