在宏观经济预测最近的网上课程,有模型的锻炼; Tibial模拟时间序列
y(t) = 3.0 + 0.55 y(t-1) + e(t)
其中e(t)
被定义为
et <- c(-1.2138662, -0.2854597, 0.5902700, 0.8285463, -0.9954260, -0.3716332)
现在我试图做到这一点的R(该课程使用了EViews),但是我没有找到给出的解决方案:第五个元素为5.648。 我想(与此类似blogpost):
y <- rep(NA,6)
y[1] <- 0
y[2] <- 3 + 0.55*y[1]+et[1]
y[3] <- 3 + 0.55*y[2]+et[2]
y[4] <- 3 + 0.55*y[3]+et[3]
y[5] <- 3 + 0.55*y[4]+et[4]
y[6] <- 3 + 0.55*y[5]+et[5]
然后
y <- rep(NA,6)
y[1] <- et[1]
y[2] <- 3 + 0.55*y[1]+et[2]
y[3] <- 3 + 0.55*y[2]+et[3]
y[4] <- 3 + 0.55*y[3]+et[4]
y[5] <- 3 + 0.55*y[4]+et[5]
y[6] <- 3 + 0.55*y[5]+et[6]
然后
arima.sim(list(order=c(1,0,0), ar=0.55), n=6, innov=head(et,6)+3)
然而所有这三种方法得出不同的结果。我想知道这是为什么,恐怕我不了解一些基本的东西。
是的我同意前两种方法必须导致不同的结果。然而,三种方法都没有给出y [5.6] = 5.648,这是练习中给出的解。 –
@Khashaa:这是(解决方案被舍入为三个数字)?你是怎么得到这个的? –
你的第一种方法有错字。在一行中检查'et [i]'项 – Khashaa