EWMA在R中的波动性

我想从N（0,1）分布生成500个观察值（返回）。设置种子（100）以获得相同的随机数序列。 然后使用EWMA方法计算这些系列的波动率：

λ=0.94，其中表示t-j期间的回报率和滞后回报率的平均值。 绘制获得的波动率序列。 我知道我应该得到一个像下面链接中那样的图

[1] ：！[问题图像+图表]

这是我到目前为止的代码——主要是生成随机数和创建新变量——我开始尝试使用“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}

---

当你尝试编程时，你在哪里卡住了？请显示您的代码。@fortranner-我在帖子中添加了我的代码-正如您所看到的，我正在苦苦挣扎。我真正知道的就是如何生成随机数并创建新的变量。对不起！如果我无能为力，别担心