*     G02CHF Example Program Text
*     Mark 14 Revised. NAG Copyright 1989.
*     .. Parameters ..
      INTEGER          K1, N, K, LDSSPZ, LDRZ, LDCOEF, LDRZNV, LDCZ,
     +                 LDWKZ
      PARAMETER        (K1=3,N=5,K=K1-1,LDSSPZ=K1,LDRZ=K1,LDCOEF=K,
     +                 LDRZNV=K,LDCZ=K,LDWKZ=K)
      INTEGER          NIN, NOUT
      PARAMETER        (NIN=5,NOUT=6)
*     .. Local Scalars ..
      INTEGER          I, IFAIL, J
*     .. Local Arrays ..
      DOUBLE PRECISION COEF(LDCOEF,3), CZ(LDCZ,K), RESULT(13),
     +                 RZ(LDRZ,K1), RZNV(LDRZNV,K), SSPZ(LDSSPZ,K1),
     +                 WKZ(LDWKZ,K)
*     .. External Subroutines ..
      EXTERNAL         G02CHF
*     .. Executable Statements ..
      WRITE (NOUT,*) 'G02CHF Example Program Results'
*     Skip heading in data file
      READ (NIN,*)
      READ (NIN,*) ((SSPZ(I,J),J=1,K1),I=1,K1),
     +  ((RZ(I,J),J=1,K1),I=1,K1)
      WRITE (NOUT,*)
      WRITE (NOUT,*) 'Sums of squares and cross-products about zero:'
      WRITE (NOUT,99999) (J,J=1,K1)
      WRITE (NOUT,99998) (I,(SSPZ(I,J),J=1,K1),I=1,K1)
      WRITE (NOUT,*)
      WRITE (NOUT,*) 'Correlation-like coefficients:'
      WRITE (NOUT,99999) (J,J=1,K1)
      WRITE (NOUT,99998) (I,(RZ(I,J),J=1,K1),I=1,K1)
      WRITE (NOUT,*)
      IFAIL = 1
*
      CALL G02CHF(N,K1,K,SSPZ,LDSSPZ,RZ,LDRZ,RESULT,COEF,LDCOEF,RZNV,
     +            LDRZNV,CZ,LDCZ,WKZ,LDWKZ,IFAIL)
*
      IF (IFAIL.NE.0) THEN
         WRITE (NOUT,*)
         WRITE (NOUT,99997) ' ** G02CHF returned with IFAIL = ', IFAIL
      ELSE
         WRITE (NOUT,*) 'Vble     Coef       Std err      t-value'
         WRITE (NOUT,99996) (I,(COEF(I,J),J=1,3),I=1,K)
         WRITE (NOUT,*)
         WRITE (NOUT,*) 'Analysis of regression table :-'
         WRITE (NOUT,*)
         WRITE (NOUT,*)
     +'      Source        Sum of squares  D.F.    Mean square     F-val
     +ue'
         WRITE (NOUT,*)
         WRITE (NOUT,99995) 'Due to regression', (RESULT(I),I=1,4)
         WRITE (NOUT,99995) 'About  regression', (RESULT(I),I=5,7)
         WRITE (NOUT,99995) 'Total            ', (RESULT(I),I=8,9)
         WRITE (NOUT,*)
         WRITE (NOUT,99994) 'Standard error of estimate =', RESULT(10)
         WRITE (NOUT,99994) 'Multiple correlation (R)   =', RESULT(11)
         WRITE (NOUT,99994) 'Determination (R squared)  =', RESULT(12)
         WRITE (NOUT,99994) 'Corrected R squared        =', RESULT(13)
         WRITE (NOUT,*)
         WRITE (NOUT,*)
     +     'Inverse of correlation matrix of independent variables:'
         WRITE (NOUT,99993) (J,J=1,K)
         WRITE (NOUT,99992) (I,(RZNV(I,J),J=1,K),I=1,K)
         WRITE (NOUT,*)
         WRITE (NOUT,*) 'Modified inverse matrix:'
         WRITE (NOUT,99993) (J,J=1,K)
         WRITE (NOUT,99992) (I,(CZ(I,J),J=1,K),I=1,K)
      END IF
*
99999 FORMAT (1X,3I10)
99998 FORMAT (1X,I4,3F10.4)
99997 FORMAT (1X,A,I5)
99996 FORMAT (1X,I3,3F12.4)
99995 FORMAT (1X,A,F14.4,F8.0,2F14.4)
99994 FORMAT (1X,A,F8.4)
99993 FORMAT (1X,2I10)
99992 FORMAT (1X,I4,2F10.4)
      END