使用指定的非零均值和 AR 系数在 R 中模拟 AR(1) 过程

机器算法验证 r 时间序列 有马 模拟 自回归的
2022-04-12 16:06:55

我需要在 R 中使用以下等式模拟 AR(1) 过程:

Xt=5+0.5Xt1+Zt

在哪里Zt〜白噪声(0,1)和T=500.

我知道我应该使用带有参数arima.sim的包中的函数但我不知道如何解释forecastn = 500ar = c(0.5)δ=5或者事实上Zt是平均值为 0 和标准差为 1 的白噪声。我一直无法找到有关如何执行此操作的明确文档。

1个回答

请注意,这是相同的模型

Xt10=.5(Xt110)+Zt.
10是平均值,而5是拦截。这意味着我们可以将 10 添加到平均零序列。

换句话说,如果你定义Yt=Xt10, 然后

Yt=.5Yt1+Zt
arima.sim用来模拟它。在你拥有之后{Yt}, 然后加10到每个值。

代码如下所示:

yt <- arima.sim(list(order=c(1,0,0), ar=.5), n=500)
xt <- yt + 10   

关于标准偏差问题,文档中对此进行了介绍。