/* nag_robust_trimmed_1var(g07ddc) Example Program.
 *
 * Copyright 1996 Numerical Algorithms Group.
 *
 * Mark 4, 1996.
 *
 */

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

#define NMAX 1000
int main(void)
{

  /* Local variables */
  double tvar, wvar;
  double alpha, x[1000], tmean, wmean, propn, sx[1000];

  Integer i, k;
  Integer n;

  Vprintf("g07ddc Example Program Results\n\n");
  /* Skip heading in data file */
  Vscanf("%*[^\n] ");

  Vscanf("%ld ", &n);
  for (i = 1; i <= n; ++i)
    Vscanf("%lf ", &x[i - 1]);
  Vscanf("%lf ", &alpha);

  g07ddc(n, x, alpha, &tmean, &wmean, &tvar, &wvar, &k, sx, NAGERR_DEFAULT);

  propn = (double) k / n;
  propn = 100.0 - propn * 200.0;
  Vprintf("Statistics from middle %6.2f%% of data\n\n", propn);
  Vprintf("               Trimmed-mean = %11.4f\n", tmean);
  Vprintf("   Variance of Trimmed-mean = %11.4f\n\n", tvar);
  Vprintf("            Winsorized-mean = %11.4f\n", wmean);
  Vprintf("Variance of Winsorized-mean = %11.4f\n", wvar);
  return EXIT_SUCCESS;
}