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 + ,1 + ,1595 + ,17 + ,60720 + ,319369 + ,143 + ,7022 + ,75 + ,0 + ,61 + ,0 + ,5565 + ,47 + ,94355 + ,493408 + ,38 + ,6243 + ,53 + ,39 + ,287 + ,0 + ,601 + ,6 + ,60720 + ,319210 + ,144 + ,19868 + ,198 + ,0 + ,268 + ,1 + ,188 + ,11 + ,77655 + ,381180 + ,67 + ,16471 + ,964 + ,0 + ,25 + ,1 + ,7146 + ,105 + ,134028 + ,297978 + ,367 + ,933 + ,14 + ,305 + ,418 + ,0 + ,1135 + ,17 + ,62285 + ,290476 + ,384 + ,5322 + ,80 + ,-9 + ,392 + ,1 + ,450 + ,8 + ,59325 + ,292136 + ,380 + ,11517 + ,205 + ,-1 + ,131 + ,1 + ,34 + ,8 + ,60630 + ,314353 + ,309 + ,14294 + ,3340 + ,0 + ,316 + ,1 + ,133 + ,14 + ,65990 + ,339445 + ,109 + ,9960 + ,1052 + ,0 + ,347 + ,1 + ,119 + ,6 + ,59118 + ,303677 + ,354 + ,17280 + ,872 + ,0 + ,68 + ,0 + ,2053 + ,53 + ,100423 + ,397144 + ,60 + ,3720 + ,96 + ,33 + ,70 + ,0 + ,4036 + ,63 + ,100269 + ,424898 + ,53 + ,3570 + ,56 + ,32 + ,147 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,0 + ,169 + ,1 + ,4655 + ,393 + ,60720 + ,341570 + ,105 + ,360 + ,30 + ,0 + ,248 + ,1 + ,131 + ,11 + ,61808 + ,308989 + ,338 + ,9908 + ,834 + ,0 + ,152 + ,1 + ,1766 + ,73 + ,60438 + ,305959 + ,349 + ,1451 + ,60 + ,0 + ,351 + ,1 + ,312 + ,14 + ,58598 + ,318690 + ,147 + ,8478 + ,381 + ,0 + ,372 + ,1 + ,448 + ,40 + ,59781 + ,323361 + ,132 + ,3084 + ,275 + ,0 + ,137 + ,1 + ,115 + ,13 + ,60945 + ,318903 + ,145 + ,9146 + ,1037 + ,0 + ,352 + ,1 + ,60 + ,10 + ,61124 + ,314049 + ,314 + ,11405 + ,1916 + ,0 + ,338 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,180 + ,0 + ,258 + ,1 + ,364 + ,35 + ,47705 + ,298466 + ,365 + ,2813 + ,271 + ,0 + ,343 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,183 + ,0 + ,398 + ,0 + ,1442 + ,118 + ,121173 + ,438493 + ,49 + ,2021 + ,165 + ,-3 + ,115 + ,0 + ,1389 + ,9 + ,57530 + ,378049 + ,71 + ,19783 + ,128 + ,6 + ,188 + ,1 + ,149 + ,5 + ,62041 + ,313332 + ,318 + ,22666 + ,761 + ,0 + ,379 + ,1 + ,2212 + ,14 + ,60720 + ,319864 + ,142 + ,8562 + ,54 + ,0 + ,48 + ,0 + ,7489 + ,322 + ,136996 + ,430866 + ,50 + ,717 + ,31 + ,70 + ,366 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,287 + ,0 + ,117 + ,1 + ,402 + ,26 + ,84990 + ,332743 + ,115 + ,5106 + ,330 + ,5 + ,412 + ,0 + ,7419 + ,29 + ,63255 + ,286963 + ,387 + ,2999 + ,12 + ,-5 + ,157 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,219 + ,0 + ,264 + ,1 + ,307 + ,15 + ,58856 + ,331955 + ,116 + ,8797 + ,430 + ,0 + ,54 + ,0 + ,1134 + ,162 + ,146216 + ,527021 + ,34 + ,2019 + ,288 + ,58 + ,64 + ,0 + ,7561 + ,87 + ,82425 + ,364304 + ,86 + ,1889 + ,22 + ,37 + ,113 + ,1 + ,5131 + ,71 + ,86111 + ,292154 + ,379 + ,1298 + ,18 + ,9 + ,148 + ,1 + ,9 + ,6 + ,60835 + ,314987 + ,300 + ,19165 + ,12458 + ,0 + ,425 + ,1 + ,2264 + ,24 + ,55174 + ,272713 + ,394 + ,3030 + ,32 + ,-31 + ,7 + ,1 + ,40949 + ,481 + ,129352 + ,1398893 + ,9 + ,2493 + ,29 + ,2865 + ,197 + ,0 + ,4980 + ,91 + ,113521 + ,409642 + ,55 + ,2304 + ,42 + ,0 + ,97 + ,1 + ,272 + ,8 + ,112995 + ,429112 + ,51 + ,28639 + ,844 + ,14 + ,168 + ,1 + ,757 + ,19 + ,45689 + ,382712 + ,66 + ,9616 + ,242 + ,0 + ,69 + ,1 + ,1424 + ,94 + ,131116 + ,374943 + ,77 + ,1861 + ,123 + ,32 + ,127 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,191 + ,0 + ,184 + ,0 + ,378 + ,52 + ,71873 + ,297413 + ,369 + ,1873 + ,258 + ,0 + ,401 + ,1 + ,1061 + ,5 + ,70415 + ,304252 + ,353 + ,20850 + ,98 + ,-4 + ,371 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,290 + ,0 + ,174 + ,1 + ,203 + ,9 + ,61938 + ,309596 + ,334 + ,12177 + ,541 + ,0 + ,46 + ,0 + ,3689 + ,25 + ,101792 + ,377305 + ,74 + ,7092 + ,48 + ,74 + ,90 + ,1 + ,567 + ,11 + ,103345 + ,386212 + ,65 + ,16928 + ,329 + ,20 + ,27 + ,0 + ,27664 + ,282 + ,74996 + ,657954 + ,24 + ,1624 + ,17 + ,247 + ,105 + ,1 + ,1686 + ,40 + ,68696 + ,368186 + ,81 + ,4205 + ,100 + ,11 + ,413 + ,0 + ,1164 + ,17 + ,57320 + ,269753 + ,395 + ,4103 + ,60 + ,-6 + ,26 + ,1 + ,15824 + ,64 + ,63838 + ,585715 + ,30 + ,6027 + ,24 + ,284 + ,16 + ,1 + ,6575 + ,159 + ,151352 + ,1071292 + ,11 + ,5480 + ,133 + ,819 + ,19 + ,0 + ,37597 + ,248 + ,37884 + ,992426 + ,12 + ,3195 + ,21 + ,652 + ,336 + ,1 + ,314 + ,5 + ,60720 + ,315371 + ,294 + ,23074 + ,367 + ,0 + ,110 + ,0 + ,1932 + ,36 + ,72954 + ,249898 + ,403 + ,1386 + ,26 + ,10 + ,86 + ,0 + ,5692 + ,72 + ,64239 + ,357312 + ,91 + ,2185 + ,28 + ,22 + ,294 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,257 + ,0 + ,322 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,177 + ,0 + ,94 + ,0 + ,6533 + ,75 + ,126910 + ,364839 + ,85 + ,2198 + ,25 + ,16 + ,400 + ,0 + ,2126 + ,138 + ,21001 + ,230621 + ,413 + ,222 + ,14 + ,-3 + ,150 + ,1 + ,14 + ,1 + ,61016 + ,315877 + ,161 + ,115877 + ,8041 + ,0 + ,107 + ,0 + ,3192 + ,60 + ,59831 + ,125390 + ,427 + ,-1243 + ,-23 + ,10 + ,166 + ,1 + ,47 + ,3 + ,60720 + ,314882 + ,305 + ,38294 + ,2432 + ,0 + ,47 + ,0 + ,4891 + ,217 + ,210568 + ,370837 + ,79 + ,787 + ,35 + ,73 + ,182 + ,0 + ,531 + ,10 + ,71641 + ,330068 + ,120 + ,13007 + ,245 + ,0 + ,245 + ,1 + ,2332 + ,4 + ,72680 + ,324385 + ,131 + ,31096 + ,53 + ,0 + ,330 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,277 + ,0 + ,10 + ,0 + ,16275 + ,109 + ,47818 + ,1443586 + ,7 + ,11409 + ,76 + ,1974 + ,429 + ,1 + ,9252 + ,48 + ,79510 + ,253537 + ,401 + ,1115 + ,6 + ,-129 + ,2 + ,0 + ,29790 + ,309 + ,81767 + ,4321023 + ,2 + ,13337 + ,138 + ,28328 + ,235 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,204 + ,0 + ,187 + ,0 + ,2300 + ,29 + ,64873 + ,352108 + ,97 + ,5245 + ,66 + ,0 + ,213 + ,1 + ,267 + ,32 + ,57640 + ,330059 + ,121 + ,4064 + ,487 + ,0 + ,296 + ,0 + ,382 + ,21 + ,59365 + ,340968 + ,107 + ,6713 + ,369 + ,0 + ,218 + ,1 + ,397 + ,28 + ,60646 + ,317736 + ,152 + ,4205 + ,297 + ,0 + ,28 + ,1 + ,11920 + ,182 + ,31080 + ,209458 + ,419 + ,52 + ,1 + ,238 + ,5 + ,0 + ,54660 + ,302 + ,126942 + ,1491348 + ,6 + ,4276 + ,24 + ,3040 + ,138 + ,1 + ,5 + ,2 + ,60720 + ,314887 + ,304 + ,57443 + ,24814 + ,0 + ,163 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,0 + ,219 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,199 + ,0 + ,125 + ,0 + ,2647 + ,37 + ,62920 + ,353058 + ,95 + ,4137 + ,58 + ,1 + ,216 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,201 + ,0 + ,196 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,263 + ,0 + ,128 + ,1 + ,94 + ,7 + ,60793 + ,314533 + ,308 + ,16362 + ,1212 + ,0 + ,223 + ,1 + ,3 + ,2 + ,60698 + ,315354 + ,296 + ,57677 + ,36505 + ,0 + ,165 + ,1 + ,454 + ,51 + ,63261 + ,302187 + ,357 + ,2004 + ,225 + ,0 + ,373 + ,1 + ,227 + ,9 + ,76644 + ,336639 + ,110 + ,15182 + ,603 + ,0 + ,324 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,245 + ,0 + ,251 + ,1 + ,206 + ,9 + ,53110 + ,296702 + ,372 + ,10745 + ,469 + ,0 + ,11 + ,1 + ,25830 + ,115 + ,245546 + ,1073089 + ,10 + ,7592 + ,34 + ,1910 + ,62 + ,1 + ,1528 + ,38 + ,60326 + ,146494 + ,426 + ,-1408 + ,-35 + ,38 + ,149 + ,1 + ,271 + ,2 + ,60720 + ,325249 + ,129 + ,62625 + ,463 + ,0 + ,270 + ,1 + ,138 + ,6 + ,69817 + ,331420 + ,117 + ,21903 + ,954 + ,0 + ,171 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,211 + ,0 + ,226 + ,1 + ,278 + ,26 + ,61404 + ,314922 + ,302 + ,4420 + ,413 + ,0 + ,276 + ,1 + ,282 + ,30 + ,62452 + ,320016 + ,141 + ,4001 + ,425 + ,0 + ,383 + ,1 + ,571 + ,80 + ,41477 + ,280398 + ,390 + ,1005 + ,141 + ,0 + ,93 + ,1 + ,2253 + ,99 + ,63593 + ,452469 + ,43 + ,2550 + ,112 + ,17 + ,302 + ,1 + ,290 + ,2 + ,58790 + ,301164 + ,360 + ,50582 + ,348 + ,0 + ,385 + ,1 + ,78 + ,16 + ,62700 + ,317330 + ,154 + ,7333 + ,1511 + ,0 + ,225 + ,1 + ,20 + ,2 + ,60805 + ,315576 + ,169 + ,57788 + ,5924 + ,0 + ,275 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,284 + ,0 + ,431 + ,1 + ,18073 + ,180 + ,27284 + ,-7170 + ,429 + ,-1151 + ,-11 + ,-366 + ,220 + ,1 + ,1866 + ,29 + ,56225 + ,322331 + ,135 + ,4218 + ,66 + ,0 + ,6 + ,0 + ,42634 + ,156 + ,157214 + ,1629616 + ,5 + ,9164 + ,34 + ,2964 + ,290 + ,1 + ,249 + ,62 + ,54323 + ,292754 + ,378 + ,1496 + ,373 + ,0 + ,129 + ,1 + ,422 + ,28 + ,57935 + ,318056 + ,150 + ,4216 + ,280 + ,0 + ,65 + ,1 + ,2675 + ,79 + ,59017 + ,355178 + ,93 + ,1964 + ,58 + ,36 + ,416 + ,1 + ,965 + ,3 + ,73490 + ,204325 + ,421 + ,1442 + ,4 + ,-8 + ,299 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,224 + ,0 + ,407 + ,1 + ,621 + ,16 + ,68005 + ,317046 + ,156 + ,7315 + ,188 + ,-5 + ,252 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,240 + ,0 + ,257 + ,1 + ,365 + ,8 + ,54820 + ,309560 + ,335 + ,13695 + ,300 + ,0 + ,63 + ,0 + ,3122 + ,34 + ,94670 + ,414462 + ,54 + ,6308 + ,69 + ,38 + ,21 + ,0 + ,12988 + ,119 + ,82340 + ,857217 + ,16 + ,5523 + ,51 + ,601 + ,38 + ,0 + ,5336 + ,81 + ,112477 + ,697458 + ,22 + ,6141 + ,93 + ,132 + ,51 + ,0 + ,3160 + ,108 + ,108094 + ,530670 + ,33 + ,3062 + ,105 + ,66 + ,386 + ,1 + ,2489 + ,40 + ,67804 + ,238125 + ,408 + ,953 + ,15 + ,0 + ,325 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,273 + ,0 + ,33 + ,0 + ,6984 + ,45 + ,80570 + ,741409 + ,19 + ,12031 + ,78 + ,195 + ,91 + ,1 + ,2001 + ,52 + ,95551 + ,393343 + ,62 + ,3718 + ,97 + ,19 + ,430 + ,1 + ,15236 + ,90 + ,77440 + ,372631 + ,78 + ,1918 + ,11 + ,-158 + ,259 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,239 + ,0 + ,382 + ,1 + ,530 + ,60 + ,73433 + ,317291 + ,155 + ,1955 + ,221 + ,0 + ,369 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,292 + ,0 + ,428 + ,1 + ,35624 + ,223 + ,157278 + ,306275 + ,347 + ,477 + ,3 + ,-70 + ,100 + ,0 + ,1383 + ,11 + ,73221 + ,317892 + ,151 + ,10717 + ,85 + ,13 + ,122 + ,1 + ,875 + ,14 + ,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 + ,1 + ,72 + ,3 + ,60398 + ,314210 + ,312 + ,38070 + ,1591 + ,0 + ,210 + ,0 + ,265 + ,31 + ,60720 + ,306948 + ,346 + ,3450 + ,404 + ,194 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,186 + ,0 + ,277 + ,1 + ,335 + ,33 + ,64175 + ,320398 + ,140 + ,3648 + ,359 + ,0 + ,298 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,230 + ,0 + ,356 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,269 + ,0 + ,43 + ,1 + ,4525 + ,226 + ,93811 + ,501749 + ,37 + ,1335 + ,67 + ,101 + ,71 + ,0 + ,3045 + ,58 + ,27330 + ,202055 + ,422 + ,35 + ,1 + ,31 + ,133 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,196 + ,0 + ,212 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,0 + ,186 + ,0 + ,638 + ,14 + ,60370 + ,333210 + ,114 + ,9515 + ,209 + ,0 + ,262 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,283 + ,0 + ,378 + ,1 + ,607 + ,11 + ,60436 + ,322340 + ,134 + ,11122 + ,201 + ,0 + ,132 + ,1 + ,1558 + ,63 + ,55637 + ,369448 + ,80 + ,2690 + ,109 + ,0 + ,279 + ,0 + ,1324 + ,117 + ,67440 + ,291841 + ,381 + ,785 + ,69 + ,0 + ,376 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,280 + ,0 + ,267 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,0 + ,403 + ,1 + ,611 + ,8 + ,59190 + ,296919 + ,371 + ,12115 + ,159 + ,-4 + ,303 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,225 + ,0 + ,388 + ,1 + ,923 + ,8 + ,58620 + ,309038 + ,337 + ,13630 + ,118 + ,0 + ,420 + ,1 + ,661 + ,3 + ,65920 + ,246541 + ,405 + ,15514 + ,70 + ,-12 + ,114 + ,0 + ,2397 + ,14 + ,74020 + ,289513 + ,385 + ,6394 + ,37 + ,6 + ,305 + ,1 + ,366 + ,53 + ,61808 + ,344425 + ,101 + ,2725 + ,394 + ,0 + ,346 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,181 + ,0 + ,354 + ,1 + ,135 + ,2 + ,62065 + ,314210 + ,311 + ,57105 + ,847 + ,0 + ,124 + ,0 + ,1659 + ,29 + ,107577 + ,480382 + ,40 + ,9668 + ,169 + ,2 + ,136 + ,1 + ,316 + ,9 + ,60505 + ,315009 + ,299 + ,12779 + ,364 + ,0 + ,130 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,197 + ,0 + ,146 + ,1 + ,309 + ,11 + ,56535 + ,312878 + ,322 + ,10262 + ,365 + ,0 + ,314 + ,0 + ,49 + ,8 + ,64107 + ,322031 + ,137 + ,15254 + ,2500 + ,0 + ,241 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,200 + ,0 + ,34 + ,0 + ,4519 + ,49 + ,102129 + ,597793 + ,29 + ,8118 + ,88 + ,163 + ,384 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,0 + ,120 + ,1 + ,837 + ,69 + ,61262 + ,315688 + ,165 + ,1677 + ,138 + ,4 + ,72 + ,1 + ,5119 + ,49 + ,39039 + ,378525 + ,69 + ,3643 + ,35 + ,31 + ,79 + ,1 + ,1280 + ,117 + ,69465 + ,312378 + ,326 + ,960 + ,88 + ,26 + ,99 + ,1 + ,2564 + ,22 + ,130140 + ,403560 + ,59 + ,9253 + ,79 + ,13 + ,108 + ,0 + ,2045 + ,20 + ,97890 + ,510834 + ,36 + ,15542 + ,152 + ,10 + ,119 + ,0 + ,2234 + ,16 + ,77200 + ,214215 + ,418 + ,888 + ,6 + ,5 + ,102 + ,1 + ,975 + ,23 + ,90534 + ,235133 + ,410 + ,1528 + ,36 + ,13 + ,209 + ,0 + ,1136 + ,32 + ,48522 + ,343613 + ,103 + ,4488 + ,126 + ,0 + ,293 + ,1 + ,453 + ,21 + ,81125 + ,365959 + ,84 + ,7903 + ,366 + ,0 + ,211 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,210 + ,0 + ,227 + ,1 + ,61 + ,14 + ,60720 + ,314551 + ,307 + ,8182 + ,1874 + ,0 + ,234 + ,1 + ,368 + ,17 + ,61686 + ,303230 + ,355 + ,6072 + ,281 + ,0 + ,247 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,233 + ,0 + ,49 + ,1 + ,4901 + ,84 + ,121920 + ,469107 + ,41 + ,3204 + ,55 + ,67 + ,88 + ,1 + ,540 + ,14 + ,103960 + ,354228 + ,94 + ,11016 + ,286 + ,22 + ,329 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,253 + ,0 + ,160 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,213 + ,0 + ,271 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,286 + ,0 + ,306 + ,1 + ,2 + ,9 + ,60735 + ,315394 + ,174 + ,12822 + ,63057 + ,0 + ,172 + ,1 + ,36 + ,4 + ,61564 + ,312412 + ,325 + ,28103 + ,3146 + ,0 + ,402 + ,1 + ,776 + ,9 + ,64230 + ,333505 + ,113 + ,14834 + ,172 + ,-4 + ,4 + ,1 + ,84738 + ,588 + ,-26007 + ,223193 + ,415 + ,39 + ,0 + ,4145 + ,156 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,217 + ,0 + ,315 + ,1 + ,3 + ,4 + ,60761 + ,315656 + ,167 + ,28914 + ,37920 + ,0 + ,399 + ,1 + ,529 + ,9 + ,63870 + ,296261 + ,374 + ,10696 + ,182 + ,-3 + ,381 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,264 + ,0 + ,159 + ,1 + ,405 + ,7 + ,60845 + ,336425 + ,111 + ,19489 + ,337 + ,0 + ,181 + ,1 + ,972 + ,34 + ,71642 + ,359335 + ,87 + ,4686 + ,164 + ,0 + ,278 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,232 + ,0 + ,89 + ,1 + ,2099 + ,64 + ,106611 + ,308636 + ,339 + ,1697 + ,52 + ,20 + ,112 + ,0 + ,3437 + ,47 + ,48022 + ,158492 + ,425 + ,-883 + ,-12 + ,9 + ,201 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,243 + ,0 + ,344 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,182 + ,0 + ,22 + ,0 + ,22330 + ,84 + ,79801 + ,711969 + ,20 + ,6095 + ,23 + ,446 + ,291 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,258 + ,0 + ,269 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,221 + ,0 + ,164 + ,1 + ,483 + ,27 + ,60830 + ,306268 + ,348 + ,3936 + ,220 + ,0 + ,358 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,279 + ,0 + ,118 + ,1 + ,2239 + ,21 + ,88590 + ,442882 + ,47 + ,11566 + ,108 + ,5 + ,92 + ,0 + ,2949 + ,41 + ,82903 + ,378509 + ,70 + ,4354 + ,61 + ,18 + ,300 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,223 + ,0 + ,312 + ,1 + ,365 + ,83 + ,87192 + ,346611 + ,99 + ,1766 + ,402 + ,0 + ,183 + ,1 + ,2461 + ,57 + ,55792 + ,314289 + ,310 + ,2005 + ,46 + ,0 + ,18 + ,0 + ,21950 + ,519 + ,114337 + ,856956 + ,17 + ,1266 + ,30 + ,745 + ,126 + ,0 + ,3294 + ,23 + ,75832 + ,217193 + ,416 + ,748 + ,5 + ,1 + ,193 + ,1 + ,141 + ,6 + ,61630 + ,315366 + ,295 + ,19228 + ,819 + ,0 + ,390 + ,1 + ,572 + ,16 + ,58580 + ,307930 + ,342 + ,6746 + ,189 + ,-1 + ,23 + ,1 + ,13326 + ,102 + ,165548 + ,702380 + ,21 + ,4925 + ,38 + ,436 + ,50 + ,1 + ,2284 + ,33 + ,76403 + ,194493 + ,423 + ,-167 + ,-2 + ,66 + ,206 + ,0 + ,10 + ,2 + ,61656 + ,316155 + ,160 + ,58077 + ,11639 + ,0 + ,162 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,207 + ,0 + ,98 + ,0 + ,1414 + ,198 + ,70184 + ,330546 + ,119 + ,659 + ,92 + ,13 + ,75 + ,1 + ,1975 + ,35 + ,118881 + ,394510 + ,61 + ,5557 + ,98 + ,30 + ,134 + ,1 + ,43 + ,3 + ,60887 + ,312846 + ,323 + ,37615 + ,2613 + ,0 + ,180 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,190 + ,0 + ,337 + ,1 + ,844 + ,78 + ,60925 + ,296139 + ,375 + ,1233 + ,114 + ,0 + ,341 + ,1 + ,304 + ,15 + ,62969 + ,295580 + ,376 + ,6372 + ,314 + ,0 + ,154 + ,1 + ,458 + ,11 + ,58625 + ,297765 + ,368 + ,8888 + ,213 + ,0 + ,36 + ,1 + ,18562 + ,155 + ,102313 + ,377934 + ,72 + ,1148 + ,10 + ,143 + ,192 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,179 + ,0 + ,17 + ,1 + ,7123 + ,109 + ,288170 + ,638830 + ,28 + ,4026 + ,62 + ,769 + ,254 + ,1 + ,622 + ,73 + ,73007 + ,304376 + ,352 + ,1430 + ,168 + ,0 + ,374 + ,1 + ,174 + ,1 + ,64820 + ,307424 + ,344 + ,107424 + ,618 + ,0 + ,35 + ,1 + ,2220 + ,22 + ,301670 + ,644190 + ,26 + ,20190 + ,200 + ,154 + ,377 + ,1 + ,121 + ,24 + ,56178 + ,295370 + ,377 + ,3974 + ,788 + ,0 + ,32 + ,0 + ,11819 + ,85 + ,106113 + ,574339 + ,31 + ,4404 + ,32 + ,199 + ,260 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,235 + ,0 + ,274 + ,1 + ,125 + ,12 + ,60798 + ,310201 + ,332 + ,9183 + ,883 + ,0 + ,406 + ,1 + ,1182 + ,13 + ,70694 + ,327007 + ,124 + ,9770 + ,107 + ,-5 + ,116 + ,1 + ,1503 + ,39 + ,56364 + ,343466 + ,104 + ,3679 + ,95 + ,6 + ,320 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,274 + ,0 + ,367 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,248 + ,0 + ,288 + ,1 + ,30 + ,4 + ,62045 + ,318098 + ,149 + ,29525 + ,3912 + ,0 + ,414 + ,0 + ,3310 + ,33 + ,75230 + ,448243 + ,44 + ,7523 + ,75 + ,-6 + ,74 + ,1 + ,554 + ,11 + ,79285 + ,325738 + ,125 + ,11431 + ,227 + ,30 + ,250 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,234 + ,0 + ,224 + ,1 + ,468 + ,21 + ,60720 + ,312161 + ,328 + ,5341 + ,239 + ,0 + ,397 + ,1 + ,4917 + ,42 + ,52811 + ,243650 + ,406 + ,1039 + ,9 + ,-3 + ,417 + ,1 + ,3256 + ,66 + ,35250 + ,407159 + ,57 + ,3139 + ,64 + ,-8 + ,323 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,185 + ,0 + ,340 + ,1 + ,125 + ,13 + ,59734 + ,317698 + ,153 + ,9054 + ,944 + ,0 + ,239 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,195 + ,0 + ,326 + ,1 + ,22 + ,2 + ,60722 + ,312502 + ,324 + ,56251 + ,5154 + ,0 + ,261 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,262 + ,0 + ,357 + ,1 + ,514 + ,52 + ,78780 + ,322378 + ,133 + ,2353 + ,238 + ,0 + ,313 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,0 + ,360 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,270 + ,0 + ,272 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,285 + ,0 + ,333 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,276 + ,0 + ,355 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,272 + ,0 + ,362 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,268 + ,0 + ,37 + ,1 + ,10327 + ,91 + ,88577 + ,640273 + ,27 + ,4838 + ,43 + ,135 + ,42 + ,1 + ,13718 + ,61 + ,96448 + ,345783 + ,100 + ,2390 + ,11 + ,104 + ,57 + ,1 + ,3748 + ,9 + ,50350 + ,652925 + ,25 + ,50325 + ,121 + ,53 + ,30 + ,1 + ,14416 + ,361 + ,49857 + ,439798 + ,48 + ,664 + ,17 + ,220 + ,404 + ,0 + ,1526 + ,25 + ,69351 + ,278990 + ,391 + ,3160 + ,52 + ,-4 + ,80 + ,0 + ,666 + ,171 + ,117869 + ,339836 + ,108 + ,818 + ,210 + ,25 + ,422 + ,1 + ,2844 + ,67 + ,72683 + ,240897 + ,407 + ,610 + ,14 + ,-15 + ,297 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,229 + ,0 + ,318 + ,1 + ,368 + ,17 + ,61167 + ,297141 + ,370 + ,5714 + ,264 + ,0 + ,175 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,205 + ,0 + ,285 + ,1 + ,333 + ,13 + ,70811 + ,331323 + ,118 + ,10102 + ,394 + ,0 + ,339 + ,1 + ,26 + ,1 + ,60896 + ,313880 + ,315 + ,113880 + ,4328 + ,0 + ,243 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,237 + ,0 + ,167 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,193 + ,0 + ,411 + ,1 + ,1303 + ,61 + ,69863 + ,309422 + ,336 + ,1794 + ,84 + ,-5 + ,155 + ,1 + ,20 + ,3 + ,60938 + ,315245 + ,297 + ,38415 + ,5674 + ,0 + ,56 + ,1 + ,2384 + ,97 + ,61348 + ,405972 + ,58 + ,2123 + ,86 + ,54 + ,204 + ,0 + ,203 + ,27 + ,50804 + ,300962 + ,361 + ,3739 + ,498 + ,0 + ,178 + ,0 + ,71 + ,1 + ,60745 + ,316176 + ,159 + ,116176 + ,1646 + ,0 + ,295 + ,1 + ,53 + ,14 + ,59506 + ,302409 + ,356 + ,7315 + ,1932 + ,0 + ,396 + ,1 + ,562 + ,13 + ,58456 + ,283587 + ,389 + ,6430 + ,149 + ,-3 + ,393 + ,1 + ,622 + ,7 + ,60950 + ,263276 + ,398 + ,9039 + ,102 + ,-2 + ,135 + ,1 + ,645 + ,5 + ,60720 + ,312075 + ,329 + ,22415 + ,174 + ,0 + ,364 + ,0 + ,1763 + ,24 + ,61600 + ,308336 + ,340 + ,4514 + ,61 + ,0 + ,179 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,0 + ,231 + ,1 + ,317 + ,9 + ,63915 + ,298700 + ,363 + ,10967 + ,311 + ,0 + ,319 + ,1 + ,1 + ,6 + ,60719 + ,315372 + ,293 + ,19229 + ,84213 + ,0 + ,191 + ,1 + ,275 + ,4 + ,59500 + ,318745 + ,146 + ,29686 + ,431 + ,0 + ,321 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,246 + ,0 + ,111 + ,1 + ,936 + ,30 + ,67939 + ,408881 + ,56 + ,6963 + ,223 + ,9 + ,29 + ,1 + ,8568 + ,37 + ,32168 + ,786690 + ,18 + ,15856 + ,68 + ,223 + ,423 + ,0 + ,11528 + ,70 + ,-14545 + ,-83265 + ,431 + ,-4047 + ,-25 + ,-25 + ,202 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,250 + ,0 + ,232 + ,1 + ,738 + ,10 + ,60720 + ,321376 + ,139 + ,12138 + ,164 + ,0 + ,266 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,281 + ,0 + ,405 + ,1 + ,592 + ,7 + ,64270 + ,276898 + ,392 + ,10985 + ,130 + ,-4 + ,176 + ,1 + ,126 + ,10 + ,60951 + ,315547 + ,171 + ,11555 + ,921 + ,0 + ,236 + ,1 + ,2 + ,1 + ,60743 + ,315487 + ,172 + ,115487 + ,52257 + ,0 + ,142 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,220 + ,0 + ,13 + ,0 + ,18014 + ,239 + ,-1710 + ,1405225 + ,8 + ,5043 + ,67 + ,943 + ,281 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,256 + ,0 + ,328 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,0 + ,9 + ,0 + ,37704 + ,452 + ,60448 + ,983660 + ,14 + ,1734 + ,21 + ,2167 + ,207 + ,0 + ,63 + ,4 + ,65688 + ,318574 + ,148 + ,29643 + ,1891 + ,0 + ,73 + ,1 + ,1431 + ,34 + ,106885 + ,310768 + ,330 + ,3258 + ,77 + ,30 + ,242 + ,1 + ,94 + ,3 + ,61360 + ,312887 + ,321 + ,37629 + ,1200 + ,0 + ,229 + ,1 + ,192 + ,4 + ,65276 + ,312339 + ,327 + ,28085 + ,586 + ,0 + ,317 + ,1 + ,32 + ,1 + ,59988 + ,314964 + ,301 + ,114964 + ,3603 + ,0 + ,44 + ,1 + ,6869 + ,105 + ,117520 + ,379983 + ,68 + ,1714 + ,26 + ,90 + ,349 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,176 + ,0 + ,140 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,0 + ,361 + ,1 + ,1 + ,4 + ,60722 + ,315398 + ,173 + ,28849 + ,172235 + ,0 + ,203 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,247 + ,0 + ,255 + ,1 + ,2328 + ,94 + ,82732 + ,253588 + ,400 + ,570 + ,23 + ,0 + ,308 + ,0 + ,209 + ,13 + ,64016 + ,316647 + ,157 + ,8973 + ,558 + ,0 + ,195 + ,1 + ,28 + ,2 + ,60890 + ,315688 + ,166 + ,57844 + ,4148 + ,0 + ,280 + ,1 + ,176 + ,28 + ,68136 + ,310670 + ,331 + ,3952 + ,629 + ,0 + ,415 + ,1 + ,1920 + ,15 + ,79420 + ,165404 + ,424 + ,-2306 + ,-18 + ,-6 + ,3 + ,0 + ,87550 + ,458 + ,153198 + ,4111912 + ,3 + ,8541 + ,45 + ,17936 + ,151 + ,1 + ,520 + ,10 + ,58650 + ,291650 + ,383 + ,9165 + ,176 + ,0 + ,304 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,231 + ,0 + ,421 + ,1 + ,1013 + ,2 + ,69770 + ,253468 + ,402 + ,26734 + ,53 + ,-14 + ,217 + ,1 + ,15 + ,8 + ,60831 + ,315688 + ,164 + ,14461 + ,7687 + ,0 + ,353 + ,1 + ,587 + ,12 + ,59595 + ,325699 + ,126 + ,10475 + ,214 + ,0 + ,58 + ,1 + ,5371 + ,66 + ,87720 + ,446211 + ,45 + ,3730 + ,46 + ,49 + ,359 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,251 + ,0 + ,77 + ,1 + ,1012 + ,47 + ,114768 + ,368078 + ,82 + ,3576 + ,166 + ,26 + ,228 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,215 + ,0 + ,83 + ,1 + ,876 + ,59 + ,138971 + ,352850 + ,96 + ,2591 + ,174 + ,24 + ,1 + ,0 + ,162556 + ,1081 + ,213118 + ,6282154 + ,1 + ,5626 + ,37 + ,56289 + ,14 + ,1 + ,43556 + ,247 + ,32648 + ,227132 + ,414 + ,110 + ,1 + ,929 + ,395 + ,0 + ,3425 + ,12 + ,83620 + ,283910 + ,388 + ,6993 + ,25 + ,-3 + ,78 + ,1 + ,810 + ,43 + ,74015 + ,236761 + ,409 + ,855 + ,45 + ,26 + ,190 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,184 + ,0 + ,311 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,244 + ,0 + ,39 + ,0 + ,2365 + ,79 + ,191778 + ,550608 + ,32 + ,4438 + ,148 + ,120 + ,101 + ,0 + ,1261 + ,33 + ,76114 + ,307528 + ,343 + ,3258 + ,85 + ,13 + ,256 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,238 + ,0 + ,123 + ,1 + ,1585 + ,107 + ,93099 + ,355864 + ,92 + ,1457 + ,98 + ,2 + ,24 + ,1 + ,16189 + ,295 + ,116384 + ,358589 + ,89 + ,538 + ,10 + ,379 + ,284 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,242 + ,0 + ,273 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,265 + ,0 + ,121 + ,1 + ,10579 + ,92 + ,110309 + ,375195 + ,76 + ,1904 + ,17 + ,4 + ,214 + ,1 + ,474 + ,20 + ,61977 + ,288985 + ,386 + ,4449 + ,188 + ,0 + ,348 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,178 + ,0 + ,87 + ,1 + ,3642 + ,86 + ,90262 + ,458343 + ,42 + ,3004 + ,71 + ,22 + ,170 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,0 + ,200 + ,0 + ,472 + ,20 + ,80045 + ,269587 + ,397 + ,3479 + ,148 + ,0 + ,158 + ,1 + ,98 + ,2 + ,61490 + ,315236 + ,298 + ,57618 + ,1179 + ,0 + ,427 + ,1 + ,3999 + ,30 + ,51252 + ,42754 + ,428 + ,-5242 + ,-39 + ,-58 + ,233 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,216 + ,0 + ,345 + ,1 + ,621 + ,48 + ,60720 + ,308256 + ,341 + ,2255 + ,174 + ,0 + ,375 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,266 + ,0 + ,173 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,189 + ,0 + ,301 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,320 + ,28291 + ,3808 + ,0 + ,368 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,291 + ,0 + ,106 + ,1 + ,746 + ,23 + ,71561 + ,269661 + ,396 + ,3029 + ,93 + ,10 + ,208 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,228 + ,0 + ,282 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,255 + ,0 + ,249 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,260 + ,0 + ,240 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,194 + ,0 + ,52 + ,0 + ,4150 + ,150 + ,134759 + ,518365 + ,35 + ,2122 + ,77 + ,66 + ,55 + ,0 + ,4658 + ,158 + ,156608 + ,233773 + ,411 + ,214 + ,7 + ,58 + ,84 + ,1 + ,814 + ,18 + ,39625 + ,301881 + ,358 + ,5660 + ,125 + ,23 + ,143 + ,1 + ,1002 + ,15 + ,56750 + ,298568 + ,364 + ,6571 + ,98 + ,0 + ,205 + ,0 + ,496 + ,17 + ,87390 + ,325479 + ,128 + ,7381 + ,253 + ,0 + ,141 + ,1 + ,389 + ,14 + ,58990 + ,325506 + ,127 + ,8965 + ,322 + ,0 + ,12 + ,1 + ,12679 + ,110 + ,48020 + ,984885 + ,13 + ,7135 + ,62 + ,1871 + ,177 + ,1 + ,400 + ,152 + ,60720 + ,313267 + ,319 + ,745 + ,283 + ,0 + ,238 + ,1 + ,53 + ,6 + ,60349 + ,315793 + ,163 + ,19299 + ,2197 + ,0 + ,331 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,278 + ,0 + ,426 + ,0 + ,5109 + ,97 + ,67038 + ,215362 + ,417 + ,158 + ,3 + ,-31 + ,104 + ,1 + ,576 + ,18 + ,113761 + ,314073 + ,313 + ,6337 + ,198 + ,12 + ,409 + ,1 + ,438 + ,11 + ,58320 + ,298096 + ,366 + ,8918 + ,224 + ,-5 + ,292 + ,1 + ,165 + ,24 + ,62841 + ,325176 + ,130 + ,5216 + ,759 + ,0 + ,40 + ,1 + ,4069 + ,145 + ,79804 + ,207393 + ,420 + ,51 + ,2 + ,119 + ,185 + ,0 + ,23 + ,5 + ,62555 + ,314806 + ,306 + ,22961 + ,4974 + ,0 + ,59 + ,1 + ,1285 + ,107 + ,99489 + ,341340 + ,106 + ,1321 + ,110 + ,49 + ,82 + ,0 + ,4677 + ,75 + ,90131 + ,426280 + ,52 + ,3017 + ,48 + ,24 + ,15 + ,0 + ,24811 + ,505 + ,95350 + ,929118 + ,15 + ,1444 + ,29 + ,822 + ,253 + ,1 + ,167 + ,2 + ,64245 + ,307322 + ,345 + ,53661 + ,642 + ,0 + ,342 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,187 + ,0 + ,45 + ,1 + ,4628 + ,62 + ,69159 + ,387475 + ,64 + ,3024 + ,41 + ,85 + ,189 + ,1 + ,226 + ,5 + ,65745 + ,291787 + ,382 + ,18357 + ,406 + ,0 + ,76 + ,0 + ,1765 + ,20 + ,77623 + ,247060 + ,404 + ,2353 + ,27 + ,28 + ,145 + ,1 + ,460 + ,15 + ,63346 + ,329784 + ,122 + ,8652 + ,282 + ,0 + ,144 + ,1 + ,36 + ,3 + ,60894 + ,315834 + ,162 + ,38611 + ,3215 + ,0 + ,408 + ,1 + ,989 + ,11 + ,58930 + ,304555 + ,350 + ,9505 + ,106 + ,-5 + ,96 + ,0 + ,3055 + ,43 + ,60247 + ,376641 + ,75 + ,4108 + ,58 + ,14 + ,265 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,288 + ,0 + ,161 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,212 + ,0 + ,246 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,236 + ,0 + ,237 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,202 + ,0 + ,230 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,214 + ,0 + ,66 + ,1 + ,13253 + ,947 + ,90829 + ,357760 + ,90 + ,167 + ,12 + ,35 + ,244 + ,1 + ,24 + ,8 + ,59818 + ,315637 + ,168 + ,14455 + ,4772 + ,0 + ,153 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,208 + ,0 + ,199 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,261 + ,0 + ,222 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,198 + ,0 + ,363 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,289 + ,0 + ,387 + ,1 + ,131 + ,6 + ,59661 + ,327071 + ,123 + ,21179 + ,969 + ,0 + ,85 + ,1 + ,514 + ,15 + ,102725 + ,377516 + ,73 + ,11834 + ,346 + ,23 + ,103 + ,1 + ,3366 + ,80 + ,108479 + ,299243 + ,362 + ,1241 + ,29 + ,12 + ,60 + ,0 + ,9327 + ,101 + ,87419 + ,387699 + ,63 + ,1858 + ,20 + ,44 + ,309 + ,1 + ,384 + ,22 + ,54683 + ,309836 + ,333 + ,4993 + ,286 + ,0 + ,20 + ,1 + ,17821 + ,373 + ,122844 + ,444477 + ,46 + ,655 + ,14 + ,643 + ,391 + ,1 + ,397 + ,10 + ,62710 + ,322327 + ,136 + ,12233 + ,308 + ,-1 + ,139 + ,1 + ,897 + ,5 + ,60720 + ,314913 + ,303 + ,22983 + ,128 + ,0 + ,370 + ,1 + ,218 + ,21 + ,60720 + ,315553 + ,170 + ,5503 + ,531 + ,0 + ,31 + ,0 + ,3369 + ,28 + ,87161 + ,688779 + ,23 + ,17456 + ,145 + ,217 + ,109 + ,1 + ,5702 + ,61 + ,101481 + ,321896 + ,138 + ,1998 + ,21 + ,10 + ,41 + ,1 + ,8636 + ,855 + ,128294 + ,301607 + ,359 + ,119 + ,12 + ,108 + ,215 + ,1 + ,534 + ,7 + ,62620 + ,304485 + ,351 + ,14926 + ,196 + ,0 + ,332 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,0 + ,310 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,227 + ,0 + ,410 + ,1 + ,726 + ,3 + ,69980 + ,231861 + ,412 + ,10620 + ,44 + ,-5 + ,81 + ,1 + ,1380 + ,26 + ,60982 + ,347385 + ,98 + ,5669 + ,107 + ,25 + ,365 + ,1 + ,180 + ,11 + ,59635 + ,316386 + ,158 + ,10581 + ,647 + ,0 + ,53 + ,0 + ,7285 + ,115 + ,188873 + ,491303 + ,39 + ,2533 + ,40 + ,59 + ,67 + ,1 + ,880 + ,74 + ,80791 + ,261216 + ,399 + ,827 + ,70 + ,34 + ,335 + ,1 + ,1 + ,2 + ,60727 + ,315388 + ,175 + ,57694 + ,117743 + ,0 + ,380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,271 + ,0 + ,334 + ,1 + ,96 + ,5 + ,60379 + ,313729 + ,316 + ,22746 + ,1188 + ,0 + ,95 + ,1 + ,1889 + ,45 + ,37527 + ,358649 + ,88 + ,3526 + ,84 + ,14 + ,8 + ,0 + ,45187 + ,353 + ,234817 + ,1926517 + ,4 + ,4891 + ,38 + ,2854 + ,221 + ,1 + ,288 + ,4 + ,60510 + ,296656 + ,373 + ,24164 + ,335 + ,0 + ,419 + ,1 + ,1270 + ,26 + ,69206 + ,275311 + ,393 + ,2897 + ,59 + ,-10 + ,424 + ,1 + ,6526 + ,26 + ,55830 + ,-42143 + ,430 + ,-9313 + ,-37 + ,-28 + ,327 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,252 + ,0 + ,263 + ,1 + ,226 + ,12 + ,72835 + ,343929 + ,102 + ,11994 + ,636 + ,0 + ,394 + ,1 + ,694 + ,8 + ,68060 + ,367655 + ,83 + ,20957 + ,241 + ,-2 + ,286 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,241 + ,0 + ,283 + ,1 + ,249 + ,72 + ,56726 + ,313491 + ,317 + ,1576 + ,455 + ,0) + ,dim=c(10 + ,431) + ,dimnames=list(c('Rank' + ,'group' + ,'Costs' + ,'Trades' + ,'Dividends' + ,'Wealth' + ,'WRank' + ,'Profit/Trades' + ,'Profit/Cost' + ,'CCScore ') + ,1:431)) > y <- array(NA,dim=c(10,431),dimnames=list(c('Rank','group','Costs','Trades','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]) + } + } > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '10' > ylab = '' > xlab = '' > main = '' > #'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 group Costs Trades Dividends Wealth WRank 1 0 350 1 1595 17 60720 319369 143 2 39 61 0 5565 47 94355 493408 38 3 0 287 0 601 6 60720 319210 144 4 0 268 1 188 11 77655 381180 67 5 305 25 1 7146 105 134028 297978 367 6 -9 418 0 1135 17 62285 290476 384 7 -1 392 1 450 8 59325 292136 380 8 0 131 1 34 8 60630 314353 309 9 0 316 1 133 14 65990 339445 109 10 0 347 1 119 6 59118 303677 354 11 33 68 0 2053 53 100423 397144 60 12 32 70 0 4036 63 100269 424898 53 13 1 147 1 0 0 60720 315380 203 14 1 4655 393 60720 341570 105 360 30 15 1 131 11 61808 308989 338 9908 834 16 1 1766 73 60438 305959 349 1451 60 17 1 312 14 58598 318690 147 8478 381 18 1 448 40 59781 323361 132 3084 275 19 1 115 13 60945 318903 145 9146 1037 20 1 60 10 61124 314049 314 11405 1916 21 35 0 0 60720 315380 180 0 258 22 0 47705 298466 365 2813 271 0 343 23 438493 60720 315380 183 0 398 0 1442 24 378049 49 2021 165 -3 115 0 1389 25 313332 71 19783 128 6 188 1 149 26 319864 318 22666 761 0 379 1 2212 27 430866 142 8562 54 0 48 0 7489 28 315380 50 717 31 70 366 1 0 29 5106 287 0 117 1 402 26 84990 30 2999 330 5 412 0 7419 29 63255 31 0 12 -5 157 1 0 0 60720 32 0 264 1 307 15 58856 331955 116 33 58 54 0 1134 162 146216 527021 34 34 37 64 0 7561 87 82425 364304 86 35 9 113 1 5131 71 86111 292154 379 36 0 148 1 9 6 60835 314987 300 37 -31 425 1 2264 24 55174 272713 394 38 2865 7 1 40949 481 129352 1398893 9 39 0 197 0 4980 91 113521 409642 55 40 14 97 1 272 8 112995 429112 51 41 0 168 1 757 19 45689 382712 66 42 32 69 1 1424 94 131116 374943 77 43 0 127 1 0 0 60720 315380 191 44 1 378 52 71873 297413 369 1873 258 45 1 1061 5 70415 304252 353 20850 98 46 9 0 0 60720 315380 290 0 174 47 25 61938 309596 334 12177 541 0 46 48 11 101792 377305 74 7092 48 74 90 49 282 103345 386212 65 16928 329 20 27 50 40 74996 657954 24 1624 17 247 105 51 17 68696 368186 81 4205 100 11 413 52 64 57320 269753 395 4103 60 -6 26 53 159 63838 585715 30 6027 24 284 16 54 248 151352 1071292 11 5480 133 819 19 55 5 37884 992426 12 3195 21 652 336 56 36 60720 315371 294 23074 367 0 110 57 72 72954 249898 403 1386 26 10 86 58 0 64239 357312 91 2185 28 22 294 59 315380 60720 315380 257 0 322 1 0 60 2198 177 0 94 0 6533 75 126910 61 222 25 16 400 0 2126 138 21001 62 115877 14 -3 150 1 14 1 61016 63 -1243 8041 0 107 0 3192 60 59831 64 38294 -23 10 166 1 47 3 60720 65 787 2432 0 47 0 4891 217 210568 66 13007 35 73 182 0 531 10 71641 67 31096 245 0 245 1 2332 4 72680 68 0 53 0 330 1 0 0 60720 69 1974 10 0 16275 109 47818 1443586 7 70 -129 429 1 9252 48 79510 253537 401 71 28328 2 0 29790 309 81767 4321023 2 72 0 235 1 0 0 60720 315380 204 73 1 2300 29 64873 352108 97 5245 66 74 0 267 32 57640 330059 121 4064 487 75 1 382 21 59365 340968 107 6713 369 76 1 397 28 60646 317736 152 4205 297 77 0 11920 182 31080 209458 419 52 1 78 1 54660 302 126942 1491348 6 4276 24 79 1 5 2 60720 314887 304 57443 24814 80 0 0 0 60720 315380 218 0 219 81 353058 60720 315380 199 0 125 0 2647 82 315380 95 4137 58 1 216 1 0 83 0 201 0 196 0 0 0 60720 84 0 128 1 94 7 60793 314533 308 85 0 223 1 3 2 60698 315354 296 86 0 165 1 454 51 63261 302187 357 87 0 373 1 227 9 76644 336639 110 88 1 324 0 0 0 60720 315380 245 89 1 206 9 53110 296702 372 10745 469 90 1 25830 115 245546 1073089 10 7592 34 91 1 1528 38 60326 146494 426 -1408 -35 92 1 271 2 60720 325249 129 62625 463 93 1 138 6 69817 331420 117 21903 954 94 26 0 0 60720 315380 211 0 226 95 30 61404 314922 302 4420 413 0 276 96 80 62452 320016 141 4001 425 0 383 97 99 41477 280398 390 1005 141 0 93 98 2 63593 452469 43 2550 112 17 302 99 16 58790 301164 360 50582 348 0 385 100 2 62700 317330 154 7333 1511 0 225 101 0 60805 315576 169 57788 5924 0 275 102 -7170 60720 315380 284 0 431 1 18073 103 322331 429 -1151 -11 -366 220 1 1866 104 1629616 135 4218 66 0 6 0 42634 105 292754 5 9164 34 2964 290 1 249 106 318056 378 1496 373 0 129 1 422 107 355178 150 4216 280 0 65 1 2675 108 204325 93 1964 58 36 416 1 965 109 315380 421 1442 4 -8 299 1 0 110 7315 224 0 407 1 621 16 68005 111 0 188 -5 252 0 0 0 60720 112 0 257 1 365 8 54820 309560 335 113 38 63 0 3122 34 94670 414462 54 114 601 21 0 12988 119 82340 857217 16 115 132 38 0 5336 81 112477 697458 22 116 66 51 0 3160 108 108094 530670 33 117 0 386 1 2489 40 67804 238125 408 118 0 325 1 0 0 60720 315380 273 119 1 6984 45 80570 741409 19 12031 78 120 1 2001 52 95551 393343 62 3718 97 121 0 15236 90 77440 372631 78 1918 11 122 60 0 0 60720 315380 239 0 382 123 0 73433 317291 155 1955 221 0 369 124 306275 60720 315380 292 0 428 1 35624 125 317892 347 477 3 -70 100 0 1383 126 334280 151 10717 85 13 122 1 875 127 315380 112 9591 153 3 198 0 0 128 0 259 0 289 0 0 0 60720 129 1 307 0 0 0 60720 315380 222 130 0 72 3 60398 314210 312 38070 1591 131 0 265 31 60720 306948 346 3450 404 132 64175 0 60720 315380 186 0 277 1 133 60720 320398 140 3648 359 0 298 0 134 269 315380 230 0 356 1 0 0 135 67 0 43 1 4525 226 93811 501749 136 1 101 71 0 3045 58 27330 202055 137 212 31 133 1 0 0 60720 315380 138 638 1 0 0 60720 315380 206 0 139 0 14 60370 333210 114 9515 209 0 140 60436 0 60720 315380 283 0 378 1 141 55637 322340 134 11122 201 0 132 1 142 67440 369448 80 2690 109 0 279 0 143 60720 291841 381 785 69 0 376 1 144 282 315380 280 0 267 0 0 0 145 159 0 403 1 611 8 59190 296919 146 388 -4 303 1 0 0 60720 315380 147 420 1 923 8 58620 309038 337 13630 148 114 1 661 3 65920 246541 405 15514 149 305 0 2397 14 74020 289513 385 6394 150 346 1 366 53 61808 344425 101 2725 151 135 1 0 0 60720 315380 181 0 152 1659 2 62065 314210 311 57105 847 0 153 316 29 107577 480382 40 9668 169 2 154 0 9 60505 315009 299 12779 364 0 155 56535 0 60720 315380 197 0 146 1 156 64107 312878 322 10262 365 0 314 0 157 60720 322031 137 15254 2500 0 241 1 158 29 315380 200 0 34 0 4519 49 159 267 8118 88 163 384 1 0 0 160 138 0 120 1 837 69 61262 315688 161 35 4 72 1 5119 49 39039 378525 162 88 31 79 1 1280 117 69465 312378 163 79 26 99 1 2564 22 130140 403560 164 152 13 108 0 2045 20 97890 510834 165 6 10 119 0 2234 16 77200 214215 166 36 5 102 1 975 23 90534 235133 167 126 13 209 0 1136 32 48522 343613 168 366 0 293 1 453 21 81125 365959 169 227 0 211 1 0 0 60720 315380 170 234 1 61 14 60720 314551 307 8182 171 247 1 368 17 61686 303230 355 6072 172 4901 0 0 0 60720 315380 233 0 173 540 84 121920 469107 41 3204 55 67 174 0 14 103960 354228 94 11016 286 22 175 60720 0 60720 315380 253 0 160 1 176 286 315380 213 0 271 0 0 0 177 63057 0 306 1 2 9 60735 315394 178 3146 0 172 1 36 4 61564 312412 179 172 0 402 1 776 9 64230 333505 180 0 -4 4 1 84738 588 -26007 223193 181 315 4145 156 1 0 0 60720 315380 182 399 1 3 4 60761 315656 167 28914 183 381 1 529 9 63870 296261 374 10696 184 405 1 0 0 60720 315380 264 0 185 972 7 60845 336425 111 19489 337 0 186 0 34 71642 359335 87 4686 164 0 187 106611 0 60720 315380 232 0 89 1 188 48022 308636 339 1697 52 20 112 0 189 60720 158492 425 -883 -12 9 201 0 190 182 315380 243 0 344 1 0 0 191 23 0 22 0 22330 84 79801 711969 192 269 446 291 0 0 0 60720 315380 193 483 0 0 0 60720 315380 221 0 194 0 27 60830 306268 348 3936 220 0 195 88590 0 60720 315380 279 0 118 1 196 82903 442882 47 11566 108 5 92 0 197 60720 378509 70 4354 61 18 300 1 198 99 315380 223 0 312 1 365 83 199 310 1766 402 0 183 1 2461 57 200 17 2005 46 0 18 0 21950 519 201 416 1266 30 745 126 0 3294 23 202 295 748 5 1 193 1 141 6 203 342 19228 819 0 390 1 572 16 204 21 6746 189 -1 23 1 13326 102 205 423 4925 38 436 50 1 2284 33 206 160 -167 -2 66 206 0 10 2 207 207 58077 11639 0 162 1 0 0 208 92 0 98 0 1414 198 70184 330546 209 98 13 75 1 1975 35 118881 394510 210 2613 30 134 1 43 3 60887 312846 211 337 0 180 1 0 0 60720 315380 212 341 1 844 78 60925 296139 375 1233 213 154 1 304 15 62969 295580 376 6372 214 36 1 458 11 58625 297765 368 8888 215 192 1 18562 155 102313 377934 72 1148 216 7123 1 0 0 60720 315380 179 0 217 622 109 288170 638830 28 4026 62 769 218 174 73 73007 304376 352 1430 168 0 219 2220 1 64820 307424 344 107424 618 0 220 121 22 301670 644190 26 20190 200 154 221 11819 24 56178 295370 377 3974 788 0 222 0 85 106113 574339 31 4404 32 199 223 60798 0 60720 315380 235 0 274 1 224 70694 310201 332 9183 883 0 406 1 225 56364 327007 124 9770 107 -5 116 1 226 60720 343466 104 3679 95 6 320 1 227 248 315380 274 0 367 1 0 0 228 3912 0 288 1 30 4 62045 318098 229 75 0 414 0 3310 33 75230 448243 230 227 -6 74 1 554 11 79285 325738 231 224 30 250 0 0 0 60720 315380 232 397 1 468 21 60720 312161 328 5341 233 417 1 4917 42 52811 243650 406 1039 234 323 1 3256 66 35250 407159 57 3139 235 125 1 0 0 60720 315380 185 0 236 0 13 59734 317698 153 9054 944 0 237 60722 0 60720 315380 195 0 326 1 238 60720 312502 324 56251 5154 0 261 0 239 133 315380 262 0 357 1 514 52 240 226 2353 238 0 313 0 0 0 241 272 0 360 1 0 0 60720 315380 242 0 0 0 0 60720 315380 285 0 243 60720 0 60720 315380 276 0 355 1 244 268 315380 272 0 362 1 0 0 245 43 0 37 1 10327 91 88577 640273 246 11 135 42 1 13718 61 96448 345783 247 121 104 57 1 3748 9 50350 652925 248 17 53 30 1 14416 361 49857 439798 249 52 220 404 0 1526 25 69351 278990 250 210 -4 80 0 666 171 117869 339836 251 14 25 422 1 2844 67 72683 240897 252 318 -15 297 1 0 0 60720 315380 253 175 1 368 17 61167 297141 370 5714 254 333 1 0 0 60720 315380 205 0 255 26 13 70811 331323 118 10102 394 0 256 0 1 60896 313880 315 113880 4328 0 257 60720 0 60720 315380 237 0 167 1 258 336 315380 193 0 411 1 1303 61 259 297 1794 84 -5 155 1 20 3 260 58 38415 5674 0 56 1 2384 97 261 361 2123 86 54 204 0 203 27 262 159 3739 498 0 178 0 71 1 263 356 116176 1646 0 295 1 53 14 264 389 7315 1932 0 396 1 562 13 265 398 6430 149 -3 393 1 622 7 266 329 9039 102 -2 135 1 645 5 267 340 22415 174 0 364 0 1763 24 268 209 4514 61 0 179 1 0 0 269 311 0 231 1 317 9 63915 298700 270 84213 0 319 1 1 6 60719 315372 271 431 0 191 1 275 4 59500 318745 272 111 0 321 1 0 0 60720 315380 273 29 1 936 30 67939 408881 56 6963 274 423 1 8568 37 32168 786690 18 15856 275 202 0 11528 70 -14545 -83265 431 -4047 276 738 0 0 0 60720 315380 250 0 277 0 10 60720 321376 139 12138 164 0 278 64270 0 60720 315380 281 0 405 1 279 60951 276898 392 10985 130 -4 176 1 280 60743 315547 171 11555 921 0 236 1 281 60720 315487 172 115487 52257 0 142 1 282 8 315380 220 0 13 0 18014 239 283 256 5043 67 943 281 0 0 0 284 9 0 328 0 0 0 60720 315380 285 207 0 37704 452 60448 983660 14 1734 286 73 0 63 4 65688 318574 148 29643 287 242 1 1431 34 106885 310768 330 3258 288 229 1 94 3 61360 312887 321 37629 289 317 1 192 4 65276 312339 327 28085 290 44 1 32 1 59988 314964 301 114964 291 349 1 6869 105 117520 379983 68 1714 292 0 1 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 1 2328 94 296 157 570 23 0 308 0 209 13 297 166 8973 558 0 195 1 28 2 298 331 57844 4148 0 280 1 176 28 299 424 3952 629 0 415 1 1920 15 300 3 -2306 -18 -6 3 0 87550 458 301 383 8541 45 17936 151 1 520 10 302 231 9165 176 0 304 0 0 0 303 53 0 421 1 1013 2 69770 253468 304 7687 -14 217 1 15 8 60831 315688 305 214 0 353 1 587 12 59595 325699 306 46 0 58 1 5371 66 87720 446211 307 77 49 359 1 0 0 60720 315380 308 228 1 1012 47 114768 368078 82 3576 309 876 1 0 0 60720 315380 215 0 310 162556 59 138971 352850 96 2591 174 24 311 43556 1081 213118 6282154 1 5626 37 56289 312 3425 247 32648 227132 414 110 1 929 313 810 12 83620 283910 388 6993 25 -3 314 0 43 74015 236761 409 855 45 26 315 60720 0 60720 315380 184 0 311 0 316 32 315380 244 0 39 0 2365 79 317 343 4438 148 120 101 0 1261 33 318 238 3258 85 13 256 1 0 0 319 98 0 123 1 1585 107 93099 355864 320 10 2 24 1 16189 295 116384 358589 321 273 379 284 1 0 0 60720 315380 322 10579 0 0 0 60720 315380 265 0 323 474 92 110309 375195 76 1904 17 4 324 0 20 61977 288985 386 4449 188 0 325 90262 0 60720 315380 178 0 87 1 326 60720 458343 42 3004 71 22 170 1 327 397 315380 192 0 200 0 472 20 328 298 3479 148 0 158 1 98 2 329 428 57618 1179 0 427 1 3999 30 330 216 -5242 -39 -58 233 1 0 0 331 174 0 345 1 621 48 60720 308256 332 173 0 375 1 0 0 60720 315380 333 30 1 0 0 60720 315380 189 0 334 0 4 60798 313164 320 28291 3808 0 335 71561 0 60720 315380 291 0 106 1 336 60720 269661 396 3029 93 10 208 0 337 255 315380 228 0 282 1 0 0 338 240 0 249 0 0 0 60720 315380 339 4150 1 0 0 60720 315380 194 0 340 4658 150 134759 518365 35 2122 77 66 341 814 158 156608 233773 411 214 7 58 342 1002 18 39625 301881 358 5660 125 23 343 496 15 56750 298568 364 6571 98 0 344 389 17 87390 325479 128 7381 253 0 345 12679 14 58990 325506 127 8965 322 0 346 400 110 48020 984885 13 7135 62 1871 347 53 152 60720 313267 319 745 283 0 348 0 6 60349 315793 163 19299 2197 0 349 67038 0 60720 315380 278 0 426 0 350 113761 215362 417 158 3 -31 104 1 351 58320 314073 313 6337 198 12 409 1 352 62841 298096 366 8918 224 -5 292 1 353 79804 325176 130 5216 759 0 40 1 354 62555 207393 420 51 2 119 185 0 355 99489 314806 306 22961 4974 0 59 1 356 90131 341340 106 1321 110 49 82 0 357 95350 426280 52 3017 48 24 15 0 358 64245 929118 15 1444 29 822 253 1 359 60720 307322 345 53661 642 0 342 1 360 64 315380 187 0 45 1 4628 62 361 382 3024 41 85 189 1 226 5 362 404 18357 406 0 76 0 1765 20 363 122 2353 27 28 145 1 460 15 364 162 8652 282 0 144 1 36 3 365 350 38611 3215 0 408 1 989 11 366 75 9505 106 -5 96 0 3055 43 367 288 4108 58 14 265 0 0 0 368 246 0 161 1 0 0 60720 315380 369 0 0 0 0 60720 315380 236 0 370 60720 0 60720 315380 202 0 230 1 371 90 315380 214 0 66 1 13253 947 372 168 167 12 35 244 1 24 8 373 208 14455 4772 0 153 1 0 0 374 222 0 199 0 0 0 60720 315380 375 0 1 0 0 60720 315380 198 0 376 59661 0 60720 315380 289 0 387 1 377 102725 327071 123 21179 969 0 85 1 378 108479 377516 73 11834 346 23 103 1 379 87419 299243 362 1241 29 12 60 0 380 54683 387699 63 1858 20 44 309 1 381 122844 309836 333 4993 286 0 20 1 382 62710 444477 46 655 14 643 391 1 383 60720 322327 136 12233 308 -1 139 1 384 60720 314913 303 22983 128 0 370 1 385 87161 315553 170 5503 531 0 31 0 386 101481 688779 23 17456 145 217 109 1 387 128294 321896 138 1998 21 10 41 1 388 62620 301607 359 119 12 108 215 1 389 60720 304485 351 14926 196 0 332 1 390 227 315380 275 0 310 0 0 0 391 44 0 410 1 726 3 69980 231861 392 107 -5 81 1 1380 26 60982 347385 393 647 25 365 1 180 11 59635 316386 394 40 0 53 0 7285 115 188873 491303 395 70 59 67 1 880 74 80791 261216 396 117743 34 335 1 1 2 60727 315388 397 334 0 380 1 0 0 60720 315380 398 95 1 96 5 60379 313729 316 22746 399 8 1 1889 45 37527 358649 88 3526 400 221 0 45187 353 234817 1926517 4 4891 401 419 1 288 4 60510 296656 373 24164 402 424 1 1270 26 69206 275311 393 2897 403 327 1 6526 26 55830 -42143 430 -9313 404 226 0 0 0 60720 315380 252 0 405 694 12 72835 343929 102 11994 636 0 406 0 8 68060 367655 83 20957 241 -2 407 56726 0 60720 315380 241 0 283 1 408 60720 313491 317 1576 455 0 350 1 409 94355 319369 143 7022 75 0 61 0 410 60720 493408 38 6243 53 39 287 0 411 77655 319210 144 19868 198 0 268 1 412 134028 381180 67 16471 964 0 25 1 413 62285 297978 367 933 14 305 418 0 414 59325 290476 384 5322 80 -9 392 1 415 60630 292136 380 11517 205 -1 131 1 416 65990 314353 309 14294 3340 0 316 1 417 59118 339445 109 9960 1052 0 347 1 418 100423 303677 354 17280 872 0 68 0 419 100269 397144 60 3720 96 33 70 0 420 60720 424898 53 3570 56 32 147 1 421 105 315380 203 0 169 1 4655 393 422 338 360 30 0 248 1 131 11 423 349 9908 834 0 152 1 1766 73 424 147 1451 60 0 351 1 312 14 425 132 8478 381 0 372 1 448 40 426 145 3084 275 0 137 1 115 13 427 314 9146 1037 0 352 1 60 10 428 180 11405 1916 0 338 1 0 0 429 271 0 258 1 364 35 47705 298466 430 398 0 343 1 0 0 60720 315380 431 115 0 1442 118 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 1 364 22 1 0 23 118 121173 24 9 57530 25 5 62041 26 14 60720 27 322 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 1 203 47 0 3689 48 1 567 49 0 27664 50 1 1686 51 0 1164 52 1 15824 53 1 6575 54 0 37597 55 1 314 56 0 1932 57 0 5692 58 1 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 1 0 81 37 62920 82 0 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 1 278 95 1 282 96 1 571 97 1 2253 98 1 290 99 1 78 100 1 20 101 1 0 102 180 27284 103 29 56225 104 156 157214 105 62 54323 106 28 57935 107 79 59017 108 3 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 1 530 123 1 0 124 223 157278 125 11 73221 126 14 67000 127 0 60720 128 315380 254 129 0 389 130 0 210 131 194 1 132 335 33 133 0 0 134 60720 315380 135 37 1335 136 422 35 137 196 0 138 186 0 139 262 0 140 607 11 141 1558 63 142 1324 117 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 1 152 124 0 153 136 1 154 130 1 155 309 11 156 49 8 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 1 173 88 1 174 329 1 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 1 185 181 1 186 278 1 187 2099 64 188 3437 47 189 0 0 190 60720 315380 191 20 6095 192 258 0 193 164 1 194 358 1 195 2239 21 196 2949 41 197 0 0 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 1 217 254 1 218 374 1 219 35 1 220 377 1 221 32 0 222 260 1 223 125 12 224 1182 13 225 1503 39 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 1 236 239 1 237 22 2 238 0 0 239 78780 322378 240 60720 315380 241 270 0 242 333 1 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 1 255 339 1 256 243 1 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 1 277 266 1 278 592 7 279 126 10 280 2 1 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 1 293 1 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 1 310 1 0 311 14 1 312 395 0 313 78 1 314 190 1 315 0 0 316 191778 550608 317 76114 307528 318 60720 315380 319 92 1457 320 89 538 321 242 0 322 121 1 323 214 1 324 348 1 325 3642 86 326 0 0 327 80045 269587 328 61490 315236 329 51252 42754 330 60720 315380 331 341 2255 332 266 0 333 301 0 334 368 1 335 746 23 336 0 0 337 60720 315380 338 260 0 339 52 0 340 55 0 341 84 1 342 143 1 343 205 0 344 141 1 345 12 1 346 177 1 347 238 1 348 331 1 349 5109 97 350 576 18 351 438 11 352 165 24 353 4069 145 354 23 5 355 1285 107 356 4677 75 357 24811 505 358 167 2 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 1 370 0 0 371 90829 357760 372 59818 315637 373 60720 315380 374 261 0 375 363 1 376 131 6 377 514 15 378 3366 80 379 9327 101 380 384 22 381 17821 373 382 397 10 383 897 5 384 218 21 385 3369 28 386 5702 61 387 8636 855 388 534 7 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 1 405 394 1 406 286 1 407 249 72 408 1595 17 409 5565 47 410 601 6 411 188 11 412 7146 105 413 1135 17 414 450 8 415 34 8 416 133 14 417 119 6 418 2053 53 419 4036 63 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 group Costs 6.002e+04 8.094e-03 -4.399e-02 -1.298e-02 Trades Dividends Wealth WRank -1.080e-01 -1.013e-01 -3.407e-02 -1.288e-01 `Profit/Trades` `Profit/Cost` -2.211e-01 -1.892e-02 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -69824 -40287 -21000 -1048 1578282 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 6.002e+04 9.178e+03 6.540 1.79e-10 *** Rank 8.094e-03 4.095e-02 0.198 0.84339 group -4.399e-02 4.766e-02 -0.923 0.35657 Costs -1.298e-02 1.573e-02 -0.825 0.40973 Trades -1.080e-01 4.264e-02 -2.533 0.01166 * Dividends -1.013e-01 3.609e-02 -2.806 0.00525 ** Wealth -3.407e-02 2.027e-02 -1.681 0.09346 . WRank -1.288e-01 4.271e-02 -3.015 0.00273 ** `Profit/Trades` -2.211e-01 8.480e-02 -2.607 0.00945 ** `Profit/Cost` -1.892e-02 2.171e-02 -0.871 0.38399 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 103300 on 421 degrees of freedom Multiple R-squared: 0.06127, Adjusted R-squared: 0.04121 F-statistic: 3.053 on 9 and 421 DF, p-value: 0.001486 > 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,] 3.315981e-12 6.631962e-12 1.000000e+00 [2,] 1.838217e-15 3.676433e-15 1.000000e+00 [3,] 1.886238e-18 3.772477e-18 1.000000e+00 [4,] 3.931924e-22 7.863848e-22 1.000000e+00 [5,] 3.996943e-25 7.993885e-25 1.000000e+00 [6,] 9.288183e-29 1.857637e-28 1.000000e+00 [7,] 2.568138e-32 5.136275e-32 1.000000e+00 [8,] 6.263582e-36 1.252716e-35 1.000000e+00 [9,] 1.402571e-39 2.805142e-39 1.000000e+00 [10,] 1.074412e-42 2.148823e-42 1.000000e+00 [11,] 3.329469e-19 6.658937e-19 1.000000e+00 [12,] 2.819928e-09 5.639857e-09 1.000000e+00 [13,] 8.808239e-09 1.761648e-08 1.000000e+00 [14,] 2.600855e-09 5.201709e-09 1.000000e+00 [15,] 4.582505e-08 9.165010e-08 1.000000e+00 [16,] 3.559604e-08 7.119209e-08 1.000000e+00 [17,] 4.677024e-07 9.354048e-07 9.999995e-01 [18,] 2.340538e-07 4.681076e-07 9.999998e-01 [19,] 7.414042e-08 1.482808e-07 9.999999e-01 [20,] 2.417730e-08 4.835460e-08 1.000000e+00 [21,] 1.120352e-07 2.240703e-07 9.999999e-01 [22,] 3.916224e-08 7.832448e-08 1.000000e+00 [23,] 1.797756e-08 3.595512e-08 1.000000e+00 [24,] 1.346892e-08 2.693784e-08 1.000000e+00 [25,] 6.266699e-09 1.253340e-08 1.000000e+00 [26,] 8.417751e-09 1.683550e-08 1.000000e+00 [27,] 2.986874e-09 5.973748e-09 1.000000e+00 [28,] 1.536508e-09 3.073016e-09 1.000000e+00 [29,] 5.265065e-10 1.053013e-09 1.000000e+00 [30,] 2.070278e-10 4.140556e-10 1.000000e+00 [31,] 6.976151e-11 1.395230e-10 1.000000e+00 [32,] 8.137766e-11 1.627553e-10 1.000000e+00 [33,] 3.508641e-11 7.017281e-11 1.000000e+00 [34,] 1.162761e-11 2.325521e-11 1.000000e+00 [35,] 3.823849e-12 7.647697e-12 1.000000e+00 [36,] 1.413595e-12 2.827191e-12 1.000000e+00 [37,] 1.774920e-12 3.549841e-12 1.000000e+00 [38,] 5.803650e-13 1.160730e-12 1.000000e+00 [39,] 1.883585e-13 3.767169e-13 1.000000e+00 [40,] 2.191561e-13 4.383122e-13 1.000000e+00 [41,] 8.022748e-14 1.604550e-13 1.000000e+00 [42,] 2.973266e-14 5.946532e-14 1.000000e+00 [43,] 9.562899e-15 1.912580e-14 1.000000e+00 [44,] 3.227526e-15 6.455051e-15 1.000000e+00 [45,] 1.035481e-15 2.070961e-15 1.000000e+00 [46,] 3.373050e-16 6.746100e-16 1.000000e+00 [47,] 1.214614e-14 2.429227e-14 1.000000e+00 [48,] 4.522923e-15 9.045846e-15 1.000000e+00 [49,] 1.584802e-15 3.169603e-15 1.000000e+00 [50,] 2.673733e-14 5.347466e-14 1.000000e+00 [51,] 1.273826e-14 2.547651e-14 1.000000e+00 [52,] 4.865133e-15 9.730266e-15 1.000000e+00 [53,] 2.061610e-15 4.123220e-15 1.000000e+00 [54,] 7.611461e-16 1.522292e-15 1.000000e+00 [55,] 2.906335e-16 5.812670e-16 1.000000e+00 [56,] 1.277609e-16 2.555219e-16 1.000000e+00 [57,] 5.624898e-17 1.124980e-16 1.000000e+00 [58,] 2.261135e-17 4.522270e-17 1.000000e+00 [59,] 2.459134e-17 4.918268e-17 1.000000e+00 [60,] 9.022847e-18 1.804569e-17 1.000000e+00 [61,] 3.200309e-18 6.400617e-18 1.000000e+00 [62,] 1.099609e-18 2.199217e-18 1.000000e+00 [63,] 3.774806e-19 7.549612e-19 1.000000e+00 [64,] 1.259363e-19 2.518726e-19 1.000000e+00 [65,] 4.194010e-20 8.388020e-20 1.000000e+00 [66,] 9.420983e-20 1.884197e-19 1.000000e+00 [67,] 3.144649e-20 6.289299e-20 1.000000e+00 [68,] 1.049495e-20 2.098989e-20 1.000000e+00 [69,] 1.644710e-18 3.289419e-18 1.000000e+00 [70,] 1.748506e-18 3.497012e-18 1.000000e+00 [71,] 7.881604e-19 1.576321e-18 1.000000e+00 [72,] 3.139831e-19 6.279662e-19 1.000000e+00 [73,] 6.991183e-17 1.398237e-16 1.000000e+00 [74,] 2.838995e-17 5.677990e-17 1.000000e+00 [75,] 1.116764e-17 2.233527e-17 1.000000e+00 [76,] 4.508142e-18 9.016284e-18 1.000000e+00 [77,] 1.788354e-18 3.576708e-18 1.000000e+00 [78,] 2.434505e-18 4.869009e-18 1.000000e+00 [79,] 9.798040e-19 1.959608e-18 1.000000e+00 [80,] 3.918616e-19 7.837233e-19 1.000000e+00 [81,] 1.529137e-19 3.058274e-19 1.000000e+00 [82,] 6.000621e-20 1.200124e-19 1.000000e+00 [83,] 2.386344e-20 4.772688e-20 1.000000e+00 [84,] 9.347271e-21 1.869454e-20 1.000000e+00 [85,] 3.518313e-21 7.036627e-21 1.000000e+00 [86,] 1.404198e-21 2.808396e-21 1.000000e+00 [87,] 5.333938e-22 1.066788e-21 1.000000e+00 [88,] 2.056359e-22 4.112718e-22 1.000000e+00 [89,] 7.844992e-23 1.568998e-22 1.000000e+00 [90,] 4.772385e-22 9.544770e-22 1.000000e+00 [91,] 1.870210e-21 3.740421e-21 1.000000e+00 [92,] 1.000000e+00 2.423132e-13 1.211566e-13 [93,] 1.000000e+00 8.000754e-15 4.000377e-15 [94,] 1.000000e+00 4.678788e-17 2.339394e-17 [95,] 1.000000e+00 7.759251e-21 3.879626e-21 [96,] 1.000000e+00 9.802799e-24 4.901400e-24 [97,] 1.000000e+00 3.943937e-28 1.971969e-28 [98,] 1.000000e+00 4.470860e-28 2.235430e-28 [99,] 1.000000e+00 3.734466e-28 1.867233e-28 [100,] 1.000000e+00 8.361844e-28 4.180922e-28 [101,] 1.000000e+00 1.863083e-27 9.315415e-28 [102,] 1.000000e+00 3.970594e-27 1.985297e-27 [103,] 1.000000e+00 8.231260e-27 4.115630e-27 [104,] 1.000000e+00 1.703613e-26 8.518063e-27 [105,] 1.000000e+00 3.728014e-26 1.864007e-26 [106,] 1.000000e+00 8.121383e-26 4.060691e-26 [107,] 1.000000e+00 1.479668e-25 7.398338e-26 [108,] 1.000000e+00 3.134139e-25 1.567069e-25 [109,] 1.000000e+00 6.572327e-25 3.286163e-25 [110,] 1.000000e+00 1.382758e-24 6.913789e-25 [111,] 1.000000e+00 9.263832e-25 4.631916e-25 [112,] 1.000000e+00 1.424736e-38 7.123678e-39 [113,] 1.000000e+00 1.481674e-46 7.408372e-47 [114,] 1.000000e+00 1.578519e-59 7.892596e-60 [115,] 1.000000e+00 9.061812e-78 4.530906e-78 [116,] 1.000000e+00 5.868714e-80 2.934357e-80 [117,] 1.000000e+00 2.269263e-79 1.134632e-79 [118,] 1.000000e+00 7.908528e-79 3.954264e-79 [119,] 1.000000e+00 2.406353e-78 1.203176e-78 [120,] 1.000000e+00 1.695212e-78 8.476062e-79 [121,] 1.000000e+00 6.643561e-79 3.321780e-79 [122,] 1.000000e+00 6.494594e-105 3.247297e-105 [123,] 1.000000e+00 4.038238e-107 2.019119e-107 [124,] 1.000000e+00 1.097650e-106 5.488252e-107 [125,] 1.000000e+00 4.190755e-106 2.095378e-106 [126,] 1.000000e+00 2.091701e-106 1.045851e-106 [127,] 1.000000e+00 6.674694e-106 3.337347e-106 [128,] 1.000000e+00 1.195148e-105 5.975741e-106 [129,] 1.000000e+00 8.206565e-109 4.103283e-109 [130,] 1.000000e+00 2.911104e-110 1.455552e-110 [131,] 1.000000e+00 3.568408e-110 1.784204e-110 [132,] 1.000000e+00 5.528180e-117 2.764090e-117 [133,] 1.000000e+00 1.981787e-116 9.908934e-117 [134,] 1.000000e+00 9.509512e-116 4.754756e-116 [135,] 1.000000e+00 1.407899e-115 7.039493e-116 [136,] 1.000000e+00 5.385378e-115 2.692689e-115 [137,] 1.000000e+00 1.885173e-114 9.425864e-115 [138,] 1.000000e+00 5.907076e-114 2.953538e-114 [139,] 1.000000e+00 2.419498e-113 1.209749e-113 [140,] 1.000000e+00 9.235663e-113 4.617832e-113 [141,] 1.000000e+00 3.291313e-112 1.645656e-112 [142,] 1.000000e+00 1.091492e-111 5.457462e-112 [143,] 1.000000e+00 2.896456e-111 1.448228e-111 [144,] 1.000000e+00 2.748161e-111 1.374080e-111 [145,] 1.000000e+00 4.146433e-111 2.073216e-111 [146,] 1.000000e+00 1.112474e-118 5.562368e-119 [147,] 1.000000e+00 3.769695e-120 1.884848e-120 [148,] 1.000000e+00 1.900260e-119 9.501301e-120 [149,] 1.000000e+00 1.049001e-118 5.245004e-119 [150,] 1.000000e+00 5.555063e-118 2.777532e-118 [151,] 1.000000e+00 3.255175e-117 1.627587e-117 [152,] 1.000000e+00 1.914690e-116 9.573448e-117 [153,] 1.000000e+00 8.716401e-116 4.358200e-116 [154,] 1.000000e+00 4.184182e-115 2.092091e-115 [155,] 1.000000e+00 2.271830e-114 1.135915e-114 [156,] 1.000000e+00 1.278765e-113 6.393823e-114 [157,] 1.000000e+00 6.791273e-113 3.395637e-113 [158,] 1.000000e+00 2.963744e-112 1.481872e-112 [159,] 1.000000e+00 1.380898e-111 6.904490e-112 [160,] 1.000000e+00 6.647754e-111 3.323877e-111 [161,] 1.000000e+00 2.488656e-110 1.244328e-110 [162,] 1.000000e+00 9.008918e-110 4.504459e-110 [163,] 1.000000e+00 2.162044e-109 1.081022e-109 [164,] 1.000000e+00 2.279207e-109 1.139603e-109 [165,] 1.000000e+00 1.476462e-109 7.382312e-110 [166,] 1.000000e+00 7.637701e-109 3.818851e-109 [167,] 1.000000e+00 4.034567e-108 2.017284e-108 [168,] 1.000000e+00 2.101744e-107 1.050872e-107 [169,] 1.000000e+00 1.101557e-106 5.507784e-107 [170,] 1.000000e+00 4.805016e-106 2.402508e-106 [171,] 1.000000e+00 2.363316e-105 1.181658e-105 [172,] 1.000000e+00 1.136833e-104 5.684164e-105 [173,] 1.000000e+00 3.987541e-104 1.993770e-104 [174,] 1.000000e+00 1.325963e-103 6.629816e-104 [175,] 1.000000e+00 5.191885e-105 2.595943e-105 [176,] 1.000000e+00 1.969993e-104 9.849966e-105 [177,] 1.000000e+00 8.833305e-104 4.416653e-104 [178,] 1.000000e+00 1.043742e-103 5.218711e-104 [179,] 1.000000e+00 3.739996e-103 1.869998e-103 [180,] 1.000000e+00 1.923785e-102 9.618925e-103 [181,] 1.000000e+00 9.308332e-102 4.654166e-102 [182,] 1.000000e+00 3.034985e-101 1.517492e-101 [183,] 1.000000e+00 7.421390e-102 3.710695e-102 [184,] 1.000000e+00 1.312126e-101 6.560628e-102 [185,] 1.000000e+00 4.830983e-101 2.415491e-101 [186,] 1.000000e+00 9.001885e-101 4.500943e-101 [187,] 1.000000e+00 6.530081e-101 3.265041e-101 [188,] 1.000000e+00 7.461425e-102 3.730712e-102 [189,] 1.000000e+00 2.230172e-101 1.115086e-101 [190,] 1.000000e+00 6.117082e-101 3.058541e-101 [191,] 1.000000e+00 1.904527e-100 9.522635e-101 [192,] 1.000000e+00 9.894227e-101 4.947113e-101 [193,] 1.000000e+00 3.582890e-100 1.791445e-100 [194,] 1.000000e+00 1.292876e-99 6.464381e-100 [195,] 1.000000e+00 5.311498e-99 2.655749e-99 [196,] 1.000000e+00 2.807260e-98 1.403630e-98 [197,] 1.000000e+00 1.556455e-97 7.782275e-98 [198,] 1.000000e+00 8.150116e-97 4.075058e-97 [199,] 1.000000e+00 4.168357e-96 2.084178e-96 [200,] 1.000000e+00 2.054488e-95 1.027244e-95 [201,] 1.000000e+00 1.013310e-94 5.066552e-95 [202,] 1.000000e+00 4.911929e-94 2.455965e-94 [203,] 1.000000e+00 2.335083e-93 1.167541e-93 [204,] 1.000000e+00 1.203448e-92 6.017242e-93 [205,] 1.000000e+00 4.324463e-92 2.162231e-92 [206,] 1.000000e+00 1.371882e-91 6.859409e-92 [207,] 1.000000e+00 5.085453e-91 2.542726e-91 [208,] 1.000000e+00 1.056373e-90 5.281865e-91 [209,] 1.000000e+00 4.460692e-90 2.230346e-90 [210,] 1.000000e+00 1.076217e-89 5.381085e-90 [211,] 1.000000e+00 2.872804e-89 1.436402e-89 [212,] 1.000000e+00 1.393194e-88 6.965970e-89 [213,] 1.000000e+00 7.096044e-88 3.548022e-88 [214,] 1.000000e+00 3.625282e-87 1.812641e-87 [215,] 1.000000e+00 7.814022e-87 3.907011e-87 [216,] 1.000000e+00 3.956665e-86 1.978332e-86 [217,] 1.000000e+00 2.049884e-85 1.024942e-85 [218,] 1.000000e+00 1.013420e-84 5.067099e-85 [219,] 1.000000e+00 4.884536e-84 2.442268e-84 [220,] 1.000000e+00 2.322474e-83 1.161237e-83 [221,] 1.000000e+00 1.064204e-82 5.321019e-83 [222,] 1.000000e+00 4.526494e-82 2.263247e-82 [223,] 1.000000e+00 2.109304e-81 1.054652e-81 [224,] 1.000000e+00 5.696432e-81 2.848216e-81 [225,] 1.000000e+00 1.499135e-80 7.495677e-81 [226,] 1.000000e+00 7.412441e-80 3.706220e-80 [227,] 1.000000e+00 2.211977e-79 1.105988e-79 [228,] 1.000000e+00 8.304417e-79 4.152208e-79 [229,] 1.000000e+00 3.877927e-78 1.938964e-78 [230,] 1.000000e+00 1.770281e-77 8.851406e-78 [231,] 1.000000e+00 4.476907e-77 2.238454e-77 [232,] 1.000000e+00 8.931691e-77 4.465846e-77 [233,] 1.000000e+00 3.951733e-76 1.975867e-76 [234,] 1.000000e+00 1.911086e-75 9.555428e-76 [235,] 1.000000e+00 8.572965e-75 4.286483e-75 [236,] 1.000000e+00 4.167981e-74 2.083991e-74 [237,] 1.000000e+00 1.844081e-73 9.220405e-74 [238,] 1.000000e+00 8.618582e-73 4.309291e-73 [239,] 1.000000e+00 3.609757e-72 1.804879e-72 [240,] 1.000000e+00 1.621781e-71 8.108904e-72 [241,] 1.000000e+00 7.146826e-71 3.573413e-71 [242,] 1.000000e+00 3.120290e-70 1.560145e-70 [243,] 1.000000e+00 8.239169e-70 4.119584e-70 [244,] 1.000000e+00 2.365168e-69 1.182584e-69 [245,] 1.000000e+00 6.060742e-69 3.030371e-69 [246,] 1.000000e+00 1.426062e-68 7.130310e-69 [247,] 1.000000e+00 5.339105e-68 2.669553e-68 [248,] 1.000000e+00 2.071792e-67 1.035896e-67 [249,] 1.000000e+00 7.858181e-67 3.929091e-67 [250,] 1.000000e+00 3.057437e-66 1.528718e-66 [251,] 1.000000e+00 1.159205e-65 5.796025e-66 [252,] 1.000000e+00 4.568381e-65 2.284191e-65 [253,] 1.000000e+00 1.808247e-64 9.041237e-65 [254,] 1.000000e+00 7.218058e-64 3.609029e-64 [255,] 1.000000e+00 2.918375e-63 1.459187e-63 [256,] 1.000000e+00 1.169392e-62 5.846958e-63 [257,] 1.000000e+00 4.916432e-62 2.458216e-62 [258,] 1.000000e+00 5.792168e-63 2.896084e-63 [259,] 1.000000e+00 2.576641e-62 1.288320e-62 [260,] 1.000000e+00 1.132539e-61 5.662695e-62 [261,] 1.000000e+00 4.671460e-61 2.335730e-61 [262,] 1.000000e+00 1.250198e-60 6.250990e-61 [263,] 1.000000e+00 2.465842e-60 1.232921e-60 [264,] 1.000000e+00 1.025144e-59 5.125718e-60 [265,] 1.000000e+00 2.614859e-59 1.307430e-59 [266,] 1.000000e+00 5.738920e-59 2.869460e-59 [267,] 1.000000e+00 2.567571e-58 1.283785e-58 [268,] 1.000000e+00 1.133928e-57 5.669641e-58 [269,] 1.000000e+00 4.843907e-57 2.421954e-57 [270,] 1.000000e+00 7.110498e-58 3.555249e-58 [271,] 1.000000e+00 2.923686e-57 1.461843e-57 [272,] 1.000000e+00 1.257164e-56 6.285820e-57 [273,] 1.000000e+00 3.320846e-56 1.660423e-56 [274,] 1.000000e+00 1.400374e-55 7.001871e-56 [275,] 1.000000e+00 6.124195e-55 3.062097e-55 [276,] 1.000000e+00 2.525420e-54 1.262710e-54 [277,] 1.000000e+00 1.030206e-53 5.151031e-54 [278,] 1.000000e+00 4.396829e-53 2.198415e-53 [279,] 1.000000e+00 1.888888e-52 9.444438e-53 [280,] 1.000000e+00 7.166955e-52 3.583478e-52 [281,] 1.000000e+00 1.788863e-51 8.944316e-52 [282,] 1.000000e+00 2.728026e-51 1.364013e-51 [283,] 1.000000e+00 8.670489e-51 4.335244e-51 [284,] 1.000000e+00 3.460445e-50 1.730223e-50 [285,] 1.000000e+00 1.381566e-49 6.907831e-50 [286,] 1.000000e+00 5.581726e-49 2.790863e-49 [287,] 1.000000e+00 2.039335e-48 1.019667e-48 [288,] 1.000000e+00 7.066107e-49 3.533053e-49 [289,] 1.000000e+00 2.829218e-48 1.414609e-48 [290,] 1.000000e+00 1.108913e-47 5.544565e-48 [291,] 1.000000e+00 4.362960e-47 2.181480e-47 [292,] 1.000000e+00 1.830824e-46 9.154121e-47 [293,] 1.000000e+00 7.373535e-46 3.686768e-46 [294,] 1.000000e+00 3.045476e-45 1.522738e-45 [295,] 1.000000e+00 1.178187e-44 5.890934e-45 [296,] 1.000000e+00 4.188908e-44 2.094454e-44 [297,] 1.000000e+00 1.563332e-43 7.816660e-44 [298,] 1.000000e+00 4.707442e-49 2.353721e-49 [299,] 1.000000e+00 3.326489e-49 1.663245e-49 [300,] 1.000000e+00 7.723036e-49 3.861518e-49 [301,] 1.000000e+00 2.577760e-48 1.288880e-48 [302,] 1.000000e+00 7.213919e-48 3.606959e-48 [303,] 1.000000e+00 1.362487e-47 6.812433e-48 [304,] 1.000000e+00 3.972601e-47 1.986300e-47 [305,] 1.000000e+00 1.758436e-46 8.792182e-47 [306,] 1.000000e+00 7.817304e-46 3.908652e-46 [307,] 1.000000e+00 3.542915e-45 1.771457e-45 [308,] 1.000000e+00 1.575713e-44 7.878567e-45 [309,] 1.000000e+00 6.639463e-44 3.319731e-44 [310,] 1.000000e+00 2.854827e-43 1.427414e-43 [311,] 1.000000e+00 1.115517e-42 5.577583e-43 [312,] 1.000000e+00 2.845697e-42 1.422849e-42 [313,] 1.000000e+00 6.204047e-43 3.102024e-43 [314,] 1.000000e+00 2.009379e-42 1.004690e-42 [315,] 1.000000e+00 6.896150e-42 3.448075e-42 [316,] 1.000000e+00 3.037191e-41 1.518596e-41 [317,] 1.000000e+00 2.213316e-41 1.106658e-41 [318,] 1.000000e+00 9.803514e-41 4.901757e-41 [319,] 1.000000e+00 4.028475e-40 2.014238e-40 [320,] 1.000000e+00 1.606988e-39 8.034939e-40 [321,] 1.000000e+00 5.445244e-39 2.722622e-39 [322,] 1.000000e+00 1.390943e-38 6.954716e-39 [323,] 1.000000e+00 1.276686e-38 6.383429e-39 [324,] 1.000000e+00 5.576252e-38 2.788126e-38 [325,] 1.000000e+00 2.030890e-37 1.015445e-37 [326,] 1.000000e+00 7.847700e-37 3.923850e-37 [327,] 1.000000e+00 2.897574e-36 1.448787e-36 [328,] 1.000000e+00 1.210634e-35 6.053170e-36 [329,] 1.000000e+00 5.221051e-35 2.610526e-35 [330,] 1.000000e+00 1.252510e-34 6.262551e-35 [331,] 1.000000e+00 3.222515e-34 1.611257e-34 [332,] 1.000000e+00 9.476425e-34 4.738212e-34 [333,] 1.000000e+00 3.252099e-33 1.626050e-33 [334,] 1.000000e+00 1.076623e-32 5.383114e-33 [335,] 1.000000e+00 2.153740e-32 1.076870e-32 [336,] 1.000000e+00 3.324066e-32 1.662033e-32 [337,] 1.000000e+00 8.543189e-32 4.271595e-32 [338,] 1.000000e+00 8.374841e-32 4.187421e-32 [339,] 1.000000e+00 3.256933e-31 1.628467e-31 [340,] 1.000000e+00 1.380673e-30 6.903365e-31 [341,] 1.000000e+00 5.938681e-30 2.969340e-30 [342,] 1.000000e+00 2.569822e-29 1.284911e-29 [343,] 1.000000e+00 4.932369e-29 2.466185e-29 [344,] 1.000000e+00 1.848786e-28 9.243929e-29 [345,] 1.000000e+00 7.368034e-28 3.684017e-28 [346,] 1.000000e+00 3.371869e-28 1.685935e-28 [347,] 1.000000e+00 1.481362e-27 7.406809e-28 [348,] 1.000000e+00 5.574067e-27 2.787034e-27 [349,] 1.000000e+00 2.395027e-26 1.197514e-26 [350,] 1.000000e+00 9.030480e-26 4.515240e-26 [351,] 1.000000e+00 3.715089e-25 1.857544e-25 [352,] 1.000000e+00 1.542702e-24 7.713508e-25 [353,] 1.000000e+00 6.398811e-24 3.199406e-24 [354,] 1.000000e+00 2.053043e-23 1.026522e-23 [355,] 1.000000e+00 8.183285e-23 4.091642e-23 [356,] 1.000000e+00 2.798573e-22 1.399287e-22 [357,] 1.000000e+00 1.013929e-21 5.069645e-22 [358,] 1.000000e+00 2.454080e-21 1.227040e-21 [359,] 1.000000e+00 2.962882e-21 1.481441e-21 [360,] 1.000000e+00 1.179028e-20 5.895142e-21 [361,] 1.000000e+00 4.738108e-20 2.369054e-20 [362,] 1.000000e+00 1.457410e-19 7.287050e-20 [363,] 1.000000e+00 5.219423e-19 2.609711e-19 [364,] 1.000000e+00 9.780683e-19 4.890341e-19 [365,] 1.000000e+00 1.844550e-18 9.222749e-19 [366,] 1.000000e+00 4.261338e-18 2.130669e-18 [367,] 1.000000e+00 1.685929e-17 8.429644e-18 [368,] 1.000000e+00 5.188040e-17 2.594020e-17 [369,] 1.000000e+00 1.048218e-16 5.241088e-17 [370,] 1.000000e+00 3.114614e-16 1.557307e-16 [371,] 1.000000e+00 1.163363e-15 5.816813e-16 [372,] 1.000000e+00 4.379777e-15 2.189888e-15 [373,] 1.000000e+00 1.448556e-14 7.242778e-15 [374,] 1.000000e+00 4.869545e-14 2.434772e-14 [375,] 1.000000e+00 2.113912e-14 1.056956e-14 [376,] 1.000000e+00 8.240725e-14 4.120362e-14 [377,] 1.000000e+00 3.099388e-13 1.549694e-13 [378,] 1.000000e+00 5.640518e-13 2.820259e-13 [379,] 1.000000e+00 1.905713e-12 9.528567e-13 [380,] 1.000000e+00 5.993590e-12 2.996795e-12 [381,] 1.000000e+00 1.667702e-11 8.338510e-12 [382,] 1.000000e+00 5.902598e-11 2.951299e-11 [383,] 1.000000e+00 1.685058e-10 8.425292e-11 [384,] 1.000000e+00 1.698365e-16 8.491824e-17 [385,] 1.000000e+00 1.140675e-15 5.703373e-16 [386,] 1.000000e+00 5.280488e-15 2.640244e-15 [387,] 1.000000e+00 1.861105e-14 9.305527e-15 [388,] 1.000000e+00 7.073184e-15 3.536592e-15 [389,] 1.000000e+00 4.918897e-14 2.459448e-14 [390,] 1.000000e+00 3.456099e-13 1.728050e-13 [391,] 1.000000e+00 7.764342e-13 3.882171e-13 [392,] 1.000000e+00 2.633420e-13 1.316710e-13 [393,] 1.000000e+00 2.160003e-12 1.080001e-12 [394,] 1.000000e+00 9.199528e-13 4.599764e-13 [395,] 1.000000e+00 7.490755e-12 3.745378e-12 [396,] 1.000000e+00 7.194675e-11 3.597337e-11 [397,] 1.000000e+00 6.668261e-10 3.334131e-10 [398,] 1.000000e+00 1.536262e-09 7.681310e-10 [399,] 1.000000e+00 9.513533e-09 4.756767e-09 [400,] 1.000000e+00 6.340780e-08 3.170390e-08 [401,] 9.999997e-01 5.383670e-07 2.691835e-07 [402,] 9.999975e-01 4.948649e-06 2.474324e-06 [403,] 9.999875e-01 2.502287e-05 1.251144e-05 [404,] 9.999219e-01 1.562587e-04 7.812934e-05 [405,] 9.996912e-01 6.175382e-04 3.087691e-04 [406,] 9.997464e-01 5.072852e-04 2.536426e-04 > postscript(file="/var/www/html/rcomp/tmp/16our1291404004.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/26our1291404004.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/3zxbu1291404004.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/4zxbu1291404004.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/5zxbu1291404004.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 -41397.85532 -32152.42206 -38574.25939 -35499.16935 -35633.16404 6 7 8 9 10 -42584.11511 -41457.35194 -39906.37242 -39535.55271 -39804.79856 11 12 13 14 15 -35423.03451 -34502.02605 -43096.68105 -22325.97826 -25360.53912 16 17 18 19 20 -26099.04367 -24477.94971 -24161.42920 -24316.27212 -24630.91034 21 22 23 24 25 -25073.30631 -46895.85716 394402.77345 319400.77652 255397.65285 26 27 28 29 30 262324.05408 374855.31196 256585.49518 29642.20343 15333.78628 31 32 33 34 35 17534.71356 -40778.32474 -26711.92725 -38688.04034 -40927.70619 36 37 38 39 40 -38616.36570 -44422.33489 4741.88689 -33977.64053 -27585.55739 41 42 43 44 45 -40204.06214 -33484.61608 -43098.78070 -26820.54122 -25485.46813 46 47 48 49 50 -25102.03227 -45426.76260 -43438.03531 -41194.48366 -31411.76089 51 52 53 54 55 -43821.23451 -47796.01473 -33821.56951 -12523.27034 -16245.38865 56 57 58 59 60 -44017.99956 -49268.46308 -44542.60454 269927.37023 39851.36266 61 62 63 64 65 -5870.60226 133560.05198 -25564.85624 55726.15704 50359.69324 66 67 68 69 70 35251.26893 52395.97354 17537.94444 -1270.74574 -43039.25347 71 72 73 74 75 127183.62523 -43097.92398 -20961.19686 -23403.05424 -22130.61521 76 77 78 79 80 -24716.89345 -36983.80722 103106.87586 -20032.44129 -25116.36878 81 82 83 84 85 307975.33025 256713.28618 17534.63440 -39466.21178 -29649.54931 86 87 88 89 90 -42814.67540 -37405.96750 -43091.19730 -26818.97490 59557.86747 91 92 93 94 95 -43421.46506 -21890.47266 -22431.93445 -25084.91570 -46068.99932 96 97 98 99 100 -45830.28009 -47738.94494 -40295.75770 -41676.67164 -45589.19061 101 102 103 104 105 -39750.18103 -50876.66268 263550.49466 1578282.07571 234561.10141 106 107 108 109 110 259273.82647 296831.53219 145953.31889 256598.82544 26220.21231 111 112 113 114 115 17534.81127 -40841.31125 -34827.73757 -20469.06966 -23295.11129 116 117 118 119 120 -30192.19198 -44743.70650 -43092.68059 21520.04323 -16149.01550 121 122 123 124 125 -18840.82230 -25023.22218 -46373.11245 267297.67551 259459.16491 126 127 128 129 130 276129.02732 256952.66276 17535.20163 -43091.41039 -23760.30674 131 132 133 134 135 -25830.88212 11023.97240 -1790.67677 -42862.24687 8406.62905 136 137 138 139 140 -32635.48870 -17074.75935 -20840.87838 -51998.66225 7358.61218 141 142 143 144 145 -6469.59332 4784.50947 -1614.61640 -42856.76166 -19210.45639 146 147 148 149 150 -16884.58599 -20142.44158 -25765.65290 -21454.18416 -17664.24751 151 152 153 154 155 -21307.56737 -45680.30928 -47717.49487 -51902.30965 3374.53240 156 157 158 159 160 1763.21884 -1423.95701 -28478.96027 -40378.23179 -16968.67751 161 162 163 164 165 -9262.29181 -17092.77102 -3064.03770 9783.20902 -29439.09896 166 167 168 169 170 -26386.56629 -13746.69660 -9528.76052 -17052.98187 -19895.61527 171 172 173 174 175 -21533.35059 -16607.22138 -47669.86262 -49638.23526 7497.53171 176 177 178 179 180 -42855.27695 46019.30092 -13927.67182 -14301.84712 -22857.67343 181 182 183 184 185 -16999.40523 -8982.27827 -21287.94074 -21077.85285 -49967.90729 186 187 188 189 190 -51653.75996 53849.13278 -13686.71760 -568.86924 -42949.97113 191 192 193 194 195 36951.88741 -17000.47324 -21000.20437 -52846.21767 35864.33697 196 197 198 199 200 20118.71466 -2285.13861 -36570.40600 -41312.83589 -17705.34678 201 202 203 204 205 -38598.28489 -40110.45165 -40955.03596 -9682.80271 -38969.08454 206 207 208 209 210 -40220.50699 -40360.58797 -14752.28345 -4726.87087 -14249.55998 211 212 213 214 215 -16948.76746 -22877.23943 -22218.57668 -22276.40091 -9534.74369 216 217 218 219 220 -14394.14436 -37862.76577 -52413.02865 -40014.96792 -36110.44505 221 222 223 224 225 -41419.03955 -47364.78119 7605.33554 8668.13178 -5822.17609 226 227 228 229 230 -2005.61041 -42880.12317 -12420.58182 875.12824 -14834.04084 231 232 233 234 235 -17049.21581 -20683.26515 -28860.97728 -14097.61248 -21320.52641 236 237 238 239 240 -52250.44696 7503.82484 -519.29718 -38847.13662 -40375.85429 241 242 243 244 245 -16988.16168 -21443.65881 7506.65994 -42860.75121 26720.45761 246 247 248 249 250 -10636.06753 27266.74849 -48.35322 -21346.40406 -11896.90698 251 252 253 254 255 -26072.27190 -16953.87651 -22327.30894 -21124.00517 -51454.62861 256 257 258 259 260 -41500.82378 7496.04202 -40829.95487 -40277.42158 -38677.50926 261 262 263 264 265 -42711.55277 -40434.08000 -41617.33684 -41250.76601 -41145.04931 266 267 268 269 270 -40392.02544 -40297.11876 -40432.50499 -18730.82007 67319.62812 271 272 273 274 275 -15880.39341 -17156.18356 -10261.07682 25974.29917 -69824.25556 276 277 278 279 280 -20729.18186 -51869.15736 11189.92399 -1102.37815 -1565.24459 281 282 283 284 285 5302.40713 -35696.84169 -40366.36821 -17257.22533 48251.99514 286 287 288 289 290 -16350.57118 -16253.85723 -16354.97962 -17259.80889 -5990.68620 291 292 293 294 295 -7967.47203 -21490.05195 7499.65613 17139.70425 -38952.49878 296 297 298 299 300 -39678.92260 -40443.38125 -38991.37094 -38798.55385 54721.50152 301 302 303 304 305 -40949.73133 -40429.69329 -24224.58157 -9283.26261 -15524.91192 306 307 308 309 310 1148.75359 -17187.80310 -9579.56736 -20625.32555 113510.68869 311 312 313 314 315 82259.69522 -51949.74488 -51078.81527 -53513.58493 7495.09493 316 317 318 319 320 -9614.95895 -36998.86639 -40383.79695 -10685.54208 -8055.23967 321 322 323 324 325 -16999.76338 -10912.21228 -49574.95005 -52967.26681 37835.79712 326 327 328 329 330 -2953.04928 -39328.06989 -40162.45379 -47679.30493 -40345.85452 331 332 333 334 335 -17873.95821 -17087.38621 -21424.03173 -50170.40827 18506.16800 336 337 338 339 340 -1407.54517 -42884.32771 -17027.26850 -17358.91390 -42465.20732 341 342 343 344 345 -49191.32261 -52706.05608 -52399.00054 -50761.21556 -39586.04746 346 347 348 349 350 -43719.58465 -53058.85647 -51146.57851 14958.57357 52147.04213 351 352 353 354 355 -4012.18593 610.95649 18211.47877 899.49037 38058.13734 356 357 358 359 360 28425.29729 37424.53215 -3143.89891 -994.51897 -39707.04202 361 362 363 364 365 -39573.23068 -37838.50594 -39636.06178 -40458.38440 -40969.50230 366 367 368 369 370 -39449.80800 -40340.98116 -17035.73919 -21466.55339 7494.40813 371 372 373 374 375 -35040.90494 -40625.71509 -40309.54313 -17047.24691 -21439.99830 376 377 378 379 380 6479.23130 40559.77486 46349.40689 27079.60963 -8345.40199 381 382 383 384 385 64374.71861 -728.90531 -1507.52730 -1461.93879 25469.81413 386 387 388 389 390 37416.92062 67631.12268 314.21187 -1522.53307 -42907.33711 391 392 393 394 395 -27343.09452 -12813.88312 -16324.54732 10583.70460 -23348.28738 396 397 398 399 400 101553.74211 -16925.06079 -18429.04672 -19083.28176 163320.86876 401 402 403 404 405 -19815.25031 -23788.13724 -58835.60030 -21234.25979 -50323.72166 406 407 408 409 410 -50051.89743 3562.84132 -1388.63677 33089.10437 -3058.39376 411 412 413 414 415 15388.11982 72826.64188 179.52151 -2839.08546 -1553.85198 416 417 418 419 420 4026.25476 -3363.22032 38736.49310 37995.21619 -2675.89132 421 422 423 424 425 -42342.84306 -40138.29311 -40476.12267 -40845.27404 -40542.84185 426 427 428 429 430 -40358.05058 -40235.90951 -40418.50850 -19498.20932 -16882.14474 431 -2053.10198 > postscript(file="/var/www/html/rcomp/tmp/6aosf1291404004.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 -41397.85532 NA 1 -32152.42206 -41397.85532 2 -38574.25939 -32152.42206 3 -35499.16935 -38574.25939 4 -35633.16404 -35499.16935 5 -42584.11511 -35633.16404 6 -41457.35194 -42584.11511 7 -39906.37242 -41457.35194 8 -39535.55271 -39906.37242 9 -39804.79856 -39535.55271 10 -35423.03451 -39804.79856 11 -34502.02605 -35423.03451 12 -43096.68105 -34502.02605 13 -22325.97826 -43096.68105 14 -25360.53912 -22325.97826 15 -26099.04367 -25360.53912 16 -24477.94971 -26099.04367 17 -24161.42920 -24477.94971 18 -24316.27212 -24161.42920 19 -24630.91034 -24316.27212 20 -25073.30631 -24630.91034 21 -46895.85716 -25073.30631 22 394402.77345 -46895.85716 23 319400.77652 394402.77345 24 255397.65285 319400.77652 25 262324.05408 255397.65285 26 374855.31196 262324.05408 27 256585.49518 374855.31196 28 29642.20343 256585.49518 29 15333.78628 29642.20343 30 17534.71356 15333.78628 31 -40778.32474 17534.71356 32 -26711.92725 -40778.32474 33 -38688.04034 -26711.92725 34 -40927.70619 -38688.04034 35 -38616.36570 -40927.70619 36 -44422.33489 -38616.36570 37 4741.88689 -44422.33489 38 -33977.64053 4741.88689 39 -27585.55739 -33977.64053 40 -40204.06214 -27585.55739 41 -33484.61608 -40204.06214 42 -43098.78070 -33484.61608 43 -26820.54122 -43098.78070 44 -25485.46813 -26820.54122 45 -25102.03227 -25485.46813 46 -45426.76260 -25102.03227 47 -43438.03531 -45426.76260 48 -41194.48366 -43438.03531 49 -31411.76089 -41194.48366 50 -43821.23451 -31411.76089 51 -47796.01473 -43821.23451 52 -33821.56951 -47796.01473 53 -12523.27034 -33821.56951 54 -16245.38865 -12523.27034 55 -44017.99956 -16245.38865 56 -49268.46308 -44017.99956 57 -44542.60454 -49268.46308 58 269927.37023 -44542.60454 59 39851.36266 269927.37023 60 -5870.60226 39851.36266 61 133560.05198 -5870.60226 62 -25564.85624 133560.05198 63 55726.15704 -25564.85624 64 50359.69324 55726.15704 65 35251.26893 50359.69324 66 52395.97354 35251.26893 67 17537.94444 52395.97354 68 -1270.74574 17537.94444 69 -43039.25347 -1270.74574 70 127183.62523 -43039.25347 71 -43097.92398 127183.62523 72 -20961.19686 -43097.92398 73 -23403.05424 -20961.19686 74 -22130.61521 -23403.05424 75 -24716.89345 -22130.61521 76 -36983.80722 -24716.89345 77 103106.87586 -36983.80722 78 -20032.44129 103106.87586 79 -25116.36878 -20032.44129 80 307975.33025 -25116.36878 81 256713.28618 307975.33025 82 17534.63440 256713.28618 83 -39466.21178 17534.63440 84 -29649.54931 -39466.21178 85 -42814.67540 -29649.54931 86 -37405.96750 -42814.67540 87 -43091.19730 -37405.96750 88 -26818.97490 -43091.19730 89 59557.86747 -26818.97490 90 -43421.46506 59557.86747 91 -21890.47266 -43421.46506 92 -22431.93445 -21890.47266 93 -25084.91570 -22431.93445 94 -46068.99932 -25084.91570 95 -45830.28009 -46068.99932 96 -47738.94494 -45830.28009 97 -40295.75770 -47738.94494 98 -41676.67164 -40295.75770 99 -45589.19061 -41676.67164 100 -39750.18103 -45589.19061 101 -50876.66268 -39750.18103 102 263550.49466 -50876.66268 103 1578282.07571 263550.49466 104 234561.10141 1578282.07571 105 259273.82647 234561.10141 106 296831.53219 259273.82647 107 145953.31889 296831.53219 108 256598.82544 145953.31889 109 26220.21231 256598.82544 110 17534.81127 26220.21231 111 -40841.31125 17534.81127 112 -34827.73757 -40841.31125 113 -20469.06966 -34827.73757 114 -23295.11129 -20469.06966 115 -30192.19198 -23295.11129 116 -44743.70650 -30192.19198 117 -43092.68059 -44743.70650 118 21520.04323 -43092.68059 119 -16149.01550 21520.04323 120 -18840.82230 -16149.01550 121 -25023.22218 -18840.82230 122 -46373.11245 -25023.22218 123 267297.67551 -46373.11245 124 259459.16491 267297.67551 125 276129.02732 259459.16491 126 256952.66276 276129.02732 127 17535.20163 256952.66276 128 -43091.41039 17535.20163 129 -23760.30674 -43091.41039 130 -25830.88212 -23760.30674 131 11023.97240 -25830.88212 132 -1790.67677 11023.97240 133 -42862.24687 -1790.67677 134 8406.62905 -42862.24687 135 -32635.48870 8406.62905 136 -17074.75935 -32635.48870 137 -20840.87838 -17074.75935 138 -51998.66225 -20840.87838 139 7358.61218 -51998.66225 140 -6469.59332 7358.61218 141 4784.50947 -6469.59332 142 -1614.61640 4784.50947 143 -42856.76166 -1614.61640 144 -19210.45639 -42856.76166 145 -16884.58599 -19210.45639 146 -20142.44158 -16884.58599 147 -25765.65290 -20142.44158 148 -21454.18416 -25765.65290 149 -17664.24751 -21454.18416 150 -21307.56737 -17664.24751 151 -45680.30928 -21307.56737 152 -47717.49487 -45680.30928 153 -51902.30965 -47717.49487 154 3374.53240 -51902.30965 155 1763.21884 3374.53240 156 -1423.95701 1763.21884 157 -28478.96027 -1423.95701 158 -40378.23179 -28478.96027 159 -16968.67751 -40378.23179 160 -9262.29181 -16968.67751 161 -17092.77102 -9262.29181 162 -3064.03770 -17092.77102 163 9783.20902 -3064.03770 164 -29439.09896 9783.20902 165 -26386.56629 -29439.09896 166 -13746.69660 -26386.56629 167 -9528.76052 -13746.69660 168 -17052.98187 -9528.76052 169 -19895.61527 -17052.98187 170 -21533.35059 -19895.61527 171 -16607.22138 -21533.35059 172 -47669.86262 -16607.22138 173 -49638.23526 -47669.86262 174 7497.53171 -49638.23526 175 -42855.27695 7497.53171 176 46019.30092 -42855.27695 177 -13927.67182 46019.30092 178 -14301.84712 -13927.67182 179 -22857.67343 -14301.84712 180 -16999.40523 -22857.67343 181 -8982.27827 -16999.40523 182 -21287.94074 -8982.27827 183 -21077.85285 -21287.94074 184 -49967.90729 -21077.85285 185 -51653.75996 -49967.90729 186 53849.13278 -51653.75996 187 -13686.71760 53849.13278 188 -568.86924 -13686.71760 189 -42949.97113 -568.86924 190 36951.88741 -42949.97113 191 -17000.47324 36951.88741 192 -21000.20437 -17000.47324 193 -52846.21767 -21000.20437 194 35864.33697 -52846.21767 195 20118.71466 35864.33697 196 -2285.13861 20118.71466 197 -36570.40600 -2285.13861 198 -41312.83589 -36570.40600 199 -17705.34678 -41312.83589 200 -38598.28489 -17705.34678 201 -40110.45165 -38598.28489 202 -40955.03596 -40110.45165 203 -9682.80271 -40955.03596 204 -38969.08454 -9682.80271 205 -40220.50699 -38969.08454 206 -40360.58797 -40220.50699 207 -14752.28345 -40360.58797 208 -4726.87087 -14752.28345 209 -14249.55998 -4726.87087 210 -16948.76746 -14249.55998 211 -22877.23943 -16948.76746 212 -22218.57668 -22877.23943 213 -22276.40091 -22218.57668 214 -9534.74369 -22276.40091 215 -14394.14436 -9534.74369 216 -37862.76577 -14394.14436 217 -52413.02865 -37862.76577 218 -40014.96792 -52413.02865 219 -36110.44505 -40014.96792 220 -41419.03955 -36110.44505 221 -47364.78119 -41419.03955 222 7605.33554 -47364.78119 223 8668.13178 7605.33554 224 -5822.17609 8668.13178 225 -2005.61041 -5822.17609 226 -42880.12317 -2005.61041 227 -12420.58182 -42880.12317 228 875.12824 -12420.58182 229 -14834.04084 875.12824 230 -17049.21581 -14834.04084 231 -20683.26515 -17049.21581 232 -28860.97728 -20683.26515 233 -14097.61248 -28860.97728 234 -21320.52641 -14097.61248 235 -52250.44696 -21320.52641 236 7503.82484 -52250.44696 237 -519.29718 7503.82484 238 -38847.13662 -519.29718 239 -40375.85429 -38847.13662 240 -16988.16168 -40375.85429 241 -21443.65881 -16988.16168 242 7506.65994 -21443.65881 243 -42860.75121 7506.65994 244 26720.45761 -42860.75121 245 -10636.06753 26720.45761 246 27266.74849 -10636.06753 247 -48.35322 27266.74849 248 -21346.40406 -48.35322 249 -11896.90698 -21346.40406 250 -26072.27190 -11896.90698 251 -16953.87651 -26072.27190 252 -22327.30894 -16953.87651 253 -21124.00517 -22327.30894 254 -51454.62861 -21124.00517 255 -41500.82378 -51454.62861 256 7496.04202 -41500.82378 257 -40829.95487 7496.04202 258 -40277.42158 -40829.95487 259 -38677.50926 -40277.42158 260 -42711.55277 -38677.50926 261 -40434.08000 -42711.55277 262 -41617.33684 -40434.08000 263 -41250.76601 -41617.33684 264 -41145.04931 -41250.76601 265 -40392.02544 -41145.04931 266 -40297.11876 -40392.02544 267 -40432.50499 -40297.11876 268 -18730.82007 -40432.50499 269 67319.62812 -18730.82007 270 -15880.39341 67319.62812 271 -17156.18356 -15880.39341 272 -10261.07682 -17156.18356 273 25974.29917 -10261.07682 274 -69824.25556 25974.29917 275 -20729.18186 -69824.25556 276 -51869.15736 -20729.18186 277 11189.92399 -51869.15736 278 -1102.37815 11189.92399 279 -1565.24459 -1102.37815 280 5302.40713 -1565.24459 281 -35696.84169 5302.40713 282 -40366.36821 -35696.84169 283 -17257.22533 -40366.36821 284 48251.99514 -17257.22533 285 -16350.57118 48251.99514 286 -16253.85723 -16350.57118 287 -16354.97962 -16253.85723 288 -17259.80889 -16354.97962 289 -5990.68620 -17259.80889 290 -7967.47203 -5990.68620 291 -21490.05195 -7967.47203 292 7499.65613 -21490.05195 293 17139.70425 7499.65613 294 -38952.49878 17139.70425 295 -39678.92260 -38952.49878 296 -40443.38125 -39678.92260 297 -38991.37094 -40443.38125 298 -38798.55385 -38991.37094 299 54721.50152 -38798.55385 300 -40949.73133 54721.50152 301 -40429.69329 -40949.73133 302 -24224.58157 -40429.69329 303 -9283.26261 -24224.58157 304 -15524.91192 -9283.26261 305 1148.75359 -15524.91192 306 -17187.80310 1148.75359 307 -9579.56736 -17187.80310 308 -20625.32555 -9579.56736 309 113510.68869 -20625.32555 310 82259.69522 113510.68869 311 -51949.74488 82259.69522 312 -51078.81527 -51949.74488 313 -53513.58493 -51078.81527 314 7495.09493 -53513.58493 315 -9614.95895 7495.09493 316 -36998.86639 -9614.95895 317 -40383.79695 -36998.86639 318 -10685.54208 -40383.79695 319 -8055.23967 -10685.54208 320 -16999.76338 -8055.23967 321 -10912.21228 -16999.76338 322 -49574.95005 -10912.21228 323 -52967.26681 -49574.95005 324 37835.79712 -52967.26681 325 -2953.04928 37835.79712 326 -39328.06989 -2953.04928 327 -40162.45379 -39328.06989 328 -47679.30493 -40162.45379 329 -40345.85452 -47679.30493 330 -17873.95821 -40345.85452 331 -17087.38621 -17873.95821 332 -21424.03173 -17087.38621 333 -50170.40827 -21424.03173 334 18506.16800 -50170.40827 335 -1407.54517 18506.16800 336 -42884.32771 -1407.54517 337 -17027.26850 -42884.32771 338 -17358.91390 -17027.26850 339 -42465.20732 -17358.91390 340 -49191.32261 -42465.20732 341 -52706.05608 -49191.32261 342 -52399.00054 -52706.05608 343 -50761.21556 -52399.00054 344 -39586.04746 -50761.21556 345 -43719.58465 -39586.04746 346 -53058.85647 -43719.58465 347 -51146.57851 -53058.85647 348 14958.57357 -51146.57851 349 52147.04213 14958.57357 350 -4012.18593 52147.04213 351 610.95649 -4012.18593 352 18211.47877 610.95649 353 899.49037 18211.47877 354 38058.13734 899.49037 355 28425.29729 38058.13734 356 37424.53215 28425.29729 357 -3143.89891 37424.53215 358 -994.51897 -3143.89891 359 -39707.04202 -994.51897 360 -39573.23068 -39707.04202 361 -37838.50594 -39573.23068 362 -39636.06178 -37838.50594 363 -40458.38440 -39636.06178 364 -40969.50230 -40458.38440 365 -39449.80800 -40969.50230 366 -40340.98116 -39449.80800 367 -17035.73919 -40340.98116 368 -21466.55339 -17035.73919 369 7494.40813 -21466.55339 370 -35040.90494 7494.40813 371 -40625.71509 -35040.90494 372 -40309.54313 -40625.71509 373 -17047.24691 -40309.54313 374 -21439.99830 -17047.24691 375 6479.23130 -21439.99830 376 40559.77486 6479.23130 377 46349.40689 40559.77486 378 27079.60963 46349.40689 379 -8345.40199 27079.60963 380 64374.71861 -8345.40199 381 -728.90531 64374.71861 382 -1507.52730 -728.90531 383 -1461.93879 -1507.52730 384 25469.81413 -1461.93879 385 37416.92062 25469.81413 386 67631.12268 37416.92062 387 314.21187 67631.12268 388 -1522.53307 314.21187 389 -42907.33711 -1522.53307 390 -27343.09452 -42907.33711 391 -12813.88312 -27343.09452 392 -16324.54732 -12813.88312 393 10583.70460 -16324.54732 394 -23348.28738 10583.70460 395 101553.74211 -23348.28738 396 -16925.06079 101553.74211 397 -18429.04672 -16925.06079 398 -19083.28176 -18429.04672 399 163320.86876 -19083.28176 400 -19815.25031 163320.86876 401 -23788.13724 -19815.25031 402 -58835.60030 -23788.13724 403 -21234.25979 -58835.60030 404 -50323.72166 -21234.25979 405 -50051.89743 -50323.72166 406 3562.84132 -50051.89743 407 -1388.63677 3562.84132 408 33089.10437 -1388.63677 409 -3058.39376 33089.10437 410 15388.11982 -3058.39376 411 72826.64188 15388.11982 412 179.52151 72826.64188 413 -2839.08546 179.52151 414 -1553.85198 -2839.08546 415 4026.25476 -1553.85198 416 -3363.22032 4026.25476 417 38736.49310 -3363.22032 418 37995.21619 38736.49310 419 -2675.89132 37995.21619 420 -42342.84306 -2675.89132 421 -40138.29311 -42342.84306 422 -40476.12267 -40138.29311 423 -40845.27404 -40476.12267 424 -40542.84185 -40845.27404 425 -40358.05058 -40542.84185 426 -40235.90951 -40358.05058 427 -40418.50850 -40235.90951 428 -19498.20932 -40418.50850 429 -16882.14474 -19498.20932 430 -2053.10198 -16882.14474 431 NA -2053.10198 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -32152.42206 -41397.85532 [2,] -38574.25939 -32152.42206 [3,] -35499.16935 -38574.25939 [4,] -35633.16404 -35499.16935 [5,] -42584.11511 -35633.16404 [6,] -41457.35194 -42584.11511 [7,] -39906.37242 -41457.35194 [8,] -39535.55271 -39906.37242 [9,] -39804.79856 -39535.55271 [10,] -35423.03451 -39804.79856 [11,] -34502.02605 -35423.03451 [12,] -43096.68105 -34502.02605 [13,] -22325.97826 -43096.68105 [14,] -25360.53912 -22325.97826 [15,] -26099.04367 -25360.53912 [16,] -24477.94971 -26099.04367 [17,] -24161.42920 -24477.94971 [18,] -24316.27212 -24161.42920 [19,] -24630.91034 -24316.27212 [20,] -25073.30631 -24630.91034 [21,] -46895.85716 -25073.30631 [22,] 394402.77345 -46895.85716 [23,] 319400.77652 394402.77345 [24,] 255397.65285 319400.77652 [25,] 262324.05408 255397.65285 [26,] 374855.31196 262324.05408 [27,] 256585.49518 374855.31196 [28,] 29642.20343 256585.49518 [29,] 15333.78628 29642.20343 [30,] 17534.71356 15333.78628 [31,] -40778.32474 17534.71356 [32,] -26711.92725 -40778.32474 [33,] -38688.04034 -26711.92725 [34,] -40927.70619 -38688.04034 [35,] -38616.36570 -40927.70619 [36,] -44422.33489 -38616.36570 [37,] 4741.88689 -44422.33489 [38,] -33977.64053 4741.88689 [39,] -27585.55739 -33977.64053 [40,] -40204.06214 -27585.55739 [41,] -33484.61608 -40204.06214 [42,] -43098.78070 -33484.61608 [43,] -26820.54122 -43098.78070 [44,] -25485.46813 -26820.54122 [45,] -25102.03227 -25485.46813 [46,] -45426.76260 -25102.03227 [47,] -43438.03531 -45426.76260 [48,] -41194.48366 -43438.03531 [49,] -31411.76089 -41194.48366 [50,] -43821.23451 -31411.76089 [51,] -47796.01473 -43821.23451 [52,] -33821.56951 -47796.01473 [53,] -12523.27034 -33821.56951 [54,] -16245.38865 -12523.27034 [55,] -44017.99956 -16245.38865 [56,] -49268.46308 -44017.99956 [57,] -44542.60454 -49268.46308 [58,] 269927.37023 -44542.60454 [59,] 39851.36266 269927.37023 [60,] -5870.60226 39851.36266 [61,] 133560.05198 -5870.60226 [62,] -25564.85624 133560.05198 [63,] 55726.15704 -25564.85624 [64,] 50359.69324 55726.15704 [65,] 35251.26893 50359.69324 [66,] 52395.97354 35251.26893 [67,] 17537.94444 52395.97354 [68,] -1270.74574 17537.94444 [69,] -43039.25347 -1270.74574 [70,] 127183.62523 -43039.25347 [71,] -43097.92398 127183.62523 [72,] -20961.19686 -43097.92398 [73,] -23403.05424 -20961.19686 [74,] -22130.61521 -23403.05424 [75,] -24716.89345 -22130.61521 [76,] -36983.80722 -24716.89345 [77,] 103106.87586 -36983.80722 [78,] -20032.44129 103106.87586 [79,] -25116.36878 -20032.44129 [80,] 307975.33025 -25116.36878 [81,] 256713.28618 307975.33025 [82,] 17534.63440 256713.28618 [83,] -39466.21178 17534.63440 [84,] -29649.54931 -39466.21178 [85,] -42814.67540 -29649.54931 [86,] -37405.96750 -42814.67540 [87,] -43091.19730 -37405.96750 [88,] -26818.97490 -43091.19730 [89,] 59557.86747 -26818.97490 [90,] -43421.46506 59557.86747 [91,] -21890.47266 -43421.46506 [92,] -22431.93445 -21890.47266 [93,] -25084.91570 -22431.93445 [94,] -46068.99932 -25084.91570 [95,] -45830.28009 -46068.99932 [96,] -47738.94494 -45830.28009 [97,] -40295.75770 -47738.94494 [98,] -41676.67164 -40295.75770 [99,] -45589.19061 -41676.67164 [100,] -39750.18103 -45589.19061 [101,] -50876.66268 -39750.18103 [102,] 263550.49466 -50876.66268 [103,] 1578282.07571 263550.49466 [104,] 234561.10141 1578282.07571 [105,] 259273.82647 234561.10141 [106,] 296831.53219 259273.82647 [107,] 145953.31889 296831.53219 [108,] 256598.82544 145953.31889 [109,] 26220.21231 256598.82544 [110,] 17534.81127 26220.21231 [111,] -40841.31125 17534.81127 [112,] -34827.73757 -40841.31125 [113,] -20469.06966 -34827.73757 [114,] -23295.11129 -20469.06966 [115,] -30192.19198 -23295.11129 [116,] -44743.70650 -30192.19198 [117,] -43092.68059 -44743.70650 [118,] 21520.04323 -43092.68059 [119,] -16149.01550 21520.04323 [120,] -18840.82230 -16149.01550 [121,] -25023.22218 -18840.82230 [122,] -46373.11245 -25023.22218 [123,] 267297.67551 -46373.11245 [124,] 259459.16491 267297.67551 [125,] 276129.02732 259459.16491 [126,] 256952.66276 276129.02732 [127,] 17535.20163 256952.66276 [128,] -43091.41039 17535.20163 [129,] -23760.30674 -43091.41039 [130,] -25830.88212 -23760.30674 [131,] 11023.97240 -25830.88212 [132,] -1790.67677 11023.97240 [133,] -42862.24687 -1790.67677 [134,] 8406.62905 -42862.24687 [135,] -32635.48870 8406.62905 [136,] -17074.75935 -32635.48870 [137,] -20840.87838 -17074.75935 [138,] -51998.66225 -20840.87838 [139,] 7358.61218 -51998.66225 [140,] -6469.59332 7358.61218 [141,] 4784.50947 -6469.59332 [142,] -1614.61640 4784.50947 [143,] -42856.76166 -1614.61640 [144,] -19210.45639 -42856.76166 [145,] -16884.58599 -19210.45639 [146,] -20142.44158 -16884.58599 [147,] -25765.65290 -20142.44158 [148,] -21454.18416 -25765.65290 [149,] -17664.24751 -21454.18416 [150,] -21307.56737 -17664.24751 [151,] -45680.30928 -21307.56737 [152,] -47717.49487 -45680.30928 [153,] -51902.30965 -47717.49487 [154,] 3374.53240 -51902.30965 [155,] 1763.21884 3374.53240 [156,] -1423.95701 1763.21884 [157,] -28478.96027 -1423.95701 [158,] -40378.23179 -28478.96027 [159,] -16968.67751 -40378.23179 [160,] -9262.29181 -16968.67751 [161,] -17092.77102 -9262.29181 [162,] -3064.03770 -17092.77102 [163,] 9783.20902 -3064.03770 [164,] -29439.09896 9783.20902 [165,] -26386.56629 -29439.09896 [166,] -13746.69660 -26386.56629 [167,] -9528.76052 -13746.69660 [168,] -17052.98187 -9528.76052 [169,] -19895.61527 -17052.98187 [170,] -21533.35059 -19895.61527 [171,] -16607.22138 -21533.35059 [172,] -47669.86262 -16607.22138 [173,] -49638.23526 -47669.86262 [174,] 7497.53171 -49638.23526 [175,] -42855.27695 7497.53171 [176,] 46019.30092 -42855.27695 [177,] -13927.67182 46019.30092 [178,] -14301.84712 -13927.67182 [179,] -22857.67343 -14301.84712 [180,] -16999.40523 -22857.67343 [181,] -8982.27827 -16999.40523 [182,] -21287.94074 -8982.27827 [183,] -21077.85285 -21287.94074 [184,] -49967.90729 -21077.85285 [185,] -51653.75996 -49967.90729 [186,] 53849.13278 -51653.75996 [187,] -13686.71760 53849.13278 [188,] -568.86924 -13686.71760 [189,] -42949.97113 -568.86924 [190,] 36951.88741 -42949.97113 [191,] -17000.47324 36951.88741 [192,] -21000.20437 -17000.47324 [193,] -52846.21767 -21000.20437 [194,] 35864.33697 -52846.21767 [195,] 20118.71466 35864.33697 [196,] -2285.13861 20118.71466 [197,] -36570.40600 -2285.13861 [198,] -41312.83589 -36570.40600 [199,] -17705.34678 -41312.83589 [200,] -38598.28489 -17705.34678 [201,] -40110.45165 -38598.28489 [202,] -40955.03596 -40110.45165 [203,] -9682.80271 -40955.03596 [204,] -38969.08454 -9682.80271 [205,] -40220.50699 -38969.08454 [206,] -40360.58797 -40220.50699 [207,] -14752.28345 -40360.58797 [208,] -4726.87087 -14752.28345 [209,] -14249.55998 -4726.87087 [210,] -16948.76746 -14249.55998 [211,] -22877.23943 -16948.76746 [212,] -22218.57668 -22877.23943 [213,] -22276.40091 -22218.57668 [214,] -9534.74369 -22276.40091 [215,] -14394.14436 -9534.74369 [216,] -37862.76577 -14394.14436 [217,] -52413.02865 -37862.76577 [218,] -40014.96792 -52413.02865 [219,] -36110.44505 -40014.96792 [220,] -41419.03955 -36110.44505 [221,] -47364.78119 -41419.03955 [222,] 7605.33554 -47364.78119 [223,] 8668.13178 7605.33554 [224,] -5822.17609 8668.13178 [225,] -2005.61041 -5822.17609 [226,] -42880.12317 -2005.61041 [227,] -12420.58182 -42880.12317 [228,] 875.12824 -12420.58182 [229,] -14834.04084 875.12824 [230,] -17049.21581 -14834.04084 [231,] -20683.26515 -17049.21581 [232,] -28860.97728 -20683.26515 [233,] -14097.61248 -28860.97728 [234,] -21320.52641 -14097.61248 [235,] -52250.44696 -21320.52641 [236,] 7503.82484 -52250.44696 [237,] -519.29718 7503.82484 [238,] -38847.13662 -519.29718 [239,] -40375.85429 -38847.13662 [240,] -16988.16168 -40375.85429 [241,] -21443.65881 -16988.16168 [242,] 7506.65994 -21443.65881 [243,] -42860.75121 7506.65994 [244,] 26720.45761 -42860.75121 [245,] -10636.06753 26720.45761 [246,] 27266.74849 -10636.06753 [247,] -48.35322 27266.74849 [248,] -21346.40406 -48.35322 [249,] -11896.90698 -21346.40406 [250,] -26072.27190 -11896.90698 [251,] -16953.87651 -26072.27190 [252,] -22327.30894 -16953.87651 [253,] -21124.00517 -22327.30894 [254,] -51454.62861 -21124.00517 [255,] -41500.82378 -51454.62861 [256,] 7496.04202 -41500.82378 [257,] -40829.95487 7496.04202 [258,] -40277.42158 -40829.95487 [259,] -38677.50926 -40277.42158 [260,] -42711.55277 -38677.50926 [261,] -40434.08000 -42711.55277 [262,] -41617.33684 -40434.08000 [263,] -41250.76601 -41617.33684 [264,] -41145.04931 -41250.76601 [265,] -40392.02544 -41145.04931 [266,] -40297.11876 -40392.02544 [267,] -40432.50499 -40297.11876 [268,] -18730.82007 -40432.50499 [269,] 67319.62812 -18730.82007 [270,] -15880.39341 67319.62812 [271,] -17156.18356 -15880.39341 [272,] -10261.07682 -17156.18356 [273,] 25974.29917 -10261.07682 [274,] -69824.25556 25974.29917 [275,] -20729.18186 -69824.25556 [276,] -51869.15736 -20729.18186 [277,] 11189.92399 -51869.15736 [278,] -1102.37815 11189.92399 [279,] -1565.24459 -1102.37815 [280,] 5302.40713 -1565.24459 [281,] -35696.84169 5302.40713 [282,] -40366.36821 -35696.84169 [283,] -17257.22533 -40366.36821 [284,] 48251.99514 -17257.22533 [285,] -16350.57118 48251.99514 [286,] -16253.85723 -16350.57118 [287,] -16354.97962 -16253.85723 [288,] -17259.80889 -16354.97962 [289,] -5990.68620 -17259.80889 [290,] -7967.47203 -5990.68620 [291,] -21490.05195 -7967.47203 [292,] 7499.65613 -21490.05195 [293,] 17139.70425 7499.65613 [294,] -38952.49878 17139.70425 [295,] -39678.92260 -38952.49878 [296,] -40443.38125 -39678.92260 [297,] -38991.37094 -40443.38125 [298,] -38798.55385 -38991.37094 [299,] 54721.50152 -38798.55385 [300,] -40949.73133 54721.50152 [301,] -40429.69329 -40949.73133 [302,] -24224.58157 -40429.69329 [303,] -9283.26261 -24224.58157 [304,] -15524.91192 -9283.26261 [305,] 1148.75359 -15524.91192 [306,] -17187.80310 1148.75359 [307,] -9579.56736 -17187.80310 [308,] -20625.32555 -9579.56736 [309,] 113510.68869 -20625.32555 [310,] 82259.69522 113510.68869 [311,] -51949.74488 82259.69522 [312,] -51078.81527 -51949.74488 [313,] -53513.58493 -51078.81527 [314,] 7495.09493 -53513.58493 [315,] -9614.95895 7495.09493 [316,] -36998.86639 -9614.95895 [317,] -40383.79695 -36998.86639 [318,] -10685.54208 -40383.79695 [319,] -8055.23967 -10685.54208 [320,] -16999.76338 -8055.23967 [321,] -10912.21228 -16999.76338 [322,] -49574.95005 -10912.21228 [323,] -52967.26681 -49574.95005 [324,] 37835.79712 -52967.26681 [325,] -2953.04928 37835.79712 [326,] -39328.06989 -2953.04928 [327,] -40162.45379 -39328.06989 [328,] -47679.30493 -40162.45379 [329,] -40345.85452 -47679.30493 [330,] -17873.95821 -40345.85452 [331,] -17087.38621 -17873.95821 [332,] -21424.03173 -17087.38621 [333,] -50170.40827 -21424.03173 [334,] 18506.16800 -50170.40827 [335,] -1407.54517 18506.16800 [336,] -42884.32771 -1407.54517 [337,] -17027.26850 -42884.32771 [338,] -17358.91390 -17027.26850 [339,] -42465.20732 -17358.91390 [340,] -49191.32261 -42465.20732 [341,] -52706.05608 -49191.32261 [342,] -52399.00054 -52706.05608 [343,] -50761.21556 -52399.00054 [344,] -39586.04746 -50761.21556 [345,] -43719.58465 -39586.04746 [346,] -53058.85647 -43719.58465 [347,] -51146.57851 -53058.85647 [348,] 14958.57357 -51146.57851 [349,] 52147.04213 14958.57357 [350,] -4012.18593 52147.04213 [351,] 610.95649 -4012.18593 [352,] 18211.47877 610.95649 [353,] 899.49037 18211.47877 [354,] 38058.13734 899.49037 [355,] 28425.29729 38058.13734 [356,] 37424.53215 28425.29729 [357,] -3143.89891 37424.53215 [358,] -994.51897 -3143.89891 [359,] -39707.04202 -994.51897 [360,] -39573.23068 -39707.04202 [361,] -37838.50594 -39573.23068 [362,] -39636.06178 -37838.50594 [363,] -40458.38440 -39636.06178 [364,] -40969.50230 -40458.38440 [365,] -39449.80800 -40969.50230 [366,] -40340.98116 -39449.80800 [367,] -17035.73919 -40340.98116 [368,] -21466.55339 -17035.73919 [369,] 7494.40813 -21466.55339 [370,] -35040.90494 7494.40813 [371,] -40625.71509 -35040.90494 [372,] -40309.54313 -40625.71509 [373,] -17047.24691 -40309.54313 [374,] -21439.99830 -17047.24691 [375,] 6479.23130 -21439.99830 [376,] 40559.77486 6479.23130 [377,] 46349.40689 40559.77486 [378,] 27079.60963 46349.40689 [379,] -8345.40199 27079.60963 [380,] 64374.71861 -8345.40199 [381,] -728.90531 64374.71861 [382,] -1507.52730 -728.90531 [383,] -1461.93879 -1507.52730 [384,] 25469.81413 -1461.93879 [385,] 37416.92062 25469.81413 [386,] 67631.12268 37416.92062 [387,] 314.21187 67631.12268 [388,] -1522.53307 314.21187 [389,] -42907.33711 -1522.53307 [390,] -27343.09452 -42907.33711 [391,] -12813.88312 -27343.09452 [392,] -16324.54732 -12813.88312 [393,] 10583.70460 -16324.54732 [394,] -23348.28738 10583.70460 [395,] 101553.74211 -23348.28738 [396,] -16925.06079 101553.74211 [397,] -18429.04672 -16925.06079 [398,] -19083.28176 -18429.04672 [399,] 163320.86876 -19083.28176 [400,] -19815.25031 163320.86876 [401,] -23788.13724 -19815.25031 [402,] -58835.60030 -23788.13724 [403,] -21234.25979 -58835.60030 [404,] -50323.72166 -21234.25979 [405,] -50051.89743 -50323.72166 [406,] 3562.84132 -50051.89743 [407,] -1388.63677 3562.84132 [408,] 33089.10437 -1388.63677 [409,] -3058.39376 33089.10437 [410,] 15388.11982 -3058.39376 [411,] 72826.64188 15388.11982 [412,] 179.52151 72826.64188 [413,] -2839.08546 179.52151 [414,] -1553.85198 -2839.08546 [415,] 4026.25476 -1553.85198 [416,] -3363.22032 4026.25476 [417,] 38736.49310 -3363.22032 [418,] 37995.21619 38736.49310 [419,] -2675.89132 37995.21619 [420,] -42342.84306 -2675.89132 [421,] -40138.29311 -42342.84306 [422,] -40476.12267 -40138.29311 [423,] -40845.27404 -40476.12267 [424,] -40542.84185 -40845.27404 [425,] -40358.05058 -40542.84185 [426,] -40235.90951 -40358.05058 [427,] -40418.50850 -40235.90951 [428,] -19498.20932 -40418.50850 [429,] -16882.14474 -19498.20932 [430,] -2053.10198 -16882.14474 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -32152.42206 -41397.85532 2 -38574.25939 -32152.42206 3 -35499.16935 -38574.25939 4 -35633.16404 -35499.16935 5 -42584.11511 -35633.16404 6 -41457.35194 -42584.11511 7 -39906.37242 -41457.35194 8 -39535.55271 -39906.37242 9 -39804.79856 -39535.55271 10 -35423.03451 -39804.79856 11 -34502.02605 -35423.03451 12 -43096.68105 -34502.02605 13 -22325.97826 -43096.68105 14 -25360.53912 -22325.97826 15 -26099.04367 -25360.53912 16 -24477.94971 -26099.04367 17 -24161.42920 -24477.94971 18 -24316.27212 -24161.42920 19 -24630.91034 -24316.27212 20 -25073.30631 -24630.91034 21 -46895.85716 -25073.30631 22 394402.77345 -46895.85716 23 319400.77652 394402.77345 24 255397.65285 319400.77652 25 262324.05408 255397.65285 26 374855.31196 262324.05408 27 256585.49518 374855.31196 28 29642.20343 256585.49518 29 15333.78628 29642.20343 30 17534.71356 15333.78628 31 -40778.32474 17534.71356 32 -26711.92725 -40778.32474 33 -38688.04034 -26711.92725 34 -40927.70619 -38688.04034 35 -38616.36570 -40927.70619 36 -44422.33489 -38616.36570 37 4741.88689 -44422.33489 38 -33977.64053 4741.88689 39 -27585.55739 -33977.64053 40 -40204.06214 -27585.55739 41 -33484.61608 -40204.06214 42 -43098.78070 -33484.61608 43 -26820.54122 -43098.78070 44 -25485.46813 -26820.54122 45 -25102.03227 -25485.46813 46 -45426.76260 -25102.03227 47 -43438.03531 -45426.76260 48 -41194.48366 -43438.03531 49 -31411.76089 -41194.48366 50 -43821.23451 -31411.76089 51 -47796.01473 -43821.23451 52 -33821.56951 -47796.01473 53 -12523.27034 -33821.56951 54 -16245.38865 -12523.27034 55 -44017.99956 -16245.38865 56 -49268.46308 -44017.99956 57 -44542.60454 -49268.46308 58 269927.37023 -44542.60454 59 39851.36266 269927.37023 60 -5870.60226 39851.36266 61 133560.05198 -5870.60226 62 -25564.85624 133560.05198 63 55726.15704 -25564.85624 64 50359.69324 55726.15704 65 35251.26893 50359.69324 66 52395.97354 35251.26893 67 17537.94444 52395.97354 68 -1270.74574 17537.94444 69 -43039.25347 -1270.74574 70 127183.62523 -43039.25347 71 -43097.92398 127183.62523 72 -20961.19686 -43097.92398 73 -23403.05424 -20961.19686 74 -22130.61521 -23403.05424 75 -24716.89345 -22130.61521 76 -36983.80722 -24716.89345 77 103106.87586 -36983.80722 78 -20032.44129 103106.87586 79 -25116.36878 -20032.44129 80 307975.33025 -25116.36878 81 256713.28618 307975.33025 82 17534.63440 256713.28618 83 -39466.21178 17534.63440 84 -29649.54931 -39466.21178 85 -42814.67540 -29649.54931 86 -37405.96750 -42814.67540 87 -43091.19730 -37405.96750 88 -26818.97490 -43091.19730 89 59557.86747 -26818.97490 90 -43421.46506 59557.86747 91 -21890.47266 -43421.46506 92 -22431.93445 -21890.47266 93 -25084.91570 -22431.93445 94 -46068.99932 -25084.91570 95 -45830.28009 -46068.99932 96 -47738.94494 -45830.28009 97 -40295.75770 -47738.94494 98 -41676.67164 -40295.75770 99 -45589.19061 -41676.67164 100 -39750.18103 -45589.19061 101 -50876.66268 -39750.18103 102 263550.49466 -50876.66268 103 1578282.07571 263550.49466 104 234561.10141 1578282.07571 105 259273.82647 234561.10141 106 296831.53219 259273.82647 107 145953.31889 296831.53219 108 256598.82544 145953.31889 109 26220.21231 256598.82544 110 17534.81127 26220.21231 111 -40841.31125 17534.81127 112 -34827.73757 -40841.31125 113 -20469.06966 -34827.73757 114 -23295.11129 -20469.06966 115 -30192.19198 -23295.11129 116 -44743.70650 -30192.19198 117 -43092.68059 -44743.70650 118 21520.04323 -43092.68059 119 -16149.01550 21520.04323 120 -18840.82230 -16149.01550 121 -25023.22218 -18840.82230 122 -46373.11245 -25023.22218 123 267297.67551 -46373.11245 124 259459.16491 267297.67551 125 276129.02732 259459.16491 126 256952.66276 276129.02732 127 17535.20163 256952.66276 128 -43091.41039 17535.20163 129 -23760.30674 -43091.41039 130 -25830.88212 -23760.30674 131 11023.97240 -25830.88212 132 -1790.67677 11023.97240 133 -42862.24687 -1790.67677 134 8406.62905 -42862.24687 135 -32635.48870 8406.62905 136 -17074.75935 -32635.48870 137 -20840.87838 -17074.75935 138 -51998.66225 -20840.87838 139 7358.61218 -51998.66225 140 -6469.59332 7358.61218 141 4784.50947 -6469.59332 142 -1614.61640 4784.50947 143 -42856.76166 -1614.61640 144 -19210.45639 -42856.76166 145 -16884.58599 -19210.45639 146 -20142.44158 -16884.58599 147 -25765.65290 -20142.44158 148 -21454.18416 -25765.65290 149 -17664.24751 -21454.18416 150 -21307.56737 -17664.24751 151 -45680.30928 -21307.56737 152 -47717.49487 -45680.30928 153 -51902.30965 -47717.49487 154 3374.53240 -51902.30965 155 1763.21884 3374.53240 156 -1423.95701 1763.21884 157 -28478.96027 -1423.95701 158 -40378.23179 -28478.96027 159 -16968.67751 -40378.23179 160 -9262.29181 -16968.67751 161 -17092.77102 -9262.29181 162 -3064.03770 -17092.77102 163 9783.20902 -3064.03770 164 -29439.09896 9783.20902 165 -26386.56629 -29439.09896 166 -13746.69660 -26386.56629 167 -9528.76052 -13746.69660 168 -17052.98187 -9528.76052 169 -19895.61527 -17052.98187 170 -21533.35059 -19895.61527 171 -16607.22138 -21533.35059 172 -47669.86262 -16607.22138 173 -49638.23526 -47669.86262 174 7497.53171 -49638.23526 175 -42855.27695 7497.53171 176 46019.30092 -42855.27695 177 -13927.67182 46019.30092 178 -14301.84712 -13927.67182 179 -22857.67343 -14301.84712 180 -16999.40523 -22857.67343 181 -8982.27827 -16999.40523 182 -21287.94074 -8982.27827 183 -21077.85285 -21287.94074 184 -49967.90729 -21077.85285 185 -51653.75996 -49967.90729 186 53849.13278 -51653.75996 187 -13686.71760 53849.13278 188 -568.86924 -13686.71760 189 -42949.97113 -568.86924 190 36951.88741 -42949.97113 191 -17000.47324 36951.88741 192 -21000.20437 -17000.47324 193 -52846.21767 -21000.20437 194 35864.33697 -52846.21767 195 20118.71466 35864.33697 196 -2285.13861 20118.71466 197 -36570.40600 -2285.13861 198 -41312.83589 -36570.40600 199 -17705.34678 -41312.83589 200 -38598.28489 -17705.34678 201 -40110.45165 -38598.28489 202 -40955.03596 -40110.45165 203 -9682.80271 -40955.03596 204 -38969.08454 -9682.80271 205 -40220.50699 -38969.08454 206 -40360.58797 -40220.50699 207 -14752.28345 -40360.58797 208 -4726.87087 -14752.28345 209 -14249.55998 -4726.87087 210 -16948.76746 -14249.55998 211 -22877.23943 -16948.76746 212 -22218.57668 -22877.23943 213 -22276.40091 -22218.57668 214 -9534.74369 -22276.40091 215 -14394.14436 -9534.74369 216 -37862.76577 -14394.14436 217 -52413.02865 -37862.76577 218 -40014.96792 -52413.02865 219 -36110.44505 -40014.96792 220 -41419.03955 -36110.44505 221 -47364.78119 -41419.03955 222 7605.33554 -47364.78119 223 8668.13178 7605.33554 224 -5822.17609 8668.13178 225 -2005.61041 -5822.17609 226 -42880.12317 -2005.61041 227 -12420.58182 -42880.12317 228 875.12824 -12420.58182 229 -14834.04084 875.12824 230 -17049.21581 -14834.04084 231 -20683.26515 -17049.21581 232 -28860.97728 -20683.26515 233 -14097.61248 -28860.97728 234 -21320.52641 -14097.61248 235 -52250.44696 -21320.52641 236 7503.82484 -52250.44696 237 -519.29718 7503.82484 238 -38847.13662 -519.29718 239 -40375.85429 -38847.13662 240 -16988.16168 -40375.85429 241 -21443.65881 -16988.16168 242 7506.65994 -21443.65881 243 -42860.75121 7506.65994 244 26720.45761 -42860.75121 245 -10636.06753 26720.45761 246 27266.74849 -10636.06753 247 -48.35322 27266.74849 248 -21346.40406 -48.35322 249 -11896.90698 -21346.40406 250 -26072.27190 -11896.90698 251 -16953.87651 -26072.27190 252 -22327.30894 -16953.87651 253 -21124.00517 -22327.30894 254 -51454.62861 -21124.00517 255 -41500.82378 -51454.62861 256 7496.04202 -41500.82378 257 -40829.95487 7496.04202 258 -40277.42158 -40829.95487 259 -38677.50926 -40277.42158 260 -42711.55277 -38677.50926 261 -40434.08000 -42711.55277 262 -41617.33684 -40434.08000 263 -41250.76601 -41617.33684 264 -41145.04931 -41250.76601 265 -40392.02544 -41145.04931 266 -40297.11876 -40392.02544 267 -40432.50499 -40297.11876 268 -18730.82007 -40432.50499 269 67319.62812 -18730.82007 270 -15880.39341 67319.62812 271 -17156.18356 -15880.39341 272 -10261.07682 -17156.18356 273 25974.29917 -10261.07682 274 -69824.25556 25974.29917 275 -20729.18186 -69824.25556 276 -51869.15736 -20729.18186 277 11189.92399 -51869.15736 278 -1102.37815 11189.92399 279 -1565.24459 -1102.37815 280 5302.40713 -1565.24459 281 -35696.84169 5302.40713 282 -40366.36821 -35696.84169 283 -17257.22533 -40366.36821 284 48251.99514 -17257.22533 285 -16350.57118 48251.99514 286 -16253.85723 -16350.57118 287 -16354.97962 -16253.85723 288 -17259.80889 -16354.97962 289 -5990.68620 -17259.80889 290 -7967.47203 -5990.68620 291 -21490.05195 -7967.47203 292 7499.65613 -21490.05195 293 17139.70425 7499.65613 294 -38952.49878 17139.70425 295 -39678.92260 -38952.49878 296 -40443.38125 -39678.92260 297 -38991.37094 -40443.38125 298 -38798.55385 -38991.37094 299 54721.50152 -38798.55385 300 -40949.73133 54721.50152 301 -40429.69329 -40949.73133 302 -24224.58157 -40429.69329 303 -9283.26261 -24224.58157 304 -15524.91192 -9283.26261 305 1148.75359 -15524.91192 306 -17187.80310 1148.75359 307 -9579.56736 -17187.80310 308 -20625.32555 -9579.56736 309 113510.68869 -20625.32555 310 82259.69522 113510.68869 311 -51949.74488 82259.69522 312 -51078.81527 -51949.74488 313 -53513.58493 -51078.81527 314 7495.09493 -53513.58493 315 -9614.95895 7495.09493 316 -36998.86639 -9614.95895 317 -40383.79695 -36998.86639 318 -10685.54208 -40383.79695 319 -8055.23967 -10685.54208 320 -16999.76338 -8055.23967 321 -10912.21228 -16999.76338 322 -49574.95005 -10912.21228 323 -52967.26681 -49574.95005 324 37835.79712 -52967.26681 325 -2953.04928 37835.79712 326 -39328.06989 -2953.04928 327 -40162.45379 -39328.06989 328 -47679.30493 -40162.45379 329 -40345.85452 -47679.30493 330 -17873.95821 -40345.85452 331 -17087.38621 -17873.95821 332 -21424.03173 -17087.38621 333 -50170.40827 -21424.03173 334 18506.16800 -50170.40827 335 -1407.54517 18506.16800 336 -42884.32771 -1407.54517 337 -17027.26850 -42884.32771 338 -17358.91390 -17027.26850 339 -42465.20732 -17358.91390 340 -49191.32261 -42465.20732 341 -52706.05608 -49191.32261 342 -52399.00054 -52706.05608 343 -50761.21556 -52399.00054 344 -39586.04746 -50761.21556 345 -43719.58465 -39586.04746 346 -53058.85647 -43719.58465 347 -51146.57851 -53058.85647 348 14958.57357 -51146.57851 349 52147.04213 14958.57357 350 -4012.18593 52147.04213 351 610.95649 -4012.18593 352 18211.47877 610.95649 353 899.49037 18211.47877 354 38058.13734 899.49037 355 28425.29729 38058.13734 356 37424.53215 28425.29729 357 -3143.89891 37424.53215 358 -994.51897 -3143.89891 359 -39707.04202 -994.51897 360 -39573.23068 -39707.04202 361 -37838.50594 -39573.23068 362 -39636.06178 -37838.50594 363 -40458.38440 -39636.06178 364 -40969.50230 -40458.38440 365 -39449.80800 -40969.50230 366 -40340.98116 -39449.80800 367 -17035.73919 -40340.98116 368 -21466.55339 -17035.73919 369 7494.40813 -21466.55339 370 -35040.90494 7494.40813 371 -40625.71509 -35040.90494 372 -40309.54313 -40625.71509 373 -17047.24691 -40309.54313 374 -21439.99830 -17047.24691 375 6479.23130 -21439.99830 376 40559.77486 6479.23130 377 46349.40689 40559.77486 378 27079.60963 46349.40689 379 -8345.40199 27079.60963 380 64374.71861 -8345.40199 381 -728.90531 64374.71861 382 -1507.52730 -728.90531 383 -1461.93879 -1507.52730 384 25469.81413 -1461.93879 385 37416.92062 25469.81413 386 67631.12268 37416.92062 387 314.21187 67631.12268 388 -1522.53307 314.21187 389 -42907.33711 -1522.53307 390 -27343.09452 -42907.33711 391 -12813.88312 -27343.09452 392 -16324.54732 -12813.88312 393 10583.70460 -16324.54732 394 -23348.28738 10583.70460 395 101553.74211 -23348.28738 396 -16925.06079 101553.74211 397 -18429.04672 -16925.06079 398 -19083.28176 -18429.04672 399 163320.86876 -19083.28176 400 -19815.25031 163320.86876 401 -23788.13724 -19815.25031 402 -58835.60030 -23788.13724 403 -21234.25979 -58835.60030 404 -50323.72166 -21234.25979 405 -50051.89743 -50323.72166 406 3562.84132 -50051.89743 407 -1388.63677 3562.84132 408 33089.10437 -1388.63677 409 -3058.39376 33089.10437 410 15388.11982 -3058.39376 411 72826.64188 15388.11982 412 179.52151 72826.64188 413 -2839.08546 179.52151 414 -1553.85198 -2839.08546 415 4026.25476 -1553.85198 416 -3363.22032 4026.25476 417 38736.49310 -3363.22032 418 37995.21619 38736.49310 419 -2675.89132 37995.21619 420 -42342.84306 -2675.89132 421 -40138.29311 -42342.84306 422 -40476.12267 -40138.29311 423 -40845.27404 -40476.12267 424 -40542.84185 -40845.27404 425 -40358.05058 -40542.84185 426 -40235.90951 -40358.05058 427 -40418.50850 -40235.90951 428 -19498.20932 -40418.50850 429 -16882.14474 -19498.20932 430 -2053.10198 -16882.14474 > 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/72ga01291404004.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/82ga01291404004.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/9dprl1291404004.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/106y8o1291404004.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/11rhou1291404004.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/12chni1291404004.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/13102b1291404004.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/14cs1e1291404004.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/158jz51291404004.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/16t2xb1291404004.tab") + } > > try(system("convert tmp/16our1291404004.ps tmp/16our1291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/26our1291404004.ps tmp/26our1291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/3zxbu1291404004.ps tmp/3zxbu1291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/4zxbu1291404004.ps tmp/4zxbu1291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/5zxbu1291404004.ps tmp/5zxbu1291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/6aosf1291404004.ps tmp/6aosf1291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/72ga01291404004.ps tmp/72ga01291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/82ga01291404004.ps tmp/82ga01291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/9dprl1291404004.ps tmp/9dprl1291404004.png",intern=TRUE)) character(0) > try(system("convert tmp/106y8o1291404004.ps tmp/106y8o1291404004.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 15.190 2.250 31.618