R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-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(1418 + ,210907 + ,56 + ,396 + ,81 + ,3 + ,79 + ,30 + ,115 + ,94 + ,112285 + ,24188 + ,144 + ,145 + ,869 + ,120982 + ,56 + ,297 + ,55 + ,4 + ,58 + ,28 + ,109 + ,103 + ,84786 + ,18273 + ,103 + ,101 + ,1530 + ,176508 + ,54 + ,559 + ,50 + ,12 + ,60 + ,38 + ,146 + ,93 + ,83123 + ,14130 + ,98 + ,98 + ,2172 + ,179321 + ,89 + ,967 + ,125 + ,2 + ,108 + ,30 + ,116 + ,103 + ,101193 + ,32287 + ,135 + ,132 + ,901 + ,123185 + ,40 + ,270 + ,40 + ,1 + ,49 + ,22 + ,68 + ,51 + ,38361 + ,8654 + ,61 + ,60 + ,463 + ,52746 + ,25 + ,143 + ,37 + ,3 + ,0 + ,26 + ,101 + ,70 + ,68504 + ,9245 + ,39 + ,38 + ,3201 + ,385534 + ,92 + ,1562 + ,63 + ,0 + ,121 + ,25 + ,96 + ,91 + ,119182 + ,33251 + ,150 + ,144 + ,371 + ,33170 + ,18 + ,109 + ,44 + ,0 + ,1 + ,18 + ,67 + ,22 + ,22807 + ,1271 + ,5 + ,5 + ,1192 + ,101645 + ,63 + ,371 + ,88 + ,0 + ,20 + ,11 + ,44 + ,38 + ,17140 + ,5279 + ,28 + ,28 + ,1583 + ,149061 + ,44 + ,656 + ,66 + ,5 + ,43 + ,26 + ,100 + ,93 + ,116174 + ,27101 + ,84 + ,84 + ,1439 + ,165446 + ,33 + ,511 + ,57 + ,0 + ,69 + ,25 + ,93 + ,60 + ,57635 + ,16373 + ,80 + ,79 + ,1764 + ,237213 + ,84 + ,655 + ,74 + ,0 + ,78 + ,38 + ,140 + ,123 + ,66198 + ,19716 + ,130 + ,127 + ,1495 + ,173326 + ,88 + ,465 + ,49 + ,7 + ,86 + ,44 + ,166 + ,148 + ,71701 + ,17753 + ,82 + ,78 + ,1373 + ,133131 + ,55 + ,525 + ,52 + ,7 + ,44 + ,30 + ,99 + ,90 + ,57793 + ,9028 + ,60 + ,60 + ,2187 + ,258873 + ,60 + ,885 + ,88 + ,3 + ,104 + ,40 + ,139 + ,124 + ,80444 + ,18653 + ,131 + ,131 + ,1491 + ,180083 + ,66 + ,497 + ,36 + ,9 + ,63 + ,34 + ,130 + ,70 + ,53855 + ,8828 + ,84 + ,84 + ,4041 + ,324799 + ,154 + ,1436 + ,108 + ,0 + ,158 + ,47 + ,181 + ,168 + ,97668 + ,29498 + ,140 + ,133 + ,1706 + ,230964 + ,53 + ,612 + ,43 + ,4 + ,102 + ,30 + ,116 + ,115 + ,133824 + ,27563 + ,151 + ,150 + ,2152 + ,236785 + ,119 + ,865 + ,75 + ,3 + ,77 + ,31 + ,116 + ,71 + ,101481 + ,18293 + ,91 + ,91 + ,1036 + ,135473 + ,41 + ,385 + ,32 + ,0 + ,82 + ,23 + ,88 + ,66 + ,99645 + ,22530 + ,138 + ,132 + ,1882 + ,202925 + ,61 + ,567 + ,44 + ,7 + ,115 + ,36 + ,139 + ,134 + ,114789 + ,15977 + ,150 + ,136 + ,1929 + ,215147 + ,58 + ,639 + ,85 + ,0 + ,101 + ,36 + ,135 + ,117 + ,99052 + ,35082 + ,124 + ,124 + ,2242 + ,344297 + ,75 + ,963 + ,86 + ,1 + ,80 + ,30 + ,108 + ,108 + ,67654 + ,16116 + ,119 + ,118 + ,1220 + ,153935 + ,33 + ,398 + ,56 + ,5 + ,50 + ,25 + ,89 + ,84 + ,65553 + ,15849 + ,73 + ,70 + ,1289 + ,132943 + ,40 + ,410 + ,50 + ,7 + ,83 + ,39 + ,156 + ,156 + ,97500 + ,16026 + ,110 + ,107 + ,2515 + ,174724 + ,92 + ,966 + ,135 + ,0 + ,123 + ,34 + ,129 + ,120 + ,69112 + ,26569 + ,123 + ,119 + ,2147 + ,174415 + ,100 + ,801 + ,63 + ,0 + ,73 + ,31 + ,118 + ,114 + ,82753 + ,24785 + ,90 + ,89 + ,2352 + ,225548 + ,112 + ,892 + ,81 + ,5 + ,81 + ,31 + ,118 + ,94 + ,85323 + ,17569 + ,116 + ,112 + ,1638 + ,223632 + ,73 + ,513 + ,52 + ,0 + ,105 + ,33 + ,125 + ,120 + ,72654 + ,23825 + ,113 + ,108 + ,1222 + ,124817 + ,40 + ,469 + ,44 + ,0 + ,47 + ,25 + ,95 + ,81 + ,30727 + ,7869 + ,56 + ,52 + ,1812 + ,221698 + ,45 + ,683 + ,113 + ,0 + ,105 + ,33 + ,126 + ,110 + ,77873 + ,14975 + ,115 + ,112 + ,1677 + ,210767 + ,60 + ,643 + ,39 + ,3 + ,94 + ,35 + ,135 + ,133 + ,117478 + ,37791 + ,119 + ,116 + ,1579 + ,170266 + ,62 + ,535 + ,73 + ,4 + ,44 + ,42 + ,154 + ,122 + ,74007 + ,9605 + ,129 + ,123 + ,1731 + ,260561 + ,75 + ,625 + ,48 + ,1 + ,114 + ,43 + ,165 + ,158 + ,90183 + ,27295 + ,127 + ,125 + ,807 + ,84853 + ,31 + ,264 + ,33 + ,4 + ,38 + ,30 + ,113 + ,109 + ,61542 + ,2746 + ,27 + ,27 + ,2452 + ,294424 + ,77 + ,992 + ,59 + ,2 + ,107 + ,33 + ,127 + ,124 + ,101494 + ,34461 + ,175 + ,162 + ,829 + ,101011 + ,34 + ,238 + ,41 + ,0 + ,30 + ,13 + ,52 + ,39 + ,27570 + ,8098 + ,35 + ,32 + ,1940 + ,215641 + ,46 + ,818 + ,69 + ,0 + ,71 + ,32 + ,121 + ,92 + ,55813 + ,4787 + ,64 + ,64 + ,2662 + ,325107 + ,99 + ,937 + ,64 + ,0 + ,84 + ,36 + ,136 + ,126 + ,79215 + ,24919 + ,96 + ,92 + ,186 + ,7176 + ,17 + ,70 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1423 + ,603 + ,0 + ,0 + ,1499 + ,167542 + ,66 + ,507 + ,59 + ,2 + ,59 + ,28 + ,108 + ,70 + ,55461 + ,16329 + ,84 + ,83 + ,865 + ,106408 + ,30 + ,260 + ,32 + ,1 + ,33 + ,14 + ,46 + ,37 + ,31081 + ,12558 + ,41 + ,41 + ,1793 + ,96560 + ,76 + ,503 + ,129 + ,0 + ,42 + ,17 + ,54 + ,38 + ,22996 + ,7784 + ,47 + ,47 + ,2527 + ,265769 + ,146 + ,927 + ,37 + ,2 + ,96 + ,32 + ,124 + ,120 + ,83122 + ,28522 + ,126 + ,120 + ,2747 + ,269651 + ,67 + ,1269 + ,31 + ,10 + ,106 + ,30 + ,115 + ,93 + ,70106 + ,22265 + ,105 + ,105 + ,1324 + ,149112 + ,56 + ,537 + ,65 + ,6 + ,56 + ,35 + ,128 + ,95 + ,60578 + ,14459 + ,80 + ,79 + ,2702 + ,175824 + ,107 + ,910 + ,107 + ,0 + ,57 + ,20 + ,80 + ,77 + ,39992 + ,14526 + ,70 + ,65 + ,1383 + ,152871 + ,58 + ,532 + ,74 + ,5 + ,59 + ,28 + ,97 + ,90 + ,79892 + ,22240 + ,73 + ,70 + ,1179 + ,111665 + ,34 + ,345 + ,54 + ,4 + ,39 + ,28 + ,104 + ,80 + ,49810 + ,11802 + ,57 + ,55 + ,2099 + ,116408 + ,61 + ,918 + ,76 + ,1 + ,34 + ,39 + ,59 + ,31 + ,71570 + ,7623 + ,40 + ,39 + ,4308 + ,362301 + ,119 + ,1635 + ,715 + ,2 + ,76 + ,34 + ,125 + ,110 + ,100708 + ,11912 + ,68 + ,67 + ,918 + ,78800 + ,42 + ,330 + ,57 + ,2 + ,20 + ,26 + ,82 + ,66 + ,33032 + ,7935 + ,21 + ,21 + ,1831 + ,183167 + ,66 + ,557 + ,66 + ,0 + ,91 + ,39 + ,149 + ,138 + ,82875 + ,18220 + ,127 + ,127 + ,3373 + ,277965 + ,89 + ,1178 + ,106 + ,8 + ,115 + ,39 + ,149 + ,133 + ,139077 + ,19199 + ,154 + ,152 + ,1713 + ,150629 + ,44 + ,740 + ,54 + ,3 + ,85 + ,33 + ,122 + ,113 + ,71595 + ,19918 + ,116 + ,113 + ,1438 + ,168809 + ,66 + ,452 + ,32 + ,0 + ,76 + ,28 + ,118 + ,100 + ,72260 + ,21884 + ,102 + ,99 + ,496 + ,24188 + ,24 + ,218 + ,20 + ,0 + ,8 + ,4 + ,12 + ,7 + ,5950 + ,2694 + ,7 + ,7 + ,2253 + ,329267 + ,259 + ,764 + ,71 + ,8 + ,79 + ,39 + ,144 + ,140 + ,115762 + ,15808 + ,148 + ,141 + ,744 + ,65029 + ,17 + ,255 + ,21 + ,5 + ,21 + ,18 + ,67 + ,61 + ,32551 + ,3597 + ,21 + ,21 + ,1161 + ,101097 + ,64 + ,454 + ,70 + ,3 + ,30 + ,14 + ,52 + ,41 + ,31701 + ,5296 + ,35 + ,35 + ,2352 + ,218946 + ,41 + ,866 + ,112 + ,1 + ,76 + ,29 + ,108 + ,96 + ,80670 + ,25239 + ,112 + ,109 + ,2144 + ,244052 + ,68 + ,574 + ,66 + ,5 + ,101 + ,44 + ,166 + ,164 + ,143558 + ,29801 + ,137 + ,133 + ,4691 + ,341570 + ,168 + ,1276 + ,190 + ,1 + ,94 + ,21 + ,80 + ,78 + ,117105 + ,18450 + ,135 + ,123 + ,1112 + ,103597 + ,43 + ,379 + ,66 + ,1 + ,27 + ,16 + ,60 + ,49 + ,23789 + ,7132 + ,26 + ,26 + ,2694 + ,233328 + ,132 + ,825 + ,165 + ,5 + ,92 + ,28 + ,107 + ,102 + ,120733 + ,34861 + ,230 + ,230 + ,1973 + ,256462 + ,105 + ,798 + ,56 + ,0 + ,123 + ,35 + ,127 + ,124 + ,105195 + ,35940 + ,181 + ,166 + ,1769 + ,206161 + ,71 + ,663 + ,61 + ,12 + ,75 + ,28 + ,107 + ,99 + ,73107 + ,16688 + ,71 + ,68 + ,3148 + ,311473 + ,112 + ,1069 + ,53 + ,8 + ,128 + ,38 + ,146 + ,129 + ,132068 + ,24683 + ,147 + ,147 + ,2474 + ,235800 + ,94 + ,921 + ,127 + ,8 + ,105 + ,23 + ,84 + ,62 + ,149193 + ,46230 + ,190 + ,179 + ,2084 + ,177939 + ,82 + ,858 + ,63 + ,8 + ,55 + ,36 + ,141 + ,73 + ,46821 + ,10387 + ,64 + ,61 + ,1954 + ,207176 + ,70 + ,711 + ,38 + ,8 + ,56 + ,32 + ,123 + ,114 + ,87011 + ,21436 + ,105 + ,101 + ,1226 + ,196553 + ,57 + ,503 + ,50 + ,2 + ,41 + ,29 + ,111 + ,99 + ,95260 + ,30546 + ,107 + ,108 + ,1389 + ,174184 + ,53 + ,382 + ,52 + ,0 + ,72 + ,25 + ,98 + ,70 + ,55183 + ,19746 + ,94 + ,90 + ,1496 + ,143246 + ,103 + ,464 + ,42 + ,5 + ,67 + ,27 + ,105 + ,104 + ,106671 + ,15977 + ,116 + ,114 + ,2269 + ,187559 + ,121 + ,717 + ,76 + ,8 + ,75 + ,36 + ,135 + ,116 + ,73511 + ,22583 + ,106 + ,103 + ,1833 + ,187681 + ,62 + ,690 + ,67 + ,2 + ,114 + ,28 + ,107 + ,91 + ,92945 + ,17274 + ,143 + ,142 + ,1268 + ,119016 + ,52 + ,462 + ,50 + ,5 + ,118 + ,23 + ,85 + ,74 + ,78664 + ,16469 + ,81 + ,79 + ,1943 + ,182192 + ,52 + ,657 + ,53 + ,12 + ,77 + ,40 + ,155 + ,138 + ,70054 + ,14251 + ,89 + ,88 + ,893 + ,73566 + ,32 + ,385 + ,39 + ,6 + ,22 + ,23 + ,88 + ,67 + ,22618 + ,3007 + ,26 + ,25 + ,1762 + ,194979 + ,62 + ,577 + ,50 + ,7 + ,66 + ,40 + ,155 + ,151 + ,74011 + ,16851 + ,84 + ,83 + ,1403 + ,167488 + ,45 + ,619 + ,77 + ,2 + ,69 + ,28 + ,104 + ,72 + ,83737 + ,21113 + ,113 + ,113 + ,1425 + ,143756 + ,46 + ,479 + ,57 + ,0 + ,105 + ,34 + ,132 + ,120 + ,69094 + ,17401 + ,120 + ,118 + ,1857 + ,275541 + ,63 + ,817 + ,73 + ,4 + ,116 + ,33 + ,127 + ,115 + ,93133 + ,23958 + ,110 + ,110 + ,1840 + ,243199 + ,75 + ,752 + ,34 + ,3 + ,88 + ,28 + ,108 + ,105 + ,95536 + ,23567 + ,134 + ,129 + ,1502 + ,182999 + ,88 + ,430 + ,39 + ,6 + ,73 + ,34 + ,129 + ,104 + ,225920 + ,13065 + ,54 + ,51 + ,1441 + ,135649 + ,46 + ,451 + ,46 + ,2 + ,99 + ,30 + ,116 + ,108 + ,62133 + ,15358 + ,96 + ,93 + ,1420 + ,152299 + ,53 + ,537 + ,63 + ,0 + ,62 + ,33 + ,122 + ,98 + ,61370 + ,14587 + ,78 + ,76 + ,1416 + ,120221 + ,37 + ,519 + ,35 + ,1 + ,53 + ,22 + ,85 + ,69 + ,43836 + ,12770 + ,51 + ,49 + ,2970 + ,346485 + ,90 + ,1000 + ,106 + ,0 + ,118 + ,38 + ,147 + ,111 + ,106117 + ,24021 + ,121 + ,118 + ,1317 + ,145790 + ,63 + ,637 + ,43 + ,5 + ,30 + ,26 + ,99 + ,99 + ,38692 + ,9648 + ,38 + ,38 + ,1644 + ,193339 + ,78 + ,465 + ,47 + ,2 + ,100 + ,35 + ,87 + ,71 + ,84651 + ,20537 + ,145 + ,141 + ,870 + ,80953 + ,25 + ,437 + ,31 + ,0 + ,49 + ,8 + ,28 + ,27 + ,56622 + ,7905 + ,59 + ,58 + ,1654 + ,122774 + ,45 + ,711 + ,162 + ,0 + ,24 + ,24 + ,90 + ,69 + ,15986 + ,4527 + ,27 + ,27 + ,1054 + ,130585 + ,46 + ,299 + ,57 + ,5 + ,67 + ,29 + ,109 + ,107 + ,95364 + ,30495 + ,91 + ,91 + ,937 + ,112611 + ,41 + ,248 + ,36 + ,0 + ,46 + ,20 + ,78 + ,73 + ,26706 + ,7117 + ,48 + ,48 + ,3004 + ,286468 + ,144 + ,1162 + ,263 + ,1 + ,57 + ,29 + ,111 + ,107 + ,89691 + ,17719 + ,68 + ,63 + ,2008 + ,241066 + ,82 + ,714 + ,78 + ,0 + ,75 + ,45 + ,158 + ,93 + ,67267 + ,27056 + ,58 + ,56 + ,2547 + ,148446 + ,91 + ,905 + ,63 + ,1 + ,135 + ,37 + ,141 + ,129 + ,126846 + ,33473 + ,150 + ,144 + ,1885 + ,204713 + ,71 + ,649 + ,54 + ,1 + ,68 + ,33 + ,122 + ,69 + ,41140 + ,9758 + ,74 + ,73 + ,1626 + ,182079 + ,63 + ,512 + ,63 + ,2 + ,124 + ,33 + ,124 + ,118 + ,102860 + ,21115 + ,181 + ,168 + ,1468 + ,140344 + ,53 + ,472 + ,77 + ,6 + ,33 + ,25 + ,93 + ,73 + ,51715 + ,7236 + ,65 + ,64 + ,2445 + ,220516 + ,62 + ,905 + ,79 + ,1 + ,98 + ,32 + ,124 + ,119 + ,55801 + ,13790 + ,97 + ,97 + ,1964 + ,243060 + ,63 + ,786 + ,110 + ,4 + ,58 + ,29 + ,112 + ,104 + ,111813 + ,32902 + ,121 + ,117 + ,1381 + ,162765 + ,32 + ,489 + ,56 + ,2 + ,68 + ,28 + ,108 + ,107 + ,120293 + ,25131 + ,99 + ,100 + ,1369 + ,182613 + ,39 + ,479 + ,56 + ,3 + ,81 + ,28 + ,99 + ,99 + ,138599 + ,30910 + ,152 + ,149 + ,1659 + ,232138 + ,62 + ,617 + ,43 + ,0 + ,131 + ,31 + ,117 + ,90 + ,161647 + ,35947 + ,188 + ,187 + ,2888 + ,265318 + ,117 + ,925 + ,111 + ,10 + ,110 + ,52 + ,199 + ,197 + ,115929 + ,29848 + ,138 + ,127 + ,1290 + ,85574 + ,34 + ,351 + ,71 + ,0 + ,37 + ,21 + ,78 + ,36 + ,24266 + ,6943 + ,40 + ,37 + ,2845 + ,310839 + ,92 + ,1144 + ,62 + ,9 + ,130 + ,24 + ,91 + ,85 + ,162901 + ,42705 + ,254 + ,245 + ,1982 + ,225060 + ,93 + ,669 + ,56 + ,7 + ,93 + ,41 + ,158 + ,139 + ,109825 + ,31808 + ,87 + ,87 + ,1904 + ,232317 + ,54 + ,707 + ,74 + ,0 + ,118 + ,33 + ,126 + ,106 + ,129838 + ,26675 + ,178 + ,177 + ,1391 + ,144966 + ,144 + ,458 + ,60 + ,0 + ,39 + ,32 + ,122 + ,50 + ,37510 + ,8435 + ,51 + ,49 + ,602 + ,43287 + ,14 + ,214 + ,43 + ,4 + ,13 + ,19 + ,71 + ,64 + ,43750 + ,7409 + ,49 + ,49 + ,1743 + ,155754 + ,61 + ,599 + ,68 + ,4 + ,74 + ,20 + ,75 + ,31 + ,40652 + ,14993 + ,73 + ,73 + ,1559 + ,164709 + ,109 + ,572 + ,53 + ,0 + ,81 + ,31 + ,115 + ,63 + ,87771 + ,36867 + ,176 + ,177 + ,2014 + ,201940 + ,38 + ,897 + ,87 + ,0 + ,109 + ,31 + ,119 + ,92 + ,85872 + ,33835 + ,94 + ,94 + ,2143 + ,235454 + ,73 + ,819 + ,46 + ,0 + ,151 + ,32 + ,124 + ,106 + ,89275 + ,24164 + ,120 + ,117 + ,2146 + ,220801 + ,75 + ,720 + ,105 + ,1 + ,51 + ,18 + ,72 + ,63 + ,44418 + ,12607 + ,66 + ,60 + ,874 + ,99466 + ,50 + ,273 + ,32 + ,0 + ,28 + ,23 + ,91 + ,69 + ,192565 + ,22609 + ,56 + ,55 + ,1590 + ,92661 + ,61 + ,508 + ,133 + ,1 + ,40 + ,17 + ,45 + ,41 + ,35232 + ,5892 + ,39 + ,39 + ,1590 + ,133328 + ,55 + ,506 + ,79 + ,0 + ,56 + ,20 + ,78 + ,56 + ,40909 + ,17014 + ,66 + ,64 + ,1210 + ,61361 + ,77 + ,451 + ,51 + ,0 + ,27 + ,12 + ,39 + ,25 + ,13294 + ,5394 + ,27 + ,26 + ,2072 + ,125930 + ,75 + ,699 + ,207 + ,4 + ,37 + ,17 + ,68 + ,65 + ,32387 + ,9178 + ,65 + ,64 + ,1281 + ,100750 + ,72 + ,407 + ,67 + ,0 + ,83 + ,30 + ,119 + ,93 + ,140867 + ,6440 + ,58 + ,58 + ,1401 + ,224549 + ,50 + ,465 + ,47 + ,4 + ,54 + ,31 + ,117 + ,114 + ,120662 + ,21916 + ,98 + ,95 + ,834 + ,82316 + ,32 + ,245 + ,34 + ,4 + ,27 + ,10 + ,39 + ,38 + ,21233 + ,4011 + ,25 + ,25 + ,1105 + ,102010 + ,53 + ,370 + ,66 + ,3 + ,28 + ,13 + ,50 + ,44 + ,44332 + ,5818 + ,26 + ,26 + ,1272 + ,101523 + ,42 + ,316 + ,76 + ,0 + ,59 + ,22 + ,88 + ,87 + ,61056 + ,18647 + ,77 + ,76 + ,1944 + ,243511 + ,71 + ,603 + ,65 + ,0 + ,133 + ,42 + ,155 + ,110 + ,101338 + ,20556 + ,130 + ,129 + ,391 + ,22938 + ,10 + ,154 + ,9 + ,0 + ,12 + ,1 + ,0 + ,0 + ,1168 + ,238 + ,11 + ,11 + ,761 + ,41566 + ,35 + ,229 + ,42 + ,5 + ,0 + ,9 + ,36 + ,27 + ,13497 + ,70 + ,2 + ,2 + ,1605 + ,152474 + ,65 + ,577 + ,45 + ,0 + ,106 + ,32 + ,123 + ,83 + ,65567 + ,22392 + ,101 + ,101 + ,530 + ,61857 + ,25 + ,192 + ,25 + ,4 + ,23 + ,11 + ,32 + ,30 + ,25162 + ,3913 + ,31 + ,28 + ,1988 + ,99923 + ,66 + ,617 + ,115 + ,0 + ,44 + ,25 + ,99 + ,80 + ,32334 + ,12237 + ,36 + ,36 + ,1386 + ,132487 + ,41 + ,411 + ,97 + ,0 + ,71 + ,36 + ,136 + ,98 + ,40735 + ,8388 + ,120 + ,89 + ,2395 + ,317394 + ,86 + ,975 + ,53 + ,1 + ,116 + ,31 + ,117 + ,82 + ,91413 + ,22120 + ,195 + ,193 + ,387 + ,21054 + ,16 + ,146 + ,2 + ,0 + ,4 + ,0 + ,0 + ,0 + ,855 + ,338 + ,4 + ,4 + ,1742 + ,209641 + ,42 + ,705 + ,52 + ,5 + ,62 + ,24 + ,88 + ,60 + ,97068 + ,11727 + ,89 + ,84 + ,620 + ,22648 + ,19 + ,184 + ,44 + ,0 + ,12 + ,13 + ,39 + ,28 + ,44339 + ,3704 + ,24 + ,23 + ,449 + ,31414 + ,19 + ,200 + ,22 + ,0 + ,18 + ,8 + ,25 + ,9 + ,14116 + ,3988 + ,39 + ,39 + ,800 + ,46698 + ,45 + ,274 + ,35 + ,0 + ,14 + ,13 + ,52 + ,33 + ,10288 + ,3030 + ,14 + ,14 + ,1684 + ,131698 + ,65 + ,502 + ,74 + ,0 + ,60 + ,19 + ,75 + ,59 + ,65622 + ,13520 + ,78 + ,78 + ,1050 + ,91735 + ,35 + ,382 + ,103 + ,0 + ,7 + ,18 + ,71 + ,49 + ,16563 + ,1421 + ,15 + ,14 + ,2699 + ,244749 + ,95 + ,964 + ,144 + ,2 + ,98 + ,33 + ,124 + ,115 + ,76643 + ,20923 + ,106 + ,101 + ,1606 + ,184510 + ,49 + ,537 + ,60 + ,7 + ,64 + ,40 + ,151 + ,140 + ,110681 + ,20237 + ,83 + ,82 + ,1502 + ,79863 + ,37 + ,438 + ,134 + ,1 + ,29 + ,22 + ,71 + ,49 + ,29011 + ,3219 + ,24 + ,24 + ,1204 + ,128423 + ,64 + ,369 + ,89 + ,8 + ,32 + ,38 + ,145 + ,120 + ,92696 + ,3769 + ,37 + ,36 + ,1138 + ,97839 + ,38 + ,417 + ,42 + ,2 + ,25 + ,24 + ,87 + ,66 + ,94785 + ,12252 + ,77 + ,75 + ,568 + ,38214 + ,34 + ,276 + ,52 + ,0 + ,16 + ,8 + ,27 + ,21 + ,8773 + ,1888 + ,16 + ,16 + ,1459 + ,151101 + ,32 + ,514 + ,98 + ,2 + ,48 + ,35 + ,131 + ,124 + ,83209 + ,14497 + ,56 + ,55 + ,2158 + ,272458 + ,65 + ,822 + ,99 + ,0 + ,100 + ,43 + ,162 + ,152 + ,93815 + ,28864 + ,132 + ,131 + ,1111 + ,172494 + ,52 + ,389 + ,52 + ,0 + ,46 + ,43 + ,165 + ,139 + ,86687 + ,21721 + ,144 + ,131 + ,1421 + ,108043 + ,62 + ,466 + ,29 + ,1 + ,45 + ,14 + ,54 + ,38 + ,34553 + ,4821 + ,40 + ,39 + ,2833 + ,328107 + ,65 + ,1255 + ,125 + ,3 + ,129 + ,41 + ,159 + ,144 + ,105547 + ,33644 + ,153 + ,144 + ,1955 + ,250579 + ,83 + ,694 + ,106 + ,0 + ,130 + ,38 + ,147 + ,120 + ,103487 + ,15923 + ,143 + ,139 + ,2922 + ,351067 + ,95 + ,1024 + ,95 + ,3 + ,136 + ,45 + ,170 + ,160 + ,213688 + ,42935 + ,220 + ,211 + ,1002 + ,158015 + ,29 + ,400 + ,40 + ,0 + ,59 + ,31 + ,119 + ,114 + ,71220 + ,18864 + ,79 + ,78 + ,1060 + ,98866 + ,18 + ,397 + ,140 + ,0 + ,25 + ,13 + ,49 + ,39 + ,23517 + ,4977 + ,50 + ,50 + ,956 + ,85439 + ,33 + ,350 + ,43 + ,0 + ,32 + ,28 + ,104 + ,78 + ,56926 + ,7785 + ,39 + ,39 + ,2186 + ,229242 + ,247 + ,719 + ,128 + ,4 + ,63 + ,31 + ,120 + ,119 + ,91721 + ,17939 + ,95 + ,90 + ,3604 + ,351619 + ,139 + ,1277 + ,142 + ,4 + ,95 + ,40 + ,150 + ,141 + ,115168 + ,23436 + ,169 + ,166 + ,1035 + ,84207 + ,29 + ,356 + ,73 + ,11 + ,14 + ,30 + ,112 + ,101 + ,111194 + ,325 + ,12 + ,12 + ,1417 + ,120445 + ,118 + ,457 + ,72 + ,0 + ,36 + ,16 + ,59 + ,56 + ,51009 + ,13539 + ,63 + ,57 + ,3261 + ,324598 + ,110 + ,1402 + ,128 + ,0 + ,113 + ,37 + ,136 + ,133 + ,135777 + ,34538 + ,134 + ,133 + ,1587 + ,131069 + ,67 + ,600 + ,61 + ,4 + ,47 + ,30 + ,107 + ,83 + ,51513 + ,12198 + ,69 + ,69 + ,1424 + ,204271 + ,42 + ,480 + ,73 + ,0 + ,92 + ,35 + ,130 + ,116 + ,74163 + ,26924 + ,119 + ,119 + ,1701 + ,165543 + ,65 + ,595 + ,148 + ,1 + ,70 + ,32 + ,115 + ,90 + ,51633 + ,12716 + ,119 + ,119 + ,1249 + ,141722 + ,94 + ,436 + ,64 + ,0 + ,19 + ,27 + ,107 + ,36 + ,75345 + ,8172 + ,75 + ,65 + ,946 + ,116048 + ,64 + ,230 + ,45 + ,0 + ,50 + ,20 + ,75 + ,50 + ,33416 + ,10855 + ,63 + ,61 + ,1926 + ,250047 + ,81 + ,651 + ,58 + ,0 + ,41 + ,18 + ,71 + ,61 + ,83305 + ,11932 + ,55 + ,49 + ,3352 + ,299775 + ,95 + ,1367 + ,97 + ,9 + ,91 + ,31 + ,120 + ,97 + ,98952 + ,14300 + ,103 + ,101 + ,1641 + ,195838 + ,67 + ,564 + ,50 + ,1 + ,111 + ,31 + ,116 + ,98 + ,102372 + ,25515 + ,197 + ,196 + ,2035 + ,173260 + ,63 + ,716 + ,37 + ,3 + ,41 + ,21 + ,79 + ,78 + ,37238 + ,2805 + ,16 + ,15 + ,2312 + ,254488 + ,83 + ,747 + ,50 + ,10 + ,120 + ,39 + ,150 + ,117 + ,103772 + ,29402 + ,140 + ,136 + ,1369 + ,104389 + ,45 + ,467 + ,105 + ,5 + ,135 + ,41 + ,156 + ,148 + ,123969 + ,16440 + ,89 + ,89 + ,1577 + ,136084 + ,30 + ,671 + ,69 + ,0 + ,27 + ,13 + ,51 + ,41 + ,27142 + ,11221 + ,40 + ,40 + ,2201 + ,199476 + ,70 + ,861 + ,46 + ,2 + ,87 + ,32 + ,118 + ,105 + ,135400 + ,28732 + ,125 + ,123 + ,961 + ,92499 + ,32 + ,319 + ,57 + ,0 + ,25 + ,18 + ,71 + ,55 + ,21399 + ,5250 + ,21 + ,21 + ,1900 + ,224330 + ,83 + ,612 + ,52 + ,1 + ,131 + ,39 + ,144 + ,132 + ,130115 + ,28608 + ,167 + ,163 + ,1254 + ,135781 + ,31 + ,433 + ,98 + ,2 + ,45 + ,14 + ,47 + ,44 + ,24874 + ,8092 + ,32 + ,29 + ,1335 + ,74408 + ,67 + ,434 + ,61 + ,4 + ,29 + ,7 + ,28 + ,21 + ,34988 + ,4473 + ,36 + ,35 + ,1597 + ,81240 + ,66 + ,503 + ,89 + ,0 + ,58 + ,17 + ,68 + ,50 + ,45549 + ,1572 + ,13 + ,13 + ,207 + ,14688 + ,10 + ,85 + ,0 + ,0 + ,4 + ,0 + ,0 + ,0 + ,6023 + ,2065 + ,5 + ,5 + ,1645 + ,181633 + ,70 + ,564 + ,48 + ,2 + ,47 + ,30 + ,110 + ,73 + ,64466 + ,14817 + ,96 + ,96 + ,2429 + ,271856 + ,103 + ,824 + ,91 + ,1 + ,109 + ,37 + ,147 + ,86 + ,54990 + ,16714 + ,151 + ,151 + ,151 + ,7199 + ,5 + ,74 + ,0 + ,0 + ,7 + ,0 + ,0 + ,0 + ,1644 + ,556 + ,6 + ,6 + ,474 + ,46660 + ,20 + ,259 + ,7 + ,0 + ,12 + ,5 + ,15 + ,13 + ,6179 + ,2089 + ,13 + ,13 + ,141 + ,17547 + ,5 + ,69 + ,3 + ,0 + ,0 + ,1 + ,4 + ,4 + ,3926 + ,2658 + ,3 + ,3 + ,1639 + ,133368 + ,36 + ,535 + ,54 + ,1 + ,37 + ,16 + ,64 + ,57 + ,32755 + ,10695 + ,57 + ,56 + ,872 + ,95227 + ,34 + ,239 + ,70 + ,0 + ,37 + ,32 + ,111 + ,48 + ,34777 + ,1669 + ,23 + ,23 + ,1318 + ,152601 + ,48 + ,438 + ,36 + ,2 + ,46 + ,24 + ,85 + ,46 + ,73224 + ,16267 + ,61 + ,57 + ,1018 + ,98146 + ,40 + ,459 + ,37 + ,0 + ,15 + ,17 + ,68 + ,48 + ,27114 + ,7768 + ,21 + ,14 + ,1383 + ,79619 + ,43 + ,426 + ,123 + ,3 + ,42 + ,11 + ,40 + ,32 + ,20760 + ,7252 + ,43 + ,43 + ,1314 + ,59194 + ,31 + ,288 + ,247 + ,6 + ,7 + ,24 + ,80 + ,68 + ,37636 + ,6387 + ,20 + ,20 + ,1335 + ,139942 + ,42 + ,498 + ,46 + ,0 + ,54 + ,22 + ,88 + ,87 + ,65461 + ,18715 + ,82 + ,72 + ,1403 + ,118612 + ,46 + ,454 + ,72 + ,2 + ,54 + ,12 + ,48 + ,43 + ,30080 + ,7936 + ,90 + ,87 + ,910 + ,72880 + ,33 + ,376 + ,41 + ,0 + ,14 + ,19 + ,76 + ,67 + ,24094 + ,8643 + ,25 + ,21 + ,616 + ,65475 + ,18 + ,225 + ,24 + ,2 + ,16 + ,13 + ,51 + ,46 + ,69008 + ,7294 + ,60 + ,56 + ,1407 + ,99643 + ,55 + ,555 + ,45 + ,1 + ,33 + ,17 + ,67 + ,46 + ,54968 + ,4570 + ,61 + ,59 + ,771 + ,71965 + ,35 + ,252 + ,33 + ,1 + ,32 + ,15 + ,59 + ,56 + ,46090 + ,7185 + ,85 + ,82 + ,766 + ,77272 + ,59 + ,208 + ,27 + ,2 + ,21 + ,16 + ,61 + ,48 + ,27507 + ,10058 + ,43 + ,43 + ,473 + ,49289 + ,19 + ,130 + ,36 + ,1 + ,15 + ,24 + ,76 + ,44 + ,10672 + ,2342 + ,25 + ,25 + ,1376 + ,135131 + ,66 + ,481 + ,87 + ,0 + ,38 + ,15 + ,60 + ,60 + ,34029 + ,8509 + ,41 + ,38 + ,1232 + ,108446 + ,60 + ,389 + ,90 + ,1 + ,22 + ,17 + ,68 + ,65 + ,46300 + ,13275 + ,26 + ,25 + ,1521 + ,89746 + ,36 + ,565 + ,114 + ,3 + ,28 + ,18 + ,71 + ,55 + ,24760 + ,6816 + ,38 + ,38 + ,572 + ,44296 + ,25 + ,173 + ,31 + ,0 + ,10 + ,20 + ,76 + ,38 + ,18779 + ,1930 + ,12 + ,12 + ,1059 + ,77648 + ,47 + ,278 + ,45 + ,0 + ,31 + ,16 + ,62 + ,52 + ,21280 + ,8086 + ,29 + ,29 + ,1544 + ,181528 + ,54 + ,609 + ,69 + ,0 + ,32 + ,16 + ,61 + ,60 + ,40662 + ,10737 + ,49 + ,47 + ,1230 + ,134019 + ,53 + ,422 + ,51 + ,0 + ,32 + ,18 + ,67 + ,54 + ,28987 + ,8033 + ,46 + ,45 + ,1206 + ,124064 + ,40 + ,445 + ,34 + ,1 + ,43 + ,22 + ,88 + ,86 + ,22827 + ,7058 + ,41 + ,40 + ,1205 + ,92630 + ,40 + ,387 + ,60 + ,4 + ,27 + ,8 + ,30 + ,24 + ,18513 + ,6782 + ,31 + ,30 + ,1255 + ,121848 + ,39 + ,339 + ,45 + ,0 + ,37 + ,17 + ,64 + ,52 + ,30594 + ,5401 + ,41 + ,41 + ,613 + ,52915 + ,14 + ,181 + ,54 + ,0 + ,20 + ,18 + ,68 + ,49 + ,24006 + ,6521 + ,26 + ,25 + ,721 + ,81872 + ,45 + ,245 + ,25 + ,0 + ,32 + ,16 + ,64 + ,61 + ,27913 + ,10856 + ,23 + ,23 + ,1109 + ,58981 + ,36 + ,384 + ,38 + ,7 + ,0 + ,23 + ,91 + ,61 + ,42744 + ,2154 + ,14 + ,14 + ,740 + ,53515 + ,28 + ,212 + ,52 + ,2 + ,5 + ,22 + ,88 + ,81 + ,12934 + ,6117 + ,16 + ,16 + ,1126 + ,60812 + ,44 + ,399 + ,67 + ,0 + ,26 + ,13 + ,52 + ,43 + ,22574 + ,5238 + ,25 + ,26 + ,728 + ,56375 + ,30 + ,229 + ,74 + ,7 + ,10 + ,13 + ,49 + ,40 + ,41385 + ,4820 + ,21 + ,21 + ,689 + ,65490 + ,22 + ,224 + ,38 + ,3 + ,27 + ,16 + ,62 + ,40 + ,18653 + ,5615 + ,32 + ,27 + ,592 + ,80949 + ,17 + ,203 + ,30 + ,0 + ,11 + ,16 + ,61 + ,56 + ,18472 + ,4272 + ,9 + ,9 + ,995 + ,76302 + ,31 + ,333 + ,26 + ,0 + ,29 + ,20 + ,76 + ,68 + ,30976 + ,8702 + ,35 + ,33 + ,1613 + ,104011 + ,55 + ,384 + ,67 + ,6 + ,25 + ,22 + ,88 + ,79 + ,63339 + ,15340 + ,42 + ,42 + ,2048 + ,98104 + ,54 + ,636 + ,132 + ,2 + ,55 + ,17 + ,66 + ,47 + ,25568 + ,8030 + ,68 + ,68 + ,705 + ,67989 + ,21 + ,185 + ,42 + ,0 + ,23 + ,18 + ,71 + ,57 + ,33747 + ,9526 + ,32 + ,32 + ,301 + ,30989 + ,14 + ,93 + ,35 + ,0 + ,5 + ,17 + ,68 + ,41 + ,4154 + ,1278 + ,6 + ,6 + ,1803 + ,135458 + ,81 + ,581 + ,118 + ,3 + ,43 + ,12 + ,48 + ,29 + ,19474 + ,4236 + ,68 + ,67 + ,799 + ,73504 + ,35 + ,248 + ,68 + ,0 + ,23 + ,7 + ,25 + ,3 + ,35130 + ,3023 + ,33 + ,33 + ,861 + ,63123 + ,43 + ,304 + ,43 + ,1 + ,34 + ,17 + ,68 + ,60 + ,39067 + ,7196 + ,84 + ,77 + ,1186 + ,61254 + ,46 + ,344 + ,76 + ,1 + ,36 + ,14 + ,41 + ,30 + ,13310 + ,3394 + ,46 + ,46 + ,1451 + ,74914 + ,30 + ,407 + ,64 + ,0 + ,35 + ,23 + ,90 + ,79 + ,65892 + ,6371 + ,30 + ,30 + ,628 + ,31774 + ,23 + ,170 + ,48 + ,1 + ,0 + ,17 + ,66 + ,47 + ,4143 + ,1574 + ,0 + ,0 + ,1161 + ,81437 + ,38 + ,312 + ,64 + ,0 + ,37 + ,14 + ,54 + ,40 + ,28579 + ,9620 + ,36 + ,36 + ,1463 + ,87186 + ,54 + ,507 + ,56 + ,0 + ,28 + ,15 + ,59 + ,48 + ,51776 + ,6978 + ,47 + ,46 + ,742 + ,50090 + ,20 + ,224 + ,71 + ,0 + ,16 + ,17 + ,60 + ,36 + ,21152 + ,4911 + ,20 + ,18 + ,979 + ,65745 + ,53 + ,340 + ,75 + ,0 + ,26 + ,21 + ,77 + ,42 + ,38084 + ,8645 + ,50 + ,48 + ,675 + ,56653 + ,45 + ,168 + ,39 + ,0 + ,38 + ,18 + ,68 + ,49 + ,27717 + ,8987 + ,30 + ,29 + ,1241 + ,158399 + ,39 + ,443 + ,42 + ,0 + ,23 + ,18 + ,72 + ,57 + ,32928 + ,5544 + ,30 + ,28 + ,676 + ,46455 + ,20 + ,204 + ,39 + ,0 + ,22 + ,17 + ,67 + ,12 + ,11342 + ,3083 + ,34 + ,34 + ,1049 + ,73624 + ,24 + ,367 + ,93 + ,0 + ,30 + ,17 + ,64 + ,40 + ,19499 + ,6909 + ,33 + ,33 + ,620 + ,38395 + ,31 + ,210 + ,38 + ,0 + ,16 + ,16 + ,63 + ,43 + ,16380 + ,3189 + ,34 + ,34 + ,1081 + ,91899 + ,35 + ,335 + ,60 + ,0 + ,18 + ,15 + ,59 + ,33 + ,36874 + ,6745 + ,37 + ,33 + ,1688 + ,139526 + ,151 + ,364 + ,71 + ,0 + ,28 + ,21 + ,84 + ,77 + ,48259 + ,16724 + ,83 + ,80 + ,736 + ,52164 + ,52 + ,178 + ,52 + ,0 + ,32 + ,16 + ,64 + ,43 + ,16734 + ,4850 + ,32 + ,32 + ,617 + ,51567 + ,30 + ,206 + ,27 + ,2 + ,21 + ,14 + ,56 + ,45 + ,28207 + ,7025 + ,30 + ,30 + ,812 + ,70551 + ,31 + ,279 + ,59 + ,0 + ,23 + ,15 + ,54 + ,47 + ,30143 + ,6047 + ,43 + ,41 + ,1051 + ,84856 + ,29 + ,387 + ,40 + ,1 + ,29 + ,17 + ,67 + ,43 + ,41369 + ,7377 + ,41 + ,41 + ,1656 + ,102538 + ,57 + ,490 + ,79 + ,1 + ,50 + ,15 + ,58 + ,45 + ,45833 + ,9078 + ,51 + ,51 + ,705 + ,86678 + ,40 + ,238 + ,44 + ,0 + ,12 + ,15 + ,59 + ,50 + ,29156 + ,4605 + ,19 + ,18 + ,945 + ,85709 + ,44 + ,343 + ,65 + ,0 + ,21 + ,10 + ,40 + ,35 + ,35944 + ,3238 + ,37 + ,34 + ,554 + ,34662 + ,25 + ,232 + ,10 + ,0 + ,18 + ,6 + ,22 + ,7 + ,36278 + ,8100 + ,33 + ,31 + ,1597 + ,150580 + ,77 + ,530 + ,124 + ,0 + ,27 + ,22 + ,83 + ,71 + ,45588 + ,9653 + ,41 + ,39 + ,982 + ,99611 + ,35 + ,291 + ,81 + ,0 + ,41 + ,21 + ,81 + ,67 + ,45097 + ,8914 + ,54 + ,54 + ,222 + ,19349 + ,11 + ,67 + ,15 + ,0 + ,13 + ,1 + ,2 + ,0 + ,3895 + ,786 + ,14 + ,14 + ,1212 + ,99373 + ,63 + ,397 + ,92 + ,1 + ,12 + ,18 + ,72 + ,62 + ,28394 + ,6700 + ,25 + ,24 + ,1143 + ,86230 + ,44 + ,467 + ,42 + ,0 + ,21 + ,17 + ,61 + ,54 + ,18632 + ,5788 + ,25 + ,24 + ,435 + ,30837 + ,19 + ,178 + ,10 + ,0 + ,8 + ,4 + ,15 + ,4 + ,2325 + ,593 + ,8 + ,8 + ,532 + ,31706 + ,13 + ,175 + ,24 + ,0 + ,26 + ,10 + ,32 + ,25 + ,25139 + ,4506 + ,26 + ,26 + ,882 + ,89806 + ,42 + ,299 + ,64 + ,0 + ,27 + ,16 + ,62 + ,40 + ,27975 + ,6382 + ,20 + ,19 + ,608 + ,62088 + ,38 + ,154 + ,45 + ,1 + ,13 + ,16 + ,58 + ,38 + ,14483 + ,5621 + ,11 + ,11 + ,459 + ,40151 + ,29 + ,106 + ,22 + ,0 + ,16 + ,9 + ,36 + ,19 + ,13127 + ,3997 + ,14 + ,14 + ,578 + ,27634 + ,20 + ,189 + ,56 + ,0 + ,2 + ,16 + ,59 + ,17 + ,5839 + ,520 + ,3 + ,1 + ,826 + ,76990 + ,27 + ,194 + ,94 + ,0 + ,42 + ,17 + ,68 + ,67 + ,24069 + ,8891 + ,40 + ,39 + ,509 + ,37460 + ,20 + ,135 + ,19 + ,0 + ,5 + ,7 + ,21 + ,14 + ,3738 + ,999 + ,5 + ,5 + ,717 + ,54157 + ,19 + ,201 + ,35 + ,0 + ,37 + ,15 + ,55 + ,30 + ,18625 + ,7067 + ,38 + ,37 + ,637 + ,49862 + ,37 + ,207 + ,32 + ,0 + ,17 + ,14 + ,54 + ,54 + ,36341 + ,4639 + ,32 + ,32 + ,857 + ,84337 + ,26 + ,280 + ,35 + ,0 + ,38 + ,14 + ,55 + ,35 + ,24548 + ,5654 + ,41 + ,38 + ,830 + ,64175 + ,42 + ,260 + ,48 + ,0 + ,37 + ,18 + ,72 + ,59 + ,21792 + ,6928 + ,46 + ,47 + ,652 + ,59382 + ,49 + ,227 + ,49 + ,0 + ,29 + ,12 + ,41 + ,24 + ,26263 + ,1514 + ,47 + ,47 + ,707 + ,119308 + ,30 + ,239 + ,48 + ,0 + ,32 + ,16 + ,61 + ,58 + ,23686 + ,9238 + ,37 + ,37 + ,954 + ,76702 + ,49 + ,333 + ,62 + ,0 + ,35 + ,21 + ,67 + ,42 + ,49303 + ,8204 + ,51 + ,51 + ,1461 + ,103425 + ,67 + ,428 + ,96 + ,1 + ,17 + ,19 + ,76 + ,46 + ,25659 + ,5926 + ,49 + ,45 + ,672 + ,70344 + ,28 + ,230 + ,45 + ,0 + ,20 + ,16 + ,64 + ,61 + ,28904 + ,5785 + ,21 + ,21 + ,778 + ,43410 + ,19 + ,292 + ,63 + ,0 + ,7 + ,1 + ,3 + ,3 + ,2781 + ,4 + ,1 + ,1 + ,1141 + ,104838 + ,49 + ,350 + ,71 + ,1 + ,46 + ,16 + ,63 + ,52 + ,29236 + ,5930 + ,44 + ,42 + ,680 + ,62215 + ,27 + ,186 + ,26 + ,0 + ,24 + ,10 + ,40 + ,25 + ,19546 + ,3710 + ,26 + ,26 + ,1090 + ,69304 + ,30 + ,326 + ,48 + ,6 + ,40 + ,19 + ,69 + ,40 + ,22818 + ,705 + ,21 + ,21 + ,616 + ,53117 + ,22 + ,155 + ,29 + ,3 + ,3 + ,12 + ,48 + ,32 + ,32689 + ,443 + ,4 + ,4 + ,285 + ,19764 + ,12 + ,75 + ,19 + ,1 + ,10 + ,2 + ,8 + ,4 + ,5752 + ,2416 + ,10 + ,10 + ,1145 + ,86680 + ,31 + ,361 + ,45 + ,2 + ,37 + ,14 + ,52 + ,49 + ,22197 + ,7747 + ,43 + ,43 + ,733 + ,84105 + ,20 + ,261 + ,45 + ,0 + ,17 + ,17 + ,66 + ,63 + ,20055 + ,5432 + ,34 + ,34 + ,888 + ,77945 + ,20 + ,299 + ,67 + ,0 + ,28 + ,19 + ,76 + ,67 + ,25272 + ,4913 + ,32 + ,31 + ,849 + ,89113 + ,39 + ,300 + ,30 + ,0 + ,19 + ,14 + ,43 + ,32 + ,82206 + ,2650 + ,20 + ,19 + ,1182 + ,91005 + ,29 + ,450 + ,36 + ,3 + ,29 + ,11 + ,39 + ,23 + ,32073 + ,2370 + ,34 + ,34 + ,528 + ,40248 + ,16 + ,183 + ,34 + ,1 + ,8 + ,4 + ,14 + ,7 + ,5444 + ,775 + ,6 + ,6 + ,642 + ,64187 + ,27 + ,238 + ,36 + ,0 + ,10 + ,16 + ,61 + ,54 + ,20154 + ,5576 + ,12 + ,11 + ,947 + ,50857 + ,21 + ,165 + ,34 + ,0 + ,15 + ,20 + ,71 + ,37 + ,36944 + ,1352 + ,24 + ,24 + ,819 + ,56613 + ,19 + ,234 + ,37 + ,1 + ,15 + ,12 + ,44 + ,35 + ,8019 + ,3080 + ,16 + ,16 + ,757 + ,62792 + ,35 + ,176 + ,46 + ,0 + ,28 + ,15 + ,60 + ,51 + ,30884 + ,10205 + ,72 + ,72 + ,894 + ,72535 + ,14 + ,329 + ,44 + ,0 + ,17 + ,16 + ,64 + ,39 + ,19540 + ,6095 + ,27 + ,21) + ,dim=c(14 + ,289) + ,dimnames=list(c('pageviews' + ,'time_in_rfc' + ,'logins' + ,'compendium_views_info' + ,'compendium_views_pr' + ,'shared_compendiums' + ,'blogged_computations' + ,'compendiums_reviewed' + ,'feedback_messages_p1' + ,'feedback_messages_p120' + ,'totsize' + ,'totrevisions' + ,'tothyperlinks' + ,'totblogs ') + ,1:289)) > y <- array(NA,dim=c(14,289),dimnames=list(c('pageviews','time_in_rfc','logins','compendium_views_info','compendium_views_pr','shared_compendiums','blogged_computations','compendiums_reviewed','feedback_messages_p1','feedback_messages_p120','totsize','totrevisions','tothyperlinks','totblogs '),1:289)) > 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 = '2' > 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 Attaching package: 'zoo' The following object(s) are masked from 'package:base': as.Date, as.Date.numeric Loading required package: sandwich Loading required package: strucchange Loading required package: vcd Loading required package: MASS Loading required package: colorspace > library(Hmisc) Hmisc library by Frank E Harrell Jr Type library(help='Hmisc'), ?Overview, or ?Hmisc.Overview') to see overall documentation. NOTE:Hmisc no longer redefines [.factor to drop unused levels when subsetting. To get the old behavior of Hmisc type dropUnusedLevels(). 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] "time_in_rfc" > x[,par1] [1] 210907 120982 176508 179321 123185 52746 385534 33170 101645 149061 [11] 165446 237213 173326 133131 258873 180083 324799 230964 236785 135473 [21] 202925 215147 344297 153935 132943 174724 174415 225548 223632 124817 [31] 221698 210767 170266 260561 84853 294424 101011 215641 325107 7176 [41] 167542 106408 96560 265769 269651 149112 175824 152871 111665 116408 [51] 362301 78800 183167 277965 150629 168809 24188 329267 65029 101097 [61] 218946 244052 341570 103597 233328 256462 206161 311473 235800 177939 [71] 207176 196553 174184 143246 187559 187681 119016 182192 73566 194979 [81] 167488 143756 275541 243199 182999 135649 152299 120221 346485 145790 [91] 193339 80953 122774 130585 112611 286468 241066 148446 204713 182079 [101] 140344 220516 243060 162765 182613 232138 265318 85574 310839 225060 [111] 232317 144966 43287 155754 164709 201940 235454 220801 99466 92661 [121] 133328 61361 125930 100750 224549 82316 102010 101523 243511 22938 [131] 41566 152474 61857 99923 132487 317394 21054 209641 22648 31414 [141] 46698 131698 91735 244749 184510 79863 128423 97839 38214 151101 [151] 272458 172494 108043 328107 250579 351067 158015 98866 85439 229242 [161] 351619 84207 120445 324598 131069 204271 165543 141722 116048 250047 [171] 299775 195838 173260 254488 104389 136084 199476 92499 224330 135781 [181] 74408 81240 14688 181633 271856 7199 46660 17547 133368 95227 [191] 152601 98146 79619 59194 139942 118612 72880 65475 99643 71965 [201] 77272 49289 135131 108446 89746 44296 77648 181528 134019 124064 [211] 92630 121848 52915 81872 58981 53515 60812 56375 65490 80949 [221] 76302 104011 98104 67989 30989 135458 73504 63123 61254 74914 [231] 31774 81437 87186 50090 65745 56653 158399 46455 73624 38395 [241] 91899 139526 52164 51567 70551 84856 102538 86678 85709 34662 [251] 150580 99611 19349 99373 86230 30837 31706 89806 62088 40151 [261] 27634 76990 37460 54157 49862 84337 64175 59382 119308 76702 [271] 103425 70344 43410 104838 62215 69304 53117 19764 86680 84105 [281] 77945 89113 91005 40248 64187 50857 56613 62792 72535 > 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]) 7176 7199 14688 17547 19349 19764 21054 22648 22938 24188 27634 1 1 1 1 1 1 1 1 1 1 1 30837 30989 31414 31706 31774 33170 34662 37460 38214 38395 40151 1 1 1 1 1 1 1 1 1 1 1 40248 41566 43287 43410 44296 46455 46660 46698 49289 49862 50090 1 1 1 1 1 1 1 1 1 1 1 50857 51567 52164 52746 52915 53117 53515 54157 56375 56613 56653 1 1 1 1 1 1 1 1 1 1 1 58981 59194 59382 60812 61254 61361 61857 62088 62215 62792 63123 1 1 1 1 1 1 1 1 1 1 1 64175 64187 65029 65475 65490 65745 67989 69304 70344 70551 71965 1 1 1 1 1 1 1 1 1 1 1 72535 72880 73504 73566 73624 74408 74914 76302 76702 76990 77272 1 1 1 1 1 1 1 1 1 1 1 77648 77945 78800 79619 79863 80949 80953 81240 81437 81872 82316 1 1 1 1 1 1 1 1 1 1 1 84105 84207 84337 84853 84856 85439 85574 85709 86230 86678 86680 1 1 1 1 1 1 1 1 1 1 1 87186 89113 89746 89806 91005 91735 91899 92499 92630 92661 95227 1 1 1 1 1 1 1 1 1 1 1 96560 97839 98104 98146 98866 99373 99466 99611 99643 99923 100750 1 1 1 1 1 1 1 1 1 1 1 101011 101097 101523 101645 102010 102538 103425 103597 104011 104389 104838 1 1 1 1 1 1 1 1 1 1 1 106408 108043 108446 111665 112611 116048 116408 118612 119016 119308 120221 1 1 1 1 1 1 1 1 1 1 1 120445 120982 121848 122774 123185 124064 124817 125930 128423 130585 131069 1 1 1 1 1 1 1 1 1 1 1 131698 132487 132943 133131 133328 133368 134019 135131 135458 135473 135649 1 1 1 1 1 1 1 1 1 1 1 135781 136084 139526 139942 140344 141722 143246 143756 144966 145790 148446 1 1 1 1 1 1 1 1 1 1 1 149061 149112 150580 150629 151101 152299 152474 152601 152871 153935 155754 1 1 1 1 1 1 1 1 1 1 1 158015 158399 162765 164709 165446 165543 167488 167542 168809 170266 172494 1 1 1 1 1 1 1 1 1 1 1 173260 173326 174184 174415 174724 175824 176508 177939 179321 180083 181528 1 1 1 1 1 1 1 1 1 1 1 181633 182079 182192 182613 182999 183167 184510 187559 187681 193339 194979 1 1 1 1 1 1 1 1 1 1 1 195838 196553 199476 201940 202925 204271 204713 206161 207176 209641 210767 1 1 1 1 1 1 1 1 1 1 1 210907 215147 215641 218946 220516 220801 221698 223632 224330 224549 225060 1 1 1 1 1 1 1 1 1 1 1 225548 229242 230964 232138 232317 233328 235454 235800 236785 237213 241066 1 1 1 1 1 1 1 1 1 1 1 243060 243199 243511 244052 244749 250047 250579 254488 256462 258873 260561 1 1 1 1 1 1 1 1 1 1 1 265318 265769 269651 271856 272458 275541 277965 286468 294424 299775 310839 1 1 1 1 1 1 1 1 1 1 1 311473 317394 324598 324799 325107 328107 329267 341570 344297 346485 351067 1 1 1 1 1 1 1 1 1 1 1 351619 362301 385534 1 1 1 > colnames(x) [1] "pageviews" "time_in_rfc" "logins" [4] "compendium_views_info" "compendium_views_pr" "shared_compendiums" [7] "blogged_computations" "compendiums_reviewed" "feedback_messages_p1" [10] "feedback_messages_p120" "totsize" "totrevisions" [13] "tothyperlinks" "totblogs." > colnames(x)[par1] [1] "time_in_rfc" > x[,par1] [1] 210907 120982 176508 179321 123185 52746 385534 33170 101645 149061 [11] 165446 237213 173326 133131 258873 180083 324799 230964 236785 135473 [21] 202925 215147 344297 153935 132943 174724 174415 225548 223632 124817 [31] 221698 210767 170266 260561 84853 294424 101011 215641 325107 7176 [41] 167542 106408 96560 265769 269651 149112 175824 152871 111665 116408 [51] 362301 78800 183167 277965 150629 168809 24188 329267 65029 101097 [61] 218946 244052 341570 103597 233328 256462 206161 311473 235800 177939 [71] 207176 196553 174184 143246 187559 187681 119016 182192 73566 194979 [81] 167488 143756 275541 243199 182999 135649 152299 120221 346485 145790 [91] 193339 80953 122774 130585 112611 286468 241066 148446 204713 182079 [101] 140344 220516 243060 162765 182613 232138 265318 85574 310839 225060 [111] 232317 144966 43287 155754 164709 201940 235454 220801 99466 92661 [121] 133328 61361 125930 100750 224549 82316 102010 101523 243511 22938 [131] 41566 152474 61857 99923 132487 317394 21054 209641 22648 31414 [141] 46698 131698 91735 244749 184510 79863 128423 97839 38214 151101 [151] 272458 172494 108043 328107 250579 351067 158015 98866 85439 229242 [161] 351619 84207 120445 324598 131069 204271 165543 141722 116048 250047 [171] 299775 195838 173260 254488 104389 136084 199476 92499 224330 135781 [181] 74408 81240 14688 181633 271856 7199 46660 17547 133368 95227 [191] 152601 98146 79619 59194 139942 118612 72880 65475 99643 71965 [201] 77272 49289 135131 108446 89746 44296 77648 181528 134019 124064 [211] 92630 121848 52915 81872 58981 53515 60812 56375 65490 80949 [221] 76302 104011 98104 67989 30989 135458 73504 63123 61254 74914 [231] 31774 81437 87186 50090 65745 56653 158399 46455 73624 38395 [241] 91899 139526 52164 51567 70551 84856 102538 86678 85709 34662 [251] 150580 99611 19349 99373 86230 30837 31706 89806 62088 40151 [261] 27634 76990 37460 54157 49862 84337 64175 59382 119308 76702 [271] 103425 70344 43410 104838 62215 69304 53117 19764 86680 84105 [281] 77945 89113 91005 40248 64187 50857 56613 62792 72535 > 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/1x44o1355151582.tab") + } + } > m Conditional inference tree with 12 terminal nodes Response: time_in_rfc Inputs: pageviews, logins, compendium_views_info, compendium_views_pr, shared_compendiums, blogged_computations, compendiums_reviewed, feedback_messages_p1, feedback_messages_p120, totsize, totrevisions, tothyperlinks, totblogs. Number of observations: 289 1) pageviews <= 1613; criterion = 1, statistic = 229.566 2) pageviews <= 1212; criterion = 1, statistic = 126.863 3) pageviews <= 800; criterion = 1, statistic = 72.955 4) pageviews <= 578; criterion = 1, statistic = 31.889 5) pageviews <= 391; criterion = 0.996, statistic = 13.203 6)* weights = 10 5) pageviews > 391 7)* weights = 15 4) pageviews > 578 8) totrevisions <= 3704; criterion = 0.952, statistic = 8.395 9)* weights = 11 8) totrevisions > 3704 10)* weights = 25 3) pageviews > 800 11) tothyperlinks <= 59; criterion = 1, statistic = 24.861 12)* weights = 55 11) tothyperlinks > 59 13)* weights = 9 2) pageviews > 1212 14) totrevisions <= 7936; criterion = 1, statistic = 29.588 15)* weights = 19 14) totrevisions > 7936 16) totblogs. <= 93; criterion = 0.992, statistic = 11.694 17)* weights = 36 16) totblogs. > 93 18)* weights = 15 1) pageviews > 1613 19) compendium_views_info <= 967; criterion = 1, statistic = 46.779 20) feedback_messages_p1 <= 99; criterion = 1, statistic = 26.511 21)* weights = 19 20) feedback_messages_p1 > 99 22)* weights = 58 19) compendium_views_info > 967 23)* weights = 17 > postscript(file="/var/wessaorg/rcomp/tmp/2khps1355151582.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/3qhik1355151582.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 210907 173457.13 37449.8667 2 120982 124193.78 -3211.7778 3 176508 173457.13 3050.8667 4 179321 224956.47 -45635.4655 5 123185 124193.78 -1008.7778 6 52746 39424.13 13321.8667 7 385534 322592.29 62941.7059 8 33170 19387.40 13782.6000 9 101645 86007.40 15637.6000 10 149061 143842.44 5218.5556 11 165446 143842.44 21603.5556 12 237213 224956.47 12256.5345 13 173326 143842.44 29483.5556 14 133131 143842.44 -10711.4444 15 258873 224956.47 33916.5345 16 180083 143842.44 36240.5556 17 324799 322592.29 2206.7059 18 230964 224956.47 6007.5345 19 236785 224956.47 11828.5345 20 135473 124193.78 11279.2222 21 202925 224956.47 -22031.4655 22 215147 224956.47 -9809.4655 23 344297 224956.47 119340.5345 24 153935 143842.44 10092.5556 25 132943 173457.13 -40514.1333 26 174724 224956.47 -50232.4655 27 174415 224956.47 -50541.4655 28 225548 224956.47 591.5345 29 223632 224956.47 -1324.4655 30 124817 98962.42 25854.5789 31 221698 224956.47 -3258.4655 32 210767 224956.47 -14189.4655 33 170266 173457.13 -3191.1333 34 260561 224956.47 35604.5345 35 84853 86007.40 -1154.4000 36 294424 322592.29 -28168.2941 37 101011 86007.40 15003.6000 38 215641 224956.47 -9315.4655 39 325107 224956.47 100150.5345 40 7176 19387.40 -12211.4000 41 167542 143842.44 23699.5556 42 106408 86007.40 20400.6000 43 96560 153513.32 -56953.3158 44 265769 224956.47 40812.5345 45 269651 322592.29 -52941.2941 46 149112 143842.44 5269.5556 47 175824 153513.32 22310.6842 48 152871 143842.44 9028.5556 49 111665 86007.40 25657.6000 50 116408 153513.32 -37105.3158 51 362301 322592.29 39708.7059 52 78800 86007.40 -7207.4000 53 183167 224956.47 -41789.4655 54 277965 322592.29 -44627.2941 55 150629 224956.47 -74327.4655 56 168809 173457.13 -4648.1333 57 24188 39424.13 -15236.1333 58 329267 224956.47 104310.5345 59 65029 47452.55 17576.4545 60 101097 86007.40 15089.6000 61 218946 224956.47 -6010.4655 62 244052 224956.47 19095.5345 63 341570 322592.29 18977.7059 64 103597 86007.40 17589.6000 65 233328 224956.47 8371.5345 66 256462 224956.47 31505.5345 67 206161 224956.47 -18795.4655 68 311473 322592.29 -11119.2941 69 235800 153513.32 82286.6842 70 177939 224956.47 -47017.4655 71 207176 224956.47 -17780.4655 72 196553 173457.13 23095.8667 73 174184 143842.44 30341.5556 74 143246 173457.13 -30211.1333 75 187559 224956.47 -37397.4655 76 187681 224956.47 -37275.4655 77 119016 143842.44 -24826.4444 78 182192 224956.47 -42764.4655 79 73566 86007.40 -12441.4000 80 194979 224956.47 -29977.4655 81 167488 173457.13 -5969.1333 82 143756 173457.13 -29701.1333 83 275541 224956.47 50584.5345 84 243199 224956.47 18242.5345 85 182999 143842.44 39156.5556 86 135649 143842.44 -8193.4444 87 152299 143842.44 8456.5556 88 120221 143842.44 -23621.4444 89 346485 322592.29 23892.7059 90 145790 143842.44 1947.5556 91 193339 153513.32 39825.6842 92 80953 86007.40 -5054.4000 93 122774 153513.32 -30739.3158 94 130585 124193.78 6391.2222 95 112611 86007.40 26603.6000 96 286468 322592.29 -36124.2941 97 241066 224956.47 16109.5345 98 148446 224956.47 -76510.4655 99 204713 224956.47 -20243.4655 100 182079 224956.47 -42877.4655 101 140344 98962.42 41381.5789 102 220516 224956.47 -4440.4655 103 243060 224956.47 18103.5345 104 162765 173457.13 -10692.1333 105 182613 173457.13 9155.8667 106 232138 224956.47 7181.5345 107 265318 224956.47 40361.5345 108 85574 98962.42 -13388.4211 109 310839 322592.29 -11753.2941 110 225060 224956.47 103.5345 111 232317 224956.47 7360.5345 112 144966 143842.44 1123.5556 113 43287 65732.56 -22445.5600 114 155754 153513.32 2240.6842 115 164709 173457.13 -8748.1333 116 201940 224956.47 -23016.4655 117 235454 224956.47 10497.5345 118 220801 153513.32 67287.6842 119 99466 86007.40 13458.6000 120 92661 98962.42 -6301.4211 121 133328 143842.44 -10514.4444 122 61361 86007.40 -24646.4000 123 125930 153513.32 -27583.3158 124 100750 98962.42 1787.5789 125 224549 173457.13 51091.8667 126 82316 86007.40 -3691.4000 127 102010 86007.40 16002.6000 128 101523 143842.44 -42319.4444 129 243511 224956.47 18554.5345 130 22938 19387.40 3550.6000 131 41566 47452.55 -5886.5455 132 152474 173457.13 -20983.1333 133 61857 39424.13 22432.8667 134 99923 153513.32 -53590.3158 135 132487 143842.44 -11355.4444 136 317394 322592.29 -5198.2941 137 21054 19387.40 1666.6000 138 209641 153513.32 56127.6842 139 22648 47452.55 -24804.5455 140 31414 39424.13 -8010.1333 141 46698 47452.55 -754.5455 142 131698 153513.32 -21815.3158 143 91735 86007.40 5727.6000 144 244749 224956.47 19792.5345 145 184510 143842.44 40667.5556 146 79863 98962.42 -19099.4211 147 128423 86007.40 42415.6000 148 97839 124193.78 -26354.7778 149 38214 39424.13 -1210.1333 150 151101 143842.44 7258.5556 151 272458 224956.47 47501.5345 152 172494 124193.78 48300.2222 153 108043 98962.42 9080.5789 154 328107 322592.29 5514.7059 155 250579 224956.47 25622.5345 156 351067 322592.29 28474.7059 157 158015 124193.78 33821.2222 158 98866 86007.40 12858.6000 159 85439 86007.40 -568.4000 160 229242 224956.47 4285.5345 161 351619 322592.29 29026.7059 162 84207 86007.40 -1800.4000 163 120445 143842.44 -23397.4444 164 324598 322592.29 2005.7059 165 131069 143842.44 -12773.4444 166 204271 173457.13 30813.8667 167 165543 224956.47 -59413.4655 168 141722 143842.44 -2120.4444 169 116048 124193.78 -8145.7778 170 250047 153513.32 96533.6842 171 299775 322592.29 -22817.2941 172 195838 224956.47 -29118.4655 173 173260 153513.32 19746.6842 174 254488 224956.47 29531.5345 175 104389 143842.44 -39453.4444 176 136084 143842.44 -7758.4444 177 199476 224956.47 -25480.4655 178 92499 86007.40 6491.6000 179 224330 224956.47 -626.4655 180 135781 143842.44 -8061.4444 181 74408 98962.42 -24554.4211 182 81240 98962.42 -17722.4211 183 14688 19387.40 -4699.4000 184 181633 224956.47 -43323.4655 185 271856 224956.47 46899.5345 186 7199 19387.40 -12188.4000 187 46660 39424.13 7235.8667 188 17547 19387.40 -1840.4000 189 133368 153513.32 -20145.3158 190 95227 86007.40 9219.6000 191 152601 143842.44 8758.5556 192 98146 86007.40 12138.6000 193 79619 98962.42 -19343.4211 194 59194 98962.42 -39768.4211 195 139942 143842.44 -3900.4444 196 118612 98962.42 19649.5789 197 72880 86007.40 -13127.4000 198 65475 65732.56 -257.5600 199 99643 98962.42 680.5789 200 71965 65732.56 6232.4400 201 77272 65732.56 11539.4400 202 49289 39424.13 9864.8667 203 135131 143842.44 -8711.4444 204 108446 143842.44 -35396.4444 205 89746 98962.42 -9216.4211 206 44296 39424.13 4871.8667 207 77648 86007.40 -8359.4000 208 181528 143842.44 37685.5556 209 134019 143842.44 -9823.4444 210 124064 86007.40 38056.6000 211 92630 86007.40 6622.6000 212 121848 98962.42 22885.5789 213 52915 65732.56 -12817.5600 214 81872 65732.56 16139.4400 215 58981 86007.40 -27026.4000 216 53515 65732.56 -12217.5600 217 60812 86007.40 -25195.4000 218 56375 65732.56 -9357.5600 219 65490 65732.56 -242.5600 220 80949 65732.56 15216.4400 221 76302 86007.40 -9705.4000 222 104011 143842.44 -39831.4444 223 98104 153513.32 -55409.3158 224 67989 65732.56 2256.4400 225 30989 19387.40 11601.6000 226 135458 153513.32 -18055.3158 227 73504 47452.55 26051.4545 228 63123 124193.78 -61070.7778 229 61254 86007.40 -24753.4000 230 74914 98962.42 -24048.4211 231 31774 47452.55 -15678.5455 232 81437 86007.40 -4570.4000 233 87186 98962.42 -11776.4211 234 50090 65732.56 -15642.5600 235 65745 86007.40 -20262.4000 236 56653 65732.56 -9079.5600 237 158399 98962.42 59436.5789 238 46455 47452.55 -997.5455 239 73624 86007.40 -12383.4000 240 38395 47452.55 -9057.5455 241 91899 86007.40 5891.6000 242 139526 153513.32 -13987.3158 243 52164 65732.56 -13568.5600 244 51567 65732.56 -14165.5600 245 70551 86007.40 -15456.4000 246 84856 86007.40 -1151.4000 247 102538 153513.32 -50975.3158 248 86678 65732.56 20945.4400 249 85709 86007.40 -298.4000 250 34662 39424.13 -4762.1333 251 150580 143842.44 6737.5556 252 99611 86007.40 13603.6000 253 19349 19387.40 -38.4000 254 99373 86007.40 13365.6000 255 86230 86007.40 222.6000 256 30837 39424.13 -8587.1333 257 31706 39424.13 -7718.1333 258 89806 86007.40 3798.6000 259 62088 65732.56 -3644.5600 260 40151 39424.13 726.8667 261 27634 39424.13 -11790.1333 262 76990 86007.40 -9017.4000 263 37460 39424.13 -1964.1333 264 54157 65732.56 -11575.5600 265 49862 65732.56 -15870.5600 266 84337 86007.40 -1670.4000 267 64175 86007.40 -21832.4000 268 59382 47452.55 11929.4545 269 119308 65732.56 53575.4400 270 76702 86007.40 -9305.4000 271 103425 98962.42 4462.5789 272 70344 65732.56 4611.4400 273 43410 47452.55 -4042.5455 274 104838 86007.40 18830.6000 275 62215 65732.56 -3517.5600 276 69304 86007.40 -16703.4000 277 53117 47452.55 5664.4545 278 19764 19387.40 376.6000 279 86680 86007.40 672.6000 280 84105 65732.56 18372.4400 281 77945 86007.40 -8062.4000 282 89113 86007.40 3105.6000 283 91005 86007.40 4997.6000 284 40248 39424.13 823.8667 285 64187 65732.56 -1545.5600 286 50857 86007.40 -35150.4000 287 56613 86007.40 -29394.4000 288 62792 65732.56 -2940.5600 289 72535 86007.40 -13472.4000 > 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/4fq8y1355151583.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/5zq7n1355151583.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/6ngqe1355151583.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/7eaua1355151583.tab") + } > > try(system("convert tmp/2khps1355151582.ps tmp/2khps1355151582.png",intern=TRUE)) character(0) > try(system("convert tmp/3qhik1355151582.ps tmp/3qhik1355151582.png",intern=TRUE)) character(0) > try(system("convert tmp/4fq8y1355151583.ps tmp/4fq8y1355151583.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 10.911 0.494 11.399