1000 REM ***** MODIFIED TO RUN ON PDP-11-20 *****
1100 REM ***** PHIL ZUMSTEG ***** 4/20/72 *****
9000 REM ***** STAT1 ***** STATISTICS *****
9001 REM ***** VERSION 1 ***** 7/31/69/ *****
9002 REM GROUP STATISTICS
9003 LET P=1
9004 LET S,S2,N=0
9005 READ X
9006 IF X=999999. THEN 9011
9007 LET S=S+X
9008 LET S2=S2+X*X
9009 LET N=N+1
9010 GOTO 9005
9011 LET S(P)=S
9012 LET Z(P)=S2
9013 LET N(P)=N
9014 IF P=2 THEN 9017
9015 LET P=2
9016 GOTO 9004
9017 REM NOW WE PRINT THE ANSWERS
9018 PRINT "GROUP","NUMBER","MEAN","VARIANCE","STD. DEV."
9019 PRINT
9020 FOR I=1 TO 2
9021 LET M(I)=S(I)/N(I)
9022 LET V(I)=(N(I)*Z(I)-S(I)*S(I))/N(I)/(N(I)-1)
9023 LET D(I)=SQR(V(I))
9024 PRINT I,N(I),M(I),V(I),D(I)
9025 NEXT I
9026 LET Q=V(1)/N(1)+V(2)/N(2)
9027 LET W=SQR(Q)
9028 LET R=M(1)-M(2)
9029 LET D=V(1)/N(1)/Q
9030 LET D=1/(K*K/(N(1)-1)+(1-K)*(1-K)/(N(2)-1))
9031 PRINT
9032 PRINT "MEAN DIFF.","VAR. DIFF.","STD. DEV. DIFF."
9033 PRINT R,Q,W
9034 PRINT
9035 PRINT "T RATIO",R/W,"ON";D;"DEGREES OF FREEDOM."
9036 GO TO 9999
9900 DATA 5,5,5,5,5,5,999999.,4,6,4,6,4,6,999999.
9999 END