PROGRAM APM
      IMPLICIT REAL*8 (A-H,O-Z)
      PARAMETER(LG=20,LF=250)
      DIMENSION G(LG,LF), SDNUM(4,LG,LF), GVAR(LG,LG), GINV(LG,LG)
      DIMENSION NTAPER(3)
      DATA NTAPER/4,8,15/
      READ(5,*)NF,NG
      DO 20 J=1,NF
      CALL FOPEN('Moments',1,.FALSE.)
      DO 10 I=1,NG
      READ(1,*,END=1110,ERR=1120)JF,NUSE,G(J,I),SDG
   10 READ(1,*,END=1110,ERR=1120)(SDNUM(K,J,I),K=1,4)
   20 CLOSE(1)
      DO 450 I=1,NG
      WRITE(6,100)I
  100 FORMAT(/,'Moment ',I3,/,
     1  'NSE estimate ',11X,'Mean',12X,'N.S.E.',11X,'Equality chi sq',/)
      DO 450 K=1,4
      EG=0.0D0
      WTSUM=0.0D0
      CALL UM0SET(NF-1,NF-1,GVAR,LG) 
      DO 200 J=1,NF
      EG=EG+G(J,I)/(SDNUM(K,J,I)**2)
  200 WTSUM=WTSUM+1/(SDNUM(K,J,I)**2)
      DO 300 J=1,NF-2
      GVAR(J,J)=(SDNUM(K,J,I)**2) + (SDNUM(K,J+1,I)**2)
      GVAR(J,J+1)=-(SDNUM(K,J+1,I)**2)
  300 GVAR(J+1,J)=GVAR(J,J+1)
      J=NF-1
      GVAR(J,J)=SDNUM(K,J,I)**2 + SDNUM(K,J+1,I)**2
      CALL DLINDS(NF-1,GVAR,LG,GINV,LG)
      CSTAT=((G(2,I)-G(1,I))**2)*GINV(1,1)
      DO 400 J=2,NF-1
      CSTAT=CSTAT+((G(J+1,I)-G(J,I))**2)*GINV(J,J) 
  400 CSTAT=CSTAT+2*(G(J+1,I)-G(J,I))*(G(J,I)-G(J-1,I))*GINV(J-1,J)
      DF=(NF-1)*1.0
      P=CUMCHI(CSTAT,DF)
      EG=EG/WTSUM
      ONSE=1/DSQRT(WTSUM)
      IF(K.EQ.1)WRITE(6,500)EG,ONSE,1-P
  450 IF(K.NE.1)WRITE(6,600)NTAPER(K-1),EG,ONSE,1-P
  500 FORMAT(' i.i.d.',7X,1PE16.7,3X,1PE16.7,5X,1PE16.7)
  600 FORMAT(I2,'% taper',5X,1PE16.7,3X,1PE16.7,5X,1PE16.7)
      STOP
 1110 WRITE(6,1115)J,I
 1115 FORMAT(/,' *** Error in file', I3, ' function ',I3,/)
      CALL TERM
 1120 WRITE(6,1125)J,I
 1125 FORMAT(/,' *** Error in file', I3, ' function ',I3,/)
      CALL TERM
      END