dm 'log;clear;output;clear;'; options ps=50 ls=70 pageno=1; goptions reset=global border ftext=swiss gunit=cm htext=0.3 htitle=0.5; goptions display noprompt; *********************************************************************; ** **; ** AUTHOR: Chris Bilder **; ** COURSE: STAT 4043 **; ** DATE: 9-22-00 **; ** UPDATE: **; ** PURPOSE: Basic matrix calculations in PROC IML **; ** NOTES: **; ** **; *********************************************************************; title1 'Chris Bilder, STAT 4043'; proc iml; A = {1 2 3, 4 5 6}; B = {-1 10 -1, 5 5 8}; A_trans = t(A); AplusB = A+B; AminusB = A-B; print A B A_trans AplusB AminusB; *Create a new B; B={3 0, 1 2, 0 1}; C=A*B; D=B*A; print A B C D; *redefine A; A = {1 2, 3 4}; A_inv = inv(A); check = A*inv(A); print A_inv check; quit; title2 'PROC REG output for the GPA example'; proc reg data=set1; model Y=X / XPX; run; quit;