蒙特卡洛仿真实验，为什么只有第一次的参数的估计值？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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估计值，怎么画？

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：蒙特卡洛 蒙特卡 估计值 simulate Uniform 蒙特卡洛

剩下全是空白

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

要画在图上与beta真值进行比较，怎么画？纵坐标是数值，横坐标是实验的次数，在同一幅图上还要画reg y x1的beta估计值，怎么画？

g b1=1
g b2=.5
g bc=5
g n=_n
sc _b* n || line b* n

sungmoo 发表于 2012-4-22 02:23
g b1=1
g b2=.5
g bc=5

谢谢！！！

学习了