Matrices

[arbongard]

proc iml;

    a = {1 2, 3 4};

    c= {1 0, 0 1};

     b = a || c;

    d = echelon(b);

    x1 = d[1,3];

    x2= d[1,4] ;

    x3 = d[2,3];

    x4 = d[2,4];

    e = {0 0, 0 0};

    e[1,1] = x1;

    e[1,2] = x2;

    e[2,1]= x3;

    e[2,2]= x4;

print d, e;

Is the and easier way to put grab the take the 2x2 matrix out of the 2 by 4 matrix?

[IanWakeling]

See the section here on "extracting rows and columns".

http://support.sas.com/publishing/pubcat/chaps/63119.pdf