R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(162556 + ,1081 + ,213118 + ,6282154 + ,29790 + ,309 + ,81767 + ,4321023 + ,87550 + ,458 + ,153198 + ,4111912 + ,84738 + ,588 + ,-26007 + ,223193 + ,54660 + ,302 + ,126942 + ,1491348 + ,42634 + ,156 + ,157214 + ,1629616 + ,40949 + ,481 + ,129352 + ,1398893 + ,45187 + ,353 + ,234817 + ,1926517 + ,37704 + ,452 + ,60448 + ,983660 + ,16275 + ,109 + ,47818 + ,1443586 + ,25830 + ,115 + ,245546 + ,1073089 + ,12679 + ,110 + ,48020 + ,984885 + ,18014 + ,239 + ,-1710 + ,1405225 + ,43556 + ,247 + ,32648 + ,227132 + ,24811 + ,505 + ,95350 + ,929118 + ,6575 + ,159 + ,151352 + ,1071292 + ,7123 + ,109 + ,288170 + ,638830 + ,21950 + ,519 + ,114337 + ,856956 + ,37597 + ,248 + ,37884 + ,992426 + ,17821 + ,373 + ,122844 + ,444477 + ,12988 + ,119 + ,82340 + ,857217 + ,22330 + ,84 + ,79801 + ,711969 + ,13326 + ,102 + ,165548 + ,702380 + ,16189 + ,295 + ,116384 + ,358589 + ,7146 + ,105 + ,134028 + ,297978 + ,15824 + ,64 + ,63838 + ,585715 + ,27664 + ,282 + ,74996 + ,657954 + ,11920 + ,182 + ,31080 + ,209458 + ,8568 + ,37 + ,32168 + ,786690 + ,14416 + ,361 + ,49857 + ,439798 + ,3369 + ,28 + ,87161 + ,688779 + ,11819 + ,85 + ,106113 + ,574339 + ,6984 + ,45 + ,80570 + ,741409 + ,4519 + ,49 + ,102129 + ,597793 + ,2220 + ,22 + ,301670 + ,644190 + ,18562 + ,155 + ,102313 + ,377934 + ,10327 + ,91 + ,88577 + ,640273 + ,5336 + ,81 + ,112477 + ,697458 + ,2365 + ,79 + ,191778 + ,550608 + ,4069 + ,145 + ,79804 + ,207393 + ,8636 + ,855 + ,128294 + ,301607 + ,13718 + ,61 + ,96448 + ,345783 + ,4525 + ,226 + ,93811 + ,501749 + ,6869 + ,105 + ,117520 + ,379983 + ,4628 + ,62 + ,69159 + ,387475 + ,3689 + ,25 + ,101792 + ,377305 + ,4891 + ,217 + ,210568 + ,370837 + ,7489 + ,322 + ,136996 + ,430866 + ,4901 + ,84 + ,121920 + ,469107 + ,2284 + ,33 + ,76403 + ,194493 + ,3160 + ,108 + ,108094 + ,530670 + ,4150 + ,150 + ,134759 + ,518365 + ,7285 + ,115 + ,188873 + ,491303 + ,1134 + ,162 + ,146216 + ,527021 + ,4658 + ,158 + ,156608 + ,233773 + ,2384 + ,97 + ,61348 + ,405972 + ,3748 + ,9 + ,50350 + ,652925 + ,5371 + ,66 + ,87720 + ,446211 + ,1285 + ,107 + ,99489 + ,341340 + ,9327 + ,101 + ,87419 + ,387699 + ,5565 + ,47 + ,94355 + ,493408 + ,1528 + ,38 + ,60326 + ,146494 + ,3122 + ,34 + ,94670 + ,414462 + ,7561 + ,87 + ,82425 + ,364304 + ,2675 + ,79 + ,59017 + ,355178 + ,13253 + ,947 + ,90829 + ,357760 + ,880 + ,74 + ,80791 + ,261216 + ,2053 + ,53 + ,100423 + ,397144 + ,1424 + ,94 + ,131116 + ,374943 + ,4036 + ,63 + ,100269 + ,424898 + ,3045 + ,58 + ,27330 + ,202055 + ,5119 + ,49 + ,39039 + ,378525 + ,1431 + ,34 + ,106885 + ,310768 + ,554 + ,11 + ,79285 + ,325738 + ,1975 + ,35 + ,118881 + ,394510 + ,1765 + ,20 + ,77623 + ,247060 + ,1012 + ,47 + ,114768 + ,368078 + ,810 + ,43 + ,74015 + ,236761 + ,1280 + ,117 + ,69465 + ,312378 + ,666 + ,171 + ,117869 + ,339836 + ,1380 + ,26 + ,60982 + ,347385 + ,4677 + ,75 + ,90131 + ,426280 + ,876 + ,59 + ,138971 + ,352850 + ,814 + ,18 + ,39625 + ,301881 + ,514 + ,15 + ,102725 + ,377516 + ,5692 + ,72 + ,64239 + ,357312 + ,3642 + ,86 + ,90262 + ,458343 + ,540 + ,14 + ,103960 + ,354228 + ,2099 + ,64 + ,106611 + ,308636 + ,567 + ,11 + ,103345 + ,386212 + ,2001 + ,52 + ,95551 + ,393343 + ,2949 + ,41 + ,82903 + ,378509 + ,2253 + ,99 + ,63593 + ,452469 + ,6533 + ,75 + ,126910 + ,364839 + ,1889 + ,45 + ,37527 + ,358649 + ,3055 + ,43 + ,60247 + ,376641 + ,272 + ,8 + ,112995 + ,429112 + ,1414 + ,198 + ,70184 + ,330546 + ,2564 + ,22 + ,130140 + ,403560 + ,1383 + ,11 + ,73221 + ,317892 + ,1261 + ,33 + ,76114 + ,307528 + ,975 + ,23 + ,90534 + ,235133 + ,3366 + ,80 + ,108479 + ,299243 + ,576 + ,18 + ,113761 + ,314073 + ,1686 + ,40 + ,68696 + ,368186 + ,746 + ,23 + ,71561 + ,269661 + ,3192 + ,60 + ,59831 + ,125390 + ,2045 + ,20 + ,97890 + ,510834 + ,5702 + ,61 + ,101481 + ,321896 + ,1932 + ,36 + ,72954 + ,249898 + ,936 + ,30 + ,67939 + ,408881 + ,3437 + ,47 + ,48022 + ,158492 + ,5131 + ,71 + ,86111 + ,292154 + ,2397 + ,14 + ,74020 + ,289513 + ,1389 + ,9 + ,57530 + ,378049 + ,1503 + ,39 + ,56364 + ,343466 + ,402 + ,26 + ,84990 + ,332743 + ,2239 + ,21 + ,88590 + ,442882 + ,2234 + ,16 + ,77200 + ,214215 + ,837 + ,69 + ,61262 + ,315688 + ,10579 + ,92 + ,110309 + ,375195 + ,875 + ,14 + ,67000 + ,334280 + ,1585 + ,107 + ,93099 + ,355864 + ,1659 + ,29 + ,107577 + ,480382 + ,2647 + ,37 + ,62920 + ,353058 + ,3294 + ,23 + ,75832 + ,217193 + ,0 + ,0 + ,60720 + ,315380 + ,94 + ,7 + ,60793 + ,314533 + ,422 + ,28 + ,57935 + ,318056 + ,0 + ,0 + ,60720 + ,315380 + ,34 + ,8 + ,60630 + ,314353 + ,1558 + ,63 + ,55637 + ,369448 + ,0 + ,0 + ,60720 + ,315380 + ,43 + ,3 + ,60887 + ,312846 + ,645 + ,5 + ,60720 + ,312075 + ,316 + ,9 + ,60505 + ,315009 + ,115 + ,13 + ,60945 + ,318903 + ,5 + ,2 + ,60720 + ,314887 + ,897 + ,5 + ,60720 + ,314913 + ,0 + ,0 + ,60720 + ,315380 + ,389 + ,14 + ,58990 + ,325506 + ,0 + ,0 + ,60720 + ,315380 + ,1002 + ,15 + ,56750 + ,298568 + ,36 + ,3 + ,60894 + ,315834 + ,460 + ,15 + ,63346 + ,329784 + ,309 + ,11 + ,56535 + ,312878 + ,0 + ,0 + ,60720 + ,315380 + ,9 + ,6 + ,60835 + ,314987 + ,271 + ,2 + ,60720 + ,325249 + ,14 + ,1 + ,61016 + ,315877 + ,520 + ,10 + ,58650 + ,291650 + ,1766 + ,73 + ,60438 + ,305959 + ,0 + ,0 + ,60720 + ,315380 + ,458 + ,11 + ,58625 + ,297765 + ,20 + ,3 + ,60938 + ,315245 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,98 + ,2 + ,61490 + ,315236 + ,405 + ,7 + ,60845 + ,336425 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,483 + ,27 + ,60830 + ,306268 + ,454 + ,51 + ,63261 + ,302187 + ,47 + ,3 + ,60720 + ,314882 + ,0 + ,0 + ,60720 + ,315380 + ,757 + ,19 + ,45689 + ,382712 + ,4655 + ,393 + ,60720 + ,341570 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,36 + ,4 + ,61564 + ,312412 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,0 + ,60720 + ,315380 + ,126 + ,10 + ,60951 + ,315547 + ,400 + ,152 + ,60720 + ,313267 + ,71 + ,1 + ,60745 + ,316176 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,972 + ,34 + ,71642 + ,359335 + ,531 + ,10 + ,71641 + ,330068 + ,2461 + ,57 + ,55792 + ,314289 + ,378 + ,52 + ,71873 + ,297413 + ,23 + ,5 + ,62555 + ,314806 + ,638 + ,14 + ,60370 + ,333210 + ,2300 + ,29 + ,64873 + ,352108 + ,149 + ,5 + ,62041 + ,313332 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,4 + ,59500 + ,318745 + ,0 + ,0 + ,60720 + ,315380 + ,141 + ,6 + ,61630 + ,315366 + ,0 + ,0 + ,60720 + ,315380 + ,28 + ,2 + ,60890 + ,315688 + ,0 + ,0 + ,60720 + ,315380 + ,4980 + ,91 + ,113521 + ,409642 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,472 + ,20 + ,80045 + ,269587 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,27 + ,50804 + ,300962 + ,496 + ,17 + ,87390 + ,325479 + ,10 + ,2 + ,61656 + ,316155 + ,63 + ,4 + ,65688 + ,318574 + ,0 + ,0 + ,60720 + ,315380 + ,1136 + ,32 + ,48522 + ,343613 + ,265 + ,31 + ,60720 + ,306948 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,32 + ,57640 + ,330059 + ,474 + ,20 + ,61977 + ,288985 + ,534 + ,7 + ,62620 + ,304485 + ,0 + ,0 + ,60720 + ,315380 + ,15 + ,8 + ,60831 + ,315688 + ,397 + ,28 + ,60646 + ,317736 + ,0 + ,0 + ,60720 + ,315380 + ,1866 + ,29 + ,56225 + ,322331 + ,288 + ,4 + ,60510 + ,296656 + ,0 + ,0 + ,60720 + ,315380 + ,3 + ,2 + ,60698 + ,315354 + ,468 + ,21 + ,60720 + ,312161 + ,20 + ,2 + ,60805 + ,315576 + ,278 + ,26 + ,61404 + ,314922 + ,61 + ,14 + ,60720 + ,314551 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,4 + ,65276 + ,312339 + ,0 + ,0 + ,60720 + ,315380 + ,317 + ,9 + ,63915 + ,298700 + ,738 + ,10 + ,60720 + ,321376 + ,0 + ,0 + ,60720 + ,315380 + ,368 + ,17 + ,61686 + ,303230 + ,0 + ,0 + ,60720 + ,315380 + ,2 + ,1 + ,60743 + ,315487 + ,0 + ,0 + ,60720 + ,315380 + ,53 + ,6 + ,60349 + ,315793 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,94 + ,3 + ,61360 + ,312887 + ,0 + ,0 + ,60720 + ,315380 + ,24 + ,8 + ,59818 + ,315637 + ,2332 + ,4 + ,72680 + ,324385 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,131 + ,11 + ,61808 + ,308989 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,0 + ,60720 + ,315380 + ,167 + ,2 + ,64245 + ,307322 + ,622 + ,73 + ,73007 + ,304376 + ,2328 + ,94 + ,82732 + ,253588 + ,0 + ,0 + ,60720 + ,315380 + ,365 + ,8 + ,54820 + ,309560 + ,364 + ,35 + ,47705 + ,298466 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,12 + ,72835 + ,343929 + ,307 + ,15 + ,58856 + ,331955 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,11 + ,77655 + ,381180 + ,0 + ,0 + ,60720 + ,315380 + ,138 + ,6 + ,69817 + ,331420 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,125 + ,12 + ,60798 + ,310201 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,30 + ,62452 + ,320016 + ,335 + ,33 + ,64175 + ,320398 + ,0 + ,0 + ,60720 + ,315380 + ,1324 + ,117 + ,67440 + ,291841 + ,176 + ,28 + ,68136 + ,310670 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,72 + ,56726 + ,313491 + ,0 + ,0 + ,60720 + ,315380 + ,333 + ,13 + ,70811 + ,331323 + ,0 + ,0 + ,60720 + ,315380 + ,601 + ,6 + ,60720 + ,319210 + ,30 + ,4 + ,62045 + ,318098 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,62 + ,54323 + ,292754 + ,0 + ,0 + ,60720 + ,315380 + ,165 + ,24 + ,62841 + ,325176 + ,453 + ,21 + ,81125 + ,365959 + ,0 + ,0 + ,60720 + ,315380 + ,53 + ,14 + ,59506 + ,302409 + ,382 + ,21 + ,59365 + ,340968 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,30 + ,4 + ,60798 + ,313164 + ,290 + ,2 + ,58790 + ,301164 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,366 + ,53 + ,61808 + ,344425 + ,2 + ,9 + ,60735 + ,315394 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,13 + ,64016 + ,316647 + ,384 + ,22 + ,54683 + ,309836 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,365 + ,83 + ,87192 + ,346611 + ,0 + ,0 + ,60720 + ,315380 + ,49 + ,8 + ,64107 + ,322031 + ,3 + ,4 + ,60761 + ,315656 + ,133 + ,14 + ,65990 + ,339445 + ,32 + ,1 + ,59988 + ,314964 + ,368 + ,17 + ,61167 + ,297141 + ,1 + ,6 + ,60719 + ,315372 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,22 + ,2 + ,60722 + ,312502 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,96 + ,5 + ,60379 + ,313729 + ,1 + ,2 + ,60727 + ,315388 + ,314 + ,5 + ,60720 + ,315371 + ,844 + ,78 + ,60925 + ,296139 + ,0 + ,0 + ,60720 + ,315380 + ,26 + ,1 + ,60896 + ,313880 + ,125 + ,13 + ,59734 + ,317698 + ,304 + ,15 + ,62969 + ,295580 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,621 + ,48 + ,60720 + ,308256 + ,0 + ,0 + ,60720 + ,315380 + ,119 + ,6 + ,59118 + ,303677 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,1595 + ,17 + ,60720 + ,319369 + ,312 + ,14 + ,58598 + ,318690 + ,60 + ,10 + ,61124 + ,314049 + ,587 + ,12 + ,59595 + ,325699 + ,135 + ,2 + ,62065 + ,314210 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,514 + ,52 + ,78780 + ,322378 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,1763 + ,24 + ,61600 + ,308336 + ,180 + ,11 + ,59635 + ,316386 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,21 + ,60720 + ,315553 + ,0 + ,0 + ,60720 + ,315380 + ,448 + ,40 + ,59781 + ,323361 + ,227 + ,9 + ,76644 + ,336639 + ,174 + ,1 + ,64820 + ,307424 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,121 + ,24 + ,56178 + ,295370 + ,607 + ,11 + ,60436 + ,322340 + ,2212 + ,14 + ,60720 + ,319864 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,530 + ,60 + ,73433 + ,317291 + ,571 + ,80 + ,41477 + ,280398 + ,0 + ,0 + ,60720 + ,315380 + ,78 + ,16 + ,62700 + ,317330 + ,2489 + ,40 + ,67804 + ,238125 + ,131 + ,6 + ,59661 + ,327071 + ,923 + ,8 + ,58620 + ,309038 + ,72 + ,3 + ,60398 + ,314210 + ,572 + ,16 + ,58580 + ,307930 + ,397 + ,10 + ,62710 + ,322327 + ,450 + ,8 + ,59325 + ,292136 + ,622 + ,7 + ,60950 + ,263276 + ,694 + ,8 + ,68060 + ,367655 + ,3425 + ,12 + ,83620 + ,283910 + ,562 + ,13 + ,58456 + ,283587 + ,4917 + ,42 + ,52811 + ,243650 + ,1442 + ,118 + ,121173 + ,438493 + ,529 + ,9 + ,63870 + ,296261 + ,2126 + ,138 + ,21001 + ,230621 + ,1061 + ,5 + ,70415 + ,304252 + ,776 + ,9 + ,64230 + ,333505 + ,611 + ,8 + ,59190 + ,296919 + ,1526 + ,25 + ,69351 + ,278990 + ,592 + ,7 + ,64270 + ,276898 + ,1182 + ,13 + ,70694 + ,327007 + ,621 + ,16 + ,68005 + ,317046 + ,989 + ,11 + ,58930 + ,304555 + ,438 + ,11 + ,58320 + ,298096 + ,726 + ,3 + ,69980 + ,231861 + ,1303 + ,61 + ,69863 + ,309422 + ,7419 + ,29 + ,63255 + ,286963 + ,1164 + ,17 + ,57320 + ,269753 + ,3310 + ,33 + ,75230 + ,448243 + ,1920 + ,15 + ,79420 + ,165404 + ,965 + ,3 + ,73490 + ,204325 + ,3256 + ,66 + ,35250 + ,407159 + ,1135 + ,17 + ,62285 + ,290476 + ,1270 + ,26 + ,69206 + ,275311 + ,661 + ,3 + ,65920 + ,246541 + ,1013 + ,2 + ,69770 + ,253468 + ,2844 + ,67 + ,72683 + ,240897 + ,11528 + ,70 + ,-14545 + ,-83265 + ,6526 + ,26 + ,55830 + ,-42143 + ,2264 + ,24 + ,55174 + ,272713 + ,5109 + ,97 + ,67038 + ,215362 + ,3999 + ,30 + ,51252 + ,42754 + ,35624 + ,223 + ,157278 + ,306275 + ,9252 + ,48 + ,79510 + ,253537 + ,15236 + ,90 + ,77440 + ,372631 + ,18073 + ,180 + ,27284 + ,-7170) + ,dim=c(4 + ,431) + ,dimnames=list(c('Costs' + ,'Trades' + ,'Dividends' + ,'Wealth ') + ,1:431)) > y <- array(NA,dim=c(4,431),dimnames=list(c('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 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '4' > #'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 Costs Trades Dividends 1 6282154 162556 1081 213118 2 4321023 29790 309 81767 3 4111912 87550 458 153198 4 223193 84738 588 -26007 5 1491348 54660 302 126942 6 1629616 42634 156 157214 7 1398893 40949 481 129352 8 1926517 45187 353 234817 9 983660 37704 452 60448 10 1443586 16275 109 47818 11 1073089 25830 115 245546 12 984885 12679 110 48020 13 1405225 18014 239 -1710 14 227132 43556 247 32648 15 929118 24811 505 95350 16 1071292 6575 159 151352 17 638830 7123 109 288170 18 856956 21950 519 114337 19 992426 37597 248 37884 20 444477 17821 373 122844 21 857217 12988 119 82340 22 711969 22330 84 79801 23 702380 13326 102 165548 24 358589 16189 295 116384 25 297978 7146 105 134028 26 585715 15824 64 63838 27 657954 27664 282 74996 28 209458 11920 182 31080 29 786690 8568 37 32168 30 439798 14416 361 49857 31 688779 3369 28 87161 32 574339 11819 85 106113 33 741409 6984 45 80570 34 597793 4519 49 102129 35 644190 2220 22 301670 36 377934 18562 155 102313 37 640273 10327 91 88577 38 697458 5336 81 112477 39 550608 2365 79 191778 40 207393 4069 145 79804 41 301607 8636 855 128294 42 345783 13718 61 96448 43 501749 4525 226 93811 44 379983 6869 105 117520 45 387475 4628 62 69159 46 377305 3689 25 101792 47 370837 4891 217 210568 48 430866 7489 322 136996 49 469107 4901 84 121920 50 194493 2284 33 76403 51 530670 3160 108 108094 52 518365 4150 150 134759 53 491303 7285 115 188873 54 527021 1134 162 146216 55 233773 4658 158 156608 56 405972 2384 97 61348 57 652925 3748 9 50350 58 446211 5371 66 87720 59 341340 1285 107 99489 60 387699 9327 101 87419 61 493408 5565 47 94355 62 146494 1528 38 60326 63 414462 3122 34 94670 64 364304 7561 87 82425 65 355178 2675 79 59017 66 357760 13253 947 90829 67 261216 880 74 80791 68 397144 2053 53 100423 69 374943 1424 94 131116 70 424898 4036 63 100269 71 202055 3045 58 27330 72 378525 5119 49 39039 73 310768 1431 34 106885 74 325738 554 11 79285 75 394510 1975 35 118881 76 247060 1765 20 77623 77 368078 1012 47 114768 78 236761 810 43 74015 79 312378 1280 117 69465 80 339836 666 171 117869 81 347385 1380 26 60982 82 426280 4677 75 90131 83 352850 876 59 138971 84 301881 814 18 39625 85 377516 514 15 102725 86 357312 5692 72 64239 87 458343 3642 86 90262 88 354228 540 14 103960 89 308636 2099 64 106611 90 386212 567 11 103345 91 393343 2001 52 95551 92 378509 2949 41 82903 93 452469 2253 99 63593 94 364839 6533 75 126910 95 358649 1889 45 37527 96 376641 3055 43 60247 97 429112 272 8 112995 98 330546 1414 198 70184 99 403560 2564 22 130140 100 317892 1383 11 73221 101 307528 1261 33 76114 102 235133 975 23 90534 103 299243 3366 80 108479 104 314073 576 18 113761 105 368186 1686 40 68696 106 269661 746 23 71561 107 125390 3192 60 59831 108 510834 2045 20 97890 109 321896 5702 61 101481 110 249898 1932 36 72954 111 408881 936 30 67939 112 158492 3437 47 48022 113 292154 5131 71 86111 114 289513 2397 14 74020 115 378049 1389 9 57530 116 343466 1503 39 56364 117 332743 402 26 84990 118 442882 2239 21 88590 119 214215 2234 16 77200 120 315688 837 69 61262 121 375195 10579 92 110309 122 334280 875 14 67000 123 355864 1585 107 93099 124 480382 1659 29 107577 125 353058 2647 37 62920 126 217193 3294 23 75832 127 315380 0 0 60720 128 314533 94 7 60793 129 318056 422 28 57935 130 315380 0 0 60720 131 314353 34 8 60630 132 369448 1558 63 55637 133 315380 0 0 60720 134 312846 43 3 60887 135 312075 645 5 60720 136 315009 316 9 60505 137 318903 115 13 60945 138 314887 5 2 60720 139 314913 897 5 60720 140 315380 0 0 60720 141 325506 389 14 58990 142 315380 0 0 60720 143 298568 1002 15 56750 144 315834 36 3 60894 145 329784 460 15 63346 146 312878 309 11 56535 147 315380 0 0 60720 148 314987 9 6 60835 149 325249 271 2 60720 150 315877 14 1 61016 151 291650 520 10 58650 152 305959 1766 73 60438 153 315380 0 0 60720 154 297765 458 11 58625 155 315245 20 3 60938 156 315380 0 0 60720 157 315380 0 0 60720 158 315236 98 2 61490 159 336425 405 7 60845 160 315380 0 0 60720 161 315380 0 0 60720 162 315380 0 0 60720 163 315380 0 0 60720 164 306268 483 27 60830 165 302187 454 51 63261 166 314882 47 3 60720 167 315380 0 0 60720 168 382712 757 19 45689 169 341570 4655 393 60720 170 315380 0 0 60720 171 315380 0 0 60720 172 312412 36 4 61564 173 315380 0 0 60720 174 309596 203 9 61938 175 315380 0 0 60720 176 315547 126 10 60951 177 313267 400 152 60720 178 316176 71 1 60745 179 315380 0 0 60720 180 315380 0 0 60720 181 359335 972 34 71642 182 330068 531 10 71641 183 314289 2461 57 55792 184 297413 378 52 71873 185 314806 23 5 62555 186 333210 638 14 60370 187 352108 2300 29 64873 188 313332 149 5 62041 189 291787 226 5 65745 190 315380 0 0 60720 191 318745 275 4 59500 192 315380 0 0 60720 193 315366 141 6 61630 194 315380 0 0 60720 195 315688 28 2 60890 196 315380 0 0 60720 197 409642 4980 91 113521 198 315380 0 0 60720 199 315380 0 0 60720 200 269587 472 20 80045 201 315380 0 0 60720 202 315380 0 0 60720 203 315380 0 0 60720 204 300962 203 27 50804 205 325479 496 17 87390 206 316155 10 2 61656 207 318574 63 4 65688 208 315380 0 0 60720 209 343613 1136 32 48522 210 306948 265 31 60720 211 315380 0 0 60720 212 315380 0 0 60720 213 330059 267 32 57640 214 288985 474 20 61977 215 304485 534 7 62620 216 315380 0 0 60720 217 315688 15 8 60831 218 317736 397 28 60646 219 315380 0 0 60720 220 322331 1866 29 56225 221 296656 288 4 60510 222 315380 0 0 60720 223 315354 3 2 60698 224 312161 468 21 60720 225 315576 20 2 60805 226 314922 278 26 61404 227 314551 61 14 60720 228 315380 0 0 60720 229 312339 192 4 65276 230 315380 0 0 60720 231 298700 317 9 63915 232 321376 738 10 60720 233 315380 0 0 60720 234 303230 368 17 61686 235 315380 0 0 60720 236 315487 2 1 60743 237 315380 0 0 60720 238 315793 53 6 60349 239 315380 0 0 60720 240 315380 0 0 60720 241 315380 0 0 60720 242 312887 94 3 61360 243 315380 0 0 60720 244 315637 24 8 59818 245 324385 2332 4 72680 246 315380 0 0 60720 247 315380 0 0 60720 248 308989 131 11 61808 249 315380 0 0 60720 250 315380 0 0 60720 251 296702 206 9 53110 252 315380 0 0 60720 253 307322 167 2 64245 254 304376 622 73 73007 255 253588 2328 94 82732 256 315380 0 0 60720 257 309560 365 8 54820 258 298466 364 35 47705 259 315380 0 0 60720 260 315380 0 0 60720 261 315380 0 0 60720 262 315380 0 0 60720 263 343929 226 12 72835 264 331955 307 15 58856 265 315380 0 0 60720 266 315380 0 0 60720 267 315380 0 0 60720 268 381180 188 11 77655 269 315380 0 0 60720 270 331420 138 6 69817 271 315380 0 0 60720 272 315380 0 0 60720 273 315380 0 0 60720 274 310201 125 12 60798 275 315380 0 0 60720 276 320016 282 30 62452 277 320398 335 33 64175 278 315380 0 0 60720 279 291841 1324 117 67440 280 310670 176 28 68136 281 315380 0 0 60720 282 315380 0 0 60720 283 313491 249 72 56726 284 315380 0 0 60720 285 331323 333 13 70811 286 315380 0 0 60720 287 319210 601 6 60720 288 318098 30 4 62045 289 315380 0 0 60720 290 292754 249 62 54323 291 315380 0 0 60720 292 325176 165 24 62841 293 365959 453 21 81125 294 315380 0 0 60720 295 302409 53 14 59506 296 340968 382 21 59365 297 315380 0 0 60720 298 315380 0 0 60720 299 315380 0 0 60720 300 315380 0 0 60720 301 313164 30 4 60798 302 301164 290 2 58790 303 315380 0 0 60720 304 315380 0 0 60720 305 344425 366 53 61808 306 315394 2 9 60735 307 315380 0 0 60720 308 316647 209 13 64016 309 309836 384 22 54683 310 315380 0 0 60720 311 315380 0 0 60720 312 346611 365 83 87192 313 315380 0 0 60720 314 322031 49 8 64107 315 315656 3 4 60761 316 339445 133 14 65990 317 314964 32 1 59988 318 297141 368 17 61167 319 315372 1 6 60719 320 315380 0 0 60720 321 315380 0 0 60720 322 315380 0 0 60720 323 315380 0 0 60720 324 315380 0 0 60720 325 315380 0 0 60720 326 312502 22 2 60722 327 315380 0 0 60720 328 315380 0 0 60720 329 315380 0 0 60720 330 315380 0 0 60720 331 315380 0 0 60720 332 315380 0 0 60720 333 315380 0 0 60720 334 313729 96 5 60379 335 315388 1 2 60727 336 315371 314 5 60720 337 296139 844 78 60925 338 315380 0 0 60720 339 313880 26 1 60896 340 317698 125 13 59734 341 295580 304 15 62969 342 315380 0 0 60720 343 315380 0 0 60720 344 315380 0 0 60720 345 308256 621 48 60720 346 315380 0 0 60720 347 303677 119 6 59118 348 315380 0 0 60720 349 315380 0 0 60720 350 319369 1595 17 60720 351 318690 312 14 58598 352 314049 60 10 61124 353 325699 587 12 59595 354 314210 135 2 62065 355 315380 0 0 60720 356 315380 0 0 60720 357 322378 514 52 78780 358 315380 0 0 60720 359 315380 0 0 60720 360 315380 0 0 60720 361 315398 1 4 60722 362 315380 0 0 60720 363 315380 0 0 60720 364 308336 1763 24 61600 365 316386 180 11 59635 366 315380 0 0 60720 367 315380 0 0 60720 368 315380 0 0 60720 369 315380 0 0 60720 370 315553 218 21 60720 371 315380 0 0 60720 372 323361 448 40 59781 373 336639 227 9 76644 374 307424 174 1 64820 375 315380 0 0 60720 376 315380 0 0 60720 377 295370 121 24 56178 378 322340 607 11 60436 379 319864 2212 14 60720 380 315380 0 0 60720 381 315380 0 0 60720 382 317291 530 60 73433 383 280398 571 80 41477 384 315380 0 0 60720 385 317330 78 16 62700 386 238125 2489 40 67804 387 327071 131 6 59661 388 309038 923 8 58620 389 314210 72 3 60398 390 307930 572 16 58580 391 322327 397 10 62710 392 292136 450 8 59325 393 263276 622 7 60950 394 367655 694 8 68060 395 283910 3425 12 83620 396 283587 562 13 58456 397 243650 4917 42 52811 398 438493 1442 118 121173 399 296261 529 9 63870 400 230621 2126 138 21001 401 304252 1061 5 70415 402 333505 776 9 64230 403 296919 611 8 59190 404 278990 1526 25 69351 405 276898 592 7 64270 406 327007 1182 13 70694 407 317046 621 16 68005 408 304555 989 11 58930 409 298096 438 11 58320 410 231861 726 3 69980 411 309422 1303 61 69863 412 286963 7419 29 63255 413 269753 1164 17 57320 414 448243 3310 33 75230 415 165404 1920 15 79420 416 204325 965 3 73490 417 407159 3256 66 35250 418 290476 1135 17 62285 419 275311 1270 26 69206 420 246541 661 3 65920 421 253468 1013 2 69770 422 240897 2844 67 72683 423 -83265 11528 70 -14545 424 -42143 6526 26 55830 425 272713 2264 24 55174 426 215362 5109 97 67038 427 42754 3999 30 51252 428 306275 35624 223 157278 429 253537 9252 48 79510 430 372631 15236 90 77440 431 -7170 18073 180 27284 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Costs Trades Dividends 159337.317 30.122 -404.577 2.000 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2198710 -4394 34580 34580 3225801 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.593e+05 3.072e+04 5.186 3.32e-07 *** Costs 3.012e+01 1.576e+00 19.114 < 2e-16 *** Trades -4.046e+02 1.669e+02 -2.424 0.0158 * Dividends 2.000e+00 4.103e-01 4.876 1.53e-06 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 254200 on 427 degrees of freedom Multiple R-squared: 0.6537, Adjusted R-squared: 0.6512 F-statistic: 268.7 on 3 and 427 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 9.875609e-93 4.937805e-93 [2,] 1.0000000 6.214917e-110 3.107459e-110 [3,] 1.0000000 4.427687e-110 2.213844e-110 [4,] 1.0000000 2.773469e-123 1.386734e-123 [5,] 1.0000000 6.380867e-130 3.190434e-130 [6,] 1.0000000 5.015475e-134 2.507738e-134 [7,] 1.0000000 1.772949e-153 8.864743e-154 [8,] 1.0000000 1.107371e-160 5.536855e-161 [9,] 1.0000000 9.758343e-165 4.879171e-165 [10,] 1.0000000 1.633175e-172 8.165877e-173 [11,] 1.0000000 3.984119e-175 1.992060e-175 [12,] 1.0000000 6.034183e-178 3.017092e-178 [13,] 1.0000000 3.954370e-183 1.977185e-183 [14,] 1.0000000 2.586062e-183 1.293031e-183 [15,] 1.0000000 7.113834e-189 3.556917e-189 [16,] 1.0000000 1.139728e-190 5.698639e-191 [17,] 1.0000000 2.389425e-191 1.194713e-191 [18,] 1.0000000 2.283247e-191 1.141623e-191 [19,] 1.0000000 2.697203e-191 1.348602e-191 [20,] 1.0000000 5.183394e-192 2.591697e-192 [21,] 1.0000000 1.198780e-194 5.993899e-195 [22,] 1.0000000 9.930032e-194 4.965016e-194 [23,] 1.0000000 1.866795e-204 9.333976e-205 [24,] 1.0000000 1.217792e-204 6.088962e-205 [25,] 1.0000000 4.154427e-209 2.077214e-209 [26,] 1.0000000 3.191013e-210 1.595506e-210 [27,] 1.0000000 2.590078e-218 1.295039e-218 [28,] 1.0000000 2.928903e-220 1.464452e-220 [29,] 1.0000000 6.016902e-220 3.008451e-220 [30,] 1.0000000 7.060014e-220 3.530007e-220 [31,] 1.0000000 4.325746e-225 2.162873e-225 [32,] 1.0000000 5.370087e-231 2.685043e-231 [33,] 1.0000000 5.293847e-230 2.646923e-230 [34,] 1.0000000 7.734743e-230 3.867372e-230 [35,] 1.0000000 1.359907e-229 6.799535e-230 [36,] 1.0000000 4.296520e-229 2.148260e-229 [37,] 1.0000000 1.600355e-229 8.001777e-230 [38,] 1.0000000 1.782341e-228 8.911703e-229 [39,] 1.0000000 7.587601e-228 3.793800e-228 [40,] 1.0000000 9.591632e-227 4.795816e-227 [41,] 1.0000000 5.476649e-228 2.738324e-228 [42,] 1.0000000 6.657804e-227 3.328902e-227 [43,] 1.0000000 3.104703e-226 1.552351e-226 [44,] 1.0000000 2.190141e-226 1.095071e-226 [45,] 1.0000000 4.576708e-227 2.288354e-227 [46,] 1.0000000 9.139549e-227 4.569775e-227 [47,] 1.0000000 6.981042e-226 3.490521e-226 [48,] 1.0000000 2.494295e-225 1.247148e-225 [49,] 1.0000000 6.357183e-228 3.178591e-228 [50,] 1.0000000 1.013824e-227 5.069121e-228 [51,] 1.0000000 1.108036e-239 5.540181e-240 [52,] 1.0000000 6.622115e-240 3.311057e-240 [53,] 1.0000000 6.961421e-239 3.480710e-239 [54,] 1.0000000 1.452938e-238 7.264689e-239 [55,] 1.0000000 5.350273e-240 2.675137e-240 [56,] 1.0000000 2.628926e-241 1.314463e-241 [57,] 1.0000000 1.262599e-240 6.312994e-241 [58,] 1.0000000 5.672154e-240 2.836077e-240 [59,] 1.0000000 3.047876e-239 1.523938e-239 [60,] 1.0000000 3.783268e-238 1.891634e-238 [61,] 1.0000000 1.046125e-237 5.230625e-238 [62,] 1.0000000 1.109933e-236 5.549664e-237 [63,] 1.0000000 7.987114e-236 3.993557e-236 [64,] 1.0000000 2.963754e-235 1.481877e-235 [65,] 1.0000000 2.974734e-234 1.487367e-234 [66,] 1.0000000 3.048998e-235 1.524499e-235 [67,] 1.0000000 1.135725e-234 5.678627e-235 [68,] 1.0000000 1.387209e-233 6.936044e-234 [69,] 1.0000000 1.802432e-232 9.012162e-233 [70,] 1.0000000 4.573414e-232 2.286707e-232 [71,] 1.0000000 4.359277e-231 2.179639e-231 [72,] 1.0000000 5.873097e-231 2.936548e-231 [73,] 1.0000000 6.755034e-230 3.377517e-230 [74,] 1.0000000 1.169430e-229 5.847150e-230 [75,] 1.0000000 8.288844e-229 4.144422e-229 [76,] 1.0000000 1.217363e-228 6.086817e-229 [77,] 1.0000000 8.132851e-229 4.066425e-229 [78,] 1.0000000 6.296147e-228 3.148074e-228 [79,] 1.0000000 7.385872e-227 3.692936e-227 [80,] 1.0000000 1.583384e-226 7.916921e-227 [81,] 1.0000000 6.919061e-227 3.459531e-227 [82,] 1.0000000 5.805052e-226 2.902526e-226 [83,] 1.0000000 1.156444e-225 5.782221e-226 [84,] 1.0000000 1.378359e-224 6.891795e-225 [85,] 1.0000000 1.430831e-223 7.154154e-224 [86,] 1.0000000 1.003700e-222 5.018501e-223 [87,] 1.0000000 2.856121e-224 1.428060e-224 [88,] 1.0000000 2.159437e-223 1.079718e-223 [89,] 1.0000000 7.728519e-224 3.864259e-224 [90,] 1.0000000 8.297363e-224 4.148681e-224 [91,] 1.0000000 9.551591e-223 4.775795e-223 [92,] 1.0000000 1.062431e-221 5.312154e-222 [93,] 1.0000000 9.105336e-221 4.552668e-221 [94,] 1.0000000 1.071484e-219 5.357421e-220 [95,] 1.0000000 1.042746e-218 5.213732e-219 [96,] 1.0000000 7.487027e-220 3.743514e-220 [97,] 1.0000000 5.823251e-220 2.911625e-220 [98,] 1.0000000 9.018533e-221 4.509266e-221 [99,] 1.0000000 4.792609e-220 2.396305e-220 [100,] 1.0000000 1.773933e-219 8.869667e-220 [101,] 1.0000000 9.041595e-222 4.520797e-222 [102,] 1.0000000 5.532620e-223 2.766310e-223 [103,] 1.0000000 3.713151e-222 1.856576e-222 [104,] 1.0000000 7.065588e-222 3.532794e-222 [105,] 1.0000000 6.374601e-222 3.187300e-222 [106,] 1.0000000 3.840291e-222 1.920146e-222 [107,] 1.0000000 2.653352e-221 1.326676e-221 [108,] 1.0000000 2.305012e-220 1.152506e-220 [109,] 1.0000000 1.950256e-220 9.751279e-221 [110,] 1.0000000 8.449564e-220 4.224782e-220 [111,] 1.0000000 6.552957e-219 3.276479e-219 [112,] 1.0000000 6.635318e-219 3.317659e-219 [113,] 1.0000000 6.463719e-220 3.231859e-220 [114,] 1.0000000 7.504225e-219 3.752112e-219 [115,] 1.0000000 3.615875e-218 1.807938e-218 [116,] 1.0000000 3.915776e-217 1.957888e-217 [117,] 1.0000000 3.478489e-216 1.739245e-216 [118,] 1.0000000 6.400805e-216 3.200402e-216 [119,] 1.0000000 2.201908e-215 1.100954e-215 [120,] 1.0000000 6.604227e-216 3.302114e-216 [121,] 1.0000000 7.549642e-215 3.774821e-215 [122,] 1.0000000 8.636219e-214 4.318109e-214 [123,] 1.0000000 9.041269e-213 4.520635e-213 [124,] 1.0000000 1.022146e-211 5.110732e-212 [125,] 1.0000000 1.154369e-210 5.771843e-211 [126,] 1.0000000 1.257251e-210 6.286256e-211 [127,] 1.0000000 1.416033e-209 7.080166e-210 [128,] 1.0000000 1.604869e-208 8.024343e-209 [129,] 1.0000000 1.808183e-207 9.040913e-208 [130,] 1.0000000 2.003049e-206 1.001525e-206 [131,] 1.0000000 2.180228e-205 1.090114e-205 [132,] 1.0000000 2.409455e-204 1.204728e-204 [133,] 1.0000000 2.575562e-203 1.287781e-203 [134,] 1.0000000 2.815797e-202 1.407898e-202 [135,] 1.0000000 2.599239e-201 1.299620e-201 [136,] 1.0000000 2.818858e-200 1.409429e-200 [137,] 1.0000000 3.105819e-199 1.552909e-199 [138,] 1.0000000 3.336564e-198 1.668282e-198 [139,] 1.0000000 3.284344e-197 1.642172e-197 [140,] 1.0000000 3.274545e-196 1.637273e-196 [141,] 1.0000000 3.474687e-195 1.737344e-195 [142,] 1.0000000 3.680834e-194 1.840417e-194 [143,] 1.0000000 3.505095e-193 1.752548e-193 [144,] 1.0000000 3.674480e-192 1.837240e-192 [145,] 1.0000000 3.829009e-191 1.914504e-191 [146,] 1.0000000 4.041399e-190 2.020699e-190 [147,] 1.0000000 4.175814e-189 2.087907e-189 [148,] 1.0000000 4.380274e-188 2.190137e-188 [149,] 1.0000000 4.497351e-187 2.248676e-187 [150,] 1.0000000 4.585045e-186 2.292523e-186 [151,] 1.0000000 4.653728e-185 2.326864e-185 [152,] 1.0000000 4.735252e-184 2.367626e-184 [153,] 1.0000000 3.569242e-183 1.784621e-183 [154,] 1.0000000 3.577261e-182 1.788630e-182 [155,] 1.0000000 3.569299e-181 1.784650e-181 [156,] 1.0000000 3.545443e-180 1.772722e-180 [157,] 1.0000000 3.505994e-179 1.752997e-179 [158,] 1.0000000 3.513381e-178 1.756691e-178 [159,] 1.0000000 3.277433e-177 1.638717e-177 [160,] 1.0000000 3.203972e-176 1.601986e-176 [161,] 1.0000000 3.112640e-175 1.556320e-175 [162,] 1.0000000 3.607545e-176 1.803772e-176 [163,] 1.0000000 1.341746e-175 6.708729e-176 [164,] 1.0000000 1.316268e-174 6.581342e-175 [165,] 1.0000000 1.285459e-173 6.427295e-174 [166,] 1.0000000 1.264426e-172 6.322130e-173 [167,] 1.0000000 1.223736e-171 6.118678e-172 [168,] 1.0000000 1.196037e-170 5.980186e-171 [169,] 1.0000000 1.147125e-169 5.735624e-170 [170,] 1.0000000 1.092220e-168 5.461100e-169 [171,] 1.0000000 9.880837e-168 4.940418e-168 [172,] 1.0000000 9.309706e-167 4.654853e-167 [173,] 1.0000000 8.774024e-166 4.387012e-166 [174,] 1.0000000 8.231790e-165 4.115895e-165 [175,] 1.0000000 5.620470e-164 2.810235e-164 [176,] 1.0000000 5.322883e-163 2.661441e-163 [177,] 1.0000000 3.053001e-162 1.526500e-162 [178,] 1.0000000 1.923482e-161 9.617409e-162 [179,] 1.0000000 1.791673e-160 8.958367e-161 [180,] 1.0000000 1.240046e-159 6.200232e-160 [181,] 1.0000000 4.710561e-159 2.355280e-159 [182,] 1.0000000 4.353831e-158 2.176916e-158 [183,] 1.0000000 3.141849e-157 1.570925e-157 [184,] 1.0000000 2.844780e-156 1.422390e-156 [185,] 1.0000000 2.422479e-155 1.211240e-155 [186,] 1.0000000 2.173839e-154 1.086920e-154 [187,] 1.0000000 1.951216e-153 9.756079e-154 [188,] 1.0000000 1.735142e-152 8.675712e-153 [189,] 1.0000000 1.535140e-151 7.675702e-152 [190,] 1.0000000 1.352829e-150 6.764143e-151 [191,] 1.0000000 1.085638e-149 5.428188e-150 [192,] 1.0000000 9.484256e-149 4.742128e-149 [193,] 1.0000000 8.248111e-148 4.124056e-148 [194,] 1.0000000 8.572229e-148 4.286114e-148 [195,] 1.0000000 7.463794e-147 3.731897e-147 [196,] 1.0000000 6.469330e-146 3.234665e-146 [197,] 1.0000000 5.582023e-145 2.791012e-145 [198,] 1.0000000 4.336191e-144 2.168095e-144 [199,] 1.0000000 1.913698e-143 9.568490e-144 [200,] 1.0000000 1.643158e-142 8.215791e-143 [201,] 1.0000000 1.410649e-141 7.053244e-142 [202,] 1.0000000 1.194740e-140 5.973702e-141 [203,] 1.0000000 1.708399e-140 8.541995e-141 [204,] 1.0000000 1.461851e-139 7.309254e-140 [205,] 1.0000000 1.230659e-138 6.153293e-139 [206,] 1.0000000 1.031322e-137 5.156612e-138 [207,] 1.0000000 6.314833e-137 3.157417e-137 [208,] 1.0000000 4.820284e-136 2.410142e-136 [209,] 1.0000000 4.036397e-135 2.018198e-135 [210,] 1.0000000 3.328507e-134 1.664253e-134 [211,] 1.0000000 2.725239e-133 1.362620e-133 [212,] 1.0000000 2.139781e-132 1.069890e-132 [213,] 1.0000000 1.740816e-131 8.704079e-132 [214,] 1.0000000 8.436053e-131 4.218026e-131 [215,] 1.0000000 6.768944e-130 3.384472e-130 [216,] 1.0000000 5.444065e-129 2.722032e-129 [217,] 1.0000000 4.355243e-128 2.177622e-128 [218,] 1.0000000 3.455383e-127 1.727691e-127 [219,] 1.0000000 2.738613e-126 1.369306e-126 [220,] 1.0000000 2.152252e-125 1.076126e-125 [221,] 1.0000000 1.688557e-124 8.442784e-125 [222,] 1.0000000 1.321247e-123 6.606237e-124 [223,] 1.0000000 1.032621e-122 5.163105e-123 [224,] 1.0000000 8.007525e-122 4.003763e-122 [225,] 1.0000000 5.887986e-121 2.943993e-121 [226,] 1.0000000 4.140144e-120 2.070072e-120 [227,] 1.0000000 3.168814e-119 1.584407e-119 [228,] 1.0000000 2.440999e-118 1.220499e-118 [229,] 1.0000000 1.851262e-117 9.256310e-118 [230,] 1.0000000 1.396814e-116 6.984068e-117 [231,] 1.0000000 1.049519e-115 5.247597e-116 [232,] 1.0000000 7.777962e-115 3.888981e-115 [233,] 1.0000000 5.789815e-114 2.894907e-114 [234,] 1.0000000 4.289704e-113 2.144852e-113 [235,] 1.0000000 3.163373e-112 1.581687e-112 [236,] 1.0000000 2.345837e-111 1.172918e-111 [237,] 1.0000000 1.713720e-110 8.568598e-111 [238,] 1.0000000 1.227663e-109 6.138316e-110 [239,] 1.0000000 8.887127e-109 4.443563e-109 [240,] 1.0000000 6.400576e-108 3.200288e-108 [241,] 1.0000000 4.587926e-107 2.293963e-107 [242,] 1.0000000 3.321162e-106 1.660581e-106 [243,] 1.0000000 2.358307e-105 1.179153e-105 [244,] 1.0000000 1.666627e-104 8.333135e-105 [245,] 1.0000000 1.155720e-103 5.778598e-104 [246,] 1.0000000 8.089828e-103 4.044914e-103 [247,] 1.0000000 5.631279e-102 2.815640e-102 [248,] 1.0000000 3.210088e-101 1.605044e-101 [249,] 1.0000000 2.432782e-101 1.216391e-101 [250,] 1.0000000 1.692684e-100 8.463420e-101 [251,] 1.0000000 1.066985e-99 5.334926e-100 [252,] 1.0000000 6.157806e-99 3.078903e-99 [253,] 1.0000000 4.227291e-98 2.113646e-98 [254,] 1.0000000 2.887882e-97 1.443941e-97 [255,] 1.0000000 1.963235e-96 9.816177e-97 [256,] 1.0000000 1.328114e-95 6.640571e-96 [257,] 1.0000000 9.029732e-95 4.514866e-95 [258,] 1.0000000 4.800266e-94 2.400133e-94 [259,] 1.0000000 3.202844e-93 1.601422e-93 [260,] 1.0000000 2.126437e-92 1.063219e-92 [261,] 1.0000000 1.404784e-91 7.023921e-92 [262,] 1.0000000 7.637153e-91 3.818576e-91 [263,] 1.0000000 4.998906e-90 2.499453e-90 [264,] 1.0000000 3.317701e-89 1.658851e-89 [265,] 1.0000000 2.150197e-88 1.075099e-88 [266,] 1.0000000 1.386521e-87 6.932604e-88 [267,] 1.0000000 8.895584e-87 4.447792e-87 [268,] 1.0000000 5.769418e-86 2.884709e-86 [269,] 1.0000000 3.664388e-85 1.832194e-85 [270,] 1.0000000 2.309817e-84 1.154908e-84 [271,] 1.0000000 1.469383e-83 7.346913e-84 [272,] 1.0000000 9.193088e-83 4.596544e-83 [273,] 1.0000000 4.780492e-82 2.390246e-82 [274,] 1.0000000 2.846048e-81 1.423024e-81 [275,] 1.0000000 1.755835e-80 8.779173e-81 [276,] 1.0000000 1.077606e-79 5.388029e-80 [277,] 1.0000000 6.430182e-79 3.215091e-79 [278,] 1.0000000 3.907637e-78 1.953819e-78 [279,] 1.0000000 2.413419e-77 1.206709e-77 [280,] 1.0000000 1.451536e-76 7.257680e-77 [281,] 1.0000000 8.238338e-76 4.119169e-76 [282,] 1.0000000 4.922932e-75 2.461466e-75 [283,] 1.0000000 2.914186e-74 1.457093e-74 [284,] 1.0000000 1.749147e-73 8.745736e-74 [285,] 1.0000000 1.024604e-72 5.123021e-73 [286,] 1.0000000 5.881046e-72 2.940523e-72 [287,] 1.0000000 3.444738e-71 1.722369e-71 [288,] 1.0000000 1.985892e-70 9.929459e-71 [289,] 1.0000000 1.164554e-69 5.822770e-70 [290,] 1.0000000 4.858391e-69 2.429195e-69 [291,] 1.0000000 2.760285e-68 1.380143e-68 [292,] 1.0000000 1.559553e-67 7.797767e-68 [293,] 1.0000000 8.762331e-67 4.381165e-67 [294,] 1.0000000 4.895535e-66 2.447767e-66 [295,] 1.0000000 2.747026e-65 1.373513e-65 [296,] 1.0000000 1.549794e-64 7.748968e-65 [297,] 1.0000000 8.513971e-64 4.256986e-64 [298,] 1.0000000 4.650528e-63 2.325264e-63 [299,] 1.0000000 2.050656e-62 1.025328e-62 [300,] 1.0000000 1.113027e-61 5.565137e-62 [301,] 1.0000000 5.987148e-61 2.993574e-61 [302,] 1.0000000 3.272556e-60 1.636278e-60 [303,] 1.0000000 1.605672e-59 8.028362e-60 [304,] 1.0000000 8.491598e-59 4.245799e-59 [305,] 1.0000000 4.464165e-58 2.232082e-58 [306,] 1.0000000 1.905455e-57 9.527277e-58 [307,] 1.0000000 9.915633e-57 4.957817e-57 [308,] 1.0000000 5.195221e-56 2.597611e-56 [309,] 1.0000000 2.675413e-55 1.337706e-55 [310,] 1.0000000 1.298822e-54 6.494109e-55 [311,] 1.0000000 6.543833e-54 3.271916e-54 [312,] 1.0000000 3.391570e-53 1.695785e-53 [313,] 1.0000000 1.708983e-52 8.544914e-53 [314,] 1.0000000 8.527669e-52 4.263835e-52 [315,] 1.0000000 4.228446e-51 2.114223e-51 [316,] 1.0000000 2.083387e-50 1.041693e-50 [317,] 1.0000000 1.019948e-49 5.099738e-50 [318,] 1.0000000 4.961176e-49 2.480588e-49 [319,] 1.0000000 2.397560e-48 1.198780e-48 [320,] 1.0000000 1.165094e-47 5.825471e-48 [321,] 1.0000000 5.557371e-47 2.778685e-47 [322,] 1.0000000 2.633245e-46 1.316622e-46 [323,] 1.0000000 1.239377e-45 6.196887e-46 [324,] 1.0000000 5.794059e-45 2.897030e-45 [325,] 1.0000000 2.690326e-44 1.345163e-44 [326,] 1.0000000 1.240637e-43 6.203183e-44 [327,] 1.0000000 5.681673e-43 2.840836e-43 [328,] 1.0000000 2.592458e-42 1.296229e-42 [329,] 1.0000000 1.172813e-41 5.864067e-42 [330,] 1.0000000 5.207948e-41 2.603974e-41 [331,] 1.0000000 2.372217e-40 1.186109e-40 [332,] 1.0000000 1.049442e-39 5.247209e-40 [333,] 1.0000000 4.648204e-39 2.324102e-39 [334,] 1.0000000 1.991062e-38 9.955308e-39 [335,] 1.0000000 8.885298e-38 4.442649e-38 [336,] 1.0000000 3.822618e-37 1.911309e-37 [337,] 1.0000000 1.632186e-36 8.160930e-37 [338,] 1.0000000 6.916143e-36 3.458071e-36 [339,] 1.0000000 3.017103e-35 1.508552e-35 [340,] 1.0000000 1.259771e-34 6.298855e-35 [341,] 1.0000000 5.363560e-34 2.681780e-34 [342,] 1.0000000 2.205094e-33 1.102547e-33 [343,] 1.0000000 8.993246e-33 4.496623e-33 [344,] 1.0000000 3.282821e-32 1.641410e-32 [345,] 1.0000000 1.252869e-31 6.264343e-32 [346,] 1.0000000 5.081274e-31 2.540637e-31 [347,] 1.0000000 1.803898e-30 9.019490e-31 [348,] 1.0000000 7.196742e-30 3.598371e-30 [349,] 1.0000000 2.796494e-29 1.398247e-29 [350,] 1.0000000 1.077198e-28 5.385992e-29 [351,] 1.0000000 3.915609e-28 1.957804e-28 [352,] 1.0000000 1.486898e-27 7.434490e-28 [353,] 1.0000000 5.595423e-27 2.797712e-27 [354,] 1.0000000 2.086393e-26 1.043196e-26 [355,] 1.0000000 7.752073e-26 3.876036e-26 [356,] 1.0000000 2.837219e-25 1.418610e-25 [357,] 1.0000000 1.028468e-24 5.142340e-25 [358,] 1.0000000 3.672638e-24 1.836319e-24 [359,] 1.0000000 1.291124e-23 6.455618e-24 [360,] 1.0000000 4.548003e-23 2.274001e-23 [361,] 1.0000000 1.585700e-22 7.928500e-23 [362,] 1.0000000 5.471218e-22 2.735609e-22 [363,] 1.0000000 1.867755e-21 9.338776e-22 [364,] 1.0000000 6.464886e-21 3.232443e-21 [365,] 1.0000000 2.161510e-20 1.080755e-20 [366,] 1.0000000 7.127188e-20 3.563594e-20 [367,] 1.0000000 2.517037e-19 1.258518e-19 [368,] 1.0000000 8.712157e-19 4.356079e-19 [369,] 1.0000000 2.797258e-18 1.398629e-18 [370,] 1.0000000 8.873494e-18 4.436747e-18 [371,] 1.0000000 2.951247e-17 1.475623e-17 [372,] 1.0000000 8.645670e-17 4.322835e-17 [373,] 1.0000000 2.275641e-16 1.137821e-16 [374,] 1.0000000 6.877425e-16 3.438713e-16 [375,] 1.0000000 2.050173e-15 1.025086e-15 [376,] 1.0000000 6.521844e-15 3.260922e-15 [377,] 1.0000000 1.983022e-14 9.915108e-15 [378,] 1.0000000 5.728861e-14 2.864430e-14 [379,] 1.0000000 1.707799e-13 8.538995e-14 [380,] 1.0000000 4.468576e-13 2.234288e-13 [381,] 1.0000000 1.156656e-12 5.783280e-13 [382,] 1.0000000 3.095127e-12 1.547563e-12 [383,] 1.0000000 8.453329e-12 4.226665e-12 [384,] 1.0000000 2.298339e-11 1.149169e-11 [385,] 1.0000000 6.040574e-11 3.020287e-11 [386,] 1.0000000 1.696765e-10 8.483827e-11 [387,] 1.0000000 4.916376e-10 2.458188e-10 [388,] 1.0000000 9.655498e-10 4.827749e-10 [389,] 1.0000000 2.666488e-09 1.333244e-09 [390,] 1.0000000 7.261224e-09 3.630612e-09 [391,] 1.0000000 1.882463e-08 9.412317e-09 [392,] 1.0000000 4.058986e-08 2.029493e-08 [393,] 0.9999999 1.067754e-07 5.338771e-08 [394,] 0.9999999 2.798909e-07 1.399455e-07 [395,] 0.9999996 7.160232e-07 3.580116e-07 [396,] 0.9999992 1.529035e-06 7.645177e-07 [397,] 0.9999982 3.539000e-06 1.769500e-06 [398,] 0.9999956 8.791658e-06 4.395829e-06 [399,] 0.9999894 2.129955e-05 1.064977e-05 [400,] 0.9999764 4.714256e-05 2.357128e-05 [401,] 0.9999475 1.050786e-04 5.253932e-05 [402,] 0.9998952 2.095866e-04 1.047933e-04 [403,] 0.9997876 4.248413e-04 2.124206e-04 [404,] 0.9995456 9.087123e-04 4.543562e-04 [405,] 0.9990200 1.960069e-03 9.800345e-04 [406,] 0.9984653 3.069301e-03 1.534651e-03 [407,] 0.9969555 6.088941e-03 3.044471e-03 [408,] 0.9973612 5.277579e-03 2.638790e-03 [409,] 0.9958550 8.289994e-03 4.144997e-03 [410,] 0.9920121 1.597584e-02 7.987918e-03 [411,] 0.9976614 4.677224e-03 2.338612e-03 [412,] 0.9951718 9.656323e-03 4.828162e-03 [413,] 0.9889162 2.216767e-02 1.108384e-02 [414,] 0.9750855 4.982890e-02 2.491445e-02 [415,] 0.9470880 1.058241e-01 5.291205e-02 [416,] 0.8920542 2.158916e-01 1.079458e-01 [417,] 0.8048247 3.903505e-01 1.951753e-01 [418,] 0.8344600 3.310800e-01 1.655400e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1ec8y1291211846.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/263811291211846.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/363811291211846.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/463811291211846.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/5hvpl1291211846.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 1237333.7716 3225800.6189 1194235.0281 -2198709.5189 -446207.8858 6 7 8 9 10 -65315.8652 -58061.4111 79150.5670 -249447.6804 742457.9366 11 12 13 14 15 -308956.7059 392076.2392 803384.0750 -1209577.7725 35999.2506 16 17 18 19 20 475470.1181 -267414.2259 17699.4160 -274855.9885 -346491.1382 21 22 23 24 25 190089.0283 -245640.1815 -148253.8127 -401854.6120 -302236.3074 26 27 28 29 30 -152079.9239 -370608.0313 -297472.3304 319888.6245 -107458.5634 31 32 33 34 35 264934.4893 -118886.7400 228735.6187 117862.7404 -176569.1359 36 37 38 39 40 -482482.7536 29495.2712 185164.6824 -31633.2784 -175484.6937 41 42 43 44 45 -28586.2323 -395020.7371 109887.2197 -178865.3985 -24526.8599 46 47 48 49 50 -86659.7032 -269247.8351 -97824.0452 -47758.9998 -173126.2485 51 52 53 54 55 103613.4092 25139.9739 -218762.7914 106580.6851 -315223.4826 56 57 58 59 60 91349.0492 283612.8987 -23682.0837 -12429.1879 -186594.5055 61 62 63 64 65 -3288.1749 -164170.1642 -14535.6869 -152468.2407 29170.0161 66 67 68 69 70 658.3290 -56301.8915 -3474.4077 -51538.3527 -31098.5552 71 72 73 74 75 -80208.5263 6724.8692 -91727.8629 -4435.9069 -47964.1058 76 77 78 79 80 -112625.8705 -32306.1501 -77635.7304 22864.3015 -6161.5256 81 82 83 84 85 35011.7284 -23890.0975 -86997.5856 46042.0802 3276.6093 86 87 88 89 90 -72852.1034 43537.4967 -23669.5973 -101295.9632 7517.6045 91 92 93 94 95 3632.1632 -18907.1836 138110.3578 -214809.1834 85549.2598 96 97 98 99 100 22161.4187 38786.3001 68328.4234 -84437.6455 -25122.8216 101 102 103 104 105 -28698.3104 -125369.5171 -146117.0270 -82896.3085 36828.6130 106 107 108 109 110 -45990.5534 -225506.3353 102172.4638 -187517.4506 -99005.2886 111 112 113 114 115 97583.6413 -181421.3633 -165268.3173 -84430.1190 65432.0980 116 117 118 119 120 41884.9294 1804.1853 47384.8091 -160370.2516 36507.6963 121 122 123 124 125 -286240.9429 20225.1998 5840.5691 67611.1743 3093.7535 126 127 128 129 130 -183753.0522 34580.2165 33587.7584 41443.9094 34580.2165 131 132 133 134 135 35945.7172 77374.3357 34580.2165 31630.6387 13869.3958 136 137 138 139 140 28761.9295 39448.6018 34745.7604 9116.6456 34580.2165 141 142 143 144 145 42113.4674 34580.2165 1596.0722 34815.4903 35943.7690 146 147 148 149 150 35592.4067 34580.2165 36113.5379 37095.3019 34467.9757 151 152 153 154 155 3373.2999 2061.9484 34580.2165 11810.4517 34620.4264 156 157 158 159 160 34580.2165 34580.2165 30753.1273 46007.7895 34580.2165 161 162 163 164 165 34580.2165 34580.2165 34580.2165 21622.8177 23262.3055 166 167 168 169 170 33880.2124 34580.2165 116864.3690 79550.9637 34580.2165 171 172 173 174 175 34580.2165 30457.8195 34580.2165 23886.1881 34580.2165 176 177 178 179 180 34535.5262 81914.1147 33592.1205 34580.2165 34580.2165 181 182 183 184 185 41164.1864 15473.1509 -7722.3720 3954.9700 31665.6161 186 187 188 189 190 39556.5727 5452.7572 27424.4312 -3849.3352 34580.2165 191 192 193 194 195 33720.4193 34580.2165 30926.1365 34580.2165 34513.8910 196 197 198 199 200 34580.2165 -89970.4929 34580.2165 34580.2165 -55995.9889 201 202 203 204 205 34580.2165 34580.2165 34580.2165 44806.6942 -16733.3662 206 207 208 209 210 33990.8040 27556.9989 34580.2165 65941.5748 30707.7678 211 212 213 214 215 34580.2165 34580.2165 60324.2404 -515.5477 6631.3916 216 217 218 219 220 34580.2165 37450.9613 36453.9569 34580.2165 6047.9156 221 222 223 224 225 9219.4592 34580.2165 35317.0126 25760.2266 34812.8986 226 227 228 229 230 34899.0433 37577.8514 34580.2165 18260.4102 34580.2165 231 232 233 234 235 5601.5458 22391.9327 34580.2165 16290.7635 34580.2165 236 237 238 239 240 34985.5410 34580.2165 36566.3486 34580.2165 34580.2165 241 242 243 244 245 34580.2165 29189.2405 34580.2165 39155.2379 -48965.4619 246 247 248 249 250 34580.2165 34580.2165 26517.1760 34580.2165 34580.2165 251 252 253 254 255 28561.0884 34580.2165 15249.6882 9795.8940 -103337.7607 256 257 258 259 260 34580.2165 32804.4768 46896.8111 34580.2165 34580.2165 261 262 263 264 265 34580.2165 34580.2165 36942.0807 51705.1001 34580.2165 266 267 268 269 270 34580.2165 34580.2165 65291.3572 34580.2165 30693.4734 271 272 273 274 275 34580.2165 34580.2165 34580.2165 30334.8589 34580.2165 276 277 278 279 280 39394.4755 35947.1018 34580.2165 5052.6817 21062.1530 281 282 283 284 285 34580.2165 34580.2165 62309.8561 34580.2165 25566.3500 286 287 288 289 290 34580.2165 22734.3419 35362.3736 34580.2165 42333.9750 291 292 293 294 295 34580.2165 44873.1463 39192.5070 34580.2165 28105.2768 296 297 298 299 300 59868.2221 34580.2165 34580.2165 34580.2165 34580.2165 301 302 303 304 305 32922.8350 16298.6976 34580.2165 34580.2165 71866.7354 306 307 308 309 310 38145.1602 34580.2165 28217.9952 38446.2874 34580.2165 311 312 313 314 315 34580.2165 35442.7760 34580.2165 36216.6003 36302.1434 316 317 318 319 320 49761.1157 35069.1596 11239.9555 36971.5570 34580.2165 321 322 323 324 325 34580.2165 34580.2165 34580.2165 34580.2165 34580.2165 326 327 328 329 330 31844.6853 34580.2165 34580.2165 34580.2165 34580.2165 331 332 333 334 335 34580.2165 34580.2165 34580.2165 32742.5135 35353.2460 336 337 338 339 340 27135.7859 21063.1600 34580.2165 32349.5558 40364.8296 341 342 343 344 345 7192.9443 34580.2165 34580.2165 34580.2165 28170.1365 346 347 348 349 350 34580.2165 24924.7505 34580.2165 34580.2165 -2597.6032 351 352 353 354 355 38401.0084 34679.5158 34322.9287 27462.3996 34580.2165 356 357 358 359 360 34580.2165 11006.8190 34580.2165 34580.2165 34580.2165 361 362 363 364 365 36182.4019 34580.2165 34580.2165 -17619.3897 36785.0008 366 367 368 369 370 34580.2165 34580.2165 34580.2165 34580.2165 36682.7327 371 372 373 374 375 34580.2165 47127.9780 20788.8182 13586.0443 34580.2165 376 377 378 379 380 34580.2165 29720.9807 28274.6000 -21901.6234 34580.2165 381 382 383 384 385 34580.2165 19370.4612 53257.8226 34580.2165 36693.1984 386 387 388 389 390 -115636.0426 46871.0853 7872.9811 33099.2810 20654.4427 391 392 393 394 395 29633.8067 3808.4378 -33887.7287 54504.4319 -141011.2663 396 397 398 399 400 -4353.0222 -152446.6167 41068.9784 -3133.3067 21066.1182 401 402 403 404 405 -25877.9536 25950.4200 4011.8418 -54926.7608 -26003.2963 406 407 408 409 410 -4089.2056 9440.9762 1995.5438 13354.0051 -88117.0685 411 412 413 414 415 -4236.9658 -210650.2878 -32429.7526 52064.9884 -204568.3343 416 417 418 419 420 -129873.5310 105933.4123 -20765.0510 -50199.8918 -63357.6347 421 422 423 424 425 -75139.5888 -122393.5869 -532433.2423 -499218.3037 -55479.1464 426 427 428 429 430 -192725.5730 -327426.9457 -1150522.8360 -324119.0097 -364142.2036 431 -692656.8969 > postscript(file="/var/www/html/freestat/rcomp/tmp/6hvpl1291211846.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 1237333.7716 NA 1 3225800.6189 1237333.7716 2 1194235.0281 3225800.6189 3 -2198709.5189 1194235.0281 4 -446207.8858 -2198709.5189 5 -65315.8652 -446207.8858 6 -58061.4111 -65315.8652 7 79150.5670 -58061.4111 8 -249447.6804 79150.5670 9 742457.9366 -249447.6804 10 -308956.7059 742457.9366 11 392076.2392 -308956.7059 12 803384.0750 392076.2392 13 -1209577.7725 803384.0750 14 35999.2506 -1209577.7725 15 475470.1181 35999.2506 16 -267414.2259 475470.1181 17 17699.4160 -267414.2259 18 -274855.9885 17699.4160 19 -346491.1382 -274855.9885 20 190089.0283 -346491.1382 21 -245640.1815 190089.0283 22 -148253.8127 -245640.1815 23 -401854.6120 -148253.8127 24 -302236.3074 -401854.6120 25 -152079.9239 -302236.3074 26 -370608.0313 -152079.9239 27 -297472.3304 -370608.0313 28 319888.6245 -297472.3304 29 -107458.5634 319888.6245 30 264934.4893 -107458.5634 31 -118886.7400 264934.4893 32 228735.6187 -118886.7400 33 117862.7404 228735.6187 34 -176569.1359 117862.7404 35 -482482.7536 -176569.1359 36 29495.2712 -482482.7536 37 185164.6824 29495.2712 38 -31633.2784 185164.6824 39 -175484.6937 -31633.2784 40 -28586.2323 -175484.6937 41 -395020.7371 -28586.2323 42 109887.2197 -395020.7371 43 -178865.3985 109887.2197 44 -24526.8599 -178865.3985 45 -86659.7032 -24526.8599 46 -269247.8351 -86659.7032 47 -97824.0452 -269247.8351 48 -47758.9998 -97824.0452 49 -173126.2485 -47758.9998 50 103613.4092 -173126.2485 51 25139.9739 103613.4092 52 -218762.7914 25139.9739 53 106580.6851 -218762.7914 54 -315223.4826 106580.6851 55 91349.0492 -315223.4826 56 283612.8987 91349.0492 57 -23682.0837 283612.8987 58 -12429.1879 -23682.0837 59 -186594.5055 -12429.1879 60 -3288.1749 -186594.5055 61 -164170.1642 -3288.1749 62 -14535.6869 -164170.1642 63 -152468.2407 -14535.6869 64 29170.0161 -152468.2407 65 658.3290 29170.0161 66 -56301.8915 658.3290 67 -3474.4077 -56301.8915 68 -51538.3527 -3474.4077 69 -31098.5552 -51538.3527 70 -80208.5263 -31098.5552 71 6724.8692 -80208.5263 72 -91727.8629 6724.8692 73 -4435.9069 -91727.8629 74 -47964.1058 -4435.9069 75 -112625.8705 -47964.1058 76 -32306.1501 -112625.8705 77 -77635.7304 -32306.1501 78 22864.3015 -77635.7304 79 -6161.5256 22864.3015 80 35011.7284 -6161.5256 81 -23890.0975 35011.7284 82 -86997.5856 -23890.0975 83 46042.0802 -86997.5856 84 3276.6093 46042.0802 85 -72852.1034 3276.6093 86 43537.4967 -72852.1034 87 -23669.5973 43537.4967 88 -101295.9632 -23669.5973 89 7517.6045 -101295.9632 90 3632.1632 7517.6045 91 -18907.1836 3632.1632 92 138110.3578 -18907.1836 93 -214809.1834 138110.3578 94 85549.2598 -214809.1834 95 22161.4187 85549.2598 96 38786.3001 22161.4187 97 68328.4234 38786.3001 98 -84437.6455 68328.4234 99 -25122.8216 -84437.6455 100 -28698.3104 -25122.8216 101 -125369.5171 -28698.3104 102 -146117.0270 -125369.5171 103 -82896.3085 -146117.0270 104 36828.6130 -82896.3085 105 -45990.5534 36828.6130 106 -225506.3353 -45990.5534 107 102172.4638 -225506.3353 108 -187517.4506 102172.4638 109 -99005.2886 -187517.4506 110 97583.6413 -99005.2886 111 -181421.3633 97583.6413 112 -165268.3173 -181421.3633 113 -84430.1190 -165268.3173 114 65432.0980 -84430.1190 115 41884.9294 65432.0980 116 1804.1853 41884.9294 117 47384.8091 1804.1853 118 -160370.2516 47384.8091 119 36507.6963 -160370.2516 120 -286240.9429 36507.6963 121 20225.1998 -286240.9429 122 5840.5691 20225.1998 123 67611.1743 5840.5691 124 3093.7535 67611.1743 125 -183753.0522 3093.7535 126 34580.2165 -183753.0522 127 33587.7584 34580.2165 128 41443.9094 33587.7584 129 34580.2165 41443.9094 130 35945.7172 34580.2165 131 77374.3357 35945.7172 132 34580.2165 77374.3357 133 31630.6387 34580.2165 134 13869.3958 31630.6387 135 28761.9295 13869.3958 136 39448.6018 28761.9295 137 34745.7604 39448.6018 138 9116.6456 34745.7604 139 34580.2165 9116.6456 140 42113.4674 34580.2165 141 34580.2165 42113.4674 142 1596.0722 34580.2165 143 34815.4903 1596.0722 144 35943.7690 34815.4903 145 35592.4067 35943.7690 146 34580.2165 35592.4067 147 36113.5379 34580.2165 148 37095.3019 36113.5379 149 34467.9757 37095.3019 150 3373.2999 34467.9757 151 2061.9484 3373.2999 152 34580.2165 2061.9484 153 11810.4517 34580.2165 154 34620.4264 11810.4517 155 34580.2165 34620.4264 156 34580.2165 34580.2165 157 30753.1273 34580.2165 158 46007.7895 30753.1273 159 34580.2165 46007.7895 160 34580.2165 34580.2165 161 34580.2165 34580.2165 162 34580.2165 34580.2165 163 21622.8177 34580.2165 164 23262.3055 21622.8177 165 33880.2124 23262.3055 166 34580.2165 33880.2124 167 116864.3690 34580.2165 168 79550.9637 116864.3690 169 34580.2165 79550.9637 170 34580.2165 34580.2165 171 30457.8195 34580.2165 172 34580.2165 30457.8195 173 23886.1881 34580.2165 174 34580.2165 23886.1881 175 34535.5262 34580.2165 176 81914.1147 34535.5262 177 33592.1205 81914.1147 178 34580.2165 33592.1205 179 34580.2165 34580.2165 180 41164.1864 34580.2165 181 15473.1509 41164.1864 182 -7722.3720 15473.1509 183 3954.9700 -7722.3720 184 31665.6161 3954.9700 185 39556.5727 31665.6161 186 5452.7572 39556.5727 187 27424.4312 5452.7572 188 -3849.3352 27424.4312 189 34580.2165 -3849.3352 190 33720.4193 34580.2165 191 34580.2165 33720.4193 192 30926.1365 34580.2165 193 34580.2165 30926.1365 194 34513.8910 34580.2165 195 34580.2165 34513.8910 196 -89970.4929 34580.2165 197 34580.2165 -89970.4929 198 34580.2165 34580.2165 199 -55995.9889 34580.2165 200 34580.2165 -55995.9889 201 34580.2165 34580.2165 202 34580.2165 34580.2165 203 44806.6942 34580.2165 204 -16733.3662 44806.6942 205 33990.8040 -16733.3662 206 27556.9989 33990.8040 207 34580.2165 27556.9989 208 65941.5748 34580.2165 209 30707.7678 65941.5748 210 34580.2165 30707.7678 211 34580.2165 34580.2165 212 60324.2404 34580.2165 213 -515.5477 60324.2404 214 6631.3916 -515.5477 215 34580.2165 6631.3916 216 37450.9613 34580.2165 217 36453.9569 37450.9613 218 34580.2165 36453.9569 219 6047.9156 34580.2165 220 9219.4592 6047.9156 221 34580.2165 9219.4592 222 35317.0126 34580.2165 223 25760.2266 35317.0126 224 34812.8986 25760.2266 225 34899.0433 34812.8986 226 37577.8514 34899.0433 227 34580.2165 37577.8514 228 18260.4102 34580.2165 229 34580.2165 18260.4102 230 5601.5458 34580.2165 231 22391.9327 5601.5458 232 34580.2165 22391.9327 233 16290.7635 34580.2165 234 34580.2165 16290.7635 235 34985.5410 34580.2165 236 34580.2165 34985.5410 237 36566.3486 34580.2165 238 34580.2165 36566.3486 239 34580.2165 34580.2165 240 34580.2165 34580.2165 241 29189.2405 34580.2165 242 34580.2165 29189.2405 243 39155.2379 34580.2165 244 -48965.4619 39155.2379 245 34580.2165 -48965.4619 246 34580.2165 34580.2165 247 26517.1760 34580.2165 248 34580.2165 26517.1760 249 34580.2165 34580.2165 250 28561.0884 34580.2165 251 34580.2165 28561.0884 252 15249.6882 34580.2165 253 9795.8940 15249.6882 254 -103337.7607 9795.8940 255 34580.2165 -103337.7607 256 32804.4768 34580.2165 257 46896.8111 32804.4768 258 34580.2165 46896.8111 259 34580.2165 34580.2165 260 34580.2165 34580.2165 261 34580.2165 34580.2165 262 36942.0807 34580.2165 263 51705.1001 36942.0807 264 34580.2165 51705.1001 265 34580.2165 34580.2165 266 34580.2165 34580.2165 267 65291.3572 34580.2165 268 34580.2165 65291.3572 269 30693.4734 34580.2165 270 34580.2165 30693.4734 271 34580.2165 34580.2165 272 34580.2165 34580.2165 273 30334.8589 34580.2165 274 34580.2165 30334.8589 275 39394.4755 34580.2165 276 35947.1018 39394.4755 277 34580.2165 35947.1018 278 5052.6817 34580.2165 279 21062.1530 5052.6817 280 34580.2165 21062.1530 281 34580.2165 34580.2165 282 62309.8561 34580.2165 283 34580.2165 62309.8561 284 25566.3500 34580.2165 285 34580.2165 25566.3500 286 22734.3419 34580.2165 287 35362.3736 22734.3419 288 34580.2165 35362.3736 289 42333.9750 34580.2165 290 34580.2165 42333.9750 291 44873.1463 34580.2165 292 39192.5070 44873.1463 293 34580.2165 39192.5070 294 28105.2768 34580.2165 295 59868.2221 28105.2768 296 34580.2165 59868.2221 297 34580.2165 34580.2165 298 34580.2165 34580.2165 299 34580.2165 34580.2165 300 32922.8350 34580.2165 301 16298.6976 32922.8350 302 34580.2165 16298.6976 303 34580.2165 34580.2165 304 71866.7354 34580.2165 305 38145.1602 71866.7354 306 34580.2165 38145.1602 307 28217.9952 34580.2165 308 38446.2874 28217.9952 309 34580.2165 38446.2874 310 34580.2165 34580.2165 311 35442.7760 34580.2165 312 34580.2165 35442.7760 313 36216.6003 34580.2165 314 36302.1434 36216.6003 315 49761.1157 36302.1434 316 35069.1596 49761.1157 317 11239.9555 35069.1596 318 36971.5570 11239.9555 319 34580.2165 36971.5570 320 34580.2165 34580.2165 321 34580.2165 34580.2165 322 34580.2165 34580.2165 323 34580.2165 34580.2165 324 34580.2165 34580.2165 325 31844.6853 34580.2165 326 34580.2165 31844.6853 327 34580.2165 34580.2165 328 34580.2165 34580.2165 329 34580.2165 34580.2165 330 34580.2165 34580.2165 331 34580.2165 34580.2165 332 34580.2165 34580.2165 333 32742.5135 34580.2165 334 35353.2460 32742.5135 335 27135.7859 35353.2460 336 21063.1600 27135.7859 337 34580.2165 21063.1600 338 32349.5558 34580.2165 339 40364.8296 32349.5558 340 7192.9443 40364.8296 341 34580.2165 7192.9443 342 34580.2165 34580.2165 343 34580.2165 34580.2165 344 28170.1365 34580.2165 345 34580.2165 28170.1365 346 24924.7505 34580.2165 347 34580.2165 24924.7505 348 34580.2165 34580.2165 349 -2597.6032 34580.2165 350 38401.0084 -2597.6032 351 34679.5158 38401.0084 352 34322.9287 34679.5158 353 27462.3996 34322.9287 354 34580.2165 27462.3996 355 34580.2165 34580.2165 356 11006.8190 34580.2165 357 34580.2165 11006.8190 358 34580.2165 34580.2165 359 34580.2165 34580.2165 360 36182.4019 34580.2165 361 34580.2165 36182.4019 362 34580.2165 34580.2165 363 -17619.3897 34580.2165 364 36785.0008 -17619.3897 365 34580.2165 36785.0008 366 34580.2165 34580.2165 367 34580.2165 34580.2165 368 34580.2165 34580.2165 369 36682.7327 34580.2165 370 34580.2165 36682.7327 371 47127.9780 34580.2165 372 20788.8182 47127.9780 373 13586.0443 20788.8182 374 34580.2165 13586.0443 375 34580.2165 34580.2165 376 29720.9807 34580.2165 377 28274.6000 29720.9807 378 -21901.6234 28274.6000 379 34580.2165 -21901.6234 380 34580.2165 34580.2165 381 19370.4612 34580.2165 382 53257.8226 19370.4612 383 34580.2165 53257.8226 384 36693.1984 34580.2165 385 -115636.0426 36693.1984 386 46871.0853 -115636.0426 387 7872.9811 46871.0853 388 33099.2810 7872.9811 389 20654.4427 33099.2810 390 29633.8067 20654.4427 391 3808.4378 29633.8067 392 -33887.7287 3808.4378 393 54504.4319 -33887.7287 394 -141011.2663 54504.4319 395 -4353.0222 -141011.2663 396 -152446.6167 -4353.0222 397 41068.9784 -152446.6167 398 -3133.3067 41068.9784 399 21066.1182 -3133.3067 400 -25877.9536 21066.1182 401 25950.4200 -25877.9536 402 4011.8418 25950.4200 403 -54926.7608 4011.8418 404 -26003.2963 -54926.7608 405 -4089.2056 -26003.2963 406 9440.9762 -4089.2056 407 1995.5438 9440.9762 408 13354.0051 1995.5438 409 -88117.0685 13354.0051 410 -4236.9658 -88117.0685 411 -210650.2878 -4236.9658 412 -32429.7526 -210650.2878 413 52064.9884 -32429.7526 414 -204568.3343 52064.9884 415 -129873.5310 -204568.3343 416 105933.4123 -129873.5310 417 -20765.0510 105933.4123 418 -50199.8918 -20765.0510 419 -63357.6347 -50199.8918 420 -75139.5888 -63357.6347 421 -122393.5869 -75139.5888 422 -532433.2423 -122393.5869 423 -499218.3037 -532433.2423 424 -55479.1464 -499218.3037 425 -192725.5730 -55479.1464 426 -327426.9457 -192725.5730 427 -1150522.8360 -327426.9457 428 -324119.0097 -1150522.8360 429 -364142.2036 -324119.0097 430 -692656.8969 -364142.2036 431 NA -692656.8969 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 3225800.6189 1237333.7716 [2,] 1194235.0281 3225800.6189 [3,] -2198709.5189 1194235.0281 [4,] -446207.8858 -2198709.5189 [5,] -65315.8652 -446207.8858 [6,] -58061.4111 -65315.8652 [7,] 79150.5670 -58061.4111 [8,] -249447.6804 79150.5670 [9,] 742457.9366 -249447.6804 [10,] -308956.7059 742457.9366 [11,] 392076.2392 -308956.7059 [12,] 803384.0750 392076.2392 [13,] -1209577.7725 803384.0750 [14,] 35999.2506 -1209577.7725 [15,] 475470.1181 35999.2506 [16,] -267414.2259 475470.1181 [17,] 17699.4160 -267414.2259 [18,] -274855.9885 17699.4160 [19,] -346491.1382 -274855.9885 [20,] 190089.0283 -346491.1382 [21,] -245640.1815 190089.0283 [22,] -148253.8127 -245640.1815 [23,] -401854.6120 -148253.8127 [24,] -302236.3074 -401854.6120 [25,] -152079.9239 -302236.3074 [26,] -370608.0313 -152079.9239 [27,] -297472.3304 -370608.0313 [28,] 319888.6245 -297472.3304 [29,] -107458.5634 319888.6245 [30,] 264934.4893 -107458.5634 [31,] -118886.7400 264934.4893 [32,] 228735.6187 -118886.7400 [33,] 117862.7404 228735.6187 [34,] -176569.1359 117862.7404 [35,] -482482.7536 -176569.1359 [36,] 29495.2712 -482482.7536 [37,] 185164.6824 29495.2712 [38,] -31633.2784 185164.6824 [39,] -175484.6937 -31633.2784 [40,] -28586.2323 -175484.6937 [41,] -395020.7371 -28586.2323 [42,] 109887.2197 -395020.7371 [43,] -178865.3985 109887.2197 [44,] -24526.8599 -178865.3985 [45,] -86659.7032 -24526.8599 [46,] -269247.8351 -86659.7032 [47,] -97824.0452 -269247.8351 [48,] -47758.9998 -97824.0452 [49,] -173126.2485 -47758.9998 [50,] 103613.4092 -173126.2485 [51,] 25139.9739 103613.4092 [52,] -218762.7914 25139.9739 [53,] 106580.6851 -218762.7914 [54,] -315223.4826 106580.6851 [55,] 91349.0492 -315223.4826 [56,] 283612.8987 91349.0492 [57,] -23682.0837 283612.8987 [58,] -12429.1879 -23682.0837 [59,] -186594.5055 -12429.1879 [60,] -3288.1749 -186594.5055 [61,] -164170.1642 -3288.1749 [62,] -14535.6869 -164170.1642 [63,] -152468.2407 -14535.6869 [64,] 29170.0161 -152468.2407 [65,] 658.3290 29170.0161 [66,] -56301.8915 658.3290 [67,] -3474.4077 -56301.8915 [68,] -51538.3527 -3474.4077 [69,] -31098.5552 -51538.3527 [70,] -80208.5263 -31098.5552 [71,] 6724.8692 -80208.5263 [72,] -91727.8629 6724.8692 [73,] -4435.9069 -91727.8629 [74,] -47964.1058 -4435.9069 [75,] -112625.8705 -47964.1058 [76,] -32306.1501 -112625.8705 [77,] -77635.7304 -32306.1501 [78,] 22864.3015 -77635.7304 [79,] -6161.5256 22864.3015 [80,] 35011.7284 -6161.5256 [81,] -23890.0975 35011.7284 [82,] -86997.5856 -23890.0975 [83,] 46042.0802 -86997.5856 [84,] 3276.6093 46042.0802 [85,] -72852.1034 3276.6093 [86,] 43537.4967 -72852.1034 [87,] -23669.5973 43537.4967 [88,] -101295.9632 -23669.5973 [89,] 7517.6045 -101295.9632 [90,] 3632.1632 7517.6045 [91,] -18907.1836 3632.1632 [92,] 138110.3578 -18907.1836 [93,] -214809.1834 138110.3578 [94,] 85549.2598 -214809.1834 [95,] 22161.4187 85549.2598 [96,] 38786.3001 22161.4187 [97,] 68328.4234 38786.3001 [98,] -84437.6455 68328.4234 [99,] -25122.8216 -84437.6455 [100,] -28698.3104 -25122.8216 [101,] -125369.5171 -28698.3104 [102,] -146117.0270 -125369.5171 [103,] -82896.3085 -146117.0270 [104,] 36828.6130 -82896.3085 [105,] -45990.5534 36828.6130 [106,] -225506.3353 -45990.5534 [107,] 102172.4638 -225506.3353 [108,] -187517.4506 102172.4638 [109,] -99005.2886 -187517.4506 [110,] 97583.6413 -99005.2886 [111,] -181421.3633 97583.6413 [112,] -165268.3173 -181421.3633 [113,] -84430.1190 -165268.3173 [114,] 65432.0980 -84430.1190 [115,] 41884.9294 65432.0980 [116,] 1804.1853 41884.9294 [117,] 47384.8091 1804.1853 [118,] -160370.2516 47384.8091 [119,] 36507.6963 -160370.2516 [120,] -286240.9429 36507.6963 [121,] 20225.1998 -286240.9429 [122,] 5840.5691 20225.1998 [123,] 67611.1743 5840.5691 [124,] 3093.7535 67611.1743 [125,] -183753.0522 3093.7535 [126,] 34580.2165 -183753.0522 [127,] 33587.7584 34580.2165 [128,] 41443.9094 33587.7584 [129,] 34580.2165 41443.9094 [130,] 35945.7172 34580.2165 [131,] 77374.3357 35945.7172 [132,] 34580.2165 77374.3357 [133,] 31630.6387 34580.2165 [134,] 13869.3958 31630.6387 [135,] 28761.9295 13869.3958 [136,] 39448.6018 28761.9295 [137,] 34745.7604 39448.6018 [138,] 9116.6456 34745.7604 [139,] 34580.2165 9116.6456 [140,] 42113.4674 34580.2165 [141,] 34580.2165 42113.4674 [142,] 1596.0722 34580.2165 [143,] 34815.4903 1596.0722 [144,] 35943.7690 34815.4903 [145,] 35592.4067 35943.7690 [146,] 34580.2165 35592.4067 [147,] 36113.5379 34580.2165 [148,] 37095.3019 36113.5379 [149,] 34467.9757 37095.3019 [150,] 3373.2999 34467.9757 [151,] 2061.9484 3373.2999 [152,] 34580.2165 2061.9484 [153,] 11810.4517 34580.2165 [154,] 34620.4264 11810.4517 [155,] 34580.2165 34620.4264 [156,] 34580.2165 34580.2165 [157,] 30753.1273 34580.2165 [158,] 46007.7895 30753.1273 [159,] 34580.2165 46007.7895 [160,] 34580.2165 34580.2165 [161,] 34580.2165 34580.2165 [162,] 34580.2165 34580.2165 [163,] 21622.8177 34580.2165 [164,] 23262.3055 21622.8177 [165,] 33880.2124 23262.3055 [166,] 34580.2165 33880.2124 [167,] 116864.3690 34580.2165 [168,] 79550.9637 116864.3690 [169,] 34580.2165 79550.9637 [170,] 34580.2165 34580.2165 [171,] 30457.8195 34580.2165 [172,] 34580.2165 30457.8195 [173,] 23886.1881 34580.2165 [174,] 34580.2165 23886.1881 [175,] 34535.5262 34580.2165 [176,] 81914.1147 34535.5262 [177,] 33592.1205 81914.1147 [178,] 34580.2165 33592.1205 [179,] 34580.2165 34580.2165 [180,] 41164.1864 34580.2165 [181,] 15473.1509 41164.1864 [182,] -7722.3720 15473.1509 [183,] 3954.9700 -7722.3720 [184,] 31665.6161 3954.9700 [185,] 39556.5727 31665.6161 [186,] 5452.7572 39556.5727 [187,] 27424.4312 5452.7572 [188,] -3849.3352 27424.4312 [189,] 34580.2165 -3849.3352 [190,] 33720.4193 34580.2165 [191,] 34580.2165 33720.4193 [192,] 30926.1365 34580.2165 [193,] 34580.2165 30926.1365 [194,] 34513.8910 34580.2165 [195,] 34580.2165 34513.8910 [196,] -89970.4929 34580.2165 [197,] 34580.2165 -89970.4929 [198,] 34580.2165 34580.2165 [199,] -55995.9889 34580.2165 [200,] 34580.2165 -55995.9889 [201,] 34580.2165 34580.2165 [202,] 34580.2165 34580.2165 [203,] 44806.6942 34580.2165 [204,] -16733.3662 44806.6942 [205,] 33990.8040 -16733.3662 [206,] 27556.9989 33990.8040 [207,] 34580.2165 27556.9989 [208,] 65941.5748 34580.2165 [209,] 30707.7678 65941.5748 [210,] 34580.2165 30707.7678 [211,] 34580.2165 34580.2165 [212,] 60324.2404 34580.2165 [213,] -515.5477 60324.2404 [214,] 6631.3916 -515.5477 [215,] 34580.2165 6631.3916 [216,] 37450.9613 34580.2165 [217,] 36453.9569 37450.9613 [218,] 34580.2165 36453.9569 [219,] 6047.9156 34580.2165 [220,] 9219.4592 6047.9156 [221,] 34580.2165 9219.4592 [222,] 35317.0126 34580.2165 [223,] 25760.2266 35317.0126 [224,] 34812.8986 25760.2266 [225,] 34899.0433 34812.8986 [226,] 37577.8514 34899.0433 [227,] 34580.2165 37577.8514 [228,] 18260.4102 34580.2165 [229,] 34580.2165 18260.4102 [230,] 5601.5458 34580.2165 [231,] 22391.9327 5601.5458 [232,] 34580.2165 22391.9327 [233,] 16290.7635 34580.2165 [234,] 34580.2165 16290.7635 [235,] 34985.5410 34580.2165 [236,] 34580.2165 34985.5410 [237,] 36566.3486 34580.2165 [238,] 34580.2165 36566.3486 [239,] 34580.2165 34580.2165 [240,] 34580.2165 34580.2165 [241,] 29189.2405 34580.2165 [242,] 34580.2165 29189.2405 [243,] 39155.2379 34580.2165 [244,] -48965.4619 39155.2379 [245,] 34580.2165 -48965.4619 [246,] 34580.2165 34580.2165 [247,] 26517.1760 34580.2165 [248,] 34580.2165 26517.1760 [249,] 34580.2165 34580.2165 [250,] 28561.0884 34580.2165 [251,] 34580.2165 28561.0884 [252,] 15249.6882 34580.2165 [253,] 9795.8940 15249.6882 [254,] -103337.7607 9795.8940 [255,] 34580.2165 -103337.7607 [256,] 32804.4768 34580.2165 [257,] 46896.8111 32804.4768 [258,] 34580.2165 46896.8111 [259,] 34580.2165 34580.2165 [260,] 34580.2165 34580.2165 [261,] 34580.2165 34580.2165 [262,] 36942.0807 34580.2165 [263,] 51705.1001 36942.0807 [264,] 34580.2165 51705.1001 [265,] 34580.2165 34580.2165 [266,] 34580.2165 34580.2165 [267,] 65291.3572 34580.2165 [268,] 34580.2165 65291.3572 [269,] 30693.4734 34580.2165 [270,] 34580.2165 30693.4734 [271,] 34580.2165 34580.2165 [272,] 34580.2165 34580.2165 [273,] 30334.8589 34580.2165 [274,] 34580.2165 30334.8589 [275,] 39394.4755 34580.2165 [276,] 35947.1018 39394.4755 [277,] 34580.2165 35947.1018 [278,] 5052.6817 34580.2165 [279,] 21062.1530 5052.6817 [280,] 34580.2165 21062.1530 [281,] 34580.2165 34580.2165 [282,] 62309.8561 34580.2165 [283,] 34580.2165 62309.8561 [284,] 25566.3500 34580.2165 [285,] 34580.2165 25566.3500 [286,] 22734.3419 34580.2165 [287,] 35362.3736 22734.3419 [288,] 34580.2165 35362.3736 [289,] 42333.9750 34580.2165 [290,] 34580.2165 42333.9750 [291,] 44873.1463 34580.2165 [292,] 39192.5070 44873.1463 [293,] 34580.2165 39192.5070 [294,] 28105.2768 34580.2165 [295,] 59868.2221 28105.2768 [296,] 34580.2165 59868.2221 [297,] 34580.2165 34580.2165 [298,] 34580.2165 34580.2165 [299,] 34580.2165 34580.2165 [300,] 32922.8350 34580.2165 [301,] 16298.6976 32922.8350 [302,] 34580.2165 16298.6976 [303,] 34580.2165 34580.2165 [304,] 71866.7354 34580.2165 [305,] 38145.1602 71866.7354 [306,] 34580.2165 38145.1602 [307,] 28217.9952 34580.2165 [308,] 38446.2874 28217.9952 [309,] 34580.2165 38446.2874 [310,] 34580.2165 34580.2165 [311,] 35442.7760 34580.2165 [312,] 34580.2165 35442.7760 [313,] 36216.6003 34580.2165 [314,] 36302.1434 36216.6003 [315,] 49761.1157 36302.1434 [316,] 35069.1596 49761.1157 [317,] 11239.9555 35069.1596 [318,] 36971.5570 11239.9555 [319,] 34580.2165 36971.5570 [320,] 34580.2165 34580.2165 [321,] 34580.2165 34580.2165 [322,] 34580.2165 34580.2165 [323,] 34580.2165 34580.2165 [324,] 34580.2165 34580.2165 [325,] 31844.6853 34580.2165 [326,] 34580.2165 31844.6853 [327,] 34580.2165 34580.2165 [328,] 34580.2165 34580.2165 [329,] 34580.2165 34580.2165 [330,] 34580.2165 34580.2165 [331,] 34580.2165 34580.2165 [332,] 34580.2165 34580.2165 [333,] 32742.5135 34580.2165 [334,] 35353.2460 32742.5135 [335,] 27135.7859 35353.2460 [336,] 21063.1600 27135.7859 [337,] 34580.2165 21063.1600 [338,] 32349.5558 34580.2165 [339,] 40364.8296 32349.5558 [340,] 7192.9443 40364.8296 [341,] 34580.2165 7192.9443 [342,] 34580.2165 34580.2165 [343,] 34580.2165 34580.2165 [344,] 28170.1365 34580.2165 [345,] 34580.2165 28170.1365 [346,] 24924.7505 34580.2165 [347,] 34580.2165 24924.7505 [348,] 34580.2165 34580.2165 [349,] -2597.6032 34580.2165 [350,] 38401.0084 -2597.6032 [351,] 34679.5158 38401.0084 [352,] 34322.9287 34679.5158 [353,] 27462.3996 34322.9287 [354,] 34580.2165 27462.3996 [355,] 34580.2165 34580.2165 [356,] 11006.8190 34580.2165 [357,] 34580.2165 11006.8190 [358,] 34580.2165 34580.2165 [359,] 34580.2165 34580.2165 [360,] 36182.4019 34580.2165 [361,] 34580.2165 36182.4019 [362,] 34580.2165 34580.2165 [363,] -17619.3897 34580.2165 [364,] 36785.0008 -17619.3897 [365,] 34580.2165 36785.0008 [366,] 34580.2165 34580.2165 [367,] 34580.2165 34580.2165 [368,] 34580.2165 34580.2165 [369,] 36682.7327 34580.2165 [370,] 34580.2165 36682.7327 [371,] 47127.9780 34580.2165 [372,] 20788.8182 47127.9780 [373,] 13586.0443 20788.8182 [374,] 34580.2165 13586.0443 [375,] 34580.2165 34580.2165 [376,] 29720.9807 34580.2165 [377,] 28274.6000 29720.9807 [378,] -21901.6234 28274.6000 [379,] 34580.2165 -21901.6234 [380,] 34580.2165 34580.2165 [381,] 19370.4612 34580.2165 [382,] 53257.8226 19370.4612 [383,] 34580.2165 53257.8226 [384,] 36693.1984 34580.2165 [385,] -115636.0426 36693.1984 [386,] 46871.0853 -115636.0426 [387,] 7872.9811 46871.0853 [388,] 33099.2810 7872.9811 [389,] 20654.4427 33099.2810 [390,] 29633.8067 20654.4427 [391,] 3808.4378 29633.8067 [392,] -33887.7287 3808.4378 [393,] 54504.4319 -33887.7287 [394,] -141011.2663 54504.4319 [395,] -4353.0222 -141011.2663 [396,] -152446.6167 -4353.0222 [397,] 41068.9784 -152446.6167 [398,] -3133.3067 41068.9784 [399,] 21066.1182 -3133.3067 [400,] -25877.9536 21066.1182 [401,] 25950.4200 -25877.9536 [402,] 4011.8418 25950.4200 [403,] -54926.7608 4011.8418 [404,] -26003.2963 -54926.7608 [405,] -4089.2056 -26003.2963 [406,] 9440.9762 -4089.2056 [407,] 1995.5438 9440.9762 [408,] 13354.0051 1995.5438 [409,] -88117.0685 13354.0051 [410,] -4236.9658 -88117.0685 [411,] -210650.2878 -4236.9658 [412,] -32429.7526 -210650.2878 [413,] 52064.9884 -32429.7526 [414,] -204568.3343 52064.9884 [415,] -129873.5310 -204568.3343 [416,] 105933.4123 -129873.5310 [417,] -20765.0510 105933.4123 [418,] -50199.8918 -20765.0510 [419,] -63357.6347 -50199.8918 [420,] -75139.5888 -63357.6347 [421,] -122393.5869 -75139.5888 [422,] -532433.2423 -122393.5869 [423,] -499218.3037 -532433.2423 [424,] -55479.1464 -499218.3037 [425,] -192725.5730 -55479.1464 [426,] -327426.9457 -192725.5730 [427,] -1150522.8360 -327426.9457 [428,] -324119.0097 -1150522.8360 [429,] -364142.2036 -324119.0097 [430,] -692656.8969 -364142.2036 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 3225800.6189 1237333.7716 2 1194235.0281 3225800.6189 3 -2198709.5189 1194235.0281 4 -446207.8858 -2198709.5189 5 -65315.8652 -446207.8858 6 -58061.4111 -65315.8652 7 79150.5670 -58061.4111 8 -249447.6804 79150.5670 9 742457.9366 -249447.6804 10 -308956.7059 742457.9366 11 392076.2392 -308956.7059 12 803384.0750 392076.2392 13 -1209577.7725 803384.0750 14 35999.2506 -1209577.7725 15 475470.1181 35999.2506 16 -267414.2259 475470.1181 17 17699.4160 -267414.2259 18 -274855.9885 17699.4160 19 -346491.1382 -274855.9885 20 190089.0283 -346491.1382 21 -245640.1815 190089.0283 22 -148253.8127 -245640.1815 23 -401854.6120 -148253.8127 24 -302236.3074 -401854.6120 25 -152079.9239 -302236.3074 26 -370608.0313 -152079.9239 27 -297472.3304 -370608.0313 28 319888.6245 -297472.3304 29 -107458.5634 319888.6245 30 264934.4893 -107458.5634 31 -118886.7400 264934.4893 32 228735.6187 -118886.7400 33 117862.7404 228735.6187 34 -176569.1359 117862.7404 35 -482482.7536 -176569.1359 36 29495.2712 -482482.7536 37 185164.6824 29495.2712 38 -31633.2784 185164.6824 39 -175484.6937 -31633.2784 40 -28586.2323 -175484.6937 41 -395020.7371 -28586.2323 42 109887.2197 -395020.7371 43 -178865.3985 109887.2197 44 -24526.8599 -178865.3985 45 -86659.7032 -24526.8599 46 -269247.8351 -86659.7032 47 -97824.0452 -269247.8351 48 -47758.9998 -97824.0452 49 -173126.2485 -47758.9998 50 103613.4092 -173126.2485 51 25139.9739 103613.4092 52 -218762.7914 25139.9739 53 106580.6851 -218762.7914 54 -315223.4826 106580.6851 55 91349.0492 -315223.4826 56 283612.8987 91349.0492 57 -23682.0837 283612.8987 58 -12429.1879 -23682.0837 59 -186594.5055 -12429.1879 60 -3288.1749 -186594.5055 61 -164170.1642 -3288.1749 62 -14535.6869 -164170.1642 63 -152468.2407 -14535.6869 64 29170.0161 -152468.2407 65 658.3290 29170.0161 66 -56301.8915 658.3290 67 -3474.4077 -56301.8915 68 -51538.3527 -3474.4077 69 -31098.5552 -51538.3527 70 -80208.5263 -31098.5552 71 6724.8692 -80208.5263 72 -91727.8629 6724.8692 73 -4435.9069 -91727.8629 74 -47964.1058 -4435.9069 75 -112625.8705 -47964.1058 76 -32306.1501 -112625.8705 77 -77635.7304 -32306.1501 78 22864.3015 -77635.7304 79 -6161.5256 22864.3015 80 35011.7284 -6161.5256 81 -23890.0975 35011.7284 82 -86997.5856 -23890.0975 83 46042.0802 -86997.5856 84 3276.6093 46042.0802 85 -72852.1034 3276.6093 86 43537.4967 -72852.1034 87 -23669.5973 43537.4967 88 -101295.9632 -23669.5973 89 7517.6045 -101295.9632 90 3632.1632 7517.6045 91 -18907.1836 3632.1632 92 138110.3578 -18907.1836 93 -214809.1834 138110.3578 94 85549.2598 -214809.1834 95 22161.4187 85549.2598 96 38786.3001 22161.4187 97 68328.4234 38786.3001 98 -84437.6455 68328.4234 99 -25122.8216 -84437.6455 100 -28698.3104 -25122.8216 101 -125369.5171 -28698.3104 102 -146117.0270 -125369.5171 103 -82896.3085 -146117.0270 104 36828.6130 -82896.3085 105 -45990.5534 36828.6130 106 -225506.3353 -45990.5534 107 102172.4638 -225506.3353 108 -187517.4506 102172.4638 109 -99005.2886 -187517.4506 110 97583.6413 -99005.2886 111 -181421.3633 97583.6413 112 -165268.3173 -181421.3633 113 -84430.1190 -165268.3173 114 65432.0980 -84430.1190 115 41884.9294 65432.0980 116 1804.1853 41884.9294 117 47384.8091 1804.1853 118 -160370.2516 47384.8091 119 36507.6963 -160370.2516 120 -286240.9429 36507.6963 121 20225.1998 -286240.9429 122 5840.5691 20225.1998 123 67611.1743 5840.5691 124 3093.7535 67611.1743 125 -183753.0522 3093.7535 126 34580.2165 -183753.0522 127 33587.7584 34580.2165 128 41443.9094 33587.7584 129 34580.2165 41443.9094 130 35945.7172 34580.2165 131 77374.3357 35945.7172 132 34580.2165 77374.3357 133 31630.6387 34580.2165 134 13869.3958 31630.6387 135 28761.9295 13869.3958 136 39448.6018 28761.9295 137 34745.7604 39448.6018 138 9116.6456 34745.7604 139 34580.2165 9116.6456 140 42113.4674 34580.2165 141 34580.2165 42113.4674 142 1596.0722 34580.2165 143 34815.4903 1596.0722 144 35943.7690 34815.4903 145 35592.4067 35943.7690 146 34580.2165 35592.4067 147 36113.5379 34580.2165 148 37095.3019 36113.5379 149 34467.9757 37095.3019 150 3373.2999 34467.9757 151 2061.9484 3373.2999 152 34580.2165 2061.9484 153 11810.4517 34580.2165 154 34620.4264 11810.4517 155 34580.2165 34620.4264 156 34580.2165 34580.2165 157 30753.1273 34580.2165 158 46007.7895 30753.1273 159 34580.2165 46007.7895 160 34580.2165 34580.2165 161 34580.2165 34580.2165 162 34580.2165 34580.2165 163 21622.8177 34580.2165 164 23262.3055 21622.8177 165 33880.2124 23262.3055 166 34580.2165 33880.2124 167 116864.3690 34580.2165 168 79550.9637 116864.3690 169 34580.2165 79550.9637 170 34580.2165 34580.2165 171 30457.8195 34580.2165 172 34580.2165 30457.8195 173 23886.1881 34580.2165 174 34580.2165 23886.1881 175 34535.5262 34580.2165 176 81914.1147 34535.5262 177 33592.1205 81914.1147 178 34580.2165 33592.1205 179 34580.2165 34580.2165 180 41164.1864 34580.2165 181 15473.1509 41164.1864 182 -7722.3720 15473.1509 183 3954.9700 -7722.3720 184 31665.6161 3954.9700 185 39556.5727 31665.6161 186 5452.7572 39556.5727 187 27424.4312 5452.7572 188 -3849.3352 27424.4312 189 34580.2165 -3849.3352 190 33720.4193 34580.2165 191 34580.2165 33720.4193 192 30926.1365 34580.2165 193 34580.2165 30926.1365 194 34513.8910 34580.2165 195 34580.2165 34513.8910 196 -89970.4929 34580.2165 197 34580.2165 -89970.4929 198 34580.2165 34580.2165 199 -55995.9889 34580.2165 200 34580.2165 -55995.9889 201 34580.2165 34580.2165 202 34580.2165 34580.2165 203 44806.6942 34580.2165 204 -16733.3662 44806.6942 205 33990.8040 -16733.3662 206 27556.9989 33990.8040 207 34580.2165 27556.9989 208 65941.5748 34580.2165 209 30707.7678 65941.5748 210 34580.2165 30707.7678 211 34580.2165 34580.2165 212 60324.2404 34580.2165 213 -515.5477 60324.2404 214 6631.3916 -515.5477 215 34580.2165 6631.3916 216 37450.9613 34580.2165 217 36453.9569 37450.9613 218 34580.2165 36453.9569 219 6047.9156 34580.2165 220 9219.4592 6047.9156 221 34580.2165 9219.4592 222 35317.0126 34580.2165 223 25760.2266 35317.0126 224 34812.8986 25760.2266 225 34899.0433 34812.8986 226 37577.8514 34899.0433 227 34580.2165 37577.8514 228 18260.4102 34580.2165 229 34580.2165 18260.4102 230 5601.5458 34580.2165 231 22391.9327 5601.5458 232 34580.2165 22391.9327 233 16290.7635 34580.2165 234 34580.2165 16290.7635 235 34985.5410 34580.2165 236 34580.2165 34985.5410 237 36566.3486 34580.2165 238 34580.2165 36566.3486 239 34580.2165 34580.2165 240 34580.2165 34580.2165 241 29189.2405 34580.2165 242 34580.2165 29189.2405 243 39155.2379 34580.2165 244 -48965.4619 39155.2379 245 34580.2165 -48965.4619 246 34580.2165 34580.2165 247 26517.1760 34580.2165 248 34580.2165 26517.1760 249 34580.2165 34580.2165 250 28561.0884 34580.2165 251 34580.2165 28561.0884 252 15249.6882 34580.2165 253 9795.8940 15249.6882 254 -103337.7607 9795.8940 255 34580.2165 -103337.7607 256 32804.4768 34580.2165 257 46896.8111 32804.4768 258 34580.2165 46896.8111 259 34580.2165 34580.2165 260 34580.2165 34580.2165 261 34580.2165 34580.2165 262 36942.0807 34580.2165 263 51705.1001 36942.0807 264 34580.2165 51705.1001 265 34580.2165 34580.2165 266 34580.2165 34580.2165 267 65291.3572 34580.2165 268 34580.2165 65291.3572 269 30693.4734 34580.2165 270 34580.2165 30693.4734 271 34580.2165 34580.2165 272 34580.2165 34580.2165 273 30334.8589 34580.2165 274 34580.2165 30334.8589 275 39394.4755 34580.2165 276 35947.1018 39394.4755 277 34580.2165 35947.1018 278 5052.6817 34580.2165 279 21062.1530 5052.6817 280 34580.2165 21062.1530 281 34580.2165 34580.2165 282 62309.8561 34580.2165 283 34580.2165 62309.8561 284 25566.3500 34580.2165 285 34580.2165 25566.3500 286 22734.3419 34580.2165 287 35362.3736 22734.3419 288 34580.2165 35362.3736 289 42333.9750 34580.2165 290 34580.2165 42333.9750 291 44873.1463 34580.2165 292 39192.5070 44873.1463 293 34580.2165 39192.5070 294 28105.2768 34580.2165 295 59868.2221 28105.2768 296 34580.2165 59868.2221 297 34580.2165 34580.2165 298 34580.2165 34580.2165 299 34580.2165 34580.2165 300 32922.8350 34580.2165 301 16298.6976 32922.8350 302 34580.2165 16298.6976 303 34580.2165 34580.2165 304 71866.7354 34580.2165 305 38145.1602 71866.7354 306 34580.2165 38145.1602 307 28217.9952 34580.2165 308 38446.2874 28217.9952 309 34580.2165 38446.2874 310 34580.2165 34580.2165 311 35442.7760 34580.2165 312 34580.2165 35442.7760 313 36216.6003 34580.2165 314 36302.1434 36216.6003 315 49761.1157 36302.1434 316 35069.1596 49761.1157 317 11239.9555 35069.1596 318 36971.5570 11239.9555 319 34580.2165 36971.5570 320 34580.2165 34580.2165 321 34580.2165 34580.2165 322 34580.2165 34580.2165 323 34580.2165 34580.2165 324 34580.2165 34580.2165 325 31844.6853 34580.2165 326 34580.2165 31844.6853 327 34580.2165 34580.2165 328 34580.2165 34580.2165 329 34580.2165 34580.2165 330 34580.2165 34580.2165 331 34580.2165 34580.2165 332 34580.2165 34580.2165 333 32742.5135 34580.2165 334 35353.2460 32742.5135 335 27135.7859 35353.2460 336 21063.1600 27135.7859 337 34580.2165 21063.1600 338 32349.5558 34580.2165 339 40364.8296 32349.5558 340 7192.9443 40364.8296 341 34580.2165 7192.9443 342 34580.2165 34580.2165 343 34580.2165 34580.2165 344 28170.1365 34580.2165 345 34580.2165 28170.1365 346 24924.7505 34580.2165 347 34580.2165 24924.7505 348 34580.2165 34580.2165 349 -2597.6032 34580.2165 350 38401.0084 -2597.6032 351 34679.5158 38401.0084 352 34322.9287 34679.5158 353 27462.3996 34322.9287 354 34580.2165 27462.3996 355 34580.2165 34580.2165 356 11006.8190 34580.2165 357 34580.2165 11006.8190 358 34580.2165 34580.2165 359 34580.2165 34580.2165 360 36182.4019 34580.2165 361 34580.2165 36182.4019 362 34580.2165 34580.2165 363 -17619.3897 34580.2165 364 36785.0008 -17619.3897 365 34580.2165 36785.0008 366 34580.2165 34580.2165 367 34580.2165 34580.2165 368 34580.2165 34580.2165 369 36682.7327 34580.2165 370 34580.2165 36682.7327 371 47127.9780 34580.2165 372 20788.8182 47127.9780 373 13586.0443 20788.8182 374 34580.2165 13586.0443 375 34580.2165 34580.2165 376 29720.9807 34580.2165 377 28274.6000 29720.9807 378 -21901.6234 28274.6000 379 34580.2165 -21901.6234 380 34580.2165 34580.2165 381 19370.4612 34580.2165 382 53257.8226 19370.4612 383 34580.2165 53257.8226 384 36693.1984 34580.2165 385 -115636.0426 36693.1984 386 46871.0853 -115636.0426 387 7872.9811 46871.0853 388 33099.2810 7872.9811 389 20654.4427 33099.2810 390 29633.8067 20654.4427 391 3808.4378 29633.8067 392 -33887.7287 3808.4378 393 54504.4319 -33887.7287 394 -141011.2663 54504.4319 395 -4353.0222 -141011.2663 396 -152446.6167 -4353.0222 397 41068.9784 -152446.6167 398 -3133.3067 41068.9784 399 21066.1182 -3133.3067 400 -25877.9536 21066.1182 401 25950.4200 -25877.9536 402 4011.8418 25950.4200 403 -54926.7608 4011.8418 404 -26003.2963 -54926.7608 405 -4089.2056 -26003.2963 406 9440.9762 -4089.2056 407 1995.5438 9440.9762 408 13354.0051 1995.5438 409 -88117.0685 13354.0051 410 -4236.9658 -88117.0685 411 -210650.2878 -4236.9658 412 -32429.7526 -210650.2878 413 52064.9884 -32429.7526 414 -204568.3343 52064.9884 415 -129873.5310 -204568.3343 416 105933.4123 -129873.5310 417 -20765.0510 105933.4123 418 -50199.8918 -20765.0510 419 -63357.6347 -50199.8918 420 -75139.5888 -63357.6347 421 -122393.5869 -75139.5888 422 -532433.2423 -122393.5869 423 -499218.3037 -532433.2423 424 -55479.1464 -499218.3037 425 -192725.5730 -55479.1464 426 -327426.9457 -192725.5730 427 -1150522.8360 -327426.9457 428 -324119.0097 -1150522.8360 429 -364142.2036 -324119.0097 430 -692656.8969 -364142.2036 > 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/7s46o1291211846.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/8s46o1291211846.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/9kv5r1291211846.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/10kv5r1291211846.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/116vmx1291211846.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/12rwkl1291211846.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/13yxhf1291211846.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/1496z01291211846.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/15cpfo1291211846.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/168hve1291211846.tab") + } > > try(system("convert tmp/1ec8y1291211846.ps tmp/1ec8y1291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/263811291211846.ps tmp/263811291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/363811291211846.ps tmp/363811291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/463811291211846.ps tmp/463811291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/5hvpl1291211846.ps tmp/5hvpl1291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/6hvpl1291211846.ps tmp/6hvpl1291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/7s46o1291211846.ps tmp/7s46o1291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/8s46o1291211846.ps tmp/8s46o1291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/9kv5r1291211846.ps tmp/9kv5r1291211846.png",intern=TRUE)) character(0) > try(system("convert tmp/10kv5r1291211846.ps tmp/10kv5r1291211846.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.352 3.110 13.872