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(1 + ,162556 + ,162556 + ,213118 + ,213118 + ,6282929 + ,1 + ,29790 + ,29790 + ,81767 + ,81767 + ,4324047 + ,1 + ,87550 + ,87550 + ,153198 + ,153198 + ,4108272 + ,0 + ,84738 + ,0 + ,-26007 + ,0 + ,-1212617 + ,1 + ,54660 + ,54660 + ,126942 + ,126942 + ,1485329 + ,1 + ,42634 + ,42634 + ,157214 + ,157214 + ,1779876 + ,0 + ,40949 + ,0 + ,129352 + ,0 + ,1367203 + ,1 + ,42312 + ,42312 + ,234817 + ,234817 + ,2519076 + ,1 + ,37704 + ,37704 + ,60448 + ,60448 + ,912684 + ,1 + ,16275 + ,16275 + ,47818 + ,47818 + ,1443586 + ,0 + ,25830 + ,0 + ,245546 + ,0 + ,1220017 + ,0 + ,12679 + ,0 + ,48020 + ,0 + ,984885 + ,1 + ,18014 + ,18014 + ,-1710 + ,-1710 + ,1457425 + ,0 + ,43556 + ,0 + ,32648 + ,0 + ,-572920 + ,1 + ,24524 + ,24524 + ,95350 + ,95350 + ,929144 + ,0 + ,6532 + ,0 + ,151352 + ,0 + ,1151176 + ,0 + ,7123 + ,0 + ,288170 + ,0 + ,790090 + ,1 + ,20813 + ,20813 + ,114337 + ,114337 + ,774497 + ,1 + ,37597 + ,37597 + ,37884 + ,37884 + ,990576 + ,0 + ,17821 + ,0 + ,122844 + ,0 + ,454195 + ,1 + ,12988 + ,12988 + ,82340 + ,82340 + ,876607 + ,1 + ,22330 + ,22330 + ,79801 + ,79801 + ,711969 + ,0 + ,13326 + ,0 + ,165548 + ,0 + ,702380 + ,0 + ,16189 + ,0 + ,116384 + ,0 + ,264449 + ,0 + ,7146 + ,0 + ,134028 + ,0 + ,450033 + ,0 + ,15824 + ,0 + ,63838 + ,0 + ,541063 + ,1 + ,26088 + ,26088 + ,74996 + ,74996 + ,588864 + ,0 + ,11326 + ,0 + ,31080 + ,0 + ,-37216 + ,0 + ,8568 + ,0 + ,32168 + ,0 + ,783310 + ,0 + ,14416 + ,0 + ,49857 + ,0 + ,467359 + ,1 + ,3369 + ,3369 + ,87161 + ,87161 + ,688779 + ,1 + ,11819 + ,11819 + ,106113 + ,106113 + ,608419 + ,1 + ,6620 + ,6620 + ,80570 + ,80570 + ,696348 + ,1 + ,4519 + ,4519 + ,102129 + ,102129 + ,597793 + ,0 + ,2220 + ,0 + ,301670 + ,0 + ,821730 + ,0 + ,18562 + ,0 + ,102313 + ,0 + ,377934 + ,0 + ,10327 + ,0 + ,88577 + ,0 + ,651939 + ,1 + ,5336 + ,5336 + ,112477 + ,112477 + ,697458 + ,1 + ,2365 + ,2365 + ,191778 + ,191778 + ,700368 + ,0 + ,4069 + ,0 + ,79804 + ,0 + ,225986 + ,0 + ,7710 + ,0 + ,128294 + ,0 + ,348695 + ,0 + ,13718 + ,0 + ,96448 + ,0 + ,373683 + ,0 + ,4525 + ,0 + ,93811 + ,0 + ,501709 + ,0 + ,6869 + ,0 + ,117520 + ,0 + ,413743 + ,0 + ,4628 + ,0 + ,69159 + ,0 + ,379825 + ,1 + ,3653 + ,3653 + ,101792 + ,101792 + ,336260 + ,1 + ,1265 + ,1265 + ,210568 + ,210568 + ,636765 + ,1 + ,7489 + ,7489 + ,136996 + ,136996 + ,481231 + ,0 + ,4901 + ,0 + ,121920 + ,0 + ,469107 + ,0 + ,2284 + ,0 + ,76403 + ,0 + ,211928 + ,1 + ,3160 + ,3160 + ,108094 + ,108094 + ,563925 + ,1 + ,4150 + ,4150 + ,134759 + ,134759 + ,511939 + ,1 + ,7285 + ,7285 + ,188873 + ,188873 + ,521016 + ,1 + ,1134 + ,1134 + ,146216 + ,146216 + ,543856 + ,1 + ,4658 + ,4658 + ,156608 + ,156608 + ,329304 + ,0 + ,2384 + ,0 + ,61348 + ,0 + ,423262 + ,0 + ,3748 + ,0 + ,50350 + ,0 + ,509665 + ,0 + ,5371 + ,0 + ,87720 + ,0 + ,455881 + ,0 + ,1285 + ,0 + ,99489 + ,0 + ,367772 + ,1 + ,9327 + ,9327 + ,87419 + ,87419 + ,406339 + ,1 + ,5565 + ,5565 + ,94355 + ,94355 + ,493408 + ,0 + ,1528 + ,0 + ,60326 + ,0 + ,232942 + ,1 + ,3122 + ,3122 + ,94670 + ,94670 + ,416002 + ,1 + ,7317 + ,7317 + ,82425 + ,82425 + ,337430 + ,0 + ,2675 + ,0 + ,59017 + ,0 + ,361517 + ,0 + ,13253 + ,0 + ,90829 + ,0 + ,360962 + ,0 + ,880 + ,0 + ,80791 + ,0 + ,235561 + ,1 + ,2053 + ,2053 + ,100423 + ,100423 + ,408247 + ,0 + ,1424 + ,0 + ,131116 + ,0 + ,450296 + ,1 + ,4036 + ,4036 + ,100269 + ,100269 + ,418799 + ,1 + ,3045 + ,3045 + ,27330 + ,27330 + ,247405 + ,0 + ,5119 + ,0 + ,39039 + ,0 + ,378519 + ,0 + ,1431 + ,0 + ,106885 + ,0 + ,326638 + ,0 + ,554 + ,0 + ,79285 + ,0 + ,328233 + ,0 + ,1975 + ,0 + ,118881 + ,0 + ,386225 + ,1 + ,1286 + ,1286 + ,77623 + ,77623 + ,283662 + ,0 + ,1012 + ,0 + ,114768 + ,0 + ,370225 + ,0 + ,810 + ,0 + ,74015 + ,0 + ,269236 + ,0 + ,1280 + ,0 + ,69465 + ,0 + ,365732 + ,1 + ,666 + ,666 + ,117869 + ,117869 + ,420383 + ,0 + ,1380 + ,0 + ,60982 + ,0 + ,345811 + ,1 + ,4608 + ,4608 + ,90131 + ,90131 + ,431809 + ,0 + ,876 + ,0 + ,138971 + ,0 + ,418876 + ,0 + ,814 + ,0 + ,39625 + ,0 + ,297476 + ,0 + ,514 + ,0 + ,102725 + ,0 + ,416776 + ,1 + ,5692 + ,5692 + ,64239 + ,64239 + ,357257 + ,0 + ,3642 + ,0 + ,90262 + ,0 + ,458343 + ,0 + ,540 + ,0 + ,103960 + ,0 + ,388386 + ,0 + ,2099 + ,0 + ,106611 + ,0 + ,358934 + ,0 + ,567 + ,0 + ,103345 + ,0 + ,407560 + ,0 + ,2001 + ,0 + ,95551 + ,0 + ,392558 + ,1 + ,2949 + ,2949 + ,82903 + ,82903 + ,373177 + ,0 + ,2253 + ,0 + ,63593 + ,0 + ,428370 + ,1 + ,6533 + ,6533 + ,126910 + ,126910 + ,369419 + ,0 + ,1889 + ,0 + ,37527 + ,0 + ,358649 + ,1 + ,3055 + ,3055 + ,60247 + ,60247 + ,376641 + ,0 + ,272 + ,0 + ,112995 + ,0 + ,467427 + ,1 + ,1414 + ,1414 + ,70184 + ,70184 + ,364885 + ,0 + ,2564 + ,0 + ,130140 + ,0 + ,436230 + ,1 + ,1383 + ,1383 + ,73221 + ,73221 + ,329118 + ,1 + ,1261 + ,1261 + ,76114 + ,76114 + ,317365 + ,0 + ,975 + ,0 + ,90534 + ,0 + ,286849 + ,0 + ,3366 + ,0 + ,108479 + ,0 + ,376685 + ,0 + ,576 + ,0 + ,113761 + ,0 + ,407198 + ,0 + ,1306 + ,0 + ,68696 + ,0 + ,377772 + ,0 + ,746 + ,0 + ,71561 + ,0 + ,271483 + ,1 + ,3192 + ,3192 + ,59831 + ,59831 + ,153661 + ,1 + ,2045 + ,2045 + ,97890 + ,97890 + ,513294 + ,0 + ,5477 + ,0 + ,101481 + ,0 + ,324881 + ,1 + ,1932 + ,1932 + ,72954 + ,72954 + ,264512 + ,0 + ,936 + ,0 + ,67939 + ,0 + ,420968 + ,1 + ,3437 + ,3437 + ,48022 + ,48022 + ,129302 + ,0 + ,5131 + ,0 + ,86111 + ,0 + ,191521 + ,1 + ,2397 + ,2397 + ,74020 + ,74020 + ,268673 + ,1 + ,1389 + ,1389 + ,57530 + ,57530 + ,353179 + ,0 + ,1503 + ,0 + ,56364 + ,0 + ,354624 + ,0 + ,402 + ,0 + ,84990 + ,0 + ,363713 + ,0 + ,2239 + ,0 + ,88590 + ,0 + ,456657 + ,1 + ,2234 + ,2234 + ,77200 + ,77200 + ,211742 + ,0 + ,837 + ,0 + ,61262 + ,0 + ,338381 + ,0 + ,10579 + ,0 + ,110309 + ,0 + ,418530 + ,0 + ,875 + ,0 + ,67000 + ,0 + ,351483 + ,0 + ,1395 + ,0 + ,93099 + ,0 + ,372928 + ,1 + ,1659 + ,1659 + ,107577 + ,107577 + ,485538 + ,1 + ,2647 + ,2647 + ,62920 + ,62920 + ,279268 + ,1 + ,3294 + ,3294 + ,75832 + ,75832 + ,219060 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,94 + ,0 + ,60793 + ,0 + ,325314 + ,0 + ,422 + ,0 + ,57935 + ,0 + ,322046 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,34 + ,0 + ,60630 + ,0 + ,325599 + ,0 + ,1558 + ,0 + ,55637 + ,0 + ,377028 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,43 + ,0 + ,60887 + ,0 + ,323850 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,316 + ,0 + ,60505 + ,0 + ,331514 + ,0 + ,115 + ,0 + ,60945 + ,0 + ,325632 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,389 + ,0 + ,58990 + ,0 + ,322265 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,1002 + ,0 + ,56750 + ,0 + ,325906 + ,0 + ,36 + ,0 + ,60894 + ,0 + ,325985 + ,0 + ,460 + ,0 + ,63346 + ,0 + ,346145 + ,0 + ,309 + ,0 + ,56535 + ,0 + ,325898 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,9 + ,0 + ,60835 + ,0 + ,325356 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,14 + ,0 + ,61016 + ,0 + ,325930 + ,0 + ,520 + ,0 + ,58650 + ,0 + ,318020 + ,0 + ,1766 + ,0 + ,60438 + ,0 + ,326389 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,458 + ,0 + ,58625 + ,0 + ,302925 + ,0 + ,20 + ,0 + ,60938 + ,0 + ,325540 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,98 + ,0 + ,61490 + ,0 + ,326736 + ,0 + ,405 + ,0 + ,60845 + ,0 + ,340580 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,483 + ,0 + ,60830 + ,0 + ,331828 + ,0 + ,454 + ,0 + ,63261 + ,0 + ,323299 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,757 + ,0 + ,45689 + ,0 + ,387722 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,36 + ,0 + ,61564 + ,0 + ,324598 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,203 + ,0 + ,61938 + ,0 + ,328726 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,90 + ,0 + ,60951 + ,0 + ,325043 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,71 + ,71 + ,60745 + ,60745 + ,325806 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,972 + ,0 + ,71642 + ,0 + ,387732 + ,1 + ,531 + ,531 + ,71641 + ,71641 + ,349729 + ,0 + ,604 + ,0 + ,55792 + ,0 + ,332202 + ,1 + ,283 + ,283 + ,71873 + ,71873 + ,305442 + ,1 + ,23 + ,23 + ,62555 + ,62555 + ,329537 + ,1 + ,638 + ,638 + ,60370 + ,60370 + ,327055 + ,1 + ,699 + ,699 + ,64873 + ,64873 + ,356245 + ,0 + ,149 + ,0 + ,62041 + ,0 + ,328451 + ,0 + ,226 + ,0 + ,65745 + ,0 + ,307062 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,275 + ,0 + ,59500 + ,0 + ,331345 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,141 + ,0 + ,61630 + ,0 + ,331824 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,28 + ,0 + ,60890 + ,0 + ,325685 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,2566 + ,2566 + ,113521 + ,113521 + ,404480 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,472 + ,472 + ,80045 + ,80045 + ,318314 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,203 + ,203 + ,50804 + ,50804 + ,311807 + ,1 + ,496 + ,496 + ,87390 + ,87390 + ,337724 + ,1 + ,10 + ,10 + ,61656 + ,61656 + ,326431 + ,1 + ,63 + ,63 + ,65688 + ,65688 + ,327556 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,1136 + ,1136 + ,48522 + ,48522 + ,356850 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,267 + ,0 + ,57640 + ,0 + ,322741 + ,0 + ,474 + ,0 + ,61977 + ,0 + ,310902 + ,0 + ,534 + ,0 + ,62620 + ,0 + ,324295 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,15 + ,0 + ,60831 + ,0 + ,326156 + ,0 + ,397 + ,0 + ,60646 + ,0 + ,326960 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,1061 + ,0 + ,56225 + ,0 + ,333411 + ,0 + ,288 + ,0 + ,60510 + ,0 + ,297761 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,3 + ,0 + ,60698 + ,0 + ,325536 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,20 + ,0 + ,60805 + ,0 + ,325762 + ,0 + ,278 + ,0 + ,61404 + ,0 + ,327957 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,192 + ,0 + ,65276 + ,0 + ,318521 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,317 + ,0 + ,63915 + ,0 + ,319775 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,368 + ,0 + ,61686 + ,0 + ,332128 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,2 + ,0 + ,60743 + ,0 + ,325486 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,53 + ,0 + ,60349 + ,0 + ,325838 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,94 + ,0 + ,61360 + ,0 + ,331767 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,24 + ,0 + ,59818 + ,0 + ,324523 + ,0 + ,2332 + ,0 + ,72680 + ,0 + ,339995 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,131 + ,0 + ,61808 + ,0 + ,319582 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,206 + ,0 + ,53110 + ,0 + ,307245 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,167 + ,0 + ,64245 + ,0 + ,317967 + ,0 + ,622 + ,0 + ,73007 + ,0 + ,331488 + ,0 + ,885 + ,0 + ,82732 + ,0 + ,335452 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,365 + ,0 + ,54820 + ,0 + ,334184 + ,0 + ,364 + ,0 + ,47705 + ,0 + ,313213 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,226 + ,0 + ,72835 + ,0 + ,348678 + ,0 + ,307 + ,0 + ,58856 + ,0 + ,328727 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,188 + ,0 + ,77655 + ,0 + ,387978 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,138 + ,0 + ,69817 + ,0 + ,336704 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,125 + ,0 + ,60798 + ,0 + ,322076 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,282 + ,0 + ,62452 + ,0 + ,334272 + ,0 + ,335 + ,0 + ,64175 + ,0 + ,338197 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,813 + ,813 + ,67440 + ,67440 + ,321024 + ,0 + ,176 + ,0 + ,68136 + ,0 + ,322145 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,249 + ,0 + ,56726 + ,0 + ,323351 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,333 + ,0 + ,70811 + ,0 + ,327748 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,30 + ,0 + ,62045 + ,0 + ,328157 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,249 + ,0 + ,54323 + ,0 + ,311594 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,165 + ,0 + ,62841 + ,0 + ,335962 + ,0 + ,453 + ,0 + ,81125 + ,0 + ,372426 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,53 + ,0 + ,59506 + ,0 + ,319844 + ,1 + ,382 + ,382 + ,59365 + ,59365 + ,355822 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,30 + ,30 + ,60798 + ,60798 + ,324047 + ,0 + ,290 + ,0 + ,58790 + ,0 + ,311464 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,366 + ,0 + ,61808 + ,0 + ,353417 + ,0 + ,2 + ,0 + ,60735 + ,0 + ,325590 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,209 + ,209 + ,64016 + ,64016 + ,328576 + ,0 + ,384 + ,0 + ,54683 + ,0 + ,326126 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,365 + ,0 + ,87192 + ,0 + ,369376 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,49 + ,49 + ,64107 + ,64107 + ,332013 + ,0 + ,3 + ,0 + ,60761 + ,0 + ,325871 + ,0 + ,133 + ,0 + ,65990 + ,0 + ,342165 + ,0 + ,32 + ,0 + ,59988 + ,0 + ,324967 + ,0 + ,368 + ,0 + ,61167 + ,0 + ,314832 + ,0 + ,1 + ,0 + ,60719 + ,0 + ,325557 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,22 + ,0 + ,60722 + ,0 + ,322649 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,1 + ,0 + ,0 + ,60720 + ,60720 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,96 + ,0 + ,60379 + ,0 + ,324598 + ,0 + ,1 + ,0 + ,60727 + ,0 + ,325567 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,81 + ,0 + ,60925 + ,0 + ,324005 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,26 + ,0 + ,60896 + ,0 + ,325748 + ,0 + ,125 + ,0 + ,59734 + ,0 + ,323385 + ,0 + ,304 + ,0 + ,62969 + ,0 + ,315409 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,119 + ,0 + ,59118 + ,0 + ,312275 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,312 + ,0 + ,58598 + ,0 + ,320576 + ,0 + ,60 + ,0 + ,61124 + ,0 + ,325246 + ,0 + ,587 + ,0 + ,59595 + ,0 + ,332961 + ,0 + ,135 + ,0 + ,62065 + ,0 + ,323010 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,514 + ,0 + ,78780 + ,0 + ,345253 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,1 + ,0 + ,60722 + ,0 + ,325559 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,1 + ,58 + ,58 + ,61600 + ,61600 + ,319634 + ,0 + ,180 + ,0 + ,59635 + ,0 + ,319951 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,448 + ,0 + ,59781 + ,0 + ,318519 + ,0 + ,227 + ,0 + ,76644 + ,0 + ,343222 + ,0 + ,174 + ,0 + ,64820 + ,0 + ,317234 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,121 + ,0 + ,56178 + ,0 + ,314025 + ,0 + ,607 + ,0 + ,60436 + ,0 + ,320249 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,530 + ,0 + ,73433 + ,0 + ,349365 + ,0 + ,571 + ,0 + ,41477 + ,0 + ,289197 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,325560 + ,0 + ,78 + ,0 + ,62700 + ,0 + ,329245 + ,0 + ,2489 + ,0 + ,67804 + ,0 + ,240869 + ,0 + ,131 + ,0 + ,59661 + ,0 + ,327182 + ,0 + ,923 + ,0 + ,58620 + ,0 + ,322876 + ,0 + ,72 + ,0 + ,60398 + ,0 + ,323117 + ,0 + ,572 + ,0 + ,58580 + ,0 + ,306351 + ,0 + ,397 + ,0 + ,62710 + ,0 + ,335137 + ,0 + ,450 + ,0 + ,59325 + ,0 + ,308271 + ,0 + ,622 + ,0 + ,60950 + ,0 + ,301731 + ,0 + ,694 + ,0 + ,68060 + ,0 + ,382409 + ,1 + ,3425 + ,3425 + ,83620 + ,83620 + ,279230 + ,0 + ,562 + ,0 + ,58456 + ,0 + ,298731 + ,0 + ,4917 + ,0 + ,52811 + ,0 + ,243650 + ,1 + ,1442 + ,1442 + ,121173 + ,121173 + ,532682 + ,0 + ,529 + ,0 + ,63870 + ,0 + ,319771 + ,1 + ,2126 + ,2126 + ,21001 + ,21001 + ,171493 + ,0 + ,1061 + ,0 + ,70415 + ,0 + ,347262 + ,0 + ,776 + ,0 + ,64230 + ,0 + ,343945 + ,0 + ,611 + ,0 + ,59190 + ,0 + ,311874 + ,1 + ,1526 + ,1526 + ,69351 + ,69351 + ,302211 + ,0 + ,592 + ,0 + ,64270 + ,0 + ,316708 + ,0 + ,1182 + ,0 + ,70694 + ,0 + ,333463 + ,0 + ,621 + ,0 + ,68005 + ,0 + ,344282 + ,0 + ,989 + ,0 + ,58930 + ,0 + ,319635 + ,0 + ,438 + ,0 + ,58320 + ,0 + ,301186 + ,0 + ,726 + ,0 + ,69980 + ,0 + ,300381 + ,0 + ,1303 + ,0 + ,69863 + ,0 + ,318765 + ,1 + ,6341 + ,6341 + ,63255 + ,63255 + ,286146 + ,1 + ,1164 + ,1164 + ,57320 + ,57320 + ,306844 + ,1 + ,3310 + ,3310 + ,75230 + ,75230 + ,307705 + ,0 + ,1366 + ,0 + ,79420 + ,0 + ,312448 + ,0 + ,965 + ,0 + ,73490 + ,0 + ,299715 + ,0 + ,3256 + ,0 + ,35250 + ,0 + ,373399 + ,1 + ,1135 + ,1135 + ,62285 + ,62285 + ,299446 + ,0 + ,1270 + ,0 + ,69206 + ,0 + ,325586 + ,0 + ,661 + ,0 + ,65920 + ,0 + ,291221 + ,0 + ,1013 + ,0 + ,69770 + ,0 + ,261173 + ,0 + ,2844 + ,0 + ,72683 + ,0 + ,255027 + ,1 + ,11528 + ,11528 + ,-14545 + ,-14545 + ,-78375 + ,0 + ,6526 + ,0 + ,55830 + ,0 + ,-58143 + ,0 + ,2264 + ,0 + ,55174 + ,0 + ,227033 + ,1 + ,4461 + ,4461 + ,67038 + ,67038 + ,235098 + ,0 + ,3999 + ,0 + ,51252 + ,0 + ,21267 + ,0 + ,35624 + ,0 + ,157278 + ,0 + ,238675 + ,0 + ,9252 + ,0 + ,79510 + ,0 + ,197687 + ,0 + ,15236 + ,0 + ,77440 + ,0 + ,418341 + ,0 + ,18073 + ,0 + ,27284 + ,0 + ,-297706) + ,dim=c(6 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'GrCosts' + ,'Dividends' + ,'GrDiv' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','GrCosts','Dividends','GrDiv','Wealth'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '6' > #'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 GrCosts Dividends GrDiv 1 6282929 1 162556 162556 213118 213118 2 4324047 1 29790 29790 81767 81767 3 4108272 1 87550 87550 153198 153198 4 -1212617 0 84738 0 -26007 0 5 1485329 1 54660 54660 126942 126942 6 1779876 1 42634 42634 157214 157214 7 1367203 0 40949 0 129352 0 8 2519076 1 42312 42312 234817 234817 9 912684 1 37704 37704 60448 60448 10 1443586 1 16275 16275 47818 47818 11 1220017 0 25830 0 245546 0 12 984885 0 12679 0 48020 0 13 1457425 1 18014 18014 -1710 -1710 14 -572920 0 43556 0 32648 0 15 929144 1 24524 24524 95350 95350 16 1151176 0 6532 0 151352 0 17 790090 0 7123 0 288170 0 18 774497 1 20813 20813 114337 114337 19 990576 1 37597 37597 37884 37884 20 454195 0 17821 0 122844 0 21 876607 1 12988 12988 82340 82340 22 711969 1 22330 22330 79801 79801 23 702380 0 13326 0 165548 0 24 264449 0 16189 0 116384 0 25 450033 0 7146 0 134028 0 26 541063 0 15824 0 63838 0 27 588864 1 26088 26088 74996 74996 28 -37216 0 11326 0 31080 0 29 783310 0 8568 0 32168 0 30 467359 0 14416 0 49857 0 31 688779 1 3369 3369 87161 87161 32 608419 1 11819 11819 106113 106113 33 696348 1 6620 6620 80570 80570 34 597793 1 4519 4519 102129 102129 35 821730 0 2220 0 301670 0 36 377934 0 18562 0 102313 0 37 651939 0 10327 0 88577 0 38 697458 1 5336 5336 112477 112477 39 700368 1 2365 2365 191778 191778 40 225986 0 4069 0 79804 0 41 348695 0 7710 0 128294 0 42 373683 0 13718 0 96448 0 43 501709 0 4525 0 93811 0 44 413743 0 6869 0 117520 0 45 379825 0 4628 0 69159 0 46 336260 1 3653 3653 101792 101792 47 636765 1 1265 1265 210568 210568 48 481231 1 7489 7489 136996 136996 49 469107 0 4901 0 121920 0 50 211928 0 2284 0 76403 0 51 563925 1 3160 3160 108094 108094 52 511939 1 4150 4150 134759 134759 53 521016 1 7285 7285 188873 188873 54 543856 1 1134 1134 146216 146216 55 329304 1 4658 4658 156608 156608 56 423262 0 2384 0 61348 0 57 509665 0 3748 0 50350 0 58 455881 0 5371 0 87720 0 59 367772 0 1285 0 99489 0 60 406339 1 9327 9327 87419 87419 61 493408 1 5565 5565 94355 94355 62 232942 0 1528 0 60326 0 63 416002 1 3122 3122 94670 94670 64 337430 1 7317 7317 82425 82425 65 361517 0 2675 0 59017 0 66 360962 0 13253 0 90829 0 67 235561 0 880 0 80791 0 68 408247 1 2053 2053 100423 100423 69 450296 0 1424 0 131116 0 70 418799 1 4036 4036 100269 100269 71 247405 1 3045 3045 27330 27330 72 378519 0 5119 0 39039 0 73 326638 0 1431 0 106885 0 74 328233 0 554 0 79285 0 75 386225 0 1975 0 118881 0 76 283662 1 1286 1286 77623 77623 77 370225 0 1012 0 114768 0 78 269236 0 810 0 74015 0 79 365732 0 1280 0 69465 0 80 420383 1 666 666 117869 117869 81 345811 0 1380 0 60982 0 82 431809 1 4608 4608 90131 90131 83 418876 0 876 0 138971 0 84 297476 0 814 0 39625 0 85 416776 0 514 0 102725 0 86 357257 1 5692 5692 64239 64239 87 458343 0 3642 0 90262 0 88 388386 0 540 0 103960 0 89 358934 0 2099 0 106611 0 90 407560 0 567 0 103345 0 91 392558 0 2001 0 95551 0 92 373177 1 2949 2949 82903 82903 93 428370 0 2253 0 63593 0 94 369419 1 6533 6533 126910 126910 95 358649 0 1889 0 37527 0 96 376641 1 3055 3055 60247 60247 97 467427 0 272 0 112995 0 98 364885 1 1414 1414 70184 70184 99 436230 0 2564 0 130140 0 100 329118 1 1383 1383 73221 73221 101 317365 1 1261 1261 76114 76114 102 286849 0 975 0 90534 0 103 376685 0 3366 0 108479 0 104 407198 0 576 0 113761 0 105 377772 0 1306 0 68696 0 106 271483 0 746 0 71561 0 107 153661 1 3192 3192 59831 59831 108 513294 1 2045 2045 97890 97890 109 324881 0 5477 0 101481 0 110 264512 1 1932 1932 72954 72954 111 420968 0 936 0 67939 0 112 129302 1 3437 3437 48022 48022 113 191521 0 5131 0 86111 0 114 268673 1 2397 2397 74020 74020 115 353179 1 1389 1389 57530 57530 116 354624 0 1503 0 56364 0 117 363713 0 402 0 84990 0 118 456657 0 2239 0 88590 0 119 211742 1 2234 2234 77200 77200 120 338381 0 837 0 61262 0 121 418530 0 10579 0 110309 0 122 351483 0 875 0 67000 0 123 372928 0 1395 0 93099 0 124 485538 1 1659 1659 107577 107577 125 279268 1 2647 2647 62920 62920 126 219060 1 3294 3294 75832 75832 127 325560 0 0 0 60720 0 128 325314 0 94 0 60793 0 129 322046 0 422 0 57935 0 130 325560 0 0 0 60720 0 131 325599 0 34 0 60630 0 132 377028 0 1558 0 55637 0 133 325560 0 0 0 60720 0 134 323850 0 43 0 60887 0 135 325560 0 0 0 60720 0 136 331514 0 316 0 60505 0 137 325632 0 115 0 60945 0 138 325560 0 0 0 60720 0 139 325560 0 0 0 60720 0 140 325560 0 0 0 60720 0 141 322265 0 389 0 58990 0 142 325560 0 0 0 60720 0 143 325906 0 1002 0 56750 0 144 325985 0 36 0 60894 0 145 346145 0 460 0 63346 0 146 325898 0 309 0 56535 0 147 325560 0 0 0 60720 0 148 325356 0 9 0 60835 0 149 325560 0 0 0 60720 0 150 325930 0 14 0 61016 0 151 318020 0 520 0 58650 0 152 326389 0 1766 0 60438 0 153 325560 0 0 0 60720 0 154 302925 0 458 0 58625 0 155 325540 0 20 0 60938 0 156 325560 0 0 0 60720 0 157 325560 0 0 0 60720 0 158 326736 0 98 0 61490 0 159 340580 0 405 0 60845 0 160 325560 0 0 0 60720 0 161 325560 0 0 0 60720 0 162 325560 0 0 0 60720 0 163 325560 0 0 0 60720 0 164 331828 0 483 0 60830 0 165 323299 0 454 0 63261 0 166 325560 0 0 0 60720 0 167 325560 0 0 0 60720 0 168 387722 0 757 0 45689 0 169 325560 0 0 0 60720 0 170 325560 0 0 0 60720 0 171 325560 0 0 0 60720 0 172 324598 0 36 0 61564 0 173 325560 0 0 0 60720 0 174 328726 0 203 0 61938 0 175 325560 0 0 0 60720 0 176 325043 0 90 0 60951 0 177 325560 0 0 0 60720 0 178 325806 1 71 71 60745 60745 179 325560 0 0 0 60720 0 180 325560 0 0 0 60720 0 181 387732 0 972 0 71642 0 182 349729 1 531 531 71641 71641 183 332202 0 604 0 55792 0 184 305442 1 283 283 71873 71873 185 329537 1 23 23 62555 62555 186 327055 1 638 638 60370 60370 187 356245 1 699 699 64873 64873 188 328451 0 149 0 62041 0 189 307062 0 226 0 65745 0 190 325560 0 0 0 60720 0 191 331345 0 275 0 59500 0 192 325560 0 0 0 60720 0 193 331824 0 141 0 61630 0 194 325560 0 0 0 60720 0 195 325685 0 28 0 60890 0 196 325560 1 0 0 60720 60720 197 404480 1 2566 2566 113521 113521 198 325560 1 0 0 60720 60720 199 325560 1 0 0 60720 60720 200 318314 1 472 472 80045 80045 201 325560 1 0 0 60720 60720 202 325560 1 0 0 60720 60720 203 325560 1 0 0 60720 60720 204 311807 1 203 203 50804 50804 205 337724 1 496 496 87390 87390 206 326431 1 10 10 61656 61656 207 327556 1 63 63 65688 65688 208 325560 1 0 0 60720 60720 209 356850 1 1136 1136 48522 48522 210 325560 1 0 0 60720 60720 211 325560 0 0 0 60720 0 212 325560 0 0 0 60720 0 213 322741 0 267 0 57640 0 214 310902 0 474 0 61977 0 215 324295 0 534 0 62620 0 216 325560 0 0 0 60720 0 217 326156 0 15 0 60831 0 218 326960 0 397 0 60646 0 219 325560 0 0 0 60720 0 220 333411 0 1061 0 56225 0 221 297761 0 288 0 60510 0 222 325560 0 0 0 60720 0 223 325536 0 3 0 60698 0 224 325560 0 0 0 60720 0 225 325762 0 20 0 60805 0 226 327957 0 278 0 61404 0 227 325560 0 0 0 60720 0 228 325560 0 0 0 60720 0 229 318521 0 192 0 65276 0 230 325560 0 0 0 60720 0 231 319775 0 317 0 63915 0 232 325560 0 0 0 60720 0 233 325560 0 0 0 60720 0 234 332128 0 368 0 61686 0 235 325560 0 0 0 60720 0 236 325486 0 2 0 60743 0 237 325560 0 0 0 60720 0 238 325838 0 53 0 60349 0 239 325560 0 0 0 60720 0 240 325560 0 0 0 60720 0 241 325560 0 0 0 60720 0 242 331767 0 94 0 61360 0 243 325560 0 0 0 60720 0 244 324523 0 24 0 59818 0 245 339995 0 2332 0 72680 0 246 325560 1 0 0 60720 60720 247 325560 1 0 0 60720 60720 248 319582 0 131 0 61808 0 249 325560 1 0 0 60720 60720 250 325560 1 0 0 60720 60720 251 307245 0 206 0 53110 0 252 325560 1 0 0 60720 60720 253 317967 0 167 0 64245 0 254 331488 0 622 0 73007 0 255 335452 0 885 0 82732 0 256 325560 0 0 0 60720 0 257 334184 0 365 0 54820 0 258 313213 0 364 0 47705 0 259 325560 1 0 0 60720 60720 260 325560 0 0 0 60720 0 261 325560 1 0 0 60720 60720 262 325560 1 0 0 60720 60720 263 348678 0 226 0 72835 0 264 328727 0 307 0 58856 0 265 325560 1 0 0 60720 60720 266 325560 0 0 0 60720 0 267 325560 1 0 0 60720 60720 268 387978 0 188 0 77655 0 269 325560 1 0 0 60720 60720 270 336704 0 138 0 69817 0 271 325560 1 0 0 60720 60720 272 325560 1 0 0 60720 60720 273 325560 1 0 0 60720 60720 274 322076 0 125 0 60798 0 275 325560 0 0 0 60720 0 276 334272 0 282 0 62452 0 277 338197 0 335 0 64175 0 278 325560 0 0 0 60720 0 279 321024 1 813 813 67440 67440 280 322145 0 176 0 68136 0 281 325560 1 0 0 60720 60720 282 325560 0 0 0 60720 0 283 323351 0 249 0 56726 0 284 325560 0 0 0 60720 0 285 327748 0 333 0 70811 0 286 325560 0 0 0 60720 0 287 325560 1 0 0 60720 60720 288 328157 0 30 0 62045 0 289 325560 1 0 0 60720 60720 290 311594 0 249 0 54323 0 291 325560 1 0 0 60720 60720 292 335962 0 165 0 62841 0 293 372426 0 453 0 81125 0 294 325560 0 0 0 60720 0 295 319844 0 53 0 59506 0 296 355822 1 382 382 59365 59365 297 325560 0 0 0 60720 0 298 325560 1 0 0 60720 60720 299 325560 0 0 0 60720 0 300 325560 0 0 0 60720 0 301 324047 1 30 30 60798 60798 302 311464 0 290 0 58790 0 303 325560 0 0 0 60720 0 304 325560 1 0 0 60720 60720 305 353417 0 366 0 61808 0 306 325590 0 2 0 60735 0 307 325560 1 0 0 60720 60720 308 328576 1 209 209 64016 64016 309 326126 0 384 0 54683 0 310 325560 1 0 0 60720 60720 311 325560 1 0 0 60720 60720 312 369376 0 365 0 87192 0 313 325560 1 0 0 60720 60720 314 332013 1 49 49 64107 64107 315 325871 0 3 0 60761 0 316 342165 0 133 0 65990 0 317 324967 0 32 0 59988 0 318 314832 0 368 0 61167 0 319 325557 0 1 0 60719 0 320 325560 0 0 0 60720 0 321 325560 0 0 0 60720 0 322 325560 0 0 0 60720 0 323 325560 0 0 0 60720 0 324 325560 1 0 0 60720 60720 325 325560 0 0 0 60720 0 326 322649 0 22 0 60722 0 327 325560 1 0 0 60720 60720 328 325560 1 0 0 60720 60720 329 325560 0 0 0 60720 0 330 325560 0 0 0 60720 0 331 325560 0 0 0 60720 0 332 325560 0 0 0 60720 0 333 325560 0 0 0 60720 0 334 324598 0 96 0 60379 0 335 325567 0 1 0 60727 0 336 325560 0 0 0 60720 0 337 324005 0 81 0 60925 0 338 325560 0 0 0 60720 0 339 325748 0 26 0 60896 0 340 323385 0 125 0 59734 0 341 315409 0 304 0 62969 0 342 325560 0 0 0 60720 0 343 325560 0 0 0 60720 0 344 325560 0 0 0 60720 0 345 325560 0 0 0 60720 0 346 325560 0 0 0 60720 0 347 312275 0 119 0 59118 0 348 325560 0 0 0 60720 0 349 325560 0 0 0 60720 0 350 325560 0 0 0 60720 0 351 320576 0 312 0 58598 0 352 325246 0 60 0 61124 0 353 332961 0 587 0 59595 0 354 323010 0 135 0 62065 0 355 325560 0 0 0 60720 0 356 325560 0 0 0 60720 0 357 345253 0 514 0 78780 0 358 325560 0 0 0 60720 0 359 325560 0 0 0 60720 0 360 325560 0 0 0 60720 0 361 325559 0 1 0 60722 0 362 325560 0 0 0 60720 0 363 325560 0 0 0 60720 0 364 319634 1 58 58 61600 61600 365 319951 0 180 0 59635 0 366 325560 0 0 0 60720 0 367 325560 0 0 0 60720 0 368 325560 0 0 0 60720 0 369 325560 0 0 0 60720 0 370 325560 0 0 0 60720 0 371 325560 0 0 0 60720 0 372 318519 0 448 0 59781 0 373 343222 0 227 0 76644 0 374 317234 0 174 0 64820 0 375 325560 0 0 0 60720 0 376 325560 0 0 0 60720 0 377 314025 0 121 0 56178 0 378 320249 0 607 0 60436 0 379 325560 0 0 0 60720 0 380 325560 0 0 0 60720 0 381 325560 0 0 0 60720 0 382 349365 0 530 0 73433 0 383 289197 0 571 0 41477 0 384 325560 0 0 0 60720 0 385 329245 0 78 0 62700 0 386 240869 0 2489 0 67804 0 387 327182 0 131 0 59661 0 388 322876 0 923 0 58620 0 389 323117 0 72 0 60398 0 390 306351 0 572 0 58580 0 391 335137 0 397 0 62710 0 392 308271 0 450 0 59325 0 393 301731 0 622 0 60950 0 394 382409 0 694 0 68060 0 395 279230 1 3425 3425 83620 83620 396 298731 0 562 0 58456 0 397 243650 0 4917 0 52811 0 398 532682 1 1442 1442 121173 121173 399 319771 0 529 0 63870 0 400 171493 1 2126 2126 21001 21001 401 347262 0 1061 0 70415 0 402 343945 0 776 0 64230 0 403 311874 0 611 0 59190 0 404 302211 1 1526 1526 69351 69351 405 316708 0 592 0 64270 0 406 333463 0 1182 0 70694 0 407 344282 0 621 0 68005 0 408 319635 0 989 0 58930 0 409 301186 0 438 0 58320 0 410 300381 0 726 0 69980 0 411 318765 0 1303 0 69863 0 412 286146 1 6341 6341 63255 63255 413 306844 1 1164 1164 57320 57320 414 307705 1 3310 3310 75230 75230 415 312448 0 1366 0 79420 0 416 299715 0 965 0 73490 0 417 373399 0 3256 0 35250 0 418 299446 1 1135 1135 62285 62285 419 325586 0 1270 0 69206 0 420 291221 0 661 0 65920 0 421 261173 0 1013 0 69770 0 422 255027 0 2844 0 72683 0 423 -78375 1 11528 11528 -14545 -14545 424 -58143 0 6526 0 55830 0 425 227033 0 2264 0 55174 0 426 235098 1 4461 4461 67038 67038 427 21267 0 3999 0 51252 0 428 238675 0 35624 0 157278 0 429 197687 0 9252 0 79510 0 430 418341 0 15236 0 77440 0 431 -297706 0 18073 0 27284 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs GrCosts Dividends GrDiv 117378.943 45851.725 -7.650 43.285 3.424 -1.613 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -855582.7 -18397.8 269.7 25822.7 2951183.2 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.174e+05 3.105e+04 3.780 0.000179 *** Group 4.585e+04 5.561e+04 0.825 0.410105 Costs -7.650e+00 1.721e+00 -4.445 1.12e-05 *** GrCosts 4.328e+01 2.065e+00 20.965 < 2e-16 *** Dividends 3.424e+00 4.168e-01 8.215 2.59e-15 *** GrDiv -1.613e+00 7.000e-01 -2.305 0.021671 * --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 210000 on 425 degrees of freedom Multiple R-squared: 0.7822, Adjusted R-squared: 0.7796 F-statistic: 305.2 on 5 and 425 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 3.928338e-86 1.964169e-86 [2,] 1.0000000 7.940290e-92 3.970145e-92 [3,] 1.0000000 1.742086e-94 8.710430e-95 [4,] 1.0000000 1.382743e-105 6.913717e-106 [5,] 1.0000000 2.448957e-120 1.224478e-120 [6,] 1.0000000 3.577209e-129 1.788604e-129 [7,] 1.0000000 1.426528e-132 7.132642e-133 [8,] 1.0000000 2.813127e-143 1.406563e-143 [9,] 1.0000000 8.847252e-148 4.423626e-148 [10,] 1.0000000 1.496397e-149 7.481987e-150 [11,] 1.0000000 3.046286e-155 1.523143e-155 [12,] 1.0000000 6.495943e-155 3.247971e-155 [13,] 1.0000000 4.108295e-159 2.054147e-159 [14,] 1.0000000 1.390631e-161 6.953153e-162 [15,] 1.0000000 1.047230e-162 5.236151e-163 [16,] 1.0000000 6.338407e-163 3.169204e-163 [17,] 1.0000000 1.270486e-162 6.352429e-163 [18,] 1.0000000 1.959619e-164 9.798096e-165 [19,] 1.0000000 1.173600e-167 5.867998e-168 [20,] 1.0000000 9.010655e-170 4.505328e-170 [21,] 1.0000000 1.834766e-182 9.173829e-183 [22,] 1.0000000 1.405182e-184 7.025911e-185 [23,] 1.0000000 3.359450e-188 1.679725e-188 [24,] 1.0000000 1.034151e-189 5.170757e-190 [25,] 1.0000000 5.416736e-196 2.708368e-196 [26,] 1.0000000 2.345896e-197 1.172948e-197 [27,] 1.0000000 1.130549e-197 5.652746e-198 [28,] 1.0000000 2.421284e-197 1.210642e-197 [29,] 1.0000000 1.792957e-204 8.964787e-205 [30,] 1.0000000 2.213544e-209 1.106772e-209 [31,] 1.0000000 1.276771e-208 6.383857e-209 [32,] 1.0000000 4.138455e-209 2.069227e-209 [33,] 1.0000000 1.025320e-208 5.126599e-209 [34,] 1.0000000 1.991331e-208 9.956653e-209 [35,] 1.0000000 4.289176e-209 2.144588e-209 [36,] 1.0000000 2.473127e-208 1.236564e-208 [37,] 1.0000000 7.943502e-208 3.971751e-208 [38,] 1.0000000 3.095188e-207 1.547594e-207 [39,] 1.0000000 2.249467e-206 1.124734e-206 [40,] 1.0000000 9.188083e-206 4.594042e-206 [41,] 1.0000000 3.330634e-205 1.665317e-205 [42,] 1.0000000 9.456548e-206 4.728274e-206 [43,] 1.0000000 1.357956e-206 6.789779e-207 [44,] 1.0000000 8.643769e-206 4.321884e-206 [45,] 1.0000000 3.851816e-205 1.925908e-205 [46,] 1.0000000 4.092463e-204 2.046231e-204 [47,] 1.0000000 8.495437e-207 4.247718e-207 [48,] 1.0000000 6.144824e-207 3.072412e-207 [49,] 1.0000000 4.076781e-211 2.038390e-211 [50,] 1.0000000 1.229191e-211 6.145953e-212 [51,] 1.0000000 7.842217e-211 3.921108e-211 [52,] 1.0000000 2.662363e-211 1.331182e-211 [53,] 1.0000000 2.672890e-212 1.336445e-212 [54,] 1.0000000 5.689039e-212 2.844519e-212 [55,] 1.0000000 4.344623e-211 2.172312e-211 [56,] 1.0000000 1.043557e-210 5.217783e-211 [57,] 1.0000000 4.103002e-210 2.051501e-210 [58,] 1.0000000 1.922985e-210 9.614924e-211 [59,] 1.0000000 3.735838e-211 1.867919e-211 [60,] 1.0000000 4.247557e-210 2.123779e-210 [61,] 1.0000000 3.216746e-209 1.608373e-209 [62,] 1.0000000 1.922920e-208 9.614602e-209 [63,] 1.0000000 7.495868e-208 3.747934e-208 [64,] 1.0000000 1.919947e-209 9.599734e-210 [65,] 1.0000000 2.192575e-209 1.096288e-209 [66,] 1.0000000 1.566899e-208 7.834494e-209 [67,] 1.0000000 6.523269e-208 3.261634e-208 [68,] 1.0000000 2.368979e-207 1.184489e-207 [69,] 1.0000000 5.077727e-207 2.538863e-207 [70,] 1.0000000 9.815383e-207 4.907692e-207 [71,] 1.0000000 7.871039e-206 3.935519e-206 [72,] 1.0000000 4.154043e-205 2.077021e-205 [73,] 1.0000000 3.180501e-204 1.590250e-204 [74,] 1.0000000 2.828822e-204 1.414411e-204 [75,] 1.0000000 2.172224e-204 1.086112e-204 [76,] 1.0000000 1.601655e-203 8.008275e-204 [77,] 1.0000000 1.502963e-202 7.514817e-203 [78,] 1.0000000 2.546787e-203 1.273393e-203 [79,] 1.0000000 1.845281e-203 9.226405e-204 [80,] 1.0000000 1.051185e-202 5.255923e-203 [81,] 1.0000000 3.332448e-202 1.666224e-202 [82,] 1.0000000 2.673888e-201 1.336944e-201 [83,] 1.0000000 2.672146e-200 1.336073e-200 [84,] 1.0000000 1.675453e-199 8.377263e-200 [85,] 1.0000000 2.694082e-200 1.347041e-200 [86,] 1.0000000 1.583800e-199 7.919001e-200 [87,] 1.0000000 2.222162e-200 1.111081e-200 [88,] 1.0000000 7.114298e-201 3.557149e-201 [89,] 1.0000000 7.086767e-200 3.543384e-200 [90,] 1.0000000 4.437362e-199 2.218681e-199 [91,] 1.0000000 1.426962e-198 7.134811e-199 [92,] 1.0000000 1.470267e-197 7.351335e-198 [93,] 1.0000000 1.279797e-196 6.398987e-197 [94,] 1.0000000 2.264937e-197 1.132468e-197 [95,] 1.0000000 8.535786e-197 4.267893e-197 [96,] 1.0000000 7.159051e-197 3.579525e-197 [97,] 1.0000000 4.177779e-196 2.088889e-196 [98,] 1.0000000 5.938464e-196 2.969232e-196 [99,] 1.0000000 1.409756e-196 7.048780e-197 [100,] 1.0000000 8.485378e-198 4.242689e-198 [101,] 1.0000000 1.867575e-197 9.337874e-198 [102,] 1.0000000 7.653100e-197 3.826550e-197 [103,] 1.0000000 8.894456e-197 4.447228e-197 [104,] 1.0000000 2.892081e-197 1.446041e-197 [105,] 1.0000000 1.085784e-199 5.428922e-200 [106,] 1.0000000 5.269902e-199 2.634951e-199 [107,] 1.0000000 1.544519e-198 7.722594e-199 [108,] 1.0000000 3.206145e-198 1.603073e-198 [109,] 1.0000000 1.426600e-197 7.133002e-198 [110,] 1.0000000 5.069907e-197 2.534953e-197 [111,] 1.0000000 4.254588e-198 2.127294e-198 [112,] 1.0000000 3.629878e-197 1.814939e-197 [113,] 1.0000000 9.976852e-197 4.988426e-197 [114,] 1.0000000 9.432640e-196 4.716320e-196 [115,] 1.0000000 3.551342e-195 1.775671e-195 [116,] 1.0000000 6.585015e-195 3.292508e-195 [117,] 1.0000000 6.502460e-194 3.251230e-194 [118,] 1.0000000 3.583924e-194 1.791962e-194 [119,] 1.0000000 3.708065e-193 1.854033e-193 [120,] 1.0000000 3.810155e-192 1.905078e-192 [121,] 1.0000000 3.599949e-191 1.799975e-191 [122,] 1.0000000 3.663133e-190 1.831567e-190 [123,] 1.0000000 3.702144e-189 1.851072e-189 [124,] 1.0000000 1.329667e-189 6.648335e-190 [125,] 1.0000000 1.344765e-188 6.723827e-189 [126,] 1.0000000 1.353959e-187 6.769795e-188 [127,] 1.0000000 1.354422e-186 6.772108e-187 [128,] 1.0000000 1.264228e-185 6.321138e-186 [129,] 1.0000000 1.248822e-184 6.244110e-185 [130,] 1.0000000 1.229087e-183 6.145434e-184 [131,] 1.0000000 1.203124e-182 6.015622e-183 [132,] 1.0000000 1.171361e-181 5.856807e-182 [133,] 1.0000000 1.084620e-180 5.423102e-181 [134,] 1.0000000 1.044744e-179 5.223721e-180 [135,] 1.0000000 6.805478e-179 3.402739e-179 [136,] 1.0000000 6.489104e-178 3.244552e-178 [137,] 1.0000000 5.404755e-177 2.702378e-177 [138,] 1.0000000 4.191970e-176 2.095985e-176 [139,] 1.0000000 3.938757e-175 1.969378e-175 [140,] 1.0000000 3.684201e-174 1.842100e-174 [141,] 1.0000000 3.425382e-173 1.712691e-173 [142,] 1.0000000 3.170865e-172 1.585432e-172 [143,] 1.0000000 2.802378e-171 1.401189e-171 [144,] 1.0000000 1.838953e-170 9.194767e-171 [145,] 1.0000000 1.676088e-169 8.380441e-170 [146,] 1.0000000 1.508932e-168 7.544659e-169 [147,] 1.0000000 1.362555e-167 6.812776e-168 [148,] 1.0000000 1.222917e-166 6.114583e-167 [149,] 1.0000000 1.091946e-165 5.459731e-166 [150,] 1.0000000 9.705648e-165 4.852824e-165 [151,] 1.0000000 7.157104e-164 3.578552e-164 [152,] 1.0000000 6.295811e-163 3.147906e-163 [153,] 1.0000000 5.509979e-162 2.754989e-162 [154,] 1.0000000 4.797745e-161 2.398872e-161 [155,] 1.0000000 4.156410e-160 2.078205e-160 [156,] 1.0000000 3.298731e-159 1.649366e-159 [157,] 1.0000000 2.830485e-158 1.415242e-158 [158,] 1.0000000 2.415774e-157 1.207887e-157 [159,] 1.0000000 2.051471e-156 1.025735e-156 [160,] 1.0000000 3.197620e-158 1.598810e-158 [161,] 1.0000000 2.750689e-157 1.375344e-157 [162,] 1.0000000 2.354424e-156 1.177212e-156 [163,] 1.0000000 2.005213e-155 1.002606e-155 [164,] 1.0000000 1.707968e-154 8.539840e-155 [165,] 1.0000000 1.440299e-153 7.201497e-154 [166,] 1.0000000 1.203869e-152 6.019347e-153 [167,] 1.0000000 1.005179e-151 5.025894e-152 [168,] 1.0000000 8.352321e-151 4.176161e-151 [169,] 1.0000000 6.905243e-150 3.452622e-150 [170,] 1.0000000 5.522659e-149 2.761330e-149 [171,] 1.0000000 4.521798e-148 2.260899e-148 [172,] 1.0000000 3.684172e-147 1.842086e-147 [173,] 1.0000000 1.969256e-146 9.846281e-147 [174,] 1.0000000 1.554459e-145 7.772297e-146 [175,] 1.0000000 7.398860e-145 3.699430e-145 [176,] 1.0000000 3.163397e-144 1.581699e-144 [177,] 1.0000000 2.468818e-143 1.234409e-143 [178,] 1.0000000 1.858128e-142 9.290638e-143 [179,] 1.0000000 1.061147e-141 5.305733e-142 [180,] 1.0000000 8.384353e-141 4.192176e-141 [181,] 1.0000000 4.849314e-140 2.424657e-140 [182,] 1.0000000 3.799232e-139 1.899616e-139 [183,] 1.0000000 2.642958e-138 1.321479e-138 [184,] 1.0000000 2.051161e-137 1.025580e-137 [185,] 1.0000000 1.552772e-136 7.763860e-137 [186,] 1.0000000 1.193501e-135 5.967507e-136 [187,] 1.0000000 9.133105e-135 4.566552e-135 [188,] 1.0000000 6.828279e-134 3.414139e-134 [189,] 1.0000000 4.661675e-133 2.330837e-133 [190,] 1.0000000 3.458445e-132 1.729222e-132 [191,] 1.0000000 2.557236e-131 1.278618e-131 [192,] 1.0000000 8.332686e-131 4.166343e-131 [193,] 1.0000000 6.135147e-130 3.067573e-130 [194,] 1.0000000 4.501097e-129 2.250548e-129 [195,] 1.0000000 3.290254e-128 1.645127e-128 [196,] 1.0000000 2.211593e-127 1.105796e-127 [197,] 1.0000000 6.135759e-127 3.067879e-127 [198,] 1.0000000 4.446292e-126 2.223146e-126 [199,] 1.0000000 3.073861e-125 1.536930e-125 [200,] 1.0000000 2.212988e-124 1.106494e-124 [201,] 1.0000000 1.851559e-124 9.257793e-125 [202,] 1.0000000 1.336855e-123 6.684274e-124 [203,] 1.0000000 9.674995e-123 4.837498e-123 [204,] 1.0000000 6.968608e-122 3.484304e-122 [205,] 1.0000000 4.523724e-121 2.261862e-121 [206,] 1.0000000 3.189903e-120 1.594951e-120 [207,] 1.0000000 2.270371e-119 1.135186e-119 [208,] 1.0000000 1.605219e-118 8.026096e-119 [209,] 1.0000000 1.128320e-117 5.641599e-118 [210,] 1.0000000 7.588995e-117 3.794497e-117 [211,] 1.0000000 5.289653e-116 2.644827e-116 [212,] 1.0000000 1.899713e-115 9.498566e-116 [213,] 1.0000000 1.215582e-114 6.077910e-115 [214,] 1.0000000 8.379833e-114 4.189917e-114 [215,] 1.0000000 5.747219e-113 2.873609e-113 [216,] 1.0000000 3.924237e-112 1.962119e-112 [217,] 1.0000000 2.665412e-111 1.332706e-111 [218,] 1.0000000 1.774317e-110 8.871586e-111 [219,] 1.0000000 1.194242e-109 5.971211e-110 [220,] 1.0000000 7.999663e-109 3.999831e-109 [221,] 1.0000000 4.930728e-108 2.465364e-108 [222,] 1.0000000 3.274238e-107 1.637119e-107 [223,] 1.0000000 2.134144e-106 1.067072e-106 [224,] 1.0000000 1.404143e-105 7.020714e-106 [225,] 1.0000000 9.194109e-105 4.597054e-105 [226,] 1.0000000 5.741959e-104 2.870980e-104 [227,] 1.0000000 3.723910e-103 1.861955e-103 [228,] 1.0000000 2.404482e-102 1.202241e-102 [229,] 1.0000000 1.544430e-101 7.722151e-102 [230,] 1.0000000 9.775098e-101 4.887549e-101 [231,] 1.0000000 6.218180e-100 3.109090e-100 [232,] 1.0000000 3.936413e-99 1.968207e-99 [233,] 1.0000000 2.479880e-98 1.239940e-98 [234,] 1.0000000 1.521671e-97 7.608353e-98 [235,] 1.0000000 9.493580e-97 4.746790e-97 [236,] 1.0000000 5.823001e-96 2.911501e-96 [237,] 1.0000000 3.590292e-95 1.795146e-95 [238,] 1.0000000 2.203973e-94 1.101986e-94 [239,] 1.0000000 1.347161e-93 6.735807e-94 [240,] 1.0000000 8.268122e-93 4.134061e-93 [241,] 1.0000000 5.007893e-92 2.503947e-92 [242,] 1.0000000 3.019966e-91 1.509983e-91 [243,] 1.0000000 1.606507e-90 8.032535e-91 [244,] 1.0000000 9.606200e-90 4.803100e-90 [245,] 1.0000000 5.513644e-89 2.756822e-89 [246,] 1.0000000 2.641251e-88 1.320626e-88 [247,] 1.0000000 5.664043e-88 2.832021e-88 [248,] 1.0000000 3.353232e-87 1.676616e-87 [249,] 1.0000000 1.199015e-86 5.995076e-87 [250,] 1.0000000 3.241162e-86 1.620581e-86 [251,] 1.0000000 1.904104e-85 9.520519e-86 [252,] 1.0000000 1.111337e-84 5.556686e-85 [253,] 1.0000000 6.467634e-84 3.233817e-84 [254,] 1.0000000 3.746668e-83 1.873334e-83 [255,] 1.0000000 1.965936e-82 9.829681e-83 [256,] 1.0000000 1.006659e-81 5.033297e-82 [257,] 1.0000000 5.755550e-81 2.877775e-81 [258,] 1.0000000 3.265676e-80 1.632838e-80 [259,] 1.0000000 1.849259e-79 9.246295e-80 [260,] 1.0000000 1.046569e-78 5.232845e-79 [261,] 1.0000000 5.868955e-78 2.934478e-78 [262,] 1.0000000 2.986678e-77 1.493339e-77 [263,] 1.0000000 1.659591e-76 8.297957e-77 [264,] 1.0000000 9.177319e-76 4.588660e-76 [265,] 1.0000000 5.050318e-75 2.525159e-75 [266,] 1.0000000 2.768578e-74 1.384289e-74 [267,] 1.0000000 1.503327e-73 7.516636e-74 [268,] 1.0000000 7.947264e-73 3.973632e-73 [269,] 1.0000000 4.236644e-72 2.118322e-72 [270,] 1.0000000 2.264846e-71 1.132423e-71 [271,] 1.0000000 1.182934e-70 5.914672e-71 [272,] 1.0000000 5.398892e-70 2.699446e-70 [273,] 1.0000000 2.857594e-69 1.428797e-69 [274,] 1.0000000 1.499261e-68 7.496307e-69 [275,] 1.0000000 6.828905e-68 3.414452e-68 [276,] 1.0000000 3.546673e-67 1.773337e-67 [277,] 1.0000000 1.510713e-66 7.553564e-67 [278,] 1.0000000 7.780897e-66 3.890448e-66 [279,] 1.0000000 4.002323e-65 2.001161e-65 [280,] 1.0000000 2.052412e-64 1.026206e-64 [281,] 1.0000000 1.044705e-63 5.223523e-64 [282,] 1.0000000 4.649705e-63 2.324853e-63 [283,] 1.0000000 2.344256e-62 1.172128e-62 [284,] 1.0000000 1.156762e-61 5.783810e-62 [285,] 1.0000000 4.717028e-61 2.358514e-61 [286,] 1.0000000 2.333805e-60 1.166903e-60 [287,] 1.0000000 1.145721e-59 5.728604e-60 [288,] 1.0000000 4.001051e-59 2.000526e-59 [289,] 1.0000000 1.952411e-58 9.762054e-59 [290,] 1.0000000 9.501120e-58 4.750560e-58 [291,] 1.0000000 4.585141e-57 2.292571e-57 [292,] 1.0000000 2.200239e-56 1.100119e-56 [293,] 1.0000000 1.055220e-55 5.276101e-56 [294,] 1.0000000 5.011137e-55 2.505569e-55 [295,] 1.0000000 2.364138e-54 1.182069e-54 [296,] 1.0000000 1.112308e-53 5.561539e-54 [297,] 1.0000000 4.162377e-53 2.081189e-53 [298,] 1.0000000 1.932435e-52 9.662176e-53 [299,] 1.0000000 8.952636e-52 4.476318e-52 [300,] 1.0000000 4.144818e-51 2.072409e-51 [301,] 1.0000000 1.399651e-50 6.998256e-51 [302,] 1.0000000 6.387042e-50 3.193521e-50 [303,] 1.0000000 2.895279e-49 1.447640e-49 [304,] 1.0000000 6.045406e-49 3.022703e-49 [305,] 1.0000000 2.724850e-48 1.362425e-48 [306,] 1.0000000 1.226473e-47 6.132366e-48 [307,] 1.0000000 5.458510e-47 2.729255e-47 [308,] 1.0000000 2.438763e-46 1.219382e-46 [309,] 1.0000000 1.061359e-45 5.306793e-46 [310,] 1.0000000 4.681957e-45 2.340978e-45 [311,] 1.0000000 2.033932e-44 1.016966e-44 [312,] 1.0000000 8.780858e-44 4.390429e-44 [313,] 1.0000000 3.766972e-43 1.883486e-43 [314,] 1.0000000 1.605777e-42 8.028884e-43 [315,] 1.0000000 6.801385e-42 3.400693e-42 [316,] 1.0000000 2.860616e-41 1.430308e-41 [317,] 1.0000000 1.196369e-40 5.981843e-41 [318,] 1.0000000 4.997265e-40 2.498632e-40 [319,] 1.0000000 2.055771e-39 1.027886e-39 [320,] 1.0000000 8.355660e-39 4.177830e-39 [321,] 1.0000000 3.405775e-38 1.702888e-38 [322,] 1.0000000 1.378918e-37 6.894590e-38 [323,] 1.0000000 5.545323e-37 2.772661e-37 [324,] 1.0000000 2.214926e-36 1.107463e-36 [325,] 1.0000000 8.786449e-36 4.393224e-36 [326,] 1.0000000 3.440225e-35 1.720112e-35 [327,] 1.0000000 1.345946e-34 6.729730e-35 [328,] 1.0000000 5.228877e-34 2.614439e-34 [329,] 1.0000000 2.022646e-33 1.011323e-33 [330,] 1.0000000 7.746975e-33 3.873487e-33 [331,] 1.0000000 2.947632e-32 1.473816e-32 [332,] 1.0000000 1.097230e-31 5.486151e-32 [333,] 1.0000000 4.140778e-31 2.070389e-31 [334,] 1.0000000 1.541156e-30 7.705778e-31 [335,] 1.0000000 5.693466e-30 2.846733e-30 [336,] 1.0000000 2.087580e-29 1.043790e-29 [337,] 1.0000000 7.596533e-29 3.798266e-29 [338,] 1.0000000 2.743233e-28 1.371616e-28 [339,] 1.0000000 9.899384e-28 4.949692e-28 [340,] 1.0000000 3.519871e-27 1.759936e-27 [341,] 1.0000000 1.241716e-26 6.208580e-27 [342,] 1.0000000 4.345703e-26 2.172851e-26 [343,] 1.0000000 1.453869e-25 7.269346e-26 [344,] 1.0000000 5.022525e-25 2.511262e-25 [345,] 1.0000000 1.539108e-24 7.695542e-25 [346,] 1.0000000 5.284325e-24 2.642162e-24 [347,] 1.0000000 1.776639e-23 8.883194e-24 [348,] 1.0000000 5.922694e-23 2.961347e-23 [349,] 1.0000000 1.580674e-22 7.903370e-23 [350,] 1.0000000 5.210726e-22 2.605363e-22 [351,] 1.0000000 1.702699e-21 8.513497e-22 [352,] 1.0000000 5.514615e-21 2.757307e-21 [353,] 1.0000000 1.770012e-20 8.850058e-21 [354,] 1.0000000 5.629662e-20 2.814831e-20 [355,] 1.0000000 1.774095e-19 8.870477e-20 [356,] 1.0000000 5.614113e-19 2.807057e-19 [357,] 1.0000000 1.724782e-18 8.623908e-19 [358,] 1.0000000 5.285323e-18 2.642661e-18 [359,] 1.0000000 1.603918e-17 8.019591e-18 [360,] 1.0000000 4.819542e-17 2.409771e-17 [361,] 1.0000000 1.433770e-16 7.168849e-17 [362,] 1.0000000 4.222194e-16 2.111097e-16 [363,] 1.0000000 1.230593e-15 6.152963e-16 [364,] 1.0000000 3.501694e-15 1.750847e-15 [365,] 1.0000000 8.914665e-15 4.457332e-15 [366,] 1.0000000 2.551357e-14 1.275678e-14 [367,] 1.0000000 7.175477e-14 3.587739e-14 [368,] 1.0000000 1.995645e-13 9.978223e-14 [369,] 1.0000000 5.263768e-13 2.631884e-13 [370,] 1.0000000 1.408639e-12 7.043197e-13 [371,] 1.0000000 3.785000e-12 1.892500e-12 [372,] 1.0000000 1.004888e-11 5.024440e-12 [373,] 1.0000000 2.635440e-11 1.317720e-11 [374,] 1.0000000 6.900571e-11 3.450285e-11 [375,] 1.0000000 1.090952e-10 5.454762e-11 [376,] 1.0000000 2.781316e-10 1.390658e-10 [377,] 1.0000000 7.083336e-10 3.541668e-10 [378,] 1.0000000 1.463665e-09 7.318327e-10 [379,] 1.0000000 3.528912e-09 1.764456e-09 [380,] 1.0000000 7.904060e-09 3.952030e-09 [381,] 1.0000000 1.900859e-08 9.504296e-09 [382,] 1.0000000 4.492066e-08 2.246033e-08 [383,] 0.9999999 1.033325e-07 5.166626e-08 [384,] 0.9999999 2.391132e-07 1.195566e-07 [385,] 0.9999997 5.564545e-07 2.782272e-07 [386,] 0.9999994 1.103763e-06 5.518813e-07 [387,] 0.9999990 2.022141e-06 1.011070e-06 [388,] 0.9999978 4.460880e-06 2.230440e-06 [389,] 0.9999960 8.098457e-06 4.049229e-06 [390,] 0.9999916 1.682333e-05 8.411663e-06 [391,] 0.9999820 3.601403e-05 1.800701e-05 [392,] 0.9999622 7.554771e-05 3.777386e-05 [393,] 0.9999222 1.556781e-04 7.783906e-05 [394,] 0.9998525 2.949126e-04 1.474563e-04 [395,] 0.9997201 5.597361e-04 2.798680e-04 [396,] 0.9994492 1.101613e-03 5.508066e-04 [397,] 0.9989476 2.104760e-03 1.052380e-03 [398,] 0.9980163 3.967341e-03 1.983671e-03 [399,] 0.9964360 7.128077e-03 3.564038e-03 [400,] 0.9942112 1.157757e-02 5.788786e-03 [401,] 0.9902445 1.951092e-02 9.755458e-03 [402,] 0.9830383 3.392334e-02 1.696167e-02 [403,] 0.9715886 5.682286e-02 2.841143e-02 [404,] 0.9566296 8.674078e-02 4.337039e-02 [405,] 0.9306208 1.387585e-01 6.937924e-02 [406,] 0.8933875 2.132250e-01 1.066125e-01 [407,] 0.8425038 3.149924e-01 1.574962e-01 [408,] 0.7738484 4.523031e-01 2.261516e-01 [409,] 0.9108151 1.783697e-01 8.918486e-02 [410,] 0.8508713 2.982575e-01 1.491287e-01 [411,] 0.7709301 4.581398e-01 2.290699e-01 [412,] 0.6587239 6.825523e-01 3.412761e-01 [413,] 0.5139631 9.720739e-01 4.860369e-01 [414,] 0.3521063 7.042125e-01 6.478937e-01 > postscript(file="/var/www/html/rcomp/tmp/1jzi41291335824.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/2jzi41291335824.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/3jzi41291335824.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/4crz71291335824.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/5crz71291335824.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 -58916.68341 2951183.24003 547776.97946 -592723.57478 -855582.68318 6 7 8 9 10 -187303.82930 1120158.08279 422848.11045 -703591.95988 613805.40397 11 12 13 14 15 459456.04993 800071.51561 655359.62146 -468897.87437 -280661.11852 16 17 18 19 20 565520.53288 -259523.08584 -337447.59549 -581028.38399 52511.45770 21 22 23 24 25 101447.04410 -391496.71174 120088.18252 -127599.19151 -71606.42087 26 27 28 29 30 326145.45395 -639817.67921 -174375.28759 621327.36936 289543.00246 31 32 33 34 35 247663.51192 -168131.35366 151317.62041 88593.24434 -311615.01065 36 37 38 39 40 52219.10887 310262.12965 140406.31363 105590.98411 -133523.16907 41 42 43 44 45 -148996.18697 30995.22637 97726.72451 -53490.33926 61041.57451 46 47 48 49 50 -141469.50111 47161.91860 -196941.82015 -28247.02662 -149590.54896 51 52 53 54 55 92352.03082 -43197.40228 -243825.49661 75448.94087 -283498.95034 56 57 58 59 60 114058.25801 248553.73190 79226.58274 -80437.40445 -247557.75152 61 62 63 64 65 -38989.03249 -79310.46952 -29908.81807 -235797.07384 62520.90448 66 67 68 69 70 33957.12914 -151722.70630 -9987.33107 -105144.12475 -69820.92442 71 72 73 74 75 -73823.45845 166625.49732 -145779.19403 -56387.81872 -123106.25890 76 77 78 79 80 -65954.22615 -132389.61114 -95381.48058 20289.54469 19983.58345 81 82 83 84 85 30180.16223 -58836.44905 -167652.48064 50643.93943 -48411.72652 86 87 88 89 90 -125132.00610 59758.15692 -80831.59816 -107434.97899 -59345.23388 91 92 93 94 95 -36690.09821 -45261.40047 110477.03904 -256423.03336 127224.14916 96 97 98 99 100 -4549.51391 -34777.47183 24177.92868 -107147.52429 -15983.74517 101 102 103 104 105 -28627.86911 -133068.99569 -86388.00661 -95303.81873 35161.57075 106 107 108 109 110 -85221.31954 -231658.23841 99931.47560 -98081.59614 -99669.94893 111 112 113 114 115 78119.21937 -243364.25034 -181459.97158 -114009.58235 36276.51649 116 117 118 119 120 55746.57428 -41605.06999 53064.69920 -170890.36224 17637.61731 121 122 123 124 125 4368.30333 11382.81709 -52559.93239 68389.53214 -92223.61851 126 127 128 129 130 -198868.43790 269.66610 492.78022 9519.96729 269.66610 131 132 133 134 135 876.92556 81060.62993 269.66610 -1683.22084 269.66610 136 137 138 139 140 9377.15855 450.96100 269.66610 269.66610 269.66610 141 142 143 144 145 5874.10042 269.66610 21874.36284 374.26243 15381.84798 146 147 148 149 150 17301.29008 269.66610 -259.25801 269.66610 -266.77164 151 152 153 154 155 3795.40766 15573.66403 269.66610 -11688.27232 -343.79348 156 157 158 159 160 269.66610 269.66610 -441.21905 17959.78921 269.66610 161 162 163 164 165 269.66610 269.66610 269.66610 9855.82870 -7219.00178 166 167 168 169 170 269.66610 269.66610 119690.15922 269.66610 269.66610 171 172 173 174 175 269.66610 -3306.88500 269.66610 818.00422 269.66610 176 177 178 179 180 -349.82655 269.66610 50048.92944 269.66610 269.66610 181 182 183 184 185 32479.16452 37849.42693 28406.06360 1979.83537 52212.88945 186 187 188 189 190 31771.85850 50634.13806 -222.76290 -33705.60342 269.66610 191 192 193 194 195 12335.74127 269.66610 4496.34472 269.66610 26.76109 196 197 198 199 200 52378.29308 -55752.85014 52378.29308 52378.29308 -6680.96187 201 202 203 204 205 52378.29308 52378.29308 52378.29308 49347.13500 -1426.44177 206 207 208 209 210 51198.04438 43133.28603 52378.29308 65274.78044 52378.29308 211 212 213 214 215 269.66610 269.66610 10039.37072 -15066.46219 -3416.17584 216 217 218 219 220 269.66610 600.33670 4959.98750 269.66610 31628.34903 221 222 223 224 225 -24607.15405 269.66610 343.94547 269.66610 333.61190 226 227 228 229 230 2451.20280 269.66610 269.66610 -20900.79062 269.66610 231 232 233 234 235 -14030.37493 269.66610 269.66610 6345.08101 269.66610 236 237 238 239 240 132.21122 269.66610 2223.44248 269.66610 269.66610 241 242 243 244 245 269.66610 5004.31516 269.66610 2504.79811 -8408.43581 246 247 248 249 250 52378.29308 52378.29308 -8431.64237 52378.29308 52378.29308 251 252 253 254 255 9587.91386 52378.29308 -18115.78580 -31116.13407 -58439.63706 256 257 258 259 260 269.66610 31888.00688 35271.83324 52378.29308 269.66610 261 262 263 264 265 52378.29308 52378.29308 -16366.47699 12167.65305 52378.29308 266 267 268 269 270 269.66610 52378.29308 6138.66863 52378.29308 -18679.71641 271 272 273 274 275 52378.29308 52378.29308 52378.29308 -2525.19903 269.66610 276 277 278 279 280 5208.34419 3639.05382 269.66610 6702.44125 -27192.11395 281 282 283 284 285 52378.29308 269.66610 13641.30379 269.66610 -29547.57748 286 287 288 289 290 269.66610 52378.29308 -1440.77572 52378.29308 10112.41764 291 292 293 294 295 52378.29308 4671.35207 -19267.78540 269.66610 -884.04067 296 297 298 299 300 71481.29254 269.66610 52378.29308 269.66610 269.66610 301 302 303 304 305 49654.99785 -4999.40151 269.66610 52378.29308 27201.04130 306 307 308 309 310 263.60403 52378.29308 41978.19496 24444.45331 52378.29308 311 312 313 314 315 52378.29308 -43764.97916 52378.29308 50952.03341 463.22713 316 317 318 319 320 -152.93183 2427.89874 -9173.81075 277.73992 269.66610 321 322 323 324 325 269.66610 269.66610 269.66610 52378.29308 269.66610 326 327 328 329 330 -2479.88831 52378.29308 52378.29308 269.66610 269.66610 331 332 333 334 335 269.66610 269.66610 269.66610 1209.65733 260.34712 336 337 338 339 340 269.66610 -1367.64739 269.66610 53.91705 2427.04406 341 342 343 344 345 -15256.62205 269.66610 269.66610 269.66610 269.66610 346 347 348 349 350 269.66610 -6619.60825 269.66610 269.66610 269.66610 351 352 353 354 355 4938.31961 -968.68749 16013.16365 -5853.03737 269.66610 356 357 358 359 360 269.66610 -37944.63660 269.66610 269.66610 269.66610 361 362 363 364 365 269.46762 269.66610 269.66610 42791.96245 -247.23550 366 367 368 369 370 269.66610 269.66610 269.66610 269.66610 269.66610 371 372 373 374 375 269.66610 -129.02984 -34857.22662 -20764.09565 269.66610 376 377 378 379 380 269.66610 5212.54708 574.48936 269.66610 269.66610 381 382 383 384 385 269.66610 -15401.57498 34164.62363 269.66610 -2228.37519 386 387 388 389 390 -89637.51538 6519.90171 11836.96694 -519.99381 -7236.11997 391 392 393 394 395 6069.64377 -8800.34051 -19588.75261 37294.67149 -157469.05530 396 397 398 399 400 -14508.02866 -16945.45973 98646.88263 -12258.55024 -105526.30134 401 402 403 404 405 -3108.63911 12572.25382 -3503.48237 -40978.81913 -16209.25828 406 407 408 409 410 -16937.34734 -1202.43237 8039.37711 -12535.91598 -51062.81077 411 412 413 414 415 -27864.30386 -217588.38755 -1660.31473 -109703.50562 -66423.50149 416 417 418 419 420 -61919.12149 160227.99049 -17015.45784 -19046.11042 -46818.19381 421 422 423 424 425 -87356.28065 -89470.05262 -626069.48963 -316767.42372 -61948.31074 426 427 428 429 430 -208492.56134 -241012.72772 -144726.09880 -121166.99639 152350.80295 431 -370254.75590 > postscript(file="/var/www/html/rcomp/tmp/6crz71291335824.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 -58916.68341 NA 1 2951183.24003 -58916.68341 2 547776.97946 2951183.24003 3 -592723.57478 547776.97946 4 -855582.68318 -592723.57478 5 -187303.82930 -855582.68318 6 1120158.08279 -187303.82930 7 422848.11045 1120158.08279 8 -703591.95988 422848.11045 9 613805.40397 -703591.95988 10 459456.04993 613805.40397 11 800071.51561 459456.04993 12 655359.62146 800071.51561 13 -468897.87437 655359.62146 14 -280661.11852 -468897.87437 15 565520.53288 -280661.11852 16 -259523.08584 565520.53288 17 -337447.59549 -259523.08584 18 -581028.38399 -337447.59549 19 52511.45770 -581028.38399 20 101447.04410 52511.45770 21 -391496.71174 101447.04410 22 120088.18252 -391496.71174 23 -127599.19151 120088.18252 24 -71606.42087 -127599.19151 25 326145.45395 -71606.42087 26 -639817.67921 326145.45395 27 -174375.28759 -639817.67921 28 621327.36936 -174375.28759 29 289543.00246 621327.36936 30 247663.51192 289543.00246 31 -168131.35366 247663.51192 32 151317.62041 -168131.35366 33 88593.24434 151317.62041 34 -311615.01065 88593.24434 35 52219.10887 -311615.01065 36 310262.12965 52219.10887 37 140406.31363 310262.12965 38 105590.98411 140406.31363 39 -133523.16907 105590.98411 40 -148996.18697 -133523.16907 41 30995.22637 -148996.18697 42 97726.72451 30995.22637 43 -53490.33926 97726.72451 44 61041.57451 -53490.33926 45 -141469.50111 61041.57451 46 47161.91860 -141469.50111 47 -196941.82015 47161.91860 48 -28247.02662 -196941.82015 49 -149590.54896 -28247.02662 50 92352.03082 -149590.54896 51 -43197.40228 92352.03082 52 -243825.49661 -43197.40228 53 75448.94087 -243825.49661 54 -283498.95034 75448.94087 55 114058.25801 -283498.95034 56 248553.73190 114058.25801 57 79226.58274 248553.73190 58 -80437.40445 79226.58274 59 -247557.75152 -80437.40445 60 -38989.03249 -247557.75152 61 -79310.46952 -38989.03249 62 -29908.81807 -79310.46952 63 -235797.07384 -29908.81807 64 62520.90448 -235797.07384 65 33957.12914 62520.90448 66 -151722.70630 33957.12914 67 -9987.33107 -151722.70630 68 -105144.12475 -9987.33107 69 -69820.92442 -105144.12475 70 -73823.45845 -69820.92442 71 166625.49732 -73823.45845 72 -145779.19403 166625.49732 73 -56387.81872 -145779.19403 74 -123106.25890 -56387.81872 75 -65954.22615 -123106.25890 76 -132389.61114 -65954.22615 77 -95381.48058 -132389.61114 78 20289.54469 -95381.48058 79 19983.58345 20289.54469 80 30180.16223 19983.58345 81 -58836.44905 30180.16223 82 -167652.48064 -58836.44905 83 50643.93943 -167652.48064 84 -48411.72652 50643.93943 85 -125132.00610 -48411.72652 86 59758.15692 -125132.00610 87 -80831.59816 59758.15692 88 -107434.97899 -80831.59816 89 -59345.23388 -107434.97899 90 -36690.09821 -59345.23388 91 -45261.40047 -36690.09821 92 110477.03904 -45261.40047 93 -256423.03336 110477.03904 94 127224.14916 -256423.03336 95 -4549.51391 127224.14916 96 -34777.47183 -4549.51391 97 24177.92868 -34777.47183 98 -107147.52429 24177.92868 99 -15983.74517 -107147.52429 100 -28627.86911 -15983.74517 101 -133068.99569 -28627.86911 102 -86388.00661 -133068.99569 103 -95303.81873 -86388.00661 104 35161.57075 -95303.81873 105 -85221.31954 35161.57075 106 -231658.23841 -85221.31954 107 99931.47560 -231658.23841 108 -98081.59614 99931.47560 109 -99669.94893 -98081.59614 110 78119.21937 -99669.94893 111 -243364.25034 78119.21937 112 -181459.97158 -243364.25034 113 -114009.58235 -181459.97158 114 36276.51649 -114009.58235 115 55746.57428 36276.51649 116 -41605.06999 55746.57428 117 53064.69920 -41605.06999 118 -170890.36224 53064.69920 119 17637.61731 -170890.36224 120 4368.30333 17637.61731 121 11382.81709 4368.30333 122 -52559.93239 11382.81709 123 68389.53214 -52559.93239 124 -92223.61851 68389.53214 125 -198868.43790 -92223.61851 126 269.66610 -198868.43790 127 492.78022 269.66610 128 9519.96729 492.78022 129 269.66610 9519.96729 130 876.92556 269.66610 131 81060.62993 876.92556 132 269.66610 81060.62993 133 -1683.22084 269.66610 134 269.66610 -1683.22084 135 9377.15855 269.66610 136 450.96100 9377.15855 137 269.66610 450.96100 138 269.66610 269.66610 139 269.66610 269.66610 140 5874.10042 269.66610 141 269.66610 5874.10042 142 21874.36284 269.66610 143 374.26243 21874.36284 144 15381.84798 374.26243 145 17301.29008 15381.84798 146 269.66610 17301.29008 147 -259.25801 269.66610 148 269.66610 -259.25801 149 -266.77164 269.66610 150 3795.40766 -266.77164 151 15573.66403 3795.40766 152 269.66610 15573.66403 153 -11688.27232 269.66610 154 -343.79348 -11688.27232 155 269.66610 -343.79348 156 269.66610 269.66610 157 -441.21905 269.66610 158 17959.78921 -441.21905 159 269.66610 17959.78921 160 269.66610 269.66610 161 269.66610 269.66610 162 269.66610 269.66610 163 9855.82870 269.66610 164 -7219.00178 9855.82870 165 269.66610 -7219.00178 166 269.66610 269.66610 167 119690.15922 269.66610 168 269.66610 119690.15922 169 269.66610 269.66610 170 269.66610 269.66610 171 -3306.88500 269.66610 172 269.66610 -3306.88500 173 818.00422 269.66610 174 269.66610 818.00422 175 -349.82655 269.66610 176 269.66610 -349.82655 177 50048.92944 269.66610 178 269.66610 50048.92944 179 269.66610 269.66610 180 32479.16452 269.66610 181 37849.42693 32479.16452 182 28406.06360 37849.42693 183 1979.83537 28406.06360 184 52212.88945 1979.83537 185 31771.85850 52212.88945 186 50634.13806 31771.85850 187 -222.76290 50634.13806 188 -33705.60342 -222.76290 189 269.66610 -33705.60342 190 12335.74127 269.66610 191 269.66610 12335.74127 192 4496.34472 269.66610 193 269.66610 4496.34472 194 26.76109 269.66610 195 52378.29308 26.76109 196 -55752.85014 52378.29308 197 52378.29308 -55752.85014 198 52378.29308 52378.29308 199 -6680.96187 52378.29308 200 52378.29308 -6680.96187 201 52378.29308 52378.29308 202 52378.29308 52378.29308 203 49347.13500 52378.29308 204 -1426.44177 49347.13500 205 51198.04438 -1426.44177 206 43133.28603 51198.04438 207 52378.29308 43133.28603 208 65274.78044 52378.29308 209 52378.29308 65274.78044 210 269.66610 52378.29308 211 269.66610 269.66610 212 10039.37072 269.66610 213 -15066.46219 10039.37072 214 -3416.17584 -15066.46219 215 269.66610 -3416.17584 216 600.33670 269.66610 217 4959.98750 600.33670 218 269.66610 4959.98750 219 31628.34903 269.66610 220 -24607.15405 31628.34903 221 269.66610 -24607.15405 222 343.94547 269.66610 223 269.66610 343.94547 224 333.61190 269.66610 225 2451.20280 333.61190 226 269.66610 2451.20280 227 269.66610 269.66610 228 -20900.79062 269.66610 229 269.66610 -20900.79062 230 -14030.37493 269.66610 231 269.66610 -14030.37493 232 269.66610 269.66610 233 6345.08101 269.66610 234 269.66610 6345.08101 235 132.21122 269.66610 236 269.66610 132.21122 237 2223.44248 269.66610 238 269.66610 2223.44248 239 269.66610 269.66610 240 269.66610 269.66610 241 5004.31516 269.66610 242 269.66610 5004.31516 243 2504.79811 269.66610 244 -8408.43581 2504.79811 245 52378.29308 -8408.43581 246 52378.29308 52378.29308 247 -8431.64237 52378.29308 248 52378.29308 -8431.64237 249 52378.29308 52378.29308 250 9587.91386 52378.29308 251 52378.29308 9587.91386 252 -18115.78580 52378.29308 253 -31116.13407 -18115.78580 254 -58439.63706 -31116.13407 255 269.66610 -58439.63706 256 31888.00688 269.66610 257 35271.83324 31888.00688 258 52378.29308 35271.83324 259 269.66610 52378.29308 260 52378.29308 269.66610 261 52378.29308 52378.29308 262 -16366.47699 52378.29308 263 12167.65305 -16366.47699 264 52378.29308 12167.65305 265 269.66610 52378.29308 266 52378.29308 269.66610 267 6138.66863 52378.29308 268 52378.29308 6138.66863 269 -18679.71641 52378.29308 270 52378.29308 -18679.71641 271 52378.29308 52378.29308 272 52378.29308 52378.29308 273 -2525.19903 52378.29308 274 269.66610 -2525.19903 275 5208.34419 269.66610 276 3639.05382 5208.34419 277 269.66610 3639.05382 278 6702.44125 269.66610 279 -27192.11395 6702.44125 280 52378.29308 -27192.11395 281 269.66610 52378.29308 282 13641.30379 269.66610 283 269.66610 13641.30379 284 -29547.57748 269.66610 285 269.66610 -29547.57748 286 52378.29308 269.66610 287 -1440.77572 52378.29308 288 52378.29308 -1440.77572 289 10112.41764 52378.29308 290 52378.29308 10112.41764 291 4671.35207 52378.29308 292 -19267.78540 4671.35207 293 269.66610 -19267.78540 294 -884.04067 269.66610 295 71481.29254 -884.04067 296 269.66610 71481.29254 297 52378.29308 269.66610 298 269.66610 52378.29308 299 269.66610 269.66610 300 49654.99785 269.66610 301 -4999.40151 49654.99785 302 269.66610 -4999.40151 303 52378.29308 269.66610 304 27201.04130 52378.29308 305 263.60403 27201.04130 306 52378.29308 263.60403 307 41978.19496 52378.29308 308 24444.45331 41978.19496 309 52378.29308 24444.45331 310 52378.29308 52378.29308 311 -43764.97916 52378.29308 312 52378.29308 -43764.97916 313 50952.03341 52378.29308 314 463.22713 50952.03341 315 -152.93183 463.22713 316 2427.89874 -152.93183 317 -9173.81075 2427.89874 318 277.73992 -9173.81075 319 269.66610 277.73992 320 269.66610 269.66610 321 269.66610 269.66610 322 269.66610 269.66610 323 52378.29308 269.66610 324 269.66610 52378.29308 325 -2479.88831 269.66610 326 52378.29308 -2479.88831 327 52378.29308 52378.29308 328 269.66610 52378.29308 329 269.66610 269.66610 330 269.66610 269.66610 331 269.66610 269.66610 332 269.66610 269.66610 333 1209.65733 269.66610 334 260.34712 1209.65733 335 269.66610 260.34712 336 -1367.64739 269.66610 337 269.66610 -1367.64739 338 53.91705 269.66610 339 2427.04406 53.91705 340 -15256.62205 2427.04406 341 269.66610 -15256.62205 342 269.66610 269.66610 343 269.66610 269.66610 344 269.66610 269.66610 345 269.66610 269.66610 346 -6619.60825 269.66610 347 269.66610 -6619.60825 348 269.66610 269.66610 349 269.66610 269.66610 350 4938.31961 269.66610 351 -968.68749 4938.31961 352 16013.16365 -968.68749 353 -5853.03737 16013.16365 354 269.66610 -5853.03737 355 269.66610 269.66610 356 -37944.63660 269.66610 357 269.66610 -37944.63660 358 269.66610 269.66610 359 269.66610 269.66610 360 269.46762 269.66610 361 269.66610 269.46762 362 269.66610 269.66610 363 42791.96245 269.66610 364 -247.23550 42791.96245 365 269.66610 -247.23550 366 269.66610 269.66610 367 269.66610 269.66610 368 269.66610 269.66610 369 269.66610 269.66610 370 269.66610 269.66610 371 -129.02984 269.66610 372 -34857.22662 -129.02984 373 -20764.09565 -34857.22662 374 269.66610 -20764.09565 375 269.66610 269.66610 376 5212.54708 269.66610 377 574.48936 5212.54708 378 269.66610 574.48936 379 269.66610 269.66610 380 269.66610 269.66610 381 -15401.57498 269.66610 382 34164.62363 -15401.57498 383 269.66610 34164.62363 384 -2228.37519 269.66610 385 -89637.51538 -2228.37519 386 6519.90171 -89637.51538 387 11836.96694 6519.90171 388 -519.99381 11836.96694 389 -7236.11997 -519.99381 390 6069.64377 -7236.11997 391 -8800.34051 6069.64377 392 -19588.75261 -8800.34051 393 37294.67149 -19588.75261 394 -157469.05530 37294.67149 395 -14508.02866 -157469.05530 396 -16945.45973 -14508.02866 397 98646.88263 -16945.45973 398 -12258.55024 98646.88263 399 -105526.30134 -12258.55024 400 -3108.63911 -105526.30134 401 12572.25382 -3108.63911 402 -3503.48237 12572.25382 403 -40978.81913 -3503.48237 404 -16209.25828 -40978.81913 405 -16937.34734 -16209.25828 406 -1202.43237 -16937.34734 407 8039.37711 -1202.43237 408 -12535.91598 8039.37711 409 -51062.81077 -12535.91598 410 -27864.30386 -51062.81077 411 -217588.38755 -27864.30386 412 -1660.31473 -217588.38755 413 -109703.50562 -1660.31473 414 -66423.50149 -109703.50562 415 -61919.12149 -66423.50149 416 160227.99049 -61919.12149 417 -17015.45784 160227.99049 418 -19046.11042 -17015.45784 419 -46818.19381 -19046.11042 420 -87356.28065 -46818.19381 421 -89470.05262 -87356.28065 422 -626069.48963 -89470.05262 423 -316767.42372 -626069.48963 424 -61948.31074 -316767.42372 425 -208492.56134 -61948.31074 426 -241012.72772 -208492.56134 427 -144726.09880 -241012.72772 428 -121166.99639 -144726.09880 429 152350.80295 -121166.99639 430 -370254.75590 152350.80295 431 NA -370254.75590 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2951183.24003 -58916.68341 [2,] 547776.97946 2951183.24003 [3,] -592723.57478 547776.97946 [4,] -855582.68318 -592723.57478 [5,] -187303.82930 -855582.68318 [6,] 1120158.08279 -187303.82930 [7,] 422848.11045 1120158.08279 [8,] -703591.95988 422848.11045 [9,] 613805.40397 -703591.95988 [10,] 459456.04993 613805.40397 [11,] 800071.51561 459456.04993 [12,] 655359.62146 800071.51561 [13,] -468897.87437 655359.62146 [14,] -280661.11852 -468897.87437 [15,] 565520.53288 -280661.11852 [16,] -259523.08584 565520.53288 [17,] -337447.59549 -259523.08584 [18,] -581028.38399 -337447.59549 [19,] 52511.45770 -581028.38399 [20,] 101447.04410 52511.45770 [21,] -391496.71174 101447.04410 [22,] 120088.18252 -391496.71174 [23,] -127599.19151 120088.18252 [24,] -71606.42087 -127599.19151 [25,] 326145.45395 -71606.42087 [26,] -639817.67921 326145.45395 [27,] -174375.28759 -639817.67921 [28,] 621327.36936 -174375.28759 [29,] 289543.00246 621327.36936 [30,] 247663.51192 289543.00246 [31,] -168131.35366 247663.51192 [32,] 151317.62041 -168131.35366 [33,] 88593.24434 151317.62041 [34,] -311615.01065 88593.24434 [35,] 52219.10887 -311615.01065 [36,] 310262.12965 52219.10887 [37,] 140406.31363 310262.12965 [38,] 105590.98411 140406.31363 [39,] -133523.16907 105590.98411 [40,] -148996.18697 -133523.16907 [41,] 30995.22637 -148996.18697 [42,] 97726.72451 30995.22637 [43,] -53490.33926 97726.72451 [44,] 61041.57451 -53490.33926 [45,] -141469.50111 61041.57451 [46,] 47161.91860 -141469.50111 [47,] -196941.82015 47161.91860 [48,] -28247.02662 -196941.82015 [49,] -149590.54896 -28247.02662 [50,] 92352.03082 -149590.54896 [51,] -43197.40228 92352.03082 [52,] -243825.49661 -43197.40228 [53,] 75448.94087 -243825.49661 [54,] -283498.95034 75448.94087 [55,] 114058.25801 -283498.95034 [56,] 248553.73190 114058.25801 [57,] 79226.58274 248553.73190 [58,] -80437.40445 79226.58274 [59,] -247557.75152 -80437.40445 [60,] -38989.03249 -247557.75152 [61,] -79310.46952 -38989.03249 [62,] -29908.81807 -79310.46952 [63,] -235797.07384 -29908.81807 [64,] 62520.90448 -235797.07384 [65,] 33957.12914 62520.90448 [66,] -151722.70630 33957.12914 [67,] -9987.33107 -151722.70630 [68,] -105144.12475 -9987.33107 [69,] -69820.92442 -105144.12475 [70,] -73823.45845 -69820.92442 [71,] 166625.49732 -73823.45845 [72,] -145779.19403 166625.49732 [73,] -56387.81872 -145779.19403 [74,] -123106.25890 -56387.81872 [75,] -65954.22615 -123106.25890 [76,] -132389.61114 -65954.22615 [77,] -95381.48058 -132389.61114 [78,] 20289.54469 -95381.48058 [79,] 19983.58345 20289.54469 [80,] 30180.16223 19983.58345 [81,] -58836.44905 30180.16223 [82,] -167652.48064 -58836.44905 [83,] 50643.93943 -167652.48064 [84,] -48411.72652 50643.93943 [85,] -125132.00610 -48411.72652 [86,] 59758.15692 -125132.00610 [87,] -80831.59816 59758.15692 [88,] -107434.97899 -80831.59816 [89,] -59345.23388 -107434.97899 [90,] -36690.09821 -59345.23388 [91,] -45261.40047 -36690.09821 [92,] 110477.03904 -45261.40047 [93,] -256423.03336 110477.03904 [94,] 127224.14916 -256423.03336 [95,] -4549.51391 127224.14916 [96,] -34777.47183 -4549.51391 [97,] 24177.92868 -34777.47183 [98,] -107147.52429 24177.92868 [99,] -15983.74517 -107147.52429 [100,] -28627.86911 -15983.74517 [101,] -133068.99569 -28627.86911 [102,] -86388.00661 -133068.99569 [103,] -95303.81873 -86388.00661 [104,] 35161.57075 -95303.81873 [105,] -85221.31954 35161.57075 [106,] -231658.23841 -85221.31954 [107,] 99931.47560 -231658.23841 [108,] -98081.59614 99931.47560 [109,] -99669.94893 -98081.59614 [110,] 78119.21937 -99669.94893 [111,] -243364.25034 78119.21937 [112,] -181459.97158 -243364.25034 [113,] -114009.58235 -181459.97158 [114,] 36276.51649 -114009.58235 [115,] 55746.57428 36276.51649 [116,] -41605.06999 55746.57428 [117,] 53064.69920 -41605.06999 [118,] -170890.36224 53064.69920 [119,] 17637.61731 -170890.36224 [120,] 4368.30333 17637.61731 [121,] 11382.81709 4368.30333 [122,] -52559.93239 11382.81709 [123,] 68389.53214 -52559.93239 [124,] -92223.61851 68389.53214 [125,] -198868.43790 -92223.61851 [126,] 269.66610 -198868.43790 [127,] 492.78022 269.66610 [128,] 9519.96729 492.78022 [129,] 269.66610 9519.96729 [130,] 876.92556 269.66610 [131,] 81060.62993 876.92556 [132,] 269.66610 81060.62993 [133,] -1683.22084 269.66610 [134,] 269.66610 -1683.22084 [135,] 9377.15855 269.66610 [136,] 450.96100 9377.15855 [137,] 269.66610 450.96100 [138,] 269.66610 269.66610 [139,] 269.66610 269.66610 [140,] 5874.10042 269.66610 [141,] 269.66610 5874.10042 [142,] 21874.36284 269.66610 [143,] 374.26243 21874.36284 [144,] 15381.84798 374.26243 [145,] 17301.29008 15381.84798 [146,] 269.66610 17301.29008 [147,] -259.25801 269.66610 [148,] 269.66610 -259.25801 [149,] -266.77164 269.66610 [150,] 3795.40766 -266.77164 [151,] 15573.66403 3795.40766 [152,] 269.66610 15573.66403 [153,] -11688.27232 269.66610 [154,] -343.79348 -11688.27232 [155,] 269.66610 -343.79348 [156,] 269.66610 269.66610 [157,] -441.21905 269.66610 [158,] 17959.78921 -441.21905 [159,] 269.66610 17959.78921 [160,] 269.66610 269.66610 [161,] 269.66610 269.66610 [162,] 269.66610 269.66610 [163,] 9855.82870 269.66610 [164,] -7219.00178 9855.82870 [165,] 269.66610 -7219.00178 [166,] 269.66610 269.66610 [167,] 119690.15922 269.66610 [168,] 269.66610 119690.15922 [169,] 269.66610 269.66610 [170,] 269.66610 269.66610 [171,] -3306.88500 269.66610 [172,] 269.66610 -3306.88500 [173,] 818.00422 269.66610 [174,] 269.66610 818.00422 [175,] -349.82655 269.66610 [176,] 269.66610 -349.82655 [177,] 50048.92944 269.66610 [178,] 269.66610 50048.92944 [179,] 269.66610 269.66610 [180,] 32479.16452 269.66610 [181,] 37849.42693 32479.16452 [182,] 28406.06360 37849.42693 [183,] 1979.83537 28406.06360 [184,] 52212.88945 1979.83537 [185,] 31771.85850 52212.88945 [186,] 50634.13806 31771.85850 [187,] -222.76290 50634.13806 [188,] -33705.60342 -222.76290 [189,] 269.66610 -33705.60342 [190,] 12335.74127 269.66610 [191,] 269.66610 12335.74127 [192,] 4496.34472 269.66610 [193,] 269.66610 4496.34472 [194,] 26.76109 269.66610 [195,] 52378.29308 26.76109 [196,] -55752.85014 52378.29308 [197,] 52378.29308 -55752.85014 [198,] 52378.29308 52378.29308 [199,] -6680.96187 52378.29308 [200,] 52378.29308 -6680.96187 [201,] 52378.29308 52378.29308 [202,] 52378.29308 52378.29308 [203,] 49347.13500 52378.29308 [204,] -1426.44177 49347.13500 [205,] 51198.04438 -1426.44177 [206,] 43133.28603 51198.04438 [207,] 52378.29308 43133.28603 [208,] 65274.78044 52378.29308 [209,] 52378.29308 65274.78044 [210,] 269.66610 52378.29308 [211,] 269.66610 269.66610 [212,] 10039.37072 269.66610 [213,] -15066.46219 10039.37072 [214,] -3416.17584 -15066.46219 [215,] 269.66610 -3416.17584 [216,] 600.33670 269.66610 [217,] 4959.98750 600.33670 [218,] 269.66610 4959.98750 [219,] 31628.34903 269.66610 [220,] -24607.15405 31628.34903 [221,] 269.66610 -24607.15405 [222,] 343.94547 269.66610 [223,] 269.66610 343.94547 [224,] 333.61190 269.66610 [225,] 2451.20280 333.61190 [226,] 269.66610 2451.20280 [227,] 269.66610 269.66610 [228,] -20900.79062 269.66610 [229,] 269.66610 -20900.79062 [230,] -14030.37493 269.66610 [231,] 269.66610 -14030.37493 [232,] 269.66610 269.66610 [233,] 6345.08101 269.66610 [234,] 269.66610 6345.08101 [235,] 132.21122 269.66610 [236,] 269.66610 132.21122 [237,] 2223.44248 269.66610 [238,] 269.66610 2223.44248 [239,] 269.66610 269.66610 [240,] 269.66610 269.66610 [241,] 5004.31516 269.66610 [242,] 269.66610 5004.31516 [243,] 2504.79811 269.66610 [244,] -8408.43581 2504.79811 [245,] 52378.29308 -8408.43581 [246,] 52378.29308 52378.29308 [247,] -8431.64237 52378.29308 [248,] 52378.29308 -8431.64237 [249,] 52378.29308 52378.29308 [250,] 9587.91386 52378.29308 [251,] 52378.29308 9587.91386 [252,] -18115.78580 52378.29308 [253,] -31116.13407 -18115.78580 [254,] -58439.63706 -31116.13407 [255,] 269.66610 -58439.63706 [256,] 31888.00688 269.66610 [257,] 35271.83324 31888.00688 [258,] 52378.29308 35271.83324 [259,] 269.66610 52378.29308 [260,] 52378.29308 269.66610 [261,] 52378.29308 52378.29308 [262,] -16366.47699 52378.29308 [263,] 12167.65305 -16366.47699 [264,] 52378.29308 12167.65305 [265,] 269.66610 52378.29308 [266,] 52378.29308 269.66610 [267,] 6138.66863 52378.29308 [268,] 52378.29308 6138.66863 [269,] -18679.71641 52378.29308 [270,] 52378.29308 -18679.71641 [271,] 52378.29308 52378.29308 [272,] 52378.29308 52378.29308 [273,] -2525.19903 52378.29308 [274,] 269.66610 -2525.19903 [275,] 5208.34419 269.66610 [276,] 3639.05382 5208.34419 [277,] 269.66610 3639.05382 [278,] 6702.44125 269.66610 [279,] -27192.11395 6702.44125 [280,] 52378.29308 -27192.11395 [281,] 269.66610 52378.29308 [282,] 13641.30379 269.66610 [283,] 269.66610 13641.30379 [284,] -29547.57748 269.66610 [285,] 269.66610 -29547.57748 [286,] 52378.29308 269.66610 [287,] -1440.77572 52378.29308 [288,] 52378.29308 -1440.77572 [289,] 10112.41764 52378.29308 [290,] 52378.29308 10112.41764 [291,] 4671.35207 52378.29308 [292,] -19267.78540 4671.35207 [293,] 269.66610 -19267.78540 [294,] -884.04067 269.66610 [295,] 71481.29254 -884.04067 [296,] 269.66610 71481.29254 [297,] 52378.29308 269.66610 [298,] 269.66610 52378.29308 [299,] 269.66610 269.66610 [300,] 49654.99785 269.66610 [301,] -4999.40151 49654.99785 [302,] 269.66610 -4999.40151 [303,] 52378.29308 269.66610 [304,] 27201.04130 52378.29308 [305,] 263.60403 27201.04130 [306,] 52378.29308 263.60403 [307,] 41978.19496 52378.29308 [308,] 24444.45331 41978.19496 [309,] 52378.29308 24444.45331 [310,] 52378.29308 52378.29308 [311,] -43764.97916 52378.29308 [312,] 52378.29308 -43764.97916 [313,] 50952.03341 52378.29308 [314,] 463.22713 50952.03341 [315,] -152.93183 463.22713 [316,] 2427.89874 -152.93183 [317,] -9173.81075 2427.89874 [318,] 277.73992 -9173.81075 [319,] 269.66610 277.73992 [320,] 269.66610 269.66610 [321,] 269.66610 269.66610 [322,] 269.66610 269.66610 [323,] 52378.29308 269.66610 [324,] 269.66610 52378.29308 [325,] -2479.88831 269.66610 [326,] 52378.29308 -2479.88831 [327,] 52378.29308 52378.29308 [328,] 269.66610 52378.29308 [329,] 269.66610 269.66610 [330,] 269.66610 269.66610 [331,] 269.66610 269.66610 [332,] 269.66610 269.66610 [333,] 1209.65733 269.66610 [334,] 260.34712 1209.65733 [335,] 269.66610 260.34712 [336,] -1367.64739 269.66610 [337,] 269.66610 -1367.64739 [338,] 53.91705 269.66610 [339,] 2427.04406 53.91705 [340,] -15256.62205 2427.04406 [341,] 269.66610 -15256.62205 [342,] 269.66610 269.66610 [343,] 269.66610 269.66610 [344,] 269.66610 269.66610 [345,] 269.66610 269.66610 [346,] -6619.60825 269.66610 [347,] 269.66610 -6619.60825 [348,] 269.66610 269.66610 [349,] 269.66610 269.66610 [350,] 4938.31961 269.66610 [351,] -968.68749 4938.31961 [352,] 16013.16365 -968.68749 [353,] -5853.03737 16013.16365 [354,] 269.66610 -5853.03737 [355,] 269.66610 269.66610 [356,] -37944.63660 269.66610 [357,] 269.66610 -37944.63660 [358,] 269.66610 269.66610 [359,] 269.66610 269.66610 [360,] 269.46762 269.66610 [361,] 269.66610 269.46762 [362,] 269.66610 269.66610 [363,] 42791.96245 269.66610 [364,] -247.23550 42791.96245 [365,] 269.66610 -247.23550 [366,] 269.66610 269.66610 [367,] 269.66610 269.66610 [368,] 269.66610 269.66610 [369,] 269.66610 269.66610 [370,] 269.66610 269.66610 [371,] -129.02984 269.66610 [372,] -34857.22662 -129.02984 [373,] -20764.09565 -34857.22662 [374,] 269.66610 -20764.09565 [375,] 269.66610 269.66610 [376,] 5212.54708 269.66610 [377,] 574.48936 5212.54708 [378,] 269.66610 574.48936 [379,] 269.66610 269.66610 [380,] 269.66610 269.66610 [381,] -15401.57498 269.66610 [382,] 34164.62363 -15401.57498 [383,] 269.66610 34164.62363 [384,] -2228.37519 269.66610 [385,] -89637.51538 -2228.37519 [386,] 6519.90171 -89637.51538 [387,] 11836.96694 6519.90171 [388,] -519.99381 11836.96694 [389,] -7236.11997 -519.99381 [390,] 6069.64377 -7236.11997 [391,] -8800.34051 6069.64377 [392,] -19588.75261 -8800.34051 [393,] 37294.67149 -19588.75261 [394,] -157469.05530 37294.67149 [395,] -14508.02866 -157469.05530 [396,] -16945.45973 -14508.02866 [397,] 98646.88263 -16945.45973 [398,] -12258.55024 98646.88263 [399,] -105526.30134 -12258.55024 [400,] -3108.63911 -105526.30134 [401,] 12572.25382 -3108.63911 [402,] -3503.48237 12572.25382 [403,] -40978.81913 -3503.48237 [404,] -16209.25828 -40978.81913 [405,] -16937.34734 -16209.25828 [406,] -1202.43237 -16937.34734 [407,] 8039.37711 -1202.43237 [408,] -12535.91598 8039.37711 [409,] -51062.81077 -12535.91598 [410,] -27864.30386 -51062.81077 [411,] -217588.38755 -27864.30386 [412,] -1660.31473 -217588.38755 [413,] -109703.50562 -1660.31473 [414,] -66423.50149 -109703.50562 [415,] -61919.12149 -66423.50149 [416,] 160227.99049 -61919.12149 [417,] -17015.45784 160227.99049 [418,] -19046.11042 -17015.45784 [419,] -46818.19381 -19046.11042 [420,] -87356.28065 -46818.19381 [421,] -89470.05262 -87356.28065 [422,] -626069.48963 -89470.05262 [423,] -316767.42372 -626069.48963 [424,] -61948.31074 -316767.42372 [425,] -208492.56134 -61948.31074 [426,] -241012.72772 -208492.56134 [427,] -144726.09880 -241012.72772 [428,] -121166.99639 -144726.09880 [429,] 152350.80295 -121166.99639 [430,] -370254.75590 152350.80295 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2951183.24003 -58916.68341 2 547776.97946 2951183.24003 3 -592723.57478 547776.97946 4 -855582.68318 -592723.57478 5 -187303.82930 -855582.68318 6 1120158.08279 -187303.82930 7 422848.11045 1120158.08279 8 -703591.95988 422848.11045 9 613805.40397 -703591.95988 10 459456.04993 613805.40397 11 800071.51561 459456.04993 12 655359.62146 800071.51561 13 -468897.87437 655359.62146 14 -280661.11852 -468897.87437 15 565520.53288 -280661.11852 16 -259523.08584 565520.53288 17 -337447.59549 -259523.08584 18 -581028.38399 -337447.59549 19 52511.45770 -581028.38399 20 101447.04410 52511.45770 21 -391496.71174 101447.04410 22 120088.18252 -391496.71174 23 -127599.19151 120088.18252 24 -71606.42087 -127599.19151 25 326145.45395 -71606.42087 26 -639817.67921 326145.45395 27 -174375.28759 -639817.67921 28 621327.36936 -174375.28759 29 289543.00246 621327.36936 30 247663.51192 289543.00246 31 -168131.35366 247663.51192 32 151317.62041 -168131.35366 33 88593.24434 151317.62041 34 -311615.01065 88593.24434 35 52219.10887 -311615.01065 36 310262.12965 52219.10887 37 140406.31363 310262.12965 38 105590.98411 140406.31363 39 -133523.16907 105590.98411 40 -148996.18697 -133523.16907 41 30995.22637 -148996.18697 42 97726.72451 30995.22637 43 -53490.33926 97726.72451 44 61041.57451 -53490.33926 45 -141469.50111 61041.57451 46 47161.91860 -141469.50111 47 -196941.82015 47161.91860 48 -28247.02662 -196941.82015 49 -149590.54896 -28247.02662 50 92352.03082 -149590.54896 51 -43197.40228 92352.03082 52 -243825.49661 -43197.40228 53 75448.94087 -243825.49661 54 -283498.95034 75448.94087 55 114058.25801 -283498.95034 56 248553.73190 114058.25801 57 79226.58274 248553.73190 58 -80437.40445 79226.58274 59 -247557.75152 -80437.40445 60 -38989.03249 -247557.75152 61 -79310.46952 -38989.03249 62 -29908.81807 -79310.46952 63 -235797.07384 -29908.81807 64 62520.90448 -235797.07384 65 33957.12914 62520.90448 66 -151722.70630 33957.12914 67 -9987.33107 -151722.70630 68 -105144.12475 -9987.33107 69 -69820.92442 -105144.12475 70 -73823.45845 -69820.92442 71 166625.49732 -73823.45845 72 -145779.19403 166625.49732 73 -56387.81872 -145779.19403 74 -123106.25890 -56387.81872 75 -65954.22615 -123106.25890 76 -132389.61114 -65954.22615 77 -95381.48058 -132389.61114 78 20289.54469 -95381.48058 79 19983.58345 20289.54469 80 30180.16223 19983.58345 81 -58836.44905 30180.16223 82 -167652.48064 -58836.44905 83 50643.93943 -167652.48064 84 -48411.72652 50643.93943 85 -125132.00610 -48411.72652 86 59758.15692 -125132.00610 87 -80831.59816 59758.15692 88 -107434.97899 -80831.59816 89 -59345.23388 -107434.97899 90 -36690.09821 -59345.23388 91 -45261.40047 -36690.09821 92 110477.03904 -45261.40047 93 -256423.03336 110477.03904 94 127224.14916 -256423.03336 95 -4549.51391 127224.14916 96 -34777.47183 -4549.51391 97 24177.92868 -34777.47183 98 -107147.52429 24177.92868 99 -15983.74517 -107147.52429 100 -28627.86911 -15983.74517 101 -133068.99569 -28627.86911 102 -86388.00661 -133068.99569 103 -95303.81873 -86388.00661 104 35161.57075 -95303.81873 105 -85221.31954 35161.57075 106 -231658.23841 -85221.31954 107 99931.47560 -231658.23841 108 -98081.59614 99931.47560 109 -99669.94893 -98081.59614 110 78119.21937 -99669.94893 111 -243364.25034 78119.21937 112 -181459.97158 -243364.25034 113 -114009.58235 -181459.97158 114 36276.51649 -114009.58235 115 55746.57428 36276.51649 116 -41605.06999 55746.57428 117 53064.69920 -41605.06999 118 -170890.36224 53064.69920 119 17637.61731 -170890.36224 120 4368.30333 17637.61731 121 11382.81709 4368.30333 122 -52559.93239 11382.81709 123 68389.53214 -52559.93239 124 -92223.61851 68389.53214 125 -198868.43790 -92223.61851 126 269.66610 -198868.43790 127 492.78022 269.66610 128 9519.96729 492.78022 129 269.66610 9519.96729 130 876.92556 269.66610 131 81060.62993 876.92556 132 269.66610 81060.62993 133 -1683.22084 269.66610 134 269.66610 -1683.22084 135 9377.15855 269.66610 136 450.96100 9377.15855 137 269.66610 450.96100 138 269.66610 269.66610 139 269.66610 269.66610 140 5874.10042 269.66610 141 269.66610 5874.10042 142 21874.36284 269.66610 143 374.26243 21874.36284 144 15381.84798 374.26243 145 17301.29008 15381.84798 146 269.66610 17301.29008 147 -259.25801 269.66610 148 269.66610 -259.25801 149 -266.77164 269.66610 150 3795.40766 -266.77164 151 15573.66403 3795.40766 152 269.66610 15573.66403 153 -11688.27232 269.66610 154 -343.79348 -11688.27232 155 269.66610 -343.79348 156 269.66610 269.66610 157 -441.21905 269.66610 158 17959.78921 -441.21905 159 269.66610 17959.78921 160 269.66610 269.66610 161 269.66610 269.66610 162 269.66610 269.66610 163 9855.82870 269.66610 164 -7219.00178 9855.82870 165 269.66610 -7219.00178 166 269.66610 269.66610 167 119690.15922 269.66610 168 269.66610 119690.15922 169 269.66610 269.66610 170 269.66610 269.66610 171 -3306.88500 269.66610 172 269.66610 -3306.88500 173 818.00422 269.66610 174 269.66610 818.00422 175 -349.82655 269.66610 176 269.66610 -349.82655 177 50048.92944 269.66610 178 269.66610 50048.92944 179 269.66610 269.66610 180 32479.16452 269.66610 181 37849.42693 32479.16452 182 28406.06360 37849.42693 183 1979.83537 28406.06360 184 52212.88945 1979.83537 185 31771.85850 52212.88945 186 50634.13806 31771.85850 187 -222.76290 50634.13806 188 -33705.60342 -222.76290 189 269.66610 -33705.60342 190 12335.74127 269.66610 191 269.66610 12335.74127 192 4496.34472 269.66610 193 269.66610 4496.34472 194 26.76109 269.66610 195 52378.29308 26.76109 196 -55752.85014 52378.29308 197 52378.29308 -55752.85014 198 52378.29308 52378.29308 199 -6680.96187 52378.29308 200 52378.29308 -6680.96187 201 52378.29308 52378.29308 202 52378.29308 52378.29308 203 49347.13500 52378.29308 204 -1426.44177 49347.13500 205 51198.04438 -1426.44177 206 43133.28603 51198.04438 207 52378.29308 43133.28603 208 65274.78044 52378.29308 209 52378.29308 65274.78044 210 269.66610 52378.29308 211 269.66610 269.66610 212 10039.37072 269.66610 213 -15066.46219 10039.37072 214 -3416.17584 -15066.46219 215 269.66610 -3416.17584 216 600.33670 269.66610 217 4959.98750 600.33670 218 269.66610 4959.98750 219 31628.34903 269.66610 220 -24607.15405 31628.34903 221 269.66610 -24607.15405 222 343.94547 269.66610 223 269.66610 343.94547 224 333.61190 269.66610 225 2451.20280 333.61190 226 269.66610 2451.20280 227 269.66610 269.66610 228 -20900.79062 269.66610 229 269.66610 -20900.79062 230 -14030.37493 269.66610 231 269.66610 -14030.37493 232 269.66610 269.66610 233 6345.08101 269.66610 234 269.66610 6345.08101 235 132.21122 269.66610 236 269.66610 132.21122 237 2223.44248 269.66610 238 269.66610 2223.44248 239 269.66610 269.66610 240 269.66610 269.66610 241 5004.31516 269.66610 242 269.66610 5004.31516 243 2504.79811 269.66610 244 -8408.43581 2504.79811 245 52378.29308 -8408.43581 246 52378.29308 52378.29308 247 -8431.64237 52378.29308 248 52378.29308 -8431.64237 249 52378.29308 52378.29308 250 9587.91386 52378.29308 251 52378.29308 9587.91386 252 -18115.78580 52378.29308 253 -31116.13407 -18115.78580 254 -58439.63706 -31116.13407 255 269.66610 -58439.63706 256 31888.00688 269.66610 257 35271.83324 31888.00688 258 52378.29308 35271.83324 259 269.66610 52378.29308 260 52378.29308 269.66610 261 52378.29308 52378.29308 262 -16366.47699 52378.29308 263 12167.65305 -16366.47699 264 52378.29308 12167.65305 265 269.66610 52378.29308 266 52378.29308 269.66610 267 6138.66863 52378.29308 268 52378.29308 6138.66863 269 -18679.71641 52378.29308 270 52378.29308 -18679.71641 271 52378.29308 52378.29308 272 52378.29308 52378.29308 273 -2525.19903 52378.29308 274 269.66610 -2525.19903 275 5208.34419 269.66610 276 3639.05382 5208.34419 277 269.66610 3639.05382 278 6702.44125 269.66610 279 -27192.11395 6702.44125 280 52378.29308 -27192.11395 281 269.66610 52378.29308 282 13641.30379 269.66610 283 269.66610 13641.30379 284 -29547.57748 269.66610 285 269.66610 -29547.57748 286 52378.29308 269.66610 287 -1440.77572 52378.29308 288 52378.29308 -1440.77572 289 10112.41764 52378.29308 290 52378.29308 10112.41764 291 4671.35207 52378.29308 292 -19267.78540 4671.35207 293 269.66610 -19267.78540 294 -884.04067 269.66610 295 71481.29254 -884.04067 296 269.66610 71481.29254 297 52378.29308 269.66610 298 269.66610 52378.29308 299 269.66610 269.66610 300 49654.99785 269.66610 301 -4999.40151 49654.99785 302 269.66610 -4999.40151 303 52378.29308 269.66610 304 27201.04130 52378.29308 305 263.60403 27201.04130 306 52378.29308 263.60403 307 41978.19496 52378.29308 308 24444.45331 41978.19496 309 52378.29308 24444.45331 310 52378.29308 52378.29308 311 -43764.97916 52378.29308 312 52378.29308 -43764.97916 313 50952.03341 52378.29308 314 463.22713 50952.03341 315 -152.93183 463.22713 316 2427.89874 -152.93183 317 -9173.81075 2427.89874 318 277.73992 -9173.81075 319 269.66610 277.73992 320 269.66610 269.66610 321 269.66610 269.66610 322 269.66610 269.66610 323 52378.29308 269.66610 324 269.66610 52378.29308 325 -2479.88831 269.66610 326 52378.29308 -2479.88831 327 52378.29308 52378.29308 328 269.66610 52378.29308 329 269.66610 269.66610 330 269.66610 269.66610 331 269.66610 269.66610 332 269.66610 269.66610 333 1209.65733 269.66610 334 260.34712 1209.65733 335 269.66610 260.34712 336 -1367.64739 269.66610 337 269.66610 -1367.64739 338 53.91705 269.66610 339 2427.04406 53.91705 340 -15256.62205 2427.04406 341 269.66610 -15256.62205 342 269.66610 269.66610 343 269.66610 269.66610 344 269.66610 269.66610 345 269.66610 269.66610 346 -6619.60825 269.66610 347 269.66610 -6619.60825 348 269.66610 269.66610 349 269.66610 269.66610 350 4938.31961 269.66610 351 -968.68749 4938.31961 352 16013.16365 -968.68749 353 -5853.03737 16013.16365 354 269.66610 -5853.03737 355 269.66610 269.66610 356 -37944.63660 269.66610 357 269.66610 -37944.63660 358 269.66610 269.66610 359 269.66610 269.66610 360 269.46762 269.66610 361 269.66610 269.46762 362 269.66610 269.66610 363 42791.96245 269.66610 364 -247.23550 42791.96245 365 269.66610 -247.23550 366 269.66610 269.66610 367 269.66610 269.66610 368 269.66610 269.66610 369 269.66610 269.66610 370 269.66610 269.66610 371 -129.02984 269.66610 372 -34857.22662 -129.02984 373 -20764.09565 -34857.22662 374 269.66610 -20764.09565 375 269.66610 269.66610 376 5212.54708 269.66610 377 574.48936 5212.54708 378 269.66610 574.48936 379 269.66610 269.66610 380 269.66610 269.66610 381 -15401.57498 269.66610 382 34164.62363 -15401.57498 383 269.66610 34164.62363 384 -2228.37519 269.66610 385 -89637.51538 -2228.37519 386 6519.90171 -89637.51538 387 11836.96694 6519.90171 388 -519.99381 11836.96694 389 -7236.11997 -519.99381 390 6069.64377 -7236.11997 391 -8800.34051 6069.64377 392 -19588.75261 -8800.34051 393 37294.67149 -19588.75261 394 -157469.05530 37294.67149 395 -14508.02866 -157469.05530 396 -16945.45973 -14508.02866 397 98646.88263 -16945.45973 398 -12258.55024 98646.88263 399 -105526.30134 -12258.55024 400 -3108.63911 -105526.30134 401 12572.25382 -3108.63911 402 -3503.48237 12572.25382 403 -40978.81913 -3503.48237 404 -16209.25828 -40978.81913 405 -16937.34734 -16209.25828 406 -1202.43237 -16937.34734 407 8039.37711 -1202.43237 408 -12535.91598 8039.37711 409 -51062.81077 -12535.91598 410 -27864.30386 -51062.81077 411 -217588.38755 -27864.30386 412 -1660.31473 -217588.38755 413 -109703.50562 -1660.31473 414 -66423.50149 -109703.50562 415 -61919.12149 -66423.50149 416 160227.99049 -61919.12149 417 -17015.45784 160227.99049 418 -19046.11042 -17015.45784 419 -46818.19381 -19046.11042 420 -87356.28065 -46818.19381 421 -89470.05262 -87356.28065 422 -626069.48963 -89470.05262 423 -316767.42372 -626069.48963 424 -61948.31074 -316767.42372 425 -208492.56134 -61948.31074 426 -241012.72772 -208492.56134 427 -144726.09880 -241012.72772 428 -121166.99639 -144726.09880 429 152350.80295 -121166.99639 430 -370254.75590 152350.80295 > 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/75izs1291335824.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/8f9gv1291335824.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/9f9gv1291335824.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/10f9gv1291335824.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/11t1e41291335824.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/12eky11291335825.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/133lvd1291335825.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/14ecug1291335825.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/15zdbm1291335825.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/16dn8d1291335825.tab") + } > > try(system("convert tmp/1jzi41291335824.ps tmp/1jzi41291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/2jzi41291335824.ps tmp/2jzi41291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/3jzi41291335824.ps tmp/3jzi41291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/4crz71291335824.ps tmp/4crz71291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/5crz71291335824.ps tmp/5crz71291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/6crz71291335824.ps tmp/6crz71291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/75izs1291335824.ps tmp/75izs1291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/8f9gv1291335824.ps tmp/8f9gv1291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/9f9gv1291335824.ps tmp/9f9gv1291335824.png",intern=TRUE)) character(0) > try(system("convert tmp/10f9gv1291335824.ps tmp/10f9gv1291335824.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.810 2.196 26.236