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(350 + ,1595 + ,17 + ,10 + ,60720 + ,319369 + ,143 + ,7022 + ,75 + ,0 + ,61 + ,5565 + ,47 + ,43 + ,94355 + ,493408 + ,38 + ,6243 + ,53 + ,39 + ,287 + ,601 + ,6 + ,5 + ,60720 + ,319210 + ,144 + ,19868 + ,198 + ,0 + ,268 + ,188 + ,11 + ,9 + ,77655 + ,381180 + ,67 + ,16471 + ,964 + ,0 + ,25 + ,7146 + ,105 + ,83 + ,134028 + ,297978 + ,367 + ,933 + ,14 + ,305 + ,418 + ,1135 + ,17 + ,20 + ,62285 + ,290476 + ,384 + ,5322 + ,80 + ,-9 + ,392 + ,450 + ,8 + ,6 + ,59325 + ,292136 + ,380 + ,11517 + ,205 + ,-1 + ,131 + ,34 + ,8 + ,7 + ,60630 + ,314353 + ,309 + ,14294 + ,3340 + ,0 + ,316 + ,133 + ,14 + ,4 + ,65990 + ,339445 + ,109 + ,9960 + ,1052 + ,0 + ,347 + ,119 + ,6 + ,3 + ,59118 + ,303677 + ,354 + ,17280 + ,872 + ,0 + ,68 + ,2053 + ,53 + ,44 + ,100423 + ,397144 + ,60 + ,3720 + ,96 + ,33 + ,70 + ,4036 + ,63 + ,85 + ,100269 + ,424898 + ,53 + ,3570 + ,56 + ,32 + ,147 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,0 + ,169 + ,4655 + ,393 + ,55 + ,60720 + ,341570 + ,105 + ,360 + ,30 + ,0 + ,248 + ,131 + ,11 + ,15 + ,61808 + ,308989 + ,338 + ,9908 + ,834 + ,0 + ,152 + ,1766 + ,73 + ,46 + ,60438 + ,305959 + ,349 + ,1451 + ,60 + ,0 + ,351 + ,312 + ,14 + ,3 + ,58598 + ,318690 + ,147 + ,8478 + ,381 + ,0 + ,372 + ,448 + ,40 + ,7 + ,59781 + ,323361 + ,132 + ,3084 + ,275 + ,0 + ,137 + ,115 + ,13 + ,10 + ,60945 + ,318903 + ,145 + ,9146 + ,1037 + ,0 + ,352 + ,60 + ,10 + ,7 + ,61124 + ,314049 + ,314 + ,11405 + ,1916 + ,0 + ,338 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,180 + ,0 + ,258 + ,364 + ,35 + ,26 + ,47705 + ,298466 + ,365 + ,2813 + ,271 + ,0 + ,343 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,183 + ,0 + ,398 + ,1442 + ,118 + ,239 + ,121173 + ,438493 + ,49 + ,2021 + ,165 + ,-3 + ,115 + ,1389 + ,9 + ,6 + ,57530 + ,378049 + ,71 + ,19783 + ,128 + ,6 + ,188 + ,149 + ,5 + ,6 + ,62041 + ,313332 + ,318 + ,22666 + ,761 + ,0 + ,379 + ,2212 + ,14 + ,10 + ,60720 + ,319864 + ,142 + ,8562 + ,54 + ,0 + ,48 + ,7489 + ,322 + ,149 + ,136996 + ,430866 + ,50 + ,717 + ,31 + ,70 + ,366 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,287 + ,0 + ,117 + ,402 + ,26 + ,8 + ,84990 + ,332743 + ,115 + ,5106 + ,330 + ,5 + ,412 + ,7419 + ,29 + ,21 + ,63255 + ,286963 + ,387 + ,2999 + ,12 + ,-5 + ,157 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,219 + ,0 + ,264 + ,307 + ,15 + ,10 + ,58856 + ,331955 + ,116 + ,8797 + ,430 + ,0 + ,54 + ,1134 + ,162 + ,63 + ,146216 + ,527021 + ,34 + ,2019 + ,288 + ,58 + ,64 + ,7561 + ,87 + ,77 + ,82425 + ,364304 + ,86 + ,1889 + ,22 + ,37 + ,113 + ,5131 + ,71 + ,44 + ,86111 + ,292154 + ,379 + ,1298 + ,18 + ,9 + ,148 + ,9 + ,6 + ,7 + ,60835 + ,314987 + ,300 + ,19165 + ,12458 + ,0 + ,425 + ,2264 + ,24 + ,20 + ,55174 + ,272713 + ,394 + ,3030 + ,32 + ,-31 + ,7 + ,40949 + ,481 + ,263 + ,129352 + ,1398893 + ,9 + ,2493 + ,29 + ,2865 + ,197 + ,4980 + ,91 + ,89 + ,113521 + ,409642 + ,55 + ,2304 + ,42 + ,0 + ,97 + ,272 + ,8 + ,9 + ,112995 + ,429112 + ,51 + ,28639 + ,844 + ,14 + ,168 + ,757 + ,19 + ,18 + ,45689 + ,382712 + ,66 + ,9616 + ,242 + ,0 + ,69 + ,1424 + ,94 + ,73 + ,131116 + ,374943 + ,77 + ,1861 + ,123 + ,32 + ,127 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,191 + ,0 + ,184 + ,378 + ,52 + ,16 + ,71873 + ,297413 + ,369 + ,1873 + ,258 + ,0 + ,401 + ,1061 + ,5 + ,3 + ,70415 + ,304252 + ,353 + ,20850 + ,98 + ,-4 + ,371 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,290 + ,0 + ,174 + ,203 + ,9 + ,9 + ,61938 + ,309596 + ,334 + ,12177 + ,541 + ,0 + ,46 + ,3689 + ,25 + ,25 + ,101792 + ,377305 + ,74 + ,7092 + ,48 + ,74 + ,90 + ,567 + ,11 + ,14 + ,103345 + ,386212 + ,65 + ,16928 + ,329 + ,20 + ,27 + ,27664 + ,282 + ,236 + ,74996 + ,657954 + ,24 + ,1624 + ,17 + ,247 + ,105 + ,1686 + ,40 + ,28 + ,68696 + ,368186 + ,81 + ,4205 + ,100 + ,11 + ,413 + ,1164 + ,17 + ,15 + ,57320 + ,269753 + ,395 + ,4103 + ,60 + ,-6 + ,26 + ,15824 + ,64 + ,56 + ,63838 + ,585715 + ,30 + ,6027 + ,24 + ,284 + ,16 + ,6575 + ,159 + ,63 + ,151352 + ,1071292 + ,11 + ,5480 + ,133 + ,819 + ,19 + ,37597 + ,248 + ,234 + ,37884 + ,992426 + ,12 + ,3195 + ,21 + ,652 + ,336 + ,314 + ,5 + ,4 + ,60720 + ,315371 + ,294 + ,23074 + ,367 + ,0 + ,110 + ,1932 + ,36 + ,20 + ,72954 + ,249898 + ,403 + ,1386 + ,26 + ,10 + ,86 + ,5692 + ,72 + ,68 + ,64239 + ,357312 + ,91 + ,2185 + ,28 + ,22 + ,294 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,257 + ,0 + ,322 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,177 + ,0 + ,94 + ,6533 + ,75 + ,29 + ,126910 + ,364839 + ,85 + ,2198 + ,25 + ,16 + ,400 + ,2126 + ,138 + ,41 + ,21001 + ,230621 + ,413 + ,222 + ,14 + ,-3 + ,150 + ,14 + ,1 + ,0 + ,61016 + ,315877 + ,161 + ,115877 + ,8041 + ,0 + ,107 + ,3192 + ,60 + ,48 + ,59831 + ,125390 + ,427 + ,-1243 + ,-23 + ,10 + ,166 + ,47 + ,3 + ,1 + ,60720 + ,314882 + ,305 + ,38294 + ,2432 + ,0 + ,47 + ,4891 + ,217 + ,49 + ,210568 + ,370837 + ,79 + ,787 + ,35 + ,73 + ,182 + ,531 + ,10 + ,11 + ,71641 + ,330068 + ,120 + ,13007 + ,245 + ,0 + ,245 + ,2332 + ,4 + ,2 + ,72680 + ,324385 + ,131 + ,31096 + ,53 + ,0 + ,330 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,277 + ,0 + ,10 + ,16275 + ,109 + ,129 + ,47818 + ,1443586 + ,7 + ,11409 + ,76 + ,1974 + ,429 + ,9252 + ,48 + ,32 + ,79510 + ,253537 + ,401 + ,1115 + ,6 + ,-129 + ,2 + ,29790 + ,309 + ,444 + ,81767 + ,4321023 + ,2 + ,13337 + ,138 + ,28328 + ,235 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,204 + ,0 + ,187 + ,2300 + ,29 + ,23 + ,64873 + ,352108 + ,97 + ,5245 + ,66 + ,0 + ,213 + ,267 + ,32 + ,11 + ,57640 + ,330059 + ,121 + ,4064 + ,487 + ,0 + ,296 + ,382 + ,21 + ,9 + ,59365 + ,340968 + ,107 + ,6713 + ,369 + ,0 + ,218 + ,397 + ,28 + ,7 + ,60646 + ,317736 + ,152 + ,4205 + ,297 + ,0 + ,28 + ,11920 + ,182 + ,73 + ,31080 + ,209458 + ,419 + ,52 + ,1 + ,238 + ,5 + ,54660 + ,302 + ,315 + ,126942 + ,1491348 + ,6 + ,4276 + ,24 + ,3040 + ,138 + ,5 + ,2 + ,1 + ,60720 + ,314887 + ,304 + ,57443 + ,24814 + ,0 + ,163 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,0 + ,219 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,199 + ,0 + ,125 + ,2647 + ,37 + ,17 + ,62920 + ,353058 + ,95 + ,4137 + ,58 + ,1 + ,216 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,201 + ,0 + ,196 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,263 + ,0 + ,128 + ,94 + ,7 + ,10 + ,60793 + ,314533 + ,308 + ,16362 + ,1212 + ,0 + ,223 + ,3 + ,2 + ,1 + ,60698 + ,315354 + ,296 + ,57677 + ,36505 + ,0 + ,165 + ,454 + ,51 + ,24 + ,63261 + ,302187 + ,357 + ,2004 + ,225 + ,0 + ,373 + ,227 + ,9 + ,7 + ,76644 + ,336639 + ,110 + ,15182 + ,603 + ,0 + ,324 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,245 + ,0 + ,251 + ,206 + ,9 + ,9 + ,53110 + ,296702 + ,372 + ,10745 + ,469 + ,0 + ,11 + ,25830 + ,115 + ,104 + ,245546 + ,1073089 + ,10 + ,7592 + ,34 + ,1910 + ,62 + ,1528 + ,38 + ,32 + ,60326 + ,146494 + ,426 + ,-1408 + ,-35 + ,38 + ,149 + ,271 + ,2 + ,2 + ,60720 + ,325249 + ,129 + ,62625 + ,463 + ,0 + ,270 + ,138 + ,6 + ,26 + ,69817 + ,331420 + ,117 + ,21903 + ,954 + ,0 + ,171 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,211 + ,0 + ,226 + ,278 + ,26 + ,22 + ,61404 + ,314922 + ,302 + ,4420 + ,413 + ,0 + ,276 + ,282 + ,30 + ,12 + ,62452 + ,320016 + ,141 + ,4001 + ,425 + ,0 + ,383 + ,571 + ,80 + ,14 + ,41477 + ,280398 + ,390 + ,1005 + ,141 + ,0 + ,93 + ,2253 + ,99 + ,68 + ,63593 + ,452469 + ,43 + ,2550 + ,112 + ,17 + ,302 + ,290 + ,2 + ,1 + ,58790 + ,301164 + ,360 + ,50582 + ,348 + ,0 + ,385 + ,78 + ,16 + ,12 + ,62700 + ,317330 + ,154 + ,7333 + ,1511 + ,0 + ,225 + ,20 + ,2 + ,1 + ,60805 + ,315576 + ,169 + ,57788 + ,5924 + ,0 + ,275 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,284 + ,0 + ,431 + ,18073 + ,180 + ,69 + ,27284 + ,-7170 + ,429 + ,-1151 + ,-11 + ,-366 + ,220 + ,1866 + ,29 + ,28 + ,56225 + ,322331 + ,135 + ,4218 + ,66 + ,0 + ,6 + ,42634 + ,156 + ,168 + ,157214 + ,1629616 + ,5 + ,9164 + ,34 + ,2964 + ,290 + ,249 + ,62 + ,13 + ,54323 + ,292754 + ,378 + ,1496 + ,373 + ,0 + ,129 + ,422 + ,28 + ,23 + ,57935 + ,318056 + ,150 + ,4216 + ,280 + ,0 + ,65 + ,2675 + ,79 + ,54 + ,59017 + ,355178 + ,93 + ,1964 + ,58 + ,36 + ,416 + ,965 + ,3 + ,2 + ,73490 + ,204325 + ,421 + ,1442 + ,4 + ,-8 + ,299 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,224 + ,0 + ,407 + ,621 + ,16 + ,11 + ,68005 + ,317046 + ,156 + ,7315 + ,188 + ,-5 + ,252 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,240 + ,0 + ,257 + ,365 + ,8 + ,7 + ,54820 + ,309560 + ,335 + ,13695 + ,300 + ,0 + ,63 + ,3122 + ,34 + ,25 + ,94670 + ,414462 + ,54 + ,6308 + ,69 + ,38 + ,21 + ,12988 + ,119 + ,73 + ,82340 + ,857217 + ,16 + ,5523 + ,51 + ,601 + ,38 + ,5336 + ,81 + ,80 + ,112477 + ,697458 + ,22 + ,6141 + ,93 + ,132 + ,51 + ,3160 + ,108 + ,90 + ,108094 + ,530670 + ,33 + ,3062 + ,105 + ,66 + ,386 + ,2489 + ,40 + ,29 + ,67804 + ,238125 + ,408 + ,953 + ,15 + ,0 + ,325 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,273 + ,0 + ,33 + ,6984 + ,45 + ,40 + ,80570 + ,741409 + ,19 + ,12031 + ,78 + ,195 + ,91 + ,2001 + ,52 + ,36 + ,95551 + ,393343 + ,62 + ,3718 + ,97 + ,19 + ,430 + ,15236 + ,90 + ,86 + ,77440 + ,372631 + ,78 + ,1918 + ,11 + ,-158 + ,259 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,239 + ,0 + ,382 + ,530 + ,60 + ,18 + ,73433 + ,317291 + ,155 + ,1955 + ,221 + ,0 + ,369 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,292 + ,0 + ,428 + ,35624 + ,223 + ,244 + ,157278 + ,306275 + ,347 + ,477 + ,3 + ,-70 + ,100 + ,1383 + ,11 + ,9 + ,73221 + ,317892 + ,151 + ,10717 + ,85 + ,13 + ,122 + ,875 + ,14 + ,12 + ,67000 + ,334280 + ,112 + ,9591 + ,153 + ,3 + ,198 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,259 + ,0 + ,289 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,254 + ,0 + ,307 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,222 + ,0 + ,389 + ,72 + ,3 + ,3 + ,60398 + ,314210 + ,312 + ,38070 + ,1591 + ,0 + ,210 + ,265 + ,31 + ,3 + ,60720 + ,306948 + ,346 + ,3450 + ,404 + ,194 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,186 + ,0 + ,277 + ,335 + ,33 + ,23 + ,64175 + ,320398 + ,140 + ,3648 + ,359 + ,0 + ,298 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,230 + ,0 + ,356 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,269 + ,0 + ,43 + ,4525 + ,226 + ,89 + ,93811 + ,501749 + ,37 + ,1335 + ,67 + ,101 + ,71 + ,3045 + ,58 + ,49 + ,27330 + ,202055 + ,422 + ,35 + ,1 + ,31 + ,133 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,196 + ,0 + ,212 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,0 + ,186 + ,638 + ,14 + ,11 + ,60370 + ,333210 + ,114 + ,9515 + ,209 + ,0 + ,262 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,283 + ,0 + ,378 + ,607 + ,11 + ,7 + ,60436 + ,322340 + ,134 + ,11122 + ,201 + ,0 + ,132 + ,1558 + ,63 + ,25 + ,55637 + ,369448 + ,80 + ,2690 + ,109 + ,0 + ,279 + ,1324 + ,117 + ,29 + ,67440 + ,291841 + ,381 + ,785 + ,69 + ,0 + ,376 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,280 + ,0 + ,267 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,0 + ,403 + ,611 + ,8 + ,6 + ,59190 + ,296919 + ,371 + ,12115 + ,159 + ,-4 + ,303 + ,0 + ,0 + ,1 + ,60720 + ,315380 + ,225 + ,0 + ,388 + ,923 + ,8 + ,6 + ,58620 + ,309038 + ,337 + ,13630 + ,118 + ,0 + ,420 + ,661 + ,3 + ,1 + ,65920 + ,246541 + ,405 + ,15514 + ,70 + ,-12 + ,114 + ,2397 + ,14 + ,13 + ,74020 + ,289513 + ,385 + ,6394 + ,37 + ,6 + ,305 + ,366 + ,53 + ,14 + ,61808 + ,344425 + ,101 + ,2725 + ,394 + ,0 + ,346 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,181 + ,0 + ,354 + ,135 + ,2 + ,1 + ,62065 + ,314210 + ,311 + ,57105 + ,847 + ,0 + ,124 + ,1659 + ,29 + ,25 + ,107577 + ,480382 + ,40 + ,9668 + ,169 + ,2 + ,136 + ,316 + ,9 + ,4 + ,60505 + ,315009 + ,299 + ,12779 + ,364 + ,0 + ,130 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,197 + ,0 + ,146 + ,309 + ,11 + ,9 + ,56535 + ,312878 + ,322 + ,10262 + ,365 + ,0 + ,314 + ,49 + ,8 + ,14 + ,64107 + ,322031 + ,137 + ,15254 + ,2500 + ,0 + ,241 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,200 + ,0 + ,34 + ,4519 + ,49 + ,42 + ,102129 + ,597793 + ,29 + ,8118 + ,88 + ,163 + ,384 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,0 + ,120 + ,837 + ,69 + ,21 + ,61262 + ,315688 + ,165 + ,1677 + ,138 + ,4 + ,72 + ,5119 + ,49 + ,38 + ,39039 + ,378525 + ,69 + ,3643 + ,35 + ,31 + ,79 + ,1280 + ,117 + ,52 + ,69465 + ,312378 + ,326 + ,960 + ,88 + ,26 + ,99 + ,2564 + ,22 + ,19 + ,130140 + ,403560 + ,59 + ,9253 + ,79 + ,13 + ,108 + ,2045 + ,20 + ,16 + ,97890 + ,510834 + ,36 + ,15542 + ,152 + ,10 + ,119 + ,2234 + ,16 + ,11 + ,77200 + ,214215 + ,418 + ,888 + ,6 + ,5 + ,102 + ,975 + ,23 + ,19 + ,90534 + ,235133 + ,410 + ,1528 + ,36 + ,13 + ,209 + ,1136 + ,32 + ,26 + ,48522 + ,343613 + ,103 + ,4488 + ,126 + ,0 + ,293 + ,453 + ,21 + ,19 + ,81125 + ,365959 + ,84 + ,7903 + ,366 + ,0 + ,211 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,210 + ,0 + ,227 + ,61 + ,14 + ,9 + ,60720 + ,314551 + ,307 + ,8182 + ,1874 + ,0 + ,234 + ,368 + ,17 + ,13 + ,61686 + ,303230 + ,355 + ,6072 + ,281 + ,0 + ,247 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,233 + ,0 + ,49 + ,4901 + ,84 + ,58 + ,121920 + ,469107 + ,41 + ,3204 + ,55 + ,67 + ,88 + ,540 + ,14 + ,13 + ,103960 + ,354228 + ,94 + ,11016 + ,286 + ,22 + ,329 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,253 + ,0 + ,160 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,213 + ,0 + ,271 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,286 + ,0 + ,306 + ,2 + ,9 + ,12 + ,60735 + ,315394 + ,174 + ,12822 + ,63057 + ,0 + ,172 + ,36 + ,4 + ,3 + ,61564 + ,312412 + ,325 + ,28103 + ,3146 + ,0 + ,402 + ,776 + ,9 + ,8 + ,64230 + ,333505 + ,113 + ,14834 + ,172 + ,-4 + ,4 + ,84738 + ,588 + ,428 + ,-26007 + ,223193 + ,415 + ,39 + ,0 + ,4145 + ,156 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,217 + ,0 + ,315 + ,3 + ,4 + ,8 + ,60761 + ,315656 + ,167 + ,28914 + ,37920 + ,0 + ,399 + ,529 + ,9 + ,7 + ,63870 + ,296261 + ,374 + ,10696 + ,182 + ,-3 + ,381 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,264 + ,0 + ,159 + ,405 + ,7 + ,5 + ,60845 + ,336425 + ,111 + ,19489 + ,337 + ,0 + ,181 + ,972 + ,34 + ,19 + ,71642 + ,359335 + ,87 + ,4686 + ,164 + ,0 + ,278 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,232 + ,0 + ,89 + ,2099 + ,64 + ,45 + ,106611 + ,308636 + ,339 + ,1697 + ,52 + ,20 + ,112 + ,3437 + ,47 + ,33 + ,48022 + ,158492 + ,425 + ,-883 + ,-12 + ,9 + ,201 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,243 + ,0 + ,344 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,182 + ,0 + ,22 + ,22330 + ,84 + ,67 + ,79801 + ,711969 + ,20 + ,6095 + ,23 + ,446 + ,291 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,258 + ,0 + ,269 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,221 + ,0 + ,164 + ,483 + ,27 + ,7 + ,60830 + ,306268 + ,348 + ,3936 + ,220 + ,0 + ,358 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,279 + ,0 + ,118 + ,2239 + ,21 + ,17 + ,88590 + ,442882 + ,47 + ,11566 + ,108 + ,5 + ,92 + ,2949 + ,41 + ,40 + ,82903 + ,378509 + ,70 + ,4354 + ,61 + ,18 + ,300 + ,0 + ,0 + ,9 + ,60720 + ,315380 + ,223 + ,0 + ,312 + ,365 + ,83 + ,32 + ,87192 + ,346611 + ,99 + ,1766 + ,402 + ,0 + ,183 + ,2461 + ,57 + ,25 + ,55792 + ,314289 + ,310 + ,2005 + ,46 + ,0 + ,18 + ,21950 + ,519 + ,265 + ,114337 + ,856956 + ,17 + ,1266 + ,30 + ,745 + ,126 + ,3294 + ,23 + ,23 + ,75832 + ,217193 + ,416 + ,748 + ,5 + ,1 + ,193 + ,141 + ,6 + ,7 + ,61630 + ,315366 + ,295 + ,19228 + ,819 + ,0 + ,390 + ,572 + ,16 + ,8 + ,58580 + ,307930 + ,342 + ,6746 + ,189 + ,-1 + ,23 + ,13326 + ,102 + ,103 + ,165548 + ,702380 + ,21 + ,4925 + ,38 + ,436 + ,50 + ,2284 + ,33 + ,41 + ,76403 + ,194493 + ,423 + ,-167 + ,-2 + ,66 + ,206 + ,10 + ,2 + ,5 + ,61656 + ,316155 + ,160 + ,58077 + ,11639 + ,0 + ,162 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,207 + ,0 + ,98 + ,1414 + ,198 + ,22 + ,70184 + ,330546 + ,119 + ,659 + ,92 + ,13 + ,75 + ,1975 + ,35 + ,34 + ,118881 + ,394510 + ,61 + ,5557 + ,98 + ,30 + ,134 + ,43 + ,3 + ,20 + ,60887 + ,312846 + ,323 + ,37615 + ,2613 + ,0 + ,180 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,190 + ,0 + ,337 + ,844 + ,78 + ,20 + ,60925 + ,296139 + ,375 + ,1233 + ,114 + ,0 + ,341 + ,304 + ,15 + ,12 + ,62969 + ,295580 + ,376 + ,6372 + ,314 + ,0 + ,154 + ,458 + ,11 + ,7 + ,58625 + ,297765 + ,368 + ,8888 + ,213 + ,0 + ,36 + ,18562 + ,155 + ,211 + ,102313 + ,377934 + ,72 + ,1148 + ,10 + ,143 + ,192 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,179 + ,0 + ,17 + ,7123 + ,109 + ,102 + ,288170 + ,638830 + ,28 + ,4026 + ,62 + ,769 + ,254 + ,622 + ,73 + ,38 + ,73007 + ,304376 + ,352 + ,1430 + ,168 + ,0 + ,374 + ,174 + ,1 + ,3 + ,64820 + ,307424 + ,344 + ,107424 + ,618 + ,0 + ,35 + ,2220 + ,22 + ,12 + ,301670 + ,644190 + ,26 + ,20190 + ,200 + ,154 + ,377 + ,121 + ,24 + ,11 + ,56178 + ,295370 + ,377 + ,3974 + ,788 + ,0 + ,32 + ,11819 + ,85 + ,70 + ,106113 + ,574339 + ,31 + ,4404 + ,32 + ,199 + ,260 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,235 + ,0 + ,274 + ,125 + ,12 + ,19 + ,60798 + ,310201 + ,332 + ,9183 + ,883 + ,0 + ,406 + ,1182 + ,13 + ,11 + ,70694 + ,327007 + ,124 + ,9770 + ,107 + ,-5 + ,116 + ,1503 + ,39 + ,35 + ,56364 + ,343466 + ,104 + ,3679 + ,95 + ,6 + ,320 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,274 + ,0 + ,367 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,248 + ,0 + ,288 + ,30 + ,4 + ,3 + ,62045 + ,318098 + ,149 + ,29525 + ,3912 + ,0 + ,414 + ,3310 + ,33 + ,32 + ,75230 + ,448243 + ,44 + ,7523 + ,75 + ,-6 + ,74 + ,554 + ,11 + ,9 + ,79285 + ,325738 + ,125 + ,11431 + ,227 + ,30 + ,250 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,234 + ,0 + ,224 + ,468 + ,21 + ,20 + ,60720 + ,312161 + ,328 + ,5341 + ,239 + ,0 + ,397 + ,4917 + ,42 + ,26 + ,52811 + ,243650 + ,406 + ,1039 + ,9 + ,-3 + ,417 + ,3256 + ,66 + ,23 + ,35250 + ,407159 + ,57 + ,3139 + ,64 + ,-8 + ,323 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,185 + ,0 + ,340 + ,125 + ,13 + ,10 + ,59734 + ,317698 + ,153 + ,9054 + ,944 + ,0 + ,239 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,195 + ,0 + ,326 + ,22 + ,2 + ,1 + ,60722 + ,312502 + ,324 + ,56251 + ,5154 + ,0 + ,261 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,262 + ,0 + ,357 + ,514 + ,52 + ,15 + ,78780 + ,322378 + ,133 + ,2353 + ,238 + ,0 + ,313 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,0 + ,360 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,270 + ,0 + ,272 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,285 + ,0 + ,333 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,276 + ,0 + ,355 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,272 + ,0 + ,362 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,268 + ,0 + ,37 + ,10327 + ,91 + ,74 + ,88577 + ,640273 + ,27 + ,4838 + ,43 + ,135 + ,42 + ,13718 + ,61 + ,56 + ,96448 + ,345783 + ,100 + ,2390 + ,11 + ,104 + ,57 + ,3748 + ,9 + ,6 + ,50350 + ,652925 + ,25 + ,50325 + ,121 + ,53 + ,30 + ,14416 + ,361 + ,139 + ,49857 + ,439798 + ,48 + ,664 + ,17 + ,220 + ,404 + ,1526 + ,25 + ,21 + ,69351 + ,278990 + ,391 + ,3160 + ,52 + ,-4 + ,80 + ,666 + ,171 + ,13 + ,117869 + ,339836 + ,108 + ,818 + ,210 + ,25 + ,422 + ,2844 + ,67 + ,74 + ,72683 + ,240897 + ,407 + ,610 + ,14 + ,-15 + ,297 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,229 + ,0 + ,318 + ,368 + ,17 + ,20 + ,61167 + ,297141 + ,370 + ,5714 + ,264 + ,0 + ,175 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,205 + ,0 + ,285 + ,333 + ,13 + ,9 + ,70811 + ,331323 + ,118 + ,10102 + ,394 + ,0 + ,339 + ,26 + ,1 + ,1 + ,60896 + ,313880 + ,315 + ,113880 + ,4328 + ,0 + ,243 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,237 + ,0 + ,167 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,193 + ,0 + ,411 + ,1303 + ,61 + ,57 + ,69863 + ,309422 + ,336 + ,1794 + ,84 + ,-5 + ,155 + ,20 + ,3 + ,2 + ,60938 + ,315245 + ,297 + ,38415 + ,5674 + ,0 + ,56 + ,2384 + ,97 + ,77 + ,61348 + ,405972 + ,58 + ,2123 + ,86 + ,54 + ,204 + ,203 + ,27 + ,12 + ,50804 + ,300962 + ,361 + ,3739 + ,498 + ,0 + ,178 + ,71 + ,1 + ,0 + ,60745 + ,316176 + ,159 + ,116176 + ,1646 + ,0 + ,295 + ,53 + ,14 + ,10 + ,59506 + ,302409 + ,356 + ,7315 + ,1932 + ,0 + ,396 + ,562 + ,13 + ,8 + ,58456 + ,283587 + ,389 + ,6430 + ,149 + ,-3 + ,393 + ,622 + ,7 + ,8 + ,60950 + ,263276 + ,398 + ,9039 + ,102 + ,-2 + ,135 + ,645 + ,5 + ,4 + ,60720 + ,312075 + ,329 + ,22415 + ,174 + ,0 + ,364 + ,1763 + ,24 + ,28 + ,61600 + ,308336 + ,340 + ,4514 + ,61 + ,0 + ,179 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,0 + ,231 + ,317 + ,9 + ,7 + ,63915 + ,298700 + ,363 + ,10967 + ,311 + ,0 + ,319 + ,1 + ,6 + ,5 + ,60719 + ,315372 + ,293 + ,19229 + ,84213 + ,0 + ,191 + ,275 + ,4 + ,7 + ,59500 + ,318745 + ,146 + ,29686 + ,431 + ,0 + ,321 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,246 + ,0 + ,111 + ,936 + ,30 + ,22 + ,67939 + ,408881 + ,56 + ,6963 + ,223 + ,9 + ,29 + ,8568 + ,37 + ,34 + ,32168 + ,786690 + ,18 + ,15856 + ,68 + ,223 + ,423 + ,11528 + ,70 + ,68 + ,-14545 + ,-83265 + ,431 + ,-4047 + ,-25 + ,-25 + ,202 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,250 + ,0 + ,232 + ,738 + ,10 + ,9 + ,60720 + ,321376 + ,139 + ,12138 + ,164 + ,0 + ,266 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,281 + ,0 + ,405 + ,592 + ,7 + ,7 + ,64270 + ,276898 + ,392 + ,10985 + ,130 + ,-4 + ,176 + ,126 + ,10 + ,8 + ,60951 + ,315547 + ,171 + ,11555 + ,921 + ,0 + ,236 + ,2 + ,1 + ,0 + ,60743 + ,315487 + ,172 + ,115487 + ,52257 + ,0 + ,142 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,220 + ,0 + ,13 + ,18014 + ,239 + ,393 + ,-1710 + ,1405225 + ,8 + ,5043 + ,67 + ,943 + ,281 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,256 + ,0 + ,328 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,0 + ,9 + ,37704 + ,452 + ,228 + ,60448 + ,983660 + ,14 + ,1734 + ,21 + ,2167 + ,207 + ,63 + ,4 + ,2 + ,65688 + ,318574 + ,148 + ,29643 + ,1891 + ,0 + ,73 + ,1431 + ,34 + ,35 + ,106885 + ,310768 + ,330 + ,3258 + ,77 + ,30 + ,242 + ,94 + ,3 + ,3 + ,61360 + ,312887 + ,321 + ,37629 + ,1200 + ,0 + ,229 + ,192 + ,4 + ,2 + ,65276 + ,312339 + ,327 + ,28085 + ,586 + ,0 + ,317 + ,32 + ,1 + ,0 + ,59988 + ,314964 + ,301 + ,114964 + ,3603 + ,0 + ,44 + ,6869 + ,105 + ,88 + ,117520 + ,379983 + ,68 + ,1714 + ,26 + ,90 + ,349 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,176 + ,0 + ,140 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,0 + ,361 + ,1 + ,4 + ,4 + ,60722 + ,315398 + ,173 + ,28849 + ,172235 + ,0 + ,203 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,247 + ,0 + ,255 + ,2328 + ,94 + ,57 + ,82732 + ,253588 + ,400 + ,570 + ,23 + ,0 + ,308 + ,209 + ,13 + ,19 + ,64016 + ,316647 + ,157 + ,8973 + ,558 + ,0 + ,195 + ,28 + ,2 + ,3 + ,60890 + ,315688 + ,166 + ,57844 + ,4148 + ,0 + ,280 + ,176 + ,28 + ,8 + ,68136 + ,310670 + ,331 + ,3952 + ,629 + ,0 + ,415 + ,1920 + ,15 + ,11 + ,79420 + ,165404 + ,424 + ,-2306 + ,-18 + ,-6 + ,3 + ,87550 + ,458 + ,412 + ,153198 + ,4111912 + ,3 + ,8541 + ,45 + ,17936 + ,151 + ,520 + ,10 + ,7 + ,58650 + ,291650 + ,383 + ,9165 + ,176 + ,0 + ,304 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,231 + ,0 + ,421 + ,1013 + ,2 + ,1 + ,69770 + ,253468 + ,402 + ,26734 + ,53 + ,-14 + ,217 + ,15 + ,8 + ,6 + ,60831 + ,315688 + ,164 + ,14461 + ,7687 + ,0 + ,353 + ,587 + ,12 + ,10 + ,59595 + ,325699 + ,126 + ,10475 + ,214 + ,0 + ,58 + ,5371 + ,66 + ,47 + ,87720 + ,446211 + ,45 + ,3730 + ,46 + ,49 + ,359 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,251 + ,0 + ,77 + ,1012 + ,47 + ,29 + ,114768 + ,368078 + ,82 + ,3576 + ,166 + ,26 + ,228 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,215 + ,0 + ,83 + ,876 + ,59 + ,33 + ,138971 + ,352850 + ,96 + ,2591 + ,174 + ,24 + ,1 + ,162556 + ,1081 + ,807 + ,213118 + ,6282154 + ,1 + ,5626 + ,37 + ,56289 + ,14 + ,43556 + ,247 + ,190 + ,32648 + ,227132 + ,414 + ,110 + ,1 + ,929 + ,395 + ,3425 + ,12 + ,9 + ,83620 + ,283910 + ,388 + ,6993 + ,25 + ,-3 + ,78 + ,810 + ,43 + ,11 + ,74015 + ,236761 + ,409 + ,855 + ,45 + ,26 + ,190 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,184 + ,0 + ,311 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,244 + ,0 + ,39 + ,2365 + ,79 + ,83 + ,191778 + ,550608 + ,32 + ,4438 + ,148 + ,120 + ,101 + ,1261 + ,33 + ,31 + ,76114 + ,307528 + ,343 + ,3258 + ,85 + ,13 + ,256 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,238 + ,0 + ,123 + ,1585 + ,107 + ,112 + ,93099 + ,355864 + ,92 + ,1457 + ,98 + ,2 + ,24 + ,16189 + ,295 + ,290 + ,116384 + ,358589 + ,89 + ,538 + ,10 + ,379 + ,284 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,242 + ,0 + ,273 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,265 + ,0 + ,121 + ,10579 + ,92 + ,92 + ,110309 + ,375195 + ,76 + ,1904 + ,17 + ,4 + ,214 + ,474 + ,20 + ,10 + ,61977 + ,288985 + ,386 + ,4449 + ,188 + ,0 + ,348 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,178 + ,0 + ,87 + ,3642 + ,86 + ,100 + ,90262 + ,458343 + ,42 + ,3004 + ,71 + ,22 + ,170 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,0 + ,200 + ,472 + ,20 + ,19 + ,80045 + ,269587 + ,397 + ,3479 + ,148 + ,0 + ,158 + ,98 + ,2 + ,2 + ,61490 + ,315236 + ,298 + ,57618 + ,1179 + ,0 + ,427 + ,3999 + ,30 + ,21 + ,51252 + ,42754 + ,428 + ,-5242 + ,-39 + ,-58 + ,233 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,216 + ,0 + ,345 + ,621 + ,48 + ,13 + ,60720 + ,308256 + ,341 + ,2255 + ,174 + ,0 + ,375 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,266 + ,0 + ,173 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,189 + ,0 + ,301 + ,30 + ,4 + ,4 + ,60798 + ,313164 + ,320 + ,28291 + ,3808 + ,0 + ,368 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,291 + ,0 + ,106 + ,746 + ,23 + ,29 + ,71561 + ,269661 + ,396 + ,3029 + ,93 + ,10 + ,208 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,228 + ,0 + ,282 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,255 + ,0 + ,249 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,260 + ,0 + ,240 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,194 + ,0 + ,52 + ,4150 + ,150 + ,136 + ,134759 + ,518365 + ,35 + ,2122 + ,77 + ,66 + ,55 + ,4658 + ,158 + ,114 + ,156608 + ,233773 + ,411 + ,214 + ,7 + ,58 + ,84 + ,814 + ,18 + ,15 + ,39625 + ,301881 + ,358 + ,5660 + ,125 + ,23 + ,143 + ,1002 + ,15 + ,11 + ,56750 + ,298568 + ,364 + ,6571 + ,98 + ,0 + ,205 + ,496 + ,17 + ,12 + ,87390 + ,325479 + ,128 + ,7381 + ,253 + ,0 + ,141 + ,389 + ,14 + ,8 + ,58990 + ,325506 + ,127 + ,8965 + ,322 + ,0 + ,12 + ,12679 + ,110 + ,122 + ,48020 + ,984885 + ,13 + ,7135 + ,62 + ,1871 + ,177 + ,400 + ,152 + ,113 + ,60720 + ,313267 + ,319 + ,745 + ,283 + ,0 + ,238 + ,53 + ,6 + ,6 + ,60349 + ,315793 + ,163 + ,19299 + ,2197 + ,0 + ,331 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,278 + ,0 + ,426 + ,5109 + ,97 + ,82 + ,67038 + ,215362 + ,417 + ,158 + ,3 + ,-31 + ,104 + ,576 + ,18 + ,8 + ,113761 + ,314073 + ,313 + ,6337 + ,198 + ,12 + ,409 + ,438 + ,11 + ,9 + ,58320 + ,298096 + ,366 + ,8918 + ,224 + ,-5 + ,292 + ,165 + ,24 + ,12 + ,62841 + ,325176 + ,130 + ,5216 + ,759 + ,0 + ,40 + ,4069 + ,145 + ,131 + ,79804 + ,207393 + ,420 + ,51 + ,2 + ,119 + ,185 + ,23 + ,5 + ,5 + ,62555 + ,314806 + ,306 + ,22961 + ,4974 + ,0 + ,59 + ,1285 + ,107 + ,51 + ,99489 + ,341340 + ,106 + ,1321 + ,110 + ,49 + ,82 + ,4677 + ,75 + ,66 + ,90131 + ,426280 + ,52 + ,3017 + ,48 + ,24 + ,15 + ,24811 + ,505 + ,280 + ,95350 + ,929118 + ,15 + ,1444 + ,29 + ,822 + ,253 + ,167 + ,2 + ,1 + ,64245 + ,307322 + ,345 + ,53661 + ,642 + ,0 + ,342 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,187 + ,0 + ,45 + ,4628 + ,62 + ,39 + ,69159 + ,387475 + ,64 + ,3024 + ,41 + ,85 + ,189 + ,226 + ,5 + ,5 + ,65745 + ,291787 + ,382 + ,18357 + ,406 + ,0 + ,76 + ,1765 + ,20 + ,20 + ,77623 + ,247060 + ,404 + ,2353 + ,27 + ,28 + ,145 + ,460 + ,15 + ,15 + ,63346 + ,329784 + ,122 + ,8652 + ,282 + ,0 + ,144 + ,36 + ,3 + ,4 + ,60894 + ,315834 + ,162 + ,38611 + ,3215 + ,0 + ,408 + ,989 + ,11 + ,12 + ,58930 + ,304555 + ,350 + ,9505 + ,106 + ,-5 + ,96 + ,3055 + ,43 + ,30 + ,60247 + ,376641 + ,75 + ,4108 + ,58 + ,14 + ,265 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,288 + ,0 + ,161 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,212 + ,0 + ,246 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,236 + ,0 + ,237 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,202 + ,0 + ,230 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,214 + ,0 + ,66 + ,13253 + ,947 + ,251 + ,90829 + ,357760 + ,90 + ,167 + ,12 + ,35 + ,244 + ,24 + ,8 + ,7 + ,59818 + ,315637 + ,168 + ,14455 + ,4772 + ,0 + ,153 + ,0 + ,0 + ,5 + ,60720 + ,315380 + ,208 + ,0 + ,199 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,261 + ,0 + ,222 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,198 + ,0 + ,363 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,289 + ,0 + ,387 + ,131 + ,6 + ,3 + ,59661 + ,327071 + ,123 + ,21179 + ,969 + ,0 + ,85 + ,514 + ,15 + ,15 + ,102725 + ,377516 + ,73 + ,11834 + ,346 + ,23 + ,103 + ,3366 + ,80 + ,55 + ,108479 + ,299243 + ,362 + ,1241 + ,29 + ,12 + ,60 + ,9327 + ,101 + ,85 + ,87419 + ,387699 + ,63 + ,1858 + ,20 + ,44 + ,309 + ,384 + ,22 + ,17 + ,54683 + ,309836 + ,333 + ,4993 + ,286 + ,0 + ,20 + ,17821 + ,373 + ,277 + ,122844 + ,444477 + ,46 + ,655 + ,14 + ,643 + ,391 + ,397 + ,10 + ,10 + ,62710 + ,322327 + ,136 + ,12233 + ,308 + ,-1 + ,139 + ,897 + ,5 + ,4 + ,60720 + ,314913 + ,303 + ,22983 + ,128 + ,0 + ,370 + ,218 + ,21 + ,7 + ,60720 + ,315553 + ,170 + ,5503 + ,531 + ,0 + ,31 + ,3369 + ,28 + ,26 + ,87161 + ,688779 + ,23 + ,17456 + ,145 + ,217 + ,109 + ,5702 + ,61 + ,47 + ,101481 + ,321896 + ,138 + ,1998 + ,21 + ,10 + ,41 + ,8636 + ,855 + ,203 + ,128294 + ,301607 + ,359 + ,119 + ,12 + ,108 + ,215 + ,534 + ,7 + ,5 + ,62620 + ,304485 + ,351 + ,14926 + ,196 + ,0 + ,332 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,0 + ,310 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,227 + ,0 + ,410 + ,726 + ,3 + ,3 + ,69980 + ,231861 + ,412 + ,10620 + ,44 + ,-5 + ,81 + ,1380 + ,26 + ,29 + ,60982 + ,347385 + ,98 + ,5669 + ,107 + ,25 + ,365 + ,180 + ,11 + ,9 + ,59635 + ,316386 + ,158 + ,10581 + ,647 + ,0 + ,53 + ,7285 + ,115 + ,97 + ,188873 + ,491303 + ,39 + ,2533 + ,40 + ,59 + ,67 + ,880 + ,74 + ,15 + ,80791 + ,261216 + ,399 + ,827 + ,70 + ,34 + ,335 + ,1 + ,2 + ,1 + ,60727 + ,315388 + ,175 + ,57694 + ,117743 + ,0 + ,380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,271 + ,0 + ,334 + ,96 + ,5 + ,4 + ,60379 + ,313729 + ,316 + ,22746 + ,1188 + ,0 + ,95 + ,1889 + ,45 + ,43 + ,37527 + ,358649 + ,88 + ,3526 + ,84 + ,14 + ,8 + ,45187 + ,353 + ,267 + ,234817 + ,1926517 + ,4 + ,4891 + ,38 + ,2854 + ,221 + ,288 + ,4 + ,3 + ,60510 + ,296656 + ,373 + ,24164 + ,335 + ,0 + ,419 + ,1270 + ,26 + ,24 + ,69206 + ,275311 + ,393 + ,2897 + ,59 + ,-10 + ,424 + ,6526 + ,26 + ,20 + ,55830 + ,-42143 + ,430 + ,-9313 + ,-37 + ,-28 + ,327 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,252 + ,0 + ,263 + ,226 + ,12 + ,13 + ,72835 + ,343929 + ,102 + ,11994 + ,636 + ,0 + ,394 + ,694 + ,8 + ,6 + ,68060 + ,367655 + ,83 + ,20957 + ,241 + ,-2 + ,286 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,241 + ,0 + ,283 + ,249 + ,72 + ,26 + ,56726 + ,313491 + ,317 + ,1576 + ,455 + ,0) + ,dim=c(10 + ,431) + ,dimnames=list(c('Rank' + ,'Costs' + ,'Trades' + ,'Orders' + ,'Dividends' + ,'Wealth' + ,'WRank' + ,'Profit/Trades' + ,'Profit/Cost' + ,'CCScore ') + ,1:431)) > y <- array(NA,dim=c(10,431),dimnames=list(c('Rank','Costs','Trades','Orders','Dividends','Wealth','WRank','Profit/Trades','Profit/Cost','CCScore '),1:431)) > 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 = 'Do not include Seasonal Dummies' > par1 = '10' > #'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 Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > 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 CCScore\r Rank Costs Trades Orders Dividends Wealth WRank 1 0 350 1595 17 10 60720 319369 143 2 39 61 5565 47 43 94355 493408 38 3 0 287 601 6 5 60720 319210 144 4 0 268 188 11 9 77655 381180 67 5 305 25 7146 105 83 134028 297978 367 6 -9 418 1135 17 20 62285 290476 384 7 -1 392 450 8 6 59325 292136 380 8 0 131 34 8 7 60630 314353 309 9 0 316 133 14 4 65990 339445 109 10 0 347 119 6 3 59118 303677 354 11 33 68 2053 53 44 100423 397144 60 12 32 70 4036 63 85 100269 424898 53 13 4655 147 0 0 0 60720 315380 203 14 131 393 55 60720 341570 105 360 30 15 1766 11 15 61808 308989 338 9908 834 16 312 73 46 60438 305959 349 1451 60 17 448 14 3 58598 318690 147 8478 381 18 115 40 7 59781 323361 132 3084 275 19 60 13 10 60945 318903 145 9146 1037 20 0 10 7 61124 314049 314 11405 1916 21 26 0 0 60720 315380 180 0 258 22 0 47705 298466 365 2813 271 0 343 23 438493 60720 315380 183 0 398 1442 118 24 378049 49 2021 165 -3 115 1389 9 25 313332 71 19783 128 6 188 149 5 26 319864 318 22666 761 0 379 2212 14 27 430866 142 8562 54 0 48 7489 322 28 315380 50 717 31 70 366 0 0 29 5106 287 0 117 402 26 8 84990 30 2999 330 5 412 7419 29 21 63255 31 0 12 -5 157 0 0 0 60720 32 0 264 307 15 10 58856 331955 116 33 58 54 1134 162 63 146216 527021 34 34 37 64 7561 87 77 82425 364304 86 35 9 113 5131 71 44 86111 292154 379 36 0 148 9 6 7 60835 314987 300 37 -31 425 2264 24 20 55174 272713 394 38 2865 7 40949 481 263 129352 1398893 9 39 0 197 4980 91 89 113521 409642 55 40 14 97 272 8 9 112995 429112 51 41 0 168 757 19 18 45689 382712 66 42 32 69 1424 94 73 131116 374943 77 43 378 127 0 0 0 60720 315380 191 44 1061 52 16 71873 297413 369 1873 258 45 0 5 3 70415 304252 353 20850 98 46 9 0 0 60720 315380 290 0 174 47 25 61938 309596 334 12177 541 0 46 48 14 101792 377305 74 7092 48 74 90 49 236 103345 386212 65 16928 329 20 27 50 28 74996 657954 24 1624 17 247 105 51 15 68696 368186 81 4205 100 11 413 52 56 57320 269753 395 4103 60 -6 26 53 63 63838 585715 30 6027 24 284 16 54 234 151352 1071292 11 5480 133 819 19 55 4 37884 992426 12 3195 21 652 336 56 20 60720 315371 294 23074 367 0 110 57 68 72954 249898 403 1386 26 10 86 58 0 64239 357312 91 2185 28 22 294 59 315380 60720 315380 257 0 322 0 0 60 2198 177 0 94 6533 75 29 126910 61 222 25 16 400 2126 138 41 21001 62 115877 14 -3 150 14 1 0 61016 63 -1243 8041 0 107 3192 60 48 59831 64 38294 -23 10 166 47 3 1 60720 65 787 2432 0 47 4891 217 49 210568 66 13007 35 73 182 531 10 11 71641 67 31096 245 0 245 2332 4 2 72680 68 0 53 0 330 0 0 0 60720 69 1974 10 16275 109 129 47818 1443586 7 70 -129 429 9252 48 32 79510 253537 401 71 28328 2 29790 309 444 81767 4321023 2 72 2300 235 0 0 0 60720 315380 204 73 267 29 23 64873 352108 97 5245 66 74 382 32 11 57640 330059 121 4064 487 75 397 21 9 59365 340968 107 6713 369 76 11920 28 7 60646 317736 152 4205 297 77 54660 182 73 31080 209458 419 52 1 78 5 302 315 126942 1491348 6 4276 24 79 0 2 1 60720 314887 304 57443 24814 80 2 0 0 60720 315380 218 0 219 81 353058 60720 315380 199 0 125 2647 37 82 315380 95 4137 58 1 216 0 0 83 0 201 0 196 0 0 0 60720 84 0 128 94 7 10 60793 314533 308 85 0 223 3 2 1 60698 315354 296 86 0 165 454 51 24 63261 302187 357 87 0 373 227 9 7 76644 336639 110 88 206 324 0 0 0 60720 315380 245 89 25830 9 9 53110 296702 372 10745 469 90 1528 115 104 245546 1073089 10 7592 34 91 271 38 32 60326 146494 426 -1408 -35 92 138 2 2 60720 325249 129 62625 463 93 0 6 26 69817 331420 117 21903 954 94 22 0 0 60720 315380 211 0 226 95 12 61404 314922 302 4420 413 0 276 96 14 62452 320016 141 4001 425 0 383 97 68 41477 280398 390 1005 141 0 93 98 1 63593 452469 43 2550 112 17 302 99 12 58790 301164 360 50582 348 0 385 100 1 62700 317330 154 7333 1511 0 225 101 0 60805 315576 169 57788 5924 0 275 102 -7170 60720 315380 284 0 431 18073 180 103 322331 429 -1151 -11 -366 220 1866 29 104 1629616 135 4218 66 0 6 42634 156 105 292754 5 9164 34 2964 290 249 62 106 318056 378 1496 373 0 129 422 28 107 355178 150 4216 280 0 65 2675 79 108 204325 93 1964 58 36 416 965 3 109 315380 421 1442 4 -8 299 0 0 110 7315 224 0 407 621 16 11 68005 111 0 188 -5 252 0 0 0 60720 112 0 257 365 8 7 54820 309560 335 113 38 63 3122 34 25 94670 414462 54 114 601 21 12988 119 73 82340 857217 16 115 132 38 5336 81 80 112477 697458 22 116 66 51 3160 108 90 108094 530670 33 117 0 386 2489 40 29 67804 238125 408 118 6984 325 0 0 0 60720 315380 273 119 2001 45 40 80570 741409 19 12031 78 120 15236 52 36 95551 393343 62 3718 97 121 0 90 86 77440 372631 78 1918 11 122 18 0 0 60720 315380 239 0 382 123 0 73433 317291 155 1955 221 0 369 124 306275 60720 315380 292 0 428 35624 223 125 317892 347 477 3 -70 100 1383 11 126 334280 151 10717 85 13 122 875 14 127 315380 112 9591 153 3 198 0 0 128 0 259 0 289 0 0 0 60720 129 72 307 0 0 0 60720 315380 222 130 265 3 3 60398 314210 312 38070 1591 131 0 31 3 60720 306948 346 3450 404 132 64175 0 60720 315380 186 0 277 335 133 60720 320398 140 3648 359 0 298 0 134 269 315380 230 0 356 0 0 0 135 67 0 43 4525 226 89 93811 501749 136 1 101 71 3045 58 49 27330 202055 137 212 31 133 0 0 0 60720 315380 138 14 0 0 0 60720 315380 206 0 139 0 11 60370 333210 114 9515 209 0 140 60436 0 60720 315380 283 0 378 607 141 55637 322340 134 11122 201 0 132 1558 142 67440 369448 80 2690 109 0 279 1324 143 60720 291841 381 785 69 0 376 0 144 282 315380 280 0 267 0 0 0 145 159 0 403 611 8 6 59190 296919 146 388 -4 303 0 0 1 60720 315380 147 420 923 8 6 58620 309038 337 13630 148 114 661 3 1 65920 246541 405 15514 149 305 2397 14 13 74020 289513 385 6394 150 346 366 53 14 61808 344425 101 2725 151 2 0 0 0 60720 315380 181 0 152 29 1 62065 314210 311 57105 847 0 153 9 25 107577 480382 40 9668 169 2 154 0 4 60505 315009 299 12779 364 0 155 56535 0 60720 315380 197 0 146 309 156 64107 312878 322 10262 365 0 314 49 157 60720 322031 137 15254 2500 0 241 0 158 29 315380 200 0 34 4519 49 42 159 267 8118 88 163 384 0 0 0 160 138 0 120 837 69 21 61262 315688 161 35 4 72 5119 49 38 39039 378525 162 88 31 79 1280 117 52 69465 312378 163 79 26 99 2564 22 19 130140 403560 164 152 13 108 2045 20 16 97890 510834 165 6 10 119 2234 16 11 77200 214215 166 36 5 102 975 23 19 90534 235133 167 126 13 209 1136 32 26 48522 343613 168 366 0 293 453 21 19 81125 365959 169 227 0 211 0 0 0 60720 315380 170 234 61 14 9 60720 314551 307 8182 171 247 368 17 13 61686 303230 355 6072 172 84 0 0 0 60720 315380 233 0 173 14 58 121920 469107 41 3204 55 67 174 0 13 103960 354228 94 11016 286 22 175 60720 0 60720 315380 253 0 160 0 176 286 315380 213 0 271 0 0 0 177 63057 0 306 2 9 12 60735 315394 178 3146 0 172 36 4 3 61564 312412 179 172 0 402 776 9 8 64230 333505 180 0 -4 4 84738 588 428 -26007 223193 181 315 4145 156 0 0 0 60720 315380 182 399 3 4 8 60761 315656 167 28914 183 381 529 9 7 63870 296261 374 10696 184 7 0 0 0 60720 315380 264 0 185 34 5 60845 336425 111 19489 337 0 186 0 19 71642 359335 87 4686 164 0 187 106611 0 60720 315380 232 0 89 2099 188 48022 308636 339 1697 52 20 112 3437 189 60720 158492 425 -883 -12 9 201 0 190 182 315380 243 0 344 0 0 0 191 23 0 22 22330 84 67 79801 711969 192 269 446 291 0 0 0 60720 315380 193 27 0 0 0 60720 315380 221 0 194 0 7 60830 306268 348 3936 220 0 195 88590 0 60720 315380 279 0 118 2239 196 82903 442882 47 11566 108 5 92 2949 197 60720 378509 70 4354 61 18 300 0 198 99 315380 223 0 312 365 83 32 199 310 1766 402 0 183 2461 57 25 200 17 2005 46 0 18 21950 519 265 201 416 1266 30 745 126 3294 23 23 202 295 748 5 1 193 141 6 7 203 342 19228 819 0 390 572 16 8 204 21 6746 189 -1 23 13326 102 103 205 423 4925 38 436 50 2284 33 41 206 160 -167 -2 66 206 10 2 5 207 207 58077 11639 0 162 0 0 0 208 92 0 98 1414 198 22 70184 330546 209 98 13 75 1975 35 34 118881 394510 210 2613 30 134 43 3 20 60887 312846 211 337 0 180 0 0 0 60720 315380 212 341 844 78 20 60925 296139 375 1233 213 154 304 15 12 62969 295580 376 6372 214 36 458 11 7 58625 297765 368 8888 215 192 18562 155 211 102313 377934 72 1148 216 109 0 0 0 60720 315380 179 0 217 73 102 288170 638830 28 4026 62 769 218 1 38 73007 304376 352 1430 168 0 219 22 3 64820 307424 344 107424 618 0 220 24 12 301670 644190 26 20190 200 154 221 85 11 56178 295370 377 3974 788 0 222 0 70 106113 574339 31 4404 32 199 223 60798 0 60720 315380 235 0 274 125 224 70694 310201 332 9183 883 0 406 1182 225 56364 327007 124 9770 107 -5 116 1503 226 60720 343466 104 3679 95 6 320 0 227 248 315380 274 0 367 0 0 0 228 3912 0 288 30 4 3 62045 318098 229 75 0 414 3310 33 32 75230 448243 230 227 -6 74 554 11 9 79285 325738 231 224 30 250 0 0 0 60720 315380 232 397 468 21 20 60720 312161 328 5341 233 417 4917 42 26 52811 243650 406 1039 234 323 3256 66 23 35250 407159 57 3139 235 13 0 0 0 60720 315380 185 0 236 0 10 59734 317698 153 9054 944 0 237 60722 0 60720 315380 195 0 326 22 238 60720 312502 324 56251 5154 0 261 0 239 133 315380 262 0 357 514 52 15 240 226 2353 238 0 313 0 0 0 241 272 0 360 0 0 0 60720 315380 242 0 0 0 0 60720 315380 285 0 243 60720 0 60720 315380 276 0 355 0 244 268 315380 272 0 362 0 0 0 245 43 0 37 10327 91 74 88577 640273 246 11 135 42 13718 61 56 96448 345783 247 121 104 57 3748 9 6 50350 652925 248 17 53 30 14416 361 139 49857 439798 249 52 220 404 1526 25 21 69351 278990 250 210 -4 80 666 171 13 117869 339836 251 14 25 422 2844 67 74 72683 240897 252 318 -15 297 0 0 0 60720 315380 253 175 368 17 20 61167 297141 370 5714 254 13 0 0 0 60720 315380 205 0 255 1 9 70811 331323 118 10102 394 0 256 0 1 60896 313880 315 113880 4328 0 257 60720 0 60720 315380 237 0 167 0 258 336 315380 193 0 411 1303 61 57 259 297 1794 84 -5 155 20 3 2 260 58 38415 5674 0 56 2384 97 77 261 361 2123 86 54 204 203 27 12 262 159 3739 498 0 178 71 1 0 263 356 116176 1646 0 295 53 14 10 264 389 7315 1932 0 396 562 13 8 265 398 6430 149 -3 393 622 7 8 266 329 9039 102 -2 135 645 5 4 267 340 22415 174 0 364 1763 24 28 268 209 4514 61 0 179 0 0 0 269 311 0 231 317 9 7 63915 298700 270 84213 0 319 1 6 5 60719 315372 271 431 0 191 275 4 7 59500 318745 272 111 0 321 0 0 0 60720 315380 273 29 936 30 22 67939 408881 56 6963 274 423 8568 37 34 32168 786690 18 15856 275 202 11528 70 68 -14545 -83265 431 -4047 276 10 0 0 0 60720 315380 250 0 277 0 9 60720 321376 139 12138 164 0 278 64270 0 60720 315380 281 0 405 592 279 60951 276898 392 10985 130 -4 176 126 280 60743 315547 171 11555 921 0 236 2 281 60720 315487 172 115487 52257 0 142 0 282 8 315380 220 0 13 18014 239 393 283 256 5043 67 943 281 0 0 0 284 9 0 328 0 0 0 60720 315380 285 207 37704 452 228 60448 983660 14 1734 286 73 63 4 2 65688 318574 148 29643 287 242 1431 34 35 106885 310768 330 3258 288 229 94 3 3 61360 312887 321 37629 289 317 192 4 2 65276 312339 327 28085 290 44 32 1 0 59988 314964 301 114964 291 349 6869 105 88 117520 379983 68 1714 292 0 0 0 0 60720 315380 176 0 293 60722 0 60720 315380 188 0 361 1 294 60720 315398 173 28849 172235 0 203 0 295 400 315380 247 0 255 2328 94 57 296 157 570 23 0 308 209 13 19 297 166 8973 558 0 195 28 2 3 298 331 57844 4148 0 280 176 28 8 299 424 3952 629 0 415 1920 15 11 300 3 -2306 -18 -6 3 87550 458 412 301 383 8541 45 17936 151 520 10 7 302 231 9165 176 0 304 0 0 0 303 53 0 421 1013 2 1 69770 253468 304 7687 -14 217 15 8 6 60831 315688 305 214 0 353 587 12 10 59595 325699 306 46 0 58 5371 66 47 87720 446211 307 77 49 359 0 0 0 60720 315380 308 228 1012 47 29 114768 368078 82 3576 309 59 0 0 0 60720 315380 215 0 310 1081 33 138971 352850 96 2591 174 24 311 247 807 213118 6282154 1 5626 37 56289 312 12 190 32648 227132 414 110 1 929 313 43 9 83620 283910 388 6993 25 -3 314 0 11 74015 236761 409 855 45 26 315 60720 0 60720 315380 184 0 311 0 316 32 315380 244 0 39 2365 79 83 317 343 4438 148 120 101 1261 33 31 318 238 3258 85 13 256 0 0 0 319 98 0 123 1585 107 112 93099 355864 320 10 2 24 16189 295 290 116384 358589 321 273 379 284 0 0 0 60720 315380 322 92 0 0 0 60720 315380 265 0 323 20 92 110309 375195 76 1904 17 4 324 0 10 61977 288985 386 4449 188 0 325 90262 0 60720 315380 178 0 87 3642 326 60720 458343 42 3004 71 22 170 0 327 397 315380 192 0 200 472 20 19 328 298 3479 148 0 158 98 2 2 329 428 57618 1179 0 427 3999 30 21 330 216 -5242 -39 -58 233 0 0 0 331 174 0 345 621 48 13 60720 308256 332 173 0 375 0 0 0 60720 315380 333 4 0 0 0 60720 315380 189 0 334 0 4 60798 313164 320 28291 3808 0 335 71561 0 60720 315380 291 0 106 746 336 60720 269661 396 3029 93 10 208 0 337 255 315380 228 0 282 0 0 0 338 240 0 249 0 0 0 60720 315380 339 150 0 0 0 60720 315380 194 0 340 158 136 134759 518365 35 2122 77 66 341 18 114 156608 233773 411 214 7 58 342 15 15 39625 301881 358 5660 125 23 343 17 11 56750 298568 364 6571 98 0 344 14 12 87390 325479 128 7381 253 0 345 110 8 58990 325506 127 8965 322 0 346 152 122 48020 984885 13 7135 62 1871 347 6 113 60720 313267 319 745 283 0 348 0 6 60349 315793 163 19299 2197 0 349 67038 0 60720 315380 278 0 426 5109 350 113761 215362 417 158 3 -31 104 576 351 58320 314073 313 6337 198 12 409 438 352 62841 298096 366 8918 224 -5 292 165 353 79804 325176 130 5216 759 0 40 4069 354 62555 207393 420 51 2 119 185 23 355 99489 314806 306 22961 4974 0 59 1285 356 90131 341340 106 1321 110 49 82 4677 357 95350 426280 52 3017 48 24 15 24811 358 64245 929118 15 1444 29 822 253 167 359 60720 307322 345 53661 642 0 342 0 360 64 315380 187 0 45 4628 62 39 361 382 3024 41 85 189 226 5 5 362 404 18357 406 0 76 1765 20 20 363 122 2353 27 28 145 460 15 15 364 162 8652 282 0 144 36 3 4 365 350 38611 3215 0 408 989 11 12 366 75 9505 106 -5 96 3055 43 30 367 288 4108 58 14 265 0 0 0 368 246 0 161 0 0 0 60720 315380 369 0 0 0 0 60720 315380 236 0 370 60720 0 60720 315380 202 0 230 0 371 90 315380 214 0 66 13253 947 251 372 168 167 12 35 244 24 8 7 373 208 14455 4772 0 153 0 0 5 374 222 0 199 0 0 0 60720 315380 375 0 0 0 0 60720 315380 198 0 376 59661 0 60720 315380 289 0 387 131 377 102725 327071 123 21179 969 0 85 514 378 108479 377516 73 11834 346 23 103 3366 379 87419 299243 362 1241 29 12 60 9327 380 54683 387699 63 1858 20 44 309 384 381 122844 309836 333 4993 286 0 20 17821 382 62710 444477 46 655 14 643 391 397 383 60720 322327 136 12233 308 -1 139 897 384 60720 314913 303 22983 128 0 370 218 385 87161 315553 170 5503 531 0 31 3369 386 101481 688779 23 17456 145 217 109 5702 387 128294 321896 138 1998 21 10 41 8636 388 62620 301607 359 119 12 108 215 534 389 60720 304485 351 14926 196 0 332 0 390 227 315380 275 0 310 0 0 0 391 44 0 410 726 3 3 69980 231861 392 107 -5 81 1380 26 29 60982 347385 393 647 25 365 180 11 9 59635 316386 394 40 0 53 7285 115 97 188873 491303 395 70 59 67 880 74 15 80791 261216 396 117743 34 335 1 2 1 60727 315388 397 334 0 380 0 0 0 60720 315380 398 95 96 5 4 60379 313729 316 22746 399 8 1889 45 43 37527 358649 88 3526 400 221 45187 353 267 234817 1926517 4 4891 401 419 288 4 3 60510 296656 373 24164 402 424 1270 26 24 69206 275311 393 2897 403 327 6526 26 20 55830 -42143 430 -9313 404 12 0 0 0 60720 315380 252 0 405 8 13 72835 343929 102 11994 636 0 406 0 6 68060 367655 83 20957 241 -2 407 56726 0 60720 315380 241 0 283 249 408 60720 313491 317 1576 455 0 350 1595 409 94355 319369 143 7022 75 0 61 5565 410 60720 493408 38 6243 53 39 287 601 411 77655 319210 144 19868 198 0 268 188 412 134028 381180 67 16471 964 0 25 7146 413 62285 297978 367 933 14 305 418 1135 414 59325 290476 384 5322 80 -9 392 450 415 60630 292136 380 11517 205 -1 131 34 416 65990 314353 309 14294 3340 0 316 133 417 59118 339445 109 9960 1052 0 347 119 418 100423 303677 354 17280 872 0 68 2053 419 100269 397144 60 3720 96 33 70 4036 420 60720 424898 53 3570 56 32 147 0 421 105 315380 203 0 169 4655 393 55 422 338 360 30 0 248 131 11 15 423 349 9908 834 0 152 1766 73 46 424 147 1451 60 0 351 312 14 3 425 132 8478 381 0 372 448 40 7 426 145 3084 275 0 137 115 13 10 427 314 9146 1037 0 352 60 10 7 428 180 11405 1916 0 338 0 0 0 429 271 0 258 364 35 26 47705 298466 430 398 0 343 0 0 0 60720 315380 431 115 1442 118 239 121173 438493 49 2021 Profit/Trades Profit/Cost 1 7022 75 2 6243 53 3 19868 198 4 16471 964 5 933 14 6 5322 80 7 11517 205 8 14294 3340 9 9960 1052 10 17280 872 11 3720 96 12 3570 56 13 0 169 14 0 248 15 0 152 16 0 351 17 0 372 18 0 137 19 0 352 20 0 338 21 364 35 22 0 0 23 239 121173 24 6 57530 25 6 62041 26 10 60720 27 149 136996 28 0 60720 29 332743 115 30 286963 387 31 315380 219 32 8797 430 33 2019 288 34 1889 22 35 1298 18 36 19165 12458 37 3030 32 38 2493 29 39 2304 42 40 28639 844 41 9616 242 42 1861 123 43 0 184 44 0 401 45 -4 371 46 203 9 47 3689 25 48 567 11 49 27664 282 50 1686 40 51 1164 17 52 15824 64 53 6575 159 54 37597 248 55 314 5 56 1932 36 57 5692 72 58 0 0 59 0 60720 60 364839 85 61 230621 413 62 315877 161 63 125390 427 64 314882 305 65 370837 79 66 330068 120 67 324385 131 68 315380 277 69 11409 76 70 1115 6 71 13337 138 72 0 187 73 0 213 74 0 296 75 0 218 76 0 28 77 238 5 78 3040 138 79 0 163 80 0 0 81 17 62920 82 2 60720 83 315380 263 84 16362 1212 85 57677 36505 86 2004 225 87 15182 603 88 0 251 89 0 11 90 1910 62 91 38 149 92 0 270 93 0 171 94 278 26 95 282 30 96 571 80 97 2253 99 98 290 2 99 78 16 100 20 2 101 0 0 102 69 27284 103 28 56225 104 168 157214 105 13 54323 106 23 57935 107 54 59017 108 2 73490 109 0 60720 110 317046 156 111 315380 240 112 13695 300 113 6308 69 114 5523 51 115 6141 93 116 3062 105 117 953 15 118 0 33 119 195 91 120 19 430 121 -158 259 122 530 60 123 0 0 124 244 157278 125 9 73221 126 12 67000 127 0 60720 128 315380 254 129 0 389 130 0 210 131 194 0 132 33 23 133 0 0 134 60720 315380 135 37 1335 136 422 35 137 196 0 138 186 638 139 262 0 140 11 7 141 63 25 142 117 29 143 0 0 144 60720 315380 145 371 12115 146 225 0 147 118 0 148 70 -12 149 37 6 150 394 0 151 354 135 152 124 1659 153 136 316 154 130 0 155 11 9 156 8 14 157 0 0 158 102129 597793 159 60720 315380 160 165 1677 161 69 3643 162 326 960 163 59 9253 164 36 15542 165 418 888 166 410 1528 167 103 4488 168 84 7903 169 210 0 170 1874 0 171 281 0 172 49 4901 173 88 540 174 329 0 175 0 0 176 60720 315380 177 174 12822 178 325 28103 179 113 14834 180 415 39 181 217 0 182 37920 0 183 182 -3 184 159 405 185 181 972 186 278 0 187 64 45 188 47 33 189 0 0 190 60720 315380 191 20 6095 192 258 0 193 164 483 194 358 0 195 21 17 196 41 40 197 0 9 198 87192 346611 199 55792 314289 200 114337 856956 201 75832 217193 202 61630 315366 203 58580 307930 204 165548 702380 205 76403 194493 206 61656 316155 207 60720 315380 208 119 659 209 61 5557 210 323 37615 211 190 0 212 114 0 213 314 0 214 213 0 215 10 143 216 17 7123 217 254 622 218 374 174 219 35 2220 220 377 121 221 32 11819 222 260 0 223 12 19 224 13 11 225 39 35 226 0 0 227 60720 315380 228 149 29525 229 44 7523 230 125 11431 231 234 0 232 239 0 233 9 -3 234 64 -8 235 340 125 236 239 0 237 2 1 238 0 0 239 78780 322378 240 60720 315380 241 270 0 242 333 0 243 0 0 244 60720 315380 245 27 4838 246 100 2390 247 25 50325 248 48 664 249 391 3160 250 108 818 251 407 610 252 229 0 253 264 0 254 285 333 255 339 26 256 243 0 257 0 0 258 69863 309422 259 60938 315245 260 61348 405972 261 50804 300962 262 60745 316176 263 59506 302409 264 58456 283587 265 60950 263276 266 60720 312075 267 61600 308336 268 60720 315380 269 363 10967 270 293 19229 271 146 29686 272 246 0 273 223 9 274 68 223 275 -25 -25 276 232 738 277 266 0 278 7 7 279 10 8 280 1 0 281 0 0 282 -1710 1405225 283 60720 315380 284 249 0 285 21 2167 286 1891 0 287 77 30 288 1200 0 289 586 0 290 3603 0 291 26 90 292 140 0 293 4 4 294 0 0 295 82732 253588 296 64016 316647 297 60890 315688 298 68136 310670 299 79420 165404 300 153198 4111912 301 58650 291650 302 60720 315380 303 402 26734 304 164 14461 305 126 10475 306 45 3730 307 251 0 308 166 26 309 83 876 310 1 162556 311 14 43556 312 395 3425 313 78 810 314 190 0 315 0 0 316 191778 550608 317 76114 307528 318 60720 315380 319 92 1457 320 89 538 321 242 0 322 121 10579 323 214 474 324 348 0 325 86 100 326 0 0 327 80045 269587 328 61490 315236 329 51252 42754 330 60720 315380 331 341 2255 332 266 0 333 301 30 334 368 0 335 23 29 336 0 0 337 60720 315380 338 260 0 339 52 4150 340 55 4658 341 84 814 342 143 1002 343 205 496 344 141 389 345 12 12679 346 177 400 347 238 53 348 331 0 349 97 82 350 18 8 351 11 9 352 24 12 353 145 131 354 5 5 355 107 51 356 75 66 357 505 280 358 2 1 359 0 0 360 69159 387475 361 65745 291787 362 77623 247060 363 63346 329784 364 60894 315834 365 58930 304555 366 60247 376641 367 60720 315380 368 212 0 369 237 0 370 0 0 371 90829 357760 372 59818 315637 373 60720 315380 374 261 0 375 363 0 376 6 3 377 15 15 378 80 55 379 101 85 380 22 17 381 373 277 382 10 10 383 5 4 384 21 7 385 28 26 386 61 47 387 855 203 388 7 5 389 0 0 390 60720 315380 391 412 10620 392 98 5669 393 158 10581 394 39 2533 395 399 827 396 175 57694 397 271 0 398 1188 0 399 84 14 400 38 2854 401 335 0 402 59 -10 403 -37 -28 404 263 226 405 394 694 406 286 0 407 72 26 408 17 10 409 47 43 410 6 5 411 11 9 412 105 83 413 17 20 414 8 6 415 8 7 416 14 4 417 6 3 418 53 44 419 63 85 420 0 0 421 60720 341570 422 61808 308989 423 60438 305959 424 58598 318690 425 59781 323361 426 60945 318903 427 61124 314049 428 60720 315380 429 365 2813 430 183 0 431 165 -3 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Rank Costs Trades 6.141e+04 4.883e-03 -4.834e-02 -2.054e-02 Orders Dividends Wealth WRank -1.038e-01 -1.064e-01 -3.090e-02 -1.430e-01 `Profit/Trades` `Profit/Cost` -2.286e-01 -1.838e-02 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -72197.9 -41308.8 -21438.5 -954.9 1572683.1 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 6.141e+04 9.121e+03 6.733 5.48e-11 *** Rank 4.883e-03 4.073e-02 0.120 0.904629 Costs -4.834e-02 4.742e-02 -1.019 0.308658 Trades -2.054e-02 1.566e-02 -1.312 0.190362 Orders -1.038e-01 4.262e-02 -2.436 0.015263 * Dividends -1.064e-01 3.580e-02 -2.972 0.003126 ** Wealth -3.090e-02 2.017e-02 -1.532 0.126316 WRank -1.430e-01 4.258e-02 -3.358 0.000856 *** `Profit/Trades` -2.286e-01 8.420e-02 -2.715 0.006909 ** `Profit/Cost` -1.838e-02 2.160e-02 -0.851 0.395413 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 102900 on 421 degrees of freedom Multiple R-squared: 0.06598, Adjusted R-squared: 0.04602 F-statistic: 3.305 on 9 and 421 DF, p-value: 0.0006545 > 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,] 1.928095e-07 3.856190e-07 9.999998e-01 [2,] 9.365485e-10 1.873097e-09 1.000000e+00 [3,] 1.664539e-10 3.329077e-10 1.000000e+00 [4,] 1.299948e-12 2.599896e-12 1.000000e+00 [5,] 9.386946e-14 1.877389e-13 1.000000e+00 [6,] 1.376855e-15 2.753710e-15 1.000000e+00 [7,] 1.815196e-17 3.630391e-17 1.000000e+00 [8,] 1.725222e-19 3.450444e-19 1.000000e+00 [9,] 2.249187e-21 4.498374e-21 1.000000e+00 [10,] 7.212513e-23 1.442503e-22 1.000000e+00 [11,] 3.709235e-19 7.418470e-19 1.000000e+00 [12,] 1.591687e-15 3.183375e-15 1.000000e+00 [13,] 6.228989e-09 1.245798e-08 1.000000e+00 [14,] 1.974784e-09 3.949567e-09 1.000000e+00 [15,] 3.709589e-06 7.419179e-06 9.999963e-01 [16,] 1.447513e-06 2.895025e-06 9.999986e-01 [17,] 1.435949e-06 2.871897e-06 9.999986e-01 [18,] 5.587233e-07 1.117447e-06 9.999994e-01 [19,] 2.386952e-07 4.773903e-07 9.999998e-01 [20,] 8.391028e-08 1.678206e-07 9.999999e-01 [21,] 9.599232e-08 1.919846e-07 9.999999e-01 [22,] 3.203073e-08 6.406146e-08 1.000000e+00 [23,] 1.430379e-08 2.860759e-08 1.000000e+00 [24,] 1.025203e-08 2.050405e-08 1.000000e+00 [25,] 4.944960e-09 9.889921e-09 1.000000e+00 [26,] 3.621164e-08 7.242329e-08 1.000000e+00 [27,] 1.385140e-08 2.770281e-08 1.000000e+00 [28,] 6.142344e-09 1.228469e-08 1.000000e+00 [29,] 2.418477e-09 4.836953e-09 1.000000e+00 [30,] 9.149042e-10 1.829808e-09 1.000000e+00 [31,] 3.456609e-10 6.913218e-10 1.000000e+00 [32,] 1.198532e-10 2.397064e-10 1.000000e+00 [33,] 4.045375e-11 8.090750e-11 1.000000e+00 [34,] 1.344809e-11 2.689618e-11 1.000000e+00 [35,] 5.662901e-12 1.132580e-11 1.000000e+00 [36,] 3.873235e-12 7.746470e-12 1.000000e+00 [37,] 1.581795e-12 3.163590e-12 1.000000e+00 [38,] 5.312941e-13 1.062588e-12 1.000000e+00 [39,] 1.799766e-13 3.599532e-13 1.000000e+00 [40,] 5.991886e-14 1.198377e-13 1.000000e+00 [41,] 1.886099e-14 3.772199e-14 1.000000e+00 [42,] 1.509394e-14 3.018788e-14 1.000000e+00 [43,] 5.649500e-15 1.129900e-14 1.000000e+00 [44,] 1.904541e-15 3.809082e-15 1.000000e+00 [45,] 6.278545e-16 1.255709e-15 1.000000e+00 [46,] 2.085260e-16 4.170519e-16 1.000000e+00 [47,] 1.196368e-15 2.392736e-15 1.000000e+00 [48,] 4.088728e-16 8.177455e-16 1.000000e+00 [49,] 1.290948e-16 2.581895e-16 1.000000e+00 [50,] 3.586156e-15 7.172312e-15 1.000000e+00 [51,] 1.447168e-15 2.894335e-15 1.000000e+00 [52,] 5.371543e-16 1.074309e-15 1.000000e+00 [53,] 1.994807e-16 3.989614e-16 1.000000e+00 [54,] 7.110466e-17 1.422093e-16 1.000000e+00 [55,] 2.606646e-17 5.213292e-17 1.000000e+00 [56,] 1.115604e-17 2.231209e-17 1.000000e+00 [57,] 3.651416e-18 7.302832e-18 1.000000e+00 [58,] 1.205666e-18 2.411332e-18 1.000000e+00 [59,] 5.418114e-19 1.083623e-18 1.000000e+00 [60,] 1.781818e-19 3.563636e-19 1.000000e+00 [61,] 5.725905e-20 1.145181e-19 1.000000e+00 [62,] 1.813576e-20 3.627152e-20 1.000000e+00 [63,] 5.706761e-21 1.141352e-20 1.000000e+00 [64,] 1.876388e-21 3.752776e-21 1.000000e+00 [65,] 1.215467e-21 2.430934e-21 1.000000e+00 [66,] 4.021006e-22 8.042013e-22 1.000000e+00 [67,] 1.226516e-22 2.453032e-22 1.000000e+00 [68,] 3.731371e-23 7.462742e-23 1.000000e+00 [69,] 1.801338e-21 3.602676e-21 1.000000e+00 [70,] 2.327290e-21 4.654581e-21 1.000000e+00 [71,] 9.791192e-22 1.958238e-21 1.000000e+00 [72,] 3.473498e-22 6.946995e-22 1.000000e+00 [73,] 2.733197e-19 5.466394e-19 1.000000e+00 [74,] 1.016691e-19 2.033382e-19 1.000000e+00 [75,] 3.666402e-20 7.332805e-20 1.000000e+00 [76,] 1.354764e-20 2.709528e-20 1.000000e+00 [77,] 5.322724e-21 1.064545e-20 1.000000e+00 [78,] 5.555997e-21 1.111199e-20 1.000000e+00 [79,] 2.069717e-21 4.139434e-21 1.000000e+00 [80,] 7.545624e-22 1.509125e-21 1.000000e+00 [81,] 2.688038e-22 5.376077e-22 1.000000e+00 [82,] 9.575049e-23 1.915010e-22 1.000000e+00 [83,] 3.432656e-23 6.865313e-23 1.000000e+00 [84,] 1.218114e-23 2.436229e-23 1.000000e+00 [85,] 4.399510e-24 8.799020e-24 1.000000e+00 [86,] 1.515806e-24 3.031611e-24 1.000000e+00 [87,] 5.243559e-25 1.048712e-24 1.000000e+00 [88,] 1.788490e-25 3.576981e-25 1.000000e+00 [89,] 6.004071e-26 1.200814e-25 1.000000e+00 [90,] 3.207239e-24 6.414479e-24 1.000000e+00 [91,] 1.690948e-23 3.381897e-23 1.000000e+00 [92,] 1.000000e+00 4.118595e-14 2.059297e-14 [93,] 1.000000e+00 1.141019e-15 5.705096e-16 [94,] 1.000000e+00 5.133154e-18 2.566577e-18 [95,] 1.000000e+00 5.328175e-22 2.664087e-22 [96,] 1.000000e+00 5.269806e-25 2.634903e-25 [97,] 1.000000e+00 1.219346e-29 6.096732e-30 [98,] 1.000000e+00 1.244700e-29 6.223501e-30 [99,] 1.000000e+00 8.576256e-30 4.288128e-30 [100,] 1.000000e+00 1.955080e-29 9.775402e-30 [101,] 1.000000e+00 4.438954e-29 2.219477e-29 [102,] 1.000000e+00 9.946944e-29 4.973472e-29 [103,] 1.000000e+00 2.201168e-28 1.100584e-28 [104,] 1.000000e+00 4.901482e-28 2.450741e-28 [105,] 1.000000e+00 1.067803e-27 5.339016e-28 [106,] 1.000000e+00 2.337824e-27 1.168912e-27 [107,] 1.000000e+00 4.503215e-27 2.251608e-27 [108,] 1.000000e+00 8.932926e-27 4.466463e-27 [109,] 1.000000e+00 1.914793e-26 9.573964e-27 [110,] 1.000000e+00 4.127931e-26 2.063965e-26 [111,] 1.000000e+00 4.237611e-26 2.118806e-26 [112,] 1.000000e+00 2.144143e-41 1.072071e-41 [113,] 1.000000e+00 4.837381e-50 2.418691e-50 [114,] 1.000000e+00 2.226221e-64 1.113110e-64 [115,] 1.000000e+00 4.223280e-85 2.111640e-85 [116,] 1.000000e+00 2.067975e-88 1.033987e-88 [117,] 1.000000e+00 8.762669e-88 4.381335e-88 [118,] 1.000000e+00 3.017415e-87 1.508707e-87 [119,] 1.000000e+00 8.515448e-87 4.257724e-87 [120,] 1.000000e+00 6.090979e-87 3.045490e-87 [121,] 1.000000e+00 6.381082e-87 3.190541e-87 [122,] 1.000000e+00 4.733147e-115 2.366574e-115 [123,] 1.000000e+00 1.712596e-114 8.562982e-115 [124,] 1.000000e+00 5.611636e-114 2.805818e-114 [125,] 1.000000e+00 2.744593e-113 1.372297e-113 [126,] 1.000000e+00 7.594713e-114 3.797357e-114 [127,] 1.000000e+00 2.732410e-113 1.366205e-113 [128,] 1.000000e+00 3.694911e-113 1.847455e-113 [129,] 1.000000e+00 2.078773e-116 1.039386e-116 [130,] 1.000000e+00 6.454111e-118 3.227055e-118 [131,] 1.000000e+00 8.152367e-118 4.076183e-118 [132,] 1.000000e+00 1.763691e-124 8.818457e-125 [133,] 1.000000e+00 8.404352e-124 4.202176e-124 [134,] 1.000000e+00 4.726755e-123 2.363377e-123 [135,] 1.000000e+00 5.433755e-123 2.716877e-123 [136,] 1.000000e+00 2.011110e-122 1.005555e-122 [137,] 1.000000e+00 6.705104e-122 3.352552e-122 [138,] 1.000000e+00 2.023952e-121 1.011976e-121 [139,] 1.000000e+00 8.429146e-121 4.214573e-121 [140,] 1.000000e+00 3.470587e-120 1.735294e-120 [141,] 1.000000e+00 1.561740e-119 7.808702e-120 [142,] 1.000000e+00 5.864891e-119 2.932446e-119 [143,] 1.000000e+00 1.217800e-118 6.088998e-119 [144,] 1.000000e+00 1.168757e-118 5.843783e-119 [145,] 1.000000e+00 1.798572e-118 8.992861e-119 [146,] 1.000000e+00 7.179302e-126 3.589651e-126 [147,] 1.000000e+00 2.034634e-127 1.017317e-127 [148,] 1.000000e+00 1.203975e-126 6.019874e-127 [149,] 1.000000e+00 7.233841e-126 3.616920e-126 [150,] 1.000000e+00 4.327265e-125 2.163632e-125 [151,] 1.000000e+00 2.710695e-124 1.355347e-124 [152,] 1.000000e+00 1.627255e-123 8.136274e-124 [153,] 1.000000e+00 8.249502e-123 4.124751e-123 [154,] 1.000000e+00 4.477126e-122 2.238563e-122 [155,] 1.000000e+00 2.594543e-121 1.297272e-121 [156,] 1.000000e+00 1.593297e-120 7.966486e-121 [157,] 1.000000e+00 9.143051e-120 4.571526e-120 [158,] 1.000000e+00 4.048928e-119 2.024464e-119 [159,] 1.000000e+00 1.938468e-118 9.692341e-119 [160,] 1.000000e+00 9.421663e-118 4.710831e-118 [161,] 1.000000e+00 4.540633e-117 2.270317e-117 [162,] 1.000000e+00 2.079542e-116 1.039771e-116 [163,] 1.000000e+00 3.679945e-116 1.839972e-116 [164,] 1.000000e+00 4.100685e-116 2.050343e-116 [165,] 1.000000e+00 1.867926e-116 9.339631e-117 [166,] 1.000000e+00 1.061320e-115 5.306599e-116 [167,] 1.000000e+00 6.086208e-115 3.043104e-115 [168,] 1.000000e+00 2.150213e-114 1.075107e-114 [169,] 1.000000e+00 1.206555e-113 6.032776e-114 [170,] 1.000000e+00 5.234621e-113 2.617311e-113 [171,] 1.000000e+00 2.687372e-112 1.343686e-112 [172,] 1.000000e+00 1.345433e-111 6.727163e-112 [173,] 1.000000e+00 5.195519e-111 2.597759e-111 [174,] 1.000000e+00 2.038610e-110 1.019305e-110 [175,] 1.000000e+00 2.234770e-112 1.117385e-112 [176,] 1.000000e+00 8.484201e-112 4.242101e-112 [177,] 1.000000e+00 4.097925e-111 2.048962e-111 [178,] 1.000000e+00 5.018353e-111 2.509177e-111 [179,] 1.000000e+00 1.835537e-110 9.177683e-111 [180,] 1.000000e+00 1.015002e-109 5.075010e-110 [181,] 1.000000e+00 5.155973e-109 2.577986e-109 [182,] 1.000000e+00 1.925529e-108 9.627644e-109 [183,] 1.000000e+00 1.920201e-109 9.601003e-110 [184,] 1.000000e+00 3.382826e-109 1.691413e-109 [185,] 1.000000e+00 1.305704e-108 6.528518e-109 [186,] 1.000000e+00 2.763124e-108 1.381562e-108 [187,] 1.000000e+00 1.978521e-108 9.892604e-109 [188,] 1.000000e+00 2.251545e-109 1.125773e-109 [189,] 1.000000e+00 6.808941e-109 3.404470e-109 [190,] 1.000000e+00 1.896998e-108 9.484988e-109 [191,] 1.000000e+00 6.105577e-108 3.052789e-108 [192,] 1.000000e+00 1.585290e-108 7.926450e-109 [193,] 1.000000e+00 5.798807e-108 2.899404e-108 [194,] 1.000000e+00 2.171538e-107 1.085769e-107 [195,] 1.000000e+00 9.535009e-107 4.767505e-107 [196,] 1.000000e+00 5.523995e-106 2.761997e-106 [197,] 1.000000e+00 3.321955e-105 1.660978e-105 [198,] 1.000000e+00 1.900115e-104 9.500573e-105 [199,] 1.000000e+00 1.054013e-103 5.270064e-104 [200,] 1.000000e+00 5.536759e-103 2.768380e-103 [201,] 1.000000e+00 2.917389e-102 1.458695e-102 [202,] 1.000000e+00 1.503081e-101 7.515404e-102 [203,] 1.000000e+00 7.452493e-101 3.726246e-101 [204,] 1.000000e+00 3.830105e-100 1.915053e-100 [205,] 1.000000e+00 1.909399e-99 9.546994e-100 [206,] 1.000000e+00 7.444187e-99 3.722093e-99 [207,] 1.000000e+00 3.123916e-98 1.561958e-98 [208,] 1.000000e+00 1.505572e-97 7.527859e-98 [209,] 1.000000e+00 5.252815e-97 2.626408e-97 [210,] 1.000000e+00 2.339569e-96 1.169785e-96 [211,] 1.000000e+00 3.849644e-96 1.924822e-96 [212,] 1.000000e+00 2.056373e-95 1.028186e-95 [213,] 1.000000e+00 1.105656e-94 5.528280e-95 [214,] 1.000000e+00 6.061924e-94 3.030962e-94 [215,] 1.000000e+00 1.443164e-93 7.215821e-94 [216,] 1.000000e+00 7.941862e-93 3.970931e-93 [217,] 1.000000e+00 4.446944e-92 2.223472e-92 [218,] 1.000000e+00 2.411258e-91 1.205629e-91 [219,] 1.000000e+00 1.244420e-90 6.222099e-91 [220,] 1.000000e+00 6.117826e-90 3.058913e-90 [221,] 1.000000e+00 2.806698e-89 1.403349e-89 [222,] 1.000000e+00 1.162715e-88 5.813576e-89 [223,] 1.000000e+00 5.561689e-88 2.780845e-88 [224,] 1.000000e+00 1.943976e-87 9.719878e-88 [225,] 1.000000e+00 3.176378e-87 1.588189e-87 [226,] 1.000000e+00 1.707795e-86 8.538973e-87 [227,] 1.000000e+00 6.147069e-86 3.073534e-86 [228,] 1.000000e+00 2.456841e-85 1.228420e-85 [229,] 1.000000e+00 1.231484e-84 6.157419e-85 [230,] 1.000000e+00 5.792806e-84 2.896403e-84 [231,] 1.000000e+00 8.932073e-84 4.466036e-84 [232,] 1.000000e+00 1.959298e-83 9.796492e-84 [233,] 1.000000e+00 9.060522e-83 4.530261e-83 [234,] 1.000000e+00 4.827779e-82 2.413889e-82 [235,] 1.000000e+00 2.338161e-81 1.169081e-81 [236,] 1.000000e+00 1.200790e-80 6.003949e-81 [237,] 1.000000e+00 5.750673e-80 2.875336e-80 [238,] 1.000000e+00 3.017867e-79 1.508934e-79 [239,] 1.000000e+00 1.386012e-78 6.930060e-79 [240,] 1.000000e+00 6.622441e-78 3.311220e-78 [241,] 1.000000e+00 3.033368e-77 1.516684e-77 [242,] 1.000000e+00 1.369174e-76 6.845871e-77 [243,] 1.000000e+00 5.100099e-76 2.550050e-76 [244,] 1.000000e+00 1.982996e-75 9.914980e-76 [245,] 1.000000e+00 3.104218e-75 1.552109e-75 [246,] 1.000000e+00 8.365977e-75 4.182989e-75 [247,] 1.000000e+00 3.371100e-74 1.685550e-74 [248,] 1.000000e+00 1.423013e-73 7.115064e-74 [249,] 1.000000e+00 5.762128e-73 2.881064e-73 [250,] 1.000000e+00 2.424375e-72 1.212188e-72 [251,] 1.000000e+00 9.977631e-72 4.988815e-72 [252,] 1.000000e+00 4.258629e-71 2.129315e-71 [253,] 1.000000e+00 1.828667e-70 9.143333e-71 [254,] 1.000000e+00 7.930111e-70 3.965055e-70 [255,] 1.000000e+00 3.491677e-69 1.745838e-69 [256,] 1.000000e+00 1.521400e-68 7.606999e-69 [257,] 1.000000e+00 6.863856e-68 3.431928e-68 [258,] 1.000000e+00 4.507764e-69 2.253882e-69 [259,] 1.000000e+00 2.165041e-68 1.082521e-68 [260,] 1.000000e+00 1.027159e-67 5.135796e-68 [261,] 1.000000e+00 4.494335e-67 2.247167e-67 [262,] 1.000000e+00 1.173564e-66 5.867820e-67 [263,] 1.000000e+00 1.543240e-66 7.716200e-67 [264,] 1.000000e+00 6.597828e-66 3.298914e-66 [265,] 1.000000e+00 2.237721e-65 1.118860e-65 [266,] 1.000000e+00 2.787890e-65 1.393945e-65 [267,] 1.000000e+00 1.392163e-64 6.960814e-65 [268,] 1.000000e+00 6.787939e-64 3.393970e-64 [269,] 1.000000e+00 3.184681e-63 1.592341e-63 [270,] 1.000000e+00 9.625090e-65 4.812545e-65 [271,] 1.000000e+00 4.428825e-64 2.214412e-64 [272,] 1.000000e+00 2.105547e-63 1.052773e-63 [273,] 1.000000e+00 5.483550e-63 2.741775e-63 [274,] 1.000000e+00 2.491187e-62 1.245593e-62 [275,] 1.000000e+00 1.226359e-61 6.131793e-62 [276,] 1.000000e+00 5.379161e-61 2.689581e-61 [277,] 1.000000e+00 2.330648e-60 1.165324e-60 [278,] 1.000000e+00 1.091801e-59 5.459003e-60 [279,] 1.000000e+00 5.294156e-59 2.647078e-59 [280,] 1.000000e+00 2.052181e-58 1.026090e-58 [281,] 1.000000e+00 3.097690e-58 1.548845e-58 [282,] 1.000000e+00 1.456375e-58 7.281875e-59 [283,] 1.000000e+00 5.889353e-58 2.944677e-58 [284,] 1.000000e+00 2.688509e-57 1.344255e-57 [285,] 1.000000e+00 1.223109e-56 6.115543e-57 [286,] 1.000000e+00 5.675671e-56 2.837835e-56 [287,] 1.000000e+00 2.246913e-55 1.123456e-55 [288,] 1.000000e+00 3.123522e-55 1.561761e-55 [289,] 1.000000e+00 1.417194e-54 7.085968e-55 [290,] 1.000000e+00 6.366689e-54 3.183344e-54 [291,] 1.000000e+00 2.716247e-53 1.358124e-53 [292,] 1.000000e+00 1.297598e-52 6.487988e-53 [293,] 1.000000e+00 5.811145e-52 2.905573e-52 [294,] 1.000000e+00 2.740939e-51 1.370470e-51 [295,] 1.000000e+00 1.179647e-50 5.898234e-51 [296,] 1.000000e+00 3.958531e-50 1.979266e-50 [297,] 1.000000e+00 1.601906e-49 8.009528e-50 [298,] 1.000000e+00 6.344121e-49 3.172060e-49 [299,] 1.000000e+00 6.266695e-49 3.133347e-49 [300,] 1.000000e+00 1.299599e-48 6.497993e-49 [301,] 1.000000e+00 4.472631e-48 2.236315e-48 [302,] 1.000000e+00 1.337887e-47 6.689434e-48 [303,] 1.000000e+00 2.371788e-47 1.185894e-47 [304,] 1.000000e+00 5.646440e-48 2.823220e-48 [305,] 1.000000e+00 2.093141e-47 1.046571e-47 [306,] 1.000000e+00 9.550081e-47 4.775040e-47 [307,] 1.000000e+00 4.478373e-46 2.239186e-46 [308,] 1.000000e+00 2.103534e-45 1.051767e-45 [309,] 1.000000e+00 9.004246e-45 4.502123e-45 [310,] 1.000000e+00 3.232914e-44 1.616457e-44 [311,] 1.000000e+00 1.351181e-43 6.755907e-44 [312,] 1.000000e+00 3.598966e-43 1.799483e-43 [313,] 1.000000e+00 6.313216e-44 3.156608e-44 [314,] 1.000000e+00 2.035552e-43 1.017776e-43 [315,] 1.000000e+00 9.503430e-43 4.751715e-43 [316,] 1.000000e+00 4.287300e-42 2.143650e-42 [317,] 1.000000e+00 1.913763e-41 9.568813e-42 [318,] 1.000000e+00 8.535603e-41 4.267802e-41 [319,] 1.000000e+00 3.504930e-40 1.752465e-40 [320,] 1.000000e+00 1.423280e-39 7.116400e-40 [321,] 1.000000e+00 4.966261e-39 2.483131e-39 [322,] 1.000000e+00 1.347292e-38 6.736461e-39 [323,] 1.000000e+00 1.120117e-38 5.600586e-39 [324,] 1.000000e+00 4.767228e-38 2.383614e-38 [325,] 1.000000e+00 1.660928e-37 8.304640e-38 [326,] 1.000000e+00 6.560050e-37 3.280025e-37 [327,] 1.000000e+00 2.214691e-36 1.107346e-36 [328,] 1.000000e+00 9.268992e-36 4.634496e-36 [329,] 1.000000e+00 3.902956e-35 1.951478e-35 [330,] 1.000000e+00 8.807706e-35 4.403853e-35 [331,] 1.000000e+00 2.379438e-34 1.189719e-34 [332,] 1.000000e+00 8.227117e-34 4.113559e-34 [333,] 1.000000e+00 1.932856e-33 9.664281e-34 [334,] 1.000000e+00 4.731850e-33 2.365925e-33 [335,] 1.000000e+00 9.678011e-33 4.839006e-33 [336,] 1.000000e+00 1.326692e-32 6.633461e-33 [337,] 1.000000e+00 3.184672e-32 1.592336e-32 [338,] 1.000000e+00 2.200237e-32 1.100118e-32 [339,] 1.000000e+00 8.761104e-32 4.380552e-32 [340,] 1.000000e+00 3.732605e-31 1.866303e-31 [341,] 1.000000e+00 1.599937e-30 7.999683e-31 [342,] 1.000000e+00 7.176616e-30 3.588308e-30 [343,] 1.000000e+00 1.947276e-29 9.736378e-30 [344,] 1.000000e+00 6.591495e-29 3.295748e-29 [345,] 1.000000e+00 2.066479e-28 1.033240e-28 [346,] 1.000000e+00 1.279981e-28 6.399903e-29 [347,] 1.000000e+00 2.899758e-28 1.449879e-28 [348,] 1.000000e+00 8.244914e-28 4.122457e-28 [349,] 1.000000e+00 3.485153e-27 1.742576e-27 [350,] 1.000000e+00 5.000723e-27 2.500361e-27 [351,] 1.000000e+00 2.200315e-26 1.100157e-26 [352,] 1.000000e+00 9.636146e-26 4.818073e-26 [353,] 1.000000e+00 4.020653e-25 2.010327e-25 [354,] 1.000000e+00 1.319237e-24 6.596187e-25 [355,] 1.000000e+00 5.685024e-24 2.842512e-24 [356,] 1.000000e+00 2.231921e-23 1.115960e-23 [357,] 1.000000e+00 8.556084e-23 4.278042e-23 [358,] 1.000000e+00 2.212745e-22 1.106372e-22 [359,] 1.000000e+00 6.458323e-23 3.229161e-23 [360,] 1.000000e+00 2.842761e-22 1.421380e-22 [361,] 1.000000e+00 8.803432e-22 4.401716e-22 [362,] 1.000000e+00 3.822409e-21 1.911205e-21 [363,] 1.000000e+00 1.362907e-20 6.814536e-21 [364,] 1.000000e+00 2.975640e-20 1.487820e-20 [365,] 1.000000e+00 9.238124e-20 4.619062e-20 [366,] 1.000000e+00 2.135273e-19 1.067637e-19 [367,] 1.000000e+00 6.634268e-19 3.317134e-19 [368,] 1.000000e+00 2.069009e-18 1.034504e-18 [369,] 1.000000e+00 7.024797e-19 3.512398e-19 [370,] 1.000000e+00 2.594264e-18 1.297132e-18 [371,] 1.000000e+00 1.012358e-17 5.061790e-18 [372,] 1.000000e+00 3.405298e-17 1.702649e-17 [373,] 1.000000e+00 1.204029e-16 6.020147e-17 [374,] 1.000000e+00 4.963718e-16 2.481859e-16 [375,] 1.000000e+00 4.805362e-17 2.402681e-17 [376,] 1.000000e+00 2.282096e-16 1.141048e-16 [377,] 1.000000e+00 9.406411e-16 4.703205e-16 [378,] 1.000000e+00 3.227348e-15 1.613674e-15 [379,] 1.000000e+00 9.914295e-15 4.957147e-15 [380,] 1.000000e+00 4.112555e-14 2.056277e-14 [381,] 1.000000e+00 1.327781e-13 6.638905e-14 [382,] 1.000000e+00 5.544948e-13 2.772474e-13 [383,] 1.000000e+00 2.267639e-12 1.133819e-12 [384,] 1.000000e+00 3.870710e-14 1.935355e-14 [385,] 1.000000e+00 2.068878e-13 1.034439e-13 [386,] 1.000000e+00 8.408664e-13 4.204332e-13 [387,] 1.000000e+00 2.372578e-12 1.186289e-12 [388,] 1.000000e+00 7.839312e-12 3.919656e-12 [389,] 1.000000e+00 2.991604e-11 1.495802e-11 [390,] 1.000000e+00 1.401894e-10 7.009470e-11 [391,] 1.000000e+00 1.748043e-10 8.740215e-11 [392,] 1.000000e+00 3.298829e-10 1.649415e-10 [393,] 1.000000e+00 1.726861e-09 8.634303e-10 [394,] 1.000000e+00 2.218296e-09 1.109148e-09 [395,] 1.000000e+00 1.104786e-08 5.523929e-09 [396,] 1.000000e+00 6.265773e-08 3.132886e-08 [397,] 9.999999e-01 1.949240e-07 9.746202e-08 [398,] 9.999997e-01 6.579573e-07 3.289787e-07 [399,] 9.999991e-01 1.753274e-06 8.766368e-07 [400,] 9.999996e-01 7.630417e-07 3.815209e-07 [401,] 9.999972e-01 5.526710e-06 2.763355e-06 [402,] 9.999842e-01 3.161723e-05 1.580862e-05 [403,] 9.999634e-01 7.320782e-05 3.660391e-05 [404,] 9.997875e-01 4.250315e-04 2.125158e-04 [405,] 9.992909e-01 1.418267e-03 7.091334e-04 [406,] 9.997988e-01 4.024794e-04 2.012397e-04 > postscript(file="/var/www/html/rcomp/tmp/18klz1291380013.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/html/rcomp/tmp/28klz1291380013.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/html/rcomp/tmp/3ic221291380013.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/html/rcomp/tmp/4ic221291380013.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/html/rcomp/tmp/5ic221291380013.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 = 431 Frequency = 1 1 2 3 4 5 6 -43371.7534 -34371.7678 -40486.6638 -37562.3114 -37007.6557 -44482.4940 7 8 9 10 11 12 -43355.1770 -41865.3726 -41577.1678 -41710.0026 -37448.7338 -36544.2726 13 14 15 16 17 18 -40512.1430 -24528.5488 -25821.3059 -27984.0194 -26323.7630 -26335.1276 19 20 21 22 23 24 -26527.2446 -26874.5029 -27244.9860 -46834.2625 394424.2984 317859.2305 25 26 27 28 29 30 254052.6223 260797.1407 372709.1556 255171.7832 31952.3948 17014.5964 31 32 33 34 35 36 19366.7834 -42834.7503 -28966.4211 -40520.5888 -42600.3811 -40545.7029 37 38 39 40 41 42 -46278.7468 1037.7896 -35880.6942 -29524.8634 -42467.9516 -35317.8389 43 44 45 46 47 48 -44790.4854 -27844.4979 -27665.3906 -27299.5662 -44539.2310 -42763.1335 49 50 51 52 53 54 -34879.5544 -29361.1868 -43154.5830 -44528.5759 -31197.3421 -919.5405 55 56 57 58 59 60 -13142.3500 -43539.2198 -48144.7547 -44173.6566 270077.6968 43016.1017 61 62 63 64 65 66 -5216.6725 135400.1381 -25124.3138 57554.8522 54770.8361 37351.2850 67 68 69 70 71 72 54472.9698 19371.4437 -6325.0237 -44477.0637 113682.6562 -42867.0989 73 74 75 76 77 78 -23058.2195 -25354.0345 -24108.8591 -15058.4672 15745.4635 96908.0479 79 80 81 82 83 84 -22103.7665 -27354.3574 307864.9365 255314.7557 19367.7117 -41405.8465 85 86 87 88 89 90 -31306.6462 -44797.8953 -39340.1469 -44954.4951 -3237.3400 57272.3226 91 92 93 94 95 96 -44674.6226 -24228.1029 -24728.5238 -27270.0839 -45857.6396 -45577.7990 97 98 99 100 101 102 -47329.2162 -39456.8467 -41753.6741 -45409.9825 -39775.0284 -52474.6781 103 104 105 106 107 108 261954.5651 1572683.1017 233148.7364 257829.4866 295179.2298 144444.5782 109 110 111 112 113 114 255189.2310 28173.9039 19368.2609 -42805.4550 -36881.3796 -23651.2955 115 116 117 118 119 120 -26076.5600 -32568.1605 -46433.0715 -38176.5027 19667.3742 -3215.1667 121 122 123 124 125 126 -21081.3276 -27190.5768 -46144.9875 263947.5932 257903.4131 274671.0640 127 128 129 130 131 132 255578.0824 19369.1730 -45089.1646 -25832.9709 -28040.8845 12265.4063 133 134 135 136 137 138 -2121.8332 -42954.3339 13462.8850 -31495.6327 -14172.0842 -21462.6264 139 140 141 142 143 144 -50553.3393 8573.1681 -6845.0189 4526.0989 -2057.2294 -42948.1585 145 146 147 148 149 150 -16622.1098 -13980.9615 -20027.9624 -25965.7807 -21680.0988 -17503.3726 151 152 153 154 155 156 -21446.2447 -45728.8460 -45255.4285 -50579.3918 4613.4976 1456.7309 157 158 159 160 161 162 -1671.0463 -28086.2722 -41456.4567 -14135.4520 -5840.9618 -14366.3933 163 164 165 166 167 168 642.2519 15156.7720 -28218.5650 -24800.7800 -10503.7408 -6014.6166 169 170 171 172 173 174 -14149.9627 -19783.1928 -21540.2681 -21344.7612 -45477.7571 -47836.2339 175 176 177 178 179 180 8757.8835 -42946.9816 48916.7561 -11086.9200 -11227.5401 -28353.4309 181 182 183 184 185 186 -14083.2630 -8297.8835 -21283.3830 -21478.2873 -49365.8792 -49986.4445 187 188 189 190 191 192 54960.0857 -14325.5584 -1451.0215 -43041.9516 43475.0731 -14095.3030 193 194 195 196 197 198 -21457.0397 -51631.5519 36954.5370 20021.0151 -2423.2971 -36458.8656 199 200 201 202 203 204 -42270.6324 -17122.5471 -39287.5223 -41195.7133 -41966.6224 -9224.0604 205 206 207 208 209 210 -39703.5078 -41317.8069 -41228.4986 -11786.2626 -1058.3495 -11404.5677 211 212 213 214 215 216 -14046.0324 -23008.6023 -22262.8924 -22264.0542 -10276.8540 -21287.9087 217 218 219 220 221 222 -33671.3010 -51342.0042 -40400.7039 -31302.0189 -51827.9467 -43920.4483 223 224 225 226 227 228 8858.5018 8254.8441 -6193.0124 -2261.6813 -42972.0650 -9501.6653 229 230 231 232 233 234 5328.4893 -11897.7741 -14145.7387 -20655.2260 -29433.4309 -13638.7970 235 236 237 238 239 240 -21438.5047 -50930.2093 8762.6115 -497.3053 -38773.8304 -41472.7734 241 242 243 244 245 246 -14084.0471 -21452.3119 8766.2970 -42952.6808 33249.1059 -8610.9119 247 248 249 250 251 252 34640.0332 3408.5755 -19117.5643 -8886.1021 -24503.5809 -14050.3900 253 254 255 256 257 258 -22368.6266 -21446.6350 -49999.7310 -39674.3499 8756.4385 -40754.4781 259 260 261 262 263 264 -41373.1437 -39504.6980 -43861.4218 -41520.0628 -42340.7744 -42283.8891 265 266 267 268 269 270 -42154.1613 -41418.8906 -41190.4944 -41522.7957 -16106.5271 70213.6112 271 272 273 274 275 276 -12965.7509 -14252.4176 -9762.3188 28326.6824 -72197.8802 -21452.9153 277 278 279 280 281 282 -50498.1878 12404.7357 -1523.0664 -1854.3021 5584.5857 -35511.0758 283 284 285 286 287 288 -41448.1311 -14353.3936 49896.4346 -15933.3086 -16502.9445 -15842.8317 289 290 291 292 293 294 -16911.7657 -4342.8956 -8187.1311 -21499.7920 8760.4757 16265.5438 295 296 297 298 299 300 -38678.7183 -40741.9988 -41513.9005 -39824.0098 -39528.8470 58580.3827 301 302 303 304 305 306 -41856.3529 -41504.9705 -22330.2434 -6385.1142 -12527.7453 5355.5399 307 308 309 310 311 312 -14283.6774 -9538.7591 -21436.5161 -43078.4997 87610.1064 -54810.7466 313 314 315 316 317 318 -50672.2440 -52783.3964 8755.3847 -8680.0354 -37876.0514 -41478.2398 319 320 321 322 323 324 -7438.8236 -6101.6914 -14094.9711 -21214.9683 -48079.1703 -51875.8077 325 326 327 328 329 330 38832.0789 -2844.9563 -39215.4444 -41242.5948 -48228.0335 -41468.5709 331 332 333 334 335 336 -15124.5938 -14183.2363 -21458.0408 -48789.3135 19713.6160 -1903.7392 337 338 339 340 341 342 -42976.1144 -14123.6980 -21293.0822 -43749.0207 -48908.9387 -52577.5135 343 344 345 346 347 348 -51717.4612 -49636.8821 -50546.0696 -37627.9805 -51854.7547 -49788.5143 349 350 351 352 353 354 15840.8768 51414.1747 -4374.0629 241.7956 17621.4624 181.3042 355 356 357 358 359 360 37760.0776 27797.2372 35603.5300 -1544.5919 -990.1866 -39439.9441 361 362 363 364 365 366 -40600.8492 -38590.3090 -40687.3700 -41530.3466 -41872.8288 -40339.2588 367 368 369 370 371 372 -41432.7408 -14132.9224 -21475.7677 8754.7509 -34028.1065 -41735.1039 373 374 375 376 377 378 -41346.6164 -14143.8862 -21448.1433 7729.7929 40343.7547 46018.7664 379 380 381 382 383 384 25959.6806 -8497.9067 62712.9737 -709.0604 -1835.9757 -1675.7406 385 386 387 388 389 390 24880.5713 37943.2659 66802.8471 -141.0171 -1818.1783 -42998.9352 391 392 393 394 395 396 -25722.1518 -9579.2388 -13423.9607 14948.5575 -21352.8312 104425.8356 397 398 399 400 401 402 -14020.8518 -18119.8950 -18812.5436 168783.5740 -19590.6876 -24061.0828 403 404 405 406 407 408 -61123.6847 -21453.1774 -49404.0111 -48253.5544 4818.9750 -1878.3993 409 410 411 412 413 414 32358.0517 -2858.9143 15164.1944 72250.8327 -325.2517 -3285.2804 415 416 417 418 419 420 -1915.0474 3736.8753 -3596.8006 38305.7594 37612.1963 -2670.6788 421 422 423 424 425 426 -42142.3471 -41220.1957 -41415.3922 -41942.1361 -41601.6787 -41443.5895 427 428 429 430 431 -41299.7036 -41479.2746 -16823.2574 -13978.7535 -1711.4855 > postscript(file="/var/www/html/rcomp/tmp/6t3251291380013.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 = 431 Frequency = 1 lag(myerror, k = 1) myerror 0 -43371.7534 NA 1 -34371.7678 -43371.7534 2 -40486.6638 -34371.7678 3 -37562.3114 -40486.6638 4 -37007.6557 -37562.3114 5 -44482.4940 -37007.6557 6 -43355.1770 -44482.4940 7 -41865.3726 -43355.1770 8 -41577.1678 -41865.3726 9 -41710.0026 -41577.1678 10 -37448.7338 -41710.0026 11 -36544.2726 -37448.7338 12 -40512.1430 -36544.2726 13 -24528.5488 -40512.1430 14 -25821.3059 -24528.5488 15 -27984.0194 -25821.3059 16 -26323.7630 -27984.0194 17 -26335.1276 -26323.7630 18 -26527.2446 -26335.1276 19 -26874.5029 -26527.2446 20 -27244.9860 -26874.5029 21 -46834.2625 -27244.9860 22 394424.2984 -46834.2625 23 317859.2305 394424.2984 24 254052.6223 317859.2305 25 260797.1407 254052.6223 26 372709.1556 260797.1407 27 255171.7832 372709.1556 28 31952.3948 255171.7832 29 17014.5964 31952.3948 30 19366.7834 17014.5964 31 -42834.7503 19366.7834 32 -28966.4211 -42834.7503 33 -40520.5888 -28966.4211 34 -42600.3811 -40520.5888 35 -40545.7029 -42600.3811 36 -46278.7468 -40545.7029 37 1037.7896 -46278.7468 38 -35880.6942 1037.7896 39 -29524.8634 -35880.6942 40 -42467.9516 -29524.8634 41 -35317.8389 -42467.9516 42 -44790.4854 -35317.8389 43 -27844.4979 -44790.4854 44 -27665.3906 -27844.4979 45 -27299.5662 -27665.3906 46 -44539.2310 -27299.5662 47 -42763.1335 -44539.2310 48 -34879.5544 -42763.1335 49 -29361.1868 -34879.5544 50 -43154.5830 -29361.1868 51 -44528.5759 -43154.5830 52 -31197.3421 -44528.5759 53 -919.5405 -31197.3421 54 -13142.3500 -919.5405 55 -43539.2198 -13142.3500 56 -48144.7547 -43539.2198 57 -44173.6566 -48144.7547 58 270077.6968 -44173.6566 59 43016.1017 270077.6968 60 -5216.6725 43016.1017 61 135400.1381 -5216.6725 62 -25124.3138 135400.1381 63 57554.8522 -25124.3138 64 54770.8361 57554.8522 65 37351.2850 54770.8361 66 54472.9698 37351.2850 67 19371.4437 54472.9698 68 -6325.0237 19371.4437 69 -44477.0637 -6325.0237 70 113682.6562 -44477.0637 71 -42867.0989 113682.6562 72 -23058.2195 -42867.0989 73 -25354.0345 -23058.2195 74 -24108.8591 -25354.0345 75 -15058.4672 -24108.8591 76 15745.4635 -15058.4672 77 96908.0479 15745.4635 78 -22103.7665 96908.0479 79 -27354.3574 -22103.7665 80 307864.9365 -27354.3574 81 255314.7557 307864.9365 82 19367.7117 255314.7557 83 -41405.8465 19367.7117 84 -31306.6462 -41405.8465 85 -44797.8953 -31306.6462 86 -39340.1469 -44797.8953 87 -44954.4951 -39340.1469 88 -3237.3400 -44954.4951 89 57272.3226 -3237.3400 90 -44674.6226 57272.3226 91 -24228.1029 -44674.6226 92 -24728.5238 -24228.1029 93 -27270.0839 -24728.5238 94 -45857.6396 -27270.0839 95 -45577.7990 -45857.6396 96 -47329.2162 -45577.7990 97 -39456.8467 -47329.2162 98 -41753.6741 -39456.8467 99 -45409.9825 -41753.6741 100 -39775.0284 -45409.9825 101 -52474.6781 -39775.0284 102 261954.5651 -52474.6781 103 1572683.1017 261954.5651 104 233148.7364 1572683.1017 105 257829.4866 233148.7364 106 295179.2298 257829.4866 107 144444.5782 295179.2298 108 255189.2310 144444.5782 109 28173.9039 255189.2310 110 19368.2609 28173.9039 111 -42805.4550 19368.2609 112 -36881.3796 -42805.4550 113 -23651.2955 -36881.3796 114 -26076.5600 -23651.2955 115 -32568.1605 -26076.5600 116 -46433.0715 -32568.1605 117 -38176.5027 -46433.0715 118 19667.3742 -38176.5027 119 -3215.1667 19667.3742 120 -21081.3276 -3215.1667 121 -27190.5768 -21081.3276 122 -46144.9875 -27190.5768 123 263947.5932 -46144.9875 124 257903.4131 263947.5932 125 274671.0640 257903.4131 126 255578.0824 274671.0640 127 19369.1730 255578.0824 128 -45089.1646 19369.1730 129 -25832.9709 -45089.1646 130 -28040.8845 -25832.9709 131 12265.4063 -28040.8845 132 -2121.8332 12265.4063 133 -42954.3339 -2121.8332 134 13462.8850 -42954.3339 135 -31495.6327 13462.8850 136 -14172.0842 -31495.6327 137 -21462.6264 -14172.0842 138 -50553.3393 -21462.6264 139 8573.1681 -50553.3393 140 -6845.0189 8573.1681 141 4526.0989 -6845.0189 142 -2057.2294 4526.0989 143 -42948.1585 -2057.2294 144 -16622.1098 -42948.1585 145 -13980.9615 -16622.1098 146 -20027.9624 -13980.9615 147 -25965.7807 -20027.9624 148 -21680.0988 -25965.7807 149 -17503.3726 -21680.0988 150 -21446.2447 -17503.3726 151 -45728.8460 -21446.2447 152 -45255.4285 -45728.8460 153 -50579.3918 -45255.4285 154 4613.4976 -50579.3918 155 1456.7309 4613.4976 156 -1671.0463 1456.7309 157 -28086.2722 -1671.0463 158 -41456.4567 -28086.2722 159 -14135.4520 -41456.4567 160 -5840.9618 -14135.4520 161 -14366.3933 -5840.9618 162 642.2519 -14366.3933 163 15156.7720 642.2519 164 -28218.5650 15156.7720 165 -24800.7800 -28218.5650 166 -10503.7408 -24800.7800 167 -6014.6166 -10503.7408 168 -14149.9627 -6014.6166 169 -19783.1928 -14149.9627 170 -21540.2681 -19783.1928 171 -21344.7612 -21540.2681 172 -45477.7571 -21344.7612 173 -47836.2339 -45477.7571 174 8757.8835 -47836.2339 175 -42946.9816 8757.8835 176 48916.7561 -42946.9816 177 -11086.9200 48916.7561 178 -11227.5401 -11086.9200 179 -28353.4309 -11227.5401 180 -14083.2630 -28353.4309 181 -8297.8835 -14083.2630 182 -21283.3830 -8297.8835 183 -21478.2873 -21283.3830 184 -49365.8792 -21478.2873 185 -49986.4445 -49365.8792 186 54960.0857 -49986.4445 187 -14325.5584 54960.0857 188 -1451.0215 -14325.5584 189 -43041.9516 -1451.0215 190 43475.0731 -43041.9516 191 -14095.3030 43475.0731 192 -21457.0397 -14095.3030 193 -51631.5519 -21457.0397 194 36954.5370 -51631.5519 195 20021.0151 36954.5370 196 -2423.2971 20021.0151 197 -36458.8656 -2423.2971 198 -42270.6324 -36458.8656 199 -17122.5471 -42270.6324 200 -39287.5223 -17122.5471 201 -41195.7133 -39287.5223 202 -41966.6224 -41195.7133 203 -9224.0604 -41966.6224 204 -39703.5078 -9224.0604 205 -41317.8069 -39703.5078 206 -41228.4986 -41317.8069 207 -11786.2626 -41228.4986 208 -1058.3495 -11786.2626 209 -11404.5677 -1058.3495 210 -14046.0324 -11404.5677 211 -23008.6023 -14046.0324 212 -22262.8924 -23008.6023 213 -22264.0542 -22262.8924 214 -10276.8540 -22264.0542 215 -21287.9087 -10276.8540 216 -33671.3010 -21287.9087 217 -51342.0042 -33671.3010 218 -40400.7039 -51342.0042 219 -31302.0189 -40400.7039 220 -51827.9467 -31302.0189 221 -43920.4483 -51827.9467 222 8858.5018 -43920.4483 223 8254.8441 8858.5018 224 -6193.0124 8254.8441 225 -2261.6813 -6193.0124 226 -42972.0650 -2261.6813 227 -9501.6653 -42972.0650 228 5328.4893 -9501.6653 229 -11897.7741 5328.4893 230 -14145.7387 -11897.7741 231 -20655.2260 -14145.7387 232 -29433.4309 -20655.2260 233 -13638.7970 -29433.4309 234 -21438.5047 -13638.7970 235 -50930.2093 -21438.5047 236 8762.6115 -50930.2093 237 -497.3053 8762.6115 238 -38773.8304 -497.3053 239 -41472.7734 -38773.8304 240 -14084.0471 -41472.7734 241 -21452.3119 -14084.0471 242 8766.2970 -21452.3119 243 -42952.6808 8766.2970 244 33249.1059 -42952.6808 245 -8610.9119 33249.1059 246 34640.0332 -8610.9119 247 3408.5755 34640.0332 248 -19117.5643 3408.5755 249 -8886.1021 -19117.5643 250 -24503.5809 -8886.1021 251 -14050.3900 -24503.5809 252 -22368.6266 -14050.3900 253 -21446.6350 -22368.6266 254 -49999.7310 -21446.6350 255 -39674.3499 -49999.7310 256 8756.4385 -39674.3499 257 -40754.4781 8756.4385 258 -41373.1437 -40754.4781 259 -39504.6980 -41373.1437 260 -43861.4218 -39504.6980 261 -41520.0628 -43861.4218 262 -42340.7744 -41520.0628 263 -42283.8891 -42340.7744 264 -42154.1613 -42283.8891 265 -41418.8906 -42154.1613 266 -41190.4944 -41418.8906 267 -41522.7957 -41190.4944 268 -16106.5271 -41522.7957 269 70213.6112 -16106.5271 270 -12965.7509 70213.6112 271 -14252.4176 -12965.7509 272 -9762.3188 -14252.4176 273 28326.6824 -9762.3188 274 -72197.8802 28326.6824 275 -21452.9153 -72197.8802 276 -50498.1878 -21452.9153 277 12404.7357 -50498.1878 278 -1523.0664 12404.7357 279 -1854.3021 -1523.0664 280 5584.5857 -1854.3021 281 -35511.0758 5584.5857 282 -41448.1311 -35511.0758 283 -14353.3936 -41448.1311 284 49896.4346 -14353.3936 285 -15933.3086 49896.4346 286 -16502.9445 -15933.3086 287 -15842.8317 -16502.9445 288 -16911.7657 -15842.8317 289 -4342.8956 -16911.7657 290 -8187.1311 -4342.8956 291 -21499.7920 -8187.1311 292 8760.4757 -21499.7920 293 16265.5438 8760.4757 294 -38678.7183 16265.5438 295 -40741.9988 -38678.7183 296 -41513.9005 -40741.9988 297 -39824.0098 -41513.9005 298 -39528.8470 -39824.0098 299 58580.3827 -39528.8470 300 -41856.3529 58580.3827 301 -41504.9705 -41856.3529 302 -22330.2434 -41504.9705 303 -6385.1142 -22330.2434 304 -12527.7453 -6385.1142 305 5355.5399 -12527.7453 306 -14283.6774 5355.5399 307 -9538.7591 -14283.6774 308 -21436.5161 -9538.7591 309 -43078.4997 -21436.5161 310 87610.1064 -43078.4997 311 -54810.7466 87610.1064 312 -50672.2440 -54810.7466 313 -52783.3964 -50672.2440 314 8755.3847 -52783.3964 315 -8680.0354 8755.3847 316 -37876.0514 -8680.0354 317 -41478.2398 -37876.0514 318 -7438.8236 -41478.2398 319 -6101.6914 -7438.8236 320 -14094.9711 -6101.6914 321 -21214.9683 -14094.9711 322 -48079.1703 -21214.9683 323 -51875.8077 -48079.1703 324 38832.0789 -51875.8077 325 -2844.9563 38832.0789 326 -39215.4444 -2844.9563 327 -41242.5948 -39215.4444 328 -48228.0335 -41242.5948 329 -41468.5709 -48228.0335 330 -15124.5938 -41468.5709 331 -14183.2363 -15124.5938 332 -21458.0408 -14183.2363 333 -48789.3135 -21458.0408 334 19713.6160 -48789.3135 335 -1903.7392 19713.6160 336 -42976.1144 -1903.7392 337 -14123.6980 -42976.1144 338 -21293.0822 -14123.6980 339 -43749.0207 -21293.0822 340 -48908.9387 -43749.0207 341 -52577.5135 -48908.9387 342 -51717.4612 -52577.5135 343 -49636.8821 -51717.4612 344 -50546.0696 -49636.8821 345 -37627.9805 -50546.0696 346 -51854.7547 -37627.9805 347 -49788.5143 -51854.7547 348 15840.8768 -49788.5143 349 51414.1747 15840.8768 350 -4374.0629 51414.1747 351 241.7956 -4374.0629 352 17621.4624 241.7956 353 181.3042 17621.4624 354 37760.0776 181.3042 355 27797.2372 37760.0776 356 35603.5300 27797.2372 357 -1544.5919 35603.5300 358 -990.1866 -1544.5919 359 -39439.9441 -990.1866 360 -40600.8492 -39439.9441 361 -38590.3090 -40600.8492 362 -40687.3700 -38590.3090 363 -41530.3466 -40687.3700 364 -41872.8288 -41530.3466 365 -40339.2588 -41872.8288 366 -41432.7408 -40339.2588 367 -14132.9224 -41432.7408 368 -21475.7677 -14132.9224 369 8754.7509 -21475.7677 370 -34028.1065 8754.7509 371 -41735.1039 -34028.1065 372 -41346.6164 -41735.1039 373 -14143.8862 -41346.6164 374 -21448.1433 -14143.8862 375 7729.7929 -21448.1433 376 40343.7547 7729.7929 377 46018.7664 40343.7547 378 25959.6806 46018.7664 379 -8497.9067 25959.6806 380 62712.9737 -8497.9067 381 -709.0604 62712.9737 382 -1835.9757 -709.0604 383 -1675.7406 -1835.9757 384 24880.5713 -1675.7406 385 37943.2659 24880.5713 386 66802.8471 37943.2659 387 -141.0171 66802.8471 388 -1818.1783 -141.0171 389 -42998.9352 -1818.1783 390 -25722.1518 -42998.9352 391 -9579.2388 -25722.1518 392 -13423.9607 -9579.2388 393 14948.5575 -13423.9607 394 -21352.8312 14948.5575 395 104425.8356 -21352.8312 396 -14020.8518 104425.8356 397 -18119.8950 -14020.8518 398 -18812.5436 -18119.8950 399 168783.5740 -18812.5436 400 -19590.6876 168783.5740 401 -24061.0828 -19590.6876 402 -61123.6847 -24061.0828 403 -21453.1774 -61123.6847 404 -49404.0111 -21453.1774 405 -48253.5544 -49404.0111 406 4818.9750 -48253.5544 407 -1878.3993 4818.9750 408 32358.0517 -1878.3993 409 -2858.9143 32358.0517 410 15164.1944 -2858.9143 411 72250.8327 15164.1944 412 -325.2517 72250.8327 413 -3285.2804 -325.2517 414 -1915.0474 -3285.2804 415 3736.8753 -1915.0474 416 -3596.8006 3736.8753 417 38305.7594 -3596.8006 418 37612.1963 38305.7594 419 -2670.6788 37612.1963 420 -42142.3471 -2670.6788 421 -41220.1957 -42142.3471 422 -41415.3922 -41220.1957 423 -41942.1361 -41415.3922 424 -41601.6787 -41942.1361 425 -41443.5895 -41601.6787 426 -41299.7036 -41443.5895 427 -41479.2746 -41299.7036 428 -16823.2574 -41479.2746 429 -13978.7535 -16823.2574 430 -1711.4855 -13978.7535 431 NA -1711.4855 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -34371.7678 -43371.7534 [2,] -40486.6638 -34371.7678 [3,] -37562.3114 -40486.6638 [4,] -37007.6557 -37562.3114 [5,] -44482.4940 -37007.6557 [6,] -43355.1770 -44482.4940 [7,] -41865.3726 -43355.1770 [8,] -41577.1678 -41865.3726 [9,] -41710.0026 -41577.1678 [10,] -37448.7338 -41710.0026 [11,] -36544.2726 -37448.7338 [12,] -40512.1430 -36544.2726 [13,] -24528.5488 -40512.1430 [14,] -25821.3059 -24528.5488 [15,] -27984.0194 -25821.3059 [16,] -26323.7630 -27984.0194 [17,] -26335.1276 -26323.7630 [18,] -26527.2446 -26335.1276 [19,] -26874.5029 -26527.2446 [20,] -27244.9860 -26874.5029 [21,] -46834.2625 -27244.9860 [22,] 394424.2984 -46834.2625 [23,] 317859.2305 394424.2984 [24,] 254052.6223 317859.2305 [25,] 260797.1407 254052.6223 [26,] 372709.1556 260797.1407 [27,] 255171.7832 372709.1556 [28,] 31952.3948 255171.7832 [29,] 17014.5964 31952.3948 [30,] 19366.7834 17014.5964 [31,] -42834.7503 19366.7834 [32,] -28966.4211 -42834.7503 [33,] -40520.5888 -28966.4211 [34,] -42600.3811 -40520.5888 [35,] -40545.7029 -42600.3811 [36,] -46278.7468 -40545.7029 [37,] 1037.7896 -46278.7468 [38,] -35880.6942 1037.7896 [39,] -29524.8634 -35880.6942 [40,] -42467.9516 -29524.8634 [41,] -35317.8389 -42467.9516 [42,] -44790.4854 -35317.8389 [43,] -27844.4979 -44790.4854 [44,] -27665.3906 -27844.4979 [45,] -27299.5662 -27665.3906 [46,] -44539.2310 -27299.5662 [47,] -42763.1335 -44539.2310 [48,] -34879.5544 -42763.1335 [49,] -29361.1868 -34879.5544 [50,] -43154.5830 -29361.1868 [51,] -44528.5759 -43154.5830 [52,] -31197.3421 -44528.5759 [53,] -919.5405 -31197.3421 [54,] -13142.3500 -919.5405 [55,] -43539.2198 -13142.3500 [56,] -48144.7547 -43539.2198 [57,] -44173.6566 -48144.7547 [58,] 270077.6968 -44173.6566 [59,] 43016.1017 270077.6968 [60,] -5216.6725 43016.1017 [61,] 135400.1381 -5216.6725 [62,] -25124.3138 135400.1381 [63,] 57554.8522 -25124.3138 [64,] 54770.8361 57554.8522 [65,] 37351.2850 54770.8361 [66,] 54472.9698 37351.2850 [67,] 19371.4437 54472.9698 [68,] -6325.0237 19371.4437 [69,] -44477.0637 -6325.0237 [70,] 113682.6562 -44477.0637 [71,] -42867.0989 113682.6562 [72,] -23058.2195 -42867.0989 [73,] -25354.0345 -23058.2195 [74,] -24108.8591 -25354.0345 [75,] -15058.4672 -24108.8591 [76,] 15745.4635 -15058.4672 [77,] 96908.0479 15745.4635 [78,] -22103.7665 96908.0479 [79,] -27354.3574 -22103.7665 [80,] 307864.9365 -27354.3574 [81,] 255314.7557 307864.9365 [82,] 19367.7117 255314.7557 [83,] -41405.8465 19367.7117 [84,] -31306.6462 -41405.8465 [85,] -44797.8953 -31306.6462 [86,] -39340.1469 -44797.8953 [87,] -44954.4951 -39340.1469 [88,] -3237.3400 -44954.4951 [89,] 57272.3226 -3237.3400 [90,] -44674.6226 57272.3226 [91,] -24228.1029 -44674.6226 [92,] -24728.5238 -24228.1029 [93,] -27270.0839 -24728.5238 [94,] -45857.6396 -27270.0839 [95,] -45577.7990 -45857.6396 [96,] -47329.2162 -45577.7990 [97,] -39456.8467 -47329.2162 [98,] -41753.6741 -39456.8467 [99,] -45409.9825 -41753.6741 [100,] -39775.0284 -45409.9825 [101,] -52474.6781 -39775.0284 [102,] 261954.5651 -52474.6781 [103,] 1572683.1017 261954.5651 [104,] 233148.7364 1572683.1017 [105,] 257829.4866 233148.7364 [106,] 295179.2298 257829.4866 [107,] 144444.5782 295179.2298 [108,] 255189.2310 144444.5782 [109,] 28173.9039 255189.2310 [110,] 19368.2609 28173.9039 [111,] -42805.4550 19368.2609 [112,] -36881.3796 -42805.4550 [113,] -23651.2955 -36881.3796 [114,] -26076.5600 -23651.2955 [115,] -32568.1605 -26076.5600 [116,] -46433.0715 -32568.1605 [117,] -38176.5027 -46433.0715 [118,] 19667.3742 -38176.5027 [119,] -3215.1667 19667.3742 [120,] -21081.3276 -3215.1667 [121,] -27190.5768 -21081.3276 [122,] -46144.9875 -27190.5768 [123,] 263947.5932 -46144.9875 [124,] 257903.4131 263947.5932 [125,] 274671.0640 257903.4131 [126,] 255578.0824 274671.0640 [127,] 19369.1730 255578.0824 [128,] -45089.1646 19369.1730 [129,] -25832.9709 -45089.1646 [130,] -28040.8845 -25832.9709 [131,] 12265.4063 -28040.8845 [132,] -2121.8332 12265.4063 [133,] -42954.3339 -2121.8332 [134,] 13462.8850 -42954.3339 [135,] -31495.6327 13462.8850 [136,] -14172.0842 -31495.6327 [137,] -21462.6264 -14172.0842 [138,] -50553.3393 -21462.6264 [139,] 8573.1681 -50553.3393 [140,] -6845.0189 8573.1681 [141,] 4526.0989 -6845.0189 [142,] -2057.2294 4526.0989 [143,] -42948.1585 -2057.2294 [144,] -16622.1098 -42948.1585 [145,] -13980.9615 -16622.1098 [146,] -20027.9624 -13980.9615 [147,] -25965.7807 -20027.9624 [148,] -21680.0988 -25965.7807 [149,] -17503.3726 -21680.0988 [150,] -21446.2447 -17503.3726 [151,] -45728.8460 -21446.2447 [152,] -45255.4285 -45728.8460 [153,] -50579.3918 -45255.4285 [154,] 4613.4976 -50579.3918 [155,] 1456.7309 4613.4976 [156,] -1671.0463 1456.7309 [157,] -28086.2722 -1671.0463 [158,] -41456.4567 -28086.2722 [159,] -14135.4520 -41456.4567 [160,] -5840.9618 -14135.4520 [161,] -14366.3933 -5840.9618 [162,] 642.2519 -14366.3933 [163,] 15156.7720 642.2519 [164,] -28218.5650 15156.7720 [165,] -24800.7800 -28218.5650 [166,] -10503.7408 -24800.7800 [167,] -6014.6166 -10503.7408 [168,] -14149.9627 -6014.6166 [169,] -19783.1928 -14149.9627 [170,] -21540.2681 -19783.1928 [171,] -21344.7612 -21540.2681 [172,] -45477.7571 -21344.7612 [173,] -47836.2339 -45477.7571 [174,] 8757.8835 -47836.2339 [175,] -42946.9816 8757.8835 [176,] 48916.7561 -42946.9816 [177,] -11086.9200 48916.7561 [178,] -11227.5401 -11086.9200 [179,] -28353.4309 -11227.5401 [180,] -14083.2630 -28353.4309 [181,] -8297.8835 -14083.2630 [182,] -21283.3830 -8297.8835 [183,] -21478.2873 -21283.3830 [184,] -49365.8792 -21478.2873 [185,] -49986.4445 -49365.8792 [186,] 54960.0857 -49986.4445 [187,] -14325.5584 54960.0857 [188,] -1451.0215 -14325.5584 [189,] -43041.9516 -1451.0215 [190,] 43475.0731 -43041.9516 [191,] -14095.3030 43475.0731 [192,] -21457.0397 -14095.3030 [193,] -51631.5519 -21457.0397 [194,] 36954.5370 -51631.5519 [195,] 20021.0151 36954.5370 [196,] -2423.2971 20021.0151 [197,] -36458.8656 -2423.2971 [198,] -42270.6324 -36458.8656 [199,] -17122.5471 -42270.6324 [200,] -39287.5223 -17122.5471 [201,] -41195.7133 -39287.5223 [202,] -41966.6224 -41195.7133 [203,] -9224.0604 -41966.6224 [204,] -39703.5078 -9224.0604 [205,] -41317.8069 -39703.5078 [206,] -41228.4986 -41317.8069 [207,] -11786.2626 -41228.4986 [208,] -1058.3495 -11786.2626 [209,] -11404.5677 -1058.3495 [210,] -14046.0324 -11404.5677 [211,] -23008.6023 -14046.0324 [212,] -22262.8924 -23008.6023 [213,] -22264.0542 -22262.8924 [214,] -10276.8540 -22264.0542 [215,] -21287.9087 -10276.8540 [216,] -33671.3010 -21287.9087 [217,] -51342.0042 -33671.3010 [218,] -40400.7039 -51342.0042 [219,] -31302.0189 -40400.7039 [220,] -51827.9467 -31302.0189 [221,] -43920.4483 -51827.9467 [222,] 8858.5018 -43920.4483 [223,] 8254.8441 8858.5018 [224,] -6193.0124 8254.8441 [225,] -2261.6813 -6193.0124 [226,] -42972.0650 -2261.6813 [227,] -9501.6653 -42972.0650 [228,] 5328.4893 -9501.6653 [229,] -11897.7741 5328.4893 [230,] -14145.7387 -11897.7741 [231,] -20655.2260 -14145.7387 [232,] -29433.4309 -20655.2260 [233,] -13638.7970 -29433.4309 [234,] -21438.5047 -13638.7970 [235,] -50930.2093 -21438.5047 [236,] 8762.6115 -50930.2093 [237,] -497.3053 8762.6115 [238,] -38773.8304 -497.3053 [239,] -41472.7734 -38773.8304 [240,] -14084.0471 -41472.7734 [241,] -21452.3119 -14084.0471 [242,] 8766.2970 -21452.3119 [243,] -42952.6808 8766.2970 [244,] 33249.1059 -42952.6808 [245,] -8610.9119 33249.1059 [246,] 34640.0332 -8610.9119 [247,] 3408.5755 34640.0332 [248,] -19117.5643 3408.5755 [249,] -8886.1021 -19117.5643 [250,] -24503.5809 -8886.1021 [251,] -14050.3900 -24503.5809 [252,] -22368.6266 -14050.3900 [253,] -21446.6350 -22368.6266 [254,] -49999.7310 -21446.6350 [255,] -39674.3499 -49999.7310 [256,] 8756.4385 -39674.3499 [257,] -40754.4781 8756.4385 [258,] -41373.1437 -40754.4781 [259,] -39504.6980 -41373.1437 [260,] -43861.4218 -39504.6980 [261,] -41520.0628 -43861.4218 [262,] -42340.7744 -41520.0628 [263,] -42283.8891 -42340.7744 [264,] -42154.1613 -42283.8891 [265,] -41418.8906 -42154.1613 [266,] -41190.4944 -41418.8906 [267,] -41522.7957 -41190.4944 [268,] -16106.5271 -41522.7957 [269,] 70213.6112 -16106.5271 [270,] -12965.7509 70213.6112 [271,] -14252.4176 -12965.7509 [272,] -9762.3188 -14252.4176 [273,] 28326.6824 -9762.3188 [274,] -72197.8802 28326.6824 [275,] -21452.9153 -72197.8802 [276,] -50498.1878 -21452.9153 [277,] 12404.7357 -50498.1878 [278,] -1523.0664 12404.7357 [279,] -1854.3021 -1523.0664 [280,] 5584.5857 -1854.3021 [281,] -35511.0758 5584.5857 [282,] -41448.1311 -35511.0758 [283,] -14353.3936 -41448.1311 [284,] 49896.4346 -14353.3936 [285,] -15933.3086 49896.4346 [286,] -16502.9445 -15933.3086 [287,] -15842.8317 -16502.9445 [288,] -16911.7657 -15842.8317 [289,] -4342.8956 -16911.7657 [290,] -8187.1311 -4342.8956 [291,] -21499.7920 -8187.1311 [292,] 8760.4757 -21499.7920 [293,] 16265.5438 8760.4757 [294,] -38678.7183 16265.5438 [295,] -40741.9988 -38678.7183 [296,] -41513.9005 -40741.9988 [297,] -39824.0098 -41513.9005 [298,] -39528.8470 -39824.0098 [299,] 58580.3827 -39528.8470 [300,] -41856.3529 58580.3827 [301,] -41504.9705 -41856.3529 [302,] -22330.2434 -41504.9705 [303,] -6385.1142 -22330.2434 [304,] -12527.7453 -6385.1142 [305,] 5355.5399 -12527.7453 [306,] -14283.6774 5355.5399 [307,] -9538.7591 -14283.6774 [308,] -21436.5161 -9538.7591 [309,] -43078.4997 -21436.5161 [310,] 87610.1064 -43078.4997 [311,] -54810.7466 87610.1064 [312,] -50672.2440 -54810.7466 [313,] -52783.3964 -50672.2440 [314,] 8755.3847 -52783.3964 [315,] -8680.0354 8755.3847 [316,] -37876.0514 -8680.0354 [317,] -41478.2398 -37876.0514 [318,] -7438.8236 -41478.2398 [319,] -6101.6914 -7438.8236 [320,] -14094.9711 -6101.6914 [321,] -21214.9683 -14094.9711 [322,] -48079.1703 -21214.9683 [323,] -51875.8077 -48079.1703 [324,] 38832.0789 -51875.8077 [325,] -2844.9563 38832.0789 [326,] -39215.4444 -2844.9563 [327,] -41242.5948 -39215.4444 [328,] -48228.0335 -41242.5948 [329,] -41468.5709 -48228.0335 [330,] -15124.5938 -41468.5709 [331,] -14183.2363 -15124.5938 [332,] -21458.0408 -14183.2363 [333,] -48789.3135 -21458.0408 [334,] 19713.6160 -48789.3135 [335,] -1903.7392 19713.6160 [336,] -42976.1144 -1903.7392 [337,] -14123.6980 -42976.1144 [338,] -21293.0822 -14123.6980 [339,] -43749.0207 -21293.0822 [340,] -48908.9387 -43749.0207 [341,] -52577.5135 -48908.9387 [342,] -51717.4612 -52577.5135 [343,] -49636.8821 -51717.4612 [344,] -50546.0696 -49636.8821 [345,] -37627.9805 -50546.0696 [346,] -51854.7547 -37627.9805 [347,] -49788.5143 -51854.7547 [348,] 15840.8768 -49788.5143 [349,] 51414.1747 15840.8768 [350,] -4374.0629 51414.1747 [351,] 241.7956 -4374.0629 [352,] 17621.4624 241.7956 [353,] 181.3042 17621.4624 [354,] 37760.0776 181.3042 [355,] 27797.2372 37760.0776 [356,] 35603.5300 27797.2372 [357,] -1544.5919 35603.5300 [358,] -990.1866 -1544.5919 [359,] -39439.9441 -990.1866 [360,] -40600.8492 -39439.9441 [361,] -38590.3090 -40600.8492 [362,] -40687.3700 -38590.3090 [363,] -41530.3466 -40687.3700 [364,] -41872.8288 -41530.3466 [365,] -40339.2588 -41872.8288 [366,] -41432.7408 -40339.2588 [367,] -14132.9224 -41432.7408 [368,] -21475.7677 -14132.9224 [369,] 8754.7509 -21475.7677 [370,] -34028.1065 8754.7509 [371,] -41735.1039 -34028.1065 [372,] -41346.6164 -41735.1039 [373,] -14143.8862 -41346.6164 [374,] -21448.1433 -14143.8862 [375,] 7729.7929 -21448.1433 [376,] 40343.7547 7729.7929 [377,] 46018.7664 40343.7547 [378,] 25959.6806 46018.7664 [379,] -8497.9067 25959.6806 [380,] 62712.9737 -8497.9067 [381,] -709.0604 62712.9737 [382,] -1835.9757 -709.0604 [383,] -1675.7406 -1835.9757 [384,] 24880.5713 -1675.7406 [385,] 37943.2659 24880.5713 [386,] 66802.8471 37943.2659 [387,] -141.0171 66802.8471 [388,] -1818.1783 -141.0171 [389,] -42998.9352 -1818.1783 [390,] -25722.1518 -42998.9352 [391,] -9579.2388 -25722.1518 [392,] -13423.9607 -9579.2388 [393,] 14948.5575 -13423.9607 [394,] -21352.8312 14948.5575 [395,] 104425.8356 -21352.8312 [396,] -14020.8518 104425.8356 [397,] -18119.8950 -14020.8518 [398,] -18812.5436 -18119.8950 [399,] 168783.5740 -18812.5436 [400,] -19590.6876 168783.5740 [401,] -24061.0828 -19590.6876 [402,] -61123.6847 -24061.0828 [403,] -21453.1774 -61123.6847 [404,] -49404.0111 -21453.1774 [405,] -48253.5544 -49404.0111 [406,] 4818.9750 -48253.5544 [407,] -1878.3993 4818.9750 [408,] 32358.0517 -1878.3993 [409,] -2858.9143 32358.0517 [410,] 15164.1944 -2858.9143 [411,] 72250.8327 15164.1944 [412,] -325.2517 72250.8327 [413,] -3285.2804 -325.2517 [414,] -1915.0474 -3285.2804 [415,] 3736.8753 -1915.0474 [416,] -3596.8006 3736.8753 [417,] 38305.7594 -3596.8006 [418,] 37612.1963 38305.7594 [419,] -2670.6788 37612.1963 [420,] -42142.3471 -2670.6788 [421,] -41220.1957 -42142.3471 [422,] -41415.3922 -41220.1957 [423,] -41942.1361 -41415.3922 [424,] -41601.6787 -41942.1361 [425,] -41443.5895 -41601.6787 [426,] -41299.7036 -41443.5895 [427,] -41479.2746 -41299.7036 [428,] -16823.2574 -41479.2746 [429,] -13978.7535 -16823.2574 [430,] -1711.4855 -13978.7535 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -34371.7678 -43371.7534 2 -40486.6638 -34371.7678 3 -37562.3114 -40486.6638 4 -37007.6557 -37562.3114 5 -44482.4940 -37007.6557 6 -43355.1770 -44482.4940 7 -41865.3726 -43355.1770 8 -41577.1678 -41865.3726 9 -41710.0026 -41577.1678 10 -37448.7338 -41710.0026 11 -36544.2726 -37448.7338 12 -40512.1430 -36544.2726 13 -24528.5488 -40512.1430 14 -25821.3059 -24528.5488 15 -27984.0194 -25821.3059 16 -26323.7630 -27984.0194 17 -26335.1276 -26323.7630 18 -26527.2446 -26335.1276 19 -26874.5029 -26527.2446 20 -27244.9860 -26874.5029 21 -46834.2625 -27244.9860 22 394424.2984 -46834.2625 23 317859.2305 394424.2984 24 254052.6223 317859.2305 25 260797.1407 254052.6223 26 372709.1556 260797.1407 27 255171.7832 372709.1556 28 31952.3948 255171.7832 29 17014.5964 31952.3948 30 19366.7834 17014.5964 31 -42834.7503 19366.7834 32 -28966.4211 -42834.7503 33 -40520.5888 -28966.4211 34 -42600.3811 -40520.5888 35 -40545.7029 -42600.3811 36 -46278.7468 -40545.7029 37 1037.7896 -46278.7468 38 -35880.6942 1037.7896 39 -29524.8634 -35880.6942 40 -42467.9516 -29524.8634 41 -35317.8389 -42467.9516 42 -44790.4854 -35317.8389 43 -27844.4979 -44790.4854 44 -27665.3906 -27844.4979 45 -27299.5662 -27665.3906 46 -44539.2310 -27299.5662 47 -42763.1335 -44539.2310 48 -34879.5544 -42763.1335 49 -29361.1868 -34879.5544 50 -43154.5830 -29361.1868 51 -44528.5759 -43154.5830 52 -31197.3421 -44528.5759 53 -919.5405 -31197.3421 54 -13142.3500 -919.5405 55 -43539.2198 -13142.3500 56 -48144.7547 -43539.2198 57 -44173.6566 -48144.7547 58 270077.6968 -44173.6566 59 43016.1017 270077.6968 60 -5216.6725 43016.1017 61 135400.1381 -5216.6725 62 -25124.3138 135400.1381 63 57554.8522 -25124.3138 64 54770.8361 57554.8522 65 37351.2850 54770.8361 66 54472.9698 37351.2850 67 19371.4437 54472.9698 68 -6325.0237 19371.4437 69 -44477.0637 -6325.0237 70 113682.6562 -44477.0637 71 -42867.0989 113682.6562 72 -23058.2195 -42867.0989 73 -25354.0345 -23058.2195 74 -24108.8591 -25354.0345 75 -15058.4672 -24108.8591 76 15745.4635 -15058.4672 77 96908.0479 15745.4635 78 -22103.7665 96908.0479 79 -27354.3574 -22103.7665 80 307864.9365 -27354.3574 81 255314.7557 307864.9365 82 19367.7117 255314.7557 83 -41405.8465 19367.7117 84 -31306.6462 -41405.8465 85 -44797.8953 -31306.6462 86 -39340.1469 -44797.8953 87 -44954.4951 -39340.1469 88 -3237.3400 -44954.4951 89 57272.3226 -3237.3400 90 -44674.6226 57272.3226 91 -24228.1029 -44674.6226 92 -24728.5238 -24228.1029 93 -27270.0839 -24728.5238 94 -45857.6396 -27270.0839 95 -45577.7990 -45857.6396 96 -47329.2162 -45577.7990 97 -39456.8467 -47329.2162 98 -41753.6741 -39456.8467 99 -45409.9825 -41753.6741 100 -39775.0284 -45409.9825 101 -52474.6781 -39775.0284 102 261954.5651 -52474.6781 103 1572683.1017 261954.5651 104 233148.7364 1572683.1017 105 257829.4866 233148.7364 106 295179.2298 257829.4866 107 144444.5782 295179.2298 108 255189.2310 144444.5782 109 28173.9039 255189.2310 110 19368.2609 28173.9039 111 -42805.4550 19368.2609 112 -36881.3796 -42805.4550 113 -23651.2955 -36881.3796 114 -26076.5600 -23651.2955 115 -32568.1605 -26076.5600 116 -46433.0715 -32568.1605 117 -38176.5027 -46433.0715 118 19667.3742 -38176.5027 119 -3215.1667 19667.3742 120 -21081.3276 -3215.1667 121 -27190.5768 -21081.3276 122 -46144.9875 -27190.5768 123 263947.5932 -46144.9875 124 257903.4131 263947.5932 125 274671.0640 257903.4131 126 255578.0824 274671.0640 127 19369.1730 255578.0824 128 -45089.1646 19369.1730 129 -25832.9709 -45089.1646 130 -28040.8845 -25832.9709 131 12265.4063 -28040.8845 132 -2121.8332 12265.4063 133 -42954.3339 -2121.8332 134 13462.8850 -42954.3339 135 -31495.6327 13462.8850 136 -14172.0842 -31495.6327 137 -21462.6264 -14172.0842 138 -50553.3393 -21462.6264 139 8573.1681 -50553.3393 140 -6845.0189 8573.1681 141 4526.0989 -6845.0189 142 -2057.2294 4526.0989 143 -42948.1585 -2057.2294 144 -16622.1098 -42948.1585 145 -13980.9615 -16622.1098 146 -20027.9624 -13980.9615 147 -25965.7807 -20027.9624 148 -21680.0988 -25965.7807 149 -17503.3726 -21680.0988 150 -21446.2447 -17503.3726 151 -45728.8460 -21446.2447 152 -45255.4285 -45728.8460 153 -50579.3918 -45255.4285 154 4613.4976 -50579.3918 155 1456.7309 4613.4976 156 -1671.0463 1456.7309 157 -28086.2722 -1671.0463 158 -41456.4567 -28086.2722 159 -14135.4520 -41456.4567 160 -5840.9618 -14135.4520 161 -14366.3933 -5840.9618 162 642.2519 -14366.3933 163 15156.7720 642.2519 164 -28218.5650 15156.7720 165 -24800.7800 -28218.5650 166 -10503.7408 -24800.7800 167 -6014.6166 -10503.7408 168 -14149.9627 -6014.6166 169 -19783.1928 -14149.9627 170 -21540.2681 -19783.1928 171 -21344.7612 -21540.2681 172 -45477.7571 -21344.7612 173 -47836.2339 -45477.7571 174 8757.8835 -47836.2339 175 -42946.9816 8757.8835 176 48916.7561 -42946.9816 177 -11086.9200 48916.7561 178 -11227.5401 -11086.9200 179 -28353.4309 -11227.5401 180 -14083.2630 -28353.4309 181 -8297.8835 -14083.2630 182 -21283.3830 -8297.8835 183 -21478.2873 -21283.3830 184 -49365.8792 -21478.2873 185 -49986.4445 -49365.8792 186 54960.0857 -49986.4445 187 -14325.5584 54960.0857 188 -1451.0215 -14325.5584 189 -43041.9516 -1451.0215 190 43475.0731 -43041.9516 191 -14095.3030 43475.0731 192 -21457.0397 -14095.3030 193 -51631.5519 -21457.0397 194 36954.5370 -51631.5519 195 20021.0151 36954.5370 196 -2423.2971 20021.0151 197 -36458.8656 -2423.2971 198 -42270.6324 -36458.8656 199 -17122.5471 -42270.6324 200 -39287.5223 -17122.5471 201 -41195.7133 -39287.5223 202 -41966.6224 -41195.7133 203 -9224.0604 -41966.6224 204 -39703.5078 -9224.0604 205 -41317.8069 -39703.5078 206 -41228.4986 -41317.8069 207 -11786.2626 -41228.4986 208 -1058.3495 -11786.2626 209 -11404.5677 -1058.3495 210 -14046.0324 -11404.5677 211 -23008.6023 -14046.0324 212 -22262.8924 -23008.6023 213 -22264.0542 -22262.8924 214 -10276.8540 -22264.0542 215 -21287.9087 -10276.8540 216 -33671.3010 -21287.9087 217 -51342.0042 -33671.3010 218 -40400.7039 -51342.0042 219 -31302.0189 -40400.7039 220 -51827.9467 -31302.0189 221 -43920.4483 -51827.9467 222 8858.5018 -43920.4483 223 8254.8441 8858.5018 224 -6193.0124 8254.8441 225 -2261.6813 -6193.0124 226 -42972.0650 -2261.6813 227 -9501.6653 -42972.0650 228 5328.4893 -9501.6653 229 -11897.7741 5328.4893 230 -14145.7387 -11897.7741 231 -20655.2260 -14145.7387 232 -29433.4309 -20655.2260 233 -13638.7970 -29433.4309 234 -21438.5047 -13638.7970 235 -50930.2093 -21438.5047 236 8762.6115 -50930.2093 237 -497.3053 8762.6115 238 -38773.8304 -497.3053 239 -41472.7734 -38773.8304 240 -14084.0471 -41472.7734 241 -21452.3119 -14084.0471 242 8766.2970 -21452.3119 243 -42952.6808 8766.2970 244 33249.1059 -42952.6808 245 -8610.9119 33249.1059 246 34640.0332 -8610.9119 247 3408.5755 34640.0332 248 -19117.5643 3408.5755 249 -8886.1021 -19117.5643 250 -24503.5809 -8886.1021 251 -14050.3900 -24503.5809 252 -22368.6266 -14050.3900 253 -21446.6350 -22368.6266 254 -49999.7310 -21446.6350 255 -39674.3499 -49999.7310 256 8756.4385 -39674.3499 257 -40754.4781 8756.4385 258 -41373.1437 -40754.4781 259 -39504.6980 -41373.1437 260 -43861.4218 -39504.6980 261 -41520.0628 -43861.4218 262 -42340.7744 -41520.0628 263 -42283.8891 -42340.7744 264 -42154.1613 -42283.8891 265 -41418.8906 -42154.1613 266 -41190.4944 -41418.8906 267 -41522.7957 -41190.4944 268 -16106.5271 -41522.7957 269 70213.6112 -16106.5271 270 -12965.7509 70213.6112 271 -14252.4176 -12965.7509 272 -9762.3188 -14252.4176 273 28326.6824 -9762.3188 274 -72197.8802 28326.6824 275 -21452.9153 -72197.8802 276 -50498.1878 -21452.9153 277 12404.7357 -50498.1878 278 -1523.0664 12404.7357 279 -1854.3021 -1523.0664 280 5584.5857 -1854.3021 281 -35511.0758 5584.5857 282 -41448.1311 -35511.0758 283 -14353.3936 -41448.1311 284 49896.4346 -14353.3936 285 -15933.3086 49896.4346 286 -16502.9445 -15933.3086 287 -15842.8317 -16502.9445 288 -16911.7657 -15842.8317 289 -4342.8956 -16911.7657 290 -8187.1311 -4342.8956 291 -21499.7920 -8187.1311 292 8760.4757 -21499.7920 293 16265.5438 8760.4757 294 -38678.7183 16265.5438 295 -40741.9988 -38678.7183 296 -41513.9005 -40741.9988 297 -39824.0098 -41513.9005 298 -39528.8470 -39824.0098 299 58580.3827 -39528.8470 300 -41856.3529 58580.3827 301 -41504.9705 -41856.3529 302 -22330.2434 -41504.9705 303 -6385.1142 -22330.2434 304 -12527.7453 -6385.1142 305 5355.5399 -12527.7453 306 -14283.6774 5355.5399 307 -9538.7591 -14283.6774 308 -21436.5161 -9538.7591 309 -43078.4997 -21436.5161 310 87610.1064 -43078.4997 311 -54810.7466 87610.1064 312 -50672.2440 -54810.7466 313 -52783.3964 -50672.2440 314 8755.3847 -52783.3964 315 -8680.0354 8755.3847 316 -37876.0514 -8680.0354 317 -41478.2398 -37876.0514 318 -7438.8236 -41478.2398 319 -6101.6914 -7438.8236 320 -14094.9711 -6101.6914 321 -21214.9683 -14094.9711 322 -48079.1703 -21214.9683 323 -51875.8077 -48079.1703 324 38832.0789 -51875.8077 325 -2844.9563 38832.0789 326 -39215.4444 -2844.9563 327 -41242.5948 -39215.4444 328 -48228.0335 -41242.5948 329 -41468.5709 -48228.0335 330 -15124.5938 -41468.5709 331 -14183.2363 -15124.5938 332 -21458.0408 -14183.2363 333 -48789.3135 -21458.0408 334 19713.6160 -48789.3135 335 -1903.7392 19713.6160 336 -42976.1144 -1903.7392 337 -14123.6980 -42976.1144 338 -21293.0822 -14123.6980 339 -43749.0207 -21293.0822 340 -48908.9387 -43749.0207 341 -52577.5135 -48908.9387 342 -51717.4612 -52577.5135 343 -49636.8821 -51717.4612 344 -50546.0696 -49636.8821 345 -37627.9805 -50546.0696 346 -51854.7547 -37627.9805 347 -49788.5143 -51854.7547 348 15840.8768 -49788.5143 349 51414.1747 15840.8768 350 -4374.0629 51414.1747 351 241.7956 -4374.0629 352 17621.4624 241.7956 353 181.3042 17621.4624 354 37760.0776 181.3042 355 27797.2372 37760.0776 356 35603.5300 27797.2372 357 -1544.5919 35603.5300 358 -990.1866 -1544.5919 359 -39439.9441 -990.1866 360 -40600.8492 -39439.9441 361 -38590.3090 -40600.8492 362 -40687.3700 -38590.3090 363 -41530.3466 -40687.3700 364 -41872.8288 -41530.3466 365 -40339.2588 -41872.8288 366 -41432.7408 -40339.2588 367 -14132.9224 -41432.7408 368 -21475.7677 -14132.9224 369 8754.7509 -21475.7677 370 -34028.1065 8754.7509 371 -41735.1039 -34028.1065 372 -41346.6164 -41735.1039 373 -14143.8862 -41346.6164 374 -21448.1433 -14143.8862 375 7729.7929 -21448.1433 376 40343.7547 7729.7929 377 46018.7664 40343.7547 378 25959.6806 46018.7664 379 -8497.9067 25959.6806 380 62712.9737 -8497.9067 381 -709.0604 62712.9737 382 -1835.9757 -709.0604 383 -1675.7406 -1835.9757 384 24880.5713 -1675.7406 385 37943.2659 24880.5713 386 66802.8471 37943.2659 387 -141.0171 66802.8471 388 -1818.1783 -141.0171 389 -42998.9352 -1818.1783 390 -25722.1518 -42998.9352 391 -9579.2388 -25722.1518 392 -13423.9607 -9579.2388 393 14948.5575 -13423.9607 394 -21352.8312 14948.5575 395 104425.8356 -21352.8312 396 -14020.8518 104425.8356 397 -18119.8950 -14020.8518 398 -18812.5436 -18119.8950 399 168783.5740 -18812.5436 400 -19590.6876 168783.5740 401 -24061.0828 -19590.6876 402 -61123.6847 -24061.0828 403 -21453.1774 -61123.6847 404 -49404.0111 -21453.1774 405 -48253.5544 -49404.0111 406 4818.9750 -48253.5544 407 -1878.3993 4818.9750 408 32358.0517 -1878.3993 409 -2858.9143 32358.0517 410 15164.1944 -2858.9143 411 72250.8327 15164.1944 412 -325.2517 72250.8327 413 -3285.2804 -325.2517 414 -1915.0474 -3285.2804 415 3736.8753 -1915.0474 416 -3596.8006 3736.8753 417 38305.7594 -3596.8006 418 37612.1963 38305.7594 419 -2670.6788 37612.1963 420 -42142.3471 -2670.6788 421 -41220.1957 -42142.3471 422 -41415.3922 -41220.1957 423 -41942.1361 -41415.3922 424 -41601.6787 -41942.1361 425 -41443.5895 -41601.6787 426 -41299.7036 -41443.5895 427 -41479.2746 -41299.7036 428 -16823.2574 -41479.2746 429 -13978.7535 -16823.2574 430 -1711.4855 -13978.7535 > 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/html/rcomp/tmp/7mcj81291380013.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/html/rcomp/tmp/8mcj81291380013.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/html/rcomp/tmp/9mcj81291380013.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/html/rcomp/tmp/10xm0t1291380013.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/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/html/rcomp/tmp/110mhz1291380013.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/html/rcomp/tmp/123nxn1291380013.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/html/rcomp/tmp/13zedd1291380013.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/html/rcomp/tmp/14d7we1291380014.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/html/rcomp/tmp/15y7u21291380014.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/html/rcomp/tmp/16dhsb1291380014.tab") + } > > try(system("convert tmp/18klz1291380013.ps tmp/18klz1291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/28klz1291380013.ps tmp/28klz1291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/3ic221291380013.ps tmp/3ic221291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/4ic221291380013.ps tmp/4ic221291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/5ic221291380013.ps tmp/5ic221291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/6t3251291380013.ps tmp/6t3251291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/7mcj81291380013.ps tmp/7mcj81291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/8mcj81291380013.ps tmp/8mcj81291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/9mcj81291380013.ps tmp/9mcj81291380013.png",intern=TRUE)) character(0) > try(system("convert tmp/10xm0t1291380013.ps tmp/10xm0t1291380013.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 15.484 2.375 32.565