R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(6282929 + ,1 + ,162556 + ,1081 + ,213118 + ,4324047 + ,1 + ,29790 + ,309 + ,81767 + ,4108272 + ,1 + ,87550 + ,458 + ,153198 + ,-1212617 + ,0 + ,84738 + ,588 + ,-26007 + ,1485329 + ,1 + ,54660 + ,299 + ,126942 + ,1779876 + ,1 + ,42634 + ,156 + ,157214 + ,1367203 + ,0 + ,40949 + ,481 + ,129352 + ,2519076 + ,1 + ,42312 + ,323 + ,234817 + ,912684 + ,1 + ,37704 + ,452 + ,60448 + ,1443586 + ,1 + ,16275 + ,109 + ,47818 + ,1220017 + ,0 + ,25830 + ,115 + ,245546 + ,984885 + ,0 + ,12679 + ,110 + ,48020 + ,1457425 + ,1 + ,18014 + ,239 + ,-1710 + ,-572920 + ,0 + ,43556 + ,247 + ,32648 + ,929144 + ,1 + ,24524 + ,497 + ,95350 + ,1151176 + ,0 + ,6532 + ,103 + ,151352 + ,790090 + ,0 + ,7123 + ,109 + ,288170 + ,774497 + ,1 + ,20813 + ,502 + ,114337 + ,990576 + ,1 + ,37597 + ,248 + ,37884 + ,454195 + ,0 + ,17821 + ,373 + ,122844 + ,876607 + ,1 + ,12988 + ,119 + ,82340 + ,711969 + ,1 + ,22330 + ,84 + ,79801 + ,702380 + ,0 + ,13326 + ,102 + ,165548 + ,264449 + ,0 + ,16189 + ,295 + ,116384 + ,450033 + ,0 + ,7146 + ,105 + ,134028 + ,541063 + ,0 + ,15824 + ,64 + ,63838 + ,588864 + ,1 + ,26088 + ,267 + ,74996 + ,-37216 + ,0 + ,11326 + ,129 + ,31080 + ,783310 + ,0 + ,8568 + ,37 + ,32168 + ,467359 + ,0 + ,14416 + ,361 + ,49857 + ,688779 + ,1 + ,3369 + ,28 + ,87161 + ,608419 + ,1 + ,11819 + ,85 + ,106113 + ,696348 + ,1 + ,6620 + ,44 + ,80570 + ,597793 + ,1 + ,4519 + ,49 + ,102129 + ,821730 + ,0 + ,2220 + ,22 + ,301670 + ,377934 + ,0 + ,18562 + ,155 + ,102313 + ,651939 + ,0 + ,10327 + ,91 + ,88577 + ,697458 + ,1 + ,5336 + ,81 + ,112477 + ,700368 + ,1 + ,2365 + ,79 + ,191778 + ,225986 + ,0 + ,4069 + ,145 + ,79804 + ,348695 + ,0 + ,7710 + ,816 + ,128294 + ,373683 + ,0 + ,13718 + ,61 + ,96448 + ,501709 + ,0 + ,4525 + ,226 + ,93811 + ,413743 + ,0 + ,6869 + ,105 + ,117520 + ,379825 + ,0 + ,4628 + ,62 + ,69159 + ,336260 + ,1 + ,3653 + ,24 + ,101792 + ,636765 + ,1 + ,1265 + ,26 + ,210568 + ,481231 + ,1 + ,7489 + ,322 + ,136996 + ,469107 + ,0 + ,4901 + ,84 + ,121920 + ,211928 + ,0 + ,2284 + ,33 + ,76403 + ,563925 + ,1 + ,3160 + ,108 + ,108094 + ,511939 + ,1 + ,4150 + ,150 + ,134759 + ,521016 + ,1 + ,7285 + ,115 + ,188873 + ,543856 + ,1 + ,1134 + ,162 + ,146216 + ,329304 + ,1 + ,4658 + ,158 + ,156608 + ,423262 + ,0 + ,2384 + ,97 + ,61348 + ,509665 + ,0 + ,3748 + ,9 + ,50350 + ,455881 + ,0 + ,5371 + ,66 + ,87720 + ,367772 + ,0 + ,1285 + ,107 + ,99489 + ,406339 + ,1 + ,9327 + ,101 + ,87419 + ,493408 + ,1 + ,5565 + ,47 + ,94355 + ,232942 + ,0 + ,1528 + ,38 + ,60326 + ,416002 + ,1 + ,3122 + ,34 + ,94670 + ,337430 + ,1 + ,7317 + ,84 + ,82425 + ,361517 + ,0 + ,2675 + ,79 + ,59017 + ,360962 + ,0 + ,13253 + ,947 + ,90829 + ,235561 + ,0 + ,880 + ,74 + ,80791 + ,408247 + ,1 + ,2053 + ,53 + ,100423 + ,450296 + ,0 + ,1424 + ,94 + ,131116 + ,418799 + ,1 + ,4036 + ,63 + ,100269 + ,247405 + ,1 + ,3045 + ,58 + ,27330 + ,378519 + ,0 + ,5119 + ,49 + ,39039 + ,326638 + ,0 + ,1431 + ,34 + ,106885 + ,328233 + ,0 + ,554 + ,11 + ,79285 + ,386225 + ,0 + ,1975 + ,35 + ,118881 + ,283662 + ,1 + ,1286 + ,17 + ,77623 + ,370225 + ,0 + ,1012 + ,47 + ,114768 + ,269236 + ,0 + ,810 + ,43 + ,74015 + ,365732 + ,0 + ,1280 + ,117 + ,69465 + ,420383 + ,1 + ,666 + ,171 + ,117869 + ,345811 + ,0 + ,1380 + ,26 + ,60982 + ,431809 + ,1 + ,4608 + ,73 + ,90131 + ,418876 + ,0 + ,876 + ,59 + ,138971 + ,297476 + ,0 + ,814 + ,18 + ,39625 + ,416776 + ,0 + ,514 + ,15 + ,102725 + ,357257 + ,1 + ,5692 + ,72 + ,64239 + ,458343 + ,0 + ,3642 + ,86 + ,90262 + ,388386 + ,0 + ,540 + ,14 + ,103960 + ,358934 + ,0 + ,2099 + ,64 + ,106611 + ,407560 + ,0 + ,567 + ,11 + ,103345 + ,392558 + ,0 + ,2001 + ,52 + ,95551 + ,373177 + ,1 + ,2949 + ,41 + ,82903 + ,428370 + ,0 + ,2253 + ,99 + ,63593 + ,369419 + ,1 + ,6533 + ,75 + ,126910 + ,358649 + ,0 + ,1889 + ,45 + ,37527 + ,376641 + ,1 + ,3055 + ,43 + ,60247 + ,467427 + ,0 + ,272 + ,8 + ,112995 + ,364885 + ,1 + ,1414 + ,198 + ,70184 + ,436230 + ,0 + ,2564 + ,22 + ,130140 + ,329118 + ,1 + ,1383 + ,11 + ,73221 + ,317365 + ,1 + ,1261 + ,33 + ,76114 + ,286849 + ,0 + ,975 + ,23 + ,90534 + ,376685 + ,0 + ,3366 + ,80 + ,108479 + ,407198 + ,0 + ,576 + ,18 + ,113761 + ,377772 + ,0 + ,1306 + ,28 + ,68696 + ,271483 + ,0 + ,746 + ,23 + ,71561 + ,153661 + ,1 + ,3192 + ,60 + ,59831 + ,513294 + ,1 + ,2045 + ,20 + ,97890 + ,324881 + ,0 + ,5477 + ,59 + ,101481 + ,264512 + ,1 + ,1932 + ,36 + ,72954 + ,420968 + ,0 + ,936 + ,30 + ,67939 + ,129302 + ,1 + ,3437 + ,47 + ,48022 + ,191521 + ,0 + ,5131 + ,71 + ,86111 + ,268673 + ,1 + ,2397 + ,14 + ,74020 + ,353179 + ,1 + ,1389 + ,9 + ,57530 + ,354624 + ,0 + ,1503 + ,39 + ,56364 + ,363713 + ,0 + ,402 + ,26 + ,84990 + ,456657 + ,0 + ,2239 + ,21 + ,88590 + ,211742 + ,1 + ,2234 + ,16 + ,77200 + ,338381 + ,0 + ,837 + ,69 + ,61262 + ,418530 + ,0 + ,10579 + ,92 + ,110309 + ,351483 + ,0 + ,875 + ,14 + ,67000 + ,372928 + ,0 + ,1395 + ,103 + ,93099 + ,485538 + ,1 + ,1659 + ,29 + ,107577 + ,279268 + ,1 + ,2647 + ,37 + ,62920 + ,219060 + ,1 + ,3294 + ,23 + ,75832 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325314 + ,0 + ,94 + ,7 + ,60793 + ,322046 + ,0 + ,422 + ,28 + ,57935 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325599 + ,0 + ,34 + ,8 + ,60630 + ,377028 + ,0 + ,1558 + ,63 + ,55637 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,323850 + ,0 + ,43 + ,3 + ,60887 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,331514 + ,0 + ,316 + ,9 + ,60505 + ,325632 + ,0 + ,115 + ,13 + ,60945 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,322265 + ,0 + ,389 + ,14 + ,58990 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325906 + ,0 + ,1002 + ,15 + ,56750 + ,325985 + ,0 + ,36 + ,3 + ,60894 + ,346145 + ,0 + ,460 + ,15 + ,63346 + ,325898 + ,0 + ,309 + ,11 + ,56535 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325356 + ,0 + ,9 + ,6 + ,60835 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325930 + ,0 + ,14 + ,1 + ,61016 + ,318020 + ,0 + ,520 + ,10 + ,58650 + ,326389 + ,0 + ,1766 + ,73 + ,60438 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,302925 + ,0 + ,458 + ,11 + ,58625 + ,325540 + ,0 + ,20 + ,3 + ,60938 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,326736 + ,0 + ,98 + ,2 + ,61490 + ,340580 + ,0 + ,405 + ,7 + ,60845 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,331828 + ,0 + ,483 + ,27 + ,60830 + ,323299 + ,0 + ,454 + ,51 + ,63261 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,387722 + ,0 + ,757 + ,19 + ,45689 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,324598 + ,0 + ,36 + ,4 + ,61564 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,328726 + ,0 + ,203 + ,9 + ,61938 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325043 + ,0 + ,90 + ,8 + ,60951 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325806 + ,1 + ,71 + ,1 + ,60745 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,387732 + ,0 + ,972 + ,34 + ,71642 + ,349729 + ,1 + ,531 + ,10 + ,71641 + ,332202 + ,0 + ,604 + ,38 + ,55792 + ,305442 + ,1 + ,283 + ,10 + ,71873 + ,329537 + ,1 + ,23 + ,5 + ,62555 + ,327055 + ,1 + ,638 + ,14 + ,60370 + ,356245 + ,1 + ,699 + ,16 + ,64873 + ,328451 + ,0 + ,149 + ,5 + ,62041 + ,307062 + ,0 + ,226 + ,5 + ,65745 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,331345 + ,0 + ,275 + ,4 + ,59500 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,331824 + ,0 + ,141 + ,6 + ,61630 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325685 + ,0 + ,28 + ,2 + ,60890 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,404480 + ,1 + ,2566 + ,80 + ,113521 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,318314 + ,1 + ,472 + ,20 + ,80045 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,311807 + ,1 + ,203 + ,27 + ,50804 + ,337724 + ,1 + ,496 + ,17 + ,87390 + ,326431 + ,1 + ,10 + ,2 + ,61656 + ,327556 + ,1 + ,63 + ,4 + ,65688 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,356850 + ,1 + ,1136 + ,32 + ,48522 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,322741 + ,0 + ,267 + ,32 + ,57640 + ,310902 + ,0 + ,474 + ,20 + ,61977 + ,324295 + ,0 + ,534 + ,7 + ,62620 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,326156 + ,0 + ,15 + ,8 + ,60831 + ,326960 + ,0 + ,397 + ,28 + ,60646 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,333411 + ,0 + ,1061 + ,20 + ,56225 + ,297761 + ,0 + ,288 + ,4 + ,60510 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325536 + ,0 + ,3 + ,2 + ,60698 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325762 + ,0 + ,20 + ,2 + ,60805 + ,327957 + ,0 + ,278 + ,26 + ,61404 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,318521 + ,0 + ,192 + ,4 + ,65276 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,319775 + ,0 + ,317 + ,9 + ,63915 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,332128 + ,0 + ,368 + ,17 + ,61686 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325486 + ,0 + ,2 + ,1 + ,60743 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325838 + ,0 + ,53 + ,6 + ,60349 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,331767 + ,0 + ,94 + ,3 + ,61360 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,324523 + ,0 + ,24 + ,8 + ,59818 + ,339995 + ,0 + ,2332 + ,4 + ,72680 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,319582 + ,0 + ,131 + ,11 + ,61808 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,307245 + ,0 + ,206 + ,9 + ,53110 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,317967 + ,0 + ,167 + ,2 + ,64245 + ,331488 + ,0 + ,622 + ,73 + ,73007 + ,335452 + ,0 + ,885 + ,85 + ,82732 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,334184 + ,0 + ,365 + ,8 + ,54820 + ,313213 + ,0 + ,364 + ,35 + ,47705 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,348678 + ,0 + ,226 + ,12 + ,72835 + ,328727 + ,0 + ,307 + ,15 + ,58856 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,387978 + ,0 + ,188 + ,11 + ,77655 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,336704 + ,0 + ,138 + ,6 + ,69817 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,322076 + ,0 + ,125 + ,12 + ,60798 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,334272 + ,0 + ,282 + ,30 + ,62452 + ,338197 + ,0 + ,335 + ,33 + ,64175 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,321024 + ,1 + ,813 + ,82 + ,67440 + ,322145 + ,0 + ,176 + ,28 + ,68136 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,323351 + ,0 + ,249 + ,72 + ,56726 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,327748 + ,0 + ,333 + ,13 + ,70811 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,328157 + ,0 + ,30 + ,4 + ,62045 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,311594 + ,0 + ,249 + ,62 + ,54323 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,335962 + ,0 + ,165 + ,24 + ,62841 + ,372426 + ,0 + ,453 + ,21 + ,81125 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,319844 + ,0 + ,53 + ,14 + ,59506 + ,355822 + ,1 + ,382 + ,21 + ,59365 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,324047 + ,1 + ,30 + ,4 + ,60798 + ,311464 + ,0 + ,290 + ,2 + ,58790 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,353417 + ,0 + ,366 + ,53 + ,61808 + ,325590 + ,0 + ,2 + ,9 + ,60735 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,328576 + ,1 + ,209 + ,13 + ,64016 + ,326126 + ,0 + ,384 + ,22 + ,54683 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,369376 + ,0 + ,365 + ,83 + ,87192 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,332013 + ,1 + ,49 + ,8 + ,64107 + ,325871 + ,0 + ,3 + ,4 + ,60761 + ,342165 + ,0 + ,133 + ,14 + ,65990 + ,324967 + ,0 + ,32 + ,1 + ,59988 + ,314832 + ,0 + ,368 + ,17 + ,61167 + ,325557 + ,0 + ,1 + ,6 + ,60719 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,322649 + ,0 + ,22 + ,2 + ,60722 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,324598 + ,0 + ,96 + ,5 + ,60379 + ,325567 + ,0 + ,1 + ,2 + ,60727 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,324005 + ,0 + ,81 + ,7 + ,60925 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325748 + ,0 + ,26 + ,1 + ,60896 + ,323385 + ,0 + ,125 + ,13 + ,59734 + ,315409 + ,0 + ,304 + ,15 + ,62969 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,312275 + ,0 + ,119 + ,6 + ,59118 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,320576 + ,0 + ,312 + ,14 + ,58598 + ,325246 + ,0 + ,60 + ,10 + ,61124 + ,332961 + ,0 + ,587 + ,12 + ,59595 + ,323010 + ,0 + ,135 + ,2 + ,62065 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,345253 + ,0 + ,514 + ,52 + ,78780 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325559 + ,0 + ,1 + ,4 + ,60722 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,319634 + ,1 + ,58 + ,3 + ,61600 + ,319951 + ,0 + ,180 + ,11 + ,59635 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,318519 + ,0 + ,448 + ,40 + ,59781 + ,343222 + ,0 + ,227 + ,9 + ,76644 + ,317234 + ,0 + ,174 + ,1 + ,64820 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,314025 + ,0 + ,121 + ,24 + ,56178 + ,320249 + ,0 + ,607 + ,11 + ,60436 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,349365 + ,0 + ,530 + ,60 + ,73433 + ,289197 + ,0 + ,571 + ,80 + ,41477 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,329245 + ,0 + ,78 + ,16 + ,62700 + ,240869 + ,0 + ,2489 + ,40 + ,67804 + ,327182 + ,0 + ,131 + ,6 + ,59661 + ,322876 + ,0 + ,923 + ,8 + ,58620 + ,323117 + ,0 + ,72 + ,3 + ,60398 + ,306351 + ,0 + ,572 + ,16 + ,58580 + ,335137 + ,0 + ,397 + ,10 + ,62710 + ,308271 + ,0 + ,450 + ,8 + ,59325 + ,301731 + ,0 + ,622 + ,7 + ,60950 + ,382409 + ,0 + ,694 + ,8 + ,68060 + ,279230 + ,1 + ,3425 + ,12 + ,83620 + ,298731 + ,0 + ,562 + ,13 + ,58456 + ,243650 + ,0 + ,4917 + ,42 + ,52811 + ,532682 + ,1 + ,1442 + ,118 + ,121173 + ,319771 + ,0 + ,529 + ,9 + ,63870 + ,171493 + ,1 + ,2126 + ,138 + ,21001 + ,347262 + ,0 + ,1061 + ,5 + ,70415 + ,343945 + ,0 + ,776 + ,9 + ,64230 + ,311874 + ,0 + ,611 + ,8 + ,59190 + ,302211 + ,1 + ,1526 + ,25 + ,69351 + ,316708 + ,0 + ,592 + ,7 + ,64270 + ,333463 + ,0 + ,1182 + ,13 + ,70694 + ,344282 + ,0 + ,621 + ,16 + ,68005 + ,319635 + ,0 + ,989 + ,11 + ,58930 + ,301186 + ,0 + ,438 + ,11 + ,58320 + ,300381 + ,0 + ,726 + ,3 + ,69980 + ,318765 + ,0 + ,1303 + ,61 + ,69863 + ,286146 + ,1 + ,6341 + ,24 + ,63255 + ,306844 + ,1 + ,1164 + ,17 + ,57320 + ,307705 + ,1 + ,3310 + ,33 + ,75230 + ,312448 + ,0 + ,1366 + ,7 + ,79420 + ,299715 + ,0 + ,965 + ,3 + ,73490 + ,373399 + ,0 + ,3256 + ,66 + ,35250 + ,299446 + ,1 + ,1135 + ,17 + ,62285 + ,325586 + ,0 + ,1270 + ,26 + ,69206 + ,291221 + ,0 + ,661 + ,3 + ,65920 + ,261173 + ,0 + ,1013 + ,2 + ,69770 + ,255027 + ,0 + ,2844 + ,67 + ,72683 + ,-78375 + ,1 + ,11528 + ,70 + ,-14545 + ,-58143 + ,0 + ,6526 + ,26 + ,55830 + ,227033 + ,0 + ,2264 + ,24 + ,55174 + ,235098 + ,1 + ,4461 + ,94 + ,67038 + ,21267 + ,0 + ,3999 + ,30 + ,51252 + ,238675 + ,0 + ,35624 + ,223 + ,157278 + ,197687 + ,0 + ,9252 + ,48 + ,79510 + ,418341 + ,0 + ,15236 + ,90 + ,77440 + ,-297706 + ,0 + ,18073 + ,180 + ,27284) + ,dim=c(5 + ,431) + ,dimnames=list(c('Wealth' + ,'Group' + ,'Costs' + ,'Trades' + ,'Dividends ') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('Wealth','Group','Costs','Trades','Dividends '),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 = 'Include Monthly Dummies' > par1 = '1' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Wealth Group Costs Trades Dividends\r M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 M11 1 6282929 1 162556 1081 213118 1 0 0 0 0 0 0 0 0 0 0 2 4324047 1 29790 309 81767 0 1 0 0 0 0 0 0 0 0 0 3 4108272 1 87550 458 153198 0 0 1 0 0 0 0 0 0 0 0 4 -1212617 0 84738 588 -26007 0 0 0 1 0 0 0 0 0 0 0 5 1485329 1 54660 299 126942 0 0 0 0 1 0 0 0 0 0 0 6 1779876 1 42634 156 157214 0 0 0 0 0 1 0 0 0 0 0 7 1367203 0 40949 481 129352 0 0 0 0 0 0 1 0 0 0 0 8 2519076 1 42312 323 234817 0 0 0 0 0 0 0 1 0 0 0 9 912684 1 37704 452 60448 0 0 0 0 0 0 0 0 1 0 0 10 1443586 1 16275 109 47818 0 0 0 0 0 0 0 0 0 1 0 11 1220017 0 25830 115 245546 0 0 0 0 0 0 0 0 0 0 1 12 984885 0 12679 110 48020 0 0 0 0 0 0 0 0 0 0 0 13 1457425 1 18014 239 -1710 1 0 0 0 0 0 0 0 0 0 0 14 -572920 0 43556 247 32648 0 1 0 0 0 0 0 0 0 0 0 15 929144 1 24524 497 95350 0 0 1 0 0 0 0 0 0 0 0 16 1151176 0 6532 103 151352 0 0 0 1 0 0 0 0 0 0 0 17 790090 0 7123 109 288170 0 0 0 0 1 0 0 0 0 0 0 18 774497 1 20813 502 114337 0 0 0 0 0 1 0 0 0 0 0 19 990576 1 37597 248 37884 0 0 0 0 0 0 1 0 0 0 0 20 454195 0 17821 373 122844 0 0 0 0 0 0 0 1 0 0 0 21 876607 1 12988 119 82340 0 0 0 0 0 0 0 0 1 0 0 22 711969 1 22330 84 79801 0 0 0 0 0 0 0 0 0 1 0 23 702380 0 13326 102 165548 0 0 0 0 0 0 0 0 0 0 1 24 264449 0 16189 295 116384 0 0 0 0 0 0 0 0 0 0 0 25 450033 0 7146 105 134028 1 0 0 0 0 0 0 0 0 0 0 26 541063 0 15824 64 63838 0 1 0 0 0 0 0 0 0 0 0 27 588864 1 26088 267 74996 0 0 1 0 0 0 0 0 0 0 0 28 -37216 0 11326 129 31080 0 0 0 1 0 0 0 0 0 0 0 29 783310 0 8568 37 32168 0 0 0 0 1 0 0 0 0 0 0 30 467359 0 14416 361 49857 0 0 0 0 0 1 0 0 0 0 0 31 688779 1 3369 28 87161 0 0 0 0 0 0 1 0 0 0 0 32 608419 1 11819 85 106113 0 0 0 0 0 0 0 1 0 0 0 33 696348 1 6620 44 80570 0 0 0 0 0 0 0 0 1 0 0 34 597793 1 4519 49 102129 0 0 0 0 0 0 0 0 0 1 0 35 821730 0 2220 22 301670 0 0 0 0 0 0 0 0 0 0 1 36 377934 0 18562 155 102313 0 0 0 0 0 0 0 0 0 0 0 37 651939 0 10327 91 88577 1 0 0 0 0 0 0 0 0 0 0 38 697458 1 5336 81 112477 0 1 0 0 0 0 0 0 0 0 0 39 700368 1 2365 79 191778 0 0 1 0 0 0 0 0 0 0 0 40 225986 0 4069 145 79804 0 0 0 1 0 0 0 0 0 0 0 41 348695 0 7710 816 128294 0 0 0 0 1 0 0 0 0 0 0 42 373683 0 13718 61 96448 0 0 0 0 0 1 0 0 0 0 0 43 501709 0 4525 226 93811 0 0 0 0 0 0 1 0 0 0 0 44 413743 0 6869 105 117520 0 0 0 0 0 0 0 1 0 0 0 45 379825 0 4628 62 69159 0 0 0 0 0 0 0 0 1 0 0 46 336260 1 3653 24 101792 0 0 0 0 0 0 0 0 0 1 0 47 636765 1 1265 26 210568 0 0 0 0 0 0 0 0 0 0 1 48 481231 1 7489 322 136996 0 0 0 0 0 0 0 0 0 0 0 49 469107 0 4901 84 121920 1 0 0 0 0 0 0 0 0 0 0 50 211928 0 2284 33 76403 0 1 0 0 0 0 0 0 0 0 0 51 563925 1 3160 108 108094 0 0 1 0 0 0 0 0 0 0 0 52 511939 1 4150 150 134759 0 0 0 1 0 0 0 0 0 0 0 53 521016 1 7285 115 188873 0 0 0 0 1 0 0 0 0 0 0 54 543856 1 1134 162 146216 0 0 0 0 0 1 0 0 0 0 0 55 329304 1 4658 158 156608 0 0 0 0 0 0 1 0 0 0 0 56 423262 0 2384 97 61348 0 0 0 0 0 0 0 1 0 0 0 57 509665 0 3748 9 50350 0 0 0 0 0 0 0 0 1 0 0 58 455881 0 5371 66 87720 0 0 0 0 0 0 0 0 0 1 0 59 367772 0 1285 107 99489 0 0 0 0 0 0 0 0 0 0 1 60 406339 1 9327 101 87419 0 0 0 0 0 0 0 0 0 0 0 61 493408 1 5565 47 94355 1 0 0 0 0 0 0 0 0 0 0 62 232942 0 1528 38 60326 0 1 0 0 0 0 0 0 0 0 0 63 416002 1 3122 34 94670 0 0 1 0 0 0 0 0 0 0 0 64 337430 1 7317 84 82425 0 0 0 1 0 0 0 0 0 0 0 65 361517 0 2675 79 59017 0 0 0 0 1 0 0 0 0 0 0 66 360962 0 13253 947 90829 0 0 0 0 0 1 0 0 0 0 0 67 235561 0 880 74 80791 0 0 0 0 0 0 1 0 0 0 0 68 408247 1 2053 53 100423 0 0 0 0 0 0 0 1 0 0 0 69 450296 0 1424 94 131116 0 0 0 0 0 0 0 0 1 0 0 70 418799 1 4036 63 100269 0 0 0 0 0 0 0 0 0 1 0 71 247405 1 3045 58 27330 0 0 0 0 0 0 0 0 0 0 1 72 378519 0 5119 49 39039 0 0 0 0 0 0 0 0 0 0 0 73 326638 0 1431 34 106885 1 0 0 0 0 0 0 0 0 0 0 74 328233 0 554 11 79285 0 1 0 0 0 0 0 0 0 0 0 75 386225 0 1975 35 118881 0 0 1 0 0 0 0 0 0 0 0 76 283662 1 1286 17 77623 0 0 0 1 0 0 0 0 0 0 0 77 370225 0 1012 47 114768 0 0 0 0 1 0 0 0 0 0 0 78 269236 0 810 43 74015 0 0 0 0 0 1 0 0 0 0 0 79 365732 0 1280 117 69465 0 0 0 0 0 0 1 0 0 0 0 80 420383 1 666 171 117869 0 0 0 0 0 0 0 1 0 0 0 81 345811 0 1380 26 60982 0 0 0 0 0 0 0 0 1 0 0 82 431809 1 4608 73 90131 0 0 0 0 0 0 0 0 0 1 0 83 418876 0 876 59 138971 0 0 0 0 0 0 0 0 0 0 1 84 297476 0 814 18 39625 0 0 0 0 0 0 0 0 0 0 0 85 416776 0 514 15 102725 1 0 0 0 0 0 0 0 0 0 0 86 357257 1 5692 72 64239 0 1 0 0 0 0 0 0 0 0 0 87 458343 0 3642 86 90262 0 0 1 0 0 0 0 0 0 0 0 88 388386 0 540 14 103960 0 0 0 1 0 0 0 0 0 0 0 89 358934 0 2099 64 106611 0 0 0 0 1 0 0 0 0 0 0 90 407560 0 567 11 103345 0 0 0 0 0 1 0 0 0 0 0 91 392558 0 2001 52 95551 0 0 0 0 0 0 1 0 0 0 0 92 373177 1 2949 41 82903 0 0 0 0 0 0 0 1 0 0 0 93 428370 0 2253 99 63593 0 0 0 0 0 0 0 0 1 0 0 94 369419 1 6533 75 126910 0 0 0 0 0 0 0 0 0 1 0 95 358649 0 1889 45 37527 0 0 0 0 0 0 0 0 0 0 1 96 376641 1 3055 43 60247 0 0 0 0 0 0 0 0 0 0 0 97 467427 0 272 8 112995 1 0 0 0 0 0 0 0 0 0 0 98 364885 1 1414 198 70184 0 1 0 0 0 0 0 0 0 0 0 99 436230 0 2564 22 130140 0 0 1 0 0 0 0 0 0 0 0 100 329118 1 1383 11 73221 0 0 0 1 0 0 0 0 0 0 0 101 317365 1 1261 33 76114 0 0 0 0 1 0 0 0 0 0 0 102 286849 0 975 23 90534 0 0 0 0 0 1 0 0 0 0 0 103 376685 0 3366 80 108479 0 0 0 0 0 0 1 0 0 0 0 104 407198 0 576 18 113761 0 0 0 0 0 0 0 1 0 0 0 105 377772 0 1306 28 68696 0 0 0 0 0 0 0 0 1 0 0 106 271483 0 746 23 71561 0 0 0 0 0 0 0 0 0 1 0 107 153661 1 3192 60 59831 0 0 0 0 0 0 0 0 0 0 1 108 513294 1 2045 20 97890 0 0 0 0 0 0 0 0 0 0 0 109 324881 0 5477 59 101481 1 0 0 0 0 0 0 0 0 0 0 110 264512 1 1932 36 72954 0 1 0 0 0 0 0 0 0 0 0 111 420968 0 936 30 67939 0 0 1 0 0 0 0 0 0 0 0 112 129302 1 3437 47 48022 0 0 0 1 0 0 0 0 0 0 0 113 191521 0 5131 71 86111 0 0 0 0 1 0 0 0 0 0 0 114 268673 1 2397 14 74020 0 0 0 0 0 1 0 0 0 0 0 115 353179 1 1389 9 57530 0 0 0 0 0 0 1 0 0 0 0 116 354624 0 1503 39 56364 0 0 0 0 0 0 0 1 0 0 0 117 363713 0 402 26 84990 0 0 0 0 0 0 0 0 1 0 0 118 456657 0 2239 21 88590 0 0 0 0 0 0 0 0 0 1 0 119 211742 1 2234 16 77200 0 0 0 0 0 0 0 0 0 0 1 120 338381 0 837 69 61262 0 0 0 0 0 0 0 0 0 0 0 121 418530 0 10579 92 110309 1 0 0 0 0 0 0 0 0 0 0 122 351483 0 875 14 67000 0 1 0 0 0 0 0 0 0 0 0 123 372928 0 1395 103 93099 0 0 1 0 0 0 0 0 0 0 0 124 485538 1 1659 29 107577 0 0 0 1 0 0 0 0 0 0 0 125 279268 1 2647 37 62920 0 0 0 0 1 0 0 0 0 0 0 126 219060 1 3294 23 75832 0 0 0 0 0 1 0 0 0 0 0 127 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 128 325314 0 94 7 60793 0 0 0 0 0 0 0 1 0 0 0 129 322046 0 422 28 57935 0 0 0 0 0 0 0 0 1 0 0 130 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 131 325599 0 34 8 60630 0 0 0 0 0 0 0 0 0 0 1 132 377028 0 1558 63 55637 0 0 0 0 0 0 0 0 0 0 0 133 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 134 323850 0 43 3 60887 0 1 0 0 0 0 0 0 0 0 0 135 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 136 331514 0 316 9 60505 0 0 0 1 0 0 0 0 0 0 0 137 325632 0 115 13 60945 0 0 0 0 1 0 0 0 0 0 0 138 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 139 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 140 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 141 322265 0 389 14 58990 0 0 0 0 0 0 0 0 1 0 0 142 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 143 325906 0 1002 15 56750 0 0 0 0 0 0 0 0 0 0 1 144 325985 0 36 3 60894 0 0 0 0 0 0 0 0 0 0 0 145 346145 0 460 15 63346 1 0 0 0 0 0 0 0 0 0 0 146 325898 0 309 11 56535 0 1 0 0 0 0 0 0 0 0 0 147 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 148 325356 0 9 6 60835 0 0 0 1 0 0 0 0 0 0 0 149 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0 150 325930 0 14 1 61016 0 0 0 0 0 1 0 0 0 0 0 151 318020 0 520 10 58650 0 0 0 0 0 0 1 0 0 0 0 152 326389 0 1766 73 60438 0 0 0 0 0 0 0 1 0 0 0 153 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 154 302925 0 458 11 58625 0 0 0 0 0 0 0 0 0 1 0 155 325540 0 20 3 60938 0 0 0 0 0 0 0 0 0 0 1 156 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 157 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 158 326736 0 98 2 61490 0 1 0 0 0 0 0 0 0 0 0 159 340580 0 405 7 60845 0 0 1 0 0 0 0 0 0 0 0 160 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 161 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0 162 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 163 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 164 331828 0 483 27 60830 0 0 0 0 0 0 0 1 0 0 0 165 323299 0 454 51 63261 0 0 0 0 0 0 0 0 1 0 0 166 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 167 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 168 387722 0 757 19 45689 0 0 0 0 0 0 0 0 0 0 0 169 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 170 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 171 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 172 324598 0 36 4 61564 0 0 0 1 0 0 0 0 0 0 0 173 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0 174 328726 0 203 9 61938 0 0 0 0 0 1 0 0 0 0 0 175 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 176 325043 0 90 8 60951 0 0 0 0 0 0 0 1 0 0 0 177 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 178 325806 1 71 1 60745 0 0 0 0 0 0 0 0 0 1 0 179 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 180 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 181 387732 0 972 34 71642 1 0 0 0 0 0 0 0 0 0 0 182 349729 1 531 10 71641 0 1 0 0 0 0 0 0 0 0 0 183 332202 0 604 38 55792 0 0 1 0 0 0 0 0 0 0 0 184 305442 1 283 10 71873 0 0 0 1 0 0 0 0 0 0 0 185 329537 1 23 5 62555 0 0 0 0 1 0 0 0 0 0 0 186 327055 1 638 14 60370 0 0 0 0 0 1 0 0 0 0 0 187 356245 1 699 16 64873 0 0 0 0 0 0 1 0 0 0 0 188 328451 0 149 5 62041 0 0 0 0 0 0 0 1 0 0 0 189 307062 0 226 5 65745 0 0 0 0 0 0 0 0 1 0 0 190 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 191 331345 0 275 4 59500 0 0 0 0 0 0 0 0 0 0 1 192 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 193 331824 0 141 6 61630 1 0 0 0 0 0 0 0 0 0 0 194 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 195 325685 0 28 2 60890 0 0 1 0 0 0 0 0 0 0 0 196 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 197 404480 1 2566 80 113521 0 0 0 0 1 0 0 0 0 0 0 198 325560 1 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 199 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 200 318314 1 472 20 80045 0 0 0 0 0 0 0 1 0 0 0 201 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 202 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 203 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 204 311807 1 203 27 50804 0 0 0 0 0 0 0 0 0 0 0 205 337724 1 496 17 87390 1 0 0 0 0 0 0 0 0 0 0 206 326431 1 10 2 61656 0 1 0 0 0 0 0 0 0 0 0 207 327556 1 63 4 65688 0 0 1 0 0 0 0 0 0 0 0 208 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 209 356850 1 1136 32 48522 0 0 0 0 1 0 0 0 0 0 0 210 325560 1 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 211 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 212 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 213 322741 0 267 32 57640 0 0 0 0 0 0 0 0 1 0 0 214 310902 0 474 20 61977 0 0 0 0 0 0 0 0 0 1 0 215 324295 0 534 7 62620 0 0 0 0 0 0 0 0 0 0 1 216 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 217 326156 0 15 8 60831 1 0 0 0 0 0 0 0 0 0 0 218 326960 0 397 28 60646 0 1 0 0 0 0 0 0 0 0 0 219 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 220 333411 0 1061 20 56225 0 0 0 1 0 0 0 0 0 0 0 221 297761 0 288 4 60510 0 0 0 0 1 0 0 0 0 0 0 222 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 223 325536 0 3 2 60698 0 0 0 0 0 0 1 0 0 0 0 224 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 225 325762 0 20 2 60805 0 0 0 0 0 0 0 0 1 0 0 226 327957 0 278 26 61404 0 0 0 0 0 0 0 0 0 1 0 227 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 228 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 229 318521 0 192 4 65276 1 0 0 0 0 0 0 0 0 0 0 230 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 231 319775 0 317 9 63915 0 0 1 0 0 0 0 0 0 0 0 232 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 233 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0 234 332128 0 368 17 61686 0 0 0 0 0 1 0 0 0 0 0 235 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 236 325486 0 2 1 60743 0 0 0 0 0 0 0 1 0 0 0 237 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 238 325838 0 53 6 60349 0 0 0 0 0 0 0 0 0 1 0 239 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 240 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 241 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 242 331767 0 94 3 61360 0 1 0 0 0 0 0 0 0 0 0 243 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 244 324523 0 24 8 59818 0 0 0 1 0 0 0 0 0 0 0 245 339995 0 2332 4 72680 0 0 0 0 1 0 0 0 0 0 0 246 325560 1 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 247 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 248 319582 0 131 11 61808 0 0 0 0 0 0 0 1 0 0 0 249 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 250 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 251 307245 0 206 9 53110 0 0 0 0 0 0 0 0 0 0 1 252 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 253 317967 0 167 2 64245 1 0 0 0 0 0 0 0 0 0 0 254 331488 0 622 73 73007 0 1 0 0 0 0 0 0 0 0 0 255 335452 0 885 85 82732 0 0 1 0 0 0 0 0 0 0 0 256 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 257 334184 0 365 8 54820 0 0 0 0 1 0 0 0 0 0 0 258 313213 0 364 35 47705 0 0 0 0 0 1 0 0 0 0 0 259 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 260 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 261 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 262 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 263 348678 0 226 12 72835 0 0 0 0 0 0 0 0 0 0 1 264 328727 0 307 15 58856 0 0 0 0 0 0 0 0 0 0 0 265 325560 1 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 266 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 267 325560 1 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 268 387978 0 188 11 77655 0 0 0 1 0 0 0 0 0 0 0 269 325560 1 0 0 60720 0 0 0 0 1 0 0 0 0 0 0 270 336704 0 138 6 69817 0 0 0 0 0 1 0 0 0 0 0 271 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 272 325560 1 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 273 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 274 322076 0 125 12 60798 0 0 0 0 0 0 0 0 0 1 0 275 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 276 334272 0 282 30 62452 0 0 0 0 0 0 0 0 0 0 0 277 338197 0 335 33 64175 1 0 0 0 0 0 0 0 0 0 0 278 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 279 321024 1 813 82 67440 0 0 1 0 0 0 0 0 0 0 0 280 322145 0 176 28 68136 0 0 0 1 0 0 0 0 0 0 0 281 325560 1 0 0 60720 0 0 0 0 1 0 0 0 0 0 0 282 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 283 323351 0 249 72 56726 0 0 0 0 0 0 1 0 0 0 0 284 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 285 327748 0 333 13 70811 0 0 0 0 0 0 0 0 1 0 0 286 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 287 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 288 328157 0 30 4 62045 0 0 0 0 0 0 0 0 0 0 0 289 325560 1 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 290 311594 0 249 62 54323 0 1 0 0 0 0 0 0 0 0 0 291 325560 1 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 292 335962 0 165 24 62841 0 0 0 1 0 0 0 0 0 0 0 293 372426 0 453 21 81125 0 0 0 0 1 0 0 0 0 0 0 294 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 295 319844 0 53 14 59506 0 0 0 0 0 0 1 0 0 0 0 296 355822 1 382 21 59365 0 0 0 0 0 0 0 1 0 0 0 297 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 298 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 299 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 300 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 301 324047 1 30 4 60798 1 0 0 0 0 0 0 0 0 0 0 302 311464 0 290 2 58790 0 1 0 0 0 0 0 0 0 0 0 303 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 304 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 305 353417 0 366 53 61808 0 0 0 0 1 0 0 0 0 0 0 306 325590 0 2 9 60735 0 0 0 0 0 1 0 0 0 0 0 307 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 308 328576 1 209 13 64016 0 0 0 0 0 0 0 1 0 0 0 309 326126 0 384 22 54683 0 0 0 0 0 0 0 0 1 0 0 310 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 311 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 312 369376 0 365 83 87192 0 0 0 0 0 0 0 0 0 0 0 313 325560 1 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 314 332013 1 49 8 64107 0 1 0 0 0 0 0 0 0 0 0 315 325871 0 3 4 60761 0 0 1 0 0 0 0 0 0 0 0 316 342165 0 133 14 65990 0 0 0 1 0 0 0 0 0 0 0 317 324967 0 32 1 59988 0 0 0 0 1 0 0 0 0 0 0 318 314832 0 368 17 61167 0 0 0 0 0 1 0 0 0 0 0 319 325557 0 1 6 60719 0 0 0 0 0 0 1 0 0 0 0 320 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 321 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 322 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 323 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 324 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 325 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 326 322649 0 22 2 60722 0 1 0 0 0 0 0 0 0 0 0 327 325560 1 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 328 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 329 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0 330 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 331 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 332 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 333 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 334 324598 0 96 5 60379 0 0 0 0 0 0 0 0 0 1 0 335 325567 0 1 2 60727 0 0 0 0 0 0 0 0 0 0 1 336 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 337 324005 0 81 7 60925 1 0 0 0 0 0 0 0 0 0 0 338 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 339 325748 0 26 1 60896 0 0 1 0 0 0 0 0 0 0 0 340 323385 0 125 13 59734 0 0 0 1 0 0 0 0 0 0 0 341 315409 0 304 15 62969 0 0 0 0 1 0 0 0 0 0 0 342 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 343 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 344 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 345 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 346 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 347 312275 0 119 6 59118 0 0 0 0 0 0 0 0 0 0 1 348 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 349 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0 350 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 351 320576 0 312 14 58598 0 0 1 0 0 0 0 0 0 0 0 352 325246 0 60 10 61124 0 0 0 1 0 0 0 0 0 0 0 353 332961 0 587 12 59595 0 0 0 0 1 0 0 0 0 0 0 354 323010 0 135 2 62065 0 0 0 0 0 1 0 0 0 0 0 355 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 356 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 357 345253 0 514 52 78780 0 0 0 0 0 0 0 0 1 0 0 358 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 359 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 360 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 361 325559 0 1 4 60722 1 0 0 0 0 0 0 0 0 0 0 362 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0 363 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 364 319634 1 58 3 61600 0 0 0 1 0 0 0 0 0 0 0 365 319951 0 180 11 59635 0 0 0 0 1 0 0 0 0 0 0 366 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0 367 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 368 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 369 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 370 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0 371 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1 372 318519 0 448 40 59781 0 0 0 0 0 0 0 0 0 0 0 373 343222 0 227 9 76644 1 0 0 0 0 0 0 0 0 0 0 374 317234 0 174 1 64820 0 1 0 0 0 0 0 0 0 0 0 375 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0 376 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0 377 314025 0 121 24 56178 0 0 0 0 1 0 0 0 0 0 0 378 320249 0 607 11 60436 0 0 0 0 0 1 0 0 0 0 0 379 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0 380 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0 381 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0 382 349365 0 530 60 73433 0 0 0 0 0 0 0 0 0 1 0 383 289197 0 571 80 41477 0 0 0 0 0 0 0 0 0 0 1 384 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0 385 329245 0 78 16 62700 1 0 0 0 0 0 0 0 0 0 0 386 240869 0 2489 40 67804 0 1 0 0 0 0 0 0 0 0 0 387 327182 0 131 6 59661 0 0 1 0 0 0 0 0 0 0 0 388 322876 0 923 8 58620 0 0 0 1 0 0 0 0 0 0 0 389 323117 0 72 3 60398 0 0 0 0 1 0 0 0 0 0 0 390 306351 0 572 16 58580 0 0 0 0 0 1 0 0 0 0 0 391 335137 0 397 10 62710 0 0 0 0 0 0 1 0 0 0 0 392 308271 0 450 8 59325 0 0 0 0 0 0 0 1 0 0 0 393 301731 0 622 7 60950 0 0 0 0 0 0 0 0 1 0 0 394 382409 0 694 8 68060 0 0 0 0 0 0 0 0 0 1 0 395 279230 1 3425 12 83620 0 0 0 0 0 0 0 0 0 0 1 396 298731 0 562 13 58456 0 0 0 0 0 0 0 0 0 0 0 397 243650 0 4917 42 52811 1 0 0 0 0 0 0 0 0 0 0 398 532682 1 1442 118 121173 0 1 0 0 0 0 0 0 0 0 0 399 319771 0 529 9 63870 0 0 1 0 0 0 0 0 0 0 0 400 171493 1 2126 138 21001 0 0 0 1 0 0 0 0 0 0 0 401 347262 0 1061 5 70415 0 0 0 0 1 0 0 0 0 0 0 402 343945 0 776 9 64230 0 0 0 0 0 1 0 0 0 0 0 403 311874 0 611 8 59190 0 0 0 0 0 0 1 0 0 0 0 404 302211 1 1526 25 69351 0 0 0 0 0 0 0 1 0 0 0 405 316708 0 592 7 64270 0 0 0 0 0 0 0 0 1 0 0 406 333463 0 1182 13 70694 0 0 0 0 0 0 0 0 0 1 0 407 344282 0 621 16 68005 0 0 0 0 0 0 0 0 0 0 1 408 319635 0 989 11 58930 0 0 0 0 0 0 0 0 0 0 0 409 301186 0 438 11 58320 1 0 0 0 0 0 0 0 0 0 0 410 300381 0 726 3 69980 0 1 0 0 0 0 0 0 0 0 0 411 318765 0 1303 61 69863 0 0 1 0 0 0 0 0 0 0 0 412 286146 1 6341 24 63255 0 0 0 1 0 0 0 0 0 0 0 413 306844 1 1164 17 57320 0 0 0 0 1 0 0 0 0 0 0 414 307705 1 3310 33 75230 0 0 0 0 0 1 0 0 0 0 0 415 312448 0 1366 7 79420 0 0 0 0 0 0 1 0 0 0 0 416 299715 0 965 3 73490 0 0 0 0 0 0 0 1 0 0 0 417 373399 0 3256 66 35250 0 0 0 0 0 0 0 0 1 0 0 418 299446 1 1135 17 62285 0 0 0 0 0 0 0 0 0 1 0 419 325586 0 1270 26 69206 0 0 0 0 0 0 0 0 0 0 1 420 291221 0 661 3 65920 0 0 0 0 0 0 0 0 0 0 0 421 261173 0 1013 2 69770 1 0 0 0 0 0 0 0 0 0 0 422 255027 0 2844 67 72683 0 1 0 0 0 0 0 0 0 0 0 423 -78375 1 11528 70 -14545 0 0 1 0 0 0 0 0 0 0 0 424 -58143 0 6526 26 55830 0 0 0 1 0 0 0 0 0 0 0 425 227033 0 2264 24 55174 0 0 0 0 1 0 0 0 0 0 0 426 235098 1 4461 94 67038 0 0 0 0 0 1 0 0 0 0 0 427 21267 0 3999 30 51252 0 0 0 0 0 0 1 0 0 0 0 428 238675 0 35624 223 157278 0 0 0 0 0 0 0 1 0 0 0 429 197687 0 9252 48 79510 0 0 0 0 0 0 0 0 1 0 0 430 418341 0 15236 90 77440 0 0 0 0 0 0 0 0 0 1 0 431 -297706 0 18073 180 27284 0 0 0 0 0 0 0 0 0 0 1 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades `Dividends\r` 7.218e+04 7.870e+04 2.628e+01 -4.662e+02 3.513e+00 M1 M2 M3 M4 M5 4.531e+04 3.026e+04 1.364e+04 -1.349e+05 -3.977e+04 M6 M7 M8 M9 M10 -2.718e+04 -1.341e+04 -3.472e+04 1.513e+04 -1.311e+04 M11 -6.342e+04 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -3011398 -45709 24947 67249 3216805 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 7.218e+04 5.942e+04 1.215 0.2252 Group 7.870e+04 3.328e+04 2.365 0.0185 * Costs 2.628e+01 1.933e+00 13.600 < 2e-16 *** Trades -4.662e+02 2.087e+02 -2.234 0.0260 * `Dividends\r` 3.513e+00 4.921e-01 7.138 4.25e-12 *** M1 4.531e+04 7.152e+04 0.634 0.5267 M2 3.026e+04 7.117e+04 0.425 0.6710 M3 1.364e+04 7.138e+04 0.191 0.8486 M4 -1.349e+05 7.144e+04 -1.889 0.0596 . M5 -3.977e+04 7.143e+04 -0.557 0.5780 M6 -2.718e+04 7.137e+04 -0.381 0.7035 M7 -1.341e+04 7.119e+04 -0.188 0.8507 M8 -3.472e+04 7.142e+04 -0.486 0.6271 M9 1.513e+04 7.115e+04 0.213 0.8318 M10 -1.311e+04 7.157e+04 -0.183 0.8547 M11 -6.342e+04 7.166e+04 -0.885 0.3767 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 299700 on 415 degrees of freedom Multiple R-squared: 0.567, Adjusted R-squared: 0.5514 F-statistic: 36.24 on 15 and 415 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 2.941805e-69 1.470903e-69 [2,] 1.0000000 2.801495e-106 1.400747e-106 [3,] 1.0000000 1.401116e-110 7.005578e-111 [4,] 1.0000000 2.975697e-120 1.487849e-120 [5,] 1.0000000 2.954009e-123 1.477004e-123 [6,] 1.0000000 1.568911e-135 7.844553e-136 [7,] 1.0000000 8.413709e-137 4.206855e-137 [8,] 1.0000000 1.099167e-137 5.495837e-138 [9,] 1.0000000 3.869193e-143 1.934596e-143 [10,] 1.0000000 2.344841e-147 1.172420e-147 [11,] 1.0000000 1.278893e-191 6.394467e-192 [12,] 1.0000000 1.577066e-202 7.885331e-203 [13,] 1.0000000 1.943662e-207 9.718312e-208 [14,] 1.0000000 2.375517e-210 1.187758e-210 [15,] 1.0000000 2.330796e-216 1.165398e-216 [16,] 1.0000000 2.256672e-218 1.128336e-218 [17,] 1.0000000 3.221510e-218 1.610755e-218 [18,] 1.0000000 4.898575e-218 2.449288e-218 [19,] 1.0000000 4.742914e-225 2.371457e-225 [20,] 1.0000000 1.948744e-236 9.743721e-237 [21,] 1.0000000 1.092052e-237 5.460258e-238 [22,] 1.0000000 2.275289e-238 1.137645e-238 [23,] 1.0000000 3.015279e-238 1.507639e-238 [24,] 1.0000000 1.051123e-238 5.255616e-239 [25,] 1.0000000 1.294814e-239 6.474068e-240 [26,] 1.0000000 6.679025e-239 3.339513e-239 [27,] 1.0000000 2.201441e-239 1.100721e-239 [28,] 1.0000000 4.107364e-239 2.053682e-239 [29,] 1.0000000 1.101589e-238 5.507944e-239 [30,] 1.0000000 2.793677e-239 1.396839e-239 [31,] 1.0000000 5.256679e-239 2.628339e-239 [32,] 1.0000000 4.475597e-239 2.237798e-239 [33,] 1.0000000 6.485537e-240 3.242769e-240 [34,] 1.0000000 3.921773e-239 1.960887e-239 [35,] 1.0000000 1.869029e-239 9.345144e-240 [36,] 1.0000000 1.153079e-238 5.765393e-239 [37,] 1.0000000 2.102602e-241 1.051301e-241 [38,] 1.0000000 4.328715e-242 2.164357e-242 [39,] 1.0000000 4.340002e-247 2.170001e-247 [40,] 1.0000000 2.925833e-247 1.462917e-247 [41,] 1.0000000 6.758724e-247 3.379362e-247 [42,] 1.0000000 8.012793e-247 4.006396e-247 [43,] 1.0000000 4.344185e-249 2.172093e-249 [44,] 1.0000000 2.003691e-248 1.001846e-248 [45,] 1.0000000 1.628414e-247 8.142070e-248 [46,] 1.0000000 9.002596e-247 4.501298e-247 [47,] 1.0000000 6.586270e-247 3.293135e-247 [48,] 1.0000000 5.501797e-246 2.750898e-246 [49,] 1.0000000 2.022656e-246 1.011328e-246 [50,] 1.0000000 2.494422e-245 1.247211e-245 [51,] 1.0000000 3.006567e-244 1.503284e-244 [52,] 1.0000000 2.157902e-243 1.078951e-243 [53,] 1.0000000 1.233786e-242 6.168932e-243 [54,] 1.0000000 1.934497e-244 9.672486e-245 [55,] 1.0000000 2.381179e-244 1.190590e-244 [56,] 1.0000000 2.654891e-243 1.327445e-243 [57,] 1.0000000 1.529057e-242 7.645283e-243 [58,] 1.0000000 7.378846e-242 3.689423e-242 [59,] 1.0000000 2.118425e-241 1.059213e-241 [60,] 1.0000000 5.504914e-241 2.752457e-241 [61,] 1.0000000 3.285034e-240 1.642517e-240 [62,] 1.0000000 1.828522e-239 9.142609e-240 [63,] 1.0000000 1.742933e-238 8.714663e-239 [64,] 1.0000000 4.386384e-238 2.193192e-238 [65,] 1.0000000 6.644184e-238 3.322092e-238 [66,] 1.0000000 3.532648e-237 1.766324e-237 [67,] 1.0000000 4.481949e-236 2.240974e-236 [68,] 1.0000000 1.343845e-236 6.719223e-237 [69,] 1.0000000 1.260994e-236 6.304970e-237 [70,] 1.0000000 2.049242e-236 1.024621e-236 [71,] 1.0000000 9.502894e-236 4.751447e-236 [72,] 1.0000000 1.004283e-234 5.021416e-235 [73,] 1.0000000 1.321405e-233 6.607025e-234 [74,] 1.0000000 1.487256e-232 7.436278e-233 [75,] 1.0000000 4.274861e-233 2.137430e-233 [76,] 1.0000000 9.845275e-233 4.922637e-233 [77,] 1.0000000 1.366558e-234 6.832788e-235 [78,] 1.0000000 2.551564e-234 1.275782e-234 [79,] 1.0000000 3.348561e-233 1.674280e-233 [80,] 1.0000000 1.629226e-232 8.146132e-233 [81,] 1.0000000 7.759503e-232 3.879751e-232 [82,] 1.0000000 7.580770e-231 3.790385e-231 [83,] 1.0000000 8.440279e-230 4.220139e-230 [84,] 1.0000000 4.101452e-230 2.050726e-230 [85,] 1.0000000 3.608508e-229 1.804254e-229 [86,] 1.0000000 6.447989e-229 3.223994e-229 [87,] 1.0000000 4.952683e-228 2.476341e-228 [88,] 1.0000000 5.543845e-228 2.771922e-228 [89,] 1.0000000 3.807794e-228 1.903897e-228 [90,] 1.0000000 4.145838e-228 2.072919e-228 [91,] 1.0000000 1.606167e-227 8.030835e-228 [92,] 1.0000000 3.877497e-227 1.938748e-227 [93,] 1.0000000 4.175775e-227 2.087887e-227 [94,] 1.0000000 3.528736e-227 1.764368e-227 [95,] 1.0000000 8.107497e-229 4.053748e-229 [96,] 1.0000000 3.521886e-228 1.760943e-228 [97,] 1.0000000 1.113767e-227 5.568834e-228 [98,] 1.0000000 3.047058e-227 1.523529e-227 [99,] 1.0000000 2.713747e-226 1.356874e-226 [100,] 1.0000000 9.420159e-226 4.710079e-226 [101,] 1.0000000 9.621586e-227 4.810793e-227 [102,] 1.0000000 9.092470e-226 4.546235e-226 [103,] 1.0000000 5.911602e-226 2.955801e-226 [104,] 1.0000000 6.783466e-225 3.391733e-225 [105,] 1.0000000 6.426902e-224 3.213451e-224 [106,] 1.0000000 3.500811e-223 1.750406e-223 [107,] 1.0000000 4.505707e-222 2.252854e-222 [108,] 1.0000000 1.212727e-222 6.063636e-223 [109,] 1.0000000 1.320725e-221 6.603624e-222 [110,] 1.0000000 1.450489e-220 7.252443e-221 [111,] 1.0000000 1.865639e-219 9.328197e-220 [112,] 1.0000000 1.933378e-218 9.666892e-219 [113,] 1.0000000 1.441552e-217 7.207758e-218 [114,] 1.0000000 1.102007e-217 5.510033e-218 [115,] 1.0000000 1.494501e-216 7.472504e-217 [116,] 1.0000000 1.991758e-215 9.958792e-216 [117,] 1.0000000 2.463684e-214 1.231842e-214 [118,] 1.0000000 7.972680e-214 3.986340e-214 [119,] 1.0000000 7.339706e-213 3.669853e-213 [120,] 1.0000000 7.351229e-212 3.675615e-212 [121,] 1.0000000 8.384224e-211 4.192112e-211 [122,] 1.0000000 9.745841e-210 4.872920e-210 [123,] 1.0000000 1.246152e-208 6.230759e-209 [124,] 1.0000000 1.359925e-207 6.799627e-208 [125,] 1.0000000 7.536415e-207 3.768207e-207 [126,] 1.0000000 8.592287e-206 4.296143e-206 [127,] 1.0000000 8.976044e-205 4.488022e-205 [128,] 1.0000000 9.182890e-204 4.591445e-204 [129,] 1.0000000 1.127261e-202 5.636306e-203 [130,] 1.0000000 5.467021e-202 2.733511e-202 [131,] 1.0000000 5.409526e-201 2.704763e-201 [132,] 1.0000000 5.700696e-200 2.850348e-200 [133,] 1.0000000 6.142554e-199 3.071277e-199 [134,] 1.0000000 5.338507e-198 2.669254e-198 [135,] 1.0000000 6.586177e-197 3.293088e-197 [136,] 1.0000000 6.702569e-196 3.351285e-196 [137,] 1.0000000 6.185776e-195 3.092888e-195 [138,] 1.0000000 6.927063e-194 3.463531e-194 [139,] 1.0000000 8.350915e-193 4.175457e-193 [140,] 1.0000000 9.971353e-192 4.985677e-192 [141,] 1.0000000 1.004645e-190 5.023224e-191 [142,] 1.0000000 5.770491e-190 2.885246e-190 [143,] 1.0000000 5.810473e-189 2.905237e-189 [144,] 1.0000000 6.123356e-188 3.061678e-188 [145,] 1.0000000 6.685978e-187 3.342989e-187 [146,] 1.0000000 6.976678e-186 3.488339e-186 [147,] 1.0000000 8.167824e-185 4.083912e-185 [148,] 1.0000000 8.585861e-184 4.292930e-184 [149,] 1.0000000 7.971005e-183 3.985502e-183 [150,] 1.0000000 7.126498e-184 3.563249e-184 [151,] 1.0000000 8.242743e-183 4.121372e-183 [152,] 1.0000000 9.466727e-182 4.733363e-182 [153,] 1.0000000 1.090390e-180 5.451948e-181 [154,] 1.0000000 7.231917e-180 3.615959e-180 [155,] 1.0000000 7.359436e-179 3.679718e-179 [156,] 1.0000000 7.753882e-178 3.876941e-178 [157,] 1.0000000 8.265849e-177 4.132925e-177 [158,] 1.0000000 9.096596e-176 4.548298e-176 [159,] 1.0000000 1.026884e-174 5.134419e-175 [160,] 1.0000000 1.104579e-173 5.522897e-174 [161,] 1.0000000 1.028828e-172 5.144139e-173 [162,] 1.0000000 1.105544e-171 5.527721e-172 [163,] 1.0000000 5.535408e-171 2.767704e-171 [164,] 1.0000000 5.543179e-170 2.771589e-170 [165,] 1.0000000 3.887074e-169 1.943537e-169 [166,] 1.0000000 2.199133e-168 1.099567e-168 [167,] 1.0000000 2.393133e-167 1.196567e-167 [168,] 1.0000000 2.506453e-166 1.253227e-166 [169,] 1.0000000 1.748126e-165 8.740628e-166 [170,] 1.0000000 1.847810e-164 9.239051e-165 [171,] 1.0000000 1.381292e-163 6.906462e-164 [172,] 1.0000000 1.354472e-162 6.772359e-163 [173,] 1.0000000 1.026893e-161 5.134467e-162 [174,] 1.0000000 1.062020e-160 5.310102e-161 [175,] 1.0000000 1.040761e-159 5.203804e-160 [176,] 1.0000000 1.077060e-158 5.385302e-159 [177,] 1.0000000 1.123061e-157 5.615303e-158 [178,] 1.0000000 1.020016e-156 5.100078e-157 [179,] 1.0000000 3.069498e-156 1.534749e-156 [180,] 1.0000000 3.201214e-155 1.600607e-155 [181,] 1.0000000 3.184626e-154 1.592313e-154 [182,] 1.0000000 7.871464e-154 3.935732e-154 [183,] 1.0000000 7.693796e-153 3.846898e-153 [184,] 1.0000000 7.319398e-152 3.659699e-152 [185,] 1.0000000 7.182237e-151 3.591118e-151 [186,] 1.0000000 6.130241e-150 3.065120e-150 [187,] 1.0000000 1.357700e-149 6.788502e-150 [188,] 1.0000000 1.307978e-148 6.539892e-149 [189,] 1.0000000 1.195505e-147 5.977524e-148 [190,] 1.0000000 1.063572e-146 5.317858e-147 [191,] 1.0000000 1.159517e-146 5.797584e-147 [192,] 1.0000000 1.171457e-145 5.857285e-146 [193,] 1.0000000 1.083696e-144 5.418481e-145 [194,] 1.0000000 1.059788e-143 5.298939e-144 [195,] 1.0000000 1.012633e-142 5.063164e-143 [196,] 1.0000000 8.286056e-142 4.143028e-142 [197,] 1.0000000 7.177807e-141 3.588904e-141 [198,] 1.0000000 6.869895e-140 3.434948e-140 [199,] 1.0000000 6.540345e-139 3.270172e-139 [200,] 1.0000000 5.612609e-138 2.806305e-138 [201,] 1.0000000 5.398515e-137 2.699258e-137 [202,] 1.0000000 1.834842e-136 9.174211e-137 [203,] 1.0000000 1.436863e-135 7.184317e-136 [204,] 1.0000000 1.311197e-134 6.555985e-135 [205,] 1.0000000 1.154118e-133 5.770588e-134 [206,] 1.0000000 1.080232e-132 5.401160e-133 [207,] 1.0000000 1.018488e-131 5.092438e-132 [208,] 1.0000000 9.175052e-131 4.587526e-131 [209,] 1.0000000 7.495193e-130 3.747597e-130 [210,] 1.0000000 6.852965e-129 3.426483e-129 [211,] 1.0000000 6.101010e-128 3.050505e-128 [212,] 1.0000000 5.397793e-127 2.698896e-127 [213,] 1.0000000 4.868744e-126 2.434372e-126 [214,] 1.0000000 3.366144e-125 1.683072e-125 [215,] 1.0000000 2.892967e-124 1.446483e-124 [216,] 1.0000000 2.436942e-123 1.218471e-123 [217,] 1.0000000 2.035962e-122 1.017981e-122 [218,] 1.0000000 1.809788e-121 9.048938e-122 [219,] 1.0000000 1.613978e-120 8.069888e-121 [220,] 1.0000000 1.360039e-119 6.800194e-120 [221,] 1.0000000 1.058236e-118 5.291179e-119 [222,] 1.0000000 9.182995e-118 4.591498e-118 [223,] 1.0000000 7.908310e-117 3.954155e-117 [224,] 1.0000000 6.284929e-116 3.142465e-116 [225,] 1.0000000 5.444925e-115 2.722463e-115 [226,] 1.0000000 3.584225e-114 1.792113e-114 [227,] 1.0000000 3.070324e-113 1.535162e-113 [228,] 1.0000000 2.636516e-112 1.318258e-112 [229,] 1.0000000 2.153023e-111 1.076511e-111 [230,] 1.0000000 1.790512e-110 8.952560e-111 [231,] 1.0000000 1.465542e-109 7.327712e-110 [232,] 1.0000000 1.124025e-108 5.620125e-109 [233,] 1.0000000 6.997735e-108 3.498868e-108 [234,] 1.0000000 5.817358e-107 2.908679e-107 [235,] 1.0000000 4.723700e-106 2.361850e-106 [236,] 1.0000000 3.790274e-105 1.895137e-105 [237,] 1.0000000 1.602432e-104 8.012161e-105 [238,] 1.0000000 1.059120e-103 5.295602e-104 [239,] 1.0000000 6.070422e-103 3.035211e-103 [240,] 1.0000000 2.846074e-102 1.423037e-102 [241,] 1.0000000 2.212079e-101 1.106040e-101 [242,] 1.0000000 1.780672e-100 8.903361e-101 [243,] 1.0000000 1.396310e-99 6.981548e-100 [244,] 1.0000000 1.015339e-98 5.076694e-99 [245,] 1.0000000 7.921415e-98 3.960707e-98 [246,] 1.0000000 5.658439e-97 2.829219e-97 [247,] 1.0000000 4.301282e-96 2.150641e-96 [248,] 1.0000000 3.194853e-95 1.597427e-95 [249,] 1.0000000 2.456480e-94 1.228240e-94 [250,] 1.0000000 1.541662e-93 7.708311e-94 [251,] 1.0000000 1.192863e-92 5.964313e-93 [252,] 1.0000000 8.703443e-92 4.351722e-92 [253,] 1.0000000 6.393698e-91 3.196849e-91 [254,] 1.0000000 4.846484e-90 2.423242e-90 [255,] 1.0000000 3.583137e-89 1.791569e-89 [256,] 1.0000000 2.476972e-88 1.238486e-88 [257,] 1.0000000 1.633560e-87 8.167800e-88 [258,] 1.0000000 1.165308e-86 5.826541e-87 [259,] 1.0000000 8.273244e-86 4.136622e-86 [260,] 1.0000000 5.761573e-85 2.880787e-85 [261,] 1.0000000 3.746063e-84 1.873032e-84 [262,] 1.0000000 2.462775e-83 1.231387e-83 [263,] 1.0000000 1.792387e-82 8.961933e-83 [264,] 1.0000000 1.271505e-81 6.357527e-82 [265,] 1.0000000 6.760946e-81 3.380473e-81 [266,] 1.0000000 4.846771e-80 2.423386e-80 [267,] 1.0000000 2.669705e-79 1.334853e-79 [268,] 1.0000000 1.765201e-78 8.826003e-79 [269,] 1.0000000 1.176541e-77 5.882705e-78 [270,] 1.0000000 8.210523e-77 4.105262e-77 [271,] 1.0000000 5.607142e-76 2.803571e-76 [272,] 1.0000000 3.082963e-75 1.541482e-75 [273,] 1.0000000 2.104632e-74 1.052316e-74 [274,] 1.0000000 1.186672e-73 5.933358e-74 [275,] 1.0000000 7.315189e-73 3.657594e-73 [276,] 1.0000000 4.916931e-72 2.458466e-72 [277,] 1.0000000 2.996756e-71 1.498378e-71 [278,] 1.0000000 1.546938e-70 7.734689e-71 [279,] 1.0000000 1.032840e-69 5.164201e-70 [280,] 1.0000000 6.108002e-69 3.054001e-69 [281,] 1.0000000 3.577420e-68 1.788710e-68 [282,] 1.0000000 2.342832e-67 1.171416e-67 [283,] 1.0000000 1.520375e-66 7.601876e-67 [284,] 1.0000000 9.336960e-66 4.668480e-66 [285,] 1.0000000 6.100591e-65 3.050296e-65 [286,] 1.0000000 3.755335e-64 1.877668e-64 [287,] 1.0000000 1.960133e-63 9.800663e-64 [288,] 1.0000000 1.236220e-62 6.181099e-63 [289,] 1.0000000 7.419573e-62 3.709786e-62 [290,] 1.0000000 4.670125e-61 2.335063e-61 [291,] 1.0000000 2.697774e-60 1.348887e-60 [292,] 1.0000000 1.464117e-59 7.320585e-60 [293,] 1.0000000 8.575283e-59 4.287641e-59 [294,] 1.0000000 4.007466e-58 2.003733e-58 [295,] 1.0000000 2.430232e-57 1.215116e-57 [296,] 1.0000000 1.427726e-56 7.138629e-57 [297,] 1.0000000 8.727023e-56 4.363512e-56 [298,] 1.0000000 4.723454e-55 2.361727e-55 [299,] 1.0000000 2.804461e-54 1.402231e-54 [300,] 1.0000000 1.677392e-53 8.386958e-54 [301,] 1.0000000 8.851561e-53 4.425781e-53 [302,] 1.0000000 5.226507e-52 2.613254e-52 [303,] 1.0000000 3.065840e-51 1.532920e-51 [304,] 1.0000000 1.625513e-50 8.127563e-51 [305,] 1.0000000 8.152298e-50 4.076149e-50 [306,] 1.0000000 4.719095e-49 2.359548e-49 [307,] 1.0000000 2.643374e-48 1.321687e-48 [308,] 1.0000000 1.384213e-47 6.921064e-48 [309,] 1.0000000 7.756228e-47 3.878114e-47 [310,] 1.0000000 4.090284e-46 2.045142e-46 [311,] 1.0000000 2.262971e-45 1.131486e-45 [312,] 1.0000000 1.239189e-44 6.195947e-45 [313,] 1.0000000 5.981734e-44 2.990867e-44 [314,] 1.0000000 3.265608e-43 1.632804e-43 [315,] 1.0000000 1.770469e-42 8.852347e-43 [316,] 1.0000000 8.724730e-42 4.362365e-42 [317,] 1.0000000 3.971450e-41 1.985725e-41 [318,] 1.0000000 2.110905e-40 1.055452e-40 [319,] 1.0000000 1.091248e-39 5.456238e-40 [320,] 1.0000000 5.067294e-39 2.533647e-39 [321,] 1.0000000 2.645497e-38 1.322749e-38 [322,] 1.0000000 1.182055e-37 5.910277e-38 [323,] 1.0000000 5.953598e-37 2.976799e-37 [324,] 1.0000000 2.994441e-36 1.497220e-36 [325,] 1.0000000 1.308535e-35 6.542674e-36 [326,] 1.0000000 6.558209e-35 3.279105e-35 [327,] 1.0000000 3.265558e-34 1.632779e-34 [328,] 1.0000000 1.452124e-33 7.260619e-34 [329,] 1.0000000 6.231999e-33 3.115999e-33 [330,] 1.0000000 3.033429e-32 1.516714e-32 [331,] 1.0000000 1.419982e-31 7.099910e-32 [332,] 1.0000000 5.855368e-31 2.927684e-31 [333,] 1.0000000 2.796118e-30 1.398059e-30 [334,] 1.0000000 1.162091e-29 5.810454e-30 [335,] 1.0000000 5.128022e-29 2.564011e-29 [336,] 1.0000000 2.374154e-28 1.187077e-28 [337,] 1.0000000 9.256266e-28 4.628133e-28 [338,] 1.0000000 4.209962e-27 2.104981e-27 [339,] 1.0000000 1.362964e-26 6.814819e-27 [340,] 1.0000000 5.524688e-26 2.762344e-26 [341,] 1.0000000 1.976074e-25 9.880368e-26 [342,] 1.0000000 8.764058e-25 4.382029e-25 [343,] 1.0000000 3.708098e-24 1.854049e-24 [344,] 1.0000000 1.287340e-23 6.436700e-24 [345,] 1.0000000 5.602286e-23 2.801143e-23 [346,] 1.0000000 2.283529e-22 1.141765e-22 [347,] 1.0000000 9.657493e-22 4.828746e-22 [348,] 1.0000000 3.961051e-21 1.980525e-21 [349,] 1.0000000 1.350552e-20 6.752761e-21 [350,] 1.0000000 5.507475e-20 2.753738e-20 [351,] 1.0000000 2.214132e-19 1.107066e-19 [352,] 1.0000000 8.000809e-19 4.000405e-19 [353,] 1.0000000 2.408746e-18 1.204373e-18 [354,] 1.0000000 9.590255e-18 4.795127e-18 [355,] 1.0000000 3.662111e-17 1.831056e-17 [356,] 1.0000000 1.250622e-16 6.253109e-17 [357,] 1.0000000 4.819349e-16 2.409674e-16 [358,] 1.0000000 1.519787e-15 7.598933e-16 [359,] 1.0000000 5.647197e-15 2.823599e-15 [360,] 1.0000000 1.985017e-14 9.925087e-15 [361,] 1.0000000 5.721411e-14 2.860706e-14 [362,] 1.0000000 2.018040e-13 1.009020e-13 [363,] 1.0000000 7.126554e-13 3.563277e-13 [364,] 1.0000000 1.556350e-12 7.781749e-13 [365,] 1.0000000 2.807860e-12 1.403930e-12 [366,] 1.0000000 9.898953e-12 4.949476e-12 [367,] 1.0000000 3.416875e-11 1.708438e-11 [368,] 1.0000000 1.157801e-10 5.789003e-11 [369,] 1.0000000 3.915537e-10 1.957769e-10 [370,] 1.0000000 8.296856e-10 4.148428e-10 [371,] 1.0000000 2.705631e-09 1.352816e-09 [372,] 1.0000000 8.215896e-09 4.107948e-09 [373,] 1.0000000 1.847582e-08 9.237910e-09 [374,] 1.0000000 5.532005e-08 2.766003e-08 [375,] 0.9999999 1.614701e-07 8.073506e-08 [376,] 0.9999998 4.908308e-07 2.454154e-07 [377,] 0.9999993 1.421330e-06 7.106648e-07 [378,] 0.9999979 4.172868e-06 2.086434e-06 [379,] 0.9999946 1.071998e-05 5.359990e-06 [380,] 0.9999852 2.953232e-05 1.476616e-05 [381,] 0.9999598 8.031072e-05 4.015536e-05 [382,] 0.9999400 1.199593e-04 5.997965e-05 [383,] 0.9998402 3.195598e-04 1.597799e-04 [384,] 0.9996486 7.027633e-04 3.513816e-04 [385,] 0.9994059 1.188147e-03 5.940734e-04 [386,] 0.9985228 2.954388e-03 1.477194e-03 [387,] 0.9967882 6.423508e-03 3.211754e-03 [388,] 0.9932633 1.347344e-02 6.736722e-03 [389,] 0.9864646 2.707088e-02 1.353544e-02 [390,] 0.9710111 5.797781e-02 2.898890e-02 [391,] 0.9422911 1.154179e-01 5.770893e-02 [392,] 0.8841626 2.316748e-01 1.158374e-01 [393,] 0.8113009 3.773983e-01 1.886991e-01 [394,] 0.8270985 3.458031e-01 1.729015e-01 > postscript(file="/var/www/html/rcomp/tmp/13xff1291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/23xff1291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3eow01291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4eow01291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5eow01291306282.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 1569858.1477 3216804.5650 1318169.8718 -3011398.4292 -368870.5277 6 7 8 9 10 56132.3433 2088.9969 616587.3876 -245850.0480 760920.0958 11 12 13 14 15 -276567.2123 462075.4986 905234.7776 -1819604.9030 16860.5273 16 17 18 19 20 558579.8501 -391037.7520 -63811.2844 -152457.7990 -309261.7393 21 22 23 24 25 135480.0583 -253842.5020 -190611.9735 -504512.5871 -277146.2879 26 27 28 29 30 -171676.1145 -400257.9907 -321172.5938 429960.1254 36650.9527 31 32 33 34 35 169625.5040 -151503.2213 93830.9434 5325.1372 -294870.5359 36 37 38 39 40 -469237.4376 -5701.6948 18715.6838 -163183.1173 -30951.2419 41 42 43 44 45 43415.1163 -342227.2146 99832.9212 -168130.2910 -43157.5823 46 47 48 49 50 -243920.4219 -211537.0486 -197604.9207 -166326.6387 -203553.7208 51 52 53 54 55 -13020.2667 -16553.3606 -391448.4984 -47772.7110 -407078.5585 56 57 58 59 60 152861.1425 151174.2325 -21731.3504 25629.2641 -249681.1403 61 62 63 64 65 -158595.9187 -103861.8336 -147288.9587 -121221.0244 88306.9426 66 67 68 69 70 90107.6762 -95650.0029 -89941.3646 -91213.1364 -147913.3273 71 72 73 74 75 10949.3791 57506.3400 -188093.9212 -62163.5410 -152804.1801 76 77 78 79 80 -30854.7275 -70049.8374 -37018.9676 83844.8336 -47622.7565 81 82 83 84 85 20131.2434 -109659.4849 -73596.1474 73092.8848 -68100.5949 86 87 88 89 90 -165579.6510 -181.4139 78258.6304 -73327.3995 -10263.5817 91 92 93 94 95 -30223.2831 -92607.5185 104610.2651 -350911.9747 189394.8250 96 97 98 99 100 -46129.8832 -50431.2966 -7656.7252 -163892.6059 24718.5248 101 102 103 104 105 -78894.3948 -91097.9397 -114331.1095 -36644.2181 27870.6053 106 107 108 109 110 -47860.2078 -199900.1793 -25894.5879 -265545.6115 -206906.6017 111 112 113 114 115 85871.1995 -123771.6851 -245146.4930 -171532.3565 -18703.3689 116 117 118 119 120 97843.3432 -20602.5726 37320.8522 -198173.0011 61162.1655 121 122 123 124 125 -321611.0043 -2794.4570 -28581.9495 61586.0824 -105202.5728 126 127 128 129 130 -246889.1051 53482.9591 75085.1371 33180.3946 53183.0694 131 132 133 134 135 106685.9316 97823.1152 -5239.7914 7784.7346 26435.2879 136 137 138 139 140 177598.0709 82157.2142 67249.1672 53482.9591 74794.3078 141 142 143 144 145 24033.0410 53183.0694 98446.4628 40338.0453 1024.4114 146 147 148 149 150 21860.2380 26435.2879 176950.4688 79836.7782 66677.6384 151 152 153 154 155 44210.8872 64236.9466 24946.6516 30999.5221 103581.6414 156 157 158 159 160 40071.6901 -5239.7914 6640.6904 33635.8992 174997.5049 161 162 163 164 165 79836.7782 67249.1672 53482.9591 80570.6160 25606.0709 166 167 168 169 170 53183.0694 103494.3525 144000.5954 -5239.7914 9812.7590 171 172 173 174 175 26435.2879 171989.4257 79836.7782 64997.4719 53482.9591 176 177 178 179 180 74830.4635 24946.6516 -26761.3976 103494.3525 40071.6901 181 182 183 184 185 8870.4918 -92379.1625 52232.5495 34221.3870 391.3591 186 187 188 189 190 -18969.3233 -20035.0601 71459.9865 -14812.3441 53183.0694 191 192 193 194 195 108202.7436 40071.6901 -3080.7867 9812.7590 26159.6990 196 197 198 199 200 96294.5998 -135572.3399 -11453.7379 -25219.9460 -82122.4322 201 202 203 204 205 -53756.2535 -25519.8357 24791.4474 -10296.1637 -170578.7136 206 207 208 209 210 -70637.5945 -67514.7389 96294.5998 60339.1038 -11453.7379 211 212 213 214 215 53482.9591 74794.3078 40850.3406 30976.8117 84784.1493 216 217 218 219 220 40071.6901 -1697.9686 14093.9080 26435.2879 180079.5502 221 222 223 224 225 47071.4215 67249.1672 54389.8956 74794.3078 25256.9161 226 227 228 229 230 57993.3973 103494.3525 40071.6901 -31464.8866 9812.7590 231 232 233 234 235 5291.3828 174997.5049 79836.7782 68678.2748 53482.9591 236 237 238 239 240 75053.1947 24946.6516 56168.9458 103494.3525 40071.6901 241 242 243 244 245 -5239.7914 12699.7489 26435.2879 180228.4215 -7166.6580 246 247 248 249 250 -11453.7379 -25219.9460 66680.0605 -53756.2535 -25519.8357 251 252 253 254 255 110695.2132 -38631.2150 -28672.4846 -9733.9543 -24628.0072 256 257 258 259 260 174997.5049 103324.4741 107375.5412 -25219.9460 74794.3078 261 262 263 264 265 -53756.2535 -25519.8357 83708.2143 48712.1463 -83942.6965 266 267 268 269 270 9812.7590 -52267.6172 178111.3539 1133.8731 45606.4224 271 272 273 274 275 -25219.9460 -3908.5973 -53756.2535 51734.8475 103494.3525 276 277 278 279 280 49275.3029 1841.8182 9812.7590 -63545.2223 153959.8020 281 282 283 284 285 1133.8731 67249.1672 92330.4772 74794.3078 -11005.1296 286 287 288 289 290 53183.0694 24791.4474 39090.5666 -83942.6965 40682.4378 291 292 293 294 295 -52267.6172 184802.0237 52906.5441 67249.1672 57166.2433 296 297 298 299 300 30865.1386 24946.6516 -25519.8357 103494.3525 40071.6901 301 302 303 304 305 -84653.1602 -4192.3673 26435.2879 96294.5998 118963.7940 306 307 308 309 310 71370.1413 -25219.9460 -11902.9183 46885.7458 -25519.8357 311 312 313 314 315 24791.4474 19998.4788 -83942.6965 -71893.3472 28388.4039 316 317 318 319 320 176121.2446 81440.5043 53205.5005 56254.6782 74794.3078 321 322 323 324 325 24946.6516 53183.0694 103494.3525 -38631.2150 -5239.7914 326 327 328 329 330 7249.0362 -52267.6172 96294.5998 79836.7782 67249.1672 331 332 333 334 335 53482.9591 74794.3078 24946.6516 53227.2049 104382.9762 336 337 338 339 340 40071.6901 -6380.0132 9812.7590 25787.9362 179062.3190 341 342 343 344 345 60789.2793 67249.1672 53482.9591 74794.3078 24946.6516 346 347 348 349 350 53183.0694 95507.1026 40071.6901 -5239.7914 9812.7590 351 352 353 354 355 27233.4321 176349.8594 81357.5962 57358.7320 53482.9591 356 357 358 359 360 74794.3078 -8068.1788 53183.0694 103494.3525 40071.6901 361 362 363 364 365 -3409.1070 9812.7590 26435.2879 87151.6133 78437.3973 366 367 368 369 370 67249.1672 53482.9591 74794.3078 24946.6516 53183.0694 371 372 373 374 375 103494.3525 43205.1677 -45287.8150 -17023.1039 26435.2879 376 377 378 379 380 174997.5049 92267.5274 52111.7090 53482.9591 74794.3078 381 382 383 384 385 24946.6516 46373.5097 157024.3228 40071.6901 -3100.4391 386 387 388 389 390 -146528.7581 31132.1234 159162.8844 78031.4273 47984.8525 391 392 393 394 395 50297.8980 54309.1973 -12773.6816 69737.5795 -186404.4461 396 397 398 399 400 12487.0534 -169009.4846 -57017.7005 -126.2104 90226.5753 401 402 403 404 405 41927.2186 57105.4751 32843.7778 -86027.3038 -8671.2602 406 407 408 409 410 1044.3173 87763.6077 19571.2178 -27565.2565 -65577.8594 411 412 413 414 415 -18282.3546 -107485.7438 -28303.5090 -151887.3283 -57958.1541 416 417 418 419 420 -19874.0729 107460.5380 -79034.8951 52454.3350 -28508.0185 421 422 423 424 425 -127109.7247 -146251.7195 -462135.4596 -350917.6996 -47518.7003 426 427 428 429 430 -197525.0426 -308661.5966 -1183573.5658 -389711.3115 -271235.1059 431 -793373.1387 > postscript(file="/var/www/html/rcomp/tmp/6ofel1291306282.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 1569858.1477 NA 1 3216804.5650 1569858.1477 2 1318169.8718 3216804.5650 3 -3011398.4292 1318169.8718 4 -368870.5277 -3011398.4292 5 56132.3433 -368870.5277 6 2088.9969 56132.3433 7 616587.3876 2088.9969 8 -245850.0480 616587.3876 9 760920.0958 -245850.0480 10 -276567.2123 760920.0958 11 462075.4986 -276567.2123 12 905234.7776 462075.4986 13 -1819604.9030 905234.7776 14 16860.5273 -1819604.9030 15 558579.8501 16860.5273 16 -391037.7520 558579.8501 17 -63811.2844 -391037.7520 18 -152457.7990 -63811.2844 19 -309261.7393 -152457.7990 20 135480.0583 -309261.7393 21 -253842.5020 135480.0583 22 -190611.9735 -253842.5020 23 -504512.5871 -190611.9735 24 -277146.2879 -504512.5871 25 -171676.1145 -277146.2879 26 -400257.9907 -171676.1145 27 -321172.5938 -400257.9907 28 429960.1254 -321172.5938 29 36650.9527 429960.1254 30 169625.5040 36650.9527 31 -151503.2213 169625.5040 32 93830.9434 -151503.2213 33 5325.1372 93830.9434 34 -294870.5359 5325.1372 35 -469237.4376 -294870.5359 36 -5701.6948 -469237.4376 37 18715.6838 -5701.6948 38 -163183.1173 18715.6838 39 -30951.2419 -163183.1173 40 43415.1163 -30951.2419 41 -342227.2146 43415.1163 42 99832.9212 -342227.2146 43 -168130.2910 99832.9212 44 -43157.5823 -168130.2910 45 -243920.4219 -43157.5823 46 -211537.0486 -243920.4219 47 -197604.9207 -211537.0486 48 -166326.6387 -197604.9207 49 -203553.7208 -166326.6387 50 -13020.2667 -203553.7208 51 -16553.3606 -13020.2667 52 -391448.4984 -16553.3606 53 -47772.7110 -391448.4984 54 -407078.5585 -47772.7110 55 152861.1425 -407078.5585 56 151174.2325 152861.1425 57 -21731.3504 151174.2325 58 25629.2641 -21731.3504 59 -249681.1403 25629.2641 60 -158595.9187 -249681.1403 61 -103861.8336 -158595.9187 62 -147288.9587 -103861.8336 63 -121221.0244 -147288.9587 64 88306.9426 -121221.0244 65 90107.6762 88306.9426 66 -95650.0029 90107.6762 67 -89941.3646 -95650.0029 68 -91213.1364 -89941.3646 69 -147913.3273 -91213.1364 70 10949.3791 -147913.3273 71 57506.3400 10949.3791 72 -188093.9212 57506.3400 73 -62163.5410 -188093.9212 74 -152804.1801 -62163.5410 75 -30854.7275 -152804.1801 76 -70049.8374 -30854.7275 77 -37018.9676 -70049.8374 78 83844.8336 -37018.9676 79 -47622.7565 83844.8336 80 20131.2434 -47622.7565 81 -109659.4849 20131.2434 82 -73596.1474 -109659.4849 83 73092.8848 -73596.1474 84 -68100.5949 73092.8848 85 -165579.6510 -68100.5949 86 -181.4139 -165579.6510 87 78258.6304 -181.4139 88 -73327.3995 78258.6304 89 -10263.5817 -73327.3995 90 -30223.2831 -10263.5817 91 -92607.5185 -30223.2831 92 104610.2651 -92607.5185 93 -350911.9747 104610.2651 94 189394.8250 -350911.9747 95 -46129.8832 189394.8250 96 -50431.2966 -46129.8832 97 -7656.7252 -50431.2966 98 -163892.6059 -7656.7252 99 24718.5248 -163892.6059 100 -78894.3948 24718.5248 101 -91097.9397 -78894.3948 102 -114331.1095 -91097.9397 103 -36644.2181 -114331.1095 104 27870.6053 -36644.2181 105 -47860.2078 27870.6053 106 -199900.1793 -47860.2078 107 -25894.5879 -199900.1793 108 -265545.6115 -25894.5879 109 -206906.6017 -265545.6115 110 85871.1995 -206906.6017 111 -123771.6851 85871.1995 112 -245146.4930 -123771.6851 113 -171532.3565 -245146.4930 114 -18703.3689 -171532.3565 115 97843.3432 -18703.3689 116 -20602.5726 97843.3432 117 37320.8522 -20602.5726 118 -198173.0011 37320.8522 119 61162.1655 -198173.0011 120 -321611.0043 61162.1655 121 -2794.4570 -321611.0043 122 -28581.9495 -2794.4570 123 61586.0824 -28581.9495 124 -105202.5728 61586.0824 125 -246889.1051 -105202.5728 126 53482.9591 -246889.1051 127 75085.1371 53482.9591 128 33180.3946 75085.1371 129 53183.0694 33180.3946 130 106685.9316 53183.0694 131 97823.1152 106685.9316 132 -5239.7914 97823.1152 133 7784.7346 -5239.7914 134 26435.2879 7784.7346 135 177598.0709 26435.2879 136 82157.2142 177598.0709 137 67249.1672 82157.2142 138 53482.9591 67249.1672 139 74794.3078 53482.9591 140 24033.0410 74794.3078 141 53183.0694 24033.0410 142 98446.4628 53183.0694 143 40338.0453 98446.4628 144 1024.4114 40338.0453 145 21860.2380 1024.4114 146 26435.2879 21860.2380 147 176950.4688 26435.2879 148 79836.7782 176950.4688 149 66677.6384 79836.7782 150 44210.8872 66677.6384 151 64236.9466 44210.8872 152 24946.6516 64236.9466 153 30999.5221 24946.6516 154 103581.6414 30999.5221 155 40071.6901 103581.6414 156 -5239.7914 40071.6901 157 6640.6904 -5239.7914 158 33635.8992 6640.6904 159 174997.5049 33635.8992 160 79836.7782 174997.5049 161 67249.1672 79836.7782 162 53482.9591 67249.1672 163 80570.6160 53482.9591 164 25606.0709 80570.6160 165 53183.0694 25606.0709 166 103494.3525 53183.0694 167 144000.5954 103494.3525 168 -5239.7914 144000.5954 169 9812.7590 -5239.7914 170 26435.2879 9812.7590 171 171989.4257 26435.2879 172 79836.7782 171989.4257 173 64997.4719 79836.7782 174 53482.9591 64997.4719 175 74830.4635 53482.9591 176 24946.6516 74830.4635 177 -26761.3976 24946.6516 178 103494.3525 -26761.3976 179 40071.6901 103494.3525 180 8870.4918 40071.6901 181 -92379.1625 8870.4918 182 52232.5495 -92379.1625 183 34221.3870 52232.5495 184 391.3591 34221.3870 185 -18969.3233 391.3591 186 -20035.0601 -18969.3233 187 71459.9865 -20035.0601 188 -14812.3441 71459.9865 189 53183.0694 -14812.3441 190 108202.7436 53183.0694 191 40071.6901 108202.7436 192 -3080.7867 40071.6901 193 9812.7590 -3080.7867 194 26159.6990 9812.7590 195 96294.5998 26159.6990 196 -135572.3399 96294.5998 197 -11453.7379 -135572.3399 198 -25219.9460 -11453.7379 199 -82122.4322 -25219.9460 200 -53756.2535 -82122.4322 201 -25519.8357 -53756.2535 202 24791.4474 -25519.8357 203 -10296.1637 24791.4474 204 -170578.7136 -10296.1637 205 -70637.5945 -170578.7136 206 -67514.7389 -70637.5945 207 96294.5998 -67514.7389 208 60339.1038 96294.5998 209 -11453.7379 60339.1038 210 53482.9591 -11453.7379 211 74794.3078 53482.9591 212 40850.3406 74794.3078 213 30976.8117 40850.3406 214 84784.1493 30976.8117 215 40071.6901 84784.1493 216 -1697.9686 40071.6901 217 14093.9080 -1697.9686 218 26435.2879 14093.9080 219 180079.5502 26435.2879 220 47071.4215 180079.5502 221 67249.1672 47071.4215 222 54389.8956 67249.1672 223 74794.3078 54389.8956 224 25256.9161 74794.3078 225 57993.3973 25256.9161 226 103494.3525 57993.3973 227 40071.6901 103494.3525 228 -31464.8866 40071.6901 229 9812.7590 -31464.8866 230 5291.3828 9812.7590 231 174997.5049 5291.3828 232 79836.7782 174997.5049 233 68678.2748 79836.7782 234 53482.9591 68678.2748 235 75053.1947 53482.9591 236 24946.6516 75053.1947 237 56168.9458 24946.6516 238 103494.3525 56168.9458 239 40071.6901 103494.3525 240 -5239.7914 40071.6901 241 12699.7489 -5239.7914 242 26435.2879 12699.7489 243 180228.4215 26435.2879 244 -7166.6580 180228.4215 245 -11453.7379 -7166.6580 246 -25219.9460 -11453.7379 247 66680.0605 -25219.9460 248 -53756.2535 66680.0605 249 -25519.8357 -53756.2535 250 110695.2132 -25519.8357 251 -38631.2150 110695.2132 252 -28672.4846 -38631.2150 253 -9733.9543 -28672.4846 254 -24628.0072 -9733.9543 255 174997.5049 -24628.0072 256 103324.4741 174997.5049 257 107375.5412 103324.4741 258 -25219.9460 107375.5412 259 74794.3078 -25219.9460 260 -53756.2535 74794.3078 261 -25519.8357 -53756.2535 262 83708.2143 -25519.8357 263 48712.1463 83708.2143 264 -83942.6965 48712.1463 265 9812.7590 -83942.6965 266 -52267.6172 9812.7590 267 178111.3539 -52267.6172 268 1133.8731 178111.3539 269 45606.4224 1133.8731 270 -25219.9460 45606.4224 271 -3908.5973 -25219.9460 272 -53756.2535 -3908.5973 273 51734.8475 -53756.2535 274 103494.3525 51734.8475 275 49275.3029 103494.3525 276 1841.8182 49275.3029 277 9812.7590 1841.8182 278 -63545.2223 9812.7590 279 153959.8020 -63545.2223 280 1133.8731 153959.8020 281 67249.1672 1133.8731 282 92330.4772 67249.1672 283 74794.3078 92330.4772 284 -11005.1296 74794.3078 285 53183.0694 -11005.1296 286 24791.4474 53183.0694 287 39090.5666 24791.4474 288 -83942.6965 39090.5666 289 40682.4378 -83942.6965 290 -52267.6172 40682.4378 291 184802.0237 -52267.6172 292 52906.5441 184802.0237 293 67249.1672 52906.5441 294 57166.2433 67249.1672 295 30865.1386 57166.2433 296 24946.6516 30865.1386 297 -25519.8357 24946.6516 298 103494.3525 -25519.8357 299 40071.6901 103494.3525 300 -84653.1602 40071.6901 301 -4192.3673 -84653.1602 302 26435.2879 -4192.3673 303 96294.5998 26435.2879 304 118963.7940 96294.5998 305 71370.1413 118963.7940 306 -25219.9460 71370.1413 307 -11902.9183 -25219.9460 308 46885.7458 -11902.9183 309 -25519.8357 46885.7458 310 24791.4474 -25519.8357 311 19998.4788 24791.4474 312 -83942.6965 19998.4788 313 -71893.3472 -83942.6965 314 28388.4039 -71893.3472 315 176121.2446 28388.4039 316 81440.5043 176121.2446 317 53205.5005 81440.5043 318 56254.6782 53205.5005 319 74794.3078 56254.6782 320 24946.6516 74794.3078 321 53183.0694 24946.6516 322 103494.3525 53183.0694 323 -38631.2150 103494.3525 324 -5239.7914 -38631.2150 325 7249.0362 -5239.7914 326 -52267.6172 7249.0362 327 96294.5998 -52267.6172 328 79836.7782 96294.5998 329 67249.1672 79836.7782 330 53482.9591 67249.1672 331 74794.3078 53482.9591 332 24946.6516 74794.3078 333 53227.2049 24946.6516 334 104382.9762 53227.2049 335 40071.6901 104382.9762 336 -6380.0132 40071.6901 337 9812.7590 -6380.0132 338 25787.9362 9812.7590 339 179062.3190 25787.9362 340 60789.2793 179062.3190 341 67249.1672 60789.2793 342 53482.9591 67249.1672 343 74794.3078 53482.9591 344 24946.6516 74794.3078 345 53183.0694 24946.6516 346 95507.1026 53183.0694 347 40071.6901 95507.1026 348 -5239.7914 40071.6901 349 9812.7590 -5239.7914 350 27233.4321 9812.7590 351 176349.8594 27233.4321 352 81357.5962 176349.8594 353 57358.7320 81357.5962 354 53482.9591 57358.7320 355 74794.3078 53482.9591 356 -8068.1788 74794.3078 357 53183.0694 -8068.1788 358 103494.3525 53183.0694 359 40071.6901 103494.3525 360 -3409.1070 40071.6901 361 9812.7590 -3409.1070 362 26435.2879 9812.7590 363 87151.6133 26435.2879 364 78437.3973 87151.6133 365 67249.1672 78437.3973 366 53482.9591 67249.1672 367 74794.3078 53482.9591 368 24946.6516 74794.3078 369 53183.0694 24946.6516 370 103494.3525 53183.0694 371 43205.1677 103494.3525 372 -45287.8150 43205.1677 373 -17023.1039 -45287.8150 374 26435.2879 -17023.1039 375 174997.5049 26435.2879 376 92267.5274 174997.5049 377 52111.7090 92267.5274 378 53482.9591 52111.7090 379 74794.3078 53482.9591 380 24946.6516 74794.3078 381 46373.5097 24946.6516 382 157024.3228 46373.5097 383 40071.6901 157024.3228 384 -3100.4391 40071.6901 385 -146528.7581 -3100.4391 386 31132.1234 -146528.7581 387 159162.8844 31132.1234 388 78031.4273 159162.8844 389 47984.8525 78031.4273 390 50297.8980 47984.8525 391 54309.1973 50297.8980 392 -12773.6816 54309.1973 393 69737.5795 -12773.6816 394 -186404.4461 69737.5795 395 12487.0534 -186404.4461 396 -169009.4846 12487.0534 397 -57017.7005 -169009.4846 398 -126.2104 -57017.7005 399 90226.5753 -126.2104 400 41927.2186 90226.5753 401 57105.4751 41927.2186 402 32843.7778 57105.4751 403 -86027.3038 32843.7778 404 -8671.2602 -86027.3038 405 1044.3173 -8671.2602 406 87763.6077 1044.3173 407 19571.2178 87763.6077 408 -27565.2565 19571.2178 409 -65577.8594 -27565.2565 410 -18282.3546 -65577.8594 411 -107485.7438 -18282.3546 412 -28303.5090 -107485.7438 413 -151887.3283 -28303.5090 414 -57958.1541 -151887.3283 415 -19874.0729 -57958.1541 416 107460.5380 -19874.0729 417 -79034.8951 107460.5380 418 52454.3350 -79034.8951 419 -28508.0185 52454.3350 420 -127109.7247 -28508.0185 421 -146251.7195 -127109.7247 422 -462135.4596 -146251.7195 423 -350917.6996 -462135.4596 424 -47518.7003 -350917.6996 425 -197525.0426 -47518.7003 426 -308661.5966 -197525.0426 427 -1183573.5658 -308661.5966 428 -389711.3115 -1183573.5658 429 -271235.1059 -389711.3115 430 -793373.1387 -271235.1059 431 NA -793373.1387 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 3216804.5650 1569858.1477 [2,] 1318169.8718 3216804.5650 [3,] -3011398.4292 1318169.8718 [4,] -368870.5277 -3011398.4292 [5,] 56132.3433 -368870.5277 [6,] 2088.9969 56132.3433 [7,] 616587.3876 2088.9969 [8,] -245850.0480 616587.3876 [9,] 760920.0958 -245850.0480 [10,] -276567.2123 760920.0958 [11,] 462075.4986 -276567.2123 [12,] 905234.7776 462075.4986 [13,] -1819604.9030 905234.7776 [14,] 16860.5273 -1819604.9030 [15,] 558579.8501 16860.5273 [16,] -391037.7520 558579.8501 [17,] -63811.2844 -391037.7520 [18,] -152457.7990 -63811.2844 [19,] -309261.7393 -152457.7990 [20,] 135480.0583 -309261.7393 [21,] -253842.5020 135480.0583 [22,] -190611.9735 -253842.5020 [23,] -504512.5871 -190611.9735 [24,] -277146.2879 -504512.5871 [25,] -171676.1145 -277146.2879 [26,] -400257.9907 -171676.1145 [27,] -321172.5938 -400257.9907 [28,] 429960.1254 -321172.5938 [29,] 36650.9527 429960.1254 [30,] 169625.5040 36650.9527 [31,] -151503.2213 169625.5040 [32,] 93830.9434 -151503.2213 [33,] 5325.1372 93830.9434 [34,] -294870.5359 5325.1372 [35,] -469237.4376 -294870.5359 [36,] -5701.6948 -469237.4376 [37,] 18715.6838 -5701.6948 [38,] -163183.1173 18715.6838 [39,] -30951.2419 -163183.1173 [40,] 43415.1163 -30951.2419 [41,] -342227.2146 43415.1163 [42,] 99832.9212 -342227.2146 [43,] -168130.2910 99832.9212 [44,] -43157.5823 -168130.2910 [45,] -243920.4219 -43157.5823 [46,] -211537.0486 -243920.4219 [47,] -197604.9207 -211537.0486 [48,] -166326.6387 -197604.9207 [49,] -203553.7208 -166326.6387 [50,] -13020.2667 -203553.7208 [51,] -16553.3606 -13020.2667 [52,] -391448.4984 -16553.3606 [53,] -47772.7110 -391448.4984 [54,] -407078.5585 -47772.7110 [55,] 152861.1425 -407078.5585 [56,] 151174.2325 152861.1425 [57,] -21731.3504 151174.2325 [58,] 25629.2641 -21731.3504 [59,] -249681.1403 25629.2641 [60,] -158595.9187 -249681.1403 [61,] -103861.8336 -158595.9187 [62,] -147288.9587 -103861.8336 [63,] -121221.0244 -147288.9587 [64,] 88306.9426 -121221.0244 [65,] 90107.6762 88306.9426 [66,] -95650.0029 90107.6762 [67,] -89941.3646 -95650.0029 [68,] -91213.1364 -89941.3646 [69,] -147913.3273 -91213.1364 [70,] 10949.3791 -147913.3273 [71,] 57506.3400 10949.3791 [72,] -188093.9212 57506.3400 [73,] -62163.5410 -188093.9212 [74,] -152804.1801 -62163.5410 [75,] -30854.7275 -152804.1801 [76,] -70049.8374 -30854.7275 [77,] -37018.9676 -70049.8374 [78,] 83844.8336 -37018.9676 [79,] -47622.7565 83844.8336 [80,] 20131.2434 -47622.7565 [81,] -109659.4849 20131.2434 [82,] -73596.1474 -109659.4849 [83,] 73092.8848 -73596.1474 [84,] -68100.5949 73092.8848 [85,] -165579.6510 -68100.5949 [86,] -181.4139 -165579.6510 [87,] 78258.6304 -181.4139 [88,] -73327.3995 78258.6304 [89,] -10263.5817 -73327.3995 [90,] -30223.2831 -10263.5817 [91,] -92607.5185 -30223.2831 [92,] 104610.2651 -92607.5185 [93,] -350911.9747 104610.2651 [94,] 189394.8250 -350911.9747 [95,] -46129.8832 189394.8250 [96,] -50431.2966 -46129.8832 [97,] -7656.7252 -50431.2966 [98,] -163892.6059 -7656.7252 [99,] 24718.5248 -163892.6059 [100,] -78894.3948 24718.5248 [101,] -91097.9397 -78894.3948 [102,] -114331.1095 -91097.9397 [103,] -36644.2181 -114331.1095 [104,] 27870.6053 -36644.2181 [105,] -47860.2078 27870.6053 [106,] -199900.1793 -47860.2078 [107,] -25894.5879 -199900.1793 [108,] -265545.6115 -25894.5879 [109,] -206906.6017 -265545.6115 [110,] 85871.1995 -206906.6017 [111,] -123771.6851 85871.1995 [112,] -245146.4930 -123771.6851 [113,] -171532.3565 -245146.4930 [114,] -18703.3689 -171532.3565 [115,] 97843.3432 -18703.3689 [116,] -20602.5726 97843.3432 [117,] 37320.8522 -20602.5726 [118,] -198173.0011 37320.8522 [119,] 61162.1655 -198173.0011 [120,] -321611.0043 61162.1655 [121,] -2794.4570 -321611.0043 [122,] -28581.9495 -2794.4570 [123,] 61586.0824 -28581.9495 [124,] -105202.5728 61586.0824 [125,] -246889.1051 -105202.5728 [126,] 53482.9591 -246889.1051 [127,] 75085.1371 53482.9591 [128,] 33180.3946 75085.1371 [129,] 53183.0694 33180.3946 [130,] 106685.9316 53183.0694 [131,] 97823.1152 106685.9316 [132,] -5239.7914 97823.1152 [133,] 7784.7346 -5239.7914 [134,] 26435.2879 7784.7346 [135,] 177598.0709 26435.2879 [136,] 82157.2142 177598.0709 [137,] 67249.1672 82157.2142 [138,] 53482.9591 67249.1672 [139,] 74794.3078 53482.9591 [140,] 24033.0410 74794.3078 [141,] 53183.0694 24033.0410 [142,] 98446.4628 53183.0694 [143,] 40338.0453 98446.4628 [144,] 1024.4114 40338.0453 [145,] 21860.2380 1024.4114 [146,] 26435.2879 21860.2380 [147,] 176950.4688 26435.2879 [148,] 79836.7782 176950.4688 [149,] 66677.6384 79836.7782 [150,] 44210.8872 66677.6384 [151,] 64236.9466 44210.8872 [152,] 24946.6516 64236.9466 [153,] 30999.5221 24946.6516 [154,] 103581.6414 30999.5221 [155,] 40071.6901 103581.6414 [156,] -5239.7914 40071.6901 [157,] 6640.6904 -5239.7914 [158,] 33635.8992 6640.6904 [159,] 174997.5049 33635.8992 [160,] 79836.7782 174997.5049 [161,] 67249.1672 79836.7782 [162,] 53482.9591 67249.1672 [163,] 80570.6160 53482.9591 [164,] 25606.0709 80570.6160 [165,] 53183.0694 25606.0709 [166,] 103494.3525 53183.0694 [167,] 144000.5954 103494.3525 [168,] -5239.7914 144000.5954 [169,] 9812.7590 -5239.7914 [170,] 26435.2879 9812.7590 [171,] 171989.4257 26435.2879 [172,] 79836.7782 171989.4257 [173,] 64997.4719 79836.7782 [174,] 53482.9591 64997.4719 [175,] 74830.4635 53482.9591 [176,] 24946.6516 74830.4635 [177,] -26761.3976 24946.6516 [178,] 103494.3525 -26761.3976 [179,] 40071.6901 103494.3525 [180,] 8870.4918 40071.6901 [181,] -92379.1625 8870.4918 [182,] 52232.5495 -92379.1625 [183,] 34221.3870 52232.5495 [184,] 391.3591 34221.3870 [185,] -18969.3233 391.3591 [186,] -20035.0601 -18969.3233 [187,] 71459.9865 -20035.0601 [188,] -14812.3441 71459.9865 [189,] 53183.0694 -14812.3441 [190,] 108202.7436 53183.0694 [191,] 40071.6901 108202.7436 [192,] -3080.7867 40071.6901 [193,] 9812.7590 -3080.7867 [194,] 26159.6990 9812.7590 [195,] 96294.5998 26159.6990 [196,] -135572.3399 96294.5998 [197,] -11453.7379 -135572.3399 [198,] -25219.9460 -11453.7379 [199,] -82122.4322 -25219.9460 [200,] -53756.2535 -82122.4322 [201,] -25519.8357 -53756.2535 [202,] 24791.4474 -25519.8357 [203,] -10296.1637 24791.4474 [204,] -170578.7136 -10296.1637 [205,] -70637.5945 -170578.7136 [206,] -67514.7389 -70637.5945 [207,] 96294.5998 -67514.7389 [208,] 60339.1038 96294.5998 [209,] -11453.7379 60339.1038 [210,] 53482.9591 -11453.7379 [211,] 74794.3078 53482.9591 [212,] 40850.3406 74794.3078 [213,] 30976.8117 40850.3406 [214,] 84784.1493 30976.8117 [215,] 40071.6901 84784.1493 [216,] -1697.9686 40071.6901 [217,] 14093.9080 -1697.9686 [218,] 26435.2879 14093.9080 [219,] 180079.5502 26435.2879 [220,] 47071.4215 180079.5502 [221,] 67249.1672 47071.4215 [222,] 54389.8956 67249.1672 [223,] 74794.3078 54389.8956 [224,] 25256.9161 74794.3078 [225,] 57993.3973 25256.9161 [226,] 103494.3525 57993.3973 [227,] 40071.6901 103494.3525 [228,] -31464.8866 40071.6901 [229,] 9812.7590 -31464.8866 [230,] 5291.3828 9812.7590 [231,] 174997.5049 5291.3828 [232,] 79836.7782 174997.5049 [233,] 68678.2748 79836.7782 [234,] 53482.9591 68678.2748 [235,] 75053.1947 53482.9591 [236,] 24946.6516 75053.1947 [237,] 56168.9458 24946.6516 [238,] 103494.3525 56168.9458 [239,] 40071.6901 103494.3525 [240,] -5239.7914 40071.6901 [241,] 12699.7489 -5239.7914 [242,] 26435.2879 12699.7489 [243,] 180228.4215 26435.2879 [244,] -7166.6580 180228.4215 [245,] -11453.7379 -7166.6580 [246,] -25219.9460 -11453.7379 [247,] 66680.0605 -25219.9460 [248,] -53756.2535 66680.0605 [249,] -25519.8357 -53756.2535 [250,] 110695.2132 -25519.8357 [251,] -38631.2150 110695.2132 [252,] -28672.4846 -38631.2150 [253,] -9733.9543 -28672.4846 [254,] -24628.0072 -9733.9543 [255,] 174997.5049 -24628.0072 [256,] 103324.4741 174997.5049 [257,] 107375.5412 103324.4741 [258,] -25219.9460 107375.5412 [259,] 74794.3078 -25219.9460 [260,] -53756.2535 74794.3078 [261,] -25519.8357 -53756.2535 [262,] 83708.2143 -25519.8357 [263,] 48712.1463 83708.2143 [264,] -83942.6965 48712.1463 [265,] 9812.7590 -83942.6965 [266,] -52267.6172 9812.7590 [267,] 178111.3539 -52267.6172 [268,] 1133.8731 178111.3539 [269,] 45606.4224 1133.8731 [270,] -25219.9460 45606.4224 [271,] -3908.5973 -25219.9460 [272,] -53756.2535 -3908.5973 [273,] 51734.8475 -53756.2535 [274,] 103494.3525 51734.8475 [275,] 49275.3029 103494.3525 [276,] 1841.8182 49275.3029 [277,] 9812.7590 1841.8182 [278,] -63545.2223 9812.7590 [279,] 153959.8020 -63545.2223 [280,] 1133.8731 153959.8020 [281,] 67249.1672 1133.8731 [282,] 92330.4772 67249.1672 [283,] 74794.3078 92330.4772 [284,] -11005.1296 74794.3078 [285,] 53183.0694 -11005.1296 [286,] 24791.4474 53183.0694 [287,] 39090.5666 24791.4474 [288,] -83942.6965 39090.5666 [289,] 40682.4378 -83942.6965 [290,] -52267.6172 40682.4378 [291,] 184802.0237 -52267.6172 [292,] 52906.5441 184802.0237 [293,] 67249.1672 52906.5441 [294,] 57166.2433 67249.1672 [295,] 30865.1386 57166.2433 [296,] 24946.6516 30865.1386 [297,] -25519.8357 24946.6516 [298,] 103494.3525 -25519.8357 [299,] 40071.6901 103494.3525 [300,] -84653.1602 40071.6901 [301,] -4192.3673 -84653.1602 [302,] 26435.2879 -4192.3673 [303,] 96294.5998 26435.2879 [304,] 118963.7940 96294.5998 [305,] 71370.1413 118963.7940 [306,] -25219.9460 71370.1413 [307,] -11902.9183 -25219.9460 [308,] 46885.7458 -11902.9183 [309,] -25519.8357 46885.7458 [310,] 24791.4474 -25519.8357 [311,] 19998.4788 24791.4474 [312,] -83942.6965 19998.4788 [313,] -71893.3472 -83942.6965 [314,] 28388.4039 -71893.3472 [315,] 176121.2446 28388.4039 [316,] 81440.5043 176121.2446 [317,] 53205.5005 81440.5043 [318,] 56254.6782 53205.5005 [319,] 74794.3078 56254.6782 [320,] 24946.6516 74794.3078 [321,] 53183.0694 24946.6516 [322,] 103494.3525 53183.0694 [323,] -38631.2150 103494.3525 [324,] -5239.7914 -38631.2150 [325,] 7249.0362 -5239.7914 [326,] -52267.6172 7249.0362 [327,] 96294.5998 -52267.6172 [328,] 79836.7782 96294.5998 [329,] 67249.1672 79836.7782 [330,] 53482.9591 67249.1672 [331,] 74794.3078 53482.9591 [332,] 24946.6516 74794.3078 [333,] 53227.2049 24946.6516 [334,] 104382.9762 53227.2049 [335,] 40071.6901 104382.9762 [336,] -6380.0132 40071.6901 [337,] 9812.7590 -6380.0132 [338,] 25787.9362 9812.7590 [339,] 179062.3190 25787.9362 [340,] 60789.2793 179062.3190 [341,] 67249.1672 60789.2793 [342,] 53482.9591 67249.1672 [343,] 74794.3078 53482.9591 [344,] 24946.6516 74794.3078 [345,] 53183.0694 24946.6516 [346,] 95507.1026 53183.0694 [347,] 40071.6901 95507.1026 [348,] -5239.7914 40071.6901 [349,] 9812.7590 -5239.7914 [350,] 27233.4321 9812.7590 [351,] 176349.8594 27233.4321 [352,] 81357.5962 176349.8594 [353,] 57358.7320 81357.5962 [354,] 53482.9591 57358.7320 [355,] 74794.3078 53482.9591 [356,] -8068.1788 74794.3078 [357,] 53183.0694 -8068.1788 [358,] 103494.3525 53183.0694 [359,] 40071.6901 103494.3525 [360,] -3409.1070 40071.6901 [361,] 9812.7590 -3409.1070 [362,] 26435.2879 9812.7590 [363,] 87151.6133 26435.2879 [364,] 78437.3973 87151.6133 [365,] 67249.1672 78437.3973 [366,] 53482.9591 67249.1672 [367,] 74794.3078 53482.9591 [368,] 24946.6516 74794.3078 [369,] 53183.0694 24946.6516 [370,] 103494.3525 53183.0694 [371,] 43205.1677 103494.3525 [372,] -45287.8150 43205.1677 [373,] -17023.1039 -45287.8150 [374,] 26435.2879 -17023.1039 [375,] 174997.5049 26435.2879 [376,] 92267.5274 174997.5049 [377,] 52111.7090 92267.5274 [378,] 53482.9591 52111.7090 [379,] 74794.3078 53482.9591 [380,] 24946.6516 74794.3078 [381,] 46373.5097 24946.6516 [382,] 157024.3228 46373.5097 [383,] 40071.6901 157024.3228 [384,] -3100.4391 40071.6901 [385,] -146528.7581 -3100.4391 [386,] 31132.1234 -146528.7581 [387,] 159162.8844 31132.1234 [388,] 78031.4273 159162.8844 [389,] 47984.8525 78031.4273 [390,] 50297.8980 47984.8525 [391,] 54309.1973 50297.8980 [392,] -12773.6816 54309.1973 [393,] 69737.5795 -12773.6816 [394,] -186404.4461 69737.5795 [395,] 12487.0534 -186404.4461 [396,] -169009.4846 12487.0534 [397,] -57017.7005 -169009.4846 [398,] -126.2104 -57017.7005 [399,] 90226.5753 -126.2104 [400,] 41927.2186 90226.5753 [401,] 57105.4751 41927.2186 [402,] 32843.7778 57105.4751 [403,] -86027.3038 32843.7778 [404,] -8671.2602 -86027.3038 [405,] 1044.3173 -8671.2602 [406,] 87763.6077 1044.3173 [407,] 19571.2178 87763.6077 [408,] -27565.2565 19571.2178 [409,] -65577.8594 -27565.2565 [410,] -18282.3546 -65577.8594 [411,] -107485.7438 -18282.3546 [412,] -28303.5090 -107485.7438 [413,] -151887.3283 -28303.5090 [414,] -57958.1541 -151887.3283 [415,] -19874.0729 -57958.1541 [416,] 107460.5380 -19874.0729 [417,] -79034.8951 107460.5380 [418,] 52454.3350 -79034.8951 [419,] -28508.0185 52454.3350 [420,] -127109.7247 -28508.0185 [421,] -146251.7195 -127109.7247 [422,] -462135.4596 -146251.7195 [423,] -350917.6996 -462135.4596 [424,] -47518.7003 -350917.6996 [425,] -197525.0426 -47518.7003 [426,] -308661.5966 -197525.0426 [427,] -1183573.5658 -308661.5966 [428,] -389711.3115 -1183573.5658 [429,] -271235.1059 -389711.3115 [430,] -793373.1387 -271235.1059 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 3216804.5650 1569858.1477 2 1318169.8718 3216804.5650 3 -3011398.4292 1318169.8718 4 -368870.5277 -3011398.4292 5 56132.3433 -368870.5277 6 2088.9969 56132.3433 7 616587.3876 2088.9969 8 -245850.0480 616587.3876 9 760920.0958 -245850.0480 10 -276567.2123 760920.0958 11 462075.4986 -276567.2123 12 905234.7776 462075.4986 13 -1819604.9030 905234.7776 14 16860.5273 -1819604.9030 15 558579.8501 16860.5273 16 -391037.7520 558579.8501 17 -63811.2844 -391037.7520 18 -152457.7990 -63811.2844 19 -309261.7393 -152457.7990 20 135480.0583 -309261.7393 21 -253842.5020 135480.0583 22 -190611.9735 -253842.5020 23 -504512.5871 -190611.9735 24 -277146.2879 -504512.5871 25 -171676.1145 -277146.2879 26 -400257.9907 -171676.1145 27 -321172.5938 -400257.9907 28 429960.1254 -321172.5938 29 36650.9527 429960.1254 30 169625.5040 36650.9527 31 -151503.2213 169625.5040 32 93830.9434 -151503.2213 33 5325.1372 93830.9434 34 -294870.5359 5325.1372 35 -469237.4376 -294870.5359 36 -5701.6948 -469237.4376 37 18715.6838 -5701.6948 38 -163183.1173 18715.6838 39 -30951.2419 -163183.1173 40 43415.1163 -30951.2419 41 -342227.2146 43415.1163 42 99832.9212 -342227.2146 43 -168130.2910 99832.9212 44 -43157.5823 -168130.2910 45 -243920.4219 -43157.5823 46 -211537.0486 -243920.4219 47 -197604.9207 -211537.0486 48 -166326.6387 -197604.9207 49 -203553.7208 -166326.6387 50 -13020.2667 -203553.7208 51 -16553.3606 -13020.2667 52 -391448.4984 -16553.3606 53 -47772.7110 -391448.4984 54 -407078.5585 -47772.7110 55 152861.1425 -407078.5585 56 151174.2325 152861.1425 57 -21731.3504 151174.2325 58 25629.2641 -21731.3504 59 -249681.1403 25629.2641 60 -158595.9187 -249681.1403 61 -103861.8336 -158595.9187 62 -147288.9587 -103861.8336 63 -121221.0244 -147288.9587 64 88306.9426 -121221.0244 65 90107.6762 88306.9426 66 -95650.0029 90107.6762 67 -89941.3646 -95650.0029 68 -91213.1364 -89941.3646 69 -147913.3273 -91213.1364 70 10949.3791 -147913.3273 71 57506.3400 10949.3791 72 -188093.9212 57506.3400 73 -62163.5410 -188093.9212 74 -152804.1801 -62163.5410 75 -30854.7275 -152804.1801 76 -70049.8374 -30854.7275 77 -37018.9676 -70049.8374 78 83844.8336 -37018.9676 79 -47622.7565 83844.8336 80 20131.2434 -47622.7565 81 -109659.4849 20131.2434 82 -73596.1474 -109659.4849 83 73092.8848 -73596.1474 84 -68100.5949 73092.8848 85 -165579.6510 -68100.5949 86 -181.4139 -165579.6510 87 78258.6304 -181.4139 88 -73327.3995 78258.6304 89 -10263.5817 -73327.3995 90 -30223.2831 -10263.5817 91 -92607.5185 -30223.2831 92 104610.2651 -92607.5185 93 -350911.9747 104610.2651 94 189394.8250 -350911.9747 95 -46129.8832 189394.8250 96 -50431.2966 -46129.8832 97 -7656.7252 -50431.2966 98 -163892.6059 -7656.7252 99 24718.5248 -163892.6059 100 -78894.3948 24718.5248 101 -91097.9397 -78894.3948 102 -114331.1095 -91097.9397 103 -36644.2181 -114331.1095 104 27870.6053 -36644.2181 105 -47860.2078 27870.6053 106 -199900.1793 -47860.2078 107 -25894.5879 -199900.1793 108 -265545.6115 -25894.5879 109 -206906.6017 -265545.6115 110 85871.1995 -206906.6017 111 -123771.6851 85871.1995 112 -245146.4930 -123771.6851 113 -171532.3565 -245146.4930 114 -18703.3689 -171532.3565 115 97843.3432 -18703.3689 116 -20602.5726 97843.3432 117 37320.8522 -20602.5726 118 -198173.0011 37320.8522 119 61162.1655 -198173.0011 120 -321611.0043 61162.1655 121 -2794.4570 -321611.0043 122 -28581.9495 -2794.4570 123 61586.0824 -28581.9495 124 -105202.5728 61586.0824 125 -246889.1051 -105202.5728 126 53482.9591 -246889.1051 127 75085.1371 53482.9591 128 33180.3946 75085.1371 129 53183.0694 33180.3946 130 106685.9316 53183.0694 131 97823.1152 106685.9316 132 -5239.7914 97823.1152 133 7784.7346 -5239.7914 134 26435.2879 7784.7346 135 177598.0709 26435.2879 136 82157.2142 177598.0709 137 67249.1672 82157.2142 138 53482.9591 67249.1672 139 74794.3078 53482.9591 140 24033.0410 74794.3078 141 53183.0694 24033.0410 142 98446.4628 53183.0694 143 40338.0453 98446.4628 144 1024.4114 40338.0453 145 21860.2380 1024.4114 146 26435.2879 21860.2380 147 176950.4688 26435.2879 148 79836.7782 176950.4688 149 66677.6384 79836.7782 150 44210.8872 66677.6384 151 64236.9466 44210.8872 152 24946.6516 64236.9466 153 30999.5221 24946.6516 154 103581.6414 30999.5221 155 40071.6901 103581.6414 156 -5239.7914 40071.6901 157 6640.6904 -5239.7914 158 33635.8992 6640.6904 159 174997.5049 33635.8992 160 79836.7782 174997.5049 161 67249.1672 79836.7782 162 53482.9591 67249.1672 163 80570.6160 53482.9591 164 25606.0709 80570.6160 165 53183.0694 25606.0709 166 103494.3525 53183.0694 167 144000.5954 103494.3525 168 -5239.7914 144000.5954 169 9812.7590 -5239.7914 170 26435.2879 9812.7590 171 171989.4257 26435.2879 172 79836.7782 171989.4257 173 64997.4719 79836.7782 174 53482.9591 64997.4719 175 74830.4635 53482.9591 176 24946.6516 74830.4635 177 -26761.3976 24946.6516 178 103494.3525 -26761.3976 179 40071.6901 103494.3525 180 8870.4918 40071.6901 181 -92379.1625 8870.4918 182 52232.5495 -92379.1625 183 34221.3870 52232.5495 184 391.3591 34221.3870 185 -18969.3233 391.3591 186 -20035.0601 -18969.3233 187 71459.9865 -20035.0601 188 -14812.3441 71459.9865 189 53183.0694 -14812.3441 190 108202.7436 53183.0694 191 40071.6901 108202.7436 192 -3080.7867 40071.6901 193 9812.7590 -3080.7867 194 26159.6990 9812.7590 195 96294.5998 26159.6990 196 -135572.3399 96294.5998 197 -11453.7379 -135572.3399 198 -25219.9460 -11453.7379 199 -82122.4322 -25219.9460 200 -53756.2535 -82122.4322 201 -25519.8357 -53756.2535 202 24791.4474 -25519.8357 203 -10296.1637 24791.4474 204 -170578.7136 -10296.1637 205 -70637.5945 -170578.7136 206 -67514.7389 -70637.5945 207 96294.5998 -67514.7389 208 60339.1038 96294.5998 209 -11453.7379 60339.1038 210 53482.9591 -11453.7379 211 74794.3078 53482.9591 212 40850.3406 74794.3078 213 30976.8117 40850.3406 214 84784.1493 30976.8117 215 40071.6901 84784.1493 216 -1697.9686 40071.6901 217 14093.9080 -1697.9686 218 26435.2879 14093.9080 219 180079.5502 26435.2879 220 47071.4215 180079.5502 221 67249.1672 47071.4215 222 54389.8956 67249.1672 223 74794.3078 54389.8956 224 25256.9161 74794.3078 225 57993.3973 25256.9161 226 103494.3525 57993.3973 227 40071.6901 103494.3525 228 -31464.8866 40071.6901 229 9812.7590 -31464.8866 230 5291.3828 9812.7590 231 174997.5049 5291.3828 232 79836.7782 174997.5049 233 68678.2748 79836.7782 234 53482.9591 68678.2748 235 75053.1947 53482.9591 236 24946.6516 75053.1947 237 56168.9458 24946.6516 238 103494.3525 56168.9458 239 40071.6901 103494.3525 240 -5239.7914 40071.6901 241 12699.7489 -5239.7914 242 26435.2879 12699.7489 243 180228.4215 26435.2879 244 -7166.6580 180228.4215 245 -11453.7379 -7166.6580 246 -25219.9460 -11453.7379 247 66680.0605 -25219.9460 248 -53756.2535 66680.0605 249 -25519.8357 -53756.2535 250 110695.2132 -25519.8357 251 -38631.2150 110695.2132 252 -28672.4846 -38631.2150 253 -9733.9543 -28672.4846 254 -24628.0072 -9733.9543 255 174997.5049 -24628.0072 256 103324.4741 174997.5049 257 107375.5412 103324.4741 258 -25219.9460 107375.5412 259 74794.3078 -25219.9460 260 -53756.2535 74794.3078 261 -25519.8357 -53756.2535 262 83708.2143 -25519.8357 263 48712.1463 83708.2143 264 -83942.6965 48712.1463 265 9812.7590 -83942.6965 266 -52267.6172 9812.7590 267 178111.3539 -52267.6172 268 1133.8731 178111.3539 269 45606.4224 1133.8731 270 -25219.9460 45606.4224 271 -3908.5973 -25219.9460 272 -53756.2535 -3908.5973 273 51734.8475 -53756.2535 274 103494.3525 51734.8475 275 49275.3029 103494.3525 276 1841.8182 49275.3029 277 9812.7590 1841.8182 278 -63545.2223 9812.7590 279 153959.8020 -63545.2223 280 1133.8731 153959.8020 281 67249.1672 1133.8731 282 92330.4772 67249.1672 283 74794.3078 92330.4772 284 -11005.1296 74794.3078 285 53183.0694 -11005.1296 286 24791.4474 53183.0694 287 39090.5666 24791.4474 288 -83942.6965 39090.5666 289 40682.4378 -83942.6965 290 -52267.6172 40682.4378 291 184802.0237 -52267.6172 292 52906.5441 184802.0237 293 67249.1672 52906.5441 294 57166.2433 67249.1672 295 30865.1386 57166.2433 296 24946.6516 30865.1386 297 -25519.8357 24946.6516 298 103494.3525 -25519.8357 299 40071.6901 103494.3525 300 -84653.1602 40071.6901 301 -4192.3673 -84653.1602 302 26435.2879 -4192.3673 303 96294.5998 26435.2879 304 118963.7940 96294.5998 305 71370.1413 118963.7940 306 -25219.9460 71370.1413 307 -11902.9183 -25219.9460 308 46885.7458 -11902.9183 309 -25519.8357 46885.7458 310 24791.4474 -25519.8357 311 19998.4788 24791.4474 312 -83942.6965 19998.4788 313 -71893.3472 -83942.6965 314 28388.4039 -71893.3472 315 176121.2446 28388.4039 316 81440.5043 176121.2446 317 53205.5005 81440.5043 318 56254.6782 53205.5005 319 74794.3078 56254.6782 320 24946.6516 74794.3078 321 53183.0694 24946.6516 322 103494.3525 53183.0694 323 -38631.2150 103494.3525 324 -5239.7914 -38631.2150 325 7249.0362 -5239.7914 326 -52267.6172 7249.0362 327 96294.5998 -52267.6172 328 79836.7782 96294.5998 329 67249.1672 79836.7782 330 53482.9591 67249.1672 331 74794.3078 53482.9591 332 24946.6516 74794.3078 333 53227.2049 24946.6516 334 104382.9762 53227.2049 335 40071.6901 104382.9762 336 -6380.0132 40071.6901 337 9812.7590 -6380.0132 338 25787.9362 9812.7590 339 179062.3190 25787.9362 340 60789.2793 179062.3190 341 67249.1672 60789.2793 342 53482.9591 67249.1672 343 74794.3078 53482.9591 344 24946.6516 74794.3078 345 53183.0694 24946.6516 346 95507.1026 53183.0694 347 40071.6901 95507.1026 348 -5239.7914 40071.6901 349 9812.7590 -5239.7914 350 27233.4321 9812.7590 351 176349.8594 27233.4321 352 81357.5962 176349.8594 353 57358.7320 81357.5962 354 53482.9591 57358.7320 355 74794.3078 53482.9591 356 -8068.1788 74794.3078 357 53183.0694 -8068.1788 358 103494.3525 53183.0694 359 40071.6901 103494.3525 360 -3409.1070 40071.6901 361 9812.7590 -3409.1070 362 26435.2879 9812.7590 363 87151.6133 26435.2879 364 78437.3973 87151.6133 365 67249.1672 78437.3973 366 53482.9591 67249.1672 367 74794.3078 53482.9591 368 24946.6516 74794.3078 369 53183.0694 24946.6516 370 103494.3525 53183.0694 371 43205.1677 103494.3525 372 -45287.8150 43205.1677 373 -17023.1039 -45287.8150 374 26435.2879 -17023.1039 375 174997.5049 26435.2879 376 92267.5274 174997.5049 377 52111.7090 92267.5274 378 53482.9591 52111.7090 379 74794.3078 53482.9591 380 24946.6516 74794.3078 381 46373.5097 24946.6516 382 157024.3228 46373.5097 383 40071.6901 157024.3228 384 -3100.4391 40071.6901 385 -146528.7581 -3100.4391 386 31132.1234 -146528.7581 387 159162.8844 31132.1234 388 78031.4273 159162.8844 389 47984.8525 78031.4273 390 50297.8980 47984.8525 391 54309.1973 50297.8980 392 -12773.6816 54309.1973 393 69737.5795 -12773.6816 394 -186404.4461 69737.5795 395 12487.0534 -186404.4461 396 -169009.4846 12487.0534 397 -57017.7005 -169009.4846 398 -126.2104 -57017.7005 399 90226.5753 -126.2104 400 41927.2186 90226.5753 401 57105.4751 41927.2186 402 32843.7778 57105.4751 403 -86027.3038 32843.7778 404 -8671.2602 -86027.3038 405 1044.3173 -8671.2602 406 87763.6077 1044.3173 407 19571.2178 87763.6077 408 -27565.2565 19571.2178 409 -65577.8594 -27565.2565 410 -18282.3546 -65577.8594 411 -107485.7438 -18282.3546 412 -28303.5090 -107485.7438 413 -151887.3283 -28303.5090 414 -57958.1541 -151887.3283 415 -19874.0729 -57958.1541 416 107460.5380 -19874.0729 417 -79034.8951 107460.5380 418 52454.3350 -79034.8951 419 -28508.0185 52454.3350 420 -127109.7247 -28508.0185 421 -146251.7195 -127109.7247 422 -462135.4596 -146251.7195 423 -350917.6996 -462135.4596 424 -47518.7003 -350917.6996 425 -197525.0426 -47518.7003 426 -308661.5966 -197525.0426 427 -1183573.5658 -308661.5966 428 -389711.3115 -1183573.5658 429 -271235.1059 -389711.3115 430 -793373.1387 -271235.1059 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7ofel1291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/8zovo1291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/9zovo1291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10sgu91291306282.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11dgtf1291306282.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12oqa01291306282.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13c87t1291306282.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14506e1291306282.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/1590n21291306282.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/16nakt1291306282.tab") + } > > try(system("convert tmp/13xff1291306282.ps tmp/13xff1291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/23xff1291306282.ps tmp/23xff1291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/3eow01291306282.ps tmp/3eow01291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/4eow01291306282.ps tmp/4eow01291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/5eow01291306282.ps tmp/5eow01291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/6ofel1291306282.ps tmp/6ofel1291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/7ofel1291306282.ps tmp/7ofel1291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/8zovo1291306282.ps tmp/8zovo1291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/9zovo1291306282.ps tmp/9zovo1291306282.png",intern=TRUE)) character(0) > try(system("convert tmp/10sgu91291306282.ps tmp/10sgu91291306282.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.658 2.188 36.624