R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(162556 + ,807 + ,213118 + ,6282154 + ,29790 + ,444 + ,81767 + ,4321023 + ,87550 + ,412 + ,153198 + ,4111912 + ,84738 + ,428 + ,-26007 + ,223193 + ,54660 + ,315 + ,126942 + ,1491348 + ,42634 + ,168 + ,157214 + ,1629616 + ,40949 + ,263 + ,129352 + ,1398893 + ,45187 + ,267 + ,234817 + ,1926517 + ,37704 + ,228 + ,60448 + ,983660 + ,16275 + ,129 + ,47818 + ,1443586 + ,25830 + ,104 + ,245546 + ,1073089 + ,12679 + ,122 + ,48020 + ,984885 + ,18014 + ,393 + ,-1710 + ,1405225 + ,43556 + ,190 + ,32648 + ,227132 + ,24811 + ,280 + ,95350 + ,929118 + ,6575 + ,63 + ,151352 + ,1071292 + ,7123 + ,102 + ,288170 + ,638830 + ,21950 + ,265 + ,114337 + ,856956 + ,37597 + ,234 + ,37884 + ,992426 + ,17821 + ,277 + ,122844 + ,444477 + ,12988 + ,73 + ,82340 + ,857217 + ,22330 + ,67 + ,79801 + ,711969 + ,13326 + ,103 + ,165548 + ,702380 + ,16189 + ,290 + ,116384 + ,358589 + ,7146 + ,83 + ,134028 + ,297978 + ,15824 + ,56 + ,63838 + ,585715 + ,27664 + ,236 + ,74996 + ,657954 + ,11920 + ,73 + ,31080 + ,209458 + ,8568 + ,34 + ,32168 + ,786690 + ,14416 + ,139 + ,49857 + ,439798 + ,3369 + ,26 + ,87161 + ,688779 + ,11819 + ,70 + ,106113 + ,574339 + ,6984 + ,40 + ,80570 + ,741409 + ,4519 + ,42 + ,102129 + ,597793 + ,2220 + ,12 + ,301670 + ,644190 + ,18562 + ,211 + ,102313 + ,377934 + ,10327 + ,74 + ,88577 + ,640273 + ,5336 + ,80 + ,112477 + ,697458 + ,2365 + ,83 + ,191778 + ,550608 + ,4069 + ,131 + ,79804 + ,207393 + ,8636 + ,203 + ,128294 + ,301607 + ,13718 + ,56 + ,96448 + ,345783 + ,4525 + ,89 + ,93811 + ,501749 + ,6869 + ,88 + ,117520 + ,379983 + ,4628 + ,39 + ,69159 + ,387475 + ,3689 + ,25 + ,101792 + ,377305 + ,4891 + ,49 + ,210568 + ,370837 + ,7489 + ,149 + ,136996 + ,430866 + ,4901 + ,58 + ,121920 + ,469107 + ,2284 + ,41 + ,76403 + ,194493 + ,3160 + ,90 + ,108094 + ,530670 + ,4150 + ,136 + ,134759 + ,518365 + ,7285 + ,97 + ,188873 + ,491303 + ,1134 + ,63 + ,146216 + ,527021 + ,4658 + ,114 + ,156608 + ,233773 + ,2384 + ,77 + ,61348 + ,405972 + ,3748 + ,6 + ,50350 + ,652925 + ,5371 + ,47 + ,87720 + ,446211 + ,1285 + ,51 + ,99489 + ,341340 + ,9327 + ,85 + ,87419 + ,387699 + ,5565 + ,43 + ,94355 + ,493408 + ,1528 + ,32 + ,60326 + ,146494 + ,3122 + ,25 + ,94670 + ,414462 + ,7561 + ,77 + ,82425 + ,364304 + ,2675 + ,54 + ,59017 + ,355178 + ,13253 + ,251 + ,90829 + ,357760 + ,880 + ,15 + ,80791 + ,261216 + ,2053 + ,44 + ,100423 + ,397144 + ,1424 + ,73 + ,131116 + ,374943 + ,4036 + ,85 + ,100269 + ,424898 + ,3045 + ,49 + ,27330 + ,202055 + ,5119 + ,38 + ,39039 + ,378525 + ,1431 + ,35 + ,106885 + ,310768 + ,554 + ,9 + ,79285 + ,325738 + ,1975 + ,34 + ,118881 + ,394510 + ,1765 + ,20 + ,77623 + ,247060 + ,1012 + ,29 + ,114768 + ,368078 + ,810 + ,11 + ,74015 + ,236761 + ,1280 + ,52 + ,69465 + ,312378 + ,666 + ,13 + ,117869 + ,339836 + ,1380 + ,29 + ,60982 + ,347385 + ,4677 + ,66 + ,90131 + ,426280 + ,876 + ,33 + ,138971 + ,352850 + ,814 + ,15 + ,39625 + ,301881 + ,514 + ,15 + ,102725 + ,377516 + ,5692 + ,68 + ,64239 + ,357312 + ,3642 + ,100 + ,90262 + ,458343 + ,540 + ,13 + ,103960 + ,354228 + ,2099 + ,45 + ,106611 + ,308636 + ,567 + ,14 + ,103345 + ,386212 + ,2001 + ,36 + ,95551 + ,393343 + ,2949 + ,40 + ,82903 + ,378509 + ,2253 + ,68 + ,63593 + ,452469 + ,6533 + ,29 + ,126910 + ,364839 + ,1889 + ,43 + ,37527 + ,358649 + ,3055 + ,30 + ,60247 + ,376641 + ,272 + ,9 + ,112995 + ,429112 + ,1414 + ,22 + ,70184 + ,330546 + ,2564 + ,19 + ,130140 + ,403560 + ,1383 + ,9 + ,73221 + ,317892 + ,1261 + ,31 + ,76114 + ,307528 + ,975 + ,19 + ,90534 + ,235133 + ,3366 + ,55 + ,108479 + ,299243 + ,576 + ,8 + ,113761 + ,314073 + ,1686 + ,28 + ,68696 + ,368186 + ,746 + ,29 + ,71561 + ,269661 + ,3192 + ,48 + ,59831 + ,125390 + ,2045 + ,16 + ,97890 + ,510834 + ,5702 + ,47 + ,101481 + ,321896 + ,1932 + ,20 + ,72954 + ,249898 + ,936 + ,22 + ,67939 + ,408881 + ,3437 + ,33 + ,48022 + ,158492 + ,5131 + ,44 + ,86111 + ,292154 + ,2397 + ,13 + ,74020 + ,289513 + ,1389 + ,6 + ,57530 + ,378049 + ,1503 + ,35 + ,56364 + ,343466 + ,402 + ,8 + ,84990 + ,332743 + ,2239 + ,17 + ,88590 + ,442882 + ,2234 + ,11 + ,77200 + ,214215 + ,837 + ,21 + ,61262 + ,315688 + ,10579 + ,92 + ,110309 + ,375195 + ,875 + ,12 + ,67000 + ,334280 + ,1585 + ,112 + ,93099 + ,355864 + ,1659 + ,25 + ,107577 + ,480382 + ,2647 + ,17 + ,62920 + ,353058 + ,3294 + ,23 + ,75832 + ,217193 + ,0 + ,0 + ,60720 + ,315380 + ,94 + ,10 + ,60793 + ,314533 + ,422 + ,23 + ,57935 + ,318056 + ,0 + ,0 + ,60720 + ,315380 + ,34 + ,7 + ,60630 + ,314353 + ,1558 + ,25 + ,55637 + ,369448 + ,0 + ,1 + ,60720 + ,315380 + ,43 + ,20 + ,60887 + ,312846 + ,645 + ,4 + ,60720 + ,312075 + ,316 + ,4 + ,60505 + ,315009 + ,115 + ,10 + ,60945 + ,318903 + ,5 + ,1 + ,60720 + ,314887 + ,897 + ,4 + ,60720 + ,314913 + ,0 + ,0 + ,60720 + ,315380 + ,389 + ,8 + ,58990 + ,325506 + ,0 + ,0 + ,60720 + ,315380 + ,1002 + ,11 + ,56750 + ,298568 + ,36 + ,4 + ,60894 + ,315834 + ,460 + ,15 + ,63346 + ,329784 + ,309 + ,9 + ,56535 + ,312878 + ,0 + ,0 + ,60720 + ,315380 + ,9 + ,7 + ,60835 + ,314987 + ,271 + ,2 + ,60720 + ,325249 + ,14 + ,0 + ,61016 + ,315877 + ,520 + ,7 + ,58650 + ,291650 + ,1766 + ,46 + ,60438 + ,305959 + ,0 + ,5 + ,60720 + ,315380 + ,458 + ,7 + ,58625 + ,297765 + ,20 + ,2 + ,60938 + ,315245 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,98 + ,2 + ,61490 + ,315236 + ,405 + ,5 + ,60845 + ,336425 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,483 + ,7 + ,60830 + ,306268 + ,454 + ,24 + ,63261 + ,302187 + ,47 + ,1 + ,60720 + ,314882 + ,0 + ,0 + ,60720 + ,315380 + ,757 + ,18 + ,45689 + ,382712 + ,4655 + ,55 + ,60720 + ,341570 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,36 + ,3 + ,61564 + ,312412 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,0 + ,60720 + ,315380 + ,126 + ,8 + ,60951 + ,315547 + ,400 + ,113 + ,60720 + ,313267 + ,71 + ,0 + ,60745 + ,316176 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,972 + ,19 + ,71642 + ,359335 + ,531 + ,11 + ,71641 + ,330068 + ,2461 + ,25 + ,55792 + ,314289 + ,378 + ,16 + ,71873 + ,297413 + ,23 + ,5 + ,62555 + ,314806 + ,638 + ,11 + ,60370 + ,333210 + ,2300 + ,23 + ,64873 + ,352108 + ,149 + ,6 + ,62041 + ,313332 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,7 + ,59500 + ,318745 + ,0 + ,0 + ,60720 + ,315380 + ,141 + ,7 + ,61630 + ,315366 + ,0 + ,0 + ,60720 + ,315380 + ,28 + ,3 + ,60890 + ,315688 + ,0 + ,0 + ,60720 + ,315380 + ,4980 + ,89 + ,113521 + ,409642 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,472 + ,19 + ,80045 + ,269587 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,12 + ,50804 + ,300962 + ,496 + ,12 + ,87390 + ,325479 + ,10 + ,5 + ,61656 + ,316155 + ,63 + ,2 + ,65688 + ,318574 + ,0 + ,0 + ,60720 + ,315380 + ,1136 + ,26 + ,48522 + ,343613 + ,265 + ,3 + ,60720 + ,306948 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,11 + ,57640 + ,330059 + ,474 + ,10 + ,61977 + ,288985 + ,534 + ,5 + ,62620 + ,304485 + ,0 + ,2 + ,60720 + ,315380 + ,15 + ,6 + ,60831 + ,315688 + ,397 + ,7 + ,60646 + ,317736 + ,0 + ,2 + ,60720 + ,315380 + ,1866 + ,28 + ,56225 + ,322331 + ,288 + ,3 + ,60510 + ,296656 + ,0 + ,0 + ,60720 + ,315380 + ,3 + ,1 + ,60698 + ,315354 + ,468 + ,20 + ,60720 + ,312161 + ,20 + ,1 + ,60805 + ,315576 + ,278 + ,22 + ,61404 + ,314922 + ,61 + ,9 + ,60720 + ,314551 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,2 + ,65276 + ,312339 + ,0 + ,0 + ,60720 + ,315380 + ,317 + ,7 + ,63915 + ,298700 + ,738 + ,9 + ,60720 + ,321376 + ,0 + ,0 + ,60720 + ,315380 + ,368 + ,13 + ,61686 + ,303230 + ,0 + ,0 + ,60720 + ,315380 + ,2 + ,0 + ,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 + ,7 + ,59818 + ,315637 + ,2332 + ,2 + ,72680 + ,324385 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,131 + ,15 + ,61808 + ,308989 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,0 + ,60720 + ,315380 + ,167 + ,1 + ,64245 + ,307322 + ,622 + ,38 + ,73007 + ,304376 + ,2328 + ,57 + ,82732 + ,253588 + ,0 + ,0 + ,60720 + ,315380 + ,365 + ,7 + ,54820 + ,309560 + ,364 + ,26 + ,47705 + ,298466 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,13 + ,72835 + ,343929 + ,307 + ,10 + ,58856 + ,331955 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,9 + ,77655 + ,381180 + ,0 + ,0 + ,60720 + ,315380 + ,138 + ,26 + ,69817 + ,331420 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,125 + ,19 + ,60798 + ,310201 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,12 + ,62452 + ,320016 + ,335 + ,23 + ,64175 + ,320398 + ,0 + ,0 + ,60720 + ,315380 + ,1324 + ,29 + ,67440 + ,291841 + ,176 + ,8 + ,68136 + ,310670 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,26 + ,56726 + ,313491 + ,0 + ,0 + ,60720 + ,315380 + ,333 + ,9 + ,70811 + ,331323 + ,0 + ,0 + ,60720 + ,315380 + ,601 + ,5 + ,60720 + ,319210 + ,30 + ,3 + ,62045 + ,318098 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,13 + ,54323 + ,292754 + ,0 + ,0 + ,60720 + ,315380 + ,165 + ,12 + ,62841 + ,325176 + ,453 + ,19 + ,81125 + ,365959 + ,0 + ,0 + ,60720 + ,315380 + ,53 + ,10 + ,59506 + ,302409 + ,382 + ,9 + ,59365 + ,340968 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,9 + ,60720 + ,315380 + ,30 + ,4 + ,60798 + ,313164 + ,290 + ,1 + ,58790 + ,301164 + ,0 + ,1 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,366 + ,14 + ,61808 + ,344425 + ,2 + ,12 + ,60735 + ,315394 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,19 + ,64016 + ,316647 + ,384 + ,17 + ,54683 + ,309836 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,365 + ,32 + ,87192 + ,346611 + ,0 + ,0 + ,60720 + ,315380 + ,49 + ,14 + ,64107 + ,322031 + ,3 + ,8 + ,60761 + ,315656 + ,133 + ,4 + ,65990 + ,339445 + ,32 + ,0 + ,59988 + ,314964 + ,368 + ,20 + ,61167 + ,297141 + ,1 + ,5 + ,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 + ,1 + ,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 + ,4 + ,60379 + ,313729 + ,1 + ,1 + ,60727 + ,315388 + ,314 + ,4 + ,60720 + ,315371 + ,844 + ,20 + ,60925 + ,296139 + ,0 + ,0 + ,60720 + ,315380 + ,26 + ,1 + ,60896 + ,313880 + ,125 + ,10 + ,59734 + ,317698 + ,304 + ,12 + ,62969 + ,295580 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,621 + ,13 + ,60720 + ,308256 + ,0 + ,0 + ,60720 + ,315380 + ,119 + ,3 + ,59118 + ,303677 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,1595 + ,10 + ,60720 + ,319369 + ,312 + ,3 + ,58598 + ,318690 + ,60 + ,7 + ,61124 + ,314049 + ,587 + ,10 + ,59595 + ,325699 + ,135 + ,1 + ,62065 + ,314210 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,514 + ,15 + ,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 + ,28 + ,61600 + ,308336 + ,180 + ,9 + ,59635 + ,316386 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,7 + ,60720 + ,315553 + ,0 + ,0 + ,60720 + ,315380 + ,448 + ,7 + ,59781 + ,323361 + ,227 + ,7 + ,76644 + ,336639 + ,174 + ,3 + ,64820 + ,307424 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,121 + ,11 + ,56178 + ,295370 + ,607 + ,7 + ,60436 + ,322340 + ,2212 + ,10 + ,60720 + ,319864 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,530 + ,18 + ,73433 + ,317291 + ,571 + ,14 + ,41477 + ,280398 + ,0 + ,0 + ,60720 + ,315380 + ,78 + ,12 + ,62700 + ,317330 + ,2489 + ,29 + ,67804 + ,238125 + ,131 + ,3 + ,59661 + ,327071 + ,923 + ,6 + ,58620 + ,309038 + ,72 + ,3 + ,60398 + ,314210 + ,572 + ,8 + ,58580 + ,307930 + ,397 + ,10 + ,62710 + ,322327 + ,450 + ,6 + ,59325 + ,292136 + ,622 + ,8 + ,60950 + ,263276 + ,694 + ,6 + ,68060 + ,367655 + ,3425 + ,9 + ,83620 + ,283910 + ,562 + ,8 + ,58456 + ,283587 + ,4917 + ,26 + ,52811 + ,243650 + ,1442 + ,239 + ,121173 + ,438493 + ,529 + ,7 + ,63870 + ,296261 + ,2126 + ,41 + ,21001 + ,230621 + ,1061 + ,3 + ,70415 + ,304252 + ,776 + ,8 + ,64230 + ,333505 + ,611 + ,6 + ,59190 + ,296919 + ,1526 + ,21 + ,69351 + ,278990 + ,592 + ,7 + ,64270 + ,276898 + ,1182 + ,11 + ,70694 + ,327007 + ,621 + ,11 + ,68005 + ,317046 + ,989 + ,12 + ,58930 + ,304555 + ,438 + ,9 + ,58320 + ,298096 + ,726 + ,3 + ,69980 + ,231861 + ,1303 + ,57 + ,69863 + ,309422 + ,7419 + ,21 + ,63255 + ,286963 + ,1164 + ,15 + ,57320 + ,269753 + ,3310 + ,32 + ,75230 + ,448243 + ,1920 + ,11 + ,79420 + ,165404 + ,965 + ,2 + ,73490 + ,204325 + ,3256 + ,23 + ,35250 + ,407159 + ,1135 + ,20 + ,62285 + ,290476 + ,1270 + ,24 + ,69206 + ,275311 + ,661 + ,1 + ,65920 + ,246541 + ,1013 + ,1 + ,69770 + ,253468 + ,2844 + ,74 + ,72683 + ,240897 + ,11528 + ,68 + ,-14545 + ,-83265 + ,6526 + ,20 + ,55830 + ,-42143 + ,2264 + ,20 + ,55174 + ,272713 + ,5109 + ,82 + ,67038 + ,215362 + ,3999 + ,21 + ,51252 + ,42754 + ,35624 + ,244 + ,157278 + ,306275 + ,9252 + ,32 + ,79510 + ,253537 + ,15236 + ,86 + ,77440 + ,372631 + ,18073 + ,69 + ,27284 + ,-7170) + ,dim=c(4 + ,431) + ,dimnames=list(c('Costs' + ,'Orders' + ,'Dividends' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(4,431),dimnames=list(c('Costs','Orders','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 Costs Orders Dividends 1 6282154 162556 807 213118 2 4321023 29790 444 81767 3 4111912 87550 412 153198 4 223193 84738 428 -26007 5 1491348 54660 315 126942 6 1629616 42634 168 157214 7 1398893 40949 263 129352 8 1926517 45187 267 234817 9 983660 37704 228 60448 10 1443586 16275 129 47818 11 1073089 25830 104 245546 12 984885 12679 122 48020 13 1405225 18014 393 -1710 14 227132 43556 190 32648 15 929118 24811 280 95350 16 1071292 6575 63 151352 17 638830 7123 102 288170 18 856956 21950 265 114337 19 992426 37597 234 37884 20 444477 17821 277 122844 21 857217 12988 73 82340 22 711969 22330 67 79801 23 702380 13326 103 165548 24 358589 16189 290 116384 25 297978 7146 83 134028 26 585715 15824 56 63838 27 657954 27664 236 74996 28 209458 11920 73 31080 29 786690 8568 34 32168 30 439798 14416 139 49857 31 688779 3369 26 87161 32 574339 11819 70 106113 33 741409 6984 40 80570 34 597793 4519 42 102129 35 644190 2220 12 301670 36 377934 18562 211 102313 37 640273 10327 74 88577 38 697458 5336 80 112477 39 550608 2365 83 191778 40 207393 4069 131 79804 41 301607 8636 203 128294 42 345783 13718 56 96448 43 501749 4525 89 93811 44 379983 6869 88 117520 45 387475 4628 39 69159 46 377305 3689 25 101792 47 370837 4891 49 210568 48 430866 7489 149 136996 49 469107 4901 58 121920 50 194493 2284 41 76403 51 530670 3160 90 108094 52 518365 4150 136 134759 53 491303 7285 97 188873 54 527021 1134 63 146216 55 233773 4658 114 156608 56 405972 2384 77 61348 57 652925 3748 6 50350 58 446211 5371 47 87720 59 341340 1285 51 99489 60 387699 9327 85 87419 61 493408 5565 43 94355 62 146494 1528 32 60326 63 414462 3122 25 94670 64 364304 7561 77 82425 65 355178 2675 54 59017 66 357760 13253 251 90829 67 261216 880 15 80791 68 397144 2053 44 100423 69 374943 1424 73 131116 70 424898 4036 85 100269 71 202055 3045 49 27330 72 378525 5119 38 39039 73 310768 1431 35 106885 74 325738 554 9 79285 75 394510 1975 34 118881 76 247060 1765 20 77623 77 368078 1012 29 114768 78 236761 810 11 74015 79 312378 1280 52 69465 80 339836 666 13 117869 81 347385 1380 29 60982 82 426280 4677 66 90131 83 352850 876 33 138971 84 301881 814 15 39625 85 377516 514 15 102725 86 357312 5692 68 64239 87 458343 3642 100 90262 88 354228 540 13 103960 89 308636 2099 45 106611 90 386212 567 14 103345 91 393343 2001 36 95551 92 378509 2949 40 82903 93 452469 2253 68 63593 94 364839 6533 29 126910 95 358649 1889 43 37527 96 376641 3055 30 60247 97 429112 272 9 112995 98 330546 1414 22 70184 99 403560 2564 19 130140 100 317892 1383 9 73221 101 307528 1261 31 76114 102 235133 975 19 90534 103 299243 3366 55 108479 104 314073 576 8 113761 105 368186 1686 28 68696 106 269661 746 29 71561 107 125390 3192 48 59831 108 510834 2045 16 97890 109 321896 5702 47 101481 110 249898 1932 20 72954 111 408881 936 22 67939 112 158492 3437 33 48022 113 292154 5131 44 86111 114 289513 2397 13 74020 115 378049 1389 6 57530 116 343466 1503 35 56364 117 332743 402 8 84990 118 442882 2239 17 88590 119 214215 2234 11 77200 120 315688 837 21 61262 121 375195 10579 92 110309 122 334280 875 12 67000 123 355864 1585 112 93099 124 480382 1659 25 107577 125 353058 2647 17 62920 126 217193 3294 23 75832 127 315380 0 0 60720 128 314533 94 10 60793 129 318056 422 23 57935 130 315380 0 0 60720 131 314353 34 7 60630 132 369448 1558 25 55637 133 315380 0 1 60720 134 312846 43 20 60887 135 312075 645 4 60720 136 315009 316 4 60505 137 318903 115 10 60945 138 314887 5 1 60720 139 314913 897 4 60720 140 315380 0 0 60720 141 325506 389 8 58990 142 315380 0 0 60720 143 298568 1002 11 56750 144 315834 36 4 60894 145 329784 460 15 63346 146 312878 309 9 56535 147 315380 0 0 60720 148 314987 9 7 60835 149 325249 271 2 60720 150 315877 14 0 61016 151 291650 520 7 58650 152 305959 1766 46 60438 153 315380 0 5 60720 154 297765 458 7 58625 155 315245 20 2 60938 156 315380 0 0 60720 157 315380 0 0 60720 158 315236 98 2 61490 159 336425 405 5 60845 160 315380 0 0 60720 161 315380 0 0 60720 162 315380 0 0 60720 163 315380 0 0 60720 164 306268 483 7 60830 165 302187 454 24 63261 166 314882 47 1 60720 167 315380 0 0 60720 168 382712 757 18 45689 169 341570 4655 55 60720 170 315380 0 0 60720 171 315380 0 0 60720 172 312412 36 3 61564 173 315380 0 0 60720 174 309596 203 9 61938 175 315380 0 0 60720 176 315547 126 8 60951 177 313267 400 113 60720 178 316176 71 0 60745 179 315380 0 0 60720 180 315380 0 0 60720 181 359335 972 19 71642 182 330068 531 11 71641 183 314289 2461 25 55792 184 297413 378 16 71873 185 314806 23 5 62555 186 333210 638 11 60370 187 352108 2300 23 64873 188 313332 149 6 62041 189 291787 226 5 65745 190 315380 0 0 60720 191 318745 275 7 59500 192 315380 0 0 60720 193 315366 141 7 61630 194 315380 0 0 60720 195 315688 28 3 60890 196 315380 0 0 60720 197 409642 4980 89 113521 198 315380 0 0 60720 199 315380 0 0 60720 200 269587 472 19 80045 201 315380 0 0 60720 202 315380 0 0 60720 203 315380 0 0 60720 204 300962 203 12 50804 205 325479 496 12 87390 206 316155 10 5 61656 207 318574 63 2 65688 208 315380 0 0 60720 209 343613 1136 26 48522 210 306948 265 3 60720 211 315380 0 0 60720 212 315380 0 0 60720 213 330059 267 11 57640 214 288985 474 10 61977 215 304485 534 5 62620 216 315380 0 2 60720 217 315688 15 6 60831 218 317736 397 7 60646 219 315380 0 2 60720 220 322331 1866 28 56225 221 296656 288 3 60510 222 315380 0 0 60720 223 315354 3 1 60698 224 312161 468 20 60720 225 315576 20 1 60805 226 314922 278 22 61404 227 314551 61 9 60720 228 315380 0 0 60720 229 312339 192 2 65276 230 315380 0 0 60720 231 298700 317 7 63915 232 321376 738 9 60720 233 315380 0 0 60720 234 303230 368 13 61686 235 315380 0 0 60720 236 315487 2 0 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 7 59818 245 324385 2332 2 72680 246 315380 0 0 60720 247 315380 0 0 60720 248 308989 131 15 61808 249 315380 0 0 60720 250 315380 0 0 60720 251 296702 206 9 53110 252 315380 0 0 60720 253 307322 167 1 64245 254 304376 622 38 73007 255 253588 2328 57 82732 256 315380 0 0 60720 257 309560 365 7 54820 258 298466 364 26 47705 259 315380 0 0 60720 260 315380 0 0 60720 261 315380 0 0 60720 262 315380 0 0 60720 263 343929 226 13 72835 264 331955 307 10 58856 265 315380 0 0 60720 266 315380 0 0 60720 267 315380 0 0 60720 268 381180 188 9 77655 269 315380 0 0 60720 270 331420 138 26 69817 271 315380 0 0 60720 272 315380 0 0 60720 273 315380 0 0 60720 274 310201 125 19 60798 275 315380 0 0 60720 276 320016 282 12 62452 277 320398 335 23 64175 278 315380 0 0 60720 279 291841 1324 29 67440 280 310670 176 8 68136 281 315380 0 0 60720 282 315380 0 0 60720 283 313491 249 26 56726 284 315380 0 0 60720 285 331323 333 9 70811 286 315380 0 0 60720 287 319210 601 5 60720 288 318098 30 3 62045 289 315380 0 0 60720 290 292754 249 13 54323 291 315380 0 0 60720 292 325176 165 12 62841 293 365959 453 19 81125 294 315380 0 0 60720 295 302409 53 10 59506 296 340968 382 9 59365 297 315380 0 0 60720 298 315380 0 0 60720 299 315380 0 0 60720 300 315380 0 9 60720 301 313164 30 4 60798 302 301164 290 1 58790 303 315380 0 1 60720 304 315380 0 0 60720 305 344425 366 14 61808 306 315394 2 12 60735 307 315380 0 0 60720 308 316647 209 19 64016 309 309836 384 17 54683 310 315380 0 0 60720 311 315380 0 0 60720 312 346611 365 32 87192 313 315380 0 0 60720 314 322031 49 14 64107 315 315656 3 8 60761 316 339445 133 4 65990 317 314964 32 0 59988 318 297141 368 20 61167 319 315372 1 5 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 1 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 4 60379 335 315388 1 1 60727 336 315371 314 4 60720 337 296139 844 20 60925 338 315380 0 0 60720 339 313880 26 1 60896 340 317698 125 10 59734 341 295580 304 12 62969 342 315380 0 0 60720 343 315380 0 0 60720 344 315380 0 0 60720 345 308256 621 13 60720 346 315380 0 0 60720 347 303677 119 3 59118 348 315380 0 0 60720 349 315380 0 0 60720 350 319369 1595 10 60720 351 318690 312 3 58598 352 314049 60 7 61124 353 325699 587 10 59595 354 314210 135 1 62065 355 315380 0 0 60720 356 315380 0 0 60720 357 322378 514 15 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 28 61600 365 316386 180 9 59635 366 315380 0 0 60720 367 315380 0 0 60720 368 315380 0 0 60720 369 315380 0 0 60720 370 315553 218 7 60720 371 315380 0 0 60720 372 323361 448 7 59781 373 336639 227 7 76644 374 307424 174 3 64820 375 315380 0 0 60720 376 315380 0 0 60720 377 295370 121 11 56178 378 322340 607 7 60436 379 319864 2212 10 60720 380 315380 0 0 60720 381 315380 0 0 60720 382 317291 530 18 73433 383 280398 571 14 41477 384 315380 0 0 60720 385 317330 78 12 62700 386 238125 2489 29 67804 387 327071 131 3 59661 388 309038 923 6 58620 389 314210 72 3 60398 390 307930 572 8 58580 391 322327 397 10 62710 392 292136 450 6 59325 393 263276 622 8 60950 394 367655 694 6 68060 395 283910 3425 9 83620 396 283587 562 8 58456 397 243650 4917 26 52811 398 438493 1442 239 121173 399 296261 529 7 63870 400 230621 2126 41 21001 401 304252 1061 3 70415 402 333505 776 8 64230 403 296919 611 6 59190 404 278990 1526 21 69351 405 276898 592 7 64270 406 327007 1182 11 70694 407 317046 621 11 68005 408 304555 989 12 58930 409 298096 438 9 58320 410 231861 726 3 69980 411 309422 1303 57 69863 412 286963 7419 21 63255 413 269753 1164 15 57320 414 448243 3310 32 75230 415 165404 1920 11 79420 416 204325 965 2 73490 417 407159 3256 23 35250 418 290476 1135 20 62285 419 275311 1270 24 69206 420 246541 661 1 65920 421 253468 1013 1 69770 422 240897 2844 74 72683 423 -83265 11528 68 -14545 424 -42143 6526 20 55830 425 272713 2264 20 55174 426 215362 5109 82 67038 427 42754 3999 21 51252 428 306275 35624 244 157278 429 253537 9252 32 79510 430 372631 15236 86 77440 431 -7170 18073 69 27284 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Costs Orders Dividends 1.650e+05 2.201e+01 1.009e+03 1.568e+00 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2198537 -21925 41823 55157 2923782 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.650e+05 3.063e+04 5.386 1.19e-07 *** Costs 2.201e+01 2.155e+00 10.217 < 2e-16 *** Orders 1.009e+03 3.497e+02 2.887 0.004092 ** Dividends 1.568e+00 4.115e-01 3.812 0.000158 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 253500 on 427 degrees of freedom Multiple R-squared: 0.6556, Adjusted R-squared: 0.6532 F-statistic: 271 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 2.430225e-32 1.215113e-32 [2,] 1.0000000 8.737333e-53 4.368667e-53 [3,] 1.0000000 2.605460e-53 1.302730e-53 [4,] 1.0000000 3.343472e-75 1.671736e-75 [5,] 1.0000000 3.787983e-78 1.893992e-78 [6,] 1.0000000 8.070447e-84 4.035223e-84 [7,] 1.0000000 3.228478e-105 1.614239e-105 [8,] 1.0000000 2.615742e-107 1.307871e-107 [9,] 1.0000000 8.429364e-116 4.214682e-116 [10,] 1.0000000 2.058424e-123 1.029212e-123 [11,] 1.0000000 1.769716e-133 8.848580e-134 [12,] 1.0000000 1.349267e-139 6.746335e-140 [13,] 1.0000000 1.674554e-144 8.372772e-145 [14,] 1.0000000 9.980505e-153 4.990253e-153 [15,] 1.0000000 4.714535e-161 2.357268e-161 [16,] 1.0000000 2.522695e-164 1.261347e-164 [17,] 1.0000000 7.254357e-165 3.627178e-165 [18,] 1.0000000 4.523802e-171 2.261901e-171 [19,] 1.0000000 1.148182e-170 5.740911e-171 [20,] 1.0000000 7.161049e-173 3.580524e-173 [21,] 1.0000000 6.187980e-176 3.093990e-176 [22,] 1.0000000 3.089039e-175 1.544519e-175 [23,] 1.0000000 4.705676e-189 2.352838e-189 [24,] 1.0000000 4.328027e-189 2.164014e-189 [25,] 1.0000000 6.441233e-195 3.220616e-195 [26,] 1.0000000 2.904189e-196 1.452095e-196 [27,] 1.0000000 1.152662e-205 5.763311e-206 [28,] 1.0000000 4.250207e-208 2.125104e-208 [29,] 1.0000000 2.165470e-207 1.082735e-207 [30,] 1.0000000 1.983385e-208 9.916925e-209 [31,] 1.0000000 4.715109e-214 2.357554e-214 [32,] 1.0000000 6.242833e-220 3.121416e-220 [33,] 1.0000000 4.577117e-219 2.288558e-219 [34,] 1.0000000 1.143453e-219 5.717266e-220 [35,] 1.0000000 3.043080e-221 1.521540e-221 [36,] 1.0000000 1.874232e-220 9.371158e-221 [37,] 1.0000000 8.048933e-221 4.024467e-221 [38,] 1.0000000 9.401952e-220 4.700976e-220 [39,] 1.0000000 2.087604e-219 1.043802e-219 [40,] 1.0000000 1.807373e-218 9.036866e-219 [41,] 1.0000000 4.582867e-219 2.291434e-219 [42,] 1.0000000 2.172893e-218 1.086446e-218 [43,] 1.0000000 8.204378e-218 4.102189e-218 [44,] 1.0000000 4.341007e-218 2.170504e-218 [45,] 1.0000000 1.573934e-218 7.869672e-219 [46,] 1.0000000 4.482458e-218 2.241229e-218 [47,] 1.0000000 2.729742e-217 1.364871e-217 [48,] 1.0000000 8.420524e-217 4.210262e-217 [49,] 1.0000000 5.048895e-220 2.524448e-220 [50,] 1.0000000 1.097779e-219 5.488896e-220 [51,] 1.0000000 3.514781e-233 1.757391e-233 [52,] 1.0000000 1.415033e-233 7.075164e-234 [53,] 1.0000000 1.300345e-232 6.501724e-233 [54,] 1.0000000 3.363462e-232 1.681731e-232 [55,] 1.0000000 6.710222e-234 3.355111e-234 [56,] 1.0000000 1.970304e-235 9.851518e-236 [57,] 1.0000000 5.446042e-235 2.723021e-235 [58,] 1.0000000 3.000963e-234 1.500481e-234 [59,] 1.0000000 1.669619e-233 8.348097e-234 [60,] 1.0000000 8.471872e-234 4.235936e-234 [61,] 1.0000000 2.271781e-233 1.135891e-233 [62,] 1.0000000 2.140684e-232 1.070342e-232 [63,] 1.0000000 1.213507e-231 6.067534e-232 [64,] 1.0000000 6.701009e-231 3.350504e-231 [65,] 1.0000000 4.718952e-230 2.359476e-230 [66,] 1.0000000 3.813748e-231 1.906874e-231 [67,] 1.0000000 1.309411e-230 6.547056e-231 [68,] 1.0000000 1.123061e-229 5.615305e-230 [69,] 1.0000000 1.348557e-228 6.742784e-229 [70,] 1.0000000 3.161448e-228 1.580724e-228 [71,] 1.0000000 2.852140e-227 1.426070e-227 [72,] 1.0000000 3.834064e-227 1.917032e-227 [73,] 1.0000000 4.032626e-226 2.016313e-226 [74,] 1.0000000 1.524307e-225 7.621534e-226 [75,] 1.0000000 9.068393e-225 4.534197e-225 [76,] 1.0000000 1.496334e-224 7.481670e-225 [77,] 1.0000000 1.458437e-224 7.292187e-225 [78,] 1.0000000 8.173524e-224 4.086762e-224 [79,] 1.0000000 7.939794e-223 3.969897e-223 [80,] 1.0000000 2.122205e-222 1.061102e-222 [81,] 1.0000000 1.693093e-222 8.465467e-223 [82,] 1.0000000 1.266378e-221 6.331888e-222 [83,] 1.0000000 2.736905e-221 1.368453e-221 [84,] 1.0000000 2.662629e-220 1.331314e-220 [85,] 1.0000000 2.420510e-219 1.210255e-219 [86,] 1.0000000 1.576943e-218 7.884716e-219 [87,] 1.0000000 5.307071e-220 2.653536e-220 [88,] 1.0000000 5.400431e-219 2.700215e-219 [89,] 1.0000000 1.923374e-219 9.616871e-220 [90,] 1.0000000 1.707259e-219 8.536295e-220 [91,] 1.0000000 1.532807e-218 7.664036e-219 [92,] 1.0000000 1.505760e-217 7.528801e-218 [93,] 1.0000000 1.317677e-216 6.588387e-217 [94,] 1.0000000 1.295529e-215 6.477644e-216 [95,] 1.0000000 1.169988e-214 5.849940e-215 [96,] 1.0000000 9.190388e-216 4.595194e-216 [97,] 1.0000000 8.431388e-216 4.215694e-216 [98,] 1.0000000 1.403601e-216 7.018006e-217 [99,] 1.0000000 6.191384e-216 3.095692e-216 [100,] 1.0000000 2.187477e-215 1.093738e-215 [101,] 1.0000000 1.345981e-217 6.729903e-218 [102,] 1.0000000 5.648113e-219 2.824057e-219 [103,] 1.0000000 4.298096e-218 2.149048e-218 [104,] 1.0000000 8.618840e-218 4.309420e-218 [105,] 1.0000000 6.133438e-218 3.066719e-218 [106,] 1.0000000 3.810103e-218 1.905052e-218 [107,] 1.0000000 2.946051e-217 1.473026e-217 [108,] 1.0000000 2.387174e-216 1.193587e-216 [109,] 1.0000000 1.584190e-216 7.920952e-217 [110,] 1.0000000 6.364207e-216 3.182103e-216 [111,] 1.0000000 4.584798e-215 2.292399e-215 [112,] 1.0000000 3.685365e-215 1.842682e-215 [113,] 1.0000000 3.977775e-216 1.988887e-216 [114,] 1.0000000 4.087402e-215 2.043701e-215 [115,] 1.0000000 1.975149e-214 9.875747e-215 [116,] 1.0000000 1.832532e-213 9.162659e-214 [117,] 1.0000000 1.589165e-212 7.945826e-213 [118,] 1.0000000 2.197853e-212 1.098926e-212 [119,] 1.0000000 6.315757e-212 3.157878e-212 [120,] 1.0000000 2.062314e-212 1.031157e-212 [121,] 1.0000000 1.925751e-211 9.628753e-212 [122,] 1.0000000 1.901750e-210 9.508748e-211 [123,] 1.0000000 1.818254e-209 9.091270e-210 [124,] 1.0000000 1.709494e-208 8.547469e-209 [125,] 1.0000000 1.668385e-207 8.341924e-208 [126,] 1.0000000 1.434972e-207 7.174861e-208 [127,] 1.0000000 1.372971e-206 6.864856e-207 [128,] 1.0000000 1.422863e-205 7.114317e-206 [129,] 1.0000000 1.407273e-204 7.036365e-205 [130,] 1.0000000 1.363319e-203 6.816593e-204 [131,] 1.0000000 1.326767e-202 6.633836e-203 [132,] 1.0000000 1.274203e-201 6.371017e-202 [133,] 1.0000000 1.218434e-200 6.092172e-201 [134,] 1.0000000 1.162635e-199 5.813177e-200 [135,] 1.0000000 9.621111e-199 4.810556e-199 [136,] 1.0000000 9.180154e-198 4.590077e-198 [137,] 1.0000000 9.398313e-197 4.699157e-197 [138,] 1.0000000 9.064287e-196 4.532144e-196 [139,] 1.0000000 8.295840e-195 4.147920e-195 [140,] 1.0000000 7.606577e-194 3.803289e-194 [141,] 1.0000000 7.212280e-193 3.606140e-193 [142,] 1.0000000 6.995395e-192 3.497697e-192 [143,] 1.0000000 6.022484e-191 3.011242e-191 [144,] 1.0000000 5.692807e-190 2.846403e-190 [145,] 1.0000000 5.542452e-189 2.771226e-189 [146,] 1.0000000 5.721135e-188 2.860568e-188 [147,] 1.0000000 5.433775e-187 2.716888e-187 [148,] 1.0000000 5.331204e-186 2.665602e-186 [149,] 1.0000000 5.012125e-185 2.506062e-185 [150,] 1.0000000 4.663512e-184 2.331756e-184 [151,] 1.0000000 4.330508e-183 2.165254e-183 [152,] 1.0000000 4.066533e-182 2.033267e-182 [153,] 1.0000000 2.844088e-181 1.422044e-181 [154,] 1.0000000 2.625283e-180 1.312641e-180 [155,] 1.0000000 2.417622e-179 1.208811e-179 [156,] 1.0000000 2.220982e-178 1.110491e-178 [157,] 1.0000000 2.035229e-177 1.017615e-177 [158,] 1.0000000 1.943694e-176 9.718470e-177 [159,] 1.0000000 1.813072e-175 9.065358e-176 [160,] 1.0000000 1.653068e-174 8.265341e-175 [161,] 1.0000000 1.496385e-173 7.481924e-174 [162,] 1.0000000 1.854320e-174 9.271601e-175 [163,] 1.0000000 3.310271e-174 1.655135e-174 [164,] 1.0000000 3.032966e-173 1.516483e-173 [165,] 1.0000000 2.770809e-172 1.385405e-172 [166,] 1.0000000 2.578881e-171 1.289440e-171 [167,] 1.0000000 2.341563e-170 1.170782e-170 [168,] 1.0000000 2.197854e-169 1.098927e-169 [169,] 1.0000000 1.982565e-168 9.912826e-169 [170,] 1.0000000 1.808351e-167 9.041757e-168 [171,] 1.0000000 1.698804e-166 8.494022e-167 [172,] 1.0000000 1.507778e-165 7.538891e-166 [173,] 1.0000000 1.340523e-164 6.702614e-165 [174,] 1.0000000 1.187954e-163 5.939772e-164 [175,] 1.0000000 7.923520e-163 3.961760e-163 [176,] 1.0000000 7.261153e-162 3.630577e-162 [177,] 1.0000000 4.232342e-161 2.116171e-161 [178,] 1.0000000 2.572383e-160 1.286192e-160 [179,] 1.0000000 2.296461e-159 1.148231e-159 [180,] 1.0000000 1.533100e-158 7.665501e-159 [181,] 1.0000000 5.734398e-158 2.867199e-158 [182,] 1.0000000 5.100126e-157 2.550063e-157 [183,] 1.0000000 3.605273e-156 1.802636e-156 [184,] 1.0000000 3.105905e-155 1.552952e-155 [185,] 1.0000000 2.543144e-154 1.271572e-154 [186,] 1.0000000 2.175442e-153 1.087721e-153 [187,] 1.0000000 1.885944e-152 9.429722e-153 [188,] 1.0000000 1.601479e-151 8.007393e-152 [189,] 1.0000000 1.360697e-150 6.803483e-151 [190,] 1.0000000 1.147050e-149 5.735252e-150 [191,] 1.0000000 8.798543e-149 4.399272e-149 [192,] 1.0000000 7.360103e-148 3.680051e-148 [193,] 1.0000000 6.134258e-147 3.067129e-147 [194,] 1.0000000 6.352474e-147 3.176237e-147 [195,] 1.0000000 5.307375e-146 2.653687e-146 [196,] 1.0000000 4.417763e-145 2.208881e-145 [197,] 1.0000000 3.663540e-144 1.831770e-144 [198,] 1.0000000 2.811964e-143 1.405982e-143 [199,] 1.0000000 1.236683e-142 6.183415e-143 [200,] 1.0000000 1.028928e-141 5.144639e-142 [201,] 1.0000000 8.559826e-141 4.279913e-141 [202,] 1.0000000 6.988784e-140 3.494392e-140 [203,] 1.0000000 1.006209e-139 5.031044e-140 [204,] 1.0000000 8.412725e-139 4.206363e-139 [205,] 1.0000000 6.838441e-138 3.419221e-138 [206,] 1.0000000 5.537162e-137 2.768581e-137 [207,] 1.0000000 3.356023e-136 1.678011e-136 [208,] 1.0000000 2.523185e-135 1.261592e-135 [209,] 1.0000000 2.072258e-134 1.036129e-134 [210,] 1.0000000 1.656891e-133 8.284456e-134 [211,] 1.0000000 1.323388e-132 6.616938e-133 [212,] 1.0000000 1.022267e-131 5.111336e-132 [213,] 1.0000000 8.074498e-131 4.037249e-131 [214,] 1.0000000 3.870213e-130 1.935107e-130 [215,] 1.0000000 3.048905e-129 1.524452e-129 [216,] 1.0000000 2.379499e-128 1.189749e-128 [217,] 1.0000000 1.850683e-127 9.253415e-128 [218,] 1.0000000 1.452543e-126 7.262714e-127 [219,] 1.0000000 1.120077e-125 5.600387e-126 [220,] 1.0000000 8.730633e-125 4.365317e-125 [221,] 1.0000000 6.729103e-124 3.364551e-124 [222,] 1.0000000 5.121344e-123 2.560672e-123 [223,] 1.0000000 3.921146e-122 1.960573e-122 [224,] 1.0000000 2.960148e-121 1.480074e-121 [225,] 1.0000000 2.148883e-120 1.074442e-120 [226,] 1.0000000 1.482857e-119 7.414286e-120 [227,] 1.0000000 1.106131e-118 5.530656e-119 [228,] 1.0000000 8.426124e-118 4.213062e-118 [229,] 1.0000000 6.232351e-117 3.116175e-117 [230,] 1.0000000 4.589173e-116 2.294586e-116 [231,] 1.0000000 3.365783e-115 1.682891e-115 [232,] 1.0000000 2.447622e-114 1.223811e-114 [233,] 1.0000000 1.779829e-113 8.899144e-114 [234,] 1.0000000 1.288705e-112 6.443523e-113 [235,] 1.0000000 9.291018e-112 4.645509e-112 [236,] 1.0000000 6.764258e-111 3.382129e-111 [237,] 1.0000000 4.834640e-110 2.417320e-110 [238,] 1.0000000 3.407681e-109 1.703841e-109 [239,] 1.0000000 2.473491e-108 1.236746e-108 [240,] 1.0000000 1.744328e-107 8.721641e-108 [241,] 1.0000000 1.224741e-106 6.123706e-107 [242,] 1.0000000 8.784312e-106 4.392156e-106 [243,] 1.0000000 6.113228e-105 3.056614e-105 [244,] 1.0000000 4.235583e-104 2.117791e-104 [245,] 1.0000000 2.901977e-103 1.450988e-103 [246,] 1.0000000 1.992722e-102 9.963610e-103 [247,] 1.0000000 1.368284e-101 6.841422e-102 [248,] 1.0000000 7.701955e-101 3.850978e-101 [249,] 1.0000000 6.137527e-101 3.068764e-101 [250,] 1.0000000 4.201373e-100 2.100687e-100 [251,] 1.0000000 2.627559e-99 1.313780e-99 [252,] 1.0000000 1.495040e-98 7.475198e-99 [253,] 1.0000000 1.010580e-97 5.052898e-98 [254,] 1.0000000 6.800028e-97 3.400014e-97 [255,] 1.0000000 4.554762e-96 2.277381e-96 [256,] 1.0000000 3.036884e-95 1.518442e-95 [257,] 1.0000000 2.041187e-94 1.020594e-94 [258,] 1.0000000 1.071082e-93 5.355410e-94 [259,] 1.0000000 7.049114e-93 3.524557e-93 [260,] 1.0000000 4.617682e-92 2.308841e-92 [261,] 1.0000000 3.010811e-91 1.505406e-91 [262,] 1.0000000 1.605316e-90 8.026580e-91 [263,] 1.0000000 1.037604e-89 5.188020e-90 [264,] 1.0000000 6.849416e-89 3.424708e-89 [265,] 1.0000000 4.385247e-88 2.192624e-88 [266,] 1.0000000 2.794251e-87 1.397125e-87 [267,] 1.0000000 1.771981e-86 8.859904e-87 [268,] 1.0000000 1.143321e-85 5.716604e-86 [269,] 1.0000000 7.180531e-85 3.590266e-85 [270,] 1.0000000 4.472021e-84 2.236011e-84 [271,] 1.0000000 2.825713e-83 1.412856e-83 [272,] 1.0000000 1.748768e-82 8.743841e-83 [273,] 1.0000000 9.991851e-82 4.995925e-82 [274,] 1.0000000 5.989676e-81 2.994838e-81 [275,] 1.0000000 3.659053e-80 1.829526e-80 [276,] 1.0000000 2.224255e-79 1.112127e-79 [277,] 1.0000000 1.274548e-78 6.372738e-79 [278,] 1.0000000 7.672426e-78 3.836213e-78 [279,] 1.0000000 4.695320e-77 2.347660e-77 [280,] 1.0000000 2.798374e-76 1.399187e-76 [281,] 1.0000000 1.583245e-75 7.916224e-76 [282,] 1.0000000 9.375018e-75 4.687509e-75 [283,] 1.0000000 5.503239e-74 2.751619e-74 [284,] 1.0000000 3.271996e-73 1.635998e-73 [285,] 1.0000000 1.900963e-72 9.504813e-73 [286,] 1.0000000 1.075367e-71 5.376835e-72 [287,] 1.0000000 6.232573e-71 3.116286e-71 [288,] 1.0000000 3.564920e-70 1.782460e-70 [289,] 1.0000000 2.076835e-69 1.038418e-69 [290,] 1.0000000 8.469871e-69 4.234935e-69 [291,] 1.0000000 4.776692e-68 2.388346e-68 [292,] 1.0000000 2.679562e-67 1.339781e-67 [293,] 1.0000000 1.495114e-66 7.475570e-67 [294,] 1.0000000 8.328321e-66 4.164161e-66 [295,] 1.0000000 4.641899e-65 2.320950e-65 [296,] 1.0000000 2.602864e-64 1.301432e-64 [297,] 1.0000000 1.421785e-63 7.108925e-64 [298,] 1.0000000 7.720436e-63 3.860218e-63 [299,] 1.0000000 3.201798e-62 1.600899e-62 [300,] 1.0000000 1.729495e-61 8.647474e-62 [301,] 1.0000000 9.248263e-61 4.624132e-61 [302,] 1.0000000 5.033930e-60 2.516965e-60 [303,] 1.0000000 2.447381e-59 1.223691e-59 [304,] 1.0000000 1.287200e-58 6.435999e-59 [305,] 1.0000000 6.731408e-58 3.365704e-58 [306,] 1.0000000 3.172027e-57 1.586013e-57 [307,] 1.0000000 1.643278e-56 8.216392e-57 [308,] 1.0000000 8.566117e-56 4.283058e-56 [309,] 1.0000000 4.395012e-55 2.197506e-55 [310,] 1.0000000 2.102942e-54 1.051471e-54 [311,] 1.0000000 1.055943e-53 5.279714e-54 [312,] 1.0000000 5.451191e-53 2.725596e-53 [313,] 1.0000000 2.731798e-52 1.365899e-52 [314,] 1.0000000 1.358346e-51 6.791732e-52 [315,] 1.0000000 6.713189e-51 3.356595e-51 [316,] 1.0000000 3.297504e-50 1.648752e-50 [317,] 1.0000000 1.609760e-49 8.048801e-50 [318,] 1.0000000 7.809751e-49 3.904875e-49 [319,] 1.0000000 3.765249e-48 1.882624e-48 [320,] 1.0000000 1.822590e-47 9.112948e-48 [321,] 1.0000000 8.676706e-47 4.338353e-47 [322,] 1.0000000 4.104312e-46 2.052156e-46 [323,] 1.0000000 1.928963e-45 9.644817e-46 [324,] 1.0000000 9.007075e-45 4.503538e-45 [325,] 1.0000000 4.178279e-44 2.089140e-44 [326,] 1.0000000 1.925494e-43 9.627469e-44 [327,] 1.0000000 8.814431e-43 4.407215e-43 [328,] 1.0000000 4.014123e-42 2.007061e-42 [329,] 1.0000000 1.813715e-41 9.068577e-42 [330,] 1.0000000 8.059688e-41 4.029844e-41 [331,] 1.0000000 3.711068e-40 1.855534e-40 [332,] 1.0000000 1.642302e-39 8.211508e-40 [333,] 1.0000000 7.269284e-39 3.634642e-39 [334,] 1.0000000 3.095741e-38 1.547870e-38 [335,] 1.0000000 1.374910e-37 6.874549e-38 [336,] 1.0000000 5.920336e-37 2.960168e-37 [337,] 1.0000000 2.530858e-36 1.265429e-36 [338,] 1.0000000 1.074009e-35 5.370047e-36 [339,] 1.0000000 4.578785e-35 2.289392e-35 [340,] 1.0000000 1.914534e-34 9.572672e-35 [341,] 1.0000000 8.122871e-34 4.061436e-34 [342,] 1.0000000 3.345712e-33 1.672856e-33 [343,] 1.0000000 1.367492e-32 6.837462e-33 [344,] 1.0000000 5.012340e-32 2.506170e-32 [345,] 1.0000000 1.897526e-31 9.487632e-32 [346,] 1.0000000 7.659029e-31 3.829515e-31 [347,] 1.0000000 2.728585e-30 1.364293e-30 [348,] 1.0000000 1.088799e-29 5.443995e-30 [349,] 1.0000000 4.246838e-29 2.123419e-29 [350,] 1.0000000 1.642728e-28 8.213639e-29 [351,] 1.0000000 6.318902e-28 3.159451e-28 [352,] 1.0000000 2.409522e-27 1.204761e-27 [353,] 1.0000000 9.109048e-27 4.554524e-27 [354,] 1.0000000 3.413653e-26 1.706826e-26 [355,] 1.0000000 1.270354e-25 6.351768e-26 [356,] 1.0000000 4.676789e-25 2.338395e-25 [357,] 1.0000000 1.706157e-24 8.530787e-25 [358,] 1.0000000 6.082174e-24 3.041087e-24 [359,] 1.0000000 2.135228e-23 1.067614e-23 [360,] 1.0000000 7.578590e-23 3.789295e-23 [361,] 1.0000000 2.664120e-22 1.332060e-22 [362,] 1.0000000 9.274171e-22 4.637086e-22 [363,] 1.0000000 3.196572e-21 1.598286e-21 [364,] 1.0000000 1.085744e-20 5.428722e-21 [365,] 1.0000000 3.667283e-20 1.833642e-20 [366,] 1.0000000 1.139737e-19 5.698684e-20 [367,] 1.0000000 3.983713e-19 1.991857e-19 [368,] 1.0000000 1.374268e-18 6.871339e-19 [369,] 1.0000000 4.456894e-18 2.228447e-18 [370,] 1.0000000 1.429328e-17 7.146640e-18 [371,] 1.0000000 4.676969e-17 2.338484e-17 [372,] 1.0000000 1.379875e-16 6.899375e-17 [373,] 1.0000000 3.754796e-16 1.877398e-16 [374,] 1.0000000 1.150579e-15 5.752895e-16 [375,] 1.0000000 3.482185e-15 1.741092e-15 [376,] 1.0000000 1.120536e-14 5.602678e-15 [377,] 1.0000000 2.962742e-14 1.481371e-14 [378,] 1.0000000 8.659144e-14 4.329572e-14 [379,] 1.0000000 2.549300e-13 1.274650e-13 [380,] 1.0000000 6.813281e-13 3.406640e-13 [381,] 1.0000000 1.781602e-12 8.908010e-13 [382,] 1.0000000 4.844327e-12 2.422163e-12 [383,] 1.0000000 1.336749e-11 6.683747e-12 [384,] 1.0000000 3.605059e-11 1.802530e-11 [385,] 1.0000000 9.509091e-11 4.754545e-11 [386,] 1.0000000 2.668014e-10 1.334007e-10 [387,] 1.0000000 7.603558e-10 3.801779e-10 [388,] 1.0000000 1.526761e-09 7.633805e-10 [389,] 1.0000000 4.142377e-09 2.071188e-09 [390,] 1.0000000 1.119048e-08 5.595240e-09 [391,] 1.0000000 2.889729e-08 1.444865e-08 [392,] 1.0000000 5.852557e-08 2.926279e-08 [393,] 0.9999999 1.525159e-07 7.625794e-08 [394,] 0.9999998 3.216941e-07 1.608470e-07 [395,] 0.9999996 8.187020e-07 4.093510e-07 [396,] 0.9999991 1.759423e-06 8.797116e-07 [397,] 0.9999980 4.097108e-06 2.048554e-06 [398,] 0.9999949 1.011585e-05 5.057927e-06 [399,] 0.9999878 2.439913e-05 1.219956e-05 [400,] 0.9999732 5.355787e-05 2.677894e-05 [401,] 0.9999412 1.176265e-04 5.881327e-05 [402,] 0.9998805 2.390873e-04 1.195437e-04 [403,] 0.9997562 4.875105e-04 2.437553e-04 [404,] 0.9994845 1.030921e-03 5.154604e-04 [405,] 0.9989024 2.195244e-03 1.097622e-03 [406,] 0.9982032 3.593599e-03 1.796799e-03 [407,] 0.9964341 7.131734e-03 3.565867e-03 [408,] 0.9968652 6.269549e-03 3.134775e-03 [409,] 0.9949695 1.006091e-02 5.030456e-03 [410,] 0.9905243 1.895148e-02 9.475738e-03 [411,] 0.9977197 4.560506e-03 2.280253e-03 [412,] 0.9952420 9.516004e-03 4.758002e-03 [413,] 0.9890590 2.188210e-02 1.094105e-02 [414,] 0.9754770 4.904606e-02 2.452303e-02 [415,] 0.9482403 1.035194e-01 5.175969e-02 [416,] 0.8942298 2.115405e-01 1.057702e-01 [417,] 0.8019672 3.960655e-01 1.980328e-01 [418,] 0.8245660 3.508681e-01 1.754340e-01 > postscript(file="/var/www/html/rcomp/tmp/1ulcq1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2mvcb1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3mvcb1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4mvcb1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5f4bw1292936583.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 1389642.9555 2923781.5148 1363354.5309 -2198536.5879 -394043.7941 6 7 8 9 10 109880.6241 -135942.2407 128929.0956 -336335.1695 715090.7022 11 12 13 14 15 -150648.2836 342304.4900 449640.5419 -1139733.8466 -214270.3673 16 17 18 19 20 460576.1545 -237908.5126 -238086.3678 -295879.7711 -585121.5240 21 22 23 24 25 203466.5526 -137406.0789 -119601.0233 -638071.4974 -318325.0883 26 27 28 29 30 -84290.3138 -471904.3808 -340382.1733 348304.7876 -261063.8100 31 32 33 34 35 286671.4176 -87934.3094 255923.1787 130741.0876 -54935.4616 36 37 38 39 40 -569154.1520 34312.4748 157831.6796 -51019.2550 -304575.1497 41 42 43 44 45 -459635.3342 -329005.7740 168.5371 -209377.2166 -27236.7467 46 47 48 49 50 -53785.8121 -281550.6842 -264265.0929 -53546.3755 -181996.3972 51 52 53 54 55 35728.4327 -86627.1882 -228213.2479 44143.8787 -394464.4083 56 57 58 59 60 14554.8677 320398.0866 -22045.0506 -59459.7426 -205533.9695 61 62 63 64 65 14512.1235 -179051.4887 7024.0939 -174142.4494 -15772.7669 66 67 68 69 70 -494815.7464 -65001.8342 -14962.1865 -100729.1806 -72008.4810 71 72 73 74 75 -122294.9668 1254.8171 -88695.3911 15115.6276 -34735.1194 76 77 78 79 80 -98719.3209 -28468.6427 -73250.3450 -42230.1873 -37806.4900 81 82 83 84 85 27097.0836 -49658.1948 -82701.2953 41682.6321 24953.4226 86 87 88 89 90 -102380.0855 -29335.5087 1174.8496 -115197.8123 32519.6180 91 92 93 94 95 -1901.5977 -21805.7745 69499.4547 -172296.0652 49811.4065 96 97 98 99 100 19621.4688 71825.6583 2142.7165 -41169.6913 -1469.7132 101 102 103 104 105 -35892.6448 -112495.2638 -165507.7762 -50097.8850 30071.9758 106 107 108 109 110 -53262.0014 -252162.6722 131140.5689 -175230.3204 -92234.7579 111 112 113 114 115 94522.0028 -190791.2910 -165266.5841 -57462.7139 86193.8813 116 117 118 119 120 21656.7901 17528.3553 72494.7632 -132141.0857 14990.3877 121 122 123 124 125 -288568.2348 32830.8847 -103092.5235 84908.0487 13950.5787 126 127 128 129 130 -162466.3229 55156.8048 42031.7011 29693.9470 55156.8048 131 132 133 134 135 46456.5140 57662.4310 54147.3844 31225.8242 33614.4964 136 137 138 139 140 44128.6218 45700.9845 53544.3097 30904.7354 55156.8048 141 142 143 144 145 51357.0400 55156.8048 11408.9397 50507.6768 40173.9040 146 147 148 149 150 43331.3385 55156.8048 47319.3565 57040.9193 54881.3370 151 152 153 154 155 16159.7755 -39133.5926 50109.7026 23678.9119 52220.7451 156 157 158 159 160 55156.8048 55156.8048 49628.8010 62042.6027 55156.8048 161 162 163 164 165 55156.8048 55156.8048 55156.8048 28173.1243 3757.5271 166 167 168 169 170 52614.6829 55156.8048 111229.1895 -76650.7953 55156.8048 171 172 173 174 175 55156.8048 47044.2411 55156.8048 33908.6287 55156.8048 176 177 178 179 180 44112.2506 -69826.6760 54350.5341 55156.8048 55156.8048 181 182 183 184 185 41403.7863 29922.3002 -17619.1547 -4775.3971 46151.2683 186 187 188 189 190 48386.6125 11520.0678 41700.1465 13659.9087 55156.8048 191 192 193 194 195 47315.2565 55156.8048 43545.4751 55156.8048 51553.4904 196 197 198 199 200 55156.8048 -132869.2446 55156.8048 55156.8048 -50516.3690 201 202 203 204 205 55156.8048 55156.8048 55156.8048 39709.4000 393.0099 206 207 208 209 210 49196.4916 47153.0040 55156.8048 51267.7733 37862.5884 211 212 213 214 215 55156.8048 55156.8048 57684.9961 6260.9944 24478.6930 216 217 218 219 220 53137.9639 48903.9612 41823.0012 53137.9639 -185.6711 221 222 223 224 225 27393.6179 55156.8048 54089.8453 21446.4109 53769.7683 226 227 228 229 230 25298.5915 43900.1103 55156.8048 38724.2768 55156.8048 231 232 233 234 235 19420.9584 35821.0062 55156.8048 20267.7317 55156.8048 236 237 238 239 240 55183.7008 55156.8048 48928.3830 55156.8048 55156.8048 241 242 243 244 245 55156.8048 46562.3377 55156.8048 49234.2381 -7954.4067 246 247 248 249 250 55156.8048 55156.8048 29034.0780 55156.8048 55156.8048 251 252 253 254 255 34794.7895 55156.8048 36884.1340 -27169.9018 -149947.4502 256 257 258 259 260 55156.8048 43489.2216 24397.7124 55156.8048 55156.8048 261 262 263 264 265 55156.8048 55156.8048 46606.2917 57802.5942 55156.8048 266 267 268 269 270 55156.8048 55156.8048 81171.6505 55156.8048 27645.6970 271 272 273 274 275 55156.8048 55156.8048 55156.8048 27924.6122 55156.8048 276 277 278 279 280 38755.0091 24164.1677 55156.8048 -37343.0787 26865.2489 281 282 283 284 285 55156.8048 55156.8048 27805.5139 55156.8048 38856.9030 286 287 288 289 290 55156.8048 40708.7327 52107.9101 55156.8048 23959.9458 291 292 293 294 295 55156.8048 45880.6314 44579.9973 55156.8048 32828.8978 296 297 298 299 300 65375.5583 55156.8048 55156.8048 55156.8048 46072.0207 301 302 303 304 305 48120.3368 36574.1491 54147.3844 55156.8048 60305.9911 306 307 308 309 310 42990.2029 55156.8048 27474.1123 33467.6100 55156.8048 311 312 313 314 315 55156.8048 4531.1075 55156.8048 41284.8743 47227.0903 316 317 318 319 320 63990.4489 55184.4267 7926.8099 50081.2561 55156.8048 321 322 323 324 325 55156.8048 55156.8048 55156.8048 55156.8048 55156.8048 326 327 328 329 330 50781.9191 55156.8048 55156.8048 55156.8048 55156.8048 331 332 333 334 335 55156.8048 55156.8048 55156.8048 47889.5289 54122.3903 336 337 338 339 340 44197.4366 -3174.7314 55156.8048 51798.9505 46175.2185 341 342 343 344 345 13023.7963 55156.8048 55156.8048 55156.8048 21239.0706 346 347 348 349 350 55156.8048 40318.4114 55156.8048 55156.8048 13937.7951 351 352 353 354 355 51898.1207 44805.3156 44223.3366 47895.8151 55156.8048 356 357 358 359 360 55156.8048 7371.7650 55156.8048 55156.8048 55156.8048 361 362 363 364 365 51111.9712 55156.8048 55156.8048 -20343.5093 44817.0938 366 367 368 369 370 55156.8048 55156.8048 55156.8048 55156.8048 43464.6080 371 372 373 374 375 55156.8048 47681.9423 39376.6041 33911.3345 55156.8048 376 377 378 379 380 55156.8048 28503.2373 42133.2398 849.5865 55156.8048 381 382 383 384 385 55156.8048 7290.7240 23653.9248 55156.8048 40171.0802 386 387 388 389 390 -117277.3791 62596.5683 25732.2345 49878.5072 30395.3699 391 392 393 394 395 40148.4756 18137.5424 -19076.5837 74584.5603 -96716.4176 396 397 398 399 400 6467.0060 -138660.7052 -189544.2257 12385.3742 -55495.2182 401 402 403 404 405 2436.6634 42617.6282 19587.8792 -49563.0225 -8991.9428 406 407 408 409 410 14014.8986 20621.8119 13253.5098 22909.7443 -61897.0645 411 412 413 414 415 -51363.8727 -161762.7527 -25904.1719 60090.8563 -177521.3404 416 417 418 419 420 -99190.4423 91986.7555 -17377.1556 -50407.0390 -37399.3791 421 422 423 424 425 -44260.1342 -175397.0259 -547868.0565 -458554.3229 -48841.8151 426 427 428 429 430 -250017.3396 -311854.7009 -1135952.0870 -272140.7589 -336046.0989 431 -682476.5187 > postscript(file="/var/www/html/rcomp/tmp/6f4bw1292936583.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 1389642.9555 NA 1 2923781.5148 1389642.9555 2 1363354.5309 2923781.5148 3 -2198536.5879 1363354.5309 4 -394043.7941 -2198536.5879 5 109880.6241 -394043.7941 6 -135942.2407 109880.6241 7 128929.0956 -135942.2407 8 -336335.1695 128929.0956 9 715090.7022 -336335.1695 10 -150648.2836 715090.7022 11 342304.4900 -150648.2836 12 449640.5419 342304.4900 13 -1139733.8466 449640.5419 14 -214270.3673 -1139733.8466 15 460576.1545 -214270.3673 16 -237908.5126 460576.1545 17 -238086.3678 -237908.5126 18 -295879.7711 -238086.3678 19 -585121.5240 -295879.7711 20 203466.5526 -585121.5240 21 -137406.0789 203466.5526 22 -119601.0233 -137406.0789 23 -638071.4974 -119601.0233 24 -318325.0883 -638071.4974 25 -84290.3138 -318325.0883 26 -471904.3808 -84290.3138 27 -340382.1733 -471904.3808 28 348304.7876 -340382.1733 29 -261063.8100 348304.7876 30 286671.4176 -261063.8100 31 -87934.3094 286671.4176 32 255923.1787 -87934.3094 33 130741.0876 255923.1787 34 -54935.4616 130741.0876 35 -569154.1520 -54935.4616 36 34312.4748 -569154.1520 37 157831.6796 34312.4748 38 -51019.2550 157831.6796 39 -304575.1497 -51019.2550 40 -459635.3342 -304575.1497 41 -329005.7740 -459635.3342 42 168.5371 -329005.7740 43 -209377.2166 168.5371 44 -27236.7467 -209377.2166 45 -53785.8121 -27236.7467 46 -281550.6842 -53785.8121 47 -264265.0929 -281550.6842 48 -53546.3755 -264265.0929 49 -181996.3972 -53546.3755 50 35728.4327 -181996.3972 51 -86627.1882 35728.4327 52 -228213.2479 -86627.1882 53 44143.8787 -228213.2479 54 -394464.4083 44143.8787 55 14554.8677 -394464.4083 56 320398.0866 14554.8677 57 -22045.0506 320398.0866 58 -59459.7426 -22045.0506 59 -205533.9695 -59459.7426 60 14512.1235 -205533.9695 61 -179051.4887 14512.1235 62 7024.0939 -179051.4887 63 -174142.4494 7024.0939 64 -15772.7669 -174142.4494 65 -494815.7464 -15772.7669 66 -65001.8342 -494815.7464 67 -14962.1865 -65001.8342 68 -100729.1806 -14962.1865 69 -72008.4810 -100729.1806 70 -122294.9668 -72008.4810 71 1254.8171 -122294.9668 72 -88695.3911 1254.8171 73 15115.6276 -88695.3911 74 -34735.1194 15115.6276 75 -98719.3209 -34735.1194 76 -28468.6427 -98719.3209 77 -73250.3450 -28468.6427 78 -42230.1873 -73250.3450 79 -37806.4900 -42230.1873 80 27097.0836 -37806.4900 81 -49658.1948 27097.0836 82 -82701.2953 -49658.1948 83 41682.6321 -82701.2953 84 24953.4226 41682.6321 85 -102380.0855 24953.4226 86 -29335.5087 -102380.0855 87 1174.8496 -29335.5087 88 -115197.8123 1174.8496 89 32519.6180 -115197.8123 90 -1901.5977 32519.6180 91 -21805.7745 -1901.5977 92 69499.4547 -21805.7745 93 -172296.0652 69499.4547 94 49811.4065 -172296.0652 95 19621.4688 49811.4065 96 71825.6583 19621.4688 97 2142.7165 71825.6583 98 -41169.6913 2142.7165 99 -1469.7132 -41169.6913 100 -35892.6448 -1469.7132 101 -112495.2638 -35892.6448 102 -165507.7762 -112495.2638 103 -50097.8850 -165507.7762 104 30071.9758 -50097.8850 105 -53262.0014 30071.9758 106 -252162.6722 -53262.0014 107 131140.5689 -252162.6722 108 -175230.3204 131140.5689 109 -92234.7579 -175230.3204 110 94522.0028 -92234.7579 111 -190791.2910 94522.0028 112 -165266.5841 -190791.2910 113 -57462.7139 -165266.5841 114 86193.8813 -57462.7139 115 21656.7901 86193.8813 116 17528.3553 21656.7901 117 72494.7632 17528.3553 118 -132141.0857 72494.7632 119 14990.3877 -132141.0857 120 -288568.2348 14990.3877 121 32830.8847 -288568.2348 122 -103092.5235 32830.8847 123 84908.0487 -103092.5235 124 13950.5787 84908.0487 125 -162466.3229 13950.5787 126 55156.8048 -162466.3229 127 42031.7011 55156.8048 128 29693.9470 42031.7011 129 55156.8048 29693.9470 130 46456.5140 55156.8048 131 57662.4310 46456.5140 132 54147.3844 57662.4310 133 31225.8242 54147.3844 134 33614.4964 31225.8242 135 44128.6218 33614.4964 136 45700.9845 44128.6218 137 53544.3097 45700.9845 138 30904.7354 53544.3097 139 55156.8048 30904.7354 140 51357.0400 55156.8048 141 55156.8048 51357.0400 142 11408.9397 55156.8048 143 50507.6768 11408.9397 144 40173.9040 50507.6768 145 43331.3385 40173.9040 146 55156.8048 43331.3385 147 47319.3565 55156.8048 148 57040.9193 47319.3565 149 54881.3370 57040.9193 150 16159.7755 54881.3370 151 -39133.5926 16159.7755 152 50109.7026 -39133.5926 153 23678.9119 50109.7026 154 52220.7451 23678.9119 155 55156.8048 52220.7451 156 55156.8048 55156.8048 157 49628.8010 55156.8048 158 62042.6027 49628.8010 159 55156.8048 62042.6027 160 55156.8048 55156.8048 161 55156.8048 55156.8048 162 55156.8048 55156.8048 163 28173.1243 55156.8048 164 3757.5271 28173.1243 165 52614.6829 3757.5271 166 55156.8048 52614.6829 167 111229.1895 55156.8048 168 -76650.7953 111229.1895 169 55156.8048 -76650.7953 170 55156.8048 55156.8048 171 47044.2411 55156.8048 172 55156.8048 47044.2411 173 33908.6287 55156.8048 174 55156.8048 33908.6287 175 44112.2506 55156.8048 176 -69826.6760 44112.2506 177 54350.5341 -69826.6760 178 55156.8048 54350.5341 179 55156.8048 55156.8048 180 41403.7863 55156.8048 181 29922.3002 41403.7863 182 -17619.1547 29922.3002 183 -4775.3971 -17619.1547 184 46151.2683 -4775.3971 185 48386.6125 46151.2683 186 11520.0678 48386.6125 187 41700.1465 11520.0678 188 13659.9087 41700.1465 189 55156.8048 13659.9087 190 47315.2565 55156.8048 191 55156.8048 47315.2565 192 43545.4751 55156.8048 193 55156.8048 43545.4751 194 51553.4904 55156.8048 195 55156.8048 51553.4904 196 -132869.2446 55156.8048 197 55156.8048 -132869.2446 198 55156.8048 55156.8048 199 -50516.3690 55156.8048 200 55156.8048 -50516.3690 201 55156.8048 55156.8048 202 55156.8048 55156.8048 203 39709.4000 55156.8048 204 393.0099 39709.4000 205 49196.4916 393.0099 206 47153.0040 49196.4916 207 55156.8048 47153.0040 208 51267.7733 55156.8048 209 37862.5884 51267.7733 210 55156.8048 37862.5884 211 55156.8048 55156.8048 212 57684.9961 55156.8048 213 6260.9944 57684.9961 214 24478.6930 6260.9944 215 53137.9639 24478.6930 216 48903.9612 53137.9639 217 41823.0012 48903.9612 218 53137.9639 41823.0012 219 -185.6711 53137.9639 220 27393.6179 -185.6711 221 55156.8048 27393.6179 222 54089.8453 55156.8048 223 21446.4109 54089.8453 224 53769.7683 21446.4109 225 25298.5915 53769.7683 226 43900.1103 25298.5915 227 55156.8048 43900.1103 228 38724.2768 55156.8048 229 55156.8048 38724.2768 230 19420.9584 55156.8048 231 35821.0062 19420.9584 232 55156.8048 35821.0062 233 20267.7317 55156.8048 234 55156.8048 20267.7317 235 55183.7008 55156.8048 236 55156.8048 55183.7008 237 48928.3830 55156.8048 238 55156.8048 48928.3830 239 55156.8048 55156.8048 240 55156.8048 55156.8048 241 46562.3377 55156.8048 242 55156.8048 46562.3377 243 49234.2381 55156.8048 244 -7954.4067 49234.2381 245 55156.8048 -7954.4067 246 55156.8048 55156.8048 247 29034.0780 55156.8048 248 55156.8048 29034.0780 249 55156.8048 55156.8048 250 34794.7895 55156.8048 251 55156.8048 34794.7895 252 36884.1340 55156.8048 253 -27169.9018 36884.1340 254 -149947.4502 -27169.9018 255 55156.8048 -149947.4502 256 43489.2216 55156.8048 257 24397.7124 43489.2216 258 55156.8048 24397.7124 259 55156.8048 55156.8048 260 55156.8048 55156.8048 261 55156.8048 55156.8048 262 46606.2917 55156.8048 263 57802.5942 46606.2917 264 55156.8048 57802.5942 265 55156.8048 55156.8048 266 55156.8048 55156.8048 267 81171.6505 55156.8048 268 55156.8048 81171.6505 269 27645.6970 55156.8048 270 55156.8048 27645.6970 271 55156.8048 55156.8048 272 55156.8048 55156.8048 273 27924.6122 55156.8048 274 55156.8048 27924.6122 275 38755.0091 55156.8048 276 24164.1677 38755.0091 277 55156.8048 24164.1677 278 -37343.0787 55156.8048 279 26865.2489 -37343.0787 280 55156.8048 26865.2489 281 55156.8048 55156.8048 282 27805.5139 55156.8048 283 55156.8048 27805.5139 284 38856.9030 55156.8048 285 55156.8048 38856.9030 286 40708.7327 55156.8048 287 52107.9101 40708.7327 288 55156.8048 52107.9101 289 23959.9458 55156.8048 290 55156.8048 23959.9458 291 45880.6314 55156.8048 292 44579.9973 45880.6314 293 55156.8048 44579.9973 294 32828.8978 55156.8048 295 65375.5583 32828.8978 296 55156.8048 65375.5583 297 55156.8048 55156.8048 298 55156.8048 55156.8048 299 46072.0207 55156.8048 300 48120.3368 46072.0207 301 36574.1491 48120.3368 302 54147.3844 36574.1491 303 55156.8048 54147.3844 304 60305.9911 55156.8048 305 42990.2029 60305.9911 306 55156.8048 42990.2029 307 27474.1123 55156.8048 308 33467.6100 27474.1123 309 55156.8048 33467.6100 310 55156.8048 55156.8048 311 4531.1075 55156.8048 312 55156.8048 4531.1075 313 41284.8743 55156.8048 314 47227.0903 41284.8743 315 63990.4489 47227.0903 316 55184.4267 63990.4489 317 7926.8099 55184.4267 318 50081.2561 7926.8099 319 55156.8048 50081.2561 320 55156.8048 55156.8048 321 55156.8048 55156.8048 322 55156.8048 55156.8048 323 55156.8048 55156.8048 324 55156.8048 55156.8048 325 50781.9191 55156.8048 326 55156.8048 50781.9191 327 55156.8048 55156.8048 328 55156.8048 55156.8048 329 55156.8048 55156.8048 330 55156.8048 55156.8048 331 55156.8048 55156.8048 332 55156.8048 55156.8048 333 47889.5289 55156.8048 334 54122.3903 47889.5289 335 44197.4366 54122.3903 336 -3174.7314 44197.4366 337 55156.8048 -3174.7314 338 51798.9505 55156.8048 339 46175.2185 51798.9505 340 13023.7963 46175.2185 341 55156.8048 13023.7963 342 55156.8048 55156.8048 343 55156.8048 55156.8048 344 21239.0706 55156.8048 345 55156.8048 21239.0706 346 40318.4114 55156.8048 347 55156.8048 40318.4114 348 55156.8048 55156.8048 349 13937.7951 55156.8048 350 51898.1207 13937.7951 351 44805.3156 51898.1207 352 44223.3366 44805.3156 353 47895.8151 44223.3366 354 55156.8048 47895.8151 355 55156.8048 55156.8048 356 7371.7650 55156.8048 357 55156.8048 7371.7650 358 55156.8048 55156.8048 359 55156.8048 55156.8048 360 51111.9712 55156.8048 361 55156.8048 51111.9712 362 55156.8048 55156.8048 363 -20343.5093 55156.8048 364 44817.0938 -20343.5093 365 55156.8048 44817.0938 366 55156.8048 55156.8048 367 55156.8048 55156.8048 368 55156.8048 55156.8048 369 43464.6080 55156.8048 370 55156.8048 43464.6080 371 47681.9423 55156.8048 372 39376.6041 47681.9423 373 33911.3345 39376.6041 374 55156.8048 33911.3345 375 55156.8048 55156.8048 376 28503.2373 55156.8048 377 42133.2398 28503.2373 378 849.5865 42133.2398 379 55156.8048 849.5865 380 55156.8048 55156.8048 381 7290.7240 55156.8048 382 23653.9248 7290.7240 383 55156.8048 23653.9248 384 40171.0802 55156.8048 385 -117277.3791 40171.0802 386 62596.5683 -117277.3791 387 25732.2345 62596.5683 388 49878.5072 25732.2345 389 30395.3699 49878.5072 390 40148.4756 30395.3699 391 18137.5424 40148.4756 392 -19076.5837 18137.5424 393 74584.5603 -19076.5837 394 -96716.4176 74584.5603 395 6467.0060 -96716.4176 396 -138660.7052 6467.0060 397 -189544.2257 -138660.7052 398 12385.3742 -189544.2257 399 -55495.2182 12385.3742 400 2436.6634 -55495.2182 401 42617.6282 2436.6634 402 19587.8792 42617.6282 403 -49563.0225 19587.8792 404 -8991.9428 -49563.0225 405 14014.8986 -8991.9428 406 20621.8119 14014.8986 407 13253.5098 20621.8119 408 22909.7443 13253.5098 409 -61897.0645 22909.7443 410 -51363.8727 -61897.0645 411 -161762.7527 -51363.8727 412 -25904.1719 -161762.7527 413 60090.8563 -25904.1719 414 -177521.3404 60090.8563 415 -99190.4423 -177521.3404 416 91986.7555 -99190.4423 417 -17377.1556 91986.7555 418 -50407.0390 -17377.1556 419 -37399.3791 -50407.0390 420 -44260.1342 -37399.3791 421 -175397.0259 -44260.1342 422 -547868.0565 -175397.0259 423 -458554.3229 -547868.0565 424 -48841.8151 -458554.3229 425 -250017.3396 -48841.8151 426 -311854.7009 -250017.3396 427 -1135952.0870 -311854.7009 428 -272140.7589 -1135952.0870 429 -336046.0989 -272140.7589 430 -682476.5187 -336046.0989 431 NA -682476.5187 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2923781.5148 1389642.9555 [2,] 1363354.5309 2923781.5148 [3,] -2198536.5879 1363354.5309 [4,] -394043.7941 -2198536.5879 [5,] 109880.6241 -394043.7941 [6,] -135942.2407 109880.6241 [7,] 128929.0956 -135942.2407 [8,] -336335.1695 128929.0956 [9,] 715090.7022 -336335.1695 [10,] -150648.2836 715090.7022 [11,] 342304.4900 -150648.2836 [12,] 449640.5419 342304.4900 [13,] -1139733.8466 449640.5419 [14,] -214270.3673 -1139733.8466 [15,] 460576.1545 -214270.3673 [16,] -237908.5126 460576.1545 [17,] -238086.3678 -237908.5126 [18,] -295879.7711 -238086.3678 [19,] -585121.5240 -295879.7711 [20,] 203466.5526 -585121.5240 [21,] -137406.0789 203466.5526 [22,] -119601.0233 -137406.0789 [23,] -638071.4974 -119601.0233 [24,] -318325.0883 -638071.4974 [25,] -84290.3138 -318325.0883 [26,] -471904.3808 -84290.3138 [27,] -340382.1733 -471904.3808 [28,] 348304.7876 -340382.1733 [29,] -261063.8100 348304.7876 [30,] 286671.4176 -261063.8100 [31,] -87934.3094 286671.4176 [32,] 255923.1787 -87934.3094 [33,] 130741.0876 255923.1787 [34,] -54935.4616 130741.0876 [35,] -569154.1520 -54935.4616 [36,] 34312.4748 -569154.1520 [37,] 157831.6796 34312.4748 [38,] -51019.2550 157831.6796 [39,] -304575.1497 -51019.2550 [40,] -459635.3342 -304575.1497 [41,] -329005.7740 -459635.3342 [42,] 168.5371 -329005.7740 [43,] -209377.2166 168.5371 [44,] -27236.7467 -209377.2166 [45,] -53785.8121 -27236.7467 [46,] -281550.6842 -53785.8121 [47,] -264265.0929 -281550.6842 [48,] -53546.3755 -264265.0929 [49,] -181996.3972 -53546.3755 [50,] 35728.4327 -181996.3972 [51,] -86627.1882 35728.4327 [52,] -228213.2479 -86627.1882 [53,] 44143.8787 -228213.2479 [54,] -394464.4083 44143.8787 [55,] 14554.8677 -394464.4083 [56,] 320398.0866 14554.8677 [57,] -22045.0506 320398.0866 [58,] -59459.7426 -22045.0506 [59,] -205533.9695 -59459.7426 [60,] 14512.1235 -205533.9695 [61,] -179051.4887 14512.1235 [62,] 7024.0939 -179051.4887 [63,] -174142.4494 7024.0939 [64,] -15772.7669 -174142.4494 [65,] -494815.7464 -15772.7669 [66,] -65001.8342 -494815.7464 [67,] -14962.1865 -65001.8342 [68,] -100729.1806 -14962.1865 [69,] -72008.4810 -100729.1806 [70,] -122294.9668 -72008.4810 [71,] 1254.8171 -122294.9668 [72,] -88695.3911 1254.8171 [73,] 15115.6276 -88695.3911 [74,] -34735.1194 15115.6276 [75,] -98719.3209 -34735.1194 [76,] -28468.6427 -98719.3209 [77,] -73250.3450 -28468.6427 [78,] -42230.1873 -73250.3450 [79,] -37806.4900 -42230.1873 [80,] 27097.0836 -37806.4900 [81,] -49658.1948 27097.0836 [82,] -82701.2953 -49658.1948 [83,] 41682.6321 -82701.2953 [84,] 24953.4226 41682.6321 [85,] -102380.0855 24953.4226 [86,] -29335.5087 -102380.0855 [87,] 1174.8496 -29335.5087 [88,] -115197.8123 1174.8496 [89,] 32519.6180 -115197.8123 [90,] -1901.5977 32519.6180 [91,] -21805.7745 -1901.5977 [92,] 69499.4547 -21805.7745 [93,] -172296.0652 69499.4547 [94,] 49811.4065 -172296.0652 [95,] 19621.4688 49811.4065 [96,] 71825.6583 19621.4688 [97,] 2142.7165 71825.6583 [98,] -41169.6913 2142.7165 [99,] -1469.7132 -41169.6913 [100,] -35892.6448 -1469.7132 [101,] -112495.2638 -35892.6448 [102,] -165507.7762 -112495.2638 [103,] -50097.8850 -165507.7762 [104,] 30071.9758 -50097.8850 [105,] -53262.0014 30071.9758 [106,] -252162.6722 -53262.0014 [107,] 131140.5689 -252162.6722 [108,] -175230.3204 131140.5689 [109,] -92234.7579 -175230.3204 [110,] 94522.0028 -92234.7579 [111,] -190791.2910 94522.0028 [112,] -165266.5841 -190791.2910 [113,] -57462.7139 -165266.5841 [114,] 86193.8813 -57462.7139 [115,] 21656.7901 86193.8813 [116,] 17528.3553 21656.7901 [117,] 72494.7632 17528.3553 [118,] -132141.0857 72494.7632 [119,] 14990.3877 -132141.0857 [120,] -288568.2348 14990.3877 [121,] 32830.8847 -288568.2348 [122,] -103092.5235 32830.8847 [123,] 84908.0487 -103092.5235 [124,] 13950.5787 84908.0487 [125,] -162466.3229 13950.5787 [126,] 55156.8048 -162466.3229 [127,] 42031.7011 55156.8048 [128,] 29693.9470 42031.7011 [129,] 55156.8048 29693.9470 [130,] 46456.5140 55156.8048 [131,] 57662.4310 46456.5140 [132,] 54147.3844 57662.4310 [133,] 31225.8242 54147.3844 [134,] 33614.4964 31225.8242 [135,] 44128.6218 33614.4964 [136,] 45700.9845 44128.6218 [137,] 53544.3097 45700.9845 [138,] 30904.7354 53544.3097 [139,] 55156.8048 30904.7354 [140,] 51357.0400 55156.8048 [141,] 55156.8048 51357.0400 [142,] 11408.9397 55156.8048 [143,] 50507.6768 11408.9397 [144,] 40173.9040 50507.6768 [145,] 43331.3385 40173.9040 [146,] 55156.8048 43331.3385 [147,] 47319.3565 55156.8048 [148,] 57040.9193 47319.3565 [149,] 54881.3370 57040.9193 [150,] 16159.7755 54881.3370 [151,] -39133.5926 16159.7755 [152,] 50109.7026 -39133.5926 [153,] 23678.9119 50109.7026 [154,] 52220.7451 23678.9119 [155,] 55156.8048 52220.7451 [156,] 55156.8048 55156.8048 [157,] 49628.8010 55156.8048 [158,] 62042.6027 49628.8010 [159,] 55156.8048 62042.6027 [160,] 55156.8048 55156.8048 [161,] 55156.8048 55156.8048 [162,] 55156.8048 55156.8048 [163,] 28173.1243 55156.8048 [164,] 3757.5271 28173.1243 [165,] 52614.6829 3757.5271 [166,] 55156.8048 52614.6829 [167,] 111229.1895 55156.8048 [168,] -76650.7953 111229.1895 [169,] 55156.8048 -76650.7953 [170,] 55156.8048 55156.8048 [171,] 47044.2411 55156.8048 [172,] 55156.8048 47044.2411 [173,] 33908.6287 55156.8048 [174,] 55156.8048 33908.6287 [175,] 44112.2506 55156.8048 [176,] -69826.6760 44112.2506 [177,] 54350.5341 -69826.6760 [178,] 55156.8048 54350.5341 [179,] 55156.8048 55156.8048 [180,] 41403.7863 55156.8048 [181,] 29922.3002 41403.7863 [182,] -17619.1547 29922.3002 [183,] -4775.3971 -17619.1547 [184,] 46151.2683 -4775.3971 [185,] 48386.6125 46151.2683 [186,] 11520.0678 48386.6125 [187,] 41700.1465 11520.0678 [188,] 13659.9087 41700.1465 [189,] 55156.8048 13659.9087 [190,] 47315.2565 55156.8048 [191,] 55156.8048 47315.2565 [192,] 43545.4751 55156.8048 [193,] 55156.8048 43545.4751 [194,] 51553.4904 55156.8048 [195,] 55156.8048 51553.4904 [196,] -132869.2446 55156.8048 [197,] 55156.8048 -132869.2446 [198,] 55156.8048 55156.8048 [199,] -50516.3690 55156.8048 [200,] 55156.8048 -50516.3690 [201,] 55156.8048 55156.8048 [202,] 55156.8048 55156.8048 [203,] 39709.4000 55156.8048 [204,] 393.0099 39709.4000 [205,] 49196.4916 393.0099 [206,] 47153.0040 49196.4916 [207,] 55156.8048 47153.0040 [208,] 51267.7733 55156.8048 [209,] 37862.5884 51267.7733 [210,] 55156.8048 37862.5884 [211,] 55156.8048 55156.8048 [212,] 57684.9961 55156.8048 [213,] 6260.9944 57684.9961 [214,] 24478.6930 6260.9944 [215,] 53137.9639 24478.6930 [216,] 48903.9612 53137.9639 [217,] 41823.0012 48903.9612 [218,] 53137.9639 41823.0012 [219,] -185.6711 53137.9639 [220,] 27393.6179 -185.6711 [221,] 55156.8048 27393.6179 [222,] 54089.8453 55156.8048 [223,] 21446.4109 54089.8453 [224,] 53769.7683 21446.4109 [225,] 25298.5915 53769.7683 [226,] 43900.1103 25298.5915 [227,] 55156.8048 43900.1103 [228,] 38724.2768 55156.8048 [229,] 55156.8048 38724.2768 [230,] 19420.9584 55156.8048 [231,] 35821.0062 19420.9584 [232,] 55156.8048 35821.0062 [233,] 20267.7317 55156.8048 [234,] 55156.8048 20267.7317 [235,] 55183.7008 55156.8048 [236,] 55156.8048 55183.7008 [237,] 48928.3830 55156.8048 [238,] 55156.8048 48928.3830 [239,] 55156.8048 55156.8048 [240,] 55156.8048 55156.8048 [241,] 46562.3377 55156.8048 [242,] 55156.8048 46562.3377 [243,] 49234.2381 55156.8048 [244,] -7954.4067 49234.2381 [245,] 55156.8048 -7954.4067 [246,] 55156.8048 55156.8048 [247,] 29034.0780 55156.8048 [248,] 55156.8048 29034.0780 [249,] 55156.8048 55156.8048 [250,] 34794.7895 55156.8048 [251,] 55156.8048 34794.7895 [252,] 36884.1340 55156.8048 [253,] -27169.9018 36884.1340 [254,] -149947.4502 -27169.9018 [255,] 55156.8048 -149947.4502 [256,] 43489.2216 55156.8048 [257,] 24397.7124 43489.2216 [258,] 55156.8048 24397.7124 [259,] 55156.8048 55156.8048 [260,] 55156.8048 55156.8048 [261,] 55156.8048 55156.8048 [262,] 46606.2917 55156.8048 [263,] 57802.5942 46606.2917 [264,] 55156.8048 57802.5942 [265,] 55156.8048 55156.8048 [266,] 55156.8048 55156.8048 [267,] 81171.6505 55156.8048 [268,] 55156.8048 81171.6505 [269,] 27645.6970 55156.8048 [270,] 55156.8048 27645.6970 [271,] 55156.8048 55156.8048 [272,] 55156.8048 55156.8048 [273,] 27924.6122 55156.8048 [274,] 55156.8048 27924.6122 [275,] 38755.0091 55156.8048 [276,] 24164.1677 38755.0091 [277,] 55156.8048 24164.1677 [278,] -37343.0787 55156.8048 [279,] 26865.2489 -37343.0787 [280,] 55156.8048 26865.2489 [281,] 55156.8048 55156.8048 [282,] 27805.5139 55156.8048 [283,] 55156.8048 27805.5139 [284,] 38856.9030 55156.8048 [285,] 55156.8048 38856.9030 [286,] 40708.7327 55156.8048 [287,] 52107.9101 40708.7327 [288,] 55156.8048 52107.9101 [289,] 23959.9458 55156.8048 [290,] 55156.8048 23959.9458 [291,] 45880.6314 55156.8048 [292,] 44579.9973 45880.6314 [293,] 55156.8048 44579.9973 [294,] 32828.8978 55156.8048 [295,] 65375.5583 32828.8978 [296,] 55156.8048 65375.5583 [297,] 55156.8048 55156.8048 [298,] 55156.8048 55156.8048 [299,] 46072.0207 55156.8048 [300,] 48120.3368 46072.0207 [301,] 36574.1491 48120.3368 [302,] 54147.3844 36574.1491 [303,] 55156.8048 54147.3844 [304,] 60305.9911 55156.8048 [305,] 42990.2029 60305.9911 [306,] 55156.8048 42990.2029 [307,] 27474.1123 55156.8048 [308,] 33467.6100 27474.1123 [309,] 55156.8048 33467.6100 [310,] 55156.8048 55156.8048 [311,] 4531.1075 55156.8048 [312,] 55156.8048 4531.1075 [313,] 41284.8743 55156.8048 [314,] 47227.0903 41284.8743 [315,] 63990.4489 47227.0903 [316,] 55184.4267 63990.4489 [317,] 7926.8099 55184.4267 [318,] 50081.2561 7926.8099 [319,] 55156.8048 50081.2561 [320,] 55156.8048 55156.8048 [321,] 55156.8048 55156.8048 [322,] 55156.8048 55156.8048 [323,] 55156.8048 55156.8048 [324,] 55156.8048 55156.8048 [325,] 50781.9191 55156.8048 [326,] 55156.8048 50781.9191 [327,] 55156.8048 55156.8048 [328,] 55156.8048 55156.8048 [329,] 55156.8048 55156.8048 [330,] 55156.8048 55156.8048 [331,] 55156.8048 55156.8048 [332,] 55156.8048 55156.8048 [333,] 47889.5289 55156.8048 [334,] 54122.3903 47889.5289 [335,] 44197.4366 54122.3903 [336,] -3174.7314 44197.4366 [337,] 55156.8048 -3174.7314 [338,] 51798.9505 55156.8048 [339,] 46175.2185 51798.9505 [340,] 13023.7963 46175.2185 [341,] 55156.8048 13023.7963 [342,] 55156.8048 55156.8048 [343,] 55156.8048 55156.8048 [344,] 21239.0706 55156.8048 [345,] 55156.8048 21239.0706 [346,] 40318.4114 55156.8048 [347,] 55156.8048 40318.4114 [348,] 55156.8048 55156.8048 [349,] 13937.7951 55156.8048 [350,] 51898.1207 13937.7951 [351,] 44805.3156 51898.1207 [352,] 44223.3366 44805.3156 [353,] 47895.8151 44223.3366 [354,] 55156.8048 47895.8151 [355,] 55156.8048 55156.8048 [356,] 7371.7650 55156.8048 [357,] 55156.8048 7371.7650 [358,] 55156.8048 55156.8048 [359,] 55156.8048 55156.8048 [360,] 51111.9712 55156.8048 [361,] 55156.8048 51111.9712 [362,] 55156.8048 55156.8048 [363,] -20343.5093 55156.8048 [364,] 44817.0938 -20343.5093 [365,] 55156.8048 44817.0938 [366,] 55156.8048 55156.8048 [367,] 55156.8048 55156.8048 [368,] 55156.8048 55156.8048 [369,] 43464.6080 55156.8048 [370,] 55156.8048 43464.6080 [371,] 47681.9423 55156.8048 [372,] 39376.6041 47681.9423 [373,] 33911.3345 39376.6041 [374,] 55156.8048 33911.3345 [375,] 55156.8048 55156.8048 [376,] 28503.2373 55156.8048 [377,] 42133.2398 28503.2373 [378,] 849.5865 42133.2398 [379,] 55156.8048 849.5865 [380,] 55156.8048 55156.8048 [381,] 7290.7240 55156.8048 [382,] 23653.9248 7290.7240 [383,] 55156.8048 23653.9248 [384,] 40171.0802 55156.8048 [385,] -117277.3791 40171.0802 [386,] 62596.5683 -117277.3791 [387,] 25732.2345 62596.5683 [388,] 49878.5072 25732.2345 [389,] 30395.3699 49878.5072 [390,] 40148.4756 30395.3699 [391,] 18137.5424 40148.4756 [392,] -19076.5837 18137.5424 [393,] 74584.5603 -19076.5837 [394,] -96716.4176 74584.5603 [395,] 6467.0060 -96716.4176 [396,] -138660.7052 6467.0060 [397,] -189544.2257 -138660.7052 [398,] 12385.3742 -189544.2257 [399,] -55495.2182 12385.3742 [400,] 2436.6634 -55495.2182 [401,] 42617.6282 2436.6634 [402,] 19587.8792 42617.6282 [403,] -49563.0225 19587.8792 [404,] -8991.9428 -49563.0225 [405,] 14014.8986 -8991.9428 [406,] 20621.8119 14014.8986 [407,] 13253.5098 20621.8119 [408,] 22909.7443 13253.5098 [409,] -61897.0645 22909.7443 [410,] -51363.8727 -61897.0645 [411,] -161762.7527 -51363.8727 [412,] -25904.1719 -161762.7527 [413,] 60090.8563 -25904.1719 [414,] -177521.3404 60090.8563 [415,] -99190.4423 -177521.3404 [416,] 91986.7555 -99190.4423 [417,] -17377.1556 91986.7555 [418,] -50407.0390 -17377.1556 [419,] -37399.3791 -50407.0390 [420,] -44260.1342 -37399.3791 [421,] -175397.0259 -44260.1342 [422,] -547868.0565 -175397.0259 [423,] -458554.3229 -547868.0565 [424,] -48841.8151 -458554.3229 [425,] -250017.3396 -48841.8151 [426,] -311854.7009 -250017.3396 [427,] -1135952.0870 -311854.7009 [428,] -272140.7589 -1135952.0870 [429,] -336046.0989 -272140.7589 [430,] -682476.5187 -336046.0989 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2923781.5148 1389642.9555 2 1363354.5309 2923781.5148 3 -2198536.5879 1363354.5309 4 -394043.7941 -2198536.5879 5 109880.6241 -394043.7941 6 -135942.2407 109880.6241 7 128929.0956 -135942.2407 8 -336335.1695 128929.0956 9 715090.7022 -336335.1695 10 -150648.2836 715090.7022 11 342304.4900 -150648.2836 12 449640.5419 342304.4900 13 -1139733.8466 449640.5419 14 -214270.3673 -1139733.8466 15 460576.1545 -214270.3673 16 -237908.5126 460576.1545 17 -238086.3678 -237908.5126 18 -295879.7711 -238086.3678 19 -585121.5240 -295879.7711 20 203466.5526 -585121.5240 21 -137406.0789 203466.5526 22 -119601.0233 -137406.0789 23 -638071.4974 -119601.0233 24 -318325.0883 -638071.4974 25 -84290.3138 -318325.0883 26 -471904.3808 -84290.3138 27 -340382.1733 -471904.3808 28 348304.7876 -340382.1733 29 -261063.8100 348304.7876 30 286671.4176 -261063.8100 31 -87934.3094 286671.4176 32 255923.1787 -87934.3094 33 130741.0876 255923.1787 34 -54935.4616 130741.0876 35 -569154.1520 -54935.4616 36 34312.4748 -569154.1520 37 157831.6796 34312.4748 38 -51019.2550 157831.6796 39 -304575.1497 -51019.2550 40 -459635.3342 -304575.1497 41 -329005.7740 -459635.3342 42 168.5371 -329005.7740 43 -209377.2166 168.5371 44 -27236.7467 -209377.2166 45 -53785.8121 -27236.7467 46 -281550.6842 -53785.8121 47 -264265.0929 -281550.6842 48 -53546.3755 -264265.0929 49 -181996.3972 -53546.3755 50 35728.4327 -181996.3972 51 -86627.1882 35728.4327 52 -228213.2479 -86627.1882 53 44143.8787 -228213.2479 54 -394464.4083 44143.8787 55 14554.8677 -394464.4083 56 320398.0866 14554.8677 57 -22045.0506 320398.0866 58 -59459.7426 -22045.0506 59 -205533.9695 -59459.7426 60 14512.1235 -205533.9695 61 -179051.4887 14512.1235 62 7024.0939 -179051.4887 63 -174142.4494 7024.0939 64 -15772.7669 -174142.4494 65 -494815.7464 -15772.7669 66 -65001.8342 -494815.7464 67 -14962.1865 -65001.8342 68 -100729.1806 -14962.1865 69 -72008.4810 -100729.1806 70 -122294.9668 -72008.4810 71 1254.8171 -122294.9668 72 -88695.3911 1254.8171 73 15115.6276 -88695.3911 74 -34735.1194 15115.6276 75 -98719.3209 -34735.1194 76 -28468.6427 -98719.3209 77 -73250.3450 -28468.6427 78 -42230.1873 -73250.3450 79 -37806.4900 -42230.1873 80 27097.0836 -37806.4900 81 -49658.1948 27097.0836 82 -82701.2953 -49658.1948 83 41682.6321 -82701.2953 84 24953.4226 41682.6321 85 -102380.0855 24953.4226 86 -29335.5087 -102380.0855 87 1174.8496 -29335.5087 88 -115197.8123 1174.8496 89 32519.6180 -115197.8123 90 -1901.5977 32519.6180 91 -21805.7745 -1901.5977 92 69499.4547 -21805.7745 93 -172296.0652 69499.4547 94 49811.4065 -172296.0652 95 19621.4688 49811.4065 96 71825.6583 19621.4688 97 2142.7165 71825.6583 98 -41169.6913 2142.7165 99 -1469.7132 -41169.6913 100 -35892.6448 -1469.7132 101 -112495.2638 -35892.6448 102 -165507.7762 -112495.2638 103 -50097.8850 -165507.7762 104 30071.9758 -50097.8850 105 -53262.0014 30071.9758 106 -252162.6722 -53262.0014 107 131140.5689 -252162.6722 108 -175230.3204 131140.5689 109 -92234.7579 -175230.3204 110 94522.0028 -92234.7579 111 -190791.2910 94522.0028 112 -165266.5841 -190791.2910 113 -57462.7139 -165266.5841 114 86193.8813 -57462.7139 115 21656.7901 86193.8813 116 17528.3553 21656.7901 117 72494.7632 17528.3553 118 -132141.0857 72494.7632 119 14990.3877 -132141.0857 120 -288568.2348 14990.3877 121 32830.8847 -288568.2348 122 -103092.5235 32830.8847 123 84908.0487 -103092.5235 124 13950.5787 84908.0487 125 -162466.3229 13950.5787 126 55156.8048 -162466.3229 127 42031.7011 55156.8048 128 29693.9470 42031.7011 129 55156.8048 29693.9470 130 46456.5140 55156.8048 131 57662.4310 46456.5140 132 54147.3844 57662.4310 133 31225.8242 54147.3844 134 33614.4964 31225.8242 135 44128.6218 33614.4964 136 45700.9845 44128.6218 137 53544.3097 45700.9845 138 30904.7354 53544.3097 139 55156.8048 30904.7354 140 51357.0400 55156.8048 141 55156.8048 51357.0400 142 11408.9397 55156.8048 143 50507.6768 11408.9397 144 40173.9040 50507.6768 145 43331.3385 40173.9040 146 55156.8048 43331.3385 147 47319.3565 55156.8048 148 57040.9193 47319.3565 149 54881.3370 57040.9193 150 16159.7755 54881.3370 151 -39133.5926 16159.7755 152 50109.7026 -39133.5926 153 23678.9119 50109.7026 154 52220.7451 23678.9119 155 55156.8048 52220.7451 156 55156.8048 55156.8048 157 49628.8010 55156.8048 158 62042.6027 49628.8010 159 55156.8048 62042.6027 160 55156.8048 55156.8048 161 55156.8048 55156.8048 162 55156.8048 55156.8048 163 28173.1243 55156.8048 164 3757.5271 28173.1243 165 52614.6829 3757.5271 166 55156.8048 52614.6829 167 111229.1895 55156.8048 168 -76650.7953 111229.1895 169 55156.8048 -76650.7953 170 55156.8048 55156.8048 171 47044.2411 55156.8048 172 55156.8048 47044.2411 173 33908.6287 55156.8048 174 55156.8048 33908.6287 175 44112.2506 55156.8048 176 -69826.6760 44112.2506 177 54350.5341 -69826.6760 178 55156.8048 54350.5341 179 55156.8048 55156.8048 180 41403.7863 55156.8048 181 29922.3002 41403.7863 182 -17619.1547 29922.3002 183 -4775.3971 -17619.1547 184 46151.2683 -4775.3971 185 48386.6125 46151.2683 186 11520.0678 48386.6125 187 41700.1465 11520.0678 188 13659.9087 41700.1465 189 55156.8048 13659.9087 190 47315.2565 55156.8048 191 55156.8048 47315.2565 192 43545.4751 55156.8048 193 55156.8048 43545.4751 194 51553.4904 55156.8048 195 55156.8048 51553.4904 196 -132869.2446 55156.8048 197 55156.8048 -132869.2446 198 55156.8048 55156.8048 199 -50516.3690 55156.8048 200 55156.8048 -50516.3690 201 55156.8048 55156.8048 202 55156.8048 55156.8048 203 39709.4000 55156.8048 204 393.0099 39709.4000 205 49196.4916 393.0099 206 47153.0040 49196.4916 207 55156.8048 47153.0040 208 51267.7733 55156.8048 209 37862.5884 51267.7733 210 55156.8048 37862.5884 211 55156.8048 55156.8048 212 57684.9961 55156.8048 213 6260.9944 57684.9961 214 24478.6930 6260.9944 215 53137.9639 24478.6930 216 48903.9612 53137.9639 217 41823.0012 48903.9612 218 53137.9639 41823.0012 219 -185.6711 53137.9639 220 27393.6179 -185.6711 221 55156.8048 27393.6179 222 54089.8453 55156.8048 223 21446.4109 54089.8453 224 53769.7683 21446.4109 225 25298.5915 53769.7683 226 43900.1103 25298.5915 227 55156.8048 43900.1103 228 38724.2768 55156.8048 229 55156.8048 38724.2768 230 19420.9584 55156.8048 231 35821.0062 19420.9584 232 55156.8048 35821.0062 233 20267.7317 55156.8048 234 55156.8048 20267.7317 235 55183.7008 55156.8048 236 55156.8048 55183.7008 237 48928.3830 55156.8048 238 55156.8048 48928.3830 239 55156.8048 55156.8048 240 55156.8048 55156.8048 241 46562.3377 55156.8048 242 55156.8048 46562.3377 243 49234.2381 55156.8048 244 -7954.4067 49234.2381 245 55156.8048 -7954.4067 246 55156.8048 55156.8048 247 29034.0780 55156.8048 248 55156.8048 29034.0780 249 55156.8048 55156.8048 250 34794.7895 55156.8048 251 55156.8048 34794.7895 252 36884.1340 55156.8048 253 -27169.9018 36884.1340 254 -149947.4502 -27169.9018 255 55156.8048 -149947.4502 256 43489.2216 55156.8048 257 24397.7124 43489.2216 258 55156.8048 24397.7124 259 55156.8048 55156.8048 260 55156.8048 55156.8048 261 55156.8048 55156.8048 262 46606.2917 55156.8048 263 57802.5942 46606.2917 264 55156.8048 57802.5942 265 55156.8048 55156.8048 266 55156.8048 55156.8048 267 81171.6505 55156.8048 268 55156.8048 81171.6505 269 27645.6970 55156.8048 270 55156.8048 27645.6970 271 55156.8048 55156.8048 272 55156.8048 55156.8048 273 27924.6122 55156.8048 274 55156.8048 27924.6122 275 38755.0091 55156.8048 276 24164.1677 38755.0091 277 55156.8048 24164.1677 278 -37343.0787 55156.8048 279 26865.2489 -37343.0787 280 55156.8048 26865.2489 281 55156.8048 55156.8048 282 27805.5139 55156.8048 283 55156.8048 27805.5139 284 38856.9030 55156.8048 285 55156.8048 38856.9030 286 40708.7327 55156.8048 287 52107.9101 40708.7327 288 55156.8048 52107.9101 289 23959.9458 55156.8048 290 55156.8048 23959.9458 291 45880.6314 55156.8048 292 44579.9973 45880.6314 293 55156.8048 44579.9973 294 32828.8978 55156.8048 295 65375.5583 32828.8978 296 55156.8048 65375.5583 297 55156.8048 55156.8048 298 55156.8048 55156.8048 299 46072.0207 55156.8048 300 48120.3368 46072.0207 301 36574.1491 48120.3368 302 54147.3844 36574.1491 303 55156.8048 54147.3844 304 60305.9911 55156.8048 305 42990.2029 60305.9911 306 55156.8048 42990.2029 307 27474.1123 55156.8048 308 33467.6100 27474.1123 309 55156.8048 33467.6100 310 55156.8048 55156.8048 311 4531.1075 55156.8048 312 55156.8048 4531.1075 313 41284.8743 55156.8048 314 47227.0903 41284.8743 315 63990.4489 47227.0903 316 55184.4267 63990.4489 317 7926.8099 55184.4267 318 50081.2561 7926.8099 319 55156.8048 50081.2561 320 55156.8048 55156.8048 321 55156.8048 55156.8048 322 55156.8048 55156.8048 323 55156.8048 55156.8048 324 55156.8048 55156.8048 325 50781.9191 55156.8048 326 55156.8048 50781.9191 327 55156.8048 55156.8048 328 55156.8048 55156.8048 329 55156.8048 55156.8048 330 55156.8048 55156.8048 331 55156.8048 55156.8048 332 55156.8048 55156.8048 333 47889.5289 55156.8048 334 54122.3903 47889.5289 335 44197.4366 54122.3903 336 -3174.7314 44197.4366 337 55156.8048 -3174.7314 338 51798.9505 55156.8048 339 46175.2185 51798.9505 340 13023.7963 46175.2185 341 55156.8048 13023.7963 342 55156.8048 55156.8048 343 55156.8048 55156.8048 344 21239.0706 55156.8048 345 55156.8048 21239.0706 346 40318.4114 55156.8048 347 55156.8048 40318.4114 348 55156.8048 55156.8048 349 13937.7951 55156.8048 350 51898.1207 13937.7951 351 44805.3156 51898.1207 352 44223.3366 44805.3156 353 47895.8151 44223.3366 354 55156.8048 47895.8151 355 55156.8048 55156.8048 356 7371.7650 55156.8048 357 55156.8048 7371.7650 358 55156.8048 55156.8048 359 55156.8048 55156.8048 360 51111.9712 55156.8048 361 55156.8048 51111.9712 362 55156.8048 55156.8048 363 -20343.5093 55156.8048 364 44817.0938 -20343.5093 365 55156.8048 44817.0938 366 55156.8048 55156.8048 367 55156.8048 55156.8048 368 55156.8048 55156.8048 369 43464.6080 55156.8048 370 55156.8048 43464.6080 371 47681.9423 55156.8048 372 39376.6041 47681.9423 373 33911.3345 39376.6041 374 55156.8048 33911.3345 375 55156.8048 55156.8048 376 28503.2373 55156.8048 377 42133.2398 28503.2373 378 849.5865 42133.2398 379 55156.8048 849.5865 380 55156.8048 55156.8048 381 7290.7240 55156.8048 382 23653.9248 7290.7240 383 55156.8048 23653.9248 384 40171.0802 55156.8048 385 -117277.3791 40171.0802 386 62596.5683 -117277.3791 387 25732.2345 62596.5683 388 49878.5072 25732.2345 389 30395.3699 49878.5072 390 40148.4756 30395.3699 391 18137.5424 40148.4756 392 -19076.5837 18137.5424 393 74584.5603 -19076.5837 394 -96716.4176 74584.5603 395 6467.0060 -96716.4176 396 -138660.7052 6467.0060 397 -189544.2257 -138660.7052 398 12385.3742 -189544.2257 399 -55495.2182 12385.3742 400 2436.6634 -55495.2182 401 42617.6282 2436.6634 402 19587.8792 42617.6282 403 -49563.0225 19587.8792 404 -8991.9428 -49563.0225 405 14014.8986 -8991.9428 406 20621.8119 14014.8986 407 13253.5098 20621.8119 408 22909.7443 13253.5098 409 -61897.0645 22909.7443 410 -51363.8727 -61897.0645 411 -161762.7527 -51363.8727 412 -25904.1719 -161762.7527 413 60090.8563 -25904.1719 414 -177521.3404 60090.8563 415 -99190.4423 -177521.3404 416 91986.7555 -99190.4423 417 -17377.1556 91986.7555 418 -50407.0390 -17377.1556 419 -37399.3791 -50407.0390 420 -44260.1342 -37399.3791 421 -175397.0259 -44260.1342 422 -547868.0565 -175397.0259 423 -458554.3229 -547868.0565 424 -48841.8151 -458554.3229 425 -250017.3396 -48841.8151 426 -311854.7009 -250017.3396 427 -1135952.0870 -311854.7009 428 -272140.7589 -1135952.0870 429 -336046.0989 -272140.7589 430 -682476.5187 -336046.0989 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7qdsz1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/8qdsz1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/904921292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/1004921292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11m5qq1292936583.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12p5ow1292936583.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13w73p1292936583.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14hpkv1292936583.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/15sg1y1292936583.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/16oqh71292936583.tab") + } > > try(system("convert tmp/1ulcq1292936583.ps tmp/1ulcq1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/2mvcb1292936583.ps tmp/2mvcb1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/3mvcb1292936583.ps tmp/3mvcb1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/4mvcb1292936583.ps tmp/4mvcb1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/5f4bw1292936583.ps tmp/5f4bw1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/6f4bw1292936583.ps tmp/6f4bw1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/7qdsz1292936583.ps tmp/7qdsz1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/8qdsz1292936583.ps tmp/8qdsz1292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/904921292936583.ps tmp/904921292936583.png",intern=TRUE)) character(0) > try(system("convert tmp/1004921292936583.ps tmp/1004921292936583.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 11.661 2.210 23.947