全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
4040 7
2013-05-08
本人第一次用stata   求教高手  
附件列表
QQ截图20130508230534.png

原图尺寸 4.05 KB

QQ截图20130508230534.png

二维码

扫码加我 拉你入群

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

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

全部回复
2013-5-8 23:46:08
你把错误提示给截个图啊
二维码

扫码加我 拉你入群

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

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

2013-5-9 00:05:57
孙悟空奥巴马 发表于 2013-5-8 23:46
你把错误提示给截个图啊
nlsuraids @ w1 w2 w3  lnp1  lnp2  lnp3  lnp4  lnexp,  parameters(a1 a2 a3 b1 b2 b3 g11 g12 g13 g22 g32 g33) ne(3) ifg
> nls (obs=4048)
@ invalid name
r(198);

我比着手册上来做的 呀
手册上是这样写的: use http://www.stata-press.com/data/r11/food, clear
                     . nlsur aids @ w1 w2 w3 lnp1 lnp2 lnp3 lnp4 lnexp,
                       > parameters(a1 a2 a3 b1 b2 b3
                       > g11 g12 g13 g22 g32 g33)
                       > neq(3) ifgnls
                        (obs = 4048)
附件列表
QQ截图20130508230534.png

原图尺寸 4.05 KB

QQ截图20130508230534.png

QQ截图20130508235749.png

原图尺寸 13.61 KB

QQ截图20130508235749.png

二维码

扫码加我 拉你入群

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

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

2013-5-9 04:58:17
手册前面还有一段很长的程序你没有执行
只执行后面的许多都不认的
aids 是没有现成程序
所以需要先执行前面编写好的aids的程序,然后才能执行后面这个。

最好先把stata的基本都先学学,再来做这些。
二维码

扫码加我 拉你入群

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

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

2013-5-9 08:41:48
蓝色 发表于 2013-5-9 04:58
手册前面还有一段很长的程序你没有执行
只执行后面的许多都不认的
aids 是没有现成程序
aids程序 我按手册里的例子编写的 写好后也保存了 才执行这个的 想不到还有错误
二维码

扫码加我 拉你入群

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

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

2013-5-9 09:15:40
clear
program drop _all

program nlsuraids
version 12
syntax varlist(min=8 max=8) if, at(name)
tokenize `varlist'
args w1 w2 w3 lnp1 lnp2 lnp3 lnp4 lnm
tempname a1 a2 a3 a4
scalar `a1' = `at'[1,1]
scalar `a2' = `at'[1,2]
scalar `a3' = `at'[1,3]
scalar `a4' = 1 - `a1' - `a2' - `a3'
tempname b1 b2 b3
scalar `b1' = `at'[1,4]
scalar `b2' = `at'[1,5]
scalar `b3' = `at'[1,6]
tempname g11 g12 g13 g14
tempname g21 g22 g23 g24
tempname g31 g32 g33 g34
tempname g41 g42 g43 g44
scalar `g11' = `at'[1,7]
scalar `g12' = `at'[1,8]
scalar `g13' = `at'[1,9]
scalar `g14' = -`g11'-`g12'-`g13'
scalar `g21' = `g12'
scalar `g22' = `at'[1,10]
scalar `g23' = `at'[1,11]
scalar `g24' = -`g21'-`g22'-`g23'
scalar `g31' = `g13'
scalar `g32' = `g23'
scalar `g33' = `at'[1,12]
scalar `g34' = -`g31'-`g32'-`g33'
scalar `g41' = `g14'
scalar `g42' = `g24'
scalar `g43' = `g34'
scalar `g44' = -`g41'-`g42'-`g43'
quietly {
tempvar lnpindex
gen double `lnpindex' = 5 + `a1'*`lnp1' + `a2'*`lnp2' + ///
`a3'*`lnp3' + `a4'*`lnp4'
forvalues i = 1/4 {
forvalues j = 1/4 {
replace `lnpindex' = `lnpindex' + ///
0.5*`g`i'`j''*`lnp`i''*`lnp`j''
}
}
replace `w1' = `a1' + `g11'*`lnp1' + `g12'*`lnp2' + ///
`g13'*`lnp3' + `g14'*`lnp4' + ///
`b1'*(`lnm' - `lnpindex')
replace `w2' = `a2' + `g21'*`lnp1' + `g22'*`lnp2' + ///
`g23'*`lnp3' + `g24'*`lnp4' + ///
`b2'*(`lnm' - `lnpindex')
replace `w3' = `a3' + `g31'*`lnp1' + `g32'*`lnp2' + ///
`g33'*`lnp3' + `g34'*`lnp4' + ///
`b3'*(`lnm' - `lnpindex')
}
end


use http://www.stata-press.com/data/r12/food
nlsur aids @ w1 w2 w3 lnp1 lnp2 lnp3 lnp4 lnexp,     ///
           parameters(a1 a2 a3 b1 b2 b3 g11 g12 g13 g22 g32 g33) ///
               neq(3) ifgnls
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群