R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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(8 + ,350 + ,165 + ,3693 + ,11.5 + ,8 + ,318 + ,150 + ,3436 + ,11 + ,8 + ,302 + ,140 + ,3449 + ,10.5 + ,8 + ,429 + ,198 + ,4341 + ,10 + ,8 + ,440 + ,215 + ,4312 + ,8.5 + ,8 + ,455 + ,225 + ,4425 + ,10 + ,8 + ,383 + ,170 + ,3563 + ,10 + ,8 + ,340 + ,160 + ,3609 + ,8 + ,8 + ,455 + ,225 + ,3086 + ,10 + ,4 + ,113 + ,95 + ,2372 + ,15 + ,6 + ,199 + ,97 + ,2774 + ,15.5 + ,4 + ,97 + ,46 + ,1835 + ,20.5 + ,4 + ,110 + ,87 + ,2672 + ,17.5 + ,4 + ,104 + ,95 + ,2375 + ,17.5 + ,4 + ,121 + ,113 + ,2234 + ,12.5 + ,8 + ,360 + ,215 + ,4615 + ,14 + ,8 + ,307 + ,200 + ,4376 + ,15 + ,8 + ,304 + ,193 + ,4732 + ,18.5 + ,4 + ,97 + ,88 + ,2130 + ,14.5 + ,4 + ,113 + ,95 + ,2228 + ,14 + ,6 + ,250 + ,100 + ,3329 + ,15.5 + ,6 + ,232 + ,100 + ,3288 + ,15.5 + ,8 + ,350 + ,165 + ,4209 + ,12 + ,8 + ,318 + ,150 + ,4096 + ,13 + ,8 + ,400 + ,170 + ,4746 + ,12 + ,8 + ,400 + ,175 + ,5140 + ,12 + ,4 + ,140 + ,72 + ,2408 + ,19 + ,6 + ,250 + ,100 + ,3282 + ,15 + ,4 + ,122 + ,86 + ,2220 + ,14 + ,4 + ,116 + ,90 + ,2123 + ,14 + ,4 + ,88 + ,76 + ,2065 + ,14.5 + ,4 + ,71 + ,65 + ,1773 + ,19 + ,4 + ,97 + ,60 + ,1834 + ,19 + ,4 + ,91 + ,70 + ,1955 + ,20.5 + ,4 + ,97.5 + ,80 + ,2126 + ,17 + ,4 + ,122 + ,86 + ,2226 + ,16.5 + ,8 + ,350 + ,165 + ,4274 + ,12 + ,8 + ,318 + ,150 + ,4135 + ,13.5 + ,8 + ,351 + ,153 + ,4129 + ,13 + ,8 + ,429 + ,208 + ,4633 + ,11 + ,8 + ,350 + ,155 + ,4502 + ,13.5 + ,8 + ,400 + ,190 + ,4422 + ,12.5 + ,3 + ,70 + ,97 + ,2330 + ,13.5 + ,8 + ,307 + ,130 + ,4098 + ,14 + ,8 + ,302 + ,140 + ,4294 + ,16 + ,4 + ,121 + ,112 + ,2933 + ,14.5 + ,4 + ,121 + ,76 + ,2511 + ,18 + ,4 + ,122 + ,86 + ,2395 + ,16 + ,4 + ,120 + ,97 + ,2506 + ,14.5 + ,4 + ,98 + ,80 + ,2164 + ,15 + ,8 + ,350 + ,175 + ,4100 + ,13 + ,8 + ,304 + ,150 + ,3672 + ,11.5 + ,8 + ,302 + ,137 + ,4042 + ,14.5 + ,8 + ,318 + ,150 + ,3777 + ,12.5 + ,8 + ,400 + ,150 + ,4464 + ,12 + ,8 + ,351 + ,158 + ,4363 + ,13 + ,8 + ,440 + ,215 + ,4735 + ,11 + ,8 + ,455 + ,225 + ,4951 + ,11 + ,6 + ,225 + ,105 + ,3121 + ,16.5 + ,6 + ,250 + ,100 + ,3278 + ,18 + ,6 + ,250 + ,88 + ,3021 + ,16.5 + ,6 + ,198 + ,95 + ,2904 + ,16 + ,8 + ,400 + ,150 + ,4997 + ,14 + ,8 + ,350 + ,180 + ,4499 + ,12.5 + ,6 + ,232 + ,100 + ,2789 + ,15 + ,4 + ,140 + ,72 + ,2401 + ,19.5 + ,4 + ,108 + ,94 + ,2379 + ,16.5 + ,4 + ,122 + ,85 + ,2310 + ,18.5 + ,6 + ,155 + ,107 + ,2472 + ,14 + ,8 + ,350 + ,145 + ,4082 + ,13 + ,8 + ,400 + ,230 + ,4278 + ,9.5 + ,4 + ,116 + ,75 + ,2158 + ,15.5 + ,4 + ,114 + ,91 + ,2582 + ,14 + ,8 + ,318 + ,150 + ,3399 + ,11 + ,4 + ,121 + ,110 + ,2660 + ,14 + ,8 + ,350 + ,180 + ,3664 + ,11 + ,6 + ,198 + ,95 + ,3102 + ,16.5 + ,6 + ,232 + ,100 + ,2901 + ,16 + ,4 + ,122 + ,80 + ,2451 + ,16.5 + ,4 + ,71 + ,65 + ,1836 + ,21 + ,6 + ,250 + ,100 + ,3781 + ,17 + ,6 + ,258 + ,110 + ,3632 + ,18 + ,8 + ,302 + ,140 + ,4141 + ,14 + ,8 + ,350 + ,150 + ,4699 + ,14.5 + ,8 + ,302 + ,140 + ,4638 + ,16 + ,8 + ,304 + ,150 + ,4257 + ,15.5 + ,4 + ,79 + ,67 + ,1963 + ,15.5 + ,4 + ,97 + ,78 + ,2300 + ,14.5 + ,4 + ,83 + ,61 + ,2003 + ,19 + ,4 + ,90 + ,75 + ,2125 + ,14.5 + ,4 + ,116 + ,75 + ,2246 + ,14 + ,4 + ,120 + ,97 + ,2489 + ,15 + ,4 + ,79 + ,67 + ,2000 + ,16 + ,6 + ,225 + ,95 + ,3264 + ,16 + ,6 + ,250 + ,72 + ,3158 + ,19.5 + ,8 + ,400 + ,170 + ,4668 + ,11.5 + ,8 + ,350 + ,145 + ,4440 + ,14 + ,8 + ,351 + ,148 + ,4657 + ,13.5 + ,6 + ,231 + ,110 + ,3907 + ,21 + ,6 + ,258 + ,110 + ,3730 + ,19 + ,6 + ,225 + ,95 + ,3785 + ,19 + ,8 + ,262 + ,110 + ,3221 + ,13.5 + ,8 + ,302 + ,129 + ,3169 + ,12 + ,4 + ,140 + ,83 + ,2639 + ,17 + ,6 + ,232 + ,100 + ,2914 + ,16 + ,4 + ,134 + ,96 + ,2702 + ,13.5 + ,4 + ,90 + ,71 + ,2223 + ,16.5 + ,6 + ,171 + ,97 + ,2984 + ,14.5 + ,4 + ,115 + ,95 + ,2694 + ,15 + ,4 + ,120 + ,88 + ,2957 + ,17 + ,4 + ,121 + ,115 + ,2671 + ,13.5 + ,4 + ,91 + ,53 + ,1795 + ,17.5 + ,4 + ,116 + ,81 + ,2220 + ,16.9 + ,4 + ,140 + ,92 + ,2572 + ,14.9 + ,4 + ,101 + ,83 + ,2202 + ,15.3 + ,8 + ,305 + ,140 + ,4215 + ,13 + ,8 + ,304 + ,120 + ,3962 + ,13.9 + ,8 + ,351 + ,152 + ,4215 + ,12.8 + ,6 + ,250 + ,105 + ,3353 + ,14.5 + ,6 + ,200 + ,81 + ,3012 + ,17.6 + ,4 + ,85 + ,52 + ,2035 + ,22.2 + ,4 + ,98 + ,60 + ,2164 + ,22.1 + ,6 + ,225 + ,100 + ,3651 + ,17.7 + ,6 + ,250 + ,110 + ,3645 + ,16.2 + ,6 + ,258 + ,95 + ,3193 + ,17.8 + ,4 + ,85 + ,70 + ,1990 + ,17 + ,4 + ,97 + ,75 + ,2155 + ,16.4 + ,4 + ,130 + ,102 + ,3150 + ,15.7 + ,8 + ,318 + ,150 + ,3940 + ,13.2 + ,6 + ,168 + ,120 + ,3820 + ,16.7 + ,8 + ,350 + ,180 + ,4380 + ,12.1 + ,8 + ,302 + ,130 + ,3870 + ,15 + ,8 + ,318 + ,150 + ,3755 + ,14 + ,4 + ,111 + ,80 + ,2155 + ,14.8 + ,4 + ,79 + ,58 + ,1825 + ,18.6 + ,4 + ,85 + ,70 + ,1945 + ,16.8 + ,8 + ,305 + ,145 + ,3880 + ,12.5 + ,8 + ,318 + ,145 + ,4140 + ,13.7 + ,6 + ,231 + ,105 + ,3425 + ,16.9 + ,6 + ,225 + ,100 + ,3630 + ,17.7 + ,8 + ,400 + ,180 + ,4220 + ,11.1 + ,8 + ,350 + ,170 + ,4165 + ,11.4 + ,8 + ,351 + ,149 + ,4335 + ,14.5 + ,4 + ,97 + ,78 + ,1940 + ,14.5 + ,4 + ,97 + ,75 + ,2265 + ,18.2 + ,4 + ,140 + ,89 + ,2755 + ,15.8 + ,4 + ,98 + ,83 + ,2075 + ,15.9 + ,4 + ,97 + ,67 + ,1985 + ,16.4 + ,6 + ,146 + ,97 + ,2815 + ,14.5 + ,4 + ,121 + ,110 + ,2600 + ,12.8 + ,4 + ,90 + ,48 + ,1985 + ,21.5 + ,4 + ,98 + ,66 + ,1800 + ,14.4 + ,4 + ,85 + ,70 + ,2070 + ,18.6 + ,8 + ,318 + ,140 + ,3735 + ,13.2 + ,8 + ,302 + ,139 + ,3570 + ,12.8 + ,6 + ,200 + ,95 + ,3155 + ,18.2 + ,6 + ,200 + ,85 + ,2965 + ,15.8 + ,6 + ,225 + ,100 + ,3430 + ,17.2 + ,6 + ,232 + ,90 + ,3210 + ,17.2 + ,6 + ,200 + ,85 + ,3070 + ,16.7 + ,6 + ,225 + ,110 + ,3620 + ,18.7 + ,8 + ,305 + ,145 + ,3425 + ,13.2 + ,6 + ,231 + ,165 + ,3445 + ,13.4 + ,8 + ,318 + ,140 + ,4080 + ,13.7 + ,4 + ,98 + ,68 + ,2155 + ,16.5 + ,4 + ,119 + ,97 + ,2300 + ,14.7 + ,4 + ,105 + ,75 + ,2230 + ,14.5 + ,4 + ,151 + ,85 + ,2855 + ,17.6 + ,5 + ,131 + ,103 + ,2830 + ,15.9 + ,6 + ,163 + ,125 + ,3140 + ,13.6 + ,6 + ,163 + ,133 + ,3410 + ,15.8 + ,4 + ,89 + ,71 + ,1990 + ,14.9 + ,4 + ,98 + ,68 + ,2135 + ,16.6 + ,6 + ,200 + ,85 + ,2990 + ,18.2 + ,4 + ,140 + ,88 + ,2890 + ,17.3 + ,6 + ,225 + ,110 + ,3360 + ,16.6 + ,8 + ,305 + ,130 + ,3840 + ,15.4 + ,8 + ,351 + ,138 + ,3955 + ,13.2 + ,8 + ,318 + ,135 + ,3830 + ,15.2 + ,8 + ,351 + ,142 + ,4054 + ,14.3 + ,8 + ,267 + ,125 + ,3605 + ,15 + ,4 + ,89 + ,71 + ,1925 + ,14 + ,4 + ,86 + ,65 + ,1975 + ,15.2 + ,4 + ,121 + ,80 + ,2670 + ,15 + ,4 + ,141 + ,71 + ,3190 + ,24.8 + ,8 + ,260 + ,90 + ,3420 + ,22.2 + ,4 + ,105 + ,70 + ,2150 + ,14.9 + ,4 + ,85 + ,65 + ,2020 + ,19.2 + ,4 + ,151 + ,90 + ,2670 + ,16 + ,6 + ,173 + ,115 + ,2595 + ,11.3 + ,4 + ,151 + ,90 + ,2556 + ,13.2 + ,4 + ,98 + ,76 + ,2144 + ,14.7 + ,4 + ,98 + ,70 + ,2120 + ,15.5 + ,4 + ,86 + ,65 + ,2019 + ,16.4 + ,4 + ,140 + ,88 + ,2870 + ,18.1 + ,4 + ,151 + ,90 + ,3003 + ,20.1 + ,4 + ,97 + ,78 + ,2188 + ,15.8 + ,4 + ,134 + ,90 + ,2711 + ,15.5 + ,4 + ,119 + ,92 + ,2434 + ,15 + ,4 + ,108 + ,75 + ,2265 + ,15.2 + ,4 + ,156 + ,105 + ,2800 + ,14.4 + ,4 + ,85 + ,65 + ,2110 + ,19.2 + ,5 + ,121 + ,67 + ,2950 + ,19.9 + ,4 + ,91 + ,67 + ,1850 + ,13.8 + ,4 + ,89 + ,62 + ,1845 + ,15.3 + ,4 + ,122 + ,88 + ,2500 + ,15.1 + ,4 + ,135 + ,84 + ,2490 + ,15.7 + ,4 + ,151 + ,84 + ,2635 + ,16.4 + ,6 + ,173 + ,110 + ,2725 + ,12.6 + ,4 + ,135 + ,84 + ,2385 + ,12.9 + ,4 + ,86 + ,64 + ,1875 + ,16.4 + ,4 + ,81 + ,60 + ,1760 + ,16.1 + ,4 + ,85 + ,65 + ,1975 + ,19.4 + ,4 + ,89 + ,62 + ,2050 + ,17.3 + ,4 + ,105 + ,63 + ,2215 + ,14.9 + ,4 + ,98 + ,65 + ,2045 + ,16.2 + ,4 + ,105 + ,74 + ,2190 + ,14.2 + ,4 + ,119 + ,100 + ,2615 + ,14.8 + ,4 + ,141 + ,80 + ,3230 + ,20.4 + ,6 + ,146 + ,120 + ,2930 + ,13.8 + ,6 + ,231 + ,110 + ,3415 + ,15.8 + ,6 + ,200 + ,88 + ,3060 + ,17.1 + ,6 + ,225 + ,85 + ,3465 + ,16.6 + ,4 + ,112 + ,88 + ,2640 + ,18.6 + ,4 + ,112 + ,88 + ,2395 + ,18 + ,4 + ,135 + ,84 + ,2525 + ,16 + ,4 + ,151 + ,90 + ,2735 + ,18 + ,4 + ,105 + ,74 + ,1980 + ,15.3 + ,4 + ,91 + ,68 + ,1970 + ,17.6 + ,4 + ,105 + ,63 + ,2125 + ,14.7 + ,4 + ,120 + ,88 + ,2160 + ,14.5 + ,4 + ,107 + ,75 + ,2205 + ,14.5 + ,4 + ,91 + ,67 + ,1965 + ,15.7 + ,6 + ,181 + ,110 + ,2945 + ,16.4 + ,6 + ,262 + ,85 + ,3015 + ,17 + ,4 + ,144 + ,96 + ,2665 + ,13.9 + ,4 + ,151 + ,90 + ,2950 + ,17.3 + ,4 + ,140 + ,86 + ,2790 + ,15.6 + ,4 + ,135 + ,84 + ,2295 + ,11.6 + ,4 + ,120 + ,79 + ,2625 + ,18.6 + ,4 + ,72 + ,69 + ,1613 + ,18 + ,4 + ,76 + ,52 + ,1649 + ,17 + ,4 + ,79 + ,58 + ,1755 + ,17 + ,4 + ,81 + ,60 + ,1760 + ,16 + ,4 + ,71 + ,65 + ,1773 + ,19 + ,4 + ,91 + ,53 + ,1795 + ,18 + ,4 + ,91 + ,53 + ,1795 + ,17 + ,4 + ,98 + ,66 + ,1800 + ,14 + ,4 + ,91 + ,60 + ,1800 + ,16 + ,4 + ,97 + ,71 + ,1825 + ,12 + ,4 + ,79 + ,58 + ,1825 + ,19 + ,4 + ,97 + ,60 + ,1834 + ,19 + ,4 + ,97 + ,46 + ,1835 + ,21 + ,4 + ,71 + ,65 + ,1836 + ,21 + ,4 + ,89 + ,62 + ,1845 + ,15 + ,4 + ,91 + ,67 + ,1850 + ,14 + ,4 + ,68 + ,49 + ,1867 + ,20 + ,4 + ,86 + ,64 + ,1875 + ,16 + ,4 + ,98 + ,80 + ,1915 + ,14 + ,4 + ,89 + ,71 + ,1925 + ,14 + ,4 + ,90 + ,70 + ,1937 + ,14 + ,4 + ,90 + ,70 + ,1937 + ,14 + ,4 + ,97 + ,78 + ,1940 + ,15 + ,4 + ,85 + ,70 + ,1945 + ,17 + ,4 + ,97 + ,46 + ,1950 + ,21 + ,4 + ,79 + ,67 + ,1950 + ,19 + ,4 + ,91 + ,70 + ,1955 + ,21 + ,4 + ,79 + ,67 + ,1963 + ,16 + ,4 + ,91 + ,67 + ,1965 + ,15 + ,4 + ,91 + ,67 + ,1965 + ,16 + ,4 + ,89 + ,60 + ,1968 + ,19 + ,4 + ,91 + ,68 + ,1970 + ,18 + ,4 + ,86 + ,65 + ,1975 + ,15 + ,4 + ,85 + ,65 + ,1975 + ,19 + ,4 + ,105 + ,74 + ,1980 + ,15 + ,4 + ,90 + ,48 + ,1985 + ,22 + ,4 + ,97 + ,67 + ,1985 + ,16 + ,4 + ,78 + ,52 + ,1985 + ,19 + ,4 + ,91 + ,68 + ,1985 + ,16 + ,4 + ,89 + ,71 + ,1990 + ,15 + ,4 + ,85 + ,70 + ,1990 + ,17 + ,4 + ,91 + ,67 + ,1995 + ,16 + ,4 + ,79 + ,67 + ,2000 + ,16 + ,4 + ,83 + ,61 + ,2003 + ,19 + ,4 + ,86 + ,65 + ,2019 + ,16 + ,4 + ,85 + ,65 + ,2020 + ,19 + ,4 + ,91 + ,68 + ,2025 + ,18 + ,4 + ,85 + ,52 + ,2035 + ,22 + ,4 + ,98 + ,65 + ,2045 + ,16 + ,4 + ,98 + ,68 + ,2045 + ,19 + ,4 + ,89 + ,62 + ,2050 + ,17 + ,4 + ,98 + ,63 + ,2051 + ,17 + ,4 + ,88 + ,76 + ,2065 + ,15 + ,4 + ,97 + ,67 + ,2065 + ,18 + ,4 + ,85 + ,70 + ,2070 + ,19 + ,4 + ,79 + ,70 + ,2074 + ,20 + ,4 + ,98 + ,83 + ,2075 + ,16 + ,4 + ,90 + ,48 + ,2085 + ,22 + ,4 + ,97 + ,88 + ,2100 + ,17 + ,4 + ,90 + ,75 + ,2108 + ,16 + ,4 + ,86 + ,65 + ,2110 + ,18 + ,4 + ,85 + ,65 + ,2110 + ,19 + ,4 + ,98 + ,70 + ,2120 + ,16 + ,4 + ,116 + ,90 + ,2123 + ,14 + ,3 + ,70 + ,90 + ,2124 + ,14 + ,4 + ,90 + ,75 + ,2125 + ,15 + ,4 + ,105 + ,63 + ,2125 + ,15 + ,4 + ,98 + ,70 + ,2125 + ,17 + ,4 + ,97.5 + ,80 + ,2126 + ,17 + ,4 + ,91 + ,69 + ,2130 + ,15 + ,4 + ,97 + ,52 + ,2130 + ,25 + ,4 + ,97 + ,88 + ,2130 + ,15 + ,4 + ,97 + ,88 + ,2130 + ,15 + ,4 + ,98 + ,68 + ,2135 + ,17 + ,4 + ,98 + ,76 + ,2144 + ,15 + ,4 + ,97 + ,67 + ,2145 + ,18 + ,4 + ,105 + ,70 + ,2150 + ,15 + ,4 + ,111 + ,80 + ,2155 + ,15 + ,4 + ,98 + ,68 + ,2155 + ,17 + ,4 + ,97 + ,75 + ,2155 + ,16 + ,4 + ,116 + ,75 + ,2158 + ,16 + ,4 + ,120 + ,88 + ,2160 + ,15 + ,4 + ,98 + ,80 + ,2164 + ,15 + ,4 + ,98 + ,60 + ,2164 + ,22 + ,4 + ,97 + ,75 + ,2171 + ,16 + ,4 + ,97 + ,78 + ,2188 + ,16 + ,4 + ,96 + ,69 + ,2189 + ,18 + ,4 + ,97 + ,78 + ,2190 + ,14 + ,4 + ,105 + ,74 + ,2190 + ,14 + ,4 + ,105 + ,70 + ,2200 + ,13 + ,4 + ,101 + ,83 + ,2202 + ,15 + ,4 + ,107 + ,75 + ,2205 + ,15 + ,4 + ,107 + ,75 + ,2210 + ,14 + ,4 + ,105 + ,63 + ,2215 + ,15 + ,4 + ,98 + ,83 + ,2219 + ,17 + ,4 + ,122 + ,86 + ,2220 + ,14 + ,4 + ,116 + ,81 + ,2220 + ,17 + ,4 + ,90 + ,71 + ,2223 + ,17 + ,4 + ,122 + ,86 + ,2226 + ,17 + ,4 + ,113 + ,95 + ,2228 + ,14 + ,4 + ,105 + ,75 + ,2230 + ,15 + ,4 + ,121 + ,113 + ,2234 + ,13 + ,4 + ,108 + ,70 + ,2245 + ,17 + ,4 + ,116 + ,75 + ,2246 + ,14 + ,4 + ,97 + ,54 + ,2254 + ,24 + ,4 + ,98 + ,79 + ,2255 + ,18 + ,4 + ,140 + ,90 + ,2264 + ,16 + ,4 + ,98 + ,90 + ,2265 + ,16 + ,4 + ,97 + ,75 + ,2265 + ,18 + ,4 + ,108 + ,75 + ,2265 + ,15 + ,4 + ,113 + ,95 + ,2278 + ,16 + ,4 + ,97 + ,88 + ,2279 + ,19 + ,4 + ,97 + ,92 + ,2288 + ,17 + ,4 + ,107 + ,72 + ,2290 + ,17 + ,4 + ,135 + ,84 + ,2295 + ,12 + ,4 + ,122 + ,96 + ,2300 + ,16 + ,4 + ,97 + ,78 + ,2300 + ,15 + ,4 + ,119 + ,97 + ,2300 + ,15 + ,4 + ,122 + ,85 + ,2310 + ,19 + ,3 + ,70 + ,97 + ,2330 + ,14 + ,4 + ,90 + ,48 + ,2335 + ,24 + ,4 + ,108 + ,75 + ,2350 + ,17 + ,4 + ,135 + ,84 + ,2370 + ,13 + ,4 + ,113 + ,95 + ,2372 + ,15 + ,4 + ,104 + ,95 + ,2375 + ,18 + ,4 + ,108 + ,94 + ,2379 + ,17 + ,4 + ,98 + ,65 + ,2380 + ,21 + ,4 + ,135 + ,84 + ,2385 + ,13 + ,4 + ,108 + ,93 + ,2391 + ,16 + ,4 + ,122 + ,86 + ,2395 + ,16 + ,4 + ,112 + ,88 + ,2395 + ,18 + ,4 + ,140 + ,72 + ,2401 + ,20 + ,4 + ,119 + ,97 + ,2405 + ,15 + ,4 + ,140 + ,72 + ,2408 + ,19 + ,4 + ,140 + ,90 + ,2408 + ,20 + ,3 + ,70 + ,100 + ,2420 + ,13 + ,4 + ,107 + ,90 + ,2430 + ,15 + ,4 + ,119 + ,92 + ,2434 + ,15 + ,4 + ,122 + ,80 + ,2451 + ,17 + ,4 + ,107 + ,86 + ,2464 + ,16 + ,6 + ,155 + ,107 + ,2472 + ,14 + ,4 + ,120 + ,97 + ,2489 + ,15 + ,4 + ,135 + ,84 + ,2490 + ,16 + ,4 + ,122 + ,88 + ,2500 + ,15 + ,4 + ,120 + ,97 + ,2506 + ,15 + ,4 + ,121 + ,76 + ,2511 + ,18 + ,4 + ,134 + ,95 + ,2515 + ,15 + ,4 + ,135 + ,84 + ,2525 + ,16 + ,4 + ,140 + ,75 + ,2542 + ,17 + ,4 + ,120 + ,75 + ,2542 + ,18 + ,4 + ,119 + ,97 + ,2545 + ,17 + ,4 + ,151 + ,90 + ,2556 + ,13 + ,4 + ,134 + ,95 + ,2560 + ,14 + ,4 + ,140 + ,72 + ,2565 + ,14 + ,4 + ,140 + ,92 + ,2572 + ,15 + ,4 + ,112 + ,85 + ,2575 + ,16 + ,4 + ,114 + ,91 + ,2582 + ,14 + ,4 + ,156 + ,92 + ,2585 + ,15 + ,6 + ,200 + ,85 + ,2587 + ,16 + ,4 + ,140 + ,78 + ,2592 + ,19 + ,6 + ,173 + ,115 + ,2595 + ,11 + ,4 + ,121 + ,110 + ,2600 + ,13 + ,4 + ,112 + ,88 + ,2605 + ,20 + ,4 + ,119 + ,100 + ,2615 + ,15 + ,4 + ,156 + ,92 + ,2620 + ,14 + ,4 + ,120 + ,79 + ,2625 + ,19 + ,6 + ,232 + ,100 + ,2634 + ,13 + ,4 + ,151 + ,84 + ,2635 + ,16 + ,4 + ,120 + ,74 + ,2635 + ,18 + ,4 + ,140 + ,83 + ,2639 + ,17 + ,4 + ,112 + ,88 + ,2640 + ,19 + ,6 + ,199 + ,90 + ,2648 + ,15 + ,4 + ,121 + ,110 + ,2660 + ,14 + ,4 + ,144 + ,96 + ,2665 + ,14 + ,4 + ,121 + ,80 + ,2670 + ,15 + ,4 + ,151 + ,90 + ,2670 + ,16 + ,4 + ,121 + ,115 + ,2671 + ,14 + ,4 + ,110 + ,87 + ,2672 + ,18 + ,4 + ,151 + ,90 + ,2678 + ,17 + ,4 + ,115 + ,95 + ,2694 + ,15 + ,6 + ,173 + ,115 + ,2700 + ,13 + ,4 + ,134 + ,96 + ,2702 + ,14 + ,4 + ,134 + ,90 + ,2711 + ,16 + ,4 + ,140 + ,88 + ,2720 + ,15 + ,4 + ,119 + ,82 + ,2720 + ,19 + ,3 + ,80 + ,110 + ,2720 + ,14 + ,6 + ,173 + ,110 + ,2725 + ,13 + ,4 + ,151 + ,90 + ,2735 + ,18 + ,4 + ,151 + ,88 + ,2740 + ,16 + ,4 + ,156 + ,105 + ,2745 + ,17 + ,4 + ,140 + ,89 + ,2755 + ,16 + ,6 + ,199 + ,97 + ,2774 + ,16 + ,6 + ,232 + ,100 + ,2789 + ,15 + ,4 + ,140 + ,86 + ,2790 + ,16 + ,4 + ,121 + ,115 + ,2795 + ,16 + ,4 + ,156 + ,105 + ,2800 + ,14 + ,6 + ,156 + ,122 + ,2807 + ,14 + ,6 + ,146 + ,97 + ,2815 + ,15 + ,5 + ,131 + ,103 + ,2830 + ,16 + ,6 + ,198 + ,95 + ,2833 + ,16 + ,6 + ,232 + ,112 + ,2835 + ,15 + ,4 + ,151 + ,85 + ,2855 + ,18 + ,4 + ,140 + ,92 + ,2865 + ,16 + ,4 + ,121 + ,112 + ,2868 + ,16 + ,4 + ,140 + ,88 + ,2870 + ,18 + ,4 + ,140 + ,88 + ,2890 + ,17 + ,6 + ,168 + ,116 + ,2900 + ,13 + ,6 + ,232 + ,100 + ,2901 + ,16 + ,6 + ,198 + ,95 + ,2904 + ,16 + ,6 + ,168 + ,132 + ,2910 + ,11 + ,6 + ,232 + ,100 + ,2914 + ,16 + ,6 + ,156 + ,108 + ,2930 + ,16 + ,6 + ,146 + ,120 + ,2930 + ,14 + ,4 + ,121 + ,112 + ,2933 + ,15 + ,6 + ,232 + ,100 + ,2945 + ,16 + ,6 + ,181 + ,110 + ,2945 + ,16 + ,4 + ,121 + ,98 + ,2945 + ,15 + ,4 + ,151 + ,90 + ,2950 + ,17 + ,5 + ,121 + ,67 + ,2950 + ,20 + ,4 + ,120 + ,88 + ,2957 + ,17 + ,6 + ,258 + ,110 + ,2962 + ,14 + ,6 + ,200 + ,85 + ,2965 + ,16 + ,4 + ,120 + ,87 + ,2979 + ,20 + ,6 + ,171 + ,97 + ,2984 + ,15 + ,6 + ,200 + ,85 + ,2990 + ,18 + ,4 + ,151 + ,90 + ,3003 + ,20 + ,6 + ,200 + ,81 + ,3012 + ,18 + ,6 + ,262 + ,85 + ,3015 + ,17 + ,6 + ,250 + ,88 + ,3021 + ,17 + ,6 + ,231 + ,110 + ,3039 + ,15 + ,6 + ,200 + ,88 + ,3060 + ,17 + ,6 + ,200 + ,85 + ,3070 + ,17 + ,6 + ,232 + ,90 + ,3085 + ,18 + ,8 + ,455 + ,225 + ,3086 + ,10 + ,6 + ,198 + ,95 + ,3102 + ,17 + ,6 + ,225 + ,105 + ,3121 + ,17 + ,6 + ,250 + ,88 + ,3139 + ,15 + ,6 + ,163 + ,125 + ,3140 + ,14 + ,4 + ,130 + ,102 + ,3150 + ,16 + ,6 + ,200 + ,95 + ,3155 + ,18 + ,6 + ,250 + ,72 + ,3158 + ,20 + ,6 + ,145 + ,76 + ,3160 + ,20 + ,8 + ,302 + ,129 + ,3169 + ,12 + ,4 + ,141 + ,71 + ,3190 + ,25 + ,6 + ,258 + ,95 + ,3193 + ,18 + ,8 + ,302 + ,139 + ,3205 + ,11 + ,6 + ,232 + ,90 + ,3210 + ,17 + ,6 + ,232 + ,90 + ,3211 + ,17 + ,8 + ,262 + ,110 + ,3221 + ,14 + ,4 + ,141 + ,80 + ,3230 + ,20 + ,6 + ,225 + ,100 + ,3233 + ,15 + ,6 + ,231 + ,115 + ,3245 + ,15 + ,4 + ,146 + ,67 + ,3250 + ,22 + ,6 + ,225 + ,95 + ,3264 + ,16 + ,6 + ,232 + ,90 + ,3265 + ,18 + ,4 + ,120 + ,88 + ,3270 + ,22 + ,6 + ,250 + ,100 + ,3278 + ,18 + ,6 + ,250 + ,100 + ,3282 + ,15 + ,6 + ,232 + ,100 + ,3288 + ,16 + ,6 + ,250 + ,88 + ,3302 + ,16 + ,6 + ,250 + ,100 + ,3329 + ,16 + ,6 + ,250 + ,100 + ,3336 + ,17 + ,6 + ,250 + ,105 + ,3353 + ,15 + ,6 + ,225 + ,110 + ,3360 + ,17 + ,8 + ,260 + ,110 + ,3365 + ,16 + ,6 + ,231 + ,105 + ,3380 + ,16 + ,6 + ,225 + ,90 + ,3381 + ,19 + ,8 + ,318 + ,150 + ,3399 + ,11 + ,6 + ,258 + ,120 + ,3410 + ,15 + ,6 + ,163 + ,133 + ,3410 + ,16 + ,6 + ,231 + ,110 + ,3415 + ,16 + ,8 + ,260 + ,90 + ,3420 + ,22 + ,6 + ,231 + ,105 + ,3425 + ,17 + ,8 + ,305 + ,145 + ,3425 + ,13 + ,6 + ,225 + ,100 + ,3430 + ,17 + ,6 + ,250 + ,72 + ,3432 + ,21 + ,8 + ,304 + ,150 + ,3433 + ,12 + ,8 + ,318 + ,150 + ,3436 + ,11 + ,6 + ,225 + ,105 + ,3439 + ,16 + ,6 + ,231 + ,165 + ,3445 + ,13 + ,8 + ,302 + ,140 + ,3449 + ,11 + ,6 + ,250 + ,105 + ,3459 + ,16 + ,6 + ,225 + ,85 + ,3465 + ,17 + ,8 + ,307 + ,130 + ,3504 + ,12 + ,6 + ,250 + ,110 + ,3520 + ,16 + ,6 + ,250 + ,98 + ,3525 + ,19 + ,5 + ,183 + ,77 + ,3530 + ,20 + ,6 + ,231 + ,105 + ,3535 + ,19 + ,8 + ,383 + ,170 + ,3563 + ,10 + ,8 + ,302 + ,139 + ,3570 + ,13 + ,6 + ,250 + ,78 + ,3574 + ,21 + ,8 + ,267 + ,125 + ,3605 + ,15 + ,8 + ,340 + ,160 + ,3609 + ,8 + ,6 + ,225 + ,105 + ,3613 + ,17 + ,6 + ,225 + ,110 + ,3620 + ,19 + ,6 + ,225 + ,100 + ,3630 + ,18 + ,6 + ,258 + ,110 + ,3632 + ,18 + ,6 + ,250 + ,110 + ,3645 + ,16 + ,6 + ,225 + ,100 + ,3651 + ,18 + ,8 + ,350 + ,180 + ,3664 + ,11 + ,8 + ,304 + ,150 + ,3672 + ,12 + ,8 + ,304 + ,150 + ,3672 + ,12 + ,8 + ,350 + ,165 + ,3693 + ,12 + ,8 + ,302 + ,129 + ,3725 + ,13 + ,8 + ,350 + ,105 + ,3725 + ,19 + ,6 + ,258 + ,110 + ,3730 + ,19 + ,8 + ,318 + ,140 + ,3735 + ,13 + ,8 + ,318 + ,150 + ,3755 + ,14 + ,8 + ,400 + ,150 + ,3761 + ,10 + ,8 + ,318 + ,150 + ,3777 + ,13 + ,6 + ,250 + ,100 + ,3781 + ,17 + ,6 + ,225 + ,95 + ,3785 + ,19 + ,6 + ,168 + ,120 + ,3820 + ,17 + ,8 + ,360 + ,175 + ,3821 + ,11 + ,8 + ,318 + ,135 + ,3830 + ,15 + ,8 + ,305 + ,130 + ,3840 + ,15 + ,8 + ,390 + ,190 + ,3850 + ,9 + ,8 + ,302 + ,130 + ,3870 + ,15 + ,8 + ,305 + ,145 + ,3880 + ,13 + ,8 + ,304 + ,150 + ,3892 + ,13 + ,6 + ,250 + ,105 + ,3897 + ,19 + ,8 + ,350 + ,125 + ,3900 + ,17 + ,6 + ,231 + ,110 + ,3907 + ,21 + ,8 + ,318 + ,150 + ,3940 + ,13 + ,8 + ,360 + ,150 + ,3940 + ,13 + ,8 + ,351 + ,138 + ,3955 + ,13 + ,8 + ,304 + ,120 + ,3962 + ,14 + ,8 + ,350 + ,145 + ,3988 + ,13 + ,8 + ,302 + ,137 + ,4042 + ,15 + ,8 + ,351 + ,142 + ,4054 + ,14 + ,8 + ,350 + ,145 + ,4055 + ,12 + ,8 + ,260 + ,110 + ,4060 + ,19 + ,8 + ,318 + ,150 + ,4077 + ,14 + ,8 + ,318 + ,140 + ,4080 + ,14 + ,8 + ,350 + ,145 + ,4082 + ,13 + ,8 + ,318 + ,150 + ,4096 + ,13 + ,8 + ,307 + ,130 + ,4098 + ,14 + ,8 + ,350 + ,175 + ,4100 + ,13 + ,8 + ,351 + ,153 + ,4129 + ,13 + ,8 + ,318 + ,150 + ,4135 + ,14 + ,8 + ,318 + ,145 + ,4140 + ,14 + ,8 + ,302 + ,140 + ,4141 + ,14 + ,8 + ,351 + ,153 + ,4154 + ,14 + ,8 + ,350 + ,170 + ,4165 + ,11 + ,8 + ,318 + ,150 + ,4190 + ,13 + ,8 + ,350 + ,165 + ,4209 + ,12 + ,8 + ,305 + ,140 + ,4215 + ,13 + ,8 + ,351 + ,152 + ,4215 + ,13 + ,8 + ,400 + ,180 + ,4220 + ,11 + ,8 + ,318 + ,150 + ,4237 + ,15 + ,8 + ,304 + ,150 + ,4257 + ,16 + ,8 + ,350 + ,165 + ,4274 + ,12 + ,8 + ,400 + ,230 + ,4278 + ,10 + ,8 + ,302 + ,140 + ,4294 + ,16 + ,8 + ,302 + ,130 + ,4295 + ,15 + ,8 + ,440 + ,215 + ,4312 + ,9 + ,8 + ,400 + ,190 + ,4325 + ,12 + ,8 + ,351 + ,149 + ,4335 + ,15 + ,8 + ,429 + ,198 + ,4341 + ,10 + ,8 + ,454 + ,220 + ,4354 + ,9 + ,8 + ,350 + ,155 + ,4360 + ,15 + ,8 + ,351 + ,158 + ,4363 + ,13 + ,8 + ,307 + ,200 + ,4376 + ,15 + ,8 + ,350 + ,180 + ,4380 + ,12 + ,8 + ,318 + ,210 + ,4382 + ,14 + ,8 + ,400 + ,175 + ,4385 + ,12 + ,8 + ,400 + ,190 + ,4422 + ,13 + ,8 + ,455 + ,225 + ,4425 + ,10 + ,8 + ,350 + ,145 + ,4440 + ,14 + ,8 + ,350 + ,160 + ,4456 + ,14 + ,8 + ,318 + ,150 + ,4457 + ,14 + ,8 + ,400 + ,175 + ,4464 + ,12 + ,8 + ,400 + ,150 + ,4464 + ,12 + ,8 + ,318 + ,150 + ,4498 + ,15 + ,8 + ,350 + ,180 + ,4499 + ,13 + ,8 + ,350 + ,155 + ,4502 + ,14 + ,8 + ,360 + ,215 + ,4615 + ,14 + ,8 + ,429 + ,208 + ,4633 + ,11 + ,8 + ,302 + ,140 + ,4638 + ,16 + ,8 + ,360 + ,170 + ,4654 + ,13 + ,8 + ,351 + ,148 + ,4657 + ,14 + ,8 + ,400 + ,170 + ,4668 + ,12 + ,8 + ,350 + ,150 + ,4699 + ,15 + ,8 + ,304 + ,193 + ,4732 + ,19 + ,8 + ,440 + ,215 + ,4735 + ,11 + ,8 + ,400 + ,170 + ,4746 + ,12 + ,8 + ,400 + ,167 + ,4906 + ,13 + ,8 + ,455 + ,225 + ,4951 + ,11 + ,8 + ,429 + ,198 + ,4952 + ,12 + ,8 + ,383 + ,180 + ,4955 + ,12 + ,8 + ,400 + ,150 + ,4997 + ,14 + ,8 + ,400 + ,175 + ,5140 + ,12) + ,dim=c(5 + ,632) + ,dimnames=list(c('cylinders' + ,'engine.displacement' + ,'horsepower' + ,'weight' + ,'acceleration ') + ,1:632)) > y <- array(NA,dim=c(5,632),dimnames=list(c('cylinders','engine.displacement','horsepower','weight','acceleration '),1:632)) > 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: Patrick Wessa > #To cite this work: Patrick Wessa, (2010), Multivariate Correlation Matrix (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/Patrick.Wessa/rwasp_pairs.wasp#output/ > #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/html/rcomp/tmp/11ptz1292951448.ps",horizontal=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) Warning messages: 1: In strwidth(labels, "user") : font width unknown for character 0xd 2: In text.default(x, y, txt, cex = cex, font = font) : font width unknown for character 0xd 3: In text.default(x, y, txt, cex = cex, font = font) : font metrics unknown for character 0xd > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > n <- length(y[,1]) > n [1] 5 > 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/html/rcomp/tmp/2npr51292951448.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 = 76.5722, df = 630, p-value < 2.2e-16 alternative hypothesis: true correlation is not equal to 0 95 percent confidence interval: 0.9420785 0.9572958 sample estimates: cor 0.950251 > 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) + } + } Warning messages: 1: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 2: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 3: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 4: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 5: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 6: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 7: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 8: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 9: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties 10: In cor.test.default(y[i, ], y[j, ], method = "spearman") : Cannot compute exact p-values with ties > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/31h7e1292951448.tab") > > try(system("convert tmp/11ptz1292951448.ps tmp/11ptz1292951448.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.020 0.618 8.452