请问谁知道错在哪？ stata程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

program nlsurxy
version 11
syntax varlist(min=10 max=10) if, at(name)
tokenize `varlist'
args w1 w2 w3 w4 lnp1 lnp2 lnp3 lnp4 lnp5 lnm
tempname a1 a2 a3 a4 a5
scalar `a1' = `at'[1,1]
scalar `a2' = `at'[1,2]
scalar `a3' = `at'[1,3]
scalar `a4' = `at'[1,4]
scalar `a5' = 1 - `a1' - `a2' - `a3'-`a4'
tempname b1 b2 b3 b4
scalar `b1' = `at'[1,5]
scalar `b2' = `at'[1,6]
scalar `b3' = `at'[1,7]
scalar `b4' = `at'[1,8]
tempname g11 g12 g13 g14 g15
tempname g21 g22 g23 g24 g25
tempname g31 g32 g33 g34 g35
tempname g41 g42 g43 g44 g45
tempname g51 g52 g53 g54 g55
scalar `g11' = `at'[1,9]
scalar `g12' = `at'[1,10]
scalar `g13' = `at'[1,11]
scalar `g14' = `at'[1,12]
scalar `g15' = -`g11'-`g12'-`g13'-`g14'
scalar `g21' = `g12'
scalar `g22' = `at'[1,13]
scalar `g23' = `at'[1,14]
scalar `g24' = `at'[1,15]
scalar `g25' = -`g21'-`g22'-`g23'-`g24'
scalar `g31' = `g13'
scalar `g32' = `g23'
scalar `g33' = `at'[1,16]
scalar `g34' = `at'[1,17]
scalar `g35' = -`g31'-`g32'-`g33'-`g34'

scalar `g41' = `g14'
scalar `g42' = `g24'
scalar `g43' = `g34'
scalar `g44' = `at'[1,18]
scalar `g45' = -`g41'-`g42'-`g43'-`g44'
scalar `g51' = `g15'
scalar `g52' = `g25'
scalar `g53' = `g35'
scalar `g54' = `g45'
scalar `g55' = -`g51'-`g52'-`g53'-`g54'
quietly {
tempvar lnpindex
gen double `lnpindex' = 5 + `a1'*`lnp1' + `a2'*`lnp2' + ///
`a3'*`lnp3' + `a4'*`lnp4'+`a5'*`lnp5'
forvalues i = 1/5 {
forvalues j = 1/5 {
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'  ///
`b1'*(`lnm' - `lnpindex')
replace `w2' = `a2' + `g21'*`lnp1' + `g22'*`lnp2' + ///
`g23'*`lnp3' + `g24'*`lnp4' + `g25'*`lnp5' ///
`b2'*(`lnm' - `lnpindex')
replace `w3' = `a3' + `g31'*`lnp1' + `g32'*`lnp2' + ///
`g33'*`lnp3' + `g34'*`lnp4' +`g35'*`lnp5'  ///
`b3'*(`lnm' - `lnpindex')
replace `w4' = `a4' + `g41'*`lnp1' + `g42'*`lnp2' + ///
`g43'*`lnp3' + `g44'*`lnp4' +`g45'*`lnp5'  ///
`b4'*(`lnm' - `lnpindex')

}
end


use "E:\资料\课题\AIDS\AIDS1.dta", clear
nlsur xy @ w1 w2 w3 w4 lnp1 lnp2 lnp3 lnp4 lnp5 lnexp,     ///
           parameters(a1 a2 a3 a4 b1 b2 b3 b4 g11 g12 g13 g14 g22 g32 g33 g42 g43 g44 ) ///
               neq(4) ifgnls

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：stata程序 Stata tata parameters forvalues 程序

这。。。你又不给点数据，又不给出错误信息，又不好好排版。。。

不好意思，忘了 发数据。  附上数据

没看明白。。

错误提示

夏目贵志 发表于 2013-11-5 09:33
这。。。你又不给点数据，又不给出错误信息，又不好好排版。。。

，数据已上传了 谢谢

LilyChanserena 发表于 2013-11-5 09:33

错哪了 给点建议呗

巨蟹小漫 发表于 2013-11-5 09:46
没看明白。。

有语法错误吗？

或与非 发表于 2013-11-5 10:06
有语法错误吗？

http://blog.stata.com/tag/nl/
http://statalist.1588530.n2.nabble.com/st-Re-retruned-r-198-in-nlsur-td2478353.html