GRASS GIS 8 Programmer's Manual  8.2.2dev(2023)-3d2c704037
shapiroe.c
Go to the documentation of this file.
1 #include <stdio.h>
2 #include <math.h>
3 
4 
5 double *Cdhc_shapiro_wilk_exp(double *x, int n)
6 {
7  static double y[2];
8  double mean, b, s1, xs, sum1 = 0.0, sum2 = 0.0;
9  int i;
10 
11  for (i = 0; i < n; ++i)
12  if (i == 0 || xs > x[i])
13  xs = x[i];
14 
15  for (i = 0; i < n; ++i) {
16  sum1 += x[i];
17  sum2 += x[i] * x[i];
18  }
19 
20  s1 = sum2 - sum1 * sum1 / n;
21  mean = sum1 / n;
22  b = (mean - xs) * sqrt((double)n / (n - 1.0));
23  y[0] = b * b / s1;
24 
25 #ifdef NOISY
26  fprintf(stdout, " TEST15 SW(E) =%10.4f\n", y[0]);
27 #endif /* NOISY */
28 
29  return y;
30 }
#define x
double b
Definition: r_raster.c:39
float mean(IClass_statistics *statistics, int band)
Helper function for computing mean.
double * Cdhc_shapiro_wilk_exp(double *x, int n)
Definition: shapiroe.c:5