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 = '2' > #'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 Orders Costs Dividends Wealth 1 807 162556 213118 6282154 2 444 29790 81767 4321023 3 412 87550 153198 4111912 4 428 84738 -26007 223193 5 315 54660 126942 1491348 6 168 42634 157214 1629616 7 263 40949 129352 1398893 8 267 45187 234817 1926517 9 228 37704 60448 983660 10 129 16275 47818 1443586 11 104 25830 245546 1073089 12 122 12679 48020 984885 13 393 18014 -1710 1405225 14 190 43556 32648 227132 15 280 24811 95350 929118 16 63 6575 151352 1071292 17 102 7123 288170 638830 18 265 21950 114337 856956 19 234 37597 37884 992426 20 277 17821 122844 444477 21 73 12988 82340 857217 22 67 22330 79801 711969 23 103 13326 165548 702380 24 290 16189 116384 358589 25 83 7146 134028 297978 26 56 15824 63838 585715 27 236 27664 74996 657954 28 73 11920 31080 209458 29 34 8568 32168 786690 30 139 14416 49857 439798 31 26 3369 87161 688779 32 70 11819 106113 574339 33 40 6984 80570 741409 34 42 4519 102129 597793 35 12 2220 301670 644190 36 211 18562 102313 377934 37 74 10327 88577 640273 38 80 5336 112477 697458 39 83 2365 191778 550608 40 131 4069 79804 207393 41 203 8636 128294 301607 42 56 13718 96448 345783 43 89 4525 93811 501749 44 88 6869 117520 379983 45 39 4628 69159 387475 46 25 3689 101792 377305 47 49 4891 210568 370837 48 149 7489 136996 430866 49 58 4901 121920 469107 50 41 2284 76403 194493 51 90 3160 108094 530670 52 136 4150 134759 518365 53 97 7285 188873 491303 54 63 1134 146216 527021 55 114 4658 156608 233773 56 77 2384 61348 405972 57 6 3748 50350 652925 58 47 5371 87720 446211 59 51 1285 99489 341340 60 85 9327 87419 387699 61 43 5565 94355 493408 62 32 1528 60326 146494 63 25 3122 94670 414462 64 77 7561 82425 364304 65 54 2675 59017 355178 66 251 13253 90829 357760 67 15 880 80791 261216 68 44 2053 100423 397144 69 73 1424 131116 374943 70 85 4036 100269 424898 71 49 3045 27330 202055 72 38 5119 39039 378525 73 35 1431 106885 310768 74 9 554 79285 325738 75 34 1975 118881 394510 76 20 1765 77623 247060 77 29 1012 114768 368078 78 11 810 74015 236761 79 52 1280 69465 312378 80 13 666 117869 339836 81 29 1380 60982 347385 82 66 4677 90131 426280 83 33 876 138971 352850 84 15 814 39625 301881 85 15 514 102725 377516 86 68 5692 64239 357312 87 100 3642 90262 458343 88 13 540 103960 354228 89 45 2099 106611 308636 90 14 567 103345 386212 91 36 2001 95551 393343 92 40 2949 82903 378509 93 68 2253 63593 452469 94 29 6533 126910 364839 95 43 1889 37527 358649 96 30 3055 60247 376641 97 9 272 112995 429112 98 22 1414 70184 330546 99 19 2564 130140 403560 100 9 1383 73221 317892 101 31 1261 76114 307528 102 19 975 90534 235133 103 55 3366 108479 299243 104 8 576 113761 314073 105 28 1686 68696 368186 106 29 746 71561 269661 107 48 3192 59831 125390 108 16 2045 97890 510834 109 47 5702 101481 321896 110 20 1932 72954 249898 111 22 936 67939 408881 112 33 3437 48022 158492 113 44 5131 86111 292154 114 13 2397 74020 289513 115 6 1389 57530 378049 116 35 1503 56364 343466 117 8 402 84990 332743 118 17 2239 88590 442882 119 11 2234 77200 214215 120 21 837 61262 315688 121 92 10579 110309 375195 122 12 875 67000 334280 123 112 1585 93099 355864 124 25 1659 107577 480382 125 17 2647 62920 353058 126 23 3294 75832 217193 127 0 0 60720 315380 128 10 94 60793 314533 129 23 422 57935 318056 130 0 0 60720 315380 131 7 34 60630 314353 132 25 1558 55637 369448 133 1 0 60720 315380 134 20 43 60887 312846 135 4 645 60720 312075 136 4 316 60505 315009 137 10 115 60945 318903 138 1 5 60720 314887 139 4 897 60720 314913 140 0 0 60720 315380 141 8 389 58990 325506 142 0 0 60720 315380 143 11 1002 56750 298568 144 4 36 60894 315834 145 15 460 63346 329784 146 9 309 56535 312878 147 0 0 60720 315380 148 7 9 60835 314987 149 2 271 60720 325249 150 0 14 61016 315877 151 7 520 58650 291650 152 46 1766 60438 305959 153 5 0 60720 315380 154 7 458 58625 297765 155 2 20 60938 315245 156 0 0 60720 315380 157 0 0 60720 315380 158 2 98 61490 315236 159 5 405 60845 336425 160 0 0 60720 315380 161 0 0 60720 315380 162 0 0 60720 315380 163 0 0 60720 315380 164 7 483 60830 306268 165 24 454 63261 302187 166 1 47 60720 314882 167 0 0 60720 315380 168 18 757 45689 382712 169 55 4655 60720 341570 170 0 0 60720 315380 171 0 0 60720 315380 172 3 36 61564 312412 173 0 0 60720 315380 174 9 203 61938 309596 175 0 0 60720 315380 176 8 126 60951 315547 177 113 400 60720 313267 178 0 71 60745 316176 179 0 0 60720 315380 180 0 0 60720 315380 181 19 972 71642 359335 182 11 531 71641 330068 183 25 2461 55792 314289 184 16 378 71873 297413 185 5 23 62555 314806 186 11 638 60370 333210 187 23 2300 64873 352108 188 6 149 62041 313332 189 5 226 65745 291787 190 0 0 60720 315380 191 7 275 59500 318745 192 0 0 60720 315380 193 7 141 61630 315366 194 0 0 60720 315380 195 3 28 60890 315688 196 0 0 60720 315380 197 89 4980 113521 409642 198 0 0 60720 315380 199 0 0 60720 315380 200 19 472 80045 269587 201 0 0 60720 315380 202 0 0 60720 315380 203 0 0 60720 315380 204 12 203 50804 300962 205 12 496 87390 325479 206 5 10 61656 316155 207 2 63 65688 318574 208 0 0 60720 315380 209 26 1136 48522 343613 210 3 265 60720 306948 211 0 0 60720 315380 212 0 0 60720 315380 213 11 267 57640 330059 214 10 474 61977 288985 215 5 534 62620 304485 216 2 0 60720 315380 217 6 15 60831 315688 218 7 397 60646 317736 219 2 0 60720 315380 220 28 1866 56225 322331 221 3 288 60510 296656 222 0 0 60720 315380 223 1 3 60698 315354 224 20 468 60720 312161 225 1 20 60805 315576 226 22 278 61404 314922 227 9 61 60720 314551 228 0 0 60720 315380 229 2 192 65276 312339 230 0 0 60720 315380 231 7 317 63915 298700 232 9 738 60720 321376 233 0 0 60720 315380 234 13 368 61686 303230 235 0 0 60720 315380 236 0 2 60743 315487 237 0 0 60720 315380 238 6 53 60349 315793 239 0 0 60720 315380 240 0 0 60720 315380 241 0 0 60720 315380 242 3 94 61360 312887 243 0 0 60720 315380 244 7 24 59818 315637 245 2 2332 72680 324385 246 0 0 60720 315380 247 0 0 60720 315380 248 15 131 61808 308989 249 0 0 60720 315380 250 0 0 60720 315380 251 9 206 53110 296702 252 0 0 60720 315380 253 1 167 64245 307322 254 38 622 73007 304376 255 57 2328 82732 253588 256 0 0 60720 315380 257 7 365 54820 309560 258 26 364 47705 298466 259 0 0 60720 315380 260 0 0 60720 315380 261 0 0 60720 315380 262 0 0 60720 315380 263 13 226 72835 343929 264 10 307 58856 331955 265 0 0 60720 315380 266 0 0 60720 315380 267 0 0 60720 315380 268 9 188 77655 381180 269 0 0 60720 315380 270 26 138 69817 331420 271 0 0 60720 315380 272 0 0 60720 315380 273 0 0 60720 315380 274 19 125 60798 310201 275 0 0 60720 315380 276 12 282 62452 320016 277 23 335 64175 320398 278 0 0 60720 315380 279 29 1324 67440 291841 280 8 176 68136 310670 281 0 0 60720 315380 282 0 0 60720 315380 283 26 249 56726 313491 284 0 0 60720 315380 285 9 333 70811 331323 286 0 0 60720 315380 287 5 601 60720 319210 288 3 30 62045 318098 289 0 0 60720 315380 290 13 249 54323 292754 291 0 0 60720 315380 292 12 165 62841 325176 293 19 453 81125 365959 294 0 0 60720 315380 295 10 53 59506 302409 296 9 382 59365 340968 297 0 0 60720 315380 298 0 0 60720 315380 299 0 0 60720 315380 300 9 0 60720 315380 301 4 30 60798 313164 302 1 290 58790 301164 303 1 0 60720 315380 304 0 0 60720 315380 305 14 366 61808 344425 306 12 2 60735 315394 307 0 0 60720 315380 308 19 209 64016 316647 309 17 384 54683 309836 310 0 0 60720 315380 311 0 0 60720 315380 312 32 365 87192 346611 313 0 0 60720 315380 314 14 49 64107 322031 315 8 3 60761 315656 316 4 133 65990 339445 317 0 32 59988 314964 318 20 368 61167 297141 319 5 1 60719 315372 320 0 0 60720 315380 321 0 0 60720 315380 322 0 0 60720 315380 323 0 0 60720 315380 324 0 0 60720 315380 325 0 0 60720 315380 326 1 22 60722 312502 327 0 0 60720 315380 328 0 0 60720 315380 329 0 0 60720 315380 330 0 0 60720 315380 331 0 0 60720 315380 332 0 0 60720 315380 333 0 0 60720 315380 334 4 96 60379 313729 335 1 1 60727 315388 336 4 314 60720 315371 337 20 844 60925 296139 338 0 0 60720 315380 339 1 26 60896 313880 340 10 125 59734 317698 341 12 304 62969 295580 342 0 0 60720 315380 343 0 0 60720 315380 344 0 0 60720 315380 345 13 621 60720 308256 346 0 0 60720 315380 347 3 119 59118 303677 348 0 0 60720 315380 349 0 0 60720 315380 350 10 1595 60720 319369 351 3 312 58598 318690 352 7 60 61124 314049 353 10 587 59595 325699 354 1 135 62065 314210 355 0 0 60720 315380 356 0 0 60720 315380 357 15 514 78780 322378 358 0 0 60720 315380 359 0 0 60720 315380 360 0 0 60720 315380 361 4 1 60722 315398 362 0 0 60720 315380 363 0 0 60720 315380 364 28 1763 61600 308336 365 9 180 59635 316386 366 0 0 60720 315380 367 0 0 60720 315380 368 0 0 60720 315380 369 0 0 60720 315380 370 7 218 60720 315553 371 0 0 60720 315380 372 7 448 59781 323361 373 7 227 76644 336639 374 3 174 64820 307424 375 0 0 60720 315380 376 0 0 60720 315380 377 11 121 56178 295370 378 7 607 60436 322340 379 10 2212 60720 319864 380 0 0 60720 315380 381 0 0 60720 315380 382 18 530 73433 317291 383 14 571 41477 280398 384 0 0 60720 315380 385 12 78 62700 317330 386 29 2489 67804 238125 387 3 131 59661 327071 388 6 923 58620 309038 389 3 72 60398 314210 390 8 572 58580 307930 391 10 397 62710 322327 392 6 450 59325 292136 393 8 622 60950 263276 394 6 694 68060 367655 395 9 3425 83620 283910 396 8 562 58456 283587 397 26 4917 52811 243650 398 239 1442 121173 438493 399 7 529 63870 296261 400 41 2126 21001 230621 401 3 1061 70415 304252 402 8 776 64230 333505 403 6 611 59190 296919 404 21 1526 69351 278990 405 7 592 64270 276898 406 11 1182 70694 327007 407 11 621 68005 317046 408 12 989 58930 304555 409 9 438 58320 298096 410 3 726 69980 231861 411 57 1303 69863 309422 412 21 7419 63255 286963 413 15 1164 57320 269753 414 32 3310 75230 448243 415 11 1920 79420 165404 416 2 965 73490 204325 417 23 3256 35250 407159 418 20 1135 62285 290476 419 24 1270 69206 275311 420 1 661 65920 246541 421 1 1013 69770 253468 422 74 2844 72683 240897 423 68 11528 -14545 -83265 424 20 6526 55830 -42143 425 20 2264 55174 272713 426 82 5109 67038 215362 427 21 3999 51252 42754 428 244 35624 157278 306275 429 32 9252 79510 253537 430 86 15236 77440 372631 431 69 18073 27284 -7170 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Costs Dividends Wealth -5.927e+00 4.792e-03 2.134e-04 1.896e-05 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -130.625 -13.010 -9.128 1.221 286.324 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -5.927e+00 4.329e+00 -1.369 0.171719 Costs 4.792e-03 2.340e-04 20.477 < 2e-16 *** Dividends 2.134e-04 5.641e-05 3.783 0.000177 *** Wealth 1.896e-05 6.569e-06 2.887 0.004092 ** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 34.75 on 427 degrees of freedom Multiple R-squared: 0.7839, Adjusted R-squared: 0.7824 F-statistic: 516.3 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,] 0.9982348 3.530380e-03 1.765190e-03 [2,] 0.9998106 3.788738e-04 1.894369e-04 [3,] 0.9997402 5.195682e-04 2.597841e-04 [4,] 0.9999833 3.342067e-05 1.671033e-05 [5,] 0.9999848 3.035395e-05 1.517698e-05 [6,] 0.9999842 3.151213e-05 1.575607e-05 [7,] 1.0000000 1.588233e-11 7.941166e-12 [8,] 1.0000000 1.766135e-12 8.830675e-13 [9,] 1.0000000 1.601033e-14 8.005167e-15 [10,] 1.0000000 1.512065e-14 7.560323e-15 [11,] 1.0000000 2.447820e-16 1.223910e-16 [12,] 1.0000000 2.959264e-18 1.479632e-18 [13,] 1.0000000 4.413380e-18 2.206690e-18 [14,] 1.0000000 7.104719e-23 3.552359e-23 [15,] 1.0000000 5.737978e-25 2.868989e-25 [16,] 1.0000000 3.051290e-29 1.525645e-29 [17,] 1.0000000 5.699245e-29 2.849622e-29 [18,] 1.0000000 5.248810e-37 2.624405e-37 [19,] 1.0000000 1.519599e-36 7.597997e-37 [20,] 1.0000000 7.895915e-40 3.947958e-40 [21,] 1.0000000 1.375510e-39 6.877550e-40 [22,] 1.0000000 2.445891e-40 1.222946e-40 [23,] 1.0000000 3.075644e-43 1.537822e-43 [24,] 1.0000000 8.130189e-43 4.065095e-43 [25,] 1.0000000 8.139010e-44 4.069505e-44 [26,] 1.0000000 3.252055e-44 1.626028e-44 [27,] 1.0000000 3.212196e-45 1.606098e-45 [28,] 1.0000000 2.612578e-45 1.306289e-45 [29,] 1.0000000 2.005803e-47 1.002901e-47 [30,] 1.0000000 8.337355e-49 4.168678e-49 [31,] 1.0000000 7.471390e-49 3.735695e-49 [32,] 1.0000000 2.504666e-48 1.252333e-48 [33,] 1.0000000 6.909073e-48 3.454537e-48 [34,] 1.0000000 5.505645e-50 2.752823e-50 [35,] 1.0000000 1.777326e-55 8.886629e-56 [36,] 1.0000000 6.343819e-57 3.171909e-57 [37,] 1.0000000 1.399601e-56 6.998006e-57 [38,] 1.0000000 5.161566e-56 2.580783e-56 [39,] 1.0000000 3.925506e-56 1.962753e-56 [40,] 1.0000000 1.967231e-56 9.836156e-57 [41,] 1.0000000 8.463077e-57 4.231539e-57 [42,] 1.0000000 6.643086e-58 3.321543e-58 [43,] 1.0000000 1.537242e-57 7.686211e-58 [44,] 1.0000000 2.473981e-57 1.236991e-57 [45,] 1.0000000 4.066066e-57 2.033033e-57 [46,] 1.0000000 1.515954e-58 7.579770e-59 [47,] 1.0000000 3.948550e-58 1.974275e-58 [48,] 1.0000000 1.400942e-57 7.004708e-58 [49,] 1.0000000 4.211118e-58 2.105559e-58 [50,] 1.0000000 1.738842e-58 8.694210e-59 [51,] 1.0000000 2.188622e-60 1.094311e-60 [52,] 1.0000000 2.848182e-60 1.424091e-60 [53,] 1.0000000 5.834498e-60 2.917249e-60 [54,] 1.0000000 1.876161e-59 9.380807e-60 [55,] 1.0000000 1.472232e-59 7.361162e-60 [56,] 1.0000000 1.522912e-59 7.614560e-60 [57,] 1.0000000 1.149525e-59 5.747626e-60 [58,] 1.0000000 3.591374e-59 1.795687e-59 [59,] 1.0000000 5.351593e-59 2.675797e-59 [60,] 1.0000000 5.287666e-71 2.643833e-71 [61,] 1.0000000 6.685561e-71 3.342781e-71 [62,] 1.0000000 1.825038e-70 9.125191e-71 [63,] 1.0000000 2.055132e-70 1.027566e-70 [64,] 1.0000000 1.867579e-70 9.337895e-71 [65,] 1.0000000 7.429756e-71 3.714878e-71 [66,] 1.0000000 9.031350e-71 4.515675e-71 [67,] 1.0000000 2.420535e-70 1.210268e-70 [68,] 1.0000000 2.558882e-70 1.279441e-70 [69,] 1.0000000 6.313101e-70 3.156551e-70 [70,] 1.0000000 1.079854e-69 5.399269e-70 [71,] 1.0000000 2.785265e-69 1.392632e-69 [72,] 1.0000000 4.096927e-69 2.048464e-69 [73,] 1.0000000 3.331330e-69 1.665665e-69 [74,] 1.0000000 4.375153e-69 2.187576e-69 [75,] 1.0000000 8.081618e-69 4.040809e-69 [76,] 1.0000000 2.184961e-68 1.092481e-68 [77,] 1.0000000 6.634984e-68 3.317492e-68 [78,] 1.0000000 8.438484e-68 4.219242e-68 [79,] 1.0000000 1.392804e-67 6.964019e-68 [80,] 1.0000000 2.449459e-67 1.224729e-67 [81,] 1.0000000 1.497267e-68 7.486335e-69 [82,] 1.0000000 2.285478e-68 1.142739e-68 [83,] 1.0000000 6.951083e-68 3.475542e-68 [84,] 1.0000000 1.073633e-67 5.368165e-68 [85,] 1.0000000 3.054618e-67 1.527309e-67 [86,] 1.0000000 8.252023e-67 4.126011e-67 [87,] 1.0000000 2.693377e-67 1.346688e-67 [88,] 1.0000000 3.944519e-68 1.972259e-68 [89,] 1.0000000 2.088630e-68 1.044315e-68 [90,] 1.0000000 4.283686e-68 2.141843e-68 [91,] 1.0000000 3.363597e-68 1.681798e-68 [92,] 1.0000000 7.608381e-68 3.804191e-68 [93,] 1.0000000 2.037087e-68 1.018544e-68 [94,] 1.0000000 2.831503e-68 1.415752e-68 [95,] 1.0000000 7.457111e-68 3.728556e-68 [96,] 1.0000000 1.789345e-67 8.946726e-68 [97,] 1.0000000 6.237970e-67 3.118985e-67 [98,] 1.0000000 3.138240e-67 1.569120e-67 [99,] 1.0000000 7.578249e-67 3.789124e-67 [100,] 1.0000000 1.858245e-66 9.291226e-67 [101,] 1.0000000 2.382560e-66 1.191280e-66 [102,] 1.0000000 1.678983e-66 8.394917e-67 [103,] 1.0000000 3.634883e-66 1.817442e-66 [104,] 1.0000000 8.738220e-66 4.369110e-66 [105,] 1.0000000 2.093571e-65 1.046786e-65 [106,] 1.0000000 3.556209e-65 1.778104e-65 [107,] 1.0000000 1.017864e-64 5.089320e-65 [108,] 1.0000000 1.528917e-64 7.644584e-65 [109,] 1.0000000 2.329265e-64 1.164633e-64 [110,] 1.0000000 3.406101e-64 1.703050e-64 [111,] 1.0000000 5.203742e-64 2.601871e-64 [112,] 1.0000000 5.742088e-64 2.871044e-64 [113,] 1.0000000 8.516681e-64 4.258341e-64 [114,] 1.0000000 1.997791e-63 9.988954e-64 [115,] 1.0000000 6.168347e-63 3.084174e-63 [116,] 1.0000000 1.382126e-62 6.910629e-63 [117,] 1.0000000 3.821532e-66 1.910766e-66 [118,] 1.0000000 4.842284e-66 2.421142e-66 [119,] 1.0000000 1.006642e-65 5.033209e-66 [120,] 1.0000000 2.373082e-65 1.186541e-65 [121,] 1.0000000 4.433866e-65 2.216933e-65 [122,] 1.0000000 1.073637e-64 5.368186e-65 [123,] 1.0000000 2.109796e-64 1.054898e-64 [124,] 1.0000000 4.076799e-64 2.038400e-64 [125,] 1.0000000 9.721224e-64 4.860612e-64 [126,] 1.0000000 2.150096e-63 1.075048e-63 [127,] 1.0000000 4.450786e-63 2.225393e-63 [128,] 1.0000000 9.865985e-63 4.932993e-63 [129,] 1.0000000 2.076681e-62 1.038341e-62 [130,] 1.0000000 4.638756e-62 2.319378e-62 [131,] 1.0000000 1.177592e-61 5.887961e-62 [132,] 1.0000000 2.534487e-61 1.267243e-61 [133,] 1.0000000 5.246779e-61 2.623389e-61 [134,] 1.0000000 1.108694e-60 5.543470e-61 [135,] 1.0000000 2.758609e-60 1.379305e-60 [136,] 1.0000000 5.910029e-60 2.955014e-60 [137,] 1.0000000 1.496616e-59 7.483079e-60 [138,] 1.0000000 3.619837e-59 1.809918e-59 [139,] 1.0000000 9.720389e-59 4.860194e-59 [140,] 1.0000000 2.452340e-58 1.226170e-58 [141,] 1.0000000 5.377179e-58 2.688590e-58 [142,] 1.0000000 1.384565e-57 6.922824e-58 [143,] 1.0000000 3.107500e-57 1.553750e-57 [144,] 1.0000000 6.863367e-57 3.431683e-57 [145,] 1.0000000 1.743476e-56 8.717378e-57 [146,] 1.0000000 1.413780e-56 7.068898e-57 [147,] 1.0000000 3.575946e-56 1.787973e-56 [148,] 1.0000000 9.118855e-56 4.559428e-56 [149,] 1.0000000 2.165974e-55 1.082987e-55 [150,] 1.0000000 4.892129e-55 2.446065e-55 [151,] 1.0000000 1.108384e-54 5.541919e-55 [152,] 1.0000000 2.614059e-54 1.307030e-54 [153,] 1.0000000 6.368322e-54 3.184161e-54 [154,] 1.0000000 1.452670e-53 7.263350e-54 [155,] 1.0000000 3.320946e-53 1.660473e-53 [156,] 1.0000000 7.607087e-53 3.803544e-53 [157,] 1.0000000 1.745619e-52 8.728093e-53 [158,] 1.0000000 4.470249e-52 2.235124e-52 [159,] 1.0000000 9.923499e-52 4.961750e-52 [160,] 1.0000000 2.329486e-51 1.164743e-51 [161,] 1.0000000 5.356619e-51 2.678309e-51 [162,] 1.0000000 1.005883e-50 5.029415e-51 [163,] 1.0000000 1.368181e-50 6.840906e-51 [164,] 1.0000000 3.154729e-50 1.577364e-50 [165,] 1.0000000 7.279591e-50 3.639795e-50 [166,] 1.0000000 1.787881e-49 8.939406e-50 [167,] 1.0000000 4.127696e-49 2.063848e-49 [168,] 1.0000000 1.081393e-48 5.406964e-49 [169,] 1.0000000 2.494917e-48 1.247459e-48 [170,] 1.0000000 6.468277e-48 3.234139e-48 [171,] 1.0000000 5.120789e-55 2.560395e-55 [172,] 1.0000000 1.253448e-54 6.267241e-55 [173,] 1.0000000 3.103329e-54 1.551665e-54 [174,] 1.0000000 7.683765e-54 3.841883e-54 [175,] 1.0000000 2.241164e-53 1.120582e-53 [176,] 1.0000000 6.033518e-53 3.016759e-53 [177,] 1.0000000 1.502999e-52 7.514993e-53 [178,] 1.0000000 4.367340e-52 2.183670e-52 [179,] 1.0000000 1.173511e-51 5.867557e-52 [180,] 1.0000000 3.211067e-51 1.605534e-51 [181,] 1.0000000 9.027439e-51 4.513720e-51 [182,] 1.0000000 2.426183e-50 1.213092e-50 [183,] 1.0000000 6.232851e-50 3.116425e-50 [184,] 1.0000000 1.521008e-49 7.605039e-50 [185,] 1.0000000 4.076696e-49 2.038348e-49 [186,] 1.0000000 9.925991e-49 4.962995e-49 [187,] 1.0000000 2.666625e-48 1.333313e-48 [188,] 1.0000000 6.475185e-48 3.237593e-48 [189,] 1.0000000 1.666113e-47 8.330563e-48 [190,] 1.0000000 4.034962e-47 2.017481e-47 [191,] 1.0000000 5.044490e-47 2.522245e-47 [192,] 1.0000000 1.220001e-46 6.100006e-47 [193,] 1.0000000 2.946672e-46 1.473336e-46 [194,] 1.0000000 8.215387e-46 4.107694e-46 [195,] 1.0000000 1.974970e-45 9.874852e-46 [196,] 1.0000000 4.740755e-45 2.370378e-45 [197,] 1.0000000 1.136212e-44 5.681060e-45 [198,] 1.0000000 2.480295e-44 1.240148e-44 [199,] 1.0000000 4.843811e-44 2.421906e-44 [200,] 1.0000000 1.250937e-43 6.254684e-44 [201,] 1.0000000 2.914489e-43 1.457245e-43 [202,] 1.0000000 6.926240e-43 3.463120e-43 [203,] 1.0000000 8.825934e-43 4.412967e-43 [204,] 1.0000000 2.177320e-42 1.088660e-42 [205,] 1.0000000 5.165025e-42 2.582512e-42 [206,] 1.0000000 1.222805e-41 6.114024e-42 [207,] 1.0000000 3.022787e-41 1.511393e-41 [208,] 1.0000000 7.829506e-41 3.914753e-41 [209,] 1.0000000 1.900291e-40 9.501456e-41 [210,] 1.0000000 4.634540e-40 2.317270e-40 [211,] 1.0000000 1.174343e-39 5.871714e-40 [212,] 1.0000000 2.952201e-39 1.476101e-39 [213,] 1.0000000 7.138251e-39 3.569125e-39 [214,] 1.0000000 1.317443e-38 6.587213e-39 [215,] 1.0000000 3.164652e-38 1.582326e-38 [216,] 1.0000000 7.323336e-38 3.661668e-38 [217,] 1.0000000 1.724476e-37 8.622382e-38 [218,] 1.0000000 3.675468e-37 1.837734e-37 [219,] 1.0000000 8.592255e-37 4.296128e-37 [220,] 1.0000000 1.673135e-36 8.365675e-37 [221,] 1.0000000 4.110628e-36 2.055314e-36 [222,] 1.0000000 9.384328e-36 4.692164e-36 [223,] 1.0000000 2.065507e-35 1.032753e-35 [224,] 1.0000000 4.690526e-35 2.345263e-35 [225,] 1.0000000 1.129985e-34 5.649925e-35 [226,] 1.0000000 2.729701e-34 1.364850e-34 [227,] 1.0000000 6.143585e-34 3.071792e-34 [228,] 1.0000000 1.464666e-33 7.323331e-34 [229,] 1.0000000 3.277330e-33 1.638665e-33 [230,] 1.0000000 7.309171e-33 3.654585e-33 [231,] 1.0000000 1.626113e-32 8.130567e-33 [232,] 1.0000000 3.853634e-32 1.926817e-32 [233,] 1.0000000 8.521092e-32 4.260546e-32 [234,] 1.0000000 1.878538e-31 9.392691e-32 [235,] 1.0000000 4.128861e-31 2.064430e-31 [236,] 1.0000000 9.375418e-31 4.687709e-31 [237,] 1.0000000 2.047621e-30 1.023811e-30 [238,] 1.0000000 4.755825e-30 2.377913e-30 [239,] 1.0000000 5.007747e-30 2.503873e-30 [240,] 1.0000000 1.087612e-29 5.438061e-30 [241,] 1.0000000 2.354563e-29 1.177282e-29 [242,] 1.0000000 5.181283e-29 2.590641e-29 [243,] 1.0000000 1.114833e-28 5.574166e-29 [244,] 1.0000000 2.390853e-28 1.195427e-28 [245,] 1.0000000 5.140044e-28 2.570022e-28 [246,] 1.0000000 1.095793e-27 5.478963e-28 [247,] 1.0000000 2.239302e-27 1.119651e-27 [248,] 1.0000000 3.268411e-27 1.634206e-27 [249,] 1.0000000 3.782429e-27 1.891215e-27 [250,] 1.0000000 8.017740e-27 4.008870e-27 [251,] 1.0000000 1.774164e-26 8.870821e-27 [252,] 1.0000000 1.720013e-26 8.600067e-27 [253,] 1.0000000 3.636534e-26 1.818267e-26 [254,] 1.0000000 7.661660e-26 3.830830e-26 [255,] 1.0000000 1.608514e-25 8.042568e-26 [256,] 1.0000000 3.364968e-25 1.682484e-25 [257,] 1.0000000 7.448833e-25 3.724416e-25 [258,] 1.0000000 1.618162e-24 8.090809e-25 [259,] 1.0000000 3.347615e-24 1.673808e-24 [260,] 1.0000000 6.900169e-24 3.450085e-24 [261,] 1.0000000 1.417041e-23 7.085207e-24 [262,] 1.0000000 2.739323e-23 1.369661e-23 [263,] 1.0000000 5.583270e-23 2.791635e-23 [264,] 1.0000000 1.031040e-22 5.155202e-23 [265,] 1.0000000 2.089033e-22 1.044517e-22 [266,] 1.0000000 4.216589e-22 2.108295e-22 [267,] 1.0000000 8.478355e-22 4.239178e-22 [268,] 1.0000000 1.563417e-21 7.817083e-22 [269,] 1.0000000 3.124630e-21 1.562315e-21 [270,] 1.0000000 6.540469e-21 3.270234e-21 [271,] 1.0000000 1.179464e-20 5.897319e-21 [272,] 1.0000000 2.333991e-20 1.166995e-20 [273,] 1.0000000 4.308600e-20 2.154300e-20 [274,] 1.0000000 8.777971e-20 4.388986e-20 [275,] 1.0000000 1.718004e-19 8.590018e-20 [276,] 1.0000000 3.348860e-19 1.674430e-19 [277,] 1.0000000 4.430779e-19 2.215390e-19 [278,] 1.0000000 8.610882e-19 4.305441e-19 [279,] 1.0000000 1.693535e-18 8.467677e-19 [280,] 1.0000000 3.263284e-18 1.631642e-18 [281,] 1.0000000 6.407520e-18 3.203760e-18 [282,] 1.0000000 1.254245e-17 6.271223e-18 [283,] 1.0000000 2.385953e-17 1.192977e-17 [284,] 1.0000000 4.330813e-17 2.165407e-17 [285,] 1.0000000 8.180411e-17 4.090206e-17 [286,] 1.0000000 1.601745e-16 8.008725e-17 [287,] 1.0000000 3.118189e-16 1.559094e-16 [288,] 1.0000000 5.811872e-16 2.905936e-16 [289,] 1.0000000 1.112268e-15 5.561342e-16 [290,] 1.0000000 2.140147e-15 1.070074e-15 [291,] 1.0000000 3.938176e-15 1.969088e-15 [292,] 1.0000000 7.214428e-15 3.607214e-15 [293,] 1.0000000 1.315682e-14 6.578411e-15 [294,] 1.0000000 2.483252e-14 1.241626e-14 [295,] 1.0000000 4.639858e-14 2.319929e-14 [296,] 1.0000000 8.420475e-14 4.210238e-14 [297,] 1.0000000 1.524302e-13 7.621512e-14 [298,] 1.0000000 2.716597e-13 1.358299e-13 [299,] 1.0000000 4.952269e-13 2.476135e-13 [300,] 1.0000000 8.924292e-13 4.462146e-13 [301,] 1.0000000 1.569698e-12 7.848492e-13 [302,] 1.0000000 2.698681e-12 1.349341e-12 [303,] 1.0000000 4.319144e-12 2.159572e-12 [304,] 1.0000000 7.514713e-12 3.757356e-12 [305,] 1.0000000 1.301089e-11 6.505444e-12 [306,] 1.0000000 2.292172e-11 1.146086e-11 [307,] 1.0000000 3.930495e-11 1.965247e-11 [308,] 1.0000000 6.858879e-11 3.429440e-11 [309,] 1.0000000 1.204359e-10 6.021794e-11 [310,] 1.0000000 2.034082e-10 1.017041e-10 [311,] 1.0000000 3.428875e-10 1.714437e-10 [312,] 1.0000000 5.454366e-10 2.727183e-10 [313,] 1.0000000 9.367951e-10 4.683976e-10 [314,] 1.0000000 1.553456e-09 7.767279e-10 [315,] 1.0000000 2.562618e-09 1.281309e-09 [316,] 1.0000000 4.205170e-09 2.102585e-09 [317,] 1.0000000 6.864068e-09 3.432034e-09 [318,] 1.0000000 1.114453e-08 5.572264e-09 [319,] 1.0000000 1.799730e-08 8.998650e-09 [320,] 1.0000000 2.912784e-08 1.456392e-08 [321,] 1.0000000 4.652786e-08 2.326393e-08 [322,] 1.0000000 7.391512e-08 3.695756e-08 [323,] 0.9999999 1.167755e-07 5.838776e-08 [324,] 0.9999999 1.834643e-07 9.173215e-08 [325,] 0.9999999 2.866257e-07 1.433129e-07 [326,] 0.9999998 4.452718e-07 2.226359e-07 [327,] 0.9999997 6.878033e-07 3.439016e-07 [328,] 0.9999995 1.081409e-06 5.407047e-07 [329,] 0.9999992 1.663716e-06 8.318581e-07 [330,] 0.9999987 2.567589e-06 1.283794e-06 [331,] 0.9999981 3.866945e-06 1.933472e-06 [332,] 0.9999971 5.806983e-06 2.903492e-06 [333,] 0.9999956 8.718161e-06 4.359080e-06 [334,] 0.9999934 1.328884e-05 6.644420e-06 [335,] 0.9999899 2.017675e-05 1.008837e-05 [336,] 0.9999852 2.958701e-05 1.479350e-05 [337,] 0.9999784 4.312090e-05 2.156045e-05 [338,] 0.9999688 6.245861e-05 3.122931e-05 [339,] 0.9999538 9.231621e-05 4.615810e-05 [340,] 0.9999340 1.320790e-04 6.603949e-05 [341,] 0.9999044 1.912683e-04 9.563415e-05 [342,] 0.9998649 2.701715e-04 1.350858e-04 [343,] 0.9998104 3.791864e-04 1.895932e-04 [344,] 0.9997326 5.347336e-04 2.673668e-04 [345,] 0.9996243 7.513796e-04 3.756898e-04 [346,] 0.9994700 1.060030e-03 5.300149e-04 [347,] 0.9992571 1.485813e-03 7.429067e-04 [348,] 0.9989944 2.011191e-03 1.005595e-03 [349,] 0.9986437 2.712575e-03 1.356288e-03 [350,] 0.9981830 3.633996e-03 1.816998e-03 [351,] 0.9975754 4.849102e-03 2.424551e-03 [352,] 0.9967959 6.408232e-03 3.204116e-03 [353,] 0.9957946 8.410789e-03 4.205395e-03 [354,] 0.9945184 1.096325e-02 5.481627e-03 [355,] 0.9927849 1.443016e-02 7.215080e-03 [356,] 0.9907272 1.854568e-02 9.272838e-03 [357,] 0.9881658 2.366836e-02 1.183418e-02 [358,] 0.9849628 3.007431e-02 1.503716e-02 [359,] 0.9806262 3.874767e-02 1.937383e-02 [360,] 0.9757661 4.846790e-02 2.423395e-02 [361,] 0.9699038 6.019245e-02 3.009622e-02 [362,] 0.9628919 7.421618e-02 3.710809e-02 [363,] 0.9545760 9.084796e-02 4.542398e-02 [364,] 0.9438182 1.123636e-01 5.618178e-02 [365,] 0.9322370 1.355261e-01 6.776304e-02 [366,] 0.9176026 1.647947e-01 8.239736e-02 [367,] 0.9043554 1.912892e-01 9.564460e-02 [368,] 0.8872373 2.255255e-01 1.127627e-01 [369,] 0.8678234 2.643532e-01 1.321766e-01 [370,] 0.8461841 3.076319e-01 1.538159e-01 [371,] 0.8189067 3.621866e-01 1.810933e-01 [372,] 0.7900137 4.199726e-01 2.099863e-01 [373,] 0.7615561 4.768878e-01 2.384439e-01 [374,] 0.7303780 5.392439e-01 2.696220e-01 [375,] 0.6972588 6.054823e-01 3.027412e-01 [376,] 0.6577735 6.844529e-01 3.422265e-01 [377,] 0.6208907 7.582186e-01 3.791093e-01 [378,] 0.5824354 8.351292e-01 4.175646e-01 [379,] 0.5369887 9.260225e-01 4.630113e-01 [380,] 0.4906021 9.812042e-01 5.093979e-01 [381,] 0.4478182 8.956365e-01 5.521818e-01 [382,] 0.4052253 8.104506e-01 5.947747e-01 [383,] 0.3637518 7.275035e-01 6.362482e-01 [384,] 0.3212151 6.424303e-01 6.787849e-01 [385,] 0.2808053 5.616107e-01 7.191947e-01 [386,] 0.2433549 4.867098e-01 7.566451e-01 [387,] 0.2080829 4.161658e-01 7.919171e-01 [388,] 0.1831176 3.662352e-01 8.168824e-01 [389,] 0.1847397 3.694794e-01 8.152603e-01 [390,] 0.1542966 3.085932e-01 8.457034e-01 [391,] 0.1263188 2.526376e-01 8.736812e-01 [392,] 0.9998227 3.546202e-04 1.773101e-04 [393,] 0.9996823 6.353814e-04 3.176907e-04 [394,] 0.9997480 5.039037e-04 2.519518e-04 [395,] 0.9996026 7.948014e-04 3.974007e-04 [396,] 0.9992983 1.403340e-03 7.016698e-04 [397,] 0.9987702 2.459577e-03 1.229788e-03 [398,] 0.9978372 4.325645e-03 2.162822e-03 [399,] 0.9963361 7.327835e-03 3.663918e-03 [400,] 0.9940502 1.189964e-02 5.949820e-03 [401,] 0.9902170 1.956593e-02 9.782967e-03 [402,] 0.9841478 3.170447e-02 1.585224e-02 [403,] 0.9748771 5.024573e-02 2.512287e-02 [404,] 0.9636439 7.271222e-02 3.635611e-02 [405,] 0.9738500 5.229995e-02 2.614998e-02 [406,] 0.9745168 5.096648e-02 2.548324e-02 [407,] 0.9585277 8.294470e-02 4.147235e-02 [408,] 0.9351713 1.296575e-01 6.482874e-02 [409,] 0.9050946 1.898109e-01 9.490545e-02 [410,] 0.8721892 2.556215e-01 1.278108e-01 [411,] 0.8136384 3.727232e-01 1.863616e-01 [412,] 0.7373297 5.253405e-01 2.626703e-01 [413,] 0.6444105 7.111789e-01 3.555895e-01 [414,] 0.5614008 8.771984e-01 4.385992e-01 [415,] 0.4995467 9.990934e-01 5.004533e-01 [416,] 0.5616004 8.767993e-01 4.383996e-01 [417,] 0.6038777 7.922447e-01 3.961223e-01 [418,] 0.4811527 9.623054e-01 5.188473e-01 > postscript(file="/var/www/html/rcomp/tmp/1plp11292936115.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/2plp11292936115.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/3plp11292936115.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/4hd6m1292936115.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/5hd6m1292936115.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 -130.62456449 207.79334272 -112.26455656 29.18643365 3.63482457 6 7 8 9 10 -94.81885653 18.57565718 -30.24135586 21.70176209 19.36167303 11 12 13 14 15 -86.59101631 38.24812170 286.32422167 -24.06338721 129.07075690 16 17 18 19 20 -15.18946340 0.19059469 125.09730397 32.86298850 162.88926644 21 22 23 24 25 -17.13496710 -64.60510293 -3.57357844 186.71670527 20.43441848 26 27 28 29 30 -38.62869905 80.88433341 11.20338030 -22.91166378 56.86844426 31 32 33 34 35 -15.87633601 -14.24201358 -18.79068589 -6.85565114 -69.29690936 36 37 38 39 40 98.98114436 -0.60087043 23.13146189 26.23152884 96.46713205 41 42 43 44 45 134.44917314 -30.94569690 43.71166625 28.72929285 0.64519469 46 47 48 49 50 -15.62559042 -20.47347369 81.63767597 5.53095440 15.99101616 51 52 53 54 55 47.65652634 83.45607185 18.39970423 22.29965975 59.75602632 56 57 58 59 60 50.71422965 -29.15787504 0.01050758 23.06755342 20.22741122 61 62 63 64 65 -7.22984443 14.95436832 -12.09344580 22.19916826 27.78031977 66 67 68 69 70 167.25442132 -5.48263787 11.12991749 37.01572608 42.13414125 71 72 73 74 75 30.67218546 3.88908437 5.36947546 -10.82257940 -2.38493941 76 77 78 79 80 -3.77907105 -1.39149175 -7.23762732 31.04721357 -15.85966015 81 82 83 84 85 8.71432865 22.19956634 -1.61548123 2.84655442 -10.61434171 86 87 88 89 90 26.16837306 60.52327270 -12.56087242 12.26736841 -12.16550127 91 92 93 94 95 4.49075950 6.92813611 40.98126276 -30.37710520 25.06648983 96 97 98 99 100 1.28995575 -18.62446858 -0.09282251 -22.78143049 -13.35236550 101 102 103 104 105 8.81145997 -3.52225122 15.97552526 -19.06330757 4.20756189 106 107 108 109 110 10.96884534 23.48641421 -18.44702831 -2.15469877 -3.63686258 111 112 113 114 115 1.19137508 9.20452636 1.42510254 -13.84374069 -14.17365641 116 117 118 119 120 15.18462428 -12.44437004 -15.10372818 -14.31341775 3.85762648 121 122 123 124 125 16.58075841 -6.90138075 83.71807617 -9.08693641 -9.87813040 126 127 128 129 130 -7.15741302 -13.01004748 -3.46000365 8.51128591 -13.01004748 131 132 133 134 135 -6.13429458 4.58352942 -12.01004748 6.79631498 -12.03816439 136 137 138 139 140 -10.47138135 -3.67593098 -12.02465886 -13.29954097 -13.01004748 141 142 143 144 145 -6.69696279 -13.01004748 -5.64564122 -9.22829329 -1.04779182 146 147 148 149 150 -4.55031080 -13.01004748 -6.07026144 -12.49579126 -13.14971888 151 152 153 154 155 -7.61018199 24.76624166 -8.01004748 -7.42370042 -11.14984283 156 157 158 159 160 -13.01004748 -13.01004748 -11.64122750 -10.37649821 -13.01004748 161 162 163 164 165 -13.01004748 -13.01004748 -13.01004748 -8.17523433 8.52238728 166 167 168 169 170 -12.22582453 -13.01004748 3.29305043 19.18697244 -13.01004748 171 172 173 174 175 -13.01004748 -10.30637057 -13.01004748 -5.13302816 -13.01004748 176 177 178 179 180 -5.66628636 98.11325251 -13.37070166 -13.01004748 -13.01004748 181 182 183 184 185 -1.83179678 -7.16339185 1.26927222 -0.86053356 -8.50093012 186 187 188 189 190 -5.33069669 -2.61405333 -7.96708434 -9.71788902 -13.01004748 191 192 193 194 195 -7.13130717 -13.01004748 -6.87961816 -13.01004748 -10.18633601 196 197 198 199 200 -13.01004748 39.07214005 -13.01004748 -13.01004748 0.47291583 201 202 203 204 205 -13.01004748 -13.01004748 -13.01004748 0.40646372 -9.26918373 206 207 208 209 210 -8.27238595 -12.43257522 -13.01004748 9.61379766 -11.12001887 211 212 213 214 215 -13.01004748 -13.01004748 -2.91061907 -5.04913646 -10.76776363 216 217 218 219 220 -11.01004748 -7.11145167 -7.94132213 -11.01004748 6.87557740 221 222 223 224 225 -10.99026770 -13.01004748 -12.01923585 4.80837403 -12.12773964 226 227 228 229 230 7.52053233 -4.28663498 -13.01004748 -12.84459258 -13.01004748 231 232 233 234 235 -7.89455090 -7.66017679 -13.01004748 -1.74921133 -13.01004748 236 237 238 239 240 -13.02656797 -13.01004748 -7.19268622 -13.01004748 -13.01004748 241 242 243 244 245 -13.01004748 -10.54977902 -13.01004748 -5.93745771 -24.90757444 246 247 248 249 250 -13.01004748 -13.01004748 1.25123907 -13.01004748 -13.01004748 251 252 253 254 255 -3.01918929 -13.01004748 -13.40966807 19.59623438 29.30913776 256 257 258 259 260 -13.01004748 -6.38979561 14.34355929 -13.01004748 -13.01004748 261 262 263 264 265 -13.01004748 -13.01004748 -4.21946369 -4.39771772 -13.01004748 266 267 268 269 270 -13.01004748 -13.01004748 -9.77221151 -13.01004748 10.08340022 271 272 273 274 275 -13.01004748 -13.01004748 -13.01004748 5.47252280 -13.01004748 276 277 278 279 280 -2.81884947 7.55228145 -13.01004748 8.65788379 -7.34654199 281 282 283 284 285 -13.01004748 -13.01004748 12.68482492 -13.01004748 -8.06128403 286 287 288 289 290 -13.01004748 -10.96261291 -10.48807205 -13.01004748 0.59078953 291 292 293 294 295 -13.01004748 -2.43904337 -1.49388020 -13.01004748 -2.75902125 296 297 298 299 300 -6.03662514 -13.01004748 -13.01004748 -13.01004748 -4.01004748 301 302 303 304 305 -9.12842915 -12.71831761 -12.01004748 -13.01004748 -1.54679346 306 307 308 309 310 -1.02309748 -13.01004748 4.26111624 3.54315760 -13.01004748 311 312 313 314 315 -13.01004748 11.00010465 -13.01004748 -0.09368542 -5.03840529 316 317 318 319 320 -11.22820302 -12.99930635 5.47699173 -8.01447432 -13.01004748 321 322 323 324 325 -13.01004748 -13.01004748 -13.01004748 -13.01004748 -13.01004748 326 327 328 329 330 -12.06132417 -13.01004748 -13.01004748 -13.01004748 -13.01004748 331 332 333 334 335 -13.01004748 -13.01004748 -13.01004748 -9.36600274 -12.01648475 336 337 338 339 340 -10.51453845 3.26667748 -13.01004748 -12.14374943 -3.44259767 341 342 343 344 345 -2.57123733 -13.01004748 -13.01004748 -13.01004748 -2.85074308 346 347 348 349 350 -13.01004748 -10.01653981 -13.01004748 -13.01004748 -10.72879238 351 352 353 354 355 -11.11509622 -6.35852982 -5.77851683 -12.92176715 -13.01004748 356 357 358 359 360 -13.01004748 -4.45943392 -13.01004748 -13.01004748 -13.01004748 361 362 363 364 365 -9.01560747 -13.01004748 -13.01004748 6.48759728 -4.66015046 366 367 368 369 370 -13.01004748 -13.01004748 -13.01004748 -13.01004748 -7.05796557 371 372 373 374 375 -13.01004748 -8.10779640 -10.89878890 -11.56783921 -13.01004748 376 377 378 379 380 -13.01004748 -1.24126971 -8.99011503 -13.69479069 -13.01004748 381 382 383 384 385 -13.01004748 -0.29870166 3.02316739 -13.01004748 -1.84328535 386 387 388 389 390 4.01618617 -10.63350187 -10.86463146 -10.26417164 -7.15312399 391 392 393 394 395 -5.46879276 -8.42799497 -7.05170825 -12.89307885 -24.71203709 396 397 398 399 400 -6.61715727 -7.52414074 203.84608517 -8.85458702 27.88477530 401 402 403 404 405 -16.95198399 -9.82122245 -9.26138167 -0.47417404 -8.87467144 406 407 408 409 410 -10.02281589 -6.57189259 -5.16203979 -5.26905751 -13.88120415 411 412 413 414 415 35.90812542 -27.56335233 -1.99717252 -2.48676517 -12.35691311 416 417 418 419 420 -16.25330300 -1.91803219 1.68941382 3.85325725 -14.98176598 421 422 423 424 425 -17.62138299 46.22141168 23.36785417 -16.45936332 -1.86649331 426 427 428 429 430 45.05644409 -3.98314145 39.85181744 -28.18161075 -4.67298230 431 -17.36361924 > postscript(file="/var/www/html/rcomp/tmp/6sm571292936115.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 -130.62456449 NA 1 207.79334272 -130.62456449 2 -112.26455656 207.79334272 3 29.18643365 -112.26455656 4 3.63482457 29.18643365 5 -94.81885653 3.63482457 6 18.57565718 -94.81885653 7 -30.24135586 18.57565718 8 21.70176209 -30.24135586 9 19.36167303 21.70176209 10 -86.59101631 19.36167303 11 38.24812170 -86.59101631 12 286.32422167 38.24812170 13 -24.06338721 286.32422167 14 129.07075690 -24.06338721 15 -15.18946340 129.07075690 16 0.19059469 -15.18946340 17 125.09730397 0.19059469 18 32.86298850 125.09730397 19 162.88926644 32.86298850 20 -17.13496710 162.88926644 21 -64.60510293 -17.13496710 22 -3.57357844 -64.60510293 23 186.71670527 -3.57357844 24 20.43441848 186.71670527 25 -38.62869905 20.43441848 26 80.88433341 -38.62869905 27 11.20338030 80.88433341 28 -22.91166378 11.20338030 29 56.86844426 -22.91166378 30 -15.87633601 56.86844426 31 -14.24201358 -15.87633601 32 -18.79068589 -14.24201358 33 -6.85565114 -18.79068589 34 -69.29690936 -6.85565114 35 98.98114436 -69.29690936 36 -0.60087043 98.98114436 37 23.13146189 -0.60087043 38 26.23152884 23.13146189 39 96.46713205 26.23152884 40 134.44917314 96.46713205 41 -30.94569690 134.44917314 42 43.71166625 -30.94569690 43 28.72929285 43.71166625 44 0.64519469 28.72929285 45 -15.62559042 0.64519469 46 -20.47347369 -15.62559042 47 81.63767597 -20.47347369 48 5.53095440 81.63767597 49 15.99101616 5.53095440 50 47.65652634 15.99101616 51 83.45607185 47.65652634 52 18.39970423 83.45607185 53 22.29965975 18.39970423 54 59.75602632 22.29965975 55 50.71422965 59.75602632 56 -29.15787504 50.71422965 57 0.01050758 -29.15787504 58 23.06755342 0.01050758 59 20.22741122 23.06755342 60 -7.22984443 20.22741122 61 14.95436832 -7.22984443 62 -12.09344580 14.95436832 63 22.19916826 -12.09344580 64 27.78031977 22.19916826 65 167.25442132 27.78031977 66 -5.48263787 167.25442132 67 11.12991749 -5.48263787 68 37.01572608 11.12991749 69 42.13414125 37.01572608 70 30.67218546 42.13414125 71 3.88908437 30.67218546 72 5.36947546 3.88908437 73 -10.82257940 5.36947546 74 -2.38493941 -10.82257940 75 -3.77907105 -2.38493941 76 -1.39149175 -3.77907105 77 -7.23762732 -1.39149175 78 31.04721357 -7.23762732 79 -15.85966015 31.04721357 80 8.71432865 -15.85966015 81 22.19956634 8.71432865 82 -1.61548123 22.19956634 83 2.84655442 -1.61548123 84 -10.61434171 2.84655442 85 26.16837306 -10.61434171 86 60.52327270 26.16837306 87 -12.56087242 60.52327270 88 12.26736841 -12.56087242 89 -12.16550127 12.26736841 90 4.49075950 -12.16550127 91 6.92813611 4.49075950 92 40.98126276 6.92813611 93 -30.37710520 40.98126276 94 25.06648983 -30.37710520 95 1.28995575 25.06648983 96 -18.62446858 1.28995575 97 -0.09282251 -18.62446858 98 -22.78143049 -0.09282251 99 -13.35236550 -22.78143049 100 8.81145997 -13.35236550 101 -3.52225122 8.81145997 102 15.97552526 -3.52225122 103 -19.06330757 15.97552526 104 4.20756189 -19.06330757 105 10.96884534 4.20756189 106 23.48641421 10.96884534 107 -18.44702831 23.48641421 108 -2.15469877 -18.44702831 109 -3.63686258 -2.15469877 110 1.19137508 -3.63686258 111 9.20452636 1.19137508 112 1.42510254 9.20452636 113 -13.84374069 1.42510254 114 -14.17365641 -13.84374069 115 15.18462428 -14.17365641 116 -12.44437004 15.18462428 117 -15.10372818 -12.44437004 118 -14.31341775 -15.10372818 119 3.85762648 -14.31341775 120 16.58075841 3.85762648 121 -6.90138075 16.58075841 122 83.71807617 -6.90138075 123 -9.08693641 83.71807617 124 -9.87813040 -9.08693641 125 -7.15741302 -9.87813040 126 -13.01004748 -7.15741302 127 -3.46000365 -13.01004748 128 8.51128591 -3.46000365 129 -13.01004748 8.51128591 130 -6.13429458 -13.01004748 131 4.58352942 -6.13429458 132 -12.01004748 4.58352942 133 6.79631498 -12.01004748 134 -12.03816439 6.79631498 135 -10.47138135 -12.03816439 136 -3.67593098 -10.47138135 137 -12.02465886 -3.67593098 138 -13.29954097 -12.02465886 139 -13.01004748 -13.29954097 140 -6.69696279 -13.01004748 141 -13.01004748 -6.69696279 142 -5.64564122 -13.01004748 143 -9.22829329 -5.64564122 144 -1.04779182 -9.22829329 145 -4.55031080 -1.04779182 146 -13.01004748 -4.55031080 147 -6.07026144 -13.01004748 148 -12.49579126 -6.07026144 149 -13.14971888 -12.49579126 150 -7.61018199 -13.14971888 151 24.76624166 -7.61018199 152 -8.01004748 24.76624166 153 -7.42370042 -8.01004748 154 -11.14984283 -7.42370042 155 -13.01004748 -11.14984283 156 -13.01004748 -13.01004748 157 -11.64122750 -13.01004748 158 -10.37649821 -11.64122750 159 -13.01004748 -10.37649821 160 -13.01004748 -13.01004748 161 -13.01004748 -13.01004748 162 -13.01004748 -13.01004748 163 -8.17523433 -13.01004748 164 8.52238728 -8.17523433 165 -12.22582453 8.52238728 166 -13.01004748 -12.22582453 167 3.29305043 -13.01004748 168 19.18697244 3.29305043 169 -13.01004748 19.18697244 170 -13.01004748 -13.01004748 171 -10.30637057 -13.01004748 172 -13.01004748 -10.30637057 173 -5.13302816 -13.01004748 174 -13.01004748 -5.13302816 175 -5.66628636 -13.01004748 176 98.11325251 -5.66628636 177 -13.37070166 98.11325251 178 -13.01004748 -13.37070166 179 -13.01004748 -13.01004748 180 -1.83179678 -13.01004748 181 -7.16339185 -1.83179678 182 1.26927222 -7.16339185 183 -0.86053356 1.26927222 184 -8.50093012 -0.86053356 185 -5.33069669 -8.50093012 186 -2.61405333 -5.33069669 187 -7.96708434 -2.61405333 188 -9.71788902 -7.96708434 189 -13.01004748 -9.71788902 190 -7.13130717 -13.01004748 191 -13.01004748 -7.13130717 192 -6.87961816 -13.01004748 193 -13.01004748 -6.87961816 194 -10.18633601 -13.01004748 195 -13.01004748 -10.18633601 196 39.07214005 -13.01004748 197 -13.01004748 39.07214005 198 -13.01004748 -13.01004748 199 0.47291583 -13.01004748 200 -13.01004748 0.47291583 201 -13.01004748 -13.01004748 202 -13.01004748 -13.01004748 203 0.40646372 -13.01004748 204 -9.26918373 0.40646372 205 -8.27238595 -9.26918373 206 -12.43257522 -8.27238595 207 -13.01004748 -12.43257522 208 9.61379766 -13.01004748 209 -11.12001887 9.61379766 210 -13.01004748 -11.12001887 211 -13.01004748 -13.01004748 212 -2.91061907 -13.01004748 213 -5.04913646 -2.91061907 214 -10.76776363 -5.04913646 215 -11.01004748 -10.76776363 216 -7.11145167 -11.01004748 217 -7.94132213 -7.11145167 218 -11.01004748 -7.94132213 219 6.87557740 -11.01004748 220 -10.99026770 6.87557740 221 -13.01004748 -10.99026770 222 -12.01923585 -13.01004748 223 4.80837403 -12.01923585 224 -12.12773964 4.80837403 225 7.52053233 -12.12773964 226 -4.28663498 7.52053233 227 -13.01004748 -4.28663498 228 -12.84459258 -13.01004748 229 -13.01004748 -12.84459258 230 -7.89455090 -13.01004748 231 -7.66017679 -7.89455090 232 -13.01004748 -7.66017679 233 -1.74921133 -13.01004748 234 -13.01004748 -1.74921133 235 -13.02656797 -13.01004748 236 -13.01004748 -13.02656797 237 -7.19268622 -13.01004748 238 -13.01004748 -7.19268622 239 -13.01004748 -13.01004748 240 -13.01004748 -13.01004748 241 -10.54977902 -13.01004748 242 -13.01004748 -10.54977902 243 -5.93745771 -13.01004748 244 -24.90757444 -5.93745771 245 -13.01004748 -24.90757444 246 -13.01004748 -13.01004748 247 1.25123907 -13.01004748 248 -13.01004748 1.25123907 249 -13.01004748 -13.01004748 250 -3.01918929 -13.01004748 251 -13.01004748 -3.01918929 252 -13.40966807 -13.01004748 253 19.59623438 -13.40966807 254 29.30913776 19.59623438 255 -13.01004748 29.30913776 256 -6.38979561 -13.01004748 257 14.34355929 -6.38979561 258 -13.01004748 14.34355929 259 -13.01004748 -13.01004748 260 -13.01004748 -13.01004748 261 -13.01004748 -13.01004748 262 -4.21946369 -13.01004748 263 -4.39771772 -4.21946369 264 -13.01004748 -4.39771772 265 -13.01004748 -13.01004748 266 -13.01004748 -13.01004748 267 -9.77221151 -13.01004748 268 -13.01004748 -9.77221151 269 10.08340022 -13.01004748 270 -13.01004748 10.08340022 271 -13.01004748 -13.01004748 272 -13.01004748 -13.01004748 273 5.47252280 -13.01004748 274 -13.01004748 5.47252280 275 -2.81884947 -13.01004748 276 7.55228145 -2.81884947 277 -13.01004748 7.55228145 278 8.65788379 -13.01004748 279 -7.34654199 8.65788379 280 -13.01004748 -7.34654199 281 -13.01004748 -13.01004748 282 12.68482492 -13.01004748 283 -13.01004748 12.68482492 284 -8.06128403 -13.01004748 285 -13.01004748 -8.06128403 286 -10.96261291 -13.01004748 287 -10.48807205 -10.96261291 288 -13.01004748 -10.48807205 289 0.59078953 -13.01004748 290 -13.01004748 0.59078953 291 -2.43904337 -13.01004748 292 -1.49388020 -2.43904337 293 -13.01004748 -1.49388020 294 -2.75902125 -13.01004748 295 -6.03662514 -2.75902125 296 -13.01004748 -6.03662514 297 -13.01004748 -13.01004748 298 -13.01004748 -13.01004748 299 -4.01004748 -13.01004748 300 -9.12842915 -4.01004748 301 -12.71831761 -9.12842915 302 -12.01004748 -12.71831761 303 -13.01004748 -12.01004748 304 -1.54679346 -13.01004748 305 -1.02309748 -1.54679346 306 -13.01004748 -1.02309748 307 4.26111624 -13.01004748 308 3.54315760 4.26111624 309 -13.01004748 3.54315760 310 -13.01004748 -13.01004748 311 11.00010465 -13.01004748 312 -13.01004748 11.00010465 313 -0.09368542 -13.01004748 314 -5.03840529 -0.09368542 315 -11.22820302 -5.03840529 316 -12.99930635 -11.22820302 317 5.47699173 -12.99930635 318 -8.01447432 5.47699173 319 -13.01004748 -8.01447432 320 -13.01004748 -13.01004748 321 -13.01004748 -13.01004748 322 -13.01004748 -13.01004748 323 -13.01004748 -13.01004748 324 -13.01004748 -13.01004748 325 -12.06132417 -13.01004748 326 -13.01004748 -12.06132417 327 -13.01004748 -13.01004748 328 -13.01004748 -13.01004748 329 -13.01004748 -13.01004748 330 -13.01004748 -13.01004748 331 -13.01004748 -13.01004748 332 -13.01004748 -13.01004748 333 -9.36600274 -13.01004748 334 -12.01648475 -9.36600274 335 -10.51453845 -12.01648475 336 3.26667748 -10.51453845 337 -13.01004748 3.26667748 338 -12.14374943 -13.01004748 339 -3.44259767 -12.14374943 340 -2.57123733 -3.44259767 341 -13.01004748 -2.57123733 342 -13.01004748 -13.01004748 343 -13.01004748 -13.01004748 344 -2.85074308 -13.01004748 345 -13.01004748 -2.85074308 346 -10.01653981 -13.01004748 347 -13.01004748 -10.01653981 348 -13.01004748 -13.01004748 349 -10.72879238 -13.01004748 350 -11.11509622 -10.72879238 351 -6.35852982 -11.11509622 352 -5.77851683 -6.35852982 353 -12.92176715 -5.77851683 354 -13.01004748 -12.92176715 355 -13.01004748 -13.01004748 356 -4.45943392 -13.01004748 357 -13.01004748 -4.45943392 358 -13.01004748 -13.01004748 359 -13.01004748 -13.01004748 360 -9.01560747 -13.01004748 361 -13.01004748 -9.01560747 362 -13.01004748 -13.01004748 363 6.48759728 -13.01004748 364 -4.66015046 6.48759728 365 -13.01004748 -4.66015046 366 -13.01004748 -13.01004748 367 -13.01004748 -13.01004748 368 -13.01004748 -13.01004748 369 -7.05796557 -13.01004748 370 -13.01004748 -7.05796557 371 -8.10779640 -13.01004748 372 -10.89878890 -8.10779640 373 -11.56783921 -10.89878890 374 -13.01004748 -11.56783921 375 -13.01004748 -13.01004748 376 -1.24126971 -13.01004748 377 -8.99011503 -1.24126971 378 -13.69479069 -8.99011503 379 -13.01004748 -13.69479069 380 -13.01004748 -13.01004748 381 -0.29870166 -13.01004748 382 3.02316739 -0.29870166 383 -13.01004748 3.02316739 384 -1.84328535 -13.01004748 385 4.01618617 -1.84328535 386 -10.63350187 4.01618617 387 -10.86463146 -10.63350187 388 -10.26417164 -10.86463146 389 -7.15312399 -10.26417164 390 -5.46879276 -7.15312399 391 -8.42799497 -5.46879276 392 -7.05170825 -8.42799497 393 -12.89307885 -7.05170825 394 -24.71203709 -12.89307885 395 -6.61715727 -24.71203709 396 -7.52414074 -6.61715727 397 203.84608517 -7.52414074 398 -8.85458702 203.84608517 399 27.88477530 -8.85458702 400 -16.95198399 27.88477530 401 -9.82122245 -16.95198399 402 -9.26138167 -9.82122245 403 -0.47417404 -9.26138167 404 -8.87467144 -0.47417404 405 -10.02281589 -8.87467144 406 -6.57189259 -10.02281589 407 -5.16203979 -6.57189259 408 -5.26905751 -5.16203979 409 -13.88120415 -5.26905751 410 35.90812542 -13.88120415 411 -27.56335233 35.90812542 412 -1.99717252 -27.56335233 413 -2.48676517 -1.99717252 414 -12.35691311 -2.48676517 415 -16.25330300 -12.35691311 416 -1.91803219 -16.25330300 417 1.68941382 -1.91803219 418 3.85325725 1.68941382 419 -14.98176598 3.85325725 420 -17.62138299 -14.98176598 421 46.22141168 -17.62138299 422 23.36785417 46.22141168 423 -16.45936332 23.36785417 424 -1.86649331 -16.45936332 425 45.05644409 -1.86649331 426 -3.98314145 45.05644409 427 39.85181744 -3.98314145 428 -28.18161075 39.85181744 429 -4.67298230 -28.18161075 430 -17.36361924 -4.67298230 431 NA -17.36361924 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 207.79334272 -130.62456449 [2,] -112.26455656 207.79334272 [3,] 29.18643365 -112.26455656 [4,] 3.63482457 29.18643365 [5,] -94.81885653 3.63482457 [6,] 18.57565718 -94.81885653 [7,] -30.24135586 18.57565718 [8,] 21.70176209 -30.24135586 [9,] 19.36167303 21.70176209 [10,] -86.59101631 19.36167303 [11,] 38.24812170 -86.59101631 [12,] 286.32422167 38.24812170 [13,] -24.06338721 286.32422167 [14,] 129.07075690 -24.06338721 [15,] -15.18946340 129.07075690 [16,] 0.19059469 -15.18946340 [17,] 125.09730397 0.19059469 [18,] 32.86298850 125.09730397 [19,] 162.88926644 32.86298850 [20,] -17.13496710 162.88926644 [21,] -64.60510293 -17.13496710 [22,] -3.57357844 -64.60510293 [23,] 186.71670527 -3.57357844 [24,] 20.43441848 186.71670527 [25,] -38.62869905 20.43441848 [26,] 80.88433341 -38.62869905 [27,] 11.20338030 80.88433341 [28,] -22.91166378 11.20338030 [29,] 56.86844426 -22.91166378 [30,] -15.87633601 56.86844426 [31,] -14.24201358 -15.87633601 [32,] -18.79068589 -14.24201358 [33,] -6.85565114 -18.79068589 [34,] -69.29690936 -6.85565114 [35,] 98.98114436 -69.29690936 [36,] -0.60087043 98.98114436 [37,] 23.13146189 -0.60087043 [38,] 26.23152884 23.13146189 [39,] 96.46713205 26.23152884 [40,] 134.44917314 96.46713205 [41,] -30.94569690 134.44917314 [42,] 43.71166625 -30.94569690 [43,] 28.72929285 43.71166625 [44,] 0.64519469 28.72929285 [45,] -15.62559042 0.64519469 [46,] -20.47347369 -15.62559042 [47,] 81.63767597 -20.47347369 [48,] 5.53095440 81.63767597 [49,] 15.99101616 5.53095440 [50,] 47.65652634 15.99101616 [51,] 83.45607185 47.65652634 [52,] 18.39970423 83.45607185 [53,] 22.29965975 18.39970423 [54,] 59.75602632 22.29965975 [55,] 50.71422965 59.75602632 [56,] -29.15787504 50.71422965 [57,] 0.01050758 -29.15787504 [58,] 23.06755342 0.01050758 [59,] 20.22741122 23.06755342 [60,] -7.22984443 20.22741122 [61,] 14.95436832 -7.22984443 [62,] -12.09344580 14.95436832 [63,] 22.19916826 -12.09344580 [64,] 27.78031977 22.19916826 [65,] 167.25442132 27.78031977 [66,] -5.48263787 167.25442132 [67,] 11.12991749 -5.48263787 [68,] 37.01572608 11.12991749 [69,] 42.13414125 37.01572608 [70,] 30.67218546 42.13414125 [71,] 3.88908437 30.67218546 [72,] 5.36947546 3.88908437 [73,] -10.82257940 5.36947546 [74,] -2.38493941 -10.82257940 [75,] -3.77907105 -2.38493941 [76,] -1.39149175 -3.77907105 [77,] -7.23762732 -1.39149175 [78,] 31.04721357 -7.23762732 [79,] -15.85966015 31.04721357 [80,] 8.71432865 -15.85966015 [81,] 22.19956634 8.71432865 [82,] -1.61548123 22.19956634 [83,] 2.84655442 -1.61548123 [84,] -10.61434171 2.84655442 [85,] 26.16837306 -10.61434171 [86,] 60.52327270 26.16837306 [87,] -12.56087242 60.52327270 [88,] 12.26736841 -12.56087242 [89,] -12.16550127 12.26736841 [90,] 4.49075950 -12.16550127 [91,] 6.92813611 4.49075950 [92,] 40.98126276 6.92813611 [93,] -30.37710520 40.98126276 [94,] 25.06648983 -30.37710520 [95,] 1.28995575 25.06648983 [96,] -18.62446858 1.28995575 [97,] -0.09282251 -18.62446858 [98,] -22.78143049 -0.09282251 [99,] -13.35236550 -22.78143049 [100,] 8.81145997 -13.35236550 [101,] -3.52225122 8.81145997 [102,] 15.97552526 -3.52225122 [103,] -19.06330757 15.97552526 [104,] 4.20756189 -19.06330757 [105,] 10.96884534 4.20756189 [106,] 23.48641421 10.96884534 [107,] -18.44702831 23.48641421 [108,] -2.15469877 -18.44702831 [109,] -3.63686258 -2.15469877 [110,] 1.19137508 -3.63686258 [111,] 9.20452636 1.19137508 [112,] 1.42510254 9.20452636 [113,] -13.84374069 1.42510254 [114,] -14.17365641 -13.84374069 [115,] 15.18462428 -14.17365641 [116,] -12.44437004 15.18462428 [117,] -15.10372818 -12.44437004 [118,] -14.31341775 -15.10372818 [119,] 3.85762648 -14.31341775 [120,] 16.58075841 3.85762648 [121,] -6.90138075 16.58075841 [122,] 83.71807617 -6.90138075 [123,] -9.08693641 83.71807617 [124,] -9.87813040 -9.08693641 [125,] -7.15741302 -9.87813040 [126,] -13.01004748 -7.15741302 [127,] -3.46000365 -13.01004748 [128,] 8.51128591 -3.46000365 [129,] -13.01004748 8.51128591 [130,] -6.13429458 -13.01004748 [131,] 4.58352942 -6.13429458 [132,] -12.01004748 4.58352942 [133,] 6.79631498 -12.01004748 [134,] -12.03816439 6.79631498 [135,] -10.47138135 -12.03816439 [136,] -3.67593098 -10.47138135 [137,] -12.02465886 -3.67593098 [138,] -13.29954097 -12.02465886 [139,] -13.01004748 -13.29954097 [140,] -6.69696279 -13.01004748 [141,] -13.01004748 -6.69696279 [142,] -5.64564122 -13.01004748 [143,] -9.22829329 -5.64564122 [144,] -1.04779182 -9.22829329 [145,] -4.55031080 -1.04779182 [146,] -13.01004748 -4.55031080 [147,] -6.07026144 -13.01004748 [148,] -12.49579126 -6.07026144 [149,] -13.14971888 -12.49579126 [150,] -7.61018199 -13.14971888 [151,] 24.76624166 -7.61018199 [152,] -8.01004748 24.76624166 [153,] -7.42370042 -8.01004748 [154,] -11.14984283 -7.42370042 [155,] -13.01004748 -11.14984283 [156,] -13.01004748 -13.01004748 [157,] -11.64122750 -13.01004748 [158,] -10.37649821 -11.64122750 [159,] -13.01004748 -10.37649821 [160,] -13.01004748 -13.01004748 [161,] -13.01004748 -13.01004748 [162,] -13.01004748 -13.01004748 [163,] -8.17523433 -13.01004748 [164,] 8.52238728 -8.17523433 [165,] -12.22582453 8.52238728 [166,] -13.01004748 -12.22582453 [167,] 3.29305043 -13.01004748 [168,] 19.18697244 3.29305043 [169,] -13.01004748 19.18697244 [170,] -13.01004748 -13.01004748 [171,] -10.30637057 -13.01004748 [172,] -13.01004748 -10.30637057 [173,] -5.13302816 -13.01004748 [174,] -13.01004748 -5.13302816 [175,] -5.66628636 -13.01004748 [176,] 98.11325251 -5.66628636 [177,] -13.37070166 98.11325251 [178,] -13.01004748 -13.37070166 [179,] -13.01004748 -13.01004748 [180,] -1.83179678 -13.01004748 [181,] -7.16339185 -1.83179678 [182,] 1.26927222 -7.16339185 [183,] -0.86053356 1.26927222 [184,] -8.50093012 -0.86053356 [185,] -5.33069669 -8.50093012 [186,] -2.61405333 -5.33069669 [187,] -7.96708434 -2.61405333 [188,] -9.71788902 -7.96708434 [189,] -13.01004748 -9.71788902 [190,] -7.13130717 -13.01004748 [191,] -13.01004748 -7.13130717 [192,] -6.87961816 -13.01004748 [193,] -13.01004748 -6.87961816 [194,] -10.18633601 -13.01004748 [195,] -13.01004748 -10.18633601 [196,] 39.07214005 -13.01004748 [197,] -13.01004748 39.07214005 [198,] -13.01004748 -13.01004748 [199,] 0.47291583 -13.01004748 [200,] -13.01004748 0.47291583 [201,] -13.01004748 -13.01004748 [202,] -13.01004748 -13.01004748 [203,] 0.40646372 -13.01004748 [204,] -9.26918373 0.40646372 [205,] -8.27238595 -9.26918373 [206,] -12.43257522 -8.27238595 [207,] -13.01004748 -12.43257522 [208,] 9.61379766 -13.01004748 [209,] -11.12001887 9.61379766 [210,] -13.01004748 -11.12001887 [211,] -13.01004748 -13.01004748 [212,] -2.91061907 -13.01004748 [213,] -5.04913646 -2.91061907 [214,] -10.76776363 -5.04913646 [215,] -11.01004748 -10.76776363 [216,] -7.11145167 -11.01004748 [217,] -7.94132213 -7.11145167 [218,] -11.01004748 -7.94132213 [219,] 6.87557740 -11.01004748 [220,] -10.99026770 6.87557740 [221,] -13.01004748 -10.99026770 [222,] -12.01923585 -13.01004748 [223,] 4.80837403 -12.01923585 [224,] -12.12773964 4.80837403 [225,] 7.52053233 -12.12773964 [226,] -4.28663498 7.52053233 [227,] -13.01004748 -4.28663498 [228,] -12.84459258 -13.01004748 [229,] -13.01004748 -12.84459258 [230,] -7.89455090 -13.01004748 [231,] -7.66017679 -7.89455090 [232,] -13.01004748 -7.66017679 [233,] -1.74921133 -13.01004748 [234,] -13.01004748 -1.74921133 [235,] -13.02656797 -13.01004748 [236,] -13.01004748 -13.02656797 [237,] -7.19268622 -13.01004748 [238,] -13.01004748 -7.19268622 [239,] -13.01004748 -13.01004748 [240,] -13.01004748 -13.01004748 [241,] -10.54977902 -13.01004748 [242,] -13.01004748 -10.54977902 [243,] -5.93745771 -13.01004748 [244,] -24.90757444 -5.93745771 [245,] -13.01004748 -24.90757444 [246,] -13.01004748 -13.01004748 [247,] 1.25123907 -13.01004748 [248,] -13.01004748 1.25123907 [249,] -13.01004748 -13.01004748 [250,] -3.01918929 -13.01004748 [251,] -13.01004748 -3.01918929 [252,] -13.40966807 -13.01004748 [253,] 19.59623438 -13.40966807 [254,] 29.30913776 19.59623438 [255,] -13.01004748 29.30913776 [256,] -6.38979561 -13.01004748 [257,] 14.34355929 -6.38979561 [258,] -13.01004748 14.34355929 [259,] -13.01004748 -13.01004748 [260,] -13.01004748 -13.01004748 [261,] -13.01004748 -13.01004748 [262,] -4.21946369 -13.01004748 [263,] -4.39771772 -4.21946369 [264,] -13.01004748 -4.39771772 [265,] -13.01004748 -13.01004748 [266,] -13.01004748 -13.01004748 [267,] -9.77221151 -13.01004748 [268,] -13.01004748 -9.77221151 [269,] 10.08340022 -13.01004748 [270,] -13.01004748 10.08340022 [271,] -13.01004748 -13.01004748 [272,] -13.01004748 -13.01004748 [273,] 5.47252280 -13.01004748 [274,] -13.01004748 5.47252280 [275,] -2.81884947 -13.01004748 [276,] 7.55228145 -2.81884947 [277,] -13.01004748 7.55228145 [278,] 8.65788379 -13.01004748 [279,] -7.34654199 8.65788379 [280,] -13.01004748 -7.34654199 [281,] -13.01004748 -13.01004748 [282,] 12.68482492 -13.01004748 [283,] -13.01004748 12.68482492 [284,] -8.06128403 -13.01004748 [285,] -13.01004748 -8.06128403 [286,] -10.96261291 -13.01004748 [287,] -10.48807205 -10.96261291 [288,] -13.01004748 -10.48807205 [289,] 0.59078953 -13.01004748 [290,] -13.01004748 0.59078953 [291,] -2.43904337 -13.01004748 [292,] -1.49388020 -2.43904337 [293,] -13.01004748 -1.49388020 [294,] -2.75902125 -13.01004748 [295,] -6.03662514 -2.75902125 [296,] -13.01004748 -6.03662514 [297,] -13.01004748 -13.01004748 [298,] -13.01004748 -13.01004748 [299,] -4.01004748 -13.01004748 [300,] -9.12842915 -4.01004748 [301,] -12.71831761 -9.12842915 [302,] -12.01004748 -12.71831761 [303,] -13.01004748 -12.01004748 [304,] -1.54679346 -13.01004748 [305,] -1.02309748 -1.54679346 [306,] -13.01004748 -1.02309748 [307,] 4.26111624 -13.01004748 [308,] 3.54315760 4.26111624 [309,] -13.01004748 3.54315760 [310,] -13.01004748 -13.01004748 [311,] 11.00010465 -13.01004748 [312,] -13.01004748 11.00010465 [313,] -0.09368542 -13.01004748 [314,] -5.03840529 -0.09368542 [315,] -11.22820302 -5.03840529 [316,] -12.99930635 -11.22820302 [317,] 5.47699173 -12.99930635 [318,] -8.01447432 5.47699173 [319,] -13.01004748 -8.01447432 [320,] -13.01004748 -13.01004748 [321,] -13.01004748 -13.01004748 [322,] -13.01004748 -13.01004748 [323,] -13.01004748 -13.01004748 [324,] -13.01004748 -13.01004748 [325,] -12.06132417 -13.01004748 [326,] -13.01004748 -12.06132417 [327,] -13.01004748 -13.01004748 [328,] -13.01004748 -13.01004748 [329,] -13.01004748 -13.01004748 [330,] -13.01004748 -13.01004748 [331,] -13.01004748 -13.01004748 [332,] -13.01004748 -13.01004748 [333,] -9.36600274 -13.01004748 [334,] -12.01648475 -9.36600274 [335,] -10.51453845 -12.01648475 [336,] 3.26667748 -10.51453845 [337,] -13.01004748 3.26667748 [338,] -12.14374943 -13.01004748 [339,] -3.44259767 -12.14374943 [340,] -2.57123733 -3.44259767 [341,] -13.01004748 -2.57123733 [342,] -13.01004748 -13.01004748 [343,] -13.01004748 -13.01004748 [344,] -2.85074308 -13.01004748 [345,] -13.01004748 -2.85074308 [346,] -10.01653981 -13.01004748 [347,] -13.01004748 -10.01653981 [348,] -13.01004748 -13.01004748 [349,] -10.72879238 -13.01004748 [350,] -11.11509622 -10.72879238 [351,] -6.35852982 -11.11509622 [352,] -5.77851683 -6.35852982 [353,] -12.92176715 -5.77851683 [354,] -13.01004748 -12.92176715 [355,] -13.01004748 -13.01004748 [356,] -4.45943392 -13.01004748 [357,] -13.01004748 -4.45943392 [358,] -13.01004748 -13.01004748 [359,] -13.01004748 -13.01004748 [360,] -9.01560747 -13.01004748 [361,] -13.01004748 -9.01560747 [362,] -13.01004748 -13.01004748 [363,] 6.48759728 -13.01004748 [364,] -4.66015046 6.48759728 [365,] -13.01004748 -4.66015046 [366,] -13.01004748 -13.01004748 [367,] -13.01004748 -13.01004748 [368,] -13.01004748 -13.01004748 [369,] -7.05796557 -13.01004748 [370,] -13.01004748 -7.05796557 [371,] -8.10779640 -13.01004748 [372,] -10.89878890 -8.10779640 [373,] -11.56783921 -10.89878890 [374,] -13.01004748 -11.56783921 [375,] -13.01004748 -13.01004748 [376,] -1.24126971 -13.01004748 [377,] -8.99011503 -1.24126971 [378,] -13.69479069 -8.99011503 [379,] -13.01004748 -13.69479069 [380,] -13.01004748 -13.01004748 [381,] -0.29870166 -13.01004748 [382,] 3.02316739 -0.29870166 [383,] -13.01004748 3.02316739 [384,] -1.84328535 -13.01004748 [385,] 4.01618617 -1.84328535 [386,] -10.63350187 4.01618617 [387,] -10.86463146 -10.63350187 [388,] -10.26417164 -10.86463146 [389,] -7.15312399 -10.26417164 [390,] -5.46879276 -7.15312399 [391,] -8.42799497 -5.46879276 [392,] -7.05170825 -8.42799497 [393,] -12.89307885 -7.05170825 [394,] -24.71203709 -12.89307885 [395,] -6.61715727 -24.71203709 [396,] -7.52414074 -6.61715727 [397,] 203.84608517 -7.52414074 [398,] -8.85458702 203.84608517 [399,] 27.88477530 -8.85458702 [400,] -16.95198399 27.88477530 [401,] -9.82122245 -16.95198399 [402,] -9.26138167 -9.82122245 [403,] -0.47417404 -9.26138167 [404,] -8.87467144 -0.47417404 [405,] -10.02281589 -8.87467144 [406,] -6.57189259 -10.02281589 [407,] -5.16203979 -6.57189259 [408,] -5.26905751 -5.16203979 [409,] -13.88120415 -5.26905751 [410,] 35.90812542 -13.88120415 [411,] -27.56335233 35.90812542 [412,] -1.99717252 -27.56335233 [413,] -2.48676517 -1.99717252 [414,] -12.35691311 -2.48676517 [415,] -16.25330300 -12.35691311 [416,] -1.91803219 -16.25330300 [417,] 1.68941382 -1.91803219 [418,] 3.85325725 1.68941382 [419,] -14.98176598 3.85325725 [420,] -17.62138299 -14.98176598 [421,] 46.22141168 -17.62138299 [422,] 23.36785417 46.22141168 [423,] -16.45936332 23.36785417 [424,] -1.86649331 -16.45936332 [425,] 45.05644409 -1.86649331 [426,] -3.98314145 45.05644409 [427,] 39.85181744 -3.98314145 [428,] -28.18161075 39.85181744 [429,] -4.67298230 -28.18161075 [430,] -17.36361924 -4.67298230 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 207.79334272 -130.62456449 2 -112.26455656 207.79334272 3 29.18643365 -112.26455656 4 3.63482457 29.18643365 5 -94.81885653 3.63482457 6 18.57565718 -94.81885653 7 -30.24135586 18.57565718 8 21.70176209 -30.24135586 9 19.36167303 21.70176209 10 -86.59101631 19.36167303 11 38.24812170 -86.59101631 12 286.32422167 38.24812170 13 -24.06338721 286.32422167 14 129.07075690 -24.06338721 15 -15.18946340 129.07075690 16 0.19059469 -15.18946340 17 125.09730397 0.19059469 18 32.86298850 125.09730397 19 162.88926644 32.86298850 20 -17.13496710 162.88926644 21 -64.60510293 -17.13496710 22 -3.57357844 -64.60510293 23 186.71670527 -3.57357844 24 20.43441848 186.71670527 25 -38.62869905 20.43441848 26 80.88433341 -38.62869905 27 11.20338030 80.88433341 28 -22.91166378 11.20338030 29 56.86844426 -22.91166378 30 -15.87633601 56.86844426 31 -14.24201358 -15.87633601 32 -18.79068589 -14.24201358 33 -6.85565114 -18.79068589 34 -69.29690936 -6.85565114 35 98.98114436 -69.29690936 36 -0.60087043 98.98114436 37 23.13146189 -0.60087043 38 26.23152884 23.13146189 39 96.46713205 26.23152884 40 134.44917314 96.46713205 41 -30.94569690 134.44917314 42 43.71166625 -30.94569690 43 28.72929285 43.71166625 44 0.64519469 28.72929285 45 -15.62559042 0.64519469 46 -20.47347369 -15.62559042 47 81.63767597 -20.47347369 48 5.53095440 81.63767597 49 15.99101616 5.53095440 50 47.65652634 15.99101616 51 83.45607185 47.65652634 52 18.39970423 83.45607185 53 22.29965975 18.39970423 54 59.75602632 22.29965975 55 50.71422965 59.75602632 56 -29.15787504 50.71422965 57 0.01050758 -29.15787504 58 23.06755342 0.01050758 59 20.22741122 23.06755342 60 -7.22984443 20.22741122 61 14.95436832 -7.22984443 62 -12.09344580 14.95436832 63 22.19916826 -12.09344580 64 27.78031977 22.19916826 65 167.25442132 27.78031977 66 -5.48263787 167.25442132 67 11.12991749 -5.48263787 68 37.01572608 11.12991749 69 42.13414125 37.01572608 70 30.67218546 42.13414125 71 3.88908437 30.67218546 72 5.36947546 3.88908437 73 -10.82257940 5.36947546 74 -2.38493941 -10.82257940 75 -3.77907105 -2.38493941 76 -1.39149175 -3.77907105 77 -7.23762732 -1.39149175 78 31.04721357 -7.23762732 79 -15.85966015 31.04721357 80 8.71432865 -15.85966015 81 22.19956634 8.71432865 82 -1.61548123 22.19956634 83 2.84655442 -1.61548123 84 -10.61434171 2.84655442 85 26.16837306 -10.61434171 86 60.52327270 26.16837306 87 -12.56087242 60.52327270 88 12.26736841 -12.56087242 89 -12.16550127 12.26736841 90 4.49075950 -12.16550127 91 6.92813611 4.49075950 92 40.98126276 6.92813611 93 -30.37710520 40.98126276 94 25.06648983 -30.37710520 95 1.28995575 25.06648983 96 -18.62446858 1.28995575 97 -0.09282251 -18.62446858 98 -22.78143049 -0.09282251 99 -13.35236550 -22.78143049 100 8.81145997 -13.35236550 101 -3.52225122 8.81145997 102 15.97552526 -3.52225122 103 -19.06330757 15.97552526 104 4.20756189 -19.06330757 105 10.96884534 4.20756189 106 23.48641421 10.96884534 107 -18.44702831 23.48641421 108 -2.15469877 -18.44702831 109 -3.63686258 -2.15469877 110 1.19137508 -3.63686258 111 9.20452636 1.19137508 112 1.42510254 9.20452636 113 -13.84374069 1.42510254 114 -14.17365641 -13.84374069 115 15.18462428 -14.17365641 116 -12.44437004 15.18462428 117 -15.10372818 -12.44437004 118 -14.31341775 -15.10372818 119 3.85762648 -14.31341775 120 16.58075841 3.85762648 121 -6.90138075 16.58075841 122 83.71807617 -6.90138075 123 -9.08693641 83.71807617 124 -9.87813040 -9.08693641 125 -7.15741302 -9.87813040 126 -13.01004748 -7.15741302 127 -3.46000365 -13.01004748 128 8.51128591 -3.46000365 129 -13.01004748 8.51128591 130 -6.13429458 -13.01004748 131 4.58352942 -6.13429458 132 -12.01004748 4.58352942 133 6.79631498 -12.01004748 134 -12.03816439 6.79631498 135 -10.47138135 -12.03816439 136 -3.67593098 -10.47138135 137 -12.02465886 -3.67593098 138 -13.29954097 -12.02465886 139 -13.01004748 -13.29954097 140 -6.69696279 -13.01004748 141 -13.01004748 -6.69696279 142 -5.64564122 -13.01004748 143 -9.22829329 -5.64564122 144 -1.04779182 -9.22829329 145 -4.55031080 -1.04779182 146 -13.01004748 -4.55031080 147 -6.07026144 -13.01004748 148 -12.49579126 -6.07026144 149 -13.14971888 -12.49579126 150 -7.61018199 -13.14971888 151 24.76624166 -7.61018199 152 -8.01004748 24.76624166 153 -7.42370042 -8.01004748 154 -11.14984283 -7.42370042 155 -13.01004748 -11.14984283 156 -13.01004748 -13.01004748 157 -11.64122750 -13.01004748 158 -10.37649821 -11.64122750 159 -13.01004748 -10.37649821 160 -13.01004748 -13.01004748 161 -13.01004748 -13.01004748 162 -13.01004748 -13.01004748 163 -8.17523433 -13.01004748 164 8.52238728 -8.17523433 165 -12.22582453 8.52238728 166 -13.01004748 -12.22582453 167 3.29305043 -13.01004748 168 19.18697244 3.29305043 169 -13.01004748 19.18697244 170 -13.01004748 -13.01004748 171 -10.30637057 -13.01004748 172 -13.01004748 -10.30637057 173 -5.13302816 -13.01004748 174 -13.01004748 -5.13302816 175 -5.66628636 -13.01004748 176 98.11325251 -5.66628636 177 -13.37070166 98.11325251 178 -13.01004748 -13.37070166 179 -13.01004748 -13.01004748 180 -1.83179678 -13.01004748 181 -7.16339185 -1.83179678 182 1.26927222 -7.16339185 183 -0.86053356 1.26927222 184 -8.50093012 -0.86053356 185 -5.33069669 -8.50093012 186 -2.61405333 -5.33069669 187 -7.96708434 -2.61405333 188 -9.71788902 -7.96708434 189 -13.01004748 -9.71788902 190 -7.13130717 -13.01004748 191 -13.01004748 -7.13130717 192 -6.87961816 -13.01004748 193 -13.01004748 -6.87961816 194 -10.18633601 -13.01004748 195 -13.01004748 -10.18633601 196 39.07214005 -13.01004748 197 -13.01004748 39.07214005 198 -13.01004748 -13.01004748 199 0.47291583 -13.01004748 200 -13.01004748 0.47291583 201 -13.01004748 -13.01004748 202 -13.01004748 -13.01004748 203 0.40646372 -13.01004748 204 -9.26918373 0.40646372 205 -8.27238595 -9.26918373 206 -12.43257522 -8.27238595 207 -13.01004748 -12.43257522 208 9.61379766 -13.01004748 209 -11.12001887 9.61379766 210 -13.01004748 -11.12001887 211 -13.01004748 -13.01004748 212 -2.91061907 -13.01004748 213 -5.04913646 -2.91061907 214 -10.76776363 -5.04913646 215 -11.01004748 -10.76776363 216 -7.11145167 -11.01004748 217 -7.94132213 -7.11145167 218 -11.01004748 -7.94132213 219 6.87557740 -11.01004748 220 -10.99026770 6.87557740 221 -13.01004748 -10.99026770 222 -12.01923585 -13.01004748 223 4.80837403 -12.01923585 224 -12.12773964 4.80837403 225 7.52053233 -12.12773964 226 -4.28663498 7.52053233 227 -13.01004748 -4.28663498 228 -12.84459258 -13.01004748 229 -13.01004748 -12.84459258 230 -7.89455090 -13.01004748 231 -7.66017679 -7.89455090 232 -13.01004748 -7.66017679 233 -1.74921133 -13.01004748 234 -13.01004748 -1.74921133 235 -13.02656797 -13.01004748 236 -13.01004748 -13.02656797 237 -7.19268622 -13.01004748 238 -13.01004748 -7.19268622 239 -13.01004748 -13.01004748 240 -13.01004748 -13.01004748 241 -10.54977902 -13.01004748 242 -13.01004748 -10.54977902 243 -5.93745771 -13.01004748 244 -24.90757444 -5.93745771 245 -13.01004748 -24.90757444 246 -13.01004748 -13.01004748 247 1.25123907 -13.01004748 248 -13.01004748 1.25123907 249 -13.01004748 -13.01004748 250 -3.01918929 -13.01004748 251 -13.01004748 -3.01918929 252 -13.40966807 -13.01004748 253 19.59623438 -13.40966807 254 29.30913776 19.59623438 255 -13.01004748 29.30913776 256 -6.38979561 -13.01004748 257 14.34355929 -6.38979561 258 -13.01004748 14.34355929 259 -13.01004748 -13.01004748 260 -13.01004748 -13.01004748 261 -13.01004748 -13.01004748 262 -4.21946369 -13.01004748 263 -4.39771772 -4.21946369 264 -13.01004748 -4.39771772 265 -13.01004748 -13.01004748 266 -13.01004748 -13.01004748 267 -9.77221151 -13.01004748 268 -13.01004748 -9.77221151 269 10.08340022 -13.01004748 270 -13.01004748 10.08340022 271 -13.01004748 -13.01004748 272 -13.01004748 -13.01004748 273 5.47252280 -13.01004748 274 -13.01004748 5.47252280 275 -2.81884947 -13.01004748 276 7.55228145 -2.81884947 277 -13.01004748 7.55228145 278 8.65788379 -13.01004748 279 -7.34654199 8.65788379 280 -13.01004748 -7.34654199 281 -13.01004748 -13.01004748 282 12.68482492 -13.01004748 283 -13.01004748 12.68482492 284 -8.06128403 -13.01004748 285 -13.01004748 -8.06128403 286 -10.96261291 -13.01004748 287 -10.48807205 -10.96261291 288 -13.01004748 -10.48807205 289 0.59078953 -13.01004748 290 -13.01004748 0.59078953 291 -2.43904337 -13.01004748 292 -1.49388020 -2.43904337 293 -13.01004748 -1.49388020 294 -2.75902125 -13.01004748 295 -6.03662514 -2.75902125 296 -13.01004748 -6.03662514 297 -13.01004748 -13.01004748 298 -13.01004748 -13.01004748 299 -4.01004748 -13.01004748 300 -9.12842915 -4.01004748 301 -12.71831761 -9.12842915 302 -12.01004748 -12.71831761 303 -13.01004748 -12.01004748 304 -1.54679346 -13.01004748 305 -1.02309748 -1.54679346 306 -13.01004748 -1.02309748 307 4.26111624 -13.01004748 308 3.54315760 4.26111624 309 -13.01004748 3.54315760 310 -13.01004748 -13.01004748 311 11.00010465 -13.01004748 312 -13.01004748 11.00010465 313 -0.09368542 -13.01004748 314 -5.03840529 -0.09368542 315 -11.22820302 -5.03840529 316 -12.99930635 -11.22820302 317 5.47699173 -12.99930635 318 -8.01447432 5.47699173 319 -13.01004748 -8.01447432 320 -13.01004748 -13.01004748 321 -13.01004748 -13.01004748 322 -13.01004748 -13.01004748 323 -13.01004748 -13.01004748 324 -13.01004748 -13.01004748 325 -12.06132417 -13.01004748 326 -13.01004748 -12.06132417 327 -13.01004748 -13.01004748 328 -13.01004748 -13.01004748 329 -13.01004748 -13.01004748 330 -13.01004748 -13.01004748 331 -13.01004748 -13.01004748 332 -13.01004748 -13.01004748 333 -9.36600274 -13.01004748 334 -12.01648475 -9.36600274 335 -10.51453845 -12.01648475 336 3.26667748 -10.51453845 337 -13.01004748 3.26667748 338 -12.14374943 -13.01004748 339 -3.44259767 -12.14374943 340 -2.57123733 -3.44259767 341 -13.01004748 -2.57123733 342 -13.01004748 -13.01004748 343 -13.01004748 -13.01004748 344 -2.85074308 -13.01004748 345 -13.01004748 -2.85074308 346 -10.01653981 -13.01004748 347 -13.01004748 -10.01653981 348 -13.01004748 -13.01004748 349 -10.72879238 -13.01004748 350 -11.11509622 -10.72879238 351 -6.35852982 -11.11509622 352 -5.77851683 -6.35852982 353 -12.92176715 -5.77851683 354 -13.01004748 -12.92176715 355 -13.01004748 -13.01004748 356 -4.45943392 -13.01004748 357 -13.01004748 -4.45943392 358 -13.01004748 -13.01004748 359 -13.01004748 -13.01004748 360 -9.01560747 -13.01004748 361 -13.01004748 -9.01560747 362 -13.01004748 -13.01004748 363 6.48759728 -13.01004748 364 -4.66015046 6.48759728 365 -13.01004748 -4.66015046 366 -13.01004748 -13.01004748 367 -13.01004748 -13.01004748 368 -13.01004748 -13.01004748 369 -7.05796557 -13.01004748 370 -13.01004748 -7.05796557 371 -8.10779640 -13.01004748 372 -10.89878890 -8.10779640 373 -11.56783921 -10.89878890 374 -13.01004748 -11.56783921 375 -13.01004748 -13.01004748 376 -1.24126971 -13.01004748 377 -8.99011503 -1.24126971 378 -13.69479069 -8.99011503 379 -13.01004748 -13.69479069 380 -13.01004748 -13.01004748 381 -0.29870166 -13.01004748 382 3.02316739 -0.29870166 383 -13.01004748 3.02316739 384 -1.84328535 -13.01004748 385 4.01618617 -1.84328535 386 -10.63350187 4.01618617 387 -10.86463146 -10.63350187 388 -10.26417164 -10.86463146 389 -7.15312399 -10.26417164 390 -5.46879276 -7.15312399 391 -8.42799497 -5.46879276 392 -7.05170825 -8.42799497 393 -12.89307885 -7.05170825 394 -24.71203709 -12.89307885 395 -6.61715727 -24.71203709 396 -7.52414074 -6.61715727 397 203.84608517 -7.52414074 398 -8.85458702 203.84608517 399 27.88477530 -8.85458702 400 -16.95198399 27.88477530 401 -9.82122245 -16.95198399 402 -9.26138167 -9.82122245 403 -0.47417404 -9.26138167 404 -8.87467144 -0.47417404 405 -10.02281589 -8.87467144 406 -6.57189259 -10.02281589 407 -5.16203979 -6.57189259 408 -5.26905751 -5.16203979 409 -13.88120415 -5.26905751 410 35.90812542 -13.88120415 411 -27.56335233 35.90812542 412 -1.99717252 -27.56335233 413 -2.48676517 -1.99717252 414 -12.35691311 -2.48676517 415 -16.25330300 -12.35691311 416 -1.91803219 -16.25330300 417 1.68941382 -1.91803219 418 3.85325725 1.68941382 419 -14.98176598 3.85325725 420 -17.62138299 -14.98176598 421 46.22141168 -17.62138299 422 23.36785417 46.22141168 423 -16.45936332 23.36785417 424 -1.86649331 -16.45936332 425 45.05644409 -1.86649331 426 -3.98314145 45.05644409 427 39.85181744 -3.98314145 428 -28.18161075 39.85181744 429 -4.67298230 -28.18161075 430 -17.36361924 -4.67298230 > 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/7lv5a1292936115.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/8lv5a1292936115.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/9lv5a1292936115.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/10w54u1292936115.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/11h5li1292936115.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/12se231292936115.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/13gfhf1292936115.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/147j4r1292936115.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/155gwr1292936115.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/169hcx1292936115.tab") + } > > try(system("convert tmp/1plp11292936115.ps tmp/1plp11292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/2plp11292936115.ps tmp/2plp11292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/3plp11292936115.ps tmp/3plp11292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/4hd6m1292936115.ps tmp/4hd6m1292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/5hd6m1292936115.ps tmp/5hd6m1292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/6sm571292936115.ps tmp/6sm571292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/7lv5a1292936115.ps tmp/7lv5a1292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/8lv5a1292936115.ps tmp/8lv5a1292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/9lv5a1292936115.ps tmp/9lv5a1292936115.png",intern=TRUE)) character(0) > try(system("convert tmp/10w54u1292936115.ps tmp/10w54u1292936115.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 11.597 2.085 26.507