data counts;
  input visits sex race marital private education
        person strata psu SamplingWeight @@;
  datalines;
5 1 1 2 1 5 71511 1 1 1002.59 1 2 1 4 2 3 307568 1 1 1002.59
2 1 1 4 4 3 457473 1 1 1002.59 9 1 1 3 1 5 849963 1 1 1002.59
3 2 1 3 2 5 892466 1 1 1002.59 0 2 1 2 3 3 249075 1 2 1002.59
3 1 1 2 4 1 835408 1 2 1002.59 1 2 1 4 2 4 159262 1 3 1002.59
2 2 1 3 4 3 539031 1 3 1002.59 0 2 1 2 2 5 953026 1 3 1002.59
4 1 1 3 4 2 31739 1 4 1002.59 1 2 1 3 2 4 427721 1 4 1002.59
3 2 2 4 4 5 591883 1 4 1002.59 1 2 1 1 4 3 693967 1 4 1002.59
2 2 1 3 2 5 935953 1 4 1002.59 3 1 2 1 2 4 388876 1 5 1002.59
1 1 1 3 2 5 478681 1 5 1002.59 5 2 1 1 1 5 747715 1 5 1002.59
2 2 1 4 3 5 791387 1 5 1002.59 1 1 1 4 3 3 803485 1 5 1002.59
3 1 1 3 3 3 61866 1 6 1002.59 0 2 2 4 2 3 319632 1 6 1002.59
2 2 1 1 2 3 329846 1 6 1002.59 4 2 1 4 1 2 465451 1 6 1002.59
3 2 2 2 2 5 719335 1 6 1002.59 3 2 1 3 2 2 835100 1 6 1002.59
1 1 1 4 4 5 890095 1 6 1002.59 3 1 2 2 2 5 949634 1 6 1002.59
0 2 1 4 2 1 33103 1 7 1002.59 8 1 3 2 1 2 120597 1 7 1002.59
1 1 1 4 3 1 245813 1 7 1002.59 1 2 1 3 2 3 390648 1 7 1002.59
3 2 1 4 2 2 512783 1 7 1002.59 1 1 1 3 4 2 735952 1 7 1002.59
3 2 1 3 3 4 762927 1 7 1002.59 1 1 1 4 3 1 904492 1 7 1002.59
3 1 1 3 2 3 50782 1 8 1002.59 1 1 3 3 2 4 432975 1 8 1002.59
2 1 1 4 4 4 482739 1 8 1002.59 1 1 1 4 2 3 498432 1 8 1002.59
1 2 1 3 4 5 534998 1 8 1002.59 1 2 1 4 2 5 723920 1 8 1002.59
3 1 1 1 2 2 414253 1 9 1002.59 4 2 3 4 3 3 516264 1 9 1002.59
4 1 1 1 2 1 547482 1 9 1002.59 0 1 1 4 2 4 566081 1 9 1002.59
2 2 1 2 4 4 589708 1 9 1002.59 3 1 1 1 2 3 608070 1 9 1002.59
3 2 3 2 2 5 692070 1 9 1002.59 1 1 1 1 2 3 906853 1 9 1002.59
2 2 1 4 2 3 370343 1 10 1002.59 2 2 1 3 4 5 69243 1 11 1002.59
0 2 2 1 2 3 109404 1 11 1002.59 4 2 1 1 2 5 223404 1 11 1002.59
2 1 1 3 2 2 296117 1 11 1002.59 2 2 1 3 2 2 444399 1 11 1002.59
3 1 2 4 4 5 513426 1 12 1002.59 13 1 1 3 1 4 690046 1 12 1002.59
0 1 1 3 2 4 823188 1 12 1002.59 1 2 1 3 4 4 128338 1 13 1002.59
4 2 1 3 2 2 241355 1 13 1002.59 3 2 1 3 2 4 482082 1 13 1002.59
2 2 1 3 4 2 730562 1 13 1002.59 1 2 1 2 4 5 969557 1 13 1002.59
2 2 1 2 4 4 134594 1 14 1002.59 3 1 1 2 2 4 171804 1 14 1002.59
1 1 1 2 2 5 527159 1 14 1002.59 3 1 3 2 2 1 539564 1 14 1002.59
1 2 1 3 4 4 159322 1 15 1002.59 3 2 1 3 2 2 299917 1 15 1002.59
2 1 1 1 2 3 271341 1 16 1002.59 2 2 2 2 2 2 490536 1 16 1002.59
2 1 1 3 2 2 625312 1 16 1002.59 3 1 1 3 2 5 18187 1 17 1002.59
3 1 1 2 1 2 36694 1 17 1002.59 3 2 1 4 2 4 983600 1 17 1002.59
0 1 1 4 2 5 156220 1 18 1002.59 1 1 1 4 2 5 362234 1 18 1002.59
1 2 1 4 3 4 688707 1 18 1002.59 1 2 1 4 2 1 768987 1 18 1002.59
7 2 1 3 1 2 874484 1 18 1002.59 1 2 2 4 2 4 588677 1 19 1002.59
6 1 1 4 1 3 583837 1 20 1002.59 1 2 3 1 2 4 651486 1 20 1002.59
0 1 2 3 2 4 705394 1 20 1002.59 1 2 1 1 2 4 824793 1 20 1002.59
5 1 1 3 2 2 860034 1 20 1002.59 2 2 1 3 4 2 862689 1 20 1002.59
1 2 1 2 4 2 862960 1 20 1002.59 1 1 1 3 4 4 885992 1 20 1002.59
2 1 2 4 3 3 162169 1 21 1002.59 0 2 1 4 2 3 491248 1 21 1002.59
1 1 1 2 2 3 108462 1 22 1002.59 4 2 2 2 2 2 325843 1 22 1002.59
3 2 2 2 2 4 442223 1 22 1002.59 1 1 1 1 3 5 461280 1 22 1002.59
0 1 2 4 2 1 693117 1 22 1002.59 1 1 2 2 2 5 180818 1 23 1002.59
2 2 1 3 2 5 188138 1 23 1002.59 2 1 1 3 2 4 283611 1 23 1002.59
2 1 1 2 2 4 328911 1 23 1002.59 3 1 3 1 2 3 568325 1 23 1002.59
1 2 1 1 3 4 582040 1 23 1002.59 0 2 1 4 4 2 754435 1 23 1002.59
2 1 1 4 2 3 174522 1 24 1002.59 2 2 2 3 2 5 502693 1 24 1002.59
1 1 1 3 3 3 608199 1 24 1002.59 1 1 1 2 2 5 897243 1 24 1002.59
3 1 1 3 2 4 988408 1 24 1002.59 3 1 2 1 2 4 157613 1 25 1002.59
0 1 1 2 2 1 240886 1 25 1002.59 5 1 1 2 1 5 600622 1 25 1002.59
2 2 1 4 2 3 767713 1 25 1002.59 1 2 1 3 3 1 849142 1 25 1002.59
1 1 1 4 3 2 921937 1 25 1002.59 2 1 3 2 2 4 966265 1 25 1002.59
3 1 1 4 4 5 100721 1 26 1002.59 1 2 1 1 2 2 171564 1 26 1002.59
9 2 1 3 1 4 537184 1 26 1002.59 2 1 1 4 2 1 679042 1 26 1002.59
1 1 2 1 2 4 718095 1 26 1002.59 2 1 1 4 4 4 748859 1 26 1002.59
1 1 2 3 2 1 924146 1 26 1002.59 2 2 1 4 2 2 946043 1 26 1002.59
2 2 1 3 2 5 39447 1 27 1002.59 10 2 1 3 1 4 275718 1 27 1002.59
1 2 2 4 2 3 280462 1 27 1002.59 0 1 1 4 3 2 400907 1 27 1002.59
0 1 3 1 2 3 489385 1 27 1002.59 0 2 1 2 2 4 681587 1 27 1002.59
3 1 3 4 4 3 699228 1 27 1002.59 2 2 3 4 2 2 153810 1 28 1002.59
0 2 3 4 2 5 336790 1 28 1002.59 1 1 1 4 4 5 454729 1 28 1002.59
1 1 1 2 2 3 762097 1 28 1002.59 3 2 1 2 2 2 179258 1 29 1002.59
1 1 1 3 4 5 249782 1 29 1002.59 2 1 1 1 4 4 409699 1 29 1002.59
1 2 1 1 2 4 460525 1 29 1002.59 4 1 2 4 4 3 623283 1 29 1002.59
4 2 1 2 1 4 751067 1 29 1002.59 2 1 1 1 2 4 822375 1 29 1002.59
2 1 1 2 2 2 985415 1 29 1002.59 7 1 1 2 2 5 200985 1 30 1002.59
0 1 1 1 4 5 418663 1 30 1002.59 1 1 1 3 3 4 586897 1 30 1002.59
2 2 1 4 2 3 833411 1 30 1002.59 1 2 1 1 2 5 30179 1 31 1002.59
0 2 2 3 3 4 168674 1 31 1002.59 2 2 1 4 2 5 223882 1 31 1002.59
1 2 1 1 2 2 548623 1 31 1002.59 1 2 1 1 2 1 570012 1 31 1002.59
2 1 3 4 2 2 649257 1 31 1002.59 1 2 1 4 2 3 680707 1 31 1002.59
2 1 1 3 2 5 794662 1 31 1002.59 6 2 1 4 1 2 935229 1 31 1002.59
3 1 1 4 3 5 944988 1 31 1002.59 3 2 1 3 2 5 326840 1 32 1002.59
0 1 1 3 3 3 344281 1 32 1002.59 2 1 1 3 2 1 575539 1 32 1002.59
3 1 3 3 2 4 648984 1 32 1002.59 1 1 1 3 3 1 697005 1 32 1002.59
2 2 1 4 2 5 720859 1 32 1002.59 0 1 1 1 2 2 803909 1 32 1002.59
1 2 1 4 4 2 433434 1 33 1002.59 1 2 1 3 3 4 83675 1 34 1002.59
3 1 2 4 2 5 514132 1 34 1002.59 8 1 1 2 1 5 763024 1 34 1002.59
2 1 1 3 2 3 125475 1 35 1002.59 0 1 3 3 2 5 176042 1 35 1002.59
4 1 1 4 2 4 361911 1 35 1002.59 3 2 2 2 2 4 547048 1 35 1002.59
0 1 1 3 2 5 759938 1 35 1002.59 1 2 1 3 2 1 763843 1 35 1002.59
3 1 1 4 2 3 918570 1 35 1002.59 3 1 1 1 3 2 957479 1 35 1002.59
2 1 1 3 4 4 52346 1 36 1002.59 1 1 1 1 2 2 69927 1 36 1002.59
1 2 1 3 2 3 111280 1 36 1002.59 2 1 1 2 2 3 518245 1 36 1002.59
2 2 3 1 2 3 23130 1 37 1002.59 3 1 1 1 4 2 85888 1 37 1002.59
0 1 1 4 2 3 90369 1 37 1002.59 1 2 1 3 2 5 569540 1 37 1002.59
2 2 1 1 2 5 469912 1 38 1002.59 2 1 1 3 2 2 536549 1 38 1002.59
2 1 1 3 2 5 711502 1 38 1002.59 6 1 1 3 4 1 180534 1 39 1002.59
1 2 1 4 2 1 348253 1 39 1002.59 2 2 1 4 2 4 479889 1 39 1002.59
2 2 1 4 3 3 592396 1 39 1002.59 5 2 1 2 2 2 619590 1 39 1002.59
4 1 1 3 2 2 785757 1 39 1002.59 3 1 1 3 2 5 265625 1 40 1002.59
6 1 1 3 4 2 382940 1 40 1002.59 0 1 1 2 3 2 383985 1 40 1002.59
1 2 1 3 2 2 600821 1 40 1002.59 8 2 2 4 1 2 723163 1 40 1002.59
1 2 1 4 4 2 988048 1 40 1002.59 3 1 1 4 2 4 166233 2 1 1000.99
0 1 3 3 2 4 222498 2 1 1000.99 2 1 1 2 2 4 366172 2 1 1000.99
1 2 1 3 4 4 562093 2 1 1000.99 2 1 1 4 4 3 915546 2 1 1000.99
10 2 1 4 1 1 955668 2 1 1000.99 4 2 1 1 2 2 316436 2 2 1000.99
3 1 1 1 2 3 998721 2 2 1000.99 1 1 1 2 4 2 11495 2 3 1000.99
6 2 1 4 1 4 449562 2 3 1000.99 4 2 1 3 2 3 722119 2 3 1000.99
2 1 1 4 2 2 346922 2 4 1000.99 2 1 1 1 2 3 453978 2 4 1000.99
3 1 1 4 2 2 717901 2 4 1000.99 1 1 3 4 2 1 625 2 5 1000.99
2 1 1 2 3 3 188183 2 5 1000.99 1 1 3 4 2 3 213840 2 5 1000.99
1 2 1 4 2 1 219077 2 5 1000.99 1 1 1 3 2 5 621749 2 5 1000.99
10 1 1 1 1 5 35023 2 6 1000.99 1 2 1 1 4 2 197591 2 6 1000.99
3 1 3 3 2 2 470108 2 6 1000.99 1 2 1 3 4 5 519986 2 6 1000.99
5 2 1 1 1 3 634148 2 6 1000.99 1 1 1 1 2 5 730430 2 6 1000.99
0 1 1 4 2 2 770625 2 6 1000.99 4 1 1 3 3 4 961252 2 6 1000.99
0 2 1 4 4 4 891820 2 7 1000.99 1 1 1 4 2 4 969999 2 7 1000.99
2 1 1 4 3 4 173053 2 8 1000.99 3 2 2 1 2 2 354931 2 8 1000.99
1 2 1 4 4 3 415848 2 8 1000.99 1 2 1 1 4 3 601372 2 8 1000.99
1 2 1 3 2 4 937558 2 8 1000.99 0 1 1 2 2 1 386038 2 9 1000.99
4 2 1 4 4 2 597841 2 9 1000.99 0 1 1 4 2 5 114197 2 10 1000.99
1 1 1 4 2 3 118060 2 10 1000.99 1 1 1 1 4 1 217419 2 10 1000.99
1 2 1 2 4 1 342760 2 10 1000.99 2 1 1 1 4 1 434970 2 10 1000.99
1 1 1 4 2 2 510121 2 10 1000.99 0 1 1 4 4 1 844822 2 10 1000.99
0 2 2 3 2 3 107733 2 11 1000.99 2 1 1 4 2 5 141489 2 11 1000.99
3 1 1 4 2 4 233177 2 11 1000.99 1 1 1 1 2 5 390433 2 11 1000.99
0 1 1 2 2 2 409285 2 11 1000.99 2 1 1 3 3 4 569731 2 11 1000.99
3 1 1 3 4 2 656062 2 11 1000.99 0 2 1 3 2 4 968459 2 11 1000.99
2 2 3 2 2 2 665243 2 12 1000.99 2 2 1 4 2 4 337880 2 13 1000.99
0 2 1 4 2 4 475625 2 13 1000.99 2 1 1 2 2 3 615428 2 13 1000.99
10 2 2 1 1 2 730881 2 13 1000.99 13 2 1 2 1 2 845624 2 13 1000.99
1 2 1 1 2 4 892901 2 13 1000.99 2 2 1 3 2 5 132975 2 14 1000.99
1 1 2 3 2 4 395501 2 14 1000.99 2 1 1 3 2 5 441404 2 14 1000.99
3 1 1 3 2 2 667836 2 14 1000.99 0 2 1 1 2 3 66196 2 15 1000.99
1 1 1 2 2 5 91356 2 15 1000.99 1 1 1 4 2 4 96501 2 15 1000.99
2 2 1 3 2 3 375625 2 15 1000.99 8 1 1 4 1 4 414371 2 15 1000.99
2 1 1 4 1 5 674621 2 15 1000.99 3 1 1 3 2 2 695960 2 15 1000.99
0 1 1 3 2 3 908548 2 15 1000.99 0 1 3 4 4 5 41061 2 16 1000.99
3 2 1 3 3 5 87349 2 16 1000.99 1 1 1 1 2 4 312337 2 16 1000.99
3 1 3 4 4 1 552496 2 16 1000.99 4 1 1 3 2 2 601344 2 16 1000.99
3 1 1 3 2 5 612685 2 16 1000.99 0 2 1 1 4 5 787362 2 16 1000.99
1 1 2 4 2 3 811432 2 16 1000.99 8 1 1 4 1 4 976761 2 16 1000.99
1 1 1 2 2 2 301987 2 17 1000.99 1 1 1 3 4 5 377285 2 17 1000.99
0 1 1 1 4 4 429968 2 17 1000.99 1 2 1 3 4 5 475580 2 17 1000.99
2 1 1 4 3 1 957214 2 18 1000.99 2 2 1 3 2 4 340995 2 19 1000.99
2 2 1 3 3 1 443842 2 19 1000.99 9 2 1 4 1 2 574864 2 19 1000.99
0 2 1 3 2 3 997681 2 19 1000.99 3 1 1 1 1 3 86209 2 20 1000.99
3 2 1 4 2 5 283918 2 20 1000.99 2 1 1 3 2 5 417057 2 20 1000.99
0 2 1 2 2 3 277011 2 21 1000.99 8 2 1 3 1 4 293375 2 21 1000.99
2 1 1 2 2 4 669282 2 21 1000.99 1 1 1 4 2 4 20886 2 22 1000.99
1 2 2 2 2 3 51659 2 22 1000.99 3 1 1 3 2 5 204168 2 22 1000.99
0 1 1 3 4 2 227469 2 22 1000.99 11 1 2 1 1 5 308086 2 22 1000.99
4 1 2 1 3 2 866568 2 22 1000.99 1 1 1 4 2 3 881079 2 22 1000.99
3 1 1 3 2 3 917849 2 22 1000.99 2 2 1 3 2 2 994768 2 22 1000.99
1 2 1 3 4 2 288147 2 23 1000.99 1 2 3 3 2 2 439337 2 23 1000.99
9 2 3 4 1 5 531604 2 23 1000.99 1 2 3 1 2 3 854974 2 23 1000.99
4 1 2 1 2 3 210498 2 24 1000.99 1 1 1 4 2 4 298144 2 24 1000.99
1 2 1 4 2 2 351473 2 24 1000.99 2 2 1 4 3 3 680146 2 24 1000.99
4 2 1 2 2 2 821825 2 24 1000.99 3 2 2 3 2 4 899954 2 24 1000.99
0 1 1 3 4 3 263496 2 25 1000.99 2 2 1 4 2 4 303866 2 25 1000.99
2 1 1 4 2 4 390812 2 25 1000.99 1 2 1 3 2 5 428283 2 25 1000.99
3 1 1 2 2 2 477432 2 25 1000.99 1 2 2 1 4 4 488217 2 25 1000.99
3 1 1 3 3 3 909124 2 25 1000.99 2 2 1 1 3 2 476969 2 26 1000.99
1 2 1 1 2 2 636489 2 26 1000.99 1 1 1 3 4 2 659865 2 26 1000.99
1 1 1 1 2 5 733715 2 26 1000.99 3 1 1 4 4 4 2826 2 27 1000.99
2 2 1 1 2 5 742421 2 27 1000.99 1 1 1 3 2 4 769330 2 27 1000.99
0 1 1 1 4 3 253147 2 28 1000.99 4 1 1 2 1 4 326795 2 28 1000.99
3 1 1 3 2 5 398311 2 28 1000.99 0 2 1 4 2 1 479384 2 28 1000.99
1 2 1 3 2 4 493199 2 28 1000.99 0 2 1 4 2 4 542943 2 28 1000.99
1 2 1 2 2 2 743440 2 28 1000.99 4 1 1 2 2 2 883725 2 28 1000.99
1 1 1 4 4 1 907484 2 28 1000.99 1 2 1 1 2 5 121357 2 29 1000.99
8 1 2 2 1 2 345491 2 29 1000.99 2 2 1 4 4 5 630604 2 29 1000.99
1 2 1 4 2 2 942561 2 29 1000.99 0 2 2 4 3 4 980123 2 29 1000.99
4 1 1 3 2 3 992975 2 29 1000.99 0 1 1 1 2 2 120707 2 30 1000.99
4 2 1 2 4 1 599662 2 30 1000.99 3 1 1 4 3 2 706806 2 30 1000.99
5 2 1 1 4 5 786072 2 30 1000.99 1 1 1 3 4 3 13148 2 31 1000.99
2 2 1 3 2 1 241523 2 31 1000.99 1 1 1 3 4 2 750924 2 31 1000.99
0 2 1 1 2 2 19016 2 32 1000.99 2 2 1 4 2 4 573331 2 32 1000.99
2 2 1 3 2 2 603617 2 32 1000.99 3 2 1 1 2 5 693503 2 32 1000.99
2 2 2 3 2 2 704454 2 32 1000.99 3 2 1 3 3 3 771110 2 32 1000.99
1 2 1 4 4 5 773754 2 32 1000.99 1 2 1 2 2 1 144583 2 33 1000.99
0 2 1 3 4 4 197703 2 33 1000.99 3 2 1 4 4 5 178606 2 34 1000.99
0 2 1 3 2 4 273101 2 34 1000.99 8 1 3 2 1 1 508459 2 34 1000.99
2 2 1 3 2 5 555590 2 34 1000.99 3 2 2 4 2 1 961202 2 34 1000.99
2 1 1 3 2 5 478258 2 35 1000.99 0 2 1 3 3 3 726325 2 35 1000.99
0 1 1 4 2 2 901093 2 35 1000.99 3 2 1 3 2 1 101453 2 36 1000.99
3 2 3 3 2 2 130407 2 36 1000.99 10 1 1 3 1 5 333408 2 36 1000.99
4 1 1 3 2 4 784856 2 36 1000.99 2 2 2 2 2 3 866724 2 36 1000.99
2 1 3 4 2 4 27578 2 37 1000.99 2 1 1 3 3 4 166230 2 37 1000.99
1 1 1 2 2 2 241280 2 37 1000.99 5 1 1 1 1 4 456428 2 37 1000.99
4 1 1 2 2 2 590297 2 37 1000.99 2 1 1 3 2 5 894552 2 37 1000.99
1 2 1 3 3 1 65763 2 38 1000.99 0 1 1 2 3 2 247288 2 38 1000.99
1 1 2 4 2 4 567506 2 38 1000.99 1 1 1 4 2 2 916432 2 38 1000.99
0 2 2 3 2 4 442871 2 39 1000.99 5 2 1 2 4 5 517566 2 39 1000.99
7 2 1 2 1 3 526107 2 39 1000.99 4 2 1 3 3 3 582509 2 39 1000.99
6 1 1 2 2 4 687269 2 39 1000.99 7 1 3 3 1 3 741437 2 39 1000.99
0 2 2 4 2 5 49356 2 40 1000.99 5 1 1 1 1 3 50118 2 40 1000.99
2 2 1 1 4 5 133987 2 40 1000.99 1 1 1 2 2 5 150096 2 40 1000.99
1 1 1 2 4 5 261146 2 40 1000.99 1 2 1 3 2 5 328658 2 40 1000.99
2 2 1 3 2 3 422440 2 40 1000.99 5 2 1 2 1 2 430482 2 40 1000.99
0 2 1 3 2 5 469440 2 40 1000.99 4 1 3 4 3 5 585573 2 40 1000.99
1 1 1 4 4 4 716656 2 40 1000.99 1 1 1 4 4 4 780490 2 40 1000.99
2 1 1 2 2 5 825369 2 40 1000.99 0 2 1 3 2 5 209080 3 1 999.96
0 1 1 2 2 2 381472 3 1 999.96 3 2 1 3 2 2 507293 3 1 999.96
0 2 1 3 2 5 531526 3 1 999.96 1 2 1 2 2 1 788791 3 1 999.96
2 2 3 2 4 3 266211 3 2 999.96 0 2 1 4 2 2 442624 3 2 999.96
1 1 2 1 3 2 448591 3 2 999.96 1 1 3 3 4 2 482688 3 2 999.96
2 2 1 4 4 3 499726 3 2 999.96 2 2 3 3 2 4 513283 3 2 999.96
0 1 1 3 2 2 689710 3 2 999.96 10 2 1 4 1 3 868515 3 2 999.96
1 2 1 3 2 3 974149 3 2 999.96 2 1 1 2 4 5 995181 3 2 999.96
8 2 2 4 1 4 256859 3 3 999.96 0 1 1 4 4 2 658908 3 3 999.96
1 1 1 1 2 3 786102 3 3 999.96 2 1 1 4 2 4 901958 3 3 999.96
2 2 1 3 4 2 142850 3 4 999.96 2 1 1 4 4 2 260754 3 4 999.96
2 2 1 1 2 3 382009 3 4 999.96 4 1 2 3 2 4 473471 3 4 999.96
0 2 1 2 2 3 538642 3 4 999.96 1 2 1 3 2 5 659341 3 4 999.96
4 2 1 2 1 3 320328 3 5 999.96 1 2 1 1 2 4 375058 3 5 999.96
0 2 1 2 2 1 898136 3 5 999.96 1 1 1 4 2 3 84803 3 6 999.96
1 2 1 1 2 4 90471 3 6 999.96 0 2 1 1 2 5 493270 3 6 999.96
0 2 2 3 2 5 886367 3 6 999.96 3 2 1 4 1 5 51055 3 7 999.96
5 1 1 3 2 2 203571 3 7 999.96 1 1 1 1 2 2 246749 3 7 999.96
1 2 1 2 2 3 362562 3 7 999.96 4 2 3 3 2 4 562049 3 7 999.96
2 1 1 2 2 2 641883 3 7 999.96 1 1 1 4 2 4 729591 3 7 999.96
0 1 1 3 2 5 552184 3 8 999.96 0 1 1 4 2 5 664339 3 8 999.96
1 1 1 3 2 3 825273 3 8 999.96 2 1 2 2 4 2 274557 3 9 999.96
3 2 1 2 3 3 377856 3 9 999.96 0 1 1 4 3 2 379195 3 9 999.96
1 1 1 4 4 2 531175 3 9 999.96 2 2 1 1 4 5 540390 3 9 999.96
0 1 1 2 4 5 552880 3 9 999.96 2 1 1 1 2 5 838951 3 9 999.96
2 1 1 2 2 2 18614 3 10 999.96 0 2 1 2 3 4 118181 3 10 999.96
2 1 1 1 2 4 254888 3 10 999.96 2 1 1 3 2 5 389657 3 10 999.96
1 1 1 4 3 1 712515 3 10 999.96 1 2 2 3 2 4 714640 3 10 999.96
1 1 1 3 2 2 172942 3 11 999.96 2 2 2 4 2 5 299444 3 11 999.96
3 2 1 1 4 5 514298 3 11 999.96 1 2 1 4 2 2 984598 3 11 999.96
0 2 1 1 2 2 80807 3 12 999.96 1 1 1 2 3 3 29608 3 13 999.96
2 1 1 4 3 1 109303 3 13 999.96 2 1 1 4 4 2 455326 3 13 999.96
1 1 1 3 4 2 526762 3 13 999.96 4 1 1 1 1 3 577313 3 13 999.96
4 1 1 1 2 4 895966 3 13 999.96 0 1 1 2 4 4 291520 3 14 999.96
2 1 3 3 4 2 351342 3 14 999.96 1 2 1 4 2 3 465132 3 14 999.96
1 1 3 2 4 2 848590 3 14 999.96 1 2 1 2 2 5 309027 3 15 999.96
0 2 1 3 3 2 661132 3 15 999.96 2 1 1 3 4 5 674171 3 15 999.96
4 2 1 3 2 4 877910 3 15 999.96 0 2 3 4 2 4 888272 3 15 999.96
2 1 2 3 2 5 50937 3 16 999.96 2 2 1 4 2 2 71514 3 16 999.96
0 2 1 4 2 4 74528 3 16 999.96 3 1 1 3 2 3 235261 3 16 999.96
1 1 1 4 2 3 263917 3 16 999.96 4 2 1 3 4 3 816932 3 16 999.96
0 2 1 2 2 5 962876 3 16 999.96 2 2 1 2 2 5 977411 3 16 999.96
1 1 1 1 4 3 256671 3 17 999.96 0 1 3 3 4 3 431807 3 17 999.96
3 2 1 4 3 2 474818 3 17 999.96 0 1 1 3 4 2 498284 3 17 999.96
3 1 1 3 4 3 569460 3 17 999.96 4 2 1 1 2 1 743559 3 17 999.96
2 2 1 1 4 5 153822 3 18 999.96 8 2 2 3 1 5 361232 3 18 999.96
0 2 1 2 3 4 479951 3 18 999.96 1 2 1 4 2 4 621026 3 18 999.96
2 2 2 3 2 3 995933 3 18 999.96 2 1 1 2 2 4 96302 3 19 999.96
2 1 1 1 3 3 322778 3 19 999.96 1 1 3 4 2 2 619519 3 19 999.96
12 2 1 2 1 2 114023 3 20 999.96 4 1 1 3 4 2 149791 3 20 999.96
0 2 3 4 2 4 858118 3 20 999.96 0 1 1 1 2 3 892410 3 20 999.96
3 2 1 2 1 5 912710 3 20 999.96 3 1 1 1 2 3 994897 3 20 999.96
2 1 1 1 2 3 292454 3 21 999.96 1 1 1 3 2 4 335751 3 21 999.96
0 2 1 4 3 3 959332 3 21 999.96 1 1 1 1 4 3 971785 3 21 999.96
2 1 1 4 2 1 997335 3 21 999.96 2 2 1 4 2 2 6008 3 22 999.96
4 1 1 4 4 1 236189 3 22 999.96 0 1 1 4 2 3 554163 3 22 999.96
2 1 1 4 2 5 817555 3 22 999.96 0 2 1 4 2 2 192867 3 23 999.96
10 1 2 4 1 5 645837 3 23 999.96 2 2 1 3 3 4 940238 3 23 999.96
5 2 1 3 3 5 949084 3 23 999.96 1 1 1 4 4 4 243003 3 24 999.96
0 2 1 1 4 3 680615 3 24 999.96 1 1 3 2 4 5 979584 3 24 999.96
2 2 1 4 2 3 61728 3 25 999.96 6 2 1 3 4 3 161546 3 25 999.96
1 1 2 4 2 2 405680 3 25 999.96 1 2 1 4 2 3 417724 3 25 999.96
1 1 1 3 4 1 505691 3 25 999.96 1 1 2 1 2 2 515576 3 25 999.96
1 2 1 3 2 2 546070 3 25 999.96 4 1 2 2 4 4 657602 3 25 999.96
2 1 1 3 2 5 862373 3 25 999.96 2 2 1 4 2 4 937347 3 25 999.96
0 1 1 2 2 5 991045 3 25 999.96 6 1 1 4 1 1 357088 3 26 999.96
5 1 1 3 3 4 541192 3 26 999.96 0 2 3 2 2 2 541673 3 26 999.96
4 2 3 4 1 4 653146 3 26 999.96 2 2 1 3 3 4 378563 3 27 999.96
1 1 1 4 3 4 616758 3 27 999.96 0 2 1 3 2 3 879276 3 27 999.96
9 1 1 1 1 3 913040 3 27 999.96 2 2 1 3 4 4 511602 3 28 999.96
1 1 2 1 2 2 537845 3 28 999.96 1 1 1 3 2 4 40622 3 29 999.96
0 1 1 1 4 2 817345 3 29 999.96 0 2 1 2 4 5 200955 3 30 999.96
4 1 1 3 2 5 584953 3 30 999.96 1 2 1 1 2 4 617022 3 30 999.96
6 1 2 3 4 2 757172 3 30 999.96 0 1 1 1 4 4 813081 3 30 999.96
2 1 1 3 2 3 822865 3 30 999.96 0 1 1 1 3 5 847941 3 30 999.96
1 1 2 3 2 4 162591 3 31 999.96 2 2 1 3 3 2 294066 3 31 999.96
3 1 1 4 2 4 310181 3 31 999.96 1 1 2 3 2 2 626979 3 31 999.96
5 2 1 2 2 5 676054 3 31 999.96 0 1 2 4 4 1 54626 3 32 999.96
3 2 1 3 2 3 165170 3 32 999.96 5 2 3 3 3 1 368968 3 32 999.96
3 2 1 2 2 1 377328 3 32 999.96 2 2 2 3 3 3 649120 3 33 999.96
2 1 1 3 2 2 738578 3 33 999.96 4 2 1 3 4 4 814182 3 33 999.96
1 1 1 3 4 5 831502 3 33 999.96 2 1 1 3 2 1 986270 3 33 999.96
1 2 1 1 3 3 13110 3 34 999.96 1 2 1 2 2 3 241330 3 34 999.96
3 2 1 3 4 4 311259 3 34 999.96 0 2 1 4 3 3 440991 3 34 999.96
1 2 1 2 3 5 602632 3 34 999.96 1 1 1 2 2 4 648852 3 34 999.96
2 2 1 4 4 4 113333 3 35 999.96 1 1 3 2 4 3 246325 3 35 999.96
4 1 3 3 2 2 408742 3 35 999.96 1 1 3 2 2 5 718665 3 35 999.96
1 2 1 3 3 5 253164 3 36 999.96 0 2 3 1 2 3 446996 3 36 999.96
4 2 1 2 1 3 690575 3 36 999.96 4 1 1 3 2 2 745715 3 36 999.96
1 2 1 3 2 5 753247 3 36 999.96 4 1 1 3 4 2 798555 3 36 999.96
4 1 1 4 2 4 814817 3 36 999.96 1 2 1 4 2 3 994277 3 36 999.96
1 1 1 3 3 5 134955 3 37 999.96 3 2 1 3 4 4 289588 3 37 999.96
0 2 1 4 2 5 652796 3 37 999.96 4 2 3 2 1 5 802367 3 37 999.96
0 2 3 3 4 5 882174 3 37 999.96 2 1 1 3 2 5 975938 3 37 999.96
0 1 1 3 2 4 27837 3 38 999.96 4 2 1 3 2 5 339484 3 38 999.96
1 2 1 3 3 1 358810 3 38 999.96 0 2 1 3 2 4 730924 3 38 999.96
2 1 2 4 4 3 755268 3 38 999.96 5 1 1 4 1 5 359709 3 39 999.96
3 2 1 4 2 4 515728 3 39 999.96 2 1 1 4 2 5 618704 3 39 999.96
1 1 1 2 2 1 724733 3 39 999.96 0 1 1 3 2 2 806377 3 39 999.96
2 1 3 3 2 4 994798 3 39 999.96 2 2 1 1 2 5 667576 3 40 999.96
6 1 1 2 1 2 732221 3 40 999.96 4 2 1 3 2 2 380844 4 1 998.22
1 1 2 3 4 3 508120 4 1 998.22 1 2 1 3 3 4 512105 4 1 998.22
1 1 1 3 2 4 644281 4 1 998.22 4 1 1 4 2 2 670794 4 1 998.22
12 1 2 3 1 4 796746 4 1 998.22 0 2 1 2 4 4 223582 4 2 998.22
2 1 1 3 4 5 653969 4 2 998.22 1 1 1 1 2 2 931106 4 2 998.22
2 2 1 3 3 1 4334 4 3 998.22 2 1 3 3 2 2 145443 4 3 998.22
1 2 1 4 4 3 327290 4 3 998.22 2 2 1 3 2 3 353675 4 3 998.22
0 1 1 4 2 2 527148 4 3 998.22 12 1 1 3 1 2 620343 4 3 998.22
0 1 1 2 2 4 794585 4 3 998.22 2 2 2 3 2 4 151189 4 4 998.22
3 1 1 3 3 3 310454 4 4 998.22 8 2 1 2 1 3 544327 4 4 998.22
0 1 1 2 2 3 570276 4 4 998.22 2 1 3 3 4 1 93200 4 5 998.22
1 2 1 4 2 3 94887 4 5 998.22 0 1 1 2 2 1 128644 4 5 998.22
2 2 1 4 2 2 222690 4 5 998.22 2 1 1 2 4 5 456110 4 5 998.22
3 1 1 1 4 1 618834 4 5 998.22 1 1 1 4 4 5 946990 4 5 998.22
11 1 1 3 1 2 9106 4 6 998.22 1 1 1 4 2 5 202885 4 6 998.22
1 2 1 3 2 5 343160 4 6 998.22 1 2 1 4 2 3 526891 4 7 998.22
1 1 1 4 4 4 682736 4 7 998.22 2 1 1 3 2 4 161504 4 8 998.22
1 1 1 3 2 4 394269 4 8 998.22 4 2 1 3 2 2 694755 4 8 998.22
2 2 1 4 4 5 752400 4 8 998.22 3 1 1 4 4 2 938446 4 8 998.22
4 1 2 3 2 3 312242 4 9 998.22 14 1 1 3 1 4 318484 4 9 998.22
3 2 2 1 2 5 321344 4 9 998.22 1 2 1 3 2 4 372362 4 9 998.22
1 1 2 4 2 4 885739 4 9 998.22 0 2 1 2 4 5 993156 4 9 998.22
0 2 1 3 2 5 85331 4 10 998.22 8 1 1 4 1 3 227885 4 10 998.22
0 2 3 1 2 1 234867 4 10 998.22 3 1 1 2 4 2 243226 4 10 998.22
2 2 1 3 4 5 458412 4 10 998.22 1 2 1 3 3 5 683785 4 10 998.22
1 2 1 3 2 5 708572 4 10 998.22 3 2 1 4 2 2 927732 4 10 998.22
3 1 1 3 2 2 72199 4 11 998.22 1 2 1 2 2 4 348589 4 11 998.22
0 1 1 1 4 3 420090 4 11 998.22 1 1 1 3 4 4 622973 4 11 998.22
3 2 1 2 4 1 701156 4 11 998.22 1 1 1 1 2 4 776993 4 11 998.22
4 2 2 3 2 3 817562 4 11 998.22 2 1 1 3 3 4 369339 4 12 998.22
3 1 1 3 2 2 935888 4 12 998.22 10 2 1 3 1 4 345837 4 13 998.22
1 1 1 3 2 1 414147 4 13 998.22 1 2 1 4 2 3 484068 4 13 998.22
7 2 2 3 1 3 642791 4 13 998.22 2 2 1 4 2 5 692030 4 13 998.22
3 2 1 3 4 5 835036 4 13 998.22 4 2 2 3 4 4 875474 4 13 998.22
0 1 1 4 4 1 916174 4 13 998.22 2 2 1 4 2 5 285559 4 14 998.22
2 2 1 2 2 1 348803 4 14 998.22 0 1 1 2 3 4 438263 4 14 998.22
2 2 1 4 2 3 498296 4 14 998.22 3 1 1 4 4 4 575216 4 14 998.22
3 2 1 4 2 4 952733 4 14 998.22 3 2 1 3 4 1 76246 4 15 998.22
8 1 3 3 1 5 616380 4 15 998.22 0 2 3 3 2 4 819444 4 15 998.22
2 1 1 4 2 2 103114 4 16 998.22 1 1 2 2 2 3 401887 4 16 998.22
4 2 1 3 1 2 563087 4 16 998.22 0 2 3 2 2 3 611940 4 16 998.22
1 2 1 4 2 2 770658 4 16 998.22 0 2 1 4 2 2 72632 4 17 998.22
3 1 1 4 4 4 88932 4 17 998.22 1 2 3 3 3 4 168300 4 17 998.22
2 1 1 2 2 2 833000 4 17 998.22 4 1 1 3 2 5 882088 4 17 998.22
3 2 3 1 3 3 119561 4 18 998.22 4 1 1 2 4 1 409373 4 18 998.22
3 2 3 3 3 4 605204 4 18 998.22 3 2 1 2 2 3 704002 4 18 998.22
4 1 1 3 2 4 45333 4 19 998.22 1 2 1 4 2 5 63301 4 19 998.22
3 1 1 2 2 2 200539 4 19 998.22 1 2 1 3 2 3 394484 4 19 998.22
0 1 1 4 2 4 822509 4 19 998.22 2 1 2 1 2 5 851337 4 19 998.22
0 2 1 4 2 4 866914 4 19 998.22 1 2 1 4 2 3 992237 4 19 998.22
2 1 1 3 2 4 46630 4 20 998.22 2 2 1 2 4 4 231608 4 20 998.22
2 2 1 3 2 1 261191 4 20 998.22 3 2 2 1 4 3 362765 4 20 998.22
1 2 2 1 4 2 726440 4 20 998.22 2 2 2 3 3 2 884638 4 20 998.22
3 1 1 1 2 4 905866 4 20 998.22 2 1 1 4 2 5 972054 4 20 998.22
2 2 1 1 2 2 624776 4 21 998.22 2 2 1 3 2 3 836630 4 21 998.22
1 1 1 4 2 2 901047 4 21 998.22 3 1 1 3 4 3 913838 4 21 998.22
0 2 1 3 2 4 914930 4 21 998.22 2 2 1 2 3 5 7311 4 22 998.22
1 2 3 3 3 2 631529 4 22 998.22 6 2 1 4 1 4 761881 4 22 998.22
3 2 1 4 2 4 303997 4 23 998.22 1 2 1 4 2 2 340764 4 23 998.22
1 1 1 2 4 4 363728 4 23 998.22 3 2 1 1 2 3 449358 4 23 998.22
1 1 1 4 2 4 519687 4 23 998.22 7 1 1 1 1 4 715036 4 23 998.22
4 2 1 3 4 3 773064 4 23 998.22 0 1 1 3 4 4 303229 4 24 998.22
3 2 1 3 2 4 439288 4 24 998.22 2 2 1 3 2 3 681268 4 24 998.22
4 2 1 1 2 2 686238 4 24 998.22 0 2 1 2 2 5 709305 4 24 998.22
1 1 1 2 2 3 8634 4 25 998.22 6 1 1 4 2 3 481326 4 25 998.22
1 1 3 1 4 2 543566 4 25 998.22 1 1 1 4 2 2 861082 4 25 998.22
1 2 3 3 2 4 83192 4 26 998.22 2 1 1 3 4 2 100973 4 26 998.22
3 2 1 3 2 4 465212 4 26 998.22 3 1 1 3 4 3 841719 4 26 998.22
4 1 2 3 2 2 357835 4 27 998.22 1 2 2 3 2 4 423481 4 27 998.22
2 2 1 2 2 2 249974 4 28 998.22 4 1 1 2 4 3 279181 4 28 998.22
4 2 1 2 2 2 406445 4 28 998.22 2 1 2 3 2 3 75287 4 29 998.22
2 1 1 4 2 4 118043 4 29 998.22 1 2 1 2 2 5 175592 4 29 998.22
3 1 1 3 2 3 354480 4 29 998.22 1 1 1 2 2 1 460511 4 29 998.22
0 2 1 4 2 2 566430 4 29 998.22 3 1 1 4 4 5 611651 4 29 998.22
4 2 1 1 2 5 647701 4 29 998.22 3 2 1 1 2 3 691686 4 29 998.22
3 2 2 3 3 2 835657 4 29 998.22 4 1 2 4 2 1 359559 4 30 998.22
1 2 2 3 4 4 464778 4 30 998.22 1 2 1 3 2 2 487872 4 30 998.22
1 1 1 2 2 4 564411 4 30 998.22 5 2 1 3 4 4 95965 4 31 998.22
8 2 1 3 1 2 430429 4 31 998.22 0 1 1 4 4 2 597691 4 31 998.22
3 1 2 2 2 2 756374 4 31 998.22 2 1 2 4 2 5 444574 4 32 998.22
1 1 1 3 2 4 626062 4 32 998.22 0 1 1 3 4 4 28450 4 33 998.22
1 2 1 3 4 2 561664 4 33 998.22 3 2 1 2 1 4 6534 4 34 998.22
6 1 1 4 4 2 57553 4 34 998.22 3 1 2 3 2 4 665516 4 34 998.22
1 2 1 4 2 4 152961 4 35 998.22 3 2 1 4 2 2 192871 4 35 998.22
0 1 1 4 2 2 741373 4 35 998.22 2 1 3 3 2 5 822312 4 35 998.22
1 2 1 2 2 4 907571 4 35 998.22 0 1 2 3 2 4 108807 4 36 998.22
2 1 1 2 2 2 114030 4 36 998.22 3 2 1 1 4 2 346461 4 36 998.22
1 1 1 2 2 2 432462 4 36 998.22 2 1 1 3 2 2 818752 4 36 998.22
0 2 1 3 3 2 913639 4 36 998.22 3 1 2 3 4 2 71492 4 37 998.22
1 2 1 1 4 5 156835 4 37 998.22 2 1 1 4 4 4 533296 4 37 998.22
3 2 1 3 2 2 589208 4 37 998.22 0 1 2 2 2 2 658644 4 37 998.22
1 1 2 3 2 2 745244 4 37 998.22 3 2 1 3 2 1 865255 4 37 998.22
1 2 1 4 2 2 900025 4 37 998.22 1 2 2 2 2 2 97378 4 38 998.22
2 2 1 3 2 4 156544 4 38 998.22 4 2 3 2 1 2 175709 4 38 998.22
2 2 1 3 2 3 398226 4 38 998.22 3 1 1 2 2 2 647754 4 39 998.22
0 2 1 1 2 4 733627 4 39 998.22 2 1 1 4 2 2 785555 4 39 998.22
2 1 2 4 3 2 951720 4 39 998.22 0 1 3 3 2 4 997381 4 39 998.22
3 1 1 3 3 5 111764 4 40 998.22 1 2 2 4 4 2 150332 4 40 998.22
2 2 1 4 2 4 282474 4 40 998.22 1 2 1 2 2 4 677474 4 40 998.22
1 2 1 3 2 5 677698 4 40 998.22 0 2 1 2 2 4 772600 4 40 998.22
0 2 1 4 2 4 948571 4 40 998.22 1 1 2 2 3 3 969195 4 40 998.22
0 1 3 4 2 2 303457 5 1 998.26 1 1 1 3 4 4 24253 5 2 998.26
1 1 1 1 2 4 153514 5 2 998.26 3 2 1 3 2 3 306193 5 2 998.26
1 2 1 1 4 2 502224 5 2 998.26 0 1 3 3 4 1 926493 5 2 998.26
1 2 1 3 3 4 36049 5 3 998.26 1 1 1 1 3 3 813091 5 3 998.26
2 2 1 1 4 5 52728 5 4 998.26 0 2 3 2 2 4 180687 5 4 998.26
0 1 1 2 2 5 502549 5 4 998.26 2 2 1 4 3 2 523768 5 4 998.26
3 2 2 4 2 3 539351 5 4 998.26 2 1 1 2 2 3 124439 5 5 998.26
2 2 1 1 4 4 195475 5 5 998.26 3 1 1 2 2 2 651509 5 5 998.26
0 2 2 4 2 4 698527 5 5 998.26 2 2 3 1 3 5 711708 5 5 998.26
2 2 1 4 2 2 884839 5 5 998.26 5 1 1 3 1 4 185376 5 6 998.26
1 2 1 2 2 5 600765 5 6 998.26 2 2 1 4 3 3 880345 5 6 998.26
4 2 1 1 2 5 125623 5 7 998.26 1 1 1 4 2 2 176695 5 7 998.26
1 1 1 1 2 4 600503 5 7 998.26 0 2 1 3 4 2 744322 5 7 998.26
1 1 1 3 2 2 963617 5 7 998.26 2 1 1 2 4 4 979370 5 7 998.26
3 2 1 1 3 2 30945 5 8 998.26 2 1 3 4 2 3 270848 5 8 998.26
1 1 1 2 3 2 399693 5 8 998.26 2 1 1 4 3 4 461061 5 8 998.26
4 1 1 2 3 3 570913 5 8 998.26 1 2 1 2 3 3 573922 5 8 998.26
4 2 1 3 4 3 605756 5 8 998.26 0 2 1 4 3 2 776019 5 8 998.26
9 1 1 2 1 2 862030 5 8 998.26 0 2 1 4 2 2 459 5 9 998.26
1 2 1 2 2 5 43429 5 9 998.26 2 2 1 2 2 4 292130 5 9 998.26
5 1 1 2 2 3 366240 5 9 998.26 7 2 1 2 1 3 488521 5 9 998.26
1 2 1 4 4 5 800418 5 9 998.26 2 2 1 4 2 2 5962 5 10 998.26
3 1 3 3 2 4 228174 5 10 998.26 2 2 1 1 1 4 343452 5 10 998.26
1 1 1 3 2 2 614309 5 10 998.26 1 1 1 3 2 4 724142 5 10 998.26
1 2 1 3 4 3 766359 5 10 998.26 1 2 1 4 4 5 822706 5 10 998.26
1 1 2 1 3 3 889308 5 10 998.26 1 1 1 2 2 1 22339 5 11 998.26
0 1 1 3 4 5 207409 5 12 998.26 1 1 2 1 3 2 398671 5 12 998.26
3 1 1 1 1 3 482282 5 12 998.26 3 2 3 3 4 2 565431 5 12 998.26
3 1 3 1 2 1 663079 5 12 998.26 1 2 1 4 2 5 780726 5 12 998.26
5 2 1 3 1 3 895 5 13 998.26 6 1 2 2 4 2 65996 5 13 998.26
5 2 1 2 1 1 350647 5 13 998.26 2 1 1 3 4 5 371762 5 13 998.26
1 2 1 3 4 2 478022 5 13 998.26 0 2 2 3 2 3 602099 5 13 998.26
4 1 2 1 2 2 126331 5 14 998.26 2 1 1 4 4 5 397371 5 14 998.26
3 1 3 3 2 2 581000 5 14 998.26 0 2 1 4 2 5 80663 5 15 998.26
4 2 1 1 2 2 99081 5 15 998.26 8 1 1 4 1 5 99349 5 15 998.26
3 1 1 3 4 5 169493 5 15 998.26 1 2 1 4 2 3 190376 5 15 998.26
1 1 1 4 4 5 214163 5 15 998.26 0 1 1 3 2 1 904528 5 15 998.26
1 2 1 4 4 2 51729 5 16 998.26 0 2 3 3 2 2 91083 5 16 998.26
4 2 3 4 2 3 96746 5 16 998.26 3 1 1 3 2 4 178267 5 16 998.26
3 2 1 3 2 5 190103 5 16 998.26 3 2 1 3 2 2 441385 5 16 998.26
2 2 1 2 3 2 804588 5 16 998.26 2 2 1 3 2 4 847192 5 16 998.26
3 1 1 2 2 2 108679 5 17 998.26 3 1 1 4 2 2 196302 5 17 998.26
5 1 3 3 1 1 243355 5 17 998.26 7 2 1 1 2 4 245947 5 17 998.26
4 2 1 3 2 3 338950 5 17 998.26 1 1 3 3 2 4 538067 5 17 998.26
1 2 1 2 4 2 29833 5 18 998.26 4 2 1 4 4 2 154421 5 18 998.26
2 1 1 3 4 2 173912 5 18 998.26 3 1 1 3 2 3 192921 5 18 998.26
2 1 1 1 3 5 202107 5 18 998.26 1 2 1 3 2 4 515032 5 18 998.26
0 1 1 3 2 4 537517 5 18 998.26 0 2 1 4 2 4 721978 5 18 998.26
0 2 1 4 4 5 728977 5 18 998.26 3 2 1 1 4 3 273536 5 19 998.26
1 2 1 1 4 1 854985 5 19 998.26 1 2 1 2 2 1 878153 5 19 998.26
2 2 1 2 2 2 93627 5 20 998.26 1 1 1 3 2 2 264296 5 20 998.26
1 1 3 4 3 5 973973 5 20 998.26 1 2 2 2 4 4 997536 5 20 998.26
2 2 1 3 2 2 506652 5 21 998.26 1 1 1 1 2 2 712321 5 21 998.26
4 1 2 4 4 4 762956 5 21 998.26 8 1 1 4 1 3 853516 5 21 998.26
1 2 1 2 2 5 875255 5 21 998.26 0 2 1 4 2 5 930897 5 21 998.26
3 2 1 4 2 4 963337 5 21 998.26 0 1 1 4 2 2 221523 5 22 998.26
4 1 1 3 2 3 257387 5 22 998.26 2 2 1 3 2 4 760053 5 22 998.26
1 2 1 4 4 3 929482 5 22 998.26 3 2 1 4 4 3 53376 5 23 998.26
10 1 1 2 1 2 95812 5 23 998.26 2 2 1 3 2 3 107964 5 23 998.26
3 1 1 2 2 5 264142 5 23 998.26 7 2 1 1 1 2 269216 5 23 998.26
3 2 1 3 4 5 477167 5 23 998.26 1 2 1 3 2 1 744712 5 23 998.26
0 2 1 4 2 2 869567 5 23 998.26 0 2 1 3 2 1 157351 5 24 998.26
0 1 1 1 3 5 266606 5 24 998.26 0 1 1 2 2 5 377528 5 24 998.26
1 1 1 3 3 4 624648 5 24 998.26 2 2 1 4 2 2 762919 5 24 998.26
2 1 1 3 2 4 268479 5 25 998.26 0 2 1 3 2 2 321092 5 25 998.26
3 2 3 4 4 1 346737 5 25 998.26 2 2 1 3 2 2 422000 5 25 998.26
0 1 3 4 3 3 424698 5 25 998.26 2 2 1 4 4 1 461497 5 25 998.26
1 2 1 3 4 1 485779 5 25 998.26 1 1 1 1 4 5 355714 5 26 998.26
3 1 1 4 2 2 734393 5 26 998.26 3 1 2 3 2 2 442499 5 27 998.26
3 1 2 1 2 3 469205 5 27 998.26 3 1 1 2 4 5 483019 5 27 998.26
1 1 2 3 2 2 425217 5 28 998.26 2 2 1 2 4 4 270953 5 29 998.26
1 2 1 2 2 3 442177 5 29 998.26 1 1 3 4 2 3 582920 5 29 998.26
6 2 2 4 1 2 669610 5 29 998.26 3 2 1 4 4 5 29960 5 30 998.26
2 1 1 3 2 2 68595 5 30 998.26 1 1 1 2 2 3 80347 5 30 998.26
1 2 2 4 2 3 485432 5 30 998.26 1 2 1 4 2 3 573849 5 30 998.26
2 2 1 1 2 3 717676 5 30 998.26 1 2 1 1 2 5 165898 5 31 998.26
2 1 1 4 4 2 172439 5 31 998.26 0 1 2 4 4 3 438961 5 31 998.26
3 2 2 4 2 4 988713 5 31 998.26 1 2 1 2 3 2 269125 5 32 998.26
3 2 2 2 2 3 667801 5 32 998.26 1 2 1 3 4 4 799141 5 32 998.26
4 2 1 3 4 5 834100 5 32 998.26 3 2 1 4 2 5 836356 5 32 998.26
0 1 3 4 2 5 897168 5 32 998.26 2 2 1 3 2 3 189498 5 33 998.26
3 2 1 3 2 2 299909 5 33 998.26 0 1 1 4 2 2 513598 5 33 998.26
0 2 1 2 2 4 981942 5 33 998.26 1 2 1 2 2 5 477049 5 34 998.26
3 1 1 4 2 2 736961 5 34 998.26 10 1 1 2 1 3 753045 5 34 998.26
1 1 1 1 2 5 983521 5 34 998.26 2 1 1 3 2 3 35829 5 35 998.26
3 2 1 3 2 2 82504 5 35 998.26 2 2 2 3 2 5 289693 5 35 998.26
4 2 2 2 4 4 443048 5 35 998.26 0 2 3 1 4 4 112296 5 36 998.26
1 2 1 4 4 4 175254 5 36 998.26 1 1 1 4 4 4 337790 5 36 998.26
0 2 1 3 2 4 514466 5 36 998.26 7 2 1 4 1 2 936565 5 36 998.26
1 1 1 2 4 5 31932 5 37 998.26 2 2 2 4 4 2 141854 5 37 998.26
0 2 2 1 2 1 297584 5 37 998.26 1 2 2 4 2 2 399094 5 37 998.26
3 2 1 3 2 4 630960 5 37 998.26 3 1 1 2 2 4 643518 5 37 998.26
5 2 1 3 2 4 736487 5 37 998.26 2 1 1 1 4 5 754663 5 37 998.26
0 2 1 3 2 2 502078 5 38 998.26 0 1 1 3 2 4 547670 5 38 998.26
1 1 1 3 3 3 615842 5 38 998.26 0 2 1 4 2 3 844728 5 38 998.26
2 2 1 4 2 3 216815 5 39 998.26 1 2 1 1 3 5 294697 5 39 998.26
0 1 1 4 3 4 306198 5 39 998.26 0 2 1 3 2 3 540641 5 39 998.26
4 2 1 4 2 1 716476 5 39 998.26 4 2 3 1 1 5 793430 5 39 998.26
2 2 1 1 4 2 244599 5 40 998.26 1 2 1 3 4 4 738928 5 40 998.26
2 2 1 3 2 2 830211 5 40 998.26 3 1 1 3 2 3 920025 5 40 998.26
;
run;

ods select none;
ods output VarianceEstimation=VE Summary=Summary;
proc surveymeans data=counts plots=none
  varmethod=jackknife(outweights=jkweights outjkcoefs=jkcoefs);
  cluster psu;
  strata strata;
  weight SamplingWeight;
  var visits;
run;

data _null_;
  set VE(where=(Label1="Number of Replicates"));
  call symput('replicates',cValue1);
run;

data _null_;
  set Summary;
  if Label1="Number of Strata" then do;
    call symput('H',cValue1);
  end;
run;

ods output ParameterEstimates=FullSample(where=(Parameter ne "Scale")
                                         keep=Parameter Estimate Level1
                                         rename=(Estimate=Estimate0))
           ParameterEstimates=parms(keep=df);
proc genmod data=jkweights;
  class sex race marital private education;
  weight SamplingWeight;
  model visits = sex race marital private education / dist=poisson;
run;

ods output Statistics=statistics;
proc surveymeans data=parms sum;
  var df;
run;

data _null_;
  set statistics;
  call symput('p',Sum);
run;

/* Begin Section 3a */

/**********************************************************************************/
/*                 S A S  W E B  E X A M P L E  M A C R O                         */
/*                                                                                */
/*        NAME: STACK                                                             */
/*       TITLE: Poisson Regressions for Complex Surveys                           */
/*        DATE: July 20, 2105                                                     */   
/*     PRODUCT: STAT                                                              */
/*       PROCS: MCMC                                                              */  
/*     SUPPORT: Allen McDowell                                                    */   
/*                                                                                */
/*   ARGUMENTS:                                                                   */
/*                                                                                */
/*         DATA - the number of jackknife replicates                              */
/*                                                                                */
/* DESCRIPTION:                                                                   */
/* The macro \%STACK reshapes the JKweights data set from wide form to long form. */
/* It creates the variable Replicate, which indexes the R copies of the original  */
/* data, the variable Repweight, which contains the replicate weights, and it     */
/* sorts the newly reshaped data set by the variable Replicate.                   */
/*                                                                                */
/**********************************************************************************/

%macro stack(data=);
  data &data;
    set &data;
    %do i=1 %to &replicates;
      Replicate=&i;
      Repweight=RepWt_&i;
      output;
    %end;
  run;

  proc sort data=&data;
    by replicate;
  run;

%mend stack;
%stack(data=jkweights)

ods output ParameterEstimates=jkparms(where=(Parameter ne "Scale")
                                      keep=Replicate Parameter Estimate Level1)
           ParameterEstimates=jkdf(where=(Parameter ne "Scale")
                                   keep= Replicate Parameter Level1 df)
           ConvergenceStatus=converge;
proc genmod data=jkweights;
  class sex race marital private education;
  weight repweight;
  model visits = sex race marital private education / dist=poisson;
  by replicate;
run;

/* End Section 3a */

/**********************************************************************************/
/*                 S A S  W E B  E X A M P L E  M A C R O                         */
/*                                                                                */
/*        NAME: JKLOOP                                                            */
/*       TITLE: Poisson Regressions for Complex Surveys                           */
/*        DATE: July 20, 2105                                                     */   
/*     PRODUCT: STAT                                                              */
/*       PROCS: MCMC                                                              */  
/*     SUPPORT: Allen McDowell                                                    */   
/*                                                                                */
/*   ARGUMENTS:                                                                   */
/*                                                                                */
/*   REPLICATES - the number of jackknife replicates                              */
/*                                                                                */
/* DESCRIPTION:                                                                   */
/* The JKLOOP macro loops through the R replicates and fits a Poisson model by    */
/* using the appropriate replicate sample and jackknife replicate weights. The    */
/* parameter estimates from each model are saved in the temporary data set Temp,  */ 
/* the degrees of freedom for each model in the temporary data set Temp2, and     */
/* the convergence status of the model is saved in the temporary data set Temp3.  */ 
/* A series of DATA steps then add the variable Replicate to Temp, Temp2, and     */
/* Temp3. The data sets Temp, Temp2, and Temp3 are then appended to the data      */
/* sets JKparms, JKDF, and Converge, respectively. Finally, the variable          */
/* Estimate in the data set FullSample is renamed as Estimate0.                   */                                         
/*                                                                                */  
/**********************************************************************************/

%macro jkloop(replicates=);
  %local _nopt; %let _nopt = %sysfunc(getoption(notes));
  options nonotes;
  ods select none;
  %do i=1 %to &replicates;
    ods output ParameterEstimates=temp(where=(Parameter ne "Scale")
                                        keep=Parameter Estimate Level1)
               ParameterEstimates=temp2(where=(Parameter ne "Scale")
                                    keep=Parameter Level1 df)
	          ConvergenceStatus=temp3;
    proc genmod data=jkweights;
      class sex race marital private education;
      weight RepWt_&i;
      model visits = sex race marital private education / dist=poisson;
    run;

    data temp;
      set temp;
      Replicate=&i;
    run;

    data temp2;
      set temp2;
      Replicate=&i;
    run;

    data temp3;
      set temp3;
      Replicate=&i;
    run;

    proc append base=jkparms data=temp;
    run;

    proc append base=jkdf data=temp2;
    run;

    proc append base=converge data=temp3;
    run;

  %end;

  data FullSample;
    set FullSample;
    rename estimate=estimate0;
  run;

  ods select all;
  options &_nopt;
%mend jkloop;

/* Before you use the macro JKLOOP, comment out the statements that are labeled Section 3a */
/* and remove the comment tags around Section 3b.                                          */

/* Begin Section 3b                                                                        */

/*
%jkloop(replicates=&replicates)
*/

/*  End Section 3b                                                                         */

ods output Statistics=statistics(keep=replicate sum);
proc surveymeans data=jkdf sum;
  var df;
  by replicate;
run;

data statistics;
  set statistics;
  full=ifn(sum=&p,0,1);
run;

data converge;
  merge converge statistics;
  by replicate;
run;

data converged;
  set converge(where=(Status=0 & full=0));
run;

data nobs;
  dsid=open("converged");
  converged_replicates=attrn(dsid, "nobs");
  call symput('R',converged_replicates);
run;

proc sort data=jkparms;
  by parameter level1;
run;

proc sort data=FullSample;
  by parameter level1;
run;

data jk;
  merge jkparms FullSample;
  by parameter level1;
run;

proc sort data=jk;
  by replicate parameter level1;
run;

data jk;
  merge jk jkcoefs converge;
  by replicate;
run;

data jkconverged;
  set jk(where=(Status=0 & full=0));
  sqrdev=JKCoefficient*(estimate-estimate0)**2;
run;

data vce;
  set jkconverged(keep= replicate parameter
                  level1 estimate estimate0);
  diff=estimate0-estimate;
run;

proc sort data=jkconverged;
  by parameter Level1;
run;

ods output Statistics=jkVariance;
proc surveymeans data=jkconverged sum plots=none;
  var sqrdev;
  by parameter Level1;
run;

data jkVariance(drop=stddev varname);
  merge jkVariance fullsample(rename=(estimate0=Estimate));
  by parameter Level1;
  StdErr=sqrt(Sum);
  rename Sum=Variance;
  DF=&R - &H;
  t=quantile('T', .975, &R-&H);
  ul=estimate+t*stderr;
  ll=estimate-t*stderr;
  label ul="Upper 95% CL";
  label ll="Lower 95% CL";
run;

ods select all;
title "Survey Poisson Regression";
title2 "with Delete-1 Jackknife Variance Estimation";
proc print data=jkVariance noobs label;
  var Parameter Level1 Estimate StdErr DF ll ul;
run;
title;title2;

proc sort data=jkdf;
  by replicate parameter level1;
run;

data temp;
  merge vce jkdf;
  by replicate parameter level1;
run;

proc transpose data=temp(where=(df=1)) out=temp2(drop=_name_) prefix=parm;
  by replicate;
  var diff;
run;

data temp3(drop=donorstratum);
  merge temp2 jkcoefs;
  by replicate;
  do i=1 to &p;
    row=i;
    output;
  end;
run;

data temp3(drop=parm: jkcoefficient i j);
  set temp3;
  array col[&p];
  array parm[*] parm:;
   do i=1 to &p;
     if row=i then do;
       do j = 1 to &p;
         col[j]=jkcoefficient*parm[i]*parm[j];
       end;
     end;
   end;
run;

proc sort data=temp3;
  by row;
run;

ods select none;
ods output Statistics=statistics(drop=StdDev);
proc surveymeans data=temp3 sum plots=none;
  var col1-col14;
  by row;
run;
ods select all;

proc transpose data=statistics out=CovB(drop=_name_ row) prefix=parm;
  var sum;
  by row;
run;

proc print data=covb noobs;
run;