1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
#include <stdio.h>
#include <math.h>
double *Cdhc_omnibus_moments(double *x, int n)
{
double diff, mean = 0., fssm, tssm, sum_cube = 0., sum_four = 0.,
sum_sq = 0.;
static double y[2];
int i;
for (i = 0; i < n; ++i)
mean += x[i];
mean /= n;
for (i = 0; i < n; ++i) {
diff = x[i] - mean;
sum_sq += diff * diff;
sum_cube += diff * diff * diff;
sum_four += diff * diff * diff * diff;
}
/*
fprintf (stdout,"n %d x-bar %g sum^2 %g sum^3 %g sum^4 %g
\n",n,mean,sum_sq,sum_cube,sum_four);
*/
tssm = sqrt((double)n) * sum_cube / pow(sum_sq, 1.5);
fssm = n * sum_four / (sum_sq * sum_sq);
#ifdef NOISY
fprintf(stdout, " TESTS OF COMPOSITE DISTRIBUTIONAL HYPOTHESES\n");
fprintf(stdout, " TEST1 TSM =%10.4f FSM =%10.4f\n", tssm, fssm);
#endif /* NOISY */
y[0] = tssm;
y[1] = fssm;
return y;
}
|
