5  REM  HP CONTRIBUTED LIBRARY, REV. B, 2/75
9000  REM ***** BITEST ***** STATISTICS PROGRAM
9010  PRINT 
9020  LET G5=6
9030  LET G0=1
9040  PRINT "THIS PROGRAM MAKES THE NECESSARY CALCULATION FOR A"
9050  PRINT "STATISTICAL TEST OF A BINOMIAL PROPORTION.  WHAT"
9060  PRINT "ARE X(SUCCESSES IN SAMPLE), N(SAMPLE SIZE), AND  P(THE"
9070  PRINT "POPULATION PROPORTION)";
9080  INPUT Y,N,P
9090  LET G5=6
9100  LET G0=1
9110  IF P<Y/N THEN 9140
9120  LET T=Y
9130  GOTO 9150
9140  LET T=Y-1
9150  GOSUB 9450
9160  PRINT 
9170  PRINT "IN SAMPLES OF SIZE ";N;"RANDOMLY SELECTED FROM A"
9180  PRINT "BINOMIAL POPULATION HAVING A TRUE PROPORTION OF";P
9190  PRINT "THE PROBABILITY OF A SAMPLE HAVING ";Y;
9200  IF P<Y/N THEN 9230
9210  PRINT "OR LESS"
9220  GOTO 9250
9230  PRINT "OR MORE"
9240  LET Z=1-Z
9250  PRINT "SUCCESSES IS";
9260  IF Z<5.E-07 THEN 9300
9270  PRINT .000001*INT(1.E+06*Z+.5)
9280  IF G0<0 THEN 9310
9290  STOP 
9300  PRINT " LESS THAN ONE IN A MILLION."
9310  PRINT 
9320  PRINT "THE PROBABILITY FOR EACH POSSIBLE NUMBER OF SUCCESSES MORE"
9330  PRINT "THAN ";G5
9340  PRINT "STANDARD DEVIATIONS AWAY FROM THE MEAN ";N*P;"  HAS NOT ";
9350  PRINT "BEEN"
9360  PRINT "INCLUDED IN COMPUTING THE RESULT ";.000001*INT(1.E+06*Z+.5)
9370  PRINT "TO INCLUDE SUCH A PROBABILTY, THE USER CAN MODIFY THIS ";
9380  PRINT "PROGRAM BY"
9390  PRINT "CHANGING ONE LINE TO READ"
9400  PRINT 
9410  PRINT "9105 LET G5=";G5+1
9420  PRINT 
9430  PRINT "OR SOME VALUE LARGER THAN";G5+1
9440  STOP 
9450  LET Q=1-P
9460  LET G=G5*SQR(P*Q*N)
9470  LET S=INT(N*P-G)
9480  LET P3=LOG(P)
9490  LET Q3=LOG(Q)
9500  IF S>0 THEN 9540
9510  LET S=0
9520  LET X=0
9530  GOTO 9580
9540  LET X=LOG(P*N/S)
9550  FOR I=1 TO S-1
9560  LET X=X+P3+LOG(N-I)-LOG(S-I)
9570  NEXT I
9580  LET F=X+(N-S)*Q3
9590  LET Z=EXP(F)
9600  FOR I=S+1 TO T
9610  LET F=F+P3+LOG(N-I+1)-Q3-LOG(I)
9620  LET Z=Z+EXP(F)
9630  IF I>N*P+G THEN 9660
9640  NEXT I
9650  GOTO 9670
9660  LET G0=-1
9670  RETURN 
9680  END 
