R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-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(210907 + ,56 + ,396 + ,3 + ,30 + ,112285 + ,1 + ,120982 + ,56 + ,297 + ,4 + ,28 + ,84786 + ,1 + ,176508 + ,54 + ,559 + ,12 + ,38 + ,83123 + ,1 + ,179321 + ,89 + ,967 + ,2 + ,30 + ,101193 + ,1 + ,123185 + ,40 + ,270 + ,1 + ,22 + ,38361 + ,1 + ,52746 + ,25 + ,143 + ,3 + ,26 + ,68504 + ,1 + ,385534 + ,92 + ,1562 + ,0 + ,25 + ,119182 + ,1 + ,33170 + ,18 + ,109 + ,0 + ,18 + ,22807 + ,1 + ,101645 + ,63 + ,371 + ,0 + ,11 + ,17140 + ,0 + ,149061 + ,44 + ,656 + ,5 + ,26 + ,116174 + ,1 + ,165446 + ,33 + ,511 + ,0 + ,25 + ,57635 + ,1 + ,237213 + ,84 + ,655 + ,0 + ,38 + ,66198 + ,1 + ,173326 + ,88 + ,465 + ,7 + ,44 + ,71701 + ,1 + ,133131 + ,55 + ,525 + ,7 + ,30 + ,57793 + ,1 + ,258873 + ,60 + ,885 + ,3 + ,40 + ,80444 + ,1 + ,180083 + ,66 + ,497 + ,9 + ,34 + ,53855 + ,1 + ,324799 + ,154 + ,1436 + ,0 + ,47 + ,97668 + ,1 + ,230964 + ,53 + ,612 + ,4 + ,30 + ,133824 + ,1 + ,236785 + ,119 + ,865 + ,3 + ,31 + ,101481 + ,1 + ,135473 + ,41 + ,385 + ,0 + ,23 + ,99645 + ,1 + ,202925 + ,61 + ,567 + ,7 + ,36 + ,114789 + ,1 + ,215147 + ,58 + ,639 + ,0 + ,36 + ,99052 + ,1 + ,344297 + ,75 + ,963 + ,1 + ,30 + ,67654 + ,1 + ,153935 + ,33 + ,398 + ,5 + ,25 + ,65553 + ,1 + ,132943 + ,40 + ,410 + ,7 + ,39 + ,97500 + ,1 + ,174724 + ,92 + ,966 + ,0 + ,34 + ,69112 + ,1 + ,174415 + ,100 + ,801 + ,0 + ,31 + ,82753 + ,1 + ,225548 + ,112 + ,892 + ,5 + ,31 + ,85323 + ,1 + ,223632 + ,73 + ,513 + ,0 + ,33 + ,72654 + ,1 + ,124817 + ,40 + ,469 + ,0 + ,25 + ,30727 + ,1 + ,221698 + ,45 + ,683 + ,0 + ,33 + ,77873 + ,1 + ,210767 + ,60 + ,643 + ,3 + ,35 + ,117478 + ,1 + ,170266 + ,62 + ,535 + ,4 + ,42 + ,74007 + ,1 + ,260561 + ,75 + ,625 + ,1 + ,43 + ,90183 + ,1 + ,84853 + ,31 + ,264 + ,4 + ,30 + ,61542 + ,1 + ,294424 + ,77 + ,992 + ,2 + ,33 + ,101494 + ,1 + ,101011 + ,34 + ,238 + ,0 + ,13 + ,27570 + ,1 + ,215641 + ,46 + ,818 + ,0 + ,32 + ,55813 + ,1 + ,325107 + ,99 + ,937 + ,0 + ,36 + ,79215 + ,1 + ,7176 + ,17 + ,70 + ,0 + ,0 + ,1423 + ,1 + ,167542 + ,66 + ,507 + ,2 + ,28 + ,55461 + ,1 + ,106408 + ,30 + ,260 + ,1 + ,14 + ,31081 + ,1 + ,96560 + ,76 + ,503 + ,0 + ,17 + ,22996 + ,1 + ,265769 + ,146 + ,927 + ,2 + ,32 + ,83122 + ,1 + ,269651 + ,67 + ,1269 + ,10 + ,30 + ,70106 + ,1 + ,149112 + ,56 + ,537 + ,6 + ,35 + ,60578 + ,1 + ,175824 + ,107 + ,910 + ,0 + ,20 + ,39992 + ,1 + ,152871 + ,58 + ,532 + ,5 + ,28 + ,79892 + ,1 + ,111665 + ,34 + ,345 + ,4 + ,28 + ,49810 + ,1 + ,116408 + ,61 + ,918 + ,1 + ,39 + ,71570 + ,1 + ,362301 + ,119 + ,1635 + ,2 + ,34 + ,100708 + ,1 + ,78800 + ,42 + ,330 + ,2 + ,26 + ,33032 + ,1 + ,183167 + ,66 + ,557 + ,0 + ,39 + ,82875 + ,1 + ,277965 + ,89 + ,1178 + ,8 + ,39 + ,139077 + ,1 + ,150629 + ,44 + ,740 + ,3 + ,33 + ,71595 + ,1 + ,168809 + ,66 + ,452 + ,0 + ,28 + ,72260 + ,1 + ,24188 + ,24 + ,218 + ,0 + ,4 + ,5950 + ,1 + ,329267 + ,259 + ,764 + ,8 + ,39 + ,115762 + ,1 + ,65029 + ,17 + ,255 + ,5 + ,18 + ,32551 + ,1 + ,101097 + ,64 + ,454 + ,3 + ,14 + ,31701 + ,1 + ,218946 + ,41 + ,866 + ,1 + ,29 + ,80670 + ,1 + ,244052 + ,68 + ,574 + ,5 + ,44 + ,143558 + ,1 + ,341570 + ,168 + ,1276 + ,1 + ,21 + ,117105 + ,1 + ,103597 + ,43 + ,379 + ,1 + ,16 + ,23789 + ,1 + ,233328 + ,132 + ,825 + ,5 + ,28 + ,120733 + ,1 + ,256462 + ,105 + ,798 + ,0 + ,35 + ,105195 + ,1 + ,206161 + ,71 + ,663 + ,12 + ,28 + ,73107 + ,1 + ,311473 + ,112 + ,1069 + ,8 + ,38 + ,132068 + ,1 + ,235800 + ,94 + ,921 + ,8 + ,23 + ,149193 + ,1 + ,177939 + ,82 + ,858 + ,8 + ,36 + ,46821 + ,1 + ,207176 + ,70 + ,711 + ,8 + ,32 + ,87011 + ,1 + ,196553 + ,57 + ,503 + ,2 + ,29 + ,95260 + ,1 + ,174184 + ,53 + ,382 + ,0 + ,25 + ,55183 + ,1 + ,143246 + ,103 + ,464 + ,5 + ,27 + ,106671 + ,1 + ,187559 + ,121 + ,717 + ,8 + ,36 + ,73511 + ,1 + ,187681 + ,62 + ,690 + ,2 + ,28 + ,92945 + ,1 + ,119016 + ,52 + ,462 + ,5 + ,23 + ,78664 + ,1 + ,182192 + ,52 + ,657 + ,12 + ,40 + ,70054 + ,1 + ,73566 + ,32 + ,385 + ,6 + ,23 + ,22618 + ,1 + ,194979 + ,62 + ,577 + ,7 + ,40 + ,74011 + ,1 + ,167488 + ,45 + ,619 + ,2 + ,28 + ,83737 + ,1 + ,143756 + ,46 + ,479 + ,0 + ,34 + ,69094 + ,1 + ,275541 + ,63 + ,817 + ,4 + ,33 + ,93133 + ,1 + ,243199 + ,75 + ,752 + ,3 + ,28 + ,95536 + ,1 + ,182999 + ,88 + ,430 + ,6 + ,34 + ,225920 + ,1 + ,135649 + ,46 + ,451 + ,2 + ,30 + ,62133 + ,1 + ,152299 + ,53 + ,537 + ,0 + ,33 + ,61370 + ,1 + ,120221 + ,37 + ,519 + ,1 + ,22 + ,43836 + ,1 + ,346485 + ,90 + ,1000 + ,0 + ,38 + ,106117 + ,1 + ,145790 + ,63 + ,637 + ,5 + ,26 + ,38692 + ,1 + ,193339 + ,78 + ,465 + ,2 + ,35 + ,84651 + ,1 + ,80953 + ,25 + ,437 + ,0 + ,8 + ,56622 + ,1 + ,122774 + ,45 + ,711 + ,0 + ,24 + ,15986 + ,1 + ,130585 + ,46 + ,299 + ,5 + ,29 + ,95364 + ,1 + ,112611 + ,41 + ,248 + ,0 + ,20 + ,26706 + ,0 + ,286468 + ,144 + ,1162 + ,1 + ,29 + ,89691 + ,1 + ,241066 + ,82 + ,714 + ,0 + ,45 + ,67267 + ,1 + ,148446 + ,91 + ,905 + ,1 + ,37 + ,126846 + ,1 + ,204713 + ,71 + ,649 + ,1 + ,33 + ,41140 + ,1 + ,182079 + ,63 + ,512 + ,2 + ,33 + ,102860 + ,1 + ,140344 + ,53 + ,472 + ,6 + ,25 + ,51715 + ,1 + ,220516 + ,62 + ,905 + ,1 + ,32 + ,55801 + ,1 + ,243060 + ,63 + ,786 + ,4 + ,29 + ,111813 + ,1 + ,162765 + ,32 + ,489 + ,2 + ,28 + ,120293 + ,1 + ,182613 + ,39 + ,479 + ,3 + ,28 + ,138599 + ,1 + ,232138 + ,62 + ,617 + ,0 + ,31 + ,161647 + ,1 + ,265318 + ,117 + ,925 + ,10 + ,52 + ,115929 + ,1 + ,85574 + ,34 + ,351 + ,0 + ,21 + ,24266 + ,0 + ,310839 + ,92 + ,1144 + ,9 + ,24 + ,162901 + ,1 + ,225060 + ,93 + ,669 + ,7 + ,41 + ,109825 + ,1 + ,232317 + ,54 + ,707 + ,0 + ,33 + ,129838 + ,1 + ,144966 + ,144 + ,458 + ,0 + ,32 + ,37510 + ,1 + ,43287 + ,14 + ,214 + ,4 + ,19 + ,43750 + ,1 + ,155754 + ,61 + ,599 + ,4 + ,20 + ,40652 + ,1 + ,164709 + ,109 + ,572 + ,0 + ,31 + ,87771 + ,1 + ,201940 + ,38 + ,897 + ,0 + ,31 + ,85872 + ,1 + ,235454 + ,73 + ,819 + ,0 + ,32 + ,89275 + ,1 + ,220801 + ,75 + ,720 + ,1 + ,18 + ,44418 + ,0 + ,99466 + ,50 + ,273 + ,0 + ,23 + ,192565 + ,1 + ,92661 + ,61 + ,508 + ,1 + ,17 + ,35232 + ,0 + ,133328 + ,55 + ,506 + ,0 + ,20 + ,40909 + ,0 + ,61361 + ,77 + ,451 + ,0 + ,12 + ,13294 + ,0 + ,125930 + ,75 + ,699 + ,4 + ,17 + ,32387 + ,0 + ,100750 + ,72 + ,407 + ,0 + ,30 + ,140867 + ,1 + ,224549 + ,50 + ,465 + ,4 + ,31 + ,120662 + ,1 + ,82316 + ,32 + ,245 + ,4 + ,10 + ,21233 + ,0 + ,102010 + ,53 + ,370 + ,3 + ,13 + ,44332 + ,0 + ,101523 + ,42 + ,316 + ,0 + ,22 + ,61056 + ,0 + ,243511 + ,71 + ,603 + ,0 + ,42 + ,101338 + ,1 + ,22938 + ,10 + ,154 + ,0 + ,1 + ,1168 + ,1 + ,41566 + ,35 + ,229 + ,5 + ,9 + ,13497 + ,0 + ,152474 + ,65 + ,577 + ,0 + ,32 + ,65567 + ,1 + ,61857 + ,25 + ,192 + ,4 + ,11 + ,25162 + ,1 + ,99923 + ,66 + ,617 + ,0 + ,25 + ,32334 + ,0 + ,132487 + ,41 + ,411 + ,0 + ,36 + ,40735 + ,1 + ,317394 + ,86 + ,975 + ,1 + ,31 + ,91413 + ,1 + ,21054 + ,16 + ,146 + ,0 + ,0 + ,855 + ,1 + ,209641 + ,42 + ,705 + ,5 + ,24 + ,97068 + ,1 + ,22648 + ,19 + ,184 + ,0 + ,13 + ,44339 + ,0 + ,31414 + ,19 + ,200 + ,0 + ,8 + ,14116 + ,1 + ,46698 + ,45 + ,274 + ,0 + ,13 + ,10288 + ,0 + ,131698 + ,65 + ,502 + ,0 + ,19 + ,65622 + ,0 + ,91735 + ,35 + ,382 + ,0 + ,18 + ,16563 + ,0 + ,244749 + ,95 + ,964 + ,2 + ,33 + ,76643 + ,1 + ,184510 + ,49 + ,537 + ,7 + ,40 + ,110681 + ,1 + ,79863 + ,37 + ,438 + ,1 + ,22 + ,29011 + ,0 + ,128423 + ,64 + ,369 + ,8 + ,38 + ,92696 + ,1 + ,97839 + ,38 + ,417 + ,2 + ,24 + ,94785 + ,1 + ,38214 + ,34 + ,276 + ,0 + ,8 + ,8773 + ,1 + ,151101 + ,32 + ,514 + ,2 + ,35 + ,83209 + ,1 + ,272458 + ,65 + ,822 + ,0 + ,43 + ,93815 + ,1 + ,172494 + ,52 + ,389 + ,0 + ,43 + ,86687 + ,1 + ,108043 + ,62 + ,466 + ,1 + ,14 + ,34553 + ,0 + ,328107 + ,65 + ,1255 + ,3 + ,41 + ,105547 + ,1 + ,250579 + ,83 + ,694 + ,0 + ,38 + ,103487 + ,1 + ,351067 + ,95 + ,1024 + ,3 + ,45 + ,213688 + ,1 + ,158015 + ,29 + ,400 + ,0 + ,31 + ,71220 + ,1 + ,98866 + ,18 + ,397 + ,0 + ,13 + ,23517 + ,0 + ,85439 + ,33 + ,350 + ,0 + ,28 + ,56926 + ,1 + ,229242 + ,247 + ,719 + ,4 + ,31 + ,91721 + ,1 + ,351619 + ,139 + ,1277 + ,4 + ,40 + ,115168 + ,1 + ,84207 + ,29 + ,356 + ,11 + ,30 + ,111194 + ,1 + ,120445 + ,118 + ,457 + ,0 + ,16 + ,51009 + ,0 + ,324598 + ,110 + ,1402 + ,0 + ,37 + ,135777 + ,1 + ,131069 + ,67 + ,600 + ,4 + ,30 + ,51513 + ,1 + ,204271 + ,42 + ,480 + ,0 + ,35 + ,74163 + ,1 + ,165543 + ,65 + ,595 + ,1 + ,32 + ,51633 + ,1 + ,141722 + ,94 + ,436 + ,0 + ,27 + ,75345 + ,1 + ,116048 + ,64 + ,230 + ,0 + ,20 + ,33416 + ,0 + ,250047 + ,81 + ,651 + ,0 + ,18 + ,83305 + ,0 + ,299775 + ,95 + ,1367 + ,9 + ,31 + ,98952 + ,1 + ,195838 + ,67 + ,564 + ,1 + ,31 + ,102372 + ,1 + ,173260 + ,63 + ,716 + ,3 + ,21 + ,37238 + ,1 + ,254488 + ,83 + ,747 + ,10 + ,39 + ,103772 + ,1 + ,104389 + ,45 + ,467 + ,5 + ,41 + ,123969 + ,1 + ,136084 + ,30 + ,671 + ,0 + ,13 + ,27142 + ,0 + ,199476 + ,70 + ,861 + ,2 + ,32 + ,135400 + ,1 + ,92499 + ,32 + ,319 + ,0 + ,18 + ,21399 + ,0 + ,224330 + ,83 + ,612 + ,1 + ,39 + ,130115 + ,1 + ,135781 + ,31 + ,433 + ,2 + ,14 + ,24874 + ,0 + ,74408 + ,67 + ,434 + ,4 + ,7 + ,34988 + ,0 + ,81240 + ,66 + ,503 + ,0 + ,17 + ,45549 + ,0 + ,14688 + ,10 + ,85 + ,0 + ,0 + ,6023 + ,1 + ,181633 + ,70 + ,564 + ,2 + ,30 + ,64466 + ,1 + ,271856 + ,103 + ,824 + ,1 + ,37 + ,54990 + ,1 + ,7199 + ,5 + ,74 + ,0 + ,0 + ,1644 + ,1 + ,46660 + ,20 + ,259 + ,0 + ,5 + ,6179 + ,1 + ,17547 + ,5 + ,69 + ,0 + ,1 + ,3926 + ,1 + ,133368 + ,36 + ,535 + ,1 + ,16 + ,32755 + ,0 + ,95227 + ,34 + ,239 + ,0 + ,32 + ,34777 + ,1 + ,152601 + ,48 + ,438 + ,2 + ,24 + ,73224 + ,1 + ,98146 + ,40 + ,459 + ,0 + ,17 + ,27114 + ,0 + ,79619 + ,43 + ,426 + ,3 + ,11 + ,20760 + ,0 + ,59194 + ,31 + ,288 + ,6 + ,24 + ,37636 + ,0 + ,139942 + ,42 + ,498 + ,0 + ,22 + ,65461 + ,0 + ,118612 + ,46 + ,454 + ,2 + ,12 + ,30080 + ,0 + ,72880 + ,33 + ,376 + ,0 + ,19 + ,24094 + ,0 + ,65475 + ,18 + ,225 + ,2 + ,13 + ,69008 + ,0 + ,99643 + ,55 + ,555 + ,1 + ,17 + ,54968 + ,0 + ,71965 + ,35 + ,252 + ,1 + ,15 + ,46090 + ,0 + ,77272 + ,59 + ,208 + ,2 + ,16 + ,27507 + ,0 + ,49289 + ,19 + ,130 + ,1 + ,24 + ,10672 + ,0 + ,135131 + ,66 + ,481 + ,0 + ,15 + ,34029 + ,0 + ,108446 + ,60 + ,389 + ,1 + ,17 + ,46300 + ,0 + ,89746 + ,36 + ,565 + ,3 + ,18 + ,24760 + ,0 + ,44296 + ,25 + ,173 + ,0 + ,20 + ,18779 + ,0 + ,77648 + ,47 + ,278 + ,0 + ,16 + ,21280 + ,0 + ,181528 + ,54 + ,609 + ,0 + ,16 + ,40662 + ,0 + ,134019 + ,53 + ,422 + ,0 + ,18 + ,28987 + ,0 + ,124064 + ,40 + ,445 + ,1 + ,22 + ,22827 + ,0 + ,92630 + ,40 + ,387 + ,4 + ,8 + ,18513 + ,0 + ,121848 + ,39 + ,339 + ,0 + ,17 + ,30594 + ,0 + ,52915 + ,14 + ,181 + ,0 + ,18 + ,24006 + ,0 + ,81872 + ,45 + ,245 + ,0 + ,16 + ,27913 + ,0 + ,58981 + ,36 + ,384 + ,7 + ,23 + ,42744 + ,0 + ,53515 + ,28 + ,212 + ,2 + ,22 + ,12934 + ,0 + ,60812 + ,44 + ,399 + ,0 + ,13 + ,22574 + ,0 + ,56375 + ,30 + ,229 + ,7 + ,13 + ,41385 + ,0 + ,65490 + ,22 + ,224 + ,3 + ,16 + ,18653 + ,0 + ,80949 + ,17 + ,203 + ,0 + ,16 + ,18472 + ,0 + ,76302 + ,31 + ,333 + ,0 + ,20 + ,30976 + ,0 + ,104011 + ,55 + ,384 + ,6 + ,22 + ,63339 + ,0 + ,98104 + ,54 + ,636 + ,2 + ,17 + ,25568 + ,0 + ,67989 + ,21 + ,185 + ,0 + ,18 + ,33747 + ,0 + ,30989 + ,14 + ,93 + ,0 + ,17 + ,4154 + ,0 + ,135458 + ,81 + ,581 + ,3 + ,12 + ,19474 + ,0 + ,73504 + ,35 + ,248 + ,0 + ,7 + ,35130 + ,0 + ,63123 + ,43 + ,304 + ,1 + ,17 + ,39067 + ,0 + ,61254 + ,46 + ,344 + ,1 + ,14 + ,13310 + ,0 + ,74914 + ,30 + ,407 + ,0 + ,23 + ,65892 + ,0 + ,31774 + ,23 + ,170 + ,1 + ,17 + ,4143 + ,0 + ,81437 + ,38 + ,312 + ,0 + ,14 + ,28579 + ,0 + ,87186 + ,54 + ,507 + ,0 + ,15 + ,51776 + ,0 + ,50090 + ,20 + ,224 + ,0 + ,17 + ,21152 + ,0 + ,65745 + ,53 + ,340 + ,0 + ,21 + ,38084 + ,0 + ,56653 + ,45 + ,168 + ,0 + ,18 + ,27717 + ,0 + ,158399 + ,39 + ,443 + ,0 + ,18 + ,32928 + ,0 + ,46455 + ,20 + ,204 + ,0 + ,17 + ,11342 + ,0 + ,73624 + ,24 + ,367 + ,0 + ,17 + ,19499 + ,0 + ,38395 + ,31 + ,210 + ,0 + ,16 + ,16380 + ,0 + ,91899 + ,35 + ,335 + ,0 + ,15 + ,36874 + ,0 + ,139526 + ,151 + ,364 + ,0 + ,21 + ,48259 + ,0 + ,52164 + ,52 + ,178 + ,0 + ,16 + ,16734 + ,0 + ,51567 + ,30 + ,206 + ,2 + ,14 + ,28207 + ,0 + ,70551 + ,31 + ,279 + ,0 + ,15 + ,30143 + ,0 + ,84856 + ,29 + ,387 + ,1 + ,17 + ,41369 + ,0 + ,102538 + ,57 + ,490 + ,1 + ,15 + ,45833 + ,0 + ,86678 + ,40 + ,238 + ,0 + ,15 + ,29156 + ,0 + ,85709 + ,44 + ,343 + ,0 + ,10 + ,35944 + ,0 + ,34662 + ,25 + ,232 + ,0 + ,6 + ,36278 + ,0 + ,150580 + ,77 + ,530 + ,0 + ,22 + ,45588 + ,0 + ,99611 + ,35 + ,291 + ,0 + ,21 + ,45097 + ,0 + ,19349 + ,11 + ,67 + ,0 + ,1 + ,3895 + ,0 + ,99373 + ,63 + ,397 + ,1 + ,18 + ,28394 + ,0 + ,86230 + ,44 + ,467 + ,0 + ,17 + ,18632 + ,0 + ,30837 + ,19 + ,178 + ,0 + ,4 + ,2325 + ,0 + ,31706 + ,13 + ,175 + ,0 + ,10 + ,25139 + ,0 + ,89806 + ,42 + ,299 + ,0 + ,16 + ,27975 + ,0 + ,62088 + ,38 + ,154 + ,1 + ,16 + ,14483 + ,0 + ,40151 + ,29 + ,106 + ,0 + ,9 + ,13127 + ,0 + ,27634 + ,20 + ,189 + ,0 + ,16 + ,5839 + ,0 + ,76990 + ,27 + ,194 + ,0 + ,17 + ,24069 + ,0 + ,37460 + ,20 + ,135 + ,0 + ,7 + ,3738 + ,0 + ,54157 + ,19 + ,201 + ,0 + ,15 + ,18625 + ,0 + ,49862 + ,37 + ,207 + ,0 + ,14 + ,36341 + ,0 + ,84337 + ,26 + ,280 + ,0 + ,14 + ,24548 + ,0 + ,64175 + ,42 + ,260 + ,0 + ,18 + ,21792 + ,0 + ,59382 + ,49 + ,227 + ,0 + ,12 + ,26263 + ,0 + ,119308 + ,30 + ,239 + ,0 + ,16 + ,23686 + ,0 + ,76702 + ,49 + ,333 + ,0 + ,21 + ,49303 + ,0 + ,103425 + ,67 + ,428 + ,1 + ,19 + ,25659 + ,0 + ,70344 + ,28 + ,230 + ,0 + ,16 + ,28904 + ,0 + ,43410 + ,19 + ,292 + ,0 + ,1 + ,2781 + ,0 + ,104838 + ,49 + ,350 + ,1 + ,16 + ,29236 + ,0 + ,62215 + ,27 + ,186 + ,0 + ,10 + ,19546 + ,0 + ,69304 + ,30 + ,326 + ,6 + ,19 + ,22818 + ,0 + ,53117 + ,22 + ,155 + ,3 + ,12 + ,32689 + ,0 + ,19764 + ,12 + ,75 + ,1 + ,2 + ,5752 + ,0 + ,86680 + ,31 + ,361 + ,2 + ,14 + ,22197 + ,0 + ,84105 + ,20 + ,261 + ,0 + ,17 + ,20055 + ,0 + ,77945 + ,20 + ,299 + ,0 + ,19 + ,25272 + ,0 + ,89113 + ,39 + ,300 + ,0 + ,14 + ,82206 + ,0 + ,91005 + ,29 + ,450 + ,3 + ,11 + ,32073 + ,0 + ,40248 + ,16 + ,183 + ,1 + ,4 + ,5444 + ,0 + ,64187 + ,27 + ,238 + ,0 + ,16 + ,20154 + ,0 + ,50857 + ,21 + ,165 + ,0 + ,20 + ,36944 + ,0 + ,56613 + ,19 + ,234 + ,1 + ,12 + ,8019 + ,0 + ,62792 + ,35 + ,176 + ,0 + ,15 + ,30884 + ,0 + ,72535 + ,14 + ,329 + ,0 + ,16 + ,19540 + ,0) + ,dim=c(7 + ,289) + ,dimnames=list(c('Time' + ,'Logins' + ,'Views' + ,'Shared' + ,'Reviewed' + ,'TotalSize' + ,'Pop') + ,1:289)) > y <- array(NA,dim=c(7,289),dimnames=list(c('Time','Logins','Views','Shared','Reviewed','TotalSize','Pop'),1:289)) > 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 = '1' > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from 'package:base': as.Date, as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Time Logins Views Shared Reviewed TotalSize Pop 1 210907 56 396 3 30 112285 1 2 120982 56 297 4 28 84786 1 3 176508 54 559 12 38 83123 1 4 179321 89 967 2 30 101193 1 5 123185 40 270 1 22 38361 1 6 52746 25 143 3 26 68504 1 7 385534 92 1562 0 25 119182 1 8 33170 18 109 0 18 22807 1 9 101645 63 371 0 11 17140 0 10 149061 44 656 5 26 116174 1 11 165446 33 511 0 25 57635 1 12 237213 84 655 0 38 66198 1 13 173326 88 465 7 44 71701 1 14 133131 55 525 7 30 57793 1 15 258873 60 885 3 40 80444 1 16 180083 66 497 9 34 53855 1 17 324799 154 1436 0 47 97668 1 18 230964 53 612 4 30 133824 1 19 236785 119 865 3 31 101481 1 20 135473 41 385 0 23 99645 1 21 202925 61 567 7 36 114789 1 22 215147 58 639 0 36 99052 1 23 344297 75 963 1 30 67654 1 24 153935 33 398 5 25 65553 1 25 132943 40 410 7 39 97500 1 26 174724 92 966 0 34 69112 1 27 174415 100 801 0 31 82753 1 28 225548 112 892 5 31 85323 1 29 223632 73 513 0 33 72654 1 30 124817 40 469 0 25 30727 1 31 221698 45 683 0 33 77873 1 32 210767 60 643 3 35 117478 1 33 170266 62 535 4 42 74007 1 34 260561 75 625 1 43 90183 1 35 84853 31 264 4 30 61542 1 36 294424 77 992 2 33 101494 1 37 101011 34 238 0 13 27570 1 38 215641 46 818 0 32 55813 1 39 325107 99 937 0 36 79215 1 40 7176 17 70 0 0 1423 1 41 167542 66 507 2 28 55461 1 42 106408 30 260 1 14 31081 1 43 96560 76 503 0 17 22996 1 44 265769 146 927 2 32 83122 1 45 269651 67 1269 10 30 70106 1 46 149112 56 537 6 35 60578 1 47 175824 107 910 0 20 39992 1 48 152871 58 532 5 28 79892 1 49 111665 34 345 4 28 49810 1 50 116408 61 918 1 39 71570 1 51 362301 119 1635 2 34 100708 1 52 78800 42 330 2 26 33032 1 53 183167 66 557 0 39 82875 1 54 277965 89 1178 8 39 139077 1 55 150629 44 740 3 33 71595 1 56 168809 66 452 0 28 72260 1 57 24188 24 218 0 4 5950 1 58 329267 259 764 8 39 115762 1 59 65029 17 255 5 18 32551 1 60 101097 64 454 3 14 31701 1 61 218946 41 866 1 29 80670 1 62 244052 68 574 5 44 143558 1 63 341570 168 1276 1 21 117105 1 64 103597 43 379 1 16 23789 1 65 233328 132 825 5 28 120733 1 66 256462 105 798 0 35 105195 1 67 206161 71 663 12 28 73107 1 68 311473 112 1069 8 38 132068 1 69 235800 94 921 8 23 149193 1 70 177939 82 858 8 36 46821 1 71 207176 70 711 8 32 87011 1 72 196553 57 503 2 29 95260 1 73 174184 53 382 0 25 55183 1 74 143246 103 464 5 27 106671 1 75 187559 121 717 8 36 73511 1 76 187681 62 690 2 28 92945 1 77 119016 52 462 5 23 78664 1 78 182192 52 657 12 40 70054 1 79 73566 32 385 6 23 22618 1 80 194979 62 577 7 40 74011 1 81 167488 45 619 2 28 83737 1 82 143756 46 479 0 34 69094 1 83 275541 63 817 4 33 93133 1 84 243199 75 752 3 28 95536 1 85 182999 88 430 6 34 225920 1 86 135649 46 451 2 30 62133 1 87 152299 53 537 0 33 61370 1 88 120221 37 519 1 22 43836 1 89 346485 90 1000 0 38 106117 1 90 145790 63 637 5 26 38692 1 91 193339 78 465 2 35 84651 1 92 80953 25 437 0 8 56622 1 93 122774 45 711 0 24 15986 1 94 130585 46 299 5 29 95364 1 95 112611 41 248 0 20 26706 0 96 286468 144 1162 1 29 89691 1 97 241066 82 714 0 45 67267 1 98 148446 91 905 1 37 126846 1 99 204713 71 649 1 33 41140 1 100 182079 63 512 2 33 102860 1 101 140344 53 472 6 25 51715 1 102 220516 62 905 1 32 55801 1 103 243060 63 786 4 29 111813 1 104 162765 32 489 2 28 120293 1 105 182613 39 479 3 28 138599 1 106 232138 62 617 0 31 161647 1 107 265318 117 925 10 52 115929 1 108 85574 34 351 0 21 24266 0 109 310839 92 1144 9 24 162901 1 110 225060 93 669 7 41 109825 1 111 232317 54 707 0 33 129838 1 112 144966 144 458 0 32 37510 1 113 43287 14 214 4 19 43750 1 114 155754 61 599 4 20 40652 1 115 164709 109 572 0 31 87771 1 116 201940 38 897 0 31 85872 1 117 235454 73 819 0 32 89275 1 118 220801 75 720 1 18 44418 0 119 99466 50 273 0 23 192565 1 120 92661 61 508 1 17 35232 0 121 133328 55 506 0 20 40909 0 122 61361 77 451 0 12 13294 0 123 125930 75 699 4 17 32387 0 124 100750 72 407 0 30 140867 1 125 224549 50 465 4 31 120662 1 126 82316 32 245 4 10 21233 0 127 102010 53 370 3 13 44332 0 128 101523 42 316 0 22 61056 0 129 243511 71 603 0 42 101338 1 130 22938 10 154 0 1 1168 1 131 41566 35 229 5 9 13497 0 132 152474 65 577 0 32 65567 1 133 61857 25 192 4 11 25162 1 134 99923 66 617 0 25 32334 0 135 132487 41 411 0 36 40735 1 136 317394 86 975 1 31 91413 1 137 21054 16 146 0 0 855 1 138 209641 42 705 5 24 97068 1 139 22648 19 184 0 13 44339 0 140 31414 19 200 0 8 14116 1 141 46698 45 274 0 13 10288 0 142 131698 65 502 0 19 65622 0 143 91735 35 382 0 18 16563 0 144 244749 95 964 2 33 76643 1 145 184510 49 537 7 40 110681 1 146 79863 37 438 1 22 29011 0 147 128423 64 369 8 38 92696 1 148 97839 38 417 2 24 94785 1 149 38214 34 276 0 8 8773 1 150 151101 32 514 2 35 83209 1 151 272458 65 822 0 43 93815 1 152 172494 52 389 0 43 86687 1 153 108043 62 466 1 14 34553 0 154 328107 65 1255 3 41 105547 1 155 250579 83 694 0 38 103487 1 156 351067 95 1024 3 45 213688 1 157 158015 29 400 0 31 71220 1 158 98866 18 397 0 13 23517 0 159 85439 33 350 0 28 56926 1 160 229242 247 719 4 31 91721 1 161 351619 139 1277 4 40 115168 1 162 84207 29 356 11 30 111194 1 163 120445 118 457 0 16 51009 0 164 324598 110 1402 0 37 135777 1 165 131069 67 600 4 30 51513 1 166 204271 42 480 0 35 74163 1 167 165543 65 595 1 32 51633 1 168 141722 94 436 0 27 75345 1 169 116048 64 230 0 20 33416 0 170 250047 81 651 0 18 83305 0 171 299775 95 1367 9 31 98952 1 172 195838 67 564 1 31 102372 1 173 173260 63 716 3 21 37238 1 174 254488 83 747 10 39 103772 1 175 104389 45 467 5 41 123969 1 176 136084 30 671 0 13 27142 0 177 199476 70 861 2 32 135400 1 178 92499 32 319 0 18 21399 0 179 224330 83 612 1 39 130115 1 180 135781 31 433 2 14 24874 0 181 74408 67 434 4 7 34988 0 182 81240 66 503 0 17 45549 0 183 14688 10 85 0 0 6023 1 184 181633 70 564 2 30 64466 1 185 271856 103 824 1 37 54990 1 186 7199 5 74 0 0 1644 1 187 46660 20 259 0 5 6179 1 188 17547 5 69 0 1 3926 1 189 133368 36 535 1 16 32755 0 190 95227 34 239 0 32 34777 1 191 152601 48 438 2 24 73224 1 192 98146 40 459 0 17 27114 0 193 79619 43 426 3 11 20760 0 194 59194 31 288 6 24 37636 0 195 139942 42 498 0 22 65461 0 196 118612 46 454 2 12 30080 0 197 72880 33 376 0 19 24094 0 198 65475 18 225 2 13 69008 0 199 99643 55 555 1 17 54968 0 200 71965 35 252 1 15 46090 0 201 77272 59 208 2 16 27507 0 202 49289 19 130 1 24 10672 0 203 135131 66 481 0 15 34029 0 204 108446 60 389 1 17 46300 0 205 89746 36 565 3 18 24760 0 206 44296 25 173 0 20 18779 0 207 77648 47 278 0 16 21280 0 208 181528 54 609 0 16 40662 0 209 134019 53 422 0 18 28987 0 210 124064 40 445 1 22 22827 0 211 92630 40 387 4 8 18513 0 212 121848 39 339 0 17 30594 0 213 52915 14 181 0 18 24006 0 214 81872 45 245 0 16 27913 0 215 58981 36 384 7 23 42744 0 216 53515 28 212 2 22 12934 0 217 60812 44 399 0 13 22574 0 218 56375 30 229 7 13 41385 0 219 65490 22 224 3 16 18653 0 220 80949 17 203 0 16 18472 0 221 76302 31 333 0 20 30976 0 222 104011 55 384 6 22 63339 0 223 98104 54 636 2 17 25568 0 224 67989 21 185 0 18 33747 0 225 30989 14 93 0 17 4154 0 226 135458 81 581 3 12 19474 0 227 73504 35 248 0 7 35130 0 228 63123 43 304 1 17 39067 0 229 61254 46 344 1 14 13310 0 230 74914 30 407 0 23 65892 0 231 31774 23 170 1 17 4143 0 232 81437 38 312 0 14 28579 0 233 87186 54 507 0 15 51776 0 234 50090 20 224 0 17 21152 0 235 65745 53 340 0 21 38084 0 236 56653 45 168 0 18 27717 0 237 158399 39 443 0 18 32928 0 238 46455 20 204 0 17 11342 0 239 73624 24 367 0 17 19499 0 240 38395 31 210 0 16 16380 0 241 91899 35 335 0 15 36874 0 242 139526 151 364 0 21 48259 0 243 52164 52 178 0 16 16734 0 244 51567 30 206 2 14 28207 0 245 70551 31 279 0 15 30143 0 246 84856 29 387 1 17 41369 0 247 102538 57 490 1 15 45833 0 248 86678 40 238 0 15 29156 0 249 85709 44 343 0 10 35944 0 250 34662 25 232 0 6 36278 0 251 150580 77 530 0 22 45588 0 252 99611 35 291 0 21 45097 0 253 19349 11 67 0 1 3895 0 254 99373 63 397 1 18 28394 0 255 86230 44 467 0 17 18632 0 256 30837 19 178 0 4 2325 0 257 31706 13 175 0 10 25139 0 258 89806 42 299 0 16 27975 0 259 62088 38 154 1 16 14483 0 260 40151 29 106 0 9 13127 0 261 27634 20 189 0 16 5839 0 262 76990 27 194 0 17 24069 0 263 37460 20 135 0 7 3738 0 264 54157 19 201 0 15 18625 0 265 49862 37 207 0 14 36341 0 266 84337 26 280 0 14 24548 0 267 64175 42 260 0 18 21792 0 268 59382 49 227 0 12 26263 0 269 119308 30 239 0 16 23686 0 270 76702 49 333 0 21 49303 0 271 103425 67 428 1 19 25659 0 272 70344 28 230 0 16 28904 0 273 43410 19 292 0 1 2781 0 274 104838 49 350 1 16 29236 0 275 62215 27 186 0 10 19546 0 276 69304 30 326 6 19 22818 0 277 53117 22 155 3 12 32689 0 278 19764 12 75 1 2 5752 0 279 86680 31 361 2 14 22197 0 280 84105 20 261 0 17 20055 0 281 77945 20 299 0 19 25272 0 282 89113 39 300 0 14 82206 0 283 91005 29 450 3 11 32073 0 284 40248 16 183 1 4 5444 0 285 64187 27 238 0 16 20154 0 286 50857 21 165 0 20 36944 0 287 56613 19 234 1 12 8019 0 288 62792 35 176 0 15 30884 0 289 72535 14 329 0 16 19540 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Logins Views Shared Reviewed TotalSize -1.516e+04 2.831e+02 1.639e+02 -1.655e+03 1.499e+03 3.693e-01 Pop 6.459e+03 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -125849 -15485 625 15770 105643 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -1.516e+04 4.286e+03 -3.538 0.000472 *** Logins 2.831e+02 7.085e+01 3.995 8.26e-05 *** Views 1.639e+02 9.291e+00 17.639 < 2e-16 *** Shared -1.655e+03 6.833e+02 -2.422 0.016087 * Reviewed 1.499e+03 2.672e+02 5.610 4.82e-08 *** TotalSize 3.693e-01 6.740e-02 5.479 9.50e-08 *** Pop 6.459e+03 4.682e+03 1.380 0.168783 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 28650 on 282 degrees of freedom Multiple R-squared: 0.8814, Adjusted R-squared: 0.8789 F-statistic: 349.4 on 6 and 282 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,] 0.9987200 2.559993e-03 1.279996e-03 [2,] 0.9966509 6.698286e-03 3.349143e-03 [3,] 0.9926542 1.469167e-02 7.345835e-03 [4,] 0.9860749 2.785019e-02 1.392510e-02 [5,] 0.9744625 5.107507e-02 2.553753e-02 [6,] 0.9608103 7.837947e-02 3.918973e-02 [7,] 0.9660898 6.782035e-02 3.391017e-02 [8,] 0.9793712 4.125769e-02 2.062884e-02 [9,] 0.9760960 4.780807e-02 2.390403e-02 [10,] 0.9671504 6.569916e-02 3.284958e-02 [11,] 0.9525469 9.490622e-02 4.745311e-02 [12,] 0.9321930 1.356140e-01 6.780701e-02 [13,] 0.9058463 1.883073e-01 9.415366e-02 [14,] 0.9967806 6.438794e-03 3.219397e-03 [15,] 0.9957175 8.564955e-03 4.282477e-03 [16,] 0.9957841 8.431722e-03 4.215861e-03 [17,] 0.9995608 8.783659e-04 4.391829e-04 [18,] 0.9996493 7.014486e-04 3.507243e-04 [19,] 0.9994300 1.140077e-03 5.700387e-04 [20,] 0.9997594 4.812001e-04 2.406001e-04 [21,] 0.9996388 7.223503e-04 3.611752e-04 [22,] 0.9994646 1.070869e-03 5.354346e-04 [23,] 0.9991516 1.696886e-03 8.484430e-04 [24,] 0.9988258 2.348395e-03 1.174198e-03 [25,] 0.9991889 1.622177e-03 8.110887e-04 [26,] 0.9991647 1.670667e-03 8.353334e-04 [27,] 0.9990961 1.807867e-03 9.039336e-04 [28,] 0.9990364 1.927238e-03 9.636188e-04 [29,] 0.9986219 2.756294e-03 1.378147e-03 [30,] 0.9997004 5.992457e-04 2.996229e-04 [31,] 0.9995379 9.241288e-04 4.620644e-04 [32,] 0.9993460 1.307967e-03 6.539833e-04 [33,] 0.9992664 1.467168e-03 7.335838e-04 [34,] 0.9992002 1.599604e-03 7.998021e-04 [35,] 0.9990570 1.885931e-03 9.429654e-04 [36,] 0.9986275 2.744955e-03 1.372477e-03 [37,] 0.9981775 3.645002e-03 1.822501e-03 [38,] 0.9981722 3.655576e-03 1.827788e-03 [39,] 0.9974453 5.109464e-03 2.554732e-03 [40,] 0.9964984 7.003108e-03 3.501554e-03 [41,] 0.9999959 8.253026e-06 4.126513e-06 [42,] 0.9999938 1.242476e-05 6.212381e-06 [43,] 0.9999928 1.439800e-05 7.199002e-06 [44,] 0.9999890 2.199257e-05 1.099629e-05 [45,] 0.9999882 2.354674e-05 1.177337e-05 [46,] 0.9999929 1.418736e-05 7.093680e-06 [47,] 0.9999895 2.090773e-05 1.045386e-05 [48,] 0.9999861 2.778799e-05 1.389399e-05 [49,] 0.9999914 1.722657e-05 8.613283e-06 [50,] 0.9999865 2.702709e-05 1.351354e-05 [51,] 0.9999799 4.019701e-05 2.009851e-05 [52,] 0.9999700 6.004433e-05 3.002216e-05 [53,] 0.9999627 7.462685e-05 3.731343e-05 [54,] 0.9999479 1.041143e-04 5.205717e-05 [55,] 0.9999247 1.506990e-04 7.534948e-05 [56,] 0.9999066 1.868516e-04 9.342581e-05 [57,] 0.9998674 2.651745e-04 1.325873e-04 [58,] 0.9998947 2.105412e-04 1.052706e-04 [59,] 0.9998636 2.728702e-04 1.364351e-04 [60,] 0.9998356 3.288187e-04 1.644094e-04 [61,] 0.9998299 3.402187e-04 1.701093e-04 [62,] 0.9997699 4.601703e-04 2.300852e-04 [63,] 0.9997509 4.982729e-04 2.491365e-04 [64,] 0.9998402 3.196801e-04 1.598400e-04 [65,] 0.9998644 2.712980e-04 1.356490e-04 [66,] 0.9998387 3.225698e-04 1.612849e-04 [67,] 0.9997793 4.413119e-04 2.206559e-04 [68,] 0.9997357 5.286151e-04 2.643075e-04 [69,] 0.9996372 7.256995e-04 3.628497e-04 [70,] 0.9995558 8.884492e-04 4.442246e-04 [71,] 0.9994421 1.115735e-03 5.578674e-04 [72,] 0.9992543 1.491439e-03 7.457197e-04 [73,] 0.9990841 1.831808e-03 9.159040e-04 [74,] 0.9995757 8.486682e-04 4.243341e-04 [75,] 0.9996069 7.862637e-04 3.931319e-04 [76,] 0.9997600 4.800976e-04 2.400488e-04 [77,] 0.9996721 6.558333e-04 3.279167e-04 [78,] 0.9995804 8.391973e-04 4.195987e-04 [79,] 0.9994594 1.081143e-03 5.405713e-04 [80,] 0.9998708 2.583371e-04 1.291685e-04 [81,] 0.9998258 3.484854e-04 1.742427e-04 [82,] 0.9997987 4.025035e-04 2.012517e-04 [83,] 0.9997640 4.719564e-04 2.359782e-04 [84,] 0.9998155 3.690959e-04 1.845480e-04 [85,] 0.9997507 4.985642e-04 2.492821e-04 [86,] 0.9997374 5.251677e-04 2.625838e-04 [87,] 0.9996536 6.927747e-04 3.463873e-04 [88,] 0.9995644 8.712606e-04 4.356303e-04 [89,] 0.9999988 2.359031e-06 1.179516e-06 [90,] 0.9999986 2.752758e-06 1.376379e-06 [91,] 0.9999979 4.110388e-06 2.055194e-06 [92,] 0.9999971 5.853512e-06 2.926756e-06 [93,] 0.9999957 8.546289e-06 4.273144e-06 [94,] 0.9999955 9.072055e-06 4.536028e-06 [95,] 0.9999933 1.346397e-05 6.731987e-06 [96,] 0.9999911 1.787574e-05 8.937869e-06 [97,] 0.9999886 2.278848e-05 1.139424e-05 [98,] 0.9999848 3.042648e-05 1.521324e-05 [99,] 0.9999808 3.832182e-05 1.916091e-05 [100,] 0.9999797 4.064590e-05 2.032295e-05 [101,] 0.9999724 5.517753e-05 2.758877e-05 [102,] 0.9999633 7.347224e-05 3.673612e-05 [103,] 0.9999592 8.165430e-05 4.082715e-05 [104,] 0.9999542 9.166666e-05 4.583333e-05 [105,] 0.9999382 1.236688e-04 6.183441e-05 [106,] 0.9999400 1.200711e-04 6.003557e-05 [107,] 0.9999388 1.224715e-04 6.123577e-05 [108,] 0.9999152 1.696650e-04 8.483248e-05 [109,] 0.9999510 9.797229e-05 4.898615e-05 [110,] 0.9999809 3.814179e-05 1.907089e-05 [111,] 0.9999861 2.786399e-05 1.393200e-05 [112,] 0.9999801 3.989023e-05 1.994512e-05 [113,] 0.9999892 2.160975e-05 1.080488e-05 [114,] 0.9999889 2.224236e-05 1.112118e-05 [115,] 0.9999989 2.181898e-06 1.090949e-06 [116,] 0.9999998 4.565769e-07 2.282885e-07 [117,] 0.9999998 3.787828e-07 1.893914e-07 [118,] 0.9999997 5.511937e-07 2.755969e-07 [119,] 0.9999996 8.456919e-07 4.228460e-07 [120,] 0.9999996 7.054938e-07 3.527469e-07 [121,] 0.9999995 1.081331e-06 5.406656e-07 [122,] 0.9999992 1.640255e-06 8.201274e-07 [123,] 0.9999991 1.784986e-06 8.924931e-07 [124,] 0.9999988 2.397076e-06 1.198538e-06 [125,] 0.9999997 5.738172e-07 2.869086e-07 [126,] 0.9999996 8.610743e-07 4.305372e-07 [127,] 0.9999999 1.484519e-07 7.422594e-08 [128,] 0.9999999 2.359066e-07 1.179533e-07 [129,] 0.9999999 2.023076e-07 1.011538e-07 [130,] 0.9999999 1.677164e-07 8.385821e-08 [131,] 0.9999999 2.296261e-07 1.148131e-07 [132,] 0.9999999 2.846984e-07 1.423492e-07 [133,] 0.9999998 4.377504e-07 2.188752e-07 [134,] 0.9999997 6.800348e-07 3.400174e-07 [135,] 0.9999995 1.024610e-06 5.123049e-07 [136,] 0.9999993 1.406416e-06 7.032081e-07 [137,] 0.9999994 1.299024e-06 6.495120e-07 [138,] 0.9999991 1.797254e-06 8.986268e-07 [139,] 0.9999994 1.248047e-06 6.240236e-07 [140,] 0.9999994 1.259733e-06 6.298666e-07 [141,] 0.9999991 1.785923e-06 8.929615e-07 [142,] 0.9999991 1.789120e-06 8.945602e-07 [143,] 0.9999987 2.640374e-06 1.320187e-06 [144,] 0.9999980 3.988206e-06 1.994103e-06 [145,] 0.9999974 5.127955e-06 2.563978e-06 [146,] 0.9999975 4.937430e-06 2.468715e-06 [147,] 0.9999984 3.145984e-06 1.572992e-06 [148,] 0.9999983 3.368484e-06 1.684242e-06 [149,] 0.9999978 4.493998e-06 2.246999e-06 [150,] 0.9999982 3.537131e-06 1.768565e-06 [151,] 0.9999983 3.479915e-06 1.739958e-06 [152,] 0.9999979 4.225194e-06 2.112597e-06 [153,] 0.9999980 4.004088e-06 2.002044e-06 [154,] 0.9999977 4.646530e-06 2.323265e-06 [155,] 0.9999979 4.255864e-06 2.127932e-06 [156,] 0.9999985 3.005035e-06 1.502517e-06 [157,] 0.9999994 1.143740e-06 5.718702e-07 [158,] 0.9999991 1.722931e-06 8.614657e-07 [159,] 0.9999990 2.016586e-06 1.008293e-06 [160,] 0.9999992 1.697484e-06 8.487420e-07 [161,] 1.0000000 2.077579e-08 1.038790e-08 [162,] 1.0000000 3.381893e-08 1.690946e-08 [163,] 1.0000000 4.266477e-08 2.133238e-08 [164,] 1.0000000 6.958600e-08 3.479300e-08 [165,] 1.0000000 1.147686e-08 5.738430e-09 [166,] 1.0000000 1.319759e-09 6.598794e-10 [167,] 1.0000000 2.311449e-09 1.155724e-09 [168,] 1.0000000 1.057180e-09 5.285900e-10 [169,] 1.0000000 1.699891e-09 8.499454e-10 [170,] 1.0000000 2.760897e-09 1.380448e-09 [171,] 1.0000000 6.575299e-10 3.287649e-10 [172,] 1.0000000 8.719129e-10 4.359564e-10 [173,] 1.0000000 1.434277e-10 7.171387e-11 [174,] 1.0000000 2.254112e-10 1.127056e-10 [175,] 1.0000000 3.779802e-10 1.889901e-10 [176,] 1.0000000 7.167197e-11 3.583598e-11 [177,] 1.0000000 1.085619e-10 5.428096e-11 [178,] 1.0000000 1.386072e-10 6.930359e-11 [179,] 1.0000000 1.704848e-10 8.524239e-11 [180,] 1.0000000 2.011115e-10 1.005557e-10 [181,] 1.0000000 2.889156e-10 1.444578e-10 [182,] 1.0000000 5.234176e-10 2.617088e-10 [183,] 1.0000000 9.348223e-10 4.674111e-10 [184,] 1.0000000 1.624339e-09 8.121695e-10 [185,] 1.0000000 2.632402e-09 1.316201e-09 [186,] 1.0000000 3.477874e-09 1.738937e-09 [187,] 1.0000000 3.987585e-09 1.993793e-09 [188,] 1.0000000 5.321290e-09 2.660645e-09 [189,] 1.0000000 9.474837e-09 4.737418e-09 [190,] 1.0000000 7.070885e-09 3.535443e-09 [191,] 1.0000000 1.295911e-08 6.479553e-09 [192,] 1.0000000 2.161659e-08 1.080830e-08 [193,] 1.0000000 3.849009e-08 1.924505e-08 [194,] 1.0000000 5.095810e-08 2.547905e-08 [195,] 1.0000000 8.845449e-08 4.422724e-08 [196,] 1.0000000 8.736445e-08 4.368223e-08 [197,] 0.9999999 1.344234e-07 6.721169e-08 [198,] 0.9999999 2.372196e-07 1.186098e-07 [199,] 1.0000000 3.573318e-08 1.786659e-08 [200,] 1.0000000 2.144854e-08 1.072427e-08 [201,] 1.0000000 2.039669e-08 1.019835e-08 [202,] 1.0000000 2.444289e-08 1.222145e-08 [203,] 1.0000000 7.845522e-09 3.922761e-09 [204,] 1.0000000 1.504052e-08 7.520259e-09 [205,] 1.0000000 2.573093e-08 1.286547e-08 [206,] 1.0000000 1.807702e-08 9.038512e-09 [207,] 1.0000000 3.116808e-08 1.558404e-08 [208,] 1.0000000 2.308020e-08 1.154010e-08 [209,] 1.0000000 4.400686e-08 2.200343e-08 [210,] 1.0000000 7.511171e-08 3.755585e-08 [211,] 1.0000000 5.667909e-08 2.833954e-08 [212,] 0.9999999 1.019756e-07 5.098781e-08 [213,] 0.9999999 1.877732e-07 9.388661e-08 [214,] 0.9999999 1.089234e-07 5.446168e-08 [215,] 0.9999999 1.750331e-07 8.751657e-08 [216,] 0.9999998 3.264709e-07 1.632354e-07 [217,] 0.9999997 5.445549e-07 2.722775e-07 [218,] 0.9999996 8.436940e-07 4.218470e-07 [219,] 0.9999995 1.052287e-06 5.261435e-07 [220,] 0.9999994 1.171083e-06 5.855415e-07 [221,] 0.9999997 6.780617e-07 3.390308e-07 [222,] 0.9999996 8.456853e-07 4.228427e-07 [223,] 0.9999992 1.582200e-06 7.910998e-07 [224,] 0.9999996 7.410761e-07 3.705381e-07 [225,] 0.9999994 1.197626e-06 5.988128e-07 [226,] 0.9999998 4.735410e-07 2.367705e-07 [227,] 0.9999996 8.613557e-07 4.306779e-07 [228,] 1.0000000 1.214814e-08 6.074068e-09 [229,] 1.0000000 2.165686e-08 1.082843e-08 [230,] 1.0000000 4.038017e-08 2.019009e-08 [231,] 1.0000000 2.592709e-08 1.296355e-08 [232,] 1.0000000 5.159209e-08 2.579604e-08 [233,] 0.9999999 1.123592e-07 5.617959e-08 [234,] 0.9999999 1.608400e-07 8.042002e-08 [235,] 0.9999998 3.083562e-07 1.541781e-07 [236,] 0.9999997 6.602484e-07 3.301242e-07 [237,] 0.9999993 1.336971e-06 6.684857e-07 [238,] 0.9999988 2.356361e-06 1.178181e-06 [239,] 0.9999983 3.360170e-06 1.680085e-06 [240,] 0.9999966 6.835874e-06 3.417937e-06 [241,] 0.9999967 6.597967e-06 3.298983e-06 [242,] 0.9999958 8.346788e-06 4.173394e-06 [243,] 0.9999932 1.354031e-05 6.770153e-06 [244,] 0.9999860 2.798781e-05 1.399390e-05 [245,] 0.9999710 5.805939e-05 2.902969e-05 [246,] 0.9999647 7.066222e-05 3.533111e-05 [247,] 0.9999357 1.285493e-04 6.427464e-05 [248,] 0.9999208 1.584448e-04 7.922238e-05 [249,] 0.9998579 2.841229e-04 1.420614e-04 [250,] 0.9997462 5.075406e-04 2.537703e-04 [251,] 0.9995041 9.917396e-04 4.958698e-04 [252,] 0.9997655 4.690350e-04 2.345175e-04 [253,] 0.9996334 7.331755e-04 3.665878e-04 [254,] 0.9992609 1.478150e-03 7.390750e-04 [255,] 0.9986414 2.717262e-03 1.358631e-03 [256,] 0.9981870 3.626054e-03 1.813027e-03 [257,] 0.9969997 6.000669e-03 3.000335e-03 [258,] 0.9957145 8.571022e-03 4.285511e-03 [259,] 0.9931683 1.366334e-02 6.831671e-03 [260,] 0.9999505 9.903193e-05 4.951597e-05 [261,] 0.9999839 3.229061e-05 1.614531e-05 [262,] 0.9999861 2.778975e-05 1.389488e-05 [263,] 0.9999459 1.082352e-04 5.411761e-05 [264,] 0.9999552 8.967855e-05 4.483928e-05 [265,] 0.9998405 3.190098e-04 1.595049e-04 [266,] 0.9994868 1.026442e-03 5.132211e-04 [267,] 0.9985044 2.991207e-03 1.495603e-03 [268,] 0.9968047 6.390590e-03 3.195295e-03 [269,] 0.9876807 2.463865e-02 1.231932e-02 [270,] 0.9553589 8.928219e-02 4.464109e-02 > postscript(file="/var/wessaorg/rcomp/tmp/16j6x1355218833.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/wessaorg/rcomp/tmp/2y2ko1355218833.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/wessaorg/rcomp/tmp/3joe61355218833.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/wessaorg/rcomp/tmp/4g21p1355218833.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/wessaorg/rcomp/tmp/5son71355218833.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 = 289 Frequency = 1 1 2 3 4 5 57381.11794 -1511.90811 10501.17082 -74681.58651 30819.96881 6 7 8 9 10 -28378.61573 30717.59077 -16496.24686 15351.12634 -35806.27187 11 12 13 14 15 22297.04832 33374.39258 47.29319 -14512.12552 20840.28937 16 17 18 19 20 32680.77750 -51961.21467 36590.20525 -8944.82804 -1805.12382 21 22 23 24 25 16657.91439 12157.59221 105642.71844 34654.50653 -19764.26410 26 27 28 29 30 -77424.81227 -53496.76887 -13349.23944 51291.41395 -3494.46935 31 32 33 34 35 27496.39450 6215.55758 -9941.44218 49487.62786 -19571.80669 36 37 38 39 40 35111.90998 31413.06822 8677.49503 68999.86930 -930.54728 41 42 43 44 45 15321.23850 33195.82855 -32664.16493 5860.64248 -2899.41017 46 47 48 49 50 -10960.83946 -39649.27269 -5239.37136 447.46891 -125849.14401 51 52 53 54 55 -15485.14617 -26339.15190 -7173.40119 -28171.72815 -45348.42941 56 57 58 59 60 16089.14399 -17823.72420 51464.29413 -3605.32157 -10452.14157 61 62 63 64 65 2505.82322 28728.51718 20529.52422 6896.81673 -8827.80248 66 67 68 69 70 13342.09226 36990.65941 20776.23197 -9379.63009 -35209.79356 71 72 73 74 75 12671.81419 31339.97577 47419.76163 -24847.57652 -23377.73801 76 77 78 79 80 -7239.02624 -17973.81838 2517.94397 -22793.42358 15850.00152 81 82 83 84 85 -7584.09708 -15554.51361 55267.99624 35136.35832 -28151.02707 86 87 88 89 90 -7195.42067 -14146.61998 -14123.02164 69668.44704 -12731.14567 91 92 93 94 95 23330.06097 -21940.41106 -39669.57886 6843.37884 35675.30501 96 97 98 99 100 -10968.58808 17234.03830 -117585.07544 23942.12876 4888.69166 101 102 103 104 105 10039.32580 -3575.20141 26967.18332 -821.09625 13579.32886 106 107 108 109 110 16005.04785 -14917.78955 -6858.40724 24773.46998 7354.72167 111 112 113 114 115 12446.17238 -23982.64205 -25068.37781 10644.38922 -30074.45043 116 117 118 119 120 -25305.67103 8327.61408 55001.97162 -56316.34712 -29540.73584 121 122 123 124 125 4903.28961 -42085.63943 -25521.43149 -74625.98305 58475.78260 126 127 128 129 130 32053.34571 10635.28023 -2522.46536 32900.92429 1641.66583 131 132 133 134 135 -913.08066 -23974.12311 12852.10903 -54136.94656 -6791.24183 136 137 138 139 140 63386.79098 985.25829 27363.47235 -33587.04807 -15245.41479 141 142 143 144 145 -19072.70354 -6527.94620 1281.75086 -5893.30000 2075.56041 146 147 148 149 150 -29273.11738 -19430.81573 -40229.48323 -23173.55634 -13384.37261 151 152 153 154 155 28934.80243 6243.54112 -2809.87691 17253.88781 26862.87247 156 157 158 159 160 23649.04916 20176.49159 15695.89678 -35561.38334 -23535.62236 161 162 163 164 165 15814.68234 -41480.04988 -15514.94985 -33213.15976 -34907.34146 166 167 168 169 170 42557.73770 -7055.04988 -15941.52180 33073.81360 77843.37210 171 172 173 174 175 -10567.50642 10517.82044 6516.62322 37028.73077 -75156.95534 176 177 178 179 180 3274.91471 -47406.71555 11433.78388 4374.93135 44340.19551 181 182 183 184 185 -17316.73907 -47022.02775 4406.14997 12614.75872 42235.77627 186 187 188 189 190 1752.18093 -2522.53629 10577.55326 16232.29357 -5684.49964 191 192 193 194 195 16221.82643 -8738.12690 -6399.47350 -21571.45769 4443.50282 196 197 198 199 200 20560.04069 -20304.07951 -2995.87050 -35850.38990 -1933.77753 201 202 203 204 205 10807.51736 -503.27315 17726.98844 1942.31810 -29043.29067 206 207 208 209 210 -12892.32868 2098.58134 42595.32910 27327.55140 15214.97874 211 212 213 214 215 20834.74317 33627.72752 -1401.91244 9847.53970 -37664.50139 216 217 218 219 220 -8444.96432 -29697.52180 2325.32435 11801.18346 27219.72710 221 222 223 224 225 -13309.64206 -5773.03200 -37869.86553 7438.13087 -76.45526 226 227 228 229 230 12256.62081 14648.21807 -21967.69165 -17232.05350 -43932.95034 231 232 233 234 235 -12799.20274 3166.83219 -37635.26986 -10419.00216 -35365.45895 236 237 238 239 240 -5678.98516 50773.80086 -7153.96261 -10841.98093 -19671.96340 241 242 243 244 245 6146.39423 2984.16950 -6733.97382 -3620.33496 -2406.48171 246 247 248 249 250 -10723.89172 -16497.28395 18256.47603 3938.06610 -17665.47305 251 252 253 254 255 7268.01690 9036.43523 17479.59684 -4178.43422 -19965.43252 256 257 258 259 260 4593.71816 -9767.78007 9758.25289 13572.10719 11390.80740 261 262 263 264 265 -19985.32693 18338.82340 12960.62004 1632.24302 -13783.69448 266 267 268 269 270 16196.37801 -10197.03482 -4218.03514 54073.75254 -26271.69019 271 272 273 274 275 -6828.52716 5224.04983 2814.00833 15639.87683 17040.77017 276 277 278 279 280 -4437.28245 11550.69999 15770.09184 8027.17270 17937.45199 281 282 283 284 285 624.75388 -7275.78206 838.97164 14537.53539 1270.06348 286 287 288 289 -10595.53443 8749.38291 5308.42024 -1388.50917 > postscript(file="/var/wessaorg/rcomp/tmp/66d1e1355218833.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 = 289 Frequency = 1 lag(myerror, k = 1) myerror 0 57381.11794 NA 1 -1511.90811 57381.11794 2 10501.17082 -1511.90811 3 -74681.58651 10501.17082 4 30819.96881 -74681.58651 5 -28378.61573 30819.96881 6 30717.59077 -28378.61573 7 -16496.24686 30717.59077 8 15351.12634 -16496.24686 9 -35806.27187 15351.12634 10 22297.04832 -35806.27187 11 33374.39258 22297.04832 12 47.29319 33374.39258 13 -14512.12552 47.29319 14 20840.28937 -14512.12552 15 32680.77750 20840.28937 16 -51961.21467 32680.77750 17 36590.20525 -51961.21467 18 -8944.82804 36590.20525 19 -1805.12382 -8944.82804 20 16657.91439 -1805.12382 21 12157.59221 16657.91439 22 105642.71844 12157.59221 23 34654.50653 105642.71844 24 -19764.26410 34654.50653 25 -77424.81227 -19764.26410 26 -53496.76887 -77424.81227 27 -13349.23944 -53496.76887 28 51291.41395 -13349.23944 29 -3494.46935 51291.41395 30 27496.39450 -3494.46935 31 6215.55758 27496.39450 32 -9941.44218 6215.55758 33 49487.62786 -9941.44218 34 -19571.80669 49487.62786 35 35111.90998 -19571.80669 36 31413.06822 35111.90998 37 8677.49503 31413.06822 38 68999.86930 8677.49503 39 -930.54728 68999.86930 40 15321.23850 -930.54728 41 33195.82855 15321.23850 42 -32664.16493 33195.82855 43 5860.64248 -32664.16493 44 -2899.41017 5860.64248 45 -10960.83946 -2899.41017 46 -39649.27269 -10960.83946 47 -5239.37136 -39649.27269 48 447.46891 -5239.37136 49 -125849.14401 447.46891 50 -15485.14617 -125849.14401 51 -26339.15190 -15485.14617 52 -7173.40119 -26339.15190 53 -28171.72815 -7173.40119 54 -45348.42941 -28171.72815 55 16089.14399 -45348.42941 56 -17823.72420 16089.14399 57 51464.29413 -17823.72420 58 -3605.32157 51464.29413 59 -10452.14157 -3605.32157 60 2505.82322 -10452.14157 61 28728.51718 2505.82322 62 20529.52422 28728.51718 63 6896.81673 20529.52422 64 -8827.80248 6896.81673 65 13342.09226 -8827.80248 66 36990.65941 13342.09226 67 20776.23197 36990.65941 68 -9379.63009 20776.23197 69 -35209.79356 -9379.63009 70 12671.81419 -35209.79356 71 31339.97577 12671.81419 72 47419.76163 31339.97577 73 -24847.57652 47419.76163 74 -23377.73801 -24847.57652 75 -7239.02624 -23377.73801 76 -17973.81838 -7239.02624 77 2517.94397 -17973.81838 78 -22793.42358 2517.94397 79 15850.00152 -22793.42358 80 -7584.09708 15850.00152 81 -15554.51361 -7584.09708 82 55267.99624 -15554.51361 83 35136.35832 55267.99624 84 -28151.02707 35136.35832 85 -7195.42067 -28151.02707 86 -14146.61998 -7195.42067 87 -14123.02164 -14146.61998 88 69668.44704 -14123.02164 89 -12731.14567 69668.44704 90 23330.06097 -12731.14567 91 -21940.41106 23330.06097 92 -39669.57886 -21940.41106 93 6843.37884 -39669.57886 94 35675.30501 6843.37884 95 -10968.58808 35675.30501 96 17234.03830 -10968.58808 97 -117585.07544 17234.03830 98 23942.12876 -117585.07544 99 4888.69166 23942.12876 100 10039.32580 4888.69166 101 -3575.20141 10039.32580 102 26967.18332 -3575.20141 103 -821.09625 26967.18332 104 13579.32886 -821.09625 105 16005.04785 13579.32886 106 -14917.78955 16005.04785 107 -6858.40724 -14917.78955 108 24773.46998 -6858.40724 109 7354.72167 24773.46998 110 12446.17238 7354.72167 111 -23982.64205 12446.17238 112 -25068.37781 -23982.64205 113 10644.38922 -25068.37781 114 -30074.45043 10644.38922 115 -25305.67103 -30074.45043 116 8327.61408 -25305.67103 117 55001.97162 8327.61408 118 -56316.34712 55001.97162 119 -29540.73584 -56316.34712 120 4903.28961 -29540.73584 121 -42085.63943 4903.28961 122 -25521.43149 -42085.63943 123 -74625.98305 -25521.43149 124 58475.78260 -74625.98305 125 32053.34571 58475.78260 126 10635.28023 32053.34571 127 -2522.46536 10635.28023 128 32900.92429 -2522.46536 129 1641.66583 32900.92429 130 -913.08066 1641.66583 131 -23974.12311 -913.08066 132 12852.10903 -23974.12311 133 -54136.94656 12852.10903 134 -6791.24183 -54136.94656 135 63386.79098 -6791.24183 136 985.25829 63386.79098 137 27363.47235 985.25829 138 -33587.04807 27363.47235 139 -15245.41479 -33587.04807 140 -19072.70354 -15245.41479 141 -6527.94620 -19072.70354 142 1281.75086 -6527.94620 143 -5893.30000 1281.75086 144 2075.56041 -5893.30000 145 -29273.11738 2075.56041 146 -19430.81573 -29273.11738 147 -40229.48323 -19430.81573 148 -23173.55634 -40229.48323 149 -13384.37261 -23173.55634 150 28934.80243 -13384.37261 151 6243.54112 28934.80243 152 -2809.87691 6243.54112 153 17253.88781 -2809.87691 154 26862.87247 17253.88781 155 23649.04916 26862.87247 156 20176.49159 23649.04916 157 15695.89678 20176.49159 158 -35561.38334 15695.89678 159 -23535.62236 -35561.38334 160 15814.68234 -23535.62236 161 -41480.04988 15814.68234 162 -15514.94985 -41480.04988 163 -33213.15976 -15514.94985 164 -34907.34146 -33213.15976 165 42557.73770 -34907.34146 166 -7055.04988 42557.73770 167 -15941.52180 -7055.04988 168 33073.81360 -15941.52180 169 77843.37210 33073.81360 170 -10567.50642 77843.37210 171 10517.82044 -10567.50642 172 6516.62322 10517.82044 173 37028.73077 6516.62322 174 -75156.95534 37028.73077 175 3274.91471 -75156.95534 176 -47406.71555 3274.91471 177 11433.78388 -47406.71555 178 4374.93135 11433.78388 179 44340.19551 4374.93135 180 -17316.73907 44340.19551 181 -47022.02775 -17316.73907 182 4406.14997 -47022.02775 183 12614.75872 4406.14997 184 42235.77627 12614.75872 185 1752.18093 42235.77627 186 -2522.53629 1752.18093 187 10577.55326 -2522.53629 188 16232.29357 10577.55326 189 -5684.49964 16232.29357 190 16221.82643 -5684.49964 191 -8738.12690 16221.82643 192 -6399.47350 -8738.12690 193 -21571.45769 -6399.47350 194 4443.50282 -21571.45769 195 20560.04069 4443.50282 196 -20304.07951 20560.04069 197 -2995.87050 -20304.07951 198 -35850.38990 -2995.87050 199 -1933.77753 -35850.38990 200 10807.51736 -1933.77753 201 -503.27315 10807.51736 202 17726.98844 -503.27315 203 1942.31810 17726.98844 204 -29043.29067 1942.31810 205 -12892.32868 -29043.29067 206 2098.58134 -12892.32868 207 42595.32910 2098.58134 208 27327.55140 42595.32910 209 15214.97874 27327.55140 210 20834.74317 15214.97874 211 33627.72752 20834.74317 212 -1401.91244 33627.72752 213 9847.53970 -1401.91244 214 -37664.50139 9847.53970 215 -8444.96432 -37664.50139 216 -29697.52180 -8444.96432 217 2325.32435 -29697.52180 218 11801.18346 2325.32435 219 27219.72710 11801.18346 220 -13309.64206 27219.72710 221 -5773.03200 -13309.64206 222 -37869.86553 -5773.03200 223 7438.13087 -37869.86553 224 -76.45526 7438.13087 225 12256.62081 -76.45526 226 14648.21807 12256.62081 227 -21967.69165 14648.21807 228 -17232.05350 -21967.69165 229 -43932.95034 -17232.05350 230 -12799.20274 -43932.95034 231 3166.83219 -12799.20274 232 -37635.26986 3166.83219 233 -10419.00216 -37635.26986 234 -35365.45895 -10419.00216 235 -5678.98516 -35365.45895 236 50773.80086 -5678.98516 237 -7153.96261 50773.80086 238 -10841.98093 -7153.96261 239 -19671.96340 -10841.98093 240 6146.39423 -19671.96340 241 2984.16950 6146.39423 242 -6733.97382 2984.16950 243 -3620.33496 -6733.97382 244 -2406.48171 -3620.33496 245 -10723.89172 -2406.48171 246 -16497.28395 -10723.89172 247 18256.47603 -16497.28395 248 3938.06610 18256.47603 249 -17665.47305 3938.06610 250 7268.01690 -17665.47305 251 9036.43523 7268.01690 252 17479.59684 9036.43523 253 -4178.43422 17479.59684 254 -19965.43252 -4178.43422 255 4593.71816 -19965.43252 256 -9767.78007 4593.71816 257 9758.25289 -9767.78007 258 13572.10719 9758.25289 259 11390.80740 13572.10719 260 -19985.32693 11390.80740 261 18338.82340 -19985.32693 262 12960.62004 18338.82340 263 1632.24302 12960.62004 264 -13783.69448 1632.24302 265 16196.37801 -13783.69448 266 -10197.03482 16196.37801 267 -4218.03514 -10197.03482 268 54073.75254 -4218.03514 269 -26271.69019 54073.75254 270 -6828.52716 -26271.69019 271 5224.04983 -6828.52716 272 2814.00833 5224.04983 273 15639.87683 2814.00833 274 17040.77017 15639.87683 275 -4437.28245 17040.77017 276 11550.69999 -4437.28245 277 15770.09184 11550.69999 278 8027.17270 15770.09184 279 17937.45199 8027.17270 280 624.75388 17937.45199 281 -7275.78206 624.75388 282 838.97164 -7275.78206 283 14537.53539 838.97164 284 1270.06348 14537.53539 285 -10595.53443 1270.06348 286 8749.38291 -10595.53443 287 5308.42024 8749.38291 288 -1388.50917 5308.42024 289 NA -1388.50917 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -1511.90811 57381.11794 [2,] 10501.17082 -1511.90811 [3,] -74681.58651 10501.17082 [4,] 30819.96881 -74681.58651 [5,] -28378.61573 30819.96881 [6,] 30717.59077 -28378.61573 [7,] -16496.24686 30717.59077 [8,] 15351.12634 -16496.24686 [9,] -35806.27187 15351.12634 [10,] 22297.04832 -35806.27187 [11,] 33374.39258 22297.04832 [12,] 47.29319 33374.39258 [13,] -14512.12552 47.29319 [14,] 20840.28937 -14512.12552 [15,] 32680.77750 20840.28937 [16,] -51961.21467 32680.77750 [17,] 36590.20525 -51961.21467 [18,] -8944.82804 36590.20525 [19,] -1805.12382 -8944.82804 [20,] 16657.91439 -1805.12382 [21,] 12157.59221 16657.91439 [22,] 105642.71844 12157.59221 [23,] 34654.50653 105642.71844 [24,] -19764.26410 34654.50653 [25,] -77424.81227 -19764.26410 [26,] -53496.76887 -77424.81227 [27,] -13349.23944 -53496.76887 [28,] 51291.41395 -13349.23944 [29,] -3494.46935 51291.41395 [30,] 27496.39450 -3494.46935 [31,] 6215.55758 27496.39450 [32,] -9941.44218 6215.55758 [33,] 49487.62786 -9941.44218 [34,] -19571.80669 49487.62786 [35,] 35111.90998 -19571.80669 [36,] 31413.06822 35111.90998 [37,] 8677.49503 31413.06822 [38,] 68999.86930 8677.49503 [39,] -930.54728 68999.86930 [40,] 15321.23850 -930.54728 [41,] 33195.82855 15321.23850 [42,] -32664.16493 33195.82855 [43,] 5860.64248 -32664.16493 [44,] -2899.41017 5860.64248 [45,] -10960.83946 -2899.41017 [46,] -39649.27269 -10960.83946 [47,] -5239.37136 -39649.27269 [48,] 447.46891 -5239.37136 [49,] -125849.14401 447.46891 [50,] -15485.14617 -125849.14401 [51,] -26339.15190 -15485.14617 [52,] -7173.40119 -26339.15190 [53,] -28171.72815 -7173.40119 [54,] -45348.42941 -28171.72815 [55,] 16089.14399 -45348.42941 [56,] -17823.72420 16089.14399 [57,] 51464.29413 -17823.72420 [58,] -3605.32157 51464.29413 [59,] -10452.14157 -3605.32157 [60,] 2505.82322 -10452.14157 [61,] 28728.51718 2505.82322 [62,] 20529.52422 28728.51718 [63,] 6896.81673 20529.52422 [64,] -8827.80248 6896.81673 [65,] 13342.09226 -8827.80248 [66,] 36990.65941 13342.09226 [67,] 20776.23197 36990.65941 [68,] -9379.63009 20776.23197 [69,] -35209.79356 -9379.63009 [70,] 12671.81419 -35209.79356 [71,] 31339.97577 12671.81419 [72,] 47419.76163 31339.97577 [73,] -24847.57652 47419.76163 [74,] -23377.73801 -24847.57652 [75,] -7239.02624 -23377.73801 [76,] -17973.81838 -7239.02624 [77,] 2517.94397 -17973.81838 [78,] -22793.42358 2517.94397 [79,] 15850.00152 -22793.42358 [80,] -7584.09708 15850.00152 [81,] -15554.51361 -7584.09708 [82,] 55267.99624 -15554.51361 [83,] 35136.35832 55267.99624 [84,] -28151.02707 35136.35832 [85,] -7195.42067 -28151.02707 [86,] -14146.61998 -7195.42067 [87,] -14123.02164 -14146.61998 [88,] 69668.44704 -14123.02164 [89,] -12731.14567 69668.44704 [90,] 23330.06097 -12731.14567 [91,] -21940.41106 23330.06097 [92,] -39669.57886 -21940.41106 [93,] 6843.37884 -39669.57886 [94,] 35675.30501 6843.37884 [95,] -10968.58808 35675.30501 [96,] 17234.03830 -10968.58808 [97,] -117585.07544 17234.03830 [98,] 23942.12876 -117585.07544 [99,] 4888.69166 23942.12876 [100,] 10039.32580 4888.69166 [101,] -3575.20141 10039.32580 [102,] 26967.18332 -3575.20141 [103,] -821.09625 26967.18332 [104,] 13579.32886 -821.09625 [105,] 16005.04785 13579.32886 [106,] -14917.78955 16005.04785 [107,] -6858.40724 -14917.78955 [108,] 24773.46998 -6858.40724 [109,] 7354.72167 24773.46998 [110,] 12446.17238 7354.72167 [111,] -23982.64205 12446.17238 [112,] -25068.37781 -23982.64205 [113,] 10644.38922 -25068.37781 [114,] -30074.45043 10644.38922 [115,] -25305.67103 -30074.45043 [116,] 8327.61408 -25305.67103 [117,] 55001.97162 8327.61408 [118,] -56316.34712 55001.97162 [119,] -29540.73584 -56316.34712 [120,] 4903.28961 -29540.73584 [121,] -42085.63943 4903.28961 [122,] -25521.43149 -42085.63943 [123,] -74625.98305 -25521.43149 [124,] 58475.78260 -74625.98305 [125,] 32053.34571 58475.78260 [126,] 10635.28023 32053.34571 [127,] -2522.46536 10635.28023 [128,] 32900.92429 -2522.46536 [129,] 1641.66583 32900.92429 [130,] -913.08066 1641.66583 [131,] -23974.12311 -913.08066 [132,] 12852.10903 -23974.12311 [133,] -54136.94656 12852.10903 [134,] -6791.24183 -54136.94656 [135,] 63386.79098 -6791.24183 [136,] 985.25829 63386.79098 [137,] 27363.47235 985.25829 [138,] -33587.04807 27363.47235 [139,] -15245.41479 -33587.04807 [140,] -19072.70354 -15245.41479 [141,] -6527.94620 -19072.70354 [142,] 1281.75086 -6527.94620 [143,] -5893.30000 1281.75086 [144,] 2075.56041 -5893.30000 [145,] -29273.11738 2075.56041 [146,] -19430.81573 -29273.11738 [147,] -40229.48323 -19430.81573 [148,] -23173.55634 -40229.48323 [149,] -13384.37261 -23173.55634 [150,] 28934.80243 -13384.37261 [151,] 6243.54112 28934.80243 [152,] -2809.87691 6243.54112 [153,] 17253.88781 -2809.87691 [154,] 26862.87247 17253.88781 [155,] 23649.04916 26862.87247 [156,] 20176.49159 23649.04916 [157,] 15695.89678 20176.49159 [158,] -35561.38334 15695.89678 [159,] -23535.62236 -35561.38334 [160,] 15814.68234 -23535.62236 [161,] -41480.04988 15814.68234 [162,] -15514.94985 -41480.04988 [163,] -33213.15976 -15514.94985 [164,] -34907.34146 -33213.15976 [165,] 42557.73770 -34907.34146 [166,] -7055.04988 42557.73770 [167,] -15941.52180 -7055.04988 [168,] 33073.81360 -15941.52180 [169,] 77843.37210 33073.81360 [170,] -10567.50642 77843.37210 [171,] 10517.82044 -10567.50642 [172,] 6516.62322 10517.82044 [173,] 37028.73077 6516.62322 [174,] -75156.95534 37028.73077 [175,] 3274.91471 -75156.95534 [176,] -47406.71555 3274.91471 [177,] 11433.78388 -47406.71555 [178,] 4374.93135 11433.78388 [179,] 44340.19551 4374.93135 [180,] -17316.73907 44340.19551 [181,] -47022.02775 -17316.73907 [182,] 4406.14997 -47022.02775 [183,] 12614.75872 4406.14997 [184,] 42235.77627 12614.75872 [185,] 1752.18093 42235.77627 [186,] -2522.53629 1752.18093 [187,] 10577.55326 -2522.53629 [188,] 16232.29357 10577.55326 [189,] -5684.49964 16232.29357 [190,] 16221.82643 -5684.49964 [191,] -8738.12690 16221.82643 [192,] -6399.47350 -8738.12690 [193,] -21571.45769 -6399.47350 [194,] 4443.50282 -21571.45769 [195,] 20560.04069 4443.50282 [196,] -20304.07951 20560.04069 [197,] -2995.87050 -20304.07951 [198,] -35850.38990 -2995.87050 [199,] -1933.77753 -35850.38990 [200,] 10807.51736 -1933.77753 [201,] -503.27315 10807.51736 [202,] 17726.98844 -503.27315 [203,] 1942.31810 17726.98844 [204,] -29043.29067 1942.31810 [205,] -12892.32868 -29043.29067 [206,] 2098.58134 -12892.32868 [207,] 42595.32910 2098.58134 [208,] 27327.55140 42595.32910 [209,] 15214.97874 27327.55140 [210,] 20834.74317 15214.97874 [211,] 33627.72752 20834.74317 [212,] -1401.91244 33627.72752 [213,] 9847.53970 -1401.91244 [214,] -37664.50139 9847.53970 [215,] -8444.96432 -37664.50139 [216,] -29697.52180 -8444.96432 [217,] 2325.32435 -29697.52180 [218,] 11801.18346 2325.32435 [219,] 27219.72710 11801.18346 [220,] -13309.64206 27219.72710 [221,] -5773.03200 -13309.64206 [222,] -37869.86553 -5773.03200 [223,] 7438.13087 -37869.86553 [224,] -76.45526 7438.13087 [225,] 12256.62081 -76.45526 [226,] 14648.21807 12256.62081 [227,] -21967.69165 14648.21807 [228,] -17232.05350 -21967.69165 [229,] -43932.95034 -17232.05350 [230,] -12799.20274 -43932.95034 [231,] 3166.83219 -12799.20274 [232,] -37635.26986 3166.83219 [233,] -10419.00216 -37635.26986 [234,] -35365.45895 -10419.00216 [235,] -5678.98516 -35365.45895 [236,] 50773.80086 -5678.98516 [237,] -7153.96261 50773.80086 [238,] -10841.98093 -7153.96261 [239,] -19671.96340 -10841.98093 [240,] 6146.39423 -19671.96340 [241,] 2984.16950 6146.39423 [242,] -6733.97382 2984.16950 [243,] -3620.33496 -6733.97382 [244,] -2406.48171 -3620.33496 [245,] -10723.89172 -2406.48171 [246,] -16497.28395 -10723.89172 [247,] 18256.47603 -16497.28395 [248,] 3938.06610 18256.47603 [249,] -17665.47305 3938.06610 [250,] 7268.01690 -17665.47305 [251,] 9036.43523 7268.01690 [252,] 17479.59684 9036.43523 [253,] -4178.43422 17479.59684 [254,] -19965.43252 -4178.43422 [255,] 4593.71816 -19965.43252 [256,] -9767.78007 4593.71816 [257,] 9758.25289 -9767.78007 [258,] 13572.10719 9758.25289 [259,] 11390.80740 13572.10719 [260,] -19985.32693 11390.80740 [261,] 18338.82340 -19985.32693 [262,] 12960.62004 18338.82340 [263,] 1632.24302 12960.62004 [264,] -13783.69448 1632.24302 [265,] 16196.37801 -13783.69448 [266,] -10197.03482 16196.37801 [267,] -4218.03514 -10197.03482 [268,] 54073.75254 -4218.03514 [269,] -26271.69019 54073.75254 [270,] -6828.52716 -26271.69019 [271,] 5224.04983 -6828.52716 [272,] 2814.00833 5224.04983 [273,] 15639.87683 2814.00833 [274,] 17040.77017 15639.87683 [275,] -4437.28245 17040.77017 [276,] 11550.69999 -4437.28245 [277,] 15770.09184 11550.69999 [278,] 8027.17270 15770.09184 [279,] 17937.45199 8027.17270 [280,] 624.75388 17937.45199 [281,] -7275.78206 624.75388 [282,] 838.97164 -7275.78206 [283,] 14537.53539 838.97164 [284,] 1270.06348 14537.53539 [285,] -10595.53443 1270.06348 [286,] 8749.38291 -10595.53443 [287,] 5308.42024 8749.38291 [288,] -1388.50917 5308.42024 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -1511.90811 57381.11794 2 10501.17082 -1511.90811 3 -74681.58651 10501.17082 4 30819.96881 -74681.58651 5 -28378.61573 30819.96881 6 30717.59077 -28378.61573 7 -16496.24686 30717.59077 8 15351.12634 -16496.24686 9 -35806.27187 15351.12634 10 22297.04832 -35806.27187 11 33374.39258 22297.04832 12 47.29319 33374.39258 13 -14512.12552 47.29319 14 20840.28937 -14512.12552 15 32680.77750 20840.28937 16 -51961.21467 32680.77750 17 36590.20525 -51961.21467 18 -8944.82804 36590.20525 19 -1805.12382 -8944.82804 20 16657.91439 -1805.12382 21 12157.59221 16657.91439 22 105642.71844 12157.59221 23 34654.50653 105642.71844 24 -19764.26410 34654.50653 25 -77424.81227 -19764.26410 26 -53496.76887 -77424.81227 27 -13349.23944 -53496.76887 28 51291.41395 -13349.23944 29 -3494.46935 51291.41395 30 27496.39450 -3494.46935 31 6215.55758 27496.39450 32 -9941.44218 6215.55758 33 49487.62786 -9941.44218 34 -19571.80669 49487.62786 35 35111.90998 -19571.80669 36 31413.06822 35111.90998 37 8677.49503 31413.06822 38 68999.86930 8677.49503 39 -930.54728 68999.86930 40 15321.23850 -930.54728 41 33195.82855 15321.23850 42 -32664.16493 33195.82855 43 5860.64248 -32664.16493 44 -2899.41017 5860.64248 45 -10960.83946 -2899.41017 46 -39649.27269 -10960.83946 47 -5239.37136 -39649.27269 48 447.46891 -5239.37136 49 -125849.14401 447.46891 50 -15485.14617 -125849.14401 51 -26339.15190 -15485.14617 52 -7173.40119 -26339.15190 53 -28171.72815 -7173.40119 54 -45348.42941 -28171.72815 55 16089.14399 -45348.42941 56 -17823.72420 16089.14399 57 51464.29413 -17823.72420 58 -3605.32157 51464.29413 59 -10452.14157 -3605.32157 60 2505.82322 -10452.14157 61 28728.51718 2505.82322 62 20529.52422 28728.51718 63 6896.81673 20529.52422 64 -8827.80248 6896.81673 65 13342.09226 -8827.80248 66 36990.65941 13342.09226 67 20776.23197 36990.65941 68 -9379.63009 20776.23197 69 -35209.79356 -9379.63009 70 12671.81419 -35209.79356 71 31339.97577 12671.81419 72 47419.76163 31339.97577 73 -24847.57652 47419.76163 74 -23377.73801 -24847.57652 75 -7239.02624 -23377.73801 76 -17973.81838 -7239.02624 77 2517.94397 -17973.81838 78 -22793.42358 2517.94397 79 15850.00152 -22793.42358 80 -7584.09708 15850.00152 81 -15554.51361 -7584.09708 82 55267.99624 -15554.51361 83 35136.35832 55267.99624 84 -28151.02707 35136.35832 85 -7195.42067 -28151.02707 86 -14146.61998 -7195.42067 87 -14123.02164 -14146.61998 88 69668.44704 -14123.02164 89 -12731.14567 69668.44704 90 23330.06097 -12731.14567 91 -21940.41106 23330.06097 92 -39669.57886 -21940.41106 93 6843.37884 -39669.57886 94 35675.30501 6843.37884 95 -10968.58808 35675.30501 96 17234.03830 -10968.58808 97 -117585.07544 17234.03830 98 23942.12876 -117585.07544 99 4888.69166 23942.12876 100 10039.32580 4888.69166 101 -3575.20141 10039.32580 102 26967.18332 -3575.20141 103 -821.09625 26967.18332 104 13579.32886 -821.09625 105 16005.04785 13579.32886 106 -14917.78955 16005.04785 107 -6858.40724 -14917.78955 108 24773.46998 -6858.40724 109 7354.72167 24773.46998 110 12446.17238 7354.72167 111 -23982.64205 12446.17238 112 -25068.37781 -23982.64205 113 10644.38922 -25068.37781 114 -30074.45043 10644.38922 115 -25305.67103 -30074.45043 116 8327.61408 -25305.67103 117 55001.97162 8327.61408 118 -56316.34712 55001.97162 119 -29540.73584 -56316.34712 120 4903.28961 -29540.73584 121 -42085.63943 4903.28961 122 -25521.43149 -42085.63943 123 -74625.98305 -25521.43149 124 58475.78260 -74625.98305 125 32053.34571 58475.78260 126 10635.28023 32053.34571 127 -2522.46536 10635.28023 128 32900.92429 -2522.46536 129 1641.66583 32900.92429 130 -913.08066 1641.66583 131 -23974.12311 -913.08066 132 12852.10903 -23974.12311 133 -54136.94656 12852.10903 134 -6791.24183 -54136.94656 135 63386.79098 -6791.24183 136 985.25829 63386.79098 137 27363.47235 985.25829 138 -33587.04807 27363.47235 139 -15245.41479 -33587.04807 140 -19072.70354 -15245.41479 141 -6527.94620 -19072.70354 142 1281.75086 -6527.94620 143 -5893.30000 1281.75086 144 2075.56041 -5893.30000 145 -29273.11738 2075.56041 146 -19430.81573 -29273.11738 147 -40229.48323 -19430.81573 148 -23173.55634 -40229.48323 149 -13384.37261 -23173.55634 150 28934.80243 -13384.37261 151 6243.54112 28934.80243 152 -2809.87691 6243.54112 153 17253.88781 -2809.87691 154 26862.87247 17253.88781 155 23649.04916 26862.87247 156 20176.49159 23649.04916 157 15695.89678 20176.49159 158 -35561.38334 15695.89678 159 -23535.62236 -35561.38334 160 15814.68234 -23535.62236 161 -41480.04988 15814.68234 162 -15514.94985 -41480.04988 163 -33213.15976 -15514.94985 164 -34907.34146 -33213.15976 165 42557.73770 -34907.34146 166 -7055.04988 42557.73770 167 -15941.52180 -7055.04988 168 33073.81360 -15941.52180 169 77843.37210 33073.81360 170 -10567.50642 77843.37210 171 10517.82044 -10567.50642 172 6516.62322 10517.82044 173 37028.73077 6516.62322 174 -75156.95534 37028.73077 175 3274.91471 -75156.95534 176 -47406.71555 3274.91471 177 11433.78388 -47406.71555 178 4374.93135 11433.78388 179 44340.19551 4374.93135 180 -17316.73907 44340.19551 181 -47022.02775 -17316.73907 182 4406.14997 -47022.02775 183 12614.75872 4406.14997 184 42235.77627 12614.75872 185 1752.18093 42235.77627 186 -2522.53629 1752.18093 187 10577.55326 -2522.53629 188 16232.29357 10577.55326 189 -5684.49964 16232.29357 190 16221.82643 -5684.49964 191 -8738.12690 16221.82643 192 -6399.47350 -8738.12690 193 -21571.45769 -6399.47350 194 4443.50282 -21571.45769 195 20560.04069 4443.50282 196 -20304.07951 20560.04069 197 -2995.87050 -20304.07951 198 -35850.38990 -2995.87050 199 -1933.77753 -35850.38990 200 10807.51736 -1933.77753 201 -503.27315 10807.51736 202 17726.98844 -503.27315 203 1942.31810 17726.98844 204 -29043.29067 1942.31810 205 -12892.32868 -29043.29067 206 2098.58134 -12892.32868 207 42595.32910 2098.58134 208 27327.55140 42595.32910 209 15214.97874 27327.55140 210 20834.74317 15214.97874 211 33627.72752 20834.74317 212 -1401.91244 33627.72752 213 9847.53970 -1401.91244 214 -37664.50139 9847.53970 215 -8444.96432 -37664.50139 216 -29697.52180 -8444.96432 217 2325.32435 -29697.52180 218 11801.18346 2325.32435 219 27219.72710 11801.18346 220 -13309.64206 27219.72710 221 -5773.03200 -13309.64206 222 -37869.86553 -5773.03200 223 7438.13087 -37869.86553 224 -76.45526 7438.13087 225 12256.62081 -76.45526 226 14648.21807 12256.62081 227 -21967.69165 14648.21807 228 -17232.05350 -21967.69165 229 -43932.95034 -17232.05350 230 -12799.20274 -43932.95034 231 3166.83219 -12799.20274 232 -37635.26986 3166.83219 233 -10419.00216 -37635.26986 234 -35365.45895 -10419.00216 235 -5678.98516 -35365.45895 236 50773.80086 -5678.98516 237 -7153.96261 50773.80086 238 -10841.98093 -7153.96261 239 -19671.96340 -10841.98093 240 6146.39423 -19671.96340 241 2984.16950 6146.39423 242 -6733.97382 2984.16950 243 -3620.33496 -6733.97382 244 -2406.48171 -3620.33496 245 -10723.89172 -2406.48171 246 -16497.28395 -10723.89172 247 18256.47603 -16497.28395 248 3938.06610 18256.47603 249 -17665.47305 3938.06610 250 7268.01690 -17665.47305 251 9036.43523 7268.01690 252 17479.59684 9036.43523 253 -4178.43422 17479.59684 254 -19965.43252 -4178.43422 255 4593.71816 -19965.43252 256 -9767.78007 4593.71816 257 9758.25289 -9767.78007 258 13572.10719 9758.25289 259 11390.80740 13572.10719 260 -19985.32693 11390.80740 261 18338.82340 -19985.32693 262 12960.62004 18338.82340 263 1632.24302 12960.62004 264 -13783.69448 1632.24302 265 16196.37801 -13783.69448 266 -10197.03482 16196.37801 267 -4218.03514 -10197.03482 268 54073.75254 -4218.03514 269 -26271.69019 54073.75254 270 -6828.52716 -26271.69019 271 5224.04983 -6828.52716 272 2814.00833 5224.04983 273 15639.87683 2814.00833 274 17040.77017 15639.87683 275 -4437.28245 17040.77017 276 11550.69999 -4437.28245 277 15770.09184 11550.69999 278 8027.17270 15770.09184 279 17937.45199 8027.17270 280 624.75388 17937.45199 281 -7275.78206 624.75388 282 838.97164 -7275.78206 283 14537.53539 838.97164 284 1270.06348 14537.53539 285 -10595.53443 1270.06348 286 8749.38291 -10595.53443 287 5308.42024 8749.38291 288 -1388.50917 5308.42024 > 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/wessaorg/rcomp/tmp/7hjkl1355218833.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/wessaorg/rcomp/tmp/82jju1355218833.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/wessaorg/rcomp/tmp/9fx901355218833.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/wessaorg/rcomp/tmp/10e5q21355218833.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/111hde1355218833.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/wessaorg/rcomp/tmp/12d6k11355218833.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/wessaorg/rcomp/tmp/13dztz1355218833.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/wessaorg/rcomp/tmp/14ycou1355218833.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/wessaorg/rcomp/tmp/15stqa1355218833.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/wessaorg/rcomp/tmp/160cpb1355218833.tab") + } > > try(system("convert tmp/16j6x1355218833.ps tmp/16j6x1355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/2y2ko1355218833.ps tmp/2y2ko1355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/3joe61355218833.ps tmp/3joe61355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/4g21p1355218833.ps tmp/4g21p1355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/5son71355218833.ps tmp/5son71355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/66d1e1355218833.ps tmp/66d1e1355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/7hjkl1355218833.ps tmp/7hjkl1355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/82jju1355218833.ps tmp/82jju1355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/9fx901355218833.ps tmp/9fx901355218833.png",intern=TRUE)) character(0) > try(system("convert tmp/10e5q21355218833.ps tmp/10e5q21355218833.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.129 1.078 13.220