R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(0 + ,0 + ,1 + ,40 + ,2898 + ,1 + ,0 + ,1 + ,26 + ,994 + ,0 + ,2 + ,0 + ,38 + ,3977 + ,0 + ,2 + ,1 + ,37 + ,3040 + ,0 + ,2 + ,0 + ,38 + ,3523 + ,0 + ,5 + ,1 + ,40 + ,3100 + ,1 + ,6 + ,0 + ,40 + ,3670 + ,1 + ,7 + ,0 + ,41 + ,3097 + ,1 + ,7 + ,1 + ,39 + ,3040 + ,1 + ,7 + ,1 + ,39 + ,3239 + ,0 + ,8 + ,0 + ,38 + ,2955 + ,0 + ,8 + ,0 + ,38 + ,2200 + ,1 + ,8 + ,0 + ,40 + ,3182 + ,0 + ,8 + ,0 + ,40 + ,3510 + ,1 + ,8 + ,1 + ,39 + ,3381 + ,1 + ,8 + ,0 + ,40 + ,3530 + ,1 + ,8 + ,1 + ,38 + ,2985 + ,1 + ,8 + ,0 + ,39 + ,3374 + ,0 + ,8 + ,0 + ,42 + ,3765 + ,1 + ,8 + ,0 + ,39 + ,2715 + ,1 + ,8 + ,0 + ,39 + ,3640 + ,1 + ,8 + ,0 + ,42 + ,3040 + ,1 + ,8 + ,0 + ,39 + ,3670 + ,1 + ,9 + ,1 + ,38 + ,2775 + ,1 + ,9 + ,0 + ,38 + ,2680 + ,1 + ,9 + ,0 + ,37 + ,2699 + ,1 + ,9 + ,1 + ,38 + ,2557 + ,1 + ,9 + ,0 + ,40 + ,3835 + ,0 + ,9 + ,1 + ,41 + ,3381 + ,1 + ,9 + ,0 + ,36 + ,2766 + ,1 + ,9 + ,0 + ,41 + ,3125 + ,1 + ,9 + ,1 + ,38 + ,2557 + ,1 + ,9 + ,0 + ,38 + ,2950 + ,0 + ,9 + ,0 + ,40 + ,3779 + ,0 + ,9 + ,0 + ,41 + ,2760 + ,1 + ,9 + ,0 + ,39 + ,3011 + ,1 + ,9 + ,1 + ,37 + ,2415 + ,1 + ,9 + ,0 + ,38 + ,3125 + ,0 + ,9 + ,1 + ,38 + ,3154 + ,0 + ,9 + ,0 + ,39 + ,3170 + ,0 + ,9 + ,0 + ,41 + ,2784 + ,1 + ,9 + ,1 + ,40 + ,2475 + ,1 + ,9 + ,0 + ,40 + ,2642 + ,0 + ,9 + ,1 + ,40 + ,4034 + ,0 + ,9 + ,0 + ,41 + ,3988 + ,1 + ,9 + ,0 + ,40 + ,3370 + ,0 + ,9 + ,0 + ,40 + ,2955 + ,1 + ,9 + ,0 + ,40 + ,2443 + ,0 + ,9 + ,0 + ,38 + ,4000 + ,1 + ,9 + ,0 + ,37 + ,3630 + ,1 + ,9 + ,0 + ,39 + ,2920 + ,1 + ,9 + ,1 + ,35 + ,2727 + ,0 + ,9 + ,0 + ,38 + ,3800 + ,0 + ,9 + ,0 + ,40 + ,3771 + ,1 + ,9 + ,0 + ,37 + ,3075 + ,1 + ,9 + ,1 + ,37 + ,2415 + ,1 + ,9 + ,0 + ,37 + ,2890 + ,0 + ,9 + ,0 + ,37 + ,3892 + ,1 + ,9 + ,1 + ,40 + ,3125 + ,0 + ,9 + ,1 + ,42 + ,3892 + ,0 + ,9 + ,0 + ,39 + ,3296 + ,0 + ,9 + ,0 + ,38 + ,3120 + ,0 + ,9 + ,1 + ,38 + ,1710 + ,1 + ,10 + ,0 + ,41 + ,3068 + ,1 + ,10 + ,0 + ,41 + ,3779 + ,1 + ,10 + ,1 + ,40 + ,2898 + ,0 + ,10 + ,0 + ,40 + ,3381 + ,1 + ,10 + ,1 + ,40 + ,3324 + ,1 + ,10 + ,1 + ,38 + ,2727 + ,0 + ,10 + ,1 + ,37 + ,2642 + ,1 + ,10 + ,1 + ,34 + ,2244 + ,0 + ,10 + ,0 + ,36 + ,2528 + ,0 + ,10 + ,1 + ,39 + ,3352 + ,0 + ,10 + ,1 + ,39 + ,3551 + ,0 + ,10 + ,0 + ,40 + ,3267 + ,1 + ,10 + ,0 + ,40 + ,3210 + ,1 + ,10 + ,0 + ,40 + ,3395 + ,0 + ,10 + ,1 + ,40 + ,3068 + ,1 + ,10 + ,1 + ,39 + ,3420 + ,1 + ,10 + ,0 + ,38 + ,2985 + ,1 + ,10 + ,0 + ,33 + ,2147 + ,1 + ,10 + ,0 + ,39 + ,2863 + ,1 + ,10 + ,0 + ,38 + ,3125 + ,0 + ,10 + ,0 + ,40 + ,4290 + ,1 + ,10 + ,1 + ,34 + ,2585 + ,1 + ,10 + ,0 + ,38 + ,2869 + ,1 + ,10 + ,1 + ,36 + ,3097 + ,1 + ,10 + ,0 + ,39 + ,2756 + ,0 + ,10 + ,1 + ,35 + ,3125 + ,0 + ,10 + ,1 + ,41 + ,3921 + ,1 + ,10 + ,0 + ,40 + ,3438 + ,1 + ,10 + ,0 + ,40 + ,2980 + ,1 + ,10 + ,0 + ,36 + ,2869 + ,0 + ,10 + ,1 + ,41 + ,4830 + ,0 + ,10 + ,1 + ,40 + ,3409 + ,1 + ,10 + ,0 + ,41 + ,3807 + ,1 + ,10 + ,1 + ,39 + ,3220 + ,0 + ,10 + ,1 + ,38 + ,2756 + ,0 + ,10 + ,1 + ,38 + ,3097 + ,1 + ,10 + ,0 + ,39 + ,3615 + ,0 + ,10 + ,0 + ,40 + ,2926 + ,1 + ,10 + ,0 + ,40 + ,4050 + ,1 + ,10 + ,0 + ,38 + ,2840 + ,1 + ,10 + ,0 + ,41 + ,2671 + ,1 + ,10 + ,0 + ,38 + ,3540 + ,1 + ,10 + ,0 + ,38 + ,2900 + ,0 + ,10 + ,0 + ,40 + ,3779 + ,1 + ,10 + ,0 + ,38 + ,2740 + ,1 + ,10 + ,0 + ,35 + ,2320 + ,1 + ,10 + ,0 + ,38 + ,2905 + ,1 + ,10 + ,1 + ,39 + ,3070 + ,0 + ,10 + ,1 + ,41 + ,3352 + ,1 + ,10 + ,0 + ,36 + ,2581 + ,1 + ,10 + ,0 + ,41 + ,4091 + ,1 + ,10 + ,1 + ,43 + ,3235 + ,1 + ,10 + ,0 + ,41 + ,3722 + ,1 + ,10 + ,1 + ,40 + ,3438 + ,0 + ,10 + ,0 + ,40 + ,3475 + ,1 + ,10 + ,0 + ,31 + ,3438 + ,1 + ,10 + ,0 + ,39 + ,3409 + ,0 + ,10 + ,0 + ,38 + ,3040 + ,0 + ,10 + ,0 + ,38 + ,3800 + ,1 + ,10 + ,1 + ,36 + ,3182 + ,1 + ,10 + ,1 + ,39 + ,3494 + ,1 + ,10 + ,1 + ,35 + ,2131 + ,1 + ,10 + ,0 + ,35 + ,1790 + ,0 + ,10 + ,0 + ,40 + ,3864 + ,1 + ,10 + ,0 + ,39 + ,2955 + ,0 + ,10 + ,0 + ,40 + ,3949 + ,1 + ,10 + ,0 + ,35 + ,2273 + ,1 + ,10 + ,0 + ,39 + ,2620 + ,1 + ,10 + ,0 + ,36 + ,2205 + ,1 + ,10 + ,0 + ,41 + ,2860 + ,1 + ,10 + ,0 + ,40 + ,3381 + ,1 + ,10 + ,0 + ,40 + ,3190 + ,1 + ,10 + ,0 + ,40 + ,3267 + ,0 + ,10 + ,1 + ,40 + ,3665 + ,1 + ,10 + ,0 + ,38 + ,3160 + ,1 + ,10 + ,0 + ,39 + ,3216 + ,1 + ,10 + ,1 + ,37 + ,2415 + ,1 + ,10 + ,0 + ,40 + ,3200 + ,0 + ,10 + ,0 + ,40 + ,3352 + ,1 + ,10 + ,1 + ,40 + ,3470 + ,1 + ,10 + ,0 + ,38 + ,3580 + ,0 + ,10 + ,0 + ,41 + ,3750 + ,1 + ,10 + ,0 + ,42 + ,3035 + ,1 + ,10 + ,0 + ,38 + ,3438 + ,1 + ,10 + ,0 + ,38 + ,3068 + ,0 + ,10 + ,0 + ,39 + ,3438 + ,1 + ,10 + ,1 + ,38 + ,3750 + ,1 + ,10 + ,1 + ,42 + ,3409 + ,1 + ,10 + ,1 + ,39 + ,1790 + ,1 + ,10 + ,0 + ,40 + ,3170 + ,0 + ,10 + ,0 + ,38 + ,3800 + ,1 + ,10 + ,1 + ,40 + ,2825 + ,1 + ,10 + ,0 + ,37 + ,2756 + ,1 + ,10 + ,0 + ,40 + ,3040 + ,0 + ,10 + ,1 + ,39 + ,2585 + ,1 + ,10 + ,1 + ,38 + ,2756 + ,0 + ,10 + ,1 + ,41 + ,3523 + ,1 + ,10 + ,1 + ,26 + ,909 + ,0 + ,10 + ,1 + ,40 + ,3779 + ,1 + ,10 + ,1 + ,38 + ,2415 + ,0 + ,10 + ,0 + ,40 + ,3665 + ,1 + ,10 + ,0 + ,39 + ,3295 + ,1 + ,10 + ,1 + ,39 + ,2699 + ,0 + ,10 + ,1 + ,39 + ,3267 + ,1 + ,10 + ,0 + ,39 + ,3460 + ,1 + ,11 + ,0 + ,28 + ,1250 + ,1 + ,11 + ,0 + ,38 + ,3730 + ,1 + ,11 + ,1 + ,40 + ,3523 + ,1 + ,11 + ,1 + ,40 + ,3011 + ,1 + ,11 + ,1 + ,39 + ,3740 + ,1 + ,11 + ,1 + ,43 + ,3949 + ,0 + ,11 + ,1 + ,37 + ,2330 + ,0 + ,11 + ,0 + ,40 + ,3182 + ,1 + ,11 + ,0 + ,39 + ,3040 + ,1 + ,11 + ,0 + ,37 + ,2557 + ,0 + ,11 + ,0 + ,38 + ,3147 + ,1 + ,11 + ,1 + ,36 + ,2671 + ,1 + ,11 + ,1 + ,40 + ,2557 + ,1 + ,11 + ,0 + ,38 + ,3200 + ,1 + ,11 + ,0 + ,41 + ,3203 + ,1 + ,11 + ,1 + ,36 + ,2756 + ,1 + ,11 + ,0 + ,32 + ,1619 + ,1 + ,11 + ,1 + ,35 + ,2330 + ,1 + ,11 + ,1 + ,39 + ,3381 + ,0 + ,11 + ,1 + ,38 + ,2983 + ,0 + ,11 + ,1 + ,38 + ,2600 + ,1 + ,11 + ,0 + ,40 + ,3585 + ,0 + ,11 + ,0 + ,40 + ,3693 + ,1 + ,11 + ,0 + ,38 + ,3020 + ,1 + ,11 + ,1 + ,34 + ,1800 + ,1 + ,11 + ,0 + ,38 + ,2290 + ,1 + ,11 + ,0 + ,40 + ,3750 + ,0 + ,11 + ,0 + ,40 + ,3551 + ,0 + ,11 + ,0 + ,37 + ,3779 + ,0 + ,11 + ,0 + ,38 + ,3125 + ,1 + ,11 + ,0 + ,38 + ,3470 + ,1 + ,11 + ,0 + ,41 + ,3210 + ,1 + ,11 + ,1 + ,40 + ,3100 + ,1 + ,11 + ,1 + ,38 + ,3182 + ,1 + ,11 + ,0 + ,37 + ,3355 + ,1 + ,11 + ,0 + ,42 + ,2750 + ,1 + ,11 + ,0 + ,39 + ,4006 + ,1 + ,11 + ,0 + ,38 + ,3438 + ,0 + ,11 + ,0 + ,40 + ,3068 + ,1 + ,11 + ,1 + ,40 + ,2650 + ,1 + ,11 + ,0 + ,39 + ,3005 + ,1 + ,11 + ,0 + ,39 + ,2993 + ,1 + ,11 + ,0 + ,40 + ,3494 + ,1 + ,11 + ,0 + ,41 + ,2795 + ,1 + ,11 + ,0 + ,41 + ,3260 + ,1 + ,11 + ,1 + ,41 + ,2959 + ,1 + ,11 + ,0 + ,38 + ,4020 + ,0 + ,11 + ,0 + ,37 + ,2949 + ,1 + ,11 + ,0 + ,38 + ,3900 + ,1 + ,11 + ,1 + ,39 + ,2926 + ,1 + ,11 + ,1 + ,38 + ,3239 + ,1 + ,11 + ,0 + ,41 + ,3409 + ,1 + ,11 + ,0 + ,42 + ,3630 + ,1 + ,11 + ,0 + ,41 + ,3466 + ,1 + ,11 + ,0 + ,40 + ,2330 + ,0 + ,11 + ,1 + ,40 + ,3949 + ,1 + ,11 + ,1 + ,38 + ,2813 + ,1 + ,11 + ,0 + ,37 + ,3137 + ,0 + ,11 + ,0 + ,40 + ,3296 + ,1 + ,11 + ,0 + ,41 + ,3267 + ,1 + ,11 + ,0 + ,38 + ,2940 + ,1 + ,11 + ,0 + ,40 + ,3323 + ,1 + ,11 + ,1 + ,37 + ,3409 + ,1 + ,11 + ,1 + ,41 + ,3296 + ,0 + ,11 + ,1 + ,40 + ,3807 + ,1 + ,11 + ,0 + ,38 + ,3296 + ,1 + ,11 + ,1 + ,38 + ,2983 + ,1 + ,11 + ,1 + ,30 + ,1680 + ,1 + ,11 + ,0 + ,39 + ,3551 + ,1 + ,11 + ,0 + ,37 + ,2385 + ,1 + ,11 + ,0 + ,40 + ,3356 + ,1 + ,11 + ,0 + ,40 + ,3466 + ,0 + ,11 + ,0 + ,41 + ,2699 + ,0 + ,11 + ,0 + ,40 + ,4233 + ,1 + ,11 + ,1 + ,40 + ,2890 + ,1 + ,11 + ,1 + ,37 + ,2443 + ,0 + ,11 + ,1 + ,39 + ,3580 + ,0 + ,11 + ,1 + ,42 + ,3551 + ,1 + ,11 + ,1 + ,31 + ,2030 + ,1 + ,11 + ,0 + ,36 + ,3125 + ,1 + ,11 + ,0 + ,37 + ,2700 + ,1 + ,11 + ,1 + ,33 + ,2216 + ,0 + ,11 + ,1 + ,40 + ,3182 + ,1 + ,11 + ,1 + ,39 + ,3040 + ,1 + ,11 + ,1 + ,39 + ,3011 + ,1 + ,11 + ,1 + ,39 + ,2800 + ,1 + ,11 + ,0 + ,38 + ,3425 + ,1 + ,11 + ,1 + ,37 + ,2760 + ,1 + ,11 + ,0 + ,41 + ,3665 + ,1 + ,11 + ,0 + ,36 + ,1733 + ,1 + ,11 + ,0 + ,40 + ,2850 + ,1 + ,11 + ,1 + ,38 + ,2485 + ,1 + ,11 + ,0 + ,40 + ,3381 + ,1 + ,11 + ,0 + ,39 + ,3065 + ,1 + ,11 + ,0 + ,39 + ,3167 + ,1 + ,11 + ,0 + ,38 + ,2635 + ,0 + ,11 + ,1 + ,40 + ,3665 + ,1 + ,11 + ,1 + ,39 + ,3300 + ,0 + ,11 + ,0 + ,40 + ,3636 + ,1 + ,11 + ,1 + ,38 + ,2813 + ,1 + ,11 + ,1 + ,42 + ,3920 + ,1 + ,11 + ,0 + ,39 + ,2671 + ,1 + ,11 + ,0 + ,40 + ,3182 + ,1 + ,11 + ,1 + ,40 + ,3409 + ,1 + ,11 + ,0 + ,37 + ,3600 + ,1 + ,11 + ,0 + ,40 + ,2557 + ,1 + ,11 + ,0 + ,40 + ,3890 + ,1 + ,11 + ,0 + ,39 + ,4148 + ,1 + ,11 + ,1 + ,32 + ,2030 + ,1 + ,11 + ,1 + ,39 + ,2890 + ,1 + ,11 + ,1 + ,36 + ,2373 + ,1 + ,11 + ,0 + ,41 + ,2977 + ,1 + ,11 + ,1 + ,38 + ,2840 + ,1 + ,11 + ,0 + ,41 + ,2658 + ,0 + ,11 + ,1 + ,40 + ,2983 + ,1 + ,11 + ,0 + ,39 + ,3700 + ,0 + ,12 + ,1 + ,40 + ,3608 + ,0 + ,12 + ,0 + ,40 + ,3559 + ,1 + ,12 + ,0 + ,40 + ,3540 + ,1 + ,12 + ,1 + ,39 + ,3440 + ,1 + ,12 + ,0 + ,39 + ,3324 + ,1 + ,12 + ,0 + ,38 + ,3494 + ,1 + ,12 + ,0 + ,41 + ,3636 + ,1 + ,12 + ,1 + ,39 + ,2813 + ,1 + ,12 + ,0 + ,39 + ,3239 + ,0 + ,12 + ,1 + ,39 + ,3068 + ,1 + ,12 + ,0 + ,41 + ,3210 + ,0 + ,12 + ,0 + ,42 + ,3892 + ,0 + ,12 + ,0 + ,39 + ,3864 + ,1 + ,12 + ,0 + ,40 + ,2784 + ,0 + ,12 + ,0 + ,40 + ,3551 + ,1 + ,12 + ,0 + ,40 + ,3608 + ,1 + ,12 + ,0 + ,39 + ,3130 + ,0 + ,12 + ,1 + ,40 + ,3466 + ,1 + ,12 + ,0 + ,40 + ,3608 + ,1 + ,12 + ,0 + ,37 + ,2225 + ,1 + ,12 + ,0 + ,36 + ,2585 + ,1 + ,12 + ,0 + ,40 + ,2917 + ,1 + ,12 + ,0 + ,36 + ,2415 + ,0 + ,12 + ,0 + ,40 + ,4489 + ,0 + ,12 + ,0 + ,40 + ,3977 + ,1 + ,12 + ,0 + ,38 + ,3080 + ,1 + ,12 + ,0 + ,40 + ,3722 + ,0 + ,12 + ,1 + ,40 + ,3097 + ,1 + ,12 + ,0 + ,40 + ,2960 + ,0 + ,12 + ,0 + ,39 + ,3807 + ,0 + ,12 + ,0 + ,40 + ,3210 + ,1 + ,12 + ,0 + ,40 + ,3210 + ,1 + ,12 + ,1 + ,20 + ,284 + ,1 + ,12 + ,0 + ,40 + ,3805 + ,0 + ,12 + ,1 + ,41 + ,3267 + ,0 + ,12 + ,1 + ,41 + ,3722 + ,0 + ,12 + ,0 + ,40 + ,3864 + ,0 + ,12 + ,1 + ,42 + ,4318 + ,0 + ,12 + ,1 + ,42 + ,3807 + ,1 + ,12 + ,1 + ,22 + ,455 + ,0 + ,12 + ,0 + ,40 + ,3722 + ,1 + ,12 + ,0 + ,38 + ,3324 + ,1 + ,12 + ,0 + ,40 + ,3665 + ,1 + ,12 + ,0 + ,41 + ,3240 + ,1 + ,12 + ,0 + ,40 + ,3182 + ,1 + ,12 + ,1 + ,40 + ,3040 + ,0 + ,12 + ,1 + ,40 + ,3665 + ,0 + ,12 + ,0 + ,40 + ,4148 + ,0 + ,12 + ,0 + ,38 + ,3636 + ,0 + ,12 + ,1 + ,33 + ,1591 + ,0 + ,12 + ,0 + ,39 + ,4135 + ,0 + ,12 + ,0 + ,40 + ,3352 + ,1 + ,12 + ,1 + ,38 + ,3011 + ,1 + ,12 + ,1 + ,34 + ,2301 + ,1 + ,12 + ,0 + ,40 + ,2740 + ,1 + ,12 + ,0 + ,41 + ,2450 + ,0 + ,12 + ,1 + ,40 + ,2472 + ,0 + ,12 + ,0 + ,40 + ,4148 + ,0 + ,12 + ,0 + ,40 + ,3352 + ,0 + ,12 + ,0 + ,38 + ,4034 + ,0 + ,12 + ,0 + ,42 + ,3807 + ,0 + ,12 + ,0 + ,39 + ,2926 + ,1 + ,12 + ,0 + ,40 + ,2699 + ,1 + ,12 + ,0 + ,39 + ,2990 + ,0 + ,12 + ,0 + ,38 + ,2898 + ,0 + ,12 + ,1 + ,41 + ,3864 + ,1 + ,12 + ,0 + ,40 + ,2970 + ,1 + ,12 + ,0 + ,42 + ,3655 + ,1 + ,12 + ,0 + ,35 + ,2890 + ,1 + ,12 + ,1 + ,38 + ,3980 + ,0 + ,12 + ,1 + ,40 + ,2756 + ,0 + ,12 + ,0 + ,40 + ,4148 + ,0 + ,12 + ,0 + ,40 + ,3494 + ,1 + ,12 + ,1 + ,37 + ,2830 + ,1 + ,12 + ,0 + ,39 + ,2983 + ,1 + ,12 + ,0 + ,40 + ,3430 + ,1 + ,12 + ,0 + ,35 + ,2159 + ,0 + ,12 + ,1 + ,42 + ,3125 + ,0 + ,12 + ,0 + ,40 + ,3239 + ,0 + ,12 + ,0 + ,38 + ,4262 + ,1 + ,12 + ,0 + ,38 + ,3459 + ,1 + ,12 + ,0 + ,28 + ,1190 + ,0 + ,12 + ,0 + ,40 + ,3182 + ,1 + ,12 + ,0 + ,35 + ,3352 + ,0 + ,12 + ,0 + ,39 + ,2898 + ,0 + ,12 + ,0 + ,40 + ,3154 + ,1 + ,12 + ,1 + ,28 + ,1165 + ,0 + ,12 + ,1 + ,39 + ,2102 + ,1 + ,12 + ,0 + ,40 + ,3110 + ,0 + ,12 + ,1 + ,34 + ,2367 + ,1 + ,12 + ,0 + ,39 + ,3892 + ,1 + ,12 + ,1 + ,40 + ,3011 + ,1 + ,12 + ,0 + ,37 + ,3125 + ,1 + ,12 + ,0 + ,41 + ,3490 + ,1 + ,12 + ,1 + ,38 + ,2955 + ,1 + ,12 + ,0 + ,41 + ,3977 + ,0 + ,12 + ,0 + ,42 + ,3665 + ,1 + ,12 + ,0 + ,39 + ,3355 + ,1 + ,12 + ,0 + ,40 + ,3182 + ,0 + ,12 + ,0 + ,40 + ,3239 + ,1 + ,12 + ,0 + ,40 + ,2955 + ,1 + ,12 + ,0 + ,39 + ,3182 + ,1 + ,12 + ,0 + ,38 + ,2784 + ,0 + ,12 + ,0 + ,40 + ,3802 + ,1 + ,12 + ,1 + ,40 + ,3270 + ,1 + ,12 + ,0 + ,39 + ,2225 + ,0 + ,12 + ,1 + ,36 + ,3182 + ,1 + ,12 + ,0 + ,41 + ,3210 + ,0 + ,12 + ,0 + ,39 + ,3779 + ,1 + ,12 + ,0 + ,36 + ,2415 + ,0 + ,12 + ,0 + ,40 + ,4675 + ,1 + ,12 + ,0 + ,38 + ,3352 + ,1 + ,12 + ,1 + ,38 + ,2900 + ,0 + ,12 + ,0 + ,41 + ,3864 + ,1 + ,12 + ,0 + ,40 + ,2841 + ,0 + ,12 + ,0 + ,40 + ,2557 + ,0 + ,12 + ,0 + ,40 + ,3494 + ,1 + ,12 + ,0 + ,40 + ,3438 + ,1 + ,12 + ,0 + ,36 + ,2445 + ,1 + ,12 + ,0 + ,42 + ,3714 + ,1 + ,12 + ,0 + ,42 + ,2935 + ,0 + ,12 + ,0 + ,41 + ,3381 + ,0 + ,12 + ,0 + ,41 + ,3864 + ,1 + ,12 + ,1 + ,40 + ,3240 + ,1 + ,12 + ,0 + ,40 + ,2980 + ,1 + ,12 + ,1 + ,38 + ,2710 + ,1 + ,12 + ,0 + ,40 + ,3570 + ,1 + ,12 + ,0 + ,42 + ,3864 + ,0 + ,12 + ,1 + ,39 + ,2500 + ,1 + ,12 + ,0 + ,28 + ,1160 + ,1 + ,12 + ,0 + ,39 + ,2443 + ,1 + ,12 + ,0 + ,43 + ,3040 + ,0 + ,12 + ,0 + ,41 + ,3864 + ,1 + ,12 + ,0 + ,42 + ,3205 + ,0 + ,12 + ,0 + ,33 + ,1875 + ,1 + ,12 + ,0 + ,39 + ,3835 + ,0 + ,12 + ,1 + ,39 + ,2841 + ,1 + ,12 + ,0 + ,36 + ,2926 + ,1 + ,12 + ,0 + ,40 + ,2983 + ,0 + ,12 + ,1 + ,40 + ,3267 + ,1 + ,12 + ,1 + ,38 + ,2890 + ,1 + ,12 + ,0 + ,38 + ,2443 + ,0 + ,12 + ,0 + ,38 + ,3352 + ,0 + ,12 + ,0 + ,39 + ,3438 + ,1 + ,12 + ,1 + ,38 + ,3120 + ,0 + ,12 + ,0 + ,40 + ,4176 + ,0 + ,12 + ,0 + ,42 + ,4148 + ,1 + ,12 + ,0 + ,40 + ,3270 + ,0 + ,12 + ,0 + ,40 + ,3523 + ,0 + ,12 + ,0 + ,42 + ,3020 + ,0 + ,12 + ,0 + ,38 + ,3620 + ,1 + ,12 + ,0 + ,42 + ,3300 + ,0 + ,12 + ,0 + ,37 + ,3779 + ,0 + ,12 + ,1 + ,41 + ,3665 + ,0 + ,12 + ,1 + ,40 + ,2869 + ,1 + ,12 + ,0 + ,38 + ,3154 + ,0 + ,12 + ,1 + ,39 + ,2810 + ,1 + ,12 + ,1 + ,34 + ,2131 + ,0 + ,12 + ,0 + ,40 + ,4233 + ,0 + ,12 + ,0 + ,38 + ,3860 + ,1 + ,12 + ,0 + ,40 + ,3920 + ,0 + ,12 + ,0 + ,42 + ,3693 + ,1 + ,12 + ,0 + ,40 + ,3580 + ,1 + ,12 + ,0 + ,40 + ,3608 + ,1 + ,12 + ,1 + ,37 + ,2330 + ,0 + ,12 + ,1 + ,40 + ,2813 + ,1 + ,12 + ,0 + ,37 + ,3068 + ,1 + ,12 + ,0 + ,39 + ,3250 + ,1 + ,12 + ,0 + ,39 + ,3175 + ,1 + ,12 + ,0 + ,39 + ,2940 + ,0 + ,12 + ,0 + ,41 + ,3267 + ,1 + ,12 + ,0 + ,40 + ,3645 + ,0 + ,12 + ,0 + ,41 + ,4205 + ,0 + ,12 + ,0 + ,34 + ,1818 + ,0 + ,12 + ,0 + ,35 + ,2983 + ,1 + ,12 + ,0 + ,40 + ,2756 + ,0 + ,12 + ,0 + ,39 + ,2642 + ,0 + ,12 + ,0 + ,41 + ,3182 + ,0 + ,12 + ,0 + ,41 + ,3125 + ,0 + ,12 + ,0 + ,38 + ,3949 + ,0 + ,12 + ,0 + ,39 + ,3466 + ,1 + ,12 + ,1 + ,38 + ,3250 + ,0 + ,12 + ,0 + ,40 + ,3750 + ,0 + ,12 + ,0 + ,40 + ,3466 + ,0 + ,12 + ,1 + ,40 + ,2557 + ,0 + ,12 + ,0 + ,40 + ,3722 + ,0 + ,12 + ,0 + ,40 + ,3551 + ,1 + ,12 + ,0 + ,39 + ,3097 + ,1 + ,12 + ,0 + ,35 + ,1763 + ,1 + ,12 + ,0 + ,40 + ,3610 + ,1 + ,12 + ,0 + ,39 + ,3551 + ,0 + ,12 + ,1 + ,40 + ,2460 + ,1 + ,12 + ,0 + ,36 + ,2358 + ,0 + ,12 + ,0 + ,40 + ,3779 + ,1 + ,12 + ,1 + ,40 + ,2671 + ,1 + ,12 + ,0 + ,38 + ,3210 + ,0 + ,12 + ,0 + ,40 + ,4375 + ,1 + ,12 + ,0 + ,40 + ,4160 + ,1 + ,12 + ,0 + ,40 + ,2841 + ,0 + ,12 + ,0 + ,40 + ,3352 + ,1 + ,12 + ,1 + ,38 + ,2790 + ,1 + ,12 + ,0 + ,42 + ,2540 + ,1 + ,12 + ,0 + ,39 + ,2890 + ,1 + ,12 + ,1 + ,40 + ,3182 + ,0 + ,12 + ,0 + ,38 + ,2557 + ,0 + ,12 + ,0 + ,40 + ,3551 + ,1 + ,12 + ,1 + ,40 + ,3324 + ,1 + ,12 + ,0 + ,39 + ,2620 + ,1 + ,12 + ,0 + ,38 + ,3068 + ,1 + ,12 + ,0 + ,38 + ,3150 + ,1 + ,12 + ,0 + ,40 + ,3636 + ,1 + ,12 + ,0 + ,40 + ,3130 + ,0 + ,12 + ,0 + ,35 + ,2784 + ,1 + ,12 + ,0 + ,40 + ,3520 + ,0 + ,12 + ,0 + ,40 + ,4205 + ,1 + ,12 + ,0 + ,38 + ,3210 + ,0 + ,12 + ,0 + ,42 + ,3523 + ,1 + ,12 + ,0 + ,38 + ,3125 + ,1 + ,12 + ,1 + ,39 + ,3097 + ,1 + ,12 + ,1 + ,39 + ,3125 + ,1 + ,12 + ,0 + ,40 + ,2930 + ,0 + ,12 + ,0 + ,39 + ,3665 + ,0 + ,12 + ,0 + ,40 + ,3693 + ,1 + ,12 + ,0 + ,41 + ,3835 + ,1 + ,12 + ,0 + ,40 + ,3636 + ,1 + ,12 + ,0 + ,43 + ,4165 + ,1 + ,12 + ,1 + ,38 + ,3035 + ,1 + ,12 + ,0 + ,37 + ,2000 + ,1 + ,12 + ,1 + ,39 + ,2740 + ,0 + ,12 + ,0 + ,38 + ,3608 + ,0 + ,12 + ,0 + ,39 + ,3750 + ,0 + ,12 + ,0 + ,40 + ,4091 + ,1 + ,12 + ,0 + ,40 + ,3267 + ,1 + ,12 + ,0 + ,41 + ,3097 + ,1 + ,12 + ,1 + ,40 + ,3381 + ,1 + ,12 + ,0 + ,39 + ,3040 + ,0 + ,12 + ,0 + ,42 + ,3779 + ,0 + ,12 + ,1 + ,40 + ,3523 + ,0 + ,12 + ,0 + ,40 + ,3551 + ,0 + ,12 + ,0 + ,40 + ,3636 + ,0 + ,12 + ,1 + ,38 + ,3154 + ,1 + ,12 + ,1 + ,38 + ,2841 + ,1 + ,12 + ,0 + ,37 + ,2960 + ,1 + ,12 + ,0 + ,40 + ,3097 + ,1 + ,12 + ,0 + ,39 + ,3600 + ,0 + ,12 + ,1 + ,38 + ,3409 + ,1 + ,12 + ,0 + ,39 + ,2880 + ,1 + ,12 + ,0 + ,39 + ,2690 + ,0 + ,12 + ,0 + ,41 + ,3090 + ,0 + ,12 + ,0 + ,40 + ,3409 + ,1 + ,12 + ,1 + ,32 + ,2159 + ,1 + ,12 + ,0 + ,40 + ,3921 + ,0 + ,12 + ,0 + ,40 + ,2940 + ,1 + ,12 + ,0 + ,40 + ,3400 + ,1 + ,12 + ,0 + ,40 + ,3239 + ,0 + ,12 + ,0 + ,40 + ,3636 + ,1 + ,12 + ,0 + ,40 + ,4347 + ,1 + ,12 + ,1 + ,38 + ,2977 + ,1 + ,12 + ,0 + ,42 + ,3636 + ,0 + ,12 + ,0 + ,39 + ,2671 + ,0 + ,12 + ,0 + ,40 + ,3949 + ,0 + ,12 + ,0 + ,37 + ,3608 + ,1 + ,12 + ,0 + ,40 + ,3710 + ,0 + ,12 + ,1 + ,35 + ,2216 + ,1 + ,12 + ,0 + ,36 + ,2500 + ,1 + ,12 + ,0 + ,29 + ,1450 + ,0 + ,12 + ,0 + ,40 + ,3296 + ,1 + ,12 + ,0 + ,38 + ,3510 + ,1 + ,12 + ,0 + ,41 + ,3950 + ,0 + ,12 + ,0 + ,39 + ,4176 + ,1 + ,12 + ,1 + ,40 + ,3636 + ,1 + ,12 + ,0 + ,40 + ,3125 + ,0 + ,12 + ,0 + ,38 + ,3845 + ,1 + ,12 + ,1 + ,41 + ,3523 + ,0 + ,12 + ,0 + ,39 + ,2784 + ,0 + ,12 + ,0 + ,40 + ,3494 + ,0 + ,12 + ,0 + ,42 + ,4233 + ,0 + ,12 + ,0 + ,40 + ,3779 + ,1 + ,12 + ,0 + ,39 + ,3720 + ,1 + ,12 + ,0 + ,38 + ,3302 + ,0 + ,12 + ,0 + ,41 + ,3693 + ,1 + ,12 + ,0 + ,35 + ,2580 + ,1 + ,12 + ,0 + ,38 + ,2597 + ,0 + ,12 + ,0 + ,38 + ,2926 + ,0 + ,12 + ,0 + ,41 + ,4505 + ,0 + ,12 + ,0 + ,41 + ,3977 + ,1 + ,12 + ,0 + ,41 + ,3182 + ,1 + ,12 + ,1 + ,36 + ,2160 + ,1 + ,12 + ,0 + ,39 + ,3402 + ,0 + ,12 + ,0 + ,40 + ,3239 + ,0 + ,12 + ,0 + ,39 + ,3608 + ,0 + ,12 + ,0 + ,38 + ,3267 + ,0 + ,12 + ,0 + ,39 + ,3995 + ,1 + ,12 + ,1 + ,39 + ,3068 + ,0 + ,12 + ,0 + ,40 + ,4065 + ,1 + ,12 + ,0 + ,38 + ,3346 + ,0 + ,12 + ,0 + ,40 + ,3154 + ,1 + ,12 + ,1 + ,33 + ,2515 + ,0 + ,12 + ,0 + ,39 + ,3182 + ,0 + ,12 + ,1 + ,36 + ,2443 + ,0 + ,12 + ,0 + ,40 + ,3807 + ,1 + ,12 + ,1 + ,38 + ,3130 + ,0 + ,12 + ,1 + ,40 + ,3722 + ,1 + ,12 + ,0 + ,40 + ,2727 + ,0 + ,12 + ,1 + ,40 + ,3011 + ,1 + ,12 + ,0 + ,29 + ,793 + ,1 + ,12 + ,0 + ,41 + ,3494 + ,1 + ,12 + ,0 + ,35 + ,2386 + ,0 + ,12 + ,0 + ,40 + ,3466 + ,0 + ,12 + ,0 + ,38 + ,2760 + ,0 + ,12 + ,0 + ,41 + ,4659 + ,0 + ,12 + ,0 + ,39 + ,3608 + ,1 + ,12 + ,0 + ,38 + ,2557 + ,0 + ,12 + ,0 + ,40 + ,2671 + ,0 + ,12 + ,0 + ,38 + ,3740 + ,1 + ,12 + ,1 + ,33 + ,2614 + ,0 + ,12 + ,0 + ,40 + ,3949 + ,1 + ,12 + ,0 + ,39 + ,3630 + ,0 + ,12 + ,0 + ,39 + ,2898 + ,1 + ,12 + ,0 + ,40 + ,3940 + ,1 + ,12 + ,0 + ,23 + ,701 + ,0 + ,12 + ,0 + ,40 + ,3977 + ,1 + ,12 + ,1 + ,34 + ,2472 + ,1 + ,12 + ,1 + ,37 + ,3011 + ,1 + ,12 + ,0 + ,39 + ,3352 + ,0 + ,12 + ,0 + ,41 + ,3494 + ,0 + ,12 + ,0 + ,40 + ,4091 + ,1 + ,12 + ,0 + ,38 + ,3530 + ,1 + ,12 + ,1 + ,36 + ,3010 + ,0 + ,12 + ,0 + ,40 + ,3068 + ,1 + ,12 + ,0 + ,42 + ,3636 + ,1 + ,12 + ,1 + ,41 + ,3665 + ,1 + ,12 + ,0 + ,42 + ,3430 + ,1 + ,12 + ,0 + ,38 + ,3182 + ,1 + ,12 + ,0 + ,39 + ,3300 + ,1 + ,12 + ,0 + ,37 + ,2557 + ,1 + ,12 + ,0 + ,37 + ,2550 + ,0 + ,12 + ,1 + ,40 + ,3097 + ,0 + ,12 + ,0 + ,40 + ,3154 + ,1 + ,12 + ,0 + ,40 + ,3625 + ,1 + ,12 + ,1 + ,40 + ,3330 + ,0 + ,12 + ,0 + ,40 + ,4347 + ,0 + ,12 + ,0 + ,39 + ,3296 + ,1 + ,12 + ,0 + ,39 + ,3154 + ,1 + ,12 + ,0 + ,41 + ,4050 + ,0 + ,12 + ,0 + ,36 + ,3239 + ,0 + ,12 + ,0 + ,40 + ,3608 + ,1 + ,12 + ,0 + ,36 + ,2500 + ,0 + ,12 + ,0 + ,41 + ,3438 + ,0 + ,12 + ,0 + ,40 + ,3977 + ,1 + ,12 + ,1 + ,38 + ,2780 + ,1 + ,12 + ,1 + ,36 + ,2189 + ,1 + ,12 + ,1 + ,34 + ,1760 + ,1 + ,12 + ,0 + ,35 + ,2360 + ,0 + ,12 + ,1 + ,40 + ,3608 + ,1 + ,12 + ,0 + ,40 + ,3892 + ,1 + ,12 + ,0 + ,39 + ,2870 + ,1 + ,12 + ,1 + ,38 + ,2926 + ,0 + ,12 + ,0 + ,40 + ,3905 + ,1 + ,12 + ,0 + ,39 + ,3020 + ,0 + ,12 + ,0 + ,39 + ,3779 + ,0 + ,12 + ,0 + ,40 + ,3381 + ,1 + ,12 + ,0 + ,39 + ,3280 + ,1 + ,12 + ,0 + ,41 + ,3154 + ,0 + ,12 + ,0 + ,40 + ,3494 + ,1 + ,12 + ,1 + ,40 + ,3523 + ,0 + ,12 + ,0 + ,32 + ,1861 + ,1 + ,12 + ,1 + ,40 + ,3360 + ,0 + ,12 + ,0 + ,40 + ,3977 + ,1 + ,12 + ,1 + ,38 + ,3296 + ,0 + ,12 + ,0 + ,38 + ,3523 + ,1 + ,12 + ,0 + ,40 + ,3693 + ,1 + ,12 + ,0 + ,39 + ,3921 + ,1 + ,12 + ,1 + ,26 + ,653 + ,1 + ,12 + ,1 + ,40 + ,3494 + ,1 + ,12 + ,0 + ,41 + ,3500 + ,1 + ,12 + ,0 + ,38 + ,2863 + ,0 + ,12 + ,0 + ,40 + ,3494 + ,1 + ,12 + ,0 + ,38 + ,3335 + ,0 + ,12 + ,0 + ,40 + ,3750 + ,0 + ,12 + ,0 + ,42 + ,3636 + ,1 + ,12 + ,1 + ,39 + ,4190 + ,1 + ,12 + ,0 + ,39 + ,3540 + ,1 + ,12 + ,1 + ,40 + ,3381 + ,1 + ,12 + ,0 + ,39 + ,3580 + ,1 + ,12 + ,0 + ,32 + ,2813 + ,0 + ,12 + ,0 + ,41 + ,3494 + ,1 + ,12 + ,0 + ,42 + ,4065 + ,0 + ,12 + ,0 + ,38 + ,3520 + ,1 + ,12 + ,0 + ,39 + ,3551 + ,0 + ,12 + ,0 + ,40 + ,4063 + ,1 + ,12 + ,0 + ,39 + ,3210 + ,1 + ,12 + ,0 + ,40 + ,3640 + ,0 + ,12 + ,1 + ,40 + ,2770 + ,1 + ,12 + ,1 + ,41 + ,3750 + ,1 + ,12 + ,0 + ,41 + ,3695 + ,1 + ,12 + ,0 + ,41 + ,3305 + ,0 + ,12 + ,0 + ,40 + ,3892 + ,0 + ,12 + ,1 + ,42 + ,3210 + ,1 + ,12 + ,1 + ,40 + ,2851 + ,1 + ,12 + ,0 + ,38 + ,3175 + ,1 + ,12 + ,0 + ,40 + ,3949 + ,0 + ,12 + ,1 + ,37 + ,2727 + ,0 + ,12 + ,0 + ,42 + ,4262 + ,1 + ,12 + ,0 + ,40 + ,3523 + ,0 + ,12 + ,0 + ,36 + ,2955 + ,0 + ,12 + ,0 + ,40 + ,3239 + ,0 + ,12 + ,0 + ,40 + ,3665 + ,0 + ,12 + ,1 + ,39 + ,2960 + ,1 + ,12 + ,0 + ,40 + ,3353 + ,0 + ,12 + ,1 + ,40 + ,3864 + ,1 + ,12 + ,0 + ,39 + ,3125 + ,1 + ,12 + ,0 + ,40 + ,3182 + ,0 + ,12 + ,0 + ,40 + ,3438 + ,0 + ,12 + ,0 + ,40 + ,3921 + ,1 + ,12 + ,0 + ,39 + ,3728 + ,0 + ,12 + ,0 + ,38 + ,2639 + ,1 + ,12 + ,1 + ,40 + ,3324 + ,1 + ,12 + ,0 + ,40 + ,2869 + ,1 + ,12 + ,1 + ,38 + ,3352 + ,1 + ,12 + ,0 + ,38 + ,3040 + ,1 + ,12 + ,0 + ,40 + ,3381 + ,1 + ,12 + ,0 + ,40 + ,4517 + ,0 + ,12 + ,0 + ,40 + ,4375 + ,0 + ,12 + ,0 + ,40 + ,4233 + ,0 + ,12 + ,0 + ,40 + ,2983 + ,1 + ,12 + ,0 + ,33 + ,1720 + ,1 + ,12 + ,1 + ,39 + ,3555 + ,0 + ,12 + ,0 + ,40 + ,3494 + ,1 + ,12 + ,0 + ,39 + ,1932 + ,1 + ,12 + ,1 + ,39 + ,3068 + ,0 + ,12 + ,1 + ,40 + ,2756 + ,1 + ,12 + ,0 + ,40 + ,3864 + ,1 + ,12 + ,1 + ,43 + ,4690 + ,0 + ,12 + ,0 + ,41 + ,3977 + ,0 + ,12 + ,0 + ,39 + ,3864 + ,0 + ,12 + ,1 + ,42 + ,3636 + ,1 + ,12 + ,1 + ,35 + ,3154 + ,0 + ,12 + ,0 + ,39 + ,3352 + ,1 + ,12 + ,0 + ,38 + ,2443 + ,1 + ,12 + ,1 + ,34 + ,2220 + ,1 + ,12 + ,0 + ,39 + ,3240 + ,1 + ,12 + ,0 + ,41 + ,3087 + ,1 + ,12 + ,1 + ,37 + ,2955 + ,1 + ,12 + ,0 + ,39 + ,2931 + ,1 + ,12 + ,0 + ,38 + ,3350 + ,0 + ,12 + ,1 + ,41 + ,3324 + ,0 + ,12 + ,0 + ,42 + ,3438 + ,1 + ,12 + ,0 + ,42 + ,3438 + ,1 + ,12 + ,0 + ,39 + ,3430 + ,0 + ,12 + ,0 + ,39 + ,3750 + ,1 + ,12 + ,0 + ,40 + ,3520 + ,1 + ,12 + ,1 + ,40 + ,3480 + ,0 + ,12 + ,1 + ,40 + ,2955 + ,0 + ,12 + ,0 + ,43 + ,3722 + ,0 + ,12 + ,0 + ,41 + ,2898 + ,0 + ,12 + ,0 + ,38 + ,3665 + ,1 + ,12 + ,0 + ,41 + ,3324 + ,1 + ,12 + ,0 + ,40 + ,2841 + ,1 + ,12 + ,1 + ,41 + ,2690 + ,1 + ,12 + ,1 + ,40 + ,3780 + ,0 + ,12 + ,0 + ,28 + ,1165 + ,1 + ,12 + ,0 + ,40 + ,3381 + ,0 + ,12 + ,1 + ,40 + ,1900 + ,1 + ,12 + ,0 + ,38 + ,2205 + ,1 + ,12 + ,0 + ,43 + ,2790 + ,1 + ,12 + ,0 + ,38 + ,3380 + ,1 + ,12 + ,0 + ,38 + ,3060 + ,0 + ,12 + ,0 + ,40 + ,3068 + ,1 + ,12 + ,0 + ,39 + ,3280 + ,1 + ,12 + ,1 + ,39 + ,2443 + ,0 + ,12 + ,0 + ,38 + ,3750 + ,0 + ,12 + ,0 + ,39 + ,3125 + ,1 + ,12 + ,0 + ,37 + ,3799 + ,1 + ,12 + ,0 + ,39 + ,3665 + ,0 + ,12 + ,0 + ,39 + ,3409 + ,0 + ,12 + ,0 + ,40 + ,3011 + ,1 + ,12 + ,0 + ,41 + ,3930 + ,0 + ,12 + ,0 + ,42 + ,3267 + ,0 + ,12 + ,1 + ,39 + ,3693 + ,1 + ,12 + ,0 + ,40 + ,3438 + ,0 + ,12 + ,0 + ,37 + ,2557 + ,1 + ,12 + ,0 + ,37 + ,2637 + ,1 + ,12 + ,0 + ,42 + ,4034 + ,0 + ,12 + ,0 + ,38 + ,3381 + ,1 + ,12 + ,1 + ,36 + ,2960 + ,0 + ,12 + ,0 + ,38 + ,3210 + ,1 + ,12 + ,0 + ,39 + ,3572 + ,1 + ,12 + ,1 + ,39 + ,3007 + ,0 + ,12 + ,0 + ,42 + ,3977 + ,0 + ,12 + ,1 + ,41 + ,3693 + ,1 + ,12 + ,0 + ,40 + ,3011 + ,0 + ,12 + ,0 + ,39 + ,2955 + ,1 + ,12 + ,0 + ,41 + ,3296 + ,1 + ,12 + ,0 + ,40 + ,3605 + ,1 + ,12 + ,0 + ,38 + ,3722 + ,1 + ,12 + ,1 + ,38 + ,3267 + ,0 + ,12 + ,0 + ,40 + ,3608 + ,1 + ,12 + ,0 + ,39 + ,2983 + ,1 + ,12 + ,0 + ,34 + ,2386 + ,0 + ,12 + ,1 + ,26 + ,682 + ,0 + ,12 + ,0 + ,35 + ,2784 + ,0 + ,12 + ,1 + ,38 + ,2528 + ,1 + ,12 + ,0 + ,30 + ,1446 + ,1 + ,12 + ,0 + ,39 + ,3480 + ,1 + ,12 + ,0 + ,38 + ,2727 + ,1 + ,12 + ,0 + ,40 + ,3267 + ,1 + ,12 + ,0 + ,38 + ,2955 + ,1 + ,12 + ,0 + ,40 + ,3068 + ,0 + ,12 + ,0 + ,34 + ,2330 + ,0 + ,12 + ,1 + ,40 + ,3154 + ,1 + ,12 + ,0 + ,37 + ,3779 + ,0 + ,12 + ,0 + ,40 + ,3714 + ,1 + ,12 + ,0 + ,39 + ,3135 + ,1 + ,12 + ,0 + ,36 + ,2188 + ,1 + ,12 + ,0 + ,38 + ,4040 + ,0 + ,12 + ,1 + ,33 + ,1750 + ,0 + ,12 + ,0 + ,42 + ,4205 + ,1 + ,12 + ,0 + ,39 + ,3335 + ,1 + ,12 + ,0 + ,27 + ,1051 + ,0 + ,12 + ,0 + ,40 + ,2898 + ,0 + ,12 + ,0 + ,42 + ,3580 + ,1 + ,12 + ,0 + ,36 + ,2860 + ,1 + ,12 + ,0 + ,41 + ,3355 + ,0 + ,12 + ,1 + ,40 + ,2699 + ,1 + ,12 + ,1 + ,39 + ,2000 + ,0 + ,13 + ,0 + ,40 + ,3154 + ,0 + ,13 + ,0 + ,38 + ,3494 + ,1 + ,13 + ,0 + ,37 + ,3210 + ,1 + ,13 + ,1 + ,24 + ,521 + ,1 + ,13 + ,0 + ,39 + ,2898 + ,0 + ,13 + ,1 + ,39 + ,2386 + ,1 + ,13 + ,0 + ,38 + ,3154 + ,1 + ,13 + ,0 + ,39 + ,3220 + ,1 + ,13 + ,0 + ,40 + ,3510 + ,0 + ,13 + ,0 + ,41 + ,3608 + ,1 + ,13 + ,0 + ,38 + ,3225 + ,0 + ,13 + ,0 + ,40 + ,3494 + ,1 + ,13 + ,0 + ,38 + ,2671 + ,0 + ,13 + ,1 + ,40 + ,2841 + ,1 + ,13 + ,0 + ,41 + ,3466 + ,0 + ,13 + ,0 + ,40 + ,3949 + ,1 + ,13 + ,1 + ,39 + ,3095 + ,1 + ,13 + ,0 + ,40 + ,2784 + ,1 + ,13 + ,0 + ,37 + ,3523 + ,1 + ,13 + ,0 + ,40 + ,4240 + ,1 + ,13 + ,0 + ,41 + ,3550 + ,1 + ,13 + ,0 + ,39 + ,3296 + ,1 + ,13 + ,1 + ,40 + ,3170 + ,1 + ,13 + ,0 + ,37 + ,3693 + ,1 + ,13 + ,0 + ,39 + ,3551 + ,1 + ,13 + ,0 + ,39 + ,3100 + ,1 + ,13 + ,1 + ,27 + ,1095 + ,1 + ,13 + ,0 + ,40 + ,3523 + ,0 + ,13 + ,0 + ,41 + ,3750 + ,1 + ,13 + ,1 + ,40 + ,3182 + ,1 + ,13 + ,1 + ,38 + ,2617 + ,1 + ,13 + ,0 + ,34 + ,1421 + ,1 + ,13 + ,1 + ,38 + ,2869 + ,0 + ,13 + ,0 + ,41 + ,3939 + ,1 + ,13 + ,0 + ,40 + ,3494 + ,1 + ,13 + ,0 + ,40 + ,3608 + ,1 + ,13 + ,0 + ,38 + ,3210 + ,1 + ,13 + ,0 + ,33 + ,1920 + ,0 + ,13 + ,0 + ,38 + ,2898 + ,1 + ,13 + ,0 + ,40 + ,2813 + ,1 + ,13 + ,0 + ,39 + ,2825 + ,1 + ,13 + ,0 + ,42 + ,4460 + ,1 + ,13 + ,0 + ,39 + ,2975 + ,1 + ,13 + ,1 + ,39 + ,2810 + ,0 + ,13 + ,0 + ,40 + ,3864 + ,1 + ,13 + ,1 + ,39 + ,2869 + ,1 + ,13 + ,0 + ,40 + ,3977 + ,0 + ,13 + ,0 + ,41 + ,3665 + ,1 + ,13 + ,0 + ,40 + ,3494 + ,1 + ,13 + ,0 + ,41 + ,3180 + ,1 + ,13 + ,0 + ,40 + ,3040 + ,0 + ,13 + ,0 + ,39 + ,4091 + ,1 + ,13 + ,0 + ,38 + ,3182 + ,1 + ,13 + ,0 + ,40 + ,3296 + ,1 + ,13 + ,0 + ,40 + ,4091 + ,0 + ,13 + ,0 + ,41 + ,3779 + ,0 + ,13 + ,1 + ,20 + ,410 + ,1 + ,13 + ,0 + ,40 + ,3750 + ,1 + ,13 + ,0 + ,40 + ,3010 + ,1 + ,13 + ,0 + ,40 + ,3324 + ,1 + ,13 + ,0 + ,37 + ,2926 + ,1 + ,13 + ,0 + ,38 + ,3580 + ,1 + ,13 + ,0 + ,40 + ,3636 + ,1 + ,13 + ,0 + ,40 + ,3324 + ,1 + ,13 + ,0 + ,39 + ,3324 + ,1 + ,14 + ,0 + ,40 + ,3352 + ,0 + ,14 + ,0 + ,40 + ,3352 + ,0 + ,14 + ,0 + ,37 + ,3125 + ,0 + ,14 + ,0 + ,41 + ,3448 + ,0 + ,14 + ,0 + ,41 + ,3494 + ,1 + ,14 + ,0 + ,40 + ,3608 + ,1 + ,14 + ,0 + ,40 + ,3779 + ,1 + ,14 + ,0 + ,40 + ,3860 + ,0 + ,14 + ,0 + ,40 + ,3921 + ,0 + ,14 + ,0 + ,36 + ,2642 + ,1 + ,14 + ,0 + ,40 + ,3285 + ,1 + ,14 + ,0 + ,40 + ,3835 + ,1 + ,14 + ,0 + ,40 + ,3250 + ,0 + ,14 + ,0 + ,40 + ,3210 + ,1 + ,14 + ,0 + ,40 + ,3892 + ,0 + ,14 + ,0 + ,42 + ,4546 + ,0 + ,14 + ,0 + ,38 + ,3864 + ,1 + ,14 + ,0 + ,38 + ,3210 + ,1 + ,14 + ,0 + ,41 + ,3420 + ,1 + ,14 + ,0 + ,40 + ,3154 + ,1 + ,14 + ,0 + ,41 + ,2472 + ,0 + ,14 + ,0 + ,41 + ,3665 + ,0 + ,14 + ,0 + ,38 + ,3580 + ,0 + ,14 + ,0 + ,40 + ,3381 + ,0 + ,14 + ,0 + ,41 + ,3125 + ,0 + ,14 + ,0 + ,36 + ,2880 + ,1 + ,14 + ,0 + ,38 + ,2600 + ,0 + ,14 + ,0 + ,40 + ,3835 + ,1 + ,14 + ,1 + ,40 + ,3612 + ,1 + ,14 + ,0 + ,21 + ,320 + ,1 + ,14 + ,1 + ,40 + ,3244 + ,1 + ,14 + ,1 + ,40 + ,3296 + ,0 + ,14 + ,0 + ,40 + ,3693 + ,0 + ,14 + ,0 + ,40 + ,3679 + ,1 + ,14 + ,0 + ,40 + ,3450 + ,1 + ,14 + ,0 + ,39 + ,3977 + ,1 + ,14 + ,0 + ,39 + ,3154 + ,1 + ,14 + ,0 + ,39 + ,4034 + ,1 + ,14 + ,0 + ,42 + ,3160 + ,1 + ,14 + ,0 + ,22 + ,682 + ,1 + ,14 + ,0 + ,38 + ,3045 + ,0 + ,14 + ,0 + ,40 + ,3409 + ,1 + ,14 + ,0 + ,40 + ,3352 + ,1 + ,14 + ,1 + ,40 + ,3665 + ,1 + ,14 + ,1 + ,40 + ,3020 + ,1 + ,14 + ,0 + ,40 + ,2557 + ,0 + ,14 + ,0 + ,40 + ,4404 + ,1 + ,14 + ,0 + ,39 + ,2840 + ,1 + ,14 + ,1 + ,40 + ,3500 + ,0 + ,14 + ,0 + ,40 + ,4063 + ,0 + ,14 + ,0 + ,40 + ,4574 + ,1 + ,14 + ,0 + ,40 + ,4006 + ,0 + ,14 + ,1 + ,40 + ,3636 + ,1 + ,14 + ,0 + ,42 + ,3494 + ,1 + ,14 + ,0 + ,41 + ,3210 + ,1 + ,14 + ,0 + ,40 + ,3409 + ,1 + ,14 + ,0 + ,37 + ,2585 + ,1 + ,14 + ,1 + ,38 + ,2555 + ,0 + ,14 + ,0 + ,41 + ,3523 + ,1 + ,14 + ,0 + ,37 + ,3155 + ,0 + ,14 + ,1 + ,40 + ,3381 + ,1 + ,14 + ,0 + ,40 + ,2983 + ,0 + ,14 + ,1 + ,42 + ,3892 + ,1 + ,14 + ,0 + ,40 + ,3324 + ,0 + ,14 + ,0 + ,40 + ,3296 + ,0 + ,14 + ,0 + ,37 + ,2443 + ,0 + ,14 + ,1 + ,39 + ,3580 + ,0 + ,14 + ,0 + ,40 + ,2784 + ,0 + ,14 + ,1 + ,39 + ,3835 + ,0 + ,14 + ,0 + ,40 + ,3721 + ,0 + ,14 + ,0 + ,39 + ,3210 + ,0 + ,14 + ,0 + ,40 + ,3892 + ,1 + ,14 + ,0 + ,39 + ,2841 + ,1 + ,14 + ,0 + ,38 + ,3068 + ,1 + ,14 + ,0 + ,36 + ,2690 + ,1 + ,14 + ,0 + ,40 + ,3025 + ,0 + ,14 + ,0 + ,39 + ,4063 + ,0 + ,14 + ,0 + ,39 + ,3551 + ,1 + ,14 + ,0 + ,38 + ,3625 + ,1 + ,14 + ,1 + ,38 + ,3210 + ,1 + ,14 + ,0 + ,38 + ,2955 + ,1 + ,14 + ,0 + ,39 + ,3230 + ,1 + ,14 + ,0 + ,39 + ,3154 + ,0 + ,15 + ,0 + ,40 + ,3438 + ,1 + ,15 + ,0 + ,40 + ,4404 + ,1 + ,15 + ,0 + ,39 + ,3615 + ,0 + ,15 + ,0 + ,40 + ,4659 + ,1 + ,15 + ,1 + ,40 + ,2841 + ,0 + ,15 + ,0 + ,40 + ,4006 + ,1 + ,15 + ,0 + ,39 + ,2869 + ,1 + ,15 + ,0 + ,22 + ,470 + ,0 + ,15 + ,0 + ,40 + ,3296 + ,1 + ,15 + ,0 + ,40 + ,3728 + ,0 + ,15 + ,0 + ,39 + ,3068 + ,0 + ,15 + ,0 + ,38 + ,2585 + ,1 + ,15 + ,0 + ,39 + ,3296 + ,0 + ,15 + ,0 + ,40 + ,3438 + ,0 + ,15 + ,0 + ,41 + ,3068 + ,1 + ,15 + ,0 + ,39 + ,3438 + ,1 + ,15 + ,0 + ,35 + ,3040 + ,0 + ,15 + ,0 + ,39 + ,3523 + ,1 + ,15 + ,0 + ,41 + ,2926 + ,1 + ,15 + ,0 + ,37 + ,2301 + ,0 + ,15 + ,0 + ,39 + ,3409 + ,1 + ,15 + ,0 + ,40 + ,3324 + ,0 + ,15 + ,0 + ,38 + ,3580 + ,1 + ,15 + ,0 + ,42 + ,3722 + ,1 + ,15 + ,0 + ,38 + ,2869 + ,0 + ,15 + ,0 + ,42 + ,4119 + ,1 + ,15 + ,1 + ,40 + ,3540 + ,0 + ,16 + ,0 + ,39 + ,3352 + ,1 + ,16 + ,1 + ,38 + ,4347 + ,0 + ,16 + ,0 + ,38 + ,3097 + ,1 + ,16 + ,0 + ,39 + ,4091 + ,1 + ,16 + ,0 + ,40 + ,3381 + ,0 + ,16 + ,0 + ,40 + ,2727 + ,0 + ,16 + ,0 + ,40 + ,3864 + ,0 + ,16 + ,0 + ,42 + ,2869 + ,1 + ,16 + ,0 + ,37 + ,3523 + ,0 + ,16 + ,0 + ,39 + ,2528 + ,1 + ,16 + ,0 + ,38 + ,3324 + ,1 + ,16 + ,1 + ,39 + ,2983 + ,0 + ,16 + ,0 + ,42 + ,4148 + ,1 + ,16 + ,0 + ,35 + ,2760 + ,1 + ,16 + ,0 + ,41 + ,4480 + ,0 + ,16 + ,0 + ,39 + ,3015 + ,1 + ,16 + ,0 + ,40 + ,2898 + ,1 + ,16 + ,0 + ,41 + ,3720 + ,1 + ,16 + ,0 + ,38 + ,4176 + ,0 + ,16 + ,0 + ,42 + ,3608 + ,0 + ,16 + ,0 + ,36 + ,2756 + ,1 + ,16 + ,1 + ,38 + ,2950 + ,0 + ,16 + ,0 + ,42 + ,3210 + ,0 + ,16 + ,0 + ,40 + ,2926 + ,1 + ,16 + ,0 + ,40 + ,3154 + ,0 + ,16 + ,0 + ,39 + ,3090 + ,1 + ,16 + ,0 + ,39 + ,3324 + ,0 + ,16 + ,0 + ,39 + ,2955 + ,1 + ,16 + ,0 + ,41 + ,4432 + ,0 + ,16 + ,0 + ,39 + ,3864 + ,0 + ,16 + ,0 + ,37 + ,3068 + ,1 + ,16 + ,0 + ,39 + ,3665 + ,0 + ,16 + ,0 + ,41 + ,3352 + ,1 + ,16 + ,0 + ,39 + ,2815 + ,0 + ,16 + ,0 + ,40 + ,3630 + ,0 + ,16 + ,0 + ,39 + ,3040 + ,1 + ,16 + ,0 + ,38 + ,3011 + ,0 + ,16 + ,0 + ,38 + ,3523 + ,0 + ,16 + ,0 + ,39 + ,3494 + ,0 + ,16 + ,1 + ,38 + ,3466 + ,0 + ,16 + ,0 + ,41 + ,3921 + ,0 + ,16 + ,0 + ,42 + ,3835 + ,1 + ,16 + ,0 + ,41 + ,3381 + ,0 + ,16 + ,0 + ,40 + ,3494 + ,1 + ,16 + ,0 + ,38 + ,4100 + ,0 + ,16 + ,0 + ,39 + ,3779 + ,0 + ,16 + ,0 + ,40 + ,3864 + ,0 + ,16 + ,0 + ,39 + ,3324 + ,1 + ,16 + ,0 + ,37 + ,2479 + ,0 + ,16 + ,0 + ,39 + ,3693 + ,0 + ,16 + ,0 + ,40 + ,4148 + ,0 + ,16 + ,0 + ,40 + ,3210 + ,1 + ,16 + ,0 + ,41 + ,3458 + ,0 + ,16 + ,0 + ,40 + ,4347 + ,1 + ,16 + ,1 + ,38 + ,3296 + ,1 + ,16 + ,0 + ,37 + ,2720 + ,0 + ,16 + ,0 + ,41 + ,3636 + ,0 + ,16 + ,0 + ,35 + ,3665 + ,0 + ,16 + ,0 + ,39 + ,3551 + ,1 + ,16 + ,0 + ,39 + ,4830 + ,0 + ,16 + ,0 + ,41 + ,3864 + ,1 + ,16 + ,0 + ,40 + ,3210 + ,1 + ,16 + ,0 + ,40 + ,3807 + ,0 + ,16 + ,0 + ,38 + ,4460 + ,0 + ,16 + ,0 + ,40 + ,3608 + ,1 + ,16 + ,0 + ,38 + ,3502 + ,0 + ,16 + ,0 + ,39 + ,3239 + ,0 + ,16 + ,0 + ,41 + ,3892 + ,1 + ,16 + ,0 + ,37 + ,2995 + ,0 + ,16 + ,0 + ,39 + ,3665 + ,0 + ,16 + ,1 + ,40 + ,3381 + ,0 + ,16 + ,1 + ,41 + ,3466 + ,0 + ,16 + ,0 + ,40 + ,3892 + ,0 + ,16 + ,0 + ,41 + ,3438 + ,1 + ,16 + ,0 + ,40 + ,2955 + ,0 + ,16 + ,0 + ,41 + ,3190 + ,0 + ,16 + ,0 + ,39 + ,3060 + ,1 + ,16 + ,0 + ,39 + ,2699 + ,0 + ,16 + ,0 + ,40 + ,4375 + ,0 + ,16 + ,0 + ,40 + ,3722 + ,0 + ,16 + ,0 + ,38 + ,3580 + ,1 + ,16 + ,0 + ,39 + ,3494 + ,0 + ,16 + ,0 + ,41 + ,3835 + ,1 + ,16 + ,0 + ,40 + ,3750 + ,0 + ,16 + ,0 + ,39 + ,3152 + ,1 + ,17 + ,0 + ,38 + ,3523 + ,0 + ,17 + ,0 + ,40 + ,3636 + ,0 + ,17 + ,0 + ,38 + ,3523 + ,1 + ,17 + ,0 + ,39 + ,2727 + ,0 + ,17 + ,0 + ,38 + ,3040 + ,0 + ,17 + ,0 + ,38 + ,3182 + ,1 + ,17 + ,0 + ,40 + ,3494 + ,1 + ,17 + ,0 + ,39 + ,2813 + ,0 + ,17 + ,0 + ,40 + ,2727 + ,1 + ,17 + ,0 + ,39 + ,3875 + ,0 + ,17 + ,0 + ,40 + ,3921 + ,0 + ,17 + ,0 + ,39 + ,2642 + ,0 + ,17 + ,0 + ,39 + ,2410 + ,1 + ,17 + ,0 + ,31 + ,997 + ,0 + ,17 + ,0 + ,38 + ,3409 + ,1 + ,17 + ,0 + ,39 + ,2557 + ,0 + ,17 + ,0 + ,41 + ,3523 + ,0 + ,17 + ,0 + ,41 + ,3921 + ,0 + ,17 + ,0 + ,42 + ,3486 + ,0 + ,17 + ,0 + ,39 + ,3739 + ,0 + ,17 + ,0 + ,38 + ,3864 + ,0 + ,17 + ,0 + ,39 + ,3267 + ,0 + ,17 + ,0 + ,39 + ,3523 + ,0 + ,17 + ,0 + ,37 + ,3011 + ,1 + ,17 + ,0 + ,33 + ,1875 + ,0 + ,17 + ,0 + ,36 + ,2330 + ,0 + ,17 + ,0 + ,41 + ,3351 + ,0 + ,17 + ,0 + ,40 + ,4006 + ,0 + ,17 + ,0 + ,38 + ,2671 + ,0 + ,17 + ,0 + ,40 + ,2784 + ,1 + ,17 + ,0 + ,38 + ,2700 + ,0 + ,17 + ,0 + ,38 + ,3324 + ,0 + ,17 + ,0 + ,41 + ,3660 + ,0 + ,17 + ,0 + ,40 + ,3125 + ,1 + ,17 + ,0 + ,40 + ,3466 + ,0 + ,17 + ,0 + ,37 + ,3197 + ,0 + ,17 + ,0 + ,39 + ,3296 + ,1 + ,17 + ,0 + ,38 + ,3415 + ,1 + ,17 + ,0 + ,40 + ,3320 + ,1 + ,17 + ,0 + ,40 + ,3315 + ,0 + ,17 + ,0 + ,40 + ,3068 + ,0 + ,17 + ,0 + ,40 + ,4318 + ,0 + ,17 + ,0 + ,40 + ,4034 + ,0 + ,17 + ,0 + ,40 + ,3608 + ,0 + ,17 + ,0 + ,40 + ,3352) + ,dim=c(5 + ,1115) + ,dimnames=list(c('black' + ,'edu' + ,'smoke' + ,'gestate' + ,'grams') + ,1:1115)) > y <- array(NA,dim=c(5,1115),dimnames=list(c('black','edu','smoke','gestate','grams'),1:1115)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par4 = 'no' > par3 = '3' > par2 = 'none' > par1 = '5' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Dr. Ian E. Holliday > #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: > #Technical description: > library(party) Loading required package: survival Loading required package: splines Loading required package: grid Loading required package: modeltools Loading required package: stats4 Loading required package: coin Loading required package: mvtnorm Loading required package: zoo Loading required package: sandwich Loading required package: strucchange Loading required package: vcd Loading required package: MASS Loading required package: colorspace > library(Hmisc) Attaching package: 'Hmisc' The following object(s) are masked from 'package:survival': untangle.specials The following object(s) are masked from 'package:base': format.pval, round.POSIXt, trunc.POSIXt, units > par1 <- as.numeric(par1) > par3 <- as.numeric(par3) > x <- data.frame(t(y)) > is.data.frame(x) [1] TRUE > x <- x[!is.na(x[,par1]),] > k <- length(x[1,]) > n <- length(x[,1]) > colnames(x)[par1] [1] "grams" > x[,par1] [1] 2898 994 3977 3040 3523 3100 3670 3097 3040 3239 2955 2200 3182 3510 [15] 3381 3530 2985 3374 3765 2715 3640 3040 3670 2775 2680 2699 2557 3835 [29] 3381 2766 3125 2557 2950 3779 2760 3011 2415 3125 3154 3170 2784 2475 [43] 2642 4034 3988 3370 2955 2443 4000 3630 2920 2727 3800 3771 3075 2415 [57] 2890 3892 3125 3892 3296 3120 1710 3068 3779 2898 3381 3324 2727 2642 [71] 2244 2528 3352 3551 3267 3210 3395 3068 3420 2985 2147 2863 3125 4290 [85] 2585 2869 3097 2756 3125 3921 3438 2980 2869 4830 3409 3807 3220 2756 [99] 3097 3615 2926 4050 2840 2671 3540 2900 3779 2740 2320 2905 3070 3352 [113] 2581 4091 3235 3722 3438 3475 3438 3409 3040 3800 3182 3494 2131 1790 [127] 3864 2955 3949 2273 2620 2205 2860 3381 3190 3267 3665 3160 3216 2415 [141] 3200 3352 3470 3580 3750 3035 3438 3068 3438 3750 3409 1790 3170 3800 [155] 2825 2756 3040 2585 2756 3523 909 3779 2415 3665 3295 2699 3267 3460 [169] 1250 3730 3523 3011 3740 3949 2330 3182 3040 2557 3147 2671 2557 3200 [183] 3203 2756 1619 2330 3381 2983 2600 3585 3693 3020 1800 2290 3750 3551 [197] 3779 3125 3470 3210 3100 3182 3355 2750 4006 3438 3068 2650 3005 2993 [211] 3494 2795 3260 2959 4020 2949 3900 2926 3239 3409 3630 3466 2330 3949 [225] 2813 3137 3296 3267 2940 3323 3409 3296 3807 3296 2983 1680 3551 2385 [239] 3356 3466 2699 4233 2890 2443 3580 3551 2030 3125 2700 2216 3182 3040 [253] 3011 2800 3425 2760 3665 1733 2850 2485 3381 3065 3167 2635 3665 3300 [267] 3636 2813 3920 2671 3182 3409 3600 2557 3890 4148 2030 2890 2373 2977 [281] 2840 2658 2983 3700 3608 3559 3540 3440 3324 3494 3636 2813 3239 3068 [295] 3210 3892 3864 2784 3551 3608 3130 3466 3608 2225 2585 2917 2415 4489 [309] 3977 3080 3722 3097 2960 3807 3210 3210 284 3805 3267 3722 3864 4318 [323] 3807 455 3722 3324 3665 3240 3182 3040 3665 4148 3636 1591 4135 3352 [337] 3011 2301 2740 2450 2472 4148 3352 4034 3807 2926 2699 2990 2898 3864 [351] 2970 3655 2890 3980 2756 4148 3494 2830 2983 3430 2159 3125 3239 4262 [365] 3459 1190 3182 3352 2898 3154 1165 2102 3110 2367 3892 3011 3125 3490 [379] 2955 3977 3665 3355 3182 3239 2955 3182 2784 3802 3270 2225 3182 3210 [393] 3779 2415 4675 3352 2900 3864 2841 2557 3494 3438 2445 3714 2935 3381 [407] 3864 3240 2980 2710 3570 3864 2500 1160 2443 3040 3864 3205 1875 3835 [421] 2841 2926 2983 3267 2890 2443 3352 3438 3120 4176 4148 3270 3523 3020 [435] 3620 3300 3779 3665 2869 3154 2810 2131 4233 3860 3920 3693 3580 3608 [449] 2330 2813 3068 3250 3175 2940 3267 3645 4205 1818 2983 2756 2642 3182 [463] 3125 3949 3466 3250 3750 3466 2557 3722 3551 3097 1763 3610 3551 2460 [477] 2358 3779 2671 3210 4375 4160 2841 3352 2790 2540 2890 3182 2557 3551 [491] 3324 2620 3068 3150 3636 3130 2784 3520 4205 3210 3523 3125 3097 3125 [505] 2930 3665 3693 3835 3636 4165 3035 2000 2740 3608 3750 4091 3267 3097 [519] 3381 3040 3779 3523 3551 3636 3154 2841 2960 3097 3600 3409 2880 2690 [533] 3090 3409 2159 3921 2940 3400 3239 3636 4347 2977 3636 2671 3949 3608 [547] 3710 2216 2500 1450 3296 3510 3950 4176 3636 3125 3845 3523 2784 3494 [561] 4233 3779 3720 3302 3693 2580 2597 2926 4505 3977 3182 2160 3402 3239 [575] 3608 3267 3995 3068 4065 3346 3154 2515 3182 2443 3807 3130 3722 2727 [589] 3011 793 3494 2386 3466 2760 4659 3608 2557 2671 3740 2614 3949 3630 [603] 2898 3940 701 3977 2472 3011 3352 3494 4091 3530 3010 3068 3636 3665 [617] 3430 3182 3300 2557 2550 3097 3154 3625 3330 4347 3296 3154 4050 3239 [631] 3608 2500 3438 3977 2780 2189 1760 2360 3608 3892 2870 2926 3905 3020 [645] 3779 3381 3280 3154 3494 3523 1861 3360 3977 3296 3523 3693 3921 653 [659] 3494 3500 2863 3494 3335 3750 3636 4190 3540 3381 3580 2813 3494 4065 [673] 3520 3551 4063 3210 3640 2770 3750 3695 3305 3892 3210 2851 3175 3949 [687] 2727 4262 3523 2955 3239 3665 2960 3353 3864 3125 3182 3438 3921 3728 [701] 2639 3324 2869 3352 3040 3381 4517 4375 4233 2983 1720 3555 3494 1932 [715] 3068 2756 3864 4690 3977 3864 3636 3154 3352 2443 2220 3240 3087 2955 [729] 2931 3350 3324 3438 3438 3430 3750 3520 3480 2955 3722 2898 3665 3324 [743] 2841 2690 3780 1165 3381 1900 2205 2790 3380 3060 3068 3280 2443 3750 [757] 3125 3799 3665 3409 3011 3930 3267 3693 3438 2557 2637 4034 3381 2960 [771] 3210 3572 3007 3977 3693 3011 2955 3296 3605 3722 3267 3608 2983 2386 [785] 682 2784 2528 1446 3480 2727 3267 2955 3068 2330 3154 3779 3714 3135 [799] 2188 4040 1750 4205 3335 1051 2898 3580 2860 3355 2699 2000 3154 3494 [813] 3210 521 2898 2386 3154 3220 3510 3608 3225 3494 2671 2841 3466 3949 [827] 3095 2784 3523 4240 3550 3296 3170 3693 3551 3100 1095 3523 3750 3182 [841] 2617 1421 2869 3939 3494 3608 3210 1920 2898 2813 2825 4460 2975 2810 [855] 3864 2869 3977 3665 3494 3180 3040 4091 3182 3296 4091 3779 410 3750 [869] 3010 3324 2926 3580 3636 3324 3324 3352 3352 3125 3448 3494 3608 3779 [883] 3860 3921 2642 3285 3835 3250 3210 3892 4546 3864 3210 3420 3154 2472 [897] 3665 3580 3381 3125 2880 2600 3835 3612 320 3244 3296 3693 3679 3450 [911] 3977 3154 4034 3160 682 3045 3409 3352 3665 3020 2557 4404 2840 3500 [925] 4063 4574 4006 3636 3494 3210 3409 2585 2555 3523 3155 3381 2983 3892 [939] 3324 3296 2443 3580 2784 3835 3721 3210 3892 2841 3068 2690 3025 4063 [953] 3551 3625 3210 2955 3230 3154 3438 4404 3615 4659 2841 4006 2869 470 [967] 3296 3728 3068 2585 3296 3438 3068 3438 3040 3523 2926 2301 3409 3324 [981] 3580 3722 2869 4119 3540 3352 4347 3097 4091 3381 2727 3864 2869 3523 [995] 2528 3324 2983 4148 2760 4480 3015 2898 3720 4176 3608 2756 2950 3210 [1009] 2926 3154 3090 3324 2955 4432 3864 3068 3665 3352 2815 3630 3040 3011 [1023] 3523 3494 3466 3921 3835 3381 3494 4100 3779 3864 3324 2479 3693 4148 [1037] 3210 3458 4347 3296 2720 3636 3665 3551 4830 3864 3210 3807 4460 3608 [1051] 3502 3239 3892 2995 3665 3381 3466 3892 3438 2955 3190 3060 2699 4375 [1065] 3722 3580 3494 3835 3750 3152 3523 3636 3523 2727 3040 3182 3494 2813 [1079] 2727 3875 3921 2642 2410 997 3409 2557 3523 3921 3486 3739 3864 3267 [1093] 3523 3011 1875 2330 3351 4006 2671 2784 2700 3324 3660 3125 3466 3197 [1107] 3296 3415 3320 3315 3068 4318 4034 3608 3352 > if (par2 == 'kmeans') { + cl <- kmeans(x[,par1], par3) + print(cl) + clm <- matrix(cbind(cl$centers,1:par3),ncol=2) + clm <- clm[sort.list(clm[,1]),] + for (i in 1:par3) { + cl$cluster[cl$cluster==clm[i,2]] <- paste('C',i,sep='') + } + cl$cluster <- as.factor(cl$cluster) + print(cl$cluster) + x[,par1] <- cl$cluster + } > if (par2 == 'quantiles') { + x[,par1] <- cut2(x[,par1],g=par3) + } > if (par2 == 'hclust') { + hc <- hclust(dist(x[,par1])^2, 'cen') + print(hc) + memb <- cutree(hc, k = par3) + dum <- c(mean(x[memb==1,par1])) + for (i in 2:par3) { + dum <- c(dum, mean(x[memb==i,par1])) + } + hcm <- matrix(cbind(dum,1:par3),ncol=2) + hcm <- hcm[sort.list(hcm[,1]),] + for (i in 1:par3) { + memb[memb==hcm[i,2]] <- paste('C',i,sep='') + } + memb <- as.factor(memb) + print(memb) + x[,par1] <- memb + } > if (par2=='equal') { + ed <- cut(as.numeric(x[,par1]),par3,labels=paste('C',1:par3,sep='')) + x[,par1] <- as.factor(ed) + } > table(x[,par1]) 284 320 410 455 470 521 653 682 701 793 909 994 997 1051 1095 1160 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1165 1190 1250 1421 1446 1450 1591 1619 1680 1710 1720 1733 1750 1760 1763 1790 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1800 1818 1861 1875 1900 1920 1932 2000 2030 2102 2131 2147 2159 2160 2188 2189 1 1 1 2 1 1 1 2 2 1 2 1 2 1 1 1 2200 2205 2216 2220 2225 2244 2273 2290 2301 2320 2330 2358 2360 2367 2373 2385 1 2 2 1 2 1 1 1 2 1 6 1 1 1 1 1 2386 2410 2415 2443 2445 2450 2460 2472 2475 2479 2485 2500 2515 2528 2540 2550 3 1 6 8 1 1 1 3 1 1 1 3 1 3 1 1 2555 2557 2580 2581 2585 2597 2600 2614 2617 2620 2635 2637 2639 2642 2650 2658 1 13 1 1 5 1 2 1 1 2 1 1 1 5 1 1 2671 2680 2690 2699 2700 2710 2715 2720 2727 2740 2750 2756 2760 2766 2770 2775 8 1 3 6 2 1 1 1 8 3 1 9 4 1 1 1 2780 2784 2790 2795 2800 2810 2813 2815 2825 2830 2840 2841 2850 2851 2860 2863 1 9 2 1 1 2 7 1 2 1 3 8 1 1 2 2 2869 2870 2880 2890 2898 2900 2905 2917 2920 2926 2930 2931 2935 2940 2949 2950 9 1 2 6 10 2 1 1 1 9 1 1 1 3 1 2 2955 2959 2960 2970 2975 2977 2980 2983 2985 2990 2993 2995 3005 3007 3010 3011 13 1 4 1 1 2 2 10 2 1 1 1 1 1 2 11 3015 3020 3025 3035 3040 3045 3060 3065 3068 3070 3075 3080 3087 3090 3095 3097 1 4 1 2 15 1 2 1 17 1 1 1 1 2 1 10 3100 3110 3120 3125 3130 3135 3137 3147 3150 3152 3154 3155 3160 3167 3170 3175 3 1 2 18 3 1 1 1 1 1 16 1 2 1 3 2 3180 3182 3190 3197 3200 3203 3205 3210 3216 3220 3225 3230 3235 3239 3240 3244 1 20 2 1 2 1 1 21 1 2 1 1 1 10 3 1 3250 3260 3267 3270 3280 3285 3295 3296 3300 3302 3305 3315 3320 3323 3324 3330 3 1 13 2 2 1 1 16 3 1 1 1 1 1 16 1 3335 3346 3350 3351 3352 3353 3355 3356 3360 3370 3374 3380 3381 3395 3400 3402 2 1 1 1 18 1 3 1 1 1 1 1 18 1 1 1 3409 3415 3420 3425 3430 3438 3440 3448 3450 3458 3459 3460 3466 3470 3475 3480 13 1 2 1 3 17 1 1 1 1 1 1 10 2 1 2 3486 3490 3494 3500 3502 3510 3520 3523 3530 3540 3550 3551 3555 3559 3570 3572 1 1 23 2 1 3 3 20 2 4 1 13 1 1 1 1 3580 3585 3600 3605 3608 3610 3612 3615 3620 3625 3630 3636 3640 3645 3655 3660 10 1 2 1 17 1 1 2 1 2 4 16 2 1 1 1 3665 3670 3679 3693 3695 3700 3710 3714 3720 3721 3722 3728 3730 3739 3740 3750 19 2 1 10 1 1 1 2 2 1 10 2 1 1 2 12 3765 3771 3779 3780 3799 3800 3802 3805 3807 3835 3845 3860 3864 3875 3890 3892 1 1 15 1 1 3 1 1 7 8 1 2 18 1 1 11 3900 3905 3920 3921 3930 3939 3940 3949 3950 3977 3980 3988 3995 4000 4006 4020 1 1 2 8 1 1 1 8 1 11 1 1 1 1 4 1 4034 4040 4050 4063 4065 4091 4100 4119 4135 4148 4160 4165 4176 4190 4205 4233 5 1 2 3 2 6 1 1 1 7 1 1 3 1 3 4 4240 4262 4290 4318 4347 4375 4404 4432 4460 4480 4489 4505 4517 4546 4574 4659 1 2 1 2 4 3 2 1 2 1 1 1 1 1 1 2 4675 4690 4830 1 1 2 > colnames(x) [1] "black" "edu" "smoke" "gestate" "grams" > colnames(x)[par1] [1] "grams" > x[,par1] [1] 2898 994 3977 3040 3523 3100 3670 3097 3040 3239 2955 2200 3182 3510 [15] 3381 3530 2985 3374 3765 2715 3640 3040 3670 2775 2680 2699 2557 3835 [29] 3381 2766 3125 2557 2950 3779 2760 3011 2415 3125 3154 3170 2784 2475 [43] 2642 4034 3988 3370 2955 2443 4000 3630 2920 2727 3800 3771 3075 2415 [57] 2890 3892 3125 3892 3296 3120 1710 3068 3779 2898 3381 3324 2727 2642 [71] 2244 2528 3352 3551 3267 3210 3395 3068 3420 2985 2147 2863 3125 4290 [85] 2585 2869 3097 2756 3125 3921 3438 2980 2869 4830 3409 3807 3220 2756 [99] 3097 3615 2926 4050 2840 2671 3540 2900 3779 2740 2320 2905 3070 3352 [113] 2581 4091 3235 3722 3438 3475 3438 3409 3040 3800 3182 3494 2131 1790 [127] 3864 2955 3949 2273 2620 2205 2860 3381 3190 3267 3665 3160 3216 2415 [141] 3200 3352 3470 3580 3750 3035 3438 3068 3438 3750 3409 1790 3170 3800 [155] 2825 2756 3040 2585 2756 3523 909 3779 2415 3665 3295 2699 3267 3460 [169] 1250 3730 3523 3011 3740 3949 2330 3182 3040 2557 3147 2671 2557 3200 [183] 3203 2756 1619 2330 3381 2983 2600 3585 3693 3020 1800 2290 3750 3551 [197] 3779 3125 3470 3210 3100 3182 3355 2750 4006 3438 3068 2650 3005 2993 [211] 3494 2795 3260 2959 4020 2949 3900 2926 3239 3409 3630 3466 2330 3949 [225] 2813 3137 3296 3267 2940 3323 3409 3296 3807 3296 2983 1680 3551 2385 [239] 3356 3466 2699 4233 2890 2443 3580 3551 2030 3125 2700 2216 3182 3040 [253] 3011 2800 3425 2760 3665 1733 2850 2485 3381 3065 3167 2635 3665 3300 [267] 3636 2813 3920 2671 3182 3409 3600 2557 3890 4148 2030 2890 2373 2977 [281] 2840 2658 2983 3700 3608 3559 3540 3440 3324 3494 3636 2813 3239 3068 [295] 3210 3892 3864 2784 3551 3608 3130 3466 3608 2225 2585 2917 2415 4489 [309] 3977 3080 3722 3097 2960 3807 3210 3210 284 3805 3267 3722 3864 4318 [323] 3807 455 3722 3324 3665 3240 3182 3040 3665 4148 3636 1591 4135 3352 [337] 3011 2301 2740 2450 2472 4148 3352 4034 3807 2926 2699 2990 2898 3864 [351] 2970 3655 2890 3980 2756 4148 3494 2830 2983 3430 2159 3125 3239 4262 [365] 3459 1190 3182 3352 2898 3154 1165 2102 3110 2367 3892 3011 3125 3490 [379] 2955 3977 3665 3355 3182 3239 2955 3182 2784 3802 3270 2225 3182 3210 [393] 3779 2415 4675 3352 2900 3864 2841 2557 3494 3438 2445 3714 2935 3381 [407] 3864 3240 2980 2710 3570 3864 2500 1160 2443 3040 3864 3205 1875 3835 [421] 2841 2926 2983 3267 2890 2443 3352 3438 3120 4176 4148 3270 3523 3020 [435] 3620 3300 3779 3665 2869 3154 2810 2131 4233 3860 3920 3693 3580 3608 [449] 2330 2813 3068 3250 3175 2940 3267 3645 4205 1818 2983 2756 2642 3182 [463] 3125 3949 3466 3250 3750 3466 2557 3722 3551 3097 1763 3610 3551 2460 [477] 2358 3779 2671 3210 4375 4160 2841 3352 2790 2540 2890 3182 2557 3551 [491] 3324 2620 3068 3150 3636 3130 2784 3520 4205 3210 3523 3125 3097 3125 [505] 2930 3665 3693 3835 3636 4165 3035 2000 2740 3608 3750 4091 3267 3097 [519] 3381 3040 3779 3523 3551 3636 3154 2841 2960 3097 3600 3409 2880 2690 [533] 3090 3409 2159 3921 2940 3400 3239 3636 4347 2977 3636 2671 3949 3608 [547] 3710 2216 2500 1450 3296 3510 3950 4176 3636 3125 3845 3523 2784 3494 [561] 4233 3779 3720 3302 3693 2580 2597 2926 4505 3977 3182 2160 3402 3239 [575] 3608 3267 3995 3068 4065 3346 3154 2515 3182 2443 3807 3130 3722 2727 [589] 3011 793 3494 2386 3466 2760 4659 3608 2557 2671 3740 2614 3949 3630 [603] 2898 3940 701 3977 2472 3011 3352 3494 4091 3530 3010 3068 3636 3665 [617] 3430 3182 3300 2557 2550 3097 3154 3625 3330 4347 3296 3154 4050 3239 [631] 3608 2500 3438 3977 2780 2189 1760 2360 3608 3892 2870 2926 3905 3020 [645] 3779 3381 3280 3154 3494 3523 1861 3360 3977 3296 3523 3693 3921 653 [659] 3494 3500 2863 3494 3335 3750 3636 4190 3540 3381 3580 2813 3494 4065 [673] 3520 3551 4063 3210 3640 2770 3750 3695 3305 3892 3210 2851 3175 3949 [687] 2727 4262 3523 2955 3239 3665 2960 3353 3864 3125 3182 3438 3921 3728 [701] 2639 3324 2869 3352 3040 3381 4517 4375 4233 2983 1720 3555 3494 1932 [715] 3068 2756 3864 4690 3977 3864 3636 3154 3352 2443 2220 3240 3087 2955 [729] 2931 3350 3324 3438 3438 3430 3750 3520 3480 2955 3722 2898 3665 3324 [743] 2841 2690 3780 1165 3381 1900 2205 2790 3380 3060 3068 3280 2443 3750 [757] 3125 3799 3665 3409 3011 3930 3267 3693 3438 2557 2637 4034 3381 2960 [771] 3210 3572 3007 3977 3693 3011 2955 3296 3605 3722 3267 3608 2983 2386 [785] 682 2784 2528 1446 3480 2727 3267 2955 3068 2330 3154 3779 3714 3135 [799] 2188 4040 1750 4205 3335 1051 2898 3580 2860 3355 2699 2000 3154 3494 [813] 3210 521 2898 2386 3154 3220 3510 3608 3225 3494 2671 2841 3466 3949 [827] 3095 2784 3523 4240 3550 3296 3170 3693 3551 3100 1095 3523 3750 3182 [841] 2617 1421 2869 3939 3494 3608 3210 1920 2898 2813 2825 4460 2975 2810 [855] 3864 2869 3977 3665 3494 3180 3040 4091 3182 3296 4091 3779 410 3750 [869] 3010 3324 2926 3580 3636 3324 3324 3352 3352 3125 3448 3494 3608 3779 [883] 3860 3921 2642 3285 3835 3250 3210 3892 4546 3864 3210 3420 3154 2472 [897] 3665 3580 3381 3125 2880 2600 3835 3612 320 3244 3296 3693 3679 3450 [911] 3977 3154 4034 3160 682 3045 3409 3352 3665 3020 2557 4404 2840 3500 [925] 4063 4574 4006 3636 3494 3210 3409 2585 2555 3523 3155 3381 2983 3892 [939] 3324 3296 2443 3580 2784 3835 3721 3210 3892 2841 3068 2690 3025 4063 [953] 3551 3625 3210 2955 3230 3154 3438 4404 3615 4659 2841 4006 2869 470 [967] 3296 3728 3068 2585 3296 3438 3068 3438 3040 3523 2926 2301 3409 3324 [981] 3580 3722 2869 4119 3540 3352 4347 3097 4091 3381 2727 3864 2869 3523 [995] 2528 3324 2983 4148 2760 4480 3015 2898 3720 4176 3608 2756 2950 3210 [1009] 2926 3154 3090 3324 2955 4432 3864 3068 3665 3352 2815 3630 3040 3011 [1023] 3523 3494 3466 3921 3835 3381 3494 4100 3779 3864 3324 2479 3693 4148 [1037] 3210 3458 4347 3296 2720 3636 3665 3551 4830 3864 3210 3807 4460 3608 [1051] 3502 3239 3892 2995 3665 3381 3466 3892 3438 2955 3190 3060 2699 4375 [1065] 3722 3580 3494 3835 3750 3152 3523 3636 3523 2727 3040 3182 3494 2813 [1079] 2727 3875 3921 2642 2410 997 3409 2557 3523 3921 3486 3739 3864 3267 [1093] 3523 3011 1875 2330 3351 4006 2671 2784 2700 3324 3660 3125 3466 3197 [1107] 3296 3415 3320 3315 3068 4318 4034 3608 3352 > if (par2 == 'none') { + m <- ctree(as.formula(paste(colnames(x)[par1],' ~ .',sep='')),data = x) + } > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > if (par2 != 'none') { + m <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data = x) + if (par4=='yes') { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'10-Fold Cross Validation',3+2*par3,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'',1,TRUE) + a<-table.element(a,'Prediction (training)',par3+1,TRUE) + a<-table.element(a,'Prediction (testing)',par3+1,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Actual',1,TRUE) + for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE) + a<-table.element(a,'CV',1,TRUE) + for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE) + a<-table.element(a,'CV',1,TRUE) + a<-table.row.end(a) + for (i in 1:10) { + ind <- sample(2, nrow(x), replace=T, prob=c(0.9,0.1)) + m.ct <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data =x[ind==1,]) + if (i==1) { + m.ct.i.pred <- predict(m.ct, newdata=x[ind==1,]) + m.ct.i.actu <- x[ind==1,par1] + m.ct.x.pred <- predict(m.ct, newdata=x[ind==2,]) + m.ct.x.actu <- x[ind==2,par1] + } else { + m.ct.i.pred <- c(m.ct.i.pred,predict(m.ct, newdata=x[ind==1,])) + m.ct.i.actu <- c(m.ct.i.actu,x[ind==1,par1]) + m.ct.x.pred <- c(m.ct.x.pred,predict(m.ct, newdata=x[ind==2,])) + m.ct.x.actu <- c(m.ct.x.actu,x[ind==2,par1]) + } + } + print(m.ct.i.tab <- table(m.ct.i.actu,m.ct.i.pred)) + numer <- 0 + for (i in 1:par3) { + print(m.ct.i.tab[i,i] / sum(m.ct.i.tab[i,])) + numer <- numer + m.ct.i.tab[i,i] + } + print(m.ct.i.cp <- numer / sum(m.ct.i.tab)) + print(m.ct.x.tab <- table(m.ct.x.actu,m.ct.x.pred)) + numer <- 0 + for (i in 1:par3) { + print(m.ct.x.tab[i,i] / sum(m.ct.x.tab[i,])) + numer <- numer + m.ct.x.tab[i,i] + } + print(m.ct.x.cp <- numer / sum(m.ct.x.tab)) + for (i in 1:par3) { + a<-table.row.start(a) + a<-table.element(a,paste('C',i,sep=''),1,TRUE) + for (jjj in 1:par3) a<-table.element(a,m.ct.i.tab[i,jjj]) + a<-table.element(a,round(m.ct.i.tab[i,i]/sum(m.ct.i.tab[i,]),4)) + for (jjj in 1:par3) a<-table.element(a,m.ct.x.tab[i,jjj]) + a<-table.element(a,round(m.ct.x.tab[i,i]/sum(m.ct.x.tab[i,]),4)) + a<-table.row.end(a) + } + a<-table.row.start(a) + a<-table.element(a,'Overall',1,TRUE) + for (jjj in 1:par3) a<-table.element(a,'-') + a<-table.element(a,round(m.ct.i.cp,4)) + for (jjj in 1:par3) a<-table.element(a,'-') + a<-table.element(a,round(m.ct.x.cp,4)) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/1nz7b1323864782.tab") + } + } > m Conditional inference tree with 14 terminal nodes Response: grams Inputs: black, edu, smoke, gestate Number of observations: 1115 1) gestate <= 34; criterion = 1, statistic = 549.039 2) gestate <= 29; criterion = 1, statistic = 38.293 3) gestate <= 26; criterion = 1, statistic = 15.81 4)* weights = 12 3) gestate > 26 5)* weights = 9 2) gestate > 29 6) edu <= 11; criterion = 0.963, statistic = 6.76 7)* weights = 10 6) edu > 11 8)* weights = 23 1) gestate > 34 9) gestate <= 37; criterion = 1, statistic = 201.281 10) gestate <= 36; criterion = 0.996, statistic = 10.616 11)* weights = 57 10) gestate > 36 12)* weights = 57 9) gestate > 37 13) gestate <= 39; criterion = 1, statistic = 74.734 14) smoke <= 0; criterion = 1, statistic = 24.623 15) black <= 0; criterion = 0.998, statistic = 12.487 16)* weights = 110 15) black > 0 17)* weights = 185 14) smoke > 0 18)* weights = 98 13) gestate > 39 19) black <= 0; criterion = 1, statistic = 29.49 20) smoke <= 0; criterion = 1, statistic = 16.207 21)* weights = 220 20) smoke > 0 22) gestate <= 40; criterion = 0.993, statistic = 9.801 23)* weights = 39 22) gestate > 40 24)* weights = 20 19) black > 0 25) edu <= 11; criterion = 0.972, statistic = 7.228 26)* weights = 73 25) edu > 11 27)* weights = 202 > postscript(file="/var/wessaorg/rcomp/tmp/2yzn41323864782.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(m) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/3lw491323864782.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,par1] ~ as.factor(where(m)),main='Response by Terminal Node',xlab='Terminal Node',ylab='Response') > dev.off() null device 1 > if (par2 == 'none') { + forec <- predict(m) + result <- as.data.frame(cbind(x[,par1],forec,x[,par1]-forec)) + colnames(result) <- c('Actuals','Forecasts','Residuals') + print(result) + } Actuals Forecasts Residuals 1 2898 3201.7179 -303.7179487 2 994 590.0833 403.9166667 3 3977 3395.2909 581.7090909 4 3040 2942.5263 97.4736842 5 3523 3395.2909 127.7090909 6 3100 3201.7179 -101.7179487 7 3670 3235.6849 434.3150685 8 3097 3235.6849 -138.6849315 9 3040 3016.7551 23.2448980 10 3239 3016.7551 222.2448980 11 2955 3395.2909 -440.2909091 12 2200 3395.2909 -1195.2909091 13 3182 3235.6849 -53.6849315 14 3510 3632.9409 -122.9409091 15 3381 3016.7551 364.2448980 16 3530 3235.6849 294.3150685 17 2985 3016.7551 -31.7551020 18 3374 3209.0541 164.9459459 19 3765 3632.9409 132.0590909 20 2715 3209.0541 -494.0540541 21 3640 3209.0541 430.9459459 22 3040 3235.6849 -195.6849315 23 3670 3209.0541 460.9459459 24 2775 3016.7551 -241.7551020 25 2680 3209.0541 -529.0540541 26 2699 2942.5263 -243.5263158 27 2557 3016.7551 -459.7551020 28 3835 3235.6849 599.3150685 29 3381 3671.9500 -290.9500000 30 2766 2633.0877 132.9122807 31 3125 3235.6849 -110.6849315 32 2557 3016.7551 -459.7551020 33 2950 3209.0541 -259.0540541 34 3779 3632.9409 146.0590909 35 2760 3632.9409 -872.9409091 36 3011 3209.0541 -198.0540541 37 2415 2942.5263 -527.5263158 38 3125 3209.0541 -84.0540541 39 3154 3016.7551 137.2448980 40 3170 3395.2909 -225.2909091 41 2784 3632.9409 -848.9409091 42 2475 3235.6849 -760.6849315 43 2642 3235.6849 -593.6849315 44 4034 3201.7179 832.2820513 45 3988 3632.9409 355.0590909 46 3370 3235.6849 134.3150685 47 2955 3632.9409 -677.9409091 48 2443 3235.6849 -792.6849315 49 4000 3395.2909 604.7090909 50 3630 2942.5263 687.4736842 51 2920 3209.0541 -289.0540541 52 2727 2633.0877 93.9122807 53 3800 3395.2909 404.7090909 54 3771 3632.9409 138.0590909 55 3075 2942.5263 132.4736842 56 2415 2942.5263 -527.5263158 57 2890 2942.5263 -52.5263158 58 3892 2942.5263 949.4736842 59 3125 3235.6849 -110.6849315 60 3892 3671.9500 220.0500000 61 3296 3395.2909 -99.2909091 62 3120 3395.2909 -275.2909091 63 1710 3016.7551 -1306.7551020 64 3068 3235.6849 -167.6849315 65 3779 3235.6849 543.3150685 66 2898 3235.6849 -337.6849315 67 3381 3632.9409 -251.9409091 68 3324 3235.6849 88.3150685 69 2727 3016.7551 -289.7551020 70 2642 2942.5263 -300.5263158 71 2244 2178.9000 65.1000000 72 2528 2633.0877 -105.0877193 73 3352 3016.7551 335.2448980 74 3551 3016.7551 534.2448980 75 3267 3632.9409 -365.9409091 76 3210 3235.6849 -25.6849315 77 3395 3235.6849 159.3150685 78 3068 3201.7179 -133.7179487 79 3420 3016.7551 403.2448980 80 2985 3209.0541 -224.0540541 81 2147 2178.9000 -31.9000000 82 2863 3209.0541 -346.0540541 83 3125 3209.0541 -84.0540541 84 4290 3632.9409 657.0590909 85 2585 2178.9000 406.1000000 86 2869 3209.0541 -340.0540541 87 3097 2633.0877 463.9122807 88 2756 3209.0541 -453.0540541 89 3125 2633.0877 491.9122807 90 3921 3671.9500 249.0500000 91 3438 3235.6849 202.3150685 92 2980 3235.6849 -255.6849315 93 2869 2633.0877 235.9122807 94 4830 3671.9500 1158.0500000 95 3409 3201.7179 207.2820513 96 3807 3235.6849 571.3150685 97 3220 3016.7551 203.2448980 98 2756 3016.7551 -260.7551020 99 3097 3016.7551 80.2448980 100 3615 3209.0541 405.9459459 101 2926 3632.9409 -706.9409091 102 4050 3235.6849 814.3150685 103 2840 3209.0541 -369.0540541 104 2671 3235.6849 -564.6849315 105 3540 3209.0541 330.9459459 106 2900 3209.0541 -309.0540541 107 3779 3632.9409 146.0590909 108 2740 3209.0541 -469.0540541 109 2320 2633.0877 -313.0877193 110 2905 3209.0541 -304.0540541 111 3070 3016.7551 53.2448980 112 3352 3671.9500 -319.9500000 113 2581 2633.0877 -52.0877193 114 4091 3235.6849 855.3150685 115 3235 3235.6849 -0.6849315 116 3722 3235.6849 486.3150685 117 3438 3235.6849 202.3150685 118 3475 3632.9409 -157.9409091 119 3438 2178.9000 1259.1000000 120 3409 3209.0541 199.9459459 121 3040 3395.2909 -355.2909091 122 3800 3395.2909 404.7090909 123 3182 2633.0877 548.9122807 124 3494 3016.7551 477.2448980 125 2131 2633.0877 -502.0877193 126 1790 2633.0877 -843.0877193 127 3864 3632.9409 231.0590909 128 2955 3209.0541 -254.0540541 129 3949 3632.9409 316.0590909 130 2273 2633.0877 -360.0877193 131 2620 3209.0541 -589.0540541 132 2205 2633.0877 -428.0877193 133 2860 3235.6849 -375.6849315 134 3381 3235.6849 145.3150685 135 3190 3235.6849 -45.6849315 136 3267 3235.6849 31.3150685 137 3665 3201.7179 463.2820513 138 3160 3209.0541 -49.0540541 139 3216 3209.0541 6.9459459 140 2415 2942.5263 -527.5263158 141 3200 3235.6849 -35.6849315 142 3352 3632.9409 -280.9409091 143 3470 3235.6849 234.3150685 144 3580 3209.0541 370.9459459 145 3750 3632.9409 117.0590909 146 3035 3235.6849 -200.6849315 147 3438 3209.0541 228.9459459 148 3068 3209.0541 -141.0540541 149 3438 3395.2909 42.7090909 150 3750 3016.7551 733.2448980 151 3409 3235.6849 173.3150685 152 1790 3016.7551 -1226.7551020 153 3170 3235.6849 -65.6849315 154 3800 3395.2909 404.7090909 155 2825 3235.6849 -410.6849315 156 2756 2942.5263 -186.5263158 157 3040 3235.6849 -195.6849315 158 2585 3016.7551 -431.7551020 159 2756 3016.7551 -260.7551020 160 3523 3671.9500 -148.9500000 161 909 590.0833 318.9166667 162 3779 3201.7179 577.2820513 163 2415 3016.7551 -601.7551020 164 3665 3632.9409 32.0590909 165 3295 3209.0541 85.9459459 166 2699 3016.7551 -317.7551020 167 3267 3016.7551 250.2448980 168 3460 3209.0541 250.9459459 169 1250 1146.5556 103.4444444 170 3730 3209.0541 520.9459459 171 3523 3235.6849 287.3150685 172 3011 3235.6849 -224.6849315 173 3740 3016.7551 723.2448980 174 3949 3235.6849 713.3150685 175 2330 2942.5263 -612.5263158 176 3182 3632.9409 -450.9409091 177 3040 3209.0541 -169.0540541 178 2557 2942.5263 -385.5263158 179 3147 3395.2909 -248.2909091 180 2671 2633.0877 37.9122807 181 2557 3235.6849 -678.6849315 182 3200 3209.0541 -9.0540541 183 3203 3235.6849 -32.6849315 184 2756 2633.0877 122.9122807 185 1619 2178.9000 -559.9000000 186 2330 2633.0877 -303.0877193 187 3381 3016.7551 364.2448980 188 2983 3016.7551 -33.7551020 189 2600 3016.7551 -416.7551020 190 3585 3235.6849 349.3150685 191 3693 3632.9409 60.0590909 192 3020 3209.0541 -189.0540541 193 1800 2178.9000 -378.9000000 194 2290 3209.0541 -919.0540541 195 3750 3235.6849 514.3150685 196 3551 3632.9409 -81.9409091 197 3779 2942.5263 836.4736842 198 3125 3395.2909 -270.2909091 199 3470 3209.0541 260.9459459 200 3210 3235.6849 -25.6849315 201 3100 3235.6849 -135.6849315 202 3182 3016.7551 165.2448980 203 3355 2942.5263 412.4736842 204 2750 3235.6849 -485.6849315 205 4006 3209.0541 796.9459459 206 3438 3209.0541 228.9459459 207 3068 3632.9409 -564.9409091 208 2650 3235.6849 -585.6849315 209 3005 3209.0541 -204.0540541 210 2993 3209.0541 -216.0540541 211 3494 3235.6849 258.3150685 212 2795 3235.6849 -440.6849315 213 3260 3235.6849 24.3150685 214 2959 3235.6849 -276.6849315 215 4020 3209.0541 810.9459459 216 2949 2942.5263 6.4736842 217 3900 3209.0541 690.9459459 218 2926 3016.7551 -90.7551020 219 3239 3016.7551 222.2448980 220 3409 3235.6849 173.3150685 221 3630 3235.6849 394.3150685 222 3466 3235.6849 230.3150685 223 2330 3235.6849 -905.6849315 224 3949 3201.7179 747.2820513 225 2813 3016.7551 -203.7551020 226 3137 2942.5263 194.4736842 227 3296 3632.9409 -336.9409091 228 3267 3235.6849 31.3150685 229 2940 3209.0541 -269.0540541 230 3323 3235.6849 87.3150685 231 3409 2942.5263 466.4736842 232 3296 3235.6849 60.3150685 233 3807 3201.7179 605.2820513 234 3296 3209.0541 86.9459459 235 2983 3016.7551 -33.7551020 236 1680 2178.9000 -498.9000000 237 3551 3209.0541 341.9459459 238 2385 2942.5263 -557.5263158 239 3356 3235.6849 120.3150685 240 3466 3235.6849 230.3150685 241 2699 3632.9409 -933.9409091 242 4233 3632.9409 600.0590909 243 2890 3235.6849 -345.6849315 244 2443 2942.5263 -499.5263158 245 3580 3016.7551 563.2448980 246 3551 3671.9500 -120.9500000 247 2030 2178.9000 -148.9000000 248 3125 2633.0877 491.9122807 249 2700 2942.5263 -242.5263158 250 2216 2178.9000 37.1000000 251 3182 3201.7179 -19.7179487 252 3040 3016.7551 23.2448980 253 3011 3016.7551 -5.7551020 254 2800 3016.7551 -216.7551020 255 3425 3209.0541 215.9459459 256 2760 2942.5263 -182.5263158 257 3665 3235.6849 429.3150685 258 1733 2633.0877 -900.0877193 259 2850 3235.6849 -385.6849315 260 2485 3016.7551 -531.7551020 261 3381 3235.6849 145.3150685 262 3065 3209.0541 -144.0540541 263 3167 3209.0541 -42.0540541 264 2635 3209.0541 -574.0540541 265 3665 3201.7179 463.2820513 266 3300 3016.7551 283.2448980 267 3636 3632.9409 3.0590909 268 2813 3016.7551 -203.7551020 269 3920 3235.6849 684.3150685 270 2671 3209.0541 -538.0540541 271 3182 3235.6849 -53.6849315 272 3409 3235.6849 173.3150685 273 3600 2942.5263 657.4736842 274 2557 3235.6849 -678.6849315 275 3890 3235.6849 654.3150685 276 4148 3209.0541 938.9459459 277 2030 2178.9000 -148.9000000 278 2890 3016.7551 -126.7551020 279 2373 2633.0877 -260.0877193 280 2977 3235.6849 -258.6849315 281 2840 3016.7551 -176.7551020 282 2658 3235.6849 -577.6849315 283 2983 3201.7179 -218.7179487 284 3700 3209.0541 490.9459459 285 3608 3201.7179 406.2820513 286 3559 3632.9409 -73.9409091 287 3540 3415.2723 124.7277228 288 3440 3016.7551 423.2448980 289 3324 3209.0541 114.9459459 290 3494 3209.0541 284.9459459 291 3636 3415.2723 220.7277228 292 2813 3016.7551 -203.7551020 293 3239 3209.0541 29.9459459 294 3068 3016.7551 51.2448980 295 3210 3415.2723 -205.2722772 296 3892 3632.9409 259.0590909 297 3864 3395.2909 468.7090909 298 2784 3415.2723 -631.2722772 299 3551 3632.9409 -81.9409091 300 3608 3415.2723 192.7277228 301 3130 3209.0541 -79.0540541 302 3466 3201.7179 264.2820513 303 3608 3415.2723 192.7277228 304 2225 2942.5263 -717.5263158 305 2585 2633.0877 -48.0877193 306 2917 3415.2723 -498.2722772 307 2415 2633.0877 -218.0877193 308 4489 3632.9409 856.0590909 309 3977 3632.9409 344.0590909 310 3080 3209.0541 -129.0540541 311 3722 3415.2723 306.7277228 312 3097 3201.7179 -104.7179487 313 2960 3415.2723 -455.2722772 314 3807 3395.2909 411.7090909 315 3210 3632.9409 -422.9409091 316 3210 3415.2723 -205.2722772 317 284 590.0833 -306.0833333 318 3805 3415.2723 389.7277228 319 3267 3671.9500 -404.9500000 320 3722 3671.9500 50.0500000 321 3864 3632.9409 231.0590909 322 4318 3671.9500 646.0500000 323 3807 3671.9500 135.0500000 324 455 590.0833 -135.0833333 325 3722 3632.9409 89.0590909 326 3324 3209.0541 114.9459459 327 3665 3415.2723 249.7277228 328 3240 3415.2723 -175.2722772 329 3182 3415.2723 -233.2722772 330 3040 3415.2723 -375.2722772 331 3665 3201.7179 463.2820513 332 4148 3632.9409 515.0590909 333 3636 3395.2909 240.7090909 334 1591 2014.8696 -423.8695652 335 4135 3395.2909 739.7090909 336 3352 3632.9409 -280.9409091 337 3011 3016.7551 -5.7551020 338 2301 2014.8696 286.1304348 339 2740 3415.2723 -675.2722772 340 2450 3415.2723 -965.2722772 341 2472 3201.7179 -729.7179487 342 4148 3632.9409 515.0590909 343 3352 3632.9409 -280.9409091 344 4034 3395.2909 638.7090909 345 3807 3632.9409 174.0590909 346 2926 3395.2909 -469.2909091 347 2699 3415.2723 -716.2722772 348 2990 3209.0541 -219.0540541 349 2898 3395.2909 -497.2909091 350 3864 3671.9500 192.0500000 351 2970 3415.2723 -445.2722772 352 3655 3415.2723 239.7277228 353 2890 2633.0877 256.9122807 354 3980 3016.7551 963.2448980 355 2756 3201.7179 -445.7179487 356 4148 3632.9409 515.0590909 357 3494 3632.9409 -138.9409091 358 2830 2942.5263 -112.5263158 359 2983 3209.0541 -226.0540541 360 3430 3415.2723 14.7277228 361 2159 2633.0877 -474.0877193 362 3125 3671.9500 -546.9500000 363 3239 3632.9409 -393.9409091 364 4262 3395.2909 866.7090909 365 3459 3209.0541 249.9459459 366 1190 1146.5556 43.4444444 367 3182 3632.9409 -450.9409091 368 3352 2633.0877 718.9122807 369 2898 3395.2909 -497.2909091 370 3154 3632.9409 -478.9409091 371 1165 1146.5556 18.4444444 372 2102 3016.7551 -914.7551020 373 3110 3415.2723 -305.2722772 374 2367 2014.8696 352.1304348 375 3892 3209.0541 682.9459459 376 3011 3415.2723 -404.2722772 377 3125 2942.5263 182.4736842 378 3490 3415.2723 74.7277228 379 2955 3016.7551 -61.7551020 380 3977 3415.2723 561.7277228 381 3665 3632.9409 32.0590909 382 3355 3209.0541 145.9459459 383 3182 3415.2723 -233.2722772 384 3239 3632.9409 -393.9409091 385 2955 3415.2723 -460.2722772 386 3182 3209.0541 -27.0540541 387 2784 3209.0541 -425.0540541 388 3802 3632.9409 169.0590909 389 3270 3415.2723 -145.2722772 390 2225 3209.0541 -984.0540541 391 3182 2633.0877 548.9122807 392 3210 3415.2723 -205.2722772 393 3779 3395.2909 383.7090909 394 2415 2633.0877 -218.0877193 395 4675 3632.9409 1042.0590909 396 3352 3209.0541 142.9459459 397 2900 3016.7551 -116.7551020 398 3864 3632.9409 231.0590909 399 2841 3415.2723 -574.2722772 400 2557 3632.9409 -1075.9409091 401 3494 3632.9409 -138.9409091 402 3438 3415.2723 22.7277228 403 2445 2633.0877 -188.0877193 404 3714 3415.2723 298.7277228 405 2935 3415.2723 -480.2722772 406 3381 3632.9409 -251.9409091 407 3864 3632.9409 231.0590909 408 3240 3415.2723 -175.2722772 409 2980 3415.2723 -435.2722772 410 2710 3016.7551 -306.7551020 411 3570 3415.2723 154.7277228 412 3864 3415.2723 448.7277228 413 2500 3016.7551 -516.7551020 414 1160 1146.5556 13.4444444 415 2443 3209.0541 -766.0540541 416 3040 3415.2723 -375.2722772 417 3864 3632.9409 231.0590909 418 3205 3415.2723 -210.2722772 419 1875 2014.8696 -139.8695652 420 3835 3209.0541 625.9459459 421 2841 3016.7551 -175.7551020 422 2926 2633.0877 292.9122807 423 2983 3415.2723 -432.2722772 424 3267 3201.7179 65.2820513 425 2890 3016.7551 -126.7551020 426 2443 3209.0541 -766.0540541 427 3352 3395.2909 -43.2909091 428 3438 3395.2909 42.7090909 429 3120 3016.7551 103.2448980 430 4176 3632.9409 543.0590909 431 4148 3632.9409 515.0590909 432 3270 3415.2723 -145.2722772 433 3523 3632.9409 -109.9409091 434 3020 3632.9409 -612.9409091 435 3620 3395.2909 224.7090909 436 3300 3415.2723 -115.2722772 437 3779 2942.5263 836.4736842 438 3665 3671.9500 -6.9500000 439 2869 3201.7179 -332.7179487 440 3154 3209.0541 -55.0540541 441 2810 3016.7551 -206.7551020 442 2131 2014.8696 116.1304348 443 4233 3632.9409 600.0590909 444 3860 3395.2909 464.7090909 445 3920 3415.2723 504.7277228 446 3693 3632.9409 60.0590909 447 3580 3415.2723 164.7277228 448 3608 3415.2723 192.7277228 449 2330 2942.5263 -612.5263158 450 2813 3201.7179 -388.7179487 451 3068 2942.5263 125.4736842 452 3250 3209.0541 40.9459459 453 3175 3209.0541 -34.0540541 454 2940 3209.0541 -269.0540541 455 3267 3632.9409 -365.9409091 456 3645 3415.2723 229.7277228 457 4205 3632.9409 572.0590909 458 1818 2014.8696 -196.8695652 459 2983 2633.0877 349.9122807 460 2756 3415.2723 -659.2722772 461 2642 3395.2909 -753.2909091 462 3182 3632.9409 -450.9409091 463 3125 3632.9409 -507.9409091 464 3949 3395.2909 553.7090909 465 3466 3395.2909 70.7090909 466 3250 3016.7551 233.2448980 467 3750 3632.9409 117.0590909 468 3466 3632.9409 -166.9409091 469 2557 3201.7179 -644.7179487 470 3722 3632.9409 89.0590909 471 3551 3632.9409 -81.9409091 472 3097 3209.0541 -112.0540541 473 1763 2633.0877 -870.0877193 474 3610 3415.2723 194.7277228 475 3551 3209.0541 341.9459459 476 2460 3201.7179 -741.7179487 477 2358 2633.0877 -275.0877193 478 3779 3632.9409 146.0590909 479 2671 3415.2723 -744.2722772 480 3210 3209.0541 0.9459459 481 4375 3632.9409 742.0590909 482 4160 3415.2723 744.7277228 483 2841 3415.2723 -574.2722772 484 3352 3632.9409 -280.9409091 485 2790 3016.7551 -226.7551020 486 2540 3415.2723 -875.2722772 487 2890 3209.0541 -319.0540541 488 3182 3415.2723 -233.2722772 489 2557 3395.2909 -838.2909091 490 3551 3632.9409 -81.9409091 491 3324 3415.2723 -91.2722772 492 2620 3209.0541 -589.0540541 493 3068 3209.0541 -141.0540541 494 3150 3209.0541 -59.0540541 495 3636 3415.2723 220.7277228 496 3130 3415.2723 -285.2722772 497 2784 2633.0877 150.9122807 498 3520 3415.2723 104.7277228 499 4205 3632.9409 572.0590909 500 3210 3209.0541 0.9459459 501 3523 3632.9409 -109.9409091 502 3125 3209.0541 -84.0540541 503 3097 3016.7551 80.2448980 504 3125 3016.7551 108.2448980 505 2930 3415.2723 -485.2722772 506 3665 3395.2909 269.7090909 507 3693 3632.9409 60.0590909 508 3835 3415.2723 419.7277228 509 3636 3415.2723 220.7277228 510 4165 3415.2723 749.7277228 511 3035 3016.7551 18.2448980 512 2000 2942.5263 -942.5263158 513 2740 3016.7551 -276.7551020 514 3608 3395.2909 212.7090909 515 3750 3395.2909 354.7090909 516 4091 3632.9409 458.0590909 517 3267 3415.2723 -148.2722772 518 3097 3415.2723 -318.2722772 519 3381 3415.2723 -34.2722772 520 3040 3209.0541 -169.0540541 521 3779 3632.9409 146.0590909 522 3523 3201.7179 321.2820513 523 3551 3632.9409 -81.9409091 524 3636 3632.9409 3.0590909 525 3154 3016.7551 137.2448980 526 2841 3016.7551 -175.7551020 527 2960 2942.5263 17.4736842 528 3097 3415.2723 -318.2722772 529 3600 3209.0541 390.9459459 530 3409 3016.7551 392.2448980 531 2880 3209.0541 -329.0540541 532 2690 3209.0541 -519.0540541 533 3090 3632.9409 -542.9409091 534 3409 3632.9409 -223.9409091 535 2159 2014.8696 144.1304348 536 3921 3415.2723 505.7277228 537 2940 3632.9409 -692.9409091 538 3400 3415.2723 -15.2722772 539 3239 3415.2723 -176.2722772 540 3636 3632.9409 3.0590909 541 4347 3415.2723 931.7277228 542 2977 3016.7551 -39.7551020 543 3636 3415.2723 220.7277228 544 2671 3395.2909 -724.2909091 545 3949 3632.9409 316.0590909 546 3608 2942.5263 665.4736842 547 3710 3415.2723 294.7277228 548 2216 2633.0877 -417.0877193 549 2500 2633.0877 -133.0877193 550 1450 1146.5556 303.4444444 551 3296 3632.9409 -336.9409091 552 3510 3209.0541 300.9459459 553 3950 3415.2723 534.7277228 554 4176 3395.2909 780.7090909 555 3636 3415.2723 220.7277228 556 3125 3415.2723 -290.2722772 557 3845 3395.2909 449.7090909 558 3523 3415.2723 107.7277228 559 2784 3395.2909 -611.2909091 560 3494 3632.9409 -138.9409091 561 4233 3632.9409 600.0590909 562 3779 3632.9409 146.0590909 563 3720 3209.0541 510.9459459 564 3302 3209.0541 92.9459459 565 3693 3632.9409 60.0590909 566 2580 2633.0877 -53.0877193 567 2597 3209.0541 -612.0540541 568 2926 3395.2909 -469.2909091 569 4505 3632.9409 872.0590909 570 3977 3632.9409 344.0590909 571 3182 3415.2723 -233.2722772 572 2160 2633.0877 -473.0877193 573 3402 3209.0541 192.9459459 574 3239 3632.9409 -393.9409091 575 3608 3395.2909 212.7090909 576 3267 3395.2909 -128.2909091 577 3995 3395.2909 599.7090909 578 3068 3016.7551 51.2448980 579 4065 3632.9409 432.0590909 580 3346 3209.0541 136.9459459 581 3154 3632.9409 -478.9409091 582 2515 2014.8696 500.1304348 583 3182 3395.2909 -213.2909091 584 2443 2633.0877 -190.0877193 585 3807 3632.9409 174.0590909 586 3130 3016.7551 113.2448980 587 3722 3201.7179 520.2820513 588 2727 3415.2723 -688.2722772 589 3011 3201.7179 -190.7179487 590 793 1146.5556 -353.5555556 591 3494 3415.2723 78.7277228 592 2386 2633.0877 -247.0877193 593 3466 3632.9409 -166.9409091 594 2760 3395.2909 -635.2909091 595 4659 3632.9409 1026.0590909 596 3608 3395.2909 212.7090909 597 2557 3209.0541 -652.0540541 598 2671 3632.9409 -961.9409091 599 3740 3395.2909 344.7090909 600 2614 2014.8696 599.1304348 601 3949 3632.9409 316.0590909 602 3630 3209.0541 420.9459459 603 2898 3395.2909 -497.2909091 604 3940 3415.2723 524.7277228 605 701 590.0833 110.9166667 606 3977 3632.9409 344.0590909 607 2472 2014.8696 457.1304348 608 3011 2942.5263 68.4736842 609 3352 3209.0541 142.9459459 610 3494 3632.9409 -138.9409091 611 4091 3632.9409 458.0590909 612 3530 3209.0541 320.9459459 613 3010 2633.0877 376.9122807 614 3068 3632.9409 -564.9409091 615 3636 3415.2723 220.7277228 616 3665 3415.2723 249.7277228 617 3430 3415.2723 14.7277228 618 3182 3209.0541 -27.0540541 619 3300 3209.0541 90.9459459 620 2557 2942.5263 -385.5263158 621 2550 2942.5263 -392.5263158 622 3097 3201.7179 -104.7179487 623 3154 3632.9409 -478.9409091 624 3625 3415.2723 209.7277228 625 3330 3415.2723 -85.2722772 626 4347 3632.9409 714.0590909 627 3296 3395.2909 -99.2909091 628 3154 3209.0541 -55.0540541 629 4050 3415.2723 634.7277228 630 3239 2633.0877 605.9122807 631 3608 3632.9409 -24.9409091 632 2500 2633.0877 -133.0877193 633 3438 3632.9409 -194.9409091 634 3977 3632.9409 344.0590909 635 2780 3016.7551 -236.7551020 636 2189 2633.0877 -444.0877193 637 1760 2014.8696 -254.8695652 638 2360 2633.0877 -273.0877193 639 3608 3201.7179 406.2820513 640 3892 3415.2723 476.7277228 641 2870 3209.0541 -339.0540541 642 2926 3016.7551 -90.7551020 643 3905 3632.9409 272.0590909 644 3020 3209.0541 -189.0540541 645 3779 3395.2909 383.7090909 646 3381 3632.9409 -251.9409091 647 3280 3209.0541 70.9459459 648 3154 3415.2723 -261.2722772 649 3494 3632.9409 -138.9409091 650 3523 3415.2723 107.7277228 651 1861 2014.8696 -153.8695652 652 3360 3415.2723 -55.2722772 653 3977 3632.9409 344.0590909 654 3296 3016.7551 279.2448980 655 3523 3395.2909 127.7090909 656 3693 3415.2723 277.7277228 657 3921 3209.0541 711.9459459 658 653 590.0833 62.9166667 659 3494 3415.2723 78.7277228 660 3500 3415.2723 84.7277228 661 2863 3209.0541 -346.0540541 662 3494 3632.9409 -138.9409091 663 3335 3209.0541 125.9459459 664 3750 3632.9409 117.0590909 665 3636 3632.9409 3.0590909 666 4190 3016.7551 1173.2448980 667 3540 3209.0541 330.9459459 668 3381 3415.2723 -34.2722772 669 3580 3209.0541 370.9459459 670 2813 2014.8696 798.1304348 671 3494 3632.9409 -138.9409091 672 4065 3415.2723 649.7277228 673 3520 3395.2909 124.7090909 674 3551 3209.0541 341.9459459 675 4063 3632.9409 430.0590909 676 3210 3209.0541 0.9459459 677 3640 3415.2723 224.7277228 678 2770 3201.7179 -431.7179487 679 3750 3415.2723 334.7277228 680 3695 3415.2723 279.7277228 681 3305 3415.2723 -110.2722772 682 3892 3632.9409 259.0590909 683 3210 3671.9500 -461.9500000 684 2851 3415.2723 -564.2722772 685 3175 3209.0541 -34.0540541 686 3949 3415.2723 533.7277228 687 2727 2942.5263 -215.5263158 688 4262 3632.9409 629.0590909 689 3523 3415.2723 107.7277228 690 2955 2633.0877 321.9122807 691 3239 3632.9409 -393.9409091 692 3665 3632.9409 32.0590909 693 2960 3016.7551 -56.7551020 694 3353 3415.2723 -62.2722772 695 3864 3201.7179 662.2820513 696 3125 3209.0541 -84.0540541 697 3182 3415.2723 -233.2722772 698 3438 3632.9409 -194.9409091 699 3921 3632.9409 288.0590909 700 3728 3209.0541 518.9459459 701 2639 3395.2909 -756.2909091 702 3324 3415.2723 -91.2722772 703 2869 3415.2723 -546.2722772 704 3352 3016.7551 335.2448980 705 3040 3209.0541 -169.0540541 706 3381 3415.2723 -34.2722772 707 4517 3415.2723 1101.7277228 708 4375 3632.9409 742.0590909 709 4233 3632.9409 600.0590909 710 2983 3632.9409 -649.9409091 711 1720 2014.8696 -294.8695652 712 3555 3016.7551 538.2448980 713 3494 3632.9409 -138.9409091 714 1932 3209.0541 -1277.0540541 715 3068 3016.7551 51.2448980 716 2756 3201.7179 -445.7179487 717 3864 3415.2723 448.7277228 718 4690 3415.2723 1274.7277228 719 3977 3632.9409 344.0590909 720 3864 3395.2909 468.7090909 721 3636 3671.9500 -35.9500000 722 3154 2633.0877 520.9122807 723 3352 3395.2909 -43.2909091 724 2443 3209.0541 -766.0540541 725 2220 2014.8696 205.1304348 726 3240 3209.0541 30.9459459 727 3087 3415.2723 -328.2722772 728 2955 2942.5263 12.4736842 729 2931 3209.0541 -278.0540541 730 3350 3209.0541 140.9459459 731 3324 3671.9500 -347.9500000 732 3438 3632.9409 -194.9409091 733 3438 3415.2723 22.7277228 734 3430 3209.0541 220.9459459 735 3750 3395.2909 354.7090909 736 3520 3415.2723 104.7277228 737 3480 3415.2723 64.7277228 738 2955 3201.7179 -246.7179487 739 3722 3632.9409 89.0590909 740 2898 3632.9409 -734.9409091 741 3665 3395.2909 269.7090909 742 3324 3415.2723 -91.2722772 743 2841 3415.2723 -574.2722772 744 2690 3415.2723 -725.2722772 745 3780 3415.2723 364.7277228 746 1165 1146.5556 18.4444444 747 3381 3415.2723 -34.2722772 748 1900 3201.7179 -1301.7179487 749 2205 3209.0541 -1004.0540541 750 2790 3415.2723 -625.2722772 751 3380 3209.0541 170.9459459 752 3060 3209.0541 -149.0540541 753 3068 3632.9409 -564.9409091 754 3280 3209.0541 70.9459459 755 2443 3016.7551 -573.7551020 756 3750 3395.2909 354.7090909 757 3125 3395.2909 -270.2909091 758 3799 2942.5263 856.4736842 759 3665 3209.0541 455.9459459 760 3409 3395.2909 13.7090909 761 3011 3632.9409 -621.9409091 762 3930 3415.2723 514.7277228 763 3267 3632.9409 -365.9409091 764 3693 3016.7551 676.2448980 765 3438 3415.2723 22.7277228 766 2557 2942.5263 -385.5263158 767 2637 2942.5263 -305.5263158 768 4034 3415.2723 618.7277228 769 3381 3395.2909 -14.2909091 770 2960 2633.0877 326.9122807 771 3210 3395.2909 -185.2909091 772 3572 3209.0541 362.9459459 773 3007 3016.7551 -9.7551020 774 3977 3632.9409 344.0590909 775 3693 3671.9500 21.0500000 776 3011 3415.2723 -404.2722772 777 2955 3395.2909 -440.2909091 778 3296 3415.2723 -119.2722772 779 3605 3415.2723 189.7277228 780 3722 3209.0541 512.9459459 781 3267 3016.7551 250.2448980 782 3608 3632.9409 -24.9409091 783 2983 3209.0541 -226.0540541 784 2386 2014.8696 371.1304348 785 682 590.0833 91.9166667 786 2784 2633.0877 150.9122807 787 2528 3016.7551 -488.7551020 788 1446 2014.8696 -568.8695652 789 3480 3209.0541 270.9459459 790 2727 3209.0541 -482.0540541 791 3267 3415.2723 -148.2722772 792 2955 3209.0541 -254.0540541 793 3068 3415.2723 -347.2722772 794 2330 2014.8696 315.1304348 795 3154 3201.7179 -47.7179487 796 3779 2942.5263 836.4736842 797 3714 3632.9409 81.0590909 798 3135 3209.0541 -74.0540541 799 2188 2633.0877 -445.0877193 800 4040 3209.0541 830.9459459 801 1750 2014.8696 -264.8695652 802 4205 3632.9409 572.0590909 803 3335 3209.0541 125.9459459 804 1051 1146.5556 -95.5555556 805 2898 3632.9409 -734.9409091 806 3580 3632.9409 -52.9409091 807 2860 2633.0877 226.9122807 808 3355 3415.2723 -60.2722772 809 2699 3201.7179 -502.7179487 810 2000 3016.7551 -1016.7551020 811 3154 3632.9409 -478.9409091 812 3494 3395.2909 98.7090909 813 3210 2942.5263 267.4736842 814 521 590.0833 -69.0833333 815 2898 3209.0541 -311.0540541 816 2386 3016.7551 -630.7551020 817 3154 3209.0541 -55.0540541 818 3220 3209.0541 10.9459459 819 3510 3415.2723 94.7277228 820 3608 3632.9409 -24.9409091 821 3225 3209.0541 15.9459459 822 3494 3632.9409 -138.9409091 823 2671 3209.0541 -538.0540541 824 2841 3201.7179 -360.7179487 825 3466 3415.2723 50.7277228 826 3949 3632.9409 316.0590909 827 3095 3016.7551 78.2448980 828 2784 3415.2723 -631.2722772 829 3523 2942.5263 580.4736842 830 4240 3415.2723 824.7277228 831 3550 3415.2723 134.7277228 832 3296 3209.0541 86.9459459 833 3170 3415.2723 -245.2722772 834 3693 2942.5263 750.4736842 835 3551 3209.0541 341.9459459 836 3100 3209.0541 -109.0540541 837 1095 1146.5556 -51.5555556 838 3523 3415.2723 107.7277228 839 3750 3632.9409 117.0590909 840 3182 3415.2723 -233.2722772 841 2617 3016.7551 -399.7551020 842 1421 2014.8696 -593.8695652 843 2869 3016.7551 -147.7551020 844 3939 3632.9409 306.0590909 845 3494 3415.2723 78.7277228 846 3608 3415.2723 192.7277228 847 3210 3209.0541 0.9459459 848 1920 2014.8696 -94.8695652 849 2898 3395.2909 -497.2909091 850 2813 3415.2723 -602.2722772 851 2825 3209.0541 -384.0540541 852 4460 3415.2723 1044.7277228 853 2975 3209.0541 -234.0540541 854 2810 3016.7551 -206.7551020 855 3864 3632.9409 231.0590909 856 2869 3016.7551 -147.7551020 857 3977 3415.2723 561.7277228 858 3665 3632.9409 32.0590909 859 3494 3415.2723 78.7277228 860 3180 3415.2723 -235.2722772 861 3040 3415.2723 -375.2722772 862 4091 3395.2909 695.7090909 863 3182 3209.0541 -27.0540541 864 3296 3415.2723 -119.2722772 865 4091 3415.2723 675.7277228 866 3779 3632.9409 146.0590909 867 410 590.0833 -180.0833333 868 3750 3415.2723 334.7277228 869 3010 3415.2723 -405.2722772 870 3324 3415.2723 -91.2722772 871 2926 2942.5263 -16.5263158 872 3580 3209.0541 370.9459459 873 3636 3415.2723 220.7277228 874 3324 3415.2723 -91.2722772 875 3324 3209.0541 114.9459459 876 3352 3415.2723 -63.2722772 877 3352 3632.9409 -280.9409091 878 3125 2942.5263 182.4736842 879 3448 3632.9409 -184.9409091 880 3494 3632.9409 -138.9409091 881 3608 3415.2723 192.7277228 882 3779 3415.2723 363.7277228 883 3860 3415.2723 444.7277228 884 3921 3632.9409 288.0590909 885 2642 2633.0877 8.9122807 886 3285 3415.2723 -130.2722772 887 3835 3415.2723 419.7277228 888 3250 3415.2723 -165.2722772 889 3210 3632.9409 -422.9409091 890 3892 3415.2723 476.7277228 891 4546 3632.9409 913.0590909 892 3864 3395.2909 468.7090909 893 3210 3209.0541 0.9459459 894 3420 3415.2723 4.7277228 895 3154 3415.2723 -261.2722772 896 2472 3415.2723 -943.2722772 897 3665 3632.9409 32.0590909 898 3580 3395.2909 184.7090909 899 3381 3632.9409 -251.9409091 900 3125 3632.9409 -507.9409091 901 2880 2633.0877 246.9122807 902 2600 3209.0541 -609.0540541 903 3835 3632.9409 202.0590909 904 3612 3415.2723 196.7277228 905 320 590.0833 -270.0833333 906 3244 3415.2723 -171.2722772 907 3296 3415.2723 -119.2722772 908 3693 3632.9409 60.0590909 909 3679 3632.9409 46.0590909 910 3450 3415.2723 34.7277228 911 3977 3209.0541 767.9459459 912 3154 3209.0541 -55.0540541 913 4034 3209.0541 824.9459459 914 3160 3415.2723 -255.2722772 915 682 590.0833 91.9166667 916 3045 3209.0541 -164.0540541 917 3409 3632.9409 -223.9409091 918 3352 3415.2723 -63.2722772 919 3665 3415.2723 249.7277228 920 3020 3415.2723 -395.2722772 921 2557 3415.2723 -858.2722772 922 4404 3632.9409 771.0590909 923 2840 3209.0541 -369.0540541 924 3500 3415.2723 84.7277228 925 4063 3632.9409 430.0590909 926 4574 3632.9409 941.0590909 927 4006 3415.2723 590.7277228 928 3636 3201.7179 434.2820513 929 3494 3415.2723 78.7277228 930 3210 3415.2723 -205.2722772 931 3409 3415.2723 -6.2722772 932 2585 2942.5263 -357.5263158 933 2555 3016.7551 -461.7551020 934 3523 3632.9409 -109.9409091 935 3155 2942.5263 212.4736842 936 3381 3201.7179 179.2820513 937 2983 3415.2723 -432.2722772 938 3892 3671.9500 220.0500000 939 3324 3415.2723 -91.2722772 940 3296 3632.9409 -336.9409091 941 2443 2942.5263 -499.5263158 942 3580 3016.7551 563.2448980 943 2784 3632.9409 -848.9409091 944 3835 3016.7551 818.2448980 945 3721 3632.9409 88.0590909 946 3210 3395.2909 -185.2909091 947 3892 3632.9409 259.0590909 948 2841 3209.0541 -368.0540541 949 3068 3209.0541 -141.0540541 950 2690 2633.0877 56.9122807 951 3025 3415.2723 -390.2722772 952 4063 3395.2909 667.7090909 953 3551 3395.2909 155.7090909 954 3625 3209.0541 415.9459459 955 3210 3016.7551 193.2448980 956 2955 3209.0541 -254.0540541 957 3230 3209.0541 20.9459459 958 3154 3209.0541 -55.0540541 959 3438 3632.9409 -194.9409091 960 4404 3415.2723 988.7277228 961 3615 3209.0541 405.9459459 962 4659 3632.9409 1026.0590909 963 2841 3415.2723 -574.2722772 964 4006 3632.9409 373.0590909 965 2869 3209.0541 -340.0540541 966 470 590.0833 -120.0833333 967 3296 3632.9409 -336.9409091 968 3728 3415.2723 312.7277228 969 3068 3395.2909 -327.2909091 970 2585 3395.2909 -810.2909091 971 3296 3209.0541 86.9459459 972 3438 3632.9409 -194.9409091 973 3068 3632.9409 -564.9409091 974 3438 3209.0541 228.9459459 975 3040 2633.0877 406.9122807 976 3523 3395.2909 127.7090909 977 2926 3415.2723 -489.2722772 978 2301 2942.5263 -641.5263158 979 3409 3395.2909 13.7090909 980 3324 3415.2723 -91.2722772 981 3580 3395.2909 184.7090909 982 3722 3415.2723 306.7277228 983 2869 3209.0541 -340.0540541 984 4119 3632.9409 486.0590909 985 3540 3415.2723 124.7277228 986 3352 3395.2909 -43.2909091 987 4347 3016.7551 1330.2448980 988 3097 3395.2909 -298.2909091 989 4091 3209.0541 881.9459459 990 3381 3415.2723 -34.2722772 991 2727 3632.9409 -905.9409091 992 3864 3632.9409 231.0590909 993 2869 3632.9409 -763.9409091 994 3523 2942.5263 580.4736842 995 2528 3395.2909 -867.2909091 996 3324 3209.0541 114.9459459 997 2983 3016.7551 -33.7551020 998 4148 3632.9409 515.0590909 999 2760 2633.0877 126.9122807 1000 4480 3415.2723 1064.7277228 1001 3015 3395.2909 -380.2909091 1002 2898 3415.2723 -517.2722772 1003 3720 3415.2723 304.7277228 1004 4176 3209.0541 966.9459459 1005 3608 3632.9409 -24.9409091 1006 2756 2633.0877 122.9122807 1007 2950 3016.7551 -66.7551020 1008 3210 3632.9409 -422.9409091 1009 2926 3632.9409 -706.9409091 1010 3154 3415.2723 -261.2722772 1011 3090 3395.2909 -305.2909091 1012 3324 3209.0541 114.9459459 1013 2955 3395.2909 -440.2909091 1014 4432 3415.2723 1016.7277228 1015 3864 3395.2909 468.7090909 1016 3068 2942.5263 125.4736842 1017 3665 3209.0541 455.9459459 1018 3352 3632.9409 -280.9409091 1019 2815 3209.0541 -394.0540541 1020 3630 3632.9409 -2.9409091 1021 3040 3395.2909 -355.2909091 1022 3011 3209.0541 -198.0540541 1023 3523 3395.2909 127.7090909 1024 3494 3395.2909 98.7090909 1025 3466 3016.7551 449.2448980 1026 3921 3632.9409 288.0590909 1027 3835 3632.9409 202.0590909 1028 3381 3415.2723 -34.2722772 1029 3494 3632.9409 -138.9409091 1030 4100 3209.0541 890.9459459 1031 3779 3395.2909 383.7090909 1032 3864 3632.9409 231.0590909 1033 3324 3395.2909 -71.2909091 1034 2479 2942.5263 -463.5263158 1035 3693 3395.2909 297.7090909 1036 4148 3632.9409 515.0590909 1037 3210 3632.9409 -422.9409091 1038 3458 3415.2723 42.7277228 1039 4347 3632.9409 714.0590909 1040 3296 3016.7551 279.2448980 1041 2720 2942.5263 -222.5263158 1042 3636 3632.9409 3.0590909 1043 3665 2633.0877 1031.9122807 1044 3551 3395.2909 155.7090909 1045 4830 3209.0541 1620.9459459 1046 3864 3632.9409 231.0590909 1047 3210 3415.2723 -205.2722772 1048 3807 3415.2723 391.7277228 1049 4460 3395.2909 1064.7090909 1050 3608 3632.9409 -24.9409091 1051 3502 3209.0541 292.9459459 1052 3239 3395.2909 -156.2909091 1053 3892 3632.9409 259.0590909 1054 2995 2942.5263 52.4736842 1055 3665 3395.2909 269.7090909 1056 3381 3201.7179 179.2820513 1057 3466 3671.9500 -205.9500000 1058 3892 3632.9409 259.0590909 1059 3438 3632.9409 -194.9409091 1060 2955 3415.2723 -460.2722772 1061 3190 3632.9409 -442.9409091 1062 3060 3395.2909 -335.2909091 1063 2699 3209.0541 -510.0540541 1064 4375 3632.9409 742.0590909 1065 3722 3632.9409 89.0590909 1066 3580 3395.2909 184.7090909 1067 3494 3209.0541 284.9459459 1068 3835 3632.9409 202.0590909 1069 3750 3415.2723 334.7277228 1070 3152 3395.2909 -243.2909091 1071 3523 3209.0541 313.9459459 1072 3636 3632.9409 3.0590909 1073 3523 3395.2909 127.7090909 1074 2727 3209.0541 -482.0540541 1075 3040 3395.2909 -355.2909091 1076 3182 3395.2909 -213.2909091 1077 3494 3415.2723 78.7277228 1078 2813 3209.0541 -396.0540541 1079 2727 3632.9409 -905.9409091 1080 3875 3209.0541 665.9459459 1081 3921 3632.9409 288.0590909 1082 2642 3395.2909 -753.2909091 1083 2410 3395.2909 -985.2909091 1084 997 2014.8696 -1017.8695652 1085 3409 3395.2909 13.7090909 1086 2557 3209.0541 -652.0540541 1087 3523 3632.9409 -109.9409091 1088 3921 3632.9409 288.0590909 1089 3486 3632.9409 -146.9409091 1090 3739 3395.2909 343.7090909 1091 3864 3395.2909 468.7090909 1092 3267 3395.2909 -128.2909091 1093 3523 3395.2909 127.7090909 1094 3011 2942.5263 68.4736842 1095 1875 2014.8696 -139.8695652 1096 2330 2633.0877 -303.0877193 1097 3351 3632.9409 -281.9409091 1098 4006 3632.9409 373.0590909 1099 2671 3395.2909 -724.2909091 1100 2784 3632.9409 -848.9409091 1101 2700 3209.0541 -509.0540541 1102 3324 3395.2909 -71.2909091 1103 3660 3632.9409 27.0590909 1104 3125 3632.9409 -507.9409091 1105 3466 3415.2723 50.7277228 1106 3197 2942.5263 254.4736842 1107 3296 3395.2909 -99.2909091 1108 3415 3209.0541 205.9459459 1109 3320 3415.2723 -95.2722772 1110 3315 3415.2723 -100.2722772 1111 3068 3632.9409 -564.9409091 1112 4318 3632.9409 685.0590909 1113 4034 3632.9409 401.0590909 1114 3608 3632.9409 -24.9409091 1115 3352 3632.9409 -280.9409091 > if (par2 != 'none') { + print(cbind(as.factor(x[,par1]),predict(m))) + myt <- table(as.factor(x[,par1]),predict(m)) + print(myt) + } > postscript(file="/var/wessaorg/rcomp/tmp/402ay1323864782.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > if(par2=='none') { + op <- par(mfrow=c(2,2)) + plot(density(result$Actuals),main='Kernel Density Plot of Actuals') + plot(density(result$Residuals),main='Kernel Density Plot of Residuals') + plot(result$Forecasts,result$Actuals,main='Actuals versus Predictions',xlab='Predictions',ylab='Actuals') + plot(density(result$Forecasts),main='Kernel Density Plot of Predictions') + par(op) + } > if(par2!='none') { + plot(myt,main='Confusion Matrix',xlab='Actual',ylab='Predicted') + } > dev.off() null device 1 > if (par2 == 'none') { + detcoef <- cor(result$Forecasts,result$Actuals) + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goodness of Fit',2,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Correlation',1,TRUE) + a<-table.element(a,round(detcoef,4)) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'R-squared',1,TRUE) + a<-table.element(a,round(detcoef*detcoef,4)) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'RMSE',1,TRUE) + a<-table.element(a,round(sqrt(mean((result$Residuals)^2)),4)) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/5tme91323864782.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Actuals, Predictions, and Residuals',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'#',header=TRUE) + a<-table.element(a,'Actuals',header=TRUE) + a<-table.element(a,'Forecasts',header=TRUE) + a<-table.element(a,'Residuals',header=TRUE) + a<-table.row.end(a) + for (i in 1:length(result$Actuals)) { + a<-table.row.start(a) + a<-table.element(a,i,header=TRUE) + a<-table.element(a,result$Actuals[i]) + a<-table.element(a,result$Forecasts[i]) + a<-table.element(a,result$Residuals[i]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/6fknu1323864782.tab") + } > if (par2 != 'none') { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Confusion Matrix (predicted in columns / actuals in rows)',par3+1,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'',1,TRUE) + for (i in 1:par3) { + a<-table.element(a,paste('C',i,sep=''),1,TRUE) + } + a<-table.row.end(a) + for (i in 1:par3) { + a<-table.row.start(a) + a<-table.element(a,paste('C',i,sep=''),1,TRUE) + for (j in 1:par3) { + a<-table.element(a,myt[i,j]) + } + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/7bn751323864782.tab") + } > > try(system("convert tmp/2yzn41323864782.ps tmp/2yzn41323864782.png",intern=TRUE)) character(0) > try(system("convert tmp/3lw491323864782.ps tmp/3lw491323864782.png",intern=TRUE)) character(0) > try(system("convert tmp/402ay1323864782.ps tmp/402ay1323864782.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 33.504 0.501 34.010