* This program finds the L matrix given the H matrix.
* Get the Reduced Design Matrix and add the H matrix to the right of the  Reduced Design Matrix.
* That is add additional columns.
* Remove the Factor columns.
* NOTE NOTE: change df  to the correct value of number of df  for h  - 5 rows down.



MATRIX.
GET W/VARIABLES = all/FILE=*/missing=omit.
* print w.

* enter in the degrees of freedom for h

compute df= 1.

compute nr = nrow(w).
compute NC = ncol(w).
* print nr.
* print nc.

compute nc1 = nc- df.
* print nr1.

compute p = make(nr,nc1,0).
* print p.

compute h = make(nr,df,0).
* print h.

loop ii = 1 to nr.
loop jj = 1 to nc1.
compute p(ii,jj) = w(ii,jj).
end loop.
end loop.

loop ii = 1 to nr.
loop jj = 1 to df.
compute h(ii,jj) = w(ii,jj+nc1).
end loop.
end loop.

 print p.
print h.

compute L = t(h)*p.

print L.

end matrix.