全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
4104 4
2012-04-22
set obs 100
set seed 500//这两个有什么用的?
prog monte
gen x1=5*runiform()
gen x2=rnormal(5,2)
gen err=invnormal(uniform())
gen y=5+x1+0.5*x2+err
reg y x1 x2
end
simulate _b _cons, reps(1000): monte//


这样会得到下面这样的结果。。
        +-------------------------------------------+
        _b_x1      _b_x2    _b_cons   _eq2_s~1
        -------------------------------------------
1.        1.104805   .5876165   4.224726          1
2.        .          .          .          .
3.        .          .          .          .
剩下全是空白。。

A.要取得100次/1000次参数(beta1 beta2 cons)的估计值,怎么做?
B.取得1000次参数估计值后,要画在图上与beta真值进行比较(line _b_x1 1???),怎么画?纵坐标是数值,横坐标是实验的次数,在同一幅图上还要画reg y x1的beta估计值,怎么画?
二维码

扫码加我 拉你入群

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

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

全部回复
2012-4-22 02:15:25
剩下全是空白
cap pr drop monte
pr monte
clear
set ob 100
g x1=5*runiform()
g x2=rnormal(5, 2)
g err=rnormal()
g y=5+x1+0.5*x2+err
reg y x1 x2
end
simulate, r(1000) nodots seed(500): monte
二维码

扫码加我 拉你入群

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

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

2012-4-22 02:23:25
要画在图上与beta真值进行比较,怎么画?纵坐标是数值,横坐标是实验的次数,在同一幅图上还要画reg y x1的beta估计值,怎么画?
g b1=1
g b2=.5
g bc=5
g n=_n
sc _b* n || line b* n
二维码

扫码加我 拉你入群

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

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

2012-4-22 12:51:10
sungmoo 发表于 2012-4-22 02:23
g b1=1
g b2=.5
g bc=5
谢谢!!!
二维码

扫码加我 拉你入群

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

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

2014-5-17 22:36:13
学习了
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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