* G02FAF Example Program Text * Mark 14 Release. NAG Copyright 1989. * .. Parameters .. INTEGER NMAX, LDSRES PARAMETER (NMAX=24,LDSRES=NMAX) INTEGER NIN, NOUT PARAMETER (NIN=5,NOUT=6) * .. Local Scalars .. DOUBLE PRECISION RMS INTEGER I, IFAIL, IP, J, N, NRES * .. Local Arrays .. DOUBLE PRECISION H(NMAX), RES(NMAX), SRES(LDSRES,4) * .. External Subroutines .. EXTERNAL G02FAF * .. Executable Statements .. WRITE (NOUT,*) 'G02FAF Example Program Results' * Skip heading in data file READ (NIN,*) READ (NIN,*) N, IP, NRES, RMS IF (NRES.LT.NMAX) THEN DO 20 I = 1, N READ (NIN,*) RES(I), H(I) 20 CONTINUE IFAIL = 1 * CALL G02FAF(N,IP,NRES,RES,H,RMS,SRES,LDSRES,IFAIL) * IF (IFAIL.EQ.0) THEN WRITE (NOUT,*) WRITE (NOUT,*) ' Internally Internally' WRITE (NOUT,*) + 'Obs. standardized standardized Cook''s D Atkinson''s T' WRITE (NOUT,*) ' residuals residuals' WRITE (NOUT,*) DO 40 I = 1, NRES WRITE (NOUT,99999) I, (SRES(I,J),J=1,4) 40 CONTINUE ELSE WRITE (NOUT,*) WRITE (NOUT,99998) ' ** G02FAF returned with IFAIL = ', + IFAIL END IF END IF * 99999 FORMAT (1X,I2,4F13.3) 99998 FORMAT (1X,A,I5) END