对于Loop OLS回归（Matlab）

Question

我有一个50x1000矩阵，其中包含来自AR（1）过程的随机生成的数据。基本上是1000个样本，每个样本50个。我必须估计每个样本的系数为50,1000倍。我必须在其1个时间段滞后时对每个列进行回归。我不知道如何使用滞后值设置for循环来执行此操作。任何人都可以建议我这样做吗？这就是我到目前为止所拥有的。

rho=0.9;

u = randn(1000,1000);
y = zeros(1000,1000);
for i=2:1000
  y(i) = rho*y(i-1) + u(i);
  for j=1:1000
      y(i,j)=rho*y(i-1,j)+u(i,j);
end;
end;

m=y(951:1000,1:1000);

此代码生成我的数据，m矩阵是蒙特卡罗过程的最后50个观察，我必须估计rho。所以m这里是51x1000，我必须为每个列在其1个时间段滞后上回归该列做OLS。我很感激帮助。

提前致谢。