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(0 + ,1979900 + ,7 + ,0 + ,1737275 + ,1742117 + ,1522522 + ,1276674 + ,1086979 + ,1149822 + ,0 + ,2061036 + ,8 + ,0 + ,1979900 + ,1737275 + ,1742117 + ,1522522 + ,1276674 + ,1086979 + ,0 + ,1867943 + ,9 + ,0 + ,2061036 + ,1979900 + ,1737275 + ,1742117 + ,1522522 + ,1276674 + ,0 + ,1707752 + ,10 + ,0 + ,1867943 + ,2061036 + ,1979900 + ,1737275 + ,1742117 + ,1522522 + ,0 + ,1298756 + ,11 + ,0 + ,1707752 + ,1867943 + ,2061036 + ,1979900 + ,1737275 + ,1742117 + ,0 + ,1281814 + ,12 + ,0 + ,1298756 + ,1707752 + ,1867943 + ,2061036 + ,1979900 + ,1737275 + ,0 + ,1281151 + ,13 + ,0 + ,1281814 + ,1298756 + ,1707752 + ,1867943 + ,2061036 + ,1979900 + ,0 + ,1164976 + ,14 + ,0 + ,1281151 + ,1281814 + ,1298756 + ,1707752 + ,1867943 + ,2061036 + ,0 + ,1454329 + ,15 + ,0 + ,1164976 + ,1281151 + ,1281814 + ,1298756 + ,1707752 + ,1867943 + ,0 + ,1645288 + ,16 + ,0 + ,1454329 + ,1164976 + ,1281151 + ,1281814 + ,1298756 + ,1707752 + ,0 + ,1817743 + ,17 + ,0 + ,1645288 + ,1454329 + ,1164976 + ,1281151 + ,1281814 + ,1298756 + ,0 + ,1895785 + ,18 + ,0 + ,1817743 + ,1645288 + ,1454329 + ,1164976 + ,1281151 + ,1281814 + ,0 + ,2236311 + ,19 + ,0 + ,1895785 + ,1817743 + ,1645288 + ,1454329 + ,1164976 + ,1281151 + ,0 + ,2295951 + ,20 + ,0 + ,2236311 + ,1895785 + ,1817743 + ,1645288 + ,1454329 + ,1164976 + ,0 + ,2087315 + ,21 + ,0 + ,2295951 + ,2236311 + ,1895785 + ,1817743 + ,1645288 + ,1454329 + ,0 + ,1980891 + ,22 + ,0 + ,2087315 + ,2295951 + ,2236311 + ,1895785 + ,1817743 + ,1645288 + ,0 + ,1465446 + ,23 + ,0 + ,1980891 + ,2087315 + ,2295951 + ,2236311 + ,1895785 + ,1817743 + ,0 + ,1445026 + ,24 + ,0 + ,1465446 + ,1980891 + ,2087315 + ,2295951 + ,2236311 + ,1895785 + ,0 + ,1488120 + ,25 + ,0 + ,1445026 + ,1465446 + ,1980891 + ,2087315 + ,2295951 + ,2236311 + ,0 + ,1338333 + ,26 + ,0 + ,1488120 + ,1445026 + ,1465446 + ,1980891 + ,2087315 + ,2295951 + ,0 + ,1715789 + ,27 + ,0 + ,1338333 + ,1488120 + ,1445026 + ,1465446 + ,1980891 + ,2087315 + ,0 + ,1806090 + ,28 + ,0 + ,1715789 + ,1338333 + ,1488120 + ,1445026 + ,1465446 + ,1980891 + ,0 + ,2083316 + ,29 + ,0 + ,1806090 + ,1715789 + ,1338333 + ,1488120 + ,1445026 + ,1465446 + ,0 + ,2092278 + ,30 + ,0 + ,2083316 + ,1806090 + ,1715789 + ,1338333 + ,1488120 + ,1445026 + ,0 + ,2430800 + ,31 + ,0 + ,2092278 + ,2083316 + ,1806090 + ,1715789 + ,1338333 + ,1488120 + ,0 + ,2424894 + ,32 + ,0 + ,2430800 + ,2092278 + ,2083316 + ,1806090 + ,1715789 + ,1338333 + ,0 + ,2299016 + ,33 + ,0 + ,2424894 + ,2430800 + ,2092278 + ,2083316 + ,1806090 + ,1715789 + ,0 + ,2130688 + ,34 + ,0 + ,2299016 + ,2424894 + ,2430800 + ,2092278 + ,2083316 + ,1806090 + ,0 + ,1652221 + ,35 + ,0 + ,2130688 + ,2299016 + ,2424894 + ,2430800 + ,2092278 + ,2083316 + ,0 + ,1608162 + ,36 + ,0 + ,1652221 + ,2130688 + ,2299016 + ,2424894 + ,2430800 + ,2092278 + ,0 + ,1647074 + ,37 + ,0 + ,1608162 + ,1652221 + ,2130688 + ,2299016 + ,2424894 + ,2430800 + ,0 + ,1479691 + ,38 + ,0 + ,1647074 + ,1608162 + ,1652221 + ,2130688 + ,2299016 + ,2424894 + ,0 + ,1884978 + ,39 + ,0 + ,1479691 + ,1647074 + ,1608162 + ,1652221 + ,2130688 + ,2299016 + ,0 + ,2007898 + ,40 + ,0 + ,1884978 + ,1479691 + ,1647074 + ,1608162 + ,1652221 + ,2130688 + ,0 + ,2208954 + ,41 + ,0 + ,2007898 + ,1884978 + ,1479691 + ,1647074 + ,1608162 + ,1652221 + ,0 + ,2217164 + ,42 + ,0 + ,2208954 + ,2007898 + ,1884978 + ,1479691 + ,1647074 + ,1608162 + ,0 + ,2534291 + ,43 + ,0 + ,2217164 + ,2208954 + ,2007898 + ,1884978 + ,1479691 + ,1647074 + ,0 + ,2560312 + ,44 + ,0 + ,2534291 + ,2217164 + ,2208954 + ,2007898 + ,1884978 + ,1479691 + ,0 + ,2429069 + ,45 + ,0 + ,2560312 + ,2534291 + ,2217164 + ,2208954 + ,2007898 + ,1884978 + ,0 + ,2315077 + ,46 + ,0 + ,2429069 + ,2560312 + ,2534291 + ,2217164 + ,2208954 + ,2007898 + ,0 + ,1799608 + ,47 + ,0 + ,2315077 + ,2429069 + ,2560312 + ,2534291 + ,2217164 + ,2208954 + ,0 + ,1772590 + ,48 + ,0 + ,1799608 + ,2315077 + ,2429069 + ,2560312 + ,2534291 + ,2217164 + ,0 + ,1744799 + ,49 + ,0 + ,1772590 + ,1799608 + ,2315077 + ,2429069 + ,2560312 + ,2534291 + ,0 + ,1659093 + ,50 + ,0 + ,1744799 + ,1772590 + ,1799608 + ,2315077 + ,2429069 + ,2560312 + ,0 + ,2099821 + ,51 + ,0 + ,1659093 + ,1744799 + ,1772590 + ,1799608 + ,2315077 + ,2429069 + ,0 + ,2135736 + ,52 + ,0 + ,2099821 + ,1659093 + ,1744799 + ,1772590 + ,1799608 + ,2315077 + ,0 + ,2427894 + ,53 + ,0 + ,2135736 + ,2099821 + ,1659093 + ,1744799 + ,1772590 + ,1799608 + ,0 + ,2468882 + ,54 + ,0 + ,2427894 + ,2135736 + ,2099821 + ,1659093 + ,1744799 + ,1772590 + ,0 + ,2703217 + ,55 + ,0 + ,2468882 + ,2427894 + ,2135736 + ,2099821 + ,1659093 + ,1744799 + ,0 + ,2766841 + ,56 + ,0 + ,2703217 + ,2468882 + ,2427894 + ,2135736 + ,2099821 + ,1659093 + ,0 + ,2655236 + ,57 + ,0 + ,2766841 + ,2703217 + ,2468882 + ,2427894 + ,2135736 + ,2099821 + ,0 + ,2550373 + ,58 + ,0 + ,2655236 + ,2766841 + ,2703217 + ,2468882 + ,2427894 + ,2135736 + ,0 + ,2052097 + ,59 + ,0 + ,2550373 + ,2655236 + ,2766841 + ,2703217 + ,2468882 + ,2427894 + ,0 + ,1998055 + ,60 + ,0 + ,2052097 + ,2550373 + ,2655236 + ,2766841 + ,2703217 + ,2468882 + ,0 + ,1920748 + ,61 + ,0 + ,1998055 + ,2052097 + ,2550373 + ,2655236 + ,2766841 + ,2703217 + ,0 + ,1876694 + ,62 + ,0 + ,1920748 + ,1998055 + ,2052097 + ,2550373 + ,2655236 + ,2766841 + ,0 + ,2380930 + ,63 + ,0 + ,1876694 + ,1920748 + ,1998055 + ,2052097 + ,2550373 + ,2655236 + ,0 + ,2467402 + ,64 + ,0 + ,2380930 + ,1876694 + ,1920748 + ,1998055 + ,2052097 + ,2550373 + ,0 + ,2770771 + ,65 + ,0 + ,2467402 + ,2380930 + ,1876694 + ,1920748 + ,1998055 + ,2052097 + ,0 + ,2781340 + ,66 + ,0 + ,2770771 + ,2467402 + ,2380930 + ,1876694 + ,1920748 + ,1998055 + ,0 + ,3143926 + ,67 + ,0 + ,2781340 + ,2770771 + ,2467402 + ,2380930 + ,1876694 + ,1920748 + ,0 + ,3172235 + ,68 + ,0 + ,3143926 + ,2781340 + ,2770771 + ,2467402 + ,2380930 + ,1876694 + ,0 + ,2952540 + ,69 + ,0 + ,3172235 + ,3143926 + ,2781340 + ,2770771 + ,2467402 + ,2380930 + ,0 + ,2920877 + ,70 + ,0 + ,2952540 + ,3172235 + ,3143926 + ,2781340 + ,2770771 + ,2467402 + ,0 + ,2384552 + ,71 + ,0 + ,2920877 + ,2952540 + ,3172235 + ,3143926 + ,2781340 + ,2770771 + ,0 + ,2248987 + ,72 + ,0 + ,2384552 + ,2920877 + ,2952540 + ,3172235 + ,3143926 + ,2781340 + ,0 + ,2208616 + ,73 + ,0 + ,2248987 + ,2384552 + ,2920877 + ,2952540 + ,3172235 + ,3143926 + ,0 + ,2178756 + ,74 + ,0 + ,2208616 + ,2248987 + ,2384552 + ,2920877 + ,2952540 + ,3172235 + ,0 + ,2632870 + ,75 + ,0 + ,2178756 + ,2208616 + ,2248987 + ,2384552 + ,2920877 + ,2952540 + ,0 + ,2706905 + ,76 + ,0 + ,2632870 + ,2178756 + ,2208616 + ,2248987 + ,2384552 + ,2920877 + ,0 + ,3029745 + ,77 + ,0 + ,2706905 + ,2632870 + ,2178756 + ,2208616 + ,2248987 + ,2384552 + ,0 + ,3015402 + ,78 + ,0 + ,3029745 + ,2706905 + ,2632870 + ,2178756 + ,2208616 + ,2248987 + ,0 + ,3391414 + ,79 + ,0 + ,3015402 + ,3029745 + ,2706905 + ,2632870 + ,2178756 + ,2208616 + ,0 + ,3507805 + ,80 + ,0 + ,3391414 + ,3015402 + ,3029745 + ,2706905 + ,2632870 + ,2178756 + ,0 + ,3177852 + ,81 + ,0 + ,3507805 + ,3391414 + ,3015402 + ,3029745 + ,2706905 + ,2632870 + ,0 + ,3142961 + ,82 + ,0 + ,3177852 + ,3507805 + ,3391414 + ,3015402 + ,3029745 + ,2706905 + ,0 + ,2545815 + ,83 + ,0 + ,3142961 + ,3177852 + ,3507805 + ,3391414 + ,3015402 + ,3029745 + ,0 + ,2414007 + ,84 + ,0 + ,2545815 + ,3142961 + ,3177852 + ,3507805 + ,3391414 + ,3015402 + ,0 + ,2372578 + ,85 + ,0 + ,2414007 + ,2545815 + ,3142961 + ,3177852 + ,3507805 + ,3391414 + ,0 + ,2332664 + ,86 + ,0 + ,2372578 + ,2414007 + ,2545815 + ,3142961 + ,3177852 + ,3507805 + ,0 + ,2825328 + ,87 + ,0 + ,2332664 + ,2372578 + ,2414007 + ,2545815 + ,3142961 + ,3177852 + ,0 + ,2901478 + ,88 + ,0 + ,2825328 + ,2332664 + ,2372578 + ,2414007 + ,2545815 + ,3142961 + ,0 + ,3263955 + ,89 + ,0 + ,2901478 + ,2825328 + ,2332664 + ,2372578 + ,2414007 + ,2545815 + ,0 + ,3226738 + ,90 + ,0 + ,3263955 + ,2901478 + ,2825328 + ,2332664 + ,2372578 + ,2414007 + ,0 + ,3610786 + ,91 + ,0 + ,3226738 + ,3263955 + ,2901478 + ,2825328 + ,2332664 + ,2372578 + ,0 + ,3709274 + ,92 + ,0 + ,3610786 + ,3226738 + ,3263955 + ,2901478 + ,2825328 + ,2332664 + ,0 + ,3467185 + ,93 + ,0 + ,3709274 + ,3610786 + ,3226738 + ,3263955 + ,2901478 + ,2825328 + ,0 + ,3449646 + ,94 + ,0 + ,3467185 + ,3709274 + ,3610786 + ,3226738 + ,3263955 + ,2901478 + ,0 + ,2802951 + ,95 + ,0 + ,3449646 + ,3467185 + ,3709274 + ,3610786 + ,3226738 + ,3263955 + ,0 + ,2462530 + ,96 + ,0 + ,2802951 + ,3449646 + ,3467185 + ,3709274 + ,3610786 + ,3226738 + ,0 + ,2490645 + ,97 + ,0 + ,2462530 + ,2802951 + ,3449646 + ,3467185 + ,3709274 + ,3610786 + ,0 + ,2561520 + ,98 + ,0 + ,2490645 + ,2462530 + ,2802951 + ,3449646 + ,3467185 + ,3709274 + ,0 + ,3067554 + ,99 + ,0 + ,2561520 + ,2490645 + ,2462530 + ,2802951 + ,3449646 + ,3467185 + ,0 + ,3226951 + ,100 + ,0 + ,3067554 + ,2561520 + ,2490645 + ,2462530 + ,2802951 + ,3449646 + ,0 + ,3546493 + ,101 + ,0 + ,3226951 + ,3067554 + ,2561520 + ,2490645 + ,2462530 + ,2802951 + ,0 + ,3492787 + ,102 + ,0 + ,3546493 + ,3226951 + ,3067554 + ,2561520 + ,2490645 + ,2462530 + ,0 + ,3952263 + ,103 + ,0 + ,3492787 + ,3546493 + ,3226951 + ,3067554 + ,2561520 + ,2490645 + ,0 + ,3932072 + ,104 + ,0 + ,3952263 + ,3492787 + ,3546493 + ,3226951 + ,3067554 + ,2561520 + ,0 + ,3720284 + ,105 + ,0 + ,3932072 + ,3952263 + ,3492787 + ,3546493 + ,3226951 + ,3067554 + ,0 + ,3651555 + ,106 + ,0 + ,3720284 + ,3932072 + ,3952263 + ,3492787 + ,3546493 + ,3226951 + ,0 + ,2914972 + ,107 + ,0 + ,3651555 + ,3720284 + ,3932072 + ,3952263 + ,3492787 + ,3546493 + ,0 + ,2713514 + ,108 + ,0 + ,2914972 + ,3651555 + ,3720284 + ,3932072 + ,3952263 + ,3492787 + ,0 + ,2703997 + ,109 + ,0 + ,2713514 + ,2914972 + ,3651555 + ,3720284 + ,3932072 + ,3952263 + ,0 + ,2591373 + ,110 + ,0 + ,2703997 + ,2713514 + ,2914972 + ,3651555 + ,3720284 + ,3932072 + ,0 + ,3163748 + ,111 + ,0 + ,2591373 + ,2703997 + ,2713514 + ,2914972 + ,3651555 + ,3720284 + ,0 + ,3355137 + ,112 + ,0 + ,3163748 + ,2591373 + ,2703997 + ,2713514 + ,2914972 + ,3651555 + ,0 + ,3613702 + ,113 + ,0 + ,3355137 + ,3163748 + ,2591373 + ,2703997 + ,2713514 + ,2914972 + ,0 + ,3686773 + ,114 + ,0 + ,3613702 + ,3355137 + ,3163748 + ,2591373 + ,2703997 + ,2713514 + ,0 + ,4098716 + ,115 + ,0 + ,3686773 + ,3613702 + ,3355137 + ,3163748 + ,2591373 + ,2703997 + ,0 + ,4063517 + ,116 + ,0 + ,4098716 + ,3686773 + ,3613702 + ,3355137 + ,3163748 + ,2591373 + ,1 + ,3551489 + ,117 + ,117 + ,4063517 + ,4098716 + ,3686773 + ,3613702 + ,3355137 + ,3163748 + ,1 + ,3226663 + ,118 + ,118 + ,3551489 + ,4063517 + ,4098716 + ,3686773 + ,3613702 + ,3355137 + ,1 + ,2656842 + ,119 + ,119 + ,3226663 + ,3551489 + ,4063517 + ,4098716 + ,3686773 + ,3613702 + ,1 + ,2597484 + ,120 + ,120 + ,2656842 + ,3226663 + ,3551489 + ,4063517 + ,4098716 + ,3686773 + ,1 + ,2572399 + ,121 + ,121 + ,2597484 + ,2656842 + ,3226663 + ,3551489 + ,4063517 + ,4098716 + ,1 + ,2596631 + ,122 + ,122 + ,2572399 + ,2597484 + ,2656842 + ,3226663 + ,3551489 + ,4063517 + ,1 + ,3165225 + ,123 + ,123 + ,2596631 + ,2572399 + ,2597484 + ,2656842 + ,3226663 + ,3551489 + ,1 + ,3303145 + ,124 + ,124 + ,3165225 + ,2596631 + ,2572399 + ,2597484 + ,2656842 + ,3226663 + ,1 + ,3698247 + ,125 + ,125 + ,3303145 + ,3165225 + ,2596631 + ,2572399 + ,2597484 + ,2656842 + ,1 + ,3668631 + ,126 + ,126 + ,3698247 + ,3303145 + ,3165225 + ,2596631 + ,2572399 + ,2597484 + ,1 + ,4130433 + ,127 + ,127 + ,3668631 + ,3698247 + ,3303145 + ,3165225 + ,2596631 + ,2572399 + ,1 + ,4131400 + ,128 + ,128 + ,4130433 + ,3668631 + ,3698247 + ,3303145 + ,3165225 + ,2596631 + ,1 + ,3864358 + ,129 + ,129 + ,4131400 + ,4130433 + ,3668631 + ,3698247 + ,3303145 + ,3165225 + ,1 + ,3721110 + ,130 + ,130 + ,3864358 + ,4131400 + ,4130433 + ,3668631 + ,3698247 + ,3303145 + ,1 + ,2892532 + ,131 + ,131 + ,3721110 + ,3864358 + ,4131400 + ,4130433 + ,3668631 + ,3698247 + ,1 + ,2843451 + ,132 + ,132 + ,2892532 + ,3721110 + ,3864358 + ,4131400 + ,4130433 + ,3668631 + ,1 + ,2747502 + ,133 + ,133 + ,2843451 + ,2892532 + ,3721110 + ,3864358 + ,4131400 + ,4130433 + ,1 + ,2668775 + ,134 + ,134 + ,2747502 + ,2843451 + ,2892532 + ,3721110 + ,3864358 + ,4131400 + ,1 + ,3018602 + ,135 + ,135 + ,2668775 + ,2747502 + ,2843451 + ,2892532 + ,3721110 + ,3864358 + ,1 + ,3013392 + ,136 + ,136 + ,3018602 + ,2668775 + ,2747502 + ,2843451 + ,2892532 + ,3721110 + ,1 + ,3393657 + ,137 + ,137 + ,3013392 + ,3018602 + ,2668775 + ,2747502 + ,2843451 + ,2892532 + ,1 + ,3544233 + ,138 + ,138 + ,3393657 + ,3013392 + ,3018602 + ,2668775 + ,2747502 + ,2843451 + ,1 + ,4075832 + ,139 + ,139 + ,3544233 + ,3393657 + ,3013392 + ,3018602 + ,2668775 + ,2747502 + ,1 + ,4032923 + ,140 + ,140 + ,4075832 + ,3544233 + ,3393657 + ,3013392 + ,3018602 + ,2668775 + ,1 + ,3734509 + ,141 + ,141 + ,4032923 + ,4075832 + ,3544233 + ,3393657 + ,3013392 + ,3018602 + ,1 + ,3761285 + ,142 + ,142 + ,3734509 + ,4032923 + ,4075832 + ,3544233 + ,3393657 + ,3013392 + ,1 + ,2970090 + ,143 + ,143 + ,3761285 + ,3734509 + ,4032923 + ,4075832 + ,3544233 + ,3393657 + ,1 + ,2847849 + ,144 + ,144 + ,2970090 + ,3761285 + ,3734509 + ,4032923 + ,4075832 + ,3544233 + ,1 + ,2741680 + ,145 + ,145 + ,2847849 + ,2970090 + ,3761285 + ,3734509 + ,4032923 + ,4075832 + ,1 + ,2830639 + ,146 + ,146 + ,2741680 + ,2847849 + ,2970090 + ,3761285 + ,3734509 + ,4032923 + ,1 + ,3257673 + ,147 + ,147 + ,2830639 + ,2741680 + ,2847849 + ,2970090 + ,3761285 + ,3734509 + ,1 + ,3480085 + ,148 + ,148 + ,3257673 + ,2830639 + ,2741680 + ,2847849 + ,2970090 + ,3761285 + ,1 + ,3843271 + ,149 + ,149 + ,3480085 + ,3257673 + ,2830639 + ,2741680 + ,2847849 + ,2970090 + ,1 + ,3796961 + ,150 + ,150 + ,3843271 + ,3480085 + ,3257673 + ,2830639 + ,2741680 + ,2847849 + ,1 + ,4337767 + ,151 + ,151 + ,3796961 + ,3843271 + ,3480085 + ,3257673 + ,2830639 + ,2741680 + ,1 + ,4243630 + ,152 + ,152 + ,4337767 + ,3796961 + ,3843271 + ,3480085 + ,3257673 + ,2830639 + ,1 + ,3927202 + ,153 + ,153 + ,4243630 + ,4337767 + ,3796961 + ,3843271 + ,3480085 + ,3257673 + ,1 + ,3915296 + ,154 + ,154 + ,3927202 + ,4243630 + ,4337767 + ,3796961 + ,3843271 + ,3480085 + ,1 + ,3087396 + ,155 + ,155 + ,3915296 + ,3927202 + ,4243630 + ,4337767 + ,3796961 + ,3843271 + ,1 + ,2963792 + ,156 + ,156 + ,3087396 + ,3915296 + ,3927202 + ,4243630 + ,4337767 + ,3796961 + ,1 + ,2955792 + ,157 + ,157 + ,2963792 + ,3087396 + ,3915296 + ,3927202 + ,4243630 + ,4337767 + ,1 + ,2829925 + ,158 + ,158 + ,2955792 + ,2963792 + ,3087396 + ,3915296 + ,3927202 + ,4243630 + ,1 + ,3281195 + ,159 + ,159 + ,2829925 + ,2955792 + ,2963792 + ,3087396 + ,3915296 + ,3927202 + ,1 + ,3548011 + ,160 + ,160 + ,3281195 + ,2829925 + ,2955792 + ,2963792 + ,3087396 + ,3915296 + ,1 + ,4059648 + ,161 + ,161 + ,3548011 + ,3281195 + ,2829925 + ,2955792 + ,2963792 + ,3087396 + ,1 + ,3941175 + ,162 + ,162 + ,4059648 + ,3548011 + ,3281195 + ,2829925 + ,2955792 + ,2963792 + ,1 + ,4528594 + ,163 + ,163 + ,3941175 + ,4059648 + ,3548011 + ,3281195 + ,2829925 + ,2955792 + ,1 + ,4433151 + ,164 + ,164 + ,4528594 + ,3941175 + ,4059648 + ,3548011 + ,3281195 + ,2829925 + ,1 + ,4145737 + ,165 + ,165 + ,4433151 + ,4528594 + ,3941175 + ,4059648 + ,3548011 + ,3281195 + ,1 + ,4077132 + ,166 + ,166 + ,4145737 + ,4433151 + ,4528594 + ,3941175 + ,4059648 + ,3548011 + ,1 + ,3198519 + ,167 + ,167 + ,4077132 + ,4145737 + ,4433151 + ,4528594 + ,3941175 + ,4059648 + ,1 + ,3078660 + ,168 + ,168 + ,3198519 + ,4077132 + ,4145737 + ,4433151 + ,4528594 + ,3941175 + ,1 + ,3028202 + ,169 + ,169 + ,3078660 + ,3198519 + ,4077132 + ,4145737 + ,4433151 + ,4528594 + ,1 + ,2858642 + ,170 + ,170 + ,3028202 + ,3078660 + ,3198519 + ,4077132 + ,4145737 + ,4433151 + ,1 + ,3398954 + ,171 + ,171 + ,2858642 + ,3028202 + ,3078660 + ,3198519 + ,4077132 + ,4145737 + ,1 + ,3808883 + ,172 + ,172 + ,3398954 + ,2858642 + ,3028202 + ,3078660 + ,3198519 + ,4077132 + ,1 + ,4175961 + ,173 + ,173 + ,3808883 + ,3398954 + ,2858642 + ,3028202 + ,3078660 + ,3198519 + ,1 + ,4227542 + ,174 + ,174 + ,4175961 + ,3808883 + ,3398954 + ,2858642 + ,3028202 + ,3078660 + ,1 + ,4744616 + ,175 + ,175 + ,4227542 + ,4175961 + ,3808883 + ,3398954 + ,2858642 + ,3028202 + ,1 + ,4608012 + ,176 + ,176 + ,4744616 + ,4227542 + ,4175961 + ,3808883 + ,3398954 + ,2858642 + ,1 + ,4295049 + ,177 + ,177 + ,4608012 + ,4744616 + ,4227542 + ,4175961 + ,3808883 + ,3398954 + ,1 + ,4201144 + ,178 + ,178 + ,4295049 + ,4608012 + ,4744616 + ,4227542 + ,4175961 + ,3808883 + ,1 + ,3353276 + ,179 + ,179 + ,4201144 + ,4295049 + ,4608012 + ,4744616 + ,4227542 + ,4175961 + ,1 + ,3286851 + ,180 + ,180 + ,3353276 + ,4201144 + ,4295049 + ,4608012 + ,4744616 + ,4227542 + ,1 + ,3169889 + ,181 + ,181 + ,3286851 + ,3353276 + ,4201144 + ,4295049 + ,4608012 + ,4744616 + ,1 + ,3051720 + ,182 + ,182 + ,3169889 + ,3286851 + ,3353276 + ,4201144 + ,4295049 + ,4608012 + ,1 + ,3695426 + ,183 + ,183 + ,3051720 + ,3169889 + ,3286851 + ,3353276 + ,4201144 + ,4295049 + ,1 + ,3905501 + ,184 + ,184 + ,3695426 + ,3051720 + ,3169889 + ,3286851 + ,3353276 + ,4201144 + ,1 + ,4296458 + ,185 + ,185 + ,3905501 + ,3695426 + ,3051720 + ,3169889 + ,3286851 + ,3353276 + ,1 + ,4246247 + ,186 + ,186 + ,4296458 + ,3905501 + ,3695426 + ,3051720 + ,3169889 + ,3286851 + ,1 + ,4921849 + ,187 + ,187 + ,4246247 + ,4296458 + ,3905501 + ,3695426 + ,3051720 + ,3169889 + ,1 + ,4821446 + ,188 + ,188 + ,4921849 + ,4246247 + ,4296458 + ,3905501 + ,3695426 + ,3051720 + ,1 + ,4425064 + ,189 + ,189 + ,4821446 + ,4921849 + ,4246247 + ,4296458 + ,3905501 + ,3695426 + ,1 + ,4379099 + ,190 + ,190 + ,4425064 + ,4821446 + ,4921849 + ,4246247 + ,4296458 + ,3905501 + ,1 + ,3472889 + ,191 + ,191 + ,4379099 + ,4425064 + ,4821446 + ,4921849 + ,4246247 + ,4296458 + ,1 + ,3359160 + ,192 + ,192 + ,3472889 + ,4379099 + ,4425064 + ,4821446 + ,4921849 + ,4246247 + ,1 + ,3200944 + ,193 + ,193 + ,3359160 + ,3472889 + ,4379099 + ,4425064 + ,4821446 + ,4921849 + ,1 + ,3153170 + ,194 + ,194 + ,3200944 + ,3359160 + ,3472889 + ,4379099 + ,4425064 + ,4821446 + ,1 + ,3741498 + ,195 + ,195 + ,3153170 + ,3200944 + ,3359160 + ,3472889 + ,4379099 + ,4425064 + ,1 + ,3918719 + ,196 + ,196 + ,3741498 + ,3153170 + ,3200944 + ,3359160 + ,3472889 + ,4379099 + ,1 + ,4403449 + ,197 + ,197 + ,3918719 + ,3741498 + ,3153170 + ,3200944 + ,3359160 + ,3472889 + ,1 + ,4400407 + ,198 + ,198 + ,4403449 + ,3918719 + ,3741498 + ,3153170 + ,3200944 + ,3359160 + ,1 + ,4847473 + ,199 + ,199 + ,4400407 + ,4403449 + ,3918719 + ,3741498 + ,3153170 + ,3200944 + ,1 + ,4716136 + ,200 + ,200 + ,4847473 + ,4400407 + ,4403449 + ,3918719 + ,3741498 + ,3153170 + ,1 + ,4297440 + ,201 + ,201 + ,4716136 + ,4847473 + ,4400407 + ,4403449 + ,3918719 + ,3741498 + ,1 + ,4272253 + ,202 + ,202 + ,4297440 + ,4716136 + ,4847473 + ,4400407 + ,4403449 + ,3918719 + ,1 + ,3271834 + ,203 + ,203 + ,4272253 + ,4297440 + ,4716136 + ,4847473 + ,4400407 + ,4403449 + ,1 + ,3168388 + ,204 + ,204 + ,3271834 + ,4272253 + ,4297440 + ,4716136 + ,4847473 + ,4400407 + ,1 + ,2911748 + ,205 + ,205 + ,3168388 + ,3271834 + ,4272253 + ,4297440 + ,4716136 + ,4847473 + ,1 + ,2720999 + ,206 + ,206 + ,2911748 + ,3168388 + ,3271834 + ,4272253 + ,4297440 + ,4716136 + ,1 + ,3199918 + ,207 + ,207 + ,2720999 + ,2911748 + ,3168388 + ,3271834 + ,4272253 + ,4297440 + ,1 + ,3672623 + ,208 + ,208 + ,3199918 + ,2720999 + ,2911748 + ,3168388 + ,3271834 + ,4272253 + ,1 + ,3892013 + ,209 + ,209 + ,3672623 + ,3199918 + ,2720999 + ,2911748 + ,3168388 + ,3271834 + ,1 + ,3850845 + ,210 + ,210 + ,3892013 + ,3672623 + ,3199918 + ,2720999 + ,2911748 + ,3168388 + ,1 + ,4532467 + ,211 + ,211 + ,3850845 + ,3892013 + ,3672623 + ,3199918 + ,2720999 + ,2911748 + ,1 + ,4484739 + ,212 + ,212 + ,4532467 + ,3850845 + ,3892013 + ,3672623 + ,3199918 + ,2720999 + ,1 + ,4014972 + ,213 + ,213 + ,4484739 + ,4532467 + ,3850845 + ,3892013 + ,3672623 + ,3199918 + ,1 + ,3983758 + ,214 + ,214 + ,4014972 + ,4484739 + ,4532467 + ,3850845 + ,3892013 + ,3672623 + ,1 + ,3158459 + ,215 + ,215 + ,3983758 + ,4014972 + ,4484739 + ,4532467 + ,3850845 + ,3892013 + ,1 + ,3100569 + ,216 + ,216 + ,3158459 + ,3983758 + ,4014972 + ,4484739 + ,4532467 + ,3850845 + ,1 + ,2935404 + ,217 + ,217 + ,3100569 + ,3158459 + ,3983758 + ,4014972 + ,4484739 + ,4532467 + ,1 + ,2855719 + ,218 + ,218 + ,2935404 + ,3100569 + ,3158459 + ,3983758 + ,4014972 + ,4484739 + ,1 + ,3465611 + ,219 + ,219 + ,2855719 + ,2935404 + ,3100569 + ,3158459 + ,3983758 + ,4014972 + ,1 + ,3006985 + ,220 + ,220 + ,3465611 + ,2855719 + ,2935404 + ,3100569 + ,3158459 + ,3983758 + ,1 + ,4095110 + ,221 + ,221 + ,3006985 + ,3465611 + ,2855719 + ,2935404 + ,3100569 + ,3158459 + ,1 + ,4104793 + ,222 + ,222 + ,4095110 + ,3006985 + ,3465611 + ,2855719 + ,2935404 + ,3100569 + ,1 + ,4730788 + ,223 + ,223 + ,4104793 + ,4095110 + ,3006985 + ,3465611 + ,2855719 + ,2935404 + ,1 + ,4642726 + ,224 + ,224 + ,4730788 + ,4104793 + ,4095110 + ,3006985 + ,3465611 + ,2855719 + ,1 + ,4246919 + ,225 + ,225 + ,4642726 + ,4730788 + ,4104793 + ,4095110 + ,3006985 + ,3465611 + ,1 + ,4308117 + ,226 + ,226 + ,4246919 + ,4642726 + ,4730788 + ,4104793 + ,4095110 + ,3006985) + ,dim=c(10 + ,220) + ,dimnames=list(c('9/11' + ,'Yt' + ,'t' + ,'T_9/11' + ,'Yt-1' + ,'Yt-2' + ,'Yt-3' + ,'Yt-4' + ,'Yt-5' + ,'Yt-6') + ,1:220)) > y <- array(NA,dim=c(10,220),dimnames=list(c('9/11','Yt','t','T_9/11','Yt-1','Yt-2','Yt-3','Yt-4','Yt-5','Yt-6'),1:220)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par1 = 'kendall' > 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/1r1pw1292020349.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) > 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] 10 > 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/2c1521292020349.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) Kendall's rank correlation tau data: y[1, ] and y[2, ] z = 9.8292, p-value < 2.2e-16 alternative hypothesis: true tau is not equal to 0 sample estimates: tau 0.5435467 > 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 17 warnings (use warnings() to see them) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/3qb3a1292020349.tab") > try(system("convert tmp/1r1pw1292020349.ps tmp/1r1pw1292020349.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.587 1.289 7.444