1  REM  ****  HP BASIC PROGRAM LIBRARY  *******************************
2  REM
3  REM        SEVPRO:   CHI-SQUARE TEST
4  REM
5  REM        36719 REV  A   10/73
6  REM
7  REM  ****  CONTRIBUTED PROGRAM  ************************************
11  DATA 5.E+06,5.39828E+06,5.7926E+06,6.17911E+06,6.55422E+06,6.91462E+06,7.25747E+06
12  DATA 7.58036E+06,7.88145E+06,8.1594E+06,8.41345E+06,8.64334E+06,8.8493E+06,9.032E+06
13  DATA 9.19243E+06,9.33193E+06,9.45201E+06,9.55435E+06,9.6407E+06,9.71283E+06,9.7725E+06
14  DATA 9.82136E+06,9.86097E+06,9.89276E+06,9.91803E+06,9.9379E+06,9.95339E+06,9.96533E+06
15  DATA 9.97445E+06,9.98134E+06,9.9865E+06,9.99032E+06,9.99313E+06,9.99517E+06,9.99663E+06
16  DATA 9.99767E+06,9.99841E+06,9.99892E+06,9.99928E+06,9.99952E+06,999683.,999793.
17  DATA 9.99987E+06,9.99992E+06,9.99995E+06,9.99997E+06,9.99998E+06,9.99999E+06,9.99999E+06
18  DIM X[49]
20  DEF FND()=X[]-X[-1]
21  FOR I=1 TO 49
22  READ X[I]
23  NEXT I
100  DATA 1217,45,948,49,1165,33,1121,44
200  DATA 1.E+37,1.E+37
205  DIM Z[100],N[100]
210  LET I=0
215  LET S1=0
220  LET S2=0
225  LET S3=0
230  LET I=I+1
235  READ N[I],Z[I]
240  IF N[I]=1.E+37 THEN 265
245  LET S1=S1+N[I]
250  LET S2=S2+Z[I]
255  LET S3=S3+(Z[I]^2)/N[I]
260  GOTO 230
265  LET M=I-1
270  LET A=(S3-(S2^2)/S1)/(M-1)
275  LET B=(S2/S1)*(1-S2/S1)
280  LET C2=(M-1)*A/B
284  LET F=M-1
285  LET B1=((C2/F)^(1/3)-(1-(2/(9*F))))/SQR(2/(9*F))
290  GOSUB 370
300  PRINT "CHI-SQUARE TEST OF SEVERAL PROPORTIONS:"
305  PRINT 
310  PRINT " SAMPLE    SUCCESSES / TOTAL    PCT SUCCESSES"
315  PRINT 
320  FOR I=1 TO M
325  PRINT "   ";I;"   ";Z[I];"/";N[I],"   ";.001*INT(.5+100000.*Z[I]/N[I])
330  NEXT I
335  PRINT 
340  PRINT "CHI-SQUARED FOR THESE DATA   = ";C2
345  PRINT "CORRESPONDING NORMAL DEVIATE = ";B1
350  PRINT 
360  PRINT " BEING EXCEEDED BY CHANCE ALONE IS";1-B2
365  STOP 
370  IF B1<-4.5 THEN 425
375  IF B1<0 THEN 410
380  IF B1<4.5 THEN 395
385  LET B2=1
390  RETURN 
395  GOSUB 435
400  LET B2=Q
405  RETURN 
410  GOSUB 435
415  LET B2=1-Q
420  RETURN 
425  LET B2=0
430  RETURN 
435  LET Z=10*ABS(B1)
440  LET K=INT(Z)
445  LET D1=Z-K
450  LET Q=X[K]+D1*FND(K+1)+(D1*(D1-1)/2)*(FND(K+2)-FND(K+1))
455  LET Q=.001*INT(.5+.0001*Q)
460  RETURN 
465  PRINT 
470  PRINT "THIS PROGRAM WILL APPLY A CHI-SQUARE TEST TO"
475  PRINT "SEVERAL SAMPLE PROPORTIONS.  TO USE, TYPE:"
480  PRINT 
485  PRINT "   100 DATA  N1,S1, N2,S2, N3,S3, ....."
490  PRINT "   RUN"
495  PRINT 
500  PRINT "WHERE S1 IS THE NUMBER OF SUCCESSES IN SAMPLE 1"
505  PRINT "OF SIZE N1, AND SO ON FOR THE M SAMPLES."
510  END 
