htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/statsort.c

   1 #include <stdio.h>
   2 #include <gsl/gsl_sort.h>
   3 #include <gsl/gsl_statistics.h>
   4
   5 int
   6 main(void)
   7 {
   8   double data[5] = {17.2, 18.1, 16.5, 18.3, 12.6};
   9   double median, upperq, lowerq;
  10
  11   printf ("Original dataset:  %g, %g, %g, %g, %g\n",
  12          data[0], data[1], data[2], data[3], data[4]);
  13
  14   gsl_sort (data, 1, 5);
  15
  16   printf ("Sorted dataset: %g, %g, %g, %g, %g\n",
  17          data[0], data[1], data[2], data[3], data[4]);
  18
  19   median
  20     = gsl_stats_median_from_sorted_data (data,
  21                                          1, 5);
  22
  23   upperq
  24     = gsl_stats_quantile_from_sorted_data (data,
  25                                            1, 5,
  26                                            0.75);
  27   lowerq
  28     = gsl_stats_quantile_from_sorted_data (data,
  29                                            1, 5,
  30                                            0.25);
  31
  32   printf ("The median is %g\n", median);
  33   printf ("The upper quartile is %g\n", upperq);
  34   printf ("The lower quartile is %g\n", lowerq);
  35   return 0;
  36 }