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 = 'Do not include Seasonal Dummies' > par1 = '1' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Wealth Group Costs Trades Dividends\r 1 6282929 1 162556 1081 213118 2 4324047 1 29790 309 81767 3 4108272 1 87550 458 153198 4 -1212617 0 84738 588 -26007 5 1485329 1 54660 299 126942 6 1779876 1 42634 156 157214 7 1367203 0 40949 481 129352 8 2519076 1 42312 323 234817 9 912684 1 37704 452 60448 10 1443586 1 16275 109 47818 11 1220017 0 25830 115 245546 12 984885 0 12679 110 48020 13 1457425 1 18014 239 -1710 14 -572920 0 43556 247 32648 15 929144 1 24524 497 95350 16 1151176 0 6532 103 151352 17 790090 0 7123 109 288170 18 774497 1 20813 502 114337 19 990576 1 37597 248 37884 20 454195 0 17821 373 122844 21 876607 1 12988 119 82340 22 711969 1 22330 84 79801 23 702380 0 13326 102 165548 24 264449 0 16189 295 116384 25 450033 0 7146 105 134028 26 541063 0 15824 64 63838 27 588864 1 26088 267 74996 28 -37216 0 11326 129 31080 29 783310 0 8568 37 32168 30 467359 0 14416 361 49857 31 688779 1 3369 28 87161 32 608419 1 11819 85 106113 33 696348 1 6620 44 80570 34 597793 1 4519 49 102129 35 821730 0 2220 22 301670 36 377934 0 18562 155 102313 37 651939 0 10327 91 88577 38 697458 1 5336 81 112477 39 700368 1 2365 79 191778 40 225986 0 4069 145 79804 41 348695 0 7710 816 128294 42 373683 0 13718 61 96448 43 501709 0 4525 226 93811 44 413743 0 6869 105 117520 45 379825 0 4628 62 69159 46 336260 1 3653 24 101792 47 636765 1 1265 26 210568 48 481231 1 7489 322 136996 49 469107 0 4901 84 121920 50 211928 0 2284 33 76403 51 563925 1 3160 108 108094 52 511939 1 4150 150 134759 53 521016 1 7285 115 188873 54 543856 1 1134 162 146216 55 329304 1 4658 158 156608 56 423262 0 2384 97 61348 57 509665 0 3748 9 50350 58 455881 0 5371 66 87720 59 367772 0 1285 107 99489 60 406339 1 9327 101 87419 61 493408 1 5565 47 94355 62 232942 0 1528 38 60326 63 416002 1 3122 34 94670 64 337430 1 7317 84 82425 65 361517 0 2675 79 59017 66 360962 0 13253 947 90829 67 235561 0 880 74 80791 68 408247 1 2053 53 100423 69 450296 0 1424 94 131116 70 418799 1 4036 63 100269 71 247405 1 3045 58 27330 72 378519 0 5119 49 39039 73 326638 0 1431 34 106885 74 328233 0 554 11 79285 75 386225 0 1975 35 118881 76 283662 1 1286 17 77623 77 370225 0 1012 47 114768 78 269236 0 810 43 74015 79 365732 0 1280 117 69465 80 420383 1 666 171 117869 81 345811 0 1380 26 60982 82 431809 1 4608 73 90131 83 418876 0 876 59 138971 84 297476 0 814 18 39625 85 416776 0 514 15 102725 86 357257 1 5692 72 64239 87 458343 0 3642 86 90262 88 388386 0 540 14 103960 89 358934 0 2099 64 106611 90 407560 0 567 11 103345 91 392558 0 2001 52 95551 92 373177 1 2949 41 82903 93 428370 0 2253 99 63593 94 369419 1 6533 75 126910 95 358649 0 1889 45 37527 96 376641 1 3055 43 60247 97 467427 0 272 8 112995 98 364885 1 1414 198 70184 99 436230 0 2564 22 130140 100 329118 1 1383 11 73221 101 317365 1 1261 33 76114 102 286849 0 975 23 90534 103 376685 0 3366 80 108479 104 407198 0 576 18 113761 105 377772 0 1306 28 68696 106 271483 0 746 23 71561 107 153661 1 3192 60 59831 108 513294 1 2045 20 97890 109 324881 0 5477 59 101481 110 264512 1 1932 36 72954 111 420968 0 936 30 67939 112 129302 1 3437 47 48022 113 191521 0 5131 71 86111 114 268673 1 2397 14 74020 115 353179 1 1389 9 57530 116 354624 0 1503 39 56364 117 363713 0 402 26 84990 118 456657 0 2239 21 88590 119 211742 1 2234 16 77200 120 338381 0 837 69 61262 121 418530 0 10579 92 110309 122 351483 0 875 14 67000 123 372928 0 1395 103 93099 124 485538 1 1659 29 107577 125 279268 1 2647 37 62920 126 219060 1 3294 23 75832 127 325560 0 0 0 60720 128 325314 0 94 7 60793 129 322046 0 422 28 57935 130 325560 0 0 0 60720 131 325599 0 34 8 60630 132 377028 0 1558 63 55637 133 325560 0 0 0 60720 134 323850 0 43 3 60887 135 325560 0 0 0 60720 136 331514 0 316 9 60505 137 325632 0 115 13 60945 138 325560 0 0 0 60720 139 325560 0 0 0 60720 140 325560 0 0 0 60720 141 322265 0 389 14 58990 142 325560 0 0 0 60720 143 325906 0 1002 15 56750 144 325985 0 36 3 60894 145 346145 0 460 15 63346 146 325898 0 309 11 56535 147 325560 0 0 0 60720 148 325356 0 9 6 60835 149 325560 0 0 0 60720 150 325930 0 14 1 61016 151 318020 0 520 10 58650 152 326389 0 1766 73 60438 153 325560 0 0 0 60720 154 302925 0 458 11 58625 155 325540 0 20 3 60938 156 325560 0 0 0 60720 157 325560 0 0 0 60720 158 326736 0 98 2 61490 159 340580 0 405 7 60845 160 325560 0 0 0 60720 161 325560 0 0 0 60720 162 325560 0 0 0 60720 163 325560 0 0 0 60720 164 331828 0 483 27 60830 165 323299 0 454 51 63261 166 325560 0 0 0 60720 167 325560 0 0 0 60720 168 387722 0 757 19 45689 169 325560 0 0 0 60720 170 325560 0 0 0 60720 171 325560 0 0 0 60720 172 324598 0 36 4 61564 173 325560 0 0 0 60720 174 328726 0 203 9 61938 175 325560 0 0 0 60720 176 325043 0 90 8 60951 177 325560 0 0 0 60720 178 325806 1 71 1 60745 179 325560 0 0 0 60720 180 325560 0 0 0 60720 181 387732 0 972 34 71642 182 349729 1 531 10 71641 183 332202 0 604 38 55792 184 305442 1 283 10 71873 185 329537 1 23 5 62555 186 327055 1 638 14 60370 187 356245 1 699 16 64873 188 328451 0 149 5 62041 189 307062 0 226 5 65745 190 325560 0 0 0 60720 191 331345 0 275 4 59500 192 325560 0 0 0 60720 193 331824 0 141 6 61630 194 325560 0 0 0 60720 195 325685 0 28 2 60890 196 325560 1 0 0 60720 197 404480 1 2566 80 113521 198 325560 1 0 0 60720 199 325560 1 0 0 60720 200 318314 1 472 20 80045 201 325560 1 0 0 60720 202 325560 1 0 0 60720 203 325560 1 0 0 60720 204 311807 1 203 27 50804 205 337724 1 496 17 87390 206 326431 1 10 2 61656 207 327556 1 63 4 65688 208 325560 1 0 0 60720 209 356850 1 1136 32 48522 210 325560 1 0 0 60720 211 325560 0 0 0 60720 212 325560 0 0 0 60720 213 322741 0 267 32 57640 214 310902 0 474 20 61977 215 324295 0 534 7 62620 216 325560 0 0 0 60720 217 326156 0 15 8 60831 218 326960 0 397 28 60646 219 325560 0 0 0 60720 220 333411 0 1061 20 56225 221 297761 0 288 4 60510 222 325560 0 0 0 60720 223 325536 0 3 2 60698 224 325560 0 0 0 60720 225 325762 0 20 2 60805 226 327957 0 278 26 61404 227 325560 0 0 0 60720 228 325560 0 0 0 60720 229 318521 0 192 4 65276 230 325560 0 0 0 60720 231 319775 0 317 9 63915 232 325560 0 0 0 60720 233 325560 0 0 0 60720 234 332128 0 368 17 61686 235 325560 0 0 0 60720 236 325486 0 2 1 60743 237 325560 0 0 0 60720 238 325838 0 53 6 60349 239 325560 0 0 0 60720 240 325560 0 0 0 60720 241 325560 0 0 0 60720 242 331767 0 94 3 61360 243 325560 0 0 0 60720 244 324523 0 24 8 59818 245 339995 0 2332 4 72680 246 325560 1 0 0 60720 247 325560 1 0 0 60720 248 319582 0 131 11 61808 249 325560 1 0 0 60720 250 325560 1 0 0 60720 251 307245 0 206 9 53110 252 325560 1 0 0 60720 253 317967 0 167 2 64245 254 331488 0 622 73 73007 255 335452 0 885 85 82732 256 325560 0 0 0 60720 257 334184 0 365 8 54820 258 313213 0 364 35 47705 259 325560 1 0 0 60720 260 325560 0 0 0 60720 261 325560 1 0 0 60720 262 325560 1 0 0 60720 263 348678 0 226 12 72835 264 328727 0 307 15 58856 265 325560 1 0 0 60720 266 325560 0 0 0 60720 267 325560 1 0 0 60720 268 387978 0 188 11 77655 269 325560 1 0 0 60720 270 336704 0 138 6 69817 271 325560 1 0 0 60720 272 325560 1 0 0 60720 273 325560 1 0 0 60720 274 322076 0 125 12 60798 275 325560 0 0 0 60720 276 334272 0 282 30 62452 277 338197 0 335 33 64175 278 325560 0 0 0 60720 279 321024 1 813 82 67440 280 322145 0 176 28 68136 281 325560 1 0 0 60720 282 325560 0 0 0 60720 283 323351 0 249 72 56726 284 325560 0 0 0 60720 285 327748 0 333 13 70811 286 325560 0 0 0 60720 287 325560 1 0 0 60720 288 328157 0 30 4 62045 289 325560 1 0 0 60720 290 311594 0 249 62 54323 291 325560 1 0 0 60720 292 335962 0 165 24 62841 293 372426 0 453 21 81125 294 325560 0 0 0 60720 295 319844 0 53 14 59506 296 355822 1 382 21 59365 297 325560 0 0 0 60720 298 325560 1 0 0 60720 299 325560 0 0 0 60720 300 325560 0 0 0 60720 301 324047 1 30 4 60798 302 311464 0 290 2 58790 303 325560 0 0 0 60720 304 325560 1 0 0 60720 305 353417 0 366 53 61808 306 325590 0 2 9 60735 307 325560 1 0 0 60720 308 328576 1 209 13 64016 309 326126 0 384 22 54683 310 325560 1 0 0 60720 311 325560 1 0 0 60720 312 369376 0 365 83 87192 313 325560 1 0 0 60720 314 332013 1 49 8 64107 315 325871 0 3 4 60761 316 342165 0 133 14 65990 317 324967 0 32 1 59988 318 314832 0 368 17 61167 319 325557 0 1 6 60719 320 325560 0 0 0 60720 321 325560 0 0 0 60720 322 325560 0 0 0 60720 323 325560 0 0 0 60720 324 325560 1 0 0 60720 325 325560 0 0 0 60720 326 322649 0 22 2 60722 327 325560 1 0 0 60720 328 325560 1 0 0 60720 329 325560 0 0 0 60720 330 325560 0 0 0 60720 331 325560 0 0 0 60720 332 325560 0 0 0 60720 333 325560 0 0 0 60720 334 324598 0 96 5 60379 335 325567 0 1 2 60727 336 325560 0 0 0 60720 337 324005 0 81 7 60925 338 325560 0 0 0 60720 339 325748 0 26 1 60896 340 323385 0 125 13 59734 341 315409 0 304 15 62969 342 325560 0 0 0 60720 343 325560 0 0 0 60720 344 325560 0 0 0 60720 345 325560 0 0 0 60720 346 325560 0 0 0 60720 347 312275 0 119 6 59118 348 325560 0 0 0 60720 349 325560 0 0 0 60720 350 325560 0 0 0 60720 351 320576 0 312 14 58598 352 325246 0 60 10 61124 353 332961 0 587 12 59595 354 323010 0 135 2 62065 355 325560 0 0 0 60720 356 325560 0 0 0 60720 357 345253 0 514 52 78780 358 325560 0 0 0 60720 359 325560 0 0 0 60720 360 325560 0 0 0 60720 361 325559 0 1 4 60722 362 325560 0 0 0 60720 363 325560 0 0 0 60720 364 319634 1 58 3 61600 365 319951 0 180 11 59635 366 325560 0 0 0 60720 367 325560 0 0 0 60720 368 325560 0 0 0 60720 369 325560 0 0 0 60720 370 325560 0 0 0 60720 371 325560 0 0 0 60720 372 318519 0 448 40 59781 373 343222 0 227 9 76644 374 317234 0 174 1 64820 375 325560 0 0 0 60720 376 325560 0 0 0 60720 377 314025 0 121 24 56178 378 320249 0 607 11 60436 379 325560 0 0 0 60720 380 325560 0 0 0 60720 381 325560 0 0 0 60720 382 349365 0 530 60 73433 383 289197 0 571 80 41477 384 325560 0 0 0 60720 385 329245 0 78 16 62700 386 240869 0 2489 40 67804 387 327182 0 131 6 59661 388 322876 0 923 8 58620 389 323117 0 72 3 60398 390 306351 0 572 16 58580 391 335137 0 397 10 62710 392 308271 0 450 8 59325 393 301731 0 622 7 60950 394 382409 0 694 8 68060 395 279230 1 3425 12 83620 396 298731 0 562 13 58456 397 243650 0 4917 42 52811 398 532682 1 1442 118 121173 399 319771 0 529 9 63870 400 171493 1 2126 138 21001 401 347262 0 1061 5 70415 402 343945 0 776 9 64230 403 311874 0 611 8 59190 404 302211 1 1526 25 69351 405 316708 0 592 7 64270 406 333463 0 1182 13 70694 407 344282 0 621 16 68005 408 319635 0 989 11 58930 409 301186 0 438 11 58320 410 300381 0 726 3 69980 411 318765 0 1303 61 69863 412 286146 1 6341 24 63255 413 306844 1 1164 17 57320 414 307705 1 3310 33 75230 415 312448 0 1366 7 79420 416 299715 0 965 3 73490 417 373399 0 3256 66 35250 418 299446 1 1135 17 62285 419 325586 0 1270 26 69206 420 291221 0 661 3 65920 421 261173 0 1013 2 69770 422 255027 0 2844 67 72683 423 -78375 1 11528 70 -14545 424 -58143 0 6526 26 55830 425 227033 0 2264 24 55174 426 235098 1 4461 94 67038 427 21267 0 3999 30 51252 428 238675 0 35624 223 157278 429 197687 0 9252 48 79510 430 418341 0 15236 90 77440 431 -297706 0 18073 180 27284 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades `Dividends\r` 56500.784 70846.109 26.482 -470.732 3.497 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -3145447 -20447 47198 56744 3267340 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 56500.7845 36533.7773 1.547 0.1227 Group 70846.1094 32881.3355 2.155 0.0318 * Costs 26.4823 1.9134 13.840 < 2e-16 *** Trades -470.7317 206.1013 -2.284 0.0229 * `Dividends\r` 3.4966 0.4833 7.235 2.18e-12 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 299400 on 426 degrees of freedom Multiple R-squared: 0.5565, Adjusted R-squared: 0.5523 F-statistic: 133.6 on 4 and 426 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 2.471451e-46 1.235726e-46 [2,] 1.0000000 4.814834e-79 2.407417e-79 [3,] 1.0000000 1.705622e-87 8.528108e-88 [4,] 1.0000000 2.115741e-89 1.057870e-89 [5,] 1.0000000 3.297289e-107 1.648645e-107 [6,] 1.0000000 3.034309e-124 1.517155e-124 [7,] 1.0000000 2.132831e-139 1.066416e-139 [8,] 1.0000000 5.694279e-151 2.847139e-151 [9,] 1.0000000 5.329587e-161 2.664794e-161 [10,] 1.0000000 2.283699e-162 1.141850e-162 [11,] 1.0000000 1.567413e-167 7.837067e-168 [12,] 1.0000000 2.995015e-177 1.497507e-177 [13,] 1.0000000 2.275749e-176 1.137875e-176 [14,] 1.0000000 1.479932e-183 7.399660e-184 [15,] 1.0000000 1.182463e-188 5.912316e-189 [16,] 1.0000000 9.793603e-189 4.896802e-189 [17,] 1.0000000 3.129126e-188 1.564563e-188 [18,] 1.0000000 2.304281e-187 1.152140e-187 [19,] 1.0000000 3.182174e-190 1.591087e-190 [20,] 1.0000000 2.314927e-195 1.157463e-195 [21,] 1.0000000 3.668875e-197 1.834437e-197 [22,] 1.0000000 2.909101e-214 1.454551e-214 [23,] 1.0000000 1.964745e-217 9.823726e-218 [24,] 1.0000000 1.014881e-221 5.074407e-222 [25,] 1.0000000 1.613932e-225 8.069660e-226 [26,] 1.0000000 1.852849e-232 9.264246e-233 [27,] 1.0000000 1.357399e-234 6.786993e-235 [28,] 1.0000000 1.918873e-234 9.594364e-235 [29,] 1.0000000 5.321350e-234 2.660675e-234 [30,] 1.0000000 4.075108e-242 2.037554e-242 [31,] 1.0000000 6.493496e-248 3.246748e-248 [32,] 1.0000000 4.203875e-249 2.101938e-249 [33,] 1.0000000 5.151715e-249 2.575858e-249 [34,] 1.0000000 6.188678e-249 3.094339e-249 [35,] 1.0000000 8.847765e-249 4.423883e-249 [36,] 1.0000000 2.189682e-249 1.094841e-249 [37,] 1.0000000 2.642796e-248 1.321398e-248 [38,] 1.0000000 3.267642e-248 1.633821e-248 [39,] 1.0000000 4.951756e-248 2.475878e-248 [40,] 1.0000000 1.665862e-248 8.329310e-249 [41,] 1.0000000 1.843172e-248 9.215858e-249 [42,] 1.0000000 1.218084e-247 6.090419e-248 [43,] 1.0000000 7.435821e-248 3.717911e-248 [44,] 1.0000000 6.502774e-249 3.251387e-249 [45,] 1.0000000 1.526112e-248 7.630560e-249 [46,] 1.0000000 6.527344e-249 3.263672e-249 [47,] 1.0000000 4.357798e-248 2.178899e-248 [48,] 1.0000000 1.944273e-250 9.721364e-251 [49,] 1.0000000 4.331229e-251 2.165614e-251 [50,] 1.0000000 1.638086e-256 8.190431e-257 [51,] 1.0000000 4.754339e-257 2.377169e-257 [52,] 1.0000000 3.496478e-256 1.748239e-256 [53,] 1.0000000 8.880228e-257 4.440114e-257 [54,] 1.0000000 4.096102e-258 2.048051e-258 [55,] 1.0000000 1.194414e-257 5.972071e-258 [56,] 1.0000000 7.649750e-257 3.824875e-257 [57,] 1.0000000 3.592465e-256 1.796233e-256 [58,] 1.0000000 8.856101e-256 4.428050e-256 [59,] 1.0000000 8.307795e-255 4.153898e-255 [60,] 1.0000000 2.358103e-255 1.179052e-255 [61,] 1.0000000 2.727739e-254 1.363870e-254 [62,] 1.0000000 3.077846e-253 1.538923e-253 [63,] 1.0000000 1.823895e-252 9.119477e-253 [64,] 1.0000000 1.887804e-251 9.439018e-252 [65,] 1.0000000 1.274720e-253 6.373598e-254 [66,] 1.0000000 2.184527e-253 1.092263e-253 [67,] 1.0000000 1.706563e-252 8.532816e-253 [68,] 1.0000000 1.049128e-251 5.245640e-252 [69,] 1.0000000 4.766407e-251 2.383203e-251 [70,] 1.0000000 1.301089e-250 6.505447e-251 [71,] 1.0000000 2.990510e-250 1.495255e-250 [72,] 1.0000000 2.050286e-249 1.025143e-249 [73,] 1.0000000 1.271508e-248 6.357538e-249 [74,] 1.0000000 7.826428e-248 3.913214e-248 [75,] 1.0000000 1.171875e-247 5.859377e-248 [76,] 1.0000000 9.946689e-248 4.973345e-248 [77,] 1.0000000 4.855968e-247 2.427984e-247 [78,] 1.0000000 5.363695e-246 2.681847e-246 [79,] 1.0000000 4.447679e-246 2.223839e-246 [80,] 1.0000000 3.891692e-246 1.945846e-246 [81,] 1.0000000 2.750131e-245 1.375066e-245 [82,] 1.0000000 1.183763e-244 5.918817e-245 [83,] 1.0000000 1.166529e-243 5.832646e-244 [84,] 1.0000000 1.450654e-242 7.253271e-243 [85,] 1.0000000 1.370590e-241 6.852949e-242 [86,] 1.0000000 1.833357e-242 9.166784e-243 [87,] 1.0000000 5.306553e-242 2.653277e-242 [88,] 1.0000000 4.338170e-243 2.169085e-243 [89,] 1.0000000 3.268875e-243 1.634438e-243 [90,] 1.0000000 3.915739e-242 1.957870e-242 [91,] 1.0000000 2.754674e-241 1.377337e-241 [92,] 1.0000000 1.327974e-240 6.639872e-241 [93,] 1.0000000 1.707697e-239 8.538486e-240 [94,] 1.0000000 1.763743e-238 8.818716e-239 [95,] 1.0000000 5.705230e-239 2.852615e-239 [96,] 1.0000000 4.298996e-238 2.149498e-238 [97,] 1.0000000 8.962346e-238 4.481173e-238 [98,] 1.0000000 4.929217e-237 2.464609e-237 [99,] 1.0000000 1.128768e-236 5.643840e-237 [100,] 1.0000000 2.713777e-237 1.356888e-237 [101,] 1.0000000 1.474139e-237 7.370694e-238 [102,] 1.0000000 8.989644e-237 4.494822e-237 [103,] 1.0000000 3.069808e-236 1.534904e-236 [104,] 1.0000000 2.801644e-236 1.400822e-236 [105,] 1.0000000 1.431704e-236 7.158518e-237 [106,] 1.0000000 3.081348e-238 1.540674e-238 [107,] 1.0000000 1.116047e-237 5.580237e-238 [108,] 1.0000000 4.254427e-237 2.127214e-237 [109,] 1.0000000 1.048865e-236 5.244324e-237 [110,] 1.0000000 8.603885e-236 4.301943e-236 [111,] 1.0000000 1.826082e-235 9.130412e-236 [112,] 1.0000000 5.788313e-237 2.894156e-237 [113,] 1.0000000 4.967795e-236 2.483897e-236 [114,] 1.0000000 6.263342e-236 3.131671e-236 [115,] 1.0000000 6.307007e-235 3.153503e-235 [116,] 1.0000000 5.573238e-234 2.786619e-234 [117,] 1.0000000 4.707670e-233 2.353835e-233 [118,] 1.0000000 5.924328e-232 2.962164e-232 [119,] 1.0000000 1.345509e-232 6.727543e-233 [120,] 1.0000000 1.487906e-231 7.439530e-232 [121,] 1.0000000 1.656035e-230 8.280173e-231 [122,] 1.0000000 1.720653e-229 8.603263e-230 [123,] 1.0000000 1.931005e-228 9.655027e-229 [124,] 1.0000000 2.164572e-227 1.082286e-227 [125,] 1.0000000 1.002074e-227 5.010368e-228 [126,] 1.0000000 1.149275e-226 5.746376e-227 [127,] 1.0000000 1.327898e-225 6.639489e-226 [128,] 1.0000000 1.529012e-224 7.645061e-225 [129,] 1.0000000 1.658907e-223 8.294537e-224 [130,] 1.0000000 1.908379e-222 9.541896e-223 [131,] 1.0000000 2.204916e-221 1.102458e-221 [132,] 1.0000000 2.548805e-220 1.274403e-220 [133,] 1.0000000 2.946973e-219 1.473487e-219 [134,] 1.0000000 3.299717e-218 1.649858e-218 [135,] 1.0000000 3.813818e-217 1.906909e-217 [136,] 1.0000000 3.287214e-216 1.643607e-216 [137,] 1.0000000 3.795535e-215 1.897768e-215 [138,] 1.0000000 3.752052e-214 1.876026e-214 [139,] 1.0000000 3.672779e-213 1.836390e-213 [140,] 1.0000000 4.233931e-212 2.116966e-212 [141,] 1.0000000 4.865316e-211 2.432658e-211 [142,] 1.0000000 5.591674e-210 2.795837e-210 [143,] 1.0000000 6.424036e-209 3.212018e-209 [144,] 1.0000000 7.251424e-208 3.625712e-208 [145,] 1.0000000 5.802227e-207 2.901113e-207 [146,] 1.0000000 6.621118e-206 3.310559e-206 [147,] 1.0000000 7.624161e-205 3.812080e-205 [148,] 1.0000000 8.662204e-204 4.331102e-204 [149,] 1.0000000 9.813242e-203 4.906621e-203 [150,] 1.0000000 1.108966e-201 5.544831e-202 [151,] 1.0000000 1.254566e-200 6.272828e-201 [152,] 1.0000000 1.185018e-199 5.925088e-200 [153,] 1.0000000 1.329196e-198 6.645981e-199 [154,] 1.0000000 1.486727e-197 7.433635e-198 [155,] 1.0000000 1.658135e-196 8.290677e-197 [156,] 1.0000000 1.843851e-195 9.219256e-196 [157,] 1.0000000 1.858139e-194 9.290694e-195 [158,] 1.0000000 2.058735e-193 1.029368e-193 [159,] 1.0000000 2.268535e-192 1.134267e-192 [160,] 1.0000000 2.491774e-191 1.245887e-191 [161,] 1.0000000 1.405656e-192 7.028278e-193 [162,] 1.0000000 1.566379e-191 7.831895e-192 [163,] 1.0000000 1.739600e-190 8.697999e-191 [164,] 1.0000000 1.925385e-189 9.626923e-190 [165,] 1.0000000 2.135505e-188 1.067753e-188 [166,] 1.0000000 2.347107e-187 1.173553e-187 [167,] 1.0000000 2.557155e-186 1.278578e-186 [168,] 1.0000000 2.790494e-185 1.395247e-185 [169,] 1.0000000 3.030016e-184 1.515008e-184 [170,] 1.0000000 3.282594e-183 1.641297e-183 [171,] 1.0000000 3.679499e-182 1.839749e-182 [172,] 1.0000000 3.952615e-181 1.976307e-181 [173,] 1.0000000 4.230357e-180 2.115178e-180 [174,] 1.0000000 2.357407e-179 1.178704e-179 [175,] 1.0000000 2.586125e-178 1.293063e-178 [176,] 1.0000000 1.685892e-177 8.429458e-178 [177,] 1.0000000 8.074023e-177 4.037011e-177 [178,] 1.0000000 8.867282e-176 4.433641e-176 [179,] 1.0000000 9.083079e-175 4.541540e-175 [180,] 1.0000000 7.528753e-174 3.764377e-174 [181,] 1.0000000 7.854207e-173 3.927104e-173 [182,] 1.0000000 6.544078e-172 3.272039e-172 [183,] 1.0000000 6.780985e-171 3.390493e-171 [184,] 1.0000000 6.387348e-170 3.193674e-170 [185,] 1.0000000 6.569784e-169 3.284892e-169 [186,] 1.0000000 6.576009e-168 3.288005e-168 [187,] 1.0000000 6.711888e-167 3.355944e-167 [188,] 1.0000000 6.822949e-166 3.411474e-166 [189,] 1.0000000 7.130121e-165 3.565060e-165 [190,] 1.0000000 2.617246e-164 1.308623e-164 [191,] 1.0000000 2.725197e-163 1.362598e-163 [192,] 1.0000000 2.824335e-162 1.412168e-162 [193,] 1.0000000 7.898688e-162 3.949344e-162 [194,] 1.0000000 8.164706e-161 4.082353e-161 [195,] 1.0000000 8.400210e-160 4.200105e-160 [196,] 1.0000000 8.602070e-159 4.301035e-159 [197,] 1.0000000 7.120288e-158 3.560144e-158 [198,] 1.0000000 1.373262e-157 6.866309e-158 [199,] 1.0000000 1.394830e-156 6.974148e-157 [200,] 1.0000000 1.312448e-155 6.562240e-156 [201,] 1.0000000 1.324345e-154 6.621726e-155 [202,] 1.0000000 9.818584e-155 4.909292e-155 [203,] 1.0000000 9.937227e-154 4.968614e-154 [204,] 1.0000000 9.680101e-153 4.840051e-153 [205,] 1.0000000 9.392578e-152 4.696289e-152 [206,] 1.0000000 8.083842e-151 4.041921e-151 [207,] 1.0000000 7.816793e-150 3.908396e-150 [208,] 1.0000000 7.576390e-149 3.788195e-149 [209,] 1.0000000 7.237398e-148 3.618699e-148 [210,] 1.0000000 6.848968e-147 3.424484e-147 [211,] 1.0000000 6.139933e-146 3.069966e-146 [212,] 1.0000000 5.795550e-145 2.897775e-145 [213,] 1.0000000 3.040793e-144 1.520396e-144 [214,] 1.0000000 2.647716e-143 1.323858e-143 [215,] 1.0000000 2.475822e-142 1.237911e-142 [216,] 1.0000000 2.302548e-141 1.151274e-141 [217,] 1.0000000 2.135280e-140 1.067640e-140 [218,] 1.0000000 1.969661e-139 9.848305e-140 [219,] 1.0000000 1.756347e-138 8.781733e-139 [220,] 1.0000000 1.608571e-137 8.042855e-138 [221,] 1.0000000 1.467017e-136 7.335087e-137 [222,] 1.0000000 1.261361e-135 6.306807e-136 [223,] 1.0000000 1.141241e-134 5.706205e-135 [224,] 1.0000000 1.026456e-133 5.132281e-134 [225,] 1.0000000 9.209166e-133 4.604583e-133 [226,] 1.0000000 8.226863e-132 4.113432e-132 [227,] 1.0000000 6.951779e-131 3.475890e-131 [228,] 1.0000000 6.157668e-130 3.078834e-130 [229,] 1.0000000 5.430229e-129 2.715114e-129 [230,] 1.0000000 4.768267e-128 2.384134e-128 [231,] 1.0000000 4.122384e-127 2.061192e-127 [232,] 1.0000000 3.588340e-126 1.794170e-126 [233,] 1.0000000 3.109754e-125 1.554877e-125 [234,] 1.0000000 2.683114e-124 1.341557e-124 [235,] 1.0000000 2.253351e-123 1.126675e-123 [236,] 1.0000000 1.927053e-122 9.635264e-123 [237,] 1.0000000 1.617826e-121 8.089131e-122 [238,] 1.0000000 1.376229e-120 6.881147e-121 [239,] 1.0000000 1.184151e-119 5.920753e-120 [240,] 1.0000000 1.014023e-118 5.070116e-119 [241,] 1.0000000 8.545465e-118 4.272732e-118 [242,] 1.0000000 7.248002e-117 3.624001e-117 [243,] 1.0000000 6.118058e-116 3.059029e-116 [244,] 1.0000000 4.654220e-115 2.327110e-115 [245,] 1.0000000 3.892974e-114 1.946487e-114 [246,] 1.0000000 3.098996e-113 1.549498e-113 [247,] 1.0000000 2.357610e-112 1.178805e-112 [248,] 1.0000000 1.130051e-111 5.650257e-112 [249,] 1.0000000 9.135237e-111 4.567619e-111 [250,] 1.0000000 4.865065e-110 2.432532e-110 [251,] 1.0000000 2.180001e-109 1.090001e-109 [252,] 1.0000000 1.780029e-108 8.900146e-109 [253,] 1.0000000 1.418969e-107 7.094844e-108 [254,] 1.0000000 1.147383e-106 5.736913e-107 [255,] 1.0000000 9.232096e-106 4.616048e-106 [256,] 1.0000000 7.093659e-105 3.546830e-105 [257,] 1.0000000 5.022811e-104 2.511405e-104 [258,] 1.0000000 3.985665e-103 1.992832e-103 [259,] 1.0000000 3.089676e-102 1.544838e-102 [260,] 1.0000000 2.427561e-101 1.213781e-101 [261,] 1.0000000 1.849520e-100 9.247598e-101 [262,] 1.0000000 1.438875e-99 7.194375e-100 [263,] 1.0000000 1.058158e-98 5.290790e-99 [264,] 1.0000000 8.152281e-98 4.076141e-98 [265,] 1.0000000 6.248897e-97 3.124449e-97 [266,] 1.0000000 4.765572e-96 2.382786e-96 [267,] 1.0000000 3.566662e-95 1.783331e-95 [268,] 1.0000000 2.647390e-94 1.323695e-94 [269,] 1.0000000 1.882924e-93 9.414622e-94 [270,] 1.0000000 1.345300e-92 6.726500e-93 [271,] 1.0000000 9.842506e-92 4.921253e-92 [272,] 1.0000000 7.072962e-91 3.536481e-91 [273,] 1.0000000 4.739310e-90 2.369655e-90 [274,] 1.0000000 3.475645e-89 1.737823e-89 [275,] 1.0000000 2.494409e-88 1.247205e-88 [276,] 1.0000000 1.523831e-87 7.619155e-88 [277,] 1.0000000 1.083975e-86 5.419873e-87 [278,] 1.0000000 6.945388e-86 3.472694e-86 [279,] 1.0000000 4.896810e-85 2.448405e-85 [280,] 1.0000000 3.488369e-84 1.744184e-84 [281,] 1.0000000 2.443689e-83 1.221844e-83 [282,] 1.0000000 1.722462e-82 8.612309e-83 [283,] 1.0000000 1.049381e-81 5.246904e-82 [284,] 1.0000000 7.326215e-81 3.663107e-81 [285,] 1.0000000 4.867411e-80 2.433705e-80 [286,] 1.0000000 3.201639e-79 1.600820e-79 [287,] 1.0000000 2.168561e-78 1.084281e-78 [288,] 1.0000000 1.456044e-77 7.280218e-78 [289,] 1.0000000 6.839108e-77 3.419554e-77 [290,] 1.0000000 4.571728e-76 2.285864e-76 [291,] 1.0000000 3.079540e-75 1.539770e-75 [292,] 1.0000000 2.036328e-74 1.018164e-74 [293,] 1.0000000 1.339189e-73 6.695944e-74 [294,] 1.0000000 8.879890e-73 4.439945e-73 [295,] 1.0000000 5.821210e-72 2.910605e-72 [296,] 1.0000000 3.765404e-71 1.882702e-71 [297,] 1.0000000 2.452459e-70 1.226230e-70 [298,] 1.0000000 1.162454e-69 5.812268e-70 [299,] 1.0000000 7.382815e-69 3.691408e-69 [300,] 1.0000000 4.741520e-68 2.370760e-68 [301,] 1.0000000 3.028420e-67 1.514210e-67 [302,] 1.0000000 1.475971e-66 7.379853e-67 [303,] 1.0000000 9.338041e-66 4.669021e-66 [304,] 1.0000000 5.872146e-65 2.936073e-65 [305,] 1.0000000 2.936485e-64 1.468243e-64 [306,] 1.0000000 1.825881e-63 9.129403e-64 [307,] 1.0000000 1.131785e-62 5.658926e-63 [308,] 1.0000000 6.863178e-62 3.431589e-62 [309,] 1.0000000 4.151007e-61 2.075504e-61 [310,] 1.0000000 2.470891e-60 1.235446e-60 [311,] 1.0000000 1.491235e-59 7.456173e-60 [312,] 1.0000000 8.828956e-59 4.414478e-59 [313,] 1.0000000 5.204148e-58 2.602074e-58 [314,] 1.0000000 3.048734e-57 1.524367e-57 [315,] 1.0000000 1.775007e-56 8.875034e-57 [316,] 1.0000000 1.027004e-55 5.135020e-56 [317,] 1.0000000 5.971449e-55 2.985725e-55 [318,] 1.0000000 3.411872e-54 1.705936e-54 [319,] 1.0000000 1.948093e-53 9.740467e-54 [320,] 1.0000000 1.110510e-52 5.552548e-53 [321,] 1.0000000 6.283688e-52 3.141844e-52 [322,] 1.0000000 3.499337e-51 1.749669e-51 [323,] 1.0000000 1.936023e-50 9.680116e-51 [324,] 1.0000000 1.064060e-49 5.320300e-50 [325,] 1.0000000 5.809375e-49 2.904688e-49 [326,] 1.0000000 3.150480e-48 1.575240e-48 [327,] 1.0000000 1.687921e-47 8.439605e-48 [328,] 1.0000000 9.027065e-47 4.513532e-47 [329,] 1.0000000 4.796323e-46 2.398162e-46 [330,] 1.0000000 2.536588e-45 1.268294e-45 [331,] 1.0000000 1.329130e-44 6.645648e-45 [332,] 1.0000000 6.917399e-44 3.458699e-44 [333,] 1.0000000 3.529195e-43 1.764597e-43 [334,] 1.0000000 1.834838e-42 9.174190e-43 [335,] 1.0000000 9.346627e-42 4.673314e-42 [336,] 1.0000000 4.726385e-41 2.363192e-41 [337,] 1.0000000 2.372407e-40 1.186204e-40 [338,] 1.0000000 1.181962e-39 5.909810e-40 [339,] 1.0000000 5.844391e-39 2.922195e-39 [340,] 1.0000000 2.899921e-38 1.449961e-38 [341,] 1.0000000 1.412231e-37 7.061154e-38 [342,] 1.0000000 6.824082e-37 3.412041e-37 [343,] 1.0000000 3.271637e-36 1.635819e-36 [344,] 1.0000000 1.513677e-35 7.568386e-36 [345,] 1.0000000 7.155764e-35 3.577882e-35 [346,] 1.0000000 3.037773e-34 1.518887e-34 [347,] 1.0000000 1.429665e-33 7.148323e-34 [348,] 1.0000000 6.586527e-33 3.293264e-33 [349,] 1.0000000 3.008965e-32 1.504482e-32 [350,] 1.0000000 1.274788e-31 6.373938e-32 [351,] 1.0000000 5.737952e-31 2.868976e-31 [352,] 1.0000000 2.560267e-30 1.280133e-30 [353,] 1.0000000 1.132330e-29 5.661650e-30 [354,] 1.0000000 4.963472e-29 2.481736e-29 [355,] 1.0000000 2.156024e-28 1.078012e-28 [356,] 1.0000000 9.279436e-28 4.639718e-28 [357,] 1.0000000 4.043815e-27 2.021908e-27 [358,] 1.0000000 1.704953e-26 8.524763e-27 [359,] 1.0000000 7.137058e-26 3.568529e-26 [360,] 1.0000000 2.958685e-25 1.479343e-25 [361,] 1.0000000 1.214464e-24 6.072322e-25 [362,] 1.0000000 4.935236e-24 2.467618e-24 [363,] 1.0000000 1.985162e-23 9.925808e-24 [364,] 1.0000000 7.902651e-23 3.951326e-23 [365,] 1.0000000 3.103054e-22 1.551527e-22 [366,] 1.0000000 1.190104e-21 5.950520e-22 [367,] 1.0000000 4.714785e-21 2.357393e-21 [368,] 1.0000000 1.806563e-20 9.032814e-21 [369,] 1.0000000 6.844780e-20 3.422390e-20 [370,] 1.0000000 2.506761e-19 1.253380e-19 [371,] 1.0000000 9.182415e-19 4.591207e-19 [372,] 1.0000000 3.360725e-18 1.680363e-18 [373,] 1.0000000 1.215114e-17 6.075571e-18 [374,] 1.0000000 4.338996e-17 2.169498e-17 [375,] 1.0000000 1.578238e-16 7.891192e-17 [376,] 1.0000000 4.187702e-16 2.093851e-16 [377,] 1.0000000 1.449459e-15 7.247296e-16 [378,] 1.0000000 4.981144e-15 2.490572e-15 [379,] 1.0000000 1.456507e-14 7.282534e-15 [380,] 1.0000000 4.766819e-14 2.383410e-14 [381,] 1.0000000 1.469654e-13 7.348271e-14 [382,] 1.0000000 4.790807e-13 2.395404e-13 [383,] 1.0000000 1.550189e-12 7.750946e-13 [384,] 1.0000000 4.786014e-12 2.393007e-12 [385,] 1.0000000 1.512044e-11 7.560220e-12 [386,] 1.0000000 4.800763e-11 2.400382e-11 [387,] 1.0000000 1.208213e-10 6.041065e-11 [388,] 1.0000000 2.583391e-10 1.291696e-10 [389,] 1.0000000 7.845136e-10 3.922568e-10 [390,] 1.0000000 2.118420e-09 1.059210e-09 [391,] 1.0000000 6.322213e-09 3.161106e-09 [392,] 1.0000000 1.832641e-08 9.163205e-09 [393,] 1.0000000 4.564867e-08 2.282434e-08 [394,] 0.9999999 1.274466e-07 6.372329e-08 [395,] 0.9999998 3.248774e-07 1.624387e-07 [396,] 0.9999996 8.473296e-07 4.236648e-07 [397,] 0.9999989 2.212691e-06 1.106345e-06 [398,] 0.9999971 5.793021e-06 2.896510e-06 [399,] 0.9999926 1.487337e-05 7.436683e-06 [400,] 0.9999821 3.578491e-05 1.789245e-05 [401,] 0.9999606 7.886807e-05 3.943403e-05 [402,] 0.9999093 1.814322e-04 9.071609e-05 [403,] 0.9997811 4.378506e-04 2.189253e-04 [404,] 0.9995127 9.746259e-04 4.873129e-04 [405,] 0.9989106 2.178836e-03 1.089418e-03 [406,] 0.9976088 4.782364e-03 2.391182e-03 [407,] 0.9949801 1.003978e-02 5.019890e-03 [408,] 0.9896587 2.068257e-02 1.034129e-02 [409,] 0.9794636 4.107281e-02 2.053640e-02 [410,] 0.9925859 1.482821e-02 7.414103e-03 [411,] 0.9843129 3.137426e-02 1.568713e-02 [412,] 0.9701958 5.960839e-02 2.980419e-02 [413,] 0.9392267 1.215467e-01 6.077333e-02 [414,] 0.8791527 2.416947e-01 1.208473e-01 [415,] 0.7984940 4.030121e-01 2.015060e-01 [416,] 0.6489841 7.020319e-01 3.510159e-01 > postscript(file="/var/www/html/rcomp/tmp/19sbw1291302152.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/29sbw1291302152.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/3kjsz1291302152.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/4kjsz1291302152.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/5kjsz1291302152.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 1614392.9083 3267339.8005 1342318.8812 -3145447.1126 -392660.3716 6 7 8 9 10 47198.0501 404.7394 602188.5762 -211744.7650 769347.7162 11 12 13 14 15 -324970.7246 476487.5596 971510.2238 -1780770.7218 52895.3919 16 17 18 19 20 440956.3216 -411358.2358 -67512.6769 -148150.4768 -328203.7754 21 22 23 24 25 173412.6685 -246220.5692 -237869.2592 -488859.5214 -214929.7252 26 27 28 29 30 -127584.6004 -365900.8069 -441606.0856 404846.4220 24692.1759 31 32 33 34 35 180623.9925 -162947.7947 132677.0574 16731.1667 -338033.7242 36 37 38 39 40 -454918.3637 55071.1870 73640.4234 -122998.4009 -149060.1825 41 42 43 44 45 23536.2045 -354630.1910 103738.8678 -186161.7635 -11873.8924 46 47 48 49 50 -232458.0978 -248121.3384 -171890.4711 -103951.1514 -156677.2080 51 52 53 54 55 25768.3763 -125901.9966 -405539.2242 -48526.5126 -394622.0142 56 57 58 59 60 134777.1577 182089.8545 -18512.2564 -20266.4483 -226136.2158 61 62 63 64 65 -89112.9730 -57073.6173 -109043.6851 -232356.0789 65003.2838 66 67 68 69 70 81678.9589 -91906.2811 -99661.3241 -58130.9141 -136377.9104 71 72 73 74 75 -28840.9276 73016.2831 -125491.3547 -14991.2288 -121785.5624 76 77 78 79 80 -141157.5865 -92252.6956 -47277.0402 87516.0975 -56250.2459 81 82 83 84 85 51772.1918 -98359.6297 -118980.2572 89337.8404 -5467.0157 86 87 88 89 90 -111554.4237 30263.8277 -39334.6250 -95805.4969 -20137.4147 91 92 93 94 95 -26562.2862 -102847.2821 136445.8672 -339389.1301 142088.0632 96 97 98 99 100 -22029.2939 12387.1874 47889.5549 -132866.7012 -85702.5943 101 102 103 104 105 -93984.4068 -101209.0829 -110607.5576 -53863.5304 59661.3416 106 107 108 109 110 -44169.0784 -239180.3303 -1079.6495 -203731.6259 -152145.4789 111 112 113 114 115 116244.2019 -234855.3016 -268536.7701 -174382.2503 -7876.3314 116 117 118 119 120 79594.8195 11626.7819 40981.2854 -237174.4564 77984.4846 121 122 123 124 125 -260529.3804 44126.2494 2437.0379 -48248.7533 -120768.4011 126 127 128 129 130 -249849.1739 56743.8463 57048.3791 64972.9222 56743.8463 131 132 133 134 135 59962.9990 114381.9051 56743.8463 54723.3658 56743.8463 136 137 138 139 140 59317.8032 59103.1536 56743.8463 56743.8463 56743.8463 141 142 143 144 145 55786.6489 56743.8463 51497.1868 57019.2654 63025.8180 146 147 148 149 150 68710.2635 56743.8463 58723.7837 56743.8463 56178.8235 151 152 153 154 155 47378.3959 46154.5855 56743.8463 34483.4455 56844.1304 156 157 158 159 160 56743.8463 56743.8463 53573.6402 63896.5616 56743.8463 161 162 163 164 165 56743.8463 56743.8463 56743.8463 62546.0270 57582.2679 166 167 168 169 170 56743.8463 56743.8463 160360.4989 56743.8463 56743.8463 171 172 173 174 175 56743.8463 53760.2551 56743.8463 54511.6315 56743.8463 176 177 178 179 180 56801.5725 56743.8463 -15353.1898 56743.8463 56743.8463 181 182 183 184 185 70989.7457 -37474.7381 82509.7404 -76005.3479 -14797.0131 186 187 188 189 190 -21688.9005 -8918.1814 53423.5954 17043.9417 56743.8463 191 192 193 194 195 61395.0315 56743.8463 58916.3004 56743.8463 56474.3785 196 197 198 199 200 -14102.2631 -150102.8416 -14102.2631 -14102.2631 -92005.6434 201 202 203 204 205 -14102.2631 -14102.2631 -14102.2631 14151.1709 -100326.1605 206 207 208 209 210 -15827.4682 -29262.9780 -14102.2631 44819.1609 -14102.2631 211 212 213 214 215 56743.8463 56743.8463 72687.1102 34552.6110 37988.8278 216 217 218 219 220 56743.8463 60320.3399 61069.6158 56743.8463 61629.1210 221 222 223 224 225 23935.1655 56743.8463 57658.7888 56743.8463 57060.4504 226 227 228 229 230 61626.0995 56743.8463 56743.8463 30572.5273 56743.8463 231 232 233 234 235 35628.8129 56743.8463 56743.8463 58191.0580 56743.8463 236 237 238 239 240 57007.1910 56743.8463 59739.9250 56743.8463 56743.8463 241 242 243 244 245 56743.8463 59635.8630 56743.8463 61991.0854 -30514.6252 246 247 248 249 250 -14102.2631 -14102.2631 48670.3819 -14102.2631 -14102.2631 251 252 253 254 255 63819.4336 -14102.2631 33344.1466 37600.1858 6243.3979 256 257 258 259 260 56743.8463 80097.7803 96741.5415 -14102.2631 56743.8463 261 262 263 264 265 -14102.2631 -14102.2631 37163.9577 65359.4770 -14102.2631 266 267 268 269 270 56743.8463 -14102.2631 60145.7968 -14102.2631 35248.8380 271 272 273 274 275 -14102.2631 -14102.2631 -14102.2631 55325.6036 56743.8463 276 277 278 279 280 66053.6292 63962.5695 56743.8463 -25065.7169 35917.4441 281 282 283 284 285 -14102.2631 56743.8463 95798.9808 56743.8463 20948.2632 286 287 288 289 290 56743.8463 -14102.2631 55796.2699 -14102.2631 87737.0651 291 292 293 294 295 -14102.2631 66657.4778 30150.0010 56743.8463 60459.4379 296 297 298 299 300 20666.8015 56743.8463 -14102.2631 56743.8463 56743.8463 301 302 303 304 305 -14799.5421 42657.9409 56743.8463 -14102.2631 96052.7763 306 307 308 309 310 60905.0179 -14102.2631 -22026.4417 78605.8994 -14102.2631 311 312 313 314 315 -14102.2631 37401.7567 -14102.2631 -17024.1271 58714.9646 316 317 318 319 320 57989.7061 58333.6778 42709.8089 59542.2511 56743.8463 321 322 323 324 325 56743.8463 56743.8463 56743.8463 -14102.2631 56743.8463 326 327 328 329 330 54184.7061 -14102.2631 -14102.2631 56743.8463 56743.8463 331 332 333 334 335 56743.8463 56743.8463 56743.8463 56785.5558 57641.3511 336 337 338 339 340 56743.8463 55622.0937 56743.8463 56098.6316 60825.7496 341 342 343 344 345 37739.2850 56743.8463 56743.8463 56743.8463 56743.8463 346 347 348 349 350 56743.8463 48733.4452 56743.8463 56743.8463 56743.8463 351 352 353 354 355 57507.4643 58135.5877 58182.2306 46857.2332 56743.8463 356 357 358 359 360 56743.8463 24153.8623 56743.8463 56743.8463 56743.8463 361 362 363 364 365 58592.2977 56743.8463 56743.8463 -23229.0752 55339.9264 366 367 368 369 370 56743.8463 56743.8463 56743.8463 56743.8463 56743.8463 371 372 373 374 375 56743.8463 59951.3846 16950.6168 29944.4766 56743.8463 376 377 378 379 380 56743.8463 69183.7439 41529.1872 56743.8463 56743.8463 381 382 383 384 385 56743.8463 50304.4797 110203.6465 56743.8463 58971.6081 386 387 388 389 390 -99802.4349 61423.9876 40725.4677 54932.2314 37401.4713 391 392 393 394 395 53556.4003 36181.4677 18933.7581 73314.7260 -225558.1571 396 397 398 399 400 29067.6811 -107954.0035 -1003.1517 30167.9153 -20626.9915 401 402 403 404 405 18801.9662 46542.0025 35992.8637 -96274.3599 23096.4155 406 407 408 409 410 4588.9031 41079.1019 36064.8764 34340.5635 -18627.8954 411 412 413 414 415 12187.3687 -219006.8707 -43752.6695 -154815.6019 -54634.8218 416 417 418 419 420 -37896.3342 138483.9631 -67743.4506 5703.9594 -11870.2290 421 422 423 424 425 -65172.7527 -99395.9490 -427200.0501 -470445.0446 -71049.1911 426 427 428 429 430 -200544.6914 -306223.7944 -1206200.7325 -359249.8705 -270057.1479 431 -843491.5827 > postscript(file="/var/www/html/rcomp/tmp/6dt9k1291302152.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 1614392.9083 NA 1 3267339.8005 1614392.9083 2 1342318.8812 3267339.8005 3 -3145447.1126 1342318.8812 4 -392660.3716 -3145447.1126 5 47198.0501 -392660.3716 6 404.7394 47198.0501 7 602188.5762 404.7394 8 -211744.7650 602188.5762 9 769347.7162 -211744.7650 10 -324970.7246 769347.7162 11 476487.5596 -324970.7246 12 971510.2238 476487.5596 13 -1780770.7218 971510.2238 14 52895.3919 -1780770.7218 15 440956.3216 52895.3919 16 -411358.2358 440956.3216 17 -67512.6769 -411358.2358 18 -148150.4768 -67512.6769 19 -328203.7754 -148150.4768 20 173412.6685 -328203.7754 21 -246220.5692 173412.6685 22 -237869.2592 -246220.5692 23 -488859.5214 -237869.2592 24 -214929.7252 -488859.5214 25 -127584.6004 -214929.7252 26 -365900.8069 -127584.6004 27 -441606.0856 -365900.8069 28 404846.4220 -441606.0856 29 24692.1759 404846.4220 30 180623.9925 24692.1759 31 -162947.7947 180623.9925 32 132677.0574 -162947.7947 33 16731.1667 132677.0574 34 -338033.7242 16731.1667 35 -454918.3637 -338033.7242 36 55071.1870 -454918.3637 37 73640.4234 55071.1870 38 -122998.4009 73640.4234 39 -149060.1825 -122998.4009 40 23536.2045 -149060.1825 41 -354630.1910 23536.2045 42 103738.8678 -354630.1910 43 -186161.7635 103738.8678 44 -11873.8924 -186161.7635 45 -232458.0978 -11873.8924 46 -248121.3384 -232458.0978 47 -171890.4711 -248121.3384 48 -103951.1514 -171890.4711 49 -156677.2080 -103951.1514 50 25768.3763 -156677.2080 51 -125901.9966 25768.3763 52 -405539.2242 -125901.9966 53 -48526.5126 -405539.2242 54 -394622.0142 -48526.5126 55 134777.1577 -394622.0142 56 182089.8545 134777.1577 57 -18512.2564 182089.8545 58 -20266.4483 -18512.2564 59 -226136.2158 -20266.4483 60 -89112.9730 -226136.2158 61 -57073.6173 -89112.9730 62 -109043.6851 -57073.6173 63 -232356.0789 -109043.6851 64 65003.2838 -232356.0789 65 81678.9589 65003.2838 66 -91906.2811 81678.9589 67 -99661.3241 -91906.2811 68 -58130.9141 -99661.3241 69 -136377.9104 -58130.9141 70 -28840.9276 -136377.9104 71 73016.2831 -28840.9276 72 -125491.3547 73016.2831 73 -14991.2288 -125491.3547 74 -121785.5624 -14991.2288 75 -141157.5865 -121785.5624 76 -92252.6956 -141157.5865 77 -47277.0402 -92252.6956 78 87516.0975 -47277.0402 79 -56250.2459 87516.0975 80 51772.1918 -56250.2459 81 -98359.6297 51772.1918 82 -118980.2572 -98359.6297 83 89337.8404 -118980.2572 84 -5467.0157 89337.8404 85 -111554.4237 -5467.0157 86 30263.8277 -111554.4237 87 -39334.6250 30263.8277 88 -95805.4969 -39334.6250 89 -20137.4147 -95805.4969 90 -26562.2862 -20137.4147 91 -102847.2821 -26562.2862 92 136445.8672 -102847.2821 93 -339389.1301 136445.8672 94 142088.0632 -339389.1301 95 -22029.2939 142088.0632 96 12387.1874 -22029.2939 97 47889.5549 12387.1874 98 -132866.7012 47889.5549 99 -85702.5943 -132866.7012 100 -93984.4068 -85702.5943 101 -101209.0829 -93984.4068 102 -110607.5576 -101209.0829 103 -53863.5304 -110607.5576 104 59661.3416 -53863.5304 105 -44169.0784 59661.3416 106 -239180.3303 -44169.0784 107 -1079.6495 -239180.3303 108 -203731.6259 -1079.6495 109 -152145.4789 -203731.6259 110 116244.2019 -152145.4789 111 -234855.3016 116244.2019 112 -268536.7701 -234855.3016 113 -174382.2503 -268536.7701 114 -7876.3314 -174382.2503 115 79594.8195 -7876.3314 116 11626.7819 79594.8195 117 40981.2854 11626.7819 118 -237174.4564 40981.2854 119 77984.4846 -237174.4564 120 -260529.3804 77984.4846 121 44126.2494 -260529.3804 122 2437.0379 44126.2494 123 -48248.7533 2437.0379 124 -120768.4011 -48248.7533 125 -249849.1739 -120768.4011 126 56743.8463 -249849.1739 127 57048.3791 56743.8463 128 64972.9222 57048.3791 129 56743.8463 64972.9222 130 59962.9990 56743.8463 131 114381.9051 59962.9990 132 56743.8463 114381.9051 133 54723.3658 56743.8463 134 56743.8463 54723.3658 135 59317.8032 56743.8463 136 59103.1536 59317.8032 137 56743.8463 59103.1536 138 56743.8463 56743.8463 139 56743.8463 56743.8463 140 55786.6489 56743.8463 141 56743.8463 55786.6489 142 51497.1868 56743.8463 143 57019.2654 51497.1868 144 63025.8180 57019.2654 145 68710.2635 63025.8180 146 56743.8463 68710.2635 147 58723.7837 56743.8463 148 56743.8463 58723.7837 149 56178.8235 56743.8463 150 47378.3959 56178.8235 151 46154.5855 47378.3959 152 56743.8463 46154.5855 153 34483.4455 56743.8463 154 56844.1304 34483.4455 155 56743.8463 56844.1304 156 56743.8463 56743.8463 157 53573.6402 56743.8463 158 63896.5616 53573.6402 159 56743.8463 63896.5616 160 56743.8463 56743.8463 161 56743.8463 56743.8463 162 56743.8463 56743.8463 163 62546.0270 56743.8463 164 57582.2679 62546.0270 165 56743.8463 57582.2679 166 56743.8463 56743.8463 167 160360.4989 56743.8463 168 56743.8463 160360.4989 169 56743.8463 56743.8463 170 56743.8463 56743.8463 171 53760.2551 56743.8463 172 56743.8463 53760.2551 173 54511.6315 56743.8463 174 56743.8463 54511.6315 175 56801.5725 56743.8463 176 56743.8463 56801.5725 177 -15353.1898 56743.8463 178 56743.8463 -15353.1898 179 56743.8463 56743.8463 180 70989.7457 56743.8463 181 -37474.7381 70989.7457 182 82509.7404 -37474.7381 183 -76005.3479 82509.7404 184 -14797.0131 -76005.3479 185 -21688.9005 -14797.0131 186 -8918.1814 -21688.9005 187 53423.5954 -8918.1814 188 17043.9417 53423.5954 189 56743.8463 17043.9417 190 61395.0315 56743.8463 191 56743.8463 61395.0315 192 58916.3004 56743.8463 193 56743.8463 58916.3004 194 56474.3785 56743.8463 195 -14102.2631 56474.3785 196 -150102.8416 -14102.2631 197 -14102.2631 -150102.8416 198 -14102.2631 -14102.2631 199 -92005.6434 -14102.2631 200 -14102.2631 -92005.6434 201 -14102.2631 -14102.2631 202 -14102.2631 -14102.2631 203 14151.1709 -14102.2631 204 -100326.1605 14151.1709 205 -15827.4682 -100326.1605 206 -29262.9780 -15827.4682 207 -14102.2631 -29262.9780 208 44819.1609 -14102.2631 209 -14102.2631 44819.1609 210 56743.8463 -14102.2631 211 56743.8463 56743.8463 212 72687.1102 56743.8463 213 34552.6110 72687.1102 214 37988.8278 34552.6110 215 56743.8463 37988.8278 216 60320.3399 56743.8463 217 61069.6158 60320.3399 218 56743.8463 61069.6158 219 61629.1210 56743.8463 220 23935.1655 61629.1210 221 56743.8463 23935.1655 222 57658.7888 56743.8463 223 56743.8463 57658.7888 224 57060.4504 56743.8463 225 61626.0995 57060.4504 226 56743.8463 61626.0995 227 56743.8463 56743.8463 228 30572.5273 56743.8463 229 56743.8463 30572.5273 230 35628.8129 56743.8463 231 56743.8463 35628.8129 232 56743.8463 56743.8463 233 58191.0580 56743.8463 234 56743.8463 58191.0580 235 57007.1910 56743.8463 236 56743.8463 57007.1910 237 59739.9250 56743.8463 238 56743.8463 59739.9250 239 56743.8463 56743.8463 240 56743.8463 56743.8463 241 59635.8630 56743.8463 242 56743.8463 59635.8630 243 61991.0854 56743.8463 244 -30514.6252 61991.0854 245 -14102.2631 -30514.6252 246 -14102.2631 -14102.2631 247 48670.3819 -14102.2631 248 -14102.2631 48670.3819 249 -14102.2631 -14102.2631 250 63819.4336 -14102.2631 251 -14102.2631 63819.4336 252 33344.1466 -14102.2631 253 37600.1858 33344.1466 254 6243.3979 37600.1858 255 56743.8463 6243.3979 256 80097.7803 56743.8463 257 96741.5415 80097.7803 258 -14102.2631 96741.5415 259 56743.8463 -14102.2631 260 -14102.2631 56743.8463 261 -14102.2631 -14102.2631 262 37163.9577 -14102.2631 263 65359.4770 37163.9577 264 -14102.2631 65359.4770 265 56743.8463 -14102.2631 266 -14102.2631 56743.8463 267 60145.7968 -14102.2631 268 -14102.2631 60145.7968 269 35248.8380 -14102.2631 270 -14102.2631 35248.8380 271 -14102.2631 -14102.2631 272 -14102.2631 -14102.2631 273 55325.6036 -14102.2631 274 56743.8463 55325.6036 275 66053.6292 56743.8463 276 63962.5695 66053.6292 277 56743.8463 63962.5695 278 -25065.7169 56743.8463 279 35917.4441 -25065.7169 280 -14102.2631 35917.4441 281 56743.8463 -14102.2631 282 95798.9808 56743.8463 283 56743.8463 95798.9808 284 20948.2632 56743.8463 285 56743.8463 20948.2632 286 -14102.2631 56743.8463 287 55796.2699 -14102.2631 288 -14102.2631 55796.2699 289 87737.0651 -14102.2631 290 -14102.2631 87737.0651 291 66657.4778 -14102.2631 292 30150.0010 66657.4778 293 56743.8463 30150.0010 294 60459.4379 56743.8463 295 20666.8015 60459.4379 296 56743.8463 20666.8015 297 -14102.2631 56743.8463 298 56743.8463 -14102.2631 299 56743.8463 56743.8463 300 -14799.5421 56743.8463 301 42657.9409 -14799.5421 302 56743.8463 42657.9409 303 -14102.2631 56743.8463 304 96052.7763 -14102.2631 305 60905.0179 96052.7763 306 -14102.2631 60905.0179 307 -22026.4417 -14102.2631 308 78605.8994 -22026.4417 309 -14102.2631 78605.8994 310 -14102.2631 -14102.2631 311 37401.7567 -14102.2631 312 -14102.2631 37401.7567 313 -17024.1271 -14102.2631 314 58714.9646 -17024.1271 315 57989.7061 58714.9646 316 58333.6778 57989.7061 317 42709.8089 58333.6778 318 59542.2511 42709.8089 319 56743.8463 59542.2511 320 56743.8463 56743.8463 321 56743.8463 56743.8463 322 56743.8463 56743.8463 323 -14102.2631 56743.8463 324 56743.8463 -14102.2631 325 54184.7061 56743.8463 326 -14102.2631 54184.7061 327 -14102.2631 -14102.2631 328 56743.8463 -14102.2631 329 56743.8463 56743.8463 330 56743.8463 56743.8463 331 56743.8463 56743.8463 332 56743.8463 56743.8463 333 56785.5558 56743.8463 334 57641.3511 56785.5558 335 56743.8463 57641.3511 336 55622.0937 56743.8463 337 56743.8463 55622.0937 338 56098.6316 56743.8463 339 60825.7496 56098.6316 340 37739.2850 60825.7496 341 56743.8463 37739.2850 342 56743.8463 56743.8463 343 56743.8463 56743.8463 344 56743.8463 56743.8463 345 56743.8463 56743.8463 346 48733.4452 56743.8463 347 56743.8463 48733.4452 348 56743.8463 56743.8463 349 56743.8463 56743.8463 350 57507.4643 56743.8463 351 58135.5877 57507.4643 352 58182.2306 58135.5877 353 46857.2332 58182.2306 354 56743.8463 46857.2332 355 56743.8463 56743.8463 356 24153.8623 56743.8463 357 56743.8463 24153.8623 358 56743.8463 56743.8463 359 56743.8463 56743.8463 360 58592.2977 56743.8463 361 56743.8463 58592.2977 362 56743.8463 56743.8463 363 -23229.0752 56743.8463 364 55339.9264 -23229.0752 365 56743.8463 55339.9264 366 56743.8463 56743.8463 367 56743.8463 56743.8463 368 56743.8463 56743.8463 369 56743.8463 56743.8463 370 56743.8463 56743.8463 371 59951.3846 56743.8463 372 16950.6168 59951.3846 373 29944.4766 16950.6168 374 56743.8463 29944.4766 375 56743.8463 56743.8463 376 69183.7439 56743.8463 377 41529.1872 69183.7439 378 56743.8463 41529.1872 379 56743.8463 56743.8463 380 56743.8463 56743.8463 381 50304.4797 56743.8463 382 110203.6465 50304.4797 383 56743.8463 110203.6465 384 58971.6081 56743.8463 385 -99802.4349 58971.6081 386 61423.9876 -99802.4349 387 40725.4677 61423.9876 388 54932.2314 40725.4677 389 37401.4713 54932.2314 390 53556.4003 37401.4713 391 36181.4677 53556.4003 392 18933.7581 36181.4677 393 73314.7260 18933.7581 394 -225558.1571 73314.7260 395 29067.6811 -225558.1571 396 -107954.0035 29067.6811 397 -1003.1517 -107954.0035 398 30167.9153 -1003.1517 399 -20626.9915 30167.9153 400 18801.9662 -20626.9915 401 46542.0025 18801.9662 402 35992.8637 46542.0025 403 -96274.3599 35992.8637 404 23096.4155 -96274.3599 405 4588.9031 23096.4155 406 41079.1019 4588.9031 407 36064.8764 41079.1019 408 34340.5635 36064.8764 409 -18627.8954 34340.5635 410 12187.3687 -18627.8954 411 -219006.8707 12187.3687 412 -43752.6695 -219006.8707 413 -154815.6019 -43752.6695 414 -54634.8218 -154815.6019 415 -37896.3342 -54634.8218 416 138483.9631 -37896.3342 417 -67743.4506 138483.9631 418 5703.9594 -67743.4506 419 -11870.2290 5703.9594 420 -65172.7527 -11870.2290 421 -99395.9490 -65172.7527 422 -427200.0501 -99395.9490 423 -470445.0446 -427200.0501 424 -71049.1911 -470445.0446 425 -200544.6914 -71049.1911 426 -306223.7944 -200544.6914 427 -1206200.7325 -306223.7944 428 -359249.8705 -1206200.7325 429 -270057.1479 -359249.8705 430 -843491.5827 -270057.1479 431 NA -843491.5827 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 3267339.8005 1614392.9083 [2,] 1342318.8812 3267339.8005 [3,] -3145447.1126 1342318.8812 [4,] -392660.3716 -3145447.1126 [5,] 47198.0501 -392660.3716 [6,] 404.7394 47198.0501 [7,] 602188.5762 404.7394 [8,] -211744.7650 602188.5762 [9,] 769347.7162 -211744.7650 [10,] -324970.7246 769347.7162 [11,] 476487.5596 -324970.7246 [12,] 971510.2238 476487.5596 [13,] -1780770.7218 971510.2238 [14,] 52895.3919 -1780770.7218 [15,] 440956.3216 52895.3919 [16,] -411358.2358 440956.3216 [17,] -67512.6769 -411358.2358 [18,] -148150.4768 -67512.6769 [19,] -328203.7754 -148150.4768 [20,] 173412.6685 -328203.7754 [21,] -246220.5692 173412.6685 [22,] -237869.2592 -246220.5692 [23,] -488859.5214 -237869.2592 [24,] -214929.7252 -488859.5214 [25,] -127584.6004 -214929.7252 [26,] -365900.8069 -127584.6004 [27,] -441606.0856 -365900.8069 [28,] 404846.4220 -441606.0856 [29,] 24692.1759 404846.4220 [30,] 180623.9925 24692.1759 [31,] -162947.7947 180623.9925 [32,] 132677.0574 -162947.7947 [33,] 16731.1667 132677.0574 [34,] -338033.7242 16731.1667 [35,] -454918.3637 -338033.7242 [36,] 55071.1870 -454918.3637 [37,] 73640.4234 55071.1870 [38,] -122998.4009 73640.4234 [39,] -149060.1825 -122998.4009 [40,] 23536.2045 -149060.1825 [41,] -354630.1910 23536.2045 [42,] 103738.8678 -354630.1910 [43,] -186161.7635 103738.8678 [44,] -11873.8924 -186161.7635 [45,] -232458.0978 -11873.8924 [46,] -248121.3384 -232458.0978 [47,] -171890.4711 -248121.3384 [48,] -103951.1514 -171890.4711 [49,] -156677.2080 -103951.1514 [50,] 25768.3763 -156677.2080 [51,] -125901.9966 25768.3763 [52,] -405539.2242 -125901.9966 [53,] -48526.5126 -405539.2242 [54,] -394622.0142 -48526.5126 [55,] 134777.1577 -394622.0142 [56,] 182089.8545 134777.1577 [57,] -18512.2564 182089.8545 [58,] -20266.4483 -18512.2564 [59,] -226136.2158 -20266.4483 [60,] -89112.9730 -226136.2158 [61,] -57073.6173 -89112.9730 [62,] -109043.6851 -57073.6173 [63,] -232356.0789 -109043.6851 [64,] 65003.2838 -232356.0789 [65,] 81678.9589 65003.2838 [66,] -91906.2811 81678.9589 [67,] -99661.3241 -91906.2811 [68,] -58130.9141 -99661.3241 [69,] -136377.9104 -58130.9141 [70,] -28840.9276 -136377.9104 [71,] 73016.2831 -28840.9276 [72,] -125491.3547 73016.2831 [73,] -14991.2288 -125491.3547 [74,] -121785.5624 -14991.2288 [75,] -141157.5865 -121785.5624 [76,] -92252.6956 -141157.5865 [77,] -47277.0402 -92252.6956 [78,] 87516.0975 -47277.0402 [79,] -56250.2459 87516.0975 [80,] 51772.1918 -56250.2459 [81,] -98359.6297 51772.1918 [82,] -118980.2572 -98359.6297 [83,] 89337.8404 -118980.2572 [84,] -5467.0157 89337.8404 [85,] -111554.4237 -5467.0157 [86,] 30263.8277 -111554.4237 [87,] -39334.6250 30263.8277 [88,] -95805.4969 -39334.6250 [89,] -20137.4147 -95805.4969 [90,] -26562.2862 -20137.4147 [91,] -102847.2821 -26562.2862 [92,] 136445.8672 -102847.2821 [93,] -339389.1301 136445.8672 [94,] 142088.0632 -339389.1301 [95,] -22029.2939 142088.0632 [96,] 12387.1874 -22029.2939 [97,] 47889.5549 12387.1874 [98,] -132866.7012 47889.5549 [99,] -85702.5943 -132866.7012 [100,] -93984.4068 -85702.5943 [101,] -101209.0829 -93984.4068 [102,] -110607.5576 -101209.0829 [103,] -53863.5304 -110607.5576 [104,] 59661.3416 -53863.5304 [105,] -44169.0784 59661.3416 [106,] -239180.3303 -44169.0784 [107,] -1079.6495 -239180.3303 [108,] -203731.6259 -1079.6495 [109,] -152145.4789 -203731.6259 [110,] 116244.2019 -152145.4789 [111,] -234855.3016 116244.2019 [112,] -268536.7701 -234855.3016 [113,] -174382.2503 -268536.7701 [114,] -7876.3314 -174382.2503 [115,] 79594.8195 -7876.3314 [116,] 11626.7819 79594.8195 [117,] 40981.2854 11626.7819 [118,] -237174.4564 40981.2854 [119,] 77984.4846 -237174.4564 [120,] -260529.3804 77984.4846 [121,] 44126.2494 -260529.3804 [122,] 2437.0379 44126.2494 [123,] -48248.7533 2437.0379 [124,] -120768.4011 -48248.7533 [125,] -249849.1739 -120768.4011 [126,] 56743.8463 -249849.1739 [127,] 57048.3791 56743.8463 [128,] 64972.9222 57048.3791 [129,] 56743.8463 64972.9222 [130,] 59962.9990 56743.8463 [131,] 114381.9051 59962.9990 [132,] 56743.8463 114381.9051 [133,] 54723.3658 56743.8463 [134,] 56743.8463 54723.3658 [135,] 59317.8032 56743.8463 [136,] 59103.1536 59317.8032 [137,] 56743.8463 59103.1536 [138,] 56743.8463 56743.8463 [139,] 56743.8463 56743.8463 [140,] 55786.6489 56743.8463 [141,] 56743.8463 55786.6489 [142,] 51497.1868 56743.8463 [143,] 57019.2654 51497.1868 [144,] 63025.8180 57019.2654 [145,] 68710.2635 63025.8180 [146,] 56743.8463 68710.2635 [147,] 58723.7837 56743.8463 [148,] 56743.8463 58723.7837 [149,] 56178.8235 56743.8463 [150,] 47378.3959 56178.8235 [151,] 46154.5855 47378.3959 [152,] 56743.8463 46154.5855 [153,] 34483.4455 56743.8463 [154,] 56844.1304 34483.4455 [155,] 56743.8463 56844.1304 [156,] 56743.8463 56743.8463 [157,] 53573.6402 56743.8463 [158,] 63896.5616 53573.6402 [159,] 56743.8463 63896.5616 [160,] 56743.8463 56743.8463 [161,] 56743.8463 56743.8463 [162,] 56743.8463 56743.8463 [163,] 62546.0270 56743.8463 [164,] 57582.2679 62546.0270 [165,] 56743.8463 57582.2679 [166,] 56743.8463 56743.8463 [167,] 160360.4989 56743.8463 [168,] 56743.8463 160360.4989 [169,] 56743.8463 56743.8463 [170,] 56743.8463 56743.8463 [171,] 53760.2551 56743.8463 [172,] 56743.8463 53760.2551 [173,] 54511.6315 56743.8463 [174,] 56743.8463 54511.6315 [175,] 56801.5725 56743.8463 [176,] 56743.8463 56801.5725 [177,] -15353.1898 56743.8463 [178,] 56743.8463 -15353.1898 [179,] 56743.8463 56743.8463 [180,] 70989.7457 56743.8463 [181,] -37474.7381 70989.7457 [182,] 82509.7404 -37474.7381 [183,] -76005.3479 82509.7404 [184,] -14797.0131 -76005.3479 [185,] -21688.9005 -14797.0131 [186,] -8918.1814 -21688.9005 [187,] 53423.5954 -8918.1814 [188,] 17043.9417 53423.5954 [189,] 56743.8463 17043.9417 [190,] 61395.0315 56743.8463 [191,] 56743.8463 61395.0315 [192,] 58916.3004 56743.8463 [193,] 56743.8463 58916.3004 [194,] 56474.3785 56743.8463 [195,] -14102.2631 56474.3785 [196,] -150102.8416 -14102.2631 [197,] -14102.2631 -150102.8416 [198,] -14102.2631 -14102.2631 [199,] -92005.6434 -14102.2631 [200,] -14102.2631 -92005.6434 [201,] -14102.2631 -14102.2631 [202,] -14102.2631 -14102.2631 [203,] 14151.1709 -14102.2631 [204,] -100326.1605 14151.1709 [205,] -15827.4682 -100326.1605 [206,] -29262.9780 -15827.4682 [207,] -14102.2631 -29262.9780 [208,] 44819.1609 -14102.2631 [209,] -14102.2631 44819.1609 [210,] 56743.8463 -14102.2631 [211,] 56743.8463 56743.8463 [212,] 72687.1102 56743.8463 [213,] 34552.6110 72687.1102 [214,] 37988.8278 34552.6110 [215,] 56743.8463 37988.8278 [216,] 60320.3399 56743.8463 [217,] 61069.6158 60320.3399 [218,] 56743.8463 61069.6158 [219,] 61629.1210 56743.8463 [220,] 23935.1655 61629.1210 [221,] 56743.8463 23935.1655 [222,] 57658.7888 56743.8463 [223,] 56743.8463 57658.7888 [224,] 57060.4504 56743.8463 [225,] 61626.0995 57060.4504 [226,] 56743.8463 61626.0995 [227,] 56743.8463 56743.8463 [228,] 30572.5273 56743.8463 [229,] 56743.8463 30572.5273 [230,] 35628.8129 56743.8463 [231,] 56743.8463 35628.8129 [232,] 56743.8463 56743.8463 [233,] 58191.0580 56743.8463 [234,] 56743.8463 58191.0580 [235,] 57007.1910 56743.8463 [236,] 56743.8463 57007.1910 [237,] 59739.9250 56743.8463 [238,] 56743.8463 59739.9250 [239,] 56743.8463 56743.8463 [240,] 56743.8463 56743.8463 [241,] 59635.8630 56743.8463 [242,] 56743.8463 59635.8630 [243,] 61991.0854 56743.8463 [244,] -30514.6252 61991.0854 [245,] -14102.2631 -30514.6252 [246,] -14102.2631 -14102.2631 [247,] 48670.3819 -14102.2631 [248,] -14102.2631 48670.3819 [249,] -14102.2631 -14102.2631 [250,] 63819.4336 -14102.2631 [251,] -14102.2631 63819.4336 [252,] 33344.1466 -14102.2631 [253,] 37600.1858 33344.1466 [254,] 6243.3979 37600.1858 [255,] 56743.8463 6243.3979 [256,] 80097.7803 56743.8463 [257,] 96741.5415 80097.7803 [258,] -14102.2631 96741.5415 [259,] 56743.8463 -14102.2631 [260,] -14102.2631 56743.8463 [261,] -14102.2631 -14102.2631 [262,] 37163.9577 -14102.2631 [263,] 65359.4770 37163.9577 [264,] -14102.2631 65359.4770 [265,] 56743.8463 -14102.2631 [266,] -14102.2631 56743.8463 [267,] 60145.7968 -14102.2631 [268,] -14102.2631 60145.7968 [269,] 35248.8380 -14102.2631 [270,] -14102.2631 35248.8380 [271,] -14102.2631 -14102.2631 [272,] -14102.2631 -14102.2631 [273,] 55325.6036 -14102.2631 [274,] 56743.8463 55325.6036 [275,] 66053.6292 56743.8463 [276,] 63962.5695 66053.6292 [277,] 56743.8463 63962.5695 [278,] -25065.7169 56743.8463 [279,] 35917.4441 -25065.7169 [280,] -14102.2631 35917.4441 [281,] 56743.8463 -14102.2631 [282,] 95798.9808 56743.8463 [283,] 56743.8463 95798.9808 [284,] 20948.2632 56743.8463 [285,] 56743.8463 20948.2632 [286,] -14102.2631 56743.8463 [287,] 55796.2699 -14102.2631 [288,] -14102.2631 55796.2699 [289,] 87737.0651 -14102.2631 [290,] -14102.2631 87737.0651 [291,] 66657.4778 -14102.2631 [292,] 30150.0010 66657.4778 [293,] 56743.8463 30150.0010 [294,] 60459.4379 56743.8463 [295,] 20666.8015 60459.4379 [296,] 56743.8463 20666.8015 [297,] -14102.2631 56743.8463 [298,] 56743.8463 -14102.2631 [299,] 56743.8463 56743.8463 [300,] -14799.5421 56743.8463 [301,] 42657.9409 -14799.5421 [302,] 56743.8463 42657.9409 [303,] -14102.2631 56743.8463 [304,] 96052.7763 -14102.2631 [305,] 60905.0179 96052.7763 [306,] -14102.2631 60905.0179 [307,] -22026.4417 -14102.2631 [308,] 78605.8994 -22026.4417 [309,] -14102.2631 78605.8994 [310,] -14102.2631 -14102.2631 [311,] 37401.7567 -14102.2631 [312,] -14102.2631 37401.7567 [313,] -17024.1271 -14102.2631 [314,] 58714.9646 -17024.1271 [315,] 57989.7061 58714.9646 [316,] 58333.6778 57989.7061 [317,] 42709.8089 58333.6778 [318,] 59542.2511 42709.8089 [319,] 56743.8463 59542.2511 [320,] 56743.8463 56743.8463 [321,] 56743.8463 56743.8463 [322,] 56743.8463 56743.8463 [323,] -14102.2631 56743.8463 [324,] 56743.8463 -14102.2631 [325,] 54184.7061 56743.8463 [326,] -14102.2631 54184.7061 [327,] -14102.2631 -14102.2631 [328,] 56743.8463 -14102.2631 [329,] 56743.8463 56743.8463 [330,] 56743.8463 56743.8463 [331,] 56743.8463 56743.8463 [332,] 56743.8463 56743.8463 [333,] 56785.5558 56743.8463 [334,] 57641.3511 56785.5558 [335,] 56743.8463 57641.3511 [336,] 55622.0937 56743.8463 [337,] 56743.8463 55622.0937 [338,] 56098.6316 56743.8463 [339,] 60825.7496 56098.6316 [340,] 37739.2850 60825.7496 [341,] 56743.8463 37739.2850 [342,] 56743.8463 56743.8463 [343,] 56743.8463 56743.8463 [344,] 56743.8463 56743.8463 [345,] 56743.8463 56743.8463 [346,] 48733.4452 56743.8463 [347,] 56743.8463 48733.4452 [348,] 56743.8463 56743.8463 [349,] 56743.8463 56743.8463 [350,] 57507.4643 56743.8463 [351,] 58135.5877 57507.4643 [352,] 58182.2306 58135.5877 [353,] 46857.2332 58182.2306 [354,] 56743.8463 46857.2332 [355,] 56743.8463 56743.8463 [356,] 24153.8623 56743.8463 [357,] 56743.8463 24153.8623 [358,] 56743.8463 56743.8463 [359,] 56743.8463 56743.8463 [360,] 58592.2977 56743.8463 [361,] 56743.8463 58592.2977 [362,] 56743.8463 56743.8463 [363,] -23229.0752 56743.8463 [364,] 55339.9264 -23229.0752 [365,] 56743.8463 55339.9264 [366,] 56743.8463 56743.8463 [367,] 56743.8463 56743.8463 [368,] 56743.8463 56743.8463 [369,] 56743.8463 56743.8463 [370,] 56743.8463 56743.8463 [371,] 59951.3846 56743.8463 [372,] 16950.6168 59951.3846 [373,] 29944.4766 16950.6168 [374,] 56743.8463 29944.4766 [375,] 56743.8463 56743.8463 [376,] 69183.7439 56743.8463 [377,] 41529.1872 69183.7439 [378,] 56743.8463 41529.1872 [379,] 56743.8463 56743.8463 [380,] 56743.8463 56743.8463 [381,] 50304.4797 56743.8463 [382,] 110203.6465 50304.4797 [383,] 56743.8463 110203.6465 [384,] 58971.6081 56743.8463 [385,] -99802.4349 58971.6081 [386,] 61423.9876 -99802.4349 [387,] 40725.4677 61423.9876 [388,] 54932.2314 40725.4677 [389,] 37401.4713 54932.2314 [390,] 53556.4003 37401.4713 [391,] 36181.4677 53556.4003 [392,] 18933.7581 36181.4677 [393,] 73314.7260 18933.7581 [394,] -225558.1571 73314.7260 [395,] 29067.6811 -225558.1571 [396,] -107954.0035 29067.6811 [397,] -1003.1517 -107954.0035 [398,] 30167.9153 -1003.1517 [399,] -20626.9915 30167.9153 [400,] 18801.9662 -20626.9915 [401,] 46542.0025 18801.9662 [402,] 35992.8637 46542.0025 [403,] -96274.3599 35992.8637 [404,] 23096.4155 -96274.3599 [405,] 4588.9031 23096.4155 [406,] 41079.1019 4588.9031 [407,] 36064.8764 41079.1019 [408,] 34340.5635 36064.8764 [409,] -18627.8954 34340.5635 [410,] 12187.3687 -18627.8954 [411,] -219006.8707 12187.3687 [412,] -43752.6695 -219006.8707 [413,] -154815.6019 -43752.6695 [414,] -54634.8218 -154815.6019 [415,] -37896.3342 -54634.8218 [416,] 138483.9631 -37896.3342 [417,] -67743.4506 138483.9631 [418,] 5703.9594 -67743.4506 [419,] -11870.2290 5703.9594 [420,] -65172.7527 -11870.2290 [421,] -99395.9490 -65172.7527 [422,] -427200.0501 -99395.9490 [423,] -470445.0446 -427200.0501 [424,] -71049.1911 -470445.0446 [425,] -200544.6914 -71049.1911 [426,] -306223.7944 -200544.6914 [427,] -1206200.7325 -306223.7944 [428,] -359249.8705 -1206200.7325 [429,] -270057.1479 -359249.8705 [430,] -843491.5827 -270057.1479 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 3267339.8005 1614392.9083 2 1342318.8812 3267339.8005 3 -3145447.1126 1342318.8812 4 -392660.3716 -3145447.1126 5 47198.0501 -392660.3716 6 404.7394 47198.0501 7 602188.5762 404.7394 8 -211744.7650 602188.5762 9 769347.7162 -211744.7650 10 -324970.7246 769347.7162 11 476487.5596 -324970.7246 12 971510.2238 476487.5596 13 -1780770.7218 971510.2238 14 52895.3919 -1780770.7218 15 440956.3216 52895.3919 16 -411358.2358 440956.3216 17 -67512.6769 -411358.2358 18 -148150.4768 -67512.6769 19 -328203.7754 -148150.4768 20 173412.6685 -328203.7754 21 -246220.5692 173412.6685 22 -237869.2592 -246220.5692 23 -488859.5214 -237869.2592 24 -214929.7252 -488859.5214 25 -127584.6004 -214929.7252 26 -365900.8069 -127584.6004 27 -441606.0856 -365900.8069 28 404846.4220 -441606.0856 29 24692.1759 404846.4220 30 180623.9925 24692.1759 31 -162947.7947 180623.9925 32 132677.0574 -162947.7947 33 16731.1667 132677.0574 34 -338033.7242 16731.1667 35 -454918.3637 -338033.7242 36 55071.1870 -454918.3637 37 73640.4234 55071.1870 38 -122998.4009 73640.4234 39 -149060.1825 -122998.4009 40 23536.2045 -149060.1825 41 -354630.1910 23536.2045 42 103738.8678 -354630.1910 43 -186161.7635 103738.8678 44 -11873.8924 -186161.7635 45 -232458.0978 -11873.8924 46 -248121.3384 -232458.0978 47 -171890.4711 -248121.3384 48 -103951.1514 -171890.4711 49 -156677.2080 -103951.1514 50 25768.3763 -156677.2080 51 -125901.9966 25768.3763 52 -405539.2242 -125901.9966 53 -48526.5126 -405539.2242 54 -394622.0142 -48526.5126 55 134777.1577 -394622.0142 56 182089.8545 134777.1577 57 -18512.2564 182089.8545 58 -20266.4483 -18512.2564 59 -226136.2158 -20266.4483 60 -89112.9730 -226136.2158 61 -57073.6173 -89112.9730 62 -109043.6851 -57073.6173 63 -232356.0789 -109043.6851 64 65003.2838 -232356.0789 65 81678.9589 65003.2838 66 -91906.2811 81678.9589 67 -99661.3241 -91906.2811 68 -58130.9141 -99661.3241 69 -136377.9104 -58130.9141 70 -28840.9276 -136377.9104 71 73016.2831 -28840.9276 72 -125491.3547 73016.2831 73 -14991.2288 -125491.3547 74 -121785.5624 -14991.2288 75 -141157.5865 -121785.5624 76 -92252.6956 -141157.5865 77 -47277.0402 -92252.6956 78 87516.0975 -47277.0402 79 -56250.2459 87516.0975 80 51772.1918 -56250.2459 81 -98359.6297 51772.1918 82 -118980.2572 -98359.6297 83 89337.8404 -118980.2572 84 -5467.0157 89337.8404 85 -111554.4237 -5467.0157 86 30263.8277 -111554.4237 87 -39334.6250 30263.8277 88 -95805.4969 -39334.6250 89 -20137.4147 -95805.4969 90 -26562.2862 -20137.4147 91 -102847.2821 -26562.2862 92 136445.8672 -102847.2821 93 -339389.1301 136445.8672 94 142088.0632 -339389.1301 95 -22029.2939 142088.0632 96 12387.1874 -22029.2939 97 47889.5549 12387.1874 98 -132866.7012 47889.5549 99 -85702.5943 -132866.7012 100 -93984.4068 -85702.5943 101 -101209.0829 -93984.4068 102 -110607.5576 -101209.0829 103 -53863.5304 -110607.5576 104 59661.3416 -53863.5304 105 -44169.0784 59661.3416 106 -239180.3303 -44169.0784 107 -1079.6495 -239180.3303 108 -203731.6259 -1079.6495 109 -152145.4789 -203731.6259 110 116244.2019 -152145.4789 111 -234855.3016 116244.2019 112 -268536.7701 -234855.3016 113 -174382.2503 -268536.7701 114 -7876.3314 -174382.2503 115 79594.8195 -7876.3314 116 11626.7819 79594.8195 117 40981.2854 11626.7819 118 -237174.4564 40981.2854 119 77984.4846 -237174.4564 120 -260529.3804 77984.4846 121 44126.2494 -260529.3804 122 2437.0379 44126.2494 123 -48248.7533 2437.0379 124 -120768.4011 -48248.7533 125 -249849.1739 -120768.4011 126 56743.8463 -249849.1739 127 57048.3791 56743.8463 128 64972.9222 57048.3791 129 56743.8463 64972.9222 130 59962.9990 56743.8463 131 114381.9051 59962.9990 132 56743.8463 114381.9051 133 54723.3658 56743.8463 134 56743.8463 54723.3658 135 59317.8032 56743.8463 136 59103.1536 59317.8032 137 56743.8463 59103.1536 138 56743.8463 56743.8463 139 56743.8463 56743.8463 140 55786.6489 56743.8463 141 56743.8463 55786.6489 142 51497.1868 56743.8463 143 57019.2654 51497.1868 144 63025.8180 57019.2654 145 68710.2635 63025.8180 146 56743.8463 68710.2635 147 58723.7837 56743.8463 148 56743.8463 58723.7837 149 56178.8235 56743.8463 150 47378.3959 56178.8235 151 46154.5855 47378.3959 152 56743.8463 46154.5855 153 34483.4455 56743.8463 154 56844.1304 34483.4455 155 56743.8463 56844.1304 156 56743.8463 56743.8463 157 53573.6402 56743.8463 158 63896.5616 53573.6402 159 56743.8463 63896.5616 160 56743.8463 56743.8463 161 56743.8463 56743.8463 162 56743.8463 56743.8463 163 62546.0270 56743.8463 164 57582.2679 62546.0270 165 56743.8463 57582.2679 166 56743.8463 56743.8463 167 160360.4989 56743.8463 168 56743.8463 160360.4989 169 56743.8463 56743.8463 170 56743.8463 56743.8463 171 53760.2551 56743.8463 172 56743.8463 53760.2551 173 54511.6315 56743.8463 174 56743.8463 54511.6315 175 56801.5725 56743.8463 176 56743.8463 56801.5725 177 -15353.1898 56743.8463 178 56743.8463 -15353.1898 179 56743.8463 56743.8463 180 70989.7457 56743.8463 181 -37474.7381 70989.7457 182 82509.7404 -37474.7381 183 -76005.3479 82509.7404 184 -14797.0131 -76005.3479 185 -21688.9005 -14797.0131 186 -8918.1814 -21688.9005 187 53423.5954 -8918.1814 188 17043.9417 53423.5954 189 56743.8463 17043.9417 190 61395.0315 56743.8463 191 56743.8463 61395.0315 192 58916.3004 56743.8463 193 56743.8463 58916.3004 194 56474.3785 56743.8463 195 -14102.2631 56474.3785 196 -150102.8416 -14102.2631 197 -14102.2631 -150102.8416 198 -14102.2631 -14102.2631 199 -92005.6434 -14102.2631 200 -14102.2631 -92005.6434 201 -14102.2631 -14102.2631 202 -14102.2631 -14102.2631 203 14151.1709 -14102.2631 204 -100326.1605 14151.1709 205 -15827.4682 -100326.1605 206 -29262.9780 -15827.4682 207 -14102.2631 -29262.9780 208 44819.1609 -14102.2631 209 -14102.2631 44819.1609 210 56743.8463 -14102.2631 211 56743.8463 56743.8463 212 72687.1102 56743.8463 213 34552.6110 72687.1102 214 37988.8278 34552.6110 215 56743.8463 37988.8278 216 60320.3399 56743.8463 217 61069.6158 60320.3399 218 56743.8463 61069.6158 219 61629.1210 56743.8463 220 23935.1655 61629.1210 221 56743.8463 23935.1655 222 57658.7888 56743.8463 223 56743.8463 57658.7888 224 57060.4504 56743.8463 225 61626.0995 57060.4504 226 56743.8463 61626.0995 227 56743.8463 56743.8463 228 30572.5273 56743.8463 229 56743.8463 30572.5273 230 35628.8129 56743.8463 231 56743.8463 35628.8129 232 56743.8463 56743.8463 233 58191.0580 56743.8463 234 56743.8463 58191.0580 235 57007.1910 56743.8463 236 56743.8463 57007.1910 237 59739.9250 56743.8463 238 56743.8463 59739.9250 239 56743.8463 56743.8463 240 56743.8463 56743.8463 241 59635.8630 56743.8463 242 56743.8463 59635.8630 243 61991.0854 56743.8463 244 -30514.6252 61991.0854 245 -14102.2631 -30514.6252 246 -14102.2631 -14102.2631 247 48670.3819 -14102.2631 248 -14102.2631 48670.3819 249 -14102.2631 -14102.2631 250 63819.4336 -14102.2631 251 -14102.2631 63819.4336 252 33344.1466 -14102.2631 253 37600.1858 33344.1466 254 6243.3979 37600.1858 255 56743.8463 6243.3979 256 80097.7803 56743.8463 257 96741.5415 80097.7803 258 -14102.2631 96741.5415 259 56743.8463 -14102.2631 260 -14102.2631 56743.8463 261 -14102.2631 -14102.2631 262 37163.9577 -14102.2631 263 65359.4770 37163.9577 264 -14102.2631 65359.4770 265 56743.8463 -14102.2631 266 -14102.2631 56743.8463 267 60145.7968 -14102.2631 268 -14102.2631 60145.7968 269 35248.8380 -14102.2631 270 -14102.2631 35248.8380 271 -14102.2631 -14102.2631 272 -14102.2631 -14102.2631 273 55325.6036 -14102.2631 274 56743.8463 55325.6036 275 66053.6292 56743.8463 276 63962.5695 66053.6292 277 56743.8463 63962.5695 278 -25065.7169 56743.8463 279 35917.4441 -25065.7169 280 -14102.2631 35917.4441 281 56743.8463 -14102.2631 282 95798.9808 56743.8463 283 56743.8463 95798.9808 284 20948.2632 56743.8463 285 56743.8463 20948.2632 286 -14102.2631 56743.8463 287 55796.2699 -14102.2631 288 -14102.2631 55796.2699 289 87737.0651 -14102.2631 290 -14102.2631 87737.0651 291 66657.4778 -14102.2631 292 30150.0010 66657.4778 293 56743.8463 30150.0010 294 60459.4379 56743.8463 295 20666.8015 60459.4379 296 56743.8463 20666.8015 297 -14102.2631 56743.8463 298 56743.8463 -14102.2631 299 56743.8463 56743.8463 300 -14799.5421 56743.8463 301 42657.9409 -14799.5421 302 56743.8463 42657.9409 303 -14102.2631 56743.8463 304 96052.7763 -14102.2631 305 60905.0179 96052.7763 306 -14102.2631 60905.0179 307 -22026.4417 -14102.2631 308 78605.8994 -22026.4417 309 -14102.2631 78605.8994 310 -14102.2631 -14102.2631 311 37401.7567 -14102.2631 312 -14102.2631 37401.7567 313 -17024.1271 -14102.2631 314 58714.9646 -17024.1271 315 57989.7061 58714.9646 316 58333.6778 57989.7061 317 42709.8089 58333.6778 318 59542.2511 42709.8089 319 56743.8463 59542.2511 320 56743.8463 56743.8463 321 56743.8463 56743.8463 322 56743.8463 56743.8463 323 -14102.2631 56743.8463 324 56743.8463 -14102.2631 325 54184.7061 56743.8463 326 -14102.2631 54184.7061 327 -14102.2631 -14102.2631 328 56743.8463 -14102.2631 329 56743.8463 56743.8463 330 56743.8463 56743.8463 331 56743.8463 56743.8463 332 56743.8463 56743.8463 333 56785.5558 56743.8463 334 57641.3511 56785.5558 335 56743.8463 57641.3511 336 55622.0937 56743.8463 337 56743.8463 55622.0937 338 56098.6316 56743.8463 339 60825.7496 56098.6316 340 37739.2850 60825.7496 341 56743.8463 37739.2850 342 56743.8463 56743.8463 343 56743.8463 56743.8463 344 56743.8463 56743.8463 345 56743.8463 56743.8463 346 48733.4452 56743.8463 347 56743.8463 48733.4452 348 56743.8463 56743.8463 349 56743.8463 56743.8463 350 57507.4643 56743.8463 351 58135.5877 57507.4643 352 58182.2306 58135.5877 353 46857.2332 58182.2306 354 56743.8463 46857.2332 355 56743.8463 56743.8463 356 24153.8623 56743.8463 357 56743.8463 24153.8623 358 56743.8463 56743.8463 359 56743.8463 56743.8463 360 58592.2977 56743.8463 361 56743.8463 58592.2977 362 56743.8463 56743.8463 363 -23229.0752 56743.8463 364 55339.9264 -23229.0752 365 56743.8463 55339.9264 366 56743.8463 56743.8463 367 56743.8463 56743.8463 368 56743.8463 56743.8463 369 56743.8463 56743.8463 370 56743.8463 56743.8463 371 59951.3846 56743.8463 372 16950.6168 59951.3846 373 29944.4766 16950.6168 374 56743.8463 29944.4766 375 56743.8463 56743.8463 376 69183.7439 56743.8463 377 41529.1872 69183.7439 378 56743.8463 41529.1872 379 56743.8463 56743.8463 380 56743.8463 56743.8463 381 50304.4797 56743.8463 382 110203.6465 50304.4797 383 56743.8463 110203.6465 384 58971.6081 56743.8463 385 -99802.4349 58971.6081 386 61423.9876 -99802.4349 387 40725.4677 61423.9876 388 54932.2314 40725.4677 389 37401.4713 54932.2314 390 53556.4003 37401.4713 391 36181.4677 53556.4003 392 18933.7581 36181.4677 393 73314.7260 18933.7581 394 -225558.1571 73314.7260 395 29067.6811 -225558.1571 396 -107954.0035 29067.6811 397 -1003.1517 -107954.0035 398 30167.9153 -1003.1517 399 -20626.9915 30167.9153 400 18801.9662 -20626.9915 401 46542.0025 18801.9662 402 35992.8637 46542.0025 403 -96274.3599 35992.8637 404 23096.4155 -96274.3599 405 4588.9031 23096.4155 406 41079.1019 4588.9031 407 36064.8764 41079.1019 408 34340.5635 36064.8764 409 -18627.8954 34340.5635 410 12187.3687 -18627.8954 411 -219006.8707 12187.3687 412 -43752.6695 -219006.8707 413 -154815.6019 -43752.6695 414 -54634.8218 -154815.6019 415 -37896.3342 -54634.8218 416 138483.9631 -37896.3342 417 -67743.4506 138483.9631 418 5703.9594 -67743.4506 419 -11870.2290 5703.9594 420 -65172.7527 -11870.2290 421 -99395.9490 -65172.7527 422 -427200.0501 -99395.9490 423 -470445.0446 -427200.0501 424 -71049.1911 -470445.0446 425 -200544.6914 -71049.1911 426 -306223.7944 -200544.6914 427 -1206200.7325 -306223.7944 428 -359249.8705 -1206200.7325 429 -270057.1479 -359249.8705 430 -843491.5827 -270057.1479 > 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/7okqn1291302152.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/8okqn1291302152.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/9okqn1291302152.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/10ztqq1291302152.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/112c6e1291302152.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/12ncnk1291302152.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/1314kt1291302152.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/14xf4u1291302153.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/15863w1291302153.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/164f1n1291302153.tab") + } > > try(system("convert tmp/19sbw1291302152.ps tmp/19sbw1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/29sbw1291302152.ps tmp/29sbw1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/3kjsz1291302152.ps tmp/3kjsz1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/4kjsz1291302152.ps tmp/4kjsz1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/5kjsz1291302152.ps tmp/5kjsz1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/6dt9k1291302152.ps tmp/6dt9k1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/7okqn1291302152.ps tmp/7okqn1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/8okqn1291302152.ps tmp/8okqn1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/9okqn1291302152.ps tmp/9okqn1291302152.png",intern=TRUE)) character(0) > try(system("convert tmp/10ztqq1291302152.ps tmp/10ztqq1291302152.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.018 2.217 37.326