48249       SUBROUTINE LOGISE(X,N1,Y,N2,PSTAMV,IWRITE,XIDTEM,STAT,  subroutine
48327    51   FORMAT('***** AT THE BEGINNING OF LOGISE--')  argument
48330    52   FORMAT('IBUGA3 = ',A4)
48333    53   FORMAT('N1,N2 = ',2I8)
48337    56     FORMAT('I,X(I),Y(I) = ',I8,2G15.7)
48339    55   CONTINUE
48354  1201   FORMAT('***** ERROR IN STANDARD ERROR OF BIAS CORRECTED ',
48355      1         'LOG ODDS RATIO')
48358  1203   FORMAT('      THE NUMBER OF OBSERVATIONS FOR THE FIRST ',
48359      1         'RESPONSE VARIABLES IS LESS THAN TWO')
48362  1205   FORMAT('SAMPLE SIZE = ',I8)
48374  1213   FORMAT('      THE NUMBER OF OBSERVATIONS FOR THE SECOND ',
48375      1         'RESPONSE VARIABLES IS LESS THAN TWO')
48397      1     (N12.EQ.0 .OR. N12.EQ.1) .AND.
48398      1     (N21.EQ.0 .OR. N21.EQ.1) .AND.
48399      1     (N22.EQ.0 .OR. N22.EQ.1)) GOTO1349
48407  1311     FORMAT('      ROW 1 COLUMN 1 OF THE COUNTS TABLE IS ',
48408      1           'NEGATIVE.')
48416  1321     FORMAT('      ROW 2 COLUMN 1 OF THE COUNTS TABLE IS ',
48417      1           'NEGATIVE.')
48425  1331     FORMAT('      ROW 1 COLUMN 2 OF THE COUNTS TABLE IS ',
48426      1           'NEGATIVE.')
48434  1341     FORMAT('      ROW 2 COLUMN 2 OF THE COUNTS TABLE IS ',
48435      1           'NEGATIVE.')
48454  1349 CONTINUE
48460      1            IBUGA3,IERROR)
48463      1            IBUGA3,IERROR)
48473  3000 CONTINUE
48484      1         (1.0/(AN12+0.5)) + (1.0/(AN22+0.5)))
48513  9011   FORMAT('***** AT THE END OF LOGISE--')  argument
48516  9012   FORMAT('IBUGA3,IERROR = ',A4,2X,A4)
48519  9013   FORMAT('N,N11,N12,N21,N22 = ',5I10)
48522  9014   FORMAT('P11,P12,P21,P22 = ',4G15.7)
48525  9015   FORMAT('STAT = ',G15.7)