R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1 + ,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 = '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 Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x CCScore\r Group Costs Trades Dividends Wealth t 1 0 1 1595 17 60720 319369 1 2 39 0 5565 47 94355 493408 2 3 0 0 601 6 60720 319210 3 4 0 1 188 11 77655 381180 4 5 305 1 7146 105 134028 297978 5 6 -9 0 1135 17 62285 290476 6 7 -1 1 450 8 59325 292136 7 8 0 1 34 8 60630 314353 8 9 0 1 133 14 65990 339445 9 10 0 1 119 6 59118 303677 10 11 33 0 2053 53 100423 397144 11 12 32 0 4036 63 100269 424898 12 13 0 1 0 0 60720 315380 13 14 0 1 4655 393 60720 341570 14 15 0 1 131 11 61808 308989 15 16 0 1 1766 73 60438 305959 16 17 0 1 312 14 58598 318690 17 18 0 1 448 40 59781 323361 18 19 0 1 115 13 60945 318903 19 20 0 1 60 10 61124 314049 20 21 0 1 0 0 60720 315380 21 22 0 1 364 35 47705 298466 22 23 0 1 0 0 60720 315380 23 24 -3 0 1442 118 121173 438493 24 25 6 0 1389 9 57530 378049 25 26 0 1 149 5 62041 313332 26 27 0 1 2212 14 60720 319864 27 28 70 0 7489 322 136996 430866 28 29 0 1 0 0 60720 315380 29 30 5 1 402 26 84990 332743 30 31 -5 0 7419 29 63255 286963 31 32 0 1 0 0 60720 315380 32 33 0 1 307 15 58856 331955 33 34 58 0 1134 162 146216 527021 34 35 37 0 7561 87 82425 364304 35 36 9 1 5131 71 86111 292154 36 37 0 1 9 6 60835 314987 37 38 -31 1 2264 24 55174 272713 38 39 2865 1 40949 481 129352 1398893 39 40 0 0 4980 91 113521 409642 40 41 14 1 272 8 112995 429112 41 42 0 1 757 19 45689 382712 42 43 32 1 1424 94 131116 374943 43 44 0 1 0 0 60720 315380 44 45 0 0 378 52 71873 297413 45 46 -4 1 1061 5 70415 304252 46 47 0 1 0 0 60720 315380 47 48 0 1 203 9 61938 309596 48 49 74 0 3689 25 101792 377305 49 50 20 1 567 11 103345 386212 50 51 247 0 27664 282 74996 657954 51 52 11 1 1686 40 68696 368186 52 53 -6 0 1164 17 57320 269753 53 54 284 1 15824 64 63838 585715 54 55 819 1 6575 159 151352 1071292 55 56 652 0 37597 248 37884 992426 56 57 0 1 314 5 60720 315371 57 58 10 0 1932 36 72954 249898 58 59 22 0 5692 72 64239 357312 59 60 0 1 0 0 60720 315380 60 61 0 1 0 0 60720 315380 61 62 16 0 6533 75 126910 364839 62 63 -3 0 2126 138 21001 230621 63 64 0 1 14 1 61016 315877 64 65 10 0 3192 60 59831 125390 65 66 0 1 47 3 60720 314882 66 67 73 0 4891 217 210568 370837 67 68 0 0 531 10 71641 330068 68 69 0 1 2332 4 72680 324385 69 70 0 1 0 0 60720 315380 70 71 1974 0 16275 109 47818 1443586 71 72 -129 1 9252 48 79510 253537 72 73 28328 0 29790 309 81767 4321023 73 74 0 1 0 0 60720 315380 74 75 0 0 2300 29 64873 352108 75 76 0 1 267 32 57640 330059 76 77 0 0 382 21 59365 340968 77 78 0 1 397 28 60646 317736 78 79 238 1 11920 182 31080 209458 79 80 3040 0 54660 302 126942 1491348 80 81 0 1 5 2 60720 314887 81 82 0 1 0 0 60720 315380 82 83 0 1 0 0 60720 315380 83 84 1 0 2647 37 62920 353058 84 85 0 1 0 0 60720 315380 85 86 0 0 0 0 60720 315380 86 87 0 1 94 7 60793 314533 87 88 0 1 3 2 60698 315354 88 89 0 1 454 51 63261 302187 89 90 0 1 227 9 76644 336639 90 91 0 0 0 0 60720 315380 91 92 0 1 206 9 53110 296702 92 93 1910 1 25830 115 245546 1073089 93 94 38 1 1528 38 60326 146494 94 95 0 1 271 2 60720 325249 95 96 0 1 138 6 69817 331420 96 97 0 1 0 0 60720 315380 97 98 0 1 278 26 61404 314922 98 99 0 1 282 30 62452 320016 99 100 0 1 571 80 41477 280398 100 101 17 1 2253 99 63593 452469 101 102 0 1 290 2 58790 301164 102 103 0 1 78 16 62700 317330 103 104 0 1 20 2 60805 315576 104 105 0 1 0 0 60720 315380 105 106 -366 1 18073 180 27284 -7170 106 107 0 1 1866 29 56225 322331 107 108 2964 0 42634 156 157214 1629616 108 109 0 1 249 62 54323 292754 109 110 0 1 422 28 57935 318056 110 111 36 1 2675 79 59017 355178 111 112 -8 1 965 3 73490 204325 112 113 0 1 0 0 60720 315380 113 114 -5 1 621 16 68005 317046 114 115 0 0 0 0 60720 315380 115 116 0 1 365 8 54820 309560 116 117 38 0 3122 34 94670 414462 117 118 601 0 12988 119 82340 857217 118 119 132 0 5336 81 112477 697458 119 120 66 0 3160 108 108094 530670 120 121 0 1 2489 40 67804 238125 121 122 0 1 0 0 60720 315380 122 123 195 0 6984 45 80570 741409 123 124 19 1 2001 52 95551 393343 124 125 -158 1 15236 90 77440 372631 125 126 0 0 0 0 60720 315380 126 127 0 1 530 60 73433 317291 127 128 0 1 0 0 60720 315380 128 129 -70 1 35624 223 157278 306275 129 130 13 0 1383 11 73221 317892 130 131 3 1 875 14 67000 334280 131 132 0 0 0 0 60720 315380 132 133 0 0 0 0 60720 315380 133 134 0 0 0 0 60720 315380 134 135 0 1 72 3 60398 314210 135 136 1 0 265 31 60720 306948 136 137 1 0 0 60720 315380 0 137 138 0 335 33 64175 320398 0 138 139 1 0 0 60720 315380 0 139 140 1 0 0 60720 315380 0 140 141 0 4525 226 93811 501749 101 141 142 1 3045 58 27330 202055 31 142 143 1 0 0 60720 315380 0 143 144 0 0 0 60720 315380 0 144 145 0 638 14 60370 333210 0 145 146 1 0 0 60720 315380 0 146 147 1 607 11 60436 322340 0 147 148 0 1558 63 55637 369448 0 148 149 1 1324 117 67440 291841 0 149 150 0 0 0 60720 315380 0 150 151 1 0 0 60720 315380 0 151 152 1 611 8 59190 296919 -4 152 153 1 0 0 60720 315380 0 153 154 1 923 8 58620 309038 0 154 155 0 661 3 65920 246541 -12 155 156 1 2397 14 74020 289513 6 156 157 1 366 53 61808 344425 0 157 158 1 0 0 60720 315380 0 158 159 0 135 2 62065 314210 0 159 160 1 1659 29 107577 480382 2 160 161 1 316 9 60505 315009 0 161 162 1 0 0 60720 315380 0 162 163 0 309 11 56535 312878 0 163 164 1 49 8 64107 322031 0 164 165 0 0 0 60720 315380 0 165 166 1 4519 49 102129 597793 163 166 167 1 0 0 60720 315380 0 167 168 1 837 69 61262 315688 4 168 169 1 5119 49 39039 378525 31 169 170 1 1280 117 69465 312378 26 170 171 0 2564 22 130140 403560 13 171 172 0 2045 20 97890 510834 10 172 173 1 2234 16 77200 214215 5 173 174 0 975 23 90534 235133 13 174 175 1 1136 32 48522 343613 0 175 176 1 453 21 81125 365959 0 176 177 1 0 0 60720 315380 0 177 178 1 61 14 60720 314551 0 178 179 0 368 17 61686 303230 0 179 180 1 0 0 60720 315380 0 180 181 1 4901 84 121920 469107 67 181 182 1 540 14 103960 354228 22 182 183 1 0 0 60720 315380 0 183 184 0 0 0 60720 315380 0 184 185 1 0 0 60720 315380 0 185 186 1 2 9 60735 315394 0 186 187 1 36 4 61564 312412 0 187 188 1 776 9 64230 333505 -4 188 189 1 84738 588 -26007 223193 4145 189 190 1 0 0 60720 315380 0 190 191 1 3 4 60761 315656 0 191 192 1 529 9 63870 296261 -3 192 193 1 0 0 60720 315380 0 193 194 1 405 7 60845 336425 0 194 195 1 972 34 71642 359335 0 195 196 1 0 0 60720 315380 0 196 197 0 2099 64 106611 308636 20 197 198 0 3437 47 48022 158492 9 198 199 1 0 0 60720 315380 0 199 200 0 0 0 60720 315380 0 200 201 0 22330 84 79801 711969 446 201 202 0 0 0 60720 315380 0 202 203 1 0 0 60720 315380 0 203 204 1 483 27 60830 306268 0 204 205 1 0 0 60720 315380 0 205 206 0 2239 21 88590 442882 5 206 207 1 2949 41 82903 378509 18 207 208 1 0 0 60720 315380 0 208 209 1 365 83 87192 346611 0 209 210 0 2461 57 55792 314289 0 210 211 0 21950 519 114337 856956 745 211 212 1 3294 23 75832 217193 1 212 213 1 141 6 61630 315366 0 213 214 1 572 16 58580 307930 -1 214 215 1 13326 102 165548 702380 436 215 216 0 2284 33 76403 194493 66 216 217 1 10 2 61656 316155 0 217 218 0 0 0 60720 315380 0 218 219 1 1414 198 70184 330546 13 219 220 1 1975 35 118881 394510 30 220 221 1 43 3 60887 312846 0 221 222 1 0 0 60720 315380 0 222 223 1 844 78 60925 296139 0 223 224 1 304 15 62969 295580 0 224 225 1 458 11 58625 297765 0 225 226 1 18562 155 102313 377934 143 226 227 1 0 0 60720 315380 0 227 228 1 7123 109 288170 638830 769 228 229 1 622 73 73007 304376 0 229 230 1 174 1 64820 307424 0 230 231 1 2220 22 301670 644190 154 231 232 0 121 24 56178 295370 0 232 233 1 11819 85 106113 574339 199 233 234 1 0 0 60720 315380 0 234 235 1 125 12 60798 310201 0 235 236 1 1182 13 70694 327007 -5 236 237 1 1503 39 56364 343466 6 237 238 1 0 0 60720 315380 0 238 239 1 0 0 60720 315380 0 239 240 0 30 4 62045 318098 0 240 241 1 3310 33 75230 448243 -6 241 242 0 554 11 79285 325738 30 242 243 1 0 0 60720 315380 0 243 244 1 468 21 60720 312161 0 244 245 1 4917 42 52811 243650 -3 245 246 1 3256 66 35250 407159 -8 246 247 1 0 0 60720 315380 0 247 248 1 125 13 59734 317698 0 248 249 1 0 0 60720 315380 0 249 250 0 22 2 60722 312502 0 250 251 1 0 0 60720 315380 0 251 252 0 514 52 78780 322378 0 252 253 1 0 0 60720 315380 0 253 254 0 0 0 60720 315380 0 254 255 1 0 0 60720 315380 0 255 256 1 0 0 60720 315380 0 256 257 1 0 0 60720 315380 0 257 258 1 0 0 60720 315380 0 258 259 1 10327 91 88577 640273 135 259 260 1 13718 61 96448 345783 104 260 261 1 3748 9 50350 652925 53 261 262 0 14416 361 49857 439798 220 262 263 0 1526 25 69351 278990 -4 263 264 1 666 171 117869 339836 25 264 265 1 2844 67 72683 240897 -15 265 266 1 0 0 60720 315380 0 266 267 1 368 17 61167 297141 0 267 268 1 0 0 60720 315380 0 268 269 1 333 13 70811 331323 0 269 270 1 26 1 60896 313880 0 270 271 1 0 0 60720 315380 0 271 272 1 0 0 60720 315380 0 272 273 1 1303 61 69863 309422 -5 273 274 1 20 3 60938 315245 0 274 275 0 2384 97 61348 405972 54 275 276 0 203 27 50804 300962 0 276 277 1 71 1 60745 316176 0 277 278 1 53 14 59506 302409 0 278 279 1 562 13 58456 283587 -3 279 280 1 622 7 60950 263276 -2 280 281 0 645 5 60720 312075 0 281 282 1 1763 24 61600 308336 0 282 283 1 0 0 60720 315380 0 283 284 1 317 9 63915 298700 0 284 285 1 1 6 60719 315372 0 285 286 1 275 4 59500 318745 0 286 287 1 0 0 60720 315380 0 287 288 1 936 30 67939 408881 9 288 289 0 8568 37 32168 786690 223 289 290 0 11528 70 -14545 -83265 -25 290 291 1 0 0 60720 315380 0 291 292 1 738 10 60720 321376 0 292 293 1 0 0 60720 315380 0 293 294 1 592 7 64270 276898 -4 294 295 1 126 10 60951 315547 0 295 296 1 2 1 60743 315487 0 296 297 0 0 0 60720 315380 0 297 298 0 18014 239 -1710 1405225 943 298 299 0 0 0 60720 315380 0 299 300 0 0 0 60720 315380 0 300 301 0 37704 452 60448 983660 2167 301 302 1 63 4 65688 318574 0 302 303 1 1431 34 106885 310768 30 303 304 1 94 3 61360 312887 0 304 305 1 192 4 65276 312339 0 305 306 1 32 1 59988 314964 0 306 307 1 6869 105 117520 379983 90 307 308 1 0 0 60720 315380 0 308 309 1 0 0 60720 315380 0 309 310 0 1 4 60722 315398 0 310 311 1 0 0 60720 315380 0 311 312 0 2328 94 82732 253588 0 312 313 1 209 13 64016 316647 0 313 314 1 28 2 60890 315688 0 314 315 1 176 28 68136 310670 0 315 316 0 1920 15 79420 165404 -6 316 317 1 87550 458 153198 4111912 17936 317 318 0 520 10 58650 291650 0 318 319 1 0 0 60720 315380 0 319 320 1 1013 2 69770 253468 -14 320 321 1 15 8 60831 315688 0 321 322 1 587 12 59595 325699 0 322 323 1 5371 66 87720 446211 49 323 324 1 0 0 60720 315380 0 324 325 1 1012 47 114768 368078 26 325 326 1 0 0 60720 315380 0 326 327 0 876 59 138971 352850 24 327 328 1 162556 1081 213118 6282154 56289 328 329 0 43556 247 32648 227132 929 329 330 1 3425 12 83620 283910 -3 330 331 1 810 43 74015 236761 26 331 332 0 0 0 60720 315380 0 332 333 0 0 0 60720 315380 0 333 334 0 2365 79 191778 550608 120 334 335 1 1261 33 76114 307528 13 335 336 1 0 0 60720 315380 0 336 337 1 1585 107 93099 355864 2 337 338 1 16189 295 116384 358589 379 338 339 0 0 0 60720 315380 0 339 340 1 0 0 60720 315380 0 340 341 1 10579 92 110309 375195 4 341 342 1 474 20 61977 288985 0 342 343 1 0 0 60720 315380 0 343 344 1 3642 86 90262 458343 22 344 345 0 0 0 60720 315380 0 345 346 1 472 20 80045 269587 0 346 347 1 98 2 61490 315236 0 347 348 1 3999 30 51252 42754 -58 348 349 1 0 0 60720 315380 0 349 350 1 621 48 60720 308256 0 350 351 1 0 0 60720 315380 0 351 352 0 0 0 60720 315380 0 352 353 1 30 4 60798 313164 0 353 354 1 0 0 60720 315380 0 354 355 0 746 23 71561 269661 10 355 356 1 0 0 60720 315380 0 356 357 0 0 0 60720 315380 0 357 358 1 0 0 60720 315380 0 358 359 0 0 0 60720 315380 0 359 360 0 4150 150 134759 518365 66 360 361 1 4658 158 156608 233773 58 361 362 1 814 18 39625 301881 23 362 363 0 1002 15 56750 298568 0 363 364 1 496 17 87390 325479 0 364 365 1 389 14 58990 325506 0 365 366 1 12679 110 48020 984885 1871 366 367 1 400 152 60720 313267 0 367 368 1 53 6 60349 315793 0 368 369 0 0 0 60720 315380 0 369 370 1 5109 97 67038 215362 -31 370 371 1 576 18 113761 314073 12 371 372 1 438 11 58320 298096 -5 372 373 1 165 24 62841 325176 0 373 374 0 4069 145 79804 207393 119 374 375 1 23 5 62555 314806 0 375 376 0 1285 107 99489 341340 49 376 377 0 4677 75 90131 426280 24 377 378 1 24811 505 95350 929118 822 378 379 1 167 2 64245 307322 0 379 380 1 0 0 60720 315380 0 380 381 1 4628 62 69159 387475 85 381 382 0 226 5 65745 291787 0 382 383 1 1765 20 77623 247060 28 383 384 1 460 15 63346 329784 0 384 385 1 36 3 60894 315834 0 385 386 0 989 11 58930 304555 -5 386 387 0 3055 43 60247 376641 14 387 388 1 0 0 60720 315380 0 388 389 0 0 0 60720 315380 0 389 390 1 0 0 60720 315380 0 390 391 1 0 0 60720 315380 0 391 392 1 0 0 60720 315380 0 392 393 1 13253 947 90829 357760 35 393 394 1 24 8 59818 315637 0 394 395 0 0 0 60720 315380 0 395 396 1 0 0 60720 315380 0 396 397 1 0 0 60720 315380 0 397 398 1 0 0 60720 315380 0 398 399 1 131 6 59661 327071 0 399 400 1 514 15 102725 377516 23 400 401 0 3366 80 108479 299243 12 401 402 1 9327 101 87419 387699 44 402 403 1 384 22 54683 309836 0 403 404 1 17821 373 122844 444477 643 404 405 1 397 10 62710 322327 -1 405 406 1 897 5 60720 314913 0 406 407 0 218 21 60720 315553 0 407 408 1 3369 28 87161 688779 217 408 409 1 5702 61 101481 321896 10 409 410 1 8636 855 128294 301607 108 410 411 1 534 7 62620 304485 0 411 412 0 0 0 60720 315380 0 412 413 1 0 0 60720 315380 0 413 414 1 726 3 69980 231861 -5 414 415 1 1380 26 60982 347385 25 415 416 0 180 11 59635 316386 0 416 417 1 7285 115 188873 491303 59 417 418 1 880 74 80791 261216 34 418 419 1 1 2 60727 315388 0 419 420 1 0 0 60720 315380 0 420 421 1 96 5 60379 313729 0 421 422 0 1889 45 37527 358649 14 422 423 1 45187 353 234817 1926517 2854 423 424 1 288 4 60510 296656 0 424 425 1 1270 26 69206 275311 -10 425 426 0 6526 26 55830 -42143 -28 426 427 1 0 0 60720 315380 0 427 428 1 226 12 72835 343929 0 428 429 1 694 8 68060 367655 -2 429 430 1 0 0 60720 315380 0 430 431 1 249 72 56726 313491 0 431 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Dividends Wealth -1.704e+03 -7.311e-03 -6.814e-02 8.654e-03 5.518e-05 5.647e-03 t 3.382e+00 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -3608.24 -282.40 -18.04 320.28 7407.69 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -1.704e+03 8.394e+01 -20.303 < 2e-16 *** Group -7.311e-03 5.787e-03 -1.263 0.207 Costs -6.814e-02 7.924e-03 -8.600 < 2e-16 *** Trades 8.654e-03 1.156e-03 7.483 4.23e-13 *** Dividends 5.518e-05 1.788e-04 0.309 0.758 Wealth 5.647e-03 1.654e-04 34.149 < 2e-16 *** t 3.382e+00 3.546e-01 9.539 < 2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 662.3 on 424 degrees of freedom Multiple R-squared: 0.7769, Adjusted R-squared: 0.7737 F-statistic: 246 on 6 and 424 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,] 2.272059e-05 4.544119e-05 9.999773e-01 [2,] 7.578351e-07 1.515670e-06 9.999992e-01 [3,] 2.085640e-08 4.171279e-08 1.000000e+00 [4,] 4.619261e-10 9.238523e-10 1.000000e+00 [5,] 5.626539e-11 1.125308e-10 1.000000e+00 [6,] 1.397903e-12 2.795805e-12 1.000000e+00 [7,] 3.701185e-14 7.402369e-14 1.000000e+00 [8,] 1.655761e-15 3.311522e-15 1.000000e+00 [9,] 5.977628e-17 1.195526e-16 1.000000e+00 [10,] 1.560074e-18 3.120148e-18 1.000000e+00 [11,] 3.533265e-20 7.066530e-20 1.000000e+00 [12,] 7.835376e-22 1.567075e-21 1.000000e+00 [13,] 5.362653e-23 1.072531e-22 1.000000e+00 [14,] 1.146912e-24 2.293824e-24 1.000000e+00 [15,] 9.885209e-26 1.977042e-25 1.000000e+00 [16,] 7.726178e-26 1.545236e-25 1.000000e+00 [17,] 2.576333e-27 5.152665e-27 1.000000e+00 [18,] 1.556371e-27 3.112743e-27 1.000000e+00 [19,] 5.875372e-28 1.175074e-27 1.000000e+00 [20,] 1.967372e-29 3.934745e-29 1.000000e+00 [21,] 9.187748e-31 1.837550e-30 1.000000e+00 [22,] 9.944343e-30 1.988869e-29 1.000000e+00 [23,] 4.989103e-31 9.978206e-31 1.000000e+00 [24,] 3.170604e-32 6.341208e-32 1.000000e+00 [25,] 2.531664e-33 5.063328e-33 1.000000e+00 [26,] 9.910268e-35 1.982054e-34 1.000000e+00 [27,] 1.048047e-34 2.096094e-34 1.000000e+00 [28,] 6.064536e-36 1.212907e-35 1.000000e+00 [29,] 3.774455e-37 7.548910e-37 1.000000e+00 [30,] 6.168944e-26 1.233789e-25 1.000000e+00 [31,] 8.648461e-27 1.729692e-26 1.000000e+00 [32,] 1.787419e-27 3.574839e-27 1.000000e+00 [33,] 3.008830e-28 6.017659e-28 1.000000e+00 [34,] 3.431045e-29 6.862091e-29 1.000000e+00 [35,] 4.929986e-30 9.859972e-30 1.000000e+00 [36,] 7.298154e-29 1.459631e-28 1.000000e+00 [37,] 8.418735e-30 1.683747e-29 1.000000e+00 [38,] 1.037285e-30 2.074569e-30 1.000000e+00 [39,] 1.274431e-31 2.548861e-31 1.000000e+00 [40,] 1.844527e-32 3.689053e-32 1.000000e+00 [41,] 2.344464e-33 4.688928e-33 1.000000e+00 [42,] 9.692575e-24 1.938515e-23 1.000000e+00 [43,] 2.027655e-24 4.055309e-24 1.000000e+00 [44,] 3.674027e-24 7.348055e-24 1.000000e+00 [45,] 2.344964e-23 4.689928e-23 1.000000e+00 [46,] 1.218621e-21 2.437242e-21 1.000000e+00 [47,] 2.867984e-19 5.735968e-19 1.000000e+00 [48,] 8.993969e-20 1.798794e-19 1.000000e+00 [49,] 1.018241e-19 2.036481e-19 1.000000e+00 [50,] 3.376849e-20 6.753697e-20 1.000000e+00 [51,] 9.972092e-21 1.994418e-20 1.000000e+00 [52,] 2.902458e-21 5.804916e-21 1.000000e+00 [53,] 7.762575e-22 1.552515e-21 1.000000e+00 [54,] 6.043655e-22 1.208731e-21 1.000000e+00 [55,] 1.691754e-22 3.383508e-22 1.000000e+00 [56,] 3.105629e-22 6.211258e-22 1.000000e+00 [57,] 8.688141e-23 1.737628e-22 1.000000e+00 [58,] 2.397275e-23 4.794550e-23 1.000000e+00 [59,] 8.262567e-24 1.652513e-23 1.000000e+00 [60,] 2.242724e-24 4.485449e-24 1.000000e+00 [61,] 6.141088e-25 1.228218e-24 1.000000e+00 [62,] 2.408162e-23 4.816325e-23 1.000000e+00 [63,] 7.502575e-24 1.500515e-23 1.000000e+00 [64,] 1.000000e+00 7.910006e-151 3.955003e-151 [65,] 1.000000e+00 5.659588e-150 2.829794e-150 [66,] 1.000000e+00 2.189305e-149 1.094653e-149 [67,] 1.000000e+00 1.352871e-148 6.764356e-149 [68,] 1.000000e+00 6.197943e-148 3.098972e-148 [69,] 1.000000e+00 4.056313e-147 2.028156e-147 [70,] 1.000000e+00 4.305415e-150 2.152708e-150 [71,] 1.000000e+00 6.168338e-157 3.084169e-157 [72,] 1.000000e+00 4.165760e-156 2.082880e-156 [73,] 1.000000e+00 2.778689e-155 1.389344e-155 [74,] 1.000000e+00 1.832344e-154 9.161719e-155 [75,] 1.000000e+00 6.000646e-154 3.000323e-154 [76,] 1.000000e+00 3.851203e-153 1.925602e-153 [77,] 1.000000e+00 2.322593e-152 1.161297e-152 [78,] 1.000000e+00 1.445592e-151 7.227958e-152 [79,] 1.000000e+00 8.900082e-151 4.450041e-151 [80,] 1.000000e+00 5.462174e-150 2.731087e-150 [81,] 1.000000e+00 2.881241e-149 1.440621e-149 [82,] 1.000000e+00 1.622212e-148 8.111059e-149 [83,] 1.000000e+00 9.936588e-148 4.968294e-148 [84,] 1.000000e+00 5.649209e-157 2.824605e-157 [85,] 1.000000e+00 7.711590e-158 3.855795e-158 [86,] 1.000000e+00 4.172124e-157 2.086062e-157 [87,] 1.000000e+00 2.206284e-156 1.103142e-156 [88,] 1.000000e+00 1.289439e-155 6.447195e-156 [89,] 1.000000e+00 7.208794e-155 3.604397e-155 [90,] 1.000000e+00 3.777758e-154 1.888879e-154 [91,] 1.000000e+00 2.343686e-153 1.171843e-153 [92,] 1.000000e+00 4.053710e-154 2.026855e-154 [93,] 1.000000e+00 2.268261e-153 1.134131e-153 [94,] 1.000000e+00 1.151534e-152 5.757670e-153 [95,] 1.000000e+00 5.810277e-152 2.905138e-152 [96,] 1.000000e+00 2.861782e-151 1.430891e-151 [97,] 1.000000e+00 1.253665e-164 6.268326e-165 [98,] 1.000000e+00 4.733908e-164 2.366954e-164 [99,] 1.000000e+00 0.000000e+00 0.000000e+00 [100,] 1.000000e+00 0.000000e+00 0.000000e+00 [101,] 1.000000e+00 0.000000e+00 0.000000e+00 [102,] 1.000000e+00 6.916919e-323 3.458460e-323 [103,] 1.000000e+00 2.371515e-322 1.185758e-322 [104,] 1.000000e+00 1.324096e-321 6.620480e-322 [105,] 1.000000e+00 5.138283e-321 2.569141e-321 [106,] 1.000000e+00 2.176853e-320 1.088427e-320 [107,] 1.000000e+00 9.061164e-320 4.530582e-320 [108,] 1.000000e+00 7.862857e-319 3.931429e-319 [109,] 1.000000e+00 0.000000e+00 0.000000e+00 [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 1.976263e-323 9.881313e-324 [307,] 1.000000e+00 1.262832e-320 6.314159e-321 [308,] 1.000000e+00 3.978809e-318 1.989405e-318 [309,] 1.000000e+00 1.759635e-315 8.798177e-316 [310,] 1.000000e+00 2.267308e-312 1.133654e-312 [311,] 1.000000e+00 2.591487e-309 1.295744e-309 [312,] 1.000000e+00 3.229788e-306 1.614894e-306 [313,] 1.000000e+00 3.961815e-303 1.980908e-303 [314,] 1.000000e+00 4.419303e-300 2.209651e-300 [315,] 1.000000e+00 5.139431e-297 2.569715e-297 [316,] 1.000000e+00 4.682010e-294 2.341005e-294 [317,] 1.000000e+00 5.205358e-291 2.602679e-291 [318,] 1.000000e+00 2.296224e-288 1.148112e-288 [319,] 1.000000e+00 2.503991e-285 1.251995e-285 [320,] 1.000000e+00 1.984290e-282 9.921452e-283 [321,] 1.000000e+00 2.106762e-279 1.053381e-279 [322,] 1.000000e+00 2.261731e-276 1.130865e-276 [323,] 1.000000e+00 1.205015e-273 6.025074e-274 [324,] 1.000000e+00 5.741619e-271 2.870810e-271 [325,] 1.000000e+00 1.458253e-268 7.291264e-269 [326,] 1.000000e+00 1.819026e-265 9.095128e-266 [327,] 1.000000e+00 2.270729e-262 1.135364e-262 [328,] 1.000000e+00 2.745144e-259 1.372572e-259 [329,] 1.000000e+00 3.018687e-256 1.509343e-256 [330,] 1.000000e+00 1.376102e-253 6.880510e-254 [331,] 1.000000e+00 1.723811e-250 8.619056e-251 [332,] 1.000000e+00 1.689219e-247 8.446094e-248 [333,] 1.000000e+00 1.999757e-244 9.998783e-245 [334,] 1.000000e+00 2.339448e-241 1.169724e-241 [335,] 1.000000e+00 2.625487e-238 1.312744e-238 [336,] 1.000000e+00 1.351028e-235 6.755138e-236 [337,] 1.000000e+00 1.514450e-232 7.572250e-233 [338,] 1.000000e+00 1.740538e-229 8.702692e-230 [339,] 1.000000e+00 1.334936e-226 6.674680e-227 [340,] 1.000000e+00 1.438906e-223 7.194532e-224 [341,] 1.000000e+00 1.522195e-220 7.610974e-221 [342,] 1.000000e+00 1.532065e-217 7.660323e-218 [343,] 1.000000e+00 9.616569e-215 4.808285e-215 [344,] 1.000000e+00 1.002621e-211 5.013105e-212 [345,] 1.000000e+00 1.005556e-208 5.027779e-209 [346,] 1.000000e+00 6.611857e-206 3.305928e-206 [347,] 1.000000e+00 6.835684e-203 3.417842e-203 [348,] 1.000000e+00 4.053978e-200 2.026989e-200 [349,] 1.000000e+00 4.360112e-197 2.180056e-197 [350,] 1.000000e+00 2.399277e-194 1.199639e-194 [351,] 1.000000e+00 5.242699e-192 2.621349e-192 [352,] 1.000000e+00 5.089703e-189 2.544851e-189 [353,] 1.000000e+00 5.827155e-186 2.913578e-186 [354,] 1.000000e+00 2.918213e-183 1.459106e-183 [355,] 1.000000e+00 3.657216e-180 1.828608e-180 [356,] 1.000000e+00 4.531390e-177 2.265695e-177 [357,] 1.000000e+00 6.812222e-174 3.406111e-174 [358,] 1.000000e+00 8.737520e-171 4.368760e-171 [359,] 1.000000e+00 9.808639e-168 4.904320e-168 [360,] 1.000000e+00 4.934027e-165 2.467014e-165 [361,] 1.000000e+00 4.215718e-162 2.107859e-162 [362,] 1.000000e+00 4.571114e-159 2.285557e-159 [363,] 1.000000e+00 4.793523e-156 2.396761e-156 [364,] 1.000000e+00 5.176718e-153 2.588359e-153 [365,] 1.000000e+00 4.213046e-150 2.106523e-150 [366,] 1.000000e+00 4.572068e-147 2.286034e-147 [367,] 1.000000e+00 1.644322e-144 8.221609e-145 [368,] 1.000000e+00 5.221508e-142 2.610754e-142 [369,] 1.000000e+00 7.020048e-139 3.510024e-139 [370,] 1.000000e+00 9.032079e-136 4.516040e-136 [371,] 1.000000e+00 1.115760e-132 5.578801e-133 [372,] 1.000000e+00 1.169964e-129 5.849820e-130 [373,] 1.000000e+00 4.955620e-127 2.477810e-127 [374,] 1.000000e+00 5.873417e-124 2.936709e-124 [375,] 1.000000e+00 7.239121e-121 3.619561e-121 [376,] 1.000000e+00 8.381749e-118 4.190874e-118 [377,] 1.000000e+00 4.429685e-115 2.214842e-115 [378,] 1.000000e+00 1.677709e-112 8.388545e-113 [379,] 1.000000e+00 2.339444e-109 1.169722e-109 [380,] 1.000000e+00 5.796437e-107 2.898219e-107 [381,] 1.000000e+00 9.015827e-104 4.507913e-104 [382,] 1.000000e+00 1.351119e-100 6.755594e-101 [383,] 1.000000e+00 1.938316e-97 9.691582e-98 [384,] 1.000000e+00 3.204998e-94 1.602499e-94 [385,] 1.000000e+00 4.358726e-91 2.179363e-91 [386,] 1.000000e+00 1.146473e-88 5.732367e-89 [387,] 1.000000e+00 1.769303e-85 8.846517e-86 [388,] 1.000000e+00 2.611961e-82 1.305980e-82 [389,] 1.000000e+00 3.650428e-79 1.825214e-79 [390,] 1.000000e+00 4.825812e-76 2.412906e-76 [391,] 1.000000e+00 6.994267e-73 3.497134e-73 [392,] 1.000000e+00 1.382151e-70 6.910756e-71 [393,] 1.000000e+00 1.210544e-67 6.052718e-68 [394,] 1.000000e+00 1.656306e-64 8.281532e-65 [395,] 1.000000e+00 1.024960e-61 5.124800e-62 [396,] 1.000000e+00 1.293015e-58 6.465077e-59 [397,] 1.000000e+00 1.159636e-55 5.798181e-56 [398,] 1.000000e+00 4.194382e-53 2.097191e-53 [399,] 1.000000e+00 6.932129e-50 3.466065e-50 [400,] 1.000000e+00 2.027471e-47 1.013735e-47 [401,] 1.000000e+00 3.356941e-44 1.678471e-44 [402,] 1.000000e+00 3.241516e-41 1.620758e-41 [403,] 1.000000e+00 5.115093e-39 2.557547e-39 [404,] 1.000000e+00 9.622435e-36 4.811218e-36 [405,] 1.000000e+00 2.001148e-32 1.000574e-32 [406,] 1.000000e+00 5.298584e-30 2.649292e-30 [407,] 1.000000e+00 2.053586e-29 1.026793e-29 [408,] 1.000000e+00 1.154518e-25 5.772589e-26 [409,] 1.000000e+00 5.226577e-23 2.613288e-23 [410,] 1.000000e+00 4.068235e-19 2.034118e-19 [411,] 1.000000e+00 2.907728e-15 1.453864e-15 [412,] 1.000000e+00 1.783946e-11 8.919728e-12 > postscript(file="/var/www/html/freestat/rcomp/tmp/1gi6c1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/299ox1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/399ox1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/4205z1293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5205z1293543676.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 6 2.659729 -676.106104 -70.852260 -453.286768 788.370749 108.464831 7 8 9 10 11 12 57.279986 -98.978853 -237.655282 -39.564559 -408.654166 -434.714144 13 14 15 16 17 18 -123.943076 38.580355 -85.847247 38.831116 -134.908052 -155.690244 19 20 21 22 23 24 -156.420520 -136.124575 -151.002743 -33.654495 -157.767660 -765.462298 25 26 27 28 29 30 -417.679304 -146.313213 -46.010046 -253.508974 -178.062410 -248.663445 31 32 33 34 35 36 475.785799 -188.209786 -264.297271 -1260.946942 75.634570 286.034949 37 38 39 40 41 42 -202.347834 155.803306 -683.006057 -411.921502 -831.305677 -550.004027 43 44 45 46 47 48 -437.427953 -228.799287 -106.038849 -105.005296 -238.946662 -195.979717 49 50 51 52 53 54 -272.511980 -592.886976 -58.126186 -428.949200 71.760313 -427.634042 55 56 57 58 59 60 -3273.920292 -879.579239 -251.367211 234.273161 -107.284695 -282.918621 61 62 63 64 65 66 -286.301080 -112.113328 328.419908 -298.325974 1000.056573 -297.224513 67 68 69 70 71 72 -223.629406 -357.433368 -205.999392 -316.743205 -3608.238642 525.708299 73 74 75 76 77 78 7407.693803 -330.273039 -384.816886 -401.842120 -458.997753 -330.292979 79 80 81 82 83 84 1301.255990 -232.807091 -350.842919 -357.332706 -360.715165 -395.939886 85 86 87 88 89 90 -367.480081 -370.869851 -363.121352 -377.292306 -276.155443 -489.928471 91 92 93 94 95 96 -387.782143 -271.305059 -1014.350311 697.575507 -438.584956 -486.413928 97 98 99 100 101 102 -408.069582 -390.185113 -422.152817 -181.398286 -1026.221006 -324.855106 103 104 105 106 107 108 -434.308476 -431.512751 -435.129250 2248.711243 -353.996441 -2004.181683 109 110 111 112 113 114 -304.108333 -438.485528 -462.470532 225.337444 -462.188917 -438.203599 115 116 117 118 119 120 -468.961145 -414.343232 -786.662676 -2055.013187 -2148.004728 -1423.818104 121 122 123 124 125 126 115.871148 -492.631043 -2232.352396 -786.666624 52.435141 -506.168188 127 128 129 130 131 132 -485.440219 -512.925793 1885.327634 -427.428132 -567.643527 -526.462938 133 134 135 136 137 138 -529.845396 -533.227855 -525.098109 -473.588737 699.041048 669.180981 139 140 141 142 143 144 692.276132 688.893673 435.779774 1003.370068 678.746298 674.363840 145 146 147 148 149 150 678.644899 668.598923 672.477558 717.521359 619.249950 654.069089 151 152 153 154 155 156 651.686631 667.597937 644.921714 667.355453 601.060755 539.552096 157 158 159 160 161 162 626.661359 628.009422 613.175548 220.745869 622.666697 614.479588 163 164 165 166 167 168 649.459806 578.940853 603.332213 262.481864 597.567296 600.276210 169 170 171 172 173 174 815.529869 529.092780 -2.398903 263.468489 457.636493 327.937844 175 176 177 178 179 180 684.994135 392.497738 563.742712 561.805966 552.137655 553.595337 181 182 183 184 185 186 53.299615 175.277073 543.447962 539.065503 536.683045 533.797905 187 188 189 190 191 192 523.313917 501.470316 1914.953066 519.770753 516.312762 491.299394 193 194 195 196 197 198 509.623377 507.436001 415.340292 499.476002 117.931460 638.533724 199 200 201 202 203 204 489.328627 484.946169 461.023649 478.181252 475.798793 477.338368 205 206 207 208 209 210 469.033877 234.208937 291.072723 458.886501 233.023390 515.704428 211 212 213 214 215 216 145.508295 345.644094 435.539816 462.799772 -391.380967 320.355560 217 218 219 220 221 222 420.511110 424.061917 362.700378 -72.724049 414.128000 411.532084 223 224 225 226 227 228 418.923043 389.642200 424.584371 180.081343 394.619791 -1539.739817 229 230 231 232 233 234 291.655534 350.771398 -1705.310865 419.636964 58.296675 370.942583 235 236 237 238 239 240 368.902510 286.779688 410.553557 357.412749 354.030291 338.523558 241 242 243 244 245 246 240.851150 186.285327 340.500457 342.148254 444.964180 574.047385 247 248 249 250 251 252 326.970623 333.792590 320.205706 316.261876 313.440790 159.687047 253 254 255 256 257 258 306.675873 302.293414 299.910956 296.528498 293.146039 289.763581 259 260 261 262 263 264 108.328391 76.005266 378.446539 491.129629 212.051931 -210.052686 265 266 267 268 269 270 192.114977 262.703913 260.308589 255.938997 167.672446 247.992025 271 272 273 274 275 276 245.791621 242.409163 173.945766 234.115830 244.563256 317.809359 277 278 279 280 281 282 225.823844 234.677255 245.090019 221.270061 216.205860 215.883083 283 284 285 286 287 288 205.202120 178.022298 198.862460 207.709783 191.672287 129.495951 289 290 291 292 293 294 468.887067 943.040286 178.142453 180.506251 171.377536 144.225546 295 296 297 298 299 300 164.207024 161.107985 156.847703 776.246826 150.082786 146.700327 301 302 303 304 305 306 403.022668 98.500580 -249.082522 129.661345 93.205735 134.065168 307 308 309 310 311 312 -314.206314 120.640660 117.258202 113.137321 110.493285 -57.539722 313 314 315 316 317 318 77.549631 99.198777 36.241978 -45.874763 -349.546161 109.521970 319 320 321 322 323 324 83.433618 12.772713 76.345943 87.561699 -127.476932 66.521326 325 326 327 328 329 330 -397.031304 59.756409 -613.567675 -650.774152 626.439418 -124.331880 331 332 333 334 335 336 -59.163814 38.461658 35.079200 -1093.427693 -92.073337 25.931825 337 338 339 340 341 342 -241.015666 -328.595692 14.784449 12.401991 -339.819126 1.044138 343 344 345 346 347 348 2.254616 -232.301619 -5.510301 -167.785667 -17.077865 113.928559 349 350 351 352 353 354 -18.040135 -13.218420 -24.805052 -29.187510 -31.630779 -34.952427 355 356 357 358 359 360 -123.662340 -41.717344 -46.099802 -48.482260 -52.864719 -667.971518 361 362 363 364 365 366 -839.421527 128.331359 -22.763563 -295.344637 -53.949194 87.047247 367 368 369 370 371 372 -65.525766 -78.322752 -86.689303 -94.089263 -546.014883 -70.133881 373 374 375 376 377 378 -115.272492 -223.832900 -121.323121 -430.886276 -335.213622 -238.506941 379 380 381 382 383 384 -148.216416 -122.896345 -165.704622 -170.851401 -261.438943 -155.560397 385 386 387 388 389 390 -140.871809 -120.094988 -121.673798 -149.956013 -154.338471 -156.720930 391 392 393 394 395 396 -160.103388 -163.485846 -268.533814 -161.738671 -174.633222 -177.015680 397 398 399 400 401 402 -180.398138 -183.780597 -177.277225 -552.824104 -577.338159 -357.521288 403 404 405 406 407 408 -143.837737 -596.724102 -221.472446 -203.915608 -212.207444 -441.708296 409 410 411 412 413 414 -528.293298 -687.586616 -239.212288 -232.135015 -234.517473 -307.884414 415 416 417 418 419 420 -233.595658 -234.265473 -1305.991394 -410.845816 -254.729647 -258.194682 421 422 423 424 425 426 -257.492528 -50.842785 -1525.522093 -266.495855 -335.218189 -167.801944 427 428 429 430 431 -281.871891 -389.199636 -349.409295 -292.019266 -254.007770 > postscript(file="/var/www/html/freestat/rcomp/tmp/6205z1293543676.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 2.659729 NA 1 -676.106104 2.659729 2 -70.852260 -676.106104 3 -453.286768 -70.852260 4 788.370749 -453.286768 5 108.464831 788.370749 6 57.279986 108.464831 7 -98.978853 57.279986 8 -237.655282 -98.978853 9 -39.564559 -237.655282 10 -408.654166 -39.564559 11 -434.714144 -408.654166 12 -123.943076 -434.714144 13 38.580355 -123.943076 14 -85.847247 38.580355 15 38.831116 -85.847247 16 -134.908052 38.831116 17 -155.690244 -134.908052 18 -156.420520 -155.690244 19 -136.124575 -156.420520 20 -151.002743 -136.124575 21 -33.654495 -151.002743 22 -157.767660 -33.654495 23 -765.462298 -157.767660 24 -417.679304 -765.462298 25 -146.313213 -417.679304 26 -46.010046 -146.313213 27 -253.508974 -46.010046 28 -178.062410 -253.508974 29 -248.663445 -178.062410 30 475.785799 -248.663445 31 -188.209786 475.785799 32 -264.297271 -188.209786 33 -1260.946942 -264.297271 34 75.634570 -1260.946942 35 286.034949 75.634570 36 -202.347834 286.034949 37 155.803306 -202.347834 38 -683.006057 155.803306 39 -411.921502 -683.006057 40 -831.305677 -411.921502 41 -550.004027 -831.305677 42 -437.427953 -550.004027 43 -228.799287 -437.427953 44 -106.038849 -228.799287 45 -105.005296 -106.038849 46 -238.946662 -105.005296 47 -195.979717 -238.946662 48 -272.511980 -195.979717 49 -592.886976 -272.511980 50 -58.126186 -592.886976 51 -428.949200 -58.126186 52 71.760313 -428.949200 53 -427.634042 71.760313 54 -3273.920292 -427.634042 55 -879.579239 -3273.920292 56 -251.367211 -879.579239 57 234.273161 -251.367211 58 -107.284695 234.273161 59 -282.918621 -107.284695 60 -286.301080 -282.918621 61 -112.113328 -286.301080 62 328.419908 -112.113328 63 -298.325974 328.419908 64 1000.056573 -298.325974 65 -297.224513 1000.056573 66 -223.629406 -297.224513 67 -357.433368 -223.629406 68 -205.999392 -357.433368 69 -316.743205 -205.999392 70 -3608.238642 -316.743205 71 525.708299 -3608.238642 72 7407.693803 525.708299 73 -330.273039 7407.693803 74 -384.816886 -330.273039 75 -401.842120 -384.816886 76 -458.997753 -401.842120 77 -330.292979 -458.997753 78 1301.255990 -330.292979 79 -232.807091 1301.255990 80 -350.842919 -232.807091 81 -357.332706 -350.842919 82 -360.715165 -357.332706 83 -395.939886 -360.715165 84 -367.480081 -395.939886 85 -370.869851 -367.480081 86 -363.121352 -370.869851 87 -377.292306 -363.121352 88 -276.155443 -377.292306 89 -489.928471 -276.155443 90 -387.782143 -489.928471 91 -271.305059 -387.782143 92 -1014.350311 -271.305059 93 697.575507 -1014.350311 94 -438.584956 697.575507 95 -486.413928 -438.584956 96 -408.069582 -486.413928 97 -390.185113 -408.069582 98 -422.152817 -390.185113 99 -181.398286 -422.152817 100 -1026.221006 -181.398286 101 -324.855106 -1026.221006 102 -434.308476 -324.855106 103 -431.512751 -434.308476 104 -435.129250 -431.512751 105 2248.711243 -435.129250 106 -353.996441 2248.711243 107 -2004.181683 -353.996441 108 -304.108333 -2004.181683 109 -438.485528 -304.108333 110 -462.470532 -438.485528 111 225.337444 -462.470532 112 -462.188917 225.337444 113 -438.203599 -462.188917 114 -468.961145 -438.203599 115 -414.343232 -468.961145 116 -786.662676 -414.343232 117 -2055.013187 -786.662676 118 -2148.004728 -2055.013187 119 -1423.818104 -2148.004728 120 115.871148 -1423.818104 121 -492.631043 115.871148 122 -2232.352396 -492.631043 123 -786.666624 -2232.352396 124 52.435141 -786.666624 125 -506.168188 52.435141 126 -485.440219 -506.168188 127 -512.925793 -485.440219 128 1885.327634 -512.925793 129 -427.428132 1885.327634 130 -567.643527 -427.428132 131 -526.462938 -567.643527 132 -529.845396 -526.462938 133 -533.227855 -529.845396 134 -525.098109 -533.227855 135 -473.588737 -525.098109 136 699.041048 -473.588737 137 669.180981 699.041048 138 692.276132 669.180981 139 688.893673 692.276132 140 435.779774 688.893673 141 1003.370068 435.779774 142 678.746298 1003.370068 143 674.363840 678.746298 144 678.644899 674.363840 145 668.598923 678.644899 146 672.477558 668.598923 147 717.521359 672.477558 148 619.249950 717.521359 149 654.069089 619.249950 150 651.686631 654.069089 151 667.597937 651.686631 152 644.921714 667.597937 153 667.355453 644.921714 154 601.060755 667.355453 155 539.552096 601.060755 156 626.661359 539.552096 157 628.009422 626.661359 158 613.175548 628.009422 159 220.745869 613.175548 160 622.666697 220.745869 161 614.479588 622.666697 162 649.459806 614.479588 163 578.940853 649.459806 164 603.332213 578.940853 165 262.481864 603.332213 166 597.567296 262.481864 167 600.276210 597.567296 168 815.529869 600.276210 169 529.092780 815.529869 170 -2.398903 529.092780 171 263.468489 -2.398903 172 457.636493 263.468489 173 327.937844 457.636493 174 684.994135 327.937844 175 392.497738 684.994135 176 563.742712 392.497738 177 561.805966 563.742712 178 552.137655 561.805966 179 553.595337 552.137655 180 53.299615 553.595337 181 175.277073 53.299615 182 543.447962 175.277073 183 539.065503 543.447962 184 536.683045 539.065503 185 533.797905 536.683045 186 523.313917 533.797905 187 501.470316 523.313917 188 1914.953066 501.470316 189 519.770753 1914.953066 190 516.312762 519.770753 191 491.299394 516.312762 192 509.623377 491.299394 193 507.436001 509.623377 194 415.340292 507.436001 195 499.476002 415.340292 196 117.931460 499.476002 197 638.533724 117.931460 198 489.328627 638.533724 199 484.946169 489.328627 200 461.023649 484.946169 201 478.181252 461.023649 202 475.798793 478.181252 203 477.338368 475.798793 204 469.033877 477.338368 205 234.208937 469.033877 206 291.072723 234.208937 207 458.886501 291.072723 208 233.023390 458.886501 209 515.704428 233.023390 210 145.508295 515.704428 211 345.644094 145.508295 212 435.539816 345.644094 213 462.799772 435.539816 214 -391.380967 462.799772 215 320.355560 -391.380967 216 420.511110 320.355560 217 424.061917 420.511110 218 362.700378 424.061917 219 -72.724049 362.700378 220 414.128000 -72.724049 221 411.532084 414.128000 222 418.923043 411.532084 223 389.642200 418.923043 224 424.584371 389.642200 225 180.081343 424.584371 226 394.619791 180.081343 227 -1539.739817 394.619791 228 291.655534 -1539.739817 229 350.771398 291.655534 230 -1705.310865 350.771398 231 419.636964 -1705.310865 232 58.296675 419.636964 233 370.942583 58.296675 234 368.902510 370.942583 235 286.779688 368.902510 236 410.553557 286.779688 237 357.412749 410.553557 238 354.030291 357.412749 239 338.523558 354.030291 240 240.851150 338.523558 241 186.285327 240.851150 242 340.500457 186.285327 243 342.148254 340.500457 244 444.964180 342.148254 245 574.047385 444.964180 246 326.970623 574.047385 247 333.792590 326.970623 248 320.205706 333.792590 249 316.261876 320.205706 250 313.440790 316.261876 251 159.687047 313.440790 252 306.675873 159.687047 253 302.293414 306.675873 254 299.910956 302.293414 255 296.528498 299.910956 256 293.146039 296.528498 257 289.763581 293.146039 258 108.328391 289.763581 259 76.005266 108.328391 260 378.446539 76.005266 261 491.129629 378.446539 262 212.051931 491.129629 263 -210.052686 212.051931 264 192.114977 -210.052686 265 262.703913 192.114977 266 260.308589 262.703913 267 255.938997 260.308589 268 167.672446 255.938997 269 247.992025 167.672446 270 245.791621 247.992025 271 242.409163 245.791621 272 173.945766 242.409163 273 234.115830 173.945766 274 244.563256 234.115830 275 317.809359 244.563256 276 225.823844 317.809359 277 234.677255 225.823844 278 245.090019 234.677255 279 221.270061 245.090019 280 216.205860 221.270061 281 215.883083 216.205860 282 205.202120 215.883083 283 178.022298 205.202120 284 198.862460 178.022298 285 207.709783 198.862460 286 191.672287 207.709783 287 129.495951 191.672287 288 468.887067 129.495951 289 943.040286 468.887067 290 178.142453 943.040286 291 180.506251 178.142453 292 171.377536 180.506251 293 144.225546 171.377536 294 164.207024 144.225546 295 161.107985 164.207024 296 156.847703 161.107985 297 776.246826 156.847703 298 150.082786 776.246826 299 146.700327 150.082786 300 403.022668 146.700327 301 98.500580 403.022668 302 -249.082522 98.500580 303 129.661345 -249.082522 304 93.205735 129.661345 305 134.065168 93.205735 306 -314.206314 134.065168 307 120.640660 -314.206314 308 117.258202 120.640660 309 113.137321 117.258202 310 110.493285 113.137321 311 -57.539722 110.493285 312 77.549631 -57.539722 313 99.198777 77.549631 314 36.241978 99.198777 315 -45.874763 36.241978 316 -349.546161 -45.874763 317 109.521970 -349.546161 318 83.433618 109.521970 319 12.772713 83.433618 320 76.345943 12.772713 321 87.561699 76.345943 322 -127.476932 87.561699 323 66.521326 -127.476932 324 -397.031304 66.521326 325 59.756409 -397.031304 326 -613.567675 59.756409 327 -650.774152 -613.567675 328 626.439418 -650.774152 329 -124.331880 626.439418 330 -59.163814 -124.331880 331 38.461658 -59.163814 332 35.079200 38.461658 333 -1093.427693 35.079200 334 -92.073337 -1093.427693 335 25.931825 -92.073337 336 -241.015666 25.931825 337 -328.595692 -241.015666 338 14.784449 -328.595692 339 12.401991 14.784449 340 -339.819126 12.401991 341 1.044138 -339.819126 342 2.254616 1.044138 343 -232.301619 2.254616 344 -5.510301 -232.301619 345 -167.785667 -5.510301 346 -17.077865 -167.785667 347 113.928559 -17.077865 348 -18.040135 113.928559 349 -13.218420 -18.040135 350 -24.805052 -13.218420 351 -29.187510 -24.805052 352 -31.630779 -29.187510 353 -34.952427 -31.630779 354 -123.662340 -34.952427 355 -41.717344 -123.662340 356 -46.099802 -41.717344 357 -48.482260 -46.099802 358 -52.864719 -48.482260 359 -667.971518 -52.864719 360 -839.421527 -667.971518 361 128.331359 -839.421527 362 -22.763563 128.331359 363 -295.344637 -22.763563 364 -53.949194 -295.344637 365 87.047247 -53.949194 366 -65.525766 87.047247 367 -78.322752 -65.525766 368 -86.689303 -78.322752 369 -94.089263 -86.689303 370 -546.014883 -94.089263 371 -70.133881 -546.014883 372 -115.272492 -70.133881 373 -223.832900 -115.272492 374 -121.323121 -223.832900 375 -430.886276 -121.323121 376 -335.213622 -430.886276 377 -238.506941 -335.213622 378 -148.216416 -238.506941 379 -122.896345 -148.216416 380 -165.704622 -122.896345 381 -170.851401 -165.704622 382 -261.438943 -170.851401 383 -155.560397 -261.438943 384 -140.871809 -155.560397 385 -120.094988 -140.871809 386 -121.673798 -120.094988 387 -149.956013 -121.673798 388 -154.338471 -149.956013 389 -156.720930 -154.338471 390 -160.103388 -156.720930 391 -163.485846 -160.103388 392 -268.533814 -163.485846 393 -161.738671 -268.533814 394 -174.633222 -161.738671 395 -177.015680 -174.633222 396 -180.398138 -177.015680 397 -183.780597 -180.398138 398 -177.277225 -183.780597 399 -552.824104 -177.277225 400 -577.338159 -552.824104 401 -357.521288 -577.338159 402 -143.837737 -357.521288 403 -596.724102 -143.837737 404 -221.472446 -596.724102 405 -203.915608 -221.472446 406 -212.207444 -203.915608 407 -441.708296 -212.207444 408 -528.293298 -441.708296 409 -687.586616 -528.293298 410 -239.212288 -687.586616 411 -232.135015 -239.212288 412 -234.517473 -232.135015 413 -307.884414 -234.517473 414 -233.595658 -307.884414 415 -234.265473 -233.595658 416 -1305.991394 -234.265473 417 -410.845816 -1305.991394 418 -254.729647 -410.845816 419 -258.194682 -254.729647 420 -257.492528 -258.194682 421 -50.842785 -257.492528 422 -1525.522093 -50.842785 423 -266.495855 -1525.522093 424 -335.218189 -266.495855 425 -167.801944 -335.218189 426 -281.871891 -167.801944 427 -389.199636 -281.871891 428 -349.409295 -389.199636 429 -292.019266 -349.409295 430 -254.007770 -292.019266 431 NA -254.007770 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -676.106104 2.659729 [2,] -70.852260 -676.106104 [3,] -453.286768 -70.852260 [4,] 788.370749 -453.286768 [5,] 108.464831 788.370749 [6,] 57.279986 108.464831 [7,] -98.978853 57.279986 [8,] -237.655282 -98.978853 [9,] -39.564559 -237.655282 [10,] -408.654166 -39.564559 [11,] -434.714144 -408.654166 [12,] -123.943076 -434.714144 [13,] 38.580355 -123.943076 [14,] -85.847247 38.580355 [15,] 38.831116 -85.847247 [16,] -134.908052 38.831116 [17,] -155.690244 -134.908052 [18,] -156.420520 -155.690244 [19,] -136.124575 -156.420520 [20,] -151.002743 -136.124575 [21,] -33.654495 -151.002743 [22,] -157.767660 -33.654495 [23,] -765.462298 -157.767660 [24,] -417.679304 -765.462298 [25,] -146.313213 -417.679304 [26,] -46.010046 -146.313213 [27,] -253.508974 -46.010046 [28,] -178.062410 -253.508974 [29,] -248.663445 -178.062410 [30,] 475.785799 -248.663445 [31,] -188.209786 475.785799 [32,] -264.297271 -188.209786 [33,] -1260.946942 -264.297271 [34,] 75.634570 -1260.946942 [35,] 286.034949 75.634570 [36,] -202.347834 286.034949 [37,] 155.803306 -202.347834 [38,] -683.006057 155.803306 [39,] -411.921502 -683.006057 [40,] -831.305677 -411.921502 [41,] -550.004027 -831.305677 [42,] -437.427953 -550.004027 [43,] -228.799287 -437.427953 [44,] -106.038849 -228.799287 [45,] -105.005296 -106.038849 [46,] -238.946662 -105.005296 [47,] -195.979717 -238.946662 [48,] -272.511980 -195.979717 [49,] -592.886976 -272.511980 [50,] -58.126186 -592.886976 [51,] -428.949200 -58.126186 [52,] 71.760313 -428.949200 [53,] -427.634042 71.760313 [54,] -3273.920292 -427.634042 [55,] -879.579239 -3273.920292 [56,] -251.367211 -879.579239 [57,] 234.273161 -251.367211 [58,] -107.284695 234.273161 [59,] -282.918621 -107.284695 [60,] -286.301080 -282.918621 [61,] -112.113328 -286.301080 [62,] 328.419908 -112.113328 [63,] -298.325974 328.419908 [64,] 1000.056573 -298.325974 [65,] -297.224513 1000.056573 [66,] -223.629406 -297.224513 [67,] -357.433368 -223.629406 [68,] -205.999392 -357.433368 [69,] -316.743205 -205.999392 [70,] -3608.238642 -316.743205 [71,] 525.708299 -3608.238642 [72,] 7407.693803 525.708299 [73,] -330.273039 7407.693803 [74,] -384.816886 -330.273039 [75,] -401.842120 -384.816886 [76,] -458.997753 -401.842120 [77,] -330.292979 -458.997753 [78,] 1301.255990 -330.292979 [79,] -232.807091 1301.255990 [80,] -350.842919 -232.807091 [81,] -357.332706 -350.842919 [82,] -360.715165 -357.332706 [83,] -395.939886 -360.715165 [84,] -367.480081 -395.939886 [85,] -370.869851 -367.480081 [86,] -363.121352 -370.869851 [87,] -377.292306 -363.121352 [88,] -276.155443 -377.292306 [89,] -489.928471 -276.155443 [90,] -387.782143 -489.928471 [91,] -271.305059 -387.782143 [92,] -1014.350311 -271.305059 [93,] 697.575507 -1014.350311 [94,] -438.584956 697.575507 [95,] -486.413928 -438.584956 [96,] -408.069582 -486.413928 [97,] -390.185113 -408.069582 [98,] -422.152817 -390.185113 [99,] -181.398286 -422.152817 [100,] -1026.221006 -181.398286 [101,] -324.855106 -1026.221006 [102,] -434.308476 -324.855106 [103,] -431.512751 -434.308476 [104,] -435.129250 -431.512751 [105,] 2248.711243 -435.129250 [106,] -353.996441 2248.711243 [107,] -2004.181683 -353.996441 [108,] -304.108333 -2004.181683 [109,] -438.485528 -304.108333 [110,] -462.470532 -438.485528 [111,] 225.337444 -462.470532 [112,] -462.188917 225.337444 [113,] -438.203599 -462.188917 [114,] -468.961145 -438.203599 [115,] -414.343232 -468.961145 [116,] -786.662676 -414.343232 [117,] -2055.013187 -786.662676 [118,] -2148.004728 -2055.013187 [119,] -1423.818104 -2148.004728 [120,] 115.871148 -1423.818104 [121,] -492.631043 115.871148 [122,] -2232.352396 -492.631043 [123,] -786.666624 -2232.352396 [124,] 52.435141 -786.666624 [125,] -506.168188 52.435141 [126,] -485.440219 -506.168188 [127,] -512.925793 -485.440219 [128,] 1885.327634 -512.925793 [129,] -427.428132 1885.327634 [130,] -567.643527 -427.428132 [131,] -526.462938 -567.643527 [132,] -529.845396 -526.462938 [133,] -533.227855 -529.845396 [134,] -525.098109 -533.227855 [135,] -473.588737 -525.098109 [136,] 699.041048 -473.588737 [137,] 669.180981 699.041048 [138,] 692.276132 669.180981 [139,] 688.893673 692.276132 [140,] 435.779774 688.893673 [141,] 1003.370068 435.779774 [142,] 678.746298 1003.370068 [143,] 674.363840 678.746298 [144,] 678.644899 674.363840 [145,] 668.598923 678.644899 [146,] 672.477558 668.598923 [147,] 717.521359 672.477558 [148,] 619.249950 717.521359 [149,] 654.069089 619.249950 [150,] 651.686631 654.069089 [151,] 667.597937 651.686631 [152,] 644.921714 667.597937 [153,] 667.355453 644.921714 [154,] 601.060755 667.355453 [155,] 539.552096 601.060755 [156,] 626.661359 539.552096 [157,] 628.009422 626.661359 [158,] 613.175548 628.009422 [159,] 220.745869 613.175548 [160,] 622.666697 220.745869 [161,] 614.479588 622.666697 [162,] 649.459806 614.479588 [163,] 578.940853 649.459806 [164,] 603.332213 578.940853 [165,] 262.481864 603.332213 [166,] 597.567296 262.481864 [167,] 600.276210 597.567296 [168,] 815.529869 600.276210 [169,] 529.092780 815.529869 [170,] -2.398903 529.092780 [171,] 263.468489 -2.398903 [172,] 457.636493 263.468489 [173,] 327.937844 457.636493 [174,] 684.994135 327.937844 [175,] 392.497738 684.994135 [176,] 563.742712 392.497738 [177,] 561.805966 563.742712 [178,] 552.137655 561.805966 [179,] 553.595337 552.137655 [180,] 53.299615 553.595337 [181,] 175.277073 53.299615 [182,] 543.447962 175.277073 [183,] 539.065503 543.447962 [184,] 536.683045 539.065503 [185,] 533.797905 536.683045 [186,] 523.313917 533.797905 [187,] 501.470316 523.313917 [188,] 1914.953066 501.470316 [189,] 519.770753 1914.953066 [190,] 516.312762 519.770753 [191,] 491.299394 516.312762 [192,] 509.623377 491.299394 [193,] 507.436001 509.623377 [194,] 415.340292 507.436001 [195,] 499.476002 415.340292 [196,] 117.931460 499.476002 [197,] 638.533724 117.931460 [198,] 489.328627 638.533724 [199,] 484.946169 489.328627 [200,] 461.023649 484.946169 [201,] 478.181252 461.023649 [202,] 475.798793 478.181252 [203,] 477.338368 475.798793 [204,] 469.033877 477.338368 [205,] 234.208937 469.033877 [206,] 291.072723 234.208937 [207,] 458.886501 291.072723 [208,] 233.023390 458.886501 [209,] 515.704428 233.023390 [210,] 145.508295 515.704428 [211,] 345.644094 145.508295 [212,] 435.539816 345.644094 [213,] 462.799772 435.539816 [214,] -391.380967 462.799772 [215,] 320.355560 -391.380967 [216,] 420.511110 320.355560 [217,] 424.061917 420.511110 [218,] 362.700378 424.061917 [219,] -72.724049 362.700378 [220,] 414.128000 -72.724049 [221,] 411.532084 414.128000 [222,] 418.923043 411.532084 [223,] 389.642200 418.923043 [224,] 424.584371 389.642200 [225,] 180.081343 424.584371 [226,] 394.619791 180.081343 [227,] -1539.739817 394.619791 [228,] 291.655534 -1539.739817 [229,] 350.771398 291.655534 [230,] -1705.310865 350.771398 [231,] 419.636964 -1705.310865 [232,] 58.296675 419.636964 [233,] 370.942583 58.296675 [234,] 368.902510 370.942583 [235,] 286.779688 368.902510 [236,] 410.553557 286.779688 [237,] 357.412749 410.553557 [238,] 354.030291 357.412749 [239,] 338.523558 354.030291 [240,] 240.851150 338.523558 [241,] 186.285327 240.851150 [242,] 340.500457 186.285327 [243,] 342.148254 340.500457 [244,] 444.964180 342.148254 [245,] 574.047385 444.964180 [246,] 326.970623 574.047385 [247,] 333.792590 326.970623 [248,] 320.205706 333.792590 [249,] 316.261876 320.205706 [250,] 313.440790 316.261876 [251,] 159.687047 313.440790 [252,] 306.675873 159.687047 [253,] 302.293414 306.675873 [254,] 299.910956 302.293414 [255,] 296.528498 299.910956 [256,] 293.146039 296.528498 [257,] 289.763581 293.146039 [258,] 108.328391 289.763581 [259,] 76.005266 108.328391 [260,] 378.446539 76.005266 [261,] 491.129629 378.446539 [262,] 212.051931 491.129629 [263,] -210.052686 212.051931 [264,] 192.114977 -210.052686 [265,] 262.703913 192.114977 [266,] 260.308589 262.703913 [267,] 255.938997 260.308589 [268,] 167.672446 255.938997 [269,] 247.992025 167.672446 [270,] 245.791621 247.992025 [271,] 242.409163 245.791621 [272,] 173.945766 242.409163 [273,] 234.115830 173.945766 [274,] 244.563256 234.115830 [275,] 317.809359 244.563256 [276,] 225.823844 317.809359 [277,] 234.677255 225.823844 [278,] 245.090019 234.677255 [279,] 221.270061 245.090019 [280,] 216.205860 221.270061 [281,] 215.883083 216.205860 [282,] 205.202120 215.883083 [283,] 178.022298 205.202120 [284,] 198.862460 178.022298 [285,] 207.709783 198.862460 [286,] 191.672287 207.709783 [287,] 129.495951 191.672287 [288,] 468.887067 129.495951 [289,] 943.040286 468.887067 [290,] 178.142453 943.040286 [291,] 180.506251 178.142453 [292,] 171.377536 180.506251 [293,] 144.225546 171.377536 [294,] 164.207024 144.225546 [295,] 161.107985 164.207024 [296,] 156.847703 161.107985 [297,] 776.246826 156.847703 [298,] 150.082786 776.246826 [299,] 146.700327 150.082786 [300,] 403.022668 146.700327 [301,] 98.500580 403.022668 [302,] -249.082522 98.500580 [303,] 129.661345 -249.082522 [304,] 93.205735 129.661345 [305,] 134.065168 93.205735 [306,] -314.206314 134.065168 [307,] 120.640660 -314.206314 [308,] 117.258202 120.640660 [309,] 113.137321 117.258202 [310,] 110.493285 113.137321 [311,] -57.539722 110.493285 [312,] 77.549631 -57.539722 [313,] 99.198777 77.549631 [314,] 36.241978 99.198777 [315,] -45.874763 36.241978 [316,] -349.546161 -45.874763 [317,] 109.521970 -349.546161 [318,] 83.433618 109.521970 [319,] 12.772713 83.433618 [320,] 76.345943 12.772713 [321,] 87.561699 76.345943 [322,] -127.476932 87.561699 [323,] 66.521326 -127.476932 [324,] -397.031304 66.521326 [325,] 59.756409 -397.031304 [326,] -613.567675 59.756409 [327,] -650.774152 -613.567675 [328,] 626.439418 -650.774152 [329,] -124.331880 626.439418 [330,] -59.163814 -124.331880 [331,] 38.461658 -59.163814 [332,] 35.079200 38.461658 [333,] -1093.427693 35.079200 [334,] -92.073337 -1093.427693 [335,] 25.931825 -92.073337 [336,] -241.015666 25.931825 [337,] -328.595692 -241.015666 [338,] 14.784449 -328.595692 [339,] 12.401991 14.784449 [340,] -339.819126 12.401991 [341,] 1.044138 -339.819126 [342,] 2.254616 1.044138 [343,] -232.301619 2.254616 [344,] -5.510301 -232.301619 [345,] -167.785667 -5.510301 [346,] -17.077865 -167.785667 [347,] 113.928559 -17.077865 [348,] -18.040135 113.928559 [349,] -13.218420 -18.040135 [350,] -24.805052 -13.218420 [351,] -29.187510 -24.805052 [352,] -31.630779 -29.187510 [353,] -34.952427 -31.630779 [354,] -123.662340 -34.952427 [355,] -41.717344 -123.662340 [356,] -46.099802 -41.717344 [357,] -48.482260 -46.099802 [358,] -52.864719 -48.482260 [359,] -667.971518 -52.864719 [360,] -839.421527 -667.971518 [361,] 128.331359 -839.421527 [362,] -22.763563 128.331359 [363,] -295.344637 -22.763563 [364,] -53.949194 -295.344637 [365,] 87.047247 -53.949194 [366,] -65.525766 87.047247 [367,] -78.322752 -65.525766 [368,] -86.689303 -78.322752 [369,] -94.089263 -86.689303 [370,] -546.014883 -94.089263 [371,] -70.133881 -546.014883 [372,] -115.272492 -70.133881 [373,] -223.832900 -115.272492 [374,] -121.323121 -223.832900 [375,] -430.886276 -121.323121 [376,] -335.213622 -430.886276 [377,] -238.506941 -335.213622 [378,] -148.216416 -238.506941 [379,] -122.896345 -148.216416 [380,] -165.704622 -122.896345 [381,] -170.851401 -165.704622 [382,] -261.438943 -170.851401 [383,] -155.560397 -261.438943 [384,] -140.871809 -155.560397 [385,] -120.094988 -140.871809 [386,] -121.673798 -120.094988 [387,] -149.956013 -121.673798 [388,] -154.338471 -149.956013 [389,] -156.720930 -154.338471 [390,] -160.103388 -156.720930 [391,] -163.485846 -160.103388 [392,] -268.533814 -163.485846 [393,] -161.738671 -268.533814 [394,] -174.633222 -161.738671 [395,] -177.015680 -174.633222 [396,] -180.398138 -177.015680 [397,] -183.780597 -180.398138 [398,] -177.277225 -183.780597 [399,] -552.824104 -177.277225 [400,] -577.338159 -552.824104 [401,] -357.521288 -577.338159 [402,] -143.837737 -357.521288 [403,] -596.724102 -143.837737 [404,] -221.472446 -596.724102 [405,] -203.915608 -221.472446 [406,] -212.207444 -203.915608 [407,] -441.708296 -212.207444 [408,] -528.293298 -441.708296 [409,] -687.586616 -528.293298 [410,] -239.212288 -687.586616 [411,] -232.135015 -239.212288 [412,] -234.517473 -232.135015 [413,] -307.884414 -234.517473 [414,] -233.595658 -307.884414 [415,] -234.265473 -233.595658 [416,] -1305.991394 -234.265473 [417,] -410.845816 -1305.991394 [418,] -254.729647 -410.845816 [419,] -258.194682 -254.729647 [420,] -257.492528 -258.194682 [421,] -50.842785 -257.492528 [422,] -1525.522093 -50.842785 [423,] -266.495855 -1525.522093 [424,] -335.218189 -266.495855 [425,] -167.801944 -335.218189 [426,] -281.871891 -167.801944 [427,] -389.199636 -281.871891 [428,] -349.409295 -389.199636 [429,] -292.019266 -349.409295 [430,] -254.007770 -292.019266 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -676.106104 2.659729 2 -70.852260 -676.106104 3 -453.286768 -70.852260 4 788.370749 -453.286768 5 108.464831 788.370749 6 57.279986 108.464831 7 -98.978853 57.279986 8 -237.655282 -98.978853 9 -39.564559 -237.655282 10 -408.654166 -39.564559 11 -434.714144 -408.654166 12 -123.943076 -434.714144 13 38.580355 -123.943076 14 -85.847247 38.580355 15 38.831116 -85.847247 16 -134.908052 38.831116 17 -155.690244 -134.908052 18 -156.420520 -155.690244 19 -136.124575 -156.420520 20 -151.002743 -136.124575 21 -33.654495 -151.002743 22 -157.767660 -33.654495 23 -765.462298 -157.767660 24 -417.679304 -765.462298 25 -146.313213 -417.679304 26 -46.010046 -146.313213 27 -253.508974 -46.010046 28 -178.062410 -253.508974 29 -248.663445 -178.062410 30 475.785799 -248.663445 31 -188.209786 475.785799 32 -264.297271 -188.209786 33 -1260.946942 -264.297271 34 75.634570 -1260.946942 35 286.034949 75.634570 36 -202.347834 286.034949 37 155.803306 -202.347834 38 -683.006057 155.803306 39 -411.921502 -683.006057 40 -831.305677 -411.921502 41 -550.004027 -831.305677 42 -437.427953 -550.004027 43 -228.799287 -437.427953 44 -106.038849 -228.799287 45 -105.005296 -106.038849 46 -238.946662 -105.005296 47 -195.979717 -238.946662 48 -272.511980 -195.979717 49 -592.886976 -272.511980 50 -58.126186 -592.886976 51 -428.949200 -58.126186 52 71.760313 -428.949200 53 -427.634042 71.760313 54 -3273.920292 -427.634042 55 -879.579239 -3273.920292 56 -251.367211 -879.579239 57 234.273161 -251.367211 58 -107.284695 234.273161 59 -282.918621 -107.284695 60 -286.301080 -282.918621 61 -112.113328 -286.301080 62 328.419908 -112.113328 63 -298.325974 328.419908 64 1000.056573 -298.325974 65 -297.224513 1000.056573 66 -223.629406 -297.224513 67 -357.433368 -223.629406 68 -205.999392 -357.433368 69 -316.743205 -205.999392 70 -3608.238642 -316.743205 71 525.708299 -3608.238642 72 7407.693803 525.708299 73 -330.273039 7407.693803 74 -384.816886 -330.273039 75 -401.842120 -384.816886 76 -458.997753 -401.842120 77 -330.292979 -458.997753 78 1301.255990 -330.292979 79 -232.807091 1301.255990 80 -350.842919 -232.807091 81 -357.332706 -350.842919 82 -360.715165 -357.332706 83 -395.939886 -360.715165 84 -367.480081 -395.939886 85 -370.869851 -367.480081 86 -363.121352 -370.869851 87 -377.292306 -363.121352 88 -276.155443 -377.292306 89 -489.928471 -276.155443 90 -387.782143 -489.928471 91 -271.305059 -387.782143 92 -1014.350311 -271.305059 93 697.575507 -1014.350311 94 -438.584956 697.575507 95 -486.413928 -438.584956 96 -408.069582 -486.413928 97 -390.185113 -408.069582 98 -422.152817 -390.185113 99 -181.398286 -422.152817 100 -1026.221006 -181.398286 101 -324.855106 -1026.221006 102 -434.308476 -324.855106 103 -431.512751 -434.308476 104 -435.129250 -431.512751 105 2248.711243 -435.129250 106 -353.996441 2248.711243 107 -2004.181683 -353.996441 108 -304.108333 -2004.181683 109 -438.485528 -304.108333 110 -462.470532 -438.485528 111 225.337444 -462.470532 112 -462.188917 225.337444 113 -438.203599 -462.188917 114 -468.961145 -438.203599 115 -414.343232 -468.961145 116 -786.662676 -414.343232 117 -2055.013187 -786.662676 118 -2148.004728 -2055.013187 119 -1423.818104 -2148.004728 120 115.871148 -1423.818104 121 -492.631043 115.871148 122 -2232.352396 -492.631043 123 -786.666624 -2232.352396 124 52.435141 -786.666624 125 -506.168188 52.435141 126 -485.440219 -506.168188 127 -512.925793 -485.440219 128 1885.327634 -512.925793 129 -427.428132 1885.327634 130 -567.643527 -427.428132 131 -526.462938 -567.643527 132 -529.845396 -526.462938 133 -533.227855 -529.845396 134 -525.098109 -533.227855 135 -473.588737 -525.098109 136 699.041048 -473.588737 137 669.180981 699.041048 138 692.276132 669.180981 139 688.893673 692.276132 140 435.779774 688.893673 141 1003.370068 435.779774 142 678.746298 1003.370068 143 674.363840 678.746298 144 678.644899 674.363840 145 668.598923 678.644899 146 672.477558 668.598923 147 717.521359 672.477558 148 619.249950 717.521359 149 654.069089 619.249950 150 651.686631 654.069089 151 667.597937 651.686631 152 644.921714 667.597937 153 667.355453 644.921714 154 601.060755 667.355453 155 539.552096 601.060755 156 626.661359 539.552096 157 628.009422 626.661359 158 613.175548 628.009422 159 220.745869 613.175548 160 622.666697 220.745869 161 614.479588 622.666697 162 649.459806 614.479588 163 578.940853 649.459806 164 603.332213 578.940853 165 262.481864 603.332213 166 597.567296 262.481864 167 600.276210 597.567296 168 815.529869 600.276210 169 529.092780 815.529869 170 -2.398903 529.092780 171 263.468489 -2.398903 172 457.636493 263.468489 173 327.937844 457.636493 174 684.994135 327.937844 175 392.497738 684.994135 176 563.742712 392.497738 177 561.805966 563.742712 178 552.137655 561.805966 179 553.595337 552.137655 180 53.299615 553.595337 181 175.277073 53.299615 182 543.447962 175.277073 183 539.065503 543.447962 184 536.683045 539.065503 185 533.797905 536.683045 186 523.313917 533.797905 187 501.470316 523.313917 188 1914.953066 501.470316 189 519.770753 1914.953066 190 516.312762 519.770753 191 491.299394 516.312762 192 509.623377 491.299394 193 507.436001 509.623377 194 415.340292 507.436001 195 499.476002 415.340292 196 117.931460 499.476002 197 638.533724 117.931460 198 489.328627 638.533724 199 484.946169 489.328627 200 461.023649 484.946169 201 478.181252 461.023649 202 475.798793 478.181252 203 477.338368 475.798793 204 469.033877 477.338368 205 234.208937 469.033877 206 291.072723 234.208937 207 458.886501 291.072723 208 233.023390 458.886501 209 515.704428 233.023390 210 145.508295 515.704428 211 345.644094 145.508295 212 435.539816 345.644094 213 462.799772 435.539816 214 -391.380967 462.799772 215 320.355560 -391.380967 216 420.511110 320.355560 217 424.061917 420.511110 218 362.700378 424.061917 219 -72.724049 362.700378 220 414.128000 -72.724049 221 411.532084 414.128000 222 418.923043 411.532084 223 389.642200 418.923043 224 424.584371 389.642200 225 180.081343 424.584371 226 394.619791 180.081343 227 -1539.739817 394.619791 228 291.655534 -1539.739817 229 350.771398 291.655534 230 -1705.310865 350.771398 231 419.636964 -1705.310865 232 58.296675 419.636964 233 370.942583 58.296675 234 368.902510 370.942583 235 286.779688 368.902510 236 410.553557 286.779688 237 357.412749 410.553557 238 354.030291 357.412749 239 338.523558 354.030291 240 240.851150 338.523558 241 186.285327 240.851150 242 340.500457 186.285327 243 342.148254 340.500457 244 444.964180 342.148254 245 574.047385 444.964180 246 326.970623 574.047385 247 333.792590 326.970623 248 320.205706 333.792590 249 316.261876 320.205706 250 313.440790 316.261876 251 159.687047 313.440790 252 306.675873 159.687047 253 302.293414 306.675873 254 299.910956 302.293414 255 296.528498 299.910956 256 293.146039 296.528498 257 289.763581 293.146039 258 108.328391 289.763581 259 76.005266 108.328391 260 378.446539 76.005266 261 491.129629 378.446539 262 212.051931 491.129629 263 -210.052686 212.051931 264 192.114977 -210.052686 265 262.703913 192.114977 266 260.308589 262.703913 267 255.938997 260.308589 268 167.672446 255.938997 269 247.992025 167.672446 270 245.791621 247.992025 271 242.409163 245.791621 272 173.945766 242.409163 273 234.115830 173.945766 274 244.563256 234.115830 275 317.809359 244.563256 276 225.823844 317.809359 277 234.677255 225.823844 278 245.090019 234.677255 279 221.270061 245.090019 280 216.205860 221.270061 281 215.883083 216.205860 282 205.202120 215.883083 283 178.022298 205.202120 284 198.862460 178.022298 285 207.709783 198.862460 286 191.672287 207.709783 287 129.495951 191.672287 288 468.887067 129.495951 289 943.040286 468.887067 290 178.142453 943.040286 291 180.506251 178.142453 292 171.377536 180.506251 293 144.225546 171.377536 294 164.207024 144.225546 295 161.107985 164.207024 296 156.847703 161.107985 297 776.246826 156.847703 298 150.082786 776.246826 299 146.700327 150.082786 300 403.022668 146.700327 301 98.500580 403.022668 302 -249.082522 98.500580 303 129.661345 -249.082522 304 93.205735 129.661345 305 134.065168 93.205735 306 -314.206314 134.065168 307 120.640660 -314.206314 308 117.258202 120.640660 309 113.137321 117.258202 310 110.493285 113.137321 311 -57.539722 110.493285 312 77.549631 -57.539722 313 99.198777 77.549631 314 36.241978 99.198777 315 -45.874763 36.241978 316 -349.546161 -45.874763 317 109.521970 -349.546161 318 83.433618 109.521970 319 12.772713 83.433618 320 76.345943 12.772713 321 87.561699 76.345943 322 -127.476932 87.561699 323 66.521326 -127.476932 324 -397.031304 66.521326 325 59.756409 -397.031304 326 -613.567675 59.756409 327 -650.774152 -613.567675 328 626.439418 -650.774152 329 -124.331880 626.439418 330 -59.163814 -124.331880 331 38.461658 -59.163814 332 35.079200 38.461658 333 -1093.427693 35.079200 334 -92.073337 -1093.427693 335 25.931825 -92.073337 336 -241.015666 25.931825 337 -328.595692 -241.015666 338 14.784449 -328.595692 339 12.401991 14.784449 340 -339.819126 12.401991 341 1.044138 -339.819126 342 2.254616 1.044138 343 -232.301619 2.254616 344 -5.510301 -232.301619 345 -167.785667 -5.510301 346 -17.077865 -167.785667 347 113.928559 -17.077865 348 -18.040135 113.928559 349 -13.218420 -18.040135 350 -24.805052 -13.218420 351 -29.187510 -24.805052 352 -31.630779 -29.187510 353 -34.952427 -31.630779 354 -123.662340 -34.952427 355 -41.717344 -123.662340 356 -46.099802 -41.717344 357 -48.482260 -46.099802 358 -52.864719 -48.482260 359 -667.971518 -52.864719 360 -839.421527 -667.971518 361 128.331359 -839.421527 362 -22.763563 128.331359 363 -295.344637 -22.763563 364 -53.949194 -295.344637 365 87.047247 -53.949194 366 -65.525766 87.047247 367 -78.322752 -65.525766 368 -86.689303 -78.322752 369 -94.089263 -86.689303 370 -546.014883 -94.089263 371 -70.133881 -546.014883 372 -115.272492 -70.133881 373 -223.832900 -115.272492 374 -121.323121 -223.832900 375 -430.886276 -121.323121 376 -335.213622 -430.886276 377 -238.506941 -335.213622 378 -148.216416 -238.506941 379 -122.896345 -148.216416 380 -165.704622 -122.896345 381 -170.851401 -165.704622 382 -261.438943 -170.851401 383 -155.560397 -261.438943 384 -140.871809 -155.560397 385 -120.094988 -140.871809 386 -121.673798 -120.094988 387 -149.956013 -121.673798 388 -154.338471 -149.956013 389 -156.720930 -154.338471 390 -160.103388 -156.720930 391 -163.485846 -160.103388 392 -268.533814 -163.485846 393 -161.738671 -268.533814 394 -174.633222 -161.738671 395 -177.015680 -174.633222 396 -180.398138 -177.015680 397 -183.780597 -180.398138 398 -177.277225 -183.780597 399 -552.824104 -177.277225 400 -577.338159 -552.824104 401 -357.521288 -577.338159 402 -143.837737 -357.521288 403 -596.724102 -143.837737 404 -221.472446 -596.724102 405 -203.915608 -221.472446 406 -212.207444 -203.915608 407 -441.708296 -212.207444 408 -528.293298 -441.708296 409 -687.586616 -528.293298 410 -239.212288 -687.586616 411 -232.135015 -239.212288 412 -234.517473 -232.135015 413 -307.884414 -234.517473 414 -233.595658 -307.884414 415 -234.265473 -233.595658 416 -1305.991394 -234.265473 417 -410.845816 -1305.991394 418 -254.729647 -410.845816 419 -258.194682 -254.729647 420 -257.492528 -258.194682 421 -50.842785 -257.492528 422 -1525.522093 -50.842785 423 -266.495855 -1525.522093 424 -335.218189 -266.495855 425 -167.801944 -335.218189 426 -281.871891 -167.801944 427 -389.199636 -281.871891 428 -349.409295 -389.199636 429 -292.019266 -349.409295 430 -254.007770 -292.019266 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7ds431293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/8nj351293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/9nj351293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/freestat/rcomp/tmp/10nj351293543676.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/11r2kb1293543676.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/12u20h1293543676.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/1313fb1293543676.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/14tcfe1293543676.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/15p4u41293543676.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/16bnbs1293543676.tab") + } > > try(system("convert tmp/1gi6c1293543676.ps tmp/1gi6c1293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/299ox1293543676.ps tmp/299ox1293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/399ox1293543676.ps tmp/399ox1293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/4205z1293543676.ps tmp/4205z1293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/5205z1293543676.ps tmp/5205z1293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/6205z1293543676.ps tmp/6205z1293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/7ds431293543676.ps tmp/7ds431293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/8nj351293543676.ps tmp/8nj351293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/9nj351293543676.ps tmp/9nj351293543676.png",intern=TRUE)) character(0) > try(system("convert tmp/10nj351293543676.ps tmp/10nj351293543676.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 14.836 3.143 15.517