SUBROUTINE LRTEST(TITLE,S,SY,R,RY,N,YR)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      PARAMETER (IDIM=50)
      INTEGER YR(IDIM),SUMDF
      DIMENSION S(IDIM),SY(IDIM),R(IDIM),RY(IDIM)
      CHARACTER*128 TITLE
      WRITE(7,10)CHAR(12),TITLE
   10 FORMAT(A1,A128/' Test of the hypothesis that adult and young ',
     ,'survival and capture probabilities are the same'/
     /' (from Pollock 1981; Biometrics 37:521-529)'/)
      SUMCHI=0.
      SUMDF=0
      DO 90 I=1,N
        WRITE(7,20)YR(I)
   20   FORMAT(/43X,'Adults     Young'/' i=',I4,35X,'Released   1st'
     /     ,' Cap'/41X,'   in i       in i')
        WRITE(7,30)
   30   FORMAT(40X,'+----------+----------+')
        WRITE(7,40)R(I),RY(I)
   40   FORMAT(' Released, recaptured after i',11X,2('!',F10.2),'!')
        R1=R(I)+RY(I)
        TOTL=DMAX1(S(I)+SY(I),.1D1)
        R2=TOTL-R1
        EX1=S(I)*R1/TOTL
        EX2=S(I)*R2/TOTL
        EX3=SY(I)*R1/TOTL
        EX4=SY(I)*R2/TOTL
        WRITE(7,50)EX1,EX3
   50   FORMAT(' Expected value             ',12X,2('!',F10.2),'!')
        WRITE(7,30)
        WRITE(7,60)S(I)-R(I),SY(I)-RY(I)
   60   FORMAT(' Released, not recaptured after i',7X,2('!',F10.2),'!')
        WRITE(7,50)EX2,EX4
        WRITE(7,30)
        CHI=0.
        DF=0.
        CPRB=.0D0
        IF(EX1.LT..2D1.OR.EX2.LT..2D1.OR.EX3.LT..2D1.OR.
     .     EX4.LT..2D1)GO TO 70
        CHI=(R(I)-EX1)**2/EX1+(S(I)-R(I)-EX2)**2/EX2+
     +      (RY(I)-EX3)**2/EX3+(SY(I)-RY(I)-EX4)**2/EX4
        DF=1.
        SUMCHI=SUMCHI+CHI
        SUMDF=SUMDF+1
        IF(CHI.GT..0D0)CPRB=CHIPRB(1,CHI)
   70   WRITE(7,80)CHI,INT(DF),CPRB
   80   FORMAT(' Chi-square = ',F10.4,' with ',I1,' degrees of freedom',
     ,   ', Probability = ',F10.4)
   90   CONTINUE
      CPRB=.0D0
      IF(SUMCHI.GT..0D0)CPRB=CHIPRB(SUMDF,SUMCHI)
      WRITE(7,100)SUMCHI,FLOAT(SUMDF),CPRB
  100 FORMAT(/20X,' Total chi-square  =',F12.2
     /       /20X,' Degrees of freedom=',F12.0
     /       /20X,' Overall Prob.     =',F12.4/)
      RETURN
      END
      FUNCTION CHIPRB(IDF,CHSQ)
C   FUNCTION COMPUTES THE PROB. OF THE CHI SQUARE VALUE 'CHSQ' WITH
C   'DF' Degrees of freedom.
C
C   SUBROUTINES NEEDED:
C        FUNCTION ALGAMA
C
      DOUBLE PRECISION CHIPRB,CHSQ
      CHIPRB=1.
      IF(CHSQ .LE. 0.002 .OR. IDF .LT. 1) GO TO 50
      A=FLOAT(IDF)/2.
      X=.5*CHSQ
      TERM=1.
      SUM=0.
      COFN=A
      IF(X.GE.A.AND.X.GE.13.)GO TO 20
C              CONVERGENT SERIES FOR X .LT. A OR .LT. 13.
      FACT=-A
   10 TEMP=SUM
      SUM=SUM+TERM
      COFN=COFN+1.
      TERM=TERM*X/COFN
      IF(SUM-TEMP)40,40,10
C                     ASYMPTOTIC SERIES FOR X .GTE. A AND X .GTE. 13.
   20 CHIPRB=0.
      FACT=X
   30 TEMP=SUM
      SUM=SUM+TERM
      COFN=COFN-1.
      RATIO=COFN/X
      TERM=TERM*RATIO
      IF(SUM.GT.TEMP)GO TO 30
C
   40 ARGMNT=ALOG(SUM)-X+A*ALOG(X)-ALGAMA(A)
      IF(ARGMNT .LT. -175.)GO TO 50
      CHIPRB=CHIPRB+EXP(ARGMNT)/FACT
   50 RETURN
      END
C$ SEGMENT 'WRMSEG'
      FUNCTION ALGAMA(A)
C   FUNCTION COMPUTES THE LOG OF THE GAMMA FUNCTION GIVEN
C   THAT THE ARGUMENT IS AN INTEGER OR A MULTIPLE OF 1/2.
C
      ALGAMA=0.0
      I=INT(A)
      B=FLOAT(I)
      IF(B.NE.A)GO TO 20
      IF(A.LE. 2.0)GO TO 40
      DO 10 J=2,I-1
   10      ALGAMA=ALGAMA+ALOG(FLOAT(J))
      GO TO 40
C                  A IS A MULTIPLE OF 1/2
   20 ALGAMA=0.5723649429
      IF(A.LE. 0.5)GO TO 40
      DO 30 J=1,I*2-1,2
   30      ALGAMA=ALGAMA+ALOG(FLOAT(J)/2.)
   40 RETURN
      END