108.g
APPENDIX C
,
FORTRAN IV Program for Calculation of Statistical Pafameters Necessary for— Pfedicting the
Stability of Suspensions of Viruses
,
_
C
CALCULATIONS FOR PREDICTING TITER. DR. DONALD GREIFF (ROBINSON)
C
DEPT OF PATHOLOGY (COMPUTING CENTER)
'
'
,
DIMENSION FK (50), SD (50), SE (50), TP (50), ALPH (13)
20 READ 900, N
,
'
.
900 FORMAT (12)
,
,
.
.
PRINT 902
‘
,
902 FORMAT (1 HO, 10 X, 31 H STATISTICS FOR PREDICTING TITER)
IC = 1
-
,
_
4 READ 901, TP (IC), M, (ALPH (I), I = 1, 13)
_
901 FORMAT (F 10.0, 8X, 12, 13 A4)
PRINT 903, (ALPH (I), I = 1, 13), TP (IC)
_
‘
…
'
903 FORMAT (1 H1, 13 A 4, 2X, 14 H TEMPE.
,
RATURE IS, F 7.2/14 HO INPUT DATA IS)
904 FORMAT (F 10.0, F 10.0)
ST
=0.
PRINT9OS,T,Y
_…
SST = 0.
905 FORMAT (1 x, 2F 10.4)
SY
=
O.
-
ST aj=.-ST-+T
,!
'
SSY = o.
SST = SST +Œ*Tv »
,
STY=O.
'
DO1
I = 1, M
SSY =_SSY+,-,Y%Y.
READ 904, T, Y
1 STY
,
"
C
COMPUTE KI, SD, SE, OUTPUT
'
'
f [
:
’
FK (IC) = (STY—(ST*SY) / FM) / (SST—(ST*ST) / FM)
(
,
SD (IC) = SQRT ((SSY —— (SY*SY) / FM — (STY —— (ST*SY) / FM),**2 / (SST — ST**2
*
/FM))/((FM—2-)*(SST—ST**ZIFMD)
_
SE (IC) = SD (IC) / SQRT (FM— 1.)
,
_
PRINT 906, FK (IC), SD (IC), SE (IC)
_
'
1
'
.
_
906 FORMAT (11—10, 3 HKI=,F8.4, 3 X, 7 HSD (Kl)=, F 10.6, 3 X, 7 HSE (K1) =, F10.‘6)
IF(N—1) 2,2, 3
,
_
,
.
'.
,g
;
_
,
3m=œ+1
,
IF(IC-N)4, 4,5
'
;
5PRINT 907, (TP(I), FK(I), SD(I),
N)
,
__
907 FORMAT (1 H1, 35 H SUMMARY OF STATISTICAL CALCULATIONS /
_
*
12 H0 TEMPERATURE, 9X, 2HK1, 13 X, 2HSD, 14X, 2HSE//
_
…
*
*
(3 X, F7.2, 6X, F8.4, 6X, F10.6, 6X, F10.6))
-
—
PRINT
908
_
;
908 FORMAT (/// 52H STANDARD ERROR OFDIEŒRENŒ;
AND SLOPE“DŒFERENCES
:
,.
-
*
/4X, 43 H FOR TEMPERATURES STANDARDERRORSLOPEU
__; ,
»
Nl = N—— 1
,
51
—
…
DO 8 I = 1, N1
-'
Ï‘
=
L = 1 +1
,
…
**
,
DO 8 I = L, N
,
.
'
SEDIF = SQRT ( SE (I)**2, +, SE
"
SLOPE = (FK(I)-FK(J)) [SEDIF {jj .—îil
f "
"
'
,
_
_
8 PRINT 909, TP (I), TP (D, SEDIF, SLOPE
-
*
909 FORMAT (3 X, F7.2, 3 X, F 7.2,
2 PRINT 910
'
ZÎ ‘
:
—>
'
.
GO TO 20
=
'i'
|
1 '
'
END
—'
'