geary.c

Go to the documentation of this file.

#include <stdio.h>
#include <math.h>
 
double *Cdhc_geary_test(double *x, int n)
{
    int i;
    static double y[2];
    double diff, s = 0.0, mean = 0.0;
 
    y[0] = 0.0;
    for (i = 0; i < n; ++i)
        mean += x[i];
 
    mean /= n;
 
    for (i = 0; i < n; ++i) {
        diff = x[i] - mean;
        y[0] += fabs(diff);
        s += diff * diff;
    }
 
    s *= n;
    y[0] /= sqrt(s);
    y[1] = (y[0] - 0.7979) * sqrt((double)n) / 0.2123;
 
#ifdef NOISY
    fprintf(stdout, "  TEST2  GTN    =%10.4f   Z(GTN) =%10.4f\n", y[0], y[1]);
#endif /* NOISY */
 
    return y;
}