data neuse; input year mo $ day gold kins; t+1; group=floor((t-1)/90)+1; ** groups of 90 days (quarters) **; season = 1 + mod((group-1),4); lgold=log(gold); lkins=log(kins); date=intnx('day','31dec70'd,_n_); format date monyy. ; A=2*4*atan(1)/365.25; s1 = sin(A*t); c1=cos(A*t); s2=sin(2*A*t); c2=cos(2*A*t); cards; 1971 JAN 1 1130 1550 1971 JAN 2 1110 1450 1971 JAN 3 1190 1380 1971 JAN 4 1350 1410 1971 JAN 5 1490 1520 1971 JAN 6 1730 1640 1971 JAN 7 2770 1830 1971 JAN 8 3500 2620 1971 JAN 9 4000 3080 1971 JAN 10 4400 3610 1971 JAN 11 4900 4110 1971 JAN 12 5420 4560 1971 JAN 13 5480 5000 1971 JAN 14 5130 5380 1971 JAN 15 4560 5820 1971 JAN 16 3870 5990 1971 JAN 17 3230 5800 1971 JAN 18 2750 5230 1971 JAN 19 2470 4490 1971 JAN 20 2200 3790 1971 JAN 21 1980 3210 1971 JAN 22 1800 2760 1971 JAN 23 1690 2460 1971 JAN 24 1900 2310 1971 JAN 25 2540 2580 1971 JAN 26 3130 3040 1971 JAN 27 3560 3420 1971 JAN 28 3670 3780 1971 JAN 29 3370 4040 1971 JAN 30 2880 4150 1971 JAN 31 2750 4220 1971 FEB 1 2820 4140 1971 FEB 2 3130 4050 1971 FEB 3 3320 3970 1971 FEB 4 3250 4050 1971 FEB 5 3300 4160 1971 FEB 6 3740 4270 1971 FEB 7 4660 4430 1971 FEB 8 5560 4740 1971 FEB 9 6480 5310 1971 FEB 10 7200 5950 1971 FEB 11 7870 6480 1971 FEB 12 8630 7000 1971 FEB 13 9510 7830 1971 FEB 14 10300 8670 1971 FEB 15 10700 9350 1971 FEB 16 10600 10100 1971 FEB 17 10300 10700 1971 FEB 18 9770 11000 1971 FEB 19 8970 11100 1971 FEB 20 7570 10800 1971 FEB 21 5990 10200 1971 FEB 22 4290 9380 1971 FEB 23 4630 8450 1971 FEB 24 5540 7460 1971 FEB 25 6440 6860 1971 FEB 26 7040 6670 1971 FEB 27 7070 6860 1971 FEB 28 6070 7200 1971 MAR 1 4750 7460 1971 MAR 2 4230 7390 1971 MAR 3 5560 7370 1971 MAR 4 7790 8220 1971 MAR 5 8410 9170 1971 MAR 6 9400 9300 1971 MAR 7 10400 9510 1971 MAR 8 11400 10000 1971 MAR 9 12200 10700 1971 MAR 10 12300 11500 1971 MAR 11 11800 12400 1971 MAR 12 10700 12900 1971 MAR 13 8570 13000 1971 MAR 14 6260 12500 1971 MAR 15 4080 11300 1971 MAR 16 3500 9070 1971 MAR 17 3410 7010 1971 MAR 18 3530 5530 1971 MAR 19 3560 4930 1971 MAR 20 3280 4520 1971 MAR 21 2810 4280 1971 MAR 22 2520 3940 1971 MAR 23 2350 3480 1971 MAR 24 2230 3090 1971 MAR 25 2040 2850 1971 MAR 26 2080 2710 1971 MAR 27 2470 2760 1971 MAR 28 2940 2980 1971 MAR 29 3300 3240 1971 MAR 30 3900 3560 1971 MAR 31 4510 3920 1971 APR 1 4930 4320 1971 APR 2 5190 4770 1971 APR 3 5360 5200 1971 APR 4 5080 5560 1971 APR 5 3930 5770 1971 APR 6 4450 5980 1971 APR 7 5280 5950 1971 APR 8 6020 6110 1971 APR 9 6580 6390 1971 APR 10 6880 6670 1971 APR 11 6930 6940 1971 APR 12 6850 7130 1971 APR 13 6640 7250 1971 APR 14 5700 7260 1971 APR 15 3930 7140 1971 APR 16 2700 6790 1971 APR 17 2350 5530 1971 APR 18 2060 4160 1971 APR 19 1840 3160 1971 APR 20 1690 2660 1971 APR 21 1560 2260 1971 APR 22 1600 2060 1971 APR 23 1630 1990 1971 APR 24 1620 2020 1971 APR 25 1670 2010 1971 APR 26 1620 1990 1971 APR 27 1530 1970 1971 APR 28 1700 1980 1971 APR 29 2030 2220 1971 APR 30 2670 2390 1971 MAY 1 3070 2660 1971 MAY 2 3040 3000 1971 MAY 3 2630 3230 1971 MAY 4 2130 3190 1971 MAY 5 1850 2830 1971 MAY 6 1600 2350 1971 MAY 7 1360 2000 1971 MAY 8 1240 1730 1971 MAY 9 1180 1540 1971 MAY 10 1100 1440 1971 MAY 11 1110 1350 1971 MAY 12 999 1310 1971 MAY 13 1120 1300 1971 MAY 14 1220 1460 1971 MAY 15 2160 1560 1971 MAY 16 2760 2210 1971 MAY 17 3020 2750 1971 MAY 18 3550 3070 1971 MAY 19 3910 3330 1971 MAY 20 4110 3620 1971 MAY 21 4280 4100 1971 MAY 22 4240 4430 1971 MAY 23 3620 4450 1971 MAY 24 1900 4370 1971 MAY 25 1370 3520 1971 MAY 26 1180 2110 1971 MAY 27 1030 1520 1971 MAY 28 902 1310 1971 MAY 29 874 1410 1971 MAY 30 1070 1270 1971 MAY 31 1910 1270 1971 JUN 1 2580 1760 1971 JUN 2 3000 2280 1971 JUN 3 3170 2700 1971 JUN 4 3350 2940 1971 JUN 5 2880 3110 1971 JUN 6 1490 3180 1971 JUN 7 1020 2320 1971 JUN 8 846 1360 1971 JUN 9 782 1060 1971 JUN 10 711 946 1971 JUN 11 525 872 1971 JUN 12 630 732 1971 JUN 13 753 764 1971 JUN 14 634 856 1971 JUN 15 634 804 1971 JUN 16 651 753 1971 JUN 17 623 773 1971 JUN 18 623 763 1971 JUN 19 641 757 1971 JUN 20 683 756 1971 JUN 21 634 779 1971 JUN 22 578 766 1971 JUN 23 606 844 1971 JUN 24 571 975 1971 JUN 25 588 826 1971 JUN 26 585 758 1971 JUN 27 634 747 1971 JUN 28 585 762 1971 JUN 29 477 729 1971 JUN 30 407 649 1971 JUL 1 394 580 1971 JUL 2 519 561 1971 JUL 3 842 675 1971 JUL 4 606 972 1971 JUL 5 560 821 1971 JUL 6 439 722 1971 JUL 7 616 680 1971 JUL 8 749 897 1971 JUL 9 854 1030 1971 JUL 10 918 999 1971 JUL 11 854 1110 1971 JUL 12 1030 1080 1971 JUL 13 1190 1140 1971 JUL 14 942 1260 1971 JUL 15 627 1160 1971 JUL 16 532 855 1971 JUL 17 442 708 1971 JUL 18 400 627 1971 JUL 19 314 552 1971 JUL 20 362 506 1971 JUL 21 445 524 1971 JUL 22 455 588 1971 JUL 23 690 596 1971 JUL 24 752 690 1971 JUL 25 602 825 1971 JUL 26 458 768 1971 JUL 27 388 722 1971 JUL 28 506 688 1971 JUL 29 602 699 1971 JUL 30 1350 720 1971 JUL 31 1480 1250 1971 AUG 1 1370 1550 1971 AUG 2 1140 1580 1971 AUG 3 914 1410 1971 AUG 4 814 1170 1971 AUG 5 763 987 1971 AUG 6 588 921 1971 AUG 7 850 784 1971 AUG 8 1090 828 1971 AUG 9 926 1080 1971 AUG 10 606 1050 1971 AUG 11 529 830 1971 AUG 12 483 670 1971 AUG 13 423 625 1971 AUG 14 350 574 1971 AUG 15 347 511 1971 AUG 16 445 535 1971 AUG 17 543 718 1971 AUG 18 602 995 1971 AUG 19 609 1090 1971 AUG 20 806 954 1971 AUG 21 874 968 1971 AUG 22 763 1060 1971 AUG 23 676 997 1971 AUG 24 496 920 1971 AUG 25 648 773 1971 AUG 26 858 779 1971 AUG 27 958 1030 1971 AUG 28 822 1330 1971 AUG 29 954 1360 1971 AUG 30 1040 1300 1971 AUG 31 1430 1300 1971 SEP 1 1520 1480 1971 SEP 2 1000 1670 1971 SEP 3 641 1370 1971 SEP 4 557 998 1971 SEP 5 435 807 1971 SEP 6 338 715 1971 SEP 7 296 603 1971 SEP 8 376 553 1971 SEP 9 474 570 1971 SEP 10 423 654 1971 SEP 11 359 672 1971 SEP 12 376 703 1971 SEP 13 407 939 1971 SEP 14 439 926 1971 SEP 15 1620 728 1971 SEP 16 2310 1340 1971 SEP 17 2590 2030 1971 SEP 18 2370 2420 1971 SEP 19 1090 2480 1971 SEP 20 651 1660 1971 SEP 21 451 1060 1971 SEP 22 496 883 1971 SEP 23 464 697 1971 SEP 24 557 648 1971 SEP 25 539 634 1971 SEP 26 529 658 1971 SEP 27 416 641 1971 SEP 28 376 578 1971 SEP 29 362 504 1971 SEP 30 567 655 1971 OCT 1 790 1840 1971 OCT 2 906 2360 1971 OCT 3 2400 2680 1971 OCT 4 3470 3080 1971 OCT 5 4040 3580 1971 OCT 6 4420 3990 1971 OCT 7 4490 4350 1971 OCT 8 4420 4560 1971 OCT 9 4400 4700 1971 OCT 10 4630 4890 1971 OCT 11 4920 5130 1971 OCT 12 5110 5470 1971 OCT 13 5240 5680 1971 OCT 14 5180 5710 1971 OCT 15 5060 5710 1971 OCT 16 5040 5660 1971 OCT 17 5020 5600 1971 OCT 18 4820 5530 1971 OCT 19 3420 5440 1971 OCT 20 2070 5240 1971 OCT 21 1650 4250 1971 OCT 22 1530 3410 1971 OCT 23 2330 3740 1971 OCT 24 3470 4440 1971 OCT 25 4600 5130 1971 OCT 26 5330 5380 1971 OCT 27 6140 5630 1971 OCT 28 6870 5990 1971 OCT 29 7470 6420 1971 OCT 30 7830 6860 1971 OCT 31 7890 7330 1971 NOV 1 7720 7780 1971 NOV 2 7580 8100 1971 NOV 3 7630 8320 1971 NOV 4 7580 8280 1971 NOV 5 6600 8170 1971 NOV 6 5120 8050 1971 NOV 7 3950 7760 1971 NOV 8 3160 7030 1971 NOV 9 2630 5650 1971 NOV 10 2090 4360 1971 NOV 11 1710 3410 1971 NOV 12 1590 2640 1971 NOV 13 1490 2310 1971 NOV 14 1490 2070 1971 NOV 15 1370 1980 1971 NOV 16 1270 1880 1971 NOV 17 1200 1750 1971 NOV 18 1180 1650 1971 NOV 19 1140 1600 1971 NOV 20 1150 1580 1971 NOV 21 1120 1580 1971 NOV 22 1100 1530 1971 NOV 23 1030 1490 1971 NOV 24 1010 1440 1971 NOV 25 1120 1490 1971 NOV 26 1180 1610 1971 NOV 27 1790 1640 1971 NOV 28 2800 1970 1971 NOV 29 3210 2730 1971 NOV 30 3480 3050 1971 DEC 1 3610 3370 1971 DEC 2 3480 3590 1971 DEC 3 2640 3730 1971 DEC 4 2210 3700 1971 DEC 5 1940 3340 1971 DEC 6 1910 2890 1971 DEC 7 1940 2600 1971 DEC 8 2000 2510 1971 DEC 9 1990 2490 1971 DEC 10 1930 2460 1971 DEC 11 2130 2090 1971 DEC 12 2130 2450 1971 DEC 13 1890 3360 1971 DEC 14 1690 4220 1971 DEC 15 1730 5840 1971 DEC 16 1850 2100 1971 DEC 17 1860 2200 1971 DEC 18 1550 2280 1971 DEC 19 1520 2070 1971 DEC 20 1630 1930 1971 DEC 21 1790 2030 1971 DEC 22 1960 2180 1971 DEC 23 2340 5270 1971 DEC 24 2380 4300 1971 DEC 25 2200 2650 1971 DEC 26 1960 2600 1971 DEC 27 1750 2430 1971 DEC 28 1620 2200 1971 DEC 29 1540 2030 1971 DEC 30 1440 1910 1971 DEC 31 1380 1810 1972 JAN 1 1390 1720 1972 JAN 2 1330 1690 1972 JAN 3 1260 3430 1972 JAN 4 1180 1600 1972 JAN 5 1180 1540 1972 JAN 6 1180 1530 1972 JAN 7 1180 1510 1972 JAN 8 1190 1490 1972 JAN 9 1190 1490 1972 JAN 10 1260 1540 1972 JAN 11 2150 1970 1972 JAN 12 3380 2650 1972 JAN 13 5260 3200 1972 JAN 14 6700 3980 1972 JAN 15 7290 5410 1972 JAN 16 7550 6740 1972 JAN 17 7600 7450 1972 JAN 18 7400 7820 1972 JAN 19 6000 7970 1972 JAN 20 4500 7910 1972 JAN 21 3600 7580 1972 JAN 22 3100 6790 1972 JAN 23 2800 5620 1972 JAN 24 2700 4540 1972 JAN 25 2600 3920 1972 JAN 26 2540 3560 1972 JAN 27 2310 3300 1972 JAN 28 2110 3040 1972 JAN 29 1960 2790 1972 JAN 30 1890 2610 1972 JAN 31 1860 2530 1972 FEB 1 1850 2470 1972 FEB 2 2670 2740 1972 FEB 3 4080 3270 1972 FEB 4 5340 3950 1972 FEB 5 5850 4720 1972 FEB 6 6420 5440 1972 FEB 7 6920 6010 1972 FEB 8 7320 6460 1972 FEB 9 7660 6840 1972 FEB 10 7910 7150 1972 FEB 11 7990 7480 1972 FEB 12 7820 7800 1972 FEB 13 7660 8460 1972 FEB 14 7310 8880 1972 FEB 15 7090 8900 1972 FEB 16 7130 8540 1972 FEB 17 7280 8160 1972 FEB 18 7400 7960 1972 FEB 19 7610 8080 1972 FEB 20 7740 8170 1972 FEB 21 7780 8270 1972 FEB 22 7700 8280 1972 FEB 23 7600 8240 1972 FEB 24 7400 8180 1972 FEB 25 6800 8050 1972 FEB 26 5800 7850 1972 FEB 27 5000 7490 1972 FEB 28 4000 6880 1972 FEB 29 3800 6090 1972 MAR 1 3580 5320 1972 MAR 2 3360 4750 1972 MAR 3 2990 4430 1972 MAR 4 2660 4140 1972 MAR 5 2570 3750 1972 MAR 6 2340 3400 1972 MAR 7 2270 3130 1972 MAR 8 2190 2910 1972 MAR 9 2030 2780 1972 MAR 10 1900 2620 1972 MAR 11 1790 2460 1972 MAR 12 1740 2320 1972 MAR 13 1640 2230 1972 MAR 14 1570 2150 1972 MAR 15 1610 2060 1972 MAR 16 1580 2060 1972 MAR 17 1910 2230 1972 MAR 18 2990 2510 1972 MAR 19 3960 2860 1972 MAR 20 4580 3360 1972 MAR 21 4610 3920 1972 MAR 22 4140 4530 1972 MAR 23 3650 4990 1972 MAR 24 3740 4990 1972 MAR 25 3800 4660 1972 MAR 26 3510 4430 1972 MAR 27 3040 4290 1972 MAR 28 2710 4050 1972 MAR 29 2600 3650 1972 MAR 30 2550 3260 1972 MAR 31 2500 3120 1972 APR 1 2550 3140 1972 APR 2 2820 3150 1972 APR 3 3040 3190 1972 APR 4 3090 3310 1972 APR 5 2890 3400 1972 APR 6 2630 3370 1972 APR 7 2320 3210 1972 APR 8 2210 3000 1972 APR 9 2300 2800 1972 APR 10 2570 2720 1972 APR 11 2840 2800 1972 APR 12 3080 2980 1972 APR 13 2840 3210 1972 APR 14 2840 3330 1972 APR 15 2940 3250 1972 APR 16 3010 3210 1972 APR 17 2920 3250 1972 APR 18 2600 3230 1972 APR 19 2200 3160 1972 APR 20 1850 2950 1972 APR 21 1650 2510 1972 APR 22 1560 2150 1972 APR 23 1590 2000 1972 APR 24 1730 1940 1972 APR 25 2140 2080 1972 APR 26 3040 2390 1972 APR 27 3500 2800 1972 APR 28 3220 3200 1972 APR 29 2480 3390 1972 APR 30 2030 3190 1972 MAY 1 1660 2550 1972 MAY 2 1440 2170 1972 MAY 3 1950 1890 1972 MAY 4 3000 2370 1972 MAY 5 3850 3400 1972 MAY 6 4460 3690 1972 MAY 7 4620 4030 1972 MAY 8 4300 4390 1972 MAY 9 3940 4620 1972 MAY 10 2410 4590 1972 MAY 11 1940 3940 1972 MAY 12 2010 2910 1972 MAY 13 1890 2530 1972 MAY 14 1480 2350 1972 MAY 15 1320 2210 1972 MAY 16 1620 1970 1972 MAY 17 2540 1980 1972 MAY 18 3240 2490 1972 MAY 19 3680 2940 1972 MAY 20 3850 3370 1972 MAY 21 4250 3670 1972 MAY 22 4810 4020 1972 MAY 23 5350 4380 1972 MAY 24 5700 4780 1972 MAY 25 5710 5200 1972 MAY 26 5690 5640 1972 MAY 27 5630 5920 1972 MAY 28 5370 6010 1972 MAY 29 4830 5980 1972 MAY 30 3230 5860 1972 MAY 31 2280 5530 1972 JUN 1 2300 4520 1972 JUN 2 2310 3990 1972 JUN 3 2240 3420 1972 JUN 4 1900 3000 1972 JUN 5 1490 2590 1972 JUN 6 1140 2110 1972 JUN 7 1130 1670 1972 JUN 8 1080 1510 1972 JUN 9 982 1430 1972 JUN 10 866 1330 1972 JUN 11 806 1210 1972 JUN 12 690 1130 1972 JUN 13 613 1020 1972 JUN 14 571 932 1972 JUN 15 592 856 1972 JUN 16 567 868 1972 JUN 17 550 844 1972 JUN 18 506 816 1972 JUN 19 432 779 1972 JUN 20 515 786 1972 JUN 21 958 1000 1972 JUN 22 1300 1530 1972 JUN 23 2300 1790 1972 JUN 24 2980 2240 1972 JUN 25 3330 2680 1972 JUN 26 3470 3030 1972 JUN 27 3520 3260 1972 JUN 28 3640 3400 1972 JUN 29 3800 3520 1972 JUN 30 2570 3700 1972 JUL 1 2080 3660 1972 JUL 2 2080 2800 1972 JUL 3 1510 2480 1972 JUL 4 1180 1990 1972 JUL 5 978 1620 1972 JUL 6 794 1420 1972 JUL 7 890 1190 1972 JUL 8 1010 1140 1972 JUL 9 862 1220 1972 JUL 10 728 1140 1972 JUL 11 686 1040 1972 JUL 12 978 1440 1972 JUL 13 1390 1820 1972 JUL 14 1750 1960 1972 JUL 15 2280 2070 1972 JUL 16 2330 2400 1972 JUL 17 2220 2540 1972 JUL 18 1690 2580 1972 JUL 19 1300 2400 1972 JUL 20 1140 1880 1972 JUL 21 1090 1540 1972 JUL 22 1010 1420 1972 JUL 23 866 1320 1972 JUL 24 739 1210 1972 JUL 25 595 1120 1972 JUL 26 509 960 1972 JUL 27 480 936 1972 JUL 28 539 876 1972 JUL 29 700 1090 1972 JUL 30 1420 1440 1972 JUL 31 1340 1590 1972 AUG 1 2480 1680 1972 AUG 2 2400 2480 1972 AUG 3 2110 2650 1972 AUG 4 1850 2490 1972 AUG 5 1390 2250 1972 AUG 6 1210 1860 1972 AUG 7 1900 1590 1972 AUG 8 2220 2020 1972 AUG 9 1990 2390 1972 AUG 10 1610 2370 1972 AUG 11 1260 2090 1972 AUG 12 1070 1690 1972 AUG 13 760 1420 1972 AUG 14 637 1160 1972 AUG 15 770 996 1972 AUG 16 942 1010 1972 AUG 17 826 1120 1972 AUG 18 862 1220 1972 AUG 19 902 1280 1972 AUG 20 2140 1270 1972 AUG 21 2160 2000 1972 AUG 22 1920 2290 1972 AUG 23 1390 2210 1972 AUG 24 982 1830 1972 AUG 25 756 1400 1972 AUG 26 648 1130 1972 AUG 27 588 976 1972 AUG 28 648 1260 1972 AUG 29 866 1200 1972 AUG 30 990 1460 1972 AUG 31 1030 1410 1972 SEP 1 958 1370 1972 SEP 2 1360 1700 1972 SEP 3 1190 2180 1972 SEP 4 1110 2030 1972 SEP 5 966 1850 1972 SEP 6 846 1740 1972 SEP 7 914 1520 1972 SEP 8 1300 1370 1972 SEP 9 1120 1550 1972 SEP 10 886 1500 1972 SEP 11 616 1290 1972 SEP 12 461 1040 1972 SEP 13 451 844 1972 SEP 14 439 776 1972 SEP 15 407 737 1972 SEP 16 400 716 1972 SEP 17 451 748 1972 SEP 18 461 786 1972 SEP 19 419 779 1972 SEP 20 529 730 1972 SEP 21 644 751 1972 SEP 22 564 872 1972 SEP 23 630 844 1972 SEP 24 669 836 1972 SEP 25 543 884 1972 SEP 26 407 836 1972 SEP 27 400 706 1972 SEP 28 480 720 1972 SEP 29 525 1140 1972 SEP 30 543 1360 1972 OCT 1 799 1140 1972 OCT 2 1810 1120 1972 OCT 3 2170 1760 1972 OCT 4 2150 2170 1972 OCT 5 1510 2330 1972 OCT 6 1470 2040 1972 OCT 7 2270 1850 1972 OCT 8 3300 2200 1972 OCT 9 3630 2730 1972 OCT 10 3710 3120 1972 OCT 11 3970 3360 1972 OCT 12 4330 3550 1972 OCT 13 4650 3740 1972 OCT 14 3810 3980 1972 OCT 15 1890 4180 1972 OCT 16 1210 3530 1972 OCT 17 961 2120 1972 OCT 18 839 1470 1972 OCT 19 790 1250 1972 OCT 20 741 1180 1972 OCT 21 720 1090 1972 OCT 22 705 1030 1972 OCT 23 731 1000 1972 OCT 24 688 1050 1972 OCT 25 749 1270 1972 OCT 26 818 1270 1972 OCT 27 954 1180 1972 OCT 28 985 1370 1972 OCT 29 923 1630 1972 OCT 30 908 1500 1972 OCT 31 1250 1350 1972 NOV 1 1410 1480 1972 NOV 2 1180 1660 1972 NOV 3 989 1560 1972 NOV 4 874 1370 1972 NOV 5 832 1230 1972 NOV 6 756 1150 1972 NOV 7 721 1080 1972 NOV 8 970 1180 1972 NOV 9 1460 1640 1972 NOV 10 2930 1910 1972 NOV 11 3510 2690 1972 NOV 12 3800 3070 1972 NOV 13 3980 3430 1972 NOV 14 4130 3720 1972 NOV 15 3170 3960 1972 NOV 16 2300 4030 1972 NOV 17 2850 3700 1972 NOV 18 3380 3360 1972 NOV 19 3930 3470 1972 NOV 20 4660 3750 1972 NOV 21 5230 4160 1972 NOV 22 5760 4610 1972 NOV 23 6300 5050 1972 NOV 24 6990 5480 1972 NOV 25 7760 5920 1972 NOV 26 8600 6670 1972 NOV 27 9120 7340 1972 NOV 28 8900 8050 1972 NOV 29 8820 8660 1972 NOV 30 8320 9210 1972 DEC 1 7860 9500 1972 DEC 2 7490 9380 1972 DEC 3 6660 8940 1972 DEC 4 6260 8390 1972 DEC 5 5680 7000 1972 DEC 6 4660 6400 1972 DEC 7 3630 5800 1972 DEC 8 3210 5400 1972 DEC 9 3060 4800 1972 DEC 10 3030 4200 1972 DEC 11 2910 3800 1972 DEC 12 2690 3600 1972 DEC 13 2640 3440 1972 DEC 14 2670 3250 1972 DEC 15 2890 3210 1972 DEC 16 4130 3340 1972 DEC 17 4770 3710 1972 DEC 18 5470 4200 1972 DEC 19 6320 4690 1972 DEC 20 8120 5260 1972 DEC 21 9680 5920 1972 DEC 22 10800 6900 1972 DEC 23 11800 8670 1972 DEC 24 12100 9740 1972 DEC 25 12100 11000 1972 DEC 26 11800 11900 1972 DEC 27 11400 12500 1972 DEC 28 10700 12600 1972 DEC 29 9970 12400 1972 DEC 30 9340 12000 1972 DEC 31 8410 11500 1973 JAN 1 6640 10900 1973 JAN 2 5180 9990 1973 JAN 3 4660 8750 1973 JAN 4 4530 7490 1973 JAN 5 4520 6430 1973 JAN 6 4460 5780 1973 JAN 7 4410 5420 1973 JAN 8 4500 5270 1973 JAN 9 4550 5180 1973 JAN 10 4410 5120 1973 JAN 11 4220 5080 1973 JAN 12 4220 4970 1973 JAN 13 3970 4860 1973 JAN 14 3670 4720 1973 JAN 15 3420 4540 1973 JAN 16 3410 4350 1973 JAN 17 3630 4220 1973 JAN 18 4020 4210 1973 JAN 19 4390 4390 1973 JAN 20 4650 4660 1973 JAN 21 4670 4920 1973 JAN 22 4580 5160 1973 JAN 23 4480 5470 1973 JAN 24 4400 5640 1973 JAN 25 4520 5530 1973 JAN 26 4580 5320 1973 JAN 27 4470 5220 1973 JAN 28 4070 5200 1973 JAN 29 4030 5180 1973 JAN 30 4300 5090 1973 JAN 31 4750 5000 1973 FEB 1 5070 5010 1973 FEB 2 5860 5380 1973 FEB 3 6760 6230 1973 FEB 4 7570 7100 1973 FEB 5 8370 7580 1973 FEB 6 10200 7980 1973 FEB 7 13400 8520 1973 FEB 8 16300 9500 1973 FEB 9 18500 11300 1973 FEB 10 18500 14800 1973 FEB 11 17500 16700 1973 FEB 12 16200 18400 1973 FEB 13 14800 18800 1973 FEB 14 13200 18600 1973 FEB 15 11700 18600 1973 FEB 16 10000 18200 1973 FEB 17 9040 16500 1973 FEB 18 8490 14200 1973 FEB 19 8050 12100 1973 FEB 20 7410 10500 1973 FEB 21 5980 9450 1973 FEB 22 4780 8560 1973 FEB 23 4130 7650 1973 FEB 24 3570 6600 1973 FEB 25 3220 5300 1973 FEB 26 2990 4760 1973 FEB 27 3040 4420 1973 FEB 28 3410 4400 1973 MAR 1 4050 4490 1973 MAR 2 4580 4610 1973 MAR 3 4920 4910 1973 MAR 4 5660 5280 1973 MAR 5 6030 5790 1973 MAR 6 6500 6410 1973 MAR 7 6980 6830 1973 MAR 8 7400 7130 1973 MAR 9 7740 7430 1973 MAR 10 7980 7760 1973 MAR 11 8150 8050 1973 MAR 12 8330 8300 1973 MAR 13 8260 8570 1973 MAR 14 8020 8770 1973 MAR 15 7730 8770 1973 MAR 16 7370 8590 1973 MAR 17 6860 8350 1973 MAR 18 5740 8020 1973 MAR 19 4640 7570 1973 MAR 20 4180 6970 1973 MAR 21 4060 6220 1973 MAR 22 4230 5470 1973 MAR 23 4150 5030 1973 MAR 24 4120 4850 1973 MAR 25 4260 4720 1973 MAR 26 4450 4730 1973 MAR 27 4570 4880 1973 MAR 28 4450 5050 1973 MAR 29 4050 5090 1973 MAR 30 3880 4990 1973 MAR 31 4100 4860 1973 APR 1 7470 6860 1973 APR 2 7660 9370 1973 APR 3 8340 9920 1973 APR 4 9130 9890 1973 APR 5 9870 9830 1973 APR 6 10800 10100 1973 APR 7 11600 10700 1973 APR 8 12600 11900 1973 APR 9 12800 13400 1973 APR 10 12900 14100 1973 APR 11 12700 14300 1973 APR 12 12300 14200 1973 APR 13 11600 14100 1973 APR 14 10900 13800 1973 APR 15 10000 13300 1973 APR 16 9120 12700 1973 APR 17 7360 11800 1973 APR 18 5170 10800 1973 APR 19 3550 9450 1973 APR 20 3000 7330 1973 APR 21 2660 5560 1973 APR 22 2440 4340 1973 APR 23 2270 3550 1973 APR 24 2100 3220 1973 APR 25 1970 2910 1973 APR 26 1880 2700 1973 APR 27 2060 2750 1973 APR 28 2870 2900 1973 APR 29 3670 3080 1973 APR 30 4230 3400 1973 MAY 1 4510 3860 1973 MAY 2 4910 4260 1973 MAY 3 5280 4590 1973 MAY 4 5670 5090 1973 MAY 5 5640 5640 1973 MAY 6 4550 5990 1973 MAY 7 3870 6000 1973 MAY 8 2990 5590 1973 MAY 9 2480 4800 1973 MAY 10 2250 3880 1973 MAY 11 2350 3380 1973 MAY 12 2550 2970 1973 MAY 13 2240 2920 1973 MAY 14 1730 2830 1973 MAY 15 1450 2450 1973 MAY 16 1300 2050 1973 MAY 17 1200 1810 1973 MAY 18 1140 1670 1973 MAY 19 1100 1580 1973 MAY 20 1110 1510 1973 MAY 21 1100 1520 1973 MAY 22 1260 1520 1973 MAY 23 1690 1570 1973 MAY 24 1610 1910 1973 MAY 25 1500 2060 1973 MAY 26 1600 1990 1973 MAY 27 1810 1970 1973 MAY 28 1630 2090 1973 MAY 29 1440 2050 1973 MAY 30 2250 1880 1973 MAY 31 3080 2270 1973 JUN 1 3520 2740 1973 JUN 2 3420 3150 1973 JUN 3 2250 3400 1973 JUN 4 1600 3230 1973 JUN 5 1210 2240 1973 JUN 6 1030 1740 1973 JUN 7 946 1450 1973 JUN 8 958 1330 1973 JUN 9 1030 1280 1973 JUN 10 882 1340 1973 JUN 11 842 1250 1973 JUN 12 826 1180 1973 JUN 13 820 1210 1973 JUN 14 840 1260 1973 JUN 15 790 1130 1973 JUN 16 700 1070 1973 JUN 17 850 1090 1973 JUN 18 1100 1240 1973 JUN 19 1300 1430 1973 JUN 20 1500 1720 1973 JUN 21 1900 2380 1973 JUN 22 2100 2580 1973 JUN 23 2350 2840 1973 JUN 24 3400 2970 1973 JUN 25 4040 3140 1973 JUN 26 4780 3540 1973 JUN 27 5610 3970 1973 JUN 28 6180 4490 1973 JUN 29 6390 5080 1973 JUN 30 5120 5610 1973 JUL 1 3930 5990 1973 JUL 2 3810 5880 1973 JUL 3 4290 5210 1973 JUL 4 5860 4840 1973 JUL 5 7370 4740 1973 JUL 6 8060 5350 1973 JUL 7 8800 6700 1973 JUL 8 9000 7270 1973 JUL 9 7500 7930 1973 JUL 10 6400 8250 1973 JUL 11 4500 8090 1973 JUL 12 3000 6310 1973 JUL 13 2200 3820 1973 JUL 14 1700 2230 1973 JUL 15 1600 1860 1973 JUL 16 1500 1770 1973 JUL 17 1800 1700 1973 JUL 18 2200 1940 1973 JUL 19 2500 2360 1973 JUL 20 2410 2550 1973 JUL 21 2010 2580 1973 JUL 22 1580 2320 1973 JUL 23 1260 1900 1973 JUL 24 1010 1530 1973 JUL 25 864 1260 1973 JUL 26 798 1080 1973 JUL 27 780 983 1973 JUL 28 741 1070 1973 JUL 29 737 1170 1973 JUL 30 830 1060 1973 JUL 31 850 1010 1973 AUG 1 918 1060 1973 AUG 2 875 1150 1973 AUG 3 800 1130 1973 AUG 4 1520 1050 1973 AUG 5 1690 1460 1973 AUG 6 2000 1740 1973 AUG 7 2200 2050 1973 AUG 8 2500 2380 1973 AUG 9 1800 2480 1973 AUG 10 1200 1790 1973 AUG 11 1000 1300 1973 AUG 12 900 1110 1973 AUG 13 800 1030 1973 AUG 14 850 931 1973 AUG 15 900 979 1973 AUG 16 880 1120 1973 AUG 17 920 1110 1973 AUG 18 860 1140 1973 AUG 19 800 987 1973 AUG 20 900 899 1973 AUG 21 1800 987 1973 AUG 22 2300 1730 1973 AUG 23 2490 2260 1973 AUG 24 1860 2560 1973 AUG 25 1630 2400 1973 AUG 26 1270 1940 1973 AUG 27 925 1570 1973 AUG 28 761 1220 1973 AUG 29 674 1000 1973 AUG 30 614 883 1973 AUG 31 566 807 1973 SEP 1 537 753 1973 SEP 2 501 711 1973 SEP 3 470 676 1973 SEP 4 451 648 1973 SEP 5 446 624 1973 SEP 6 416 613 1973 SEP 7 394 592 1973 SEP 8 382 564 1973 SEP 9 389 557 1973 SEP 10 398 634 1973 SEP 11 388 624 1973 SEP 12 383 582 1973 SEP 13 375 564 1973 SEP 14 404 575 1973 SEP 15 367 627 1973 SEP 16 381 599 1973 SEP 17 558 543 1973 SEP 18 542 655 1973 SEP 19 486 704 1973 SEP 20 432 673 1973 SEP 21 392 617 1973 SEP 22 366 568 1973 SEP 23 344 536 1973 SEP 24 323 508 1973 SEP 25 301 484 1973 SEP 26 294 466 1973 SEP 27 295 459 1973 SEP 28 280 445 1973 SEP 29 267 430 1973 SEP 30 256 415 1973 OCT 1 275 409 1973 OCT 2 279 608 1973 OCT 3 273 778 1973 OCT 4 263 641 1973 OCT 5 271 529 1973 OCT 6 275 480 1973 OCT 7 268 456 1973 OCT 8 257 442 1973 OCT 9 242 439 1973 OCT 10 234 424 1973 OCT 11 226 403 1973 OCT 12 226 391 1973 OCT 13 220 385 1973 OCT 14 221 382 1973 OCT 15 216 379 1973 OCT 16 208 370 1973 OCT 17 200 364 1973 OCT 18 205 352 1973 OCT 19 203 343 1973 OCT 20 193 340 1973 OCT 21 188 331 1973 OCT 22 189 322 1973 OCT 23 188 313 1973 OCT 24 185 313 1973 OCT 25 185 301 1973 OCT 26 190 295 1973 OCT 27 188 292 1973 OCT 28 183 292 1973 OCT 29 192 313 1973 OCT 30 173 355 1973 OCT 31 177 328 1973 NOV 1 202 304 1973 NOV 2 237 304 1973 NOV 3 222 325 1973 NOV 4 206 327 1973 NOV 5 201 322 1973 NOV 6 183 310 1973 NOV 7 177 298 1973 NOV 8 180 286 1973 NOV 9 182 283 1973 NOV 10 181 280 1973 NOV 11 174 280 1973 NOV 12 174 277 1973 NOV 13 167 271 1973 NOV 14 167 277 1973 NOV 15 168 274 1973 NOV 16 176 274 1973 NOV 17 177 271 1973 NOV 18 177 277 1973 NOV 19 196 286 1973 NOV 20 202 301 1973 NOV 21 207 316 1973 NOV 22 234 316 1973 NOV 23 236 331 1973 NOV 24 250 334 1973 NOV 25 286 340 1973 NOV 26 283 361 1973 NOV 27 271 385 1973 NOV 28 262 385 1973 NOV 29 262 379 1973 NOV 30 253 376 1973 DEC 1 256 364 1973 DEC 2 271 355 1973 DEC 3 271 358 1973 DEC 4 265 364 1973 DEC 5 262 367 1973 DEC 6 256 364 1973 DEC 7 253 358 1973 DEC 8 368 442 1973 DEC 9 597 863 1973 DEC 10 549 1120 1973 DEC 11 1030 915 1973 DEC 12 1130 1060 1973 DEC 13 1020 1260 1973 DEC 14 888 1320 1973 DEC 15 702 1260 1973 DEC 16 778 1110 1973 DEC 17 981 1390 1973 DEC 18 1380 1630 1973 DEC 19 1980 1650 1973 DEC 20 2110 2080 1973 DEC 21 2240 2500 1973 DEC 22 2130 2850 1973 DEC 23 2800 2840 1973 DEC 24 3230 2920 1973 DEC 25 3450 3150 1973 DEC 26 3600 3390 1973 DEC 27 3750 3590 1973 DEC 28 3750 3750 1973 DEC 29 2320 3930 1973 DEC 30 1520 3760 1973 DEC 31 1330 2580 1974 JAN 1 1450 2070 1974 JAN 2 2490 2130 1974 JAN 3 3250 2670 1974 JAN 4 3710 3060 1974 JAN 5 3980 3490 1974 JAN 6 4190 3830 1974 JAN 7 4280 4100 1974 JAN 8 4060 4320 1974 JAN 9 3480 4420 1974 JAN 10 2820 4320 1974 JAN 11 2500 3920 1974 JAN 12 2730 3360 1974 JAN 13 3060 3090 1974 JAN 14 3050 3140 1974 JAN 15 2440 3250 1974 JAN 16 1880 3090 1974 JAN 17 1620 2620 1974 JAN 18 1480 2160 1974 JAN 19 1420 1920 1974 JAN 20 1350 1820 1974 JAN 21 1320 1810 1974 JAN 22 1300 1880 1974 JAN 23 1360 1840 1974 JAN 24 1930 1780 1974 JAN 25 2740 2060 1974 JAN 26 3190 2680 1974 JAN 27 3330 3080 1974 JAN 28 3520 3350 1974 JAN 29 3700 3620 1974 JAN 30 3870 3920 1974 JAN 31 4100 4130 1974 FEB 1 4370 4270 1974 FEB 2 4760 4430 1974 FEB 3 5550 4710 1974 FEB 4 6360 5160 1974 FEB 5 6840 5730 1974 FEB 6 7050 6250 1974 FEB 7 7010 6700 1974 FEB 8 6580 7040 1974 FEB 9 5910 7180 1974 FEB 10 5220 7080 1974 FEB 11 4830 6710 1974 FEB 12 4670 6180 1974 FEB 13 4600 5660 1974 FEB 14 4390 5300 1974 FEB 15 3310 5130 1974 FEB 16 2950 4910 1974 FEB 17 3510 4500 1974 FEB 18 4390 4280 1974 FEB 19 5210 4390 1974 FEB 20 5790 4800 1974 FEB 21 6070 5290 1974 FEB 22 6200 5740 1974 FEB 23 6240 6120 1974 FEB 24 6130 6340 1974 FEB 25 5970 6440 1974 FEB 26 5840 6410 1974 FEB 27 5630 6320 1974 FEB 28 5070 6200 1974 MAR 1 3670 6040 1974 MAR 2 2620 5670 1974 MAR 3 2360 4540 1974 MAR 4 2150 3560 1974 MAR 5 2010 3020 1974 MAR 6 1900 2620 1974 MAR 7 1800 2400 1974 MAR 8 1710 2240 1974 MAR 9 1610 2120 1974 MAR 10 1540 2010 1974 MAR 11 1480 1920 1974 MAR 12 1490 1930 1974 MAR 13 1550 1990 1974 MAR 14 1720 1990 1974 MAR 15 1800 2050 1974 MAR 16 1930 2260 1974 MAR 17 2340 2780 1974 MAR 18 2660 3100 1974 MAR 19 3180 3200 1974 MAR 20 3620 3580 1974 MAR 21 3840 3950 1974 MAR 22 3930 4190 1974 MAR 23 3980 4380 1974 MAR 24 4090 4600 1974 MAR 25 4200 4790 1974 MAR 26 4460 4950 1974 MAR 27 4570 5130 1974 MAR 28 4660 5230 1974 MAR 29 4580 5290 1974 MAR 30 5340 5560 1974 MAR 31 5920 5980 1974 APR 1 6520 6350 1974 APR 2 7040 6570 1974 APR 3 7230 6810 1974 APR 4 7070 7090 1974 APR 5 7110 7650 1974 APR 6 7140 7940 1974 APR 7 6730 7950 1974 APR 8 5900 7810 1974 APR 9 5080 7610 1974 APR 10 4160 7120 1974 APR 11 3370 6360 1974 APR 12 3170 5420 1974 APR 13 3220 4590 1974 APR 14 4430 4350 1974 APR 15 4790 4500 1974 APR 16 4850 5000 1974 APR 17 4370 5260 1974 APR 18 3660 5330 1974 APR 19 3040 5120 1974 APR 20 2590 4600 1974 APR 21 2250 3910 1974 APR 22 1930 3290 1974 APR 23 1690 2800 1974 APR 24 1500 2390 1974 APR 25 1360 2020 1974 APR 26 1260 1870 1974 APR 27 1180 1710 1974 APR 28 1080 1600 1974 APR 29 1000 1490 1974 APR 30 933 1400 1974 MAY 1 862 1310 1974 MAY 2 819 1230 1974 MAY 3 840 1190 1974 MAY 4 1050 1170 1974 MAY 5 1410 1280 1974 MAY 6 1600 1660 1974 MAY 7 1650 1900 1974 MAY 8 1930 1930 1974 MAY 9 2380 2080 1974 MAY 10 2290 2400 1974 MAY 11 1970 2550 1974 MAY 12 2700 2380 1974 MAY 13 3270 2540 1974 MAY 14 3520 2930 1974 MAY 15 3730 3240 1974 MAY 16 3960 3470 1974 MAY 17 4160 3660 1974 MAY 18 4340 3830 1974 MAY 19 3820 4020 1974 MAY 20 2280 4180 1974 MAY 21 1890 3670 1974 MAY 22 1920 2670 1974 MAY 23 2080 2330 1974 MAY 24 2110 2340 1974 MAY 25 2130 2440 1974 MAY 26 2300 2400 1974 MAY 27 2380 2760 1974 MAY 28 2930 2890 1974 MAY 29 3670 2970 1974 MAY 30 4110 3290 1974 MAY 31 4350 3690 1974 JUN 1 3890 4070 1974 JUN 2 2490 4470 1974 JUN 3 2050 4290 1974 JUN 4 2330 3160 1974 JUN 5 2390 2910 1974 JUN 6 2160 2910 1974 JUN 7 1760 2690 1974 JUN 8 1370 2270 1974 JUN 9 1120 1820 1974 JUN 10 970 1510 1974 JUN 11 874 1320 1974 JUN 12 773 1190 1974 JUN 13 724 1090 1974 JUN 14 676 1020 1974 JUN 15 629 973 1974 JUN 16 609 923 1974 JUN 17 1490 929 1974 JUN 18 2930 1610 1974 JUN 19 3260 2550 1974 JUN 20 2810 2890 1974 JUN 21 2210 2960 1974 JUN 22 1660 2670 1974 JUN 23 1150 1780 1974 JUN 24 1060 1390 1974 JUN 25 1180 1390 1974 JUN 26 1010 1340 1974 JUN 27 919 1310 1974 JUN 28 821 1250 1974 JUN 29 861 1200 1974 JUN 30 851 1120 1974 JUL 1 862 1100 1974 JUL 2 859 1050 1974 JUL 3 681 1090 1974 JUL 4 572 961 1974 JUL 5 514 849 1974 JUL 6 486 822 1974 JUL 7 550 792 1974 JUL 8 722 857 1974 JUL 9 836 994 1974 JUL 10 828 1060 1974 JUL 11 709 1080 1974 JUL 12 600 982 1974 JUL 13 539 866 1974 JUL 14 484 789 1974 JUL 15 420 737 1974 JUL 16 371 744 1974 JUL 17 339 706 1974 JUL 18 316 614 1974 JUL 19 300 572 1974 JUL 20 299 544 1974 JUL 21 484 597 1974 JUL 22 497 785 1974 JUL 23 686 784 1974 JUL 24 726 874 1974 JUL 25 762 975 1974 JUL 26 738 1040 1974 JUL 27 912 1070 1974 JUL 28 2040 1210 1974 JUL 29 2580 2080 1974 JUL 30 2640 2480 1974 JUL 31 1910 2700 1974 AUG 1 1600 2470 1974 AUG 2 2360 2210 1974 AUG 3 2850 3000 1974 AUG 4 2880 3430 1974 AUG 5 4370 3940 1974 AUG 6 7350 4440 1974 AUG 7 8240 5690 1974 AUG 8 9100 7020 1974 AUG 9 9660 7940 1974 AUG 10 10000 9170 1974 AUG 11 10100 10000 1974 AUG 12 10000 10400 1974 AUG 13 9800 10600 1974 AUG 14 8040 10700 1974 AUG 15 5960 10600 1974 AUG 16 4670 9600 1974 AUG 17 4200 9000 1974 AUG 18 3700 8000 1974 AUG 19 3210 6600 1974 AUG 20 2770 5440 1974 AUG 21 2450 4700 1974 AUG 22 2590 4010 1974 AUG 23 2380 3620 1974 AUG 24 1950 3510 1974 AUG 25 1910 3400 1974 AUG 26 1920 3290 1974 AUG 27 1820 3000 1974 AUG 28 1590 2760 1974 AUG 29 1530 2320 1974 AUG 30 1500 2010 1974 AUG 31 1310 1870 1974 SEP 1 1130 1770 1974 SEP 2 1200 1830 1974 SEP 3 978 1810 1974 SEP 4 1120 1580 1974 SEP 5 2290 1550 1974 SEP 6 3050 2620 1974 SEP 7 3840 3410 1974 SEP 8 4480 4010 1974 SEP 9 4930 4440 1974 SEP 10 5160 4800 1974 SEP 11 5220 5100 1974 SEP 12 5200 5310 1974 SEP 13 5180 5430 1974 SEP 14 5170 5480 1974 SEP 15 5280 5460 1974 SEP 16 5370 5440 1974 SEP 17 3490 5470 1974 SEP 18 1850 5440 1974 SEP 19 1570 4210 1974 SEP 20 1540 2560 1974 SEP 21 1330 2050 1974 SEP 22 1200 1700 1974 SEP 23 1260 1520 1974 SEP 24 1470 1490 1974 SEP 25 1500 1600 1974 SEP 26 1300 1670 1974 SEP 27 1150 1560 1974 SEP 28 1040 1380 1974 SEP 29 952 1200 1974 SEP 30 940 1100 1974 OCT 1 944 1150 1974 OCT 2 903 1130 1974 OCT 3 811 1080 1974 OCT 4 737 1020 1974 OCT 5 687 948 1974 OCT 6 640 892 1974 OCT 7 610 845 1974 OCT 8 584 810 1974 OCT 9 567 783 1974 OCT 10 542 760 1974 OCT 11 524 734 1974 OCT 12 515 714 1974 OCT 13 499 698 1974 OCT 14 477 689 1974 OCT 15 459 668 1974 OCT 16 493 660 1974 OCT 17 615 725 1974 OCT 18 701 806 1974 OCT 19 915 883 1974 OCT 20 980 1090 1974 OCT 21 1110 1160 1974 OCT 22 1190 1220 1974 OCT 23 1100 1310 1974 OCT 24 988 1270 1974 OCT 25 867 1180 1974 OCT 26 768 1080 1974 OCT 27 704 981 1974 OCT 28 667 907 1974 OCT 29 637 863 1974 OCT 30 617 833 1974 OCT 31 594 795 1974 NOV 1 577 769 1974 NOV 2 570 743 1974 NOV 3 567 729 1974 NOV 4 560 718 1974 NOV 5 554 711 1974 NOV 6 551 701 1974 NOV 7 538 694 1974 NOV 8 532 680 1974 NOV 9 535 673 1974 NOV 10 533 666 1974 NOV 11 518 666 1974 NOV 12 514 666 1974 NOV 13 510 655 1974 NOV 14 509 645 1974 NOV 15 518 634 1974 NOV 16 525 634 1974 NOV 17 521 638 1974 NOV 18 603 669 1974 NOV 19 705 825 1974 NOV 20 994 1060 1974 NOV 21 1350 1380 1974 NOV 22 1610 1600 1974 NOV 23 1790 1740 1974 NOV 24 1690 1890 1974 NOV 25 1490 1880 1974 NOV 26 1310 1750 1974 NOV 27 1130 1580 1974 NOV 28 1040 1410 1974 NOV 29 975 1280 1974 NOV 30 918 1200 1974 DEC 1 1150 1540 1974 DEC 2 1440 1960 1974 DEC 3 2470 2100 1974 DEC 4 2820 2450 1974 DEC 5 2730 2830 1974 DEC 6 2220 2980 1974 DEC 7 1770 2800 1974 DEC 8 1650 2400 1974 DEC 9 1830 2190 1974 DEC 10 2690 2200 1974 DEC 11 3260 2520 1974 DEC 12 3580 2950 1974 DEC 13 3420 3310 1974 DEC 14 2460 3510 1974 DEC 15 2040 3350 1974 DEC 16 1810 2800 1974 DEC 17 1750 2300 1974 DEC 18 1910 2080 1974 DEC 19 2100 2100 1974 DEC 20 2110 2270 1974 DEC 21 2090 2480 1974 DEC 22 2110 2570 1974 DEC 23 2240 2550 1974 DEC 24 2240 2560 1974 DEC 25 2130 2570 1974 DEC 26 1970 2490 1974 DEC 27 1810 2320 1974 DEC 28 1960 2260 1974 DEC 29 2270 2440 1974 DEC 30 2920 2650 1974 DEC 31 3580 2930 ; goptions reset=all; proc gplot data=neuse; plot lgold*date lkins*date/overlay legend; symbol1 v=none i=join c=red w=1; symbol2 v=none i=join c=blue w=1; proc arima data=neuse; i var =lgold(1) outcov=gold; by group season; proc arima data=neuse; i var=lgold(1) outcov=outno; proc print; run; data gold; set gold; nlag=-1*lag; lags = lag+.1*(season-1); *** separate seasons ***; seas=season+.0001*ranuni(1827655) ; ** jitter **; lb=.; ub=.; if group=1 then do; lb=-1.96*stderr; ub=1.96*stderr; lbmn=-1.96*stderr/4; ubmn=1.96*stderr/4; end; color="green"; if season=1 then color="red"; if season=2 then color="cyan"; if season=3 then color="black"; proc g3d data=gold; scatter group*nlag=corr/ rotate=80 tilt=80 color=color shape="cylinder" noneedle size=.75; plot group*nlag=corr/rotate=80 tilt=80; plot lag*group=corr/rotate=80 tilt=80; scatter seas*nlag=corr/ rotate=80 tilt=80 color=color shape="balloon" size=.75; run; proc gplot data=gold; plot corr*lags=season; symbol1 v=dot i=none c=red; symbol2 v=dot i=none c=cyan; symbol3 v=dot i=none c=black; symbol4 v=dot i=none c=green; where lags < 10; run; goptions reset=all; proc gplot data=gold; plot corr*lag=2 (ub lb)*lag = 1/overlay vref=0; symbol1 v=none i=join c=red width=3; symbol2 v=none i=boxt; proc gplot data=gold; plot invcorr*lag=2/overlay vref=-.22 -.055 0 .055 .22; proc gplot data=gold; plot partcorr*lag=2/overlay vref=-.22 -.055 0 .055 .22; data gold1; set gold; if lag in (1 9 10 11 22); proc sort data=gold1; by lag; proc glm data=gold1; class season lag; model corr = season/solution; means season/tukey; by lag; run; ************ crosscorrelations ****************; proc arima data=neuse; i var=lgold(1) noprint; e p=1 q=1 ml noprint; i var=lkins(1) crosscor=lgold(1) noprint outcov=cross; by group; data cross; set cross; where crossvar="LGOLD"; season = 1 + mod(group-1,4); lags = lag + .1*(season-2.5); proc print; where corr>.8; proc gplot data=cross; plot corr*lags=season; symbol1 v=dot i=none c=red; symbol2 v=dot i=none c=cyan; symbol3 v=dot i=none c=black; symbol4 v=dot i=none c=green; where abs(lag) < 8; proc gplot data=cross; plot corr*lags=season/ vref=0; symbol1 v=none i=boxt c=red; symbol2 v=none i=boxt c=cyan; symbol3 v=none i=boxt c=black; symbol4 v=none i=boxt c=green; proc gplot data=cross; plot corr*lags=season/ vref=0; where abs(lag) < 8; run; proc greplay tc=tempcat nofs igout=work.gseg; tdef three des="three panel stack" 1/ llx= 5 lly=67 ulx= 5 uly=98 urx=95 ury=98 lrx=95 lry=67 2/ llx= 5 lly=34 ulx= 5 uly=65 urx=95 ury=65 lrx=95 lry=34 3/ llx= 5 lly= 2 ulx= 5 uly=32 urx=95 ury=32 lrx=95 lry= 2 ; template=three ; goptions lfactor=3 targetdevice=pscolor; treplay 1:gplot 2:g3d3 3:gplot1; treplay 1:gplot2 2:gplot3 3:gplot4; treplay 1:gplot5 2:gplot6 3:gplot7; run;