options ps=54 ls=76; FILENAME website HTTP "http://www.stat.ufl.edu/~winner/data/wnba_ht_wt.dat"; data wnba; infile website; input id height weight; run; /* Matrix Notes: A` is the transpose of A A*B does matrix multiplication inv(A) gives the inverse of A */ proc iml; use wnba; read all; Y=weight; X0=j(nrow(height),1,1); X=X0||height; * print X Y; run;