蒙特卡洛模拟程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这是我从王群勇上照搬的蒙特卡洛模拟程序，怎么不出结果，问题在哪里？谢谢

program rdwalk5,rclass
version 10
syntax [,n(integer 1000)
drop _all
quietly set obs 'n'
capture drop y x1 x2 x3 x4 x5 rndy rndx1 rndx2 rndx3 rndx4 rndx5
quietly {
    gen rndy=100+(1000-100)*invnorm(uniform())
    gen rndx1=10+(100-10)*invnorm(uniform())
    gen rndx2=9+(22-9)*invnorm(uniform())
    gen rndx3=2+(12-2)*invnorm(uniform())
    gen rndx4=6+(24-6)*invnorm(uniform())
    gen rndx5=0+(100-0)*invnorm(uniform())
    gen y=rndy
    gen x1=rndx1
    gen x2=rndx2
    gen x3=rndx3
    gen x4=rndx4
    gen x5=rndx5
    replace y=y[_n-1]+rndy[-n] in 2/'n'
    replace x1=x1[_n-1]+rndx1[-n] in 2/'n'
    replace x3=x3[_n-1]+rndx2[-n] in 2/'n'
    replace x4=x4[_n-1]+rndx3[-n] in 2/'n'
    replace x5=x5[_n-1]+rndx4[-n] in 2/'n'
    regress y x1 x2 x3 x4 x5
    matrix matb=e(b)
    matrix matv=e(V)
    return scalar b=el(matb,1,2)
    return scalar se=sqrt(el(matv,2,2))
  }
end
simulate beta=(r(b)) se=(r(se)),reps(10000) nodots:rdwalk

扫码加我 拉你入群

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

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

分享0 收藏5 回帖

关键词：蒙特卡洛模拟 模拟程序 蒙特卡洛 蒙特卡 simulate 蒙特卡洛 模拟程序

模拟程序Stata里simulate命令帮助中就有

program define lnsim, rclass
        version 10.1
        syntax [, obs(integer 1) mu(real 0) sigma(real 1) ]
        drop _all
        set obs `obs'
        tempvar z
        gen `z' = exp(rnormal(`mu',`sigma'))
        summarize `z'
        return scalar mean = r(mean)
        return scalar Var  = r(Var)
end
simulate mean=r(mean) var=r(Var), reps(10000): lnsim, obs(100)

谢谢，这个我也试过了。那么，我再问一个菜鸟的问题：

syntax [, obs(integer 1) mu(real 0) sigma(real 1) ]

能具体给我解释一下各项的含义吗？谢谢！

以下是引用wangjie20012006在2008-9-21 8:06:00的发言：

谢谢，这个我也试过了。那么，我再问一个菜鸟的问题：

syntax [, obs(integer 1) mu(real 0) sigma(real 1) ]

能具体给我解释一下各项的含义吗？谢谢！

一言难尽，参见help syntax。

已经解决了，谢谢！

2# eblog

为什么我用这个程度执行不了？结果显示
Unknown function rnormal()
r(133);

请博主赐教!谢谢

程序中有几个明显的小错误
1.第三行 “syntax [,n(integer 1000) ”换成 syntax [,n(integer 1000)]
2.程序中的所有的'n'换成 `n'，左边的撇是Esc下面的那个符号
3.所有的[-n]换成[_n]
另：修改后的程序，我跑了一遍没有问题

rnd函数需要自己安装,可能是LZ没有装

怎么安装rnd函数？谢谢！

好东西。。