R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(162556 + ,807 + ,213118 + ,6282154 + ,29790 + ,444 + ,81767 + ,4321023 + ,87550 + ,412 + ,153198 + ,4111912 + ,84738 + ,428 + ,-26007 + ,223193 + ,54660 + ,315 + ,126942 + ,1491348 + ,42634 + ,168 + ,157214 + ,1629616 + ,40949 + ,263 + ,129352 + ,1398893 + ,45187 + ,267 + ,234817 + ,1926517 + ,37704 + ,228 + ,60448 + ,983660 + ,16275 + ,129 + ,47818 + ,1443586 + ,25830 + ,104 + ,245546 + ,1073089 + ,12679 + ,122 + ,48020 + ,984885 + ,18014 + ,393 + ,-1710 + ,1405225 + ,43556 + ,190 + ,32648 + ,227132 + ,24811 + ,280 + ,95350 + ,929118 + ,6575 + ,63 + ,151352 + ,1071292 + ,7123 + ,102 + ,288170 + ,638830 + ,21950 + ,265 + ,114337 + ,856956 + ,37597 + ,234 + ,37884 + ,992426 + ,17821 + ,277 + ,122844 + ,444477 + ,12988 + ,73 + ,82340 + ,857217 + ,22330 + ,67 + ,79801 + ,711969 + ,13326 + ,103 + ,165548 + ,702380 + ,16189 + ,290 + ,116384 + ,358589 + ,7146 + ,83 + ,134028 + ,297978 + ,15824 + ,56 + ,63838 + ,585715 + ,27664 + ,236 + ,74996 + ,657954 + ,11920 + ,73 + ,31080 + ,209458 + ,8568 + ,34 + ,32168 + ,786690 + ,14416 + ,139 + ,49857 + ,439798 + ,3369 + ,26 + ,87161 + ,688779 + ,11819 + ,70 + ,106113 + ,574339 + ,6984 + ,40 + ,80570 + ,741409 + ,4519 + ,42 + ,102129 + ,597793 + ,2220 + ,12 + ,301670 + ,644190 + ,18562 + ,211 + ,102313 + ,377934 + ,10327 + ,74 + ,88577 + ,640273 + ,5336 + ,80 + ,112477 + ,697458 + ,2365 + ,83 + ,191778 + ,550608 + ,4069 + ,131 + ,79804 + ,207393 + ,8636 + ,203 + ,128294 + ,301607 + ,13718 + ,56 + ,96448 + ,345783 + ,4525 + ,89 + ,93811 + ,501749 + ,6869 + ,88 + ,117520 + ,379983 + ,4628 + ,39 + ,69159 + ,387475 + ,3689 + ,25 + ,101792 + ,377305 + ,4891 + ,49 + ,210568 + ,370837 + ,7489 + ,149 + ,136996 + ,430866 + ,4901 + ,58 + ,121920 + ,469107 + ,2284 + ,41 + ,76403 + ,194493 + ,3160 + ,90 + ,108094 + ,530670 + ,4150 + ,136 + ,134759 + ,518365 + ,7285 + ,97 + ,188873 + ,491303 + ,1134 + ,63 + ,146216 + ,527021 + ,4658 + ,114 + ,156608 + ,233773 + ,2384 + ,77 + ,61348 + ,405972 + ,3748 + ,6 + ,50350 + ,652925 + ,5371 + ,47 + ,87720 + ,446211 + ,1285 + ,51 + ,99489 + ,341340 + ,9327 + ,85 + ,87419 + ,387699 + ,5565 + ,43 + ,94355 + ,493408 + ,1528 + ,32 + ,60326 + ,146494 + ,3122 + ,25 + ,94670 + ,414462 + ,7561 + ,77 + ,82425 + ,364304 + ,2675 + ,54 + ,59017 + ,355178 + ,13253 + ,251 + ,90829 + ,357760 + ,880 + ,15 + ,80791 + ,261216 + ,2053 + ,44 + ,100423 + ,397144 + ,1424 + ,73 + ,131116 + ,374943 + ,4036 + ,85 + ,100269 + ,424898 + ,3045 + ,49 + ,27330 + ,202055 + ,5119 + ,38 + ,39039 + ,378525 + ,1431 + ,35 + ,106885 + ,310768 + ,554 + ,9 + ,79285 + ,325738 + ,1975 + ,34 + ,118881 + ,394510 + ,1765 + ,20 + ,77623 + ,247060 + ,1012 + ,29 + ,114768 + ,368078 + ,810 + ,11 + ,74015 + ,236761 + ,1280 + ,52 + ,69465 + ,312378 + ,666 + ,13 + ,117869 + ,339836 + ,1380 + ,29 + ,60982 + ,347385 + ,4677 + ,66 + ,90131 + ,426280 + ,876 + ,33 + ,138971 + ,352850 + ,814 + ,15 + ,39625 + ,301881 + ,514 + ,15 + ,102725 + ,377516 + ,5692 + ,68 + ,64239 + ,357312 + ,3642 + ,100 + ,90262 + ,458343 + ,540 + ,13 + ,103960 + ,354228 + ,2099 + ,45 + ,106611 + ,308636 + ,567 + ,14 + ,103345 + ,386212 + ,2001 + ,36 + ,95551 + ,393343 + ,2949 + ,40 + ,82903 + ,378509 + ,2253 + ,68 + ,63593 + ,452469 + ,6533 + ,29 + ,126910 + ,364839 + ,1889 + ,43 + ,37527 + ,358649 + ,3055 + ,30 + ,60247 + ,376641 + ,272 + ,9 + ,112995 + ,429112 + ,1414 + ,22 + ,70184 + ,330546 + ,2564 + ,19 + ,130140 + ,403560 + ,1383 + ,9 + ,73221 + ,317892 + ,1261 + ,31 + ,76114 + ,307528 + ,975 + ,19 + ,90534 + ,235133 + ,3366 + ,55 + ,108479 + ,299243 + ,576 + ,8 + ,113761 + ,314073 + ,1686 + ,28 + ,68696 + ,368186 + ,746 + ,29 + ,71561 + ,269661 + ,3192 + ,48 + ,59831 + ,125390 + ,2045 + ,16 + ,97890 + ,510834 + ,5702 + ,47 + ,101481 + ,321896 + ,1932 + ,20 + ,72954 + ,249898 + ,936 + ,22 + ,67939 + ,408881 + ,3437 + ,33 + ,48022 + ,158492 + ,5131 + ,44 + ,86111 + ,292154 + ,2397 + ,13 + ,74020 + ,289513 + ,1389 + ,6 + ,57530 + ,378049 + ,1503 + ,35 + ,56364 + ,343466 + ,402 + ,8 + ,84990 + ,332743 + ,2239 + ,17 + ,88590 + ,442882 + ,2234 + ,11 + ,77200 + ,214215 + ,837 + ,21 + ,61262 + ,315688 + ,10579 + ,92 + ,110309 + ,375195 + ,875 + ,12 + ,67000 + ,334280 + ,1585 + ,112 + ,93099 + ,355864 + ,1659 + ,25 + ,107577 + ,480382 + ,2647 + ,17 + ,62920 + ,353058 + ,3294 + ,23 + ,75832 + ,217193 + ,0 + ,0 + ,60720 + ,315380 + ,94 + ,10 + ,60793 + ,314533 + ,422 + ,23 + ,57935 + ,318056 + ,0 + ,0 + ,60720 + ,315380 + ,34 + ,7 + ,60630 + ,314353 + ,1558 + ,25 + ,55637 + ,369448 + ,0 + ,1 + ,60720 + ,315380 + ,43 + ,20 + ,60887 + ,312846 + ,645 + ,4 + ,60720 + ,312075 + ,316 + ,4 + ,60505 + ,315009 + ,115 + ,10 + ,60945 + ,318903 + ,5 + ,1 + ,60720 + ,314887 + ,897 + ,4 + ,60720 + ,314913 + ,0 + ,0 + ,60720 + ,315380 + ,389 + ,8 + ,58990 + ,325506 + ,0 + ,0 + ,60720 + ,315380 + ,1002 + ,11 + ,56750 + ,298568 + ,36 + ,4 + ,60894 + ,315834 + ,460 + ,15 + ,63346 + ,329784 + ,309 + ,9 + ,56535 + ,312878 + ,0 + ,0 + ,60720 + ,315380 + ,9 + ,7 + ,60835 + ,314987 + ,271 + ,2 + ,60720 + ,325249 + ,14 + ,0 + ,61016 + ,315877 + ,520 + ,7 + ,58650 + ,291650 + ,1766 + ,46 + ,60438 + ,305959 + ,0 + ,5 + ,60720 + ,315380 + ,458 + ,7 + ,58625 + ,297765 + ,20 + ,2 + ,60938 + ,315245 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,98 + ,2 + ,61490 + ,315236 + ,405 + ,5 + ,60845 + ,336425 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,483 + ,7 + ,60830 + ,306268 + ,454 + ,24 + ,63261 + ,302187 + ,47 + ,1 + ,60720 + ,314882 + ,0 + ,0 + ,60720 + ,315380 + ,757 + ,18 + ,45689 + ,382712 + ,4655 + ,55 + ,60720 + ,341570 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,36 + ,3 + ,61564 + ,312412 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,0 + ,60720 + ,315380 + ,126 + ,8 + ,60951 + ,315547 + ,400 + ,113 + ,60720 + ,313267 + ,71 + ,0 + ,60745 + ,316176 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,972 + ,19 + ,71642 + ,359335 + ,531 + ,11 + ,71641 + ,330068 + ,2461 + ,25 + ,55792 + ,314289 + ,378 + ,16 + ,71873 + ,297413 + ,23 + ,5 + ,62555 + ,314806 + ,638 + ,11 + ,60370 + ,333210 + ,2300 + ,23 + ,64873 + ,352108 + ,149 + ,6 + ,62041 + ,313332 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,7 + ,59500 + ,318745 + ,0 + ,0 + ,60720 + ,315380 + ,141 + ,7 + ,61630 + ,315366 + ,0 + ,0 + ,60720 + ,315380 + ,28 + ,3 + ,60890 + ,315688 + ,0 + ,0 + ,60720 + ,315380 + ,4980 + ,89 + ,113521 + ,409642 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,472 + ,19 + ,80045 + ,269587 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,12 + ,50804 + ,300962 + ,496 + ,12 + ,87390 + ,325479 + ,10 + ,5 + ,61656 + ,316155 + ,63 + ,2 + ,65688 + ,318574 + ,0 + ,0 + ,60720 + ,315380 + ,1136 + ,26 + ,48522 + ,343613 + ,265 + ,3 + ,60720 + ,306948 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,11 + ,57640 + ,330059 + ,474 + ,10 + ,61977 + ,288985 + ,534 + ,5 + ,62620 + ,304485 + ,0 + ,2 + ,60720 + ,315380 + ,15 + ,6 + ,60831 + ,315688 + ,397 + ,7 + ,60646 + ,317736 + ,0 + ,2 + ,60720 + ,315380 + ,1866 + ,28 + ,56225 + ,322331 + ,288 + ,3 + ,60510 + ,296656 + ,0 + ,0 + ,60720 + ,315380 + ,3 + ,1 + ,60698 + ,315354 + ,468 + ,20 + ,60720 + ,312161 + ,20 + ,1 + ,60805 + ,315576 + ,278 + ,22 + ,61404 + ,314922 + ,61 + ,9 + ,60720 + ,314551 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,2 + ,65276 + ,312339 + ,0 + ,0 + ,60720 + ,315380 + ,317 + ,7 + ,63915 + ,298700 + ,738 + ,9 + ,60720 + ,321376 + ,0 + ,0 + ,60720 + ,315380 + ,368 + ,13 + ,61686 + ,303230 + ,0 + ,0 + ,60720 + ,315380 + ,2 + ,0 + ,60743 + ,315487 + ,0 + ,0 + ,60720 + ,315380 + ,53 + ,6 + ,60349 + ,315793 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,94 + ,3 + ,61360 + ,312887 + ,0 + ,0 + ,60720 + ,315380 + ,24 + ,7 + ,59818 + ,315637 + ,2332 + ,2 + ,72680 + ,324385 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,131 + ,15 + ,61808 + ,308989 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,0 + ,60720 + ,315380 + ,167 + ,1 + ,64245 + ,307322 + ,622 + ,38 + ,73007 + ,304376 + ,2328 + ,57 + ,82732 + ,253588 + ,0 + ,0 + ,60720 + ,315380 + ,365 + ,7 + ,54820 + ,309560 + ,364 + ,26 + ,47705 + ,298466 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,13 + ,72835 + ,343929 + ,307 + ,10 + ,58856 + ,331955 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,9 + ,77655 + ,381180 + ,0 + ,0 + ,60720 + ,315380 + ,138 + ,26 + ,69817 + ,331420 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,125 + ,19 + ,60798 + ,310201 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,12 + ,62452 + ,320016 + ,335 + ,23 + ,64175 + ,320398 + ,0 + ,0 + ,60720 + ,315380 + ,1324 + ,29 + ,67440 + ,291841 + ,176 + ,8 + ,68136 + ,310670 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,26 + ,56726 + ,313491 + ,0 + ,0 + ,60720 + ,315380 + ,333 + ,9 + ,70811 + ,331323 + ,0 + ,0 + ,60720 + ,315380 + ,601 + ,5 + ,60720 + ,319210 + ,30 + ,3 + ,62045 + ,318098 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,13 + ,54323 + ,292754 + ,0 + ,0 + ,60720 + ,315380 + ,165 + ,12 + ,62841 + ,325176 + ,453 + ,19 + ,81125 + ,365959 + ,0 + ,0 + ,60720 + ,315380 + ,53 + ,10 + ,59506 + ,302409 + ,382 + ,9 + ,59365 + ,340968 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,9 + ,60720 + ,315380 + ,30 + ,4 + ,60798 + ,313164 + ,290 + ,1 + ,58790 + ,301164 + ,0 + ,1 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,366 + ,14 + ,61808 + ,344425 + ,2 + ,12 + ,60735 + ,315394 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,19 + ,64016 + ,316647 + ,384 + ,17 + ,54683 + ,309836 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,365 + ,32 + ,87192 + ,346611 + ,0 + ,0 + ,60720 + ,315380 + ,49 + ,14 + ,64107 + ,322031 + ,3 + ,8 + ,60761 + ,315656 + ,133 + ,4 + ,65990 + ,339445 + ,32 + ,0 + ,59988 + ,314964 + ,368 + ,20 + ,61167 + ,297141 + ,1 + ,5 + ,60719 + ,315372 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,22 + ,1 + ,60722 + ,312502 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,96 + ,4 + ,60379 + ,313729 + ,1 + ,1 + ,60727 + ,315388 + ,314 + ,4 + ,60720 + ,315371 + ,844 + ,20 + ,60925 + ,296139 + ,0 + ,0 + ,60720 + ,315380 + ,26 + ,1 + ,60896 + ,313880 + ,125 + ,10 + ,59734 + ,317698 + ,304 + ,12 + ,62969 + ,295580 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,621 + ,13 + ,60720 + ,308256 + ,0 + ,0 + ,60720 + ,315380 + ,119 + ,3 + ,59118 + ,303677 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,1595 + ,10 + ,60720 + ,319369 + ,312 + ,3 + ,58598 + ,318690 + ,60 + ,7 + ,61124 + ,314049 + ,587 + ,10 + ,59595 + ,325699 + ,135 + ,1 + ,62065 + ,314210 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,514 + ,15 + ,78780 + ,322378 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,1763 + ,28 + ,61600 + ,308336 + ,180 + ,9 + ,59635 + ,316386 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,7 + ,60720 + ,315553 + ,0 + ,0 + ,60720 + ,315380 + ,448 + ,7 + ,59781 + ,323361 + ,227 + ,7 + ,76644 + ,336639 + ,174 + ,3 + ,64820 + ,307424 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,121 + ,11 + ,56178 + ,295370 + ,607 + ,7 + ,60436 + ,322340 + ,2212 + ,10 + ,60720 + ,319864 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,60720 + ,315380 + ,530 + ,18 + ,73433 + ,317291 + ,571 + ,14 + ,41477 + ,280398 + ,0 + ,0 + ,60720 + ,315380 + ,78 + ,12 + ,62700 + ,317330 + ,2489 + ,29 + ,67804 + ,238125 + ,131 + ,3 + ,59661 + ,327071 + ,923 + ,6 + ,58620 + ,309038 + ,72 + ,3 + ,60398 + ,314210 + ,572 + ,8 + ,58580 + ,307930 + ,397 + ,10 + ,62710 + ,322327 + ,450 + ,6 + ,59325 + ,292136 + ,622 + ,8 + ,60950 + ,263276 + ,694 + ,6 + ,68060 + ,367655 + ,3425 + ,9 + ,83620 + ,283910 + ,562 + ,8 + ,58456 + ,283587 + ,4917 + ,26 + ,52811 + ,243650 + ,1442 + ,239 + ,121173 + ,438493 + ,529 + ,7 + ,63870 + ,296261 + ,2126 + ,41 + ,21001 + ,230621 + ,1061 + ,3 + ,70415 + ,304252 + ,776 + ,8 + ,64230 + ,333505 + ,611 + ,6 + ,59190 + ,296919 + ,1526 + ,21 + ,69351 + ,278990 + ,592 + ,7 + ,64270 + ,276898 + ,1182 + ,11 + ,70694 + ,327007 + ,621 + ,11 + ,68005 + ,317046 + ,989 + ,12 + ,58930 + ,304555 + ,438 + ,9 + ,58320 + ,298096 + ,726 + ,3 + ,69980 + ,231861 + ,1303 + ,57 + ,69863 + ,309422 + ,7419 + ,21 + ,63255 + ,286963 + ,1164 + ,15 + ,57320 + ,269753 + ,3310 + ,32 + ,75230 + ,448243 + ,1920 + ,11 + ,79420 + ,165404 + ,965 + ,2 + ,73490 + ,204325 + ,3256 + ,23 + ,35250 + ,407159 + ,1135 + ,20 + ,62285 + ,290476 + ,1270 + ,24 + ,69206 + ,275311 + ,661 + ,1 + ,65920 + ,246541 + ,1013 + ,1 + ,69770 + ,253468 + ,2844 + ,74 + ,72683 + ,240897 + ,11528 + ,68 + ,-14545 + ,-83265 + ,6526 + ,20 + ,55830 + ,-42143 + ,2264 + ,20 + ,55174 + ,272713 + ,5109 + ,82 + ,67038 + ,215362 + ,3999 + ,21 + ,51252 + ,42754 + ,35624 + ,244 + ,157278 + ,306275 + ,9252 + ,32 + ,79510 + ,253537 + ,15236 + ,86 + ,77440 + ,372631 + ,18073 + ,69 + ,27284 + ,-7170) + ,dim=c(4 + ,431) + ,dimnames=list(c('Costs' + ,'Orders' + ,'Dividends' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(4,431),dimnames=list(c('Costs','Orders','Dividends','Wealth'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '1' > #'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 Costs Orders Dividends Wealth 1 162556 807 213118 6282154 2 29790 444 81767 4321023 3 87550 412 153198 4111912 4 84738 428 -26007 223193 5 54660 315 126942 1491348 6 42634 168 157214 1629616 7 40949 263 129352 1398893 8 45187 267 234817 1926517 9 37704 228 60448 983660 10 16275 129 47818 1443586 11 25830 104 245546 1073089 12 12679 122 48020 984885 13 18014 393 -1710 1405225 14 43556 190 32648 227132 15 24811 280 95350 929118 16 6575 63 151352 1071292 17 7123 102 288170 638830 18 21950 265 114337 856956 19 37597 234 37884 992426 20 17821 277 122844 444477 21 12988 73 82340 857217 22 22330 67 79801 711969 23 13326 103 165548 702380 24 16189 290 116384 358589 25 7146 83 134028 297978 26 15824 56 63838 585715 27 27664 236 74996 657954 28 11920 73 31080 209458 29 8568 34 32168 786690 30 14416 139 49857 439798 31 3369 26 87161 688779 32 11819 70 106113 574339 33 6984 40 80570 741409 34 4519 42 102129 597793 35 2220 12 301670 644190 36 18562 211 102313 377934 37 10327 74 88577 640273 38 5336 80 112477 697458 39 2365 83 191778 550608 40 4069 131 79804 207393 41 8636 203 128294 301607 42 13718 56 96448 345783 43 4525 89 93811 501749 44 6869 88 117520 379983 45 4628 39 69159 387475 46 3689 25 101792 377305 47 4891 49 210568 370837 48 7489 149 136996 430866 49 4901 58 121920 469107 50 2284 41 76403 194493 51 3160 90 108094 530670 52 4150 136 134759 518365 53 7285 97 188873 491303 54 1134 63 146216 527021 55 4658 114 156608 233773 56 2384 77 61348 405972 57 3748 6 50350 652925 58 5371 47 87720 446211 59 1285 51 99489 341340 60 9327 85 87419 387699 61 5565 43 94355 493408 62 1528 32 60326 146494 63 3122 25 94670 414462 64 7561 77 82425 364304 65 2675 54 59017 355178 66 13253 251 90829 357760 67 880 15 80791 261216 68 2053 44 100423 397144 69 1424 73 131116 374943 70 4036 85 100269 424898 71 3045 49 27330 202055 72 5119 38 39039 378525 73 1431 35 106885 310768 74 554 9 79285 325738 75 1975 34 118881 394510 76 1765 20 77623 247060 77 1012 29 114768 368078 78 810 11 74015 236761 79 1280 52 69465 312378 80 666 13 117869 339836 81 1380 29 60982 347385 82 4677 66 90131 426280 83 876 33 138971 352850 84 814 15 39625 301881 85 514 15 102725 377516 86 5692 68 64239 357312 87 3642 100 90262 458343 88 540 13 103960 354228 89 2099 45 106611 308636 90 567 14 103345 386212 91 2001 36 95551 393343 92 2949 40 82903 378509 93 2253 68 63593 452469 94 6533 29 126910 364839 95 1889 43 37527 358649 96 3055 30 60247 376641 97 272 9 112995 429112 98 1414 22 70184 330546 99 2564 19 130140 403560 100 1383 9 73221 317892 101 1261 31 76114 307528 102 975 19 90534 235133 103 3366 55 108479 299243 104 576 8 113761 314073 105 1686 28 68696 368186 106 746 29 71561 269661 107 3192 48 59831 125390 108 2045 16 97890 510834 109 5702 47 101481 321896 110 1932 20 72954 249898 111 936 22 67939 408881 112 3437 33 48022 158492 113 5131 44 86111 292154 114 2397 13 74020 289513 115 1389 6 57530 378049 116 1503 35 56364 343466 117 402 8 84990 332743 118 2239 17 88590 442882 119 2234 11 77200 214215 120 837 21 61262 315688 121 10579 92 110309 375195 122 875 12 67000 334280 123 1585 112 93099 355864 124 1659 25 107577 480382 125 2647 17 62920 353058 126 3294 23 75832 217193 127 0 0 60720 315380 128 94 10 60793 314533 129 422 23 57935 318056 130 0 0 60720 315380 131 34 7 60630 314353 132 1558 25 55637 369448 133 0 1 60720 315380 134 43 20 60887 312846 135 645 4 60720 312075 136 316 4 60505 315009 137 115 10 60945 318903 138 5 1 60720 314887 139 897 4 60720 314913 140 0 0 60720 315380 141 389 8 58990 325506 142 0 0 60720 315380 143 1002 11 56750 298568 144 36 4 60894 315834 145 460 15 63346 329784 146 309 9 56535 312878 147 0 0 60720 315380 148 9 7 60835 314987 149 271 2 60720 325249 150 14 0 61016 315877 151 520 7 58650 291650 152 1766 46 60438 305959 153 0 5 60720 315380 154 458 7 58625 297765 155 20 2 60938 315245 156 0 0 60720 315380 157 0 0 60720 315380 158 98 2 61490 315236 159 405 5 60845 336425 160 0 0 60720 315380 161 0 0 60720 315380 162 0 0 60720 315380 163 0 0 60720 315380 164 483 7 60830 306268 165 454 24 63261 302187 166 47 1 60720 314882 167 0 0 60720 315380 168 757 18 45689 382712 169 4655 55 60720 341570 170 0 0 60720 315380 171 0 0 60720 315380 172 36 3 61564 312412 173 0 0 60720 315380 174 203 9 61938 309596 175 0 0 60720 315380 176 126 8 60951 315547 177 400 113 60720 313267 178 71 0 60745 316176 179 0 0 60720 315380 180 0 0 60720 315380 181 972 19 71642 359335 182 531 11 71641 330068 183 2461 25 55792 314289 184 378 16 71873 297413 185 23 5 62555 314806 186 638 11 60370 333210 187 2300 23 64873 352108 188 149 6 62041 313332 189 226 5 65745 291787 190 0 0 60720 315380 191 275 7 59500 318745 192 0 0 60720 315380 193 141 7 61630 315366 194 0 0 60720 315380 195 28 3 60890 315688 196 0 0 60720 315380 197 4980 89 113521 409642 198 0 0 60720 315380 199 0 0 60720 315380 200 472 19 80045 269587 201 0 0 60720 315380 202 0 0 60720 315380 203 0 0 60720 315380 204 203 12 50804 300962 205 496 12 87390 325479 206 10 5 61656 316155 207 63 2 65688 318574 208 0 0 60720 315380 209 1136 26 48522 343613 210 265 3 60720 306948 211 0 0 60720 315380 212 0 0 60720 315380 213 267 11 57640 330059 214 474 10 61977 288985 215 534 5 62620 304485 216 0 2 60720 315380 217 15 6 60831 315688 218 397 7 60646 317736 219 0 2 60720 315380 220 1866 28 56225 322331 221 288 3 60510 296656 222 0 0 60720 315380 223 3 1 60698 315354 224 468 20 60720 312161 225 20 1 60805 315576 226 278 22 61404 314922 227 61 9 60720 314551 228 0 0 60720 315380 229 192 2 65276 312339 230 0 0 60720 315380 231 317 7 63915 298700 232 738 9 60720 321376 233 0 0 60720 315380 234 368 13 61686 303230 235 0 0 60720 315380 236 2 0 60743 315487 237 0 0 60720 315380 238 53 6 60349 315793 239 0 0 60720 315380 240 0 0 60720 315380 241 0 0 60720 315380 242 94 3 61360 312887 243 0 0 60720 315380 244 24 7 59818 315637 245 2332 2 72680 324385 246 0 0 60720 315380 247 0 0 60720 315380 248 131 15 61808 308989 249 0 0 60720 315380 250 0 0 60720 315380 251 206 9 53110 296702 252 0 0 60720 315380 253 167 1 64245 307322 254 622 38 73007 304376 255 2328 57 82732 253588 256 0 0 60720 315380 257 365 7 54820 309560 258 364 26 47705 298466 259 0 0 60720 315380 260 0 0 60720 315380 261 0 0 60720 315380 262 0 0 60720 315380 263 226 13 72835 343929 264 307 10 58856 331955 265 0 0 60720 315380 266 0 0 60720 315380 267 0 0 60720 315380 268 188 9 77655 381180 269 0 0 60720 315380 270 138 26 69817 331420 271 0 0 60720 315380 272 0 0 60720 315380 273 0 0 60720 315380 274 125 19 60798 310201 275 0 0 60720 315380 276 282 12 62452 320016 277 335 23 64175 320398 278 0 0 60720 315380 279 1324 29 67440 291841 280 176 8 68136 310670 281 0 0 60720 315380 282 0 0 60720 315380 283 249 26 56726 313491 284 0 0 60720 315380 285 333 9 70811 331323 286 0 0 60720 315380 287 601 5 60720 319210 288 30 3 62045 318098 289 0 0 60720 315380 290 249 13 54323 292754 291 0 0 60720 315380 292 165 12 62841 325176 293 453 19 81125 365959 294 0 0 60720 315380 295 53 10 59506 302409 296 382 9 59365 340968 297 0 0 60720 315380 298 0 0 60720 315380 299 0 0 60720 315380 300 0 9 60720 315380 301 30 4 60798 313164 302 290 1 58790 301164 303 0 1 60720 315380 304 0 0 60720 315380 305 366 14 61808 344425 306 2 12 60735 315394 307 0 0 60720 315380 308 209 19 64016 316647 309 384 17 54683 309836 310 0 0 60720 315380 311 0 0 60720 315380 312 365 32 87192 346611 313 0 0 60720 315380 314 49 14 64107 322031 315 3 8 60761 315656 316 133 4 65990 339445 317 32 0 59988 314964 318 368 20 61167 297141 319 1 5 60719 315372 320 0 0 60720 315380 321 0 0 60720 315380 322 0 0 60720 315380 323 0 0 60720 315380 324 0 0 60720 315380 325 0 0 60720 315380 326 22 1 60722 312502 327 0 0 60720 315380 328 0 0 60720 315380 329 0 0 60720 315380 330 0 0 60720 315380 331 0 0 60720 315380 332 0 0 60720 315380 333 0 0 60720 315380 334 96 4 60379 313729 335 1 1 60727 315388 336 314 4 60720 315371 337 844 20 60925 296139 338 0 0 60720 315380 339 26 1 60896 313880 340 125 10 59734 317698 341 304 12 62969 295580 342 0 0 60720 315380 343 0 0 60720 315380 344 0 0 60720 315380 345 621 13 60720 308256 346 0 0 60720 315380 347 119 3 59118 303677 348 0 0 60720 315380 349 0 0 60720 315380 350 1595 10 60720 319369 351 312 3 58598 318690 352 60 7 61124 314049 353 587 10 59595 325699 354 135 1 62065 314210 355 0 0 60720 315380 356 0 0 60720 315380 357 514 15 78780 322378 358 0 0 60720 315380 359 0 0 60720 315380 360 0 0 60720 315380 361 1 4 60722 315398 362 0 0 60720 315380 363 0 0 60720 315380 364 1763 28 61600 308336 365 180 9 59635 316386 366 0 0 60720 315380 367 0 0 60720 315380 368 0 0 60720 315380 369 0 0 60720 315380 370 218 7 60720 315553 371 0 0 60720 315380 372 448 7 59781 323361 373 227 7 76644 336639 374 174 3 64820 307424 375 0 0 60720 315380 376 0 0 60720 315380 377 121 11 56178 295370 378 607 7 60436 322340 379 2212 10 60720 319864 380 0 0 60720 315380 381 0 0 60720 315380 382 530 18 73433 317291 383 571 14 41477 280398 384 0 0 60720 315380 385 78 12 62700 317330 386 2489 29 67804 238125 387 131 3 59661 327071 388 923 6 58620 309038 389 72 3 60398 314210 390 572 8 58580 307930 391 397 10 62710 322327 392 450 6 59325 292136 393 622 8 60950 263276 394 694 6 68060 367655 395 3425 9 83620 283910 396 562 8 58456 283587 397 4917 26 52811 243650 398 1442 239 121173 438493 399 529 7 63870 296261 400 2126 41 21001 230621 401 1061 3 70415 304252 402 776 8 64230 333505 403 611 6 59190 296919 404 1526 21 69351 278990 405 592 7 64270 276898 406 1182 11 70694 327007 407 621 11 68005 317046 408 989 12 58930 304555 409 438 9 58320 298096 410 726 3 69980 231861 411 1303 57 69863 309422 412 7419 21 63255 286963 413 1164 15 57320 269753 414 3310 32 75230 448243 415 1920 11 79420 165404 416 965 2 73490 204325 417 3256 23 35250 407159 418 1135 20 62285 290476 419 1270 24 69206 275311 420 661 1 65920 246541 421 1013 1 69770 253468 422 2844 74 72683 240897 423 11528 68 -14545 -83265 424 6526 20 55830 -42143 425 2264 20 55174 272713 426 5109 82 67038 215362 427 3999 21 51252 42754 428 35624 244 157278 306275 429 9252 32 79510 253537 430 15236 86 77440 372631 431 18073 69 27284 -7170 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Orders Dividends Wealth -2.185e+03 1.034e+02 -1.737e-02 8.923e-03 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -51069.2 -521.0 339.6 425.6 40226.2 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -2.185e+03 6.285e+02 -3.476 0.00056 *** Orders 1.034e+02 5.049e+00 20.477 < 2e-16 *** Dividends -1.737e-02 8.381e-03 -2.073 0.03878 * Wealth 8.923e-03 8.734e-04 10.217 < 2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 5104 on 427 degrees of freedom Multiple R-squared: 0.8148, Adjusted R-squared: 0.8135 F-statistic: 626.3 on 3 and 427 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 4.085643e-25 2.042821e-25 [2,] 1.0000000 8.013174e-34 4.006587e-34 [3,] 1.0000000 1.743977e-35 8.719885e-36 [4,] 1.0000000 1.437569e-37 7.187845e-38 [5,] 1.0000000 4.368345e-39 2.184173e-39 [6,] 1.0000000 1.418758e-38 7.093788e-39 [7,] 1.0000000 4.755889e-64 2.377945e-64 [8,] 1.0000000 2.825377e-81 1.412688e-81 [9,] 1.0000000 3.521898e-87 1.760949e-87 [10,] 1.0000000 1.807901e-86 9.039504e-87 [11,] 1.0000000 1.508309e-92 7.541544e-93 [12,] 1.0000000 5.621955e-96 2.810978e-96 [13,] 1.0000000 9.309984e-101 4.654992e-101 [14,] 1.0000000 2.977978e-106 1.488989e-106 [15,] 1.0000000 8.854312e-107 4.427156e-107 [16,] 1.0000000 3.168395e-116 1.584198e-116 [17,] 1.0000000 2.236986e-116 1.118493e-116 [18,] 1.0000000 2.888298e-121 1.444149e-121 [19,] 1.0000000 1.990625e-120 9.953126e-121 [20,] 1.0000000 5.992246e-125 2.996123e-125 [21,] 1.0000000 3.305700e-128 1.652850e-128 [22,] 1.0000000 2.263402e-129 1.131701e-129 [23,] 1.0000000 1.523253e-130 7.616263e-131 [24,] 1.0000000 7.901529e-131 3.950764e-131 [25,] 1.0000000 3.505588e-130 1.752794e-130 [26,] 1.0000000 6.312617e-132 3.156309e-132 [27,] 1.0000000 9.084338e-133 4.542169e-133 [28,] 1.0000000 2.714646e-132 1.357323e-132 [29,] 1.0000000 1.374016e-131 6.870080e-132 [30,] 1.0000000 2.601006e-132 1.300503e-132 [31,] 1.0000000 3.693725e-134 1.846863e-134 [32,] 1.0000000 5.972390e-134 2.986195e-134 [33,] 1.0000000 6.518888e-134 3.259444e-134 [34,] 1.0000000 4.089708e-135 2.044854e-135 [35,] 1.0000000 2.163515e-137 1.081757e-137 [36,] 1.0000000 3.641977e-142 1.820988e-142 [37,] 1.0000000 1.725256e-141 8.626279e-142 [38,] 1.0000000 1.066638e-140 5.333191e-141 [39,] 1.0000000 3.175310e-140 1.587655e-140 [40,] 1.0000000 1.175778e-139 5.878890e-140 [41,] 1.0000000 8.340194e-139 4.170097e-139 [42,] 1.0000000 8.269485e-139 4.134743e-139 [43,] 1.0000000 3.690466e-138 1.845233e-138 [44,] 1.0000000 1.882976e-137 9.414880e-138 [45,] 1.0000000 5.297066e-137 2.648533e-137 [46,] 1.0000000 8.547439e-138 4.273719e-138 [47,] 1.0000000 4.204441e-137 2.102221e-137 [48,] 1.0000000 1.299299e-136 6.496496e-137 [49,] 1.0000000 2.964563e-137 1.482281e-137 [50,] 1.0000000 9.859533e-137 4.929766e-137 [51,] 1.0000000 5.217121e-139 2.608561e-139 [52,] 1.0000000 7.227137e-139 3.613569e-139 [53,] 1.0000000 2.186093e-138 1.093046e-138 [54,] 1.0000000 1.427059e-138 7.135294e-139 [55,] 1.0000000 5.381793e-139 2.690897e-139 [56,] 1.0000000 1.110500e-138 5.552500e-139 [57,] 1.0000000 3.366352e-138 1.683176e-138 [58,] 1.0000000 7.258759e-138 3.629380e-138 [59,] 1.0000000 4.859801e-137 2.429900e-137 [60,] 1.0000000 2.585264e-138 1.292632e-138 [61,] 1.0000000 1.324930e-137 6.624649e-138 [62,] 1.0000000 8.739328e-137 4.369664e-137 [63,] 1.0000000 4.580119e-137 2.290059e-137 [64,] 1.0000000 2.162010e-136 1.081005e-136 [65,] 1.0000000 1.170902e-135 5.854512e-136 [66,] 1.0000000 6.704424e-136 3.352212e-136 [67,] 1.0000000 2.492221e-135 1.246110e-135 [68,] 1.0000000 1.510686e-134 7.553429e-135 [69,] 1.0000000 9.598418e-134 4.799209e-134 [70,] 1.0000000 4.851001e-133 2.425500e-133 [71,] 1.0000000 2.351729e-132 1.175865e-132 [72,] 1.0000000 1.000073e-131 5.000365e-132 [73,] 1.0000000 2.606310e-131 1.303155e-131 [74,] 1.0000000 1.235611e-130 6.178053e-131 [75,] 1.0000000 7.734243e-130 3.867122e-130 [76,] 1.0000000 3.714373e-129 1.857186e-129 [77,] 1.0000000 6.587755e-129 3.293878e-129 [78,] 1.0000000 3.826878e-128 1.913439e-128 [79,] 1.0000000 2.138458e-127 1.069229e-127 [80,] 1.0000000 7.873420e-127 3.936710e-127 [81,] 1.0000000 1.834445e-126 9.172225e-127 [82,] 1.0000000 9.145845e-126 4.572922e-126 [83,] 1.0000000 2.269572e-125 1.134786e-125 [84,] 1.0000000 1.249982e-124 6.249911e-125 [85,] 1.0000000 7.289162e-124 3.644581e-124 [86,] 1.0000000 4.004996e-123 2.002498e-123 [87,] 1.0000000 1.767328e-122 8.836642e-123 [88,] 1.0000000 1.498555e-122 7.492776e-123 [89,] 1.0000000 8.049316e-122 4.024658e-122 [90,] 1.0000000 2.241169e-121 1.120584e-121 [91,] 1.0000000 1.258239e-120 6.291194e-121 [92,] 1.0000000 6.965492e-120 3.482746e-120 [93,] 1.0000000 3.449690e-119 1.724845e-119 [94,] 1.0000000 1.662040e-118 8.310201e-119 [95,] 1.0000000 7.269682e-118 3.634841e-118 [96,] 1.0000000 1.528420e-117 7.642098e-118 [97,] 1.0000000 4.307527e-117 2.153763e-117 [98,] 1.0000000 1.240698e-116 6.203492e-117 [99,] 1.0000000 6.426506e-116 3.213253e-116 [100,] 1.0000000 1.665630e-115 8.328150e-116 [101,] 1.0000000 1.914778e-115 9.573891e-116 [102,] 1.0000000 2.280451e-115 1.140225e-115 [103,] 1.0000000 9.157633e-115 4.578816e-115 [104,] 1.0000000 3.861827e-114 1.930913e-114 [105,] 1.0000000 1.695488e-113 8.477442e-114 [106,] 1.0000000 6.454186e-113 3.227093e-113 [107,] 1.0000000 2.810095e-112 1.405047e-112 [108,] 1.0000000 1.213132e-111 6.065659e-112 [109,] 1.0000000 2.473181e-111 1.236591e-111 [110,] 1.0000000 1.271882e-110 6.359409e-111 [111,] 1.0000000 5.999493e-110 2.999747e-110 [112,] 1.0000000 1.090656e-109 5.453282e-110 [113,] 1.0000000 3.781119e-109 1.890560e-109 [114,] 1.0000000 1.885474e-108 9.427370e-109 [115,] 1.0000000 9.571150e-109 4.785575e-109 [116,] 1.0000000 4.673687e-108 2.336844e-108 [117,] 1.0000000 4.412795e-110 2.206398e-110 [118,] 1.0000000 1.450072e-109 7.250360e-110 [119,] 1.0000000 3.083309e-109 1.541654e-109 [120,] 1.0000000 1.234498e-108 6.172490e-109 [121,] 1.0000000 6.019137e-108 3.009569e-108 [122,] 1.0000000 2.944444e-107 1.472222e-107 [123,] 1.0000000 1.348614e-106 6.743069e-107 [124,] 1.0000000 6.482020e-106 3.241010e-106 [125,] 1.0000000 3.150290e-105 1.575145e-105 [126,] 1.0000000 1.182313e-104 5.911564e-105 [127,] 1.0000000 5.645668e-104 2.822834e-104 [128,] 1.0000000 2.253763e-103 1.126882e-103 [129,] 1.0000000 1.020912e-102 5.104559e-103 [130,] 1.0000000 4.790598e-102 2.395299e-102 [131,] 1.0000000 2.268329e-101 1.134164e-101 [132,] 1.0000000 1.058653e-100 5.293263e-101 [133,] 1.0000000 4.421614e-100 2.210807e-100 [134,] 1.0000000 2.034251e-99 1.017125e-99 [135,] 1.0000000 9.403457e-99 4.701728e-99 [136,] 1.0000000 4.287094e-98 2.143547e-98 [137,] 1.0000000 1.946287e-97 9.731437e-98 [138,] 1.0000000 8.930099e-97 4.465049e-97 [139,] 1.0000000 4.114978e-96 2.057489e-96 [140,] 1.0000000 1.881974e-95 9.409869e-96 [141,] 1.0000000 8.381252e-95 4.190626e-95 [142,] 1.0000000 3.749645e-94 1.874822e-94 [143,] 1.0000000 1.605834e-93 8.029171e-94 [144,] 1.0000000 7.047948e-93 3.523974e-93 [145,] 1.0000000 3.079251e-92 1.539625e-92 [146,] 1.0000000 1.038682e-91 5.193409e-92 [147,] 1.0000000 4.562480e-91 2.281240e-91 [148,] 1.0000000 1.979023e-90 9.895115e-91 [149,] 1.0000000 8.575381e-90 4.287690e-90 [150,] 1.0000000 3.666720e-89 1.833360e-89 [151,] 1.0000000 1.560848e-88 7.804241e-89 [152,] 1.0000000 6.650047e-88 3.325024e-88 [153,] 1.0000000 2.676778e-87 1.338389e-87 [154,] 1.0000000 1.124140e-86 5.620701e-87 [155,] 1.0000000 4.699754e-86 2.349877e-86 [156,] 1.0000000 1.956020e-85 9.780101e-86 [157,] 1.0000000 8.104243e-85 4.052121e-85 [158,] 1.0000000 3.369956e-84 1.684978e-84 [159,] 1.0000000 1.102936e-83 5.514679e-84 [160,] 1.0000000 4.522433e-83 2.261216e-83 [161,] 1.0000000 1.841137e-82 9.205687e-83 [162,] 1.0000000 5.106366e-82 2.553183e-82 [163,] 1.0000000 1.435321e-81 7.176606e-82 [164,] 1.0000000 5.776310e-81 2.888155e-81 [165,] 1.0000000 2.314197e-80 1.157098e-80 [166,] 1.0000000 9.307069e-80 4.653534e-80 [167,] 1.0000000 3.695118e-79 1.847559e-79 [168,] 1.0000000 1.455127e-78 7.275637e-79 [169,] 1.0000000 5.724983e-78 2.862492e-78 [170,] 1.0000000 2.262161e-77 1.131081e-77 [171,] 1.0000000 3.520369e-80 1.760185e-80 [172,] 1.0000000 1.403849e-79 7.019247e-80 [173,] 1.0000000 5.611591e-79 2.805795e-79 [174,] 1.0000000 2.233102e-78 1.116551e-78 [175,] 1.0000000 8.866462e-78 4.433231e-78 [176,] 1.0000000 3.544959e-77 1.772480e-77 [177,] 1.0000000 1.251536e-76 6.257681e-77 [178,] 1.0000000 3.988130e-76 1.994065e-76 [179,] 1.0000000 1.564439e-75 7.822195e-76 [180,] 1.0000000 6.029570e-75 3.014785e-75 [181,] 1.0000000 1.809786e-74 9.048929e-75 [182,] 1.0000000 7.023035e-74 3.511518e-74 [183,] 1.0000000 2.582580e-73 1.291290e-73 [184,] 1.0000000 9.845236e-73 4.922618e-73 [185,] 1.0000000 3.783150e-72 1.891575e-72 [186,] 1.0000000 1.429175e-71 7.145877e-72 [187,] 1.0000000 5.422773e-71 2.711386e-71 [188,] 1.0000000 2.030199e-70 1.015099e-70 [189,] 1.0000000 7.641352e-70 3.820676e-70 [190,] 1.0000000 2.835141e-69 1.417570e-69 [191,] 1.0000000 8.096974e-69 4.048487e-69 [192,] 1.0000000 2.980593e-68 1.490297e-68 [193,] 1.0000000 1.092294e-67 5.461468e-68 [194,] 1.0000000 2.074603e-67 1.037301e-67 [195,] 1.0000000 7.561976e-67 3.780988e-67 [196,] 1.0000000 2.744040e-66 1.372020e-66 [197,] 1.0000000 9.912855e-66 4.956427e-66 [198,] 1.0000000 3.544907e-65 1.772453e-65 [199,] 1.0000000 1.151289e-64 5.756443e-65 [200,] 1.0000000 4.136243e-64 2.068121e-64 [201,] 1.0000000 1.477256e-63 7.386280e-64 [202,] 1.0000000 5.221301e-63 2.610650e-63 [203,] 1.0000000 1.780844e-62 8.904221e-63 [204,] 1.0000000 6.257096e-62 3.128548e-62 [205,] 1.0000000 2.181629e-61 1.090814e-61 [206,] 1.0000000 7.572442e-61 3.786221e-61 [207,] 1.0000000 2.645128e-60 1.322564e-60 [208,] 1.0000000 8.726672e-60 4.363336e-60 [209,] 1.0000000 2.992599e-59 1.496299e-59 [210,] 1.0000000 1.027373e-58 5.136865e-59 [211,] 1.0000000 3.506199e-58 1.753099e-58 [212,] 1.0000000 1.192575e-57 5.962873e-58 [213,] 1.0000000 4.037924e-57 2.018962e-57 [214,] 1.0000000 1.345856e-56 6.729281e-57 [215,] 1.0000000 4.482498e-56 2.241249e-56 [216,] 1.0000000 1.487238e-55 7.436188e-56 [217,] 1.0000000 4.930612e-55 2.465306e-55 [218,] 1.0000000 1.516240e-54 7.581202e-55 [219,] 1.0000000 4.979406e-54 2.489703e-54 [220,] 1.0000000 1.417900e-53 7.089499e-54 [221,] 1.0000000 4.569938e-53 2.284969e-53 [222,] 1.0000000 1.477357e-52 7.386784e-53 [223,] 1.0000000 4.765218e-52 2.382609e-52 [224,] 1.0000000 1.526565e-51 7.632825e-52 [225,] 1.0000000 4.781891e-51 2.390946e-51 [226,] 1.0000000 1.503842e-50 7.519209e-51 [227,] 1.0000000 4.753015e-50 2.376507e-50 [228,] 1.0000000 1.441938e-49 7.209690e-50 [229,] 1.0000000 4.517130e-49 2.258565e-49 [230,] 1.0000000 1.408390e-48 7.041948e-49 [231,] 1.0000000 4.371662e-48 2.185831e-48 [232,] 1.0000000 1.358515e-47 6.792577e-48 [233,] 1.0000000 4.178088e-47 2.089044e-47 [234,] 1.0000000 1.279040e-46 6.395198e-47 [235,] 1.0000000 3.897451e-46 1.948725e-46 [236,] 1.0000000 1.188549e-45 5.942746e-46 [237,] 1.0000000 3.588234e-45 1.794117e-45 [238,] 1.0000000 1.081169e-44 5.405846e-45 [239,] 1.0000000 2.027337e-44 1.013668e-44 [240,] 1.0000000 6.052591e-44 3.026296e-44 [241,] 1.0000000 1.798560e-43 8.992801e-44 [242,] 1.0000000 4.948368e-43 2.474184e-43 [243,] 1.0000000 1.457604e-42 7.288020e-43 [244,] 1.0000000 4.273413e-42 2.136706e-42 [245,] 1.0000000 1.242267e-41 6.211337e-42 [246,] 1.0000000 3.608257e-41 1.804128e-41 [247,] 1.0000000 1.044309e-40 5.221547e-41 [248,] 1.0000000 1.586334e-40 7.931668e-41 [249,] 1.0000000 1.184849e-40 5.924243e-41 [250,] 1.0000000 3.425570e-40 1.712785e-40 [251,] 1.0000000 9.897958e-40 4.948979e-40 [252,] 1.0000000 2.479625e-39 1.239812e-39 [253,] 1.0000000 7.074025e-39 3.537013e-39 [254,] 1.0000000 2.008473e-38 1.004237e-38 [255,] 1.0000000 5.675176e-38 2.837588e-38 [256,] 1.0000000 1.595884e-37 7.979421e-38 [257,] 1.0000000 4.431481e-37 2.215740e-37 [258,] 1.0000000 1.238607e-36 6.193037e-37 [259,] 1.0000000 3.432938e-36 1.716469e-36 [260,] 1.0000000 9.468629e-36 4.734314e-36 [261,] 1.0000000 2.598911e-35 1.299456e-35 [262,] 1.0000000 7.001339e-35 3.500669e-35 [263,] 1.0000000 1.902925e-34 9.514625e-35 [264,] 1.0000000 4.138979e-34 2.069489e-34 [265,] 1.0000000 1.116072e-33 5.580361e-34 [266,] 1.0000000 2.994655e-33 1.497328e-33 [267,] 1.0000000 7.995592e-33 3.997796e-33 [268,] 1.0000000 1.896396e-32 9.481981e-33 [269,] 1.0000000 5.018627e-32 2.509314e-32 [270,] 1.0000000 1.314031e-31 6.570153e-32 [271,] 1.0000000 3.065727e-31 1.532863e-31 [272,] 1.0000000 8.001928e-31 4.000964e-31 [273,] 1.0000000 1.782208e-30 8.911040e-31 [274,] 1.0000000 4.490005e-30 2.245003e-30 [275,] 1.0000000 1.156403e-29 5.782016e-30 [276,] 1.0000000 2.963207e-29 1.481603e-29 [277,] 1.0000000 6.474778e-29 3.237389e-29 [278,] 1.0000000 1.644660e-28 8.223302e-29 [279,] 1.0000000 4.161496e-28 2.080748e-28 [280,] 1.0000000 1.046268e-27 5.231341e-28 [281,] 1.0000000 2.579966e-27 1.289983e-27 [282,] 1.0000000 6.444166e-27 3.222083e-27 [283,] 1.0000000 1.595215e-26 7.976073e-27 [284,] 1.0000000 3.775190e-26 1.887595e-26 [285,] 1.0000000 9.252249e-26 4.626124e-26 [286,] 1.0000000 2.235179e-25 1.117589e-25 [287,] 1.0000000 5.330000e-25 2.665000e-25 [288,] 1.0000000 1.286164e-24 6.430820e-25 [289,] 1.0000000 2.980022e-24 1.490011e-24 [290,] 1.0000000 7.044267e-24 3.522133e-24 [291,] 1.0000000 1.673516e-23 8.367580e-24 [292,] 1.0000000 3.954373e-23 1.977186e-23 [293,] 1.0000000 9.293273e-23 4.646636e-23 [294,] 1.0000000 2.144975e-22 1.072487e-22 [295,] 1.0000000 4.995305e-22 2.497653e-22 [296,] 1.0000000 1.152722e-21 5.763610e-22 [297,] 1.0000000 2.655463e-21 1.327732e-21 [298,] 1.0000000 6.073007e-21 3.036504e-21 [299,] 1.0000000 1.383771e-20 6.918855e-21 [300,] 1.0000000 3.038388e-20 1.519194e-20 [301,] 1.0000000 6.835534e-20 3.417767e-20 [302,] 1.0000000 1.413329e-19 7.066646e-20 [303,] 1.0000000 3.080122e-19 1.540061e-19 [304,] 1.0000000 6.821276e-19 3.410638e-19 [305,] 1.0000000 1.501999e-18 7.509997e-19 [306,] 1.0000000 2.432095e-18 1.216048e-18 [307,] 1.0000000 5.313499e-18 2.656750e-18 [308,] 1.0000000 1.108769e-17 5.543843e-18 [309,] 1.0000000 2.373140e-17 1.186570e-17 [310,] 1.0000000 5.093454e-17 2.546727e-17 [311,] 1.0000000 1.086684e-16 5.433421e-17 [312,] 1.0000000 2.073913e-16 1.036956e-16 [313,] 1.0000000 4.376292e-16 2.188146e-16 [314,] 1.0000000 9.192408e-16 4.596204e-16 [315,] 1.0000000 1.919187e-15 9.595933e-16 [316,] 1.0000000 3.982484e-15 1.991242e-15 [317,] 1.0000000 8.213410e-15 4.106705e-15 [318,] 1.0000000 1.683485e-14 8.417423e-15 [319,] 1.0000000 3.429210e-14 1.714605e-14 [320,] 1.0000000 6.946940e-14 3.473470e-14 [321,] 1.0000000 1.397416e-13 6.987078e-14 [322,] 1.0000000 2.793214e-13 1.396607e-13 [323,] 1.0000000 5.547674e-13 2.773837e-13 [324,] 1.0000000 1.094781e-12 5.473906e-13 [325,] 1.0000000 2.146516e-12 1.073258e-12 [326,] 1.0000000 4.181308e-12 2.090654e-12 [327,] 1.0000000 8.091732e-12 4.045866e-12 [328,] 1.0000000 1.557247e-11 7.786235e-12 [329,] 1.0000000 2.976498e-11 1.488249e-11 [330,] 1.0000000 5.647814e-11 2.823907e-11 [331,] 1.0000000 1.026319e-10 5.131593e-11 [332,] 1.0000000 1.921970e-10 9.609848e-11 [333,] 1.0000000 3.576974e-10 1.788487e-10 [334,] 1.0000000 6.571180e-10 3.285590e-10 [335,] 1.0000000 1.161354e-09 5.806772e-10 [336,] 1.0000000 2.117193e-09 1.058596e-09 [337,] 1.0000000 3.832523e-09 1.916262e-09 [338,] 1.0000000 6.888317e-09 3.444158e-09 [339,] 1.0000000 1.223454e-08 6.117269e-09 [340,] 1.0000000 2.167700e-08 1.083850e-08 [341,] 1.0000000 3.809290e-08 1.904645e-08 [342,] 1.0000000 6.651171e-08 3.325586e-08 [343,] 0.9999999 1.152718e-07 5.763589e-08 [344,] 0.9999999 1.915526e-07 9.577629e-08 [345,] 0.9999998 3.261596e-07 1.630798e-07 [346,] 0.9999997 5.509601e-07 2.754800e-07 [347,] 0.9999995 9.269720e-07 4.634860e-07 [348,] 0.9999992 1.547445e-06 7.737227e-07 [349,] 0.9999987 2.563232e-06 1.281616e-06 [350,] 0.9999979 4.212304e-06 2.106152e-06 [351,] 0.9999967 6.652101e-06 3.326051e-06 [352,] 0.9999946 1.076666e-05 5.383329e-06 [353,] 0.9999914 1.728466e-05 8.642328e-06 [354,] 0.9999862 2.752094e-05 1.376047e-05 [355,] 0.9999783 4.345261e-05 2.172631e-05 [356,] 0.9999660 6.804043e-05 3.402022e-05 [357,] 0.9999472 1.056407e-04 5.282037e-05 [358,] 0.9999190 1.620307e-04 8.101535e-05 [359,] 0.9998766 2.468922e-04 1.234461e-04 [360,] 0.9998132 3.735697e-04 1.867849e-04 [361,] 0.9997199 5.602646e-04 2.801323e-04 [362,] 0.9995836 8.327819e-04 4.163910e-04 [363,] 0.9993866 1.226718e-03 6.133588e-04 [364,] 0.9991048 1.790309e-03 8.951545e-04 [365,] 0.9987054 2.589264e-03 1.294632e-03 [366,] 0.9981456 3.708726e-03 1.854363e-03 [367,] 0.9973738 5.252402e-03 2.626201e-03 [368,] 0.9963130 7.374070e-03 3.687035e-03 [369,] 0.9948650 1.027003e-02 5.135017e-03 [370,] 0.9929176 1.416473e-02 7.082363e-03 [371,] 0.9904297 1.914052e-02 9.570261e-03 [372,] 0.9870664 2.586718e-02 1.293359e-02 [373,] 0.9833306 3.333880e-02 1.666940e-02 [374,] 0.9778618 4.427642e-02 2.213821e-02 [375,] 0.9708986 5.820288e-02 2.910144e-02 [376,] 0.9627698 7.446040e-02 3.723020e-02 [377,] 0.9520631 9.587390e-02 4.793695e-02 [378,] 0.9388658 1.222685e-01 6.113424e-02 [379,] 0.9233200 1.533599e-01 7.667996e-02 [380,] 0.9066587 1.866826e-01 9.334128e-02 [381,] 0.8845274 2.309452e-01 1.154726e-01 [382,] 0.8589360 2.821281e-01 1.410640e-01 [383,] 0.8291497 3.417006e-01 1.708503e-01 [384,] 0.7953074 4.093852e-01 2.046926e-01 [385,] 0.7575010 4.849981e-01 2.424990e-01 [386,] 0.7160851 5.678298e-01 2.839149e-01 [387,] 0.6732332 6.535336e-01 3.267668e-01 [388,] 0.6271913 7.456174e-01 3.728087e-01 [389,] 0.5808355 8.383290e-01 4.191645e-01 [390,] 0.5289826 9.420348e-01 4.710174e-01 [391,] 0.4832403 9.664806e-01 5.167597e-01 [392,] 0.9993777 1.244643e-03 6.223214e-04 [393,] 0.9989266 2.146779e-03 1.073389e-03 [394,] 0.9986825 2.634923e-03 1.317461e-03 [395,] 0.9978588 4.282482e-03 2.141241e-03 [396,] 0.9964592 7.081687e-03 3.540843e-03 [397,] 0.9942167 1.156654e-02 5.783268e-03 [398,] 0.9910051 1.798971e-02 8.994853e-03 [399,] 0.9858042 2.839162e-02 1.419581e-02 [400,] 0.9781250 4.375003e-02 2.187502e-02 [401,] 0.9668825 6.623494e-02 3.311747e-02 [402,] 0.9508920 9.821603e-02 4.910802e-02 [403,] 0.9288966 1.422068e-01 7.110342e-02 [404,] 0.8993719 2.012562e-01 1.006281e-01 [405,] 0.9293866 1.412267e-01 7.061337e-02 [406,] 0.9344830 1.310339e-01 6.551697e-02 [407,] 0.9039307 1.921386e-01 9.606931e-02 [408,] 0.8627064 2.745871e-01 1.372936e-01 [409,] 0.8089671 3.820657e-01 1.910329e-01 [410,] 0.7405855 5.188291e-01 2.594145e-01 [411,] 0.6806729 6.386543e-01 3.193271e-01 [412,] 0.5903046 8.193908e-01 4.096954e-01 [413,] 0.5013504 9.972992e-01 4.986496e-01 [414,] 0.3948549 7.897097e-01 6.051451e-01 [415,] 0.2925943 5.851887e-01 7.074057e-01 [416,] 0.4531570 9.063140e-01 5.468430e-01 [417,] 0.3410712 6.821425e-01 6.589288e-01 [418,] 0.2134217 4.268434e-01 7.865783e-01 > postscript(file="/var/www/html/rcomp/tmp/1e5sx1292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2e5sx1292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3e5sx1292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/47fr01292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/57fr01292938441.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 2.894692e+04 -5.106920e+04 1.310643e+04 4.022619e+04 1.317328e+04 6 7 8 9 10 1.563856e+04 5.705697e+03 6.654402e+03 8.587580e+03 -6.928696e+03 11 12 13 14 15 1.195252e+04 -5.704364e+03 -3.300428e+04 2.463621e+04 -8.588906e+03 16 17 18 19 20 -4.683766e+03 -1.932186e+03 -8.925183e+03 7.389960e+03 -1.046652e+04 21 22 23 24 25 1.406486e+03 1.262081e+04 1.469915e+03 -1.278850e+04 4.187309e+02 26 27 28 29 30 8.101405e+03 8.794954e+02 5.227941e+03 7.765397e+02 -8.292800e+02 31 32 33 34 35 -1.766174e+03 3.484866e+03 -1.828464e+02 -1.198564e+03 2.657120e+03 36 37 38 39 40 -2.664355e+03 6.862752e+02 -5.020128e+03 -5.613173e+03 -7.755034e+03 41 42 43 44 45 -1.063072e+04 8.702920e+03 -5.339650e+03 -1.393799e+03 5.244980e+02 46 47 48 49 50 1.690748e+03 2.358851e+03 -7.196577e+03 -9.787577e+02 -1.784293e+02 51 52 53 54 55 -6.817959e+03 -1.001107e+04 -1.661993e+03 -5.357390e+03 -4.309310e+03 56 57 58 59 60 -5.949287e+03 3.611094e+02 2.387043e+02 -3.120623e+03 7.825609e+02 61 62 63 64 65 5.404169e+02 1.451091e+02 6.684533e+02 -3.428489e+01 -2.867449e+03 66 67 68 69 70 -1.212868e+04 5.867092e+02 -2.110575e+03 -5.006678e+03 -4.617115e+03 71 72 73 74 75 -1.164665e+03 6.754499e+02 -9.190005e+02 2.791772e+02 -8.104282e+02 76 77 78 79 80 1.026007e+03 -1.092054e+03 1.030780e+03 -3.492224e+03 5.221564e+02 81 82 83 84 85 -1.473887e+03 -2.200072e+03 -1.085251e+03 -5.573695e+02 -4.359706e+02 86 87 88 89 90 -1.226300e+03 -7.034342e+03 2.607924e+01 -1.270691e+03 -3.463990e+02 91 92 93 94 95 -1.386142e+03 -9.391039e+02 -5.525623e+03 4.668803e+03 -2.920440e+03 96 97 98 99 100 -1.761075e+02 -3.395656e+02 -4.059864e+02 1.444362e+03 1.072832e+03 101 102 103 104 105 -1.181124e+03 6.701451e+02 -9.213749e+02 1.107648e+03 -1.116078e+03 106 107 108 109 110 -1.230544e+03 3.344967e+02 -2.819851e+02 1.918069e+03 1.086564e+03 111 112 113 114 115 -1.621986e+03 1.629884e+03 1.655608e+03 1.940362e+03 5.796129e+02 116 117 118 119 120 -2.016522e+03 2.671847e+02 2.533866e+02 2.711297e+03 -9.020219e+02 121 122 123 124 125 1.820060e+03 3.296227e-01 -9.368362e+03 -1.158514e+03 1.016908e+03 126 127 128 129 130 2.480212e+03 4.256129e+02 -5.055148e+02 -1.602746e+03 4.256129e+02 131 132 133 134 135 -2.565545e+02 -1.172041e+03 3.222175e+02 -1.573782e+03 6.865223e+02 136 137 138 139 140 3.276064e+02 -5.208681e+02 3.316166e+02 9.131984e+02 4.256129e+02 141 142 143 144 145 -1.329631e+02 4.256129e+02 3.713046e+02 4.700327e+01 -7.482228e+02 146 147 148 149 150 -2.463302e+02 4.256129e+02 -2.836501e+02 4.017596e+02 4.403208e+02 151 152 153 154 155 3.976271e+02 -2.485410e+03 -9.136409e+01 2.806277e+02 2.438142e+02 156 157 158 159 160 4.256129e+02 4.256129e+02 3.314850e+02 1.280201e+02 4.256129e+02 161 162 163 164 165 4.256129e+02 4.256129e+02 4.256129e+02 2.680640e+02 -1.440006e+03 166 167 168 169 170 3.736612e+02 4.256129e+02 -1.540466e+03 -8.398307e+02 4.256129e+02 171 172 173 174 175 4.256129e+02 1.925743e+02 4.256129e+02 -2.291726e+02 4.256129e+02 176 177 178 179 180 -2.730270e+02 -1.083921e+04 4.899444e+02 4.256129e+02 4.256129e+02 181 182 183 184 185 -7.693573e+02 -1.220579e+02 2.258441e+02 -4.966188e+02 -3.136088e+01 186 187 188 189 190 -2.389170e+02 9.194398e+01 -4.533803e+00 4.324641e+02 4.256129e+02 191 192 193 194 195 -7.437756e+01 4.256129e+02 -1.412196e+02 4.256129e+02 1.436319e+02 196 197 198 199 200 4.256129e+02 -3.720324e+03 4.256129e+02 4.256129e+02 -3.225290e+02 201 202 203 204 205 4.256129e+02 4.256129e+02 4.256129e+02 -6.557587e+02 5.411853e+01 206 207 208 209 210 -7.201745e+01 3.396357e+02 4.256129e+02 -1.590523e+03 4.556666e+02 211 212 213 214 215 4.256129e+02 4.256129e+02 -6.292313e+02 1.230245e+02 5.728642e+02 216 217 218 219 220 2.188221e+02 -1.805793e+02 7.653652e+01 2.188221e+02 -7.435791e+02 221 222 223 224 225 5.668551e+02 4.256129e+02 3.250672e+02 -1.145571e+03 3.419453e+02 226 227 228 229 230 -1.555115e+03 -4.365484e+02 4.256129e+02 5.171134e+02 4.256129e+02 231 232 233 234 235 2.231932e+02 1.795512e+02 4.256129e+02 -4.253277e+02 4.256129e+02 236 237 238 239 240 4.270577e+02 4.256129e+02 -1.518905e+02 4.256129e+02 4.256129e+02 241 242 243 244 245 4.256129e+02 2.427915e+02 4.256129e+02 -2.921195e+02 2.678262e+03 246 247 248 249 250 4.256129e+02 4.256129e+02 -9.183872e+02 4.256129e+02 4.256129e+02 251 252 253 254 255 -2.644955e+02 4.256129e+02 6.223636e+02 -2.569747e+03 -2.206108e+03 256 257 258 259 260 4.256129e+02 1.627109e+01 -1.973864e+03 4.256129e+02 4.256129e+02 261 262 263 264 265 4.256129e+02 4.256129e+02 -7.367878e+02 -4.816273e+02 4.256129e+02 266 267 268 269 270 4.256129e+02 4.256129e+02 -6.098594e+02 4.256129e+02 -2.109743e+03 271 272 273 274 275 4.256129e+02 4.256129e+02 4.256129e+02 -1.366331e+03 4.256129e+02 276 277 278 279 280 -5.444077e+02 -1.602230e+03 4.256129e+02 -9.220581e+02 -5.467651e+01 281 282 283 284 285 4.256129e+02 4.256129e+02 -2.066203e+03 4.256129e+02 -1.388861e+02 286 287 288 289 290 4.256129e+02 4.754603e+02 1.441942e+02 4.256129e+02 -5.787738e+02 291 292 293 294 295 4.256129e+02 -7.006926e+02 -1.182706e+03 4.256129e+02 -4.606910e+02 296 297 298 299 300 -3.748129e+02 4.256129e+02 4.256129e+02 4.256129e+02 -5.049457e+02 301 302 303 304 305 6.316016e+01 7.055370e+02 3.222175e+02 4.256129e+02 -8.961924e+02 306 307 308 309 310 -8.129961e+02 4.256129e+02 -1.283940e+03 -1.003526e+03 4.256129e+02 311 312 313 314 315 4.256129e+02 -2.336793e+03 4.256129e+02 -9.734246e+02 -4.003007e+02 316 317 318 319 320 2.185694e+01 4.486071e+02 -1.103780e+03 -9.031008e+01 4.256129e+02 321 322 323 324 325 4.256129e+02 4.256129e+02 4.256129e+02 4.256129e+02 4.256129e+02 326 327 328 329 330 3.699330e+02 4.256129e+02 4.256129e+02 4.256129e+02 4.256129e+02 331 332 333 334 335 4.256129e+02 4.256129e+02 4.256129e+02 1.168389e+02 3.232677e+02 336 337 338 339 340 3.261116e+02 -6.230431e+02 4.256129e+02 3.646600e+02 -5.211557e+02 341 342 343 344 345 -2.953794e+02 4.256129e+02 4.256129e+02 4.256129e+02 -2.339587e+02 346 347 348 349 350 4.256129e+02 3.110211e+02 4.256129e+02 4.256129e+02 9.510645e+02 351 352 353 354 355 3.610234e+02 -2.192591e+02 -1.329647e+02 4.910256e+02 4.256129e+02 356 357 358 359 360 4.256129e+02 -3.599873e+02 4.256129e+02 4.256129e+02 4.256129e+02 361 362 363 364 365 1.290543e+01 4.256129e+02 4.256129e+02 -6.283143e+02 -3.527732e+02 366 367 368 369 370 4.256129e+02 4.256129e+02 4.256129e+02 4.256129e+02 -8.169858e+01 371 372 373 374 375 4.256129e+02 6.231531e+01 1.581195e+01 4.316527e+02 4.256129e+02 376 377 378 379 380 4.256129e+02 -4.910971e+02 2.418058e+02 1.563648e+03 4.256129e+02 381 382 383 384 385 4.256129e+02 -7.016803e+02 -4.731015e+02 4.256129e+02 -7.201314e+02 386 387 388 389 390 7.285816e+02 1.237071e+02 7.483456e+02 1.922723e+02 1.997467e+02 391 392 393 394 395 -2.387558e+02 4.384133e+02 6.893772e+02 1.603085e+02 3.598731e+03 396 397 398 399 400 4.048085e+02 3.156979e+03 -2.289213e+04 4.561749e+02 -1.621360e+03 401 402 403 404 405 1.444165e+03 2.737006e+02 5.543884e+02 2.549781e+02 6.989033e+02 406 407 408 409 410 5.398027e+02 2.096736e+01 2.393617e+02 4.558431e+01 1.747562e+03 411 412 413 414 415 -3.952910e+03 5.970922e+03 3.867466e+02 -5.064982e+02 2.871415e+03 416 417 418 419 420 2.396648e+03 4.204533e+01 -2.578826e+02 -2.808991e+02 1.687823e+03 421 422 423 424 425 2.044902e+03 -3.509178e+03 7.172237e+03 7.988975e+03 9.060725e+02 426 427 428 429 430 -1.941565e+03 4.521493e+03 1.257998e+04 7.247252e+03 6.549244e+03 431 1.366157e+04 > postscript(file="/var/www/html/rcomp/tmp/67fr01292938441.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 2.894692e+04 NA 1 -5.106920e+04 2.894692e+04 2 1.310643e+04 -5.106920e+04 3 4.022619e+04 1.310643e+04 4 1.317328e+04 4.022619e+04 5 1.563856e+04 1.317328e+04 6 5.705697e+03 1.563856e+04 7 6.654402e+03 5.705697e+03 8 8.587580e+03 6.654402e+03 9 -6.928696e+03 8.587580e+03 10 1.195252e+04 -6.928696e+03 11 -5.704364e+03 1.195252e+04 12 -3.300428e+04 -5.704364e+03 13 2.463621e+04 -3.300428e+04 14 -8.588906e+03 2.463621e+04 15 -4.683766e+03 -8.588906e+03 16 -1.932186e+03 -4.683766e+03 17 -8.925183e+03 -1.932186e+03 18 7.389960e+03 -8.925183e+03 19 -1.046652e+04 7.389960e+03 20 1.406486e+03 -1.046652e+04 21 1.262081e+04 1.406486e+03 22 1.469915e+03 1.262081e+04 23 -1.278850e+04 1.469915e+03 24 4.187309e+02 -1.278850e+04 25 8.101405e+03 4.187309e+02 26 8.794954e+02 8.101405e+03 27 5.227941e+03 8.794954e+02 28 7.765397e+02 5.227941e+03 29 -8.292800e+02 7.765397e+02 30 -1.766174e+03 -8.292800e+02 31 3.484866e+03 -1.766174e+03 32 -1.828464e+02 3.484866e+03 33 -1.198564e+03 -1.828464e+02 34 2.657120e+03 -1.198564e+03 35 -2.664355e+03 2.657120e+03 36 6.862752e+02 -2.664355e+03 37 -5.020128e+03 6.862752e+02 38 -5.613173e+03 -5.020128e+03 39 -7.755034e+03 -5.613173e+03 40 -1.063072e+04 -7.755034e+03 41 8.702920e+03 -1.063072e+04 42 -5.339650e+03 8.702920e+03 43 -1.393799e+03 -5.339650e+03 44 5.244980e+02 -1.393799e+03 45 1.690748e+03 5.244980e+02 46 2.358851e+03 1.690748e+03 47 -7.196577e+03 2.358851e+03 48 -9.787577e+02 -7.196577e+03 49 -1.784293e+02 -9.787577e+02 50 -6.817959e+03 -1.784293e+02 51 -1.001107e+04 -6.817959e+03 52 -1.661993e+03 -1.001107e+04 53 -5.357390e+03 -1.661993e+03 54 -4.309310e+03 -5.357390e+03 55 -5.949287e+03 -4.309310e+03 56 3.611094e+02 -5.949287e+03 57 2.387043e+02 3.611094e+02 58 -3.120623e+03 2.387043e+02 59 7.825609e+02 -3.120623e+03 60 5.404169e+02 7.825609e+02 61 1.451091e+02 5.404169e+02 62 6.684533e+02 1.451091e+02 63 -3.428489e+01 6.684533e+02 64 -2.867449e+03 -3.428489e+01 65 -1.212868e+04 -2.867449e+03 66 5.867092e+02 -1.212868e+04 67 -2.110575e+03 5.867092e+02 68 -5.006678e+03 -2.110575e+03 69 -4.617115e+03 -5.006678e+03 70 -1.164665e+03 -4.617115e+03 71 6.754499e+02 -1.164665e+03 72 -9.190005e+02 6.754499e+02 73 2.791772e+02 -9.190005e+02 74 -8.104282e+02 2.791772e+02 75 1.026007e+03 -8.104282e+02 76 -1.092054e+03 1.026007e+03 77 1.030780e+03 -1.092054e+03 78 -3.492224e+03 1.030780e+03 79 5.221564e+02 -3.492224e+03 80 -1.473887e+03 5.221564e+02 81 -2.200072e+03 -1.473887e+03 82 -1.085251e+03 -2.200072e+03 83 -5.573695e+02 -1.085251e+03 84 -4.359706e+02 -5.573695e+02 85 -1.226300e+03 -4.359706e+02 86 -7.034342e+03 -1.226300e+03 87 2.607924e+01 -7.034342e+03 88 -1.270691e+03 2.607924e+01 89 -3.463990e+02 -1.270691e+03 90 -1.386142e+03 -3.463990e+02 91 -9.391039e+02 -1.386142e+03 92 -5.525623e+03 -9.391039e+02 93 4.668803e+03 -5.525623e+03 94 -2.920440e+03 4.668803e+03 95 -1.761075e+02 -2.920440e+03 96 -3.395656e+02 -1.761075e+02 97 -4.059864e+02 -3.395656e+02 98 1.444362e+03 -4.059864e+02 99 1.072832e+03 1.444362e+03 100 -1.181124e+03 1.072832e+03 101 6.701451e+02 -1.181124e+03 102 -9.213749e+02 6.701451e+02 103 1.107648e+03 -9.213749e+02 104 -1.116078e+03 1.107648e+03 105 -1.230544e+03 -1.116078e+03 106 3.344967e+02 -1.230544e+03 107 -2.819851e+02 3.344967e+02 108 1.918069e+03 -2.819851e+02 109 1.086564e+03 1.918069e+03 110 -1.621986e+03 1.086564e+03 111 1.629884e+03 -1.621986e+03 112 1.655608e+03 1.629884e+03 113 1.940362e+03 1.655608e+03 114 5.796129e+02 1.940362e+03 115 -2.016522e+03 5.796129e+02 116 2.671847e+02 -2.016522e+03 117 2.533866e+02 2.671847e+02 118 2.711297e+03 2.533866e+02 119 -9.020219e+02 2.711297e+03 120 1.820060e+03 -9.020219e+02 121 3.296227e-01 1.820060e+03 122 -9.368362e+03 3.296227e-01 123 -1.158514e+03 -9.368362e+03 124 1.016908e+03 -1.158514e+03 125 2.480212e+03 1.016908e+03 126 4.256129e+02 2.480212e+03 127 -5.055148e+02 4.256129e+02 128 -1.602746e+03 -5.055148e+02 129 4.256129e+02 -1.602746e+03 130 -2.565545e+02 4.256129e+02 131 -1.172041e+03 -2.565545e+02 132 3.222175e+02 -1.172041e+03 133 -1.573782e+03 3.222175e+02 134 6.865223e+02 -1.573782e+03 135 3.276064e+02 6.865223e+02 136 -5.208681e+02 3.276064e+02 137 3.316166e+02 -5.208681e+02 138 9.131984e+02 3.316166e+02 139 4.256129e+02 9.131984e+02 140 -1.329631e+02 4.256129e+02 141 4.256129e+02 -1.329631e+02 142 3.713046e+02 4.256129e+02 143 4.700327e+01 3.713046e+02 144 -7.482228e+02 4.700327e+01 145 -2.463302e+02 -7.482228e+02 146 4.256129e+02 -2.463302e+02 147 -2.836501e+02 4.256129e+02 148 4.017596e+02 -2.836501e+02 149 4.403208e+02 4.017596e+02 150 3.976271e+02 4.403208e+02 151 -2.485410e+03 3.976271e+02 152 -9.136409e+01 -2.485410e+03 153 2.806277e+02 -9.136409e+01 154 2.438142e+02 2.806277e+02 155 4.256129e+02 2.438142e+02 156 4.256129e+02 4.256129e+02 157 3.314850e+02 4.256129e+02 158 1.280201e+02 3.314850e+02 159 4.256129e+02 1.280201e+02 160 4.256129e+02 4.256129e+02 161 4.256129e+02 4.256129e+02 162 4.256129e+02 4.256129e+02 163 2.680640e+02 4.256129e+02 164 -1.440006e+03 2.680640e+02 165 3.736612e+02 -1.440006e+03 166 4.256129e+02 3.736612e+02 167 -1.540466e+03 4.256129e+02 168 -8.398307e+02 -1.540466e+03 169 4.256129e+02 -8.398307e+02 170 4.256129e+02 4.256129e+02 171 1.925743e+02 4.256129e+02 172 4.256129e+02 1.925743e+02 173 -2.291726e+02 4.256129e+02 174 4.256129e+02 -2.291726e+02 175 -2.730270e+02 4.256129e+02 176 -1.083921e+04 -2.730270e+02 177 4.899444e+02 -1.083921e+04 178 4.256129e+02 4.899444e+02 179 4.256129e+02 4.256129e+02 180 -7.693573e+02 4.256129e+02 181 -1.220579e+02 -7.693573e+02 182 2.258441e+02 -1.220579e+02 183 -4.966188e+02 2.258441e+02 184 -3.136088e+01 -4.966188e+02 185 -2.389170e+02 -3.136088e+01 186 9.194398e+01 -2.389170e+02 187 -4.533803e+00 9.194398e+01 188 4.324641e+02 -4.533803e+00 189 4.256129e+02 4.324641e+02 190 -7.437756e+01 4.256129e+02 191 4.256129e+02 -7.437756e+01 192 -1.412196e+02 4.256129e+02 193 4.256129e+02 -1.412196e+02 194 1.436319e+02 4.256129e+02 195 4.256129e+02 1.436319e+02 196 -3.720324e+03 4.256129e+02 197 4.256129e+02 -3.720324e+03 198 4.256129e+02 4.256129e+02 199 -3.225290e+02 4.256129e+02 200 4.256129e+02 -3.225290e+02 201 4.256129e+02 4.256129e+02 202 4.256129e+02 4.256129e+02 203 -6.557587e+02 4.256129e+02 204 5.411853e+01 -6.557587e+02 205 -7.201745e+01 5.411853e+01 206 3.396357e+02 -7.201745e+01 207 4.256129e+02 3.396357e+02 208 -1.590523e+03 4.256129e+02 209 4.556666e+02 -1.590523e+03 210 4.256129e+02 4.556666e+02 211 4.256129e+02 4.256129e+02 212 -6.292313e+02 4.256129e+02 213 1.230245e+02 -6.292313e+02 214 5.728642e+02 1.230245e+02 215 2.188221e+02 5.728642e+02 216 -1.805793e+02 2.188221e+02 217 7.653652e+01 -1.805793e+02 218 2.188221e+02 7.653652e+01 219 -7.435791e+02 2.188221e+02 220 5.668551e+02 -7.435791e+02 221 4.256129e+02 5.668551e+02 222 3.250672e+02 4.256129e+02 223 -1.145571e+03 3.250672e+02 224 3.419453e+02 -1.145571e+03 225 -1.555115e+03 3.419453e+02 226 -4.365484e+02 -1.555115e+03 227 4.256129e+02 -4.365484e+02 228 5.171134e+02 4.256129e+02 229 4.256129e+02 5.171134e+02 230 2.231932e+02 4.256129e+02 231 1.795512e+02 2.231932e+02 232 4.256129e+02 1.795512e+02 233 -4.253277e+02 4.256129e+02 234 4.256129e+02 -4.253277e+02 235 4.270577e+02 4.256129e+02 236 4.256129e+02 4.270577e+02 237 -1.518905e+02 4.256129e+02 238 4.256129e+02 -1.518905e+02 239 4.256129e+02 4.256129e+02 240 4.256129e+02 4.256129e+02 241 2.427915e+02 4.256129e+02 242 4.256129e+02 2.427915e+02 243 -2.921195e+02 4.256129e+02 244 2.678262e+03 -2.921195e+02 245 4.256129e+02 2.678262e+03 246 4.256129e+02 4.256129e+02 247 -9.183872e+02 4.256129e+02 248 4.256129e+02 -9.183872e+02 249 4.256129e+02 4.256129e+02 250 -2.644955e+02 4.256129e+02 251 4.256129e+02 -2.644955e+02 252 6.223636e+02 4.256129e+02 253 -2.569747e+03 6.223636e+02 254 -2.206108e+03 -2.569747e+03 255 4.256129e+02 -2.206108e+03 256 1.627109e+01 4.256129e+02 257 -1.973864e+03 1.627109e+01 258 4.256129e+02 -1.973864e+03 259 4.256129e+02 4.256129e+02 260 4.256129e+02 4.256129e+02 261 4.256129e+02 4.256129e+02 262 -7.367878e+02 4.256129e+02 263 -4.816273e+02 -7.367878e+02 264 4.256129e+02 -4.816273e+02 265 4.256129e+02 4.256129e+02 266 4.256129e+02 4.256129e+02 267 -6.098594e+02 4.256129e+02 268 4.256129e+02 -6.098594e+02 269 -2.109743e+03 4.256129e+02 270 4.256129e+02 -2.109743e+03 271 4.256129e+02 4.256129e+02 272 4.256129e+02 4.256129e+02 273 -1.366331e+03 4.256129e+02 274 4.256129e+02 -1.366331e+03 275 -5.444077e+02 4.256129e+02 276 -1.602230e+03 -5.444077e+02 277 4.256129e+02 -1.602230e+03 278 -9.220581e+02 4.256129e+02 279 -5.467651e+01 -9.220581e+02 280 4.256129e+02 -5.467651e+01 281 4.256129e+02 4.256129e+02 282 -2.066203e+03 4.256129e+02 283 4.256129e+02 -2.066203e+03 284 -1.388861e+02 4.256129e+02 285 4.256129e+02 -1.388861e+02 286 4.754603e+02 4.256129e+02 287 1.441942e+02 4.754603e+02 288 4.256129e+02 1.441942e+02 289 -5.787738e+02 4.256129e+02 290 4.256129e+02 -5.787738e+02 291 -7.006926e+02 4.256129e+02 292 -1.182706e+03 -7.006926e+02 293 4.256129e+02 -1.182706e+03 294 -4.606910e+02 4.256129e+02 295 -3.748129e+02 -4.606910e+02 296 4.256129e+02 -3.748129e+02 297 4.256129e+02 4.256129e+02 298 4.256129e+02 4.256129e+02 299 -5.049457e+02 4.256129e+02 300 6.316016e+01 -5.049457e+02 301 7.055370e+02 6.316016e+01 302 3.222175e+02 7.055370e+02 303 4.256129e+02 3.222175e+02 304 -8.961924e+02 4.256129e+02 305 -8.129961e+02 -8.961924e+02 306 4.256129e+02 -8.129961e+02 307 -1.283940e+03 4.256129e+02 308 -1.003526e+03 -1.283940e+03 309 4.256129e+02 -1.003526e+03 310 4.256129e+02 4.256129e+02 311 -2.336793e+03 4.256129e+02 312 4.256129e+02 -2.336793e+03 313 -9.734246e+02 4.256129e+02 314 -4.003007e+02 -9.734246e+02 315 2.185694e+01 -4.003007e+02 316 4.486071e+02 2.185694e+01 317 -1.103780e+03 4.486071e+02 318 -9.031008e+01 -1.103780e+03 319 4.256129e+02 -9.031008e+01 320 4.256129e+02 4.256129e+02 321 4.256129e+02 4.256129e+02 322 4.256129e+02 4.256129e+02 323 4.256129e+02 4.256129e+02 324 4.256129e+02 4.256129e+02 325 3.699330e+02 4.256129e+02 326 4.256129e+02 3.699330e+02 327 4.256129e+02 4.256129e+02 328 4.256129e+02 4.256129e+02 329 4.256129e+02 4.256129e+02 330 4.256129e+02 4.256129e+02 331 4.256129e+02 4.256129e+02 332 4.256129e+02 4.256129e+02 333 1.168389e+02 4.256129e+02 334 3.232677e+02 1.168389e+02 335 3.261116e+02 3.232677e+02 336 -6.230431e+02 3.261116e+02 337 4.256129e+02 -6.230431e+02 338 3.646600e+02 4.256129e+02 339 -5.211557e+02 3.646600e+02 340 -2.953794e+02 -5.211557e+02 341 4.256129e+02 -2.953794e+02 342 4.256129e+02 4.256129e+02 343 4.256129e+02 4.256129e+02 344 -2.339587e+02 4.256129e+02 345 4.256129e+02 -2.339587e+02 346 3.110211e+02 4.256129e+02 347 4.256129e+02 3.110211e+02 348 4.256129e+02 4.256129e+02 349 9.510645e+02 4.256129e+02 350 3.610234e+02 9.510645e+02 351 -2.192591e+02 3.610234e+02 352 -1.329647e+02 -2.192591e+02 353 4.910256e+02 -1.329647e+02 354 4.256129e+02 4.910256e+02 355 4.256129e+02 4.256129e+02 356 -3.599873e+02 4.256129e+02 357 4.256129e+02 -3.599873e+02 358 4.256129e+02 4.256129e+02 359 4.256129e+02 4.256129e+02 360 1.290543e+01 4.256129e+02 361 4.256129e+02 1.290543e+01 362 4.256129e+02 4.256129e+02 363 -6.283143e+02 4.256129e+02 364 -3.527732e+02 -6.283143e+02 365 4.256129e+02 -3.527732e+02 366 4.256129e+02 4.256129e+02 367 4.256129e+02 4.256129e+02 368 4.256129e+02 4.256129e+02 369 -8.169858e+01 4.256129e+02 370 4.256129e+02 -8.169858e+01 371 6.231531e+01 4.256129e+02 372 1.581195e+01 6.231531e+01 373 4.316527e+02 1.581195e+01 374 4.256129e+02 4.316527e+02 375 4.256129e+02 4.256129e+02 376 -4.910971e+02 4.256129e+02 377 2.418058e+02 -4.910971e+02 378 1.563648e+03 2.418058e+02 379 4.256129e+02 1.563648e+03 380 4.256129e+02 4.256129e+02 381 -7.016803e+02 4.256129e+02 382 -4.731015e+02 -7.016803e+02 383 4.256129e+02 -4.731015e+02 384 -7.201314e+02 4.256129e+02 385 7.285816e+02 -7.201314e+02 386 1.237071e+02 7.285816e+02 387 7.483456e+02 1.237071e+02 388 1.922723e+02 7.483456e+02 389 1.997467e+02 1.922723e+02 390 -2.387558e+02 1.997467e+02 391 4.384133e+02 -2.387558e+02 392 6.893772e+02 4.384133e+02 393 1.603085e+02 6.893772e+02 394 3.598731e+03 1.603085e+02 395 4.048085e+02 3.598731e+03 396 3.156979e+03 4.048085e+02 397 -2.289213e+04 3.156979e+03 398 4.561749e+02 -2.289213e+04 399 -1.621360e+03 4.561749e+02 400 1.444165e+03 -1.621360e+03 401 2.737006e+02 1.444165e+03 402 5.543884e+02 2.737006e+02 403 2.549781e+02 5.543884e+02 404 6.989033e+02 2.549781e+02 405 5.398027e+02 6.989033e+02 406 2.096736e+01 5.398027e+02 407 2.393617e+02 2.096736e+01 408 4.558431e+01 2.393617e+02 409 1.747562e+03 4.558431e+01 410 -3.952910e+03 1.747562e+03 411 5.970922e+03 -3.952910e+03 412 3.867466e+02 5.970922e+03 413 -5.064982e+02 3.867466e+02 414 2.871415e+03 -5.064982e+02 415 2.396648e+03 2.871415e+03 416 4.204533e+01 2.396648e+03 417 -2.578826e+02 4.204533e+01 418 -2.808991e+02 -2.578826e+02 419 1.687823e+03 -2.808991e+02 420 2.044902e+03 1.687823e+03 421 -3.509178e+03 2.044902e+03 422 7.172237e+03 -3.509178e+03 423 7.988975e+03 7.172237e+03 424 9.060725e+02 7.988975e+03 425 -1.941565e+03 9.060725e+02 426 4.521493e+03 -1.941565e+03 427 1.257998e+04 4.521493e+03 428 7.247252e+03 1.257998e+04 429 6.549244e+03 7.247252e+03 430 1.366157e+04 6.549244e+03 431 NA 1.366157e+04 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -5.106920e+04 2.894692e+04 [2,] 1.310643e+04 -5.106920e+04 [3,] 4.022619e+04 1.310643e+04 [4,] 1.317328e+04 4.022619e+04 [5,] 1.563856e+04 1.317328e+04 [6,] 5.705697e+03 1.563856e+04 [7,] 6.654402e+03 5.705697e+03 [8,] 8.587580e+03 6.654402e+03 [9,] -6.928696e+03 8.587580e+03 [10,] 1.195252e+04 -6.928696e+03 [11,] -5.704364e+03 1.195252e+04 [12,] -3.300428e+04 -5.704364e+03 [13,] 2.463621e+04 -3.300428e+04 [14,] -8.588906e+03 2.463621e+04 [15,] -4.683766e+03 -8.588906e+03 [16,] -1.932186e+03 -4.683766e+03 [17,] -8.925183e+03 -1.932186e+03 [18,] 7.389960e+03 -8.925183e+03 [19,] -1.046652e+04 7.389960e+03 [20,] 1.406486e+03 -1.046652e+04 [21,] 1.262081e+04 1.406486e+03 [22,] 1.469915e+03 1.262081e+04 [23,] -1.278850e+04 1.469915e+03 [24,] 4.187309e+02 -1.278850e+04 [25,] 8.101405e+03 4.187309e+02 [26,] 8.794954e+02 8.101405e+03 [27,] 5.227941e+03 8.794954e+02 [28,] 7.765397e+02 5.227941e+03 [29,] -8.292800e+02 7.765397e+02 [30,] -1.766174e+03 -8.292800e+02 [31,] 3.484866e+03 -1.766174e+03 [32,] -1.828464e+02 3.484866e+03 [33,] -1.198564e+03 -1.828464e+02 [34,] 2.657120e+03 -1.198564e+03 [35,] -2.664355e+03 2.657120e+03 [36,] 6.862752e+02 -2.664355e+03 [37,] -5.020128e+03 6.862752e+02 [38,] -5.613173e+03 -5.020128e+03 [39,] -7.755034e+03 -5.613173e+03 [40,] -1.063072e+04 -7.755034e+03 [41,] 8.702920e+03 -1.063072e+04 [42,] -5.339650e+03 8.702920e+03 [43,] -1.393799e+03 -5.339650e+03 [44,] 5.244980e+02 -1.393799e+03 [45,] 1.690748e+03 5.244980e+02 [46,] 2.358851e+03 1.690748e+03 [47,] -7.196577e+03 2.358851e+03 [48,] -9.787577e+02 -7.196577e+03 [49,] -1.784293e+02 -9.787577e+02 [50,] -6.817959e+03 -1.784293e+02 [51,] -1.001107e+04 -6.817959e+03 [52,] -1.661993e+03 -1.001107e+04 [53,] -5.357390e+03 -1.661993e+03 [54,] -4.309310e+03 -5.357390e+03 [55,] -5.949287e+03 -4.309310e+03 [56,] 3.611094e+02 -5.949287e+03 [57,] 2.387043e+02 3.611094e+02 [58,] -3.120623e+03 2.387043e+02 [59,] 7.825609e+02 -3.120623e+03 [60,] 5.404169e+02 7.825609e+02 [61,] 1.451091e+02 5.404169e+02 [62,] 6.684533e+02 1.451091e+02 [63,] -3.428489e+01 6.684533e+02 [64,] -2.867449e+03 -3.428489e+01 [65,] -1.212868e+04 -2.867449e+03 [66,] 5.867092e+02 -1.212868e+04 [67,] -2.110575e+03 5.867092e+02 [68,] -5.006678e+03 -2.110575e+03 [69,] -4.617115e+03 -5.006678e+03 [70,] -1.164665e+03 -4.617115e+03 [71,] 6.754499e+02 -1.164665e+03 [72,] -9.190005e+02 6.754499e+02 [73,] 2.791772e+02 -9.190005e+02 [74,] -8.104282e+02 2.791772e+02 [75,] 1.026007e+03 -8.104282e+02 [76,] -1.092054e+03 1.026007e+03 [77,] 1.030780e+03 -1.092054e+03 [78,] -3.492224e+03 1.030780e+03 [79,] 5.221564e+02 -3.492224e+03 [80,] -1.473887e+03 5.221564e+02 [81,] -2.200072e+03 -1.473887e+03 [82,] -1.085251e+03 -2.200072e+03 [83,] -5.573695e+02 -1.085251e+03 [84,] -4.359706e+02 -5.573695e+02 [85,] -1.226300e+03 -4.359706e+02 [86,] -7.034342e+03 -1.226300e+03 [87,] 2.607924e+01 -7.034342e+03 [88,] -1.270691e+03 2.607924e+01 [89,] -3.463990e+02 -1.270691e+03 [90,] -1.386142e+03 -3.463990e+02 [91,] -9.391039e+02 -1.386142e+03 [92,] -5.525623e+03 -9.391039e+02 [93,] 4.668803e+03 -5.525623e+03 [94,] -2.920440e+03 4.668803e+03 [95,] -1.761075e+02 -2.920440e+03 [96,] -3.395656e+02 -1.761075e+02 [97,] -4.059864e+02 -3.395656e+02 [98,] 1.444362e+03 -4.059864e+02 [99,] 1.072832e+03 1.444362e+03 [100,] -1.181124e+03 1.072832e+03 [101,] 6.701451e+02 -1.181124e+03 [102,] -9.213749e+02 6.701451e+02 [103,] 1.107648e+03 -9.213749e+02 [104,] -1.116078e+03 1.107648e+03 [105,] -1.230544e+03 -1.116078e+03 [106,] 3.344967e+02 -1.230544e+03 [107,] -2.819851e+02 3.344967e+02 [108,] 1.918069e+03 -2.819851e+02 [109,] 1.086564e+03 1.918069e+03 [110,] -1.621986e+03 1.086564e+03 [111,] 1.629884e+03 -1.621986e+03 [112,] 1.655608e+03 1.629884e+03 [113,] 1.940362e+03 1.655608e+03 [114,] 5.796129e+02 1.940362e+03 [115,] -2.016522e+03 5.796129e+02 [116,] 2.671847e+02 -2.016522e+03 [117,] 2.533866e+02 2.671847e+02 [118,] 2.711297e+03 2.533866e+02 [119,] -9.020219e+02 2.711297e+03 [120,] 1.820060e+03 -9.020219e+02 [121,] 3.296227e-01 1.820060e+03 [122,] -9.368362e+03 3.296227e-01 [123,] -1.158514e+03 -9.368362e+03 [124,] 1.016908e+03 -1.158514e+03 [125,] 2.480212e+03 1.016908e+03 [126,] 4.256129e+02 2.480212e+03 [127,] -5.055148e+02 4.256129e+02 [128,] -1.602746e+03 -5.055148e+02 [129,] 4.256129e+02 -1.602746e+03 [130,] -2.565545e+02 4.256129e+02 [131,] -1.172041e+03 -2.565545e+02 [132,] 3.222175e+02 -1.172041e+03 [133,] -1.573782e+03 3.222175e+02 [134,] 6.865223e+02 -1.573782e+03 [135,] 3.276064e+02 6.865223e+02 [136,] -5.208681e+02 3.276064e+02 [137,] 3.316166e+02 -5.208681e+02 [138,] 9.131984e+02 3.316166e+02 [139,] 4.256129e+02 9.131984e+02 [140,] -1.329631e+02 4.256129e+02 [141,] 4.256129e+02 -1.329631e+02 [142,] 3.713046e+02 4.256129e+02 [143,] 4.700327e+01 3.713046e+02 [144,] -7.482228e+02 4.700327e+01 [145,] -2.463302e+02 -7.482228e+02 [146,] 4.256129e+02 -2.463302e+02 [147,] -2.836501e+02 4.256129e+02 [148,] 4.017596e+02 -2.836501e+02 [149,] 4.403208e+02 4.017596e+02 [150,] 3.976271e+02 4.403208e+02 [151,] -2.485410e+03 3.976271e+02 [152,] -9.136409e+01 -2.485410e+03 [153,] 2.806277e+02 -9.136409e+01 [154,] 2.438142e+02 2.806277e+02 [155,] 4.256129e+02 2.438142e+02 [156,] 4.256129e+02 4.256129e+02 [157,] 3.314850e+02 4.256129e+02 [158,] 1.280201e+02 3.314850e+02 [159,] 4.256129e+02 1.280201e+02 [160,] 4.256129e+02 4.256129e+02 [161,] 4.256129e+02 4.256129e+02 [162,] 4.256129e+02 4.256129e+02 [163,] 2.680640e+02 4.256129e+02 [164,] -1.440006e+03 2.680640e+02 [165,] 3.736612e+02 -1.440006e+03 [166,] 4.256129e+02 3.736612e+02 [167,] -1.540466e+03 4.256129e+02 [168,] -8.398307e+02 -1.540466e+03 [169,] 4.256129e+02 -8.398307e+02 [170,] 4.256129e+02 4.256129e+02 [171,] 1.925743e+02 4.256129e+02 [172,] 4.256129e+02 1.925743e+02 [173,] -2.291726e+02 4.256129e+02 [174,] 4.256129e+02 -2.291726e+02 [175,] -2.730270e+02 4.256129e+02 [176,] -1.083921e+04 -2.730270e+02 [177,] 4.899444e+02 -1.083921e+04 [178,] 4.256129e+02 4.899444e+02 [179,] 4.256129e+02 4.256129e+02 [180,] -7.693573e+02 4.256129e+02 [181,] -1.220579e+02 -7.693573e+02 [182,] 2.258441e+02 -1.220579e+02 [183,] -4.966188e+02 2.258441e+02 [184,] -3.136088e+01 -4.966188e+02 [185,] -2.389170e+02 -3.136088e+01 [186,] 9.194398e+01 -2.389170e+02 [187,] -4.533803e+00 9.194398e+01 [188,] 4.324641e+02 -4.533803e+00 [189,] 4.256129e+02 4.324641e+02 [190,] -7.437756e+01 4.256129e+02 [191,] 4.256129e+02 -7.437756e+01 [192,] -1.412196e+02 4.256129e+02 [193,] 4.256129e+02 -1.412196e+02 [194,] 1.436319e+02 4.256129e+02 [195,] 4.256129e+02 1.436319e+02 [196,] -3.720324e+03 4.256129e+02 [197,] 4.256129e+02 -3.720324e+03 [198,] 4.256129e+02 4.256129e+02 [199,] -3.225290e+02 4.256129e+02 [200,] 4.256129e+02 -3.225290e+02 [201,] 4.256129e+02 4.256129e+02 [202,] 4.256129e+02 4.256129e+02 [203,] -6.557587e+02 4.256129e+02 [204,] 5.411853e+01 -6.557587e+02 [205,] -7.201745e+01 5.411853e+01 [206,] 3.396357e+02 -7.201745e+01 [207,] 4.256129e+02 3.396357e+02 [208,] -1.590523e+03 4.256129e+02 [209,] 4.556666e+02 -1.590523e+03 [210,] 4.256129e+02 4.556666e+02 [211,] 4.256129e+02 4.256129e+02 [212,] -6.292313e+02 4.256129e+02 [213,] 1.230245e+02 -6.292313e+02 [214,] 5.728642e+02 1.230245e+02 [215,] 2.188221e+02 5.728642e+02 [216,] -1.805793e+02 2.188221e+02 [217,] 7.653652e+01 -1.805793e+02 [218,] 2.188221e+02 7.653652e+01 [219,] -7.435791e+02 2.188221e+02 [220,] 5.668551e+02 -7.435791e+02 [221,] 4.256129e+02 5.668551e+02 [222,] 3.250672e+02 4.256129e+02 [223,] -1.145571e+03 3.250672e+02 [224,] 3.419453e+02 -1.145571e+03 [225,] -1.555115e+03 3.419453e+02 [226,] -4.365484e+02 -1.555115e+03 [227,] 4.256129e+02 -4.365484e+02 [228,] 5.171134e+02 4.256129e+02 [229,] 4.256129e+02 5.171134e+02 [230,] 2.231932e+02 4.256129e+02 [231,] 1.795512e+02 2.231932e+02 [232,] 4.256129e+02 1.795512e+02 [233,] -4.253277e+02 4.256129e+02 [234,] 4.256129e+02 -4.253277e+02 [235,] 4.270577e+02 4.256129e+02 [236,] 4.256129e+02 4.270577e+02 [237,] -1.518905e+02 4.256129e+02 [238,] 4.256129e+02 -1.518905e+02 [239,] 4.256129e+02 4.256129e+02 [240,] 4.256129e+02 4.256129e+02 [241,] 2.427915e+02 4.256129e+02 [242,] 4.256129e+02 2.427915e+02 [243,] -2.921195e+02 4.256129e+02 [244,] 2.678262e+03 -2.921195e+02 [245,] 4.256129e+02 2.678262e+03 [246,] 4.256129e+02 4.256129e+02 [247,] -9.183872e+02 4.256129e+02 [248,] 4.256129e+02 -9.183872e+02 [249,] 4.256129e+02 4.256129e+02 [250,] -2.644955e+02 4.256129e+02 [251,] 4.256129e+02 -2.644955e+02 [252,] 6.223636e+02 4.256129e+02 [253,] -2.569747e+03 6.223636e+02 [254,] -2.206108e+03 -2.569747e+03 [255,] 4.256129e+02 -2.206108e+03 [256,] 1.627109e+01 4.256129e+02 [257,] -1.973864e+03 1.627109e+01 [258,] 4.256129e+02 -1.973864e+03 [259,] 4.256129e+02 4.256129e+02 [260,] 4.256129e+02 4.256129e+02 [261,] 4.256129e+02 4.256129e+02 [262,] -7.367878e+02 4.256129e+02 [263,] -4.816273e+02 -7.367878e+02 [264,] 4.256129e+02 -4.816273e+02 [265,] 4.256129e+02 4.256129e+02 [266,] 4.256129e+02 4.256129e+02 [267,] -6.098594e+02 4.256129e+02 [268,] 4.256129e+02 -6.098594e+02 [269,] -2.109743e+03 4.256129e+02 [270,] 4.256129e+02 -2.109743e+03 [271,] 4.256129e+02 4.256129e+02 [272,] 4.256129e+02 4.256129e+02 [273,] -1.366331e+03 4.256129e+02 [274,] 4.256129e+02 -1.366331e+03 [275,] -5.444077e+02 4.256129e+02 [276,] -1.602230e+03 -5.444077e+02 [277,] 4.256129e+02 -1.602230e+03 [278,] -9.220581e+02 4.256129e+02 [279,] -5.467651e+01 -9.220581e+02 [280,] 4.256129e+02 -5.467651e+01 [281,] 4.256129e+02 4.256129e+02 [282,] -2.066203e+03 4.256129e+02 [283,] 4.256129e+02 -2.066203e+03 [284,] -1.388861e+02 4.256129e+02 [285,] 4.256129e+02 -1.388861e+02 [286,] 4.754603e+02 4.256129e+02 [287,] 1.441942e+02 4.754603e+02 [288,] 4.256129e+02 1.441942e+02 [289,] -5.787738e+02 4.256129e+02 [290,] 4.256129e+02 -5.787738e+02 [291,] -7.006926e+02 4.256129e+02 [292,] -1.182706e+03 -7.006926e+02 [293,] 4.256129e+02 -1.182706e+03 [294,] -4.606910e+02 4.256129e+02 [295,] -3.748129e+02 -4.606910e+02 [296,] 4.256129e+02 -3.748129e+02 [297,] 4.256129e+02 4.256129e+02 [298,] 4.256129e+02 4.256129e+02 [299,] -5.049457e+02 4.256129e+02 [300,] 6.316016e+01 -5.049457e+02 [301,] 7.055370e+02 6.316016e+01 [302,] 3.222175e+02 7.055370e+02 [303,] 4.256129e+02 3.222175e+02 [304,] -8.961924e+02 4.256129e+02 [305,] -8.129961e+02 -8.961924e+02 [306,] 4.256129e+02 -8.129961e+02 [307,] -1.283940e+03 4.256129e+02 [308,] -1.003526e+03 -1.283940e+03 [309,] 4.256129e+02 -1.003526e+03 [310,] 4.256129e+02 4.256129e+02 [311,] -2.336793e+03 4.256129e+02 [312,] 4.256129e+02 -2.336793e+03 [313,] -9.734246e+02 4.256129e+02 [314,] -4.003007e+02 -9.734246e+02 [315,] 2.185694e+01 -4.003007e+02 [316,] 4.486071e+02 2.185694e+01 [317,] -1.103780e+03 4.486071e+02 [318,] -9.031008e+01 -1.103780e+03 [319,] 4.256129e+02 -9.031008e+01 [320,] 4.256129e+02 4.256129e+02 [321,] 4.256129e+02 4.256129e+02 [322,] 4.256129e+02 4.256129e+02 [323,] 4.256129e+02 4.256129e+02 [324,] 4.256129e+02 4.256129e+02 [325,] 3.699330e+02 4.256129e+02 [326,] 4.256129e+02 3.699330e+02 [327,] 4.256129e+02 4.256129e+02 [328,] 4.256129e+02 4.256129e+02 [329,] 4.256129e+02 4.256129e+02 [330,] 4.256129e+02 4.256129e+02 [331,] 4.256129e+02 4.256129e+02 [332,] 4.256129e+02 4.256129e+02 [333,] 1.168389e+02 4.256129e+02 [334,] 3.232677e+02 1.168389e+02 [335,] 3.261116e+02 3.232677e+02 [336,] -6.230431e+02 3.261116e+02 [337,] 4.256129e+02 -6.230431e+02 [338,] 3.646600e+02 4.256129e+02 [339,] -5.211557e+02 3.646600e+02 [340,] -2.953794e+02 -5.211557e+02 [341,] 4.256129e+02 -2.953794e+02 [342,] 4.256129e+02 4.256129e+02 [343,] 4.256129e+02 4.256129e+02 [344,] -2.339587e+02 4.256129e+02 [345,] 4.256129e+02 -2.339587e+02 [346,] 3.110211e+02 4.256129e+02 [347,] 4.256129e+02 3.110211e+02 [348,] 4.256129e+02 4.256129e+02 [349,] 9.510645e+02 4.256129e+02 [350,] 3.610234e+02 9.510645e+02 [351,] -2.192591e+02 3.610234e+02 [352,] -1.329647e+02 -2.192591e+02 [353,] 4.910256e+02 -1.329647e+02 [354,] 4.256129e+02 4.910256e+02 [355,] 4.256129e+02 4.256129e+02 [356,] -3.599873e+02 4.256129e+02 [357,] 4.256129e+02 -3.599873e+02 [358,] 4.256129e+02 4.256129e+02 [359,] 4.256129e+02 4.256129e+02 [360,] 1.290543e+01 4.256129e+02 [361,] 4.256129e+02 1.290543e+01 [362,] 4.256129e+02 4.256129e+02 [363,] -6.283143e+02 4.256129e+02 [364,] -3.527732e+02 -6.283143e+02 [365,] 4.256129e+02 -3.527732e+02 [366,] 4.256129e+02 4.256129e+02 [367,] 4.256129e+02 4.256129e+02 [368,] 4.256129e+02 4.256129e+02 [369,] -8.169858e+01 4.256129e+02 [370,] 4.256129e+02 -8.169858e+01 [371,] 6.231531e+01 4.256129e+02 [372,] 1.581195e+01 6.231531e+01 [373,] 4.316527e+02 1.581195e+01 [374,] 4.256129e+02 4.316527e+02 [375,] 4.256129e+02 4.256129e+02 [376,] -4.910971e+02 4.256129e+02 [377,] 2.418058e+02 -4.910971e+02 [378,] 1.563648e+03 2.418058e+02 [379,] 4.256129e+02 1.563648e+03 [380,] 4.256129e+02 4.256129e+02 [381,] -7.016803e+02 4.256129e+02 [382,] -4.731015e+02 -7.016803e+02 [383,] 4.256129e+02 -4.731015e+02 [384,] -7.201314e+02 4.256129e+02 [385,] 7.285816e+02 -7.201314e+02 [386,] 1.237071e+02 7.285816e+02 [387,] 7.483456e+02 1.237071e+02 [388,] 1.922723e+02 7.483456e+02 [389,] 1.997467e+02 1.922723e+02 [390,] -2.387558e+02 1.997467e+02 [391,] 4.384133e+02 -2.387558e+02 [392,] 6.893772e+02 4.384133e+02 [393,] 1.603085e+02 6.893772e+02 [394,] 3.598731e+03 1.603085e+02 [395,] 4.048085e+02 3.598731e+03 [396,] 3.156979e+03 4.048085e+02 [397,] -2.289213e+04 3.156979e+03 [398,] 4.561749e+02 -2.289213e+04 [399,] -1.621360e+03 4.561749e+02 [400,] 1.444165e+03 -1.621360e+03 [401,] 2.737006e+02 1.444165e+03 [402,] 5.543884e+02 2.737006e+02 [403,] 2.549781e+02 5.543884e+02 [404,] 6.989033e+02 2.549781e+02 [405,] 5.398027e+02 6.989033e+02 [406,] 2.096736e+01 5.398027e+02 [407,] 2.393617e+02 2.096736e+01 [408,] 4.558431e+01 2.393617e+02 [409,] 1.747562e+03 4.558431e+01 [410,] -3.952910e+03 1.747562e+03 [411,] 5.970922e+03 -3.952910e+03 [412,] 3.867466e+02 5.970922e+03 [413,] -5.064982e+02 3.867466e+02 [414,] 2.871415e+03 -5.064982e+02 [415,] 2.396648e+03 2.871415e+03 [416,] 4.204533e+01 2.396648e+03 [417,] -2.578826e+02 4.204533e+01 [418,] -2.808991e+02 -2.578826e+02 [419,] 1.687823e+03 -2.808991e+02 [420,] 2.044902e+03 1.687823e+03 [421,] -3.509178e+03 2.044902e+03 [422,] 7.172237e+03 -3.509178e+03 [423,] 7.988975e+03 7.172237e+03 [424,] 9.060725e+02 7.988975e+03 [425,] -1.941565e+03 9.060725e+02 [426,] 4.521493e+03 -1.941565e+03 [427,] 1.257998e+04 4.521493e+03 [428,] 7.247252e+03 1.257998e+04 [429,] 6.549244e+03 7.247252e+03 [430,] 1.366157e+04 6.549244e+03 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -5.106920e+04 2.894692e+04 2 1.310643e+04 -5.106920e+04 3 4.022619e+04 1.310643e+04 4 1.317328e+04 4.022619e+04 5 1.563856e+04 1.317328e+04 6 5.705697e+03 1.563856e+04 7 6.654402e+03 5.705697e+03 8 8.587580e+03 6.654402e+03 9 -6.928696e+03 8.587580e+03 10 1.195252e+04 -6.928696e+03 11 -5.704364e+03 1.195252e+04 12 -3.300428e+04 -5.704364e+03 13 2.463621e+04 -3.300428e+04 14 -8.588906e+03 2.463621e+04 15 -4.683766e+03 -8.588906e+03 16 -1.932186e+03 -4.683766e+03 17 -8.925183e+03 -1.932186e+03 18 7.389960e+03 -8.925183e+03 19 -1.046652e+04 7.389960e+03 20 1.406486e+03 -1.046652e+04 21 1.262081e+04 1.406486e+03 22 1.469915e+03 1.262081e+04 23 -1.278850e+04 1.469915e+03 24 4.187309e+02 -1.278850e+04 25 8.101405e+03 4.187309e+02 26 8.794954e+02 8.101405e+03 27 5.227941e+03 8.794954e+02 28 7.765397e+02 5.227941e+03 29 -8.292800e+02 7.765397e+02 30 -1.766174e+03 -8.292800e+02 31 3.484866e+03 -1.766174e+03 32 -1.828464e+02 3.484866e+03 33 -1.198564e+03 -1.828464e+02 34 2.657120e+03 -1.198564e+03 35 -2.664355e+03 2.657120e+03 36 6.862752e+02 -2.664355e+03 37 -5.020128e+03 6.862752e+02 38 -5.613173e+03 -5.020128e+03 39 -7.755034e+03 -5.613173e+03 40 -1.063072e+04 -7.755034e+03 41 8.702920e+03 -1.063072e+04 42 -5.339650e+03 8.702920e+03 43 -1.393799e+03 -5.339650e+03 44 5.244980e+02 -1.393799e+03 45 1.690748e+03 5.244980e+02 46 2.358851e+03 1.690748e+03 47 -7.196577e+03 2.358851e+03 48 -9.787577e+02 -7.196577e+03 49 -1.784293e+02 -9.787577e+02 50 -6.817959e+03 -1.784293e+02 51 -1.001107e+04 -6.817959e+03 52 -1.661993e+03 -1.001107e+04 53 -5.357390e+03 -1.661993e+03 54 -4.309310e+03 -5.357390e+03 55 -5.949287e+03 -4.309310e+03 56 3.611094e+02 -5.949287e+03 57 2.387043e+02 3.611094e+02 58 -3.120623e+03 2.387043e+02 59 7.825609e+02 -3.120623e+03 60 5.404169e+02 7.825609e+02 61 1.451091e+02 5.404169e+02 62 6.684533e+02 1.451091e+02 63 -3.428489e+01 6.684533e+02 64 -2.867449e+03 -3.428489e+01 65 -1.212868e+04 -2.867449e+03 66 5.867092e+02 -1.212868e+04 67 -2.110575e+03 5.867092e+02 68 -5.006678e+03 -2.110575e+03 69 -4.617115e+03 -5.006678e+03 70 -1.164665e+03 -4.617115e+03 71 6.754499e+02 -1.164665e+03 72 -9.190005e+02 6.754499e+02 73 2.791772e+02 -9.190005e+02 74 -8.104282e+02 2.791772e+02 75 1.026007e+03 -8.104282e+02 76 -1.092054e+03 1.026007e+03 77 1.030780e+03 -1.092054e+03 78 -3.492224e+03 1.030780e+03 79 5.221564e+02 -3.492224e+03 80 -1.473887e+03 5.221564e+02 81 -2.200072e+03 -1.473887e+03 82 -1.085251e+03 -2.200072e+03 83 -5.573695e+02 -1.085251e+03 84 -4.359706e+02 -5.573695e+02 85 -1.226300e+03 -4.359706e+02 86 -7.034342e+03 -1.226300e+03 87 2.607924e+01 -7.034342e+03 88 -1.270691e+03 2.607924e+01 89 -3.463990e+02 -1.270691e+03 90 -1.386142e+03 -3.463990e+02 91 -9.391039e+02 -1.386142e+03 92 -5.525623e+03 -9.391039e+02 93 4.668803e+03 -5.525623e+03 94 -2.920440e+03 4.668803e+03 95 -1.761075e+02 -2.920440e+03 96 -3.395656e+02 -1.761075e+02 97 -4.059864e+02 -3.395656e+02 98 1.444362e+03 -4.059864e+02 99 1.072832e+03 1.444362e+03 100 -1.181124e+03 1.072832e+03 101 6.701451e+02 -1.181124e+03 102 -9.213749e+02 6.701451e+02 103 1.107648e+03 -9.213749e+02 104 -1.116078e+03 1.107648e+03 105 -1.230544e+03 -1.116078e+03 106 3.344967e+02 -1.230544e+03 107 -2.819851e+02 3.344967e+02 108 1.918069e+03 -2.819851e+02 109 1.086564e+03 1.918069e+03 110 -1.621986e+03 1.086564e+03 111 1.629884e+03 -1.621986e+03 112 1.655608e+03 1.629884e+03 113 1.940362e+03 1.655608e+03 114 5.796129e+02 1.940362e+03 115 -2.016522e+03 5.796129e+02 116 2.671847e+02 -2.016522e+03 117 2.533866e+02 2.671847e+02 118 2.711297e+03 2.533866e+02 119 -9.020219e+02 2.711297e+03 120 1.820060e+03 -9.020219e+02 121 3.296227e-01 1.820060e+03 122 -9.368362e+03 3.296227e-01 123 -1.158514e+03 -9.368362e+03 124 1.016908e+03 -1.158514e+03 125 2.480212e+03 1.016908e+03 126 4.256129e+02 2.480212e+03 127 -5.055148e+02 4.256129e+02 128 -1.602746e+03 -5.055148e+02 129 4.256129e+02 -1.602746e+03 130 -2.565545e+02 4.256129e+02 131 -1.172041e+03 -2.565545e+02 132 3.222175e+02 -1.172041e+03 133 -1.573782e+03 3.222175e+02 134 6.865223e+02 -1.573782e+03 135 3.276064e+02 6.865223e+02 136 -5.208681e+02 3.276064e+02 137 3.316166e+02 -5.208681e+02 138 9.131984e+02 3.316166e+02 139 4.256129e+02 9.131984e+02 140 -1.329631e+02 4.256129e+02 141 4.256129e+02 -1.329631e+02 142 3.713046e+02 4.256129e+02 143 4.700327e+01 3.713046e+02 144 -7.482228e+02 4.700327e+01 145 -2.463302e+02 -7.482228e+02 146 4.256129e+02 -2.463302e+02 147 -2.836501e+02 4.256129e+02 148 4.017596e+02 -2.836501e+02 149 4.403208e+02 4.017596e+02 150 3.976271e+02 4.403208e+02 151 -2.485410e+03 3.976271e+02 152 -9.136409e+01 -2.485410e+03 153 2.806277e+02 -9.136409e+01 154 2.438142e+02 2.806277e+02 155 4.256129e+02 2.438142e+02 156 4.256129e+02 4.256129e+02 157 3.314850e+02 4.256129e+02 158 1.280201e+02 3.314850e+02 159 4.256129e+02 1.280201e+02 160 4.256129e+02 4.256129e+02 161 4.256129e+02 4.256129e+02 162 4.256129e+02 4.256129e+02 163 2.680640e+02 4.256129e+02 164 -1.440006e+03 2.680640e+02 165 3.736612e+02 -1.440006e+03 166 4.256129e+02 3.736612e+02 167 -1.540466e+03 4.256129e+02 168 -8.398307e+02 -1.540466e+03 169 4.256129e+02 -8.398307e+02 170 4.256129e+02 4.256129e+02 171 1.925743e+02 4.256129e+02 172 4.256129e+02 1.925743e+02 173 -2.291726e+02 4.256129e+02 174 4.256129e+02 -2.291726e+02 175 -2.730270e+02 4.256129e+02 176 -1.083921e+04 -2.730270e+02 177 4.899444e+02 -1.083921e+04 178 4.256129e+02 4.899444e+02 179 4.256129e+02 4.256129e+02 180 -7.693573e+02 4.256129e+02 181 -1.220579e+02 -7.693573e+02 182 2.258441e+02 -1.220579e+02 183 -4.966188e+02 2.258441e+02 184 -3.136088e+01 -4.966188e+02 185 -2.389170e+02 -3.136088e+01 186 9.194398e+01 -2.389170e+02 187 -4.533803e+00 9.194398e+01 188 4.324641e+02 -4.533803e+00 189 4.256129e+02 4.324641e+02 190 -7.437756e+01 4.256129e+02 191 4.256129e+02 -7.437756e+01 192 -1.412196e+02 4.256129e+02 193 4.256129e+02 -1.412196e+02 194 1.436319e+02 4.256129e+02 195 4.256129e+02 1.436319e+02 196 -3.720324e+03 4.256129e+02 197 4.256129e+02 -3.720324e+03 198 4.256129e+02 4.256129e+02 199 -3.225290e+02 4.256129e+02 200 4.256129e+02 -3.225290e+02 201 4.256129e+02 4.256129e+02 202 4.256129e+02 4.256129e+02 203 -6.557587e+02 4.256129e+02 204 5.411853e+01 -6.557587e+02 205 -7.201745e+01 5.411853e+01 206 3.396357e+02 -7.201745e+01 207 4.256129e+02 3.396357e+02 208 -1.590523e+03 4.256129e+02 209 4.556666e+02 -1.590523e+03 210 4.256129e+02 4.556666e+02 211 4.256129e+02 4.256129e+02 212 -6.292313e+02 4.256129e+02 213 1.230245e+02 -6.292313e+02 214 5.728642e+02 1.230245e+02 215 2.188221e+02 5.728642e+02 216 -1.805793e+02 2.188221e+02 217 7.653652e+01 -1.805793e+02 218 2.188221e+02 7.653652e+01 219 -7.435791e+02 2.188221e+02 220 5.668551e+02 -7.435791e+02 221 4.256129e+02 5.668551e+02 222 3.250672e+02 4.256129e+02 223 -1.145571e+03 3.250672e+02 224 3.419453e+02 -1.145571e+03 225 -1.555115e+03 3.419453e+02 226 -4.365484e+02 -1.555115e+03 227 4.256129e+02 -4.365484e+02 228 5.171134e+02 4.256129e+02 229 4.256129e+02 5.171134e+02 230 2.231932e+02 4.256129e+02 231 1.795512e+02 2.231932e+02 232 4.256129e+02 1.795512e+02 233 -4.253277e+02 4.256129e+02 234 4.256129e+02 -4.253277e+02 235 4.270577e+02 4.256129e+02 236 4.256129e+02 4.270577e+02 237 -1.518905e+02 4.256129e+02 238 4.256129e+02 -1.518905e+02 239 4.256129e+02 4.256129e+02 240 4.256129e+02 4.256129e+02 241 2.427915e+02 4.256129e+02 242 4.256129e+02 2.427915e+02 243 -2.921195e+02 4.256129e+02 244 2.678262e+03 -2.921195e+02 245 4.256129e+02 2.678262e+03 246 4.256129e+02 4.256129e+02 247 -9.183872e+02 4.256129e+02 248 4.256129e+02 -9.183872e+02 249 4.256129e+02 4.256129e+02 250 -2.644955e+02 4.256129e+02 251 4.256129e+02 -2.644955e+02 252 6.223636e+02 4.256129e+02 253 -2.569747e+03 6.223636e+02 254 -2.206108e+03 -2.569747e+03 255 4.256129e+02 -2.206108e+03 256 1.627109e+01 4.256129e+02 257 -1.973864e+03 1.627109e+01 258 4.256129e+02 -1.973864e+03 259 4.256129e+02 4.256129e+02 260 4.256129e+02 4.256129e+02 261 4.256129e+02 4.256129e+02 262 -7.367878e+02 4.256129e+02 263 -4.816273e+02 -7.367878e+02 264 4.256129e+02 -4.816273e+02 265 4.256129e+02 4.256129e+02 266 4.256129e+02 4.256129e+02 267 -6.098594e+02 4.256129e+02 268 4.256129e+02 -6.098594e+02 269 -2.109743e+03 4.256129e+02 270 4.256129e+02 -2.109743e+03 271 4.256129e+02 4.256129e+02 272 4.256129e+02 4.256129e+02 273 -1.366331e+03 4.256129e+02 274 4.256129e+02 -1.366331e+03 275 -5.444077e+02 4.256129e+02 276 -1.602230e+03 -5.444077e+02 277 4.256129e+02 -1.602230e+03 278 -9.220581e+02 4.256129e+02 279 -5.467651e+01 -9.220581e+02 280 4.256129e+02 -5.467651e+01 281 4.256129e+02 4.256129e+02 282 -2.066203e+03 4.256129e+02 283 4.256129e+02 -2.066203e+03 284 -1.388861e+02 4.256129e+02 285 4.256129e+02 -1.388861e+02 286 4.754603e+02 4.256129e+02 287 1.441942e+02 4.754603e+02 288 4.256129e+02 1.441942e+02 289 -5.787738e+02 4.256129e+02 290 4.256129e+02 -5.787738e+02 291 -7.006926e+02 4.256129e+02 292 -1.182706e+03 -7.006926e+02 293 4.256129e+02 -1.182706e+03 294 -4.606910e+02 4.256129e+02 295 -3.748129e+02 -4.606910e+02 296 4.256129e+02 -3.748129e+02 297 4.256129e+02 4.256129e+02 298 4.256129e+02 4.256129e+02 299 -5.049457e+02 4.256129e+02 300 6.316016e+01 -5.049457e+02 301 7.055370e+02 6.316016e+01 302 3.222175e+02 7.055370e+02 303 4.256129e+02 3.222175e+02 304 -8.961924e+02 4.256129e+02 305 -8.129961e+02 -8.961924e+02 306 4.256129e+02 -8.129961e+02 307 -1.283940e+03 4.256129e+02 308 -1.003526e+03 -1.283940e+03 309 4.256129e+02 -1.003526e+03 310 4.256129e+02 4.256129e+02 311 -2.336793e+03 4.256129e+02 312 4.256129e+02 -2.336793e+03 313 -9.734246e+02 4.256129e+02 314 -4.003007e+02 -9.734246e+02 315 2.185694e+01 -4.003007e+02 316 4.486071e+02 2.185694e+01 317 -1.103780e+03 4.486071e+02 318 -9.031008e+01 -1.103780e+03 319 4.256129e+02 -9.031008e+01 320 4.256129e+02 4.256129e+02 321 4.256129e+02 4.256129e+02 322 4.256129e+02 4.256129e+02 323 4.256129e+02 4.256129e+02 324 4.256129e+02 4.256129e+02 325 3.699330e+02 4.256129e+02 326 4.256129e+02 3.699330e+02 327 4.256129e+02 4.256129e+02 328 4.256129e+02 4.256129e+02 329 4.256129e+02 4.256129e+02 330 4.256129e+02 4.256129e+02 331 4.256129e+02 4.256129e+02 332 4.256129e+02 4.256129e+02 333 1.168389e+02 4.256129e+02 334 3.232677e+02 1.168389e+02 335 3.261116e+02 3.232677e+02 336 -6.230431e+02 3.261116e+02 337 4.256129e+02 -6.230431e+02 338 3.646600e+02 4.256129e+02 339 -5.211557e+02 3.646600e+02 340 -2.953794e+02 -5.211557e+02 341 4.256129e+02 -2.953794e+02 342 4.256129e+02 4.256129e+02 343 4.256129e+02 4.256129e+02 344 -2.339587e+02 4.256129e+02 345 4.256129e+02 -2.339587e+02 346 3.110211e+02 4.256129e+02 347 4.256129e+02 3.110211e+02 348 4.256129e+02 4.256129e+02 349 9.510645e+02 4.256129e+02 350 3.610234e+02 9.510645e+02 351 -2.192591e+02 3.610234e+02 352 -1.329647e+02 -2.192591e+02 353 4.910256e+02 -1.329647e+02 354 4.256129e+02 4.910256e+02 355 4.256129e+02 4.256129e+02 356 -3.599873e+02 4.256129e+02 357 4.256129e+02 -3.599873e+02 358 4.256129e+02 4.256129e+02 359 4.256129e+02 4.256129e+02 360 1.290543e+01 4.256129e+02 361 4.256129e+02 1.290543e+01 362 4.256129e+02 4.256129e+02 363 -6.283143e+02 4.256129e+02 364 -3.527732e+02 -6.283143e+02 365 4.256129e+02 -3.527732e+02 366 4.256129e+02 4.256129e+02 367 4.256129e+02 4.256129e+02 368 4.256129e+02 4.256129e+02 369 -8.169858e+01 4.256129e+02 370 4.256129e+02 -8.169858e+01 371 6.231531e+01 4.256129e+02 372 1.581195e+01 6.231531e+01 373 4.316527e+02 1.581195e+01 374 4.256129e+02 4.316527e+02 375 4.256129e+02 4.256129e+02 376 -4.910971e+02 4.256129e+02 377 2.418058e+02 -4.910971e+02 378 1.563648e+03 2.418058e+02 379 4.256129e+02 1.563648e+03 380 4.256129e+02 4.256129e+02 381 -7.016803e+02 4.256129e+02 382 -4.731015e+02 -7.016803e+02 383 4.256129e+02 -4.731015e+02 384 -7.201314e+02 4.256129e+02 385 7.285816e+02 -7.201314e+02 386 1.237071e+02 7.285816e+02 387 7.483456e+02 1.237071e+02 388 1.922723e+02 7.483456e+02 389 1.997467e+02 1.922723e+02 390 -2.387558e+02 1.997467e+02 391 4.384133e+02 -2.387558e+02 392 6.893772e+02 4.384133e+02 393 1.603085e+02 6.893772e+02 394 3.598731e+03 1.603085e+02 395 4.048085e+02 3.598731e+03 396 3.156979e+03 4.048085e+02 397 -2.289213e+04 3.156979e+03 398 4.561749e+02 -2.289213e+04 399 -1.621360e+03 4.561749e+02 400 1.444165e+03 -1.621360e+03 401 2.737006e+02 1.444165e+03 402 5.543884e+02 2.737006e+02 403 2.549781e+02 5.543884e+02 404 6.989033e+02 2.549781e+02 405 5.398027e+02 6.989033e+02 406 2.096736e+01 5.398027e+02 407 2.393617e+02 2.096736e+01 408 4.558431e+01 2.393617e+02 409 1.747562e+03 4.558431e+01 410 -3.952910e+03 1.747562e+03 411 5.970922e+03 -3.952910e+03 412 3.867466e+02 5.970922e+03 413 -5.064982e+02 3.867466e+02 414 2.871415e+03 -5.064982e+02 415 2.396648e+03 2.871415e+03 416 4.204533e+01 2.396648e+03 417 -2.578826e+02 4.204533e+01 418 -2.808991e+02 -2.578826e+02 419 1.687823e+03 -2.808991e+02 420 2.044902e+03 1.687823e+03 421 -3.509178e+03 2.044902e+03 422 7.172237e+03 -3.509178e+03 423 7.988975e+03 7.172237e+03 424 9.060725e+02 7.988975e+03 425 -1.941565e+03 9.060725e+02 426 4.521493e+03 -1.941565e+03 427 1.257998e+04 4.521493e+03 428 7.247252e+03 1.257998e+04 429 6.549244e+03 7.247252e+03 430 1.366157e+04 6.549244e+03 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/706831292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/8tfpo1292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/9tfpo1292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10tfpo1292938441.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/116p5f1292938441.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12apml1292938441.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13h9m61292938442.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14r03r1292938442.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/15uj2f1292938442.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/16rth61292938442.tab") + } > > try(system("convert tmp/1e5sx1292938441.ps tmp/1e5sx1292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/2e5sx1292938441.ps tmp/2e5sx1292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/3e5sx1292938441.ps tmp/3e5sx1292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/47fr01292938441.ps tmp/47fr01292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/57fr01292938441.ps tmp/57fr01292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/67fr01292938441.ps tmp/67fr01292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/706831292938441.ps tmp/706831292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/8tfpo1292938441.ps tmp/8tfpo1292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/9tfpo1292938441.ps tmp/9tfpo1292938441.png",intern=TRUE)) character(0) > try(system("convert tmp/10tfpo1292938441.ps tmp/10tfpo1292938441.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 11.658 2.172 23.914