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.

302 APPENDIX IV<br />

C UPON CONVERGENCE OF MAXIMUM ITERATIONS WRITE OUT RESULTS<br />

GO TO 1<br />

9 wRITEC~, 108)(B(J>, J=l ,NP)<br />

108 FORMAT(' ESTIMATED PARNETER VECTOR',(/lX,lOF12.6))<br />

WRITE(6,109>(SCORE(J), J=l ,NP)<br />

109 FORMAT(' FIRST DERIVATIVE LOG-LIKELIHOOD', (/1X, 10F12.6))<br />

WRITE(6,llO)<br />

110 FOR~T(' INFORMATION MATRIX')<br />

DO 10 J=l,NP<br />

K=J::( J-1 )/2+1<br />

JJ=J::(J+l)/?<br />

10 WRITE(^, 111) (COVI(I), I=K, JJ)<br />

111 FORMAT(lX, 10F12.6)<br />

C INVERT INFORMATION M4TRIX WID WRITE OUT<br />

C CALL SYM'[NV(COV I, NP, COV, W, NULLTY, I FAULT, NM1)<br />

DO 600 I=l,NMl<br />

COv(I)=CovI(I)<br />

600 CONTINUE<br />

CALL SINV(COV,IVP, EPS, IER)<br />

IF (IER.NE.0) WRITE (6,501 )IER<br />

WRITE(6,112)<br />

112 FORM4T(' ESTIMATED COVARIANCE MATRIX')<br />

DO 11 J=l,NP<br />

K=J::( J- 1 )/ 2+ 1<br />

JJ=J::(J+I )/2<br />

SCORE(J~=B(J~/SQRT(COV(JJ))<br />

11 WRITE(6,ll l)(COV(I), I=K, JJ)<br />

WRITE(~,~~~)(SCORE(J), J=l ,NP)<br />

113 FORMAT(' STANDARDIZED REGRESSION COEFFICIENTS',(/lX,1OOF12.6))<br />

RETURN<br />

END

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

Saved successfully!

Ooh no, something went wrong!