R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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 + ,10 + ,60720 + ,319369 + ,143 + ,7022 + ,75 + ,0 + ,61 + ,0 + ,5565 + ,47 + ,43 + ,94355 + ,493408 + ,38 + ,6243 + ,53 + ,39 + ,287 + ,0 + ,601 + ,6 + ,5 + ,60720 + ,319210 + ,144 + ,19868 + ,198 + ,0 + ,268 + ,1 + ,188 + ,11 + ,9 + ,77655 + ,381180 + ,67 + ,16471 + ,964 + ,0 + ,25 + ,1 + ,7146 + ,105 + ,83 + ,134028 + ,297978 + ,367 + ,933 + ,14 + ,305 + ,418 + ,0 + ,1135 + ,17 + ,20 + ,62285 + ,290476 + ,384 + ,5322 + ,80 + ,-9 + ,392 + ,1 + ,450 + ,8 + ,6 + ,59325 + ,292136 + ,380 + ,11517 + ,205 + ,-1 + ,131 + ,1 + ,34 + ,8 + ,7 + ,60630 + ,314353 + ,309 + ,14294 + ,3340 + ,0 + ,316 + ,1 + ,133 + ,14 + ,4 + ,65990 + ,339445 + ,109 + ,9960 + ,1052 + ,0 + ,347 + ,1 + ,119 + ,6 + ,3 + ,59118 + ,303677 + ,354 + ,17280 + ,872 + ,0 + ,68 + ,0 + ,2053 + ,53 + ,44 + ,100423 + ,397144 + ,60 + ,3720 + ,96 + ,33 + ,70 + ,0 + ,4036 + ,63 + ,85 + ,100269 + ,424898 + ,53 + ,3570 + ,56 + ,32 + ,147 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,0 + ,169 + ,1 + ,4655 + ,393 + ,55 + ,60720 + ,341570 + ,105 + ,360 + ,30 + ,0 + ,248 + ,1 + ,131 + ,11 + ,15 + ,61808 + ,308989 + ,338 + ,9908 + ,834 + ,0 + ,152 + ,1 + ,1766 + ,73 + ,46 + ,60438 + ,305959 + ,349 + ,1451 + ,60 + ,0 + ,351 + ,1 + ,312 + ,14 + ,3 + ,58598 + ,318690 + ,147 + ,8478 + ,381 + ,0 + ,372 + ,1 + ,448 + ,40 + ,7 + ,59781 + ,323361 + ,132 + ,3084 + ,275 + ,0 + ,137 + ,1 + ,115 + ,13 + ,10 + ,60945 + ,318903 + ,145 + ,9146 + ,1037 + ,0 + ,352 + ,1 + ,60 + ,10 + ,7 + ,61124 + ,314049 + ,314 + ,11405 + ,1916 + ,0 + ,338 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,180 + ,0 + ,258 + ,1 + ,364 + ,35 + ,26 + ,47705 + ,298466 + ,365 + ,2813 + ,271 + ,0 + ,343 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,183 + ,0 + ,398 + ,0 + ,1442 + ,118 + ,239 + ,121173 + ,438493 + ,49 + ,2021 + ,165 + ,-3 + ,115 + ,0 + ,1389 + ,9 + ,6 + ,57530 + ,378049 + ,71 + ,19783 + ,128 + ,6 + ,188 + ,1 + ,149 + ,5 + ,6 + ,62041 + ,313332 + ,318 + ,22666 + ,761 + ,0 + ,379 + ,1 + ,2212 + ,14 + ,10 + ,60720 + ,319864 + ,142 + ,8562 + ,54 + ,0 + ,48 + ,0 + ,7489 + ,322 + ,149 + ,136996 + ,430866 + ,50 + ,717 + ,31 + ,70 + ,366 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,287 + ,0 + ,117 + ,1 + ,402 + ,26 + ,8 + ,84990 + ,332743 + ,115 + ,5106 + ,330 + ,5 + ,412 + ,0 + ,7419 + ,29 + ,21 + ,63255 + ,286963 + ,387 + ,2999 + ,12 + ,-5 + ,157 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,219 + ,0 + ,264 + ,1 + ,307 + ,15 + ,10 + ,58856 + ,331955 + ,116 + ,8797 + ,430 + ,0 + ,54 + ,0 + ,1134 + ,162 + ,63 + ,146216 + ,527021 + ,34 + ,2019 + ,288 + ,58 + ,64 + ,0 + ,7561 + ,87 + ,77 + ,82425 + ,364304 + ,86 + ,1889 + ,22 + ,37 + ,113 + ,1 + ,5131 + ,71 + ,44 + ,86111 + ,292154 + ,379 + ,1298 + ,18 + ,9 + ,148 + ,1 + ,9 + ,6 + ,7 + ,60835 + ,314987 + ,300 + ,19165 + ,12458 + ,0 + ,425 + ,1 + ,2264 + ,24 + ,20 + ,55174 + ,272713 + ,394 + ,3030 + ,32 + ,-31 + ,7 + ,1 + ,40949 + ,481 + ,263 + ,129352 + ,1398893 + ,9 + ,2493 + ,29 + ,2865 + ,197 + ,0 + ,4980 + ,91 + ,89 + ,113521 + ,409642 + ,55 + ,2304 + ,42 + ,0 + ,97 + ,1 + ,272 + ,8 + ,9 + ,112995 + ,429112 + ,51 + ,28639 + ,844 + ,14 + ,168 + ,1 + ,757 + ,19 + ,18 + ,45689 + ,382712 + ,66 + ,9616 + ,242 + ,0 + ,69 + ,1 + ,1424 + ,94 + ,73 + ,131116 + ,374943 + ,77 + ,1861 + ,123 + ,32 + ,127 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,191 + ,0 + ,184 + ,0 + ,378 + ,52 + ,16 + ,71873 + ,297413 + ,369 + ,1873 + ,258 + ,0 + ,401 + ,1 + ,1061 + ,5 + ,3 + ,70415 + ,304252 + ,353 + ,20850 + ,98 + ,-4 + ,371 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,290 + ,0 + ,174 + ,1 + ,203 + ,9 + ,9 + ,61938 + ,309596 + ,334 + ,12177 + ,541 + ,0 + ,46 + ,0 + ,3689 + ,25 + ,25 + ,101792 + ,377305 + ,74 + ,7092 + ,48 + ,74 + ,90 + ,1 + ,567 + ,11 + ,14 + ,103345 + ,386212 + ,65 + ,16928 + ,329 + ,20 + ,27 + ,0 + ,27664 + ,282 + ,236 + ,74996 + ,657954 + ,24 + ,1624 + ,17 + ,247 + ,105 + ,1 + ,1686 + ,40 + ,28 + ,68696 + ,368186 + ,81 + ,4205 + ,100 + ,11 + ,413 + ,0 + ,1164 + ,17 + ,15 + ,57320 + ,269753 + ,395 + ,4103 + ,60 + ,-6 + ,26 + ,1 + ,15824 + ,64 + ,56 + ,63838 + ,585715 + ,30 + ,6027 + ,24 + ,284 + ,16 + ,1 + ,6575 + ,159 + ,63 + ,151352 + ,1071292 + ,11 + ,5480 + ,133 + ,819 + ,19 + ,0 + ,37597 + ,248 + ,234 + ,37884 + ,992426 + ,12 + ,3195 + ,21 + ,652 + ,336 + ,1 + ,314 + ,5 + ,4 + ,60720 + ,315371 + ,294 + ,23074 + ,367 + ,0 + ,110 + ,0 + ,1932 + ,36 + ,20 + ,72954 + ,249898 + ,403 + ,1386 + ,26 + ,10 + ,86 + ,0 + ,5692 + ,72 + ,68 + ,64239 + ,357312 + ,91 + ,2185 + ,28 + ,22 + ,294 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,257 + ,0 + ,322 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,177 + ,0 + ,94 + ,0 + ,6533 + ,75 + ,29 + ,126910 + ,364839 + ,85 + ,2198 + ,25 + ,16 + ,400 + ,0 + ,2126 + ,138 + ,41 + ,21001 + ,230621 + ,413 + ,222 + ,14 + ,-3 + ,150 + ,1 + ,14 + ,1 + ,0 + ,61016 + ,315877 + ,161 + ,115877 + ,8041 + ,0 + ,107 + ,0 + ,3192 + ,60 + ,48 + ,59831 + ,125390 + ,427 + ,-1243 + ,-23 + ,10 + ,166 + ,1 + ,47 + ,3 + ,1 + ,60720 + ,314882 + ,305 + ,38294 + ,2432 + ,0 + ,47 + ,0 + ,4891 + ,217 + ,49 + ,210568 + ,370837 + ,79 + ,787 + ,35 + ,73 + ,182 + ,0 + ,531 + ,10 + ,11 + ,71641 + ,330068 + ,120 + ,13007 + ,245 + ,0 + ,245 + ,1 + ,2332 + ,4 + ,2 + ,72680 + ,324385 + ,131 + ,31096 + ,53 + ,0 + ,330 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,277 + ,0 + ,10 + ,0 + ,16275 + ,109 + ,129 + ,47818 + ,1443586 + ,7 + ,11409 + ,76 + ,1974 + ,429 + ,1 + ,9252 + ,48 + ,32 + ,79510 + ,253537 + ,401 + ,1115 + ,6 + ,-129 + ,2 + ,0 + ,29790 + ,309 + ,444 + ,81767 + ,4321023 + ,2 + ,13337 + ,138 + ,28328 + ,235 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,204 + ,0 + ,187 + ,0 + ,2300 + ,29 + ,23 + ,64873 + ,352108 + ,97 + ,5245 + ,66 + ,0 + ,213 + ,1 + ,267 + ,32 + ,11 + ,57640 + ,330059 + ,121 + ,4064 + ,487 + ,0 + ,296 + ,0 + ,382 + ,21 + ,9 + ,59365 + ,340968 + ,107 + ,6713 + ,369 + ,0 + ,218 + ,1 + ,397 + ,28 + ,7 + ,60646 + ,317736 + ,152 + ,4205 + ,297 + ,0 + ,28 + ,1 + ,11920 + ,182 + ,73 + ,31080 + ,209458 + ,419 + ,52 + ,1 + ,238 + ,5 + ,0 + ,54660 + ,302 + ,315 + ,126942 + ,1491348 + ,6 + ,4276 + ,24 + ,3040 + ,138 + ,1 + ,5 + ,2 + ,1 + ,60720 + ,314887 + ,304 + ,57443 + ,24814 + ,0 + ,163 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,0 + ,219 + ,1 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,199 + ,0 + ,125 + ,0 + ,2647 + ,37 + ,17 + ,62920 + ,353058 + ,95 + ,4137 + ,58 + ,1 + ,216 + ,1 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,201 + ,0 + ,196 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,263 + ,0 + ,128 + ,1 + ,94 + ,7 + ,10 + ,60793 + ,314533 + ,308 + ,16362 + ,1212 + ,0 + ,223 + ,1 + ,3 + ,2 + ,1 + ,60698 + ,315354 + ,296 + ,57677 + ,36505 + ,0 + ,165 + ,1 + ,454 + ,51 + ,24 + ,63261 + ,302187 + ,357 + ,2004 + ,225 + ,0 + ,373 + ,1 + ,227 + ,9 + ,7 + ,76644 + ,336639 + ,110 + ,15182 + ,603 + ,0 + ,324 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,245 + ,0 + ,251 + ,1 + ,206 + ,9 + ,9 + ,53110 + ,296702 + ,372 + ,10745 + ,469 + ,0 + ,11 + ,1 + ,25830 + ,115 + ,104 + ,245546 + ,1073089 + ,10 + ,7592 + ,34 + ,1910 + ,62 + ,1 + ,1528 + ,38 + ,32 + ,60326 + ,146494 + ,426 + ,-1408 + ,-35 + ,38 + ,149 + ,1 + ,271 + ,2 + ,2 + ,60720 + ,325249 + ,129 + ,62625 + ,463 + ,0 + ,270 + ,1 + ,138 + ,6 + ,26 + ,69817 + ,331420 + ,117 + ,21903 + ,954 + ,0 + ,171 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,211 + ,0 + ,226 + ,1 + ,278 + ,26 + ,22 + ,61404 + ,314922 + ,302 + ,4420 + ,413 + ,0 + ,276 + ,1 + ,282 + ,30 + ,12 + ,62452 + ,320016 + ,141 + ,4001 + ,425 + ,0 + ,383 + ,1 + ,571 + ,80 + ,14 + ,41477 + ,280398 + ,390 + ,1005 + ,141 + ,0 + ,93 + ,1 + ,2253 + ,99 + ,68 + ,63593 + ,452469 + ,43 + ,2550 + ,112 + ,17 + ,302 + ,1 + ,290 + ,2 + ,1 + ,58790 + ,301164 + ,360 + ,50582 + ,348 + ,0 + ,385 + ,1 + ,78 + ,16 + ,12 + ,62700 + ,317330 + ,154 + ,7333 + ,1511 + ,0 + ,225 + ,1 + ,20 + ,2 + ,1 + ,60805 + ,315576 + ,169 + ,57788 + ,5924 + ,0 + ,275 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,284 + ,0 + ,431 + ,1 + ,18073 + ,180 + ,69 + ,27284 + ,-7170 + ,429 + ,-1151 + ,-11 + ,-366 + ,220 + ,1 + ,1866 + ,29 + ,28 + ,56225 + ,322331 + ,135 + ,4218 + ,66 + ,0 + ,6 + ,0 + ,42634 + ,156 + ,168 + ,157214 + ,1629616 + ,5 + ,9164 + ,34 + ,2964 + ,290 + ,1 + ,249 + ,62 + ,13 + ,54323 + ,292754 + ,378 + ,1496 + ,373 + ,0 + ,129 + ,1 + ,422 + ,28 + ,23 + ,57935 + ,318056 + ,150 + ,4216 + ,280 + ,0 + ,65 + ,1 + ,2675 + ,79 + ,54 + ,59017 + ,355178 + ,93 + ,1964 + ,58 + ,36 + ,416 + ,1 + ,965 + ,3 + ,2 + ,73490 + ,204325 + ,421 + ,1442 + ,4 + ,-8 + ,299 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,224 + ,0 + ,407 + ,1 + ,621 + ,16 + ,11 + ,68005 + ,317046 + ,156 + ,7315 + ,188 + ,-5 + ,252 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,240 + ,0 + ,257 + ,1 + ,365 + ,8 + ,7 + ,54820 + ,309560 + ,335 + ,13695 + ,300 + ,0 + ,63 + ,0 + ,3122 + ,34 + ,25 + ,94670 + ,414462 + ,54 + ,6308 + ,69 + ,38 + ,21 + ,0 + ,12988 + ,119 + ,73 + ,82340 + ,857217 + ,16 + ,5523 + ,51 + ,601 + ,38 + ,0 + ,5336 + ,81 + ,80 + ,112477 + ,697458 + ,22 + ,6141 + ,93 + ,132 + ,51 + ,0 + ,3160 + ,108 + ,90 + ,108094 + ,530670 + ,33 + ,3062 + ,105 + ,66 + ,386 + ,1 + ,2489 + ,40 + ,29 + ,67804 + ,238125 + ,408 + ,953 + ,15 + ,0 + ,325 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,273 + ,0 + ,33 + ,0 + ,6984 + ,45 + ,40 + ,80570 + ,741409 + ,19 + ,12031 + ,78 + ,195 + ,91 + ,1 + ,2001 + ,52 + ,36 + ,95551 + ,393343 + ,62 + ,3718 + ,97 + ,19 + ,430 + ,1 + ,15236 + ,90 + ,86 + ,77440 + ,372631 + ,78 + ,1918 + ,11 + ,-158 + ,259 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,239 + ,0 + ,382 + ,1 + ,530 + ,60 + ,18 + ,73433 + ,317291 + ,155 + ,1955 + ,221 + ,0 + ,369 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,292 + ,0 + ,428 + ,1 + ,35624 + ,223 + ,244 + ,157278 + ,306275 + ,347 + ,477 + ,3 + ,-70 + ,100 + ,0 + ,1383 + ,11 + ,9 + ,73221 + ,317892 + ,151 + ,10717 + ,85 + ,13 + ,122 + ,1 + ,875 + ,14 + ,12 + ,67000 + ,334280 + ,112 + ,9591 + ,153 + ,3 + ,198 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,259 + ,0 + ,289 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,254 + ,0 + ,307 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,222 + ,0 + ,389 + ,1 + ,72 + ,3 + ,3 + ,60398 + ,314210 + ,312 + ,38070 + ,1591 + ,0 + ,210 + ,0 + ,265 + ,31 + ,3 + ,60720 + ,306948 + ,346 + ,3450 + ,404 + ,194 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,186 + ,0 + ,277 + ,1 + ,335 + ,33 + ,23 + ,64175 + ,320398 + ,140 + ,3648 + ,359 + ,0 + ,298 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,230 + ,0 + ,356 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,269 + ,0 + ,43 + ,1 + ,4525 + ,226 + ,89 + ,93811 + ,501749 + ,37 + ,1335 + ,67 + ,101 + ,71 + ,0 + ,3045 + ,58 + ,49 + ,27330 + ,202055 + ,422 + ,35 + ,1 + ,31 + ,133 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,196 + ,0 + ,212 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,0 + ,186 + ,0 + ,638 + ,14 + ,11 + ,60370 + ,333210 + ,114 + ,9515 + ,209 + ,0 + ,262 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,283 + ,0 + ,378 + ,1 + ,607 + ,11 + ,7 + ,60436 + ,322340 + ,134 + ,11122 + ,201 + ,0 + ,132 + ,1 + ,1558 + ,63 + ,25 + ,55637 + ,369448 + ,80 + ,2690 + ,109 + ,0 + ,279 + ,0 + ,1324 + ,117 + ,29 + ,67440 + ,291841 + ,381 + ,785 + ,69 + ,0 + ,376 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,280 + ,0 + ,267 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,0 + ,403 + ,1 + ,611 + ,8 + ,6 + ,59190 + ,296919 + ,371 + ,12115 + ,159 + ,-4 + ,303 + ,1 + ,0 + ,0 + ,1 + ,60720 + ,315380 + ,225 + ,0 + ,388 + ,1 + ,923 + ,8 + ,6 + ,58620 + ,309038 + ,337 + ,13630 + ,118 + ,0 + ,420 + ,1 + ,661 + ,3 + ,1 + ,65920 + ,246541 + ,405 + ,15514 + ,70 + ,-12 + ,114 + ,0 + ,2397 + ,14 + ,13 + ,74020 + ,289513 + ,385 + ,6394 + ,37 + ,6 + ,305 + ,1 + ,366 + ,53 + ,14 + ,61808 + ,344425 + ,101 + ,2725 + ,394 + ,0 + ,346 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,181 + ,0 + ,354 + ,1 + ,135 + ,2 + ,1 + ,62065 + ,314210 + ,311 + ,57105 + ,847 + ,0 + ,124 + ,0 + ,1659 + ,29 + ,25 + ,107577 + ,480382 + ,40 + ,9668 + ,169 + ,2 + ,136 + ,1 + ,316 + ,9 + ,4 + ,60505 + ,315009 + ,299 + ,12779 + ,364 + ,0 + ,130 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,197 + ,0 + ,146 + ,1 + ,309 + ,11 + ,9 + ,56535 + ,312878 + ,322 + ,10262 + ,365 + ,0 + ,314 + ,0 + ,49 + ,8 + ,14 + ,64107 + ,322031 + ,137 + ,15254 + ,2500 + ,0 + ,241 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,200 + ,0 + ,34 + ,0 + ,4519 + ,49 + ,42 + ,102129 + ,597793 + ,29 + ,8118 + ,88 + ,163 + ,384 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,0 + ,120 + ,1 + ,837 + ,69 + ,21 + ,61262 + ,315688 + ,165 + ,1677 + ,138 + ,4 + ,72 + ,1 + ,5119 + ,49 + ,38 + ,39039 + ,378525 + ,69 + ,3643 + ,35 + ,31 + ,79 + ,1 + ,1280 + ,117 + ,52 + ,69465 + ,312378 + ,326 + ,960 + ,88 + ,26 + ,99 + ,1 + ,2564 + ,22 + ,19 + ,130140 + ,403560 + ,59 + ,9253 + ,79 + ,13 + ,108 + ,0 + ,2045 + ,20 + ,16 + ,97890 + ,510834 + ,36 + ,15542 + ,152 + ,10 + ,119 + ,0 + ,2234 + ,16 + ,11 + ,77200 + ,214215 + ,418 + ,888 + ,6 + ,5 + ,102 + ,1 + ,975 + ,23 + ,19 + ,90534 + ,235133 + ,410 + ,1528 + ,36 + ,13 + ,209 + ,0 + ,1136 + ,32 + ,26 + ,48522 + ,343613 + ,103 + ,4488 + ,126 + ,0 + ,293 + ,1 + ,453 + ,21 + ,19 + ,81125 + ,365959 + ,84 + ,7903 + ,366 + ,0 + ,211 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,210 + ,0 + ,227 + ,1 + ,61 + ,14 + ,9 + ,60720 + ,314551 + ,307 + ,8182 + ,1874 + ,0 + ,234 + ,1 + ,368 + ,17 + ,13 + ,61686 + ,303230 + ,355 + ,6072 + ,281 + ,0 + ,247 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,233 + ,0 + ,49 + ,1 + ,4901 + ,84 + ,58 + ,121920 + ,469107 + ,41 + ,3204 + ,55 + ,67 + ,88 + ,1 + ,540 + ,14 + ,13 + ,103960 + ,354228 + ,94 + ,11016 + ,286 + ,22 + ,329 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,253 + ,0 + ,160 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,213 + ,0 + ,271 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,286 + ,0 + ,306 + ,1 + ,2 + ,9 + ,12 + ,60735 + ,315394 + ,174 + ,12822 + ,63057 + ,0 + ,172 + ,1 + ,36 + ,4 + ,3 + ,61564 + ,312412 + ,325 + ,28103 + ,3146 + ,0 + ,402 + ,1 + ,776 + ,9 + ,8 + ,64230 + ,333505 + ,113 + ,14834 + ,172 + ,-4 + ,4 + ,1 + ,84738 + ,588 + ,428 + ,-26007 + ,223193 + ,415 + ,39 + ,0 + ,4145 + ,156 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,217 + ,0 + ,315 + ,1 + ,3 + ,4 + ,8 + ,60761 + ,315656 + ,167 + ,28914 + ,37920 + ,0 + ,399 + ,1 + ,529 + ,9 + ,7 + ,63870 + ,296261 + ,374 + ,10696 + ,182 + ,-3 + ,381 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,264 + ,0 + ,159 + ,1 + ,405 + ,7 + ,5 + ,60845 + ,336425 + ,111 + ,19489 + ,337 + ,0 + ,181 + ,1 + ,972 + ,34 + ,19 + ,71642 + ,359335 + ,87 + ,4686 + ,164 + ,0 + ,278 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,232 + ,0 + ,89 + ,1 + ,2099 + ,64 + ,45 + ,106611 + ,308636 + ,339 + ,1697 + ,52 + ,20 + ,112 + ,0 + ,3437 + ,47 + ,33 + ,48022 + ,158492 + ,425 + ,-883 + ,-12 + ,9 + ,201 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,243 + ,0 + ,344 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,182 + ,0 + ,22 + ,0 + ,22330 + ,84 + ,67 + ,79801 + ,711969 + ,20 + ,6095 + ,23 + ,446 + ,291 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,258 + ,0 + ,269 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,221 + ,0 + ,164 + ,1 + ,483 + ,27 + ,7 + ,60830 + ,306268 + ,348 + ,3936 + ,220 + ,0 + ,358 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,279 + ,0 + ,118 + ,1 + ,2239 + ,21 + ,17 + ,88590 + ,442882 + ,47 + ,11566 + ,108 + ,5 + ,92 + ,0 + ,2949 + ,41 + ,40 + ,82903 + ,378509 + ,70 + ,4354 + ,61 + ,18 + ,300 + ,1 + ,0 + ,0 + ,9 + ,60720 + ,315380 + ,223 + ,0 + ,312 + ,1 + ,365 + ,83 + ,32 + ,87192 + ,346611 + ,99 + ,1766 + ,402 + ,0 + ,183 + ,1 + ,2461 + ,57 + ,25 + ,55792 + ,314289 + ,310 + ,2005 + ,46 + ,0 + ,18 + ,0 + ,21950 + ,519 + ,265 + ,114337 + ,856956 + ,17 + ,1266 + ,30 + ,745 + ,126 + ,0 + ,3294 + ,23 + ,23 + ,75832 + ,217193 + ,416 + ,748 + ,5 + ,1 + ,193 + ,1 + ,141 + ,6 + ,7 + ,61630 + ,315366 + ,295 + ,19228 + ,819 + ,0 + ,390 + ,1 + ,572 + ,16 + ,8 + ,58580 + ,307930 + ,342 + ,6746 + ,189 + ,-1 + ,23 + ,1 + ,13326 + ,102 + ,103 + ,165548 + ,702380 + ,21 + ,4925 + ,38 + ,436 + ,50 + ,1 + ,2284 + ,33 + ,41 + ,76403 + ,194493 + ,423 + ,-167 + ,-2 + ,66 + ,206 + ,0 + ,10 + ,2 + ,5 + ,61656 + ,316155 + ,160 + ,58077 + ,11639 + ,0 + ,162 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,207 + ,0 + ,98 + ,0 + ,1414 + ,198 + ,22 + ,70184 + ,330546 + ,119 + ,659 + ,92 + ,13 + ,75 + ,1 + ,1975 + ,35 + ,34 + ,118881 + ,394510 + ,61 + ,5557 + ,98 + ,30 + ,134 + ,1 + ,43 + ,3 + ,20 + ,60887 + ,312846 + ,323 + ,37615 + ,2613 + ,0 + ,180 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,190 + ,0 + ,337 + ,1 + ,844 + ,78 + ,20 + ,60925 + ,296139 + ,375 + ,1233 + ,114 + ,0 + ,341 + ,1 + ,304 + ,15 + ,12 + ,62969 + ,295580 + ,376 + ,6372 + ,314 + ,0 + ,154 + ,1 + ,458 + ,11 + ,7 + ,58625 + ,297765 + ,368 + ,8888 + ,213 + ,0 + ,36 + ,1 + ,18562 + ,155 + ,211 + ,102313 + ,377934 + ,72 + ,1148 + ,10 + ,143 + ,192 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,179 + ,0 + ,17 + ,1 + ,7123 + ,109 + ,102 + ,288170 + ,638830 + ,28 + ,4026 + ,62 + ,769 + ,254 + ,1 + ,622 + ,73 + ,38 + ,73007 + ,304376 + ,352 + ,1430 + ,168 + ,0 + ,374 + ,1 + ,174 + ,1 + ,3 + ,64820 + ,307424 + ,344 + ,107424 + ,618 + ,0 + ,35 + ,1 + ,2220 + ,22 + ,12 + ,301670 + ,644190 + ,26 + ,20190 + ,200 + ,154 + ,377 + ,1 + ,121 + ,24 + ,11 + ,56178 + ,295370 + ,377 + ,3974 + ,788 + ,0 + ,32 + ,0 + ,11819 + ,85 + ,70 + ,106113 + ,574339 + ,31 + ,4404 + ,32 + ,199 + ,260 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,235 + ,0 + ,274 + ,1 + ,125 + ,12 + ,19 + ,60798 + ,310201 + ,332 + ,9183 + ,883 + ,0 + ,406 + ,1 + ,1182 + ,13 + ,11 + ,70694 + ,327007 + ,124 + ,9770 + ,107 + ,-5 + ,116 + ,1 + ,1503 + ,39 + ,35 + ,56364 + ,343466 + ,104 + ,3679 + ,95 + ,6 + ,320 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,274 + ,0 + ,367 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,248 + ,0 + ,288 + ,1 + ,30 + ,4 + ,3 + ,62045 + ,318098 + ,149 + ,29525 + ,3912 + ,0 + ,414 + ,0 + ,3310 + ,33 + ,32 + ,75230 + ,448243 + ,44 + ,7523 + ,75 + ,-6 + ,74 + ,1 + ,554 + ,11 + ,9 + ,79285 + ,325738 + ,125 + ,11431 + ,227 + ,30 + ,250 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,234 + ,0 + ,224 + ,1 + ,468 + ,21 + ,20 + ,60720 + ,312161 + ,328 + ,5341 + ,239 + ,0 + ,397 + ,1 + ,4917 + ,42 + ,26 + ,52811 + ,243650 + ,406 + ,1039 + ,9 + ,-3 + ,417 + ,1 + ,3256 + ,66 + ,23 + ,35250 + ,407159 + ,57 + ,3139 + ,64 + ,-8 + ,323 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,185 + ,0 + ,340 + ,1 + ,125 + ,13 + ,10 + ,59734 + ,317698 + ,153 + ,9054 + ,944 + ,0 + ,239 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,195 + ,0 + ,326 + ,1 + ,22 + ,2 + ,1 + ,60722 + ,312502 + ,324 + ,56251 + ,5154 + ,0 + ,261 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,262 + ,0 + ,357 + ,1 + ,514 + ,52 + ,15 + ,78780 + ,322378 + ,133 + ,2353 + ,238 + ,0 + ,313 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,0 + ,360 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,270 + ,0 + ,272 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,285 + ,0 + ,333 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,276 + ,0 + ,355 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,272 + ,0 + ,362 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,268 + ,0 + ,37 + ,1 + ,10327 + ,91 + ,74 + ,88577 + ,640273 + ,27 + ,4838 + ,43 + ,135 + ,42 + ,1 + ,13718 + ,61 + ,56 + ,96448 + ,345783 + ,100 + ,2390 + ,11 + ,104 + ,57 + ,1 + ,3748 + ,9 + ,6 + ,50350 + ,652925 + ,25 + ,50325 + ,121 + ,53 + ,30 + ,1 + ,14416 + ,361 + ,139 + ,49857 + ,439798 + ,48 + ,664 + ,17 + ,220 + ,404 + ,0 + ,1526 + ,25 + ,21 + ,69351 + ,278990 + ,391 + ,3160 + ,52 + ,-4 + ,80 + ,0 + ,666 + ,171 + ,13 + ,117869 + ,339836 + ,108 + ,818 + ,210 + ,25 + ,422 + ,1 + ,2844 + ,67 + ,74 + ,72683 + ,240897 + ,407 + ,610 + ,14 + ,-15 + ,297 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,229 + ,0 + ,318 + ,1 + ,368 + ,17 + ,20 + ,61167 + ,297141 + ,370 + ,5714 + ,264 + ,0 + ,175 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,205 + ,0 + ,285 + ,1 + ,333 + ,13 + ,9 + ,70811 + ,331323 + ,118 + ,10102 + ,394 + ,0 + ,339 + ,1 + ,26 + ,1 + ,1 + ,60896 + ,313880 + ,315 + ,113880 + ,4328 + ,0 + ,243 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,237 + ,0 + ,167 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,193 + ,0 + ,411 + ,1 + ,1303 + ,61 + ,57 + ,69863 + ,309422 + ,336 + ,1794 + ,84 + ,-5 + ,155 + ,1 + ,20 + ,3 + ,2 + ,60938 + ,315245 + ,297 + ,38415 + ,5674 + ,0 + ,56 + ,1 + ,2384 + ,97 + ,77 + ,61348 + ,405972 + ,58 + ,2123 + ,86 + ,54 + ,204 + ,0 + ,203 + ,27 + ,12 + ,50804 + ,300962 + ,361 + ,3739 + ,498 + ,0 + ,178 + ,0 + ,71 + ,1 + ,0 + ,60745 + ,316176 + ,159 + ,116176 + ,1646 + ,0 + ,295 + ,1 + ,53 + ,14 + ,10 + ,59506 + ,302409 + ,356 + ,7315 + ,1932 + ,0 + ,396 + ,1 + ,562 + ,13 + ,8 + ,58456 + ,283587 + ,389 + ,6430 + ,149 + ,-3 + ,393 + ,1 + ,622 + ,7 + ,8 + ,60950 + ,263276 + ,398 + ,9039 + ,102 + ,-2 + ,135 + ,1 + ,645 + ,5 + ,4 + ,60720 + ,312075 + ,329 + ,22415 + ,174 + ,0 + ,364 + ,0 + ,1763 + ,24 + ,28 + ,61600 + ,308336 + ,340 + ,4514 + ,61 + ,0 + ,179 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,0 + ,231 + ,1 + ,317 + ,9 + ,7 + ,63915 + ,298700 + ,363 + ,10967 + ,311 + ,0 + ,319 + ,1 + ,1 + ,6 + ,5 + ,60719 + ,315372 + ,293 + ,19229 + ,84213 + ,0 + ,191 + ,1 + ,275 + ,4 + ,7 + ,59500 + ,318745 + ,146 + ,29686 + ,431 + ,0 + ,321 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,246 + ,0 + ,111 + ,1 + ,936 + ,30 + ,22 + ,67939 + ,408881 + ,56 + ,6963 + ,223 + ,9 + ,29 + ,1 + ,8568 + ,37 + ,34 + ,32168 + ,786690 + ,18 + ,15856 + ,68 + ,223 + ,423 + ,0 + ,11528 + ,70 + ,68 + ,-14545 + ,-83265 + ,431 + ,-4047 + ,-25 + ,-25 + ,202 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,250 + ,0 + ,232 + ,1 + ,738 + ,10 + ,9 + ,60720 + ,321376 + ,139 + ,12138 + ,164 + ,0 + ,266 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,281 + ,0 + ,405 + ,1 + ,592 + ,7 + ,7 + ,64270 + ,276898 + ,392 + ,10985 + ,130 + ,-4 + ,176 + ,1 + ,126 + ,10 + ,8 + ,60951 + ,315547 + ,171 + ,11555 + ,921 + ,0 + ,236 + ,1 + ,2 + ,1 + ,0 + ,60743 + ,315487 + ,172 + ,115487 + ,52257 + ,0 + ,142 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,220 + ,0 + ,13 + ,0 + ,18014 + ,239 + ,393 + ,-1710 + ,1405225 + ,8 + ,5043 + ,67 + ,943 + ,281 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,256 + ,0 + ,328 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,0 + ,9 + ,0 + ,37704 + ,452 + ,228 + ,60448 + ,983660 + ,14 + ,1734 + ,21 + ,2167 + ,207 + ,0 + ,63 + ,4 + ,2 + ,65688 + ,318574 + ,148 + ,29643 + ,1891 + ,0 + ,73 + ,1 + ,1431 + ,34 + ,35 + ,106885 + ,310768 + ,330 + ,3258 + ,77 + ,30 + ,242 + ,1 + ,94 + ,3 + ,3 + ,61360 + ,312887 + ,321 + ,37629 + ,1200 + ,0 + ,229 + ,1 + ,192 + ,4 + ,2 + ,65276 + ,312339 + ,327 + ,28085 + ,586 + ,0 + ,317 + ,1 + ,32 + ,1 + ,0 + ,59988 + ,314964 + ,301 + ,114964 + ,3603 + ,0 + ,44 + ,1 + ,6869 + ,105 + ,88 + ,117520 + ,379983 + ,68 + ,1714 + ,26 + ,90 + ,349 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,176 + ,0 + ,140 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,0 + ,361 + ,1 + ,1 + ,4 + ,4 + ,60722 + ,315398 + ,173 + ,28849 + ,172235 + ,0 + ,203 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,247 + ,0 + ,255 + ,1 + ,2328 + ,94 + ,57 + ,82732 + ,253588 + ,400 + ,570 + ,23 + ,0 + ,308 + ,0 + ,209 + ,13 + ,19 + ,64016 + ,316647 + ,157 + ,8973 + ,558 + ,0 + ,195 + ,1 + ,28 + ,2 + ,3 + ,60890 + ,315688 + ,166 + ,57844 + ,4148 + ,0 + ,280 + ,1 + ,176 + ,28 + ,8 + ,68136 + ,310670 + ,331 + ,3952 + ,629 + ,0 + ,415 + ,1 + ,1920 + ,15 + ,11 + ,79420 + ,165404 + ,424 + ,-2306 + ,-18 + ,-6 + ,3 + ,0 + ,87550 + ,458 + ,412 + ,153198 + ,4111912 + ,3 + ,8541 + ,45 + ,17936 + ,151 + ,1 + ,520 + ,10 + ,7 + ,58650 + ,291650 + ,383 + ,9165 + ,176 + ,0 + ,304 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,231 + ,0 + ,421 + ,1 + ,1013 + ,2 + ,1 + ,69770 + ,253468 + ,402 + ,26734 + ,53 + ,-14 + ,217 + ,1 + ,15 + ,8 + ,6 + ,60831 + ,315688 + ,164 + ,14461 + ,7687 + ,0 + ,353 + ,1 + ,587 + ,12 + ,10 + ,59595 + ,325699 + ,126 + ,10475 + ,214 + ,0 + ,58 + ,1 + ,5371 + ,66 + ,47 + ,87720 + ,446211 + ,45 + ,3730 + ,46 + ,49 + ,359 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,251 + ,0 + ,77 + ,1 + ,1012 + ,47 + ,29 + ,114768 + ,368078 + ,82 + ,3576 + ,166 + ,26 + ,228 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,215 + ,0 + ,83 + ,1 + ,876 + ,59 + ,33 + ,138971 + ,352850 + ,96 + ,2591 + ,174 + ,24 + ,1 + ,0 + ,162556 + ,1081 + ,807 + ,213118 + ,6282154 + ,1 + ,5626 + ,37 + ,56289 + ,14 + ,1 + ,43556 + ,247 + ,190 + ,32648 + ,227132 + ,414 + ,110 + ,1 + ,929 + ,395 + ,0 + ,3425 + ,12 + ,9 + ,83620 + ,283910 + ,388 + ,6993 + ,25 + ,-3 + ,78 + ,1 + ,810 + ,43 + ,11 + ,74015 + ,236761 + ,409 + ,855 + ,45 + ,26 + ,190 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,184 + ,0 + ,311 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,244 + ,0 + ,39 + ,0 + ,2365 + ,79 + ,83 + ,191778 + ,550608 + ,32 + ,4438 + ,148 + ,120 + ,101 + ,0 + ,1261 + ,33 + ,31 + ,76114 + ,307528 + ,343 + ,3258 + ,85 + ,13 + ,256 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,238 + ,0 + ,123 + ,1 + ,1585 + ,107 + ,112 + ,93099 + ,355864 + ,92 + ,1457 + ,98 + ,2 + ,24 + ,1 + ,16189 + ,295 + ,290 + ,116384 + ,358589 + ,89 + ,538 + ,10 + ,379 + ,284 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,242 + ,0 + ,273 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,265 + ,0 + ,121 + ,1 + ,10579 + ,92 + ,92 + ,110309 + ,375195 + ,76 + ,1904 + ,17 + ,4 + ,214 + ,1 + ,474 + ,20 + ,10 + ,61977 + ,288985 + ,386 + ,4449 + ,188 + ,0 + ,348 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,178 + ,0 + ,87 + ,1 + ,3642 + ,86 + ,100 + ,90262 + ,458343 + ,42 + ,3004 + ,71 + ,22 + ,170 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,0 + ,200 + ,0 + ,472 + ,20 + ,19 + ,80045 + ,269587 + ,397 + ,3479 + ,148 + ,0 + ,158 + ,1 + ,98 + ,2 + ,2 + ,61490 + ,315236 + ,298 + ,57618 + ,1179 + ,0 + ,427 + ,1 + ,3999 + ,30 + ,21 + ,51252 + ,42754 + ,428 + ,-5242 + ,-39 + ,-58 + ,233 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,216 + ,0 + ,345 + ,1 + ,621 + ,48 + ,13 + ,60720 + ,308256 + ,341 + ,2255 + ,174 + ,0 + ,375 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,266 + ,0 + ,173 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,189 + ,0 + ,301 + ,0 + ,30 + ,4 + ,4 + ,60798 + ,313164 + ,320 + ,28291 + ,3808 + ,0 + ,368 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,291 + ,0 + ,106 + ,1 + ,746 + ,23 + ,29 + ,71561 + ,269661 + ,396 + ,3029 + ,93 + ,10 + ,208 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,228 + ,0 + ,282 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,255 + ,0 + ,249 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,260 + ,0 + ,240 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,194 + ,0 + ,52 + ,0 + ,4150 + ,150 + ,136 + ,134759 + ,518365 + ,35 + ,2122 + ,77 + ,66 + ,55 + ,0 + ,4658 + ,158 + ,114 + ,156608 + ,233773 + ,411 + ,214 + ,7 + ,58 + ,84 + ,1 + ,814 + ,18 + ,15 + ,39625 + ,301881 + ,358 + ,5660 + ,125 + ,23 + ,143 + ,1 + ,1002 + ,15 + ,11 + ,56750 + ,298568 + ,364 + ,6571 + ,98 + ,0 + ,205 + ,0 + ,496 + ,17 + ,12 + ,87390 + ,325479 + ,128 + ,7381 + ,253 + ,0 + ,141 + ,1 + ,389 + ,14 + ,8 + ,58990 + ,325506 + ,127 + ,8965 + ,322 + ,0 + ,12 + ,1 + ,12679 + ,110 + ,122 + ,48020 + ,984885 + ,13 + ,7135 + ,62 + ,1871 + ,177 + ,1 + ,400 + ,152 + ,113 + ,60720 + ,313267 + ,319 + ,745 + ,283 + ,0 + ,238 + ,1 + ,53 + ,6 + ,6 + ,60349 + ,315793 + ,163 + ,19299 + ,2197 + ,0 + ,331 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,278 + ,0 + ,426 + ,0 + ,5109 + ,97 + ,82 + ,67038 + ,215362 + ,417 + ,158 + ,3 + ,-31 + ,104 + ,1 + ,576 + ,18 + ,8 + ,113761 + ,314073 + ,313 + ,6337 + ,198 + ,12 + ,409 + ,1 + ,438 + ,11 + ,9 + ,58320 + ,298096 + ,366 + ,8918 + ,224 + ,-5 + ,292 + ,1 + ,165 + ,24 + ,12 + ,62841 + ,325176 + ,130 + ,5216 + ,759 + ,0 + ,40 + ,1 + ,4069 + ,145 + ,131 + ,79804 + ,207393 + ,420 + ,51 + ,2 + ,119 + ,185 + ,0 + ,23 + ,5 + ,5 + ,62555 + ,314806 + ,306 + ,22961 + ,4974 + ,0 + ,59 + ,1 + ,1285 + ,107 + ,51 + ,99489 + ,341340 + ,106 + ,1321 + ,110 + ,49 + ,82 + ,0 + ,4677 + ,75 + ,66 + ,90131 + ,426280 + ,52 + ,3017 + ,48 + ,24 + ,15 + ,0 + ,24811 + ,505 + ,280 + ,95350 + ,929118 + ,15 + ,1444 + ,29 + ,822 + ,253 + ,1 + ,167 + ,2 + ,1 + ,64245 + ,307322 + ,345 + ,53661 + ,642 + ,0 + ,342 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,187 + ,0 + ,45 + ,1 + ,4628 + ,62 + ,39 + ,69159 + ,387475 + ,64 + ,3024 + ,41 + ,85 + ,189 + ,1 + ,226 + ,5 + ,5 + ,65745 + ,291787 + ,382 + ,18357 + ,406 + ,0 + ,76 + ,0 + ,1765 + ,20 + ,20 + ,77623 + ,247060 + ,404 + ,2353 + ,27 + ,28 + ,145 + ,1 + ,460 + ,15 + ,15 + ,63346 + ,329784 + ,122 + ,8652 + ,282 + ,0 + ,144 + ,1 + ,36 + ,3 + ,4 + ,60894 + ,315834 + ,162 + ,38611 + ,3215 + ,0 + ,408 + ,1 + ,989 + ,11 + ,12 + ,58930 + ,304555 + ,350 + ,9505 + ,106 + ,-5 + ,96 + ,0 + ,3055 + ,43 + ,30 + ,60247 + ,376641 + ,75 + ,4108 + ,58 + ,14 + ,265 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,288 + ,0 + ,161 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,212 + ,0 + ,246 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,236 + ,0 + ,237 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,202 + ,0 + ,230 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,214 + ,0 + ,66 + ,1 + ,13253 + ,947 + ,251 + ,90829 + ,357760 + ,90 + ,167 + ,12 + ,35 + ,244 + ,1 + ,24 + ,8 + ,7 + ,59818 + ,315637 + ,168 + ,14455 + ,4772 + ,0 + ,153 + ,1 + ,0 + ,0 + ,5 + ,60720 + ,315380 + ,208 + ,0 + ,199 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,261 + ,0 + ,222 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,198 + ,0 + ,363 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,289 + ,0 + ,387 + ,1 + ,131 + ,6 + ,3 + ,59661 + ,327071 + ,123 + ,21179 + ,969 + ,0 + ,85 + ,1 + ,514 + ,15 + ,15 + ,102725 + ,377516 + ,73 + ,11834 + ,346 + ,23 + ,103 + ,1 + ,3366 + ,80 + ,55 + ,108479 + ,299243 + ,362 + ,1241 + ,29 + ,12 + ,60 + ,0 + ,9327 + ,101 + ,85 + ,87419 + ,387699 + ,63 + ,1858 + ,20 + ,44 + ,309 + ,1 + ,384 + ,22 + ,17 + ,54683 + ,309836 + ,333 + ,4993 + ,286 + ,0 + ,20 + ,1 + ,17821 + ,373 + ,277 + ,122844 + ,444477 + ,46 + ,655 + ,14 + ,643 + ,391 + ,1 + ,397 + ,10 + ,10 + ,62710 + ,322327 + ,136 + ,12233 + ,308 + ,-1 + ,139 + ,1 + ,897 + ,5 + ,4 + ,60720 + ,314913 + ,303 + ,22983 + ,128 + ,0 + ,370 + ,1 + ,218 + ,21 + ,7 + ,60720 + ,315553 + ,170 + ,5503 + ,531 + ,0 + ,31 + ,0 + ,3369 + ,28 + ,26 + ,87161 + ,688779 + ,23 + ,17456 + ,145 + ,217 + ,109 + ,1 + ,5702 + ,61 + ,47 + ,101481 + ,321896 + ,138 + ,1998 + ,21 + ,10 + ,41 + ,1 + ,8636 + ,855 + ,203 + ,128294 + ,301607 + ,359 + ,119 + ,12 + ,108 + ,215 + ,1 + ,534 + ,7 + ,5 + ,62620 + ,304485 + ,351 + ,14926 + ,196 + ,0 + ,332 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,0 + ,310 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,227 + ,0 + ,410 + ,1 + ,726 + ,3 + ,3 + ,69980 + ,231861 + ,412 + ,10620 + ,44 + ,-5 + ,81 + ,1 + ,1380 + ,26 + ,29 + ,60982 + ,347385 + ,98 + ,5669 + ,107 + ,25 + ,365 + ,1 + ,180 + ,11 + ,9 + ,59635 + ,316386 + ,158 + ,10581 + ,647 + ,0 + ,53 + ,0 + ,7285 + ,115 + ,97 + ,188873 + ,491303 + ,39 + ,2533 + ,40 + ,59 + ,67 + ,1 + ,880 + ,74 + ,15 + ,80791 + ,261216 + ,399 + ,827 + ,70 + ,34 + ,335 + ,1 + ,1 + ,2 + ,1 + ,60727 + ,315388 + ,175 + ,57694 + ,117743 + ,0 + ,380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,271 + ,0 + ,334 + ,1 + ,96 + ,5 + ,4 + ,60379 + ,313729 + ,316 + ,22746 + ,1188 + ,0 + ,95 + ,1 + ,1889 + ,45 + ,43 + ,37527 + ,358649 + ,88 + ,3526 + ,84 + ,14 + ,8 + ,0 + ,45187 + ,353 + ,267 + ,234817 + ,1926517 + ,4 + ,4891 + ,38 + ,2854 + ,221 + ,1 + ,288 + ,4 + ,3 + ,60510 + ,296656 + ,373 + ,24164 + ,335 + ,0 + ,419 + ,1 + ,1270 + ,26 + ,24 + ,69206 + ,275311 + ,393 + ,2897 + ,59 + ,-10 + ,424 + ,1 + ,6526 + ,26 + ,20 + ,55830 + ,-42143 + ,430 + ,-9313 + ,-37 + ,-28 + ,327 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,252 + ,0 + ,263 + ,1 + ,226 + ,12 + ,13 + ,72835 + ,343929 + ,102 + ,11994 + ,636 + ,0 + ,394 + ,1 + ,694 + ,8 + ,6 + ,68060 + ,367655 + ,83 + ,20957 + ,241 + ,-2 + ,286 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,241 + ,0 + ,283 + ,1 + ,249 + ,72 + ,26 + ,56726 + ,313491 + ,317 + ,1576 + ,455 + ,0) + ,dim=c(11 + ,431) + ,dimnames=list(c('Rank' + ,'Group' + ,'Costs' + ,'Trades' + ,'Orders' + ,'Dividends' + ,'Wealth' + ,'WRank' + ,'Profit/Trades' + ,'Profit/Cost' + ,'CCScore ') + ,1:431)) > y <- array(NA,dim=c(11,431),dimnames=list(c('Rank','Group','Costs','Trades','Orders','Dividends','Wealth','WRank','Profit/Trades','Profit/Cost','CCScore '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > 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 = '11' > 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 Orders Dividends Wealth WRank 1 0 350 1 1595 17 10 60720 319369 143 2 39 61 0 5565 47 43 94355 493408 38 3 0 287 0 601 6 5 60720 319210 144 4 0 268 1 188 11 9 77655 381180 67 5 305 25 1 7146 105 83 134028 297978 367 6 -9 418 0 1135 17 20 62285 290476 384 7 -1 392 1 450 8 6 59325 292136 380 8 0 131 1 34 8 7 60630 314353 309 9 0 316 1 133 14 4 65990 339445 109 10 0 347 1 119 6 3 59118 303677 354 11 33 68 0 2053 53 44 100423 397144 60 12 32 70 0 4036 63 85 100269 424898 53 13 1 147 1 0 0 0 60720 315380 203 14 1 4655 393 55 60720 341570 105 360 30 15 1 131 11 15 61808 308989 338 9908 834 16 1 1766 73 46 60438 305959 349 1451 60 17 1 312 14 3 58598 318690 147 8478 381 18 1 448 40 7 59781 323361 132 3084 275 19 1 115 13 10 60945 318903 145 9146 1037 20 1 60 10 7 61124 314049 314 11405 1916 21 35 0 0 0 60720 315380 180 0 258 22 0 26 47705 298466 365 2813 271 0 343 23 121173 0 60720 315380 183 0 398 0 1442 24 57530 438493 49 2021 165 -3 115 0 1389 25 62041 378049 71 19783 128 6 188 1 149 26 60720 313332 318 22666 761 0 379 1 2212 27 136996 319864 142 8562 54 0 48 0 7489 28 60720 430866 50 717 31 70 366 1 0 29 115 315380 287 0 117 1 402 26 8 30 387 5106 330 5 412 0 7419 29 21 31 219 2999 12 -5 157 1 0 0 0 32 430 0 264 1 307 15 10 58856 331955 33 288 0 54 0 1134 162 63 146216 527021 34 22 58 64 0 7561 87 77 82425 364304 35 18 37 113 1 5131 71 44 86111 292154 36 12458 9 148 1 9 6 7 60835 314987 37 32 0 425 1 2264 24 20 55174 272713 38 29 -31 7 1 40949 481 263 129352 1398893 39 42 2865 197 0 4980 91 89 113521 409642 40 844 0 97 1 272 8 9 112995 429112 41 242 14 168 1 757 19 18 45689 382712 42 123 0 69 1 1424 94 73 131116 374943 43 184 32 127 1 0 0 0 60720 315380 44 401 0 378 52 16 71873 297413 369 1873 45 371 1 1061 5 3 70415 304252 353 20850 46 203 1 0 0 0 60720 315380 290 0 47 3689 9 9 61938 309596 334 12177 541 0 48 567 25 25 101792 377305 74 7092 48 74 49 27664 11 14 103345 386212 65 16928 329 20 50 1686 282 236 74996 657954 24 1624 17 247 51 1164 40 28 68696 368186 81 4205 100 11 52 15824 17 15 57320 269753 395 4103 60 -6 53 6575 64 56 63838 585715 30 6027 24 284 54 37597 159 63 151352 1071292 11 5480 133 819 55 314 248 234 37884 992426 12 3195 21 652 56 1932 5 4 60720 315371 294 23074 367 0 57 5692 36 20 72954 249898 403 1386 26 10 58 0 72 68 64239 357312 91 2185 28 22 59 0 0 0 60720 315380 257 0 322 1 60 364839 60720 315380 177 0 94 0 6533 75 61 230621 85 2198 25 16 400 0 2126 138 62 315877 413 222 14 -3 150 1 14 1 63 125390 161 115877 8041 0 107 0 3192 60 64 314882 427 -1243 -23 10 166 1 47 3 65 370837 305 38294 2432 0 47 0 4891 217 66 330068 79 787 35 73 182 0 531 10 67 324385 120 13007 245 0 245 1 2332 4 68 315380 131 31096 53 0 330 1 0 0 69 11409 277 0 10 0 16275 109 129 47818 70 1115 76 1974 429 1 9252 48 32 79510 71 13337 6 -129 2 0 29790 309 444 81767 72 0 138 28328 235 1 0 0 0 60720 73 0 187 0 2300 29 23 64873 352108 97 74 0 213 1 267 32 11 57640 330059 121 75 0 296 0 382 21 9 59365 340968 107 76 0 218 1 397 28 7 60646 317736 152 77 238 28 1 11920 182 73 31080 209458 419 78 3040 5 0 54660 302 315 126942 1491348 6 79 0 138 1 5 2 1 60720 314887 304 80 1 163 1 0 0 0 60720 315380 218 81 37 0 0 2 60720 315380 199 0 125 82 0 17 62920 353058 95 4137 58 1 216 83 60720 2 60720 315380 201 0 196 0 0 84 308 315380 263 0 128 1 94 7 10 85 296 16362 1212 0 223 1 3 2 1 86 357 57677 36505 0 165 1 454 51 24 87 110 2004 225 0 373 1 227 9 7 88 245 15182 603 0 324 0 0 0 0 89 469 0 251 1 206 9 9 53110 296702 90 34 0 11 1 25830 115 104 245546 1073089 91 -35 1910 62 1 1528 38 32 60326 146494 92 463 38 149 1 271 2 2 60720 325249 93 954 0 270 1 138 6 26 69817 331420 94 226 0 171 1 0 0 0 60720 315380 95 276 1 278 26 22 61404 314922 302 4420 96 383 1 282 30 12 62452 320016 141 4001 97 93 1 571 80 14 41477 280398 390 1005 98 302 1 2253 99 68 63593 452469 43 2550 99 385 1 290 2 1 58790 301164 360 50582 100 225 1 78 16 12 62700 317330 154 7333 101 275 1 20 2 1 60805 315576 169 57788 102 18073 1 0 0 0 60720 315380 284 0 103 1866 180 69 27284 -7170 429 -1151 -11 -366 104 42634 29 28 56225 322331 135 4218 66 0 105 249 156 168 157214 1629616 5 9164 34 2964 106 422 62 13 54323 292754 378 1496 373 0 107 2675 28 23 57935 318056 150 4216 280 0 108 965 79 54 59017 355178 93 1964 58 36 109 0 3 2 73490 204325 421 1442 4 -8 110 11 0 0 60720 315380 224 0 407 1 111 0 68005 317046 156 7315 188 -5 252 0 112 309560 60720 315380 240 0 257 1 365 8 113 414462 335 13695 300 0 63 0 3122 34 114 857217 54 6308 69 38 21 0 12988 119 115 697458 16 5523 51 601 38 0 5336 81 116 530670 22 6141 93 132 51 0 3160 108 117 238125 33 3062 105 66 386 1 2489 40 118 315380 408 953 15 0 325 1 0 0 119 12031 273 0 33 0 6984 45 40 80570 120 3718 78 195 91 1 2001 52 36 95551 121 1918 97 19 430 1 15236 90 86 77440 122 0 11 -158 259 0 0 0 0 60720 123 0 382 1 530 60 18 73433 317291 155 124 1 369 1 0 0 0 60720 315380 292 125 0 35624 223 244 157278 306275 347 477 3 126 1 1383 11 9 73221 317892 151 10717 85 127 0 875 14 12 67000 334280 112 9591 153 128 0 0 0 0 60720 315380 259 0 289 129 60720 0 60720 315380 254 0 307 0 0 130 312 315380 222 0 389 1 72 3 3 131 346 38070 1591 0 210 0 265 31 3 132 0 3450 404 194 1 0 0 0 60720 133 0 277 1 335 33 23 64175 320398 140 134 1 298 0 0 0 0 60720 315380 230 135 226 0 0 0 60720 315380 269 0 43 136 58 89 93811 501749 37 1335 67 101 71 137 0 49 27330 202055 422 35 1 31 133 138 60720 0 60720 315380 196 0 212 1 0 139 114 315380 206 0 186 0 638 14 11 140 283 9515 209 0 262 0 0 0 0 141 201 0 378 1 607 11 7 60436 322340 142 109 0 132 1 1558 63 25 55637 369448 143 69 0 279 0 1324 117 29 67440 291841 144 267 0 376 1 0 0 0 60720 315380 145 611 0 0 0 0 60720 315380 282 0 146 0 8 6 59190 296919 371 12115 159 -4 147 6 0 1 60720 315380 225 0 388 1 148 1 58620 309038 337 13630 118 0 420 1 149 13 65920 246541 405 15514 70 -12 114 0 150 14 74020 289513 385 6394 37 6 305 1 151 0 61808 344425 101 2725 394 0 346 1 152 314210 60720 315380 181 0 354 1 135 2 153 480382 311 57105 847 0 124 0 1659 29 154 315009 40 9668 169 2 136 1 316 9 155 315380 299 12779 364 0 130 1 0 0 156 10262 197 0 146 1 309 11 9 56535 157 15254 365 0 314 0 49 8 14 64107 158 0 2500 0 241 1 0 0 0 60720 159 163 34 0 4519 49 42 102129 597793 29 160 1 384 1 0 0 0 60720 315380 267 161 1 837 69 21 61262 315688 165 1677 138 162 1 5119 49 38 39039 378525 69 3643 35 163 1 1280 117 52 69465 312378 326 960 88 164 0 2564 22 19 130140 403560 59 9253 79 165 0 2045 20 16 97890 510834 36 15542 152 166 1 2234 16 11 77200 214215 418 888 6 167 0 975 23 19 90534 235133 410 1528 36 168 1 1136 32 26 48522 343613 103 4488 126 169 1 453 21 19 81125 365959 84 7903 366 170 14 0 0 0 60720 315380 210 0 227 171 17 9 60720 314551 307 8182 1874 0 234 172 0 13 61686 303230 355 6072 281 0 247 173 121920 0 60720 315380 233 0 49 1 4901 174 103960 469107 41 3204 55 67 88 1 540 175 60720 354228 94 11016 286 22 329 1 0 176 213 315380 253 0 160 1 0 0 0 177 306 0 271 0 0 0 0 60720 315380 178 172 1 2 9 12 60735 315394 174 12822 179 402 1 36 4 3 61564 312412 325 28103 180 4 1 776 9 8 64230 333505 113 14834 181 156 1 84738 588 428 -26007 223193 415 39 182 3 1 0 0 0 60720 315380 217 0 183 529 4 8 60761 315656 167 28914 37920 0 184 0 9 7 63870 296261 374 10696 182 -3 185 5 0 0 60720 315380 264 0 159 1 186 19 60845 336425 111 19489 337 0 181 1 187 0 71642 359335 87 4686 164 0 278 1 188 308636 60720 315380 232 0 89 1 2099 64 189 158492 339 1697 52 20 112 0 3437 47 190 315380 425 -883 -12 9 201 0 0 0 191 0 243 0 344 1 0 0 0 60720 192 446 22 0 22330 84 67 79801 711969 20 193 0 291 0 0 0 0 60720 315380 258 194 27 0 0 0 60720 315380 221 0 164 195 0 7 60830 306268 348 3936 220 0 358 196 88590 0 60720 315380 279 0 118 1 2239 197 82903 442882 47 11566 108 5 92 0 2949 198 60720 378509 70 4354 61 18 300 1 0 199 99 315380 223 0 312 1 365 83 32 200 310 1766 402 0 183 1 2461 57 25 201 17 2005 46 0 18 0 21950 519 265 202 416 1266 30 745 126 0 3294 23 23 203 295 748 5 1 193 1 141 6 7 204 342 19228 819 0 390 1 572 16 8 205 21 6746 189 -1 23 1 13326 102 103 206 423 4925 38 436 50 1 2284 33 41 207 160 -167 -2 66 206 0 10 2 5 208 207 58077 11639 0 162 1 0 0 0 209 92 0 98 0 1414 198 22 70184 330546 210 98 13 75 1 1975 35 34 118881 394510 211 2613 30 134 1 43 3 20 60887 312846 212 337 0 180 1 0 0 0 60720 315380 213 341 1 844 78 20 60925 296139 375 1233 214 154 1 304 15 12 62969 295580 376 6372 215 36 1 458 11 7 58625 297765 368 8888 216 192 1 18562 155 211 102313 377934 72 1148 217 7123 1 0 0 0 60720 315380 179 0 218 622 109 102 288170 638830 28 4026 62 769 219 174 73 38 73007 304376 352 1430 168 0 220 2220 1 3 64820 307424 344 107424 618 0 221 121 22 12 301670 644190 26 20190 200 154 222 11819 24 11 56178 295370 377 3974 788 0 223 0 85 70 106113 574339 31 4404 32 199 224 19 0 0 60720 315380 235 0 274 1 225 11 60798 310201 332 9183 883 0 406 1 226 35 70694 327007 124 9770 107 -5 116 1 227 0 56364 343466 104 3679 95 6 320 1 228 315380 60720 315380 274 0 367 1 0 0 229 29525 248 0 288 1 30 4 3 62045 230 7523 3912 0 414 0 3310 33 32 75230 231 11431 75 -6 74 1 554 11 9 79285 232 0 227 30 250 0 0 0 0 60720 233 0 224 1 468 21 20 60720 312161 328 234 -3 397 1 4917 42 26 52811 243650 406 235 -8 417 1 3256 66 23 35250 407159 57 236 1 323 1 0 0 0 60720 315380 185 237 1 125 13 10 59734 317698 153 9054 944 238 2 0 0 0 60720 315380 195 0 326 239 0 1 60722 312502 324 56251 5154 0 261 240 78780 0 60720 315380 262 0 357 1 514 241 60720 322378 133 2353 238 0 313 0 0 242 270 315380 226 0 360 1 0 0 0 243 333 0 272 0 0 0 0 60720 315380 244 0 1 0 0 0 60720 315380 276 0 245 0 0 0 60720 315380 272 0 362 1 246 640273 60720 315380 268 0 37 1 10327 91 247 345783 27 4838 43 135 42 1 13718 61 248 652925 100 2390 11 104 57 1 3748 9 249 439798 25 50325 121 53 30 1 14416 361 250 278990 48 664 17 220 404 0 1526 25 251 339836 391 3160 52 -4 80 0 666 171 252 240897 108 818 210 25 422 1 2844 67 253 315380 407 610 14 -15 297 1 0 0 254 5714 229 0 318 1 368 17 20 61167 255 0 264 0 175 1 0 0 0 60720 256 0 285 1 333 13 9 70811 331323 118 257 0 339 1 26 1 1 60896 313880 315 258 1 243 1 0 0 0 60720 315380 237 259 61 0 0 0 60720 315380 193 0 411 260 3 57 69863 309422 336 1794 84 -5 155 261 97 2 60938 315245 297 38415 5674 0 56 262 27 77 61348 405972 58 2123 86 54 204 263 1 12 50804 300962 361 3739 498 0 178 264 14 0 60745 316176 159 116176 1646 0 295 265 13 10 59506 302409 356 7315 1932 0 396 266 7 8 58456 283587 389 6430 149 -3 393 267 5 8 60950 263276 398 9039 102 -2 135 268 24 4 60720 312075 329 22415 174 0 364 269 0 28 61600 308336 340 4514 61 0 179 270 63915 0 60720 315380 209 0 231 1 317 271 60719 298700 363 10967 311 0 319 1 1 272 59500 315372 293 19229 84213 0 191 1 275 273 60720 318745 146 29686 431 0 321 1 0 274 56 315380 246 0 111 1 936 30 22 275 18 6963 223 9 29 1 8568 37 34 276 431 15856 68 223 423 0 11528 70 68 277 250 -4047 -25 -25 202 0 0 0 0 278 164 0 232 1 738 10 9 60720 321376 279 405 0 266 1 0 0 0 60720 315380 280 176 1 592 7 7 64270 276898 392 10985 281 236 1 126 10 8 60951 315547 171 11555 282 142 1 2 1 0 60743 315487 172 115487 283 18014 1 0 0 0 60720 315380 220 0 284 0 239 393 -1710 1405225 8 5043 67 943 285 0 0 0 60720 315380 256 0 328 0 286 983660 60720 315380 249 0 9 0 37704 452 287 318574 14 1734 21 2167 207 0 63 4 288 310768 148 29643 1891 0 73 1 1431 34 289 312887 330 3258 77 30 242 1 94 3 290 312339 321 37629 1200 0 229 1 192 4 291 314964 327 28085 586 0 317 1 32 1 292 379983 301 114964 3603 0 44 1 6869 105 293 315380 68 1714 26 90 349 1 0 0 294 0 176 0 140 1 0 0 0 60720 295 0 361 1 1 4 4 60722 315398 173 296 1 203 0 0 0 0 60720 315380 247 297 0 2328 94 57 82732 253588 400 570 23 298 1 209 13 19 64016 316647 157 8973 558 299 1 28 2 3 60890 315688 166 57844 4148 300 1 176 28 8 68136 310670 331 3952 629 301 0 1920 15 11 79420 165404 424 -2306 -18 302 1 87550 458 412 153198 4111912 3 8541 45 303 0 520 10 7 58650 291650 383 9165 176 304 2 0 0 0 60720 315380 231 0 421 305 8 1 69770 253468 402 26734 53 -14 217 306 12 6 60831 315688 164 14461 7687 0 353 307 66 10 59595 325699 126 10475 214 0 58 308 0 47 87720 446211 45 3730 46 49 359 309 114768 0 60720 315380 251 0 77 1 1012 310 60720 368078 82 3576 166 26 228 1 0 311 96 315380 215 0 83 1 876 59 33 312 1 2591 174 24 1 0 162556 1081 807 313 414 5626 37 56289 14 1 43556 247 190 314 388 110 1 929 395 0 3425 12 9 315 409 6993 25 -3 78 1 810 43 11 316 184 855 45 26 190 1 0 0 0 317 39 0 311 0 0 0 0 60720 315380 318 101 0 2365 79 83 191778 550608 32 4438 319 256 0 1261 33 31 76114 307528 343 3258 320 1585 1 0 0 0 60720 315380 238 0 321 16189 107 112 93099 355864 92 1457 98 2 322 0 295 290 116384 358589 89 538 10 379 323 0 0 0 60720 315380 242 0 273 0 324 375195 60720 315380 265 0 121 1 10579 92 325 288985 76 1904 17 4 214 1 474 20 326 315380 386 4449 188 0 348 1 0 0 327 3004 178 0 87 1 3642 86 100 90262 328 0 71 22 170 1 0 0 0 60720 329 0 200 0 472 20 19 80045 269587 397 330 0 158 1 98 2 2 61490 315236 298 331 -58 427 1 3999 30 21 51252 42754 428 332 1 233 1 0 0 0 60720 315380 216 333 1 621 48 13 60720 308256 341 2255 174 334 0 0 0 0 60720 315380 266 0 173 335 60798 0 60720 315380 189 0 301 0 30 336 60720 313164 320 28291 3808 0 368 1 0 337 396 315380 291 0 106 1 746 23 29 338 228 3029 93 10 208 0 0 0 0 339 249 0 282 1 0 0 0 60720 315380 340 0 0 0 0 0 60720 315380 260 0 341 136 0 0 60720 315380 194 0 52 0 342 114 134759 518365 35 2122 77 66 55 0 343 15 156608 233773 411 214 7 58 84 1 344 11 39625 301881 358 5660 125 23 143 1 345 12 56750 298568 364 6571 98 0 205 0 346 8 87390 325479 128 7381 253 0 141 1 347 122 58990 325506 127 8965 322 0 12 1 348 113 48020 984885 13 7135 62 1871 177 1 349 6 60720 313267 319 745 283 0 238 1 350 0 60349 315793 163 19299 2197 0 331 1 351 215362 60720 315380 278 0 426 0 5109 97 352 314073 417 158 3 -31 104 1 576 18 353 298096 313 6337 198 12 409 1 438 11 354 325176 366 8918 224 -5 292 1 165 24 355 207393 130 5216 759 0 40 1 4069 145 356 314806 420 51 2 119 185 0 23 5 357 341340 306 22961 4974 0 59 1 1285 107 358 426280 106 1321 110 49 82 0 4677 75 359 929118 52 3017 48 24 15 0 24811 505 360 307322 15 1444 29 822 253 1 167 2 361 315380 345 53661 642 0 342 1 0 0 362 3024 187 0 45 1 4628 62 39 69159 363 18357 41 85 189 1 226 5 5 65745 364 2353 406 0 76 0 1765 20 20 77623 365 8652 27 28 145 1 460 15 15 63346 366 38611 282 0 144 1 36 3 4 60894 367 9505 3215 0 408 1 989 11 12 58930 368 4108 106 -5 96 0 3055 43 30 60247 369 0 58 14 265 0 0 0 0 60720 370 0 161 1 0 0 0 60720 315380 212 371 0 0 0 0 60720 315380 236 0 237 372 60720 0 60720 315380 202 0 230 1 0 373 90 315380 214 0 66 1 13253 947 251 374 168 167 12 35 244 1 24 8 7 375 208 14455 4772 0 153 1 0 0 5 376 222 0 199 0 0 0 0 60720 315380 377 0 1 0 0 0 60720 315380 198 0 378 3 0 0 60720 315380 289 0 387 1 379 15 59661 327071 123 21179 969 0 85 1 380 55 102725 377516 73 11834 346 23 103 1 381 85 108479 299243 362 1241 29 12 60 0 382 17 87419 387699 63 1858 20 44 309 1 383 277 54683 309836 333 4993 286 0 20 1 384 10 122844 444477 46 655 14 643 391 1 385 4 62710 322327 136 12233 308 -1 139 1 386 7 60720 314913 303 22983 128 0 370 1 387 26 60720 315553 170 5503 531 0 31 0 388 47 87161 688779 23 17456 145 217 109 1 389 203 101481 321896 138 1998 21 10 41 1 390 5 128294 301607 359 119 12 108 215 1 391 0 62620 304485 351 14926 196 0 332 1 392 315380 60720 315380 275 0 310 0 0 0 393 10620 227 0 410 1 726 3 3 69980 394 5669 44 -5 81 1 1380 26 29 60982 395 10581 107 25 365 1 180 11 9 59635 396 2533 647 0 53 0 7285 115 97 188873 397 827 40 59 67 1 880 74 15 80791 398 57694 70 34 335 1 1 2 1 60727 399 0 117743 0 380 1 0 0 0 60720 400 0 334 1 96 5 4 60379 313729 316 401 14 95 1 1889 45 43 37527 358649 88 402 2854 8 0 45187 353 267 234817 1926517 4 403 0 221 1 288 4 3 60510 296656 373 404 -10 419 1 1270 26 24 69206 275311 393 405 -28 424 1 6526 26 20 55830 -42143 430 406 1 327 0 0 0 0 60720 315380 252 407 1 226 12 13 72835 343929 102 11994 636 408 1 694 8 6 68060 367655 83 20957 241 409 72 0 0 0 60720 315380 241 0 283 410 17 26 56726 313491 317 1576 455 0 350 411 47 10 60720 319369 143 7022 75 0 61 412 6 43 94355 493408 38 6243 53 39 287 413 11 5 60720 319210 144 19868 198 0 268 414 105 9 77655 381180 67 16471 964 0 25 415 17 83 134028 297978 367 933 14 305 418 416 8 20 62285 290476 384 5322 80 -9 392 417 8 6 59325 292136 380 11517 205 -1 131 418 14 7 60630 314353 309 14294 3340 0 316 419 6 4 65990 339445 109 9960 1052 0 347 420 53 3 59118 303677 354 17280 872 0 68 421 63 44 100423 397144 60 3720 96 33 70 422 0 85 100269 424898 53 3570 56 32 147 423 60720 0 60720 315380 203 0 169 1 4655 424 61808 341570 105 360 30 0 248 1 131 425 60438 308989 338 9908 834 0 152 1 1766 426 58598 305959 349 1451 60 0 351 1 312 427 59781 318690 147 8478 381 0 372 1 448 428 60945 323361 132 3084 275 0 137 1 115 429 61124 318903 145 9146 1037 0 352 1 60 430 60720 314049 314 11405 1916 0 338 1 0 431 365 315380 180 0 258 1 364 35 26 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 239 24 9 6 25 5 6 26 14 10 27 322 149 28 0 0 29 84990 332743 30 63255 286963 31 60720 315380 32 116 8797 33 34 2019 34 86 1889 35 379 1298 36 300 19165 37 394 3030 38 9 2493 39 55 2304 40 51 28639 41 66 9616 42 77 1861 43 191 0 44 258 0 45 98 -4 46 174 1 47 46 0 48 90 1 49 27 0 50 105 1 51 413 0 52 26 1 53 16 1 54 19 0 55 336 1 56 110 0 57 86 0 58 294 1 59 0 0 60 29 126910 61 41 21001 62 0 61016 63 48 59831 64 1 60720 65 49 210568 66 11 71641 67 2 72680 68 0 60720 69 1443586 7 70 253537 401 71 4321023 2 72 315380 204 73 5245 66 74 4064 487 75 6713 369 76 4205 297 77 52 1 78 4276 24 79 57443 24814 80 0 219 81 0 2647 82 1 0 83 0 0 84 60793 314533 85 60698 315354 86 63261 302187 87 76644 336639 88 60720 315380 89 372 10745 90 10 7592 91 426 -1408 92 129 62625 93 117 21903 94 211 0 95 413 0 96 425 0 97 141 0 98 112 17 99 348 0 100 1511 0 101 5924 0 102 431 1 103 220 1 104 6 0 105 290 1 106 129 1 107 65 1 108 416 1 109 299 1 110 621 16 111 0 0 112 7 54820 113 25 94670 114 73 82340 115 80 112477 116 90 108094 117 29 67804 118 0 60720 119 741409 19 120 393343 62 121 372631 78 122 315380 239 123 1955 221 124 0 428 125 -70 100 126 13 122 127 3 198 128 0 0 129 0 0 130 60398 314210 131 60720 306948 132 315380 186 133 3648 359 134 0 356 135 1 4525 136 0 3045 137 1 0 138 0 0 139 60370 333210 140 60720 315380 141 134 11122 142 80 2690 143 381 785 144 280 0 145 403 1 146 303 1 147 923 8 148 661 3 149 2397 14 150 366 53 151 0 0 152 1 62065 153 25 107577 154 4 60505 155 0 60720 156 312878 322 157 322031 137 158 315380 200 159 8118 88 160 0 120 161 4 72 162 31 79 163 26 99 164 13 108 165 10 119 166 5 102 167 13 209 168 0 293 169 0 211 170 1 61 171 1 368 172 0 0 173 84 58 174 14 13 175 0 0 176 60720 315380 177 286 0 178 63057 0 179 3146 0 180 172 -4 181 0 4145 182 315 1 183 399 1 184 381 1 185 405 7 186 972 34 187 0 0 188 45 106611 189 33 48022 190 0 60720 191 315380 182 192 6095 23 193 0 269 194 1 483 195 1 0 196 21 17 197 41 40 198 0 9 199 87192 346611 200 55792 314289 201 114337 856956 202 75832 217193 203 61630 315366 204 58580 307930 205 165548 702380 206 76403 194493 207 61656 316155 208 60720 315380 209 119 659 210 61 5557 211 323 37615 212 190 0 213 114 0 214 314 0 215 213 0 216 10 143 217 17 1 218 254 1 219 374 1 220 35 1 221 377 1 222 32 0 223 260 1 224 125 12 225 1182 13 226 1503 39 227 0 0 228 0 60720 229 318098 149 230 448243 44 231 325738 125 232 315380 234 233 5341 239 234 1039 9 235 3139 64 236 0 340 237 0 239 238 1 22 239 0 0 240 52 15 241 0 0 242 60720 315380 243 285 0 244 355 1 245 0 0 246 74 88577 247 56 96448 248 6 50350 249 139 49857 250 21 69351 251 13 117869 252 74 72683 253 0 60720 254 297141 370 255 315380 205 256 10102 394 257 113880 4328 258 0 167 259 1 1303 260 1 20 261 1 2384 262 0 203 263 0 71 264 1 53 265 1 562 266 1 622 267 1 645 268 0 1763 269 1 0 270 9 7 271 6 5 272 4 7 273 0 0 274 67939 408881 275 32168 786690 276 -14545 -83265 277 60720 315380 278 139 12138 279 281 0 280 130 -4 281 921 0 282 52257 0 283 13 0 284 281 0 285 0 0 286 228 60448 287 2 65688 288 35 106885 289 3 61360 290 2 65276 291 0 59988 292 88 117520 293 0 60720 294 315380 188 295 28849 172235 296 0 255 297 0 308 298 0 195 299 0 280 300 0 415 301 -6 3 302 17936 151 303 0 304 304 1 1013 305 1 15 306 1 587 307 1 5371 308 1 0 309 47 29 310 0 0 311 138971 352850 312 213118 6282154 313 32648 227132 314 83620 283910 315 74015 236761 316 60720 315380 317 244 0 318 148 120 319 85 13 320 123 1 321 24 1 322 284 1 323 0 0 324 92 110309 325 10 61977 326 0 60720 327 458343 42 328 315380 192 329 3479 148 330 57618 1179 331 -5242 -39 332 0 345 333 0 375 334 1 0 335 4 4 336 0 0 337 71561 269661 338 60720 315380 339 255 0 340 240 1 341 4150 150 342 4658 158 343 814 18 344 1002 15 345 496 17 346 389 14 347 12679 110 348 400 152 349 53 6 350 0 0 351 82 67038 352 8 113761 353 9 58320 354 12 62841 355 131 79804 356 5 62555 357 51 99489 358 66 90131 359 280 95350 360 1 64245 361 0 60720 362 387475 64 363 291787 382 364 247060 404 365 329784 122 366 315834 162 367 304555 350 368 376641 75 369 315380 288 370 0 246 371 1 0 372 0 0 373 90829 357760 374 59818 315637 375 60720 315380 376 261 0 377 363 1 378 131 6 379 514 15 380 3366 80 381 9327 101 382 384 22 383 17821 373 384 397 10 385 897 5 386 218 21 387 3369 28 388 5702 61 389 8636 855 390 534 7 391 0 0 392 0 60720 393 231861 412 394 347385 98 395 316386 158 396 491303 39 397 261216 399 398 315388 175 399 315380 271 400 22746 1188 401 3526 84 402 4891 38 403 24164 335 404 2897 59 405 -9313 -37 406 0 263 407 0 394 408 -2 286 409 1 249 410 1 1595 411 0 5565 412 0 601 413 1 188 414 1 7146 415 0 1135 416 1 450 417 1 34 418 1 133 419 1 119 420 0 2053 421 0 4036 422 1 0 423 393 55 424 11 15 425 73 46 426 14 3 427 40 7 428 13 10 429 10 7 430 0 0 431 47705 298466 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Rank Group Costs 1.125e+05 -1.981e-01 5.507e-03 -2.735e-01 Trades Orders Dividends Wealth -1.375e-01 -7.563e-02 -3.408e-01 -1.215e-01 WRank `Profit/Trades` `Profit/Cost` -2.103e-01 -7.038e-02 -5.011e-03 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -107203 -76101 -38779 12429 886489 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.125e+05 1.081e+04 10.409 < 2e-16 *** Rank -1.981e-01 7.110e-02 -2.786 0.005581 ** Group 5.507e-03 5.921e-02 0.093 0.925933 Costs -2.735e-01 6.319e-02 -4.329 1.88e-05 *** Trades -1.375e-01 4.037e-02 -3.405 0.000726 *** Orders -7.563e-02 2.899e-02 -2.609 0.009403 ** Dividends -3.408e-01 7.638e-02 -4.462 1.04e-05 *** Wealth -1.215e-01 4.128e-02 -2.944 0.003416 ** WRank -2.103e-01 5.387e-02 -3.905 0.000110 *** `Profit/Trades` -7.037e-02 2.668e-02 -2.638 0.008652 ** `Profit/Cost` -5.011e-03 1.970e-02 -0.254 0.799327 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 128800 on 420 degrees of freedom Multiple R-squared: 0.1322, Adjusted R-squared: 0.1115 F-statistic: 6.397 on 10 and 420 DF, p-value: 3.461e-09 > 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.688394e-12 7.376788e-12 1.000000e+00 [2,] 1.733192e-15 3.466385e-15 1.000000e+00 [3,] 5.276523e-19 1.055305e-18 1.000000e+00 [4,] 1.116093e-22 2.232186e-22 1.000000e+00 [5,] 2.065889e-26 4.131779e-26 1.000000e+00 [6,] 3.602514e-30 7.205029e-30 1.000000e+00 [7,] 5.827073e-34 1.165415e-33 1.000000e+00 [8,] 1.151468e-37 2.302935e-37 1.000000e+00 [9,] 4.326492e-41 8.652983e-41 1.000000e+00 [10,] 5.491528e-21 1.098306e-20 1.000000e+00 [11,] 2.783435e-22 5.566870e-22 1.000000e+00 [12,] 4.717505e-15 9.435011e-15 1.000000e+00 [13,] 2.518676e-14 5.037352e-14 1.000000e+00 [14,] 5.688787e-15 1.137757e-14 1.000000e+00 [15,] 6.198193e-16 1.239639e-15 1.000000e+00 [16,] 5.176667e-17 1.035333e-16 1.000000e+00 [17,] 1.375800e-17 2.751599e-17 1.000000e+00 [18,] 1.273658e-18 2.547316e-18 1.000000e+00 [19,] 2.333264e-14 4.666528e-14 1.000000e+00 [20,] 5.607291e-15 1.121458e-14 1.000000e+00 [21,] 8.837687e-16 1.767537e-15 1.000000e+00 [22,] 1.487347e-16 2.974695e-16 1.000000e+00 [23,] 2.103504e-17 4.207008e-17 1.000000e+00 [24,] 4.385373e-18 8.770746e-18 1.000000e+00 [25,] 1.841784e-18 3.683569e-18 1.000000e+00 [26,] 2.712490e-19 5.424980e-19 1.000000e+00 [27,] 3.778591e-20 7.557183e-20 1.000000e+00 [28,] 5.250583e-21 1.050117e-20 1.000000e+00 [29,] 6.992542e-22 1.398508e-21 1.000000e+00 [30,] 9.513714e-23 1.902743e-22 1.000000e+00 [31,] 1.226987e-23 2.453974e-23 1.000000e+00 [32,] 1.591617e-24 3.183234e-24 1.000000e+00 [33,] 1.992949e-25 3.985897e-25 1.000000e+00 [34,] 7.437804e-26 1.487561e-25 1.000000e+00 [35,] 1.607432e-26 3.214863e-26 1.000000e+00 [36,] 8.121190e-27 1.624238e-26 1.000000e+00 [37,] 1.837952e-26 3.675903e-26 1.000000e+00 [38,] 2.675354e-27 5.350707e-27 1.000000e+00 [39,] 4.025613e-28 8.051225e-28 1.000000e+00 [40,] 6.253353e-29 1.250671e-28 1.000000e+00 [41,] 1.810591e-29 3.621183e-29 1.000000e+00 [42,] 8.650288e-30 1.730058e-29 1.000000e+00 [43,] 1.261012e-30 2.522025e-30 1.000000e+00 [44,] 1.714500e-31 3.429001e-31 1.000000e+00 [45,] 2.563186e-32 5.126372e-32 1.000000e+00 [46,] 3.744217e-33 7.488433e-33 1.000000e+00 [47,] 6.842242e-34 1.368448e-33 1.000000e+00 [48,] 3.205729e-21 6.411457e-21 1.000000e+00 [49,] 6.975445e-15 1.395089e-14 1.000000e+00 [50,] 3.990555e-15 7.981111e-15 1.000000e+00 [51,] 8.545785e-13 1.709157e-12 1.000000e+00 [52,] 4.478064e-13 8.956128e-13 1.000000e+00 [53,] 6.534394e-12 1.306879e-11 1.000000e+00 [54,] 2.676902e-11 5.353804e-11 1.000000e+00 [55,] 7.871734e-11 1.574347e-10 1.000000e+00 [56,] 1.849966e-06 3.699932e-06 9.999982e-01 [57,] 1.907635e-06 3.815271e-06 9.999981e-01 [58,] 4.820426e-06 9.640853e-06 9.999952e-01 [59,] 6.951496e-06 1.390299e-05 9.999930e-01 [60,] 4.236270e-06 8.472540e-06 9.999958e-01 [61,] 2.582358e-06 5.164716e-06 9.999974e-01 [62,] 1.554771e-06 3.109543e-06 9.999984e-01 [63,] 9.358234e-07 1.871647e-06 9.999991e-01 [64,] 6.339984e-07 1.267997e-06 9.999994e-01 [65,] 4.391504e-06 8.783008e-06 9.999956e-01 [66,] 2.882146e-06 5.764291e-06 9.999971e-01 [67,] 1.879807e-06 3.759614e-06 9.999981e-01 [68,] 1.174792e-06 2.349585e-06 9.999988e-01 [69,] 7.713549e-07 1.542710e-06 9.999992e-01 [70,] 4.627911e-07 9.255822e-07 9.999995e-01 [71,] 3.888454e-07 7.776907e-07 9.999996e-01 [72,] 3.747272e-07 7.494543e-07 9.999996e-01 [73,] 3.965556e-07 7.931112e-07 9.999996e-01 [74,] 2.967781e-07 5.935561e-07 9.999997e-01 [75,] 2.137654e-07 4.275309e-07 9.999998e-01 [76,] 1.418771e-07 2.837542e-07 9.999999e-01 [77,] 1.351079e-07 2.702158e-07 9.999999e-01 [78,] 1.007435e-07 2.014871e-07 9.999999e-01 [79,] 6.337771e-08 1.267554e-07 9.999999e-01 [80,] 3.984466e-08 7.968933e-08 1.000000e+00 [81,] 2.540229e-08 5.080458e-08 1.000000e+00 [82,] 1.494061e-08 2.988122e-08 1.000000e+00 [83,] 8.678179e-09 1.735636e-08 1.000000e+00 [84,] 4.943773e-09 9.887546e-09 1.000000e+00 [85,] 3.095847e-09 6.191695e-09 1.000000e+00 [86,] 1.741076e-09 3.482151e-09 1.000000e+00 [87,] 9.717370e-10 1.943474e-09 1.000000e+00 [88,] 5.414763e-10 1.082953e-09 1.000000e+00 [89,] 3.044255e-10 6.088511e-10 1.000000e+00 [90,] 2.478290e-10 4.956580e-10 1.000000e+00 [91,] 1.358189e-10 2.716379e-10 1.000000e+00 [92,] 1.375844e-10 2.751688e-10 1.000000e+00 [93,] 8.598749e-11 1.719750e-10 1.000000e+00 [94,] 5.164316e-11 1.032863e-10 1.000000e+00 [95,] 3.065901e-11 6.131802e-11 1.000000e+00 [96,] 1.882414e-11 3.764827e-11 1.000000e+00 [97,] 1.120190e-11 2.240380e-11 1.000000e+00 [98,] 7.360563e-10 1.472113e-09 1.000000e+00 [99,] 7.240750e-10 1.448150e-09 1.000000e+00 [100,] 7.450263e-08 1.490053e-07 9.999999e-01 [101,] 5.541553e-02 1.108311e-01 9.445845e-01 [102,] 5.364233e-01 9.271534e-01 4.635767e-01 [103,] 7.875479e-01 4.249043e-01 2.124521e-01 [104,] 7.796291e-01 4.407419e-01 2.203709e-01 [105,] 8.014342e-01 3.971316e-01 1.985658e-01 [106,] 7.937186e-01 4.125628e-01 2.062814e-01 [107,] 7.808217e-01 4.383565e-01 2.191783e-01 [108,] 7.666996e-01 4.666008e-01 2.333004e-01 [109,] 7.540077e-01 4.919845e-01 2.459923e-01 [110,] 7.330802e-01 5.338395e-01 2.669198e-01 [111,] 7.122219e-01 5.755562e-01 2.877781e-01 [112,] 6.877260e-01 6.245481e-01 3.122740e-01 [113,] 6.645081e-01 6.709838e-01 3.354919e-01 [114,] 6.409022e-01 7.181956e-01 3.590978e-01 [115,] 6.174117e-01 7.651766e-01 3.825883e-01 [116,] 5.889579e-01 8.220842e-01 4.110421e-01 [117,] 5.662459e-01 8.675083e-01 4.337541e-01 [118,] 5.638005e-01 8.723989e-01 4.361995e-01 [119,] 5.473113e-01 9.053775e-01 4.526887e-01 [120,] 5.226303e-01 9.547394e-01 4.773697e-01 [121,] 4.984627e-01 9.969254e-01 5.015373e-01 [122,] 4.744810e-01 9.489621e-01 5.255190e-01 [123,] 4.452757e-01 8.905514e-01 5.547243e-01 [124,] 4.235338e-01 8.470677e-01 5.764662e-01 [125,] 3.952909e-01 7.905818e-01 6.047091e-01 [126,] 3.708983e-01 7.417966e-01 6.291017e-01 [127,] 3.677158e-01 7.354315e-01 6.322842e-01 [128,] 3.438152e-01 6.876304e-01 6.561848e-01 [129,] 3.193868e-01 6.387736e-01 6.806132e-01 [130,] 2.972580e-01 5.945160e-01 7.027420e-01 [131,] 2.752237e-01 5.504475e-01 7.247763e-01 [132,] 2.507305e-01 5.014611e-01 7.492695e-01 [133,] 2.327052e-01 4.654104e-01 7.672948e-01 [134,] 2.155124e-01 4.310248e-01 7.844876e-01 [135,] 2.438836e-01 4.877673e-01 7.561164e-01 [136,] 2.508666e-01 5.017331e-01 7.491334e-01 [137,] 2.528364e-01 5.056729e-01 7.471636e-01 [138,] 2.508619e-01 5.017239e-01 7.491381e-01 [139,] 2.834218e-01 5.668436e-01 7.165782e-01 [140,] 4.655945e-01 9.311889e-01 5.344055e-01 [141,] 4.978454e-01 9.956908e-01 5.021546e-01 [142,] 5.288026e-01 9.423947e-01 4.711974e-01 [143,] 5.099241e-01 9.801518e-01 4.900759e-01 [144,] 4.892660e-01 9.785319e-01 5.107340e-01 [145,] 4.714301e-01 9.428601e-01 5.285699e-01 [146,] 4.431017e-01 8.862033e-01 5.568983e-01 [147,] 4.205812e-01 8.411624e-01 5.794188e-01 [148,] 3.990024e-01 7.980047e-01 6.009976e-01 [149,] 3.781101e-01 7.562203e-01 6.218899e-01 [150,] 3.574679e-01 7.149358e-01 6.425321e-01 [151,] 3.363933e-01 6.727866e-01 6.636067e-01 [152,] 3.186686e-01 6.373373e-01 6.813314e-01 [153,] 3.022865e-01 6.045729e-01 6.977135e-01 [154,] 2.853743e-01 5.707486e-01 7.146257e-01 [155,] 2.683662e-01 5.367324e-01 7.316338e-01 [156,] 2.512386e-01 5.024771e-01 7.487614e-01 [157,] 2.359147e-01 4.718295e-01 7.640853e-01 [158,] 2.155392e-01 4.310783e-01 7.844608e-01 [159,] 1.964694e-01 3.929389e-01 8.035306e-01 [160,] 1.872942e-01 3.745884e-01 8.127058e-01 [161,] 1.794983e-01 3.589967e-01 8.205017e-01 [162,] 1.618014e-01 3.236027e-01 8.381986e-01 [163,] 1.480625e-01 2.961250e-01 8.519375e-01 [164,] 1.340211e-01 2.680422e-01 8.659789e-01 [165,] 1.193389e-01 2.386778e-01 8.806611e-01 [166,] 1.056876e-01 2.113752e-01 8.943124e-01 [167,] 9.325826e-02 1.865165e-01 9.067417e-01 [168,] 8.437847e-02 1.687569e-01 9.156215e-01 [169,] 7.386283e-02 1.477257e-01 9.261372e-01 [170,] 6.561576e-02 1.312315e-01 9.343842e-01 [171,] 5.896315e-02 1.179263e-01 9.410368e-01 [172,] 5.296953e-02 1.059391e-01 9.470305e-01 [173,] 5.292991e-02 1.058598e-01 9.470701e-01 [174,] 5.159874e-02 1.031975e-01 9.484013e-01 [175,] 6.198250e-02 1.239650e-01 9.380175e-01 [176,] 5.419346e-02 1.083869e-01 9.458065e-01 [177,] 6.350326e-02 1.270065e-01 9.364967e-01 [178,] 5.811810e-02 1.162362e-01 9.418819e-01 [179,] 5.071795e-02 1.014359e-01 9.492820e-01 [180,] 4.545221e-02 9.090442e-02 9.545478e-01 [181,] 4.125984e-02 8.251967e-02 9.587402e-01 [182,] 3.579667e-02 7.159334e-02 9.642033e-01 [183,] 3.161143e-02 6.322285e-02 9.683886e-01 [184,] 2.828141e-02 5.656282e-02 9.717186e-01 [185,] 2.417456e-02 4.834912e-02 9.758254e-01 [186,] 2.115348e-02 4.230696e-02 9.788465e-01 [187,] 2.141254e-02 4.282508e-02 9.785875e-01 [188,] 2.023789e-02 4.047578e-02 9.797621e-01 [189,] 1.973622e-02 3.947245e-02 9.802638e-01 [190,] 1.927808e-02 3.855617e-02 9.807219e-01 [191,] 1.856732e-02 3.713465e-02 9.814327e-01 [192,] 1.649810e-02 3.299620e-02 9.835019e-01 [193,] 1.594242e-02 3.188485e-02 9.840576e-01 [194,] 1.539827e-02 3.079655e-02 9.846017e-01 [195,] 1.433758e-02 2.867517e-02 9.856624e-01 [196,] 1.216814e-02 2.433628e-02 9.878319e-01 [197,] 1.020224e-02 2.040448e-02 9.897978e-01 [198,] 8.581904e-03 1.716381e-02 9.914181e-01 [199,] 7.219514e-03 1.443903e-02 9.927805e-01 [200,] 5.935395e-03 1.187079e-02 9.940646e-01 [201,] 4.857280e-03 9.714561e-03 9.951427e-01 [202,] 3.956890e-03 7.913780e-03 9.960431e-01 [203,] 3.250733e-03 6.501466e-03 9.967493e-01 [204,] 2.623991e-03 5.247981e-03 9.973760e-01 [205,] 2.268646e-03 4.537291e-03 9.977314e-01 [206,] 1.918568e-03 3.837135e-03 9.980814e-01 [207,] 1.541219e-03 3.082439e-03 9.984588e-01 [208,] 1.377354e-03 2.754708e-03 9.986226e-01 [209,] 1.140833e-03 2.281666e-03 9.988592e-01 [210,] 9.062064e-04 1.812413e-03 9.990938e-01 [211,] 7.564474e-04 1.512895e-03 9.992436e-01 [212,] 7.458737e-04 1.491747e-03 9.992541e-01 [213,] 7.141789e-04 1.428358e-03 9.992858e-01 [214,] 6.825969e-04 1.365194e-03 9.993174e-01 [215,] 1.025697e-03 2.051393e-03 9.989743e-01 [216,] 8.371816e-04 1.674363e-03 9.991628e-01 [217,] 7.027261e-04 1.405452e-03 9.992973e-01 [218,] 5.779794e-04 1.155959e-03 9.994220e-01 [219,] 4.838227e-04 9.676454e-04 9.995162e-01 [220,] 4.063374e-04 8.126748e-04 9.995937e-01 [221,] 3.514075e-04 7.028150e-04 9.996486e-01 [222,] 2.940221e-04 5.880442e-04 9.997060e-01 [223,] 2.483281e-04 4.966561e-04 9.997517e-01 [224,] 2.186839e-04 4.373679e-04 9.997813e-01 [225,] 1.940512e-04 3.881024e-04 9.998059e-01 [226,] 1.494455e-04 2.988909e-04 9.998506e-01 [227,] 1.182776e-04 2.365552e-04 9.998817e-01 [228,] 8.969758e-05 1.793952e-04 9.999103e-01 [229,] 6.900170e-05 1.380034e-04 9.999310e-01 [230,] 5.308981e-05 1.061796e-04 9.999469e-01 [231,] 3.972095e-05 7.944190e-05 9.999603e-01 [232,] 3.188302e-05 6.376604e-05 9.999681e-01 [233,] 2.265463e-03 4.530926e-03 9.977345e-01 [234,] 3.421373e-03 6.842746e-03 9.965786e-01 [235,] 5.496936e-02 1.099387e-01 9.450306e-01 [236,] 1.030390e-01 2.060780e-01 8.969610e-01 [237,] 1.061002e-01 2.122004e-01 8.938998e-01 [238,] 1.279744e-01 2.559488e-01 8.720256e-01 [239,] 1.230748e-01 2.461495e-01 8.769252e-01 [240,] 1.372939e-01 2.745877e-01 8.627061e-01 [241,] 1.263267e-01 2.526535e-01 8.736733e-01 [242,] 1.162265e-01 2.324529e-01 8.837735e-01 [243,] 1.063360e-01 2.126720e-01 8.936640e-01 [244,] 9.550874e-02 1.910175e-01 9.044913e-01 [245,] 8.798599e-02 1.759720e-01 9.120140e-01 [246,] 8.243884e-02 1.648777e-01 9.175612e-01 [247,] 7.281352e-02 1.456270e-01 9.271865e-01 [248,] 6.367137e-02 1.273427e-01 9.363286e-01 [249,] 5.531894e-02 1.106379e-01 9.446811e-01 [250,] 4.827372e-02 9.654743e-02 9.517263e-01 [251,] 4.156922e-02 8.313845e-02 9.584308e-01 [252,] 3.593465e-02 7.186930e-02 9.640653e-01 [253,] 3.109063e-02 6.218127e-02 9.689094e-01 [254,] 2.694702e-02 5.389404e-02 9.730530e-01 [255,] 2.293281e-02 4.586561e-02 9.770672e-01 [256,] 1.948837e-02 3.897675e-02 9.805116e-01 [257,] 1.656132e-02 3.312264e-02 9.834387e-01 [258,] 1.382159e-02 2.764318e-02 9.861784e-01 [259,] 1.156986e-02 2.313972e-02 9.884301e-01 [260,] 9.601201e-03 1.920240e-02 9.903988e-01 [261,] 7.962442e-03 1.592488e-02 9.920376e-01 [262,] 7.512445e-03 1.502489e-02 9.924876e-01 [263,] 7.838881e-03 1.567776e-02 9.921611e-01 [264,] 7.685302e-03 1.537060e-02 9.923147e-01 [265,] 6.376080e-03 1.275216e-02 9.936239e-01 [266,] 5.282910e-03 1.056582e-02 9.947171e-01 [267,] 4.299503e-03 8.599006e-03 9.957005e-01 [268,] 3.467569e-03 6.935138e-03 9.965324e-01 [269,] 2.857971e-03 5.715941e-03 9.971420e-01 [270,] 2.287994e-03 4.575988e-03 9.977120e-01 [271,] 3.841158e-03 7.682316e-03 9.961588e-01 [272,] 3.194689e-03 6.389378e-03 9.968053e-01 [273,] 7.782784e-01 4.434433e-01 2.217216e-01 [274,] 7.979951e-01 4.040099e-01 2.020049e-01 [275,] 8.149919e-01 3.700161e-01 1.850081e-01 [276,] 8.314027e-01 3.371947e-01 1.685973e-01 [277,] 8.485124e-01 3.029753e-01 1.514876e-01 [278,] 8.661040e-01 2.677921e-01 1.338960e-01 [279,] 9.127235e-01 1.745530e-01 8.727651e-02 [280,] 9.262354e-01 1.475292e-01 7.376462e-02 [281,] 9.174684e-01 1.650633e-01 8.253165e-02 [282,] 9.091628e-01 1.816744e-01 9.083720e-02 [283,] 9.018736e-01 1.962528e-01 9.812641e-02 [284,] 8.942037e-01 2.115926e-01 1.057963e-01 [285,] 8.871781e-01 2.256437e-01 1.128219e-01 [286,] 8.796592e-01 2.406816e-01 1.203408e-01 [287,] 8.726483e-01 2.547035e-01 1.273517e-01 [288,] 8.669464e-01 2.661072e-01 1.330536e-01 [289,] 9.697658e-01 6.046848e-02 3.023424e-02 [290,] 9.647301e-01 7.053990e-02 3.526995e-02 [291,] 9.588838e-01 8.223243e-02 4.111622e-02 [292,] 9.525581e-01 9.488372e-02 4.744186e-02 [293,] 9.445503e-01 1.108993e-01 5.544966e-02 [294,] 9.355460e-01 1.289080e-01 6.445402e-02 [295,] 9.251497e-01 1.497006e-01 7.485029e-02 [296,] 9.188305e-01 1.623390e-01 8.116951e-02 [297,] 9.069063e-01 1.861873e-01 9.309366e-02 [298,] 8.932598e-01 2.134805e-01 1.067402e-01 [299,] 8.808791e-01 2.382419e-01 1.191209e-01 [300,] 8.757023e-01 2.485954e-01 1.242977e-01 [301,] 8.776129e-01 2.447742e-01 1.223871e-01 [302,] 8.834752e-01 2.330497e-01 1.165248e-01 [303,] 8.977678e-01 2.044644e-01 1.022322e-01 [304,] 8.832631e-01 2.334738e-01 1.167369e-01 [305,] 8.968036e-01 2.063928e-01 1.031964e-01 [306,] 8.829158e-01 2.341685e-01 1.170842e-01 [307,] 8.690630e-01 2.618741e-01 1.309370e-01 [308,] 8.523468e-01 2.953065e-01 1.476532e-01 [309,] 8.352443e-01 3.295114e-01 1.647557e-01 [310,] 8.145924e-01 3.708152e-01 1.854076e-01 [311,] 8.914646e-01 2.170709e-01 1.085354e-01 [312,] 8.908391e-01 2.183217e-01 1.091609e-01 [313,] 8.981984e-01 2.036032e-01 1.018016e-01 [314,] 8.860156e-01 2.279689e-01 1.139844e-01 [315,] 8.705655e-01 2.588690e-01 1.294345e-01 [316,] 8.645661e-01 2.708679e-01 1.354339e-01 [317,] 8.562870e-01 2.874260e-01 1.437130e-01 [318,] 8.654160e-01 2.691680e-01 1.345840e-01 [319,] 8.669992e-01 2.660016e-01 1.330008e-01 [320,] 8.491583e-01 3.016834e-01 1.508417e-01 [321,] 8.296825e-01 3.406351e-01 1.703175e-01 [322,] 8.076739e-01 3.846522e-01 1.923261e-01 [323,] 7.839812e-01 4.320376e-01 2.160188e-01 [324,] 7.582048e-01 4.835905e-01 2.417952e-01 [325,] 7.922081e-01 4.155839e-01 2.077919e-01 [326,] 7.681773e-01 4.636455e-01 2.318227e-01 [327,] 7.466176e-01 5.067647e-01 2.533824e-01 [328,] 7.184060e-01 5.631880e-01 2.815940e-01 [329,] 7.159737e-01 5.680525e-01 2.840263e-01 [330,] 6.971270e-01 6.057460e-01 3.028730e-01 [331,] 6.913423e-01 6.173154e-01 3.086577e-01 [332,] 6.820579e-01 6.358843e-01 3.179421e-01 [333,] 6.658530e-01 6.682940e-01 3.341470e-01 [334,] 6.511514e-01 6.976973e-01 3.488486e-01 [335,] 6.827562e-01 6.344877e-01 3.172438e-01 [336,] 6.650359e-01 6.699283e-01 3.349641e-01 [337,] 6.434724e-01 7.130552e-01 3.565276e-01 [338,] 6.473736e-01 7.052529e-01 3.526264e-01 [339,] 6.413151e-01 7.173697e-01 3.586849e-01 [340,] 6.303541e-01 7.392918e-01 3.696459e-01 [341,] 6.354571e-01 7.290858e-01 3.645429e-01 [342,] 6.019631e-01 7.960737e-01 3.980369e-01 [343,] 6.023992e-01 7.952016e-01 3.976008e-01 [344,] 6.300597e-01 7.398806e-01 3.699403e-01 [345,] 7.423377e-01 5.153246e-01 2.576623e-01 [346,] 1.000000e+00 1.394802e-10 6.974011e-11 [347,] 1.000000e+00 4.648572e-13 2.324286e-13 [348,] 1.000000e+00 9.070194e-19 4.535097e-19 [349,] 1.000000e+00 2.797657e-18 1.398828e-18 [350,] 1.000000e+00 8.563706e-18 4.281853e-18 [351,] 1.000000e+00 2.557067e-17 1.278534e-17 [352,] 1.000000e+00 7.749200e-17 3.874600e-17 [353,] 1.000000e+00 1.977538e-16 9.887689e-17 [354,] 1.000000e+00 5.909194e-16 2.954597e-16 [355,] 1.000000e+00 1.745141e-15 8.725705e-16 [356,] 1.000000e+00 4.991345e-15 2.495673e-15 [357,] 1.000000e+00 1.411695e-14 7.058474e-15 [358,] 1.000000e+00 4.038559e-14 2.019280e-14 [359,] 1.000000e+00 6.818108e-14 3.409054e-14 [360,] 1.000000e+00 1.635830e-13 8.179151e-14 [361,] 1.000000e+00 4.456747e-13 2.228374e-13 [362,] 1.000000e+00 1.130368e-12 5.651841e-13 [363,] 1.000000e+00 3.094730e-12 1.547365e-12 [364,] 1.000000e+00 6.241849e-12 3.120925e-12 [365,] 1.000000e+00 7.485619e-12 3.742809e-12 [366,] 1.000000e+00 1.974381e-11 9.871907e-12 [367,] 1.000000e+00 5.152911e-11 2.576456e-11 [368,] 1.000000e+00 1.026822e-10 5.134108e-11 [369,] 1.000000e+00 1.886690e-10 9.433451e-11 [370,] 1.000000e+00 3.999594e-10 1.999797e-10 [371,] 1.000000e+00 4.465540e-10 2.232770e-10 [372,] 1.000000e+00 1.119007e-09 5.595035e-10 [373,] 1.000000e+00 2.240344e-09 1.120172e-09 [374,] 1.000000e+00 3.789802e-09 1.894901e-09 [375,] 1.000000e+00 6.655373e-09 3.327687e-09 [376,] 1.000000e+00 1.843343e-09 9.216715e-10 [377,] 1.000000e+00 1.469799e-13 7.348994e-14 [378,] 1.000000e+00 3.115920e-18 1.557960e-18 [379,] 1.000000e+00 9.632580e-22 4.816290e-22 [380,] 1.000000e+00 8.523750e-21 4.261875e-21 [381,] 1.000000e+00 6.615324e-20 3.307662e-20 [382,] 1.000000e+00 5.435512e-19 2.717756e-19 [383,] 1.000000e+00 4.605667e-18 2.302833e-18 [384,] 1.000000e+00 1.144106e-18 5.720528e-19 [385,] 1.000000e+00 2.347027e-21 1.173513e-21 [386,] 1.000000e+00 5.563664e-28 2.781832e-28 [387,] 1.000000e+00 1.721774e-26 8.608871e-27 [388,] 1.000000e+00 5.240890e-25 2.620445e-25 [389,] 1.000000e+00 1.418878e-23 7.094389e-24 [390,] 1.000000e+00 3.906285e-22 1.953142e-22 [391,] 1.000000e+00 9.814741e-21 4.907371e-21 [392,] 1.000000e+00 2.691137e-19 1.345568e-19 [393,] 1.000000e+00 7.054742e-18 3.527371e-18 [394,] 1.000000e+00 1.414344e-16 7.071718e-17 [395,] 1.000000e+00 3.373651e-15 1.686825e-15 [396,] 1.000000e+00 7.946003e-14 3.973002e-14 [397,] 1.000000e+00 1.560026e-12 7.800130e-13 [398,] 1.000000e+00 3.402505e-11 1.701253e-11 [399,] 1.000000e+00 3.696942e-10 1.848471e-10 [400,] 1.000000e+00 6.993274e-09 3.496637e-09 [401,] 9.999999e-01 1.425640e-07 7.128199e-08 [402,] 9.999992e-01 1.580879e-06 7.904396e-07 [403,] 9.999828e-01 3.431282e-05 1.715641e-05 [404,] 9.997198e-01 5.603863e-04 2.801932e-04 > postscript(file="/var/www/html/freestat/rcomp/tmp/1rhou1291403799.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/freestat/rcomp/tmp/2k86x1291403799.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/freestat/rcomp/tmp/3k86x1291403799.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/freestat/rcomp/tmp/4k86x1291403799.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/freestat/rcomp/tmp/5k86x1291403799.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') > qqline(mysum$resid) > grid() > dev.off() null device 1 > (myerror <- as.ts(mysum$resid)) Time Series: Start = 1 End = 431 Frequency = 1 1 2 3 4 5 6 -51984.0821 -18372.0453 -51384.7569 -38448.7631 -28204.3908 -55151.6551 7 8 9 10 11 12 -55709.7680 -52533.0890 -47951.0220 -54073.2014 -29140.1596 -25285.6329 13 14 15 16 17 18 -53428.4879 -77322.8544 -79133.9065 -80405.2747 -79141.9387 -79282.3001 19 20 21 22 23 24 -78621.9433 -78459.2013 -80174.0945 -30726.2388 95041.9604 32759.1670 25 26 27 28 29 30 29915.8614 17150.3681 91785.8572 33861.8705 -42141.6666 -102645.2920 31 32 33 34 35 36 -105839.0484 -35023.5739 16584.5218 -24762.1198 -39826.6969 -26299.3480 37 38 39 40 41 42 -48066.7863 203220.1216 -11222.0719 -7504.4171 -26067.1170 -17362.2156 43 44 45 46 47 48 -38607.7672 -4853.7732 1296.8835 -195.4289 -45095.7820 -29801.9203 49 50 51 52 53 54 2311.3611 780.5657 -40474.6696 -42507.0539 -5850.3256 115807.4229 55 56 57 58 59 60 35862.4758 -42698.6495 -52010.5455 -45048.6944 -52510.2117 264105.7978 61 62 63 64 65 66 118534.8649 203754.6431 15169.0314 202766.3805 260528.1083 218013.3857 67 68 69 70 71 72 212545.9028 203053.0516 11875.1888 -76000.6433 224517.6363 -77622.5047 73 74 75 76 77 78 -46561.4052 -52333.5476 -50190.7121 -52755.5158 -72851.0447 130344.3494 79 80 81 82 83 84 -49302.1727 -53421.9134 -80181.6776 -15917.2363 34210.1823 -43847.9955 85 86 87 88 89 90 -103111.7453 -94791.3445 -104808.1266 -103379.0923 -43084.6284 146690.6464 91 92 93 94 95 96 -73791.4766 -35904.4095 -33232.7310 -38570.9399 729.7612 2544.9957 97 98 99 100 101 102 -13441.4433 47369.3150 5654.4717 2266.4497 12495.0906 17691.9284 103 104 105 106 107 108 -104570.5179 -8747.1120 158518.1673 -56401.4898 -48795.5967 -45863.7494 109 110 111 112 113 114 -63797.5544 -52447.5941 -99711.4355 207730.2291 302876.7551 746714.1728 115 116 117 118 119 120 586239.2520 419117.0070 126309.0107 203263.1535 -30759.6203 -60815.3023 121 122 123 124 125 126 -66764.7204 -77484.1720 -48535.5660 -53364.5011 -60449.4067 -76767.6583 127 128 129 130 131 132 -76620.2602 -80177.5438 34254.9020 -43846.7676 -98713.5407 -76824.0443 133 134 135 136 137 138 -51274.0349 -53391.9593 -79977.1345 24431.4312 -57312.2709 34214.6736 139 140 141 142 143 144 -43783.5052 -104469.8502 -37030.7118 -27701.4535 -42646.5360 -38526.2130 145 146 147 148 149 150 227.5167 -51326.5332 -52433.6198 -100544.0350 -98387.3333 -98393.4294 151 152 153 154 155 156 -101706.6340 212378.0938 368592.9577 202838.8708 203257.1363 -68244.6674 157 158 159 160 161 162 -60957.7973 -76997.0518 -3073.2985 -53368.3318 -79766.9541 -77029.7275 163 164 165 166 167 168 -78834.9155 -62439.7263 -58091.1291 -85013.8242 -81764.6926 -79026.4824 169 170 171 172 173 174 -72533.5203 -80192.9086 -25461.7417 -29274.6169 96401.2540 85379.0730 175 176 177 178 179 180 21518.2627 -43974.4263 -38486.4860 6891.7816 5181.5796 9146.0670 181 182 183 184 185 186 -38435.9034 -394.3779 -37486.4399 -50611.7184 -52495.9551 -99483.0970 187 188 189 190 191 192 -99600.9712 207276.0284 46720.1004 203260.7656 -77415.9932 8207.0563 193 194 195 196 197 198 -53388.8923 -80187.2961 -28600.0927 62536.5457 61927.4653 24464.3847 199 200 201 202 203 204 -41906.6927 -105490.8002 -92169.2391 -104082.7742 -106088.0254 -102462.1405 205 206 207 208 209 210 -91421.0283 -103858.7786 -106424.4940 -95004.2818 -34151.0419 -14680.6530 211 212 213 214 215 216 -36481.1758 -38461.4673 -6314.7240 -5457.9028 -4641.2721 24432.1750 217 218 219 220 221 222 6700.0320 56301.5034 -49972.0866 -13605.4337 65617.6260 -43257.0152 223 224 225 226 227 228 -2975.6361 -52489.8522 -100628.8139 -98786.4913 -102671.3692 213550.8750 229 230 231 232 233 234 -47431.5150 -56480.3045 -61411.8271 -77444.9130 -53270.5395 -63328.5218 235 236 237 238 239 240 -49819.4800 -53396.5585 -78906.1370 -80189.3932 -21281.1783 52445.0103 241 242 243 244 245 246 12826.6985 -43889.7870 -38459.5619 -387.3924 -52504.2159 539836.2739 247 248 249 250 251 252 235437.0846 541138.0297 329130.6501 167075.3061 228098.5024 129210.4649 253 254 255 256 257 258 203260.3915 -72864.4329 -77457.9396 -47237.6400 -45447.7206 -53402.3328 259 260 261 262 263 264 -80106.7779 -28027.5395 -21641.5347 -1540.8844 -29949.9993 -16941.1908 265 266 267 268 269 270 -28783.6810 -34602.8242 -40045.0405 -25598.5392 -28083.8066 37485.2802 271 272 273 274 275 276 10503.5281 26393.0639 19611.9413 -42834.3117 -101985.9126 -106324.5296 277 278 279 280 281 282 -107202.7353 -37211.0991 -38387.5368 -10749.8817 2383.6092 27723.8676 283 284 285 286 287 288 17595.7284 82149.9237 -52509.7685 886488.6528 206698.5841 199350.5118 289 290 291 292 293 294 200772.2861 200366.9096 202837.6765 269324.7964 203208.8164 -77485.0265 295 296 297 298 299 300 -50497.0133 -53407.7054 -81286.3200 -78467.8443 -72312.1892 -78898.0977 301 302 303 304 305 306 -88856.0022 239288.1022 -81018.1094 -80152.1766 -41439.6638 -22694.4401 307 308 309 310 311 312 -22786.5675 9423.3411 88440.5338 22176.0282 -38094.4522 -9827.2335 313 314 315 316 317 318 -77250.4057 -103329.5351 -104043.0651 -106285.8564 -38756.6620 90705.3731 319 320 321 322 323 324 -962.5984 1176.6623 -21421.6578 -31058.0088 -52517.5116 274904.7996 325 326 327 328 329 330 176858.6770 203288.5967 -57904.2013 -77497.7190 -51979.2262 -49075.2270 331 332 333 334 335 336 -89015.1393 -53407.8385 -80309.3676 -80199.4864 34328.5342 18604.2797 337 338 339 340 341 342 -43002.0530 -105813.5158 -38545.4547 -397.6282 -52119.1919 -87911.3164 343 344 345 346 347 348 -82550.4707 -105347.7632 -101848.0355 -95888.2868 -100326.9524 -106648.7371 349 350 351 352 353 354 -101974.8185 -99408.2689 114216.8324 202278.1743 186032.8914 213097.5189 355 356 357 358 359 360 96009.3980 202711.7632 230794.4162 314851.5059 820222.6083 195274.4309 361 362 363 364 365 366 203133.1586 -67260.8431 -59724.8910 -76213.0596 -67253.9069 -38779.2886 367 368 369 370 371 372 -68362.3153 -68942.4717 -77473.9239 -53424.4362 -80196.2496 34221.6330 373 374 375 376 377 378 -37094.1100 -106479.6225 -103605.4195 -38571.8488 -396.3090 -52487.6426 379 380 381 382 383 384 -99430.7005 -92273.4987 -91663.3273 -96975.6662 -101066.8817 -90234.5815 385 386 387 388 389 390 -100054.6555 -98914.5719 -101127.2940 -96103.3282 -93061.5587 -88555.6314 391 392 393 394 395 396 -99597.0803 213546.4967 -70653.6435 -69434.5041 -66995.4019 -34945.9212 397 398 399 400 401 402 -76200.7920 -19756.0272 -54133.6556 -52052.4458 -55322.1652 217267.5293 403 404 405 406 407 408 -53945.3063 -54767.1197 -97338.8189 -53382.0542 -74824.7040 -72596.7299 409 410 411 412 413 414 -80111.6225 -26677.6112 -24847.6612 22478.3077 -23897.8391 -6970.4024 415 416 417 418 419 420 -31477.8782 -32846.2212 -31924.1019 -24540.5460 -18844.6989 -28065.5915 421 422 423 424 425 426 -4028.1168 3478.4659 35208.0174 17148.7214 12362.0026 7260.0986 427 428 429 430 431 12969.3097 13416.6999 14536.8355 13891.0421 -44675.4567 > postscript(file="/var/www/html/freestat/rcomp/tmp/6chni1291403799.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 -51984.0821 NA 1 -18372.0453 -51984.0821 2 -51384.7569 -18372.0453 3 -38448.7631 -51384.7569 4 -28204.3908 -38448.7631 5 -55151.6551 -28204.3908 6 -55709.7680 -55151.6551 7 -52533.0890 -55709.7680 8 -47951.0220 -52533.0890 9 -54073.2014 -47951.0220 10 -29140.1596 -54073.2014 11 -25285.6329 -29140.1596 12 -53428.4879 -25285.6329 13 -77322.8544 -53428.4879 14 -79133.9065 -77322.8544 15 -80405.2747 -79133.9065 16 -79141.9387 -80405.2747 17 -79282.3001 -79141.9387 18 -78621.9433 -79282.3001 19 -78459.2013 -78621.9433 20 -80174.0945 -78459.2013 21 -30726.2388 -80174.0945 22 95041.9604 -30726.2388 23 32759.1670 95041.9604 24 29915.8614 32759.1670 25 17150.3681 29915.8614 26 91785.8572 17150.3681 27 33861.8705 91785.8572 28 -42141.6666 33861.8705 29 -102645.2920 -42141.6666 30 -105839.0484 -102645.2920 31 -35023.5739 -105839.0484 32 16584.5218 -35023.5739 33 -24762.1198 16584.5218 34 -39826.6969 -24762.1198 35 -26299.3480 -39826.6969 36 -48066.7863 -26299.3480 37 203220.1216 -48066.7863 38 -11222.0719 203220.1216 39 -7504.4171 -11222.0719 40 -26067.1170 -7504.4171 41 -17362.2156 -26067.1170 42 -38607.7672 -17362.2156 43 -4853.7732 -38607.7672 44 1296.8835 -4853.7732 45 -195.4289 1296.8835 46 -45095.7820 -195.4289 47 -29801.9203 -45095.7820 48 2311.3611 -29801.9203 49 780.5657 2311.3611 50 -40474.6696 780.5657 51 -42507.0539 -40474.6696 52 -5850.3256 -42507.0539 53 115807.4229 -5850.3256 54 35862.4758 115807.4229 55 -42698.6495 35862.4758 56 -52010.5455 -42698.6495 57 -45048.6944 -52010.5455 58 -52510.2117 -45048.6944 59 264105.7978 -52510.2117 60 118534.8649 264105.7978 61 203754.6431 118534.8649 62 15169.0314 203754.6431 63 202766.3805 15169.0314 64 260528.1083 202766.3805 65 218013.3857 260528.1083 66 212545.9028 218013.3857 67 203053.0516 212545.9028 68 11875.1888 203053.0516 69 -76000.6433 11875.1888 70 224517.6363 -76000.6433 71 -77622.5047 224517.6363 72 -46561.4052 -77622.5047 73 -52333.5476 -46561.4052 74 -50190.7121 -52333.5476 75 -52755.5158 -50190.7121 76 -72851.0447 -52755.5158 77 130344.3494 -72851.0447 78 -49302.1727 130344.3494 79 -53421.9134 -49302.1727 80 -80181.6776 -53421.9134 81 -15917.2363 -80181.6776 82 34210.1823 -15917.2363 83 -43847.9955 34210.1823 84 -103111.7453 -43847.9955 85 -94791.3445 -103111.7453 86 -104808.1266 -94791.3445 87 -103379.0923 -104808.1266 88 -43084.6284 -103379.0923 89 146690.6464 -43084.6284 90 -73791.4766 146690.6464 91 -35904.4095 -73791.4766 92 -33232.7310 -35904.4095 93 -38570.9399 -33232.7310 94 729.7612 -38570.9399 95 2544.9957 729.7612 96 -13441.4433 2544.9957 97 47369.3150 -13441.4433 98 5654.4717 47369.3150 99 2266.4497 5654.4717 100 12495.0906 2266.4497 101 17691.9284 12495.0906 102 -104570.5179 17691.9284 103 -8747.1120 -104570.5179 104 158518.1673 -8747.1120 105 -56401.4898 158518.1673 106 -48795.5967 -56401.4898 107 -45863.7494 -48795.5967 108 -63797.5544 -45863.7494 109 -52447.5941 -63797.5544 110 -99711.4355 -52447.5941 111 207730.2291 -99711.4355 112 302876.7551 207730.2291 113 746714.1728 302876.7551 114 586239.2520 746714.1728 115 419117.0070 586239.2520 116 126309.0107 419117.0070 117 203263.1535 126309.0107 118 -30759.6203 203263.1535 119 -60815.3023 -30759.6203 120 -66764.7204 -60815.3023 121 -77484.1720 -66764.7204 122 -48535.5660 -77484.1720 123 -53364.5011 -48535.5660 124 -60449.4067 -53364.5011 125 -76767.6583 -60449.4067 126 -76620.2602 -76767.6583 127 -80177.5438 -76620.2602 128 34254.9020 -80177.5438 129 -43846.7676 34254.9020 130 -98713.5407 -43846.7676 131 -76824.0443 -98713.5407 132 -51274.0349 -76824.0443 133 -53391.9593 -51274.0349 134 -79977.1345 -53391.9593 135 24431.4312 -79977.1345 136 -57312.2709 24431.4312 137 34214.6736 -57312.2709 138 -43783.5052 34214.6736 139 -104469.8502 -43783.5052 140 -37030.7118 -104469.8502 141 -27701.4535 -37030.7118 142 -42646.5360 -27701.4535 143 -38526.2130 -42646.5360 144 227.5167 -38526.2130 145 -51326.5332 227.5167 146 -52433.6198 -51326.5332 147 -100544.0350 -52433.6198 148 -98387.3333 -100544.0350 149 -98393.4294 -98387.3333 150 -101706.6340 -98393.4294 151 212378.0938 -101706.6340 152 368592.9577 212378.0938 153 202838.8708 368592.9577 154 203257.1363 202838.8708 155 -68244.6674 203257.1363 156 -60957.7973 -68244.6674 157 -76997.0518 -60957.7973 158 -3073.2985 -76997.0518 159 -53368.3318 -3073.2985 160 -79766.9541 -53368.3318 161 -77029.7275 -79766.9541 162 -78834.9155 -77029.7275 163 -62439.7263 -78834.9155 164 -58091.1291 -62439.7263 165 -85013.8242 -58091.1291 166 -81764.6926 -85013.8242 167 -79026.4824 -81764.6926 168 -72533.5203 -79026.4824 169 -80192.9086 -72533.5203 170 -25461.7417 -80192.9086 171 -29274.6169 -25461.7417 172 96401.2540 -29274.6169 173 85379.0730 96401.2540 174 21518.2627 85379.0730 175 -43974.4263 21518.2627 176 -38486.4860 -43974.4263 177 6891.7816 -38486.4860 178 5181.5796 6891.7816 179 9146.0670 5181.5796 180 -38435.9034 9146.0670 181 -394.3779 -38435.9034 182 -37486.4399 -394.3779 183 -50611.7184 -37486.4399 184 -52495.9551 -50611.7184 185 -99483.0970 -52495.9551 186 -99600.9712 -99483.0970 187 207276.0284 -99600.9712 188 46720.1004 207276.0284 189 203260.7656 46720.1004 190 -77415.9932 203260.7656 191 8207.0563 -77415.9932 192 -53388.8923 8207.0563 193 -80187.2961 -53388.8923 194 -28600.0927 -80187.2961 195 62536.5457 -28600.0927 196 61927.4653 62536.5457 197 24464.3847 61927.4653 198 -41906.6927 24464.3847 199 -105490.8002 -41906.6927 200 -92169.2391 -105490.8002 201 -104082.7742 -92169.2391 202 -106088.0254 -104082.7742 203 -102462.1405 -106088.0254 204 -91421.0283 -102462.1405 205 -103858.7786 -91421.0283 206 -106424.4940 -103858.7786 207 -95004.2818 -106424.4940 208 -34151.0419 -95004.2818 209 -14680.6530 -34151.0419 210 -36481.1758 -14680.6530 211 -38461.4673 -36481.1758 212 -6314.7240 -38461.4673 213 -5457.9028 -6314.7240 214 -4641.2721 -5457.9028 215 24432.1750 -4641.2721 216 6700.0320 24432.1750 217 56301.5034 6700.0320 218 -49972.0866 56301.5034 219 -13605.4337 -49972.0866 220 65617.6260 -13605.4337 221 -43257.0152 65617.6260 222 -2975.6361 -43257.0152 223 -52489.8522 -2975.6361 224 -100628.8139 -52489.8522 225 -98786.4913 -100628.8139 226 -102671.3692 -98786.4913 227 213550.8750 -102671.3692 228 -47431.5150 213550.8750 229 -56480.3045 -47431.5150 230 -61411.8271 -56480.3045 231 -77444.9130 -61411.8271 232 -53270.5395 -77444.9130 233 -63328.5218 -53270.5395 234 -49819.4800 -63328.5218 235 -53396.5585 -49819.4800 236 -78906.1370 -53396.5585 237 -80189.3932 -78906.1370 238 -21281.1783 -80189.3932 239 52445.0103 -21281.1783 240 12826.6985 52445.0103 241 -43889.7870 12826.6985 242 -38459.5619 -43889.7870 243 -387.3924 -38459.5619 244 -52504.2159 -387.3924 245 539836.2739 -52504.2159 246 235437.0846 539836.2739 247 541138.0297 235437.0846 248 329130.6501 541138.0297 249 167075.3061 329130.6501 250 228098.5024 167075.3061 251 129210.4649 228098.5024 252 203260.3915 129210.4649 253 -72864.4329 203260.3915 254 -77457.9396 -72864.4329 255 -47237.6400 -77457.9396 256 -45447.7206 -47237.6400 257 -53402.3328 -45447.7206 258 -80106.7779 -53402.3328 259 -28027.5395 -80106.7779 260 -21641.5347 -28027.5395 261 -1540.8844 -21641.5347 262 -29949.9993 -1540.8844 263 -16941.1908 -29949.9993 264 -28783.6810 -16941.1908 265 -34602.8242 -28783.6810 266 -40045.0405 -34602.8242 267 -25598.5392 -40045.0405 268 -28083.8066 -25598.5392 269 37485.2802 -28083.8066 270 10503.5281 37485.2802 271 26393.0639 10503.5281 272 19611.9413 26393.0639 273 -42834.3117 19611.9413 274 -101985.9126 -42834.3117 275 -106324.5296 -101985.9126 276 -107202.7353 -106324.5296 277 -37211.0991 -107202.7353 278 -38387.5368 -37211.0991 279 -10749.8817 -38387.5368 280 2383.6092 -10749.8817 281 27723.8676 2383.6092 282 17595.7284 27723.8676 283 82149.9237 17595.7284 284 -52509.7685 82149.9237 285 886488.6528 -52509.7685 286 206698.5841 886488.6528 287 199350.5118 206698.5841 288 200772.2861 199350.5118 289 200366.9096 200772.2861 290 202837.6765 200366.9096 291 269324.7964 202837.6765 292 203208.8164 269324.7964 293 -77485.0265 203208.8164 294 -50497.0133 -77485.0265 295 -53407.7054 -50497.0133 296 -81286.3200 -53407.7054 297 -78467.8443 -81286.3200 298 -72312.1892 -78467.8443 299 -78898.0977 -72312.1892 300 -88856.0022 -78898.0977 301 239288.1022 -88856.0022 302 -81018.1094 239288.1022 303 -80152.1766 -81018.1094 304 -41439.6638 -80152.1766 305 -22694.4401 -41439.6638 306 -22786.5675 -22694.4401 307 9423.3411 -22786.5675 308 88440.5338 9423.3411 309 22176.0282 88440.5338 310 -38094.4522 22176.0282 311 -9827.2335 -38094.4522 312 -77250.4057 -9827.2335 313 -103329.5351 -77250.4057 314 -104043.0651 -103329.5351 315 -106285.8564 -104043.0651 316 -38756.6620 -106285.8564 317 90705.3731 -38756.6620 318 -962.5984 90705.3731 319 1176.6623 -962.5984 320 -21421.6578 1176.6623 321 -31058.0088 -21421.6578 322 -52517.5116 -31058.0088 323 274904.7996 -52517.5116 324 176858.6770 274904.7996 325 203288.5967 176858.6770 326 -57904.2013 203288.5967 327 -77497.7190 -57904.2013 328 -51979.2262 -77497.7190 329 -49075.2270 -51979.2262 330 -89015.1393 -49075.2270 331 -53407.8385 -89015.1393 332 -80309.3676 -53407.8385 333 -80199.4864 -80309.3676 334 34328.5342 -80199.4864 335 18604.2797 34328.5342 336 -43002.0530 18604.2797 337 -105813.5158 -43002.0530 338 -38545.4547 -105813.5158 339 -397.6282 -38545.4547 340 -52119.1919 -397.6282 341 -87911.3164 -52119.1919 342 -82550.4707 -87911.3164 343 -105347.7632 -82550.4707 344 -101848.0355 -105347.7632 345 -95888.2868 -101848.0355 346 -100326.9524 -95888.2868 347 -106648.7371 -100326.9524 348 -101974.8185 -106648.7371 349 -99408.2689 -101974.8185 350 114216.8324 -99408.2689 351 202278.1743 114216.8324 352 186032.8914 202278.1743 353 213097.5189 186032.8914 354 96009.3980 213097.5189 355 202711.7632 96009.3980 356 230794.4162 202711.7632 357 314851.5059 230794.4162 358 820222.6083 314851.5059 359 195274.4309 820222.6083 360 203133.1586 195274.4309 361 -67260.8431 203133.1586 362 -59724.8910 -67260.8431 363 -76213.0596 -59724.8910 364 -67253.9069 -76213.0596 365 -38779.2886 -67253.9069 366 -68362.3153 -38779.2886 367 -68942.4717 -68362.3153 368 -77473.9239 -68942.4717 369 -53424.4362 -77473.9239 370 -80196.2496 -53424.4362 371 34221.6330 -80196.2496 372 -37094.1100 34221.6330 373 -106479.6225 -37094.1100 374 -103605.4195 -106479.6225 375 -38571.8488 -103605.4195 376 -396.3090 -38571.8488 377 -52487.6426 -396.3090 378 -99430.7005 -52487.6426 379 -92273.4987 -99430.7005 380 -91663.3273 -92273.4987 381 -96975.6662 -91663.3273 382 -101066.8817 -96975.6662 383 -90234.5815 -101066.8817 384 -100054.6555 -90234.5815 385 -98914.5719 -100054.6555 386 -101127.2940 -98914.5719 387 -96103.3282 -101127.2940 388 -93061.5587 -96103.3282 389 -88555.6314 -93061.5587 390 -99597.0803 -88555.6314 391 213546.4967 -99597.0803 392 -70653.6435 213546.4967 393 -69434.5041 -70653.6435 394 -66995.4019 -69434.5041 395 -34945.9212 -66995.4019 396 -76200.7920 -34945.9212 397 -19756.0272 -76200.7920 398 -54133.6556 -19756.0272 399 -52052.4458 -54133.6556 400 -55322.1652 -52052.4458 401 217267.5293 -55322.1652 402 -53945.3063 217267.5293 403 -54767.1197 -53945.3063 404 -97338.8189 -54767.1197 405 -53382.0542 -97338.8189 406 -74824.7040 -53382.0542 407 -72596.7299 -74824.7040 408 -80111.6225 -72596.7299 409 -26677.6112 -80111.6225 410 -24847.6612 -26677.6112 411 22478.3077 -24847.6612 412 -23897.8391 22478.3077 413 -6970.4024 -23897.8391 414 -31477.8782 -6970.4024 415 -32846.2212 -31477.8782 416 -31924.1019 -32846.2212 417 -24540.5460 -31924.1019 418 -18844.6989 -24540.5460 419 -28065.5915 -18844.6989 420 -4028.1168 -28065.5915 421 3478.4659 -4028.1168 422 35208.0174 3478.4659 423 17148.7214 35208.0174 424 12362.0026 17148.7214 425 7260.0986 12362.0026 426 12969.3097 7260.0986 427 13416.6999 12969.3097 428 14536.8355 13416.6999 429 13891.0421 14536.8355 430 -44675.4567 13891.0421 431 NA -44675.4567 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -18372.0453 -51984.0821 [2,] -51384.7569 -18372.0453 [3,] -38448.7631 -51384.7569 [4,] -28204.3908 -38448.7631 [5,] -55151.6551 -28204.3908 [6,] -55709.7680 -55151.6551 [7,] -52533.0890 -55709.7680 [8,] -47951.0220 -52533.0890 [9,] -54073.2014 -47951.0220 [10,] -29140.1596 -54073.2014 [11,] -25285.6329 -29140.1596 [12,] -53428.4879 -25285.6329 [13,] -77322.8544 -53428.4879 [14,] -79133.9065 -77322.8544 [15,] -80405.2747 -79133.9065 [16,] -79141.9387 -80405.2747 [17,] -79282.3001 -79141.9387 [18,] -78621.9433 -79282.3001 [19,] -78459.2013 -78621.9433 [20,] -80174.0945 -78459.2013 [21,] -30726.2388 -80174.0945 [22,] 95041.9604 -30726.2388 [23,] 32759.1670 95041.9604 [24,] 29915.8614 32759.1670 [25,] 17150.3681 29915.8614 [26,] 91785.8572 17150.3681 [27,] 33861.8705 91785.8572 [28,] -42141.6666 33861.8705 [29,] -102645.2920 -42141.6666 [30,] -105839.0484 -102645.2920 [31,] -35023.5739 -105839.0484 [32,] 16584.5218 -35023.5739 [33,] -24762.1198 16584.5218 [34,] -39826.6969 -24762.1198 [35,] -26299.3480 -39826.6969 [36,] -48066.7863 -26299.3480 [37,] 203220.1216 -48066.7863 [38,] -11222.0719 203220.1216 [39,] -7504.4171 -11222.0719 [40,] -26067.1170 -7504.4171 [41,] -17362.2156 -26067.1170 [42,] -38607.7672 -17362.2156 [43,] -4853.7732 -38607.7672 [44,] 1296.8835 -4853.7732 [45,] -195.4289 1296.8835 [46,] -45095.7820 -195.4289 [47,] -29801.9203 -45095.7820 [48,] 2311.3611 -29801.9203 [49,] 780.5657 2311.3611 [50,] -40474.6696 780.5657 [51,] -42507.0539 -40474.6696 [52,] -5850.3256 -42507.0539 [53,] 115807.4229 -5850.3256 [54,] 35862.4758 115807.4229 [55,] -42698.6495 35862.4758 [56,] -52010.5455 -42698.6495 [57,] -45048.6944 -52010.5455 [58,] -52510.2117 -45048.6944 [59,] 264105.7978 -52510.2117 [60,] 118534.8649 264105.7978 [61,] 203754.6431 118534.8649 [62,] 15169.0314 203754.6431 [63,] 202766.3805 15169.0314 [64,] 260528.1083 202766.3805 [65,] 218013.3857 260528.1083 [66,] 212545.9028 218013.3857 [67,] 203053.0516 212545.9028 [68,] 11875.1888 203053.0516 [69,] -76000.6433 11875.1888 [70,] 224517.6363 -76000.6433 [71,] -77622.5047 224517.6363 [72,] -46561.4052 -77622.5047 [73,] -52333.5476 -46561.4052 [74,] -50190.7121 -52333.5476 [75,] -52755.5158 -50190.7121 [76,] -72851.0447 -52755.5158 [77,] 130344.3494 -72851.0447 [78,] -49302.1727 130344.3494 [79,] -53421.9134 -49302.1727 [80,] -80181.6776 -53421.9134 [81,] -15917.2363 -80181.6776 [82,] 34210.1823 -15917.2363 [83,] -43847.9955 34210.1823 [84,] -103111.7453 -43847.9955 [85,] -94791.3445 -103111.7453 [86,] -104808.1266 -94791.3445 [87,] -103379.0923 -104808.1266 [88,] -43084.6284 -103379.0923 [89,] 146690.6464 -43084.6284 [90,] -73791.4766 146690.6464 [91,] -35904.4095 -73791.4766 [92,] -33232.7310 -35904.4095 [93,] -38570.9399 -33232.7310 [94,] 729.7612 -38570.9399 [95,] 2544.9957 729.7612 [96,] -13441.4433 2544.9957 [97,] 47369.3150 -13441.4433 [98,] 5654.4717 47369.3150 [99,] 2266.4497 5654.4717 [100,] 12495.0906 2266.4497 [101,] 17691.9284 12495.0906 [102,] -104570.5179 17691.9284 [103,] -8747.1120 -104570.5179 [104,] 158518.1673 -8747.1120 [105,] -56401.4898 158518.1673 [106,] -48795.5967 -56401.4898 [107,] -45863.7494 -48795.5967 [108,] -63797.5544 -45863.7494 [109,] -52447.5941 -63797.5544 [110,] -99711.4355 -52447.5941 [111,] 207730.2291 -99711.4355 [112,] 302876.7551 207730.2291 [113,] 746714.1728 302876.7551 [114,] 586239.2520 746714.1728 [115,] 419117.0070 586239.2520 [116,] 126309.0107 419117.0070 [117,] 203263.1535 126309.0107 [118,] -30759.6203 203263.1535 [119,] -60815.3023 -30759.6203 [120,] -66764.7204 -60815.3023 [121,] -77484.1720 -66764.7204 [122,] -48535.5660 -77484.1720 [123,] -53364.5011 -48535.5660 [124,] -60449.4067 -53364.5011 [125,] -76767.6583 -60449.4067 [126,] -76620.2602 -76767.6583 [127,] -80177.5438 -76620.2602 [128,] 34254.9020 -80177.5438 [129,] -43846.7676 34254.9020 [130,] -98713.5407 -43846.7676 [131,] -76824.0443 -98713.5407 [132,] -51274.0349 -76824.0443 [133,] -53391.9593 -51274.0349 [134,] -79977.1345 -53391.9593 [135,] 24431.4312 -79977.1345 [136,] -57312.2709 24431.4312 [137,] 34214.6736 -57312.2709 [138,] -43783.5052 34214.6736 [139,] -104469.8502 -43783.5052 [140,] -37030.7118 -104469.8502 [141,] -27701.4535 -37030.7118 [142,] -42646.5360 -27701.4535 [143,] -38526.2130 -42646.5360 [144,] 227.5167 -38526.2130 [145,] -51326.5332 227.5167 [146,] -52433.6198 -51326.5332 [147,] -100544.0350 -52433.6198 [148,] -98387.3333 -100544.0350 [149,] -98393.4294 -98387.3333 [150,] -101706.6340 -98393.4294 [151,] 212378.0938 -101706.6340 [152,] 368592.9577 212378.0938 [153,] 202838.8708 368592.9577 [154,] 203257.1363 202838.8708 [155,] -68244.6674 203257.1363 [156,] -60957.7973 -68244.6674 [157,] -76997.0518 -60957.7973 [158,] -3073.2985 -76997.0518 [159,] -53368.3318 -3073.2985 [160,] -79766.9541 -53368.3318 [161,] -77029.7275 -79766.9541 [162,] -78834.9155 -77029.7275 [163,] -62439.7263 -78834.9155 [164,] -58091.1291 -62439.7263 [165,] -85013.8242 -58091.1291 [166,] -81764.6926 -85013.8242 [167,] -79026.4824 -81764.6926 [168,] -72533.5203 -79026.4824 [169,] -80192.9086 -72533.5203 [170,] -25461.7417 -80192.9086 [171,] -29274.6169 -25461.7417 [172,] 96401.2540 -29274.6169 [173,] 85379.0730 96401.2540 [174,] 21518.2627 85379.0730 [175,] -43974.4263 21518.2627 [176,] -38486.4860 -43974.4263 [177,] 6891.7816 -38486.4860 [178,] 5181.5796 6891.7816 [179,] 9146.0670 5181.5796 [180,] -38435.9034 9146.0670 [181,] -394.3779 -38435.9034 [182,] -37486.4399 -394.3779 [183,] -50611.7184 -37486.4399 [184,] -52495.9551 -50611.7184 [185,] -99483.0970 -52495.9551 [186,] -99600.9712 -99483.0970 [187,] 207276.0284 -99600.9712 [188,] 46720.1004 207276.0284 [189,] 203260.7656 46720.1004 [190,] -77415.9932 203260.7656 [191,] 8207.0563 -77415.9932 [192,] -53388.8923 8207.0563 [193,] -80187.2961 -53388.8923 [194,] -28600.0927 -80187.2961 [195,] 62536.5457 -28600.0927 [196,] 61927.4653 62536.5457 [197,] 24464.3847 61927.4653 [198,] -41906.6927 24464.3847 [199,] -105490.8002 -41906.6927 [200,] -92169.2391 -105490.8002 [201,] -104082.7742 -92169.2391 [202,] -106088.0254 -104082.7742 [203,] -102462.1405 -106088.0254 [204,] -91421.0283 -102462.1405 [205,] -103858.7786 -91421.0283 [206,] -106424.4940 -103858.7786 [207,] -95004.2818 -106424.4940 [208,] -34151.0419 -95004.2818 [209,] -14680.6530 -34151.0419 [210,] -36481.1758 -14680.6530 [211,] -38461.4673 -36481.1758 [212,] -6314.7240 -38461.4673 [213,] -5457.9028 -6314.7240 [214,] -4641.2721 -5457.9028 [215,] 24432.1750 -4641.2721 [216,] 6700.0320 24432.1750 [217,] 56301.5034 6700.0320 [218,] -49972.0866 56301.5034 [219,] -13605.4337 -49972.0866 [220,] 65617.6260 -13605.4337 [221,] -43257.0152 65617.6260 [222,] -2975.6361 -43257.0152 [223,] -52489.8522 -2975.6361 [224,] -100628.8139 -52489.8522 [225,] -98786.4913 -100628.8139 [226,] -102671.3692 -98786.4913 [227,] 213550.8750 -102671.3692 [228,] -47431.5150 213550.8750 [229,] -56480.3045 -47431.5150 [230,] -61411.8271 -56480.3045 [231,] -77444.9130 -61411.8271 [232,] -53270.5395 -77444.9130 [233,] -63328.5218 -53270.5395 [234,] -49819.4800 -63328.5218 [235,] -53396.5585 -49819.4800 [236,] -78906.1370 -53396.5585 [237,] -80189.3932 -78906.1370 [238,] -21281.1783 -80189.3932 [239,] 52445.0103 -21281.1783 [240,] 12826.6985 52445.0103 [241,] -43889.7870 12826.6985 [242,] -38459.5619 -43889.7870 [243,] -387.3924 -38459.5619 [244,] -52504.2159 -387.3924 [245,] 539836.2739 -52504.2159 [246,] 235437.0846 539836.2739 [247,] 541138.0297 235437.0846 [248,] 329130.6501 541138.0297 [249,] 167075.3061 329130.6501 [250,] 228098.5024 167075.3061 [251,] 129210.4649 228098.5024 [252,] 203260.3915 129210.4649 [253,] -72864.4329 203260.3915 [254,] -77457.9396 -72864.4329 [255,] -47237.6400 -77457.9396 [256,] -45447.7206 -47237.6400 [257,] -53402.3328 -45447.7206 [258,] -80106.7779 -53402.3328 [259,] -28027.5395 -80106.7779 [260,] -21641.5347 -28027.5395 [261,] -1540.8844 -21641.5347 [262,] -29949.9993 -1540.8844 [263,] -16941.1908 -29949.9993 [264,] -28783.6810 -16941.1908 [265,] -34602.8242 -28783.6810 [266,] -40045.0405 -34602.8242 [267,] -25598.5392 -40045.0405 [268,] -28083.8066 -25598.5392 [269,] 37485.2802 -28083.8066 [270,] 10503.5281 37485.2802 [271,] 26393.0639 10503.5281 [272,] 19611.9413 26393.0639 [273,] -42834.3117 19611.9413 [274,] -101985.9126 -42834.3117 [275,] -106324.5296 -101985.9126 [276,] -107202.7353 -106324.5296 [277,] -37211.0991 -107202.7353 [278,] -38387.5368 -37211.0991 [279,] -10749.8817 -38387.5368 [280,] 2383.6092 -10749.8817 [281,] 27723.8676 2383.6092 [282,] 17595.7284 27723.8676 [283,] 82149.9237 17595.7284 [284,] -52509.7685 82149.9237 [285,] 886488.6528 -52509.7685 [286,] 206698.5841 886488.6528 [287,] 199350.5118 206698.5841 [288,] 200772.2861 199350.5118 [289,] 200366.9096 200772.2861 [290,] 202837.6765 200366.9096 [291,] 269324.7964 202837.6765 [292,] 203208.8164 269324.7964 [293,] -77485.0265 203208.8164 [294,] -50497.0133 -77485.0265 [295,] -53407.7054 -50497.0133 [296,] -81286.3200 -53407.7054 [297,] -78467.8443 -81286.3200 [298,] -72312.1892 -78467.8443 [299,] -78898.0977 -72312.1892 [300,] -88856.0022 -78898.0977 [301,] 239288.1022 -88856.0022 [302,] -81018.1094 239288.1022 [303,] -80152.1766 -81018.1094 [304,] -41439.6638 -80152.1766 [305,] -22694.4401 -41439.6638 [306,] -22786.5675 -22694.4401 [307,] 9423.3411 -22786.5675 [308,] 88440.5338 9423.3411 [309,] 22176.0282 88440.5338 [310,] -38094.4522 22176.0282 [311,] -9827.2335 -38094.4522 [312,] -77250.4057 -9827.2335 [313,] -103329.5351 -77250.4057 [314,] -104043.0651 -103329.5351 [315,] -106285.8564 -104043.0651 [316,] -38756.6620 -106285.8564 [317,] 90705.3731 -38756.6620 [318,] -962.5984 90705.3731 [319,] 1176.6623 -962.5984 [320,] -21421.6578 1176.6623 [321,] -31058.0088 -21421.6578 [322,] -52517.5116 -31058.0088 [323,] 274904.7996 -52517.5116 [324,] 176858.6770 274904.7996 [325,] 203288.5967 176858.6770 [326,] -57904.2013 203288.5967 [327,] -77497.7190 -57904.2013 [328,] -51979.2262 -77497.7190 [329,] -49075.2270 -51979.2262 [330,] -89015.1393 -49075.2270 [331,] -53407.8385 -89015.1393 [332,] -80309.3676 -53407.8385 [333,] -80199.4864 -80309.3676 [334,] 34328.5342 -80199.4864 [335,] 18604.2797 34328.5342 [336,] -43002.0530 18604.2797 [337,] -105813.5158 -43002.0530 [338,] -38545.4547 -105813.5158 [339,] -397.6282 -38545.4547 [340,] -52119.1919 -397.6282 [341,] -87911.3164 -52119.1919 [342,] -82550.4707 -87911.3164 [343,] -105347.7632 -82550.4707 [344,] -101848.0355 -105347.7632 [345,] -95888.2868 -101848.0355 [346,] -100326.9524 -95888.2868 [347,] -106648.7371 -100326.9524 [348,] -101974.8185 -106648.7371 [349,] -99408.2689 -101974.8185 [350,] 114216.8324 -99408.2689 [351,] 202278.1743 114216.8324 [352,] 186032.8914 202278.1743 [353,] 213097.5189 186032.8914 [354,] 96009.3980 213097.5189 [355,] 202711.7632 96009.3980 [356,] 230794.4162 202711.7632 [357,] 314851.5059 230794.4162 [358,] 820222.6083 314851.5059 [359,] 195274.4309 820222.6083 [360,] 203133.1586 195274.4309 [361,] -67260.8431 203133.1586 [362,] -59724.8910 -67260.8431 [363,] -76213.0596 -59724.8910 [364,] -67253.9069 -76213.0596 [365,] -38779.2886 -67253.9069 [366,] -68362.3153 -38779.2886 [367,] -68942.4717 -68362.3153 [368,] -77473.9239 -68942.4717 [369,] -53424.4362 -77473.9239 [370,] -80196.2496 -53424.4362 [371,] 34221.6330 -80196.2496 [372,] -37094.1100 34221.6330 [373,] -106479.6225 -37094.1100 [374,] -103605.4195 -106479.6225 [375,] -38571.8488 -103605.4195 [376,] -396.3090 -38571.8488 [377,] -52487.6426 -396.3090 [378,] -99430.7005 -52487.6426 [379,] -92273.4987 -99430.7005 [380,] -91663.3273 -92273.4987 [381,] -96975.6662 -91663.3273 [382,] -101066.8817 -96975.6662 [383,] -90234.5815 -101066.8817 [384,] -100054.6555 -90234.5815 [385,] -98914.5719 -100054.6555 [386,] -101127.2940 -98914.5719 [387,] -96103.3282 -101127.2940 [388,] -93061.5587 -96103.3282 [389,] -88555.6314 -93061.5587 [390,] -99597.0803 -88555.6314 [391,] 213546.4967 -99597.0803 [392,] -70653.6435 213546.4967 [393,] -69434.5041 -70653.6435 [394,] -66995.4019 -69434.5041 [395,] -34945.9212 -66995.4019 [396,] -76200.7920 -34945.9212 [397,] -19756.0272 -76200.7920 [398,] -54133.6556 -19756.0272 [399,] -52052.4458 -54133.6556 [400,] -55322.1652 -52052.4458 [401,] 217267.5293 -55322.1652 [402,] -53945.3063 217267.5293 [403,] -54767.1197 -53945.3063 [404,] -97338.8189 -54767.1197 [405,] -53382.0542 -97338.8189 [406,] -74824.7040 -53382.0542 [407,] -72596.7299 -74824.7040 [408,] -80111.6225 -72596.7299 [409,] -26677.6112 -80111.6225 [410,] -24847.6612 -26677.6112 [411,] 22478.3077 -24847.6612 [412,] -23897.8391 22478.3077 [413,] -6970.4024 -23897.8391 [414,] -31477.8782 -6970.4024 [415,] -32846.2212 -31477.8782 [416,] -31924.1019 -32846.2212 [417,] -24540.5460 -31924.1019 [418,] -18844.6989 -24540.5460 [419,] -28065.5915 -18844.6989 [420,] -4028.1168 -28065.5915 [421,] 3478.4659 -4028.1168 [422,] 35208.0174 3478.4659 [423,] 17148.7214 35208.0174 [424,] 12362.0026 17148.7214 [425,] 7260.0986 12362.0026 [426,] 12969.3097 7260.0986 [427,] 13416.6999 12969.3097 [428,] 14536.8355 13416.6999 [429,] 13891.0421 14536.8355 [430,] -44675.4567 13891.0421 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -18372.0453 -51984.0821 2 -51384.7569 -18372.0453 3 -38448.7631 -51384.7569 4 -28204.3908 -38448.7631 5 -55151.6551 -28204.3908 6 -55709.7680 -55151.6551 7 -52533.0890 -55709.7680 8 -47951.0220 -52533.0890 9 -54073.2014 -47951.0220 10 -29140.1596 -54073.2014 11 -25285.6329 -29140.1596 12 -53428.4879 -25285.6329 13 -77322.8544 -53428.4879 14 -79133.9065 -77322.8544 15 -80405.2747 -79133.9065 16 -79141.9387 -80405.2747 17 -79282.3001 -79141.9387 18 -78621.9433 -79282.3001 19 -78459.2013 -78621.9433 20 -80174.0945 -78459.2013 21 -30726.2388 -80174.0945 22 95041.9604 -30726.2388 23 32759.1670 95041.9604 24 29915.8614 32759.1670 25 17150.3681 29915.8614 26 91785.8572 17150.3681 27 33861.8705 91785.8572 28 -42141.6666 33861.8705 29 -102645.2920 -42141.6666 30 -105839.0484 -102645.2920 31 -35023.5739 -105839.0484 32 16584.5218 -35023.5739 33 -24762.1198 16584.5218 34 -39826.6969 -24762.1198 35 -26299.3480 -39826.6969 36 -48066.7863 -26299.3480 37 203220.1216 -48066.7863 38 -11222.0719 203220.1216 39 -7504.4171 -11222.0719 40 -26067.1170 -7504.4171 41 -17362.2156 -26067.1170 42 -38607.7672 -17362.2156 43 -4853.7732 -38607.7672 44 1296.8835 -4853.7732 45 -195.4289 1296.8835 46 -45095.7820 -195.4289 47 -29801.9203 -45095.7820 48 2311.3611 -29801.9203 49 780.5657 2311.3611 50 -40474.6696 780.5657 51 -42507.0539 -40474.6696 52 -5850.3256 -42507.0539 53 115807.4229 -5850.3256 54 35862.4758 115807.4229 55 -42698.6495 35862.4758 56 -52010.5455 -42698.6495 57 -45048.6944 -52010.5455 58 -52510.2117 -45048.6944 59 264105.7978 -52510.2117 60 118534.8649 264105.7978 61 203754.6431 118534.8649 62 15169.0314 203754.6431 63 202766.3805 15169.0314 64 260528.1083 202766.3805 65 218013.3857 260528.1083 66 212545.9028 218013.3857 67 203053.0516 212545.9028 68 11875.1888 203053.0516 69 -76000.6433 11875.1888 70 224517.6363 -76000.6433 71 -77622.5047 224517.6363 72 -46561.4052 -77622.5047 73 -52333.5476 -46561.4052 74 -50190.7121 -52333.5476 75 -52755.5158 -50190.7121 76 -72851.0447 -52755.5158 77 130344.3494 -72851.0447 78 -49302.1727 130344.3494 79 -53421.9134 -49302.1727 80 -80181.6776 -53421.9134 81 -15917.2363 -80181.6776 82 34210.1823 -15917.2363 83 -43847.9955 34210.1823 84 -103111.7453 -43847.9955 85 -94791.3445 -103111.7453 86 -104808.1266 -94791.3445 87 -103379.0923 -104808.1266 88 -43084.6284 -103379.0923 89 146690.6464 -43084.6284 90 -73791.4766 146690.6464 91 -35904.4095 -73791.4766 92 -33232.7310 -35904.4095 93 -38570.9399 -33232.7310 94 729.7612 -38570.9399 95 2544.9957 729.7612 96 -13441.4433 2544.9957 97 47369.3150 -13441.4433 98 5654.4717 47369.3150 99 2266.4497 5654.4717 100 12495.0906 2266.4497 101 17691.9284 12495.0906 102 -104570.5179 17691.9284 103 -8747.1120 -104570.5179 104 158518.1673 -8747.1120 105 -56401.4898 158518.1673 106 -48795.5967 -56401.4898 107 -45863.7494 -48795.5967 108 -63797.5544 -45863.7494 109 -52447.5941 -63797.5544 110 -99711.4355 -52447.5941 111 207730.2291 -99711.4355 112 302876.7551 207730.2291 113 746714.1728 302876.7551 114 586239.2520 746714.1728 115 419117.0070 586239.2520 116 126309.0107 419117.0070 117 203263.1535 126309.0107 118 -30759.6203 203263.1535 119 -60815.3023 -30759.6203 120 -66764.7204 -60815.3023 121 -77484.1720 -66764.7204 122 -48535.5660 -77484.1720 123 -53364.5011 -48535.5660 124 -60449.4067 -53364.5011 125 -76767.6583 -60449.4067 126 -76620.2602 -76767.6583 127 -80177.5438 -76620.2602 128 34254.9020 -80177.5438 129 -43846.7676 34254.9020 130 -98713.5407 -43846.7676 131 -76824.0443 -98713.5407 132 -51274.0349 -76824.0443 133 -53391.9593 -51274.0349 134 -79977.1345 -53391.9593 135 24431.4312 -79977.1345 136 -57312.2709 24431.4312 137 34214.6736 -57312.2709 138 -43783.5052 34214.6736 139 -104469.8502 -43783.5052 140 -37030.7118 -104469.8502 141 -27701.4535 -37030.7118 142 -42646.5360 -27701.4535 143 -38526.2130 -42646.5360 144 227.5167 -38526.2130 145 -51326.5332 227.5167 146 -52433.6198 -51326.5332 147 -100544.0350 -52433.6198 148 -98387.3333 -100544.0350 149 -98393.4294 -98387.3333 150 -101706.6340 -98393.4294 151 212378.0938 -101706.6340 152 368592.9577 212378.0938 153 202838.8708 368592.9577 154 203257.1363 202838.8708 155 -68244.6674 203257.1363 156 -60957.7973 -68244.6674 157 -76997.0518 -60957.7973 158 -3073.2985 -76997.0518 159 -53368.3318 -3073.2985 160 -79766.9541 -53368.3318 161 -77029.7275 -79766.9541 162 -78834.9155 -77029.7275 163 -62439.7263 -78834.9155 164 -58091.1291 -62439.7263 165 -85013.8242 -58091.1291 166 -81764.6926 -85013.8242 167 -79026.4824 -81764.6926 168 -72533.5203 -79026.4824 169 -80192.9086 -72533.5203 170 -25461.7417 -80192.9086 171 -29274.6169 -25461.7417 172 96401.2540 -29274.6169 173 85379.0730 96401.2540 174 21518.2627 85379.0730 175 -43974.4263 21518.2627 176 -38486.4860 -43974.4263 177 6891.7816 -38486.4860 178 5181.5796 6891.7816 179 9146.0670 5181.5796 180 -38435.9034 9146.0670 181 -394.3779 -38435.9034 182 -37486.4399 -394.3779 183 -50611.7184 -37486.4399 184 -52495.9551 -50611.7184 185 -99483.0970 -52495.9551 186 -99600.9712 -99483.0970 187 207276.0284 -99600.9712 188 46720.1004 207276.0284 189 203260.7656 46720.1004 190 -77415.9932 203260.7656 191 8207.0563 -77415.9932 192 -53388.8923 8207.0563 193 -80187.2961 -53388.8923 194 -28600.0927 -80187.2961 195 62536.5457 -28600.0927 196 61927.4653 62536.5457 197 24464.3847 61927.4653 198 -41906.6927 24464.3847 199 -105490.8002 -41906.6927 200 -92169.2391 -105490.8002 201 -104082.7742 -92169.2391 202 -106088.0254 -104082.7742 203 -102462.1405 -106088.0254 204 -91421.0283 -102462.1405 205 -103858.7786 -91421.0283 206 -106424.4940 -103858.7786 207 -95004.2818 -106424.4940 208 -34151.0419 -95004.2818 209 -14680.6530 -34151.0419 210 -36481.1758 -14680.6530 211 -38461.4673 -36481.1758 212 -6314.7240 -38461.4673 213 -5457.9028 -6314.7240 214 -4641.2721 -5457.9028 215 24432.1750 -4641.2721 216 6700.0320 24432.1750 217 56301.5034 6700.0320 218 -49972.0866 56301.5034 219 -13605.4337 -49972.0866 220 65617.6260 -13605.4337 221 -43257.0152 65617.6260 222 -2975.6361 -43257.0152 223 -52489.8522 -2975.6361 224 -100628.8139 -52489.8522 225 -98786.4913 -100628.8139 226 -102671.3692 -98786.4913 227 213550.8750 -102671.3692 228 -47431.5150 213550.8750 229 -56480.3045 -47431.5150 230 -61411.8271 -56480.3045 231 -77444.9130 -61411.8271 232 -53270.5395 -77444.9130 233 -63328.5218 -53270.5395 234 -49819.4800 -63328.5218 235 -53396.5585 -49819.4800 236 -78906.1370 -53396.5585 237 -80189.3932 -78906.1370 238 -21281.1783 -80189.3932 239 52445.0103 -21281.1783 240 12826.6985 52445.0103 241 -43889.7870 12826.6985 242 -38459.5619 -43889.7870 243 -387.3924 -38459.5619 244 -52504.2159 -387.3924 245 539836.2739 -52504.2159 246 235437.0846 539836.2739 247 541138.0297 235437.0846 248 329130.6501 541138.0297 249 167075.3061 329130.6501 250 228098.5024 167075.3061 251 129210.4649 228098.5024 252 203260.3915 129210.4649 253 -72864.4329 203260.3915 254 -77457.9396 -72864.4329 255 -47237.6400 -77457.9396 256 -45447.7206 -47237.6400 257 -53402.3328 -45447.7206 258 -80106.7779 -53402.3328 259 -28027.5395 -80106.7779 260 -21641.5347 -28027.5395 261 -1540.8844 -21641.5347 262 -29949.9993 -1540.8844 263 -16941.1908 -29949.9993 264 -28783.6810 -16941.1908 265 -34602.8242 -28783.6810 266 -40045.0405 -34602.8242 267 -25598.5392 -40045.0405 268 -28083.8066 -25598.5392 269 37485.2802 -28083.8066 270 10503.5281 37485.2802 271 26393.0639 10503.5281 272 19611.9413 26393.0639 273 -42834.3117 19611.9413 274 -101985.9126 -42834.3117 275 -106324.5296 -101985.9126 276 -107202.7353 -106324.5296 277 -37211.0991 -107202.7353 278 -38387.5368 -37211.0991 279 -10749.8817 -38387.5368 280 2383.6092 -10749.8817 281 27723.8676 2383.6092 282 17595.7284 27723.8676 283 82149.9237 17595.7284 284 -52509.7685 82149.9237 285 886488.6528 -52509.7685 286 206698.5841 886488.6528 287 199350.5118 206698.5841 288 200772.2861 199350.5118 289 200366.9096 200772.2861 290 202837.6765 200366.9096 291 269324.7964 202837.6765 292 203208.8164 269324.7964 293 -77485.0265 203208.8164 294 -50497.0133 -77485.0265 295 -53407.7054 -50497.0133 296 -81286.3200 -53407.7054 297 -78467.8443 -81286.3200 298 -72312.1892 -78467.8443 299 -78898.0977 -72312.1892 300 -88856.0022 -78898.0977 301 239288.1022 -88856.0022 302 -81018.1094 239288.1022 303 -80152.1766 -81018.1094 304 -41439.6638 -80152.1766 305 -22694.4401 -41439.6638 306 -22786.5675 -22694.4401 307 9423.3411 -22786.5675 308 88440.5338 9423.3411 309 22176.0282 88440.5338 310 -38094.4522 22176.0282 311 -9827.2335 -38094.4522 312 -77250.4057 -9827.2335 313 -103329.5351 -77250.4057 314 -104043.0651 -103329.5351 315 -106285.8564 -104043.0651 316 -38756.6620 -106285.8564 317 90705.3731 -38756.6620 318 -962.5984 90705.3731 319 1176.6623 -962.5984 320 -21421.6578 1176.6623 321 -31058.0088 -21421.6578 322 -52517.5116 -31058.0088 323 274904.7996 -52517.5116 324 176858.6770 274904.7996 325 203288.5967 176858.6770 326 -57904.2013 203288.5967 327 -77497.7190 -57904.2013 328 -51979.2262 -77497.7190 329 -49075.2270 -51979.2262 330 -89015.1393 -49075.2270 331 -53407.8385 -89015.1393 332 -80309.3676 -53407.8385 333 -80199.4864 -80309.3676 334 34328.5342 -80199.4864 335 18604.2797 34328.5342 336 -43002.0530 18604.2797 337 -105813.5158 -43002.0530 338 -38545.4547 -105813.5158 339 -397.6282 -38545.4547 340 -52119.1919 -397.6282 341 -87911.3164 -52119.1919 342 -82550.4707 -87911.3164 343 -105347.7632 -82550.4707 344 -101848.0355 -105347.7632 345 -95888.2868 -101848.0355 346 -100326.9524 -95888.2868 347 -106648.7371 -100326.9524 348 -101974.8185 -106648.7371 349 -99408.2689 -101974.8185 350 114216.8324 -99408.2689 351 202278.1743 114216.8324 352 186032.8914 202278.1743 353 213097.5189 186032.8914 354 96009.3980 213097.5189 355 202711.7632 96009.3980 356 230794.4162 202711.7632 357 314851.5059 230794.4162 358 820222.6083 314851.5059 359 195274.4309 820222.6083 360 203133.1586 195274.4309 361 -67260.8431 203133.1586 362 -59724.8910 -67260.8431 363 -76213.0596 -59724.8910 364 -67253.9069 -76213.0596 365 -38779.2886 -67253.9069 366 -68362.3153 -38779.2886 367 -68942.4717 -68362.3153 368 -77473.9239 -68942.4717 369 -53424.4362 -77473.9239 370 -80196.2496 -53424.4362 371 34221.6330 -80196.2496 372 -37094.1100 34221.6330 373 -106479.6225 -37094.1100 374 -103605.4195 -106479.6225 375 -38571.8488 -103605.4195 376 -396.3090 -38571.8488 377 -52487.6426 -396.3090 378 -99430.7005 -52487.6426 379 -92273.4987 -99430.7005 380 -91663.3273 -92273.4987 381 -96975.6662 -91663.3273 382 -101066.8817 -96975.6662 383 -90234.5815 -101066.8817 384 -100054.6555 -90234.5815 385 -98914.5719 -100054.6555 386 -101127.2940 -98914.5719 387 -96103.3282 -101127.2940 388 -93061.5587 -96103.3282 389 -88555.6314 -93061.5587 390 -99597.0803 -88555.6314 391 213546.4967 -99597.0803 392 -70653.6435 213546.4967 393 -69434.5041 -70653.6435 394 -66995.4019 -69434.5041 395 -34945.9212 -66995.4019 396 -76200.7920 -34945.9212 397 -19756.0272 -76200.7920 398 -54133.6556 -19756.0272 399 -52052.4458 -54133.6556 400 -55322.1652 -52052.4458 401 217267.5293 -55322.1652 402 -53945.3063 217267.5293 403 -54767.1197 -53945.3063 404 -97338.8189 -54767.1197 405 -53382.0542 -97338.8189 406 -74824.7040 -53382.0542 407 -72596.7299 -74824.7040 408 -80111.6225 -72596.7299 409 -26677.6112 -80111.6225 410 -24847.6612 -26677.6112 411 22478.3077 -24847.6612 412 -23897.8391 22478.3077 413 -6970.4024 -23897.8391 414 -31477.8782 -6970.4024 415 -32846.2212 -31477.8782 416 -31924.1019 -32846.2212 417 -24540.5460 -31924.1019 418 -18844.6989 -24540.5460 419 -28065.5915 -18844.6989 420 -4028.1168 -28065.5915 421 3478.4659 -4028.1168 422 35208.0174 3478.4659 423 17148.7214 35208.0174 424 12362.0026 17148.7214 425 7260.0986 12362.0026 426 12969.3097 7260.0986 427 13416.6999 12969.3097 428 14536.8355 13416.6999 429 13891.0421 14536.8355 430 -44675.4567 13891.0421 > 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/freestat/rcomp/tmp/75q431291403799.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/freestat/rcomp/tmp/85q431291403799.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/freestat/rcomp/tmp/9g0mo1291403799.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/freestat/rcomp/tmp/10g0mo1291403799.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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/11102b1291403799.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/freestat/rcomp/tmp/12511z1291403799.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/freestat/rcomp/tmp/13t2xb1291403799.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/freestat/rcomp/tmp/144bxe1291403799.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/freestat/rcomp/tmp/15pcv21291403799.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/freestat/rcomp/tmp/16tuuq1291403799.tab") + } > > try(system("convert tmp/1rhou1291403799.ps tmp/1rhou1291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/2k86x1291403799.ps tmp/2k86x1291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/3k86x1291403799.ps tmp/3k86x1291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/4k86x1291403799.ps tmp/4k86x1291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/5k86x1291403799.ps tmp/5k86x1291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/6chni1291403799.ps tmp/6chni1291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/75q431291403799.ps tmp/75q431291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/85q431291403799.ps tmp/85q431291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/9g0mo1291403799.ps tmp/9g0mo1291403799.png",intern=TRUE)) character(0) > try(system("convert tmp/10g0mo1291403799.ps tmp/10g0mo1291403799.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 20.090 3.187 20.609