Запустите y с процессом AR 1-го порядка в OLS

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

Первая часть, создание х в порядке:

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