如题，其实就是想问如果不用这个函数，自己写程序来模拟生成一列符合arma(p,q)模型的数据，是怎么生成的。我主要是对产生时间序列的逻辑有些疑惑，求大神解答。

例如：y(t)=a*y(t-1)+error(t)-b*error(t-1)

下面是我猜想的模拟过程
error=rnorm(1000,mean=0,sd=1)


for (i in 1:1000){
if i=1
y(i)=error（i）
else
y(i)=a*y(i-1)+error(i)-b*(y(i-1)-y(target))
}


我主要有两个疑问，标红的第一个地方，表示这列数据的初始值，是直接白噪声生成的么？
标红的第二个地方表示arma模型中的error(t-1)的部分，是不是像我写的这样来的？还是说应该直接写error(i-1)
R是像我这样产生模拟数据的么？