R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1 + ,1595 + ,17 + ,60720 + ,319369 + ,0 + ,0 + ,5565 + ,47 + ,94355 + ,493408 + ,39 + ,0 + ,601 + ,6 + ,60720 + ,319210 + ,0 + ,1 + ,188 + ,11 + ,77655 + ,381180 + ,0 + ,1 + ,7146 + ,105 + ,134028 + ,297978 + ,305 + ,0 + ,1135 + ,17 + ,62285 + ,290476 + ,-9 + ,1 + ,450 + ,8 + ,59325 + ,292136 + ,-1 + ,1 + ,34 + ,8 + ,60630 + ,314353 + ,0 + ,1 + ,133 + ,14 + ,65990 + ,339445 + ,0 + ,1 + ,119 + ,6 + ,59118 + ,303677 + ,0 + ,0 + ,2053 + ,53 + ,100423 + ,397144 + ,33 + ,0 + ,4036 + ,63 + ,100269 + ,424898 + ,32 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4655 + ,393 + ,60720 + ,341570 + ,0 + ,1 + ,131 + ,11 + ,61808 + ,308989 + ,0 + ,1 + ,1766 + ,73 + ,60438 + ,305959 + ,0 + ,1 + ,312 + ,14 + ,58598 + ,318690 + ,0 + ,1 + ,448 + ,40 + ,59781 + ,323361 + ,0 + ,1 + ,115 + ,13 + ,60945 + ,318903 + ,0 + ,1 + ,60 + ,10 + ,61124 + ,314049 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,364 + ,35 + ,47705 + ,298466 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,1442 + ,118 + ,121173 + ,438493 + ,-3 + ,0 + ,1389 + ,9 + ,57530 + ,378049 + ,6 + ,1 + ,149 + ,5 + ,62041 + ,313332 + ,0 + ,1 + ,2212 + ,14 + ,60720 + ,319864 + ,0 + ,0 + ,7489 + ,322 + ,136996 + ,430866 + ,70 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,402 + ,26 + ,84990 + ,332743 + ,5 + ,0 + ,7419 + ,29 + ,63255 + ,286963 + ,-5 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,307 + ,15 + ,58856 + ,331955 + ,0 + ,0 + ,1134 + ,162 + ,146216 + ,527021 + ,58 + ,0 + ,7561 + ,87 + ,82425 + ,364304 + ,37 + ,1 + ,5131 + ,71 + ,86111 + ,292154 + ,9 + ,1 + ,9 + ,6 + ,60835 + ,314987 + ,0 + ,1 + ,2264 + ,24 + ,55174 + ,272713 + ,-31 + ,1 + ,40949 + ,481 + ,129352 + ,1398893 + ,2865 + ,0 + ,4980 + ,91 + ,113521 + ,409642 + ,0 + ,1 + ,272 + ,8 + ,112995 + ,429112 + ,14 + ,1 + ,757 + ,19 + ,45689 + ,382712 + ,0 + ,1 + ,1424 + ,94 + ,131116 + ,374943 + ,32 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,378 + ,52 + ,71873 + ,297413 + ,0 + ,1 + ,1061 + ,5 + ,70415 + ,304252 + ,-4 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,0 + ,3689 + ,25 + ,101792 + ,377305 + ,74 + ,1 + ,567 + ,11 + ,103345 + ,386212 + ,20 + ,0 + ,27664 + ,282 + ,74996 + ,657954 + ,247 + ,1 + ,1686 + ,40 + ,68696 + ,368186 + ,11 + ,0 + ,1164 + ,17 + ,57320 + ,269753 + ,-6 + ,1 + ,15824 + ,64 + ,63838 + ,585715 + ,284 + ,1 + ,6575 + ,159 + ,151352 + ,1071292 + ,819 + ,0 + ,37597 + ,248 + ,37884 + ,992426 + ,652 + ,1 + ,314 + ,5 + ,60720 + ,315371 + ,0 + ,0 + ,1932 + ,36 + ,72954 + ,249898 + ,10 + ,0 + ,5692 + ,72 + ,64239 + ,357312 + ,22 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,6533 + ,75 + ,126910 + ,364839 + ,16 + ,0 + ,2126 + ,138 + ,21001 + ,230621 + ,-3 + ,1 + ,14 + ,1 + ,61016 + ,315877 + ,0 + ,0 + ,3192 + ,60 + ,59831 + ,125390 + ,10 + ,1 + ,47 + ,3 + ,60720 + ,314882 + ,0 + ,0 + ,4891 + ,217 + ,210568 + ,370837 + ,73 + ,0 + ,531 + ,10 + ,71641 + ,330068 + ,0 + ,1 + ,2332 + ,4 + ,72680 + ,324385 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,16275 + ,109 + ,47818 + ,1443586 + ,1974 + ,1 + ,9252 + ,48 + ,79510 + ,253537 + ,-129 + ,0 + ,29790 + ,309 + ,81767 + ,4321023 + ,28328 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,2300 + ,29 + ,64873 + ,352108 + ,0 + ,1 + ,267 + ,32 + ,57640 + ,330059 + ,0 + ,0 + ,382 + ,21 + ,59365 + ,340968 + ,0 + ,1 + ,397 + ,28 + ,60646 + ,317736 + ,0 + ,1 + ,11920 + ,182 + ,31080 + ,209458 + ,238 + ,0 + ,54660 + ,302 + ,126942 + ,1491348 + ,3040 + ,1 + ,5 + ,2 + ,60720 + ,314887 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,2647 + ,37 + ,62920 + ,353058 + ,1 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,94 + ,7 + ,60793 + ,314533 + ,0 + ,1 + ,3 + ,2 + ,60698 + ,315354 + ,0 + ,1 + ,454 + ,51 + ,63261 + ,302187 + ,0 + ,1 + ,227 + ,9 + ,76644 + ,336639 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,1 + ,25830 + ,115 + ,245546 + ,1073089 + ,1910 + ,1 + ,1528 + ,38 + ,60326 + ,146494 + ,38 + ,1 + ,271 + ,2 + ,60720 + ,325249 + ,0 + ,1 + ,138 + ,6 + ,69817 + ,331420 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,278 + ,26 + ,61404 + ,314922 + ,0 + ,1 + ,282 + ,30 + ,62452 + ,320016 + ,0 + ,1 + ,571 + ,80 + ,41477 + ,280398 + ,0 + ,1 + ,2253 + ,99 + ,63593 + ,452469 + ,17 + ,1 + ,290 + ,2 + ,58790 + ,301164 + ,0 + ,1 + ,78 + ,16 + ,62700 + ,317330 + ,0 + ,1 + ,20 + ,2 + ,60805 + ,315576 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,18073 + ,180 + ,27284 + ,-7170 + ,-366 + ,1 + ,1866 + ,29 + ,56225 + ,322331 + ,0 + ,0 + ,42634 + ,156 + ,157214 + ,1629616 + ,2964 + ,1 + ,249 + ,62 + ,54323 + ,292754 + ,0 + ,1 + ,422 + ,28 + ,57935 + ,318056 + ,0 + ,1 + ,2675 + ,79 + ,59017 + ,355178 + ,36 + ,1 + ,965 + ,3 + ,73490 + ,204325 + ,-8 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,621 + ,16 + ,68005 + ,317046 + ,-5 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,365 + ,8 + ,54820 + ,309560 + ,0 + ,0 + ,3122 + ,34 + ,94670 + ,414462 + ,38 + ,0 + ,12988 + ,119 + ,82340 + ,857217 + ,601 + ,0 + ,5336 + ,81 + ,112477 + ,697458 + ,132 + ,0 + ,3160 + ,108 + ,108094 + ,530670 + ,66 + ,1 + ,2489 + ,40 + ,67804 + ,238125 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,6984 + ,45 + ,80570 + ,741409 + ,195 + ,1 + ,2001 + ,52 + ,95551 + ,393343 + ,19 + ,1 + ,15236 + ,90 + ,77440 + ,372631 + ,-158 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,530 + ,60 + ,73433 + ,317291 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,35624 + ,223 + ,157278 + ,306275 + ,-70 + ,0 + ,1383 + ,11 + ,73221 + ,317892 + ,13 + ,1 + ,875 + ,14 + ,67000 + ,334280 + ,3 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,72 + ,3 + ,60398 + ,314210 + ,0 + ,0 + ,265 + ,31 + ,60720 + ,306948 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,335 + ,33 + ,64175 + ,320398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4525 + ,226 + ,93811 + ,501749 + ,101 + ,0 + ,3045 + ,58 + ,27330 + ,202055 + ,31 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,638 + ,14 + ,60370 + ,333210 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,607 + ,11 + ,60436 + ,322340 + ,0 + ,1 + ,1558 + ,63 + ,55637 + ,369448 + ,0 + ,0 + ,1324 + ,117 + ,67440 + ,291841 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,611 + ,8 + ,59190 + ,296919 + ,-4 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,923 + ,8 + ,58620 + ,309038 + ,0 + ,1 + ,661 + ,3 + ,65920 + ,246541 + ,-12 + ,0 + ,2397 + ,14 + ,74020 + ,289513 + ,6 + ,1 + ,366 + ,53 + ,61808 + ,344425 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,135 + ,2 + ,62065 + ,314210 + ,0 + ,0 + ,1659 + ,29 + ,107577 + ,480382 + ,2 + ,1 + ,316 + ,9 + ,60505 + ,315009 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,309 + ,11 + ,56535 + ,312878 + ,0 + ,0 + ,49 + ,8 + ,64107 + ,322031 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,4519 + ,49 + ,102129 + ,597793 + ,163 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,837 + ,69 + ,61262 + ,315688 + ,4 + ,1 + ,5119 + ,49 + ,39039 + ,378525 + ,31 + ,1 + ,1280 + ,117 + ,69465 + ,312378 + ,26 + ,1 + ,2564 + ,22 + ,130140 + ,403560 + ,13 + ,0 + ,2045 + ,20 + ,97890 + ,510834 + ,10 + ,0 + ,2234 + ,16 + ,77200 + ,214215 + ,5 + ,1 + ,975 + ,23 + ,90534 + ,235133 + ,13 + ,0 + ,1136 + ,32 + ,48522 + ,343613 + ,0 + ,1 + ,453 + ,21 + ,81125 + ,365959 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,61 + ,14 + ,60720 + ,314551 + ,0 + ,1 + ,368 + ,17 + ,61686 + ,303230 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4901 + ,84 + ,121920 + ,469107 + ,67 + ,1 + ,540 + ,14 + ,103960 + ,354228 + ,22 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,2 + ,9 + ,60735 + ,315394 + ,0 + ,1 + ,36 + ,4 + ,61564 + ,312412 + ,0 + ,1 + ,776 + ,9 + ,64230 + ,333505 + ,-4 + ,1 + ,84738 + ,588 + ,-26007 + ,223193 + ,4145 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,3 + ,4 + ,60761 + ,315656 + ,0 + ,1 + ,529 + ,9 + ,63870 + ,296261 + ,-3 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,405 + ,7 + ,60845 + ,336425 + ,0 + ,1 + ,972 + ,34 + ,71642 + ,359335 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,2099 + ,64 + ,106611 + ,308636 + ,20 + ,0 + ,3437 + ,47 + ,48022 + ,158492 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,22330 + ,84 + ,79801 + ,711969 + ,446 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,483 + ,27 + ,60830 + ,306268 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,2239 + ,21 + ,88590 + ,442882 + ,5 + ,0 + ,2949 + ,41 + ,82903 + ,378509 + ,18 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,365 + ,83 + ,87192 + ,346611 + ,0 + ,1 + ,2461 + ,57 + ,55792 + ,314289 + ,0 + ,0 + ,21950 + ,519 + ,114337 + ,856956 + ,745 + ,0 + ,3294 + ,23 + ,75832 + ,217193 + ,1 + ,1 + ,141 + ,6 + ,61630 + ,315366 + ,0 + ,1 + ,572 + ,16 + ,58580 + ,307930 + ,-1 + ,1 + ,13326 + ,102 + ,165548 + ,702380 + ,436 + ,1 + ,2284 + ,33 + ,76403 + ,194493 + ,66 + ,0 + ,10 + ,2 + ,61656 + ,316155 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,1414 + ,198 + ,70184 + ,330546 + ,13 + ,1 + ,1975 + ,35 + ,118881 + ,394510 + ,30 + ,1 + ,43 + ,3 + ,60887 + ,312846 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,844 + ,78 + ,60925 + ,296139 + ,0 + ,1 + ,304 + ,15 + ,62969 + ,295580 + ,0 + ,1 + ,458 + ,11 + ,58625 + ,297765 + ,0 + ,1 + ,18562 + ,155 + ,102313 + ,377934 + ,143 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,7123 + ,109 + ,288170 + ,638830 + ,769 + ,1 + ,622 + ,73 + ,73007 + ,304376 + ,0 + ,1 + ,174 + ,1 + ,64820 + ,307424 + ,0 + ,1 + ,2220 + ,22 + ,301670 + ,644190 + ,154 + ,1 + ,121 + ,24 + ,56178 + ,295370 + ,0 + ,0 + ,11819 + ,85 + ,106113 + ,574339 + ,199 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,125 + ,12 + ,60798 + ,310201 + ,0 + ,1 + ,1182 + ,13 + ,70694 + ,327007 + ,-5 + ,1 + ,1503 + ,39 + ,56364 + ,343466 + ,6 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,30 + ,4 + ,62045 + ,318098 + ,0 + ,0 + ,3310 + ,33 + ,75230 + ,448243 + ,-6 + ,1 + ,554 + ,11 + ,79285 + ,325738 + ,30 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,468 + ,21 + ,60720 + ,312161 + ,0 + ,1 + ,4917 + ,42 + ,52811 + ,243650 + ,-3 + ,1 + ,3256 + ,66 + ,35250 + ,407159 + ,-8 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,125 + ,13 + ,59734 + ,317698 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,22 + ,2 + ,60722 + ,312502 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,514 + ,52 + ,78780 + ,322378 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,10327 + ,91 + ,88577 + ,640273 + ,135 + ,1 + ,13718 + ,61 + ,96448 + ,345783 + ,104 + ,1 + ,3748 + ,9 + ,50350 + ,652925 + ,53 + ,1 + ,14416 + ,361 + ,49857 + ,439798 + ,220 + ,0 + ,1526 + ,25 + ,69351 + ,278990 + ,-4 + ,0 + ,666 + ,171 + ,117869 + ,339836 + ,25 + ,1 + ,2844 + ,67 + ,72683 + ,240897 + ,-15 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,368 + ,17 + ,61167 + ,297141 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,333 + ,13 + ,70811 + ,331323 + ,0 + ,1 + ,26 + ,1 + ,60896 + ,313880 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1303 + ,61 + ,69863 + ,309422 + ,-5 + ,1 + ,20 + ,3 + ,60938 + ,315245 + ,0 + ,1 + ,2384 + ,97 + ,61348 + ,405972 + ,54 + ,0 + ,203 + ,27 + ,50804 + ,300962 + ,0 + ,0 + ,71 + ,1 + ,60745 + ,316176 + ,0 + ,1 + ,53 + ,14 + ,59506 + ,302409 + ,0 + ,1 + ,562 + ,13 + ,58456 + ,283587 + ,-3 + ,1 + ,622 + ,7 + ,60950 + ,263276 + ,-2 + ,1 + ,645 + ,5 + ,60720 + ,312075 + ,0 + ,0 + ,1763 + ,24 + ,61600 + ,308336 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,317 + ,9 + ,63915 + ,298700 + ,0 + ,1 + ,1 + ,6 + ,60719 + ,315372 + ,0 + ,1 + ,275 + ,4 + ,59500 + ,318745 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,936 + ,30 + ,67939 + ,408881 + ,9 + ,1 + ,8568 + ,37 + ,32168 + ,786690 + ,223 + ,0 + ,11528 + ,70 + ,-14545 + ,-83265 + ,-25 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,738 + ,10 + ,60720 + ,321376 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,592 + ,7 + ,64270 + ,276898 + ,-4 + ,1 + ,126 + ,10 + ,60951 + ,315547 + ,0 + ,1 + ,2 + ,1 + ,60743 + ,315487 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,18014 + ,239 + ,-1710 + ,1405225 + ,943 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,37704 + ,452 + ,60448 + ,983660 + ,2167 + ,0 + ,63 + ,4 + ,65688 + ,318574 + ,0 + ,1 + ,1431 + ,34 + ,106885 + ,310768 + ,30 + ,1 + ,94 + ,3 + ,61360 + ,312887 + ,0 + ,1 + ,192 + ,4 + ,65276 + ,312339 + ,0 + ,1 + ,32 + ,1 + ,59988 + ,314964 + ,0 + ,1 + ,6869 + ,105 + ,117520 + ,379983 + ,90 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,2328 + ,94 + ,82732 + ,253588 + ,0 + ,0 + ,209 + ,13 + ,64016 + ,316647 + ,0 + ,1 + ,28 + ,2 + ,60890 + ,315688 + ,0 + ,1 + ,176 + ,28 + ,68136 + ,310670 + ,0 + ,1 + ,1920 + ,15 + ,79420 + ,165404 + ,-6 + ,0 + ,87550 + ,458 + ,153198 + ,4111912 + ,17936 + ,1 + ,520 + ,10 + ,58650 + ,291650 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1013 + ,2 + ,69770 + ,253468 + ,-14 + ,1 + ,15 + ,8 + ,60831 + ,315688 + ,0 + ,1 + ,587 + ,12 + ,59595 + ,325699 + ,0 + ,1 + ,5371 + ,66 + ,87720 + ,446211 + ,49 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1012 + ,47 + ,114768 + ,368078 + ,26 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,876 + ,59 + ,138971 + ,352850 + ,24 + ,0 + ,162556 + ,1081 + ,213118 + ,6282154 + ,56289 + ,1 + ,43556 + ,247 + ,32648 + ,227132 + ,929 + ,0 + ,3425 + ,12 + ,83620 + ,283910 + ,-3 + ,1 + ,810 + ,43 + ,74015 + ,236761 + ,26 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,2365 + ,79 + ,191778 + ,550608 + ,120 + ,0 + ,1261 + ,33 + ,76114 + ,307528 + ,13 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1585 + ,107 + ,93099 + ,355864 + ,2 + ,1 + ,16189 + ,295 + ,116384 + ,358589 + ,379 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,10579 + ,92 + ,110309 + ,375195 + ,4 + ,1 + ,474 + ,20 + ,61977 + ,288985 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,3642 + ,86 + ,90262 + ,458343 + ,22 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,472 + ,20 + ,80045 + ,269587 + ,0 + ,1 + ,98 + ,2 + ,61490 + ,315236 + ,0 + ,1 + ,3999 + ,30 + ,51252 + ,42754 + ,-58 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,621 + ,48 + ,60720 + ,308256 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,746 + ,23 + ,71561 + ,269661 + ,10 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,4150 + ,150 + ,134759 + ,518365 + ,66 + ,0 + ,4658 + ,158 + ,156608 + ,233773 + ,58 + ,1 + ,814 + ,18 + ,39625 + ,301881 + ,23 + ,1 + ,1002 + ,15 + ,56750 + ,298568 + ,0 + ,0 + ,496 + ,17 + ,87390 + ,325479 + ,0 + ,1 + ,389 + ,14 + ,58990 + ,325506 + ,0 + ,1 + ,12679 + ,110 + ,48020 + ,984885 + ,1871 + ,1 + ,400 + ,152 + ,60720 + ,313267 + ,0 + ,1 + ,53 + ,6 + ,60349 + ,315793 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,5109 + ,97 + ,67038 + ,215362 + ,-31 + ,1 + ,576 + ,18 + ,113761 + ,314073 + ,12 + ,1 + ,438 + ,11 + ,58320 + ,298096 + ,-5 + ,1 + ,165 + ,24 + ,62841 + ,325176 + ,0 + ,1 + ,4069 + ,145 + ,79804 + ,207393 + ,119 + ,0 + ,23 + ,5 + ,62555 + ,314806 + ,0 + ,1 + ,1285 + ,107 + ,99489 + ,341340 + ,49 + ,0 + ,4677 + ,75 + ,90131 + ,426280 + ,24 + ,0 + ,24811 + ,505 + ,95350 + ,929118 + ,822 + ,1 + ,167 + ,2 + ,64245 + ,307322 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4628 + ,62 + ,69159 + ,387475 + ,85 + ,1 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,0 + ,1765 + ,20 + ,77623 + ,247060 + ,28 + ,1 + ,460 + ,15 + ,63346 + ,329784 + ,0 + ,1 + ,36 + ,3 + ,60894 + ,315834 + ,0 + ,1 + ,989 + ,11 + ,58930 + ,304555 + ,-5 + ,0 + ,3055 + ,43 + ,60247 + ,376641 + ,14 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,13253 + ,947 + ,90829 + ,357760 + ,35 + ,1 + ,24 + ,8 + ,59818 + ,315637 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,131 + ,6 + ,59661 + ,327071 + ,0 + ,1 + ,514 + ,15 + ,102725 + ,377516 + ,23 + ,1 + ,3366 + ,80 + ,108479 + ,299243 + ,12 + ,0 + ,9327 + ,101 + ,87419 + ,387699 + ,44 + ,1 + ,384 + ,22 + ,54683 + ,309836 + ,0 + ,1 + ,17821 + ,373 + ,122844 + ,444477 + ,643 + ,1 + ,397 + ,10 + ,62710 + ,322327 + ,-1 + ,1 + ,897 + ,5 + ,60720 + ,314913 + ,0 + ,1 + ,218 + ,21 + ,60720 + ,315553 + ,0 + ,0 + ,3369 + ,28 + ,87161 + ,688779 + ,217 + ,1 + ,5702 + ,61 + ,101481 + ,321896 + ,10 + ,1 + ,8636 + ,855 + ,128294 + ,301607 + ,108 + ,1 + ,534 + ,7 + ,62620 + ,304485 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,726 + ,3 + ,69980 + ,231861 + ,-5 + ,1 + ,1380 + ,26 + ,60982 + ,347385 + ,25 + ,1 + ,180 + ,11 + ,59635 + ,316386 + ,0 + ,0 + ,7285 + ,115 + ,188873 + ,491303 + ,59 + ,1 + ,880 + ,74 + ,80791 + ,261216 + ,34 + ,1 + ,1 + ,2 + ,60727 + ,315388 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,96 + ,5 + ,60379 + ,313729 + ,0 + ,1 + ,1889 + ,45 + ,37527 + ,358649 + ,14 + ,0 + ,45187 + ,353 + ,234817 + ,1926517 + ,2854 + ,1 + ,288 + ,4 + ,60510 + ,296656 + ,0 + ,1 + ,1270 + ,26 + ,69206 + ,275311 + ,-10 + ,1 + ,6526 + ,26 + ,55830 + ,-42143 + ,-28 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,226 + ,12 + ,72835 + ,343929 + ,0 + ,1 + ,694 + ,8 + ,68060 + ,367655 + ,-2 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,249 + ,72 + ,56726 + ,313491 + ,0) + ,dim=c(6 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Trades' + ,'Dividends' + ,'Wealth' + ,'CCScore ') + ,1:431)) > y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','Trades','Dividends','Wealth','CCScore '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '6' > ylab = '' > xlab = '' > main = '' > #'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 > 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 CCScore\r Group Costs Trades Dividends Wealth 1 0 1 1595 17 60720 319369 2 39 0 5565 47 94355 493408 3 0 0 601 6 60720 319210 4 0 1 188 11 77655 381180 5 305 1 7146 105 134028 297978 6 -9 0 1135 17 62285 290476 7 -1 1 450 8 59325 292136 8 0 1 34 8 60630 314353 9 0 1 133 14 65990 339445 10 0 1 119 6 59118 303677 11 33 0 2053 53 100423 397144 12 32 0 4036 63 100269 424898 13 0 1 0 0 60720 315380 14 0 1 4655 393 60720 341570 15 0 1 131 11 61808 308989 16 0 1 1766 73 60438 305959 17 0 1 312 14 58598 318690 18 0 1 448 40 59781 323361 19 0 1 115 13 60945 318903 20 0 1 60 10 61124 314049 21 0 1 0 0 60720 315380 22 0 1 364 35 47705 298466 23 0 1 0 0 60720 315380 24 -3 0 1442 118 121173 438493 25 6 0 1389 9 57530 378049 26 0 1 149 5 62041 313332 27 0 1 2212 14 60720 319864 28 70 0 7489 322 136996 430866 29 0 1 0 0 60720 315380 30 5 1 402 26 84990 332743 31 -5 0 7419 29 63255 286963 32 0 1 0 0 60720 315380 33 0 1 307 15 58856 331955 34 58 0 1134 162 146216 527021 35 37 0 7561 87 82425 364304 36 9 1 5131 71 86111 292154 37 0 1 9 6 60835 314987 38 -31 1 2264 24 55174 272713 39 2865 1 40949 481 129352 1398893 40 0 0 4980 91 113521 409642 41 14 1 272 8 112995 429112 42 0 1 757 19 45689 382712 43 32 1 1424 94 131116 374943 44 0 1 0 0 60720 315380 45 0 0 378 52 71873 297413 46 -4 1 1061 5 70415 304252 47 0 1 0 0 60720 315380 48 0 1 203 9 61938 309596 49 74 0 3689 25 101792 377305 50 20 1 567 11 103345 386212 51 247 0 27664 282 74996 657954 52 11 1 1686 40 68696 368186 53 -6 0 1164 17 57320 269753 54 284 1 15824 64 63838 585715 55 819 1 6575 159 151352 1071292 56 652 0 37597 248 37884 992426 57 0 1 314 5 60720 315371 58 10 0 1932 36 72954 249898 59 22 0 5692 72 64239 357312 60 0 1 0 0 60720 315380 61 0 1 0 0 60720 315380 62 16 0 6533 75 126910 364839 63 -3 0 2126 138 21001 230621 64 0 1 14 1 61016 315877 65 10 0 3192 60 59831 125390 66 0 1 47 3 60720 314882 67 73 0 4891 217 210568 370837 68 0 0 531 10 71641 330068 69 0 1 2332 4 72680 324385 70 0 1 0 0 60720 315380 71 1974 0 16275 109 47818 1443586 72 -129 1 9252 48 79510 253537 73 28328 0 29790 309 81767 4321023 74 0 1 0 0 60720 315380 75 0 0 2300 29 64873 352108 76 0 1 267 32 57640 330059 77 0 0 382 21 59365 340968 78 0 1 397 28 60646 317736 79 238 1 11920 182 31080 209458 80 3040 0 54660 302 126942 1491348 81 0 1 5 2 60720 314887 82 0 1 0 0 60720 315380 83 0 1 0 0 60720 315380 84 1 0 2647 37 62920 353058 85 0 1 0 0 60720 315380 86 0 0 0 0 60720 315380 87 0 1 94 7 60793 314533 88 0 1 3 2 60698 315354 89 0 1 454 51 63261 302187 90 0 1 227 9 76644 336639 91 0 0 0 0 60720 315380 92 0 1 206 9 53110 296702 93 1910 1 25830 115 245546 1073089 94 38 1 1528 38 60326 146494 95 0 1 271 2 60720 325249 96 0 1 138 6 69817 331420 97 0 1 0 0 60720 315380 98 0 1 278 26 61404 314922 99 0 1 282 30 62452 320016 100 0 1 571 80 41477 280398 101 17 1 2253 99 63593 452469 102 0 1 290 2 58790 301164 103 0 1 78 16 62700 317330 104 0 1 20 2 60805 315576 105 0 1 0 0 60720 315380 106 -366 1 18073 180 27284 -7170 107 0 1 1866 29 56225 322331 108 2964 0 42634 156 157214 1629616 109 0 1 249 62 54323 292754 110 0 1 422 28 57935 318056 111 36 1 2675 79 59017 355178 112 -8 1 965 3 73490 204325 113 0 1 0 0 60720 315380 114 -5 1 621 16 68005 317046 115 0 0 0 0 60720 315380 116 0 1 365 8 54820 309560 117 38 0 3122 34 94670 414462 118 601 0 12988 119 82340 857217 119 132 0 5336 81 112477 697458 120 66 0 3160 108 108094 530670 121 0 1 2489 40 67804 238125 122 0 1 0 0 60720 315380 123 195 0 6984 45 80570 741409 124 19 1 2001 52 95551 393343 125 -158 1 15236 90 77440 372631 126 0 0 0 0 60720 315380 127 0 1 530 60 73433 317291 128 0 1 0 0 60720 315380 129 -70 1 35624 223 157278 306275 130 13 0 1383 11 73221 317892 131 3 1 875 14 67000 334280 132 0 0 0 0 60720 315380 133 0 0 0 0 60720 315380 134 0 0 0 0 60720 315380 135 0 1 72 3 60398 314210 136 1 0 265 31 60720 306948 137 1 0 0 60720 315380 0 138 0 335 33 64175 320398 0 139 1 0 0 60720 315380 0 140 1 0 0 60720 315380 0 141 0 4525 226 93811 501749 101 142 1 3045 58 27330 202055 31 143 1 0 0 60720 315380 0 144 0 0 0 60720 315380 0 145 0 638 14 60370 333210 0 146 1 0 0 60720 315380 0 147 1 607 11 60436 322340 0 148 0 1558 63 55637 369448 0 149 1 1324 117 67440 291841 0 150 0 0 0 60720 315380 0 151 1 0 0 60720 315380 0 152 1 611 8 59190 296919 -4 153 1 0 0 60720 315380 0 154 1 923 8 58620 309038 0 155 0 661 3 65920 246541 -12 156 1 2397 14 74020 289513 6 157 1 366 53 61808 344425 0 158 1 0 0 60720 315380 0 159 0 135 2 62065 314210 0 160 1 1659 29 107577 480382 2 161 1 316 9 60505 315009 0 162 1 0 0 60720 315380 0 163 0 309 11 56535 312878 0 164 1 49 8 64107 322031 0 165 0 0 0 60720 315380 0 166 1 4519 49 102129 597793 163 167 1 0 0 60720 315380 0 168 1 837 69 61262 315688 4 169 1 5119 49 39039 378525 31 170 1 1280 117 69465 312378 26 171 0 2564 22 130140 403560 13 172 0 2045 20 97890 510834 10 173 1 2234 16 77200 214215 5 174 0 975 23 90534 235133 13 175 1 1136 32 48522 343613 0 176 1 453 21 81125 365959 0 177 1 0 0 60720 315380 0 178 1 61 14 60720 314551 0 179 0 368 17 61686 303230 0 180 1 0 0 60720 315380 0 181 1 4901 84 121920 469107 67 182 1 540 14 103960 354228 22 183 1 0 0 60720 315380 0 184 0 0 0 60720 315380 0 185 1 0 0 60720 315380 0 186 1 2 9 60735 315394 0 187 1 36 4 61564 312412 0 188 1 776 9 64230 333505 -4 189 1 84738 588 -26007 223193 4145 190 1 0 0 60720 315380 0 191 1 3 4 60761 315656 0 192 1 529 9 63870 296261 -3 193 1 0 0 60720 315380 0 194 1 405 7 60845 336425 0 195 1 972 34 71642 359335 0 196 1 0 0 60720 315380 0 197 0 2099 64 106611 308636 20 198 0 3437 47 48022 158492 9 199 1 0 0 60720 315380 0 200 0 0 0 60720 315380 0 201 0 22330 84 79801 711969 446 202 0 0 0 60720 315380 0 203 1 0 0 60720 315380 0 204 1 483 27 60830 306268 0 205 1 0 0 60720 315380 0 206 0 2239 21 88590 442882 5 207 1 2949 41 82903 378509 18 208 1 0 0 60720 315380 0 209 1 365 83 87192 346611 0 210 0 2461 57 55792 314289 0 211 0 21950 519 114337 856956 745 212 1 3294 23 75832 217193 1 213 1 141 6 61630 315366 0 214 1 572 16 58580 307930 -1 215 1 13326 102 165548 702380 436 216 0 2284 33 76403 194493 66 217 1 10 2 61656 316155 0 218 0 0 0 60720 315380 0 219 1 1414 198 70184 330546 13 220 1 1975 35 118881 394510 30 221 1 43 3 60887 312846 0 222 1 0 0 60720 315380 0 223 1 844 78 60925 296139 0 224 1 304 15 62969 295580 0 225 1 458 11 58625 297765 0 226 1 18562 155 102313 377934 143 227 1 0 0 60720 315380 0 228 1 7123 109 288170 638830 769 229 1 622 73 73007 304376 0 230 1 174 1 64820 307424 0 231 1 2220 22 301670 644190 154 232 0 121 24 56178 295370 0 233 1 11819 85 106113 574339 199 234 1 0 0 60720 315380 0 235 1 125 12 60798 310201 0 236 1 1182 13 70694 327007 -5 237 1 1503 39 56364 343466 6 238 1 0 0 60720 315380 0 239 1 0 0 60720 315380 0 240 0 30 4 62045 318098 0 241 1 3310 33 75230 448243 -6 242 0 554 11 79285 325738 30 243 1 0 0 60720 315380 0 244 1 468 21 60720 312161 0 245 1 4917 42 52811 243650 -3 246 1 3256 66 35250 407159 -8 247 1 0 0 60720 315380 0 248 1 125 13 59734 317698 0 249 1 0 0 60720 315380 0 250 0 22 2 60722 312502 0 251 1 0 0 60720 315380 0 252 0 514 52 78780 322378 0 253 1 0 0 60720 315380 0 254 0 0 0 60720 315380 0 255 1 0 0 60720 315380 0 256 1 0 0 60720 315380 0 257 1 0 0 60720 315380 0 258 1 0 0 60720 315380 0 259 1 10327 91 88577 640273 135 260 1 13718 61 96448 345783 104 261 1 3748 9 50350 652925 53 262 0 14416 361 49857 439798 220 263 0 1526 25 69351 278990 -4 264 1 666 171 117869 339836 25 265 1 2844 67 72683 240897 -15 266 1 0 0 60720 315380 0 267 1 368 17 61167 297141 0 268 1 0 0 60720 315380 0 269 1 333 13 70811 331323 0 270 1 26 1 60896 313880 0 271 1 0 0 60720 315380 0 272 1 0 0 60720 315380 0 273 1 1303 61 69863 309422 -5 274 1 20 3 60938 315245 0 275 0 2384 97 61348 405972 54 276 0 203 27 50804 300962 0 277 1 71 1 60745 316176 0 278 1 53 14 59506 302409 0 279 1 562 13 58456 283587 -3 280 1 622 7 60950 263276 -2 281 0 645 5 60720 312075 0 282 1 1763 24 61600 308336 0 283 1 0 0 60720 315380 0 284 1 317 9 63915 298700 0 285 1 1 6 60719 315372 0 286 1 275 4 59500 318745 0 287 1 0 0 60720 315380 0 288 1 936 30 67939 408881 9 289 0 8568 37 32168 786690 223 290 0 11528 70 -14545 -83265 -25 291 1 0 0 60720 315380 0 292 1 738 10 60720 321376 0 293 1 0 0 60720 315380 0 294 1 592 7 64270 276898 -4 295 1 126 10 60951 315547 0 296 1 2 1 60743 315487 0 297 0 0 0 60720 315380 0 298 0 18014 239 -1710 1405225 943 299 0 0 0 60720 315380 0 300 0 0 0 60720 315380 0 301 0 37704 452 60448 983660 2167 302 1 63 4 65688 318574 0 303 1 1431 34 106885 310768 30 304 1 94 3 61360 312887 0 305 1 192 4 65276 312339 0 306 1 32 1 59988 314964 0 307 1 6869 105 117520 379983 90 308 1 0 0 60720 315380 0 309 1 0 0 60720 315380 0 310 0 1 4 60722 315398 0 311 1 0 0 60720 315380 0 312 0 2328 94 82732 253588 0 313 1 209 13 64016 316647 0 314 1 28 2 60890 315688 0 315 1 176 28 68136 310670 0 316 0 1920 15 79420 165404 -6 317 1 87550 458 153198 4111912 17936 318 0 520 10 58650 291650 0 319 1 0 0 60720 315380 0 320 1 1013 2 69770 253468 -14 321 1 15 8 60831 315688 0 322 1 587 12 59595 325699 0 323 1 5371 66 87720 446211 49 324 1 0 0 60720 315380 0 325 1 1012 47 114768 368078 26 326 1 0 0 60720 315380 0 327 0 876 59 138971 352850 24 328 1 162556 1081 213118 6282154 56289 329 0 43556 247 32648 227132 929 330 1 3425 12 83620 283910 -3 331 1 810 43 74015 236761 26 332 0 0 0 60720 315380 0 333 0 0 0 60720 315380 0 334 0 2365 79 191778 550608 120 335 1 1261 33 76114 307528 13 336 1 0 0 60720 315380 0 337 1 1585 107 93099 355864 2 338 1 16189 295 116384 358589 379 339 0 0 0 60720 315380 0 340 1 0 0 60720 315380 0 341 1 10579 92 110309 375195 4 342 1 474 20 61977 288985 0 343 1 0 0 60720 315380 0 344 1 3642 86 90262 458343 22 345 0 0 0 60720 315380 0 346 1 472 20 80045 269587 0 347 1 98 2 61490 315236 0 348 1 3999 30 51252 42754 -58 349 1 0 0 60720 315380 0 350 1 621 48 60720 308256 0 351 1 0 0 60720 315380 0 352 0 0 0 60720 315380 0 353 1 30 4 60798 313164 0 354 1 0 0 60720 315380 0 355 0 746 23 71561 269661 10 356 1 0 0 60720 315380 0 357 0 0 0 60720 315380 0 358 1 0 0 60720 315380 0 359 0 0 0 60720 315380 0 360 0 4150 150 134759 518365 66 361 1 4658 158 156608 233773 58 362 1 814 18 39625 301881 23 363 0 1002 15 56750 298568 0 364 1 496 17 87390 325479 0 365 1 389 14 58990 325506 0 366 1 12679 110 48020 984885 1871 367 1 400 152 60720 313267 0 368 1 53 6 60349 315793 0 369 0 0 0 60720 315380 0 370 1 5109 97 67038 215362 -31 371 1 576 18 113761 314073 12 372 1 438 11 58320 298096 -5 373 1 165 24 62841 325176 0 374 0 4069 145 79804 207393 119 375 1 23 5 62555 314806 0 376 0 1285 107 99489 341340 49 377 0 4677 75 90131 426280 24 378 1 24811 505 95350 929118 822 379 1 167 2 64245 307322 0 380 1 0 0 60720 315380 0 381 1 4628 62 69159 387475 85 382 0 226 5 65745 291787 0 383 1 1765 20 77623 247060 28 384 1 460 15 63346 329784 0 385 1 36 3 60894 315834 0 386 0 989 11 58930 304555 -5 387 0 3055 43 60247 376641 14 388 1 0 0 60720 315380 0 389 0 0 0 60720 315380 0 390 1 0 0 60720 315380 0 391 1 0 0 60720 315380 0 392 1 0 0 60720 315380 0 393 1 13253 947 90829 357760 35 394 1 24 8 59818 315637 0 395 0 0 0 60720 315380 0 396 1 0 0 60720 315380 0 397 1 0 0 60720 315380 0 398 1 0 0 60720 315380 0 399 1 131 6 59661 327071 0 400 1 514 15 102725 377516 23 401 0 3366 80 108479 299243 12 402 1 9327 101 87419 387699 44 403 1 384 22 54683 309836 0 404 1 17821 373 122844 444477 643 405 1 397 10 62710 322327 -1 406 1 897 5 60720 314913 0 407 0 218 21 60720 315553 0 408 1 3369 28 87161 688779 217 409 1 5702 61 101481 321896 10 410 1 8636 855 128294 301607 108 411 1 534 7 62620 304485 0 412 0 0 0 60720 315380 0 413 1 0 0 60720 315380 0 414 1 726 3 69980 231861 -5 415 1 1380 26 60982 347385 25 416 0 180 11 59635 316386 0 417 1 7285 115 188873 491303 59 418 1 880 74 80791 261216 34 419 1 1 2 60727 315388 0 420 1 0 0 60720 315380 0 421 1 96 5 60379 313729 0 422 0 1889 45 37527 358649 14 423 1 45187 353 234817 1926517 2854 424 1 288 4 60510 296656 0 425 1 1270 26 69206 275311 -10 426 0 6526 26 55830 -42143 -28 427 1 0 0 60720 315380 0 428 1 226 12 72835 343929 0 429 1 694 8 68060 367655 -2 430 1 0 0 60720 315380 0 431 1 249 72 56726 313491 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Dividends Wealth -1.185e+03 -1.050e-02 -6.235e-02 1.372e-02 1.202e-04 5.227e-03 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -3378.6 -428.6 247.2 315.4 8771.3 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -1.185e+03 7.036e+01 -16.846 < 2e-16 *** Group -1.050e-02 6.360e-03 -1.651 0.0995 . Costs -6.235e-02 8.697e-03 -7.170 3.35e-12 *** Trades 1.372e-02 1.131e-03 12.128 < 2e-16 *** Dividends 1.202e-04 1.966e-04 0.611 0.5413 Wealth 5.227e-03 1.755e-04 29.790 < 2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 729.1 on 425 degrees of freedom Multiple R-squared: 0.729, Adjusted R-squared: 0.7258 F-statistic: 228.6 on 5 and 425 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.677346e-05 3.354693e-05 9.999832e-01 [2,] 3.451887e-07 6.903774e-07 9.999997e-01 [3,] 1.078783e-08 2.157566e-08 1.000000e+00 [4,] 2.509298e-10 5.018595e-10 1.000000e+00 [5,] 4.963145e-12 9.926291e-12 1.000000e+00 [6,] 9.039728e-13 1.807946e-12 1.000000e+00 [7,] 1.938604e-14 3.877208e-14 1.000000e+00 [8,] 7.691098e-16 1.538220e-15 1.000000e+00 [9,] 2.015417e-17 4.030833e-17 1.000000e+00 [10,] 5.345730e-19 1.069146e-18 1.000000e+00 [11,] 1.163285e-20 2.326571e-20 1.000000e+00 [12,] 2.300441e-22 4.600882e-22 1.000000e+00 [13,] 4.488800e-24 8.977599e-24 1.000000e+00 [14,] 3.298070e-25 6.596140e-25 1.000000e+00 [15,] 6.454870e-27 1.290974e-26 1.000000e+00 [16,] 4.879198e-28 9.758395e-28 1.000000e+00 [17,] 1.189891e-27 2.379781e-27 1.000000e+00 [18,] 3.035039e-29 6.070078e-29 1.000000e+00 [19,] 8.414731e-30 1.682946e-29 1.000000e+00 [20,] 4.047575e-30 8.095150e-30 1.000000e+00 [21,] 1.212880e-31 2.425760e-31 1.000000e+00 [22,] 6.472157e-33 1.294431e-32 1.000000e+00 [23,] 6.714550e-31 1.342910e-30 1.000000e+00 [24,] 2.696599e-32 5.393198e-32 1.000000e+00 [25,] 1.118978e-33 2.237957e-33 1.000000e+00 [26,] 5.807141e-35 1.161428e-34 1.000000e+00 [27,] 2.688470e-36 5.376939e-36 1.000000e+00 [28,] 9.209374e-36 1.841875e-35 1.000000e+00 [29,] 4.032605e-37 8.065210e-37 1.000000e+00 [30,] 5.127136e-38 1.025427e-37 1.000000e+00 [31,] 9.240865e-28 1.848173e-27 1.000000e+00 [32,] 1.665227e-28 3.330454e-28 1.000000e+00 [33,] 5.400613e-29 1.080123e-28 1.000000e+00 [34,] 1.689924e-29 3.379849e-29 1.000000e+00 [35,] 1.781298e-30 3.562596e-30 1.000000e+00 [36,] 1.869708e-31 3.739416e-31 1.000000e+00 [37,] 7.357510e-31 1.471502e-30 1.000000e+00 [38,] 7.924913e-32 1.584983e-31 1.000000e+00 [39,] 8.775655e-33 1.755131e-32 1.000000e+00 [40,] 9.692653e-34 1.938531e-33 1.000000e+00 [41,] 1.234657e-34 2.469314e-34 1.000000e+00 [42,] 1.602581e-35 3.205162e-35 1.000000e+00 [43,] 2.748169e-25 5.496338e-25 1.000000e+00 [44,] 7.024019e-26 1.404804e-25 1.000000e+00 [45,] 6.052559e-26 1.210512e-25 1.000000e+00 [46,] 6.074398e-25 1.214880e-24 1.000000e+00 [47,] 5.207021e-23 1.041404e-22 1.000000e+00 [48,] 1.457187e-20 2.914374e-20 1.000000e+00 [49,] 3.944977e-21 7.889954e-21 1.000000e+00 [50,] 2.362866e-21 4.725732e-21 1.000000e+00 [51,] 6.645107e-22 1.329021e-21 1.000000e+00 [52,] 1.802758e-22 3.605516e-22 1.000000e+00 [53,] 4.863480e-23 9.726961e-23 1.000000e+00 [54,] 1.258580e-23 2.517160e-23 1.000000e+00 [55,] 6.884702e-24 1.376940e-23 1.000000e+00 [56,] 1.846280e-24 3.692560e-24 1.000000e+00 [57,] 2.616074e-24 5.232149e-24 1.000000e+00 [58,] 7.130378e-25 1.426076e-24 1.000000e+00 [59,] 1.928661e-25 3.857323e-25 1.000000e+00 [60,] 6.478062e-26 1.295612e-25 1.000000e+00 [61,] 1.721957e-26 3.443913e-26 1.000000e+00 [62,] 4.646390e-27 9.292781e-27 1.000000e+00 [63,] 1.300471e-25 2.600943e-25 1.000000e+00 [64,] 3.647360e-26 7.294719e-26 1.000000e+00 [65,] 1.000000e+00 2.762772e-146 1.381386e-146 [66,] 1.000000e+00 1.779885e-145 8.899427e-146 [67,] 1.000000e+00 5.727666e-145 2.863833e-145 [68,] 1.000000e+00 3.144854e-144 1.572427e-144 [69,] 1.000000e+00 1.257698e-143 6.288489e-144 [70,] 1.000000e+00 7.241033e-143 3.620516e-143 [71,] 1.000000e+00 4.231106e-146 2.115553e-146 [72,] 1.000000e+00 8.348916e-153 4.174458e-153 [73,] 1.000000e+00 5.101400e-152 2.550700e-152 [74,] 1.000000e+00 3.070059e-151 1.535030e-151 [75,] 1.000000e+00 1.821170e-150 9.105848e-151 [76,] 1.000000e+00 5.135912e-150 2.567956e-150 [77,] 1.000000e+00 2.943525e-149 1.471763e-149 [78,] 1.000000e+00 1.611638e-148 8.058191e-149 [79,] 1.000000e+00 8.874433e-148 4.437216e-148 [80,] 1.000000e+00 4.846123e-147 2.423062e-147 [81,] 1.000000e+00 2.551639e-146 1.275820e-146 [82,] 1.000000e+00 1.182848e-145 5.914241e-146 [83,] 1.000000e+00 5.960437e-145 2.980218e-145 [84,] 1.000000e+00 3.160320e-144 1.580160e-144 [85,] 1.000000e+00 1.810311e-153 9.051554e-154 [86,] 1.000000e+00 1.632360e-154 8.161801e-155 [87,] 1.000000e+00 8.211026e-154 4.105513e-154 [88,] 1.000000e+00 4.007385e-153 2.003693e-153 [89,] 1.000000e+00 2.120810e-152 1.060405e-152 [90,] 1.000000e+00 1.062634e-151 5.313168e-152 [91,] 1.000000e+00 5.031077e-151 2.515538e-151 [92,] 1.000000e+00 2.727164e-150 1.363582e-150 [93,] 1.000000e+00 5.755644e-151 2.877822e-151 [94,] 1.000000e+00 2.726361e-150 1.363180e-150 [95,] 1.000000e+00 1.203205e-149 6.016024e-150 [96,] 1.000000e+00 5.254052e-149 2.627026e-149 [97,] 1.000000e+00 2.231271e-148 1.115635e-148 [98,] 1.000000e+00 5.965881e-163 2.982940e-163 [99,] 1.000000e+00 1.995278e-162 9.976392e-163 [100,] 1.000000e+00 0.000000e+00 0.000000e+00 [101,] 1.000000e+00 0.000000e+00 0.000000e+00 [102,] 1.000000e+00 9.881313e-324 4.940656e-324 [103,] 1.000000e+00 1.284571e-322 6.422853e-323 [104,] 1.000000e+00 2.667954e-322 1.333977e-322 [105,] 1.000000e+00 1.333977e-321 6.669886e-322 [106,] 1.000000e+00 4.555285e-321 2.277643e-321 [107,] 1.000000e+00 1.747016e-320 8.735081e-321 [108,] 1.000000e+00 6.540441e-320 3.270221e-320 [109,] 1.000000e+00 5.744995e-319 2.872498e-319 [110,] 1.000000e+00 0.000000e+00 0.000000e+00 [111,] 1.000000e+00 0.000000e+00 0.000000e+00 [112,] 1.000000e+00 0.000000e+00 0.000000e+00 [113,] 1.000000e+00 0.000000e+00 0.000000e+00 [114,] 1.000000e+00 0.000000e+00 0.000000e+00 [115,] 1.000000e+00 0.000000e+00 0.000000e+00 [116,] 1.000000e+00 0.000000e+00 0.000000e+00 [117,] 1.000000e+00 0.000000e+00 0.000000e+00 [118,] 1.000000e+00 0.000000e+00 0.000000e+00 [119,] 1.000000e+00 0.000000e+00 0.000000e+00 [120,] 1.000000e+00 0.000000e+00 0.000000e+00 [121,] 1.000000e+00 0.000000e+00 0.000000e+00 [122,] 1.000000e+00 0.000000e+00 0.000000e+00 [123,] 1.000000e+00 0.000000e+00 0.000000e+00 [124,] 1.000000e+00 0.000000e+00 0.000000e+00 [125,] 1.000000e+00 0.000000e+00 0.000000e+00 [126,] 1.000000e+00 0.000000e+00 0.000000e+00 [127,] 1.000000e+00 0.000000e+00 0.000000e+00 [128,] 1.000000e+00 0.000000e+00 0.000000e+00 [129,] 1.000000e+00 0.000000e+00 0.000000e+00 [130,] 1.000000e+00 0.000000e+00 0.000000e+00 [131,] 1.000000e+00 0.000000e+00 0.000000e+00 [132,] 1.000000e+00 0.000000e+00 0.000000e+00 [133,] 1.000000e+00 0.000000e+00 0.000000e+00 [134,] 1.000000e+00 0.000000e+00 0.000000e+00 [135,] 1.000000e+00 0.000000e+00 0.000000e+00 [136,] 1.000000e+00 0.000000e+00 0.000000e+00 [137,] 1.000000e+00 0.000000e+00 0.000000e+00 [138,] 1.000000e+00 0.000000e+00 0.000000e+00 [139,] 1.000000e+00 0.000000e+00 0.000000e+00 [140,] 1.000000e+00 0.000000e+00 0.000000e+00 [141,] 1.000000e+00 0.000000e+00 0.000000e+00 [142,] 1.000000e+00 0.000000e+00 0.000000e+00 [143,] 1.000000e+00 0.000000e+00 0.000000e+00 [144,] 1.000000e+00 0.000000e+00 0.000000e+00 [145,] 1.000000e+00 0.000000e+00 0.000000e+00 [146,] 1.000000e+00 0.000000e+00 0.000000e+00 [147,] 1.000000e+00 0.000000e+00 0.000000e+00 [148,] 1.000000e+00 0.000000e+00 0.000000e+00 [149,] 1.000000e+00 0.000000e+00 0.000000e+00 [150,] 1.000000e+00 0.000000e+00 0.000000e+00 [151,] 1.000000e+00 0.000000e+00 0.000000e+00 [152,] 1.000000e+00 0.000000e+00 0.000000e+00 [153,] 1.000000e+00 0.000000e+00 0.000000e+00 [154,] 1.000000e+00 0.000000e+00 0.000000e+00 [155,] 1.000000e+00 0.000000e+00 0.000000e+00 [156,] 1.000000e+00 0.000000e+00 0.000000e+00 [157,] 1.000000e+00 0.000000e+00 0.000000e+00 [158,] 1.000000e+00 0.000000e+00 0.000000e+00 [159,] 1.000000e+00 0.000000e+00 0.000000e+00 [160,] 1.000000e+00 0.000000e+00 0.000000e+00 [161,] 1.000000e+00 0.000000e+00 0.000000e+00 [162,] 1.000000e+00 0.000000e+00 0.000000e+00 [163,] 1.000000e+00 0.000000e+00 0.000000e+00 [164,] 1.000000e+00 0.000000e+00 0.000000e+00 [165,] 1.000000e+00 0.000000e+00 0.000000e+00 [166,] 1.000000e+00 0.000000e+00 0.000000e+00 [167,] 1.000000e+00 0.000000e+00 0.000000e+00 [168,] 1.000000e+00 0.000000e+00 0.000000e+00 [169,] 1.000000e+00 0.000000e+00 0.000000e+00 [170,] 1.000000e+00 0.000000e+00 0.000000e+00 [171,] 1.000000e+00 0.000000e+00 0.000000e+00 [172,] 1.000000e+00 0.000000e+00 0.000000e+00 [173,] 1.000000e+00 0.000000e+00 0.000000e+00 [174,] 1.000000e+00 0.000000e+00 0.000000e+00 [175,] 1.000000e+00 0.000000e+00 0.000000e+00 [176,] 1.000000e+00 0.000000e+00 0.000000e+00 [177,] 1.000000e+00 0.000000e+00 0.000000e+00 [178,] 1.000000e+00 0.000000e+00 0.000000e+00 [179,] 1.000000e+00 0.000000e+00 0.000000e+00 [180,] 1.000000e+00 0.000000e+00 0.000000e+00 [181,] 1.000000e+00 0.000000e+00 0.000000e+00 [182,] 1.000000e+00 0.000000e+00 0.000000e+00 [183,] 1.000000e+00 0.000000e+00 0.000000e+00 [184,] 1.000000e+00 0.000000e+00 0.000000e+00 [185,] 1.000000e+00 0.000000e+00 0.000000e+00 [186,] 1.000000e+00 0.000000e+00 0.000000e+00 [187,] 1.000000e+00 0.000000e+00 0.000000e+00 [188,] 1.000000e+00 0.000000e+00 0.000000e+00 [189,] 1.000000e+00 0.000000e+00 0.000000e+00 [190,] 1.000000e+00 0.000000e+00 0.000000e+00 [191,] 1.000000e+00 0.000000e+00 0.000000e+00 [192,] 1.000000e+00 0.000000e+00 0.000000e+00 [193,] 1.000000e+00 0.000000e+00 0.000000e+00 [194,] 1.000000e+00 0.000000e+00 0.000000e+00 [195,] 1.000000e+00 0.000000e+00 0.000000e+00 [196,] 1.000000e+00 0.000000e+00 0.000000e+00 [197,] 1.000000e+00 0.000000e+00 0.000000e+00 [198,] 1.000000e+00 0.000000e+00 0.000000e+00 [199,] 1.000000e+00 0.000000e+00 0.000000e+00 [200,] 1.000000e+00 0.000000e+00 0.000000e+00 [201,] 1.000000e+00 0.000000e+00 0.000000e+00 [202,] 1.000000e+00 0.000000e+00 0.000000e+00 [203,] 1.000000e+00 0.000000e+00 0.000000e+00 [204,] 1.000000e+00 0.000000e+00 0.000000e+00 [205,] 1.000000e+00 0.000000e+00 0.000000e+00 [206,] 1.000000e+00 0.000000e+00 0.000000e+00 [207,] 1.000000e+00 0.000000e+00 0.000000e+00 [208,] 1.000000e+00 0.000000e+00 0.000000e+00 [209,] 1.000000e+00 0.000000e+00 0.000000e+00 [210,] 1.000000e+00 0.000000e+00 0.000000e+00 [211,] 1.000000e+00 0.000000e+00 0.000000e+00 [212,] 1.000000e+00 0.000000e+00 0.000000e+00 [213,] 1.000000e+00 0.000000e+00 0.000000e+00 [214,] 1.000000e+00 0.000000e+00 0.000000e+00 [215,] 1.000000e+00 0.000000e+00 0.000000e+00 [216,] 1.000000e+00 0.000000e+00 0.000000e+00 [217,] 1.000000e+00 0.000000e+00 0.000000e+00 [218,] 1.000000e+00 0.000000e+00 0.000000e+00 [219,] 1.000000e+00 0.000000e+00 0.000000e+00 [220,] 1.000000e+00 0.000000e+00 0.000000e+00 [221,] 1.000000e+00 0.000000e+00 0.000000e+00 [222,] 1.000000e+00 0.000000e+00 0.000000e+00 [223,] 1.000000e+00 0.000000e+00 0.000000e+00 [224,] 1.000000e+00 0.000000e+00 0.000000e+00 [225,] 1.000000e+00 0.000000e+00 0.000000e+00 [226,] 1.000000e+00 0.000000e+00 0.000000e+00 [227,] 1.000000e+00 0.000000e+00 0.000000e+00 [228,] 1.000000e+00 0.000000e+00 0.000000e+00 [229,] 1.000000e+00 0.000000e+00 0.000000e+00 [230,] 1.000000e+00 0.000000e+00 0.000000e+00 [231,] 1.000000e+00 0.000000e+00 0.000000e+00 [232,] 1.000000e+00 0.000000e+00 0.000000e+00 [233,] 1.000000e+00 0.000000e+00 0.000000e+00 [234,] 1.000000e+00 0.000000e+00 0.000000e+00 [235,] 1.000000e+00 0.000000e+00 0.000000e+00 [236,] 1.000000e+00 0.000000e+00 0.000000e+00 [237,] 1.000000e+00 0.000000e+00 0.000000e+00 [238,] 1.000000e+00 0.000000e+00 0.000000e+00 [239,] 1.000000e+00 0.000000e+00 0.000000e+00 [240,] 1.000000e+00 0.000000e+00 0.000000e+00 [241,] 1.000000e+00 0.000000e+00 0.000000e+00 [242,] 1.000000e+00 0.000000e+00 0.000000e+00 [243,] 1.000000e+00 0.000000e+00 0.000000e+00 [244,] 1.000000e+00 0.000000e+00 0.000000e+00 [245,] 1.000000e+00 0.000000e+00 0.000000e+00 [246,] 1.000000e+00 0.000000e+00 0.000000e+00 [247,] 1.000000e+00 0.000000e+00 0.000000e+00 [248,] 1.000000e+00 0.000000e+00 0.000000e+00 [249,] 1.000000e+00 0.000000e+00 0.000000e+00 [250,] 1.000000e+00 0.000000e+00 0.000000e+00 [251,] 1.000000e+00 0.000000e+00 0.000000e+00 [252,] 1.000000e+00 0.000000e+00 0.000000e+00 [253,] 1.000000e+00 0.000000e+00 0.000000e+00 [254,] 1.000000e+00 0.000000e+00 0.000000e+00 [255,] 1.000000e+00 0.000000e+00 0.000000e+00 [256,] 1.000000e+00 0.000000e+00 0.000000e+00 [257,] 1.000000e+00 0.000000e+00 0.000000e+00 [258,] 1.000000e+00 0.000000e+00 0.000000e+00 [259,] 1.000000e+00 0.000000e+00 0.000000e+00 [260,] 1.000000e+00 0.000000e+00 0.000000e+00 [261,] 1.000000e+00 0.000000e+00 0.000000e+00 [262,] 1.000000e+00 0.000000e+00 0.000000e+00 [263,] 1.000000e+00 0.000000e+00 0.000000e+00 [264,] 1.000000e+00 0.000000e+00 0.000000e+00 [265,] 1.000000e+00 0.000000e+00 0.000000e+00 [266,] 1.000000e+00 0.000000e+00 0.000000e+00 [267,] 1.000000e+00 0.000000e+00 0.000000e+00 [268,] 1.000000e+00 0.000000e+00 0.000000e+00 [269,] 1.000000e+00 0.000000e+00 0.000000e+00 [270,] 1.000000e+00 0.000000e+00 0.000000e+00 [271,] 1.000000e+00 0.000000e+00 0.000000e+00 [272,] 1.000000e+00 0.000000e+00 0.000000e+00 [273,] 1.000000e+00 0.000000e+00 0.000000e+00 [274,] 1.000000e+00 0.000000e+00 0.000000e+00 [275,] 1.000000e+00 0.000000e+00 0.000000e+00 [276,] 1.000000e+00 0.000000e+00 0.000000e+00 [277,] 1.000000e+00 0.000000e+00 0.000000e+00 [278,] 1.000000e+00 0.000000e+00 0.000000e+00 [279,] 1.000000e+00 0.000000e+00 0.000000e+00 [280,] 1.000000e+00 0.000000e+00 0.000000e+00 [281,] 1.000000e+00 0.000000e+00 0.000000e+00 [282,] 1.000000e+00 0.000000e+00 0.000000e+00 [283,] 1.000000e+00 0.000000e+00 0.000000e+00 [284,] 1.000000e+00 0.000000e+00 0.000000e+00 [285,] 1.000000e+00 0.000000e+00 0.000000e+00 [286,] 1.000000e+00 0.000000e+00 0.000000e+00 [287,] 1.000000e+00 0.000000e+00 0.000000e+00 [288,] 1.000000e+00 0.000000e+00 0.000000e+00 [289,] 1.000000e+00 0.000000e+00 0.000000e+00 [290,] 1.000000e+00 0.000000e+00 0.000000e+00 [291,] 1.000000e+00 0.000000e+00 0.000000e+00 [292,] 1.000000e+00 0.000000e+00 0.000000e+00 [293,] 1.000000e+00 0.000000e+00 0.000000e+00 [294,] 1.000000e+00 0.000000e+00 0.000000e+00 [295,] 1.000000e+00 0.000000e+00 0.000000e+00 [296,] 1.000000e+00 0.000000e+00 0.000000e+00 [297,] 1.000000e+00 0.000000e+00 0.000000e+00 [298,] 1.000000e+00 0.000000e+00 0.000000e+00 [299,] 1.000000e+00 0.000000e+00 0.000000e+00 [300,] 1.000000e+00 0.000000e+00 0.000000e+00 [301,] 1.000000e+00 0.000000e+00 0.000000e+00 [302,] 1.000000e+00 0.000000e+00 0.000000e+00 [303,] 1.000000e+00 0.000000e+00 0.000000e+00 [304,] 1.000000e+00 0.000000e+00 0.000000e+00 [305,] 1.000000e+00 0.000000e+00 0.000000e+00 [306,] 1.000000e+00 0.000000e+00 0.000000e+00 [307,] 1.000000e+00 0.000000e+00 0.000000e+00 [308,] 1.000000e+00 0.000000e+00 0.000000e+00 [309,] 1.000000e+00 0.000000e+00 0.000000e+00 [310,] 1.000000e+00 9.881313e-324 4.940656e-324 [311,] 1.000000e+00 1.077063e-320 5.385316e-321 [312,] 1.000000e+00 1.531626e-317 7.658131e-318 [313,] 1.000000e+00 2.291963e-314 1.145981e-314 [314,] 1.000000e+00 3.403938e-311 1.701969e-311 [315,] 1.000000e+00 5.208096e-308 2.604048e-308 [316,] 1.000000e+00 7.605405e-305 3.802703e-305 [317,] 1.000000e+00 1.036104e-301 5.180519e-302 [318,] 1.000000e+00 1.494995e-298 7.474976e-299 [319,] 1.000000e+00 6.476543e-296 3.238271e-296 [320,] 1.000000e+00 4.149540e-293 2.074770e-293 [321,] 1.000000e+00 2.511571e-290 1.255786e-290 [322,] 1.000000e+00 3.692146e-287 1.846073e-287 [323,] 1.000000e+00 5.378287e-284 2.689143e-284 [324,] 1.000000e+00 2.445436e-281 1.222718e-281 [325,] 1.000000e+00 1.054342e-278 5.271711e-279 [326,] 1.000000e+00 2.622121e-276 1.311060e-276 [327,] 1.000000e+00 4.068665e-273 2.034333e-273 [328,] 1.000000e+00 6.048531e-270 3.024265e-270 [329,] 1.000000e+00 9.656588e-267 4.828294e-267 [330,] 1.000000e+00 1.458358e-263 7.291788e-264 [331,] 1.000000e+00 6.126811e-261 3.063406e-261 [332,] 1.000000e+00 9.155059e-258 4.577530e-258 [333,] 1.000000e+00 1.261405e-254 6.307025e-255 [334,] 1.000000e+00 1.854940e-251 9.274701e-252 [335,] 1.000000e+00 2.721525e-248 1.360763e-248 [336,] 1.000000e+00 4.233826e-245 2.116913e-245 [337,] 1.000000e+00 1.816616e-242 9.083079e-243 [338,] 1.000000e+00 2.710570e-239 1.355285e-239 [339,] 1.000000e+00 3.972915e-236 1.986457e-236 [340,] 1.000000e+00 4.288895e-233 2.144448e-233 [341,] 1.000000e+00 6.207579e-230 3.103789e-230 [342,] 1.000000e+00 8.989598e-227 4.494799e-227 [343,] 1.000000e+00 1.278874e-223 6.394368e-224 [344,] 1.000000e+00 5.768274e-221 2.884137e-221 [345,] 1.000000e+00 8.321561e-218 4.160781e-218 [346,] 1.000000e+00 1.188932e-214 5.944661e-215 [347,] 1.000000e+00 5.773763e-212 2.886881e-212 [348,] 1.000000e+00 8.328352e-209 4.164176e-209 [349,] 1.000000e+00 3.568128e-206 1.784064e-206 [350,] 1.000000e+00 5.227315e-203 2.613657e-203 [351,] 1.000000e+00 2.129162e-200 1.064581e-200 [352,] 1.000000e+00 3.741876e-198 1.870938e-198 [353,] 1.000000e+00 4.974826e-195 2.487413e-195 [354,] 1.000000e+00 6.819293e-192 3.409647e-192 [355,] 1.000000e+00 2.788510e-189 1.394255e-189 [356,] 1.000000e+00 4.557584e-186 2.278792e-186 [357,] 1.000000e+00 7.057526e-183 3.528763e-183 [358,] 1.000000e+00 1.185826e-179 5.929129e-180 [359,] 1.000000e+00 1.911288e-176 9.556440e-177 [360,] 1.000000e+00 2.864863e-173 1.432432e-173 [361,] 1.000000e+00 1.079954e-170 5.399771e-171 [362,] 1.000000e+00 1.293567e-167 6.467834e-168 [363,] 1.000000e+00 2.026949e-164 1.013474e-164 [364,] 1.000000e+00 2.965199e-161 1.482599e-161 [365,] 1.000000e+00 4.538975e-158 2.269488e-158 [366,] 1.000000e+00 2.529751e-155 1.264875e-155 [367,] 1.000000e+00 3.833139e-152 1.916570e-152 [368,] 1.000000e+00 8.938501e-150 4.469250e-150 [369,] 1.000000e+00 2.452597e-147 1.226299e-147 [370,] 1.000000e+00 3.957988e-144 1.978994e-144 [371,] 1.000000e+00 6.325725e-141 3.162863e-141 [372,] 1.000000e+00 9.909567e-138 4.954784e-138 [373,] 1.000000e+00 1.410438e-134 7.052189e-135 [374,] 1.000000e+00 4.460902e-132 2.230451e-132 [375,] 1.000000e+00 7.006714e-129 3.503357e-129 [376,] 1.000000e+00 1.118503e-125 5.592516e-126 [377,] 1.000000e+00 1.747088e-122 8.735442e-123 [378,] 1.000000e+00 6.685272e-120 3.342636e-120 [379,] 1.000000e+00 2.138947e-117 1.069474e-117 [380,] 1.000000e+00 3.571440e-114 1.785720e-114 [381,] 1.000000e+00 8.163125e-112 4.081563e-112 [382,] 1.000000e+00 1.433081e-108 7.165406e-109 [383,] 1.000000e+00 2.474733e-105 1.237367e-105 [384,] 1.000000e+00 4.201036e-102 2.100518e-102 [385,] 1.000000e+00 7.982663e-99 3.991332e-99 [386,] 1.000000e+00 1.316488e-95 6.582439e-96 [387,] 1.000000e+00 3.080963e-93 1.540482e-93 [388,] 1.000000e+00 5.390582e-90 2.695291e-90 [389,] 1.000000e+00 9.248397e-87 4.624199e-87 [390,] 1.000000e+00 1.554471e-83 7.772355e-84 [391,] 1.000000e+00 2.573945e-80 1.286972e-80 [392,] 1.000000e+00 4.630680e-77 2.315340e-77 [393,] 1.000000e+00 7.223479e-75 3.611739e-75 [394,] 1.000000e+00 7.106914e-72 3.553457e-72 [395,] 1.000000e+00 1.148832e-68 5.744161e-69 [396,] 1.000000e+00 1.021952e-65 5.109762e-66 [397,] 1.000000e+00 1.780568e-62 8.902840e-63 [398,] 1.000000e+00 2.697331e-59 1.348666e-59 [399,] 1.000000e+00 4.678282e-57 2.339141e-57 [400,] 1.000000e+00 9.482406e-54 4.741203e-54 [401,] 1.000000e+00 6.771965e-51 3.385983e-51 [402,] 1.000000e+00 1.394498e-47 6.972492e-48 [403,] 1.000000e+00 2.547762e-44 1.273881e-44 [404,] 1.000000e+00 1.309126e-42 6.545628e-43 [405,] 1.000000e+00 3.321918e-39 1.660959e-39 [406,] 1.000000e+00 8.204901e-36 4.102451e-36 [407,] 1.000000e+00 1.161418e-32 5.807089e-33 [408,] 1.000000e+00 9.278476e-33 4.639238e-33 [409,] 1.000000e+00 4.304074e-29 2.152037e-29 [410,] 1.000000e+00 3.303227e-27 1.651613e-27 [411,] 1.000000e+00 3.025130e-23 1.512565e-23 [412,] 1.000000e+00 2.580880e-19 1.290440e-19 [413,] 1.000000e+00 1.989474e-15 9.947370e-16 [414,] 1.000000e+00 3.302574e-12 1.651287e-12 > postscript(file="/var/www/wessaorg/rcomp/tmp/1n3yc1293639134.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/wessaorg/rcomp/tmp/2n3yc1293639134.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/wessaorg/rcomp/tmp/3n3yc1293639134.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/wessaorg/rcomp/tmp/4gcff1293639134.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/wessaorg/rcomp/tmp/5gcff1293639134.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 -392.1615896 -1019.7563837 -453.1708724 -804.9260504 360.7594783 6 7 8 9 10 -279.0235781 -321.9235182 -463.1453074 -588.8512942 -401.8339358 11 12 13 14 15 -742.3992506 -764.9346651 -470.5344475 -322.5549404 -429.2427177 16 17 18 19 20 -312.1408758 -468.3176302 -484.7509013 -481.9832552 -460.0219261 21 22 23 24 25 -470.5344475 -358.3457368 -470.5344475 -1036.0095528 -705.2365682 26 27 28 29 30 -450.7663177 -356.2347827 -550.7845271 -470.5344475 -534.4957791 31 32 33 34 35 134.8952139 -470.5344475 -538.0085226 -1460.5526455 -221.6013291 36 37 38 39 40 -24.2183041 -468.0152199 -137.0108679 -730.3905299 -660.3075539 41 42 43 44 45 -1040.4291487 -773.7211062 -670.8199921 -470.5344475 -355.1175642 46 47 48 49 50 -351.4452806 -470.5344475 -427.9140673 -495.4713487 -790.6830912 51 52 53 54 55 -294.7326741 -631.9214648 -165.3021303 -614.0887094 -3205.6775652 56 57 58 59 60 -1013.6751239 -450.9765612 0.2260205 -314.2046138 -470.5344475 61 62 63 64 65 -470.5344475 -314.6810109 104.9268068 -472.3085284 730.8280805 66 67 68 69 70 -465.0419398 -403.4219727 -515.6566081 -373.6832878 -470.5344475 71 72 73 74 75 -3378.6463520 297.7009682 8771.2754479 -470.5344475 -519.9980963 76 77 78 79 80 -530.6796029 -580.5956424 -458.4692428 1065.4420485 -180.9583579 81 82 83 84 85 -467.6732617 -470.5344475 -470.5344475 -502.2014647 -470.5344475 86 87 88 89 90 -470.5449476 -460.3507269 -470.2362739 -374.2722757 -569.5353756 91 92 93 94 95 -470.5449476 -359.2706537 -934.1458012 545.0145513 -505.2476610 96 97 98 99 100 -546.9441833 -470.5344475 -451.2447467 -477.8018249 -250.8679047 101 102 103 104 105 -1031.2981959 -377.9418429 -476.3206260 -470.3494670 -470.5344475 106 107 108 109 110 1977.8810030 -390.3697034 -1731.1824846 -336.8265263 -458.2571058 111 112 113 114 115 -476.6329577 160.5315987 -470.5344475 -446.6151630 -470.5449476 116 117 118 119 120 -416.7551092 -760.3086408 -1896.0208159 -2010.2224842 -1339.9715274 121 122 123 124 125 87.0682504 -470.5344475 -2069.8586465 -739.1645505 19.0167286 126 127 128 129 130 -470.5449476 -449.8260321 -470.5344475 1713.7209097 -386.0915604 131 132 133 134 135 -512.7086149 -470.5449476 -470.5449476 -470.5449476 -459.9319090 136 137 138 139 140 -409.3731914 315.3994283 271.9793364 315.3994283 315.3994283 141 142 143 144 145 -100.8336478 822.4580612 315.3994283 314.3994283 324.6291894 146 147 148 149 150 315.3994283 325.5178793 397.9112373 247.2483103 314.3994283 151 152 153 154 155 315.3994283 345.5406820 315.3994283 355.1577545 258.5359260 156 157 158 159 160 162.0830619 304.1319653 315.3994283 297.6329451 -327.9527531 161 162 163 164 165 322.2723763 315.3994283 376.0361236 269.1539227 314.3994283 166 167 168 169 170 -236.8998305 315.3994283 320.9979119 661.8501481 216.4049847 171 172 173 174 175 -620.2062952 -196.2867470 125.9323451 -73.3090202 493.2488721 176 177 178 179 180 35.4909132 315.3994283 317.0125444 307.5333072 315.3994283 181 182 183 184 185 -486.2091934 -275.9641179 315.3994283 314.3994283 315.3994283 186 187 188 189 190 315.7741848 304.8064760 273.8043439 2420.8634789 315.3994283 191 192 193 194 195 315.0847766 280.6204554 315.3994283 315.8441910 172.6253245 196 197 198 199 200 315.3994283 -288.3499046 546.4077924 315.3994283 314.3994283 201 202 203 204 205 242.3681110 314.3994283 315.3994283 321.7409222 315.3994283 206 207 208 209 210 -58.4255405 36.9547142 315.3994283 -42.4623428 411.5226277 211 212 213 214 215 -227.2172816 155.9266370 304.7732975 352.6581697 -1025.0363920 216 217 218 219 220 139.5015577 302.6968945 314.3994283 210.8843809 -469.1418232 221 222 223 224 225 314.0518467 315.3994283 328.6259273 291.0572476 351.7487357 226 227 228 229 230 -58.8305535 315.3994283 -2765.8363590 159.2645634 262.0055034 231 232 233 234 235 -3005.3518062 381.8741556 -210.0222391 315.3994283 317.0127834 236 237 238 239 240 190.4365370 389.9566574 315.3994283 315.3994283 296.4621526 241 242 243 244 245 137.2401429 64.8448672 315.3994283 322.0098201 486.7721881 246 247 248 249 250 692.0842751 315.3994283 330.7688591 315.3994283 315.0736409 251 252 253 254 255 315.3994283 74.4691667 315.3994283 314.3994283 315.3994283 256 257 258 259 260 315.3994283 315.3994283 315.3994283 7.6359586 -31.0356061 261 262 263 264 265 456.7097840 621.1805299 217.9851292 -453.9265658 194.3743963 266 267 268 269 270 315.3994283 316.3843162 315.3994283 179.3722050 313.5008966 271 272 273 274 275 315.3994283 315.3994283 208.2126011 312.8224200 325.6942428 276 277 278 279 280 455.9650879 315.7686817 335.0404225 357.0034963 325.4853920 281 282 283 284 285 321.8809852 324.1832544 315.3994283 277.4683358 315.7987361 286 287 288 289 290 334.8665928 315.3994283 216.7893830 740.5007750 1520.2646473 291 292 293 294 295 315.3994283 323.0512807 315.3994283 278.0031205 314.1572864 296 297 298 299 300 315.1544314 314.3994283 1238.8720897 314.3994283 314.3994283 301 302 303 304 305 650.5540791 247.7804980 -300.3015591 308.0942909 255.5358441 306 307 308 309 310 325.8886133 -393.2884841 315.3994283 315.3994283 314.6297502 311 312 313 314 315 315.3994283 50.1942145 273.0411048 313.4492315 217.8344257 316 317 318 319 320 97.0458816 -555.3755167 351.7294826 315.3994283 209.5370554 321 322 323 324 325 314.4961618 336.5024572 -10.4298297 315.3994283 -418.8891389 326 327 328 329 330 315.3994283 -750.7199719 -1012.2136382 1177.9550887 41.7899478 331 332 333 334 335 153.5326488 314.3994283 314.3994283 -1482.4629570 120.4144274 336 337 338 339 340 315.3994283 -110.3046243 -266.9372740 314.3994283 315.3994283 341 342 343 344 345 -255.2058402 307.5539870 315.3994283 -63.5229823 314.3994283 346 347 348 349 350 62.0263233 306.0083792 522.2048761 315.3994283 325.7692883 351 352 353 354 355 315.3994283 314.3994283 315.1602883 315.3994283 180.4040504 356 357 358 359 360 315.3994283 314.3994283 315.3994283 314.3994283 -673.0071334 361 362 363 364 365 -931.6266442 615.9307969 382.3329805 -45.3781510 342.8701291 366 367 368 369 370 539.3398741 329.3313646 321.3694143 314.3994283 300.6128820 371 372 373 374 375 -404.8978638 355.7087606 288.3572654 116.7488226 290.8510489 376 377 378 379 380 -200.6059805 -48.7009041 54.3181962 269.8938791 315.3994283 381 382 383 384 385 242.9919123 250.9922927 111.3867762 283.4126263 313.5231695 386 387 388 389 390 351.3505507 348.2096572 315.3994283 314.3994283 315.3994283 391 392 393 394 395 315.3994283 315.3994283 95.3249937 328.4920881 314.3994283 396 397 398 399 400 315.3994283 315.3994283 315.3994283 330.2700782 -262.0387481 401 402 403 404 405 -298.5011399 44.4788159 404.2790165 -345.2528458 292.0649056 406 407 408 409 410 325.1858653 317.9770930 -56.1874205 -180.8788554 -466.4287385 411 412 413 414 415 296.6902579 314.3994283 315.3994283 206.2554596 323.9391740 416 417 418 419 420 331.7373343 -1380.2628987 60.2731184 315.4376576 315.3994283 421 422 423 424 425 321.5951448 649.9036455 -1784.7849797 323.8040500 218.8220148 426 427 428 429 430 494.7405549 315.3994283 148.9079456 216.2296820 315.3994283 431 377.5161438 > postscript(file="/var/www/wessaorg/rcomp/tmp/6gcff1293639134.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 -392.1615896 NA 1 -1019.7563837 -392.1615896 2 -453.1708724 -1019.7563837 3 -804.9260504 -453.1708724 4 360.7594783 -804.9260504 5 -279.0235781 360.7594783 6 -321.9235182 -279.0235781 7 -463.1453074 -321.9235182 8 -588.8512942 -463.1453074 9 -401.8339358 -588.8512942 10 -742.3992506 -401.8339358 11 -764.9346651 -742.3992506 12 -470.5344475 -764.9346651 13 -322.5549404 -470.5344475 14 -429.2427177 -322.5549404 15 -312.1408758 -429.2427177 16 -468.3176302 -312.1408758 17 -484.7509013 -468.3176302 18 -481.9832552 -484.7509013 19 -460.0219261 -481.9832552 20 -470.5344475 -460.0219261 21 -358.3457368 -470.5344475 22 -470.5344475 -358.3457368 23 -1036.0095528 -470.5344475 24 -705.2365682 -1036.0095528 25 -450.7663177 -705.2365682 26 -356.2347827 -450.7663177 27 -550.7845271 -356.2347827 28 -470.5344475 -550.7845271 29 -534.4957791 -470.5344475 30 134.8952139 -534.4957791 31 -470.5344475 134.8952139 32 -538.0085226 -470.5344475 33 -1460.5526455 -538.0085226 34 -221.6013291 -1460.5526455 35 -24.2183041 -221.6013291 36 -468.0152199 -24.2183041 37 -137.0108679 -468.0152199 38 -730.3905299 -137.0108679 39 -660.3075539 -730.3905299 40 -1040.4291487 -660.3075539 41 -773.7211062 -1040.4291487 42 -670.8199921 -773.7211062 43 -470.5344475 -670.8199921 44 -355.1175642 -470.5344475 45 -351.4452806 -355.1175642 46 -470.5344475 -351.4452806 47 -427.9140673 -470.5344475 48 -495.4713487 -427.9140673 49 -790.6830912 -495.4713487 50 -294.7326741 -790.6830912 51 -631.9214648 -294.7326741 52 -165.3021303 -631.9214648 53 -614.0887094 -165.3021303 54 -3205.6775652 -614.0887094 55 -1013.6751239 -3205.6775652 56 -450.9765612 -1013.6751239 57 0.2260205 -450.9765612 58 -314.2046138 0.2260205 59 -470.5344475 -314.2046138 60 -470.5344475 -470.5344475 61 -314.6810109 -470.5344475 62 104.9268068 -314.6810109 63 -472.3085284 104.9268068 64 730.8280805 -472.3085284 65 -465.0419398 730.8280805 66 -403.4219727 -465.0419398 67 -515.6566081 -403.4219727 68 -373.6832878 -515.6566081 69 -470.5344475 -373.6832878 70 -3378.6463520 -470.5344475 71 297.7009682 -3378.6463520 72 8771.2754479 297.7009682 73 -470.5344475 8771.2754479 74 -519.9980963 -470.5344475 75 -530.6796029 -519.9980963 76 -580.5956424 -530.6796029 77 -458.4692428 -580.5956424 78 1065.4420485 -458.4692428 79 -180.9583579 1065.4420485 80 -467.6732617 -180.9583579 81 -470.5344475 -467.6732617 82 -470.5344475 -470.5344475 83 -502.2014647 -470.5344475 84 -470.5344475 -502.2014647 85 -470.5449476 -470.5344475 86 -460.3507269 -470.5449476 87 -470.2362739 -460.3507269 88 -374.2722757 -470.2362739 89 -569.5353756 -374.2722757 90 -470.5449476 -569.5353756 91 -359.2706537 -470.5449476 92 -934.1458012 -359.2706537 93 545.0145513 -934.1458012 94 -505.2476610 545.0145513 95 -546.9441833 -505.2476610 96 -470.5344475 -546.9441833 97 -451.2447467 -470.5344475 98 -477.8018249 -451.2447467 99 -250.8679047 -477.8018249 100 -1031.2981959 -250.8679047 101 -377.9418429 -1031.2981959 102 -476.3206260 -377.9418429 103 -470.3494670 -476.3206260 104 -470.5344475 -470.3494670 105 1977.8810030 -470.5344475 106 -390.3697034 1977.8810030 107 -1731.1824846 -390.3697034 108 -336.8265263 -1731.1824846 109 -458.2571058 -336.8265263 110 -476.6329577 -458.2571058 111 160.5315987 -476.6329577 112 -470.5344475 160.5315987 113 -446.6151630 -470.5344475 114 -470.5449476 -446.6151630 115 -416.7551092 -470.5449476 116 -760.3086408 -416.7551092 117 -1896.0208159 -760.3086408 118 -2010.2224842 -1896.0208159 119 -1339.9715274 -2010.2224842 120 87.0682504 -1339.9715274 121 -470.5344475 87.0682504 122 -2069.8586465 -470.5344475 123 -739.1645505 -2069.8586465 124 19.0167286 -739.1645505 125 -470.5449476 19.0167286 126 -449.8260321 -470.5449476 127 -470.5344475 -449.8260321 128 1713.7209097 -470.5344475 129 -386.0915604 1713.7209097 130 -512.7086149 -386.0915604 131 -470.5449476 -512.7086149 132 -470.5449476 -470.5449476 133 -470.5449476 -470.5449476 134 -459.9319090 -470.5449476 135 -409.3731914 -459.9319090 136 315.3994283 -409.3731914 137 271.9793364 315.3994283 138 315.3994283 271.9793364 139 315.3994283 315.3994283 140 -100.8336478 315.3994283 141 822.4580612 -100.8336478 142 315.3994283 822.4580612 143 314.3994283 315.3994283 144 324.6291894 314.3994283 145 315.3994283 324.6291894 146 325.5178793 315.3994283 147 397.9112373 325.5178793 148 247.2483103 397.9112373 149 314.3994283 247.2483103 150 315.3994283 314.3994283 151 345.5406820 315.3994283 152 315.3994283 345.5406820 153 355.1577545 315.3994283 154 258.5359260 355.1577545 155 162.0830619 258.5359260 156 304.1319653 162.0830619 157 315.3994283 304.1319653 158 297.6329451 315.3994283 159 -327.9527531 297.6329451 160 322.2723763 -327.9527531 161 315.3994283 322.2723763 162 376.0361236 315.3994283 163 269.1539227 376.0361236 164 314.3994283 269.1539227 165 -236.8998305 314.3994283 166 315.3994283 -236.8998305 167 320.9979119 315.3994283 168 661.8501481 320.9979119 169 216.4049847 661.8501481 170 -620.2062952 216.4049847 171 -196.2867470 -620.2062952 172 125.9323451 -196.2867470 173 -73.3090202 125.9323451 174 493.2488721 -73.3090202 175 35.4909132 493.2488721 176 315.3994283 35.4909132 177 317.0125444 315.3994283 178 307.5333072 317.0125444 179 315.3994283 307.5333072 180 -486.2091934 315.3994283 181 -275.9641179 -486.2091934 182 315.3994283 -275.9641179 183 314.3994283 315.3994283 184 315.3994283 314.3994283 185 315.7741848 315.3994283 186 304.8064760 315.7741848 187 273.8043439 304.8064760 188 2420.8634789 273.8043439 189 315.3994283 2420.8634789 190 315.0847766 315.3994283 191 280.6204554 315.0847766 192 315.3994283 280.6204554 193 315.8441910 315.3994283 194 172.6253245 315.8441910 195 315.3994283 172.6253245 196 -288.3499046 315.3994283 197 546.4077924 -288.3499046 198 315.3994283 546.4077924 199 314.3994283 315.3994283 200 242.3681110 314.3994283 201 314.3994283 242.3681110 202 315.3994283 314.3994283 203 321.7409222 315.3994283 204 315.3994283 321.7409222 205 -58.4255405 315.3994283 206 36.9547142 -58.4255405 207 315.3994283 36.9547142 208 -42.4623428 315.3994283 209 411.5226277 -42.4623428 210 -227.2172816 411.5226277 211 155.9266370 -227.2172816 212 304.7732975 155.9266370 213 352.6581697 304.7732975 214 -1025.0363920 352.6581697 215 139.5015577 -1025.0363920 216 302.6968945 139.5015577 217 314.3994283 302.6968945 218 210.8843809 314.3994283 219 -469.1418232 210.8843809 220 314.0518467 -469.1418232 221 315.3994283 314.0518467 222 328.6259273 315.3994283 223 291.0572476 328.6259273 224 351.7487357 291.0572476 225 -58.8305535 351.7487357 226 315.3994283 -58.8305535 227 -2765.8363590 315.3994283 228 159.2645634 -2765.8363590 229 262.0055034 159.2645634 230 -3005.3518062 262.0055034 231 381.8741556 -3005.3518062 232 -210.0222391 381.8741556 233 315.3994283 -210.0222391 234 317.0127834 315.3994283 235 190.4365370 317.0127834 236 389.9566574 190.4365370 237 315.3994283 389.9566574 238 315.3994283 315.3994283 239 296.4621526 315.3994283 240 137.2401429 296.4621526 241 64.8448672 137.2401429 242 315.3994283 64.8448672 243 322.0098201 315.3994283 244 486.7721881 322.0098201 245 692.0842751 486.7721881 246 315.3994283 692.0842751 247 330.7688591 315.3994283 248 315.3994283 330.7688591 249 315.0736409 315.3994283 250 315.3994283 315.0736409 251 74.4691667 315.3994283 252 315.3994283 74.4691667 253 314.3994283 315.3994283 254 315.3994283 314.3994283 255 315.3994283 315.3994283 256 315.3994283 315.3994283 257 315.3994283 315.3994283 258 7.6359586 315.3994283 259 -31.0356061 7.6359586 260 456.7097840 -31.0356061 261 621.1805299 456.7097840 262 217.9851292 621.1805299 263 -453.9265658 217.9851292 264 194.3743963 -453.9265658 265 315.3994283 194.3743963 266 316.3843162 315.3994283 267 315.3994283 316.3843162 268 179.3722050 315.3994283 269 313.5008966 179.3722050 270 315.3994283 313.5008966 271 315.3994283 315.3994283 272 208.2126011 315.3994283 273 312.8224200 208.2126011 274 325.6942428 312.8224200 275 455.9650879 325.6942428 276 315.7686817 455.9650879 277 335.0404225 315.7686817 278 357.0034963 335.0404225 279 325.4853920 357.0034963 280 321.8809852 325.4853920 281 324.1832544 321.8809852 282 315.3994283 324.1832544 283 277.4683358 315.3994283 284 315.7987361 277.4683358 285 334.8665928 315.7987361 286 315.3994283 334.8665928 287 216.7893830 315.3994283 288 740.5007750 216.7893830 289 1520.2646473 740.5007750 290 315.3994283 1520.2646473 291 323.0512807 315.3994283 292 315.3994283 323.0512807 293 278.0031205 315.3994283 294 314.1572864 278.0031205 295 315.1544314 314.1572864 296 314.3994283 315.1544314 297 1238.8720897 314.3994283 298 314.3994283 1238.8720897 299 314.3994283 314.3994283 300 650.5540791 314.3994283 301 247.7804980 650.5540791 302 -300.3015591 247.7804980 303 308.0942909 -300.3015591 304 255.5358441 308.0942909 305 325.8886133 255.5358441 306 -393.2884841 325.8886133 307 315.3994283 -393.2884841 308 315.3994283 315.3994283 309 314.6297502 315.3994283 310 315.3994283 314.6297502 311 50.1942145 315.3994283 312 273.0411048 50.1942145 313 313.4492315 273.0411048 314 217.8344257 313.4492315 315 97.0458816 217.8344257 316 -555.3755167 97.0458816 317 351.7294826 -555.3755167 318 315.3994283 351.7294826 319 209.5370554 315.3994283 320 314.4961618 209.5370554 321 336.5024572 314.4961618 322 -10.4298297 336.5024572 323 315.3994283 -10.4298297 324 -418.8891389 315.3994283 325 315.3994283 -418.8891389 326 -750.7199719 315.3994283 327 -1012.2136382 -750.7199719 328 1177.9550887 -1012.2136382 329 41.7899478 1177.9550887 330 153.5326488 41.7899478 331 314.3994283 153.5326488 332 314.3994283 314.3994283 333 -1482.4629570 314.3994283 334 120.4144274 -1482.4629570 335 315.3994283 120.4144274 336 -110.3046243 315.3994283 337 -266.9372740 -110.3046243 338 314.3994283 -266.9372740 339 315.3994283 314.3994283 340 -255.2058402 315.3994283 341 307.5539870 -255.2058402 342 315.3994283 307.5539870 343 -63.5229823 315.3994283 344 314.3994283 -63.5229823 345 62.0263233 314.3994283 346 306.0083792 62.0263233 347 522.2048761 306.0083792 348 315.3994283 522.2048761 349 325.7692883 315.3994283 350 315.3994283 325.7692883 351 314.3994283 315.3994283 352 315.1602883 314.3994283 353 315.3994283 315.1602883 354 180.4040504 315.3994283 355 315.3994283 180.4040504 356 314.3994283 315.3994283 357 315.3994283 314.3994283 358 314.3994283 315.3994283 359 -673.0071334 314.3994283 360 -931.6266442 -673.0071334 361 615.9307969 -931.6266442 362 382.3329805 615.9307969 363 -45.3781510 382.3329805 364 342.8701291 -45.3781510 365 539.3398741 342.8701291 366 329.3313646 539.3398741 367 321.3694143 329.3313646 368 314.3994283 321.3694143 369 300.6128820 314.3994283 370 -404.8978638 300.6128820 371 355.7087606 -404.8978638 372 288.3572654 355.7087606 373 116.7488226 288.3572654 374 290.8510489 116.7488226 375 -200.6059805 290.8510489 376 -48.7009041 -200.6059805 377 54.3181962 -48.7009041 378 269.8938791 54.3181962 379 315.3994283 269.8938791 380 242.9919123 315.3994283 381 250.9922927 242.9919123 382 111.3867762 250.9922927 383 283.4126263 111.3867762 384 313.5231695 283.4126263 385 351.3505507 313.5231695 386 348.2096572 351.3505507 387 315.3994283 348.2096572 388 314.3994283 315.3994283 389 315.3994283 314.3994283 390 315.3994283 315.3994283 391 315.3994283 315.3994283 392 95.3249937 315.3994283 393 328.4920881 95.3249937 394 314.3994283 328.4920881 395 315.3994283 314.3994283 396 315.3994283 315.3994283 397 315.3994283 315.3994283 398 330.2700782 315.3994283 399 -262.0387481 330.2700782 400 -298.5011399 -262.0387481 401 44.4788159 -298.5011399 402 404.2790165 44.4788159 403 -345.2528458 404.2790165 404 292.0649056 -345.2528458 405 325.1858653 292.0649056 406 317.9770930 325.1858653 407 -56.1874205 317.9770930 408 -180.8788554 -56.1874205 409 -466.4287385 -180.8788554 410 296.6902579 -466.4287385 411 314.3994283 296.6902579 412 315.3994283 314.3994283 413 206.2554596 315.3994283 414 323.9391740 206.2554596 415 331.7373343 323.9391740 416 -1380.2628987 331.7373343 417 60.2731184 -1380.2628987 418 315.4376576 60.2731184 419 315.3994283 315.4376576 420 321.5951448 315.3994283 421 649.9036455 321.5951448 422 -1784.7849797 649.9036455 423 323.8040500 -1784.7849797 424 218.8220148 323.8040500 425 494.7405549 218.8220148 426 315.3994283 494.7405549 427 148.9079456 315.3994283 428 216.2296820 148.9079456 429 315.3994283 216.2296820 430 377.5161438 315.3994283 431 NA 377.5161438 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -1019.7563837 -392.1615896 [2,] -453.1708724 -1019.7563837 [3,] -804.9260504 -453.1708724 [4,] 360.7594783 -804.9260504 [5,] -279.0235781 360.7594783 [6,] -321.9235182 -279.0235781 [7,] -463.1453074 -321.9235182 [8,] -588.8512942 -463.1453074 [9,] -401.8339358 -588.8512942 [10,] -742.3992506 -401.8339358 [11,] -764.9346651 -742.3992506 [12,] -470.5344475 -764.9346651 [13,] -322.5549404 -470.5344475 [14,] -429.2427177 -322.5549404 [15,] -312.1408758 -429.2427177 [16,] -468.3176302 -312.1408758 [17,] -484.7509013 -468.3176302 [18,] -481.9832552 -484.7509013 [19,] -460.0219261 -481.9832552 [20,] -470.5344475 -460.0219261 [21,] -358.3457368 -470.5344475 [22,] -470.5344475 -358.3457368 [23,] -1036.0095528 -470.5344475 [24,] -705.2365682 -1036.0095528 [25,] -450.7663177 -705.2365682 [26,] -356.2347827 -450.7663177 [27,] -550.7845271 -356.2347827 [28,] -470.5344475 -550.7845271 [29,] -534.4957791 -470.5344475 [30,] 134.8952139 -534.4957791 [31,] -470.5344475 134.8952139 [32,] -538.0085226 -470.5344475 [33,] -1460.5526455 -538.0085226 [34,] -221.6013291 -1460.5526455 [35,] -24.2183041 -221.6013291 [36,] -468.0152199 -24.2183041 [37,] -137.0108679 -468.0152199 [38,] -730.3905299 -137.0108679 [39,] -660.3075539 -730.3905299 [40,] -1040.4291487 -660.3075539 [41,] -773.7211062 -1040.4291487 [42,] -670.8199921 -773.7211062 [43,] -470.5344475 -670.8199921 [44,] -355.1175642 -470.5344475 [45,] -351.4452806 -355.1175642 [46,] -470.5344475 -351.4452806 [47,] -427.9140673 -470.5344475 [48,] -495.4713487 -427.9140673 [49,] -790.6830912 -495.4713487 [50,] -294.7326741 -790.6830912 [51,] -631.9214648 -294.7326741 [52,] -165.3021303 -631.9214648 [53,] -614.0887094 -165.3021303 [54,] -3205.6775652 -614.0887094 [55,] -1013.6751239 -3205.6775652 [56,] -450.9765612 -1013.6751239 [57,] 0.2260205 -450.9765612 [58,] -314.2046138 0.2260205 [59,] -470.5344475 -314.2046138 [60,] -470.5344475 -470.5344475 [61,] -314.6810109 -470.5344475 [62,] 104.9268068 -314.6810109 [63,] -472.3085284 104.9268068 [64,] 730.8280805 -472.3085284 [65,] -465.0419398 730.8280805 [66,] -403.4219727 -465.0419398 [67,] -515.6566081 -403.4219727 [68,] -373.6832878 -515.6566081 [69,] -470.5344475 -373.6832878 [70,] -3378.6463520 -470.5344475 [71,] 297.7009682 -3378.6463520 [72,] 8771.2754479 297.7009682 [73,] -470.5344475 8771.2754479 [74,] -519.9980963 -470.5344475 [75,] -530.6796029 -519.9980963 [76,] -580.5956424 -530.6796029 [77,] -458.4692428 -580.5956424 [78,] 1065.4420485 -458.4692428 [79,] -180.9583579 1065.4420485 [80,] -467.6732617 -180.9583579 [81,] -470.5344475 -467.6732617 [82,] -470.5344475 -470.5344475 [83,] -502.2014647 -470.5344475 [84,] -470.5344475 -502.2014647 [85,] -470.5449476 -470.5344475 [86,] -460.3507269 -470.5449476 [87,] -470.2362739 -460.3507269 [88,] -374.2722757 -470.2362739 [89,] -569.5353756 -374.2722757 [90,] -470.5449476 -569.5353756 [91,] -359.2706537 -470.5449476 [92,] -934.1458012 -359.2706537 [93,] 545.0145513 -934.1458012 [94,] -505.2476610 545.0145513 [95,] -546.9441833 -505.2476610 [96,] -470.5344475 -546.9441833 [97,] -451.2447467 -470.5344475 [98,] -477.8018249 -451.2447467 [99,] -250.8679047 -477.8018249 [100,] -1031.2981959 -250.8679047 [101,] -377.9418429 -1031.2981959 [102,] -476.3206260 -377.9418429 [103,] -470.3494670 -476.3206260 [104,] -470.5344475 -470.3494670 [105,] 1977.8810030 -470.5344475 [106,] -390.3697034 1977.8810030 [107,] -1731.1824846 -390.3697034 [108,] -336.8265263 -1731.1824846 [109,] -458.2571058 -336.8265263 [110,] -476.6329577 -458.2571058 [111,] 160.5315987 -476.6329577 [112,] -470.5344475 160.5315987 [113,] -446.6151630 -470.5344475 [114,] -470.5449476 -446.6151630 [115,] -416.7551092 -470.5449476 [116,] -760.3086408 -416.7551092 [117,] -1896.0208159 -760.3086408 [118,] -2010.2224842 -1896.0208159 [119,] -1339.9715274 -2010.2224842 [120,] 87.0682504 -1339.9715274 [121,] -470.5344475 87.0682504 [122,] -2069.8586465 -470.5344475 [123,] -739.1645505 -2069.8586465 [124,] 19.0167286 -739.1645505 [125,] -470.5449476 19.0167286 [126,] -449.8260321 -470.5449476 [127,] -470.5344475 -449.8260321 [128,] 1713.7209097 -470.5344475 [129,] -386.0915604 1713.7209097 [130,] -512.7086149 -386.0915604 [131,] -470.5449476 -512.7086149 [132,] -470.5449476 -470.5449476 [133,] -470.5449476 -470.5449476 [134,] -459.9319090 -470.5449476 [135,] -409.3731914 -459.9319090 [136,] 315.3994283 -409.3731914 [137,] 271.9793364 315.3994283 [138,] 315.3994283 271.9793364 [139,] 315.3994283 315.3994283 [140,] -100.8336478 315.3994283 [141,] 822.4580612 -100.8336478 [142,] 315.3994283 822.4580612 [143,] 314.3994283 315.3994283 [144,] 324.6291894 314.3994283 [145,] 315.3994283 324.6291894 [146,] 325.5178793 315.3994283 [147,] 397.9112373 325.5178793 [148,] 247.2483103 397.9112373 [149,] 314.3994283 247.2483103 [150,] 315.3994283 314.3994283 [151,] 345.5406820 315.3994283 [152,] 315.3994283 345.5406820 [153,] 355.1577545 315.3994283 [154,] 258.5359260 355.1577545 [155,] 162.0830619 258.5359260 [156,] 304.1319653 162.0830619 [157,] 315.3994283 304.1319653 [158,] 297.6329451 315.3994283 [159,] -327.9527531 297.6329451 [160,] 322.2723763 -327.9527531 [161,] 315.3994283 322.2723763 [162,] 376.0361236 315.3994283 [163,] 269.1539227 376.0361236 [164,] 314.3994283 269.1539227 [165,] -236.8998305 314.3994283 [166,] 315.3994283 -236.8998305 [167,] 320.9979119 315.3994283 [168,] 661.8501481 320.9979119 [169,] 216.4049847 661.8501481 [170,] -620.2062952 216.4049847 [171,] -196.2867470 -620.2062952 [172,] 125.9323451 -196.2867470 [173,] -73.3090202 125.9323451 [174,] 493.2488721 -73.3090202 [175,] 35.4909132 493.2488721 [176,] 315.3994283 35.4909132 [177,] 317.0125444 315.3994283 [178,] 307.5333072 317.0125444 [179,] 315.3994283 307.5333072 [180,] -486.2091934 315.3994283 [181,] -275.9641179 -486.2091934 [182,] 315.3994283 -275.9641179 [183,] 314.3994283 315.3994283 [184,] 315.3994283 314.3994283 [185,] 315.7741848 315.3994283 [186,] 304.8064760 315.7741848 [187,] 273.8043439 304.8064760 [188,] 2420.8634789 273.8043439 [189,] 315.3994283 2420.8634789 [190,] 315.0847766 315.3994283 [191,] 280.6204554 315.0847766 [192,] 315.3994283 280.6204554 [193,] 315.8441910 315.3994283 [194,] 172.6253245 315.8441910 [195,] 315.3994283 172.6253245 [196,] -288.3499046 315.3994283 [197,] 546.4077924 -288.3499046 [198,] 315.3994283 546.4077924 [199,] 314.3994283 315.3994283 [200,] 242.3681110 314.3994283 [201,] 314.3994283 242.3681110 [202,] 315.3994283 314.3994283 [203,] 321.7409222 315.3994283 [204,] 315.3994283 321.7409222 [205,] -58.4255405 315.3994283 [206,] 36.9547142 -58.4255405 [207,] 315.3994283 36.9547142 [208,] -42.4623428 315.3994283 [209,] 411.5226277 -42.4623428 [210,] -227.2172816 411.5226277 [211,] 155.9266370 -227.2172816 [212,] 304.7732975 155.9266370 [213,] 352.6581697 304.7732975 [214,] -1025.0363920 352.6581697 [215,] 139.5015577 -1025.0363920 [216,] 302.6968945 139.5015577 [217,] 314.3994283 302.6968945 [218,] 210.8843809 314.3994283 [219,] -469.1418232 210.8843809 [220,] 314.0518467 -469.1418232 [221,] 315.3994283 314.0518467 [222,] 328.6259273 315.3994283 [223,] 291.0572476 328.6259273 [224,] 351.7487357 291.0572476 [225,] -58.8305535 351.7487357 [226,] 315.3994283 -58.8305535 [227,] -2765.8363590 315.3994283 [228,] 159.2645634 -2765.8363590 [229,] 262.0055034 159.2645634 [230,] -3005.3518062 262.0055034 [231,] 381.8741556 -3005.3518062 [232,] -210.0222391 381.8741556 [233,] 315.3994283 -210.0222391 [234,] 317.0127834 315.3994283 [235,] 190.4365370 317.0127834 [236,] 389.9566574 190.4365370 [237,] 315.3994283 389.9566574 [238,] 315.3994283 315.3994283 [239,] 296.4621526 315.3994283 [240,] 137.2401429 296.4621526 [241,] 64.8448672 137.2401429 [242,] 315.3994283 64.8448672 [243,] 322.0098201 315.3994283 [244,] 486.7721881 322.0098201 [245,] 692.0842751 486.7721881 [246,] 315.3994283 692.0842751 [247,] 330.7688591 315.3994283 [248,] 315.3994283 330.7688591 [249,] 315.0736409 315.3994283 [250,] 315.3994283 315.0736409 [251,] 74.4691667 315.3994283 [252,] 315.3994283 74.4691667 [253,] 314.3994283 315.3994283 [254,] 315.3994283 314.3994283 [255,] 315.3994283 315.3994283 [256,] 315.3994283 315.3994283 [257,] 315.3994283 315.3994283 [258,] 7.6359586 315.3994283 [259,] -31.0356061 7.6359586 [260,] 456.7097840 -31.0356061 [261,] 621.1805299 456.7097840 [262,] 217.9851292 621.1805299 [263,] -453.9265658 217.9851292 [264,] 194.3743963 -453.9265658 [265,] 315.3994283 194.3743963 [266,] 316.3843162 315.3994283 [267,] 315.3994283 316.3843162 [268,] 179.3722050 315.3994283 [269,] 313.5008966 179.3722050 [270,] 315.3994283 313.5008966 [271,] 315.3994283 315.3994283 [272,] 208.2126011 315.3994283 [273,] 312.8224200 208.2126011 [274,] 325.6942428 312.8224200 [275,] 455.9650879 325.6942428 [276,] 315.7686817 455.9650879 [277,] 335.0404225 315.7686817 [278,] 357.0034963 335.0404225 [279,] 325.4853920 357.0034963 [280,] 321.8809852 325.4853920 [281,] 324.1832544 321.8809852 [282,] 315.3994283 324.1832544 [283,] 277.4683358 315.3994283 [284,] 315.7987361 277.4683358 [285,] 334.8665928 315.7987361 [286,] 315.3994283 334.8665928 [287,] 216.7893830 315.3994283 [288,] 740.5007750 216.7893830 [289,] 1520.2646473 740.5007750 [290,] 315.3994283 1520.2646473 [291,] 323.0512807 315.3994283 [292,] 315.3994283 323.0512807 [293,] 278.0031205 315.3994283 [294,] 314.1572864 278.0031205 [295,] 315.1544314 314.1572864 [296,] 314.3994283 315.1544314 [297,] 1238.8720897 314.3994283 [298,] 314.3994283 1238.8720897 [299,] 314.3994283 314.3994283 [300,] 650.5540791 314.3994283 [301,] 247.7804980 650.5540791 [302,] -300.3015591 247.7804980 [303,] 308.0942909 -300.3015591 [304,] 255.5358441 308.0942909 [305,] 325.8886133 255.5358441 [306,] -393.2884841 325.8886133 [307,] 315.3994283 -393.2884841 [308,] 315.3994283 315.3994283 [309,] 314.6297502 315.3994283 [310,] 315.3994283 314.6297502 [311,] 50.1942145 315.3994283 [312,] 273.0411048 50.1942145 [313,] 313.4492315 273.0411048 [314,] 217.8344257 313.4492315 [315,] 97.0458816 217.8344257 [316,] -555.3755167 97.0458816 [317,] 351.7294826 -555.3755167 [318,] 315.3994283 351.7294826 [319,] 209.5370554 315.3994283 [320,] 314.4961618 209.5370554 [321,] 336.5024572 314.4961618 [322,] -10.4298297 336.5024572 [323,] 315.3994283 -10.4298297 [324,] -418.8891389 315.3994283 [325,] 315.3994283 -418.8891389 [326,] -750.7199719 315.3994283 [327,] -1012.2136382 -750.7199719 [328,] 1177.9550887 -1012.2136382 [329,] 41.7899478 1177.9550887 [330,] 153.5326488 41.7899478 [331,] 314.3994283 153.5326488 [332,] 314.3994283 314.3994283 [333,] -1482.4629570 314.3994283 [334,] 120.4144274 -1482.4629570 [335,] 315.3994283 120.4144274 [336,] -110.3046243 315.3994283 [337,] -266.9372740 -110.3046243 [338,] 314.3994283 -266.9372740 [339,] 315.3994283 314.3994283 [340,] -255.2058402 315.3994283 [341,] 307.5539870 -255.2058402 [342,] 315.3994283 307.5539870 [343,] -63.5229823 315.3994283 [344,] 314.3994283 -63.5229823 [345,] 62.0263233 314.3994283 [346,] 306.0083792 62.0263233 [347,] 522.2048761 306.0083792 [348,] 315.3994283 522.2048761 [349,] 325.7692883 315.3994283 [350,] 315.3994283 325.7692883 [351,] 314.3994283 315.3994283 [352,] 315.1602883 314.3994283 [353,] 315.3994283 315.1602883 [354,] 180.4040504 315.3994283 [355,] 315.3994283 180.4040504 [356,] 314.3994283 315.3994283 [357,] 315.3994283 314.3994283 [358,] 314.3994283 315.3994283 [359,] -673.0071334 314.3994283 [360,] -931.6266442 -673.0071334 [361,] 615.9307969 -931.6266442 [362,] 382.3329805 615.9307969 [363,] -45.3781510 382.3329805 [364,] 342.8701291 -45.3781510 [365,] 539.3398741 342.8701291 [366,] 329.3313646 539.3398741 [367,] 321.3694143 329.3313646 [368,] 314.3994283 321.3694143 [369,] 300.6128820 314.3994283 [370,] -404.8978638 300.6128820 [371,] 355.7087606 -404.8978638 [372,] 288.3572654 355.7087606 [373,] 116.7488226 288.3572654 [374,] 290.8510489 116.7488226 [375,] -200.6059805 290.8510489 [376,] -48.7009041 -200.6059805 [377,] 54.3181962 -48.7009041 [378,] 269.8938791 54.3181962 [379,] 315.3994283 269.8938791 [380,] 242.9919123 315.3994283 [381,] 250.9922927 242.9919123 [382,] 111.3867762 250.9922927 [383,] 283.4126263 111.3867762 [384,] 313.5231695 283.4126263 [385,] 351.3505507 313.5231695 [386,] 348.2096572 351.3505507 [387,] 315.3994283 348.2096572 [388,] 314.3994283 315.3994283 [389,] 315.3994283 314.3994283 [390,] 315.3994283 315.3994283 [391,] 315.3994283 315.3994283 [392,] 95.3249937 315.3994283 [393,] 328.4920881 95.3249937 [394,] 314.3994283 328.4920881 [395,] 315.3994283 314.3994283 [396,] 315.3994283 315.3994283 [397,] 315.3994283 315.3994283 [398,] 330.2700782 315.3994283 [399,] -262.0387481 330.2700782 [400,] -298.5011399 -262.0387481 [401,] 44.4788159 -298.5011399 [402,] 404.2790165 44.4788159 [403,] -345.2528458 404.2790165 [404,] 292.0649056 -345.2528458 [405,] 325.1858653 292.0649056 [406,] 317.9770930 325.1858653 [407,] -56.1874205 317.9770930 [408,] -180.8788554 -56.1874205 [409,] -466.4287385 -180.8788554 [410,] 296.6902579 -466.4287385 [411,] 314.3994283 296.6902579 [412,] 315.3994283 314.3994283 [413,] 206.2554596 315.3994283 [414,] 323.9391740 206.2554596 [415,] 331.7373343 323.9391740 [416,] -1380.2628987 331.7373343 [417,] 60.2731184 -1380.2628987 [418,] 315.4376576 60.2731184 [419,] 315.3994283 315.4376576 [420,] 321.5951448 315.3994283 [421,] 649.9036455 321.5951448 [422,] -1784.7849797 649.9036455 [423,] 323.8040500 -1784.7849797 [424,] 218.8220148 323.8040500 [425,] 494.7405549 218.8220148 [426,] 315.3994283 494.7405549 [427,] 148.9079456 315.3994283 [428,] 216.2296820 148.9079456 [429,] 315.3994283 216.2296820 [430,] 377.5161438 315.3994283 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -1019.7563837 -392.1615896 2 -453.1708724 -1019.7563837 3 -804.9260504 -453.1708724 4 360.7594783 -804.9260504 5 -279.0235781 360.7594783 6 -321.9235182 -279.0235781 7 -463.1453074 -321.9235182 8 -588.8512942 -463.1453074 9 -401.8339358 -588.8512942 10 -742.3992506 -401.8339358 11 -764.9346651 -742.3992506 12 -470.5344475 -764.9346651 13 -322.5549404 -470.5344475 14 -429.2427177 -322.5549404 15 -312.1408758 -429.2427177 16 -468.3176302 -312.1408758 17 -484.7509013 -468.3176302 18 -481.9832552 -484.7509013 19 -460.0219261 -481.9832552 20 -470.5344475 -460.0219261 21 -358.3457368 -470.5344475 22 -470.5344475 -358.3457368 23 -1036.0095528 -470.5344475 24 -705.2365682 -1036.0095528 25 -450.7663177 -705.2365682 26 -356.2347827 -450.7663177 27 -550.7845271 -356.2347827 28 -470.5344475 -550.7845271 29 -534.4957791 -470.5344475 30 134.8952139 -534.4957791 31 -470.5344475 134.8952139 32 -538.0085226 -470.5344475 33 -1460.5526455 -538.0085226 34 -221.6013291 -1460.5526455 35 -24.2183041 -221.6013291 36 -468.0152199 -24.2183041 37 -137.0108679 -468.0152199 38 -730.3905299 -137.0108679 39 -660.3075539 -730.3905299 40 -1040.4291487 -660.3075539 41 -773.7211062 -1040.4291487 42 -670.8199921 -773.7211062 43 -470.5344475 -670.8199921 44 -355.1175642 -470.5344475 45 -351.4452806 -355.1175642 46 -470.5344475 -351.4452806 47 -427.9140673 -470.5344475 48 -495.4713487 -427.9140673 49 -790.6830912 -495.4713487 50 -294.7326741 -790.6830912 51 -631.9214648 -294.7326741 52 -165.3021303 -631.9214648 53 -614.0887094 -165.3021303 54 -3205.6775652 -614.0887094 55 -1013.6751239 -3205.6775652 56 -450.9765612 -1013.6751239 57 0.2260205 -450.9765612 58 -314.2046138 0.2260205 59 -470.5344475 -314.2046138 60 -470.5344475 -470.5344475 61 -314.6810109 -470.5344475 62 104.9268068 -314.6810109 63 -472.3085284 104.9268068 64 730.8280805 -472.3085284 65 -465.0419398 730.8280805 66 -403.4219727 -465.0419398 67 -515.6566081 -403.4219727 68 -373.6832878 -515.6566081 69 -470.5344475 -373.6832878 70 -3378.6463520 -470.5344475 71 297.7009682 -3378.6463520 72 8771.2754479 297.7009682 73 -470.5344475 8771.2754479 74 -519.9980963 -470.5344475 75 -530.6796029 -519.9980963 76 -580.5956424 -530.6796029 77 -458.4692428 -580.5956424 78 1065.4420485 -458.4692428 79 -180.9583579 1065.4420485 80 -467.6732617 -180.9583579 81 -470.5344475 -467.6732617 82 -470.5344475 -470.5344475 83 -502.2014647 -470.5344475 84 -470.5344475 -502.2014647 85 -470.5449476 -470.5344475 86 -460.3507269 -470.5449476 87 -470.2362739 -460.3507269 88 -374.2722757 -470.2362739 89 -569.5353756 -374.2722757 90 -470.5449476 -569.5353756 91 -359.2706537 -470.5449476 92 -934.1458012 -359.2706537 93 545.0145513 -934.1458012 94 -505.2476610 545.0145513 95 -546.9441833 -505.2476610 96 -470.5344475 -546.9441833 97 -451.2447467 -470.5344475 98 -477.8018249 -451.2447467 99 -250.8679047 -477.8018249 100 -1031.2981959 -250.8679047 101 -377.9418429 -1031.2981959 102 -476.3206260 -377.9418429 103 -470.3494670 -476.3206260 104 -470.5344475 -470.3494670 105 1977.8810030 -470.5344475 106 -390.3697034 1977.8810030 107 -1731.1824846 -390.3697034 108 -336.8265263 -1731.1824846 109 -458.2571058 -336.8265263 110 -476.6329577 -458.2571058 111 160.5315987 -476.6329577 112 -470.5344475 160.5315987 113 -446.6151630 -470.5344475 114 -470.5449476 -446.6151630 115 -416.7551092 -470.5449476 116 -760.3086408 -416.7551092 117 -1896.0208159 -760.3086408 118 -2010.2224842 -1896.0208159 119 -1339.9715274 -2010.2224842 120 87.0682504 -1339.9715274 121 -470.5344475 87.0682504 122 -2069.8586465 -470.5344475 123 -739.1645505 -2069.8586465 124 19.0167286 -739.1645505 125 -470.5449476 19.0167286 126 -449.8260321 -470.5449476 127 -470.5344475 -449.8260321 128 1713.7209097 -470.5344475 129 -386.0915604 1713.7209097 130 -512.7086149 -386.0915604 131 -470.5449476 -512.7086149 132 -470.5449476 -470.5449476 133 -470.5449476 -470.5449476 134 -459.9319090 -470.5449476 135 -409.3731914 -459.9319090 136 315.3994283 -409.3731914 137 271.9793364 315.3994283 138 315.3994283 271.9793364 139 315.3994283 315.3994283 140 -100.8336478 315.3994283 141 822.4580612 -100.8336478 142 315.3994283 822.4580612 143 314.3994283 315.3994283 144 324.6291894 314.3994283 145 315.3994283 324.6291894 146 325.5178793 315.3994283 147 397.9112373 325.5178793 148 247.2483103 397.9112373 149 314.3994283 247.2483103 150 315.3994283 314.3994283 151 345.5406820 315.3994283 152 315.3994283 345.5406820 153 355.1577545 315.3994283 154 258.5359260 355.1577545 155 162.0830619 258.5359260 156 304.1319653 162.0830619 157 315.3994283 304.1319653 158 297.6329451 315.3994283 159 -327.9527531 297.6329451 160 322.2723763 -327.9527531 161 315.3994283 322.2723763 162 376.0361236 315.3994283 163 269.1539227 376.0361236 164 314.3994283 269.1539227 165 -236.8998305 314.3994283 166 315.3994283 -236.8998305 167 320.9979119 315.3994283 168 661.8501481 320.9979119 169 216.4049847 661.8501481 170 -620.2062952 216.4049847 171 -196.2867470 -620.2062952 172 125.9323451 -196.2867470 173 -73.3090202 125.9323451 174 493.2488721 -73.3090202 175 35.4909132 493.2488721 176 315.3994283 35.4909132 177 317.0125444 315.3994283 178 307.5333072 317.0125444 179 315.3994283 307.5333072 180 -486.2091934 315.3994283 181 -275.9641179 -486.2091934 182 315.3994283 -275.9641179 183 314.3994283 315.3994283 184 315.3994283 314.3994283 185 315.7741848 315.3994283 186 304.8064760 315.7741848 187 273.8043439 304.8064760 188 2420.8634789 273.8043439 189 315.3994283 2420.8634789 190 315.0847766 315.3994283 191 280.6204554 315.0847766 192 315.3994283 280.6204554 193 315.8441910 315.3994283 194 172.6253245 315.8441910 195 315.3994283 172.6253245 196 -288.3499046 315.3994283 197 546.4077924 -288.3499046 198 315.3994283 546.4077924 199 314.3994283 315.3994283 200 242.3681110 314.3994283 201 314.3994283 242.3681110 202 315.3994283 314.3994283 203 321.7409222 315.3994283 204 315.3994283 321.7409222 205 -58.4255405 315.3994283 206 36.9547142 -58.4255405 207 315.3994283 36.9547142 208 -42.4623428 315.3994283 209 411.5226277 -42.4623428 210 -227.2172816 411.5226277 211 155.9266370 -227.2172816 212 304.7732975 155.9266370 213 352.6581697 304.7732975 214 -1025.0363920 352.6581697 215 139.5015577 -1025.0363920 216 302.6968945 139.5015577 217 314.3994283 302.6968945 218 210.8843809 314.3994283 219 -469.1418232 210.8843809 220 314.0518467 -469.1418232 221 315.3994283 314.0518467 222 328.6259273 315.3994283 223 291.0572476 328.6259273 224 351.7487357 291.0572476 225 -58.8305535 351.7487357 226 315.3994283 -58.8305535 227 -2765.8363590 315.3994283 228 159.2645634 -2765.8363590 229 262.0055034 159.2645634 230 -3005.3518062 262.0055034 231 381.8741556 -3005.3518062 232 -210.0222391 381.8741556 233 315.3994283 -210.0222391 234 317.0127834 315.3994283 235 190.4365370 317.0127834 236 389.9566574 190.4365370 237 315.3994283 389.9566574 238 315.3994283 315.3994283 239 296.4621526 315.3994283 240 137.2401429 296.4621526 241 64.8448672 137.2401429 242 315.3994283 64.8448672 243 322.0098201 315.3994283 244 486.7721881 322.0098201 245 692.0842751 486.7721881 246 315.3994283 692.0842751 247 330.7688591 315.3994283 248 315.3994283 330.7688591 249 315.0736409 315.3994283 250 315.3994283 315.0736409 251 74.4691667 315.3994283 252 315.3994283 74.4691667 253 314.3994283 315.3994283 254 315.3994283 314.3994283 255 315.3994283 315.3994283 256 315.3994283 315.3994283 257 315.3994283 315.3994283 258 7.6359586 315.3994283 259 -31.0356061 7.6359586 260 456.7097840 -31.0356061 261 621.1805299 456.7097840 262 217.9851292 621.1805299 263 -453.9265658 217.9851292 264 194.3743963 -453.9265658 265 315.3994283 194.3743963 266 316.3843162 315.3994283 267 315.3994283 316.3843162 268 179.3722050 315.3994283 269 313.5008966 179.3722050 270 315.3994283 313.5008966 271 315.3994283 315.3994283 272 208.2126011 315.3994283 273 312.8224200 208.2126011 274 325.6942428 312.8224200 275 455.9650879 325.6942428 276 315.7686817 455.9650879 277 335.0404225 315.7686817 278 357.0034963 335.0404225 279 325.4853920 357.0034963 280 321.8809852 325.4853920 281 324.1832544 321.8809852 282 315.3994283 324.1832544 283 277.4683358 315.3994283 284 315.7987361 277.4683358 285 334.8665928 315.7987361 286 315.3994283 334.8665928 287 216.7893830 315.3994283 288 740.5007750 216.7893830 289 1520.2646473 740.5007750 290 315.3994283 1520.2646473 291 323.0512807 315.3994283 292 315.3994283 323.0512807 293 278.0031205 315.3994283 294 314.1572864 278.0031205 295 315.1544314 314.1572864 296 314.3994283 315.1544314 297 1238.8720897 314.3994283 298 314.3994283 1238.8720897 299 314.3994283 314.3994283 300 650.5540791 314.3994283 301 247.7804980 650.5540791 302 -300.3015591 247.7804980 303 308.0942909 -300.3015591 304 255.5358441 308.0942909 305 325.8886133 255.5358441 306 -393.2884841 325.8886133 307 315.3994283 -393.2884841 308 315.3994283 315.3994283 309 314.6297502 315.3994283 310 315.3994283 314.6297502 311 50.1942145 315.3994283 312 273.0411048 50.1942145 313 313.4492315 273.0411048 314 217.8344257 313.4492315 315 97.0458816 217.8344257 316 -555.3755167 97.0458816 317 351.7294826 -555.3755167 318 315.3994283 351.7294826 319 209.5370554 315.3994283 320 314.4961618 209.5370554 321 336.5024572 314.4961618 322 -10.4298297 336.5024572 323 315.3994283 -10.4298297 324 -418.8891389 315.3994283 325 315.3994283 -418.8891389 326 -750.7199719 315.3994283 327 -1012.2136382 -750.7199719 328 1177.9550887 -1012.2136382 329 41.7899478 1177.9550887 330 153.5326488 41.7899478 331 314.3994283 153.5326488 332 314.3994283 314.3994283 333 -1482.4629570 314.3994283 334 120.4144274 -1482.4629570 335 315.3994283 120.4144274 336 -110.3046243 315.3994283 337 -266.9372740 -110.3046243 338 314.3994283 -266.9372740 339 315.3994283 314.3994283 340 -255.2058402 315.3994283 341 307.5539870 -255.2058402 342 315.3994283 307.5539870 343 -63.5229823 315.3994283 344 314.3994283 -63.5229823 345 62.0263233 314.3994283 346 306.0083792 62.0263233 347 522.2048761 306.0083792 348 315.3994283 522.2048761 349 325.7692883 315.3994283 350 315.3994283 325.7692883 351 314.3994283 315.3994283 352 315.1602883 314.3994283 353 315.3994283 315.1602883 354 180.4040504 315.3994283 355 315.3994283 180.4040504 356 314.3994283 315.3994283 357 315.3994283 314.3994283 358 314.3994283 315.3994283 359 -673.0071334 314.3994283 360 -931.6266442 -673.0071334 361 615.9307969 -931.6266442 362 382.3329805 615.9307969 363 -45.3781510 382.3329805 364 342.8701291 -45.3781510 365 539.3398741 342.8701291 366 329.3313646 539.3398741 367 321.3694143 329.3313646 368 314.3994283 321.3694143 369 300.6128820 314.3994283 370 -404.8978638 300.6128820 371 355.7087606 -404.8978638 372 288.3572654 355.7087606 373 116.7488226 288.3572654 374 290.8510489 116.7488226 375 -200.6059805 290.8510489 376 -48.7009041 -200.6059805 377 54.3181962 -48.7009041 378 269.8938791 54.3181962 379 315.3994283 269.8938791 380 242.9919123 315.3994283 381 250.9922927 242.9919123 382 111.3867762 250.9922927 383 283.4126263 111.3867762 384 313.5231695 283.4126263 385 351.3505507 313.5231695 386 348.2096572 351.3505507 387 315.3994283 348.2096572 388 314.3994283 315.3994283 389 315.3994283 314.3994283 390 315.3994283 315.3994283 391 315.3994283 315.3994283 392 95.3249937 315.3994283 393 328.4920881 95.3249937 394 314.3994283 328.4920881 395 315.3994283 314.3994283 396 315.3994283 315.3994283 397 315.3994283 315.3994283 398 330.2700782 315.3994283 399 -262.0387481 330.2700782 400 -298.5011399 -262.0387481 401 44.4788159 -298.5011399 402 404.2790165 44.4788159 403 -345.2528458 404.2790165 404 292.0649056 -345.2528458 405 325.1858653 292.0649056 406 317.9770930 325.1858653 407 -56.1874205 317.9770930 408 -180.8788554 -56.1874205 409 -466.4287385 -180.8788554 410 296.6902579 -466.4287385 411 314.3994283 296.6902579 412 315.3994283 314.3994283 413 206.2554596 315.3994283 414 323.9391740 206.2554596 415 331.7373343 323.9391740 416 -1380.2628987 331.7373343 417 60.2731184 -1380.2628987 418 315.4376576 60.2731184 419 315.3994283 315.4376576 420 321.5951448 315.3994283 421 649.9036455 321.5951448 422 -1784.7849797 649.9036455 423 323.8040500 -1784.7849797 424 218.8220148 323.8040500 425 494.7405549 218.8220148 426 315.3994283 494.7405549 427 148.9079456 315.3994283 428 216.2296820 148.9079456 429 315.3994283 216.2296820 430 377.5161438 315.3994283 > 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/wessaorg/rcomp/tmp/7q3wi1293639134.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/wessaorg/rcomp/tmp/81ve31293639134.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/wessaorg/rcomp/tmp/91ve31293639134.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/wessaorg/rcomp/tmp/10c4dn1293639134.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/wessaorg/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/wessaorg/rcomp/tmp/11f4bb1293639134.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/wessaorg/rcomp/tmp/128ebw1293639134.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/wessaorg/rcomp/tmp/13fx8q1293639134.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/wessaorg/rcomp/tmp/14popt1293639134.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/wessaorg/rcomp/tmp/154yn21293639134.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/wessaorg/rcomp/tmp/167gl81293639134.tab") + } > > try(system("convert tmp/1n3yc1293639134.ps tmp/1n3yc1293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/2n3yc1293639134.ps tmp/2n3yc1293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/3n3yc1293639134.ps tmp/3n3yc1293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/4gcff1293639134.ps tmp/4gcff1293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/5gcff1293639134.ps tmp/5gcff1293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/6gcff1293639134.ps tmp/6gcff1293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/7q3wi1293639134.ps tmp/7q3wi1293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/81ve31293639134.ps tmp/81ve31293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/91ve31293639134.ps tmp/91ve31293639134.png",intern=TRUE)) character(0) > try(system("convert tmp/10c4dn1293639134.ps tmp/10c4dn1293639134.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 16.48 0.50 17.27