在OLS中使用一阶AR流程运行y

Question

我有一个回归：* yt =β1+β2* xi + ei *，n = 23，“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即可：

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

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

如何使用此回归创建y： * yt =β1+β2* xi + ei *

且β1= 0.913且β2= 0.015