/* Geary statistic */
new; cls;
x=rndn(100,1);
call G1(x);

print;
call G2(x);



proc G1(x);
   local n,stat;
   if ismiss(x);
      errorlog "Warning: missing data found.";
      x=packr(x);
   endif;
   n=rows(x);
   stat=sumc(abs(x-meanc(x)))/sqrt(n*sumc((x-meanc(x))^2));
/* results */
   print "Geary Test:";
   print/rz "n          =" n;
   print    "stat ( G1 )=" stat;
   retp(stat);
endp;

proc G2(x);
   local n,stat;
   if ismiss(x);
      errorlog "Warning: missing data found.";
      x=packr(x);
   endif;
   n=rows(x);
   stat=sumc(abs(x-median(x)))/sqrt(n*sumc((x-meanc(x))^2));
/* results */
   print "Geary Test(median version):";
   print/rz "n          =" n;
   print    "stat ( G2 )=" stat;
   retp(stat);
endp;