10 X$=SYS(CHR$(7)):R$=LEFT(X$,7):L$=RIGHT(X$,8) 20 REM *************************************************************** 30 REM CMOD4 CMOD4 CMOD4 CMOD4 CMOD4 CMOD4 40 REM*************************************************************** 50 REM 60 REM PRIOR FOR STANDARD DEVIATION OF TWO PARAMETER NORMAL 70 REM 80 REM************************************************************ 90 S8=0 100 DIM #1,A1(111),A2(111),A3(111) 110 OPEN R$ AS FILE 1 150 A1(0)=4 160 DATA 1,0,0,0,1,0,0,0,1,1,1,1 170 DIM J(4),S(4,5),W(4,3),V(76),C(72,2) 180 MAT READ W 190 REM****************************************************************** 200 REM N1=DEGREES OF FREEDOM 210 REM l1 =scale factor 220 REM****************************************************************** 230 D9=0 235 S9=0 240 MAT READ C 250 FOR K5=1 TO 72 260 K6=INT(C(K5,1)*.1) 270 V(K5+4)=K6/1000 280 C(K5,1)=C(K5,1)-10*K6 290 NEXT K5 300 PRINT L$ 310 PRINT 320 PRINT "THIS MODULE WILL ASSIST YOU IN FITTING AN INVERSE CHI DISTRIBUTION" 330 PRINT "TO YOUR PRIOR BELIEFS ABOUT THE STANDARD DEVIATION OF A NORMAL" 340 PRINT "DISTRIBUTION" 350 PRINT 360 PRINT "WE BEGIN BY ASKING YOU TO SPECIFY THE 25TH, 50TH AND" 370 PRINT "75TH PERCENTILES OF YOUR PRIOR DISTRIBUTION." 380 PRINT 390 PRINT "SPECIFY 25TH PERCENTILE. YOUR BETTING ODDS ARE 3 TO 1 THAT THE " 400 PRINT "STANDARD DEVIATION IS GREATER THAN THIS VALUE. INPUT 25TH."; 410 GOSUB 9000 420 Q1=O1 430 IF Q1>0 THEN 460 440 PRINT "STANDARD DEVIATION MUST BE POSITIVE." 450 GOTO 390 460 PRINT 470 PRINT "SPECIFY 50TH PERCENTILE. YOUR BETTING ODDS ARE EVEN THAT THE" 480 PRINT "STANDARD DEVIATION IS LESS THAN THIS VALUE. INPUT 50TH."; 490 GOSUB 9000 500 Q2=O1 510 IF Q2>Q1 THEN 560 520 PRINT L$ 530 PRINT "THE 50TH PERCENTILE MUST BE LARGER THAN THE 25TH,AND THE" 540 PRINT "75TH MUST BE LARGER THAN THE 50TH. PLEASE RESPECIFY." 550 GOTO 380 560 Q8=Q2 570 PRINT 580 PRINT "SPECIFY 75TH PERCENTILE. YOUR BETTING ODDS ARE 1 TO 3 THAT THE" 590 PRINT "STANDARD DEVIATION IS GREATER THAN THIS VALUE. INPUT 75TH."; 600 GOSUB 9000 610 Q3=O1 620 IF Q3>Q2 THEN 640 630 GOTO 520 640 PRINT 650A$="FOUR POSSIBLE APPROXIMATE PRIOR DISTRIBUTIONS ARE NOW BEING" 651PRINT A$ 660A$="COMPUTED FOR YOUR CONSIDERATION." 661PRINT A$ 670 FOR K5=1 TO 4 680 R9=1.E+30 690 FOR K6=1 TO 72 700 R=W(K5,1)*(Q2/Q1-C(K6,1))^2 710 R=R+W(K5,2)*(Q3/Q1-C(K6,2))^2 720 R=R+W(K5,3)*(Q2/Q3-C(K6,1)/C(K6,2))^2 730 IF R>R9 THEN 900 740 R9=R 750 IF K6 <= 16 THEN 820 760 IF K6 <= 28 THEN 860 770 IF K6 <= 32 THEN 800 780 V(K5)=K6-22 790 GOTO 870 800 V(K5)=8+(K6-28)*.5 810 GOTO 870 820 V(K5)=2+(K6-1)*.2 830 IF V(K5) <> 2 THEN 850 840 V(K5)=2.01 850 GOTO 870 860 V(K5)=5+(K6-16)*.25 870 J(K5)=K6 880 NEXT K6 890 GOTO 910 900 IF K5=4 THEN 880 910 IF K5=1 THEN 980 920 IF K5=2 THEN 1070 930 IF K5=3 THEN 1200 940 IF K5=4 THEN 1290 950 NEXT K5 960 GOTO 1460 970 REM ASSUME 25TH AND 50TH ARE MORE ACCURATELY ESTIMATED 980 IF J(1) <> 72 THEN 1040 990 K=(Q2-Q1)/.64 1000 GOSUB 2290 1010 J(1)=L1 1020 V(1)=N1 1030 GOTO 1330 1040 J(1)=Q1/V(J(1)+4) 1050 GOTO 1330 1060 REM ASSUME 25TH AND 75TH ARE MOST ACCURATELY ESTIMATED 1070 IF J(2) <> 72 THEN 1160 1080 Q6=Q2 1090 Q2=Q1+(.64/1.34)*(Q3-Q1) 1100 K=(Q2-Q1)/.64 1110 GOSUB 2290 1120 J(2)=L1 1130 V(2)=N1 1140 Q2=Q6 1150 GOTO 1330 1160 J(2)=Q3/(C(J(2),2)*V(J(2)+4)) 1170 GOTO 1330 1180 REM ASSUME 50TH AND 75TH ARE MOST ACCURATELY ESTIMATED 1190 Q8=J(3) 1200 IF J(3) <> 72 THEN 1260 1210 K=(Q3-Q2)/.695 1220 GOSUB 2290 1230 V(3)=N1 1240 J(3)=L1 1250 GOTO 1330 1260 J(3)=Q2/(C(J(3),1)*V(J(3)+4)) 1270 GOTO 1330 1280 REM ASSUME ALL ARE EQUALLY WELL ESTIMATED 1290 IF Q8 <> J(4) THEN 1320 1300 J(4)=Q1/V(J(1)+4) 1310 GOTO 1330 1320 J(4)=Q2/(C(J(4),1)*V(J(4)+4)) 1330 N1=V(K5) 1340 L1=J(K5) 1350 GOSUB 1680 1360 IF K5 <> 1 THEN 1430 1370 PRINT L$ 1380 PRINT "HERE ARE THE PERCENTILES OF FOUR INVERSE CHI DISTRIBUTIONS" 1390 PRINT "FITTED TO YOUR PERCENTILE SPECIFICATIONS." 1400 PRINT 1410A$=" 10TH 25TH 50TH 75TH 90TH" 1411PRINT A$ 1430A$=" # #####.## #####.## #####.## #####.## #####.## 1431 PRINT USING A$ ,K5,S(K5,1),S(K5,2),S(K5,3),S(K5,4),S(K5,5) 1450 GOTO 950 1460 PRINT 1470 PRINT "COMPARE THE PERCENTILES OF THESE DISTRIBUTIONS AND DECIDE" 1480 PRINT "WHICH MOST CLOSELY CORRESPONDS TO YOUR PRIOR BELIEFS. YOU" 1490 PRINT "CAN EITHER TENTATIVELY ACCEPT THIS DISTRIBUTION OR SPECIFY" 1500 PRINT "NEW VALUES FOR THE PERCENTILES." 1510 PRINT 1520 PRINT "IF YOU WANT ONE OF THE DISTRIBUTION TYPE ITS NUMBER."; 1530 PRINT "IF YOU WANT TO RESPECIFY THE PERCENTILES TYPE '0'." 1540 GOSUB 9000 1550 IF O1=0 THEN 1970 1560 GOTO 1600 1570 PRINT "REENTER. INPUT MUST BE 0 OR THE NUMBER OF A DISTRIBUTION."; 1580 GOSUB 9000 1590 GOTO 1550 1600 IF O1=1 THEN 1660 1610 IF O1=2 THEN 1660 1620 IF O1=3 THEN 1660 1630 IF O1=4 THEN 1660 1640 PRINT 1650 GOTO 1570 1660 K5=O1 1670 GOTO 2400 1680 E1=0 1690 E2=.5*L1 1700 P0=.1 1710 GOSUB 1890 1720 E1=X3 1730 E2=L1*.65 1740 P0=.25 1750 GOSUB 1890 1760 P0=.5 1770 E1=X3 1780 E2=L1*.87 1790 GOSUB 1890 1800 E1=X3 1810 E2=1.4*L1 1820 P0=.75 1830 GOSUB 1890 1840 E1=X3 1850 P0=.9 1860 E2=2.2*L1 1870 GOSUB 1890 1880 RETURN 1890 GOSUB 3930 1900 S(K5,INT(P0/.25+1.5))=X3 1910 RETURN 1920 REM **** END OF ROUTINE TO FIND PERCENTILE POINTS ******************* 1930 PRINT 1940 PRINT "RESPECIFY YOUR PERCENTILES." 1950 PRINT 1960 GOTO 1980 1970 PRINT 1980 PRINT "INPUT 25TH"; 1990 GOSUB 9000 2000 C1=O1 2010 PRINT 2020 IF C1>0 THEN 2050 2030 PRINT "REENTER. STANDARD DEVIATION MUST BE POSITIVE." 2040 GOTO 1980 2050 PRINT "INPUT 50TH"; 2060 GOSUB 9000 2070 C2=O1 2080 PRINT 2090 IF C2>C1 THEN 2130 2100 PRINT "REENTER. 50TH MUST BE LARGER THAN 25TH." 2110 GOTO 2050 2120 GOTO 1980 2130 PRINT "INPUT 75TH"; 2140 GOSUB 9000 2150 C3=O1 2160 IF C3>C2 THEN 2190 2170 PRINT "REENTER. 75TH MUST BE LARGER THAN 50TH." 2180 GOTO 2130 2190 Q1=C1 2200 Q2=C2 2210 Q3=C3 2220 GOTO 640 2250 REM ***************************************************************** 2260 REM ***** ROUTINE TO SOLVE FOR INVERSE CHI PARAMETERS ASSUMING" 2270 REM ***** THE VARIANCE AND MEDIAN ARE KNOWN ******* 2280 REM ****************************************************************** 2290 K=K^2 2300 B=-22*K/3-Q2*Q2 2310 A=2*K 2320 C=20*K/3+2*Q2*Q2/3 2330 F=(-B+SQR(B^2-4*A*C))/(2*A) 2333 IF F<500 THEN 2340 2336 F=500 2340 L=Q2^2*(F-2/3) 2350 L=SQR(L) 2360 N1=F 2370 L1=L 2380 RETURN 2390 REM ****** SWITCH FROM PERCENTILE TO HYPOTHETICAL SAMPLE SIZE ******* 2400 M6=S(K5,3) 2410 K7=K5 2420 N1=V(K5) 2430 L1=J(K5) 2440 GOTO 2850 2450 GOSUB 7300 2460 Q5=L1*J1 2470 Q6=L1*J2 2480 RETURN 2490 PRINT "IF YOU DO NOT FEEL THAT THIS HYPOTHETICAL SAMPLE SIZE ( M )" 2500 PRINT "REFLECTS YOUR PRIOR INFORMATION ABOUT THE STANDARD DEVIATION" 2510 PRINT "YOU CAN SPECIFY A DIFFERENT ONE. A DIFFERENT M WILL NOT AFFECT" 2520 PRINT "THE MEDIAN, BUT WILL CHANGE THE HDRS AND OTHER PERCENTILES. A" 2530 PRINT "LARGER M WILL SHORTEN THE HDR INTERVALS, AND A SMALLER M WILL" 2540 PRINT "LENGTHEN THEM." 2550 PRINT 2560 REM 2570 PRINT "IF YOU WANT TO CHANGE M TYPE THE NEW VALUE (GREATER THAN 3)." 2580 PRINT "IF YOU DO NOT WANT TO CHANGE M TYPE '0'." 2590 REM 2600 GOSUB 9000 2610 IF O1=0 THEN 3220 2612 IF O1>9999 THEN 2642 2620 IF O1>3 THEN 2650 2630 PRINT "REENTER. M MUST BE GREATER THAN 3." 2640 GOTO 2590 2642 PRINT"REENTER M MUST NOT BE GREATER THAN 9999." 2643 GOTO 2590 2650 N1=O1-1 2660 S9=1 2670 GOTO 3380 2680 E1=0 2690 E2=.65*L1 2700 P0=.25 2710 GOSUB 3930 2720 S(K7,2)=X3 2730 IF S9=1 THEN 2790 2740 E1=X3 2750 E2=.87*L1 2760 P0=.5 2770 GOSUB 3930 2780 S(K7,3)=X3 2790 E2=1.4*L1 2800 E1=X3 2810 P0=.75 2820 GOSUB 3930 2830 S(K7,4)=X3 2840 K5=K7 2850 PRINT L$ 2860 PRINT "HERE ARE SOME OF THE CHARACTERISTICS OF THE INVERSE CHI" 2870 PRINT "DISTRIBUTION YOU ARE NOW CONSIDERING." 2880 PRINT 2890A$=" HYPOTHETICAL SAMPLE SIZE(M) #####.##" 2891 PRINT USING A$ ,N1+1 2900 M9=L1/SQR(N1+1) 2920 IF D9=0 THEN 2990 2930A$=" DEGREES OF FREEDOM #####.##" 2931 PRINT USING A$ ,N1 2940A$=" SCALE PARAMETER ########.##" 2941 PRINT USING A$ ,L1 2980A$=" MODE #####.##" 2981 PRINT USING A$ ,M9 2990A$=" 25TH PERCENTILE #####.##" 2991 PRINT USING A$ ,S(K5,2) 3030A$=" 50TH (MEDIAN) #####.##" 3031 PRINT USING A$ ,S(K5,3) 3040A$=" 75TH PERCENTILE #####.##" 3041 PRINT USING A$ ,S(K5,4) 3080 G=N1 3090 J5=.5 3100 GOSUB 2450 3110A$=" 50% HDR #####.## - #####.##" 3111 PRINT USING A$ ,Q5,Q6 3120 J5=.75 3130 GOSUB 2450 3140A$=" 75% HDR #####.## - #####.##" 3141 PRINT USING A$ ,Q5,Q6 3150 J5=.95 3160 GOSUB 2450 3170A$=" 95% HDR #####.## - #####.##" 3171 PRINT USING A$ ,Q5,Q6 3180 IF D9=1 THEN 3600 3190 PRINT 3200 IF S9=0 THEN 2490 3210 GOTO 3460 3220 PRINT 3230 IF S8=1 THEN 3280 3240 PRINT "YOU CAN CHANGE THE CENTERING OF THE DISTRIBUTION BY" 3250 PRINT "SPECIFYING A DIFFERENT MEDIAN. THIS WILL NOT AFFECT" 3260 PRINT "THE HYPOTHETICAL SAMPLE SIZE." 3270 PRINT 3280 PRINT "IF YOU WANT TO CHANGE THE MEDIAN TYPE THE NEW VALUE." 3290 PRINT "IF YOU DO NOT TYPE '0'." 3300 GOSUB 9000 3310 IF O1=0 THEN 3550 3320 M6=O1 3330 IF O1>0 THEN 3370 3340 PRINT "REENTER. INPUT MUST BE '0' OR POSITIVE VALUE FOR MEDIAN."; 3350 GOSUB 9000 3360 GOTO 3310 3370 S8=1 3380 L1=1 3390 P0=.5 3400 E1=0 3410 E2=.86 3420 GOSUB 3930 3430 L1=M6/X3 3440 S(K5,3)=L1*X3 3450 GOTO 2680 3460 PRINT "IF YOU WANT TO CHANGE M TYPE NEW VALUE ELSE '0'."; 3470 GOSUB 9000 3480 IF O1=0 THEN 3220 3485 IF O1>9999 THEN 3512 3490 IF O1>3 THEN 3520 3500 PRINT "REENTER. M MUST BE GREATER THAN 3." 3510 GOTO 3470 3512PRINT"REENTER. M MUST BE LESS THAN OR EQUAL TO 9999." 3513 GOTO 3470 3520 N1=O1-1 3530 GOTO 3380 3540 D9=1 3550 PRINT L$ 3560 D9=1 3570 PRINT "HERE ARE SOME OF THE CHARACTERISTICS OF THE INVERSE CHI" 3580 PRINT "DISTRIBUTION FITTED TO YOUR PRIOR BELIEFS ABOUT SIGMA." 3590 GOTO 2880 3600 PRINT 3610 PRINT "THIS COMPLETES THE SPECIFICATION OF A PRIOR DISTRIBUTION" 3620 PRINT "ON SIGMA. IF YOU DO NOT WANT TO FIT A PRIOR DISTRIBTION" 3630 PRINT "ON THE MEAN YOU SHOULD RECORD THE PARAMETERS OF YOUR" 3640 PRINT "PRIOR DISTRIBUTION ON SIGMA (DEGREES AND SCALE)." 3650 PRINT 3660 PRINT "IF YOU WANT TO SPECIFY THE PRIOR ON THE MEAN TYPE '1'." 3670 PRINT "IF YOU DO NOT TYPE '0'." 3680 GOSUB 9000 3690 IF O1=0 THEN 3760 3700 IF O1=1 THEN 3740 3710 PRINT "REENTER. MUST BE '0' OR '1'." 3720 GOSUB 9000 3730 GOTO 3680 3740A2(0)=N1:A2(1)=L1:A2(2)=S(K5,3):A2(3)=0:A2(4)=0:A2(5)=0:A2(6)=0 3750 CLOSE 1: CHAIN "CMOD5" 3760 CLOSE 1: CHAIN "RSTRT" 3930 REM********************************************************* 3945 IF E1 <> 0 THEN 3950 3946 P8=0 3947 GOTO 3980 3950 X3=E1 3960 GOSUB 4020 3970 P8=P 3980 X3=E2 3990 GOSUB 4020 4000 P9=P 4004 X3=.5*(E1+E2) 4006 GOSUB 4020 4010 GOTO 4080 4020 X=L1*L1/X3/X3 4030 G=N1 4040 GOSUB 5775 4050 GOSUB 5500 4060 P=1-P 4070 RETURN 4080 IF ABS(P-P0)<.001 THEN 4200 4090 IF P.15 THEN 4140 4138 P=.15 4140 X3=P*(E2-E1)+E1 4150 GOSUB 4020 4160 GOTO 4080 4170 E1=X3 4180 P8=P 4190 GOTO 4120 4200 RETURN 4210 DATA 6011.41,2.1944 4220 DATA 5731.39,2.0971 4230 DATA 5491.36,2.0152 4240 DATA 5281.35,1.9474 4250 DATA 5091.33,1.8907 4260 DATA 4931.32,1.8392 4270 DATA 4781.31,1.7965 4280 DATA 4641.29,1.7605 4290 DATA 4521.29,1.7292 4300 DATA 4411.27,1.7003 4310 DATA 4301.27,1.6737 4320 DATA 4211.26,1.6487 4330 DATA 4121.25,1.627 4340 DATA 4031.25,1.6116 4350 DATA 3961.24,1.5896 4360 DATA 3881.23,1.5736 4370 DATA 3801.23,1.553 4380 DATA 3721.22,1.5354 4390 DATA 3641.22,1.5211 4400 DATA 3571.21,1.5041 4410 DATA 3501.21,1.4922 4420 DATA 3441.2,1.4809 4430 DATA 3381.2,1.4675 4440 DATA 3321.19,1.4581 4450 DATA 3271.19,1.4478 4460 DATA 3221.19,1.4364 4470 DATA 3171.18,1.4277 4480 DATA 3121.18,1.4197 4490 DATA 3041.17,1.4029 4500 DATA 2961.17,1.3888 4510 DATA 2891.16,1.3767 4520 DATA 2821.16,1.3651 4530 DATA 2701.15,1.3454 4540 DATA 2591.14,1.3264 4550 DATA 2501.14,1.3125 4560 DATA 2411.13,1.297 4570 DATA 2341.13,1.2857 4580 DATA 2271.12,1.275 4590 DATA 2201.12,1.2652 4600 DATA 2151.12,1.2577 4610 DATA 2091.11,1.2491 4620 DATA 2041.11,1.242 4630 DATA 2001.11,1.2354 4640 DATA 1961.1,1.2279 4650 DATA 1911.1,1.2239 4660 DATA 1881.1,1.2179 4670 DATA 1841.1,1.213 4680 DATA 1811.1,1.2081 4690 DATA 1781.09,1.2036 4700 DATA 1751.09,1.2001 4710 DATA 1721.09,1.1963 4720 DATA 1691.09,1.1922 4730 DATA 1661.09,1.1886 4740 DATA 1641.09,1.1856 4750 DATA 1621.09,1.1822 4760 DATA 1591.08,1.1791 4770 DATA 1571.08,1.1765 4780 DATA 1551.08,1.1733 4790 DATA 1531.08,1.1717 4800 DATA 1511.08,1.1681 4810 DATA 1491.08,1.1666 4820 DATA 1481.08,1.1641 4830 DATA 1461.08,1.1618 4840 DATA 1441.08,1.1603 4850 DATA 1431.07,1.1578 4860 DATA 1411.07,1.1557 4870 DATA 1401.07,1.1543 4880 DATA 1381.07,1.1519 4890 DATA 1371.07,1.1504 4900 DATA 1351.07,1.1482 4910 DATA 1341.07,1.147 4920 DATA 1331.06,1.1452 5500 REM ******************************************************** 5501 REM CHI-SQUARE CDF ROUTINE-LOWER TAIL 5502 REM INPUT G X 5503 REM OUTPUT P 5504 REM PRIOR GOSUB 5775 5505 REM 5506 REM 5510 IF G>30 THEN 5795 5515 T0=.5*G-1 5520 X=X*.5 5525 T=81.4983 5530 P=5.575E-35*(T+X)^T0 5535 T=69.9622 5540 P=P+4.0883E-30*(T+X)^T0 5545 T=61.0585 5550 P=P+2.45182E-26*(T+X)^T0 5555 T=53.6086 5560 P=P+3.60577E-23*(T+X)^T0 5565 T=47.1531 5570 P=P+2.01052E-20*(T+X)^T0 5575 T=41.4517 5580 P=P+5.35019E-18*(T+X)^T0 5585 T=36.3584 5590 P=P+7.8198E-16*(T+X)^T0 5595 T=31.776 5600 P=P+6.89418E-14*(T+X)^T0 5605 T=27.6359 5610 P=P+3.91774E-12*(T+X)^T0 5615 T=23.8873 5620 P=P+1.50701E-10*(T+X)^T0 5625 T=20.4915 5630 P=P+4.07286E-09*(T+X)^T0 5635 T=17.418 5640 P=P+7.96081E-08*(T+X)^T0 5645 T=14.6427 5650 P=P+1.15132E-06*(T+X)^T0 5655 T=12.1461 5660 P=P+1.25447E-05*(T+X)^T0 5665 T=9.9121 5670 P=P+1.04461E-04*(T+X)^T0 5675 T=7.92754 5680 P=P+6.72163E-04*(T+X)^T0 5685 T=6.18154 5690 P=P+3.36935E-03*(T+X)^T0 5695 T=4.66508 5700 P=P+.013226*(T+X)^T0 5705 T=3.37077 5710 P=P+4.07325E-02*(T+X)^T0 5715 T=2.29256 5720 P=P+9.81663E-02*(T+X)^T0 5725 T=1.4256 5730 P=P+.183323*(T+X)^T0 5735 T=.766097 5740 P=P+.258807*(T+X)^T0 5745 T=.311239 5750 P=P+.258774*(T+X)^T0 5755 T=5.90199E-02 5760 P=P+.142812*(T+X)^T0 5761 IF X<80 THEN 5765 5762 P=1 5763 RETURN 5765 P=1-(T1*P*EXP(-X)) 5770 RETURN 5775 G9=.5*G 5778 IF G>30 THEN 5790 5780 GOSUB 5850 5785 T1=1/EXP(G0) 5790 RETURN 5795 X2=((X/G)^(1/3)-(1-2/9/G))/SQR(2/9/G) 5800 Y3=X2 5805 GOSUB 8000 5810 RETURN 5812 REM 5813 REM END OF CHI-SQUARE CDF ROUTINE 5815 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 1830 THEN 7307 7306 GOSUB 5775 7307 J8=1 7308 GOSUB 7337 7309 X=1/(J1*J1) 7310 GOSUB 5500 7311 J3=P 7312 X=1/(J2*J2) 7313 GOSUB 5500 7314 J3=J3-P 7315 IF ABS(J3-J5)>.0001 THEN 7317 7316 RETURN 7317 IF J3>J5 THEN 7320 7318 J8=J8+1 7319 GOTO 7308 7320 J9=J8-1 7321 J0=J8 7322 J8=(J0+J9)/2 7323 GOSUB 7337 7324 X=1/(J1*J1) 7325 GOSUB 5500 7326 J3=P 7327 X=1/(J2*J2) 7328 GOSUB 5500 7329 J3=J3-P 7330 IF ABS(J3-J5)>.0001 THEN 7332 7331 RETURN 7332 IF J3>J5 THEN 7335 7333 J9=J8 7334 GOTO 7322 7335 J0=J8 7336 GOTO 7322 7337 J=J8*(1+EXP(-2*J8/(G+1)))/(1-EXP(-2*J8/(G+1))) 7338 J1=1/SQR(J+J8) 7339 J2=1/SQR(J-J8) 7340 RETURN 7341 REM 7342 REM END OF ROUTINE FOR INVERSE CHI HDR'S 7343 REM************************************************** 8000 REM ********************************************************** 8001 REM ROUTINE CALCULATES THE CDF FOR NORMAL DISTRIBUTION 8002 REM INPUT Y3 8003 REM OUTPUT P 8004 REM 8005 Y4=ABS(Y3) 8010 X1=X 8015 X=Y3 8020 T=1/(1+.231642*Y4) 8021 IF X*X/2<80 THEN 8025 8022 D=0 8023 GOTO 8030 8025 D=.398942*EXP(-X*X/2) 8030 C1=1.33027 8035 C2=1.82126 8040 C3=1.78148 8045 C4=.356564 8050 C5=.319382 8055 P=1-D*T*((((C1*T-C2)*T+C3)*T-C4)*T+C5) 8060 IF X >= 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