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(0 + ,1281151 + ,13 + ,0 + ,1281814 + ,1298756 + ,1707752 + ,1867943 + ,0 + ,1164976 + ,14 + ,0 + ,1281151 + ,1281814 + ,1298756 + ,1707752 + ,0 + ,1454329 + ,15 + ,0 + ,1164976 + ,1281151 + ,1281814 + ,1298756 + ,0 + ,1645288 + ,16 + ,0 + ,1454329 + ,1164976 + ,1281151 + ,1281814 + ,0 + ,1817743 + ,17 + ,0 + ,1645288 + ,1454329 + ,1164976 + ,1281151 + ,0 + ,1895785 + ,18 + ,0 + ,1817743 + ,1645288 + ,1454329 + ,1164976 + ,0 + ,2236311 + ,19 + ,0 + ,1895785 + ,1817743 + ,1645288 + ,1454329 + ,0 + ,2295951 + ,20 + ,0 + ,2236311 + ,1895785 + ,1817743 + ,1645288 + ,0 + ,2087315 + ,21 + ,0 + ,2295951 + ,2236311 + ,1895785 + ,1817743 + ,0 + ,1980891 + ,22 + ,0 + ,2087315 + ,2295951 + ,2236311 + ,1895785 + ,0 + ,1465446 + ,23 + ,0 + ,1980891 + ,2087315 + ,2295951 + ,2236311 + ,0 + ,1445026 + ,24 + ,0 + ,1465446 + ,1980891 + ,2087315 + ,2295951 + ,0 + ,1488120 + ,25 + ,0 + ,1445026 + ,1465446 + ,1980891 + ,2087315 + ,0 + ,1338333 + ,26 + ,0 + ,1488120 + ,1445026 + ,1465446 + ,1980891 + ,0 + ,1715789 + ,27 + ,0 + ,1338333 + ,1488120 + ,1445026 + ,1465446 + ,0 + ,1806090 + ,28 + ,0 + ,1715789 + ,1338333 + ,1488120 + ,1445026 + ,0 + ,2083316 + ,29 + ,0 + ,1806090 + ,1715789 + ,1338333 + ,1488120 + ,0 + ,2092278 + ,30 + ,0 + ,2083316 + ,1806090 + ,1715789 + ,1338333 + ,0 + ,2430800 + ,31 + ,0 + ,2092278 + ,2083316 + ,1806090 + ,1715789 + ,0 + ,2424894 + ,32 + ,0 + ,2430800 + ,2092278 + ,2083316 + ,1806090 + ,0 + ,2299016 + ,33 + ,0 + ,2424894 + ,2430800 + ,2092278 + ,2083316 + ,0 + ,2130688 + ,34 + ,0 + ,2299016 + ,2424894 + ,2430800 + ,2092278 + ,0 + ,1652221 + ,35 + ,0 + ,2130688 + ,2299016 + ,2424894 + ,2430800 + ,0 + ,1608162 + ,36 + ,0 + ,1652221 + ,2130688 + ,2299016 + ,2424894 + ,0 + ,1647074 + ,37 + ,0 + ,1608162 + ,1652221 + ,2130688 + ,2299016 + ,0 + ,1479691 + ,38 + ,0 + ,1647074 + ,1608162 + ,1652221 + ,2130688 + ,0 + ,1884978 + ,39 + ,0 + ,1479691 + ,1647074 + ,1608162 + ,1652221 + ,0 + ,2007898 + ,40 + ,0 + ,1884978 + ,1479691 + ,1647074 + ,1608162 + ,0 + ,2208954 + ,41 + ,0 + ,2007898 + ,1884978 + ,1479691 + ,1647074 + ,0 + ,2217164 + ,42 + ,0 + ,2208954 + ,2007898 + ,1884978 + ,1479691 + ,0 + ,2534291 + ,43 + ,0 + ,2217164 + ,2208954 + ,2007898 + ,1884978 + ,0 + ,2560312 + ,44 + ,0 + ,2534291 + ,2217164 + ,2208954 + ,2007898 + ,0 + ,2429069 + ,45 + ,0 + ,2560312 + ,2534291 + ,2217164 + ,2208954 + ,0 + ,2315077 + ,46 + ,0 + ,2429069 + ,2560312 + ,2534291 + ,2217164 + ,0 + ,1799608 + ,47 + ,0 + ,2315077 + ,2429069 + ,2560312 + ,2534291 + ,0 + ,1772590 + ,48 + ,0 + ,1799608 + ,2315077 + ,2429069 + ,2560312 + ,0 + ,1744799 + ,49 + ,0 + ,1772590 + ,1799608 + ,2315077 + ,2429069 + ,0 + ,1659093 + ,50 + ,0 + ,1744799 + ,1772590 + ,1799608 + ,2315077 + ,0 + ,2099821 + ,51 + ,0 + ,1659093 + ,1744799 + ,1772590 + ,1799608 + ,0 + ,2135736 + ,52 + ,0 + ,2099821 + ,1659093 + ,1744799 + ,1772590 + ,0 + ,2427894 + ,53 + ,0 + ,2135736 + ,2099821 + ,1659093 + ,1744799 + ,0 + ,2468882 + ,54 + ,0 + ,2427894 + ,2135736 + ,2099821 + ,1659093 + ,0 + ,2703217 + ,55 + ,0 + ,2468882 + ,2427894 + ,2135736 + ,2099821 + ,0 + ,2766841 + ,56 + ,0 + ,2703217 + ,2468882 + ,2427894 + ,2135736 + ,0 + ,2655236 + ,57 + ,0 + ,2766841 + ,2703217 + ,2468882 + ,2427894 + ,0 + ,2550373 + ,58 + ,0 + ,2655236 + ,2766841 + ,2703217 + ,2468882 + ,0 + ,2052097 + ,59 + ,0 + ,2550373 + ,2655236 + ,2766841 + ,2703217 + ,0 + ,1998055 + ,60 + ,0 + ,2052097 + ,2550373 + ,2655236 + ,2766841 + ,0 + ,1920748 + ,61 + ,0 + ,1998055 + ,2052097 + ,2550373 + ,2655236 + ,0 + ,1876694 + ,62 + ,0 + ,1920748 + ,1998055 + ,2052097 + ,2550373 + ,0 + ,2380930 + ,63 + ,0 + ,1876694 + ,1920748 + ,1998055 + ,2052097 + ,0 + ,2467402 + ,64 + ,0 + ,2380930 + ,1876694 + ,1920748 + ,1998055 + ,0 + ,2770771 + ,65 + ,0 + ,2467402 + ,2380930 + ,1876694 + ,1920748 + ,0 + ,2781340 + ,66 + ,0 + ,2770771 + ,2467402 + ,2380930 + ,1876694 + ,0 + ,3143926 + ,67 + ,0 + ,2781340 + ,2770771 + ,2467402 + ,2380930 + ,0 + ,3172235 + ,68 + ,0 + ,3143926 + ,2781340 + ,2770771 + ,2467402 + ,0 + ,2952540 + ,69 + ,0 + ,3172235 + ,3143926 + ,2781340 + ,2770771 + ,0 + ,2920877 + ,70 + ,0 + ,2952540 + ,3172235 + ,3143926 + ,2781340 + ,0 + ,2384552 + ,71 + ,0 + ,2920877 + ,2952540 + ,3172235 + ,3143926 + ,0 + ,2248987 + ,72 + ,0 + ,2384552 + ,2920877 + ,2952540 + ,3172235 + ,0 + ,2208616 + ,73 + ,0 + ,2248987 + ,2384552 + ,2920877 + ,2952540 + ,0 + ,2178756 + ,74 + ,0 + ,2208616 + ,2248987 + ,2384552 + ,2920877 + ,0 + ,2632870 + ,75 + ,0 + ,2178756 + ,2208616 + ,2248987 + ,2384552 + ,0 + ,2706905 + ,76 + ,0 + ,2632870 + ,2178756 + ,2208616 + ,2248987 + ,0 + ,3029745 + ,77 + ,0 + ,2706905 + ,2632870 + ,2178756 + ,2208616 + ,0 + ,3015402 + ,78 + ,0 + ,3029745 + ,2706905 + ,2632870 + ,2178756 + ,0 + ,3391414 + ,79 + ,0 + ,3015402 + ,3029745 + ,2706905 + ,2632870 + ,0 + ,3507805 + ,80 + ,0 + ,3391414 + ,3015402 + ,3029745 + ,2706905 + ,0 + ,3177852 + ,81 + ,0 + ,3507805 + ,3391414 + ,3015402 + ,3029745 + ,0 + ,3142961 + ,82 + ,0 + ,3177852 + ,3507805 + ,3391414 + ,3015402 + ,0 + ,2545815 + ,83 + ,0 + ,3142961 + ,3177852 + ,3507805 + ,3391414 + ,0 + ,2414007 + ,84 + ,0 + ,2545815 + ,3142961 + ,3177852 + ,3507805 + ,0 + ,2372578 + ,85 + ,0 + ,2414007 + ,2545815 + ,3142961 + ,3177852 + ,0 + ,2332664 + ,86 + ,0 + ,2372578 + ,2414007 + ,2545815 + ,3142961 + ,0 + ,2825328 + ,87 + ,0 + ,2332664 + ,2372578 + ,2414007 + ,2545815 + ,0 + ,2901478 + ,88 + ,0 + ,2825328 + ,2332664 + ,2372578 + ,2414007 + ,0 + ,3263955 + ,89 + ,0 + ,2901478 + ,2825328 + ,2332664 + ,2372578 + ,0 + ,3226738 + ,90 + ,0 + ,3263955 + ,2901478 + ,2825328 + ,2332664 + ,0 + ,3610786 + ,91 + ,0 + ,3226738 + ,3263955 + ,2901478 + ,2825328 + ,0 + ,3709274 + ,92 + ,0 + ,3610786 + ,3226738 + ,3263955 + ,2901478 + ,0 + ,3467185 + ,93 + ,0 + ,3709274 + ,3610786 + ,3226738 + ,3263955 + ,0 + ,3449646 + ,94 + ,0 + ,3467185 + ,3709274 + ,3610786 + ,3226738 + ,0 + ,2802951 + ,95 + ,0 + ,3449646 + ,3467185 + ,3709274 + ,3610786 + ,0 + ,2462530 + ,96 + ,0 + ,2802951 + ,3449646 + ,3467185 + ,3709274 + ,0 + ,2490645 + ,97 + ,0 + ,2462530 + ,2802951 + ,3449646 + ,3467185 + ,0 + ,2561520 + ,98 + ,0 + ,2490645 + ,2462530 + ,2802951 + ,3449646 + ,0 + ,3067554 + ,99 + ,0 + ,2561520 + ,2490645 + ,2462530 + ,2802951 + ,0 + ,3226951 + ,100 + ,0 + ,3067554 + ,2561520 + ,2490645 + ,2462530 + ,0 + ,3546493 + ,101 + ,0 + ,3226951 + ,3067554 + ,2561520 + ,2490645 + ,0 + ,3492787 + ,102 + ,0 + ,3546493 + ,3226951 + ,3067554 + ,2561520 + ,0 + ,3952263 + ,103 + ,0 + ,3492787 + ,3546493 + ,3226951 + ,3067554 + ,0 + ,3932072 + ,104 + ,0 + ,3952263 + ,3492787 + ,3546493 + ,3226951 + ,0 + ,3720284 + ,105 + ,0 + ,3932072 + ,3952263 + ,3492787 + ,3546493 + ,0 + ,3651555 + ,106 + ,0 + ,3720284 + ,3932072 + ,3952263 + ,3492787 + ,0 + ,2914972 + ,107 + ,0 + ,3651555 + ,3720284 + ,3932072 + ,3952263 + ,0 + ,2713514 + ,108 + ,0 + ,2914972 + ,3651555 + ,3720284 + ,3932072 + ,0 + ,2703997 + ,109 + ,0 + ,2713514 + ,2914972 + ,3651555 + ,3720284 + ,0 + ,2591373 + ,110 + ,0 + ,2703997 + ,2713514 + ,2914972 + ,3651555 + ,0 + ,3163748 + ,111 + ,0 + ,2591373 + ,2703997 + ,2713514 + ,2914972 + ,0 + ,3355137 + ,112 + ,0 + ,3163748 + ,2591373 + ,2703997 + ,2713514 + ,0 + ,3613702 + ,113 + ,0 + ,3355137 + ,3163748 + ,2591373 + ,2703997 + ,0 + ,3686773 + ,114 + ,0 + ,3613702 + ,3355137 + ,3163748 + ,2591373 + ,0 + ,4098716 + ,115 + ,0 + ,3686773 + ,3613702 + ,3355137 + ,3163748 + ,0 + ,4063517 + ,116 + ,0 + ,4098716 + ,3686773 + ,3613702 + ,3355137 + ,1 + ,3551489 + ,117 + ,117 + ,4063517 + ,4098716 + ,3686773 + ,3613702 + ,1 + ,3226663 + ,118 + ,118 + ,3551489 + ,4063517 + ,4098716 + ,3686773 + ,1 + ,2656842 + ,119 + ,119 + ,3226663 + ,3551489 + ,4063517 + ,4098716 + ,1 + ,2597484 + ,120 + ,120 + ,2656842 + ,3226663 + ,3551489 + ,4063517 + ,1 + ,2572399 + ,121 + ,121 + ,2597484 + ,2656842 + ,3226663 + ,3551489 + ,1 + ,2596631 + ,122 + ,122 + ,2572399 + ,2597484 + ,2656842 + ,3226663 + ,1 + ,3165225 + ,123 + ,123 + ,2596631 + ,2572399 + ,2597484 + ,2656842 + ,1 + ,3303145 + ,124 + ,124 + ,3165225 + ,2596631 + ,2572399 + ,2597484 + ,1 + ,3698247 + ,125 + ,125 + ,3303145 + ,3165225 + ,2596631 + ,2572399 + ,1 + ,3668631 + ,126 + ,126 + ,3698247 + ,3303145 + ,3165225 + ,2596631 + ,1 + ,4130433 + ,127 + ,127 + ,3668631 + ,3698247 + ,3303145 + ,3165225 + ,1 + ,4131400 + ,128 + ,128 + ,4130433 + ,3668631 + ,3698247 + ,3303145 + ,1 + ,3864358 + ,129 + ,129 + ,4131400 + ,4130433 + ,3668631 + ,3698247 + ,1 + ,3721110 + ,130 + ,130 + ,3864358 + ,4131400 + ,4130433 + ,3668631 + ,1 + ,2892532 + ,131 + ,131 + ,3721110 + ,3864358 + ,4131400 + ,4130433 + ,1 + ,2843451 + ,132 + ,132 + ,2892532 + ,3721110 + ,3864358 + ,4131400 + ,1 + ,2747502 + ,133 + ,133 + ,2843451 + ,2892532 + ,3721110 + ,3864358 + ,1 + ,2668775 + ,134 + ,134 + ,2747502 + ,2843451 + ,2892532 + ,3721110 + ,1 + ,3018602 + ,135 + ,135 + ,2668775 + ,2747502 + ,2843451 + ,2892532 + ,1 + ,3013392 + ,136 + ,136 + ,3018602 + ,2668775 + ,2747502 + ,2843451 + ,1 + ,3393657 + ,137 + ,137 + ,3013392 + ,3018602 + ,2668775 + ,2747502 + ,1 + ,3544233 + ,138 + ,138 + ,3393657 + ,3013392 + ,3018602 + ,2668775 + ,1 + ,4075832 + ,139 + ,139 + ,3544233 + ,3393657 + ,3013392 + ,3018602 + ,1 + ,4032923 + ,140 + ,140 + ,4075832 + ,3544233 + ,3393657 + ,3013392 + ,1 + ,3734509 + ,141 + ,141 + ,4032923 + ,4075832 + ,3544233 + ,3393657 + ,1 + ,3761285 + ,142 + ,142 + ,3734509 + ,4032923 + ,4075832 + ,3544233 + ,1 + ,2970090 + ,143 + ,143 + ,3761285 + ,3734509 + ,4032923 + ,4075832 + ,1 + ,2847849 + ,144 + ,144 + ,2970090 + ,3761285 + ,3734509 + ,4032923 + ,1 + ,2741680 + ,145 + ,145 + ,2847849 + ,2970090 + ,3761285 + ,3734509 + ,1 + ,2830639 + ,146 + ,146 + ,2741680 + ,2847849 + ,2970090 + ,3761285 + ,1 + ,3257673 + ,147 + ,147 + ,2830639 + ,2741680 + ,2847849 + ,2970090 + ,1 + ,3480085 + ,148 + ,148 + ,3257673 + ,2830639 + ,2741680 + ,2847849 + ,1 + ,3843271 + ,149 + ,149 + ,3480085 + ,3257673 + ,2830639 + ,2741680 + ,1 + ,3796961 + ,150 + ,150 + ,3843271 + ,3480085 + ,3257673 + ,2830639 + ,1 + ,4337767 + ,151 + ,151 + ,3796961 + ,3843271 + ,3480085 + ,3257673 + ,1 + ,4243630 + ,152 + ,152 + ,4337767 + ,3796961 + ,3843271 + ,3480085 + ,1 + ,3927202 + ,153 + ,153 + ,4243630 + ,4337767 + ,3796961 + ,3843271 + ,1 + ,3915296 + ,154 + ,154 + ,3927202 + ,4243630 + ,4337767 + ,3796961 + ,1 + ,3087396 + ,155 + ,155 + ,3915296 + ,3927202 + ,4243630 + ,4337767 + ,1 + ,2963792 + ,156 + ,156 + ,3087396 + ,3915296 + ,3927202 + ,4243630 + ,1 + ,2955792 + ,157 + ,157 + ,2963792 + ,3087396 + ,3915296 + ,3927202 + ,1 + ,2829925 + ,158 + ,158 + ,2955792 + ,2963792 + ,3087396 + ,3915296 + ,1 + ,3281195 + ,159 + ,159 + ,2829925 + ,2955792 + ,2963792 + ,3087396 + ,1 + ,3548011 + ,160 + ,160 + ,3281195 + ,2829925 + ,2955792 + ,2963792 + ,1 + ,4059648 + ,161 + ,161 + ,3548011 + ,3281195 + ,2829925 + ,2955792 + ,1 + ,3941175 + ,162 + ,162 + ,4059648 + ,3548011 + ,3281195 + ,2829925 + ,1 + ,4528594 + ,163 + ,163 + ,3941175 + ,4059648 + ,3548011 + ,3281195 + ,1 + ,4433151 + ,164 + ,164 + ,4528594 + ,3941175 + ,4059648 + ,3548011 + ,1 + ,4145737 + ,165 + ,165 + ,4433151 + ,4528594 + ,3941175 + ,4059648 + ,1 + ,4077132 + ,166 + ,166 + ,4145737 + ,4433151 + ,4528594 + ,3941175 + ,1 + ,3198519 + ,167 + ,167 + ,4077132 + ,4145737 + ,4433151 + ,4528594 + ,1 + ,3078660 + ,168 + ,168 + ,3198519 + ,4077132 + ,4145737 + ,4433151 + ,1 + ,3028202 + ,169 + ,169 + ,3078660 + ,3198519 + ,4077132 + ,4145737 + ,1 + ,2858642 + ,170 + ,170 + ,3028202 + ,3078660 + ,3198519 + ,4077132 + ,1 + ,3398954 + ,171 + ,171 + ,2858642 + ,3028202 + ,3078660 + ,3198519 + ,1 + ,3808883 + ,172 + ,172 + ,3398954 + ,2858642 + ,3028202 + ,3078660 + ,1 + ,4175961 + ,173 + ,173 + ,3808883 + ,3398954 + ,2858642 + ,3028202 + ,1 + ,4227542 + ,174 + ,174 + ,4175961 + ,3808883 + ,3398954 + ,2858642 + ,1 + ,4744616 + ,175 + ,175 + ,4227542 + ,4175961 + ,3808883 + ,3398954 + ,1 + ,4608012 + ,176 + ,176 + ,4744616 + ,4227542 + ,4175961 + ,3808883 + ,1 + ,4295049 + ,177 + ,177 + ,4608012 + ,4744616 + ,4227542 + ,4175961 + ,1 + ,4201144 + ,178 + ,178 + ,4295049 + ,4608012 + ,4744616 + ,4227542 + ,1 + ,3353276 + ,179 + ,179 + ,4201144 + ,4295049 + ,4608012 + ,4744616 + ,1 + ,3286851 + ,180 + ,180 + ,3353276 + ,4201144 + ,4295049 + ,4608012 + ,1 + ,3169889 + ,181 + ,181 + ,3286851 + ,3353276 + ,4201144 + ,4295049 + ,1 + ,3051720 + ,182 + ,182 + ,3169889 + ,3286851 + ,3353276 + ,4201144 + ,1 + ,3695426 + ,183 + ,183 + ,3051720 + ,3169889 + ,3286851 + ,3353276 + ,1 + ,3905501 + ,184 + ,184 + ,3695426 + ,3051720 + ,3169889 + ,3286851 + ,1 + ,4296458 + ,185 + ,185 + ,3905501 + ,3695426 + ,3051720 + ,3169889 + ,1 + ,4246247 + ,186 + ,186 + ,4296458 + ,3905501 + ,3695426 + ,3051720 + ,1 + ,4921849 + ,187 + ,187 + ,4246247 + ,4296458 + ,3905501 + ,3695426 + ,1 + ,4821446 + ,188 + ,188 + ,4921849 + ,4246247 + ,4296458 + ,3905501 + ,1 + ,4425064 + ,189 + ,189 + ,4821446 + ,4921849 + ,4246247 + ,4296458 + ,1 + ,4379099 + ,190 + ,190 + ,4425064 + ,4821446 + ,4921849 + ,4246247 + ,1 + ,3472889 + ,191 + ,191 + ,4379099 + ,4425064 + ,4821446 + ,4921849 + ,1 + ,3359160 + ,192 + ,192 + ,3472889 + ,4379099 + ,4425064 + ,4821446 + ,1 + ,3200944 + ,193 + ,193 + ,3359160 + ,3472889 + ,4379099 + ,4425064 + ,1 + ,3153170 + ,194 + ,194 + ,3200944 + ,3359160 + ,3472889 + ,4379099 + ,1 + ,3741498 + ,195 + ,195 + ,3153170 + ,3200944 + ,3359160 + ,3472889 + ,1 + ,3918719 + ,196 + ,196 + ,3741498 + ,3153170 + ,3200944 + ,3359160 + ,1 + ,4403449 + ,197 + ,197 + ,3918719 + ,3741498 + ,3153170 + ,3200944 + ,1 + ,4400407 + ,198 + ,198 + ,4403449 + ,3918719 + ,3741498 + ,3153170 + ,1 + ,4847473 + ,199 + ,199 + ,4400407 + ,4403449 + ,3918719 + ,3741498 + ,1 + ,4716136 + ,200 + ,200 + ,4847473 + ,4400407 + ,4403449 + ,3918719 + ,1 + ,4297440 + ,201 + ,201 + ,4716136 + ,4847473 + ,4400407 + ,4403449 + ,1 + ,4272253 + ,202 + ,202 + ,4297440 + ,4716136 + ,4847473 + ,4400407 + ,1 + ,3271834 + ,203 + ,203 + ,4272253 + ,4297440 + ,4716136 + ,4847473 + ,1 + ,3168388 + ,204 + ,204 + ,3271834 + ,4272253 + ,4297440 + ,4716136 + ,1 + ,2911748 + ,205 + ,205 + ,3168388 + ,3271834 + ,4272253 + ,4297440 + ,1 + ,2720999 + ,206 + ,206 + ,2911748 + ,3168388 + ,3271834 + ,4272253 + ,1 + ,3199918 + ,207 + ,207 + ,2720999 + ,2911748 + ,3168388 + ,3271834 + ,1 + ,3672623 + ,208 + ,208 + ,3199918 + ,2720999 + ,2911748 + ,3168388 + ,1 + ,3892013 + ,209 + ,209 + ,3672623 + ,3199918 + ,2720999 + ,2911748 + ,1 + ,3850845 + ,210 + ,210 + ,3892013 + ,3672623 + ,3199918 + ,2720999 + ,1 + ,4532467 + ,211 + ,211 + ,3850845 + ,3892013 + ,3672623 + ,3199918 + ,1 + ,4484739 + ,212 + ,212 + ,4532467 + ,3850845 + ,3892013 + ,3672623 + ,1 + ,4014972 + ,213 + ,213 + ,4484739 + ,4532467 + ,3850845 + ,3892013 + ,1 + ,3983758 + ,214 + ,214 + ,4014972 + ,4484739 + ,4532467 + ,3850845 + ,1 + ,3158459 + ,215 + ,215 + ,3983758 + ,4014972 + ,4484739 + ,4532467 + ,1 + ,3100569 + ,216 + ,216 + ,3158459 + ,3983758 + ,4014972 + ,4484739 + ,1 + ,2935404 + ,217 + ,217 + ,3100569 + ,3158459 + ,3983758 + ,4014972 + ,1 + ,2855719 + ,218 + ,218 + ,2935404 + ,3100569 + ,3158459 + ,3983758 + ,1 + ,3465611 + ,219 + ,219 + ,2855719 + ,2935404 + ,3100569 + ,3158459 + ,1 + ,3006985 + ,220 + ,220 + ,3465611 + ,2855719 + ,2935404 + ,3100569 + ,1 + ,4095110 + ,221 + ,221 + ,3006985 + ,3465611 + ,2855719 + ,2935404 + ,1 + ,4104793 + ,222 + ,222 + ,4095110 + ,3006985 + ,3465611 + ,2855719 + ,1 + ,4730788 + ,223 + ,223 + ,4104793 + ,4095110 + ,3006985 + ,3465611 + ,1 + ,4642726 + ,224 + ,224 + ,4730788 + ,4104793 + ,4095110 + ,3006985 + ,1 + ,4246919 + ,225 + ,225 + ,4642726 + ,4730788 + ,4104793 + ,4095110 + ,1 + ,4308117 + ,226 + ,226 + ,4246919 + ,4642726 + ,4730788 + ,4104793) + ,dim=c(8 + ,214) + ,dimnames=list(c('9/11' + ,'Passagiers' + ,'t' + ,'T_9/11' + ,'Passagiers-1' + ,'Passagiers-2' + ,'Passagiers-3' + ,'Passagiers-4') + ,1:214)) > y <- array(NA,dim=c(8,214),dimnames=list(c('9/11','Passagiers','t','T_9/11','Passagiers-1','Passagiers-2','Passagiers-3','Passagiers-4'),1:214)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Include Monthly Dummies' > par1 = '2' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), 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: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Passagiers 9/11 t T_9/11 Passagiers-1 Passagiers-2 Passagiers-3 1 1281151 0 13 0 1281814 1298756 1707752 2 1164976 0 14 0 1281151 1281814 1298756 3 1454329 0 15 0 1164976 1281151 1281814 4 1645288 0 16 0 1454329 1164976 1281151 5 1817743 0 17 0 1645288 1454329 1164976 6 1895785 0 18 0 1817743 1645288 1454329 7 2236311 0 19 0 1895785 1817743 1645288 8 2295951 0 20 0 2236311 1895785 1817743 9 2087315 0 21 0 2295951 2236311 1895785 10 1980891 0 22 0 2087315 2295951 2236311 11 1465446 0 23 0 1980891 2087315 2295951 12 1445026 0 24 0 1465446 1980891 2087315 13 1488120 0 25 0 1445026 1465446 1980891 14 1338333 0 26 0 1488120 1445026 1465446 15 1715789 0 27 0 1338333 1488120 1445026 16 1806090 0 28 0 1715789 1338333 1488120 17 2083316 0 29 0 1806090 1715789 1338333 18 2092278 0 30 0 2083316 1806090 1715789 19 2430800 0 31 0 2092278 2083316 1806090 20 2424894 0 32 0 2430800 2092278 2083316 21 2299016 0 33 0 2424894 2430800 2092278 22 2130688 0 34 0 2299016 2424894 2430800 23 1652221 0 35 0 2130688 2299016 2424894 24 1608162 0 36 0 1652221 2130688 2299016 25 1647074 0 37 0 1608162 1652221 2130688 26 1479691 0 38 0 1647074 1608162 1652221 27 1884978 0 39 0 1479691 1647074 1608162 28 2007898 0 40 0 1884978 1479691 1647074 29 2208954 0 41 0 2007898 1884978 1479691 30 2217164 0 42 0 2208954 2007898 1884978 31 2534291 0 43 0 2217164 2208954 2007898 32 2560312 0 44 0 2534291 2217164 2208954 33 2429069 0 45 0 2560312 2534291 2217164 34 2315077 0 46 0 2429069 2560312 2534291 35 1799608 0 47 0 2315077 2429069 2560312 36 1772590 0 48 0 1799608 2315077 2429069 37 1744799 0 49 0 1772590 1799608 2315077 38 1659093 0 50 0 1744799 1772590 1799608 39 2099821 0 51 0 1659093 1744799 1772590 40 2135736 0 52 0 2099821 1659093 1744799 41 2427894 0 53 0 2135736 2099821 1659093 42 2468882 0 54 0 2427894 2135736 2099821 43 2703217 0 55 0 2468882 2427894 2135736 44 2766841 0 56 0 2703217 2468882 2427894 45 2655236 0 57 0 2766841 2703217 2468882 46 2550373 0 58 0 2655236 2766841 2703217 47 2052097 0 59 0 2550373 2655236 2766841 48 1998055 0 60 0 2052097 2550373 2655236 49 1920748 0 61 0 1998055 2052097 2550373 50 1876694 0 62 0 1920748 1998055 2052097 51 2380930 0 63 0 1876694 1920748 1998055 52 2467402 0 64 0 2380930 1876694 1920748 53 2770771 0 65 0 2467402 2380930 1876694 54 2781340 0 66 0 2770771 2467402 2380930 55 3143926 0 67 0 2781340 2770771 2467402 56 3172235 0 68 0 3143926 2781340 2770771 57 2952540 0 69 0 3172235 3143926 2781340 58 2920877 0 70 0 2952540 3172235 3143926 59 2384552 0 71 0 2920877 2952540 3172235 60 2248987 0 72 0 2384552 2920877 2952540 61 2208616 0 73 0 2248987 2384552 2920877 62 2178756 0 74 0 2208616 2248987 2384552 63 2632870 0 75 0 2178756 2208616 2248987 64 2706905 0 76 0 2632870 2178756 2208616 65 3029745 0 77 0 2706905 2632870 2178756 66 3015402 0 78 0 3029745 2706905 2632870 67 3391414 0 79 0 3015402 3029745 2706905 68 3507805 0 80 0 3391414 3015402 3029745 69 3177852 0 81 0 3507805 3391414 3015402 70 3142961 0 82 0 3177852 3507805 3391414 71 2545815 0 83 0 3142961 3177852 3507805 72 2414007 0 84 0 2545815 3142961 3177852 73 2372578 0 85 0 2414007 2545815 3142961 74 2332664 0 86 0 2372578 2414007 2545815 75 2825328 0 87 0 2332664 2372578 2414007 76 2901478 0 88 0 2825328 2332664 2372578 77 3263955 0 89 0 2901478 2825328 2332664 78 3226738 0 90 0 3263955 2901478 2825328 79 3610786 0 91 0 3226738 3263955 2901478 80 3709274 0 92 0 3610786 3226738 3263955 81 3467185 0 93 0 3709274 3610786 3226738 82 3449646 0 94 0 3467185 3709274 3610786 83 2802951 0 95 0 3449646 3467185 3709274 84 2462530 0 96 0 2802951 3449646 3467185 85 2490645 0 97 0 2462530 2802951 3449646 86 2561520 0 98 0 2490645 2462530 2802951 87 3067554 0 99 0 2561520 2490645 2462530 88 3226951 0 100 0 3067554 2561520 2490645 89 3546493 0 101 0 3226951 3067554 2561520 90 3492787 0 102 0 3546493 3226951 3067554 91 3952263 0 103 0 3492787 3546493 3226951 92 3932072 0 104 0 3952263 3492787 3546493 93 3720284 0 105 0 3932072 3952263 3492787 94 3651555 0 106 0 3720284 3932072 3952263 95 2914972 0 107 0 3651555 3720284 3932072 96 2713514 0 108 0 2914972 3651555 3720284 97 2703997 0 109 0 2713514 2914972 3651555 98 2591373 0 110 0 2703997 2713514 2914972 99 3163748 0 111 0 2591373 2703997 2713514 100 3355137 0 112 0 3163748 2591373 2703997 101 3613702 0 113 0 3355137 3163748 2591373 102 3686773 0 114 0 3613702 3355137 3163748 103 4098716 0 115 0 3686773 3613702 3355137 104 4063517 0 116 0 4098716 3686773 3613702 105 3551489 1 117 117 4063517 4098716 3686773 106 3226663 1 118 118 3551489 4063517 4098716 107 2656842 1 119 119 3226663 3551489 4063517 108 2597484 1 120 120 2656842 3226663 3551489 109 2572399 1 121 121 2597484 2656842 3226663 110 2596631 1 122 122 2572399 2597484 2656842 111 3165225 1 123 123 2596631 2572399 2597484 112 3303145 1 124 124 3165225 2596631 2572399 113 3698247 1 125 125 3303145 3165225 2596631 114 3668631 1 126 126 3698247 3303145 3165225 115 4130433 1 127 127 3668631 3698247 3303145 116 4131400 1 128 128 4130433 3668631 3698247 117 3864358 1 129 129 4131400 4130433 3668631 118 3721110 1 130 130 3864358 4131400 4130433 119 2892532 1 131 131 3721110 3864358 4131400 120 2843451 1 132 132 2892532 3721110 3864358 121 2747502 1 133 133 2843451 2892532 3721110 122 2668775 1 134 134 2747502 2843451 2892532 123 3018602 1 135 135 2668775 2747502 2843451 124 3013392 1 136 136 3018602 2668775 2747502 125 3393657 1 137 137 3013392 3018602 2668775 126 3544233 1 138 138 3393657 3013392 3018602 127 4075832 1 139 139 3544233 3393657 3013392 128 4032923 1 140 140 4075832 3544233 3393657 129 3734509 1 141 141 4032923 4075832 3544233 130 3761285 1 142 142 3734509 4032923 4075832 131 2970090 1 143 143 3761285 3734509 4032923 132 2847849 1 144 144 2970090 3761285 3734509 133 2741680 1 145 145 2847849 2970090 3761285 134 2830639 1 146 146 2741680 2847849 2970090 135 3257673 1 147 147 2830639 2741680 2847849 136 3480085 1 148 148 3257673 2830639 2741680 137 3843271 1 149 149 3480085 3257673 2830639 138 3796961 1 150 150 3843271 3480085 3257673 139 4337767 1 151 151 3796961 3843271 3480085 140 4243630 1 152 152 4337767 3796961 3843271 141 3927202 1 153 153 4243630 4337767 3796961 142 3915296 1 154 154 3927202 4243630 4337767 143 3087396 1 155 155 3915296 3927202 4243630 144 2963792 1 156 156 3087396 3915296 3927202 145 2955792 1 157 157 2963792 3087396 3915296 146 2829925 1 158 158 2955792 2963792 3087396 147 3281195 1 159 159 2829925 2955792 2963792 148 3548011 1 160 160 3281195 2829925 2955792 149 4059648 1 161 161 3548011 3281195 2829925 150 3941175 1 162 162 4059648 3548011 3281195 151 4528594 1 163 163 3941175 4059648 3548011 152 4433151 1 164 164 4528594 3941175 4059648 153 4145737 1 165 165 4433151 4528594 3941175 154 4077132 1 166 166 4145737 4433151 4528594 155 3198519 1 167 167 4077132 4145737 4433151 156 3078660 1 168 168 3198519 4077132 4145737 157 3028202 1 169 169 3078660 3198519 4077132 158 2858642 1 170 170 3028202 3078660 3198519 159 3398954 1 171 171 2858642 3028202 3078660 160 3808883 1 172 172 3398954 2858642 3028202 161 4175961 1 173 173 3808883 3398954 2858642 162 4227542 1 174 174 4175961 3808883 3398954 163 4744616 1 175 175 4227542 4175961 3808883 164 4608012 1 176 176 4744616 4227542 4175961 165 4295049 1 177 177 4608012 4744616 4227542 166 4201144 1 178 178 4295049 4608012 4744616 167 3353276 1 179 179 4201144 4295049 4608012 168 3286851 1 180 180 3353276 4201144 4295049 169 3169889 1 181 181 3286851 3353276 4201144 170 3051720 1 182 182 3169889 3286851 3353276 171 3695426 1 183 183 3051720 3169889 3286851 172 3905501 1 184 184 3695426 3051720 3169889 173 4296458 1 185 185 3905501 3695426 3051720 174 4246247 1 186 186 4296458 3905501 3695426 175 4921849 1 187 187 4246247 4296458 3905501 176 4821446 1 188 188 4921849 4246247 4296458 177 4425064 1 189 189 4821446 4921849 4246247 178 4379099 1 190 190 4425064 4821446 4921849 179 3472889 1 191 191 4379099 4425064 4821446 180 3359160 1 192 192 3472889 4379099 4425064 181 3200944 1 193 193 3359160 3472889 4379099 182 3153170 1 194 194 3200944 3359160 3472889 183 3741498 1 195 195 3153170 3200944 3359160 184 3918719 1 196 196 3741498 3153170 3200944 185 4403449 1 197 197 3918719 3741498 3153170 186 4400407 1 198 198 4403449 3918719 3741498 187 4847473 1 199 199 4400407 4403449 3918719 188 4716136 1 200 200 4847473 4400407 4403449 189 4297440 1 201 201 4716136 4847473 4400407 190 4272253 1 202 202 4297440 4716136 4847473 191 3271834 1 203 203 4272253 4297440 4716136 192 3168388 1 204 204 3271834 4272253 4297440 193 2911748 1 205 205 3168388 3271834 4272253 194 2720999 1 206 206 2911748 3168388 3271834 195 3199918 1 207 207 2720999 2911748 3168388 196 3672623 1 208 208 3199918 2720999 2911748 197 3892013 1 209 209 3672623 3199918 2720999 198 3850845 1 210 210 3892013 3672623 3199918 199 4532467 1 211 211 3850845 3892013 3672623 200 4484739 1 212 212 4532467 3850845 3892013 201 4014972 1 213 213 4484739 4532467 3850845 202 3983758 1 214 214 4014972 4484739 4532467 203 3158459 1 215 215 3983758 4014972 4484739 204 3100569 1 216 216 3158459 3983758 4014972 205 2935404 1 217 217 3100569 3158459 3983758 206 2855719 1 218 218 2935404 3100569 3158459 207 3465611 1 219 219 2855719 2935404 3100569 208 3006985 1 220 220 3465611 2855719 2935404 209 4095110 1 221 221 3006985 3465611 2855719 210 4104793 1 222 222 4095110 3006985 3465611 211 4730788 1 223 223 4104793 4095110 3006985 212 4642726 1 224 224 4730788 4104793 4095110 213 4246919 1 225 225 4642726 4730788 4104793 214 4308117 1 226 226 4246919 4642726 4730788 Passagiers-4 M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 M11 1 1867943 1 0 0 0 0 0 0 0 0 0 0 2 1707752 0 1 0 0 0 0 0 0 0 0 0 3 1298756 0 0 1 0 0 0 0 0 0 0 0 4 1281814 0 0 0 1 0 0 0 0 0 0 0 5 1281151 0 0 0 0 1 0 0 0 0 0 0 6 1164976 0 0 0 0 0 1 0 0 0 0 0 7 1454329 0 0 0 0 0 0 1 0 0 0 0 8 1645288 0 0 0 0 0 0 0 1 0 0 0 9 1817743 0 0 0 0 0 0 0 0 1 0 0 10 1895785 0 0 0 0 0 0 0 0 0 1 0 11 2236311 0 0 0 0 0 0 0 0 0 0 1 12 2295951 0 0 0 0 0 0 0 0 0 0 0 13 2087315 1 0 0 0 0 0 0 0 0 0 0 14 1980891 0 1 0 0 0 0 0 0 0 0 0 15 1465446 0 0 1 0 0 0 0 0 0 0 0 16 1445026 0 0 0 1 0 0 0 0 0 0 0 17 1488120 0 0 0 0 1 0 0 0 0 0 0 18 1338333 0 0 0 0 0 1 0 0 0 0 0 19 1715789 0 0 0 0 0 0 1 0 0 0 0 20 1806090 0 0 0 0 0 0 0 1 0 0 0 21 2083316 0 0 0 0 0 0 0 0 1 0 0 22 2092278 0 0 0 0 0 0 0 0 0 1 0 23 2430800 0 0 0 0 0 0 0 0 0 0 1 24 2424894 0 0 0 0 0 0 0 0 0 0 0 25 2299016 1 0 0 0 0 0 0 0 0 0 0 26 2130688 0 1 0 0 0 0 0 0 0 0 0 27 1652221 0 0 1 0 0 0 0 0 0 0 0 28 1608162 0 0 0 1 0 0 0 0 0 0 0 29 1647074 0 0 0 0 1 0 0 0 0 0 0 30 1479691 0 0 0 0 0 1 0 0 0 0 0 31 1884978 0 0 0 0 0 0 1 0 0 0 0 32 2007898 0 0 0 0 0 0 0 1 0 0 0 33 2208954 0 0 0 0 0 0 0 0 1 0 0 34 2217164 0 0 0 0 0 0 0 0 0 1 0 35 2534291 0 0 0 0 0 0 0 0 0 0 1 36 2560312 0 0 0 0 0 0 0 0 0 0 0 37 2429069 1 0 0 0 0 0 0 0 0 0 0 38 2315077 0 1 0 0 0 0 0 0 0 0 0 39 1799608 0 0 1 0 0 0 0 0 0 0 0 40 1772590 0 0 0 1 0 0 0 0 0 0 0 41 1744799 0 0 0 0 1 0 0 0 0 0 0 42 1659093 0 0 0 0 0 1 0 0 0 0 0 43 2099821 0 0 0 0 0 0 1 0 0 0 0 44 2135736 0 0 0 0 0 0 0 1 0 0 0 45 2427894 0 0 0 0 0 0 0 0 1 0 0 46 2468882 0 0 0 0 0 0 0 0 0 1 0 47 2703217 0 0 0 0 0 0 0 0 0 0 1 48 2766841 0 0 0 0 0 0 0 0 0 0 0 49 2655236 1 0 0 0 0 0 0 0 0 0 0 50 2550373 0 1 0 0 0 0 0 0 0 0 0 51 2052097 0 0 1 0 0 0 0 0 0 0 0 52 1998055 0 0 0 1 0 0 0 0 0 0 0 53 1920748 0 0 0 0 1 0 0 0 0 0 0 54 1876694 0 0 0 0 0 1 0 0 0 0 0 55 2380930 0 0 0 0 0 0 1 0 0 0 0 56 2467402 0 0 0 0 0 0 0 1 0 0 0 57 2770771 0 0 0 0 0 0 0 0 1 0 0 58 2781340 0 0 0 0 0 0 0 0 0 1 0 59 3143926 0 0 0 0 0 0 0 0 0 0 1 60 3172235 0 0 0 0 0 0 0 0 0 0 0 61 2952540 1 0 0 0 0 0 0 0 0 0 0 62 2920877 0 1 0 0 0 0 0 0 0 0 0 63 2384552 0 0 1 0 0 0 0 0 0 0 0 64 2248987 0 0 0 1 0 0 0 0 0 0 0 65 2208616 0 0 0 0 1 0 0 0 0 0 0 66 2178756 0 0 0 0 0 1 0 0 0 0 0 67 2632870 0 0 0 0 0 0 1 0 0 0 0 68 2706905 0 0 0 0 0 0 0 1 0 0 0 69 3029745 0 0 0 0 0 0 0 0 1 0 0 70 3015402 0 0 0 0 0 0 0 0 0 1 0 71 3391414 0 0 0 0 0 0 0 0 0 0 1 72 3507805 0 0 0 0 0 0 0 0 0 0 0 73 3177852 1 0 0 0 0 0 0 0 0 0 0 74 3142961 0 1 0 0 0 0 0 0 0 0 0 75 2545815 0 0 1 0 0 0 0 0 0 0 0 76 2414007 0 0 0 1 0 0 0 0 0 0 0 77 2372578 0 0 0 0 1 0 0 0 0 0 0 78 2332664 0 0 0 0 0 1 0 0 0 0 0 79 2825328 0 0 0 0 0 0 1 0 0 0 0 80 2901478 0 0 0 0 0 0 0 1 0 0 0 81 3263955 0 0 0 0 0 0 0 0 1 0 0 82 3226738 0 0 0 0 0 0 0 0 0 1 0 83 3610786 0 0 0 0 0 0 0 0 0 0 1 84 3709274 0 0 0 0 0 0 0 0 0 0 0 85 3467185 1 0 0 0 0 0 0 0 0 0 0 86 3449646 0 1 0 0 0 0 0 0 0 0 0 87 2802951 0 0 1 0 0 0 0 0 0 0 0 88 2462530 0 0 0 1 0 0 0 0 0 0 0 89 2490645 0 0 0 0 1 0 0 0 0 0 0 90 2561520 0 0 0 0 0 1 0 0 0 0 0 91 3067554 0 0 0 0 0 0 1 0 0 0 0 92 3226951 0 0 0 0 0 0 0 1 0 0 0 93 3546493 0 0 0 0 0 0 0 0 1 0 0 94 3492787 0 0 0 0 0 0 0 0 0 1 0 95 3952263 0 0 0 0 0 0 0 0 0 0 1 96 3932072 0 0 0 0 0 0 0 0 0 0 0 97 3720284 1 0 0 0 0 0 0 0 0 0 0 98 3651555 0 1 0 0 0 0 0 0 0 0 0 99 2914972 0 0 1 0 0 0 0 0 0 0 0 100 2713514 0 0 0 1 0 0 0 0 0 0 0 101 2703997 0 0 0 0 1 0 0 0 0 0 0 102 2591373 0 0 0 0 0 1 0 0 0 0 0 103 3163748 0 0 0 0 0 0 1 0 0 0 0 104 3355137 0 0 0 0 0 0 0 1 0 0 0 105 3613702 0 0 0 0 0 0 0 0 1 0 0 106 3686773 0 0 0 0 0 0 0 0 0 1 0 107 4098716 0 0 0 0 0 0 0 0 0 0 1 108 4063517 0 0 0 0 0 0 0 0 0 0 0 109 3551489 1 0 0 0 0 0 0 0 0 0 0 110 3226663 0 1 0 0 0 0 0 0 0 0 0 111 2656842 0 0 1 0 0 0 0 0 0 0 0 112 2597484 0 0 0 1 0 0 0 0 0 0 0 113 2572399 0 0 0 0 1 0 0 0 0 0 0 114 2596631 0 0 0 0 0 1 0 0 0 0 0 115 3165225 0 0 0 0 0 0 1 0 0 0 0 116 3303145 0 0 0 0 0 0 0 1 0 0 0 117 3698247 0 0 0 0 0 0 0 0 1 0 0 118 3668631 0 0 0 0 0 0 0 0 0 1 0 119 4130433 0 0 0 0 0 0 0 0 0 0 1 120 4131400 0 0 0 0 0 0 0 0 0 0 0 121 3864358 1 0 0 0 0 0 0 0 0 0 0 122 3721110 0 1 0 0 0 0 0 0 0 0 0 123 2892532 0 0 1 0 0 0 0 0 0 0 0 124 2843451 0 0 0 1 0 0 0 0 0 0 0 125 2747502 0 0 0 0 1 0 0 0 0 0 0 126 2668775 0 0 0 0 0 1 0 0 0 0 0 127 3018602 0 0 0 0 0 0 1 0 0 0 0 128 3013392 0 0 0 0 0 0 0 1 0 0 0 129 3393657 0 0 0 0 0 0 0 0 1 0 0 130 3544233 0 0 0 0 0 0 0 0 0 1 0 131 4075832 0 0 0 0 0 0 0 0 0 0 1 132 4032923 0 0 0 0 0 0 0 0 0 0 0 133 3734509 1 0 0 0 0 0 0 0 0 0 0 134 3761285 0 1 0 0 0 0 0 0 0 0 0 135 2970090 0 0 1 0 0 0 0 0 0 0 0 136 2847849 0 0 0 1 0 0 0 0 0 0 0 137 2741680 0 0 0 0 1 0 0 0 0 0 0 138 2830639 0 0 0 0 0 1 0 0 0 0 0 139 3257673 0 0 0 0 0 0 1 0 0 0 0 140 3480085 0 0 0 0 0 0 0 1 0 0 0 141 3843271 0 0 0 0 0 0 0 0 1 0 0 142 3796961 0 0 0 0 0 0 0 0 0 1 0 143 4337767 0 0 0 0 0 0 0 0 0 0 1 144 4243630 0 0 0 0 0 0 0 0 0 0 0 145 3927202 1 0 0 0 0 0 0 0 0 0 0 146 3915296 0 1 0 0 0 0 0 0 0 0 0 147 3087396 0 0 1 0 0 0 0 0 0 0 0 148 2963792 0 0 0 1 0 0 0 0 0 0 0 149 2955792 0 0 0 0 1 0 0 0 0 0 0 150 2829925 0 0 0 0 0 1 0 0 0 0 0 151 3281195 0 0 0 0 0 0 1 0 0 0 0 152 3548011 0 0 0 0 0 0 0 1 0 0 0 153 4059648 0 0 0 0 0 0 0 0 1 0 0 154 3941175 0 0 0 0 0 0 0 0 0 1 0 155 4528594 0 0 0 0 0 0 0 0 0 0 1 156 4433151 0 0 0 0 0 0 0 0 0 0 0 157 4145737 1 0 0 0 0 0 0 0 0 0 0 158 4077132 0 1 0 0 0 0 0 0 0 0 0 159 3198519 0 0 1 0 0 0 0 0 0 0 0 160 3078660 0 0 0 1 0 0 0 0 0 0 0 161 3028202 0 0 0 0 1 0 0 0 0 0 0 162 2858642 0 0 0 0 0 1 0 0 0 0 0 163 3398954 0 0 0 0 0 0 1 0 0 0 0 164 3808883 0 0 0 0 0 0 0 1 0 0 0 165 4175961 0 0 0 0 0 0 0 0 1 0 0 166 4227542 0 0 0 0 0 0 0 0 0 1 0 167 4744616 0 0 0 0 0 0 0 0 0 0 1 168 4608012 0 0 0 0 0 0 0 0 0 0 0 169 4295049 1 0 0 0 0 0 0 0 0 0 0 170 4201144 0 1 0 0 0 0 0 0 0 0 0 171 3353276 0 0 1 0 0 0 0 0 0 0 0 172 3286851 0 0 0 1 0 0 0 0 0 0 0 173 3169889 0 0 0 0 1 0 0 0 0 0 0 174 3051720 0 0 0 0 0 1 0 0 0 0 0 175 3695426 0 0 0 0 0 0 1 0 0 0 0 176 3905501 0 0 0 0 0 0 0 1 0 0 0 177 4296458 0 0 0 0 0 0 0 0 1 0 0 178 4246247 0 0 0 0 0 0 0 0 0 1 0 179 4921849 0 0 0 0 0 0 0 0 0 0 1 180 4821446 0 0 0 0 0 0 0 0 0 0 0 181 4425064 1 0 0 0 0 0 0 0 0 0 0 182 4379099 0 1 0 0 0 0 0 0 0 0 0 183 3472889 0 0 1 0 0 0 0 0 0 0 0 184 3359160 0 0 0 1 0 0 0 0 0 0 0 185 3200944 0 0 0 0 1 0 0 0 0 0 0 186 3153170 0 0 0 0 0 1 0 0 0 0 0 187 3741498 0 0 0 0 0 0 1 0 0 0 0 188 3918719 0 0 0 0 0 0 0 1 0 0 0 189 4403449 0 0 0 0 0 0 0 0 1 0 0 190 4400407 0 0 0 0 0 0 0 0 0 1 0 191 4847473 0 0 0 0 0 0 0 0 0 0 1 192 4716136 0 0 0 0 0 0 0 0 0 0 0 193 4297440 1 0 0 0 0 0 0 0 0 0 0 194 4272253 0 1 0 0 0 0 0 0 0 0 0 195 3271834 0 0 1 0 0 0 0 0 0 0 0 196 3168388 0 0 0 1 0 0 0 0 0 0 0 197 2911748 0 0 0 0 1 0 0 0 0 0 0 198 2720999 0 0 0 0 0 1 0 0 0 0 0 199 3199918 0 0 0 0 0 0 1 0 0 0 0 200 3672623 0 0 0 0 0 0 0 1 0 0 0 201 3892013 0 0 0 0 0 0 0 0 1 0 0 202 3850845 0 0 0 0 0 0 0 0 0 1 0 203 4532467 0 0 0 0 0 0 0 0 0 0 1 204 4484739 0 0 0 0 0 0 0 0 0 0 0 205 4014972 1 0 0 0 0 0 0 0 0 0 0 206 3983758 0 1 0 0 0 0 0 0 0 0 0 207 3158459 0 0 1 0 0 0 0 0 0 0 0 208 3100569 0 0 0 1 0 0 0 0 0 0 0 209 2935404 0 0 0 0 1 0 0 0 0 0 0 210 2855719 0 0 0 0 0 1 0 0 0 0 0 211 3465611 0 0 0 0 0 0 1 0 0 0 0 212 3006985 0 0 0 0 0 0 0 1 0 0 0 213 4095110 0 0 0 0 0 0 0 0 1 0 0 214 4104793 0 0 0 0 0 0 0 0 0 1 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) `9/11` t `T_9/11` `Passagiers-1` 1.940e+05 3.801e+05 5.405e+03 -3.722e+03 6.587e-01 `Passagiers-2` `Passagiers-3` `Passagiers-4` M1 M2 2.647e-01 -2.761e-03 -2.097e-01 1.343e+05 1.003e+05 M3 M4 M5 M6 M7 5.012e+05 3.151e+05 4.537e+05 1.558e+05 6.148e+05 M8 M9 M10 M11 3.144e+05 -7.045e+03 1.172e+05 -3.740e+05 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -632906 -50188 -2490 62624 420778 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.940e+05 9.745e+04 1.991 0.047893 * `9/11` 3.801e+05 9.458e+04 4.019 8.36e-05 *** t 5.405e+03 1.010e+03 5.352 2.43e-07 *** `T_9/11` -3.722e+03 8.183e+02 -4.548 9.50e-06 *** `Passagiers-1` 6.587e-01 6.991e-02 9.422 < 2e-16 *** `Passagiers-2` 2.647e-01 8.454e-02 3.131 0.002013 ** `Passagiers-3` -2.761e-03 8.369e-02 -0.033 0.973716 `Passagiers-4` -2.097e-01 6.920e-02 -3.031 0.002768 ** M1 1.343e+05 6.410e+04 2.095 0.037451 * M2 1.003e+05 6.629e+04 1.513 0.131935 M3 5.012e+05 6.780e+04 7.392 4.15e-12 *** M4 3.151e+05 9.234e+04 3.413 0.000782 *** M5 4.537e+05 8.790e+04 5.162 6.00e-07 *** M6 1.558e+05 9.935e+04 1.569 0.118379 M7 6.148e+05 7.784e+04 7.898 2.02e-13 *** M8 3.144e+05 9.604e+04 3.274 0.001257 ** M9 -7.045e+03 7.884e+04 -0.089 0.928892 M10 1.172e+05 6.842e+04 1.713 0.088297 . M11 -3.740e+05 6.017e+04 -6.215 3.04e-09 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 111900 on 195 degrees of freedom Multiple R-squared: 0.9848, Adjusted R-squared: 0.9833 F-statistic: 699.7 on 18 and 195 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 5.503968e-02 1.100794e-01 0.9449603 [2,] 2.268071e-02 4.536142e-02 0.9773193 [3,] 6.910215e-03 1.382043e-02 0.9930898 [4,] 2.244675e-03 4.489349e-03 0.9977553 [5,] 8.240686e-04 1.648137e-03 0.9991759 [6,] 2.592132e-04 5.184264e-04 0.9997408 [7,] 8.320187e-05 1.664037e-04 0.9999168 [8,] 3.276403e-05 6.552806e-05 0.9999672 [9,] 2.100580e-05 4.201159e-05 0.9999790 [10,] 1.997461e-05 3.994922e-05 0.9999800 [11,] 7.475639e-06 1.495128e-05 0.9999925 [12,] 2.297465e-06 4.594930e-06 0.9999977 [13,] 7.538132e-07 1.507626e-06 0.9999992 [14,] 2.227786e-07 4.455572e-07 0.9999998 [15,] 5.942688e-08 1.188538e-07 0.9999999 [16,] 4.056417e-08 8.112834e-08 1.0000000 [17,] 1.061834e-08 2.123668e-08 1.0000000 [18,] 7.000681e-08 1.400136e-07 0.9999999 [19,] 3.927747e-08 7.855494e-08 1.0000000 [20,] 1.721160e-08 3.442320e-08 1.0000000 [21,] 1.618468e-08 3.236937e-08 1.0000000 [22,] 1.635081e-08 3.270162e-08 1.0000000 [23,] 5.746189e-09 1.149238e-08 1.0000000 [24,] 3.444419e-09 6.888838e-09 1.0000000 [25,] 2.783449e-09 5.566898e-09 1.0000000 [26,] 1.854487e-09 3.708973e-09 1.0000000 [27,] 6.421296e-10 1.284259e-09 1.0000000 [28,] 7.131577e-10 1.426315e-09 1.0000000 [29,] 3.194761e-10 6.389521e-10 1.0000000 [30,] 9.305914e-09 1.861183e-08 1.0000000 [31,] 7.979599e-09 1.595920e-08 1.0000000 [32,] 3.219159e-08 6.438319e-08 1.0000000 [33,] 1.535458e-08 3.070916e-08 1.0000000 [34,] 1.264126e-08 2.528253e-08 1.0000000 [35,] 5.135807e-09 1.027161e-08 1.0000000 [36,] 4.393649e-09 8.787298e-09 1.0000000 [37,] 3.208311e-09 6.416621e-09 1.0000000 [38,] 1.809115e-09 3.618231e-09 1.0000000 [39,] 2.702841e-09 5.405682e-09 1.0000000 [40,] 2.452507e-09 4.905014e-09 1.0000000 [41,] 1.528883e-09 3.057766e-09 1.0000000 [42,] 1.089229e-09 2.178459e-09 1.0000000 [43,] 4.700153e-10 9.400306e-10 1.0000000 [44,] 5.599262e-10 1.119852e-09 1.0000000 [45,] 2.424227e-10 4.848454e-10 1.0000000 [46,] 2.206127e-10 4.412254e-10 1.0000000 [47,] 3.785197e-10 7.570395e-10 1.0000000 [48,] 1.062501e-09 2.125001e-09 1.0000000 [49,] 4.803793e-10 9.607586e-10 1.0000000 [50,] 5.200693e-10 1.040139e-09 1.0000000 [51,] 6.021107e-10 1.204221e-09 1.0000000 [52,] 5.524302e-10 1.104860e-09 1.0000000 [53,] 2.812170e-10 5.624341e-10 1.0000000 [54,] 2.667208e-10 5.334417e-10 1.0000000 [55,] 1.264864e-10 2.529728e-10 1.0000000 [56,] 1.606567e-10 3.213133e-10 1.0000000 [57,] 8.500874e-11 1.700175e-10 1.0000000 [58,] 1.017013e-10 2.034027e-10 1.0000000 [59,] 8.019629e-11 1.603926e-10 1.0000000 [60,] 3.762807e-11 7.525614e-11 1.0000000 [61,] 2.849166e-11 5.698331e-11 1.0000000 [62,] 3.560369e-11 7.120738e-11 1.0000000 [63,] 2.926880e-08 5.853760e-08 1.0000000 [64,] 1.715960e-08 3.431920e-08 1.0000000 [65,] 3.547913e-08 7.095826e-08 1.0000000 [66,] 2.926445e-08 5.852890e-08 1.0000000 [67,] 1.972926e-08 3.945853e-08 1.0000000 [68,] 1.480418e-08 2.960836e-08 1.0000000 [69,] 9.973404e-09 1.994681e-08 1.0000000 [70,] 3.350715e-08 6.701430e-08 1.0000000 [71,] 2.081733e-08 4.163467e-08 1.0000000 [72,] 1.154438e-08 2.308876e-08 1.0000000 [73,] 6.461862e-09 1.292372e-08 1.0000000 [74,] 1.644439e-08 3.288878e-08 1.0000000 [75,] 2.264896e-08 4.529793e-08 1.0000000 [76,] 1.381039e-08 2.762079e-08 1.0000000 [77,] 1.007923e-08 2.015847e-08 1.0000000 [78,] 7.304676e-09 1.460935e-08 1.0000000 [79,] 5.413549e-09 1.082710e-08 1.0000000 [80,] 3.725478e-09 7.450956e-09 1.0000000 [81,] 2.432209e-09 4.864419e-09 1.0000000 [82,] 1.446065e-09 2.892131e-09 1.0000000 [83,] 9.020298e-10 1.804060e-09 1.0000000 [84,] 7.721606e-10 1.544321e-09 1.0000000 [85,] 4.881690e-09 9.763379e-09 1.0000000 [86,] 4.504987e-09 9.009974e-09 1.0000000 [87,] 2.498152e-09 4.996305e-09 1.0000000 [88,] 3.898742e-09 7.797485e-09 1.0000000 [89,] 2.861289e-09 5.722578e-09 1.0000000 [90,] 1.872157e-09 3.744314e-09 1.0000000 [91,] 1.292028e-09 2.584057e-09 1.0000000 [92,] 7.587223e-10 1.517445e-09 1.0000000 [93,] 8.640420e-10 1.728084e-09 1.0000000 [94,] 8.080094e-10 1.616019e-09 1.0000000 [95,] 6.342273e-10 1.268455e-09 1.0000000 [96,] 5.907779e-10 1.181556e-09 1.0000000 [97,] 5.575577e-10 1.115115e-09 1.0000000 [98,] 1.218475e-08 2.436950e-08 1.0000000 [99,] 6.830575e-09 1.366115e-08 1.0000000 [100,] 8.468299e-09 1.693660e-08 1.0000000 [101,] 5.355428e-09 1.071086e-08 1.0000000 [102,] 3.806004e-08 7.612008e-08 1.0000000 [103,] 4.703977e-07 9.407955e-07 0.9999995 [104,] 1.688519e-06 3.377038e-06 0.9999983 [105,] 1.214998e-06 2.429995e-06 0.9999988 [106,] 1.073698e-06 2.147396e-06 0.9999989 [107,] 1.177217e-06 2.354433e-06 0.9999988 [108,] 9.138129e-07 1.827626e-06 0.9999991 [109,] 9.617738e-07 1.923548e-06 0.9999990 [110,] 1.452408e-06 2.904815e-06 0.9999985 [111,] 8.506275e-07 1.701255e-06 0.9999991 [112,] 6.795914e-07 1.359183e-06 0.9999993 [113,] 1.806171e-06 3.612342e-06 0.9999982 [114,] 1.115100e-06 2.230199e-06 0.9999989 [115,] 8.652097e-07 1.730419e-06 0.9999991 [116,] 1.031221e-06 2.062443e-06 0.9999990 [117,] 9.103058e-07 1.820612e-06 0.9999991 [118,] 2.387738e-06 4.775477e-06 0.9999976 [119,] 1.895158e-06 3.790316e-06 0.9999981 [120,] 1.147474e-06 2.294948e-06 0.9999989 [121,] 9.688313e-07 1.937663e-06 0.9999990 [122,] 1.002739e-06 2.005477e-06 0.9999990 [123,] 5.902345e-07 1.180469e-06 0.9999994 [124,] 4.166478e-07 8.332955e-07 0.9999996 [125,] 2.653059e-07 5.306119e-07 0.9999997 [126,] 2.167876e-07 4.335751e-07 0.9999998 [127,] 1.625023e-07 3.250047e-07 0.9999998 [128,] 2.830914e-07 5.661829e-07 0.9999997 [129,] 3.353270e-07 6.706540e-07 0.9999997 [130,] 5.903468e-07 1.180694e-06 0.9999994 [131,] 4.385780e-07 8.771560e-07 0.9999996 [132,] 2.534046e-07 5.068092e-07 0.9999997 [133,] 1.602647e-07 3.205293e-07 0.9999998 [134,] 1.576951e-07 3.153902e-07 0.9999998 [135,] 1.028323e-07 2.056646e-07 0.9999999 [136,] 5.671141e-08 1.134228e-07 0.9999999 [137,] 4.142025e-08 8.284050e-08 1.0000000 [138,] 2.926407e-08 5.852813e-08 1.0000000 [139,] 1.828126e-07 3.656252e-07 0.9999998 [140,] 1.096685e-07 2.193370e-07 0.9999999 [141,] 5.846766e-08 1.169353e-07 0.9999999 [142,] 5.294914e-08 1.058983e-07 0.9999999 [143,] 4.315926e-08 8.631852e-08 1.0000000 [144,] 2.153324e-08 4.306647e-08 1.0000000 [145,] 1.673689e-08 3.347378e-08 1.0000000 [146,] 9.525291e-09 1.905058e-08 1.0000000 [147,] 4.921374e-09 9.842748e-09 1.0000000 [148,] 2.928608e-09 5.857215e-09 1.0000000 [149,] 1.433057e-09 2.866113e-09 1.0000000 [150,] 1.109205e-09 2.218410e-09 1.0000000 [151,] 2.290698e-09 4.581396e-09 1.0000000 [152,] 1.124213e-09 2.248427e-09 1.0000000 [153,] 6.853806e-10 1.370761e-09 1.0000000 [154,] 1.809688e-09 3.619375e-09 1.0000000 [155,] 7.887702e-10 1.577540e-09 1.0000000 [156,] 4.459495e-10 8.918991e-10 1.0000000 [157,] 1.927033e-10 3.854065e-10 1.0000000 [158,] 1.137755e-10 2.275510e-10 1.0000000 [159,] 4.399068e-11 8.798135e-11 1.0000000 [160,] 3.241201e-11 6.482401e-11 1.0000000 [161,] 3.058442e-11 6.116883e-11 1.0000000 [162,] 4.442718e-11 8.885435e-11 1.0000000 [163,] 2.284143e-09 4.568286e-09 1.0000000 [164,] 4.847850e-09 9.695700e-09 1.0000000 [165,] 2.040623e-08 4.081245e-08 1.0000000 [166,] 8.609384e-09 1.721877e-08 1.0000000 [167,] 3.584265e-09 7.168530e-09 1.0000000 [168,] 2.127695e-09 4.255390e-09 1.0000000 [169,] 1.149321e-09 2.298643e-09 1.0000000 [170,] 2.347826e-09 4.695651e-09 1.0000000 [171,] 1.289481e-09 2.578961e-09 1.0000000 > postscript(file="/var/www/rcomp/tmp/191aa1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/2karv1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/3karv1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/4karv1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5vj8x1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') > qqline(mysum$resid) > grid() > dev.off() null device 1 > (myerror <- as.ts(mysum$resid)) Time Series: Start = 1 End = 214 Frequency = 1 1 2 3 4 5 6 90945.6890 -26431.6816 -152473.9555 55710.4893 -118660.0363 64123.7765 7 8 9 10 11 12 -95510.1935 54659.5163 69004.7594 -28123.8249 139124.4708 118966.5259 13 14 15 16 17 18 128188.4557 -39721.7390 -89436.6614 -31658.9483 -49185.9588 15337.1862 19 20 21 22 23 24 -110325.6568 -26920.0787 135674.2027 -75019.3021 147479.8459 82190.0877 25 26 27 28 29 30 110200.8293 -79175.1075 -80670.8124 -8919.5067 -132395.6533 -30695.8873 31 32 33 34 35 36 -151170.7823 -14914.5713 110968.5848 -50521.9376 96203.1530 64628.7423 37 38 39 40 41 42 23529.3743 -33446.4849 -43363.3649 -100191.3284 -98397.1124 16325.1713 43 44 45 46 47 48 -225453.3180 -23728.1457 138133.4213 -30470.8909 104975.4532 40576.0430 49 50 51 52 53 54 -32644.3365 -6235.5921 36555.0189 -28379.6699 -75757.2988 -3309.1513 55 56 57 58 59 60 -86306.5507 14306.0401 59663.0798 38782.5595 143370.5497 -4567.2274 61 62 63 64 65 66 453.0063 53549.4485 18886.6483 -46220.9533 -44882.8762 -4042.5790 67 68 69 70 71 72 -72903.7388 110985.6216 -11476.5180 8544.5494 86679.7326 1587.4990 73 74 75 76 77 78 -3960.9673 37939.0716 35983.3656 -48957.3089 -19828.3947 -30534.3341 79 80 81 82 83 84 -78691.5335 88608.4850 71926.7697 51384.1134 46920.3948 -222250.1779 85 86 87 88 89 90 10756.1608 176351.7003 85410.7913 2025.4305 -55265.3548 -52937.9954 91 92 93 94 95 96 -407.4535 20232.0154 83014.6718 19483.2485 -33675.7533 -115929.6101 97 98 99 100 101 102 17910.1839 -22968.5277 64799.9742 47314.7810 -117985.4785 4516.1914 103 104 105 106 107 108 -43857.0717 -33927.4872 -206026.0388 -293724.3211 61755.6975 79277.8501 109 110 111 112 113 114 -150.4262 18941.2861 55987.1307 -15217.2609 -6923.6792 -30484.1118 115 116 117 118 119 120 5287.7493 38602.6563 51213.7185 -47258.6252 -124432.2239 34021.0692 121 122 123 124 125 126 -2670.9085 -5209.2059 -154593.4596 -195608.0717 -65111.2965 116979.0506 127 128 129 130 131 132 61498.3824 -72793.0986 -83739.3648 58076.0874 -70893.9582 -64516.3751 133 134 135 136 137 138 -79240.1214 147766.5629 -24527.6230 51466.0230 -7173.8251 -35582.5990 139 140 141 142 143 144 69180.9646 -22599.5075 -24376.7115 62913.8086 -70722.6475 -42085.0054 145 146 147 148 149 150 48088.1086 -12251.9365 -52482.0974 108795.0843 182937.3208 -72166.1508 151 152 153 154 155 156 92660.8092 -2308.4495 44387.2132 41205.5935 -103700.9420 -23098.8464 157 158 159 160 161 162 41500.5171 -47586.9472 30618.6918 288586.6540 91303.7015 54683.3950 163 164 165 166 167 168 94476.7178 -10707.7616 26322.3494 61036.8960 -44563.8256 67212.4307 169 170 171 172 173 174 16543.1990 3291.0507 175235.5645 162661.7182 79736.1003 -10445.5213 175 176 177 178 179 180 269742.7892 81423.3615 -26043.8208 81068.5625 -59026.1680 38551.5690 181 182 183 184 185 186 -24132.5768 82599.4785 151346.5140 113730.9565 152419.6456 70942.5853 187 188 189 190 191 192 54998.2251 -32827.0687 -61945.2561 98093.5329 -192008.4879 -34138.0828 193 194 195 196 197 198 -181707.6964 -151739.3894 -91898.8012 277768.0944 -135608.0500 -188919.2089 199 200 201 202 203 204 102905.7018 15510.8803 -237591.8634 -79415.9118 -127485.2911 -20426.4917 205 206 207 208 209 210 -163608.4908 -95671.9868 34623.0762 -632906.1831 420778.2462 116210.1826 211 212 213 214 113874.9593 -183602.4075 -139109.1972 83945.8618 > postscript(file="/var/www/rcomp/tmp/6vj8x1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > dum <- cbind(lag(myerror,k=1),myerror) > dum Time Series: Start = 0 End = 214 Frequency = 1 lag(myerror, k = 1) myerror 0 90945.6890 NA 1 -26431.6816 90945.6890 2 -152473.9555 -26431.6816 3 55710.4893 -152473.9555 4 -118660.0363 55710.4893 5 64123.7765 -118660.0363 6 -95510.1935 64123.7765 7 54659.5163 -95510.1935 8 69004.7594 54659.5163 9 -28123.8249 69004.7594 10 139124.4708 -28123.8249 11 118966.5259 139124.4708 12 128188.4557 118966.5259 13 -39721.7390 128188.4557 14 -89436.6614 -39721.7390 15 -31658.9483 -89436.6614 16 -49185.9588 -31658.9483 17 15337.1862 -49185.9588 18 -110325.6568 15337.1862 19 -26920.0787 -110325.6568 20 135674.2027 -26920.0787 21 -75019.3021 135674.2027 22 147479.8459 -75019.3021 23 82190.0877 147479.8459 24 110200.8293 82190.0877 25 -79175.1075 110200.8293 26 -80670.8124 -79175.1075 27 -8919.5067 -80670.8124 28 -132395.6533 -8919.5067 29 -30695.8873 -132395.6533 30 -151170.7823 -30695.8873 31 -14914.5713 -151170.7823 32 110968.5848 -14914.5713 33 -50521.9376 110968.5848 34 96203.1530 -50521.9376 35 64628.7423 96203.1530 36 23529.3743 64628.7423 37 -33446.4849 23529.3743 38 -43363.3649 -33446.4849 39 -100191.3284 -43363.3649 40 -98397.1124 -100191.3284 41 16325.1713 -98397.1124 42 -225453.3180 16325.1713 43 -23728.1457 -225453.3180 44 138133.4213 -23728.1457 45 -30470.8909 138133.4213 46 104975.4532 -30470.8909 47 40576.0430 104975.4532 48 -32644.3365 40576.0430 49 -6235.5921 -32644.3365 50 36555.0189 -6235.5921 51 -28379.6699 36555.0189 52 -75757.2988 -28379.6699 53 -3309.1513 -75757.2988 54 -86306.5507 -3309.1513 55 14306.0401 -86306.5507 56 59663.0798 14306.0401 57 38782.5595 59663.0798 58 143370.5497 38782.5595 59 -4567.2274 143370.5497 60 453.0063 -4567.2274 61 53549.4485 453.0063 62 18886.6483 53549.4485 63 -46220.9533 18886.6483 64 -44882.8762 -46220.9533 65 -4042.5790 -44882.8762 66 -72903.7388 -4042.5790 67 110985.6216 -72903.7388 68 -11476.5180 110985.6216 69 8544.5494 -11476.5180 70 86679.7326 8544.5494 71 1587.4990 86679.7326 72 -3960.9673 1587.4990 73 37939.0716 -3960.9673 74 35983.3656 37939.0716 75 -48957.3089 35983.3656 76 -19828.3947 -48957.3089 77 -30534.3341 -19828.3947 78 -78691.5335 -30534.3341 79 88608.4850 -78691.5335 80 71926.7697 88608.4850 81 51384.1134 71926.7697 82 46920.3948 51384.1134 83 -222250.1779 46920.3948 84 10756.1608 -222250.1779 85 176351.7003 10756.1608 86 85410.7913 176351.7003 87 2025.4305 85410.7913 88 -55265.3548 2025.4305 89 -52937.9954 -55265.3548 90 -407.4535 -52937.9954 91 20232.0154 -407.4535 92 83014.6718 20232.0154 93 19483.2485 83014.6718 94 -33675.7533 19483.2485 95 -115929.6101 -33675.7533 96 17910.1839 -115929.6101 97 -22968.5277 17910.1839 98 64799.9742 -22968.5277 99 47314.7810 64799.9742 100 -117985.4785 47314.7810 101 4516.1914 -117985.4785 102 -43857.0717 4516.1914 103 -33927.4872 -43857.0717 104 -206026.0388 -33927.4872 105 -293724.3211 -206026.0388 106 61755.6975 -293724.3211 107 79277.8501 61755.6975 108 -150.4262 79277.8501 109 18941.2861 -150.4262 110 55987.1307 18941.2861 111 -15217.2609 55987.1307 112 -6923.6792 -15217.2609 113 -30484.1118 -6923.6792 114 5287.7493 -30484.1118 115 38602.6563 5287.7493 116 51213.7185 38602.6563 117 -47258.6252 51213.7185 118 -124432.2239 -47258.6252 119 34021.0692 -124432.2239 120 -2670.9085 34021.0692 121 -5209.2059 -2670.9085 122 -154593.4596 -5209.2059 123 -195608.0717 -154593.4596 124 -65111.2965 -195608.0717 125 116979.0506 -65111.2965 126 61498.3824 116979.0506 127 -72793.0986 61498.3824 128 -83739.3648 -72793.0986 129 58076.0874 -83739.3648 130 -70893.9582 58076.0874 131 -64516.3751 -70893.9582 132 -79240.1214 -64516.3751 133 147766.5629 -79240.1214 134 -24527.6230 147766.5629 135 51466.0230 -24527.6230 136 -7173.8251 51466.0230 137 -35582.5990 -7173.8251 138 69180.9646 -35582.5990 139 -22599.5075 69180.9646 140 -24376.7115 -22599.5075 141 62913.8086 -24376.7115 142 -70722.6475 62913.8086 143 -42085.0054 -70722.6475 144 48088.1086 -42085.0054 145 -12251.9365 48088.1086 146 -52482.0974 -12251.9365 147 108795.0843 -52482.0974 148 182937.3208 108795.0843 149 -72166.1508 182937.3208 150 92660.8092 -72166.1508 151 -2308.4495 92660.8092 152 44387.2132 -2308.4495 153 41205.5935 44387.2132 154 -103700.9420 41205.5935 155 -23098.8464 -103700.9420 156 41500.5171 -23098.8464 157 -47586.9472 41500.5171 158 30618.6918 -47586.9472 159 288586.6540 30618.6918 160 91303.7015 288586.6540 161 54683.3950 91303.7015 162 94476.7178 54683.3950 163 -10707.7616 94476.7178 164 26322.3494 -10707.7616 165 61036.8960 26322.3494 166 -44563.8256 61036.8960 167 67212.4307 -44563.8256 168 16543.1990 67212.4307 169 3291.0507 16543.1990 170 175235.5645 3291.0507 171 162661.7182 175235.5645 172 79736.1003 162661.7182 173 -10445.5213 79736.1003 174 269742.7892 -10445.5213 175 81423.3615 269742.7892 176 -26043.8208 81423.3615 177 81068.5625 -26043.8208 178 -59026.1680 81068.5625 179 38551.5690 -59026.1680 180 -24132.5768 38551.5690 181 82599.4785 -24132.5768 182 151346.5140 82599.4785 183 113730.9565 151346.5140 184 152419.6456 113730.9565 185 70942.5853 152419.6456 186 54998.2251 70942.5853 187 -32827.0687 54998.2251 188 -61945.2561 -32827.0687 189 98093.5329 -61945.2561 190 -192008.4879 98093.5329 191 -34138.0828 -192008.4879 192 -181707.6964 -34138.0828 193 -151739.3894 -181707.6964 194 -91898.8012 -151739.3894 195 277768.0944 -91898.8012 196 -135608.0500 277768.0944 197 -188919.2089 -135608.0500 198 102905.7018 -188919.2089 199 15510.8803 102905.7018 200 -237591.8634 15510.8803 201 -79415.9118 -237591.8634 202 -127485.2911 -79415.9118 203 -20426.4917 -127485.2911 204 -163608.4908 -20426.4917 205 -95671.9868 -163608.4908 206 34623.0762 -95671.9868 207 -632906.1831 34623.0762 208 420778.2462 -632906.1831 209 116210.1826 420778.2462 210 113874.9593 116210.1826 211 -183602.4075 113874.9593 212 -139109.1972 -183602.4075 213 83945.8618 -139109.1972 214 NA 83945.8618 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -26431.6816 90945.6890 [2,] -152473.9555 -26431.6816 [3,] 55710.4893 -152473.9555 [4,] -118660.0363 55710.4893 [5,] 64123.7765 -118660.0363 [6,] -95510.1935 64123.7765 [7,] 54659.5163 -95510.1935 [8,] 69004.7594 54659.5163 [9,] -28123.8249 69004.7594 [10,] 139124.4708 -28123.8249 [11,] 118966.5259 139124.4708 [12,] 128188.4557 118966.5259 [13,] -39721.7390 128188.4557 [14,] -89436.6614 -39721.7390 [15,] -31658.9483 -89436.6614 [16,] -49185.9588 -31658.9483 [17,] 15337.1862 -49185.9588 [18,] -110325.6568 15337.1862 [19,] -26920.0787 -110325.6568 [20,] 135674.2027 -26920.0787 [21,] -75019.3021 135674.2027 [22,] 147479.8459 -75019.3021 [23,] 82190.0877 147479.8459 [24,] 110200.8293 82190.0877 [25,] -79175.1075 110200.8293 [26,] -80670.8124 -79175.1075 [27,] -8919.5067 -80670.8124 [28,] -132395.6533 -8919.5067 [29,] -30695.8873 -132395.6533 [30,] -151170.7823 -30695.8873 [31,] -14914.5713 -151170.7823 [32,] 110968.5848 -14914.5713 [33,] -50521.9376 110968.5848 [34,] 96203.1530 -50521.9376 [35,] 64628.7423 96203.1530 [36,] 23529.3743 64628.7423 [37,] -33446.4849 23529.3743 [38,] -43363.3649 -33446.4849 [39,] -100191.3284 -43363.3649 [40,] -98397.1124 -100191.3284 [41,] 16325.1713 -98397.1124 [42,] -225453.3180 16325.1713 [43,] -23728.1457 -225453.3180 [44,] 138133.4213 -23728.1457 [45,] -30470.8909 138133.4213 [46,] 104975.4532 -30470.8909 [47,] 40576.0430 104975.4532 [48,] -32644.3365 40576.0430 [49,] -6235.5921 -32644.3365 [50,] 36555.0189 -6235.5921 [51,] -28379.6699 36555.0189 [52,] -75757.2988 -28379.6699 [53,] -3309.1513 -75757.2988 [54,] -86306.5507 -3309.1513 [55,] 14306.0401 -86306.5507 [56,] 59663.0798 14306.0401 [57,] 38782.5595 59663.0798 [58,] 143370.5497 38782.5595 [59,] -4567.2274 143370.5497 [60,] 453.0063 -4567.2274 [61,] 53549.4485 453.0063 [62,] 18886.6483 53549.4485 [63,] -46220.9533 18886.6483 [64,] -44882.8762 -46220.9533 [65,] -4042.5790 -44882.8762 [66,] -72903.7388 -4042.5790 [67,] 110985.6216 -72903.7388 [68,] -11476.5180 110985.6216 [69,] 8544.5494 -11476.5180 [70,] 86679.7326 8544.5494 [71,] 1587.4990 86679.7326 [72,] -3960.9673 1587.4990 [73,] 37939.0716 -3960.9673 [74,] 35983.3656 37939.0716 [75,] -48957.3089 35983.3656 [76,] -19828.3947 -48957.3089 [77,] -30534.3341 -19828.3947 [78,] -78691.5335 -30534.3341 [79,] 88608.4850 -78691.5335 [80,] 71926.7697 88608.4850 [81,] 51384.1134 71926.7697 [82,] 46920.3948 51384.1134 [83,] -222250.1779 46920.3948 [84,] 10756.1608 -222250.1779 [85,] 176351.7003 10756.1608 [86,] 85410.7913 176351.7003 [87,] 2025.4305 85410.7913 [88,] -55265.3548 2025.4305 [89,] -52937.9954 -55265.3548 [90,] -407.4535 -52937.9954 [91,] 20232.0154 -407.4535 [92,] 83014.6718 20232.0154 [93,] 19483.2485 83014.6718 [94,] -33675.7533 19483.2485 [95,] -115929.6101 -33675.7533 [96,] 17910.1839 -115929.6101 [97,] -22968.5277 17910.1839 [98,] 64799.9742 -22968.5277 [99,] 47314.7810 64799.9742 [100,] -117985.4785 47314.7810 [101,] 4516.1914 -117985.4785 [102,] -43857.0717 4516.1914 [103,] -33927.4872 -43857.0717 [104,] -206026.0388 -33927.4872 [105,] -293724.3211 -206026.0388 [106,] 61755.6975 -293724.3211 [107,] 79277.8501 61755.6975 [108,] -150.4262 79277.8501 [109,] 18941.2861 -150.4262 [110,] 55987.1307 18941.2861 [111,] -15217.2609 55987.1307 [112,] -6923.6792 -15217.2609 [113,] -30484.1118 -6923.6792 [114,] 5287.7493 -30484.1118 [115,] 38602.6563 5287.7493 [116,] 51213.7185 38602.6563 [117,] -47258.6252 51213.7185 [118,] -124432.2239 -47258.6252 [119,] 34021.0692 -124432.2239 [120,] -2670.9085 34021.0692 [121,] -5209.2059 -2670.9085 [122,] -154593.4596 -5209.2059 [123,] -195608.0717 -154593.4596 [124,] -65111.2965 -195608.0717 [125,] 116979.0506 -65111.2965 [126,] 61498.3824 116979.0506 [127,] -72793.0986 61498.3824 [128,] -83739.3648 -72793.0986 [129,] 58076.0874 -83739.3648 [130,] -70893.9582 58076.0874 [131,] -64516.3751 -70893.9582 [132,] -79240.1214 -64516.3751 [133,] 147766.5629 -79240.1214 [134,] -24527.6230 147766.5629 [135,] 51466.0230 -24527.6230 [136,] -7173.8251 51466.0230 [137,] -35582.5990 -7173.8251 [138,] 69180.9646 -35582.5990 [139,] -22599.5075 69180.9646 [140,] -24376.7115 -22599.5075 [141,] 62913.8086 -24376.7115 [142,] -70722.6475 62913.8086 [143,] -42085.0054 -70722.6475 [144,] 48088.1086 -42085.0054 [145,] -12251.9365 48088.1086 [146,] -52482.0974 -12251.9365 [147,] 108795.0843 -52482.0974 [148,] 182937.3208 108795.0843 [149,] -72166.1508 182937.3208 [150,] 92660.8092 -72166.1508 [151,] -2308.4495 92660.8092 [152,] 44387.2132 -2308.4495 [153,] 41205.5935 44387.2132 [154,] -103700.9420 41205.5935 [155,] -23098.8464 -103700.9420 [156,] 41500.5171 -23098.8464 [157,] -47586.9472 41500.5171 [158,] 30618.6918 -47586.9472 [159,] 288586.6540 30618.6918 [160,] 91303.7015 288586.6540 [161,] 54683.3950 91303.7015 [162,] 94476.7178 54683.3950 [163,] -10707.7616 94476.7178 [164,] 26322.3494 -10707.7616 [165,] 61036.8960 26322.3494 [166,] -44563.8256 61036.8960 [167,] 67212.4307 -44563.8256 [168,] 16543.1990 67212.4307 [169,] 3291.0507 16543.1990 [170,] 175235.5645 3291.0507 [171,] 162661.7182 175235.5645 [172,] 79736.1003 162661.7182 [173,] -10445.5213 79736.1003 [174,] 269742.7892 -10445.5213 [175,] 81423.3615 269742.7892 [176,] -26043.8208 81423.3615 [177,] 81068.5625 -26043.8208 [178,] -59026.1680 81068.5625 [179,] 38551.5690 -59026.1680 [180,] -24132.5768 38551.5690 [181,] 82599.4785 -24132.5768 [182,] 151346.5140 82599.4785 [183,] 113730.9565 151346.5140 [184,] 152419.6456 113730.9565 [185,] 70942.5853 152419.6456 [186,] 54998.2251 70942.5853 [187,] -32827.0687 54998.2251 [188,] -61945.2561 -32827.0687 [189,] 98093.5329 -61945.2561 [190,] -192008.4879 98093.5329 [191,] -34138.0828 -192008.4879 [192,] -181707.6964 -34138.0828 [193,] -151739.3894 -181707.6964 [194,] -91898.8012 -151739.3894 [195,] 277768.0944 -91898.8012 [196,] -135608.0500 277768.0944 [197,] -188919.2089 -135608.0500 [198,] 102905.7018 -188919.2089 [199,] 15510.8803 102905.7018 [200,] -237591.8634 15510.8803 [201,] -79415.9118 -237591.8634 [202,] -127485.2911 -79415.9118 [203,] -20426.4917 -127485.2911 [204,] -163608.4908 -20426.4917 [205,] -95671.9868 -163608.4908 [206,] 34623.0762 -95671.9868 [207,] -632906.1831 34623.0762 [208,] 420778.2462 -632906.1831 [209,] 116210.1826 420778.2462 [210,] 113874.9593 116210.1826 [211,] -183602.4075 113874.9593 [212,] -139109.1972 -183602.4075 [213,] 83945.8618 -139109.1972 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -26431.6816 90945.6890 2 -152473.9555 -26431.6816 3 55710.4893 -152473.9555 4 -118660.0363 55710.4893 5 64123.7765 -118660.0363 6 -95510.1935 64123.7765 7 54659.5163 -95510.1935 8 69004.7594 54659.5163 9 -28123.8249 69004.7594 10 139124.4708 -28123.8249 11 118966.5259 139124.4708 12 128188.4557 118966.5259 13 -39721.7390 128188.4557 14 -89436.6614 -39721.7390 15 -31658.9483 -89436.6614 16 -49185.9588 -31658.9483 17 15337.1862 -49185.9588 18 -110325.6568 15337.1862 19 -26920.0787 -110325.6568 20 135674.2027 -26920.0787 21 -75019.3021 135674.2027 22 147479.8459 -75019.3021 23 82190.0877 147479.8459 24 110200.8293 82190.0877 25 -79175.1075 110200.8293 26 -80670.8124 -79175.1075 27 -8919.5067 -80670.8124 28 -132395.6533 -8919.5067 29 -30695.8873 -132395.6533 30 -151170.7823 -30695.8873 31 -14914.5713 -151170.7823 32 110968.5848 -14914.5713 33 -50521.9376 110968.5848 34 96203.1530 -50521.9376 35 64628.7423 96203.1530 36 23529.3743 64628.7423 37 -33446.4849 23529.3743 38 -43363.3649 -33446.4849 39 -100191.3284 -43363.3649 40 -98397.1124 -100191.3284 41 16325.1713 -98397.1124 42 -225453.3180 16325.1713 43 -23728.1457 -225453.3180 44 138133.4213 -23728.1457 45 -30470.8909 138133.4213 46 104975.4532 -30470.8909 47 40576.0430 104975.4532 48 -32644.3365 40576.0430 49 -6235.5921 -32644.3365 50 36555.0189 -6235.5921 51 -28379.6699 36555.0189 52 -75757.2988 -28379.6699 53 -3309.1513 -75757.2988 54 -86306.5507 -3309.1513 55 14306.0401 -86306.5507 56 59663.0798 14306.0401 57 38782.5595 59663.0798 58 143370.5497 38782.5595 59 -4567.2274 143370.5497 60 453.0063 -4567.2274 61 53549.4485 453.0063 62 18886.6483 53549.4485 63 -46220.9533 18886.6483 64 -44882.8762 -46220.9533 65 -4042.5790 -44882.8762 66 -72903.7388 -4042.5790 67 110985.6216 -72903.7388 68 -11476.5180 110985.6216 69 8544.5494 -11476.5180 70 86679.7326 8544.5494 71 1587.4990 86679.7326 72 -3960.9673 1587.4990 73 37939.0716 -3960.9673 74 35983.3656 37939.0716 75 -48957.3089 35983.3656 76 -19828.3947 -48957.3089 77 -30534.3341 -19828.3947 78 -78691.5335 -30534.3341 79 88608.4850 -78691.5335 80 71926.7697 88608.4850 81 51384.1134 71926.7697 82 46920.3948 51384.1134 83 -222250.1779 46920.3948 84 10756.1608 -222250.1779 85 176351.7003 10756.1608 86 85410.7913 176351.7003 87 2025.4305 85410.7913 88 -55265.3548 2025.4305 89 -52937.9954 -55265.3548 90 -407.4535 -52937.9954 91 20232.0154 -407.4535 92 83014.6718 20232.0154 93 19483.2485 83014.6718 94 -33675.7533 19483.2485 95 -115929.6101 -33675.7533 96 17910.1839 -115929.6101 97 -22968.5277 17910.1839 98 64799.9742 -22968.5277 99 47314.7810 64799.9742 100 -117985.4785 47314.7810 101 4516.1914 -117985.4785 102 -43857.0717 4516.1914 103 -33927.4872 -43857.0717 104 -206026.0388 -33927.4872 105 -293724.3211 -206026.0388 106 61755.6975 -293724.3211 107 79277.8501 61755.6975 108 -150.4262 79277.8501 109 18941.2861 -150.4262 110 55987.1307 18941.2861 111 -15217.2609 55987.1307 112 -6923.6792 -15217.2609 113 -30484.1118 -6923.6792 114 5287.7493 -30484.1118 115 38602.6563 5287.7493 116 51213.7185 38602.6563 117 -47258.6252 51213.7185 118 -124432.2239 -47258.6252 119 34021.0692 -124432.2239 120 -2670.9085 34021.0692 121 -5209.2059 -2670.9085 122 -154593.4596 -5209.2059 123 -195608.0717 -154593.4596 124 -65111.2965 -195608.0717 125 116979.0506 -65111.2965 126 61498.3824 116979.0506 127 -72793.0986 61498.3824 128 -83739.3648 -72793.0986 129 58076.0874 -83739.3648 130 -70893.9582 58076.0874 131 -64516.3751 -70893.9582 132 -79240.1214 -64516.3751 133 147766.5629 -79240.1214 134 -24527.6230 147766.5629 135 51466.0230 -24527.6230 136 -7173.8251 51466.0230 137 -35582.5990 -7173.8251 138 69180.9646 -35582.5990 139 -22599.5075 69180.9646 140 -24376.7115 -22599.5075 141 62913.8086 -24376.7115 142 -70722.6475 62913.8086 143 -42085.0054 -70722.6475 144 48088.1086 -42085.0054 145 -12251.9365 48088.1086 146 -52482.0974 -12251.9365 147 108795.0843 -52482.0974 148 182937.3208 108795.0843 149 -72166.1508 182937.3208 150 92660.8092 -72166.1508 151 -2308.4495 92660.8092 152 44387.2132 -2308.4495 153 41205.5935 44387.2132 154 -103700.9420 41205.5935 155 -23098.8464 -103700.9420 156 41500.5171 -23098.8464 157 -47586.9472 41500.5171 158 30618.6918 -47586.9472 159 288586.6540 30618.6918 160 91303.7015 288586.6540 161 54683.3950 91303.7015 162 94476.7178 54683.3950 163 -10707.7616 94476.7178 164 26322.3494 -10707.7616 165 61036.8960 26322.3494 166 -44563.8256 61036.8960 167 67212.4307 -44563.8256 168 16543.1990 67212.4307 169 3291.0507 16543.1990 170 175235.5645 3291.0507 171 162661.7182 175235.5645 172 79736.1003 162661.7182 173 -10445.5213 79736.1003 174 269742.7892 -10445.5213 175 81423.3615 269742.7892 176 -26043.8208 81423.3615 177 81068.5625 -26043.8208 178 -59026.1680 81068.5625 179 38551.5690 -59026.1680 180 -24132.5768 38551.5690 181 82599.4785 -24132.5768 182 151346.5140 82599.4785 183 113730.9565 151346.5140 184 152419.6456 113730.9565 185 70942.5853 152419.6456 186 54998.2251 70942.5853 187 -32827.0687 54998.2251 188 -61945.2561 -32827.0687 189 98093.5329 -61945.2561 190 -192008.4879 98093.5329 191 -34138.0828 -192008.4879 192 -181707.6964 -34138.0828 193 -151739.3894 -181707.6964 194 -91898.8012 -151739.3894 195 277768.0944 -91898.8012 196 -135608.0500 277768.0944 197 -188919.2089 -135608.0500 198 102905.7018 -188919.2089 199 15510.8803 102905.7018 200 -237591.8634 15510.8803 201 -79415.9118 -237591.8634 202 -127485.2911 -79415.9118 203 -20426.4917 -127485.2911 204 -163608.4908 -20426.4917 205 -95671.9868 -163608.4908 206 34623.0762 -95671.9868 207 -632906.1831 34623.0762 208 420778.2462 -632906.1831 209 116210.1826 420778.2462 210 113874.9593 116210.1826 211 -183602.4075 113874.9593 212 -139109.1972 -183602.4075 213 83945.8618 -139109.1972 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7ns801292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/8ns801292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/9g2p41292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/rcomp/tmp/10g2p41292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + 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") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/11uc841292011608.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/12xvps1292011608.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/134em41292011608.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/14f5l71292011608.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/rcomp/tmp/150n1u1292011608.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/rcomp/tmp/16exz31292011608.tab") + } > try(system("convert tmp/191aa1292011607.ps tmp/191aa1292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/2karv1292011607.ps tmp/2karv1292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/3karv1292011607.ps tmp/3karv1292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/4karv1292011607.ps tmp/4karv1292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/5vj8x1292011607.ps tmp/5vj8x1292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/6vj8x1292011607.ps tmp/6vj8x1292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/7ns801292011607.ps tmp/7ns801292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/8ns801292011607.ps tmp/8ns801292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/9g2p41292011607.ps tmp/9g2p41292011607.png",intern=TRUE)) character(0) > try(system("convert tmp/10g2p41292011607.ps tmp/10g2p41292011607.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 6.60 0.89 7.49