10X$=SYS(CHR$(7)):R$=LEFT(X$,7):L$=RIGHT(X$,8) 15 REM****************************************************************** 16 REM CMODU CMODU CMODU CMODU CMODU CMODU 17 REM****************************************************************** 20 DIM #1,A1(111),A2(111),A3(111) 25 DIM L(10) 30 OPEN R$ AS FILE 1 70 A1(0)=30 80 X=0 190 GOSUB 5160 200 PRINT L$ 205 PRINT " COMPARISON OF TWO STANDARD DEVIATIONS" 206 PRINT 210 PRINT "THIS MODULE COMPUTES THE PROBABILITY THAT THE RATIO OF TWO" 220 PRINT "STANDARD DEVIATIONS DISTRIBUTED AS INDEPENDENT INVERSE" 230 PRINT "CHIS IS GREATER THAN K, WHERE K CAN TAKE ON ANY POSITIVE" 240 PRINT "VALUE." 242 PRINT 244 PRINT 250 REM***************************************************************** 260 PRINT "INPUT THE PARAMETERS OF THE INVERSE CHI DISTRIBUTION ON THE" 270 PRINT "FIRST STANDARD DEVIATION." 280 REM***************************************************************** 290 PRINT 300PRINT"INPUT THE DEGREES OF FREEDOM."; 310 GOSUB 9000 312 IF O1>2 THEN 320 313PRINT"REENTER. THE DEGREES OF FREEDOM MUST BE GREATER THAN 2." 314 GOTO 310 320 V1=O1 330 PRINT "INPUT THE SCALE PARAMETER."; 340 GOSUB 9000 342 IF O1>0 THEN 350 344 PRINT "REENTER. SCALE PARAMETER MUST BE POSITIVE." 346 GOTO 340 350 L1=O1 360 PRINT 370 REM***************************************************************** 380 PRINT "INPUT THE PARAMETERS OF THE SECOND INVERSE CHI DISTRIBUTION." 390 REM****************************************************************** 400 PRINT 410 PRINT "INPUT THE DEGREES OF FREEDOM"; 420 GOSUB 9000 422 IF O1>2 THEN 430 424 PRINT "REENTER. DEGREES OF FREEEDOM MUST BE GREATER THAN 2." 426 GOTO 420 430 V2=O1 440 PRINT "INPUT THE SCALE PARAMETER."; 450 GOSUB 9000 452 IF O1>0 THEN 460 454 PRINT "REENTER. SCALE PARAMETER MUST BE POSITIVE." 456 GOTO 450 460 L2=O1 465 PRINT L$ 470 PRINT "THIS MODULE COMPUTES AND PRINTS THE PROBABILITY THAT THE RATIO" 480 IF L1/SQR(V1+1)= 1 THEN 650 630 PRINT "REENTER. N MUST BE BETWEEN 1 AND 6 (INCLUSIVE)." 640 GOTO 610 650 IF O1>6 THEN 630 660 N=O1 670 FOR K5=1 TO N 690A$="INPUT VALUE ## " 691PRINT USING A$,K5; 700 GOSUB 9000 710 IF O1>0 THEN 740 720 PRINT "REENTER. K MUST BE POSITIVE." 730 GOTO 700 740 L(K5)=O1 750 NEXT K5 760 PRINT L$ 770 PRINT " COMPARISON OF TWO STANDARD DEVIATIONS" 780 PRINT " SIGMA-1 INVERSE CHI DISTRIBUTIONS SIGMA-2" 790 PRINT "============================================================" 810A$="######.## DEGREES OF FREEDOM #######.##" 811 PRINT USING A$,V1,V2 830A$="######.## SCALE PARAMETER #######.##" 831 PRINT USING A$,L1,L2 834A$="######.## MEAN #######.##" 835 PRINT USING A$,L1/SQR(V1-1.5),L2/SQR(V2-1.5) 836A$="######.## MODE #######.##" 837 PRINT USING A$,L1/SQR(V1+1),L2/SQR(V2+1) 840 PRINT "============================================================" 850 IF B2=2 THEN 900 860 F7=(V1/(L1*L1))/(V2/(L2*L2)) 870 A=V2/2 880 B=V1/2 890 GOTO 910 900 F7=(V2/(L2*L2))/(V1/(L1*L1)) 910 A=V1/2 920 B=V2/2 930 GOSUB 5220 940 FOR K5=1 TO N 950 K=L(K5)^2 960 J1=0 970 J2=(K*F7*A)/(B+K*F7*A) 980 GOSUB 5000 990 P=1-P 1000 IF B2=1 THEN 1050 1020A$=" PROB( SIGMA-2/SIGMA-1 >######.##)=#.##" 1021 PRINT USING A$ ,L(K5),P 1030 GOTO 1060 1050A$=" PROB( SIGMA-1/SIGMA-2 >######.##)=#.##" 1051 PRINT USING A$ ,L(K5),P 1060 NEXT K5 1070 PRINT "============================================================" 1080 PRINT "WHEN YOU ARE READY TO CONTINUE TYPE '1'."; 1090 GOSUB 9000 1100 PRINT L$ 1110 PRINT "TYPE THE NUMBER OF THE OPTION YOU WANT." 1120 PRINT " 1. FURTHER EVALUATE THIS RATIO" 1130 PRINT " 2. EVALUATE A DIFFERENT RATIO" 1140 PRINT " 3. EXIT MODULE" 1150 GOSUB 9000 1160 IF O1=3 THEN 1210 1170 IF O1=2 THEN 1220 1180 IF O1=1 THEN 1240 1190 PRINT "REENTER. INPUT MUST BE THE NUMBER OF AN OPTION." 1200 GOTO 1150 1210 CLOSE 1: CHAIN "RSTRT" 1220 B2=0 1222 PRINT L$ 1230 GOTO 260 1240 PRINT L$ 1250 GOTO 590 5000 REM *************************************************** 5005 REM BETA CDF ROUTINE 5010 REM INPUT A B J2 5015 REM OUTPUT P 5020 REM GOSUB'S TO BE CALLED PRIOR 5160 AND 5220 5025 DIM W(16),O(16) 5030 GOSUB 5340 5035 IF A+B>85 THEN 5280 5040 P=0 5045 C6=0 5050 IF A>1 THEN 5080 5055 C6=A 5060 C7=B 5065 A=C7 5070 B=C6 5075 J2=1-J2 5080 D0=(J2-J1)*.5 5085 D1=(J1+J2)*.5 5090 FOR I1=1 TO 16 5095 D9=D0*O(I1)+D1 5100 IF D9=0 THEN 5115 5105 IF D9=1 THEN 5115 5107 D9=LOG(D9)*(A-1)+LOG(1-D9)*(B-1) 5108 IF D9<-80 THEN 5115 5110 P=P+W(I1)*EXP(D9) 5115 NEXT I1 5120 P=P*F0 5125 P=P*D0 5130 IF C6=0 THEN 5155 5135 A=C6 5140 B=C7 5145 P=1-P 5150 J2=1-J2 5155 RETURN 5160 FOR I1=1 TO 16 5165 READ W(I1),O(I1) 5170 NEXT I1 5175 DATA 2.71525E-02,-.989401 5180 DATA 6.22535E-02,-.944575,9.51585E-02,-.865631 5185 DATA .124629,-.755404,.149596,-.617876 5190 DATA .169156,-.458017,.182603,-.281604,.189451,-9.50125E-02 5195 DATA .189451,9.50125E-02,.182603,.281604,.169156,.458017 5200 DATA .149596,.617876,.124629,.755404 5205 DATA 9.51585E-02,.865631,6.22535E-02,.944575,2.71525E-02 5210 DATA .989401 5215 RETURN 5220 G9=A+B 5225 GOSUB 5850 5230 F0=G0 5235 G9=A 5240 GOSUB 5850 5245 F0=F0-G0 5250 G9=B 5255 GOSUB 5850 5260 F0=F0-G0 5265 IF A+B>85 THEN 5275 5270 F0=EXP(F0) 5275 RETURN 5280 W1=(B*J2)^(1/3) 5285 W2=(A*(1-J2))^(1/3) 5290 GOSUB 5325 5295 I1=P 5300 W1=(B*J1)^(1/3) 5305 W2=(A*(1-J1))^(1/3) 5310 GOSUB 5325 5315 P=I1-P 5320 RETURN 5325 Y3=3*(W1*(1-1/9/B)-W2*(1-1/9/A))/SQR(W1*W1/B+W2*W2/A) 5330 GOSUB 8000 5335 RETURN 5340 REM 2/16/76 CHANGED TO ALL LOG 5345 D2=F0+(A-1)*LOG(J2)+(B-1)*LOG(1-J2) 5350 IF D2<-80 THEN 5370 5355 IF D2>85 THEN 5380 5360 D2=EXP(D2) 5365 RETURN 5370 D2=1.E-37 5375 RETURN 5380 D2=1.E+37 5385 RETURN 5390 REM 5395 REM END OF BETA CDF ROUTINE 5400 REM ******************************************************* 5850 REM **************************************************** 5852 REM LOG GAMMA ROUTINE 5853 REM INPUT G9 5854 REM OUTPUT G0 5860 G5=G9 5863 IF G9 <= 1.E+30 THEN 5872 5866 G0=1.E+38 5869 RETURN 5872 IF G9>1.E-09 THEN 5881 5875 G0=0 5878 RETURN 5881 IF G9<1.E+10 THEN 5890 5884 G0=G9*(LOG(G9)-1) 5887 RETURN 5890 G6=1 5893 IF 18= 0 THEN 8070 8065 P=1-P 8070 X=X1 8075 RETURN 8076 REM 8077 REM END OF NORMAL CDF ROUTINE 8078 REM ********************************************************** 9000 REM--SUBROUTINE THAT DETERMINES IF RESTART HAS BEEN REQUESTED. 9005 INPUT O1 9015 IF O1=-9999 THEN 9025 9020 RETURN 9025 CLOSE 1: CHAIN "RSTRT" 9035 REM*************END ROUTINE 9999 END