17.01.2014 Views

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

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

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

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

APPENDIX V<br />

SUBROUTINE SYMINVCA, N, C,W, NULLTY, IFAULT, NFAC)<br />

DIMENSION A(NFAC),C(NFAC),W(N)<br />

ALGORITHM AS7 J .R. STATIST.SOC .C .(1968) VOL. 17,N0.2<br />

NROW=N<br />

I FAULT= 1<br />

IF(NROW.LE.O)GOTO loo<br />

I FAlJLT=O<br />

CALL CHOLCA, NRCW, C, NLILLTY, I FAULT, MFAC)<br />

IF(IFAULT.NE.O)GOTO loo<br />

IUI\I=(NROW::(NROW+I I)/ 2<br />

I ROW =NROW<br />

NDIAG=NN<br />

IF(C(NDIAG> .EQ .090 )GOTO 11<br />

L=ND IAG<br />

DO 10 I=IRCW,NROW<br />

W(I )=C(L)<br />

L=L+l<br />

CON T I NUE<br />

ICOL=NROW<br />

JCOL=NN<br />

MD I AG=NN<br />

L=JCOL<br />

K= NROW<br />

IF(K .EQ. 1ROW)GOTOlZ<br />

x=x-w(K)::c(L)<br />

K=K-1<br />

L=L-1<br />

I F( L .GT .MD IAG)L=L-K+I<br />

GOT0 13<br />

C(L)=X/W(IROW)<br />

IF( ICOL .EQ. IROW)GOTO 14<br />

MD IAG=MD IAG -I COL<br />

ICOL=ICOL-1<br />

JCOL=JCOL-1<br />

GOT0 15<br />

L=ND I AG<br />

DO 1 7 J= I ROW, IUROW<br />

C(L)=O .o<br />

L=L+J<br />

CONTINUE<br />

ND I AG=ND I AG - I ROW<br />

I RCWIROW-1<br />

IF(IRCW.NE.O)GOTO 16<br />

RE TLIRN<br />

END

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!