/* nag_summary_stats_1var (g01aac) Example Program.
 *
 * Copyright 1990 Numerical Algorithms Group.
 *
 * Mark 1, 1990.
 *
 * Mark 5 revised, 1998.
 * Mark 8 revised, 2004.
 *
 */

#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nagg01.h>

int main(void)
{
  Integer exit_status=0, i, j, n, nprob, nvalid, weight;
  NagError fail;
  double wsum, *wt=0, *x=0, xkurt, xmax, xmean, xmin, xsd, xskew;

  INIT_FAIL(fail);
  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vprintf("nag_summary_stats_1var (g01aac) Example Program Results\n");
  Vscanf("%ld", &nprob);
  for (j=1; j<=nprob; j++)
    {
      Vscanf("%ld %ld", &n, &weight);
      Vprintf("Problem  %5ld\n", j);
      Vprintf("Number of cases %ld\n", n);
      if (n>=1)
        {
          if ( !( wt = NAG_ALLOC(n, double)) ||
               !( x = NAG_ALLOC(n, double)) )
            {
              Vprintf("Allocation failure\n");
              exit_status = -1;
              goto END;
            }
        }
      else
        {
          Vprintf("Invalid n.\n");
          exit_status = 1;
          return exit_status;
        }
      for (i=0; i<n; i++)
        Vscanf("%lf", &x[i]);
      Vprintf("Data as input -\n");
      for (i=0; i<n; i++)
        Vprintf("%12.1f%c", x[i], (i%5==4 || i==n-1) ? '\n' : ' ');
      if (weight)
        {
          Vprintf("Weights as input -\n");
          for (i=0; i<n; i++)
            Vscanf("%lf", &wt[i]);
          for (i=0; i<n; i++)
            Vprintf("%12.1f%c", wt[i], (i%5==4 || i==n-1) ? '\n' : ' ');
          /* nag_summary_stats_1var (g01aac).
           * Mean, variance, skewness, kurtosis, etc., one variable,
           * from raw data
           */
          nag_summary_stats_1var(n, x, wt, &nvalid, &xmean, &xsd, &xskew,
                                 &xkurt, &xmin, &xmax, &wsum, &fail);
        }
      else
        /* nag_summary_stats_1var (g01aac), see above. */
        nag_summary_stats_1var(n, x, (double *)0, &nvalid, &xmean, &xsd,
                               &xskew, &xkurt, &xmin, &xmax, &wsum, &fail);

      if (fail.code==NE_NOERROR)
        {
          Vprintf("\n");
          Vprintf("Successful call of g01aac\n");
          Vprintf("No. of valid cases %5ld\n", nvalid);
          Vprintf("Mean           %13.1f\n", xmean);
          Vprintf("Std devn       %13.1f\n", xsd);
          Vprintf("Skewness       %13.1f\n", xskew);
          Vprintf("Kurtosis       %13.1f\n", xkurt);
          Vprintf("Minimum        %13.1f\n", xmin);
          Vprintf("Maximum        %13.1f\n", xmax);
          Vprintf("Sum of weights %13.1f\n", wsum);
        }
      else
        {
          Vprintf("Unsuccessful call of g01aac\n");
          Vprintf("%s \n", fail.message);
          if (fail.code==NE_CASES_ONE)
            {
              Vprintf("No. of valid cases %5ld\n", nvalid);
              Vprintf("Mean           %13.1f\n", xmean);
              Vprintf("Minimum        %13.1f\n", xmin);
              Vprintf("Maximum        %13.1f\n", xmax);
              Vprintf("Sum of weights %13.1f\n", wsum);
              Vprintf("Std devn and coeffts of skewness\n");
              Vprintf("and kurtosis not defined\n");
              exit_status = 2;
            }
          else
            {
              Vprintf("Error from nag_summary_stats_1var (g01aac).\n%s\n",
                      fail.message);
              exit_status = 1;
              goto END;
            }
        }
    }
 END:
  if (wt) NAG_FREE(wt);
  if (x) NAG_FREE(x);
  return exit_status;
}