Hi,
I developed a script for a 9-good AIDS demand system based on the example given in the nlsur manual. But, when I run it it keeps keeps giving me error messages. I was searching the internet for an exact solution for the problem and trying different suggestions but without a success. I really appreciate if somebody can have a look on my script and identify where i have gone wrong. I am using STATA/IC 13.1. The error messages I am getting are like "nlsuraids returned 102 verify that nlsuraids is a function evaluator program" (Earlier I got number 199 as well). Quick response is much appreciated.
Manoj
* function evaluator program ***************************** program nlsuraids
version 13
syntax varlist(min=18 max=18) if, at(name)
tokenize `varlist'
args w1 w2 w3 w4 w5 w6 w7 w8 lnp1 lnp2 lnp3 lnp4 lnp5 lnp6 lnp7 lnp8 lnp9 lnm
tempname a1 a2 a3 a4 a5 a6 a7 a8 a9
scalar `a1' = `at'[1,1]
scalar `a2' = `at'[1,2]
scalar `a3' = `at'[1,3]
scalar `a4' = `at'[1,4]
scalar `a5' = `at'[1,5]
scalar `a6' = `at'[1,6]
scalar `a7' = `at'[1,7]
scalar `a8' = `at'[1,8]
scalar `a9' = 1 - `a1' - `a2' - `a3'- `a4' - `a5' - `a6'- `a7' - `a8'
tempname b1 b2 b3 b4 b5 b6 b7 b8
scalar `b1' = `at'[1,9]
scalar `b2' = `at'[1,10]
scalar `b3' = `at'[1,11]
scalar `b4' = `at'[1,12]
scalar `b5' = `at'[1,13]
scalar `b6' = `at'[1,14]
scalar `b7' = `at'[1,15]
scalar `b8' = `at'[1,16]
tempname g11 g12 g13 g14 g15 g16 g17 g18 g19
tempname g21 g22 g23 g24 g25 g26 g27 g28 g29
tempname g31 g32 g33 g34 g35 g36 g37 g38 g39
tempname g41 g42 g43 g44 g45 g46 g47 g48 g49
tempname g51 g52 g53 g54 g55 g56 g57 g58 g59
tempname g61 g62 g63 g64 g65 g66 g67 g68 g69
tempname g71 g72 g73 g74 g75 g76 g77 g78 g79
tempname g81 g82 g83 g84 g85 g86 g87 g88 g89
tempname g91 g92 g93 g94 g95 g96 g97 g98 g99
scalar `g11' = `at'[1,17]
scalar `g12' = `at'[1,18]
scalar `g13' = `at'[1,19]
scalar `g14' = `at'[1,20]
scalar `g15' = `at'[1,21]
scalar `g16' = `at'[1,22]
scalar `g17' = `at'[1,23]
scalar `g18' = `at'[1,24]
scalar `g19' = -`g11'-`g12'-`g13'-`g14'-`g15'-`g16'-`g17'-`g18'
scalar `g21' = `g12'
scalar `g22' = `at'[1,25]
scalar `g23' = `at'[1,26]
scalar `g24' = `at'[1,27]
scalar `g25' = `at'[1,28]
scalar `g26' = `at'[1,29]
scalar `g27' = `at'[1,30]
scalar `g28' = `at'[1,31]
scalar `g29' = -`g21'-`g22'-`g23'-`g24'-`g25'-`g26'-`g27'-`g28'
scalar `g31' = `g13'
scalar `g32' = `g23'
scalar `g33' = `at'[1,32]
scalar `g34' = `at'[1,33]
scalar `g35' = `at'[1,34]
scalar `g36' = `at'[1,35]
scalar `g37' = `at'[1,36]
scalar `g38' = `at'[1,37]
scalar `g39' = -`g31'-`g32'-`g33'-`g34'-`g35'-`g36'-`g37'-`g38'
scalar `g41' = `g14'
scalar `g42' = `g24'
scalar `g43' = `g34'
scalar `g44' = `at'[1,38]
scalar `g45' = `at'[1,39]
scalar `g46' = `at'[1,40]
scalar `g47' = `at'[1,41]
scalar `g48' = `at'[1,42]
scalar `g49' = -`g41'-`g42'-`g43'-`g44'-`g45'-`g46'-`g47'-`g48'
scalar `g51' = `g15'
scalar `g52' = `g25'
scalar `g53' = `g35'
scalar `g54' = `g45'
scalar `g55' = `at'[1,43]
scalar `g56' = `at'[1,44]
scalar `g57' = `at'[1,45]
scalar `g58' = `at'[1,46]
scalar `g59' = -`g51'-`g52'-`g53'-`g54'-`g55'-`g56'-`g57'-`g58'
scalar `g61' = `g16'
scalar `g62' = `g26'
scalar `g63' = `g36'
scalar `g64' = `g46'
scalar `g65' = `g56'
scalar `g66' = `at'[1,47]
scalar `g67' = `at'[1,48]
scalar `g68' = `at'[1,49]
scalar `g69' = -`g61'-`g62'-`g63'-`g64'-`g65'-`g66'-`g67'-`g68'
scalar `g71' = `g17'
scalar `g72' = `g27'
scalar `g73' = `g37'
scalar `g74' = `g47'
scalar `g75' = `g57'
scalar `g76' = `g67'
scalar `g77' = `at'[1,50]
scalar `g78' = `at'[1,51]
scalar `g79' = -`g71'-`g72'-`g73'-`g74'-`g75'-`g76'-`g77'-`g78'
scalar `g81' = `g18'
scalar `g82' = `g28'
scalar `g83' = `g38'
scalar `g84' = `g48'
scalar `g85' = `g58'
scalar `g86' = `g68'
scalar `g87' = `g78'
scalar `g88' = `at'[1,52]
scalar `g89' = -`g81'-`g82'-`g83'-`g84'-`g85'-`g86'-`g87'-`g88'
scalar `g91' = `g19'
scalar `g92' = `g29'
scalar `g93' = `g39'
scalar `g94' = `g49'
scalar `g95' = `g59'
scalar `g96' = `g69'
scalar `g97' = `g79'
scalar `g98' = `g89'
scalar `g99' = -`g91'-`g92'-`g93'-`g94'-`g95'-`g96'-`g97'-`g98'
quietly {
tempvar
gen double `lnpindex' = 5 + `a1'*`lnp1' + `a2'*`lnp2' + `a3'*`lnp3' + ///
`a4'*`lnp4' + `a5'*`lnp5' + `a6'*`lnp6' + ///
`a7'*`lnp7' + `a8'*`lnp8' + `a9'*`lnp9'
forvalues i = 1/9 {
forvalues j = 1/9 {
replace `lnpindex' = `lnpindex' + 0.5*`g`i'`j''*`lnp`i''*`lnp`j''
}
}
replace `w1' = `a1' + `g11'*`lnp1' + `g12'*`lnp2' + `g13'*`lnp3' +///
`g14'*`lnp4' + `g15'*`lnp5'+ `g16'*`lnp6' +///
`g17'*`lnp7' + `g18'*`lnp8'+ `g19'*`lnp9' +///
`b1'*(`lnm' - `lnpindex')
replace `w2' = `a2' + `g21'*`lnp1' + `g22'*`lnp2' + `g23'*`lnp3' +///
`g24'*`lnp4' + `g25'*`lnp5'+ `g26'*`lnp6' +///
`g27'*`lnp7' + `g28'*`lnp8'+ `g29'*`lnp9' +///
`b2'*(`lnm' - `lnpindex')
replace `w3' = `a3' + `g31'*`lnp1' + `g32'*`lnp2' + `g33'*`lnp3' +///
`g34'*`lnp4' + `g35'*`lnp5'+ `g36'*`lnp6' +///
`g37'*`lnp7' + `g38'*`lnp8'+ `g39'*`lnp9' +///
`b3'*(`lnm' - `lnpindex')
replace `w4' = `a4' + `g41'*`lnp1' + `g42'*`lnp2' + `g43'*`lnp3' +///
`g44'*`lnp4' + `g45'*`lnp5'+ `g46'*`lnp6' +///
`g47'*`lnp7' + `g48'*`lnp8'+ `g49'*`lnp9' +///
`b4'*(`lnm' - `lnpindex')
replace `w5' = `a5' + `g51'*`lnp1' + `g52'*`lnp2' + `g53'*`lnp3' +///
`g54'*`lnp4' + `g55'*`lnp5'+ `g56'*`lnp6' +///
`g57'*`lnp7' + `g58'*`lnp8'+ `g59'*`lnp9' +///
`b5'*(`lnm' - `lnpindex')
replace `w6' = `a6' + `g61'*`lnp1' + `g62'*`lnp2' + `g63'*`lnp3' +///
`g64'*`lnp4' + `g65'*`lnp5'+ `g66'*`lnp6' +///
`g67'*`lnp7' + `g68'*`lnp8'+ `g69'*`lnp9' +///
`b6'*(`lnm' - `lnpindex')
replace `w7' = `a7' + `g71'*`lnp1' + `g72'*`lnp2' + `g73'*`lnp3' +///
`g74'*`lnp4' + `g75'*`lnp5'+ `g76'*`lnp6' +///
`g77'*`lnp7' + `g78'*`lnp8'+ `g79'*`lnp9' +///
`b7'*(`lnm' - `lnpindex')
replace `w8' = `a8' + `g81'*`lnp1' + `g82'*`lnp2' + `g83'*`lnp3' +///
`g84'*`lnp4' + `g85'*`lnp5'+ `g86'*`lnp6' +///
`g87'*`lnp7' + `g88'*`lnp8'+ `g89'*`lnp9' +///
`b8'*(`lnm' - `lnpindex')
}
end
nlsur aids @ w1 w2 w3 w4 w5 w6 w7 w8 lnp1 lnp2 lnp3 lnp4 lnp5 lnp6 lnp7 lnp8 lnp9 lnm,parameters(a1 a2 a3 a4 a5 a6 a7 a8 b1 b2 b3 b4 b5 b6 b7 b8 g11 g12 g13 g14 g15 g16 g17 g18 g22 g23 g24 g25 g26 g27 g28 g33 g34 g35 g36 g37 g38 g44 g45 g46 g47 g48 g55 g56 g57 g58 g66 g67 g68 g77 g78 g88) nequations(8) ifgnls
I developed a script for a 9-good AIDS demand system based on the example given in the nlsur manual. But, when I run it it keeps keeps giving me error messages. I was searching the internet for an exact solution for the problem and trying different suggestions but without a success. I really appreciate if somebody can have a look on my script and identify where i have gone wrong. I am using STATA/IC 13.1. The error messages I am getting are like "nlsuraids returned 102 verify that nlsuraids is a function evaluator program" (Earlier I got number 199 as well). Quick response is much appreciated.
Manoj
* function evaluator program ***************************** program nlsuraids
version 13
syntax varlist(min=18 max=18) if, at(name)
tokenize `varlist'
args w1 w2 w3 w4 w5 w6 w7 w8 lnp1 lnp2 lnp3 lnp4 lnp5 lnp6 lnp7 lnp8 lnp9 lnm
tempname a1 a2 a3 a4 a5 a6 a7 a8 a9
scalar `a1' = `at'[1,1]
scalar `a2' = `at'[1,2]
scalar `a3' = `at'[1,3]
scalar `a4' = `at'[1,4]
scalar `a5' = `at'[1,5]
scalar `a6' = `at'[1,6]
scalar `a7' = `at'[1,7]
scalar `a8' = `at'[1,8]
scalar `a9' = 1 - `a1' - `a2' - `a3'- `a4' - `a5' - `a6'- `a7' - `a8'
tempname b1 b2 b3 b4 b5 b6 b7 b8
scalar `b1' = `at'[1,9]
scalar `b2' = `at'[1,10]
scalar `b3' = `at'[1,11]
scalar `b4' = `at'[1,12]
scalar `b5' = `at'[1,13]
scalar `b6' = `at'[1,14]
scalar `b7' = `at'[1,15]
scalar `b8' = `at'[1,16]
tempname g11 g12 g13 g14 g15 g16 g17 g18 g19
tempname g21 g22 g23 g24 g25 g26 g27 g28 g29
tempname g31 g32 g33 g34 g35 g36 g37 g38 g39
tempname g41 g42 g43 g44 g45 g46 g47 g48 g49
tempname g51 g52 g53 g54 g55 g56 g57 g58 g59
tempname g61 g62 g63 g64 g65 g66 g67 g68 g69
tempname g71 g72 g73 g74 g75 g76 g77 g78 g79
tempname g81 g82 g83 g84 g85 g86 g87 g88 g89
tempname g91 g92 g93 g94 g95 g96 g97 g98 g99
scalar `g11' = `at'[1,17]
scalar `g12' = `at'[1,18]
scalar `g13' = `at'[1,19]
scalar `g14' = `at'[1,20]
scalar `g15' = `at'[1,21]
scalar `g16' = `at'[1,22]
scalar `g17' = `at'[1,23]
scalar `g18' = `at'[1,24]
scalar `g19' = -`g11'-`g12'-`g13'-`g14'-`g15'-`g16'-`g17'-`g18'
scalar `g21' = `g12'
scalar `g22' = `at'[1,25]
scalar `g23' = `at'[1,26]
scalar `g24' = `at'[1,27]
scalar `g25' = `at'[1,28]
scalar `g26' = `at'[1,29]
scalar `g27' = `at'[1,30]
scalar `g28' = `at'[1,31]
scalar `g29' = -`g21'-`g22'-`g23'-`g24'-`g25'-`g26'-`g27'-`g28'
scalar `g31' = `g13'
scalar `g32' = `g23'
scalar `g33' = `at'[1,32]
scalar `g34' = `at'[1,33]
scalar `g35' = `at'[1,34]
scalar `g36' = `at'[1,35]
scalar `g37' = `at'[1,36]
scalar `g38' = `at'[1,37]
scalar `g39' = -`g31'-`g32'-`g33'-`g34'-`g35'-`g36'-`g37'-`g38'
scalar `g41' = `g14'
scalar `g42' = `g24'
scalar `g43' = `g34'
scalar `g44' = `at'[1,38]
scalar `g45' = `at'[1,39]
scalar `g46' = `at'[1,40]
scalar `g47' = `at'[1,41]
scalar `g48' = `at'[1,42]
scalar `g49' = -`g41'-`g42'-`g43'-`g44'-`g45'-`g46'-`g47'-`g48'
scalar `g51' = `g15'
scalar `g52' = `g25'
scalar `g53' = `g35'
scalar `g54' = `g45'
scalar `g55' = `at'[1,43]
scalar `g56' = `at'[1,44]
scalar `g57' = `at'[1,45]
scalar `g58' = `at'[1,46]
scalar `g59' = -`g51'-`g52'-`g53'-`g54'-`g55'-`g56'-`g57'-`g58'
scalar `g61' = `g16'
scalar `g62' = `g26'
scalar `g63' = `g36'
scalar `g64' = `g46'
scalar `g65' = `g56'
scalar `g66' = `at'[1,47]
scalar `g67' = `at'[1,48]
scalar `g68' = `at'[1,49]
scalar `g69' = -`g61'-`g62'-`g63'-`g64'-`g65'-`g66'-`g67'-`g68'
scalar `g71' = `g17'
scalar `g72' = `g27'
scalar `g73' = `g37'
scalar `g74' = `g47'
scalar `g75' = `g57'
scalar `g76' = `g67'
scalar `g77' = `at'[1,50]
scalar `g78' = `at'[1,51]
scalar `g79' = -`g71'-`g72'-`g73'-`g74'-`g75'-`g76'-`g77'-`g78'
scalar `g81' = `g18'
scalar `g82' = `g28'
scalar `g83' = `g38'
scalar `g84' = `g48'
scalar `g85' = `g58'
scalar `g86' = `g68'
scalar `g87' = `g78'
scalar `g88' = `at'[1,52]
scalar `g89' = -`g81'-`g82'-`g83'-`g84'-`g85'-`g86'-`g87'-`g88'
scalar `g91' = `g19'
scalar `g92' = `g29'
scalar `g93' = `g39'
scalar `g94' = `g49'
scalar `g95' = `g59'
scalar `g96' = `g69'
scalar `g97' = `g79'
scalar `g98' = `g89'
scalar `g99' = -`g91'-`g92'-`g93'-`g94'-`g95'-`g96'-`g97'-`g98'
quietly {
tempvar
gen double `lnpindex' = 5 + `a1'*`lnp1' + `a2'*`lnp2' + `a3'*`lnp3' + ///
`a4'*`lnp4' + `a5'*`lnp5' + `a6'*`lnp6' + ///
`a7'*`lnp7' + `a8'*`lnp8' + `a9'*`lnp9'
forvalues i = 1/9 {
forvalues j = 1/9 {
replace `lnpindex' = `lnpindex' + 0.5*`g`i'`j''*`lnp`i''*`lnp`j''
}
}
replace `w1' = `a1' + `g11'*`lnp1' + `g12'*`lnp2' + `g13'*`lnp3' +///
`g14'*`lnp4' + `g15'*`lnp5'+ `g16'*`lnp6' +///
`g17'*`lnp7' + `g18'*`lnp8'+ `g19'*`lnp9' +///
`b1'*(`lnm' - `lnpindex')
replace `w2' = `a2' + `g21'*`lnp1' + `g22'*`lnp2' + `g23'*`lnp3' +///
`g24'*`lnp4' + `g25'*`lnp5'+ `g26'*`lnp6' +///
`g27'*`lnp7' + `g28'*`lnp8'+ `g29'*`lnp9' +///
`b2'*(`lnm' - `lnpindex')
replace `w3' = `a3' + `g31'*`lnp1' + `g32'*`lnp2' + `g33'*`lnp3' +///
`g34'*`lnp4' + `g35'*`lnp5'+ `g36'*`lnp6' +///
`g37'*`lnp7' + `g38'*`lnp8'+ `g39'*`lnp9' +///
`b3'*(`lnm' - `lnpindex')
replace `w4' = `a4' + `g41'*`lnp1' + `g42'*`lnp2' + `g43'*`lnp3' +///
`g44'*`lnp4' + `g45'*`lnp5'+ `g46'*`lnp6' +///
`g47'*`lnp7' + `g48'*`lnp8'+ `g49'*`lnp9' +///
`b4'*(`lnm' - `lnpindex')
replace `w5' = `a5' + `g51'*`lnp1' + `g52'*`lnp2' + `g53'*`lnp3' +///
`g54'*`lnp4' + `g55'*`lnp5'+ `g56'*`lnp6' +///
`g57'*`lnp7' + `g58'*`lnp8'+ `g59'*`lnp9' +///
`b5'*(`lnm' - `lnpindex')
replace `w6' = `a6' + `g61'*`lnp1' + `g62'*`lnp2' + `g63'*`lnp3' +///
`g64'*`lnp4' + `g65'*`lnp5'+ `g66'*`lnp6' +///
`g67'*`lnp7' + `g68'*`lnp8'+ `g69'*`lnp9' +///
`b6'*(`lnm' - `lnpindex')
replace `w7' = `a7' + `g71'*`lnp1' + `g72'*`lnp2' + `g73'*`lnp3' +///
`g74'*`lnp4' + `g75'*`lnp5'+ `g76'*`lnp6' +///
`g77'*`lnp7' + `g78'*`lnp8'+ `g79'*`lnp9' +///
`b7'*(`lnm' - `lnpindex')
replace `w8' = `a8' + `g81'*`lnp1' + `g82'*`lnp2' + `g83'*`lnp3' +///
`g84'*`lnp4' + `g85'*`lnp5'+ `g86'*`lnp6' +///
`g87'*`lnp7' + `g88'*`lnp8'+ `g89'*`lnp9' +///
`b8'*(`lnm' - `lnpindex')
}
end
nlsur aids @ w1 w2 w3 w4 w5 w6 w7 w8 lnp1 lnp2 lnp3 lnp4 lnp5 lnp6 lnp7 lnp8 lnp9 lnm,parameters(a1 a2 a3 a4 a5 a6 a7 a8 b1 b2 b3 b4 b5 b6 b7 b8 g11 g12 g13 g14 g15 g16 g17 g18 g22 g23 g24 g25 g26 g27 g28 g33 g34 g35 g36 g37 g38 g44 g45 g46 g47 g48 g55 g56 g57 g58 g66 g67 g68 g77 g78 g88) nequations(8) ifgnls
Comment