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(30 + ,210907 + ,56 + ,79 + ,112285 + ,1418 + ,396 + ,81 + ,24.188 + ,144 + ,28 + ,120982 + ,56 + ,58 + ,84786 + ,869 + ,297 + ,55 + ,18.273 + ,103 + ,38 + ,176508 + ,54 + ,60 + ,83123 + ,1530 + ,559 + ,50 + ,14.130 + ,98 + ,30 + ,179321 + ,89 + ,108 + ,101193 + ,2172 + ,967 + ,125 + ,32287 + ,135 + ,22 + ,123185 + ,40 + ,49 + ,38361 + ,901 + ,270 + ,40 + ,8654 + ,61 + ,26 + ,52746 + ,25 + ,0 + ,68504 + ,463 + ,143 + ,37 + ,9245 + ,39 + ,25 + ,385534 + ,92 + ,121 + ,119182 + ,3201 + ,1562 + ,63 + ,33251 + ,150 + ,18 + ,33170 + ,18 + ,1 + ,22807 + ,371 + ,109 + ,44 + ,1271 + ,5 + ,11 + ,101645 + ,63 + ,20 + ,17140 + ,1192 + ,371 + ,88 + ,5279 + ,28 + ,26 + ,149061 + ,44 + ,43 + ,116174 + ,1583 + ,656 + ,66 + ,27.101 + ,84 + ,25 + ,165446 + ,33 + ,69 + ,57635 + ,1439 + ,511 + ,57 + ,16.373 + ,80 + ,38 + ,237213 + ,84 + ,78 + ,66198 + ,1764 + ,655 + ,74 + ,19.716 + ,130 + ,44 + ,173326 + ,88 + ,86 + ,71701 + ,1495 + ,465 + ,49 + ,17753 + ,82 + ,30 + ,133131 + ,55 + ,44 + ,57793 + ,1373 + ,525 + ,52 + ,9.028 + ,60 + ,40 + ,258873 + ,60 + ,104 + ,80444 + ,2187 + ,885 + ,88 + ,18653 + ,131 + ,34 + ,180083 + ,66 + ,63 + ,53855 + ,1491 + ,497 + ,36 + ,8828 + ,84 + ,47 + ,324799 + ,154 + ,158 + ,97668 + ,4041 + ,1436 + ,108 + ,29498 + ,140 + ,30 + ,230964 + ,53 + ,102 + ,133824 + ,1706 + ,612 + ,43 + ,27563 + ,151 + ,31 + ,236785 + ,119 + ,77 + ,101481 + ,2152 + ,865 + ,75 + ,18.293 + ,91 + ,23 + ,135473 + ,41 + ,82 + ,99645 + ,1036 + ,385 + ,32 + ,22.530 + ,138 + ,36 + ,202925 + ,61 + ,115 + ,114789 + ,1882 + ,567 + ,44 + ,15977 + ,150 + ,36 + ,215147 + ,58 + ,101 + ,99052 + ,1929 + ,639 + ,85 + ,35.082 + ,124 + ,30 + ,344297 + ,75 + ,80 + ,67654 + ,2242 + ,963 + ,86 + ,16.116 + ,119 + ,25 + ,153935 + ,33 + ,50 + ,65553 + ,1220 + ,398 + ,56 + ,15.849 + ,73 + ,39 + ,132943 + ,40 + ,83 + ,97500 + ,1289 + ,410 + ,50 + ,16026 + ,110 + ,34 + ,174724 + ,92 + ,123 + ,69112 + ,2515 + ,966 + ,135 + ,26569 + ,123 + ,31 + ,174415 + ,100 + ,73 + ,82753 + ,2147 + ,801 + ,63 + ,24.785 + ,90 + ,31 + ,225548 + ,112 + ,81 + ,85323 + ,2352 + ,892 + ,81 + ,17.569 + ,116 + ,33 + ,223632 + ,73 + ,105 + ,72654 + ,1638 + ,513 + ,52 + ,23.825 + ,113 + ,25 + ,124817 + ,40 + ,47 + ,30727 + ,1222 + ,469 + ,44 + ,7869 + ,56 + ,33 + ,221698 + ,45 + ,105 + ,77873 + ,1812 + ,683 + ,113 + ,14975 + ,115 + ,35 + ,210767 + ,60 + ,94 + ,117478 + ,1677 + ,643 + ,39 + ,37.791 + ,119 + ,42 + ,170266 + ,62 + ,44 + ,74007 + ,1579 + ,535 + ,73 + ,9.605 + ,129 + ,43 + ,260561 + ,75 + ,114 + ,90183 + ,1731 + ,625 + ,48 + ,27295 + ,127 + ,30 + ,84853 + ,31 + ,38 + ,61542 + ,807 + ,264 + ,33 + ,2746 + ,27 + ,33 + ,294424 + ,77 + ,107 + ,101494 + ,2452 + ,992 + ,59 + ,34.461 + ,175 + ,13 + ,101011 + ,34 + ,30 + ,27570 + ,829 + ,238 + ,41 + ,8098 + ,35 + ,32 + ,215641 + ,46 + ,71 + ,55813 + ,1940 + ,818 + ,69 + ,4787 + ,64 + ,36 + ,325107 + ,99 + ,84 + ,79215 + ,2662 + ,937 + ,64 + ,24.919 + ,96 + ,0 + ,7176 + ,17 + ,0 + ,1423 + ,186 + ,70 + ,1 + ,603 + ,0 + ,28 + ,167542 + ,66 + ,59 + ,55461 + ,1499 + ,507 + ,59 + ,16329 + ,84 + ,14 + ,106408 + ,30 + ,33 + ,31081 + ,865 + ,260 + ,32 + ,12558 + ,41 + ,17 + ,96560 + ,76 + ,42 + ,22996 + ,1793 + ,503 + ,129 + ,7.784 + ,47 + ,32 + ,265769 + ,146 + ,96 + ,83122 + ,2527 + ,927 + ,37 + ,28.522 + ,126 + ,30 + ,269651 + ,67 + ,106 + ,70106 + ,2747 + ,1269 + ,31 + ,22265 + ,105 + ,35 + ,149112 + ,56 + ,56 + ,60578 + ,1324 + ,537 + ,65 + ,14.459 + ,80 + ,20 + ,175824 + ,107 + ,57 + ,39992 + ,2702 + ,910 + ,107 + ,14.526 + ,70 + ,28 + ,152871 + ,58 + ,59 + ,79892 + ,1383 + ,532 + ,74 + ,22.240 + ,73 + ,28 + ,111665 + ,34 + ,39 + ,49810 + ,1179 + ,345 + ,54 + ,11.802 + ,57 + ,39 + ,116408 + ,61 + ,34 + ,71570 + ,2099 + ,918 + ,76 + ,7623 + ,40 + ,34 + ,362301 + ,119 + ,76 + ,100708 + ,4308 + ,1635 + ,715 + ,11.912 + ,68 + ,26 + ,78800 + ,42 + ,20 + ,33032 + ,918 + ,330 + ,57 + ,7935 + ,21 + ,39 + ,183167 + ,66 + ,91 + ,82875 + ,1831 + ,557 + ,66 + ,18.220 + ,127 + ,39 + ,277965 + ,89 + ,115 + ,139077 + ,3373 + ,1178 + ,106 + ,19199 + ,154 + ,33 + ,150629 + ,44 + ,85 + ,71595 + ,1713 + ,740 + ,54 + ,19918 + ,116 + ,28 + ,168809 + ,66 + ,76 + ,72260 + ,1438 + ,452 + ,32 + ,21.884 + ,102 + ,4 + ,24188 + ,24 + ,8 + ,5950 + ,496 + ,218 + ,20 + ,2.694 + ,7 + ,39 + ,329267 + ,259 + ,79 + ,115762 + ,2253 + ,764 + ,71 + ,15.808 + ,148 + ,18 + ,65029 + ,17 + ,21 + ,32551 + ,744 + ,255 + ,21 + ,3597 + ,21 + ,14 + ,101097 + ,64 + ,30 + ,31701 + ,1161 + ,454 + ,70 + ,5296 + ,35 + ,29 + ,218946 + ,41 + ,76 + ,80670 + ,2352 + ,866 + ,112 + ,25.239 + ,112 + ,44 + ,244052 + ,68 + ,101 + ,143558 + ,2144 + ,574 + ,66 + ,29.801 + ,137 + ,21 + ,341570 + ,168 + ,94 + ,117105 + ,4691 + ,1276 + ,190 + ,18.450 + ,135 + ,16 + ,103597 + ,43 + ,27 + ,23789 + ,1112 + ,379 + ,66 + ,7.132 + ,26 + ,28 + ,233328 + ,132 + ,92 + ,120733 + ,2694 + ,825 + ,165 + ,34861 + ,230 + ,35 + ,256462 + ,105 + ,123 + ,105195 + ,1973 + ,798 + ,56 + ,35.940 + ,181 + ,28 + ,206161 + ,71 + ,75 + ,73107 + ,1769 + ,663 + ,61 + ,16688 + ,71 + ,38 + ,311473 + ,112 + ,128 + ,132068 + ,3148 + ,1069 + ,53 + ,24683 + ,147 + ,23 + ,235800 + ,94 + ,105 + ,149193 + ,2474 + ,921 + ,127 + ,46.230 + ,190 + ,36 + ,177939 + ,82 + ,55 + ,46821 + ,2084 + ,858 + ,63 + ,10387 + ,64 + ,32 + ,207176 + ,70 + ,56 + ,87011 + ,1954 + ,711 + ,38 + ,21436 + ,105 + ,29 + ,196553 + ,57 + ,41 + ,95260 + ,1226 + ,503 + ,50 + ,30.546 + ,107 + ,25 + ,174184 + ,53 + ,72 + ,55183 + ,1389 + ,382 + ,52 + ,19.746 + ,94 + ,27 + ,143246 + ,103 + ,67 + ,106671 + ,1496 + ,464 + ,42 + ,15.977 + ,116 + ,36 + ,187559 + ,121 + ,75 + ,73511 + ,2269 + ,717 + ,76 + ,22.583 + ,106 + ,28 + ,187681 + ,62 + ,114 + ,92945 + ,1833 + ,690 + ,67 + ,17.274 + ,143 + ,23 + ,119016 + ,52 + ,118 + ,78664 + ,1268 + ,462 + ,50 + ,16469 + ,81 + ,40 + ,182192 + ,52 + ,77 + ,70054 + ,1943 + ,657 + ,53 + ,14251 + ,89 + ,23 + ,73566 + ,32 + ,22 + ,22618 + ,893 + ,385 + ,39 + ,3.007 + ,26 + ,40 + ,194979 + ,62 + ,66 + ,74011 + ,1762 + ,577 + ,50 + ,16851 + ,84 + ,28 + ,167488 + ,45 + ,69 + ,83737 + ,1403 + ,619 + ,77 + ,21.113 + ,113 + ,34 + ,143756 + ,46 + ,105 + ,69094 + ,1425 + ,479 + ,57 + ,17401 + ,120 + ,33 + ,275541 + ,63 + ,116 + ,93133 + ,1857 + ,817 + ,73 + ,23958 + ,110 + ,28 + ,243199 + ,75 + ,88 + ,95536 + ,1840 + ,752 + ,34 + ,23567 + ,134 + ,34 + ,182999 + ,88 + ,73 + ,225920 + ,1502 + ,430 + ,39 + ,13.065 + ,54 + ,30 + ,135649 + ,46 + ,99 + ,62133 + ,1441 + ,451 + ,46 + ,15358 + ,96 + ,33 + ,152299 + ,53 + ,62 + ,61370 + ,1420 + ,537 + ,63 + ,14.587 + ,78 + ,22 + ,120221 + ,37 + ,53 + ,43836 + ,1416 + ,519 + ,35 + ,12770 + ,51 + ,38 + ,346485 + ,90 + ,118 + ,106117 + ,2970 + ,1000 + ,106 + ,24.021 + ,121 + ,26 + ,145790 + ,63 + ,30 + ,38692 + ,1317 + ,637 + ,43 + ,9648 + ,38 + ,35 + ,193339 + ,78 + ,100 + ,84651 + ,1644 + ,465 + ,47 + ,20.537 + ,145 + ,8 + ,80953 + ,25 + ,49 + ,56622 + ,870 + ,437 + ,31 + ,7905 + ,59 + ,24 + ,122774 + ,45 + ,24 + ,15986 + ,1654 + ,711 + ,162 + ,4.527 + ,27 + ,29 + ,130585 + ,46 + ,67 + ,95364 + ,1054 + ,299 + ,57 + ,30.495 + ,91 + ,20 + ,112611 + ,41 + ,46 + ,26706 + ,937 + ,248 + ,36 + ,7.117 + ,48 + ,29 + ,286468 + ,144 + ,57 + ,89691 + ,3004 + ,1162 + ,263 + ,17719 + ,68 + ,45 + ,241066 + ,82 + ,75 + ,67267 + ,2008 + ,714 + ,78 + ,27056 + ,58 + ,37 + ,148446 + ,91 + ,135 + ,126846 + ,2547 + ,905 + ,63 + ,33.473 + ,150 + ,33 + ,204713 + ,71 + ,68 + ,41140 + ,1885 + ,649 + ,54 + ,9758 + ,74 + ,33 + ,182079 + ,63 + ,124 + ,102860 + ,1626 + ,512 + ,63 + ,21.115 + ,181 + ,25 + ,140344 + ,53 + ,33 + ,51715 + ,1468 + ,472 + ,77 + ,7.236 + ,65 + ,32 + ,220516 + ,62 + ,98 + ,55801 + ,2445 + ,905 + ,79 + ,13.790 + ,97 + ,29 + ,243060 + ,63 + ,58 + ,111813 + ,1964 + ,786 + ,110 + ,32.902 + ,121 + ,28 + ,162765 + ,32 + ,68 + ,120293 + ,1381 + ,489 + ,56 + ,25.131 + ,99 + ,28 + ,182613 + ,39 + ,81 + ,138599 + ,1369 + ,479 + ,56 + ,30910 + ,152 + ,31 + ,232138 + ,62 + ,131 + ,161647 + ,1659 + ,617 + ,43 + ,35.947 + ,188 + ,52 + ,265318 + ,117 + ,110 + ,115929 + ,2888 + ,925 + ,111 + ,29.848 + ,138 + ,21 + ,85574 + ,34 + ,37 + ,24266 + ,1290 + ,351 + ,71 + ,6943 + ,40 + ,24 + ,310839 + ,92 + ,130 + ,162901 + ,2845 + ,1144 + ,62 + ,42.705 + ,254 + ,41 + ,225060 + ,93 + ,93 + ,109825 + ,1982 + ,669 + ,56 + ,31.808 + ,87 + ,33 + ,232317 + ,54 + ,118 + ,129838 + ,1904 + ,707 + ,74 + ,26.675 + ,178 + ,32 + ,144966 + ,144 + ,39 + ,37510 + ,1391 + ,458 + ,60 + ,8.435 + ,51 + ,19 + ,43287 + ,14 + ,13 + ,43750 + ,602 + ,214 + ,43 + ,7409 + ,49 + ,20 + ,155754 + ,61 + ,74 + ,40652 + ,1743 + ,599 + ,68 + ,14993 + ,73 + ,31 + ,164709 + ,109 + ,81 + ,87771 + ,1559 + ,572 + ,53 + ,36.867 + ,176 + ,31 + ,201940 + ,38 + ,109 + ,85872 + ,2014 + ,897 + ,87 + ,33835 + ,94 + ,32 + ,235454 + ,73 + ,151 + ,89275 + ,2143 + ,819 + ,46 + ,24164 + ,120 + ,18 + ,220801 + ,75 + ,51 + ,44418 + ,2146 + ,720 + ,105 + ,12.607 + ,66 + ,23 + ,99466 + ,50 + ,28 + ,192565 + ,874 + ,273 + ,32 + ,22609 + ,56 + ,17 + ,92661 + ,61 + ,40 + ,35232 + ,1590 + ,508 + ,133 + ,5.892 + ,39 + ,20 + ,133328 + ,55 + ,56 + ,40909 + ,1590 + ,506 + ,79 + ,17.014 + ,66 + ,12 + ,61361 + ,77 + ,27 + ,13294 + ,1210 + ,451 + ,51 + ,5.394 + ,27 + ,17 + ,125930 + ,75 + ,37 + ,32387 + ,2072 + ,699 + ,207 + ,9178 + ,65 + ,30 + ,100750 + ,72 + ,83 + ,140867 + ,1281 + ,407 + ,67 + ,6.440 + ,58 + ,31 + ,224549 + ,50 + ,54 + ,120662 + ,1401 + ,465 + ,47 + ,21916 + ,98 + ,10 + ,82316 + ,32 + ,27 + ,21233 + ,834 + ,245 + ,34 + ,4011 + ,25 + ,13 + ,102010 + ,53 + ,28 + ,44332 + ,1105 + ,370 + ,66 + ,5.818 + ,26 + ,22 + ,101523 + ,42 + ,59 + ,61056 + ,1272 + ,316 + ,76 + ,18.647 + ,77 + ,42 + ,243511 + ,71 + ,133 + ,101338 + ,1944 + ,603 + ,65 + ,20.556 + ,130 + ,1 + ,22938 + ,10 + ,12 + ,1168 + ,391 + ,154 + ,9 + ,238 + ,11 + ,9 + ,41566 + ,35 + ,0 + ,13497 + ,761 + ,229 + ,42 + ,70 + ,2 + ,32 + ,152474 + ,65 + ,106 + ,65567 + ,1605 + ,577 + ,45 + ,22.392 + ,101 + ,11 + ,61857 + ,25 + ,23 + ,25162 + ,530 + ,192 + ,25 + ,3913 + ,31 + ,25 + ,99923 + ,66 + ,44 + ,32334 + ,1988 + ,617 + ,115 + ,12.237 + ,36 + ,36 + ,132487 + ,41 + ,71 + ,40735 + ,1386 + ,411 + ,97 + ,8.388 + ,120 + ,31 + ,317394 + ,86 + ,116 + ,91413 + ,2395 + ,975 + ,53 + ,22.120 + ,195 + ,0 + ,21054 + ,16 + ,4 + ,855 + ,387 + ,146 + ,2 + ,338 + ,4 + ,24 + ,209641 + ,42 + ,62 + ,97068 + ,1742 + ,705 + ,52 + ,11727 + ,89 + ,13 + ,22648 + ,19 + ,12 + ,44339 + ,620 + ,184 + ,44 + ,3.704 + ,24 + ,8 + ,31414 + ,19 + ,18 + ,14116 + ,449 + ,200 + ,22 + ,3.988 + ,39 + ,13 + ,46698 + ,45 + ,14 + ,10288 + ,800 + ,274 + ,35 + ,3.030 + ,14 + ,19 + ,131698 + ,65 + ,60 + ,65622 + ,1684 + ,502 + ,74 + ,13.520 + ,78 + ,18 + ,91735 + ,35 + ,7 + ,16563 + ,1050 + ,382 + ,103 + ,1421 + ,15 + ,33 + ,244749 + ,95 + ,98 + ,76643 + ,2699 + ,964 + ,144 + ,20.923 + ,106 + ,40 + ,184510 + ,49 + ,64 + ,110681 + ,1606 + ,537 + ,60 + ,20237 + ,83 + ,22 + ,79863 + ,37 + ,29 + ,29011 + ,1502 + ,438 + ,134 + ,3219 + ,24 + ,38 + ,128423 + ,64 + ,32 + ,92696 + ,1204 + ,369 + ,89 + ,3.769 + ,37 + ,24 + ,97839 + ,38 + ,25 + ,94785 + ,1138 + ,417 + ,42 + ,12.252 + ,77 + ,8 + ,38214 + ,34 + ,16 + ,8773 + ,568 + ,276 + ,52 + ,1888 + ,16 + ,35 + ,151101 + ,32 + ,48 + ,83209 + ,1459 + ,514 + ,98 + ,14497 + ,56 + ,43 + ,272458 + ,65 + ,100 + ,93815 + ,2158 + ,822 + ,99 + ,28864 + ,132 + ,43 + ,172494 + ,52 + ,46 + ,86687 + ,1111 + ,389 + ,52 + ,21.721 + ,144 + ,14 + ,108043 + ,62 + ,45 + ,34553 + ,1421 + ,466 + ,29 + ,4821 + ,40 + ,41 + ,328107 + ,65 + ,129 + ,105547 + ,2833 + ,1255 + ,125 + ,33.644 + ,153 + ,38 + ,250579 + ,83 + ,130 + ,103487 + ,1955 + ,694 + ,106 + ,15923 + ,143 + ,45 + ,351067 + ,95 + ,136 + ,213688 + ,2922 + ,1024 + ,95 + ,42935 + ,220 + ,31 + ,158015 + ,29 + ,59 + ,71220 + ,1002 + ,400 + ,40 + ,18.864 + ,79 + ,13 + ,98866 + ,18 + ,25 + ,23517 + ,1060 + ,397 + ,140 + ,4977 + ,50 + ,28 + ,85439 + ,33 + ,32 + ,56926 + ,956 + ,350 + ,43 + ,7785 + ,39 + ,31 + ,229242 + ,247 + ,63 + ,91721 + ,2186 + ,719 + ,128 + ,17.939 + ,95 + ,40 + ,351619 + ,139 + ,95 + ,115168 + ,3604 + ,1277 + ,142 + ,23436 + ,169 + ,30 + ,84207 + ,29 + ,14 + ,111194 + ,1035 + ,356 + ,73 + ,325 + ,12 + ,16 + ,120445 + ,118 + ,36 + ,51009 + ,1417 + ,457 + ,72 + ,13.539 + ,63 + ,37 + ,324598 + ,110 + ,113 + ,135777 + ,3261 + ,1402 + ,128 + ,34.538 + ,134 + ,30 + ,131069 + ,67 + ,47 + ,51513 + ,1587 + ,600 + ,61 + ,12.198 + ,69 + ,35 + ,204271 + ,42 + ,92 + ,74163 + ,1424 + ,480 + ,73 + ,26924 + ,119 + ,32 + ,165543 + ,65 + ,70 + ,51633 + ,1701 + ,595 + ,148 + ,12716 + ,119 + ,27 + ,141722 + ,94 + ,19 + ,75345 + ,1249 + ,436 + ,64 + ,8172 + ,75 + ,20 + ,116048 + ,64 + ,50 + ,33416 + ,946 + ,230 + ,45 + ,10.855 + ,63 + ,18 + ,250047 + ,81 + ,41 + ,83305 + ,1926 + ,651 + ,58 + ,11.932 + ,55 + ,31 + ,299775 + ,95 + ,91 + ,98952 + ,3352 + ,1367 + ,97 + ,14.300 + ,103 + ,31 + ,195838 + ,67 + ,111 + ,102372 + ,1641 + ,564 + ,50 + ,25.515 + ,197 + ,21 + ,173260 + ,63 + ,41 + ,37238 + ,2035 + ,716 + ,37 + ,2.805 + ,16 + ,39 + ,254488 + ,83 + ,120 + ,103772 + ,2312 + ,747 + ,50 + ,29.402 + ,140 + ,41 + ,104389 + ,45 + ,135 + ,123969 + ,1369 + ,467 + ,105 + ,16440 + ,89 + ,13 + ,136084 + ,30 + ,27 + ,27142 + ,1577 + ,671 + ,69 + ,11221 + ,40 + ,32 + ,199476 + ,70 + ,87 + ,135400 + ,2201 + ,861 + ,46 + ,28732 + ,125 + ,18 + ,92499 + ,32 + ,25 + ,21399 + ,961 + ,319 + ,57 + ,5.250 + ,21 + ,39 + ,224330 + ,83 + ,131 + ,130115 + ,1900 + ,612 + ,52 + ,28.608 + ,167 + ,14 + ,135781 + ,31 + ,45 + ,24874 + ,1254 + ,433 + ,98 + ,8.092 + ,32 + ,7 + ,74408 + ,67 + ,29 + ,34988 + ,1335 + ,434 + ,61 + ,4.473 + ,36 + ,17 + ,81240 + ,66 + ,58 + ,45549 + ,1597 + ,503 + ,89 + ,1.572 + ,13 + ,0 + ,14688 + ,10 + ,4 + ,6023 + ,207 + ,85 + ,0 + ,2065 + ,5 + ,30 + ,181633 + ,70 + ,47 + ,64466 + ,1645 + ,564 + ,48 + ,14.817 + ,96 + ,37 + ,271856 + ,103 + ,109 + ,54990 + ,2429 + ,824 + ,91 + ,16714 + ,151 + ,0 + ,7199 + ,5 + ,7 + ,1644 + ,151 + ,74 + ,0 + ,556 + ,6 + ,5 + ,46660 + ,20 + ,12 + ,6179 + ,474 + ,259 + ,7 + ,2089 + ,13 + ,1 + ,17547 + ,5 + ,0 + ,3926 + ,141 + ,69 + ,3 + ,2658 + ,3 + ,16 + ,133368 + ,36 + ,37 + ,32755 + ,1639 + ,535 + ,54 + ,10695 + ,57 + ,32 + ,95227 + ,34 + ,37 + ,34777 + ,872 + ,239 + ,70 + ,1.669 + ,23 + ,24 + ,152601 + ,48 + ,46 + ,73224 + ,1318 + ,438 + ,36 + ,16267 + ,61 + ,17 + ,98146 + ,40 + ,15 + ,27114 + ,1018 + ,459 + ,37 + ,7.768 + ,21 + ,11 + ,79619 + ,43 + ,42 + ,20760 + ,1383 + ,426 + ,123 + ,7252 + ,43 + ,24 + ,59194 + ,31 + ,7 + ,37636 + ,1314 + ,288 + ,247 + ,6.387 + ,20 + ,22 + ,139942 + ,42 + ,54 + ,65461 + ,1335 + ,498 + ,46 + ,18.715 + ,82 + ,12 + ,118612 + ,46 + ,54 + ,30080 + ,1403 + ,454 + ,72 + ,7.936 + ,90 + ,19 + ,72880 + ,33 + ,14 + ,24094 + ,910 + ,376 + ,41 + ,8.643 + ,25 + ,13 + ,65475 + ,18 + ,16 + ,69008 + ,616 + ,225 + ,24 + ,7.294 + ,60 + ,17 + ,99643 + ,55 + ,33 + ,54968 + ,1407 + ,555 + ,45 + ,4570 + ,61 + ,15 + ,71965 + ,35 + ,32 + ,46090 + ,771 + ,252 + ,33 + ,7.185 + ,85 + ,16 + ,77272 + ,59 + ,21 + ,27507 + ,766 + ,208 + ,27 + ,10058 + ,43 + ,24 + ,49289 + ,19 + ,15 + ,10672 + ,473 + ,130 + ,36 + ,2342 + ,25 + ,15 + ,135131 + ,66 + ,38 + ,34029 + ,1376 + ,481 + ,87 + ,8.509 + ,41 + ,17 + ,108446 + ,60 + ,22 + ,46300 + ,1232 + ,389 + ,90 + ,13.275 + ,26 + ,18 + ,89746 + ,36 + ,28 + ,24760 + ,1521 + ,565 + ,114 + ,6816 + ,38 + ,20 + ,44296 + ,25 + ,10 + ,18779 + ,572 + ,173 + ,31 + ,1930 + ,12 + ,16 + ,77648 + ,47 + ,31 + ,21280 + ,1059 + ,278 + ,45 + ,8086 + ,29 + ,16 + ,181528 + ,54 + ,32 + ,40662 + ,1544 + ,609 + ,69 + ,10.737 + ,49 + ,18 + ,134019 + ,53 + ,32 + ,28987 + ,1230 + ,422 + ,51 + ,8.033 + ,46 + ,22 + ,124064 + ,40 + ,43 + ,22827 + ,1206 + ,445 + ,34 + ,7058 + ,41 + ,8 + ,92630 + ,40 + ,27 + ,18513 + ,1205 + ,387 + ,60 + ,6782 + ,31 + ,17 + ,121848 + ,39 + ,37 + ,30594 + ,1255 + ,339 + ,45 + ,5.401 + ,41 + ,18 + ,52915 + ,14 + ,20 + ,24006 + ,613 + ,181 + ,54 + ,6521 + ,26 + ,16 + ,81872 + ,45 + ,32 + ,27913 + ,721 + ,245 + ,25 + ,10.856 + ,23 + ,23 + ,58981 + ,36 + ,0 + ,42744 + ,1109 + ,384 + ,38 + ,2.154 + ,14 + ,22 + ,53515 + ,28 + ,5 + ,12934 + ,740 + ,212 + ,52 + ,6.117 + ,16 + ,13 + ,60812 + ,44 + ,26 + ,22574 + ,1126 + ,399 + ,67 + ,5238 + ,25 + ,13 + ,56375 + ,30 + ,10 + ,41385 + ,728 + ,229 + ,74 + ,4.820 + ,21 + ,16 + ,65490 + ,22 + ,27 + ,18653 + ,689 + ,224 + ,38 + ,5.615 + ,32 + ,16 + ,80949 + ,17 + ,11 + ,18472 + ,592 + ,203 + ,30 + ,4272 + ,9 + ,20 + ,76302 + ,31 + ,29 + ,30976 + ,995 + ,333 + ,26 + ,8.702 + ,35 + ,22 + ,104011 + ,55 + ,25 + ,63339 + ,1613 + ,384 + ,67 + ,15.340 + ,42 + ,17 + ,98104 + ,54 + ,55 + ,25568 + ,2048 + ,636 + ,132 + ,8.030 + ,68 + ,18 + ,67989 + ,21 + ,23 + ,33747 + ,705 + ,185 + ,42 + ,9526 + ,32 + ,17 + ,30989 + ,14 + ,5 + ,4154 + ,301 + ,93 + ,35 + ,1278 + ,6 + ,12 + ,135458 + ,81 + ,43 + ,19474 + ,1803 + ,581 + ,118 + ,4.236 + ,68 + ,7 + ,73504 + ,35 + ,23 + ,35130 + ,799 + ,248 + ,68 + ,3023 + ,33 + ,17 + ,63123 + ,43 + ,34 + ,39067 + ,861 + ,304 + ,43 + ,7.196 + ,84 + ,14 + ,61254 + ,46 + ,36 + ,13310 + ,1186 + ,344 + ,76 + ,3394 + ,46 + ,23 + ,74914 + ,30 + ,35 + ,65892 + ,1451 + ,407 + ,64 + ,6.371 + ,30 + ,17 + ,31774 + ,23 + ,0 + ,4143 + ,628 + ,170 + ,48 + ,1.574 + ,0 + ,14 + ,81437 + ,38 + ,37 + ,28579 + ,1161 + ,312 + ,64 + ,9.620 + ,36 + ,15 + ,87186 + ,54 + ,28 + ,51776 + ,1463 + ,507 + ,56 + ,6978 + ,47 + ,17 + ,50090 + ,20 + ,16 + ,21152 + ,742 + ,224 + ,71 + ,4911 + ,20 + ,21 + ,65745 + ,53 + ,26 + ,38084 + ,979 + ,340 + ,75 + ,8.645 + ,50 + ,18 + ,56653 + ,45 + ,38 + ,27717 + ,675 + ,168 + ,39 + ,8.987 + ,30 + ,18 + ,158399 + ,39 + ,23 + ,32928 + ,1241 + ,443 + ,42 + ,5.544 + ,30 + ,17 + ,46455 + ,20 + ,22 + ,11342 + ,676 + ,204 + ,39 + ,3083 + ,34 + ,17 + ,73624 + ,24 + ,30 + ,19499 + ,1049 + ,367 + ,93 + ,6.909 + ,33 + ,16 + ,38395 + ,31 + ,16 + ,16380 + ,620 + ,210 + ,38 + ,3189 + ,34 + ,15 + ,91899 + ,35 + ,18 + ,36874 + ,1081 + ,335 + ,60 + ,6.745 + ,37 + ,21 + ,139526 + ,151 + ,28 + ,48259 + ,1688 + ,364 + ,71 + ,16.724 + ,83 + ,16 + ,52164 + ,52 + ,32 + ,16734 + ,736 + ,178 + ,52 + ,4850 + ,32 + ,14 + ,51567 + ,30 + ,21 + ,28207 + ,617 + ,206 + ,27 + ,7.025 + ,30 + ,15 + ,70551 + ,31 + ,23 + ,30143 + ,812 + ,279 + ,59 + ,6047 + ,43 + ,17 + ,84856 + ,29 + ,29 + ,41369 + ,1051 + ,387 + ,40 + ,7377 + ,41 + ,15 + ,102538 + ,57 + ,50 + ,45833 + ,1656 + ,490 + ,79 + ,9.078 + ,51 + ,15 + ,86678 + ,40 + ,12 + ,29156 + ,705 + ,238 + ,44 + ,4.605 + ,19 + ,10 + ,85709 + ,44 + ,21 + ,35944 + ,945 + ,343 + ,65 + ,3238 + ,37 + ,6 + ,34662 + ,25 + ,18 + ,36278 + ,554 + ,232 + ,10 + ,8100 + ,33 + ,22 + ,150580 + ,77 + ,27 + ,45588 + ,1597 + ,530 + ,124 + ,9.653 + ,41 + ,21 + ,99611 + ,35 + ,41 + ,45097 + ,982 + ,291 + ,81 + ,8.914 + ,54 + ,1 + ,19349 + ,11 + ,13 + ,3895 + ,222 + ,67 + ,15 + ,786 + ,14 + ,18 + ,99373 + ,63 + ,12 + ,28394 + ,1212 + ,397 + ,92 + ,6.700 + ,25 + ,17 + ,86230 + ,44 + ,21 + ,18632 + ,1143 + ,467 + ,42 + ,5.788 + ,25 + ,4 + ,30837 + ,19 + ,8 + ,2325 + ,435 + ,178 + ,10 + ,593 + ,8 + ,10 + ,31706 + ,13 + ,26 + ,25139 + ,532 + ,175 + ,24 + ,4.506 + ,26 + ,16 + ,89806 + ,42 + ,27 + ,27975 + ,882 + ,299 + ,64 + ,6.382 + ,20 + ,16 + ,62088 + ,38 + ,13 + ,14483 + ,608 + ,154 + ,45 + ,5621 + ,11 + ,9 + ,40151 + ,29 + ,16 + ,13127 + ,459 + ,106 + ,22 + ,3997 + ,14 + ,16 + ,27634 + ,20 + ,2 + ,5839 + ,578 + ,189 + ,56 + ,520 + ,3 + ,17 + ,76990 + ,27 + ,42 + ,24069 + ,826 + ,194 + ,94 + ,8891 + ,40 + ,7 + ,37460 + ,20 + ,5 + ,3738 + ,509 + ,135 + ,19 + ,999 + ,5 + ,15 + ,54157 + ,19 + ,37 + ,18625 + ,717 + ,201 + ,35 + ,7067 + ,38 + ,14 + ,49862 + ,37 + ,17 + ,36341 + ,637 + ,207 + ,32 + ,4639 + ,32 + ,14 + ,84337 + ,26 + ,38 + ,24548 + ,857 + ,280 + ,35 + ,5654 + ,41 + ,18 + ,64175 + ,42 + ,37 + ,21792 + ,830 + ,260 + ,48 + ,6.928 + ,46 + ,12 + ,59382 + ,49 + ,29 + ,26263 + ,652 + ,227 + ,49 + ,1.514 + ,47 + ,16 + ,119308 + ,30 + ,32 + ,23686 + ,707 + ,239 + ,48 + ,9.238 + ,37 + ,21 + ,76702 + ,49 + ,35 + ,49303 + ,954 + ,333 + ,62 + ,8.204 + ,51 + ,19 + ,103425 + ,67 + ,17 + ,25659 + ,1461 + ,428 + ,96 + ,5926 + ,49 + ,16 + ,70344 + ,28 + ,20 + ,28904 + ,672 + ,230 + ,45 + ,5785 + ,21 + ,1 + ,43410 + ,19 + ,7 + ,2781 + ,778 + ,292 + ,63 + ,4 + ,1 + ,16 + ,104838 + ,49 + ,46 + ,29236 + ,1141 + ,350 + ,71 + ,5930 + ,44 + ,10 + ,62215 + ,27 + ,24 + ,19546 + ,680 + ,186 + ,26 + ,3710 + ,26 + ,19 + ,69304 + ,30 + ,40 + ,22818 + ,1090 + ,326 + ,48 + ,705 + ,21 + ,12 + ,53117 + ,22 + ,3 + ,32689 + ,616 + ,155 + ,29 + ,443 + ,4 + ,2 + ,19764 + ,12 + ,10 + ,5752 + ,285 + ,75 + ,19 + ,2.416 + ,10 + ,14 + ,86680 + ,31 + ,37 + ,22197 + ,1145 + ,361 + ,45 + ,7747 + ,43 + ,17 + ,84105 + ,20 + ,17 + ,20055 + ,733 + ,261 + ,45 + ,5.432 + ,34 + ,19 + ,77945 + ,20 + ,28 + ,25272 + ,888 + ,299 + ,67 + ,4913 + ,32 + ,14 + ,89113 + ,39 + ,19 + ,82206 + ,849 + ,300 + ,30 + ,2650 + ,20 + ,11 + ,91005 + ,29 + ,29 + ,32073 + ,1182 + ,450 + ,36 + ,2370 + ,34 + ,4 + ,40248 + ,16 + ,8 + ,5444 + ,528 + ,183 + ,34 + ,775 + ,6 + ,16 + ,64187 + ,27 + ,10 + ,20154 + ,642 + ,238 + ,36 + ,5576 + ,12 + ,20 + ,50857 + ,21 + ,15 + ,36944 + ,947 + ,165 + ,34 + ,1352 + ,24 + ,12 + ,56613 + ,19 + ,15 + ,8019 + ,819 + ,234 + ,37 + ,3080 + ,16 + ,15 + ,62792 + ,35 + ,28 + ,30884 + ,757 + ,176 + ,46 + ,10205 + ,72 + ,16 + ,72535 + ,14 + ,17 + ,19540 + ,894 + ,329 + ,44 + ,6095 + ,27) + ,dim=c(10 + ,289) + ,dimnames=list(c('compendiums_reviewed' + ,'time_in_rfc' + ,'logins' + ,'blogged_computations' + ,'totsize' + ,'pageviews' + ,'compendium_views_info' + ,'compendium_views_pr' + ,'totrevisions' + ,'tothyperlinks') + ,1:289)) > y <- array(NA,dim=c(10,289),dimnames=list(c('compendiums_reviewed','time_in_rfc','logins','blogged_computations','totsize','pageviews','compendium_views_info','compendium_views_pr','totrevisions','tothyperlinks'),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 = 'none' > par2 = 'none' > par1 = '1' > 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) Warning message: NAs introduced by coercion > 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] "compendiums_reviewed" > x[,par1] [1] 30 28 38 30 22 26 25 18 11 26 25 38 44 30 40 34 47 30 31 23 36 36 30 25 39 [26] 34 31 31 33 25 33 35 42 43 30 33 13 32 36 0 28 14 17 32 30 35 20 28 28 39 [51] 34 26 39 39 33 28 4 39 18 14 29 44 21 16 28 35 28 38 23 36 32 29 25 27 36 [76] 28 23 40 23 40 28 34 33 28 34 30 33 22 38 26 35 8 24 29 20 29 45 37 33 33 [101] 25 32 29 28 28 31 52 21 24 41 33 32 19 20 31 31 32 18 23 17 20 12 17 30 31 [126] 10 13 22 42 1 9 32 11 25 36 31 0 24 13 8 13 19 18 33 40 22 38 24 8 35 [151] 43 43 14 41 38 45 31 13 28 31 40 30 16 37 30 35 32 27 20 18 31 31 21 39 41 [176] 13 32 18 39 14 7 17 0 30 37 0 5 1 16 32 24 17 11 24 22 12 19 13 17 15 [201] 16 24 15 17 18 20 16 16 18 22 8 17 18 16 23 22 13 13 16 16 20 22 17 18 17 [226] 12 7 17 14 23 17 14 15 17 21 18 18 17 17 16 15 21 16 14 15 17 15 15 10 6 [251] 22 21 1 18 17 4 10 16 16 9 16 17 7 15 14 14 18 12 16 21 19 16 1 16 10 [276] 19 12 2 14 17 19 14 11 4 16 20 12 15 16 > 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]) 0 1 2 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 4 4 1 3 1 1 3 4 2 4 4 6 9 11 9 19 19 14 6 8 7 9 7 7 7 4 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 47 52 2 13 5 12 12 10 10 5 6 6 3 6 7 5 3 2 3 2 2 1 1 > colnames(x) [1] "compendiums_reviewed" "time_in_rfc" "logins" [4] "blogged_computations" "totsize" "pageviews" [7] "compendium_views_info" "compendium_views_pr" "totrevisions" [10] "tothyperlinks" > colnames(x)[par1] [1] "compendiums_reviewed" > x[,par1] [1] 30 28 38 30 22 26 25 18 11 26 25 38 44 30 40 34 47 30 31 23 36 36 30 25 39 [26] 34 31 31 33 25 33 35 42 43 30 33 13 32 36 0 28 14 17 32 30 35 20 28 28 39 [51] 34 26 39 39 33 28 4 39 18 14 29 44 21 16 28 35 28 38 23 36 32 29 25 27 36 [76] 28 23 40 23 40 28 34 33 28 34 30 33 22 38 26 35 8 24 29 20 29 45 37 33 33 [101] 25 32 29 28 28 31 52 21 24 41 33 32 19 20 31 31 32 18 23 17 20 12 17 30 31 [126] 10 13 22 42 1 9 32 11 25 36 31 0 24 13 8 13 19 18 33 40 22 38 24 8 35 [151] 43 43 14 41 38 45 31 13 28 31 40 30 16 37 30 35 32 27 20 18 31 31 21 39 41 [176] 13 32 18 39 14 7 17 0 30 37 0 5 1 16 32 24 17 11 24 22 12 19 13 17 15 [201] 16 24 15 17 18 20 16 16 18 22 8 17 18 16 23 22 13 13 16 16 20 22 17 18 17 [226] 12 7 17 14 23 17 14 15 17 21 18 18 17 17 16 15 21 16 14 15 17 15 15 10 6 [251] 22 21 1 18 17 4 10 16 16 9 16 17 7 15 14 14 18 12 16 21 19 16 1 16 10 [276] 19 12 2 14 17 19 14 11 4 16 20 12 15 16 > 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/1nu8q1324133728.tab") + } + } > m Conditional inference tree with 8 terminal nodes Response: compendiums_reviewed Inputs: time_in_rfc, logins, blogged_computations, totsize, pageviews, compendium_views_info, compendium_views_pr, totrevisions, tothyperlinks Number of observations: 289 1) blogged_computations <= 54; criterion = 1, statistic = 167.343 2) totsize <= 56622; criterion = 1, statistic = 66.676 3) totsize <= 8773; criterion = 1, statistic = 42.98 4)* weights = 19 3) totsize > 8773 5) time_in_rfc <= 108446; criterion = 0.997, statistic = 13.063 6) compendium_views_pr <= 35; criterion = 0.956, statistic = 7.872 7)* weights = 24 6) compendium_views_pr > 35 8)* weights = 77 5) time_in_rfc > 108446 9)* weights = 29 2) totsize > 56622 10) compendium_views_pr <= 46; criterion = 0.959, statistic = 8.029 11)* weights = 9 10) compendium_views_pr > 46 12)* weights = 16 1) blogged_computations > 54 13) blogged_computations <= 74; criterion = 1, statistic = 19.205 14)* weights = 35 13) blogged_computations > 74 15)* weights = 80 > postscript(file="/var/wessaorg/rcomp/tmp/2noso1324133728.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/3vzel1324133728.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 30 34.362500 -4.36250000 2 28 28.885714 -0.88571429 3 38 28.885714 9.11428571 4 30 34.362500 -4.36250000 5 22 19.931034 2.06896552 6 26 22.666667 3.33333333 7 25 34.362500 -9.36250000 8 18 16.636364 1.36363636 9 11 16.636364 -5.63636364 10 26 30.500000 -4.50000000 11 25 28.885714 -3.88571429 12 38 34.362500 3.63750000 13 44 34.362500 9.63750000 14 30 30.500000 -0.50000000 15 40 34.362500 5.63750000 16 34 28.885714 5.11428571 17 47 34.362500 12.63750000 18 30 34.362500 -4.36250000 19 31 34.362500 -3.36250000 20 23 34.362500 -11.36250000 21 36 34.362500 1.63750000 22 36 34.362500 1.63750000 23 30 34.362500 -4.36250000 24 25 30.500000 -5.50000000 25 39 34.362500 4.63750000 26 34 34.362500 -0.36250000 27 31 28.885714 2.11428571 28 31 34.362500 -3.36250000 29 33 34.362500 -1.36250000 30 25 19.931034 5.06896552 31 33 34.362500 -1.36250000 32 35 34.362500 0.63750000 33 42 30.500000 11.50000000 34 43 34.362500 8.63750000 35 30 22.666667 7.33333333 36 33 34.362500 -1.36250000 37 13 16.636364 -3.63636364 38 32 28.885714 3.11428571 39 36 34.362500 1.63750000 40 0 5.263158 -5.26315789 41 28 28.885714 -0.88571429 42 14 13.458333 0.54166667 43 17 16.636364 0.36363636 44 32 34.362500 -2.36250000 45 30 34.362500 -4.36250000 46 35 28.885714 6.11428571 47 20 28.885714 -8.88571429 48 28 28.885714 -0.88571429 49 28 19.931034 8.06896552 50 39 30.500000 8.50000000 51 34 34.362500 -0.36250000 52 26 16.636364 9.36363636 53 39 34.362500 4.63750000 54 39 34.362500 4.63750000 55 33 34.362500 -1.36250000 56 28 34.362500 -6.36250000 57 4 5.263158 -1.26315789 58 39 34.362500 4.63750000 59 18 13.458333 4.54166667 60 14 16.636364 -2.63636364 61 29 34.362500 -5.36250000 62 44 34.362500 9.63750000 63 21 34.362500 -13.36250000 64 16 16.636364 -0.63636364 65 28 34.362500 -6.36250000 66 35 34.362500 0.63750000 67 28 34.362500 -6.36250000 68 38 34.362500 3.63750000 69 23 34.362500 -11.36250000 70 36 28.885714 7.11428571 71 32 28.885714 3.11428571 72 29 30.500000 -1.50000000 73 25 28.885714 -3.88571429 74 27 28.885714 -1.88571429 75 36 34.362500 1.63750000 76 28 34.362500 -6.36250000 77 23 34.362500 -11.36250000 78 40 34.362500 5.63750000 79 23 16.636364 6.36363636 80 40 28.885714 11.11428571 81 28 28.885714 -0.88571429 82 34 34.362500 -0.36250000 83 33 34.362500 -1.36250000 84 28 34.362500 -6.36250000 85 34 28.885714 5.11428571 86 30 34.362500 -4.36250000 87 33 28.885714 4.11428571 88 22 19.931034 2.06896552 89 38 34.362500 3.63750000 90 26 19.931034 6.06896552 91 35 34.362500 0.63750000 92 8 13.458333 -5.45833333 93 24 19.931034 4.06896552 94 29 28.885714 0.11428571 95 20 19.931034 0.06896552 96 29 28.885714 0.11428571 97 45 34.362500 10.63750000 98 37 34.362500 2.63750000 99 33 28.885714 4.11428571 100 33 34.362500 -1.36250000 101 25 19.931034 5.06896552 102 32 34.362500 -2.36250000 103 29 28.885714 0.11428571 104 28 28.885714 -0.88571429 105 28 34.362500 -6.36250000 106 31 34.362500 -3.36250000 107 52 34.362500 17.63750000 108 21 16.636364 4.36363636 109 24 34.362500 -10.36250000 110 41 34.362500 6.63750000 111 33 34.362500 -1.36250000 112 32 19.931034 12.06896552 113 19 16.636364 2.36363636 114 20 28.885714 -8.88571429 115 31 34.362500 -3.36250000 116 31 34.362500 -3.36250000 117 32 34.362500 -2.36250000 118 18 19.931034 -1.93103448 119 23 22.666667 0.33333333 120 17 16.636364 0.36363636 121 20 28.885714 -8.88571429 122 12 16.636364 -4.63636364 123 17 19.931034 -2.93103448 124 30 34.362500 -4.36250000 125 31 30.500000 0.50000000 126 10 13.458333 -3.45833333 127 13 16.636364 -3.63636364 128 22 28.885714 -6.88571429 129 42 34.362500 7.63750000 130 1 5.263158 -4.26315789 131 9 16.636364 -7.63636364 132 32 34.362500 -2.36250000 133 11 13.458333 -2.45833333 134 25 16.636364 8.36363636 135 36 28.885714 7.11428571 136 31 34.362500 -3.36250000 137 0 5.263158 -5.26315789 138 24 28.885714 -4.88571429 139 13 16.636364 -3.63636364 140 8 13.458333 -5.45833333 141 13 13.458333 -0.45833333 142 19 28.885714 -9.88571429 143 18 16.636364 1.36363636 144 33 34.362500 -1.36250000 145 40 28.885714 11.11428571 146 22 16.636364 5.36363636 147 38 30.500000 7.50000000 148 24 22.666667 1.33333333 149 8 5.263158 2.73684211 150 35 30.500000 4.50000000 151 43 34.362500 8.63750000 152 43 30.500000 12.50000000 153 14 13.458333 0.54166667 154 41 34.362500 6.63750000 155 38 34.362500 3.63750000 156 45 34.362500 10.63750000 157 31 28.885714 2.11428571 158 13 16.636364 -3.63636364 159 28 22.666667 5.33333333 160 31 28.885714 2.11428571 161 40 34.362500 5.63750000 162 30 30.500000 -0.50000000 163 16 19.931034 -3.93103448 164 37 34.362500 2.63750000 165 30 19.931034 10.06896552 166 35 34.362500 0.63750000 167 32 28.885714 3.11428571 168 27 30.500000 -3.50000000 169 20 19.931034 0.06896552 170 18 30.500000 -12.50000000 171 31 34.362500 -3.36250000 172 31 34.362500 -3.36250000 173 21 19.931034 1.06896552 174 39 34.362500 4.63750000 175 41 34.362500 6.63750000 176 13 19.931034 -6.93103448 177 32 34.362500 -2.36250000 178 18 16.636364 1.36363636 179 39 34.362500 4.63750000 180 14 19.931034 -5.93103448 181 7 16.636364 -9.63636364 182 17 28.885714 -11.88571429 183 0 5.263158 -5.26315789 184 30 30.500000 -0.50000000 185 37 34.362500 2.63750000 186 0 5.263158 -5.26315789 187 5 5.263158 -0.26315789 188 1 5.263158 -4.26315789 189 16 19.931034 -3.93103448 190 32 16.636364 15.36363636 191 24 22.666667 1.33333333 192 17 16.636364 0.36363636 193 11 16.636364 -5.63636364 194 24 16.636364 7.36363636 195 22 22.666667 -0.66666667 196 12 19.931034 -7.93103448 197 19 16.636364 2.36363636 198 13 22.666667 -9.66666667 199 17 16.636364 0.36363636 200 15 13.458333 1.54166667 201 16 13.458333 2.54166667 202 24 16.636364 7.36363636 203 15 19.931034 -4.93103448 204 17 16.636364 0.36363636 205 18 16.636364 1.36363636 206 20 13.458333 6.54166667 207 16 16.636364 -0.63636364 208 16 19.931034 -3.93103448 209 18 19.931034 -1.93103448 210 22 19.931034 2.06896552 211 8 16.636364 -8.63636364 212 17 19.931034 -2.93103448 213 18 16.636364 1.36363636 214 16 13.458333 2.54166667 215 23 16.636364 6.36363636 216 22 16.636364 5.36363636 217 13 16.636364 -3.63636364 218 13 16.636364 -3.63636364 219 16 16.636364 -0.63636364 220 16 13.458333 2.54166667 221 20 13.458333 6.54166667 222 22 30.500000 -8.50000000 223 17 28.885714 -11.88571429 224 18 16.636364 1.36363636 225 17 5.263158 11.73684211 226 12 19.931034 -7.93103448 227 7 16.636364 -9.63636364 228 17 16.636364 0.36363636 229 14 16.636364 -2.63636364 230 23 30.500000 -7.50000000 231 17 5.263158 11.73684211 232 14 16.636364 -2.63636364 233 15 16.636364 -1.63636364 234 17 16.636364 0.36363636 235 21 16.636364 4.36363636 236 18 16.636364 1.36363636 237 18 19.931034 -1.93103448 238 17 16.636364 0.36363636 239 17 16.636364 0.36363636 240 16 16.636364 -0.63636364 241 15 16.636364 -1.63636364 242 21 19.931034 1.06896552 243 16 16.636364 -0.63636364 244 14 13.458333 0.54166667 245 15 16.636364 -1.63636364 246 17 16.636364 0.36363636 247 15 16.636364 -1.63636364 248 15 16.636364 -1.63636364 249 10 16.636364 -6.63636364 250 6 13.458333 -7.45833333 251 22 19.931034 2.06896552 252 21 16.636364 4.36363636 253 1 5.263158 -4.26315789 254 18 16.636364 1.36363636 255 17 16.636364 0.36363636 256 4 5.263158 -1.26315789 257 10 13.458333 -3.45833333 258 16 16.636364 -0.63636364 259 16 16.636364 -0.63636364 260 9 13.458333 -4.45833333 261 16 5.263158 10.73684211 262 17 16.636364 0.36363636 263 7 5.263158 1.73684211 264 15 13.458333 1.54166667 265 14 13.458333 0.54166667 266 14 13.458333 0.54166667 267 18 16.636364 1.36363636 268 12 16.636364 -4.63636364 269 16 19.931034 -3.93103448 270 21 16.636364 4.36363636 271 19 16.636364 2.36363636 272 16 16.636364 -0.63636364 273 1 5.263158 -4.26315789 274 16 16.636364 -0.63636364 275 10 13.458333 -3.45833333 276 19 16.636364 2.36363636 277 12 13.458333 -1.45833333 278 2 5.263158 -3.26315789 279 14 16.636364 -2.63636364 280 17 16.636364 0.36363636 281 19 16.636364 2.36363636 282 14 22.666667 -8.66666667 283 11 16.636364 -5.63636364 284 4 5.263158 -1.26315789 285 16 16.636364 -0.63636364 286 20 13.458333 6.54166667 287 12 5.263158 6.73684211 288 15 16.636364 -1.63636364 289 16 16.636364 -0.63636364 > 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/491gq1324133728.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/5wab61324133728.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/6bsdq1324133728.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/7r2lx1324133728.tab") + } > > try(system("convert tmp/2noso1324133728.ps tmp/2noso1324133728.png",intern=TRUE)) character(0) > try(system("convert tmp/3vzel1324133728.ps tmp/3vzel1324133728.png",intern=TRUE)) character(0) > try(system("convert tmp/491gq1324133728.ps tmp/491gq1324133728.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 6.325 0.267 6.595