1  REM  ****  HP BASIC PROGRAM LIBRARY  *******************************
2  REM
3  REM        BICONF:   CONFIDENCE LIMITS
4  REM
5  REM        36691 REV  A   10/73
6  REM
7  REM  ****  CONTRIBUTED PROGRAM  ************************************
100  PRINT 
110  PRINT "CONFIDENCE LIMITS FOR A POPULATION PROPORTION BASED ON"
120  PRINT "THE EXACT BINOMIAL DISTRIBUTION.  WHAT ARE THE VALUES"
130  PRINT "OF X(SUCCESSES), N(SAMPLE SIZE), C(CONFIDENCE COEFF-"
140  PRINT "ICIENT IN PERCENT)";
150  DEF FNA()=.0001*INT(.5+1.E+06*)
160  INPUT Y,N,C
164  PRINT "PLEASE WAIT..........."
170  LET P0=Y/N
180  LET Q0=1-P0
190  LET G0=6*SQR(P0*Q0/N)
200  LET P1=P0-G0
210  IF P1>0 THEN 230
220  LET P1=0
230  LET P1=.1*INT(10*P1)
240  LET P2=P0+G0
250  IF P2<1 THEN 270
260  LET P2=1
270  LET P2=1-.1*INT(10-10*P2)
280  LET C0=.005*(100-C)
290  LET P=P1
300  LET T=Y-1
310  LET J=0
320  LET J=J-1
330  IF J=-7 THEN 390
340  LET P=P+10^J
350  GOSUB 600
360  IF 1-Z<C0 THEN 340
370  LET P=P-10^J
380  GOTO 320
390  LET L=P
400  LET P=P2
410  LET T=Y
420  LET J=0
430  LET J=J-1
440  IF J=-7 THEN 500
450  LET P=P-10^J
460  GOSUB 600
470  IF Z<C0 THEN 450
480  LET P=P+10^J
490  GOTO 430
500  LET U=P
510  PRINT 
520  PRINT 
530  PRINT "   BEST ESTIMATE OF POPULATION PROPORTION (PCT) = ";
540  PRINT FNA(Y/N)
550  PRINT 
560  PRINT "   THE ";C;"PERCENT CONFIDENCE LIMITS ON THE POPULATION"
570  PRINT "   PROPORTION (PCT) ARE ";FNA(L);"AND ";FNA(U)
580  PRINT 
590  STOP 
600  LET Q=1-P
610  LET G=6*SQR(P*Q*N)
620  LET S=INT(N*P-G)
630  IF S>0 THEN 670
640  LET S=0
650  LET X=1
652  X1=LOG(X)
660  GOTO 710
670  LET X=P*N/S
672  X1=LOG(X)
680  FOR I=1 TO S-1
690  X1=X1+LOG(P)+LOG(N-I)-LOG(S-I)
700  NEXT I
710  F1=X1+(N-S)*LOG(Q)
712  F=EXP(F1)
720  LET Z=F
730  FOR I=S+1 TO T
740  LET F=F*P*(N-I+1)/(Q*I)
750  LET Z=Z+F
760  NEXT I
770  RETURN 
780  END 
