我想从 N(0,1) 分布中生成 500 个观察值(返回)。设置种子(100)以获得相同的随机数序列。然后使用 EWMA 方法计算这些系列的波动率:
其中 λ=0.94,其中表示 tj 期收益的均值和滞后收益。绘制获得的波动率序列。我知道我应该得到如下链接中的图表
[1]:![问题图片+图表] http://imgur.com/mSiUZxA
到目前为止,这是我的代码——主要是生成随机数和创建新变量——我开始尝试使用“for”命令,但不知道该去哪里:
set.seed(100)
data=rnorm(n=500,mean=0,sd=1)
lambda=0.94
rbar=mean(data)
dsquared=NULL
for (j in 0:100){dsquared[j]=lambda^j*(data[j]-rbar)^2}