R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(210907 + ,79 + ,81 + ,94 + ,56 + ,3 + ,1418 + ,112285 + ,120982 + ,58 + ,55 + ,103 + ,56 + ,4 + ,869 + ,84786 + ,176508 + ,60 + ,50 + ,93 + ,54 + ,12 + ,1530 + ,83123 + ,179321 + ,108 + ,125 + ,103 + ,89 + ,2 + ,2172 + ,101193 + ,123185 + ,49 + ,40 + ,51 + ,40 + ,1 + ,901 + ,38361 + ,52746 + ,0 + ,37 + ,70 + ,25 + ,3 + ,463 + ,68504 + ,385534 + ,121 + ,63 + ,91 + ,92 + ,0 + ,3201 + ,119182 + ,33170 + ,1 + ,44 + ,22 + ,18 + ,0 + ,371 + ,22807 + ,101645 + ,20 + ,88 + ,38 + ,63 + ,0 + ,1192 + ,17140 + ,149061 + ,43 + ,66 + ,93 + ,44 + ,5 + ,1583 + ,116174 + ,165446 + ,69 + ,57 + ,60 + ,33 + ,0 + ,1439 + ,57635 + ,237213 + ,78 + ,74 + ,123 + ,84 + ,0 + ,1764 + ,66198 + ,173326 + ,86 + ,49 + ,148 + ,88 + ,7 + ,1495 + ,71701 + ,133131 + ,44 + ,52 + ,90 + ,55 + ,7 + ,1373 + ,57793 + ,258873 + ,104 + ,88 + ,124 + ,60 + ,3 + ,2187 + ,80444 + ,180083 + ,63 + ,36 + ,70 + ,66 + ,9 + ,1491 + ,53855 + ,324799 + ,158 + ,108 + ,168 + ,154 + ,0 + ,4041 + ,97668 + ,230964 + ,102 + ,43 + ,115 + ,53 + ,4 + ,1706 + ,133824 + ,236785 + ,77 + ,75 + ,71 + ,119 + ,3 + ,2152 + ,101481 + ,135473 + ,82 + ,32 + ,66 + ,41 + ,0 + ,1036 + ,99645 + ,202925 + ,115 + ,44 + ,134 + ,61 + ,7 + ,1882 + ,114789 + ,215147 + ,101 + ,85 + ,117 + ,58 + ,0 + ,1929 + ,99052 + ,344297 + ,80 + ,86 + ,108 + ,75 + ,1 + ,2242 + ,67654 + ,153935 + ,50 + ,56 + ,84 + ,33 + ,5 + ,1220 + ,65553 + ,132943 + ,83 + ,50 + ,156 + ,40 + ,7 + ,1289 + ,97500 + ,174724 + ,123 + ,135 + ,120 + ,92 + ,0 + ,2515 + ,69112 + ,174415 + ,73 + ,63 + ,114 + ,100 + ,0 + ,2147 + ,82753 + ,225548 + ,81 + ,81 + ,94 + ,112 + ,5 + ,2352 + ,85323 + ,223632 + ,105 + ,52 + ,120 + ,73 + ,0 + ,1638 + ,72654 + ,124817 + ,47 + ,44 + ,81 + ,40 + ,0 + ,1222 + ,30727 + ,221698 + ,105 + ,113 + ,110 + ,45 + ,0 + ,1812 + ,77873 + ,210767 + ,94 + ,39 + ,133 + ,60 + ,3 + ,1677 + ,117478 + ,170266 + ,44 + ,73 + ,122 + ,62 + ,4 + ,1579 + ,74007 + ,260561 + ,114 + ,48 + ,158 + ,75 + ,1 + ,1731 + ,90183 + ,84853 + ,38 + ,33 + ,109 + ,31 + ,4 + ,807 + ,61542 + ,294424 + ,107 + ,59 + ,124 + ,77 + ,2 + ,2452 + ,101494 + ,101011 + ,30 + ,41 + ,39 + ,34 + ,0 + ,829 + ,27570 + ,215641 + ,71 + ,69 + ,92 + ,46 + ,0 + ,1940 + ,55813 + ,325107 + ,84 + ,64 + ,126 + ,99 + ,0 + ,2662 + ,79215 + ,7176 + ,0 + ,1 + ,0 + ,17 + ,0 + ,186 + ,1423 + ,167542 + ,59 + ,59 + ,70 + ,66 + ,2 + ,1499 + ,55461 + ,106408 + ,33 + ,32 + ,37 + ,30 + ,1 + ,865 + ,31081 + ,96560 + ,42 + ,129 + ,38 + ,76 + ,0 + ,1793 + ,22996 + ,265769 + ,96 + ,37 + ,120 + ,146 + ,2 + ,2527 + ,83122 + ,269651 + ,106 + ,31 + ,93 + ,67 + ,10 + ,2747 + ,70106 + ,149112 + ,56 + ,65 + ,95 + ,56 + ,6 + ,1324 + ,60578 + ,175824 + ,57 + ,107 + ,77 + ,107 + ,0 + ,2702 + ,39992 + ,152871 + ,59 + ,74 + ,90 + ,58 + ,5 + ,1383 + ,79892 + ,111665 + ,39 + ,54 + ,80 + ,34 + ,4 + ,1179 + ,49810 + ,116408 + ,34 + ,76 + ,31 + ,61 + ,1 + ,2099 + ,71570 + ,362301 + ,76 + ,715 + ,110 + ,119 + ,2 + ,4308 + ,100708 + ,78800 + ,20 + ,57 + ,66 + ,42 + ,2 + ,918 + ,33032 + ,183167 + ,91 + ,66 + ,138 + ,66 + ,0 + ,1831 + ,82875 + ,277965 + ,115 + ,106 + ,133 + ,89 + ,8 + ,3373 + ,139077 + ,150629 + ,85 + ,54 + ,113 + ,44 + ,3 + ,1713 + ,71595 + ,168809 + ,76 + ,32 + ,100 + ,66 + ,0 + ,1438 + ,72260 + ,24188 + ,8 + ,20 + ,7 + ,24 + ,0 + ,496 + ,5950 + ,329267 + ,79 + ,71 + ,140 + ,259 + ,8 + ,2253 + ,115762 + ,65029 + ,21 + ,21 + ,61 + ,17 + ,5 + ,744 + ,32551 + ,101097 + ,30 + ,70 + ,41 + ,64 + ,3 + ,1161 + ,31701 + ,218946 + ,76 + ,112 + ,96 + ,41 + ,1 + ,2352 + ,80670 + ,244052 + ,101 + ,66 + ,164 + ,68 + ,5 + ,2144 + ,143558 + ,341570 + ,94 + ,190 + ,78 + ,168 + ,1 + ,4691 + ,117105 + ,103597 + ,27 + ,66 + ,49 + ,43 + ,1 + ,1112 + ,23789 + ,233328 + ,92 + ,165 + ,102 + ,132 + ,5 + ,2694 + ,120733 + ,256462 + ,123 + ,56 + ,124 + ,105 + ,0 + ,1973 + ,105195 + ,206161 + ,75 + ,61 + ,99 + ,71 + ,12 + ,1769 + ,73107 + ,311473 + ,128 + ,53 + ,129 + ,112 + ,8 + ,3148 + ,132068 + ,235800 + ,105 + ,127 + ,62 + ,94 + ,8 + ,2474 + ,149193 + ,177939 + ,55 + ,63 + ,73 + ,82 + ,8 + ,2084 + ,46821 + ,207176 + ,56 + ,38 + ,114 + ,70 + ,8 + ,1954 + ,87011 + ,196553 + ,41 + ,50 + ,99 + ,57 + ,2 + ,1226 + ,95260 + ,174184 + ,72 + ,52 + ,70 + ,53 + ,0 + ,1389 + ,55183 + ,143246 + ,67 + ,42 + ,104 + ,103 + ,5 + ,1496 + ,106671 + ,187559 + ,75 + ,76 + ,116 + ,121 + ,8 + ,2269 + ,73511 + ,187681 + ,114 + ,67 + ,91 + ,62 + ,2 + ,1833 + ,92945 + ,119016 + ,118 + ,50 + ,74 + ,52 + ,5 + ,1268 + ,78664 + ,182192 + ,77 + ,53 + ,138 + ,52 + ,12 + ,1943 + ,70054 + ,73566 + ,22 + ,39 + ,67 + ,32 + ,6 + ,893 + ,22618 + ,194979 + ,66 + ,50 + ,151 + ,62 + ,7 + ,1762 + ,74011 + ,167488 + ,69 + ,77 + ,72 + ,45 + ,2 + ,1403 + ,83737 + ,143756 + ,105 + ,57 + ,120 + ,46 + ,0 + ,1425 + ,69094 + ,275541 + ,116 + ,73 + ,115 + ,63 + ,4 + ,1857 + ,93133 + ,243199 + ,88 + ,34 + ,105 + ,75 + ,3 + ,1840 + ,95536 + ,182999 + ,73 + ,39 + ,104 + ,88 + ,6 + ,1502 + ,225920 + ,135649 + ,99 + ,46 + ,108 + ,46 + ,2 + ,1441 + ,62133 + ,152299 + ,62 + ,63 + ,98 + ,53 + ,0 + ,1420 + ,61370 + ,120221 + ,53 + ,35 + ,69 + ,37 + ,1 + ,1416 + ,43836 + ,346485 + ,118 + ,106 + ,111 + ,90 + ,0 + ,2970 + ,106117 + ,145790 + ,30 + ,43 + ,99 + ,63 + ,5 + ,1317 + ,38692 + ,193339 + ,100 + ,47 + ,71 + ,78 + ,2 + ,1644 + ,84651 + ,80953 + ,49 + ,31 + ,27 + ,25 + ,0 + ,870 + ,56622 + ,122774 + ,24 + ,162 + ,69 + ,45 + ,0 + ,1654 + ,15986 + ,130585 + ,67 + ,57 + ,107 + ,46 + ,5 + ,1054 + ,95364 + ,112611 + ,46 + ,36 + ,73 + ,41 + ,0 + ,937 + ,26706 + ,286468 + ,57 + ,263 + ,107 + ,144 + ,1 + ,3004 + ,89691 + ,241066 + ,75 + ,78 + ,93 + ,82 + ,0 + ,2008 + ,67267 + ,148446 + ,135 + ,63 + ,129 + ,91 + ,1 + ,2547 + ,126846 + ,204713 + ,68 + ,54 + ,69 + ,71 + ,1 + ,1885 + ,41140 + ,182079 + ,124 + ,63 + ,118 + ,63 + ,2 + ,1626 + ,102860 + ,140344 + ,33 + ,77 + ,73 + ,53 + ,6 + ,1468 + ,51715 + ,220516 + ,98 + ,79 + ,119 + ,62 + ,1 + ,2445 + ,55801 + ,243060 + ,58 + ,110 + ,104 + ,63 + ,4 + ,1964 + ,111813 + ,162765 + ,68 + ,56 + ,107 + ,32 + ,2 + ,1381 + ,120293 + ,182613 + ,81 + ,56 + ,99 + ,39 + ,3 + ,1369 + ,138599 + ,232138 + ,131 + ,43 + ,90 + ,62 + ,0 + ,1659 + ,161647 + ,265318 + ,110 + ,111 + ,197 + ,117 + ,10 + ,2888 + ,115929 + ,85574 + ,37 + ,71 + ,36 + ,34 + ,0 + ,1290 + ,24266 + ,310839 + ,130 + ,62 + ,85 + ,92 + ,9 + ,2845 + ,162901 + ,225060 + ,93 + ,56 + ,139 + ,93 + ,7 + ,1982 + ,109825 + ,232317 + ,118 + ,74 + ,106 + ,54 + ,0 + ,1904 + ,129838 + ,144966 + ,39 + ,60 + ,50 + ,144 + ,0 + ,1391 + ,37510 + ,43287 + ,13 + ,43 + ,64 + ,14 + ,4 + ,602 + ,43750 + ,155754 + ,74 + ,68 + ,31 + ,61 + ,4 + ,1743 + ,40652 + ,164709 + ,81 + ,53 + ,63 + ,109 + ,0 + ,1559 + ,87771 + ,201940 + ,109 + ,87 + ,92 + ,38 + ,0 + ,2014 + ,85872 + ,235454 + ,151 + ,46 + ,106 + ,73 + ,0 + ,2143 + ,89275 + ,220801 + ,51 + ,105 + ,63 + ,75 + ,1 + ,2146 + ,44418 + ,99466 + ,28 + ,32 + ,69 + ,50 + ,0 + ,874 + ,192565 + ,92661 + ,40 + ,133 + ,41 + ,61 + ,1 + ,1590 + ,35232 + ,133328 + ,56 + ,79 + ,56 + ,55 + ,0 + ,1590 + ,40909 + ,61361 + ,27 + ,51 + ,25 + ,77 + ,0 + ,1210 + ,13294 + ,125930 + ,37 + ,207 + ,65 + ,75 + ,4 + ,2072 + ,32387 + ,100750 + ,83 + ,67 + ,93 + ,72 + ,0 + ,1281 + ,140867 + ,224549 + ,54 + ,47 + ,114 + ,50 + ,4 + ,1401 + ,120662 + ,82316 + ,27 + ,34 + ,38 + ,32 + ,4 + ,834 + ,21233 + ,102010 + ,28 + ,66 + ,44 + ,53 + ,3 + ,1105 + ,44332 + ,101523 + ,59 + ,76 + ,87 + ,42 + ,0 + ,1272 + ,61056 + ,243511 + ,133 + ,65 + ,110 + ,71 + ,0 + ,1944 + ,101338 + ,22938 + ,12 + ,9 + ,0 + ,10 + ,0 + ,391 + ,1168 + ,41566 + ,0 + ,42 + ,27 + ,35 + ,5 + ,761 + ,13497 + ,152474 + ,106 + ,45 + ,83 + ,65 + ,0 + ,1605 + ,65567 + ,61857 + ,23 + ,25 + ,30 + ,25 + ,4 + ,530 + ,25162 + ,99923 + ,44 + ,115 + ,80 + ,66 + ,0 + ,1988 + ,32334 + ,132487 + ,71 + ,97 + ,98 + ,41 + ,0 + ,1386 + ,40735 + ,317394 + ,116 + ,53 + ,82 + ,86 + ,1 + ,2395 + ,91413 + ,21054 + ,4 + ,2 + ,0 + ,16 + ,0 + ,387 + ,855 + ,209641 + ,62 + ,52 + ,60 + ,42 + ,5 + ,1742 + ,97068 + ,22648 + ,12 + ,44 + ,28 + ,19 + ,0 + ,620 + ,44339 + ,31414 + ,18 + ,22 + ,9 + ,19 + ,0 + ,449 + ,14116 + ,46698 + ,14 + ,35 + ,33 + ,45 + ,0 + ,800 + ,10288 + ,131698 + ,60 + ,74 + ,59 + ,65 + ,0 + ,1684 + ,65622 + ,91735 + ,7 + ,103 + ,49 + ,35 + ,0 + ,1050 + ,16563 + ,244749 + ,98 + ,144 + ,115 + ,95 + ,2 + ,2699 + ,76643 + ,184510 + ,64 + ,60 + ,140 + ,49 + ,7 + ,1606 + ,110681 + ,79863 + ,29 + ,134 + ,49 + ,37 + ,1 + ,1502 + ,29011 + ,128423 + ,32 + ,89 + ,120 + ,64 + ,8 + ,1204 + ,92696 + ,97839 + ,25 + ,42 + ,66 + ,38 + ,2 + ,1138 + ,94785 + ,38214 + ,16 + ,52 + ,21 + ,34 + ,0 + ,568 + ,8773 + ,151101 + ,48 + ,98 + ,124 + ,32 + ,2 + ,1459 + ,83209 + ,272458 + ,100 + ,99 + ,152 + ,65 + ,0 + ,2158 + ,93815 + ,172494 + ,46 + ,52 + ,139 + ,52 + ,0 + ,1111 + ,86687 + ,108043 + ,45 + ,29 + ,38 + ,62 + ,1 + ,1421 + ,34553 + ,328107 + ,129 + ,125 + ,144 + ,65 + ,3 + ,2833 + ,105547 + ,250579 + ,130 + ,106 + ,120 + ,83 + ,0 + ,1955 + ,103487 + ,351067 + ,136 + ,95 + ,160 + ,95 + ,3 + ,2922 + ,213688 + ,158015 + ,59 + ,40 + ,114 + ,29 + ,0 + ,1002 + ,71220 + ,98866 + ,25 + ,140 + ,39 + ,18 + ,0 + ,1060 + ,23517 + ,85439 + ,32 + ,43 + ,78 + ,33 + ,0 + ,956 + ,56926 + ,229242 + ,63 + ,128 + ,119 + ,247 + ,4 + ,2186 + ,91721 + ,351619 + ,95 + ,142 + ,141 + ,139 + ,4 + ,3604 + ,115168 + ,84207 + ,14 + ,73 + ,101 + ,29 + ,11 + ,1035 + ,111194 + ,120445 + ,36 + ,72 + ,56 + ,118 + ,0 + ,1417 + ,51009 + ,324598 + ,113 + ,128 + ,133 + ,110 + ,0 + ,3261 + ,135777 + ,131069 + ,47 + ,61 + ,83 + ,67 + ,4 + ,1587 + ,51513 + ,204271 + ,92 + ,73 + ,116 + ,42 + ,0 + ,1424 + ,74163 + ,165543 + ,70 + ,148 + ,90 + ,65 + ,1 + ,1701 + ,51633 + ,141722 + ,19 + ,64 + ,36 + ,94 + ,0 + ,1249 + ,75345 + ,116048 + ,50 + ,45 + ,50 + ,64 + ,0 + ,946 + ,33416 + ,250047 + ,41 + ,58 + ,61 + ,81 + ,0 + ,1926 + ,83305 + ,299775 + ,91 + ,97 + ,97 + ,95 + ,9 + ,3352 + ,98952 + ,195838 + ,111 + ,50 + ,98 + ,67 + ,1 + ,1641 + ,102372 + ,173260 + ,41 + ,37 + ,78 + ,63 + ,3 + ,2035 + ,37238 + ,254488 + ,120 + ,50 + ,117 + ,83 + ,10 + ,2312 + ,103772 + ,104389 + ,135 + ,105 + ,148 + ,45 + ,5 + ,1369 + ,123969 + ,136084 + ,27 + ,69 + ,41 + ,30 + ,0 + ,1577 + ,27142 + ,199476 + ,87 + ,46 + ,105 + ,70 + ,2 + ,2201 + ,135400 + ,92499 + ,25 + ,57 + ,55 + ,32 + ,0 + ,961 + ,21399 + ,224330 + ,131 + ,52 + ,132 + ,83 + ,1 + ,1900 + ,130115 + ,135781 + ,45 + ,98 + ,44 + ,31 + ,2 + ,1254 + ,24874 + ,74408 + ,29 + ,61 + ,21 + ,67 + ,4 + ,1335 + ,34988 + ,81240 + ,58 + ,89 + ,50 + ,66 + ,0 + ,1597 + ,45549 + ,14688 + ,4 + ,0 + ,0 + ,10 + ,0 + ,207 + ,6023 + ,181633 + ,47 + ,48 + ,73 + ,70 + ,2 + ,1645 + ,64466 + ,271856 + ,109 + ,91 + ,86 + ,103 + ,1 + ,2429 + ,54990 + ,7199 + ,7 + ,0 + ,0 + ,5 + ,0 + ,151 + ,1644 + ,46660 + ,12 + ,7 + ,13 + ,20 + ,0 + ,474 + ,6179 + ,17547 + ,0 + ,3 + ,4 + ,5 + ,0 + ,141 + ,3926 + ,133368 + ,37 + ,54 + ,57 + ,36 + ,1 + ,1639 + ,32755 + ,95227 + ,37 + ,70 + ,48 + ,34 + ,0 + ,872 + ,34777 + ,152601 + ,46 + ,36 + ,46 + ,48 + ,2 + ,1318 + ,73224 + ,98146 + ,15 + ,37 + ,48 + ,40 + ,0 + ,1018 + ,27114 + ,79619 + ,42 + ,123 + ,32 + ,43 + ,3 + ,1383 + ,20760 + ,59194 + ,7 + ,247 + ,68 + ,31 + ,6 + ,1314 + ,37636 + ,139942 + ,54 + ,46 + ,87 + ,42 + ,0 + ,1335 + ,65461 + ,118612 + ,54 + ,72 + ,43 + ,46 + ,2 + ,1403 + ,30080 + ,72880 + ,14 + ,41 + ,67 + ,33 + ,0 + ,910 + ,24094 + ,65475 + ,16 + ,24 + ,46 + ,18 + ,2 + ,616 + ,69008 + ,99643 + ,33 + ,45 + ,46 + ,55 + ,1 + ,1407 + ,54968 + ,71965 + ,32 + ,33 + ,56 + ,35 + ,1 + ,771 + ,46090 + ,77272 + ,21 + ,27 + ,48 + ,59 + ,2 + ,766 + ,27507 + ,49289 + ,15 + ,36 + ,44 + ,19 + ,1 + ,473 + ,10672 + ,135131 + ,38 + ,87 + ,60 + ,66 + ,0 + ,1376 + ,34029 + ,108446 + ,22 + ,90 + ,65 + ,60 + ,1 + ,1232 + ,46300 + ,89746 + ,28 + ,114 + ,55 + ,36 + ,3 + ,1521 + ,24760 + ,44296 + ,10 + ,31 + ,38 + ,25 + ,0 + ,572 + ,18779 + ,77648 + ,31 + ,45 + ,52 + ,47 + ,0 + ,1059 + ,21280 + ,181528 + ,32 + ,69 + ,60 + ,54 + ,0 + ,1544 + ,40662 + ,134019 + ,32 + ,51 + ,54 + ,53 + ,0 + ,1230 + ,28987 + ,124064 + ,43 + ,34 + ,86 + ,40 + ,1 + ,1206 + ,22827 + ,92630 + ,27 + ,60 + ,24 + ,40 + ,4 + ,1205 + ,18513 + ,121848 + ,37 + ,45 + ,52 + ,39 + ,0 + ,1255 + ,30594 + ,52915 + ,20 + ,54 + ,49 + ,14 + ,0 + ,613 + ,24006 + ,81872 + ,32 + ,25 + ,61 + ,45 + ,0 + ,721 + ,27913 + ,58981 + ,0 + ,38 + ,61 + ,36 + ,7 + ,1109 + ,42744 + ,53515 + ,5 + ,52 + ,81 + ,28 + ,2 + ,740 + ,12934 + ,60812 + ,26 + ,67 + ,43 + ,44 + ,0 + ,1126 + ,22574 + ,56375 + ,10 + ,74 + ,40 + ,30 + ,7 + ,728 + ,41385 + ,65490 + ,27 + ,38 + ,40 + ,22 + ,3 + ,689 + ,18653 + ,80949 + ,11 + ,30 + ,56 + ,17 + ,0 + ,592 + ,18472 + ,76302 + ,29 + ,26 + ,68 + ,31 + ,0 + ,995 + ,30976 + ,104011 + ,25 + ,67 + ,79 + ,55 + ,6 + ,1613 + ,63339 + ,98104 + ,55 + ,132 + ,47 + ,54 + ,2 + ,2048 + ,25568 + ,67989 + ,23 + ,42 + ,57 + ,21 + ,0 + ,705 + ,33747 + ,30989 + ,5 + ,35 + ,41 + ,14 + ,0 + ,301 + ,4154 + ,135458 + ,43 + ,118 + ,29 + ,81 + ,3 + ,1803 + ,19474 + ,73504 + ,23 + ,68 + ,3 + ,35 + ,0 + ,799 + ,35130 + ,63123 + ,34 + ,43 + ,60 + ,43 + ,1 + ,861 + ,39067 + ,61254 + ,36 + ,76 + ,30 + ,46 + ,1 + ,1186 + ,13310 + ,74914 + ,35 + ,64 + ,79 + ,30 + ,0 + ,1451 + ,65892 + ,31774 + ,0 + ,48 + ,47 + ,23 + ,1 + ,628 + ,4143 + ,81437 + ,37 + ,64 + ,40 + ,38 + ,0 + ,1161 + ,28579 + ,87186 + ,28 + ,56 + ,48 + ,54 + ,0 + ,1463 + ,51776 + ,50090 + ,16 + ,71 + ,36 + ,20 + ,0 + ,742 + ,21152 + ,65745 + ,26 + ,75 + ,42 + ,53 + ,0 + ,979 + ,38084 + ,56653 + ,38 + ,39 + ,49 + ,45 + ,0 + ,675 + ,27717 + ,158399 + ,23 + ,42 + ,57 + ,39 + ,0 + ,1241 + ,32928 + ,46455 + ,22 + ,39 + ,12 + ,20 + ,0 + ,676 + ,11342 + ,73624 + ,30 + ,93 + ,40 + ,24 + ,0 + ,1049 + ,19499 + ,38395 + ,16 + ,38 + ,43 + ,31 + ,0 + ,620 + ,16380 + ,91899 + ,18 + ,60 + ,33 + ,35 + ,0 + ,1081 + ,36874 + ,139526 + ,28 + ,71 + ,77 + ,151 + ,0 + ,1688 + ,48259 + ,52164 + ,32 + ,52 + ,43 + ,52 + ,0 + ,736 + ,16734 + ,51567 + ,21 + ,27 + ,45 + ,30 + ,2 + ,617 + ,28207 + ,70551 + ,23 + ,59 + ,47 + ,31 + ,0 + ,812 + ,30143 + ,84856 + ,29 + ,40 + ,43 + ,29 + ,1 + ,1051 + ,41369 + ,102538 + ,50 + ,79 + ,45 + ,57 + ,1 + ,1656 + ,45833 + ,86678 + ,12 + ,44 + ,50 + ,40 + ,0 + ,705 + ,29156 + ,85709 + ,21 + ,65 + ,35 + ,44 + ,0 + ,945 + ,35944 + ,34662 + ,18 + ,10 + ,7 + ,25 + ,0 + ,554 + ,36278 + ,150580 + ,27 + ,124 + ,71 + ,77 + ,0 + ,1597 + ,45588 + ,99611 + ,41 + ,81 + ,67 + ,35 + ,0 + ,982 + ,45097 + ,19349 + ,13 + ,15 + ,0 + ,11 + ,0 + ,222 + ,3895 + ,99373 + ,12 + ,92 + ,62 + ,63 + ,1 + ,1212 + ,28394 + ,86230 + ,21 + ,42 + ,54 + ,44 + ,0 + ,1143 + ,18632 + ,30837 + ,8 + ,10 + ,4 + ,19 + ,0 + ,435 + ,2325 + ,31706 + ,26 + ,24 + ,25 + ,13 + ,0 + ,532 + ,25139 + ,89806 + ,27 + ,64 + ,40 + ,42 + ,0 + ,882 + ,27975 + ,62088 + ,13 + ,45 + ,38 + ,38 + ,1 + ,608 + ,14483 + ,40151 + ,16 + ,22 + ,19 + ,29 + ,0 + ,459 + ,13127 + ,27634 + ,2 + ,56 + ,17 + ,20 + ,0 + ,578 + ,5839 + ,76990 + ,42 + ,94 + ,67 + ,27 + ,0 + ,826 + ,24069 + ,37460 + ,5 + ,19 + ,14 + ,20 + ,0 + ,509 + ,3738 + ,54157 + ,37 + ,35 + ,30 + ,19 + ,0 + ,717 + ,18625 + ,49862 + ,17 + ,32 + ,54 + ,37 + ,0 + ,637 + ,36341 + ,84337 + ,38 + ,35 + ,35 + ,26 + ,0 + ,857 + ,24548 + ,64175 + ,37 + ,48 + ,59 + ,42 + ,0 + ,830 + ,21792 + ,59382 + ,29 + ,49 + ,24 + ,49 + ,0 + ,652 + ,26263 + ,119308 + ,32 + ,48 + ,58 + ,30 + ,0 + ,707 + ,23686 + ,76702 + ,35 + ,62 + ,42 + ,49 + ,0 + ,954 + ,49303 + ,103425 + ,17 + ,96 + ,46 + ,67 + ,1 + ,1461 + ,25659 + ,70344 + ,20 + ,45 + ,61 + ,28 + ,0 + ,672 + ,28904 + ,43410 + ,7 + ,63 + ,3 + ,19 + ,0 + ,778 + ,2781 + ,104838 + ,46 + ,71 + ,52 + ,49 + ,1 + ,1141 + ,29236 + ,62215 + ,24 + ,26 + ,25 + ,27 + ,0 + ,680 + ,19546 + ,69304 + ,40 + ,48 + ,40 + ,30 + ,6 + ,1090 + ,22818 + ,53117 + ,3 + ,29 + ,32 + ,22 + ,3 + ,616 + ,32689 + ,19764 + ,10 + ,19 + ,4 + ,12 + ,1 + ,285 + ,5752 + ,86680 + ,37 + ,45 + ,49 + ,31 + ,2 + ,1145 + ,22197 + ,84105 + ,17 + ,45 + ,63 + ,20 + ,0 + ,733 + ,20055 + ,77945 + ,28 + ,67 + ,67 + ,20 + ,0 + ,888 + ,25272 + ,89113 + ,19 + ,30 + ,32 + ,39 + ,0 + ,849 + ,82206 + ,91005 + ,29 + ,36 + ,23 + ,29 + ,3 + ,1182 + ,32073 + ,40248 + ,8 + ,34 + ,7 + ,16 + ,1 + ,528 + ,5444 + ,64187 + ,10 + ,36 + ,54 + ,27 + ,0 + ,642 + ,20154 + ,50857 + ,15 + ,34 + ,37 + ,21 + ,0 + ,947 + ,36944 + ,56613 + ,15 + ,37 + ,35 + ,19 + ,1 + ,819 + ,8019 + ,62792 + ,28 + ,46 + ,51 + ,35 + ,0 + ,757 + ,30884 + ,72535 + ,17 + ,44 + ,39 + ,14 + ,0 + ,894 + ,19540) + ,dim=c(8 + ,289) + ,dimnames=list(c('time_in_rfc' + ,'blogged_computations' + ,'compendium_views_pr' + ,'feedback_messages_p120' + ,'logins' + ,'shared_compendiums' + ,'pageviews' + ,'totsize ') + ,1:289)) > y <- array(NA,dim=c(8,289),dimnames=list(c('time_in_rfc','blogged_computations','compendium_views_pr','feedback_messages_p120','logins','shared_compendiums','pageviews','totsize '),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 = '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 > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x time_in_rfc blogged_computations compendium_views_pr feedback_messages_p120 1 210907 79 81 94 2 120982 58 55 103 3 176508 60 50 93 4 179321 108 125 103 5 123185 49 40 51 6 52746 0 37 70 7 385534 121 63 91 8 33170 1 44 22 9 101645 20 88 38 10 149061 43 66 93 11 165446 69 57 60 12 237213 78 74 123 13 173326 86 49 148 14 133131 44 52 90 15 258873 104 88 124 16 180083 63 36 70 17 324799 158 108 168 18 230964 102 43 115 19 236785 77 75 71 20 135473 82 32 66 21 202925 115 44 134 22 215147 101 85 117 23 344297 80 86 108 24 153935 50 56 84 25 132943 83 50 156 26 174724 123 135 120 27 174415 73 63 114 28 225548 81 81 94 29 223632 105 52 120 30 124817 47 44 81 31 221698 105 113 110 32 210767 94 39 133 33 170266 44 73 122 34 260561 114 48 158 35 84853 38 33 109 36 294424 107 59 124 37 101011 30 41 39 38 215641 71 69 92 39 325107 84 64 126 40 7176 0 1 0 41 167542 59 59 70 42 106408 33 32 37 43 96560 42 129 38 44 265769 96 37 120 45 269651 106 31 93 46 149112 56 65 95 47 175824 57 107 77 48 152871 59 74 90 49 111665 39 54 80 50 116408 34 76 31 51 362301 76 715 110 52 78800 20 57 66 53 183167 91 66 138 54 277965 115 106 133 55 150629 85 54 113 56 168809 76 32 100 57 24188 8 20 7 58 329267 79 71 140 59 65029 21 21 61 60 101097 30 70 41 61 218946 76 112 96 62 244052 101 66 164 63 341570 94 190 78 64 103597 27 66 49 65 233328 92 165 102 66 256462 123 56 124 67 206161 75 61 99 68 311473 128 53 129 69 235800 105 127 62 70 177939 55 63 73 71 207176 56 38 114 72 196553 41 50 99 73 174184 72 52 70 74 143246 67 42 104 75 187559 75 76 116 76 187681 114 67 91 77 119016 118 50 74 78 182192 77 53 138 79 73566 22 39 67 80 194979 66 50 151 81 167488 69 77 72 82 143756 105 57 120 83 275541 116 73 115 84 243199 88 34 105 85 182999 73 39 104 86 135649 99 46 108 87 152299 62 63 98 88 120221 53 35 69 89 346485 118 106 111 90 145790 30 43 99 91 193339 100 47 71 92 80953 49 31 27 93 122774 24 162 69 94 130585 67 57 107 95 112611 46 36 73 96 286468 57 263 107 97 241066 75 78 93 98 148446 135 63 129 99 204713 68 54 69 100 182079 124 63 118 101 140344 33 77 73 102 220516 98 79 119 103 243060 58 110 104 104 162765 68 56 107 105 182613 81 56 99 106 232138 131 43 90 107 265318 110 111 197 108 85574 37 71 36 109 310839 130 62 85 110 225060 93 56 139 111 232317 118 74 106 112 144966 39 60 50 113 43287 13 43 64 114 155754 74 68 31 115 164709 81 53 63 116 201940 109 87 92 117 235454 151 46 106 118 220801 51 105 63 119 99466 28 32 69 120 92661 40 133 41 121 133328 56 79 56 122 61361 27 51 25 123 125930 37 207 65 124 100750 83 67 93 125 224549 54 47 114 126 82316 27 34 38 127 102010 28 66 44 128 101523 59 76 87 129 243511 133 65 110 130 22938 12 9 0 131 41566 0 42 27 132 152474 106 45 83 133 61857 23 25 30 134 99923 44 115 80 135 132487 71 97 98 136 317394 116 53 82 137 21054 4 2 0 138 209641 62 52 60 139 22648 12 44 28 140 31414 18 22 9 141 46698 14 35 33 142 131698 60 74 59 143 91735 7 103 49 144 244749 98 144 115 145 184510 64 60 140 146 79863 29 134 49 147 128423 32 89 120 148 97839 25 42 66 149 38214 16 52 21 150 151101 48 98 124 151 272458 100 99 152 152 172494 46 52 139 153 108043 45 29 38 154 328107 129 125 144 155 250579 130 106 120 156 351067 136 95 160 157 158015 59 40 114 158 98866 25 140 39 159 85439 32 43 78 160 229242 63 128 119 161 351619 95 142 141 162 84207 14 73 101 163 120445 36 72 56 164 324598 113 128 133 165 131069 47 61 83 166 204271 92 73 116 167 165543 70 148 90 168 141722 19 64 36 169 116048 50 45 50 170 250047 41 58 61 171 299775 91 97 97 172 195838 111 50 98 173 173260 41 37 78 174 254488 120 50 117 175 104389 135 105 148 176 136084 27 69 41 177 199476 87 46 105 178 92499 25 57 55 179 224330 131 52 132 180 135781 45 98 44 181 74408 29 61 21 182 81240 58 89 50 183 14688 4 0 0 184 181633 47 48 73 185 271856 109 91 86 186 7199 7 0 0 187 46660 12 7 13 188 17547 0 3 4 189 133368 37 54 57 190 95227 37 70 48 191 152601 46 36 46 192 98146 15 37 48 193 79619 42 123 32 194 59194 7 247 68 195 139942 54 46 87 196 118612 54 72 43 197 72880 14 41 67 198 65475 16 24 46 199 99643 33 45 46 200 71965 32 33 56 201 77272 21 27 48 202 49289 15 36 44 203 135131 38 87 60 204 108446 22 90 65 205 89746 28 114 55 206 44296 10 31 38 207 77648 31 45 52 208 181528 32 69 60 209 134019 32 51 54 210 124064 43 34 86 211 92630 27 60 24 212 121848 37 45 52 213 52915 20 54 49 214 81872 32 25 61 215 58981 0 38 61 216 53515 5 52 81 217 60812 26 67 43 218 56375 10 74 40 219 65490 27 38 40 220 80949 11 30 56 221 76302 29 26 68 222 104011 25 67 79 223 98104 55 132 47 224 67989 23 42 57 225 30989 5 35 41 226 135458 43 118 29 227 73504 23 68 3 228 63123 34 43 60 229 61254 36 76 30 230 74914 35 64 79 231 31774 0 48 47 232 81437 37 64 40 233 87186 28 56 48 234 50090 16 71 36 235 65745 26 75 42 236 56653 38 39 49 237 158399 23 42 57 238 46455 22 39 12 239 73624 30 93 40 240 38395 16 38 43 241 91899 18 60 33 242 139526 28 71 77 243 52164 32 52 43 244 51567 21 27 45 245 70551 23 59 47 246 84856 29 40 43 247 102538 50 79 45 248 86678 12 44 50 249 85709 21 65 35 250 34662 18 10 7 251 150580 27 124 71 252 99611 41 81 67 253 19349 13 15 0 254 99373 12 92 62 255 86230 21 42 54 256 30837 8 10 4 257 31706 26 24 25 258 89806 27 64 40 259 62088 13 45 38 260 40151 16 22 19 261 27634 2 56 17 262 76990 42 94 67 263 37460 5 19 14 264 54157 37 35 30 265 49862 17 32 54 266 84337 38 35 35 267 64175 37 48 59 268 59382 29 49 24 269 119308 32 48 58 270 76702 35 62 42 271 103425 17 96 46 272 70344 20 45 61 273 43410 7 63 3 274 104838 46 71 52 275 62215 24 26 25 276 69304 40 48 40 277 53117 3 29 32 278 19764 10 19 4 279 86680 37 45 49 280 84105 17 45 63 281 77945 28 67 67 282 89113 19 30 32 283 91005 29 36 23 284 40248 8 34 7 285 64187 10 36 54 286 50857 15 34 37 287 56613 15 37 35 288 62792 28 46 51 289 72535 17 44 39 logins shared_compendiums pageviews totsize\r\r\r t 1 56 3 1418 112285 1 2 56 4 869 84786 2 3 54 12 1530 83123 3 4 89 2 2172 101193 4 5 40 1 901 38361 5 6 25 3 463 68504 6 7 92 0 3201 119182 7 8 18 0 371 22807 8 9 63 0 1192 17140 9 10 44 5 1583 116174 10 11 33 0 1439 57635 11 12 84 0 1764 66198 12 13 88 7 1495 71701 13 14 55 7 1373 57793 14 15 60 3 2187 80444 15 16 66 9 1491 53855 16 17 154 0 4041 97668 17 18 53 4 1706 133824 18 19 119 3 2152 101481 19 20 41 0 1036 99645 20 21 61 7 1882 114789 21 22 58 0 1929 99052 22 23 75 1 2242 67654 23 24 33 5 1220 65553 24 25 40 7 1289 97500 25 26 92 0 2515 69112 26 27 100 0 2147 82753 27 28 112 5 2352 85323 28 29 73 0 1638 72654 29 30 40 0 1222 30727 30 31 45 0 1812 77873 31 32 60 3 1677 117478 32 33 62 4 1579 74007 33 34 75 1 1731 90183 34 35 31 4 807 61542 35 36 77 2 2452 101494 36 37 34 0 829 27570 37 38 46 0 1940 55813 38 39 99 0 2662 79215 39 40 17 0 186 1423 40 41 66 2 1499 55461 41 42 30 1 865 31081 42 43 76 0 1793 22996 43 44 146 2 2527 83122 44 45 67 10 2747 70106 45 46 56 6 1324 60578 46 47 107 0 2702 39992 47 48 58 5 1383 79892 48 49 34 4 1179 49810 49 50 61 1 2099 71570 50 51 119 2 4308 100708 51 52 42 2 918 33032 52 53 66 0 1831 82875 53 54 89 8 3373 139077 54 55 44 3 1713 71595 55 56 66 0 1438 72260 56 57 24 0 496 5950 57 58 259 8 2253 115762 58 59 17 5 744 32551 59 60 64 3 1161 31701 60 61 41 1 2352 80670 61 62 68 5 2144 143558 62 63 168 1 4691 117105 63 64 43 1 1112 23789 64 65 132 5 2694 120733 65 66 105 0 1973 105195 66 67 71 12 1769 73107 67 68 112 8 3148 132068 68 69 94 8 2474 149193 69 70 82 8 2084 46821 70 71 70 8 1954 87011 71 72 57 2 1226 95260 72 73 53 0 1389 55183 73 74 103 5 1496 106671 74 75 121 8 2269 73511 75 76 62 2 1833 92945 76 77 52 5 1268 78664 77 78 52 12 1943 70054 78 79 32 6 893 22618 79 80 62 7 1762 74011 80 81 45 2 1403 83737 81 82 46 0 1425 69094 82 83 63 4 1857 93133 83 84 75 3 1840 95536 84 85 88 6 1502 225920 85 86 46 2 1441 62133 86 87 53 0 1420 61370 87 88 37 1 1416 43836 88 89 90 0 2970 106117 89 90 63 5 1317 38692 90 91 78 2 1644 84651 91 92 25 0 870 56622 92 93 45 0 1654 15986 93 94 46 5 1054 95364 94 95 41 0 937 26706 95 96 144 1 3004 89691 96 97 82 0 2008 67267 97 98 91 1 2547 126846 98 99 71 1 1885 41140 99 100 63 2 1626 102860 100 101 53 6 1468 51715 101 102 62 1 2445 55801 102 103 63 4 1964 111813 103 104 32 2 1381 120293 104 105 39 3 1369 138599 105 106 62 0 1659 161647 106 107 117 10 2888 115929 107 108 34 0 1290 24266 108 109 92 9 2845 162901 109 110 93 7 1982 109825 110 111 54 0 1904 129838 111 112 144 0 1391 37510 112 113 14 4 602 43750 113 114 61 4 1743 40652 114 115 109 0 1559 87771 115 116 38 0 2014 85872 116 117 73 0 2143 89275 117 118 75 1 2146 44418 118 119 50 0 874 192565 119 120 61 1 1590 35232 120 121 55 0 1590 40909 121 122 77 0 1210 13294 122 123 75 4 2072 32387 123 124 72 0 1281 140867 124 125 50 4 1401 120662 125 126 32 4 834 21233 126 127 53 3 1105 44332 127 128 42 0 1272 61056 128 129 71 0 1944 101338 129 130 10 0 391 1168 130 131 35 5 761 13497 131 132 65 0 1605 65567 132 133 25 4 530 25162 133 134 66 0 1988 32334 134 135 41 0 1386 40735 135 136 86 1 2395 91413 136 137 16 0 387 855 137 138 42 5 1742 97068 138 139 19 0 620 44339 139 140 19 0 449 14116 140 141 45 0 800 10288 141 142 65 0 1684 65622 142 143 35 0 1050 16563 143 144 95 2 2699 76643 144 145 49 7 1606 110681 145 146 37 1 1502 29011 146 147 64 8 1204 92696 147 148 38 2 1138 94785 148 149 34 0 568 8773 149 150 32 2 1459 83209 150 151 65 0 2158 93815 151 152 52 0 1111 86687 152 153 62 1 1421 34553 153 154 65 3 2833 105547 154 155 83 0 1955 103487 155 156 95 3 2922 213688 156 157 29 0 1002 71220 157 158 18 0 1060 23517 158 159 33 0 956 56926 159 160 247 4 2186 91721 160 161 139 4 3604 115168 161 162 29 11 1035 111194 162 163 118 0 1417 51009 163 164 110 0 3261 135777 164 165 67 4 1587 51513 165 166 42 0 1424 74163 166 167 65 1 1701 51633 167 168 94 0 1249 75345 168 169 64 0 946 33416 169 170 81 0 1926 83305 170 171 95 9 3352 98952 171 172 67 1 1641 102372 172 173 63 3 2035 37238 173 174 83 10 2312 103772 174 175 45 5 1369 123969 175 176 30 0 1577 27142 176 177 70 2 2201 135400 177 178 32 0 961 21399 178 179 83 1 1900 130115 179 180 31 2 1254 24874 180 181 67 4 1335 34988 181 182 66 0 1597 45549 182 183 10 0 207 6023 183 184 70 2 1645 64466 184 185 103 1 2429 54990 185 186 5 0 151 1644 186 187 20 0 474 6179 187 188 5 0 141 3926 188 189 36 1 1639 32755 189 190 34 0 872 34777 190 191 48 2 1318 73224 191 192 40 0 1018 27114 192 193 43 3 1383 20760 193 194 31 6 1314 37636 194 195 42 0 1335 65461 195 196 46 2 1403 30080 196 197 33 0 910 24094 197 198 18 2 616 69008 198 199 55 1 1407 54968 199 200 35 1 771 46090 200 201 59 2 766 27507 201 202 19 1 473 10672 202 203 66 0 1376 34029 203 204 60 1 1232 46300 204 205 36 3 1521 24760 205 206 25 0 572 18779 206 207 47 0 1059 21280 207 208 54 0 1544 40662 208 209 53 0 1230 28987 209 210 40 1 1206 22827 210 211 40 4 1205 18513 211 212 39 0 1255 30594 212 213 14 0 613 24006 213 214 45 0 721 27913 214 215 36 7 1109 42744 215 216 28 2 740 12934 216 217 44 0 1126 22574 217 218 30 7 728 41385 218 219 22 3 689 18653 219 220 17 0 592 18472 220 221 31 0 995 30976 221 222 55 6 1613 63339 222 223 54 2 2048 25568 223 224 21 0 705 33747 224 225 14 0 301 4154 225 226 81 3 1803 19474 226 227 35 0 799 35130 227 228 43 1 861 39067 228 229 46 1 1186 13310 229 230 30 0 1451 65892 230 231 23 1 628 4143 231 232 38 0 1161 28579 232 233 54 0 1463 51776 233 234 20 0 742 21152 234 235 53 0 979 38084 235 236 45 0 675 27717 236 237 39 0 1241 32928 237 238 20 0 676 11342 238 239 24 0 1049 19499 239 240 31 0 620 16380 240 241 35 0 1081 36874 241 242 151 0 1688 48259 242 243 52 0 736 16734 243 244 30 2 617 28207 244 245 31 0 812 30143 245 246 29 1 1051 41369 246 247 57 1 1656 45833 247 248 40 0 705 29156 248 249 44 0 945 35944 249 250 25 0 554 36278 250 251 77 0 1597 45588 251 252 35 0 982 45097 252 253 11 0 222 3895 253 254 63 1 1212 28394 254 255 44 0 1143 18632 255 256 19 0 435 2325 256 257 13 0 532 25139 257 258 42 0 882 27975 258 259 38 1 608 14483 259 260 29 0 459 13127 260 261 20 0 578 5839 261 262 27 0 826 24069 262 263 20 0 509 3738 263 264 19 0 717 18625 264 265 37 0 637 36341 265 266 26 0 857 24548 266 267 42 0 830 21792 267 268 49 0 652 26263 268 269 30 0 707 23686 269 270 49 0 954 49303 270 271 67 1 1461 25659 271 272 28 0 672 28904 272 273 19 0 778 2781 273 274 49 1 1141 29236 274 275 27 0 680 19546 275 276 30 6 1090 22818 276 277 22 3 616 32689 277 278 12 1 285 5752 278 279 31 2 1145 22197 279 280 20 0 733 20055 280 281 20 0 888 25272 281 282 39 0 849 82206 282 283 29 3 1182 32073 283 284 16 1 528 5444 284 285 27 0 642 20154 285 286 21 0 947 36944 286 287 19 1 819 8019 287 288 35 0 757 30884 288 289 14 0 894 19540 289 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) blogged_computations compendium_views_pr 857.0455 303.7885 -105.3452 feedback_messages_p120 logins shared_compendiums 315.8860 145.5679 -1527.7978 pageviews `totsize\r\r\r` t 66.2646 0.2503 -53.5537 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -134526 -13771 341 15298 120427 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 8.570e+02 7.068e+03 0.121 0.903579 blogged_computations 3.038e+02 1.001e+02 3.036 0.002622 ** compendium_views_pr -1.053e+02 4.302e+01 -2.449 0.014949 * feedback_messages_p120 3.159e+02 8.001e+01 3.948 9.97e-05 *** logins 1.456e+02 7.874e+01 1.849 0.065560 . shared_compendiums -1.528e+03 7.263e+02 -2.103 0.036321 * pageviews 6.626e+01 5.548e+00 11.944 < 2e-16 *** `totsize\r\r\r` 2.503e-01 7.348e-02 3.406 0.000757 *** t -5.355e+01 2.521e+01 -2.125 0.034488 * --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 29010 on 280 degrees of freedom Multiple R-squared: 0.8793, Adjusted R-squared: 0.8759 F-statistic: 255 on 8 and 280 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.8026480 3.947039e-01 1.973520e-01 [2,] 0.6813588 6.372824e-01 3.186412e-01 [3,] 0.5664261 8.671477e-01 4.335739e-01 [4,] 0.4426078 8.852155e-01 5.573922e-01 [5,] 0.5426066 9.147868e-01 4.573934e-01 [6,] 0.9097224 1.805553e-01 9.027764e-02 [7,] 0.8708865 2.582269e-01 1.291135e-01 [8,] 0.8579714 2.840571e-01 1.420286e-01 [9,] 0.8895884 2.208232e-01 1.104116e-01 [10,] 0.8722575 2.554851e-01 1.277425e-01 [11,] 0.8457848 3.084305e-01 1.542152e-01 [12,] 0.9992363 1.527413e-03 7.637066e-04 [13,] 0.9987298 2.540438e-03 1.270219e-03 [14,] 0.9987319 2.536156e-03 1.268078e-03 [15,] 0.9995995 8.010495e-04 4.005247e-04 [16,] 0.9997623 4.753800e-04 2.376900e-04 [17,] 0.9996037 7.925020e-04 3.962510e-04 [18,] 0.9995992 8.016203e-04 4.008102e-04 [19,] 0.9994042 1.191692e-03 5.958458e-04 [20,] 0.9994487 1.102540e-03 5.512700e-04 [21,] 0.9991086 1.782775e-03 8.913877e-04 [22,] 0.9987348 2.530480e-03 1.265240e-03 [23,] 0.9991637 1.672566e-03 8.362831e-04 [24,] 0.9990330 1.933975e-03 9.669873e-04 [25,] 0.9988027 2.394662e-03 1.197331e-03 [26,] 0.9982586 3.482897e-03 1.741449e-03 [27,] 0.9975478 4.904351e-03 2.452176e-03 [28,] 0.9985692 2.861645e-03 1.430822e-03 [29,] 0.9987660 2.467985e-03 1.233993e-03 [30,] 0.9983005 3.398958e-03 1.699479e-03 [31,] 0.9976672 4.665542e-03 2.332771e-03 [32,] 0.9975619 4.876225e-03 2.438112e-03 [33,] 0.9964846 7.030890e-03 3.515445e-03 [34,] 0.9958735 8.253051e-03 4.126525e-03 [35,] 0.9950342 9.931613e-03 4.965807e-03 [36,] 0.9967258 6.548342e-03 3.274171e-03 [37,] 0.9958564 8.287150e-03 4.143575e-03 [38,] 0.9943510 1.129800e-02 5.649001e-03 [39,] 0.9974769 5.046204e-03 2.523102e-03 [40,] 0.9998609 2.782776e-04 1.391388e-04 [41,] 0.9997920 4.159889e-04 2.079944e-04 [42,] 0.9997657 4.686139e-04 2.343070e-04 [43,] 0.9998531 2.937282e-04 1.468641e-04 [44,] 0.9998607 2.785973e-04 1.392986e-04 [45,] 0.9998016 3.967474e-04 1.983737e-04 [46,] 0.9997141 5.718778e-04 2.859389e-04 [47,] 0.9999529 9.422995e-05 4.711498e-05 [48,] 0.9999289 1.421343e-04 7.106713e-05 [49,] 0.9998931 2.138768e-04 1.069384e-04 [50,] 0.9998481 3.038168e-04 1.519084e-04 [51,] 0.9997793 4.414870e-04 2.207435e-04 [52,] 0.9998595 2.810431e-04 1.405216e-04 [53,] 0.9998038 3.923056e-04 1.961528e-04 [54,] 0.9997866 4.267453e-04 2.133727e-04 [55,] 0.9997214 5.571724e-04 2.785862e-04 [56,] 0.9997459 5.081415e-04 2.540707e-04 [57,] 0.9996424 7.151755e-04 3.575877e-04 [58,] 0.9995081 9.838168e-04 4.919084e-04 [59,] 0.9993089 1.382137e-03 6.910685e-04 [60,] 0.9991501 1.699775e-03 8.498877e-04 [61,] 0.9995896 8.208572e-04 4.104286e-04 [62,] 0.9995477 9.046319e-04 4.523159e-04 [63,] 0.9996263 7.473914e-04 3.736957e-04 [64,] 0.9996474 7.051547e-04 3.525773e-04 [65,] 0.9995450 9.100231e-04 4.550115e-04 [66,] 0.9995827 8.346358e-04 4.173179e-04 [67,] 0.9994288 1.142332e-03 5.711659e-04 [68,] 0.9992091 1.581849e-03 7.909244e-04 [69,] 0.9989326 2.134753e-03 1.067377e-03 [70,] 0.9987713 2.457411e-03 1.228706e-03 [71,] 0.9987509 2.498203e-03 1.249102e-03 [72,] 0.9996941 6.118130e-04 3.059065e-04 [73,] 0.9997663 4.674123e-04 2.337062e-04 [74,] 0.9997411 5.177935e-04 2.588967e-04 [75,] 0.9997490 5.019089e-04 2.509545e-04 [76,] 0.9996488 7.024251e-04 3.512126e-04 [77,] 0.9995535 8.930740e-04 4.465370e-04 [78,] 0.9998254 3.492482e-04 1.746241e-04 [79,] 0.9997824 4.352234e-04 2.176117e-04 [80,] 0.9997125 5.750584e-04 2.875292e-04 [81,] 0.9996078 7.844997e-04 3.922498e-04 [82,] 0.9994588 1.082489e-03 5.412447e-04 [83,] 0.9992682 1.463690e-03 7.318452e-04 [84,] 0.9990499 1.900273e-03 9.501365e-04 [85,] 0.9991119 1.776229e-03 8.881146e-04 [86,] 0.9993359 1.328119e-03 6.640594e-04 [87,] 0.9999996 8.626920e-07 4.313460e-07 [88,] 0.9999996 8.429039e-07 4.214520e-07 [89,] 0.9999994 1.109019e-06 5.545095e-07 [90,] 0.9999992 1.545155e-06 7.725777e-07 [91,] 0.9999990 1.964914e-06 9.824572e-07 [92,] 0.9999996 7.124137e-07 3.562068e-07 [93,] 0.9999994 1.114142e-06 5.570709e-07 [94,] 0.9999992 1.539787e-06 7.698937e-07 [95,] 0.9999990 2.042600e-06 1.021300e-06 [96,] 0.9999992 1.574929e-06 7.874643e-07 [97,] 0.9999990 2.051676e-06 1.025838e-06 [98,] 0.9999990 2.082021e-06 1.041010e-06 [99,] 0.9999984 3.172007e-06 1.586003e-06 [100,] 0.9999978 4.473893e-06 2.236946e-06 [101,] 0.9999967 6.506268e-06 3.253134e-06 [102,] 0.9999956 8.718464e-06 4.359232e-06 [103,] 0.9999938 1.246910e-05 6.234552e-06 [104,] 0.9999910 1.792110e-05 8.960551e-06 [105,] 0.9999870 2.596574e-05 1.298287e-05 [106,] 0.9999830 3.408323e-05 1.704162e-05 [107,] 0.9999888 2.234907e-05 1.117453e-05 [108,] 0.9999887 2.259257e-05 1.129628e-05 [109,] 0.9999885 2.294458e-05 1.147229e-05 [110,] 0.9999838 3.249072e-05 1.624536e-05 [111,] 0.9999871 2.573186e-05 1.286593e-05 [112,] 0.9999843 3.141846e-05 1.570923e-05 [113,] 0.9999976 4.778673e-06 2.389337e-06 [114,] 0.9999995 1.043523e-06 5.217616e-07 [115,] 0.9999993 1.416762e-06 7.083811e-07 [116,] 0.9999990 2.073775e-06 1.036887e-06 [117,] 0.9999991 1.759447e-06 8.797236e-07 [118,] 0.9999988 2.310344e-06 1.155172e-06 [119,] 0.9999983 3.494644e-06 1.747322e-06 [120,] 0.9999974 5.266982e-06 2.633491e-06 [121,] 0.9999977 4.626703e-06 2.313352e-06 [122,] 0.9999971 5.853882e-06 2.926941e-06 [123,] 0.9999997 6.952090e-07 3.476045e-07 [124,] 0.9999995 9.791189e-07 4.895595e-07 [125,] 1.0000000 7.289551e-08 3.644776e-08 [126,] 0.9999999 1.160952e-07 5.804759e-08 [127,] 1.0000000 3.933056e-08 1.966528e-08 [128,] 1.0000000 3.485228e-08 1.742614e-08 [129,] 1.0000000 5.633664e-08 2.816832e-08 [130,] 1.0000000 6.698680e-08 3.349340e-08 [131,] 1.0000000 6.009510e-08 3.004755e-08 [132,] 1.0000000 8.399503e-08 4.199752e-08 [133,] 0.9999999 1.268164e-07 6.340821e-08 [134,] 0.9999999 1.954728e-07 9.773639e-08 [135,] 0.9999999 1.602550e-07 8.012750e-08 [136,] 0.9999999 2.421748e-07 1.210874e-07 [137,] 0.9999998 3.018233e-07 1.509116e-07 [138,] 0.9999998 4.562874e-07 2.281437e-07 [139,] 0.9999997 6.957098e-07 3.478549e-07 [140,] 0.9999997 5.602047e-07 2.801024e-07 [141,] 0.9999997 6.202345e-07 3.101172e-07 [142,] 0.9999997 6.627639e-07 3.313820e-07 [143,] 0.9999999 2.894049e-07 1.447025e-07 [144,] 0.9999998 3.013428e-07 1.506714e-07 [145,] 0.9999999 2.879696e-07 1.439848e-07 [146,] 0.9999999 2.675885e-07 1.337942e-07 [147,] 0.9999999 2.702299e-07 1.351150e-07 [148,] 0.9999998 3.339899e-07 1.669950e-07 [149,] 0.9999997 5.338862e-07 2.669431e-07 [150,] 0.9999997 6.050862e-07 3.025431e-07 [151,] 0.9999996 8.648947e-07 4.324474e-07 [152,] 0.9999995 1.007450e-06 5.037248e-07 [153,] 0.9999992 1.520466e-06 7.602328e-07 [154,] 0.9999990 2.010977e-06 1.005488e-06 [155,] 0.9999995 9.632174e-07 4.816087e-07 [156,] 0.9999994 1.251632e-06 6.258162e-07 [157,] 0.9999993 1.452232e-06 7.261162e-07 [158,] 0.9999990 1.962459e-06 9.812295e-07 [159,] 1.0000000 6.857280e-08 3.428640e-08 [160,] 1.0000000 7.153306e-08 3.576653e-08 [161,] 0.9999999 1.004740e-07 5.023702e-08 [162,] 0.9999999 1.639423e-07 8.197115e-08 [163,] 1.0000000 7.391583e-08 3.695792e-08 [164,] 1.0000000 5.291284e-09 2.645642e-09 [165,] 1.0000000 6.609671e-09 3.304836e-09 [166,] 1.0000000 8.109397e-09 4.054698e-09 [167,] 1.0000000 1.346114e-08 6.730571e-09 [168,] 1.0000000 2.120307e-08 1.060154e-08 [169,] 1.0000000 8.352994e-09 4.176497e-09 [170,] 1.0000000 9.696424e-09 4.848212e-09 [171,] 1.0000000 6.381423e-10 3.190712e-10 [172,] 1.0000000 1.134803e-09 5.674017e-10 [173,] 1.0000000 8.634431e-10 4.317215e-10 [174,] 1.0000000 1.033885e-10 5.169427e-11 [175,] 1.0000000 1.752413e-10 8.762067e-11 [176,] 1.0000000 3.188824e-10 1.594412e-10 [177,] 1.0000000 5.757686e-10 2.878843e-10 [178,] 1.0000000 1.003423e-09 5.017114e-10 [179,] 1.0000000 1.482370e-09 7.411850e-10 [180,] 1.0000000 5.329343e-10 2.664671e-10 [181,] 1.0000000 8.454807e-10 4.227404e-10 [182,] 1.0000000 1.425859e-09 7.129293e-10 [183,] 1.0000000 2.407034e-09 1.203517e-09 [184,] 1.0000000 3.723410e-09 1.861705e-09 [185,] 1.0000000 5.503356e-09 2.751678e-09 [186,] 1.0000000 9.123764e-09 4.561882e-09 [187,] 1.0000000 1.556371e-08 7.781855e-09 [188,] 1.0000000 2.387216e-08 1.193608e-08 [189,] 1.0000000 4.237646e-08 2.118823e-08 [190,] 1.0000000 7.402676e-08 3.701338e-08 [191,] 0.9999999 1.276420e-07 6.382100e-08 [192,] 0.9999999 1.426127e-07 7.130637e-08 [193,] 0.9999999 2.358677e-07 1.179338e-07 [194,] 0.9999998 3.666916e-07 1.833458e-07 [195,] 0.9999997 6.026002e-07 3.013001e-07 [196,] 0.9999996 8.779492e-07 4.389746e-07 [197,] 1.0000000 2.889177e-08 1.444589e-08 [198,] 1.0000000 1.109052e-08 5.545262e-09 [199,] 1.0000000 1.132810e-08 5.664050e-09 [200,] 1.0000000 1.138417e-08 5.692084e-09 [201,] 1.0000000 5.360419e-09 2.680209e-09 [202,] 1.0000000 1.043743e-08 5.218714e-09 [203,] 1.0000000 1.551091e-08 7.755454e-09 [204,] 1.0000000 2.396202e-08 1.198101e-08 [205,] 1.0000000 3.809559e-08 1.904780e-08 [206,] 1.0000000 4.551524e-08 2.275762e-08 [207,] 1.0000000 8.436733e-08 4.218367e-08 [208,] 0.9999999 1.084483e-07 5.422415e-08 [209,] 1.0000000 6.103934e-08 3.051967e-08 [210,] 0.9999999 1.135204e-07 5.676019e-08 [211,] 0.9999999 1.918536e-07 9.592681e-08 [212,] 1.0000000 9.082595e-08 4.541297e-08 [213,] 0.9999999 1.575523e-07 7.877615e-08 [214,] 0.9999999 2.927092e-07 1.463546e-07 [215,] 0.9999998 3.425875e-07 1.712938e-07 [216,] 0.9999999 2.978407e-07 1.489204e-07 [217,] 0.9999997 5.062589e-07 2.531294e-07 [218,] 0.9999996 7.100118e-07 3.550059e-07 [219,] 1.0000000 7.400682e-08 3.700341e-08 [220,] 1.0000000 7.932638e-08 3.966319e-08 [221,] 0.9999999 1.550363e-07 7.751813e-08 [222,] 1.0000000 7.501654e-08 3.750827e-08 [223,] 0.9999999 1.274761e-07 6.373807e-08 [224,] 0.9999999 1.674849e-07 8.374247e-08 [225,] 0.9999999 2.909156e-07 1.454578e-07 [226,] 1.0000000 1.191653e-09 5.958263e-10 [227,] 1.0000000 2.828302e-09 1.414151e-09 [228,] 1.0000000 6.159047e-09 3.079524e-09 [229,] 1.0000000 6.812864e-09 3.406432e-09 [230,] 1.0000000 1.232355e-08 6.161773e-09 [231,] 1.0000000 2.698261e-08 1.349131e-08 [232,] 1.0000000 2.874906e-08 1.437453e-08 [233,] 1.0000000 5.809263e-08 2.904631e-08 [234,] 0.9999999 1.340533e-07 6.702665e-08 [235,] 0.9999998 3.034887e-07 1.517444e-07 [236,] 0.9999998 3.094724e-07 1.547362e-07 [237,] 0.9999998 3.819331e-07 1.909666e-07 [238,] 0.9999996 8.078442e-07 4.039221e-07 [239,] 0.9999992 1.649223e-06 8.246115e-07 [240,] 0.9999991 1.894189e-06 9.470945e-07 [241,] 0.9999979 4.112367e-06 2.056183e-06 [242,] 0.9999954 9.179102e-06 4.589551e-06 [243,] 0.9999899 2.027985e-05 1.013992e-05 [244,] 0.9999779 4.417161e-05 2.208581e-05 [245,] 0.9999531 9.379897e-05 4.689949e-05 [246,] 0.9999403 1.193483e-04 5.967417e-05 [247,] 0.9999150 1.700955e-04 8.504774e-05 [248,] 0.9998470 3.060888e-04 1.530444e-04 [249,] 0.9996777 6.445862e-04 3.222931e-04 [250,] 0.9994319 1.136196e-03 5.680981e-04 [251,] 0.9990279 1.944137e-03 9.720683e-04 [252,] 0.9980765 3.846987e-03 1.923494e-03 [253,] 0.9971994 5.601195e-03 2.800598e-03 [254,] 0.9974298 5.140481e-03 2.570241e-03 [255,] 0.9949007 1.019857e-02 5.099283e-03 [256,] 0.9970569 5.886180e-03 2.943090e-03 [257,] 0.9945902 1.081958e-02 5.409790e-03 [258,] 0.9994968 1.006438e-03 5.032188e-04 [259,] 0.9991689 1.662265e-03 8.311323e-04 [260,] 0.9978952 4.209536e-03 2.104768e-03 [261,] 0.9948092 1.038153e-02 5.190763e-03 [262,] 0.9911931 1.761377e-02 8.806885e-03 [263,] 0.9790116 4.197670e-02 2.098835e-02 [264,] 0.9497541 1.004917e-01 5.024587e-02 [265,] 0.8925233 2.149535e-01 1.074767e-01 [266,] 0.8759831 2.480338e-01 1.240169e-01 > postscript(file="/var/wessaorg/rcomp/tmp/1de2p1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/264xe1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/3q5ap1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/4795c1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5juzk1324466480.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 39311.89863 -4972.98917 21759.92248 -52650.25936 22214.18449 6 7 8 9 10 -12884.05949 70852.89528 -2788.76485 13.24984 -19484.53103 11 12 13 14 15 16685.95385 36559.53426 -13674.68602 -6050.96424 28122.54285 16 17 18 19 20 34489.00676 -79472.81624 20145.41417 18289.56003 -6256.57013 21 22 23 24 25 -21058.69191 -4274.86990 120426.53329 24125.90655 -40744.09666 26 27 28 29 30 -83134.99463 -54083.51388 -5448.57143 22651.63412 -4150.84535 31 32 33 34 35 21649.41276 486.78876 894.87108 35376.51515 -23907.22038 36 37 38 39 40 34001.26317 18239.07291 24240.14611 57128.89000 -6589.53964 41 42 43 44 45 15298.32456 21522.02447 -48798.04363 -2353.94603 18841.98645 46 47 48 49 50 8664.97771 -57514.33608 3585.19881 -7427.40490 -58238.70703 51 52 53 54 55 56723.59676 -12347.54869 -30814.14673 -44831.05397 -36361.26423 56 57 58 59 60 -3334.23601 -14001.00292 67027.86879 -8387.45094 -837.26567 61 62 63 64 65 -741.83980 -9277.54637 -52174.27831 5068.19553 -27142.39072 66 67 68 69 70 16152.25866 33741.72169 -5526.75993 -2202.30423 -1827.58717 71 72 73 74 75 11877.50789 50769.15196 25162.84293 -35611.31504 -34844.72236 76 77 78 79 80 -16118.33614 -35313.21680 -11408.28481 -7125.58562 2313.89018 81 82 83 84 85 17955.01964 -34923.67530 65832.39974 38350.94122 -23937.45809 86 87 88 89 90 -34624.47622 -4223.08798 -18792.41659 54186.68826 15410.96423 91 92 93 94 95 11077.39931 -10585.54936 -5276.53501 -6152.71367 8858.39833 96 97 98 99 100 26403.29203 39628.97342 -134525.99604 28380.69994 -21333.35190 101 102 103 104 105 11155.85303 -17397.07988 47649.56301 -4298.99524 10903.14281 106 107 108 109 110 13849.58930 -35899.70759 -21135.24029 27075.53204 2168.30495 111 112 113 114 115 9347.35468 6262.56626 -17921.48319 7451.54542 -10052.69102 116 117 118 119 120 -6192.99932 -8620.75945 39220.86161 -35334.19029 -34391.25173 121 122 123 124 125 -11033.35493 -38405.16575 -28518.05125 -71619.23916 58721.09947 126 127 128 129 130 12457.08959 5053.40988 -35560.79787 16744.00235 -1311.85449 131 132 133 134 135 -7640.92877 -27218.89093 15347.24417 -69714.11314 -11452.95326 136 137 138 139 140 75689.71901 -1658.10440 45664.29384 -33566.25443 -3990.35253 141 142 143 144 145 -19735.06652 -28097.47854 12964.00702 -8126.90361 3514.30084 146 147 148 149 150 -33995.14486 -2888.06724 -20716.72490 -5462.93781 -4257.33307 151 152 153 154 155 35784.21577 24488.90835 -17546.08380 44968.81845 24262.52367 156 157 158 159 160 20366.97621 27403.10064 22558.45456 -19132.46386 -3943.91461 161 162 163 164 165 19181.04977 -20268.57120 -16563.36847 3586.21046 -16717.60936 166 167 168 169 170 36368.12609 5953.73709 24156.52882 17633.01975 72415.25487 171 172 173 174 175 13047.05478 2197.08783 -283.91364 18825.45159 -84449.68761 176 177 178 179 180 15108.76463 -33521.45373 8517.00026 -11975.80454 36540.21971 181 182 183 184 185 -26634.41095 -60739.68987 5736.37107 26076.02674 42029.49517 186 187 188 189 190 3031.21024 12935.80988 14756.88250 -1441.99545 14081.31360 191 192 193 194 195 27666.71469 11684.07754 -19327.53844 -20696.39693 -472.17236 196 197 198 199 200 1711.10049 -9659.48215 704.98200 -23841.08339 -8307.29664 201 202 203 204 205 5306.01056 9334.50984 14509.62400 347.77189 -21644.91521 206 207 208 209 210 -3546.99712 -15567.82447 50055.30204 26473.38992 7883.18184 211 212 213 214 215 9416.53186 12922.97848 -1066.52060 4805.98484 -24358.30973 216 217 218 219 220 -10695.15857 -29515.63779 7044.81888 10583.11348 27676.75548 221 222 223 224 225 -18468.72380 -32024.12292 -55373.60902 340.87940 8375.18731 226 227 228 229 230 5355.00092 17200.42427 -21838.24561 -26835.06367 -59478.25015 231 232 233 234 235 -10973.46033 -13746.04441 -36725.50802 -4361.65297 -17910.76376 236 237 238 239 240 -12694.87786 53513.99412 1433.55821 -4270.69738 -13745.65009 241 242 243 244 245 8421.71025 -19633.49856 -14033.98363 -3228.80810 1333.17937 246 247 248 249 250 -3696.73029 -34147.18670 24461.88553 9578.34396 -8861.25266 251 252 253 254 255 17155.42935 5710.10336 12385.22920 3518.42688 -6792.00279 256 257 258 259 260 8876.56821 -12091.52252 17109.94531 15971.82214 6750.99521 261 262 263 264 265 -1997.38968 1454.61676 9172.29191 -4530.12013 -12345.74724 266 267 268 269 270 12095.46304 -13771.19778 4738.22989 52727.10689 -9751.74114 271 272 273 274 275 -3960.24266 7610.07914 5720.07428 7204.32321 9754.30241 276 277 278 279 280 -9641.70626 11510.80817 10951.10329 -4098.58943 21416.05188 281 282 283 284 285 1445.36551 8129.50701 7031.44229 16389.03335 10773.53438 286 287 288 289 -22401.63529 1895.28429 -5398.14417 8137.64420 > postscript(file="/var/wessaorg/rcomp/tmp/643dw1324466480.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 39311.89863 NA 1 -4972.98917 39311.89863 2 21759.92248 -4972.98917 3 -52650.25936 21759.92248 4 22214.18449 -52650.25936 5 -12884.05949 22214.18449 6 70852.89528 -12884.05949 7 -2788.76485 70852.89528 8 13.24984 -2788.76485 9 -19484.53103 13.24984 10 16685.95385 -19484.53103 11 36559.53426 16685.95385 12 -13674.68602 36559.53426 13 -6050.96424 -13674.68602 14 28122.54285 -6050.96424 15 34489.00676 28122.54285 16 -79472.81624 34489.00676 17 20145.41417 -79472.81624 18 18289.56003 20145.41417 19 -6256.57013 18289.56003 20 -21058.69191 -6256.57013 21 -4274.86990 -21058.69191 22 120426.53329 -4274.86990 23 24125.90655 120426.53329 24 -40744.09666 24125.90655 25 -83134.99463 -40744.09666 26 -54083.51388 -83134.99463 27 -5448.57143 -54083.51388 28 22651.63412 -5448.57143 29 -4150.84535 22651.63412 30 21649.41276 -4150.84535 31 486.78876 21649.41276 32 894.87108 486.78876 33 35376.51515 894.87108 34 -23907.22038 35376.51515 35 34001.26317 -23907.22038 36 18239.07291 34001.26317 37 24240.14611 18239.07291 38 57128.89000 24240.14611 39 -6589.53964 57128.89000 40 15298.32456 -6589.53964 41 21522.02447 15298.32456 42 -48798.04363 21522.02447 43 -2353.94603 -48798.04363 44 18841.98645 -2353.94603 45 8664.97771 18841.98645 46 -57514.33608 8664.97771 47 3585.19881 -57514.33608 48 -7427.40490 3585.19881 49 -58238.70703 -7427.40490 50 56723.59676 -58238.70703 51 -12347.54869 56723.59676 52 -30814.14673 -12347.54869 53 -44831.05397 -30814.14673 54 -36361.26423 -44831.05397 55 -3334.23601 -36361.26423 56 -14001.00292 -3334.23601 57 67027.86879 -14001.00292 58 -8387.45094 67027.86879 59 -837.26567 -8387.45094 60 -741.83980 -837.26567 61 -9277.54637 -741.83980 62 -52174.27831 -9277.54637 63 5068.19553 -52174.27831 64 -27142.39072 5068.19553 65 16152.25866 -27142.39072 66 33741.72169 16152.25866 67 -5526.75993 33741.72169 68 -2202.30423 -5526.75993 69 -1827.58717 -2202.30423 70 11877.50789 -1827.58717 71 50769.15196 11877.50789 72 25162.84293 50769.15196 73 -35611.31504 25162.84293 74 -34844.72236 -35611.31504 75 -16118.33614 -34844.72236 76 -35313.21680 -16118.33614 77 -11408.28481 -35313.21680 78 -7125.58562 -11408.28481 79 2313.89018 -7125.58562 80 17955.01964 2313.89018 81 -34923.67530 17955.01964 82 65832.39974 -34923.67530 83 38350.94122 65832.39974 84 -23937.45809 38350.94122 85 -34624.47622 -23937.45809 86 -4223.08798 -34624.47622 87 -18792.41659 -4223.08798 88 54186.68826 -18792.41659 89 15410.96423 54186.68826 90 11077.39931 15410.96423 91 -10585.54936 11077.39931 92 -5276.53501 -10585.54936 93 -6152.71367 -5276.53501 94 8858.39833 -6152.71367 95 26403.29203 8858.39833 96 39628.97342 26403.29203 97 -134525.99604 39628.97342 98 28380.69994 -134525.99604 99 -21333.35190 28380.69994 100 11155.85303 -21333.35190 101 -17397.07988 11155.85303 102 47649.56301 -17397.07988 103 -4298.99524 47649.56301 104 10903.14281 -4298.99524 105 13849.58930 10903.14281 106 -35899.70759 13849.58930 107 -21135.24029 -35899.70759 108 27075.53204 -21135.24029 109 2168.30495 27075.53204 110 9347.35468 2168.30495 111 6262.56626 9347.35468 112 -17921.48319 6262.56626 113 7451.54542 -17921.48319 114 -10052.69102 7451.54542 115 -6192.99932 -10052.69102 116 -8620.75945 -6192.99932 117 39220.86161 -8620.75945 118 -35334.19029 39220.86161 119 -34391.25173 -35334.19029 120 -11033.35493 -34391.25173 121 -38405.16575 -11033.35493 122 -28518.05125 -38405.16575 123 -71619.23916 -28518.05125 124 58721.09947 -71619.23916 125 12457.08959 58721.09947 126 5053.40988 12457.08959 127 -35560.79787 5053.40988 128 16744.00235 -35560.79787 129 -1311.85449 16744.00235 130 -7640.92877 -1311.85449 131 -27218.89093 -7640.92877 132 15347.24417 -27218.89093 133 -69714.11314 15347.24417 134 -11452.95326 -69714.11314 135 75689.71901 -11452.95326 136 -1658.10440 75689.71901 137 45664.29384 -1658.10440 138 -33566.25443 45664.29384 139 -3990.35253 -33566.25443 140 -19735.06652 -3990.35253 141 -28097.47854 -19735.06652 142 12964.00702 -28097.47854 143 -8126.90361 12964.00702 144 3514.30084 -8126.90361 145 -33995.14486 3514.30084 146 -2888.06724 -33995.14486 147 -20716.72490 -2888.06724 148 -5462.93781 -20716.72490 149 -4257.33307 -5462.93781 150 35784.21577 -4257.33307 151 24488.90835 35784.21577 152 -17546.08380 24488.90835 153 44968.81845 -17546.08380 154 24262.52367 44968.81845 155 20366.97621 24262.52367 156 27403.10064 20366.97621 157 22558.45456 27403.10064 158 -19132.46386 22558.45456 159 -3943.91461 -19132.46386 160 19181.04977 -3943.91461 161 -20268.57120 19181.04977 162 -16563.36847 -20268.57120 163 3586.21046 -16563.36847 164 -16717.60936 3586.21046 165 36368.12609 -16717.60936 166 5953.73709 36368.12609 167 24156.52882 5953.73709 168 17633.01975 24156.52882 169 72415.25487 17633.01975 170 13047.05478 72415.25487 171 2197.08783 13047.05478 172 -283.91364 2197.08783 173 18825.45159 -283.91364 174 -84449.68761 18825.45159 175 15108.76463 -84449.68761 176 -33521.45373 15108.76463 177 8517.00026 -33521.45373 178 -11975.80454 8517.00026 179 36540.21971 -11975.80454 180 -26634.41095 36540.21971 181 -60739.68987 -26634.41095 182 5736.37107 -60739.68987 183 26076.02674 5736.37107 184 42029.49517 26076.02674 185 3031.21024 42029.49517 186 12935.80988 3031.21024 187 14756.88250 12935.80988 188 -1441.99545 14756.88250 189 14081.31360 -1441.99545 190 27666.71469 14081.31360 191 11684.07754 27666.71469 192 -19327.53844 11684.07754 193 -20696.39693 -19327.53844 194 -472.17236 -20696.39693 195 1711.10049 -472.17236 196 -9659.48215 1711.10049 197 704.98200 -9659.48215 198 -23841.08339 704.98200 199 -8307.29664 -23841.08339 200 5306.01056 -8307.29664 201 9334.50984 5306.01056 202 14509.62400 9334.50984 203 347.77189 14509.62400 204 -21644.91521 347.77189 205 -3546.99712 -21644.91521 206 -15567.82447 -3546.99712 207 50055.30204 -15567.82447 208 26473.38992 50055.30204 209 7883.18184 26473.38992 210 9416.53186 7883.18184 211 12922.97848 9416.53186 212 -1066.52060 12922.97848 213 4805.98484 -1066.52060 214 -24358.30973 4805.98484 215 -10695.15857 -24358.30973 216 -29515.63779 -10695.15857 217 7044.81888 -29515.63779 218 10583.11348 7044.81888 219 27676.75548 10583.11348 220 -18468.72380 27676.75548 221 -32024.12292 -18468.72380 222 -55373.60902 -32024.12292 223 340.87940 -55373.60902 224 8375.18731 340.87940 225 5355.00092 8375.18731 226 17200.42427 5355.00092 227 -21838.24561 17200.42427 228 -26835.06367 -21838.24561 229 -59478.25015 -26835.06367 230 -10973.46033 -59478.25015 231 -13746.04441 -10973.46033 232 -36725.50802 -13746.04441 233 -4361.65297 -36725.50802 234 -17910.76376 -4361.65297 235 -12694.87786 -17910.76376 236 53513.99412 -12694.87786 237 1433.55821 53513.99412 238 -4270.69738 1433.55821 239 -13745.65009 -4270.69738 240 8421.71025 -13745.65009 241 -19633.49856 8421.71025 242 -14033.98363 -19633.49856 243 -3228.80810 -14033.98363 244 1333.17937 -3228.80810 245 -3696.73029 1333.17937 246 -34147.18670 -3696.73029 247 24461.88553 -34147.18670 248 9578.34396 24461.88553 249 -8861.25266 9578.34396 250 17155.42935 -8861.25266 251 5710.10336 17155.42935 252 12385.22920 5710.10336 253 3518.42688 12385.22920 254 -6792.00279 3518.42688 255 8876.56821 -6792.00279 256 -12091.52252 8876.56821 257 17109.94531 -12091.52252 258 15971.82214 17109.94531 259 6750.99521 15971.82214 260 -1997.38968 6750.99521 261 1454.61676 -1997.38968 262 9172.29191 1454.61676 263 -4530.12013 9172.29191 264 -12345.74724 -4530.12013 265 12095.46304 -12345.74724 266 -13771.19778 12095.46304 267 4738.22989 -13771.19778 268 52727.10689 4738.22989 269 -9751.74114 52727.10689 270 -3960.24266 -9751.74114 271 7610.07914 -3960.24266 272 5720.07428 7610.07914 273 7204.32321 5720.07428 274 9754.30241 7204.32321 275 -9641.70626 9754.30241 276 11510.80817 -9641.70626 277 10951.10329 11510.80817 278 -4098.58943 10951.10329 279 21416.05188 -4098.58943 280 1445.36551 21416.05188 281 8129.50701 1445.36551 282 7031.44229 8129.50701 283 16389.03335 7031.44229 284 10773.53438 16389.03335 285 -22401.63529 10773.53438 286 1895.28429 -22401.63529 287 -5398.14417 1895.28429 288 8137.64420 -5398.14417 289 NA 8137.64420 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -4972.98917 39311.89863 [2,] 21759.92248 -4972.98917 [3,] -52650.25936 21759.92248 [4,] 22214.18449 -52650.25936 [5,] -12884.05949 22214.18449 [6,] 70852.89528 -12884.05949 [7,] -2788.76485 70852.89528 [8,] 13.24984 -2788.76485 [9,] -19484.53103 13.24984 [10,] 16685.95385 -19484.53103 [11,] 36559.53426 16685.95385 [12,] -13674.68602 36559.53426 [13,] -6050.96424 -13674.68602 [14,] 28122.54285 -6050.96424 [15,] 34489.00676 28122.54285 [16,] -79472.81624 34489.00676 [17,] 20145.41417 -79472.81624 [18,] 18289.56003 20145.41417 [19,] -6256.57013 18289.56003 [20,] -21058.69191 -6256.57013 [21,] -4274.86990 -21058.69191 [22,] 120426.53329 -4274.86990 [23,] 24125.90655 120426.53329 [24,] -40744.09666 24125.90655 [25,] -83134.99463 -40744.09666 [26,] -54083.51388 -83134.99463 [27,] -5448.57143 -54083.51388 [28,] 22651.63412 -5448.57143 [29,] -4150.84535 22651.63412 [30,] 21649.41276 -4150.84535 [31,] 486.78876 21649.41276 [32,] 894.87108 486.78876 [33,] 35376.51515 894.87108 [34,] -23907.22038 35376.51515 [35,] 34001.26317 -23907.22038 [36,] 18239.07291 34001.26317 [37,] 24240.14611 18239.07291 [38,] 57128.89000 24240.14611 [39,] -6589.53964 57128.89000 [40,] 15298.32456 -6589.53964 [41,] 21522.02447 15298.32456 [42,] -48798.04363 21522.02447 [43,] -2353.94603 -48798.04363 [44,] 18841.98645 -2353.94603 [45,] 8664.97771 18841.98645 [46,] -57514.33608 8664.97771 [47,] 3585.19881 -57514.33608 [48,] -7427.40490 3585.19881 [49,] -58238.70703 -7427.40490 [50,] 56723.59676 -58238.70703 [51,] -12347.54869 56723.59676 [52,] -30814.14673 -12347.54869 [53,] -44831.05397 -30814.14673 [54,] -36361.26423 -44831.05397 [55,] -3334.23601 -36361.26423 [56,] -14001.00292 -3334.23601 [57,] 67027.86879 -14001.00292 [58,] -8387.45094 67027.86879 [59,] -837.26567 -8387.45094 [60,] -741.83980 -837.26567 [61,] -9277.54637 -741.83980 [62,] -52174.27831 -9277.54637 [63,] 5068.19553 -52174.27831 [64,] -27142.39072 5068.19553 [65,] 16152.25866 -27142.39072 [66,] 33741.72169 16152.25866 [67,] -5526.75993 33741.72169 [68,] -2202.30423 -5526.75993 [69,] -1827.58717 -2202.30423 [70,] 11877.50789 -1827.58717 [71,] 50769.15196 11877.50789 [72,] 25162.84293 50769.15196 [73,] -35611.31504 25162.84293 [74,] -34844.72236 -35611.31504 [75,] -16118.33614 -34844.72236 [76,] -35313.21680 -16118.33614 [77,] -11408.28481 -35313.21680 [78,] -7125.58562 -11408.28481 [79,] 2313.89018 -7125.58562 [80,] 17955.01964 2313.89018 [81,] -34923.67530 17955.01964 [82,] 65832.39974 -34923.67530 [83,] 38350.94122 65832.39974 [84,] -23937.45809 38350.94122 [85,] -34624.47622 -23937.45809 [86,] -4223.08798 -34624.47622 [87,] -18792.41659 -4223.08798 [88,] 54186.68826 -18792.41659 [89,] 15410.96423 54186.68826 [90,] 11077.39931 15410.96423 [91,] -10585.54936 11077.39931 [92,] -5276.53501 -10585.54936 [93,] -6152.71367 -5276.53501 [94,] 8858.39833 -6152.71367 [95,] 26403.29203 8858.39833 [96,] 39628.97342 26403.29203 [97,] -134525.99604 39628.97342 [98,] 28380.69994 -134525.99604 [99,] -21333.35190 28380.69994 [100,] 11155.85303 -21333.35190 [101,] -17397.07988 11155.85303 [102,] 47649.56301 -17397.07988 [103,] -4298.99524 47649.56301 [104,] 10903.14281 -4298.99524 [105,] 13849.58930 10903.14281 [106,] -35899.70759 13849.58930 [107,] -21135.24029 -35899.70759 [108,] 27075.53204 -21135.24029 [109,] 2168.30495 27075.53204 [110,] 9347.35468 2168.30495 [111,] 6262.56626 9347.35468 [112,] -17921.48319 6262.56626 [113,] 7451.54542 -17921.48319 [114,] -10052.69102 7451.54542 [115,] -6192.99932 -10052.69102 [116,] -8620.75945 -6192.99932 [117,] 39220.86161 -8620.75945 [118,] -35334.19029 39220.86161 [119,] -34391.25173 -35334.19029 [120,] -11033.35493 -34391.25173 [121,] -38405.16575 -11033.35493 [122,] -28518.05125 -38405.16575 [123,] -71619.23916 -28518.05125 [124,] 58721.09947 -71619.23916 [125,] 12457.08959 58721.09947 [126,] 5053.40988 12457.08959 [127,] -35560.79787 5053.40988 [128,] 16744.00235 -35560.79787 [129,] -1311.85449 16744.00235 [130,] -7640.92877 -1311.85449 [131,] -27218.89093 -7640.92877 [132,] 15347.24417 -27218.89093 [133,] -69714.11314 15347.24417 [134,] -11452.95326 -69714.11314 [135,] 75689.71901 -11452.95326 [136,] -1658.10440 75689.71901 [137,] 45664.29384 -1658.10440 [138,] -33566.25443 45664.29384 [139,] -3990.35253 -33566.25443 [140,] -19735.06652 -3990.35253 [141,] -28097.47854 -19735.06652 [142,] 12964.00702 -28097.47854 [143,] -8126.90361 12964.00702 [144,] 3514.30084 -8126.90361 [145,] -33995.14486 3514.30084 [146,] -2888.06724 -33995.14486 [147,] -20716.72490 -2888.06724 [148,] -5462.93781 -20716.72490 [149,] -4257.33307 -5462.93781 [150,] 35784.21577 -4257.33307 [151,] 24488.90835 35784.21577 [152,] -17546.08380 24488.90835 [153,] 44968.81845 -17546.08380 [154,] 24262.52367 44968.81845 [155,] 20366.97621 24262.52367 [156,] 27403.10064 20366.97621 [157,] 22558.45456 27403.10064 [158,] -19132.46386 22558.45456 [159,] -3943.91461 -19132.46386 [160,] 19181.04977 -3943.91461 [161,] -20268.57120 19181.04977 [162,] -16563.36847 -20268.57120 [163,] 3586.21046 -16563.36847 [164,] -16717.60936 3586.21046 [165,] 36368.12609 -16717.60936 [166,] 5953.73709 36368.12609 [167,] 24156.52882 5953.73709 [168,] 17633.01975 24156.52882 [169,] 72415.25487 17633.01975 [170,] 13047.05478 72415.25487 [171,] 2197.08783 13047.05478 [172,] -283.91364 2197.08783 [173,] 18825.45159 -283.91364 [174,] -84449.68761 18825.45159 [175,] 15108.76463 -84449.68761 [176,] -33521.45373 15108.76463 [177,] 8517.00026 -33521.45373 [178,] -11975.80454 8517.00026 [179,] 36540.21971 -11975.80454 [180,] -26634.41095 36540.21971 [181,] -60739.68987 -26634.41095 [182,] 5736.37107 -60739.68987 [183,] 26076.02674 5736.37107 [184,] 42029.49517 26076.02674 [185,] 3031.21024 42029.49517 [186,] 12935.80988 3031.21024 [187,] 14756.88250 12935.80988 [188,] -1441.99545 14756.88250 [189,] 14081.31360 -1441.99545 [190,] 27666.71469 14081.31360 [191,] 11684.07754 27666.71469 [192,] -19327.53844 11684.07754 [193,] -20696.39693 -19327.53844 [194,] -472.17236 -20696.39693 [195,] 1711.10049 -472.17236 [196,] -9659.48215 1711.10049 [197,] 704.98200 -9659.48215 [198,] -23841.08339 704.98200 [199,] -8307.29664 -23841.08339 [200,] 5306.01056 -8307.29664 [201,] 9334.50984 5306.01056 [202,] 14509.62400 9334.50984 [203,] 347.77189 14509.62400 [204,] -21644.91521 347.77189 [205,] -3546.99712 -21644.91521 [206,] -15567.82447 -3546.99712 [207,] 50055.30204 -15567.82447 [208,] 26473.38992 50055.30204 [209,] 7883.18184 26473.38992 [210,] 9416.53186 7883.18184 [211,] 12922.97848 9416.53186 [212,] -1066.52060 12922.97848 [213,] 4805.98484 -1066.52060 [214,] -24358.30973 4805.98484 [215,] -10695.15857 -24358.30973 [216,] -29515.63779 -10695.15857 [217,] 7044.81888 -29515.63779 [218,] 10583.11348 7044.81888 [219,] 27676.75548 10583.11348 [220,] -18468.72380 27676.75548 [221,] -32024.12292 -18468.72380 [222,] -55373.60902 -32024.12292 [223,] 340.87940 -55373.60902 [224,] 8375.18731 340.87940 [225,] 5355.00092 8375.18731 [226,] 17200.42427 5355.00092 [227,] -21838.24561 17200.42427 [228,] -26835.06367 -21838.24561 [229,] -59478.25015 -26835.06367 [230,] -10973.46033 -59478.25015 [231,] -13746.04441 -10973.46033 [232,] -36725.50802 -13746.04441 [233,] -4361.65297 -36725.50802 [234,] -17910.76376 -4361.65297 [235,] -12694.87786 -17910.76376 [236,] 53513.99412 -12694.87786 [237,] 1433.55821 53513.99412 [238,] -4270.69738 1433.55821 [239,] -13745.65009 -4270.69738 [240,] 8421.71025 -13745.65009 [241,] -19633.49856 8421.71025 [242,] -14033.98363 -19633.49856 [243,] -3228.80810 -14033.98363 [244,] 1333.17937 -3228.80810 [245,] -3696.73029 1333.17937 [246,] -34147.18670 -3696.73029 [247,] 24461.88553 -34147.18670 [248,] 9578.34396 24461.88553 [249,] -8861.25266 9578.34396 [250,] 17155.42935 -8861.25266 [251,] 5710.10336 17155.42935 [252,] 12385.22920 5710.10336 [253,] 3518.42688 12385.22920 [254,] -6792.00279 3518.42688 [255,] 8876.56821 -6792.00279 [256,] -12091.52252 8876.56821 [257,] 17109.94531 -12091.52252 [258,] 15971.82214 17109.94531 [259,] 6750.99521 15971.82214 [260,] -1997.38968 6750.99521 [261,] 1454.61676 -1997.38968 [262,] 9172.29191 1454.61676 [263,] -4530.12013 9172.29191 [264,] -12345.74724 -4530.12013 [265,] 12095.46304 -12345.74724 [266,] -13771.19778 12095.46304 [267,] 4738.22989 -13771.19778 [268,] 52727.10689 4738.22989 [269,] -9751.74114 52727.10689 [270,] -3960.24266 -9751.74114 [271,] 7610.07914 -3960.24266 [272,] 5720.07428 7610.07914 [273,] 7204.32321 5720.07428 [274,] 9754.30241 7204.32321 [275,] -9641.70626 9754.30241 [276,] 11510.80817 -9641.70626 [277,] 10951.10329 11510.80817 [278,] -4098.58943 10951.10329 [279,] 21416.05188 -4098.58943 [280,] 1445.36551 21416.05188 [281,] 8129.50701 1445.36551 [282,] 7031.44229 8129.50701 [283,] 16389.03335 7031.44229 [284,] 10773.53438 16389.03335 [285,] -22401.63529 10773.53438 [286,] 1895.28429 -22401.63529 [287,] -5398.14417 1895.28429 [288,] 8137.64420 -5398.14417 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -4972.98917 39311.89863 2 21759.92248 -4972.98917 3 -52650.25936 21759.92248 4 22214.18449 -52650.25936 5 -12884.05949 22214.18449 6 70852.89528 -12884.05949 7 -2788.76485 70852.89528 8 13.24984 -2788.76485 9 -19484.53103 13.24984 10 16685.95385 -19484.53103 11 36559.53426 16685.95385 12 -13674.68602 36559.53426 13 -6050.96424 -13674.68602 14 28122.54285 -6050.96424 15 34489.00676 28122.54285 16 -79472.81624 34489.00676 17 20145.41417 -79472.81624 18 18289.56003 20145.41417 19 -6256.57013 18289.56003 20 -21058.69191 -6256.57013 21 -4274.86990 -21058.69191 22 120426.53329 -4274.86990 23 24125.90655 120426.53329 24 -40744.09666 24125.90655 25 -83134.99463 -40744.09666 26 -54083.51388 -83134.99463 27 -5448.57143 -54083.51388 28 22651.63412 -5448.57143 29 -4150.84535 22651.63412 30 21649.41276 -4150.84535 31 486.78876 21649.41276 32 894.87108 486.78876 33 35376.51515 894.87108 34 -23907.22038 35376.51515 35 34001.26317 -23907.22038 36 18239.07291 34001.26317 37 24240.14611 18239.07291 38 57128.89000 24240.14611 39 -6589.53964 57128.89000 40 15298.32456 -6589.53964 41 21522.02447 15298.32456 42 -48798.04363 21522.02447 43 -2353.94603 -48798.04363 44 18841.98645 -2353.94603 45 8664.97771 18841.98645 46 -57514.33608 8664.97771 47 3585.19881 -57514.33608 48 -7427.40490 3585.19881 49 -58238.70703 -7427.40490 50 56723.59676 -58238.70703 51 -12347.54869 56723.59676 52 -30814.14673 -12347.54869 53 -44831.05397 -30814.14673 54 -36361.26423 -44831.05397 55 -3334.23601 -36361.26423 56 -14001.00292 -3334.23601 57 67027.86879 -14001.00292 58 -8387.45094 67027.86879 59 -837.26567 -8387.45094 60 -741.83980 -837.26567 61 -9277.54637 -741.83980 62 -52174.27831 -9277.54637 63 5068.19553 -52174.27831 64 -27142.39072 5068.19553 65 16152.25866 -27142.39072 66 33741.72169 16152.25866 67 -5526.75993 33741.72169 68 -2202.30423 -5526.75993 69 -1827.58717 -2202.30423 70 11877.50789 -1827.58717 71 50769.15196 11877.50789 72 25162.84293 50769.15196 73 -35611.31504 25162.84293 74 -34844.72236 -35611.31504 75 -16118.33614 -34844.72236 76 -35313.21680 -16118.33614 77 -11408.28481 -35313.21680 78 -7125.58562 -11408.28481 79 2313.89018 -7125.58562 80 17955.01964 2313.89018 81 -34923.67530 17955.01964 82 65832.39974 -34923.67530 83 38350.94122 65832.39974 84 -23937.45809 38350.94122 85 -34624.47622 -23937.45809 86 -4223.08798 -34624.47622 87 -18792.41659 -4223.08798 88 54186.68826 -18792.41659 89 15410.96423 54186.68826 90 11077.39931 15410.96423 91 -10585.54936 11077.39931 92 -5276.53501 -10585.54936 93 -6152.71367 -5276.53501 94 8858.39833 -6152.71367 95 26403.29203 8858.39833 96 39628.97342 26403.29203 97 -134525.99604 39628.97342 98 28380.69994 -134525.99604 99 -21333.35190 28380.69994 100 11155.85303 -21333.35190 101 -17397.07988 11155.85303 102 47649.56301 -17397.07988 103 -4298.99524 47649.56301 104 10903.14281 -4298.99524 105 13849.58930 10903.14281 106 -35899.70759 13849.58930 107 -21135.24029 -35899.70759 108 27075.53204 -21135.24029 109 2168.30495 27075.53204 110 9347.35468 2168.30495 111 6262.56626 9347.35468 112 -17921.48319 6262.56626 113 7451.54542 -17921.48319 114 -10052.69102 7451.54542 115 -6192.99932 -10052.69102 116 -8620.75945 -6192.99932 117 39220.86161 -8620.75945 118 -35334.19029 39220.86161 119 -34391.25173 -35334.19029 120 -11033.35493 -34391.25173 121 -38405.16575 -11033.35493 122 -28518.05125 -38405.16575 123 -71619.23916 -28518.05125 124 58721.09947 -71619.23916 125 12457.08959 58721.09947 126 5053.40988 12457.08959 127 -35560.79787 5053.40988 128 16744.00235 -35560.79787 129 -1311.85449 16744.00235 130 -7640.92877 -1311.85449 131 -27218.89093 -7640.92877 132 15347.24417 -27218.89093 133 -69714.11314 15347.24417 134 -11452.95326 -69714.11314 135 75689.71901 -11452.95326 136 -1658.10440 75689.71901 137 45664.29384 -1658.10440 138 -33566.25443 45664.29384 139 -3990.35253 -33566.25443 140 -19735.06652 -3990.35253 141 -28097.47854 -19735.06652 142 12964.00702 -28097.47854 143 -8126.90361 12964.00702 144 3514.30084 -8126.90361 145 -33995.14486 3514.30084 146 -2888.06724 -33995.14486 147 -20716.72490 -2888.06724 148 -5462.93781 -20716.72490 149 -4257.33307 -5462.93781 150 35784.21577 -4257.33307 151 24488.90835 35784.21577 152 -17546.08380 24488.90835 153 44968.81845 -17546.08380 154 24262.52367 44968.81845 155 20366.97621 24262.52367 156 27403.10064 20366.97621 157 22558.45456 27403.10064 158 -19132.46386 22558.45456 159 -3943.91461 -19132.46386 160 19181.04977 -3943.91461 161 -20268.57120 19181.04977 162 -16563.36847 -20268.57120 163 3586.21046 -16563.36847 164 -16717.60936 3586.21046 165 36368.12609 -16717.60936 166 5953.73709 36368.12609 167 24156.52882 5953.73709 168 17633.01975 24156.52882 169 72415.25487 17633.01975 170 13047.05478 72415.25487 171 2197.08783 13047.05478 172 -283.91364 2197.08783 173 18825.45159 -283.91364 174 -84449.68761 18825.45159 175 15108.76463 -84449.68761 176 -33521.45373 15108.76463 177 8517.00026 -33521.45373 178 -11975.80454 8517.00026 179 36540.21971 -11975.80454 180 -26634.41095 36540.21971 181 -60739.68987 -26634.41095 182 5736.37107 -60739.68987 183 26076.02674 5736.37107 184 42029.49517 26076.02674 185 3031.21024 42029.49517 186 12935.80988 3031.21024 187 14756.88250 12935.80988 188 -1441.99545 14756.88250 189 14081.31360 -1441.99545 190 27666.71469 14081.31360 191 11684.07754 27666.71469 192 -19327.53844 11684.07754 193 -20696.39693 -19327.53844 194 -472.17236 -20696.39693 195 1711.10049 -472.17236 196 -9659.48215 1711.10049 197 704.98200 -9659.48215 198 -23841.08339 704.98200 199 -8307.29664 -23841.08339 200 5306.01056 -8307.29664 201 9334.50984 5306.01056 202 14509.62400 9334.50984 203 347.77189 14509.62400 204 -21644.91521 347.77189 205 -3546.99712 -21644.91521 206 -15567.82447 -3546.99712 207 50055.30204 -15567.82447 208 26473.38992 50055.30204 209 7883.18184 26473.38992 210 9416.53186 7883.18184 211 12922.97848 9416.53186 212 -1066.52060 12922.97848 213 4805.98484 -1066.52060 214 -24358.30973 4805.98484 215 -10695.15857 -24358.30973 216 -29515.63779 -10695.15857 217 7044.81888 -29515.63779 218 10583.11348 7044.81888 219 27676.75548 10583.11348 220 -18468.72380 27676.75548 221 -32024.12292 -18468.72380 222 -55373.60902 -32024.12292 223 340.87940 -55373.60902 224 8375.18731 340.87940 225 5355.00092 8375.18731 226 17200.42427 5355.00092 227 -21838.24561 17200.42427 228 -26835.06367 -21838.24561 229 -59478.25015 -26835.06367 230 -10973.46033 -59478.25015 231 -13746.04441 -10973.46033 232 -36725.50802 -13746.04441 233 -4361.65297 -36725.50802 234 -17910.76376 -4361.65297 235 -12694.87786 -17910.76376 236 53513.99412 -12694.87786 237 1433.55821 53513.99412 238 -4270.69738 1433.55821 239 -13745.65009 -4270.69738 240 8421.71025 -13745.65009 241 -19633.49856 8421.71025 242 -14033.98363 -19633.49856 243 -3228.80810 -14033.98363 244 1333.17937 -3228.80810 245 -3696.73029 1333.17937 246 -34147.18670 -3696.73029 247 24461.88553 -34147.18670 248 9578.34396 24461.88553 249 -8861.25266 9578.34396 250 17155.42935 -8861.25266 251 5710.10336 17155.42935 252 12385.22920 5710.10336 253 3518.42688 12385.22920 254 -6792.00279 3518.42688 255 8876.56821 -6792.00279 256 -12091.52252 8876.56821 257 17109.94531 -12091.52252 258 15971.82214 17109.94531 259 6750.99521 15971.82214 260 -1997.38968 6750.99521 261 1454.61676 -1997.38968 262 9172.29191 1454.61676 263 -4530.12013 9172.29191 264 -12345.74724 -4530.12013 265 12095.46304 -12345.74724 266 -13771.19778 12095.46304 267 4738.22989 -13771.19778 268 52727.10689 4738.22989 269 -9751.74114 52727.10689 270 -3960.24266 -9751.74114 271 7610.07914 -3960.24266 272 5720.07428 7610.07914 273 7204.32321 5720.07428 274 9754.30241 7204.32321 275 -9641.70626 9754.30241 276 11510.80817 -9641.70626 277 10951.10329 11510.80817 278 -4098.58943 10951.10329 279 21416.05188 -4098.58943 280 1445.36551 21416.05188 281 8129.50701 1445.36551 282 7031.44229 8129.50701 283 16389.03335 7031.44229 284 10773.53438 16389.03335 285 -22401.63529 10773.53438 286 1895.28429 -22401.63529 287 -5398.14417 1895.28429 288 8137.64420 -5398.14417 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7l93k1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/880xl1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/9e0ak1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/wessaorg/rcomp/tmp/1048zv1324466480.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/11d9kl1324466480.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/12iqx11324466480.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/13kwmp1324466480.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/14zzml1324466480.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/159p6h1324466480.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/16q4241324466480.tab") + } > > try(system("convert tmp/1de2p1324466480.ps tmp/1de2p1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/264xe1324466480.ps tmp/264xe1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/3q5ap1324466480.ps tmp/3q5ap1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/4795c1324466480.ps tmp/4795c1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/5juzk1324466480.ps tmp/5juzk1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/643dw1324466480.ps tmp/643dw1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/7l93k1324466480.ps tmp/7l93k1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/880xl1324466480.ps tmp/880xl1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/9e0ak1324466480.ps tmp/9e0ak1324466480.png",intern=TRUE)) character(0) > try(system("convert tmp/1048zv1324466480.ps tmp/1048zv1324466480.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 10.162 0.704 10.878