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('1' + ,'2' + ,'3' + ,'4') + ,1:431)) > y <- array(NA,dim=c(4,431),dimnames=list(c('1','2','3','4'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '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 4 1 2 3 t 1 6282154 162556 807 213118 1 2 4321023 29790 444 81767 2 3 4111912 87550 412 153198 3 4 223193 84738 428 -26007 4 5 1491348 54660 315 126942 5 6 1629616 42634 168 157214 6 7 1398893 40949 263 129352 7 8 1926517 45187 267 234817 8 9 983660 37704 228 60448 9 10 1443586 16275 129 47818 10 11 1073089 25830 104 245546 11 12 984885 12679 122 48020 12 13 1405225 18014 393 -1710 13 14 227132 43556 190 32648 14 15 929118 24811 280 95350 15 16 1071292 6575 63 151352 16 17 638830 7123 102 288170 17 18 856956 21950 265 114337 18 19 992426 37597 234 37884 19 20 444477 17821 277 122844 20 21 857217 12988 73 82340 21 22 711969 22330 67 79801 22 23 702380 13326 103 165548 23 24 358589 16189 290 116384 24 25 297978 7146 83 134028 25 26 585715 15824 56 63838 26 27 657954 27664 236 74996 27 28 209458 11920 73 31080 28 29 786690 8568 34 32168 29 30 439798 14416 139 49857 30 31 688779 3369 26 87161 31 32 574339 11819 70 106113 32 33 741409 6984 40 80570 33 34 597793 4519 42 102129 34 35 644190 2220 12 301670 35 36 377934 18562 211 102313 36 37 640273 10327 74 88577 37 38 697458 5336 80 112477 38 39 550608 2365 83 191778 39 40 207393 4069 131 79804 40 41 301607 8636 203 128294 41 42 345783 13718 56 96448 42 43 501749 4525 89 93811 43 44 379983 6869 88 117520 44 45 387475 4628 39 69159 45 46 377305 3689 25 101792 46 47 370837 4891 49 210568 47 48 430866 7489 149 136996 48 49 469107 4901 58 121920 49 50 194493 2284 41 76403 50 51 530670 3160 90 108094 51 52 518365 4150 136 134759 52 53 491303 7285 97 188873 53 54 527021 1134 63 146216 54 55 233773 4658 114 156608 55 56 405972 2384 77 61348 56 57 652925 3748 6 50350 57 58 446211 5371 47 87720 58 59 341340 1285 51 99489 59 60 387699 9327 85 87419 60 61 493408 5565 43 94355 61 62 146494 1528 32 60326 62 63 414462 3122 25 94670 63 64 364304 7561 77 82425 64 65 355178 2675 54 59017 65 66 357760 13253 251 90829 66 67 261216 880 15 80791 67 68 397144 2053 44 100423 68 69 374943 1424 73 131116 69 70 424898 4036 85 100269 70 71 202055 3045 49 27330 71 72 378525 5119 38 39039 72 73 310768 1431 35 106885 73 74 325738 554 9 79285 74 75 394510 1975 34 118881 75 76 247060 1765 20 77623 76 77 368078 1012 29 114768 77 78 236761 810 11 74015 78 79 312378 1280 52 69465 79 80 339836 666 13 117869 80 81 347385 1380 29 60982 81 82 426280 4677 66 90131 82 83 352850 876 33 138971 83 84 301881 814 15 39625 84 85 377516 514 15 102725 85 86 357312 5692 68 64239 86 87 458343 3642 100 90262 87 88 354228 540 13 103960 88 89 308636 2099 45 106611 89 90 386212 567 14 103345 90 91 393343 2001 36 95551 91 92 378509 2949 40 82903 92 93 452469 2253 68 63593 93 94 364839 6533 29 126910 94 95 358649 1889 43 37527 95 96 376641 3055 30 60247 96 97 429112 272 9 112995 97 98 330546 1414 22 70184 98 99 403560 2564 19 130140 99 100 317892 1383 9 73221 100 101 307528 1261 31 76114 101 102 235133 975 19 90534 102 103 299243 3366 55 108479 103 104 314073 576 8 113761 104 105 368186 1686 28 68696 105 106 269661 746 29 71561 106 107 125390 3192 48 59831 107 108 510834 2045 16 97890 108 109 321896 5702 47 101481 109 110 249898 1932 20 72954 110 111 408881 936 22 67939 111 112 158492 3437 33 48022 112 113 292154 5131 44 86111 113 114 289513 2397 13 74020 114 115 378049 1389 6 57530 115 116 343466 1503 35 56364 116 117 332743 402 8 84990 117 118 442882 2239 17 88590 118 119 214215 2234 11 77200 119 120 315688 837 21 61262 120 121 375195 10579 92 110309 121 122 334280 875 12 67000 122 123 355864 1585 112 93099 123 124 480382 1659 25 107577 124 125 353058 2647 17 62920 125 126 217193 3294 23 75832 126 127 315380 0 0 60720 127 128 314533 94 10 60793 128 129 318056 422 23 57935 129 130 315380 0 0 60720 130 131 314353 34 7 60630 131 132 369448 1558 25 55637 132 133 315380 0 0 60720 133 134 312846 43 20 60887 134 135 312075 645 4 60720 135 136 315009 316 4 60505 136 137 318903 115 10 60945 137 138 314887 5 1 60720 138 139 314913 897 4 60720 139 140 315380 0 0 60720 140 141 325506 389 8 58990 141 142 315380 0 0 60720 142 143 298568 1002 11 56750 143 144 315834 36 4 60894 144 145 329784 460 15 63346 145 146 312878 309 9 56535 146 147 315380 0 0 60720 147 148 314987 9 7 60835 148 149 325249 271 2 60720 149 150 315877 14 0 61016 150 151 291650 520 7 58650 151 152 305959 1766 46 60438 152 153 315380 0 5 60720 153 154 297765 458 7 58625 154 155 315245 20 2 60938 155 156 315380 0 0 60720 156 157 315380 0 0 60720 157 158 315236 98 2 61490 158 159 336425 405 5 60845 159 160 315380 0 0 60720 160 161 315380 0 0 60720 161 162 315380 0 0 60720 162 163 315380 0 0 60720 163 164 306268 483 7 60830 164 165 302187 454 24 63261 165 166 314882 47 1 60720 166 167 315380 0 0 60720 167 168 382712 757 18 45689 168 169 341570 4655 55 60720 169 170 315380 0 0 60720 170 171 315380 0 0 60720 171 172 312412 36 3 61564 172 173 315380 0 0 60720 173 174 309596 203 9 61938 174 175 315380 0 0 60720 175 176 315547 126 8 60951 176 177 313267 400 113 60720 177 178 316176 71 0 60745 178 179 315380 0 0 60720 179 180 315380 0 0 60720 180 181 359335 972 19 71642 181 182 330068 531 11 71641 182 183 314289 2461 25 55792 183 184 297413 378 16 71873 184 185 314806 23 5 62555 185 186 333210 638 11 60370 186 187 352108 2300 23 64873 187 188 313332 149 6 62041 188 189 291787 226 5 65745 189 190 315380 0 0 60720 190 191 318745 275 7 59500 191 192 315380 0 0 60720 192 193 315366 141 7 61630 193 194 315380 0 0 60720 194 195 315688 28 3 60890 195 196 315380 0 0 60720 196 197 409642 4980 89 113521 197 198 315380 0 0 60720 198 199 315380 0 0 60720 199 200 269587 472 19 80045 200 201 315380 0 0 60720 201 202 315380 0 0 60720 202 203 315380 0 0 60720 203 204 300962 203 12 50804 204 205 325479 496 12 87390 205 206 316155 10 5 61656 206 207 318574 63 2 65688 207 208 315380 0 0 60720 208 209 343613 1136 26 48522 209 210 306948 265 3 60720 210 211 315380 0 0 60720 211 212 315380 0 0 60720 212 213 330059 267 11 57640 213 214 288985 474 10 61977 214 215 304485 534 5 62620 215 216 315380 0 2 60720 216 217 315688 15 6 60831 217 218 317736 397 7 60646 218 219 315380 0 2 60720 219 220 322331 1866 28 56225 220 221 296656 288 3 60510 221 222 315380 0 0 60720 222 223 315354 3 1 60698 223 224 312161 468 20 60720 224 225 315576 20 1 60805 225 226 314922 278 22 61404 226 227 314551 61 9 60720 227 228 315380 0 0 60720 228 229 312339 192 2 65276 229 230 315380 0 0 60720 230 231 298700 317 7 63915 231 232 321376 738 9 60720 232 233 315380 0 0 60720 233 234 303230 368 13 61686 234 235 315380 0 0 60720 235 236 315487 2 0 60743 236 237 315380 0 0 60720 237 238 315793 53 6 60349 238 239 315380 0 0 60720 239 240 315380 0 0 60720 240 241 315380 0 0 60720 241 242 312887 94 3 61360 242 243 315380 0 0 60720 243 244 315637 24 7 59818 244 245 324385 2332 2 72680 245 246 315380 0 0 60720 246 247 315380 0 0 60720 247 248 308989 131 15 61808 248 249 315380 0 0 60720 249 250 315380 0 0 60720 250 251 296702 206 9 53110 251 252 315380 0 0 60720 252 253 307322 167 1 64245 253 254 304376 622 38 73007 254 255 253588 2328 57 82732 255 256 315380 0 0 60720 256 257 309560 365 7 54820 257 258 298466 364 26 47705 258 259 315380 0 0 60720 259 260 315380 0 0 60720 260 261 315380 0 0 60720 261 262 315380 0 0 60720 262 263 343929 226 13 72835 263 264 331955 307 10 58856 264 265 315380 0 0 60720 265 266 315380 0 0 60720 266 267 315380 0 0 60720 267 268 381180 188 9 77655 268 269 315380 0 0 60720 269 270 331420 138 26 69817 270 271 315380 0 0 60720 271 272 315380 0 0 60720 272 273 315380 0 0 60720 273 274 310201 125 19 60798 274 275 315380 0 0 60720 275 276 320016 282 12 62452 276 277 320398 335 23 64175 277 278 315380 0 0 60720 278 279 291841 1324 29 67440 279 280 310670 176 8 68136 280 281 315380 0 0 60720 281 282 315380 0 0 60720 282 283 313491 249 26 56726 283 284 315380 0 0 60720 284 285 331323 333 9 70811 285 286 315380 0 0 60720 286 287 319210 601 5 60720 287 288 318098 30 3 62045 288 289 315380 0 0 60720 289 290 292754 249 13 54323 290 291 315380 0 0 60720 291 292 325176 165 12 62841 292 293 365959 453 19 81125 293 294 315380 0 0 60720 294 295 302409 53 10 59506 295 296 340968 382 9 59365 296 297 315380 0 0 60720 297 298 315380 0 0 60720 298 299 315380 0 0 60720 299 300 315380 0 9 60720 300 301 313164 30 4 60798 301 302 301164 290 1 58790 302 303 315380 0 1 60720 303 304 315380 0 0 60720 304 305 344425 366 14 61808 305 306 315394 2 12 60735 306 307 315380 0 0 60720 307 308 316647 209 19 64016 308 309 309836 384 17 54683 309 310 315380 0 0 60720 310 311 315380 0 0 60720 311 312 346611 365 32 87192 312 313 315380 0 0 60720 313 314 322031 49 14 64107 314 315 315656 3 8 60761 315 316 339445 133 4 65990 316 317 314964 32 0 59988 317 318 297141 368 20 61167 318 319 315372 1 5 60719 319 320 315380 0 0 60720 320 321 315380 0 0 60720 321 322 315380 0 0 60720 322 323 315380 0 0 60720 323 324 315380 0 0 60720 324 325 315380 0 0 60720 325 326 312502 22 1 60722 326 327 315380 0 0 60720 327 328 315380 0 0 60720 328 329 315380 0 0 60720 329 330 315380 0 0 60720 330 331 315380 0 0 60720 331 332 315380 0 0 60720 332 333 315380 0 0 60720 333 334 313729 96 4 60379 334 335 315388 1 1 60727 335 336 315371 314 4 60720 336 337 296139 844 20 60925 337 338 315380 0 0 60720 338 339 313880 26 1 60896 339 340 317698 125 10 59734 340 341 295580 304 12 62969 341 342 315380 0 0 60720 342 343 315380 0 0 60720 343 344 315380 0 0 60720 344 345 308256 621 13 60720 345 346 315380 0 0 60720 346 347 303677 119 3 59118 347 348 315380 0 0 60720 348 349 315380 0 0 60720 349 350 319369 1595 10 60720 350 351 318690 312 3 58598 351 352 314049 60 7 61124 352 353 325699 587 10 59595 353 354 314210 135 1 62065 354 355 315380 0 0 60720 355 356 315380 0 0 60720 356 357 322378 514 15 78780 357 358 315380 0 0 60720 358 359 315380 0 0 60720 359 360 315380 0 0 60720 360 361 315398 1 4 60722 361 362 315380 0 0 60720 362 363 315380 0 0 60720 363 364 308336 1763 28 61600 364 365 316386 180 9 59635 365 366 315380 0 0 60720 366 367 315380 0 0 60720 367 368 315380 0 0 60720 368 369 315380 0 0 60720 369 370 315553 218 7 60720 370 371 315380 0 0 60720 371 372 323361 448 7 59781 372 373 336639 227 7 76644 373 374 307424 174 3 64820 374 375 315380 0 0 60720 375 376 315380 0 0 60720 376 377 295370 121 11 56178 377 378 322340 607 7 60436 378 379 319864 2212 10 60720 379 380 315380 0 0 60720 380 381 315380 0 0 60720 381 382 317291 530 18 73433 382 383 280398 571 14 41477 383 384 315380 0 0 60720 384 385 317330 78 12 62700 385 386 238125 2489 29 67804 386 387 327071 131 3 59661 387 388 309038 923 6 58620 388 389 314210 72 3 60398 389 390 307930 572 8 58580 390 391 322327 397 10 62710 391 392 292136 450 6 59325 392 393 263276 622 8 60950 393 394 367655 694 6 68060 394 395 283910 3425 9 83620 395 396 283587 562 8 58456 396 397 243650 4917 26 52811 397 398 438493 1442 239 121173 398 399 296261 529 7 63870 399 400 230621 2126 41 21001 400 401 304252 1061 3 70415 401 402 333505 776 8 64230 402 403 296919 611 6 59190 403 404 278990 1526 21 69351 404 405 276898 592 7 64270 405 406 327007 1182 11 70694 406 407 317046 621 11 68005 407 408 304555 989 12 58930 408 409 298096 438 9 58320 409 410 231861 726 3 69980 410 411 309422 1303 57 69863 411 412 286963 7419 21 63255 412 413 269753 1164 15 57320 413 414 448243 3310 32 75230 414 415 165404 1920 11 79420 415 416 204325 965 2 73490 416 417 407159 3256 23 35250 417 418 290476 1135 20 62285 418 419 275311 1270 24 69206 419 420 246541 661 1 65920 420 421 253468 1013 1 69770 421 422 240897 2844 74 72683 422 423 -83265 11528 68 -14545 423 424 -42143 6526 20 55830 424 425 272713 2264 20 55174 425 426 215362 5109 82 67038 426 427 42754 3999 21 51252 427 428 306275 35624 244 157278 428 429 253537 9252 32 79510 429 430 372631 15236 86 77440 430 431 -7170 18073 69 27284 431 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) `1` `2` `3` t 1.521e+05 2.194e+01 1.042e+03 1.614e+00 4.063e+01 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2192330 -24812 40310 54069 2920559 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.521e+05 4.768e+04 3.190 0.001528 ** `1` 2.194e+01 2.167e+00 10.122 < 2e-16 *** `2` 1.042e+03 3.623e+02 2.877 0.004216 ** `3` 1.614e+00 4.313e-01 3.741 0.000208 *** t 4.063e+01 1.151e+02 0.353 0.724132 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 253800 on 426 degrees of freedom Multiple R-squared: 0.6557, Adjusted R-squared: 0.6525 F-statistic: 202.9 on 4 and 426 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.049545e-12 1.024773e-12 [2,] 1.0000000 3.330072e-43 1.665036e-43 [3,] 1.0000000 1.626109e-77 8.130544e-78 [4,] 1.0000000 1.016329e-79 5.081643e-80 [5,] 1.0000000 2.684890e-88 1.342445e-88 [6,] 1.0000000 4.986150e-103 2.493075e-103 [7,] 1.0000000 3.954754e-107 1.977377e-107 [8,] 1.0000000 4.642465e-110 2.321233e-110 [9,] 1.0000000 7.044137e-122 3.522069e-122 [10,] 1.0000000 2.804502e-126 1.402251e-126 [11,] 1.0000000 2.382863e-128 1.191432e-128 [12,] 1.0000000 1.032210e-139 5.161048e-140 [13,] 1.0000000 3.557656e-141 1.778828e-141 [14,] 1.0000000 1.784468e-157 8.922340e-158 [15,] 1.0000000 7.160686e-165 3.580343e-165 [16,] 1.0000000 7.548287e-166 3.774144e-166 [17,] 1.0000000 3.907770e-167 1.953885e-167 [18,] 1.0000000 2.415530e-167 1.207765e-167 [19,] 1.0000000 1.099518e-172 5.497589e-173 [20,] 1.0000000 8.340245e-175 4.170122e-175 [21,] 1.0000000 5.046338e-176 2.523169e-176 [22,] 1.0000000 1.552772e-191 7.763858e-192 [23,] 1.0000000 1.760173e-191 8.800864e-192 [24,] 1.0000000 3.673260e-198 1.836630e-198 [25,] 1.0000000 1.422199e-199 7.110994e-200 [26,] 1.0000000 3.881845e-209 1.940923e-209 [27,] 1.0000000 1.286468e-211 6.432339e-212 [28,] 1.0000000 1.327234e-210 6.636171e-211 [29,] 1.0000000 2.135419e-210 1.067710e-210 [30,] 1.0000000 5.494673e-216 2.747336e-216 [31,] 1.0000000 3.961285e-222 1.980643e-222 [32,] 1.0000000 3.526740e-221 1.763370e-221 [33,] 1.0000000 8.043914e-222 4.021957e-222 [34,] 1.0000000 2.137789e-222 1.068895e-222 [35,] 1.0000000 1.699589e-221 8.497944e-222 [36,] 1.0000000 4.932411e-222 2.466206e-222 [37,] 1.0000000 6.345028e-221 3.172514e-221 [38,] 1.0000000 8.041461e-221 4.020731e-221 [39,] 1.0000000 3.975388e-220 1.987694e-220 [40,] 1.0000000 2.555460e-220 1.277730e-220 [41,] 1.0000000 2.889349e-219 1.444674e-219 [42,] 1.0000000 7.725807e-219 3.862904e-219 [43,] 1.0000000 9.978300e-220 4.989150e-220 [44,] 1.0000000 2.655385e-220 1.327692e-220 [45,] 1.0000000 8.671152e-220 4.335576e-220 [46,] 1.0000000 8.761674e-219 4.380837e-219 [47,] 1.0000000 1.575460e-218 7.877298e-219 [48,] 1.0000000 2.328749e-221 1.164374e-221 [49,] 1.0000000 3.402870e-221 1.701435e-221 [50,] 1.0000000 2.513238e-234 1.256619e-234 [51,] 1.0000000 1.171262e-234 5.856308e-235 [52,] 1.0000000 7.689298e-234 3.844649e-234 [53,] 1.0000000 3.111290e-233 1.555645e-233 [54,] 1.0000000 8.405666e-235 4.202833e-235 [55,] 1.0000000 6.666983e-237 3.333491e-237 [56,] 1.0000000 1.970578e-236 9.852890e-237 [57,] 1.0000000 1.501413e-235 7.507064e-236 [58,] 1.0000000 9.260166e-235 4.630083e-235 [59,] 1.0000000 1.332603e-234 6.663017e-235 [60,] 1.0000000 1.829533e-234 9.147664e-235 [61,] 1.0000000 1.564282e-233 7.821408e-234 [62,] 1.0000000 8.407455e-233 4.203728e-233 [63,] 1.0000000 5.312013e-232 2.656007e-232 [64,] 1.0000000 1.922456e-231 9.612281e-232 [65,] 1.0000000 3.788180e-232 1.894090e-232 [66,] 1.0000000 9.843170e-232 4.921585e-232 [67,] 1.0000000 6.800964e-231 3.400482e-231 [68,] 1.0000000 7.712337e-230 3.856169e-230 [69,] 1.0000000 1.039340e-229 5.196698e-230 [70,] 1.0000000 8.276236e-229 4.138118e-229 [71,] 1.0000000 5.609776e-229 2.804888e-229 [72,] 1.0000000 4.928620e-228 2.464310e-228 [73,] 1.0000000 1.491341e-227 7.456706e-228 [74,] 1.0000000 1.116697e-226 5.583485e-227 [75,] 1.0000000 2.867995e-226 1.433997e-226 [76,] 1.0000000 2.458527e-226 1.229264e-226 [77,] 1.0000000 1.757355e-225 8.786776e-226 [78,] 1.0000000 1.679467e-224 8.397333e-225 [79,] 1.0000000 7.687452e-224 3.843726e-224 [80,] 1.0000000 9.529757e-224 4.764878e-224 [81,] 1.0000000 6.501092e-223 3.250546e-223 [82,] 1.0000000 1.010119e-222 5.050594e-223 [83,] 1.0000000 1.003948e-221 5.019741e-222 [84,] 1.0000000 1.054173e-220 5.270864e-221 [85,] 1.0000000 8.916359e-220 4.458179e-220 [86,] 1.0000000 7.195645e-221 3.597822e-221 [87,] 1.0000000 6.541842e-220 3.270921e-220 [88,] 1.0000000 5.170201e-220 2.585100e-220 [89,] 1.0000000 8.063827e-220 4.031914e-220 [90,] 1.0000000 8.165304e-219 4.082652e-219 [91,] 1.0000000 9.001268e-218 4.500634e-218 [92,] 1.0000000 7.597372e-217 3.798686e-217 [93,] 1.0000000 7.753580e-216 3.876790e-216 [94,] 1.0000000 6.430493e-215 3.215246e-215 [95,] 1.0000000 2.752225e-216 1.376112e-216 [96,] 1.0000000 1.486490e-216 7.432450e-217 [97,] 1.0000000 1.517389e-217 7.586946e-218 [98,] 1.0000000 9.892778e-217 4.946389e-217 [99,] 1.0000000 2.342522e-216 1.171261e-216 [100,] 1.0000000 2.522600e-219 1.261300e-219 [101,] 1.0000000 2.341966e-220 1.170983e-220 [102,] 1.0000000 1.178566e-219 5.892829e-220 [103,] 1.0000000 1.144424e-219 5.722122e-220 [104,] 1.0000000 1.736820e-219 8.684101e-220 [105,] 1.0000000 2.925303e-220 1.462651e-220 [106,] 1.0000000 1.350341e-219 6.751704e-220 [107,] 1.0000000 7.573311e-219 3.786655e-219 [108,] 1.0000000 1.228767e-218 6.143834e-219 [109,] 1.0000000 8.561574e-218 4.280787e-218 [110,] 1.0000000 4.899369e-217 2.449684e-217 [111,] 1.0000000 8.307913e-217 4.153957e-217 [112,] 1.0000000 2.081188e-218 1.040594e-218 [113,] 1.0000000 2.424848e-217 1.212424e-217 [114,] 1.0000000 1.269171e-216 6.345855e-217 [115,] 1.0000000 1.429915e-215 7.149574e-216 [116,] 1.0000000 7.834940e-215 3.917470e-215 [117,] 1.0000000 2.189057e-214 1.094528e-214 [118,] 1.0000000 1.185418e-213 5.927091e-214 [119,] 1.0000000 8.834189e-215 4.417094e-215 [120,] 1.0000000 9.101181e-214 4.550591e-214 [121,] 1.0000000 9.682872e-213 4.841436e-213 [122,] 1.0000000 1.114808e-211 5.574041e-212 [123,] 1.0000000 1.161532e-210 5.807660e-211 [124,] 1.0000000 1.229455e-209 6.147273e-210 [125,] 1.0000000 2.558451e-209 1.279226e-209 [126,] 1.0000000 2.718409e-208 1.359205e-208 [127,] 1.0000000 2.915887e-207 1.457944e-207 [128,] 1.0000000 3.222851e-206 1.611426e-206 [129,] 1.0000000 3.534030e-205 1.767015e-205 [130,] 1.0000000 3.904286e-204 1.952143e-204 [131,] 1.0000000 4.145103e-203 2.072551e-203 [132,] 1.0000000 4.642596e-202 2.321298e-202 [133,] 1.0000000 4.925795e-201 2.462897e-201 [134,] 1.0000000 5.321324e-200 2.660662e-200 [135,] 1.0000000 5.639162e-199 2.819581e-199 [136,] 1.0000000 6.118033e-198 3.059016e-198 [137,] 1.0000000 6.494678e-197 3.247339e-197 [138,] 1.0000000 7.165212e-196 3.582606e-196 [139,] 1.0000000 7.825957e-195 3.912979e-195 [140,] 1.0000000 8.195410e-194 4.097705e-194 [141,] 1.0000000 8.558754e-193 4.279377e-193 [142,] 1.0000000 9.128375e-192 4.564188e-192 [143,] 1.0000000 9.480667e-191 4.740333e-191 [144,] 1.0000000 8.289324e-190 4.144662e-190 [145,] 1.0000000 8.239061e-189 4.119531e-189 [146,] 1.0000000 8.453780e-188 4.226890e-188 [147,] 1.0000000 7.916688e-187 3.958344e-187 [148,] 1.0000000 8.024859e-186 4.012429e-186 [149,] 1.0000000 8.120972e-185 4.060486e-185 [150,] 1.0000000 8.186254e-184 4.093127e-184 [151,] 1.0000000 8.167932e-183 4.083966e-183 [152,] 1.0000000 7.820314e-182 3.910157e-182 [153,] 1.0000000 7.804298e-181 3.902149e-181 [154,] 1.0000000 7.754825e-180 3.877412e-180 [155,] 1.0000000 7.671860e-179 3.835930e-179 [156,] 1.0000000 7.555880e-178 3.777940e-178 [157,] 1.0000000 7.044219e-177 3.522110e-177 [158,] 1.0000000 5.401004e-176 2.700502e-176 [159,] 1.0000000 5.233808e-175 2.616904e-175 [160,] 1.0000000 5.052368e-174 2.526184e-174 [161,] 1.0000000 1.885060e-174 9.425301e-175 [162,] 1.0000000 5.776604e-174 2.888302e-174 [163,] 1.0000000 5.673396e-173 2.836698e-173 [164,] 1.0000000 5.545939e-172 2.772969e-172 [165,] 1.0000000 5.227816e-171 2.613908e-171 [166,] 1.0000000 5.060671e-170 2.530336e-170 [167,] 1.0000000 4.589327e-169 2.294663e-169 [168,] 1.0000000 4.398487e-168 2.199244e-168 [169,] 1.0000000 4.190331e-167 2.095166e-167 [170,] 1.0000000 2.764189e-166 1.382094e-166 [171,] 1.0000000 2.630388e-165 1.315194e-165 [172,] 1.0000000 2.476354e-164 1.238177e-164 [173,] 1.0000000 2.320174e-163 1.160087e-163 [174,] 1.0000000 1.957354e-162 9.786769e-163 [175,] 1.0000000 1.710260e-161 8.551298e-162 [176,] 1.0000000 1.296522e-160 6.482609e-161 [177,] 1.0000000 5.275363e-160 2.637682e-160 [178,] 1.0000000 4.687470e-159 2.343735e-159 [179,] 1.0000000 3.997488e-158 1.998744e-158 [180,] 1.0000000 2.228521e-157 1.114261e-157 [181,] 1.0000000 1.982023e-156 9.910117e-157 [182,] 1.0000000 1.054603e-155 5.273013e-156 [183,] 1.0000000 9.538764e-155 4.769382e-155 [184,] 1.0000000 8.691757e-154 4.345879e-154 [185,] 1.0000000 7.790837e-153 3.895418e-153 [186,] 1.0000000 6.876207e-152 3.438104e-152 [187,] 1.0000000 6.104203e-151 3.052101e-151 [188,] 1.0000000 5.385718e-150 2.692859e-150 [189,] 1.0000000 4.735138e-149 2.367569e-149 [190,] 1.0000000 2.541282e-148 1.270641e-148 [191,] 1.0000000 2.220579e-147 1.110290e-147 [192,] 1.0000000 1.931129e-146 9.655643e-147 [193,] 1.0000000 7.854352e-147 3.927176e-147 [194,] 1.0000000 6.810993e-146 3.405496e-146 [195,] 1.0000000 5.878131e-145 2.939065e-145 [196,] 1.0000000 5.048869e-144 2.524434e-144 [197,] 1.0000000 4.380650e-143 2.190325e-143 [198,] 1.0000000 1.139059e-142 5.695293e-143 [199,] 1.0000000 9.536621e-142 4.768310e-142 [200,] 1.0000000 7.448065e-141 3.724033e-141 [201,] 1.0000000 6.251751e-140 3.125876e-140 [202,] 1.0000000 1.787956e-139 8.939778e-140 [203,] 1.0000000 1.432400e-138 7.161999e-139 [204,] 1.0000000 1.198620e-137 5.993102e-138 [205,] 1.0000000 9.982143e-137 4.991072e-137 [206,] 1.0000000 7.738719e-136 3.869359e-136 [207,] 1.0000000 4.470903e-135 2.235451e-135 [208,] 1.0000000 3.293740e-134 1.646870e-134 [209,] 1.0000000 2.692748e-133 1.346374e-133 [210,] 1.0000000 2.187244e-132 1.093622e-132 [211,] 1.0000000 1.798440e-131 8.992202e-132 [212,] 1.0000000 1.450250e-130 7.251249e-131 [213,] 1.0000000 9.423269e-130 4.711634e-130 [214,] 1.0000000 6.453251e-129 3.226626e-129 [215,] 1.0000000 5.155729e-128 2.577864e-128 [216,] 1.0000000 4.098718e-127 2.049359e-127 [217,] 1.0000000 3.181582e-126 1.590791e-126 [218,] 1.0000000 2.507079e-125 1.253540e-125 [219,] 1.0000000 1.924120e-124 9.620602e-125 [220,] 1.0000000 1.490570e-123 7.452851e-124 [221,] 1.0000000 1.158274e-122 5.791372e-123 [222,] 1.0000000 7.950795e-122 3.975397e-122 [223,] 1.0000000 6.116264e-121 3.058132e-121 [224,] 1.0000000 3.567815e-120 1.783907e-120 [225,] 1.0000000 2.746564e-119 1.373282e-119 [226,] 1.0000000 2.084118e-118 1.042059e-118 [227,] 1.0000000 1.399589e-117 6.997946e-118 [228,] 1.0000000 1.051524e-116 5.257619e-117 [229,] 1.0000000 7.865248e-116 3.932624e-116 [230,] 1.0000000 5.853286e-115 2.926643e-115 [231,] 1.0000000 4.355257e-114 2.177628e-114 [232,] 1.0000000 3.211141e-113 1.605570e-113 [233,] 1.0000000 2.356340e-112 1.178170e-112 [234,] 1.0000000 1.720876e-111 8.604378e-112 [235,] 1.0000000 1.222055e-110 6.110276e-111 [236,] 1.0000000 8.838614e-110 4.419307e-110 [237,] 1.0000000 6.411972e-109 3.205986e-109 [238,] 1.0000000 4.367455e-108 2.183728e-108 [239,] 1.0000000 3.115592e-107 1.557796e-107 [240,] 1.0000000 2.211972e-106 1.105986e-106 [241,] 1.0000000 1.450680e-105 7.253401e-106 [242,] 1.0000000 1.019867e-104 5.099335e-105 [243,] 1.0000000 7.135728e-104 3.567864e-104 [244,] 1.0000000 4.981615e-103 2.490807e-103 [245,] 1.0000000 3.453495e-102 1.726747e-102 [246,] 1.0000000 2.055096e-101 1.027548e-101 [247,] 1.0000000 7.298345e-101 3.649172e-101 [248,] 1.0000000 1.232723e-101 6.163615e-102 [249,] 1.0000000 8.441044e-101 4.220522e-101 [250,] 1.0000000 5.993491e-100 2.996746e-100 [251,] 1.0000000 4.159178e-99 2.079589e-99 [252,] 1.0000000 2.819596e-98 1.409798e-98 [253,] 1.0000000 1.902144e-97 9.510722e-98 [254,] 1.0000000 1.276952e-96 6.384760e-97 [255,] 1.0000000 8.530561e-96 4.265281e-96 [256,] 1.0000000 5.541413e-95 2.770707e-95 [257,] 1.0000000 3.707030e-94 1.853515e-94 [258,] 1.0000000 2.449153e-93 1.224576e-93 [259,] 1.0000000 1.610200e-92 8.051000e-93 [260,] 1.0000000 1.053453e-91 5.267264e-92 [261,] 1.0000000 6.940826e-91 3.470413e-91 [262,] 1.0000000 4.511618e-90 2.255809e-90 [263,] 1.0000000 2.670130e-89 1.335065e-89 [264,] 1.0000000 1.718038e-88 8.590190e-89 [265,] 1.0000000 1.100023e-87 5.500114e-88 [266,] 1.0000000 7.008676e-87 3.504338e-87 [267,] 1.0000000 4.192747e-86 2.096373e-86 [268,] 1.0000000 2.644247e-85 1.322123e-85 [269,] 1.0000000 1.659495e-84 8.297473e-85 [270,] 1.0000000 9.959005e-84 4.979502e-84 [271,] 1.0000000 6.190178e-83 3.095089e-83 [272,] 1.0000000 2.273109e-82 1.136554e-82 [273,] 1.0000000 1.045523e-81 5.227616e-82 [274,] 1.0000000 6.376613e-81 3.188306e-81 [275,] 1.0000000 3.869417e-80 1.934709e-80 [276,] 1.0000000 2.391071e-79 1.195535e-79 [277,] 1.0000000 1.438320e-78 7.191599e-79 [278,] 1.0000000 7.949586e-78 3.974793e-78 [279,] 1.0000000 4.729177e-77 2.364589e-77 [280,] 1.0000000 2.896300e-76 1.448150e-76 [281,] 1.0000000 1.696594e-75 8.482968e-76 [282,] 1.0000000 9.951494e-75 4.975747e-75 [283,] 1.0000000 5.505402e-74 2.752701e-74 [284,] 1.0000000 3.194072e-73 1.597036e-73 [285,] 1.0000000 1.879396e-72 9.396979e-73 [286,] 1.0000000 1.064516e-71 5.322580e-72 [287,] 1.0000000 6.086496e-71 3.043248e-71 [288,] 1.0000000 3.134627e-70 1.567313e-70 [289,] 1.0000000 1.765606e-69 8.828028e-70 [290,] 1.0000000 9.986554e-69 4.993277e-69 [291,] 1.0000000 5.619210e-68 2.809605e-68 [292,] 1.0000000 3.145333e-67 1.572666e-67 [293,] 1.0000000 1.734036e-66 8.670179e-67 [294,] 1.0000000 9.419122e-66 4.709561e-66 [295,] 1.0000000 4.898885e-65 2.449442e-65 [296,] 1.0000000 2.678501e-64 1.339251e-64 [297,] 1.0000000 1.458233e-63 7.291166e-64 [298,] 1.0000000 8.009358e-63 4.004679e-63 [299,] 1.0000000 4.280786e-62 2.140393e-62 [300,] 1.0000000 2.305605e-61 1.152802e-61 [301,] 1.0000000 1.153463e-60 5.767314e-61 [302,] 1.0000000 6.339271e-60 3.169635e-60 [303,] 1.0000000 3.363436e-59 1.681718e-59 [304,] 1.0000000 1.774965e-58 8.874826e-59 [305,] 1.0000000 5.502149e-58 2.751074e-58 [306,] 1.0000000 2.864765e-57 1.432383e-57 [307,] 1.0000000 1.439110e-56 7.195548e-57 [308,] 1.0000000 7.345919e-56 3.672960e-56 [309,] 1.0000000 3.925280e-55 1.962640e-55 [310,] 1.0000000 2.021098e-54 1.010549e-54 [311,] 1.0000000 8.396060e-54 4.198030e-54 [312,] 1.0000000 4.199148e-53 2.099574e-53 [313,] 1.0000000 2.101321e-52 1.050661e-52 [314,] 1.0000000 1.045648e-51 5.228238e-52 [315,] 1.0000000 5.174065e-51 2.587032e-51 [316,] 1.0000000 2.545798e-50 1.272899e-50 [317,] 1.0000000 1.245528e-49 6.227639e-50 [318,] 1.0000000 6.059150e-49 3.029575e-49 [319,] 1.0000000 2.876549e-48 1.438275e-48 [320,] 1.0000000 1.382779e-47 6.913896e-48 [321,] 1.0000000 6.608976e-47 3.304488e-47 [322,] 1.0000000 3.140556e-46 1.570278e-46 [323,] 1.0000000 1.483750e-45 7.418749e-46 [324,] 1.0000000 6.969266e-45 3.484633e-45 [325,] 1.0000000 3.254437e-44 1.627218e-44 [326,] 1.0000000 1.510835e-43 7.554175e-44 [327,] 1.0000000 6.934608e-43 3.467304e-43 [328,] 1.0000000 3.177623e-42 1.588811e-42 [329,] 1.0000000 1.457985e-41 7.289924e-42 [330,] 1.0000000 5.678135e-41 2.839067e-41 [331,] 1.0000000 2.553022e-40 1.276511e-40 [332,] 1.0000000 1.127867e-39 5.639333e-40 [333,] 1.0000000 5.092999e-39 2.546500e-39 [334,] 1.0000000 1.766471e-38 8.832354e-39 [335,] 1.0000000 7.726309e-38 3.863154e-38 [336,] 1.0000000 3.358570e-37 1.679285e-37 [337,] 1.0000000 1.450917e-36 7.254584e-37 [338,] 1.0000000 5.988279e-36 2.994140e-36 [339,] 1.0000000 2.552443e-35 1.276221e-35 [340,] 1.0000000 1.028623e-34 5.143117e-35 [341,] 1.0000000 4.323169e-34 2.161584e-34 [342,] 1.0000000 1.805438e-33 9.027189e-34 [343,] 1.0000000 7.780755e-33 3.890377e-33 [344,] 1.0000000 3.339604e-32 1.669802e-32 [345,] 1.0000000 1.349280e-31 6.746398e-32 [346,] 1.0000000 5.746229e-31 2.873115e-31 [347,] 1.0000000 2.295931e-30 1.147966e-30 [348,] 1.0000000 9.290337e-30 4.645169e-30 [349,] 1.0000000 3.734444e-29 1.867222e-29 [350,] 1.0000000 1.049958e-28 5.249791e-29 [351,] 1.0000000 4.146631e-28 2.073316e-28 [352,] 1.0000000 1.626632e-27 8.133160e-28 [353,] 1.0000000 6.337733e-27 3.168866e-27 [354,] 1.0000000 2.442982e-26 1.221491e-26 [355,] 1.0000000 9.387696e-26 4.693848e-26 [356,] 1.0000000 3.582463e-25 1.791232e-25 [357,] 1.0000000 1.314851e-24 6.574254e-25 [358,] 1.0000000 4.997395e-24 2.498698e-24 [359,] 1.0000000 1.867743e-23 9.338717e-24 [360,] 1.0000000 6.930392e-23 3.465196e-23 [361,] 1.0000000 2.552871e-22 1.276436e-22 [362,] 1.0000000 9.334479e-22 4.667240e-22 [363,] 1.0000000 3.385332e-21 1.692666e-21 [364,] 1.0000000 1.219398e-20 6.096991e-21 [365,] 1.0000000 4.454049e-20 2.227025e-20 [366,] 1.0000000 1.450759e-19 7.253795e-20 [367,] 1.0000000 4.695025e-19 2.347512e-19 [368,] 1.0000000 1.634997e-18 8.174985e-19 [369,] 1.0000000 5.647176e-18 2.823588e-18 [370,] 1.0000000 1.864363e-17 9.321815e-18 [371,] 1.0000000 6.422569e-17 3.211285e-17 [372,] 1.0000000 2.177250e-16 1.088625e-16 [373,] 1.0000000 7.284053e-16 3.642026e-16 [374,] 1.0000000 2.413466e-15 1.206733e-15 [375,] 1.0000000 7.012941e-15 3.506470e-15 [376,] 1.0000000 2.282268e-14 1.141134e-14 [377,] 1.0000000 7.340541e-14 3.670270e-14 [378,] 1.0000000 2.318502e-13 1.159251e-13 [379,] 1.0000000 3.771027e-13 1.885514e-13 [380,] 1.0000000 1.192594e-12 5.962968e-13 [381,] 1.0000000 3.706694e-12 1.853347e-12 [382,] 1.0000000 1.137015e-11 5.685077e-12 [383,] 1.0000000 3.439899e-11 1.719949e-11 [384,] 1.0000000 1.028068e-10 5.140339e-11 [385,] 1.0000000 2.956548e-10 1.478274e-10 [386,] 1.0000000 7.250624e-10 3.625312e-10 [387,] 1.0000000 1.941714e-09 9.708568e-10 [388,] 1.0000000 3.970079e-09 1.985039e-09 [389,] 1.0000000 1.069071e-08 5.345357e-09 [390,] 1.0000000 2.672203e-08 1.336102e-08 [391,] 1.0000000 4.588624e-08 2.294312e-08 [392,] 0.9999999 1.200968e-07 6.004840e-08 [393,] 0.9999998 3.221606e-07 1.610803e-07 [394,] 0.9999996 8.220961e-07 4.110480e-07 [395,] 0.9999989 2.122023e-06 1.061011e-06 [396,] 0.9999973 5.407133e-06 2.703566e-06 [397,] 0.9999940 1.205979e-05 6.029895e-06 [398,] 0.9999861 2.774546e-05 1.387273e-05 [399,] 0.9999666 6.678666e-05 3.339333e-05 [400,] 0.9999216 1.568585e-04 7.842925e-05 [401,] 0.9998208 3.584182e-04 1.792091e-04 [402,] 0.9996009 7.981064e-04 3.990532e-04 [403,] 0.9992772 1.445526e-03 7.227630e-04 [404,] 0.9985135 2.972919e-03 1.486460e-03 [405,] 0.9969994 6.001128e-03 3.000564e-03 [406,] 0.9939295 1.214096e-02 6.070480e-03 [407,] 0.9931498 1.370039e-02 6.850194e-03 [408,] 0.9904063 1.918734e-02 9.593669e-03 [409,] 0.9850379 2.992422e-02 1.496211e-02 [410,] 0.9942056 1.158886e-02 5.794429e-03 [411,] 0.9892813 2.143749e-02 1.071875e-02 [412,] 0.9783294 4.334119e-02 2.167060e-02 [413,] 0.9591700 8.166006e-02 4.083003e-02 [414,] 0.9437181 1.125637e-01 5.628187e-02 [415,] 0.8785829 2.428342e-01 1.214171e-01 [416,] 0.8364419 3.271163e-01 1.635581e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/11lt11291213940.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/21lt11291213940.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/31lt11291213940.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/4cual1291213940.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/5cual1291213940.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 1378556.5473 2920559.1318 1362238.5911 -2192329.7703 -393334.6417 6 7 8 9 10 113103.7594 -134743.8218 125506.6345 -331198.5858 722401.0795 11 12 13 14 15 -150778.8045 349484.8648 450535.3685 -1131858.9834 -213626.6275 16 17 18 19 20 464402.8104 -241544.4240 -238144.6973 -290334.6560 -586348.5220 21 22 23 24 25 210361.0732 -129541.0024 -117507.2879 -639717.8407 -314694.6713 26 27 28 29 30 -75992.6191 -469171.0226 -331535.4341 358090.3559 -255125.5876 31 32 33 34 35 293763.5421 -82550.2598 263043.2554 136595.9837 -57327.0272 36 37 38 39 40 -567882.3723 40042.8987 161870.7193 -50925.6784 -300909.0262 41 42 43 44 45 -460221.7187 -322989.6990 4491.9591 -205957.5863 -20232.1883 46 47 48 49 50 -47907.6524 -281326.6209 -263843.2642 -49691.7965 -175763.5452 51 52 53 54 55 38946.2512 -86090.2880 -228648.1861 46251.1114 -394276.9548 56 57 58 59 60 20050.5834 328781.9677 -16614.4277 -55038.8645 -201121.6032 61 62 63 64 65 19666.7573 -172340.9734 12490.8868 -169536.6759 -9760.9545 66 67 68 69 70 -495954.0026 -58910.1095 -10662.0331 -98855.4293 -68979.4098 71 72 73 74 75 -114903.1380 8593.0754 -84641.7246 21163.4329 -31231.3166 76 77 78 79 80 -92947.9122 -24768.2825 -67173.3567 -37298.2761 -33869.1040 81 82 83 84 85 33093.2582 -45987.1372 -80479.9940 48939.0844 29295.0480 86 87 88 89 90 -97691.0212 -27066.6374 5406.3288 -112059.9978 36666.8338 91 92 93 94 95 1942.7568 -17490.6982 73675.7476 -169421.7860 55877.4815 96 97 98 99 100 25133.9722 75394.2829 7264.4673 -38613.6608 3857.7440 101 102 103 104 105 -31467.5557 -108390.1822 -163256.4947 -46792.7400 34800.0371 106 107 108 109 110 -48807.2575 -247658.6891 134847.9625 -172469.1537 -87621.3880 111 112 113 114 115 99181.1773 -185446.2654 -161917.6497 -52795.5546 91719.8485 116 117 118 119 120 26251.8232 21592.4992 76197.8356 -127767.3956 19611.0134 121 122 123 124 125 -287804.9497 37409.0695 -102962.4649 87203.2981 18559.8907 126 127 128 129 130 -158629.6053 60143.8504 46653.7379 34002.5266 60021.9485 131 132 133 134 135 51057.9569 61972.4103 59900.0467 35267.8752 38193.5252 136 137 138 139 140 48652.0961 49952.0189 58051.9448 35340.0993 59615.6090 141 142 143 144 145 55620.0338 59534.3411 15639.3696 54667.5314 43853.1344 146 147 148 149 150 47463.3110 59331.1713 51218.8854 61088.6887 58921.4727 151 152 153 154 155 20074.4171 -36526.7328 53876.2080 27468.1387 55996.0625 156 157 158 159 160 58965.4658 58924.8318 53263.1106 65589.9836 58802.9300 161 162 163 164 165 58762.2960 58721.6621 58681.0281 31458.2314 6332.1673 166 167 168 169 170 55987.7124 58518.4923 114695.6319 -74826.4184 58396.5905 171 172 173 174 175 58355.9565 50068.8768 58274.6886 36650.7308 58193.4207 176 177 178 179 180 46844.7364 -70549.0715 57269.4351 58030.8849 57990.2509 181 182 183 184 185 43152.3619 31859.7547 -15320.8909 -3105.2102 48536.2790 186 187 188 189 190 50678.9229 13298.9575 43963.1090 15753.4171 57583.9114 191 192 193 194 195 49547.7785 57502.6435 45590.4319 57421.3756 53673.4067 196 197 198 199 200 57340.1077 -135660.1478 57258.8398 57218.2059 -49957.0442 201 202 203 204 205 57136.9380 57096.3040 57055.6701 41637.2367 648.6942 206 207 208 209 210 50767.8451 48603.9056 56852.5003 52706.0792 39398.4260 211 212 213 214 215 56730.5985 56689.9645 58975.7867 7363.4582 25680.0137 216 217 218 219 220 54442.9648 50033.1890 42915.7504 54321.0630 446.5766 221 222 223 224 225 28493.6951 56283.6250 55144.4390 21870.7797 54739.5314 226 227 228 229 230 25530.9287 44533.0254 56039.8213 39309.5114 55958.5534 231 232 233 234 235 19831.7371 36301.4485 55836.6515 20464.2956 55755.3836 236 237 238 239 240 55740.7560 55674.1157 49228.9269 55592.8478 55552.2139 241 242 243 244 245 55511.5799 46756.1595 55430.3120 49279.9930 -8193.6765 246 247 248 249 250 55308.4102 55267.7762 28572.8824 55186.5083 55145.8744 251 252 253 254 255 34807.2727 55064.6065 36571.6879 -29098.9220 -152852.2784 256 257 258 259 260 54902.0706 43258.1493 23824.0670 54780.1688 54739.5348 261 262 263 264 265 54698.9009 54658.2669 45109.9720 57001.9044 54536.3651 266 267 268 269 270 54495.7311 54455.0972 79382.7272 54373.8293 25568.1954 271 272 273 274 275 54292.5614 54251.9274 54211.2935 26320.8842 54130.0256 276 277 278 279 280 37236.7001 22170.3984 54008.1237 -39688.4814 25050.8307 281 282 283 284 285 53886.2219 53845.5879 25799.7149 53764.3200 36697.3672 286 287 288 289 290 53683.0521 39075.2944 50396.8199 53561.1503 22204.8590 291 292 293 294 295 53479.8824 43685.8372 41310.1326 53357.9805 30720.1589 296 297 298 299 300 63290.0063 53236.0786 53195.4447 53154.8107 43734.0894 301 302 303 304 305 45904.5503 34526.3823 51950.0430 52951.6410 57579.0697 306 307 308 309 310 40309.5053 52829.7391 24349.6867 30803.0719 52707.8373 311 312 313 314 315 52667.2033 -218.1321 52585.9354 38064.6083 44310.8329 316 317 318 319 320 60938.2202 52486.4997 4503.8971 47102.6457 52301.4978 321 322 323 324 325 52260.8638 52220.2299 52179.5959 52138.9620 52098.3280 326 327 328 329 330 47651.5540 52017.0601 51976.4262 51935.7922 51895.1583 331 332 333 334 335 51854.5243 51813.8903 51773.2564 44356.7007 50624.5211 336 337 338 339 340 40584.2542 -7323.1072 51570.0866 48132.7798 42233.0392 341 342 343 344 345 8842.5627 51407.5508 51366.9169 51326.2829 16987.8690 346 347 348 349 350 51245.0150 36348.8650 51163.7471 51123.1132 9654.7960 351 352 353 354 355 47803.9920 40406.2766 39793.7947 43575.4660 50879.3095 356 357 358 359 360 50838.6755 1743.1207 50757.4076 50716.7737 50676.1397 361 362 363 364 365 46459.4107 50594.8718 50554.2379 -25813.1803 39900.4701 366 367 368 369 370 50432.3360 50391.7021 50351.0681 50310.4342 38364.2479 371 372 373 374 375 50229.1663 42559.9748 33435.3504 28591.0951 50066.6304 376 377 378 379 380 50025.9965 23185.1999 36749.7733 -4565.5933 49863.4607 381 382 383 384 385 49822.8267 790.6458 18691.7412 49700.9249 34197.1859 386 387 388 389 390 -123899.8243 56978.0222 20080.9755 44140.9672 24572.7433 391 392 393 394 395 34019.8369 12256.4655 -25124.4773 68245.7197 -103693.0314 396 397 398 399 400 405.4308 -144772.2788 -206035.0405 5987.5583 -60992.2552 401 402 403 404 405 -4167.1290 36067.3258 13277.9864 -56796.4158 -15646.9222 406 407 408 409 410 6941.8966 13587.6916 6583.6352 16283.9769 -68871.9896 411 412 413 414 415 -60102.7556 -168604.4202 -32789.7838 51958.0975 -185299.1439 416 417 418 419 420 -106517.0924 85830.2016 -24856.5535 -58360.9935 -44536.4016 421 422 423 424 425 -51585.4272 -185152.7264 -552875.0594 -465582.0946 -56199.7327 426 427 428 429 430 -259773.3347 -319019.5228 -1152897.9866 -280631.4909 -345807.6229 431 -689241.6941 > postscript(file="/var/www/html/freestat/rcomp/tmp/653ro1291213940.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 1378556.5473 NA 1 2920559.1318 1378556.5473 2 1362238.5911 2920559.1318 3 -2192329.7703 1362238.5911 4 -393334.6417 -2192329.7703 5 113103.7594 -393334.6417 6 -134743.8218 113103.7594 7 125506.6345 -134743.8218 8 -331198.5858 125506.6345 9 722401.0795 -331198.5858 10 -150778.8045 722401.0795 11 349484.8648 -150778.8045 12 450535.3685 349484.8648 13 -1131858.9834 450535.3685 14 -213626.6275 -1131858.9834 15 464402.8104 -213626.6275 16 -241544.4240 464402.8104 17 -238144.6973 -241544.4240 18 -290334.6560 -238144.6973 19 -586348.5220 -290334.6560 20 210361.0732 -586348.5220 21 -129541.0024 210361.0732 22 -117507.2879 -129541.0024 23 -639717.8407 -117507.2879 24 -314694.6713 -639717.8407 25 -75992.6191 -314694.6713 26 -469171.0226 -75992.6191 27 -331535.4341 -469171.0226 28 358090.3559 -331535.4341 29 -255125.5876 358090.3559 30 293763.5421 -255125.5876 31 -82550.2598 293763.5421 32 263043.2554 -82550.2598 33 136595.9837 263043.2554 34 -57327.0272 136595.9837 35 -567882.3723 -57327.0272 36 40042.8987 -567882.3723 37 161870.7193 40042.8987 38 -50925.6784 161870.7193 39 -300909.0262 -50925.6784 40 -460221.7187 -300909.0262 41 -322989.6990 -460221.7187 42 4491.9591 -322989.6990 43 -205957.5863 4491.9591 44 -20232.1883 -205957.5863 45 -47907.6524 -20232.1883 46 -281326.6209 -47907.6524 47 -263843.2642 -281326.6209 48 -49691.7965 -263843.2642 49 -175763.5452 -49691.7965 50 38946.2512 -175763.5452 51 -86090.2880 38946.2512 52 -228648.1861 -86090.2880 53 46251.1114 -228648.1861 54 -394276.9548 46251.1114 55 20050.5834 -394276.9548 56 328781.9677 20050.5834 57 -16614.4277 328781.9677 58 -55038.8645 -16614.4277 59 -201121.6032 -55038.8645 60 19666.7573 -201121.6032 61 -172340.9734 19666.7573 62 12490.8868 -172340.9734 63 -169536.6759 12490.8868 64 -9760.9545 -169536.6759 65 -495954.0026 -9760.9545 66 -58910.1095 -495954.0026 67 -10662.0331 -58910.1095 68 -98855.4293 -10662.0331 69 -68979.4098 -98855.4293 70 -114903.1380 -68979.4098 71 8593.0754 -114903.1380 72 -84641.7246 8593.0754 73 21163.4329 -84641.7246 74 -31231.3166 21163.4329 75 -92947.9122 -31231.3166 76 -24768.2825 -92947.9122 77 -67173.3567 -24768.2825 78 -37298.2761 -67173.3567 79 -33869.1040 -37298.2761 80 33093.2582 -33869.1040 81 -45987.1372 33093.2582 82 -80479.9940 -45987.1372 83 48939.0844 -80479.9940 84 29295.0480 48939.0844 85 -97691.0212 29295.0480 86 -27066.6374 -97691.0212 87 5406.3288 -27066.6374 88 -112059.9978 5406.3288 89 36666.8338 -112059.9978 90 1942.7568 36666.8338 91 -17490.6982 1942.7568 92 73675.7476 -17490.6982 93 -169421.7860 73675.7476 94 55877.4815 -169421.7860 95 25133.9722 55877.4815 96 75394.2829 25133.9722 97 7264.4673 75394.2829 98 -38613.6608 7264.4673 99 3857.7440 -38613.6608 100 -31467.5557 3857.7440 101 -108390.1822 -31467.5557 102 -163256.4947 -108390.1822 103 -46792.7400 -163256.4947 104 34800.0371 -46792.7400 105 -48807.2575 34800.0371 106 -247658.6891 -48807.2575 107 134847.9625 -247658.6891 108 -172469.1537 134847.9625 109 -87621.3880 -172469.1537 110 99181.1773 -87621.3880 111 -185446.2654 99181.1773 112 -161917.6497 -185446.2654 113 -52795.5546 -161917.6497 114 91719.8485 -52795.5546 115 26251.8232 91719.8485 116 21592.4992 26251.8232 117 76197.8356 21592.4992 118 -127767.3956 76197.8356 119 19611.0134 -127767.3956 120 -287804.9497 19611.0134 121 37409.0695 -287804.9497 122 -102962.4649 37409.0695 123 87203.2981 -102962.4649 124 18559.8907 87203.2981 125 -158629.6053 18559.8907 126 60143.8504 -158629.6053 127 46653.7379 60143.8504 128 34002.5266 46653.7379 129 60021.9485 34002.5266 130 51057.9569 60021.9485 131 61972.4103 51057.9569 132 59900.0467 61972.4103 133 35267.8752 59900.0467 134 38193.5252 35267.8752 135 48652.0961 38193.5252 136 49952.0189 48652.0961 137 58051.9448 49952.0189 138 35340.0993 58051.9448 139 59615.6090 35340.0993 140 55620.0338 59615.6090 141 59534.3411 55620.0338 142 15639.3696 59534.3411 143 54667.5314 15639.3696 144 43853.1344 54667.5314 145 47463.3110 43853.1344 146 59331.1713 47463.3110 147 51218.8854 59331.1713 148 61088.6887 51218.8854 149 58921.4727 61088.6887 150 20074.4171 58921.4727 151 -36526.7328 20074.4171 152 53876.2080 -36526.7328 153 27468.1387 53876.2080 154 55996.0625 27468.1387 155 58965.4658 55996.0625 156 58924.8318 58965.4658 157 53263.1106 58924.8318 158 65589.9836 53263.1106 159 58802.9300 65589.9836 160 58762.2960 58802.9300 161 58721.6621 58762.2960 162 58681.0281 58721.6621 163 31458.2314 58681.0281 164 6332.1673 31458.2314 165 55987.7124 6332.1673 166 58518.4923 55987.7124 167 114695.6319 58518.4923 168 -74826.4184 114695.6319 169 58396.5905 -74826.4184 170 58355.9565 58396.5905 171 50068.8768 58355.9565 172 58274.6886 50068.8768 173 36650.7308 58274.6886 174 58193.4207 36650.7308 175 46844.7364 58193.4207 176 -70549.0715 46844.7364 177 57269.4351 -70549.0715 178 58030.8849 57269.4351 179 57990.2509 58030.8849 180 43152.3619 57990.2509 181 31859.7547 43152.3619 182 -15320.8909 31859.7547 183 -3105.2102 -15320.8909 184 48536.2790 -3105.2102 185 50678.9229 48536.2790 186 13298.9575 50678.9229 187 43963.1090 13298.9575 188 15753.4171 43963.1090 189 57583.9114 15753.4171 190 49547.7785 57583.9114 191 57502.6435 49547.7785 192 45590.4319 57502.6435 193 57421.3756 45590.4319 194 53673.4067 57421.3756 195 57340.1077 53673.4067 196 -135660.1478 57340.1077 197 57258.8398 -135660.1478 198 57218.2059 57258.8398 199 -49957.0442 57218.2059 200 57136.9380 -49957.0442 201 57096.3040 57136.9380 202 57055.6701 57096.3040 203 41637.2367 57055.6701 204 648.6942 41637.2367 205 50767.8451 648.6942 206 48603.9056 50767.8451 207 56852.5003 48603.9056 208 52706.0792 56852.5003 209 39398.4260 52706.0792 210 56730.5985 39398.4260 211 56689.9645 56730.5985 212 58975.7867 56689.9645 213 7363.4582 58975.7867 214 25680.0137 7363.4582 215 54442.9648 25680.0137 216 50033.1890 54442.9648 217 42915.7504 50033.1890 218 54321.0630 42915.7504 219 446.5766 54321.0630 220 28493.6951 446.5766 221 56283.6250 28493.6951 222 55144.4390 56283.6250 223 21870.7797 55144.4390 224 54739.5314 21870.7797 225 25530.9287 54739.5314 226 44533.0254 25530.9287 227 56039.8213 44533.0254 228 39309.5114 56039.8213 229 55958.5534 39309.5114 230 19831.7371 55958.5534 231 36301.4485 19831.7371 232 55836.6515 36301.4485 233 20464.2956 55836.6515 234 55755.3836 20464.2956 235 55740.7560 55755.3836 236 55674.1157 55740.7560 237 49228.9269 55674.1157 238 55592.8478 49228.9269 239 55552.2139 55592.8478 240 55511.5799 55552.2139 241 46756.1595 55511.5799 242 55430.3120 46756.1595 243 49279.9930 55430.3120 244 -8193.6765 49279.9930 245 55308.4102 -8193.6765 246 55267.7762 55308.4102 247 28572.8824 55267.7762 248 55186.5083 28572.8824 249 55145.8744 55186.5083 250 34807.2727 55145.8744 251 55064.6065 34807.2727 252 36571.6879 55064.6065 253 -29098.9220 36571.6879 254 -152852.2784 -29098.9220 255 54902.0706 -152852.2784 256 43258.1493 54902.0706 257 23824.0670 43258.1493 258 54780.1688 23824.0670 259 54739.5348 54780.1688 260 54698.9009 54739.5348 261 54658.2669 54698.9009 262 45109.9720 54658.2669 263 57001.9044 45109.9720 264 54536.3651 57001.9044 265 54495.7311 54536.3651 266 54455.0972 54495.7311 267 79382.7272 54455.0972 268 54373.8293 79382.7272 269 25568.1954 54373.8293 270 54292.5614 25568.1954 271 54251.9274 54292.5614 272 54211.2935 54251.9274 273 26320.8842 54211.2935 274 54130.0256 26320.8842 275 37236.7001 54130.0256 276 22170.3984 37236.7001 277 54008.1237 22170.3984 278 -39688.4814 54008.1237 279 25050.8307 -39688.4814 280 53886.2219 25050.8307 281 53845.5879 53886.2219 282 25799.7149 53845.5879 283 53764.3200 25799.7149 284 36697.3672 53764.3200 285 53683.0521 36697.3672 286 39075.2944 53683.0521 287 50396.8199 39075.2944 288 53561.1503 50396.8199 289 22204.8590 53561.1503 290 53479.8824 22204.8590 291 43685.8372 53479.8824 292 41310.1326 43685.8372 293 53357.9805 41310.1326 294 30720.1589 53357.9805 295 63290.0063 30720.1589 296 53236.0786 63290.0063 297 53195.4447 53236.0786 298 53154.8107 53195.4447 299 43734.0894 53154.8107 300 45904.5503 43734.0894 301 34526.3823 45904.5503 302 51950.0430 34526.3823 303 52951.6410 51950.0430 304 57579.0697 52951.6410 305 40309.5053 57579.0697 306 52829.7391 40309.5053 307 24349.6867 52829.7391 308 30803.0719 24349.6867 309 52707.8373 30803.0719 310 52667.2033 52707.8373 311 -218.1321 52667.2033 312 52585.9354 -218.1321 313 38064.6083 52585.9354 314 44310.8329 38064.6083 315 60938.2202 44310.8329 316 52486.4997 60938.2202 317 4503.8971 52486.4997 318 47102.6457 4503.8971 319 52301.4978 47102.6457 320 52260.8638 52301.4978 321 52220.2299 52260.8638 322 52179.5959 52220.2299 323 52138.9620 52179.5959 324 52098.3280 52138.9620 325 47651.5540 52098.3280 326 52017.0601 47651.5540 327 51976.4262 52017.0601 328 51935.7922 51976.4262 329 51895.1583 51935.7922 330 51854.5243 51895.1583 331 51813.8903 51854.5243 332 51773.2564 51813.8903 333 44356.7007 51773.2564 334 50624.5211 44356.7007 335 40584.2542 50624.5211 336 -7323.1072 40584.2542 337 51570.0866 -7323.1072 338 48132.7798 51570.0866 339 42233.0392 48132.7798 340 8842.5627 42233.0392 341 51407.5508 8842.5627 342 51366.9169 51407.5508 343 51326.2829 51366.9169 344 16987.8690 51326.2829 345 51245.0150 16987.8690 346 36348.8650 51245.0150 347 51163.7471 36348.8650 348 51123.1132 51163.7471 349 9654.7960 51123.1132 350 47803.9920 9654.7960 351 40406.2766 47803.9920 352 39793.7947 40406.2766 353 43575.4660 39793.7947 354 50879.3095 43575.4660 355 50838.6755 50879.3095 356 1743.1207 50838.6755 357 50757.4076 1743.1207 358 50716.7737 50757.4076 359 50676.1397 50716.7737 360 46459.4107 50676.1397 361 50594.8718 46459.4107 362 50554.2379 50594.8718 363 -25813.1803 50554.2379 364 39900.4701 -25813.1803 365 50432.3360 39900.4701 366 50391.7021 50432.3360 367 50351.0681 50391.7021 368 50310.4342 50351.0681 369 38364.2479 50310.4342 370 50229.1663 38364.2479 371 42559.9748 50229.1663 372 33435.3504 42559.9748 373 28591.0951 33435.3504 374 50066.6304 28591.0951 375 50025.9965 50066.6304 376 23185.1999 50025.9965 377 36749.7733 23185.1999 378 -4565.5933 36749.7733 379 49863.4607 -4565.5933 380 49822.8267 49863.4607 381 790.6458 49822.8267 382 18691.7412 790.6458 383 49700.9249 18691.7412 384 34197.1859 49700.9249 385 -123899.8243 34197.1859 386 56978.0222 -123899.8243 387 20080.9755 56978.0222 388 44140.9672 20080.9755 389 24572.7433 44140.9672 390 34019.8369 24572.7433 391 12256.4655 34019.8369 392 -25124.4773 12256.4655 393 68245.7197 -25124.4773 394 -103693.0314 68245.7197 395 405.4308 -103693.0314 396 -144772.2788 405.4308 397 -206035.0405 -144772.2788 398 5987.5583 -206035.0405 399 -60992.2552 5987.5583 400 -4167.1290 -60992.2552 401 36067.3258 -4167.1290 402 13277.9864 36067.3258 403 -56796.4158 13277.9864 404 -15646.9222 -56796.4158 405 6941.8966 -15646.9222 406 13587.6916 6941.8966 407 6583.6352 13587.6916 408 16283.9769 6583.6352 409 -68871.9896 16283.9769 410 -60102.7556 -68871.9896 411 -168604.4202 -60102.7556 412 -32789.7838 -168604.4202 413 51958.0975 -32789.7838 414 -185299.1439 51958.0975 415 -106517.0924 -185299.1439 416 85830.2016 -106517.0924 417 -24856.5535 85830.2016 418 -58360.9935 -24856.5535 419 -44536.4016 -58360.9935 420 -51585.4272 -44536.4016 421 -185152.7264 -51585.4272 422 -552875.0594 -185152.7264 423 -465582.0946 -552875.0594 424 -56199.7327 -465582.0946 425 -259773.3347 -56199.7327 426 -319019.5228 -259773.3347 427 -1152897.9866 -319019.5228 428 -280631.4909 -1152897.9866 429 -345807.6229 -280631.4909 430 -689241.6941 -345807.6229 431 NA -689241.6941 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2920559.1318 1378556.5473 [2,] 1362238.5911 2920559.1318 [3,] -2192329.7703 1362238.5911 [4,] -393334.6417 -2192329.7703 [5,] 113103.7594 -393334.6417 [6,] -134743.8218 113103.7594 [7,] 125506.6345 -134743.8218 [8,] -331198.5858 125506.6345 [9,] 722401.0795 -331198.5858 [10,] -150778.8045 722401.0795 [11,] 349484.8648 -150778.8045 [12,] 450535.3685 349484.8648 [13,] -1131858.9834 450535.3685 [14,] -213626.6275 -1131858.9834 [15,] 464402.8104 -213626.6275 [16,] -241544.4240 464402.8104 [17,] -238144.6973 -241544.4240 [18,] -290334.6560 -238144.6973 [19,] -586348.5220 -290334.6560 [20,] 210361.0732 -586348.5220 [21,] -129541.0024 210361.0732 [22,] -117507.2879 -129541.0024 [23,] -639717.8407 -117507.2879 [24,] -314694.6713 -639717.8407 [25,] -75992.6191 -314694.6713 [26,] -469171.0226 -75992.6191 [27,] -331535.4341 -469171.0226 [28,] 358090.3559 -331535.4341 [29,] -255125.5876 358090.3559 [30,] 293763.5421 -255125.5876 [31,] -82550.2598 293763.5421 [32,] 263043.2554 -82550.2598 [33,] 136595.9837 263043.2554 [34,] -57327.0272 136595.9837 [35,] -567882.3723 -57327.0272 [36,] 40042.8987 -567882.3723 [37,] 161870.7193 40042.8987 [38,] -50925.6784 161870.7193 [39,] -300909.0262 -50925.6784 [40,] -460221.7187 -300909.0262 [41,] -322989.6990 -460221.7187 [42,] 4491.9591 -322989.6990 [43,] -205957.5863 4491.9591 [44,] -20232.1883 -205957.5863 [45,] -47907.6524 -20232.1883 [46,] -281326.6209 -47907.6524 [47,] -263843.2642 -281326.6209 [48,] -49691.7965 -263843.2642 [49,] -175763.5452 -49691.7965 [50,] 38946.2512 -175763.5452 [51,] -86090.2880 38946.2512 [52,] -228648.1861 -86090.2880 [53,] 46251.1114 -228648.1861 [54,] -394276.9548 46251.1114 [55,] 20050.5834 -394276.9548 [56,] 328781.9677 20050.5834 [57,] -16614.4277 328781.9677 [58,] -55038.8645 -16614.4277 [59,] -201121.6032 -55038.8645 [60,] 19666.7573 -201121.6032 [61,] -172340.9734 19666.7573 [62,] 12490.8868 -172340.9734 [63,] -169536.6759 12490.8868 [64,] -9760.9545 -169536.6759 [65,] -495954.0026 -9760.9545 [66,] -58910.1095 -495954.0026 [67,] -10662.0331 -58910.1095 [68,] -98855.4293 -10662.0331 [69,] -68979.4098 -98855.4293 [70,] -114903.1380 -68979.4098 [71,] 8593.0754 -114903.1380 [72,] -84641.7246 8593.0754 [73,] 21163.4329 -84641.7246 [74,] -31231.3166 21163.4329 [75,] -92947.9122 -31231.3166 [76,] -24768.2825 -92947.9122 [77,] -67173.3567 -24768.2825 [78,] -37298.2761 -67173.3567 [79,] -33869.1040 -37298.2761 [80,] 33093.2582 -33869.1040 [81,] -45987.1372 33093.2582 [82,] -80479.9940 -45987.1372 [83,] 48939.0844 -80479.9940 [84,] 29295.0480 48939.0844 [85,] -97691.0212 29295.0480 [86,] -27066.6374 -97691.0212 [87,] 5406.3288 -27066.6374 [88,] -112059.9978 5406.3288 [89,] 36666.8338 -112059.9978 [90,] 1942.7568 36666.8338 [91,] -17490.6982 1942.7568 [92,] 73675.7476 -17490.6982 [93,] -169421.7860 73675.7476 [94,] 55877.4815 -169421.7860 [95,] 25133.9722 55877.4815 [96,] 75394.2829 25133.9722 [97,] 7264.4673 75394.2829 [98,] -38613.6608 7264.4673 [99,] 3857.7440 -38613.6608 [100,] -31467.5557 3857.7440 [101,] -108390.1822 -31467.5557 [102,] -163256.4947 -108390.1822 [103,] -46792.7400 -163256.4947 [104,] 34800.0371 -46792.7400 [105,] -48807.2575 34800.0371 [106,] -247658.6891 -48807.2575 [107,] 134847.9625 -247658.6891 [108,] -172469.1537 134847.9625 [109,] -87621.3880 -172469.1537 [110,] 99181.1773 -87621.3880 [111,] -185446.2654 99181.1773 [112,] -161917.6497 -185446.2654 [113,] -52795.5546 -161917.6497 [114,] 91719.8485 -52795.5546 [115,] 26251.8232 91719.8485 [116,] 21592.4992 26251.8232 [117,] 76197.8356 21592.4992 [118,] -127767.3956 76197.8356 [119,] 19611.0134 -127767.3956 [120,] -287804.9497 19611.0134 [121,] 37409.0695 -287804.9497 [122,] -102962.4649 37409.0695 [123,] 87203.2981 -102962.4649 [124,] 18559.8907 87203.2981 [125,] -158629.6053 18559.8907 [126,] 60143.8504 -158629.6053 [127,] 46653.7379 60143.8504 [128,] 34002.5266 46653.7379 [129,] 60021.9485 34002.5266 [130,] 51057.9569 60021.9485 [131,] 61972.4103 51057.9569 [132,] 59900.0467 61972.4103 [133,] 35267.8752 59900.0467 [134,] 38193.5252 35267.8752 [135,] 48652.0961 38193.5252 [136,] 49952.0189 48652.0961 [137,] 58051.9448 49952.0189 [138,] 35340.0993 58051.9448 [139,] 59615.6090 35340.0993 [140,] 55620.0338 59615.6090 [141,] 59534.3411 55620.0338 [142,] 15639.3696 59534.3411 [143,] 54667.5314 15639.3696 [144,] 43853.1344 54667.5314 [145,] 47463.3110 43853.1344 [146,] 59331.1713 47463.3110 [147,] 51218.8854 59331.1713 [148,] 61088.6887 51218.8854 [149,] 58921.4727 61088.6887 [150,] 20074.4171 58921.4727 [151,] -36526.7328 20074.4171 [152,] 53876.2080 -36526.7328 [153,] 27468.1387 53876.2080 [154,] 55996.0625 27468.1387 [155,] 58965.4658 55996.0625 [156,] 58924.8318 58965.4658 [157,] 53263.1106 58924.8318 [158,] 65589.9836 53263.1106 [159,] 58802.9300 65589.9836 [160,] 58762.2960 58802.9300 [161,] 58721.6621 58762.2960 [162,] 58681.0281 58721.6621 [163,] 31458.2314 58681.0281 [164,] 6332.1673 31458.2314 [165,] 55987.7124 6332.1673 [166,] 58518.4923 55987.7124 [167,] 114695.6319 58518.4923 [168,] -74826.4184 114695.6319 [169,] 58396.5905 -74826.4184 [170,] 58355.9565 58396.5905 [171,] 50068.8768 58355.9565 [172,] 58274.6886 50068.8768 [173,] 36650.7308 58274.6886 [174,] 58193.4207 36650.7308 [175,] 46844.7364 58193.4207 [176,] -70549.0715 46844.7364 [177,] 57269.4351 -70549.0715 [178,] 58030.8849 57269.4351 [179,] 57990.2509 58030.8849 [180,] 43152.3619 57990.2509 [181,] 31859.7547 43152.3619 [182,] -15320.8909 31859.7547 [183,] -3105.2102 -15320.8909 [184,] 48536.2790 -3105.2102 [185,] 50678.9229 48536.2790 [186,] 13298.9575 50678.9229 [187,] 43963.1090 13298.9575 [188,] 15753.4171 43963.1090 [189,] 57583.9114 15753.4171 [190,] 49547.7785 57583.9114 [191,] 57502.6435 49547.7785 [192,] 45590.4319 57502.6435 [193,] 57421.3756 45590.4319 [194,] 53673.4067 57421.3756 [195,] 57340.1077 53673.4067 [196,] -135660.1478 57340.1077 [197,] 57258.8398 -135660.1478 [198,] 57218.2059 57258.8398 [199,] -49957.0442 57218.2059 [200,] 57136.9380 -49957.0442 [201,] 57096.3040 57136.9380 [202,] 57055.6701 57096.3040 [203,] 41637.2367 57055.6701 [204,] 648.6942 41637.2367 [205,] 50767.8451 648.6942 [206,] 48603.9056 50767.8451 [207,] 56852.5003 48603.9056 [208,] 52706.0792 56852.5003 [209,] 39398.4260 52706.0792 [210,] 56730.5985 39398.4260 [211,] 56689.9645 56730.5985 [212,] 58975.7867 56689.9645 [213,] 7363.4582 58975.7867 [214,] 25680.0137 7363.4582 [215,] 54442.9648 25680.0137 [216,] 50033.1890 54442.9648 [217,] 42915.7504 50033.1890 [218,] 54321.0630 42915.7504 [219,] 446.5766 54321.0630 [220,] 28493.6951 446.5766 [221,] 56283.6250 28493.6951 [222,] 55144.4390 56283.6250 [223,] 21870.7797 55144.4390 [224,] 54739.5314 21870.7797 [225,] 25530.9287 54739.5314 [226,] 44533.0254 25530.9287 [227,] 56039.8213 44533.0254 [228,] 39309.5114 56039.8213 [229,] 55958.5534 39309.5114 [230,] 19831.7371 55958.5534 [231,] 36301.4485 19831.7371 [232,] 55836.6515 36301.4485 [233,] 20464.2956 55836.6515 [234,] 55755.3836 20464.2956 [235,] 55740.7560 55755.3836 [236,] 55674.1157 55740.7560 [237,] 49228.9269 55674.1157 [238,] 55592.8478 49228.9269 [239,] 55552.2139 55592.8478 [240,] 55511.5799 55552.2139 [241,] 46756.1595 55511.5799 [242,] 55430.3120 46756.1595 [243,] 49279.9930 55430.3120 [244,] -8193.6765 49279.9930 [245,] 55308.4102 -8193.6765 [246,] 55267.7762 55308.4102 [247,] 28572.8824 55267.7762 [248,] 55186.5083 28572.8824 [249,] 55145.8744 55186.5083 [250,] 34807.2727 55145.8744 [251,] 55064.6065 34807.2727 [252,] 36571.6879 55064.6065 [253,] -29098.9220 36571.6879 [254,] -152852.2784 -29098.9220 [255,] 54902.0706 -152852.2784 [256,] 43258.1493 54902.0706 [257,] 23824.0670 43258.1493 [258,] 54780.1688 23824.0670 [259,] 54739.5348 54780.1688 [260,] 54698.9009 54739.5348 [261,] 54658.2669 54698.9009 [262,] 45109.9720 54658.2669 [263,] 57001.9044 45109.9720 [264,] 54536.3651 57001.9044 [265,] 54495.7311 54536.3651 [266,] 54455.0972 54495.7311 [267,] 79382.7272 54455.0972 [268,] 54373.8293 79382.7272 [269,] 25568.1954 54373.8293 [270,] 54292.5614 25568.1954 [271,] 54251.9274 54292.5614 [272,] 54211.2935 54251.9274 [273,] 26320.8842 54211.2935 [274,] 54130.0256 26320.8842 [275,] 37236.7001 54130.0256 [276,] 22170.3984 37236.7001 [277,] 54008.1237 22170.3984 [278,] -39688.4814 54008.1237 [279,] 25050.8307 -39688.4814 [280,] 53886.2219 25050.8307 [281,] 53845.5879 53886.2219 [282,] 25799.7149 53845.5879 [283,] 53764.3200 25799.7149 [284,] 36697.3672 53764.3200 [285,] 53683.0521 36697.3672 [286,] 39075.2944 53683.0521 [287,] 50396.8199 39075.2944 [288,] 53561.1503 50396.8199 [289,] 22204.8590 53561.1503 [290,] 53479.8824 22204.8590 [291,] 43685.8372 53479.8824 [292,] 41310.1326 43685.8372 [293,] 53357.9805 41310.1326 [294,] 30720.1589 53357.9805 [295,] 63290.0063 30720.1589 [296,] 53236.0786 63290.0063 [297,] 53195.4447 53236.0786 [298,] 53154.8107 53195.4447 [299,] 43734.0894 53154.8107 [300,] 45904.5503 43734.0894 [301,] 34526.3823 45904.5503 [302,] 51950.0430 34526.3823 [303,] 52951.6410 51950.0430 [304,] 57579.0697 52951.6410 [305,] 40309.5053 57579.0697 [306,] 52829.7391 40309.5053 [307,] 24349.6867 52829.7391 [308,] 30803.0719 24349.6867 [309,] 52707.8373 30803.0719 [310,] 52667.2033 52707.8373 [311,] -218.1321 52667.2033 [312,] 52585.9354 -218.1321 [313,] 38064.6083 52585.9354 [314,] 44310.8329 38064.6083 [315,] 60938.2202 44310.8329 [316,] 52486.4997 60938.2202 [317,] 4503.8971 52486.4997 [318,] 47102.6457 4503.8971 [319,] 52301.4978 47102.6457 [320,] 52260.8638 52301.4978 [321,] 52220.2299 52260.8638 [322,] 52179.5959 52220.2299 [323,] 52138.9620 52179.5959 [324,] 52098.3280 52138.9620 [325,] 47651.5540 52098.3280 [326,] 52017.0601 47651.5540 [327,] 51976.4262 52017.0601 [328,] 51935.7922 51976.4262 [329,] 51895.1583 51935.7922 [330,] 51854.5243 51895.1583 [331,] 51813.8903 51854.5243 [332,] 51773.2564 51813.8903 [333,] 44356.7007 51773.2564 [334,] 50624.5211 44356.7007 [335,] 40584.2542 50624.5211 [336,] -7323.1072 40584.2542 [337,] 51570.0866 -7323.1072 [338,] 48132.7798 51570.0866 [339,] 42233.0392 48132.7798 [340,] 8842.5627 42233.0392 [341,] 51407.5508 8842.5627 [342,] 51366.9169 51407.5508 [343,] 51326.2829 51366.9169 [344,] 16987.8690 51326.2829 [345,] 51245.0150 16987.8690 [346,] 36348.8650 51245.0150 [347,] 51163.7471 36348.8650 [348,] 51123.1132 51163.7471 [349,] 9654.7960 51123.1132 [350,] 47803.9920 9654.7960 [351,] 40406.2766 47803.9920 [352,] 39793.7947 40406.2766 [353,] 43575.4660 39793.7947 [354,] 50879.3095 43575.4660 [355,] 50838.6755 50879.3095 [356,] 1743.1207 50838.6755 [357,] 50757.4076 1743.1207 [358,] 50716.7737 50757.4076 [359,] 50676.1397 50716.7737 [360,] 46459.4107 50676.1397 [361,] 50594.8718 46459.4107 [362,] 50554.2379 50594.8718 [363,] -25813.1803 50554.2379 [364,] 39900.4701 -25813.1803 [365,] 50432.3360 39900.4701 [366,] 50391.7021 50432.3360 [367,] 50351.0681 50391.7021 [368,] 50310.4342 50351.0681 [369,] 38364.2479 50310.4342 [370,] 50229.1663 38364.2479 [371,] 42559.9748 50229.1663 [372,] 33435.3504 42559.9748 [373,] 28591.0951 33435.3504 [374,] 50066.6304 28591.0951 [375,] 50025.9965 50066.6304 [376,] 23185.1999 50025.9965 [377,] 36749.7733 23185.1999 [378,] -4565.5933 36749.7733 [379,] 49863.4607 -4565.5933 [380,] 49822.8267 49863.4607 [381,] 790.6458 49822.8267 [382,] 18691.7412 790.6458 [383,] 49700.9249 18691.7412 [384,] 34197.1859 49700.9249 [385,] -123899.8243 34197.1859 [386,] 56978.0222 -123899.8243 [387,] 20080.9755 56978.0222 [388,] 44140.9672 20080.9755 [389,] 24572.7433 44140.9672 [390,] 34019.8369 24572.7433 [391,] 12256.4655 34019.8369 [392,] -25124.4773 12256.4655 [393,] 68245.7197 -25124.4773 [394,] -103693.0314 68245.7197 [395,] 405.4308 -103693.0314 [396,] -144772.2788 405.4308 [397,] -206035.0405 -144772.2788 [398,] 5987.5583 -206035.0405 [399,] -60992.2552 5987.5583 [400,] -4167.1290 -60992.2552 [401,] 36067.3258 -4167.1290 [402,] 13277.9864 36067.3258 [403,] -56796.4158 13277.9864 [404,] -15646.9222 -56796.4158 [405,] 6941.8966 -15646.9222 [406,] 13587.6916 6941.8966 [407,] 6583.6352 13587.6916 [408,] 16283.9769 6583.6352 [409,] -68871.9896 16283.9769 [410,] -60102.7556 -68871.9896 [411,] -168604.4202 -60102.7556 [412,] -32789.7838 -168604.4202 [413,] 51958.0975 -32789.7838 [414,] -185299.1439 51958.0975 [415,] -106517.0924 -185299.1439 [416,] 85830.2016 -106517.0924 [417,] -24856.5535 85830.2016 [418,] -58360.9935 -24856.5535 [419,] -44536.4016 -58360.9935 [420,] -51585.4272 -44536.4016 [421,] -185152.7264 -51585.4272 [422,] -552875.0594 -185152.7264 [423,] -465582.0946 -552875.0594 [424,] -56199.7327 -465582.0946 [425,] -259773.3347 -56199.7327 [426,] -319019.5228 -259773.3347 [427,] -1152897.9866 -319019.5228 [428,] -280631.4909 -1152897.9866 [429,] -345807.6229 -280631.4909 [430,] -689241.6941 -345807.6229 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2920559.1318 1378556.5473 2 1362238.5911 2920559.1318 3 -2192329.7703 1362238.5911 4 -393334.6417 -2192329.7703 5 113103.7594 -393334.6417 6 -134743.8218 113103.7594 7 125506.6345 -134743.8218 8 -331198.5858 125506.6345 9 722401.0795 -331198.5858 10 -150778.8045 722401.0795 11 349484.8648 -150778.8045 12 450535.3685 349484.8648 13 -1131858.9834 450535.3685 14 -213626.6275 -1131858.9834 15 464402.8104 -213626.6275 16 -241544.4240 464402.8104 17 -238144.6973 -241544.4240 18 -290334.6560 -238144.6973 19 -586348.5220 -290334.6560 20 210361.0732 -586348.5220 21 -129541.0024 210361.0732 22 -117507.2879 -129541.0024 23 -639717.8407 -117507.2879 24 -314694.6713 -639717.8407 25 -75992.6191 -314694.6713 26 -469171.0226 -75992.6191 27 -331535.4341 -469171.0226 28 358090.3559 -331535.4341 29 -255125.5876 358090.3559 30 293763.5421 -255125.5876 31 -82550.2598 293763.5421 32 263043.2554 -82550.2598 33 136595.9837 263043.2554 34 -57327.0272 136595.9837 35 -567882.3723 -57327.0272 36 40042.8987 -567882.3723 37 161870.7193 40042.8987 38 -50925.6784 161870.7193 39 -300909.0262 -50925.6784 40 -460221.7187 -300909.0262 41 -322989.6990 -460221.7187 42 4491.9591 -322989.6990 43 -205957.5863 4491.9591 44 -20232.1883 -205957.5863 45 -47907.6524 -20232.1883 46 -281326.6209 -47907.6524 47 -263843.2642 -281326.6209 48 -49691.7965 -263843.2642 49 -175763.5452 -49691.7965 50 38946.2512 -175763.5452 51 -86090.2880 38946.2512 52 -228648.1861 -86090.2880 53 46251.1114 -228648.1861 54 -394276.9548 46251.1114 55 20050.5834 -394276.9548 56 328781.9677 20050.5834 57 -16614.4277 328781.9677 58 -55038.8645 -16614.4277 59 -201121.6032 -55038.8645 60 19666.7573 -201121.6032 61 -172340.9734 19666.7573 62 12490.8868 -172340.9734 63 -169536.6759 12490.8868 64 -9760.9545 -169536.6759 65 -495954.0026 -9760.9545 66 -58910.1095 -495954.0026 67 -10662.0331 -58910.1095 68 -98855.4293 -10662.0331 69 -68979.4098 -98855.4293 70 -114903.1380 -68979.4098 71 8593.0754 -114903.1380 72 -84641.7246 8593.0754 73 21163.4329 -84641.7246 74 -31231.3166 21163.4329 75 -92947.9122 -31231.3166 76 -24768.2825 -92947.9122 77 -67173.3567 -24768.2825 78 -37298.2761 -67173.3567 79 -33869.1040 -37298.2761 80 33093.2582 -33869.1040 81 -45987.1372 33093.2582 82 -80479.9940 -45987.1372 83 48939.0844 -80479.9940 84 29295.0480 48939.0844 85 -97691.0212 29295.0480 86 -27066.6374 -97691.0212 87 5406.3288 -27066.6374 88 -112059.9978 5406.3288 89 36666.8338 -112059.9978 90 1942.7568 36666.8338 91 -17490.6982 1942.7568 92 73675.7476 -17490.6982 93 -169421.7860 73675.7476 94 55877.4815 -169421.7860 95 25133.9722 55877.4815 96 75394.2829 25133.9722 97 7264.4673 75394.2829 98 -38613.6608 7264.4673 99 3857.7440 -38613.6608 100 -31467.5557 3857.7440 101 -108390.1822 -31467.5557 102 -163256.4947 -108390.1822 103 -46792.7400 -163256.4947 104 34800.0371 -46792.7400 105 -48807.2575 34800.0371 106 -247658.6891 -48807.2575 107 134847.9625 -247658.6891 108 -172469.1537 134847.9625 109 -87621.3880 -172469.1537 110 99181.1773 -87621.3880 111 -185446.2654 99181.1773 112 -161917.6497 -185446.2654 113 -52795.5546 -161917.6497 114 91719.8485 -52795.5546 115 26251.8232 91719.8485 116 21592.4992 26251.8232 117 76197.8356 21592.4992 118 -127767.3956 76197.8356 119 19611.0134 -127767.3956 120 -287804.9497 19611.0134 121 37409.0695 -287804.9497 122 -102962.4649 37409.0695 123 87203.2981 -102962.4649 124 18559.8907 87203.2981 125 -158629.6053 18559.8907 126 60143.8504 -158629.6053 127 46653.7379 60143.8504 128 34002.5266 46653.7379 129 60021.9485 34002.5266 130 51057.9569 60021.9485 131 61972.4103 51057.9569 132 59900.0467 61972.4103 133 35267.8752 59900.0467 134 38193.5252 35267.8752 135 48652.0961 38193.5252 136 49952.0189 48652.0961 137 58051.9448 49952.0189 138 35340.0993 58051.9448 139 59615.6090 35340.0993 140 55620.0338 59615.6090 141 59534.3411 55620.0338 142 15639.3696 59534.3411 143 54667.5314 15639.3696 144 43853.1344 54667.5314 145 47463.3110 43853.1344 146 59331.1713 47463.3110 147 51218.8854 59331.1713 148 61088.6887 51218.8854 149 58921.4727 61088.6887 150 20074.4171 58921.4727 151 -36526.7328 20074.4171 152 53876.2080 -36526.7328 153 27468.1387 53876.2080 154 55996.0625 27468.1387 155 58965.4658 55996.0625 156 58924.8318 58965.4658 157 53263.1106 58924.8318 158 65589.9836 53263.1106 159 58802.9300 65589.9836 160 58762.2960 58802.9300 161 58721.6621 58762.2960 162 58681.0281 58721.6621 163 31458.2314 58681.0281 164 6332.1673 31458.2314 165 55987.7124 6332.1673 166 58518.4923 55987.7124 167 114695.6319 58518.4923 168 -74826.4184 114695.6319 169 58396.5905 -74826.4184 170 58355.9565 58396.5905 171 50068.8768 58355.9565 172 58274.6886 50068.8768 173 36650.7308 58274.6886 174 58193.4207 36650.7308 175 46844.7364 58193.4207 176 -70549.0715 46844.7364 177 57269.4351 -70549.0715 178 58030.8849 57269.4351 179 57990.2509 58030.8849 180 43152.3619 57990.2509 181 31859.7547 43152.3619 182 -15320.8909 31859.7547 183 -3105.2102 -15320.8909 184 48536.2790 -3105.2102 185 50678.9229 48536.2790 186 13298.9575 50678.9229 187 43963.1090 13298.9575 188 15753.4171 43963.1090 189 57583.9114 15753.4171 190 49547.7785 57583.9114 191 57502.6435 49547.7785 192 45590.4319 57502.6435 193 57421.3756 45590.4319 194 53673.4067 57421.3756 195 57340.1077 53673.4067 196 -135660.1478 57340.1077 197 57258.8398 -135660.1478 198 57218.2059 57258.8398 199 -49957.0442 57218.2059 200 57136.9380 -49957.0442 201 57096.3040 57136.9380 202 57055.6701 57096.3040 203 41637.2367 57055.6701 204 648.6942 41637.2367 205 50767.8451 648.6942 206 48603.9056 50767.8451 207 56852.5003 48603.9056 208 52706.0792 56852.5003 209 39398.4260 52706.0792 210 56730.5985 39398.4260 211 56689.9645 56730.5985 212 58975.7867 56689.9645 213 7363.4582 58975.7867 214 25680.0137 7363.4582 215 54442.9648 25680.0137 216 50033.1890 54442.9648 217 42915.7504 50033.1890 218 54321.0630 42915.7504 219 446.5766 54321.0630 220 28493.6951 446.5766 221 56283.6250 28493.6951 222 55144.4390 56283.6250 223 21870.7797 55144.4390 224 54739.5314 21870.7797 225 25530.9287 54739.5314 226 44533.0254 25530.9287 227 56039.8213 44533.0254 228 39309.5114 56039.8213 229 55958.5534 39309.5114 230 19831.7371 55958.5534 231 36301.4485 19831.7371 232 55836.6515 36301.4485 233 20464.2956 55836.6515 234 55755.3836 20464.2956 235 55740.7560 55755.3836 236 55674.1157 55740.7560 237 49228.9269 55674.1157 238 55592.8478 49228.9269 239 55552.2139 55592.8478 240 55511.5799 55552.2139 241 46756.1595 55511.5799 242 55430.3120 46756.1595 243 49279.9930 55430.3120 244 -8193.6765 49279.9930 245 55308.4102 -8193.6765 246 55267.7762 55308.4102 247 28572.8824 55267.7762 248 55186.5083 28572.8824 249 55145.8744 55186.5083 250 34807.2727 55145.8744 251 55064.6065 34807.2727 252 36571.6879 55064.6065 253 -29098.9220 36571.6879 254 -152852.2784 -29098.9220 255 54902.0706 -152852.2784 256 43258.1493 54902.0706 257 23824.0670 43258.1493 258 54780.1688 23824.0670 259 54739.5348 54780.1688 260 54698.9009 54739.5348 261 54658.2669 54698.9009 262 45109.9720 54658.2669 263 57001.9044 45109.9720 264 54536.3651 57001.9044 265 54495.7311 54536.3651 266 54455.0972 54495.7311 267 79382.7272 54455.0972 268 54373.8293 79382.7272 269 25568.1954 54373.8293 270 54292.5614 25568.1954 271 54251.9274 54292.5614 272 54211.2935 54251.9274 273 26320.8842 54211.2935 274 54130.0256 26320.8842 275 37236.7001 54130.0256 276 22170.3984 37236.7001 277 54008.1237 22170.3984 278 -39688.4814 54008.1237 279 25050.8307 -39688.4814 280 53886.2219 25050.8307 281 53845.5879 53886.2219 282 25799.7149 53845.5879 283 53764.3200 25799.7149 284 36697.3672 53764.3200 285 53683.0521 36697.3672 286 39075.2944 53683.0521 287 50396.8199 39075.2944 288 53561.1503 50396.8199 289 22204.8590 53561.1503 290 53479.8824 22204.8590 291 43685.8372 53479.8824 292 41310.1326 43685.8372 293 53357.9805 41310.1326 294 30720.1589 53357.9805 295 63290.0063 30720.1589 296 53236.0786 63290.0063 297 53195.4447 53236.0786 298 53154.8107 53195.4447 299 43734.0894 53154.8107 300 45904.5503 43734.0894 301 34526.3823 45904.5503 302 51950.0430 34526.3823 303 52951.6410 51950.0430 304 57579.0697 52951.6410 305 40309.5053 57579.0697 306 52829.7391 40309.5053 307 24349.6867 52829.7391 308 30803.0719 24349.6867 309 52707.8373 30803.0719 310 52667.2033 52707.8373 311 -218.1321 52667.2033 312 52585.9354 -218.1321 313 38064.6083 52585.9354 314 44310.8329 38064.6083 315 60938.2202 44310.8329 316 52486.4997 60938.2202 317 4503.8971 52486.4997 318 47102.6457 4503.8971 319 52301.4978 47102.6457 320 52260.8638 52301.4978 321 52220.2299 52260.8638 322 52179.5959 52220.2299 323 52138.9620 52179.5959 324 52098.3280 52138.9620 325 47651.5540 52098.3280 326 52017.0601 47651.5540 327 51976.4262 52017.0601 328 51935.7922 51976.4262 329 51895.1583 51935.7922 330 51854.5243 51895.1583 331 51813.8903 51854.5243 332 51773.2564 51813.8903 333 44356.7007 51773.2564 334 50624.5211 44356.7007 335 40584.2542 50624.5211 336 -7323.1072 40584.2542 337 51570.0866 -7323.1072 338 48132.7798 51570.0866 339 42233.0392 48132.7798 340 8842.5627 42233.0392 341 51407.5508 8842.5627 342 51366.9169 51407.5508 343 51326.2829 51366.9169 344 16987.8690 51326.2829 345 51245.0150 16987.8690 346 36348.8650 51245.0150 347 51163.7471 36348.8650 348 51123.1132 51163.7471 349 9654.7960 51123.1132 350 47803.9920 9654.7960 351 40406.2766 47803.9920 352 39793.7947 40406.2766 353 43575.4660 39793.7947 354 50879.3095 43575.4660 355 50838.6755 50879.3095 356 1743.1207 50838.6755 357 50757.4076 1743.1207 358 50716.7737 50757.4076 359 50676.1397 50716.7737 360 46459.4107 50676.1397 361 50594.8718 46459.4107 362 50554.2379 50594.8718 363 -25813.1803 50554.2379 364 39900.4701 -25813.1803 365 50432.3360 39900.4701 366 50391.7021 50432.3360 367 50351.0681 50391.7021 368 50310.4342 50351.0681 369 38364.2479 50310.4342 370 50229.1663 38364.2479 371 42559.9748 50229.1663 372 33435.3504 42559.9748 373 28591.0951 33435.3504 374 50066.6304 28591.0951 375 50025.9965 50066.6304 376 23185.1999 50025.9965 377 36749.7733 23185.1999 378 -4565.5933 36749.7733 379 49863.4607 -4565.5933 380 49822.8267 49863.4607 381 790.6458 49822.8267 382 18691.7412 790.6458 383 49700.9249 18691.7412 384 34197.1859 49700.9249 385 -123899.8243 34197.1859 386 56978.0222 -123899.8243 387 20080.9755 56978.0222 388 44140.9672 20080.9755 389 24572.7433 44140.9672 390 34019.8369 24572.7433 391 12256.4655 34019.8369 392 -25124.4773 12256.4655 393 68245.7197 -25124.4773 394 -103693.0314 68245.7197 395 405.4308 -103693.0314 396 -144772.2788 405.4308 397 -206035.0405 -144772.2788 398 5987.5583 -206035.0405 399 -60992.2552 5987.5583 400 -4167.1290 -60992.2552 401 36067.3258 -4167.1290 402 13277.9864 36067.3258 403 -56796.4158 13277.9864 404 -15646.9222 -56796.4158 405 6941.8966 -15646.9222 406 13587.6916 6941.8966 407 6583.6352 13587.6916 408 16283.9769 6583.6352 409 -68871.9896 16283.9769 410 -60102.7556 -68871.9896 411 -168604.4202 -60102.7556 412 -32789.7838 -168604.4202 413 51958.0975 -32789.7838 414 -185299.1439 51958.0975 415 -106517.0924 -185299.1439 416 85830.2016 -106517.0924 417 -24856.5535 85830.2016 418 -58360.9935 -24856.5535 419 -44536.4016 -58360.9935 420 -51585.4272 -44536.4016 421 -185152.7264 -51585.4272 422 -552875.0594 -185152.7264 423 -465582.0946 -552875.0594 424 -56199.7327 -465582.0946 425 -259773.3347 -56199.7327 426 -319019.5228 -259773.3347 427 -1152897.9866 -319019.5228 428 -280631.4909 -1152897.9866 429 -345807.6229 -280631.4909 430 -689241.6941 -345807.6229 > 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/753ro1291213940.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/8fd8r1291213940.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/9fd8r1291213940.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/1084pu1291213940.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/11tmo01291213940.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/12f5m61291213940.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/13tfkf1291213940.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/14ef1k1291213940.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/1576051291213940.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/163yye1291213940.tab") + } > > try(system("convert tmp/11lt11291213940.ps tmp/11lt11291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/21lt11291213940.ps tmp/21lt11291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/31lt11291213940.ps tmp/31lt11291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/4cual1291213940.ps tmp/4cual1291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/5cual1291213940.ps tmp/5cual1291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/653ro1291213940.ps tmp/653ro1291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/753ro1291213940.ps tmp/753ro1291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/8fd8r1291213940.ps tmp/8fd8r1291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/9fd8r1291213940.ps tmp/9fd8r1291213940.png",intern=TRUE)) character(0) > try(system("convert tmp/1084pu1291213940.ps tmp/1084pu1291213940.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.656 3.204 14.140