|
Packit |
67cb25 |
#include <stdio.h>
|
|
Packit |
67cb25 |
#include <gsl/gsl_rstat.h>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
int
|
|
Packit |
67cb25 |
main(void)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double data[5] = {17.2, 18.1, 16.5, 18.3, 12.6};
|
|
Packit |
67cb25 |
double mean, variance, largest, smallest, sd,
|
|
Packit |
67cb25 |
rms, sd_mean, median, skew, kurtosis;
|
|
Packit |
67cb25 |
gsl_rstat_workspace *rstat_p = gsl_rstat_alloc();
|
|
Packit |
67cb25 |
size_t i, n;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* add data to rstat accumulator */
|
|
Packit |
67cb25 |
for (i = 0; i < 5; ++i)
|
|
Packit |
67cb25 |
gsl_rstat_add(data[i], rstat_p);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
mean = gsl_rstat_mean(rstat_p);
|
|
Packit |
67cb25 |
variance = gsl_rstat_variance(rstat_p);
|
|
Packit |
67cb25 |
largest = gsl_rstat_max(rstat_p);
|
|
Packit |
67cb25 |
smallest = gsl_rstat_min(rstat_p);
|
|
Packit |
67cb25 |
median = gsl_rstat_median(rstat_p);
|
|
Packit |
67cb25 |
sd = gsl_rstat_sd(rstat_p);
|
|
Packit |
67cb25 |
sd_mean = gsl_rstat_sd_mean(rstat_p);
|
|
Packit |
67cb25 |
skew = gsl_rstat_skew(rstat_p);
|
|
Packit |
67cb25 |
rms = gsl_rstat_rms(rstat_p);
|
|
Packit |
67cb25 |
kurtosis = gsl_rstat_kurtosis(rstat_p);
|
|
Packit |
67cb25 |
n = gsl_rstat_n(rstat_p);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
printf ("The dataset is %g, %g, %g, %g, %g\n",
|
|
Packit |
67cb25 |
data[0], data[1], data[2], data[3], data[4]);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
printf ("The sample mean is %g\n", mean);
|
|
Packit |
67cb25 |
printf ("The estimated variance is %g\n", variance);
|
|
Packit |
67cb25 |
printf ("The largest value is %g\n", largest);
|
|
Packit |
67cb25 |
printf ("The smallest value is %g\n", smallest);
|
|
Packit |
67cb25 |
printf( "The median is %g\n", median);
|
|
Packit |
67cb25 |
printf( "The standard deviation is %g\n", sd);
|
|
Packit |
67cb25 |
printf( "The root mean square is %g\n", rms);
|
|
Packit |
67cb25 |
printf( "The standard devation of the mean is %g\n", sd_mean);
|
|
Packit |
67cb25 |
printf( "The skew is %g\n", skew);
|
|
Packit |
67cb25 |
printf( "The kurtosis %g\n", kurtosis);
|
|
Packit |
67cb25 |
printf( "There are %zu items in the accumulator\n", n);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_rstat_reset(rstat_p);
|
|
Packit |
67cb25 |
n = gsl_rstat_n(rstat_p);
|
|
Packit |
67cb25 |
printf( "There are %zu items in the accumulator\n", n);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_rstat_free(rstat_p);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return 0;
|
|
Packit |
67cb25 |
}
|