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(162556 + ,1081 + ,213118 + ,6282154 + ,1 + ,5626 + ,37 + ,56289 + ,29790 + ,309 + ,81767 + ,4321023 + ,2 + ,13337 + ,138 + ,28328 + ,87550 + ,458 + ,153198 + ,4111912 + ,3 + ,8541 + ,45 + ,17936 + ,84738 + ,588 + ,-26007 + ,223193 + ,415 + ,39 + ,0 + ,4145 + ,54660 + ,302 + ,126942 + ,1491348 + ,6 + ,4276 + ,24 + ,3040 + ,42634 + ,156 + ,157214 + ,1629616 + ,5 + ,9164 + ,34 + ,2964 + ,40949 + ,481 + ,129352 + ,1398893 + ,9 + ,2493 + ,29 + ,2865 + ,45187 + ,353 + ,234817 + ,1926517 + ,4 + ,4891 + ,38 + ,2854 + ,37704 + ,452 + ,60448 + ,983660 + ,14 + ,1734 + ,21 + ,2167 + ,16275 + ,109 + ,47818 + ,1443586 + ,7 + ,11409 + ,76 + ,1974 + ,25830 + ,115 + ,245546 + ,1073089 + ,10 + ,7592 + ,34 + ,1910 + ,12679 + ,110 + ,48020 + ,984885 + ,13 + ,7135 + ,62 + ,1871 + ,18014 + ,239 + ,-1710 + ,1405225 + ,8 + ,5043 + ,67 + ,943 + ,43556 + ,247 + ,32648 + ,227132 + ,414 + ,110 + ,1 + ,929 + ,24811 + ,505 + ,95350 + ,929118 + ,15 + ,1444 + ,29 + ,822 + ,6575 + ,159 + ,151352 + ,1071292 + ,11 + ,5480 + ,133 + ,819 + ,7123 + ,109 + ,288170 + ,638830 + ,28 + ,4026 + ,62 + ,769 + ,21950 + ,519 + ,114337 + ,856956 + ,17 + ,1266 + ,30 + ,745 + ,37597 + ,248 + ,37884 + ,992426 + ,12 + ,3195 + ,21 + ,652 + ,17821 + ,373 + ,122844 + ,444477 + ,46 + ,655 + ,14 + ,643 + ,12988 + ,119 + ,82340 + ,857217 + ,16 + ,5523 + ,51 + ,601 + ,22330 + ,84 + ,79801 + ,711969 + ,20 + ,6095 + ,23 + ,446 + ,13326 + ,102 + ,165548 + ,702380 + ,21 + ,4925 + ,38 + ,436 + ,16189 + ,295 + ,116384 + ,358589 + ,89 + ,538 + ,10 + ,379 + ,7146 + ,105 + ,134028 + ,297978 + ,367 + ,933 + ,14 + ,305 + ,15824 + ,64 + ,63838 + ,585715 + ,30 + ,6027 + ,24 + ,284 + ,27664 + ,282 + ,74996 + ,657954 + ,24 + ,1624 + ,17 + ,247 + ,11920 + ,182 + ,31080 + ,209458 + ,419 + ,52 + ,1 + ,238 + ,8568 + ,37 + ,32168 + ,786690 + ,18 + ,15856 + ,68 + ,223 + ,14416 + ,361 + ,49857 + ,439798 + ,48 + ,664 + ,17 + ,220 + ,3369 + ,28 + ,87161 + ,688779 + ,23 + ,17456 + ,145 + ,217 + ,11819 + ,85 + ,106113 + ,574339 + ,31 + ,4404 + ,32 + ,199 + ,6984 + ,45 + ,80570 + ,741409 + ,19 + ,12031 + ,78 + ,195 + ,4519 + ,49 + ,102129 + ,597793 + ,29 + ,8118 + ,88 + ,163 + ,2220 + ,22 + ,301670 + ,644190 + ,26 + ,20190 + ,200 + ,154 + ,18562 + ,155 + ,102313 + ,377934 + ,72 + ,1148 + ,10 + ,143 + ,10327 + ,91 + ,88577 + ,640273 + ,27 + ,4838 + ,43 + ,135 + ,5336 + ,81 + ,112477 + ,697458 + ,22 + ,6141 + ,93 + ,132 + ,2365 + ,79 + ,191778 + ,550608 + ,32 + ,4438 + ,148 + ,120 + ,4069 + ,145 + ,79804 + ,207393 + ,420 + ,51 + ,2 + ,119 + ,8636 + ,855 + ,128294 + ,301607 + ,359 + ,119 + ,12 + ,108 + ,13718 + ,61 + ,96448 + ,345783 + ,100 + ,2390 + ,11 + ,104 + ,4525 + ,226 + ,93811 + ,501749 + ,37 + ,1335 + ,67 + ,101 + ,6869 + ,105 + ,117520 + ,379983 + ,68 + ,1714 + ,26 + ,90 + ,4628 + ,62 + ,69159 + ,387475 + ,64 + ,3024 + ,41 + ,85 + ,3689 + ,25 + ,101792 + ,377305 + ,74 + ,7092 + ,48 + ,74 + ,4891 + ,217 + ,210568 + ,370837 + ,79 + ,787 + ,35 + ,73 + ,7489 + ,322 + ,136996 + ,430866 + ,50 + ,717 + ,31 + ,70 + ,4901 + ,84 + ,121920 + ,469107 + ,41 + ,3204 + ,55 + ,67 + ,2284 + ,33 + ,76403 + ,194493 + ,423 + ,-167 + ,-2 + ,66 + ,3160 + ,108 + ,108094 + ,530670 + ,33 + ,3062 + ,105 + ,66 + ,4150 + ,150 + ,134759 + ,518365 + ,35 + ,2122 + ,77 + ,66 + ,7285 + ,115 + ,188873 + ,491303 + ,39 + ,2533 + ,40 + ,59 + ,1134 + ,162 + ,146216 + ,527021 + ,34 + ,2019 + ,288 + ,58 + ,4658 + ,158 + ,156608 + ,233773 + ,411 + ,214 + ,7 + ,58 + ,2384 + ,97 + ,61348 + ,405972 + ,58 + ,2123 + ,86 + ,54 + ,3748 + ,9 + ,50350 + ,652925 + ,25 + ,50325 + ,121 + ,53 + ,5371 + ,66 + ,87720 + ,446211 + ,45 + ,3730 + ,46 + ,49 + ,1285 + ,107 + ,99489 + ,341340 + ,106 + ,1321 + ,110 + ,49 + ,9327 + ,101 + ,87419 + ,387699 + ,63 + ,1858 + ,20 + ,44 + ,5565 + ,47 + ,94355 + ,493408 + ,38 + ,6243 + ,53 + ,39 + ,1528 + ,38 + ,60326 + ,146494 + ,426 + ,-1408 + ,-35 + ,38 + ,3122 + ,34 + ,94670 + ,414462 + ,54 + ,6308 + ,69 + ,38 + ,7561 + ,87 + ,82425 + ,364304 + ,86 + ,1889 + ,22 + ,37 + ,2675 + ,79 + ,59017 + ,355178 + ,93 + ,1964 + ,58 + ,36 + ,13253 + ,947 + ,90829 + ,357760 + ,90 + ,167 + ,12 + ,35 + ,880 + ,74 + ,80791 + ,261216 + ,399 + ,827 + ,70 + ,34 + ,2053 + ,53 + ,100423 + ,397144 + ,60 + ,3720 + ,96 + ,33 + ,1424 + ,94 + ,131116 + ,374943 + ,77 + ,1861 + ,123 + ,32 + ,4036 + ,63 + ,100269 + ,424898 + ,53 + ,3570 + ,56 + ,32 + ,3045 + ,58 + ,27330 + ,202055 + ,422 + ,35 + ,1 + ,31 + ,5119 + ,49 + ,39039 + ,378525 + ,69 + ,3643 + ,35 + ,31 + ,1431 + ,34 + ,106885 + ,310768 + ,330 + ,3258 + ,77 + ,30 + ,554 + ,11 + ,79285 + ,325738 + ,125 + ,11431 + ,227 + ,30 + ,1975 + ,35 + ,118881 + ,394510 + ,61 + ,5557 + ,98 + ,30 + ,1765 + ,20 + ,77623 + ,247060 + ,404 + ,2353 + ,27 + ,28 + ,1012 + ,47 + ,114768 + ,368078 + ,82 + ,3576 + ,166 + ,26 + ,810 + ,43 + ,74015 + ,236761 + ,409 + ,855 + ,45 + ,26 + ,1280 + ,117 + ,69465 + ,312378 + ,326 + ,960 + ,88 + ,26 + ,666 + ,171 + ,117869 + ,339836 + ,108 + ,818 + ,210 + ,25 + ,1380 + ,26 + ,60982 + ,347385 + ,98 + ,5669 + ,107 + ,25 + ,4677 + ,75 + ,90131 + ,426280 + ,52 + ,3017 + ,48 + ,24 + ,876 + ,59 + ,138971 + ,352850 + ,96 + ,2591 + ,174 + ,24 + ,814 + ,18 + ,39625 + ,301881 + ,358 + ,5660 + ,125 + ,23 + ,514 + ,15 + ,102725 + ,377516 + ,73 + ,11834 + ,346 + ,23 + ,5692 + ,72 + ,64239 + ,357312 + ,91 + ,2185 + ,28 + ,22 + ,3642 + ,86 + ,90262 + ,458343 + ,42 + ,3004 + ,71 + ,22 + ,540 + ,14 + ,103960 + ,354228 + ,94 + ,11016 + ,286 + ,22 + ,2099 + ,64 + ,106611 + ,308636 + ,339 + ,1697 + ,52 + ,20 + ,567 + ,11 + ,103345 + ,386212 + ,65 + ,16928 + ,329 + ,20 + ,2001 + ,52 + ,95551 + ,393343 + ,62 + ,3718 + ,97 + ,19 + ,2949 + ,41 + ,82903 + ,378509 + ,70 + ,4354 + ,61 + ,18 + ,2253 + ,99 + ,63593 + ,452469 + ,43 + ,2550 + ,112 + ,17 + ,6533 + ,75 + ,126910 + ,364839 + ,85 + ,2198 + ,25 + ,16 + ,1889 + ,45 + ,37527 + ,358649 + ,88 + ,3526 + ,84 + ,14 + ,3055 + ,43 + ,60247 + ,376641 + ,75 + ,4108 + ,58 + ,14 + ,272 + ,8 + ,112995 + ,429112 + ,51 + ,28639 + ,844 + ,14 + ,1414 + ,198 + ,70184 + ,330546 + ,119 + ,659 + ,92 + ,13 + ,2564 + ,22 + ,130140 + ,403560 + ,59 + ,9253 + ,79 + ,13 + ,1383 + ,11 + ,73221 + ,317892 + ,151 + ,10717 + ,85 + ,13 + ,1261 + ,33 + ,76114 + ,307528 + ,343 + ,3258 + ,85 + ,13 + ,975 + ,23 + ,90534 + ,235133 + ,410 + ,1528 + ,36 + ,13 + ,3366 + ,80 + ,108479 + ,299243 + ,362 + ,1241 + ,29 + ,12 + ,576 + ,18 + ,113761 + ,314073 + ,313 + ,6337 + ,198 + ,12 + ,1686 + ,40 + ,68696 + ,368186 + ,81 + ,4205 + ,100 + ,11 + ,746 + ,23 + ,71561 + ,269661 + ,396 + ,3029 + ,93 + ,10 + ,3192 + ,60 + ,59831 + ,125390 + ,427 + ,-1243 + ,-23 + ,10 + ,2045 + ,20 + ,97890 + ,510834 + ,36 + ,15542 + ,152 + ,10 + ,5702 + ,61 + ,101481 + ,321896 + ,138 + ,1998 + ,21 + ,10 + ,1932 + ,36 + ,72954 + ,249898 + ,403 + ,1386 + ,26 + ,10 + ,936 + ,30 + ,67939 + ,408881 + ,56 + ,6963 + ,223 + ,9 + ,3437 + ,47 + ,48022 + ,158492 + ,425 + ,-883 + ,-12 + ,9 + ,5131 + ,71 + ,86111 + ,292154 + ,379 + ,1298 + ,18 + ,9 + ,2397 + ,14 + ,74020 + ,289513 + ,385 + ,6394 + ,37 + ,6 + ,1389 + ,9 + ,57530 + ,378049 + ,71 + ,19783 + ,128 + ,6 + ,1503 + ,39 + ,56364 + ,343466 + ,104 + ,3679 + ,95 + ,6 + ,402 + ,26 + ,84990 + ,332743 + ,115 + ,5106 + ,330 + ,5 + ,2239 + ,21 + ,88590 + ,442882 + ,47 + ,11566 + ,108 + ,5 + ,2234 + ,16 + ,77200 + ,214215 + ,418 + ,888 + ,6 + ,5 + ,837 + ,69 + ,61262 + ,315688 + ,165 + ,1677 + ,138 + ,4 + ,10579 + ,92 + ,110309 + ,375195 + ,76 + ,1904 + ,17 + ,4 + ,875 + ,14 + ,67000 + ,334280 + ,112 + ,9591 + ,153 + ,3 + ,1585 + ,107 + ,93099 + ,355864 + ,92 + ,1457 + ,98 + ,2 + ,1659 + ,29 + ,107577 + ,480382 + ,40 + ,9668 + ,169 + ,2 + ,2647 + ,37 + ,62920 + ,353058 + ,95 + ,4137 + ,58 + ,1 + ,3294 + ,23 + ,75832 + ,217193 + ,416 + ,748 + ,5 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,191 + ,0 + ,94 + ,7 + ,60793 + ,314533 + ,308 + ,16362 + ,1212 + ,0 + ,422 + ,28 + ,57935 + ,318056 + ,150 + ,4216 + ,280 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,197 + ,0 + ,34 + ,8 + ,60630 + ,314353 + ,309 + ,14294 + ,3340 + ,0 + ,1558 + ,63 + ,55637 + ,369448 + ,80 + ,2690 + ,109 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,196 + ,0 + ,43 + ,3 + ,60887 + ,312846 + ,323 + ,37615 + ,2613 + ,0 + ,645 + ,5 + ,60720 + ,312075 + ,329 + ,22415 + ,174 + ,0 + ,316 + ,9 + ,60505 + ,315009 + ,299 + ,12779 + ,364 + ,0 + ,115 + ,13 + ,60945 + ,318903 + ,145 + ,9146 + ,1037 + ,0 + ,5 + ,2 + ,60720 + ,314887 + ,304 + ,57443 + ,24814 + ,0 + ,897 + ,5 + ,60720 + ,314913 + ,303 + ,22983 + ,128 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,0 + ,389 + ,14 + ,58990 + ,325506 + ,127 + ,8965 + ,322 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,220 + ,0 + ,1002 + ,15 + ,56750 + ,298568 + ,364 + ,6571 + ,98 + ,0 + ,36 + ,3 + ,60894 + ,315834 + ,162 + ,38611 + ,3215 + ,0 + ,460 + ,15 + ,63346 + ,329784 + ,122 + ,8652 + ,282 + ,0 + ,309 + ,11 + ,56535 + ,312878 + ,322 + ,10262 + ,365 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,0 + ,9 + ,6 + ,60835 + ,314987 + ,300 + ,19165 + ,12458 + ,0 + ,271 + ,2 + ,60720 + ,325249 + ,129 + ,62625 + ,463 + ,0 + ,14 + ,1 + ,61016 + ,315877 + ,161 + ,115877 + ,8041 + ,0 + ,520 + ,10 + ,58650 + ,291650 + ,383 + ,9165 + ,176 + ,0 + ,1766 + ,73 + ,60438 + ,305959 + ,349 + ,1451 + ,60 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,208 + ,0 + ,458 + ,11 + ,58625 + ,297765 + ,368 + ,8888 + ,213 + ,0 + ,20 + ,3 + ,60938 + ,315245 + ,297 + ,38415 + ,5674 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,217 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,219 + ,0 + ,98 + ,2 + ,61490 + ,315236 + ,298 + ,57618 + ,1179 + ,0 + ,405 + ,7 + ,60845 + ,336425 + ,111 + ,19489 + ,337 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,213 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,212 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,207 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,0 + ,483 + ,27 + ,60830 + ,306268 + ,348 + ,3936 + ,220 + ,0 + ,454 + ,51 + ,63261 + ,302187 + ,357 + ,2004 + ,225 + ,0 + ,47 + ,3 + ,60720 + ,314882 + ,305 + ,38294 + ,2432 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,193 + ,0 + ,757 + ,19 + ,45689 + ,382712 + ,66 + ,9616 + ,242 + ,0 + ,4655 + ,393 + ,60720 + ,341570 + ,105 + ,360 + ,30 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,211 + ,0 + ,36 + ,4 + ,61564 + ,312412 + ,325 + ,28103 + ,3146 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,189 + ,0 + ,203 + ,9 + ,61938 + ,309596 + ,334 + ,12177 + ,541 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,205 + ,0 + ,126 + ,10 + ,60951 + ,315547 + ,171 + ,11555 + ,921 + ,0 + ,400 + ,152 + ,60720 + ,313267 + ,319 + ,745 + ,283 + ,0 + ,71 + ,1 + ,60745 + ,316176 + ,159 + ,116176 + ,1646 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,190 + ,0 + ,972 + ,34 + ,71642 + ,359335 + ,87 + ,4686 + ,164 + ,0 + ,531 + ,10 + ,71641 + ,330068 + ,120 + ,13007 + ,245 + ,0 + ,2461 + ,57 + ,55792 + ,314289 + ,310 + ,2005 + ,46 + ,0 + ,378 + ,52 + ,71873 + ,297413 + ,369 + ,1873 + ,258 + ,0 + ,23 + ,5 + ,62555 + ,314806 + ,306 + ,22961 + ,4974 + ,0 + ,638 + ,14 + ,60370 + ,333210 + ,114 + ,9515 + ,209 + ,0 + ,2300 + ,29 + ,64873 + ,352108 + ,97 + ,5245 + ,66 + ,0 + ,149 + ,5 + ,62041 + ,313332 + ,318 + ,22666 + ,761 + ,0 + ,226 + ,5 + ,65745 + ,291787 + ,382 + ,18357 + ,406 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,184 + ,0 + ,275 + ,4 + ,59500 + ,318745 + ,146 + ,29686 + ,431 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,179 + ,0 + ,141 + ,6 + ,61630 + ,315366 + ,295 + ,19228 + ,819 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,186 + ,0 + ,28 + ,2 + ,60890 + ,315688 + ,166 + ,57844 + ,4148 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,263 + ,0 + ,4980 + ,91 + ,113521 + ,409642 + ,55 + ,2304 + ,42 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,259 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,261 + ,0 + ,472 + ,20 + ,80045 + ,269587 + ,397 + ,3479 + ,148 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,243 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,250 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,247 + ,0 + ,203 + ,27 + ,50804 + ,300962 + ,361 + ,3739 + ,498 + ,0 + ,496 + ,17 + ,87390 + ,325479 + ,128 + ,7381 + ,253 + ,0 + ,10 + ,2 + ,61656 + ,316155 + ,160 + ,58077 + ,11639 + ,0 + ,63 + ,4 + ,65688 + ,318574 + ,148 + ,29643 + ,1891 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,228 + ,0 + ,1136 + ,32 + ,48522 + ,343613 + ,103 + ,4488 + ,126 + ,0 + ,265 + ,31 + ,60720 + ,306948 + ,346 + ,3450 + ,404 + ,0 + ,0 + ,60720 + ,315380 + ,210 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,0 + ,267 + ,32 + ,57640 + ,330059 + ,121 + ,4064 + ,487 + ,0 + ,474 + ,20 + ,61977 + ,288985 + ,386 + ,4449 + ,188 + ,0 + ,534 + ,7 + ,62620 + ,304485 + ,351 + ,14926 + ,196 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,201 + ,0 + ,15 + ,8 + ,60831 + ,315688 + ,164 + ,14461 + ,7687 + ,0 + ,397 + ,28 + ,60646 + ,317736 + ,152 + ,4205 + ,297 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,199 + ,0 + ,1866 + ,29 + ,56225 + ,322331 + ,135 + ,4218 + ,66 + ,0 + ,288 + ,4 + ,60510 + ,296656 + ,373 + ,24164 + ,335 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,198 + ,0 + ,3 + ,2 + ,60698 + ,315354 + ,296 + ,57677 + ,36505 + ,0 + ,468 + ,21 + ,60720 + ,312161 + ,328 + ,5341 + ,239 + ,0 + ,20 + ,2 + ,60805 + ,315576 + ,169 + ,57788 + ,5924 + ,0 + ,278 + ,26 + ,61404 + ,314922 + ,302 + ,4420 + ,413 + ,0 + ,61 + ,14 + ,60720 + ,314551 + ,307 + ,8182 + ,1874 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,215 + ,0 + ,192 + ,4 + ,65276 + ,312339 + ,327 + ,28085 + ,586 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,214 + ,0 + ,317 + ,9 + ,63915 + ,298700 + ,363 + ,10967 + ,311 + ,0 + ,738 + ,10 + ,60720 + ,321376 + ,139 + ,12138 + ,164 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,216 + ,0 + ,368 + ,17 + ,61686 + ,303230 + ,355 + ,6072 + ,281 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,204 + ,0 + ,2 + ,1 + ,60743 + ,315487 + ,172 + ,115487 + ,52257 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,202 + ,0 + ,53 + ,6 + ,60349 + ,315793 + ,163 + ,19299 + ,2197 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,195 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,194 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,200 + ,0 + ,94 + ,3 + ,61360 + ,312887 + ,321 + ,37629 + ,1200 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,237 + ,0 + ,24 + ,8 + ,59818 + ,315637 + ,168 + ,14455 + ,4772 + ,0 + ,2332 + ,4 + ,72680 + ,324385 + ,131 + ,31096 + ,53 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,236 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,233 + ,0 + ,131 + ,11 + ,61808 + ,308989 + ,338 + ,9908 + ,834 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,260 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,234 + ,0 + ,206 + ,9 + ,53110 + ,296702 + ,372 + ,10745 + ,469 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,240 + ,0 + ,167 + ,2 + ,64245 + ,307322 + ,345 + ,53661 + ,642 + ,0 + ,622 + ,73 + ,73007 + ,304376 + ,352 + ,1430 + ,168 + ,0 + ,2328 + ,94 + ,82732 + ,253588 + ,400 + ,570 + ,23 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,238 + ,0 + ,365 + ,8 + ,54820 + ,309560 + ,335 + ,13695 + ,300 + ,0 + ,364 + ,35 + ,47705 + ,298466 + ,365 + ,2813 + ,271 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,239 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,235 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,262 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,283 + ,0 + ,226 + ,12 + ,72835 + ,343929 + ,102 + ,11994 + ,636 + ,0 + ,307 + ,15 + ,58856 + ,331955 + ,116 + ,8797 + ,430 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,288 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,281 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,0 + ,188 + ,11 + ,77655 + ,381180 + ,67 + ,16471 + ,964 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,221 + ,0 + ,138 + ,6 + ,69817 + ,331420 + ,117 + ,21903 + ,954 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,286 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,285 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,265 + ,0 + ,125 + ,12 + ,60798 + ,310201 + ,332 + ,9183 + ,883 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,284 + ,0 + ,282 + ,30 + ,62452 + ,320016 + ,141 + ,4001 + ,425 + ,0 + ,335 + ,33 + ,64175 + ,320398 + ,140 + ,3648 + ,359 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,232 + ,0 + ,1324 + ,117 + ,67440 + ,291841 + ,381 + ,785 + ,69 + ,0 + ,176 + ,28 + ,68136 + ,310670 + ,331 + ,3952 + ,629 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,256 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,255 + ,0 + ,249 + ,72 + ,56726 + ,313491 + ,317 + ,1576 + ,455 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,242 + ,0 + ,333 + ,13 + ,70811 + ,331323 + ,118 + ,10102 + ,394 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,241 + ,0 + ,601 + ,6 + ,60720 + ,319210 + ,144 + ,19868 + ,198 + ,0 + ,30 + ,4 + ,62045 + ,318098 + ,149 + ,29525 + ,3912 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,254 + ,0 + ,249 + ,62 + ,54323 + ,292754 + ,378 + ,1496 + ,373 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,258 + ,0 + ,165 + ,24 + ,62841 + ,325176 + ,130 + ,5216 + ,759 + ,0 + ,453 + ,21 + ,81125 + ,365959 + ,84 + ,7903 + ,366 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,257 + ,0 + ,53 + ,14 + ,59506 + ,302409 + ,356 + ,7315 + ,1932 + ,0 + ,382 + ,21 + ,59365 + ,340968 + ,107 + ,6713 + ,369 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,229 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,230 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,224 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,223 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,320 + ,28291 + ,3808 + ,0 + ,290 + ,2 + ,58790 + ,301164 + ,360 + ,50582 + ,348 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,225 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,231 + ,0 + ,366 + ,53 + ,61808 + ,344425 + ,101 + ,2725 + ,394 + ,0 + ,2 + ,9 + ,60735 + ,315394 + ,174 + ,12822 + ,63057 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,222 + ,0 + ,209 + ,13 + ,64016 + ,316647 + ,157 + ,8973 + ,558 + ,0 + ,384 + ,22 + ,54683 + ,309836 + ,333 + ,4993 + ,286 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,227 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,244 + ,0 + ,365 + ,83 + ,87192 + ,346611 + ,99 + ,1766 + ,402 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,0 + ,49 + ,8 + ,64107 + ,322031 + ,137 + ,15254 + ,2500 + ,0 + ,3 + ,4 + ,60761 + ,315656 + ,167 + ,28914 + ,37920 + ,0 + ,133 + ,14 + ,65990 + ,339445 + ,109 + ,9960 + ,1052 + ,0 + ,32 + ,1 + ,59988 + ,314964 + ,301 + ,114964 + ,3603 + ,0 + ,368 + ,17 + ,61167 + ,297141 + ,370 + ,5714 + ,264 + ,0 + ,1 + ,6 + ,60719 + ,315372 + ,293 + ,19229 + ,84213 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,274 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,246 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,177 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,185 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,245 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,273 + ,0 + ,22 + ,2 + ,60722 + ,312502 + ,324 + ,56251 + ,5154 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,252 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,253 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,277 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,278 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,276 + ,0 + ,96 + ,5 + ,60379 + ,313729 + ,316 + ,22746 + ,1188 + ,0 + ,1 + ,2 + ,60727 + ,315388 + ,175 + ,57694 + ,117743 + ,0 + ,314 + ,5 + ,60720 + ,315371 + ,294 + ,23074 + ,367 + ,0 + ,844 + ,78 + ,60925 + ,296139 + ,375 + ,1233 + ,114 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,180 + ,0 + ,26 + ,1 + ,60896 + ,313880 + ,315 + ,113880 + ,4328 + ,0 + ,125 + ,13 + ,59734 + ,317698 + ,153 + ,9054 + ,944 + ,0 + ,304 + ,15 + ,62969 + ,295580 + ,376 + ,6372 + ,314 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,187 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,183 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,182 + ,0 + ,621 + ,48 + ,60720 + ,308256 + ,341 + ,2255 + ,174 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,181 + ,0 + ,119 + ,6 + ,59118 + ,303677 + ,354 + ,17280 + ,872 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,178 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,176 + ,0 + ,1595 + ,17 + ,60720 + ,319369 + ,143 + ,7022 + ,75 + ,0 + ,312 + ,14 + ,58598 + ,318690 + ,147 + ,8478 + ,381 + ,0 + ,60 + ,10 + ,61124 + ,314049 + ,314 + ,11405 + ,1916 + ,0 + ,587 + ,12 + ,59595 + ,325699 + ,126 + ,10475 + ,214 + ,0 + ,135 + ,2 + ,62065 + ,314210 + ,311 + ,57105 + ,847 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,272 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,269 + ,0 + ,514 + ,52 + ,78780 + ,322378 + ,133 + ,2353 + ,238 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,279 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,251 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,270 + ,0 + ,1 + ,4 + ,60722 + ,315398 + ,173 + ,28849 + ,172235 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,268 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,289 + ,0 + ,1763 + ,24 + ,61600 + ,308336 + ,340 + ,4514 + ,61 + ,0 + ,180 + ,11 + ,59635 + ,316386 + ,158 + ,10581 + ,647 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,287 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,248 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,291 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,292 + ,0 + ,218 + ,21 + ,60720 + ,315553 + ,170 + ,5503 + ,531 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,290 + ,0 + ,448 + ,40 + ,59781 + ,323361 + ,132 + ,3084 + ,275 + ,0 + ,227 + ,9 + ,76644 + ,336639 + ,110 + ,15182 + ,603 + ,0 + ,174 + ,1 + ,64820 + ,307424 + ,344 + ,107424 + ,618 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,266 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,280 + ,0 + ,121 + ,24 + ,56178 + ,295370 + ,377 + ,3974 + ,788 + ,0 + ,607 + ,11 + ,60436 + ,322340 + ,134 + ,11122 + ,201 + ,0 + ,2212 + ,14 + ,60720 + ,319864 + ,142 + ,8562 + ,54 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,271 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,264 + ,0 + ,530 + ,60 + ,73433 + ,317291 + ,155 + ,1955 + ,221 + ,0 + ,571 + ,80 + ,41477 + ,280398 + ,390 + ,1005 + ,141 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,0 + ,78 + ,16 + ,62700 + ,317330 + ,154 + ,7333 + ,1511 + ,0 + ,2489 + ,40 + ,67804 + ,238125 + ,408 + ,953 + ,15 + ,0 + ,131 + ,6 + ,59661 + ,327071 + ,123 + ,21179 + ,969 + ,0 + ,923 + ,8 + ,58620 + ,309038 + ,337 + ,13630 + ,118 + ,0 + ,72 + ,3 + ,60398 + ,314210 + ,312 + ,38070 + ,1591 + ,0 + ,572 + ,16 + ,58580 + ,307930 + ,342 + ,6746 + ,189 + ,-1 + ,397 + ,10 + ,62710 + ,322327 + ,136 + ,12233 + ,308 + ,-1 + ,450 + ,8 + ,59325 + ,292136 + ,380 + ,11517 + ,205 + ,-1 + ,622 + ,7 + ,60950 + ,263276 + ,398 + ,9039 + ,102 + ,-2 + ,694 + ,8 + ,68060 + ,367655 + ,83 + ,20957 + ,241 + ,-2 + ,3425 + ,12 + ,83620 + ,283910 + ,388 + ,6993 + ,25 + ,-3 + ,562 + ,13 + ,58456 + ,283587 + ,389 + ,6430 + ,149 + ,-3 + ,4917 + ,42 + ,52811 + ,243650 + ,406 + ,1039 + ,9 + ,-3 + ,1442 + ,118 + ,121173 + ,438493 + ,49 + ,2021 + ,165 + ,-3 + ,529 + ,9 + ,63870 + ,296261 + ,374 + ,10696 + ,182 + ,-3 + ,2126 + ,138 + ,21001 + ,230621 + ,413 + ,222 + ,14 + ,-3 + ,1061 + ,5 + ,70415 + ,304252 + ,353 + ,20850 + ,98 + ,-4 + ,776 + ,9 + ,64230 + ,333505 + ,113 + ,14834 + ,172 + ,-4 + ,611 + ,8 + ,59190 + ,296919 + ,371 + ,12115 + ,159 + ,-4 + ,1526 + ,25 + ,69351 + ,278990 + ,391 + ,3160 + ,52 + ,-4 + ,592 + ,7 + ,64270 + ,276898 + ,392 + ,10985 + ,130 + ,-4 + ,1182 + ,13 + ,70694 + ,327007 + ,124 + ,9770 + ,107 + ,-5 + ,621 + ,16 + ,68005 + ,317046 + ,156 + ,7315 + ,188 + ,-5 + ,989 + ,11 + ,58930 + ,304555 + ,350 + ,9505 + ,106 + ,-5 + ,438 + ,11 + ,58320 + ,298096 + ,366 + ,8918 + ,224 + ,-5 + ,726 + ,3 + ,69980 + ,231861 + ,412 + ,10620 + ,44 + ,-5 + ,1303 + ,61 + ,69863 + ,309422 + ,336 + ,1794 + ,84 + ,-5 + ,7419 + ,29 + ,63255 + ,286963 + ,387 + ,2999 + ,12 + ,-5 + ,1164 + ,17 + ,57320 + ,269753 + ,395 + ,4103 + ,60 + ,-6 + ,3310 + ,33 + ,75230 + ,448243 + ,44 + ,7523 + ,75 + ,-6 + ,1920 + ,15 + ,79420 + ,165404 + ,424 + ,-2306 + ,-18 + ,-6 + ,965 + ,3 + ,73490 + ,204325 + ,421 + ,1442 + ,4 + ,-8 + ,3256 + ,66 + ,35250 + ,407159 + ,57 + ,3139 + ,64 + ,-8 + ,1135 + ,17 + ,62285 + ,290476 + ,384 + ,5322 + ,80 + ,-9 + ,1270 + ,26 + ,69206 + ,275311 + ,393 + ,2897 + ,59 + ,-10 + ,661 + ,3 + ,65920 + ,246541 + ,405 + ,15514 + ,70 + ,-12 + ,1013 + ,2 + ,69770 + ,253468 + ,402 + ,26734 + ,53 + ,-14 + ,2844 + ,67 + ,72683 + ,240897 + ,407 + ,610 + ,14 + ,-15 + ,11528 + ,70 + ,-14545 + ,-83265 + ,431 + ,-4047 + ,-25 + ,-25 + ,6526 + ,26 + ,55830 + ,-42143 + ,430 + ,-9313 + ,-37 + ,-28 + ,2264 + ,24 + ,55174 + ,272713 + ,394 + ,3030 + ,32 + ,-31 + ,5109 + ,97 + ,67038 + ,215362 + ,417 + ,158 + ,3 + ,-31 + ,3999 + ,30 + ,51252 + ,42754 + ,428 + ,-5242 + ,-39 + ,-58 + ,35624 + ,223 + ,157278 + ,306275 + ,347 + ,477 + ,3 + ,-70 + ,9252 + ,48 + ,79510 + ,253537 + ,401 + ,1115 + ,6 + ,-129 + ,15236 + ,90 + ,77440 + ,372631 + ,78 + ,1918 + ,11 + ,-158 + ,18073 + ,180 + ,27284 + ,-7170 + ,429 + ,-1151 + ,-11 + ,-366) + ,dim=c(8 + ,431) + ,dimnames=list(c('Costs' + ,'Trades' + ,'Dividends' + ,'Wealth' + ,'WRank' + ,'Profit/Trades' + ,'Profit/Cost' + ,'CCScore ') + ,1:431)) > y <- array(NA,dim=c(8,431),dimnames=list(c('Costs','Trades','Dividends','Wealth','WRank','Profit/Trades','Profit/Cost','CCScore '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '8' > #'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 Costs Trades Dividends Wealth WRank Profit/Trades Profit/Cost 1 56289 162556 1081 213118 6282154 1 5626 37 2 28328 29790 309 81767 4321023 2 13337 138 3 17936 87550 458 153198 4111912 3 8541 45 4 4145 84738 588 -26007 223193 415 39 0 5 3040 54660 302 126942 1491348 6 4276 24 6 2964 42634 156 157214 1629616 5 9164 34 7 2865 40949 481 129352 1398893 9 2493 29 8 2854 45187 353 234817 1926517 4 4891 38 9 2167 37704 452 60448 983660 14 1734 21 10 1974 16275 109 47818 1443586 7 11409 76 11 1910 25830 115 245546 1073089 10 7592 34 12 1871 12679 110 48020 984885 13 7135 62 13 943 18014 239 -1710 1405225 8 5043 67 14 929 43556 247 32648 227132 414 110 1 15 822 24811 505 95350 929118 15 1444 29 16 819 6575 159 151352 1071292 11 5480 133 17 769 7123 109 288170 638830 28 4026 62 18 745 21950 519 114337 856956 17 1266 30 19 652 37597 248 37884 992426 12 3195 21 20 643 17821 373 122844 444477 46 655 14 21 601 12988 119 82340 857217 16 5523 51 22 446 22330 84 79801 711969 20 6095 23 23 436 13326 102 165548 702380 21 4925 38 24 379 16189 295 116384 358589 89 538 10 25 305 7146 105 134028 297978 367 933 14 26 284 15824 64 63838 585715 30 6027 24 27 247 27664 282 74996 657954 24 1624 17 28 238 11920 182 31080 209458 419 52 1 29 223 8568 37 32168 786690 18 15856 68 30 220 14416 361 49857 439798 48 664 17 31 217 3369 28 87161 688779 23 17456 145 32 199 11819 85 106113 574339 31 4404 32 33 195 6984 45 80570 741409 19 12031 78 34 163 4519 49 102129 597793 29 8118 88 35 154 2220 22 301670 644190 26 20190 200 36 143 18562 155 102313 377934 72 1148 10 37 135 10327 91 88577 640273 27 4838 43 38 132 5336 81 112477 697458 22 6141 93 39 120 2365 79 191778 550608 32 4438 148 40 119 4069 145 79804 207393 420 51 2 41 108 8636 855 128294 301607 359 119 12 42 104 13718 61 96448 345783 100 2390 11 43 101 4525 226 93811 501749 37 1335 67 44 90 6869 105 117520 379983 68 1714 26 45 85 4628 62 69159 387475 64 3024 41 46 74 3689 25 101792 377305 74 7092 48 47 73 4891 217 210568 370837 79 787 35 48 70 7489 322 136996 430866 50 717 31 49 67 4901 84 121920 469107 41 3204 55 50 66 2284 33 76403 194493 423 -167 -2 51 66 3160 108 108094 530670 33 3062 105 52 66 4150 150 134759 518365 35 2122 77 53 59 7285 115 188873 491303 39 2533 40 54 58 1134 162 146216 527021 34 2019 288 55 58 4658 158 156608 233773 411 214 7 56 54 2384 97 61348 405972 58 2123 86 57 53 3748 9 50350 652925 25 50325 121 58 49 5371 66 87720 446211 45 3730 46 59 49 1285 107 99489 341340 106 1321 110 60 44 9327 101 87419 387699 63 1858 20 61 39 5565 47 94355 493408 38 6243 53 62 38 1528 38 60326 146494 426 -1408 -35 63 38 3122 34 94670 414462 54 6308 69 64 37 7561 87 82425 364304 86 1889 22 65 36 2675 79 59017 355178 93 1964 58 66 35 13253 947 90829 357760 90 167 12 67 34 880 74 80791 261216 399 827 70 68 33 2053 53 100423 397144 60 3720 96 69 32 1424 94 131116 374943 77 1861 123 70 32 4036 63 100269 424898 53 3570 56 71 31 3045 58 27330 202055 422 35 1 72 31 5119 49 39039 378525 69 3643 35 73 30 1431 34 106885 310768 330 3258 77 74 30 554 11 79285 325738 125 11431 227 75 30 1975 35 118881 394510 61 5557 98 76 28 1765 20 77623 247060 404 2353 27 77 26 1012 47 114768 368078 82 3576 166 78 26 810 43 74015 236761 409 855 45 79 26 1280 117 69465 312378 326 960 88 80 25 666 171 117869 339836 108 818 210 81 25 1380 26 60982 347385 98 5669 107 82 24 4677 75 90131 426280 52 3017 48 83 24 876 59 138971 352850 96 2591 174 84 23 814 18 39625 301881 358 5660 125 85 23 514 15 102725 377516 73 11834 346 86 22 5692 72 64239 357312 91 2185 28 87 22 3642 86 90262 458343 42 3004 71 88 22 540 14 103960 354228 94 11016 286 89 20 2099 64 106611 308636 339 1697 52 90 20 567 11 103345 386212 65 16928 329 91 19 2001 52 95551 393343 62 3718 97 92 18 2949 41 82903 378509 70 4354 61 93 17 2253 99 63593 452469 43 2550 112 94 16 6533 75 126910 364839 85 2198 25 95 14 1889 45 37527 358649 88 3526 84 96 14 3055 43 60247 376641 75 4108 58 97 14 272 8 112995 429112 51 28639 844 98 13 1414 198 70184 330546 119 659 92 99 13 2564 22 130140 403560 59 9253 79 100 13 1383 11 73221 317892 151 10717 85 101 13 1261 33 76114 307528 343 3258 85 102 13 975 23 90534 235133 410 1528 36 103 12 3366 80 108479 299243 362 1241 29 104 12 576 18 113761 314073 313 6337 198 105 11 1686 40 68696 368186 81 4205 100 106 10 746 23 71561 269661 396 3029 93 107 10 3192 60 59831 125390 427 -1243 -23 108 10 2045 20 97890 510834 36 15542 152 109 10 5702 61 101481 321896 138 1998 21 110 10 1932 36 72954 249898 403 1386 26 111 9 936 30 67939 408881 56 6963 223 112 9 3437 47 48022 158492 425 -883 -12 113 9 5131 71 86111 292154 379 1298 18 114 6 2397 14 74020 289513 385 6394 37 115 6 1389 9 57530 378049 71 19783 128 116 6 1503 39 56364 343466 104 3679 95 117 5 402 26 84990 332743 115 5106 330 118 5 2239 21 88590 442882 47 11566 108 119 5 2234 16 77200 214215 418 888 6 120 4 837 69 61262 315688 165 1677 138 121 4 10579 92 110309 375195 76 1904 17 122 3 875 14 67000 334280 112 9591 153 123 2 1585 107 93099 355864 92 1457 98 124 2 1659 29 107577 480382 40 9668 169 125 1 2647 37 62920 353058 95 4137 58 126 1 3294 23 75832 217193 416 748 5 127 7 0 0 60720 315380 191 0 94 128 28 60793 314533 308 16362 1212 0 422 129 0 57935 318056 150 4216 280 0 0 130 314353 60720 315380 197 0 34 8 60630 131 369448 309 14294 3340 0 1558 63 55637 132 315380 80 2690 109 0 0 0 60720 133 37615 196 0 43 3 60887 312846 323 134 22415 2613 0 645 5 60720 312075 329 135 12779 174 0 316 9 60505 315009 299 136 9146 364 0 115 13 60945 318903 145 137 57443 1037 0 5 2 60720 314887 304 138 22983 24814 0 897 5 60720 314913 303 139 0 128 0 0 0 60720 315380 188 140 0 389 14 58990 325506 127 8965 322 141 15 0 0 60720 315380 220 0 1002 142 3 56750 298568 364 6571 98 0 36 143 15 60894 315834 162 38611 3215 0 460 144 11 63346 329784 122 8652 282 0 309 145 0 56535 312878 322 10262 365 0 0 146 314987 60720 315380 203 0 9 6 60835 147 325249 300 19165 12458 0 271 2 60720 148 315877 129 62625 463 0 14 1 61016 149 291650 161 115877 8041 0 520 10 58650 150 305959 383 9165 176 0 1766 73 60438 151 315380 349 1451 60 0 0 0 60720 152 8888 208 0 458 11 58625 297765 368 153 38415 213 0 20 3 60938 315245 297 154 0 5674 0 0 0 60720 315380 217 155 2 0 0 60720 315380 219 0 98 156 7 61490 315236 298 57618 1179 0 405 157 0 60845 336425 111 19489 337 0 0 158 315380 60720 315380 213 0 0 0 60720 159 0 212 0 0 0 60720 315380 207 160 27 0 0 60720 315380 218 0 483 161 51 60830 306268 348 3936 220 0 454 162 3 63261 302187 357 2004 225 0 47 163 0 60720 314882 305 38294 2432 0 0 164 382712 60720 315380 193 0 757 19 45689 165 341570 66 9616 242 0 4655 393 60720 166 315380 105 360 30 0 0 0 60720 167 0 192 0 0 0 60720 315380 211 168 0 36 4 61564 312412 325 28103 3146 169 9 0 0 60720 315380 189 0 203 170 0 61938 309596 334 12177 541 0 0 171 315547 60720 315380 205 0 126 10 60951 172 313267 171 11555 921 0 400 152 60720 173 316176 319 745 283 0 71 1 60745 174 315380 159 116176 1646 0 0 0 60720 175 0 209 0 0 0 60720 315380 190 176 0 972 34 71642 359335 87 4686 164 177 0 531 10 71641 330068 120 13007 245 178 0 2461 57 55792 314289 310 2005 46 179 0 378 52 71873 297413 369 1873 258 180 0 23 5 62555 314806 306 22961 4974 181 0 638 14 60370 333210 114 9515 209 182 0 2300 29 64873 352108 97 5245 66 183 0 149 5 62041 313332 318 22666 761 184 0 226 5 65745 291787 382 18357 406 185 4 0 0 60720 315380 184 0 275 186 0 59500 318745 146 29686 431 0 0 187 315366 60720 315380 179 0 141 6 61630 188 315380 295 19228 819 0 0 0 60720 189 57844 186 0 28 2 60890 315688 166 190 0 4148 0 0 0 60720 315380 263 191 0 4980 91 113521 409642 55 2304 42 192 0 0 0 60720 315380 259 0 0 193 269587 60720 315380 261 0 472 20 80045 194 315380 397 3479 148 0 0 0 60720 195 0 243 0 0 0 60720 315380 250 196 27 0 0 60720 315380 247 0 203 197 17 50804 300962 361 3739 498 0 496 198 2 87390 325479 128 7381 253 0 10 199 4 61656 316155 160 58077 11639 0 63 200 0 65688 318574 148 29643 1891 0 0 201 343613 60720 315380 228 0 1136 32 48522 202 306948 103 4488 126 0 265 31 60720 203 210 346 3450 404 0 0 60720 315380 204 267 0 0 0 60720 315380 206 0 205 474 32 57640 330059 121 4064 487 0 206 534 20 61977 288985 386 4449 188 0 207 0 7 62620 304485 351 14926 196 0 208 60831 0 60720 315380 201 0 15 8 209 60646 315688 164 14461 7687 0 397 28 210 60720 317736 152 4205 297 0 0 0 211 135 315380 199 0 1866 29 56225 322331 212 373 4218 66 0 288 4 60510 296656 213 198 24164 335 0 0 0 60720 315380 214 36505 0 3 2 60698 315354 296 57677 215 239 0 468 21 60720 312161 328 5341 216 5924 0 20 2 60805 315576 169 57788 217 413 0 278 26 61404 314922 302 4420 218 1874 0 61 14 60720 314551 307 8182 219 192 0 0 0 60720 315380 215 0 220 0 4 65276 312339 327 28085 586 0 221 63915 0 60720 315380 214 0 317 9 222 60720 298700 363 10967 311 0 738 10 223 60720 321376 139 12138 164 0 0 0 224 355 315380 216 0 368 17 61686 303230 225 204 6072 281 0 0 0 60720 315380 226 52257 0 2 1 60743 315487 172 115487 227 53 0 0 0 60720 315380 202 0 228 0 6 60349 315793 163 19299 2197 0 229 60720 0 60720 315380 195 0 0 0 230 200 315380 194 0 0 0 60720 315380 231 1200 0 94 3 61360 312887 321 37629 232 24 0 0 0 60720 315380 237 0 233 2332 8 59818 315637 168 14455 4772 0 234 0 4 72680 324385 131 31096 53 0 235 60720 0 60720 315380 236 0 0 0 236 338 315380 233 0 131 11 61808 308989 237 260 9908 834 0 0 0 60720 315380 238 206 0 0 0 60720 315380 234 0 239 0 9 53110 296702 372 10745 469 0 240 64245 0 60720 315380 240 0 167 2 241 73007 307322 345 53661 642 0 622 73 242 82732 304376 352 1430 168 0 2328 94 243 60720 253588 400 570 23 0 0 0 244 335 315380 238 0 365 8 54820 309560 245 365 13695 300 0 364 35 47705 298466 246 239 2813 271 0 0 0 60720 315380 247 0 0 0 0 60720 315380 235 0 248 60720 0 60720 315380 262 0 0 0 249 102 315380 283 0 226 12 72835 343929 250 116 11994 636 0 307 15 58856 331955 251 288 8797 430 0 0 0 60720 315380 252 0 0 0 0 60720 315380 281 0 253 77655 0 60720 315380 282 0 188 11 254 60720 381180 67 16471 964 0 0 0 255 117 315380 221 0 138 6 69817 331420 256 286 21903 954 0 0 0 60720 315380 257 0 0 0 0 60720 315380 285 0 258 60798 0 60720 315380 265 0 125 12 259 60720 310201 332 9183 883 0 0 0 260 141 315380 284 0 282 30 62452 320016 261 140 4001 425 0 335 33 64175 320398 262 232 3648 359 0 0 0 60720 315380 263 69 0 1324 117 67440 291841 381 785 264 629 0 176 28 68136 310670 331 3952 265 0 0 0 0 60720 315380 256 0 266 56726 0 60720 315380 255 0 249 72 267 60720 313491 317 1576 455 0 0 0 268 118 315380 242 0 333 13 70811 331323 269 241 10102 394 0 0 0 60720 315380 270 198 0 601 6 60720 319210 144 19868 271 3912 0 30 4 62045 318098 149 29525 272 249 0 0 0 60720 315380 254 0 273 0 62 54323 292754 378 1496 373 0 274 62841 0 60720 315380 258 0 165 24 275 81125 325176 130 5216 759 0 453 21 276 60720 365959 84 7903 366 0 0 0 277 356 315380 257 0 53 14 59506 302409 278 107 7315 1932 0 382 21 59365 340968 279 229 6713 369 0 0 0 60720 315380 280 0 0 0 0 60720 315380 230 0 281 60720 0 60720 315380 224 0 0 0 282 320 315380 223 0 30 4 60798 313164 283 360 28291 3808 0 290 2 58790 301164 284 225 50582 348 0 0 0 60720 315380 285 366 0 0 0 60720 315380 231 0 286 2 53 61808 344425 101 2725 394 0 287 0 9 60735 315394 174 12822 63057 0 288 64016 0 60720 315380 222 0 209 13 289 54683 316647 157 8973 558 0 384 22 290 60720 309836 333 4993 286 0 0 0 291 244 315380 227 0 0 0 60720 315380 292 402 0 365 83 87192 346611 99 1766 293 49 0 0 0 60720 315380 226 0 294 3 8 64107 322031 137 15254 2500 0 295 133 4 60761 315656 167 28914 37920 0 296 32 14 65990 339445 109 9960 1052 0 297 368 1 59988 314964 301 114964 3603 0 298 1 17 61167 297141 370 5714 264 0 299 0 6 60719 315372 293 19229 84213 0 300 60720 0 60720 315380 274 0 0 0 301 177 315380 246 0 0 0 60720 315380 302 0 0 0 0 60720 315380 185 0 303 60720 0 60720 315380 245 0 0 0 304 324 315380 273 0 22 2 60722 312502 305 252 56251 5154 0 0 0 60720 315380 306 0 0 0 0 60720 315380 249 0 307 60720 0 60720 315380 253 0 0 0 308 278 315380 277 0 0 0 60720 315380 309 0 0 0 0 60720 315380 275 0 310 60379 0 60720 315380 276 0 96 5 311 60727 313729 316 22746 1188 0 1 2 312 60720 315388 175 57694 117743 0 314 5 313 60925 315371 294 23074 367 0 844 78 314 60720 296139 375 1233 114 0 0 0 315 315 315380 180 0 26 1 60896 313880 316 153 113880 4328 0 125 13 59734 317698 317 376 9054 944 0 304 15 62969 295580 318 187 6372 314 0 0 0 60720 315380 319 0 0 0 0 60720 315380 183 0 320 60720 0 60720 315380 182 0 621 48 321 60720 308256 341 2255 174 0 0 0 322 354 315380 181 0 119 6 59118 303677 323 178 17280 872 0 0 0 60720 315380 324 1595 0 0 0 60720 315380 176 0 325 312 17 60720 319369 143 7022 75 0 326 60 14 58598 318690 147 8478 381 0 327 587 10 61124 314049 314 11405 1916 0 328 135 12 59595 325699 126 10475 214 0 329 0 2 62065 314210 311 57105 847 0 330 60720 0 60720 315380 272 0 0 0 331 133 315380 269 0 514 52 78780 322378 332 279 2353 238 0 0 0 60720 315380 333 0 0 0 0 60720 315380 251 0 334 60722 0 60720 315380 270 0 1 4 335 60720 315398 173 28849 172235 0 0 0 336 289 315380 268 0 0 0 60720 315380 337 61 0 1763 24 61600 308336 340 4514 338 647 0 180 11 59635 316386 158 10581 339 0 0 0 0 60720 315380 287 0 340 60720 0 60720 315380 248 0 0 0 341 292 315380 291 0 0 0 60720 315380 342 531 0 218 21 60720 315553 170 5503 343 448 0 0 0 60720 315380 290 0 344 227 40 59781 323361 132 3084 275 0 345 174 9 76644 336639 110 15182 603 0 346 0 1 64820 307424 344 107424 618 0 347 60720 0 60720 315380 266 0 0 0 348 377 315380 280 0 121 24 56178 295370 349 134 3974 788 0 607 11 60436 322340 350 142 11122 201 0 2212 14 60720 319864 351 271 8562 54 0 0 0 60720 315380 352 530 0 0 0 60720 315380 264 0 353 571 60 73433 317291 155 1955 221 0 354 0 80 41477 280398 390 1005 141 0 355 62700 0 60720 315380 267 0 78 16 356 67804 317330 154 7333 1511 0 2489 40 357 59661 238125 408 953 15 0 131 6 358 58620 327071 123 21179 969 0 923 8 359 60398 309038 337 13630 118 0 72 3 360 58580 314210 312 38070 1591 0 572 16 361 62710 307930 342 6746 189 -1 397 10 362 59325 322327 136 12233 308 -1 450 8 363 60950 292136 380 11517 205 -1 622 7 364 68060 263276 398 9039 102 -2 694 8 365 83620 367655 83 20957 241 -2 3425 12 366 58456 283910 388 6993 25 -3 562 13 367 52811 283587 389 6430 149 -3 4917 42 368 121173 243650 406 1039 9 -3 1442 118 369 63870 438493 49 2021 165 -3 529 9 370 21001 296261 374 10696 182 -3 2126 138 371 70415 230621 413 222 14 -3 1061 5 372 64230 304252 353 20850 98 -4 776 9 373 59190 333505 113 14834 172 -4 611 8 374 69351 296919 371 12115 159 -4 1526 25 375 64270 278990 391 3160 52 -4 592 7 376 70694 276898 392 10985 130 -4 1182 13 377 68005 327007 124 9770 107 -5 621 16 378 58930 317046 156 7315 188 -5 989 11 379 58320 304555 350 9505 106 -5 438 11 380 69980 298096 366 8918 224 -5 726 3 381 69863 231861 412 10620 44 -5 1303 61 382 63255 309422 336 1794 84 -5 7419 29 383 57320 286963 387 2999 12 -5 1164 17 384 75230 269753 395 4103 60 -6 3310 33 385 79420 448243 44 7523 75 -6 1920 15 386 73490 165404 424 -2306 -18 -6 965 3 387 35250 204325 421 1442 4 -8 3256 66 388 62285 407159 57 3139 64 -8 1135 17 389 69206 290476 384 5322 80 -9 1270 26 390 65920 275311 393 2897 59 -10 661 3 391 69770 246541 405 15514 70 -12 1013 2 392 72683 253468 402 26734 53 -14 2844 67 393 -14545 240897 407 610 14 -15 11528 70 394 55830 -83265 431 -4047 -25 -25 6526 26 395 55174 -42143 430 -9313 -37 -28 2264 24 396 67038 272713 394 3030 32 -31 5109 97 397 51252 215362 417 158 3 -31 3999 30 398 157278 42754 428 -5242 -39 -58 35624 223 399 79510 306275 347 477 3 -70 9252 48 400 77440 253537 401 1115 6 -129 15236 90 401 27284 372631 78 1918 11 -158 18073 180 402 213118 -7170 429 -1151 -11 -366 162556 1081 403 81767 6282154 1 5626 37 56289 29790 309 404 153198 4321023 2 13337 138 28328 87550 458 405 -26007 4111912 3 8541 45 17936 84738 588 406 126942 223193 415 39 0 4145 54660 302 407 157214 1491348 6 4276 24 3040 42634 156 408 129352 1629616 5 9164 34 2964 40949 481 409 234817 1398893 9 2493 29 2865 45187 353 410 60448 1926517 4 4891 38 2854 37704 452 411 47818 983660 14 1734 21 2167 16275 109 412 245546 1443586 7 11409 76 1974 25830 115 413 48020 1073089 10 7592 34 1910 12679 110 414 -1710 984885 13 7135 62 1871 18014 239 415 32648 1405225 8 5043 67 943 43556 247 416 95350 227132 414 110 1 929 24811 505 417 151352 929118 15 1444 29 822 6575 159 418 288170 1071292 11 5480 133 819 7123 109 419 114337 638830 28 4026 62 769 21950 519 420 37884 856956 17 1266 30 745 37597 248 421 122844 992426 12 3195 21 652 17821 373 422 82340 444477 46 655 14 643 12988 119 423 79801 857217 16 5523 51 601 22330 84 424 165548 711969 20 6095 23 446 13326 102 425 116384 702380 21 4925 38 436 16189 295 426 134028 358589 89 538 10 379 7146 105 427 63838 297978 367 933 14 305 15824 64 428 74996 585715 30 6027 24 284 27664 282 429 31080 657954 24 1624 17 247 11920 182 430 32168 209458 419 52 1 238 8568 37 431 49857 786690 18 15856 68 223 14416 361 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Costs Trades Dividends 6.150e+04 1.494e-02 1.962e-01 -1.944e-01 Wealth WRank `Profit/Trades` `Profit/Cost` -2.252e-02 -1.769e-01 -1.341e-01 -1.034e-01 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -132663 -32371 -13268 2920 311831 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 6.150e+04 6.383e+03 9.634 < 2e-16 *** Costs 1.494e-02 7.505e-03 1.990 0.04722 * Trades 1.962e-01 4.786e-02 4.099 4.98e-05 *** Dividends -1.944e-01 3.778e-02 -5.146 4.09e-07 *** Wealth -2.252e-02 7.422e-03 -3.035 0.00255 ** WRank -1.769e-01 4.136e-02 -4.278 2.33e-05 *** `Profit/Trades` -1.341e-01 5.882e-02 -2.279 0.02316 * `Profit/Cost` -1.034e-01 3.806e-02 -2.717 0.00685 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 71370 on 423 degrees of freedom Multiple R-squared: 0.1742, Adjusted R-squared: 0.1606 F-statistic: 12.75 on 7 and 423 DF, p-value: 6.96e-15 > 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,] 4.872156e-04 9.744313e-04 9.995128e-01 [2,] 4.739735e-05 9.479471e-05 9.999526e-01 [3,] 2.890755e-06 5.781511e-06 9.999971e-01 [4,] 6.794851e-07 1.358970e-06 9.999993e-01 [5,] 4.474491e-08 8.948983e-08 1.000000e+00 [6,] 3.021524e-09 6.043048e-09 1.000000e+00 [7,] 3.025228e-10 6.050457e-10 1.000000e+00 [8,] 1.861662e-11 3.723323e-11 1.000000e+00 [9,] 1.499580e-12 2.999159e-12 1.000000e+00 [10,] 1.354864e-13 2.709728e-13 1.000000e+00 [11,] 1.140312e-14 2.280624e-14 1.000000e+00 [12,] 1.352415e-15 2.704830e-15 1.000000e+00 [13,] 1.101070e-16 2.202139e-16 1.000000e+00 [14,] 8.082624e-18 1.616525e-17 1.000000e+00 [15,] 4.544945e-19 9.089889e-19 1.000000e+00 [16,] 4.814673e-20 9.629346e-20 1.000000e+00 [17,] 2.796444e-21 5.592889e-21 1.000000e+00 [18,] 1.540525e-22 3.081050e-22 1.000000e+00 [19,] 1.077094e-23 2.154188e-23 1.000000e+00 [20,] 6.525240e-25 1.305048e-24 1.000000e+00 [21,] 3.650639e-26 7.301278e-26 1.000000e+00 [22,] 2.927226e-27 5.854452e-27 1.000000e+00 [23,] 1.642894e-28 3.285787e-28 1.000000e+00 [24,] 1.067032e-29 2.134064e-29 1.000000e+00 [25,] 5.730272e-31 1.146054e-30 1.000000e+00 [26,] 4.593154e-32 9.186309e-32 1.000000e+00 [27,] 2.757153e-33 5.514305e-33 1.000000e+00 [28,] 1.463202e-34 2.926404e-34 1.000000e+00 [29,] 7.399926e-36 1.479985e-35 1.000000e+00 [30,] 3.459506e-37 6.919012e-37 1.000000e+00 [31,] 2.379736e-38 4.759472e-38 1.000000e+00 [32,] 1.702657e-39 3.405315e-39 1.000000e+00 [33,] 8.546066e-41 1.709213e-40 1.000000e+00 [34,] 5.615994e-42 1.123199e-41 1.000000e+00 [35,] 3.577194e-43 7.154388e-43 1.000000e+00 [36,] 2.550311e-44 5.100622e-44 1.000000e+00 [37,] 1.357952e-45 2.715905e-45 1.000000e+00 [38,] 6.807868e-47 1.361574e-46 1.000000e+00 [39,] 3.444232e-48 6.888464e-48 1.000000e+00 [40,] 1.506731e-49 3.013462e-49 1.000000e+00 [41,] 6.539457e-51 1.307891e-50 1.000000e+00 [42,] 2.861489e-52 5.722978e-52 1.000000e+00 [43,] 1.363264e-53 2.726529e-53 1.000000e+00 [44,] 8.176284e-55 1.635257e-54 1.000000e+00 [45,] 3.484477e-56 6.968953e-56 1.000000e+00 [46,] 1.651285e-57 3.302569e-57 1.000000e+00 [47,] 1.312822e-58 2.625644e-58 1.000000e+00 [48,] 6.885214e-60 1.377043e-59 1.000000e+00 [49,] 3.100648e-61 6.201297e-61 1.000000e+00 [50,] 1.640081e-62 3.280162e-62 1.000000e+00 [51,] 7.863243e-64 1.572649e-63 1.000000e+00 [52,] 3.270914e-65 6.541829e-65 1.000000e+00 [53,] 1.626312e-66 3.252625e-66 1.000000e+00 [54,] 7.994801e-68 1.598960e-67 1.000000e+00 [55,] 3.731800e-69 7.463601e-69 1.000000e+00 [56,] 1.433508e-70 2.867015e-70 1.000000e+00 [57,] 5.866120e-72 1.173224e-71 1.000000e+00 [58,] 2.592281e-73 5.184563e-73 1.000000e+00 [59,] 1.048853e-74 2.097706e-74 1.000000e+00 [60,] 4.693888e-76 9.387775e-76 1.000000e+00 [61,] 1.804904e-77 3.609808e-77 1.000000e+00 [62,] 8.573290e-79 1.714658e-78 1.000000e+00 [63,] 3.172206e-80 6.344412e-80 1.000000e+00 [64,] 1.150913e-81 2.301827e-81 1.000000e+00 [65,] 4.835013e-83 9.670026e-83 1.000000e+00 [66,] 1.751944e-84 3.503888e-84 1.000000e+00 [67,] 6.333092e-86 1.266618e-85 1.000000e+00 [68,] 2.307364e-87 4.614727e-87 1.000000e+00 [69,] 8.433607e-89 1.686721e-88 1.000000e+00 [70,] 2.887486e-90 5.774971e-90 1.000000e+00 [71,] 1.126416e-91 2.252831e-91 1.000000e+00 [72,] 4.596143e-93 9.192285e-93 1.000000e+00 [73,] 1.569603e-94 3.139205e-94 1.000000e+00 [74,] 5.801305e-96 1.160261e-95 1.000000e+00 [75,] 1.986152e-97 3.972304e-97 1.000000e+00 [76,] 8.223569e-99 1.644714e-98 1.000000e+00 [77,] 3.035027e-100 6.070054e-100 1.000000e+00 [78,] 9.860300e-102 1.972060e-101 1.000000e+00 [79,] 3.227620e-103 6.455240e-103 1.000000e+00 [80,] 1.041265e-104 2.082531e-104 1.000000e+00 [81,] 3.812831e-106 7.625662e-106 1.000000e+00 [82,] 1.507286e-107 3.014573e-107 1.000000e+00 [83,] 4.994613e-109 9.989225e-109 1.000000e+00 [84,] 1.971231e-110 3.942462e-110 1.000000e+00 [85,] 7.010784e-112 1.402157e-111 1.000000e+00 [86,] 2.613754e-113 5.227507e-113 1.000000e+00 [87,] 1.580557e-114 3.161113e-114 1.000000e+00 [88,] 5.166052e-116 1.033210e-115 1.000000e+00 [89,] 2.048326e-117 4.096651e-117 1.000000e+00 [90,] 7.460484e-119 1.492097e-118 1.000000e+00 [91,] 2.349852e-120 4.699703e-120 1.000000e+00 [92,] 7.306830e-122 1.461366e-121 1.000000e+00 [93,] 2.250549e-123 4.501099e-123 1.000000e+00 [94,] 7.048306e-125 1.409661e-124 1.000000e+00 [95,] 2.366799e-126 4.733597e-126 1.000000e+00 [96,] 7.407112e-128 1.481422e-127 1.000000e+00 [97,] 2.416171e-129 4.832342e-129 1.000000e+00 [98,] 8.339251e-131 1.667850e-130 1.000000e+00 [99,] 2.781414e-132 5.562828e-132 1.000000e+00 [100,] 8.258268e-134 1.651654e-133 1.000000e+00 [101,] 2.463403e-135 4.926806e-135 1.000000e+00 [102,] 7.617107e-137 1.523421e-136 1.000000e+00 [103,] 2.233555e-138 4.467110e-138 1.000000e+00 [104,] 6.345982e-140 1.269196e-139 1.000000e+00 [105,] 2.502184e-141 5.004368e-141 1.000000e+00 [106,] 7.829014e-143 1.565803e-142 1.000000e+00 [107,] 2.165002e-144 4.330004e-144 1.000000e+00 [108,] 7.218269e-146 1.443654e-145 1.000000e+00 [109,] 2.036855e-147 4.073709e-147 1.000000e+00 [110,] 5.647086e-149 1.129417e-148 1.000000e+00 [111,] 1.795480e-150 3.590961e-150 1.000000e+00 [112,] 5.462506e-152 1.092501e-151 1.000000e+00 [113,] 1.611037e-153 3.222074e-153 1.000000e+00 [114,] 5.000010e-155 1.000002e-154 1.000000e+00 [115,] 1.559521e-156 3.119042e-156 1.000000e+00 [116,] 4.237275e-158 8.474549e-158 1.000000e+00 [117,] 1.130146e-159 2.260291e-159 1.000000e+00 [118,] 5.370568e-161 1.074114e-160 1.000000e+00 [119,] 1.346883e-160 2.693766e-160 1.000000e+00 [120,] 2.418311e-161 4.836621e-161 1.000000e+00 [121,] 1.078027e-127 2.156054e-127 1.000000e+00 [122,] 2.262106e-120 4.524212e-120 1.000000e+00 [123,] 3.046019e-104 6.092037e-104 1.000000e+00 [124,] 1.585279e-104 3.170557e-104 1.000000e+00 [125,] 1.461846e-104 2.923693e-104 1.000000e+00 [126,] 7.198913e-105 1.439783e-104 1.000000e+00 [127,] 1.036071e-100 2.072142e-100 1.000000e+00 [128,] 1.266942e-101 2.533885e-101 1.000000e+00 [129,] 4.570257e-100 9.140515e-100 1.000000e+00 [130,] 4.476994e-101 8.953988e-101 1.000000e+00 [131,] 6.216593e-102 1.243319e-101 1.000000e+00 [132,] 2.429357e-102 4.858715e-102 1.000000e+00 [133,] 4.434745e-103 8.869490e-103 1.000000e+00 [134,] 1.360759e-103 2.721518e-103 1.000000e+00 [135,] 4.679531e-104 9.359062e-104 1.000000e+00 [136,] 7.791464e-103 1.558293e-102 1.000000e+00 [137,] 3.071137e-102 6.142273e-102 1.000000e+00 [138,] 5.115950e-101 1.023190e-100 1.000000e+00 [139,] 6.048062e-99 1.209612e-98 1.000000e+00 [140,] 2.878542e-97 5.757084e-97 1.000000e+00 [141,] 1.245421e-96 2.490842e-96 1.000000e+00 [142,] 7.260497e-97 1.452099e-96 1.000000e+00 [143,] 6.165313e-97 1.233063e-96 1.000000e+00 [144,] 2.606179e-96 5.212357e-96 1.000000e+00 [145,] 3.205350e-97 6.410699e-97 1.000000e+00 [146,] 9.322094e-98 1.864419e-97 1.000000e+00 [147,] 4.492477e-98 8.984953e-98 1.000000e+00 [148,] 6.253268e-98 1.250654e-97 1.000000e+00 [149,] 1.261658e-97 2.523316e-97 1.000000e+00 [150,] 1.667532e-98 3.335063e-98 1.000000e+00 [151,] 7.455330e-99 1.491066e-98 1.000000e+00 [152,] 4.483959e-99 8.967918e-99 1.000000e+00 [153,] 2.869111e-99 5.738222e-99 1.000000e+00 [154,] 2.588696e-65 5.177393e-65 1.000000e+00 [155,] 1.048779e-63 2.097559e-63 1.000000e+00 [156,] 3.538764e-62 7.077528e-62 1.000000e+00 [157,] 1.846823e-62 3.693645e-62 1.000000e+00 [158,] 5.573191e-63 1.114638e-62 1.000000e+00 [159,] 1.279660e-63 2.559320e-63 1.000000e+00 [160,] 1.993731e-63 3.987463e-63 1.000000e+00 [161,] 2.045364e-62 4.090728e-62 1.000000e+00 [162,] 1.548511e-60 3.097021e-60 1.000000e+00 [163,] 2.023644e-58 4.047288e-58 1.000000e+00 [164,] 1.985119e-56 3.970238e-56 1.000000e+00 [165,] 1.089270e-56 2.178541e-56 1.000000e+00 [166,] 2.871774e-57 5.743548e-57 1.000000e+00 [167,] 7.608622e-58 1.521724e-57 1.000000e+00 [168,] 1.972603e-58 3.945206e-58 1.000000e+00 [169,] 5.088585e-59 1.017717e-58 1.000000e+00 [170,] 3.541256e-59 7.082512e-59 1.000000e+00 [171,] 9.217331e-60 1.843466e-59 1.000000e+00 [172,] 2.379340e-60 4.758680e-60 1.000000e+00 [173,] 6.188236e-61 1.237647e-60 1.000000e+00 [174,] 1.620832e-61 3.241665e-61 1.000000e+00 [175,] 4.087979e-62 8.175959e-62 1.000000e+00 [176,] 8.587464e-62 1.717493e-61 1.000000e+00 [177,] 9.253684e-60 1.850737e-59 1.000000e+00 [178,] 1.504369e-56 3.008737e-56 1.000000e+00 [179,] 2.660566e-55 5.321131e-55 1.000000e+00 [180,] 2.083487e-55 4.166975e-55 1.000000e+00 [181,] 6.431813e-56 1.286363e-55 1.000000e+00 [182,] 1.800246e-56 3.600491e-56 1.000000e+00 [183,] 6.754679e-37 1.350936e-36 1.000000e+00 [184,] 4.311226e-33 8.622452e-33 1.000000e+00 [185,] 1.269967e-32 2.539933e-32 1.000000e+00 [186,] 5.233391e-33 1.046678e-32 1.000000e+00 [187,] 6.720969e-33 1.344194e-32 1.000000e+00 [188,] 2.065099e-32 4.130198e-32 1.000000e+00 [189,] 4.246037e-31 8.492073e-31 1.000000e+00 [190,] 4.257662e-26 8.515324e-26 1.000000e+00 [191,] 8.325134e-23 1.665027e-22 1.000000e+00 [192,] 1.998463e-17 3.996927e-17 1.000000e+00 [193,] 1.000000e+00 1.143217e-23 5.716084e-24 [194,] 1.000000e+00 2.775034e-28 1.387517e-28 [195,] 1.000000e+00 3.291978e-28 1.645989e-28 [196,] 1.000000e+00 4.796450e-28 2.398225e-28 [197,] 1.000000e+00 7.266861e-28 3.633431e-28 [198,] 1.000000e+00 6.071084e-28 3.035542e-28 [199,] 1.000000e+00 9.873087e-28 4.936544e-28 [200,] 1.000000e+00 2.073386e-27 1.036693e-27 [201,] 1.000000e+00 3.883645e-32 1.941823e-32 [202,] 1.000000e+00 1.114063e-34 5.570313e-35 [203,] 1.000000e+00 5.915626e-36 2.957813e-36 [204,] 1.000000e+00 6.629735e-36 3.314868e-36 [205,] 1.000000e+00 1.636604e-35 8.183020e-36 [206,] 1.000000e+00 3.373933e-35 1.686967e-35 [207,] 1.000000e+00 8.444152e-35 4.222076e-35 [208,] 1.000000e+00 2.104272e-34 1.052136e-34 [209,] 1.000000e+00 5.210313e-34 2.605156e-34 [210,] 1.000000e+00 9.120719e-34 4.560359e-34 [211,] 1.000000e+00 8.779980e-34 4.389990e-34 [212,] 1.000000e+00 1.499905e-33 7.499527e-34 [213,] 1.000000e+00 2.845509e-33 1.422754e-33 [214,] 1.000000e+00 2.379355e-33 1.189678e-33 [215,] 1.000000e+00 1.037226e-33 5.186131e-34 [216,] 1.000000e+00 2.631325e-34 1.315662e-34 [217,] 1.000000e+00 6.590263e-34 3.295132e-34 [218,] 1.000000e+00 1.208060e-33 6.040300e-34 [219,] 1.000000e+00 1.456349e-33 7.281747e-34 [220,] 1.000000e+00 2.287165e-33 1.143583e-33 [221,] 1.000000e+00 5.351535e-33 2.675768e-33 [222,] 1.000000e+00 1.323622e-32 6.618112e-33 [223,] 1.000000e+00 2.496505e-32 1.248252e-32 [224,] 1.000000e+00 4.946949e-32 2.473475e-32 [225,] 1.000000e+00 6.644781e-32 3.322390e-32 [226,] 1.000000e+00 1.219173e-31 6.095865e-32 [227,] 1.000000e+00 1.182558e-31 5.912792e-32 [228,] 1.000000e+00 2.887092e-31 1.443546e-31 [229,] 1.000000e+00 5.378568e-31 2.689284e-31 [230,] 1.000000e+00 6.936849e-31 3.468424e-31 [231,] 1.000000e+00 1.134481e-30 5.672406e-31 [232,] 1.000000e+00 2.087920e-30 1.043960e-30 [233,] 1.000000e+00 4.524782e-30 2.262391e-30 [234,] 1.000000e+00 8.661524e-30 4.330762e-30 [235,] 1.000000e+00 1.250741e-29 6.253706e-30 [236,] 1.000000e+00 1.938435e-29 9.692176e-30 [237,] 1.000000e+00 4.622061e-29 2.311031e-29 [238,] 1.000000e+00 6.621565e-29 3.310782e-29 [239,] 1.000000e+00 1.381801e-28 6.909007e-29 [240,] 1.000000e+00 2.315820e-28 1.157910e-28 [241,] 1.000000e+00 4.205106e-28 2.102553e-28 [242,] 1.000000e+00 9.877361e-28 4.938681e-28 [243,] 1.000000e+00 8.809676e-28 4.404838e-28 [244,] 1.000000e+00 1.851332e-27 9.256660e-28 [245,] 1.000000e+00 4.027591e-27 2.013795e-27 [246,] 1.000000e+00 7.769302e-27 3.884651e-27 [247,] 1.000000e+00 1.799996e-26 8.999978e-27 [248,] 1.000000e+00 2.578680e-26 1.289340e-26 [249,] 1.000000e+00 5.535568e-26 2.767784e-26 [250,] 1.000000e+00 1.204203e-25 6.021013e-26 [251,] 1.000000e+00 2.361229e-25 1.180614e-25 [252,] 1.000000e+00 4.701850e-25 2.350925e-25 [253,] 1.000000e+00 1.052763e-24 5.263815e-25 [254,] 1.000000e+00 2.374294e-24 1.187147e-24 [255,] 1.000000e+00 5.328741e-24 2.664371e-24 [256,] 1.000000e+00 8.244192e-24 4.122096e-24 [257,] 1.000000e+00 1.727087e-23 8.635436e-24 [258,] 1.000000e+00 3.704405e-23 1.852202e-23 [259,] 1.000000e+00 7.383002e-23 3.691501e-23 [260,] 1.000000e+00 1.606382e-22 8.031910e-23 [261,] 1.000000e+00 3.356024e-22 1.678012e-22 [262,] 1.000000e+00 7.314347e-22 3.657174e-22 [263,] 1.000000e+00 1.320704e-21 6.603522e-22 [264,] 1.000000e+00 1.798964e-21 8.994821e-22 [265,] 1.000000e+00 3.546363e-21 1.773182e-21 [266,] 1.000000e+00 7.241486e-21 3.620743e-21 [267,] 1.000000e+00 1.498788e-20 7.493940e-21 [268,] 1.000000e+00 2.881520e-20 1.440760e-20 [269,] 1.000000e+00 5.703927e-20 2.851964e-20 [270,] 1.000000e+00 1.204005e-19 6.020024e-20 [271,] 1.000000e+00 1.683924e-19 8.419622e-20 [272,] 1.000000e+00 3.455356e-19 1.727678e-19 [273,] 1.000000e+00 6.799029e-19 3.399514e-19 [274,] 1.000000e+00 1.346980e-18 6.734901e-19 [275,] 1.000000e+00 2.784629e-18 1.392315e-18 [276,] 1.000000e+00 5.344313e-18 2.672157e-18 [277,] 1.000000e+00 5.495617e-18 2.747809e-18 [278,] 1.000000e+00 7.373131e-18 3.686565e-18 [279,] 1.000000e+00 1.418686e-17 7.093432e-18 [280,] 1.000000e+00 2.793233e-17 1.396616e-17 [281,] 1.000000e+00 5.572520e-17 2.786260e-17 [282,] 1.000000e+00 1.108098e-16 5.540488e-17 [283,] 1.000000e+00 2.224757e-16 1.112378e-16 [284,] 1.000000e+00 4.050663e-16 2.025332e-16 [285,] 1.000000e+00 5.659672e-16 2.829836e-16 [286,] 1.000000e+00 1.021060e-15 5.105299e-16 [287,] 1.000000e+00 2.003960e-15 1.001980e-15 [288,] 1.000000e+00 3.308861e-15 1.654431e-15 [289,] 1.000000e+00 1.521078e-15 7.605390e-16 [290,] 1.000000e+00 2.376830e-15 1.188415e-15 [291,] 1.000000e+00 4.628653e-15 2.314327e-15 [292,] 1.000000e+00 9.095276e-15 4.547638e-15 [293,] 1.000000e+00 1.391484e-14 6.957422e-15 [294,] 1.000000e+00 2.678143e-14 1.339072e-14 [295,] 1.000000e+00 5.021128e-14 2.510564e-14 [296,] 1.000000e+00 9.677686e-14 4.838843e-14 [297,] 1.000000e+00 1.442590e-13 7.212952e-14 [298,] 1.000000e+00 2.728597e-13 1.364299e-13 [299,] 1.000000e+00 5.185559e-13 2.592779e-13 [300,] 1.000000e+00 7.570203e-13 3.785102e-13 [301,] 1.000000e+00 1.396716e-12 6.983580e-13 [302,] 1.000000e+00 2.558573e-12 1.279286e-12 [303,] 1.000000e+00 4.673584e-12 2.336792e-12 [304,] 1.000000e+00 8.475785e-12 4.237893e-12 [305,] 1.000000e+00 1.545233e-11 7.726167e-12 [306,] 1.000000e+00 2.767275e-11 1.383638e-11 [307,] 1.000000e+00 4.873261e-11 2.436630e-11 [308,] 1.000000e+00 8.633495e-11 4.316748e-11 [309,] 1.000000e+00 1.553674e-10 7.768369e-11 [310,] 1.000000e+00 2.170818e-10 1.085409e-10 [311,] 1.000000e+00 3.798850e-10 1.899425e-10 [312,] 1.000000e+00 6.682702e-10 3.341351e-10 [313,] 1.000000e+00 1.157098e-09 5.785491e-10 [314,] 1.000000e+00 2.026669e-09 1.013334e-09 [315,] 1.000000e+00 3.277419e-09 1.638710e-09 [316,] 1.000000e+00 5.188701e-09 2.594351e-09 [317,] 1.000000e+00 8.057891e-09 4.028945e-09 [318,] 1.000000e+00 1.231741e-08 6.158707e-09 [319,] 1.000000e+00 1.979156e-08 9.895778e-09 [320,] 1.000000e+00 2.890040e-08 1.445020e-08 [321,] 1.000000e+00 4.785377e-08 2.392688e-08 [322,] 1.000000e+00 7.923415e-08 3.961707e-08 [323,] 9.999999e-01 1.322164e-07 6.610821e-08 [324,] 9.999999e-01 1.854130e-07 9.270649e-08 [325,] 9.999999e-01 2.965083e-07 1.482541e-07 [326,] 9.999998e-01 4.842640e-07 2.421320e-07 [327,] 9.999996e-01 7.884213e-07 3.942106e-07 [328,] 9.999994e-01 1.273053e-06 6.365267e-07 [329,] 9.999990e-01 2.045738e-06 1.022869e-06 [330,] 9.999986e-01 2.752494e-06 1.376247e-06 [331,] 9.999978e-01 4.358926e-06 2.179463e-06 [332,] 9.999966e-01 6.879301e-06 3.439650e-06 [333,] 9.999946e-01 1.078278e-05 5.391391e-06 [334,] 9.999921e-01 1.573744e-05 7.868721e-06 [335,] 9.999884e-01 2.322320e-05 1.161160e-05 [336,] 9.999829e-01 3.414019e-05 1.707010e-05 [337,] 9.999775e-01 4.497611e-05 2.248805e-05 [338,] 9.999660e-01 6.808067e-05 3.404033e-05 [339,] 9.999491e-01 1.017672e-04 5.088360e-05 [340,] 9.999246e-01 1.507739e-04 7.538693e-05 [341,] 9.998907e-01 2.185507e-04 1.092753e-04 [342,] 9.998431e-01 3.137548e-04 1.568774e-04 [343,] 9.997966e-01 4.068960e-04 2.034480e-04 [344,] 9.997356e-01 5.287440e-04 2.643720e-04 [345,] 9.998205e-01 3.590187e-04 1.795094e-04 [346,] 9.997372e-01 5.255761e-04 2.627880e-04 [347,] 9.996167e-01 7.665325e-04 3.832662e-04 [348,] 9.994609e-01 1.078298e-03 5.391492e-04 [349,] 9.992275e-01 1.545053e-03 7.725266e-04 [350,] 9.989026e-01 2.194818e-03 1.097409e-03 [351,] 9.984520e-01 3.095974e-03 1.547987e-03 [352,] 9.979264e-01 4.147163e-03 2.073581e-03 [353,] 9.971202e-01 5.759612e-03 2.879806e-03 [354,] 9.960531e-01 7.893805e-03 3.946902e-03 [355,] 9.946816e-01 1.063686e-02 5.318429e-03 [356,] 9.927929e-01 1.441421e-02 7.207104e-03 [357,] 9.904381e-01 1.912371e-02 9.561857e-03 [358,] 9.904595e-01 1.908098e-02 9.540489e-03 [359,] 9.881386e-01 2.372276e-02 1.186138e-02 [360,] 9.866002e-01 2.679957e-02 1.339979e-02 [361,] 9.824744e-01 3.505124e-02 1.752562e-02 [362,] 9.770409e-01 4.591821e-02 2.295911e-02 [363,] 9.729759e-01 5.404829e-02 2.702414e-02 [364,] 9.652402e-01 6.951953e-02 3.475976e-02 [365,] 9.555266e-01 8.894681e-02 4.447340e-02 [366,] 9.436890e-01 1.126219e-01 5.631095e-02 [367,] 9.320979e-01 1.358041e-01 6.790207e-02 [368,] 9.280124e-01 1.439752e-01 7.198760e-02 [369,] 9.123836e-01 1.752328e-01 8.761642e-02 [370,] 9.108697e-01 1.782605e-01 8.913025e-02 [371,] 8.903547e-01 2.192907e-01 1.096453e-01 [372,] 8.770575e-01 2.458849e-01 1.229425e-01 [373,] 8.504501e-01 2.990998e-01 1.495499e-01 [374,] 8.209589e-01 3.580821e-01 1.790411e-01 [375,] 7.980734e-01 4.038532e-01 2.019266e-01 [376,] 7.684359e-01 4.631281e-01 2.315641e-01 [377,] 7.305300e-01 5.389400e-01 2.694700e-01 [378,] 7.210710e-01 5.578580e-01 2.789290e-01 [379,] 6.894066e-01 6.211868e-01 3.105934e-01 [380,] 6.534525e-01 6.930950e-01 3.465475e-01 [381,] 6.067458e-01 7.865083e-01 3.932542e-01 [382,] 5.766948e-01 8.466104e-01 4.233052e-01 [383,] 6.141980e-01 7.716041e-01 3.858020e-01 [384,] 5.625375e-01 8.749249e-01 4.374625e-01 [385,] 5.127025e-01 9.745950e-01 4.872975e-01 [386,] 4.625908e-01 9.251815e-01 5.374092e-01 [387,] 4.132957e-01 8.265914e-01 5.867043e-01 [388,] 4.306196e-01 8.612392e-01 5.693804e-01 [389,] 3.748533e-01 7.497066e-01 6.251467e-01 [390,] 3.201809e-01 6.403617e-01 6.798191e-01 [391,] 3.149654e-01 6.299308e-01 6.850346e-01 [392,] 3.228041e-01 6.456083e-01 6.771959e-01 [393,] 2.867308e-01 5.734615e-01 7.132692e-01 [394,] 2.522553e-01 5.045106e-01 7.477447e-01 [395,] 3.699063e-01 7.398126e-01 6.300937e-01 [396,] 3.098995e-01 6.197990e-01 6.901005e-01 [397,] 2.755246e-01 5.510491e-01 7.244754e-01 [398,] 2.336092e-01 4.672183e-01 7.663908e-01 [399,] 5.611366e-01 8.777269e-01 4.388634e-01 [400,] 4.843312e-01 9.686624e-01 5.156688e-01 [401,] 4.173255e-01 8.346511e-01 5.826745e-01 [402,] 8.475032e-01 3.049935e-01 1.524968e-01 [403,] 8.011517e-01 3.976967e-01 1.988483e-01 [404,] 9.256544e-01 1.486912e-01 7.434561e-02 [405,] 8.945130e-01 2.109740e-01 1.054870e-01 [406,] 8.495060e-01 3.009879e-01 1.504940e-01 [407,] 7.795463e-01 4.409074e-01 2.204537e-01 [408,] 8.115482e-01 3.769035e-01 1.884518e-01 [409,] 7.338995e-01 5.322010e-01 2.661005e-01 [410,] 6.080603e-01 7.838793e-01 3.919397e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1nyif1291299426.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/2nyif1291299426.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/3nyif1291299426.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/4f7z01291299426.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/5f7z01291299426.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 175839.7180 81350.2769 78590.8368 -58681.6234 -485.3187 9302.4862 7 8 9 10 11 12 -2342.5076 30317.0338 -25836.0535 -16436.4824 12935.2203 -27351.5918 13 14 15 16 17 18 -28865.8114 -49714.4406 -21479.9857 -6502.1158 10108.2758 -19474.5117 19 20 21 22 23 24 -31303.1213 -27201.2617 -25048.1118 -29026.3780 -12606.5837 -30624.4492 25 26 27 28 29 30 -28358.5798 -35041.5445 -32093.9150 -50629.9749 -35299.0443 -41863.8545 31 32 33 34 35 36 -26516.4798 -27325.0288 -27426.6463 -26987.8611 14509.6491 -33089.4648 37 38 39 40 41 42 -29233.3932 -23046.5555 -11125.4232 -41198.4622 -29868.7112 -34730.4736 43 44 45 46 47 48 -31775.0815 -29878.6655 -38898.3383 -32225.1461 -12126.1734 -25154.4592 49 50 51 52 53 54 -26807.4039 -42183.6733 -28103.6677 -23348.5010 -13432.7929 -20883.6733 55 56 57 58 59 60 -25724.5201 -40121.7901 -30242.1033 -33923.0665 -34249.6826 -35620.8912 61 62 63 64 65 66 -31242.9225 -46577.2192 -32908.5839 -37088.1477 -41755.8466 -36088.6335 67 68 69 70 71 72 -39710.3953 -32516.0638 -27292.1925 -31978.9497 -51577.6128 -44931.0421 73 74 75 76 77 78 -33211.3911 -37147.3803 -28738.6549 -40452.7374 -30380.3899 -41576.4536 79 80 81 82 83 84 -40775.5548 -30794.3775 -41027.9823 -34013.8013 -26149.1762 -46150.1871 85 86 87 88 89 90 -31374.2837 -40724.0988 -33255.9965 -31772.0565 -33548.7993 -30381.1736 91 92 93 94 95 96 -33561.6713 -36284.7098 -38616.2274 -28389.5734 -45647.8919 -40769.6802 97 98 99 100 101 102 -25919.1456 -40334.0261 -25876.0592 -38638.1751 -39277.6458 -38323.7626 103 104 105 106 107 108 -33486.7093 -31380.0197 -39281.1150 -41029.5240 -47182.4188 -28877.8311 109 110 111 112 113 114 -34309.0997 -41450.0565 -38122.4971 -48685.6801 -38012.9905 -39687.8624 115 116 117 118 119 120 -39134.7264 -42304.3726 -36745.2130 -32760.0580 -41500.0659 -42228.8996 121 122 123 124 125 126 -31501.2936 -39632.2950 -35201.7884 -28469.4748 -40779.5001 -41739.3726 127 128 129 130 131 132 -42536.9338 -123392.1292 -124581.4028 196396.3015 311830.6417 259655.8262 133 134 135 136 137 138 28871.8957 13620.5855 4309.3256 1218.4567 48922.0226 14283.7498 139 140 141 142 143 144 -8454.3222 -41446.6589 -42429.9080 -120671.4974 -122830.7648 -126825.1745 145 146 147 148 149 150 -123360.0469 197047.9748 268738.5417 248496.4379 215141.6004 249266.2120 151 152 153 154 155 156 259885.3422 -2191.7581 29995.1120 -8534.1588 -42536.5657 -122641.7626 157 158 159 160 161 162 -127881.8844 197428.6301 -8453.6121 -42471.9306 -122192.6253 -121559.3162 163 164 165 166 167 168 -122821.8647 263338.9146 285389.5626 260097.1803 -8452.8997 -38341.2725 169 170 171 172 173 174 -42524.0163 -122720.2949 197641.5980 256051.3986 260878.9248 237690.4336 175 176 177 178 179 180 -8455.3252 -38834.7399 -38353.1194 -43289.6112 -40496.6762 -38598.3549 181 182 183 184 185 186 -40948.8268 -40265.7478 -39206.4434 -39575.8156 -42522.4556 -124140.3273 187 188 189 190 191 192 197527.8752 256546.3170 49463.3960 -8506.6096 -29968.6716 -42541.6218 193 194 195 196 197 198 153729.5184 259503.8920 -8449.6286 -42495.7533 -120984.3518 -126412.2789 199 200 201 202 203 204 -121028.5164 -123941.4431 224608.4790 250925.1133 -21136.0257 -4027.6097 205 206 207 208 209 210 -7375.9673 -16117.7170 -11910.1747 48743.8142 -2556.0992 -4726.7520 211 212 213 214 215 216 -25193.8165 -22402.7120 -20971.2314 38181.4139 -4144.2725 7633.2241 217 218 219 220 221 222 -3526.7948 -1716.9207 -4101.4031 -8524.0974 51868.6960 -3068.9476 223 224 225 226 227 228 -3239.3226 -26256.1157 -20684.4148 59919.3474 -4242.1459 -8228.4725 229 230 231 232 233 234 48629.8409 -25291.1916 367.3699 -4266.4539 -6334.4799 -7177.8073 235 236 237 238 239 240 48630.7644 -25670.9695 -20794.1938 -4084.8560 -12260.4136 52178.4490 241 242 243 244 245 246 17391.2971 17224.9782 -4530.1306 -26547.9644 -24120.4174 -20598.7764 247 248 249 250 251 252 -4290.7220 48631.3500 -20823.1303 -19456.8388 -20670.3455 -4284.5553 253 254 255 256 257 258 65593.1409 -3258.0440 -22497.1280 -20970.8928 -4284.0191 48727.4158 259 260 261 262 263 264 -3668.5562 -24644.5995 -19750.2746 -20635.5113 -8371.7599 -3935.6539 265 266 267 268 269 270 -4287.9067 44678.0182 -5203.2582 -22371.3936 -20729.7748 -1489.4416 271 272 273 274 275 276 3168.5435 -4039.1749 -14916.0059 50776.8613 15841.2923 -4713.1975 277 278 279 280 281 282 -26647.9296 -18647.1931 -20686.2521 -4291.3923 48630.4941 -25394.1924 283 284 285 286 287 288 -23274.0781 -21341.3626 -3925.2582 -6122.9133 -1368.8576 51955.8116 289 290 291 292 293 294 -9761.8961 -4491.3185 -25253.6654 2342.6689 -4242.9285 -8425.7215 295 296 297 298 299 300 -1712.4997 -6512.2831 9168.4821 -14672.5101 2602.5664 48631.6203 301 302 303 304 305 306 -25324.3927 -4297.4249 48630.9671 -25479.1796 -22341.8481 -4288.8452 307 308 309 310 311 312 48631.1473 -25229.4741 -4285.3596 48304.0520 -1067.1399 8390.4218 313 314 315 316 317 318 -723.2047 -5029.8283 -25304.2001 -23026.9251 -22423.4944 -20712.3693 319 320 321 322 323 324 -4297.6930 48717.7619 -5004.1018 -26555.8428 -20993.7571 -2703.6314 325 326 327 328 329 330 -9751.5130 -9420.4433 -9563.7042 -7847.8881 -2360.7844 48631.5752 331 332 333 334 335 336 -22207.3850 -20545.4321 -4288.5770 48634.0779 3967.6920 -25216.7085 337 338 339 340 341 342 -5316.6458 -2439.4883 -4283.7509 48631.0347 -25218.2205 -3207.4540 343 344 345 346 347 348 -3835.3488 -9546.8239 -8141.9104 4653.3535 48631.4401 -27802.1744 349 350 351 352 353 354 -20125.2750 -20290.1680 -20610.0740 -3756.8343 -13267.8453 -14916.0444 355 356 357 358 359 360 50623.5737 3336.1802 -5267.2833 -3520.9037 -3116.9002 -154.3528 361 362 363 364 365 366 -2081.9245 -4565.2495 -2655.7431 4407.2953 21156.6960 -5919.7833 367 368 369 370 371 372 -11079.9913 56365.9726 -3716.4247 -42610.4884 5579.7210 2280.9130 373 374 375 376 377 378 -4339.0472 5913.3227 -774.2213 7283.5463 3586.8921 -5772.2392 379 380 381 382 383 384 -5883.6893 5795.9656 7069.4099 -1580.4930 -7797.1841 10873.1710 385 386 387 388 389 390 12943.1023 9120.8558 -28657.4555 -4538.7319 4504.5095 887.1391 391 392 393 394 395 396 7664.3256 12907.1600 -78049.3679 -4420.4893 -7286.6132 2671.2908 397 398 399 400 401 402 -12977.3502 98828.5992 14697.6847 14324.9978 -37005.7122 173260.5421 403 404 405 406 407 408 -58478.9410 46555.8013 -132662.9390 70131.4619 80543.8928 51361.6742 409 410 411 412 413 414 159512.3570 -23264.7612 -25458.0482 168531.3781 -25979.2891 -73759.0312 415 416 417 418 419 420 -42824.1270 33945.0739 77301.0111 212851.0134 47210.9221 -30969.8114 421 422 423 424 425 426 49687.9245 16191.7387 9682.4484 96476.0414 47629.9444 68299.9656 427 428 429 430 431 183.6247 9706.6965 -38270.7328 -31333.3773 -18298.5227 > postscript(file="/var/www/html/freestat/rcomp/tmp/6qgg31291299426.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 175839.7180 NA 1 81350.2769 175839.7180 2 78590.8368 81350.2769 3 -58681.6234 78590.8368 4 -485.3187 -58681.6234 5 9302.4862 -485.3187 6 -2342.5076 9302.4862 7 30317.0338 -2342.5076 8 -25836.0535 30317.0338 9 -16436.4824 -25836.0535 10 12935.2203 -16436.4824 11 -27351.5918 12935.2203 12 -28865.8114 -27351.5918 13 -49714.4406 -28865.8114 14 -21479.9857 -49714.4406 15 -6502.1158 -21479.9857 16 10108.2758 -6502.1158 17 -19474.5117 10108.2758 18 -31303.1213 -19474.5117 19 -27201.2617 -31303.1213 20 -25048.1118 -27201.2617 21 -29026.3780 -25048.1118 22 -12606.5837 -29026.3780 23 -30624.4492 -12606.5837 24 -28358.5798 -30624.4492 25 -35041.5445 -28358.5798 26 -32093.9150 -35041.5445 27 -50629.9749 -32093.9150 28 -35299.0443 -50629.9749 29 -41863.8545 -35299.0443 30 -26516.4798 -41863.8545 31 -27325.0288 -26516.4798 32 -27426.6463 -27325.0288 33 -26987.8611 -27426.6463 34 14509.6491 -26987.8611 35 -33089.4648 14509.6491 36 -29233.3932 -33089.4648 37 -23046.5555 -29233.3932 38 -11125.4232 -23046.5555 39 -41198.4622 -11125.4232 40 -29868.7112 -41198.4622 41 -34730.4736 -29868.7112 42 -31775.0815 -34730.4736 43 -29878.6655 -31775.0815 44 -38898.3383 -29878.6655 45 -32225.1461 -38898.3383 46 -12126.1734 -32225.1461 47 -25154.4592 -12126.1734 48 -26807.4039 -25154.4592 49 -42183.6733 -26807.4039 50 -28103.6677 -42183.6733 51 -23348.5010 -28103.6677 52 -13432.7929 -23348.5010 53 -20883.6733 -13432.7929 54 -25724.5201 -20883.6733 55 -40121.7901 -25724.5201 56 -30242.1033 -40121.7901 57 -33923.0665 -30242.1033 58 -34249.6826 -33923.0665 59 -35620.8912 -34249.6826 60 -31242.9225 -35620.8912 61 -46577.2192 -31242.9225 62 -32908.5839 -46577.2192 63 -37088.1477 -32908.5839 64 -41755.8466 -37088.1477 65 -36088.6335 -41755.8466 66 -39710.3953 -36088.6335 67 -32516.0638 -39710.3953 68 -27292.1925 -32516.0638 69 -31978.9497 -27292.1925 70 -51577.6128 -31978.9497 71 -44931.0421 -51577.6128 72 -33211.3911 -44931.0421 73 -37147.3803 -33211.3911 74 -28738.6549 -37147.3803 75 -40452.7374 -28738.6549 76 -30380.3899 -40452.7374 77 -41576.4536 -30380.3899 78 -40775.5548 -41576.4536 79 -30794.3775 -40775.5548 80 -41027.9823 -30794.3775 81 -34013.8013 -41027.9823 82 -26149.1762 -34013.8013 83 -46150.1871 -26149.1762 84 -31374.2837 -46150.1871 85 -40724.0988 -31374.2837 86 -33255.9965 -40724.0988 87 -31772.0565 -33255.9965 88 -33548.7993 -31772.0565 89 -30381.1736 -33548.7993 90 -33561.6713 -30381.1736 91 -36284.7098 -33561.6713 92 -38616.2274 -36284.7098 93 -28389.5734 -38616.2274 94 -45647.8919 -28389.5734 95 -40769.6802 -45647.8919 96 -25919.1456 -40769.6802 97 -40334.0261 -25919.1456 98 -25876.0592 -40334.0261 99 -38638.1751 -25876.0592 100 -39277.6458 -38638.1751 101 -38323.7626 -39277.6458 102 -33486.7093 -38323.7626 103 -31380.0197 -33486.7093 104 -39281.1150 -31380.0197 105 -41029.5240 -39281.1150 106 -47182.4188 -41029.5240 107 -28877.8311 -47182.4188 108 -34309.0997 -28877.8311 109 -41450.0565 -34309.0997 110 -38122.4971 -41450.0565 111 -48685.6801 -38122.4971 112 -38012.9905 -48685.6801 113 -39687.8624 -38012.9905 114 -39134.7264 -39687.8624 115 -42304.3726 -39134.7264 116 -36745.2130 -42304.3726 117 -32760.0580 -36745.2130 118 -41500.0659 -32760.0580 119 -42228.8996 -41500.0659 120 -31501.2936 -42228.8996 121 -39632.2950 -31501.2936 122 -35201.7884 -39632.2950 123 -28469.4748 -35201.7884 124 -40779.5001 -28469.4748 125 -41739.3726 -40779.5001 126 -42536.9338 -41739.3726 127 -123392.1292 -42536.9338 128 -124581.4028 -123392.1292 129 196396.3015 -124581.4028 130 311830.6417 196396.3015 131 259655.8262 311830.6417 132 28871.8957 259655.8262 133 13620.5855 28871.8957 134 4309.3256 13620.5855 135 1218.4567 4309.3256 136 48922.0226 1218.4567 137 14283.7498 48922.0226 138 -8454.3222 14283.7498 139 -41446.6589 -8454.3222 140 -42429.9080 -41446.6589 141 -120671.4974 -42429.9080 142 -122830.7648 -120671.4974 143 -126825.1745 -122830.7648 144 -123360.0469 -126825.1745 145 197047.9748 -123360.0469 146 268738.5417 197047.9748 147 248496.4379 268738.5417 148 215141.6004 248496.4379 149 249266.2120 215141.6004 150 259885.3422 249266.2120 151 -2191.7581 259885.3422 152 29995.1120 -2191.7581 153 -8534.1588 29995.1120 154 -42536.5657 -8534.1588 155 -122641.7626 -42536.5657 156 -127881.8844 -122641.7626 157 197428.6301 -127881.8844 158 -8453.6121 197428.6301 159 -42471.9306 -8453.6121 160 -122192.6253 -42471.9306 161 -121559.3162 -122192.6253 162 -122821.8647 -121559.3162 163 263338.9146 -122821.8647 164 285389.5626 263338.9146 165 260097.1803 285389.5626 166 -8452.8997 260097.1803 167 -38341.2725 -8452.8997 168 -42524.0163 -38341.2725 169 -122720.2949 -42524.0163 170 197641.5980 -122720.2949 171 256051.3986 197641.5980 172 260878.9248 256051.3986 173 237690.4336 260878.9248 174 -8455.3252 237690.4336 175 -38834.7399 -8455.3252 176 -38353.1194 -38834.7399 177 -43289.6112 -38353.1194 178 -40496.6762 -43289.6112 179 -38598.3549 -40496.6762 180 -40948.8268 -38598.3549 181 -40265.7478 -40948.8268 182 -39206.4434 -40265.7478 183 -39575.8156 -39206.4434 184 -42522.4556 -39575.8156 185 -124140.3273 -42522.4556 186 197527.8752 -124140.3273 187 256546.3170 197527.8752 188 49463.3960 256546.3170 189 -8506.6096 49463.3960 190 -29968.6716 -8506.6096 191 -42541.6218 -29968.6716 192 153729.5184 -42541.6218 193 259503.8920 153729.5184 194 -8449.6286 259503.8920 195 -42495.7533 -8449.6286 196 -120984.3518 -42495.7533 197 -126412.2789 -120984.3518 198 -121028.5164 -126412.2789 199 -123941.4431 -121028.5164 200 224608.4790 -123941.4431 201 250925.1133 224608.4790 202 -21136.0257 250925.1133 203 -4027.6097 -21136.0257 204 -7375.9673 -4027.6097 205 -16117.7170 -7375.9673 206 -11910.1747 -16117.7170 207 48743.8142 -11910.1747 208 -2556.0992 48743.8142 209 -4726.7520 -2556.0992 210 -25193.8165 -4726.7520 211 -22402.7120 -25193.8165 212 -20971.2314 -22402.7120 213 38181.4139 -20971.2314 214 -4144.2725 38181.4139 215 7633.2241 -4144.2725 216 -3526.7948 7633.2241 217 -1716.9207 -3526.7948 218 -4101.4031 -1716.9207 219 -8524.0974 -4101.4031 220 51868.6960 -8524.0974 221 -3068.9476 51868.6960 222 -3239.3226 -3068.9476 223 -26256.1157 -3239.3226 224 -20684.4148 -26256.1157 225 59919.3474 -20684.4148 226 -4242.1459 59919.3474 227 -8228.4725 -4242.1459 228 48629.8409 -8228.4725 229 -25291.1916 48629.8409 230 367.3699 -25291.1916 231 -4266.4539 367.3699 232 -6334.4799 -4266.4539 233 -7177.8073 -6334.4799 234 48630.7644 -7177.8073 235 -25670.9695 48630.7644 236 -20794.1938 -25670.9695 237 -4084.8560 -20794.1938 238 -12260.4136 -4084.8560 239 52178.4490 -12260.4136 240 17391.2971 52178.4490 241 17224.9782 17391.2971 242 -4530.1306 17224.9782 243 -26547.9644 -4530.1306 244 -24120.4174 -26547.9644 245 -20598.7764 -24120.4174 246 -4290.7220 -20598.7764 247 48631.3500 -4290.7220 248 -20823.1303 48631.3500 249 -19456.8388 -20823.1303 250 -20670.3455 -19456.8388 251 -4284.5553 -20670.3455 252 65593.1409 -4284.5553 253 -3258.0440 65593.1409 254 -22497.1280 -3258.0440 255 -20970.8928 -22497.1280 256 -4284.0191 -20970.8928 257 48727.4158 -4284.0191 258 -3668.5562 48727.4158 259 -24644.5995 -3668.5562 260 -19750.2746 -24644.5995 261 -20635.5113 -19750.2746 262 -8371.7599 -20635.5113 263 -3935.6539 -8371.7599 264 -4287.9067 -3935.6539 265 44678.0182 -4287.9067 266 -5203.2582 44678.0182 267 -22371.3936 -5203.2582 268 -20729.7748 -22371.3936 269 -1489.4416 -20729.7748 270 3168.5435 -1489.4416 271 -4039.1749 3168.5435 272 -14916.0059 -4039.1749 273 50776.8613 -14916.0059 274 15841.2923 50776.8613 275 -4713.1975 15841.2923 276 -26647.9296 -4713.1975 277 -18647.1931 -26647.9296 278 -20686.2521 -18647.1931 279 -4291.3923 -20686.2521 280 48630.4941 -4291.3923 281 -25394.1924 48630.4941 282 -23274.0781 -25394.1924 283 -21341.3626 -23274.0781 284 -3925.2582 -21341.3626 285 -6122.9133 -3925.2582 286 -1368.8576 -6122.9133 287 51955.8116 -1368.8576 288 -9761.8961 51955.8116 289 -4491.3185 -9761.8961 290 -25253.6654 -4491.3185 291 2342.6689 -25253.6654 292 -4242.9285 2342.6689 293 -8425.7215 -4242.9285 294 -1712.4997 -8425.7215 295 -6512.2831 -1712.4997 296 9168.4821 -6512.2831 297 -14672.5101 9168.4821 298 2602.5664 -14672.5101 299 48631.6203 2602.5664 300 -25324.3927 48631.6203 301 -4297.4249 -25324.3927 302 48630.9671 -4297.4249 303 -25479.1796 48630.9671 304 -22341.8481 -25479.1796 305 -4288.8452 -22341.8481 306 48631.1473 -4288.8452 307 -25229.4741 48631.1473 308 -4285.3596 -25229.4741 309 48304.0520 -4285.3596 310 -1067.1399 48304.0520 311 8390.4218 -1067.1399 312 -723.2047 8390.4218 313 -5029.8283 -723.2047 314 -25304.2001 -5029.8283 315 -23026.9251 -25304.2001 316 -22423.4944 -23026.9251 317 -20712.3693 -22423.4944 318 -4297.6930 -20712.3693 319 48717.7619 -4297.6930 320 -5004.1018 48717.7619 321 -26555.8428 -5004.1018 322 -20993.7571 -26555.8428 323 -2703.6314 -20993.7571 324 -9751.5130 -2703.6314 325 -9420.4433 -9751.5130 326 -9563.7042 -9420.4433 327 -7847.8881 -9563.7042 328 -2360.7844 -7847.8881 329 48631.5752 -2360.7844 330 -22207.3850 48631.5752 331 -20545.4321 -22207.3850 332 -4288.5770 -20545.4321 333 48634.0779 -4288.5770 334 3967.6920 48634.0779 335 -25216.7085 3967.6920 336 -5316.6458 -25216.7085 337 -2439.4883 -5316.6458 338 -4283.7509 -2439.4883 339 48631.0347 -4283.7509 340 -25218.2205 48631.0347 341 -3207.4540 -25218.2205 342 -3835.3488 -3207.4540 343 -9546.8239 -3835.3488 344 -8141.9104 -9546.8239 345 4653.3535 -8141.9104 346 48631.4401 4653.3535 347 -27802.1744 48631.4401 348 -20125.2750 -27802.1744 349 -20290.1680 -20125.2750 350 -20610.0740 -20290.1680 351 -3756.8343 -20610.0740 352 -13267.8453 -3756.8343 353 -14916.0444 -13267.8453 354 50623.5737 -14916.0444 355 3336.1802 50623.5737 356 -5267.2833 3336.1802 357 -3520.9037 -5267.2833 358 -3116.9002 -3520.9037 359 -154.3528 -3116.9002 360 -2081.9245 -154.3528 361 -4565.2495 -2081.9245 362 -2655.7431 -4565.2495 363 4407.2953 -2655.7431 364 21156.6960 4407.2953 365 -5919.7833 21156.6960 366 -11079.9913 -5919.7833 367 56365.9726 -11079.9913 368 -3716.4247 56365.9726 369 -42610.4884 -3716.4247 370 5579.7210 -42610.4884 371 2280.9130 5579.7210 372 -4339.0472 2280.9130 373 5913.3227 -4339.0472 374 -774.2213 5913.3227 375 7283.5463 -774.2213 376 3586.8921 7283.5463 377 -5772.2392 3586.8921 378 -5883.6893 -5772.2392 379 5795.9656 -5883.6893 380 7069.4099 5795.9656 381 -1580.4930 7069.4099 382 -7797.1841 -1580.4930 383 10873.1710 -7797.1841 384 12943.1023 10873.1710 385 9120.8558 12943.1023 386 -28657.4555 9120.8558 387 -4538.7319 -28657.4555 388 4504.5095 -4538.7319 389 887.1391 4504.5095 390 7664.3256 887.1391 391 12907.1600 7664.3256 392 -78049.3679 12907.1600 393 -4420.4893 -78049.3679 394 -7286.6132 -4420.4893 395 2671.2908 -7286.6132 396 -12977.3502 2671.2908 397 98828.5992 -12977.3502 398 14697.6847 98828.5992 399 14324.9978 14697.6847 400 -37005.7122 14324.9978 401 173260.5421 -37005.7122 402 -58478.9410 173260.5421 403 46555.8013 -58478.9410 404 -132662.9390 46555.8013 405 70131.4619 -132662.9390 406 80543.8928 70131.4619 407 51361.6742 80543.8928 408 159512.3570 51361.6742 409 -23264.7612 159512.3570 410 -25458.0482 -23264.7612 411 168531.3781 -25458.0482 412 -25979.2891 168531.3781 413 -73759.0312 -25979.2891 414 -42824.1270 -73759.0312 415 33945.0739 -42824.1270 416 77301.0111 33945.0739 417 212851.0134 77301.0111 418 47210.9221 212851.0134 419 -30969.8114 47210.9221 420 49687.9245 -30969.8114 421 16191.7387 49687.9245 422 9682.4484 16191.7387 423 96476.0414 9682.4484 424 47629.9444 96476.0414 425 68299.9656 47629.9444 426 183.6247 68299.9656 427 9706.6965 183.6247 428 -38270.7328 9706.6965 429 -31333.3773 -38270.7328 430 -18298.5227 -31333.3773 431 NA -18298.5227 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 81350.2769 175839.7180 [2,] 78590.8368 81350.2769 [3,] -58681.6234 78590.8368 [4,] -485.3187 -58681.6234 [5,] 9302.4862 -485.3187 [6,] -2342.5076 9302.4862 [7,] 30317.0338 -2342.5076 [8,] -25836.0535 30317.0338 [9,] -16436.4824 -25836.0535 [10,] 12935.2203 -16436.4824 [11,] -27351.5918 12935.2203 [12,] -28865.8114 -27351.5918 [13,] -49714.4406 -28865.8114 [14,] -21479.9857 -49714.4406 [15,] -6502.1158 -21479.9857 [16,] 10108.2758 -6502.1158 [17,] -19474.5117 10108.2758 [18,] -31303.1213 -19474.5117 [19,] -27201.2617 -31303.1213 [20,] -25048.1118 -27201.2617 [21,] -29026.3780 -25048.1118 [22,] -12606.5837 -29026.3780 [23,] -30624.4492 -12606.5837 [24,] -28358.5798 -30624.4492 [25,] -35041.5445 -28358.5798 [26,] -32093.9150 -35041.5445 [27,] -50629.9749 -32093.9150 [28,] -35299.0443 -50629.9749 [29,] -41863.8545 -35299.0443 [30,] -26516.4798 -41863.8545 [31,] -27325.0288 -26516.4798 [32,] -27426.6463 -27325.0288 [33,] -26987.8611 -27426.6463 [34,] 14509.6491 -26987.8611 [35,] -33089.4648 14509.6491 [36,] -29233.3932 -33089.4648 [37,] -23046.5555 -29233.3932 [38,] -11125.4232 -23046.5555 [39,] -41198.4622 -11125.4232 [40,] -29868.7112 -41198.4622 [41,] -34730.4736 -29868.7112 [42,] -31775.0815 -34730.4736 [43,] -29878.6655 -31775.0815 [44,] -38898.3383 -29878.6655 [45,] -32225.1461 -38898.3383 [46,] -12126.1734 -32225.1461 [47,] -25154.4592 -12126.1734 [48,] -26807.4039 -25154.4592 [49,] -42183.6733 -26807.4039 [50,] -28103.6677 -42183.6733 [51,] -23348.5010 -28103.6677 [52,] -13432.7929 -23348.5010 [53,] -20883.6733 -13432.7929 [54,] -25724.5201 -20883.6733 [55,] -40121.7901 -25724.5201 [56,] -30242.1033 -40121.7901 [57,] -33923.0665 -30242.1033 [58,] -34249.6826 -33923.0665 [59,] -35620.8912 -34249.6826 [60,] -31242.9225 -35620.8912 [61,] -46577.2192 -31242.9225 [62,] -32908.5839 -46577.2192 [63,] -37088.1477 -32908.5839 [64,] -41755.8466 -37088.1477 [65,] -36088.6335 -41755.8466 [66,] -39710.3953 -36088.6335 [67,] -32516.0638 -39710.3953 [68,] -27292.1925 -32516.0638 [69,] -31978.9497 -27292.1925 [70,] -51577.6128 -31978.9497 [71,] -44931.0421 -51577.6128 [72,] -33211.3911 -44931.0421 [73,] -37147.3803 -33211.3911 [74,] -28738.6549 -37147.3803 [75,] -40452.7374 -28738.6549 [76,] -30380.3899 -40452.7374 [77,] -41576.4536 -30380.3899 [78,] -40775.5548 -41576.4536 [79,] -30794.3775 -40775.5548 [80,] -41027.9823 -30794.3775 [81,] -34013.8013 -41027.9823 [82,] -26149.1762 -34013.8013 [83,] -46150.1871 -26149.1762 [84,] -31374.2837 -46150.1871 [85,] -40724.0988 -31374.2837 [86,] -33255.9965 -40724.0988 [87,] -31772.0565 -33255.9965 [88,] -33548.7993 -31772.0565 [89,] -30381.1736 -33548.7993 [90,] -33561.6713 -30381.1736 [91,] -36284.7098 -33561.6713 [92,] -38616.2274 -36284.7098 [93,] -28389.5734 -38616.2274 [94,] -45647.8919 -28389.5734 [95,] -40769.6802 -45647.8919 [96,] -25919.1456 -40769.6802 [97,] -40334.0261 -25919.1456 [98,] -25876.0592 -40334.0261 [99,] -38638.1751 -25876.0592 [100,] -39277.6458 -38638.1751 [101,] -38323.7626 -39277.6458 [102,] -33486.7093 -38323.7626 [103,] -31380.0197 -33486.7093 [104,] -39281.1150 -31380.0197 [105,] -41029.5240 -39281.1150 [106,] -47182.4188 -41029.5240 [107,] -28877.8311 -47182.4188 [108,] -34309.0997 -28877.8311 [109,] -41450.0565 -34309.0997 [110,] -38122.4971 -41450.0565 [111,] -48685.6801 -38122.4971 [112,] -38012.9905 -48685.6801 [113,] -39687.8624 -38012.9905 [114,] -39134.7264 -39687.8624 [115,] -42304.3726 -39134.7264 [116,] -36745.2130 -42304.3726 [117,] -32760.0580 -36745.2130 [118,] -41500.0659 -32760.0580 [119,] -42228.8996 -41500.0659 [120,] -31501.2936 -42228.8996 [121,] -39632.2950 -31501.2936 [122,] -35201.7884 -39632.2950 [123,] -28469.4748 -35201.7884 [124,] -40779.5001 -28469.4748 [125,] -41739.3726 -40779.5001 [126,] -42536.9338 -41739.3726 [127,] -123392.1292 -42536.9338 [128,] -124581.4028 -123392.1292 [129,] 196396.3015 -124581.4028 [130,] 311830.6417 196396.3015 [131,] 259655.8262 311830.6417 [132,] 28871.8957 259655.8262 [133,] 13620.5855 28871.8957 [134,] 4309.3256 13620.5855 [135,] 1218.4567 4309.3256 [136,] 48922.0226 1218.4567 [137,] 14283.7498 48922.0226 [138,] -8454.3222 14283.7498 [139,] -41446.6589 -8454.3222 [140,] -42429.9080 -41446.6589 [141,] -120671.4974 -42429.9080 [142,] -122830.7648 -120671.4974 [143,] -126825.1745 -122830.7648 [144,] -123360.0469 -126825.1745 [145,] 197047.9748 -123360.0469 [146,] 268738.5417 197047.9748 [147,] 248496.4379 268738.5417 [148,] 215141.6004 248496.4379 [149,] 249266.2120 215141.6004 [150,] 259885.3422 249266.2120 [151,] -2191.7581 259885.3422 [152,] 29995.1120 -2191.7581 [153,] -8534.1588 29995.1120 [154,] -42536.5657 -8534.1588 [155,] -122641.7626 -42536.5657 [156,] -127881.8844 -122641.7626 [157,] 197428.6301 -127881.8844 [158,] -8453.6121 197428.6301 [159,] -42471.9306 -8453.6121 [160,] -122192.6253 -42471.9306 [161,] -121559.3162 -122192.6253 [162,] -122821.8647 -121559.3162 [163,] 263338.9146 -122821.8647 [164,] 285389.5626 263338.9146 [165,] 260097.1803 285389.5626 [166,] -8452.8997 260097.1803 [167,] -38341.2725 -8452.8997 [168,] -42524.0163 -38341.2725 [169,] -122720.2949 -42524.0163 [170,] 197641.5980 -122720.2949 [171,] 256051.3986 197641.5980 [172,] 260878.9248 256051.3986 [173,] 237690.4336 260878.9248 [174,] -8455.3252 237690.4336 [175,] -38834.7399 -8455.3252 [176,] -38353.1194 -38834.7399 [177,] -43289.6112 -38353.1194 [178,] -40496.6762 -43289.6112 [179,] -38598.3549 -40496.6762 [180,] -40948.8268 -38598.3549 [181,] -40265.7478 -40948.8268 [182,] -39206.4434 -40265.7478 [183,] -39575.8156 -39206.4434 [184,] -42522.4556 -39575.8156 [185,] -124140.3273 -42522.4556 [186,] 197527.8752 -124140.3273 [187,] 256546.3170 197527.8752 [188,] 49463.3960 256546.3170 [189,] -8506.6096 49463.3960 [190,] -29968.6716 -8506.6096 [191,] -42541.6218 -29968.6716 [192,] 153729.5184 -42541.6218 [193,] 259503.8920 153729.5184 [194,] -8449.6286 259503.8920 [195,] -42495.7533 -8449.6286 [196,] -120984.3518 -42495.7533 [197,] -126412.2789 -120984.3518 [198,] -121028.5164 -126412.2789 [199,] -123941.4431 -121028.5164 [200,] 224608.4790 -123941.4431 [201,] 250925.1133 224608.4790 [202,] -21136.0257 250925.1133 [203,] -4027.6097 -21136.0257 [204,] -7375.9673 -4027.6097 [205,] -16117.7170 -7375.9673 [206,] -11910.1747 -16117.7170 [207,] 48743.8142 -11910.1747 [208,] -2556.0992 48743.8142 [209,] -4726.7520 -2556.0992 [210,] -25193.8165 -4726.7520 [211,] -22402.7120 -25193.8165 [212,] -20971.2314 -22402.7120 [213,] 38181.4139 -20971.2314 [214,] -4144.2725 38181.4139 [215,] 7633.2241 -4144.2725 [216,] -3526.7948 7633.2241 [217,] -1716.9207 -3526.7948 [218,] -4101.4031 -1716.9207 [219,] -8524.0974 -4101.4031 [220,] 51868.6960 -8524.0974 [221,] -3068.9476 51868.6960 [222,] -3239.3226 -3068.9476 [223,] -26256.1157 -3239.3226 [224,] -20684.4148 -26256.1157 [225,] 59919.3474 -20684.4148 [226,] -4242.1459 59919.3474 [227,] -8228.4725 -4242.1459 [228,] 48629.8409 -8228.4725 [229,] -25291.1916 48629.8409 [230,] 367.3699 -25291.1916 [231,] -4266.4539 367.3699 [232,] -6334.4799 -4266.4539 [233,] -7177.8073 -6334.4799 [234,] 48630.7644 -7177.8073 [235,] -25670.9695 48630.7644 [236,] -20794.1938 -25670.9695 [237,] -4084.8560 -20794.1938 [238,] -12260.4136 -4084.8560 [239,] 52178.4490 -12260.4136 [240,] 17391.2971 52178.4490 [241,] 17224.9782 17391.2971 [242,] -4530.1306 17224.9782 [243,] -26547.9644 -4530.1306 [244,] -24120.4174 -26547.9644 [245,] -20598.7764 -24120.4174 [246,] -4290.7220 -20598.7764 [247,] 48631.3500 -4290.7220 [248,] -20823.1303 48631.3500 [249,] -19456.8388 -20823.1303 [250,] -20670.3455 -19456.8388 [251,] -4284.5553 -20670.3455 [252,] 65593.1409 -4284.5553 [253,] -3258.0440 65593.1409 [254,] -22497.1280 -3258.0440 [255,] -20970.8928 -22497.1280 [256,] -4284.0191 -20970.8928 [257,] 48727.4158 -4284.0191 [258,] -3668.5562 48727.4158 [259,] -24644.5995 -3668.5562 [260,] -19750.2746 -24644.5995 [261,] -20635.5113 -19750.2746 [262,] -8371.7599 -20635.5113 [263,] -3935.6539 -8371.7599 [264,] -4287.9067 -3935.6539 [265,] 44678.0182 -4287.9067 [266,] -5203.2582 44678.0182 [267,] -22371.3936 -5203.2582 [268,] -20729.7748 -22371.3936 [269,] -1489.4416 -20729.7748 [270,] 3168.5435 -1489.4416 [271,] -4039.1749 3168.5435 [272,] -14916.0059 -4039.1749 [273,] 50776.8613 -14916.0059 [274,] 15841.2923 50776.8613 [275,] -4713.1975 15841.2923 [276,] -26647.9296 -4713.1975 [277,] -18647.1931 -26647.9296 [278,] -20686.2521 -18647.1931 [279,] -4291.3923 -20686.2521 [280,] 48630.4941 -4291.3923 [281,] -25394.1924 48630.4941 [282,] -23274.0781 -25394.1924 [283,] -21341.3626 -23274.0781 [284,] -3925.2582 -21341.3626 [285,] -6122.9133 -3925.2582 [286,] -1368.8576 -6122.9133 [287,] 51955.8116 -1368.8576 [288,] -9761.8961 51955.8116 [289,] -4491.3185 -9761.8961 [290,] -25253.6654 -4491.3185 [291,] 2342.6689 -25253.6654 [292,] -4242.9285 2342.6689 [293,] -8425.7215 -4242.9285 [294,] -1712.4997 -8425.7215 [295,] -6512.2831 -1712.4997 [296,] 9168.4821 -6512.2831 [297,] -14672.5101 9168.4821 [298,] 2602.5664 -14672.5101 [299,] 48631.6203 2602.5664 [300,] -25324.3927 48631.6203 [301,] -4297.4249 -25324.3927 [302,] 48630.9671 -4297.4249 [303,] -25479.1796 48630.9671 [304,] -22341.8481 -25479.1796 [305,] -4288.8452 -22341.8481 [306,] 48631.1473 -4288.8452 [307,] -25229.4741 48631.1473 [308,] -4285.3596 -25229.4741 [309,] 48304.0520 -4285.3596 [310,] -1067.1399 48304.0520 [311,] 8390.4218 -1067.1399 [312,] -723.2047 8390.4218 [313,] -5029.8283 -723.2047 [314,] -25304.2001 -5029.8283 [315,] -23026.9251 -25304.2001 [316,] -22423.4944 -23026.9251 [317,] -20712.3693 -22423.4944 [318,] -4297.6930 -20712.3693 [319,] 48717.7619 -4297.6930 [320,] -5004.1018 48717.7619 [321,] -26555.8428 -5004.1018 [322,] -20993.7571 -26555.8428 [323,] -2703.6314 -20993.7571 [324,] -9751.5130 -2703.6314 [325,] -9420.4433 -9751.5130 [326,] -9563.7042 -9420.4433 [327,] -7847.8881 -9563.7042 [328,] -2360.7844 -7847.8881 [329,] 48631.5752 -2360.7844 [330,] -22207.3850 48631.5752 [331,] -20545.4321 -22207.3850 [332,] -4288.5770 -20545.4321 [333,] 48634.0779 -4288.5770 [334,] 3967.6920 48634.0779 [335,] -25216.7085 3967.6920 [336,] -5316.6458 -25216.7085 [337,] -2439.4883 -5316.6458 [338,] -4283.7509 -2439.4883 [339,] 48631.0347 -4283.7509 [340,] -25218.2205 48631.0347 [341,] -3207.4540 -25218.2205 [342,] -3835.3488 -3207.4540 [343,] -9546.8239 -3835.3488 [344,] -8141.9104 -9546.8239 [345,] 4653.3535 -8141.9104 [346,] 48631.4401 4653.3535 [347,] -27802.1744 48631.4401 [348,] -20125.2750 -27802.1744 [349,] -20290.1680 -20125.2750 [350,] -20610.0740 -20290.1680 [351,] -3756.8343 -20610.0740 [352,] -13267.8453 -3756.8343 [353,] -14916.0444 -13267.8453 [354,] 50623.5737 -14916.0444 [355,] 3336.1802 50623.5737 [356,] -5267.2833 3336.1802 [357,] -3520.9037 -5267.2833 [358,] -3116.9002 -3520.9037 [359,] -154.3528 -3116.9002 [360,] -2081.9245 -154.3528 [361,] -4565.2495 -2081.9245 [362,] -2655.7431 -4565.2495 [363,] 4407.2953 -2655.7431 [364,] 21156.6960 4407.2953 [365,] -5919.7833 21156.6960 [366,] -11079.9913 -5919.7833 [367,] 56365.9726 -11079.9913 [368,] -3716.4247 56365.9726 [369,] -42610.4884 -3716.4247 [370,] 5579.7210 -42610.4884 [371,] 2280.9130 5579.7210 [372,] -4339.0472 2280.9130 [373,] 5913.3227 -4339.0472 [374,] -774.2213 5913.3227 [375,] 7283.5463 -774.2213 [376,] 3586.8921 7283.5463 [377,] -5772.2392 3586.8921 [378,] -5883.6893 -5772.2392 [379,] 5795.9656 -5883.6893 [380,] 7069.4099 5795.9656 [381,] -1580.4930 7069.4099 [382,] -7797.1841 -1580.4930 [383,] 10873.1710 -7797.1841 [384,] 12943.1023 10873.1710 [385,] 9120.8558 12943.1023 [386,] -28657.4555 9120.8558 [387,] -4538.7319 -28657.4555 [388,] 4504.5095 -4538.7319 [389,] 887.1391 4504.5095 [390,] 7664.3256 887.1391 [391,] 12907.1600 7664.3256 [392,] -78049.3679 12907.1600 [393,] -4420.4893 -78049.3679 [394,] -7286.6132 -4420.4893 [395,] 2671.2908 -7286.6132 [396,] -12977.3502 2671.2908 [397,] 98828.5992 -12977.3502 [398,] 14697.6847 98828.5992 [399,] 14324.9978 14697.6847 [400,] -37005.7122 14324.9978 [401,] 173260.5421 -37005.7122 [402,] -58478.9410 173260.5421 [403,] 46555.8013 -58478.9410 [404,] -132662.9390 46555.8013 [405,] 70131.4619 -132662.9390 [406,] 80543.8928 70131.4619 [407,] 51361.6742 80543.8928 [408,] 159512.3570 51361.6742 [409,] -23264.7612 159512.3570 [410,] -25458.0482 -23264.7612 [411,] 168531.3781 -25458.0482 [412,] -25979.2891 168531.3781 [413,] -73759.0312 -25979.2891 [414,] -42824.1270 -73759.0312 [415,] 33945.0739 -42824.1270 [416,] 77301.0111 33945.0739 [417,] 212851.0134 77301.0111 [418,] 47210.9221 212851.0134 [419,] -30969.8114 47210.9221 [420,] 49687.9245 -30969.8114 [421,] 16191.7387 49687.9245 [422,] 9682.4484 16191.7387 [423,] 96476.0414 9682.4484 [424,] 47629.9444 96476.0414 [425,] 68299.9656 47629.9444 [426,] 183.6247 68299.9656 [427,] 9706.6965 183.6247 [428,] -38270.7328 9706.6965 [429,] -31333.3773 -38270.7328 [430,] -18298.5227 -31333.3773 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 81350.2769 175839.7180 2 78590.8368 81350.2769 3 -58681.6234 78590.8368 4 -485.3187 -58681.6234 5 9302.4862 -485.3187 6 -2342.5076 9302.4862 7 30317.0338 -2342.5076 8 -25836.0535 30317.0338 9 -16436.4824 -25836.0535 10 12935.2203 -16436.4824 11 -27351.5918 12935.2203 12 -28865.8114 -27351.5918 13 -49714.4406 -28865.8114 14 -21479.9857 -49714.4406 15 -6502.1158 -21479.9857 16 10108.2758 -6502.1158 17 -19474.5117 10108.2758 18 -31303.1213 -19474.5117 19 -27201.2617 -31303.1213 20 -25048.1118 -27201.2617 21 -29026.3780 -25048.1118 22 -12606.5837 -29026.3780 23 -30624.4492 -12606.5837 24 -28358.5798 -30624.4492 25 -35041.5445 -28358.5798 26 -32093.9150 -35041.5445 27 -50629.9749 -32093.9150 28 -35299.0443 -50629.9749 29 -41863.8545 -35299.0443 30 -26516.4798 -41863.8545 31 -27325.0288 -26516.4798 32 -27426.6463 -27325.0288 33 -26987.8611 -27426.6463 34 14509.6491 -26987.8611 35 -33089.4648 14509.6491 36 -29233.3932 -33089.4648 37 -23046.5555 -29233.3932 38 -11125.4232 -23046.5555 39 -41198.4622 -11125.4232 40 -29868.7112 -41198.4622 41 -34730.4736 -29868.7112 42 -31775.0815 -34730.4736 43 -29878.6655 -31775.0815 44 -38898.3383 -29878.6655 45 -32225.1461 -38898.3383 46 -12126.1734 -32225.1461 47 -25154.4592 -12126.1734 48 -26807.4039 -25154.4592 49 -42183.6733 -26807.4039 50 -28103.6677 -42183.6733 51 -23348.5010 -28103.6677 52 -13432.7929 -23348.5010 53 -20883.6733 -13432.7929 54 -25724.5201 -20883.6733 55 -40121.7901 -25724.5201 56 -30242.1033 -40121.7901 57 -33923.0665 -30242.1033 58 -34249.6826 -33923.0665 59 -35620.8912 -34249.6826 60 -31242.9225 -35620.8912 61 -46577.2192 -31242.9225 62 -32908.5839 -46577.2192 63 -37088.1477 -32908.5839 64 -41755.8466 -37088.1477 65 -36088.6335 -41755.8466 66 -39710.3953 -36088.6335 67 -32516.0638 -39710.3953 68 -27292.1925 -32516.0638 69 -31978.9497 -27292.1925 70 -51577.6128 -31978.9497 71 -44931.0421 -51577.6128 72 -33211.3911 -44931.0421 73 -37147.3803 -33211.3911 74 -28738.6549 -37147.3803 75 -40452.7374 -28738.6549 76 -30380.3899 -40452.7374 77 -41576.4536 -30380.3899 78 -40775.5548 -41576.4536 79 -30794.3775 -40775.5548 80 -41027.9823 -30794.3775 81 -34013.8013 -41027.9823 82 -26149.1762 -34013.8013 83 -46150.1871 -26149.1762 84 -31374.2837 -46150.1871 85 -40724.0988 -31374.2837 86 -33255.9965 -40724.0988 87 -31772.0565 -33255.9965 88 -33548.7993 -31772.0565 89 -30381.1736 -33548.7993 90 -33561.6713 -30381.1736 91 -36284.7098 -33561.6713 92 -38616.2274 -36284.7098 93 -28389.5734 -38616.2274 94 -45647.8919 -28389.5734 95 -40769.6802 -45647.8919 96 -25919.1456 -40769.6802 97 -40334.0261 -25919.1456 98 -25876.0592 -40334.0261 99 -38638.1751 -25876.0592 100 -39277.6458 -38638.1751 101 -38323.7626 -39277.6458 102 -33486.7093 -38323.7626 103 -31380.0197 -33486.7093 104 -39281.1150 -31380.0197 105 -41029.5240 -39281.1150 106 -47182.4188 -41029.5240 107 -28877.8311 -47182.4188 108 -34309.0997 -28877.8311 109 -41450.0565 -34309.0997 110 -38122.4971 -41450.0565 111 -48685.6801 -38122.4971 112 -38012.9905 -48685.6801 113 -39687.8624 -38012.9905 114 -39134.7264 -39687.8624 115 -42304.3726 -39134.7264 116 -36745.2130 -42304.3726 117 -32760.0580 -36745.2130 118 -41500.0659 -32760.0580 119 -42228.8996 -41500.0659 120 -31501.2936 -42228.8996 121 -39632.2950 -31501.2936 122 -35201.7884 -39632.2950 123 -28469.4748 -35201.7884 124 -40779.5001 -28469.4748 125 -41739.3726 -40779.5001 126 -42536.9338 -41739.3726 127 -123392.1292 -42536.9338 128 -124581.4028 -123392.1292 129 196396.3015 -124581.4028 130 311830.6417 196396.3015 131 259655.8262 311830.6417 132 28871.8957 259655.8262 133 13620.5855 28871.8957 134 4309.3256 13620.5855 135 1218.4567 4309.3256 136 48922.0226 1218.4567 137 14283.7498 48922.0226 138 -8454.3222 14283.7498 139 -41446.6589 -8454.3222 140 -42429.9080 -41446.6589 141 -120671.4974 -42429.9080 142 -122830.7648 -120671.4974 143 -126825.1745 -122830.7648 144 -123360.0469 -126825.1745 145 197047.9748 -123360.0469 146 268738.5417 197047.9748 147 248496.4379 268738.5417 148 215141.6004 248496.4379 149 249266.2120 215141.6004 150 259885.3422 249266.2120 151 -2191.7581 259885.3422 152 29995.1120 -2191.7581 153 -8534.1588 29995.1120 154 -42536.5657 -8534.1588 155 -122641.7626 -42536.5657 156 -127881.8844 -122641.7626 157 197428.6301 -127881.8844 158 -8453.6121 197428.6301 159 -42471.9306 -8453.6121 160 -122192.6253 -42471.9306 161 -121559.3162 -122192.6253 162 -122821.8647 -121559.3162 163 263338.9146 -122821.8647 164 285389.5626 263338.9146 165 260097.1803 285389.5626 166 -8452.8997 260097.1803 167 -38341.2725 -8452.8997 168 -42524.0163 -38341.2725 169 -122720.2949 -42524.0163 170 197641.5980 -122720.2949 171 256051.3986 197641.5980 172 260878.9248 256051.3986 173 237690.4336 260878.9248 174 -8455.3252 237690.4336 175 -38834.7399 -8455.3252 176 -38353.1194 -38834.7399 177 -43289.6112 -38353.1194 178 -40496.6762 -43289.6112 179 -38598.3549 -40496.6762 180 -40948.8268 -38598.3549 181 -40265.7478 -40948.8268 182 -39206.4434 -40265.7478 183 -39575.8156 -39206.4434 184 -42522.4556 -39575.8156 185 -124140.3273 -42522.4556 186 197527.8752 -124140.3273 187 256546.3170 197527.8752 188 49463.3960 256546.3170 189 -8506.6096 49463.3960 190 -29968.6716 -8506.6096 191 -42541.6218 -29968.6716 192 153729.5184 -42541.6218 193 259503.8920 153729.5184 194 -8449.6286 259503.8920 195 -42495.7533 -8449.6286 196 -120984.3518 -42495.7533 197 -126412.2789 -120984.3518 198 -121028.5164 -126412.2789 199 -123941.4431 -121028.5164 200 224608.4790 -123941.4431 201 250925.1133 224608.4790 202 -21136.0257 250925.1133 203 -4027.6097 -21136.0257 204 -7375.9673 -4027.6097 205 -16117.7170 -7375.9673 206 -11910.1747 -16117.7170 207 48743.8142 -11910.1747 208 -2556.0992 48743.8142 209 -4726.7520 -2556.0992 210 -25193.8165 -4726.7520 211 -22402.7120 -25193.8165 212 -20971.2314 -22402.7120 213 38181.4139 -20971.2314 214 -4144.2725 38181.4139 215 7633.2241 -4144.2725 216 -3526.7948 7633.2241 217 -1716.9207 -3526.7948 218 -4101.4031 -1716.9207 219 -8524.0974 -4101.4031 220 51868.6960 -8524.0974 221 -3068.9476 51868.6960 222 -3239.3226 -3068.9476 223 -26256.1157 -3239.3226 224 -20684.4148 -26256.1157 225 59919.3474 -20684.4148 226 -4242.1459 59919.3474 227 -8228.4725 -4242.1459 228 48629.8409 -8228.4725 229 -25291.1916 48629.8409 230 367.3699 -25291.1916 231 -4266.4539 367.3699 232 -6334.4799 -4266.4539 233 -7177.8073 -6334.4799 234 48630.7644 -7177.8073 235 -25670.9695 48630.7644 236 -20794.1938 -25670.9695 237 -4084.8560 -20794.1938 238 -12260.4136 -4084.8560 239 52178.4490 -12260.4136 240 17391.2971 52178.4490 241 17224.9782 17391.2971 242 -4530.1306 17224.9782 243 -26547.9644 -4530.1306 244 -24120.4174 -26547.9644 245 -20598.7764 -24120.4174 246 -4290.7220 -20598.7764 247 48631.3500 -4290.7220 248 -20823.1303 48631.3500 249 -19456.8388 -20823.1303 250 -20670.3455 -19456.8388 251 -4284.5553 -20670.3455 252 65593.1409 -4284.5553 253 -3258.0440 65593.1409 254 -22497.1280 -3258.0440 255 -20970.8928 -22497.1280 256 -4284.0191 -20970.8928 257 48727.4158 -4284.0191 258 -3668.5562 48727.4158 259 -24644.5995 -3668.5562 260 -19750.2746 -24644.5995 261 -20635.5113 -19750.2746 262 -8371.7599 -20635.5113 263 -3935.6539 -8371.7599 264 -4287.9067 -3935.6539 265 44678.0182 -4287.9067 266 -5203.2582 44678.0182 267 -22371.3936 -5203.2582 268 -20729.7748 -22371.3936 269 -1489.4416 -20729.7748 270 3168.5435 -1489.4416 271 -4039.1749 3168.5435 272 -14916.0059 -4039.1749 273 50776.8613 -14916.0059 274 15841.2923 50776.8613 275 -4713.1975 15841.2923 276 -26647.9296 -4713.1975 277 -18647.1931 -26647.9296 278 -20686.2521 -18647.1931 279 -4291.3923 -20686.2521 280 48630.4941 -4291.3923 281 -25394.1924 48630.4941 282 -23274.0781 -25394.1924 283 -21341.3626 -23274.0781 284 -3925.2582 -21341.3626 285 -6122.9133 -3925.2582 286 -1368.8576 -6122.9133 287 51955.8116 -1368.8576 288 -9761.8961 51955.8116 289 -4491.3185 -9761.8961 290 -25253.6654 -4491.3185 291 2342.6689 -25253.6654 292 -4242.9285 2342.6689 293 -8425.7215 -4242.9285 294 -1712.4997 -8425.7215 295 -6512.2831 -1712.4997 296 9168.4821 -6512.2831 297 -14672.5101 9168.4821 298 2602.5664 -14672.5101 299 48631.6203 2602.5664 300 -25324.3927 48631.6203 301 -4297.4249 -25324.3927 302 48630.9671 -4297.4249 303 -25479.1796 48630.9671 304 -22341.8481 -25479.1796 305 -4288.8452 -22341.8481 306 48631.1473 -4288.8452 307 -25229.4741 48631.1473 308 -4285.3596 -25229.4741 309 48304.0520 -4285.3596 310 -1067.1399 48304.0520 311 8390.4218 -1067.1399 312 -723.2047 8390.4218 313 -5029.8283 -723.2047 314 -25304.2001 -5029.8283 315 -23026.9251 -25304.2001 316 -22423.4944 -23026.9251 317 -20712.3693 -22423.4944 318 -4297.6930 -20712.3693 319 48717.7619 -4297.6930 320 -5004.1018 48717.7619 321 -26555.8428 -5004.1018 322 -20993.7571 -26555.8428 323 -2703.6314 -20993.7571 324 -9751.5130 -2703.6314 325 -9420.4433 -9751.5130 326 -9563.7042 -9420.4433 327 -7847.8881 -9563.7042 328 -2360.7844 -7847.8881 329 48631.5752 -2360.7844 330 -22207.3850 48631.5752 331 -20545.4321 -22207.3850 332 -4288.5770 -20545.4321 333 48634.0779 -4288.5770 334 3967.6920 48634.0779 335 -25216.7085 3967.6920 336 -5316.6458 -25216.7085 337 -2439.4883 -5316.6458 338 -4283.7509 -2439.4883 339 48631.0347 -4283.7509 340 -25218.2205 48631.0347 341 -3207.4540 -25218.2205 342 -3835.3488 -3207.4540 343 -9546.8239 -3835.3488 344 -8141.9104 -9546.8239 345 4653.3535 -8141.9104 346 48631.4401 4653.3535 347 -27802.1744 48631.4401 348 -20125.2750 -27802.1744 349 -20290.1680 -20125.2750 350 -20610.0740 -20290.1680 351 -3756.8343 -20610.0740 352 -13267.8453 -3756.8343 353 -14916.0444 -13267.8453 354 50623.5737 -14916.0444 355 3336.1802 50623.5737 356 -5267.2833 3336.1802 357 -3520.9037 -5267.2833 358 -3116.9002 -3520.9037 359 -154.3528 -3116.9002 360 -2081.9245 -154.3528 361 -4565.2495 -2081.9245 362 -2655.7431 -4565.2495 363 4407.2953 -2655.7431 364 21156.6960 4407.2953 365 -5919.7833 21156.6960 366 -11079.9913 -5919.7833 367 56365.9726 -11079.9913 368 -3716.4247 56365.9726 369 -42610.4884 -3716.4247 370 5579.7210 -42610.4884 371 2280.9130 5579.7210 372 -4339.0472 2280.9130 373 5913.3227 -4339.0472 374 -774.2213 5913.3227 375 7283.5463 -774.2213 376 3586.8921 7283.5463 377 -5772.2392 3586.8921 378 -5883.6893 -5772.2392 379 5795.9656 -5883.6893 380 7069.4099 5795.9656 381 -1580.4930 7069.4099 382 -7797.1841 -1580.4930 383 10873.1710 -7797.1841 384 12943.1023 10873.1710 385 9120.8558 12943.1023 386 -28657.4555 9120.8558 387 -4538.7319 -28657.4555 388 4504.5095 -4538.7319 389 887.1391 4504.5095 390 7664.3256 887.1391 391 12907.1600 7664.3256 392 -78049.3679 12907.1600 393 -4420.4893 -78049.3679 394 -7286.6132 -4420.4893 395 2671.2908 -7286.6132 396 -12977.3502 2671.2908 397 98828.5992 -12977.3502 398 14697.6847 98828.5992 399 14324.9978 14697.6847 400 -37005.7122 14324.9978 401 173260.5421 -37005.7122 402 -58478.9410 173260.5421 403 46555.8013 -58478.9410 404 -132662.9390 46555.8013 405 70131.4619 -132662.9390 406 80543.8928 70131.4619 407 51361.6742 80543.8928 408 159512.3570 51361.6742 409 -23264.7612 159512.3570 410 -25458.0482 -23264.7612 411 168531.3781 -25458.0482 412 -25979.2891 168531.3781 413 -73759.0312 -25979.2891 414 -42824.1270 -73759.0312 415 33945.0739 -42824.1270 416 77301.0111 33945.0739 417 212851.0134 77301.0111 418 47210.9221 212851.0134 419 -30969.8114 47210.9221 420 49687.9245 -30969.8114 421 16191.7387 49687.9245 422 9682.4484 16191.7387 423 96476.0414 9682.4484 424 47629.9444 96476.0414 425 68299.9656 47629.9444 426 183.6247 68299.9656 427 9706.6965 183.6247 428 -38270.7328 9706.6965 429 -31333.3773 -38270.7328 430 -18298.5227 -31333.3773 > 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/7qgg31291299426.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/8jqxo1291299426.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/9jqxo1291299426.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/10jqxo1291299426.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/11fzvx1291299426.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/1200u31291299426.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/13essc1291299426.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/14ia8z1291299426.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/153bo51291299426.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/166b5t1291299426.tab") + } > > try(system("convert tmp/1nyif1291299426.ps tmp/1nyif1291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/2nyif1291299426.ps tmp/2nyif1291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/3nyif1291299426.ps tmp/3nyif1291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/4f7z01291299426.ps tmp/4f7z01291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/5f7z01291299426.ps tmp/5f7z01291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/6qgg31291299426.ps tmp/6qgg31291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/7qgg31291299426.ps tmp/7qgg31291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/8jqxo1291299426.ps tmp/8jqxo1291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/9jqxo1291299426.ps tmp/9jqxo1291299426.png",intern=TRUE)) character(0) > try(system("convert tmp/10jqxo1291299426.ps tmp/10jqxo1291299426.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 16.425 3.184 16.918