/* skewness */
new; cls;
x=rndn(10000,1);
print skew(x);





proc skew(x);
   local n,sum,i,xbar,s,z,sk;
   n=rows(x);
   sum=0;
   i=1;
   do while i<=n;
      sum=sum+x[i,.];
      i=i+1;
   endo;
   xbar=sum/n;
   s=sqrt(sumc((x-xbar)^2)/(n-1));
   z=(x-xbar)./s';
   sk=sumc((z)^3)/n;
   retp(sk);
endp;