proc iml;

use Sasuser.paving;

Read all var {x1  } into x;

print x;

nr = nrow(x);

Print nr;

sum1 = 0;

do i = 1 to nr;

sum1 = sum1 + x[i,1];

end;

mean = sum1/nr;

print mean;

sum2 = 0;

	do i = 1 to nr;

		sum2 = sum2 + (x[i,1] - mean)**2;

	end;

svar = sum2/(nr-1);

print svar;

tst = tinv(.975,nr-1);

print tst;

 a =tst*sqrt(svar/nr);

lower = mean - a;

upper = mean + a;

print lower upper;

 

m1 ={11 12 13,14 15 16, 17 18 19,20 21 22};

m2 = {11 12 13 14 15 16 17 18 19 20 21 22};

print m1 ;

print m2;

 

r = 4;

c= 3;

col = 3;

ind1 = col*(r-1) + c;

print ind1;

aa1= m2[1, ind1];

print aa1;

w= m1[r,c];

print w;
 
ident = i(4);
print ident;

j1 = j(1,nr);
j2 = j(nr,1);

m11 = j1*x/nr;
print m11;

diff1 = x-m11*j2;

var1 = t(diff1)*diff1/(nr-1);

print var1;