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(1418 + ,210907 + ,396 + ,81 + ,79 + ,30 + ,115 + ,94 + ,112285 + ,146283 + ,144 + ,869 + ,120982 + ,297 + ,55 + ,58 + ,28 + ,109 + ,103 + ,84786 + ,98364 + ,103 + ,1530 + ,176508 + ,559 + ,50 + ,60 + ,38 + ,146 + ,93 + ,83123 + ,86146 + ,98 + ,2172 + ,179321 + ,967 + ,125 + ,108 + ,30 + ,116 + ,103 + ,101193 + ,96933 + ,135 + ,901 + ,123185 + ,270 + ,40 + ,49 + ,22 + ,68 + ,51 + ,38361 + ,79234 + ,61 + ,463 + ,52746 + ,143 + ,37 + ,0 + ,26 + ,101 + ,70 + ,68504 + ,42551 + ,39 + ,3201 + ,385534 + ,1562 + ,63 + ,121 + ,25 + ,96 + ,91 + ,119182 + ,195663 + ,150 + ,371 + ,33170 + ,109 + ,44 + ,1 + ,18 + ,67 + ,22 + ,22807 + ,6853 + ,5 + ,1192 + ,101645 + ,371 + ,88 + ,20 + ,11 + ,44 + ,38 + ,17140 + ,21529 + ,28 + ,1583 + ,149061 + ,656 + ,66 + ,43 + ,26 + ,100 + ,93 + ,116174 + ,95757 + ,84 + ,1439 + ,165446 + ,511 + ,57 + ,69 + ,25 + ,93 + ,60 + ,57635 + ,85584 + ,80 + ,1764 + ,237213 + ,655 + ,74 + ,78 + ,38 + ,140 + ,123 + ,66198 + ,143983 + ,130 + ,1495 + ,173326 + ,465 + ,49 + ,86 + ,44 + ,166 + ,148 + ,71701 + ,75851 + ,82 + ,1373 + ,133131 + ,525 + ,52 + ,44 + ,30 + ,99 + ,90 + ,57793 + ,59238 + ,60 + ,2187 + ,258873 + ,885 + ,88 + ,104 + ,40 + ,139 + ,124 + ,80444 + ,93163 + ,131 + ,1491 + ,180083 + ,497 + ,36 + ,63 + ,34 + ,130 + ,70 + ,53855 + ,96037 + ,84 + ,4041 + ,324799 + ,1436 + ,108 + ,158 + ,47 + ,181 + ,168 + ,97668 + ,151511 + ,140 + ,1706 + ,230964 + ,612 + ,43 + ,102 + ,30 + ,116 + ,115 + ,133824 + ,136368 + ,151 + ,2152 + ,236785 + ,865 + ,75 + ,77 + ,31 + ,116 + ,71 + ,101481 + ,112642 + ,91 + ,1036 + ,135473 + ,385 + ,32 + ,82 + ,23 + ,88 + ,66 + ,99645 + ,94728 + ,138 + ,1882 + ,202925 + ,567 + ,44 + ,115 + ,36 + ,139 + ,134 + ,114789 + ,105499 + ,150 + ,1929 + ,215147 + ,639 + ,85 + ,101 + ,36 + ,135 + ,117 + ,99052 + ,121527 + ,124 + ,2242 + ,344297 + ,963 + ,86 + ,80 + ,30 + ,108 + ,108 + ,67654 + ,127766 + ,119 + ,1220 + ,153935 + ,398 + ,56 + ,50 + ,25 + ,89 + ,84 + ,65553 + ,98958 + ,73 + ,1289 + ,132943 + ,410 + ,50 + ,83 + ,39 + ,156 + ,156 + ,97500 + ,77900 + ,110 + ,2515 + ,174724 + ,966 + ,135 + ,123 + ,34 + ,129 + ,120 + ,69112 + ,85646 + ,123 + ,2147 + ,174415 + ,801 + ,63 + ,73 + ,31 + ,118 + ,114 + ,82753 + ,98579 + ,90 + ,2352 + ,225548 + ,892 + ,81 + ,81 + ,31 + ,118 + ,94 + ,85323 + ,130767 + ,116 + ,1638 + ,223632 + ,513 + ,52 + ,105 + ,33 + ,125 + ,120 + ,72654 + ,131741 + ,113 + ,1222 + ,124817 + ,469 + ,44 + ,47 + ,25 + ,95 + ,81 + ,30727 + ,53907 + ,56 + ,1812 + ,221698 + ,683 + ,113 + ,105 + ,33 + ,126 + ,110 + ,77873 + ,178812 + ,115 + ,1677 + ,210767 + ,643 + ,39 + ,94 + ,35 + ,135 + ,133 + ,117478 + ,146761 + ,119 + ,1579 + ,170266 + ,535 + ,73 + ,44 + ,42 + ,154 + ,122 + ,74007 + ,82036 + ,129 + ,1731 + ,260561 + ,625 + ,48 + ,114 + ,43 + ,165 + ,158 + ,90183 + ,163253 + ,127 + ,807 + ,84853 + ,264 + ,33 + ,38 + ,30 + ,113 + ,109 + ,61542 + ,27032 + ,27 + ,2452 + ,294424 + ,992 + ,59 + ,107 + ,33 + ,127 + ,124 + ,101494 + ,171975 + ,175 + ,829 + ,101011 + ,238 + ,41 + ,30 + ,13 + ,52 + ,39 + ,27570 + ,65990 + ,35 + ,1940 + ,215641 + ,818 + ,69 + ,71 + ,32 + ,121 + ,92 + ,55813 + ,86572 + ,64 + ,2662 + ,325107 + ,937 + ,64 + ,84 + ,36 + ,136 + ,126 + ,79215 + ,159676 + ,96 + ,186 + ,7176 + ,70 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1423 + ,1929 + ,0 + ,1499 + ,167542 + ,507 + ,59 + ,59 + ,28 + ,108 + ,70 + ,55461 + ,85371 + ,84 + ,865 + ,106408 + ,260 + ,32 + ,33 + ,14 + ,46 + ,37 + ,31081 + ,58391 + ,41 + ,1793 + ,96560 + ,503 + ,129 + ,42 + ,17 + ,54 + ,38 + ,22996 + ,31580 + ,47 + ,2527 + ,265769 + ,927 + ,37 + ,96 + ,32 + ,124 + ,120 + ,83122 + ,136815 + ,126 + ,2747 + ,269651 + ,1269 + ,31 + ,106 + ,30 + ,115 + ,93 + ,70106 + ,120642 + ,105 + ,1324 + ,149112 + ,537 + ,65 + ,56 + ,35 + ,128 + ,95 + ,60578 + ,69107 + ,80 + ,2702 + ,175824 + ,910 + ,107 + ,57 + ,20 + ,80 + ,77 + ,39992 + ,50495 + ,70 + ,1383 + ,152871 + ,532 + ,74 + ,59 + ,28 + ,97 + ,90 + ,79892 + ,108016 + ,73 + ,1179 + ,111665 + ,345 + ,54 + ,39 + ,28 + ,104 + ,80 + ,49810 + ,46341 + ,57 + ,2099 + ,116408 + ,918 + ,76 + ,34 + ,39 + ,59 + ,31 + ,71570 + ,78348 + ,40 + ,4308 + ,362301 + ,1635 + ,715 + ,76 + ,34 + ,125 + ,110 + ,100708 + ,79336 + ,68 + ,918 + ,78800 + ,330 + ,57 + ,20 + ,26 + ,82 + ,66 + ,33032 + ,56968 + ,21 + ,1831 + ,183167 + ,557 + ,66 + ,91 + ,39 + ,149 + ,138 + ,82875 + ,93176 + ,127 + ,3373 + ,277965 + ,1178 + ,106 + ,115 + ,39 + ,149 + ,133 + ,139077 + ,161632 + ,154 + ,1713 + ,150629 + ,740 + ,54 + ,85 + ,33 + ,122 + ,113 + ,71595 + ,87850 + ,116 + ,1438 + ,168809 + ,452 + ,32 + ,76 + ,28 + ,118 + ,100 + ,72260 + ,127969 + ,102 + ,496 + ,24188 + ,218 + ,20 + ,8 + ,4 + ,12 + ,7 + ,5950 + ,15049 + ,7 + ,2253 + ,329267 + ,764 + ,71 + ,79 + ,39 + ,144 + ,140 + ,115762 + ,155135 + ,148 + ,744 + ,65029 + ,255 + ,21 + ,21 + ,18 + ,67 + ,61 + ,32551 + ,25109 + ,21 + ,1161 + ,101097 + ,454 + ,70 + ,30 + ,14 + ,52 + ,41 + ,31701 + ,45824 + ,35 + ,2352 + ,218946 + ,866 + ,112 + ,76 + ,29 + ,108 + ,96 + ,80670 + ,102996 + ,112 + ,2144 + ,244052 + ,574 + ,66 + ,101 + ,44 + ,166 + ,164 + ,143558 + ,160604 + ,137 + ,4691 + ,341570 + ,1276 + ,190 + ,94 + ,21 + ,80 + ,78 + ,117105 + ,158051 + ,135 + ,1112 + ,103597 + ,379 + ,66 + ,27 + ,16 + ,60 + ,49 + ,23789 + ,44547 + ,26 + ,2694 + ,233328 + ,825 + ,165 + ,92 + ,28 + ,107 + ,102 + ,120733 + ,162647 + ,230 + ,1973 + ,256462 + ,798 + ,56 + ,123 + ,35 + ,127 + ,124 + ,105195 + ,174141 + ,181 + ,1769 + ,206161 + ,663 + ,61 + ,75 + ,28 + ,107 + ,99 + ,73107 + ,60622 + ,71 + ,3148 + ,311473 + ,1069 + ,53 + ,128 + ,38 + ,146 + ,129 + ,132068 + ,179566 + ,147 + ,2474 + ,235800 + ,921 + ,127 + ,105 + ,23 + ,84 + ,62 + ,149193 + ,184301 + ,190 + ,2084 + ,177939 + ,858 + ,63 + ,55 + ,36 + ,141 + ,73 + ,46821 + ,75661 + ,64 + ,1954 + ,207176 + ,711 + ,38 + ,56 + ,32 + ,123 + ,114 + ,87011 + ,96144 + ,105 + ,1226 + ,196553 + ,503 + ,50 + ,41 + ,29 + ,111 + ,99 + ,95260 + ,129847 + ,107 + ,1389 + ,174184 + ,382 + ,52 + ,72 + ,25 + ,98 + ,70 + ,55183 + ,117286 + ,94 + ,1496 + ,143246 + ,464 + ,42 + ,67 + ,27 + ,105 + ,104 + ,106671 + ,71180 + ,116 + ,2269 + ,187559 + ,717 + ,76 + ,75 + ,36 + ,135 + ,116 + ,73511 + ,109377 + ,106 + ,1833 + ,187681 + ,690 + ,67 + ,114 + ,28 + ,107 + ,91 + ,92945 + ,85298 + ,143 + ,1268 + ,119016 + ,462 + ,50 + ,118 + ,23 + ,85 + ,74 + ,78664 + ,73631 + ,81 + ,1943 + ,182192 + ,657 + ,53 + ,77 + ,40 + ,155 + ,138 + ,70054 + ,86767 + ,89 + ,893 + ,73566 + ,385 + ,39 + ,22 + ,23 + ,88 + ,67 + ,22618 + ,23824 + ,26 + ,1762 + ,194979 + ,577 + ,50 + ,66 + ,40 + ,155 + ,151 + ,74011 + ,93487 + ,84 + ,1403 + ,167488 + ,619 + ,77 + ,69 + ,28 + ,104 + ,72 + ,83737 + ,82981 + ,113 + ,1425 + ,143756 + ,479 + ,57 + ,105 + ,34 + ,132 + ,120 + ,69094 + ,73815 + ,120 + ,1857 + ,275541 + ,817 + ,73 + ,116 + ,33 + ,127 + ,115 + ,93133 + ,94552 + ,110 + ,1840 + ,243199 + ,752 + ,34 + ,88 + ,28 + ,108 + ,105 + ,95536 + ,132190 + ,134 + ,1502 + ,182999 + ,430 + ,39 + ,73 + ,34 + ,129 + ,104 + ,225920 + ,128754 + ,54 + ,1441 + ,135649 + ,451 + ,46 + ,99 + ,30 + ,116 + ,108 + ,62133 + ,66363 + ,96 + ,1420 + ,152299 + ,537 + ,63 + ,62 + ,33 + ,122 + ,98 + ,61370 + ,67808 + ,78 + ,1416 + ,120221 + ,519 + ,35 + ,53 + ,22 + ,85 + ,69 + ,43836 + ,61724 + ,51 + ,2970 + ,346485 + ,1000 + ,106 + ,118 + ,38 + ,147 + ,111 + ,106117 + ,131722 + ,121 + ,1317 + ,145790 + ,637 + ,43 + ,30 + ,26 + ,99 + ,99 + ,38692 + ,68580 + ,38 + ,1644 + ,193339 + ,465 + ,47 + ,100 + ,35 + ,87 + ,71 + ,84651 + ,106175 + ,145 + ,870 + ,80953 + ,437 + ,31 + ,49 + ,8 + ,28 + ,27 + ,56622 + ,55792 + ,59 + ,1654 + ,122774 + ,711 + ,162 + ,24 + ,24 + ,90 + ,69 + ,15986 + ,25157 + ,27 + ,1054 + ,130585 + ,299 + ,57 + ,67 + ,29 + ,109 + ,107 + ,95364 + ,76669 + ,91 + ,937 + ,112611 + ,248 + ,36 + ,46 + ,20 + ,78 + ,73 + ,26706 + ,57283 + ,48 + ,3004 + ,286468 + ,1162 + ,263 + ,57 + ,29 + ,111 + ,107 + ,89691 + ,105805 + ,68 + ,2008 + ,241066 + ,714 + ,78 + ,75 + ,45 + ,158 + ,93 + ,67267 + ,129484 + ,58 + ,2547 + ,148446 + ,905 + ,63 + ,135 + ,37 + ,141 + ,129 + ,126846 + ,72413 + ,150 + ,1885 + ,204713 + ,649 + ,54 + ,68 + ,33 + ,122 + ,69 + ,41140 + ,87831 + ,74 + ,1626 + ,182079 + ,512 + ,63 + ,124 + ,33 + ,124 + ,118 + ,102860 + ,96971 + ,181 + ,1468 + ,140344 + ,472 + ,77 + ,33 + ,25 + ,93 + ,73 + ,51715 + ,71299 + ,65 + ,2445 + ,220516 + ,905 + ,79 + ,98 + ,32 + ,124 + ,119 + ,55801 + ,77494 + ,97 + ,1964 + ,243060 + ,786 + ,110 + ,58 + ,29 + ,112 + ,104 + ,111813 + ,120336 + ,121 + ,1381 + ,162765 + ,489 + ,56 + ,68 + ,28 + ,108 + ,107 + ,120293 + ,93913 + ,99 + ,1369 + ,182613 + ,479 + ,56 + ,81 + ,28 + ,99 + ,99 + ,138599 + ,136048 + ,152 + ,1659 + ,232138 + ,617 + ,43 + ,131 + ,31 + ,117 + ,90 + ,161647 + ,181248 + ,188 + ,2888 + ,265318 + ,925 + ,111 + ,110 + ,52 + ,199 + ,197 + ,115929 + ,146123 + ,138 + ,1290 + ,85574 + ,351 + ,71 + ,37 + ,21 + ,78 + ,36 + ,24266 + ,32036 + ,40 + ,2845 + ,310839 + ,1144 + ,62 + ,130 + ,24 + ,91 + ,85 + ,162901 + ,186646 + ,254 + ,1982 + ,225060 + ,669 + ,56 + ,93 + ,41 + ,158 + ,139 + ,109825 + ,102255 + ,87 + ,1904 + ,232317 + ,707 + ,74 + ,118 + ,33 + ,126 + ,106 + ,129838 + ,168237 + ,178 + ,1391 + ,144966 + ,458 + ,60 + ,39 + ,32 + ,122 + ,50 + ,37510 + ,64219 + ,51 + ,602 + ,43287 + ,214 + ,43 + ,13 + ,19 + ,71 + ,64 + ,43750 + ,19630 + ,49 + ,1743 + ,155754 + ,599 + ,68 + ,74 + ,20 + ,75 + ,31 + ,40652 + ,76825 + ,73 + ,1559 + ,164709 + ,572 + ,53 + ,81 + ,31 + ,115 + ,63 + ,87771 + ,115338 + ,176 + ,2014 + ,201940 + ,897 + ,87 + ,109 + ,31 + ,119 + ,92 + ,85872 + ,109427 + ,94 + ,2143 + ,235454 + ,819 + ,46 + ,151 + ,32 + ,124 + ,106 + ,89275 + ,118168 + ,120 + ,2146 + ,220801 + ,720 + ,105 + ,51 + ,18 + ,72 + ,63 + ,44418 + ,84845 + ,66 + ,874 + ,99466 + ,273 + ,32 + ,28 + ,23 + ,91 + ,69 + ,192565 + ,153197 + ,56 + ,1590 + ,92661 + ,508 + ,133 + ,40 + ,17 + ,45 + ,41 + ,35232 + ,29877 + ,39 + ,1590 + ,133328 + ,506 + ,79 + ,56 + ,20 + ,78 + ,56 + ,40909 + ,63506 + ,66 + ,1210 + ,61361 + ,451 + ,51 + ,27 + ,12 + ,39 + ,25 + ,13294 + ,22445 + ,27 + ,2072 + ,125930 + ,699 + ,207 + ,37 + ,17 + ,68 + ,65 + ,32387 + ,47695 + ,65 + ,1281 + ,100750 + ,407 + ,67 + ,83 + ,30 + ,119 + ,93 + ,140867 + ,68370 + ,58 + ,1401 + ,224549 + ,465 + ,47 + ,54 + ,31 + ,117 + ,114 + ,120662 + ,146304 + ,98 + ,834 + ,82316 + ,245 + ,34 + ,27 + ,10 + ,39 + ,38 + ,21233 + ,38233 + ,25 + ,1105 + ,102010 + ,370 + ,66 + ,28 + ,13 + ,50 + ,44 + ,44332 + ,42071 + ,26 + ,1272 + ,101523 + ,316 + ,76 + ,59 + ,22 + ,88 + ,87 + ,61056 + ,50517 + ,77 + ,1944 + ,243511 + ,603 + ,65 + ,133 + ,42 + ,155 + ,110 + ,101338 + ,103950 + ,130 + ,391 + ,22938 + ,154 + ,9 + ,12 + ,1 + ,0 + ,0 + ,1168 + ,5841 + ,11 + ,761 + ,41566 + ,229 + ,42 + ,0 + ,9 + ,36 + ,27 + ,13497 + ,2341 + ,2 + ,1605 + ,152474 + ,577 + ,45 + ,106 + ,32 + ,123 + ,83 + ,65567 + ,84396 + ,101 + ,530 + ,61857 + ,192 + ,25 + ,23 + ,11 + ,32 + ,30 + ,25162 + ,24610 + ,31 + ,1988 + ,99923 + ,617 + ,115 + ,44 + ,25 + ,99 + ,80 + ,32334 + ,35753 + ,36 + ,1386 + ,132487 + ,411 + ,97 + ,71 + ,36 + ,136 + ,98 + ,40735 + ,55515 + ,120 + ,2395 + ,317394 + ,975 + ,53 + ,116 + ,31 + ,117 + ,82 + ,91413 + ,209056 + ,195 + ,387 + ,21054 + ,146 + ,2 + ,4 + ,0 + ,0 + ,0 + ,855 + ,6622 + ,4 + ,1742 + ,209641 + ,705 + ,52 + ,62 + ,24 + ,88 + ,60 + ,97068 + ,115814 + ,89 + ,620 + ,22648 + ,184 + ,44 + ,12 + ,13 + ,39 + ,28 + ,44339 + ,11609 + ,24 + ,449 + ,31414 + ,200 + ,22 + ,18 + ,8 + ,25 + ,9 + ,14116 + ,13155 + ,39 + ,800 + ,46698 + ,274 + ,35 + ,14 + ,13 + ,52 + ,33 + ,10288 + ,18274 + ,14 + ,1684 + ,131698 + ,502 + ,74 + ,60 + ,19 + ,75 + ,59 + ,65622 + ,72875 + ,78 + ,1050 + ,91735 + ,382 + ,103 + ,7 + ,18 + ,71 + ,49 + ,16563 + ,10112 + ,15 + ,2699 + ,244749 + ,964 + ,144 + ,98 + ,33 + ,124 + ,115 + ,76643 + ,142775 + ,106 + ,1606 + ,184510 + ,537 + ,60 + ,64 + ,40 + ,151 + ,140 + ,110681 + ,68847 + ,83 + ,1502 + ,79863 + ,438 + ,134 + ,29 + ,22 + ,71 + ,49 + ,29011 + ,17659 + ,24 + ,1204 + ,128423 + ,369 + ,89 + ,32 + ,38 + ,145 + ,120 + ,92696 + ,20112 + ,37 + ,1138 + ,97839 + ,417 + ,42 + ,25 + ,24 + ,87 + ,66 + ,94785 + ,61023 + ,77 + ,568 + ,38214 + ,276 + ,52 + ,16 + ,8 + ,27 + ,21 + ,8773 + ,13983 + ,16 + ,1459 + ,151101 + ,514 + ,98 + ,48 + ,35 + ,131 + ,124 + ,83209 + ,65176 + ,56 + ,2158 + ,272458 + ,822 + ,99 + ,100 + ,43 + ,162 + ,152 + ,93815 + ,132432 + ,132 + ,1111 + ,172494 + ,389 + ,52 + ,46 + ,43 + ,165 + ,139 + ,86687 + ,112494 + ,144 + ,1421 + ,108043 + ,466 + ,29 + ,45 + ,14 + ,54 + ,38 + ,34553 + ,45109 + ,40 + ,2833 + ,328107 + ,1255 + ,125 + ,129 + ,41 + ,159 + ,144 + ,105547 + ,170875 + ,153 + ,1955 + ,250579 + ,694 + ,106 + ,130 + ,38 + ,147 + ,120 + ,103487 + ,180759 + ,143 + ,2922 + ,351067 + ,1024 + ,95 + ,136 + ,45 + ,170 + ,160 + ,213688 + ,214921 + ,220 + ,1002 + ,158015 + ,400 + ,40 + ,59 + ,31 + ,119 + ,114 + ,71220 + ,100226 + ,79 + ,1060 + ,98866 + ,397 + ,140 + ,25 + ,13 + ,49 + ,39 + ,23517 + ,32043 + ,50 + ,956 + ,85439 + ,350 + ,43 + ,32 + ,28 + ,104 + ,78 + ,56926 + ,54454 + ,39 + ,2186 + ,229242 + ,719 + ,128 + ,63 + ,31 + ,120 + ,119 + ,91721 + ,78876 + ,95 + ,3604 + ,351619 + ,1277 + ,142 + ,95 + ,40 + ,150 + ,141 + ,115168 + ,170745 + ,169 + ,1035 + ,84207 + ,356 + ,73 + ,14 + ,30 + ,112 + ,101 + ,111194 + ,6940 + ,12 + ,1417 + ,120445 + ,457 + ,72 + ,36 + ,16 + ,59 + ,56 + ,51009 + ,49025 + ,63 + ,3261 + ,324598 + ,1402 + ,128 + ,113 + ,37 + ,136 + ,133 + ,135777 + ,122037 + ,134 + ,1587 + ,131069 + ,600 + ,61 + ,47 + ,30 + ,107 + ,83 + ,51513 + ,53782 + ,69 + ,1424 + ,204271 + ,480 + ,73 + ,92 + ,35 + ,130 + ,116 + ,74163 + ,127748 + ,119 + ,1701 + ,165543 + ,595 + ,148 + ,70 + ,32 + ,115 + ,90 + ,51633 + ,86839 + ,119 + ,1249 + ,141722 + ,436 + ,64 + ,19 + ,27 + ,107 + ,36 + ,75345 + ,44830 + ,75 + ,946 + ,116048 + ,230 + ,45 + ,50 + ,20 + ,75 + ,50 + ,33416 + ,77395 + ,63 + ,1926 + ,250047 + ,651 + ,58 + ,41 + ,18 + ,71 + ,61 + ,83305 + ,89324 + ,55 + ,3352 + ,299775 + ,1367 + ,97 + ,91 + ,31 + ,120 + ,97 + ,98952 + ,103300 + ,103 + ,1641 + ,195838 + ,564 + ,50 + ,111 + ,31 + ,116 + ,98 + ,102372 + ,112283 + ,197 + ,2035 + ,173260 + ,716 + ,37 + ,41 + ,21 + ,79 + ,78 + ,37238 + ,10901 + ,16 + ,2312 + ,254488 + ,747 + ,50 + ,120 + ,39 + ,150 + ,117 + ,103772 + ,120691 + ,140 + ,1369 + ,104389 + ,467 + ,105 + ,135 + ,41 + ,156 + ,148 + ,123969 + ,58106 + ,89 + ,1577 + ,136084 + ,671 + ,69 + ,27 + ,13 + ,51 + ,41 + ,27142 + ,57140 + ,40 + ,2201 + ,199476 + ,861 + ,46 + ,87 + ,32 + ,118 + ,105 + ,135400 + ,122422 + ,125 + ,961 + ,92499 + ,319 + ,57 + ,25 + ,18 + ,71 + ,55 + ,21399 + ,25899 + ,21 + ,1900 + ,224330 + ,612 + ,52 + ,131 + ,39 + ,144 + ,132 + ,130115 + ,139296 + ,167 + ,1254 + ,135781 + ,433 + ,98 + ,45 + ,14 + ,47 + ,44 + ,24874 + ,52678 + ,32 + ,1335 + ,74408 + ,434 + ,61 + ,29 + ,7 + ,28 + ,21 + ,34988 + ,23853 + ,36 + ,1597 + ,81240 + ,503 + ,89 + ,58 + ,17 + ,68 + ,50 + ,45549 + ,17306 + ,13 + ,207 + ,14688 + ,85 + ,0 + ,4 + ,0 + ,0 + ,0 + ,6023 + ,7953 + ,5 + ,1645 + ,181633 + ,564 + ,48 + ,47 + ,30 + ,110 + ,73 + ,64466 + ,89455 + ,96 + ,2429 + ,271856 + ,824 + ,91 + ,109 + ,37 + ,147 + ,86 + ,54990 + ,147866 + ,151 + ,151 + ,7199 + ,74 + ,0 + ,7 + ,0 + ,0 + ,0 + ,1644 + ,4245 + ,6 + ,474 + ,46660 + ,259 + ,7 + ,12 + ,5 + ,15 + ,13 + ,6179 + ,21509 + ,13 + ,141 + ,17547 + ,69 + ,3 + ,0 + ,1 + ,4 + ,4 + ,3926 + ,7670 + ,3 + ,1639 + ,133368 + ,535 + ,54 + ,37 + ,16 + ,64 + ,57 + ,32755 + ,66675 + ,57 + ,872 + ,95227 + ,239 + ,70 + ,37 + ,32 + ,111 + ,48 + ,34777 + ,14336 + ,23 + ,1318 + ,152601 + ,438 + ,36 + ,46 + ,24 + ,85 + ,46 + ,73224 + ,53608 + ,61 + ,1018 + ,98146 + ,459 + ,37 + ,15 + ,17 + ,68 + ,48 + ,27114 + ,30059 + ,21 + ,1383 + ,79619 + ,426 + ,123 + ,42 + ,11 + ,40 + ,32 + ,20760 + ,29668 + ,43 + ,1314 + ,59194 + ,288 + ,247 + ,7 + ,24 + ,80 + ,68 + ,37636 + ,22097 + ,20 + ,1335 + ,139942 + ,498 + ,46 + ,54 + ,22 + ,88 + ,87 + ,65461 + ,96841 + ,82 + ,1403 + ,118612 + ,454 + ,72 + ,54 + ,12 + ,48 + ,43 + ,30080 + ,41907 + ,90 + ,910 + ,72880 + ,376 + ,41 + ,14 + ,19 + ,76 + ,67 + ,24094 + ,27080 + ,25 + ,616 + ,65475 + ,225 + ,24 + ,16 + ,13 + ,51 + ,46 + ,69008 + ,35885 + ,60 + ,1407 + ,99643 + ,555 + ,45 + ,33 + ,17 + ,67 + ,46 + ,54968 + ,41247 + ,61 + ,771 + ,71965 + ,252 + ,33 + ,32 + ,15 + ,59 + ,56 + ,46090 + ,28313 + ,85 + ,766 + ,77272 + ,208 + ,27 + ,21 + ,16 + ,61 + ,48 + ,27507 + ,36845 + ,43 + ,473 + ,49289 + ,130 + ,36 + ,15 + ,24 + ,76 + ,44 + ,10672 + ,16548 + ,25 + ,1376 + ,135131 + ,481 + ,87 + ,38 + ,15 + ,60 + ,60 + ,34029 + ,36134 + ,41 + ,1232 + ,108446 + ,389 + ,90 + ,22 + ,17 + ,68 + ,65 + ,46300 + ,55764 + ,26 + ,1521 + ,89746 + ,565 + ,114 + ,28 + ,18 + ,71 + ,55 + ,24760 + ,28910 + ,38 + ,572 + ,44296 + ,173 + ,31 + ,10 + ,20 + ,76 + ,38 + ,18779 + ,13339 + ,12 + ,1059 + ,77648 + ,278 + ,45 + ,31 + ,16 + ,62 + ,52 + ,21280 + ,25319 + ,29 + ,1544 + ,181528 + ,609 + ,69 + ,32 + ,16 + ,61 + ,60 + ,40662 + ,66956 + ,49 + ,1230 + ,134019 + ,422 + ,51 + ,32 + ,18 + ,67 + ,54 + ,28987 + ,47487 + ,46 + ,1206 + ,124064 + ,445 + ,34 + ,43 + ,22 + ,88 + ,86 + ,22827 + ,52785 + ,41 + ,1205 + ,92630 + ,387 + ,60 + ,27 + ,8 + ,30 + ,24 + ,18513 + ,44683 + ,31 + ,1255 + ,121848 + ,339 + ,45 + ,37 + ,17 + ,64 + ,52 + ,30594 + ,35619 + ,41 + ,613 + ,52915 + ,181 + ,54 + ,20 + ,18 + ,68 + ,49 + ,24006 + ,21920 + ,26 + ,721 + ,81872 + ,245 + ,25 + ,32 + ,16 + ,64 + ,61 + ,27913 + ,45608 + ,23 + ,1109 + ,58981 + ,384 + ,38 + ,0 + ,23 + ,91 + ,61 + ,42744 + ,7721 + ,14 + ,740 + ,53515 + ,212 + ,52 + ,5 + ,22 + ,88 + ,81 + ,12934 + ,20634 + ,16 + ,1126 + ,60812 + ,399 + ,67 + ,26 + ,13 + ,52 + ,43 + ,22574 + ,29788 + ,25 + ,728 + ,56375 + ,229 + ,74 + ,10 + ,13 + ,49 + ,40 + ,41385 + ,31931 + ,21 + ,689 + ,65490 + ,224 + ,38 + ,27 + ,16 + ,62 + ,40 + ,18653 + ,37754 + ,32 + ,592 + ,80949 + ,203 + ,30 + ,11 + ,16 + ,61 + ,56 + ,18472 + ,32505 + ,9 + ,995 + ,76302 + ,333 + ,26 + ,29 + ,20 + ,76 + ,68 + ,30976 + ,40557 + ,35 + ,1613 + ,104011 + ,384 + ,67 + ,25 + ,22 + ,88 + ,79 + ,63339 + ,94238 + ,42 + ,2048 + ,98104 + ,636 + ,132 + ,55 + ,17 + ,66 + ,47 + ,25568 + ,44197 + ,68 + ,705 + ,67989 + ,185 + ,42 + ,23 + ,18 + ,71 + ,57 + ,33747 + ,43228 + ,32 + ,301 + ,30989 + ,93 + ,35 + ,5 + ,17 + ,68 + ,41 + ,4154 + ,4103 + ,6 + ,1803 + ,135458 + ,581 + ,118 + ,43 + ,12 + ,48 + ,29 + ,19474 + ,44144 + ,68 + ,799 + ,73504 + ,248 + ,68 + ,23 + ,7 + ,25 + ,3 + ,35130 + ,32868 + ,33 + ,861 + ,63123 + ,304 + ,43 + ,34 + ,17 + ,68 + ,60 + ,39067 + ,27640 + ,84 + ,1186 + ,61254 + ,344 + ,76 + ,36 + ,14 + ,41 + ,30 + ,13310 + ,14063 + ,46 + ,1451 + ,74914 + ,407 + ,64 + ,35 + ,23 + ,90 + ,79 + ,65892 + ,28990 + ,30 + ,628 + ,31774 + ,170 + ,48 + ,0 + ,17 + ,66 + ,47 + ,4143 + ,4694 + ,0 + ,1161 + ,81437 + ,312 + ,64 + ,37 + ,14 + ,54 + ,40 + ,28579 + ,42648 + ,36 + ,1463 + ,87186 + ,507 + ,56 + ,28 + ,15 + ,59 + ,48 + ,51776 + ,64329 + ,47 + ,742 + ,50090 + ,224 + ,71 + ,16 + ,17 + ,60 + ,36 + ,21152 + ,21928 + ,20 + ,979 + ,65745 + ,340 + ,75 + ,26 + ,21 + ,77 + ,42 + ,38084 + ,25836 + ,50 + ,675 + ,56653 + ,168 + ,39 + ,38 + ,18 + ,68 + ,49 + ,27717 + ,22779 + ,30 + ,1241 + ,158399 + ,443 + ,42 + ,23 + ,18 + ,72 + ,57 + ,32928 + ,40820 + ,30 + ,676 + ,46455 + ,204 + ,39 + ,22 + ,17 + ,67 + ,12 + ,11342 + ,27530 + ,34 + ,1049 + ,73624 + ,367 + ,93 + ,30 + ,17 + ,64 + ,40 + ,19499 + ,32378 + ,33 + ,620 + ,38395 + ,210 + ,38 + ,16 + ,16 + ,63 + ,43 + ,16380 + ,10824 + ,34 + ,1081 + ,91899 + ,335 + ,60 + ,18 + ,15 + ,59 + ,33 + ,36874 + ,39613 + ,37 + ,1688 + ,139526 + ,364 + ,71 + ,28 + ,21 + ,84 + ,77 + ,48259 + ,60865 + ,83 + ,736 + ,52164 + ,178 + ,52 + ,32 + ,16 + ,64 + ,43 + ,16734 + ,19787 + ,32 + ,617 + ,51567 + ,206 + ,27 + ,21 + ,14 + ,56 + ,45 + ,28207 + ,20107 + ,30 + ,812 + ,70551 + ,279 + ,59 + ,23 + ,15 + ,54 + ,47 + ,30143 + ,36605 + ,43 + ,1051 + ,84856 + ,387 + ,40 + ,29 + ,17 + ,67 + ,43 + ,41369 + ,40961 + ,41 + ,1656 + ,102538 + ,490 + ,79 + ,50 + ,15 + ,58 + ,45 + ,45833 + ,48231 + ,51 + ,705 + ,86678 + ,238 + ,44 + ,12 + ,15 + ,59 + ,50 + ,29156 + ,39725 + ,19 + ,945 + ,85709 + ,343 + ,65 + ,21 + ,10 + ,40 + ,35 + ,35944 + ,21455 + ,37 + ,554 + ,34662 + ,232 + ,10 + ,18 + ,6 + ,22 + ,7 + ,36278 + ,23430 + ,33 + ,1597 + ,150580 + ,530 + ,124 + ,27 + ,22 + ,83 + ,71 + ,45588 + ,62991 + ,41 + ,982 + ,99611 + ,291 + ,81 + ,41 + ,21 + ,81 + ,67 + ,45097 + ,49363 + ,54 + ,222 + ,19349 + ,67 + ,15 + ,13 + ,1 + ,2 + ,0 + ,3895 + ,9604 + ,14 + ,1212 + ,99373 + ,397 + ,92 + ,12 + ,18 + ,72 + ,62 + ,28394 + ,24552 + ,25 + ,1143 + ,86230 + ,467 + ,42 + ,21 + ,17 + ,61 + ,54 + ,18632 + ,31493 + ,25 + ,435 + ,30837 + ,178 + ,10 + ,8 + ,4 + ,15 + ,4 + ,2325 + ,3439 + ,8 + ,532 + ,31706 + ,175 + ,24 + ,26 + ,10 + ,32 + ,25 + ,25139 + ,19555 + ,26 + ,882 + ,89806 + ,299 + ,64 + ,27 + ,16 + ,62 + ,40 + ,27975 + ,21228 + ,20 + ,608 + ,62088 + ,154 + ,45 + ,13 + ,16 + ,58 + ,38 + ,14483 + ,23177 + ,11 + ,459 + ,40151 + ,106 + ,22 + ,16 + ,9 + ,36 + ,19 + ,13127 + ,22094 + ,14 + ,578 + ,27634 + ,189 + ,56 + ,2 + ,16 + ,59 + ,17 + ,5839 + ,2342 + ,3 + ,826 + ,76990 + ,194 + ,94 + ,42 + ,17 + ,68 + ,67 + ,24069 + ,38798 + ,40 + ,509 + ,37460 + ,135 + ,19 + ,5 + ,7 + ,21 + ,14 + ,3738 + ,3255 + ,5 + ,717 + ,54157 + ,201 + ,35 + ,37 + ,15 + ,55 + ,30 + ,18625 + ,24261 + ,38 + ,637 + ,49862 + ,207 + ,32 + ,17 + ,14 + ,54 + ,54 + ,36341 + ,18511 + ,32 + ,857 + ,84337 + ,280 + ,35 + ,38 + ,14 + ,55 + ,35 + ,24548 + ,40798 + ,41 + ,830 + ,64175 + ,260 + ,48 + ,37 + ,18 + ,72 + ,59 + ,21792 + ,28893 + ,46 + ,652 + ,59382 + ,227 + ,49 + ,29 + ,12 + ,41 + ,24 + ,26263 + ,21425 + ,47 + ,707 + ,119308 + ,239 + ,48 + ,32 + ,16 + ,61 + ,58 + ,23686 + ,50276 + ,37 + ,954 + ,76702 + ,333 + ,62 + ,35 + ,21 + ,67 + ,42 + ,49303 + ,37643 + ,51 + ,1461 + ,103425 + ,428 + ,96 + ,17 + ,19 + ,76 + ,46 + ,25659 + ,30377 + ,49 + ,672 + ,70344 + ,230 + ,45 + ,20 + ,16 + ,64 + ,61 + ,28904 + ,27126 + ,21 + ,778 + ,43410 + ,292 + ,63 + ,7 + ,1 + ,3 + ,3 + ,2781 + ,13 + ,1 + ,1141 + ,104838 + ,350 + ,71 + ,46 + ,16 + ,63 + ,52 + ,29236 + ,42097 + ,44 + ,680 + ,62215 + ,186 + ,26 + ,24 + ,10 + ,40 + ,25 + ,19546 + ,24451 + ,26 + ,1090 + ,69304 + ,326 + ,48 + ,40 + ,19 + ,69 + ,40 + ,22818 + ,14335 + ,21 + ,616 + ,53117 + ,155 + ,29 + ,3 + ,12 + ,48 + ,32 + ,32689 + ,5084 + ,4 + ,285 + ,19764 + ,75 + ,19 + ,10 + ,2 + ,8 + ,4 + ,5752 + ,9927 + ,10 + ,1145 + ,86680 + ,361 + ,45 + ,37 + ,14 + ,52 + ,49 + ,22197 + ,43527 + ,43 + ,733 + ,84105 + ,261 + ,45 + ,17 + ,17 + ,66 + ,63 + ,20055 + ,27184 + ,34 + ,888 + ,77945 + ,299 + ,67 + ,28 + ,19 + ,76 + ,67 + ,25272 + ,21610 + ,32 + ,849 + ,89113 + ,300 + ,30 + ,19 + ,14 + ,43 + ,32 + ,82206 + ,20484 + ,20 + ,1182 + ,91005 + ,450 + ,36 + ,29 + ,11 + ,39 + ,23 + ,32073 + ,20156 + ,34 + ,528 + ,40248 + ,183 + ,34 + ,8 + ,4 + ,14 + ,7 + ,5444 + ,6012 + ,6 + ,642 + ,64187 + ,238 + ,36 + ,10 + ,16 + ,61 + ,54 + ,20154 + ,18475 + ,12 + ,947 + ,50857 + ,165 + ,34 + ,15 + ,20 + ,71 + ,37 + ,36944 + ,12645 + ,24 + ,819 + ,56613 + ,234 + ,37 + ,15 + ,12 + ,44 + ,35 + ,8019 + ,11017 + ,16 + ,757 + ,62792 + ,176 + ,46 + ,28 + ,15 + ,60 + ,51 + ,30884 + ,37623 + ,72 + ,894 + ,72535 + ,329 + ,44 + ,17 + ,16 + ,64 + ,39 + ,19540 + ,35873 + ,27) + ,dim=c(11 + ,289) + ,dimnames=list(c('pageviews' + ,'time_in_rfc' + ,'compendium_views_info' + ,'compendium_views_pr' + ,'blogged_computations' + ,'compendiums_reviewed' + ,'feedback_messages_p1' + ,'feedback_messages_p120' + ,'totsize' + ,'totseconds' + ,'tothyperlinks') + ,1:289)) > y <- array(NA,dim=c(11,289),dimnames=list(c('pageviews','time_in_rfc','compendium_views_info','compendium_views_pr','blogged_computations','compendiums_reviewed','feedback_messages_p1','feedback_messages_p120','totsize','totseconds','tothyperlinks'),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' > par3 <- 'No Linear Trend' > par2 <- 'Do not include Seasonal Dummies' > par1 <- '1' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from 'package:base': as.Date, 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 pageviews time_in_rfc compendium_views_info compendium_views_pr 1 1418 210907 396 81 2 869 120982 297 55 3 1530 176508 559 50 4 2172 179321 967 125 5 901 123185 270 40 6 463 52746 143 37 7 3201 385534 1562 63 8 371 33170 109 44 9 1192 101645 371 88 10 1583 149061 656 66 11 1439 165446 511 57 12 1764 237213 655 74 13 1495 173326 465 49 14 1373 133131 525 52 15 2187 258873 885 88 16 1491 180083 497 36 17 4041 324799 1436 108 18 1706 230964 612 43 19 2152 236785 865 75 20 1036 135473 385 32 21 1882 202925 567 44 22 1929 215147 639 85 23 2242 344297 963 86 24 1220 153935 398 56 25 1289 132943 410 50 26 2515 174724 966 135 27 2147 174415 801 63 28 2352 225548 892 81 29 1638 223632 513 52 30 1222 124817 469 44 31 1812 221698 683 113 32 1677 210767 643 39 33 1579 170266 535 73 34 1731 260561 625 48 35 807 84853 264 33 36 2452 294424 992 59 37 829 101011 238 41 38 1940 215641 818 69 39 2662 325107 937 64 40 186 7176 70 1 41 1499 167542 507 59 42 865 106408 260 32 43 1793 96560 503 129 44 2527 265769 927 37 45 2747 269651 1269 31 46 1324 149112 537 65 47 2702 175824 910 107 48 1383 152871 532 74 49 1179 111665 345 54 50 2099 116408 918 76 51 4308 362301 1635 715 52 918 78800 330 57 53 1831 183167 557 66 54 3373 277965 1178 106 55 1713 150629 740 54 56 1438 168809 452 32 57 496 24188 218 20 58 2253 329267 764 71 59 744 65029 255 21 60 1161 101097 454 70 61 2352 218946 866 112 62 2144 244052 574 66 63 4691 341570 1276 190 64 1112 103597 379 66 65 2694 233328 825 165 66 1973 256462 798 56 67 1769 206161 663 61 68 3148 311473 1069 53 69 2474 235800 921 127 70 2084 177939 858 63 71 1954 207176 711 38 72 1226 196553 503 50 73 1389 174184 382 52 74 1496 143246 464 42 75 2269 187559 717 76 76 1833 187681 690 67 77 1268 119016 462 50 78 1943 182192 657 53 79 893 73566 385 39 80 1762 194979 577 50 81 1403 167488 619 77 82 1425 143756 479 57 83 1857 275541 817 73 84 1840 243199 752 34 85 1502 182999 430 39 86 1441 135649 451 46 87 1420 152299 537 63 88 1416 120221 519 35 89 2970 346485 1000 106 90 1317 145790 637 43 91 1644 193339 465 47 92 870 80953 437 31 93 1654 122774 711 162 94 1054 130585 299 57 95 937 112611 248 36 96 3004 286468 1162 263 97 2008 241066 714 78 98 2547 148446 905 63 99 1885 204713 649 54 100 1626 182079 512 63 101 1468 140344 472 77 102 2445 220516 905 79 103 1964 243060 786 110 104 1381 162765 489 56 105 1369 182613 479 56 106 1659 232138 617 43 107 2888 265318 925 111 108 1290 85574 351 71 109 2845 310839 1144 62 110 1982 225060 669 56 111 1904 232317 707 74 112 1391 144966 458 60 113 602 43287 214 43 114 1743 155754 599 68 115 1559 164709 572 53 116 2014 201940 897 87 117 2143 235454 819 46 118 2146 220801 720 105 119 874 99466 273 32 120 1590 92661 508 133 121 1590 133328 506 79 122 1210 61361 451 51 123 2072 125930 699 207 124 1281 100750 407 67 125 1401 224549 465 47 126 834 82316 245 34 127 1105 102010 370 66 128 1272 101523 316 76 129 1944 243511 603 65 130 391 22938 154 9 131 761 41566 229 42 132 1605 152474 577 45 133 530 61857 192 25 134 1988 99923 617 115 135 1386 132487 411 97 136 2395 317394 975 53 137 387 21054 146 2 138 1742 209641 705 52 139 620 22648 184 44 140 449 31414 200 22 141 800 46698 274 35 142 1684 131698 502 74 143 1050 91735 382 103 144 2699 244749 964 144 145 1606 184510 537 60 146 1502 79863 438 134 147 1204 128423 369 89 148 1138 97839 417 42 149 568 38214 276 52 150 1459 151101 514 98 151 2158 272458 822 99 152 1111 172494 389 52 153 1421 108043 466 29 154 2833 328107 1255 125 155 1955 250579 694 106 156 2922 351067 1024 95 157 1002 158015 400 40 158 1060 98866 397 140 159 956 85439 350 43 160 2186 229242 719 128 161 3604 351619 1277 142 162 1035 84207 356 73 163 1417 120445 457 72 164 3261 324598 1402 128 165 1587 131069 600 61 166 1424 204271 480 73 167 1701 165543 595 148 168 1249 141722 436 64 169 946 116048 230 45 170 1926 250047 651 58 171 3352 299775 1367 97 172 1641 195838 564 50 173 2035 173260 716 37 174 2312 254488 747 50 175 1369 104389 467 105 176 1577 136084 671 69 177 2201 199476 861 46 178 961 92499 319 57 179 1900 224330 612 52 180 1254 135781 433 98 181 1335 74408 434 61 182 1597 81240 503 89 183 207 14688 85 0 184 1645 181633 564 48 185 2429 271856 824 91 186 151 7199 74 0 187 474 46660 259 7 188 141 17547 69 3 189 1639 133368 535 54 190 872 95227 239 70 191 1318 152601 438 36 192 1018 98146 459 37 193 1383 79619 426 123 194 1314 59194 288 247 195 1335 139942 498 46 196 1403 118612 454 72 197 910 72880 376 41 198 616 65475 225 24 199 1407 99643 555 45 200 771 71965 252 33 201 766 77272 208 27 202 473 49289 130 36 203 1376 135131 481 87 204 1232 108446 389 90 205 1521 89746 565 114 206 572 44296 173 31 207 1059 77648 278 45 208 1544 181528 609 69 209 1230 134019 422 51 210 1206 124064 445 34 211 1205 92630 387 60 212 1255 121848 339 45 213 613 52915 181 54 214 721 81872 245 25 215 1109 58981 384 38 216 740 53515 212 52 217 1126 60812 399 67 218 728 56375 229 74 219 689 65490 224 38 220 592 80949 203 30 221 995 76302 333 26 222 1613 104011 384 67 223 2048 98104 636 132 224 705 67989 185 42 225 301 30989 93 35 226 1803 135458 581 118 227 799 73504 248 68 228 861 63123 304 43 229 1186 61254 344 76 230 1451 74914 407 64 231 628 31774 170 48 232 1161 81437 312 64 233 1463 87186 507 56 234 742 50090 224 71 235 979 65745 340 75 236 675 56653 168 39 237 1241 158399 443 42 238 676 46455 204 39 239 1049 73624 367 93 240 620 38395 210 38 241 1081 91899 335 60 242 1688 139526 364 71 243 736 52164 178 52 244 617 51567 206 27 245 812 70551 279 59 246 1051 84856 387 40 247 1656 102538 490 79 248 705 86678 238 44 249 945 85709 343 65 250 554 34662 232 10 251 1597 150580 530 124 252 982 99611 291 81 253 222 19349 67 15 254 1212 99373 397 92 255 1143 86230 467 42 256 435 30837 178 10 257 532 31706 175 24 258 882 89806 299 64 259 608 62088 154 45 260 459 40151 106 22 261 578 27634 189 56 262 826 76990 194 94 263 509 37460 135 19 264 717 54157 201 35 265 637 49862 207 32 266 857 84337 280 35 267 830 64175 260 48 268 652 59382 227 49 269 707 119308 239 48 270 954 76702 333 62 271 1461 103425 428 96 272 672 70344 230 45 273 778 43410 292 63 274 1141 104838 350 71 275 680 62215 186 26 276 1090 69304 326 48 277 616 53117 155 29 278 285 19764 75 19 279 1145 86680 361 45 280 733 84105 261 45 281 888 77945 299 67 282 849 89113 300 30 283 1182 91005 450 36 284 528 40248 183 34 285 642 64187 238 36 286 947 50857 165 34 287 819 56613 234 37 288 757 62792 176 46 289 894 72535 329 44 blogged_computations compendiums_reviewed feedback_messages_p1 1 79 30 115 2 58 28 109 3 60 38 146 4 108 30 116 5 49 22 68 6 0 26 101 7 121 25 96 8 1 18 67 9 20 11 44 10 43 26 100 11 69 25 93 12 78 38 140 13 86 44 166 14 44 30 99 15 104 40 139 16 63 34 130 17 158 47 181 18 102 30 116 19 77 31 116 20 82 23 88 21 115 36 139 22 101 36 135 23 80 30 108 24 50 25 89 25 83 39 156 26 123 34 129 27 73 31 118 28 81 31 118 29 105 33 125 30 47 25 95 31 105 33 126 32 94 35 135 33 44 42 154 34 114 43 165 35 38 30 113 36 107 33 127 37 30 13 52 38 71 32 121 39 84 36 136 40 0 0 0 41 59 28 108 42 33 14 46 43 42 17 54 44 96 32 124 45 106 30 115 46 56 35 128 47 57 20 80 48 59 28 97 49 39 28 104 50 34 39 59 51 76 34 125 52 20 26 82 53 91 39 149 54 115 39 149 55 85 33 122 56 76 28 118 57 8 4 12 58 79 39 144 59 21 18 67 60 30 14 52 61 76 29 108 62 101 44 166 63 94 21 80 64 27 16 60 65 92 28 107 66 123 35 127 67 75 28 107 68 128 38 146 69 105 23 84 70 55 36 141 71 56 32 123 72 41 29 111 73 72 25 98 74 67 27 105 75 75 36 135 76 114 28 107 77 118 23 85 78 77 40 155 79 22 23 88 80 66 40 155 81 69 28 104 82 105 34 132 83 116 33 127 84 88 28 108 85 73 34 129 86 99 30 116 87 62 33 122 88 53 22 85 89 118 38 147 90 30 26 99 91 100 35 87 92 49 8 28 93 24 24 90 94 67 29 109 95 46 20 78 96 57 29 111 97 75 45 158 98 135 37 141 99 68 33 122 100 124 33 124 101 33 25 93 102 98 32 124 103 58 29 112 104 68 28 108 105 81 28 99 106 131 31 117 107 110 52 199 108 37 21 78 109 130 24 91 110 93 41 158 111 118 33 126 112 39 32 122 113 13 19 71 114 74 20 75 115 81 31 115 116 109 31 119 117 151 32 124 118 51 18 72 119 28 23 91 120 40 17 45 121 56 20 78 122 27 12 39 123 37 17 68 124 83 30 119 125 54 31 117 126 27 10 39 127 28 13 50 128 59 22 88 129 133 42 155 130 12 1 0 131 0 9 36 132 106 32 123 133 23 11 32 134 44 25 99 135 71 36 136 136 116 31 117 137 4 0 0 138 62 24 88 139 12 13 39 140 18 8 25 141 14 13 52 142 60 19 75 143 7 18 71 144 98 33 124 145 64 40 151 146 29 22 71 147 32 38 145 148 25 24 87 149 16 8 27 150 48 35 131 151 100 43 162 152 46 43 165 153 45 14 54 154 129 41 159 155 130 38 147 156 136 45 170 157 59 31 119 158 25 13 49 159 32 28 104 160 63 31 120 161 95 40 150 162 14 30 112 163 36 16 59 164 113 37 136 165 47 30 107 166 92 35 130 167 70 32 115 168 19 27 107 169 50 20 75 170 41 18 71 171 91 31 120 172 111 31 116 173 41 21 79 174 120 39 150 175 135 41 156 176 27 13 51 177 87 32 118 178 25 18 71 179 131 39 144 180 45 14 47 181 29 7 28 182 58 17 68 183 4 0 0 184 47 30 110 185 109 37 147 186 7 0 0 187 12 5 15 188 0 1 4 189 37 16 64 190 37 32 111 191 46 24 85 192 15 17 68 193 42 11 40 194 7 24 80 195 54 22 88 196 54 12 48 197 14 19 76 198 16 13 51 199 33 17 67 200 32 15 59 201 21 16 61 202 15 24 76 203 38 15 60 204 22 17 68 205 28 18 71 206 10 20 76 207 31 16 62 208 32 16 61 209 32 18 67 210 43 22 88 211 27 8 30 212 37 17 64 213 20 18 68 214 32 16 64 215 0 23 91 216 5 22 88 217 26 13 52 218 10 13 49 219 27 16 62 220 11 16 61 221 29 20 76 222 25 22 88 223 55 17 66 224 23 18 71 225 5 17 68 226 43 12 48 227 23 7 25 228 34 17 68 229 36 14 41 230 35 23 90 231 0 17 66 232 37 14 54 233 28 15 59 234 16 17 60 235 26 21 77 236 38 18 68 237 23 18 72 238 22 17 67 239 30 17 64 240 16 16 63 241 18 15 59 242 28 21 84 243 32 16 64 244 21 14 56 245 23 15 54 246 29 17 67 247 50 15 58 248 12 15 59 249 21 10 40 250 18 6 22 251 27 22 83 252 41 21 81 253 13 1 2 254 12 18 72 255 21 17 61 256 8 4 15 257 26 10 32 258 27 16 62 259 13 16 58 260 16 9 36 261 2 16 59 262 42 17 68 263 5 7 21 264 37 15 55 265 17 14 54 266 38 14 55 267 37 18 72 268 29 12 41 269 32 16 61 270 35 21 67 271 17 19 76 272 20 16 64 273 7 1 3 274 46 16 63 275 24 10 40 276 40 19 69 277 3 12 48 278 10 2 8 279 37 14 52 280 17 17 66 281 28 19 76 282 19 14 43 283 29 11 39 284 8 4 14 285 10 16 61 286 15 20 71 287 15 12 44 288 28 15 60 289 17 16 64 feedback_messages_p120 totsize totseconds tothyperlinks 1 94 112285 146283 144 2 103 84786 98364 103 3 93 83123 86146 98 4 103 101193 96933 135 5 51 38361 79234 61 6 70 68504 42551 39 7 91 119182 195663 150 8 22 22807 6853 5 9 38 17140 21529 28 10 93 116174 95757 84 11 60 57635 85584 80 12 123 66198 143983 130 13 148 71701 75851 82 14 90 57793 59238 60 15 124 80444 93163 131 16 70 53855 96037 84 17 168 97668 151511 140 18 115 133824 136368 151 19 71 101481 112642 91 20 66 99645 94728 138 21 134 114789 105499 150 22 117 99052 121527 124 23 108 67654 127766 119 24 84 65553 98958 73 25 156 97500 77900 110 26 120 69112 85646 123 27 114 82753 98579 90 28 94 85323 130767 116 29 120 72654 131741 113 30 81 30727 53907 56 31 110 77873 178812 115 32 133 117478 146761 119 33 122 74007 82036 129 34 158 90183 163253 127 35 109 61542 27032 27 36 124 101494 171975 175 37 39 27570 65990 35 38 92 55813 86572 64 39 126 79215 159676 96 40 0 1423 1929 0 41 70 55461 85371 84 42 37 31081 58391 41 43 38 22996 31580 47 44 120 83122 136815 126 45 93 70106 120642 105 46 95 60578 69107 80 47 77 39992 50495 70 48 90 79892 108016 73 49 80 49810 46341 57 50 31 71570 78348 40 51 110 100708 79336 68 52 66 33032 56968 21 53 138 82875 93176 127 54 133 139077 161632 154 55 113 71595 87850 116 56 100 72260 127969 102 57 7 5950 15049 7 58 140 115762 155135 148 59 61 32551 25109 21 60 41 31701 45824 35 61 96 80670 102996 112 62 164 143558 160604 137 63 78 117105 158051 135 64 49 23789 44547 26 65 102 120733 162647 230 66 124 105195 174141 181 67 99 73107 60622 71 68 129 132068 179566 147 69 62 149193 184301 190 70 73 46821 75661 64 71 114 87011 96144 105 72 99 95260 129847 107 73 70 55183 117286 94 74 104 106671 71180 116 75 116 73511 109377 106 76 91 92945 85298 143 77 74 78664 73631 81 78 138 70054 86767 89 79 67 22618 23824 26 80 151 74011 93487 84 81 72 83737 82981 113 82 120 69094 73815 120 83 115 93133 94552 110 84 105 95536 132190 134 85 104 225920 128754 54 86 108 62133 66363 96 87 98 61370 67808 78 88 69 43836 61724 51 89 111 106117 131722 121 90 99 38692 68580 38 91 71 84651 106175 145 92 27 56622 55792 59 93 69 15986 25157 27 94 107 95364 76669 91 95 73 26706 57283 48 96 107 89691 105805 68 97 93 67267 129484 58 98 129 126846 72413 150 99 69 41140 87831 74 100 118 102860 96971 181 101 73 51715 71299 65 102 119 55801 77494 97 103 104 111813 120336 121 104 107 120293 93913 99 105 99 138599 136048 152 106 90 161647 181248 188 107 197 115929 146123 138 108 36 24266 32036 40 109 85 162901 186646 254 110 139 109825 102255 87 111 106 129838 168237 178 112 50 37510 64219 51 113 64 43750 19630 49 114 31 40652 76825 73 115 63 87771 115338 176 116 92 85872 109427 94 117 106 89275 118168 120 118 63 44418 84845 66 119 69 192565 153197 56 120 41 35232 29877 39 121 56 40909 63506 66 122 25 13294 22445 27 123 65 32387 47695 65 124 93 140867 68370 58 125 114 120662 146304 98 126 38 21233 38233 25 127 44 44332 42071 26 128 87 61056 50517 77 129 110 101338 103950 130 130 0 1168 5841 11 131 27 13497 2341 2 132 83 65567 84396 101 133 30 25162 24610 31 134 80 32334 35753 36 135 98 40735 55515 120 136 82 91413 209056 195 137 0 855 6622 4 138 60 97068 115814 89 139 28 44339 11609 24 140 9 14116 13155 39 141 33 10288 18274 14 142 59 65622 72875 78 143 49 16563 10112 15 144 115 76643 142775 106 145 140 110681 68847 83 146 49 29011 17659 24 147 120 92696 20112 37 148 66 94785 61023 77 149 21 8773 13983 16 150 124 83209 65176 56 151 152 93815 132432 132 152 139 86687 112494 144 153 38 34553 45109 40 154 144 105547 170875 153 155 120 103487 180759 143 156 160 213688 214921 220 157 114 71220 100226 79 158 39 23517 32043 50 159 78 56926 54454 39 160 119 91721 78876 95 161 141 115168 170745 169 162 101 111194 6940 12 163 56 51009 49025 63 164 133 135777 122037 134 165 83 51513 53782 69 166 116 74163 127748 119 167 90 51633 86839 119 168 36 75345 44830 75 169 50 33416 77395 63 170 61 83305 89324 55 171 97 98952 103300 103 172 98 102372 112283 197 173 78 37238 10901 16 174 117 103772 120691 140 175 148 123969 58106 89 176 41 27142 57140 40 177 105 135400 122422 125 178 55 21399 25899 21 179 132 130115 139296 167 180 44 24874 52678 32 181 21 34988 23853 36 182 50 45549 17306 13 183 0 6023 7953 5 184 73 64466 89455 96 185 86 54990 147866 151 186 0 1644 4245 6 187 13 6179 21509 13 188 4 3926 7670 3 189 57 32755 66675 57 190 48 34777 14336 23 191 46 73224 53608 61 192 48 27114 30059 21 193 32 20760 29668 43 194 68 37636 22097 20 195 87 65461 96841 82 196 43 30080 41907 90 197 67 24094 27080 25 198 46 69008 35885 60 199 46 54968 41247 61 200 56 46090 28313 85 201 48 27507 36845 43 202 44 10672 16548 25 203 60 34029 36134 41 204 65 46300 55764 26 205 55 24760 28910 38 206 38 18779 13339 12 207 52 21280 25319 29 208 60 40662 66956 49 209 54 28987 47487 46 210 86 22827 52785 41 211 24 18513 44683 31 212 52 30594 35619 41 213 49 24006 21920 26 214 61 27913 45608 23 215 61 42744 7721 14 216 81 12934 20634 16 217 43 22574 29788 25 218 40 41385 31931 21 219 40 18653 37754 32 220 56 18472 32505 9 221 68 30976 40557 35 222 79 63339 94238 42 223 47 25568 44197 68 224 57 33747 43228 32 225 41 4154 4103 6 226 29 19474 44144 68 227 3 35130 32868 33 228 60 39067 27640 84 229 30 13310 14063 46 230 79 65892 28990 30 231 47 4143 4694 0 232 40 28579 42648 36 233 48 51776 64329 47 234 36 21152 21928 20 235 42 38084 25836 50 236 49 27717 22779 30 237 57 32928 40820 30 238 12 11342 27530 34 239 40 19499 32378 33 240 43 16380 10824 34 241 33 36874 39613 37 242 77 48259 60865 83 243 43 16734 19787 32 244 45 28207 20107 30 245 47 30143 36605 43 246 43 41369 40961 41 247 45 45833 48231 51 248 50 29156 39725 19 249 35 35944 21455 37 250 7 36278 23430 33 251 71 45588 62991 41 252 67 45097 49363 54 253 0 3895 9604 14 254 62 28394 24552 25 255 54 18632 31493 25 256 4 2325 3439 8 257 25 25139 19555 26 258 40 27975 21228 20 259 38 14483 23177 11 260 19 13127 22094 14 261 17 5839 2342 3 262 67 24069 38798 40 263 14 3738 3255 5 264 30 18625 24261 38 265 54 36341 18511 32 266 35 24548 40798 41 267 59 21792 28893 46 268 24 26263 21425 47 269 58 23686 50276 37 270 42 49303 37643 51 271 46 25659 30377 49 272 61 28904 27126 21 273 3 2781 13 1 274 52 29236 42097 44 275 25 19546 24451 26 276 40 22818 14335 21 277 32 32689 5084 4 278 4 5752 9927 10 279 49 22197 43527 43 280 63 20055 27184 34 281 67 25272 21610 32 282 32 82206 20484 20 283 23 32073 20156 34 284 7 5444 6012 6 285 54 20154 18475 12 286 37 36944 12645 24 287 35 8019 11017 16 288 51 30884 37623 72 289 39 19540 35873 27 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) time_in_rfc compendium_views_info 1.195e+02 1.264e-03 1.850e+00 compendium_views_pr blogged_computations compendiums_reviewed 1.612e+00 1.217e+00 -4.987e+00 feedback_messages_p1 feedback_messages_p120 totsize 1.864e+00 -1.999e-02 2.081e-04 totseconds tothyperlinks -5.580e-04 2.649e-01 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -594.68 -92.35 -17.83 75.95 1344.02 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.195e+02 2.814e+01 4.246 2.98e-05 *** time_in_rfc 1.264e-03 4.832e-04 2.616 0.00938 ** compendium_views_info 1.850e+00 9.560e-02 19.349 < 2e-16 *** compendium_views_pr 1.612e+00 2.513e-01 6.413 6.10e-10 *** blogged_computations 1.217e+00 7.204e-01 1.690 0.09224 . compendiums_reviewed -4.987e+00 5.726e+00 -0.871 0.38448 feedback_messages_p1 1.864e+00 1.741e+00 1.071 0.28516 feedback_messages_p120 -1.999e-02 8.446e-01 -0.024 0.98113 totsize 2.081e-04 5.125e-04 0.406 0.68500 totseconds -5.580e-04 6.821e-04 -0.818 0.41405 tothyperlinks 2.649e-01 5.575e-01 0.475 0.63509 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 179.1 on 278 degrees of freedom Multiple R-squared: 0.9429, Adjusted R-squared: 0.9408 F-statistic: 458.9 on 10 and 278 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.1104966 2.209931e-01 8.895034e-01 [2,] 0.1410375 2.820751e-01 8.589625e-01 [3,] 0.1273157 2.546314e-01 8.726843e-01 [4,] 0.3677743 7.355485e-01 6.322257e-01 [5,] 0.5599303 8.801394e-01 4.400697e-01 [6,] 0.5323272 9.353457e-01 4.676728e-01 [7,] 0.5487421 9.025158e-01 4.512579e-01 [8,] 0.5366051 9.267899e-01 4.633949e-01 [9,] 0.4510291 9.020581e-01 5.489709e-01 [10,] 0.4087136 8.174272e-01 5.912864e-01 [11,] 0.3277835 6.555670e-01 6.722165e-01 [12,] 0.3142443 6.284885e-01 6.857557e-01 [13,] 0.2585127 5.170253e-01 7.414873e-01 [14,] 0.2693458 5.386915e-01 7.306542e-01 [15,] 0.2830831 5.661662e-01 7.169169e-01 [16,] 0.2712327 5.424655e-01 7.287673e-01 [17,] 0.2200424 4.400849e-01 7.799576e-01 [18,] 0.3306789 6.613577e-01 6.693211e-01 [19,] 0.2998765 5.997529e-01 7.001235e-01 [20,] 0.2927649 5.855298e-01 7.072351e-01 [21,] 0.3109951 6.219902e-01 6.890049e-01 [22,] 0.2854918 5.709836e-01 7.145082e-01 [23,] 0.2459535 4.919070e-01 7.540465e-01 [24,] 0.2170617 4.341233e-01 7.829383e-01 [25,] 0.2154932 4.309864e-01 7.845068e-01 [26,] 0.2837172 5.674344e-01 7.162828e-01 [27,] 0.2396977 4.793955e-01 7.603023e-01 [28,] 0.2013531 4.027062e-01 7.986469e-01 [29,] 0.1815025 3.630051e-01 8.184975e-01 [30,] 0.3337543 6.675087e-01 6.662457e-01 [31,] 0.3851528 7.703056e-01 6.148472e-01 [32,] 0.3836516 7.673032e-01 6.163484e-01 [33,] 0.3927238 7.854477e-01 6.072762e-01 [34,] 0.5673488 8.653023e-01 4.326512e-01 [35,] 0.5220932 9.558135e-01 4.779068e-01 [36,] 0.4837602 9.675204e-01 5.162398e-01 [37,] 0.5091440 9.817119e-01 4.908560e-01 [38,] 0.7283221 5.433559e-01 2.716779e-01 [39,] 0.6932610 6.134780e-01 3.067390e-01 [40,] 0.6755117 6.489767e-01 3.244883e-01 [41,] 0.8808640 2.382721e-01 1.191360e-01 [42,] 0.9006628 1.986743e-01 9.933716e-02 [43,] 0.8857198 2.285605e-01 1.142802e-01 [44,] 0.8684847 2.630305e-01 1.315153e-01 [45,] 0.8639698 2.720603e-01 1.360302e-01 [46,] 0.8389281 3.221438e-01 1.610719e-01 [47,] 0.8138527 3.722946e-01 1.861473e-01 [48,] 0.7951524 4.096952e-01 2.048476e-01 [49,] 0.8811365 2.377269e-01 1.188635e-01 [50,] 1.0000000 1.986683e-10 9.933417e-11 [51,] 1.0000000 3.711154e-10 1.855577e-10 [52,] 1.0000000 1.926756e-10 9.633780e-11 [53,] 1.0000000 1.287199e-10 6.435997e-11 [54,] 1.0000000 2.161715e-10 1.080857e-10 [55,] 1.0000000 4.144731e-11 2.072366e-11 [56,] 1.0000000 7.056773e-11 3.528387e-11 [57,] 1.0000000 1.203706e-10 6.018531e-11 [58,] 1.0000000 2.025365e-10 1.012682e-10 [59,] 1.0000000 1.423049e-10 7.115245e-11 [60,] 1.0000000 1.699735e-10 8.498676e-11 [61,] 1.0000000 2.545604e-10 1.272802e-10 [62,] 1.0000000 6.371928e-11 3.185964e-11 [63,] 1.0000000 9.115108e-11 4.557554e-11 [64,] 1.0000000 1.359941e-10 6.799703e-11 [65,] 1.0000000 1.733734e-10 8.668669e-11 [66,] 1.0000000 2.240582e-10 1.120291e-10 [67,] 1.0000000 3.259493e-10 1.629747e-10 [68,] 1.0000000 9.354169e-11 4.677085e-11 [69,] 1.0000000 1.660450e-10 8.302248e-11 [70,] 1.0000000 1.030659e-11 5.153297e-12 [71,] 1.0000000 1.070380e-11 5.351899e-12 [72,] 1.0000000 1.428935e-11 7.144674e-12 [73,] 1.0000000 2.337813e-11 1.168906e-11 [74,] 1.0000000 3.523967e-11 1.761983e-11 [75,] 1.0000000 6.203278e-11 3.101639e-11 [76,] 1.0000000 5.404181e-11 2.702091e-11 [77,] 1.0000000 2.174891e-11 1.087446e-11 [78,] 1.0000000 1.105233e-11 5.526165e-12 [79,] 1.0000000 5.857353e-12 2.928676e-12 [80,] 1.0000000 2.127021e-12 1.063510e-12 [81,] 1.0000000 4.010938e-12 2.005469e-12 [82,] 1.0000000 6.136939e-12 3.068469e-12 [83,] 1.0000000 3.546278e-12 1.773139e-12 [84,] 1.0000000 6.439136e-12 3.219568e-12 [85,] 1.0000000 5.108403e-12 2.554202e-12 [86,] 1.0000000 7.267881e-12 3.633941e-12 [87,] 1.0000000 1.299870e-11 6.499350e-12 [88,] 1.0000000 1.858168e-11 9.290838e-12 [89,] 1.0000000 3.021845e-11 1.510923e-11 [90,] 1.0000000 1.852521e-11 9.262603e-12 [91,] 1.0000000 3.051216e-11 1.525608e-11 [92,] 1.0000000 4.723738e-11 2.361869e-11 [93,] 1.0000000 4.924260e-11 2.462130e-11 [94,] 1.0000000 1.000064e-11 5.000318e-12 [95,] 1.0000000 6.313658e-12 3.156829e-12 [96,] 1.0000000 1.022354e-11 5.111768e-12 [97,] 1.0000000 1.637676e-11 8.188378e-12 [98,] 1.0000000 2.205511e-11 1.102755e-11 [99,] 1.0000000 3.753229e-11 1.876614e-11 [100,] 1.0000000 6.082306e-11 3.041153e-11 [101,] 1.0000000 9.008190e-11 4.504095e-11 [102,] 1.0000000 1.499761e-10 7.498806e-11 [103,] 1.0000000 2.612226e-11 1.306113e-11 [104,] 1.0000000 4.039068e-11 2.019534e-11 [105,] 1.0000000 4.738130e-11 2.369065e-11 [106,] 1.0000000 8.478535e-11 4.239267e-11 [107,] 1.0000000 8.862275e-11 4.431138e-11 [108,] 1.0000000 1.122758e-10 5.613790e-11 [109,] 1.0000000 1.745518e-10 8.727591e-11 [110,] 1.0000000 2.474178e-10 1.237089e-10 [111,] 1.0000000 4.010473e-10 2.005236e-10 [112,] 1.0000000 6.849792e-10 3.424896e-10 [113,] 1.0000000 1.065453e-09 5.327264e-10 [114,] 1.0000000 1.807220e-09 9.036098e-10 [115,] 1.0000000 1.712446e-09 8.562231e-10 [116,] 1.0000000 2.656484e-09 1.328242e-09 [117,] 1.0000000 4.424096e-09 2.212048e-09 [118,] 1.0000000 6.686089e-09 3.343045e-09 [119,] 1.0000000 1.108285e-08 5.541425e-09 [120,] 1.0000000 1.756071e-08 8.780354e-09 [121,] 1.0000000 7.111565e-09 3.555783e-09 [122,] 1.0000000 1.185909e-08 5.929544e-09 [123,] 1.0000000 1.256243e-08 6.281215e-09 [124,] 1.0000000 2.085094e-08 1.042547e-08 [125,] 1.0000000 2.408126e-08 1.204063e-08 [126,] 1.0000000 3.822247e-08 1.911123e-08 [127,] 1.0000000 4.637215e-08 2.318607e-08 [128,] 1.0000000 7.503511e-08 3.751756e-08 [129,] 1.0000000 5.342098e-08 2.671049e-08 [130,] 1.0000000 5.961257e-08 2.980629e-08 [131,] 1.0000000 8.443720e-08 4.221860e-08 [132,] 0.9999999 1.338164e-07 6.690818e-08 [133,] 0.9999999 1.195813e-07 5.979066e-08 [134,] 0.9999999 1.865187e-07 9.325936e-08 [135,] 0.9999999 2.979297e-07 1.489649e-07 [136,] 0.9999999 2.228283e-07 1.114141e-07 [137,] 0.9999998 3.302268e-07 1.651134e-07 [138,] 0.9999998 3.585011e-07 1.792506e-07 [139,] 0.9999998 3.789317e-07 1.894659e-07 [140,] 0.9999998 3.240757e-07 1.620378e-07 [141,] 1.0000000 9.028324e-09 4.514162e-09 [142,] 1.0000000 6.589968e-09 3.294984e-09 [143,] 1.0000000 1.091961e-08 5.459804e-09 [144,] 1.0000000 5.357354e-09 2.678677e-09 [145,] 1.0000000 2.030655e-09 1.015328e-09 [146,] 1.0000000 3.181484e-09 1.590742e-09 [147,] 1.0000000 5.167728e-09 2.583864e-09 [148,] 1.0000000 3.653967e-09 1.826984e-09 [149,] 1.0000000 5.887299e-09 2.943650e-09 [150,] 1.0000000 8.052971e-09 4.026485e-09 [151,] 1.0000000 1.219992e-09 6.099959e-10 [152,] 1.0000000 2.148884e-09 1.074442e-09 [153,] 1.0000000 2.198327e-09 1.099164e-09 [154,] 1.0000000 1.531985e-09 7.659923e-10 [155,] 1.0000000 1.660050e-09 8.300251e-10 [156,] 1.0000000 2.518146e-09 1.259073e-09 [157,] 1.0000000 4.068150e-09 2.034075e-09 [158,] 1.0000000 4.097095e-09 2.048548e-09 [159,] 1.0000000 5.289029e-09 2.644514e-09 [160,] 1.0000000 2.052217e-09 1.026108e-09 [161,] 1.0000000 2.385416e-09 1.192708e-09 [162,] 1.0000000 6.055665e-10 3.027833e-10 [163,] 1.0000000 6.559333e-10 3.279666e-10 [164,] 1.0000000 7.095716e-10 3.547858e-10 [165,] 1.0000000 1.269009e-09 6.345047e-10 [166,] 1.0000000 1.337894e-09 6.689468e-10 [167,] 1.0000000 1.905008e-09 9.525040e-10 [168,] 1.0000000 2.171499e-09 1.085749e-09 [169,] 1.0000000 3.029891e-09 1.514945e-09 [170,] 1.0000000 5.284160e-09 2.642080e-09 [171,] 1.0000000 9.157148e-09 4.578574e-09 [172,] 1.0000000 7.845505e-09 3.922752e-09 [173,] 1.0000000 1.328919e-08 6.644595e-09 [174,] 1.0000000 1.455663e-08 7.278315e-09 [175,] 1.0000000 2.416396e-08 1.208198e-08 [176,] 1.0000000 2.321516e-08 1.160758e-08 [177,] 1.0000000 3.542272e-08 1.771136e-08 [178,] 1.0000000 5.909283e-08 2.954641e-08 [179,] 1.0000000 4.604267e-08 2.302134e-08 [180,] 1.0000000 7.528390e-08 3.764195e-08 [181,] 0.9999999 1.077262e-07 5.386312e-08 [182,] 1.0000000 7.351597e-08 3.675799e-08 [183,] 0.9999999 1.251707e-07 6.258534e-08 [184,] 0.9999999 1.543904e-07 7.719518e-08 [185,] 0.9999999 1.719451e-07 8.597253e-08 [186,] 0.9999999 2.062701e-07 1.031350e-07 [187,] 0.9999999 2.817200e-07 1.408600e-07 [188,] 0.9999998 4.529300e-07 2.264650e-07 [189,] 0.9999996 7.729316e-07 3.864658e-07 [190,] 0.9999994 1.215229e-06 6.076143e-07 [191,] 0.9999990 1.959033e-06 9.795167e-07 [192,] 0.9999988 2.362973e-06 1.181487e-06 [193,] 0.9999980 3.947090e-06 1.973545e-06 [194,] 0.9999985 2.900734e-06 1.450367e-06 [195,] 0.9999989 2.289176e-06 1.144588e-06 [196,] 0.9999981 3.743496e-06 1.871748e-06 [197,] 0.9999974 5.190323e-06 2.595162e-06 [198,] 0.9999963 7.405645e-06 3.702822e-06 [199,] 0.9999975 5.089018e-06 2.544509e-06 [200,] 0.9999961 7.866169e-06 3.933084e-06 [201,] 0.9999941 1.181329e-05 5.906643e-06 [202,] 0.9999907 1.858523e-05 9.292613e-06 [203,] 0.9999851 2.989558e-05 1.494779e-05 [204,] 0.9999764 4.715450e-05 2.357725e-05 [205,] 0.9999693 6.147584e-05 3.073792e-05 [206,] 0.9999556 8.870611e-05 4.435305e-05 [207,] 0.9999335 1.330751e-04 6.653753e-05 [208,] 0.9998953 2.093456e-04 1.046728e-04 [209,] 0.9999855 2.897179e-05 1.448589e-05 [210,] 0.9999870 2.595819e-05 1.297910e-05 [211,] 0.9999784 4.326535e-05 2.163267e-05 [212,] 0.9999704 5.929837e-05 2.964918e-05 [213,] 0.9999565 8.694929e-05 4.347464e-05 [214,] 0.9999324 1.352693e-04 6.763463e-05 [215,] 0.9999450 1.099235e-04 5.496175e-05 [216,] 0.9999519 9.622257e-05 4.811129e-05 [217,] 0.9999755 4.907962e-05 2.453981e-05 [218,] 0.9999691 6.180704e-05 3.090352e-05 [219,] 0.9999765 4.696233e-05 2.348116e-05 [220,] 0.9999786 4.274539e-05 2.137269e-05 [221,] 0.9999627 7.459777e-05 3.729889e-05 [222,] 0.9999631 7.381096e-05 3.690548e-05 [223,] 0.9999373 1.254564e-04 6.272819e-05 [224,] 0.9999283 1.433651e-04 7.168255e-05 [225,] 0.9999043 1.913065e-04 9.565326e-05 [226,] 0.9998537 2.926655e-04 1.463327e-04 [227,] 0.9998025 3.950552e-04 1.975276e-04 [228,] 0.9996692 6.616341e-04 3.308171e-04 [229,] 0.9999971 5.770292e-06 2.885146e-06 [230,] 0.9999942 1.160960e-05 5.804802e-06 [231,] 0.9999891 2.179439e-05 1.089720e-05 [232,] 0.9999796 4.078485e-05 2.039242e-05 [233,] 0.9999666 6.684154e-05 3.342077e-05 [234,] 0.9999963 7.319417e-06 3.659708e-06 [235,] 0.9999924 1.521345e-05 7.606723e-06 [236,] 0.9999869 2.614594e-05 1.307297e-05 [237,] 0.9999770 4.595719e-05 2.297859e-05 [238,] 0.9999656 6.888021e-05 3.444011e-05 [239,] 0.9999316 1.367969e-04 6.839843e-05 [240,] 0.9998681 2.637703e-04 1.318851e-04 [241,] 0.9997569 4.861051e-04 2.430526e-04 [242,] 0.9995322 9.356907e-04 4.678454e-04 [243,] 0.9992997 1.400614e-03 7.003072e-04 [244,] 0.9986774 2.645264e-03 1.322632e-03 [245,] 0.9984062 3.187678e-03 1.593839e-03 [246,] 0.9970894 5.821222e-03 2.910611e-03 [247,] 0.9948766 1.024683e-02 5.123417e-03 [248,] 0.9969234 6.153223e-03 3.076611e-03 [249,] 0.9954139 9.172215e-03 4.586107e-03 [250,] 0.9923648 1.527045e-02 7.635227e-03 [251,] 0.9875677 2.486452e-02 1.243226e-02 [252,] 0.9780946 4.381084e-02 2.190542e-02 [253,] 0.9652838 6.943238e-02 3.471619e-02 [254,] 0.9462546 1.074909e-01 5.374544e-02 [255,] 0.9656826 6.863490e-02 3.431745e-02 [256,] 0.9508718 9.825632e-02 4.912816e-02 [257,] 0.9831372 3.372551e-02 1.686275e-02 [258,] 0.9660882 6.782360e-02 3.391180e-02 [259,] 0.9319767 1.360466e-01 6.802330e-02 [260,] 0.8891883 2.216234e-01 1.108117e-01 [261,] 0.8012676 3.974648e-01 1.987324e-01 [262,] 0.6540443 6.919115e-01 3.459557e-01 > postscript(file="/var/fisher/rcomp/tmp/150q01355434873.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/fisher/rcomp/tmp/2bwx51355434873.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/fisher/rcomp/tmp/3o3vp1355434873.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/fisher/rcomp/tmp/4nf2f1355434873.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/fisher/rcomp/tmp/5kl4c1355434873.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 29.89680087 -163.59692548 -76.25268655 -363.10568715 6.28096365 6 7 8 9 10 -105.39449048 -551.79617950 -101.12593542 66.20578365 -144.96806083 11 12 13 14 15 -43.63991998 -118.05314830 31.21731667 -51.28301505 -221.80806034 16 17 18 19 20 38.56372192 415.76548871 -87.06096291 -124.05592383 -170.74778310 21 22 23 24 25 164.61426552 40.16418142 -353.89403103 1.27380576 -37.73155458 26 27 28 29 30 -47.20884833 85.83032764 29.25439589 37.67056955 -92.88673947 31 32 33 34 35 -176.06937912 -124.13642303 4.44330430 -141.40785122 -71.27927800 36 37 38 39 40 -141.18157533 29.46387318 -207.22088665 168.61904164 -72.87228618 41 42 43 44 45 16.52662332 38.37566722 347.25072540 136.83149694 -277.67037259 46 47 48 49 50 -207.64239921 388.60083516 -119.97281924 93.42932904 73.78423823 51 52 53 54 55 -594.68330623 -30.26776893 153.09780218 352.14406014 -213.41915731 56 57 58 59 60 75.94595716 -96.25686273 45.57958136 -21.03338795 -92.01013079 61 62 63 64 65 37.05701620 361.61699692 1344.01533560 3.19419509 322.50788534 66 67 68 69 70 -219.15798211 -85.14415187 369.84845906 1.13112897 -82.33865044 71 72 73 74 75 68.32089590 -238.92124371 143.70438508 115.67585584 320.19651842 76 77 78 79 80 -114.45165761 -119.75341421 122.30761104 -167.58830894 95.85386093 81 82 83 84 85 -335.15826960 -61.01667857 -446.55060581 -181.06369354 145.73515580 86 87 88 89 90 55.34478929 -118.83876352 28.03993600 185.33387497 -303.66557525 91 92 93 94 95 239.50475926 -277.80566380 -269.32680846 -14.67504960 71.88188087 96 97 98 99 100 -158.09954777 20.65591723 198.50235739 95.79609323 -2.67603829 101 102 103 104 105 98.30836672 62.85447118 -214.80666435 -80.21857559 -92.35235205 106 107 108 109 110 -168.01513443 322.43573848 215.87065084 -86.98525529 60.91661051 111 112 113 114 115 -128.33400893 44.55575911 -100.64422564 94.56657781 -69.82907444 116 117 118 119 120 -339.99754357 -100.89651931 161.68030524 15.13390985 151.38039564 121 122 123 124 125 135.39585461 53.89541762 83.65278006 -4.61720773 -34.64379998 126 127 128 129 130 57.81855603 11.55571969 187.83115236 59.80751362 -66.36453109 131 132 133 134 135 73.95882290 -37.45162666 -95.03733133 307.39829755 14.69171425 136 137 138 139 140 -171.73710929 -34.79860527 -128.16782123 29.58559852 -149.97110263 141 142 143 144 145 14.09463919 239.61730198 -109.98177680 106.93380452 -0.48353932 146 147 148 149 150 195.88856867 -39.07247206 -21.89171677 -221.88780808 -81.62438678 151 152 153 154 155 -172.80696037 -169.70003193 178.73097238 -437.38503520 -134.74980470 156 157 158 159 160 73.79214278 -238.21965337 -200.88441598 -71.64254670 96.79039968 161 162 163 164 165 282.59922845 -63.76403276 111.03018983 -268.01760679 -10.73902051 166 167 168 169 170 -112.39500895 -98.90997586 -57.05457772 101.44396203 119.40457097 171 172 173 174 175 0.45067722 -55.53247389 215.63869598 188.25417871 -180.17957373 176 177 178 179 180 -113.55204823 5.40874321 -24.84993763 56.71478374 -51.97627456 181 182 183 184 185 164.67344274 186.09797615 -91.29190464 75.14522380 105.57721018 186 187 188 189 190 -122.55109946 -204.93856765 -132.85952400 206.14304209 -19.53770827 191 192 193 194 195 42.06237228 -187.94692948 107.28202769 151.44238765 -56.34364369 196 197 198 199 200 76.52117269 -122.36230799 -100.16582292 -21.53917262 -48.10930633 201 202 203 204 205 65.47896641 -46.27031455 -24.14332043 57.93600607 -15.44406481 206 207 208 209 210 -26.45687240 184.12894892 -98.42745231 13.57249013 -39.37595079 211 212 213 214 215 120.31477434 206.00465263 -55.23902181 -59.17211053 81.37862810 216 217 218 219 220 22.67718765 26.00939770 -39.86338357 -48.01458086 -87.96114960 221 222 223 224 225 52.21565858 488.82174542 312.37712708 28.83450586 -133.52494238 226 227 228 229 230 168.56038488 -19.08587316 -67.04258633 173.33742243 281.22518499 231 232 233 234 235 40.95688555 191.51941196 149.58425832 -12.88109822 -49.77304663 236 237 238 239 240 26.99750720 -29.19814137 -5.08975899 -57.28693376 -60.36432890 241 242 243 244 245 77.16212825 521.90875213 58.97492157 -54.08015067 -58.00343153 246 247 248 249 250 -27.25677854 283.77516860 -72.94769495 -76.93131152 -90.63952589 251 252 253 254 255 45.27324730 -23.22625175 -83.78021211 27.60167212 -63.36462173 256 257 258 259 260 -87.19796946 -32.09382940 -74.37228689 16.23829297 21.78533103 261 262 263 264 265 -49.27826544 13.10222606 51.48661329 28.27204725 -36.16903697 266 267 268 269 270 -14.91430371 -17.82581426 -98.78961121 -141.16023858 -43.04548120 271 272 273 274 275 196.27796433 -93.47114167 -47.92774506 40.28219379 45.18093852 276 277 278 279 280 118.44276851 58.23228362 -44.16439501 112.75136398 -103.82826083 281 282 283 284 285 -72.43252949 -30.22288612 -0.05604489 -50.78493517 -98.97742086 286 287 288 289 346.08184338 95.96999499 83.82039405 -47.30008121 > postscript(file="/var/fisher/rcomp/tmp/6r8381355434873.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 29.89680087 NA 1 -163.59692548 29.89680087 2 -76.25268655 -163.59692548 3 -363.10568715 -76.25268655 4 6.28096365 -363.10568715 5 -105.39449048 6.28096365 6 -551.79617950 -105.39449048 7 -101.12593542 -551.79617950 8 66.20578365 -101.12593542 9 -144.96806083 66.20578365 10 -43.63991998 -144.96806083 11 -118.05314830 -43.63991998 12 31.21731667 -118.05314830 13 -51.28301505 31.21731667 14 -221.80806034 -51.28301505 15 38.56372192 -221.80806034 16 415.76548871 38.56372192 17 -87.06096291 415.76548871 18 -124.05592383 -87.06096291 19 -170.74778310 -124.05592383 20 164.61426552 -170.74778310 21 40.16418142 164.61426552 22 -353.89403103 40.16418142 23 1.27380576 -353.89403103 24 -37.73155458 1.27380576 25 -47.20884833 -37.73155458 26 85.83032764 -47.20884833 27 29.25439589 85.83032764 28 37.67056955 29.25439589 29 -92.88673947 37.67056955 30 -176.06937912 -92.88673947 31 -124.13642303 -176.06937912 32 4.44330430 -124.13642303 33 -141.40785122 4.44330430 34 -71.27927800 -141.40785122 35 -141.18157533 -71.27927800 36 29.46387318 -141.18157533 37 -207.22088665 29.46387318 38 168.61904164 -207.22088665 39 -72.87228618 168.61904164 40 16.52662332 -72.87228618 41 38.37566722 16.52662332 42 347.25072540 38.37566722 43 136.83149694 347.25072540 44 -277.67037259 136.83149694 45 -207.64239921 -277.67037259 46 388.60083516 -207.64239921 47 -119.97281924 388.60083516 48 93.42932904 -119.97281924 49 73.78423823 93.42932904 50 -594.68330623 73.78423823 51 -30.26776893 -594.68330623 52 153.09780218 -30.26776893 53 352.14406014 153.09780218 54 -213.41915731 352.14406014 55 75.94595716 -213.41915731 56 -96.25686273 75.94595716 57 45.57958136 -96.25686273 58 -21.03338795 45.57958136 59 -92.01013079 -21.03338795 60 37.05701620 -92.01013079 61 361.61699692 37.05701620 62 1344.01533560 361.61699692 63 3.19419509 1344.01533560 64 322.50788534 3.19419509 65 -219.15798211 322.50788534 66 -85.14415187 -219.15798211 67 369.84845906 -85.14415187 68 1.13112897 369.84845906 69 -82.33865044 1.13112897 70 68.32089590 -82.33865044 71 -238.92124371 68.32089590 72 143.70438508 -238.92124371 73 115.67585584 143.70438508 74 320.19651842 115.67585584 75 -114.45165761 320.19651842 76 -119.75341421 -114.45165761 77 122.30761104 -119.75341421 78 -167.58830894 122.30761104 79 95.85386093 -167.58830894 80 -335.15826960 95.85386093 81 -61.01667857 -335.15826960 82 -446.55060581 -61.01667857 83 -181.06369354 -446.55060581 84 145.73515580 -181.06369354 85 55.34478929 145.73515580 86 -118.83876352 55.34478929 87 28.03993600 -118.83876352 88 185.33387497 28.03993600 89 -303.66557525 185.33387497 90 239.50475926 -303.66557525 91 -277.80566380 239.50475926 92 -269.32680846 -277.80566380 93 -14.67504960 -269.32680846 94 71.88188087 -14.67504960 95 -158.09954777 71.88188087 96 20.65591723 -158.09954777 97 198.50235739 20.65591723 98 95.79609323 198.50235739 99 -2.67603829 95.79609323 100 98.30836672 -2.67603829 101 62.85447118 98.30836672 102 -214.80666435 62.85447118 103 -80.21857559 -214.80666435 104 -92.35235205 -80.21857559 105 -168.01513443 -92.35235205 106 322.43573848 -168.01513443 107 215.87065084 322.43573848 108 -86.98525529 215.87065084 109 60.91661051 -86.98525529 110 -128.33400893 60.91661051 111 44.55575911 -128.33400893 112 -100.64422564 44.55575911 113 94.56657781 -100.64422564 114 -69.82907444 94.56657781 115 -339.99754357 -69.82907444 116 -100.89651931 -339.99754357 117 161.68030524 -100.89651931 118 15.13390985 161.68030524 119 151.38039564 15.13390985 120 135.39585461 151.38039564 121 53.89541762 135.39585461 122 83.65278006 53.89541762 123 -4.61720773 83.65278006 124 -34.64379998 -4.61720773 125 57.81855603 -34.64379998 126 11.55571969 57.81855603 127 187.83115236 11.55571969 128 59.80751362 187.83115236 129 -66.36453109 59.80751362 130 73.95882290 -66.36453109 131 -37.45162666 73.95882290 132 -95.03733133 -37.45162666 133 307.39829755 -95.03733133 134 14.69171425 307.39829755 135 -171.73710929 14.69171425 136 -34.79860527 -171.73710929 137 -128.16782123 -34.79860527 138 29.58559852 -128.16782123 139 -149.97110263 29.58559852 140 14.09463919 -149.97110263 141 239.61730198 14.09463919 142 -109.98177680 239.61730198 143 106.93380452 -109.98177680 144 -0.48353932 106.93380452 145 195.88856867 -0.48353932 146 -39.07247206 195.88856867 147 -21.89171677 -39.07247206 148 -221.88780808 -21.89171677 149 -81.62438678 -221.88780808 150 -172.80696037 -81.62438678 151 -169.70003193 -172.80696037 152 178.73097238 -169.70003193 153 -437.38503520 178.73097238 154 -134.74980470 -437.38503520 155 73.79214278 -134.74980470 156 -238.21965337 73.79214278 157 -200.88441598 -238.21965337 158 -71.64254670 -200.88441598 159 96.79039968 -71.64254670 160 282.59922845 96.79039968 161 -63.76403276 282.59922845 162 111.03018983 -63.76403276 163 -268.01760679 111.03018983 164 -10.73902051 -268.01760679 165 -112.39500895 -10.73902051 166 -98.90997586 -112.39500895 167 -57.05457772 -98.90997586 168 101.44396203 -57.05457772 169 119.40457097 101.44396203 170 0.45067722 119.40457097 171 -55.53247389 0.45067722 172 215.63869598 -55.53247389 173 188.25417871 215.63869598 174 -180.17957373 188.25417871 175 -113.55204823 -180.17957373 176 5.40874321 -113.55204823 177 -24.84993763 5.40874321 178 56.71478374 -24.84993763 179 -51.97627456 56.71478374 180 164.67344274 -51.97627456 181 186.09797615 164.67344274 182 -91.29190464 186.09797615 183 75.14522380 -91.29190464 184 105.57721018 75.14522380 185 -122.55109946 105.57721018 186 -204.93856765 -122.55109946 187 -132.85952400 -204.93856765 188 206.14304209 -132.85952400 189 -19.53770827 206.14304209 190 42.06237228 -19.53770827 191 -187.94692948 42.06237228 192 107.28202769 -187.94692948 193 151.44238765 107.28202769 194 -56.34364369 151.44238765 195 76.52117269 -56.34364369 196 -122.36230799 76.52117269 197 -100.16582292 -122.36230799 198 -21.53917262 -100.16582292 199 -48.10930633 -21.53917262 200 65.47896641 -48.10930633 201 -46.27031455 65.47896641 202 -24.14332043 -46.27031455 203 57.93600607 -24.14332043 204 -15.44406481 57.93600607 205 -26.45687240 -15.44406481 206 184.12894892 -26.45687240 207 -98.42745231 184.12894892 208 13.57249013 -98.42745231 209 -39.37595079 13.57249013 210 120.31477434 -39.37595079 211 206.00465263 120.31477434 212 -55.23902181 206.00465263 213 -59.17211053 -55.23902181 214 81.37862810 -59.17211053 215 22.67718765 81.37862810 216 26.00939770 22.67718765 217 -39.86338357 26.00939770 218 -48.01458086 -39.86338357 219 -87.96114960 -48.01458086 220 52.21565858 -87.96114960 221 488.82174542 52.21565858 222 312.37712708 488.82174542 223 28.83450586 312.37712708 224 -133.52494238 28.83450586 225 168.56038488 -133.52494238 226 -19.08587316 168.56038488 227 -67.04258633 -19.08587316 228 173.33742243 -67.04258633 229 281.22518499 173.33742243 230 40.95688555 281.22518499 231 191.51941196 40.95688555 232 149.58425832 191.51941196 233 -12.88109822 149.58425832 234 -49.77304663 -12.88109822 235 26.99750720 -49.77304663 236 -29.19814137 26.99750720 237 -5.08975899 -29.19814137 238 -57.28693376 -5.08975899 239 -60.36432890 -57.28693376 240 77.16212825 -60.36432890 241 521.90875213 77.16212825 242 58.97492157 521.90875213 243 -54.08015067 58.97492157 244 -58.00343153 -54.08015067 245 -27.25677854 -58.00343153 246 283.77516860 -27.25677854 247 -72.94769495 283.77516860 248 -76.93131152 -72.94769495 249 -90.63952589 -76.93131152 250 45.27324730 -90.63952589 251 -23.22625175 45.27324730 252 -83.78021211 -23.22625175 253 27.60167212 -83.78021211 254 -63.36462173 27.60167212 255 -87.19796946 -63.36462173 256 -32.09382940 -87.19796946 257 -74.37228689 -32.09382940 258 16.23829297 -74.37228689 259 21.78533103 16.23829297 260 -49.27826544 21.78533103 261 13.10222606 -49.27826544 262 51.48661329 13.10222606 263 28.27204725 51.48661329 264 -36.16903697 28.27204725 265 -14.91430371 -36.16903697 266 -17.82581426 -14.91430371 267 -98.78961121 -17.82581426 268 -141.16023858 -98.78961121 269 -43.04548120 -141.16023858 270 196.27796433 -43.04548120 271 -93.47114167 196.27796433 272 -47.92774506 -93.47114167 273 40.28219379 -47.92774506 274 45.18093852 40.28219379 275 118.44276851 45.18093852 276 58.23228362 118.44276851 277 -44.16439501 58.23228362 278 112.75136398 -44.16439501 279 -103.82826083 112.75136398 280 -72.43252949 -103.82826083 281 -30.22288612 -72.43252949 282 -0.05604489 -30.22288612 283 -50.78493517 -0.05604489 284 -98.97742086 -50.78493517 285 346.08184338 -98.97742086 286 95.96999499 346.08184338 287 83.82039405 95.96999499 288 -47.30008121 83.82039405 289 NA -47.30008121 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -163.59692548 29.89680087 [2,] -76.25268655 -163.59692548 [3,] -363.10568715 -76.25268655 [4,] 6.28096365 -363.10568715 [5,] -105.39449048 6.28096365 [6,] -551.79617950 -105.39449048 [7,] -101.12593542 -551.79617950 [8,] 66.20578365 -101.12593542 [9,] -144.96806083 66.20578365 [10,] -43.63991998 -144.96806083 [11,] -118.05314830 -43.63991998 [12,] 31.21731667 -118.05314830 [13,] -51.28301505 31.21731667 [14,] -221.80806034 -51.28301505 [15,] 38.56372192 -221.80806034 [16,] 415.76548871 38.56372192 [17,] -87.06096291 415.76548871 [18,] -124.05592383 -87.06096291 [19,] -170.74778310 -124.05592383 [20,] 164.61426552 -170.74778310 [21,] 40.16418142 164.61426552 [22,] -353.89403103 40.16418142 [23,] 1.27380576 -353.89403103 [24,] -37.73155458 1.27380576 [25,] -47.20884833 -37.73155458 [26,] 85.83032764 -47.20884833 [27,] 29.25439589 85.83032764 [28,] 37.67056955 29.25439589 [29,] -92.88673947 37.67056955 [30,] -176.06937912 -92.88673947 [31,] -124.13642303 -176.06937912 [32,] 4.44330430 -124.13642303 [33,] -141.40785122 4.44330430 [34,] -71.27927800 -141.40785122 [35,] -141.18157533 -71.27927800 [36,] 29.46387318 -141.18157533 [37,] -207.22088665 29.46387318 [38,] 168.61904164 -207.22088665 [39,] -72.87228618 168.61904164 [40,] 16.52662332 -72.87228618 [41,] 38.37566722 16.52662332 [42,] 347.25072540 38.37566722 [43,] 136.83149694 347.25072540 [44,] -277.67037259 136.83149694 [45,] -207.64239921 -277.67037259 [46,] 388.60083516 -207.64239921 [47,] -119.97281924 388.60083516 [48,] 93.42932904 -119.97281924 [49,] 73.78423823 93.42932904 [50,] -594.68330623 73.78423823 [51,] -30.26776893 -594.68330623 [52,] 153.09780218 -30.26776893 [53,] 352.14406014 153.09780218 [54,] -213.41915731 352.14406014 [55,] 75.94595716 -213.41915731 [56,] -96.25686273 75.94595716 [57,] 45.57958136 -96.25686273 [58,] -21.03338795 45.57958136 [59,] -92.01013079 -21.03338795 [60,] 37.05701620 -92.01013079 [61,] 361.61699692 37.05701620 [62,] 1344.01533560 361.61699692 [63,] 3.19419509 1344.01533560 [64,] 322.50788534 3.19419509 [65,] -219.15798211 322.50788534 [66,] -85.14415187 -219.15798211 [67,] 369.84845906 -85.14415187 [68,] 1.13112897 369.84845906 [69,] -82.33865044 1.13112897 [70,] 68.32089590 -82.33865044 [71,] -238.92124371 68.32089590 [72,] 143.70438508 -238.92124371 [73,] 115.67585584 143.70438508 [74,] 320.19651842 115.67585584 [75,] -114.45165761 320.19651842 [76,] -119.75341421 -114.45165761 [77,] 122.30761104 -119.75341421 [78,] -167.58830894 122.30761104 [79,] 95.85386093 -167.58830894 [80,] -335.15826960 95.85386093 [81,] -61.01667857 -335.15826960 [82,] -446.55060581 -61.01667857 [83,] -181.06369354 -446.55060581 [84,] 145.73515580 -181.06369354 [85,] 55.34478929 145.73515580 [86,] -118.83876352 55.34478929 [87,] 28.03993600 -118.83876352 [88,] 185.33387497 28.03993600 [89,] -303.66557525 185.33387497 [90,] 239.50475926 -303.66557525 [91,] -277.80566380 239.50475926 [92,] -269.32680846 -277.80566380 [93,] -14.67504960 -269.32680846 [94,] 71.88188087 -14.67504960 [95,] -158.09954777 71.88188087 [96,] 20.65591723 -158.09954777 [97,] 198.50235739 20.65591723 [98,] 95.79609323 198.50235739 [99,] -2.67603829 95.79609323 [100,] 98.30836672 -2.67603829 [101,] 62.85447118 98.30836672 [102,] -214.80666435 62.85447118 [103,] -80.21857559 -214.80666435 [104,] -92.35235205 -80.21857559 [105,] -168.01513443 -92.35235205 [106,] 322.43573848 -168.01513443 [107,] 215.87065084 322.43573848 [108,] -86.98525529 215.87065084 [109,] 60.91661051 -86.98525529 [110,] -128.33400893 60.91661051 [111,] 44.55575911 -128.33400893 [112,] -100.64422564 44.55575911 [113,] 94.56657781 -100.64422564 [114,] -69.82907444 94.56657781 [115,] -339.99754357 -69.82907444 [116,] -100.89651931 -339.99754357 [117,] 161.68030524 -100.89651931 [118,] 15.13390985 161.68030524 [119,] 151.38039564 15.13390985 [120,] 135.39585461 151.38039564 [121,] 53.89541762 135.39585461 [122,] 83.65278006 53.89541762 [123,] -4.61720773 83.65278006 [124,] -34.64379998 -4.61720773 [125,] 57.81855603 -34.64379998 [126,] 11.55571969 57.81855603 [127,] 187.83115236 11.55571969 [128,] 59.80751362 187.83115236 [129,] -66.36453109 59.80751362 [130,] 73.95882290 -66.36453109 [131,] -37.45162666 73.95882290 [132,] -95.03733133 -37.45162666 [133,] 307.39829755 -95.03733133 [134,] 14.69171425 307.39829755 [135,] -171.73710929 14.69171425 [136,] -34.79860527 -171.73710929 [137,] -128.16782123 -34.79860527 [138,] 29.58559852 -128.16782123 [139,] -149.97110263 29.58559852 [140,] 14.09463919 -149.97110263 [141,] 239.61730198 14.09463919 [142,] -109.98177680 239.61730198 [143,] 106.93380452 -109.98177680 [144,] -0.48353932 106.93380452 [145,] 195.88856867 -0.48353932 [146,] -39.07247206 195.88856867 [147,] -21.89171677 -39.07247206 [148,] -221.88780808 -21.89171677 [149,] -81.62438678 -221.88780808 [150,] -172.80696037 -81.62438678 [151,] -169.70003193 -172.80696037 [152,] 178.73097238 -169.70003193 [153,] -437.38503520 178.73097238 [154,] -134.74980470 -437.38503520 [155,] 73.79214278 -134.74980470 [156,] -238.21965337 73.79214278 [157,] -200.88441598 -238.21965337 [158,] -71.64254670 -200.88441598 [159,] 96.79039968 -71.64254670 [160,] 282.59922845 96.79039968 [161,] -63.76403276 282.59922845 [162,] 111.03018983 -63.76403276 [163,] -268.01760679 111.03018983 [164,] -10.73902051 -268.01760679 [165,] -112.39500895 -10.73902051 [166,] -98.90997586 -112.39500895 [167,] -57.05457772 -98.90997586 [168,] 101.44396203 -57.05457772 [169,] 119.40457097 101.44396203 [170,] 0.45067722 119.40457097 [171,] -55.53247389 0.45067722 [172,] 215.63869598 -55.53247389 [173,] 188.25417871 215.63869598 [174,] -180.17957373 188.25417871 [175,] -113.55204823 -180.17957373 [176,] 5.40874321 -113.55204823 [177,] -24.84993763 5.40874321 [178,] 56.71478374 -24.84993763 [179,] -51.97627456 56.71478374 [180,] 164.67344274 -51.97627456 [181,] 186.09797615 164.67344274 [182,] -91.29190464 186.09797615 [183,] 75.14522380 -91.29190464 [184,] 105.57721018 75.14522380 [185,] -122.55109946 105.57721018 [186,] -204.93856765 -122.55109946 [187,] -132.85952400 -204.93856765 [188,] 206.14304209 -132.85952400 [189,] -19.53770827 206.14304209 [190,] 42.06237228 -19.53770827 [191,] -187.94692948 42.06237228 [192,] 107.28202769 -187.94692948 [193,] 151.44238765 107.28202769 [194,] -56.34364369 151.44238765 [195,] 76.52117269 -56.34364369 [196,] -122.36230799 76.52117269 [197,] -100.16582292 -122.36230799 [198,] -21.53917262 -100.16582292 [199,] -48.10930633 -21.53917262 [200,] 65.47896641 -48.10930633 [201,] -46.27031455 65.47896641 [202,] -24.14332043 -46.27031455 [203,] 57.93600607 -24.14332043 [204,] -15.44406481 57.93600607 [205,] -26.45687240 -15.44406481 [206,] 184.12894892 -26.45687240 [207,] -98.42745231 184.12894892 [208,] 13.57249013 -98.42745231 [209,] -39.37595079 13.57249013 [210,] 120.31477434 -39.37595079 [211,] 206.00465263 120.31477434 [212,] -55.23902181 206.00465263 [213,] -59.17211053 -55.23902181 [214,] 81.37862810 -59.17211053 [215,] 22.67718765 81.37862810 [216,] 26.00939770 22.67718765 [217,] -39.86338357 26.00939770 [218,] -48.01458086 -39.86338357 [219,] -87.96114960 -48.01458086 [220,] 52.21565858 -87.96114960 [221,] 488.82174542 52.21565858 [222,] 312.37712708 488.82174542 [223,] 28.83450586 312.37712708 [224,] -133.52494238 28.83450586 [225,] 168.56038488 -133.52494238 [226,] -19.08587316 168.56038488 [227,] -67.04258633 -19.08587316 [228,] 173.33742243 -67.04258633 [229,] 281.22518499 173.33742243 [230,] 40.95688555 281.22518499 [231,] 191.51941196 40.95688555 [232,] 149.58425832 191.51941196 [233,] -12.88109822 149.58425832 [234,] -49.77304663 -12.88109822 [235,] 26.99750720 -49.77304663 [236,] -29.19814137 26.99750720 [237,] -5.08975899 -29.19814137 [238,] -57.28693376 -5.08975899 [239,] -60.36432890 -57.28693376 [240,] 77.16212825 -60.36432890 [241,] 521.90875213 77.16212825 [242,] 58.97492157 521.90875213 [243,] -54.08015067 58.97492157 [244,] -58.00343153 -54.08015067 [245,] -27.25677854 -58.00343153 [246,] 283.77516860 -27.25677854 [247,] -72.94769495 283.77516860 [248,] -76.93131152 -72.94769495 [249,] -90.63952589 -76.93131152 [250,] 45.27324730 -90.63952589 [251,] -23.22625175 45.27324730 [252,] -83.78021211 -23.22625175 [253,] 27.60167212 -83.78021211 [254,] -63.36462173 27.60167212 [255,] -87.19796946 -63.36462173 [256,] -32.09382940 -87.19796946 [257,] -74.37228689 -32.09382940 [258,] 16.23829297 -74.37228689 [259,] 21.78533103 16.23829297 [260,] -49.27826544 21.78533103 [261,] 13.10222606 -49.27826544 [262,] 51.48661329 13.10222606 [263,] 28.27204725 51.48661329 [264,] -36.16903697 28.27204725 [265,] -14.91430371 -36.16903697 [266,] -17.82581426 -14.91430371 [267,] -98.78961121 -17.82581426 [268,] -141.16023858 -98.78961121 [269,] -43.04548120 -141.16023858 [270,] 196.27796433 -43.04548120 [271,] -93.47114167 196.27796433 [272,] -47.92774506 -93.47114167 [273,] 40.28219379 -47.92774506 [274,] 45.18093852 40.28219379 [275,] 118.44276851 45.18093852 [276,] 58.23228362 118.44276851 [277,] -44.16439501 58.23228362 [278,] 112.75136398 -44.16439501 [279,] -103.82826083 112.75136398 [280,] -72.43252949 -103.82826083 [281,] -30.22288612 -72.43252949 [282,] -0.05604489 -30.22288612 [283,] -50.78493517 -0.05604489 [284,] -98.97742086 -50.78493517 [285,] 346.08184338 -98.97742086 [286,] 95.96999499 346.08184338 [287,] 83.82039405 95.96999499 [288,] -47.30008121 83.82039405 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -163.59692548 29.89680087 2 -76.25268655 -163.59692548 3 -363.10568715 -76.25268655 4 6.28096365 -363.10568715 5 -105.39449048 6.28096365 6 -551.79617950 -105.39449048 7 -101.12593542 -551.79617950 8 66.20578365 -101.12593542 9 -144.96806083 66.20578365 10 -43.63991998 -144.96806083 11 -118.05314830 -43.63991998 12 31.21731667 -118.05314830 13 -51.28301505 31.21731667 14 -221.80806034 -51.28301505 15 38.56372192 -221.80806034 16 415.76548871 38.56372192 17 -87.06096291 415.76548871 18 -124.05592383 -87.06096291 19 -170.74778310 -124.05592383 20 164.61426552 -170.74778310 21 40.16418142 164.61426552 22 -353.89403103 40.16418142 23 1.27380576 -353.89403103 24 -37.73155458 1.27380576 25 -47.20884833 -37.73155458 26 85.83032764 -47.20884833 27 29.25439589 85.83032764 28 37.67056955 29.25439589 29 -92.88673947 37.67056955 30 -176.06937912 -92.88673947 31 -124.13642303 -176.06937912 32 4.44330430 -124.13642303 33 -141.40785122 4.44330430 34 -71.27927800 -141.40785122 35 -141.18157533 -71.27927800 36 29.46387318 -141.18157533 37 -207.22088665 29.46387318 38 168.61904164 -207.22088665 39 -72.87228618 168.61904164 40 16.52662332 -72.87228618 41 38.37566722 16.52662332 42 347.25072540 38.37566722 43 136.83149694 347.25072540 44 -277.67037259 136.83149694 45 -207.64239921 -277.67037259 46 388.60083516 -207.64239921 47 -119.97281924 388.60083516 48 93.42932904 -119.97281924 49 73.78423823 93.42932904 50 -594.68330623 73.78423823 51 -30.26776893 -594.68330623 52 153.09780218 -30.26776893 53 352.14406014 153.09780218 54 -213.41915731 352.14406014 55 75.94595716 -213.41915731 56 -96.25686273 75.94595716 57 45.57958136 -96.25686273 58 -21.03338795 45.57958136 59 -92.01013079 -21.03338795 60 37.05701620 -92.01013079 61 361.61699692 37.05701620 62 1344.01533560 361.61699692 63 3.19419509 1344.01533560 64 322.50788534 3.19419509 65 -219.15798211 322.50788534 66 -85.14415187 -219.15798211 67 369.84845906 -85.14415187 68 1.13112897 369.84845906 69 -82.33865044 1.13112897 70 68.32089590 -82.33865044 71 -238.92124371 68.32089590 72 143.70438508 -238.92124371 73 115.67585584 143.70438508 74 320.19651842 115.67585584 75 -114.45165761 320.19651842 76 -119.75341421 -114.45165761 77 122.30761104 -119.75341421 78 -167.58830894 122.30761104 79 95.85386093 -167.58830894 80 -335.15826960 95.85386093 81 -61.01667857 -335.15826960 82 -446.55060581 -61.01667857 83 -181.06369354 -446.55060581 84 145.73515580 -181.06369354 85 55.34478929 145.73515580 86 -118.83876352 55.34478929 87 28.03993600 -118.83876352 88 185.33387497 28.03993600 89 -303.66557525 185.33387497 90 239.50475926 -303.66557525 91 -277.80566380 239.50475926 92 -269.32680846 -277.80566380 93 -14.67504960 -269.32680846 94 71.88188087 -14.67504960 95 -158.09954777 71.88188087 96 20.65591723 -158.09954777 97 198.50235739 20.65591723 98 95.79609323 198.50235739 99 -2.67603829 95.79609323 100 98.30836672 -2.67603829 101 62.85447118 98.30836672 102 -214.80666435 62.85447118 103 -80.21857559 -214.80666435 104 -92.35235205 -80.21857559 105 -168.01513443 -92.35235205 106 322.43573848 -168.01513443 107 215.87065084 322.43573848 108 -86.98525529 215.87065084 109 60.91661051 -86.98525529 110 -128.33400893 60.91661051 111 44.55575911 -128.33400893 112 -100.64422564 44.55575911 113 94.56657781 -100.64422564 114 -69.82907444 94.56657781 115 -339.99754357 -69.82907444 116 -100.89651931 -339.99754357 117 161.68030524 -100.89651931 118 15.13390985 161.68030524 119 151.38039564 15.13390985 120 135.39585461 151.38039564 121 53.89541762 135.39585461 122 83.65278006 53.89541762 123 -4.61720773 83.65278006 124 -34.64379998 -4.61720773 125 57.81855603 -34.64379998 126 11.55571969 57.81855603 127 187.83115236 11.55571969 128 59.80751362 187.83115236 129 -66.36453109 59.80751362 130 73.95882290 -66.36453109 131 -37.45162666 73.95882290 132 -95.03733133 -37.45162666 133 307.39829755 -95.03733133 134 14.69171425 307.39829755 135 -171.73710929 14.69171425 136 -34.79860527 -171.73710929 137 -128.16782123 -34.79860527 138 29.58559852 -128.16782123 139 -149.97110263 29.58559852 140 14.09463919 -149.97110263 141 239.61730198 14.09463919 142 -109.98177680 239.61730198 143 106.93380452 -109.98177680 144 -0.48353932 106.93380452 145 195.88856867 -0.48353932 146 -39.07247206 195.88856867 147 -21.89171677 -39.07247206 148 -221.88780808 -21.89171677 149 -81.62438678 -221.88780808 150 -172.80696037 -81.62438678 151 -169.70003193 -172.80696037 152 178.73097238 -169.70003193 153 -437.38503520 178.73097238 154 -134.74980470 -437.38503520 155 73.79214278 -134.74980470 156 -238.21965337 73.79214278 157 -200.88441598 -238.21965337 158 -71.64254670 -200.88441598 159 96.79039968 -71.64254670 160 282.59922845 96.79039968 161 -63.76403276 282.59922845 162 111.03018983 -63.76403276 163 -268.01760679 111.03018983 164 -10.73902051 -268.01760679 165 -112.39500895 -10.73902051 166 -98.90997586 -112.39500895 167 -57.05457772 -98.90997586 168 101.44396203 -57.05457772 169 119.40457097 101.44396203 170 0.45067722 119.40457097 171 -55.53247389 0.45067722 172 215.63869598 -55.53247389 173 188.25417871 215.63869598 174 -180.17957373 188.25417871 175 -113.55204823 -180.17957373 176 5.40874321 -113.55204823 177 -24.84993763 5.40874321 178 56.71478374 -24.84993763 179 -51.97627456 56.71478374 180 164.67344274 -51.97627456 181 186.09797615 164.67344274 182 -91.29190464 186.09797615 183 75.14522380 -91.29190464 184 105.57721018 75.14522380 185 -122.55109946 105.57721018 186 -204.93856765 -122.55109946 187 -132.85952400 -204.93856765 188 206.14304209 -132.85952400 189 -19.53770827 206.14304209 190 42.06237228 -19.53770827 191 -187.94692948 42.06237228 192 107.28202769 -187.94692948 193 151.44238765 107.28202769 194 -56.34364369 151.44238765 195 76.52117269 -56.34364369 196 -122.36230799 76.52117269 197 -100.16582292 -122.36230799 198 -21.53917262 -100.16582292 199 -48.10930633 -21.53917262 200 65.47896641 -48.10930633 201 -46.27031455 65.47896641 202 -24.14332043 -46.27031455 203 57.93600607 -24.14332043 204 -15.44406481 57.93600607 205 -26.45687240 -15.44406481 206 184.12894892 -26.45687240 207 -98.42745231 184.12894892 208 13.57249013 -98.42745231 209 -39.37595079 13.57249013 210 120.31477434 -39.37595079 211 206.00465263 120.31477434 212 -55.23902181 206.00465263 213 -59.17211053 -55.23902181 214 81.37862810 -59.17211053 215 22.67718765 81.37862810 216 26.00939770 22.67718765 217 -39.86338357 26.00939770 218 -48.01458086 -39.86338357 219 -87.96114960 -48.01458086 220 52.21565858 -87.96114960 221 488.82174542 52.21565858 222 312.37712708 488.82174542 223 28.83450586 312.37712708 224 -133.52494238 28.83450586 225 168.56038488 -133.52494238 226 -19.08587316 168.56038488 227 -67.04258633 -19.08587316 228 173.33742243 -67.04258633 229 281.22518499 173.33742243 230 40.95688555 281.22518499 231 191.51941196 40.95688555 232 149.58425832 191.51941196 233 -12.88109822 149.58425832 234 -49.77304663 -12.88109822 235 26.99750720 -49.77304663 236 -29.19814137 26.99750720 237 -5.08975899 -29.19814137 238 -57.28693376 -5.08975899 239 -60.36432890 -57.28693376 240 77.16212825 -60.36432890 241 521.90875213 77.16212825 242 58.97492157 521.90875213 243 -54.08015067 58.97492157 244 -58.00343153 -54.08015067 245 -27.25677854 -58.00343153 246 283.77516860 -27.25677854 247 -72.94769495 283.77516860 248 -76.93131152 -72.94769495 249 -90.63952589 -76.93131152 250 45.27324730 -90.63952589 251 -23.22625175 45.27324730 252 -83.78021211 -23.22625175 253 27.60167212 -83.78021211 254 -63.36462173 27.60167212 255 -87.19796946 -63.36462173 256 -32.09382940 -87.19796946 257 -74.37228689 -32.09382940 258 16.23829297 -74.37228689 259 21.78533103 16.23829297 260 -49.27826544 21.78533103 261 13.10222606 -49.27826544 262 51.48661329 13.10222606 263 28.27204725 51.48661329 264 -36.16903697 28.27204725 265 -14.91430371 -36.16903697 266 -17.82581426 -14.91430371 267 -98.78961121 -17.82581426 268 -141.16023858 -98.78961121 269 -43.04548120 -141.16023858 270 196.27796433 -43.04548120 271 -93.47114167 196.27796433 272 -47.92774506 -93.47114167 273 40.28219379 -47.92774506 274 45.18093852 40.28219379 275 118.44276851 45.18093852 276 58.23228362 118.44276851 277 -44.16439501 58.23228362 278 112.75136398 -44.16439501 279 -103.82826083 112.75136398 280 -72.43252949 -103.82826083 281 -30.22288612 -72.43252949 282 -0.05604489 -30.22288612 283 -50.78493517 -0.05604489 284 -98.97742086 -50.78493517 285 346.08184338 -98.97742086 286 95.96999499 346.08184338 287 83.82039405 95.96999499 288 -47.30008121 83.82039405 > 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/fisher/rcomp/tmp/7k4t81355434874.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/fisher/rcomp/tmp/8muz41355434874.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/fisher/rcomp/tmp/959kq1355434874.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/fisher/rcomp/tmp/109evb1355434874.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/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/fisher/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/fisher/rcomp/tmp/110p201355434874.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/fisher/rcomp/tmp/121qr81355434874.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/fisher/rcomp/tmp/1305ba1355434874.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/fisher/rcomp/tmp/144i861355434874.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/fisher/rcomp/tmp/15k1x91355434874.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/fisher/rcomp/tmp/16750l1355434875.tab") + } > > try(system("convert tmp/150q01355434873.ps tmp/150q01355434873.png",intern=TRUE)) character(0) > try(system("convert tmp/2bwx51355434873.ps tmp/2bwx51355434873.png",intern=TRUE)) character(0) > try(system("convert tmp/3o3vp1355434873.ps tmp/3o3vp1355434873.png",intern=TRUE)) character(0) > try(system("convert tmp/4nf2f1355434873.ps tmp/4nf2f1355434873.png",intern=TRUE)) character(0) > try(system("convert tmp/5kl4c1355434873.ps tmp/5kl4c1355434873.png",intern=TRUE)) character(0) > try(system("convert tmp/6r8381355434873.ps tmp/6r8381355434873.png",intern=TRUE)) character(0) > try(system("convert tmp/7k4t81355434874.ps tmp/7k4t81355434874.png",intern=TRUE)) character(0) > try(system("convert tmp/8muz41355434874.ps tmp/8muz41355434874.png",intern=TRUE)) character(0) > try(system("convert tmp/959kq1355434874.ps tmp/959kq1355434874.png",intern=TRUE)) character(0) > try(system("convert tmp/109evb1355434874.ps tmp/109evb1355434874.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.970 1.636 15.612