matrix.

* L to H.

*NOTE NOTE: Remove all columns form the design matrices except the p columns.

* P1 is the design matrix for the population parameter model.
* P2 is the design matrix for the computer parameter model.
* l is the matrix of coefficients for what is being tested in the computer model  .
* h is the matrix of coefficients for what is being tested in the population parameter model.


* Note: save the design matrices as C:\p1.sav  c:\P2.sav.

* Note:  save L as C:\l.sav.

* Note the output is h ans it is saved as C:\h.sav.


get p1 /file = 'c:\p1.sav'.
get p2 /file = 'c:\p2.sav'.
get l1 /file ='C:\l.sav'.

compute l = t(l1).

print l.

compute d = ginv(T(p2)*p2)*T(p2)*p1.

compute h = l*d.

print h.

save h /outfile = 'c:\h.sav'.
end matrix.