从DGP自回归过程生成观测值？R或Stata

我想知道如何从AR进程yt=0.2+0.5yt-1+et生成数据 式中，et~N0,1，其中y1=0.2

我需要从t2到t100生成100个观测值

---

我需要在R或Stata中执行此操作

这将是R语言中的一个直接翻译

y <- numeric( length=100)
y[1] <- 0.2
et <- rnorm(100)
for( t in 2:100)
  { y[t] = 0.2 + 0.5*y[t-1] +et[t] }

---

嗨，BondedDust，我在R中尝试了该代码，但它似乎不起作用。我的R技能有点生疏，你有什么建议吗？我的错。纠正了三个语法错误。现在应该成功了。

plot(y)