我需要一个矩阵形式的OLS循环

Question

这是我的回归：Yt= B0 + B1*t + B2*Yt-1

如果我写我的矩阵：

n=1
round(solve(t(X[,c(1, n + 1, n + 4)])%*%X[,c(1, n + 1, n + 4)])%*%t(X[,c(1, n + 1, n + 4)])%*%Y[,n], digits=4)

它告诉我这个

#      [,1]
#[1,]  0.0920
#[2,] -0.0007
#[3,]  0.4000

对于n = 1是正确的，但是当n = 1,2,3

时我需要一个代码

直到现在我都试过了：

bh<-array (0,dim=(c(3,20,1)))
for (n in 1:3)
bh[,n] = round(solve(t(X[,c(1, n + 1, n + 4)])%*%X[,c(1, n + 1, n + 4)])%*%t(X[,c(1, n + 1, n + 4)])%*%Y[,n], digits=4)

但它说：

  

矩阵中的订阅者数量不正确