SUBROUTINE TAURE(NT,NU,ALPH,CRTAU) TAURE001 C***********************************************************************TAURE002 C* TAURE003 C* TAURE COMPUTES THE CRITICAL VALUE FOR REJECTION OF STANDARDIZED TAURE004 C* RESIDUALS WITH CONTROL OF TYPE I ERROR. TAURE005 C* TAURE006 C* TAURE007 C* INPUT: TAURE008 C* NT- NUMBER OF OBSERVATIONS TAURE009 C* NU- DEGREES OF FREEDOM TAURE010 C* ALPH- DESIRED PROBABILITY OF TYPE I ERROR TAURE011 C* TAURE012 C* OUTPUT: TAURE013 C* CRTAU- CRITICAL VALUE (TAU-MAX) TAURE014 C* TAURE015 C* TAURE016 C* REFERENCE: TAURE017 C* A.J. POPE (1976)- THE STATISTTAURE018 C* OF RESIDUALS AND THE DETECTIOTAURE019 C* OF OUTLIERS; U.S. DEPT OR COMTAURE020 C* NOAA TECHNICAL REPORT TAURE021 C* NO. 65 NGS1. TAURE022 C* TAURE023 C***********************************************************************TAURE024 IMPLICIT REAL*8(A-H,O-Z) TAURE025 DATA PI/ 3.1415926535898 / TAURE026 PD = 2. /PI TAURE027 S = 1. TAURE028 WNU = NU TAURE029 U = WNU -1. TAURE030 IF( U.EQ.0. ) GO TO 72 TAURE031 IF ( ALPH.EQ.0. ) GO TO 72 TAURE032 IF ( ALPH.LT.1. ) GO TO 10 TAURE033 CRTAU = 0. TAURE034 C TAURE035 RETURN TAURE036 C TAURE037 10 Q = NT TAURE038 IF ( ALPH.GT.0.5 ) GO TO 19 TAURE039 AA = ALPH / Q TAURE040 DELT = AA TAURE041 DO 18 I = 1,100 TAURE042 XI = I TAURE043 DELT = DELT * ALPH * (( XI*Q - 1.)/(( XI+1.)*Q)) TAURE044 IF ( DELT.LT.1.D-14 ) GO TO 20 TAURE045 18 AA = AA + DELT TAURE046 19 AA = 1. - (1.-ALPH)**(1./Q) TAURE047 20 P = 1. - AA TAURE048 IF(U.EQ.1. ) GO TO 71 TAURE049 F = 1.3862943611199 - 2.*DLOG(AA) TAURE050 G = DSQRT(F) TAURE051 X = G - (2.515517 + 0.802853*G + 0.010328*F) TAURE052 $ / (1. + 1.432788*G + F*(0.189269 + 0.001308*G)) TAURE053 Y = X*X TAURE054 A = X*(1. + Y)/4. TAURE055 B = X*(3. + Y*(16. + 5.*Y))/96. TAURE056 C = X*(-15. + Y*(17. + Y*(19. + 3.*Y)))/384. TAURE057 E = X*(-945. + Y*(-1920. + Y*(1482. + Y*(776. + 79.*Y))))/92160. TAURE058 V = 1./U TAURE059 T = X + V*(A + V*(B + V*(C + E*V))) TAURE060 S = T/DSQRT(U + T*T) TAURE061 UM = U - 1. TAURE062 DELL = 1. TAURE063 DO 70 M = 1,50 TAURE064 H = 1. - S*S TAURE065 R = 0.0 TAURE066 IF ( DMOD(U,2.D0).EQ.0.0 ) GO TO 49 TAURE067 DD = DSQRT(H) TAURE068 N = 0.5*UM TAURE069 DO 45 I = 1,N TAURE070 Z = 2*I TAURE071 R = R + DD TAURE072 D = DD TAURE073 45 DD = DD * H * (Z/(Z+1.)) TAURE074 R = PD*(R*S + DARSIN(S)) TAURE075 D = PD*D*UM TAURE076 GO TO 61 TAURE077 49 DD = 1. TAURE078 N = 0.5*U TAURE079 DO 55 I = 1,N TAURE080 Z = 2*I TAURE081 R = R + DD TAURE082 D = DD TAURE083 55 DD = DD*H*((Z-1.)/Z) TAURE084 R = R*S TAURE085 D = D*UM TAURE086 61 CONTINUE TAURE087 DEL = (P-R)/D TAURE088 IF( DABS( DEL/DELL ) .GT.0.5) GO TO 72 TAURE089 DELL = DEL TAURE090 S = S + DEL TAURE091 IF( DABS(DEL) .LT. 1.D-8 ) GO TO 72 TAURE092 70 CONTINUE TAURE093 GO TO 72 TAURE094 71 S =DSIN(P/PD) TAURE095 72 CRTAU = S*DSQRT(WNU) TAURE096 RETURN TAURE097 END TAURE098