matrix.

* H to L.

*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 H as C:\H.sav.

* Note the output is L and it is saved as C:\L.sav.


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

print h.

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

compute l = h*d.

print l.

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