wolffd@0: function yy = convert_to_lagged_form(y, k)
wolffd@0: % Create an observation vector yy(:,t) containing the last k values of y, newest first
wolffd@0: % e.g., k=2, y = (a1 a2 a3)     yy  = a2 a3
wolffd@0: %                (b1 b2 b3)           b2 b2
wolffd@0: %                                     a1 a2
wolffd@0: %                                     b1 b2
wolffd@0: 
wolffd@0: [s T] = size(y);
wolffd@0: bs = s*ones(1,k);
wolffd@0: yy = zeros(k*s, T-k+1);
wolffd@0: for i=1:k
wolffd@0:   yy(block(i,bs), :) = y(:, k-i+1:end-i+1);
wolffd@0: end
wolffd@0: