* This program find the H matrix given the L matrix.
* Get the Reduced Design Matrix and add the Lmatix at the bottom of the Reduced Design Matrix.
* That is add additional rows.
* Remove the Factor columns.
* NOTE NOTE: change n_pop to the correct value of number of populations - 5 rows down.

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

* enter in the number of populations.
compute n_pop = 5.

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

compute nr1 = nr- n_pop.
* print nr1.

compute p = make(n_pop,nc,0).
* print p.

compute h = make(nr1,nc,0).
* print h.

loop ii = 1 to n_pop.
loop jj = 1 to nc.
compute p(ii,jj) = w(ii,jj).
end loop.
end loop.

loop ii = 1 to nr1.
loop jj = 1 to nc.
compute h(ii,jj) = w(ii+n_pop,jj).
end loop.
end loop.

 print p.
print h.

compute h = t(H*t(p)*ginv(p*t(p))).

print h.

* compute h1 = t(hy)*p.
* print h1.

end matrix.