R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1595 + ,17 + ,60720 + ,319369 + ,0 + ,5565 + ,47 + ,94355 + ,493408 + ,39 + ,601 + ,6 + ,60720 + ,319210 + ,0 + ,188 + ,11 + ,77655 + ,381180 + ,0 + ,7146 + ,105 + ,134028 + ,297978 + ,305 + ,1135 + ,17 + ,62285 + ,290476 + ,-9 + ,450 + ,8 + ,59325 + ,292136 + ,-1 + ,34 + ,8 + ,60630 + ,314353 + ,0 + ,133 + ,14 + ,65990 + ,339445 + ,0 + ,119 + ,6 + ,59118 + ,303677 + ,0 + ,2053 + ,53 + ,100423 + ,397144 + ,33 + ,4036 + ,63 + ,100269 + ,424898 + ,32 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4655 + ,393 + ,60720 + ,341570 + ,0 + ,131 + ,11 + ,61808 + ,308989 + ,0 + ,1766 + ,73 + ,60438 + ,305959 + ,0 + ,312 + ,14 + ,58598 + ,318690 + ,0 + ,448 + ,40 + ,59781 + ,323361 + ,0 + ,115 + ,13 + ,60945 + ,318903 + ,0 + ,60 + ,10 + ,61124 + ,314049 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,364 + ,35 + ,47705 + ,298466 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1442 + ,118 + ,121173 + ,438493 + ,-3 + ,1389 + ,9 + ,57530 + ,378049 + ,6 + ,149 + ,5 + ,62041 + ,313332 + ,0 + ,2212 + ,14 + ,60720 + ,319864 + ,0 + ,7489 + ,322 + ,136996 + ,430866 + ,70 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,402 + ,26 + ,84990 + ,332743 + ,5 + ,7419 + ,29 + ,63255 + ,286963 + ,-5 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,307 + ,15 + ,58856 + ,331955 + ,0 + ,1134 + ,162 + ,146216 + ,527021 + ,58 + ,7561 + ,87 + ,82425 + ,364304 + ,37 + ,5131 + ,71 + ,86111 + ,292154 + ,9 + ,9 + ,6 + ,60835 + ,314987 + ,0 + ,2264 + ,24 + ,55174 + ,272713 + ,-31 + ,40949 + ,481 + ,129352 + ,1398893 + ,2865 + ,4980 + ,91 + ,113521 + ,409642 + ,0 + ,272 + ,8 + ,112995 + ,429112 + ,14 + ,757 + ,19 + ,45689 + ,382712 + ,0 + ,1424 + ,94 + ,131116 + ,374943 + ,32 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,378 + ,52 + ,71873 + ,297413 + ,0 + ,1061 + ,5 + ,70415 + ,304252 + ,-4 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,3689 + ,25 + ,101792 + ,377305 + ,74 + ,567 + ,11 + ,103345 + ,386212 + ,20 + ,27664 + ,282 + ,74996 + ,657954 + ,247 + ,1686 + ,40 + ,68696 + ,368186 + ,11 + ,1164 + ,17 + ,57320 + ,269753 + ,-6 + ,15824 + ,64 + ,63838 + ,585715 + ,284 + ,6575 + ,159 + ,151352 + ,1071292 + ,819 + ,37597 + ,248 + ,37884 + ,992426 + ,652 + ,314 + ,5 + ,60720 + ,315371 + ,0 + ,1932 + ,36 + ,72954 + ,249898 + ,10 + ,5692 + ,72 + ,64239 + ,357312 + ,22 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,6533 + ,75 + ,126910 + ,364839 + ,16 + ,2126 + ,138 + ,21001 + ,230621 + ,-3 + ,14 + ,1 + ,61016 + ,315877 + ,0 + ,3192 + ,60 + ,59831 + ,125390 + ,10 + ,47 + ,3 + ,60720 + ,314882 + ,0 + ,4891 + ,217 + ,210568 + ,370837 + ,73 + ,531 + ,10 + ,71641 + ,330068 + ,0 + ,2332 + ,4 + ,72680 + ,324385 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,16275 + ,109 + ,47818 + ,1443586 + ,1974 + ,9252 + ,48 + ,79510 + ,253537 + ,-129 + ,29790 + ,309 + ,81767 + ,4321023 + ,28328 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2300 + ,29 + ,64873 + ,352108 + ,0 + ,267 + ,32 + ,57640 + ,330059 + ,0 + ,382 + ,21 + ,59365 + ,340968 + ,0 + ,397 + ,28 + ,60646 + ,317736 + ,0 + ,11920 + ,182 + ,31080 + ,209458 + ,238 + ,54660 + ,302 + ,126942 + ,1491348 + ,3040 + ,5 + ,2 + ,60720 + ,314887 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2647 + ,37 + ,62920 + ,353058 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,94 + ,7 + ,60793 + ,314533 + ,0 + ,3 + ,2 + ,60698 + ,315354 + ,0 + ,454 + ,51 + ,63261 + ,302187 + ,0 + ,227 + ,9 + ,76644 + ,336639 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,25830 + ,115 + ,245546 + ,1073089 + ,1910 + ,1528 + ,38 + ,60326 + ,146494 + ,38 + ,271 + ,2 + ,60720 + ,325249 + ,0 + ,138 + ,6 + ,69817 + ,331420 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,278 + ,26 + ,61404 + ,314922 + ,0 + ,282 + ,30 + ,62452 + ,320016 + ,0 + ,571 + ,80 + ,41477 + ,280398 + ,0 + ,2253 + ,99 + ,63593 + ,452469 + ,17 + ,290 + ,2 + ,58790 + ,301164 + ,0 + ,78 + ,16 + ,62700 + ,317330 + ,0 + ,20 + ,2 + ,60805 + ,315576 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,18073 + ,180 + ,27284 + ,-7170 + ,-366 + ,1866 + ,29 + ,56225 + ,322331 + ,0 + ,42634 + ,156 + ,157214 + ,1629616 + ,2964 + ,249 + ,62 + ,54323 + ,292754 + ,0 + ,422 + ,28 + ,57935 + ,318056 + ,0 + ,2675 + ,79 + ,59017 + ,355178 + ,36 + ,965 + ,3 + ,73490 + ,204325 + ,-8 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,621 + ,16 + ,68005 + ,317046 + ,-5 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,8 + ,54820 + ,309560 + ,0 + ,3122 + ,34 + ,94670 + ,414462 + ,38 + ,12988 + ,119 + ,82340 + ,857217 + ,601 + ,5336 + ,81 + ,112477 + ,697458 + ,132 + ,3160 + ,108 + ,108094 + ,530670 + ,66 + ,2489 + ,40 + ,67804 + ,238125 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,6984 + ,45 + ,80570 + ,741409 + ,195 + ,2001 + ,52 + ,95551 + ,393343 + ,19 + ,15236 + ,90 + ,77440 + ,372631 + ,-158 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,530 + ,60 + ,73433 + ,317291 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,35624 + ,223 + ,157278 + ,306275 + ,-70 + ,1383 + ,11 + ,73221 + ,317892 + ,13 + ,875 + ,14 + ,67000 + ,334280 + ,3 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,72 + ,3 + ,60398 + ,314210 + ,0 + ,265 + ,31 + ,60720 + ,306948 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,335 + ,33 + ,64175 + ,320398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4525 + ,226 + ,93811 + ,501749 + ,101 + ,3045 + ,58 + ,27330 + ,202055 + ,31 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,638 + ,14 + ,60370 + ,333210 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,607 + ,11 + ,60436 + ,322340 + ,0 + ,1558 + ,63 + ,55637 + ,369448 + ,0 + ,1324 + ,117 + ,67440 + ,291841 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,611 + ,8 + ,59190 + ,296919 + ,-4 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,923 + ,8 + ,58620 + ,309038 + ,0 + ,661 + ,3 + ,65920 + ,246541 + ,-12 + ,2397 + ,14 + ,74020 + ,289513 + ,6 + ,366 + ,53 + ,61808 + ,344425 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,135 + ,2 + ,62065 + ,314210 + ,0 + ,1659 + ,29 + ,107577 + ,480382 + ,2 + ,316 + ,9 + ,60505 + ,315009 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,309 + ,11 + ,56535 + ,312878 + ,0 + ,49 + ,8 + ,64107 + ,322031 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4519 + ,49 + ,102129 + ,597793 + ,163 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,837 + ,69 + ,61262 + ,315688 + ,4 + ,5119 + ,49 + ,39039 + ,378525 + ,31 + ,1280 + ,117 + ,69465 + ,312378 + ,26 + ,2564 + ,22 + ,130140 + ,403560 + ,13 + ,2045 + ,20 + ,97890 + ,510834 + ,10 + ,2234 + ,16 + ,77200 + ,214215 + ,5 + ,975 + ,23 + ,90534 + ,235133 + ,13 + ,1136 + ,32 + ,48522 + ,343613 + ,0 + ,453 + ,21 + ,81125 + ,365959 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,61 + ,14 + ,60720 + ,314551 + ,0 + ,368 + ,17 + ,61686 + ,303230 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4901 + ,84 + ,121920 + ,469107 + ,67 + ,540 + ,14 + ,103960 + ,354228 + ,22 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2 + ,9 + ,60735 + ,315394 + ,0 + ,36 + ,4 + ,61564 + ,312412 + ,0 + ,776 + ,9 + ,64230 + ,333505 + ,-4 + ,84738 + ,588 + ,-26007 + ,223193 + ,4145 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,3 + ,4 + ,60761 + ,315656 + ,0 + ,529 + ,9 + ,63870 + ,296261 + ,-3 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,405 + ,7 + ,60845 + ,336425 + ,0 + ,972 + ,34 + ,71642 + ,359335 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2099 + ,64 + ,106611 + ,308636 + ,20 + ,3437 + ,47 + ,48022 + ,158492 + ,9 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,22330 + ,84 + ,79801 + ,711969 + ,446 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,483 + ,27 + ,60830 + ,306268 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2239 + ,21 + ,88590 + ,442882 + ,5 + ,2949 + ,41 + ,82903 + ,378509 + ,18 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,83 + ,87192 + ,346611 + ,0 + ,2461 + ,57 + ,55792 + ,314289 + ,0 + ,21950 + ,519 + ,114337 + ,856956 + ,745 + ,3294 + ,23 + ,75832 + ,217193 + ,1 + ,141 + ,6 + ,61630 + ,315366 + ,0 + ,572 + ,16 + ,58580 + ,307930 + ,-1 + ,13326 + ,102 + ,165548 + ,702380 + ,436 + ,2284 + ,33 + ,76403 + ,194493 + ,66 + ,10 + ,2 + ,61656 + ,316155 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1414 + ,198 + ,70184 + ,330546 + ,13 + ,1975 + ,35 + ,118881 + ,394510 + ,30 + ,43 + ,3 + ,60887 + ,312846 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,844 + ,78 + ,60925 + ,296139 + ,0 + ,304 + ,15 + ,62969 + ,295580 + ,0 + ,458 + ,11 + ,58625 + ,297765 + ,0 + ,18562 + ,155 + ,102313 + ,377934 + ,143 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,7123 + ,109 + ,288170 + ,638830 + ,769 + ,622 + ,73 + ,73007 + ,304376 + ,0 + ,174 + ,1 + ,64820 + ,307424 + ,0 + ,2220 + ,22 + ,301670 + ,644190 + ,154 + ,121 + ,24 + ,56178 + ,295370 + ,0 + ,11819 + ,85 + ,106113 + ,574339 + ,199 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,125 + ,12 + ,60798 + ,310201 + ,0 + ,1182 + ,13 + ,70694 + ,327007 + ,-5 + ,1503 + ,39 + ,56364 + ,343466 + ,6 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,30 + ,4 + ,62045 + ,318098 + ,0 + ,3310 + ,33 + ,75230 + ,448243 + ,-6 + ,554 + ,11 + ,79285 + ,325738 + ,30 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,468 + ,21 + ,60720 + ,312161 + ,0 + ,4917 + ,42 + ,52811 + ,243650 + ,-3 + ,3256 + ,66 + ,35250 + ,407159 + ,-8 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,125 + ,13 + ,59734 + ,317698 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,22 + ,2 + ,60722 + ,312502 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,514 + ,52 + ,78780 + ,322378 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,10327 + ,91 + ,88577 + ,640273 + ,135 + ,13718 + ,61 + ,96448 + ,345783 + ,104 + ,3748 + ,9 + ,50350 + ,652925 + ,53 + ,14416 + ,361 + ,49857 + ,439798 + ,220 + ,1526 + ,25 + ,69351 + ,278990 + ,-4 + ,666 + ,171 + ,117869 + ,339836 + ,25 + ,2844 + ,67 + ,72683 + ,240897 + ,-15 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,368 + ,17 + ,61167 + ,297141 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,333 + ,13 + ,70811 + ,331323 + ,0 + ,26 + ,1 + ,60896 + ,313880 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1303 + ,61 + ,69863 + ,309422 + ,-5 + ,20 + ,3 + ,60938 + ,315245 + ,0 + ,2384 + ,97 + ,61348 + ,405972 + ,54 + ,203 + ,27 + ,50804 + ,300962 + ,0 + ,71 + ,1 + ,60745 + ,316176 + ,0 + ,53 + ,14 + ,59506 + ,302409 + ,0 + ,562 + ,13 + ,58456 + ,283587 + ,-3 + ,622 + ,7 + ,60950 + ,263276 + ,-2 + ,645 + ,5 + ,60720 + ,312075 + ,0 + ,1763 + ,24 + ,61600 + ,308336 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,317 + ,9 + ,63915 + ,298700 + ,0 + ,1 + ,6 + ,60719 + ,315372 + ,0 + ,275 + ,4 + ,59500 + ,318745 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,936 + ,30 + ,67939 + ,408881 + ,9 + ,8568 + ,37 + ,32168 + ,786690 + ,223 + ,11528 + ,70 + ,-14545 + ,-83265 + ,-25 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,738 + ,10 + ,60720 + ,321376 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,592 + ,7 + ,64270 + ,276898 + ,-4 + ,126 + ,10 + ,60951 + ,315547 + ,0 + ,2 + ,1 + ,60743 + ,315487 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,18014 + ,239 + ,-1710 + ,1405225 + ,943 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,37704 + ,452 + ,60448 + ,983660 + ,2167 + ,63 + ,4 + ,65688 + ,318574 + ,0 + ,1431 + ,34 + ,106885 + ,310768 + ,30 + ,94 + ,3 + ,61360 + ,312887 + ,0 + ,192 + ,4 + ,65276 + ,312339 + ,0 + ,32 + ,1 + ,59988 + ,314964 + ,0 + ,6869 + ,105 + ,117520 + ,379983 + ,90 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2328 + ,94 + ,82732 + ,253588 + ,0 + ,209 + ,13 + ,64016 + ,316647 + ,0 + ,28 + ,2 + ,60890 + ,315688 + ,0 + ,176 + ,28 + ,68136 + ,310670 + ,0 + ,1920 + ,15 + ,79420 + ,165404 + ,-6 + ,87550 + ,458 + ,153198 + ,4111912 + ,17936 + ,520 + ,10 + ,58650 + ,291650 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1013 + ,2 + ,69770 + ,253468 + ,-14 + ,15 + ,8 + ,60831 + ,315688 + ,0 + ,587 + ,12 + ,59595 + ,325699 + ,0 + ,5371 + ,66 + ,87720 + ,446211 + ,49 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1012 + ,47 + ,114768 + ,368078 + ,26 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,876 + ,59 + ,138971 + ,352850 + ,24 + ,162556 + ,1081 + ,213118 + ,6282154 + ,56289 + ,43556 + ,247 + ,32648 + ,227132 + ,929 + ,3425 + ,12 + ,83620 + ,283910 + ,-3 + ,810 + ,43 + ,74015 + ,236761 + ,26 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2365 + ,79 + ,191778 + ,550608 + ,120 + ,1261 + ,33 + ,76114 + ,307528 + ,13 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1585 + ,107 + ,93099 + ,355864 + ,2 + ,16189 + ,295 + ,116384 + ,358589 + ,379 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,10579 + ,92 + ,110309 + ,375195 + ,4 + ,474 + ,20 + ,61977 + ,288985 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,3642 + ,86 + ,90262 + ,458343 + ,22 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,472 + ,20 + ,80045 + ,269587 + ,0 + ,98 + ,2 + ,61490 + ,315236 + ,0 + ,3999 + ,30 + ,51252 + ,42754 + ,-58 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,621 + ,48 + ,60720 + ,308256 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,746 + ,23 + ,71561 + ,269661 + ,10 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4150 + ,150 + ,134759 + ,518365 + ,66 + ,4658 + ,158 + ,156608 + ,233773 + ,58 + ,814 + ,18 + ,39625 + ,301881 + ,23 + ,1002 + ,15 + ,56750 + ,298568 + ,0 + ,496 + ,17 + ,87390 + ,325479 + ,0 + ,389 + ,14 + ,58990 + ,325506 + ,0 + ,12679 + ,110 + ,48020 + ,984885 + ,1871 + ,400 + ,152 + ,60720 + ,313267 + ,0 + ,53 + ,6 + ,60349 + ,315793 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,5109 + ,97 + ,67038 + ,215362 + ,-31 + ,576 + ,18 + ,113761 + ,314073 + ,12 + ,438 + ,11 + ,58320 + ,298096 + ,-5 + ,165 + ,24 + ,62841 + ,325176 + ,0 + ,4069 + ,145 + ,79804 + ,207393 + ,119 + ,23 + ,5 + ,62555 + ,314806 + ,0 + ,1285 + ,107 + ,99489 + ,341340 + ,49 + ,4677 + ,75 + ,90131 + ,426280 + ,24 + ,24811 + ,505 + ,95350 + ,929118 + ,822 + ,167 + ,2 + ,64245 + ,307322 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4628 + ,62 + ,69159 + ,387475 + ,85 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,1765 + ,20 + ,77623 + ,247060 + ,28 + ,460 + ,15 + ,63346 + ,329784 + ,0 + ,36 + ,3 + ,60894 + ,315834 + ,0 + ,989 + ,11 + ,58930 + ,304555 + ,-5 + ,3055 + ,43 + ,60247 + ,376641 + ,14 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,13253 + ,947 + ,90829 + ,357760 + ,35 + ,24 + ,8 + ,59818 + ,315637 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,131 + ,6 + ,59661 + ,327071 + ,0 + ,514 + ,15 + ,102725 + ,377516 + ,23 + ,3366 + ,80 + ,108479 + ,299243 + ,12 + ,9327 + ,101 + ,87419 + ,387699 + ,44 + ,384 + ,22 + ,54683 + ,309836 + ,0 + ,17821 + ,373 + ,122844 + ,444477 + ,643 + ,397 + ,10 + ,62710 + ,322327 + ,-1 + ,897 + ,5 + ,60720 + ,314913 + ,0 + ,218 + ,21 + ,60720 + ,315553 + ,0 + ,3369 + ,28 + ,87161 + ,688779 + ,217 + ,5702 + ,61 + ,101481 + ,321896 + ,10 + ,8636 + ,855 + ,128294 + ,301607 + ,108 + ,534 + ,7 + ,62620 + ,304485 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,726 + ,3 + ,69980 + ,231861 + ,-5 + ,1380 + ,26 + ,60982 + ,347385 + ,25 + ,180 + ,11 + ,59635 + ,316386 + ,0 + ,7285 + ,115 + ,188873 + ,491303 + ,59 + ,880 + ,74 + ,80791 + ,261216 + ,34 + ,1 + ,2 + ,60727 + ,315388 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,96 + ,5 + ,60379 + ,313729 + ,0 + ,1889 + ,45 + ,37527 + ,358649 + ,14 + ,45187 + ,353 + ,234817 + ,1926517 + ,2854 + ,288 + ,4 + ,60510 + ,296656 + ,0 + ,1270 + ,26 + ,69206 + ,275311 + ,-10 + ,6526 + ,26 + ,55830 + ,-42143 + ,-28 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,226 + ,12 + ,72835 + ,343929 + ,0 + ,694 + ,8 + ,68060 + ,367655 + ,-2 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,249 + ,72 + ,56726 + ,313491 + ,0) + ,dim=c(5 + ,431) + ,dimnames=list(c('Costs' + ,'Trades' + ,'Dividends' + ,'Wealth' + ,'CCScore ') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('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]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '5' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x CCScore\r Costs Trades Dividends Wealth 1 0 1595 17 60720 319369 2 39 5565 47 94355 493408 3 0 601 6 60720 319210 4 0 188 11 77655 381180 5 305 7146 105 134028 297978 6 -9 1135 17 62285 290476 7 -1 450 8 59325 292136 8 0 34 8 60630 314353 9 0 133 14 65990 339445 10 0 119 6 59118 303677 11 33 2053 53 100423 397144 12 32 4036 63 100269 424898 13 0 0 0 60720 315380 14 0 4655 393 60720 341570 15 0 131 11 61808 308989 16 0 1766 73 60438 305959 17 0 312 14 58598 318690 18 0 448 40 59781 323361 19 0 115 13 60945 318903 20 0 60 10 61124 314049 21 0 0 0 60720 315380 22 0 364 35 47705 298466 23 0 0 0 60720 315380 24 -3 1442 118 121173 438493 25 6 1389 9 57530 378049 26 0 149 5 62041 313332 27 0 2212 14 60720 319864 28 70 7489 322 136996 430866 29 0 0 0 60720 315380 30 5 402 26 84990 332743 31 -5 7419 29 63255 286963 32 0 0 0 60720 315380 33 0 307 15 58856 331955 34 58 1134 162 146216 527021 35 37 7561 87 82425 364304 36 9 5131 71 86111 292154 37 0 9 6 60835 314987 38 -31 2264 24 55174 272713 39 2865 40949 481 129352 1398893 40 0 4980 91 113521 409642 41 14 272 8 112995 429112 42 0 757 19 45689 382712 43 32 1424 94 131116 374943 44 0 0 0 60720 315380 45 0 378 52 71873 297413 46 -4 1061 5 70415 304252 47 0 0 0 60720 315380 48 0 203 9 61938 309596 49 74 3689 25 101792 377305 50 20 567 11 103345 386212 51 247 27664 282 74996 657954 52 11 1686 40 68696 368186 53 -6 1164 17 57320 269753 54 284 15824 64 63838 585715 55 819 6575 159 151352 1071292 56 652 37597 248 37884 992426 57 0 314 5 60720 315371 58 10 1932 36 72954 249898 59 22 5692 72 64239 357312 60 0 0 0 60720 315380 61 0 0 0 60720 315380 62 16 6533 75 126910 364839 63 -3 2126 138 21001 230621 64 0 14 1 61016 315877 65 10 3192 60 59831 125390 66 0 47 3 60720 314882 67 73 4891 217 210568 370837 68 0 531 10 71641 330068 69 0 2332 4 72680 324385 70 0 0 0 60720 315380 71 1974 16275 109 47818 1443586 72 -129 9252 48 79510 253537 73 28328 29790 309 81767 4321023 74 0 0 0 60720 315380 75 0 2300 29 64873 352108 76 0 267 32 57640 330059 77 0 382 21 59365 340968 78 0 397 28 60646 317736 79 238 11920 182 31080 209458 80 3040 54660 302 126942 1491348 81 0 5 2 60720 314887 82 0 0 0 60720 315380 83 0 0 0 60720 315380 84 1 2647 37 62920 353058 85 0 0 0 60720 315380 86 0 0 0 60720 315380 87 0 94 7 60793 314533 88 0 3 2 60698 315354 89 0 454 51 63261 302187 90 0 227 9 76644 336639 91 0 0 0 60720 315380 92 0 206 9 53110 296702 93 1910 25830 115 245546 1073089 94 38 1528 38 60326 146494 95 0 271 2 60720 325249 96 0 138 6 69817 331420 97 0 0 0 60720 315380 98 0 278 26 61404 314922 99 0 282 30 62452 320016 100 0 571 80 41477 280398 101 17 2253 99 63593 452469 102 0 290 2 58790 301164 103 0 78 16 62700 317330 104 0 20 2 60805 315576 105 0 0 0 60720 315380 106 -366 18073 180 27284 -7170 107 0 1866 29 56225 322331 108 2964 42634 156 157214 1629616 109 0 249 62 54323 292754 110 0 422 28 57935 318056 111 36 2675 79 59017 355178 112 -8 965 3 73490 204325 113 0 0 0 60720 315380 114 -5 621 16 68005 317046 115 0 0 0 60720 315380 116 0 365 8 54820 309560 117 38 3122 34 94670 414462 118 601 12988 119 82340 857217 119 132 5336 81 112477 697458 120 66 3160 108 108094 530670 121 0 2489 40 67804 238125 122 0 0 0 60720 315380 123 195 6984 45 80570 741409 124 19 2001 52 95551 393343 125 -158 15236 90 77440 372631 126 0 0 0 60720 315380 127 0 530 60 73433 317291 128 0 0 0 60720 315380 129 -70 35624 223 157278 306275 130 13 1383 11 73221 317892 131 3 875 14 67000 334280 132 0 0 0 60720 315380 133 0 0 0 60720 315380 134 0 0 0 60720 315380 135 0 72 3 60398 314210 136 0 265 31 60720 306948 137 335 0 60720 315380 0 138 0 33 64175 320398 0 139 0 0 60720 315380 0 140 4525 0 60720 315380 0 141 3045 226 93811 501749 101 142 0 58 27330 202055 31 143 0 0 60720 315380 0 144 638 0 60720 315380 0 145 0 14 60370 333210 0 146 607 0 60720 315380 0 147 1558 11 60436 322340 0 148 1324 63 55637 369448 0 149 0 117 67440 291841 0 150 0 0 60720 315380 0 151 611 0 60720 315380 0 152 0 8 59190 296919 -4 153 923 0 60720 315380 0 154 661 8 58620 309038 0 155 2397 3 65920 246541 -12 156 366 14 74020 289513 6 157 0 53 61808 344425 0 158 135 0 60720 315380 0 159 1659 2 62065 314210 0 160 316 29 107577 480382 2 161 0 9 60505 315009 0 162 309 0 60720 315380 0 163 49 11 56535 312878 0 164 0 8 64107 322031 0 165 4519 0 60720 315380 0 166 0 49 102129 597793 163 167 837 0 60720 315380 0 168 5119 69 61262 315688 4 169 1280 49 39039 378525 31 170 2564 117 69465 312378 26 171 2045 22 130140 403560 13 172 2234 20 97890 510834 10 173 975 16 77200 214215 5 174 1136 23 90534 235133 13 175 453 32 48522 343613 0 176 0 21 81125 365959 0 177 61 0 60720 315380 0 178 368 14 60720 314551 0 179 0 17 61686 303230 0 180 4901 0 60720 315380 0 181 540 84 121920 469107 67 182 0 14 103960 354228 22 183 0 0 60720 315380 0 184 0 0 60720 315380 0 185 2 0 60720 315380 0 186 36 9 60735 315394 0 187 776 4 61564 312412 0 188 84738 9 64230 333505 -4 189 0 588 -26007 223193 4145 190 3 0 60720 315380 0 191 529 4 60761 315656 0 192 0 9 63870 296261 -3 193 405 0 60720 315380 0 194 972 7 60845 336425 0 195 0 34 71642 359335 0 196 2099 0 60720 315380 0 197 3437 64 106611 308636 20 198 0 47 48022 158492 9 199 0 0 60720 315380 0 200 22330 0 60720 315380 0 201 0 84 79801 711969 446 202 0 0 60720 315380 0 203 483 0 60720 315380 0 204 0 27 60830 306268 0 205 2239 0 60720 315380 0 206 2949 21 88590 442882 5 207 0 41 82903 378509 18 208 365 0 60720 315380 0 209 2461 83 87192 346611 0 210 21950 57 55792 314289 0 211 3294 519 114337 856956 745 212 141 23 75832 217193 1 213 572 6 61630 315366 0 214 13326 16 58580 307930 -1 215 2284 102 165548 702380 436 216 10 33 76403 194493 66 217 0 2 61656 316155 0 218 1414 0 60720 315380 0 219 1975 198 70184 330546 13 220 43 35 118881 394510 30 221 0 3 60887 312846 0 222 844 0 60720 315380 0 223 304 78 60925 296139 0 224 458 15 62969 295580 0 225 18562 11 58625 297765 0 226 0 155 102313 377934 143 227 7123 0 60720 315380 0 228 622 109 288170 638830 769 229 174 73 73007 304376 0 230 2220 1 64820 307424 0 231 121 22 301670 644190 154 232 11819 24 56178 295370 0 233 0 85 106113 574339 199 234 125 0 60720 315380 0 235 1182 12 60798 310201 0 236 1503 13 70694 327007 -5 237 0 39 56364 343466 6 238 0 0 60720 315380 0 239 30 0 60720 315380 0 240 3310 4 62045 318098 0 241 554 33 75230 448243 -6 242 0 11 79285 325738 30 243 468 0 60720 315380 0 244 4917 21 60720 312161 0 245 3256 42 52811 243650 -3 246 0 66 35250 407159 -8 247 125 0 60720 315380 0 248 0 13 59734 317698 0 249 22 0 60720 315380 0 250 0 2 60722 312502 0 251 514 0 60720 315380 0 252 0 52 78780 322378 0 253 0 0 60720 315380 0 254 0 0 60720 315380 0 255 0 0 60720 315380 0 256 0 0 60720 315380 0 257 0 0 60720 315380 0 258 10327 0 60720 315380 0 259 13718 91 88577 640273 135 260 3748 61 96448 345783 104 261 14416 9 50350 652925 53 262 1526 361 49857 439798 220 263 666 25 69351 278990 -4 264 2844 171 117869 339836 25 265 0 67 72683 240897 -15 266 368 0 60720 315380 0 267 0 17 61167 297141 0 268 333 0 60720 315380 0 269 26 13 70811 331323 0 270 0 1 60896 313880 0 271 0 0 60720 315380 0 272 1303 0 60720 315380 0 273 20 61 69863 309422 -5 274 2384 3 60938 315245 0 275 203 97 61348 405972 54 276 71 27 50804 300962 0 277 53 1 60745 316176 0 278 562 14 59506 302409 0 279 622 13 58456 283587 -3 280 645 7 60950 263276 -2 281 1763 5 60720 312075 0 282 0 24 61600 308336 0 283 317 0 60720 315380 0 284 1 9 63915 298700 0 285 275 6 60719 315372 0 286 0 4 59500 318745 0 287 936 0 60720 315380 0 288 8568 30 67939 408881 9 289 11528 37 32168 786690 223 290 0 70 -14545 -83265 -25 291 738 0 60720 315380 0 292 0 10 60720 321376 0 293 592 0 60720 315380 0 294 126 7 64270 276898 -4 295 2 10 60951 315547 0 296 0 1 60743 315487 0 297 18014 0 60720 315380 0 298 0 239 -1710 1405225 943 299 0 0 60720 315380 0 300 37704 0 60720 315380 0 301 63 452 60448 983660 2167 302 1431 4 65688 318574 0 303 94 34 106885 310768 30 304 192 3 61360 312887 0 305 32 4 65276 312339 0 306 6869 1 59988 314964 0 307 0 105 117520 379983 90 308 0 0 60720 315380 0 309 1 0 60720 315380 0 310 0 4 60722 315398 0 311 2328 0 60720 315380 0 312 209 94 82732 253588 0 313 28 13 64016 316647 0 314 176 2 60890 315688 0 315 1920 28 68136 310670 0 316 87550 15 79420 165404 -6 317 520 458 153198 4111912 17936 318 0 10 58650 291650 0 319 1013 0 60720 315380 0 320 15 2 69770 253468 -14 321 587 8 60831 315688 0 322 5371 12 59595 325699 0 323 0 66 87720 446211 49 324 1012 0 60720 315380 0 325 0 47 114768 368078 26 326 876 0 60720 315380 0 327 162556 59 138971 352850 24 328 43556 1081 213118 6282154 56289 329 3425 247 32648 227132 929 330 810 12 83620 283910 -3 331 0 43 74015 236761 26 332 0 0 60720 315380 0 333 2365 0 60720 315380 0 334 1261 79 191778 550608 120 335 0 33 76114 307528 13 336 1585 0 60720 315380 0 337 16189 107 93099 355864 2 338 0 295 116384 358589 379 339 0 0 60720 315380 0 340 10579 0 60720 315380 0 341 474 92 110309 375195 4 342 0 20 61977 288985 0 343 3642 0 60720 315380 0 344 0 86 90262 458343 22 345 472 0 60720 315380 0 346 98 20 80045 269587 0 347 3999 2 61490 315236 0 348 0 30 51252 42754 -58 349 621 0 60720 315380 0 350 0 48 60720 308256 0 351 0 0 60720 315380 0 352 30 0 60720 315380 0 353 0 4 60798 313164 0 354 746 0 60720 315380 0 355 0 23 71561 269661 10 356 0 0 60720 315380 0 357 0 0 60720 315380 0 358 0 0 60720 315380 0 359 4150 0 60720 315380 0 360 4658 150 134759 518365 66 361 814 158 156608 233773 58 362 1002 18 39625 301881 23 363 496 15 56750 298568 0 364 389 17 87390 325479 0 365 12679 14 58990 325506 0 366 400 110 48020 984885 1871 367 53 152 60720 313267 0 368 0 6 60349 315793 0 369 5109 0 60720 315380 0 370 576 97 67038 215362 -31 371 438 18 113761 314073 12 372 165 11 58320 298096 -5 373 4069 24 62841 325176 0 374 23 145 79804 207393 119 375 1285 5 62555 314806 0 376 4677 107 99489 341340 49 377 24811 75 90131 426280 24 378 167 505 95350 929118 822 379 0 2 64245 307322 0 380 4628 0 60720 315380 0 381 226 62 69159 387475 85 382 1765 5 65745 291787 0 383 460 20 77623 247060 28 384 36 15 63346 329784 0 385 989 3 60894 315834 0 386 3055 11 58930 304555 -5 387 0 43 60247 376641 14 388 0 0 60720 315380 0 389 0 0 60720 315380 0 390 0 0 60720 315380 0 391 0 0 60720 315380 0 392 13253 0 60720 315380 0 393 24 947 90829 357760 35 394 0 8 59818 315637 0 395 0 0 60720 315380 0 396 0 0 60720 315380 0 397 0 0 60720 315380 0 398 131 0 60720 315380 0 399 514 6 59661 327071 0 400 3366 15 102725 377516 23 401 9327 80 108479 299243 12 402 384 101 87419 387699 44 403 17821 22 54683 309836 0 404 397 373 122844 444477 643 405 897 10 62710 322327 -1 406 218 5 60720 314913 0 407 3369 21 60720 315553 0 408 5702 28 87161 688779 217 409 8636 61 101481 321896 10 410 534 855 128294 301607 108 411 0 7 62620 304485 0 412 0 0 60720 315380 0 413 726 0 60720 315380 0 414 1380 3 69980 231861 -5 415 180 26 60982 347385 25 416 7285 11 59635 316386 0 417 880 115 188873 491303 59 418 1 74 80791 261216 34 419 0 2 60727 315388 0 420 96 0 60720 315380 0 421 1889 5 60379 313729 0 422 45187 45 37527 358649 14 423 288 353 234817 1926517 2854 424 1270 4 60510 296656 0 425 6526 26 69206 275311 -10 426 0 26 55830 -42143 -28 427 226 0 60720 315380 0 428 694 12 72835 343929 0 429 0 8 68060 367655 -2 430 249 0 60720 315380 0 431 1595 72 56726 313491 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Costs Trades Dividends Wealth -1.245e+03 -6.645e-02 4.547e-02 2.878e-03 5.080e-03 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -17095 -2424 -1405 -514 156470 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -1.245e+03 9.850e+02 -1.264 0.20697 Costs -6.645e-02 1.263e-01 -0.526 0.59912 Trades 4.547e-02 1.587e-02 2.866 0.00436 ** Dividends 2.878e-03 1.542e-03 1.866 0.06268 . Wealth 5.080e-03 2.546e-03 1.995 0.04666 * --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 10590 on 426 degrees of freedom Multiple R-squared: 0.04961, Adjusted R-squared: 0.04069 F-statistic: 5.56 on 4 and 426 DF, p-value: 0.000227 > 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,] 6.355009e-09 1.271002e-08 1.000000e+00 [2,] 9.135221e-12 1.827044e-11 1.000000e+00 [3,] 1.449072e-14 2.898144e-14 1.000000e+00 [4,] 1.080923e-16 2.161846e-16 1.000000e+00 [5,] 1.974032e-19 3.948064e-19 1.000000e+00 [6,] 3.066813e-22 6.133625e-22 1.000000e+00 [7,] 9.776785e-24 1.955357e-23 1.000000e+00 [8,] 1.855129e-26 3.710259e-26 1.000000e+00 [9,] 4.133112e-29 8.266224e-29 1.000000e+00 [10,] 1.039992e-31 2.079983e-31 1.000000e+00 [11,] 2.587061e-34 5.174123e-34 1.000000e+00 [12,] 5.256491e-37 1.051298e-36 1.000000e+00 [13,] 9.422029e-40 1.884406e-39 1.000000e+00 [14,] 1.695782e-42 3.391564e-42 1.000000e+00 [15,] 7.203974e-45 1.440795e-44 1.000000e+00 [16,] 1.245824e-47 2.491649e-47 1.000000e+00 [17,] 1.039243e-49 2.078485e-49 1.000000e+00 [18,] 9.860932e-52 1.972186e-51 1.000000e+00 [19,] 1.780336e-54 3.560673e-54 1.000000e+00 [20,] 1.703993e-56 3.407987e-56 1.000000e+00 [21,] 6.968396e-58 1.393679e-57 1.000000e+00 [22,] 1.471697e-60 2.943393e-60 1.000000e+00 [23,] 4.126217e-63 8.252433e-63 1.000000e+00 [24,] 1.463468e-61 2.926935e-61 1.000000e+00 [25,] 4.246066e-64 8.492133e-64 1.000000e+00 [26,] 1.377975e-66 2.755950e-66 1.000000e+00 [27,] 4.181158e-69 8.362316e-69 1.000000e+00 [28,] 1.429607e-71 2.859213e-71 1.000000e+00 [29,] 7.942876e-73 1.588575e-72 1.000000e+00 [30,] 2.354922e-75 4.709845e-75 1.000000e+00 [31,] 1.818977e-77 3.637954e-77 1.000000e+00 [32,] 1.026964e-63 2.053929e-63 1.000000e+00 [33,] 2.208274e-65 4.416549e-65 1.000000e+00 [34,] 1.950868e-67 3.901736e-67 1.000000e+00 [35,] 3.239172e-69 6.478344e-69 1.000000e+00 [36,] 5.639839e-71 1.127968e-70 1.000000e+00 [37,] 5.491757e-73 1.098351e-72 1.000000e+00 [38,] 7.647239e-75 1.529448e-74 1.000000e+00 [39,] 7.058672e-77 1.411734e-76 1.000000e+00 [40,] 6.598836e-79 1.319767e-78 1.000000e+00 [41,] 6.262627e-81 1.252525e-80 1.000000e+00 [42,] 5.401203e-83 1.080241e-82 1.000000e+00 [43,] 4.617989e-85 9.235977e-85 1.000000e+00 [44,] 2.738523e-77 5.477046e-77 1.000000e+00 [45,] 4.018980e-79 8.037961e-79 1.000000e+00 [46,] 7.058372e-81 1.411674e-80 1.000000e+00 [47,] 1.374041e-81 2.748083e-81 1.000000e+00 [48,] 1.255381e-81 2.510762e-81 1.000000e+00 [49,] 1.036583e-79 2.073166e-79 1.000000e+00 [50,] 1.838095e-81 3.676190e-81 1.000000e+00 [51,] 4.103905e-83 8.207810e-83 1.000000e+00 [52,] 7.080039e-85 1.416008e-84 1.000000e+00 [53,] 1.225312e-86 2.450624e-86 1.000000e+00 [54,] 2.097743e-88 4.195487e-88 1.000000e+00 [55,] 3.521624e-90 7.043248e-90 1.000000e+00 [56,] 6.550513e-92 1.310103e-91 1.000000e+00 [57,] 1.086812e-93 2.173624e-93 1.000000e+00 [58,] 5.969403e-95 1.193881e-94 1.000000e+00 [59,] 9.883463e-97 1.976693e-96 1.000000e+00 [60,] 1.585118e-98 3.170236e-98 1.000000e+00 [61,] 2.530176e-100 5.060352e-100 1.000000e+00 [62,] 4.003172e-102 8.006344e-102 1.000000e+00 [63,] 6.393274e-104 1.278655e-103 1.000000e+00 [64,] 1.699581e-105 3.399162e-105 1.000000e+00 [65,] 2.655375e-107 5.310749e-107 1.000000e+00 [66,] 1.621016e-61 3.242033e-61 1.000000e+00 [67,] 1.239384e-62 2.478768e-62 1.000000e+00 [68,] 9.404972e-64 1.880994e-63 1.000000e+00 [69,] 7.211589e-65 1.442318e-64 1.000000e+00 [70,] 5.541259e-66 1.108252e-65 1.000000e+00 [71,] 4.098929e-67 8.197857e-67 1.000000e+00 [72,] 2.242646e-67 4.485291e-67 1.000000e+00 [73,] 1.898967e-68 3.797935e-68 1.000000e+00 [74,] 1.411180e-69 2.822360e-69 1.000000e+00 [75,] 1.040439e-70 2.080878e-70 1.000000e+00 [76,] 7.610784e-72 1.522157e-71 1.000000e+00 [77,] 5.539715e-73 1.107943e-72 1.000000e+00 [78,] 3.989917e-74 7.979833e-74 1.000000e+00 [79,] 2.851788e-75 5.703576e-75 1.000000e+00 [80,] 2.022319e-76 4.044638e-76 1.000000e+00 [81,] 1.424042e-77 2.848084e-77 1.000000e+00 [82,] 9.952282e-79 1.990456e-78 1.000000e+00 [83,] 6.908787e-80 1.381757e-79 1.000000e+00 [84,] 4.757988e-81 9.515975e-81 1.000000e+00 [85,] 3.252465e-82 6.504930e-82 1.000000e+00 [86,] 2.458891e-83 4.917783e-83 1.000000e+00 [87,] 5.774142e-84 1.154828e-83 1.000000e+00 [88,] 3.948514e-85 7.897029e-85 1.000000e+00 [89,] 2.678534e-86 5.357068e-86 1.000000e+00 [90,] 1.796914e-87 3.593829e-87 1.000000e+00 [91,] 1.197050e-88 2.394100e-88 1.000000e+00 [92,] 7.946404e-90 1.589281e-89 1.000000e+00 [93,] 5.212504e-91 1.042501e-90 1.000000e+00 [94,] 6.265803e-92 1.253161e-91 1.000000e+00 [95,] 4.097152e-93 8.194303e-93 1.000000e+00 [96,] 2.657281e-94 5.314563e-94 1.000000e+00 [97,] 1.710868e-95 3.421737e-95 1.000000e+00 [98,] 1.094237e-96 2.188473e-96 1.000000e+00 [99,] 6.466360e-95 1.293272e-94 1.000000e+00 [100,] 4.365693e-96 8.731386e-96 1.000000e+00 [101,] 1.431757e-96 2.863514e-96 1.000000e+00 [102,] 9.743408e-98 1.948682e-97 1.000000e+00 [103,] 6.621878e-99 1.324376e-98 1.000000e+00 [104,] 4.603028e-100 9.206056e-100 1.000000e+00 [105,] 5.735959e-101 1.147192e-100 1.000000e+00 [106,] 3.817379e-102 7.634759e-102 1.000000e+00 [107,] 2.521017e-103 5.042035e-103 1.000000e+00 [108,] 1.657272e-104 3.314544e-104 1.000000e+00 [109,] 1.083129e-105 2.166257e-105 1.000000e+00 [110,] 7.098872e-107 1.419774e-106 1.000000e+00 [111,] 7.726397e-107 1.545279e-106 1.000000e+00 [112,] 3.583620e-107 7.167240e-107 1.000000e+00 [113,] 4.060135e-108 8.120271e-108 1.000000e+00 [114,] 3.907457e-109 7.814914e-109 1.000000e+00 [115,] 2.587052e-110 5.174105e-110 1.000000e+00 [116,] 2.422664e-110 4.845329e-110 1.000000e+00 [117,] 1.626297e-111 3.252594e-111 1.000000e+00 [118,] 1.830911e-112 3.661822e-112 1.000000e+00 [119,] 1.216073e-113 2.432145e-113 1.000000e+00 [120,] 8.029455e-115 1.605891e-114 1.000000e+00 [121,] 5.272320e-116 1.054464e-115 1.000000e+00 [122,] 1.435193e-113 2.870386e-113 1.000000e+00 [123,] 1.025980e-114 2.051960e-114 1.000000e+00 [124,] 7.101983e-116 1.420397e-115 1.000000e+00 [125,] 4.879584e-117 9.759168e-117 1.000000e+00 [126,] 3.335184e-118 6.670369e-118 1.000000e+00 [127,] 2.269188e-119 4.538377e-119 1.000000e+00 [128,] 1.539602e-120 3.079203e-120 1.000000e+00 [129,] 3.164801e-121 6.329602e-121 1.000000e+00 [130,] 5.549521e-122 1.109904e-121 1.000000e+00 [131,] 4.191066e-123 8.382132e-123 1.000000e+00 [132,] 2.858683e-124 5.717366e-124 1.000000e+00 [133,] 2.757753e-121 5.515506e-121 1.000000e+00 [134,] 2.272874e-122 4.545748e-122 1.000000e+00 [135,] 2.256224e-123 4.512448e-123 1.000000e+00 [136,] 4.026316e-124 8.052632e-124 1.000000e+00 [137,] 3.316278e-125 6.632557e-125 1.000000e+00 [138,] 3.722789e-126 7.445579e-126 1.000000e+00 [139,] 2.882557e-127 5.765115e-127 1.000000e+00 [140,] 2.603447e-128 5.206894e-128 1.000000e+00 [141,] 3.845593e-129 7.691187e-129 1.000000e+00 [142,] 1.552153e-129 3.104307e-129 1.000000e+00 [143,] 1.807332e-130 3.614665e-130 1.000000e+00 [144,] 1.335573e-131 2.671146e-131 1.000000e+00 [145,] 1.443281e-132 2.886563e-132 1.000000e+00 [146,] 1.022703e-133 2.045406e-133 1.000000e+00 [147,] 7.177108e-135 1.435422e-134 1.000000e+00 [148,] 6.724100e-136 1.344820e-135 1.000000e+00 [149,] 2.133752e-136 4.267505e-136 1.000000e+00 [150,] 1.917652e-137 3.835305e-137 1.000000e+00 [151,] 1.648861e-138 3.297722e-138 1.000000e+00 [152,] 1.718317e-139 3.436634e-139 1.000000e+00 [153,] 5.827789e-139 1.165558e-138 1.000000e+00 [154,] 4.927695e-140 9.855390e-140 1.000000e+00 [155,] 3.601478e-141 7.202955e-141 1.000000e+00 [156,] 2.614367e-142 5.228733e-142 1.000000e+00 [157,] 2.462159e-143 4.924318e-143 1.000000e+00 [158,] 9.304644e-141 1.860929e-140 1.000000e+00 [159,] 2.826516e-141 5.653031e-141 1.000000e+00 [160,] 2.159855e-142 4.319710e-142 1.000000e+00 [161,] 2.867877e-139 5.735753e-139 1.000000e+00 [162,] 3.459159e-139 6.918318e-139 1.000000e+00 [163,] 5.380062e-140 1.076012e-139 1.000000e+00 [164,] 6.236011e-140 1.247202e-139 1.000000e+00 [165,] 4.926119e-141 9.852237e-141 1.000000e+00 [166,] 5.098702e-142 1.019740e-141 1.000000e+00 [167,] 6.894762e-143 1.378952e-142 1.000000e+00 [168,] 5.921565e-144 1.184313e-143 1.000000e+00 [169,] 1.251246e-144 2.502493e-144 1.000000e+00 [170,] 1.104826e-145 2.209652e-145 1.000000e+00 [171,] 8.776002e-147 1.755200e-146 1.000000e+00 [172,] 8.116032e-148 1.623206e-147 1.000000e+00 [173,] 1.812266e-145 3.624531e-145 1.000000e+00 [174,] 2.603590e-145 5.207180e-145 1.000000e+00 [175,] 1.668377e-145 3.336754e-145 1.000000e+00 [176,] 1.523634e-146 3.047268e-146 1.000000e+00 [177,] 1.382216e-147 2.764433e-147 1.000000e+00 [178,] 1.244788e-148 2.489575e-148 1.000000e+00 [179,] 1.099044e-149 2.198088e-149 1.000000e+00 [180,] 9.141012e-151 1.828202e-150 1.000000e+00 [181,] 1.320131e-22 2.640262e-22 1.000000e+00 [182,] 6.539563e-23 1.307913e-22 1.000000e+00 [183,] 3.268426e-23 6.536852e-23 1.000000e+00 [184,] 1.595432e-23 3.190864e-23 1.000000e+00 [185,] 7.905886e-24 1.581177e-23 1.000000e+00 [186,] 3.838458e-24 7.676916e-24 1.000000e+00 [187,] 1.829378e-24 3.658756e-24 1.000000e+00 [188,] 9.145385e-25 1.829077e-24 1.000000e+00 [189,] 4.269519e-25 8.539039e-25 1.000000e+00 [190,] 1.989307e-25 3.978614e-25 1.000000e+00 [191,] 9.358103e-26 1.871621e-25 1.000000e+00 [192,] 4.490088e-26 8.980177e-26 1.000000e+00 [193,] 1.446296e-24 2.892591e-24 1.000000e+00 [194,] 8.391581e-25 1.678316e-24 1.000000e+00 [195,] 4.113225e-25 8.226450e-25 1.000000e+00 [196,] 1.971736e-25 3.943472e-25 1.000000e+00 [197,] 9.562663e-26 1.912533e-25 1.000000e+00 [198,] 4.434294e-26 8.868587e-26 1.000000e+00 [199,] 2.046302e-26 4.092604e-26 1.000000e+00 [200,] 1.034614e-26 2.069227e-26 1.000000e+00 [201,] 4.867161e-27 9.734322e-27 1.000000e+00 [202,] 2.222400e-27 4.444801e-27 1.000000e+00 [203,] 6.585074e-26 1.317015e-25 1.000000e+00 [204,] 3.335289e-26 6.670577e-26 1.000000e+00 [205,] 1.645543e-26 3.291085e-26 1.000000e+00 [206,] 7.774521e-27 1.554904e-26 1.000000e+00 [207,] 1.309976e-26 2.619953e-26 1.000000e+00 [208,] 7.863099e-27 1.572620e-26 1.000000e+00 [209,] 3.861431e-27 7.722861e-27 1.000000e+00 [210,] 1.849486e-27 3.698973e-27 1.000000e+00 [211,] 8.471627e-28 1.694325e-27 1.000000e+00 [212,] 3.840444e-28 7.680889e-28 1.000000e+00 [213,] 2.109667e-28 4.219334e-28 1.000000e+00 [214,] 9.911567e-29 1.982313e-28 1.000000e+00 [215,] 4.509720e-29 9.019439e-29 1.000000e+00 [216,] 2.065923e-29 4.131845e-29 1.000000e+00 [217,] 9.429323e-30 1.885865e-29 1.000000e+00 [218,] 7.720242e-29 1.544048e-28 1.000000e+00 [219,] 3.970073e-29 7.940146e-29 1.000000e+00 [220,] 2.313574e-29 4.627148e-29 1.000000e+00 [221,] 2.555871e-29 5.111742e-29 1.000000e+00 [222,] 1.188713e-29 2.377427e-29 1.000000e+00 [223,] 5.302085e-30 1.060417e-29 1.000000e+00 [224,] 7.120146e-30 1.424029e-29 1.000000e+00 [225,] 9.244596e-30 1.848919e-29 1.000000e+00 [226,] 5.250926e-30 1.050185e-29 1.000000e+00 [227,] 2.411001e-30 4.822002e-30 1.000000e+00 [228,] 1.070399e-30 2.140798e-30 1.000000e+00 [229,] 4.749662e-31 9.499323e-31 1.000000e+00 [230,] 2.164159e-31 4.328318e-31 1.000000e+00 [231,] 9.811165e-32 1.962233e-31 1.000000e+00 [232,] 4.423177e-32 8.846354e-32 1.000000e+00 [233,] 1.949868e-32 3.899737e-32 1.000000e+00 [234,] 9.007438e-33 1.801488e-32 1.000000e+00 [235,] 4.148043e-33 8.296087e-33 1.000000e+00 [236,] 1.809456e-33 3.618911e-33 1.000000e+00 [237,] 8.526697e-34 1.705339e-33 1.000000e+00 [238,] 3.793147e-34 7.586295e-34 1.000000e+00 [239,] 1.697779e-34 3.395557e-34 1.000000e+00 [240,] 7.371153e-35 1.474231e-34 1.000000e+00 [241,] 3.195337e-35 6.390675e-35 1.000000e+00 [242,] 1.381074e-35 2.762148e-35 1.000000e+00 [243,] 5.943969e-36 1.188794e-35 1.000000e+00 [244,] 2.500815e-36 5.001629e-36 1.000000e+00 [245,] 1.096252e-36 2.192504e-36 1.000000e+00 [246,] 4.662055e-37 9.324110e-37 1.000000e+00 [247,] 1.974129e-37 3.948258e-37 1.000000e+00 [248,] 8.323521e-38 1.664704e-37 1.000000e+00 [249,] 3.494415e-38 6.988830e-38 1.000000e+00 [250,] 1.460764e-38 2.921528e-38 1.000000e+00 [251,] 1.366975e-38 2.733951e-38 1.000000e+00 [252,] 1.365018e-38 2.730037e-38 1.000000e+00 [253,] 5.697156e-39 1.139431e-38 1.000000e+00 [254,] 6.882827e-39 1.376565e-38 1.000000e+00 [255,] 2.923155e-39 5.846311e-39 1.000000e+00 [256,] 1.187891e-39 2.375782e-39 1.000000e+00 [257,] 5.025793e-40 1.005159e-39 1.000000e+00 [258,] 2.066367e-40 4.132735e-40 1.000000e+00 [259,] 8.337316e-41 1.667463e-40 1.000000e+00 [260,] 3.379082e-41 6.758164e-41 1.000000e+00 [261,] 1.353180e-41 2.706361e-41 1.000000e+00 [262,] 5.582541e-42 1.116508e-41 1.000000e+00 [263,] 2.249545e-42 4.499091e-42 1.000000e+00 [264,] 9.028984e-43 1.805797e-42 1.000000e+00 [265,] 3.448727e-43 6.897454e-43 1.000000e+00 [266,] 1.374539e-43 2.749077e-43 1.000000e+00 [267,] 5.174104e-44 1.034821e-43 1.000000e+00 [268,] 2.095197e-44 4.190393e-44 1.000000e+00 [269,] 8.023418e-45 1.604684e-44 1.000000e+00 [270,] 3.125512e-45 6.251024e-45 1.000000e+00 [271,] 1.178087e-45 2.356174e-45 1.000000e+00 [272,] 4.397364e-46 8.794728e-46 1.000000e+00 [273,] 1.641481e-46 3.282962e-46 1.000000e+00 [274,] 5.972512e-47 1.194502e-46 1.000000e+00 [275,] 2.269851e-47 4.539702e-47 1.000000e+00 [276,] 8.497050e-48 1.699410e-47 1.000000e+00 [277,] 3.215841e-48 6.431682e-48 1.000000e+00 [278,] 1.194283e-48 2.388566e-48 1.000000e+00 [279,] 4.474488e-49 8.948977e-49 1.000000e+00 [280,] 1.607086e-49 3.214172e-49 1.000000e+00 [281,] 8.523482e-50 1.704696e-49 1.000000e+00 [282,] 5.994889e-50 1.198978e-49 1.000000e+00 [283,] 2.380717e-50 4.761434e-50 1.000000e+00 [284,] 8.464673e-51 1.692935e-50 1.000000e+00 [285,] 3.099942e-51 6.199884e-51 1.000000e+00 [286,] 1.097849e-51 2.195698e-51 1.000000e+00 [287,] 3.949921e-52 7.899843e-52 1.000000e+00 [288,] 1.424379e-52 2.848757e-52 1.000000e+00 [289,] 5.121351e-53 1.024270e-52 1.000000e+00 [290,] 5.812597e-52 1.162519e-51 1.000000e+00 [291,] 3.611454e-51 7.222909e-51 1.000000e+00 [292,] 1.303537e-51 2.607073e-51 1.000000e+00 [293,] 3.882168e-45 7.764337e-45 1.000000e+00 [294,] 2.444750e-45 4.889500e-45 1.000000e+00 [295,] 9.029282e-46 1.805856e-45 1.000000e+00 [296,] 4.445364e-46 8.890727e-46 1.000000e+00 [297,] 1.674072e-46 3.348143e-46 1.000000e+00 [298,] 6.403981e-47 1.280796e-46 1.000000e+00 [299,] 3.047112e-47 6.094225e-47 1.000000e+00 [300,] 1.549263e-47 3.098526e-47 1.000000e+00 [301,] 5.773357e-48 1.154671e-47 1.000000e+00 [302,] 2.140656e-48 4.281312e-48 1.000000e+00 [303,] 7.890016e-49 1.578003e-48 1.000000e+00 [304,] 2.759074e-49 5.518148e-49 1.000000e+00 [305,] 1.082943e-49 2.165887e-49 1.000000e+00 [306,] 3.963922e-50 7.927844e-50 1.000000e+00 [307,] 1.421729e-50 2.843457e-50 1.000000e+00 [308,] 4.907726e-51 9.815453e-51 1.000000e+00 [309,] 9.296435e-24 1.859287e-23 1.000000e+00 [310,] 5.790618e-24 1.158124e-23 1.000000e+00 [311,] 2.849208e-24 5.698415e-24 1.000000e+00 [312,] 1.366308e-24 2.732615e-24 1.000000e+00 [313,] 6.972594e-25 1.394519e-24 1.000000e+00 [314,] 3.330515e-25 6.661031e-25 1.000000e+00 [315,] 1.670830e-25 3.341661e-25 1.000000e+00 [316,] 8.429982e-26 1.685996e-25 1.000000e+00 [317,] 3.927732e-26 7.855463e-26 1.000000e+00 [318,] 2.456025e-26 4.912050e-26 1.000000e+00 [319,] 1.135525e-26 2.271050e-26 1.000000e+00 [320,] 1.000000e+00 1.686519e-12 8.432593e-13 [321,] 1.000000e+00 7.470124e-13 3.735062e-13 [322,] 1.000000e+00 9.618237e-13 4.809118e-13 [323,] 1.000000e+00 1.856430e-12 9.282149e-13 [324,] 1.000000e+00 3.547471e-12 1.773736e-12 [325,] 1.000000e+00 6.329640e-12 3.164820e-12 [326,] 1.000000e+00 1.220358e-11 6.101791e-12 [327,] 1.000000e+00 2.149302e-11 1.074651e-11 [328,] 1.000000e+00 3.863128e-11 1.931564e-11 [329,] 1.000000e+00 7.207940e-11 3.603970e-11 [330,] 1.000000e+00 2.378167e-11 1.189084e-11 [331,] 1.000000e+00 4.443872e-11 2.221936e-11 [332,] 1.000000e+00 7.841180e-11 3.920590e-11 [333,] 1.000000e+00 8.300183e-11 4.150091e-11 [334,] 1.000000e+00 1.505770e-10 7.528850e-11 [335,] 1.000000e+00 2.696785e-10 1.348393e-10 [336,] 1.000000e+00 5.057881e-10 2.528940e-10 [337,] 1.000000e+00 8.080858e-10 4.040429e-10 [338,] 1.000000e+00 1.427065e-09 7.135325e-10 [339,] 1.000000e+00 2.562581e-09 1.281291e-09 [340,] 1.000000e+00 4.665974e-09 2.332987e-09 [341,] 1.000000e+00 8.498698e-09 4.249349e-09 [342,] 1.000000e+00 1.471871e-08 7.359356e-09 [343,] 1.000000e+00 2.472924e-08 1.236462e-08 [344,] 1.000000e+00 4.104065e-08 2.052032e-08 [345,] 1.000000e+00 6.771729e-08 3.385864e-08 [346,] 9.999999e-01 1.109401e-07 5.547004e-08 [347,] 9.999999e-01 1.861573e-07 9.307863e-08 [348,] 9.999998e-01 3.126103e-07 1.563051e-07 [349,] 9.999997e-01 5.001518e-07 2.500759e-07 [350,] 9.999996e-01 7.941141e-07 3.970570e-07 [351,] 9.999994e-01 1.251163e-06 6.255816e-07 [352,] 9.999989e-01 2.113199e-06 1.056600e-06 [353,] 9.999982e-01 3.526443e-06 1.763222e-06 [354,] 9.999972e-01 5.561596e-06 2.780798e-06 [355,] 9.999956e-01 8.784262e-06 4.392131e-06 [356,] 9.999931e-01 1.372988e-05 6.864942e-06 [357,] 9.999892e-01 2.161115e-05 1.080557e-05 [358,] 9.999898e-01 2.047832e-05 1.023916e-05 [359,] 9.999836e-01 3.281802e-05 1.640901e-05 [360,] 9.999751e-01 4.970085e-05 2.485043e-05 [361,] 9.999629e-01 7.428604e-05 3.714302e-05 [362,] 9.999420e-01 1.160242e-04 5.801209e-05 [363,] 9.999094e-01 1.811884e-04 9.059422e-05 [364,] 9.998607e-01 2.785229e-04 1.392614e-04 [365,] 9.997954e-01 4.092862e-04 2.046431e-04 [366,] 9.996864e-01 6.272569e-04 3.136285e-04 [367,] 9.995247e-01 9.505048e-04 4.752524e-04 [368,] 9.993004e-01 1.399169e-03 6.995846e-04 [369,] 9.989851e-01 2.029748e-03 1.014874e-03 [370,] 9.998718e-01 2.564812e-04 1.282406e-04 [371,] 9.998283e-01 3.434906e-04 1.717453e-04 [372,] 9.997414e-01 5.171018e-04 2.585509e-04 [373,] 9.995979e-01 8.041495e-04 4.020747e-04 [374,] 9.994111e-01 1.177769e-03 5.888844e-04 [375,] 9.990972e-01 1.805523e-03 9.027615e-04 [376,] 9.986319e-01 2.736240e-03 1.368120e-03 [377,] 9.980594e-01 3.881248e-03 1.940624e-03 [378,] 9.971915e-01 5.617053e-03 2.808526e-03 [379,] 9.958603e-01 8.279328e-03 4.139664e-03 [380,] 9.944932e-01 1.101366e-02 5.506832e-03 [381,] 9.924994e-01 1.500127e-02 7.500634e-03 [382,] 9.898955e-01 2.020898e-02 1.010449e-02 [383,] 9.865383e-01 2.692345e-02 1.346172e-02 [384,] 9.822660e-01 3.546797e-02 1.773398e-02 [385,] 9.822820e-01 3.543609e-02 1.771805e-02 [386,] 9.787556e-01 4.248875e-02 2.124437e-02 [387,] 9.721992e-01 5.560150e-02 2.780075e-02 [388,] 9.638944e-01 7.221111e-02 3.610555e-02 [389,] 9.536783e-01 9.264345e-02 4.632173e-02 [390,] 9.412996e-01 1.174009e-01 5.870043e-02 [391,] 9.262661e-01 1.474677e-01 7.373386e-02 [392,] 9.085606e-01 1.828788e-01 9.143940e-02 [393,] 8.817873e-01 2.364254e-01 1.182127e-01 [394,] 8.809948e-01 2.380103e-01 1.190052e-01 [395,] 8.511728e-01 2.976543e-01 1.488272e-01 [396,] 8.754873e-01 2.490254e-01 1.245127e-01 [397,] 8.413464e-01 3.173071e-01 1.586536e-01 [398,] 8.029967e-01 3.940065e-01 1.970033e-01 [399,] 7.619043e-01 4.761915e-01 2.380957e-01 [400,] 7.064266e-01 5.871468e-01 2.935734e-01 [401,] 6.542459e-01 6.915081e-01 3.457541e-01 [402,] 6.338070e-01 7.323861e-01 3.661930e-01 [403,] 8.200569e-01 3.598863e-01 1.799431e-01 [404,] 7.666385e-01 4.667230e-01 2.333615e-01 [405,] 7.027285e-01 5.945430e-01 2.972715e-01 [406,] 6.285612e-01 7.428777e-01 3.714388e-01 [407,] 5.651762e-01 8.696475e-01 4.348238e-01 [408,] 5.252115e-01 9.495770e-01 4.747885e-01 [409,] 4.402295e-01 8.804589e-01 5.597705e-01 [410,] 7.995580e-01 4.008840e-01 2.004420e-01 [411,] 7.717690e-01 4.564620e-01 2.282310e-01 [412,] 6.967081e-01 6.065838e-01 3.032919e-01 [413,] 6.112946e-01 7.774108e-01 3.887054e-01 [414,] 5.061653e-01 9.876693e-01 4.938347e-01 [415,] 9.951747e-01 9.650621e-03 4.825311e-03 [416,] 9.811266e-01 3.774679e-02 1.887339e-02 > postscript(file="/var/www/rcomp/tmp/1uinc1292675864.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/rcomp/tmp/2nr4f1292675864.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/rcomp/tmp/3nr4f1292675864.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/rcomp/tmp/4nr4f1292675864.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/rcomp/tmp/5gi401292675864.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 -446.799014 -1126.182058 -511.543395 -902.728119 120.760443 6 7 8 9 10 -344.107520 -381.130939 -524.381751 -660.956188 -460.062309 11 12 13 14 15 -894.317669 -904.534116 -531.752996 -373.326204 -494.215890 16 17 18 19 20 -369.054268 -522.363459 -541.639199 -543.244901 -522.622415 21 22 23 24 25 -531.752996 -385.786972 -531.752996 -1243.621731 -743.010377 26 27 28 29 30 -515.477844 -408.176768 -784.862499 -531.752996 -659.261369 31 32 33 34 35 91.977173 -531.752996 -590.863508 -1726.837388 -307.246345 36 37 38 39 40 -140.114489 -529.762446 -180.711314 -668.761908 -835.720751 41 42 43 44 45 -1228.182183 -781.072067 -914.536718 -531.752996 -449.831058 46 47 48 49 50 -436.851090 -531.752996 -492.797868 -646.497911 -957.033521 51 52 53 54 55 -240.473167 -701.717963 -219.629339 -581.289588 -3383.558440 56 57 58 59 60 -766.028236 -511.069047 -97.596255 -357.909723 -531.752996 61 62 63 64 65 -531.752996 -526.747763 145.084534 -534.244498 655.246817 66 67 68 69 70 -526.236610 -856.532957 -602.958505 -457.130721 -531.752996 71 72 73 74 75 -3174.838918 211.927550 9354.448906 -531.752996 -578.746650 76 77 78 79 80 -581.164370 -633.399028 -518.399530 1113.396719 -37.206470 81 82 83 84 85 -529.007476 -531.752996 -531.752996 -554.257189 -531.752996 86 87 88 89 90 -531.752996 -521.732620 -531.449203 -444.201528 -670.889173 91 92 93 94 95 -531.752996 -401.698129 -1291.243244 465.050125 -563.965545 96 97 98 99 100 -630.510274 -531.752996 -514.103829 -542.910900 -264.378147 101 102 103 104 105 -1074.155039 -434.808579 -542.900405 -531.755114 -531.752996 106 107 108 109 110 2029.650731 -431.446382 -1695.157848 -384.687406 -510.562047 111 112 113 114 115 -518.842890 51.592784 -531.752996 -525.641602 -531.752996 116 117 118 119 120 -461.320462 -888.829048 -1887.590163 -2138.572831 -1490.581032 121 122 123 124 125 3.856788 -531.752996 -2095.853648 -878.399845 -20.323604 126 127 128 129 130 -531.752996 -545.554349 -531.752996 1523.730439 -476.086445 131 132 133 134 135 -585.320451 -531.752996 -531.752996 -531.752996 -520.235254 136 137 138 139 140 -472.722529 -2088.524181 -2592.868414 -2423.524181 2101.475819 141 142 143 144 145 -1404.972548 -575.483408 -2423.524181 -1785.524181 -2457.990409 146 147 148 149 150 -1816.524181 -871.909281 -1019.813082 -2653.563065 -2423.524181 151 152 153 154 155 -1812.524181 -2300.277558 -1500.524181 -1648.256178 -64.600804 156 157 158 159 160 -2586.926486 -2553.057962 -2288.524181 -822.180481 -4711.001339 161 162 163 164 165 -2412.082568 -2114.524181 -2176.304013 -2596.137187 2095.475819 166 167 168 169 170 -5116.656277 -1586.524181 2674.509898 -336.323074 -240.871400 171 172 173 174 175 -3787.369387 -2440.813076 -1905.689883 -2410.750338 -1495.011601 176 177 178 179 180 -3495.484878 -2362.524181 -2052.208190 -2431.352844 2477.475819 181 182 183 184 185 -5103.395872 -4500.593947 -2423.524181 -2423.524181 -2421.524181 186 187 188 189 190 -2387.648451 -1677.093216 82103.337225 1803.203792 -2420.524181 191 192 193 194 195 -1896.916885 -2511.118916 -2018.524181 -1517.305540 -3044.373271 196 197 198 199 200 -324.524181 -1049.596659 -1391.588625 -2423.524181 19906.475819 201 202 203 204 205 -4429.102815 -2423.524181 -1940.524181 -2400.509309 -184.524181 206 207 208 209 210 -1107.305949 -3611.207824 -2058.524181 -1250.547881 19757.475870 211 212 213 214 215 -3095.283585 -2685.572436 -1892.462246 11022.287859 -6015.116649 216 217 218 219 220 -2776.875487 -2468.180819 -1009.524181 -909.398752 -5250.609171 221 222 223 224 225 -2423.625909 -1579.524181 -2068.290876 -2009.807846 16285.157028 226 227 228 229 230 -4485.171685 4699.475819 -13070.996208 -2771.687434 -366.986238 231 232 233 234 235 -14203.913927 9661.176538 -5228.101162 -2298.524181 -1229.369359 236 237 238 239 240 -1406.605618 -2303.724111 -2423.524181 -2393.524181 818.673002 241 242 243 244 245 -2909.410634 -3296.891128 -1955.524181 2504.134858 1401.321731 246 247 248 249 250 -1525.114460 -2298.524181 -2384.498284 -2401.524181 -2415.199970 251 252 253 254 255 -1909.524181 -3261.382924 -2423.524181 -2423.524181 -2423.524181 256 257 258 259 260 -2423.524181 -2423.524181 7903.475819 9098.228752 -384.019079 261 262 263 264 265 11492.941439 -738.767412 -2043.565720 -2237.191848 -2748.599377 266 267 268 269 270 -2055.524181 -2390.231484 -2090.524181 -2901.371416 -2427.143658 271 272 273 274 275 -2423.524181 -1120.524181 -2798.025251 -48.848634 -2503.611189 276 277 278 279 280 -1858.364649 -2373.885174 -1768.066474 -1606.209373 -1638.563020 281 282 283 284 285 -650.680867 -2441.671267 -2106.524181 -2519.199175 -2148.056984 286 287 288 289 290 -2377.469561 -1487.524181 5549.105588 9047.715220 2150.713080 291 292 293 294 295 -1685.524181 -2440.114833 -1831.524181 -2347.712027 -2431.843668 296 297 298 299 300 -2424.811446 15590.475819 -2710.117007 -2423.524181 35280.475819 301 302 303 304 305 -4252.290035 -1227.341467 -4413.234993 -2253.250875 -2589.665167 306 307 308 309 310 4480.022951 -5185.573175 -2423.524181 -2422.524181 -2423.401116 311 312 313 314 315 -95.524181 -3031.324247 -2548.173266 -2256.007415 -825.309520 316 317 318 319 320 84708.824961 -17094.684688 -2260.448538 -1410.524181 -2641.648261 321 322 323 324 325 -1841.926021 2969.730390 -4023.560995 -1411.524181 -5029.710777 326 327 328 329 330 -1547.524181 156470.426224 16817.902759 2543.542624 -2563.394724 331 332 333 334 335 -2799.065069 -2423.524181 -58.524181 -7793.933302 -3098.755366 336 337 338 339 340 -838.524181 12183.821560 -5061.195168 -2423.524181 8155.475819 341 342 343 344 345 -4370.342296 -2403.391098 1218.475819 -4172.590994 -1951.524181 346 347 348 349 350 -3071.107195 1541.011763 -1206.175433 -1802.524181 -2399.833242 351 352 353 354 355 -2423.524181 -2393.524181 -2420.427824 -1677.524181 -2783.410082 356 357 358 359 360 -2423.524181 -2423.524181 -2423.524181 1726.475819 -1706.538957 361 362 363 364 365 -5724.432563 -422.422949 -1697.633094 -3275.123340 10305.927627 366 367 368 369 370 -3374.945564 -2354.342898 -2407.444889 2685.475819 -1840.366319 371 372 373 374 375 -4392.364717 -2098.901984 1522.439542 -2948.466625 -1219.976240 376 377 378 379 380 422.830876 20735.894586 -5567.945758 -2560.481408 2204.475819 381 382 383 384 385 -2785.024837 -818.779707 -2534.295064 -2547.381315 -1443.542997 386 387 388 389 390 744.774172 -2575.526532 -2423.524181 -2423.524181 -2423.524181 391 392 393 394 395 -2423.524181 10829.475819 -3827.768280 -2382.718857 -2423.524181 396 397 398 399 400 -2423.524181 -2423.524181 -2292.524181 -1894.617602 -1145.395706 401 402 403 404 405 4783.601521 -3455.138968 15689.390287 -5201.251309 -1636.330421 406 407 408 409 410 -2203.848003 946.373432 1002.420813 4344.352928 -4866.163850 411 412 413 414 415 -2478.097288 -2423.524181 -1697.524181 -1223.995912 -2345.939710 416 417 418 419 420 4908.645920 -7869.476244 -3174.521760 -2423.732586 -2327.524181 421 422 423 424 425 -513.935683 43696.445781 -14679.127902 -1089.826295 3833.711605 426 427 428 429 430 -1170.437279 -2197.524181 -2361.745002 -2907.162935 -2174.524181 431 -636.699222 > postscript(file="/var/www/rcomp/tmp/6gi401292675864.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 -446.799014 NA 1 -1126.182058 -446.799014 2 -511.543395 -1126.182058 3 -902.728119 -511.543395 4 120.760443 -902.728119 5 -344.107520 120.760443 6 -381.130939 -344.107520 7 -524.381751 -381.130939 8 -660.956188 -524.381751 9 -460.062309 -660.956188 10 -894.317669 -460.062309 11 -904.534116 -894.317669 12 -531.752996 -904.534116 13 -373.326204 -531.752996 14 -494.215890 -373.326204 15 -369.054268 -494.215890 16 -522.363459 -369.054268 17 -541.639199 -522.363459 18 -543.244901 -541.639199 19 -522.622415 -543.244901 20 -531.752996 -522.622415 21 -385.786972 -531.752996 22 -531.752996 -385.786972 23 -1243.621731 -531.752996 24 -743.010377 -1243.621731 25 -515.477844 -743.010377 26 -408.176768 -515.477844 27 -784.862499 -408.176768 28 -531.752996 -784.862499 29 -659.261369 -531.752996 30 91.977173 -659.261369 31 -531.752996 91.977173 32 -590.863508 -531.752996 33 -1726.837388 -590.863508 34 -307.246345 -1726.837388 35 -140.114489 -307.246345 36 -529.762446 -140.114489 37 -180.711314 -529.762446 38 -668.761908 -180.711314 39 -835.720751 -668.761908 40 -1228.182183 -835.720751 41 -781.072067 -1228.182183 42 -914.536718 -781.072067 43 -531.752996 -914.536718 44 -449.831058 -531.752996 45 -436.851090 -449.831058 46 -531.752996 -436.851090 47 -492.797868 -531.752996 48 -646.497911 -492.797868 49 -957.033521 -646.497911 50 -240.473167 -957.033521 51 -701.717963 -240.473167 52 -219.629339 -701.717963 53 -581.289588 -219.629339 54 -3383.558440 -581.289588 55 -766.028236 -3383.558440 56 -511.069047 -766.028236 57 -97.596255 -511.069047 58 -357.909723 -97.596255 59 -531.752996 -357.909723 60 -531.752996 -531.752996 61 -526.747763 -531.752996 62 145.084534 -526.747763 63 -534.244498 145.084534 64 655.246817 -534.244498 65 -526.236610 655.246817 66 -856.532957 -526.236610 67 -602.958505 -856.532957 68 -457.130721 -602.958505 69 -531.752996 -457.130721 70 -3174.838918 -531.752996 71 211.927550 -3174.838918 72 9354.448906 211.927550 73 -531.752996 9354.448906 74 -578.746650 -531.752996 75 -581.164370 -578.746650 76 -633.399028 -581.164370 77 -518.399530 -633.399028 78 1113.396719 -518.399530 79 -37.206470 1113.396719 80 -529.007476 -37.206470 81 -531.752996 -529.007476 82 -531.752996 -531.752996 83 -554.257189 -531.752996 84 -531.752996 -554.257189 85 -531.752996 -531.752996 86 -521.732620 -531.752996 87 -531.449203 -521.732620 88 -444.201528 -531.449203 89 -670.889173 -444.201528 90 -531.752996 -670.889173 91 -401.698129 -531.752996 92 -1291.243244 -401.698129 93 465.050125 -1291.243244 94 -563.965545 465.050125 95 -630.510274 -563.965545 96 -531.752996 -630.510274 97 -514.103829 -531.752996 98 -542.910900 -514.103829 99 -264.378147 -542.910900 100 -1074.155039 -264.378147 101 -434.808579 -1074.155039 102 -542.900405 -434.808579 103 -531.755114 -542.900405 104 -531.752996 -531.755114 105 2029.650731 -531.752996 106 -431.446382 2029.650731 107 -1695.157848 -431.446382 108 -384.687406 -1695.157848 109 -510.562047 -384.687406 110 -518.842890 -510.562047 111 51.592784 -518.842890 112 -531.752996 51.592784 113 -525.641602 -531.752996 114 -531.752996 -525.641602 115 -461.320462 -531.752996 116 -888.829048 -461.320462 117 -1887.590163 -888.829048 118 -2138.572831 -1887.590163 119 -1490.581032 -2138.572831 120 3.856788 -1490.581032 121 -531.752996 3.856788 122 -2095.853648 -531.752996 123 -878.399845 -2095.853648 124 -20.323604 -878.399845 125 -531.752996 -20.323604 126 -545.554349 -531.752996 127 -531.752996 -545.554349 128 1523.730439 -531.752996 129 -476.086445 1523.730439 130 -585.320451 -476.086445 131 -531.752996 -585.320451 132 -531.752996 -531.752996 133 -531.752996 -531.752996 134 -520.235254 -531.752996 135 -472.722529 -520.235254 136 -2088.524181 -472.722529 137 -2592.868414 -2088.524181 138 -2423.524181 -2592.868414 139 2101.475819 -2423.524181 140 -1404.972548 2101.475819 141 -575.483408 -1404.972548 142 -2423.524181 -575.483408 143 -1785.524181 -2423.524181 144 -2457.990409 -1785.524181 145 -1816.524181 -2457.990409 146 -871.909281 -1816.524181 147 -1019.813082 -871.909281 148 -2653.563065 -1019.813082 149 -2423.524181 -2653.563065 150 -1812.524181 -2423.524181 151 -2300.277558 -1812.524181 152 -1500.524181 -2300.277558 153 -1648.256178 -1500.524181 154 -64.600804 -1648.256178 155 -2586.926486 -64.600804 156 -2553.057962 -2586.926486 157 -2288.524181 -2553.057962 158 -822.180481 -2288.524181 159 -4711.001339 -822.180481 160 -2412.082568 -4711.001339 161 -2114.524181 -2412.082568 162 -2176.304013 -2114.524181 163 -2596.137187 -2176.304013 164 2095.475819 -2596.137187 165 -5116.656277 2095.475819 166 -1586.524181 -5116.656277 167 2674.509898 -1586.524181 168 -336.323074 2674.509898 169 -240.871400 -336.323074 170 -3787.369387 -240.871400 171 -2440.813076 -3787.369387 172 -1905.689883 -2440.813076 173 -2410.750338 -1905.689883 174 -1495.011601 -2410.750338 175 -3495.484878 -1495.011601 176 -2362.524181 -3495.484878 177 -2052.208190 -2362.524181 178 -2431.352844 -2052.208190 179 2477.475819 -2431.352844 180 -5103.395872 2477.475819 181 -4500.593947 -5103.395872 182 -2423.524181 -4500.593947 183 -2423.524181 -2423.524181 184 -2421.524181 -2423.524181 185 -2387.648451 -2421.524181 186 -1677.093216 -2387.648451 187 82103.337225 -1677.093216 188 1803.203792 82103.337225 189 -2420.524181 1803.203792 190 -1896.916885 -2420.524181 191 -2511.118916 -1896.916885 192 -2018.524181 -2511.118916 193 -1517.305540 -2018.524181 194 -3044.373271 -1517.305540 195 -324.524181 -3044.373271 196 -1049.596659 -324.524181 197 -1391.588625 -1049.596659 198 -2423.524181 -1391.588625 199 19906.475819 -2423.524181 200 -4429.102815 19906.475819 201 -2423.524181 -4429.102815 202 -1940.524181 -2423.524181 203 -2400.509309 -1940.524181 204 -184.524181 -2400.509309 205 -1107.305949 -184.524181 206 -3611.207824 -1107.305949 207 -2058.524181 -3611.207824 208 -1250.547881 -2058.524181 209 19757.475870 -1250.547881 210 -3095.283585 19757.475870 211 -2685.572436 -3095.283585 212 -1892.462246 -2685.572436 213 11022.287859 -1892.462246 214 -6015.116649 11022.287859 215 -2776.875487 -6015.116649 216 -2468.180819 -2776.875487 217 -1009.524181 -2468.180819 218 -909.398752 -1009.524181 219 -5250.609171 -909.398752 220 -2423.625909 -5250.609171 221 -1579.524181 -2423.625909 222 -2068.290876 -1579.524181 223 -2009.807846 -2068.290876 224 16285.157028 -2009.807846 225 -4485.171685 16285.157028 226 4699.475819 -4485.171685 227 -13070.996208 4699.475819 228 -2771.687434 -13070.996208 229 -366.986238 -2771.687434 230 -14203.913927 -366.986238 231 9661.176538 -14203.913927 232 -5228.101162 9661.176538 233 -2298.524181 -5228.101162 234 -1229.369359 -2298.524181 235 -1406.605618 -1229.369359 236 -2303.724111 -1406.605618 237 -2423.524181 -2303.724111 238 -2393.524181 -2423.524181 239 818.673002 -2393.524181 240 -2909.410634 818.673002 241 -3296.891128 -2909.410634 242 -1955.524181 -3296.891128 243 2504.134858 -1955.524181 244 1401.321731 2504.134858 245 -1525.114460 1401.321731 246 -2298.524181 -1525.114460 247 -2384.498284 -2298.524181 248 -2401.524181 -2384.498284 249 -2415.199970 -2401.524181 250 -1909.524181 -2415.199970 251 -3261.382924 -1909.524181 252 -2423.524181 -3261.382924 253 -2423.524181 -2423.524181 254 -2423.524181 -2423.524181 255 -2423.524181 -2423.524181 256 -2423.524181 -2423.524181 257 7903.475819 -2423.524181 258 9098.228752 7903.475819 259 -384.019079 9098.228752 260 11492.941439 -384.019079 261 -738.767412 11492.941439 262 -2043.565720 -738.767412 263 -2237.191848 -2043.565720 264 -2748.599377 -2237.191848 265 -2055.524181 -2748.599377 266 -2390.231484 -2055.524181 267 -2090.524181 -2390.231484 268 -2901.371416 -2090.524181 269 -2427.143658 -2901.371416 270 -2423.524181 -2427.143658 271 -1120.524181 -2423.524181 272 -2798.025251 -1120.524181 273 -48.848634 -2798.025251 274 -2503.611189 -48.848634 275 -1858.364649 -2503.611189 276 -2373.885174 -1858.364649 277 -1768.066474 -2373.885174 278 -1606.209373 -1768.066474 279 -1638.563020 -1606.209373 280 -650.680867 -1638.563020 281 -2441.671267 -650.680867 282 -2106.524181 -2441.671267 283 -2519.199175 -2106.524181 284 -2148.056984 -2519.199175 285 -2377.469561 -2148.056984 286 -1487.524181 -2377.469561 287 5549.105588 -1487.524181 288 9047.715220 5549.105588 289 2150.713080 9047.715220 290 -1685.524181 2150.713080 291 -2440.114833 -1685.524181 292 -1831.524181 -2440.114833 293 -2347.712027 -1831.524181 294 -2431.843668 -2347.712027 295 -2424.811446 -2431.843668 296 15590.475819 -2424.811446 297 -2710.117007 15590.475819 298 -2423.524181 -2710.117007 299 35280.475819 -2423.524181 300 -4252.290035 35280.475819 301 -1227.341467 -4252.290035 302 -4413.234993 -1227.341467 303 -2253.250875 -4413.234993 304 -2589.665167 -2253.250875 305 4480.022951 -2589.665167 306 -5185.573175 4480.022951 307 -2423.524181 -5185.573175 308 -2422.524181 -2423.524181 309 -2423.401116 -2422.524181 310 -95.524181 -2423.401116 311 -3031.324247 -95.524181 312 -2548.173266 -3031.324247 313 -2256.007415 -2548.173266 314 -825.309520 -2256.007415 315 84708.824961 -825.309520 316 -17094.684688 84708.824961 317 -2260.448538 -17094.684688 318 -1410.524181 -2260.448538 319 -2641.648261 -1410.524181 320 -1841.926021 -2641.648261 321 2969.730390 -1841.926021 322 -4023.560995 2969.730390 323 -1411.524181 -4023.560995 324 -5029.710777 -1411.524181 325 -1547.524181 -5029.710777 326 156470.426224 -1547.524181 327 16817.902759 156470.426224 328 2543.542624 16817.902759 329 -2563.394724 2543.542624 330 -2799.065069 -2563.394724 331 -2423.524181 -2799.065069 332 -58.524181 -2423.524181 333 -7793.933302 -58.524181 334 -3098.755366 -7793.933302 335 -838.524181 -3098.755366 336 12183.821560 -838.524181 337 -5061.195168 12183.821560 338 -2423.524181 -5061.195168 339 8155.475819 -2423.524181 340 -4370.342296 8155.475819 341 -2403.391098 -4370.342296 342 1218.475819 -2403.391098 343 -4172.590994 1218.475819 344 -1951.524181 -4172.590994 345 -3071.107195 -1951.524181 346 1541.011763 -3071.107195 347 -1206.175433 1541.011763 348 -1802.524181 -1206.175433 349 -2399.833242 -1802.524181 350 -2423.524181 -2399.833242 351 -2393.524181 -2423.524181 352 -2420.427824 -2393.524181 353 -1677.524181 -2420.427824 354 -2783.410082 -1677.524181 355 -2423.524181 -2783.410082 356 -2423.524181 -2423.524181 357 -2423.524181 -2423.524181 358 1726.475819 -2423.524181 359 -1706.538957 1726.475819 360 -5724.432563 -1706.538957 361 -422.422949 -5724.432563 362 -1697.633094 -422.422949 363 -3275.123340 -1697.633094 364 10305.927627 -3275.123340 365 -3374.945564 10305.927627 366 -2354.342898 -3374.945564 367 -2407.444889 -2354.342898 368 2685.475819 -2407.444889 369 -1840.366319 2685.475819 370 -4392.364717 -1840.366319 371 -2098.901984 -4392.364717 372 1522.439542 -2098.901984 373 -2948.466625 1522.439542 374 -1219.976240 -2948.466625 375 422.830876 -1219.976240 376 20735.894586 422.830876 377 -5567.945758 20735.894586 378 -2560.481408 -5567.945758 379 2204.475819 -2560.481408 380 -2785.024837 2204.475819 381 -818.779707 -2785.024837 382 -2534.295064 -818.779707 383 -2547.381315 -2534.295064 384 -1443.542997 -2547.381315 385 744.774172 -1443.542997 386 -2575.526532 744.774172 387 -2423.524181 -2575.526532 388 -2423.524181 -2423.524181 389 -2423.524181 -2423.524181 390 -2423.524181 -2423.524181 391 10829.475819 -2423.524181 392 -3827.768280 10829.475819 393 -2382.718857 -3827.768280 394 -2423.524181 -2382.718857 395 -2423.524181 -2423.524181 396 -2423.524181 -2423.524181 397 -2292.524181 -2423.524181 398 -1894.617602 -2292.524181 399 -1145.395706 -1894.617602 400 4783.601521 -1145.395706 401 -3455.138968 4783.601521 402 15689.390287 -3455.138968 403 -5201.251309 15689.390287 404 -1636.330421 -5201.251309 405 -2203.848003 -1636.330421 406 946.373432 -2203.848003 407 1002.420813 946.373432 408 4344.352928 1002.420813 409 -4866.163850 4344.352928 410 -2478.097288 -4866.163850 411 -2423.524181 -2478.097288 412 -1697.524181 -2423.524181 413 -1223.995912 -1697.524181 414 -2345.939710 -1223.995912 415 4908.645920 -2345.939710 416 -7869.476244 4908.645920 417 -3174.521760 -7869.476244 418 -2423.732586 -3174.521760 419 -2327.524181 -2423.732586 420 -513.935683 -2327.524181 421 43696.445781 -513.935683 422 -14679.127902 43696.445781 423 -1089.826295 -14679.127902 424 3833.711605 -1089.826295 425 -1170.437279 3833.711605 426 -2197.524181 -1170.437279 427 -2361.745002 -2197.524181 428 -2907.162935 -2361.745002 429 -2174.524181 -2907.162935 430 -636.699222 -2174.524181 431 NA -636.699222 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -1126.182058 -446.799014 [2,] -511.543395 -1126.182058 [3,] -902.728119 -511.543395 [4,] 120.760443 -902.728119 [5,] -344.107520 120.760443 [6,] -381.130939 -344.107520 [7,] -524.381751 -381.130939 [8,] -660.956188 -524.381751 [9,] -460.062309 -660.956188 [10,] -894.317669 -460.062309 [11,] -904.534116 -894.317669 [12,] -531.752996 -904.534116 [13,] -373.326204 -531.752996 [14,] -494.215890 -373.326204 [15,] -369.054268 -494.215890 [16,] -522.363459 -369.054268 [17,] -541.639199 -522.363459 [18,] -543.244901 -541.639199 [19,] -522.622415 -543.244901 [20,] -531.752996 -522.622415 [21,] -385.786972 -531.752996 [22,] -531.752996 -385.786972 [23,] -1243.621731 -531.752996 [24,] -743.010377 -1243.621731 [25,] -515.477844 -743.010377 [26,] -408.176768 -515.477844 [27,] -784.862499 -408.176768 [28,] -531.752996 -784.862499 [29,] -659.261369 -531.752996 [30,] 91.977173 -659.261369 [31,] -531.752996 91.977173 [32,] -590.863508 -531.752996 [33,] -1726.837388 -590.863508 [34,] -307.246345 -1726.837388 [35,] -140.114489 -307.246345 [36,] -529.762446 -140.114489 [37,] -180.711314 -529.762446 [38,] -668.761908 -180.711314 [39,] -835.720751 -668.761908 [40,] -1228.182183 -835.720751 [41,] -781.072067 -1228.182183 [42,] -914.536718 -781.072067 [43,] -531.752996 -914.536718 [44,] -449.831058 -531.752996 [45,] -436.851090 -449.831058 [46,] -531.752996 -436.851090 [47,] -492.797868 -531.752996 [48,] -646.497911 -492.797868 [49,] -957.033521 -646.497911 [50,] -240.473167 -957.033521 [51,] -701.717963 -240.473167 [52,] -219.629339 -701.717963 [53,] -581.289588 -219.629339 [54,] -3383.558440 -581.289588 [55,] -766.028236 -3383.558440 [56,] -511.069047 -766.028236 [57,] -97.596255 -511.069047 [58,] -357.909723 -97.596255 [59,] -531.752996 -357.909723 [60,] -531.752996 -531.752996 [61,] -526.747763 -531.752996 [62,] 145.084534 -526.747763 [63,] -534.244498 145.084534 [64,] 655.246817 -534.244498 [65,] -526.236610 655.246817 [66,] -856.532957 -526.236610 [67,] -602.958505 -856.532957 [68,] -457.130721 -602.958505 [69,] -531.752996 -457.130721 [70,] -3174.838918 -531.752996 [71,] 211.927550 -3174.838918 [72,] 9354.448906 211.927550 [73,] -531.752996 9354.448906 [74,] -578.746650 -531.752996 [75,] -581.164370 -578.746650 [76,] -633.399028 -581.164370 [77,] -518.399530 -633.399028 [78,] 1113.396719 -518.399530 [79,] -37.206470 1113.396719 [80,] -529.007476 -37.206470 [81,] -531.752996 -529.007476 [82,] -531.752996 -531.752996 [83,] -554.257189 -531.752996 [84,] -531.752996 -554.257189 [85,] -531.752996 -531.752996 [86,] -521.732620 -531.752996 [87,] -531.449203 -521.732620 [88,] -444.201528 -531.449203 [89,] -670.889173 -444.201528 [90,] -531.752996 -670.889173 [91,] -401.698129 -531.752996 [92,] -1291.243244 -401.698129 [93,] 465.050125 -1291.243244 [94,] -563.965545 465.050125 [95,] -630.510274 -563.965545 [96,] -531.752996 -630.510274 [97,] -514.103829 -531.752996 [98,] -542.910900 -514.103829 [99,] -264.378147 -542.910900 [100,] -1074.155039 -264.378147 [101,] -434.808579 -1074.155039 [102,] -542.900405 -434.808579 [103,] -531.755114 -542.900405 [104,] -531.752996 -531.755114 [105,] 2029.650731 -531.752996 [106,] -431.446382 2029.650731 [107,] -1695.157848 -431.446382 [108,] -384.687406 -1695.157848 [109,] -510.562047 -384.687406 [110,] -518.842890 -510.562047 [111,] 51.592784 -518.842890 [112,] -531.752996 51.592784 [113,] -525.641602 -531.752996 [114,] -531.752996 -525.641602 [115,] -461.320462 -531.752996 [116,] -888.829048 -461.320462 [117,] -1887.590163 -888.829048 [118,] -2138.572831 -1887.590163 [119,] -1490.581032 -2138.572831 [120,] 3.856788 -1490.581032 [121,] -531.752996 3.856788 [122,] -2095.853648 -531.752996 [123,] -878.399845 -2095.853648 [124,] -20.323604 -878.399845 [125,] -531.752996 -20.323604 [126,] -545.554349 -531.752996 [127,] -531.752996 -545.554349 [128,] 1523.730439 -531.752996 [129,] -476.086445 1523.730439 [130,] -585.320451 -476.086445 [131,] -531.752996 -585.320451 [132,] -531.752996 -531.752996 [133,] -531.752996 -531.752996 [134,] -520.235254 -531.752996 [135,] -472.722529 -520.235254 [136,] -2088.524181 -472.722529 [137,] -2592.868414 -2088.524181 [138,] -2423.524181 -2592.868414 [139,] 2101.475819 -2423.524181 [140,] -1404.972548 2101.475819 [141,] -575.483408 -1404.972548 [142,] -2423.524181 -575.483408 [143,] -1785.524181 -2423.524181 [144,] -2457.990409 -1785.524181 [145,] -1816.524181 -2457.990409 [146,] -871.909281 -1816.524181 [147,] -1019.813082 -871.909281 [148,] -2653.563065 -1019.813082 [149,] -2423.524181 -2653.563065 [150,] -1812.524181 -2423.524181 [151,] -2300.277558 -1812.524181 [152,] -1500.524181 -2300.277558 [153,] -1648.256178 -1500.524181 [154,] -64.600804 -1648.256178 [155,] -2586.926486 -64.600804 [156,] -2553.057962 -2586.926486 [157,] -2288.524181 -2553.057962 [158,] -822.180481 -2288.524181 [159,] -4711.001339 -822.180481 [160,] -2412.082568 -4711.001339 [161,] -2114.524181 -2412.082568 [162,] -2176.304013 -2114.524181 [163,] -2596.137187 -2176.304013 [164,] 2095.475819 -2596.137187 [165,] -5116.656277 2095.475819 [166,] -1586.524181 -5116.656277 [167,] 2674.509898 -1586.524181 [168,] -336.323074 2674.509898 [169,] -240.871400 -336.323074 [170,] -3787.369387 -240.871400 [171,] -2440.813076 -3787.369387 [172,] -1905.689883 -2440.813076 [173,] -2410.750338 -1905.689883 [174,] -1495.011601 -2410.750338 [175,] -3495.484878 -1495.011601 [176,] -2362.524181 -3495.484878 [177,] -2052.208190 -2362.524181 [178,] -2431.352844 -2052.208190 [179,] 2477.475819 -2431.352844 [180,] -5103.395872 2477.475819 [181,] -4500.593947 -5103.395872 [182,] -2423.524181 -4500.593947 [183,] -2423.524181 -2423.524181 [184,] -2421.524181 -2423.524181 [185,] -2387.648451 -2421.524181 [186,] -1677.093216 -2387.648451 [187,] 82103.337225 -1677.093216 [188,] 1803.203792 82103.337225 [189,] -2420.524181 1803.203792 [190,] -1896.916885 -2420.524181 [191,] -2511.118916 -1896.916885 [192,] -2018.524181 -2511.118916 [193,] -1517.305540 -2018.524181 [194,] -3044.373271 -1517.305540 [195,] -324.524181 -3044.373271 [196,] -1049.596659 -324.524181 [197,] -1391.588625 -1049.596659 [198,] -2423.524181 -1391.588625 [199,] 19906.475819 -2423.524181 [200,] -4429.102815 19906.475819 [201,] -2423.524181 -4429.102815 [202,] -1940.524181 -2423.524181 [203,] -2400.509309 -1940.524181 [204,] -184.524181 -2400.509309 [205,] -1107.305949 -184.524181 [206,] -3611.207824 -1107.305949 [207,] -2058.524181 -3611.207824 [208,] -1250.547881 -2058.524181 [209,] 19757.475870 -1250.547881 [210,] -3095.283585 19757.475870 [211,] -2685.572436 -3095.283585 [212,] -1892.462246 -2685.572436 [213,] 11022.287859 -1892.462246 [214,] -6015.116649 11022.287859 [215,] -2776.875487 -6015.116649 [216,] -2468.180819 -2776.875487 [217,] -1009.524181 -2468.180819 [218,] -909.398752 -1009.524181 [219,] -5250.609171 -909.398752 [220,] -2423.625909 -5250.609171 [221,] -1579.524181 -2423.625909 [222,] -2068.290876 -1579.524181 [223,] -2009.807846 -2068.290876 [224,] 16285.157028 -2009.807846 [225,] -4485.171685 16285.157028 [226,] 4699.475819 -4485.171685 [227,] -13070.996208 4699.475819 [228,] -2771.687434 -13070.996208 [229,] -366.986238 -2771.687434 [230,] -14203.913927 -366.986238 [231,] 9661.176538 -14203.913927 [232,] -5228.101162 9661.176538 [233,] -2298.524181 -5228.101162 [234,] -1229.369359 -2298.524181 [235,] -1406.605618 -1229.369359 [236,] -2303.724111 -1406.605618 [237,] -2423.524181 -2303.724111 [238,] -2393.524181 -2423.524181 [239,] 818.673002 -2393.524181 [240,] -2909.410634 818.673002 [241,] -3296.891128 -2909.410634 [242,] -1955.524181 -3296.891128 [243,] 2504.134858 -1955.524181 [244,] 1401.321731 2504.134858 [245,] -1525.114460 1401.321731 [246,] -2298.524181 -1525.114460 [247,] -2384.498284 -2298.524181 [248,] -2401.524181 -2384.498284 [249,] -2415.199970 -2401.524181 [250,] -1909.524181 -2415.199970 [251,] -3261.382924 -1909.524181 [252,] -2423.524181 -3261.382924 [253,] -2423.524181 -2423.524181 [254,] -2423.524181 -2423.524181 [255,] -2423.524181 -2423.524181 [256,] -2423.524181 -2423.524181 [257,] 7903.475819 -2423.524181 [258,] 9098.228752 7903.475819 [259,] -384.019079 9098.228752 [260,] 11492.941439 -384.019079 [261,] -738.767412 11492.941439 [262,] -2043.565720 -738.767412 [263,] -2237.191848 -2043.565720 [264,] -2748.599377 -2237.191848 [265,] -2055.524181 -2748.599377 [266,] -2390.231484 -2055.524181 [267,] -2090.524181 -2390.231484 [268,] -2901.371416 -2090.524181 [269,] -2427.143658 -2901.371416 [270,] -2423.524181 -2427.143658 [271,] -1120.524181 -2423.524181 [272,] -2798.025251 -1120.524181 [273,] -48.848634 -2798.025251 [274,] -2503.611189 -48.848634 [275,] -1858.364649 -2503.611189 [276,] -2373.885174 -1858.364649 [277,] -1768.066474 -2373.885174 [278,] -1606.209373 -1768.066474 [279,] -1638.563020 -1606.209373 [280,] -650.680867 -1638.563020 [281,] -2441.671267 -650.680867 [282,] -2106.524181 -2441.671267 [283,] -2519.199175 -2106.524181 [284,] -2148.056984 -2519.199175 [285,] -2377.469561 -2148.056984 [286,] -1487.524181 -2377.469561 [287,] 5549.105588 -1487.524181 [288,] 9047.715220 5549.105588 [289,] 2150.713080 9047.715220 [290,] -1685.524181 2150.713080 [291,] -2440.114833 -1685.524181 [292,] -1831.524181 -2440.114833 [293,] -2347.712027 -1831.524181 [294,] -2431.843668 -2347.712027 [295,] -2424.811446 -2431.843668 [296,] 15590.475819 -2424.811446 [297,] -2710.117007 15590.475819 [298,] -2423.524181 -2710.117007 [299,] 35280.475819 -2423.524181 [300,] -4252.290035 35280.475819 [301,] -1227.341467 -4252.290035 [302,] -4413.234993 -1227.341467 [303,] -2253.250875 -4413.234993 [304,] -2589.665167 -2253.250875 [305,] 4480.022951 -2589.665167 [306,] -5185.573175 4480.022951 [307,] -2423.524181 -5185.573175 [308,] -2422.524181 -2423.524181 [309,] -2423.401116 -2422.524181 [310,] -95.524181 -2423.401116 [311,] -3031.324247 -95.524181 [312,] -2548.173266 -3031.324247 [313,] -2256.007415 -2548.173266 [314,] -825.309520 -2256.007415 [315,] 84708.824961 -825.309520 [316,] -17094.684688 84708.824961 [317,] -2260.448538 -17094.684688 [318,] -1410.524181 -2260.448538 [319,] -2641.648261 -1410.524181 [320,] -1841.926021 -2641.648261 [321,] 2969.730390 -1841.926021 [322,] -4023.560995 2969.730390 [323,] -1411.524181 -4023.560995 [324,] -5029.710777 -1411.524181 [325,] -1547.524181 -5029.710777 [326,] 156470.426224 -1547.524181 [327,] 16817.902759 156470.426224 [328,] 2543.542624 16817.902759 [329,] -2563.394724 2543.542624 [330,] -2799.065069 -2563.394724 [331,] -2423.524181 -2799.065069 [332,] -58.524181 -2423.524181 [333,] -7793.933302 -58.524181 [334,] -3098.755366 -7793.933302 [335,] -838.524181 -3098.755366 [336,] 12183.821560 -838.524181 [337,] -5061.195168 12183.821560 [338,] -2423.524181 -5061.195168 [339,] 8155.475819 -2423.524181 [340,] -4370.342296 8155.475819 [341,] -2403.391098 -4370.342296 [342,] 1218.475819 -2403.391098 [343,] -4172.590994 1218.475819 [344,] -1951.524181 -4172.590994 [345,] -3071.107195 -1951.524181 [346,] 1541.011763 -3071.107195 [347,] -1206.175433 1541.011763 [348,] -1802.524181 -1206.175433 [349,] -2399.833242 -1802.524181 [350,] -2423.524181 -2399.833242 [351,] -2393.524181 -2423.524181 [352,] -2420.427824 -2393.524181 [353,] -1677.524181 -2420.427824 [354,] -2783.410082 -1677.524181 [355,] -2423.524181 -2783.410082 [356,] -2423.524181 -2423.524181 [357,] -2423.524181 -2423.524181 [358,] 1726.475819 -2423.524181 [359,] -1706.538957 1726.475819 [360,] -5724.432563 -1706.538957 [361,] -422.422949 -5724.432563 [362,] -1697.633094 -422.422949 [363,] -3275.123340 -1697.633094 [364,] 10305.927627 -3275.123340 [365,] -3374.945564 10305.927627 [366,] -2354.342898 -3374.945564 [367,] -2407.444889 -2354.342898 [368,] 2685.475819 -2407.444889 [369,] -1840.366319 2685.475819 [370,] -4392.364717 -1840.366319 [371,] -2098.901984 -4392.364717 [372,] 1522.439542 -2098.901984 [373,] -2948.466625 1522.439542 [374,] -1219.976240 -2948.466625 [375,] 422.830876 -1219.976240 [376,] 20735.894586 422.830876 [377,] -5567.945758 20735.894586 [378,] -2560.481408 -5567.945758 [379,] 2204.475819 -2560.481408 [380,] -2785.024837 2204.475819 [381,] -818.779707 -2785.024837 [382,] -2534.295064 -818.779707 [383,] -2547.381315 -2534.295064 [384,] -1443.542997 -2547.381315 [385,] 744.774172 -1443.542997 [386,] -2575.526532 744.774172 [387,] -2423.524181 -2575.526532 [388,] -2423.524181 -2423.524181 [389,] -2423.524181 -2423.524181 [390,] -2423.524181 -2423.524181 [391,] 10829.475819 -2423.524181 [392,] -3827.768280 10829.475819 [393,] -2382.718857 -3827.768280 [394,] -2423.524181 -2382.718857 [395,] -2423.524181 -2423.524181 [396,] -2423.524181 -2423.524181 [397,] -2292.524181 -2423.524181 [398,] -1894.617602 -2292.524181 [399,] -1145.395706 -1894.617602 [400,] 4783.601521 -1145.395706 [401,] -3455.138968 4783.601521 [402,] 15689.390287 -3455.138968 [403,] -5201.251309 15689.390287 [404,] -1636.330421 -5201.251309 [405,] -2203.848003 -1636.330421 [406,] 946.373432 -2203.848003 [407,] 1002.420813 946.373432 [408,] 4344.352928 1002.420813 [409,] -4866.163850 4344.352928 [410,] -2478.097288 -4866.163850 [411,] -2423.524181 -2478.097288 [412,] -1697.524181 -2423.524181 [413,] -1223.995912 -1697.524181 [414,] -2345.939710 -1223.995912 [415,] 4908.645920 -2345.939710 [416,] -7869.476244 4908.645920 [417,] -3174.521760 -7869.476244 [418,] -2423.732586 -3174.521760 [419,] -2327.524181 -2423.732586 [420,] -513.935683 -2327.524181 [421,] 43696.445781 -513.935683 [422,] -14679.127902 43696.445781 [423,] -1089.826295 -14679.127902 [424,] 3833.711605 -1089.826295 [425,] -1170.437279 3833.711605 [426,] -2197.524181 -1170.437279 [427,] -2361.745002 -2197.524181 [428,] -2907.162935 -2361.745002 [429,] -2174.524181 -2907.162935 [430,] -636.699222 -2174.524181 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -1126.182058 -446.799014 2 -511.543395 -1126.182058 3 -902.728119 -511.543395 4 120.760443 -902.728119 5 -344.107520 120.760443 6 -381.130939 -344.107520 7 -524.381751 -381.130939 8 -660.956188 -524.381751 9 -460.062309 -660.956188 10 -894.317669 -460.062309 11 -904.534116 -894.317669 12 -531.752996 -904.534116 13 -373.326204 -531.752996 14 -494.215890 -373.326204 15 -369.054268 -494.215890 16 -522.363459 -369.054268 17 -541.639199 -522.363459 18 -543.244901 -541.639199 19 -522.622415 -543.244901 20 -531.752996 -522.622415 21 -385.786972 -531.752996 22 -531.752996 -385.786972 23 -1243.621731 -531.752996 24 -743.010377 -1243.621731 25 -515.477844 -743.010377 26 -408.176768 -515.477844 27 -784.862499 -408.176768 28 -531.752996 -784.862499 29 -659.261369 -531.752996 30 91.977173 -659.261369 31 -531.752996 91.977173 32 -590.863508 -531.752996 33 -1726.837388 -590.863508 34 -307.246345 -1726.837388 35 -140.114489 -307.246345 36 -529.762446 -140.114489 37 -180.711314 -529.762446 38 -668.761908 -180.711314 39 -835.720751 -668.761908 40 -1228.182183 -835.720751 41 -781.072067 -1228.182183 42 -914.536718 -781.072067 43 -531.752996 -914.536718 44 -449.831058 -531.752996 45 -436.851090 -449.831058 46 -531.752996 -436.851090 47 -492.797868 -531.752996 48 -646.497911 -492.797868 49 -957.033521 -646.497911 50 -240.473167 -957.033521 51 -701.717963 -240.473167 52 -219.629339 -701.717963 53 -581.289588 -219.629339 54 -3383.558440 -581.289588 55 -766.028236 -3383.558440 56 -511.069047 -766.028236 57 -97.596255 -511.069047 58 -357.909723 -97.596255 59 -531.752996 -357.909723 60 -531.752996 -531.752996 61 -526.747763 -531.752996 62 145.084534 -526.747763 63 -534.244498 145.084534 64 655.246817 -534.244498 65 -526.236610 655.246817 66 -856.532957 -526.236610 67 -602.958505 -856.532957 68 -457.130721 -602.958505 69 -531.752996 -457.130721 70 -3174.838918 -531.752996 71 211.927550 -3174.838918 72 9354.448906 211.927550 73 -531.752996 9354.448906 74 -578.746650 -531.752996 75 -581.164370 -578.746650 76 -633.399028 -581.164370 77 -518.399530 -633.399028 78 1113.396719 -518.399530 79 -37.206470 1113.396719 80 -529.007476 -37.206470 81 -531.752996 -529.007476 82 -531.752996 -531.752996 83 -554.257189 -531.752996 84 -531.752996 -554.257189 85 -531.752996 -531.752996 86 -521.732620 -531.752996 87 -531.449203 -521.732620 88 -444.201528 -531.449203 89 -670.889173 -444.201528 90 -531.752996 -670.889173 91 -401.698129 -531.752996 92 -1291.243244 -401.698129 93 465.050125 -1291.243244 94 -563.965545 465.050125 95 -630.510274 -563.965545 96 -531.752996 -630.510274 97 -514.103829 -531.752996 98 -542.910900 -514.103829 99 -264.378147 -542.910900 100 -1074.155039 -264.378147 101 -434.808579 -1074.155039 102 -542.900405 -434.808579 103 -531.755114 -542.900405 104 -531.752996 -531.755114 105 2029.650731 -531.752996 106 -431.446382 2029.650731 107 -1695.157848 -431.446382 108 -384.687406 -1695.157848 109 -510.562047 -384.687406 110 -518.842890 -510.562047 111 51.592784 -518.842890 112 -531.752996 51.592784 113 -525.641602 -531.752996 114 -531.752996 -525.641602 115 -461.320462 -531.752996 116 -888.829048 -461.320462 117 -1887.590163 -888.829048 118 -2138.572831 -1887.590163 119 -1490.581032 -2138.572831 120 3.856788 -1490.581032 121 -531.752996 3.856788 122 -2095.853648 -531.752996 123 -878.399845 -2095.853648 124 -20.323604 -878.399845 125 -531.752996 -20.323604 126 -545.554349 -531.752996 127 -531.752996 -545.554349 128 1523.730439 -531.752996 129 -476.086445 1523.730439 130 -585.320451 -476.086445 131 -531.752996 -585.320451 132 -531.752996 -531.752996 133 -531.752996 -531.752996 134 -520.235254 -531.752996 135 -472.722529 -520.235254 136 -2088.524181 -472.722529 137 -2592.868414 -2088.524181 138 -2423.524181 -2592.868414 139 2101.475819 -2423.524181 140 -1404.972548 2101.475819 141 -575.483408 -1404.972548 142 -2423.524181 -575.483408 143 -1785.524181 -2423.524181 144 -2457.990409 -1785.524181 145 -1816.524181 -2457.990409 146 -871.909281 -1816.524181 147 -1019.813082 -871.909281 148 -2653.563065 -1019.813082 149 -2423.524181 -2653.563065 150 -1812.524181 -2423.524181 151 -2300.277558 -1812.524181 152 -1500.524181 -2300.277558 153 -1648.256178 -1500.524181 154 -64.600804 -1648.256178 155 -2586.926486 -64.600804 156 -2553.057962 -2586.926486 157 -2288.524181 -2553.057962 158 -822.180481 -2288.524181 159 -4711.001339 -822.180481 160 -2412.082568 -4711.001339 161 -2114.524181 -2412.082568 162 -2176.304013 -2114.524181 163 -2596.137187 -2176.304013 164 2095.475819 -2596.137187 165 -5116.656277 2095.475819 166 -1586.524181 -5116.656277 167 2674.509898 -1586.524181 168 -336.323074 2674.509898 169 -240.871400 -336.323074 170 -3787.369387 -240.871400 171 -2440.813076 -3787.369387 172 -1905.689883 -2440.813076 173 -2410.750338 -1905.689883 174 -1495.011601 -2410.750338 175 -3495.484878 -1495.011601 176 -2362.524181 -3495.484878 177 -2052.208190 -2362.524181 178 -2431.352844 -2052.208190 179 2477.475819 -2431.352844 180 -5103.395872 2477.475819 181 -4500.593947 -5103.395872 182 -2423.524181 -4500.593947 183 -2423.524181 -2423.524181 184 -2421.524181 -2423.524181 185 -2387.648451 -2421.524181 186 -1677.093216 -2387.648451 187 82103.337225 -1677.093216 188 1803.203792 82103.337225 189 -2420.524181 1803.203792 190 -1896.916885 -2420.524181 191 -2511.118916 -1896.916885 192 -2018.524181 -2511.118916 193 -1517.305540 -2018.524181 194 -3044.373271 -1517.305540 195 -324.524181 -3044.373271 196 -1049.596659 -324.524181 197 -1391.588625 -1049.596659 198 -2423.524181 -1391.588625 199 19906.475819 -2423.524181 200 -4429.102815 19906.475819 201 -2423.524181 -4429.102815 202 -1940.524181 -2423.524181 203 -2400.509309 -1940.524181 204 -184.524181 -2400.509309 205 -1107.305949 -184.524181 206 -3611.207824 -1107.305949 207 -2058.524181 -3611.207824 208 -1250.547881 -2058.524181 209 19757.475870 -1250.547881 210 -3095.283585 19757.475870 211 -2685.572436 -3095.283585 212 -1892.462246 -2685.572436 213 11022.287859 -1892.462246 214 -6015.116649 11022.287859 215 -2776.875487 -6015.116649 216 -2468.180819 -2776.875487 217 -1009.524181 -2468.180819 218 -909.398752 -1009.524181 219 -5250.609171 -909.398752 220 -2423.625909 -5250.609171 221 -1579.524181 -2423.625909 222 -2068.290876 -1579.524181 223 -2009.807846 -2068.290876 224 16285.157028 -2009.807846 225 -4485.171685 16285.157028 226 4699.475819 -4485.171685 227 -13070.996208 4699.475819 228 -2771.687434 -13070.996208 229 -366.986238 -2771.687434 230 -14203.913927 -366.986238 231 9661.176538 -14203.913927 232 -5228.101162 9661.176538 233 -2298.524181 -5228.101162 234 -1229.369359 -2298.524181 235 -1406.605618 -1229.369359 236 -2303.724111 -1406.605618 237 -2423.524181 -2303.724111 238 -2393.524181 -2423.524181 239 818.673002 -2393.524181 240 -2909.410634 818.673002 241 -3296.891128 -2909.410634 242 -1955.524181 -3296.891128 243 2504.134858 -1955.524181 244 1401.321731 2504.134858 245 -1525.114460 1401.321731 246 -2298.524181 -1525.114460 247 -2384.498284 -2298.524181 248 -2401.524181 -2384.498284 249 -2415.199970 -2401.524181 250 -1909.524181 -2415.199970 251 -3261.382924 -1909.524181 252 -2423.524181 -3261.382924 253 -2423.524181 -2423.524181 254 -2423.524181 -2423.524181 255 -2423.524181 -2423.524181 256 -2423.524181 -2423.524181 257 7903.475819 -2423.524181 258 9098.228752 7903.475819 259 -384.019079 9098.228752 260 11492.941439 -384.019079 261 -738.767412 11492.941439 262 -2043.565720 -738.767412 263 -2237.191848 -2043.565720 264 -2748.599377 -2237.191848 265 -2055.524181 -2748.599377 266 -2390.231484 -2055.524181 267 -2090.524181 -2390.231484 268 -2901.371416 -2090.524181 269 -2427.143658 -2901.371416 270 -2423.524181 -2427.143658 271 -1120.524181 -2423.524181 272 -2798.025251 -1120.524181 273 -48.848634 -2798.025251 274 -2503.611189 -48.848634 275 -1858.364649 -2503.611189 276 -2373.885174 -1858.364649 277 -1768.066474 -2373.885174 278 -1606.209373 -1768.066474 279 -1638.563020 -1606.209373 280 -650.680867 -1638.563020 281 -2441.671267 -650.680867 282 -2106.524181 -2441.671267 283 -2519.199175 -2106.524181 284 -2148.056984 -2519.199175 285 -2377.469561 -2148.056984 286 -1487.524181 -2377.469561 287 5549.105588 -1487.524181 288 9047.715220 5549.105588 289 2150.713080 9047.715220 290 -1685.524181 2150.713080 291 -2440.114833 -1685.524181 292 -1831.524181 -2440.114833 293 -2347.712027 -1831.524181 294 -2431.843668 -2347.712027 295 -2424.811446 -2431.843668 296 15590.475819 -2424.811446 297 -2710.117007 15590.475819 298 -2423.524181 -2710.117007 299 35280.475819 -2423.524181 300 -4252.290035 35280.475819 301 -1227.341467 -4252.290035 302 -4413.234993 -1227.341467 303 -2253.250875 -4413.234993 304 -2589.665167 -2253.250875 305 4480.022951 -2589.665167 306 -5185.573175 4480.022951 307 -2423.524181 -5185.573175 308 -2422.524181 -2423.524181 309 -2423.401116 -2422.524181 310 -95.524181 -2423.401116 311 -3031.324247 -95.524181 312 -2548.173266 -3031.324247 313 -2256.007415 -2548.173266 314 -825.309520 -2256.007415 315 84708.824961 -825.309520 316 -17094.684688 84708.824961 317 -2260.448538 -17094.684688 318 -1410.524181 -2260.448538 319 -2641.648261 -1410.524181 320 -1841.926021 -2641.648261 321 2969.730390 -1841.926021 322 -4023.560995 2969.730390 323 -1411.524181 -4023.560995 324 -5029.710777 -1411.524181 325 -1547.524181 -5029.710777 326 156470.426224 -1547.524181 327 16817.902759 156470.426224 328 2543.542624 16817.902759 329 -2563.394724 2543.542624 330 -2799.065069 -2563.394724 331 -2423.524181 -2799.065069 332 -58.524181 -2423.524181 333 -7793.933302 -58.524181 334 -3098.755366 -7793.933302 335 -838.524181 -3098.755366 336 12183.821560 -838.524181 337 -5061.195168 12183.821560 338 -2423.524181 -5061.195168 339 8155.475819 -2423.524181 340 -4370.342296 8155.475819 341 -2403.391098 -4370.342296 342 1218.475819 -2403.391098 343 -4172.590994 1218.475819 344 -1951.524181 -4172.590994 345 -3071.107195 -1951.524181 346 1541.011763 -3071.107195 347 -1206.175433 1541.011763 348 -1802.524181 -1206.175433 349 -2399.833242 -1802.524181 350 -2423.524181 -2399.833242 351 -2393.524181 -2423.524181 352 -2420.427824 -2393.524181 353 -1677.524181 -2420.427824 354 -2783.410082 -1677.524181 355 -2423.524181 -2783.410082 356 -2423.524181 -2423.524181 357 -2423.524181 -2423.524181 358 1726.475819 -2423.524181 359 -1706.538957 1726.475819 360 -5724.432563 -1706.538957 361 -422.422949 -5724.432563 362 -1697.633094 -422.422949 363 -3275.123340 -1697.633094 364 10305.927627 -3275.123340 365 -3374.945564 10305.927627 366 -2354.342898 -3374.945564 367 -2407.444889 -2354.342898 368 2685.475819 -2407.444889 369 -1840.366319 2685.475819 370 -4392.364717 -1840.366319 371 -2098.901984 -4392.364717 372 1522.439542 -2098.901984 373 -2948.466625 1522.439542 374 -1219.976240 -2948.466625 375 422.830876 -1219.976240 376 20735.894586 422.830876 377 -5567.945758 20735.894586 378 -2560.481408 -5567.945758 379 2204.475819 -2560.481408 380 -2785.024837 2204.475819 381 -818.779707 -2785.024837 382 -2534.295064 -818.779707 383 -2547.381315 -2534.295064 384 -1443.542997 -2547.381315 385 744.774172 -1443.542997 386 -2575.526532 744.774172 387 -2423.524181 -2575.526532 388 -2423.524181 -2423.524181 389 -2423.524181 -2423.524181 390 -2423.524181 -2423.524181 391 10829.475819 -2423.524181 392 -3827.768280 10829.475819 393 -2382.718857 -3827.768280 394 -2423.524181 -2382.718857 395 -2423.524181 -2423.524181 396 -2423.524181 -2423.524181 397 -2292.524181 -2423.524181 398 -1894.617602 -2292.524181 399 -1145.395706 -1894.617602 400 4783.601521 -1145.395706 401 -3455.138968 4783.601521 402 15689.390287 -3455.138968 403 -5201.251309 15689.390287 404 -1636.330421 -5201.251309 405 -2203.848003 -1636.330421 406 946.373432 -2203.848003 407 1002.420813 946.373432 408 4344.352928 1002.420813 409 -4866.163850 4344.352928 410 -2478.097288 -4866.163850 411 -2423.524181 -2478.097288 412 -1697.524181 -2423.524181 413 -1223.995912 -1697.524181 414 -2345.939710 -1223.995912 415 4908.645920 -2345.939710 416 -7869.476244 4908.645920 417 -3174.521760 -7869.476244 418 -2423.732586 -3174.521760 419 -2327.524181 -2423.732586 420 -513.935683 -2327.524181 421 43696.445781 -513.935683 422 -14679.127902 43696.445781 423 -1089.826295 -14679.127902 424 3833.711605 -1089.826295 425 -1170.437279 3833.711605 426 -2197.524181 -1170.437279 427 -2361.745002 -2197.524181 428 -2907.162935 -2361.745002 429 -2174.524181 -2907.162935 430 -636.699222 -2174.524181 > 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/rcomp/tmp/7qr331292675864.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/rcomp/tmp/8qr331292675864.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/rcomp/tmp/9ctnj1292675865.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/rcomp/tmp/10ctnj1292675865.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/rcomp/tmp/11ftl71292675865.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/rcomp/tmp/120c1v1292675865.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/rcomp/tmp/13pdhp1292675865.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/rcomp/tmp/14i4yr1292675865.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/rcomp/tmp/15lnwx1292675865.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/rcomp/tmp/16zwuo1292675865.tab") + } > > try(system("convert tmp/1uinc1292675864.ps tmp/1uinc1292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/2nr4f1292675864.ps tmp/2nr4f1292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/3nr4f1292675864.ps tmp/3nr4f1292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/4nr4f1292675864.ps tmp/4nr4f1292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/5gi401292675864.ps tmp/5gi401292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/6gi401292675864.ps tmp/6gi401292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/7qr331292675864.ps tmp/7qr331292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/8qr331292675864.ps tmp/8qr331292675864.png",intern=TRUE)) character(0) > try(system("convert tmp/9ctnj1292675865.ps tmp/9ctnj1292675865.png",intern=TRUE)) character(0) > try(system("convert tmp/10ctnj1292675865.ps tmp/10ctnj1292675865.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 11.790 1.860 13.702