R version 2.12.0 (2010-10-15) Copyright (C) 2010 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(1536 + ,127476 + ,78 + ,490 + ,107 + ,0 + ,20 + ,17 + ,66 + ,59 + ,18158 + ,5636 + ,22622 + ,30 + ,28 + ,1134 + ,130358 + ,46 + ,329 + ,68 + ,1 + ,38 + ,17 + ,68 + ,50 + ,30461 + ,9079 + ,73570 + ,42 + ,39 + ,192 + ,7215 + ,18 + ,72 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1423 + ,603 + ,1929 + ,0 + ,0 + ,2032 + ,112861 + ,84 + ,584 + ,146 + ,0 + ,49 + ,22 + ,68 + ,51 + ,25629 + ,8874 + ,36294 + ,54 + ,54 + ,3231 + ,210171 + ,125 + ,1077 + ,124 + ,0 + ,74 + ,30 + ,120 + ,112 + ,48758 + ,17988 + ,62378 + ,86 + ,80 + ,5777 + ,393802 + ,215 + ,1578 + ,267 + ,1 + ,104 + ,31 + ,120 + ,118 + ,129230 + ,21325 + ,167760 + ,157 + ,144 + ,1322 + ,117604 + ,50 + ,442 + ,83 + ,1 + ,37 + ,19 + ,72 + ,59 + ,27376 + ,8325 + ,52443 + ,36 + ,36 + ,1181 + ,126029 + ,48 + ,319 + ,48 + ,0 + ,53 + ,25 + ,96 + ,90 + ,26706 + ,7117 + ,57283 + ,48 + ,48 + ,1462 + ,99729 + ,37 + ,406 + ,87 + ,0 + ,42 + ,30 + ,109 + ,50 + ,26505 + ,7996 + ,36614 + ,45 + ,42 + ,2568 + ,256310 + ,86 + ,818 + ,129 + ,1 + ,62 + ,26 + ,104 + ,79 + ,49801 + ,14218 + ,93268 + ,77 + ,71 + ,1810 + ,113066 + ,69 + ,568 + ,146 + ,2 + ,50 + ,20 + ,54 + ,49 + ,46580 + ,6321 + ,35439 + ,49 + ,49 + ,1789 + ,156212 + ,59 + ,551 + ,94 + ,0 + ,65 + ,25 + ,98 + ,74 + ,48352 + ,19690 + ,72405 + ,77 + ,74 + ,1334 + ,69952 + ,85 + ,494 + ,57 + ,0 + ,28 + ,15 + ,49 + ,32 + ,13899 + ,5659 + ,24044 + ,28 + ,27 + ,2415 + ,152673 + ,84 + ,818 + ,240 + ,4 + ,48 + ,22 + ,88 + ,82 + ,39342 + ,11370 + ,55909 + ,84 + ,83 + ,1156 + ,125841 + ,44 + ,331 + ,40 + ,4 + ,42 + ,12 + ,45 + ,43 + ,27465 + ,4778 + ,44689 + ,31 + ,31 + ,1374 + ,125769 + ,67 + ,419 + ,81 + ,3 + ,47 + ,19 + ,74 + ,65 + ,55211 + ,5954 + ,49319 + ,28 + ,28 + ,1504 + ,123467 + ,50 + ,364 + ,85 + ,0 + ,71 + ,28 + ,112 + ,111 + ,74098 + ,22924 + ,62075 + ,99 + ,98 + ,999 + ,56232 + ,47 + ,284 + ,62 + ,5 + ,0 + ,12 + ,45 + ,36 + ,13497 + ,70 + ,2341 + ,2 + ,2 + ,2190 + ,108244 + ,77 + ,667 + ,126 + ,0 + ,50 + ,28 + ,110 + ,89 + ,38338 + ,14369 + ,40551 + ,41 + ,43 + ,633 + ,22762 + ,20 + ,188 + ,44 + ,0 + ,12 + ,13 + ,39 + ,28 + ,52505 + ,3706 + ,11621 + ,25 + ,24 + ,838 + ,48554 + ,49 + ,286 + ,37 + ,0 + ,16 + ,14 + ,55 + ,35 + ,10663 + ,3147 + ,18741 + ,16 + ,16 + ,2167 + ,178697 + ,81 + ,633 + ,94 + ,0 + ,76 + ,27 + ,102 + ,78 + ,74484 + ,16801 + ,84202 + ,96 + ,95 + ,1452 + ,139115 + ,58 + ,514 + ,127 + ,0 + ,29 + ,25 + ,96 + ,67 + ,28895 + ,2162 + ,15334 + ,23 + ,22 + ,1790 + ,93773 + ,45 + ,532 + ,159 + ,1 + ,38 + ,30 + ,86 + ,61 + ,32827 + ,4721 + ,28024 + ,33 + ,33 + ,1718 + ,132796 + ,76 + ,540 + ,41 + ,1 + ,50 + ,20 + ,74 + ,55 + ,36188 + ,5290 + ,53306 + ,46 + ,45 + ,1179 + ,113933 + ,22 + ,428 + ,153 + ,0 + ,33 + ,17 + ,64 + ,49 + ,28173 + ,6446 + ,37918 + ,59 + ,59 + ,1688 + ,144781 + ,138 + ,539 + ,86 + ,0 + ,45 + ,22 + ,82 + ,77 + ,54926 + ,14711 + ,54819 + ,72 + ,66 + ,1101 + ,140711 + ,75 + ,266 + ,55 + ,0 + ,59 + ,28 + ,100 + ,71 + ,38900 + ,13311 + ,89058 + ,72 + ,70 + ,2259 + ,283337 + ,102 + ,745 + ,73 + ,0 + ,49 + ,25 + ,95 + ,82 + ,88530 + ,13577 + ,103354 + ,62 + ,56 + ,1768 + ,158146 + ,36 + ,733 + ,79 + ,0 + ,40 + ,16 + ,63 + ,53 + ,35482 + ,14634 + ,70239 + ,55 + ,55 + ,1300 + ,123344 + ,39 + ,394 + ,71 + ,0 + ,40 + ,23 + ,87 + ,71 + ,26730 + ,6931 + ,33045 + ,27 + ,27 + ,1432 + ,157640 + ,38 + ,482 + ,111 + ,2 + ,51 + ,20 + ,65 + ,58 + ,29806 + ,9992 + ,63852 + ,41 + ,37 + ,1791 + ,91279 + ,88 + ,567 + ,71 + ,4 + ,41 + ,11 + ,43 + ,25 + ,41799 + ,6185 + ,30905 + ,51 + ,48 + ,2475 + ,189374 + ,102 + ,746 + ,243 + ,0 + ,73 + ,20 + ,80 + ,59 + ,54289 + ,3445 + ,24242 + ,26 + ,26 + ,1930 + ,167915 + ,42 + ,626 + ,66 + ,1 + ,43 + ,21 + ,84 + ,77 + ,36805 + ,12327 + ,78907 + ,65 + ,64 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1782 + ,175403 + ,54 + ,835 + ,58 + ,0 + ,46 + ,27 + ,105 + ,75 + ,33146 + ,9898 + ,36005 + ,28 + ,21 + ,1505 + ,92342 + ,46 + ,464 + ,131 + ,3 + ,44 + ,14 + ,51 + ,39 + ,23333 + ,8022 + ,31972 + ,44 + ,44 + ,1820 + ,100023 + ,41 + ,418 + ,258 + ,9 + ,31 + ,29 + ,98 + ,83 + ,47686 + ,10765 + ,35853 + ,36 + ,36 + ,1648 + ,178277 + ,49 + ,607 + ,56 + ,0 + ,71 + ,31 + ,124 + ,123 + ,77783 + ,22717 + ,115301 + ,100 + ,89 + ,1668 + ,145062 + ,56 + ,539 + ,90 + ,2 + ,61 + ,19 + ,75 + ,67 + ,36042 + ,10090 + ,47689 + ,104 + ,101 + ,1366 + ,110980 + ,47 + ,519 + ,57 + ,0 + ,28 + ,30 + ,120 + ,105 + ,34541 + ,12385 + ,34223 + ,35 + ,31 + ,864 + ,86039 + ,25 + ,309 + ,35 + ,2 + ,21 + ,23 + ,84 + ,76 + ,75620 + ,8513 + ,43431 + ,69 + ,65 + ,1602 + ,119514 + ,62 + ,609 + ,50 + ,1 + ,42 + ,20 + ,78 + ,54 + ,60610 + ,5508 + ,52220 + ,73 + ,71 + ,1023 + ,95535 + ,41 + ,321 + ,46 + ,2 + ,44 + ,22 + ,87 + ,82 + ,55041 + ,9628 + ,33863 + ,106 + ,102 + ,963 + ,109894 + ,73 + ,246 + ,32 + ,2 + ,34 + ,19 + ,70 + ,57 + ,32087 + ,11872 + ,46879 + ,53 + ,53 + ,629 + ,61554 + ,26 + ,180 + ,45 + ,1 + ,15 + ,32 + ,97 + ,57 + ,16356 + ,4186 + ,23228 + ,43 + ,41 + ,1568 + ,156520 + ,77 + ,544 + ,96 + ,0 + ,46 + ,18 + ,72 + ,72 + ,40161 + ,10877 + ,42827 + ,49 + ,46 + ,1715 + ,159121 + ,75 + ,544 + ,104 + ,1 + ,43 + ,26 + ,104 + ,94 + ,55459 + ,17066 + ,65765 + ,38 + ,37 + ,2093 + ,129362 + ,51 + ,758 + ,150 + ,4 + ,47 + ,25 + ,93 + ,72 + ,36679 + ,9175 + ,38167 + ,51 + ,51 + ,658 + ,48188 + ,28 + ,205 + ,37 + ,0 + ,12 + ,22 + ,82 + ,39 + ,22346 + ,2102 + ,14812 + ,14 + ,14 + ,1199 + ,91198 + ,54 + ,309 + ,49 + ,0 + ,42 + ,19 + ,73 + ,60 + ,27377 + ,10807 + ,32615 + ,40 + ,40 + ,2059 + ,229864 + ,64 + ,709 + ,83 + ,0 + ,56 + ,24 + ,87 + ,84 + ,50273 + ,13662 + ,82188 + ,79 + ,77 + ,1592 + ,180317 + ,67 + ,542 + ,67 + ,0 + ,41 + ,26 + ,95 + ,69 + ,32104 + ,9224 + ,51763 + ,52 + ,51 + ,1447 + ,150640 + ,48 + ,526 + ,39 + ,1 + ,48 + ,27 + ,105 + ,102 + ,27016 + ,9001 + ,59325 + ,44 + ,43 + ,1342 + ,104416 + ,44 + ,418 + ,68 + ,5 + ,30 + ,10 + ,37 + ,28 + ,19715 + ,7204 + ,48976 + ,34 + ,33 + ,1527 + ,159645 + ,55 + ,409 + ,58 + ,0 + ,44 + ,26 + ,96 + ,65 + ,33629 + ,6572 + ,43384 + ,47 + ,47 + ,670 + ,60368 + ,17 + ,189 + ,59 + ,0 + ,25 + ,21 + ,80 + ,59 + ,27084 + ,7509 + ,26692 + ,32 + ,31 + ,859 + ,100056 + ,55 + ,293 + ,30 + ,0 + ,42 + ,21 + ,83 + ,80 + ,32352 + ,12920 + ,53279 + ,31 + ,31 + ,2329 + ,137214 + ,73 + ,781 + ,54 + ,10 + ,28 + ,34 + ,124 + ,79 + ,51845 + ,5438 + ,20652 + ,40 + ,40 + ,1326 + ,99630 + ,47 + ,383 + ,65 + ,6 + ,33 + ,29 + ,116 + ,107 + ,26591 + ,11489 + ,38338 + ,42 + ,42 + ,1567 + ,84557 + ,62 + ,572 + ,81 + ,0 + ,32 + ,18 + ,72 + ,57 + ,29677 + ,6661 + ,36735 + ,34 + ,35 + ,1081 + ,91199 + ,45 + ,308 + ,81 + ,11 + ,28 + ,16 + ,55 + ,44 + ,54237 + ,7941 + ,42764 + ,40 + ,40 + ,897 + ,83419 + ,29 + ,288 + ,45 + ,3 + ,31 + ,23 + ,86 + ,59 + ,20284 + ,6173 + ,44331 + ,35 + ,30 + ,855 + ,101723 + ,25 + ,285 + ,52 + ,0 + ,13 + ,22 + ,85 + ,80 + ,22741 + ,5562 + ,41354 + ,11 + ,11 + ,1229 + ,94982 + ,37 + ,391 + ,36 + ,0 + ,38 + ,29 + ,107 + ,89 + ,34178 + ,9492 + ,47879 + ,43 + ,41 + ,1939 + ,129700 + ,60 + ,446 + ,80 + ,8 + ,39 + ,31 + ,124 + ,115 + ,69551 + ,17456 + ,103793 + ,53 + ,53 + ,2293 + ,110708 + ,57 + ,690 + ,137 + ,2 + ,68 + ,21 + ,78 + ,59 + ,29653 + ,9422 + ,52235 + ,82 + ,82 + ,820 + ,81518 + ,32 + ,208 + ,45 + ,0 + ,32 + ,21 + ,83 + ,66 + ,38071 + ,10913 + ,49825 + ,41 + ,41 + ,340 + ,31970 + ,15 + ,101 + ,40 + ,0 + ,5 + ,21 + ,78 + ,42 + ,4157 + ,1283 + ,4105 + ,6 + ,6 + ,2443 + ,192268 + ,102 + ,858 + ,126 + ,3 + ,53 + ,15 + ,59 + ,35 + ,28321 + ,6198 + ,58687 + ,82 + ,81 + ,993 + ,87611 + ,52 + ,293 + ,74 + ,1 + ,33 + ,9 + ,33 + ,3 + ,40195 + ,4501 + ,40745 + ,47 + ,47 + ,1038 + ,77890 + ,53 + ,349 + ,48 + ,2 + ,48 + ,21 + ,84 + ,68 + ,48158 + ,9560 + ,33187 + ,108 + ,100 + ,1380 + ,83261 + ,58 + ,411 + ,82 + ,1 + ,36 + ,18 + ,52 + ,38 + ,13310 + ,3394 + ,14063 + ,46 + ,46 + ,2186 + ,116290 + ,51 + ,561 + ,86 + ,0 + ,52 + ,31 + ,121 + ,107 + ,78474 + ,9871 + ,37407 + ,38 + ,38 + ,1069 + ,55254 + ,31 + ,289 + ,60 + ,2 + ,0 + ,24 + ,88 + ,69 + ,6386 + ,2419 + ,7190 + ,0 + ,0 + ,1763 + ,116173 + ,50 + ,492 + ,99 + ,1 + ,52 + ,24 + ,99 + ,80 + ,31588 + ,10630 + ,49562 + ,45 + ,45 + ,1995 + ,111488 + ,78 + ,669 + ,63 + ,0 + ,45 + ,22 + ,86 + ,69 + ,61254 + ,8536 + ,76324 + ,57 + ,56 + ,816 + ,60138 + ,23 + ,253 + ,76 + ,0 + ,16 + ,21 + ,75 + ,46 + ,21152 + ,4911 + ,21928 + ,20 + ,18 + ,1121 + ,73422 + ,66 + ,366 + ,92 + ,0 + ,33 + ,26 + ,96 + ,52 + ,41272 + ,9775 + ,27860 + ,56 + ,54 + ,808 + ,67751 + ,56 + ,192 + ,45 + ,0 + ,48 + ,22 + ,81 + ,58 + ,34165 + ,11227 + ,28078 + ,38 + ,37 + ,1690 + ,213351 + ,51 + ,616 + ,57 + ,0 + ,33 + ,26 + ,104 + ,85 + ,37054 + ,6916 + ,49577 + ,42 + ,40 + ,751 + ,51185 + ,24 + ,221 + ,44 + ,0 + ,24 + ,20 + ,76 + ,13 + ,12368 + ,3424 + ,28145 + ,37 + ,37 + ,1309 + ,97181 + ,32 + ,438 + ,132 + ,0 + ,37 + ,25 + ,90 + ,61 + ,23168 + ,8637 + ,36241 + ,36 + ,36 + ,685 + ,42311 + ,37 + ,229 + ,43 + ,0 + ,16 + ,19 + ,75 + ,49 + ,16380 + ,3189 + ,10824 + ,34 + ,34 + ,1326 + ,115801 + ,42 + ,388 + ,67 + ,0 + ,32 + ,22 + ,86 + ,47 + ,41242 + ,8178 + ,46892 + ,53 + ,49 + ,2224 + ,183637 + ,182 + ,536 + ,82 + ,0 + ,55 + ,25 + ,100 + ,93 + ,48450 + ,16739 + ,61264 + ,85 + ,82 + ,923 + ,68161 + ,84 + ,220 + ,71 + ,0 + ,36 + ,22 + ,88 + ,65 + ,20790 + ,6094 + ,22933 + ,36 + ,36 + ,967 + ,76441 + ,46 + ,313 + ,44 + ,4 + ,29 + ,21 + ,80 + ,64 + ,34585 + ,7237 + ,20787 + ,33 + ,33 + ,1099 + ,103613 + ,40 + ,422 + ,68 + ,0 + ,26 + ,20 + ,73 + ,64 + ,35672 + ,7355 + ,43978 + ,57 + ,55 + ,1300 + ,98707 + ,33 + ,452 + ,54 + ,3 + ,37 + ,23 + ,88 + ,57 + ,52168 + ,9734 + ,51305 + ,50 + ,50 + ,1872 + ,126527 + ,66 + ,556 + ,86 + ,1 + ,58 + ,22 + ,79 + ,61 + ,53933 + ,11225 + ,55593 + ,71 + ,71 + ,1091 + ,136781 + ,52 + ,366 + ,59 + ,0 + ,35 + ,21 + ,81 + ,71 + ,34474 + ,6213 + ,51648 + ,32 + ,31 + ,1106 + ,105863 + ,51 + ,406 + ,74 + ,0 + ,24 + ,12 + ,48 + ,43 + ,43753 + ,4875 + ,30552 + ,45 + ,42 + ,632 + ,38775 + ,30 + ,254 + ,18 + ,0 + ,18 + ,9 + ,33 + ,18 + ,36456 + ,8159 + ,23470 + ,33 + ,31 + ,1901 + ,179984 + ,89 + ,606 + ,156 + ,0 + ,37 + ,32 + ,120 + ,103 + ,51183 + ,11893 + ,77530 + ,53 + ,51 + ,1580 + ,164808 + ,49 + ,479 + ,87 + ,0 + ,86 + ,24 + ,90 + ,76 + ,52742 + ,10754 + ,57299 + ,64 + ,64 + ,223 + ,19349 + ,12 + ,67 + ,15 + ,0 + ,13 + ,1 + ,2 + ,0 + ,3895 + ,786 + ,9604 + ,14 + ,14 + ,1698 + ,143902 + ,83 + ,578 + ,104 + ,1 + ,20 + ,24 + ,96 + ,83 + ,37076 + ,9706 + ,34684 + ,38 + ,37 + ,1420 + ,108660 + ,52 + ,581 + ,49 + ,0 + ,32 + ,22 + ,79 + ,70 + ,24079 + ,7796 + ,41094 + ,39 + ,37 + ,552 + ,43803 + ,24 + ,240 + ,11 + ,0 + ,8 + ,4 + ,15 + ,4 + ,2325 + ,593 + ,3439 + ,8 + ,8 + ,708 + ,47062 + ,19 + ,219 + ,37 + ,0 + ,38 + ,15 + ,48 + ,41 + ,29354 + ,5600 + ,25171 + ,38 + ,38 + ,1079 + ,110845 + ,44 + ,349 + ,80 + ,0 + ,45 + ,21 + ,81 + ,57 + ,30341 + ,7245 + ,23437 + ,24 + ,23 + ,957 + ,92517 + ,52 + ,241 + ,66 + ,1 + ,24 + ,23 + ,84 + ,52 + ,18992 + ,7360 + ,34086 + ,22 + ,22 + ,584 + ,58660 + ,35 + ,136 + ,27 + ,0 + ,23 + ,12 + ,46 + ,24 + ,15292 + ,4574 + ,24649 + ,18 + ,18 + ,596 + ,27676 + ,22 + ,194 + ,59 + ,0 + ,2 + ,16 + ,59 + ,17 + ,5842 + ,522 + ,2342 + ,3 + ,1 + ,980 + ,98550 + ,32 + ,222 + ,113 + ,0 + ,52 + ,24 + ,96 + ,89 + ,28918 + ,10905 + ,45571 + ,49 + ,48 + ,576 + ,43284 + ,22 + ,151 + ,24 + ,0 + ,5 + ,9 + ,29 + ,20 + ,3738 + ,999 + ,3255 + ,5 + ,5 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,880 + ,66016 + ,26 + ,239 + ,54 + ,0 + ,43 + ,22 + ,79 + ,45 + ,95352 + ,9016 + ,30002 + ,47 + ,46 + ,750 + ,57359 + ,48 + ,240 + ,43 + ,0 + ,18 + ,17 + ,63 + ,63 + ,37478 + ,5134 + ,19360 + ,33 + ,33 + ,999 + ,96933 + ,35 + ,323 + ,45 + ,0 + ,41 + ,18 + ,68 + ,48 + ,26839 + ,6608 + ,43320 + ,44 + ,41 + ,931 + ,70369 + ,47 + ,302 + ,55 + ,0 + ,45 + ,21 + ,84 + ,70 + ,26783 + ,8577 + ,35513 + ,56 + ,57 + ,782 + ,65494 + ,55 + ,267 + ,66 + ,0 + ,29 + ,17 + ,54 + ,32 + ,33392 + ,1543 + ,23536 + ,49 + ,49 + ,78 + ,3616 + ,5 + ,14 + ,5 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,874 + ,143931 + ,37 + ,287 + ,67 + ,0 + ,32 + ,20 + ,75 + ,72 + ,25446 + ,9803 + ,54438 + ,45 + ,45 + ,1262 + ,109894 + ,65 + ,442 + ,67 + ,0 + ,58 + ,26 + ,87 + ,56 + ,59847 + ,12140 + ,56812 + ,78 + ,78 + ,1711 + ,122973 + ,81 + ,490 + ,115 + ,1 + ,17 + ,26 + ,104 + ,64 + ,28162 + ,6678 + ,33838 + ,51 + ,46 + ,749 + ,84336 + ,32 + ,243 + ,51 + ,0 + ,24 + ,20 + ,80 + ,77 + ,33298 + ,6420 + ,32366 + ,25 + ,25 + ,778 + ,43410 + ,19 + ,292 + ,63 + ,0 + ,7 + ,1 + ,3 + ,3 + ,2781 + ,4 + ,13 + ,1 + ,1 + ,1373 + ,136250 + ,58 + ,410 + ,84 + ,1 + ,62 + ,24 + ,93 + ,73 + ,37121 + ,7979 + ,55082 + ,62 + ,59 + ,806 + ,79015 + ,33 + ,217 + ,35 + ,0 + ,30 + ,14 + ,55 + ,37 + ,22698 + ,5141 + ,31334 + ,29 + ,29 + ,1448 + ,92937 + ,42 + ,422 + ,57 + ,8 + ,49 + ,26 + ,96 + ,54 + ,27615 + ,1311 + ,16612 + ,26 + ,26 + ,684 + ,57586 + ,37 + ,160 + ,29 + ,3 + ,3 + ,12 + ,48 + ,32 + ,32689 + ,443 + ,5084 + ,4 + ,4 + ,285 + ,19764 + ,12 + ,75 + ,19 + ,1 + ,10 + ,2 + ,8 + ,4 + ,5752 + ,2416 + ,9927 + ,10 + ,10 + ,1336 + ,105757 + ,42 + ,412 + ,51 + ,2 + ,42 + ,16 + ,60 + ,55 + ,23164 + ,8396 + ,47413 + ,43 + ,43 + ,841 + ,96410 + ,23 + ,293 + ,51 + ,0 + ,18 + ,22 + ,84 + ,81 + ,20304 + ,5462 + ,27389 + ,36 + ,36 + ,1283 + ,113402 + ,35 + ,417 + ,96 + ,0 + ,40 + ,28 + ,112 + ,90 + ,34409 + ,7271 + ,30425 + ,43 + ,41 + ,256 + ,11796 + ,9 + ,79 + ,22 + ,0 + ,1 + ,2 + ,8 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,81 + ,7627 + ,9 + ,25 + ,7 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1214 + ,121085 + ,49 + ,431 + ,34 + ,0 + ,29 + ,17 + ,52 + ,38 + ,92538 + ,4423 + ,33510 + ,33 + ,32 + ,41 + ,6836 + ,3 + ,11 + ,5 + ,0 + ,0 + ,1 + ,4 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1629 + ,139563 + ,41 + ,564 + ,43 + ,4 + ,46 + ,17 + ,57 + ,36 + ,46037 + ,5331 + ,40389 + ,53 + ,53 + ,42 + ,5118 + ,3 + ,6 + ,1 + ,0 + ,5 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,528 + ,40248 + ,16 + ,183 + ,34 + ,1 + ,8 + ,4 + ,14 + ,7 + ,5444 + ,775 + ,6012 + ,6 + ,6 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,889 + ,95079 + ,41 + ,295 + ,49 + ,0 + ,21 + ,25 + ,91 + ,75 + ,23924 + ,6676 + ,22205 + ,19 + ,18 + ,1197 + ,80750 + ,31 + ,228 + ,44 + ,0 + ,21 + ,26 + ,89 + ,52 + ,52230 + ,1489 + ,17231 + ,26 + ,26 + ,81 + ,7131 + ,4 + ,27 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,61 + ,4194 + ,11 + ,14 + ,4 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,849 + ,60378 + ,20 + ,240 + ,40 + ,1 + ,15 + ,15 + ,54 + ,45 + ,8019 + ,3080 + ,11017 + ,16 + ,16 + ,970 + ,96971 + ,40 + ,233 + ,49 + ,0 + ,40 + ,18 + ,69 + ,60 + ,34542 + ,11409 + ,46741 + ,84 + ,84 + ,964 + ,83484 + ,16 + ,347 + ,47 + ,0 + ,17 + ,19 + ,76 + ,48 + ,21157 + ,6769 + ,39869 + ,28 + ,22) + ,dim=c(15 + ,144) + ,dimnames=list(c('Pageviews' + ,'Time' + ,'Logins' + ,'CompendiumViews' + ,'CompendiumViews(PRonly)' + ,'Shared' + ,'Blogs' + ,'Reviews' + ,'Submits' + ,'Submits(+120)' + ,'Characters' + ,'CW:Revisions' + ,'CW:seconds' + ,'CW:Hyperlinks' + ,'CW:blogs') + ,1:144)) > y <- array(NA,dim=c(15,144),dimnames=list(c('Pageviews','Time','Logins','CompendiumViews','CompendiumViews(PRonly)','Shared','Blogs','Reviews','Submits','Submits(+120)','Characters','CW:Revisions','CW:seconds','CW:Hyperlinks','CW:blogs'),1:144)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par1 = 'pearson' > main = 'Correlation Matrix' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Dr. Ian E. Holliday > #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: > #Technical description: > panel.tau <- function(x, y, digits=2, prefix='', cex.cor) + { + usr <- par('usr'); on.exit(par(usr)) + par(usr = c(0, 1, 0, 1)) + rr <- cor.test(x, y, method=par1) + r <- round(rr$p.value,2) + txt <- format(c(r, 0.123456789), digits=digits)[1] + txt <- paste(prefix, txt, sep='') + if(missing(cex.cor)) cex <- 0.5/strwidth(txt) + text(0.5, 0.5, txt, cex = cex) + } > panel.hist <- function(x, ...) + { + usr <- par('usr'); on.exit(par(usr)) + par(usr = c(usr[1:2], 0, 1.5) ) + h <- hist(x, plot = FALSE) + breaks <- h$breaks; nB <- length(breaks) + y <- h$counts; y <- y/max(y) + rect(breaks[-nB], 0, breaks[-1], y, col='grey', ...) + } > postscript(file="/var/www/rcomp/tmp/1pbxw1323871646.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pairs(t(y),diag.panel=panel.hist, upper.panel=panel.smooth, lower.panel=panel.tau, main=main) > dev.off() null device 1 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/rcomp/createtable") > > n <- length(y[,1]) > n [1] 15 > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,paste('Correlations for all pairs of data series (method=',par1,')',sep=''),n+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,' ',header=TRUE) > for (i in 1:n) { + a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE) + for (j in 1:n) { + r <- cor.test(y[i,],y[j,],method=par1) + a<-table.element(a,round(r$estimate,3)) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/2t8c31323871646.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Correlations for all pairs of data series with p-values',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'pair',1,TRUE) > a<-table.element(a,'Pearson r',1,TRUE) > a<-table.element(a,'Spearman rho',1,TRUE) > a<-table.element(a,'Kendall tau',1,TRUE) > a<-table.row.end(a) > cor.test(y[1,],y[2,],method=par1) Pearson's product-moment correlation data: y[1, ] and y[2, ] t = 22.3588, df = 142, p-value < 2.2e-16 alternative hypothesis: true correlation is not equal to 0 95 percent confidence interval: 0.8401993 0.9141097 sample estimates: cor 0.8824891 > for (i in 1:(n-1)) + { + for (j in (i+1):n) + { + a<-table.row.start(a) + dum <- paste(dimnames(t(x))[[2]][i],';',dimnames(t(x))[[2]][j],sep='') + a<-table.element(a,dum,header=TRUE) + rp <- cor.test(y[i,],y[j,],method='pearson') + a<-table.element(a,round(rp$estimate,4)) + rs <- cor.test(y[i,],y[j,],method='spearman') + a<-table.element(a,round(rs$estimate,4)) + rk <- cor.test(y[i,],y[j,],method='kendall') + a<-table.element(a,round(rk$estimate,4)) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-value',header=T) + a<-table.element(a,paste('(',round(rp$p.value,4),')',sep='')) + a<-table.element(a,paste('(',round(rs$p.value,4),')',sep='')) + a<-table.element(a,paste('(',round(rk$p.value,4),')',sep='')) + a<-table.row.end(a) + } + } There were 50 or more warnings (use warnings() to see the first 50) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/3tpjb1323871646.tab") > > try(system("convert tmp/1pbxw1323871646.ps tmp/1pbxw1323871646.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.110 0.060 4.137