R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1 + ,162556 + ,1081 + ,807 + ,213118 + ,230380558 + ,6282929 + ,1 + ,29790 + ,309 + ,444 + ,81767 + ,25266003 + ,4324047 + ,1 + ,87550 + ,458 + ,412 + ,153198 + ,70164684 + ,4108272 + ,0 + ,84738 + ,588 + ,428 + ,-26007 + ,-15292116 + ,-1212617 + ,1 + ,54660 + ,299 + ,312 + ,126942 + ,37955658 + ,1485329 + ,1 + ,42634 + ,156 + ,166 + ,157214 + ,24525384 + ,1779876 + ,0 + ,40949 + ,481 + ,263 + ,129352 + ,62218312 + ,1367203 + ,1 + ,42312 + ,323 + ,237 + ,234817 + ,75845891 + ,2519076 + ,1 + ,37704 + ,452 + ,228 + ,60448 + ,27322496 + ,912684 + ,1 + ,16275 + ,109 + ,129 + ,47818 + ,5212162 + ,1443586 + ,0 + ,25830 + ,115 + ,104 + ,245546 + ,28237790 + ,1220017 + ,0 + ,12679 + ,110 + ,122 + ,48020 + ,5282200 + ,984885 + ,1 + ,18014 + ,239 + ,393 + ,-1710 + ,-408690 + ,1457425 + ,0 + ,43556 + ,247 + ,190 + ,32648 + ,8064056 + ,-572920 + ,1 + ,24524 + ,497 + ,275 + ,95350 + ,47388950 + ,929144 + ,0 + ,6532 + ,103 + ,62 + ,151352 + ,15589256 + ,1151176 + ,0 + ,7123 + ,109 + ,102 + ,288170 + ,31410530 + ,790090 + ,1 + ,20813 + ,502 + ,255 + ,114337 + ,57397174 + ,774497 + ,1 + ,37597 + ,248 + ,234 + ,37884 + ,9395232 + ,990576 + ,0 + ,17821 + ,373 + ,277 + ,122844 + ,45820812 + ,454195 + ,1 + ,12988 + ,119 + ,73 + ,82340 + ,9798460 + ,876607 + ,1 + ,22330 + ,84 + ,67 + ,79801 + ,6703284 + ,711969 + ,0 + ,13326 + ,102 + ,103 + ,165548 + ,16885896 + ,702380 + ,0 + ,16189 + ,295 + ,290 + ,116384 + ,34333280 + ,264449 + ,0 + ,7146 + ,105 + ,83 + ,134028 + ,14072940 + ,450033 + ,0 + ,15824 + ,64 + ,56 + ,63838 + ,4085632 + ,541063 + ,1 + ,26088 + ,267 + ,224 + ,74996 + ,20023932 + ,588864 + ,0 + ,11326 + ,129 + ,64 + ,31080 + ,4009320 + ,-37216 + ,0 + ,8568 + ,37 + ,34 + ,32168 + ,1190216 + ,783310 + ,0 + ,14416 + ,361 + ,139 + ,49857 + ,17998377 + ,467359 + ,1 + ,3369 + ,28 + ,26 + ,87161 + ,2440508 + ,688779 + ,1 + ,11819 + ,85 + ,70 + ,106113 + ,9019605 + ,608419 + ,1 + ,6620 + ,44 + ,40 + ,80570 + ,3545080 + ,696348 + ,1 + ,4519 + ,49 + ,42 + ,102129 + ,5004321 + ,597793 + ,0 + ,2220 + ,22 + ,12 + ,301670 + ,6636740 + ,821730 + ,0 + ,18562 + ,155 + ,211 + ,102313 + ,15858515 + ,377934 + ,0 + ,10327 + ,91 + ,74 + ,88577 + ,8060507 + ,651939 + ,1 + ,5336 + ,81 + ,80 + ,112477 + ,9110637 + ,697458 + ,1 + ,2365 + ,79 + ,83 + ,191778 + ,15150462 + ,700368 + ,0 + ,4069 + ,145 + ,131 + ,79804 + ,11571580 + ,225986 + ,0 + ,7710 + ,816 + ,187 + ,128294 + ,104687904 + ,348695 + ,0 + ,13718 + ,61 + ,56 + ,96448 + ,5883328 + ,373683 + ,0 + ,4525 + ,226 + ,89 + ,93811 + ,21201286 + ,501709 + ,0 + ,6869 + ,105 + ,88 + ,117520 + ,12339600 + ,413743 + ,0 + ,4628 + ,62 + ,39 + ,69159 + ,4287858 + ,379825 + ,1 + ,3653 + ,24 + ,24 + ,101792 + ,2443008 + ,336260 + ,1 + ,1265 + ,26 + ,19 + ,210568 + ,5474768 + ,636765 + ,1 + ,7489 + ,322 + ,149 + ,136996 + ,44112712 + ,481231 + ,0 + ,4901 + ,84 + ,58 + ,121920 + ,10241280 + ,469107 + ,0 + ,2284 + ,33 + ,41 + ,76403 + ,2521299 + ,211928 + ,1 + ,3160 + ,108 + ,90 + ,108094 + ,11674152 + ,563925 + ,1 + ,4150 + ,150 + ,136 + ,134759 + ,20213850 + ,511939 + ,1 + ,7285 + ,115 + ,97 + ,188873 + ,21720395 + ,521016 + ,1 + ,1134 + ,162 + ,63 + ,146216 + ,23686992 + ,543856 + ,1 + ,4658 + ,158 + ,114 + ,156608 + ,24744064 + ,329304 + ,0 + ,2384 + ,97 + ,77 + ,61348 + ,5950756 + ,423262 + ,0 + ,3748 + ,9 + ,6 + ,50350 + ,453150 + ,509665 + ,0 + ,5371 + ,66 + ,47 + ,87720 + ,5789520 + ,455881 + ,0 + ,1285 + ,107 + ,51 + ,99489 + ,10645323 + ,367772 + ,1 + ,9327 + ,101 + ,85 + ,87419 + ,8829319 + ,406339 + ,1 + ,5565 + ,47 + ,43 + ,94355 + ,4434685 + ,493408 + ,0 + ,1528 + ,38 + ,32 + ,60326 + ,2292388 + ,232942 + ,1 + ,3122 + ,34 + ,25 + ,94670 + ,3218780 + ,416002 + ,1 + ,7317 + ,84 + ,74 + ,82425 + ,6923700 + ,337430 + ,0 + ,2675 + ,79 + ,54 + ,59017 + ,4662343 + ,361517 + ,0 + ,13253 + ,947 + ,251 + ,90829 + ,86015063 + ,360962 + ,0 + ,880 + ,74 + ,15 + ,80791 + ,5978534 + ,235561 + ,1 + ,2053 + ,53 + ,44 + ,100423 + ,5322419 + ,408247 + ,0 + ,1424 + ,94 + ,73 + ,131116 + ,12324904 + ,450296 + ,1 + ,4036 + ,63 + ,85 + ,100269 + ,6316947 + ,418799 + ,1 + ,3045 + ,58 + ,49 + ,27330 + ,1585140 + ,247405 + ,0 + ,5119 + ,49 + ,38 + ,39039 + ,1912911 + ,378519 + ,0 + ,1431 + ,34 + ,35 + ,106885 + ,3634090 + ,326638 + ,0 + ,554 + ,11 + ,9 + ,79285 + ,872135 + ,328233 + ,0 + ,1975 + ,35 + ,34 + ,118881 + ,4160835 + ,386225 + ,1 + ,1286 + ,17 + ,17 + ,77623 + ,1319591 + ,283662 + ,0 + ,1012 + ,47 + ,29 + ,114768 + ,5394096 + ,370225 + ,0 + ,810 + ,43 + ,11 + ,74015 + ,3182645 + ,269236 + ,0 + ,1280 + ,117 + ,52 + ,69465 + ,8127405 + ,365732 + ,1 + ,666 + ,171 + ,13 + ,117869 + ,20155599 + ,420383 + ,0 + ,1380 + ,26 + ,29 + ,60982 + ,1585532 + ,345811 + ,1 + ,4608 + ,73 + ,64 + ,90131 + ,6579563 + ,431809 + ,0 + ,876 + ,59 + ,33 + ,138971 + ,8199289 + ,418876 + ,0 + ,814 + ,18 + ,15 + ,39625 + ,713250 + ,297476 + ,0 + ,514 + ,15 + ,15 + ,102725 + ,1540875 + ,416776 + ,1 + ,5692 + ,72 + ,68 + ,64239 + ,4625208 + ,357257 + ,0 + ,3642 + ,86 + ,100 + ,90262 + ,7762532 + ,458343 + ,0 + ,540 + ,14 + ,13 + ,103960 + ,1455440 + ,388386 + ,0 + ,2099 + ,64 + ,45 + ,106611 + ,6823104 + ,358934 + ,0 + ,567 + ,11 + ,14 + ,103345 + ,1136795 + ,407560 + ,0 + ,2001 + ,52 + ,36 + ,95551 + ,4968652 + ,392558 + ,1 + ,2949 + ,41 + ,40 + ,82903 + ,3399023 + ,373177 + ,0 + ,2253 + ,99 + ,68 + ,63593 + ,6295707 + ,428370 + ,1 + ,6533 + ,75 + ,29 + ,126910 + ,9518250 + ,369419 + ,0 + ,1889 + ,45 + ,43 + ,37527 + ,1688715 + ,358649 + ,1 + ,3055 + ,43 + ,30 + ,60247 + ,2590621 + ,376641 + ,0 + ,272 + ,8 + ,9 + ,112995 + ,903960 + ,467427 + ,1 + ,1414 + ,198 + ,22 + ,70184 + ,13896432 + ,364885 + ,0 + ,2564 + ,22 + ,19 + ,130140 + ,2863080 + ,436230 + ,1 + ,1383 + ,11 + ,9 + ,73221 + ,805431 + ,329118 + ,1 + ,1261 + ,33 + ,31 + ,76114 + ,2511762 + ,317365 + ,0 + ,975 + ,23 + ,19 + ,90534 + ,2082282 + ,286849 + ,0 + ,3366 + ,80 + ,55 + ,108479 + ,8678320 + ,376685 + ,0 + ,576 + ,18 + ,8 + ,113761 + ,2047698 + ,407198 + ,0 + ,1306 + ,28 + ,26 + ,68696 + ,1923488 + ,377772 + ,0 + ,746 + ,23 + ,29 + ,71561 + ,1645903 + ,271483 + ,1 + ,3192 + ,60 + ,48 + ,59831 + ,3589860 + ,153661 + ,1 + ,2045 + ,20 + ,16 + ,97890 + ,1957800 + ,513294 + ,0 + ,5477 + ,59 + ,45 + ,101481 + ,5987379 + ,324881 + ,1 + ,1932 + ,36 + ,20 + ,72954 + ,2626344 + ,264512 + ,0 + ,936 + ,30 + ,22 + ,67939 + ,2038170 + ,420968 + ,1 + ,3437 + ,47 + ,33 + ,48022 + ,2257034 + ,129302 + ,0 + ,5131 + ,71 + ,44 + ,86111 + ,6113881 + ,191521 + ,1 + ,2397 + ,14 + ,13 + ,74020 + ,1036280 + ,268673 + ,1 + ,1389 + ,9 + ,6 + ,57530 + ,517770 + ,353179 + ,0 + ,1503 + ,39 + ,35 + ,56364 + ,2198196 + ,354624 + ,0 + ,402 + ,26 + ,8 + ,84990 + ,2209740 + ,363713 + ,0 + ,2239 + ,21 + ,17 + ,88590 + ,1860390 + ,456657 + ,1 + ,2234 + ,16 + ,11 + ,77200 + ,1235200 + ,211742 + ,0 + ,837 + ,69 + ,21 + ,61262 + ,4227078 + ,338381 + ,0 + ,10579 + ,92 + ,92 + ,110309 + ,10148428 + ,418530 + ,0 + ,875 + ,14 + ,12 + ,67000 + ,938000 + ,351483 + ,0 + ,1395 + ,103 + ,108 + ,93099 + ,9589197 + ,372928 + ,1 + ,1659 + ,29 + ,25 + ,107577 + ,3119733 + ,485538 + ,1 + ,2647 + ,37 + ,17 + ,62920 + ,2328040 + ,279268 + ,1 + ,3294 + ,23 + ,23 + ,75832 + ,1744136 + ,219060 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,94 + ,7 + ,10 + ,60793 + ,425551 + ,325314 + ,0 + ,422 + ,28 + ,23 + ,57935 + ,1622180 + ,322046 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,34 + ,8 + ,7 + ,60630 + ,485040 + ,325599 + ,0 + ,1558 + ,63 + ,25 + ,55637 + ,3505131 + ,377028 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,43 + ,3 + ,20 + ,60887 + ,182661 + ,323850 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,316 + ,9 + ,4 + ,60505 + ,544545 + ,331514 + ,0 + ,115 + ,13 + ,10 + ,60945 + ,792285 + ,325632 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,389 + ,14 + ,7 + ,58990 + ,825860 + ,322265 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,1002 + ,15 + ,11 + ,56750 + ,851250 + ,325906 + ,0 + ,36 + ,3 + ,4 + ,60894 + ,182682 + ,325985 + ,0 + ,460 + ,15 + ,15 + ,63346 + ,950190 + ,346145 + ,0 + ,309 + ,11 + ,9 + ,56535 + ,621885 + ,325898 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,9 + ,6 + ,7 + ,60835 + ,365010 + ,325356 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,14 + ,1 + ,0 + ,61016 + ,61016 + ,325930 + ,0 + ,520 + ,10 + ,7 + ,58650 + ,586500 + ,318020 + ,0 + ,1766 + ,73 + ,46 + ,60438 + ,4411974 + ,326389 + ,0 + ,0 + ,0 + ,5 + ,60720 + ,0 + ,325560 + ,0 + ,458 + ,11 + ,7 + ,58625 + ,644875 + ,302925 + ,0 + ,20 + ,3 + ,2 + ,60938 + ,182814 + ,325540 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,98 + ,2 + ,2 + ,61490 + ,122980 + ,326736 + ,0 + ,405 + ,7 + ,5 + ,60845 + ,425915 + ,340580 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,483 + ,27 + ,7 + ,60830 + ,1642410 + ,331828 + ,0 + ,454 + ,51 + ,24 + ,63261 + ,3226311 + ,323299 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,757 + ,19 + ,18 + ,45689 + ,868091 + ,387722 + ,0 + ,0 + ,0 + ,4 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,36 + ,4 + ,3 + ,61564 + ,246256 + ,324598 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,203 + ,9 + ,9 + ,61938 + ,557442 + ,328726 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,90 + ,8 + ,6 + ,60951 + ,487608 + ,325043 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,71 + ,1 + ,0 + ,60745 + ,60745 + ,325806 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,972 + ,34 + ,19 + ,71642 + ,2435828 + ,387732 + ,1 + ,531 + ,10 + ,11 + ,71641 + ,716410 + ,349729 + ,0 + ,604 + ,38 + ,8 + ,55792 + ,2120096 + ,332202 + ,1 + ,283 + ,10 + ,15 + ,71873 + ,718730 + ,305442 + ,1 + ,23 + ,5 + ,5 + ,62555 + ,312775 + ,329537 + ,1 + ,638 + ,14 + ,11 + ,60370 + ,845180 + ,327055 + ,1 + ,699 + ,16 + ,10 + ,64873 + ,1037968 + ,356245 + ,0 + ,149 + ,5 + ,6 + ,62041 + ,310205 + ,328451 + ,0 + ,226 + ,5 + ,5 + ,65745 + ,328725 + ,307062 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,275 + ,4 + ,7 + ,59500 + ,238000 + ,331345 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,141 + ,6 + ,7 + ,61630 + ,369780 + ,331824 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,28 + ,2 + ,3 + ,60890 + ,121780 + ,325685 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,2566 + ,80 + ,74 + ,113521 + ,9081680 + ,404480 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,472 + ,20 + ,19 + ,80045 + ,1600900 + ,318314 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,203 + ,27 + ,12 + ,50804 + ,1371708 + ,311807 + ,1 + ,496 + ,17 + ,12 + ,87390 + ,1485630 + ,337724 + ,1 + ,10 + ,2 + ,5 + ,61656 + ,123312 + ,326431 + ,1 + ,63 + ,4 + ,2 + ,65688 + ,262752 + ,327556 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,1136 + ,32 + ,26 + ,48522 + ,1552704 + ,356850 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,267 + ,32 + ,11 + ,57640 + ,1844480 + ,322741 + ,0 + ,474 + ,20 + ,10 + ,61977 + ,1239540 + ,310902 + ,0 + ,534 + ,7 + ,5 + ,62620 + ,438340 + ,324295 + ,0 + ,0 + ,0 + ,2 + ,60720 + ,0 + ,325560 + ,0 + ,15 + ,8 + ,6 + ,60831 + ,486648 + ,326156 + ,0 + ,397 + ,28 + ,7 + ,60646 + ,1698088 + ,326960 + ,0 + ,0 + ,0 + ,2 + ,60720 + ,0 + ,325560 + ,0 + ,1061 + ,20 + ,22 + ,56225 + ,1124500 + ,333411 + ,0 + ,288 + ,4 + ,3 + ,60510 + ,242040 + ,297761 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,3 + ,2 + ,1 + ,60698 + ,121396 + ,325536 + ,0 + ,0 + ,0 + ,6 + ,60720 + ,0 + ,325560 + ,0 + ,20 + ,2 + ,1 + ,60805 + ,121610 + ,325762 + ,0 + ,278 + ,26 + ,22 + ,61404 + ,1596504 + ,327957 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,192 + ,4 + ,2 + ,65276 + ,261104 + ,318521 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,317 + ,9 + ,7 + ,63915 + ,575235 + ,319775 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,368 + ,17 + ,13 + ,61686 + ,1048662 + ,332128 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,2 + ,1 + ,0 + ,60743 + ,60743 + ,325486 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,53 + ,6 + ,6 + ,60349 + ,362094 + ,325838 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,94 + ,3 + ,3 + ,61360 + ,184080 + ,331767 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,24 + ,8 + ,7 + ,59818 + ,478544 + ,324523 + ,0 + ,2332 + ,4 + ,2 + ,72680 + ,290720 + ,339995 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,131 + ,11 + ,15 + ,61808 + ,679888 + ,319582 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,206 + ,9 + ,9 + ,53110 + ,477990 + ,307245 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,167 + ,2 + ,1 + ,64245 + ,128490 + ,317967 + ,0 + ,622 + ,73 + ,38 + ,73007 + ,5329511 + ,331488 + ,0 + ,885 + ,85 + ,49 + ,82732 + ,7032220 + ,335452 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,365 + ,8 + ,6 + ,54820 + ,438560 + ,334184 + ,0 + ,364 + ,35 + ,26 + ,47705 + ,1669675 + ,313213 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,226 + ,12 + ,13 + ,72835 + ,874020 + ,348678 + ,0 + ,307 + ,15 + ,10 + ,58856 + ,882840 + ,328727 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,188 + ,11 + ,9 + ,77655 + ,854205 + ,387978 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,138 + ,6 + ,26 + ,69817 + ,418902 + ,336704 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,125 + ,12 + ,19 + ,60798 + ,729576 + ,322076 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,282 + ,30 + ,12 + ,62452 + ,1873560 + ,334272 + ,0 + ,335 + ,33 + ,23 + ,64175 + ,2117775 + ,338197 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,813 + ,82 + ,24 + ,67440 + ,5530080 + ,321024 + ,0 + ,176 + ,28 + ,8 + ,68136 + ,1907808 + ,322145 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,249 + ,72 + ,26 + ,56726 + ,4084272 + ,323351 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,333 + ,13 + ,9 + ,70811 + ,920543 + ,327748 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,30 + ,4 + ,3 + ,62045 + ,248180 + ,328157 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,249 + ,62 + ,13 + ,54323 + ,3368026 + ,311594 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,165 + ,24 + ,12 + ,62841 + ,1508184 + ,335962 + ,0 + ,453 + ,21 + ,19 + ,81125 + ,1703625 + ,372426 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,53 + ,14 + ,10 + ,59506 + ,833084 + ,319844 + ,1 + ,382 + ,21 + ,9 + ,59365 + ,1246665 + ,355822 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,9 + ,60720 + ,0 + ,325560 + ,1 + ,30 + ,4 + ,4 + ,60798 + ,243192 + ,324047 + ,0 + ,290 + ,2 + ,1 + ,58790 + ,117580 + ,311464 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,366 + ,53 + ,14 + ,61808 + ,3275824 + ,353417 + ,0 + ,2 + ,9 + ,12 + ,60735 + ,546615 + ,325590 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,209 + ,13 + ,19 + ,64016 + ,832208 + ,328576 + ,0 + ,384 + ,22 + ,17 + ,54683 + ,1203026 + ,326126 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,365 + ,83 + ,32 + ,87192 + ,7236936 + ,369376 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,49 + ,8 + ,14 + ,64107 + ,512856 + ,332013 + ,0 + ,3 + ,4 + ,8 + ,60761 + ,243044 + ,325871 + ,0 + ,133 + ,14 + ,4 + ,65990 + ,923860 + ,342165 + ,0 + ,32 + ,1 + ,0 + ,59988 + ,59988 + ,324967 + ,0 + ,368 + ,17 + ,20 + ,61167 + ,1039839 + ,314832 + ,0 + ,1 + ,6 + ,5 + ,60719 + ,364314 + ,325557 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,22 + ,2 + ,1 + ,60722 + ,121444 + ,322649 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,96 + ,5 + ,4 + ,60379 + ,301895 + ,324598 + ,0 + ,1 + ,2 + ,1 + ,60727 + ,121454 + ,325567 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,81 + ,7 + ,4 + ,60925 + ,426475 + ,324005 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,26 + ,1 + ,1 + ,60896 + ,60896 + ,325748 + ,0 + ,125 + ,13 + ,10 + ,59734 + ,776542 + ,323385 + ,0 + ,304 + ,15 + ,12 + ,62969 + ,944535 + ,315409 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,119 + ,6 + ,3 + ,59118 + ,354708 + ,312275 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,312 + ,14 + ,3 + ,58598 + ,820372 + ,320576 + ,0 + ,60 + ,10 + ,7 + ,61124 + ,611240 + ,325246 + ,0 + ,587 + ,12 + ,10 + ,59595 + ,715140 + ,332961 + ,0 + ,135 + ,2 + ,1 + ,62065 + ,124130 + ,323010 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,514 + ,52 + ,15 + ,78780 + ,4096560 + ,345253 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,1 + ,4 + ,4 + ,60722 + ,242888 + ,325559 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,58 + ,3 + ,2 + ,61600 + ,184800 + ,319634 + ,0 + ,180 + ,11 + ,9 + ,59635 + ,655985 + ,319951 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,448 + ,40 + ,7 + ,59781 + ,2391240 + ,318519 + ,0 + ,227 + ,9 + ,7 + ,76644 + ,689796 + ,343222 + ,0 + ,174 + ,1 + ,3 + ,64820 + ,64820 + ,317234 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,121 + ,24 + ,11 + ,56178 + ,1348272 + ,314025 + ,0 + ,607 + ,11 + ,7 + ,60436 + ,664796 + ,320249 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,530 + ,60 + ,18 + ,73433 + ,4405980 + ,349365 + ,0 + ,571 + ,80 + ,14 + ,41477 + ,3318160 + ,289197 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,78 + ,16 + ,12 + ,62700 + ,1003200 + ,329245 + ,0 + ,2489 + ,40 + ,29 + ,67804 + ,2712160 + ,240869 + ,0 + ,131 + ,6 + ,3 + ,59661 + ,357966 + ,327182 + ,0 + ,923 + ,8 + ,6 + ,58620 + ,468960 + ,322876 + ,0 + ,72 + ,3 + ,3 + ,60398 + ,181194 + ,323117 + ,0 + ,572 + ,16 + ,8 + ,58580 + ,937280 + ,306351 + ,0 + ,397 + ,10 + ,10 + ,62710 + ,627100 + ,335137 + ,0 + ,450 + ,8 + ,6 + ,59325 + ,474600 + ,308271 + ,0 + ,622 + ,7 + ,8 + ,60950 + ,426650 + ,301731 + ,0 + ,694 + ,8 + ,6 + ,68060 + ,544480 + ,382409 + ,1 + ,3425 + ,12 + ,9 + ,83620 + ,1003440 + ,279230 + ,0 + ,562 + ,13 + ,8 + ,58456 + ,759928 + ,298731 + ,0 + ,4917 + ,42 + ,26 + ,52811 + ,2218062 + ,243650 + ,1 + ,1442 + ,118 + ,239 + ,121173 + ,14298414 + ,532682 + ,0 + ,529 + ,9 + ,7 + ,63870 + ,574830 + ,319771 + ,1 + ,2126 + ,138 + ,41 + ,21001 + ,2898138 + ,171493 + ,0 + ,1061 + ,5 + ,3 + ,70415 + ,352075 + ,347262 + ,0 + ,776 + ,9 + ,8 + ,64230 + ,578070 + ,343945 + ,0 + ,611 + ,8 + ,6 + ,59190 + ,473520 + ,311874 + ,1 + ,1526 + ,25 + ,21 + ,69351 + ,1733775 + ,302211 + ,0 + ,592 + ,7 + ,7 + ,64270 + ,449890 + ,316708 + ,0 + ,1182 + ,13 + ,11 + ,70694 + ,919022 + ,333463 + ,0 + ,621 + ,16 + ,11 + ,68005 + ,1088080 + ,344282 + ,0 + ,989 + ,11 + ,12 + ,58930 + ,648230 + ,319635 + ,0 + ,438 + ,11 + ,9 + ,58320 + ,641520 + ,301186 + ,0 + ,726 + ,3 + ,3 + ,69980 + ,209940 + ,300381 + ,0 + ,1303 + ,61 + ,57 + ,69863 + ,4261643 + ,318765 + ,1 + ,6341 + ,24 + ,16 + ,63255 + ,1518120 + ,286146 + ,1 + ,1164 + ,17 + ,15 + ,57320 + ,974440 + ,306844 + ,1 + ,3310 + ,33 + ,32 + ,75230 + ,2482590 + ,307705 + ,0 + ,1366 + ,7 + ,5 + ,79420 + ,555940 + ,312448 + ,0 + ,965 + ,3 + ,2 + ,73490 + ,220470 + ,299715 + ,0 + ,3256 + ,66 + ,23 + ,35250 + ,2326500 + ,373399 + ,1 + ,1135 + ,17 + ,20 + ,62285 + ,1058845 + ,299446 + ,0 + ,1270 + ,26 + ,24 + ,69206 + ,1799356 + ,325586 + ,0 + ,661 + ,3 + ,1 + ,65920 + ,197760 + ,291221 + ,0 + ,1013 + ,2 + ,1 + ,69770 + ,139540 + ,261173 + ,0 + ,2844 + ,67 + ,74 + ,72683 + ,4869761 + ,255027 + ,1 + ,11528 + ,70 + ,68 + ,-14545 + ,-1018150 + ,-78375 + ,0 + ,6526 + ,26 + ,20 + ,55830 + ,1451580 + ,-58143 + ,0 + ,2264 + ,24 + ,20 + ,55174 + ,1324176 + ,227033 + ,1 + ,4461 + ,94 + ,79 + ,67038 + ,6301572 + ,235098 + ,0 + ,3999 + ,30 + ,21 + ,51252 + ,1537560 + ,21267 + ,0 + ,35624 + ,223 + ,244 + ,157278 + ,35072994 + ,238675 + ,0 + ,9252 + ,48 + ,32 + ,79510 + ,3816480 + ,197687 + ,0 + ,15236 + ,90 + ,86 + ,77440 + ,6969600 + ,418341 + ,0 + ,18073 + ,180 + ,69 + ,27284 + ,4911120 + ,-297706) + ,dim=c(7 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Trades' + ,'Orders' + ,'Dividends' + ,'TrDiv' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(7,431),dimnames=list(c('Group','Costs','Trades','Orders','Dividends','TrDiv','Wealth'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '7' > #'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 Group Costs Trades Orders Dividends TrDiv 1 6282929 1 162556 1081 807 213118 230380558 2 4324047 1 29790 309 444 81767 25266003 3 4108272 1 87550 458 412 153198 70164684 4 -1212617 0 84738 588 428 -26007 -15292116 5 1485329 1 54660 299 312 126942 37955658 6 1779876 1 42634 156 166 157214 24525384 7 1367203 0 40949 481 263 129352 62218312 8 2519076 1 42312 323 237 234817 75845891 9 912684 1 37704 452 228 60448 27322496 10 1443586 1 16275 109 129 47818 5212162 11 1220017 0 25830 115 104 245546 28237790 12 984885 0 12679 110 122 48020 5282200 13 1457425 1 18014 239 393 -1710 -408690 14 -572920 0 43556 247 190 32648 8064056 15 929144 1 24524 497 275 95350 47388950 16 1151176 0 6532 103 62 151352 15589256 17 790090 0 7123 109 102 288170 31410530 18 774497 1 20813 502 255 114337 57397174 19 990576 1 37597 248 234 37884 9395232 20 454195 0 17821 373 277 122844 45820812 21 876607 1 12988 119 73 82340 9798460 22 711969 1 22330 84 67 79801 6703284 23 702380 0 13326 102 103 165548 16885896 24 264449 0 16189 295 290 116384 34333280 25 450033 0 7146 105 83 134028 14072940 26 541063 0 15824 64 56 63838 4085632 27 588864 1 26088 267 224 74996 20023932 28 -37216 0 11326 129 64 31080 4009320 29 783310 0 8568 37 34 32168 1190216 30 467359 0 14416 361 139 49857 17998377 31 688779 1 3369 28 26 87161 2440508 32 608419 1 11819 85 70 106113 9019605 33 696348 1 6620 44 40 80570 3545080 34 597793 1 4519 49 42 102129 5004321 35 821730 0 2220 22 12 301670 6636740 36 377934 0 18562 155 211 102313 15858515 37 651939 0 10327 91 74 88577 8060507 38 697458 1 5336 81 80 112477 9110637 39 700368 1 2365 79 83 191778 15150462 40 225986 0 4069 145 131 79804 11571580 41 348695 0 7710 816 187 128294 104687904 42 373683 0 13718 61 56 96448 5883328 43 501709 0 4525 226 89 93811 21201286 44 413743 0 6869 105 88 117520 12339600 45 379825 0 4628 62 39 69159 4287858 46 336260 1 3653 24 24 101792 2443008 47 636765 1 1265 26 19 210568 5474768 48 481231 1 7489 322 149 136996 44112712 49 469107 0 4901 84 58 121920 10241280 50 211928 0 2284 33 41 76403 2521299 51 563925 1 3160 108 90 108094 11674152 52 511939 1 4150 150 136 134759 20213850 53 521016 1 7285 115 97 188873 21720395 54 543856 1 1134 162 63 146216 23686992 55 329304 1 4658 158 114 156608 24744064 56 423262 0 2384 97 77 61348 5950756 57 509665 0 3748 9 6 50350 453150 58 455881 0 5371 66 47 87720 5789520 59 367772 0 1285 107 51 99489 10645323 60 406339 1 9327 101 85 87419 8829319 61 493408 1 5565 47 43 94355 4434685 62 232942 0 1528 38 32 60326 2292388 63 416002 1 3122 34 25 94670 3218780 64 337430 1 7317 84 74 82425 6923700 65 361517 0 2675 79 54 59017 4662343 66 360962 0 13253 947 251 90829 86015063 67 235561 0 880 74 15 80791 5978534 68 408247 1 2053 53 44 100423 5322419 69 450296 0 1424 94 73 131116 12324904 70 418799 1 4036 63 85 100269 6316947 71 247405 1 3045 58 49 27330 1585140 72 378519 0 5119 49 38 39039 1912911 73 326638 0 1431 34 35 106885 3634090 74 328233 0 554 11 9 79285 872135 75 386225 0 1975 35 34 118881 4160835 76 283662 1 1286 17 17 77623 1319591 77 370225 0 1012 47 29 114768 5394096 78 269236 0 810 43 11 74015 3182645 79 365732 0 1280 117 52 69465 8127405 80 420383 1 666 171 13 117869 20155599 81 345811 0 1380 26 29 60982 1585532 82 431809 1 4608 73 64 90131 6579563 83 418876 0 876 59 33 138971 8199289 84 297476 0 814 18 15 39625 713250 85 416776 0 514 15 15 102725 1540875 86 357257 1 5692 72 68 64239 4625208 87 458343 0 3642 86 100 90262 7762532 88 388386 0 540 14 13 103960 1455440 89 358934 0 2099 64 45 106611 6823104 90 407560 0 567 11 14 103345 1136795 91 392558 0 2001 52 36 95551 4968652 92 373177 1 2949 41 40 82903 3399023 93 428370 0 2253 99 68 63593 6295707 94 369419 1 6533 75 29 126910 9518250 95 358649 0 1889 45 43 37527 1688715 96 376641 1 3055 43 30 60247 2590621 97 467427 0 272 8 9 112995 903960 98 364885 1 1414 198 22 70184 13896432 99 436230 0 2564 22 19 130140 2863080 100 329118 1 1383 11 9 73221 805431 101 317365 1 1261 33 31 76114 2511762 102 286849 0 975 23 19 90534 2082282 103 376685 0 3366 80 55 108479 8678320 104 407198 0 576 18 8 113761 2047698 105 377772 0 1306 28 26 68696 1923488 106 271483 0 746 23 29 71561 1645903 107 153661 1 3192 60 48 59831 3589860 108 513294 1 2045 20 16 97890 1957800 109 324881 0 5477 59 45 101481 5987379 110 264512 1 1932 36 20 72954 2626344 111 420968 0 936 30 22 67939 2038170 112 129302 1 3437 47 33 48022 2257034 113 191521 0 5131 71 44 86111 6113881 114 268673 1 2397 14 13 74020 1036280 115 353179 1 1389 9 6 57530 517770 116 354624 0 1503 39 35 56364 2198196 117 363713 0 402 26 8 84990 2209740 118 456657 0 2239 21 17 88590 1860390 119 211742 1 2234 16 11 77200 1235200 120 338381 0 837 69 21 61262 4227078 121 418530 0 10579 92 92 110309 10148428 122 351483 0 875 14 12 67000 938000 123 372928 0 1395 103 108 93099 9589197 124 485538 1 1659 29 25 107577 3119733 125 279268 1 2647 37 17 62920 2328040 126 219060 1 3294 23 23 75832 1744136 127 325560 0 0 0 0 60720 0 128 325314 0 94 7 10 60793 425551 129 322046 0 422 28 23 57935 1622180 130 325560 0 0 0 0 60720 0 131 325599 0 34 8 7 60630 485040 132 377028 0 1558 63 25 55637 3505131 133 325560 0 0 0 0 60720 0 134 323850 0 43 3 20 60887 182661 135 325560 0 0 0 0 60720 0 136 331514 0 316 9 4 60505 544545 137 325632 0 115 13 10 60945 792285 138 325560 0 0 0 0 60720 0 139 325560 0 0 0 0 60720 0 140 325560 0 0 0 0 60720 0 141 322265 0 389 14 7 58990 825860 142 325560 0 0 0 0 60720 0 143 325906 0 1002 15 11 56750 851250 144 325985 0 36 3 4 60894 182682 145 346145 0 460 15 15 63346 950190 146 325898 0 309 11 9 56535 621885 147 325560 0 0 0 0 60720 0 148 325356 0 9 6 7 60835 365010 149 325560 0 0 0 0 60720 0 150 325930 0 14 1 0 61016 61016 151 318020 0 520 10 7 58650 586500 152 326389 0 1766 73 46 60438 4411974 153 325560 0 0 0 5 60720 0 154 302925 0 458 11 7 58625 644875 155 325540 0 20 3 2 60938 182814 156 325560 0 0 0 0 60720 0 157 325560 0 0 0 0 60720 0 158 326736 0 98 2 2 61490 122980 159 340580 0 405 7 5 60845 425915 160 325560 0 0 0 0 60720 0 161 325560 0 0 0 0 60720 0 162 325560 0 0 0 0 60720 0 163 325560 0 0 0 0 60720 0 164 331828 0 483 27 7 60830 1642410 165 323299 0 454 51 24 63261 3226311 166 325560 0 0 0 0 60720 0 167 325560 0 0 0 0 60720 0 168 387722 0 757 19 18 45689 868091 169 325560 0 0 0 4 60720 0 170 325560 0 0 0 0 60720 0 171 325560 0 0 0 0 60720 0 172 324598 0 36 4 3 61564 246256 173 325560 0 0 0 0 60720 0 174 328726 0 203 9 9 61938 557442 175 325560 0 0 0 0 60720 0 176 325043 0 90 8 6 60951 487608 177 325560 0 0 0 0 60720 0 178 325806 1 71 1 0 60745 60745 179 325560 0 0 0 0 60720 0 180 325560 0 0 0 0 60720 0 181 387732 0 972 34 19 71642 2435828 182 349729 1 531 10 11 71641 716410 183 332202 0 604 38 8 55792 2120096 184 305442 1 283 10 15 71873 718730 185 329537 1 23 5 5 62555 312775 186 327055 1 638 14 11 60370 845180 187 356245 1 699 16 10 64873 1037968 188 328451 0 149 5 6 62041 310205 189 307062 0 226 5 5 65745 328725 190 325560 0 0 0 0 60720 0 191 331345 0 275 4 7 59500 238000 192 325560 0 0 0 0 60720 0 193 331824 0 141 6 7 61630 369780 194 325560 0 0 0 0 60720 0 195 325685 0 28 2 3 60890 121780 196 325560 1 0 0 0 60720 0 197 404480 1 2566 80 74 113521 9081680 198 325560 1 0 0 0 60720 0 199 325560 1 0 0 0 60720 0 200 318314 1 472 20 19 80045 1600900 201 325560 1 0 0 0 60720 0 202 325560 1 0 0 0 60720 0 203 325560 1 0 0 0 60720 0 204 311807 1 203 27 12 50804 1371708 205 337724 1 496 17 12 87390 1485630 206 326431 1 10 2 5 61656 123312 207 327556 1 63 4 2 65688 262752 208 325560 1 0 0 0 60720 0 209 356850 1 1136 32 26 48522 1552704 210 325560 1 0 0 0 60720 0 211 325560 0 0 0 0 60720 0 212 325560 0 0 0 0 60720 0 213 322741 0 267 32 11 57640 1844480 214 310902 0 474 20 10 61977 1239540 215 324295 0 534 7 5 62620 438340 216 325560 0 0 0 2 60720 0 217 326156 0 15 8 6 60831 486648 218 326960 0 397 28 7 60646 1698088 219 325560 0 0 0 2 60720 0 220 333411 0 1061 20 22 56225 1124500 221 297761 0 288 4 3 60510 242040 222 325560 0 0 0 0 60720 0 223 325536 0 3 2 1 60698 121396 224 325560 0 0 0 6 60720 0 225 325762 0 20 2 1 60805 121610 226 327957 0 278 26 22 61404 1596504 227 325560 0 0 0 0 60720 0 228 325560 0 0 0 0 60720 0 229 318521 0 192 4 2 65276 261104 230 325560 0 0 0 0 60720 0 231 319775 0 317 9 7 63915 575235 232 325560 0 0 0 0 60720 0 233 325560 0 0 0 0 60720 0 234 332128 0 368 17 13 61686 1048662 235 325560 0 0 0 0 60720 0 236 325486 0 2 1 0 60743 60743 237 325560 0 0 0 0 60720 0 238 325838 0 53 6 6 60349 362094 239 325560 0 0 0 0 60720 0 240 325560 0 0 0 0 60720 0 241 325560 0 0 0 0 60720 0 242 331767 0 94 3 3 61360 184080 243 325560 0 0 0 0 60720 0 244 324523 0 24 8 7 59818 478544 245 339995 0 2332 4 2 72680 290720 246 325560 1 0 0 0 60720 0 247 325560 1 0 0 0 60720 0 248 319582 0 131 11 15 61808 679888 249 325560 1 0 0 0 60720 0 250 325560 1 0 0 0 60720 0 251 307245 0 206 9 9 53110 477990 252 325560 1 0 0 0 60720 0 253 317967 0 167 2 1 64245 128490 254 331488 0 622 73 38 73007 5329511 255 335452 0 885 85 49 82732 7032220 256 325560 0 0 0 0 60720 0 257 334184 0 365 8 6 54820 438560 258 313213 0 364 35 26 47705 1669675 259 325560 1 0 0 0 60720 0 260 325560 0 0 0 0 60720 0 261 325560 1 0 0 0 60720 0 262 325560 1 0 0 0 60720 0 263 348678 0 226 12 13 72835 874020 264 328727 0 307 15 10 58856 882840 265 325560 1 0 0 0 60720 0 266 325560 0 0 0 0 60720 0 267 325560 1 0 0 0 60720 0 268 387978 0 188 11 9 77655 854205 269 325560 1 0 0 0 60720 0 270 336704 0 138 6 26 69817 418902 271 325560 1 0 0 0 60720 0 272 325560 1 0 0 0 60720 0 273 325560 1 0 0 0 60720 0 274 322076 0 125 12 19 60798 729576 275 325560 0 0 0 0 60720 0 276 334272 0 282 30 12 62452 1873560 277 338197 0 335 33 23 64175 2117775 278 325560 0 0 0 0 60720 0 279 321024 1 813 82 24 67440 5530080 280 322145 0 176 28 8 68136 1907808 281 325560 1 0 0 0 60720 0 282 325560 0 0 0 0 60720 0 283 323351 0 249 72 26 56726 4084272 284 325560 0 0 0 0 60720 0 285 327748 0 333 13 9 70811 920543 286 325560 0 0 0 0 60720 0 287 325560 1 0 0 0 60720 0 288 328157 0 30 4 3 62045 248180 289 325560 1 0 0 0 60720 0 290 311594 0 249 62 13 54323 3368026 291 325560 1 0 0 0 60720 0 292 335962 0 165 24 12 62841 1508184 293 372426 0 453 21 19 81125 1703625 294 325560 0 0 0 0 60720 0 295 319844 0 53 14 10 59506 833084 296 355822 1 382 21 9 59365 1246665 297 325560 0 0 0 0 60720 0 298 325560 1 0 0 0 60720 0 299 325560 0 0 0 0 60720 0 300 325560 0 0 0 9 60720 0 301 324047 1 30 4 4 60798 243192 302 311464 0 290 2 1 58790 117580 303 325560 0 0 0 0 60720 0 304 325560 1 0 0 0 60720 0 305 353417 0 366 53 14 61808 3275824 306 325590 0 2 9 12 60735 546615 307 325560 1 0 0 0 60720 0 308 328576 1 209 13 19 64016 832208 309 326126 0 384 22 17 54683 1203026 310 325560 1 0 0 0 60720 0 311 325560 1 0 0 0 60720 0 312 369376 0 365 83 32 87192 7236936 313 325560 1 0 0 0 60720 0 314 332013 1 49 8 14 64107 512856 315 325871 0 3 4 8 60761 243044 316 342165 0 133 14 4 65990 923860 317 324967 0 32 1 0 59988 59988 318 314832 0 368 17 20 61167 1039839 319 325557 0 1 6 5 60719 364314 320 325560 0 0 0 0 60720 0 321 325560 0 0 0 0 60720 0 322 325560 0 0 0 0 60720 0 323 325560 0 0 0 0 60720 0 324 325560 1 0 0 0 60720 0 325 325560 0 0 0 0 60720 0 326 322649 0 22 2 1 60722 121444 327 325560 1 0 0 0 60720 0 328 325560 1 0 0 0 60720 0 329 325560 0 0 0 0 60720 0 330 325560 0 0 0 0 60720 0 331 325560 0 0 0 0 60720 0 332 325560 0 0 0 0 60720 0 333 325560 0 0 0 0 60720 0 334 324598 0 96 5 4 60379 301895 335 325567 0 1 2 1 60727 121454 336 325560 0 0 0 0 60720 0 337 324005 0 81 7 4 60925 426475 338 325560 0 0 0 0 60720 0 339 325748 0 26 1 1 60896 60896 340 323385 0 125 13 10 59734 776542 341 315409 0 304 15 12 62969 944535 342 325560 0 0 0 0 60720 0 343 325560 0 0 0 0 60720 0 344 325560 0 0 0 0 60720 0 345 325560 0 0 0 0 60720 0 346 325560 0 0 0 0 60720 0 347 312275 0 119 6 3 59118 354708 348 325560 0 0 0 0 60720 0 349 325560 0 0 0 0 60720 0 350 325560 0 0 0 0 60720 0 351 320576 0 312 14 3 58598 820372 352 325246 0 60 10 7 61124 611240 353 332961 0 587 12 10 59595 715140 354 323010 0 135 2 1 62065 124130 355 325560 0 0 0 0 60720 0 356 325560 0 0 0 0 60720 0 357 345253 0 514 52 15 78780 4096560 358 325560 0 0 0 0 60720 0 359 325560 0 0 0 0 60720 0 360 325560 0 0 0 0 60720 0 361 325559 0 1 4 4 60722 242888 362 325560 0 0 0 0 60720 0 363 325560 0 0 0 0 60720 0 364 319634 1 58 3 2 61600 184800 365 319951 0 180 11 9 59635 655985 366 325560 0 0 0 0 60720 0 367 325560 0 0 0 0 60720 0 368 325560 0 0 0 0 60720 0 369 325560 0 0 0 0 60720 0 370 325560 0 0 0 0 60720 0 371 325560 0 0 0 0 60720 0 372 318519 0 448 40 7 59781 2391240 373 343222 0 227 9 7 76644 689796 374 317234 0 174 1 3 64820 64820 375 325560 0 0 0 0 60720 0 376 325560 0 0 0 0 60720 0 377 314025 0 121 24 11 56178 1348272 378 320249 0 607 11 7 60436 664796 379 325560 0 0 0 0 60720 0 380 325560 0 0 0 0 60720 0 381 325560 0 0 0 0 60720 0 382 349365 0 530 60 18 73433 4405980 383 289197 0 571 80 14 41477 3318160 384 325560 0 0 0 0 60720 0 385 329245 0 78 16 12 62700 1003200 386 240869 0 2489 40 29 67804 2712160 387 327182 0 131 6 3 59661 357966 388 322876 0 923 8 6 58620 468960 389 323117 0 72 3 3 60398 181194 390 306351 0 572 16 8 58580 937280 391 335137 0 397 10 10 62710 627100 392 308271 0 450 8 6 59325 474600 393 301731 0 622 7 8 60950 426650 394 382409 0 694 8 6 68060 544480 395 279230 1 3425 12 9 83620 1003440 396 298731 0 562 13 8 58456 759928 397 243650 0 4917 42 26 52811 2218062 398 532682 1 1442 118 239 121173 14298414 399 319771 0 529 9 7 63870 574830 400 171493 1 2126 138 41 21001 2898138 401 347262 0 1061 5 3 70415 352075 402 343945 0 776 9 8 64230 578070 403 311874 0 611 8 6 59190 473520 404 302211 1 1526 25 21 69351 1733775 405 316708 0 592 7 7 64270 449890 406 333463 0 1182 13 11 70694 919022 407 344282 0 621 16 11 68005 1088080 408 319635 0 989 11 12 58930 648230 409 301186 0 438 11 9 58320 641520 410 300381 0 726 3 3 69980 209940 411 318765 0 1303 61 57 69863 4261643 412 286146 1 6341 24 16 63255 1518120 413 306844 1 1164 17 15 57320 974440 414 307705 1 3310 33 32 75230 2482590 415 312448 0 1366 7 5 79420 555940 416 299715 0 965 3 2 73490 220470 417 373399 0 3256 66 23 35250 2326500 418 299446 1 1135 17 20 62285 1058845 419 325586 0 1270 26 24 69206 1799356 420 291221 0 661 3 1 65920 197760 421 261173 0 1013 2 1 69770 139540 422 255027 0 2844 67 74 72683 4869761 423 -78375 1 11528 70 68 -14545 -1018150 424 -58143 0 6526 26 20 55830 1451580 425 227033 0 2264 24 20 55174 1324176 426 235098 1 4461 94 79 67038 6301572 427 21267 0 3999 30 21 51252 1537560 428 238675 0 35624 223 244 157278 35072994 429 197687 0 9252 48 32 79510 3816480 430 418341 0 15236 90 86 77440 6969600 431 -297706 0 18073 180 69 27284 4911120 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Orders Dividends 3.354e+05 5.511e+04 5.927e+00 -4.315e+03 4.170e+03 -5.567e-01 TrDiv 2.875e-02 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -1284038 -31164 18420 23944 2557838 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 3.354e+05 3.252e+04 10.314 < 2e-16 *** Group 5.511e+04 2.528e+04 2.180 0.02979 * Costs 5.927e+00 2.091e+00 2.834 0.00482 ** Trades -4.315e+03 2.752e+02 -15.680 < 2e-16 *** Orders 4.170e+03 4.287e+02 9.726 < 2e-16 *** Dividends -5.567e-01 4.374e-01 -1.273 0.20380 TrDiv 2.875e-02 1.767e-03 16.273 < 2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 228400 on 424 degrees of freedom Multiple R-squared: 0.7431, Adjusted R-squared: 0.7395 F-statistic: 204.4 on 6 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 6.002798e-36 3.001399e-36 [2,] 1.0000000 2.824924e-37 1.412462e-37 [3,] 1.0000000 1.595201e-38 7.976005e-39 [4,] 1.0000000 1.467970e-62 7.339852e-63 [5,] 1.0000000 2.653330e-76 1.326665e-76 [6,] 1.0000000 1.195983e-80 5.979917e-81 [7,] 1.0000000 2.246340e-88 1.123170e-88 [8,] 1.0000000 5.632362e-94 2.816181e-94 [9,] 1.0000000 3.075679e-97 1.537839e-97 [10,] 1.0000000 5.127542e-104 2.563771e-104 [11,] 1.0000000 6.942451e-108 3.471226e-108 [12,] 1.0000000 3.627250e-113 1.813625e-113 [13,] 1.0000000 8.517863e-117 4.258931e-117 [14,] 1.0000000 5.687746e-117 2.843873e-117 [15,] 1.0000000 1.112129e-125 5.560644e-126 [16,] 1.0000000 9.704956e-125 4.852478e-125 [17,] 1.0000000 4.084167e-127 2.042084e-127 [18,] 1.0000000 2.740490e-133 1.370245e-133 [19,] 1.0000000 1.323807e-134 6.619034e-135 [20,] 1.0000000 9.832949e-149 4.916474e-149 [21,] 1.0000000 2.362293e-159 1.181146e-159 [22,] 1.0000000 7.670821e-164 3.835410e-164 [23,] 1.0000000 7.720804e-168 3.860402e-168 [24,] 1.0000000 3.601832e-175 1.800916e-175 [25,] 1.0000000 1.510464e-177 7.552322e-178 [26,] 1.0000000 8.266993e-178 4.133496e-178 [27,] 1.0000000 1.566351e-180 7.831755e-181 [28,] 1.0000000 2.428257e-188 1.214128e-188 [29,] 1.0000000 3.309918e-194 1.654959e-194 [30,] 1.0000000 1.556186e-195 7.780932e-196 [31,] 1.0000000 1.292405e-196 6.462024e-197 [32,] 1.0000000 7.452755e-199 3.726378e-199 [33,] 1.0000000 2.325587e-199 1.162793e-199 [34,] 1.0000000 2.289329e-200 1.144665e-200 [35,] 1.0000000 2.268201e-199 1.134101e-199 [36,] 1.0000000 4.373332e-199 2.186666e-199 [37,] 1.0000000 3.439152e-199 1.719576e-199 [38,] 1.0000000 1.426529e-198 7.132645e-199 [39,] 1.0000000 9.397522e-200 4.698761e-200 [40,] 1.0000000 4.364562e-199 2.182281e-199 [41,] 1.0000000 1.730090e-199 8.650450e-200 [42,] 1.0000000 9.464975e-201 4.732487e-201 [43,] 1.0000000 6.352046e-202 3.176023e-202 [44,] 1.0000000 7.366366e-203 3.683183e-203 [45,] 1.0000000 6.049162e-202 3.024581e-202 [46,] 1.0000000 2.334303e-208 1.167151e-208 [47,] 1.0000000 2.294666e-208 1.147333e-208 [48,] 1.0000000 7.090789e-214 3.545395e-214 [49,] 1.0000000 6.656088e-215 3.328044e-215 [50,] 1.0000000 4.429154e-214 2.214577e-214 [51,] 1.0000000 2.545440e-215 1.272720e-215 [52,] 1.0000000 1.111902e-217 5.559511e-218 [53,] 1.0000000 2.703935e-217 1.351968e-217 [54,] 1.0000000 5.974204e-217 2.987102e-217 [55,] 1.0000000 7.994260e-217 3.997130e-217 [56,] 1.0000000 2.594442e-216 1.297221e-216 [57,] 1.0000000 1.454873e-219 7.274365e-220 [58,] 1.0000000 6.123754e-220 3.061877e-220 [59,] 1.0000000 4.614312e-219 2.307156e-219 [60,] 1.0000000 2.513912e-218 1.256956e-218 [61,] 1.0000000 5.599035e-218 2.799518e-218 [62,] 1.0000000 3.015624e-217 1.507812e-217 [63,] 1.0000000 1.260718e-219 6.303592e-220 [64,] 1.0000000 3.505790e-219 1.752895e-219 [65,] 1.0000000 3.397266e-218 1.698633e-218 [66,] 1.0000000 3.137009e-217 1.568505e-217 [67,] 1.0000000 1.194170e-216 5.970851e-217 [68,] 1.0000000 4.991704e-216 2.495852e-216 [69,] 1.0000000 1.469281e-215 7.346403e-216 [70,] 1.0000000 9.691553e-215 4.845776e-215 [71,] 1.0000000 7.731077e-215 3.865538e-215 [72,] 1.0000000 6.471011e-214 3.235505e-214 [73,] 1.0000000 4.484858e-214 2.242429e-214 [74,] 1.0000000 4.768882e-214 2.384441e-214 [75,] 1.0000000 4.877578e-213 2.438789e-213 [76,] 1.0000000 3.971417e-212 1.985708e-212 [77,] 1.0000000 8.543522e-213 4.271761e-213 [78,] 1.0000000 9.944656e-213 4.972328e-213 [79,] 1.0000000 9.597108e-212 4.798554e-212 [80,] 1.0000000 3.929530e-211 1.964765e-211 [81,] 1.0000000 3.507644e-210 1.753822e-210 [82,] 1.0000000 3.620712e-209 1.810356e-209 [83,] 1.0000000 1.668600e-208 8.343002e-209 [84,] 1.0000000 1.980987e-209 9.904936e-210 [85,] 1.0000000 1.479132e-208 7.395658e-209 [86,] 1.0000000 2.510024e-209 1.255012e-209 [87,] 1.0000000 7.771629e-210 3.885815e-210 [88,] 1.0000000 2.226946e-209 1.113473e-209 [89,] 1.0000000 7.982705e-209 3.991352e-209 [90,] 1.0000000 6.597699e-208 3.298850e-208 [91,] 1.0000000 6.318079e-207 3.159040e-207 [92,] 1.0000000 4.705927e-206 2.352964e-206 [93,] 1.0000000 6.327862e-206 3.163931e-206 [94,] 1.0000000 4.403444e-205 2.201722e-205 [95,] 1.0000000 3.590693e-204 1.795346e-204 [96,] 1.0000000 1.417922e-203 7.089608e-204 [97,] 1.0000000 3.653709e-203 1.826854e-203 [98,] 1.0000000 2.092618e-204 1.046309e-204 [99,] 1.0000000 6.477906e-207 3.238953e-207 [100,] 1.0000000 6.875801e-206 3.437901e-206 [101,] 1.0000000 2.510838e-205 1.255419e-205 [102,] 1.0000000 1.055984e-205 5.279918e-206 [103,] 1.0000000 7.895318e-207 3.947659e-207 [104,] 1.0000000 4.953772e-208 2.476886e-208 [105,] 1.0000000 2.710585e-207 1.355292e-207 [106,] 1.0000000 8.919880e-207 4.459940e-207 [107,] 1.0000000 2.636587e-206 1.318293e-206 [108,] 1.0000000 2.469933e-205 1.234967e-205 [109,] 1.0000000 9.358925e-207 4.679463e-207 [110,] 1.0000000 1.396874e-207 6.984371e-208 [111,] 1.0000000 1.016909e-206 5.084546e-207 [112,] 1.0000000 2.280916e-208 1.140458e-208 [113,] 1.0000000 1.660128e-207 8.300642e-208 [114,] 1.0000000 5.227016e-207 2.613508e-207 [115,] 1.0000000 2.302634e-208 1.151317e-208 [116,] 1.0000000 2.451444e-207 1.225722e-207 [117,] 1.0000000 2.666786e-207 1.333393e-207 [118,] 1.0000000 2.910706e-206 1.455353e-206 [119,] 1.0000000 3.280003e-205 1.640001e-205 [120,] 1.0000000 3.693380e-204 1.846690e-204 [121,] 1.0000000 3.990237e-203 1.995119e-203 [122,] 1.0000000 4.360552e-202 2.180276e-202 [123,] 1.0000000 1.544943e-202 7.724715e-203 [124,] 1.0000000 1.677675e-201 8.388377e-202 [125,] 1.0000000 1.868701e-200 9.343507e-201 [126,] 1.0000000 2.010425e-199 1.005212e-199 [127,] 1.0000000 2.030963e-198 1.015481e-198 [128,] 1.0000000 2.200412e-197 1.100206e-197 [129,] 1.0000000 2.341951e-196 1.170975e-196 [130,] 1.0000000 2.483175e-195 1.241587e-195 [131,] 1.0000000 2.622741e-194 1.311370e-194 [132,] 1.0000000 2.741870e-193 1.370935e-193 [133,] 1.0000000 2.873026e-192 1.436513e-192 [134,] 1.0000000 2.521973e-191 1.260986e-191 [135,] 1.0000000 2.647161e-190 1.323581e-190 [136,] 1.0000000 2.216655e-189 1.108328e-189 [137,] 1.0000000 2.246036e-188 1.123018e-188 [138,] 1.0000000 2.306315e-187 1.153157e-187 [139,] 1.0000000 2.386137e-186 1.193068e-186 [140,] 1.0000000 2.428429e-185 1.214214e-185 [141,] 1.0000000 2.454408e-184 1.227204e-184 [142,] 1.0000000 2.498241e-183 1.249120e-183 [143,] 1.0000000 1.846172e-182 9.230859e-183 [144,] 1.0000000 1.874281e-181 9.371405e-182 [145,] 1.0000000 1.766246e-180 8.831229e-181 [146,] 1.0000000 1.752562e-179 8.762812e-180 [147,] 1.0000000 1.729442e-178 8.647208e-179 [148,] 1.0000000 1.699032e-177 8.495161e-178 [149,] 1.0000000 1.664753e-176 8.323765e-177 [150,] 1.0000000 1.350890e-175 6.754449e-176 [151,] 1.0000000 1.310423e-174 6.552114e-175 [152,] 1.0000000 1.265427e-173 6.327134e-174 [153,] 1.0000000 1.216441e-172 6.082205e-173 [154,] 1.0000000 1.164041e-171 5.820206e-172 [155,] 1.0000000 9.803631e-171 4.901816e-171 [156,] 1.0000000 9.153033e-170 4.576516e-170 [157,] 1.0000000 8.648284e-169 4.324142e-169 [158,] 1.0000000 8.133850e-168 4.066925e-168 [159,] 1.0000000 2.607154e-168 1.303577e-168 [160,] 1.0000000 2.485401e-167 1.242701e-167 [161,] 1.0000000 2.337313e-166 1.168657e-166 [162,] 1.0000000 2.187902e-165 1.093951e-165 [163,] 1.0000000 2.041177e-164 1.020588e-164 [164,] 1.0000000 1.892998e-163 9.464989e-164 [165,] 1.0000000 1.746163e-162 8.730814e-163 [166,] 1.0000000 1.604275e-161 8.021376e-162 [167,] 1.0000000 1.472332e-160 7.361659e-161 [168,] 1.0000000 1.340108e-159 6.700541e-160 [169,] 1.0000000 1.184561e-158 5.922804e-159 [170,] 1.0000000 1.067843e-157 5.339216e-158 [171,] 1.0000000 9.581788e-157 4.790894e-157 [172,] 1.0000000 2.085174e-156 1.042587e-156 [173,] 1.0000000 1.594143e-155 7.970716e-156 [174,] 1.0000000 9.994779e-155 4.997389e-155 [175,] 1.0000000 5.513397e-154 2.756699e-154 [176,] 1.0000000 4.752891e-153 2.376446e-153 [177,] 1.0000000 3.883291e-152 1.941645e-152 [178,] 1.0000000 2.119812e-151 1.059906e-151 [179,] 1.0000000 1.853598e-150 9.267988e-151 [180,] 1.0000000 1.407746e-149 7.038730e-150 [181,] 1.0000000 1.220464e-148 6.102320e-149 [182,] 1.0000000 1.011774e-147 5.058870e-148 [183,] 1.0000000 8.692031e-147 4.346015e-147 [184,] 1.0000000 7.321939e-146 3.660970e-146 [185,] 1.0000000 6.232571e-145 3.116286e-145 [186,] 1.0000000 5.299201e-144 2.649601e-144 [187,] 1.0000000 4.385964e-143 2.192982e-143 [188,] 1.0000000 2.157530e-142 1.078765e-142 [189,] 1.0000000 1.779330e-141 8.896650e-142 [190,] 1.0000000 1.462099e-140 7.310494e-141 [191,] 1.0000000 8.171960e-140 4.085980e-140 [192,] 1.0000000 6.674793e-139 3.337396e-139 [193,] 1.0000000 5.430631e-138 2.715316e-138 [194,] 1.0000000 4.400840e-137 2.200420e-137 [195,] 1.0000000 3.558390e-136 1.779195e-136 [196,] 1.0000000 2.442963e-135 1.221482e-135 [197,] 1.0000000 1.929561e-134 9.647805e-135 [198,] 1.0000000 1.527932e-133 7.639658e-134 [199,] 1.0000000 1.213008e-132 6.065040e-133 [200,] 1.0000000 1.783082e-132 8.915412e-133 [201,] 1.0000000 1.418570e-131 7.092849e-132 [202,] 1.0000000 1.128402e-130 5.642008e-131 [203,] 1.0000000 8.935136e-130 4.467568e-130 [204,] 1.0000000 6.712429e-129 3.356214e-129 [205,] 1.0000000 5.185937e-128 2.592968e-128 [206,] 1.0000000 4.013507e-127 2.006754e-127 [207,] 1.0000000 3.132851e-126 1.566426e-126 [208,] 1.0000000 2.430512e-125 1.215256e-125 [209,] 1.0000000 1.759180e-124 8.795901e-125 [210,] 1.0000000 1.354653e-123 6.773263e-124 [211,] 1.0000000 7.173323e-123 3.586661e-123 [212,] 1.0000000 4.729795e-122 2.364897e-122 [213,] 1.0000000 3.588150e-121 1.794075e-121 [214,] 1.0000000 2.708732e-120 1.354366e-120 [215,] 1.0000000 2.049289e-119 1.024644e-119 [216,] 1.0000000 1.532110e-118 7.660552e-119 [217,] 1.0000000 1.135773e-117 5.678867e-118 [218,] 1.0000000 8.417690e-117 4.208845e-117 [219,] 1.0000000 6.209841e-116 3.104920e-116 [220,] 1.0000000 4.453989e-115 2.226995e-115 [221,] 1.0000000 3.255831e-114 1.627916e-114 [222,] 1.0000000 2.365443e-113 1.182721e-113 [223,] 1.0000000 1.713044e-112 8.565219e-113 [224,] 1.0000000 1.234800e-111 6.174000e-112 [225,] 1.0000000 8.511066e-111 4.255533e-111 [226,] 1.0000000 6.078732e-110 3.039366e-110 [227,] 1.0000000 4.317244e-109 2.158622e-109 [228,] 1.0000000 3.054662e-108 1.527331e-108 [229,] 1.0000000 2.155014e-107 1.077507e-107 [230,] 1.0000000 1.510410e-106 7.552051e-107 [231,] 1.0000000 1.053633e-105 5.268163e-106 [232,] 1.0000000 7.315254e-105 3.657627e-105 [233,] 1.0000000 4.967653e-104 2.483826e-104 [234,] 1.0000000 3.416655e-103 1.708327e-103 [235,] 1.0000000 2.345060e-102 1.172530e-102 [236,] 1.0000000 9.635195e-102 4.817598e-102 [237,] 1.0000000 6.502954e-101 3.251477e-101 [238,] 1.0000000 4.369610e-100 2.184805e-100 [239,] 1.0000000 2.928495e-99 1.464248e-99 [240,] 1.0000000 1.949721e-98 9.748605e-99 [241,] 1.0000000 1.292247e-97 6.461236e-98 [242,] 1.0000000 8.585482e-97 4.292741e-97 [243,] 1.0000000 5.637233e-96 2.818617e-96 [244,] 1.0000000 3.649198e-95 1.824599e-95 [245,] 1.0000000 2.210509e-94 1.105254e-94 [246,] 1.0000000 8.656227e-94 4.328113e-94 [247,] 1.0000000 5.626883e-93 2.813441e-93 [248,] 1.0000000 3.095274e-92 1.547637e-92 [249,] 1.0000000 1.809838e-91 9.049189e-92 [250,] 1.0000000 1.155237e-90 5.776184e-91 [251,] 1.0000000 7.373300e-90 3.686650e-90 [252,] 1.0000000 4.661453e-89 2.330727e-89 [253,] 1.0000000 2.933174e-88 1.466587e-88 [254,] 1.0000000 1.813141e-87 9.065706e-88 [255,] 1.0000000 1.091623e-86 5.458116e-87 [256,] 1.0000000 6.771084e-86 3.385542e-86 [257,] 1.0000000 4.195986e-85 2.097993e-85 [258,] 1.0000000 2.577269e-84 1.288634e-84 [259,] 1.0000000 1.004948e-83 5.024741e-84 [260,] 1.0000000 6.129898e-83 3.064949e-83 [261,] 1.0000000 3.681522e-82 1.840761e-82 [262,] 1.0000000 2.223619e-81 1.111810e-81 [263,] 1.0000000 1.336425e-80 6.682127e-81 [264,] 1.0000000 7.992157e-80 3.996078e-80 [265,] 1.0000000 4.770203e-79 2.385101e-79 [266,] 1.0000000 2.832675e-78 1.416338e-78 [267,] 1.0000000 1.606535e-77 8.032676e-78 [268,] 1.0000000 9.158524e-77 4.579262e-77 [269,] 1.0000000 5.357910e-76 2.678955e-76 [270,] 1.0000000 2.908977e-75 1.454488e-75 [271,] 1.0000000 1.599946e-74 7.999729e-75 [272,] 1.0000000 9.184860e-74 4.592430e-74 [273,] 1.0000000 5.267360e-73 2.633680e-73 [274,] 1.0000000 2.809018e-72 1.404509e-72 [275,] 1.0000000 1.595049e-71 7.975247e-72 [276,] 1.0000000 8.959682e-71 4.479841e-71 [277,] 1.0000000 5.034532e-70 2.517266e-70 [278,] 1.0000000 2.801642e-69 1.400821e-69 [279,] 1.0000000 1.555724e-68 7.778620e-69 [280,] 1.0000000 8.565057e-68 4.282529e-68 [281,] 1.0000000 4.406702e-67 2.203351e-67 [282,] 1.0000000 2.401065e-66 1.200533e-66 [283,] 1.0000000 1.271454e-65 6.357270e-66 [284,] 1.0000000 6.229709e-65 3.114854e-65 [285,] 1.0000000 3.358146e-64 1.679073e-64 [286,] 1.0000000 1.805286e-63 9.026428e-64 [287,] 1.0000000 7.432713e-63 3.716357e-63 [288,] 1.0000000 3.949071e-62 1.974536e-62 [289,] 1.0000000 2.081487e-61 1.040744e-61 [290,] 1.0000000 1.093651e-60 5.468254e-61 [291,] 1.0000000 5.724377e-60 2.862188e-60 [292,] 1.0000000 2.953599e-59 1.476799e-59 [293,] 1.0000000 1.525815e-58 7.629074e-59 [294,] 1.0000000 7.837902e-58 3.918951e-58 [295,] 1.0000000 3.993526e-57 1.996763e-57 [296,] 1.0000000 1.537549e-56 7.687743e-57 [297,] 1.0000000 7.801712e-56 3.900856e-56 [298,] 1.0000000 3.916165e-55 1.958082e-55 [299,] 1.0000000 1.915961e-54 9.579807e-55 [300,] 1.0000000 8.885438e-54 4.442719e-54 [301,] 1.0000000 4.385613e-53 2.192806e-53 [302,] 1.0000000 2.151387e-52 1.075694e-52 [303,] 1.0000000 8.878233e-52 4.439116e-52 [304,] 1.0000000 4.309386e-51 2.154693e-51 [305,] 1.0000000 2.058336e-50 1.029168e-50 [306,] 1.0000000 9.918849e-50 4.959425e-50 [307,] 1.0000000 4.570749e-49 2.285375e-49 [308,] 1.0000000 2.169033e-48 1.084516e-48 [309,] 1.0000000 1.022335e-47 5.111677e-48 [310,] 1.0000000 4.804484e-47 2.402242e-47 [311,] 1.0000000 2.240942e-46 1.120471e-46 [312,] 1.0000000 1.038664e-45 5.193319e-46 [313,] 1.0000000 4.783678e-45 2.391839e-45 [314,] 1.0000000 2.189137e-44 1.094569e-44 [315,] 1.0000000 9.927348e-44 4.963674e-44 [316,] 1.0000000 4.485069e-43 2.242534e-43 [317,] 1.0000000 2.016889e-42 1.008444e-42 [318,] 1.0000000 8.965039e-42 4.482520e-42 [319,] 1.0000000 3.955446e-41 1.977723e-41 [320,] 1.0000000 1.740893e-40 8.704467e-41 [321,] 1.0000000 7.610872e-40 3.805436e-40 [322,] 1.0000000 3.304914e-39 1.652457e-39 [323,] 1.0000000 1.425367e-38 7.126836e-39 [324,] 1.0000000 6.105360e-38 3.052680e-38 [325,] 1.0000000 2.595886e-37 1.297943e-37 [326,] 1.0000000 1.096565e-36 5.482826e-37 [327,] 1.0000000 4.599661e-36 2.299831e-36 [328,] 1.0000000 1.917371e-35 9.586853e-36 [329,] 1.0000000 7.928914e-35 3.964457e-35 [330,] 1.0000000 3.254633e-34 1.627316e-34 [331,] 1.0000000 1.326541e-33 6.632707e-34 [332,] 1.0000000 5.362006e-33 2.681003e-33 [333,] 1.0000000 2.153727e-32 1.076864e-32 [334,] 1.0000000 8.586370e-32 4.293185e-32 [335,] 1.0000000 3.397456e-31 1.698728e-31 [336,] 1.0000000 1.334113e-30 6.670565e-31 [337,] 1.0000000 5.198691e-30 2.599346e-30 [338,] 1.0000000 2.009276e-29 1.004638e-29 [339,] 1.0000000 7.709729e-29 3.854865e-29 [340,] 1.0000000 2.934955e-28 1.467478e-28 [341,] 1.0000000 1.108387e-27 5.541934e-28 [342,] 1.0000000 4.123429e-27 2.061714e-27 [343,] 1.0000000 1.533604e-26 7.668019e-27 [344,] 1.0000000 5.318953e-26 2.659476e-26 [345,] 1.0000000 1.948710e-25 9.743552e-26 [346,] 1.0000000 7.066181e-25 3.533090e-25 [347,] 1.0000000 2.540468e-24 1.270234e-24 [348,] 1.0000000 8.778946e-24 4.389473e-24 [349,] 1.0000000 3.104047e-23 1.552024e-23 [350,] 1.0000000 1.087860e-22 5.439302e-23 [351,] 1.0000000 3.778587e-22 1.889294e-22 [352,] 1.0000000 1.301628e-21 6.508142e-22 [353,] 1.0000000 4.439317e-21 2.219658e-21 [354,] 1.0000000 1.500047e-20 7.500233e-21 [355,] 1.0000000 4.963784e-20 2.481892e-20 [356,] 1.0000000 1.649745e-19 8.248726e-20 [357,] 1.0000000 5.419878e-19 2.709939e-19 [358,] 1.0000000 1.763185e-18 8.815926e-19 [359,] 1.0000000 5.679108e-18 2.839554e-18 [360,] 1.0000000 1.810806e-17 9.054031e-18 [361,] 1.0000000 5.714874e-17 2.857437e-17 [362,] 1.0000000 1.784904e-16 8.924522e-17 [363,] 1.0000000 5.465338e-16 2.732669e-16 [364,] 1.0000000 1.675998e-15 8.379989e-16 [365,] 1.0000000 5.084768e-15 2.542384e-15 [366,] 1.0000000 1.521808e-14 7.609042e-15 [367,] 1.0000000 4.503430e-14 2.251715e-14 [368,] 1.0000000 1.320324e-13 6.601618e-14 [369,] 1.0000000 3.789784e-13 1.894892e-13 [370,] 1.0000000 1.082890e-12 5.414452e-13 [371,] 1.0000000 3.056873e-12 1.528436e-12 [372,] 1.0000000 8.522921e-12 4.261461e-12 [373,] 1.0000000 2.332347e-11 1.166173e-11 [374,] 1.0000000 5.600820e-11 2.800410e-11 [375,] 1.0000000 1.507013e-10 7.535066e-11 [376,] 1.0000000 3.996184e-10 1.998092e-10 [377,] 1.0000000 9.117798e-10 4.558899e-10 [378,] 1.0000000 2.348823e-09 1.174412e-09 [379,] 1.0000000 5.804373e-09 2.902187e-09 [380,] 1.0000000 1.466598e-08 7.332992e-09 [381,] 1.0000000 3.665646e-08 1.832823e-08 [382,] 1.0000000 8.777207e-08 4.388603e-08 [383,] 0.9999999 2.130467e-07 1.065234e-07 [384,] 0.9999997 5.097429e-07 2.548714e-07 [385,] 0.9999995 9.650960e-07 4.825480e-07 [386,] 0.9999990 2.098654e-06 1.049327e-06 [387,] 0.9999976 4.823996e-06 2.411998e-06 [388,] 0.9999949 1.017128e-05 5.085642e-06 [389,] 0.9999932 1.369560e-05 6.847801e-06 [390,] 0.9999849 3.021501e-05 1.510750e-05 [391,] 0.9999750 5.006669e-05 2.503334e-05 [392,] 0.9999484 1.032830e-04 5.164149e-05 [393,] 0.9998980 2.040522e-04 1.020261e-04 [394,] 0.9997933 4.133326e-04 2.066663e-04 [395,] 0.9995787 8.426896e-04 4.213448e-04 [396,] 0.9991633 1.673403e-03 8.367016e-04 [397,] 0.9983938 3.212325e-03 1.606163e-03 [398,] 0.9970859 5.828218e-03 2.914109e-03 [399,] 0.9948777 1.024464e-02 5.122321e-03 [400,] 0.9909115 1.817702e-02 9.088511e-03 [401,] 0.9837571 3.248589e-02 1.624295e-02 [402,] 0.9719053 5.618932e-02 2.809466e-02 [403,] 0.9572517 8.549669e-02 4.274834e-02 [404,] 0.9325706 1.348587e-01 6.742936e-02 [405,] 0.8936672 2.126657e-01 1.063328e-01 [406,] 0.8372058 3.255884e-01 1.627942e-01 [407,] 0.7608735 4.782530e-01 2.391265e-01 [408,] 0.9735282 5.294361e-02 2.647180e-02 [409,] 0.9434052 1.131895e-01 5.659476e-02 [410,] 0.9051731 1.896538e-01 9.482689e-02 [411,] 0.8616411 2.767177e-01 1.383589e-01 [412,] 0.7500628 4.998745e-01 2.499372e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/198301291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/2kh231291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/3kh231291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/4kh231291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5kh231291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') > qqline(mysum$resid) > grid() > dev.off() null device 1 > (myerror <- as.ts(mysum$resid)) Time Series: Start = 1 End = 431 Frequency = 1 1 2 3 4 5 -277397.7292 2557837.5403 1524849.5589 -872756.8463 -260723.9355 6 7 8 9 10 499912.6898 50803.3342 233096.1568 546213.5869 765744.5259 11 12 13 14 15 118815.3978 415065.5112 363390.6666 -1106733.5215 81428.2065 16 17 18 19 20 598941.0692 -285286.3107 -223449.8448 222456.0823 -781630.3797 21 22 23 24 25 382247.6700 123835.2468 -94773.1470 -1025742.4420 -150822.8397 26 27 28 29 30 72550.9570 -272316.0688 -248027.7430 398664.0991 534699.5439 31 32 33 34 35 269029.6805 22432.4533 232558.4609 129733.1617 495152.7680 36 37 38 39 40 -677577.9129 156925.7774 91863.3885 -38270.7135 -342467.4789 41 42 43 44 45 -230150.3517 -128825.6111 186076.1321 -165671.9290 37071.3337 46 47 48 49 50 -86017.8086 231502.1696 -377790.7035 -1374.9665 -195568.0272 51 52 53 54 55 -30126.9217 -329273.4078 -340355.4879 -16799.7902 -506760.5780 56 57 58 59 60 34203.7748 180846.5192 59780.1536 23037.7797 -163328.4059 61 62 63 64 65 18400.9273 -113341.5416 9576.2456 -195793.0162 24724.8221 66 67 68 69 70 563636.9711 24732.4057 -46375.2413 -73770.5888 -204071.2810 71 72 73 74 75 -145607.0407 32454.6517 -61495.9718 18524.9487 -5115.8936 76 77 78 79 80 -106754.6017 19465.0654 18368.4564 115684.9807 195574.6939 81 82 83 84 85 -18171.0815 -76936.2201 36834.4012 -26100.2435 93365.5534 86 87 88 89 90 -137127.4478 -117531.4463 71990.4549 -37266.1893 82710.6706 91 92 93 94 95 29855.4842 -76301.2867 77576.4644 -60184.6443 -774.7735 96 97 98 99 100 -12507.9392 164297.9901 368028.3966 91437.6969 -42069.6757 101 102 103 104 105 -97365.9034 -43810.1211 -51991.9092 117123.1199 29944.0083 106 107 108 109 110 -97531.6477 -266971.0776 128424.8455 -91738.3138 -100438.5848 111 112 113 114 115 96922.0217 -254576.2834 -179297.6440 -118452.4114 -14628.1591 116 117 118 119 120 796.5017 88509.4743 123514.9171 -161398.0342 120703.3278 121 122 123 124 125 -196652.6095 31573.9512 -200582.8122 76243.4774 -70105.2464 126 127 128 129 130 -195592.4510 23944.4897 -549.6192 -5360.9040 23944.4897 131 132 133 134 135 15113.5220 130137.6496 23944.4897 -53630.8243 23944.4897 136 137 138 139 140 34400.4677 15071.0163 23944.4897 23944.4897 23944.4897 141 142 143 144 145 24850.1705 23944.4897 10516.6532 15264.9363 18116.0680 146 147 148 149 150 12172.3994 23944.4897 9954.9738 23944.4897 26956.4369 151 152 153 154 155 9263.7127 10444.8738 3095.7455 -2841.6848 23274.9586 156 157 158 159 160 23944.4897 23944.4897 21721.8609 33740.4414 23944.4897 161 162 163 164 165 23944.4897 23944.4897 23944.4897 67491.0308 47608.3635 166 167 168 169 170 23944.4897 23944.4897 55212.5717 7265.4943 23944.4897 171 172 173 174 175 23944.4897 20907.2606 23944.4897 11860.3540 23944.4897 176 177 178 179 180 18500.2461 23944.4897 -28757.0335 23944.4897 23944.4897 181 182 183 184 185 83867.6908 -27383.2329 93898.3554 -86816.9886 -34571.1742 186 187 188 189 190 -43410.3807 -4819.4832 14322.6612 -1823.4228 23944.4897 191 192 193 194 195 8647.1684 23944.4897 15946.0910 23944.4897 16616.4565 196 197 198 199 200 -31164.0924 -162580.2171 -31164.0924 -31164.0924 -69414.3865 201 202 203 204 205 -31164.0924 -31164.0924 -31164.0924 -24625.9356 -26498.6333 206 207 208 209 210 -45596.5627 -25412.0680 -31164.0924 -28390.8870 -31164.0924 211 212 213 214 215 23944.4897 23944.4897 56991.7930 16129.6807 17321.8019 216 217 218 219 220 15604.9920 20018.5381 65743.9018 15604.9920 -14771.6417 221 222 223 224 225 -7888.7766 23944.4897 24859.2736 -1074.0034 25037.9354 226 227 228 229 230 -386.2961 23944.4897 23944.4897 19715.0099 23944.4897 231 232 233 234 235 11162.2202 23944.4897 23944.4897 17857.3750 23944.4897 236 237 238 239 240 26439.4266 23944.4897 14159.2402 23944.4897 23944.4897 241 242 243 244 245 23944.4897 25092.1689 23944.4897 13831.5253 31776.3908 246 247 248 249 250 -31164.0924 -31164.0924 -16839.4265 -31164.0924 -31164.0924 251 252 253 254 255 -12268.5111 -31164.0924 18088.9535 36297.4716 1066.2229 256 257 258 259 260 23944.4897 24008.5743 -3218.5157 -31164.0924 23944.4897 261 262 263 264 265 -31164.0924 -31164.0924 24904.4705 21890.4878 -31164.0924 266 267 268 269 270 23944.4897 -31164.0924 80047.1915 -31164.0924 -55236.0503 271 272 273 274 275 -31164.0924 -31164.0924 -31164.0924 -28665.2822 23944.4897 276 277 278 279 280 57477.8478 22102.3950 23944.4897 57933.9111 56208.6106 281 282 283 284 285 -31164.0924 23944.4897 102833.7758 23944.4897 21869.4117 286 287 288 289 290 23944.4897 -31164.0924 24714.2721 -31164.0924 121394.6428 291 292 293 294 295 -31164.0924 44696.2917 41880.9308 23944.4897 11990.8023 296 297 298 299 300 13311.6374 23944.4897 -31164.0924 23944.4897 -13583.2499 301 302 303 304 305 -39224.8440 8133.8845 23944.4897 -31164.0924 126341.8884 306 307 308 309 310 -2952.0544 -31164.0924 -74616.6713 8317.1496 -31164.0924 311 312 313 314 315 -31164.0924 96925.5331 -31164.0924 -61722.3414 1172.4179 316 317 318 319 320 59855.7592 25344.0301 -28662.1095 18498.3189 23944.4897 321 322 323 324 325 23944.4897 23944.4897 23944.4897 -31164.0924 23944.4897 326 327 328 329 330 21871.6492 -31164.0924 -31164.0924 23944.4897 23944.4897 331 332 333 334 335 23944.4897 23944.4897 23944.4897 18436.9685 24916.6033 336 337 338 339 340 23944.4897 23283.8362 23944.4897 22470.2156 12543.2403 341 342 343 344 345 766.5812 23944.4897 23944.4897 23944.4897 23944.4897 346 347 348 349 350 23944.4897 12241.4031 23944.4897 23944.4897 23944.4897 351 352 353 354 355 40236.0807 19881.8725 16754.7417 22233.3625 23944.4897 356 357 358 359 360 23944.4897 94666.5242 23944.4897 23944.4897 23944.4897 361 362 363 364 365 17534.0406 23944.4897 23944.4897 -37653.4017 7735.2309 366 367 368 369 370 23944.4897 23944.4897 23944.4897 23944.4897 23944.4897 371 372 373 374 375 23944.4897 88363.2415 38934.9047 6811.2441 23944.4897 376 377 378 379 380 23944.4897 28078.4082 14034.6561 23944.4897 23944.4897 381 382 383 384 385 23944.4897 108817.1919 164863.0184 23944.4897 18419.9700 386 387 388 389 390 -97878.2616 27285.8956 10634.9196 16119.9896 8879.0068 391 392 393 394 395 15693.2771 -936.5007 -18866.5525 74608.9461 -99649.4056 396 397 398 399 400 -6594.9425 -82488.7376 -697508.8098 9888.3789 121170.3199 401 402 403 404 405 43695.7911 28536.0182 1668.2186 -88304.8133 1638.0231 406 407 408 409 410 14191.8476 34921.3543 -10054.1141 -12874.9790 -5984.2293 411 412 413 414 415 -82507.1098 -113565.0487 -75888.7341 -122992.5639 6515.1212 416 417 418 419 420 -2245.6856 160267.5784 -103626.4598 -18465.1137 -8329.4936 421 422 423 424 425 -40960.8889 -216291.7420 -497573.8832 -414112.1634 -109008.0932 426 427 428 429 430 -249580.5640 -311658.1463 -1284038.3256 -184370.2558 -134950.8637 431 -377355.7461 > postscript(file="/var/www/html/freestat/rcomp/tmp/6icqu1291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > dum <- cbind(lag(myerror,k=1),myerror) > dum Time Series: Start = 0 End = 431 Frequency = 1 lag(myerror, k = 1) myerror 0 -277397.7292 NA 1 2557837.5403 -277397.7292 2 1524849.5589 2557837.5403 3 -872756.8463 1524849.5589 4 -260723.9355 -872756.8463 5 499912.6898 -260723.9355 6 50803.3342 499912.6898 7 233096.1568 50803.3342 8 546213.5869 233096.1568 9 765744.5259 546213.5869 10 118815.3978 765744.5259 11 415065.5112 118815.3978 12 363390.6666 415065.5112 13 -1106733.5215 363390.6666 14 81428.2065 -1106733.5215 15 598941.0692 81428.2065 16 -285286.3107 598941.0692 17 -223449.8448 -285286.3107 18 222456.0823 -223449.8448 19 -781630.3797 222456.0823 20 382247.6700 -781630.3797 21 123835.2468 382247.6700 22 -94773.1470 123835.2468 23 -1025742.4420 -94773.1470 24 -150822.8397 -1025742.4420 25 72550.9570 -150822.8397 26 -272316.0688 72550.9570 27 -248027.7430 -272316.0688 28 398664.0991 -248027.7430 29 534699.5439 398664.0991 30 269029.6805 534699.5439 31 22432.4533 269029.6805 32 232558.4609 22432.4533 33 129733.1617 232558.4609 34 495152.7680 129733.1617 35 -677577.9129 495152.7680 36 156925.7774 -677577.9129 37 91863.3885 156925.7774 38 -38270.7135 91863.3885 39 -342467.4789 -38270.7135 40 -230150.3517 -342467.4789 41 -128825.6111 -230150.3517 42 186076.1321 -128825.6111 43 -165671.9290 186076.1321 44 37071.3337 -165671.9290 45 -86017.8086 37071.3337 46 231502.1696 -86017.8086 47 -377790.7035 231502.1696 48 -1374.9665 -377790.7035 49 -195568.0272 -1374.9665 50 -30126.9217 -195568.0272 51 -329273.4078 -30126.9217 52 -340355.4879 -329273.4078 53 -16799.7902 -340355.4879 54 -506760.5780 -16799.7902 55 34203.7748 -506760.5780 56 180846.5192 34203.7748 57 59780.1536 180846.5192 58 23037.7797 59780.1536 59 -163328.4059 23037.7797 60 18400.9273 -163328.4059 61 -113341.5416 18400.9273 62 9576.2456 -113341.5416 63 -195793.0162 9576.2456 64 24724.8221 -195793.0162 65 563636.9711 24724.8221 66 24732.4057 563636.9711 67 -46375.2413 24732.4057 68 -73770.5888 -46375.2413 69 -204071.2810 -73770.5888 70 -145607.0407 -204071.2810 71 32454.6517 -145607.0407 72 -61495.9718 32454.6517 73 18524.9487 -61495.9718 74 -5115.8936 18524.9487 75 -106754.6017 -5115.8936 76 19465.0654 -106754.6017 77 18368.4564 19465.0654 78 115684.9807 18368.4564 79 195574.6939 115684.9807 80 -18171.0815 195574.6939 81 -76936.2201 -18171.0815 82 36834.4012 -76936.2201 83 -26100.2435 36834.4012 84 93365.5534 -26100.2435 85 -137127.4478 93365.5534 86 -117531.4463 -137127.4478 87 71990.4549 -117531.4463 88 -37266.1893 71990.4549 89 82710.6706 -37266.1893 90 29855.4842 82710.6706 91 -76301.2867 29855.4842 92 77576.4644 -76301.2867 93 -60184.6443 77576.4644 94 -774.7735 -60184.6443 95 -12507.9392 -774.7735 96 164297.9901 -12507.9392 97 368028.3966 164297.9901 98 91437.6969 368028.3966 99 -42069.6757 91437.6969 100 -97365.9034 -42069.6757 101 -43810.1211 -97365.9034 102 -51991.9092 -43810.1211 103 117123.1199 -51991.9092 104 29944.0083 117123.1199 105 -97531.6477 29944.0083 106 -266971.0776 -97531.6477 107 128424.8455 -266971.0776 108 -91738.3138 128424.8455 109 -100438.5848 -91738.3138 110 96922.0217 -100438.5848 111 -254576.2834 96922.0217 112 -179297.6440 -254576.2834 113 -118452.4114 -179297.6440 114 -14628.1591 -118452.4114 115 796.5017 -14628.1591 116 88509.4743 796.5017 117 123514.9171 88509.4743 118 -161398.0342 123514.9171 119 120703.3278 -161398.0342 120 -196652.6095 120703.3278 121 31573.9512 -196652.6095 122 -200582.8122 31573.9512 123 76243.4774 -200582.8122 124 -70105.2464 76243.4774 125 -195592.4510 -70105.2464 126 23944.4897 -195592.4510 127 -549.6192 23944.4897 128 -5360.9040 -549.6192 129 23944.4897 -5360.9040 130 15113.5220 23944.4897 131 130137.6496 15113.5220 132 23944.4897 130137.6496 133 -53630.8243 23944.4897 134 23944.4897 -53630.8243 135 34400.4677 23944.4897 136 15071.0163 34400.4677 137 23944.4897 15071.0163 138 23944.4897 23944.4897 139 23944.4897 23944.4897 140 24850.1705 23944.4897 141 23944.4897 24850.1705 142 10516.6532 23944.4897 143 15264.9363 10516.6532 144 18116.0680 15264.9363 145 12172.3994 18116.0680 146 23944.4897 12172.3994 147 9954.9738 23944.4897 148 23944.4897 9954.9738 149 26956.4369 23944.4897 150 9263.7127 26956.4369 151 10444.8738 9263.7127 152 3095.7455 10444.8738 153 -2841.6848 3095.7455 154 23274.9586 -2841.6848 155 23944.4897 23274.9586 156 23944.4897 23944.4897 157 21721.8609 23944.4897 158 33740.4414 21721.8609 159 23944.4897 33740.4414 160 23944.4897 23944.4897 161 23944.4897 23944.4897 162 23944.4897 23944.4897 163 67491.0308 23944.4897 164 47608.3635 67491.0308 165 23944.4897 47608.3635 166 23944.4897 23944.4897 167 55212.5717 23944.4897 168 7265.4943 55212.5717 169 23944.4897 7265.4943 170 23944.4897 23944.4897 171 20907.2606 23944.4897 172 23944.4897 20907.2606 173 11860.3540 23944.4897 174 23944.4897 11860.3540 175 18500.2461 23944.4897 176 23944.4897 18500.2461 177 -28757.0335 23944.4897 178 23944.4897 -28757.0335 179 23944.4897 23944.4897 180 83867.6908 23944.4897 181 -27383.2329 83867.6908 182 93898.3554 -27383.2329 183 -86816.9886 93898.3554 184 -34571.1742 -86816.9886 185 -43410.3807 -34571.1742 186 -4819.4832 -43410.3807 187 14322.6612 -4819.4832 188 -1823.4228 14322.6612 189 23944.4897 -1823.4228 190 8647.1684 23944.4897 191 23944.4897 8647.1684 192 15946.0910 23944.4897 193 23944.4897 15946.0910 194 16616.4565 23944.4897 195 -31164.0924 16616.4565 196 -162580.2171 -31164.0924 197 -31164.0924 -162580.2171 198 -31164.0924 -31164.0924 199 -69414.3865 -31164.0924 200 -31164.0924 -69414.3865 201 -31164.0924 -31164.0924 202 -31164.0924 -31164.0924 203 -24625.9356 -31164.0924 204 -26498.6333 -24625.9356 205 -45596.5627 -26498.6333 206 -25412.0680 -45596.5627 207 -31164.0924 -25412.0680 208 -28390.8870 -31164.0924 209 -31164.0924 -28390.8870 210 23944.4897 -31164.0924 211 23944.4897 23944.4897 212 56991.7930 23944.4897 213 16129.6807 56991.7930 214 17321.8019 16129.6807 215 15604.9920 17321.8019 216 20018.5381 15604.9920 217 65743.9018 20018.5381 218 15604.9920 65743.9018 219 -14771.6417 15604.9920 220 -7888.7766 -14771.6417 221 23944.4897 -7888.7766 222 24859.2736 23944.4897 223 -1074.0034 24859.2736 224 25037.9354 -1074.0034 225 -386.2961 25037.9354 226 23944.4897 -386.2961 227 23944.4897 23944.4897 228 19715.0099 23944.4897 229 23944.4897 19715.0099 230 11162.2202 23944.4897 231 23944.4897 11162.2202 232 23944.4897 23944.4897 233 17857.3750 23944.4897 234 23944.4897 17857.3750 235 26439.4266 23944.4897 236 23944.4897 26439.4266 237 14159.2402 23944.4897 238 23944.4897 14159.2402 239 23944.4897 23944.4897 240 23944.4897 23944.4897 241 25092.1689 23944.4897 242 23944.4897 25092.1689 243 13831.5253 23944.4897 244 31776.3908 13831.5253 245 -31164.0924 31776.3908 246 -31164.0924 -31164.0924 247 -16839.4265 -31164.0924 248 -31164.0924 -16839.4265 249 -31164.0924 -31164.0924 250 -12268.5111 -31164.0924 251 -31164.0924 -12268.5111 252 18088.9535 -31164.0924 253 36297.4716 18088.9535 254 1066.2229 36297.4716 255 23944.4897 1066.2229 256 24008.5743 23944.4897 257 -3218.5157 24008.5743 258 -31164.0924 -3218.5157 259 23944.4897 -31164.0924 260 -31164.0924 23944.4897 261 -31164.0924 -31164.0924 262 24904.4705 -31164.0924 263 21890.4878 24904.4705 264 -31164.0924 21890.4878 265 23944.4897 -31164.0924 266 -31164.0924 23944.4897 267 80047.1915 -31164.0924 268 -31164.0924 80047.1915 269 -55236.0503 -31164.0924 270 -31164.0924 -55236.0503 271 -31164.0924 -31164.0924 272 -31164.0924 -31164.0924 273 -28665.2822 -31164.0924 274 23944.4897 -28665.2822 275 57477.8478 23944.4897 276 22102.3950 57477.8478 277 23944.4897 22102.3950 278 57933.9111 23944.4897 279 56208.6106 57933.9111 280 -31164.0924 56208.6106 281 23944.4897 -31164.0924 282 102833.7758 23944.4897 283 23944.4897 102833.7758 284 21869.4117 23944.4897 285 23944.4897 21869.4117 286 -31164.0924 23944.4897 287 24714.2721 -31164.0924 288 -31164.0924 24714.2721 289 121394.6428 -31164.0924 290 -31164.0924 121394.6428 291 44696.2917 -31164.0924 292 41880.9308 44696.2917 293 23944.4897 41880.9308 294 11990.8023 23944.4897 295 13311.6374 11990.8023 296 23944.4897 13311.6374 297 -31164.0924 23944.4897 298 23944.4897 -31164.0924 299 -13583.2499 23944.4897 300 -39224.8440 -13583.2499 301 8133.8845 -39224.8440 302 23944.4897 8133.8845 303 -31164.0924 23944.4897 304 126341.8884 -31164.0924 305 -2952.0544 126341.8884 306 -31164.0924 -2952.0544 307 -74616.6713 -31164.0924 308 8317.1496 -74616.6713 309 -31164.0924 8317.1496 310 -31164.0924 -31164.0924 311 96925.5331 -31164.0924 312 -31164.0924 96925.5331 313 -61722.3414 -31164.0924 314 1172.4179 -61722.3414 315 59855.7592 1172.4179 316 25344.0301 59855.7592 317 -28662.1095 25344.0301 318 18498.3189 -28662.1095 319 23944.4897 18498.3189 320 23944.4897 23944.4897 321 23944.4897 23944.4897 322 23944.4897 23944.4897 323 -31164.0924 23944.4897 324 23944.4897 -31164.0924 325 21871.6492 23944.4897 326 -31164.0924 21871.6492 327 -31164.0924 -31164.0924 328 23944.4897 -31164.0924 329 23944.4897 23944.4897 330 23944.4897 23944.4897 331 23944.4897 23944.4897 332 23944.4897 23944.4897 333 18436.9685 23944.4897 334 24916.6033 18436.9685 335 23944.4897 24916.6033 336 23283.8362 23944.4897 337 23944.4897 23283.8362 338 22470.2156 23944.4897 339 12543.2403 22470.2156 340 766.5812 12543.2403 341 23944.4897 766.5812 342 23944.4897 23944.4897 343 23944.4897 23944.4897 344 23944.4897 23944.4897 345 23944.4897 23944.4897 346 12241.4031 23944.4897 347 23944.4897 12241.4031 348 23944.4897 23944.4897 349 23944.4897 23944.4897 350 40236.0807 23944.4897 351 19881.8725 40236.0807 352 16754.7417 19881.8725 353 22233.3625 16754.7417 354 23944.4897 22233.3625 355 23944.4897 23944.4897 356 94666.5242 23944.4897 357 23944.4897 94666.5242 358 23944.4897 23944.4897 359 23944.4897 23944.4897 360 17534.0406 23944.4897 361 23944.4897 17534.0406 362 23944.4897 23944.4897 363 -37653.4017 23944.4897 364 7735.2309 -37653.4017 365 23944.4897 7735.2309 366 23944.4897 23944.4897 367 23944.4897 23944.4897 368 23944.4897 23944.4897 369 23944.4897 23944.4897 370 23944.4897 23944.4897 371 88363.2415 23944.4897 372 38934.9047 88363.2415 373 6811.2441 38934.9047 374 23944.4897 6811.2441 375 23944.4897 23944.4897 376 28078.4082 23944.4897 377 14034.6561 28078.4082 378 23944.4897 14034.6561 379 23944.4897 23944.4897 380 23944.4897 23944.4897 381 108817.1919 23944.4897 382 164863.0184 108817.1919 383 23944.4897 164863.0184 384 18419.9700 23944.4897 385 -97878.2616 18419.9700 386 27285.8956 -97878.2616 387 10634.9196 27285.8956 388 16119.9896 10634.9196 389 8879.0068 16119.9896 390 15693.2771 8879.0068 391 -936.5007 15693.2771 392 -18866.5525 -936.5007 393 74608.9461 -18866.5525 394 -99649.4056 74608.9461 395 -6594.9425 -99649.4056 396 -82488.7376 -6594.9425 397 -697508.8098 -82488.7376 398 9888.3789 -697508.8098 399 121170.3199 9888.3789 400 43695.7911 121170.3199 401 28536.0182 43695.7911 402 1668.2186 28536.0182 403 -88304.8133 1668.2186 404 1638.0231 -88304.8133 405 14191.8476 1638.0231 406 34921.3543 14191.8476 407 -10054.1141 34921.3543 408 -12874.9790 -10054.1141 409 -5984.2293 -12874.9790 410 -82507.1098 -5984.2293 411 -113565.0487 -82507.1098 412 -75888.7341 -113565.0487 413 -122992.5639 -75888.7341 414 6515.1212 -122992.5639 415 -2245.6856 6515.1212 416 160267.5784 -2245.6856 417 -103626.4598 160267.5784 418 -18465.1137 -103626.4598 419 -8329.4936 -18465.1137 420 -40960.8889 -8329.4936 421 -216291.7420 -40960.8889 422 -497573.8832 -216291.7420 423 -414112.1634 -497573.8832 424 -109008.0932 -414112.1634 425 -249580.5640 -109008.0932 426 -311658.1463 -249580.5640 427 -1284038.3256 -311658.1463 428 -184370.2558 -1284038.3256 429 -134950.8637 -184370.2558 430 -377355.7461 -134950.8637 431 NA -377355.7461 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2557837.5403 -277397.7292 [2,] 1524849.5589 2557837.5403 [3,] -872756.8463 1524849.5589 [4,] -260723.9355 -872756.8463 [5,] 499912.6898 -260723.9355 [6,] 50803.3342 499912.6898 [7,] 233096.1568 50803.3342 [8,] 546213.5869 233096.1568 [9,] 765744.5259 546213.5869 [10,] 118815.3978 765744.5259 [11,] 415065.5112 118815.3978 [12,] 363390.6666 415065.5112 [13,] -1106733.5215 363390.6666 [14,] 81428.2065 -1106733.5215 [15,] 598941.0692 81428.2065 [16,] -285286.3107 598941.0692 [17,] -223449.8448 -285286.3107 [18,] 222456.0823 -223449.8448 [19,] -781630.3797 222456.0823 [20,] 382247.6700 -781630.3797 [21,] 123835.2468 382247.6700 [22,] -94773.1470 123835.2468 [23,] -1025742.4420 -94773.1470 [24,] -150822.8397 -1025742.4420 [25,] 72550.9570 -150822.8397 [26,] -272316.0688 72550.9570 [27,] -248027.7430 -272316.0688 [28,] 398664.0991 -248027.7430 [29,] 534699.5439 398664.0991 [30,] 269029.6805 534699.5439 [31,] 22432.4533 269029.6805 [32,] 232558.4609 22432.4533 [33,] 129733.1617 232558.4609 [34,] 495152.7680 129733.1617 [35,] -677577.9129 495152.7680 [36,] 156925.7774 -677577.9129 [37,] 91863.3885 156925.7774 [38,] -38270.7135 91863.3885 [39,] -342467.4789 -38270.7135 [40,] -230150.3517 -342467.4789 [41,] -128825.6111 -230150.3517 [42,] 186076.1321 -128825.6111 [43,] -165671.9290 186076.1321 [44,] 37071.3337 -165671.9290 [45,] -86017.8086 37071.3337 [46,] 231502.1696 -86017.8086 [47,] -377790.7035 231502.1696 [48,] -1374.9665 -377790.7035 [49,] -195568.0272 -1374.9665 [50,] -30126.9217 -195568.0272 [51,] -329273.4078 -30126.9217 [52,] -340355.4879 -329273.4078 [53,] -16799.7902 -340355.4879 [54,] -506760.5780 -16799.7902 [55,] 34203.7748 -506760.5780 [56,] 180846.5192 34203.7748 [57,] 59780.1536 180846.5192 [58,] 23037.7797 59780.1536 [59,] -163328.4059 23037.7797 [60,] 18400.9273 -163328.4059 [61,] -113341.5416 18400.9273 [62,] 9576.2456 -113341.5416 [63,] -195793.0162 9576.2456 [64,] 24724.8221 -195793.0162 [65,] 563636.9711 24724.8221 [66,] 24732.4057 563636.9711 [67,] -46375.2413 24732.4057 [68,] -73770.5888 -46375.2413 [69,] -204071.2810 -73770.5888 [70,] -145607.0407 -204071.2810 [71,] 32454.6517 -145607.0407 [72,] -61495.9718 32454.6517 [73,] 18524.9487 -61495.9718 [74,] -5115.8936 18524.9487 [75,] -106754.6017 -5115.8936 [76,] 19465.0654 -106754.6017 [77,] 18368.4564 19465.0654 [78,] 115684.9807 18368.4564 [79,] 195574.6939 115684.9807 [80,] -18171.0815 195574.6939 [81,] -76936.2201 -18171.0815 [82,] 36834.4012 -76936.2201 [83,] -26100.2435 36834.4012 [84,] 93365.5534 -26100.2435 [85,] -137127.4478 93365.5534 [86,] -117531.4463 -137127.4478 [87,] 71990.4549 -117531.4463 [88,] -37266.1893 71990.4549 [89,] 82710.6706 -37266.1893 [90,] 29855.4842 82710.6706 [91,] -76301.2867 29855.4842 [92,] 77576.4644 -76301.2867 [93,] -60184.6443 77576.4644 [94,] -774.7735 -60184.6443 [95,] -12507.9392 -774.7735 [96,] 164297.9901 -12507.9392 [97,] 368028.3966 164297.9901 [98,] 91437.6969 368028.3966 [99,] -42069.6757 91437.6969 [100,] -97365.9034 -42069.6757 [101,] -43810.1211 -97365.9034 [102,] -51991.9092 -43810.1211 [103,] 117123.1199 -51991.9092 [104,] 29944.0083 117123.1199 [105,] -97531.6477 29944.0083 [106,] -266971.0776 -97531.6477 [107,] 128424.8455 -266971.0776 [108,] -91738.3138 128424.8455 [109,] -100438.5848 -91738.3138 [110,] 96922.0217 -100438.5848 [111,] -254576.2834 96922.0217 [112,] -179297.6440 -254576.2834 [113,] -118452.4114 -179297.6440 [114,] -14628.1591 -118452.4114 [115,] 796.5017 -14628.1591 [116,] 88509.4743 796.5017 [117,] 123514.9171 88509.4743 [118,] -161398.0342 123514.9171 [119,] 120703.3278 -161398.0342 [120,] -196652.6095 120703.3278 [121,] 31573.9512 -196652.6095 [122,] -200582.8122 31573.9512 [123,] 76243.4774 -200582.8122 [124,] -70105.2464 76243.4774 [125,] -195592.4510 -70105.2464 [126,] 23944.4897 -195592.4510 [127,] -549.6192 23944.4897 [128,] -5360.9040 -549.6192 [129,] 23944.4897 -5360.9040 [130,] 15113.5220 23944.4897 [131,] 130137.6496 15113.5220 [132,] 23944.4897 130137.6496 [133,] -53630.8243 23944.4897 [134,] 23944.4897 -53630.8243 [135,] 34400.4677 23944.4897 [136,] 15071.0163 34400.4677 [137,] 23944.4897 15071.0163 [138,] 23944.4897 23944.4897 [139,] 23944.4897 23944.4897 [140,] 24850.1705 23944.4897 [141,] 23944.4897 24850.1705 [142,] 10516.6532 23944.4897 [143,] 15264.9363 10516.6532 [144,] 18116.0680 15264.9363 [145,] 12172.3994 18116.0680 [146,] 23944.4897 12172.3994 [147,] 9954.9738 23944.4897 [148,] 23944.4897 9954.9738 [149,] 26956.4369 23944.4897 [150,] 9263.7127 26956.4369 [151,] 10444.8738 9263.7127 [152,] 3095.7455 10444.8738 [153,] -2841.6848 3095.7455 [154,] 23274.9586 -2841.6848 [155,] 23944.4897 23274.9586 [156,] 23944.4897 23944.4897 [157,] 21721.8609 23944.4897 [158,] 33740.4414 21721.8609 [159,] 23944.4897 33740.4414 [160,] 23944.4897 23944.4897 [161,] 23944.4897 23944.4897 [162,] 23944.4897 23944.4897 [163,] 67491.0308 23944.4897 [164,] 47608.3635 67491.0308 [165,] 23944.4897 47608.3635 [166,] 23944.4897 23944.4897 [167,] 55212.5717 23944.4897 [168,] 7265.4943 55212.5717 [169,] 23944.4897 7265.4943 [170,] 23944.4897 23944.4897 [171,] 20907.2606 23944.4897 [172,] 23944.4897 20907.2606 [173,] 11860.3540 23944.4897 [174,] 23944.4897 11860.3540 [175,] 18500.2461 23944.4897 [176,] 23944.4897 18500.2461 [177,] -28757.0335 23944.4897 [178,] 23944.4897 -28757.0335 [179,] 23944.4897 23944.4897 [180,] 83867.6908 23944.4897 [181,] -27383.2329 83867.6908 [182,] 93898.3554 -27383.2329 [183,] -86816.9886 93898.3554 [184,] -34571.1742 -86816.9886 [185,] -43410.3807 -34571.1742 [186,] -4819.4832 -43410.3807 [187,] 14322.6612 -4819.4832 [188,] -1823.4228 14322.6612 [189,] 23944.4897 -1823.4228 [190,] 8647.1684 23944.4897 [191,] 23944.4897 8647.1684 [192,] 15946.0910 23944.4897 [193,] 23944.4897 15946.0910 [194,] 16616.4565 23944.4897 [195,] -31164.0924 16616.4565 [196,] -162580.2171 -31164.0924 [197,] -31164.0924 -162580.2171 [198,] -31164.0924 -31164.0924 [199,] -69414.3865 -31164.0924 [200,] -31164.0924 -69414.3865 [201,] -31164.0924 -31164.0924 [202,] -31164.0924 -31164.0924 [203,] -24625.9356 -31164.0924 [204,] -26498.6333 -24625.9356 [205,] -45596.5627 -26498.6333 [206,] -25412.0680 -45596.5627 [207,] -31164.0924 -25412.0680 [208,] -28390.8870 -31164.0924 [209,] -31164.0924 -28390.8870 [210,] 23944.4897 -31164.0924 [211,] 23944.4897 23944.4897 [212,] 56991.7930 23944.4897 [213,] 16129.6807 56991.7930 [214,] 17321.8019 16129.6807 [215,] 15604.9920 17321.8019 [216,] 20018.5381 15604.9920 [217,] 65743.9018 20018.5381 [218,] 15604.9920 65743.9018 [219,] -14771.6417 15604.9920 [220,] -7888.7766 -14771.6417 [221,] 23944.4897 -7888.7766 [222,] 24859.2736 23944.4897 [223,] -1074.0034 24859.2736 [224,] 25037.9354 -1074.0034 [225,] -386.2961 25037.9354 [226,] 23944.4897 -386.2961 [227,] 23944.4897 23944.4897 [228,] 19715.0099 23944.4897 [229,] 23944.4897 19715.0099 [230,] 11162.2202 23944.4897 [231,] 23944.4897 11162.2202 [232,] 23944.4897 23944.4897 [233,] 17857.3750 23944.4897 [234,] 23944.4897 17857.3750 [235,] 26439.4266 23944.4897 [236,] 23944.4897 26439.4266 [237,] 14159.2402 23944.4897 [238,] 23944.4897 14159.2402 [239,] 23944.4897 23944.4897 [240,] 23944.4897 23944.4897 [241,] 25092.1689 23944.4897 [242,] 23944.4897 25092.1689 [243,] 13831.5253 23944.4897 [244,] 31776.3908 13831.5253 [245,] -31164.0924 31776.3908 [246,] -31164.0924 -31164.0924 [247,] -16839.4265 -31164.0924 [248,] -31164.0924 -16839.4265 [249,] -31164.0924 -31164.0924 [250,] -12268.5111 -31164.0924 [251,] -31164.0924 -12268.5111 [252,] 18088.9535 -31164.0924 [253,] 36297.4716 18088.9535 [254,] 1066.2229 36297.4716 [255,] 23944.4897 1066.2229 [256,] 24008.5743 23944.4897 [257,] -3218.5157 24008.5743 [258,] -31164.0924 -3218.5157 [259,] 23944.4897 -31164.0924 [260,] -31164.0924 23944.4897 [261,] -31164.0924 -31164.0924 [262,] 24904.4705 -31164.0924 [263,] 21890.4878 24904.4705 [264,] -31164.0924 21890.4878 [265,] 23944.4897 -31164.0924 [266,] -31164.0924 23944.4897 [267,] 80047.1915 -31164.0924 [268,] -31164.0924 80047.1915 [269,] -55236.0503 -31164.0924 [270,] -31164.0924 -55236.0503 [271,] -31164.0924 -31164.0924 [272,] -31164.0924 -31164.0924 [273,] -28665.2822 -31164.0924 [274,] 23944.4897 -28665.2822 [275,] 57477.8478 23944.4897 [276,] 22102.3950 57477.8478 [277,] 23944.4897 22102.3950 [278,] 57933.9111 23944.4897 [279,] 56208.6106 57933.9111 [280,] -31164.0924 56208.6106 [281,] 23944.4897 -31164.0924 [282,] 102833.7758 23944.4897 [283,] 23944.4897 102833.7758 [284,] 21869.4117 23944.4897 [285,] 23944.4897 21869.4117 [286,] -31164.0924 23944.4897 [287,] 24714.2721 -31164.0924 [288,] -31164.0924 24714.2721 [289,] 121394.6428 -31164.0924 [290,] -31164.0924 121394.6428 [291,] 44696.2917 -31164.0924 [292,] 41880.9308 44696.2917 [293,] 23944.4897 41880.9308 [294,] 11990.8023 23944.4897 [295,] 13311.6374 11990.8023 [296,] 23944.4897 13311.6374 [297,] -31164.0924 23944.4897 [298,] 23944.4897 -31164.0924 [299,] -13583.2499 23944.4897 [300,] -39224.8440 -13583.2499 [301,] 8133.8845 -39224.8440 [302,] 23944.4897 8133.8845 [303,] -31164.0924 23944.4897 [304,] 126341.8884 -31164.0924 [305,] -2952.0544 126341.8884 [306,] -31164.0924 -2952.0544 [307,] -74616.6713 -31164.0924 [308,] 8317.1496 -74616.6713 [309,] -31164.0924 8317.1496 [310,] -31164.0924 -31164.0924 [311,] 96925.5331 -31164.0924 [312,] -31164.0924 96925.5331 [313,] -61722.3414 -31164.0924 [314,] 1172.4179 -61722.3414 [315,] 59855.7592 1172.4179 [316,] 25344.0301 59855.7592 [317,] -28662.1095 25344.0301 [318,] 18498.3189 -28662.1095 [319,] 23944.4897 18498.3189 [320,] 23944.4897 23944.4897 [321,] 23944.4897 23944.4897 [322,] 23944.4897 23944.4897 [323,] -31164.0924 23944.4897 [324,] 23944.4897 -31164.0924 [325,] 21871.6492 23944.4897 [326,] -31164.0924 21871.6492 [327,] -31164.0924 -31164.0924 [328,] 23944.4897 -31164.0924 [329,] 23944.4897 23944.4897 [330,] 23944.4897 23944.4897 [331,] 23944.4897 23944.4897 [332,] 23944.4897 23944.4897 [333,] 18436.9685 23944.4897 [334,] 24916.6033 18436.9685 [335,] 23944.4897 24916.6033 [336,] 23283.8362 23944.4897 [337,] 23944.4897 23283.8362 [338,] 22470.2156 23944.4897 [339,] 12543.2403 22470.2156 [340,] 766.5812 12543.2403 [341,] 23944.4897 766.5812 [342,] 23944.4897 23944.4897 [343,] 23944.4897 23944.4897 [344,] 23944.4897 23944.4897 [345,] 23944.4897 23944.4897 [346,] 12241.4031 23944.4897 [347,] 23944.4897 12241.4031 [348,] 23944.4897 23944.4897 [349,] 23944.4897 23944.4897 [350,] 40236.0807 23944.4897 [351,] 19881.8725 40236.0807 [352,] 16754.7417 19881.8725 [353,] 22233.3625 16754.7417 [354,] 23944.4897 22233.3625 [355,] 23944.4897 23944.4897 [356,] 94666.5242 23944.4897 [357,] 23944.4897 94666.5242 [358,] 23944.4897 23944.4897 [359,] 23944.4897 23944.4897 [360,] 17534.0406 23944.4897 [361,] 23944.4897 17534.0406 [362,] 23944.4897 23944.4897 [363,] -37653.4017 23944.4897 [364,] 7735.2309 -37653.4017 [365,] 23944.4897 7735.2309 [366,] 23944.4897 23944.4897 [367,] 23944.4897 23944.4897 [368,] 23944.4897 23944.4897 [369,] 23944.4897 23944.4897 [370,] 23944.4897 23944.4897 [371,] 88363.2415 23944.4897 [372,] 38934.9047 88363.2415 [373,] 6811.2441 38934.9047 [374,] 23944.4897 6811.2441 [375,] 23944.4897 23944.4897 [376,] 28078.4082 23944.4897 [377,] 14034.6561 28078.4082 [378,] 23944.4897 14034.6561 [379,] 23944.4897 23944.4897 [380,] 23944.4897 23944.4897 [381,] 108817.1919 23944.4897 [382,] 164863.0184 108817.1919 [383,] 23944.4897 164863.0184 [384,] 18419.9700 23944.4897 [385,] -97878.2616 18419.9700 [386,] 27285.8956 -97878.2616 [387,] 10634.9196 27285.8956 [388,] 16119.9896 10634.9196 [389,] 8879.0068 16119.9896 [390,] 15693.2771 8879.0068 [391,] -936.5007 15693.2771 [392,] -18866.5525 -936.5007 [393,] 74608.9461 -18866.5525 [394,] -99649.4056 74608.9461 [395,] -6594.9425 -99649.4056 [396,] -82488.7376 -6594.9425 [397,] -697508.8098 -82488.7376 [398,] 9888.3789 -697508.8098 [399,] 121170.3199 9888.3789 [400,] 43695.7911 121170.3199 [401,] 28536.0182 43695.7911 [402,] 1668.2186 28536.0182 [403,] -88304.8133 1668.2186 [404,] 1638.0231 -88304.8133 [405,] 14191.8476 1638.0231 [406,] 34921.3543 14191.8476 [407,] -10054.1141 34921.3543 [408,] -12874.9790 -10054.1141 [409,] -5984.2293 -12874.9790 [410,] -82507.1098 -5984.2293 [411,] -113565.0487 -82507.1098 [412,] -75888.7341 -113565.0487 [413,] -122992.5639 -75888.7341 [414,] 6515.1212 -122992.5639 [415,] -2245.6856 6515.1212 [416,] 160267.5784 -2245.6856 [417,] -103626.4598 160267.5784 [418,] -18465.1137 -103626.4598 [419,] -8329.4936 -18465.1137 [420,] -40960.8889 -8329.4936 [421,] -216291.7420 -40960.8889 [422,] -497573.8832 -216291.7420 [423,] -414112.1634 -497573.8832 [424,] -109008.0932 -414112.1634 [425,] -249580.5640 -109008.0932 [426,] -311658.1463 -249580.5640 [427,] -1284038.3256 -311658.1463 [428,] -184370.2558 -1284038.3256 [429,] -134950.8637 -184370.2558 [430,] -377355.7461 -134950.8637 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2557837.5403 -277397.7292 2 1524849.5589 2557837.5403 3 -872756.8463 1524849.5589 4 -260723.9355 -872756.8463 5 499912.6898 -260723.9355 6 50803.3342 499912.6898 7 233096.1568 50803.3342 8 546213.5869 233096.1568 9 765744.5259 546213.5869 10 118815.3978 765744.5259 11 415065.5112 118815.3978 12 363390.6666 415065.5112 13 -1106733.5215 363390.6666 14 81428.2065 -1106733.5215 15 598941.0692 81428.2065 16 -285286.3107 598941.0692 17 -223449.8448 -285286.3107 18 222456.0823 -223449.8448 19 -781630.3797 222456.0823 20 382247.6700 -781630.3797 21 123835.2468 382247.6700 22 -94773.1470 123835.2468 23 -1025742.4420 -94773.1470 24 -150822.8397 -1025742.4420 25 72550.9570 -150822.8397 26 -272316.0688 72550.9570 27 -248027.7430 -272316.0688 28 398664.0991 -248027.7430 29 534699.5439 398664.0991 30 269029.6805 534699.5439 31 22432.4533 269029.6805 32 232558.4609 22432.4533 33 129733.1617 232558.4609 34 495152.7680 129733.1617 35 -677577.9129 495152.7680 36 156925.7774 -677577.9129 37 91863.3885 156925.7774 38 -38270.7135 91863.3885 39 -342467.4789 -38270.7135 40 -230150.3517 -342467.4789 41 -128825.6111 -230150.3517 42 186076.1321 -128825.6111 43 -165671.9290 186076.1321 44 37071.3337 -165671.9290 45 -86017.8086 37071.3337 46 231502.1696 -86017.8086 47 -377790.7035 231502.1696 48 -1374.9665 -377790.7035 49 -195568.0272 -1374.9665 50 -30126.9217 -195568.0272 51 -329273.4078 -30126.9217 52 -340355.4879 -329273.4078 53 -16799.7902 -340355.4879 54 -506760.5780 -16799.7902 55 34203.7748 -506760.5780 56 180846.5192 34203.7748 57 59780.1536 180846.5192 58 23037.7797 59780.1536 59 -163328.4059 23037.7797 60 18400.9273 -163328.4059 61 -113341.5416 18400.9273 62 9576.2456 -113341.5416 63 -195793.0162 9576.2456 64 24724.8221 -195793.0162 65 563636.9711 24724.8221 66 24732.4057 563636.9711 67 -46375.2413 24732.4057 68 -73770.5888 -46375.2413 69 -204071.2810 -73770.5888 70 -145607.0407 -204071.2810 71 32454.6517 -145607.0407 72 -61495.9718 32454.6517 73 18524.9487 -61495.9718 74 -5115.8936 18524.9487 75 -106754.6017 -5115.8936 76 19465.0654 -106754.6017 77 18368.4564 19465.0654 78 115684.9807 18368.4564 79 195574.6939 115684.9807 80 -18171.0815 195574.6939 81 -76936.2201 -18171.0815 82 36834.4012 -76936.2201 83 -26100.2435 36834.4012 84 93365.5534 -26100.2435 85 -137127.4478 93365.5534 86 -117531.4463 -137127.4478 87 71990.4549 -117531.4463 88 -37266.1893 71990.4549 89 82710.6706 -37266.1893 90 29855.4842 82710.6706 91 -76301.2867 29855.4842 92 77576.4644 -76301.2867 93 -60184.6443 77576.4644 94 -774.7735 -60184.6443 95 -12507.9392 -774.7735 96 164297.9901 -12507.9392 97 368028.3966 164297.9901 98 91437.6969 368028.3966 99 -42069.6757 91437.6969 100 -97365.9034 -42069.6757 101 -43810.1211 -97365.9034 102 -51991.9092 -43810.1211 103 117123.1199 -51991.9092 104 29944.0083 117123.1199 105 -97531.6477 29944.0083 106 -266971.0776 -97531.6477 107 128424.8455 -266971.0776 108 -91738.3138 128424.8455 109 -100438.5848 -91738.3138 110 96922.0217 -100438.5848 111 -254576.2834 96922.0217 112 -179297.6440 -254576.2834 113 -118452.4114 -179297.6440 114 -14628.1591 -118452.4114 115 796.5017 -14628.1591 116 88509.4743 796.5017 117 123514.9171 88509.4743 118 -161398.0342 123514.9171 119 120703.3278 -161398.0342 120 -196652.6095 120703.3278 121 31573.9512 -196652.6095 122 -200582.8122 31573.9512 123 76243.4774 -200582.8122 124 -70105.2464 76243.4774 125 -195592.4510 -70105.2464 126 23944.4897 -195592.4510 127 -549.6192 23944.4897 128 -5360.9040 -549.6192 129 23944.4897 -5360.9040 130 15113.5220 23944.4897 131 130137.6496 15113.5220 132 23944.4897 130137.6496 133 -53630.8243 23944.4897 134 23944.4897 -53630.8243 135 34400.4677 23944.4897 136 15071.0163 34400.4677 137 23944.4897 15071.0163 138 23944.4897 23944.4897 139 23944.4897 23944.4897 140 24850.1705 23944.4897 141 23944.4897 24850.1705 142 10516.6532 23944.4897 143 15264.9363 10516.6532 144 18116.0680 15264.9363 145 12172.3994 18116.0680 146 23944.4897 12172.3994 147 9954.9738 23944.4897 148 23944.4897 9954.9738 149 26956.4369 23944.4897 150 9263.7127 26956.4369 151 10444.8738 9263.7127 152 3095.7455 10444.8738 153 -2841.6848 3095.7455 154 23274.9586 -2841.6848 155 23944.4897 23274.9586 156 23944.4897 23944.4897 157 21721.8609 23944.4897 158 33740.4414 21721.8609 159 23944.4897 33740.4414 160 23944.4897 23944.4897 161 23944.4897 23944.4897 162 23944.4897 23944.4897 163 67491.0308 23944.4897 164 47608.3635 67491.0308 165 23944.4897 47608.3635 166 23944.4897 23944.4897 167 55212.5717 23944.4897 168 7265.4943 55212.5717 169 23944.4897 7265.4943 170 23944.4897 23944.4897 171 20907.2606 23944.4897 172 23944.4897 20907.2606 173 11860.3540 23944.4897 174 23944.4897 11860.3540 175 18500.2461 23944.4897 176 23944.4897 18500.2461 177 -28757.0335 23944.4897 178 23944.4897 -28757.0335 179 23944.4897 23944.4897 180 83867.6908 23944.4897 181 -27383.2329 83867.6908 182 93898.3554 -27383.2329 183 -86816.9886 93898.3554 184 -34571.1742 -86816.9886 185 -43410.3807 -34571.1742 186 -4819.4832 -43410.3807 187 14322.6612 -4819.4832 188 -1823.4228 14322.6612 189 23944.4897 -1823.4228 190 8647.1684 23944.4897 191 23944.4897 8647.1684 192 15946.0910 23944.4897 193 23944.4897 15946.0910 194 16616.4565 23944.4897 195 -31164.0924 16616.4565 196 -162580.2171 -31164.0924 197 -31164.0924 -162580.2171 198 -31164.0924 -31164.0924 199 -69414.3865 -31164.0924 200 -31164.0924 -69414.3865 201 -31164.0924 -31164.0924 202 -31164.0924 -31164.0924 203 -24625.9356 -31164.0924 204 -26498.6333 -24625.9356 205 -45596.5627 -26498.6333 206 -25412.0680 -45596.5627 207 -31164.0924 -25412.0680 208 -28390.8870 -31164.0924 209 -31164.0924 -28390.8870 210 23944.4897 -31164.0924 211 23944.4897 23944.4897 212 56991.7930 23944.4897 213 16129.6807 56991.7930 214 17321.8019 16129.6807 215 15604.9920 17321.8019 216 20018.5381 15604.9920 217 65743.9018 20018.5381 218 15604.9920 65743.9018 219 -14771.6417 15604.9920 220 -7888.7766 -14771.6417 221 23944.4897 -7888.7766 222 24859.2736 23944.4897 223 -1074.0034 24859.2736 224 25037.9354 -1074.0034 225 -386.2961 25037.9354 226 23944.4897 -386.2961 227 23944.4897 23944.4897 228 19715.0099 23944.4897 229 23944.4897 19715.0099 230 11162.2202 23944.4897 231 23944.4897 11162.2202 232 23944.4897 23944.4897 233 17857.3750 23944.4897 234 23944.4897 17857.3750 235 26439.4266 23944.4897 236 23944.4897 26439.4266 237 14159.2402 23944.4897 238 23944.4897 14159.2402 239 23944.4897 23944.4897 240 23944.4897 23944.4897 241 25092.1689 23944.4897 242 23944.4897 25092.1689 243 13831.5253 23944.4897 244 31776.3908 13831.5253 245 -31164.0924 31776.3908 246 -31164.0924 -31164.0924 247 -16839.4265 -31164.0924 248 -31164.0924 -16839.4265 249 -31164.0924 -31164.0924 250 -12268.5111 -31164.0924 251 -31164.0924 -12268.5111 252 18088.9535 -31164.0924 253 36297.4716 18088.9535 254 1066.2229 36297.4716 255 23944.4897 1066.2229 256 24008.5743 23944.4897 257 -3218.5157 24008.5743 258 -31164.0924 -3218.5157 259 23944.4897 -31164.0924 260 -31164.0924 23944.4897 261 -31164.0924 -31164.0924 262 24904.4705 -31164.0924 263 21890.4878 24904.4705 264 -31164.0924 21890.4878 265 23944.4897 -31164.0924 266 -31164.0924 23944.4897 267 80047.1915 -31164.0924 268 -31164.0924 80047.1915 269 -55236.0503 -31164.0924 270 -31164.0924 -55236.0503 271 -31164.0924 -31164.0924 272 -31164.0924 -31164.0924 273 -28665.2822 -31164.0924 274 23944.4897 -28665.2822 275 57477.8478 23944.4897 276 22102.3950 57477.8478 277 23944.4897 22102.3950 278 57933.9111 23944.4897 279 56208.6106 57933.9111 280 -31164.0924 56208.6106 281 23944.4897 -31164.0924 282 102833.7758 23944.4897 283 23944.4897 102833.7758 284 21869.4117 23944.4897 285 23944.4897 21869.4117 286 -31164.0924 23944.4897 287 24714.2721 -31164.0924 288 -31164.0924 24714.2721 289 121394.6428 -31164.0924 290 -31164.0924 121394.6428 291 44696.2917 -31164.0924 292 41880.9308 44696.2917 293 23944.4897 41880.9308 294 11990.8023 23944.4897 295 13311.6374 11990.8023 296 23944.4897 13311.6374 297 -31164.0924 23944.4897 298 23944.4897 -31164.0924 299 -13583.2499 23944.4897 300 -39224.8440 -13583.2499 301 8133.8845 -39224.8440 302 23944.4897 8133.8845 303 -31164.0924 23944.4897 304 126341.8884 -31164.0924 305 -2952.0544 126341.8884 306 -31164.0924 -2952.0544 307 -74616.6713 -31164.0924 308 8317.1496 -74616.6713 309 -31164.0924 8317.1496 310 -31164.0924 -31164.0924 311 96925.5331 -31164.0924 312 -31164.0924 96925.5331 313 -61722.3414 -31164.0924 314 1172.4179 -61722.3414 315 59855.7592 1172.4179 316 25344.0301 59855.7592 317 -28662.1095 25344.0301 318 18498.3189 -28662.1095 319 23944.4897 18498.3189 320 23944.4897 23944.4897 321 23944.4897 23944.4897 322 23944.4897 23944.4897 323 -31164.0924 23944.4897 324 23944.4897 -31164.0924 325 21871.6492 23944.4897 326 -31164.0924 21871.6492 327 -31164.0924 -31164.0924 328 23944.4897 -31164.0924 329 23944.4897 23944.4897 330 23944.4897 23944.4897 331 23944.4897 23944.4897 332 23944.4897 23944.4897 333 18436.9685 23944.4897 334 24916.6033 18436.9685 335 23944.4897 24916.6033 336 23283.8362 23944.4897 337 23944.4897 23283.8362 338 22470.2156 23944.4897 339 12543.2403 22470.2156 340 766.5812 12543.2403 341 23944.4897 766.5812 342 23944.4897 23944.4897 343 23944.4897 23944.4897 344 23944.4897 23944.4897 345 23944.4897 23944.4897 346 12241.4031 23944.4897 347 23944.4897 12241.4031 348 23944.4897 23944.4897 349 23944.4897 23944.4897 350 40236.0807 23944.4897 351 19881.8725 40236.0807 352 16754.7417 19881.8725 353 22233.3625 16754.7417 354 23944.4897 22233.3625 355 23944.4897 23944.4897 356 94666.5242 23944.4897 357 23944.4897 94666.5242 358 23944.4897 23944.4897 359 23944.4897 23944.4897 360 17534.0406 23944.4897 361 23944.4897 17534.0406 362 23944.4897 23944.4897 363 -37653.4017 23944.4897 364 7735.2309 -37653.4017 365 23944.4897 7735.2309 366 23944.4897 23944.4897 367 23944.4897 23944.4897 368 23944.4897 23944.4897 369 23944.4897 23944.4897 370 23944.4897 23944.4897 371 88363.2415 23944.4897 372 38934.9047 88363.2415 373 6811.2441 38934.9047 374 23944.4897 6811.2441 375 23944.4897 23944.4897 376 28078.4082 23944.4897 377 14034.6561 28078.4082 378 23944.4897 14034.6561 379 23944.4897 23944.4897 380 23944.4897 23944.4897 381 108817.1919 23944.4897 382 164863.0184 108817.1919 383 23944.4897 164863.0184 384 18419.9700 23944.4897 385 -97878.2616 18419.9700 386 27285.8956 -97878.2616 387 10634.9196 27285.8956 388 16119.9896 10634.9196 389 8879.0068 16119.9896 390 15693.2771 8879.0068 391 -936.5007 15693.2771 392 -18866.5525 -936.5007 393 74608.9461 -18866.5525 394 -99649.4056 74608.9461 395 -6594.9425 -99649.4056 396 -82488.7376 -6594.9425 397 -697508.8098 -82488.7376 398 9888.3789 -697508.8098 399 121170.3199 9888.3789 400 43695.7911 121170.3199 401 28536.0182 43695.7911 402 1668.2186 28536.0182 403 -88304.8133 1668.2186 404 1638.0231 -88304.8133 405 14191.8476 1638.0231 406 34921.3543 14191.8476 407 -10054.1141 34921.3543 408 -12874.9790 -10054.1141 409 -5984.2293 -12874.9790 410 -82507.1098 -5984.2293 411 -113565.0487 -82507.1098 412 -75888.7341 -113565.0487 413 -122992.5639 -75888.7341 414 6515.1212 -122992.5639 415 -2245.6856 6515.1212 416 160267.5784 -2245.6856 417 -103626.4598 160267.5784 418 -18465.1137 -103626.4598 419 -8329.4936 -18465.1137 420 -40960.8889 -8329.4936 421 -216291.7420 -40960.8889 422 -497573.8832 -216291.7420 423 -414112.1634 -497573.8832 424 -109008.0932 -414112.1634 425 -249580.5640 -109008.0932 426 -311658.1463 -249580.5640 427 -1284038.3256 -311658.1463 428 -184370.2558 -1284038.3256 429 -134950.8637 -184370.2558 430 -377355.7461 -134950.8637 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7ni0r1291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/8ni0r1291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/9ni0r1291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/freestat/rcomp/tmp/10gric1291358215.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/111ayh1291358215.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/125sf51291358215.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/13jkue1291358215.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/14m2t21291358215.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/1583981291358215.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/16blqw1291358215.tab") + } > > try(system("convert tmp/198301291358215.ps tmp/198301291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/2kh231291358215.ps tmp/2kh231291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/3kh231291358215.ps tmp/3kh231291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/4kh231291358215.ps tmp/4kh231291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/5kh231291358215.ps tmp/5kh231291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/6icqu1291358215.ps tmp/6icqu1291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/7ni0r1291358215.ps tmp/7ni0r1291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/8ni0r1291358215.ps tmp/8ni0r1291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/9ni0r1291358215.ps tmp/9ni0r1291358215.png",intern=TRUE)) character(0) > try(system("convert tmp/10gric1291358215.ps tmp/10gric1291358215.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 15.740 3.059 16.373