r语言 - 在 OLS 中使用一阶自回归过程创建 x

Question

我有一个简单的回归：yt=β1+β2xi+ei，n=27，"x"是AR（1）：

xi = c + ∅

x（i-1） + ηi ， 其中 ηi~N（0,1） ， x0~N（c/（1-∅），1/（1-∅^2） ， c=2 ， ∅=0.6

我需要创建"x"，为此我已经设置了包括"x0"在内的所有内容，但是我卡住了：

phi <- 0.6
c <- 2
ni <- rnorm(27)
x0 <- rnorm(1,(c/(1-phi)),(1/(1-(phi)^2)))

有了for我无法让它工作：

xa<- vector(mode="numeric",length=0)
x<- rep(0,27)
for(i in 1:27){
  w<-c+phi*x[i-1]+ni[i]
  xa <- c(xa,w)
}

我应该如何创建"x"？

Answer 1

看起来你计算的 x[i-1] 不用于计算 x[i]。 如果我了解您的AR模型，您可以尝试

x <- rep(0,27)
x[1] <- x0
for(i in 2:27){
  x[i] <-c + phi*x[i-1] + ni[i]
}

或使用R函数生成更通用的AR模拟

mu <- c/(1-phi)
x <- arima.sim(model=list(ar=c(phi,0), sd=1), n=27) + mu

更新

要让 x[1] 包含 x0 之后的第一个计算值，请尝试

x <- rep(0,28)
x[1] <- x0
for(i in 2:28){
  x[i] <-c + phi*x[i-1] + ni[i-1]
}
# to change indexing so x[1] is not x0 but is first computed value
for( i in 1:27) { x[i] <- x[i+1] }
x <- x[-28]

Answer 2

好吧，我设法提出了一个解决方案：

设置后

phi <- 0.6
c <- 2
e <- as.vector(rnorm(27))
ni<- as.vector(rnorm(27))
x0 <- rnorm(1,(c/(1-phi)),(1/(1-(phi)^2)))

我们需要做这样的for：

x<- rep(0,27)
x[1]<- c+phi*x0+ni[1]
for(i in 2:27){
  x[i]<-c+phi*x[i-1]+ni[i]
}

然后它起作用了!!!