For earlier releases of SAS/IML, you can use the following module, which appears in my book Statistical Programming with SAS/IML Software: /* Var: return sample variance of each column of a data matrix */ /* x is a matrix that contains the data */ start Var(x); mean = x[:,]; countn = j(1, ncol(x)); /* allocate vector for counts */ do i = 1 to ncol(x); /* count nonmissing values */ countn = sum(x[,i]^=.);/* in each column */ end; var = (x-mean)[##,] / (countn-1); return ( var ); finish; For STD, use sqrt(var(x)).
... View more