R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1 + ,162556 + ,162556 + ,1081 + ,1081 + ,213118 + ,213118 + ,230380558 + ,6282929 + ,1 + ,29790 + ,29790 + ,309 + ,309 + ,81767 + ,81767 + ,25266003 + ,4324047 + ,1 + ,87550 + ,87550 + ,458 + ,458 + ,153198 + ,153198 + ,70164684 + ,4108272 + ,0 + ,84738 + ,0 + ,588 + ,0 + ,-26007 + ,0 + ,-15292116 + ,-1212617 + ,1 + ,54660 + ,54660 + ,299 + ,299 + ,126942 + ,126942 + ,37955658 + ,1485329 + ,1 + ,42634 + ,42634 + ,156 + ,156 + ,157214 + ,157214 + ,24525384 + ,1779876 + ,0 + ,40949 + ,0 + ,481 + ,0 + ,129352 + ,0 + ,62218312 + ,1367203 + ,1 + ,42312 + ,42312 + ,323 + ,323 + ,234817 + ,234817 + ,75845891 + ,2519076 + ,1 + ,37704 + ,37704 + ,452 + ,452 + ,60448 + ,60448 + ,27322496 + ,912684 + ,1 + ,16275 + ,16275 + ,109 + ,109 + ,47818 + ,47818 + ,5212162 + ,1443586 + ,0 + ,25830 + ,0 + ,115 + ,0 + ,245546 + ,0 + ,28237790 + ,1220017 + ,0 + ,12679 + ,0 + ,110 + ,0 + ,48020 + ,0 + ,5282200 + ,984885 + ,1 + ,18014 + ,18014 + ,239 + ,239 + ,-1710 + ,-1710 + ,-408690 + ,1457425 + ,0 + ,43556 + ,0 + ,247 + ,0 + ,32648 + ,0 + ,8064056 + ,-572920 + ,1 + ,24524 + ,24524 + ,497 + ,497 + ,95350 + ,95350 + ,47388950 + ,929144 + ,0 + ,6532 + ,0 + ,103 + ,0 + ,151352 + ,0 + ,15589256 + ,1151176 + ,0 + ,7123 + ,0 + ,109 + ,0 + ,288170 + ,0 + ,31410530 + ,790090 + ,1 + ,20813 + ,20813 + ,502 + ,502 + ,114337 + ,114337 + ,57397174 + ,774497 + ,1 + ,37597 + ,37597 + ,248 + ,248 + ,37884 + ,37884 + ,9395232 + ,990576 + ,0 + ,17821 + ,0 + ,373 + ,0 + ,122844 + ,0 + ,45820812 + ,454195 + ,1 + ,12988 + ,12988 + ,119 + ,119 + ,82340 + ,82340 + ,9798460 + ,876607 + ,1 + ,22330 + ,22330 + ,84 + ,84 + ,79801 + ,79801 + ,6703284 + ,711969 + ,0 + ,13326 + ,0 + ,102 + ,0 + ,165548 + ,0 + ,16885896 + ,702380 + ,0 + ,16189 + ,0 + ,295 + ,0 + ,116384 + ,0 + ,34333280 + ,264449 + ,0 + ,7146 + ,0 + ,105 + ,0 + ,134028 + ,0 + ,14072940 + ,450033 + ,0 + ,15824 + ,0 + ,64 + ,0 + ,63838 + ,0 + ,4085632 + ,541063 + ,1 + ,26088 + ,26088 + ,267 + ,267 + ,74996 + ,74996 + ,20023932 + ,588864 + ,0 + ,11326 + ,0 + ,129 + ,0 + ,31080 + ,0 + ,4009320 + ,-37216 + ,0 + ,8568 + ,0 + ,37 + ,0 + ,32168 + ,0 + ,1190216 + ,783310 + ,0 + ,14416 + ,0 + ,361 + ,0 + ,49857 + ,0 + ,17998377 + ,467359 + ,1 + ,3369 + ,3369 + ,28 + ,28 + ,87161 + ,87161 + ,2440508 + ,688779 + ,1 + ,11819 + ,11819 + ,85 + ,85 + ,106113 + ,106113 + ,9019605 + ,608419 + ,1 + ,6620 + ,6620 + ,44 + ,44 + ,80570 + ,80570 + ,3545080 + ,696348 + ,1 + ,4519 + ,4519 + ,49 + ,49 + ,102129 + ,102129 + ,5004321 + ,597793 + ,0 + ,2220 + ,0 + ,22 + ,0 + ,301670 + ,0 + ,6636740 + ,821730 + ,0 + ,18562 + ,0 + ,155 + ,0 + ,102313 + ,0 + ,15858515 + ,377934 + ,0 + ,10327 + ,0 + ,91 + ,0 + ,88577 + ,0 + ,8060507 + ,651939 + ,1 + ,5336 + ,5336 + ,81 + ,81 + ,112477 + ,112477 + ,9110637 + ,697458 + ,1 + ,2365 + ,2365 + ,79 + ,79 + ,191778 + ,191778 + ,15150462 + ,700368 + ,0 + ,4069 + ,0 + ,145 + ,0 + ,79804 + ,0 + ,11571580 + ,225986 + ,0 + ,7710 + ,0 + ,816 + ,0 + ,128294 + ,0 + ,104687904 + ,348695 + ,0 + ,13718 + ,0 + ,61 + ,0 + ,96448 + ,0 + ,5883328 + ,373683 + ,0 + ,4525 + ,0 + ,226 + ,0 + ,93811 + ,0 + ,21201286 + ,501709 + ,0 + ,6869 + ,0 + ,105 + ,0 + ,117520 + ,0 + ,12339600 + ,413743 + ,0 + ,4628 + ,0 + ,62 + ,0 + ,69159 + ,0 + ,4287858 + ,379825 + ,1 + ,3653 + ,3653 + ,24 + ,24 + ,101792 + ,101792 + ,2443008 + ,336260 + ,1 + ,1265 + ,1265 + ,26 + ,26 + ,210568 + ,210568 + ,5474768 + ,636765 + ,1 + ,7489 + ,7489 + ,322 + ,322 + ,136996 + ,136996 + ,44112712 + ,481231 + ,0 + ,4901 + ,0 + ,84 + ,0 + ,121920 + ,0 + ,10241280 + ,469107 + ,0 + ,2284 + ,0 + ,33 + ,0 + ,76403 + ,0 + ,2521299 + ,211928 + ,1 + ,3160 + ,3160 + ,108 + ,108 + ,108094 + ,108094 + ,11674152 + ,563925 + ,1 + ,4150 + ,4150 + ,150 + ,150 + ,134759 + ,134759 + ,20213850 + ,511939 + ,1 + ,7285 + ,7285 + ,115 + ,115 + ,188873 + ,188873 + ,21720395 + ,521016 + ,1 + ,1134 + ,1134 + ,162 + ,162 + ,146216 + ,146216 + ,23686992 + ,543856 + ,1 + ,4658 + ,4658 + ,158 + ,158 + ,156608 + ,156608 + ,24744064 + ,329304 + ,0 + ,2384 + ,0 + ,97 + ,0 + ,61348 + ,0 + ,5950756 + ,423262 + ,0 + ,3748 + ,0 + ,9 + ,0 + ,50350 + ,0 + ,453150 + ,509665 + ,0 + ,5371 + ,0 + ,66 + ,0 + ,87720 + ,0 + ,5789520 + ,455881 + ,0 + ,1285 + ,0 + ,107 + ,0 + ,99489 + ,0 + ,10645323 + ,367772 + ,1 + ,9327 + ,9327 + ,101 + ,101 + ,87419 + ,87419 + ,8829319 + ,406339 + ,1 + ,5565 + ,5565 + ,47 + ,47 + ,94355 + ,94355 + ,4434685 + ,493408 + ,0 + ,1528 + ,0 + ,38 + ,0 + ,60326 + ,0 + ,2292388 + ,232942 + ,1 + ,3122 + ,3122 + ,34 + ,34 + ,94670 + ,94670 + ,3218780 + ,416002 + ,1 + ,7317 + ,7317 + ,84 + ,84 + ,82425 + ,82425 + ,6923700 + ,337430 + ,0 + ,2675 + ,0 + ,79 + ,0 + ,59017 + ,0 + ,4662343 + ,361517 + ,0 + ,13253 + ,0 + ,947 + ,0 + ,90829 + ,0 + ,86015063 + ,360962 + ,0 + ,880 + ,0 + ,74 + ,0 + ,80791 + ,0 + ,5978534 + ,235561 + ,1 + ,2053 + ,2053 + ,53 + ,53 + ,100423 + ,100423 + ,5322419 + ,408247 + ,0 + ,1424 + ,0 + ,94 + ,0 + ,131116 + ,0 + ,12324904 + ,450296 + ,1 + ,4036 + ,4036 + ,63 + ,63 + ,100269 + ,100269 + ,6316947 + ,418799 + ,1 + ,3045 + ,3045 + ,58 + ,58 + ,27330 + ,27330 + ,1585140 + ,247405 + ,0 + ,5119 + ,0 + ,49 + ,0 + ,39039 + ,0 + ,1912911 + ,378519 + ,0 + ,1431 + ,0 + ,34 + ,0 + ,106885 + ,0 + ,3634090 + ,326638 + ,0 + ,554 + ,0 + ,11 + ,0 + ,79285 + ,0 + ,872135 + ,328233 + ,0 + ,1975 + ,0 + ,35 + ,0 + ,118881 + ,0 + ,4160835 + ,386225 + ,1 + ,1286 + ,1286 + ,17 + ,17 + ,77623 + ,77623 + ,1319591 + ,283662 + ,0 + ,1012 + ,0 + ,47 + ,0 + ,114768 + ,0 + ,5394096 + ,370225 + ,0 + ,810 + ,0 + ,43 + ,0 + ,74015 + ,0 + ,3182645 + ,269236 + ,0 + ,1280 + ,0 + ,117 + ,0 + ,69465 + ,0 + ,8127405 + ,365732 + ,1 + ,666 + ,666 + ,171 + ,171 + ,117869 + ,117869 + ,20155599 + ,420383 + ,0 + ,1380 + ,0 + ,26 + ,0 + ,60982 + ,0 + ,1585532 + ,345811 + ,1 + ,4608 + ,4608 + ,73 + ,73 + ,90131 + ,90131 + ,6579563 + ,431809 + ,0 + ,876 + ,0 + ,59 + ,0 + ,138971 + ,0 + ,8199289 + ,418876 + ,0 + ,814 + ,0 + ,18 + ,0 + ,39625 + ,0 + ,713250 + ,297476 + ,0 + ,514 + ,0 + ,15 + ,0 + ,102725 + ,0 + ,1540875 + ,416776 + ,1 + ,5692 + ,5692 + ,72 + ,72 + ,64239 + ,64239 + ,4625208 + ,357257 + ,0 + ,3642 + ,0 + ,86 + ,0 + ,90262 + ,0 + ,7762532 + ,458343 + ,0 + ,540 + ,0 + ,14 + ,0 + ,103960 + ,0 + ,1455440 + ,388386 + ,0 + ,2099 + ,0 + ,64 + ,0 + ,106611 + ,0 + ,6823104 + ,358934 + ,0 + ,567 + ,0 + ,11 + ,0 + ,103345 + ,0 + ,1136795 + ,407560 + ,0 + ,2001 + ,0 + ,52 + ,0 + ,95551 + ,0 + ,4968652 + ,392558 + ,1 + ,2949 + ,2949 + ,41 + ,41 + ,82903 + ,82903 + ,3399023 + ,373177 + ,0 + ,2253 + ,0 + ,99 + ,0 + ,63593 + ,0 + ,6295707 + ,428370 + ,1 + ,6533 + ,6533 + ,75 + ,75 + ,126910 + ,126910 + ,9518250 + ,369419 + ,0 + ,1889 + ,0 + ,45 + ,0 + ,37527 + ,0 + ,1688715 + ,358649 + ,1 + ,3055 + ,3055 + ,43 + ,43 + ,60247 + ,60247 + ,2590621 + ,376641 + ,0 + ,272 + ,0 + ,8 + ,0 + ,112995 + ,0 + ,903960 + ,467427 + ,1 + ,1414 + ,1414 + ,198 + ,198 + ,70184 + ,70184 + ,13896432 + ,364885 + ,0 + ,2564 + ,0 + ,22 + ,0 + ,130140 + ,0 + ,2863080 + ,436230 + ,1 + ,1383 + ,1383 + ,11 + ,11 + ,73221 + ,73221 + ,805431 + ,329118 + ,1 + ,1261 + ,1261 + ,33 + ,33 + ,76114 + ,76114 + ,2511762 + ,317365 + ,0 + ,975 + ,0 + ,23 + ,0 + ,90534 + ,0 + ,2082282 + ,286849 + ,0 + ,3366 + ,0 + ,80 + ,0 + ,108479 + ,0 + ,8678320 + ,376685 + ,0 + ,576 + ,0 + ,18 + ,0 + ,113761 + ,0 + ,2047698 + ,407198 + ,0 + ,1306 + ,0 + ,28 + ,0 + ,68696 + ,0 + ,1923488 + ,377772 + ,0 + ,746 + ,0 + ,23 + ,0 + ,71561 + ,0 + ,1645903 + ,271483 + ,1 + ,3192 + ,3192 + ,60 + ,60 + ,59831 + ,59831 + ,3589860 + ,153661 + ,1 + ,2045 + ,2045 + ,20 + ,20 + ,97890 + ,97890 + ,1957800 + ,513294 + ,0 + ,5477 + ,0 + ,59 + ,0 + ,101481 + ,0 + ,5987379 + ,324881 + ,1 + ,1932 + ,1932 + ,36 + ,36 + ,72954 + ,72954 + ,2626344 + ,264512 + ,0 + ,936 + ,0 + ,30 + ,0 + ,67939 + ,0 + ,2038170 + ,420968 + ,1 + ,3437 + ,3437 + ,47 + ,47 + ,48022 + ,48022 + ,2257034 + ,129302 + ,0 + ,5131 + ,0 + ,71 + ,0 + ,86111 + ,0 + ,6113881 + ,191521 + ,1 + ,2397 + ,2397 + ,14 + ,14 + ,74020 + ,74020 + ,1036280 + ,268673 + ,1 + ,1389 + ,1389 + ,9 + ,9 + ,57530 + ,57530 + ,517770 + ,353179 + ,0 + ,1503 + ,0 + ,39 + ,0 + ,56364 + ,0 + ,2198196 + ,354624 + ,0 + ,402 + ,0 + ,26 + ,0 + ,84990 + ,0 + ,2209740 + ,363713 + ,0 + ,2239 + ,0 + ,21 + ,0 + ,88590 + ,0 + ,1860390 + ,456657 + ,1 + ,2234 + ,2234 + ,16 + ,16 + ,77200 + ,77200 + ,1235200 + ,211742 + ,0 + ,837 + ,0 + ,69 + ,0 + ,61262 + ,0 + ,4227078 + ,338381 + ,0 + ,10579 + ,0 + ,92 + ,0 + ,110309 + ,0 + ,10148428 + ,418530 + ,0 + ,875 + ,0 + ,14 + ,0 + ,67000 + ,0 + ,938000 + ,351483 + ,0 + ,1395 + ,0 + ,103 + ,0 + ,93099 + ,0 + ,9589197 + ,372928 + ,1 + ,1659 + ,1659 + ,29 + ,29 + ,107577 + ,107577 + ,3119733 + ,485538 + ,1 + ,2647 + ,2647 + ,37 + ,37 + ,62920 + ,62920 + ,2328040 + ,279268 + ,1 + ,3294 + ,3294 + ,23 + ,23 + ,75832 + ,75832 + ,1744136 + ,219060 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,94 + ,0 + ,7 + ,0 + ,60793 + ,0 + ,425551 + ,325314 + ,0 + ,422 + ,0 + ,28 + ,0 + ,57935 + ,0 + ,1622180 + ,322046 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,34 + ,0 + ,8 + ,0 + ,60630 + ,0 + ,485040 + ,325599 + ,0 + ,1558 + ,0 + ,63 + ,0 + ,55637 + ,0 + ,3505131 + ,377028 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,43 + ,0 + ,3 + ,0 + ,60887 + ,0 + ,182661 + ,323850 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,316 + ,0 + ,9 + ,0 + ,60505 + ,0 + ,544545 + ,331514 + ,0 + ,115 + ,0 + ,13 + ,0 + ,60945 + ,0 + ,792285 + ,325632 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,389 + ,0 + ,14 + ,0 + ,58990 + ,0 + ,825860 + ,322265 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,1002 + ,0 + ,15 + ,0 + ,56750 + ,0 + ,851250 + ,325906 + ,0 + ,36 + ,0 + ,3 + ,0 + ,60894 + ,0 + ,182682 + ,325985 + ,0 + ,460 + ,0 + ,15 + ,0 + ,63346 + ,0 + ,950190 + ,346145 + ,0 + ,309 + ,0 + ,11 + ,0 + ,56535 + ,0 + ,621885 + ,325898 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,9 + ,0 + ,6 + ,0 + ,60835 + ,0 + ,365010 + ,325356 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,14 + ,0 + ,1 + ,0 + ,61016 + ,0 + ,61016 + ,325930 + ,0 + ,520 + ,0 + ,10 + ,0 + ,58650 + ,0 + ,586500 + ,318020 + ,0 + ,1766 + ,0 + ,73 + ,0 + ,60438 + ,0 + ,4411974 + ,326389 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,458 + ,0 + ,11 + ,0 + ,58625 + ,0 + ,644875 + ,302925 + ,0 + ,20 + ,0 + ,3 + ,0 + ,60938 + ,0 + ,182814 + ,325540 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,98 + ,0 + ,2 + ,0 + ,61490 + ,0 + ,122980 + ,326736 + ,0 + ,405 + ,0 + ,7 + ,0 + ,60845 + ,0 + ,425915 + ,340580 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,483 + ,0 + ,27 + ,0 + ,60830 + ,0 + ,1642410 + ,331828 + ,0 + ,454 + ,0 + ,51 + ,0 + ,63261 + ,0 + ,3226311 + ,323299 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,757 + ,0 + ,19 + ,0 + ,45689 + ,0 + ,868091 + ,387722 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,36 + ,0 + ,4 + ,0 + ,61564 + ,0 + ,246256 + ,324598 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,203 + ,0 + ,9 + ,0 + ,61938 + ,0 + ,557442 + ,328726 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,90 + ,0 + ,8 + ,0 + ,60951 + ,0 + ,487608 + ,325043 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,71 + ,71 + ,1 + ,1 + ,60745 + ,60745 + ,60745 + ,325806 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,972 + ,0 + ,34 + ,0 + ,71642 + ,0 + ,2435828 + ,387732 + ,1 + ,531 + ,531 + ,10 + ,10 + ,71641 + ,71641 + ,716410 + ,349729 + ,0 + ,604 + ,0 + ,38 + ,0 + ,55792 + ,0 + ,2120096 + ,332202 + ,1 + ,283 + ,283 + ,10 + ,10 + ,71873 + ,71873 + ,718730 + ,305442 + ,1 + ,23 + ,23 + ,5 + ,5 + ,62555 + ,62555 + ,312775 + ,329537 + ,1 + ,638 + ,638 + ,14 + ,14 + ,60370 + ,60370 + ,845180 + ,327055 + ,1 + ,699 + ,699 + ,16 + ,16 + ,64873 + ,64873 + ,1037968 + ,356245 + ,0 + ,149 + ,0 + ,5 + ,0 + ,62041 + ,0 + ,310205 + ,328451 + ,0 + ,226 + ,0 + ,5 + ,0 + ,65745 + ,0 + ,328725 + ,307062 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,275 + ,0 + ,4 + ,0 + ,59500 + ,0 + ,238000 + ,331345 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,141 + ,0 + ,6 + ,0 + ,61630 + ,0 + ,369780 + ,331824 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,28 + ,0 + ,2 + ,0 + ,60890 + ,0 + ,121780 + ,325685 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,2566 + ,2566 + ,80 + ,80 + ,113521 + ,113521 + ,9081680 + ,404480 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,472 + ,472 + ,20 + ,20 + ,80045 + ,80045 + ,1600900 + ,318314 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,203 + ,203 + ,27 + ,27 + ,50804 + ,50804 + ,1371708 + ,311807 + ,1 + ,496 + ,496 + ,17 + ,17 + ,87390 + ,87390 + ,1485630 + ,337724 + ,1 + ,10 + ,10 + ,2 + ,2 + ,61656 + ,61656 + ,123312 + ,326431 + ,1 + ,63 + ,63 + ,4 + ,4 + ,65688 + ,65688 + ,262752 + ,327556 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,1136 + ,1136 + ,32 + ,32 + ,48522 + ,48522 + ,1552704 + ,356850 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,267 + ,0 + ,32 + ,0 + ,57640 + ,0 + ,1844480 + ,322741 + ,0 + ,474 + ,0 + ,20 + ,0 + ,61977 + ,0 + ,1239540 + ,310902 + ,0 + ,534 + ,0 + ,7 + ,0 + ,62620 + ,0 + ,438340 + ,324295 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,15 + ,0 + ,8 + ,0 + ,60831 + ,0 + ,486648 + ,326156 + ,0 + ,397 + ,0 + ,28 + ,0 + ,60646 + ,0 + ,1698088 + ,326960 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,1061 + ,0 + ,20 + ,0 + ,56225 + ,0 + ,1124500 + ,333411 + ,0 + ,288 + ,0 + ,4 + ,0 + ,60510 + ,0 + ,242040 + ,297761 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,3 + ,0 + ,2 + ,0 + ,60698 + ,0 + ,121396 + ,325536 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,20 + ,0 + ,2 + ,0 + ,60805 + ,0 + ,121610 + ,325762 + ,0 + ,278 + ,0 + ,26 + ,0 + ,61404 + ,0 + ,1596504 + ,327957 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,192 + ,0 + ,4 + ,0 + ,65276 + ,0 + ,261104 + ,318521 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,317 + ,0 + ,9 + ,0 + ,63915 + ,0 + ,575235 + ,319775 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,368 + ,0 + ,17 + ,0 + ,61686 + ,0 + ,1048662 + ,332128 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,2 + ,0 + ,1 + ,0 + ,60743 + ,0 + ,60743 + ,325486 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,53 + ,0 + ,6 + ,0 + ,60349 + ,0 + ,362094 + ,325838 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,94 + ,0 + ,3 + ,0 + ,61360 + ,0 + ,184080 + ,331767 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,24 + ,0 + ,8 + ,0 + ,59818 + ,0 + ,478544 + ,324523 + ,0 + ,2332 + ,0 + ,4 + ,0 + ,72680 + ,0 + ,290720 + ,339995 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,131 + ,0 + ,11 + ,0 + ,61808 + ,0 + ,679888 + ,319582 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,206 + ,0 + ,9 + ,0 + ,53110 + ,0 + ,477990 + ,307245 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,167 + ,0 + ,2 + ,0 + ,64245 + ,0 + ,128490 + ,317967 + ,0 + ,622 + ,0 + ,73 + ,0 + ,73007 + ,0 + ,5329511 + ,331488 + ,0 + ,885 + ,0 + ,85 + ,0 + ,82732 + ,0 + ,7032220 + ,335452 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,365 + ,0 + ,8 + ,0 + ,54820 + ,0 + ,438560 + ,334184 + ,0 + ,364 + ,0 + ,35 + ,0 + ,47705 + ,0 + ,1669675 + ,313213 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,226 + ,0 + ,12 + ,0 + ,72835 + ,0 + ,874020 + ,348678 + ,0 + ,307 + ,0 + ,15 + ,0 + ,58856 + ,0 + ,882840 + ,328727 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,188 + ,0 + ,11 + ,0 + ,77655 + ,0 + ,854205 + ,387978 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,138 + ,0 + ,6 + ,0 + ,69817 + ,0 + ,418902 + ,336704 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,125 + ,0 + ,12 + ,0 + ,60798 + ,0 + ,729576 + ,322076 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,282 + ,0 + ,30 + ,0 + ,62452 + ,0 + ,1873560 + ,334272 + ,0 + ,335 + ,0 + ,33 + ,0 + ,64175 + ,0 + ,2117775 + ,338197 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,813 + ,813 + ,82 + ,82 + ,67440 + ,67440 + ,5530080 + ,321024 + ,0 + ,176 + ,0 + ,28 + ,0 + ,68136 + ,0 + ,1907808 + ,322145 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,249 + ,0 + ,72 + ,0 + ,56726 + ,0 + ,4084272 + ,323351 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,333 + ,0 + ,13 + ,0 + ,70811 + ,0 + ,920543 + ,327748 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,30 + ,0 + ,4 + ,0 + ,62045 + ,0 + ,248180 + ,328157 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,249 + ,0 + ,62 + ,0 + ,54323 + ,0 + ,3368026 + ,311594 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,165 + ,0 + ,24 + ,0 + ,62841 + ,0 + ,1508184 + ,335962 + ,0 + ,453 + ,0 + ,21 + ,0 + ,81125 + ,0 + ,1703625 + ,372426 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,53 + ,0 + ,14 + ,0 + ,59506 + ,0 + ,833084 + ,319844 + ,1 + ,382 + ,382 + ,21 + ,21 + ,59365 + ,59365 + ,1246665 + ,355822 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,30 + ,30 + ,4 + ,4 + ,60798 + ,60798 + ,243192 + ,324047 + ,0 + ,290 + ,0 + ,2 + ,0 + ,58790 + ,0 + ,117580 + ,311464 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,366 + ,0 + ,53 + ,0 + ,61808 + ,0 + ,3275824 + ,353417 + ,0 + ,2 + ,0 + ,9 + ,0 + ,60735 + ,0 + ,546615 + ,325590 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,209 + ,209 + ,13 + ,13 + ,64016 + ,64016 + ,832208 + ,328576 + ,0 + ,384 + ,0 + ,22 + ,0 + ,54683 + ,0 + ,1203026 + ,326126 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,365 + ,0 + ,83 + ,0 + ,87192 + ,0 + ,7236936 + ,369376 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,49 + ,49 + ,8 + ,8 + ,64107 + ,64107 + ,512856 + ,332013 + ,0 + ,3 + ,0 + ,4 + ,0 + ,60761 + ,0 + ,243044 + ,325871 + ,0 + ,133 + ,0 + ,14 + ,0 + ,65990 + ,0 + ,923860 + ,342165 + ,0 + ,32 + ,0 + ,1 + ,0 + ,59988 + ,0 + ,59988 + ,324967 + ,0 + ,368 + ,0 + ,17 + ,0 + ,61167 + ,0 + ,1039839 + ,314832 + ,0 + ,1 + ,0 + ,6 + ,0 + ,60719 + ,0 + ,364314 + ,325557 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,22 + ,0 + ,2 + ,0 + ,60722 + ,0 + ,121444 + ,322649 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,96 + ,0 + ,5 + ,0 + ,60379 + ,0 + ,301895 + ,324598 + ,0 + ,1 + ,0 + ,2 + ,0 + ,60727 + ,0 + ,121454 + ,325567 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,81 + ,0 + ,7 + ,0 + ,60925 + ,0 + ,426475 + ,324005 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,26 + ,0 + ,1 + ,0 + ,60896 + ,0 + ,60896 + ,325748 + ,0 + ,125 + ,0 + ,13 + ,0 + ,59734 + ,0 + ,776542 + ,323385 + ,0 + ,304 + ,0 + ,15 + ,0 + ,62969 + ,0 + ,944535 + ,315409 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,119 + ,0 + ,6 + ,0 + ,59118 + ,0 + ,354708 + ,312275 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,312 + ,0 + ,14 + ,0 + ,58598 + ,0 + ,820372 + ,320576 + ,0 + ,60 + ,0 + ,10 + ,0 + ,61124 + ,0 + ,611240 + ,325246 + ,0 + ,587 + ,0 + ,12 + ,0 + ,59595 + ,0 + ,715140 + ,332961 + ,0 + ,135 + ,0 + ,2 + ,0 + ,62065 + ,0 + ,124130 + ,323010 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,514 + ,0 + ,52 + ,0 + ,78780 + ,0 + ,4096560 + ,345253 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,1 + ,0 + ,4 + ,0 + ,60722 + ,0 + ,242888 + ,325559 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,1 + ,58 + ,58 + ,3 + ,3 + ,61600 + ,61600 + ,184800 + ,319634 + ,0 + ,180 + ,0 + ,11 + ,0 + ,59635 + ,0 + ,655985 + ,319951 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,448 + ,0 + ,40 + ,0 + ,59781 + ,0 + ,2391240 + ,318519 + ,0 + ,227 + ,0 + ,9 + ,0 + ,76644 + ,0 + ,689796 + ,343222 + ,0 + ,174 + ,0 + ,1 + ,0 + ,64820 + ,0 + ,64820 + ,317234 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,121 + ,0 + ,24 + ,0 + ,56178 + ,0 + ,1348272 + ,314025 + ,0 + ,607 + ,0 + ,11 + ,0 + ,60436 + ,0 + ,664796 + ,320249 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,530 + ,0 + ,60 + ,0 + ,73433 + ,0 + ,4405980 + ,349365 + ,0 + ,571 + ,0 + ,80 + ,0 + ,41477 + ,0 + ,3318160 + ,289197 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,0 + ,325560 + ,0 + ,78 + ,0 + ,16 + ,0 + ,62700 + ,0 + ,1003200 + ,329245 + ,0 + ,2489 + ,0 + ,40 + ,0 + ,67804 + ,0 + ,2712160 + ,240869 + ,0 + ,131 + ,0 + ,6 + ,0 + ,59661 + ,0 + ,357966 + ,327182 + ,0 + ,923 + ,0 + ,8 + ,0 + ,58620 + ,0 + ,468960 + ,322876 + ,0 + ,72 + ,0 + ,3 + ,0 + ,60398 + ,0 + ,181194 + ,323117 + ,0 + ,572 + ,0 + ,16 + ,0 + ,58580 + ,0 + ,937280 + ,306351 + ,0 + ,397 + ,0 + ,10 + ,0 + ,62710 + ,0 + ,627100 + ,335137 + ,0 + ,450 + ,0 + ,8 + ,0 + ,59325 + ,0 + ,474600 + ,308271 + ,0 + ,622 + ,0 + ,7 + ,0 + ,60950 + ,0 + ,426650 + ,301731 + ,0 + ,694 + ,0 + ,8 + ,0 + ,68060 + ,0 + ,544480 + ,382409 + ,1 + ,3425 + ,3425 + ,12 + ,12 + ,83620 + ,83620 + ,1003440 + ,279230 + ,0 + ,562 + ,0 + ,13 + ,0 + ,58456 + ,0 + ,759928 + ,298731 + ,0 + ,4917 + ,0 + ,42 + ,0 + ,52811 + ,0 + ,2218062 + ,243650 + ,1 + ,1442 + ,1442 + ,118 + ,118 + ,121173 + ,121173 + ,14298414 + ,532682 + ,0 + ,529 + ,0 + ,9 + ,0 + ,63870 + ,0 + ,574830 + ,319771 + ,1 + ,2126 + ,2126 + ,138 + ,138 + ,21001 + ,21001 + ,2898138 + ,171493 + ,0 + ,1061 + ,0 + ,5 + ,0 + ,70415 + ,0 + ,352075 + ,347262 + ,0 + ,776 + ,0 + ,9 + ,0 + ,64230 + ,0 + ,578070 + ,343945 + ,0 + ,611 + ,0 + ,8 + ,0 + ,59190 + ,0 + ,473520 + ,311874 + ,1 + ,1526 + ,1526 + ,25 + ,25 + ,69351 + ,69351 + ,1733775 + ,302211 + ,0 + ,592 + ,0 + ,7 + ,0 + ,64270 + ,0 + ,449890 + ,316708 + ,0 + ,1182 + ,0 + ,13 + ,0 + ,70694 + ,0 + ,919022 + ,333463 + ,0 + ,621 + ,0 + ,16 + ,0 + ,68005 + ,0 + ,1088080 + ,344282 + ,0 + ,989 + ,0 + ,11 + ,0 + ,58930 + ,0 + ,648230 + ,319635 + ,0 + ,438 + ,0 + ,11 + ,0 + ,58320 + ,0 + ,641520 + ,301186 + ,0 + ,726 + ,0 + ,3 + ,0 + ,69980 + ,0 + ,209940 + ,300381 + ,0 + ,1303 + ,0 + ,61 + ,0 + ,69863 + ,0 + ,4261643 + ,318765 + ,1 + ,6341 + ,6341 + ,24 + ,24 + ,63255 + ,63255 + ,1518120 + ,286146 + ,1 + ,1164 + ,1164 + ,17 + ,17 + ,57320 + ,57320 + ,974440 + ,306844 + ,1 + ,3310 + ,3310 + ,33 + ,33 + ,75230 + ,75230 + ,2482590 + ,307705 + ,0 + ,1366 + ,0 + ,7 + ,0 + ,79420 + ,0 + ,555940 + ,312448 + ,0 + ,965 + ,0 + ,3 + ,0 + ,73490 + ,0 + ,220470 + ,299715 + ,0 + ,3256 + ,0 + ,66 + ,0 + ,35250 + ,0 + ,2326500 + ,373399 + ,1 + ,1135 + ,1135 + ,17 + ,17 + ,62285 + ,62285 + ,1058845 + ,299446 + ,0 + ,1270 + ,0 + ,26 + ,0 + ,69206 + ,0 + ,1799356 + ,325586 + ,0 + ,661 + ,0 + ,3 + ,0 + ,65920 + ,0 + ,197760 + ,291221 + ,0 + ,1013 + ,0 + ,2 + ,0 + ,69770 + ,0 + ,139540 + ,261173 + ,0 + ,2844 + ,0 + ,67 + ,0 + ,72683 + ,0 + ,4869761 + ,255027 + ,1 + ,11528 + ,11528 + ,70 + ,70 + ,-14545 + ,-14545 + ,-1018150 + ,-78375 + ,0 + ,6526 + ,0 + ,26 + ,0 + ,55830 + ,0 + ,1451580 + ,-58143 + ,0 + ,2264 + ,0 + ,24 + ,0 + ,55174 + ,0 + ,1324176 + ,227033 + ,1 + ,4461 + ,4461 + ,94 + ,94 + ,67038 + ,67038 + ,6301572 + ,235098 + ,0 + ,3999 + ,0 + ,30 + ,0 + ,51252 + ,0 + ,1537560 + ,21267 + ,0 + ,35624 + ,0 + ,223 + ,0 + ,157278 + ,0 + ,35072994 + ,238675 + ,0 + ,9252 + ,0 + ,48 + ,0 + ,79510 + ,0 + ,3816480 + ,197687 + ,0 + ,15236 + ,0 + ,90 + ,0 + ,77440 + ,0 + ,6969600 + ,418341 + ,0 + ,18073 + ,0 + ,180 + ,0 + ,27284 + ,0 + ,4911120 + ,-297706) + ,dim=c(9 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'GrCosts' + ,'Trades' + ,'GrTrades' + ,'Dividends' + ,'GrDiv' + ,'TrDiv' + ,'Wealth ') + ,1:431)) > y <- array(NA,dim=c(9,431),dimnames=list(c('Group','Costs','GrCosts','Trades','GrTrades','Dividends','GrDiv','TrDiv','Wealth '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '9' > #'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\r Group Costs GrCosts Trades GrTrades Dividends GrDiv TrDiv 1 6282929 1 162556 162556 1081 1081 213118 213118 230380558 2 4324047 1 29790 29790 309 309 81767 81767 25266003 3 4108272 1 87550 87550 458 458 153198 153198 70164684 4 -1212617 0 84738 0 588 0 -26007 0 -15292116 5 1485329 1 54660 54660 299 299 126942 126942 37955658 6 1779876 1 42634 42634 156 156 157214 157214 24525384 7 1367203 0 40949 0 481 0 129352 0 62218312 8 2519076 1 42312 42312 323 323 234817 234817 75845891 9 912684 1 37704 37704 452 452 60448 60448 27322496 10 1443586 1 16275 16275 109 109 47818 47818 5212162 11 1220017 0 25830 0 115 0 245546 0 28237790 12 984885 0 12679 0 110 0 48020 0 5282200 13 1457425 1 18014 18014 239 239 -1710 -1710 -408690 14 -572920 0 43556 0 247 0 32648 0 8064056 15 929144 1 24524 24524 497 497 95350 95350 47388950 16 1151176 0 6532 0 103 0 151352 0 15589256 17 790090 0 7123 0 109 0 288170 0 31410530 18 774497 1 20813 20813 502 502 114337 114337 57397174 19 990576 1 37597 37597 248 248 37884 37884 9395232 20 454195 0 17821 0 373 0 122844 0 45820812 21 876607 1 12988 12988 119 119 82340 82340 9798460 22 711969 1 22330 22330 84 84 79801 79801 6703284 23 702380 0 13326 0 102 0 165548 0 16885896 24 264449 0 16189 0 295 0 116384 0 34333280 25 450033 0 7146 0 105 0 134028 0 14072940 26 541063 0 15824 0 64 0 63838 0 4085632 27 588864 1 26088 26088 267 267 74996 74996 20023932 28 -37216 0 11326 0 129 0 31080 0 4009320 29 783310 0 8568 0 37 0 32168 0 1190216 30 467359 0 14416 0 361 0 49857 0 17998377 31 688779 1 3369 3369 28 28 87161 87161 2440508 32 608419 1 11819 11819 85 85 106113 106113 9019605 33 696348 1 6620 6620 44 44 80570 80570 3545080 34 597793 1 4519 4519 49 49 102129 102129 5004321 35 821730 0 2220 0 22 0 301670 0 6636740 36 377934 0 18562 0 155 0 102313 0 15858515 37 651939 0 10327 0 91 0 88577 0 8060507 38 697458 1 5336 5336 81 81 112477 112477 9110637 39 700368 1 2365 2365 79 79 191778 191778 15150462 40 225986 0 4069 0 145 0 79804 0 11571580 41 348695 0 7710 0 816 0 128294 0 104687904 42 373683 0 13718 0 61 0 96448 0 5883328 43 501709 0 4525 0 226 0 93811 0 21201286 44 413743 0 6869 0 105 0 117520 0 12339600 45 379825 0 4628 0 62 0 69159 0 4287858 46 336260 1 3653 3653 24 24 101792 101792 2443008 47 636765 1 1265 1265 26 26 210568 210568 5474768 48 481231 1 7489 7489 322 322 136996 136996 44112712 49 469107 0 4901 0 84 0 121920 0 10241280 50 211928 0 2284 0 33 0 76403 0 2521299 51 563925 1 3160 3160 108 108 108094 108094 11674152 52 511939 1 4150 4150 150 150 134759 134759 20213850 53 521016 1 7285 7285 115 115 188873 188873 21720395 54 543856 1 1134 1134 162 162 146216 146216 23686992 55 329304 1 4658 4658 158 158 156608 156608 24744064 56 423262 0 2384 0 97 0 61348 0 5950756 57 509665 0 3748 0 9 0 50350 0 453150 58 455881 0 5371 0 66 0 87720 0 5789520 59 367772 0 1285 0 107 0 99489 0 10645323 60 406339 1 9327 9327 101 101 87419 87419 8829319 61 493408 1 5565 5565 47 47 94355 94355 4434685 62 232942 0 1528 0 38 0 60326 0 2292388 63 416002 1 3122 3122 34 34 94670 94670 3218780 64 337430 1 7317 7317 84 84 82425 82425 6923700 65 361517 0 2675 0 79 0 59017 0 4662343 66 360962 0 13253 0 947 0 90829 0 86015063 67 235561 0 880 0 74 0 80791 0 5978534 68 408247 1 2053 2053 53 53 100423 100423 5322419 69 450296 0 1424 0 94 0 131116 0 12324904 70 418799 1 4036 4036 63 63 100269 100269 6316947 71 247405 1 3045 3045 58 58 27330 27330 1585140 72 378519 0 5119 0 49 0 39039 0 1912911 73 326638 0 1431 0 34 0 106885 0 3634090 74 328233 0 554 0 11 0 79285 0 872135 75 386225 0 1975 0 35 0 118881 0 4160835 76 283662 1 1286 1286 17 17 77623 77623 1319591 77 370225 0 1012 0 47 0 114768 0 5394096 78 269236 0 810 0 43 0 74015 0 3182645 79 365732 0 1280 0 117 0 69465 0 8127405 80 420383 1 666 666 171 171 117869 117869 20155599 81 345811 0 1380 0 26 0 60982 0 1585532 82 431809 1 4608 4608 73 73 90131 90131 6579563 83 418876 0 876 0 59 0 138971 0 8199289 84 297476 0 814 0 18 0 39625 0 713250 85 416776 0 514 0 15 0 102725 0 1540875 86 357257 1 5692 5692 72 72 64239 64239 4625208 87 458343 0 3642 0 86 0 90262 0 7762532 88 388386 0 540 0 14 0 103960 0 1455440 89 358934 0 2099 0 64 0 106611 0 6823104 90 407560 0 567 0 11 0 103345 0 1136795 91 392558 0 2001 0 52 0 95551 0 4968652 92 373177 1 2949 2949 41 41 82903 82903 3399023 93 428370 0 2253 0 99 0 63593 0 6295707 94 369419 1 6533 6533 75 75 126910 126910 9518250 95 358649 0 1889 0 45 0 37527 0 1688715 96 376641 1 3055 3055 43 43 60247 60247 2590621 97 467427 0 272 0 8 0 112995 0 903960 98 364885 1 1414 1414 198 198 70184 70184 13896432 99 436230 0 2564 0 22 0 130140 0 2863080 100 329118 1 1383 1383 11 11 73221 73221 805431 101 317365 1 1261 1261 33 33 76114 76114 2511762 102 286849 0 975 0 23 0 90534 0 2082282 103 376685 0 3366 0 80 0 108479 0 8678320 104 407198 0 576 0 18 0 113761 0 2047698 105 377772 0 1306 0 28 0 68696 0 1923488 106 271483 0 746 0 23 0 71561 0 1645903 107 153661 1 3192 3192 60 60 59831 59831 3589860 108 513294 1 2045 2045 20 20 97890 97890 1957800 109 324881 0 5477 0 59 0 101481 0 5987379 110 264512 1 1932 1932 36 36 72954 72954 2626344 111 420968 0 936 0 30 0 67939 0 2038170 112 129302 1 3437 3437 47 47 48022 48022 2257034 113 191521 0 5131 0 71 0 86111 0 6113881 114 268673 1 2397 2397 14 14 74020 74020 1036280 115 353179 1 1389 1389 9 9 57530 57530 517770 116 354624 0 1503 0 39 0 56364 0 2198196 117 363713 0 402 0 26 0 84990 0 2209740 118 456657 0 2239 0 21 0 88590 0 1860390 119 211742 1 2234 2234 16 16 77200 77200 1235200 120 338381 0 837 0 69 0 61262 0 4227078 121 418530 0 10579 0 92 0 110309 0 10148428 122 351483 0 875 0 14 0 67000 0 938000 123 372928 0 1395 0 103 0 93099 0 9589197 124 485538 1 1659 1659 29 29 107577 107577 3119733 125 279268 1 2647 2647 37 37 62920 62920 2328040 126 219060 1 3294 3294 23 23 75832 75832 1744136 127 325560 0 0 0 0 0 60720 0 0 128 325314 0 94 0 7 0 60793 0 425551 129 322046 0 422 0 28 0 57935 0 1622180 130 325560 0 0 0 0 0 60720 0 0 131 325599 0 34 0 8 0 60630 0 485040 132 377028 0 1558 0 63 0 55637 0 3505131 133 325560 0 0 0 0 0 60720 0 0 134 323850 0 43 0 3 0 60887 0 182661 135 325560 0 0 0 0 0 60720 0 0 136 331514 0 316 0 9 0 60505 0 544545 137 325632 0 115 0 13 0 60945 0 792285 138 325560 0 0 0 0 0 60720 0 0 139 325560 0 0 0 0 0 60720 0 0 140 325560 0 0 0 0 0 60720 0 0 141 322265 0 389 0 14 0 58990 0 825860 142 325560 0 0 0 0 0 60720 0 0 143 325906 0 1002 0 15 0 56750 0 851250 144 325985 0 36 0 3 0 60894 0 182682 145 346145 0 460 0 15 0 63346 0 950190 146 325898 0 309 0 11 0 56535 0 621885 147 325560 0 0 0 0 0 60720 0 0 148 325356 0 9 0 6 0 60835 0 365010 149 325560 0 0 0 0 0 60720 0 0 150 325930 0 14 0 1 0 61016 0 61016 151 318020 0 520 0 10 0 58650 0 586500 152 326389 0 1766 0 73 0 60438 0 4411974 153 325560 0 0 0 0 0 60720 0 0 154 302925 0 458 0 11 0 58625 0 644875 155 325540 0 20 0 3 0 60938 0 182814 156 325560 0 0 0 0 0 60720 0 0 157 325560 0 0 0 0 0 60720 0 0 158 326736 0 98 0 2 0 61490 0 122980 159 340580 0 405 0 7 0 60845 0 425915 160 325560 0 0 0 0 0 60720 0 0 161 325560 0 0 0 0 0 60720 0 0 162 325560 0 0 0 0 0 60720 0 0 163 325560 0 0 0 0 0 60720 0 0 164 331828 0 483 0 27 0 60830 0 1642410 165 323299 0 454 0 51 0 63261 0 3226311 166 325560 0 0 0 0 0 60720 0 0 167 325560 0 0 0 0 0 60720 0 0 168 387722 0 757 0 19 0 45689 0 868091 169 325560 0 0 0 0 0 60720 0 0 170 325560 0 0 0 0 0 60720 0 0 171 325560 0 0 0 0 0 60720 0 0 172 324598 0 36 0 4 0 61564 0 246256 173 325560 0 0 0 0 0 60720 0 0 174 328726 0 203 0 9 0 61938 0 557442 175 325560 0 0 0 0 0 60720 0 0 176 325043 0 90 0 8 0 60951 0 487608 177 325560 0 0 0 0 0 60720 0 0 178 325806 1 71 71 1 1 60745 60745 60745 179 325560 0 0 0 0 0 60720 0 0 180 325560 0 0 0 0 0 60720 0 0 181 387732 0 972 0 34 0 71642 0 2435828 182 349729 1 531 531 10 10 71641 71641 716410 183 332202 0 604 0 38 0 55792 0 2120096 184 305442 1 283 283 10 10 71873 71873 718730 185 329537 1 23 23 5 5 62555 62555 312775 186 327055 1 638 638 14 14 60370 60370 845180 187 356245 1 699 699 16 16 64873 64873 1037968 188 328451 0 149 0 5 0 62041 0 310205 189 307062 0 226 0 5 0 65745 0 328725 190 325560 0 0 0 0 0 60720 0 0 191 331345 0 275 0 4 0 59500 0 238000 192 325560 0 0 0 0 0 60720 0 0 193 331824 0 141 0 6 0 61630 0 369780 194 325560 0 0 0 0 0 60720 0 0 195 325685 0 28 0 2 0 60890 0 121780 196 325560 1 0 0 0 0 60720 60720 0 197 404480 1 2566 2566 80 80 113521 113521 9081680 198 325560 1 0 0 0 0 60720 60720 0 199 325560 1 0 0 0 0 60720 60720 0 200 318314 1 472 472 20 20 80045 80045 1600900 201 325560 1 0 0 0 0 60720 60720 0 202 325560 1 0 0 0 0 60720 60720 0 203 325560 1 0 0 0 0 60720 60720 0 204 311807 1 203 203 27 27 50804 50804 1371708 205 337724 1 496 496 17 17 87390 87390 1485630 206 326431 1 10 10 2 2 61656 61656 123312 207 327556 1 63 63 4 4 65688 65688 262752 208 325560 1 0 0 0 0 60720 60720 0 209 356850 1 1136 1136 32 32 48522 48522 1552704 210 325560 1 0 0 0 0 60720 60720 0 211 325560 0 0 0 0 0 60720 0 0 212 325560 0 0 0 0 0 60720 0 0 213 322741 0 267 0 32 0 57640 0 1844480 214 310902 0 474 0 20 0 61977 0 1239540 215 324295 0 534 0 7 0 62620 0 438340 216 325560 0 0 0 0 0 60720 0 0 217 326156 0 15 0 8 0 60831 0 486648 218 326960 0 397 0 28 0 60646 0 1698088 219 325560 0 0 0 0 0 60720 0 0 220 333411 0 1061 0 20 0 56225 0 1124500 221 297761 0 288 0 4 0 60510 0 242040 222 325560 0 0 0 0 0 60720 0 0 223 325536 0 3 0 2 0 60698 0 121396 224 325560 0 0 0 0 0 60720 0 0 225 325762 0 20 0 2 0 60805 0 121610 226 327957 0 278 0 26 0 61404 0 1596504 227 325560 0 0 0 0 0 60720 0 0 228 325560 0 0 0 0 0 60720 0 0 229 318521 0 192 0 4 0 65276 0 261104 230 325560 0 0 0 0 0 60720 0 0 231 319775 0 317 0 9 0 63915 0 575235 232 325560 0 0 0 0 0 60720 0 0 233 325560 0 0 0 0 0 60720 0 0 234 332128 0 368 0 17 0 61686 0 1048662 235 325560 0 0 0 0 0 60720 0 0 236 325486 0 2 0 1 0 60743 0 60743 237 325560 0 0 0 0 0 60720 0 0 238 325838 0 53 0 6 0 60349 0 362094 239 325560 0 0 0 0 0 60720 0 0 240 325560 0 0 0 0 0 60720 0 0 241 325560 0 0 0 0 0 60720 0 0 242 331767 0 94 0 3 0 61360 0 184080 243 325560 0 0 0 0 0 60720 0 0 244 324523 0 24 0 8 0 59818 0 478544 245 339995 0 2332 0 4 0 72680 0 290720 246 325560 1 0 0 0 0 60720 60720 0 247 325560 1 0 0 0 0 60720 60720 0 248 319582 0 131 0 11 0 61808 0 679888 249 325560 1 0 0 0 0 60720 60720 0 250 325560 1 0 0 0 0 60720 60720 0 251 307245 0 206 0 9 0 53110 0 477990 252 325560 1 0 0 0 0 60720 60720 0 253 317967 0 167 0 2 0 64245 0 128490 254 331488 0 622 0 73 0 73007 0 5329511 255 335452 0 885 0 85 0 82732 0 7032220 256 325560 0 0 0 0 0 60720 0 0 257 334184 0 365 0 8 0 54820 0 438560 258 313213 0 364 0 35 0 47705 0 1669675 259 325560 1 0 0 0 0 60720 60720 0 260 325560 0 0 0 0 0 60720 0 0 261 325560 1 0 0 0 0 60720 60720 0 262 325560 1 0 0 0 0 60720 60720 0 263 348678 0 226 0 12 0 72835 0 874020 264 328727 0 307 0 15 0 58856 0 882840 265 325560 1 0 0 0 0 60720 60720 0 266 325560 0 0 0 0 0 60720 0 0 267 325560 1 0 0 0 0 60720 60720 0 268 387978 0 188 0 11 0 77655 0 854205 269 325560 1 0 0 0 0 60720 60720 0 270 336704 0 138 0 6 0 69817 0 418902 271 325560 1 0 0 0 0 60720 60720 0 272 325560 1 0 0 0 0 60720 60720 0 273 325560 1 0 0 0 0 60720 60720 0 274 322076 0 125 0 12 0 60798 0 729576 275 325560 0 0 0 0 0 60720 0 0 276 334272 0 282 0 30 0 62452 0 1873560 277 338197 0 335 0 33 0 64175 0 2117775 278 325560 0 0 0 0 0 60720 0 0 279 321024 1 813 813 82 82 67440 67440 5530080 280 322145 0 176 0 28 0 68136 0 1907808 281 325560 1 0 0 0 0 60720 60720 0 282 325560 0 0 0 0 0 60720 0 0 283 323351 0 249 0 72 0 56726 0 4084272 284 325560 0 0 0 0 0 60720 0 0 285 327748 0 333 0 13 0 70811 0 920543 286 325560 0 0 0 0 0 60720 0 0 287 325560 1 0 0 0 0 60720 60720 0 288 328157 0 30 0 4 0 62045 0 248180 289 325560 1 0 0 0 0 60720 60720 0 290 311594 0 249 0 62 0 54323 0 3368026 291 325560 1 0 0 0 0 60720 60720 0 292 335962 0 165 0 24 0 62841 0 1508184 293 372426 0 453 0 21 0 81125 0 1703625 294 325560 0 0 0 0 0 60720 0 0 295 319844 0 53 0 14 0 59506 0 833084 296 355822 1 382 382 21 21 59365 59365 1246665 297 325560 0 0 0 0 0 60720 0 0 298 325560 1 0 0 0 0 60720 60720 0 299 325560 0 0 0 0 0 60720 0 0 300 325560 0 0 0 0 0 60720 0 0 301 324047 1 30 30 4 4 60798 60798 243192 302 311464 0 290 0 2 0 58790 0 117580 303 325560 0 0 0 0 0 60720 0 0 304 325560 1 0 0 0 0 60720 60720 0 305 353417 0 366 0 53 0 61808 0 3275824 306 325590 0 2 0 9 0 60735 0 546615 307 325560 1 0 0 0 0 60720 60720 0 308 328576 1 209 209 13 13 64016 64016 832208 309 326126 0 384 0 22 0 54683 0 1203026 310 325560 1 0 0 0 0 60720 60720 0 311 325560 1 0 0 0 0 60720 60720 0 312 369376 0 365 0 83 0 87192 0 7236936 313 325560 1 0 0 0 0 60720 60720 0 314 332013 1 49 49 8 8 64107 64107 512856 315 325871 0 3 0 4 0 60761 0 243044 316 342165 0 133 0 14 0 65990 0 923860 317 324967 0 32 0 1 0 59988 0 59988 318 314832 0 368 0 17 0 61167 0 1039839 319 325557 0 1 0 6 0 60719 0 364314 320 325560 0 0 0 0 0 60720 0 0 321 325560 0 0 0 0 0 60720 0 0 322 325560 0 0 0 0 0 60720 0 0 323 325560 0 0 0 0 0 60720 0 0 324 325560 1 0 0 0 0 60720 60720 0 325 325560 0 0 0 0 0 60720 0 0 326 322649 0 22 0 2 0 60722 0 121444 327 325560 1 0 0 0 0 60720 60720 0 328 325560 1 0 0 0 0 60720 60720 0 329 325560 0 0 0 0 0 60720 0 0 330 325560 0 0 0 0 0 60720 0 0 331 325560 0 0 0 0 0 60720 0 0 332 325560 0 0 0 0 0 60720 0 0 333 325560 0 0 0 0 0 60720 0 0 334 324598 0 96 0 5 0 60379 0 301895 335 325567 0 1 0 2 0 60727 0 121454 336 325560 0 0 0 0 0 60720 0 0 337 324005 0 81 0 7 0 60925 0 426475 338 325560 0 0 0 0 0 60720 0 0 339 325748 0 26 0 1 0 60896 0 60896 340 323385 0 125 0 13 0 59734 0 776542 341 315409 0 304 0 15 0 62969 0 944535 342 325560 0 0 0 0 0 60720 0 0 343 325560 0 0 0 0 0 60720 0 0 344 325560 0 0 0 0 0 60720 0 0 345 325560 0 0 0 0 0 60720 0 0 346 325560 0 0 0 0 0 60720 0 0 347 312275 0 119 0 6 0 59118 0 354708 348 325560 0 0 0 0 0 60720 0 0 349 325560 0 0 0 0 0 60720 0 0 350 325560 0 0 0 0 0 60720 0 0 351 320576 0 312 0 14 0 58598 0 820372 352 325246 0 60 0 10 0 61124 0 611240 353 332961 0 587 0 12 0 59595 0 715140 354 323010 0 135 0 2 0 62065 0 124130 355 325560 0 0 0 0 0 60720 0 0 356 325560 0 0 0 0 0 60720 0 0 357 345253 0 514 0 52 0 78780 0 4096560 358 325560 0 0 0 0 0 60720 0 0 359 325560 0 0 0 0 0 60720 0 0 360 325560 0 0 0 0 0 60720 0 0 361 325559 0 1 0 4 0 60722 0 242888 362 325560 0 0 0 0 0 60720 0 0 363 325560 0 0 0 0 0 60720 0 0 364 319634 1 58 58 3 3 61600 61600 184800 365 319951 0 180 0 11 0 59635 0 655985 366 325560 0 0 0 0 0 60720 0 0 367 325560 0 0 0 0 0 60720 0 0 368 325560 0 0 0 0 0 60720 0 0 369 325560 0 0 0 0 0 60720 0 0 370 325560 0 0 0 0 0 60720 0 0 371 325560 0 0 0 0 0 60720 0 0 372 318519 0 448 0 40 0 59781 0 2391240 373 343222 0 227 0 9 0 76644 0 689796 374 317234 0 174 0 1 0 64820 0 64820 375 325560 0 0 0 0 0 60720 0 0 376 325560 0 0 0 0 0 60720 0 0 377 314025 0 121 0 24 0 56178 0 1348272 378 320249 0 607 0 11 0 60436 0 664796 379 325560 0 0 0 0 0 60720 0 0 380 325560 0 0 0 0 0 60720 0 0 381 325560 0 0 0 0 0 60720 0 0 382 349365 0 530 0 60 0 73433 0 4405980 383 289197 0 571 0 80 0 41477 0 3318160 384 325560 0 0 0 0 0 60720 0 0 385 329245 0 78 0 16 0 62700 0 1003200 386 240869 0 2489 0 40 0 67804 0 2712160 387 327182 0 131 0 6 0 59661 0 357966 388 322876 0 923 0 8 0 58620 0 468960 389 323117 0 72 0 3 0 60398 0 181194 390 306351 0 572 0 16 0 58580 0 937280 391 335137 0 397 0 10 0 62710 0 627100 392 308271 0 450 0 8 0 59325 0 474600 393 301731 0 622 0 7 0 60950 0 426650 394 382409 0 694 0 8 0 68060 0 544480 395 279230 1 3425 3425 12 12 83620 83620 1003440 396 298731 0 562 0 13 0 58456 0 759928 397 243650 0 4917 0 42 0 52811 0 2218062 398 532682 1 1442 1442 118 118 121173 121173 14298414 399 319771 0 529 0 9 0 63870 0 574830 400 171493 1 2126 2126 138 138 21001 21001 2898138 401 347262 0 1061 0 5 0 70415 0 352075 402 343945 0 776 0 9 0 64230 0 578070 403 311874 0 611 0 8 0 59190 0 473520 404 302211 1 1526 1526 25 25 69351 69351 1733775 405 316708 0 592 0 7 0 64270 0 449890 406 333463 0 1182 0 13 0 70694 0 919022 407 344282 0 621 0 16 0 68005 0 1088080 408 319635 0 989 0 11 0 58930 0 648230 409 301186 0 438 0 11 0 58320 0 641520 410 300381 0 726 0 3 0 69980 0 209940 411 318765 0 1303 0 61 0 69863 0 4261643 412 286146 1 6341 6341 24 24 63255 63255 1518120 413 306844 1 1164 1164 17 17 57320 57320 974440 414 307705 1 3310 3310 33 33 75230 75230 2482590 415 312448 0 1366 0 7 0 79420 0 555940 416 299715 0 965 0 3 0 73490 0 220470 417 373399 0 3256 0 66 0 35250 0 2326500 418 299446 1 1135 1135 17 17 62285 62285 1058845 419 325586 0 1270 0 26 0 69206 0 1799356 420 291221 0 661 0 3 0 65920 0 197760 421 261173 0 1013 0 2 0 69770 0 139540 422 255027 0 2844 0 67 0 72683 0 4869761 423 -78375 1 11528 11528 70 70 -14545 -14545 -1018150 424 -58143 0 6526 0 26 0 55830 0 1451580 425 227033 0 2264 0 24 0 55174 0 1324176 426 235098 1 4461 4461 94 94 67038 67038 6301572 427 21267 0 3999 0 30 0 51252 0 1537560 428 238675 0 35624 0 223 0 157278 0 35072994 429 197687 0 9252 0 48 0 79510 0 3816480 430 418341 0 15236 0 90 0 77440 0 6969600 431 -297706 0 18073 0 180 0 27284 0 4911120 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs GrCosts Trades GrTrades 1.638e+05 7.066e+04 -5.434e+00 3.857e+01 -5.819e+02 -7.036e+01 Dividends GrDiv TrDiv 2.704e+00 -1.643e+00 6.407e-03 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -743164 -18133 -2401 21979 3055410 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.638e+05 3.431e+04 4.774 2.50e-06 *** Group 7.066e+04 5.587e+04 1.265 0.20671 Costs -5.434e+00 2.425e+00 -2.241 0.02553 * GrCosts 3.857e+01 3.744e+00 10.302 < 2e-16 *** Trades -5.819e+02 2.781e+02 -2.092 0.03703 * GrTrades -7.036e+01 3.405e+02 -0.207 0.83636 Dividends 2.704e+00 4.765e-01 5.675 2.58e-08 *** GrDiv -1.643e+00 7.088e-01 -2.318 0.02091 * TrDiv 6.407e-03 2.115e-03 3.029 0.00261 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 208400 on 422 degrees of freedom Multiple R-squared: 0.7871, Adjusted R-squared: 0.7831 F-statistic: 195 on 8 and 422 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 7.851251e-100 3.925626e-100 [2,] 1.0000000 2.191165e-113 1.095583e-113 [3,] 1.0000000 6.072675e-121 3.036338e-121 [4,] 1.0000000 2.671786e-125 1.335893e-125 [5,] 1.0000000 2.731368e-135 1.365684e-135 [6,] 1.0000000 1.148300e-137 5.741501e-138 [7,] 1.0000000 7.993861e-139 3.996931e-139 [8,] 1.0000000 7.659853e-146 3.829927e-146 [9,] 1.0000000 6.118071e-148 3.059036e-148 [10,] 1.0000000 8.791256e-154 4.395628e-154 [11,] 1.0000000 3.525189e-158 1.762594e-158 [12,] 1.0000000 4.274488e-159 2.137244e-159 [13,] 1.0000000 1.199647e-159 5.998233e-160 [14,] 1.0000000 6.364360e-159 3.182180e-159 [15,] 1.0000000 9.999843e-161 4.999921e-161 [16,] 1.0000000 4.209709e-164 2.104855e-164 [17,] 1.0000000 2.856631e-166 1.428315e-166 [18,] 1.0000000 5.112184e-179 2.556092e-179 [19,] 1.0000000 1.071187e-180 5.355933e-181 [20,] 1.0000000 2.749321e-184 1.374661e-184 [21,] 1.0000000 4.370172e-186 2.185086e-186 [22,] 1.0000000 2.256602e-192 1.128301e-192 [23,] 1.0000000 7.310422e-194 3.655211e-194 [24,] 1.0000000 7.134161e-193 3.567081e-193 [25,] 1.0000000 1.271116e-192 6.355578e-193 [26,] 1.0000000 7.610557e-200 3.805278e-200 [27,] 1.0000000 8.528374e-205 4.264187e-205 [28,] 1.0000000 4.029553e-204 2.014777e-204 [29,] 1.0000000 2.044370e-204 1.022185e-204 [30,] 1.0000000 2.069842e-206 1.034921e-206 [31,] 1.0000000 1.001833e-206 5.009163e-207 [32,] 1.0000000 7.168054e-207 3.584027e-207 [33,] 1.0000000 4.672260e-206 2.336130e-206 [34,] 1.0000000 1.134090e-205 5.670451e-206 [35,] 1.0000000 1.986613e-205 9.933063e-206 [36,] 1.0000000 2.011350e-204 1.005675e-204 [37,] 1.0000000 5.774158e-204 2.887079e-204 [38,] 1.0000000 2.333566e-203 1.166783e-203 [39,] 1.0000000 8.272180e-204 4.136090e-204 [40,] 1.0000000 7.636110e-205 3.818055e-205 [41,] 1.0000000 5.339899e-204 2.669949e-204 [42,] 1.0000000 2.070578e-203 1.035289e-203 [43,] 1.0000000 2.474863e-202 1.237432e-202 [44,] 1.0000000 2.854662e-207 1.427331e-207 [45,] 1.0000000 7.652691e-208 3.826346e-208 [46,] 1.0000000 2.359341e-212 1.179670e-212 [47,] 1.0000000 2.011848e-213 1.005924e-213 [48,] 1.0000000 1.527366e-212 7.636832e-213 [49,] 1.0000000 3.523830e-213 1.761915e-213 [50,] 1.0000000 1.206829e-214 6.034144e-215 [51,] 1.0000000 2.840292e-214 1.420146e-214 [52,] 1.0000000 1.118817e-213 5.594083e-214 [53,] 1.0000000 1.511618e-213 7.558088e-214 [54,] 1.0000000 3.172995e-213 1.586497e-213 [55,] 1.0000000 3.768960e-213 1.884480e-213 [56,] 1.0000000 3.646202e-214 1.823101e-214 [57,] 1.0000000 3.452764e-213 1.726382e-213 [58,] 1.0000000 1.532771e-212 7.663855e-213 [59,] 1.0000000 6.261774e-212 3.130887e-212 [60,] 1.0000000 1.366611e-211 6.833056e-212 [61,] 1.0000000 6.446728e-214 3.223364e-214 [62,] 1.0000000 1.523090e-213 7.615452e-214 [63,] 1.0000000 1.363732e-212 6.818661e-213 [64,] 1.0000000 1.055586e-211 5.277929e-212 [65,] 1.0000000 3.963795e-211 1.981897e-211 [66,] 1.0000000 1.178119e-210 5.890596e-211 [67,] 1.0000000 2.056105e-210 1.028052e-210 [68,] 1.0000000 2.275088e-209 1.137544e-209 [69,] 1.0000000 1.168665e-210 5.843323e-211 [70,] 1.0000000 8.805335e-210 4.402668e-210 [71,] 1.0000000 7.342751e-210 3.671376e-210 [72,] 1.0000000 3.768693e-210 1.884347e-210 [73,] 1.0000000 3.739645e-209 1.869822e-209 [74,] 1.0000000 3.589547e-208 1.794774e-208 [75,] 1.0000000 4.122509e-209 2.061254e-209 [76,] 1.0000000 3.914147e-209 1.957074e-209 [77,] 1.0000000 4.008049e-208 2.004024e-208 [78,] 1.0000000 1.199480e-207 5.997398e-208 [79,] 1.0000000 1.205720e-206 6.028602e-207 [80,] 1.0000000 1.288189e-205 6.440946e-206 [81,] 1.0000000 5.512733e-205 2.756367e-205 [82,] 1.0000000 7.561750e-206 3.780875e-206 [83,] 1.0000000 4.113166e-205 2.056583e-205 [84,] 1.0000000 5.673904e-206 2.836952e-206 [85,] 1.0000000 1.464725e-206 7.323624e-207 [86,] 1.0000000 6.771670e-206 3.385835e-206 [87,] 1.0000000 6.783791e-205 3.391896e-205 [88,] 1.0000000 5.970192e-204 2.985096e-204 [89,] 1.0000000 5.683514e-203 2.841757e-203 [90,] 1.0000000 5.157830e-202 2.578915e-202 [91,] 1.0000000 5.349010e-202 2.674505e-202 [92,] 1.0000000 2.878033e-201 1.439016e-201 [93,] 1.0000000 2.505846e-200 1.252923e-200 [94,] 1.0000000 1.052635e-199 5.263176e-200 [95,] 1.0000000 2.389309e-199 1.194655e-199 [96,] 1.0000000 2.079708e-200 1.039854e-200 [97,] 1.0000000 2.845376e-203 1.422688e-203 [98,] 1.0000000 2.353012e-202 1.176506e-202 [99,] 1.0000000 9.211686e-202 4.605843e-202 [100,] 1.0000000 4.993942e-202 2.496971e-202 [101,] 1.0000000 2.918974e-203 1.459487e-203 [102,] 1.0000000 6.757189e-205 3.378595e-205 [103,] 1.0000000 4.220176e-204 2.110088e-204 [104,] 1.0000000 1.434609e-203 7.173043e-204 [105,] 1.0000000 3.413438e-203 1.706719e-203 [106,] 1.0000000 3.329291e-202 1.664645e-202 [107,] 1.0000000 4.451413e-203 2.225707e-203 [108,] 1.0000000 8.239628e-204 4.119814e-204 [109,] 1.0000000 7.876436e-203 3.938218e-203 [110,] 1.0000000 4.385686e-204 2.192843e-204 [111,] 1.0000000 3.354392e-203 1.677196e-203 [112,] 1.0000000 1.351314e-202 6.756568e-203 [113,] 1.0000000 3.493354e-204 1.746677e-204 [114,] 1.0000000 3.420620e-203 1.710310e-203 [115,] 1.0000000 4.386036e-203 2.193018e-203 [116,] 1.0000000 4.941792e-202 2.470896e-202 [117,] 1.0000000 5.546052e-201 2.773026e-201 [118,] 1.0000000 6.158798e-200 3.079399e-200 [119,] 1.0000000 6.827413e-199 3.413707e-199 [120,] 1.0000000 7.536911e-198 3.768456e-198 [121,] 1.0000000 3.194940e-198 1.597470e-198 [122,] 1.0000000 3.521270e-197 1.760635e-197 [123,] 1.0000000 3.851471e-196 1.925736e-196 [124,] 1.0000000 4.200029e-195 2.100014e-195 [125,] 1.0000000 4.357628e-194 2.178814e-194 [126,] 1.0000000 4.711477e-193 2.355739e-193 [127,] 1.0000000 5.056614e-192 2.528307e-192 [128,] 1.0000000 5.398251e-191 2.699126e-191 [129,] 1.0000000 5.732416e-190 2.866208e-190 [130,] 1.0000000 6.013174e-189 3.006587e-189 [131,] 1.0000000 6.317827e-188 3.158913e-188 [132,] 1.0000000 5.250876e-187 2.625438e-187 [133,] 1.0000000 5.466402e-186 2.733201e-186 [134,] 1.0000000 4.527111e-185 2.263556e-185 [135,] 1.0000000 4.460852e-184 2.230426e-184 [136,] 1.0000000 4.571938e-183 2.285969e-183 [137,] 1.0000000 4.661299e-182 2.330649e-182 [138,] 1.0000000 4.727525e-181 2.363763e-181 [139,] 1.0000000 4.770023e-180 2.385012e-180 [140,] 1.0000000 4.760345e-179 2.380172e-179 [141,] 1.0000000 3.289288e-178 1.644644e-178 [142,] 1.0000000 3.271132e-177 1.635566e-177 [143,] 1.0000000 3.039135e-176 1.519567e-176 [144,] 1.0000000 2.992293e-175 1.496147e-175 [145,] 1.0000000 2.930747e-174 1.465373e-174 [146,] 1.0000000 2.855705e-173 1.427853e-173 [147,] 1.0000000 2.767403e-172 1.383701e-172 [148,] 1.0000000 2.232514e-171 1.116257e-171 [149,] 1.0000000 2.143563e-170 1.071782e-170 [150,] 1.0000000 2.047659e-169 1.023829e-169 [151,] 1.0000000 1.946088e-168 9.730439e-169 [152,] 1.0000000 1.840156e-167 9.200782e-168 [153,] 1.0000000 1.619920e-166 8.099599e-167 [154,] 1.0000000 1.515315e-165 7.576574e-166 [155,] 1.0000000 1.411462e-164 7.057310e-165 [156,] 1.0000000 1.308099e-163 6.540497e-164 [157,] 1.0000000 1.999802e-164 9.999011e-165 [158,] 1.0000000 1.864968e-163 9.324838e-164 [159,] 1.0000000 1.730530e-162 8.652649e-163 [160,] 1.0000000 1.597768e-161 7.988842e-162 [161,] 1.0000000 1.465500e-160 7.327501e-161 [162,] 1.0000000 1.339647e-159 6.698234e-160 [163,] 1.0000000 1.208928e-158 6.044638e-159 [164,] 1.0000000 1.094183e-157 5.470916e-158 [165,] 1.0000000 9.856067e-157 4.928033e-157 [166,] 1.0000000 8.832444e-156 4.416222e-156 [167,] 1.0000000 7.892415e-155 3.946208e-155 [168,] 1.0000000 7.002654e-154 3.501327e-154 [169,] 1.0000000 6.182615e-153 3.091308e-153 [170,] 1.0000000 1.699175e-152 8.495876e-153 [171,] 1.0000000 1.273459e-151 6.367296e-152 [172,] 1.0000000 8.111266e-151 4.055633e-151 [173,] 1.0000000 5.208978e-150 2.604489e-150 [174,] 1.0000000 4.530619e-149 2.265309e-149 [175,] 1.0000000 3.778162e-148 1.889081e-148 [176,] 1.0000000 2.040727e-147 1.020363e-147 [177,] 1.0000000 1.741004e-146 8.705019e-147 [178,] 1.0000000 1.228971e-145 6.144856e-146 [179,] 1.0000000 1.044676e-144 5.223381e-145 [180,] 1.0000000 8.320041e-144 4.160021e-144 [181,] 1.0000000 7.005587e-143 3.502793e-143 [182,] 1.0000000 5.752095e-142 2.876047e-142 [183,] 1.0000000 4.796909e-141 2.398455e-141 [184,] 1.0000000 3.980739e-140 1.990370e-140 [185,] 1.0000000 3.290844e-139 1.645422e-139 [186,] 1.0000000 2.390340e-138 1.195170e-138 [187,] 1.0000000 1.958266e-137 9.791331e-138 [188,] 1.0000000 1.596482e-136 7.982410e-137 [189,] 1.0000000 9.867813e-136 4.933906e-136 [190,] 1.0000000 7.973688e-135 3.986844e-135 [191,] 1.0000000 6.411728e-134 3.205864e-134 [192,] 1.0000000 5.130600e-133 2.565300e-133 [193,] 1.0000000 4.089697e-132 2.044848e-132 [194,] 1.0000000 2.919224e-131 1.459612e-131 [195,] 1.0000000 2.301929e-130 1.150965e-130 [196,] 1.0000000 1.796011e-129 8.980055e-130 [197,] 1.0000000 1.402682e-128 7.013411e-129 [198,] 1.0000000 2.338485e-128 1.169242e-128 [199,] 1.0000000 1.824807e-127 9.124033e-128 [200,] 1.0000000 1.415115e-126 7.075575e-127 [201,] 1.0000000 1.092154e-125 5.460771e-126 [202,] 1.0000000 8.122477e-125 4.061239e-125 [203,] 1.0000000 6.072738e-124 3.036369e-124 [204,] 1.0000000 4.567646e-123 2.283823e-123 [205,] 1.0000000 3.459292e-122 1.729646e-122 [206,] 1.0000000 2.607234e-121 1.303617e-121 [207,] 1.0000000 1.897200e-120 9.486001e-121 [208,] 1.0000000 1.416463e-119 7.082315e-120 [209,] 1.0000000 6.630758e-119 3.315379e-119 [210,] 1.0000000 4.210428e-118 2.105214e-118 [211,] 1.0000000 3.108270e-117 1.554135e-117 [212,] 1.0000000 2.283718e-116 1.141859e-116 [213,] 1.0000000 1.669771e-115 8.348853e-116 [214,] 1.0000000 1.214805e-114 6.074023e-115 [215,] 1.0000000 8.659823e-114 4.329911e-114 [216,] 1.0000000 6.241285e-113 3.120642e-113 [217,] 1.0000000 4.476565e-112 2.238283e-112 [218,] 1.0000000 3.070093e-111 1.535047e-111 [219,] 1.0000000 2.181458e-110 1.090729e-110 [220,] 1.0000000 1.527339e-109 7.636694e-110 [221,] 1.0000000 1.074895e-108 5.374477e-109 [222,] 1.0000000 7.528258e-108 3.764129e-108 [223,] 1.0000000 5.010585e-107 2.505293e-107 [224,] 1.0000000 3.476160e-106 1.738080e-106 [225,] 1.0000000 2.400047e-105 1.200023e-105 [226,] 1.0000000 1.648967e-104 8.244835e-105 [227,] 1.0000000 1.125397e-103 5.626986e-104 [228,] 1.0000000 7.657198e-103 3.828599e-103 [229,] 1.0000000 5.184612e-102 2.592306e-102 [230,] 1.0000000 3.493353e-101 1.746676e-101 [231,] 1.0000000 2.304431e-100 1.152215e-100 [232,] 1.0000000 1.537686e-99 7.688429e-100 [233,] 1.0000000 1.020333e-98 5.101666e-99 [234,] 1.0000000 4.464562e-98 2.232281e-98 [235,] 1.0000000 2.944635e-97 1.472318e-97 [236,] 1.0000000 1.932582e-96 9.662911e-97 [237,] 1.0000000 1.253793e-95 6.268963e-96 [238,] 1.0000000 8.147958e-95 4.073979e-95 [239,] 1.0000000 5.268838e-94 2.634419e-94 [240,] 1.0000000 3.390350e-93 1.695175e-93 [241,] 1.0000000 2.170645e-92 1.085322e-92 [242,] 1.0000000 1.346710e-91 6.733550e-92 [243,] 1.0000000 7.624759e-91 3.812380e-91 [244,] 1.0000000 2.223465e-90 1.111733e-90 [245,] 1.0000000 1.401791e-89 7.008956e-90 [246,] 1.0000000 7.100697e-89 3.550348e-89 [247,] 1.0000000 3.679878e-88 1.839939e-88 [248,] 1.0000000 2.293719e-87 1.146859e-87 [249,] 1.0000000 1.420145e-86 7.100723e-87 [250,] 1.0000000 8.762918e-86 4.381459e-86 [251,] 1.0000000 5.379695e-85 2.689847e-85 [252,] 1.0000000 3.275713e-84 1.637856e-84 [253,] 1.0000000 1.892728e-83 9.463641e-84 [254,] 1.0000000 1.144745e-82 5.723726e-83 [255,] 1.0000000 6.876357e-82 3.438178e-82 [256,] 1.0000000 4.116488e-81 2.058244e-81 [257,] 1.0000000 1.874660e-80 9.373301e-81 [258,] 1.0000000 1.112719e-79 5.563593e-80 [259,] 1.0000000 6.554998e-79 3.277499e-79 [260,] 1.0000000 3.850680e-78 1.925340e-78 [261,] 1.0000000 2.250268e-77 1.125134e-77 [262,] 1.0000000 1.308141e-76 6.540707e-77 [263,] 1.0000000 7.563932e-76 3.781966e-76 [264,] 1.0000000 4.344533e-75 2.172267e-75 [265,] 1.0000000 2.432268e-74 1.216134e-74 [266,] 1.0000000 1.351404e-73 6.757022e-74 [267,] 1.0000000 7.642411e-73 3.821205e-73 [268,] 1.0000000 3.979334e-72 1.989667e-72 [269,] 1.0000000 2.073574e-71 1.036787e-71 [270,] 1.0000000 1.157067e-70 5.785337e-71 [271,] 1.0000000 6.413274e-70 3.206637e-70 [272,] 1.0000000 3.424981e-69 1.712490e-69 [273,] 1.0000000 1.878607e-68 9.393035e-69 [274,] 1.0000000 9.962060e-68 4.981030e-68 [275,] 1.0000000 5.406997e-67 2.703499e-67 [276,] 1.0000000 2.922362e-66 1.461181e-66 [277,] 1.0000000 1.568347e-65 7.841735e-66 [278,] 1.0000000 8.383805e-65 4.191902e-65 [279,] 1.0000000 4.385742e-64 2.192871e-64 [280,] 1.0000000 2.318883e-63 1.159441e-63 [281,] 1.0000000 1.196050e-62 5.980249e-63 [282,] 1.0000000 6.092224e-62 3.046112e-62 [283,] 1.0000000 3.165327e-61 1.582663e-61 [284,] 1.0000000 1.637804e-60 8.189018e-61 [285,] 1.0000000 6.441852e-60 3.220926e-60 [286,] 1.0000000 3.298147e-59 1.649074e-59 [287,] 1.0000000 1.681981e-58 8.409904e-59 [288,] 1.0000000 8.513785e-58 4.256892e-58 [289,] 1.0000000 4.284748e-57 2.142374e-57 [290,] 1.0000000 2.147118e-56 1.073559e-56 [291,] 1.0000000 1.067749e-55 5.338743e-56 [292,] 1.0000000 5.281349e-55 2.640674e-55 [293,] 1.0000000 2.601522e-54 1.300761e-54 [294,] 1.0000000 1.036493e-53 5.182465e-54 [295,] 1.0000000 5.047364e-53 2.523682e-53 [296,] 1.0000000 2.447067e-52 1.223534e-52 [297,] 1.0000000 1.179226e-51 5.896132e-52 [298,] 1.0000000 5.062402e-51 2.531201e-51 [299,] 1.0000000 2.412903e-50 1.206452e-50 [300,] 1.0000000 1.143088e-49 5.715441e-50 [301,] 1.0000000 3.602734e-49 1.801367e-49 [302,] 1.0000000 1.692441e-48 8.462204e-49 [303,] 1.0000000 7.891565e-48 3.945783e-48 [304,] 1.0000000 3.653582e-47 1.826791e-47 [305,] 1.0000000 1.664452e-46 8.322259e-47 [306,] 1.0000000 7.589775e-46 3.794887e-46 [307,] 1.0000000 3.447869e-45 1.723934e-45 [308,] 1.0000000 1.557706e-44 7.788532e-45 [309,] 1.0000000 6.984510e-44 3.492255e-44 [310,] 1.0000000 3.111567e-43 1.555783e-43 [311,] 1.0000000 1.377199e-42 6.885994e-43 [312,] 1.0000000 6.055769e-42 3.027884e-42 [313,] 1.0000000 2.652684e-41 1.326342e-41 [314,] 1.0000000 1.151144e-40 5.755722e-41 [315,] 1.0000000 4.973864e-40 2.486932e-40 [316,] 1.0000000 2.135430e-39 1.067715e-39 [317,] 1.0000000 9.104121e-39 4.552060e-39 [318,] 1.0000000 3.845887e-38 1.922944e-38 [319,] 1.0000000 1.613490e-37 8.067452e-38 [320,] 1.0000000 6.722403e-37 3.361201e-37 [321,] 1.0000000 2.781298e-36 1.390649e-36 [322,] 1.0000000 1.142645e-35 5.713223e-36 [323,] 1.0000000 4.655590e-35 2.327795e-35 [324,] 1.0000000 1.886435e-34 9.432175e-35 [325,] 1.0000000 7.585461e-34 3.792731e-34 [326,] 1.0000000 3.033512e-33 1.516756e-33 [327,] 1.0000000 1.202116e-32 6.010580e-33 [328,] 1.0000000 4.727634e-32 2.363817e-32 [329,] 1.0000000 1.845013e-31 9.225067e-32 [330,] 1.0000000 7.120755e-31 3.560377e-31 [331,] 1.0000000 2.739196e-30 1.369598e-30 [332,] 1.0000000 1.045648e-29 5.228241e-30 [333,] 1.0000000 3.960789e-29 1.980395e-29 [334,] 1.0000000 1.488602e-28 7.443009e-29 [335,] 1.0000000 5.550629e-28 2.775315e-28 [336,] 1.0000000 2.058112e-27 1.029056e-27 [337,] 1.0000000 7.553423e-27 3.776711e-27 [338,] 1.0000000 2.749663e-26 1.374831e-26 [339,] 1.0000000 9.927456e-26 4.963728e-26 [340,] 1.0000000 3.531216e-25 1.765608e-25 [341,] 1.0000000 1.256723e-24 6.283614e-25 [342,] 1.0000000 4.127511e-24 2.063755e-24 [343,] 1.0000000 1.446441e-23 7.232207e-24 [344,] 1.0000000 5.008317e-23 2.504158e-23 [345,] 1.0000000 1.718906e-22 8.594531e-23 [346,] 1.0000000 5.164669e-22 2.582335e-22 [347,] 1.0000000 1.743900e-21 8.719502e-22 [348,] 1.0000000 5.834753e-21 2.917376e-21 [349,] 1.0000000 1.934160e-20 9.670800e-21 [350,] 1.0000000 6.363140e-20 3.181570e-20 [351,] 1.0000000 2.069828e-19 1.034914e-19 [352,] 1.0000000 6.668038e-19 3.334019e-19 [353,] 1.0000000 2.139356e-18 1.069678e-18 [354,] 1.0000000 6.782242e-18 3.391121e-18 [355,] 1.0000000 2.121598e-17 1.060799e-17 [356,] 1.0000000 6.569150e-17 3.284575e-17 [357,] 1.0000000 2.012988e-16 1.006494e-16 [358,] 1.0000000 6.103593e-16 3.051796e-16 [359,] 1.0000000 1.830904e-15 9.154522e-16 [360,] 1.0000000 5.432516e-15 2.716258e-15 [361,] 1.0000000 1.608125e-14 8.040626e-15 [362,] 1.0000000 4.629544e-14 2.314772e-14 [363,] 1.0000000 1.336757e-13 6.683787e-14 [364,] 1.0000000 3.797814e-13 1.898907e-13 [365,] 1.0000000 1.066256e-12 5.331282e-13 [366,] 1.0000000 2.972788e-12 1.486394e-12 [367,] 1.0000000 8.084023e-12 4.042012e-12 [368,] 1.0000000 2.187892e-11 1.093946e-11 [369,] 1.0000000 5.845773e-11 2.922886e-11 [370,] 1.0000000 1.541517e-10 7.707586e-11 [371,] 1.0000000 3.873563e-10 1.936781e-10 [372,] 1.0000000 9.260512e-10 4.630256e-10 [373,] 1.0000000 2.354000e-09 1.177000e-09 [374,] 1.0000000 5.951081e-09 2.975541e-09 [375,] 1.0000000 1.213660e-08 6.068299e-09 [376,] 1.0000000 2.948997e-08 1.474499e-08 [377,] 1.0000000 6.727029e-08 3.363515e-08 [378,] 0.9999999 1.603641e-07 8.018205e-08 [379,] 0.9999998 3.792423e-07 1.896212e-07 [380,] 0.9999996 8.630042e-07 4.315021e-07 [381,] 0.9999990 1.973840e-06 9.869201e-07 [382,] 0.9999978 4.469783e-06 2.234892e-06 [383,] 0.9999958 8.448100e-06 4.224050e-06 [384,] 0.9999914 1.722933e-05 8.614666e-06 [385,] 0.9999814 3.722305e-05 1.861153e-05 [386,] 0.9999656 6.872351e-05 3.436175e-05 [387,] 0.9999299 1.402851e-04 7.014253e-05 [388,] 0.9998565 2.869441e-04 1.434721e-04 [389,] 0.9997498 5.003118e-04 2.501559e-04 [390,] 0.9995169 9.661836e-04 4.830918e-04 [391,] 0.9991177 1.764696e-03 8.823479e-04 [392,] 0.9983704 3.259159e-03 1.629580e-03 [393,] 0.9969592 6.081585e-03 3.040793e-03 [394,] 0.9945000 1.099991e-02 5.499953e-03 [395,] 0.9902795 1.944102e-02 9.720508e-03 [396,] 0.9834526 3.309474e-02 1.654737e-02 [397,] 0.9740392 5.192155e-02 2.596078e-02 [398,] 0.9583871 8.322573e-02 4.161287e-02 [399,] 0.9330291 1.339418e-01 6.697091e-02 [400,] 0.8963089 2.073822e-01 1.036911e-01 [401,] 0.8463392 3.073216e-01 1.536608e-01 [402,] 0.7773996 4.452009e-01 2.226004e-01 [403,] 0.6896033 6.207933e-01 3.103967e-01 [404,] 0.5858268 8.283464e-01 4.141732e-01 [405,] 0.4683597 9.367194e-01 5.316403e-01 [406,] 0.7841718 4.316563e-01 2.158282e-01 [407,] 0.6533324 6.933352e-01 3.466676e-01 [408,] 0.4907134 9.814268e-01 5.092866e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1frax1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/2frax1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/3frax1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/4q0si1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5q0si1291330094.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 -335057.38899 3055410.06433 659410.54377 -405429.22905 -743164.38986 6 7 8 9 10 -89456.37662 957456.07911 358218.87563 -515485.62638 656837.66595 11 12 13 14 15 418629.65698 790336.11651 786396.36249 -496204.59220 -198517.34823 16 17 18 19 20 573684.82775 -252041.74304 -311194.45538 -428313.88915 -21431.49731 21 22 23 24 25 139295.64762 -335211.15220 114533.03501 -174367.24805 -66395.70200 26 27 28 29 30 301730.34172 -543727.48346 -174098.17192 593025.71736 341865.12044 31 32 33 34 35 252881.74560 -132561.70348 163074.51332 105181.03944 -175459.84008 36 37 38 39 40 26963.16687 306078.95930 161359.42712 138591.64288 -121232.08493 41 42 43 44 45 -315987.14980 21458.74625 104529.24113 -48445.33886 62799.44298 46 47 48 49 50 -127194.57210 118925.53513 -219278.61704 -14451.06244 -142981.86240 51 52 53 54 55 105768.72434 -34627.40188 -219322.05982 70651.42675 -281080.27993 56 57 58 59 60 124875.96679 232448.43092 85407.63172 -63984.83679 -220578.16036 61 62 63 64 65 -23272.87970 -78225.09783 -20749.17848 -216463.75382 68796.51174 66 67 68 69 70 23572.03863 -137137.24569 -268.85031 -84555.82611 -55112.87409 71 72 73 74 75 -89236.14587 153263.11954 -121884.06595 -46107.67637 -94572.99222 76 77 78 79 80 -73086.17103 -105602.61775 -85644.93657 37088.85778 21253.31059 81 82 83 84 85 29611.70914 -45462.15309 -134129.65691 36887.62176 -23125.20612 86 87 88 89 90 -116597.88333 70596.47129 -54748.11133 -88191.48526 -33468.48369 91 92 93 94 95 -20293.38577 -41947.40635 122154.90850 -228180.60334 119037.29664 96 97 98 99 100 -11477.33974 -1554.05235 49267.09848 -71063.71119 -26793.91038 101 102 103 104 105 -34155.89201 -116393.68974 -71182.11744 -63711.14507 39309.00263 106 107 108 109 110 -78901.81531 -233869.45225 107762.43360 -87569.45515 -104668.69729 111 112 113 114 115 82970.34063 -253759.38582 -175074.28537 -121209.46808 14252.74836 116 117 118 119 120 55220.38255 -26722.55853 65798.28707 -176082.65391 26570.83341 121 122 123 124 125 2477.23170 13431.45499 -36512.52694 80946.69568 -100397.40476 126 127 128 129 130 -201133.25705 -2401.22685 -987.01213 9809.10476 -2401.22685 131 132 133 134 135 -386.54690 85481.18719 -2401.22685 -3753.75972 -2401.22685 136 137 138 139 140 7599.65197 175.79713 -2401.22685 -2401.22685 -2401.22685 141 142 143 144 145 3951.31681 -2401.22685 17400.26177 -1675.86504 16223.04061 146 147 148 149 150 13349.51735 -2401.22685 -1714.53127 -2401.22685 -2564.61669 151 152 153 154 155 543.65038 22998.67080 -2401.22685 -14612.79492 -2327.64702 156 157 158 159 160 -2401.22685 -2401.22685 -2399.01202 15826.16411 -2401.22685 161 162 163 164 165 -2401.22685 -2401.22685 -2401.22685 11382.48822 -60.46908 166 167 168 169 170 -2401.22685 -2401.22685 110015.49676 -2401.22685 -2401.22685 171 172 173 174 175 -2401.22685 -4700.09393 -2401.22685 239.82022 -2401.22685 176 177 178 179 180 -1522.71355 -2401.22685 24853.40506 -2401.22685 -2401.22685 181 182 183 184 185 39696.20723 23643.86380 29378.05747 -12686.20265 29249.09785 186 187 188 189 190 11165.20288 33626.14463 -1350.71885 -32456.24329 -2401.22685 191 192 193 194 195 8980.09568 -2401.22685 3290.42860 -2401.22685 -2200.22282 196 197 198 199 200 26723.56888 -41404.00191 26723.56888 26723.56888 -13875.78642 201 202 203 204 205 26723.56888 26723.56888 26723.56888 25585.75864 -4271.33958 206 207 208 209 210 26784.69571 22287.18057 26723.56888 44234.67155 26723.56888 211 212 213 214 215 -2401.22685 -2401.22685 11362.80316 -14186.21899 -4637.33104 216 217 218 219 220 -2401.22685 -486.64678 6769.86053 -2401.22685 27804.36286 221 222 223 224 225 -27290.37324 -2401.22685 -1963.41995 -2401.22685 -1935.75325 226 227 228 229 230 4557.43857 -2401.22685 -2401.22685 -20062.40201 -2401.22685 231 232 233 234 235 -13551.83006 -2401.22685 -2401.22685 6727.90438 -2401.22685 236 237 238 239 240 -2333.83700 -2401.22685 339.50561 -2401.22685 -2401.22685 241 242 243 244 245 -2401.22685 3152.22517 -2401.22685 720.53012 -7170.17477 246 247 248 249 250 26723.56888 26723.56888 -8564.63380 26723.56888 26723.56888 251 252 253 254 255 3156.75624 26723.56888 -18278.38317 -17986.72143 -42817.86848 256 257 258 259 260 -2401.22685 26006.40867 32093.64563 26723.56888 -2401.22685 261 262 263 264 265 26723.56888 26723.56888 -9433.34743 10546.87598 26723.56888 266 267 268 269 270 -2401.22685 26723.56888 16171.00545 26723.56888 -14299.80626 271 272 273 274 275 26723.56888 26723.56888 26723.56888 -3108.46350 -2401.22685 276 277 278 279 280 8612.68765 8347.42108 -2401.22685 6172.76214 -20844.16653 281 282 283 284 285 26723.56888 -2401.22685 23272.10884 -2401.22685 -24024.75157 286 287 288 289 290 -2401.22685 26723.56888 -2486.74385 26723.56888 16783.25672 291 292 293 294 295 26723.56888 7464.49902 -6947.57562 -2401.22685 -1737.31668 296 297 298 299 300 51474.88847 -2401.22685 26723.56888 -2401.22685 -2401.22685 301 302 303 304 305 25184.63237 -9291.67440 -2401.22685 26723.56888 34355.02811 306 307 308 309 310 -666.00381 26723.56888 22464.83173 21670.74259 26723.56888 311 312 313 314 315 26723.56888 -26255.79903 26723.56888 29891.95209 -1414.38736 316 317 318 319 320 2902.88817 -643.29997 -9108.09383 -1238.86183 -2401.22685 321 322 323 324 325 -2401.22685 -2401.22685 -2401.22685 26723.56888 -2401.22685 326 327 328 329 330 -4812.37279 26723.56888 26723.56888 -2401.22685 -2401.22685 331 332 333 334 335 -2401.22685 -2401.22685 -2401.22685 -944.13962 -2022.08210 336 337 338 339 340 -2401.22685 -2729.53365 -2401.22685 -2356.13091 1358.77905 341 342 343 344 345 -14305.02958 -2401.22685 -2401.22685 -2401.22685 -2401.22685 346 347 348 349 350 -2401.22685 -9488.63760 -2401.22685 -2401.22685 -2401.22685 351 352 353 354 355 2939.06427 -1578.88762 13632.90069 -7486.21423 -2401.22685 356 357 358 359 360 -2401.22685 -24740.26924 -2401.22685 -2401.22685 -2401.22685 361 362 363 364 365 -1630.79340 -2401.22685 -2401.22685 18714.84729 -1899.99534 366 367 368 369 370 -2401.22685 -2401.22685 -2401.22685 -2401.22685 -2401.22685 371 372 373 374 375 -2401.22685 3486.88107 -25749.55270 -20702.21744 -2401.22685 376 377 378 379 380 -2401.22685 4330.91660 -1503.97547 -2401.22685 -2401.22685 381 382 383 384 385 -2401.22685 -3409.33390 41667.47187 -2401.22685 -763.76613 386 387 388 389 390 -86823.08487 3994.32689 7260.16022 -2997.41084 -9409.52872 391 392 393 394 395 5753.07660 -11857.94258 -22132.20075 39537.25471 -155995.44039 396 397 398 399 400 -17357.95088 -25973.91882 107289.31148 -12279.43491 -84213.96533 401 402 403 404 405 -496.61868 12242.61718 -7008.00459 -51150.14027 -16445.04507 406 407 408 409 410 -13369.73486 2334.62395 4136.64604 -15614.21169 -48274.99259 411 412 413 414 415 -18647.91267 -219572.24477 -22111.50353 -110588.42235 -58146.69165 416 417 418 419 420 -57201.32722 155507.03181 -34356.48016 -14820.25975 -46731.17993 421 422 423 424 425 -85486.26190 -82041.92819 -627186.38066 -331586.22377 -68145.58823 426 427 428 429 430 -197325.48460 -251752.72344 -251741.39555 -127326.89068 135682.64604 431 -363758.36870 > postscript(file="/var/www/html/freestat/rcomp/tmp/6q0si1291330094.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 -335057.38899 NA 1 3055410.06433 -335057.38899 2 659410.54377 3055410.06433 3 -405429.22905 659410.54377 4 -743164.38986 -405429.22905 5 -89456.37662 -743164.38986 6 957456.07911 -89456.37662 7 358218.87563 957456.07911 8 -515485.62638 358218.87563 9 656837.66595 -515485.62638 10 418629.65698 656837.66595 11 790336.11651 418629.65698 12 786396.36249 790336.11651 13 -496204.59220 786396.36249 14 -198517.34823 -496204.59220 15 573684.82775 -198517.34823 16 -252041.74304 573684.82775 17 -311194.45538 -252041.74304 18 -428313.88915 -311194.45538 19 -21431.49731 -428313.88915 20 139295.64762 -21431.49731 21 -335211.15220 139295.64762 22 114533.03501 -335211.15220 23 -174367.24805 114533.03501 24 -66395.70200 -174367.24805 25 301730.34172 -66395.70200 26 -543727.48346 301730.34172 27 -174098.17192 -543727.48346 28 593025.71736 -174098.17192 29 341865.12044 593025.71736 30 252881.74560 341865.12044 31 -132561.70348 252881.74560 32 163074.51332 -132561.70348 33 105181.03944 163074.51332 34 -175459.84008 105181.03944 35 26963.16687 -175459.84008 36 306078.95930 26963.16687 37 161359.42712 306078.95930 38 138591.64288 161359.42712 39 -121232.08493 138591.64288 40 -315987.14980 -121232.08493 41 21458.74625 -315987.14980 42 104529.24113 21458.74625 43 -48445.33886 104529.24113 44 62799.44298 -48445.33886 45 -127194.57210 62799.44298 46 118925.53513 -127194.57210 47 -219278.61704 118925.53513 48 -14451.06244 -219278.61704 49 -142981.86240 -14451.06244 50 105768.72434 -142981.86240 51 -34627.40188 105768.72434 52 -219322.05982 -34627.40188 53 70651.42675 -219322.05982 54 -281080.27993 70651.42675 55 124875.96679 -281080.27993 56 232448.43092 124875.96679 57 85407.63172 232448.43092 58 -63984.83679 85407.63172 59 -220578.16036 -63984.83679 60 -23272.87970 -220578.16036 61 -78225.09783 -23272.87970 62 -20749.17848 -78225.09783 63 -216463.75382 -20749.17848 64 68796.51174 -216463.75382 65 23572.03863 68796.51174 66 -137137.24569 23572.03863 67 -268.85031 -137137.24569 68 -84555.82611 -268.85031 69 -55112.87409 -84555.82611 70 -89236.14587 -55112.87409 71 153263.11954 -89236.14587 72 -121884.06595 153263.11954 73 -46107.67637 -121884.06595 74 -94572.99222 -46107.67637 75 -73086.17103 -94572.99222 76 -105602.61775 -73086.17103 77 -85644.93657 -105602.61775 78 37088.85778 -85644.93657 79 21253.31059 37088.85778 80 29611.70914 21253.31059 81 -45462.15309 29611.70914 82 -134129.65691 -45462.15309 83 36887.62176 -134129.65691 84 -23125.20612 36887.62176 85 -116597.88333 -23125.20612 86 70596.47129 -116597.88333 87 -54748.11133 70596.47129 88 -88191.48526 -54748.11133 89 -33468.48369 -88191.48526 90 -20293.38577 -33468.48369 91 -41947.40635 -20293.38577 92 122154.90850 -41947.40635 93 -228180.60334 122154.90850 94 119037.29664 -228180.60334 95 -11477.33974 119037.29664 96 -1554.05235 -11477.33974 97 49267.09848 -1554.05235 98 -71063.71119 49267.09848 99 -26793.91038 -71063.71119 100 -34155.89201 -26793.91038 101 -116393.68974 -34155.89201 102 -71182.11744 -116393.68974 103 -63711.14507 -71182.11744 104 39309.00263 -63711.14507 105 -78901.81531 39309.00263 106 -233869.45225 -78901.81531 107 107762.43360 -233869.45225 108 -87569.45515 107762.43360 109 -104668.69729 -87569.45515 110 82970.34063 -104668.69729 111 -253759.38582 82970.34063 112 -175074.28537 -253759.38582 113 -121209.46808 -175074.28537 114 14252.74836 -121209.46808 115 55220.38255 14252.74836 116 -26722.55853 55220.38255 117 65798.28707 -26722.55853 118 -176082.65391 65798.28707 119 26570.83341 -176082.65391 120 2477.23170 26570.83341 121 13431.45499 2477.23170 122 -36512.52694 13431.45499 123 80946.69568 -36512.52694 124 -100397.40476 80946.69568 125 -201133.25705 -100397.40476 126 -2401.22685 -201133.25705 127 -987.01213 -2401.22685 128 9809.10476 -987.01213 129 -2401.22685 9809.10476 130 -386.54690 -2401.22685 131 85481.18719 -386.54690 132 -2401.22685 85481.18719 133 -3753.75972 -2401.22685 134 -2401.22685 -3753.75972 135 7599.65197 -2401.22685 136 175.79713 7599.65197 137 -2401.22685 175.79713 138 -2401.22685 -2401.22685 139 -2401.22685 -2401.22685 140 3951.31681 -2401.22685 141 -2401.22685 3951.31681 142 17400.26177 -2401.22685 143 -1675.86504 17400.26177 144 16223.04061 -1675.86504 145 13349.51735 16223.04061 146 -2401.22685 13349.51735 147 -1714.53127 -2401.22685 148 -2401.22685 -1714.53127 149 -2564.61669 -2401.22685 150 543.65038 -2564.61669 151 22998.67080 543.65038 152 -2401.22685 22998.67080 153 -14612.79492 -2401.22685 154 -2327.64702 -14612.79492 155 -2401.22685 -2327.64702 156 -2401.22685 -2401.22685 157 -2399.01202 -2401.22685 158 15826.16411 -2399.01202 159 -2401.22685 15826.16411 160 -2401.22685 -2401.22685 161 -2401.22685 -2401.22685 162 -2401.22685 -2401.22685 163 11382.48822 -2401.22685 164 -60.46908 11382.48822 165 -2401.22685 -60.46908 166 -2401.22685 -2401.22685 167 110015.49676 -2401.22685 168 -2401.22685 110015.49676 169 -2401.22685 -2401.22685 170 -2401.22685 -2401.22685 171 -4700.09393 -2401.22685 172 -2401.22685 -4700.09393 173 239.82022 -2401.22685 174 -2401.22685 239.82022 175 -1522.71355 -2401.22685 176 -2401.22685 -1522.71355 177 24853.40506 -2401.22685 178 -2401.22685 24853.40506 179 -2401.22685 -2401.22685 180 39696.20723 -2401.22685 181 23643.86380 39696.20723 182 29378.05747 23643.86380 183 -12686.20265 29378.05747 184 29249.09785 -12686.20265 185 11165.20288 29249.09785 186 33626.14463 11165.20288 187 -1350.71885 33626.14463 188 -32456.24329 -1350.71885 189 -2401.22685 -32456.24329 190 8980.09568 -2401.22685 191 -2401.22685 8980.09568 192 3290.42860 -2401.22685 193 -2401.22685 3290.42860 194 -2200.22282 -2401.22685 195 26723.56888 -2200.22282 196 -41404.00191 26723.56888 197 26723.56888 -41404.00191 198 26723.56888 26723.56888 199 -13875.78642 26723.56888 200 26723.56888 -13875.78642 201 26723.56888 26723.56888 202 26723.56888 26723.56888 203 25585.75864 26723.56888 204 -4271.33958 25585.75864 205 26784.69571 -4271.33958 206 22287.18057 26784.69571 207 26723.56888 22287.18057 208 44234.67155 26723.56888 209 26723.56888 44234.67155 210 -2401.22685 26723.56888 211 -2401.22685 -2401.22685 212 11362.80316 -2401.22685 213 -14186.21899 11362.80316 214 -4637.33104 -14186.21899 215 -2401.22685 -4637.33104 216 -486.64678 -2401.22685 217 6769.86053 -486.64678 218 -2401.22685 6769.86053 219 27804.36286 -2401.22685 220 -27290.37324 27804.36286 221 -2401.22685 -27290.37324 222 -1963.41995 -2401.22685 223 -2401.22685 -1963.41995 224 -1935.75325 -2401.22685 225 4557.43857 -1935.75325 226 -2401.22685 4557.43857 227 -2401.22685 -2401.22685 228 -20062.40201 -2401.22685 229 -2401.22685 -20062.40201 230 -13551.83006 -2401.22685 231 -2401.22685 -13551.83006 232 -2401.22685 -2401.22685 233 6727.90438 -2401.22685 234 -2401.22685 6727.90438 235 -2333.83700 -2401.22685 236 -2401.22685 -2333.83700 237 339.50561 -2401.22685 238 -2401.22685 339.50561 239 -2401.22685 -2401.22685 240 -2401.22685 -2401.22685 241 3152.22517 -2401.22685 242 -2401.22685 3152.22517 243 720.53012 -2401.22685 244 -7170.17477 720.53012 245 26723.56888 -7170.17477 246 26723.56888 26723.56888 247 -8564.63380 26723.56888 248 26723.56888 -8564.63380 249 26723.56888 26723.56888 250 3156.75624 26723.56888 251 26723.56888 3156.75624 252 -18278.38317 26723.56888 253 -17986.72143 -18278.38317 254 -42817.86848 -17986.72143 255 -2401.22685 -42817.86848 256 26006.40867 -2401.22685 257 32093.64563 26006.40867 258 26723.56888 32093.64563 259 -2401.22685 26723.56888 260 26723.56888 -2401.22685 261 26723.56888 26723.56888 262 -9433.34743 26723.56888 263 10546.87598 -9433.34743 264 26723.56888 10546.87598 265 -2401.22685 26723.56888 266 26723.56888 -2401.22685 267 16171.00545 26723.56888 268 26723.56888 16171.00545 269 -14299.80626 26723.56888 270 26723.56888 -14299.80626 271 26723.56888 26723.56888 272 26723.56888 26723.56888 273 -3108.46350 26723.56888 274 -2401.22685 -3108.46350 275 8612.68765 -2401.22685 276 8347.42108 8612.68765 277 -2401.22685 8347.42108 278 6172.76214 -2401.22685 279 -20844.16653 6172.76214 280 26723.56888 -20844.16653 281 -2401.22685 26723.56888 282 23272.10884 -2401.22685 283 -2401.22685 23272.10884 284 -24024.75157 -2401.22685 285 -2401.22685 -24024.75157 286 26723.56888 -2401.22685 287 -2486.74385 26723.56888 288 26723.56888 -2486.74385 289 16783.25672 26723.56888 290 26723.56888 16783.25672 291 7464.49902 26723.56888 292 -6947.57562 7464.49902 293 -2401.22685 -6947.57562 294 -1737.31668 -2401.22685 295 51474.88847 -1737.31668 296 -2401.22685 51474.88847 297 26723.56888 -2401.22685 298 -2401.22685 26723.56888 299 -2401.22685 -2401.22685 300 25184.63237 -2401.22685 301 -9291.67440 25184.63237 302 -2401.22685 -9291.67440 303 26723.56888 -2401.22685 304 34355.02811 26723.56888 305 -666.00381 34355.02811 306 26723.56888 -666.00381 307 22464.83173 26723.56888 308 21670.74259 22464.83173 309 26723.56888 21670.74259 310 26723.56888 26723.56888 311 -26255.79903 26723.56888 312 26723.56888 -26255.79903 313 29891.95209 26723.56888 314 -1414.38736 29891.95209 315 2902.88817 -1414.38736 316 -643.29997 2902.88817 317 -9108.09383 -643.29997 318 -1238.86183 -9108.09383 319 -2401.22685 -1238.86183 320 -2401.22685 -2401.22685 321 -2401.22685 -2401.22685 322 -2401.22685 -2401.22685 323 26723.56888 -2401.22685 324 -2401.22685 26723.56888 325 -4812.37279 -2401.22685 326 26723.56888 -4812.37279 327 26723.56888 26723.56888 328 -2401.22685 26723.56888 329 -2401.22685 -2401.22685 330 -2401.22685 -2401.22685 331 -2401.22685 -2401.22685 332 -2401.22685 -2401.22685 333 -944.13962 -2401.22685 334 -2022.08210 -944.13962 335 -2401.22685 -2022.08210 336 -2729.53365 -2401.22685 337 -2401.22685 -2729.53365 338 -2356.13091 -2401.22685 339 1358.77905 -2356.13091 340 -14305.02958 1358.77905 341 -2401.22685 -14305.02958 342 -2401.22685 -2401.22685 343 -2401.22685 -2401.22685 344 -2401.22685 -2401.22685 345 -2401.22685 -2401.22685 346 -9488.63760 -2401.22685 347 -2401.22685 -9488.63760 348 -2401.22685 -2401.22685 349 -2401.22685 -2401.22685 350 2939.06427 -2401.22685 351 -1578.88762 2939.06427 352 13632.90069 -1578.88762 353 -7486.21423 13632.90069 354 -2401.22685 -7486.21423 355 -2401.22685 -2401.22685 356 -24740.26924 -2401.22685 357 -2401.22685 -24740.26924 358 -2401.22685 -2401.22685 359 -2401.22685 -2401.22685 360 -1630.79340 -2401.22685 361 -2401.22685 -1630.79340 362 -2401.22685 -2401.22685 363 18714.84729 -2401.22685 364 -1899.99534 18714.84729 365 -2401.22685 -1899.99534 366 -2401.22685 -2401.22685 367 -2401.22685 -2401.22685 368 -2401.22685 -2401.22685 369 -2401.22685 -2401.22685 370 -2401.22685 -2401.22685 371 3486.88107 -2401.22685 372 -25749.55270 3486.88107 373 -20702.21744 -25749.55270 374 -2401.22685 -20702.21744 375 -2401.22685 -2401.22685 376 4330.91660 -2401.22685 377 -1503.97547 4330.91660 378 -2401.22685 -1503.97547 379 -2401.22685 -2401.22685 380 -2401.22685 -2401.22685 381 -3409.33390 -2401.22685 382 41667.47187 -3409.33390 383 -2401.22685 41667.47187 384 -763.76613 -2401.22685 385 -86823.08487 -763.76613 386 3994.32689 -86823.08487 387 7260.16022 3994.32689 388 -2997.41084 7260.16022 389 -9409.52872 -2997.41084 390 5753.07660 -9409.52872 391 -11857.94258 5753.07660 392 -22132.20075 -11857.94258 393 39537.25471 -22132.20075 394 -155995.44039 39537.25471 395 -17357.95088 -155995.44039 396 -25973.91882 -17357.95088 397 107289.31148 -25973.91882 398 -12279.43491 107289.31148 399 -84213.96533 -12279.43491 400 -496.61868 -84213.96533 401 12242.61718 -496.61868 402 -7008.00459 12242.61718 403 -51150.14027 -7008.00459 404 -16445.04507 -51150.14027 405 -13369.73486 -16445.04507 406 2334.62395 -13369.73486 407 4136.64604 2334.62395 408 -15614.21169 4136.64604 409 -48274.99259 -15614.21169 410 -18647.91267 -48274.99259 411 -219572.24477 -18647.91267 412 -22111.50353 -219572.24477 413 -110588.42235 -22111.50353 414 -58146.69165 -110588.42235 415 -57201.32722 -58146.69165 416 155507.03181 -57201.32722 417 -34356.48016 155507.03181 418 -14820.25975 -34356.48016 419 -46731.17993 -14820.25975 420 -85486.26190 -46731.17993 421 -82041.92819 -85486.26190 422 -627186.38066 -82041.92819 423 -331586.22377 -627186.38066 424 -68145.58823 -331586.22377 425 -197325.48460 -68145.58823 426 -251752.72344 -197325.48460 427 -251741.39555 -251752.72344 428 -127326.89068 -251741.39555 429 135682.64604 -127326.89068 430 -363758.36870 135682.64604 431 NA -363758.36870 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 3055410.06433 -335057.38899 [2,] 659410.54377 3055410.06433 [3,] -405429.22905 659410.54377 [4,] -743164.38986 -405429.22905 [5,] -89456.37662 -743164.38986 [6,] 957456.07911 -89456.37662 [7,] 358218.87563 957456.07911 [8,] -515485.62638 358218.87563 [9,] 656837.66595 -515485.62638 [10,] 418629.65698 656837.66595 [11,] 790336.11651 418629.65698 [12,] 786396.36249 790336.11651 [13,] -496204.59220 786396.36249 [14,] -198517.34823 -496204.59220 [15,] 573684.82775 -198517.34823 [16,] -252041.74304 573684.82775 [17,] -311194.45538 -252041.74304 [18,] -428313.88915 -311194.45538 [19,] -21431.49731 -428313.88915 [20,] 139295.64762 -21431.49731 [21,] -335211.15220 139295.64762 [22,] 114533.03501 -335211.15220 [23,] -174367.24805 114533.03501 [24,] -66395.70200 -174367.24805 [25,] 301730.34172 -66395.70200 [26,] -543727.48346 301730.34172 [27,] -174098.17192 -543727.48346 [28,] 593025.71736 -174098.17192 [29,] 341865.12044 593025.71736 [30,] 252881.74560 341865.12044 [31,] -132561.70348 252881.74560 [32,] 163074.51332 -132561.70348 [33,] 105181.03944 163074.51332 [34,] -175459.84008 105181.03944 [35,] 26963.16687 -175459.84008 [36,] 306078.95930 26963.16687 [37,] 161359.42712 306078.95930 [38,] 138591.64288 161359.42712 [39,] -121232.08493 138591.64288 [40,] -315987.14980 -121232.08493 [41,] 21458.74625 -315987.14980 [42,] 104529.24113 21458.74625 [43,] -48445.33886 104529.24113 [44,] 62799.44298 -48445.33886 [45,] -127194.57210 62799.44298 [46,] 118925.53513 -127194.57210 [47,] -219278.61704 118925.53513 [48,] -14451.06244 -219278.61704 [49,] -142981.86240 -14451.06244 [50,] 105768.72434 -142981.86240 [51,] -34627.40188 105768.72434 [52,] -219322.05982 -34627.40188 [53,] 70651.42675 -219322.05982 [54,] -281080.27993 70651.42675 [55,] 124875.96679 -281080.27993 [56,] 232448.43092 124875.96679 [57,] 85407.63172 232448.43092 [58,] -63984.83679 85407.63172 [59,] -220578.16036 -63984.83679 [60,] -23272.87970 -220578.16036 [61,] -78225.09783 -23272.87970 [62,] -20749.17848 -78225.09783 [63,] -216463.75382 -20749.17848 [64,] 68796.51174 -216463.75382 [65,] 23572.03863 68796.51174 [66,] -137137.24569 23572.03863 [67,] -268.85031 -137137.24569 [68,] -84555.82611 -268.85031 [69,] -55112.87409 -84555.82611 [70,] -89236.14587 -55112.87409 [71,] 153263.11954 -89236.14587 [72,] -121884.06595 153263.11954 [73,] -46107.67637 -121884.06595 [74,] -94572.99222 -46107.67637 [75,] -73086.17103 -94572.99222 [76,] -105602.61775 -73086.17103 [77,] -85644.93657 -105602.61775 [78,] 37088.85778 -85644.93657 [79,] 21253.31059 37088.85778 [80,] 29611.70914 21253.31059 [81,] -45462.15309 29611.70914 [82,] -134129.65691 -45462.15309 [83,] 36887.62176 -134129.65691 [84,] -23125.20612 36887.62176 [85,] -116597.88333 -23125.20612 [86,] 70596.47129 -116597.88333 [87,] -54748.11133 70596.47129 [88,] -88191.48526 -54748.11133 [89,] -33468.48369 -88191.48526 [90,] -20293.38577 -33468.48369 [91,] -41947.40635 -20293.38577 [92,] 122154.90850 -41947.40635 [93,] -228180.60334 122154.90850 [94,] 119037.29664 -228180.60334 [95,] -11477.33974 119037.29664 [96,] -1554.05235 -11477.33974 [97,] 49267.09848 -1554.05235 [98,] -71063.71119 49267.09848 [99,] -26793.91038 -71063.71119 [100,] -34155.89201 -26793.91038 [101,] -116393.68974 -34155.89201 [102,] -71182.11744 -116393.68974 [103,] -63711.14507 -71182.11744 [104,] 39309.00263 -63711.14507 [105,] -78901.81531 39309.00263 [106,] -233869.45225 -78901.81531 [107,] 107762.43360 -233869.45225 [108,] -87569.45515 107762.43360 [109,] -104668.69729 -87569.45515 [110,] 82970.34063 -104668.69729 [111,] -253759.38582 82970.34063 [112,] -175074.28537 -253759.38582 [113,] -121209.46808 -175074.28537 [114,] 14252.74836 -121209.46808 [115,] 55220.38255 14252.74836 [116,] -26722.55853 55220.38255 [117,] 65798.28707 -26722.55853 [118,] -176082.65391 65798.28707 [119,] 26570.83341 -176082.65391 [120,] 2477.23170 26570.83341 [121,] 13431.45499 2477.23170 [122,] -36512.52694 13431.45499 [123,] 80946.69568 -36512.52694 [124,] -100397.40476 80946.69568 [125,] -201133.25705 -100397.40476 [126,] -2401.22685 -201133.25705 [127,] -987.01213 -2401.22685 [128,] 9809.10476 -987.01213 [129,] -2401.22685 9809.10476 [130,] -386.54690 -2401.22685 [131,] 85481.18719 -386.54690 [132,] -2401.22685 85481.18719 [133,] -3753.75972 -2401.22685 [134,] -2401.22685 -3753.75972 [135,] 7599.65197 -2401.22685 [136,] 175.79713 7599.65197 [137,] -2401.22685 175.79713 [138,] -2401.22685 -2401.22685 [139,] -2401.22685 -2401.22685 [140,] 3951.31681 -2401.22685 [141,] -2401.22685 3951.31681 [142,] 17400.26177 -2401.22685 [143,] -1675.86504 17400.26177 [144,] 16223.04061 -1675.86504 [145,] 13349.51735 16223.04061 [146,] -2401.22685 13349.51735 [147,] -1714.53127 -2401.22685 [148,] -2401.22685 -1714.53127 [149,] -2564.61669 -2401.22685 [150,] 543.65038 -2564.61669 [151,] 22998.67080 543.65038 [152,] -2401.22685 22998.67080 [153,] -14612.79492 -2401.22685 [154,] -2327.64702 -14612.79492 [155,] -2401.22685 -2327.64702 [156,] -2401.22685 -2401.22685 [157,] -2399.01202 -2401.22685 [158,] 15826.16411 -2399.01202 [159,] -2401.22685 15826.16411 [160,] -2401.22685 -2401.22685 [161,] -2401.22685 -2401.22685 [162,] -2401.22685 -2401.22685 [163,] 11382.48822 -2401.22685 [164,] -60.46908 11382.48822 [165,] -2401.22685 -60.46908 [166,] -2401.22685 -2401.22685 [167,] 110015.49676 -2401.22685 [168,] -2401.22685 110015.49676 [169,] -2401.22685 -2401.22685 [170,] -2401.22685 -2401.22685 [171,] -4700.09393 -2401.22685 [172,] -2401.22685 -4700.09393 [173,] 239.82022 -2401.22685 [174,] -2401.22685 239.82022 [175,] -1522.71355 -2401.22685 [176,] -2401.22685 -1522.71355 [177,] 24853.40506 -2401.22685 [178,] -2401.22685 24853.40506 [179,] -2401.22685 -2401.22685 [180,] 39696.20723 -2401.22685 [181,] 23643.86380 39696.20723 [182,] 29378.05747 23643.86380 [183,] -12686.20265 29378.05747 [184,] 29249.09785 -12686.20265 [185,] 11165.20288 29249.09785 [186,] 33626.14463 11165.20288 [187,] -1350.71885 33626.14463 [188,] -32456.24329 -1350.71885 [189,] -2401.22685 -32456.24329 [190,] 8980.09568 -2401.22685 [191,] -2401.22685 8980.09568 [192,] 3290.42860 -2401.22685 [193,] -2401.22685 3290.42860 [194,] -2200.22282 -2401.22685 [195,] 26723.56888 -2200.22282 [196,] -41404.00191 26723.56888 [197,] 26723.56888 -41404.00191 [198,] 26723.56888 26723.56888 [199,] -13875.78642 26723.56888 [200,] 26723.56888 -13875.78642 [201,] 26723.56888 26723.56888 [202,] 26723.56888 26723.56888 [203,] 25585.75864 26723.56888 [204,] -4271.33958 25585.75864 [205,] 26784.69571 -4271.33958 [206,] 22287.18057 26784.69571 [207,] 26723.56888 22287.18057 [208,] 44234.67155 26723.56888 [209,] 26723.56888 44234.67155 [210,] -2401.22685 26723.56888 [211,] -2401.22685 -2401.22685 [212,] 11362.80316 -2401.22685 [213,] -14186.21899 11362.80316 [214,] -4637.33104 -14186.21899 [215,] -2401.22685 -4637.33104 [216,] -486.64678 -2401.22685 [217,] 6769.86053 -486.64678 [218,] -2401.22685 6769.86053 [219,] 27804.36286 -2401.22685 [220,] -27290.37324 27804.36286 [221,] -2401.22685 -27290.37324 [222,] -1963.41995 -2401.22685 [223,] -2401.22685 -1963.41995 [224,] -1935.75325 -2401.22685 [225,] 4557.43857 -1935.75325 [226,] -2401.22685 4557.43857 [227,] -2401.22685 -2401.22685 [228,] -20062.40201 -2401.22685 [229,] -2401.22685 -20062.40201 [230,] -13551.83006 -2401.22685 [231,] -2401.22685 -13551.83006 [232,] -2401.22685 -2401.22685 [233,] 6727.90438 -2401.22685 [234,] -2401.22685 6727.90438 [235,] -2333.83700 -2401.22685 [236,] -2401.22685 -2333.83700 [237,] 339.50561 -2401.22685 [238,] -2401.22685 339.50561 [239,] -2401.22685 -2401.22685 [240,] -2401.22685 -2401.22685 [241,] 3152.22517 -2401.22685 [242,] -2401.22685 3152.22517 [243,] 720.53012 -2401.22685 [244,] -7170.17477 720.53012 [245,] 26723.56888 -7170.17477 [246,] 26723.56888 26723.56888 [247,] -8564.63380 26723.56888 [248,] 26723.56888 -8564.63380 [249,] 26723.56888 26723.56888 [250,] 3156.75624 26723.56888 [251,] 26723.56888 3156.75624 [252,] -18278.38317 26723.56888 [253,] -17986.72143 -18278.38317 [254,] -42817.86848 -17986.72143 [255,] -2401.22685 -42817.86848 [256,] 26006.40867 -2401.22685 [257,] 32093.64563 26006.40867 [258,] 26723.56888 32093.64563 [259,] -2401.22685 26723.56888 [260,] 26723.56888 -2401.22685 [261,] 26723.56888 26723.56888 [262,] -9433.34743 26723.56888 [263,] 10546.87598 -9433.34743 [264,] 26723.56888 10546.87598 [265,] -2401.22685 26723.56888 [266,] 26723.56888 -2401.22685 [267,] 16171.00545 26723.56888 [268,] 26723.56888 16171.00545 [269,] -14299.80626 26723.56888 [270,] 26723.56888 -14299.80626 [271,] 26723.56888 26723.56888 [272,] 26723.56888 26723.56888 [273,] -3108.46350 26723.56888 [274,] -2401.22685 -3108.46350 [275,] 8612.68765 -2401.22685 [276,] 8347.42108 8612.68765 [277,] -2401.22685 8347.42108 [278,] 6172.76214 -2401.22685 [279,] -20844.16653 6172.76214 [280,] 26723.56888 -20844.16653 [281,] -2401.22685 26723.56888 [282,] 23272.10884 -2401.22685 [283,] -2401.22685 23272.10884 [284,] -24024.75157 -2401.22685 [285,] -2401.22685 -24024.75157 [286,] 26723.56888 -2401.22685 [287,] -2486.74385 26723.56888 [288,] 26723.56888 -2486.74385 [289,] 16783.25672 26723.56888 [290,] 26723.56888 16783.25672 [291,] 7464.49902 26723.56888 [292,] -6947.57562 7464.49902 [293,] -2401.22685 -6947.57562 [294,] -1737.31668 -2401.22685 [295,] 51474.88847 -1737.31668 [296,] -2401.22685 51474.88847 [297,] 26723.56888 -2401.22685 [298,] -2401.22685 26723.56888 [299,] -2401.22685 -2401.22685 [300,] 25184.63237 -2401.22685 [301,] -9291.67440 25184.63237 [302,] -2401.22685 -9291.67440 [303,] 26723.56888 -2401.22685 [304,] 34355.02811 26723.56888 [305,] -666.00381 34355.02811 [306,] 26723.56888 -666.00381 [307,] 22464.83173 26723.56888 [308,] 21670.74259 22464.83173 [309,] 26723.56888 21670.74259 [310,] 26723.56888 26723.56888 [311,] -26255.79903 26723.56888 [312,] 26723.56888 -26255.79903 [313,] 29891.95209 26723.56888 [314,] -1414.38736 29891.95209 [315,] 2902.88817 -1414.38736 [316,] -643.29997 2902.88817 [317,] -9108.09383 -643.29997 [318,] -1238.86183 -9108.09383 [319,] -2401.22685 -1238.86183 [320,] -2401.22685 -2401.22685 [321,] -2401.22685 -2401.22685 [322,] -2401.22685 -2401.22685 [323,] 26723.56888 -2401.22685 [324,] -2401.22685 26723.56888 [325,] -4812.37279 -2401.22685 [326,] 26723.56888 -4812.37279 [327,] 26723.56888 26723.56888 [328,] -2401.22685 26723.56888 [329,] -2401.22685 -2401.22685 [330,] -2401.22685 -2401.22685 [331,] -2401.22685 -2401.22685 [332,] -2401.22685 -2401.22685 [333,] -944.13962 -2401.22685 [334,] -2022.08210 -944.13962 [335,] -2401.22685 -2022.08210 [336,] -2729.53365 -2401.22685 [337,] -2401.22685 -2729.53365 [338,] -2356.13091 -2401.22685 [339,] 1358.77905 -2356.13091 [340,] -14305.02958 1358.77905 [341,] -2401.22685 -14305.02958 [342,] -2401.22685 -2401.22685 [343,] -2401.22685 -2401.22685 [344,] -2401.22685 -2401.22685 [345,] -2401.22685 -2401.22685 [346,] -9488.63760 -2401.22685 [347,] -2401.22685 -9488.63760 [348,] -2401.22685 -2401.22685 [349,] -2401.22685 -2401.22685 [350,] 2939.06427 -2401.22685 [351,] -1578.88762 2939.06427 [352,] 13632.90069 -1578.88762 [353,] -7486.21423 13632.90069 [354,] -2401.22685 -7486.21423 [355,] -2401.22685 -2401.22685 [356,] -24740.26924 -2401.22685 [357,] -2401.22685 -24740.26924 [358,] -2401.22685 -2401.22685 [359,] -2401.22685 -2401.22685 [360,] -1630.79340 -2401.22685 [361,] -2401.22685 -1630.79340 [362,] -2401.22685 -2401.22685 [363,] 18714.84729 -2401.22685 [364,] -1899.99534 18714.84729 [365,] -2401.22685 -1899.99534 [366,] -2401.22685 -2401.22685 [367,] -2401.22685 -2401.22685 [368,] -2401.22685 -2401.22685 [369,] -2401.22685 -2401.22685 [370,] -2401.22685 -2401.22685 [371,] 3486.88107 -2401.22685 [372,] -25749.55270 3486.88107 [373,] -20702.21744 -25749.55270 [374,] -2401.22685 -20702.21744 [375,] -2401.22685 -2401.22685 [376,] 4330.91660 -2401.22685 [377,] -1503.97547 4330.91660 [378,] -2401.22685 -1503.97547 [379,] -2401.22685 -2401.22685 [380,] -2401.22685 -2401.22685 [381,] -3409.33390 -2401.22685 [382,] 41667.47187 -3409.33390 [383,] -2401.22685 41667.47187 [384,] -763.76613 -2401.22685 [385,] -86823.08487 -763.76613 [386,] 3994.32689 -86823.08487 [387,] 7260.16022 3994.32689 [388,] -2997.41084 7260.16022 [389,] -9409.52872 -2997.41084 [390,] 5753.07660 -9409.52872 [391,] -11857.94258 5753.07660 [392,] -22132.20075 -11857.94258 [393,] 39537.25471 -22132.20075 [394,] -155995.44039 39537.25471 [395,] -17357.95088 -155995.44039 [396,] -25973.91882 -17357.95088 [397,] 107289.31148 -25973.91882 [398,] -12279.43491 107289.31148 [399,] -84213.96533 -12279.43491 [400,] -496.61868 -84213.96533 [401,] 12242.61718 -496.61868 [402,] -7008.00459 12242.61718 [403,] -51150.14027 -7008.00459 [404,] -16445.04507 -51150.14027 [405,] -13369.73486 -16445.04507 [406,] 2334.62395 -13369.73486 [407,] 4136.64604 2334.62395 [408,] -15614.21169 4136.64604 [409,] -48274.99259 -15614.21169 [410,] -18647.91267 -48274.99259 [411,] -219572.24477 -18647.91267 [412,] -22111.50353 -219572.24477 [413,] -110588.42235 -22111.50353 [414,] -58146.69165 -110588.42235 [415,] -57201.32722 -58146.69165 [416,] 155507.03181 -57201.32722 [417,] -34356.48016 155507.03181 [418,] -14820.25975 -34356.48016 [419,] -46731.17993 -14820.25975 [420,] -85486.26190 -46731.17993 [421,] -82041.92819 -85486.26190 [422,] -627186.38066 -82041.92819 [423,] -331586.22377 -627186.38066 [424,] -68145.58823 -331586.22377 [425,] -197325.48460 -68145.58823 [426,] -251752.72344 -197325.48460 [427,] -251741.39555 -251752.72344 [428,] -127326.89068 -251741.39555 [429,] 135682.64604 -127326.89068 [430,] -363758.36870 135682.64604 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 3055410.06433 -335057.38899 2 659410.54377 3055410.06433 3 -405429.22905 659410.54377 4 -743164.38986 -405429.22905 5 -89456.37662 -743164.38986 6 957456.07911 -89456.37662 7 358218.87563 957456.07911 8 -515485.62638 358218.87563 9 656837.66595 -515485.62638 10 418629.65698 656837.66595 11 790336.11651 418629.65698 12 786396.36249 790336.11651 13 -496204.59220 786396.36249 14 -198517.34823 -496204.59220 15 573684.82775 -198517.34823 16 -252041.74304 573684.82775 17 -311194.45538 -252041.74304 18 -428313.88915 -311194.45538 19 -21431.49731 -428313.88915 20 139295.64762 -21431.49731 21 -335211.15220 139295.64762 22 114533.03501 -335211.15220 23 -174367.24805 114533.03501 24 -66395.70200 -174367.24805 25 301730.34172 -66395.70200 26 -543727.48346 301730.34172 27 -174098.17192 -543727.48346 28 593025.71736 -174098.17192 29 341865.12044 593025.71736 30 252881.74560 341865.12044 31 -132561.70348 252881.74560 32 163074.51332 -132561.70348 33 105181.03944 163074.51332 34 -175459.84008 105181.03944 35 26963.16687 -175459.84008 36 306078.95930 26963.16687 37 161359.42712 306078.95930 38 138591.64288 161359.42712 39 -121232.08493 138591.64288 40 -315987.14980 -121232.08493 41 21458.74625 -315987.14980 42 104529.24113 21458.74625 43 -48445.33886 104529.24113 44 62799.44298 -48445.33886 45 -127194.57210 62799.44298 46 118925.53513 -127194.57210 47 -219278.61704 118925.53513 48 -14451.06244 -219278.61704 49 -142981.86240 -14451.06244 50 105768.72434 -142981.86240 51 -34627.40188 105768.72434 52 -219322.05982 -34627.40188 53 70651.42675 -219322.05982 54 -281080.27993 70651.42675 55 124875.96679 -281080.27993 56 232448.43092 124875.96679 57 85407.63172 232448.43092 58 -63984.83679 85407.63172 59 -220578.16036 -63984.83679 60 -23272.87970 -220578.16036 61 -78225.09783 -23272.87970 62 -20749.17848 -78225.09783 63 -216463.75382 -20749.17848 64 68796.51174 -216463.75382 65 23572.03863 68796.51174 66 -137137.24569 23572.03863 67 -268.85031 -137137.24569 68 -84555.82611 -268.85031 69 -55112.87409 -84555.82611 70 -89236.14587 -55112.87409 71 153263.11954 -89236.14587 72 -121884.06595 153263.11954 73 -46107.67637 -121884.06595 74 -94572.99222 -46107.67637 75 -73086.17103 -94572.99222 76 -105602.61775 -73086.17103 77 -85644.93657 -105602.61775 78 37088.85778 -85644.93657 79 21253.31059 37088.85778 80 29611.70914 21253.31059 81 -45462.15309 29611.70914 82 -134129.65691 -45462.15309 83 36887.62176 -134129.65691 84 -23125.20612 36887.62176 85 -116597.88333 -23125.20612 86 70596.47129 -116597.88333 87 -54748.11133 70596.47129 88 -88191.48526 -54748.11133 89 -33468.48369 -88191.48526 90 -20293.38577 -33468.48369 91 -41947.40635 -20293.38577 92 122154.90850 -41947.40635 93 -228180.60334 122154.90850 94 119037.29664 -228180.60334 95 -11477.33974 119037.29664 96 -1554.05235 -11477.33974 97 49267.09848 -1554.05235 98 -71063.71119 49267.09848 99 -26793.91038 -71063.71119 100 -34155.89201 -26793.91038 101 -116393.68974 -34155.89201 102 -71182.11744 -116393.68974 103 -63711.14507 -71182.11744 104 39309.00263 -63711.14507 105 -78901.81531 39309.00263 106 -233869.45225 -78901.81531 107 107762.43360 -233869.45225 108 -87569.45515 107762.43360 109 -104668.69729 -87569.45515 110 82970.34063 -104668.69729 111 -253759.38582 82970.34063 112 -175074.28537 -253759.38582 113 -121209.46808 -175074.28537 114 14252.74836 -121209.46808 115 55220.38255 14252.74836 116 -26722.55853 55220.38255 117 65798.28707 -26722.55853 118 -176082.65391 65798.28707 119 26570.83341 -176082.65391 120 2477.23170 26570.83341 121 13431.45499 2477.23170 122 -36512.52694 13431.45499 123 80946.69568 -36512.52694 124 -100397.40476 80946.69568 125 -201133.25705 -100397.40476 126 -2401.22685 -201133.25705 127 -987.01213 -2401.22685 128 9809.10476 -987.01213 129 -2401.22685 9809.10476 130 -386.54690 -2401.22685 131 85481.18719 -386.54690 132 -2401.22685 85481.18719 133 -3753.75972 -2401.22685 134 -2401.22685 -3753.75972 135 7599.65197 -2401.22685 136 175.79713 7599.65197 137 -2401.22685 175.79713 138 -2401.22685 -2401.22685 139 -2401.22685 -2401.22685 140 3951.31681 -2401.22685 141 -2401.22685 3951.31681 142 17400.26177 -2401.22685 143 -1675.86504 17400.26177 144 16223.04061 -1675.86504 145 13349.51735 16223.04061 146 -2401.22685 13349.51735 147 -1714.53127 -2401.22685 148 -2401.22685 -1714.53127 149 -2564.61669 -2401.22685 150 543.65038 -2564.61669 151 22998.67080 543.65038 152 -2401.22685 22998.67080 153 -14612.79492 -2401.22685 154 -2327.64702 -14612.79492 155 -2401.22685 -2327.64702 156 -2401.22685 -2401.22685 157 -2399.01202 -2401.22685 158 15826.16411 -2399.01202 159 -2401.22685 15826.16411 160 -2401.22685 -2401.22685 161 -2401.22685 -2401.22685 162 -2401.22685 -2401.22685 163 11382.48822 -2401.22685 164 -60.46908 11382.48822 165 -2401.22685 -60.46908 166 -2401.22685 -2401.22685 167 110015.49676 -2401.22685 168 -2401.22685 110015.49676 169 -2401.22685 -2401.22685 170 -2401.22685 -2401.22685 171 -4700.09393 -2401.22685 172 -2401.22685 -4700.09393 173 239.82022 -2401.22685 174 -2401.22685 239.82022 175 -1522.71355 -2401.22685 176 -2401.22685 -1522.71355 177 24853.40506 -2401.22685 178 -2401.22685 24853.40506 179 -2401.22685 -2401.22685 180 39696.20723 -2401.22685 181 23643.86380 39696.20723 182 29378.05747 23643.86380 183 -12686.20265 29378.05747 184 29249.09785 -12686.20265 185 11165.20288 29249.09785 186 33626.14463 11165.20288 187 -1350.71885 33626.14463 188 -32456.24329 -1350.71885 189 -2401.22685 -32456.24329 190 8980.09568 -2401.22685 191 -2401.22685 8980.09568 192 3290.42860 -2401.22685 193 -2401.22685 3290.42860 194 -2200.22282 -2401.22685 195 26723.56888 -2200.22282 196 -41404.00191 26723.56888 197 26723.56888 -41404.00191 198 26723.56888 26723.56888 199 -13875.78642 26723.56888 200 26723.56888 -13875.78642 201 26723.56888 26723.56888 202 26723.56888 26723.56888 203 25585.75864 26723.56888 204 -4271.33958 25585.75864 205 26784.69571 -4271.33958 206 22287.18057 26784.69571 207 26723.56888 22287.18057 208 44234.67155 26723.56888 209 26723.56888 44234.67155 210 -2401.22685 26723.56888 211 -2401.22685 -2401.22685 212 11362.80316 -2401.22685 213 -14186.21899 11362.80316 214 -4637.33104 -14186.21899 215 -2401.22685 -4637.33104 216 -486.64678 -2401.22685 217 6769.86053 -486.64678 218 -2401.22685 6769.86053 219 27804.36286 -2401.22685 220 -27290.37324 27804.36286 221 -2401.22685 -27290.37324 222 -1963.41995 -2401.22685 223 -2401.22685 -1963.41995 224 -1935.75325 -2401.22685 225 4557.43857 -1935.75325 226 -2401.22685 4557.43857 227 -2401.22685 -2401.22685 228 -20062.40201 -2401.22685 229 -2401.22685 -20062.40201 230 -13551.83006 -2401.22685 231 -2401.22685 -13551.83006 232 -2401.22685 -2401.22685 233 6727.90438 -2401.22685 234 -2401.22685 6727.90438 235 -2333.83700 -2401.22685 236 -2401.22685 -2333.83700 237 339.50561 -2401.22685 238 -2401.22685 339.50561 239 -2401.22685 -2401.22685 240 -2401.22685 -2401.22685 241 3152.22517 -2401.22685 242 -2401.22685 3152.22517 243 720.53012 -2401.22685 244 -7170.17477 720.53012 245 26723.56888 -7170.17477 246 26723.56888 26723.56888 247 -8564.63380 26723.56888 248 26723.56888 -8564.63380 249 26723.56888 26723.56888 250 3156.75624 26723.56888 251 26723.56888 3156.75624 252 -18278.38317 26723.56888 253 -17986.72143 -18278.38317 254 -42817.86848 -17986.72143 255 -2401.22685 -42817.86848 256 26006.40867 -2401.22685 257 32093.64563 26006.40867 258 26723.56888 32093.64563 259 -2401.22685 26723.56888 260 26723.56888 -2401.22685 261 26723.56888 26723.56888 262 -9433.34743 26723.56888 263 10546.87598 -9433.34743 264 26723.56888 10546.87598 265 -2401.22685 26723.56888 266 26723.56888 -2401.22685 267 16171.00545 26723.56888 268 26723.56888 16171.00545 269 -14299.80626 26723.56888 270 26723.56888 -14299.80626 271 26723.56888 26723.56888 272 26723.56888 26723.56888 273 -3108.46350 26723.56888 274 -2401.22685 -3108.46350 275 8612.68765 -2401.22685 276 8347.42108 8612.68765 277 -2401.22685 8347.42108 278 6172.76214 -2401.22685 279 -20844.16653 6172.76214 280 26723.56888 -20844.16653 281 -2401.22685 26723.56888 282 23272.10884 -2401.22685 283 -2401.22685 23272.10884 284 -24024.75157 -2401.22685 285 -2401.22685 -24024.75157 286 26723.56888 -2401.22685 287 -2486.74385 26723.56888 288 26723.56888 -2486.74385 289 16783.25672 26723.56888 290 26723.56888 16783.25672 291 7464.49902 26723.56888 292 -6947.57562 7464.49902 293 -2401.22685 -6947.57562 294 -1737.31668 -2401.22685 295 51474.88847 -1737.31668 296 -2401.22685 51474.88847 297 26723.56888 -2401.22685 298 -2401.22685 26723.56888 299 -2401.22685 -2401.22685 300 25184.63237 -2401.22685 301 -9291.67440 25184.63237 302 -2401.22685 -9291.67440 303 26723.56888 -2401.22685 304 34355.02811 26723.56888 305 -666.00381 34355.02811 306 26723.56888 -666.00381 307 22464.83173 26723.56888 308 21670.74259 22464.83173 309 26723.56888 21670.74259 310 26723.56888 26723.56888 311 -26255.79903 26723.56888 312 26723.56888 -26255.79903 313 29891.95209 26723.56888 314 -1414.38736 29891.95209 315 2902.88817 -1414.38736 316 -643.29997 2902.88817 317 -9108.09383 -643.29997 318 -1238.86183 -9108.09383 319 -2401.22685 -1238.86183 320 -2401.22685 -2401.22685 321 -2401.22685 -2401.22685 322 -2401.22685 -2401.22685 323 26723.56888 -2401.22685 324 -2401.22685 26723.56888 325 -4812.37279 -2401.22685 326 26723.56888 -4812.37279 327 26723.56888 26723.56888 328 -2401.22685 26723.56888 329 -2401.22685 -2401.22685 330 -2401.22685 -2401.22685 331 -2401.22685 -2401.22685 332 -2401.22685 -2401.22685 333 -944.13962 -2401.22685 334 -2022.08210 -944.13962 335 -2401.22685 -2022.08210 336 -2729.53365 -2401.22685 337 -2401.22685 -2729.53365 338 -2356.13091 -2401.22685 339 1358.77905 -2356.13091 340 -14305.02958 1358.77905 341 -2401.22685 -14305.02958 342 -2401.22685 -2401.22685 343 -2401.22685 -2401.22685 344 -2401.22685 -2401.22685 345 -2401.22685 -2401.22685 346 -9488.63760 -2401.22685 347 -2401.22685 -9488.63760 348 -2401.22685 -2401.22685 349 -2401.22685 -2401.22685 350 2939.06427 -2401.22685 351 -1578.88762 2939.06427 352 13632.90069 -1578.88762 353 -7486.21423 13632.90069 354 -2401.22685 -7486.21423 355 -2401.22685 -2401.22685 356 -24740.26924 -2401.22685 357 -2401.22685 -24740.26924 358 -2401.22685 -2401.22685 359 -2401.22685 -2401.22685 360 -1630.79340 -2401.22685 361 -2401.22685 -1630.79340 362 -2401.22685 -2401.22685 363 18714.84729 -2401.22685 364 -1899.99534 18714.84729 365 -2401.22685 -1899.99534 366 -2401.22685 -2401.22685 367 -2401.22685 -2401.22685 368 -2401.22685 -2401.22685 369 -2401.22685 -2401.22685 370 -2401.22685 -2401.22685 371 3486.88107 -2401.22685 372 -25749.55270 3486.88107 373 -20702.21744 -25749.55270 374 -2401.22685 -20702.21744 375 -2401.22685 -2401.22685 376 4330.91660 -2401.22685 377 -1503.97547 4330.91660 378 -2401.22685 -1503.97547 379 -2401.22685 -2401.22685 380 -2401.22685 -2401.22685 381 -3409.33390 -2401.22685 382 41667.47187 -3409.33390 383 -2401.22685 41667.47187 384 -763.76613 -2401.22685 385 -86823.08487 -763.76613 386 3994.32689 -86823.08487 387 7260.16022 3994.32689 388 -2997.41084 7260.16022 389 -9409.52872 -2997.41084 390 5753.07660 -9409.52872 391 -11857.94258 5753.07660 392 -22132.20075 -11857.94258 393 39537.25471 -22132.20075 394 -155995.44039 39537.25471 395 -17357.95088 -155995.44039 396 -25973.91882 -17357.95088 397 107289.31148 -25973.91882 398 -12279.43491 107289.31148 399 -84213.96533 -12279.43491 400 -496.61868 -84213.96533 401 12242.61718 -496.61868 402 -7008.00459 12242.61718 403 -51150.14027 -7008.00459 404 -16445.04507 -51150.14027 405 -13369.73486 -16445.04507 406 2334.62395 -13369.73486 407 4136.64604 2334.62395 408 -15614.21169 4136.64604 409 -48274.99259 -15614.21169 410 -18647.91267 -48274.99259 411 -219572.24477 -18647.91267 412 -22111.50353 -219572.24477 413 -110588.42235 -22111.50353 414 -58146.69165 -110588.42235 415 -57201.32722 -58146.69165 416 155507.03181 -57201.32722 417 -34356.48016 155507.03181 418 -14820.25975 -34356.48016 419 -46731.17993 -14820.25975 420 -85486.26190 -46731.17993 421 -82041.92819 -85486.26190 422 -627186.38066 -82041.92819 423 -331586.22377 -627186.38066 424 -68145.58823 -331586.22377 425 -197325.48460 -68145.58823 426 -251752.72344 -197325.48460 427 -251741.39555 -251752.72344 428 -127326.89068 -251741.39555 429 135682.64604 -127326.89068 430 -363758.36870 135682.64604 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7i9rl1291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/8tj861291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/9tj861291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/freestat/rcomp/tmp/104a781291330094.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/117a6e1291330094.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/12st4k1291330094.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/137l2b1291330094.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/14s3jz1291330094.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/15v4h51291330094.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/162nh81291330095.tab") + } > > try(system("convert tmp/1frax1291330094.ps tmp/1frax1291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/2frax1291330094.ps tmp/2frax1291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/3frax1291330094.ps tmp/3frax1291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/4q0si1291330094.ps tmp/4q0si1291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/5q0si1291330094.ps tmp/5q0si1291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/6q0si1291330094.ps tmp/6q0si1291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/7i9rl1291330094.ps tmp/7i9rl1291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/8tj861291330094.ps tmp/8tj861291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/9tj861291330094.ps tmp/9tj861291330094.png",intern=TRUE)) character(0) > try(system("convert tmp/104a781291330094.ps tmp/104a781291330094.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 17.837 3.053 18.343