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(4321023 + ,29790 + ,444 + ,81767 + ,1 + ,1 + ,4111912 + ,87550 + ,412 + ,153198 + ,1 + ,1 + ,223193 + ,84738 + ,428 + ,-26007 + ,0 + ,1 + ,1491348 + ,54660 + ,315 + ,126942 + ,1 + ,1 + ,1629616 + ,42634 + ,168 + ,157214 + ,1 + ,0 + ,1398893 + ,40949 + ,263 + ,129352 + ,0 + ,1 + ,1926517 + ,45187 + ,267 + ,234817 + ,1 + ,0 + ,983660 + ,37704 + ,228 + ,60448 + ,1 + ,1 + ,1443586 + ,16275 + ,129 + ,47818 + ,1 + ,1 + ,1073089 + ,25830 + ,104 + ,245546 + ,0 + ,1 + ,984885 + ,12679 + ,122 + ,48020 + ,0 + ,1 + ,1405225 + ,18014 + ,393 + ,-1710 + ,1 + ,1 + ,227132 + ,43556 + ,190 + ,32648 + ,0 + ,1 + ,929118 + ,24811 + ,280 + ,95350 + ,1 + ,1 + ,1071292 + ,6575 + ,63 + ,151352 + ,0 + ,0 + ,638830 + ,7123 + ,102 + ,288170 + ,0 + ,0 + ,856956 + ,21950 + ,265 + ,114337 + ,1 + ,1 + ,992426 + ,37597 + ,234 + ,37884 + ,1 + ,1 + ,444477 + ,17821 + ,277 + ,122844 + ,0 + ,1 + ,857217 + ,12988 + ,73 + ,82340 + ,1 + ,1 + ,711969 + ,22330 + ,67 + ,79801 + ,1 + ,0 + ,702380 + ,13326 + ,103 + ,165548 + ,0 + ,0 + ,358589 + ,16189 + ,290 + ,116384 + ,0 + ,1 + ,297978 + ,7146 + ,83 + ,134028 + ,0 + ,1 + ,585715 + ,15824 + ,56 + ,63838 + ,0 + ,1 + ,657954 + ,27664 + ,236 + ,74996 + ,1 + ,1 + ,209458 + ,11920 + ,73 + ,31080 + ,0 + ,1 + ,786690 + ,8568 + ,34 + ,32168 + ,0 + ,1 + ,439798 + ,14416 + ,139 + ,49857 + ,0 + ,1 + ,688779 + ,3369 + ,26 + ,87161 + ,1 + ,1 + ,574339 + ,11819 + ,70 + ,106113 + ,1 + ,1 + ,741409 + ,6984 + ,40 + ,80570 + ,1 + ,0 + ,597793 + ,4519 + ,42 + ,102129 + ,1 + ,1 + ,644190 + ,2220 + ,12 + ,301670 + ,0 + ,1 + ,377934 + ,18562 + ,211 + ,102313 + ,0 + ,1 + ,640273 + ,10327 + ,74 + ,88577 + ,0 + ,1 + ,697458 + ,5336 + ,80 + ,112477 + ,1 + ,1 + ,550608 + ,2365 + ,83 + ,191778 + ,1 + ,1 + ,207393 + ,4069 + ,131 + ,79804 + ,0 + ,1 + ,301607 + ,8636 + ,203 + ,128294 + ,0 + ,0 + ,345783 + ,13718 + ,56 + ,96448 + ,0 + ,1 + ,501749 + ,4525 + ,89 + ,93811 + ,0 + ,0 + ,379983 + ,6869 + ,88 + ,117520 + ,0 + ,0 + ,387475 + ,4628 + ,39 + ,69159 + ,0 + ,0 + ,377305 + ,3689 + ,25 + ,101792 + ,1 + ,1 + ,370837 + ,4891 + ,49 + ,210568 + ,1 + ,1 + ,430866 + ,7489 + ,149 + ,136996 + ,1 + ,1 + ,469107 + ,4901 + ,58 + ,121920 + ,0 + ,1 + ,194493 + ,2284 + ,41 + ,76403 + ,0 + ,1 + ,530670 + ,3160 + ,90 + ,108094 + ,1 + ,1 + ,518365 + ,4150 + ,136 + ,134759 + ,1 + ,1 + ,491303 + ,7285 + ,97 + ,188873 + ,1 + ,1 + ,527021 + ,1134 + ,63 + ,146216 + ,1 + ,0 + ,233773 + ,4658 + ,114 + ,156608 + ,1 + ,1 + ,405972 + ,2384 + ,77 + ,61348 + ,0 + ,0 + ,652925 + ,3748 + ,6 + ,50350 + ,0 + ,1 + ,446211 + ,5371 + ,47 + ,87720 + ,0 + ,0 + ,341340 + ,1285 + ,51 + ,99489 + ,0 + ,1 + ,387699 + ,9327 + ,85 + ,87419 + ,1 + ,1 + ,493408 + ,5565 + ,43 + ,94355 + ,1 + ,0 + ,146494 + ,1528 + ,32 + ,60326 + ,0 + ,0 + ,414462 + ,3122 + ,25 + ,94670 + ,1 + ,1 + ,364304 + ,7561 + ,77 + ,82425 + ,1 + ,0 + ,355178 + ,2675 + ,54 + ,59017 + ,0 + ,0 + ,357760 + ,13253 + ,251 + ,90829 + ,0 + ,1 + ,261216 + ,880 + ,15 + ,80791 + ,0 + ,1 + ,397144 + ,2053 + ,44 + ,100423 + ,1 + ,1 + ,374943 + ,1424 + ,73 + ,131116 + ,0 + ,1 + ,424898 + ,4036 + ,85 + ,100269 + ,1 + ,1 + ,202055 + ,3045 + ,49 + ,27330 + ,1 + ,1 + ,378525 + ,5119 + ,38 + ,39039 + ,0 + ,0 + ,310768 + ,1431 + ,35 + ,106885 + ,0 + ,0 + ,325738 + ,554 + ,9 + ,79285 + ,0 + ,1 + ,394510 + ,1975 + ,34 + ,118881 + ,0 + ,1 + ,247060 + ,1765 + ,20 + ,77623 + ,1 + ,1 + ,368078 + ,1012 + ,29 + ,114768 + ,0 + ,0 + ,236761 + ,810 + ,11 + ,74015 + ,0 + ,0 + ,312378 + ,1280 + ,52 + ,69465 + ,0 + ,1 + ,339836 + ,666 + ,13 + ,117869 + ,1 + ,1 + ,347385 + ,1380 + ,29 + ,60982 + ,0 + ,1 + ,426280 + ,4677 + ,66 + ,90131 + ,1 + ,0 + ,352850 + ,876 + ,33 + ,138971 + ,0 + ,1 + ,301881 + ,814 + ,15 + ,39625 + ,0 + ,0 + ,377516 + ,514 + ,15 + ,102725 + ,0 + ,1 + ,357312 + ,5692 + ,68 + ,64239 + ,1 + ,1 + ,458343 + ,3642 + ,100 + ,90262 + ,0 + ,1 + ,354228 + ,540 + ,13 + ,103960 + ,0 + ,0 + ,308636 + ,2099 + ,45 + ,106611 + ,0 + ,0 + ,386212 + ,567 + ,14 + ,103345 + ,0 + ,0 + ,393343 + ,2001 + ,36 + ,95551 + ,0 + ,0 + ,378509 + ,2949 + ,40 + ,82903 + ,1 + ,1 + ,452469 + ,2253 + ,68 + ,63593 + ,0 + ,1 + ,364839 + ,6533 + ,29 + ,126910 + ,1 + ,1 + ,358649 + ,1889 + ,43 + ,37527 + ,0 + ,1 + ,376641 + ,3055 + ,30 + ,60247 + ,1 + ,0 + ,429112 + ,272 + ,9 + ,112995 + ,0 + ,1 + ,330546 + ,1414 + ,22 + ,70184 + ,1 + ,1 + ,403560 + ,2564 + ,19 + ,130140 + ,0 + ,0 + ,317892 + ,1383 + ,9 + ,73221 + ,1 + ,1 + ,307528 + ,1261 + ,31 + ,76114 + ,1 + ,1 + ,235133 + ,975 + ,19 + ,90534 + ,0 + ,1 + ,299243 + ,3366 + ,55 + ,108479 + ,0 + ,1 + ,314073 + ,576 + ,8 + ,113761 + ,0 + ,1 + ,368186 + ,1686 + ,28 + ,68696 + ,0 + ,1 + ,269661 + ,746 + ,29 + ,71561 + ,0 + ,1 + ,125390 + ,3192 + ,48 + ,59831 + ,1 + ,1 + ,510834 + ,2045 + ,16 + ,97890 + ,1 + ,1 + ,321896 + ,5702 + ,47 + ,101481 + ,0 + ,0 + ,249898 + ,1932 + ,20 + ,72954 + ,1 + ,1 + ,408881 + ,936 + ,22 + ,67939 + ,0 + ,1 + ,158492 + ,3437 + ,33 + ,48022 + ,1 + ,1 + ,292154 + ,5131 + ,44 + ,86111 + ,0 + ,0 + ,289513 + ,2397 + ,13 + ,74020 + ,1 + ,1 + ,378049 + ,1389 + ,6 + ,57530 + ,1 + ,1 + ,343466 + ,1503 + ,35 + ,56364 + ,0 + ,0 + ,332743 + ,402 + ,8 + ,84990 + ,0 + ,0 + ,442882 + ,2239 + ,17 + ,88590 + ,0 + ,0 + ,214215 + ,2234 + ,11 + ,77200 + ,1 + ,1 + ,315688 + ,837 + ,21 + ,61262 + ,0 + ,1 + ,375195 + ,10579 + ,92 + ,110309 + ,0 + ,1 + ,334280 + ,875 + ,12 + ,67000 + ,0 + ,0 + ,355864 + ,1585 + ,112 + ,93099 + ,0 + ,1 + ,480382 + ,1659 + ,25 + ,107577 + ,1 + ,1 + ,353058 + ,2647 + ,17 + ,62920 + ,1 + ,1 + ,217193 + ,3294 + ,23 + ,75832 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,314533 + ,94 + ,10 + ,60793 + ,0 + ,0 + ,318056 + ,422 + ,23 + ,57935 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,314353 + ,34 + ,7 + ,60630 + ,0 + ,0 + ,369448 + ,1558 + ,25 + ,55637 + ,0 + ,1 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,1 + ,312846 + ,43 + ,20 + ,60887 + ,0 + ,0 + ,312075 + ,645 + ,4 + ,60720 + ,0 + ,0 + ,315009 + ,316 + ,4 + ,60505 + ,0 + ,0 + ,318903 + ,115 + ,10 + ,60945 + ,0 + ,1 + ,314887 + ,5 + ,1 + ,60720 + ,0 + ,0 + ,314913 + ,897 + ,4 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325506 + ,389 + ,8 + ,58990 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,298568 + ,1002 + ,11 + ,56750 + ,0 + ,0 + ,315834 + ,36 + ,4 + ,60894 + ,0 + ,0 + ,329784 + ,460 + ,15 + ,63346 + ,0 + ,1 + ,312878 + ,309 + ,9 + ,56535 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,314987 + ,9 + ,7 + ,60835 + ,0 + ,0 + ,325249 + ,271 + ,2 + ,60720 + ,0 + ,1 + ,315877 + ,14 + ,0 + ,61016 + ,0 + ,1 + ,291650 + ,520 + ,7 + ,58650 + ,0 + ,0 + ,305959 + ,1766 + ,46 + ,60438 + ,0 + ,0 + ,315380 + ,0 + ,5 + ,60720 + ,0 + ,0 + ,297765 + ,458 + ,7 + ,58625 + ,0 + ,1 + ,315245 + ,20 + ,2 + ,60938 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315236 + ,98 + ,2 + ,61490 + ,0 + ,0 + ,336425 + ,405 + ,5 + ,60845 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,306268 + ,483 + ,7 + ,60830 + ,0 + ,0 + ,302187 + ,454 + ,24 + ,63261 + ,0 + ,1 + ,314882 + ,47 + ,1 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,382712 + ,757 + ,18 + ,45689 + ,0 + ,1 + ,341570 + ,4655 + ,55 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,312412 + ,36 + ,3 + ,61564 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,309596 + ,203 + ,9 + ,61938 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315547 + ,126 + ,8 + ,60951 + ,0 + ,0 + ,313267 + ,400 + ,113 + ,60720 + ,0 + ,0 + ,316176 + ,71 + ,0 + ,60745 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,359335 + ,972 + ,19 + ,71642 + ,0 + ,1 + ,330068 + ,531 + ,11 + ,71641 + ,1 + ,0 + ,314289 + ,2461 + ,25 + ,55792 + ,0 + ,0 + ,297413 + ,378 + ,16 + ,71873 + ,1 + ,0 + ,314806 + ,23 + ,5 + ,62555 + ,1 + ,1 + ,333210 + ,638 + ,11 + ,60370 + ,1 + ,1 + ,352108 + ,2300 + ,23 + ,64873 + ,1 + ,0 + ,313332 + ,149 + ,6 + ,62041 + ,0 + ,0 + ,291787 + ,226 + ,5 + ,65745 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,318745 + ,275 + ,7 + ,59500 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315366 + ,141 + ,7 + ,61630 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315688 + ,28 + ,3 + ,60890 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,409642 + ,4980 + ,89 + ,113521 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,269587 + ,472 + ,19 + ,80045 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,300962 + ,203 + ,12 + ,50804 + ,1 + ,0 + ,325479 + ,496 + ,12 + ,87390 + ,1 + ,1 + ,316155 + ,10 + ,5 + ,61656 + ,1 + ,1 + ,318574 + ,63 + ,2 + ,65688 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,343613 + ,1136 + ,26 + ,48522 + ,1 + ,0 + ,306948 + ,265 + ,3 + ,60720 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,330059 + ,267 + ,11 + ,57640 + ,0 + ,1 + ,288985 + ,474 + ,10 + ,61977 + ,0 + ,1 + ,304485 + ,534 + ,5 + ,62620 + ,0 + ,1 + ,315380 + ,0 + ,2 + ,60720 + ,0 + ,1 + ,315688 + ,15 + ,6 + ,60831 + ,0 + ,0 + ,317736 + ,397 + ,7 + ,60646 + ,0 + ,0 + ,315380 + ,0 + ,2 + ,60720 + ,0 + ,0 + ,322331 + ,1866 + ,28 + ,56225 + ,0 + ,0 + ,296656 + ,288 + ,3 + ,60510 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315354 + ,3 + ,1 + ,60698 + ,0 + ,1 + ,312161 + ,468 + ,20 + ,60720 + ,0 + ,0 + ,315576 + ,20 + ,1 + ,60805 + ,0 + ,0 + ,314922 + ,278 + ,22 + ,61404 + ,0 + ,0 + ,314551 + ,61 + ,9 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,312339 + ,192 + ,2 + ,65276 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,298700 + ,317 + ,7 + ,63915 + ,0 + ,1 + ,321376 + ,738 + ,9 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,303230 + ,368 + ,13 + ,61686 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315487 + ,2 + ,0 + ,60743 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315793 + ,53 + ,6 + ,60349 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,312887 + ,94 + ,3 + ,61360 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315637 + ,24 + ,7 + ,59818 + ,0 + ,1 + ,324385 + ,2332 + ,2 + ,72680 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,308989 + ,131 + ,15 + ,61808 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,0 + ,296702 + ,206 + ,9 + ,53110 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,307322 + ,167 + ,1 + ,64245 + ,0 + ,1 + ,304376 + ,622 + ,38 + ,73007 + ,0 + ,0 + ,253588 + ,2328 + ,57 + ,82732 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,309560 + ,365 + ,7 + ,54820 + ,0 + ,0 + ,298466 + ,364 + ,26 + ,47705 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,343929 + ,226 + ,13 + ,72835 + ,0 + ,0 + ,331955 + ,307 + ,10 + ,58856 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,0 + ,381180 + ,188 + ,9 + ,77655 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,331420 + ,138 + ,26 + ,69817 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,310201 + ,125 + ,19 + ,60798 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,320016 + ,282 + ,12 + ,62452 + ,0 + ,1 + ,320398 + ,335 + ,23 + ,64175 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,291841 + ,1324 + ,29 + ,67440 + ,1 + ,1 + ,310670 + ,176 + ,8 + ,68136 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,313491 + ,249 + ,26 + ,56726 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,331323 + ,333 + ,9 + ,70811 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,319210 + ,601 + ,5 + ,60720 + ,1 + ,1 + ,318098 + ,30 + ,3 + ,62045 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,0 + ,292754 + ,249 + ,13 + ,54323 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,325176 + ,165 + ,12 + ,62841 + ,0 + ,0 + ,365959 + ,453 + ,19 + ,81125 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,302409 + ,53 + ,10 + ,59506 + ,0 + ,1 + ,340968 + ,382 + ,9 + ,59365 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,9 + ,60720 + ,0 + ,0 + ,313164 + ,30 + ,4 + ,60798 + ,1 + ,1 + ,301164 + ,290 + ,1 + ,58790 + ,0 + ,1 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,344425 + ,366 + ,14 + ,61808 + ,0 + ,0 + ,315394 + ,2 + ,12 + ,60735 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,316647 + ,209 + ,19 + ,64016 + ,1 + ,1 + ,309836 + ,384 + ,17 + ,54683 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,346611 + ,365 + ,32 + ,87192 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,322031 + ,49 + ,14 + ,64107 + ,1 + ,0 + ,315656 + ,3 + ,8 + ,60761 + ,0 + ,0 + ,339445 + ,133 + ,4 + ,65990 + ,0 + ,0 + ,314964 + ,32 + ,0 + ,59988 + ,0 + ,0 + ,297141 + ,368 + ,20 + ,61167 + ,0 + ,0 + ,315372 + ,1 + ,5 + ,60719 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,312502 + ,22 + ,1 + ,60722 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,313729 + ,96 + ,4 + ,60379 + ,0 + ,1 + ,315388 + ,1 + ,1 + ,60727 + ,0 + ,0 + ,315371 + ,314 + ,4 + ,60720 + ,0 + ,0 + ,296139 + ,844 + ,20 + ,60925 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,313880 + ,26 + ,1 + ,60896 + ,0 + ,0 + ,317698 + ,125 + ,10 + ,59734 + ,0 + ,1 + ,295580 + ,304 + ,12 + ,62969 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,308256 + ,621 + ,13 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,303677 + ,119 + ,3 + ,59118 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,319369 + ,1595 + ,10 + ,60720 + ,0 + ,0 + ,318690 + ,312 + ,3 + ,58598 + ,0 + ,1 + ,314049 + ,60 + ,7 + ,61124 + ,0 + ,0 + ,325699 + ,587 + ,10 + ,59595 + ,0 + ,1 + ,314210 + ,135 + ,1 + ,62065 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,322378 + ,514 + ,15 + ,78780 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315398 + ,1 + ,4 + ,60722 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,308336 + ,1763 + ,28 + ,61600 + ,1 + ,1 + ,316386 + ,180 + ,9 + ,59635 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315553 + ,218 + ,7 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,323361 + ,448 + ,7 + ,59781 + ,0 + ,1 + ,336639 + ,227 + ,7 + ,76644 + ,0 + ,0 + ,307424 + ,174 + ,3 + ,64820 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,295370 + ,121 + ,11 + ,56178 + ,0 + ,1 + ,322340 + ,607 + ,7 + ,60436 + ,0 + ,0 + ,319864 + ,2212 + ,10 + ,60720 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,317291 + ,530 + ,18 + ,73433 + ,0 + ,0 + ,280398 + ,571 + ,14 + ,41477 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,1 + ,317330 + ,78 + ,12 + ,62700 + ,0 + ,0 + ,238125 + ,2489 + ,29 + ,67804 + ,0 + ,1 + ,327071 + ,131 + ,3 + ,59661 + ,0 + ,0 + ,309038 + ,923 + ,6 + ,58620 + ,0 + ,1 + ,314210 + ,72 + ,3 + ,60398 + ,0 + ,0 + ,307930 + ,572 + ,8 + ,58580 + ,0 + ,1 + ,322327 + ,397 + ,10 + ,62710 + ,0 + ,1 + ,292136 + ,450 + ,6 + ,59325 + ,0 + ,0 + ,263276 + ,622 + ,8 + ,60950 + ,0 + ,1 + ,367655 + ,694 + ,6 + ,68060 + ,0 + ,0 + ,283910 + ,3425 + ,9 + ,83620 + ,1 + ,0 + ,283587 + ,562 + ,8 + ,58456 + ,0 + ,1 + ,243650 + ,4917 + ,26 + ,52811 + ,0 + ,1 + ,438493 + ,1442 + ,239 + ,121173 + ,1 + ,1 + ,296261 + ,529 + ,7 + ,63870 + ,0 + ,1 + ,230621 + ,2126 + ,41 + ,21001 + ,1 + ,1 + ,304252 + ,1061 + ,3 + ,70415 + ,0 + ,0 + ,333505 + ,776 + ,8 + ,64230 + ,0 + ,0 + ,296919 + ,611 + ,6 + ,59190 + ,0 + ,0 + ,278990 + ,1526 + ,21 + ,69351 + ,1 + ,0 + ,276898 + ,592 + ,7 + ,64270 + ,0 + ,0 + ,327007 + ,1182 + ,11 + ,70694 + ,0 + ,0 + ,317046 + ,621 + ,11 + ,68005 + ,0 + ,1 + ,304555 + ,989 + ,12 + ,58930 + ,0 + ,1 + ,298096 + ,438 + ,9 + ,58320 + ,0 + ,0 + ,231861 + ,726 + ,3 + ,69980 + ,0 + ,1 + ,309422 + ,1303 + ,57 + ,69863 + ,0 + ,1 + ,286963 + ,7419 + ,21 + ,63255 + ,1 + ,1 + ,269753 + ,1164 + ,15 + ,57320 + ,1 + ,1 + ,448243 + ,3310 + ,32 + ,75230 + ,1 + ,1 + ,165404 + ,1920 + ,11 + ,79420 + ,0 + ,1 + ,204325 + ,965 + ,2 + ,73490 + ,0 + ,0 + ,407159 + ,3256 + ,23 + ,35250 + ,0 + ,1 + ,290476 + ,1135 + ,20 + ,62285 + ,1 + ,0 + ,275311 + ,1270 + ,24 + ,69206 + ,0 + ,0 + ,246541 + ,661 + ,1 + ,65920 + ,0 + ,0 + ,253468 + ,1013 + ,1 + ,69770 + ,0 + ,0 + ,240897 + ,2844 + ,74 + ,72683 + ,0 + ,1 + ,-83265 + ,11528 + ,68 + ,-14545 + ,1 + ,1 + ,-42143 + ,6526 + ,20 + ,55830 + ,0 + ,1 + ,272713 + ,2264 + ,20 + ,55174 + ,0 + ,1 + ,215362 + ,5109 + ,82 + ,67038 + ,1 + ,1 + ,42754 + ,3999 + ,21 + ,51252 + ,0 + ,1 + ,306275 + ,35624 + ,244 + ,157278 + ,0 + ,1 + ,253537 + ,9252 + ,32 + ,79510 + ,0 + ,1 + ,372631 + ,15236 + ,86 + ,77440 + ,0 + ,1 + ,-7170 + ,18073 + ,69 + ,27284 + ,0 + ,1) + ,dim=c(6 + ,430) + ,dimnames=list(c('Wealth' + ,'Costs' + ,'Orders' + ,'Dividends' + ,'Group' + ,'Gender ') + ,1:430)) > y <- array(NA,dim=c(6,430),dimnames=list(c('Wealth','Costs','Orders','Dividends','Group','Gender '),1:430)) > 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 Wealth Costs Orders Dividends Group Gender\r 1 4321023 29790 444 81767 1 1 2 4111912 87550 412 153198 1 1 3 223193 84738 428 -26007 0 1 4 1491348 54660 315 126942 1 1 5 1629616 42634 168 157214 1 0 6 1398893 40949 263 129352 0 1 7 1926517 45187 267 234817 1 0 8 983660 37704 228 60448 1 1 9 1443586 16275 129 47818 1 1 10 1073089 25830 104 245546 0 1 11 984885 12679 122 48020 0 1 12 1405225 18014 393 -1710 1 1 13 227132 43556 190 32648 0 1 14 929118 24811 280 95350 1 1 15 1071292 6575 63 151352 0 0 16 638830 7123 102 288170 0 0 17 856956 21950 265 114337 1 1 18 992426 37597 234 37884 1 1 19 444477 17821 277 122844 0 1 20 857217 12988 73 82340 1 1 21 711969 22330 67 79801 1 0 22 702380 13326 103 165548 0 0 23 358589 16189 290 116384 0 1 24 297978 7146 83 134028 0 1 25 585715 15824 56 63838 0 1 26 657954 27664 236 74996 1 1 27 209458 11920 73 31080 0 1 28 786690 8568 34 32168 0 1 29 439798 14416 139 49857 0 1 30 688779 3369 26 87161 1 1 31 574339 11819 70 106113 1 1 32 741409 6984 40 80570 1 0 33 597793 4519 42 102129 1 1 34 644190 2220 12 301670 0 1 35 377934 18562 211 102313 0 1 36 640273 10327 74 88577 0 1 37 697458 5336 80 112477 1 1 38 550608 2365 83 191778 1 1 39 207393 4069 131 79804 0 1 40 301607 8636 203 128294 0 0 41 345783 13718 56 96448 0 1 42 501749 4525 89 93811 0 0 43 379983 6869 88 117520 0 0 44 387475 4628 39 69159 0 0 45 377305 3689 25 101792 1 1 46 370837 4891 49 210568 1 1 47 430866 7489 149 136996 1 1 48 469107 4901 58 121920 0 1 49 194493 2284 41 76403 0 1 50 530670 3160 90 108094 1 1 51 518365 4150 136 134759 1 1 52 491303 7285 97 188873 1 1 53 527021 1134 63 146216 1 0 54 233773 4658 114 156608 1 1 55 405972 2384 77 61348 0 0 56 652925 3748 6 50350 0 1 57 446211 5371 47 87720 0 0 58 341340 1285 51 99489 0 1 59 387699 9327 85 87419 1 1 60 493408 5565 43 94355 1 0 61 146494 1528 32 60326 0 0 62 414462 3122 25 94670 1 1 63 364304 7561 77 82425 1 0 64 355178 2675 54 59017 0 0 65 357760 13253 251 90829 0 1 66 261216 880 15 80791 0 1 67 397144 2053 44 100423 1 1 68 374943 1424 73 131116 0 1 69 424898 4036 85 100269 1 1 70 202055 3045 49 27330 1 1 71 378525 5119 38 39039 0 0 72 310768 1431 35 106885 0 0 73 325738 554 9 79285 0 1 74 394510 1975 34 118881 0 1 75 247060 1765 20 77623 1 1 76 368078 1012 29 114768 0 0 77 236761 810 11 74015 0 0 78 312378 1280 52 69465 0 1 79 339836 666 13 117869 1 1 80 347385 1380 29 60982 0 1 81 426280 4677 66 90131 1 0 82 352850 876 33 138971 0 1 83 301881 814 15 39625 0 0 84 377516 514 15 102725 0 1 85 357312 5692 68 64239 1 1 86 458343 3642 100 90262 0 1 87 354228 540 13 103960 0 0 88 308636 2099 45 106611 0 0 89 386212 567 14 103345 0 0 90 393343 2001 36 95551 0 0 91 378509 2949 40 82903 1 1 92 452469 2253 68 63593 0 1 93 364839 6533 29 126910 1 1 94 358649 1889 43 37527 0 1 95 376641 3055 30 60247 1 0 96 429112 272 9 112995 0 1 97 330546 1414 22 70184 1 1 98 403560 2564 19 130140 0 0 99 317892 1383 9 73221 1 1 100 307528 1261 31 76114 1 1 101 235133 975 19 90534 0 1 102 299243 3366 55 108479 0 1 103 314073 576 8 113761 0 1 104 368186 1686 28 68696 0 1 105 269661 746 29 71561 0 1 106 125390 3192 48 59831 1 1 107 510834 2045 16 97890 1 1 108 321896 5702 47 101481 0 0 109 249898 1932 20 72954 1 1 110 408881 936 22 67939 0 1 111 158492 3437 33 48022 1 1 112 292154 5131 44 86111 0 0 113 289513 2397 13 74020 1 1 114 378049 1389 6 57530 1 1 115 343466 1503 35 56364 0 0 116 332743 402 8 84990 0 0 117 442882 2239 17 88590 0 0 118 214215 2234 11 77200 1 1 119 315688 837 21 61262 0 1 120 375195 10579 92 110309 0 1 121 334280 875 12 67000 0 0 122 355864 1585 112 93099 0 1 123 480382 1659 25 107577 1 1 124 353058 2647 17 62920 1 1 125 217193 3294 23 75832 1 0 126 315380 0 0 60720 0 0 127 314533 94 10 60793 0 0 128 318056 422 23 57935 0 0 129 315380 0 0 60720 0 1 130 314353 34 7 60630 0 0 131 369448 1558 25 55637 0 1 132 315380 0 1 60720 0 1 133 312846 43 20 60887 0 0 134 312075 645 4 60720 0 0 135 315009 316 4 60505 0 0 136 318903 115 10 60945 0 1 137 314887 5 1 60720 0 0 138 314913 897 4 60720 0 0 139 315380 0 0 60720 0 0 140 325506 389 8 58990 0 0 141 315380 0 0 60720 0 1 142 298568 1002 11 56750 0 0 143 315834 36 4 60894 0 0 144 329784 460 15 63346 0 1 145 312878 309 9 56535 0 0 146 315380 0 0 60720 0 1 147 314987 9 7 60835 0 0 148 325249 271 2 60720 0 1 149 315877 14 0 61016 0 1 150 291650 520 7 58650 0 0 151 305959 1766 46 60438 0 0 152 315380 0 5 60720 0 0 153 297765 458 7 58625 0 1 154 315245 20 2 60938 0 0 155 315380 0 0 60720 0 0 156 315380 0 0 60720 0 0 157 315236 98 2 61490 0 0 158 336425 405 5 60845 0 0 159 315380 0 0 60720 0 0 160 315380 0 0 60720 0 0 161 315380 0 0 60720 0 1 162 315380 0 0 60720 0 0 163 306268 483 7 60830 0 0 164 302187 454 24 63261 0 1 165 314882 47 1 60720 0 0 166 315380 0 0 60720 0 1 167 382712 757 18 45689 0 1 168 341570 4655 55 60720 0 1 169 315380 0 0 60720 0 0 170 315380 0 0 60720 0 0 171 312412 36 3 61564 0 1 172 315380 0 0 60720 0 1 173 309596 203 9 61938 0 1 174 315380 0 0 60720 0 1 175 315547 126 8 60951 0 0 176 313267 400 113 60720 0 0 177 316176 71 0 60745 1 1 178 315380 0 0 60720 0 0 179 315380 0 0 60720 0 0 180 359335 972 19 71642 0 1 181 330068 531 11 71641 1 0 182 314289 2461 25 55792 0 0 183 297413 378 16 71873 1 0 184 314806 23 5 62555 1 1 185 333210 638 11 60370 1 1 186 352108 2300 23 64873 1 0 187 313332 149 6 62041 0 0 188 291787 226 5 65745 0 1 189 315380 0 0 60720 0 0 190 318745 275 7 59500 0 0 191 315380 0 0 60720 0 0 192 315366 141 7 61630 0 0 193 315380 0 0 60720 0 1 194 315688 28 3 60890 0 0 195 315380 0 0 60720 1 1 196 409642 4980 89 113521 1 1 197 315380 0 0 60720 1 1 198 315380 0 0 60720 1 1 199 269587 472 19 80045 1 1 200 315380 0 0 60720 1 1 201 315380 0 0 60720 1 1 202 315380 0 0 60720 1 1 203 300962 203 12 50804 1 0 204 325479 496 12 87390 1 1 205 316155 10 5 61656 1 1 206 318574 63 2 65688 1 1 207 315380 0 0 60720 1 1 208 343613 1136 26 48522 1 0 209 306948 265 3 60720 1 0 210 315380 0 0 60720 0 0 211 315380 0 0 60720 0 1 212 330059 267 11 57640 0 1 213 288985 474 10 61977 0 1 214 304485 534 5 62620 0 1 215 315380 0 2 60720 0 1 216 315688 15 6 60831 0 0 217 317736 397 7 60646 0 0 218 315380 0 2 60720 0 0 219 322331 1866 28 56225 0 0 220 296656 288 3 60510 0 1 221 315380 0 0 60720 0 1 222 315354 3 1 60698 0 1 223 312161 468 20 60720 0 0 224 315576 20 1 60805 0 0 225 314922 278 22 61404 0 0 226 314551 61 9 60720 0 0 227 315380 0 0 60720 0 1 228 312339 192 2 65276 0 1 229 315380 0 0 60720 0 0 230 298700 317 7 63915 0 1 231 321376 738 9 60720 0 1 232 315380 0 0 60720 0 1 233 303230 368 13 61686 0 0 234 315380 0 0 60720 0 0 235 315487 2 0 60743 0 0 236 315380 0 0 60720 0 0 237 315793 53 6 60349 0 0 238 315380 0 0 60720 0 1 239 315380 0 0 60720 0 1 240 315380 0 0 60720 0 0 241 312887 94 3 61360 0 1 242 315380 0 0 60720 0 0 243 315637 24 7 59818 0 1 244 324385 2332 2 72680 0 1 245 315380 0 0 60720 1 1 246 315380 0 0 60720 1 1 247 308989 131 15 61808 0 0 248 315380 0 0 60720 1 1 249 315380 0 0 60720 1 0 250 296702 206 9 53110 0 1 251 315380 0 0 60720 1 1 252 307322 167 1 64245 0 1 253 304376 622 38 73007 0 0 254 253588 2328 57 82732 0 0 255 315380 0 0 60720 0 0 256 309560 365 7 54820 0 0 257 298466 364 26 47705 0 0 258 315380 0 0 60720 1 1 259 315380 0 0 60720 0 1 260 315380 0 0 60720 1 1 261 315380 0 0 60720 1 1 262 343929 226 13 72835 0 0 263 331955 307 10 58856 0 1 264 315380 0 0 60720 1 1 265 315380 0 0 60720 0 1 266 315380 0 0 60720 1 0 267 381180 188 9 77655 0 1 268 315380 0 0 60720 1 1 269 331420 138 26 69817 0 1 270 315380 0 0 60720 1 1 271 315380 0 0 60720 1 1 272 315380 0 0 60720 1 1 273 310201 125 19 60798 0 0 274 315380 0 0 60720 0 0 275 320016 282 12 62452 0 1 276 320398 335 23 64175 0 0 277 315380 0 0 60720 0 1 278 291841 1324 29 67440 1 1 279 310670 176 8 68136 0 1 280 315380 0 0 60720 1 1 281 315380 0 0 60720 0 0 282 313491 249 26 56726 0 1 283 315380 0 0 60720 0 0 284 331323 333 9 70811 0 1 285 315380 0 0 60720 0 1 286 319210 601 5 60720 1 1 287 318098 30 3 62045 0 0 288 315380 0 0 60720 1 0 289 292754 249 13 54323 0 0 290 315380 0 0 60720 1 1 291 325176 165 12 62841 0 0 292 365959 453 19 81125 0 1 293 315380 0 0 60720 0 0 294 302409 53 10 59506 0 1 295 340968 382 9 59365 1 1 296 315380 0 0 60720 0 0 297 315380 0 0 60720 1 0 298 315380 0 0 60720 0 1 299 315380 0 9 60720 0 0 300 313164 30 4 60798 1 1 301 301164 290 1 58790 0 1 302 315380 0 1 60720 0 1 303 315380 0 0 60720 1 1 304 344425 366 14 61808 0 0 305 315394 2 12 60735 0 0 306 315380 0 0 60720 1 1 307 316647 209 19 64016 1 1 308 309836 384 17 54683 0 0 309 315380 0 0 60720 1 1 310 315380 0 0 60720 1 1 311 346611 365 32 87192 0 1 312 315380 0 0 60720 1 1 313 322031 49 14 64107 1 0 314 315656 3 8 60761 0 0 315 339445 133 4 65990 0 0 316 314964 32 0 59988 0 0 317 297141 368 20 61167 0 0 318 315372 1 5 60719 0 0 319 315380 0 0 60720 0 0 320 315380 0 0 60720 0 1 321 315380 0 0 60720 0 1 322 315380 0 0 60720 0 1 323 315380 0 0 60720 1 0 324 315380 0 0 60720 0 1 325 312502 22 1 60722 0 1 326 315380 0 0 60720 1 0 327 315380 0 0 60720 1 1 328 315380 0 0 60720 0 1 329 315380 0 0 60720 0 1 330 315380 0 0 60720 0 1 331 315380 0 0 60720 0 0 332 315380 0 0 60720 0 1 333 313729 96 4 60379 0 1 334 315388 1 1 60727 0 0 335 315371 314 4 60720 0 0 336 296139 844 20 60925 0 1 337 315380 0 0 60720 0 0 338 313880 26 1 60896 0 0 339 317698 125 10 59734 0 1 340 295580 304 12 62969 0 0 341 315380 0 0 60720 0 0 342 315380 0 0 60720 0 0 343 315380 0 0 60720 0 1 344 308256 621 13 60720 0 0 345 315380 0 0 60720 0 0 346 303677 119 3 59118 0 0 347 315380 0 0 60720 0 0 348 315380 0 0 60720 0 1 349 319369 1595 10 60720 0 0 350 318690 312 3 58598 0 1 351 314049 60 7 61124 0 0 352 325699 587 10 59595 0 1 353 314210 135 1 62065 0 0 354 315380 0 0 60720 0 0 355 315380 0 0 60720 0 0 356 322378 514 15 78780 0 0 357 315380 0 0 60720 0 0 358 315380 0 0 60720 0 0 359 315380 0 0 60720 0 1 360 315398 1 4 60722 0 0 361 315380 0 0 60720 0 1 362 315380 0 0 60720 0 0 363 308336 1763 28 61600 1 1 364 316386 180 9 59635 0 1 365 315380 0 0 60720 0 1 366 315380 0 0 60720 0 1 367 315380 0 0 60720 0 0 368 315380 0 0 60720 0 1 369 315553 218 7 60720 0 0 370 315380 0 0 60720 0 1 371 323361 448 7 59781 0 1 372 336639 227 7 76644 0 0 373 307424 174 3 64820 0 1 374 315380 0 0 60720 0 1 375 315380 0 0 60720 0 1 376 295370 121 11 56178 0 1 377 322340 607 7 60436 0 0 378 319864 2212 10 60720 0 0 379 315380 0 0 60720 0 1 380 315380 0 0 60720 0 1 381 317291 530 18 73433 0 0 382 280398 571 14 41477 0 0 383 315380 0 0 60720 0 1 384 317330 78 12 62700 0 0 385 238125 2489 29 67804 0 1 386 327071 131 3 59661 0 0 387 309038 923 6 58620 0 1 388 314210 72 3 60398 0 0 389 307930 572 8 58580 0 1 390 322327 397 10 62710 0 1 391 292136 450 6 59325 0 0 392 263276 622 8 60950 0 1 393 367655 694 6 68060 0 0 394 283910 3425 9 83620 1 0 395 283587 562 8 58456 0 1 396 243650 4917 26 52811 0 1 397 438493 1442 239 121173 1 1 398 296261 529 7 63870 0 1 399 230621 2126 41 21001 1 1 400 304252 1061 3 70415 0 0 401 333505 776 8 64230 0 0 402 296919 611 6 59190 0 0 403 278990 1526 21 69351 1 0 404 276898 592 7 64270 0 0 405 327007 1182 11 70694 0 0 406 317046 621 11 68005 0 1 407 304555 989 12 58930 0 1 408 298096 438 9 58320 0 0 409 231861 726 3 69980 0 1 410 309422 1303 57 69863 0 1 411 286963 7419 21 63255 1 1 412 269753 1164 15 57320 1 1 413 448243 3310 32 75230 1 1 414 165404 1920 11 79420 0 1 415 204325 965 2 73490 0 0 416 407159 3256 23 35250 0 1 417 290476 1135 20 62285 1 0 418 275311 1270 24 69206 0 0 419 246541 661 1 65920 0 0 420 253468 1013 1 69770 0 0 421 240897 2844 74 72683 0 1 422 -83265 11528 68 -14545 1 1 423 -42143 6526 20 55830 0 1 424 272713 2264 20 55174 0 1 425 215362 5109 82 67038 1 1 426 42754 3999 21 51252 0 1 427 306275 35624 244 157278 0 1 428 253537 9252 32 79510 0 1 429 372631 15236 86 77440 0 1 430 -7170 18073 69 27284 0 1 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Costs Orders Dividends Group `Gender\r` 197015.803 12.618 1501.735 1.356 66885.361 -42165.685 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -1608369 -36822 20648 51704 2945709 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.970e+05 3.124e+04 6.307 7.14e-10 *** Costs 1.262e+01 2.322e+00 5.434 9.32e-08 *** Orders 1.502e+03 3.337e+02 4.500 8.78e-06 *** Dividends 1.356e+00 3.839e-01 3.533 0.000456 *** Group 6.689e+04 2.674e+04 2.501 0.012745 * `Gender\r` -4.217e+04 2.403e+04 -1.755 0.080049 . --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 235600 on 424 degrees of freedom Multiple R-squared: 0.4759, Adjusted R-squared: 0.4697 F-statistic: 77.01 on 5 and 424 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 3.119064e-53 1.559532e-53 [2,] 1.0000000 9.391202e-54 4.695601e-54 [3,] 1.0000000 2.804233e-65 1.402117e-65 [4,] 1.0000000 2.178794e-89 1.089397e-89 [5,] 1.0000000 4.181877e-91 2.090939e-91 [6,] 1.0000000 2.087251e-106 1.043626e-106 [7,] 1.0000000 2.655190e-117 1.327595e-117 [8,] 1.0000000 5.038473e-122 2.519236e-122 [9,] 1.0000000 2.112852e-132 1.056426e-132 [10,] 1.0000000 2.896492e-137 1.448246e-137 [11,] 1.0000000 1.732373e-140 8.661863e-141 [12,] 1.0000000 2.895087e-145 1.447543e-145 [13,] 1.0000000 1.463764e-146 7.318819e-147 [14,] 1.0000000 3.203861e-147 1.601930e-147 [15,] 1.0000000 1.626336e-149 8.131679e-150 [16,] 1.0000000 5.983008e-149 2.991504e-149 [17,] 1.0000000 1.149654e-152 5.748272e-153 [18,] 1.0000000 2.928685e-158 1.464343e-158 [19,] 1.0000000 6.008239e-158 3.004119e-158 [20,] 1.0000000 1.739595e-172 8.697973e-173 [21,] 1.0000000 1.161595e-172 5.807976e-173 [22,] 1.0000000 1.261577e-176 6.307883e-177 [23,] 1.0000000 1.086674e-178 5.433368e-179 [24,] 1.0000000 3.225017e-186 1.612509e-186 [25,] 1.0000000 2.193136e-188 1.096568e-188 [26,] 1.0000000 2.002572e-187 1.001286e-187 [27,] 1.0000000 3.010327e-187 1.505164e-187 [28,] 1.0000000 1.676785e-193 8.383923e-194 [29,] 1.0000000 2.675121e-199 1.337560e-199 [30,] 1.0000000 4.075137e-200 2.037568e-200 [31,] 1.0000000 2.007825e-200 1.003913e-200 [32,] 1.0000000 2.315420e-201 1.157710e-201 [33,] 1.0000000 1.145022e-200 5.725110e-201 [34,] 1.0000000 3.235417e-201 1.617709e-201 [35,] 1.0000000 3.952575e-200 1.976288e-200 [36,] 1.0000000 5.801512e-200 2.900756e-200 [37,] 1.0000000 2.549246e-199 1.274623e-199 [38,] 1.0000000 1.921825e-201 9.609124e-202 [39,] 1.0000000 3.638143e-202 1.819071e-202 [40,] 1.0000000 6.835358e-202 3.417679e-202 [41,] 1.0000000 3.434304e-202 1.717152e-202 [42,] 1.0000000 8.509666e-203 4.254833e-203 [43,] 1.0000000 4.624594e-203 2.312297e-203 [44,] 1.0000000 5.660809e-203 2.830405e-203 [45,] 1.0000000 9.425505e-203 4.712752e-203 [46,] 1.0000000 5.852011e-207 2.926006e-207 [47,] 1.0000000 1.494129e-206 7.470643e-207 [48,] 1.0000000 1.149188e-220 5.745941e-221 [49,] 1.0000000 5.270960e-221 2.635480e-221 [50,] 1.0000000 3.863063e-220 1.931532e-220 [51,] 1.0000000 4.912123e-220 2.456062e-220 [52,] 1.0000000 1.562514e-221 7.812569e-222 [53,] 1.0000000 4.198600e-223 2.099300e-223 [54,] 1.0000000 1.543670e-222 7.718348e-223 [55,] 1.0000000 2.851558e-222 1.425779e-222 [56,] 1.0000000 1.712246e-221 8.561229e-222 [57,] 1.0000000 3.167759e-221 1.583880e-221 [58,] 1.0000000 6.530479e-221 3.265240e-221 [59,] 1.0000000 6.028345e-220 3.014173e-220 [60,] 1.0000000 3.164635e-219 1.582317e-219 [61,] 1.0000000 1.180140e-218 5.900698e-219 [62,] 1.0000000 7.196561e-218 3.598280e-218 [63,] 1.0000000 6.924825e-219 3.462413e-219 [64,] 1.0000000 2.197226e-218 1.098613e-218 [65,] 1.0000000 1.189270e-217 5.946349e-218 [66,] 1.0000000 1.008801e-216 5.044007e-217 [67,] 1.0000000 2.120505e-216 1.060252e-216 [68,] 1.0000000 1.795763e-215 8.978817e-216 [69,] 1.0000000 2.357403e-215 1.178702e-215 [70,] 1.0000000 1.985358e-214 9.926792e-215 [71,] 1.0000000 7.346452e-214 3.673226e-214 [72,] 1.0000000 3.003100e-213 1.501550e-213 [73,] 1.0000000 3.503251e-213 1.751626e-213 [74,] 1.0000000 3.078370e-213 1.539185e-213 [75,] 1.0000000 1.976286e-212 9.881428e-213 [76,] 1.0000000 1.354311e-211 6.771557e-212 [77,] 1.0000000 2.780740e-211 1.390370e-211 [78,] 1.0000000 1.956036e-211 9.780180e-212 [79,] 1.0000000 1.403677e-210 7.018385e-211 [80,] 1.0000000 2.735503e-210 1.367752e-210 [81,] 1.0000000 2.592666e-209 1.296333e-209 [82,] 1.0000000 2.335285e-208 1.167643e-208 [83,] 1.0000000 1.449452e-207 7.247261e-208 [84,] 1.0000000 3.665838e-209 1.832919e-209 [85,] 1.0000000 2.861349e-208 1.430674e-208 [86,] 1.0000000 7.671012e-209 3.835506e-209 [87,] 1.0000000 8.183865e-209 4.091932e-209 [88,] 1.0000000 5.041356e-208 2.520678e-208 [89,] 1.0000000 5.299331e-207 2.649666e-207 [90,] 1.0000000 4.175145e-206 2.087573e-206 [91,] 1.0000000 4.538039e-205 2.269019e-205 [92,] 1.0000000 3.976922e-204 1.988461e-204 [93,] 1.0000000 3.148935e-205 1.574467e-205 [94,] 1.0000000 3.000841e-205 1.500421e-205 [95,] 1.0000000 4.306849e-206 2.153424e-206 [96,] 1.0000000 1.524646e-205 7.623232e-206 [97,] 1.0000000 4.739711e-205 2.369855e-205 [98,] 1.0000000 2.011683e-207 1.005841e-207 [99,] 1.0000000 1.005795e-208 5.028975e-209 [100,] 1.0000000 7.115215e-208 3.557608e-208 [101,] 1.0000000 1.389550e-207 6.947748e-208 [102,] 1.0000000 7.618880e-208 3.809440e-208 [103,] 1.0000000 3.668680e-208 1.834340e-208 [104,] 1.0000000 2.631100e-207 1.315550e-207 [105,] 1.0000000 2.076715e-206 1.038358e-206 [106,] 1.0000000 1.737609e-206 8.688046e-207 [107,] 1.0000000 7.196635e-206 3.598318e-206 [108,] 1.0000000 4.943050e-205 2.471525e-205 [109,] 1.0000000 4.028912e-205 2.014456e-205 [110,] 1.0000000 3.825979e-206 1.912989e-206 [111,] 1.0000000 3.229987e-205 1.614994e-205 [112,] 1.0000000 1.984586e-204 9.922928e-205 [113,] 1.0000000 1.832957e-203 9.164785e-204 [114,] 1.0000000 1.575958e-202 7.879789e-203 [115,] 1.0000000 2.464870e-202 1.232435e-202 [116,] 1.0000000 7.421702e-202 3.710851e-202 [117,] 1.0000000 1.415200e-202 7.076000e-203 [118,] 1.0000000 1.322493e-201 6.612466e-202 [119,] 1.0000000 1.308347e-200 6.541737e-201 [120,] 1.0000000 1.258909e-199 6.294546e-200 [121,] 1.0000000 9.268900e-199 4.634450e-199 [122,] 1.0000000 9.083554e-198 4.541777e-198 [123,] 1.0000000 5.931070e-198 2.965535e-198 [124,] 1.0000000 4.582489e-197 2.291245e-197 [125,] 1.0000000 4.730272e-196 2.365136e-196 [126,] 1.0000000 4.648341e-195 2.324170e-195 [127,] 1.0000000 4.509954e-194 2.254977e-194 [128,] 1.0000000 3.629916e-193 1.814958e-193 [129,] 1.0000000 3.501883e-192 1.750941e-192 [130,] 1.0000000 3.325466e-191 1.662733e-191 [131,] 1.0000000 3.188987e-190 1.594494e-190 [132,] 1.0000000 2.673568e-189 1.336784e-189 [133,] 1.0000000 2.085916e-188 1.042958e-188 [134,] 1.0000000 2.094039e-187 1.047019e-187 [135,] 1.0000000 2.019585e-186 1.009792e-186 [136,] 1.0000000 1.540432e-185 7.702161e-186 [137,] 1.0000000 1.415270e-184 7.076352e-185 [138,] 1.0000000 1.114503e-183 5.572515e-184 [139,] 1.0000000 1.075315e-182 5.376575e-183 [140,] 1.0000000 7.617526e-182 3.808763e-182 [141,] 1.0000000 6.075973e-181 3.037986e-181 [142,] 1.0000000 5.731942e-180 2.865971e-180 [143,] 1.0000000 5.553901e-179 2.776951e-179 [144,] 1.0000000 5.228394e-178 2.614197e-178 [145,] 1.0000000 4.500132e-177 2.250066e-177 [146,] 1.0000000 4.186493e-176 2.093246e-176 [147,] 1.0000000 3.857105e-175 1.928553e-175 [148,] 1.0000000 3.543484e-174 1.771742e-174 [149,] 1.0000000 3.276616e-173 1.638308e-173 [150,] 1.0000000 2.293475e-172 1.146737e-172 [151,] 1.0000000 2.088166e-171 1.044083e-171 [152,] 1.0000000 1.894897e-170 9.474486e-171 [153,] 1.0000000 1.476460e-169 7.382300e-170 [154,] 1.0000000 1.330815e-168 6.654075e-169 [155,] 1.0000000 1.229603e-167 6.148016e-168 [156,] 1.0000000 1.045383e-166 5.226914e-167 [157,] 1.0000000 9.340998e-166 4.670499e-166 [158,] 1.0000000 7.230059e-165 3.615029e-165 [159,] 1.0000000 6.864815e-166 3.432408e-166 [160,] 1.0000000 1.050247e-165 5.251235e-166 [161,] 1.0000000 9.443733e-165 4.721866e-165 [162,] 1.0000000 8.459355e-164 4.229678e-164 [163,] 1.0000000 6.882020e-163 3.441010e-163 [164,] 1.0000000 5.403056e-162 2.701528e-162 [165,] 1.0000000 4.533571e-161 2.266785e-161 [166,] 1.0000000 3.562971e-160 1.781486e-160 [167,] 1.0000000 3.161754e-159 1.580877e-159 [168,] 1.0000000 2.490193e-158 1.245096e-158 [169,] 1.0000000 2.216461e-157 1.108231e-157 [170,] 1.0000000 1.918934e-156 9.594671e-157 [171,] 1.0000000 1.654776e-155 8.273881e-156 [172,] 1.0000000 9.287793e-155 4.643897e-155 [173,] 1.0000000 7.685498e-154 3.842749e-154 [174,] 1.0000000 4.331669e-153 2.165834e-153 [175,] 1.0000000 1.951747e-152 9.758736e-153 [176,] 1.0000000 1.707014e-151 8.535070e-152 [177,] 1.0000000 1.163927e-150 5.819633e-151 [178,] 1.0000000 4.752783e-150 2.376391e-150 [179,] 1.0000000 4.055371e-149 2.027686e-149 [180,] 1.0000000 2.714794e-148 1.357397e-148 [181,] 1.0000000 2.255632e-147 1.127816e-147 [182,] 1.0000000 1.785160e-146 8.925800e-147 [183,] 1.0000000 1.471050e-145 7.355251e-146 [184,] 1.0000000 1.223294e-144 6.116472e-145 [185,] 1.0000000 9.133870e-144 4.566935e-144 [186,] 1.0000000 7.461621e-143 3.730811e-143 [187,] 1.0000000 6.190266e-142 3.095133e-142 [188,] 1.0000000 4.034319e-141 2.017159e-141 [189,] 1.0000000 3.314928e-140 1.657464e-140 [190,] 1.0000000 2.710295e-139 1.355147e-139 [191,] 1.0000000 2.282429e-139 1.141215e-139 [192,] 1.0000000 1.874873e-138 9.374367e-139 [193,] 1.0000000 1.532662e-137 7.663309e-138 [194,] 1.0000000 1.246856e-136 6.234280e-137 [195,] 1.0000000 9.342345e-136 4.671172e-136 [196,] 1.0000000 3.479348e-135 1.739674e-135 [197,] 1.0000000 2.812736e-134 1.406368e-134 [198,] 1.0000000 2.230367e-133 1.115183e-133 [199,] 1.0000000 1.781902e-132 8.909511e-133 [200,] 1.0000000 3.024369e-132 1.512184e-132 [201,] 1.0000000 2.350495e-131 1.175248e-131 [202,] 1.0000000 1.811408e-130 9.057038e-131 [203,] 1.0000000 1.294973e-129 6.474867e-130 [204,] 1.0000000 6.814186e-129 3.407093e-129 [205,] 1.0000000 4.776280e-128 2.388140e-128 [206,] 1.0000000 3.562670e-127 1.781335e-127 [207,] 1.0000000 2.540584e-126 1.270292e-126 [208,] 1.0000000 1.923010e-125 9.615050e-126 [209,] 1.0000000 1.403194e-124 7.015970e-125 [210,] 1.0000000 1.049727e-123 5.248633e-124 [211,] 1.0000000 4.658856e-123 2.329428e-123 [212,] 1.0000000 3.354149e-122 1.677075e-122 [213,] 1.0000000 2.341114e-121 1.170557e-121 [214,] 1.0000000 1.632779e-120 8.163896e-121 [215,] 1.0000000 1.202557e-119 6.012784e-120 [216,] 1.0000000 8.760590e-119 4.380295e-119 [217,] 1.0000000 6.409275e-118 3.204637e-118 [218,] 1.0000000 4.658212e-117 2.329106e-117 [219,] 1.0000000 3.182352e-116 1.591176e-116 [220,] 1.0000000 2.222739e-115 1.111370e-115 [221,] 1.0000000 1.583412e-114 7.917060e-115 [222,] 1.0000000 1.066319e-113 5.331593e-114 [223,] 1.0000000 6.614849e-113 3.307424e-113 [224,] 1.0000000 4.453289e-112 2.226645e-112 [225,] 1.0000000 3.197515e-111 1.598758e-111 [226,] 1.0000000 2.230605e-110 1.115302e-110 [227,] 1.0000000 1.548521e-109 7.742604e-110 [228,] 1.0000000 1.070353e-108 5.351767e-109 [229,] 1.0000000 7.320488e-108 3.660244e-108 [230,] 1.0000000 4.810673e-107 2.405337e-107 [231,] 1.0000000 3.151298e-106 1.575649e-106 [232,] 1.0000000 2.139560e-105 1.069780e-105 [233,] 1.0000000 1.418807e-104 7.094035e-105 [234,] 1.0000000 9.545832e-104 4.772916e-104 [235,] 1.0000000 6.111366e-103 3.055683e-103 [236,] 1.0000000 4.070336e-102 2.035168e-102 [237,] 1.0000000 2.769436e-101 1.384718e-101 [238,] 1.0000000 1.875004e-100 9.375022e-101 [239,] 1.0000000 1.264913e-99 6.324564e-100 [240,] 1.0000000 8.480102e-99 4.240051e-99 [241,] 1.0000000 5.616045e-98 2.808022e-98 [242,] 1.0000000 3.508506e-97 1.754253e-97 [243,] 1.0000000 2.318360e-96 1.159180e-96 [244,] 1.0000000 1.467934e-95 7.339668e-96 [245,] 1.0000000 7.787919e-95 3.893960e-95 [246,] 1.0000000 5.870566e-95 2.935283e-95 [247,] 1.0000000 3.786195e-94 1.893098e-94 [248,] 1.0000000 2.229166e-93 1.114583e-93 [249,] 1.0000000 1.181037e-92 5.905187e-93 [250,] 1.0000000 7.659204e-92 3.829602e-92 [251,] 1.0000000 4.678480e-91 2.339240e-91 [252,] 1.0000000 3.004381e-90 1.502190e-90 [253,] 1.0000000 1.919313e-89 9.596566e-90 [254,] 1.0000000 1.211357e-88 6.056786e-89 [255,] 1.0000000 5.657238e-88 2.828619e-88 [256,] 1.0000000 3.566104e-87 1.783052e-87 [257,] 1.0000000 2.120619e-86 1.060310e-86 [258,] 1.0000000 1.319610e-85 6.598048e-86 [259,] 1.0000000 6.234355e-85 3.117177e-85 [260,] 1.0000000 3.856098e-84 1.928049e-84 [261,] 1.0000000 2.356679e-83 1.178340e-83 [262,] 1.0000000 1.442532e-82 7.212658e-83 [263,] 1.0000000 8.781296e-82 4.390648e-82 [264,] 1.0000000 5.315850e-81 2.657925e-81 [265,] 1.0000000 3.213512e-80 1.606756e-80 [266,] 1.0000000 1.896571e-79 9.482857e-80 [267,] 1.0000000 1.076193e-78 5.380964e-79 [268,] 1.0000000 6.368686e-78 3.184343e-78 [269,] 1.0000000 3.581063e-77 1.790531e-77 [270,] 1.0000000 1.803531e-76 9.017655e-77 [271,] 1.0000000 1.014375e-75 5.071875e-76 [272,] 1.0000000 5.917060e-75 2.958530e-75 [273,] 1.0000000 3.376103e-74 1.688052e-74 [274,] 1.0000000 1.753947e-73 8.769736e-74 [275,] 1.0000000 9.910888e-73 4.955444e-73 [276,] 1.0000000 5.585390e-72 2.792695e-72 [277,] 1.0000000 3.030375e-71 1.515188e-71 [278,] 1.0000000 1.652394e-70 8.261968e-71 [279,] 1.0000000 9.183571e-70 4.591786e-70 [280,] 1.0000000 5.137785e-69 2.568892e-69 [281,] 1.0000000 2.856848e-68 1.428424e-68 [282,] 1.0000000 1.580626e-67 7.903132e-68 [283,] 1.0000000 8.394063e-67 4.197031e-67 [284,] 1.0000000 4.450574e-66 2.225287e-66 [285,] 1.0000000 2.387072e-65 1.193536e-65 [286,] 1.0000000 1.285730e-64 6.428651e-65 [287,] 1.0000000 5.203901e-64 2.601950e-64 [288,] 1.0000000 2.753913e-63 1.376956e-63 [289,] 1.0000000 1.467843e-62 7.339214e-63 [290,] 1.0000000 7.467417e-62 3.733709e-62 [291,] 1.0000000 3.901597e-61 1.950799e-61 [292,] 1.0000000 2.049672e-60 1.024836e-60 [293,] 1.0000000 1.056153e-59 5.280763e-60 [294,] 1.0000000 5.263601e-59 2.631801e-59 [295,] 1.0000000 2.700465e-58 1.350232e-58 [296,] 1.0000000 1.068365e-57 5.341826e-58 [297,] 1.0000000 5.405417e-57 2.702708e-57 [298,] 1.0000000 2.732371e-56 1.366186e-56 [299,] 1.0000000 1.385621e-55 6.928106e-56 [300,] 1.0000000 6.338311e-55 3.169156e-55 [301,] 1.0000000 3.147763e-54 1.573881e-54 [302,] 1.0000000 1.550330e-53 7.751648e-54 [303,] 1.0000000 6.913106e-53 3.456553e-53 [304,] 1.0000000 3.362089e-52 1.681044e-52 [305,] 1.0000000 1.631869e-51 8.159347e-52 [306,] 1.0000000 7.834841e-51 3.917420e-51 [307,] 1.0000000 3.519905e-50 1.759952e-50 [308,] 1.0000000 1.652140e-49 8.260698e-50 [309,] 1.0000000 7.913803e-49 3.956902e-49 [310,] 1.0000000 3.703998e-48 1.851999e-48 [311,] 1.0000000 1.718858e-47 8.594288e-48 [312,] 1.0000000 7.789872e-47 3.894936e-47 [313,] 1.0000000 3.508059e-46 1.754030e-46 [314,] 1.0000000 1.569654e-45 7.848270e-46 [315,] 1.0000000 7.154481e-45 3.577240e-45 [316,] 1.0000000 3.162775e-44 1.581388e-44 [317,] 1.0000000 1.405488e-43 7.027441e-44 [318,] 1.0000000 6.272830e-43 3.136415e-43 [319,] 1.0000000 2.748366e-42 1.374183e-42 [320,] 1.0000000 1.183991e-41 5.919953e-42 [321,] 1.0000000 5.064765e-41 2.532382e-41 [322,] 1.0000000 2.151038e-40 1.075519e-40 [323,] 1.0000000 9.281042e-40 4.640521e-40 [324,] 1.0000000 3.889104e-39 1.944552e-39 [325,] 1.0000000 1.622848e-38 8.114239e-39 [326,] 1.0000000 6.870085e-38 3.435043e-38 [327,] 1.0000000 2.862942e-37 1.431471e-37 [328,] 1.0000000 1.222146e-36 6.110729e-37 [329,] 1.0000000 5.059980e-36 2.529990e-36 [330,] 1.0000000 2.092263e-35 1.046132e-35 [331,] 1.0000000 8.173019e-35 4.086509e-35 [332,] 1.0000000 3.352694e-34 1.676347e-34 [333,] 1.0000000 1.350493e-33 6.752466e-34 [334,] 1.0000000 5.397024e-33 2.698512e-33 [335,] 1.0000000 2.097327e-32 1.048664e-32 [336,] 1.0000000 8.345836e-32 4.172918e-32 [337,] 1.0000000 3.259795e-31 1.629897e-31 [338,] 1.0000000 1.287244e-30 6.436219e-31 [339,] 1.0000000 4.945802e-30 2.472901e-30 [340,] 1.0000000 1.854176e-29 9.270879e-30 [341,] 1.0000000 6.385035e-29 3.192517e-29 [342,] 1.0000000 2.217132e-28 1.108566e-28 [343,] 1.0000000 8.333414e-28 4.166707e-28 [344,] 1.0000000 2.719755e-27 1.359877e-27 [345,] 1.0000000 1.010475e-26 5.052376e-27 [346,] 1.0000000 3.672330e-26 1.836165e-26 [347,] 1.0000000 1.321805e-25 6.609023e-26 [348,] 1.0000000 4.695365e-25 2.347683e-25 [349,] 1.0000000 1.665203e-24 8.326014e-25 [350,] 1.0000000 5.845746e-24 2.922873e-24 [351,] 1.0000000 2.012322e-23 1.006161e-23 [352,] 1.0000000 6.948032e-23 3.474016e-23 [353,] 1.0000000 2.349383e-22 1.174691e-22 [354,] 1.0000000 7.948465e-22 3.974233e-22 [355,] 1.0000000 2.632301e-21 1.316150e-21 [356,] 1.0000000 8.492560e-21 4.246280e-21 [357,] 1.0000000 2.762021e-20 1.381010e-20 [358,] 1.0000000 8.880549e-20 4.440274e-20 [359,] 1.0000000 2.880110e-19 1.440055e-19 [360,] 1.0000000 9.075136e-19 4.537568e-19 [361,] 1.0000000 2.871064e-18 1.435532e-18 [362,] 1.0000000 8.862034e-18 4.431017e-18 [363,] 1.0000000 2.504471e-17 1.252235e-17 [364,] 1.0000000 8.090228e-17 4.045114e-17 [365,] 1.0000000 2.556178e-16 1.278089e-16 [366,] 1.0000000 7.517475e-16 3.758738e-16 [367,] 1.0000000 2.177930e-15 1.088965e-15 [368,] 1.0000000 6.462148e-15 3.231074e-15 [369,] 1.0000000 1.804762e-14 9.023808e-15 [370,] 1.0000000 4.659063e-14 2.329531e-14 [371,] 1.0000000 1.291315e-13 6.456577e-14 [372,] 1.0000000 3.515439e-13 1.757719e-13 [373,] 1.0000000 1.041513e-12 5.207565e-13 [374,] 1.0000000 2.588360e-12 1.294180e-12 [375,] 1.0000000 6.816873e-12 3.408437e-12 [376,] 1.0000000 1.859086e-11 9.295429e-12 [377,] 1.0000000 4.744356e-11 2.372178e-11 [378,] 1.0000000 1.142402e-10 5.712012e-11 [379,] 1.0000000 2.817406e-10 1.408703e-10 [380,] 1.0000000 7.112249e-10 3.556124e-10 [381,] 1.0000000 1.732707e-09 8.663534e-10 [382,] 1.0000000 4.091044e-09 2.045522e-09 [383,] 1.0000000 1.051183e-08 5.255917e-09 [384,] 1.0000000 2.757690e-08 1.378845e-08 [385,] 1.0000000 5.064842e-08 2.532421e-08 [386,] 0.9999999 1.137411e-07 5.687056e-08 [387,] 0.9999999 2.775260e-07 1.387630e-07 [388,] 0.9999997 6.440629e-07 3.220315e-07 [389,] 0.9999994 1.206373e-06 6.031866e-07 [390,] 0.9999986 2.823280e-06 1.411640e-06 [391,] 0.9999971 5.876508e-06 2.938254e-06 [392,] 0.9999931 1.373613e-05 6.868066e-06 [393,] 0.9999865 2.694751e-05 1.347375e-05 [394,] 0.9999719 5.613207e-05 2.806604e-05 [395,] 0.9999393 1.214554e-04 6.072768e-05 [396,] 0.9998668 2.664745e-04 1.332373e-04 [397,] 0.9997393 5.214706e-04 2.607353e-04 [398,] 0.9994818 1.036438e-03 5.182188e-04 [399,] 0.9990712 1.857575e-03 9.287875e-04 [400,] 0.9983556 3.288838e-03 1.644419e-03 [401,] 0.9967621 6.475853e-03 3.237926e-03 [402,] 0.9938543 1.229134e-02 6.145668e-03 [403,] 0.9887839 2.243216e-02 1.121608e-02 [404,] 0.9793192 4.136155e-02 2.068078e-02 [405,] 0.9810969 3.780613e-02 1.890307e-02 [406,] 0.9692122 6.157567e-02 3.078784e-02 [407,] 0.9489260 1.021480e-01 5.107402e-02 [408,] 0.9843747 3.125068e-02 1.562534e-02 [409,] 0.9668674 6.626526e-02 3.313263e-02 [410,] 0.9315350 1.369300e-01 6.846500e-02 [411,] 0.8669904 2.660191e-01 1.330096e-01 [412,] 0.7595129 4.809743e-01 2.404871e-01 [413,] 0.6009843 7.980313e-01 3.990157e-01 > postscript(file="/var/www/html/rcomp/tmp/13jx21292932180.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/23jx21292932180.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/33jx21292932180.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/4esxn1292932180.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/5esxn1292932180.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 = 430 Frequency = 1 1 2 3 4 5 2.945709e+06 1.958932e+06 -1.608369e+06 -6.533468e+04 3.622083e+05 6 7 8 9 10 1.569258e+05 3.729608e+05 -1.382227e+05 7.579032e+05 1.030653e+05 11 12 13 14 15 4.217010e+05 3.683217e+05 -8.069330e+05 -1.555096e+05 4.914049e+05 16 17 18 19 20 -1.921223e+05 -1.947991e+05 -1.065107e+05 -5.178519e+05 2.502813e+05 21 22 23 24 25 -4.255778e+04 -4.201810e+04 -5.939069e+05 -2.534845e+05 6.050526e+04 26 27 28 29 30 -3.689886e+05 -2.475860e+05 4.290343e+05 -1.733252e+05 2.672605e+05 31 32 33 34 35 -4.558695e+04 2.200257e+05 1.174329e+05 3.411520e+04 -4.667818e+05 36 37 38 39 40 1.238381e+05 1.356866e+05 -8.574527e+04 -3.037759e+05 -4.832529e+05 41 42 43 44 45 -1.970856e+05 -1.326601e+04 -1.952668e+05 -2.031444e+04 -6.659535e+04 46 47 48 49 50 -2.718181e+05 -2.949503e+05 -6.060906e+01 -1.543829e+05 -1.271625e+04 51 52 53 54 55 -1.427621e+05 -2.242160e+05 -4.412895e+04 -4.303624e+05 -1.997302e+04 56 57 58 59 60 3.734755e+05 -8.144209e+03 -4.126198e+04 -1.979514e+05 -3.327342e+04 61 62 63 64 65 -1.996854e+05 -1.262337e+04 -2.224404e+05 -3.673721e+04 -4.644577e+05 66 67 68 69 70 -3.685077e+04 -5.278891e+04 -8.535061e+04 -1.114192e+05 -1.687591e+05 71 72 73 74 75 6.897109e+03 -1.018463e+05 4.283796e+04 2.427217e+03 -1.322708e+05 76 77 78 79 80 -4.093145e+04 -8.739013e+04 -3.093746e+04 -6.970571e+04 4.885414e+04 81 82 83 84 85 -1.180068e+05 -5.111407e+04 1.831975e+04 5.431579e+04 -1.254997e+05 86 87 88 89 90 -1.506956e+04 -1.013768e+04 -1.270529e+05 2.083809e+04 -1.259173e+04 91 92 93 94 95 -5.295845e+04 8.081311e+04 -1.550252e+05 6.448599e+04 -5.258125e+04 96 97 98 99 100 1.040454e+05 -3.726888e+04 -3.086653e+04 -3.412860e+04 -7.991547e+04 101 102 103 104 105 -8.335505e+04 -1.278187e+05 -1.436685e+04 5.683226e+04 -3.521946e+04 106 107 108 109 110 -2.898622e+05 1.064864e+05 -1.553016e+05 -1.252069e+05 1.170282e+05 111 112 113 114 115 -2.213077e+05 -1.524851e+05 -8.239323e+04 5.174149e+04 -1.529056e+03 116 117 118 119 120 3.358610e+03 7.191914e+04 -1.569444e+05 3.564294e+04 -2.009287e+05 121 122 123 124 125 1.732224e+04 -1.134618e+05 5.424982e+04 -1.295347e+04 -2.256727e+05 126 127 128 129 130 3.600235e+04 1.885287e+04 2.591166e+03 7.816804e+04 2.415626e+04 131 132 133 134 135 8.192810e+04 7.666630e+04 2.664539e+03 1.855164e+04 2.592867e+04 136 137 138 139 140 6.491739e+04 3.394453e+04 1.820984e+04 3.600235e+04 3.155258e+04 141 142 143 144 145 7.816804e+04 -4.587239e+03 2.975914e+04 6.067965e+04 2.176231e+04 146 147 148 149 150 7.816804e+04 2.482765e+04 8.161402e+04 7.808688e+04 -1.993497e+03 151 152 153 154 155 -6.439980e+04 2.849368e+04 4.710343e+04 3.231582e+04 3.600235e+04 156 157 158 159 160 3.600235e+04 3.057385e+04 4.425873e+04 3.600235e+04 3.600235e+04 161 162 163 164 165 7.816804e+04 3.600235e+04 1.013438e+04 1.975804e+04 3.340956e+04 166 167 168 169 170 7.816804e+04 1.293051e+05 -3.697539e+04 3.600235e+04 3.600235e+04 171 172 173 174 175 6.909576e+04 7.816804e+04 5.465479e+04 7.816804e+04 2.225224e+04 176 177 178 179 180 -1.408540e+05 1.114887e+04 3.600235e+04 3.600235e+04 6.651030e+04 181 182 183 184 185 -5.422786e+04 -2.700112e+04 -9.277563e+04 4.207503e+02 5.017889e+03 186 187 188 189 190 -6.335012e+04 2.127199e+04 3.739862e+04 3.600235e+04 2.704002e+04 191 192 193 194 195 3.600235e+04 2.246269e+04 7.816804e+04 3.122124e+04 1.128268e+04 196 197 198 199 200 -1.625690e+05 1.128268e+04 1.128268e+04 -9.521193e+04 1.128268e+04 201 202 203 204 205 1.128268e+04 1.128268e+04 -5.243307e+04 -3.907353e+04 3.153210e+03 206 207 208 209 210 3.939561e+03 1.128268e+04 -3.948385e+04 -4.716405e+04 3.600235e+04 211 212 213 214 215 7.816804e+04 7.713665e+04 2.906961e+04 5.044902e+04 7.516457e+04 216 217 218 219 220 2.696010e+04 2.293713e+04 3.299888e+04 -1.654379e+04 5.158962e+04 221 222 223 224 225 7.816804e+04 7.663229e+04 -3.156698e+03 3.432896e+04 -1.929492e+03 226 227 228 229 230 2.088802e+04 7.816804e+04 6.352101e+04 3.600235e+04 4.264214e+04 231 232 233 234 235 6.133615e+04 7.816804e+04 -1.624027e+03 3.600235e+04 3.605292e+04 236 237 238 239 240 3.600235e+04 2.723941e+04 7.816804e+04 7.816804e+04 3.600235e+04 241 242 243 244 245 6.911561e+04 3.600235e+04 6.883354e+04 3.852101e+04 1.128268e+04 246 247 248 249 250 1.128268e+04 3.956546e+03 1.128268e+04 -3.088301e+04 5.369742e+04 251 252 253 254 255 1.128268e+04 6.171967e+04 -5.658248e+04 -1.706214e+05 3.600235e+04 256 257 258 259 260 2.306742e+04 -6.896001e+03 1.128268e+04 7.816804e+04 1.128268e+04 261 262 263 264 265 1.128268e+04 2.574404e+04 7.838025e+04 1.128268e+04 7.816804e+04 266 267 268 269 270 -3.088301e+04 1.051092e+05 1.128268e+04 4.108224e+04 1.128268e+04 271 272 273 274 275 1.128268e+04 1.128268e+04 6.072985e+02 3.600235e+04 5.887555e+04 276 277 278 279 280 -2.433108e+03 7.816804e+04 -8.162836e+04 4.916413e+04 1.128268e+04 281 282 283 284 285 3.600235e+04 3.950952e+04 3.600235e+04 6.270590e+04 7.816804e+04 286 287 288 289 290 2.043001e+01 3.203934e+04 -3.088301e+04 -6.111311e+02 1.128268e+04 291 292 293 294 295 2.281855e+04 6.682024e+04 3.600235e+04 5.115761e+04 2.037284e+04 296 297 298 299 300 3.600235e+04 -3.088301e+04 7.816804e+04 2.248674e+04 2.575388e+03 301 302 303 304 305 6.140890e+04 7.666630e+04 1.128268e+04 3.792899e+04 1.794995e+04 306 307 308 309 310 1.128268e+04 -2.309127e+04 8.272152e+03 1.128268e+04 1.128268e+04 311 312 313 314 315 2.083068e+04 1.128268e+04 -5.046879e+04 2.417100e+04 4.523385e+04 316 317 318 319 320 3.617547e+04 -1.752119e+04 2.847441e+04 3.600235e+04 7.816804e+04 321 322 323 324 325 7.816804e+04 7.816804e+04 -3.088301e+04 7.816804e+04 7.350799e+04 326 327 328 329 330 -3.088301e+04 1.128268e+04 7.816804e+04 7.816804e+04 7.816804e+04 331 332 333 334 335 3.600235e+04 7.816804e+04 6.976128e+04 3.448651e+04 2.602428e+04 336 337 338 339 340 1.796446e+04 3.600235e+04 3.243381e+04 6.522883e+04 -8.705010e+03 341 342 343 344 345 3.600235e+04 3.600235e+04 7.816804e+04 1.519857e+03 3.600235e+04 346 347 348 349 350 2.046556e+04 3.600235e+04 7.816804e+04 4.847882e+03 7.591426e+04 351 352 353 354 355 2.285412e+04 6.758874e+04 2.980277e+04 3.600235e+04 3.600235e+04 356 357 358 359 360 -1.050841e+04 3.600235e+04 3.600235e+04 7.816804e+04 2.999808e+04 361 362 363 364 365 7.816804e+04 3.600235e+04 -6.124955e+04 6.485885e+04 7.816804e+04 366 367 368 369 370 7.816804e+04 3.600235e+04 7.816804e+04 2.291243e+04 7.816804e+04 371 372 373 374 375 7.125759e+04 2.228522e+04 5.794993e+04 7.816804e+04 7.816804e+04 376 377 378 379 380 4.627300e+04 2.517615e+04 -2.442581e+03 7.816804e+04 7.816804e+04 381 382 383 384 385 -1.304973e+04 -1.107369e+03 7.816804e+04 1.626159e+04 -8.365301e+04 386 387 388 389 390 4.297161e+04 5.401746e+04 2.985540e+04 5.438925e+04 6.238896e+04 391 392 393 394 395 -3.806724e+01 5.889624e+03 6.055375e+04 -1.501482e+05 3.034063e+04 396 397 398 399 400 -8.392311e+04 -3.247143e+05 3.758911e+04 -1.079982e+05 -6.169317e+03 401 402 403 404 405 2.756067e+04 2.896510e+03 -1.297722e+05 -2.527709e+04 2.666550e+03 406 407 408 409 410 4.559749e+04 3.927075e+04 2.931348e+03 -3.157747e+04 -4.223222e+04 411 412 413 414 415 -1.457241e+05 -6.694617e+04 3.464206e+04 -1.379192e+05 -1.075542e+05 416 417 418 419 420 1.288701e+05 -1.022662e+05 -6.764406e+04 -4.973243e+04 -5.246928e+04 421 422 423 424 425 -1.595566e+05 -5.328526e+05 -3.851035e+05 -1.557869e+04 -2.849142e+05 426 427 428 429 430 -2.636122e+05 -8.778464e+05 -1.739617e+05 -2.086613e+05 -5.306982e+05 > postscript(file="/var/www/html/rcomp/tmp/6esxn1292932180.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 = 430 Frequency = 1 lag(myerror, k = 1) myerror 0 2.945709e+06 NA 1 1.958932e+06 2.945709e+06 2 -1.608369e+06 1.958932e+06 3 -6.533468e+04 -1.608369e+06 4 3.622083e+05 -6.533468e+04 5 1.569258e+05 3.622083e+05 6 3.729608e+05 1.569258e+05 7 -1.382227e+05 3.729608e+05 8 7.579032e+05 -1.382227e+05 9 1.030653e+05 7.579032e+05 10 4.217010e+05 1.030653e+05 11 3.683217e+05 4.217010e+05 12 -8.069330e+05 3.683217e+05 13 -1.555096e+05 -8.069330e+05 14 4.914049e+05 -1.555096e+05 15 -1.921223e+05 4.914049e+05 16 -1.947991e+05 -1.921223e+05 17 -1.065107e+05 -1.947991e+05 18 -5.178519e+05 -1.065107e+05 19 2.502813e+05 -5.178519e+05 20 -4.255778e+04 2.502813e+05 21 -4.201810e+04 -4.255778e+04 22 -5.939069e+05 -4.201810e+04 23 -2.534845e+05 -5.939069e+05 24 6.050526e+04 -2.534845e+05 25 -3.689886e+05 6.050526e+04 26 -2.475860e+05 -3.689886e+05 27 4.290343e+05 -2.475860e+05 28 -1.733252e+05 4.290343e+05 29 2.672605e+05 -1.733252e+05 30 -4.558695e+04 2.672605e+05 31 2.200257e+05 -4.558695e+04 32 1.174329e+05 2.200257e+05 33 3.411520e+04 1.174329e+05 34 -4.667818e+05 3.411520e+04 35 1.238381e+05 -4.667818e+05 36 1.356866e+05 1.238381e+05 37 -8.574527e+04 1.356866e+05 38 -3.037759e+05 -8.574527e+04 39 -4.832529e+05 -3.037759e+05 40 -1.970856e+05 -4.832529e+05 41 -1.326601e+04 -1.970856e+05 42 -1.952668e+05 -1.326601e+04 43 -2.031444e+04 -1.952668e+05 44 -6.659535e+04 -2.031444e+04 45 -2.718181e+05 -6.659535e+04 46 -2.949503e+05 -2.718181e+05 47 -6.060906e+01 -2.949503e+05 48 -1.543829e+05 -6.060906e+01 49 -1.271625e+04 -1.543829e+05 50 -1.427621e+05 -1.271625e+04 51 -2.242160e+05 -1.427621e+05 52 -4.412895e+04 -2.242160e+05 53 -4.303624e+05 -4.412895e+04 54 -1.997302e+04 -4.303624e+05 55 3.734755e+05 -1.997302e+04 56 -8.144209e+03 3.734755e+05 57 -4.126198e+04 -8.144209e+03 58 -1.979514e+05 -4.126198e+04 59 -3.327342e+04 -1.979514e+05 60 -1.996854e+05 -3.327342e+04 61 -1.262337e+04 -1.996854e+05 62 -2.224404e+05 -1.262337e+04 63 -3.673721e+04 -2.224404e+05 64 -4.644577e+05 -3.673721e+04 65 -3.685077e+04 -4.644577e+05 66 -5.278891e+04 -3.685077e+04 67 -8.535061e+04 -5.278891e+04 68 -1.114192e+05 -8.535061e+04 69 -1.687591e+05 -1.114192e+05 70 6.897109e+03 -1.687591e+05 71 -1.018463e+05 6.897109e+03 72 4.283796e+04 -1.018463e+05 73 2.427217e+03 4.283796e+04 74 -1.322708e+05 2.427217e+03 75 -4.093145e+04 -1.322708e+05 76 -8.739013e+04 -4.093145e+04 77 -3.093746e+04 -8.739013e+04 78 -6.970571e+04 -3.093746e+04 79 4.885414e+04 -6.970571e+04 80 -1.180068e+05 4.885414e+04 81 -5.111407e+04 -1.180068e+05 82 1.831975e+04 -5.111407e+04 83 5.431579e+04 1.831975e+04 84 -1.254997e+05 5.431579e+04 85 -1.506956e+04 -1.254997e+05 86 -1.013768e+04 -1.506956e+04 87 -1.270529e+05 -1.013768e+04 88 2.083809e+04 -1.270529e+05 89 -1.259173e+04 2.083809e+04 90 -5.295845e+04 -1.259173e+04 91 8.081311e+04 -5.295845e+04 92 -1.550252e+05 8.081311e+04 93 6.448599e+04 -1.550252e+05 94 -5.258125e+04 6.448599e+04 95 1.040454e+05 -5.258125e+04 96 -3.726888e+04 1.040454e+05 97 -3.086653e+04 -3.726888e+04 98 -3.412860e+04 -3.086653e+04 99 -7.991547e+04 -3.412860e+04 100 -8.335505e+04 -7.991547e+04 101 -1.278187e+05 -8.335505e+04 102 -1.436685e+04 -1.278187e+05 103 5.683226e+04 -1.436685e+04 104 -3.521946e+04 5.683226e+04 105 -2.898622e+05 -3.521946e+04 106 1.064864e+05 -2.898622e+05 107 -1.553016e+05 1.064864e+05 108 -1.252069e+05 -1.553016e+05 109 1.170282e+05 -1.252069e+05 110 -2.213077e+05 1.170282e+05 111 -1.524851e+05 -2.213077e+05 112 -8.239323e+04 -1.524851e+05 113 5.174149e+04 -8.239323e+04 114 -1.529056e+03 5.174149e+04 115 3.358610e+03 -1.529056e+03 116 7.191914e+04 3.358610e+03 117 -1.569444e+05 7.191914e+04 118 3.564294e+04 -1.569444e+05 119 -2.009287e+05 3.564294e+04 120 1.732224e+04 -2.009287e+05 121 -1.134618e+05 1.732224e+04 122 5.424982e+04 -1.134618e+05 123 -1.295347e+04 5.424982e+04 124 -2.256727e+05 -1.295347e+04 125 3.600235e+04 -2.256727e+05 126 1.885287e+04 3.600235e+04 127 2.591166e+03 1.885287e+04 128 7.816804e+04 2.591166e+03 129 2.415626e+04 7.816804e+04 130 8.192810e+04 2.415626e+04 131 7.666630e+04 8.192810e+04 132 2.664539e+03 7.666630e+04 133 1.855164e+04 2.664539e+03 134 2.592867e+04 1.855164e+04 135 6.491739e+04 2.592867e+04 136 3.394453e+04 6.491739e+04 137 1.820984e+04 3.394453e+04 138 3.600235e+04 1.820984e+04 139 3.155258e+04 3.600235e+04 140 7.816804e+04 3.155258e+04 141 -4.587239e+03 7.816804e+04 142 2.975914e+04 -4.587239e+03 143 6.067965e+04 2.975914e+04 144 2.176231e+04 6.067965e+04 145 7.816804e+04 2.176231e+04 146 2.482765e+04 7.816804e+04 147 8.161402e+04 2.482765e+04 148 7.808688e+04 8.161402e+04 149 -1.993497e+03 7.808688e+04 150 -6.439980e+04 -1.993497e+03 151 2.849368e+04 -6.439980e+04 152 4.710343e+04 2.849368e+04 153 3.231582e+04 4.710343e+04 154 3.600235e+04 3.231582e+04 155 3.600235e+04 3.600235e+04 156 3.057385e+04 3.600235e+04 157 4.425873e+04 3.057385e+04 158 3.600235e+04 4.425873e+04 159 3.600235e+04 3.600235e+04 160 7.816804e+04 3.600235e+04 161 3.600235e+04 7.816804e+04 162 1.013438e+04 3.600235e+04 163 1.975804e+04 1.013438e+04 164 3.340956e+04 1.975804e+04 165 7.816804e+04 3.340956e+04 166 1.293051e+05 7.816804e+04 167 -3.697539e+04 1.293051e+05 168 3.600235e+04 -3.697539e+04 169 3.600235e+04 3.600235e+04 170 6.909576e+04 3.600235e+04 171 7.816804e+04 6.909576e+04 172 5.465479e+04 7.816804e+04 173 7.816804e+04 5.465479e+04 174 2.225224e+04 7.816804e+04 175 -1.408540e+05 2.225224e+04 176 1.114887e+04 -1.408540e+05 177 3.600235e+04 1.114887e+04 178 3.600235e+04 3.600235e+04 179 6.651030e+04 3.600235e+04 180 -5.422786e+04 6.651030e+04 181 -2.700112e+04 -5.422786e+04 182 -9.277563e+04 -2.700112e+04 183 4.207503e+02 -9.277563e+04 184 5.017889e+03 4.207503e+02 185 -6.335012e+04 5.017889e+03 186 2.127199e+04 -6.335012e+04 187 3.739862e+04 2.127199e+04 188 3.600235e+04 3.739862e+04 189 2.704002e+04 3.600235e+04 190 3.600235e+04 2.704002e+04 191 2.246269e+04 3.600235e+04 192 7.816804e+04 2.246269e+04 193 3.122124e+04 7.816804e+04 194 1.128268e+04 3.122124e+04 195 -1.625690e+05 1.128268e+04 196 1.128268e+04 -1.625690e+05 197 1.128268e+04 1.128268e+04 198 -9.521193e+04 1.128268e+04 199 1.128268e+04 -9.521193e+04 200 1.128268e+04 1.128268e+04 201 1.128268e+04 1.128268e+04 202 -5.243307e+04 1.128268e+04 203 -3.907353e+04 -5.243307e+04 204 3.153210e+03 -3.907353e+04 205 3.939561e+03 3.153210e+03 206 1.128268e+04 3.939561e+03 207 -3.948385e+04 1.128268e+04 208 -4.716405e+04 -3.948385e+04 209 3.600235e+04 -4.716405e+04 210 7.816804e+04 3.600235e+04 211 7.713665e+04 7.816804e+04 212 2.906961e+04 7.713665e+04 213 5.044902e+04 2.906961e+04 214 7.516457e+04 5.044902e+04 215 2.696010e+04 7.516457e+04 216 2.293713e+04 2.696010e+04 217 3.299888e+04 2.293713e+04 218 -1.654379e+04 3.299888e+04 219 5.158962e+04 -1.654379e+04 220 7.816804e+04 5.158962e+04 221 7.663229e+04 7.816804e+04 222 -3.156698e+03 7.663229e+04 223 3.432896e+04 -3.156698e+03 224 -1.929492e+03 3.432896e+04 225 2.088802e+04 -1.929492e+03 226 7.816804e+04 2.088802e+04 227 6.352101e+04 7.816804e+04 228 3.600235e+04 6.352101e+04 229 4.264214e+04 3.600235e+04 230 6.133615e+04 4.264214e+04 231 7.816804e+04 6.133615e+04 232 -1.624027e+03 7.816804e+04 233 3.600235e+04 -1.624027e+03 234 3.605292e+04 3.600235e+04 235 3.600235e+04 3.605292e+04 236 2.723941e+04 3.600235e+04 237 7.816804e+04 2.723941e+04 238 7.816804e+04 7.816804e+04 239 3.600235e+04 7.816804e+04 240 6.911561e+04 3.600235e+04 241 3.600235e+04 6.911561e+04 242 6.883354e+04 3.600235e+04 243 3.852101e+04 6.883354e+04 244 1.128268e+04 3.852101e+04 245 1.128268e+04 1.128268e+04 246 3.956546e+03 1.128268e+04 247 1.128268e+04 3.956546e+03 248 -3.088301e+04 1.128268e+04 249 5.369742e+04 -3.088301e+04 250 1.128268e+04 5.369742e+04 251 6.171967e+04 1.128268e+04 252 -5.658248e+04 6.171967e+04 253 -1.706214e+05 -5.658248e+04 254 3.600235e+04 -1.706214e+05 255 2.306742e+04 3.600235e+04 256 -6.896001e+03 2.306742e+04 257 1.128268e+04 -6.896001e+03 258 7.816804e+04 1.128268e+04 259 1.128268e+04 7.816804e+04 260 1.128268e+04 1.128268e+04 261 2.574404e+04 1.128268e+04 262 7.838025e+04 2.574404e+04 263 1.128268e+04 7.838025e+04 264 7.816804e+04 1.128268e+04 265 -3.088301e+04 7.816804e+04 266 1.051092e+05 -3.088301e+04 267 1.128268e+04 1.051092e+05 268 4.108224e+04 1.128268e+04 269 1.128268e+04 4.108224e+04 270 1.128268e+04 1.128268e+04 271 1.128268e+04 1.128268e+04 272 6.072985e+02 1.128268e+04 273 3.600235e+04 6.072985e+02 274 5.887555e+04 3.600235e+04 275 -2.433108e+03 5.887555e+04 276 7.816804e+04 -2.433108e+03 277 -8.162836e+04 7.816804e+04 278 4.916413e+04 -8.162836e+04 279 1.128268e+04 4.916413e+04 280 3.600235e+04 1.128268e+04 281 3.950952e+04 3.600235e+04 282 3.600235e+04 3.950952e+04 283 6.270590e+04 3.600235e+04 284 7.816804e+04 6.270590e+04 285 2.043001e+01 7.816804e+04 286 3.203934e+04 2.043001e+01 287 -3.088301e+04 3.203934e+04 288 -6.111311e+02 -3.088301e+04 289 1.128268e+04 -6.111311e+02 290 2.281855e+04 1.128268e+04 291 6.682024e+04 2.281855e+04 292 3.600235e+04 6.682024e+04 293 5.115761e+04 3.600235e+04 294 2.037284e+04 5.115761e+04 295 3.600235e+04 2.037284e+04 296 -3.088301e+04 3.600235e+04 297 7.816804e+04 -3.088301e+04 298 2.248674e+04 7.816804e+04 299 2.575388e+03 2.248674e+04 300 6.140890e+04 2.575388e+03 301 7.666630e+04 6.140890e+04 302 1.128268e+04 7.666630e+04 303 3.792899e+04 1.128268e+04 304 1.794995e+04 3.792899e+04 305 1.128268e+04 1.794995e+04 306 -2.309127e+04 1.128268e+04 307 8.272152e+03 -2.309127e+04 308 1.128268e+04 8.272152e+03 309 1.128268e+04 1.128268e+04 310 2.083068e+04 1.128268e+04 311 1.128268e+04 2.083068e+04 312 -5.046879e+04 1.128268e+04 313 2.417100e+04 -5.046879e+04 314 4.523385e+04 2.417100e+04 315 3.617547e+04 4.523385e+04 316 -1.752119e+04 3.617547e+04 317 2.847441e+04 -1.752119e+04 318 3.600235e+04 2.847441e+04 319 7.816804e+04 3.600235e+04 320 7.816804e+04 7.816804e+04 321 7.816804e+04 7.816804e+04 322 -3.088301e+04 7.816804e+04 323 7.816804e+04 -3.088301e+04 324 7.350799e+04 7.816804e+04 325 -3.088301e+04 7.350799e+04 326 1.128268e+04 -3.088301e+04 327 7.816804e+04 1.128268e+04 328 7.816804e+04 7.816804e+04 329 7.816804e+04 7.816804e+04 330 3.600235e+04 7.816804e+04 331 7.816804e+04 3.600235e+04 332 6.976128e+04 7.816804e+04 333 3.448651e+04 6.976128e+04 334 2.602428e+04 3.448651e+04 335 1.796446e+04 2.602428e+04 336 3.600235e+04 1.796446e+04 337 3.243381e+04 3.600235e+04 338 6.522883e+04 3.243381e+04 339 -8.705010e+03 6.522883e+04 340 3.600235e+04 -8.705010e+03 341 3.600235e+04 3.600235e+04 342 7.816804e+04 3.600235e+04 343 1.519857e+03 7.816804e+04 344 3.600235e+04 1.519857e+03 345 2.046556e+04 3.600235e+04 346 3.600235e+04 2.046556e+04 347 7.816804e+04 3.600235e+04 348 4.847882e+03 7.816804e+04 349 7.591426e+04 4.847882e+03 350 2.285412e+04 7.591426e+04 351 6.758874e+04 2.285412e+04 352 2.980277e+04 6.758874e+04 353 3.600235e+04 2.980277e+04 354 3.600235e+04 3.600235e+04 355 -1.050841e+04 3.600235e+04 356 3.600235e+04 -1.050841e+04 357 3.600235e+04 3.600235e+04 358 7.816804e+04 3.600235e+04 359 2.999808e+04 7.816804e+04 360 7.816804e+04 2.999808e+04 361 3.600235e+04 7.816804e+04 362 -6.124955e+04 3.600235e+04 363 6.485885e+04 -6.124955e+04 364 7.816804e+04 6.485885e+04 365 7.816804e+04 7.816804e+04 366 3.600235e+04 7.816804e+04 367 7.816804e+04 3.600235e+04 368 2.291243e+04 7.816804e+04 369 7.816804e+04 2.291243e+04 370 7.125759e+04 7.816804e+04 371 2.228522e+04 7.125759e+04 372 5.794993e+04 2.228522e+04 373 7.816804e+04 5.794993e+04 374 7.816804e+04 7.816804e+04 375 4.627300e+04 7.816804e+04 376 2.517615e+04 4.627300e+04 377 -2.442581e+03 2.517615e+04 378 7.816804e+04 -2.442581e+03 379 7.816804e+04 7.816804e+04 380 -1.304973e+04 7.816804e+04 381 -1.107369e+03 -1.304973e+04 382 7.816804e+04 -1.107369e+03 383 1.626159e+04 7.816804e+04 384 -8.365301e+04 1.626159e+04 385 4.297161e+04 -8.365301e+04 386 5.401746e+04 4.297161e+04 387 2.985540e+04 5.401746e+04 388 5.438925e+04 2.985540e+04 389 6.238896e+04 5.438925e+04 390 -3.806724e+01 6.238896e+04 391 5.889624e+03 -3.806724e+01 392 6.055375e+04 5.889624e+03 393 -1.501482e+05 6.055375e+04 394 3.034063e+04 -1.501482e+05 395 -8.392311e+04 3.034063e+04 396 -3.247143e+05 -8.392311e+04 397 3.758911e+04 -3.247143e+05 398 -1.079982e+05 3.758911e+04 399 -6.169317e+03 -1.079982e+05 400 2.756067e+04 -6.169317e+03 401 2.896510e+03 2.756067e+04 402 -1.297722e+05 2.896510e+03 403 -2.527709e+04 -1.297722e+05 404 2.666550e+03 -2.527709e+04 405 4.559749e+04 2.666550e+03 406 3.927075e+04 4.559749e+04 407 2.931348e+03 3.927075e+04 408 -3.157747e+04 2.931348e+03 409 -4.223222e+04 -3.157747e+04 410 -1.457241e+05 -4.223222e+04 411 -6.694617e+04 -1.457241e+05 412 3.464206e+04 -6.694617e+04 413 -1.379192e+05 3.464206e+04 414 -1.075542e+05 -1.379192e+05 415 1.288701e+05 -1.075542e+05 416 -1.022662e+05 1.288701e+05 417 -6.764406e+04 -1.022662e+05 418 -4.973243e+04 -6.764406e+04 419 -5.246928e+04 -4.973243e+04 420 -1.595566e+05 -5.246928e+04 421 -5.328526e+05 -1.595566e+05 422 -3.851035e+05 -5.328526e+05 423 -1.557869e+04 -3.851035e+05 424 -2.849142e+05 -1.557869e+04 425 -2.636122e+05 -2.849142e+05 426 -8.778464e+05 -2.636122e+05 427 -1.739617e+05 -8.778464e+05 428 -2.086613e+05 -1.739617e+05 429 -5.306982e+05 -2.086613e+05 430 NA -5.306982e+05 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 1.958932e+06 2.945709e+06 [2,] -1.608369e+06 1.958932e+06 [3,] -6.533468e+04 -1.608369e+06 [4,] 3.622083e+05 -6.533468e+04 [5,] 1.569258e+05 3.622083e+05 [6,] 3.729608e+05 1.569258e+05 [7,] -1.382227e+05 3.729608e+05 [8,] 7.579032e+05 -1.382227e+05 [9,] 1.030653e+05 7.579032e+05 [10,] 4.217010e+05 1.030653e+05 [11,] 3.683217e+05 4.217010e+05 [12,] -8.069330e+05 3.683217e+05 [13,] -1.555096e+05 -8.069330e+05 [14,] 4.914049e+05 -1.555096e+05 [15,] -1.921223e+05 4.914049e+05 [16,] -1.947991e+05 -1.921223e+05 [17,] -1.065107e+05 -1.947991e+05 [18,] -5.178519e+05 -1.065107e+05 [19,] 2.502813e+05 -5.178519e+05 [20,] -4.255778e+04 2.502813e+05 [21,] -4.201810e+04 -4.255778e+04 [22,] -5.939069e+05 -4.201810e+04 [23,] -2.534845e+05 -5.939069e+05 [24,] 6.050526e+04 -2.534845e+05 [25,] -3.689886e+05 6.050526e+04 [26,] -2.475860e+05 -3.689886e+05 [27,] 4.290343e+05 -2.475860e+05 [28,] -1.733252e+05 4.290343e+05 [29,] 2.672605e+05 -1.733252e+05 [30,] -4.558695e+04 2.672605e+05 [31,] 2.200257e+05 -4.558695e+04 [32,] 1.174329e+05 2.200257e+05 [33,] 3.411520e+04 1.174329e+05 [34,] -4.667818e+05 3.411520e+04 [35,] 1.238381e+05 -4.667818e+05 [36,] 1.356866e+05 1.238381e+05 [37,] -8.574527e+04 1.356866e+05 [38,] -3.037759e+05 -8.574527e+04 [39,] -4.832529e+05 -3.037759e+05 [40,] -1.970856e+05 -4.832529e+05 [41,] -1.326601e+04 -1.970856e+05 [42,] -1.952668e+05 -1.326601e+04 [43,] -2.031444e+04 -1.952668e+05 [44,] -6.659535e+04 -2.031444e+04 [45,] -2.718181e+05 -6.659535e+04 [46,] -2.949503e+05 -2.718181e+05 [47,] -6.060906e+01 -2.949503e+05 [48,] -1.543829e+05 -6.060906e+01 [49,] -1.271625e+04 -1.543829e+05 [50,] -1.427621e+05 -1.271625e+04 [51,] -2.242160e+05 -1.427621e+05 [52,] -4.412895e+04 -2.242160e+05 [53,] -4.303624e+05 -4.412895e+04 [54,] -1.997302e+04 -4.303624e+05 [55,] 3.734755e+05 -1.997302e+04 [56,] -8.144209e+03 3.734755e+05 [57,] -4.126198e+04 -8.144209e+03 [58,] -1.979514e+05 -4.126198e+04 [59,] -3.327342e+04 -1.979514e+05 [60,] -1.996854e+05 -3.327342e+04 [61,] -1.262337e+04 -1.996854e+05 [62,] -2.224404e+05 -1.262337e+04 [63,] -3.673721e+04 -2.224404e+05 [64,] -4.644577e+05 -3.673721e+04 [65,] -3.685077e+04 -4.644577e+05 [66,] -5.278891e+04 -3.685077e+04 [67,] -8.535061e+04 -5.278891e+04 [68,] -1.114192e+05 -8.535061e+04 [69,] -1.687591e+05 -1.114192e+05 [70,] 6.897109e+03 -1.687591e+05 [71,] -1.018463e+05 6.897109e+03 [72,] 4.283796e+04 -1.018463e+05 [73,] 2.427217e+03 4.283796e+04 [74,] -1.322708e+05 2.427217e+03 [75,] -4.093145e+04 -1.322708e+05 [76,] -8.739013e+04 -4.093145e+04 [77,] -3.093746e+04 -8.739013e+04 [78,] -6.970571e+04 -3.093746e+04 [79,] 4.885414e+04 -6.970571e+04 [80,] -1.180068e+05 4.885414e+04 [81,] -5.111407e+04 -1.180068e+05 [82,] 1.831975e+04 -5.111407e+04 [83,] 5.431579e+04 1.831975e+04 [84,] -1.254997e+05 5.431579e+04 [85,] -1.506956e+04 -1.254997e+05 [86,] -1.013768e+04 -1.506956e+04 [87,] -1.270529e+05 -1.013768e+04 [88,] 2.083809e+04 -1.270529e+05 [89,] -1.259173e+04 2.083809e+04 [90,] -5.295845e+04 -1.259173e+04 [91,] 8.081311e+04 -5.295845e+04 [92,] -1.550252e+05 8.081311e+04 [93,] 6.448599e+04 -1.550252e+05 [94,] -5.258125e+04 6.448599e+04 [95,] 1.040454e+05 -5.258125e+04 [96,] -3.726888e+04 1.040454e+05 [97,] -3.086653e+04 -3.726888e+04 [98,] -3.412860e+04 -3.086653e+04 [99,] -7.991547e+04 -3.412860e+04 [100,] -8.335505e+04 -7.991547e+04 [101,] -1.278187e+05 -8.335505e+04 [102,] -1.436685e+04 -1.278187e+05 [103,] 5.683226e+04 -1.436685e+04 [104,] -3.521946e+04 5.683226e+04 [105,] -2.898622e+05 -3.521946e+04 [106,] 1.064864e+05 -2.898622e+05 [107,] -1.553016e+05 1.064864e+05 [108,] -1.252069e+05 -1.553016e+05 [109,] 1.170282e+05 -1.252069e+05 [110,] -2.213077e+05 1.170282e+05 [111,] -1.524851e+05 -2.213077e+05 [112,] -8.239323e+04 -1.524851e+05 [113,] 5.174149e+04 -8.239323e+04 [114,] -1.529056e+03 5.174149e+04 [115,] 3.358610e+03 -1.529056e+03 [116,] 7.191914e+04 3.358610e+03 [117,] -1.569444e+05 7.191914e+04 [118,] 3.564294e+04 -1.569444e+05 [119,] -2.009287e+05 3.564294e+04 [120,] 1.732224e+04 -2.009287e+05 [121,] -1.134618e+05 1.732224e+04 [122,] 5.424982e+04 -1.134618e+05 [123,] -1.295347e+04 5.424982e+04 [124,] -2.256727e+05 -1.295347e+04 [125,] 3.600235e+04 -2.256727e+05 [126,] 1.885287e+04 3.600235e+04 [127,] 2.591166e+03 1.885287e+04 [128,] 7.816804e+04 2.591166e+03 [129,] 2.415626e+04 7.816804e+04 [130,] 8.192810e+04 2.415626e+04 [131,] 7.666630e+04 8.192810e+04 [132,] 2.664539e+03 7.666630e+04 [133,] 1.855164e+04 2.664539e+03 [134,] 2.592867e+04 1.855164e+04 [135,] 6.491739e+04 2.592867e+04 [136,] 3.394453e+04 6.491739e+04 [137,] 1.820984e+04 3.394453e+04 [138,] 3.600235e+04 1.820984e+04 [139,] 3.155258e+04 3.600235e+04 [140,] 7.816804e+04 3.155258e+04 [141,] -4.587239e+03 7.816804e+04 [142,] 2.975914e+04 -4.587239e+03 [143,] 6.067965e+04 2.975914e+04 [144,] 2.176231e+04 6.067965e+04 [145,] 7.816804e+04 2.176231e+04 [146,] 2.482765e+04 7.816804e+04 [147,] 8.161402e+04 2.482765e+04 [148,] 7.808688e+04 8.161402e+04 [149,] -1.993497e+03 7.808688e+04 [150,] -6.439980e+04 -1.993497e+03 [151,] 2.849368e+04 -6.439980e+04 [152,] 4.710343e+04 2.849368e+04 [153,] 3.231582e+04 4.710343e+04 [154,] 3.600235e+04 3.231582e+04 [155,] 3.600235e+04 3.600235e+04 [156,] 3.057385e+04 3.600235e+04 [157,] 4.425873e+04 3.057385e+04 [158,] 3.600235e+04 4.425873e+04 [159,] 3.600235e+04 3.600235e+04 [160,] 7.816804e+04 3.600235e+04 [161,] 3.600235e+04 7.816804e+04 [162,] 1.013438e+04 3.600235e+04 [163,] 1.975804e+04 1.013438e+04 [164,] 3.340956e+04 1.975804e+04 [165,] 7.816804e+04 3.340956e+04 [166,] 1.293051e+05 7.816804e+04 [167,] -3.697539e+04 1.293051e+05 [168,] 3.600235e+04 -3.697539e+04 [169,] 3.600235e+04 3.600235e+04 [170,] 6.909576e+04 3.600235e+04 [171,] 7.816804e+04 6.909576e+04 [172,] 5.465479e+04 7.816804e+04 [173,] 7.816804e+04 5.465479e+04 [174,] 2.225224e+04 7.816804e+04 [175,] -1.408540e+05 2.225224e+04 [176,] 1.114887e+04 -1.408540e+05 [177,] 3.600235e+04 1.114887e+04 [178,] 3.600235e+04 3.600235e+04 [179,] 6.651030e+04 3.600235e+04 [180,] -5.422786e+04 6.651030e+04 [181,] -2.700112e+04 -5.422786e+04 [182,] -9.277563e+04 -2.700112e+04 [183,] 4.207503e+02 -9.277563e+04 [184,] 5.017889e+03 4.207503e+02 [185,] -6.335012e+04 5.017889e+03 [186,] 2.127199e+04 -6.335012e+04 [187,] 3.739862e+04 2.127199e+04 [188,] 3.600235e+04 3.739862e+04 [189,] 2.704002e+04 3.600235e+04 [190,] 3.600235e+04 2.704002e+04 [191,] 2.246269e+04 3.600235e+04 [192,] 7.816804e+04 2.246269e+04 [193,] 3.122124e+04 7.816804e+04 [194,] 1.128268e+04 3.122124e+04 [195,] -1.625690e+05 1.128268e+04 [196,] 1.128268e+04 -1.625690e+05 [197,] 1.128268e+04 1.128268e+04 [198,] -9.521193e+04 1.128268e+04 [199,] 1.128268e+04 -9.521193e+04 [200,] 1.128268e+04 1.128268e+04 [201,] 1.128268e+04 1.128268e+04 [202,] -5.243307e+04 1.128268e+04 [203,] -3.907353e+04 -5.243307e+04 [204,] 3.153210e+03 -3.907353e+04 [205,] 3.939561e+03 3.153210e+03 [206,] 1.128268e+04 3.939561e+03 [207,] -3.948385e+04 1.128268e+04 [208,] -4.716405e+04 -3.948385e+04 [209,] 3.600235e+04 -4.716405e+04 [210,] 7.816804e+04 3.600235e+04 [211,] 7.713665e+04 7.816804e+04 [212,] 2.906961e+04 7.713665e+04 [213,] 5.044902e+04 2.906961e+04 [214,] 7.516457e+04 5.044902e+04 [215,] 2.696010e+04 7.516457e+04 [216,] 2.293713e+04 2.696010e+04 [217,] 3.299888e+04 2.293713e+04 [218,] -1.654379e+04 3.299888e+04 [219,] 5.158962e+04 -1.654379e+04 [220,] 7.816804e+04 5.158962e+04 [221,] 7.663229e+04 7.816804e+04 [222,] -3.156698e+03 7.663229e+04 [223,] 3.432896e+04 -3.156698e+03 [224,] -1.929492e+03 3.432896e+04 [225,] 2.088802e+04 -1.929492e+03 [226,] 7.816804e+04 2.088802e+04 [227,] 6.352101e+04 7.816804e+04 [228,] 3.600235e+04 6.352101e+04 [229,] 4.264214e+04 3.600235e+04 [230,] 6.133615e+04 4.264214e+04 [231,] 7.816804e+04 6.133615e+04 [232,] -1.624027e+03 7.816804e+04 [233,] 3.600235e+04 -1.624027e+03 [234,] 3.605292e+04 3.600235e+04 [235,] 3.600235e+04 3.605292e+04 [236,] 2.723941e+04 3.600235e+04 [237,] 7.816804e+04 2.723941e+04 [238,] 7.816804e+04 7.816804e+04 [239,] 3.600235e+04 7.816804e+04 [240,] 6.911561e+04 3.600235e+04 [241,] 3.600235e+04 6.911561e+04 [242,] 6.883354e+04 3.600235e+04 [243,] 3.852101e+04 6.883354e+04 [244,] 1.128268e+04 3.852101e+04 [245,] 1.128268e+04 1.128268e+04 [246,] 3.956546e+03 1.128268e+04 [247,] 1.128268e+04 3.956546e+03 [248,] -3.088301e+04 1.128268e+04 [249,] 5.369742e+04 -3.088301e+04 [250,] 1.128268e+04 5.369742e+04 [251,] 6.171967e+04 1.128268e+04 [252,] -5.658248e+04 6.171967e+04 [253,] -1.706214e+05 -5.658248e+04 [254,] 3.600235e+04 -1.706214e+05 [255,] 2.306742e+04 3.600235e+04 [256,] -6.896001e+03 2.306742e+04 [257,] 1.128268e+04 -6.896001e+03 [258,] 7.816804e+04 1.128268e+04 [259,] 1.128268e+04 7.816804e+04 [260,] 1.128268e+04 1.128268e+04 [261,] 2.574404e+04 1.128268e+04 [262,] 7.838025e+04 2.574404e+04 [263,] 1.128268e+04 7.838025e+04 [264,] 7.816804e+04 1.128268e+04 [265,] -3.088301e+04 7.816804e+04 [266,] 1.051092e+05 -3.088301e+04 [267,] 1.128268e+04 1.051092e+05 [268,] 4.108224e+04 1.128268e+04 [269,] 1.128268e+04 4.108224e+04 [270,] 1.128268e+04 1.128268e+04 [271,] 1.128268e+04 1.128268e+04 [272,] 6.072985e+02 1.128268e+04 [273,] 3.600235e+04 6.072985e+02 [274,] 5.887555e+04 3.600235e+04 [275,] -2.433108e+03 5.887555e+04 [276,] 7.816804e+04 -2.433108e+03 [277,] -8.162836e+04 7.816804e+04 [278,] 4.916413e+04 -8.162836e+04 [279,] 1.128268e+04 4.916413e+04 [280,] 3.600235e+04 1.128268e+04 [281,] 3.950952e+04 3.600235e+04 [282,] 3.600235e+04 3.950952e+04 [283,] 6.270590e+04 3.600235e+04 [284,] 7.816804e+04 6.270590e+04 [285,] 2.043001e+01 7.816804e+04 [286,] 3.203934e+04 2.043001e+01 [287,] -3.088301e+04 3.203934e+04 [288,] -6.111311e+02 -3.088301e+04 [289,] 1.128268e+04 -6.111311e+02 [290,] 2.281855e+04 1.128268e+04 [291,] 6.682024e+04 2.281855e+04 [292,] 3.600235e+04 6.682024e+04 [293,] 5.115761e+04 3.600235e+04 [294,] 2.037284e+04 5.115761e+04 [295,] 3.600235e+04 2.037284e+04 [296,] -3.088301e+04 3.600235e+04 [297,] 7.816804e+04 -3.088301e+04 [298,] 2.248674e+04 7.816804e+04 [299,] 2.575388e+03 2.248674e+04 [300,] 6.140890e+04 2.575388e+03 [301,] 7.666630e+04 6.140890e+04 [302,] 1.128268e+04 7.666630e+04 [303,] 3.792899e+04 1.128268e+04 [304,] 1.794995e+04 3.792899e+04 [305,] 1.128268e+04 1.794995e+04 [306,] -2.309127e+04 1.128268e+04 [307,] 8.272152e+03 -2.309127e+04 [308,] 1.128268e+04 8.272152e+03 [309,] 1.128268e+04 1.128268e+04 [310,] 2.083068e+04 1.128268e+04 [311,] 1.128268e+04 2.083068e+04 [312,] -5.046879e+04 1.128268e+04 [313,] 2.417100e+04 -5.046879e+04 [314,] 4.523385e+04 2.417100e+04 [315,] 3.617547e+04 4.523385e+04 [316,] -1.752119e+04 3.617547e+04 [317,] 2.847441e+04 -1.752119e+04 [318,] 3.600235e+04 2.847441e+04 [319,] 7.816804e+04 3.600235e+04 [320,] 7.816804e+04 7.816804e+04 [321,] 7.816804e+04 7.816804e+04 [322,] -3.088301e+04 7.816804e+04 [323,] 7.816804e+04 -3.088301e+04 [324,] 7.350799e+04 7.816804e+04 [325,] -3.088301e+04 7.350799e+04 [326,] 1.128268e+04 -3.088301e+04 [327,] 7.816804e+04 1.128268e+04 [328,] 7.816804e+04 7.816804e+04 [329,] 7.816804e+04 7.816804e+04 [330,] 3.600235e+04 7.816804e+04 [331,] 7.816804e+04 3.600235e+04 [332,] 6.976128e+04 7.816804e+04 [333,] 3.448651e+04 6.976128e+04 [334,] 2.602428e+04 3.448651e+04 [335,] 1.796446e+04 2.602428e+04 [336,] 3.600235e+04 1.796446e+04 [337,] 3.243381e+04 3.600235e+04 [338,] 6.522883e+04 3.243381e+04 [339,] -8.705010e+03 6.522883e+04 [340,] 3.600235e+04 -8.705010e+03 [341,] 3.600235e+04 3.600235e+04 [342,] 7.816804e+04 3.600235e+04 [343,] 1.519857e+03 7.816804e+04 [344,] 3.600235e+04 1.519857e+03 [345,] 2.046556e+04 3.600235e+04 [346,] 3.600235e+04 2.046556e+04 [347,] 7.816804e+04 3.600235e+04 [348,] 4.847882e+03 7.816804e+04 [349,] 7.591426e+04 4.847882e+03 [350,] 2.285412e+04 7.591426e+04 [351,] 6.758874e+04 2.285412e+04 [352,] 2.980277e+04 6.758874e+04 [353,] 3.600235e+04 2.980277e+04 [354,] 3.600235e+04 3.600235e+04 [355,] -1.050841e+04 3.600235e+04 [356,] 3.600235e+04 -1.050841e+04 [357,] 3.600235e+04 3.600235e+04 [358,] 7.816804e+04 3.600235e+04 [359,] 2.999808e+04 7.816804e+04 [360,] 7.816804e+04 2.999808e+04 [361,] 3.600235e+04 7.816804e+04 [362,] -6.124955e+04 3.600235e+04 [363,] 6.485885e+04 -6.124955e+04 [364,] 7.816804e+04 6.485885e+04 [365,] 7.816804e+04 7.816804e+04 [366,] 3.600235e+04 7.816804e+04 [367,] 7.816804e+04 3.600235e+04 [368,] 2.291243e+04 7.816804e+04 [369,] 7.816804e+04 2.291243e+04 [370,] 7.125759e+04 7.816804e+04 [371,] 2.228522e+04 7.125759e+04 [372,] 5.794993e+04 2.228522e+04 [373,] 7.816804e+04 5.794993e+04 [374,] 7.816804e+04 7.816804e+04 [375,] 4.627300e+04 7.816804e+04 [376,] 2.517615e+04 4.627300e+04 [377,] -2.442581e+03 2.517615e+04 [378,] 7.816804e+04 -2.442581e+03 [379,] 7.816804e+04 7.816804e+04 [380,] -1.304973e+04 7.816804e+04 [381,] -1.107369e+03 -1.304973e+04 [382,] 7.816804e+04 -1.107369e+03 [383,] 1.626159e+04 7.816804e+04 [384,] -8.365301e+04 1.626159e+04 [385,] 4.297161e+04 -8.365301e+04 [386,] 5.401746e+04 4.297161e+04 [387,] 2.985540e+04 5.401746e+04 [388,] 5.438925e+04 2.985540e+04 [389,] 6.238896e+04 5.438925e+04 [390,] -3.806724e+01 6.238896e+04 [391,] 5.889624e+03 -3.806724e+01 [392,] 6.055375e+04 5.889624e+03 [393,] -1.501482e+05 6.055375e+04 [394,] 3.034063e+04 -1.501482e+05 [395,] -8.392311e+04 3.034063e+04 [396,] -3.247143e+05 -8.392311e+04 [397,] 3.758911e+04 -3.247143e+05 [398,] -1.079982e+05 3.758911e+04 [399,] -6.169317e+03 -1.079982e+05 [400,] 2.756067e+04 -6.169317e+03 [401,] 2.896510e+03 2.756067e+04 [402,] -1.297722e+05 2.896510e+03 [403,] -2.527709e+04 -1.297722e+05 [404,] 2.666550e+03 -2.527709e+04 [405,] 4.559749e+04 2.666550e+03 [406,] 3.927075e+04 4.559749e+04 [407,] 2.931348e+03 3.927075e+04 [408,] -3.157747e+04 2.931348e+03 [409,] -4.223222e+04 -3.157747e+04 [410,] -1.457241e+05 -4.223222e+04 [411,] -6.694617e+04 -1.457241e+05 [412,] 3.464206e+04 -6.694617e+04 [413,] -1.379192e+05 3.464206e+04 [414,] -1.075542e+05 -1.379192e+05 [415,] 1.288701e+05 -1.075542e+05 [416,] -1.022662e+05 1.288701e+05 [417,] -6.764406e+04 -1.022662e+05 [418,] -4.973243e+04 -6.764406e+04 [419,] -5.246928e+04 -4.973243e+04 [420,] -1.595566e+05 -5.246928e+04 [421,] -5.328526e+05 -1.595566e+05 [422,] -3.851035e+05 -5.328526e+05 [423,] -1.557869e+04 -3.851035e+05 [424,] -2.849142e+05 -1.557869e+04 [425,] -2.636122e+05 -2.849142e+05 [426,] -8.778464e+05 -2.636122e+05 [427,] -1.739617e+05 -8.778464e+05 [428,] -2.086613e+05 -1.739617e+05 [429,] -5.306982e+05 -2.086613e+05 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 1.958932e+06 2.945709e+06 2 -1.608369e+06 1.958932e+06 3 -6.533468e+04 -1.608369e+06 4 3.622083e+05 -6.533468e+04 5 1.569258e+05 3.622083e+05 6 3.729608e+05 1.569258e+05 7 -1.382227e+05 3.729608e+05 8 7.579032e+05 -1.382227e+05 9 1.030653e+05 7.579032e+05 10 4.217010e+05 1.030653e+05 11 3.683217e+05 4.217010e+05 12 -8.069330e+05 3.683217e+05 13 -1.555096e+05 -8.069330e+05 14 4.914049e+05 -1.555096e+05 15 -1.921223e+05 4.914049e+05 16 -1.947991e+05 -1.921223e+05 17 -1.065107e+05 -1.947991e+05 18 -5.178519e+05 -1.065107e+05 19 2.502813e+05 -5.178519e+05 20 -4.255778e+04 2.502813e+05 21 -4.201810e+04 -4.255778e+04 22 -5.939069e+05 -4.201810e+04 23 -2.534845e+05 -5.939069e+05 24 6.050526e+04 -2.534845e+05 25 -3.689886e+05 6.050526e+04 26 -2.475860e+05 -3.689886e+05 27 4.290343e+05 -2.475860e+05 28 -1.733252e+05 4.290343e+05 29 2.672605e+05 -1.733252e+05 30 -4.558695e+04 2.672605e+05 31 2.200257e+05 -4.558695e+04 32 1.174329e+05 2.200257e+05 33 3.411520e+04 1.174329e+05 34 -4.667818e+05 3.411520e+04 35 1.238381e+05 -4.667818e+05 36 1.356866e+05 1.238381e+05 37 -8.574527e+04 1.356866e+05 38 -3.037759e+05 -8.574527e+04 39 -4.832529e+05 -3.037759e+05 40 -1.970856e+05 -4.832529e+05 41 -1.326601e+04 -1.970856e+05 42 -1.952668e+05 -1.326601e+04 43 -2.031444e+04 -1.952668e+05 44 -6.659535e+04 -2.031444e+04 45 -2.718181e+05 -6.659535e+04 46 -2.949503e+05 -2.718181e+05 47 -6.060906e+01 -2.949503e+05 48 -1.543829e+05 -6.060906e+01 49 -1.271625e+04 -1.543829e+05 50 -1.427621e+05 -1.271625e+04 51 -2.242160e+05 -1.427621e+05 52 -4.412895e+04 -2.242160e+05 53 -4.303624e+05 -4.412895e+04 54 -1.997302e+04 -4.303624e+05 55 3.734755e+05 -1.997302e+04 56 -8.144209e+03 3.734755e+05 57 -4.126198e+04 -8.144209e+03 58 -1.979514e+05 -4.126198e+04 59 -3.327342e+04 -1.979514e+05 60 -1.996854e+05 -3.327342e+04 61 -1.262337e+04 -1.996854e+05 62 -2.224404e+05 -1.262337e+04 63 -3.673721e+04 -2.224404e+05 64 -4.644577e+05 -3.673721e+04 65 -3.685077e+04 -4.644577e+05 66 -5.278891e+04 -3.685077e+04 67 -8.535061e+04 -5.278891e+04 68 -1.114192e+05 -8.535061e+04 69 -1.687591e+05 -1.114192e+05 70 6.897109e+03 -1.687591e+05 71 -1.018463e+05 6.897109e+03 72 4.283796e+04 -1.018463e+05 73 2.427217e+03 4.283796e+04 74 -1.322708e+05 2.427217e+03 75 -4.093145e+04 -1.322708e+05 76 -8.739013e+04 -4.093145e+04 77 -3.093746e+04 -8.739013e+04 78 -6.970571e+04 -3.093746e+04 79 4.885414e+04 -6.970571e+04 80 -1.180068e+05 4.885414e+04 81 -5.111407e+04 -1.180068e+05 82 1.831975e+04 -5.111407e+04 83 5.431579e+04 1.831975e+04 84 -1.254997e+05 5.431579e+04 85 -1.506956e+04 -1.254997e+05 86 -1.013768e+04 -1.506956e+04 87 -1.270529e+05 -1.013768e+04 88 2.083809e+04 -1.270529e+05 89 -1.259173e+04 2.083809e+04 90 -5.295845e+04 -1.259173e+04 91 8.081311e+04 -5.295845e+04 92 -1.550252e+05 8.081311e+04 93 6.448599e+04 -1.550252e+05 94 -5.258125e+04 6.448599e+04 95 1.040454e+05 -5.258125e+04 96 -3.726888e+04 1.040454e+05 97 -3.086653e+04 -3.726888e+04 98 -3.412860e+04 -3.086653e+04 99 -7.991547e+04 -3.412860e+04 100 -8.335505e+04 -7.991547e+04 101 -1.278187e+05 -8.335505e+04 102 -1.436685e+04 -1.278187e+05 103 5.683226e+04 -1.436685e+04 104 -3.521946e+04 5.683226e+04 105 -2.898622e+05 -3.521946e+04 106 1.064864e+05 -2.898622e+05 107 -1.553016e+05 1.064864e+05 108 -1.252069e+05 -1.553016e+05 109 1.170282e+05 -1.252069e+05 110 -2.213077e+05 1.170282e+05 111 -1.524851e+05 -2.213077e+05 112 -8.239323e+04 -1.524851e+05 113 5.174149e+04 -8.239323e+04 114 -1.529056e+03 5.174149e+04 115 3.358610e+03 -1.529056e+03 116 7.191914e+04 3.358610e+03 117 -1.569444e+05 7.191914e+04 118 3.564294e+04 -1.569444e+05 119 -2.009287e+05 3.564294e+04 120 1.732224e+04 -2.009287e+05 121 -1.134618e+05 1.732224e+04 122 5.424982e+04 -1.134618e+05 123 -1.295347e+04 5.424982e+04 124 -2.256727e+05 -1.295347e+04 125 3.600235e+04 -2.256727e+05 126 1.885287e+04 3.600235e+04 127 2.591166e+03 1.885287e+04 128 7.816804e+04 2.591166e+03 129 2.415626e+04 7.816804e+04 130 8.192810e+04 2.415626e+04 131 7.666630e+04 8.192810e+04 132 2.664539e+03 7.666630e+04 133 1.855164e+04 2.664539e+03 134 2.592867e+04 1.855164e+04 135 6.491739e+04 2.592867e+04 136 3.394453e+04 6.491739e+04 137 1.820984e+04 3.394453e+04 138 3.600235e+04 1.820984e+04 139 3.155258e+04 3.600235e+04 140 7.816804e+04 3.155258e+04 141 -4.587239e+03 7.816804e+04 142 2.975914e+04 -4.587239e+03 143 6.067965e+04 2.975914e+04 144 2.176231e+04 6.067965e+04 145 7.816804e+04 2.176231e+04 146 2.482765e+04 7.816804e+04 147 8.161402e+04 2.482765e+04 148 7.808688e+04 8.161402e+04 149 -1.993497e+03 7.808688e+04 150 -6.439980e+04 -1.993497e+03 151 2.849368e+04 -6.439980e+04 152 4.710343e+04 2.849368e+04 153 3.231582e+04 4.710343e+04 154 3.600235e+04 3.231582e+04 155 3.600235e+04 3.600235e+04 156 3.057385e+04 3.600235e+04 157 4.425873e+04 3.057385e+04 158 3.600235e+04 4.425873e+04 159 3.600235e+04 3.600235e+04 160 7.816804e+04 3.600235e+04 161 3.600235e+04 7.816804e+04 162 1.013438e+04 3.600235e+04 163 1.975804e+04 1.013438e+04 164 3.340956e+04 1.975804e+04 165 7.816804e+04 3.340956e+04 166 1.293051e+05 7.816804e+04 167 -3.697539e+04 1.293051e+05 168 3.600235e+04 -3.697539e+04 169 3.600235e+04 3.600235e+04 170 6.909576e+04 3.600235e+04 171 7.816804e+04 6.909576e+04 172 5.465479e+04 7.816804e+04 173 7.816804e+04 5.465479e+04 174 2.225224e+04 7.816804e+04 175 -1.408540e+05 2.225224e+04 176 1.114887e+04 -1.408540e+05 177 3.600235e+04 1.114887e+04 178 3.600235e+04 3.600235e+04 179 6.651030e+04 3.600235e+04 180 -5.422786e+04 6.651030e+04 181 -2.700112e+04 -5.422786e+04 182 -9.277563e+04 -2.700112e+04 183 4.207503e+02 -9.277563e+04 184 5.017889e+03 4.207503e+02 185 -6.335012e+04 5.017889e+03 186 2.127199e+04 -6.335012e+04 187 3.739862e+04 2.127199e+04 188 3.600235e+04 3.739862e+04 189 2.704002e+04 3.600235e+04 190 3.600235e+04 2.704002e+04 191 2.246269e+04 3.600235e+04 192 7.816804e+04 2.246269e+04 193 3.122124e+04 7.816804e+04 194 1.128268e+04 3.122124e+04 195 -1.625690e+05 1.128268e+04 196 1.128268e+04 -1.625690e+05 197 1.128268e+04 1.128268e+04 198 -9.521193e+04 1.128268e+04 199 1.128268e+04 -9.521193e+04 200 1.128268e+04 1.128268e+04 201 1.128268e+04 1.128268e+04 202 -5.243307e+04 1.128268e+04 203 -3.907353e+04 -5.243307e+04 204 3.153210e+03 -3.907353e+04 205 3.939561e+03 3.153210e+03 206 1.128268e+04 3.939561e+03 207 -3.948385e+04 1.128268e+04 208 -4.716405e+04 -3.948385e+04 209 3.600235e+04 -4.716405e+04 210 7.816804e+04 3.600235e+04 211 7.713665e+04 7.816804e+04 212 2.906961e+04 7.713665e+04 213 5.044902e+04 2.906961e+04 214 7.516457e+04 5.044902e+04 215 2.696010e+04 7.516457e+04 216 2.293713e+04 2.696010e+04 217 3.299888e+04 2.293713e+04 218 -1.654379e+04 3.299888e+04 219 5.158962e+04 -1.654379e+04 220 7.816804e+04 5.158962e+04 221 7.663229e+04 7.816804e+04 222 -3.156698e+03 7.663229e+04 223 3.432896e+04 -3.156698e+03 224 -1.929492e+03 3.432896e+04 225 2.088802e+04 -1.929492e+03 226 7.816804e+04 2.088802e+04 227 6.352101e+04 7.816804e+04 228 3.600235e+04 6.352101e+04 229 4.264214e+04 3.600235e+04 230 6.133615e+04 4.264214e+04 231 7.816804e+04 6.133615e+04 232 -1.624027e+03 7.816804e+04 233 3.600235e+04 -1.624027e+03 234 3.605292e+04 3.600235e+04 235 3.600235e+04 3.605292e+04 236 2.723941e+04 3.600235e+04 237 7.816804e+04 2.723941e+04 238 7.816804e+04 7.816804e+04 239 3.600235e+04 7.816804e+04 240 6.911561e+04 3.600235e+04 241 3.600235e+04 6.911561e+04 242 6.883354e+04 3.600235e+04 243 3.852101e+04 6.883354e+04 244 1.128268e+04 3.852101e+04 245 1.128268e+04 1.128268e+04 246 3.956546e+03 1.128268e+04 247 1.128268e+04 3.956546e+03 248 -3.088301e+04 1.128268e+04 249 5.369742e+04 -3.088301e+04 250 1.128268e+04 5.369742e+04 251 6.171967e+04 1.128268e+04 252 -5.658248e+04 6.171967e+04 253 -1.706214e+05 -5.658248e+04 254 3.600235e+04 -1.706214e+05 255 2.306742e+04 3.600235e+04 256 -6.896001e+03 2.306742e+04 257 1.128268e+04 -6.896001e+03 258 7.816804e+04 1.128268e+04 259 1.128268e+04 7.816804e+04 260 1.128268e+04 1.128268e+04 261 2.574404e+04 1.128268e+04 262 7.838025e+04 2.574404e+04 263 1.128268e+04 7.838025e+04 264 7.816804e+04 1.128268e+04 265 -3.088301e+04 7.816804e+04 266 1.051092e+05 -3.088301e+04 267 1.128268e+04 1.051092e+05 268 4.108224e+04 1.128268e+04 269 1.128268e+04 4.108224e+04 270 1.128268e+04 1.128268e+04 271 1.128268e+04 1.128268e+04 272 6.072985e+02 1.128268e+04 273 3.600235e+04 6.072985e+02 274 5.887555e+04 3.600235e+04 275 -2.433108e+03 5.887555e+04 276 7.816804e+04 -2.433108e+03 277 -8.162836e+04 7.816804e+04 278 4.916413e+04 -8.162836e+04 279 1.128268e+04 4.916413e+04 280 3.600235e+04 1.128268e+04 281 3.950952e+04 3.600235e+04 282 3.600235e+04 3.950952e+04 283 6.270590e+04 3.600235e+04 284 7.816804e+04 6.270590e+04 285 2.043001e+01 7.816804e+04 286 3.203934e+04 2.043001e+01 287 -3.088301e+04 3.203934e+04 288 -6.111311e+02 -3.088301e+04 289 1.128268e+04 -6.111311e+02 290 2.281855e+04 1.128268e+04 291 6.682024e+04 2.281855e+04 292 3.600235e+04 6.682024e+04 293 5.115761e+04 3.600235e+04 294 2.037284e+04 5.115761e+04 295 3.600235e+04 2.037284e+04 296 -3.088301e+04 3.600235e+04 297 7.816804e+04 -3.088301e+04 298 2.248674e+04 7.816804e+04 299 2.575388e+03 2.248674e+04 300 6.140890e+04 2.575388e+03 301 7.666630e+04 6.140890e+04 302 1.128268e+04 7.666630e+04 303 3.792899e+04 1.128268e+04 304 1.794995e+04 3.792899e+04 305 1.128268e+04 1.794995e+04 306 -2.309127e+04 1.128268e+04 307 8.272152e+03 -2.309127e+04 308 1.128268e+04 8.272152e+03 309 1.128268e+04 1.128268e+04 310 2.083068e+04 1.128268e+04 311 1.128268e+04 2.083068e+04 312 -5.046879e+04 1.128268e+04 313 2.417100e+04 -5.046879e+04 314 4.523385e+04 2.417100e+04 315 3.617547e+04 4.523385e+04 316 -1.752119e+04 3.617547e+04 317 2.847441e+04 -1.752119e+04 318 3.600235e+04 2.847441e+04 319 7.816804e+04 3.600235e+04 320 7.816804e+04 7.816804e+04 321 7.816804e+04 7.816804e+04 322 -3.088301e+04 7.816804e+04 323 7.816804e+04 -3.088301e+04 324 7.350799e+04 7.816804e+04 325 -3.088301e+04 7.350799e+04 326 1.128268e+04 -3.088301e+04 327 7.816804e+04 1.128268e+04 328 7.816804e+04 7.816804e+04 329 7.816804e+04 7.816804e+04 330 3.600235e+04 7.816804e+04 331 7.816804e+04 3.600235e+04 332 6.976128e+04 7.816804e+04 333 3.448651e+04 6.976128e+04 334 2.602428e+04 3.448651e+04 335 1.796446e+04 2.602428e+04 336 3.600235e+04 1.796446e+04 337 3.243381e+04 3.600235e+04 338 6.522883e+04 3.243381e+04 339 -8.705010e+03 6.522883e+04 340 3.600235e+04 -8.705010e+03 341 3.600235e+04 3.600235e+04 342 7.816804e+04 3.600235e+04 343 1.519857e+03 7.816804e+04 344 3.600235e+04 1.519857e+03 345 2.046556e+04 3.600235e+04 346 3.600235e+04 2.046556e+04 347 7.816804e+04 3.600235e+04 348 4.847882e+03 7.816804e+04 349 7.591426e+04 4.847882e+03 350 2.285412e+04 7.591426e+04 351 6.758874e+04 2.285412e+04 352 2.980277e+04 6.758874e+04 353 3.600235e+04 2.980277e+04 354 3.600235e+04 3.600235e+04 355 -1.050841e+04 3.600235e+04 356 3.600235e+04 -1.050841e+04 357 3.600235e+04 3.600235e+04 358 7.816804e+04 3.600235e+04 359 2.999808e+04 7.816804e+04 360 7.816804e+04 2.999808e+04 361 3.600235e+04 7.816804e+04 362 -6.124955e+04 3.600235e+04 363 6.485885e+04 -6.124955e+04 364 7.816804e+04 6.485885e+04 365 7.816804e+04 7.816804e+04 366 3.600235e+04 7.816804e+04 367 7.816804e+04 3.600235e+04 368 2.291243e+04 7.816804e+04 369 7.816804e+04 2.291243e+04 370 7.125759e+04 7.816804e+04 371 2.228522e+04 7.125759e+04 372 5.794993e+04 2.228522e+04 373 7.816804e+04 5.794993e+04 374 7.816804e+04 7.816804e+04 375 4.627300e+04 7.816804e+04 376 2.517615e+04 4.627300e+04 377 -2.442581e+03 2.517615e+04 378 7.816804e+04 -2.442581e+03 379 7.816804e+04 7.816804e+04 380 -1.304973e+04 7.816804e+04 381 -1.107369e+03 -1.304973e+04 382 7.816804e+04 -1.107369e+03 383 1.626159e+04 7.816804e+04 384 -8.365301e+04 1.626159e+04 385 4.297161e+04 -8.365301e+04 386 5.401746e+04 4.297161e+04 387 2.985540e+04 5.401746e+04 388 5.438925e+04 2.985540e+04 389 6.238896e+04 5.438925e+04 390 -3.806724e+01 6.238896e+04 391 5.889624e+03 -3.806724e+01 392 6.055375e+04 5.889624e+03 393 -1.501482e+05 6.055375e+04 394 3.034063e+04 -1.501482e+05 395 -8.392311e+04 3.034063e+04 396 -3.247143e+05 -8.392311e+04 397 3.758911e+04 -3.247143e+05 398 -1.079982e+05 3.758911e+04 399 -6.169317e+03 -1.079982e+05 400 2.756067e+04 -6.169317e+03 401 2.896510e+03 2.756067e+04 402 -1.297722e+05 2.896510e+03 403 -2.527709e+04 -1.297722e+05 404 2.666550e+03 -2.527709e+04 405 4.559749e+04 2.666550e+03 406 3.927075e+04 4.559749e+04 407 2.931348e+03 3.927075e+04 408 -3.157747e+04 2.931348e+03 409 -4.223222e+04 -3.157747e+04 410 -1.457241e+05 -4.223222e+04 411 -6.694617e+04 -1.457241e+05 412 3.464206e+04 -6.694617e+04 413 -1.379192e+05 3.464206e+04 414 -1.075542e+05 -1.379192e+05 415 1.288701e+05 -1.075542e+05 416 -1.022662e+05 1.288701e+05 417 -6.764406e+04 -1.022662e+05 418 -4.973243e+04 -6.764406e+04 419 -5.246928e+04 -4.973243e+04 420 -1.595566e+05 -5.246928e+04 421 -5.328526e+05 -1.595566e+05 422 -3.851035e+05 -5.328526e+05 423 -1.557869e+04 -3.851035e+05 424 -2.849142e+05 -1.557869e+04 425 -2.636122e+05 -2.849142e+05 426 -8.778464e+05 -2.636122e+05 427 -1.739617e+05 -8.778464e+05 428 -2.086613e+05 -1.739617e+05 429 -5.306982e+05 -2.086613e+05 > 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/7p2wq1292932180.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/8zbds1292932180.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/9zbds1292932180.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/10zbds1292932180.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/11lccg1292932180.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/12oca41292932180.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/13dvpg1292932180.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/1464o11292932180.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/15kema1292932180.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/16ne3x1292932180.tab") + } > > try(system("convert tmp/13jx21292932180.ps tmp/13jx21292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/23jx21292932180.ps tmp/23jx21292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/33jx21292932180.ps tmp/33jx21292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/4esxn1292932180.ps tmp/4esxn1292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/5esxn1292932180.ps tmp/5esxn1292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/6esxn1292932180.ps tmp/6esxn1292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/7p2wq1292932180.ps tmp/7p2wq1292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/8zbds1292932180.ps tmp/8zbds1292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/9zbds1292932180.ps tmp/9zbds1292932180.png",intern=TRUE)) character(0) > try(system("convert tmp/10zbds1292932180.ps tmp/10zbds1292932180.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.987 2.234 31.619