#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#define LF 250
#define LG 20


FILE *fp1;
extern double sq(double a);
main()
{
  char hfile[40], s[]=" ", c[]="%";
/* INT */
  int nf, ng, i, j, k, jf, nuse, ss;
  int ntaper[3]={4, 8, 15}, which=1, stat;
/* DOUBLE */
  double g[LG][LF], sdnum[4][LG][LF], gvar[LG][LG], inv[LG][LG];
  double sdg, wtsum, eg, nse, cstat, p, q, bound, df;

  ss=scanf("%d%d", &nf, &ng);
  if (ss!=2) { printf("read error in control file, line 1\n"); return -1;}
  for (j=0; j< nf; j++) { 
     ss=scanf("%s", &hfile);
     if (ss!=1) {printf("error reading file name\n"); return -1;}
     fp1=fopen(hfile,"r");
     if (fp1==NULL) {printf("failure to open file %s for input\n", hfile);
                     return -1;}
     for (i=0; i<ng; i++) {
        ss=0;
        ss=ss+fscanf(fp1,"%d %d", &jf, &nuse);
        ss=ss+fscanf(fp1,"%le",&g[j][i]);
        ss=ss+fscanf(fp1,"%le",&sdg);
        for (k=0; k<4; k++) 
           ss=ss+fscanf(fp1,"%le",&sdnum[k][j][i]);  /*j th file i th function*/
        if (ss!=8) {printf("error in file %s function %3d. \n",hfile,i); 
                    return -1;}

     } /* for i<=ng */
     fclose(fp1);
  } /* for j<= nf */
  for (i=0; i<ng; i++) { 
     printf("\nMoment %3d\n", i+1);
     printf("NSE estimate %10s Mean %10s N.S.E. %9s Equality chi sq\n",s,s,s);
     for (k=0; k<4; k++) { 
        eg=0; nse=0; wtsum=0; 
        umoset(nf-1,gvar,LG);
        for  (j=0; j<nf; j++) {
           eg=eg+g[j][i]/sq(sdnum[k][j][i]);
           wtsum=wtsum+1/sq(sdnum[k][j][i]); 
        }
        for  (j=0; j<nf-1; j++) {
           gvar[j][j]=sq(sdnum[k][j][i])+sq(sdnum[k][j+1][i]);
           gvar[j][j+1]=-sq(sdnum[k][j+1][i]);
           gvar[j+1][j]=gvar[j][j+1];
        }
        dlinds(nf-1,gvar,LG,inv,LG);
        cstat=sq(g[1][i]-g[0][i])*inv[0][0]; 
        for (j=1; j<nf-1; j++) {
           cstat=cstat+sq(g[j+1][i]-g[j][i])*inv[j][j] + 
                 2*(g[j+1][i]-g[j][i])*(g[j][i]-g[j-1][i])*inv[j-1][j]; 
        }  
        df=nf-1;
        cdfchi(&which,&p,&q,&cstat,&df,&stat,&bound); 
        eg=eg/wtsum;
        nse=1/sqrt(wtsum);
        if (k==0)
          printf(" i.i.d. %6s",s);
        else
          printf("%2d%s taper%4s ",ntaper[k-1],c,s);
        printf("%3s%16.7e %16.7e %3s %16.7e\n",s,eg,nse,s,q);
     }
  }
return 1;  
} /* main*/