/* nag_regsn_std_resid_influence(g02fac) Example Program
 *
 * Copyright 1990 Numerical Algorithms Group.
 *
 * Mark 1, 1990.
 */

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

#define NMAX 24

int main(void)
{
  double rms;
  Integer i, ip, j, n, nres;
  double h[NMAX], res[NMAX], sres[NMAX][4];

  Vprintf("g02fac Example Program Results\n");
  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vscanf("%ld %ld %ld %lf", &n, &ip, &nres, &rms);
  if (nres<NMAX)
    {
      for (i=0; i<n; i++)
        Vscanf("%lf%lf",&res[i],&h[i]);
      g02fac(n, ip, nres, res, h, rms, &sres[0][0], NAGERR_DEFAULT);
      Vprintf("\n");
      Vprintf("         Internally     Externally\n");
      Vprintf(" Obs.    studentized    studentized    Cook's D    Atkinson's T"
              "\n");
      Vprintf("          residuals      residuals\n\n");
      for (i=0; i<nres; i++)
        {
          Vprintf("%2ld", i+1);
          for (j=0; j<4; j++)
            Vprintf("%14.3f", sres[i][j]);
          Vprintf("\n");
        }
      return EXIT_SUCCESS;
    }
  else
    {
      Vfprintf(stderr, "nres is out of range: nres = %5ld\n", nres);
      return EXIT_FAILURE;
    }
}