Web pdp-10.trailing-edge.com

Trailing-Edge - PDP-10 Archives - decuslib10-01 - 43,50110/fvalue.sta
There are 2 other files named fvalue.sta in the archive. Click here to see a list.
```100'  NAME--FVALUE
110'
120'  DESCRIPTION--COMPUTES THE EXACT PROBABILITY OF AN F-RATIO WITH
130'  DEGREES OF FREEDOM (M,N).
140'
150'  SOURCE--PROGRAMMED BY VICTOR E MCGEE,DARTMOUTH COLLEGE,
160'  HANOVER, N.H. 03755 AFTER THE FORTRAN PROGRAM GIVEN
170'  IN DONALD J. VELDMAN'S BOOK,"FORTRAN PROGRAMMING FOR THE
180'  BEHAVIORAL SCIENCES",(HOLT,RINEHART, AND WINSTON),1967,P.131
190'
200'  INSTRUCTIONS--TYPE "RUN" AND FOLLOW DIRECTIONS AS THEY APPEAR.
210'
220'  *  *  *  *  *  *  MAIN PROGRAM  *  *  *  *  *  *  *  *  *  *  *
230'
240 PRINT "ENTER F-VALUE,NUMERATOR D.F., AND DENOMINATOR D.F."
250 INPUT G,M,N
260 LET P=1
270 IF G*M*N=0 THEN 500
280 IF G<1 THEN 330
290 LET A=M
300 LET B=N
310 LET F=G
320 GO TO 360
330 LET A=N
340 LET B=M
350 LET F=1/G
360 LET A1=2/(9*A)
370 LET B1=2/(9*B)
380 LET Z=ABS((1-B1)*F^(.333333)-1+A1)
390 LET Z=Z/SQR(B1*F^(.666667)+A1)
400 IF B<4 THEN 440
410 LET P=(1+Z*(.196854+Z*(.115194+Z*(.000344+Z*.019527))))^4
420 LET P=.5/P
430 GO TO 460
440 LET Z=Z*(1+.08*Z^4/B^3)
450 GO TO 410
460 IF G<1 THEN 480
470 GO TO 520
480 LET P=1-P
490 GO TO 520
500 PRINT "YOU HAVE ERRED - INPUT THE 3 VALUES AGAIN"
510 GO TO 250
520 PRINT
530 PRINT
540 PRINT"EXACT PROBABILITY OF F="G;"WITH   ( "M;",  "N;")  D.F."
550 PRINT
560 PRINT "    IS "INT(100000*P)/100000
570 END

```