u1 = v(:,1); a2 = v(:,1)\v(:,2); u2 = v(:,2) - v(:,1)*a2; a3 = v(:,[1 2])\v(:,3); u3 = v(:,3) - v(:,[1 2])*a3; a4 = v(:,[1 2 3])\v(:,4); u4 = v(:,4) - v(:,[1 2 3])*a4; % Structural residuals and covariance
b02 = [ 1 0 0 0; -a2(1) 1 0 0; -a3(1) -a3(2) 1 0; -a4(1) -a4(2) -a4(3) 1 ]; u = [u1 u2 u3 u4 ];
Please advise how can I estimate a , u and b with for loop…in case I have many more variables then just four.
Best Answer