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(807 + ,213118 + ,171986226 + ,6282154 + ,444 + ,81767 + ,36304548 + ,4321023 + ,412 + ,153198 + ,63117576 + ,4111912 + ,428 + ,-26007 + ,-11130996 + ,223193 + ,315 + ,126942 + ,39986730 + ,1491348 + ,168 + ,157214 + ,26411952 + ,1629616 + ,263 + ,129352 + ,34019576 + ,1398893 + ,267 + ,234817 + ,62696139 + ,1926517 + ,228 + ,60448 + ,13782144 + ,983660 + ,129 + ,47818 + ,6168522 + ,1443586 + ,104 + ,245546 + ,25536784 + ,1073089 + ,122 + ,48020 + ,5858440 + ,984885 + ,393 + ,-1710 + ,-672030 + ,1405225 + ,190 + ,32648 + ,6203120 + ,227132 + ,280 + ,95350 + ,26698000 + ,929118 + ,63 + ,151352 + ,9535176 + ,1071292 + ,102 + ,288170 + ,29393340 + ,638830 + ,265 + ,114337 + ,30299305 + ,856956 + ,234 + ,37884 + ,8864856 + ,992426 + ,277 + ,122844 + ,34027788 + ,444477 + ,73 + ,82340 + ,6010820 + ,857217 + ,67 + ,79801 + ,5346667 + ,711969 + ,103 + ,165548 + ,17051444 + ,702380 + ,290 + ,116384 + ,33751360 + ,358589 + ,83 + ,134028 + ,11124324 + ,297978 + ,56 + ,63838 + ,3574928 + ,585715 + ,236 + ,74996 + ,17699056 + ,657954 + ,73 + ,31080 + ,2268840 + ,209458 + ,34 + ,32168 + ,1093712 + ,786690 + ,139 + ,49857 + ,6930123 + ,439798 + ,26 + ,87161 + ,2266186 + ,688779 + ,70 + ,106113 + ,7427910 + ,574339 + ,40 + ,80570 + ,3222800 + ,741409 + ,42 + ,102129 + ,4289418 + ,597793 + ,12 + ,301670 + ,3620040 + ,644190 + ,211 + ,102313 + ,21588043 + ,377934 + ,74 + ,88577 + ,6554698 + ,640273 + ,80 + ,112477 + ,8998160 + ,697458 + ,83 + ,191778 + ,15917574 + ,550608 + ,131 + ,79804 + ,10454324 + ,207393 + ,203 + ,128294 + ,26043682 + ,301607 + ,56 + ,96448 + ,5401088 + ,345783 + ,89 + ,93811 + ,8349179 + ,501749 + ,88 + ,117520 + ,10341760 + ,379983 + ,39 + ,69159 + ,2697201 + ,387475 + ,25 + ,101792 + ,2544800 + ,377305 + ,49 + ,210568 + ,10317832 + ,370837 + ,149 + ,136996 + ,20412404 + ,430866 + ,58 + ,121920 + ,7071360 + ,469107 + ,41 + ,76403 + ,3132523 + ,194493 + ,90 + ,108094 + ,9728460 + ,530670 + ,136 + ,134759 + ,18327224 + ,518365 + ,97 + ,188873 + ,18320681 + ,491303 + ,63 + ,146216 + ,9211608 + ,527021 + ,114 + ,156608 + ,17853312 + ,233773 + ,77 + ,61348 + ,4723796 + ,405972 + ,6 + ,50350 + ,302100 + ,652925 + ,47 + ,87720 + ,4122840 + ,446211 + ,51 + ,99489 + ,5073939 + ,341340 + ,85 + ,87419 + ,7430615 + ,387699 + ,43 + ,94355 + ,4057265 + ,493408 + ,32 + ,60326 + ,1930432 + ,146494 + ,25 + ,94670 + ,2366750 + ,414462 + ,77 + ,82425 + ,6346725 + ,364304 + ,54 + ,59017 + ,3186918 + ,355178 + ,251 + ,90829 + ,22798079 + ,357760 + ,15 + ,80791 + ,1211865 + ,261216 + ,44 + ,100423 + ,4418612 + ,397144 + ,73 + ,131116 + ,9571468 + ,374943 + ,85 + ,100269 + ,8522865 + ,424898 + ,49 + ,27330 + ,1339170 + ,202055 + ,38 + ,39039 + ,1483482 + ,378525 + ,35 + ,106885 + ,3740975 + ,310768 + ,9 + ,79285 + ,713565 + ,325738 + ,34 + ,118881 + ,4041954 + ,394510 + ,20 + ,77623 + ,1552460 + ,247060 + ,29 + ,114768 + ,3328272 + ,368078 + ,11 + ,74015 + ,814165 + ,236761 + ,52 + ,69465 + ,3612180 + ,312378 + ,13 + ,117869 + ,1532297 + ,339836 + ,29 + ,60982 + ,1768478 + ,347385 + ,66 + ,90131 + ,5948646 + ,426280 + ,33 + ,138971 + ,4586043 + ,352850 + ,15 + ,39625 + ,594375 + ,301881 + ,15 + ,102725 + ,1540875 + ,377516 + ,68 + ,64239 + ,4368252 + ,357312 + ,100 + ,90262 + ,9026200 + ,458343 + ,13 + ,103960 + ,1351480 + ,354228 + ,45 + ,106611 + ,4797495 + ,308636 + ,14 + ,103345 + ,1446830 + ,386212 + ,36 + ,95551 + ,3439836 + ,393343 + ,40 + ,82903 + ,3316120 + ,378509 + ,68 + ,63593 + ,4324324 + ,452469 + ,29 + ,126910 + ,3680390 + ,364839 + ,43 + ,37527 + ,1613661 + ,358649 + ,30 + ,60247 + ,1807410 + ,376641 + ,9 + ,112995 + ,1016955 + ,429112 + ,22 + ,70184 + ,1544048 + ,330546 + ,19 + ,130140 + ,2472660 + ,403560 + ,9 + ,73221 + ,658989 + ,317892 + ,31 + ,76114 + ,2359534 + ,307528 + ,19 + ,90534 + ,1720146 + ,235133 + ,55 + ,108479 + ,5966345 + ,299243 + ,8 + ,113761 + ,910088 + ,314073 + ,28 + ,68696 + ,1923488 + ,368186 + ,29 + ,71561 + ,2075269 + ,269661 + ,48 + ,59831 + ,2871888 + ,125390 + ,16 + ,97890 + ,1566240 + ,510834 + ,47 + ,101481 + ,4769607 + ,321896 + ,20 + ,72954 + ,1459080 + ,249898 + ,22 + ,67939 + ,1494658 + ,408881 + ,33 + ,48022 + ,1584726 + ,158492 + ,44 + ,86111 + ,3788884 + ,292154 + ,13 + ,74020 + ,962260 + ,289513 + ,6 + ,57530 + ,345180 + ,378049 + ,35 + ,56364 + ,1972740 + ,343466 + ,8 + ,84990 + ,679920 + ,332743 + ,17 + ,88590 + ,1506030 + ,442882 + ,11 + ,77200 + ,849200 + ,214215 + ,21 + ,61262 + ,1286502 + ,315688 + ,92 + ,110309 + ,10148428 + ,375195 + ,12 + ,67000 + ,804000 + ,334280 + ,112 + ,93099 + ,10427088 + ,355864 + ,25 + ,107577 + ,2689425 + ,480382 + ,17 + ,62920 + ,1069640 + ,353058 + ,23 + ,75832 + ,1744136 + ,217193 + ,0 + ,60720 + ,0 + ,315380 + ,10 + ,60793 + ,607930 + ,314533 + ,23 + ,57935 + ,1332505 + ,318056 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,60630 + ,424410 + ,314353 + ,25 + ,55637 + ,1390925 + ,369448 + ,0 + ,60720 + ,0 + ,315380 + ,20 + ,60887 + ,1217740 + ,312846 + ,4 + ,60720 + ,242880 + ,312075 + ,4 + ,60505 + ,242020 + ,315009 + ,10 + ,60945 + ,609450 + ,318903 + ,1 + ,60720 + ,60720 + ,314887 + ,4 + ,60720 + ,242880 + ,314913 + ,0 + ,60720 + ,0 + ,315380 + ,8 + ,58990 + ,471920 + ,325506 + ,0 + ,60720 + ,0 + ,315380 + ,11 + ,56750 + ,624250 + ,298568 + ,4 + ,60894 + ,243576 + ,315834 + ,15 + ,63346 + ,950190 + ,329784 + ,9 + ,56535 + ,508815 + ,312878 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,60835 + ,425845 + ,314987 + ,2 + ,60720 + ,121440 + ,325249 + ,0 + ,61016 + ,0 + ,315877 + ,7 + ,58650 + ,410550 + ,291650 + ,46 + ,60438 + ,2780148 + ,305959 + ,5 + ,60720 + ,303600 + ,315380 + ,7 + ,58625 + ,410375 + ,297765 + ,2 + ,60938 + ,121876 + ,315245 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,2 + ,61490 + ,122980 + ,315236 + ,5 + ,60845 + ,304225 + ,336425 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,60830 + ,425810 + ,306268 + ,24 + ,63261 + ,1518264 + ,302187 + ,1 + ,60720 + ,60720 + ,314882 + ,0 + ,60720 + ,0 + ,315380 + ,18 + ,45689 + ,822402 + ,382712 + ,55 + ,60720 + ,3339600 + ,341570 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,3 + ,61564 + ,184692 + ,312412 + ,0 + ,60720 + ,0 + ,315380 + ,9 + ,61938 + ,557442 + ,309596 + ,0 + ,60720 + ,0 + ,315380 + ,8 + ,60951 + ,487608 + ,315547 + ,113 + ,60720 + ,6861360 + ,313267 + ,0 + ,60745 + ,0 + ,316176 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,19 + ,71642 + ,1361198 + ,359335 + ,11 + ,71641 + ,788051 + ,330068 + ,25 + ,55792 + ,1394800 + ,314289 + ,16 + ,71873 + ,1149968 + ,297413 + ,5 + ,62555 + ,312775 + ,314806 + ,11 + ,60370 + ,664070 + ,333210 + ,23 + ,64873 + ,1492079 + ,352108 + ,6 + ,62041 + ,372246 + ,313332 + ,5 + ,65745 + ,328725 + ,291787 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,59500 + ,416500 + ,318745 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,61630 + ,431410 + ,315366 + ,0 + ,60720 + ,0 + ,315380 + ,3 + ,60890 + ,182670 + ,315688 + ,0 + ,60720 + ,0 + ,315380 + ,89 + ,113521 + ,10103369 + ,409642 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,19 + ,80045 + ,1520855 + ,269587 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,12 + ,50804 + ,609648 + ,300962 + ,12 + ,87390 + ,1048680 + ,325479 + ,5 + ,61656 + ,308280 + ,316155 + ,2 + ,65688 + ,131376 + ,318574 + ,0 + ,60720 + ,0 + ,315380 + ,26 + ,48522 + ,1261572 + ,343613 + ,3 + ,60720 + ,182160 + ,306948 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,11 + ,57640 + ,634040 + ,330059 + ,10 + ,61977 + ,619770 + ,288985 + ,5 + ,62620 + ,313100 + ,304485 + ,2 + ,60720 + ,121440 + ,315380 + ,6 + ,60831 + ,364986 + ,315688 + ,7 + ,60646 + ,424522 + ,317736 + ,2 + ,60720 + ,121440 + ,315380 + ,28 + ,56225 + ,1574300 + ,322331 + ,3 + ,60510 + ,181530 + ,296656 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,60698 + ,60698 + ,315354 + ,20 + ,60720 + ,1214400 + ,312161 + ,1 + ,60805 + ,60805 + ,315576 + ,22 + ,61404 + ,1350888 + ,314922 + ,9 + ,60720 + ,546480 + ,314551 + ,0 + ,60720 + ,0 + ,315380 + ,2 + ,65276 + ,130552 + ,312339 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,63915 + ,447405 + ,298700 + ,9 + ,60720 + ,546480 + ,321376 + ,0 + ,60720 + ,0 + ,315380 + ,13 + ,61686 + ,801918 + ,303230 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60743 + ,0 + ,315487 + ,0 + ,60720 + ,0 + ,315380 + ,6 + ,60349 + ,362094 + ,315793 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,3 + ,61360 + ,184080 + ,312887 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,59818 + ,418726 + ,315637 + ,2 + ,72680 + ,145360 + ,324385 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,15 + ,61808 + ,927120 + ,308989 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,9 + ,53110 + ,477990 + ,296702 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,64245 + ,64245 + ,307322 + ,38 + ,73007 + ,2774266 + ,304376 + ,57 + ,82732 + ,4715724 + ,253588 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,54820 + ,383740 + ,309560 + ,26 + ,47705 + ,1240330 + ,298466 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,13 + ,72835 + ,946855 + ,343929 + ,10 + ,58856 + ,588560 + ,331955 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,9 + ,77655 + ,698895 + ,381180 + ,0 + ,60720 + ,0 + ,315380 + ,26 + ,69817 + ,1815242 + ,331420 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,19 + ,60798 + ,1155162 + ,310201 + ,0 + ,60720 + ,0 + ,315380 + ,12 + ,62452 + ,749424 + ,320016 + ,23 + ,64175 + ,1476025 + ,320398 + ,0 + ,60720 + ,0 + ,315380 + ,29 + ,67440 + ,1955760 + ,291841 + ,8 + ,68136 + ,545088 + ,310670 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,26 + ,56726 + ,1474876 + ,313491 + ,0 + ,60720 + ,0 + ,315380 + ,9 + ,70811 + ,637299 + ,331323 + ,0 + ,60720 + ,0 + ,315380 + ,5 + ,60720 + ,303600 + ,319210 + ,3 + ,62045 + ,186135 + ,318098 + ,0 + ,60720 + ,0 + ,315380 + ,13 + ,54323 + ,706199 + ,292754 + ,0 + ,60720 + ,0 + ,315380 + ,12 + ,62841 + ,754092 + ,325176 + ,19 + ,81125 + ,1541375 + ,365959 + ,0 + ,60720 + ,0 + ,315380 + ,10 + ,59506 + ,595060 + ,302409 + ,9 + ,59365 + ,534285 + ,340968 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,9 + ,60720 + ,546480 + ,315380 + ,4 + ,60798 + ,243192 + ,313164 + ,1 + ,58790 + ,58790 + ,301164 + ,1 + ,60720 + ,60720 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,14 + ,61808 + ,865312 + ,344425 + ,12 + ,60735 + ,728820 + ,315394 + ,0 + ,60720 + ,0 + ,315380 + ,19 + ,64016 + ,1216304 + ,316647 + ,17 + ,54683 + ,929611 + ,309836 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,32 + ,87192 + ,2790144 + ,346611 + ,0 + ,60720 + ,0 + ,315380 + ,14 + ,64107 + ,897498 + ,322031 + ,8 + ,60761 + ,486088 + ,315656 + ,4 + ,65990 + ,263960 + ,339445 + ,0 + ,59988 + ,0 + ,314964 + ,20 + ,61167 + ,1223340 + ,297141 + ,5 + ,60719 + ,303595 + ,315372 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,60722 + ,60722 + ,312502 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,4 + ,60379 + ,241516 + ,313729 + ,1 + ,60727 + ,60727 + ,315388 + ,4 + ,60720 + ,242880 + ,315371 + ,20 + ,60925 + ,1218500 + ,296139 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,60896 + ,60896 + ,313880 + ,10 + ,59734 + ,597340 + ,317698 + ,12 + ,62969 + ,755628 + ,295580 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,13 + ,60720 + ,789360 + ,308256 + ,0 + ,60720 + ,0 + ,315380 + ,3 + ,59118 + ,177354 + ,303677 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,10 + ,60720 + ,607200 + ,319369 + ,3 + ,58598 + ,175794 + ,318690 + ,7 + ,61124 + ,427868 + ,314049 + ,10 + ,59595 + ,595950 + ,325699 + ,1 + ,62065 + ,62065 + ,314210 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,15 + ,78780 + ,1181700 + ,322378 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,4 + ,60722 + ,242888 + ,315398 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,28 + ,61600 + ,1724800 + ,308336 + ,9 + ,59635 + ,536715 + ,316386 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,60720 + ,425040 + ,315553 + ,0 + ,60720 + ,0 + ,315380 + ,7 + ,59781 + ,418467 + ,323361 + ,7 + ,76644 + ,536508 + ,336639 + ,3 + ,64820 + ,194460 + ,307424 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,11 + ,56178 + ,617958 + ,295370 + ,7 + ,60436 + ,423052 + ,322340 + ,10 + ,60720 + ,607200 + ,319864 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,60720 + ,0 + ,315380 + ,18 + ,73433 + ,1321794 + ,317291 + ,14 + ,41477 + ,580678 + ,280398 + ,0 + ,60720 + ,0 + ,315380 + ,12 + ,62700 + ,752400 + ,317330 + ,29 + ,67804 + ,1966316 + ,238125 + ,3 + ,59661 + ,178983 + ,327071 + ,6 + ,58620 + ,351720 + ,309038 + ,3 + ,60398 + ,181194 + ,314210 + ,8 + ,58580 + ,468640 + ,307930 + ,10 + ,62710 + ,627100 + ,322327 + ,6 + ,59325 + ,355950 + ,292136 + ,8 + ,60950 + ,487600 + ,263276 + ,6 + ,68060 + ,408360 + ,367655 + ,9 + ,83620 + ,752580 + ,283910 + ,8 + ,58456 + ,467648 + ,283587 + ,26 + ,52811 + ,1373086 + ,243650 + ,239 + ,121173 + ,28960347 + ,438493 + ,7 + ,63870 + ,447090 + ,296261 + ,41 + ,21001 + ,861041 + ,230621 + ,3 + ,70415 + ,211245 + ,304252 + ,8 + ,64230 + ,513840 + ,333505 + ,6 + ,59190 + ,355140 + ,296919 + ,21 + ,69351 + ,1456371 + ,278990 + ,7 + ,64270 + ,449890 + ,276898 + ,11 + ,70694 + ,777634 + ,327007 + ,11 + ,68005 + ,748055 + ,317046 + ,12 + ,58930 + ,707160 + ,304555 + ,9 + ,58320 + ,524880 + ,298096 + ,3 + ,69980 + ,209940 + ,231861 + ,57 + ,69863 + ,3982191 + ,309422 + ,21 + ,63255 + ,1328355 + ,286963 + ,15 + ,57320 + ,859800 + ,269753 + ,32 + ,75230 + ,2407360 + ,448243 + ,11 + ,79420 + ,873620 + ,165404 + ,2 + ,73490 + ,146980 + ,204325 + ,23 + ,35250 + ,810750 + ,407159 + ,20 + ,62285 + ,1245700 + ,290476 + ,24 + ,69206 + ,1660944 + ,275311 + ,1 + ,65920 + ,65920 + ,246541 + ,1 + ,69770 + ,69770 + ,253468 + ,74 + ,72683 + ,5378542 + ,240897 + ,68 + ,-14545 + ,-989060 + ,-83265 + ,20 + ,55830 + ,1116600 + ,-42143 + ,20 + ,55174 + ,1103480 + ,272713 + ,82 + ,67038 + ,5497116 + ,215362 + ,21 + ,51252 + ,1076292 + ,42754 + ,244 + ,157278 + ,38375832 + ,306275 + ,32 + ,79510 + ,2544320 + ,253537 + ,86 + ,77440 + ,6659840 + ,372631 + ,69 + ,27284 + ,1882596 + ,-7170) + ,dim=c(4 + ,431) + ,dimnames=list(c('Orders' + ,'Dividends' + ,'O_D' + ,'Wealth ') + ,1:431)) > y <- array(NA,dim=c(4,431),dimnames=list(c('Orders','Dividends','O_D','Wealth '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '4' > #'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 Wealth\r\r Orders Dividends O_D 1 6282154 807 213118 171986226 2 4321023 444 81767 36304548 3 4111912 412 153198 63117576 4 223193 428 -26007 -11130996 5 1491348 315 126942 39986730 6 1629616 168 157214 26411952 7 1398893 263 129352 34019576 8 1926517 267 234817 62696139 9 983660 228 60448 13782144 10 1443586 129 47818 6168522 11 1073089 104 245546 25536784 12 984885 122 48020 5858440 13 1405225 393 -1710 -672030 14 227132 190 32648 6203120 15 929118 280 95350 26698000 16 1071292 63 151352 9535176 17 638830 102 288170 29393340 18 856956 265 114337 30299305 19 992426 234 37884 8864856 20 444477 277 122844 34027788 21 857217 73 82340 6010820 22 711969 67 79801 5346667 23 702380 103 165548 17051444 24 358589 290 116384 33751360 25 297978 83 134028 11124324 26 585715 56 63838 3574928 27 657954 236 74996 17699056 28 209458 73 31080 2268840 29 786690 34 32168 1093712 30 439798 139 49857 6930123 31 688779 26 87161 2266186 32 574339 70 106113 7427910 33 741409 40 80570 3222800 34 597793 42 102129 4289418 35 644190 12 301670 3620040 36 377934 211 102313 21588043 37 640273 74 88577 6554698 38 697458 80 112477 8998160 39 550608 83 191778 15917574 40 207393 131 79804 10454324 41 301607 203 128294 26043682 42 345783 56 96448 5401088 43 501749 89 93811 8349179 44 379983 88 117520 10341760 45 387475 39 69159 2697201 46 377305 25 101792 2544800 47 370837 49 210568 10317832 48 430866 149 136996 20412404 49 469107 58 121920 7071360 50 194493 41 76403 3132523 51 530670 90 108094 9728460 52 518365 136 134759 18327224 53 491303 97 188873 18320681 54 527021 63 146216 9211608 55 233773 114 156608 17853312 56 405972 77 61348 4723796 57 652925 6 50350 302100 58 446211 47 87720 4122840 59 341340 51 99489 5073939 60 387699 85 87419 7430615 61 493408 43 94355 4057265 62 146494 32 60326 1930432 63 414462 25 94670 2366750 64 364304 77 82425 6346725 65 355178 54 59017 3186918 66 357760 251 90829 22798079 67 261216 15 80791 1211865 68 397144 44 100423 4418612 69 374943 73 131116 9571468 70 424898 85 100269 8522865 71 202055 49 27330 1339170 72 378525 38 39039 1483482 73 310768 35 106885 3740975 74 325738 9 79285 713565 75 394510 34 118881 4041954 76 247060 20 77623 1552460 77 368078 29 114768 3328272 78 236761 11 74015 814165 79 312378 52 69465 3612180 80 339836 13 117869 1532297 81 347385 29 60982 1768478 82 426280 66 90131 5948646 83 352850 33 138971 4586043 84 301881 15 39625 594375 85 377516 15 102725 1540875 86 357312 68 64239 4368252 87 458343 100 90262 9026200 88 354228 13 103960 1351480 89 308636 45 106611 4797495 90 386212 14 103345 1446830 91 393343 36 95551 3439836 92 378509 40 82903 3316120 93 452469 68 63593 4324324 94 364839 29 126910 3680390 95 358649 43 37527 1613661 96 376641 30 60247 1807410 97 429112 9 112995 1016955 98 330546 22 70184 1544048 99 403560 19 130140 2472660 100 317892 9 73221 658989 101 307528 31 76114 2359534 102 235133 19 90534 1720146 103 299243 55 108479 5966345 104 314073 8 113761 910088 105 368186 28 68696 1923488 106 269661 29 71561 2075269 107 125390 48 59831 2871888 108 510834 16 97890 1566240 109 321896 47 101481 4769607 110 249898 20 72954 1459080 111 408881 22 67939 1494658 112 158492 33 48022 1584726 113 292154 44 86111 3788884 114 289513 13 74020 962260 115 378049 6 57530 345180 116 343466 35 56364 1972740 117 332743 8 84990 679920 118 442882 17 88590 1506030 119 214215 11 77200 849200 120 315688 21 61262 1286502 121 375195 92 110309 10148428 122 334280 12 67000 804000 123 355864 112 93099 10427088 124 480382 25 107577 2689425 125 353058 17 62920 1069640 126 217193 23 75832 1744136 127 315380 0 60720 0 128 314533 10 60793 607930 129 318056 23 57935 1332505 130 315380 0 60720 0 131 314353 7 60630 424410 132 369448 25 55637 1390925 133 315380 0 60720 0 134 312846 20 60887 1217740 135 312075 4 60720 242880 136 315009 4 60505 242020 137 318903 10 60945 609450 138 314887 1 60720 60720 139 314913 4 60720 242880 140 315380 0 60720 0 141 325506 8 58990 471920 142 315380 0 60720 0 143 298568 11 56750 624250 144 315834 4 60894 243576 145 329784 15 63346 950190 146 312878 9 56535 508815 147 315380 0 60720 0 148 314987 7 60835 425845 149 325249 2 60720 121440 150 315877 0 61016 0 151 291650 7 58650 410550 152 305959 46 60438 2780148 153 315380 5 60720 303600 154 297765 7 58625 410375 155 315245 2 60938 121876 156 315380 0 60720 0 157 315380 0 60720 0 158 315236 2 61490 122980 159 336425 5 60845 304225 160 315380 0 60720 0 161 315380 0 60720 0 162 315380 0 60720 0 163 315380 0 60720 0 164 306268 7 60830 425810 165 302187 24 63261 1518264 166 314882 1 60720 60720 167 315380 0 60720 0 168 382712 18 45689 822402 169 341570 55 60720 3339600 170 315380 0 60720 0 171 315380 0 60720 0 172 312412 3 61564 184692 173 315380 0 60720 0 174 309596 9 61938 557442 175 315380 0 60720 0 176 315547 8 60951 487608 177 313267 113 60720 6861360 178 316176 0 60745 0 179 315380 0 60720 0 180 315380 0 60720 0 181 359335 19 71642 1361198 182 330068 11 71641 788051 183 314289 25 55792 1394800 184 297413 16 71873 1149968 185 314806 5 62555 312775 186 333210 11 60370 664070 187 352108 23 64873 1492079 188 313332 6 62041 372246 189 291787 5 65745 328725 190 315380 0 60720 0 191 318745 7 59500 416500 192 315380 0 60720 0 193 315366 7 61630 431410 194 315380 0 60720 0 195 315688 3 60890 182670 196 315380 0 60720 0 197 409642 89 113521 10103369 198 315380 0 60720 0 199 315380 0 60720 0 200 269587 19 80045 1520855 201 315380 0 60720 0 202 315380 0 60720 0 203 315380 0 60720 0 204 300962 12 50804 609648 205 325479 12 87390 1048680 206 316155 5 61656 308280 207 318574 2 65688 131376 208 315380 0 60720 0 209 343613 26 48522 1261572 210 306948 3 60720 182160 211 315380 0 60720 0 212 315380 0 60720 0 213 330059 11 57640 634040 214 288985 10 61977 619770 215 304485 5 62620 313100 216 315380 2 60720 121440 217 315688 6 60831 364986 218 317736 7 60646 424522 219 315380 2 60720 121440 220 322331 28 56225 1574300 221 296656 3 60510 181530 222 315380 0 60720 0 223 315354 1 60698 60698 224 312161 20 60720 1214400 225 315576 1 60805 60805 226 314922 22 61404 1350888 227 314551 9 60720 546480 228 315380 0 60720 0 229 312339 2 65276 130552 230 315380 0 60720 0 231 298700 7 63915 447405 232 321376 9 60720 546480 233 315380 0 60720 0 234 303230 13 61686 801918 235 315380 0 60720 0 236 315487 0 60743 0 237 315380 0 60720 0 238 315793 6 60349 362094 239 315380 0 60720 0 240 315380 0 60720 0 241 315380 0 60720 0 242 312887 3 61360 184080 243 315380 0 60720 0 244 315637 7 59818 418726 245 324385 2 72680 145360 246 315380 0 60720 0 247 315380 0 60720 0 248 308989 15 61808 927120 249 315380 0 60720 0 250 315380 0 60720 0 251 296702 9 53110 477990 252 315380 0 60720 0 253 307322 1 64245 64245 254 304376 38 73007 2774266 255 253588 57 82732 4715724 256 315380 0 60720 0 257 309560 7 54820 383740 258 298466 26 47705 1240330 259 315380 0 60720 0 260 315380 0 60720 0 261 315380 0 60720 0 262 315380 0 60720 0 263 343929 13 72835 946855 264 331955 10 58856 588560 265 315380 0 60720 0 266 315380 0 60720 0 267 315380 0 60720 0 268 381180 9 77655 698895 269 315380 0 60720 0 270 331420 26 69817 1815242 271 315380 0 60720 0 272 315380 0 60720 0 273 315380 0 60720 0 274 310201 19 60798 1155162 275 315380 0 60720 0 276 320016 12 62452 749424 277 320398 23 64175 1476025 278 315380 0 60720 0 279 291841 29 67440 1955760 280 310670 8 68136 545088 281 315380 0 60720 0 282 315380 0 60720 0 283 313491 26 56726 1474876 284 315380 0 60720 0 285 331323 9 70811 637299 286 315380 0 60720 0 287 319210 5 60720 303600 288 318098 3 62045 186135 289 315380 0 60720 0 290 292754 13 54323 706199 291 315380 0 60720 0 292 325176 12 62841 754092 293 365959 19 81125 1541375 294 315380 0 60720 0 295 302409 10 59506 595060 296 340968 9 59365 534285 297 315380 0 60720 0 298 315380 0 60720 0 299 315380 0 60720 0 300 315380 9 60720 546480 301 313164 4 60798 243192 302 301164 1 58790 58790 303 315380 1 60720 60720 304 315380 0 60720 0 305 344425 14 61808 865312 306 315394 12 60735 728820 307 315380 0 60720 0 308 316647 19 64016 1216304 309 309836 17 54683 929611 310 315380 0 60720 0 311 315380 0 60720 0 312 346611 32 87192 2790144 313 315380 0 60720 0 314 322031 14 64107 897498 315 315656 8 60761 486088 316 339445 4 65990 263960 317 314964 0 59988 0 318 297141 20 61167 1223340 319 315372 5 60719 303595 320 315380 0 60720 0 321 315380 0 60720 0 322 315380 0 60720 0 323 315380 0 60720 0 324 315380 0 60720 0 325 315380 0 60720 0 326 312502 1 60722 60722 327 315380 0 60720 0 328 315380 0 60720 0 329 315380 0 60720 0 330 315380 0 60720 0 331 315380 0 60720 0 332 315380 0 60720 0 333 315380 0 60720 0 334 313729 4 60379 241516 335 315388 1 60727 60727 336 315371 4 60720 242880 337 296139 20 60925 1218500 338 315380 0 60720 0 339 313880 1 60896 60896 340 317698 10 59734 597340 341 295580 12 62969 755628 342 315380 0 60720 0 343 315380 0 60720 0 344 315380 0 60720 0 345 308256 13 60720 789360 346 315380 0 60720 0 347 303677 3 59118 177354 348 315380 0 60720 0 349 315380 0 60720 0 350 319369 10 60720 607200 351 318690 3 58598 175794 352 314049 7 61124 427868 353 325699 10 59595 595950 354 314210 1 62065 62065 355 315380 0 60720 0 356 315380 0 60720 0 357 322378 15 78780 1181700 358 315380 0 60720 0 359 315380 0 60720 0 360 315380 0 60720 0 361 315398 4 60722 242888 362 315380 0 60720 0 363 315380 0 60720 0 364 308336 28 61600 1724800 365 316386 9 59635 536715 366 315380 0 60720 0 367 315380 0 60720 0 368 315380 0 60720 0 369 315380 0 60720 0 370 315553 7 60720 425040 371 315380 0 60720 0 372 323361 7 59781 418467 373 336639 7 76644 536508 374 307424 3 64820 194460 375 315380 0 60720 0 376 315380 0 60720 0 377 295370 11 56178 617958 378 322340 7 60436 423052 379 319864 10 60720 607200 380 315380 0 60720 0 381 315380 0 60720 0 382 317291 18 73433 1321794 383 280398 14 41477 580678 384 315380 0 60720 0 385 317330 12 62700 752400 386 238125 29 67804 1966316 387 327071 3 59661 178983 388 309038 6 58620 351720 389 314210 3 60398 181194 390 307930 8 58580 468640 391 322327 10 62710 627100 392 292136 6 59325 355950 393 263276 8 60950 487600 394 367655 6 68060 408360 395 283910 9 83620 752580 396 283587 8 58456 467648 397 243650 26 52811 1373086 398 438493 239 121173 28960347 399 296261 7 63870 447090 400 230621 41 21001 861041 401 304252 3 70415 211245 402 333505 8 64230 513840 403 296919 6 59190 355140 404 278990 21 69351 1456371 405 276898 7 64270 449890 406 327007 11 70694 777634 407 317046 11 68005 748055 408 304555 12 58930 707160 409 298096 9 58320 524880 410 231861 3 69980 209940 411 309422 57 69863 3982191 412 286963 21 63255 1328355 413 269753 15 57320 859800 414 448243 32 75230 2407360 415 165404 11 79420 873620 416 204325 2 73490 146980 417 407159 23 35250 810750 418 290476 20 62285 1245700 419 275311 24 69206 1660944 420 246541 1 65920 65920 421 253468 1 69770 69770 422 240897 74 72683 5378542 423 -83265 68 -14545 -989060 424 -42143 20 55830 1116600 425 272713 20 55174 1103480 426 215362 82 67038 5497116 427 42754 21 51252 1076292 428 306275 244 157278 38375832 429 253537 32 79510 2544320 430 372631 86 77440 6659840 431 -7170 69 27284 1882596 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Orders Dividends O_D 3.917e+05 8.432e+02 -1.823e+00 3.057e-02 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -1177595 -23719 23639 34415 2594168 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 3.917e+05 3.341e+04 11.725 < 2e-16 *** Orders 8.432e+02 2.857e+02 2.951 0.00334 ** Dividends -1.823e+00 4.490e-01 -4.061 5.81e-05 *** O_D 3.057e-02 2.202e-03 13.883 < 2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 234800 on 427 degrees of freedom Multiple R-squared: 0.7047, Adjusted R-squared: 0.7026 F-statistic: 339.7 on 3 and 427 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 2.412205e-13 1.206102e-13 [2,] 1.0000000 1.766274e-72 8.831368e-73 [3,] 1.0000000 4.413854e-73 2.206927e-73 [4,] 1.0000000 2.307902e-87 1.153951e-87 [5,] 1.0000000 8.967041e-91 4.483520e-91 [6,] 1.0000000 1.097964e-93 5.489818e-94 [7,] 1.0000000 6.203572e-95 3.101786e-95 [8,] 1.0000000 2.514953e-100 1.257476e-100 [9,] 1.0000000 4.047570e-106 2.023785e-106 [10,] 1.0000000 9.648272e-114 4.824136e-114 [11,] 1.0000000 3.830036e-115 1.915018e-115 [12,] 1.0000000 2.375063e-122 1.187532e-122 [13,] 1.0000000 3.431091e-127 1.715546e-127 [14,] 1.0000000 5.530630e-137 2.765315e-137 [15,] 1.0000000 2.300575e-143 1.150288e-143 [16,] 1.0000000 2.623895e-146 1.311948e-146 [17,] 1.0000000 2.403232e-147 1.201616e-147 [18,] 1.0000000 4.039002e-157 2.019501e-157 [19,] 1.0000000 4.866210e-157 2.433105e-157 [20,] 1.0000000 3.271720e-158 1.635860e-158 [21,] 1.0000000 2.096207e-161 1.048103e-161 [22,] 1.0000000 5.982313e-161 2.991157e-161 [23,] 1.0000000 5.428947e-171 2.714473e-171 [24,] 1.0000000 3.620383e-171 1.810192e-171 [25,] 1.0000000 8.639791e-176 4.319896e-176 [26,] 1.0000000 1.617980e-176 8.089899e-177 [27,] 1.0000000 5.849597e-184 2.924798e-184 [28,] 1.0000000 8.062963e-186 4.031482e-186 [29,] 1.0000000 1.128703e-186 5.643516e-187 [30,] 1.0000000 3.965101e-189 1.982550e-189 [31,] 1.0000000 5.263843e-193 2.631921e-193 [32,] 1.0000000 3.451767e-198 1.725883e-198 [33,] 1.0000000 1.813116e-197 9.065582e-198 [34,] 1.0000000 1.141767e-198 5.708837e-199 [35,] 1.0000000 2.206972e-202 1.103486e-202 [36,] 1.0000000 2.400054e-201 1.200027e-201 [37,] 1.0000000 7.448426e-202 3.724213e-202 [38,] 1.0000000 5.198009e-201 2.599005e-201 [39,] 1.0000000 2.803767e-200 1.401883e-200 [40,] 1.0000000 2.917671e-199 1.458836e-199 [41,] 1.0000000 9.883787e-200 4.941893e-200 [42,] 1.0000000 5.845799e-200 2.922900e-200 [43,] 1.0000000 2.838734e-199 1.419367e-199 [44,] 1.0000000 1.771905e-199 8.859527e-200 [45,] 1.0000000 3.040286e-200 1.520143e-200 [46,] 1.0000000 1.067648e-200 5.338238e-201 [47,] 1.0000000 6.883826e-200 3.441913e-200 [48,] 1.0000000 1.191128e-199 5.955641e-200 [49,] 1.0000000 8.534312e-203 4.267156e-203 [50,] 1.0000000 9.299268e-203 4.649634e-203 [51,] 1.0000000 1.796220e-212 8.981098e-213 [52,] 1.0000000 2.208819e-212 1.104409e-212 [53,] 1.0000000 2.421322e-211 1.210661e-211 [54,] 1.0000000 1.167255e-210 5.836275e-211 [55,] 1.0000000 1.829768e-211 9.148838e-212 [56,] 1.0000000 1.440850e-212 7.204252e-213 [57,] 1.0000000 6.281024e-212 3.140512e-212 [58,] 1.0000000 4.107223e-211 2.053611e-211 [59,] 1.0000000 1.826123e-210 9.130615e-211 [60,] 1.0000000 4.783515e-213 2.391758e-213 [61,] 1.0000000 1.836706e-212 9.183530e-213 [62,] 1.0000000 1.414260e-211 7.071301e-212 [63,] 1.0000000 1.244691e-210 6.223456e-211 [64,] 1.0000000 3.003511e-210 1.501756e-210 [65,] 1.0000000 2.108897e-209 1.054449e-209 [66,] 1.0000000 1.150879e-209 5.754396e-210 [67,] 1.0000000 6.600556e-209 3.300278e-209 [68,] 1.0000000 6.424888e-208 3.212444e-208 [69,] 1.0000000 6.369586e-207 3.184793e-207 [70,] 1.0000000 1.767582e-206 8.837912e-207 [71,] 1.0000000 1.688796e-205 8.443982e-206 [72,] 1.0000000 3.262399e-205 1.631200e-205 [73,] 1.0000000 3.246616e-204 1.623308e-204 [74,] 1.0000000 1.513042e-203 7.565208e-204 [75,] 1.0000000 9.343020e-203 4.671510e-203 [76,] 1.0000000 1.353297e-202 6.766486e-203 [77,] 1.0000000 4.143209e-202 2.071605e-202 [78,] 1.0000000 3.422896e-201 1.711448e-201 [79,] 1.0000000 2.891949e-200 1.445975e-200 [80,] 1.0000000 7.406198e-200 3.703099e-200 [81,] 1.0000000 5.486298e-201 2.743149e-201 [82,] 1.0000000 4.335011e-200 2.167505e-200 [83,] 1.0000000 2.029235e-199 1.014617e-199 [84,] 1.0000000 1.610871e-198 8.054354e-199 [85,] 1.0000000 1.070382e-197 5.351908e-198 [86,] 1.0000000 5.763223e-197 2.881612e-197 [87,] 1.0000000 2.397773e-199 1.198886e-199 [88,] 1.0000000 1.464572e-198 7.322861e-199 [89,] 1.0000000 3.359115e-199 1.679558e-199 [90,] 1.0000000 4.821630e-199 2.410815e-199 [91,] 1.0000000 2.432947e-198 1.216473e-198 [92,] 1.0000000 2.242750e-197 1.121375e-197 [93,] 1.0000000 1.623685e-196 8.118426e-197 [94,] 1.0000000 1.529186e-195 7.645930e-196 [95,] 1.0000000 1.501032e-194 7.505160e-195 [96,] 1.0000000 5.254379e-195 2.627190e-195 [97,] 1.0000000 1.278745e-194 6.393724e-195 [98,] 1.0000000 8.344600e-195 4.172300e-195 [99,] 1.0000000 3.305689e-194 1.652844e-194 [100,] 1.0000000 2.015257e-193 1.007629e-193 [101,] 1.0000000 5.848798e-195 2.924399e-195 [102,] 1.0000000 4.088605e-196 2.044303e-196 [103,] 1.0000000 2.387243e-195 1.193622e-195 [104,] 1.0000000 5.714812e-195 2.857406e-195 [105,] 1.0000000 4.196784e-195 2.098392e-195 [106,] 1.0000000 2.103089e-195 1.051545e-195 [107,] 1.0000000 1.220652e-194 6.103260e-195 [108,] 1.0000000 7.554101e-194 3.777050e-194 [109,] 1.0000000 1.408838e-193 7.044192e-194 [110,] 1.0000000 4.118940e-193 2.059470e-193 [111,] 1.0000000 3.141841e-192 1.570920e-192 [112,] 1.0000000 3.576496e-192 1.788248e-192 [113,] 1.0000000 4.391048e-193 2.195524e-193 [114,] 1.0000000 4.032869e-192 2.016435e-192 [115,] 1.0000000 2.771426e-191 1.385713e-191 [116,] 1.0000000 2.451994e-190 1.225997e-190 [117,] 1.0000000 6.462261e-190 3.231130e-190 [118,] 1.0000000 4.241149e-190 2.120575e-190 [119,] 1.0000000 2.001628e-189 1.000814e-189 [120,] 1.0000000 7.152407e-190 3.576203e-190 [121,] 1.0000000 6.761905e-189 3.380952e-189 [122,] 1.0000000 6.418398e-188 3.209199e-188 [123,] 1.0000000 4.910069e-187 2.455034e-187 [124,] 1.0000000 4.584165e-186 2.292083e-186 [125,] 1.0000000 4.307217e-185 2.153609e-185 [126,] 1.0000000 4.452139e-185 2.226069e-185 [127,] 1.0000000 4.136262e-184 2.068131e-184 [128,] 1.0000000 3.686498e-183 1.843249e-183 [129,] 1.0000000 3.444550e-182 1.722275e-182 [130,] 1.0000000 3.175257e-181 1.587628e-181 [131,] 1.0000000 2.815349e-180 1.407675e-180 [132,] 1.0000000 2.566842e-179 1.283421e-179 [133,] 1.0000000 2.337149e-178 1.168575e-178 [134,] 1.0000000 2.106308e-177 1.053154e-177 [135,] 1.0000000 1.667822e-176 8.339112e-177 [136,] 1.0000000 1.490812e-175 7.454062e-176 [137,] 1.0000000 1.365502e-174 6.827509e-175 [138,] 1.0000000 1.211944e-173 6.059719e-174 [139,] 1.0000000 9.261667e-173 4.630834e-173 [140,] 1.0000000 7.789612e-172 3.894806e-172 [141,] 1.0000000 6.806988e-171 3.403494e-171 [142,] 1.0000000 5.928339e-170 2.964170e-170 [143,] 1.0000000 4.866655e-169 2.433328e-169 [144,] 1.0000000 4.194526e-168 2.097263e-168 [145,] 1.0000000 3.561164e-167 1.780582e-167 [146,] 1.0000000 2.158148e-166 1.079074e-166 [147,] 1.0000000 1.837582e-165 9.187910e-166 [148,] 1.0000000 1.581189e-164 7.905943e-165 [149,] 1.0000000 1.336549e-163 6.682747e-164 [150,] 1.0000000 1.121648e-162 5.608239e-163 [151,] 1.0000000 9.372090e-162 4.686045e-162 [152,] 1.0000000 7.824196e-161 3.912098e-161 [153,] 1.0000000 5.288853e-160 2.644427e-160 [154,] 1.0000000 4.364169e-159 2.182085e-159 [155,] 1.0000000 3.585304e-158 1.792652e-158 [156,] 1.0000000 2.932445e-157 1.466223e-157 [157,] 1.0000000 2.387858e-156 1.193929e-156 [158,] 1.0000000 1.967429e-155 9.837145e-156 [159,] 1.0000000 1.598483e-154 7.992413e-155 [160,] 1.0000000 1.285524e-153 6.427620e-154 [161,] 1.0000000 1.027301e-152 5.136507e-153 [162,] 1.0000000 2.312843e-153 1.156422e-153 [163,] 1.0000000 3.353002e-153 1.676501e-153 [164,] 1.0000000 2.696874e-152 1.348437e-152 [165,] 1.0000000 2.159534e-151 1.079767e-151 [166,] 1.0000000 1.735157e-150 8.675784e-151 [167,] 1.0000000 1.376993e-149 6.884964e-150 [168,] 1.0000000 1.102464e-148 5.512318e-149 [169,] 1.0000000 8.669010e-148 4.334505e-148 [170,] 1.0000000 6.700901e-147 3.350450e-147 [171,] 1.0000000 8.465111e-148 4.232556e-148 [172,] 1.0000000 6.674165e-147 3.337083e-147 [173,] 1.0000000 5.245210e-146 2.622605e-146 [174,] 1.0000000 4.103956e-145 2.051978e-145 [175,] 1.0000000 1.898971e-144 9.494854e-145 [176,] 1.0000000 1.459144e-143 7.295721e-144 [177,] 1.0000000 7.552188e-143 3.776094e-143 [178,] 1.0000000 5.285778e-142 2.642889e-142 [179,] 1.0000000 4.076854e-141 2.038427e-141 [180,] 1.0000000 2.376149e-140 1.188075e-140 [181,] 1.0000000 8.011954e-140 4.005977e-140 [182,] 1.0000000 6.123967e-139 3.061983e-139 [183,] 1.0000000 3.955153e-138 1.977577e-138 [184,] 1.0000000 2.984162e-137 1.492081e-137 [185,] 1.0000000 2.125321e-136 1.062660e-136 [186,] 1.0000000 1.589491e-135 7.947455e-136 [187,] 1.0000000 1.174327e-134 5.871633e-135 [188,] 1.0000000 8.704013e-134 4.352006e-134 [189,] 1.0000000 6.408503e-133 3.204251e-133 [190,] 1.0000000 4.707580e-132 2.353790e-132 [191,] 1.0000000 1.263677e-131 6.318384e-132 [192,] 1.0000000 9.241575e-131 4.620787e-131 [193,] 1.0000000 6.728548e-130 3.364274e-130 [194,] 1.0000000 2.026907e-129 1.013453e-129 [195,] 1.0000000 1.468098e-128 7.340491e-129 [196,] 1.0000000 1.058620e-127 5.293101e-128 [197,] 1.0000000 7.599537e-127 3.799769e-127 [198,] 1.0000000 4.930846e-126 2.465423e-126 [199,] 1.0000000 2.918631e-125 1.459316e-125 [200,] 1.0000000 2.059901e-124 1.029950e-124 [201,] 1.0000000 1.448669e-123 7.243345e-124 [202,] 1.0000000 1.017938e-122 5.089692e-123 [203,] 1.0000000 1.100163e-122 5.500815e-123 [204,] 1.0000000 7.765123e-122 3.882561e-122 [205,] 1.0000000 5.428828e-121 2.714414e-121 [206,] 1.0000000 3.778464e-120 1.889232e-120 [207,] 1.0000000 1.918485e-119 9.592423e-120 [208,] 1.0000000 1.267244e-118 6.336221e-119 [209,] 1.0000000 8.688477e-118 4.344239e-118 [210,] 1.0000000 5.939891e-117 2.969945e-117 [211,] 1.0000000 3.992126e-116 1.996063e-116 [212,] 1.0000000 2.620472e-115 1.310236e-115 [213,] 1.0000000 1.767101e-114 8.835504e-115 [214,] 1.0000000 6.117434e-114 3.058717e-114 [215,] 1.0000000 3.994508e-113 1.997254e-113 [216,] 1.0000000 2.669638e-112 1.334819e-112 [217,] 1.0000000 1.775367e-111 8.876837e-112 [218,] 1.0000000 1.042669e-110 5.213345e-111 [219,] 1.0000000 6.872403e-110 3.436201e-110 [220,] 1.0000000 3.825780e-109 1.912890e-109 [221,] 1.0000000 2.433006e-108 1.216503e-108 [222,] 1.0000000 1.583965e-107 7.919827e-108 [223,] 1.0000000 1.015766e-106 5.078831e-107 [224,] 1.0000000 6.553756e-106 3.276878e-106 [225,] 1.0000000 4.103610e-105 2.051805e-105 [226,] 1.0000000 2.428433e-104 1.214216e-104 [227,] 1.0000000 1.546124e-103 7.730620e-104 [228,] 1.0000000 9.839761e-103 4.919881e-103 [229,] 1.0000000 6.206677e-102 3.103338e-102 [230,] 1.0000000 3.896610e-101 1.948305e-101 [231,] 1.0000000 2.435400e-100 1.217700e-100 [232,] 1.0000000 1.484500e-99 7.422501e-100 [233,] 1.0000000 9.193155e-99 4.596577e-99 [234,] 1.0000000 5.666847e-98 2.833424e-98 [235,] 1.0000000 3.477022e-97 1.738511e-97 [236,] 1.0000000 2.130140e-96 1.065070e-96 [237,] 1.0000000 1.294866e-95 6.474332e-96 [238,] 1.0000000 7.594533e-95 3.797267e-95 [239,] 1.0000000 4.437576e-94 2.218788e-94 [240,] 1.0000000 2.660942e-93 1.330471e-93 [241,] 1.0000000 1.588151e-92 7.940756e-93 [242,] 1.0000000 9.203820e-92 4.601910e-92 [243,] 1.0000000 5.441907e-91 2.720954e-91 [244,] 1.0000000 3.202502e-90 1.601251e-90 [245,] 1.0000000 1.837726e-89 9.188630e-90 [246,] 1.0000000 1.071250e-88 5.356250e-89 [247,] 1.0000000 6.083797e-88 3.041898e-88 [248,] 1.0000000 3.395655e-87 1.697828e-87 [249,] 1.0000000 1.087657e-86 5.438283e-87 [250,] 1.0000000 6.241738e-86 3.120869e-86 [251,] 1.0000000 3.355179e-85 1.677590e-85 [252,] 1.0000000 1.193107e-84 5.965533e-85 [253,] 1.0000000 6.767021e-84 3.383510e-84 [254,] 1.0000000 3.819713e-83 1.909856e-83 [255,] 1.0000000 2.145722e-82 1.072861e-82 [256,] 1.0000000 1.199555e-81 5.997776e-82 [257,] 1.0000000 6.196906e-81 3.098453e-81 [258,] 1.0000000 2.685745e-80 1.342872e-80 [259,] 1.0000000 1.482240e-79 7.411200e-80 [260,] 1.0000000 8.140541e-79 4.070270e-79 [261,] 1.0000000 4.449006e-78 2.224503e-78 [262,] 1.0000000 1.672812e-77 8.364060e-78 [263,] 1.0000000 9.074498e-77 4.537249e-77 [264,] 1.0000000 4.175415e-76 2.087708e-76 [265,] 1.0000000 2.245004e-75 1.122502e-75 [266,] 1.0000000 1.201096e-74 6.005480e-75 [267,] 1.0000000 6.394035e-74 3.197018e-74 [268,] 1.0000000 3.124242e-73 1.562121e-73 [269,] 1.0000000 1.647256e-72 8.236279e-73 [270,] 1.0000000 8.123251e-72 4.061625e-72 [271,] 1.0000000 3.658849e-71 1.829425e-71 [272,] 1.0000000 1.902360e-70 9.511799e-71 [273,] 1.0000000 9.761299e-70 4.880650e-70 [274,] 1.0000000 5.026295e-69 2.513148e-69 [275,] 1.0000000 2.574047e-68 1.287024e-68 [276,] 1.0000000 1.311473e-67 6.557367e-68 [277,] 1.0000000 4.802342e-67 2.401171e-67 [278,] 1.0000000 2.426741e-66 1.213370e-66 [279,] 1.0000000 1.202440e-65 6.012198e-66 [280,] 1.0000000 6.014112e-65 3.007056e-65 [281,] 1.0000000 2.900420e-64 1.450210e-64 [282,] 1.0000000 1.424875e-63 7.124374e-64 [283,] 1.0000000 7.016138e-63 3.508069e-63 [284,] 1.0000000 3.360883e-62 1.680441e-62 [285,] 1.0000000 1.637498e-61 8.187488e-62 [286,] 1.0000000 7.226363e-61 3.613181e-61 [287,] 1.0000000 2.900484e-60 1.450242e-60 [288,] 1.0000000 1.393494e-59 6.967468e-60 [289,] 1.0000000 6.651419e-59 3.325710e-59 [290,] 1.0000000 2.310122e-58 1.155061e-58 [291,] 1.0000000 1.094015e-57 5.470076e-58 [292,] 1.0000000 5.152731e-57 2.576365e-57 [293,] 1.0000000 2.413604e-56 1.206802e-56 [294,] 1.0000000 1.083765e-55 5.418825e-56 [295,] 1.0000000 5.001869e-55 2.500934e-55 [296,] 1.0000000 2.324253e-54 1.162127e-54 [297,] 1.0000000 1.062764e-53 5.313822e-54 [298,] 1.0000000 4.842389e-53 2.421194e-53 [299,] 1.0000000 1.527104e-52 7.635522e-53 [300,] 1.0000000 6.504914e-52 3.252457e-52 [301,] 1.0000000 2.922012e-51 1.461006e-51 [302,] 1.0000000 1.208992e-50 6.044962e-51 [303,] 1.0000000 4.504451e-50 2.252225e-50 [304,] 1.0000000 1.991805e-49 9.959023e-50 [305,] 1.0000000 8.756181e-49 4.378091e-49 [306,] 1.0000000 3.726025e-48 1.863013e-48 [307,] 1.0000000 1.619686e-47 8.098431e-48 [308,] 1.0000000 6.513086e-47 3.256543e-47 [309,] 1.0000000 2.706139e-46 1.353069e-46 [310,] 1.0000000 1.055300e-45 5.276500e-46 [311,] 1.0000000 4.473746e-45 2.236873e-45 [312,] 1.0000000 1.853870e-44 9.269348e-45 [313,] 1.0000000 7.647664e-44 3.823832e-44 [314,] 1.0000000 3.190879e-43 1.595440e-43 [315,] 1.0000000 1.323121e-42 6.615607e-43 [316,] 1.0000000 5.452286e-42 2.726143e-42 [317,] 1.0000000 2.232693e-41 1.116346e-41 [318,] 1.0000000 9.085143e-41 4.542571e-41 [319,] 1.0000000 3.673397e-40 1.836698e-40 [320,] 1.0000000 1.482302e-39 7.411512e-40 [321,] 1.0000000 5.916782e-39 2.958391e-39 [322,] 1.0000000 2.346429e-38 1.173214e-38 [323,] 1.0000000 9.244446e-38 4.622223e-38 [324,] 1.0000000 3.618133e-37 1.809066e-37 [325,] 1.0000000 1.406682e-36 7.033410e-37 [326,] 1.0000000 5.432406e-36 2.716203e-36 [327,] 1.0000000 2.083771e-35 1.041886e-35 [328,] 1.0000000 7.862928e-35 3.931464e-35 [329,] 1.0000000 2.966948e-34 1.483474e-34 [330,] 1.0000000 1.100472e-33 5.502358e-34 [331,] 1.0000000 4.050514e-33 2.025257e-33 [332,] 1.0000000 1.501189e-32 7.505943e-33 [333,] 1.0000000 5.536218e-32 2.768109e-32 [334,] 1.0000000 1.893921e-31 9.469604e-32 [335,] 1.0000000 6.988680e-31 3.494340e-31 [336,] 1.0000000 2.518117e-30 1.259058e-30 [337,] 1.0000000 9.006482e-30 4.503241e-30 [338,] 1.0000000 3.197436e-29 1.598718e-29 [339,] 1.0000000 1.099070e-28 5.495349e-29 [340,] 1.0000000 3.844277e-28 1.922139e-28 [341,] 1.0000000 1.350689e-27 6.753445e-28 [342,] 1.0000000 4.652514e-27 2.326257e-27 [343,] 1.0000000 1.590083e-26 7.950415e-27 [344,] 1.0000000 5.077499e-26 2.538749e-26 [345,] 1.0000000 1.639629e-25 8.198145e-26 [346,] 1.0000000 5.382420e-25 2.691210e-25 [347,] 1.0000000 1.588019e-24 7.940097e-25 [348,] 1.0000000 5.260168e-24 2.630084e-24 [349,] 1.0000000 1.714120e-23 8.570599e-24 [350,] 1.0000000 5.538428e-23 2.769214e-23 [351,] 1.0000000 1.822020e-22 9.110101e-23 [352,] 1.0000000 5.790287e-22 2.895143e-22 [353,] 1.0000000 1.823969e-21 9.119847e-22 [354,] 1.0000000 5.694449e-21 2.847225e-21 [355,] 1.0000000 1.739242e-20 8.696212e-21 [356,] 1.0000000 5.331895e-20 2.665948e-20 [357,] 1.0000000 1.619365e-19 8.096826e-20 [358,] 1.0000000 4.459122e-19 2.229561e-19 [359,] 1.0000000 1.267156e-18 6.335782e-19 [360,] 1.0000000 3.747410e-18 1.873705e-18 [361,] 1.0000000 1.097285e-17 5.486424e-18 [362,] 1.0000000 3.180664e-17 1.590332e-17 [363,] 1.0000000 9.125245e-17 4.562623e-17 [364,] 1.0000000 2.528501e-16 1.264250e-16 [365,] 1.0000000 7.103848e-16 3.551924e-16 [366,] 1.0000000 1.834875e-15 9.174377e-16 [367,] 1.0000000 5.180395e-15 2.590198e-15 [368,] 1.0000000 1.463573e-14 7.317863e-15 [369,] 1.0000000 3.942440e-14 1.971220e-14 [370,] 1.0000000 1.049551e-13 5.247753e-14 [371,] 1.0000000 2.755193e-13 1.377596e-13 [372,] 1.0000000 6.762735e-13 3.381368e-13 [373,] 1.0000000 1.645781e-12 8.228906e-13 [374,] 1.0000000 4.173207e-12 2.086603e-12 [375,] 1.0000000 1.044065e-11 5.220326e-12 [376,] 1.0000000 2.710694e-11 1.355347e-11 [377,] 1.0000000 5.748747e-11 2.874374e-11 [378,] 1.0000000 1.379712e-10 6.898559e-11 [379,] 1.0000000 3.207677e-10 1.603838e-10 [380,] 1.0000000 7.492963e-10 3.746482e-10 [381,] 1.0000000 1.601327e-09 8.006635e-10 [382,] 1.0000000 3.597136e-09 1.798568e-09 [383,] 1.0000000 8.000891e-09 4.000446e-09 [384,] 1.0000000 1.736538e-08 8.682691e-09 [385,] 1.0000000 3.658660e-08 1.829330e-08 [386,] 1.0000000 8.269409e-08 4.134704e-08 [387,] 0.9999999 1.948187e-07 9.740936e-08 [388,] 0.9999998 3.321823e-07 1.660911e-07 [389,] 0.9999996 7.145376e-07 3.572688e-07 [390,] 0.9999992 1.553908e-06 7.769542e-07 [391,] 0.9999983 3.415269e-06 1.707634e-06 [392,] 0.9999989 2.211253e-06 1.105626e-06 [393,] 0.9999975 4.958514e-06 2.479257e-06 [394,] 0.9999969 6.285862e-06 3.142931e-06 [395,] 0.9999929 1.423416e-05 7.117079e-06 [396,] 0.9999865 2.695334e-05 1.347667e-05 [397,] 0.9999730 5.406222e-05 2.703111e-05 [398,] 0.9999415 1.170785e-04 5.853926e-05 [399,] 0.9998766 2.467473e-04 1.233736e-04 [400,] 0.9997595 4.809721e-04 2.404861e-04 [401,] 0.9995423 9.153458e-04 4.576729e-04 [402,] 0.9992143 1.571397e-03 7.856986e-04 [403,] 0.9986712 2.657539e-03 1.328770e-03 [404,] 0.9974800 5.040079e-03 2.520040e-03 [405,] 0.9956407 8.718587e-03 4.359293e-03 [406,] 0.9923942 1.521159e-02 7.605797e-03 [407,] 0.9872296 2.554074e-02 1.277037e-02 [408,] 0.9889947 2.201058e-02 1.100529e-02 [409,] 0.9843607 3.127861e-02 1.563931e-02 [410,] 0.9739277 5.214466e-02 2.607233e-02 [411,] 0.9976557 4.688672e-03 2.344336e-03 [412,] 0.9963847 7.230510e-03 3.615255e-03 [413,] 0.9920132 1.597370e-02 7.986849e-03 [414,] 0.9844017 3.119656e-02 1.559828e-02 [415,] 0.9707689 5.846214e-02 2.923107e-02 [416,] 0.9348681 1.302639e-01 6.513194e-02 [417,] 0.8855443 2.289113e-01 1.144557e-01 [418,] 0.8880029 2.239942e-01 1.119971e-01 > postscript(file="/var/www/rcomp/tmp/1nb4g1292956049.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/2nb4g1292956049.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/3nb4g1292956049.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/4xk311292956049.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/5xk311292956049.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 340675.6874 2594167.7524 1722550.3597 -236477.4865 -156944.7777 6 7 8 9 10 575467.7008 -18734.1681 -178890.7705 88611.2273 841746.1429 11 12 13 14 15 260708.2843 398795.4080 699614.9668 -454857.4185 -341002.5886 16 17 18 19 20 610942.5762 -212053.1222 -475997.8305 201509.3228 -997071.1842 21 22 23 24 25 370357.3971 245843.6206 4400.2931 -1097247.5545 -259403.5402 26 27 28 29 30 153925.6423 -337062.9511 -256461.0148 391564.1640 -190039.5441 31 32 33 34 35 364818.4903 90029.4346 364381.7852 225777.9029 681742.2419 36 37 38 39 40 -665091.9651 147314.3445 168314.3345 -48020.6220 -468840.2425 41 42 43 44 45 -823521.3136 -82380.0097 -49175.0083 -187787.8009 6554.6514 46 47 48 49 50 72345.5298 6328.8551 -460706.9652 34635.5935 -188214.8335 51 52 53 54 55 -37223.1732 -302577.5598 -197893.8234 67199.5523 -514295.7847 56 57 58 59 60 -83186.7776 338758.8715 48801.7867 -67061.1856 -143424.2107 61 62 63 64 65 113473.3182 -221187.1722 101960.8337 -136042.5409 -71852.6127 66 67 68 69 70 -776923.7604 -32851.0573 16382.4429 -131842.4345 -116188.8282 71 72 73 74 75 -222043.2780 -19362.1218 -29906.7408 49218.1317 67348.3632 76 77 78 79 80 -67411.4994 59451.8497 -54129.1308 -106918.2436 105260.6267 81 82 83 84 85 -11619.3933 -38571.8810 46526.3361 -48363.1395 113381.1226 86 87 88 89 90 -108117.6610 -129024.2839 99821.2971 -73269.7527 126925.8198 91 92 93 94 95 40367.0213 2882.4075 -12795.5059 67585.5605 -50190.5068 96 97 98 99 100 14263.1426 204777.8549 1082.3575 157549.7872 31984.5842 101 102 103 104 105 -43643.3516 -60082.1371 -123422.5483 95245.7345 19350.1987 106 107 108 109 110 -79434.6543 -285466.4373 236265.2051 -70196.6536 -70231.3401 111 112 113 114 115 76834.1596 -221896.4443 -95449.6266 -7581.8911 75656.5911 116 117 118 119 120 -35261.7376 68496.3892 152354.7718 -71939.2435 -21325.6342 121 122 123 124 125 -203185.2554 30067.5382 -279276.5402 181548.4371 29069.8249 126 127 128 129 130 -108933.2935 34415.4883 6684.3415 -28116.2424 34415.4883 131 132 133 134 135 14347.2110 15613.6376 34415.4883 -21905.9925 20312.5184 136 137 138 139 140 22880.8169 11285.0025 31222.9958 23150.5184 34415.4883 141 142 143 144 145 20214.4836 34415.4883 -17994.0601 24367.4815 11910.6657 146 147 148 149 150 1139.3479 34415.4883 15311.1015 38885.5034 35452.1628 151 152 153 154 155 -11542.0486 -99299.7376 20918.0259 -5467.2792 29265.6370 156 157 158 159 160 34415.4883 34415.4883 30229.3058 42171.8216 34415.4883 161 162 163 164 165 34415.4883 34415.4883 34415.4883 6584.0554 -40796.9083 166 167 168 169 170 31217.9958 34415.4883 34023.1804 -87866.5980 34415.4883 171 172 173 174 175 34415.4883 24810.4047 34415.4883 6221.6155 34415.4883 176 177 178 179 180 13351.2167 -272740.1617 35257.0689 34415.4883 34415.4883 181 182 183 184 185 40649.1911 35647.8860 -39381.2283 -11864.4442 23409.1469 186 187 188 189 190 22030.6439 13706.8251 18336.7011 5718.6132 34415.4883 191 192 193 194 195 16920.7901 34415.4883 16969.4330 34415.4883 26919.3675 196 197 198 199 200 34415.4883 -158975.0055 34415.4883 34415.4883 -38659.2549 201 202 203 204 205 34415.4883 34415.4883 34415.4883 -26837.7169 50961.7796 206 207 208 209 210 23256.4882 40964.5220 34415.4883 -20082.2236 17885.0109 211 212 213 214 215 34415.4883 34415.4883 14820.3119 -17066.9307 13196.7206 216 217 218 219 220 29016.5034 18708.5506 17755.9585 29016.5034 -38566.9285 221 222 223 224 225 7229.3939 34415.4883 31650.5575 -22793.3612 32064.3713 226 227 228 229 230 -24644.3051 9291.0560 34415.4883 34003.5446 34415.4883 231 232 233 234 235 3980.5053 16116.0560 34415.4883 -11450.5996 34415.4883 236 237 238 239 240 34564.4225 34415.4883 18023.1699 34415.4883 34415.4883 241 242 243 244 245 34415.4883 24932.1768 34415.4883 14324.5227 59095.9927 246 247 248 249 250 34415.4883 34415.4883 -10983.1623 34415.4883 34415.4883 251 252 253 254 255 -20338.8207 34415.4883 29977.0971 -71041.4238 -179472.7167 256 257 258 259 260 34415.4883 204.6311 -66069.3918 34415.4883 34415.4883 261 262 263 264 265 34415.4883 34415.4883 45144.5464 21166.9312 34415.4883 266 267 268 269 270 34415.4883 34415.4883 102136.7639 34415.4883 -10376.3459 271 272 273 274 275 34415.4883 34415.4883 34415.4883 -21956.9647 34415.4883 276 277 278 279 280 9179.9993 -18784.9860 34415.4883 -61114.5573 19816.8185 281 282 283 284 285 34415.4883 34415.4883 -41767.5768 34415.4883 41684.7152 286 287 288 289 290 34415.4883 24748.0259 31329.2606 34415.4883 -32424.7001 291 292 293 294 295 34415.4883 14906.5245 59054.4943 34415.4883 -7392.6895 296 297 298 299 300 33610.4098 34415.4883 34415.4883 34415.4883 10120.0560 301 302 303 304 305 21534.1915 14040.1758 31715.9958 34415.4883 27185.5923 306 307 308 309 310 2057.4240 34415.4883 -11513.0490 -24889.1373 34415.4883 311 312 313 314 315 34415.4883 1629.6336 34415.4883 7999.2017 13160.2732 316 317 318 319 320 56646.4587 32664.8879 -37271.6918 20908.3556 34415.4883 321 322 323 324 325 34415.4883 34415.4883 34415.4883 34415.4883 34415.4883 326 327 328 329 330 28841.5811 34415.4883 34415.4883 34415.4883 34415.4883 331 332 333 334 335 34415.4883 34415.4883 34415.4883 21386.4987 31736.5444 336 337 338 339 340 23608.5184 -38566.9446 34415.4883 30531.5027 8242.3020 341 342 343 344 345 -14503.0610 34415.4883 34415.4883 34415.4883 -7801.9139 346 347 348 349 350 34415.4883 11840.1331 34415.4883 34415.4883 11409.5635 351 352 353 354 355 25952.7483 14838.1667 16032.3686 32957.1139 34415.4883 356 357 358 359 360 34415.4883 25566.6364 34415.4883 34415.4883 34415.4883 361 362 363 364 365 23638.9203 34415.4883 34415.4883 -47363.1526 9446.3909 366 367 368 369 370 34415.4883 34415.4883 34415.4883 34415.4883 15692.0410 371 372 373 374 375 34415.4883 21988.9815 62403.2924 25460.1988 34415.4883 376 377 378 379 380 34415.4883 -22042.5868 22022.0219 11904.5635 34415.4883 381 382 383 384 385 34415.4883 3918.4103 -65207.6130 34415.4883 6855.1773 386 387 388 389 390 -114489.6195 36174.3426 8432.9664 24589.4649 1991.2370 391 392 393 394 395 17387.4017 -7312.9789 -38921.3620 82529.6201 14101.0516 396 397 398 399 400 -22547.5156 -115634.5973 -819140.4195 1469.0903 -183653.8205 401 402 403 404 405 31975.9930 36485.6108 -2751.3511 -48468.7720 -17250.2209 406 407 408 409 410 31178.7587 17219.3911 -11410.3112 -10879.3322 -41168.2135 411 412 413 414 415 -124676.3668 -47696.4693 -56343.7006 93154.6182 -117449.2482 416 417 418 419 420 -59536.7217 35577.8986 -42581.9122 -61195.8229 -27801.2093 421 422 423 424 425 -13972.4958 -245082.7760 -528552.8930 -383023.0073 -68961.9405 426 427 428 429 430 -291280.2801 -306083.6139 -1177595.1299 -97935.0948 -153965.3588 431 -464829.1094 > postscript(file="/var/www/rcomp/tmp/68u241292956049.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 340675.6874 NA 1 2594167.7524 340675.6874 2 1722550.3597 2594167.7524 3 -236477.4865 1722550.3597 4 -156944.7777 -236477.4865 5 575467.7008 -156944.7777 6 -18734.1681 575467.7008 7 -178890.7705 -18734.1681 8 88611.2273 -178890.7705 9 841746.1429 88611.2273 10 260708.2843 841746.1429 11 398795.4080 260708.2843 12 699614.9668 398795.4080 13 -454857.4185 699614.9668 14 -341002.5886 -454857.4185 15 610942.5762 -341002.5886 16 -212053.1222 610942.5762 17 -475997.8305 -212053.1222 18 201509.3228 -475997.8305 19 -997071.1842 201509.3228 20 370357.3971 -997071.1842 21 245843.6206 370357.3971 22 4400.2931 245843.6206 23 -1097247.5545 4400.2931 24 -259403.5402 -1097247.5545 25 153925.6423 -259403.5402 26 -337062.9511 153925.6423 27 -256461.0148 -337062.9511 28 391564.1640 -256461.0148 29 -190039.5441 391564.1640 30 364818.4903 -190039.5441 31 90029.4346 364818.4903 32 364381.7852 90029.4346 33 225777.9029 364381.7852 34 681742.2419 225777.9029 35 -665091.9651 681742.2419 36 147314.3445 -665091.9651 37 168314.3345 147314.3445 38 -48020.6220 168314.3345 39 -468840.2425 -48020.6220 40 -823521.3136 -468840.2425 41 -82380.0097 -823521.3136 42 -49175.0083 -82380.0097 43 -187787.8009 -49175.0083 44 6554.6514 -187787.8009 45 72345.5298 6554.6514 46 6328.8551 72345.5298 47 -460706.9652 6328.8551 48 34635.5935 -460706.9652 49 -188214.8335 34635.5935 50 -37223.1732 -188214.8335 51 -302577.5598 -37223.1732 52 -197893.8234 -302577.5598 53 67199.5523 -197893.8234 54 -514295.7847 67199.5523 55 -83186.7776 -514295.7847 56 338758.8715 -83186.7776 57 48801.7867 338758.8715 58 -67061.1856 48801.7867 59 -143424.2107 -67061.1856 60 113473.3182 -143424.2107 61 -221187.1722 113473.3182 62 101960.8337 -221187.1722 63 -136042.5409 101960.8337 64 -71852.6127 -136042.5409 65 -776923.7604 -71852.6127 66 -32851.0573 -776923.7604 67 16382.4429 -32851.0573 68 -131842.4345 16382.4429 69 -116188.8282 -131842.4345 70 -222043.2780 -116188.8282 71 -19362.1218 -222043.2780 72 -29906.7408 -19362.1218 73 49218.1317 -29906.7408 74 67348.3632 49218.1317 75 -67411.4994 67348.3632 76 59451.8497 -67411.4994 77 -54129.1308 59451.8497 78 -106918.2436 -54129.1308 79 105260.6267 -106918.2436 80 -11619.3933 105260.6267 81 -38571.8810 -11619.3933 82 46526.3361 -38571.8810 83 -48363.1395 46526.3361 84 113381.1226 -48363.1395 85 -108117.6610 113381.1226 86 -129024.2839 -108117.6610 87 99821.2971 -129024.2839 88 -73269.7527 99821.2971 89 126925.8198 -73269.7527 90 40367.0213 126925.8198 91 2882.4075 40367.0213 92 -12795.5059 2882.4075 93 67585.5605 -12795.5059 94 -50190.5068 67585.5605 95 14263.1426 -50190.5068 96 204777.8549 14263.1426 97 1082.3575 204777.8549 98 157549.7872 1082.3575 99 31984.5842 157549.7872 100 -43643.3516 31984.5842 101 -60082.1371 -43643.3516 102 -123422.5483 -60082.1371 103 95245.7345 -123422.5483 104 19350.1987 95245.7345 105 -79434.6543 19350.1987 106 -285466.4373 -79434.6543 107 236265.2051 -285466.4373 108 -70196.6536 236265.2051 109 -70231.3401 -70196.6536 110 76834.1596 -70231.3401 111 -221896.4443 76834.1596 112 -95449.6266 -221896.4443 113 -7581.8911 -95449.6266 114 75656.5911 -7581.8911 115 -35261.7376 75656.5911 116 68496.3892 -35261.7376 117 152354.7718 68496.3892 118 -71939.2435 152354.7718 119 -21325.6342 -71939.2435 120 -203185.2554 -21325.6342 121 30067.5382 -203185.2554 122 -279276.5402 30067.5382 123 181548.4371 -279276.5402 124 29069.8249 181548.4371 125 -108933.2935 29069.8249 126 34415.4883 -108933.2935 127 6684.3415 34415.4883 128 -28116.2424 6684.3415 129 34415.4883 -28116.2424 130 14347.2110 34415.4883 131 15613.6376 14347.2110 132 34415.4883 15613.6376 133 -21905.9925 34415.4883 134 20312.5184 -21905.9925 135 22880.8169 20312.5184 136 11285.0025 22880.8169 137 31222.9958 11285.0025 138 23150.5184 31222.9958 139 34415.4883 23150.5184 140 20214.4836 34415.4883 141 34415.4883 20214.4836 142 -17994.0601 34415.4883 143 24367.4815 -17994.0601 144 11910.6657 24367.4815 145 1139.3479 11910.6657 146 34415.4883 1139.3479 147 15311.1015 34415.4883 148 38885.5034 15311.1015 149 35452.1628 38885.5034 150 -11542.0486 35452.1628 151 -99299.7376 -11542.0486 152 20918.0259 -99299.7376 153 -5467.2792 20918.0259 154 29265.6370 -5467.2792 155 34415.4883 29265.6370 156 34415.4883 34415.4883 157 30229.3058 34415.4883 158 42171.8216 30229.3058 159 34415.4883 42171.8216 160 34415.4883 34415.4883 161 34415.4883 34415.4883 162 34415.4883 34415.4883 163 6584.0554 34415.4883 164 -40796.9083 6584.0554 165 31217.9958 -40796.9083 166 34415.4883 31217.9958 167 34023.1804 34415.4883 168 -87866.5980 34023.1804 169 34415.4883 -87866.5980 170 34415.4883 34415.4883 171 24810.4047 34415.4883 172 34415.4883 24810.4047 173 6221.6155 34415.4883 174 34415.4883 6221.6155 175 13351.2167 34415.4883 176 -272740.1617 13351.2167 177 35257.0689 -272740.1617 178 34415.4883 35257.0689 179 34415.4883 34415.4883 180 40649.1911 34415.4883 181 35647.8860 40649.1911 182 -39381.2283 35647.8860 183 -11864.4442 -39381.2283 184 23409.1469 -11864.4442 185 22030.6439 23409.1469 186 13706.8251 22030.6439 187 18336.7011 13706.8251 188 5718.6132 18336.7011 189 34415.4883 5718.6132 190 16920.7901 34415.4883 191 34415.4883 16920.7901 192 16969.4330 34415.4883 193 34415.4883 16969.4330 194 26919.3675 34415.4883 195 34415.4883 26919.3675 196 -158975.0055 34415.4883 197 34415.4883 -158975.0055 198 34415.4883 34415.4883 199 -38659.2549 34415.4883 200 34415.4883 -38659.2549 201 34415.4883 34415.4883 202 34415.4883 34415.4883 203 -26837.7169 34415.4883 204 50961.7796 -26837.7169 205 23256.4882 50961.7796 206 40964.5220 23256.4882 207 34415.4883 40964.5220 208 -20082.2236 34415.4883 209 17885.0109 -20082.2236 210 34415.4883 17885.0109 211 34415.4883 34415.4883 212 14820.3119 34415.4883 213 -17066.9307 14820.3119 214 13196.7206 -17066.9307 215 29016.5034 13196.7206 216 18708.5506 29016.5034 217 17755.9585 18708.5506 218 29016.5034 17755.9585 219 -38566.9285 29016.5034 220 7229.3939 -38566.9285 221 34415.4883 7229.3939 222 31650.5575 34415.4883 223 -22793.3612 31650.5575 224 32064.3713 -22793.3612 225 -24644.3051 32064.3713 226 9291.0560 -24644.3051 227 34415.4883 9291.0560 228 34003.5446 34415.4883 229 34415.4883 34003.5446 230 3980.5053 34415.4883 231 16116.0560 3980.5053 232 34415.4883 16116.0560 233 -11450.5996 34415.4883 234 34415.4883 -11450.5996 235 34564.4225 34415.4883 236 34415.4883 34564.4225 237 18023.1699 34415.4883 238 34415.4883 18023.1699 239 34415.4883 34415.4883 240 34415.4883 34415.4883 241 24932.1768 34415.4883 242 34415.4883 24932.1768 243 14324.5227 34415.4883 244 59095.9927 14324.5227 245 34415.4883 59095.9927 246 34415.4883 34415.4883 247 -10983.1623 34415.4883 248 34415.4883 -10983.1623 249 34415.4883 34415.4883 250 -20338.8207 34415.4883 251 34415.4883 -20338.8207 252 29977.0971 34415.4883 253 -71041.4238 29977.0971 254 -179472.7167 -71041.4238 255 34415.4883 -179472.7167 256 204.6311 34415.4883 257 -66069.3918 204.6311 258 34415.4883 -66069.3918 259 34415.4883 34415.4883 260 34415.4883 34415.4883 261 34415.4883 34415.4883 262 45144.5464 34415.4883 263 21166.9312 45144.5464 264 34415.4883 21166.9312 265 34415.4883 34415.4883 266 34415.4883 34415.4883 267 102136.7639 34415.4883 268 34415.4883 102136.7639 269 -10376.3459 34415.4883 270 34415.4883 -10376.3459 271 34415.4883 34415.4883 272 34415.4883 34415.4883 273 -21956.9647 34415.4883 274 34415.4883 -21956.9647 275 9179.9993 34415.4883 276 -18784.9860 9179.9993 277 34415.4883 -18784.9860 278 -61114.5573 34415.4883 279 19816.8185 -61114.5573 280 34415.4883 19816.8185 281 34415.4883 34415.4883 282 -41767.5768 34415.4883 283 34415.4883 -41767.5768 284 41684.7152 34415.4883 285 34415.4883 41684.7152 286 24748.0259 34415.4883 287 31329.2606 24748.0259 288 34415.4883 31329.2606 289 -32424.7001 34415.4883 290 34415.4883 -32424.7001 291 14906.5245 34415.4883 292 59054.4943 14906.5245 293 34415.4883 59054.4943 294 -7392.6895 34415.4883 295 33610.4098 -7392.6895 296 34415.4883 33610.4098 297 34415.4883 34415.4883 298 34415.4883 34415.4883 299 10120.0560 34415.4883 300 21534.1915 10120.0560 301 14040.1758 21534.1915 302 31715.9958 14040.1758 303 34415.4883 31715.9958 304 27185.5923 34415.4883 305 2057.4240 27185.5923 306 34415.4883 2057.4240 307 -11513.0490 34415.4883 308 -24889.1373 -11513.0490 309 34415.4883 -24889.1373 310 34415.4883 34415.4883 311 1629.6336 34415.4883 312 34415.4883 1629.6336 313 7999.2017 34415.4883 314 13160.2732 7999.2017 315 56646.4587 13160.2732 316 32664.8879 56646.4587 317 -37271.6918 32664.8879 318 20908.3556 -37271.6918 319 34415.4883 20908.3556 320 34415.4883 34415.4883 321 34415.4883 34415.4883 322 34415.4883 34415.4883 323 34415.4883 34415.4883 324 34415.4883 34415.4883 325 28841.5811 34415.4883 326 34415.4883 28841.5811 327 34415.4883 34415.4883 328 34415.4883 34415.4883 329 34415.4883 34415.4883 330 34415.4883 34415.4883 331 34415.4883 34415.4883 332 34415.4883 34415.4883 333 21386.4987 34415.4883 334 31736.5444 21386.4987 335 23608.5184 31736.5444 336 -38566.9446 23608.5184 337 34415.4883 -38566.9446 338 30531.5027 34415.4883 339 8242.3020 30531.5027 340 -14503.0610 8242.3020 341 34415.4883 -14503.0610 342 34415.4883 34415.4883 343 34415.4883 34415.4883 344 -7801.9139 34415.4883 345 34415.4883 -7801.9139 346 11840.1331 34415.4883 347 34415.4883 11840.1331 348 34415.4883 34415.4883 349 11409.5635 34415.4883 350 25952.7483 11409.5635 351 14838.1667 25952.7483 352 16032.3686 14838.1667 353 32957.1139 16032.3686 354 34415.4883 32957.1139 355 34415.4883 34415.4883 356 25566.6364 34415.4883 357 34415.4883 25566.6364 358 34415.4883 34415.4883 359 34415.4883 34415.4883 360 23638.9203 34415.4883 361 34415.4883 23638.9203 362 34415.4883 34415.4883 363 -47363.1526 34415.4883 364 9446.3909 -47363.1526 365 34415.4883 9446.3909 366 34415.4883 34415.4883 367 34415.4883 34415.4883 368 34415.4883 34415.4883 369 15692.0410 34415.4883 370 34415.4883 15692.0410 371 21988.9815 34415.4883 372 62403.2924 21988.9815 373 25460.1988 62403.2924 374 34415.4883 25460.1988 375 34415.4883 34415.4883 376 -22042.5868 34415.4883 377 22022.0219 -22042.5868 378 11904.5635 22022.0219 379 34415.4883 11904.5635 380 34415.4883 34415.4883 381 3918.4103 34415.4883 382 -65207.6130 3918.4103 383 34415.4883 -65207.6130 384 6855.1773 34415.4883 385 -114489.6195 6855.1773 386 36174.3426 -114489.6195 387 8432.9664 36174.3426 388 24589.4649 8432.9664 389 1991.2370 24589.4649 390 17387.4017 1991.2370 391 -7312.9789 17387.4017 392 -38921.3620 -7312.9789 393 82529.6201 -38921.3620 394 14101.0516 82529.6201 395 -22547.5156 14101.0516 396 -115634.5973 -22547.5156 397 -819140.4195 -115634.5973 398 1469.0903 -819140.4195 399 -183653.8205 1469.0903 400 31975.9930 -183653.8205 401 36485.6108 31975.9930 402 -2751.3511 36485.6108 403 -48468.7720 -2751.3511 404 -17250.2209 -48468.7720 405 31178.7587 -17250.2209 406 17219.3911 31178.7587 407 -11410.3112 17219.3911 408 -10879.3322 -11410.3112 409 -41168.2135 -10879.3322 410 -124676.3668 -41168.2135 411 -47696.4693 -124676.3668 412 -56343.7006 -47696.4693 413 93154.6182 -56343.7006 414 -117449.2482 93154.6182 415 -59536.7217 -117449.2482 416 35577.8986 -59536.7217 417 -42581.9122 35577.8986 418 -61195.8229 -42581.9122 419 -27801.2093 -61195.8229 420 -13972.4958 -27801.2093 421 -245082.7760 -13972.4958 422 -528552.8930 -245082.7760 423 -383023.0073 -528552.8930 424 -68961.9405 -383023.0073 425 -291280.2801 -68961.9405 426 -306083.6139 -291280.2801 427 -1177595.1299 -306083.6139 428 -97935.0948 -1177595.1299 429 -153965.3588 -97935.0948 430 -464829.1094 -153965.3588 431 NA -464829.1094 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2594167.7524 340675.6874 [2,] 1722550.3597 2594167.7524 [3,] -236477.4865 1722550.3597 [4,] -156944.7777 -236477.4865 [5,] 575467.7008 -156944.7777 [6,] -18734.1681 575467.7008 [7,] -178890.7705 -18734.1681 [8,] 88611.2273 -178890.7705 [9,] 841746.1429 88611.2273 [10,] 260708.2843 841746.1429 [11,] 398795.4080 260708.2843 [12,] 699614.9668 398795.4080 [13,] -454857.4185 699614.9668 [14,] -341002.5886 -454857.4185 [15,] 610942.5762 -341002.5886 [16,] -212053.1222 610942.5762 [17,] -475997.8305 -212053.1222 [18,] 201509.3228 -475997.8305 [19,] -997071.1842 201509.3228 [20,] 370357.3971 -997071.1842 [21,] 245843.6206 370357.3971 [22,] 4400.2931 245843.6206 [23,] -1097247.5545 4400.2931 [24,] -259403.5402 -1097247.5545 [25,] 153925.6423 -259403.5402 [26,] -337062.9511 153925.6423 [27,] -256461.0148 -337062.9511 [28,] 391564.1640 -256461.0148 [29,] -190039.5441 391564.1640 [30,] 364818.4903 -190039.5441 [31,] 90029.4346 364818.4903 [32,] 364381.7852 90029.4346 [33,] 225777.9029 364381.7852 [34,] 681742.2419 225777.9029 [35,] -665091.9651 681742.2419 [36,] 147314.3445 -665091.9651 [37,] 168314.3345 147314.3445 [38,] -48020.6220 168314.3345 [39,] -468840.2425 -48020.6220 [40,] -823521.3136 -468840.2425 [41,] -82380.0097 -823521.3136 [42,] -49175.0083 -82380.0097 [43,] -187787.8009 -49175.0083 [44,] 6554.6514 -187787.8009 [45,] 72345.5298 6554.6514 [46,] 6328.8551 72345.5298 [47,] -460706.9652 6328.8551 [48,] 34635.5935 -460706.9652 [49,] -188214.8335 34635.5935 [50,] -37223.1732 -188214.8335 [51,] -302577.5598 -37223.1732 [52,] -197893.8234 -302577.5598 [53,] 67199.5523 -197893.8234 [54,] -514295.7847 67199.5523 [55,] -83186.7776 -514295.7847 [56,] 338758.8715 -83186.7776 [57,] 48801.7867 338758.8715 [58,] -67061.1856 48801.7867 [59,] -143424.2107 -67061.1856 [60,] 113473.3182 -143424.2107 [61,] -221187.1722 113473.3182 [62,] 101960.8337 -221187.1722 [63,] -136042.5409 101960.8337 [64,] -71852.6127 -136042.5409 [65,] -776923.7604 -71852.6127 [66,] -32851.0573 -776923.7604 [67,] 16382.4429 -32851.0573 [68,] -131842.4345 16382.4429 [69,] -116188.8282 -131842.4345 [70,] -222043.2780 -116188.8282 [71,] -19362.1218 -222043.2780 [72,] -29906.7408 -19362.1218 [73,] 49218.1317 -29906.7408 [74,] 67348.3632 49218.1317 [75,] -67411.4994 67348.3632 [76,] 59451.8497 -67411.4994 [77,] -54129.1308 59451.8497 [78,] -106918.2436 -54129.1308 [79,] 105260.6267 -106918.2436 [80,] -11619.3933 105260.6267 [81,] -38571.8810 -11619.3933 [82,] 46526.3361 -38571.8810 [83,] -48363.1395 46526.3361 [84,] 113381.1226 -48363.1395 [85,] -108117.6610 113381.1226 [86,] -129024.2839 -108117.6610 [87,] 99821.2971 -129024.2839 [88,] -73269.7527 99821.2971 [89,] 126925.8198 -73269.7527 [90,] 40367.0213 126925.8198 [91,] 2882.4075 40367.0213 [92,] -12795.5059 2882.4075 [93,] 67585.5605 -12795.5059 [94,] -50190.5068 67585.5605 [95,] 14263.1426 -50190.5068 [96,] 204777.8549 14263.1426 [97,] 1082.3575 204777.8549 [98,] 157549.7872 1082.3575 [99,] 31984.5842 157549.7872 [100,] -43643.3516 31984.5842 [101,] -60082.1371 -43643.3516 [102,] -123422.5483 -60082.1371 [103,] 95245.7345 -123422.5483 [104,] 19350.1987 95245.7345 [105,] -79434.6543 19350.1987 [106,] -285466.4373 -79434.6543 [107,] 236265.2051 -285466.4373 [108,] -70196.6536 236265.2051 [109,] -70231.3401 -70196.6536 [110,] 76834.1596 -70231.3401 [111,] -221896.4443 76834.1596 [112,] -95449.6266 -221896.4443 [113,] -7581.8911 -95449.6266 [114,] 75656.5911 -7581.8911 [115,] -35261.7376 75656.5911 [116,] 68496.3892 -35261.7376 [117,] 152354.7718 68496.3892 [118,] -71939.2435 152354.7718 [119,] -21325.6342 -71939.2435 [120,] -203185.2554 -21325.6342 [121,] 30067.5382 -203185.2554 [122,] -279276.5402 30067.5382 [123,] 181548.4371 -279276.5402 [124,] 29069.8249 181548.4371 [125,] -108933.2935 29069.8249 [126,] 34415.4883 -108933.2935 [127,] 6684.3415 34415.4883 [128,] -28116.2424 6684.3415 [129,] 34415.4883 -28116.2424 [130,] 14347.2110 34415.4883 [131,] 15613.6376 14347.2110 [132,] 34415.4883 15613.6376 [133,] -21905.9925 34415.4883 [134,] 20312.5184 -21905.9925 [135,] 22880.8169 20312.5184 [136,] 11285.0025 22880.8169 [137,] 31222.9958 11285.0025 [138,] 23150.5184 31222.9958 [139,] 34415.4883 23150.5184 [140,] 20214.4836 34415.4883 [141,] 34415.4883 20214.4836 [142,] -17994.0601 34415.4883 [143,] 24367.4815 -17994.0601 [144,] 11910.6657 24367.4815 [145,] 1139.3479 11910.6657 [146,] 34415.4883 1139.3479 [147,] 15311.1015 34415.4883 [148,] 38885.5034 15311.1015 [149,] 35452.1628 38885.5034 [150,] -11542.0486 35452.1628 [151,] -99299.7376 -11542.0486 [152,] 20918.0259 -99299.7376 [153,] -5467.2792 20918.0259 [154,] 29265.6370 -5467.2792 [155,] 34415.4883 29265.6370 [156,] 34415.4883 34415.4883 [157,] 30229.3058 34415.4883 [158,] 42171.8216 30229.3058 [159,] 34415.4883 42171.8216 [160,] 34415.4883 34415.4883 [161,] 34415.4883 34415.4883 [162,] 34415.4883 34415.4883 [163,] 6584.0554 34415.4883 [164,] -40796.9083 6584.0554 [165,] 31217.9958 -40796.9083 [166,] 34415.4883 31217.9958 [167,] 34023.1804 34415.4883 [168,] -87866.5980 34023.1804 [169,] 34415.4883 -87866.5980 [170,] 34415.4883 34415.4883 [171,] 24810.4047 34415.4883 [172,] 34415.4883 24810.4047 [173,] 6221.6155 34415.4883 [174,] 34415.4883 6221.6155 [175,] 13351.2167 34415.4883 [176,] -272740.1617 13351.2167 [177,] 35257.0689 -272740.1617 [178,] 34415.4883 35257.0689 [179,] 34415.4883 34415.4883 [180,] 40649.1911 34415.4883 [181,] 35647.8860 40649.1911 [182,] -39381.2283 35647.8860 [183,] -11864.4442 -39381.2283 [184,] 23409.1469 -11864.4442 [185,] 22030.6439 23409.1469 [186,] 13706.8251 22030.6439 [187,] 18336.7011 13706.8251 [188,] 5718.6132 18336.7011 [189,] 34415.4883 5718.6132 [190,] 16920.7901 34415.4883 [191,] 34415.4883 16920.7901 [192,] 16969.4330 34415.4883 [193,] 34415.4883 16969.4330 [194,] 26919.3675 34415.4883 [195,] 34415.4883 26919.3675 [196,] -158975.0055 34415.4883 [197,] 34415.4883 -158975.0055 [198,] 34415.4883 34415.4883 [199,] -38659.2549 34415.4883 [200,] 34415.4883 -38659.2549 [201,] 34415.4883 34415.4883 [202,] 34415.4883 34415.4883 [203,] -26837.7169 34415.4883 [204,] 50961.7796 -26837.7169 [205,] 23256.4882 50961.7796 [206,] 40964.5220 23256.4882 [207,] 34415.4883 40964.5220 [208,] -20082.2236 34415.4883 [209,] 17885.0109 -20082.2236 [210,] 34415.4883 17885.0109 [211,] 34415.4883 34415.4883 [212,] 14820.3119 34415.4883 [213,] -17066.9307 14820.3119 [214,] 13196.7206 -17066.9307 [215,] 29016.5034 13196.7206 [216,] 18708.5506 29016.5034 [217,] 17755.9585 18708.5506 [218,] 29016.5034 17755.9585 [219,] -38566.9285 29016.5034 [220,] 7229.3939 -38566.9285 [221,] 34415.4883 7229.3939 [222,] 31650.5575 34415.4883 [223,] -22793.3612 31650.5575 [224,] 32064.3713 -22793.3612 [225,] -24644.3051 32064.3713 [226,] 9291.0560 -24644.3051 [227,] 34415.4883 9291.0560 [228,] 34003.5446 34415.4883 [229,] 34415.4883 34003.5446 [230,] 3980.5053 34415.4883 [231,] 16116.0560 3980.5053 [232,] 34415.4883 16116.0560 [233,] -11450.5996 34415.4883 [234,] 34415.4883 -11450.5996 [235,] 34564.4225 34415.4883 [236,] 34415.4883 34564.4225 [237,] 18023.1699 34415.4883 [238,] 34415.4883 18023.1699 [239,] 34415.4883 34415.4883 [240,] 34415.4883 34415.4883 [241,] 24932.1768 34415.4883 [242,] 34415.4883 24932.1768 [243,] 14324.5227 34415.4883 [244,] 59095.9927 14324.5227 [245,] 34415.4883 59095.9927 [246,] 34415.4883 34415.4883 [247,] -10983.1623 34415.4883 [248,] 34415.4883 -10983.1623 [249,] 34415.4883 34415.4883 [250,] -20338.8207 34415.4883 [251,] 34415.4883 -20338.8207 [252,] 29977.0971 34415.4883 [253,] -71041.4238 29977.0971 [254,] -179472.7167 -71041.4238 [255,] 34415.4883 -179472.7167 [256,] 204.6311 34415.4883 [257,] -66069.3918 204.6311 [258,] 34415.4883 -66069.3918 [259,] 34415.4883 34415.4883 [260,] 34415.4883 34415.4883 [261,] 34415.4883 34415.4883 [262,] 45144.5464 34415.4883 [263,] 21166.9312 45144.5464 [264,] 34415.4883 21166.9312 [265,] 34415.4883 34415.4883 [266,] 34415.4883 34415.4883 [267,] 102136.7639 34415.4883 [268,] 34415.4883 102136.7639 [269,] -10376.3459 34415.4883 [270,] 34415.4883 -10376.3459 [271,] 34415.4883 34415.4883 [272,] 34415.4883 34415.4883 [273,] -21956.9647 34415.4883 [274,] 34415.4883 -21956.9647 [275,] 9179.9993 34415.4883 [276,] -18784.9860 9179.9993 [277,] 34415.4883 -18784.9860 [278,] -61114.5573 34415.4883 [279,] 19816.8185 -61114.5573 [280,] 34415.4883 19816.8185 [281,] 34415.4883 34415.4883 [282,] -41767.5768 34415.4883 [283,] 34415.4883 -41767.5768 [284,] 41684.7152 34415.4883 [285,] 34415.4883 41684.7152 [286,] 24748.0259 34415.4883 [287,] 31329.2606 24748.0259 [288,] 34415.4883 31329.2606 [289,] -32424.7001 34415.4883 [290,] 34415.4883 -32424.7001 [291,] 14906.5245 34415.4883 [292,] 59054.4943 14906.5245 [293,] 34415.4883 59054.4943 [294,] -7392.6895 34415.4883 [295,] 33610.4098 -7392.6895 [296,] 34415.4883 33610.4098 [297,] 34415.4883 34415.4883 [298,] 34415.4883 34415.4883 [299,] 10120.0560 34415.4883 [300,] 21534.1915 10120.0560 [301,] 14040.1758 21534.1915 [302,] 31715.9958 14040.1758 [303,] 34415.4883 31715.9958 [304,] 27185.5923 34415.4883 [305,] 2057.4240 27185.5923 [306,] 34415.4883 2057.4240 [307,] -11513.0490 34415.4883 [308,] -24889.1373 -11513.0490 [309,] 34415.4883 -24889.1373 [310,] 34415.4883 34415.4883 [311,] 1629.6336 34415.4883 [312,] 34415.4883 1629.6336 [313,] 7999.2017 34415.4883 [314,] 13160.2732 7999.2017 [315,] 56646.4587 13160.2732 [316,] 32664.8879 56646.4587 [317,] -37271.6918 32664.8879 [318,] 20908.3556 -37271.6918 [319,] 34415.4883 20908.3556 [320,] 34415.4883 34415.4883 [321,] 34415.4883 34415.4883 [322,] 34415.4883 34415.4883 [323,] 34415.4883 34415.4883 [324,] 34415.4883 34415.4883 [325,] 28841.5811 34415.4883 [326,] 34415.4883 28841.5811 [327,] 34415.4883 34415.4883 [328,] 34415.4883 34415.4883 [329,] 34415.4883 34415.4883 [330,] 34415.4883 34415.4883 [331,] 34415.4883 34415.4883 [332,] 34415.4883 34415.4883 [333,] 21386.4987 34415.4883 [334,] 31736.5444 21386.4987 [335,] 23608.5184 31736.5444 [336,] -38566.9446 23608.5184 [337,] 34415.4883 -38566.9446 [338,] 30531.5027 34415.4883 [339,] 8242.3020 30531.5027 [340,] -14503.0610 8242.3020 [341,] 34415.4883 -14503.0610 [342,] 34415.4883 34415.4883 [343,] 34415.4883 34415.4883 [344,] -7801.9139 34415.4883 [345,] 34415.4883 -7801.9139 [346,] 11840.1331 34415.4883 [347,] 34415.4883 11840.1331 [348,] 34415.4883 34415.4883 [349,] 11409.5635 34415.4883 [350,] 25952.7483 11409.5635 [351,] 14838.1667 25952.7483 [352,] 16032.3686 14838.1667 [353,] 32957.1139 16032.3686 [354,] 34415.4883 32957.1139 [355,] 34415.4883 34415.4883 [356,] 25566.6364 34415.4883 [357,] 34415.4883 25566.6364 [358,] 34415.4883 34415.4883 [359,] 34415.4883 34415.4883 [360,] 23638.9203 34415.4883 [361,] 34415.4883 23638.9203 [362,] 34415.4883 34415.4883 [363,] -47363.1526 34415.4883 [364,] 9446.3909 -47363.1526 [365,] 34415.4883 9446.3909 [366,] 34415.4883 34415.4883 [367,] 34415.4883 34415.4883 [368,] 34415.4883 34415.4883 [369,] 15692.0410 34415.4883 [370,] 34415.4883 15692.0410 [371,] 21988.9815 34415.4883 [372,] 62403.2924 21988.9815 [373,] 25460.1988 62403.2924 [374,] 34415.4883 25460.1988 [375,] 34415.4883 34415.4883 [376,] -22042.5868 34415.4883 [377,] 22022.0219 -22042.5868 [378,] 11904.5635 22022.0219 [379,] 34415.4883 11904.5635 [380,] 34415.4883 34415.4883 [381,] 3918.4103 34415.4883 [382,] -65207.6130 3918.4103 [383,] 34415.4883 -65207.6130 [384,] 6855.1773 34415.4883 [385,] -114489.6195 6855.1773 [386,] 36174.3426 -114489.6195 [387,] 8432.9664 36174.3426 [388,] 24589.4649 8432.9664 [389,] 1991.2370 24589.4649 [390,] 17387.4017 1991.2370 [391,] -7312.9789 17387.4017 [392,] -38921.3620 -7312.9789 [393,] 82529.6201 -38921.3620 [394,] 14101.0516 82529.6201 [395,] -22547.5156 14101.0516 [396,] -115634.5973 -22547.5156 [397,] -819140.4195 -115634.5973 [398,] 1469.0903 -819140.4195 [399,] -183653.8205 1469.0903 [400,] 31975.9930 -183653.8205 [401,] 36485.6108 31975.9930 [402,] -2751.3511 36485.6108 [403,] -48468.7720 -2751.3511 [404,] -17250.2209 -48468.7720 [405,] 31178.7587 -17250.2209 [406,] 17219.3911 31178.7587 [407,] -11410.3112 17219.3911 [408,] -10879.3322 -11410.3112 [409,] -41168.2135 -10879.3322 [410,] -124676.3668 -41168.2135 [411,] -47696.4693 -124676.3668 [412,] -56343.7006 -47696.4693 [413,] 93154.6182 -56343.7006 [414,] -117449.2482 93154.6182 [415,] -59536.7217 -117449.2482 [416,] 35577.8986 -59536.7217 [417,] -42581.9122 35577.8986 [418,] -61195.8229 -42581.9122 [419,] -27801.2093 -61195.8229 [420,] -13972.4958 -27801.2093 [421,] -245082.7760 -13972.4958 [422,] -528552.8930 -245082.7760 [423,] -383023.0073 -528552.8930 [424,] -68961.9405 -383023.0073 [425,] -291280.2801 -68961.9405 [426,] -306083.6139 -291280.2801 [427,] -1177595.1299 -306083.6139 [428,] -97935.0948 -1177595.1299 [429,] -153965.3588 -97935.0948 [430,] -464829.1094 -153965.3588 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2594167.7524 340675.6874 2 1722550.3597 2594167.7524 3 -236477.4865 1722550.3597 4 -156944.7777 -236477.4865 5 575467.7008 -156944.7777 6 -18734.1681 575467.7008 7 -178890.7705 -18734.1681 8 88611.2273 -178890.7705 9 841746.1429 88611.2273 10 260708.2843 841746.1429 11 398795.4080 260708.2843 12 699614.9668 398795.4080 13 -454857.4185 699614.9668 14 -341002.5886 -454857.4185 15 610942.5762 -341002.5886 16 -212053.1222 610942.5762 17 -475997.8305 -212053.1222 18 201509.3228 -475997.8305 19 -997071.1842 201509.3228 20 370357.3971 -997071.1842 21 245843.6206 370357.3971 22 4400.2931 245843.6206 23 -1097247.5545 4400.2931 24 -259403.5402 -1097247.5545 25 153925.6423 -259403.5402 26 -337062.9511 153925.6423 27 -256461.0148 -337062.9511 28 391564.1640 -256461.0148 29 -190039.5441 391564.1640 30 364818.4903 -190039.5441 31 90029.4346 364818.4903 32 364381.7852 90029.4346 33 225777.9029 364381.7852 34 681742.2419 225777.9029 35 -665091.9651 681742.2419 36 147314.3445 -665091.9651 37 168314.3345 147314.3445 38 -48020.6220 168314.3345 39 -468840.2425 -48020.6220 40 -823521.3136 -468840.2425 41 -82380.0097 -823521.3136 42 -49175.0083 -82380.0097 43 -187787.8009 -49175.0083 44 6554.6514 -187787.8009 45 72345.5298 6554.6514 46 6328.8551 72345.5298 47 -460706.9652 6328.8551 48 34635.5935 -460706.9652 49 -188214.8335 34635.5935 50 -37223.1732 -188214.8335 51 -302577.5598 -37223.1732 52 -197893.8234 -302577.5598 53 67199.5523 -197893.8234 54 -514295.7847 67199.5523 55 -83186.7776 -514295.7847 56 338758.8715 -83186.7776 57 48801.7867 338758.8715 58 -67061.1856 48801.7867 59 -143424.2107 -67061.1856 60 113473.3182 -143424.2107 61 -221187.1722 113473.3182 62 101960.8337 -221187.1722 63 -136042.5409 101960.8337 64 -71852.6127 -136042.5409 65 -776923.7604 -71852.6127 66 -32851.0573 -776923.7604 67 16382.4429 -32851.0573 68 -131842.4345 16382.4429 69 -116188.8282 -131842.4345 70 -222043.2780 -116188.8282 71 -19362.1218 -222043.2780 72 -29906.7408 -19362.1218 73 49218.1317 -29906.7408 74 67348.3632 49218.1317 75 -67411.4994 67348.3632 76 59451.8497 -67411.4994 77 -54129.1308 59451.8497 78 -106918.2436 -54129.1308 79 105260.6267 -106918.2436 80 -11619.3933 105260.6267 81 -38571.8810 -11619.3933 82 46526.3361 -38571.8810 83 -48363.1395 46526.3361 84 113381.1226 -48363.1395 85 -108117.6610 113381.1226 86 -129024.2839 -108117.6610 87 99821.2971 -129024.2839 88 -73269.7527 99821.2971 89 126925.8198 -73269.7527 90 40367.0213 126925.8198 91 2882.4075 40367.0213 92 -12795.5059 2882.4075 93 67585.5605 -12795.5059 94 -50190.5068 67585.5605 95 14263.1426 -50190.5068 96 204777.8549 14263.1426 97 1082.3575 204777.8549 98 157549.7872 1082.3575 99 31984.5842 157549.7872 100 -43643.3516 31984.5842 101 -60082.1371 -43643.3516 102 -123422.5483 -60082.1371 103 95245.7345 -123422.5483 104 19350.1987 95245.7345 105 -79434.6543 19350.1987 106 -285466.4373 -79434.6543 107 236265.2051 -285466.4373 108 -70196.6536 236265.2051 109 -70231.3401 -70196.6536 110 76834.1596 -70231.3401 111 -221896.4443 76834.1596 112 -95449.6266 -221896.4443 113 -7581.8911 -95449.6266 114 75656.5911 -7581.8911 115 -35261.7376 75656.5911 116 68496.3892 -35261.7376 117 152354.7718 68496.3892 118 -71939.2435 152354.7718 119 -21325.6342 -71939.2435 120 -203185.2554 -21325.6342 121 30067.5382 -203185.2554 122 -279276.5402 30067.5382 123 181548.4371 -279276.5402 124 29069.8249 181548.4371 125 -108933.2935 29069.8249 126 34415.4883 -108933.2935 127 6684.3415 34415.4883 128 -28116.2424 6684.3415 129 34415.4883 -28116.2424 130 14347.2110 34415.4883 131 15613.6376 14347.2110 132 34415.4883 15613.6376 133 -21905.9925 34415.4883 134 20312.5184 -21905.9925 135 22880.8169 20312.5184 136 11285.0025 22880.8169 137 31222.9958 11285.0025 138 23150.5184 31222.9958 139 34415.4883 23150.5184 140 20214.4836 34415.4883 141 34415.4883 20214.4836 142 -17994.0601 34415.4883 143 24367.4815 -17994.0601 144 11910.6657 24367.4815 145 1139.3479 11910.6657 146 34415.4883 1139.3479 147 15311.1015 34415.4883 148 38885.5034 15311.1015 149 35452.1628 38885.5034 150 -11542.0486 35452.1628 151 -99299.7376 -11542.0486 152 20918.0259 -99299.7376 153 -5467.2792 20918.0259 154 29265.6370 -5467.2792 155 34415.4883 29265.6370 156 34415.4883 34415.4883 157 30229.3058 34415.4883 158 42171.8216 30229.3058 159 34415.4883 42171.8216 160 34415.4883 34415.4883 161 34415.4883 34415.4883 162 34415.4883 34415.4883 163 6584.0554 34415.4883 164 -40796.9083 6584.0554 165 31217.9958 -40796.9083 166 34415.4883 31217.9958 167 34023.1804 34415.4883 168 -87866.5980 34023.1804 169 34415.4883 -87866.5980 170 34415.4883 34415.4883 171 24810.4047 34415.4883 172 34415.4883 24810.4047 173 6221.6155 34415.4883 174 34415.4883 6221.6155 175 13351.2167 34415.4883 176 -272740.1617 13351.2167 177 35257.0689 -272740.1617 178 34415.4883 35257.0689 179 34415.4883 34415.4883 180 40649.1911 34415.4883 181 35647.8860 40649.1911 182 -39381.2283 35647.8860 183 -11864.4442 -39381.2283 184 23409.1469 -11864.4442 185 22030.6439 23409.1469 186 13706.8251 22030.6439 187 18336.7011 13706.8251 188 5718.6132 18336.7011 189 34415.4883 5718.6132 190 16920.7901 34415.4883 191 34415.4883 16920.7901 192 16969.4330 34415.4883 193 34415.4883 16969.4330 194 26919.3675 34415.4883 195 34415.4883 26919.3675 196 -158975.0055 34415.4883 197 34415.4883 -158975.0055 198 34415.4883 34415.4883 199 -38659.2549 34415.4883 200 34415.4883 -38659.2549 201 34415.4883 34415.4883 202 34415.4883 34415.4883 203 -26837.7169 34415.4883 204 50961.7796 -26837.7169 205 23256.4882 50961.7796 206 40964.5220 23256.4882 207 34415.4883 40964.5220 208 -20082.2236 34415.4883 209 17885.0109 -20082.2236 210 34415.4883 17885.0109 211 34415.4883 34415.4883 212 14820.3119 34415.4883 213 -17066.9307 14820.3119 214 13196.7206 -17066.9307 215 29016.5034 13196.7206 216 18708.5506 29016.5034 217 17755.9585 18708.5506 218 29016.5034 17755.9585 219 -38566.9285 29016.5034 220 7229.3939 -38566.9285 221 34415.4883 7229.3939 222 31650.5575 34415.4883 223 -22793.3612 31650.5575 224 32064.3713 -22793.3612 225 -24644.3051 32064.3713 226 9291.0560 -24644.3051 227 34415.4883 9291.0560 228 34003.5446 34415.4883 229 34415.4883 34003.5446 230 3980.5053 34415.4883 231 16116.0560 3980.5053 232 34415.4883 16116.0560 233 -11450.5996 34415.4883 234 34415.4883 -11450.5996 235 34564.4225 34415.4883 236 34415.4883 34564.4225 237 18023.1699 34415.4883 238 34415.4883 18023.1699 239 34415.4883 34415.4883 240 34415.4883 34415.4883 241 24932.1768 34415.4883 242 34415.4883 24932.1768 243 14324.5227 34415.4883 244 59095.9927 14324.5227 245 34415.4883 59095.9927 246 34415.4883 34415.4883 247 -10983.1623 34415.4883 248 34415.4883 -10983.1623 249 34415.4883 34415.4883 250 -20338.8207 34415.4883 251 34415.4883 -20338.8207 252 29977.0971 34415.4883 253 -71041.4238 29977.0971 254 -179472.7167 -71041.4238 255 34415.4883 -179472.7167 256 204.6311 34415.4883 257 -66069.3918 204.6311 258 34415.4883 -66069.3918 259 34415.4883 34415.4883 260 34415.4883 34415.4883 261 34415.4883 34415.4883 262 45144.5464 34415.4883 263 21166.9312 45144.5464 264 34415.4883 21166.9312 265 34415.4883 34415.4883 266 34415.4883 34415.4883 267 102136.7639 34415.4883 268 34415.4883 102136.7639 269 -10376.3459 34415.4883 270 34415.4883 -10376.3459 271 34415.4883 34415.4883 272 34415.4883 34415.4883 273 -21956.9647 34415.4883 274 34415.4883 -21956.9647 275 9179.9993 34415.4883 276 -18784.9860 9179.9993 277 34415.4883 -18784.9860 278 -61114.5573 34415.4883 279 19816.8185 -61114.5573 280 34415.4883 19816.8185 281 34415.4883 34415.4883 282 -41767.5768 34415.4883 283 34415.4883 -41767.5768 284 41684.7152 34415.4883 285 34415.4883 41684.7152 286 24748.0259 34415.4883 287 31329.2606 24748.0259 288 34415.4883 31329.2606 289 -32424.7001 34415.4883 290 34415.4883 -32424.7001 291 14906.5245 34415.4883 292 59054.4943 14906.5245 293 34415.4883 59054.4943 294 -7392.6895 34415.4883 295 33610.4098 -7392.6895 296 34415.4883 33610.4098 297 34415.4883 34415.4883 298 34415.4883 34415.4883 299 10120.0560 34415.4883 300 21534.1915 10120.0560 301 14040.1758 21534.1915 302 31715.9958 14040.1758 303 34415.4883 31715.9958 304 27185.5923 34415.4883 305 2057.4240 27185.5923 306 34415.4883 2057.4240 307 -11513.0490 34415.4883 308 -24889.1373 -11513.0490 309 34415.4883 -24889.1373 310 34415.4883 34415.4883 311 1629.6336 34415.4883 312 34415.4883 1629.6336 313 7999.2017 34415.4883 314 13160.2732 7999.2017 315 56646.4587 13160.2732 316 32664.8879 56646.4587 317 -37271.6918 32664.8879 318 20908.3556 -37271.6918 319 34415.4883 20908.3556 320 34415.4883 34415.4883 321 34415.4883 34415.4883 322 34415.4883 34415.4883 323 34415.4883 34415.4883 324 34415.4883 34415.4883 325 28841.5811 34415.4883 326 34415.4883 28841.5811 327 34415.4883 34415.4883 328 34415.4883 34415.4883 329 34415.4883 34415.4883 330 34415.4883 34415.4883 331 34415.4883 34415.4883 332 34415.4883 34415.4883 333 21386.4987 34415.4883 334 31736.5444 21386.4987 335 23608.5184 31736.5444 336 -38566.9446 23608.5184 337 34415.4883 -38566.9446 338 30531.5027 34415.4883 339 8242.3020 30531.5027 340 -14503.0610 8242.3020 341 34415.4883 -14503.0610 342 34415.4883 34415.4883 343 34415.4883 34415.4883 344 -7801.9139 34415.4883 345 34415.4883 -7801.9139 346 11840.1331 34415.4883 347 34415.4883 11840.1331 348 34415.4883 34415.4883 349 11409.5635 34415.4883 350 25952.7483 11409.5635 351 14838.1667 25952.7483 352 16032.3686 14838.1667 353 32957.1139 16032.3686 354 34415.4883 32957.1139 355 34415.4883 34415.4883 356 25566.6364 34415.4883 357 34415.4883 25566.6364 358 34415.4883 34415.4883 359 34415.4883 34415.4883 360 23638.9203 34415.4883 361 34415.4883 23638.9203 362 34415.4883 34415.4883 363 -47363.1526 34415.4883 364 9446.3909 -47363.1526 365 34415.4883 9446.3909 366 34415.4883 34415.4883 367 34415.4883 34415.4883 368 34415.4883 34415.4883 369 15692.0410 34415.4883 370 34415.4883 15692.0410 371 21988.9815 34415.4883 372 62403.2924 21988.9815 373 25460.1988 62403.2924 374 34415.4883 25460.1988 375 34415.4883 34415.4883 376 -22042.5868 34415.4883 377 22022.0219 -22042.5868 378 11904.5635 22022.0219 379 34415.4883 11904.5635 380 34415.4883 34415.4883 381 3918.4103 34415.4883 382 -65207.6130 3918.4103 383 34415.4883 -65207.6130 384 6855.1773 34415.4883 385 -114489.6195 6855.1773 386 36174.3426 -114489.6195 387 8432.9664 36174.3426 388 24589.4649 8432.9664 389 1991.2370 24589.4649 390 17387.4017 1991.2370 391 -7312.9789 17387.4017 392 -38921.3620 -7312.9789 393 82529.6201 -38921.3620 394 14101.0516 82529.6201 395 -22547.5156 14101.0516 396 -115634.5973 -22547.5156 397 -819140.4195 -115634.5973 398 1469.0903 -819140.4195 399 -183653.8205 1469.0903 400 31975.9930 -183653.8205 401 36485.6108 31975.9930 402 -2751.3511 36485.6108 403 -48468.7720 -2751.3511 404 -17250.2209 -48468.7720 405 31178.7587 -17250.2209 406 17219.3911 31178.7587 407 -11410.3112 17219.3911 408 -10879.3322 -11410.3112 409 -41168.2135 -10879.3322 410 -124676.3668 -41168.2135 411 -47696.4693 -124676.3668 412 -56343.7006 -47696.4693 413 93154.6182 -56343.7006 414 -117449.2482 93154.6182 415 -59536.7217 -117449.2482 416 35577.8986 -59536.7217 417 -42581.9122 35577.8986 418 -61195.8229 -42581.9122 419 -27801.2093 -61195.8229 420 -13972.4958 -27801.2093 421 -245082.7760 -13972.4958 422 -528552.8930 -245082.7760 423 -383023.0073 -528552.8930 424 -68961.9405 -383023.0073 425 -291280.2801 -68961.9405 426 -306083.6139 -291280.2801 427 -1177595.1299 -306083.6139 428 -97935.0948 -1177595.1299 429 -153965.3588 -97935.0948 430 -464829.1094 -153965.3588 > 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/7jljo1292956049.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/8jljo1292956049.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/9jljo1292956049.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/10jljo1292956049.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/114l0c1292956049.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/12bn0y1292956050.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/13oxg71292956050.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/14h6xs1292956050.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/15dydi1292956050.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/16hyb61292956050.tab") + } > > try(system("convert tmp/1nb4g1292956049.ps tmp/1nb4g1292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/2nb4g1292956049.ps tmp/2nb4g1292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/3nb4g1292956049.ps tmp/3nb4g1292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/4xk311292956049.ps tmp/4xk311292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/5xk311292956049.ps tmp/5xk311292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/68u241292956049.ps tmp/68u241292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/7jljo1292956049.ps tmp/7jljo1292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/8jljo1292956049.ps tmp/8jljo1292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/9jljo1292956049.ps tmp/9jljo1292956049.png",intern=TRUE)) character(0) > try(system("convert tmp/10jljo1292956049.ps tmp/10jljo1292956049.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 11.360 1.530 12.892