3. general considerations for the analysis of case-control ... - IARC

328 APPENDIX VI
SUBROUTINE MYSTCT, SUM, PROB, E,V, P, IME)
C THE EXTERNALS :QUAD : AND :B INCOF:::: ARE NECESSARY TO USE THIS
C
SUBROUTINE
INTEGER P
DOUBLE PREC IS ION T, PROB, E, V,TOT, B INCOF, PI
DIMENSION T(4)
DOUBLE PRECISION SUM,EM,EN,TEE ,ENMT,A,B,C,U,G,H,HH
b PI= 3 .I4159 26535 89793 23846 26433
u=o H=O G=O .
HH=O .
EM =T(l)+T(3)
EN = T(2) + T(4)
TEE = T(1) + T(2)
ENMT = T(3) + T(4)
C VARIABLE ::IME:: IS THE FLAT TO DETERMINE EXACT OR ASYMPTOTIC COMP
C IME = 1 GIVES EXACT
C IME = 2 GIVES ASYMPTOTIC
IF(IME.LT.2) GO TO 100
GO TO 101
C EXACT CALCULATIONS
100 SAP=SNGL(DEXP(SUM))
'[F(SAP .EQ.1 .O)GO TO 403
A=BINCOF (EM, TEE-T( 2))
B=BINCOF(EN, T(2))
C = DEXP(T(2) :: SUM)
FNUM = A:B:c
KK = ID[NT(TEE + 1.)
DO 14 J=l,KK
A=BI NCOFCEM, TEE-U)
B=B INCOF(EN, U)
C =DEXP(U::SUM)
H=A::B::C+H
G = A:B:c:u + G
HH = A:B:C:U:2 + HH
14U=U+1.
PROB= FNUM/H
E = G/H
V = HH/H - E:2
SUM=DEXP(SUM)
GO TO 677
403 A=BINCOF(TEE,T(2))
B=BINCOF(ENMT, T(4))
TOT = ENMT + TEE
c = BINCOFCTOT, EN)
PROB = (A::B)/c
EN T(2) + T(4)
E = TEE :: (EN/(TOT))
v = (EM::EN::TEE::ENMT)/( TOT::TOT::( TOT- 1 ) )
SUM=DEXP(SUM)
GO TO 677