/* kurtosis */
new; cls;
x=rndn(10000,1);
print kurt(x);





proc kurt(x);
   local n,sum,i,xbar,s,z,kur;
   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';
   kur=sumc((z)^4)/n;
   retp(kur);
endp;