我想从N（0,1）分布中生成500个观测值（返回值）。设置种子（100）以获得相同的随机数序列。 然后使用EWMA方法计算这些序列的波动率：
λ= 0.94，其中表示回报的均值和时段t-j的滞后回报。 绘制获得的波动率系列。 我知道我应该得到一个如下图所示的图表
[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}