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(1 + ,162556 + ,1081 + ,213118 + ,6282154 + ,1 + ,29790 + ,309 + ,81767 + ,4321023 + ,1 + ,87550 + ,458 + ,153198 + ,4111912 + ,0 + ,84738 + ,588 + ,-26007 + ,223193 + ,1 + ,54660 + ,302 + ,126942 + ,1491348 + ,1 + ,42634 + ,156 + ,157214 + ,1629616 + ,0 + ,40949 + ,481 + ,129352 + ,1398893 + ,1 + ,45187 + ,353 + ,234817 + ,1926517 + ,1 + ,37704 + ,452 + ,60448 + ,983660 + ,1 + ,16275 + ,109 + ,47818 + ,1443586 + ,0 + ,25830 + ,115 + ,245546 + ,1073089 + ,0 + ,12679 + ,110 + ,48020 + ,984885 + ,1 + ,18014 + ,239 + ,-1710 + ,1405225 + ,0 + ,43556 + ,247 + ,32648 + ,227132 + ,1 + ,24811 + ,505 + ,95350 + ,929118 + ,0 + ,6575 + ,159 + ,151352 + ,1071292 + ,0 + ,7123 + ,109 + ,288170 + ,638830 + ,1 + ,21950 + ,519 + ,114337 + ,856956 + ,1 + ,37597 + ,248 + ,37884 + ,992426 + ,0 + ,17821 + ,373 + ,122844 + ,444477 + ,1 + ,12988 + ,119 + ,82340 + ,857217 + ,1 + ,22330 + ,84 + ,79801 + ,711969 + ,0 + ,13326 + ,102 + ,165548 + ,702380 + ,0 + ,16189 + ,295 + ,116384 + ,358589 + ,0 + ,7146 + ,105 + ,134028 + ,297978 + ,0 + ,15824 + ,64 + ,63838 + ,585715 + ,1 + ,27664 + ,282 + ,74996 + ,657954 + ,0 + ,11920 + ,182 + ,31080 + ,209458 + ,0 + ,8568 + ,37 + ,32168 + ,786690 + ,0 + ,14416 + ,361 + ,49857 + ,439798 + ,1 + ,3369 + ,28 + ,87161 + ,688779 + ,1 + ,11819 + ,85 + ,106113 + ,574339 + ,1 + ,6984 + ,45 + ,80570 + ,741409 + ,1 + ,4519 + ,49 + ,102129 + ,597793 + ,0 + ,2220 + ,22 + ,301670 + ,644190 + ,0 + ,18562 + ,155 + ,102313 + ,377934 + ,0 + ,10327 + ,91 + ,88577 + ,640273 + ,1 + ,5336 + ,81 + ,112477 + ,697458 + ,1 + ,2365 + ,79 + ,191778 + ,550608 + ,0 + ,4069 + ,145 + ,79804 + ,207393 + ,0 + ,8636 + ,855 + ,128294 + ,301607 + ,0 + ,13718 + ,61 + ,96448 + ,345783 + ,0 + ,4525 + ,226 + ,93811 + ,501749 + ,0 + ,6869 + ,105 + ,117520 + ,379983 + ,0 + ,4628 + ,62 + ,69159 + ,387475 + ,1 + ,3689 + ,25 + ,101792 + ,377305 + ,1 + ,4891 + ,217 + ,210568 + ,370837 + ,1 + ,7489 + ,322 + ,136996 + ,430866 + ,0 + ,4901 + ,84 + ,121920 + ,469107 + ,0 + ,2284 + ,33 + ,76403 + ,194493 + ,1 + ,3160 + ,108 + ,108094 + ,530670 + ,1 + ,4150 + ,150 + ,134759 + ,518365 + ,1 + ,7285 + ,115 + ,188873 + ,491303 + ,1 + ,1134 + ,162 + ,146216 + ,527021 + ,1 + ,4658 + ,158 + ,156608 + ,233773 + ,0 + ,2384 + ,97 + ,61348 + ,405972 + ,0 + ,3748 + ,9 + ,50350 + ,652925 + ,0 + ,5371 + ,66 + ,87720 + ,446211 + ,0 + ,1285 + ,107 + ,99489 + ,341340 + ,1 + ,9327 + ,101 + ,87419 + ,387699 + ,1 + ,5565 + ,47 + ,94355 + ,493408 + ,0 + ,1528 + ,38 + ,60326 + ,146494 + ,1 + ,3122 + ,34 + ,94670 + ,414462 + ,1 + ,7561 + ,87 + ,82425 + ,364304 + ,0 + ,2675 + ,79 + ,59017 + ,355178 + ,0 + ,13253 + ,947 + ,90829 + ,357760 + ,0 + ,880 + ,74 + ,80791 + ,261216 + ,1 + ,2053 + ,53 + ,100423 + ,397144 + ,0 + ,1424 + ,94 + ,131116 + ,374943 + ,1 + ,4036 + ,63 + ,100269 + ,424898 + ,1 + ,3045 + ,58 + ,27330 + ,202055 + ,0 + ,5119 + ,49 + ,39039 + ,378525 + ,0 + ,1431 + ,34 + ,106885 + ,310768 + ,0 + ,554 + ,11 + ,79285 + ,325738 + ,0 + ,1975 + ,35 + ,118881 + ,394510 + ,1 + ,1765 + ,20 + ,77623 + ,247060 + ,0 + ,1012 + ,47 + ,114768 + ,368078 + ,0 + ,810 + ,43 + ,74015 + ,236761 + ,0 + ,1280 + ,117 + ,69465 + ,312378 + ,1 + ,666 + ,171 + ,117869 + ,339836 + ,0 + ,1380 + ,26 + ,60982 + ,347385 + ,1 + ,4677 + ,75 + ,90131 + ,426280 + ,0 + ,876 + ,59 + ,138971 + ,352850 + ,0 + ,814 + ,18 + ,39625 + ,301881 + ,0 + ,514 + ,15 + ,102725 + ,377516 + ,1 + ,5692 + ,72 + ,64239 + ,357312 + ,0 + ,3642 + ,86 + ,90262 + ,458343 + ,0 + ,540 + ,14 + ,103960 + ,354228 + ,0 + ,2099 + ,64 + ,106611 + ,308636 + ,0 + ,567 + ,11 + ,103345 + ,386212 + ,0 + ,2001 + ,52 + ,95551 + ,393343 + ,1 + ,2949 + ,41 + ,82903 + ,378509 + ,0 + ,2253 + ,99 + ,63593 + ,452469 + ,1 + ,6533 + ,75 + ,126910 + ,364839 + ,0 + ,1889 + ,45 + ,37527 + ,358649 + ,1 + ,3055 + ,43 + ,60247 + ,376641 + ,0 + ,272 + ,8 + ,112995 + ,429112 + ,1 + ,1414 + ,198 + ,70184 + ,330546 + ,0 + ,2564 + ,22 + ,130140 + ,403560 + ,1 + ,1383 + ,11 + ,73221 + ,317892 + ,1 + ,1261 + ,33 + ,76114 + ,307528 + ,0 + ,975 + ,23 + ,90534 + ,235133 + ,0 + ,3366 + ,80 + ,108479 + ,299243 + ,0 + ,576 + ,18 + ,113761 + ,314073 + ,0 + ,1686 + ,40 + ,68696 + ,368186 + ,0 + ,746 + ,23 + ,71561 + ,269661 + ,1 + ,3192 + ,60 + ,59831 + ,125390 + ,1 + ,2045 + ,20 + ,97890 + ,510834 + ,0 + ,5702 + ,61 + ,101481 + ,321896 + ,1 + ,1932 + ,36 + ,72954 + ,249898 + ,0 + ,936 + ,30 + ,67939 + ,408881 + ,1 + ,3437 + ,47 + ,48022 + ,158492 + ,0 + ,5131 + ,71 + ,86111 + ,292154 + ,1 + ,2397 + ,14 + ,74020 + ,289513 + ,1 + ,1389 + ,9 + ,57530 + ,378049 + ,0 + ,1503 + ,39 + ,56364 + ,343466 + ,0 + ,402 + ,26 + ,84990 + ,332743 + ,0 + ,2239 + ,21 + ,88590 + ,442882 + ,1 + ,2234 + ,16 + ,77200 + ,214215 + ,0 + ,837 + ,69 + ,61262 + ,315688 + ,0 + ,10579 + ,92 + ,110309 + ,375195 + ,0 + ,875 + ,14 + ,67000 + ,334280 + ,0 + ,1585 + ,107 + ,93099 + ,355864 + ,1 + ,1659 + ,29 + ,107577 + ,480382 + ,1 + ,2647 + ,37 + ,62920 + ,353058 + ,1 + ,3294 + ,23 + ,75832 + ,217193 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,94 + ,7 + ,60793 + ,314533 + ,0 + ,422 + ,28 + ,57935 + ,318056 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,34 + ,8 + ,60630 + ,314353 + ,0 + ,1558 + ,63 + ,55637 + ,369448 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,43 + ,3 + ,60887 + ,312846 + ,0 + ,645 + ,5 + ,60720 + ,312075 + ,0 + ,316 + ,9 + ,60505 + ,315009 + ,0 + ,115 + ,13 + ,60945 + ,318903 + ,0 + ,5 + ,2 + ,60720 + ,314887 + ,0 + ,897 + ,5 + ,60720 + ,314913 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,389 + ,14 + ,58990 + ,325506 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1002 + ,15 + ,56750 + ,298568 + ,0 + ,36 + ,3 + ,60894 + ,315834 + ,0 + ,460 + ,15 + ,63346 + ,329784 + ,0 + ,309 + ,11 + ,56535 + ,312878 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,9 + ,6 + ,60835 + ,314987 + ,0 + ,271 + ,2 + ,60720 + ,325249 + ,0 + ,14 + ,1 + ,61016 + ,315877 + ,0 + ,520 + ,10 + ,58650 + ,291650 + ,0 + ,1766 + ,73 + ,60438 + ,305959 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,458 + ,11 + ,58625 + ,297765 + ,0 + ,20 + ,3 + ,60938 + ,315245 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,98 + ,2 + ,61490 + ,315236 + ,0 + ,405 + ,7 + ,60845 + ,336425 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,483 + ,27 + ,60830 + ,306268 + ,0 + ,454 + ,51 + ,63261 + ,302187 + ,0 + ,47 + ,3 + ,60720 + ,314882 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,757 + ,19 + ,45689 + ,382712 + ,0 + ,4655 + ,393 + ,60720 + ,341570 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,36 + ,4 + ,61564 + ,312412 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,126 + ,10 + ,60951 + ,315547 + ,0 + ,400 + ,152 + ,60720 + ,313267 + ,1 + ,71 + ,1 + ,60745 + ,316176 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,972 + ,34 + ,71642 + ,359335 + ,1 + ,531 + ,10 + ,71641 + ,330068 + ,0 + ,2461 + ,57 + ,55792 + ,314289 + ,1 + ,378 + ,52 + ,71873 + ,297413 + ,1 + ,23 + ,5 + ,62555 + ,314806 + ,1 + ,638 + ,14 + ,60370 + ,333210 + ,1 + ,2300 + ,29 + ,64873 + ,352108 + ,0 + ,149 + ,5 + ,62041 + ,313332 + ,0 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,275 + ,4 + ,59500 + ,318745 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,141 + ,6 + ,61630 + ,315366 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,28 + ,2 + ,60890 + ,315688 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4980 + ,91 + ,113521 + ,409642 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,472 + ,20 + ,80045 + ,269587 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,203 + ,27 + ,50804 + ,300962 + ,1 + ,496 + ,17 + ,87390 + ,325479 + ,1 + ,10 + ,2 + ,61656 + ,316155 + ,1 + ,63 + ,4 + ,65688 + ,318574 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1136 + ,32 + ,48522 + ,343613 + ,1 + ,265 + ,31 + ,60720 + ,306948 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,267 + ,32 + ,57640 + ,330059 + ,0 + ,474 + ,20 + ,61977 + ,288985 + ,0 + ,534 + ,7 + ,62620 + ,304485 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,15 + ,8 + ,60831 + ,315688 + ,0 + ,397 + ,28 + ,60646 + ,317736 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1866 + ,29 + ,56225 + ,322331 + ,0 + ,288 + ,4 + ,60510 + ,296656 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,3 + ,2 + ,60698 + ,315354 + ,0 + ,468 + ,21 + ,60720 + ,312161 + ,0 + ,20 + ,2 + ,60805 + ,315576 + ,0 + ,278 + ,26 + ,61404 + ,314922 + ,0 + ,61 + ,14 + ,60720 + ,314551 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,192 + ,4 + ,65276 + ,312339 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,317 + ,9 + ,63915 + ,298700 + ,0 + ,738 + ,10 + ,60720 + ,321376 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,368 + ,17 + ,61686 + ,303230 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2 + ,1 + ,60743 + ,315487 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,53 + ,6 + ,60349 + ,315793 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,94 + ,3 + ,61360 + ,312887 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,24 + ,8 + ,59818 + ,315637 + ,0 + ,2332 + ,4 + ,72680 + ,324385 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,131 + ,11 + ,61808 + ,308989 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,206 + ,9 + ,53110 + ,296702 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,167 + ,2 + ,64245 + ,307322 + ,0 + ,622 + ,73 + ,73007 + ,304376 + ,0 + ,2328 + ,94 + ,82732 + ,253588 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,8 + ,54820 + ,309560 + ,0 + ,364 + ,35 + ,47705 + ,298466 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,226 + ,12 + ,72835 + ,343929 + ,0 + ,307 + ,15 + ,58856 + ,331955 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,188 + ,11 + ,77655 + ,381180 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,138 + ,6 + ,69817 + ,331420 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,125 + ,12 + ,60798 + ,310201 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,282 + ,30 + ,62452 + ,320016 + ,0 + ,335 + ,33 + ,64175 + ,320398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1324 + ,117 + ,67440 + ,291841 + ,0 + ,176 + ,28 + ,68136 + ,310670 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,249 + ,72 + ,56726 + ,313491 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,333 + ,13 + ,70811 + ,331323 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,601 + ,6 + ,60720 + ,319210 + ,0 + ,30 + ,4 + ,62045 + ,318098 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,249 + ,62 + ,54323 + ,292754 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,165 + ,24 + ,62841 + ,325176 + ,0 + ,453 + ,21 + ,81125 + ,365959 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,53 + ,14 + ,59506 + ,302409 + ,1 + ,382 + ,21 + ,59365 + ,340968 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,30 + ,4 + ,60798 + ,313164 + ,0 + ,290 + ,2 + ,58790 + ,301164 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,366 + ,53 + ,61808 + ,344425 + ,0 + ,2 + ,9 + ,60735 + ,315394 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,209 + ,13 + ,64016 + ,316647 + ,0 + ,384 + ,22 + ,54683 + ,309836 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,83 + ,87192 + ,346611 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,49 + ,8 + ,64107 + ,322031 + ,0 + ,3 + ,4 + ,60761 + ,315656 + ,0 + ,133 + ,14 + ,65990 + ,339445 + ,0 + ,32 + ,1 + ,59988 + ,314964 + ,0 + ,368 + ,17 + ,61167 + ,297141 + ,0 + ,1 + ,6 + ,60719 + ,315372 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,22 + ,2 + ,60722 + ,312502 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,96 + ,5 + ,60379 + ,313729 + ,0 + ,1 + ,2 + ,60727 + ,315388 + ,0 + ,314 + ,5 + ,60720 + ,315371 + ,0 + ,844 + ,78 + ,60925 + ,296139 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,26 + ,1 + ,60896 + ,313880 + ,0 + ,125 + ,13 + ,59734 + ,317698 + ,0 + ,304 + ,15 + ,62969 + ,295580 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,621 + ,48 + ,60720 + ,308256 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,119 + ,6 + ,59118 + ,303677 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1595 + ,17 + ,60720 + ,319369 + ,0 + ,312 + ,14 + ,58598 + ,318690 + ,0 + ,60 + ,10 + ,61124 + ,314049 + ,0 + ,587 + ,12 + ,59595 + ,325699 + ,0 + ,135 + ,2 + ,62065 + ,314210 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,514 + ,52 + ,78780 + ,322378 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1763 + ,24 + ,61600 + ,308336 + ,0 + ,180 + ,11 + ,59635 + ,316386 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,218 + ,21 + ,60720 + ,315553 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,448 + ,40 + ,59781 + ,323361 + ,0 + ,227 + ,9 + ,76644 + ,336639 + ,0 + ,174 + ,1 + ,64820 + ,307424 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,121 + ,24 + ,56178 + ,295370 + ,0 + ,607 + ,11 + ,60436 + ,322340 + ,0 + ,2212 + ,14 + ,60720 + ,319864 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,530 + ,60 + ,73433 + ,317291 + ,0 + ,571 + ,80 + ,41477 + ,280398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,78 + ,16 + ,62700 + ,317330 + ,0 + ,2489 + ,40 + ,67804 + ,238125 + ,0 + ,131 + ,6 + ,59661 + ,327071 + ,0 + ,923 + ,8 + ,58620 + ,309038 + ,0 + ,72 + ,3 + ,60398 + ,314210 + ,0 + ,572 + ,16 + ,58580 + ,307930 + ,0 + ,397 + ,10 + ,62710 + ,322327 + ,0 + ,450 + ,8 + ,59325 + ,292136 + ,0 + ,622 + ,7 + ,60950 + ,263276 + ,0 + ,694 + ,8 + ,68060 + ,367655 + ,1 + ,3425 + ,12 + ,83620 + ,283910 + ,0 + ,562 + ,13 + ,58456 + ,283587 + ,0 + ,4917 + ,42 + ,52811 + ,243650 + ,1 + ,1442 + ,118 + ,121173 + ,438493 + ,0 + ,529 + ,9 + ,63870 + ,296261 + ,1 + ,2126 + ,138 + ,21001 + ,230621 + ,0 + ,1061 + ,5 + ,70415 + ,304252 + ,0 + ,776 + ,9 + ,64230 + ,333505 + ,0 + ,611 + ,8 + ,59190 + ,296919 + ,1 + ,1526 + ,25 + ,69351 + ,278990 + ,0 + ,592 + ,7 + ,64270 + ,276898 + ,0 + ,1182 + ,13 + ,70694 + ,327007 + ,0 + ,621 + ,16 + ,68005 + ,317046 + ,0 + ,989 + ,11 + ,58930 + ,304555 + ,0 + ,438 + ,11 + ,58320 + ,298096 + ,0 + ,726 + ,3 + ,69980 + ,231861 + ,0 + ,1303 + ,61 + ,69863 + ,309422 + ,1 + ,7419 + ,29 + ,63255 + ,286963 + ,1 + ,1164 + ,17 + ,57320 + ,269753 + ,1 + ,3310 + ,33 + ,75230 + ,448243 + ,0 + ,1920 + ,15 + ,79420 + ,165404 + ,0 + ,965 + ,3 + ,73490 + ,204325 + ,0 + ,3256 + ,66 + ,35250 + ,407159 + ,1 + ,1135 + ,17 + ,62285 + ,290476 + ,0 + ,1270 + ,26 + ,69206 + ,275311 + ,0 + ,661 + ,3 + ,65920 + ,246541 + ,0 + ,1013 + ,2 + ,69770 + ,253468 + ,0 + ,2844 + ,67 + ,72683 + ,240897 + ,1 + ,11528 + ,70 + ,-14545 + ,-83265 + ,0 + ,6526 + ,26 + ,55830 + ,-42143 + ,0 + ,2264 + ,24 + ,55174 + ,272713 + ,1 + ,5109 + ,97 + ,67038 + ,215362 + ,0 + ,3999 + ,30 + ,51252 + ,42754 + ,0 + ,35624 + ,223 + ,157278 + ,306275 + ,0 + ,9252 + ,48 + ,79510 + ,253537 + ,0 + ,15236 + ,90 + ,77440 + ,372631 + ,0 + ,18073 + ,180 + ,27284 + ,-7170) + ,dim=c(5 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Trades' + ,'Dividends' + ,'Wealth ') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('Group','Costs','Trades','Dividends','Wealth '),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 = 'Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '5' > #'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 Wealth\r Group Costs Trades Dividends t 1 6282154 1 162556 1081 213118 1 2 4321023 1 29790 309 81767 2 3 4111912 1 87550 458 153198 3 4 223193 0 84738 588 -26007 4 5 1491348 1 54660 302 126942 5 6 1629616 1 42634 156 157214 6 7 1398893 0 40949 481 129352 7 8 1926517 1 45187 353 234817 8 9 983660 1 37704 452 60448 9 10 1443586 1 16275 109 47818 10 11 1073089 0 25830 115 245546 11 12 984885 0 12679 110 48020 12 13 1405225 1 18014 239 -1710 13 14 227132 0 43556 247 32648 14 15 929118 1 24811 505 95350 15 16 1071292 0 6575 159 151352 16 17 638830 0 7123 109 288170 17 18 856956 1 21950 519 114337 18 19 992426 1 37597 248 37884 19 20 444477 0 17821 373 122844 20 21 857217 1 12988 119 82340 21 22 711969 1 22330 84 79801 22 23 702380 0 13326 102 165548 23 24 358589 0 16189 295 116384 24 25 297978 0 7146 105 134028 25 26 585715 0 15824 64 63838 26 27 657954 1 27664 282 74996 27 28 209458 0 11920 182 31080 28 29 786690 0 8568 37 32168 29 30 439798 0 14416 361 49857 30 31 688779 1 3369 28 87161 31 32 574339 1 11819 85 106113 32 33 741409 1 6984 45 80570 33 34 597793 1 4519 49 102129 34 35 644190 0 2220 22 301670 35 36 377934 0 18562 155 102313 36 37 640273 0 10327 91 88577 37 38 697458 1 5336 81 112477 38 39 550608 1 2365 79 191778 39 40 207393 0 4069 145 79804 40 41 301607 0 8636 855 128294 41 42 345783 0 13718 61 96448 42 43 501749 0 4525 226 93811 43 44 379983 0 6869 105 117520 44 45 387475 0 4628 62 69159 45 46 377305 1 3689 25 101792 46 47 370837 1 4891 217 210568 47 48 430866 1 7489 322 136996 48 49 469107 0 4901 84 121920 49 50 194493 0 2284 33 76403 50 51 530670 1 3160 108 108094 51 52 518365 1 4150 150 134759 52 53 491303 1 7285 115 188873 53 54 527021 1 1134 162 146216 54 55 233773 1 4658 158 156608 55 56 405972 0 2384 97 61348 56 57 652925 0 3748 9 50350 57 58 446211 0 5371 66 87720 58 59 341340 0 1285 107 99489 59 60 387699 1 9327 101 87419 60 61 493408 1 5565 47 94355 61 62 146494 0 1528 38 60326 62 63 414462 1 3122 34 94670 63 64 364304 1 7561 87 82425 64 65 355178 0 2675 79 59017 65 66 357760 0 13253 947 90829 66 67 261216 0 880 74 80791 67 68 397144 1 2053 53 100423 68 69 374943 0 1424 94 131116 69 70 424898 1 4036 63 100269 70 71 202055 1 3045 58 27330 71 72 378525 0 5119 49 39039 72 73 310768 0 1431 34 106885 73 74 325738 0 554 11 79285 74 75 394510 0 1975 35 118881 75 76 247060 1 1765 20 77623 76 77 368078 0 1012 47 114768 77 78 236761 0 810 43 74015 78 79 312378 0 1280 117 69465 79 80 339836 1 666 171 117869 80 81 347385 0 1380 26 60982 81 82 426280 1 4677 75 90131 82 83 352850 0 876 59 138971 83 84 301881 0 814 18 39625 84 85 377516 0 514 15 102725 85 86 357312 1 5692 72 64239 86 87 458343 0 3642 86 90262 87 88 354228 0 540 14 103960 88 89 308636 0 2099 64 106611 89 90 386212 0 567 11 103345 90 91 393343 0 2001 52 95551 91 92 378509 1 2949 41 82903 92 93 452469 0 2253 99 63593 93 94 364839 1 6533 75 126910 94 95 358649 0 1889 45 37527 95 96 376641 1 3055 43 60247 96 97 429112 0 272 8 112995 97 98 330546 1 1414 198 70184 98 99 403560 0 2564 22 130140 99 100 317892 1 1383 11 73221 100 101 307528 1 1261 33 76114 101 102 235133 0 975 23 90534 102 103 299243 0 3366 80 108479 103 104 314073 0 576 18 113761 104 105 368186 0 1686 40 68696 105 106 269661 0 746 23 71561 106 107 125390 1 3192 60 59831 107 108 510834 1 2045 20 97890 108 109 321896 0 5702 61 101481 109 110 249898 1 1932 36 72954 110 111 408881 0 936 30 67939 111 112 158492 1 3437 47 48022 112 113 292154 0 5131 71 86111 113 114 289513 1 2397 14 74020 114 115 378049 1 1389 9 57530 115 116 343466 0 1503 39 56364 116 117 332743 0 402 26 84990 117 118 442882 0 2239 21 88590 118 119 214215 1 2234 16 77200 119 120 315688 0 837 69 61262 120 121 375195 0 10579 92 110309 121 122 334280 0 875 14 67000 122 123 355864 0 1585 107 93099 123 124 480382 1 1659 29 107577 124 125 353058 1 2647 37 62920 125 126 217193 1 3294 23 75832 126 127 315380 0 0 0 60720 127 128 314533 0 94 7 60793 128 129 318056 0 422 28 57935 129 130 315380 0 0 0 60720 130 131 314353 0 34 8 60630 131 132 369448 0 1558 63 55637 132 133 315380 0 0 0 60720 133 134 312846 0 43 3 60887 134 135 312075 0 645 5 60720 135 136 315009 0 316 9 60505 136 137 318903 0 115 13 60945 137 138 314887 0 5 2 60720 138 139 314913 0 897 5 60720 139 140 315380 0 0 0 60720 140 141 325506 0 389 14 58990 141 142 315380 0 0 0 60720 142 143 298568 0 1002 15 56750 143 144 315834 0 36 3 60894 144 145 329784 0 460 15 63346 145 146 312878 0 309 11 56535 146 147 315380 0 0 0 60720 147 148 314987 0 9 6 60835 148 149 325249 0 271 2 60720 149 150 315877 0 14 1 61016 150 151 291650 0 520 10 58650 151 152 305959 0 1766 73 60438 152 153 315380 0 0 0 60720 153 154 297765 0 458 11 58625 154 155 315245 0 20 3 60938 155 156 315380 0 0 0 60720 156 157 315380 0 0 0 60720 157 158 315236 0 98 2 61490 158 159 336425 0 405 7 60845 159 160 315380 0 0 0 60720 160 161 315380 0 0 0 60720 161 162 315380 0 0 0 60720 162 163 315380 0 0 0 60720 163 164 306268 0 483 27 60830 164 165 302187 0 454 51 63261 165 166 314882 0 47 3 60720 166 167 315380 0 0 0 60720 167 168 382712 0 757 19 45689 168 169 341570 0 4655 393 60720 169 170 315380 0 0 0 60720 170 171 315380 0 0 0 60720 171 172 312412 0 36 4 61564 172 173 315380 0 0 0 60720 173 174 309596 0 203 9 61938 174 175 315380 0 0 0 60720 175 176 315547 0 126 10 60951 176 177 313267 0 400 152 60720 177 178 316176 1 71 1 60745 178 179 315380 0 0 0 60720 179 180 315380 0 0 0 60720 180 181 359335 0 972 34 71642 181 182 330068 1 531 10 71641 182 183 314289 0 2461 57 55792 183 184 297413 1 378 52 71873 184 185 314806 1 23 5 62555 185 186 333210 1 638 14 60370 186 187 352108 1 2300 29 64873 187 188 313332 0 149 5 62041 188 189 291787 0 226 5 65745 189 190 315380 0 0 0 60720 190 191 318745 0 275 4 59500 191 192 315380 0 0 0 60720 192 193 315366 0 141 6 61630 193 194 315380 0 0 0 60720 194 195 315688 0 28 2 60890 195 196 315380 1 0 0 60720 196 197 409642 1 4980 91 113521 197 198 315380 1 0 0 60720 198 199 315380 1 0 0 60720 199 200 269587 1 472 20 80045 200 201 315380 1 0 0 60720 201 202 315380 1 0 0 60720 202 203 315380 1 0 0 60720 203 204 300962 1 203 27 50804 204 205 325479 1 496 17 87390 205 206 316155 1 10 2 61656 206 207 318574 1 63 4 65688 207 208 315380 1 0 0 60720 208 209 343613 1 1136 32 48522 209 210 306948 1 265 31 60720 210 211 315380 0 0 0 60720 211 212 315380 0 0 0 60720 212 213 330059 0 267 32 57640 213 214 288985 0 474 20 61977 214 215 304485 0 534 7 62620 215 216 315380 0 0 0 60720 216 217 315688 0 15 8 60831 217 218 317736 0 397 28 60646 218 219 315380 0 0 0 60720 219 220 322331 0 1866 29 56225 220 221 296656 0 288 4 60510 221 222 315380 0 0 0 60720 222 223 315354 0 3 2 60698 223 224 312161 0 468 21 60720 224 225 315576 0 20 2 60805 225 226 314922 0 278 26 61404 226 227 314551 0 61 14 60720 227 228 315380 0 0 0 60720 228 229 312339 0 192 4 65276 229 230 315380 0 0 0 60720 230 231 298700 0 317 9 63915 231 232 321376 0 738 10 60720 232 233 315380 0 0 0 60720 233 234 303230 0 368 17 61686 234 235 315380 0 0 0 60720 235 236 315487 0 2 1 60743 236 237 315380 0 0 0 60720 237 238 315793 0 53 6 60349 238 239 315380 0 0 0 60720 239 240 315380 0 0 0 60720 240 241 315380 0 0 0 60720 241 242 312887 0 94 3 61360 242 243 315380 0 0 0 60720 243 244 315637 0 24 8 59818 244 245 324385 0 2332 4 72680 245 246 315380 1 0 0 60720 246 247 315380 1 0 0 60720 247 248 308989 0 131 11 61808 248 249 315380 1 0 0 60720 249 250 315380 1 0 0 60720 250 251 296702 0 206 9 53110 251 252 315380 1 0 0 60720 252 253 307322 0 167 2 64245 253 254 304376 0 622 73 73007 254 255 253588 0 2328 94 82732 255 256 315380 0 0 0 60720 256 257 309560 0 365 8 54820 257 258 298466 0 364 35 47705 258 259 315380 1 0 0 60720 259 260 315380 0 0 0 60720 260 261 315380 1 0 0 60720 261 262 315380 1 0 0 60720 262 263 343929 0 226 12 72835 263 264 331955 0 307 15 58856 264 265 315380 1 0 0 60720 265 266 315380 0 0 0 60720 266 267 315380 1 0 0 60720 267 268 381180 0 188 11 77655 268 269 315380 1 0 0 60720 269 270 331420 0 138 6 69817 270 271 315380 1 0 0 60720 271 272 315380 1 0 0 60720 272 273 315380 1 0 0 60720 273 274 310201 0 125 12 60798 274 275 315380 0 0 0 60720 275 276 320016 0 282 30 62452 276 277 320398 0 335 33 64175 277 278 315380 0 0 0 60720 278 279 291841 1 1324 117 67440 279 280 310670 0 176 28 68136 280 281 315380 1 0 0 60720 281 282 315380 0 0 0 60720 282 283 313491 0 249 72 56726 283 284 315380 0 0 0 60720 284 285 331323 0 333 13 70811 285 286 315380 0 0 0 60720 286 287 319210 1 601 6 60720 287 288 318098 0 30 4 62045 288 289 315380 1 0 0 60720 289 290 292754 0 249 62 54323 290 291 315380 1 0 0 60720 291 292 325176 0 165 24 62841 292 293 365959 0 453 21 81125 293 294 315380 0 0 0 60720 294 295 302409 0 53 14 59506 295 296 340968 1 382 21 59365 296 297 315380 0 0 0 60720 297 298 315380 1 0 0 60720 298 299 315380 0 0 0 60720 299 300 315380 0 0 0 60720 300 301 313164 1 30 4 60798 301 302 301164 0 290 2 58790 302 303 315380 0 0 0 60720 303 304 315380 1 0 0 60720 304 305 344425 0 366 53 61808 305 306 315394 0 2 9 60735 306 307 315380 1 0 0 60720 307 308 316647 1 209 13 64016 308 309 309836 0 384 22 54683 309 310 315380 1 0 0 60720 310 311 315380 1 0 0 60720 311 312 346611 0 365 83 87192 312 313 315380 1 0 0 60720 313 314 322031 1 49 8 64107 314 315 315656 0 3 4 60761 315 316 339445 0 133 14 65990 316 317 314964 0 32 1 59988 317 318 297141 0 368 17 61167 318 319 315372 0 1 6 60719 319 320 315380 0 0 0 60720 320 321 315380 0 0 0 60720 321 322 315380 0 0 0 60720 322 323 315380 0 0 0 60720 323 324 315380 1 0 0 60720 324 325 315380 0 0 0 60720 325 326 312502 0 22 2 60722 326 327 315380 1 0 0 60720 327 328 315380 1 0 0 60720 328 329 315380 0 0 0 60720 329 330 315380 0 0 0 60720 330 331 315380 0 0 0 60720 331 332 315380 0 0 0 60720 332 333 315380 0 0 0 60720 333 334 313729 0 96 5 60379 334 335 315388 0 1 2 60727 335 336 315371 0 314 5 60720 336 337 296139 0 844 78 60925 337 338 315380 0 0 0 60720 338 339 313880 0 26 1 60896 339 340 317698 0 125 13 59734 340 341 295580 0 304 15 62969 341 342 315380 0 0 0 60720 342 343 315380 0 0 0 60720 343 344 315380 0 0 0 60720 344 345 308256 0 621 48 60720 345 346 315380 0 0 0 60720 346 347 303677 0 119 6 59118 347 348 315380 0 0 0 60720 348 349 315380 0 0 0 60720 349 350 319369 0 1595 17 60720 350 351 318690 0 312 14 58598 351 352 314049 0 60 10 61124 352 353 325699 0 587 12 59595 353 354 314210 0 135 2 62065 354 355 315380 0 0 0 60720 355 356 315380 0 0 0 60720 356 357 322378 0 514 52 78780 357 358 315380 0 0 0 60720 358 359 315380 0 0 0 60720 359 360 315380 0 0 0 60720 360 361 315398 0 1 4 60722 361 362 315380 0 0 0 60720 362 363 315380 0 0 0 60720 363 364 308336 1 1763 24 61600 364 365 316386 0 180 11 59635 365 366 315380 0 0 0 60720 366 367 315380 0 0 0 60720 367 368 315380 0 0 0 60720 368 369 315380 0 0 0 60720 369 370 315553 0 218 21 60720 370 371 315380 0 0 0 60720 371 372 323361 0 448 40 59781 372 373 336639 0 227 9 76644 373 374 307424 0 174 1 64820 374 375 315380 0 0 0 60720 375 376 315380 0 0 0 60720 376 377 295370 0 121 24 56178 377 378 322340 0 607 11 60436 378 379 319864 0 2212 14 60720 379 380 315380 0 0 0 60720 380 381 315380 0 0 0 60720 381 382 317291 0 530 60 73433 382 383 280398 0 571 80 41477 383 384 315380 0 0 0 60720 384 385 317330 0 78 16 62700 385 386 238125 0 2489 40 67804 386 387 327071 0 131 6 59661 387 388 309038 0 923 8 58620 388 389 314210 0 72 3 60398 389 390 307930 0 572 16 58580 390 391 322327 0 397 10 62710 391 392 292136 0 450 8 59325 392 393 263276 0 622 7 60950 393 394 367655 0 694 8 68060 394 395 283910 1 3425 12 83620 395 396 283587 0 562 13 58456 396 397 243650 0 4917 42 52811 397 398 438493 1 1442 118 121173 398 399 296261 0 529 9 63870 399 400 230621 1 2126 138 21001 400 401 304252 0 1061 5 70415 401 402 333505 0 776 9 64230 402 403 296919 0 611 8 59190 403 404 278990 1 1526 25 69351 404 405 276898 0 592 7 64270 405 406 327007 0 1182 13 70694 406 407 317046 0 621 16 68005 407 408 304555 0 989 11 58930 408 409 298096 0 438 11 58320 409 410 231861 0 726 3 69980 410 411 309422 0 1303 61 69863 411 412 286963 1 7419 29 63255 412 413 269753 1 1164 17 57320 413 414 448243 1 3310 33 75230 414 415 165404 0 1920 15 79420 415 416 204325 0 965 3 73490 416 417 407159 0 3256 66 35250 417 418 290476 1 1135 17 62285 418 419 275311 0 1270 26 69206 419 420 246541 0 661 3 65920 420 421 253468 0 1013 2 69770 421 422 240897 0 2844 67 72683 422 423 -83265 1 11528 70 -14545 423 424 -42143 0 6526 26 55830 424 425 272713 0 2264 24 55174 425 426 215362 1 5109 97 67038 426 427 42754 0 3999 30 51252 427 428 306275 0 35624 223 157278 428 429 253537 0 9252 48 79510 429 430 372631 0 15236 90 77440 430 431 -7170 0 18073 180 27284 431 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Dividends t 175770.913 53021.362 29.801 -431.749 1.838 -78.375 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2175875 -15297 35999 51627 3187769 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.758e+05 4.820e+04 3.647 0.000299 *** Group 5.302e+04 2.809e+04 1.888 0.059763 . Costs 2.980e+01 1.579e+00 18.870 < 2e-16 *** Trades -4.317e+02 1.703e+02 -2.536 0.011575 * Dividends 1.838e+00 4.307e-01 4.267 2.45e-05 *** t -7.837e+01 1.145e+02 -0.684 0.494218 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 253500 on 425 degrees of freedom Multiple R-squared: 0.6573, Adjusted R-squared: 0.6532 F-statistic: 163 on 5 and 425 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 1.515478e-18 7.577388e-19 [2,] 1.0000000 4.419952e-64 2.209976e-64 [3,] 1.0000000 4.935590e-79 2.467795e-79 [4,] 1.0000000 9.104659e-109 4.552329e-109 [5,] 1.0000000 1.680338e-129 8.401691e-130 [6,] 1.0000000 3.121883e-133 1.560941e-133 [7,] 1.0000000 3.217971e-136 1.608985e-136 [8,] 1.0000000 3.820751e-151 1.910375e-151 [9,] 1.0000000 4.497188e-150 2.248594e-150 [10,] 1.0000000 9.318449e-153 4.659224e-153 [11,] 1.0000000 1.161379e-158 5.806893e-159 [12,] 1.0000000 6.858018e-159 3.429009e-159 [13,] 1.0000000 1.208447e-164 6.042236e-165 [14,] 1.0000000 6.215520e-166 3.107760e-166 [15,] 1.0000000 1.817677e-169 9.088387e-170 [16,] 1.0000000 1.297091e-169 6.485457e-170 [17,] 1.0000000 1.694399e-170 8.471993e-171 [18,] 1.0000000 1.406974e-174 7.034868e-175 [19,] 1.0000000 1.561645e-176 7.808225e-177 [20,] 1.0000000 5.798162e-178 2.899081e-178 [21,] 1.0000000 1.993953e-192 9.969767e-193 [22,] 1.0000000 1.315769e-193 6.578846e-194 [23,] 1.0000000 4.356192e-197 2.178096e-197 [24,] 1.0000000 6.875675e-198 3.437838e-198 [25,] 1.0000000 5.143499e-205 2.571749e-205 [26,] 1.0000000 1.611490e-206 8.057452e-207 [27,] 1.0000000 2.024035e-205 1.012018e-205 [28,] 1.0000000 1.343738e-204 6.718692e-205 [29,] 1.0000000 4.749207e-211 2.374604e-211 [30,] 1.0000000 1.690935e-216 8.454676e-217 [31,] 1.0000000 8.165241e-216 4.082621e-216 [32,] 1.0000000 1.032994e-216 5.164970e-217 [33,] 1.0000000 1.818056e-216 9.090279e-217 [34,] 1.0000000 1.311054e-215 6.555271e-216 [35,] 1.0000000 3.547422e-217 1.773711e-217 [36,] 1.0000000 2.262280e-216 1.131140e-216 [37,] 1.0000000 1.882959e-216 9.414796e-217 [38,] 1.0000000 1.865443e-215 9.327213e-216 [39,] 1.0000000 5.660137e-217 2.830068e-217 [40,] 1.0000000 5.208572e-216 2.604286e-216 [41,] 1.0000000 5.039992e-216 2.519996e-216 [42,] 1.0000000 4.853354e-217 2.426677e-217 [43,] 1.0000000 1.602300e-217 8.011499e-218 [44,] 1.0000000 3.142225e-217 1.571112e-217 [45,] 1.0000000 1.847783e-216 9.238914e-217 [46,] 1.0000000 5.500221e-216 2.750110e-216 [47,] 1.0000000 1.261439e-218 6.307196e-219 [48,] 1.0000000 1.472648e-219 7.363239e-220 [49,] 1.0000000 3.768798e-232 1.884399e-232 [50,] 1.0000000 1.082690e-232 5.413451e-233 [51,] 1.0000000 3.711894e-232 1.855947e-232 [52,] 1.0000000 1.143142e-231 5.715709e-232 [53,] 1.0000000 8.744860e-233 4.372430e-233 [54,] 1.0000000 8.764647e-235 4.382323e-235 [55,] 1.0000000 5.489724e-234 2.744862e-234 [56,] 1.0000000 3.315814e-233 1.657907e-233 [57,] 1.0000000 6.575947e-233 3.287974e-233 [58,] 1.0000000 4.504210e-232 2.252105e-232 [59,] 1.0000000 4.223847e-232 2.111924e-232 [60,] 1.0000000 4.861336e-231 2.430668e-231 [61,] 1.0000000 1.962121e-230 9.810604e-231 [62,] 1.0000000 8.842435e-230 4.421218e-230 [63,] 1.0000000 5.884206e-229 2.942103e-229 [64,] 1.0000000 7.544807e-230 3.772404e-230 [65,] 1.0000000 1.554422e-229 7.772108e-230 [66,] 1.0000000 8.834832e-229 4.417416e-229 [67,] 1.0000000 7.675549e-228 3.837774e-228 [68,] 1.0000000 1.076490e-227 5.382448e-228 [69,] 1.0000000 6.193443e-227 3.096721e-227 [70,] 1.0000000 3.306115e-227 1.653058e-227 [71,] 1.0000000 2.062552e-226 1.031276e-226 [72,] 1.0000000 3.482571e-226 1.741286e-226 [73,] 1.0000000 1.898070e-225 9.490351e-226 [74,] 1.0000000 4.012425e-225 2.006213e-225 [75,] 1.0000000 2.208483e-225 1.104242e-225 [76,] 1.0000000 1.344902e-224 6.724510e-225 [77,] 1.0000000 1.083559e-223 5.417793e-224 [78,] 1.0000000 3.431825e-223 1.715913e-223 [79,] 1.0000000 1.620583e-223 8.102917e-224 [80,] 1.0000000 9.558939e-223 4.779470e-223 [81,] 1.0000000 1.350018e-222 6.750091e-223 [82,] 1.0000000 1.218524e-221 6.092618e-222 [83,] 1.0000000 1.130777e-220 5.653884e-221 [84,] 1.0000000 8.945680e-220 4.472840e-220 [85,] 1.0000000 3.543613e-221 1.771806e-221 [86,] 1.0000000 1.740925e-220 8.704624e-221 [87,] 1.0000000 9.237789e-221 4.618895e-221 [88,] 1.0000000 1.509665e-220 7.548325e-221 [89,] 1.0000000 1.426899e-219 7.134495e-220 [90,] 1.0000000 1.709761e-218 8.548805e-219 [91,] 1.0000000 1.347585e-217 6.737925e-218 [92,] 1.0000000 1.457848e-216 7.289241e-217 [93,] 1.0000000 1.225221e-215 6.126103e-216 [94,] 1.0000000 4.867287e-217 2.433643e-217 [95,] 1.0000000 2.612817e-217 1.306409e-217 [96,] 1.0000000 2.196936e-218 1.098468e-218 [97,] 1.0000000 1.293048e-217 6.465242e-218 [98,] 1.0000000 2.754812e-217 1.377406e-217 [99,] 1.0000000 2.822821e-220 1.411410e-220 [100,] 1.0000000 3.184031e-221 1.592015e-221 [101,] 1.0000000 1.683356e-220 8.416779e-221 [102,] 1.0000000 1.590035e-220 7.950175e-221 [103,] 1.0000000 2.031209e-220 1.015604e-220 [104,] 1.0000000 3.355586e-221 1.677793e-221 [105,] 1.0000000 1.581953e-220 7.909765e-221 [106,] 1.0000000 8.652811e-220 4.326406e-220 [107,] 1.0000000 1.638948e-219 8.194741e-220 [108,] 1.0000000 9.051606e-219 4.525803e-219 [109,] 1.0000000 4.368337e-218 2.184169e-218 [110,] 1.0000000 7.049943e-218 3.524971e-218 [111,] 1.0000000 1.475061e-219 7.375304e-220 [112,] 1.0000000 1.445986e-218 7.229931e-219 [113,] 1.0000000 1.334807e-217 6.674035e-218 [114,] 1.0000000 1.379293e-216 6.896466e-217 [115,] 1.0000000 8.159952e-216 4.079976e-216 [116,] 1.0000000 2.805980e-215 1.402990e-215 [117,] 1.0000000 1.621156e-214 8.105779e-215 [118,] 1.0000000 1.033328e-215 5.166640e-216 [119,] 1.0000000 1.005872e-214 5.029362e-215 [120,] 1.0000000 9.861832e-214 4.930916e-214 [121,] 1.0000000 9.986274e-213 4.993137e-213 [122,] 1.0000000 1.000083e-211 5.000417e-212 [123,] 1.0000000 9.970742e-211 4.985371e-211 [124,] 1.0000000 1.941811e-210 9.709054e-211 [125,] 1.0000000 2.007670e-209 1.003835e-209 [126,] 1.0000000 2.051664e-208 1.025832e-208 [127,] 1.0000000 2.227086e-207 1.113543e-207 [128,] 1.0000000 2.372285e-206 1.186142e-206 [129,] 1.0000000 2.501003e-205 1.250502e-205 [130,] 1.0000000 2.606333e-204 1.303166e-204 [131,] 1.0000000 2.893069e-203 1.446534e-203 [132,] 1.0000000 3.029840e-202 1.514920e-202 [133,] 1.0000000 3.152845e-201 1.576422e-201 [134,] 1.0000000 3.308396e-200 1.654198e-200 [135,] 1.0000000 3.584480e-199 1.792240e-199 [136,] 1.0000000 3.753738e-198 1.876869e-198 [137,] 1.0000000 4.017008e-197 2.008504e-197 [138,] 1.0000000 4.294421e-196 2.147211e-196 [139,] 1.0000000 4.483868e-195 2.241934e-195 [140,] 1.0000000 4.636728e-194 2.318364e-194 [141,] 1.0000000 4.902323e-193 2.451162e-193 [142,] 1.0000000 5.089050e-192 2.544525e-192 [143,] 1.0000000 4.536544e-191 2.268272e-191 [144,] 1.0000000 4.810067e-190 2.405033e-190 [145,] 1.0000000 4.943833e-189 2.471917e-189 [146,] 1.0000000 4.693623e-188 2.346812e-188 [147,] 1.0000000 4.768557e-187 2.384279e-187 [148,] 1.0000000 4.850010e-186 2.425005e-186 [149,] 1.0000000 4.916451e-185 2.458226e-185 [150,] 1.0000000 4.938567e-184 2.469284e-184 [151,] 1.0000000 4.691499e-183 2.345749e-183 [152,] 1.0000000 4.718413e-182 2.359206e-182 [153,] 1.0000000 4.727624e-181 2.363812e-181 [154,] 1.0000000 4.718569e-180 2.359285e-180 [155,] 1.0000000 4.690923e-179 2.345461e-179 [156,] 1.0000000 4.393029e-178 2.196514e-178 [157,] 1.0000000 3.414431e-177 1.707215e-177 [158,] 1.0000000 3.341140e-176 1.670570e-176 [159,] 1.0000000 3.260951e-175 1.630475e-175 [160,] 1.0000000 1.036608e-175 5.183042e-176 [161,] 1.0000000 5.933276e-175 2.966638e-175 [162,] 1.0000000 5.881166e-174 2.940583e-174 [163,] 1.0000000 5.803678e-173 2.901839e-173 [164,] 1.0000000 5.532135e-172 2.766068e-172 [165,] 1.0000000 5.408149e-171 2.704074e-171 [166,] 1.0000000 5.001847e-170 2.500923e-170 [167,] 1.0000000 4.842125e-169 2.421062e-169 [168,] 1.0000000 4.683089e-168 2.341544e-168 [169,] 1.0000000 4.458892e-167 2.229446e-167 [170,] 1.0000000 4.033280e-166 2.016640e-166 [171,] 1.0000000 3.835111e-165 1.917555e-165 [172,] 1.0000000 3.629539e-164 1.814769e-164 [173,] 1.0000000 3.132239e-163 1.566119e-163 [174,] 1.0000000 2.565970e-162 1.282985e-162 [175,] 1.0000000 2.019233e-161 1.009617e-161 [176,] 1.0000000 7.887108e-161 3.943554e-161 [177,] 1.0000000 6.739485e-160 3.369742e-160 [178,] 1.0000000 5.611098e-159 2.805549e-159 [179,] 1.0000000 3.006365e-158 1.503183e-158 [180,] 1.0000000 2.722168e-157 1.361084e-157 [181,] 1.0000000 1.485235e-156 7.426177e-157 [182,] 1.0000000 1.358335e-155 6.791673e-156 [183,] 1.0000000 1.256463e-154 6.282317e-155 [184,] 1.0000000 1.139527e-153 5.697634e-154 [185,] 1.0000000 1.024080e-152 5.120398e-153 [186,] 1.0000000 9.201480e-152 4.600740e-152 [187,] 1.0000000 8.234211e-151 4.117106e-151 [188,] 1.0000000 7.049946e-150 3.524973e-150 [189,] 1.0000000 4.862527e-149 2.431264e-149 [190,] 1.0000000 4.177865e-148 2.088933e-148 [191,] 1.0000000 3.580136e-147 1.790068e-147 [192,] 1.0000000 1.421982e-147 7.109909e-148 [193,] 1.0000000 1.216127e-146 6.080634e-147 [194,] 1.0000000 1.036096e-145 5.180481e-146 [195,] 1.0000000 8.791750e-145 4.395875e-145 [196,] 1.0000000 7.669867e-144 3.834933e-144 [197,] 1.0000000 1.804365e-143 9.021823e-144 [198,] 1.0000000 1.489763e-142 7.448815e-143 [199,] 1.0000000 1.124895e-141 5.624473e-142 [200,] 1.0000000 9.285318e-141 4.642659e-141 [201,] 1.0000000 2.830105e-140 1.415053e-140 [202,] 1.0000000 2.233113e-139 1.116556e-139 [203,] 1.0000000 1.897353e-138 9.486763e-139 [204,] 1.0000000 1.605142e-137 8.025712e-138 [205,] 1.0000000 1.236876e-136 6.184382e-137 [206,] 1.0000000 7.594920e-136 3.797460e-136 [207,] 1.0000000 5.802934e-135 2.901467e-135 [208,] 1.0000000 4.828425e-134 2.414213e-134 [209,] 1.0000000 4.001215e-133 2.000607e-133 [210,] 1.0000000 3.358999e-132 1.679500e-132 [211,] 1.0000000 2.759562e-131 1.379781e-131 [212,] 1.0000000 1.859136e-130 9.295678e-131 [213,] 1.0000000 1.320228e-129 6.601142e-130 [214,] 1.0000000 1.073671e-128 5.368357e-129 [215,] 1.0000000 8.694309e-128 4.347155e-128 [216,] 1.0000000 7.040859e-127 3.520430e-127 [217,] 1.0000000 5.653002e-126 2.826501e-126 [218,] 1.0000000 4.522004e-125 2.261002e-125 [219,] 1.0000000 3.591051e-124 1.795525e-124 [220,] 1.0000000 2.842281e-123 1.421141e-123 [221,] 1.0000000 2.009232e-122 1.004616e-122 [222,] 1.0000000 1.574298e-121 7.871490e-122 [223,] 1.0000000 9.592876e-121 4.796438e-121 [224,] 1.0000000 7.566092e-120 3.783046e-120 [225,] 1.0000000 5.845772e-119 2.922886e-119 [226,] 1.0000000 4.101848e-118 2.050924e-118 [227,] 1.0000000 3.136694e-117 1.568347e-117 [228,] 1.0000000 2.388046e-116 1.194023e-116 [229,] 1.0000000 1.808837e-115 9.044183e-116 [230,] 1.0000000 1.374516e-114 6.872582e-115 [231,] 1.0000000 1.031393e-113 5.156967e-114 [232,] 1.0000000 7.701906e-113 3.850953e-113 [233,] 1.0000000 5.723502e-112 2.861751e-112 [234,] 1.0000000 4.151705e-111 2.075853e-111 [235,] 1.0000000 3.054121e-110 1.527061e-110 [236,] 1.0000000 2.262317e-109 1.131158e-109 [237,] 1.0000000 1.577420e-108 7.887100e-109 [238,] 1.0000000 1.117470e-107 5.587352e-108 [239,] 1.0000000 7.885128e-107 3.942564e-107 [240,] 1.0000000 5.358671e-106 2.679335e-106 [241,] 1.0000000 3.747670e-105 1.873835e-105 [242,] 1.0000000 2.609784e-104 1.304892e-104 [243,] 1.0000000 1.870328e-103 9.351641e-104 [244,] 1.0000000 1.291012e-102 6.455060e-103 [245,] 1.0000000 7.857566e-102 3.928783e-102 [246,] 1.0000000 3.186289e-101 1.593144e-101 [247,] 1.0000000 6.151534e-102 3.075767e-102 [248,] 1.0000000 4.270225e-101 2.135112e-101 [249,] 1.0000000 3.086918e-100 1.543459e-100 [250,] 1.0000000 2.205365e-99 1.102682e-99 [251,] 1.0000000 1.507519e-98 7.537593e-99 [252,] 1.0000000 1.031568e-97 5.157840e-98 [253,] 1.0000000 6.987961e-97 3.493980e-97 [254,] 1.0000000 4.711436e-96 2.355718e-96 [255,] 1.0000000 3.116187e-95 1.558094e-95 [256,] 1.0000000 2.116939e-94 1.058470e-94 [257,] 1.0000000 1.410379e-93 7.051893e-94 [258,] 1.0000000 9.423939e-93 4.711970e-93 [259,] 1.0000000 6.219654e-92 3.109827e-92 [260,] 1.0000000 4.156511e-91 2.078255e-91 [261,] 1.0000000 2.721886e-90 1.360943e-90 [262,] 1.0000000 1.693425e-89 8.467125e-90 [263,] 1.0000000 1.097962e-88 5.489808e-89 [264,] 1.0000000 7.080714e-88 3.540357e-88 [265,] 1.0000000 4.540766e-87 2.270383e-87 [266,] 1.0000000 2.829782e-86 1.414891e-86 [267,] 1.0000000 1.817427e-85 9.087134e-86 [268,] 1.0000000 1.161385e-84 5.806923e-85 [269,] 1.0000000 7.189676e-84 3.594838e-84 [270,] 1.0000000 4.546009e-83 2.273005e-83 [271,] 1.0000000 1.398810e-82 6.994051e-83 [272,] 1.0000000 6.275686e-82 3.137843e-82 [273,] 1.0000000 3.863207e-81 1.931604e-81 [274,] 1.0000000 2.398864e-80 1.199432e-80 [275,] 1.0000000 1.468775e-79 7.343877e-80 [276,] 1.0000000 9.028335e-79 4.514168e-79 [277,] 1.0000000 5.052306e-78 2.526153e-78 [278,] 1.0000000 3.068809e-77 1.534405e-77 [279,] 1.0000000 1.890460e-76 9.452298e-77 [280,] 1.0000000 1.129239e-75 5.646194e-76 [281,] 1.0000000 6.704396e-75 3.352198e-75 [282,] 1.0000000 3.492248e-74 1.746124e-74 [283,] 1.0000000 2.052099e-73 1.026050e-73 [284,] 1.0000000 1.216532e-72 6.082660e-73 [285,] 1.0000000 6.933466e-72 3.466733e-72 [286,] 1.0000000 4.045770e-71 2.022885e-71 [287,] 1.0000000 2.113795e-70 1.056898e-70 [288,] 1.0000000 1.209689e-69 6.048443e-70 [289,] 1.0000000 6.968100e-69 3.484050e-69 [290,] 1.0000000 3.988206e-68 1.994103e-68 [291,] 1.0000000 2.275480e-67 1.137740e-67 [292,] 1.0000000 1.291287e-66 6.456433e-67 [293,] 1.0000000 7.145543e-66 3.572771e-66 [294,] 1.0000000 3.816662e-65 1.908331e-65 [295,] 1.0000000 2.130522e-64 1.065261e-64 [296,] 1.0000000 1.179190e-63 5.895950e-64 [297,] 1.0000000 6.724510e-63 3.362255e-63 [298,] 1.0000000 3.674526e-62 1.837263e-62 [299,] 1.0000000 2.003543e-61 1.001771e-61 [300,] 1.0000000 1.020572e-60 5.102861e-61 [301,] 1.0000000 5.729089e-60 2.864545e-60 [302,] 1.0000000 3.052548e-59 1.526274e-59 [303,] 1.0000000 1.610343e-58 8.051713e-59 [304,] 1.0000000 3.626974e-58 1.813487e-58 [305,] 1.0000000 1.893155e-57 9.465775e-58 [306,] 1.0000000 9.417253e-57 4.708626e-57 [307,] 1.0000000 4.936750e-56 2.468375e-56 [308,] 1.0000000 2.663227e-55 1.331614e-55 [309,] 1.0000000 1.401157e-54 7.005784e-55 [310,] 1.0000000 6.002863e-54 3.001431e-54 [311,] 1.0000000 3.062509e-53 1.531254e-53 [312,] 1.0000000 1.568720e-52 7.843602e-53 [313,] 1.0000000 7.992001e-52 3.996001e-52 [314,] 1.0000000 4.049453e-51 2.024727e-51 [315,] 1.0000000 2.040594e-50 1.020297e-50 [316,] 1.0000000 9.858699e-50 4.929349e-50 [317,] 1.0000000 4.911602e-49 2.455801e-49 [318,] 1.0000000 2.388471e-48 1.194235e-48 [319,] 1.0000000 1.118447e-47 5.592235e-48 [320,] 1.0000000 5.091928e-47 2.545964e-47 [321,] 1.0000000 2.471670e-46 1.235835e-46 [322,] 1.0000000 1.192980e-45 5.964902e-46 [323,] 1.0000000 5.725285e-45 2.862643e-45 [324,] 1.0000000 2.731922e-44 1.365961e-44 [325,] 1.0000000 1.296081e-43 6.480404e-44 [326,] 1.0000000 6.078515e-43 3.039257e-43 [327,] 1.0000000 2.843432e-42 1.421716e-42 [328,] 1.0000000 1.333776e-41 6.668881e-42 [329,] 1.0000000 4.466083e-41 2.233042e-41 [330,] 1.0000000 2.055124e-40 1.027562e-40 [331,] 1.0000000 9.298580e-40 4.649290e-40 [332,] 1.0000000 4.273578e-39 2.136789e-39 [333,] 1.0000000 1.500689e-38 7.503443e-39 [334,] 1.0000000 6.723012e-38 3.361506e-38 [335,] 1.0000000 2.993973e-37 1.496986e-37 [336,] 1.0000000 1.325334e-36 6.626668e-37 [337,] 1.0000000 5.176679e-36 2.588339e-36 [338,] 1.0000000 2.260424e-35 1.130212e-35 [339,] 1.0000000 9.323828e-35 4.661914e-35 [340,] 1.0000000 4.016530e-34 2.008265e-34 [341,] 1.0000000 1.719686e-33 8.598431e-34 [342,] 1.0000000 7.574332e-33 3.787166e-33 [343,] 1.0000000 3.295568e-32 1.647784e-32 [344,] 1.0000000 1.356348e-31 6.781740e-32 [345,] 1.0000000 5.881285e-31 2.940642e-31 [346,] 1.0000000 2.407009e-30 1.203504e-30 [347,] 1.0000000 9.976592e-30 4.988296e-30 [348,] 1.0000000 4.108175e-29 2.054088e-29 [349,] 1.0000000 9.444924e-29 4.722462e-29 [350,] 1.0000000 3.834291e-28 1.917146e-28 [351,] 1.0000000 1.546504e-27 7.732519e-28 [352,] 1.0000000 6.196694e-27 3.098347e-27 [353,] 1.0000000 2.453025e-26 1.226513e-26 [354,] 1.0000000 9.696166e-26 4.848083e-26 [355,] 1.0000000 3.806424e-25 1.903212e-25 [356,] 1.0000000 1.336261e-24 6.681303e-25 [357,] 1.0000000 5.190094e-24 2.595047e-24 [358,] 1.0000000 1.991091e-23 9.955454e-24 [359,] 1.0000000 7.583449e-23 3.791725e-23 [360,] 1.0000000 2.867061e-22 1.433531e-22 [361,] 1.0000000 1.075795e-21 5.378977e-22 [362,] 1.0000000 3.946850e-21 1.973425e-21 [363,] 1.0000000 1.457539e-20 7.287697e-21 [364,] 1.0000000 5.355144e-20 2.677572e-20 [365,] 1.0000000 1.762319e-19 8.811597e-20 [366,] 1.0000000 5.903226e-19 2.951613e-19 [367,] 1.0000000 2.103424e-18 1.051712e-18 [368,] 1.0000000 7.430136e-18 3.715068e-18 [369,] 1.0000000 2.492973e-17 1.246486e-17 [370,] 1.0000000 8.711215e-17 4.355607e-17 [371,] 1.0000000 2.980005e-16 1.490002e-16 [372,] 1.0000000 1.015336e-15 5.076679e-16 [373,] 1.0000000 3.421759e-15 1.710880e-15 [374,] 1.0000000 9.315506e-15 4.657753e-15 [375,] 1.0000000 3.054829e-14 1.527415e-14 [376,] 1.0000000 9.975289e-14 4.987645e-14 [377,] 1.0000000 3.197171e-13 1.598586e-13 [378,] 1.0000000 5.175394e-13 2.587697e-13 [379,] 1.0000000 1.648750e-12 8.243748e-13 [380,] 1.0000000 5.197907e-12 2.598953e-12 [381,] 1.0000000 1.616802e-11 8.084009e-12 [382,] 1.0000000 4.954103e-11 2.477052e-11 [383,] 1.0000000 1.492874e-10 7.464372e-11 [384,] 1.0000000 4.398662e-10 2.199331e-10 [385,] 1.0000000 1.144726e-09 5.723628e-10 [386,] 1.0000000 2.983455e-09 1.491728e-09 [387,] 1.0000000 5.922027e-09 2.961013e-09 [388,] 1.0000000 1.656015e-08 8.280077e-09 [389,] 1.0000000 4.396224e-08 2.198112e-08 [390,] 1.0000000 7.238622e-08 3.619311e-08 [391,] 0.9999999 1.964606e-07 9.823030e-08 [392,] 0.9999998 4.450165e-07 2.225082e-07 [393,] 0.9999994 1.161090e-06 5.805450e-07 [394,] 0.9999985 2.984062e-06 1.492031e-06 [395,] 0.9999962 7.664921e-06 3.832461e-06 [396,] 0.9999928 1.430934e-05 7.154669e-06 [397,] 0.9999827 3.460435e-05 1.730218e-05 [398,] 0.9999577 8.460541e-05 4.230270e-05 [399,] 0.9998991 2.018981e-04 1.009490e-04 [400,] 0.9997735 4.530900e-04 2.265450e-04 [401,] 0.9995066 9.868163e-04 4.934082e-04 [402,] 0.9990026 1.994881e-03 9.974404e-04 [403,] 0.9979303 4.139357e-03 2.069679e-03 [404,] 0.9958082 8.383675e-03 4.191837e-03 [405,] 0.9915924 1.681520e-02 8.407600e-03 [406,] 0.9890337 2.193266e-02 1.096633e-02 [407,] 0.9856072 2.878567e-02 1.439284e-02 [408,] 0.9775658 4.486845e-02 2.243423e-02 [409,] 0.9881766 2.364676e-02 1.182338e-02 [410,] 0.9751101 4.977978e-02 2.488989e-02 [411,] 0.9514808 9.703838e-02 4.851919e-02 [412,] 0.9068904 1.862192e-01 9.310962e-02 [413,] 0.8341798 3.316403e-01 1.658202e-01 [414,] 0.7693017 4.613966e-01 2.306983e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1fdkq1291223383.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/2fdkq1291223383.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/3fdkq1291223383.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/484kt1291223383.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/584kt1291223383.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 1284205.5414 3187769.1732 1190500.6719 -2175874.6365 -468857.7504 6 7 8 9 10 -90790.2427 -26679.9252 72632.9742 -283972.6976 689755.2837 11 12 13 14 15 -273154.5235 391457.0462 746948.1834 -1198904.2326 4923.7717 16 17 18 19 20 491349.8015 -230376.0355 -10589.6859 -317846.2350 -325510.8117 21 22 23 24 25 143082.0979 -290933.1323 -128896.5763 -384255.3887 -289754.1051 26 27 28 29 30 -149268.0248 -409198.1778 -297881.3181 314718.7937 -98990.2136 31 32 33 34 35 213934.4307 -162462.4642 178442.5490 70473.2483 -139861.8286 36 37 38 39 40 -469274.7105 36163.3018 140904.3362 -63919.8915 -170551.2589 41 42 43 44 45 5074.2454 -386406.3673 119682.3108 -167669.1703 -23009.4475 46 47 48 49 50 -134082.1738 -293289.1576 -130071.2149 -36657.4800 -171578.8841 51 52 53 54 55 59693.6577 -12903.4165 -247866.8332 69915.9152 -349096.0264 56 57 58 59 60 92688.0486 281287.5672 -17778.3211 -4730.0110 -231383.0298 61 62 63 64 65 -49544.9854 -164405.2161 -61722.1836 -198703.3558 30439.4419 66 67 68 69 70 34164.2920 -52044.9352 -49159.9099 -38217.7441 -75743.8977 71 72 73 74 75 -137098.0592 5261.7630 -83664.6369 -1691.8845 -37590.3118 76 77 78 79 80 -162383.4934 -22428.0277 -74484.0321 27515.6407 -45306.8495 81 82 83 84 85 35998.9557 -68712.6880 -72428.5580 43390.1431 10792.7087 86 87 88 89 90 -121329.8845 52116.0046 -15736.2398 -90993.6726 15434.7945 91 92 93 94 95 11934.0118 -65600.8623 142726.8062 -252110.6991 84496.8895 96 97 98 99 100 -27816.9201 48646.1016 23808.4614 -70514.9164 -74082.9420 101 102 103 104 105 -76550.7077 -118139.3378 -133571.9433 -71993.9074 41430.7063 106 107 108 109 110 -41607.6573 -274184.3298 58310.5508 -175406.9829 -146369.5250 111 112 113 114 115 102020.2593 -231903.5840 -155256.9593 -131755.8649 15041.9739 116 117 118 119 120 45256.5727 9205.4005 55904.1942 -206784.6764 41591.1305 121 122 123 124 125 -269345.1492 24916.8122 17612.2972 26700.1079 -44470.8119 126 127 128 129 130 -229310.8781 37980.4618 37298.6434 45444.0513 38215.5862 131 132 133 134 135 39873.1091 82551.4642 38450.7106 35702.0045 18239.6111 136 137 138 139 140 33178.5791 44059.3673 39064.0778 13881.2785 38999.3343 141 142 143 144 145 46834.7729 39156.0840 6333.6478 39669.4956 41737.3438 146 147 148 149 150 40198.9096 39547.9580 41344.2876 42361.1560 39750.6735 151 152 153 154 155 8756.4069 9926.2905 40018.2069 17431.8785 40338.5765 156 157 158 159 160 40253.3314 40331.7062 36794.0965 52256.6186 40566.8306 161 162 163 164 165 40645.2054 40723.5803 40801.9551 28829.5624 31585.8041 166 167 168 169 170 40433.6828 41115.4543 121791.5073 98416.2146 41350.5788 171 172 173 174 175 41428.9536 37642.5124 41585.7032 31477.9707 41742.4528 176 177 178 179 180 42125.9025 93491.6577 -11977.8421 42055.9521 42134.3269 181 182 183 184 185 51809.8333 -27618.0822 1604.0199 -37849.6760 -12967.9289 186 187 188 189 190 -4912.1108 -37263.5792 36004.1943 5436.2354 42918.0750 191 192 193 194 195 42135.1276 43074.8247 39855.4995 43231.5743 43334.6202 196 197 198 199 200 -9633.0377 -121441.9611 -9476.2880 -9397.9132 -96056.2200 201 202 203 204 205 -9241.1636 -9162.7888 -9084.4140 405.7495 -55280.3682 206 207 208 209 210 -9228.8433 -14856.8323 -8692.5399 21996.6312 -11480.8201 211 212 213 214 215 44563.9461 44642.3210 70918.7810 10603.4828 17599.4573 216 217 218 219 220 44955.8202 48145.1926 47862.5583 45190.9447 17392.7647 221 222 223 224 225 20153.9316 45426.0691 46292.9672 37483.7138 45968.4724 226 227 228 229 230 46965.4284 49215.5708 45896.3180 30566.5692 46053.0676 231 232 233 234 235 18018.9923 34530.2267 46288.1921 28814.3868 46444.9417 236 237 238 239 240 46960.1975 46601.6913 48785.8793 46758.4409 46836.8158 241 242 243 244 245 46915.1906 41818.4282 47071.9402 51803.6457 -33513.3776 246 247 248 249 250 -5714.2970 -5635.9222 39918.7656 -5479.1725 -5400.7977 251 252 253 254 255 40752.2201 -5244.0481 29206.7092 27332.3264 -83022.1079 256 257 258 259 260 48090.8128 45767.9891 59514.2771 -4695.4244 48404.3120 261 262 263 264 265 -4538.6748 -4460.2999 53371.2611 66045.5447 -4225.1755 266 267 268 269 270 48874.5609 -4068.4259 82857.3389 -3911.6762 46988.9064 271 272 273 274 275 -3754.9266 -3676.5518 -3598.1770 45635.0936 49579.9342 276 277 278 279 280 55660.1016 52669.9998 49815.0587 -27957.7533 38477.7869 281 282 283 284 285 -2971.1785 50128.5579 79323.0019 50285.3076 43451.9650 286 287 288 289 290 50442.0572 -13990.7895 51714.8857 -2344.1800 59233.0149 291 292 293 294 295 -2187.4303 62255.4666 59639.3015 51069.0557 44872.3750 296 297 298 299 300 23965.2363 51304.1802 -1638.8066 51460.9298 51539.3046 301 302 303 304 305 -2930.0513 33247.9534 51774.4290 -1168.5578 90952.3650 306 307 308 309 310 55822.1264 -933.4333 -6260.6167 55849.5040 -698.3089 311 312 313 314 315 -619.9341 60022.2974 -463.1844 2035.8088 54553.1757 316 317 318 319 320 69254.8372 53278.9571 30262.6125 55582.9558 53106.8009 321 322 323 324 325 53185.1757 53263.5505 53341.9253 398.9385 53498.6750 326 327 328 329 330 50903.2519 634.0630 712.4378 53812.1742 53890.5490 331 332 333 334 335 53968.9238 54047.2987 54125.6735 52477.5435 55111.2563 336 337 338 339 340 47153.0532 43345.8869 54517.5475 52429.4205 60691.8419 341 342 343 344 345 28236.5456 54831.0468 54909.4216 54987.7964 50159.7441 346 347 348 349 350 55144.5460 45508.0229 55301.2957 55379.6705 19254.3082 351 352 353 354 355 59492.5143 56070.8167 55767.3742 48970.2762 55849.9194 356 357 358 359 360 55928.2942 36949.9413 56085.0438 56163.4186 56241.7934 361 362 363 364 365 58031.6874 56398.5431 56476.9179 -47304.2524 59018.5915 366 367 368 369 370 56712.0423 56790.4171 56868.7920 56947.1668 59768.6666 371 372 373 374 375 57103.9164 70807.9878 46377.8269 37095.0566 57417.4157 376 377 378 379 380 57495.7905 52666.8413 51794.5121 2339.7637 57809.2897 381 382 383 384 385 57887.6645 46625.4127 75947.9325 58122.7890 61096.1226 386 387 388 389 390 -88897.9241 70681.5594 31901.0989 57085.9673 44937.4683 391 392 393 394 395 54448.0084 28112.8887 -9212.4280 80465.3240 -164475.8443 396 397 398 399 400 20295.3487 -126452.0213 26193.5371 22511.8604 -9189.2937 401 402 403 404 405 1051.0932 51968.6029 29208.1466 -80264.3333 143.0929 406 407 408 409 410 23533.3418 36605.7253 27744.3089 38904.9583 -60715.3444 411 412 413 414 415 24985.3375 -234351.8343 -59353.2269 29258.0085 -174529.2171 416 417 418 419 420 -101353.6737 130756.7812 -46498.0682 -21419.1056 -35853.6695 421 422 423 424 425 -46844.9274 -91192.4427 -565498.6666 -470534.4341 -28246.5365 426 427 428 429 430 -213608.3783 -299955.6190 -1090321.6984 -289716.9918 -326935.9312 431 -660177.0992 > postscript(file="/var/www/html/freestat/rcomp/tmp/684kt1291223383.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 1284205.5414 NA 1 3187769.1732 1284205.5414 2 1190500.6719 3187769.1732 3 -2175874.6365 1190500.6719 4 -468857.7504 -2175874.6365 5 -90790.2427 -468857.7504 6 -26679.9252 -90790.2427 7 72632.9742 -26679.9252 8 -283972.6976 72632.9742 9 689755.2837 -283972.6976 10 -273154.5235 689755.2837 11 391457.0462 -273154.5235 12 746948.1834 391457.0462 13 -1198904.2326 746948.1834 14 4923.7717 -1198904.2326 15 491349.8015 4923.7717 16 -230376.0355 491349.8015 17 -10589.6859 -230376.0355 18 -317846.2350 -10589.6859 19 -325510.8117 -317846.2350 20 143082.0979 -325510.8117 21 -290933.1323 143082.0979 22 -128896.5763 -290933.1323 23 -384255.3887 -128896.5763 24 -289754.1051 -384255.3887 25 -149268.0248 -289754.1051 26 -409198.1778 -149268.0248 27 -297881.3181 -409198.1778 28 314718.7937 -297881.3181 29 -98990.2136 314718.7937 30 213934.4307 -98990.2136 31 -162462.4642 213934.4307 32 178442.5490 -162462.4642 33 70473.2483 178442.5490 34 -139861.8286 70473.2483 35 -469274.7105 -139861.8286 36 36163.3018 -469274.7105 37 140904.3362 36163.3018 38 -63919.8915 140904.3362 39 -170551.2589 -63919.8915 40 5074.2454 -170551.2589 41 -386406.3673 5074.2454 42 119682.3108 -386406.3673 43 -167669.1703 119682.3108 44 -23009.4475 -167669.1703 45 -134082.1738 -23009.4475 46 -293289.1576 -134082.1738 47 -130071.2149 -293289.1576 48 -36657.4800 -130071.2149 49 -171578.8841 -36657.4800 50 59693.6577 -171578.8841 51 -12903.4165 59693.6577 52 -247866.8332 -12903.4165 53 69915.9152 -247866.8332 54 -349096.0264 69915.9152 55 92688.0486 -349096.0264 56 281287.5672 92688.0486 57 -17778.3211 281287.5672 58 -4730.0110 -17778.3211 59 -231383.0298 -4730.0110 60 -49544.9854 -231383.0298 61 -164405.2161 -49544.9854 62 -61722.1836 -164405.2161 63 -198703.3558 -61722.1836 64 30439.4419 -198703.3558 65 34164.2920 30439.4419 66 -52044.9352 34164.2920 67 -49159.9099 -52044.9352 68 -38217.7441 -49159.9099 69 -75743.8977 -38217.7441 70 -137098.0592 -75743.8977 71 5261.7630 -137098.0592 72 -83664.6369 5261.7630 73 -1691.8845 -83664.6369 74 -37590.3118 -1691.8845 75 -162383.4934 -37590.3118 76 -22428.0277 -162383.4934 77 -74484.0321 -22428.0277 78 27515.6407 -74484.0321 79 -45306.8495 27515.6407 80 35998.9557 -45306.8495 81 -68712.6880 35998.9557 82 -72428.5580 -68712.6880 83 43390.1431 -72428.5580 84 10792.7087 43390.1431 85 -121329.8845 10792.7087 86 52116.0046 -121329.8845 87 -15736.2398 52116.0046 88 -90993.6726 -15736.2398 89 15434.7945 -90993.6726 90 11934.0118 15434.7945 91 -65600.8623 11934.0118 92 142726.8062 -65600.8623 93 -252110.6991 142726.8062 94 84496.8895 -252110.6991 95 -27816.9201 84496.8895 96 48646.1016 -27816.9201 97 23808.4614 48646.1016 98 -70514.9164 23808.4614 99 -74082.9420 -70514.9164 100 -76550.7077 -74082.9420 101 -118139.3378 -76550.7077 102 -133571.9433 -118139.3378 103 -71993.9074 -133571.9433 104 41430.7063 -71993.9074 105 -41607.6573 41430.7063 106 -274184.3298 -41607.6573 107 58310.5508 -274184.3298 108 -175406.9829 58310.5508 109 -146369.5250 -175406.9829 110 102020.2593 -146369.5250 111 -231903.5840 102020.2593 112 -155256.9593 -231903.5840 113 -131755.8649 -155256.9593 114 15041.9739 -131755.8649 115 45256.5727 15041.9739 116 9205.4005 45256.5727 117 55904.1942 9205.4005 118 -206784.6764 55904.1942 119 41591.1305 -206784.6764 120 -269345.1492 41591.1305 121 24916.8122 -269345.1492 122 17612.2972 24916.8122 123 26700.1079 17612.2972 124 -44470.8119 26700.1079 125 -229310.8781 -44470.8119 126 37980.4618 -229310.8781 127 37298.6434 37980.4618 128 45444.0513 37298.6434 129 38215.5862 45444.0513 130 39873.1091 38215.5862 131 82551.4642 39873.1091 132 38450.7106 82551.4642 133 35702.0045 38450.7106 134 18239.6111 35702.0045 135 33178.5791 18239.6111 136 44059.3673 33178.5791 137 39064.0778 44059.3673 138 13881.2785 39064.0778 139 38999.3343 13881.2785 140 46834.7729 38999.3343 141 39156.0840 46834.7729 142 6333.6478 39156.0840 143 39669.4956 6333.6478 144 41737.3438 39669.4956 145 40198.9096 41737.3438 146 39547.9580 40198.9096 147 41344.2876 39547.9580 148 42361.1560 41344.2876 149 39750.6735 42361.1560 150 8756.4069 39750.6735 151 9926.2905 8756.4069 152 40018.2069 9926.2905 153 17431.8785 40018.2069 154 40338.5765 17431.8785 155 40253.3314 40338.5765 156 40331.7062 40253.3314 157 36794.0965 40331.7062 158 52256.6186 36794.0965 159 40566.8306 52256.6186 160 40645.2054 40566.8306 161 40723.5803 40645.2054 162 40801.9551 40723.5803 163 28829.5624 40801.9551 164 31585.8041 28829.5624 165 40433.6828 31585.8041 166 41115.4543 40433.6828 167 121791.5073 41115.4543 168 98416.2146 121791.5073 169 41350.5788 98416.2146 170 41428.9536 41350.5788 171 37642.5124 41428.9536 172 41585.7032 37642.5124 173 31477.9707 41585.7032 174 41742.4528 31477.9707 175 42125.9025 41742.4528 176 93491.6577 42125.9025 177 -11977.8421 93491.6577 178 42055.9521 -11977.8421 179 42134.3269 42055.9521 180 51809.8333 42134.3269 181 -27618.0822 51809.8333 182 1604.0199 -27618.0822 183 -37849.6760 1604.0199 184 -12967.9289 -37849.6760 185 -4912.1108 -12967.9289 186 -37263.5792 -4912.1108 187 36004.1943 -37263.5792 188 5436.2354 36004.1943 189 42918.0750 5436.2354 190 42135.1276 42918.0750 191 43074.8247 42135.1276 192 39855.4995 43074.8247 193 43231.5743 39855.4995 194 43334.6202 43231.5743 195 -9633.0377 43334.6202 196 -121441.9611 -9633.0377 197 -9476.2880 -121441.9611 198 -9397.9132 -9476.2880 199 -96056.2200 -9397.9132 200 -9241.1636 -96056.2200 201 -9162.7888 -9241.1636 202 -9084.4140 -9162.7888 203 405.7495 -9084.4140 204 -55280.3682 405.7495 205 -9228.8433 -55280.3682 206 -14856.8323 -9228.8433 207 -8692.5399 -14856.8323 208 21996.6312 -8692.5399 209 -11480.8201 21996.6312 210 44563.9461 -11480.8201 211 44642.3210 44563.9461 212 70918.7810 44642.3210 213 10603.4828 70918.7810 214 17599.4573 10603.4828 215 44955.8202 17599.4573 216 48145.1926 44955.8202 217 47862.5583 48145.1926 218 45190.9447 47862.5583 219 17392.7647 45190.9447 220 20153.9316 17392.7647 221 45426.0691 20153.9316 222 46292.9672 45426.0691 223 37483.7138 46292.9672 224 45968.4724 37483.7138 225 46965.4284 45968.4724 226 49215.5708 46965.4284 227 45896.3180 49215.5708 228 30566.5692 45896.3180 229 46053.0676 30566.5692 230 18018.9923 46053.0676 231 34530.2267 18018.9923 232 46288.1921 34530.2267 233 28814.3868 46288.1921 234 46444.9417 28814.3868 235 46960.1975 46444.9417 236 46601.6913 46960.1975 237 48785.8793 46601.6913 238 46758.4409 48785.8793 239 46836.8158 46758.4409 240 46915.1906 46836.8158 241 41818.4282 46915.1906 242 47071.9402 41818.4282 243 51803.6457 47071.9402 244 -33513.3776 51803.6457 245 -5714.2970 -33513.3776 246 -5635.9222 -5714.2970 247 39918.7656 -5635.9222 248 -5479.1725 39918.7656 249 -5400.7977 -5479.1725 250 40752.2201 -5400.7977 251 -5244.0481 40752.2201 252 29206.7092 -5244.0481 253 27332.3264 29206.7092 254 -83022.1079 27332.3264 255 48090.8128 -83022.1079 256 45767.9891 48090.8128 257 59514.2771 45767.9891 258 -4695.4244 59514.2771 259 48404.3120 -4695.4244 260 -4538.6748 48404.3120 261 -4460.2999 -4538.6748 262 53371.2611 -4460.2999 263 66045.5447 53371.2611 264 -4225.1755 66045.5447 265 48874.5609 -4225.1755 266 -4068.4259 48874.5609 267 82857.3389 -4068.4259 268 -3911.6762 82857.3389 269 46988.9064 -3911.6762 270 -3754.9266 46988.9064 271 -3676.5518 -3754.9266 272 -3598.1770 -3676.5518 273 45635.0936 -3598.1770 274 49579.9342 45635.0936 275 55660.1016 49579.9342 276 52669.9998 55660.1016 277 49815.0587 52669.9998 278 -27957.7533 49815.0587 279 38477.7869 -27957.7533 280 -2971.1785 38477.7869 281 50128.5579 -2971.1785 282 79323.0019 50128.5579 283 50285.3076 79323.0019 284 43451.9650 50285.3076 285 50442.0572 43451.9650 286 -13990.7895 50442.0572 287 51714.8857 -13990.7895 288 -2344.1800 51714.8857 289 59233.0149 -2344.1800 290 -2187.4303 59233.0149 291 62255.4666 -2187.4303 292 59639.3015 62255.4666 293 51069.0557 59639.3015 294 44872.3750 51069.0557 295 23965.2363 44872.3750 296 51304.1802 23965.2363 297 -1638.8066 51304.1802 298 51460.9298 -1638.8066 299 51539.3046 51460.9298 300 -2930.0513 51539.3046 301 33247.9534 -2930.0513 302 51774.4290 33247.9534 303 -1168.5578 51774.4290 304 90952.3650 -1168.5578 305 55822.1264 90952.3650 306 -933.4333 55822.1264 307 -6260.6167 -933.4333 308 55849.5040 -6260.6167 309 -698.3089 55849.5040 310 -619.9341 -698.3089 311 60022.2974 -619.9341 312 -463.1844 60022.2974 313 2035.8088 -463.1844 314 54553.1757 2035.8088 315 69254.8372 54553.1757 316 53278.9571 69254.8372 317 30262.6125 53278.9571 318 55582.9558 30262.6125 319 53106.8009 55582.9558 320 53185.1757 53106.8009 321 53263.5505 53185.1757 322 53341.9253 53263.5505 323 398.9385 53341.9253 324 53498.6750 398.9385 325 50903.2519 53498.6750 326 634.0630 50903.2519 327 712.4378 634.0630 328 53812.1742 712.4378 329 53890.5490 53812.1742 330 53968.9238 53890.5490 331 54047.2987 53968.9238 332 54125.6735 54047.2987 333 52477.5435 54125.6735 334 55111.2563 52477.5435 335 47153.0532 55111.2563 336 43345.8869 47153.0532 337 54517.5475 43345.8869 338 52429.4205 54517.5475 339 60691.8419 52429.4205 340 28236.5456 60691.8419 341 54831.0468 28236.5456 342 54909.4216 54831.0468 343 54987.7964 54909.4216 344 50159.7441 54987.7964 345 55144.5460 50159.7441 346 45508.0229 55144.5460 347 55301.2957 45508.0229 348 55379.6705 55301.2957 349 19254.3082 55379.6705 350 59492.5143 19254.3082 351 56070.8167 59492.5143 352 55767.3742 56070.8167 353 48970.2762 55767.3742 354 55849.9194 48970.2762 355 55928.2942 55849.9194 356 36949.9413 55928.2942 357 56085.0438 36949.9413 358 56163.4186 56085.0438 359 56241.7934 56163.4186 360 58031.6874 56241.7934 361 56398.5431 58031.6874 362 56476.9179 56398.5431 363 -47304.2524 56476.9179 364 59018.5915 -47304.2524 365 56712.0423 59018.5915 366 56790.4171 56712.0423 367 56868.7920 56790.4171 368 56947.1668 56868.7920 369 59768.6666 56947.1668 370 57103.9164 59768.6666 371 70807.9878 57103.9164 372 46377.8269 70807.9878 373 37095.0566 46377.8269 374 57417.4157 37095.0566 375 57495.7905 57417.4157 376 52666.8413 57495.7905 377 51794.5121 52666.8413 378 2339.7637 51794.5121 379 57809.2897 2339.7637 380 57887.6645 57809.2897 381 46625.4127 57887.6645 382 75947.9325 46625.4127 383 58122.7890 75947.9325 384 61096.1226 58122.7890 385 -88897.9241 61096.1226 386 70681.5594 -88897.9241 387 31901.0989 70681.5594 388 57085.9673 31901.0989 389 44937.4683 57085.9673 390 54448.0084 44937.4683 391 28112.8887 54448.0084 392 -9212.4280 28112.8887 393 80465.3240 -9212.4280 394 -164475.8443 80465.3240 395 20295.3487 -164475.8443 396 -126452.0213 20295.3487 397 26193.5371 -126452.0213 398 22511.8604 26193.5371 399 -9189.2937 22511.8604 400 1051.0932 -9189.2937 401 51968.6029 1051.0932 402 29208.1466 51968.6029 403 -80264.3333 29208.1466 404 143.0929 -80264.3333 405 23533.3418 143.0929 406 36605.7253 23533.3418 407 27744.3089 36605.7253 408 38904.9583 27744.3089 409 -60715.3444 38904.9583 410 24985.3375 -60715.3444 411 -234351.8343 24985.3375 412 -59353.2269 -234351.8343 413 29258.0085 -59353.2269 414 -174529.2171 29258.0085 415 -101353.6737 -174529.2171 416 130756.7812 -101353.6737 417 -46498.0682 130756.7812 418 -21419.1056 -46498.0682 419 -35853.6695 -21419.1056 420 -46844.9274 -35853.6695 421 -91192.4427 -46844.9274 422 -565498.6666 -91192.4427 423 -470534.4341 -565498.6666 424 -28246.5365 -470534.4341 425 -213608.3783 -28246.5365 426 -299955.6190 -213608.3783 427 -1090321.6984 -299955.6190 428 -289716.9918 -1090321.6984 429 -326935.9312 -289716.9918 430 -660177.0992 -326935.9312 431 NA -660177.0992 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 3187769.1732 1284205.5414 [2,] 1190500.6719 3187769.1732 [3,] -2175874.6365 1190500.6719 [4,] -468857.7504 -2175874.6365 [5,] -90790.2427 -468857.7504 [6,] -26679.9252 -90790.2427 [7,] 72632.9742 -26679.9252 [8,] -283972.6976 72632.9742 [9,] 689755.2837 -283972.6976 [10,] -273154.5235 689755.2837 [11,] 391457.0462 -273154.5235 [12,] 746948.1834 391457.0462 [13,] -1198904.2326 746948.1834 [14,] 4923.7717 -1198904.2326 [15,] 491349.8015 4923.7717 [16,] -230376.0355 491349.8015 [17,] -10589.6859 -230376.0355 [18,] -317846.2350 -10589.6859 [19,] -325510.8117 -317846.2350 [20,] 143082.0979 -325510.8117 [21,] -290933.1323 143082.0979 [22,] -128896.5763 -290933.1323 [23,] -384255.3887 -128896.5763 [24,] -289754.1051 -384255.3887 [25,] -149268.0248 -289754.1051 [26,] -409198.1778 -149268.0248 [27,] -297881.3181 -409198.1778 [28,] 314718.7937 -297881.3181 [29,] -98990.2136 314718.7937 [30,] 213934.4307 -98990.2136 [31,] -162462.4642 213934.4307 [32,] 178442.5490 -162462.4642 [33,] 70473.2483 178442.5490 [34,] -139861.8286 70473.2483 [35,] -469274.7105 -139861.8286 [36,] 36163.3018 -469274.7105 [37,] 140904.3362 36163.3018 [38,] -63919.8915 140904.3362 [39,] -170551.2589 -63919.8915 [40,] 5074.2454 -170551.2589 [41,] -386406.3673 5074.2454 [42,] 119682.3108 -386406.3673 [43,] -167669.1703 119682.3108 [44,] -23009.4475 -167669.1703 [45,] -134082.1738 -23009.4475 [46,] -293289.1576 -134082.1738 [47,] -130071.2149 -293289.1576 [48,] -36657.4800 -130071.2149 [49,] -171578.8841 -36657.4800 [50,] 59693.6577 -171578.8841 [51,] -12903.4165 59693.6577 [52,] -247866.8332 -12903.4165 [53,] 69915.9152 -247866.8332 [54,] -349096.0264 69915.9152 [55,] 92688.0486 -349096.0264 [56,] 281287.5672 92688.0486 [57,] -17778.3211 281287.5672 [58,] -4730.0110 -17778.3211 [59,] -231383.0298 -4730.0110 [60,] -49544.9854 -231383.0298 [61,] -164405.2161 -49544.9854 [62,] -61722.1836 -164405.2161 [63,] -198703.3558 -61722.1836 [64,] 30439.4419 -198703.3558 [65,] 34164.2920 30439.4419 [66,] -52044.9352 34164.2920 [67,] -49159.9099 -52044.9352 [68,] -38217.7441 -49159.9099 [69,] -75743.8977 -38217.7441 [70,] -137098.0592 -75743.8977 [71,] 5261.7630 -137098.0592 [72,] -83664.6369 5261.7630 [73,] -1691.8845 -83664.6369 [74,] -37590.3118 -1691.8845 [75,] -162383.4934 -37590.3118 [76,] -22428.0277 -162383.4934 [77,] -74484.0321 -22428.0277 [78,] 27515.6407 -74484.0321 [79,] -45306.8495 27515.6407 [80,] 35998.9557 -45306.8495 [81,] -68712.6880 35998.9557 [82,] -72428.5580 -68712.6880 [83,] 43390.1431 -72428.5580 [84,] 10792.7087 43390.1431 [85,] -121329.8845 10792.7087 [86,] 52116.0046 -121329.8845 [87,] -15736.2398 52116.0046 [88,] -90993.6726 -15736.2398 [89,] 15434.7945 -90993.6726 [90,] 11934.0118 15434.7945 [91,] -65600.8623 11934.0118 [92,] 142726.8062 -65600.8623 [93,] -252110.6991 142726.8062 [94,] 84496.8895 -252110.6991 [95,] -27816.9201 84496.8895 [96,] 48646.1016 -27816.9201 [97,] 23808.4614 48646.1016 [98,] -70514.9164 23808.4614 [99,] -74082.9420 -70514.9164 [100,] -76550.7077 -74082.9420 [101,] -118139.3378 -76550.7077 [102,] -133571.9433 -118139.3378 [103,] -71993.9074 -133571.9433 [104,] 41430.7063 -71993.9074 [105,] -41607.6573 41430.7063 [106,] -274184.3298 -41607.6573 [107,] 58310.5508 -274184.3298 [108,] -175406.9829 58310.5508 [109,] -146369.5250 -175406.9829 [110,] 102020.2593 -146369.5250 [111,] -231903.5840 102020.2593 [112,] -155256.9593 -231903.5840 [113,] -131755.8649 -155256.9593 [114,] 15041.9739 -131755.8649 [115,] 45256.5727 15041.9739 [116,] 9205.4005 45256.5727 [117,] 55904.1942 9205.4005 [118,] -206784.6764 55904.1942 [119,] 41591.1305 -206784.6764 [120,] -269345.1492 41591.1305 [121,] 24916.8122 -269345.1492 [122,] 17612.2972 24916.8122 [123,] 26700.1079 17612.2972 [124,] -44470.8119 26700.1079 [125,] -229310.8781 -44470.8119 [126,] 37980.4618 -229310.8781 [127,] 37298.6434 37980.4618 [128,] 45444.0513 37298.6434 [129,] 38215.5862 45444.0513 [130,] 39873.1091 38215.5862 [131,] 82551.4642 39873.1091 [132,] 38450.7106 82551.4642 [133,] 35702.0045 38450.7106 [134,] 18239.6111 35702.0045 [135,] 33178.5791 18239.6111 [136,] 44059.3673 33178.5791 [137,] 39064.0778 44059.3673 [138,] 13881.2785 39064.0778 [139,] 38999.3343 13881.2785 [140,] 46834.7729 38999.3343 [141,] 39156.0840 46834.7729 [142,] 6333.6478 39156.0840 [143,] 39669.4956 6333.6478 [144,] 41737.3438 39669.4956 [145,] 40198.9096 41737.3438 [146,] 39547.9580 40198.9096 [147,] 41344.2876 39547.9580 [148,] 42361.1560 41344.2876 [149,] 39750.6735 42361.1560 [150,] 8756.4069 39750.6735 [151,] 9926.2905 8756.4069 [152,] 40018.2069 9926.2905 [153,] 17431.8785 40018.2069 [154,] 40338.5765 17431.8785 [155,] 40253.3314 40338.5765 [156,] 40331.7062 40253.3314 [157,] 36794.0965 40331.7062 [158,] 52256.6186 36794.0965 [159,] 40566.8306 52256.6186 [160,] 40645.2054 40566.8306 [161,] 40723.5803 40645.2054 [162,] 40801.9551 40723.5803 [163,] 28829.5624 40801.9551 [164,] 31585.8041 28829.5624 [165,] 40433.6828 31585.8041 [166,] 41115.4543 40433.6828 [167,] 121791.5073 41115.4543 [168,] 98416.2146 121791.5073 [169,] 41350.5788 98416.2146 [170,] 41428.9536 41350.5788 [171,] 37642.5124 41428.9536 [172,] 41585.7032 37642.5124 [173,] 31477.9707 41585.7032 [174,] 41742.4528 31477.9707 [175,] 42125.9025 41742.4528 [176,] 93491.6577 42125.9025 [177,] -11977.8421 93491.6577 [178,] 42055.9521 -11977.8421 [179,] 42134.3269 42055.9521 [180,] 51809.8333 42134.3269 [181,] -27618.0822 51809.8333 [182,] 1604.0199 -27618.0822 [183,] -37849.6760 1604.0199 [184,] -12967.9289 -37849.6760 [185,] -4912.1108 -12967.9289 [186,] -37263.5792 -4912.1108 [187,] 36004.1943 -37263.5792 [188,] 5436.2354 36004.1943 [189,] 42918.0750 5436.2354 [190,] 42135.1276 42918.0750 [191,] 43074.8247 42135.1276 [192,] 39855.4995 43074.8247 [193,] 43231.5743 39855.4995 [194,] 43334.6202 43231.5743 [195,] -9633.0377 43334.6202 [196,] -121441.9611 -9633.0377 [197,] -9476.2880 -121441.9611 [198,] -9397.9132 -9476.2880 [199,] -96056.2200 -9397.9132 [200,] -9241.1636 -96056.2200 [201,] -9162.7888 -9241.1636 [202,] -9084.4140 -9162.7888 [203,] 405.7495 -9084.4140 [204,] -55280.3682 405.7495 [205,] -9228.8433 -55280.3682 [206,] -14856.8323 -9228.8433 [207,] -8692.5399 -14856.8323 [208,] 21996.6312 -8692.5399 [209,] -11480.8201 21996.6312 [210,] 44563.9461 -11480.8201 [211,] 44642.3210 44563.9461 [212,] 70918.7810 44642.3210 [213,] 10603.4828 70918.7810 [214,] 17599.4573 10603.4828 [215,] 44955.8202 17599.4573 [216,] 48145.1926 44955.8202 [217,] 47862.5583 48145.1926 [218,] 45190.9447 47862.5583 [219,] 17392.7647 45190.9447 [220,] 20153.9316 17392.7647 [221,] 45426.0691 20153.9316 [222,] 46292.9672 45426.0691 [223,] 37483.7138 46292.9672 [224,] 45968.4724 37483.7138 [225,] 46965.4284 45968.4724 [226,] 49215.5708 46965.4284 [227,] 45896.3180 49215.5708 [228,] 30566.5692 45896.3180 [229,] 46053.0676 30566.5692 [230,] 18018.9923 46053.0676 [231,] 34530.2267 18018.9923 [232,] 46288.1921 34530.2267 [233,] 28814.3868 46288.1921 [234,] 46444.9417 28814.3868 [235,] 46960.1975 46444.9417 [236,] 46601.6913 46960.1975 [237,] 48785.8793 46601.6913 [238,] 46758.4409 48785.8793 [239,] 46836.8158 46758.4409 [240,] 46915.1906 46836.8158 [241,] 41818.4282 46915.1906 [242,] 47071.9402 41818.4282 [243,] 51803.6457 47071.9402 [244,] -33513.3776 51803.6457 [245,] -5714.2970 -33513.3776 [246,] -5635.9222 -5714.2970 [247,] 39918.7656 -5635.9222 [248,] -5479.1725 39918.7656 [249,] -5400.7977 -5479.1725 [250,] 40752.2201 -5400.7977 [251,] -5244.0481 40752.2201 [252,] 29206.7092 -5244.0481 [253,] 27332.3264 29206.7092 [254,] -83022.1079 27332.3264 [255,] 48090.8128 -83022.1079 [256,] 45767.9891 48090.8128 [257,] 59514.2771 45767.9891 [258,] -4695.4244 59514.2771 [259,] 48404.3120 -4695.4244 [260,] -4538.6748 48404.3120 [261,] -4460.2999 -4538.6748 [262,] 53371.2611 -4460.2999 [263,] 66045.5447 53371.2611 [264,] -4225.1755 66045.5447 [265,] 48874.5609 -4225.1755 [266,] -4068.4259 48874.5609 [267,] 82857.3389 -4068.4259 [268,] -3911.6762 82857.3389 [269,] 46988.9064 -3911.6762 [270,] -3754.9266 46988.9064 [271,] -3676.5518 -3754.9266 [272,] -3598.1770 -3676.5518 [273,] 45635.0936 -3598.1770 [274,] 49579.9342 45635.0936 [275,] 55660.1016 49579.9342 [276,] 52669.9998 55660.1016 [277,] 49815.0587 52669.9998 [278,] -27957.7533 49815.0587 [279,] 38477.7869 -27957.7533 [280,] -2971.1785 38477.7869 [281,] 50128.5579 -2971.1785 [282,] 79323.0019 50128.5579 [283,] 50285.3076 79323.0019 [284,] 43451.9650 50285.3076 [285,] 50442.0572 43451.9650 [286,] -13990.7895 50442.0572 [287,] 51714.8857 -13990.7895 [288,] -2344.1800 51714.8857 [289,] 59233.0149 -2344.1800 [290,] -2187.4303 59233.0149 [291,] 62255.4666 -2187.4303 [292,] 59639.3015 62255.4666 [293,] 51069.0557 59639.3015 [294,] 44872.3750 51069.0557 [295,] 23965.2363 44872.3750 [296,] 51304.1802 23965.2363 [297,] -1638.8066 51304.1802 [298,] 51460.9298 -1638.8066 [299,] 51539.3046 51460.9298 [300,] -2930.0513 51539.3046 [301,] 33247.9534 -2930.0513 [302,] 51774.4290 33247.9534 [303,] -1168.5578 51774.4290 [304,] 90952.3650 -1168.5578 [305,] 55822.1264 90952.3650 [306,] -933.4333 55822.1264 [307,] -6260.6167 -933.4333 [308,] 55849.5040 -6260.6167 [309,] -698.3089 55849.5040 [310,] -619.9341 -698.3089 [311,] 60022.2974 -619.9341 [312,] -463.1844 60022.2974 [313,] 2035.8088 -463.1844 [314,] 54553.1757 2035.8088 [315,] 69254.8372 54553.1757 [316,] 53278.9571 69254.8372 [317,] 30262.6125 53278.9571 [318,] 55582.9558 30262.6125 [319,] 53106.8009 55582.9558 [320,] 53185.1757 53106.8009 [321,] 53263.5505 53185.1757 [322,] 53341.9253 53263.5505 [323,] 398.9385 53341.9253 [324,] 53498.6750 398.9385 [325,] 50903.2519 53498.6750 [326,] 634.0630 50903.2519 [327,] 712.4378 634.0630 [328,] 53812.1742 712.4378 [329,] 53890.5490 53812.1742 [330,] 53968.9238 53890.5490 [331,] 54047.2987 53968.9238 [332,] 54125.6735 54047.2987 [333,] 52477.5435 54125.6735 [334,] 55111.2563 52477.5435 [335,] 47153.0532 55111.2563 [336,] 43345.8869 47153.0532 [337,] 54517.5475 43345.8869 [338,] 52429.4205 54517.5475 [339,] 60691.8419 52429.4205 [340,] 28236.5456 60691.8419 [341,] 54831.0468 28236.5456 [342,] 54909.4216 54831.0468 [343,] 54987.7964 54909.4216 [344,] 50159.7441 54987.7964 [345,] 55144.5460 50159.7441 [346,] 45508.0229 55144.5460 [347,] 55301.2957 45508.0229 [348,] 55379.6705 55301.2957 [349,] 19254.3082 55379.6705 [350,] 59492.5143 19254.3082 [351,] 56070.8167 59492.5143 [352,] 55767.3742 56070.8167 [353,] 48970.2762 55767.3742 [354,] 55849.9194 48970.2762 [355,] 55928.2942 55849.9194 [356,] 36949.9413 55928.2942 [357,] 56085.0438 36949.9413 [358,] 56163.4186 56085.0438 [359,] 56241.7934 56163.4186 [360,] 58031.6874 56241.7934 [361,] 56398.5431 58031.6874 [362,] 56476.9179 56398.5431 [363,] -47304.2524 56476.9179 [364,] 59018.5915 -47304.2524 [365,] 56712.0423 59018.5915 [366,] 56790.4171 56712.0423 [367,] 56868.7920 56790.4171 [368,] 56947.1668 56868.7920 [369,] 59768.6666 56947.1668 [370,] 57103.9164 59768.6666 [371,] 70807.9878 57103.9164 [372,] 46377.8269 70807.9878 [373,] 37095.0566 46377.8269 [374,] 57417.4157 37095.0566 [375,] 57495.7905 57417.4157 [376,] 52666.8413 57495.7905 [377,] 51794.5121 52666.8413 [378,] 2339.7637 51794.5121 [379,] 57809.2897 2339.7637 [380,] 57887.6645 57809.2897 [381,] 46625.4127 57887.6645 [382,] 75947.9325 46625.4127 [383,] 58122.7890 75947.9325 [384,] 61096.1226 58122.7890 [385,] -88897.9241 61096.1226 [386,] 70681.5594 -88897.9241 [387,] 31901.0989 70681.5594 [388,] 57085.9673 31901.0989 [389,] 44937.4683 57085.9673 [390,] 54448.0084 44937.4683 [391,] 28112.8887 54448.0084 [392,] -9212.4280 28112.8887 [393,] 80465.3240 -9212.4280 [394,] -164475.8443 80465.3240 [395,] 20295.3487 -164475.8443 [396,] -126452.0213 20295.3487 [397,] 26193.5371 -126452.0213 [398,] 22511.8604 26193.5371 [399,] -9189.2937 22511.8604 [400,] 1051.0932 -9189.2937 [401,] 51968.6029 1051.0932 [402,] 29208.1466 51968.6029 [403,] -80264.3333 29208.1466 [404,] 143.0929 -80264.3333 [405,] 23533.3418 143.0929 [406,] 36605.7253 23533.3418 [407,] 27744.3089 36605.7253 [408,] 38904.9583 27744.3089 [409,] -60715.3444 38904.9583 [410,] 24985.3375 -60715.3444 [411,] -234351.8343 24985.3375 [412,] -59353.2269 -234351.8343 [413,] 29258.0085 -59353.2269 [414,] -174529.2171 29258.0085 [415,] -101353.6737 -174529.2171 [416,] 130756.7812 -101353.6737 [417,] -46498.0682 130756.7812 [418,] -21419.1056 -46498.0682 [419,] -35853.6695 -21419.1056 [420,] -46844.9274 -35853.6695 [421,] -91192.4427 -46844.9274 [422,] -565498.6666 -91192.4427 [423,] -470534.4341 -565498.6666 [424,] -28246.5365 -470534.4341 [425,] -213608.3783 -28246.5365 [426,] -299955.6190 -213608.3783 [427,] -1090321.6984 -299955.6190 [428,] -289716.9918 -1090321.6984 [429,] -326935.9312 -289716.9918 [430,] -660177.0992 -326935.9312 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 3187769.1732 1284205.5414 2 1190500.6719 3187769.1732 3 -2175874.6365 1190500.6719 4 -468857.7504 -2175874.6365 5 -90790.2427 -468857.7504 6 -26679.9252 -90790.2427 7 72632.9742 -26679.9252 8 -283972.6976 72632.9742 9 689755.2837 -283972.6976 10 -273154.5235 689755.2837 11 391457.0462 -273154.5235 12 746948.1834 391457.0462 13 -1198904.2326 746948.1834 14 4923.7717 -1198904.2326 15 491349.8015 4923.7717 16 -230376.0355 491349.8015 17 -10589.6859 -230376.0355 18 -317846.2350 -10589.6859 19 -325510.8117 -317846.2350 20 143082.0979 -325510.8117 21 -290933.1323 143082.0979 22 -128896.5763 -290933.1323 23 -384255.3887 -128896.5763 24 -289754.1051 -384255.3887 25 -149268.0248 -289754.1051 26 -409198.1778 -149268.0248 27 -297881.3181 -409198.1778 28 314718.7937 -297881.3181 29 -98990.2136 314718.7937 30 213934.4307 -98990.2136 31 -162462.4642 213934.4307 32 178442.5490 -162462.4642 33 70473.2483 178442.5490 34 -139861.8286 70473.2483 35 -469274.7105 -139861.8286 36 36163.3018 -469274.7105 37 140904.3362 36163.3018 38 -63919.8915 140904.3362 39 -170551.2589 -63919.8915 40 5074.2454 -170551.2589 41 -386406.3673 5074.2454 42 119682.3108 -386406.3673 43 -167669.1703 119682.3108 44 -23009.4475 -167669.1703 45 -134082.1738 -23009.4475 46 -293289.1576 -134082.1738 47 -130071.2149 -293289.1576 48 -36657.4800 -130071.2149 49 -171578.8841 -36657.4800 50 59693.6577 -171578.8841 51 -12903.4165 59693.6577 52 -247866.8332 -12903.4165 53 69915.9152 -247866.8332 54 -349096.0264 69915.9152 55 92688.0486 -349096.0264 56 281287.5672 92688.0486 57 -17778.3211 281287.5672 58 -4730.0110 -17778.3211 59 -231383.0298 -4730.0110 60 -49544.9854 -231383.0298 61 -164405.2161 -49544.9854 62 -61722.1836 -164405.2161 63 -198703.3558 -61722.1836 64 30439.4419 -198703.3558 65 34164.2920 30439.4419 66 -52044.9352 34164.2920 67 -49159.9099 -52044.9352 68 -38217.7441 -49159.9099 69 -75743.8977 -38217.7441 70 -137098.0592 -75743.8977 71 5261.7630 -137098.0592 72 -83664.6369 5261.7630 73 -1691.8845 -83664.6369 74 -37590.3118 -1691.8845 75 -162383.4934 -37590.3118 76 -22428.0277 -162383.4934 77 -74484.0321 -22428.0277 78 27515.6407 -74484.0321 79 -45306.8495 27515.6407 80 35998.9557 -45306.8495 81 -68712.6880 35998.9557 82 -72428.5580 -68712.6880 83 43390.1431 -72428.5580 84 10792.7087 43390.1431 85 -121329.8845 10792.7087 86 52116.0046 -121329.8845 87 -15736.2398 52116.0046 88 -90993.6726 -15736.2398 89 15434.7945 -90993.6726 90 11934.0118 15434.7945 91 -65600.8623 11934.0118 92 142726.8062 -65600.8623 93 -252110.6991 142726.8062 94 84496.8895 -252110.6991 95 -27816.9201 84496.8895 96 48646.1016 -27816.9201 97 23808.4614 48646.1016 98 -70514.9164 23808.4614 99 -74082.9420 -70514.9164 100 -76550.7077 -74082.9420 101 -118139.3378 -76550.7077 102 -133571.9433 -118139.3378 103 -71993.9074 -133571.9433 104 41430.7063 -71993.9074 105 -41607.6573 41430.7063 106 -274184.3298 -41607.6573 107 58310.5508 -274184.3298 108 -175406.9829 58310.5508 109 -146369.5250 -175406.9829 110 102020.2593 -146369.5250 111 -231903.5840 102020.2593 112 -155256.9593 -231903.5840 113 -131755.8649 -155256.9593 114 15041.9739 -131755.8649 115 45256.5727 15041.9739 116 9205.4005 45256.5727 117 55904.1942 9205.4005 118 -206784.6764 55904.1942 119 41591.1305 -206784.6764 120 -269345.1492 41591.1305 121 24916.8122 -269345.1492 122 17612.2972 24916.8122 123 26700.1079 17612.2972 124 -44470.8119 26700.1079 125 -229310.8781 -44470.8119 126 37980.4618 -229310.8781 127 37298.6434 37980.4618 128 45444.0513 37298.6434 129 38215.5862 45444.0513 130 39873.1091 38215.5862 131 82551.4642 39873.1091 132 38450.7106 82551.4642 133 35702.0045 38450.7106 134 18239.6111 35702.0045 135 33178.5791 18239.6111 136 44059.3673 33178.5791 137 39064.0778 44059.3673 138 13881.2785 39064.0778 139 38999.3343 13881.2785 140 46834.7729 38999.3343 141 39156.0840 46834.7729 142 6333.6478 39156.0840 143 39669.4956 6333.6478 144 41737.3438 39669.4956 145 40198.9096 41737.3438 146 39547.9580 40198.9096 147 41344.2876 39547.9580 148 42361.1560 41344.2876 149 39750.6735 42361.1560 150 8756.4069 39750.6735 151 9926.2905 8756.4069 152 40018.2069 9926.2905 153 17431.8785 40018.2069 154 40338.5765 17431.8785 155 40253.3314 40338.5765 156 40331.7062 40253.3314 157 36794.0965 40331.7062 158 52256.6186 36794.0965 159 40566.8306 52256.6186 160 40645.2054 40566.8306 161 40723.5803 40645.2054 162 40801.9551 40723.5803 163 28829.5624 40801.9551 164 31585.8041 28829.5624 165 40433.6828 31585.8041 166 41115.4543 40433.6828 167 121791.5073 41115.4543 168 98416.2146 121791.5073 169 41350.5788 98416.2146 170 41428.9536 41350.5788 171 37642.5124 41428.9536 172 41585.7032 37642.5124 173 31477.9707 41585.7032 174 41742.4528 31477.9707 175 42125.9025 41742.4528 176 93491.6577 42125.9025 177 -11977.8421 93491.6577 178 42055.9521 -11977.8421 179 42134.3269 42055.9521 180 51809.8333 42134.3269 181 -27618.0822 51809.8333 182 1604.0199 -27618.0822 183 -37849.6760 1604.0199 184 -12967.9289 -37849.6760 185 -4912.1108 -12967.9289 186 -37263.5792 -4912.1108 187 36004.1943 -37263.5792 188 5436.2354 36004.1943 189 42918.0750 5436.2354 190 42135.1276 42918.0750 191 43074.8247 42135.1276 192 39855.4995 43074.8247 193 43231.5743 39855.4995 194 43334.6202 43231.5743 195 -9633.0377 43334.6202 196 -121441.9611 -9633.0377 197 -9476.2880 -121441.9611 198 -9397.9132 -9476.2880 199 -96056.2200 -9397.9132 200 -9241.1636 -96056.2200 201 -9162.7888 -9241.1636 202 -9084.4140 -9162.7888 203 405.7495 -9084.4140 204 -55280.3682 405.7495 205 -9228.8433 -55280.3682 206 -14856.8323 -9228.8433 207 -8692.5399 -14856.8323 208 21996.6312 -8692.5399 209 -11480.8201 21996.6312 210 44563.9461 -11480.8201 211 44642.3210 44563.9461 212 70918.7810 44642.3210 213 10603.4828 70918.7810 214 17599.4573 10603.4828 215 44955.8202 17599.4573 216 48145.1926 44955.8202 217 47862.5583 48145.1926 218 45190.9447 47862.5583 219 17392.7647 45190.9447 220 20153.9316 17392.7647 221 45426.0691 20153.9316 222 46292.9672 45426.0691 223 37483.7138 46292.9672 224 45968.4724 37483.7138 225 46965.4284 45968.4724 226 49215.5708 46965.4284 227 45896.3180 49215.5708 228 30566.5692 45896.3180 229 46053.0676 30566.5692 230 18018.9923 46053.0676 231 34530.2267 18018.9923 232 46288.1921 34530.2267 233 28814.3868 46288.1921 234 46444.9417 28814.3868 235 46960.1975 46444.9417 236 46601.6913 46960.1975 237 48785.8793 46601.6913 238 46758.4409 48785.8793 239 46836.8158 46758.4409 240 46915.1906 46836.8158 241 41818.4282 46915.1906 242 47071.9402 41818.4282 243 51803.6457 47071.9402 244 -33513.3776 51803.6457 245 -5714.2970 -33513.3776 246 -5635.9222 -5714.2970 247 39918.7656 -5635.9222 248 -5479.1725 39918.7656 249 -5400.7977 -5479.1725 250 40752.2201 -5400.7977 251 -5244.0481 40752.2201 252 29206.7092 -5244.0481 253 27332.3264 29206.7092 254 -83022.1079 27332.3264 255 48090.8128 -83022.1079 256 45767.9891 48090.8128 257 59514.2771 45767.9891 258 -4695.4244 59514.2771 259 48404.3120 -4695.4244 260 -4538.6748 48404.3120 261 -4460.2999 -4538.6748 262 53371.2611 -4460.2999 263 66045.5447 53371.2611 264 -4225.1755 66045.5447 265 48874.5609 -4225.1755 266 -4068.4259 48874.5609 267 82857.3389 -4068.4259 268 -3911.6762 82857.3389 269 46988.9064 -3911.6762 270 -3754.9266 46988.9064 271 -3676.5518 -3754.9266 272 -3598.1770 -3676.5518 273 45635.0936 -3598.1770 274 49579.9342 45635.0936 275 55660.1016 49579.9342 276 52669.9998 55660.1016 277 49815.0587 52669.9998 278 -27957.7533 49815.0587 279 38477.7869 -27957.7533 280 -2971.1785 38477.7869 281 50128.5579 -2971.1785 282 79323.0019 50128.5579 283 50285.3076 79323.0019 284 43451.9650 50285.3076 285 50442.0572 43451.9650 286 -13990.7895 50442.0572 287 51714.8857 -13990.7895 288 -2344.1800 51714.8857 289 59233.0149 -2344.1800 290 -2187.4303 59233.0149 291 62255.4666 -2187.4303 292 59639.3015 62255.4666 293 51069.0557 59639.3015 294 44872.3750 51069.0557 295 23965.2363 44872.3750 296 51304.1802 23965.2363 297 -1638.8066 51304.1802 298 51460.9298 -1638.8066 299 51539.3046 51460.9298 300 -2930.0513 51539.3046 301 33247.9534 -2930.0513 302 51774.4290 33247.9534 303 -1168.5578 51774.4290 304 90952.3650 -1168.5578 305 55822.1264 90952.3650 306 -933.4333 55822.1264 307 -6260.6167 -933.4333 308 55849.5040 -6260.6167 309 -698.3089 55849.5040 310 -619.9341 -698.3089 311 60022.2974 -619.9341 312 -463.1844 60022.2974 313 2035.8088 -463.1844 314 54553.1757 2035.8088 315 69254.8372 54553.1757 316 53278.9571 69254.8372 317 30262.6125 53278.9571 318 55582.9558 30262.6125 319 53106.8009 55582.9558 320 53185.1757 53106.8009 321 53263.5505 53185.1757 322 53341.9253 53263.5505 323 398.9385 53341.9253 324 53498.6750 398.9385 325 50903.2519 53498.6750 326 634.0630 50903.2519 327 712.4378 634.0630 328 53812.1742 712.4378 329 53890.5490 53812.1742 330 53968.9238 53890.5490 331 54047.2987 53968.9238 332 54125.6735 54047.2987 333 52477.5435 54125.6735 334 55111.2563 52477.5435 335 47153.0532 55111.2563 336 43345.8869 47153.0532 337 54517.5475 43345.8869 338 52429.4205 54517.5475 339 60691.8419 52429.4205 340 28236.5456 60691.8419 341 54831.0468 28236.5456 342 54909.4216 54831.0468 343 54987.7964 54909.4216 344 50159.7441 54987.7964 345 55144.5460 50159.7441 346 45508.0229 55144.5460 347 55301.2957 45508.0229 348 55379.6705 55301.2957 349 19254.3082 55379.6705 350 59492.5143 19254.3082 351 56070.8167 59492.5143 352 55767.3742 56070.8167 353 48970.2762 55767.3742 354 55849.9194 48970.2762 355 55928.2942 55849.9194 356 36949.9413 55928.2942 357 56085.0438 36949.9413 358 56163.4186 56085.0438 359 56241.7934 56163.4186 360 58031.6874 56241.7934 361 56398.5431 58031.6874 362 56476.9179 56398.5431 363 -47304.2524 56476.9179 364 59018.5915 -47304.2524 365 56712.0423 59018.5915 366 56790.4171 56712.0423 367 56868.7920 56790.4171 368 56947.1668 56868.7920 369 59768.6666 56947.1668 370 57103.9164 59768.6666 371 70807.9878 57103.9164 372 46377.8269 70807.9878 373 37095.0566 46377.8269 374 57417.4157 37095.0566 375 57495.7905 57417.4157 376 52666.8413 57495.7905 377 51794.5121 52666.8413 378 2339.7637 51794.5121 379 57809.2897 2339.7637 380 57887.6645 57809.2897 381 46625.4127 57887.6645 382 75947.9325 46625.4127 383 58122.7890 75947.9325 384 61096.1226 58122.7890 385 -88897.9241 61096.1226 386 70681.5594 -88897.9241 387 31901.0989 70681.5594 388 57085.9673 31901.0989 389 44937.4683 57085.9673 390 54448.0084 44937.4683 391 28112.8887 54448.0084 392 -9212.4280 28112.8887 393 80465.3240 -9212.4280 394 -164475.8443 80465.3240 395 20295.3487 -164475.8443 396 -126452.0213 20295.3487 397 26193.5371 -126452.0213 398 22511.8604 26193.5371 399 -9189.2937 22511.8604 400 1051.0932 -9189.2937 401 51968.6029 1051.0932 402 29208.1466 51968.6029 403 -80264.3333 29208.1466 404 143.0929 -80264.3333 405 23533.3418 143.0929 406 36605.7253 23533.3418 407 27744.3089 36605.7253 408 38904.9583 27744.3089 409 -60715.3444 38904.9583 410 24985.3375 -60715.3444 411 -234351.8343 24985.3375 412 -59353.2269 -234351.8343 413 29258.0085 -59353.2269 414 -174529.2171 29258.0085 415 -101353.6737 -174529.2171 416 130756.7812 -101353.6737 417 -46498.0682 130756.7812 418 -21419.1056 -46498.0682 419 -35853.6695 -21419.1056 420 -46844.9274 -35853.6695 421 -91192.4427 -46844.9274 422 -565498.6666 -91192.4427 423 -470534.4341 -565498.6666 424 -28246.5365 -470534.4341 425 -213608.3783 -28246.5365 426 -299955.6190 -213608.3783 427 -1090321.6984 -299955.6190 428 -289716.9918 -1090321.6984 429 -326935.9312 -289716.9918 430 -660177.0992 -326935.9312 > 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/71d1w1291223383.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/8b5iz1291223383.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/9b5iz1291223383.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/10mwh21291223383.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/11qxyq1291223383.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/12ioxb1291223383.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/13p7cn1291223383.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/140gtp1291223383.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/15lzsv1291223383.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/16h9q41291223383.tab") + } > > try(system("convert tmp/1fdkq1291223383.ps tmp/1fdkq1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/2fdkq1291223383.ps tmp/2fdkq1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/3fdkq1291223383.ps tmp/3fdkq1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/484kt1291223383.ps tmp/484kt1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/584kt1291223383.ps tmp/584kt1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/684kt1291223383.ps tmp/684kt1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/71d1w1291223383.ps tmp/71d1w1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/8b5iz1291223383.ps tmp/8b5iz1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/9b5iz1291223383.ps tmp/9b5iz1291223383.png",intern=TRUE)) character(0) > try(system("convert tmp/10mwh21291223383.ps tmp/10mwh21291223383.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 14.417 3.177 14.893