monte carlo simulation编程求助！！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在通过montecarlo simulation做VaR预测。window size是250天，第251天的预测可以做出来，
mu=mean(whole[1:250,1])
> sigma=sd(whole[1:250,1])
> e=array(20)
> q=array()
> for(k in 1:100){
+ e=rnorm(20)
+ for(i in 2:20){
+ p[1]=whole[250,1]*(1+mu/400+sigma*e[1]/20)
+ p=p[i-1]*(1+mu/400+sigma*e/20)}
+ q[k]=p[20]}
> s=mean(q)
但是如果想做rolling window，加多一个循环，就出不来答案了，全都是NA
s=array()
for(m in 1:400)
+ {
+ mu=mean(whole[m:m+249,1])
+ sigma=sd(whole[m:m+249,1])
+ e=array(20)
+ q=array()
+ for(k in 1:1000){
+ e=rnorm(20)
+ for(i in 2:20){
+ p[1]=whole[m+249,1]*(1+mu/400+sigma*e[1]/20)
+ p=p[i-1]*(1+mu/400+sigma*e/20)}
+ q[k]=p[20]}
+ s[m]=mean(q)}
请大家帮忙找一下问题在哪里好吗？ 谢谢！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Monte Carlo Simulation ulation Carlo ATION window sigma whole

弱弱的问一下whole函数在那个包里，我怎么没找到