R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(162556 + ,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 + ,0 + ,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('Kosten' + ,'Oders' + ,'Divid.' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(4,431),dimnames=list(c('Kosten','Oders','Divid.','Wealth'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '4' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Wealth Kosten Oders Divid. 1 6282154 162556 807 213118 2 4321023 29790 444 81767 3 4111912 87550 412 153198 4 223193 84738 428 -26007 5 1491348 54660 315 126942 6 1629616 42634 168 157214 7 1398893 40949 263 129352 8 1926517 45187 267 234817 9 983660 37704 228 60448 10 1443586 16275 129 47818 11 1073089 25830 104 245546 12 984885 12679 122 48020 13 1405225 18014 393 -1710 14 227132 43556 190 32648 15 929118 24811 280 95350 16 1071292 6575 63 151352 17 638830 7123 102 288170 18 856956 21950 265 114337 19 992426 37597 234 37884 20 444477 17821 277 122844 21 857217 12988 73 82340 22 711969 22330 67 79801 23 702380 13326 103 165548 24 358589 16189 290 116384 25 297978 7146 83 134028 26 585715 15824 56 63838 27 657954 27664 236 74996 28 209458 11920 73 31080 29 786690 8568 34 32168 30 439798 14416 139 49857 31 688779 3369 26 87161 32 574339 11819 70 106113 33 741409 6984 40 80570 34 597793 4519 42 102129 35 644190 2220 12 301670 36 377934 18562 211 102313 37 640273 10327 74 88577 38 697458 5336 80 112477 39 550608 2365 83 191778 40 207393 4069 131 79804 41 301607 8636 203 128294 42 345783 13718 56 96448 43 501749 4525 89 93811 44 379983 6869 88 117520 45 387475 4628 39 69159 46 377305 3689 25 101792 47 370837 4891 49 210568 48 430866 7489 149 136996 49 469107 4901 58 121920 50 194493 2284 41 76403 51 530670 3160 90 108094 52 518365 4150 136 134759 53 491303 7285 97 188873 54 527021 1134 63 146216 55 233773 4658 114 156608 56 405972 2384 77 61348 57 652925 3748 6 50350 58 446211 5371 47 87720 59 341340 1285 51 99489 60 387699 9327 85 87419 61 493408 5565 43 94355 62 146494 1528 32 60326 63 414462 3122 25 94670 64 364304 7561 77 82425 65 355178 2675 54 59017 66 357760 13253 251 90829 67 261216 880 15 80791 68 397144 2053 44 100423 69 374943 1424 73 131116 70 424898 4036 85 100269 71 202055 3045 49 27330 72 378525 5119 38 39039 73 310768 1431 35 106885 74 325738 554 9 79285 75 394510 1975 34 118881 76 247060 1765 20 77623 77 368078 1012 29 114768 78 236761 810 11 74015 79 312378 1280 52 69465 80 339836 666 13 117869 81 347385 1380 29 60982 82 426280 4677 66 90131 83 352850 876 33 138971 84 301881 814 15 39625 85 377516 514 15 102725 86 357312 5692 68 64239 87 458343 3642 100 90262 88 354228 540 13 103960 89 308636 2099 45 106611 90 386212 567 14 103345 91 393343 2001 36 95551 92 378509 2949 40 82903 93 452469 2253 68 63593 94 364839 6533 29 126910 95 358649 1889 43 37527 96 376641 3055 30 60247 97 429112 272 9 112995 98 330546 1414 22 70184 99 403560 2564 19 130140 100 317892 1383 9 73221 101 307528 1261 31 76114 102 235133 975 19 90534 103 299243 3366 55 108479 104 314073 576 8 113761 105 368186 1686 28 68696 106 269661 746 29 71561 107 125390 3192 48 59831 108 510834 2045 16 97890 109 321896 5702 47 101481 110 249898 1932 20 72954 111 408881 936 22 67939 112 158492 3437 33 48022 113 292154 5131 44 86111 114 289513 2397 13 74020 115 378049 1389 6 57530 116 343466 1503 35 56364 117 332743 402 8 84990 118 442882 2239 17 88590 119 214215 2234 11 77200 120 315688 837 21 61262 121 375195 10579 92 110309 122 334280 875 12 67000 123 355864 1585 112 93099 124 480382 1659 25 107577 125 353058 2647 17 62920 126 217193 3294 23 75832 127 315380 0 0 60720 128 314533 94 10 60793 129 318056 422 23 57935 130 315380 0 0 60720 131 314353 34 7 60630 132 369448 1558 25 55637 133 315380 0 0 60720 134 312846 43 20 60887 135 312075 645 4 60720 136 315009 316 4 60505 137 318903 115 10 60945 138 314887 5 1 60720 139 314913 897 4 60720 140 315380 0 0 60720 141 325506 389 8 58990 142 315380 0 0 60720 143 298568 1002 11 56750 144 315834 36 4 60894 145 329784 460 15 63346 146 312878 309 9 56535 147 315380 0 0 60720 148 314987 9 7 60835 149 325249 271 2 60720 150 315877 14 0 61016 151 291650 520 7 58650 152 305959 1766 46 60438 153 315380 0 5 60720 154 297765 458 7 58625 155 315245 20 2 60938 156 315380 0 0 60720 157 315380 0 0 60720 158 315236 98 2 61490 159 336425 405 5 60845 160 315380 0 0 60720 161 315380 0 0 60720 162 315380 0 0 60720 163 315380 0 0 60720 164 306268 483 7 60830 165 302187 454 24 63261 166 314882 47 1 60720 167 315380 0 0 60720 168 382712 757 18 45689 169 341570 4655 55 60720 170 315380 0 0 60720 171 315380 0 0 60720 172 312412 36 3 61564 173 315380 0 0 60720 174 309596 203 9 61938 175 315380 0 0 60720 176 315547 126 8 60951 177 313267 400 113 60720 178 316176 71 0 60745 179 315380 0 0 60720 180 315380 0 0 60720 181 359335 972 19 71642 182 330068 531 11 71641 183 314289 2461 25 55792 184 297413 378 16 71873 185 314806 23 5 62555 186 333210 638 11 60370 187 352108 2300 23 64873 188 313332 149 6 62041 189 291787 226 5 65745 190 315380 0 0 60720 191 318745 275 7 59500 192 315380 0 0 60720 193 315366 141 7 61630 194 315380 0 0 60720 195 315688 28 3 60890 196 315380 0 0 60720 197 409642 4980 89 113521 198 315380 0 0 60720 199 315380 0 0 60720 200 269587 472 19 80045 201 315380 0 0 60720 202 315380 0 0 60720 203 315380 0 0 60720 204 300962 203 12 50804 205 325479 496 12 87390 206 316155 10 5 61656 207 318574 63 2 65688 208 315380 0 0 60720 209 343613 1136 26 48522 210 306948 265 3 60720 211 315380 0 0 60720 212 315380 0 0 60720 213 330059 267 11 57640 214 288985 474 10 61977 215 304485 534 5 62620 216 315380 0 2 60720 217 315688 15 6 60831 218 317736 397 7 60646 219 315380 0 2 60720 220 322331 1866 28 56225 221 296656 288 3 60510 222 315380 0 0 60720 223 315354 3 1 60698 224 312161 468 20 60720 225 315576 20 1 60805 226 314922 278 22 61404 227 314551 61 9 60720 228 315380 0 0 60720 229 312339 192 2 65276 230 315380 0 0 60720 231 298700 317 7 63915 232 321376 738 9 60720 233 315380 0 0 60720 234 303230 368 13 61686 235 315380 0 0 60720 236 315487 2 0 60743 237 315380 0 0 60720 238 315793 53 6 60349 239 315380 0 0 60720 240 315380 0 0 60720 241 315380 0 0 60720 242 312887 94 3 61360 243 315380 0 0 60720 244 315637 24 7 59818 245 324385 2332 2 72680 246 315380 0 0 60720 247 315380 0 0 60720 248 308989 131 15 61808 249 315380 0 0 60720 250 315380 0 0 60720 251 296702 206 9 53110 252 315380 0 0 60720 253 307322 167 1 64245 254 304376 622 38 73007 255 253588 2328 57 82732 256 315380 0 0 60720 257 309560 365 7 54820 258 298466 364 26 47705 259 315380 0 0 60720 260 315380 0 0 60720 261 315380 0 0 60720 262 315380 0 0 60720 263 343929 226 13 72835 264 331955 307 10 58856 265 315380 0 0 60720 266 315380 0 0 60720 267 315380 0 0 60720 268 381180 188 9 77655 269 315380 0 0 60720 270 331420 138 26 69817 271 315380 0 0 60720 272 315380 0 0 60720 273 315380 0 0 60720 274 310201 125 19 60798 275 315380 0 0 60720 276 320016 282 12 62452 277 320398 335 23 64175 278 315380 0 0 60720 279 291841 1324 29 67440 280 310670 176 8 68136 281 315380 0 0 60720 282 315380 0 0 60720 283 313491 249 26 56726 284 315380 0 0 60720 285 331323 333 9 70811 286 315380 0 0 60720 287 319210 601 5 60720 288 318098 30 3 62045 289 315380 0 0 60720 290 292754 249 13 54323 291 315380 0 0 60720 292 325176 165 12 62841 293 365959 453 19 81125 294 315380 0 0 60720 295 302409 53 10 59506 296 340968 382 9 59365 297 315380 0 0 60720 298 315380 0 0 60720 299 315380 0 0 60720 300 315380 0 9 60720 301 313164 30 4 60798 302 301164 290 1 58790 303 315380 0 1 60720 304 315380 0 0 60720 305 344425 366 14 61808 306 315394 2 12 60735 307 315380 0 0 60720 308 316647 209 19 64016 309 309836 384 17 54683 310 315380 0 0 60720 311 315380 0 0 60720 312 346611 365 32 87192 313 315380 0 0 60720 314 322031 49 14 64107 315 315656 3 8 60761 316 339445 133 4 65990 317 314964 32 0 59988 318 297141 368 20 61167 319 315372 1 5 60719 320 315380 0 0 60720 321 315380 0 0 60720 322 315380 0 0 60720 323 315380 0 0 60720 324 315380 0 0 60720 325 315380 0 0 60720 326 312502 22 1 60722 327 315380 0 0 60720 328 315380 0 0 60720 329 315380 0 0 60720 330 315380 0 0 60720 331 315380 0 0 60720 332 315380 0 0 60720 333 315380 0 0 60720 334 313729 96 4 60379 335 315388 1 1 60727 336 315371 314 4 60720 337 296139 844 20 60925 338 315380 0 0 60720 339 313880 26 1 60896 340 317698 125 10 59734 341 295580 304 12 62969 342 315380 0 0 60720 343 315380 0 0 60720 344 315380 0 0 60720 345 308256 621 13 60720 346 315380 0 0 60720 347 303677 119 3 59118 348 315380 0 0 60720 349 315380 0 0 60720 350 319369 1595 10 60720 351 318690 312 3 58598 352 314049 60 7 61124 353 325699 587 10 59595 354 314210 135 1 62065 355 315380 0 0 60720 356 315380 0 0 60720 357 322378 514 15 78780 358 315380 0 0 60720 359 315380 0 0 60720 360 315380 0 0 60720 361 315398 1 4 60722 362 315380 0 0 60720 363 315380 0 0 60720 364 308336 1763 28 61600 365 316386 180 9 59635 366 315380 0 0 60720 367 315380 0 0 60720 368 315380 0 0 60720 369 315380 0 0 60720 370 315553 218 7 60720 371 315380 0 0 60720 372 323361 448 7 59781 373 336639 227 7 76644 374 307424 174 3 64820 375 315380 0 0 60720 376 315380 0 0 60720 377 295370 121 11 56178 378 322340 607 7 60436 379 319864 2212 10 60720 380 315380 0 0 60720 381 315380 0 0 60720 382 317291 530 18 73433 383 280398 571 14 41477 384 315380 0 0 60720 385 317330 78 12 62700 386 238125 2489 29 67804 387 327071 131 3 59661 388 309038 923 6 58620 389 314210 72 3 60398 390 307930 572 8 58580 391 322327 397 10 62710 392 292136 450 6 59325 393 263276 622 8 60950 394 367655 694 6 68060 395 283910 3425 9 83620 396 283587 562 8 58456 397 243650 4917 26 52811 398 438493 1442 239 121173 399 296261 529 7 63870 400 230621 2126 41 21001 401 304252 1061 3 70415 402 333505 776 8 64230 403 296919 611 6 59190 404 278990 1526 21 69351 405 276898 592 7 64270 406 327007 1182 11 70694 407 317046 621 11 68005 408 304555 989 12 58930 409 298096 438 9 58320 410 231861 726 3 69980 411 309422 1303 57 69863 412 286963 7419 21 63255 413 269753 1164 15 57320 414 448243 3310 32 75230 415 165404 1920 11 79420 416 204325 965 2 73490 417 407159 3256 23 35250 418 290476 1135 20 62285 419 275311 1270 24 69206 420 246541 661 1 65920 421 253468 1013 1 69770 422 240897 2844 74 72683 423 -83265 11528 68 -14545 424 -42143 6526 20 55830 425 272713 2264 20 55174 426 215362 5109 82 67038 427 42754 3999 21 51252 428 306275 35624 244 157278 429 253537 9252 32 79510 430 372631 15236 86 77440 431 -7170 18073 69 27284 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Kosten Oders Divid. 1.650e+05 2.202e+01 1.009e+03 1.568e+00 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2198539 -21927 41819 55153 2923814 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.650e+05 3.063e+04 5.386 1.19e-07 *** Kosten 2.202e+01 2.155e+00 10.217 < 2e-16 *** Oders 1.009e+03 3.497e+02 2.886 0.004096 ** Divid. 1.569e+00 4.115e-01 3.812 0.000158 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 253500 on 427 degrees of freedom Multiple R-squared: 0.6556, Adjusted R-squared: 0.6532 F-statistic: 271 on 3 and 427 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 2.430592e-32 1.215296e-32 [2,] 1.0000000 8.739682e-53 4.369841e-53 [3,] 1.0000000 2.606180e-53 1.303090e-53 [4,] 1.0000000 3.344876e-75 1.672438e-75 [5,] 1.0000000 3.789859e-78 1.894929e-78 [6,] 1.0000000 8.074852e-84 4.037426e-84 [7,] 1.0000000 3.228560e-105 1.614280e-105 [8,] 1.0000000 2.615482e-107 1.307741e-107 [9,] 1.0000000 8.428079e-116 4.214039e-116 [10,] 1.0000000 2.058138e-123 1.029069e-123 [11,] 1.0000000 1.769463e-133 8.847316e-134 [12,] 1.0000000 1.349026e-139 6.745128e-140 [13,] 1.0000000 1.674390e-144 8.371948e-145 [14,] 1.0000000 9.979506e-153 4.989753e-153 [15,] 1.0000000 4.714147e-161 2.357073e-161 [16,] 1.0000000 2.522661e-164 1.261330e-164 [17,] 1.0000000 7.254403e-165 3.627202e-165 [18,] 1.0000000 4.523952e-171 2.261976e-171 [19,] 1.0000000 1.148221e-170 5.741105e-171 [20,] 1.0000000 7.161661e-173 3.580830e-173 [21,] 1.0000000 6.188206e-176 3.094103e-176 [22,] 1.0000000 3.089147e-175 1.544573e-175 [23,] 1.0000000 4.706165e-189 2.353083e-189 [24,] 1.0000000 4.328376e-189 2.164188e-189 [25,] 1.0000000 6.442047e-195 3.221023e-195 [26,] 1.0000000 2.904653e-196 1.452326e-196 [27,] 1.0000000 1.152967e-205 5.764835e-206 [28,] 1.0000000 4.251535e-208 2.125767e-208 [29,] 1.0000000 2.166124e-207 1.083062e-207 [30,] 1.0000000 1.983912e-208 9.919558e-209 [31,] 1.0000000 4.716584e-214 2.358292e-214 [32,] 1.0000000 6.244748e-220 3.122374e-220 [33,] 1.0000000 4.578503e-219 2.289252e-219 [34,] 1.0000000 1.143927e-219 5.719637e-220 [35,] 1.0000000 3.044899e-221 1.522449e-221 [36,] 1.0000000 1.875420e-220 9.377099e-221 [37,] 1.0000000 8.053626e-221 4.026813e-221 [38,] 1.0000000 9.407400e-220 4.703700e-220 [39,] 1.0000000 2.088861e-219 1.044430e-219 [40,] 1.0000000 1.808484e-218 9.042421e-219 [41,] 1.0000000 4.585386e-219 2.292693e-219 [42,] 1.0000000 2.174066e-218 1.087033e-218 [43,] 1.0000000 8.208906e-218 4.104453e-218 [44,] 1.0000000 4.343402e-218 2.171701e-218 [45,] 1.0000000 1.574703e-218 7.873517e-219 [46,] 1.0000000 4.484256e-218 2.242128e-218 [47,] 1.0000000 2.730826e-217 1.365413e-217 [48,] 1.0000000 8.423766e-217 4.211883e-217 [49,] 1.0000000 5.051767e-220 2.525883e-220 [50,] 1.0000000 1.098346e-219 5.491728e-220 [51,] 1.0000000 3.518083e-233 1.759042e-233 [52,] 1.0000000 1.416428e-233 7.082138e-234 [53,] 1.0000000 1.301640e-232 6.508202e-233 [54,] 1.0000000 3.366758e-232 1.683379e-232 [55,] 1.0000000 6.717500e-234 3.358750e-234 [56,] 1.0000000 1.972486e-235 9.862431e-236 [57,] 1.0000000 5.452362e-235 2.726181e-235 [58,] 1.0000000 3.004395e-234 1.502198e-234 [59,] 1.0000000 1.671496e-233 8.357482e-234 [60,] 1.0000000 8.480176e-234 4.240088e-234 [61,] 1.0000000 2.273922e-233 1.136961e-233 [62,] 1.0000000 2.142692e-232 1.071346e-232 [63,] 1.0000000 1.214715e-231 6.073577e-232 [64,] 1.0000000 6.707206e-231 3.353603e-231 [65,] 1.0000000 4.723453e-230 2.361727e-230 [66,] 1.0000000 3.817530e-231 1.908765e-231 [67,] 1.0000000 1.310705e-230 6.553526e-231 [68,] 1.0000000 1.124158e-229 5.620790e-230 [69,] 1.0000000 1.349869e-228 6.749343e-229 [70,] 1.0000000 3.164417e-228 1.582209e-228 [71,] 1.0000000 2.854783e-227 1.427392e-227 [72,] 1.0000000 3.837427e-227 1.918714e-227 [73,] 1.0000000 4.036207e-226 2.018104e-226 [74,] 1.0000000 1.525577e-225 7.627887e-226 [75,] 1.0000000 9.075956e-225 4.537978e-225 [76,] 1.0000000 1.497523e-224 7.487613e-225 [77,] 1.0000000 1.459550e-224 7.297752e-225 [78,] 1.0000000 8.179804e-224 4.089902e-224 [79,] 1.0000000 7.945820e-223 3.972910e-223 [80,] 1.0000000 2.123714e-222 1.061857e-222 [81,] 1.0000000 1.694008e-222 8.470042e-223 [82,] 1.0000000 1.267032e-221 6.335159e-222 [83,] 1.0000000 2.738350e-221 1.369175e-221 [84,] 1.0000000 2.664020e-220 1.332010e-220 [85,] 1.0000000 2.421778e-219 1.210889e-219 [86,] 1.0000000 1.577767e-218 7.888834e-219 [87,] 1.0000000 5.308984e-220 2.654492e-220 [88,] 1.0000000 5.402230e-219 2.701115e-219 [89,] 1.0000000 1.923829e-219 9.619147e-220 [90,] 1.0000000 1.707649e-219 8.538247e-220 [91,] 1.0000000 1.533171e-218 7.665856e-219 [92,] 1.0000000 1.506117e-217 7.530586e-218 [93,] 1.0000000 1.317935e-216 6.589677e-217 [94,] 1.0000000 1.295777e-215 6.478886e-216 [95,] 1.0000000 1.170218e-214 5.851090e-215 [96,] 1.0000000 9.191774e-216 4.595887e-216 [97,] 1.0000000 8.432706e-216 4.216353e-216 [98,] 1.0000000 1.403658e-216 7.018291e-217 [99,] 1.0000000 6.191573e-216 3.095787e-216 [100,] 1.0000000 2.187586e-215 1.093793e-215 [101,] 1.0000000 1.346190e-217 6.730948e-218 [102,] 1.0000000 5.649664e-219 2.824832e-219 [103,] 1.0000000 4.299186e-218 2.149593e-218 [104,] 1.0000000 8.620852e-218 4.310426e-218 [105,] 1.0000000 6.134868e-218 3.067434e-218 [106,] 1.0000000 3.811163e-218 1.905582e-218 [107,] 1.0000000 2.946839e-217 1.473420e-217 [108,] 1.0000000 2.387761e-216 1.193881e-216 [109,] 1.0000000 1.584638e-216 7.923190e-217 [110,] 1.0000000 6.365781e-216 3.182891e-216 [111,] 1.0000000 4.585815e-215 2.292908e-215 [112,] 1.0000000 3.686476e-215 1.843238e-215 [113,] 1.0000000 3.978613e-216 1.989306e-216 [114,] 1.0000000 4.088253e-215 2.044126e-215 [115,] 1.0000000 1.975573e-214 9.877866e-215 [116,] 1.0000000 1.832932e-213 9.164659e-214 [117,] 1.0000000 1.589629e-212 7.948143e-213 [118,] 1.0000000 2.198678e-212 1.099339e-212 [119,] 1.0000000 6.318306e-212 3.159153e-212 [120,] 1.0000000 2.063023e-212 1.031511e-212 [121,] 1.0000000 1.926414e-211 9.632071e-212 [122,] 1.0000000 1.902403e-210 9.512017e-211 [123,] 1.0000000 1.818865e-209 9.094323e-210 [124,] 1.0000000 1.710070e-208 8.550349e-209 [125,] 1.0000000 1.668946e-207 8.344731e-208 [126,] 1.0000000 1.435411e-207 7.177056e-208 [127,] 1.0000000 1.367575e-206 6.837874e-207 [128,] 1.0000000 1.417370e-205 7.086848e-206 [129,] 1.0000000 1.401973e-204 7.009863e-205 [130,] 1.0000000 1.358313e-203 6.791565e-204 [131,] 1.0000000 1.322004e-202 6.610022e-203 [132,] 1.0000000 1.269751e-201 6.348753e-202 [133,] 1.0000000 1.214274e-200 6.071372e-201 [134,] 1.0000000 1.158771e-199 5.793855e-200 [135,] 1.0000000 9.589873e-199 4.794937e-199 [136,] 1.0000000 9.151115e-198 4.575558e-198 [137,] 1.0000000 9.369098e-197 4.684549e-197 [138,] 1.0000000 9.036756e-196 4.518378e-196 [139,] 1.0000000 8.271089e-195 4.135544e-195 [140,] 1.0000000 7.584350e-194 3.792175e-194 [141,] 1.0000000 7.191716e-193 3.595858e-193 [142,] 1.0000000 6.975854e-192 3.487927e-192 [143,] 1.0000000 6.006059e-191 3.003029e-191 [144,] 1.0000000 5.677643e-190 2.838821e-190 [145,] 1.0000000 5.527930e-189 2.763965e-189 [146,] 1.0000000 5.706216e-188 2.853108e-188 [147,] 1.0000000 5.419894e-187 2.709947e-187 [148,] 1.0000000 5.317814e-186 2.658907e-186 [149,] 1.0000000 4.999803e-185 2.499902e-185 [150,] 1.0000000 4.652301e-184 2.326150e-184 [151,] 1.0000000 4.320326e-183 2.160163e-183 [152,] 1.0000000 4.057165e-182 2.028583e-182 [153,] 1.0000000 2.837675e-181 1.418837e-181 [154,] 1.0000000 2.619488e-180 1.309744e-180 [155,] 1.0000000 2.412398e-179 1.206199e-179 [156,] 1.0000000 2.216283e-178 1.108142e-178 [157,] 1.0000000 2.031012e-177 1.015506e-177 [158,] 1.0000000 1.939728e-176 9.698638e-177 [159,] 1.0000000 1.809404e-175 9.047021e-176 [160,] 1.0000000 1.649790e-174 8.248948e-175 [161,] 1.0000000 1.493476e-173 7.467381e-174 [162,] 1.0000000 1.850795e-174 9.253976e-175 [163,] 1.0000000 3.303986e-174 1.651993e-174 [164,] 1.0000000 3.027321e-173 1.513660e-173 [165,] 1.0000000 2.765752e-172 1.382876e-172 [166,] 1.0000000 2.574255e-171 1.287127e-171 [167,] 1.0000000 2.337443e-170 1.168721e-170 [168,] 1.0000000 2.194040e-169 1.097020e-169 [169,] 1.0000000 1.979190e-168 9.895948e-169 [170,] 1.0000000 1.805319e-167 9.026594e-168 [171,] 1.0000000 1.695989e-166 8.479946e-167 [172,] 1.0000000 1.505326e-165 7.526629e-166 [173,] 1.0000000 1.338383e-164 6.691913e-165 [174,] 1.0000000 1.186093e-163 5.930463e-164 [175,] 1.0000000 7.911250e-163 3.955625e-163 [176,] 1.0000000 7.250046e-162 3.625023e-162 [177,] 1.0000000 4.225910e-161 2.112955e-161 [178,] 1.0000000 2.568503e-160 1.284252e-160 [179,] 1.0000000 2.293050e-159 1.146525e-159 [180,] 1.0000000 1.530856e-158 7.654282e-159 [181,] 1.0000000 5.726082e-158 2.863041e-158 [182,] 1.0000000 5.092835e-157 2.546418e-157 [183,] 1.0000000 3.600176e-156 1.800088e-156 [184,] 1.0000000 3.101590e-155 1.550795e-155 [185,] 1.0000000 2.539664e-154 1.269832e-154 [186,] 1.0000000 2.172516e-153 1.086258e-153 [187,] 1.0000000 1.883443e-152 9.417215e-153 [188,] 1.0000000 1.599390e-151 7.996952e-152 [189,] 1.0000000 1.358950e-150 6.794751e-151 [190,] 1.0000000 1.145603e-149 5.728014e-150 [191,] 1.0000000 8.787550e-149 4.393775e-149 [192,] 1.0000000 7.351060e-148 3.675530e-148 [193,] 1.0000000 6.126847e-147 3.063423e-147 [194,] 1.0000000 6.344816e-147 3.172408e-147 [195,] 1.0000000 5.301083e-146 2.650541e-146 [196,] 1.0000000 4.412611e-145 2.206306e-145 [197,] 1.0000000 3.659338e-144 1.829669e-144 [198,] 1.0000000 2.808778e-143 1.404389e-143 [199,] 1.0000000 1.235292e-142 6.176459e-143 [200,] 1.0000000 1.027787e-141 5.138934e-142 [201,] 1.0000000 8.550471e-141 4.275235e-141 [202,] 1.0000000 6.981267e-140 3.490634e-140 [203,] 1.0000000 1.005138e-139 5.025692e-140 [204,] 1.0000000 8.403893e-139 4.201947e-139 [205,] 1.0000000 6.831375e-138 3.415688e-138 [206,] 1.0000000 5.531530e-137 2.765765e-137 [207,] 1.0000000 3.352656e-136 1.676328e-136 [208,] 1.0000000 2.520675e-135 1.260337e-135 [209,] 1.0000000 2.070219e-134 1.035110e-134 [210,] 1.0000000 1.655285e-133 8.276427e-134 [211,] 1.0000000 1.322122e-132 6.610610e-133 [212,] 1.0000000 1.021302e-131 5.106509e-132 [213,] 1.0000000 8.066985e-131 4.033492e-131 [214,] 1.0000000 3.866636e-130 1.933318e-130 [215,] 1.0000000 3.046118e-129 1.523059e-129 [216,] 1.0000000 2.377357e-128 1.188679e-128 [217,] 1.0000000 1.849042e-127 9.245209e-128 [218,] 1.0000000 1.451266e-126 7.256329e-127 [219,] 1.0000000 1.119107e-125 5.595536e-126 [220,] 1.0000000 8.723133e-125 4.361566e-125 [221,] 1.0000000 6.723391e-124 3.361695e-124 [222,] 1.0000000 5.117062e-123 2.558531e-123 [223,] 1.0000000 3.917908e-122 1.958954e-122 [224,] 1.0000000 2.957740e-121 1.478870e-121 [225,] 1.0000000 2.147152e-120 1.073576e-120 [226,] 1.0000000 1.481676e-119 7.408380e-120 [227,] 1.0000000 1.105263e-118 5.526315e-119 [228,] 1.0000000 8.419570e-118 4.209785e-118 [229,] 1.0000000 6.227574e-117 3.113787e-117 [230,] 1.0000000 4.585707e-116 2.292853e-116 [231,] 1.0000000 3.363278e-115 1.681639e-115 [232,] 1.0000000 2.445823e-114 1.222912e-114 [233,] 1.0000000 1.778540e-113 8.892699e-114 [234,] 1.0000000 1.287785e-112 6.438925e-113 [235,] 1.0000000 9.284482e-112 4.642241e-112 [236,] 1.0000000 6.759561e-111 3.379780e-111 [237,] 1.0000000 4.831331e-110 2.415665e-110 [238,] 1.0000000 3.405378e-109 1.702689e-109 [239,] 1.0000000 2.471831e-108 1.235915e-108 [240,] 1.0000000 1.743174e-107 8.715870e-108 [241,] 1.0000000 1.223942e-106 6.119712e-107 [242,] 1.0000000 8.778635e-106 4.389318e-106 [243,] 1.0000000 6.109333e-105 3.054667e-105 [244,] 1.0000000 4.232923e-104 2.116461e-104 [245,] 1.0000000 2.900174e-103 1.450087e-103 [246,] 1.0000000 1.991501e-102 9.957507e-103 [247,] 1.0000000 1.367456e-101 6.837281e-102 [248,] 1.0000000 7.697318e-101 3.848659e-101 [249,] 1.0000000 6.133853e-101 3.066926e-101 [250,] 1.0000000 4.198894e-100 2.099447e-100 [251,] 1.0000000 2.626027e-99 1.313014e-99 [252,] 1.0000000 1.494175e-98 7.470875e-99 [253,] 1.0000000 1.010003e-97 5.050017e-98 [254,] 1.0000000 6.796205e-97 3.398102e-97 [255,] 1.0000000 4.552238e-96 2.276119e-96 [256,] 1.0000000 3.035225e-95 1.517612e-95 [257,] 1.0000000 2.040084e-94 1.020042e-94 [258,] 1.0000000 1.070511e-93 5.352554e-94 [259,] 1.0000000 7.045406e-93 3.522703e-93 [260,] 1.0000000 4.615288e-92 2.307644e-92 [261,] 1.0000000 3.009272e-91 1.504636e-91 [262,] 1.0000000 1.604508e-90 8.022538e-91 [263,] 1.0000000 1.037089e-89 5.185444e-90 [264,] 1.0000000 6.846043e-89 3.423022e-89 [265,] 1.0000000 4.383118e-88 2.191559e-88 [266,] 1.0000000 2.792913e-87 1.396457e-87 [267,] 1.0000000 1.771145e-86 8.855723e-87 [268,] 1.0000000 1.142786e-85 5.713931e-86 [269,] 1.0000000 7.177221e-85 3.588611e-85 [270,] 1.0000000 4.469982e-84 2.234991e-84 [271,] 1.0000000 2.824435e-83 1.412218e-83 [272,] 1.0000000 1.747989e-82 8.739945e-83 [273,] 1.0000000 9.987421e-82 4.993710e-82 [274,] 1.0000000 5.987048e-81 2.993524e-81 [275,] 1.0000000 3.657469e-80 1.828735e-80 [276,] 1.0000000 2.223306e-79 1.111653e-79 [277,] 1.0000000 1.274009e-78 6.370044e-79 [278,] 1.0000000 7.669229e-78 3.834615e-78 [279,] 1.0000000 4.693385e-77 2.346693e-77 [280,] 1.0000000 2.797238e-76 1.398619e-76 [281,] 1.0000000 1.582610e-75 7.913049e-76 [282,] 1.0000000 9.371309e-75 4.685654e-75 [283,] 1.0000000 5.501093e-74 2.750546e-74 [284,] 1.0000000 3.270733e-73 1.635366e-73 [285,] 1.0000000 1.900239e-72 9.501195e-73 [286,] 1.0000000 1.074963e-71 5.374813e-72 [287,] 1.0000000 6.230254e-71 3.115127e-71 [288,] 1.0000000 3.563613e-70 1.781807e-70 [289,] 1.0000000 2.076082e-69 1.038041e-69 [290,] 1.0000000 8.466843e-69 4.233422e-69 [291,] 1.0000000 4.775010e-68 2.387505e-68 [292,] 1.0000000 2.678632e-67 1.339316e-67 [293,] 1.0000000 1.494603e-66 7.473014e-67 [294,] 1.0000000 8.325508e-66 4.162754e-66 [295,] 1.0000000 4.640352e-65 2.320176e-65 [296,] 1.0000000 2.602008e-64 1.301004e-64 [297,] 1.0000000 1.421324e-63 7.106620e-64 [298,] 1.0000000 7.717969e-63 3.858985e-63 [299,] 1.0000000 3.200789e-62 1.600394e-62 [300,] 1.0000000 1.728956e-61 8.644781e-62 [301,] 1.0000000 9.245426e-61 4.622713e-61 [302,] 1.0000000 5.032402e-60 2.516201e-60 [303,] 1.0000000 2.446646e-59 1.223323e-59 [304,] 1.0000000 1.286819e-58 6.434096e-59 [305,] 1.0000000 6.729447e-58 3.364724e-58 [306,] 1.0000000 3.171111e-57 1.585556e-57 [307,] 1.0000000 1.642811e-56 8.214056e-57 [308,] 1.0000000 8.563712e-56 4.281856e-56 [309,] 1.0000000 4.393795e-55 2.196897e-55 [310,] 1.0000000 2.102368e-54 1.051184e-54 [311,] 1.0000000 1.055659e-53 5.278297e-54 [312,] 1.0000000 5.449743e-53 2.724871e-53 [313,] 1.0000000 2.731083e-52 1.365541e-52 [314,] 1.0000000 1.357996e-51 6.789981e-52 [315,] 1.0000000 6.711486e-51 3.355743e-51 [316,] 1.0000000 3.296681e-50 1.648340e-50 [317,] 1.0000000 1.609365e-49 8.046823e-50 [318,] 1.0000000 7.807862e-49 3.903931e-49 [319,] 1.0000000 3.764353e-48 1.882176e-48 [320,] 1.0000000 1.822163e-47 9.110813e-48 [321,] 1.0000000 8.674707e-47 4.337353e-47 [322,] 1.0000000 4.103382e-46 2.051691e-46 [323,] 1.0000000 1.928533e-45 9.642665e-46 [324,] 1.0000000 9.005099e-45 4.502550e-45 [325,] 1.0000000 4.177378e-44 2.088689e-44 [326,] 1.0000000 1.925085e-43 9.625426e-44 [327,] 1.0000000 8.812592e-43 4.406296e-43 [328,] 1.0000000 4.013298e-42 2.006649e-42 [329,] 1.0000000 1.813349e-41 9.066746e-42 [330,] 1.0000000 8.058085e-41 4.029042e-41 [331,] 1.0000000 3.710338e-40 1.855169e-40 [332,] 1.0000000 1.641984e-39 8.209921e-40 [333,] 1.0000000 7.267902e-39 3.633951e-39 [334,] 1.0000000 3.095162e-38 1.547581e-38 [335,] 1.0000000 1.374656e-37 6.873279e-38 [336,] 1.0000000 5.919262e-37 2.959631e-37 [337,] 1.0000000 2.530407e-36 1.265204e-36 [338,] 1.0000000 1.073822e-35 5.369108e-36 [339,] 1.0000000 4.577995e-35 2.288997e-35 [340,] 1.0000000 1.914210e-34 9.571050e-35 [341,] 1.0000000 8.121518e-34 4.060759e-34 [342,] 1.0000000 3.345165e-33 1.672583e-33 [343,] 1.0000000 1.367273e-32 6.836365e-33 [344,] 1.0000000 5.011547e-32 2.505774e-32 [345,] 1.0000000 1.897232e-31 9.486159e-32 [346,] 1.0000000 7.657861e-31 3.828930e-31 [347,] 1.0000000 2.728176e-30 1.364088e-30 [348,] 1.0000000 1.088639e-29 5.443194e-30 [349,] 1.0000000 4.246226e-29 2.123113e-29 [350,] 1.0000000 1.642496e-28 8.212478e-29 [351,] 1.0000000 6.318022e-28 3.159011e-28 [352,] 1.0000000 2.409193e-27 1.204596e-27 [353,] 1.0000000 9.107831e-27 4.553916e-27 [354,] 1.0000000 3.413206e-26 1.706603e-26 [355,] 1.0000000 1.270191e-25 6.350954e-26 [356,] 1.0000000 4.676202e-25 2.338101e-25 [357,] 1.0000000 1.705948e-24 8.529739e-25 [358,] 1.0000000 6.081439e-24 3.040719e-24 [359,] 1.0000000 2.134975e-23 1.067487e-23 [360,] 1.0000000 7.577712e-23 3.788856e-23 [361,] 1.0000000 2.663818e-22 1.331909e-22 [362,] 1.0000000 9.273145e-22 4.636572e-22 [363,] 1.0000000 3.196226e-21 1.598113e-21 [364,] 1.0000000 1.085630e-20 5.428148e-21 [365,] 1.0000000 3.666904e-20 1.833452e-20 [366,] 1.0000000 1.139622e-19 5.698109e-20 [367,] 1.0000000 3.983320e-19 1.991660e-19 [368,] 1.0000000 1.374135e-18 6.870676e-19 [369,] 1.0000000 4.456475e-18 2.228238e-18 [370,] 1.0000000 1.429197e-17 7.145986e-18 [371,] 1.0000000 4.676550e-17 2.338275e-17 [372,] 1.0000000 1.379755e-16 6.898773e-17 [373,] 1.0000000 3.754474e-16 1.877237e-16 [374,] 1.0000000 1.150483e-15 5.752416e-16 [375,] 1.0000000 3.481903e-15 1.740951e-15 [376,] 1.0000000 1.120447e-14 5.602236e-15 [377,] 1.0000000 2.962514e-14 1.481257e-14 [378,] 1.0000000 8.658496e-14 4.329248e-14 [379,] 1.0000000 2.549114e-13 1.274557e-13 [380,] 1.0000000 6.812795e-13 3.406397e-13 [381,] 1.0000000 1.781479e-12 8.907395e-13 [382,] 1.0000000 4.844002e-12 2.422001e-12 [383,] 1.0000000 1.336663e-11 6.683313e-12 [384,] 1.0000000 3.604832e-11 1.802416e-11 [385,] 1.0000000 9.508511e-11 4.754255e-11 [386,] 1.0000000 2.667857e-10 1.333928e-10 [387,] 1.0000000 7.603123e-10 3.801562e-10 [388,] 1.0000000 1.526677e-09 7.633383e-10 [389,] 1.0000000 4.142153e-09 2.071077e-09 [390,] 1.0000000 1.118990e-08 5.594948e-09 [391,] 1.0000000 2.889582e-08 1.444791e-08 [392,] 1.0000000 5.852297e-08 2.926148e-08 [393,] 0.9999999 1.525094e-07 7.625468e-08 [394,] 0.9999998 3.216808e-07 1.608404e-07 [395,] 0.9999996 8.186697e-07 4.093349e-07 [396,] 0.9999991 1.759357e-06 8.796784e-07 [397,] 0.9999980 4.096960e-06 2.048480e-06 [398,] 0.9999949 1.011551e-05 5.057754e-06 [399,] 0.9999878 2.439833e-05 1.219916e-05 [400,] 0.9999732 5.355622e-05 2.677811e-05 [401,] 0.9999412 1.176231e-04 5.881155e-05 [402,] 0.9998805 2.390808e-04 1.195404e-04 [403,] 0.9997563 4.874979e-04 2.437490e-04 [404,] 0.9994846 1.030896e-03 5.154479e-04 [405,] 0.9989024 2.195194e-03 1.097597e-03 [406,] 0.9982032 3.593522e-03 1.796761e-03 [407,] 0.9964342 7.131593e-03 3.565796e-03 [408,] 0.9968653 6.269431e-03 3.134716e-03 [409,] 0.9949696 1.006073e-02 5.030367e-03 [410,] 0.9905244 1.895117e-02 9.475586e-03 [411,] 0.9977198 4.560433e-03 2.280217e-03 [412,] 0.9952421 9.515868e-03 4.757934e-03 [413,] 0.9890591 2.188182e-02 1.094091e-02 [414,] 0.9754772 4.904553e-02 2.452276e-02 [415,] 0.9482408 1.035184e-01 5.175922e-02 [416,] 0.8942305 2.115390e-01 1.057695e-01 [417,] 0.8019684 3.960632e-01 1.980316e-01 [418,] 0.8245668 3.508663e-01 1.754332e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1uzlk1291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/2m8k51291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/3m8k51291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/4m8k51291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5fzj81291211442.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 1389636.7423 2923814.3860 1363346.2455 -2198539.3196 -394043.3532 6 7 8 9 10 109869.7630 -135939.7551 128928.3808 -336334.3876 715092.5993 11 12 13 14 15 -150657.5105 342307.7668 449675.2072 -1139741.2853 -214255.2259 16 17 18 19 20 460574.8057 -237906.6182 -238071.5057 -295877.9417 -585102.6254 21 22 23 24 25 203463.1027 -137416.1514 -119601.7096 -638049.8630 -318324.0983 26 27 28 29 30 -84297.5193 -471896.4012 -340384.4471 348299.6853 -261059.4968 31 32 33 34 35 286668.0415 -87937.6349 255919.3586 130738.8627 -54941.9881 36 37 38 39 40 -569143.8636 34310.7675 157833.6426 -51015.8179 -304565.6193 41 42 43 44 45 -459620.0553 -329011.9718 172.3323 -209375.2579 -27239.0962 46 47 48 49 50 -53789.6599 -281553.3224 -264256.0037 -53547.0125 -181997.0873 51 52 53 54 55 35733.0747 -86617.6159 -228211.1136 44146.0101 -394458.1521 56 57 58 59 60 14558.8138 320392.3041 -22047.0390 -59458.7651 -205533.6443 61 62 63 64 65 14509.4419 -179052.6824 7020.6734 -174141.9735 -15771.8893 66 67 68 69 70 -494796.9414 -65004.9696 -14962.5874 -100725.8195 -72004.9465 71 72 73 74 75 -122294.6437 1252.3272 -88696.6011 15111.9538 -34736.9165 76 77 78 79 80 -98722.3509 -28470.4246 -73253.8756 -42228.7857 -37810.1071 81 82 83 84 85 27095.5975 -49657.3672 -82702.7231 41679.9422 24950.3027 86 87 88 89 90 -102379.3908 -29329.7313 1171.4484 -115198.1764 32516.3322 91 92 93 94 95 -1902.9294 -21807.0740 69502.3235 -172301.4464 49811.5996 96 97 98 99 100 19619.0604 71821.8316 2140.2336 -41173.8667 -1473.8500 101 102 103 104 105 -35893.9515 -112498.0489 -165507.6930 -50102.0373 30070.0947 106 107 108 109 110 -53263.1917 -252162.8870 131136.6578 -175232.6525 -92237.8474 111 112 113 114 115 94519.8432 -190793.4409 -165268.7848 -57466.9920 86189.5153 116 117 118 119 120 21656.0303 17524.5950 72490.9475 -132145.5451 14988.2296 121 122 123 124 125 -288568.0383 32827.5084 -103083.9599 84905.4238 13946.7576 126 127 128 129 130 -162469.9198 55152.5245 42028.6252 29692.3360 55152.5245 131 132 133 134 135 46453.0982 57660.3793 55152.5245 31224.0439 33610.3153 136 137 138 139 140 44124.6501 45697.8939 53540.1528 30900.3954 55152.5245 141 142 143 144 145 51353.5430 55152.5245 11405.4575 50503.8779 40171.2046 146 147 148 149 150 43328.0424 55152.5245 47315.9545 57036.7210 54877.0451 151 152 153 154 155 16156.0728 -39133.1674 50106.0545 23675.2486 52216.7029 156 157 158 159 160 55152.5245 55152.5245 49624.7042 62038.6981 55152.5245 161 162 163 164 165 55152.5245 55152.5245 55152.5245 28169.4235 3755.9701 166 167 168 169 170 52610.4994 55152.5245 111226.8553 -76651.0561 55152.5245 171 172 173 174 175 55152.5245 47040.3092 55152.5245 33905.3464 55152.5245 176 177 178 179 180 44108.9001 -69816.9214 54346.2088 55152.5245 55152.5245 181 182 183 184 185 41401.1884 29918.9688 -17621.7771 -4778.0021 46147.5877 186 187 188 189 190 48383.3242 11517.2049 41696.5179 13656.0689 55152.5245 191 192 193 194 195 47311.6999 55152.5245 43541.9820 55152.5245 51549.5701 196 197 198 199 200 55152.5245 -132865.9295 55152.5245 55152.5245 -50518.7341 201 202 203 204 205 55152.5245 55152.5245 55152.5245 39706.6063 389.6725 206 207 208 209 210 49192.8281 47148.8882 55152.5245 51266.1839 37858.5204 211 212 213 214 215 55152.5245 55152.5245 57681.9685 6257.6673 24474.6896 216 217 218 219 220 53133.9365 48900.4290 41819.3564 53133.9365 -187.5434 221 222 223 224 225 27389.5375 55152.5245 54085.6898 21444.3644 53765.6010 226 227 228 229 230 25296.9108 43896.9295 55152.5245 38720.0837 55152.5245 231 232 233 234 235 19417.3320 35817.3986 55152.5245 20264.8536 55152.5245 236 237 238 239 240 55179.4190 55152.5245 48924.8316 55152.5245 55152.5245 241 242 243 244 245 55152.5245 46558.3712 55152.5245 49230.8366 -7960.0214 246 247 248 249 250 55152.5245 55152.5245 29031.6010 55152.5245 55152.5245 251 252 253 254 255 34791.5919 55152.5245 36879.8403 -27169.8902 -149946.2072 256 257 258 259 260 55152.5245 43485.6542 24396.6177 55152.5245 55152.5245 261 262 263 264 265 55152.5245 55152.5245 46603.3938 57799.4030 55152.5245 266 267 268 269 270 55152.5245 55152.5245 81168.2234 55152.5245 27644.5278 271 272 273 274 275 55152.5245 55152.5245 55152.5245 27922.6546 55152.5245 276 277 278 279 280 38752.0513 24162.5503 55152.5245 -37344.5929 26861.7964 281 282 283 284 285 55152.5245 55152.5245 27804.4032 55152.5245 38853.4518 286 287 288 289 290 55152.5245 40704.7058 52103.9773 55152.5245 23957.2150 291 292 293 294 295 55152.5245 45877.7435 44577.6335 55152.5245 32825.8604 296 297 298 299 300 65372.1884 55152.5245 55152.5245 55152.5245 46068.8784 301 302 303 304 305 48116.5426 36569.8314 54143.2305 55152.5245 60303.2396 306 307 308 309 310 42987.4385 55152.5245 27472.0702 33465.2963 55152.5245 311 312 313 314 315 55152.5245 4530.3833 55152.5245 41282.3000 47223.8192 316 317 318 319 320 63986.5389 55180.1335 7924.8219 50077.6074 55152.5245 321 322 323 324 325 55152.5245 55152.5245 55152.5245 55152.5245 55152.5245 326 327 328 329 330 50777.7514 55152.5245 55152.5245 55152.5245 55152.5245 331 332 333 334 335 55152.5245 55152.5245 55152.5245 47885.6972 54118.2358 336 337 338 339 340 44193.4641 -3177.0170 55152.5245 51794.7786 46172.1334 341 342 343 344 345 13020.8195 55152.5245 55152.5245 55152.5245 21236.0425 346 347 348 349 350 55152.5245 40314.4511 55152.5245 55152.5245 13933.7738 351 352 353 354 355 51894.0439 44801.8786 44219.9617 47891.5630 55152.5245 356 357 358 359 360 55152.5245 7368.8801 55152.5245 55152.5245 55152.5245 361 362 363 364 365 51108.1960 55152.5245 55152.5245 -20345.3694 44813.8486 366 367 368 369 370 55152.5245 55152.5245 55152.5245 55152.5245 43461.0754 371 372 373 374 375 55152.5245 47678.2739 39372.9095 33907.2837 55152.5245 376 377 378 379 380 55152.5245 28500.3161 42129.4647 845.1762 55152.5245 381 382 383 384 385 55152.5245 7288.2607 23651.2435 55152.5245 40168.2486 386 387 388 389 390 -117279.6312 62592.5952 25728.1516 49874.5640 30391.7616 391 392 393 394 395 40145.1899 18133.7508 -19080.2469 74580.5292 -96721.9437 396 397 398 399 400 6463.4051 -138664.7201 -189519.7904 12381.6147 -55495.2651 401 402 403 404 405 2431.9986 42613.8359 19583.9874 -49565.6943 -8995.7460 406 407 408 409 410 14011.1662 20618.4595 13250.1409 22906.3494 -61901.5139 411 412 413 414 415 -51361.8573 -161769.0794 -25907.2560 60088.3931 -177525.6236 416 417 418 419 420 -99195.2033 91983.5806 -17379.6380 -50409.1687 -37404.0006 421 422 423 424 425 -44265.0154 -175393.8602 -547870.2678 -458560.1402 -48844.9394 426 427 428 429 430 -250014.5348 -311858.7540 -1135948.9209 -272147.0098 -336049.2741 431 -682483.1396 > postscript(file="/var/www/html/freestat/rcomp/tmp/6fzj81291211442.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 1389636.7423 NA 1 2923814.3860 1389636.7423 2 1363346.2455 2923814.3860 3 -2198539.3196 1363346.2455 4 -394043.3532 -2198539.3196 5 109869.7630 -394043.3532 6 -135939.7551 109869.7630 7 128928.3808 -135939.7551 8 -336334.3876 128928.3808 9 715092.5993 -336334.3876 10 -150657.5105 715092.5993 11 342307.7668 -150657.5105 12 449675.2072 342307.7668 13 -1139741.2853 449675.2072 14 -214255.2259 -1139741.2853 15 460574.8057 -214255.2259 16 -237906.6182 460574.8057 17 -238071.5057 -237906.6182 18 -295877.9417 -238071.5057 19 -585102.6254 -295877.9417 20 203463.1027 -585102.6254 21 -137416.1514 203463.1027 22 -119601.7096 -137416.1514 23 -638049.8630 -119601.7096 24 -318324.0983 -638049.8630 25 -84297.5193 -318324.0983 26 -471896.4012 -84297.5193 27 -340384.4471 -471896.4012 28 348299.6853 -340384.4471 29 -261059.4968 348299.6853 30 286668.0415 -261059.4968 31 -87937.6349 286668.0415 32 255919.3586 -87937.6349 33 130738.8627 255919.3586 34 -54941.9881 130738.8627 35 -569143.8636 -54941.9881 36 34310.7675 -569143.8636 37 157833.6426 34310.7675 38 -51015.8179 157833.6426 39 -304565.6193 -51015.8179 40 -459620.0553 -304565.6193 41 -329011.9718 -459620.0553 42 172.3323 -329011.9718 43 -209375.2579 172.3323 44 -27239.0962 -209375.2579 45 -53789.6599 -27239.0962 46 -281553.3224 -53789.6599 47 -264256.0037 -281553.3224 48 -53547.0125 -264256.0037 49 -181997.0873 -53547.0125 50 35733.0747 -181997.0873 51 -86617.6159 35733.0747 52 -228211.1136 -86617.6159 53 44146.0101 -228211.1136 54 -394458.1521 44146.0101 55 14558.8138 -394458.1521 56 320392.3041 14558.8138 57 -22047.0390 320392.3041 58 -59458.7651 -22047.0390 59 -205533.6443 -59458.7651 60 14509.4419 -205533.6443 61 -179052.6824 14509.4419 62 7020.6734 -179052.6824 63 -174141.9735 7020.6734 64 -15771.8893 -174141.9735 65 -494796.9414 -15771.8893 66 -65004.9696 -494796.9414 67 -14962.5874 -65004.9696 68 -100725.8195 -14962.5874 69 -72004.9465 -100725.8195 70 -122294.6437 -72004.9465 71 1252.3272 -122294.6437 72 -88696.6011 1252.3272 73 15111.9538 -88696.6011 74 -34736.9165 15111.9538 75 -98722.3509 -34736.9165 76 -28470.4246 -98722.3509 77 -73253.8756 -28470.4246 78 -42228.7857 -73253.8756 79 -37810.1071 -42228.7857 80 27095.5975 -37810.1071 81 -49657.3672 27095.5975 82 -82702.7231 -49657.3672 83 41679.9422 -82702.7231 84 24950.3027 41679.9422 85 -102379.3908 24950.3027 86 -29329.7313 -102379.3908 87 1171.4484 -29329.7313 88 -115198.1764 1171.4484 89 32516.3322 -115198.1764 90 -1902.9294 32516.3322 91 -21807.0740 -1902.9294 92 69502.3235 -21807.0740 93 -172301.4464 69502.3235 94 49811.5996 -172301.4464 95 19619.0604 49811.5996 96 71821.8316 19619.0604 97 2140.2336 71821.8316 98 -41173.8667 2140.2336 99 -1473.8500 -41173.8667 100 -35893.9515 -1473.8500 101 -112498.0489 -35893.9515 102 -165507.6930 -112498.0489 103 -50102.0373 -165507.6930 104 30070.0947 -50102.0373 105 -53263.1917 30070.0947 106 -252162.8870 -53263.1917 107 131136.6578 -252162.8870 108 -175232.6525 131136.6578 109 -92237.8474 -175232.6525 110 94519.8432 -92237.8474 111 -190793.4409 94519.8432 112 -165268.7848 -190793.4409 113 -57466.9920 -165268.7848 114 86189.5153 -57466.9920 115 21656.0303 86189.5153 116 17524.5950 21656.0303 117 72490.9475 17524.5950 118 -132145.5451 72490.9475 119 14988.2296 -132145.5451 120 -288568.0383 14988.2296 121 32827.5084 -288568.0383 122 -103083.9599 32827.5084 123 84905.4238 -103083.9599 124 13946.7576 84905.4238 125 -162469.9198 13946.7576 126 55152.5245 -162469.9198 127 42028.6252 55152.5245 128 29692.3360 42028.6252 129 55152.5245 29692.3360 130 46453.0982 55152.5245 131 57660.3793 46453.0982 132 55152.5245 57660.3793 133 31224.0439 55152.5245 134 33610.3153 31224.0439 135 44124.6501 33610.3153 136 45697.8939 44124.6501 137 53540.1528 45697.8939 138 30900.3954 53540.1528 139 55152.5245 30900.3954 140 51353.5430 55152.5245 141 55152.5245 51353.5430 142 11405.4575 55152.5245 143 50503.8779 11405.4575 144 40171.2046 50503.8779 145 43328.0424 40171.2046 146 55152.5245 43328.0424 147 47315.9545 55152.5245 148 57036.7210 47315.9545 149 54877.0451 57036.7210 150 16156.0728 54877.0451 151 -39133.1674 16156.0728 152 50106.0545 -39133.1674 153 23675.2486 50106.0545 154 52216.7029 23675.2486 155 55152.5245 52216.7029 156 55152.5245 55152.5245 157 49624.7042 55152.5245 158 62038.6981 49624.7042 159 55152.5245 62038.6981 160 55152.5245 55152.5245 161 55152.5245 55152.5245 162 55152.5245 55152.5245 163 28169.4235 55152.5245 164 3755.9701 28169.4235 165 52610.4994 3755.9701 166 55152.5245 52610.4994 167 111226.8553 55152.5245 168 -76651.0561 111226.8553 169 55152.5245 -76651.0561 170 55152.5245 55152.5245 171 47040.3092 55152.5245 172 55152.5245 47040.3092 173 33905.3464 55152.5245 174 55152.5245 33905.3464 175 44108.9001 55152.5245 176 -69816.9214 44108.9001 177 54346.2088 -69816.9214 178 55152.5245 54346.2088 179 55152.5245 55152.5245 180 41401.1884 55152.5245 181 29918.9688 41401.1884 182 -17621.7771 29918.9688 183 -4778.0021 -17621.7771 184 46147.5877 -4778.0021 185 48383.3242 46147.5877 186 11517.2049 48383.3242 187 41696.5179 11517.2049 188 13656.0689 41696.5179 189 55152.5245 13656.0689 190 47311.6999 55152.5245 191 55152.5245 47311.6999 192 43541.9820 55152.5245 193 55152.5245 43541.9820 194 51549.5701 55152.5245 195 55152.5245 51549.5701 196 -132865.9295 55152.5245 197 55152.5245 -132865.9295 198 55152.5245 55152.5245 199 -50518.7341 55152.5245 200 55152.5245 -50518.7341 201 55152.5245 55152.5245 202 55152.5245 55152.5245 203 39706.6063 55152.5245 204 389.6725 39706.6063 205 49192.8281 389.6725 206 47148.8882 49192.8281 207 55152.5245 47148.8882 208 51266.1839 55152.5245 209 37858.5204 51266.1839 210 55152.5245 37858.5204 211 55152.5245 55152.5245 212 57681.9685 55152.5245 213 6257.6673 57681.9685 214 24474.6896 6257.6673 215 53133.9365 24474.6896 216 48900.4290 53133.9365 217 41819.3564 48900.4290 218 53133.9365 41819.3564 219 -187.5434 53133.9365 220 27389.5375 -187.5434 221 55152.5245 27389.5375 222 54085.6898 55152.5245 223 21444.3644 54085.6898 224 53765.6010 21444.3644 225 25296.9108 53765.6010 226 43896.9295 25296.9108 227 55152.5245 43896.9295 228 38720.0837 55152.5245 229 55152.5245 38720.0837 230 19417.3320 55152.5245 231 35817.3986 19417.3320 232 55152.5245 35817.3986 233 20264.8536 55152.5245 234 55152.5245 20264.8536 235 55179.4190 55152.5245 236 55152.5245 55179.4190 237 48924.8316 55152.5245 238 55152.5245 48924.8316 239 55152.5245 55152.5245 240 55152.5245 55152.5245 241 46558.3712 55152.5245 242 55152.5245 46558.3712 243 49230.8366 55152.5245 244 -7960.0214 49230.8366 245 55152.5245 -7960.0214 246 55152.5245 55152.5245 247 29031.6010 55152.5245 248 55152.5245 29031.6010 249 55152.5245 55152.5245 250 34791.5919 55152.5245 251 55152.5245 34791.5919 252 36879.8403 55152.5245 253 -27169.8902 36879.8403 254 -149946.2072 -27169.8902 255 55152.5245 -149946.2072 256 43485.6542 55152.5245 257 24396.6177 43485.6542 258 55152.5245 24396.6177 259 55152.5245 55152.5245 260 55152.5245 55152.5245 261 55152.5245 55152.5245 262 46603.3938 55152.5245 263 57799.4030 46603.3938 264 55152.5245 57799.4030 265 55152.5245 55152.5245 266 55152.5245 55152.5245 267 81168.2234 55152.5245 268 55152.5245 81168.2234 269 27644.5278 55152.5245 270 55152.5245 27644.5278 271 55152.5245 55152.5245 272 55152.5245 55152.5245 273 27922.6546 55152.5245 274 55152.5245 27922.6546 275 38752.0513 55152.5245 276 24162.5503 38752.0513 277 55152.5245 24162.5503 278 -37344.5929 55152.5245 279 26861.7964 -37344.5929 280 55152.5245 26861.7964 281 55152.5245 55152.5245 282 27804.4032 55152.5245 283 55152.5245 27804.4032 284 38853.4518 55152.5245 285 55152.5245 38853.4518 286 40704.7058 55152.5245 287 52103.9773 40704.7058 288 55152.5245 52103.9773 289 23957.2150 55152.5245 290 55152.5245 23957.2150 291 45877.7435 55152.5245 292 44577.6335 45877.7435 293 55152.5245 44577.6335 294 32825.8604 55152.5245 295 65372.1884 32825.8604 296 55152.5245 65372.1884 297 55152.5245 55152.5245 298 55152.5245 55152.5245 299 46068.8784 55152.5245 300 48116.5426 46068.8784 301 36569.8314 48116.5426 302 54143.2305 36569.8314 303 55152.5245 54143.2305 304 60303.2396 55152.5245 305 42987.4385 60303.2396 306 55152.5245 42987.4385 307 27472.0702 55152.5245 308 33465.2963 27472.0702 309 55152.5245 33465.2963 310 55152.5245 55152.5245 311 4530.3833 55152.5245 312 55152.5245 4530.3833 313 41282.3000 55152.5245 314 47223.8192 41282.3000 315 63986.5389 47223.8192 316 55180.1335 63986.5389 317 7924.8219 55180.1335 318 50077.6074 7924.8219 319 55152.5245 50077.6074 320 55152.5245 55152.5245 321 55152.5245 55152.5245 322 55152.5245 55152.5245 323 55152.5245 55152.5245 324 55152.5245 55152.5245 325 50777.7514 55152.5245 326 55152.5245 50777.7514 327 55152.5245 55152.5245 328 55152.5245 55152.5245 329 55152.5245 55152.5245 330 55152.5245 55152.5245 331 55152.5245 55152.5245 332 55152.5245 55152.5245 333 47885.6972 55152.5245 334 54118.2358 47885.6972 335 44193.4641 54118.2358 336 -3177.0170 44193.4641 337 55152.5245 -3177.0170 338 51794.7786 55152.5245 339 46172.1334 51794.7786 340 13020.8195 46172.1334 341 55152.5245 13020.8195 342 55152.5245 55152.5245 343 55152.5245 55152.5245 344 21236.0425 55152.5245 345 55152.5245 21236.0425 346 40314.4511 55152.5245 347 55152.5245 40314.4511 348 55152.5245 55152.5245 349 13933.7738 55152.5245 350 51894.0439 13933.7738 351 44801.8786 51894.0439 352 44219.9617 44801.8786 353 47891.5630 44219.9617 354 55152.5245 47891.5630 355 55152.5245 55152.5245 356 7368.8801 55152.5245 357 55152.5245 7368.8801 358 55152.5245 55152.5245 359 55152.5245 55152.5245 360 51108.1960 55152.5245 361 55152.5245 51108.1960 362 55152.5245 55152.5245 363 -20345.3694 55152.5245 364 44813.8486 -20345.3694 365 55152.5245 44813.8486 366 55152.5245 55152.5245 367 55152.5245 55152.5245 368 55152.5245 55152.5245 369 43461.0754 55152.5245 370 55152.5245 43461.0754 371 47678.2739 55152.5245 372 39372.9095 47678.2739 373 33907.2837 39372.9095 374 55152.5245 33907.2837 375 55152.5245 55152.5245 376 28500.3161 55152.5245 377 42129.4647 28500.3161 378 845.1762 42129.4647 379 55152.5245 845.1762 380 55152.5245 55152.5245 381 7288.2607 55152.5245 382 23651.2435 7288.2607 383 55152.5245 23651.2435 384 40168.2486 55152.5245 385 -117279.6312 40168.2486 386 62592.5952 -117279.6312 387 25728.1516 62592.5952 388 49874.5640 25728.1516 389 30391.7616 49874.5640 390 40145.1899 30391.7616 391 18133.7508 40145.1899 392 -19080.2469 18133.7508 393 74580.5292 -19080.2469 394 -96721.9437 74580.5292 395 6463.4051 -96721.9437 396 -138664.7201 6463.4051 397 -189519.7904 -138664.7201 398 12381.6147 -189519.7904 399 -55495.2651 12381.6147 400 2431.9986 -55495.2651 401 42613.8359 2431.9986 402 19583.9874 42613.8359 403 -49565.6943 19583.9874 404 -8995.7460 -49565.6943 405 14011.1662 -8995.7460 406 20618.4595 14011.1662 407 13250.1409 20618.4595 408 22906.3494 13250.1409 409 -61901.5139 22906.3494 410 -51361.8573 -61901.5139 411 -161769.0794 -51361.8573 412 -25907.2560 -161769.0794 413 60088.3931 -25907.2560 414 -177525.6236 60088.3931 415 -99195.2033 -177525.6236 416 91983.5806 -99195.2033 417 -17379.6380 91983.5806 418 -50409.1687 -17379.6380 419 -37404.0006 -50409.1687 420 -44265.0154 -37404.0006 421 -175393.8602 -44265.0154 422 -547870.2678 -175393.8602 423 -458560.1402 -547870.2678 424 -48844.9394 -458560.1402 425 -250014.5348 -48844.9394 426 -311858.7540 -250014.5348 427 -1135948.9209 -311858.7540 428 -272147.0098 -1135948.9209 429 -336049.2741 -272147.0098 430 -682483.1396 -336049.2741 431 NA -682483.1396 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2923814.3860 1389636.7423 [2,] 1363346.2455 2923814.3860 [3,] -2198539.3196 1363346.2455 [4,] -394043.3532 -2198539.3196 [5,] 109869.7630 -394043.3532 [6,] -135939.7551 109869.7630 [7,] 128928.3808 -135939.7551 [8,] -336334.3876 128928.3808 [9,] 715092.5993 -336334.3876 [10,] -150657.5105 715092.5993 [11,] 342307.7668 -150657.5105 [12,] 449675.2072 342307.7668 [13,] -1139741.2853 449675.2072 [14,] -214255.2259 -1139741.2853 [15,] 460574.8057 -214255.2259 [16,] -237906.6182 460574.8057 [17,] -238071.5057 -237906.6182 [18,] -295877.9417 -238071.5057 [19,] -585102.6254 -295877.9417 [20,] 203463.1027 -585102.6254 [21,] -137416.1514 203463.1027 [22,] -119601.7096 -137416.1514 [23,] -638049.8630 -119601.7096 [24,] -318324.0983 -638049.8630 [25,] -84297.5193 -318324.0983 [26,] -471896.4012 -84297.5193 [27,] -340384.4471 -471896.4012 [28,] 348299.6853 -340384.4471 [29,] -261059.4968 348299.6853 [30,] 286668.0415 -261059.4968 [31,] -87937.6349 286668.0415 [32,] 255919.3586 -87937.6349 [33,] 130738.8627 255919.3586 [34,] -54941.9881 130738.8627 [35,] -569143.8636 -54941.9881 [36,] 34310.7675 -569143.8636 [37,] 157833.6426 34310.7675 [38,] -51015.8179 157833.6426 [39,] -304565.6193 -51015.8179 [40,] -459620.0553 -304565.6193 [41,] -329011.9718 -459620.0553 [42,] 172.3323 -329011.9718 [43,] -209375.2579 172.3323 [44,] -27239.0962 -209375.2579 [45,] -53789.6599 -27239.0962 [46,] -281553.3224 -53789.6599 [47,] -264256.0037 -281553.3224 [48,] -53547.0125 -264256.0037 [49,] -181997.0873 -53547.0125 [50,] 35733.0747 -181997.0873 [51,] -86617.6159 35733.0747 [52,] -228211.1136 -86617.6159 [53,] 44146.0101 -228211.1136 [54,] -394458.1521 44146.0101 [55,] 14558.8138 -394458.1521 [56,] 320392.3041 14558.8138 [57,] -22047.0390 320392.3041 [58,] -59458.7651 -22047.0390 [59,] -205533.6443 -59458.7651 [60,] 14509.4419 -205533.6443 [61,] -179052.6824 14509.4419 [62,] 7020.6734 -179052.6824 [63,] -174141.9735 7020.6734 [64,] -15771.8893 -174141.9735 [65,] -494796.9414 -15771.8893 [66,] -65004.9696 -494796.9414 [67,] -14962.5874 -65004.9696 [68,] -100725.8195 -14962.5874 [69,] -72004.9465 -100725.8195 [70,] -122294.6437 -72004.9465 [71,] 1252.3272 -122294.6437 [72,] -88696.6011 1252.3272 [73,] 15111.9538 -88696.6011 [74,] -34736.9165 15111.9538 [75,] -98722.3509 -34736.9165 [76,] -28470.4246 -98722.3509 [77,] -73253.8756 -28470.4246 [78,] -42228.7857 -73253.8756 [79,] -37810.1071 -42228.7857 [80,] 27095.5975 -37810.1071 [81,] -49657.3672 27095.5975 [82,] -82702.7231 -49657.3672 [83,] 41679.9422 -82702.7231 [84,] 24950.3027 41679.9422 [85,] -102379.3908 24950.3027 [86,] -29329.7313 -102379.3908 [87,] 1171.4484 -29329.7313 [88,] -115198.1764 1171.4484 [89,] 32516.3322 -115198.1764 [90,] -1902.9294 32516.3322 [91,] -21807.0740 -1902.9294 [92,] 69502.3235 -21807.0740 [93,] -172301.4464 69502.3235 [94,] 49811.5996 -172301.4464 [95,] 19619.0604 49811.5996 [96,] 71821.8316 19619.0604 [97,] 2140.2336 71821.8316 [98,] -41173.8667 2140.2336 [99,] -1473.8500 -41173.8667 [100,] -35893.9515 -1473.8500 [101,] -112498.0489 -35893.9515 [102,] -165507.6930 -112498.0489 [103,] -50102.0373 -165507.6930 [104,] 30070.0947 -50102.0373 [105,] -53263.1917 30070.0947 [106,] -252162.8870 -53263.1917 [107,] 131136.6578 -252162.8870 [108,] -175232.6525 131136.6578 [109,] -92237.8474 -175232.6525 [110,] 94519.8432 -92237.8474 [111,] -190793.4409 94519.8432 [112,] -165268.7848 -190793.4409 [113,] -57466.9920 -165268.7848 [114,] 86189.5153 -57466.9920 [115,] 21656.0303 86189.5153 [116,] 17524.5950 21656.0303 [117,] 72490.9475 17524.5950 [118,] -132145.5451 72490.9475 [119,] 14988.2296 -132145.5451 [120,] -288568.0383 14988.2296 [121,] 32827.5084 -288568.0383 [122,] -103083.9599 32827.5084 [123,] 84905.4238 -103083.9599 [124,] 13946.7576 84905.4238 [125,] -162469.9198 13946.7576 [126,] 55152.5245 -162469.9198 [127,] 42028.6252 55152.5245 [128,] 29692.3360 42028.6252 [129,] 55152.5245 29692.3360 [130,] 46453.0982 55152.5245 [131,] 57660.3793 46453.0982 [132,] 55152.5245 57660.3793 [133,] 31224.0439 55152.5245 [134,] 33610.3153 31224.0439 [135,] 44124.6501 33610.3153 [136,] 45697.8939 44124.6501 [137,] 53540.1528 45697.8939 [138,] 30900.3954 53540.1528 [139,] 55152.5245 30900.3954 [140,] 51353.5430 55152.5245 [141,] 55152.5245 51353.5430 [142,] 11405.4575 55152.5245 [143,] 50503.8779 11405.4575 [144,] 40171.2046 50503.8779 [145,] 43328.0424 40171.2046 [146,] 55152.5245 43328.0424 [147,] 47315.9545 55152.5245 [148,] 57036.7210 47315.9545 [149,] 54877.0451 57036.7210 [150,] 16156.0728 54877.0451 [151,] -39133.1674 16156.0728 [152,] 50106.0545 -39133.1674 [153,] 23675.2486 50106.0545 [154,] 52216.7029 23675.2486 [155,] 55152.5245 52216.7029 [156,] 55152.5245 55152.5245 [157,] 49624.7042 55152.5245 [158,] 62038.6981 49624.7042 [159,] 55152.5245 62038.6981 [160,] 55152.5245 55152.5245 [161,] 55152.5245 55152.5245 [162,] 55152.5245 55152.5245 [163,] 28169.4235 55152.5245 [164,] 3755.9701 28169.4235 [165,] 52610.4994 3755.9701 [166,] 55152.5245 52610.4994 [167,] 111226.8553 55152.5245 [168,] -76651.0561 111226.8553 [169,] 55152.5245 -76651.0561 [170,] 55152.5245 55152.5245 [171,] 47040.3092 55152.5245 [172,] 55152.5245 47040.3092 [173,] 33905.3464 55152.5245 [174,] 55152.5245 33905.3464 [175,] 44108.9001 55152.5245 [176,] -69816.9214 44108.9001 [177,] 54346.2088 -69816.9214 [178,] 55152.5245 54346.2088 [179,] 55152.5245 55152.5245 [180,] 41401.1884 55152.5245 [181,] 29918.9688 41401.1884 [182,] -17621.7771 29918.9688 [183,] -4778.0021 -17621.7771 [184,] 46147.5877 -4778.0021 [185,] 48383.3242 46147.5877 [186,] 11517.2049 48383.3242 [187,] 41696.5179 11517.2049 [188,] 13656.0689 41696.5179 [189,] 55152.5245 13656.0689 [190,] 47311.6999 55152.5245 [191,] 55152.5245 47311.6999 [192,] 43541.9820 55152.5245 [193,] 55152.5245 43541.9820 [194,] 51549.5701 55152.5245 [195,] 55152.5245 51549.5701 [196,] -132865.9295 55152.5245 [197,] 55152.5245 -132865.9295 [198,] 55152.5245 55152.5245 [199,] -50518.7341 55152.5245 [200,] 55152.5245 -50518.7341 [201,] 55152.5245 55152.5245 [202,] 55152.5245 55152.5245 [203,] 39706.6063 55152.5245 [204,] 389.6725 39706.6063 [205,] 49192.8281 389.6725 [206,] 47148.8882 49192.8281 [207,] 55152.5245 47148.8882 [208,] 51266.1839 55152.5245 [209,] 37858.5204 51266.1839 [210,] 55152.5245 37858.5204 [211,] 55152.5245 55152.5245 [212,] 57681.9685 55152.5245 [213,] 6257.6673 57681.9685 [214,] 24474.6896 6257.6673 [215,] 53133.9365 24474.6896 [216,] 48900.4290 53133.9365 [217,] 41819.3564 48900.4290 [218,] 53133.9365 41819.3564 [219,] -187.5434 53133.9365 [220,] 27389.5375 -187.5434 [221,] 55152.5245 27389.5375 [222,] 54085.6898 55152.5245 [223,] 21444.3644 54085.6898 [224,] 53765.6010 21444.3644 [225,] 25296.9108 53765.6010 [226,] 43896.9295 25296.9108 [227,] 55152.5245 43896.9295 [228,] 38720.0837 55152.5245 [229,] 55152.5245 38720.0837 [230,] 19417.3320 55152.5245 [231,] 35817.3986 19417.3320 [232,] 55152.5245 35817.3986 [233,] 20264.8536 55152.5245 [234,] 55152.5245 20264.8536 [235,] 55179.4190 55152.5245 [236,] 55152.5245 55179.4190 [237,] 48924.8316 55152.5245 [238,] 55152.5245 48924.8316 [239,] 55152.5245 55152.5245 [240,] 55152.5245 55152.5245 [241,] 46558.3712 55152.5245 [242,] 55152.5245 46558.3712 [243,] 49230.8366 55152.5245 [244,] -7960.0214 49230.8366 [245,] 55152.5245 -7960.0214 [246,] 55152.5245 55152.5245 [247,] 29031.6010 55152.5245 [248,] 55152.5245 29031.6010 [249,] 55152.5245 55152.5245 [250,] 34791.5919 55152.5245 [251,] 55152.5245 34791.5919 [252,] 36879.8403 55152.5245 [253,] -27169.8902 36879.8403 [254,] -149946.2072 -27169.8902 [255,] 55152.5245 -149946.2072 [256,] 43485.6542 55152.5245 [257,] 24396.6177 43485.6542 [258,] 55152.5245 24396.6177 [259,] 55152.5245 55152.5245 [260,] 55152.5245 55152.5245 [261,] 55152.5245 55152.5245 [262,] 46603.3938 55152.5245 [263,] 57799.4030 46603.3938 [264,] 55152.5245 57799.4030 [265,] 55152.5245 55152.5245 [266,] 55152.5245 55152.5245 [267,] 81168.2234 55152.5245 [268,] 55152.5245 81168.2234 [269,] 27644.5278 55152.5245 [270,] 55152.5245 27644.5278 [271,] 55152.5245 55152.5245 [272,] 55152.5245 55152.5245 [273,] 27922.6546 55152.5245 [274,] 55152.5245 27922.6546 [275,] 38752.0513 55152.5245 [276,] 24162.5503 38752.0513 [277,] 55152.5245 24162.5503 [278,] -37344.5929 55152.5245 [279,] 26861.7964 -37344.5929 [280,] 55152.5245 26861.7964 [281,] 55152.5245 55152.5245 [282,] 27804.4032 55152.5245 [283,] 55152.5245 27804.4032 [284,] 38853.4518 55152.5245 [285,] 55152.5245 38853.4518 [286,] 40704.7058 55152.5245 [287,] 52103.9773 40704.7058 [288,] 55152.5245 52103.9773 [289,] 23957.2150 55152.5245 [290,] 55152.5245 23957.2150 [291,] 45877.7435 55152.5245 [292,] 44577.6335 45877.7435 [293,] 55152.5245 44577.6335 [294,] 32825.8604 55152.5245 [295,] 65372.1884 32825.8604 [296,] 55152.5245 65372.1884 [297,] 55152.5245 55152.5245 [298,] 55152.5245 55152.5245 [299,] 46068.8784 55152.5245 [300,] 48116.5426 46068.8784 [301,] 36569.8314 48116.5426 [302,] 54143.2305 36569.8314 [303,] 55152.5245 54143.2305 [304,] 60303.2396 55152.5245 [305,] 42987.4385 60303.2396 [306,] 55152.5245 42987.4385 [307,] 27472.0702 55152.5245 [308,] 33465.2963 27472.0702 [309,] 55152.5245 33465.2963 [310,] 55152.5245 55152.5245 [311,] 4530.3833 55152.5245 [312,] 55152.5245 4530.3833 [313,] 41282.3000 55152.5245 [314,] 47223.8192 41282.3000 [315,] 63986.5389 47223.8192 [316,] 55180.1335 63986.5389 [317,] 7924.8219 55180.1335 [318,] 50077.6074 7924.8219 [319,] 55152.5245 50077.6074 [320,] 55152.5245 55152.5245 [321,] 55152.5245 55152.5245 [322,] 55152.5245 55152.5245 [323,] 55152.5245 55152.5245 [324,] 55152.5245 55152.5245 [325,] 50777.7514 55152.5245 [326,] 55152.5245 50777.7514 [327,] 55152.5245 55152.5245 [328,] 55152.5245 55152.5245 [329,] 55152.5245 55152.5245 [330,] 55152.5245 55152.5245 [331,] 55152.5245 55152.5245 [332,] 55152.5245 55152.5245 [333,] 47885.6972 55152.5245 [334,] 54118.2358 47885.6972 [335,] 44193.4641 54118.2358 [336,] -3177.0170 44193.4641 [337,] 55152.5245 -3177.0170 [338,] 51794.7786 55152.5245 [339,] 46172.1334 51794.7786 [340,] 13020.8195 46172.1334 [341,] 55152.5245 13020.8195 [342,] 55152.5245 55152.5245 [343,] 55152.5245 55152.5245 [344,] 21236.0425 55152.5245 [345,] 55152.5245 21236.0425 [346,] 40314.4511 55152.5245 [347,] 55152.5245 40314.4511 [348,] 55152.5245 55152.5245 [349,] 13933.7738 55152.5245 [350,] 51894.0439 13933.7738 [351,] 44801.8786 51894.0439 [352,] 44219.9617 44801.8786 [353,] 47891.5630 44219.9617 [354,] 55152.5245 47891.5630 [355,] 55152.5245 55152.5245 [356,] 7368.8801 55152.5245 [357,] 55152.5245 7368.8801 [358,] 55152.5245 55152.5245 [359,] 55152.5245 55152.5245 [360,] 51108.1960 55152.5245 [361,] 55152.5245 51108.1960 [362,] 55152.5245 55152.5245 [363,] -20345.3694 55152.5245 [364,] 44813.8486 -20345.3694 [365,] 55152.5245 44813.8486 [366,] 55152.5245 55152.5245 [367,] 55152.5245 55152.5245 [368,] 55152.5245 55152.5245 [369,] 43461.0754 55152.5245 [370,] 55152.5245 43461.0754 [371,] 47678.2739 55152.5245 [372,] 39372.9095 47678.2739 [373,] 33907.2837 39372.9095 [374,] 55152.5245 33907.2837 [375,] 55152.5245 55152.5245 [376,] 28500.3161 55152.5245 [377,] 42129.4647 28500.3161 [378,] 845.1762 42129.4647 [379,] 55152.5245 845.1762 [380,] 55152.5245 55152.5245 [381,] 7288.2607 55152.5245 [382,] 23651.2435 7288.2607 [383,] 55152.5245 23651.2435 [384,] 40168.2486 55152.5245 [385,] -117279.6312 40168.2486 [386,] 62592.5952 -117279.6312 [387,] 25728.1516 62592.5952 [388,] 49874.5640 25728.1516 [389,] 30391.7616 49874.5640 [390,] 40145.1899 30391.7616 [391,] 18133.7508 40145.1899 [392,] -19080.2469 18133.7508 [393,] 74580.5292 -19080.2469 [394,] -96721.9437 74580.5292 [395,] 6463.4051 -96721.9437 [396,] -138664.7201 6463.4051 [397,] -189519.7904 -138664.7201 [398,] 12381.6147 -189519.7904 [399,] -55495.2651 12381.6147 [400,] 2431.9986 -55495.2651 [401,] 42613.8359 2431.9986 [402,] 19583.9874 42613.8359 [403,] -49565.6943 19583.9874 [404,] -8995.7460 -49565.6943 [405,] 14011.1662 -8995.7460 [406,] 20618.4595 14011.1662 [407,] 13250.1409 20618.4595 [408,] 22906.3494 13250.1409 [409,] -61901.5139 22906.3494 [410,] -51361.8573 -61901.5139 [411,] -161769.0794 -51361.8573 [412,] -25907.2560 -161769.0794 [413,] 60088.3931 -25907.2560 [414,] -177525.6236 60088.3931 [415,] -99195.2033 -177525.6236 [416,] 91983.5806 -99195.2033 [417,] -17379.6380 91983.5806 [418,] -50409.1687 -17379.6380 [419,] -37404.0006 -50409.1687 [420,] -44265.0154 -37404.0006 [421,] -175393.8602 -44265.0154 [422,] -547870.2678 -175393.8602 [423,] -458560.1402 -547870.2678 [424,] -48844.9394 -458560.1402 [425,] -250014.5348 -48844.9394 [426,] -311858.7540 -250014.5348 [427,] -1135948.9209 -311858.7540 [428,] -272147.0098 -1135948.9209 [429,] -336049.2741 -272147.0098 [430,] -682483.1396 -336049.2741 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2923814.3860 1389636.7423 2 1363346.2455 2923814.3860 3 -2198539.3196 1363346.2455 4 -394043.3532 -2198539.3196 5 109869.7630 -394043.3532 6 -135939.7551 109869.7630 7 128928.3808 -135939.7551 8 -336334.3876 128928.3808 9 715092.5993 -336334.3876 10 -150657.5105 715092.5993 11 342307.7668 -150657.5105 12 449675.2072 342307.7668 13 -1139741.2853 449675.2072 14 -214255.2259 -1139741.2853 15 460574.8057 -214255.2259 16 -237906.6182 460574.8057 17 -238071.5057 -237906.6182 18 -295877.9417 -238071.5057 19 -585102.6254 -295877.9417 20 203463.1027 -585102.6254 21 -137416.1514 203463.1027 22 -119601.7096 -137416.1514 23 -638049.8630 -119601.7096 24 -318324.0983 -638049.8630 25 -84297.5193 -318324.0983 26 -471896.4012 -84297.5193 27 -340384.4471 -471896.4012 28 348299.6853 -340384.4471 29 -261059.4968 348299.6853 30 286668.0415 -261059.4968 31 -87937.6349 286668.0415 32 255919.3586 -87937.6349 33 130738.8627 255919.3586 34 -54941.9881 130738.8627 35 -569143.8636 -54941.9881 36 34310.7675 -569143.8636 37 157833.6426 34310.7675 38 -51015.8179 157833.6426 39 -304565.6193 -51015.8179 40 -459620.0553 -304565.6193 41 -329011.9718 -459620.0553 42 172.3323 -329011.9718 43 -209375.2579 172.3323 44 -27239.0962 -209375.2579 45 -53789.6599 -27239.0962 46 -281553.3224 -53789.6599 47 -264256.0037 -281553.3224 48 -53547.0125 -264256.0037 49 -181997.0873 -53547.0125 50 35733.0747 -181997.0873 51 -86617.6159 35733.0747 52 -228211.1136 -86617.6159 53 44146.0101 -228211.1136 54 -394458.1521 44146.0101 55 14558.8138 -394458.1521 56 320392.3041 14558.8138 57 -22047.0390 320392.3041 58 -59458.7651 -22047.0390 59 -205533.6443 -59458.7651 60 14509.4419 -205533.6443 61 -179052.6824 14509.4419 62 7020.6734 -179052.6824 63 -174141.9735 7020.6734 64 -15771.8893 -174141.9735 65 -494796.9414 -15771.8893 66 -65004.9696 -494796.9414 67 -14962.5874 -65004.9696 68 -100725.8195 -14962.5874 69 -72004.9465 -100725.8195 70 -122294.6437 -72004.9465 71 1252.3272 -122294.6437 72 -88696.6011 1252.3272 73 15111.9538 -88696.6011 74 -34736.9165 15111.9538 75 -98722.3509 -34736.9165 76 -28470.4246 -98722.3509 77 -73253.8756 -28470.4246 78 -42228.7857 -73253.8756 79 -37810.1071 -42228.7857 80 27095.5975 -37810.1071 81 -49657.3672 27095.5975 82 -82702.7231 -49657.3672 83 41679.9422 -82702.7231 84 24950.3027 41679.9422 85 -102379.3908 24950.3027 86 -29329.7313 -102379.3908 87 1171.4484 -29329.7313 88 -115198.1764 1171.4484 89 32516.3322 -115198.1764 90 -1902.9294 32516.3322 91 -21807.0740 -1902.9294 92 69502.3235 -21807.0740 93 -172301.4464 69502.3235 94 49811.5996 -172301.4464 95 19619.0604 49811.5996 96 71821.8316 19619.0604 97 2140.2336 71821.8316 98 -41173.8667 2140.2336 99 -1473.8500 -41173.8667 100 -35893.9515 -1473.8500 101 -112498.0489 -35893.9515 102 -165507.6930 -112498.0489 103 -50102.0373 -165507.6930 104 30070.0947 -50102.0373 105 -53263.1917 30070.0947 106 -252162.8870 -53263.1917 107 131136.6578 -252162.8870 108 -175232.6525 131136.6578 109 -92237.8474 -175232.6525 110 94519.8432 -92237.8474 111 -190793.4409 94519.8432 112 -165268.7848 -190793.4409 113 -57466.9920 -165268.7848 114 86189.5153 -57466.9920 115 21656.0303 86189.5153 116 17524.5950 21656.0303 117 72490.9475 17524.5950 118 -132145.5451 72490.9475 119 14988.2296 -132145.5451 120 -288568.0383 14988.2296 121 32827.5084 -288568.0383 122 -103083.9599 32827.5084 123 84905.4238 -103083.9599 124 13946.7576 84905.4238 125 -162469.9198 13946.7576 126 55152.5245 -162469.9198 127 42028.6252 55152.5245 128 29692.3360 42028.6252 129 55152.5245 29692.3360 130 46453.0982 55152.5245 131 57660.3793 46453.0982 132 55152.5245 57660.3793 133 31224.0439 55152.5245 134 33610.3153 31224.0439 135 44124.6501 33610.3153 136 45697.8939 44124.6501 137 53540.1528 45697.8939 138 30900.3954 53540.1528 139 55152.5245 30900.3954 140 51353.5430 55152.5245 141 55152.5245 51353.5430 142 11405.4575 55152.5245 143 50503.8779 11405.4575 144 40171.2046 50503.8779 145 43328.0424 40171.2046 146 55152.5245 43328.0424 147 47315.9545 55152.5245 148 57036.7210 47315.9545 149 54877.0451 57036.7210 150 16156.0728 54877.0451 151 -39133.1674 16156.0728 152 50106.0545 -39133.1674 153 23675.2486 50106.0545 154 52216.7029 23675.2486 155 55152.5245 52216.7029 156 55152.5245 55152.5245 157 49624.7042 55152.5245 158 62038.6981 49624.7042 159 55152.5245 62038.6981 160 55152.5245 55152.5245 161 55152.5245 55152.5245 162 55152.5245 55152.5245 163 28169.4235 55152.5245 164 3755.9701 28169.4235 165 52610.4994 3755.9701 166 55152.5245 52610.4994 167 111226.8553 55152.5245 168 -76651.0561 111226.8553 169 55152.5245 -76651.0561 170 55152.5245 55152.5245 171 47040.3092 55152.5245 172 55152.5245 47040.3092 173 33905.3464 55152.5245 174 55152.5245 33905.3464 175 44108.9001 55152.5245 176 -69816.9214 44108.9001 177 54346.2088 -69816.9214 178 55152.5245 54346.2088 179 55152.5245 55152.5245 180 41401.1884 55152.5245 181 29918.9688 41401.1884 182 -17621.7771 29918.9688 183 -4778.0021 -17621.7771 184 46147.5877 -4778.0021 185 48383.3242 46147.5877 186 11517.2049 48383.3242 187 41696.5179 11517.2049 188 13656.0689 41696.5179 189 55152.5245 13656.0689 190 47311.6999 55152.5245 191 55152.5245 47311.6999 192 43541.9820 55152.5245 193 55152.5245 43541.9820 194 51549.5701 55152.5245 195 55152.5245 51549.5701 196 -132865.9295 55152.5245 197 55152.5245 -132865.9295 198 55152.5245 55152.5245 199 -50518.7341 55152.5245 200 55152.5245 -50518.7341 201 55152.5245 55152.5245 202 55152.5245 55152.5245 203 39706.6063 55152.5245 204 389.6725 39706.6063 205 49192.8281 389.6725 206 47148.8882 49192.8281 207 55152.5245 47148.8882 208 51266.1839 55152.5245 209 37858.5204 51266.1839 210 55152.5245 37858.5204 211 55152.5245 55152.5245 212 57681.9685 55152.5245 213 6257.6673 57681.9685 214 24474.6896 6257.6673 215 53133.9365 24474.6896 216 48900.4290 53133.9365 217 41819.3564 48900.4290 218 53133.9365 41819.3564 219 -187.5434 53133.9365 220 27389.5375 -187.5434 221 55152.5245 27389.5375 222 54085.6898 55152.5245 223 21444.3644 54085.6898 224 53765.6010 21444.3644 225 25296.9108 53765.6010 226 43896.9295 25296.9108 227 55152.5245 43896.9295 228 38720.0837 55152.5245 229 55152.5245 38720.0837 230 19417.3320 55152.5245 231 35817.3986 19417.3320 232 55152.5245 35817.3986 233 20264.8536 55152.5245 234 55152.5245 20264.8536 235 55179.4190 55152.5245 236 55152.5245 55179.4190 237 48924.8316 55152.5245 238 55152.5245 48924.8316 239 55152.5245 55152.5245 240 55152.5245 55152.5245 241 46558.3712 55152.5245 242 55152.5245 46558.3712 243 49230.8366 55152.5245 244 -7960.0214 49230.8366 245 55152.5245 -7960.0214 246 55152.5245 55152.5245 247 29031.6010 55152.5245 248 55152.5245 29031.6010 249 55152.5245 55152.5245 250 34791.5919 55152.5245 251 55152.5245 34791.5919 252 36879.8403 55152.5245 253 -27169.8902 36879.8403 254 -149946.2072 -27169.8902 255 55152.5245 -149946.2072 256 43485.6542 55152.5245 257 24396.6177 43485.6542 258 55152.5245 24396.6177 259 55152.5245 55152.5245 260 55152.5245 55152.5245 261 55152.5245 55152.5245 262 46603.3938 55152.5245 263 57799.4030 46603.3938 264 55152.5245 57799.4030 265 55152.5245 55152.5245 266 55152.5245 55152.5245 267 81168.2234 55152.5245 268 55152.5245 81168.2234 269 27644.5278 55152.5245 270 55152.5245 27644.5278 271 55152.5245 55152.5245 272 55152.5245 55152.5245 273 27922.6546 55152.5245 274 55152.5245 27922.6546 275 38752.0513 55152.5245 276 24162.5503 38752.0513 277 55152.5245 24162.5503 278 -37344.5929 55152.5245 279 26861.7964 -37344.5929 280 55152.5245 26861.7964 281 55152.5245 55152.5245 282 27804.4032 55152.5245 283 55152.5245 27804.4032 284 38853.4518 55152.5245 285 55152.5245 38853.4518 286 40704.7058 55152.5245 287 52103.9773 40704.7058 288 55152.5245 52103.9773 289 23957.2150 55152.5245 290 55152.5245 23957.2150 291 45877.7435 55152.5245 292 44577.6335 45877.7435 293 55152.5245 44577.6335 294 32825.8604 55152.5245 295 65372.1884 32825.8604 296 55152.5245 65372.1884 297 55152.5245 55152.5245 298 55152.5245 55152.5245 299 46068.8784 55152.5245 300 48116.5426 46068.8784 301 36569.8314 48116.5426 302 54143.2305 36569.8314 303 55152.5245 54143.2305 304 60303.2396 55152.5245 305 42987.4385 60303.2396 306 55152.5245 42987.4385 307 27472.0702 55152.5245 308 33465.2963 27472.0702 309 55152.5245 33465.2963 310 55152.5245 55152.5245 311 4530.3833 55152.5245 312 55152.5245 4530.3833 313 41282.3000 55152.5245 314 47223.8192 41282.3000 315 63986.5389 47223.8192 316 55180.1335 63986.5389 317 7924.8219 55180.1335 318 50077.6074 7924.8219 319 55152.5245 50077.6074 320 55152.5245 55152.5245 321 55152.5245 55152.5245 322 55152.5245 55152.5245 323 55152.5245 55152.5245 324 55152.5245 55152.5245 325 50777.7514 55152.5245 326 55152.5245 50777.7514 327 55152.5245 55152.5245 328 55152.5245 55152.5245 329 55152.5245 55152.5245 330 55152.5245 55152.5245 331 55152.5245 55152.5245 332 55152.5245 55152.5245 333 47885.6972 55152.5245 334 54118.2358 47885.6972 335 44193.4641 54118.2358 336 -3177.0170 44193.4641 337 55152.5245 -3177.0170 338 51794.7786 55152.5245 339 46172.1334 51794.7786 340 13020.8195 46172.1334 341 55152.5245 13020.8195 342 55152.5245 55152.5245 343 55152.5245 55152.5245 344 21236.0425 55152.5245 345 55152.5245 21236.0425 346 40314.4511 55152.5245 347 55152.5245 40314.4511 348 55152.5245 55152.5245 349 13933.7738 55152.5245 350 51894.0439 13933.7738 351 44801.8786 51894.0439 352 44219.9617 44801.8786 353 47891.5630 44219.9617 354 55152.5245 47891.5630 355 55152.5245 55152.5245 356 7368.8801 55152.5245 357 55152.5245 7368.8801 358 55152.5245 55152.5245 359 55152.5245 55152.5245 360 51108.1960 55152.5245 361 55152.5245 51108.1960 362 55152.5245 55152.5245 363 -20345.3694 55152.5245 364 44813.8486 -20345.3694 365 55152.5245 44813.8486 366 55152.5245 55152.5245 367 55152.5245 55152.5245 368 55152.5245 55152.5245 369 43461.0754 55152.5245 370 55152.5245 43461.0754 371 47678.2739 55152.5245 372 39372.9095 47678.2739 373 33907.2837 39372.9095 374 55152.5245 33907.2837 375 55152.5245 55152.5245 376 28500.3161 55152.5245 377 42129.4647 28500.3161 378 845.1762 42129.4647 379 55152.5245 845.1762 380 55152.5245 55152.5245 381 7288.2607 55152.5245 382 23651.2435 7288.2607 383 55152.5245 23651.2435 384 40168.2486 55152.5245 385 -117279.6312 40168.2486 386 62592.5952 -117279.6312 387 25728.1516 62592.5952 388 49874.5640 25728.1516 389 30391.7616 49874.5640 390 40145.1899 30391.7616 391 18133.7508 40145.1899 392 -19080.2469 18133.7508 393 74580.5292 -19080.2469 394 -96721.9437 74580.5292 395 6463.4051 -96721.9437 396 -138664.7201 6463.4051 397 -189519.7904 -138664.7201 398 12381.6147 -189519.7904 399 -55495.2651 12381.6147 400 2431.9986 -55495.2651 401 42613.8359 2431.9986 402 19583.9874 42613.8359 403 -49565.6943 19583.9874 404 -8995.7460 -49565.6943 405 14011.1662 -8995.7460 406 20618.4595 14011.1662 407 13250.1409 20618.4595 408 22906.3494 13250.1409 409 -61901.5139 22906.3494 410 -51361.8573 -61901.5139 411 -161769.0794 -51361.8573 412 -25907.2560 -161769.0794 413 60088.3931 -25907.2560 414 -177525.6236 60088.3931 415 -99195.2033 -177525.6236 416 91983.5806 -99195.2033 417 -17379.6380 91983.5806 418 -50409.1687 -17379.6380 419 -37404.0006 -50409.1687 420 -44265.0154 -37404.0006 421 -175393.8602 -44265.0154 422 -547870.2678 -175393.8602 423 -458560.1402 -547870.2678 424 -48844.9394 -458560.1402 425 -250014.5348 -48844.9394 426 -311858.7540 -250014.5348 427 -1135948.9209 -311858.7540 428 -272147.0098 -1135948.9209 429 -336049.2741 -272147.0098 430 -682483.1396 -336049.2741 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7q81b1291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/8q81b1291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/9j0ie1291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/freestat/rcomp/tmp/10j0ie1291211442.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/114izj1291211442.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/12p1x71291211442.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/13wkcj1291211442.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/14pbt41291211442.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/15auas1291211442.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/1663811291211442.tab") + } > > try(system("convert tmp/1uzlk1291211442.ps tmp/1uzlk1291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/2m8k51291211442.ps tmp/2m8k51291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/3m8k51291211442.ps tmp/3m8k51291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/4m8k51291211442.ps tmp/4m8k51291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/5fzj81291211442.ps tmp/5fzj81291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/6fzj81291211442.ps tmp/6fzj81291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/7q81b1291211442.ps tmp/7q81b1291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/8q81b1291211442.ps tmp/8q81b1291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/9j0ie1291211442.ps tmp/9j0ie1291211442.png",intern=TRUE)) character(0) > try(system("convert tmp/10j0ie1291211442.ps tmp/10j0ie1291211442.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.372 3.174 13.841