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 + ,56 + ,396 + ,81 + ,3 + ,79 + ,30 + ,115 + ,94 + ,112285 + ,24188 + ,146283 + ,144 + ,145 + ,869 + ,120982 + ,56 + ,297 + ,55 + ,4 + ,58 + ,28 + ,109 + ,103 + ,84786 + ,18273 + ,98364 + ,103 + ,101 + ,1530 + ,176508 + ,54 + ,559 + ,50 + ,12 + ,60 + ,38 + ,146 + ,93 + ,83123 + ,14130 + ,86146 + ,98 + ,98 + ,2172 + ,179321 + ,89 + ,967 + ,125 + ,2 + ,108 + ,30 + ,116 + ,103 + ,101193 + ,32287 + ,96933 + ,135 + ,132 + ,901 + ,123185 + ,40 + ,270 + ,40 + ,1 + ,49 + ,22 + ,68 + ,51 + ,38361 + ,8654 + ,79234 + ,61 + ,60 + ,463 + ,52746 + ,25 + ,143 + ,37 + ,3 + ,0 + ,26 + ,101 + ,70 + ,68504 + ,9245 + ,42551 + ,39 + ,38 + ,3201 + ,385534 + ,92 + ,1562 + ,63 + ,0 + ,121 + ,25 + ,96 + ,91 + ,119182 + ,33251 + ,195663 + ,150 + ,144 + ,371 + ,33170 + ,18 + ,109 + ,44 + ,0 + ,1 + ,18 + ,67 + ,22 + ,22807 + ,1271 + ,6853 + ,5 + ,5 + ,1192 + ,101645 + ,63 + ,371 + ,88 + ,0 + ,20 + ,11 + ,44 + ,38 + ,17140 + ,5279 + ,21529 + ,28 + ,28 + ,1583 + ,149061 + ,44 + ,656 + ,66 + ,5 + ,43 + ,26 + ,100 + ,93 + ,116174 + ,27101 + ,95757 + ,84 + ,84 + ,1439 + ,165446 + ,33 + ,511 + ,57 + ,0 + ,69 + ,25 + ,93 + ,60 + ,57635 + ,16373 + ,85584 + ,80 + ,79 + ,1764 + ,237213 + ,84 + ,655 + ,74 + ,0 + ,78 + ,38 + ,140 + ,123 + ,66198 + ,19716 + ,143983 + ,130 + ,127 + ,1495 + ,173326 + ,88 + ,465 + ,49 + ,7 + ,86 + ,44 + ,166 + ,148 + ,71701 + ,17753 + ,75851 + ,82 + ,78 + ,1373 + ,133131 + ,55 + ,525 + ,52 + ,7 + ,44 + ,30 + ,99 + ,90 + ,57793 + ,9028 + ,59238 + ,60 + ,60 + ,2187 + ,258873 + ,60 + ,885 + ,88 + ,3 + ,104 + ,40 + ,139 + ,124 + ,80444 + ,18653 + ,93163 + ,131 + ,131 + ,1491 + ,180083 + ,66 + ,497 + ,36 + ,9 + ,63 + ,34 + ,130 + ,70 + ,53855 + ,8828 + ,96037 + ,84 + ,84 + ,4041 + ,324799 + ,154 + ,1436 + ,108 + ,0 + ,158 + ,47 + ,181 + ,168 + ,97668 + ,29498 + ,151511 + ,140 + ,133 + ,1706 + ,230964 + ,53 + ,612 + ,43 + ,4 + ,102 + ,30 + ,116 + ,115 + ,133824 + ,27563 + ,136368 + ,151 + ,150 + ,2152 + ,236785 + ,119 + ,865 + ,75 + ,3 + ,77 + ,31 + ,116 + ,71 + ,101481 + ,18293 + ,112642 + ,91 + ,91 + ,1036 + ,135473 + ,41 + ,385 + ,32 + ,0 + ,82 + ,23 + ,88 + ,66 + ,99645 + ,22530 + ,94728 + ,138 + ,132 + ,1882 + ,202925 + ,61 + ,567 + ,44 + ,7 + ,115 + ,36 + ,139 + ,134 + ,114789 + ,15977 + ,105499 + ,150 + ,136 + ,1929 + ,215147 + ,58 + ,639 + ,85 + ,0 + ,101 + ,36 + ,135 + ,117 + ,99052 + ,35082 + ,121527 + ,124 + ,124 + ,2242 + ,344297 + ,75 + ,963 + ,86 + ,1 + ,80 + ,30 + ,108 + ,108 + ,67654 + ,16116 + ,127766 + ,119 + ,118 + ,1220 + ,153935 + ,33 + ,398 + ,56 + ,5 + ,50 + ,25 + ,89 + ,84 + ,65553 + ,15849 + ,98958 + ,73 + ,70 + ,1289 + ,132943 + ,40 + ,410 + ,50 + ,7 + ,83 + ,39 + ,156 + ,156 + ,97500 + ,16026 + ,77900 + ,110 + ,107 + ,2515 + ,174724 + ,92 + ,966 + ,135 + ,0 + ,123 + ,34 + ,129 + ,120 + ,69112 + ,26569 + ,85646 + ,123 + ,119 + ,2147 + ,174415 + ,100 + ,801 + ,63 + ,0 + ,73 + ,31 + ,118 + ,114 + ,82753 + ,24785 + ,98579 + ,90 + ,89 + ,2352 + ,225548 + ,112 + ,892 + ,81 + ,5 + ,81 + ,31 + ,118 + ,94 + ,85323 + ,17569 + ,130767 + ,116 + ,112 + ,1638 + ,223632 + ,73 + ,513 + ,52 + ,0 + ,105 + ,33 + ,125 + ,120 + ,72654 + ,23825 + ,131741 + ,113 + ,108 + ,1222 + ,124817 + ,40 + ,469 + ,44 + ,0 + ,47 + ,25 + ,95 + ,81 + ,30727 + ,7869 + ,53907 + ,56 + ,52 + ,1812 + ,221698 + ,45 + ,683 + ,113 + ,0 + ,105 + ,33 + ,126 + ,110 + ,77873 + ,14975 + ,178812 + ,115 + ,112 + ,1677 + ,210767 + ,60 + ,643 + ,39 + ,3 + ,94 + ,35 + ,135 + ,133 + ,117478 + ,37791 + ,146761 + ,119 + ,116 + ,1579 + ,170266 + ,62 + ,535 + ,73 + ,4 + ,44 + ,42 + ,154 + ,122 + ,74007 + ,9605 + ,82036 + ,129 + ,123 + ,1731 + ,260561 + ,75 + ,625 + ,48 + ,1 + ,114 + ,43 + ,165 + ,158 + ,90183 + ,27295 + ,163253 + ,127 + ,125 + ,807 + ,84853 + ,31 + ,264 + ,33 + ,4 + ,38 + ,30 + ,113 + ,109 + ,61542 + ,2746 + ,27032 + ,27 + ,27 + ,2452 + ,294424 + ,77 + ,992 + ,59 + ,2 + ,107 + ,33 + ,127 + ,124 + ,101494 + ,34461 + ,171975 + ,175 + ,162 + ,829 + ,101011 + ,34 + ,238 + ,41 + ,0 + ,30 + ,13 + ,52 + ,39 + ,27570 + ,8098 + ,65990 + ,35 + ,32 + ,1940 + ,215641 + ,46 + ,818 + ,69 + ,0 + ,71 + ,32 + ,121 + ,92 + ,55813 + ,4787 + ,86572 + ,64 + ,64 + ,2662 + ,325107 + ,99 + ,937 + ,64 + ,0 + ,84 + ,36 + ,136 + ,126 + ,79215 + ,24919 + ,159676 + ,96 + ,92 + ,186 + ,7176 + ,17 + ,70 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1423 + ,603 + ,1929 + ,0 + ,0 + ,1499 + ,167542 + ,66 + ,507 + ,59 + ,2 + ,59 + ,28 + ,108 + ,70 + ,55461 + ,16329 + ,85371 + ,84 + ,83 + ,865 + ,106408 + ,30 + ,260 + ,32 + ,1 + ,33 + ,14 + ,46 + ,37 + ,31081 + ,12558 + ,58391 + ,41 + ,41 + ,1793 + ,96560 + ,76 + ,503 + ,129 + ,0 + ,42 + ,17 + ,54 + ,38 + ,22996 + ,7784 + ,31580 + ,47 + ,47 + ,2527 + ,265769 + ,146 + ,927 + ,37 + ,2 + ,96 + ,32 + ,124 + ,120 + ,83122 + ,28522 + ,136815 + ,126 + ,120 + ,2747 + ,269651 + ,67 + ,1269 + ,31 + ,10 + ,106 + ,30 + ,115 + ,93 + ,70106 + ,22265 + ,120642 + ,105 + ,105 + ,1324 + ,149112 + ,56 + ,537 + ,65 + ,6 + ,56 + ,35 + ,128 + ,95 + ,60578 + ,14459 + ,69107 + ,80 + ,79 + ,2702 + ,175824 + ,107 + ,910 + ,107 + ,0 + ,57 + ,20 + ,80 + ,77 + ,39992 + ,14526 + ,50495 + ,70 + ,65 + ,1383 + ,152871 + ,58 + ,532 + ,74 + ,5 + ,59 + ,28 + ,97 + ,90 + ,79892 + ,22240 + ,108016 + ,73 + ,70 + ,1179 + ,111665 + ,34 + ,345 + ,54 + ,4 + ,39 + ,28 + ,104 + ,80 + ,49810 + ,11802 + ,46341 + ,57 + ,55 + ,2099 + ,116408 + ,61 + ,918 + ,76 + ,1 + ,34 + ,39 + ,59 + ,31 + ,71570 + ,7623 + ,78348 + ,40 + ,39 + ,4308 + ,362301 + ,119 + ,1635 + ,715 + ,2 + ,76 + ,34 + ,125 + ,110 + ,100708 + ,11912 + ,79336 + ,68 + ,67 + ,918 + ,78800 + ,42 + ,330 + ,57 + ,2 + ,20 + ,26 + ,82 + ,66 + ,33032 + ,7935 + ,56968 + ,21 + ,21 + ,1831 + ,183167 + ,66 + ,557 + ,66 + ,0 + ,91 + ,39 + ,149 + ,138 + ,82875 + ,18220 + ,93176 + ,127 + ,127 + ,3373 + ,277965 + ,89 + ,1178 + ,106 + ,8 + ,115 + ,39 + ,149 + ,133 + ,139077 + ,19199 + ,161632 + ,154 + ,152 + ,1713 + ,150629 + ,44 + ,740 + ,54 + ,3 + ,85 + ,33 + ,122 + ,113 + ,71595 + ,19918 + ,87850 + ,116 + ,113 + ,1438 + ,168809 + ,66 + ,452 + ,32 + ,0 + ,76 + ,28 + ,118 + ,100 + ,72260 + ,21884 + ,127969 + ,102 + ,99 + ,496 + ,24188 + ,24 + ,218 + ,20 + ,0 + ,8 + ,4 + ,12 + ,7 + ,5950 + ,2694 + ,15049 + ,7 + ,7 + ,2253 + ,329267 + ,259 + ,764 + ,71 + ,8 + ,79 + ,39 + ,144 + ,140 + ,115762 + ,15808 + ,155135 + ,148 + ,141 + ,744 + ,65029 + ,17 + ,255 + ,21 + ,5 + ,21 + ,18 + ,67 + ,61 + ,32551 + ,3597 + ,25109 + ,21 + ,21 + ,1161 + ,101097 + ,64 + ,454 + ,70 + ,3 + ,30 + ,14 + ,52 + ,41 + ,31701 + ,5296 + ,45824 + ,35 + ,35 + ,2352 + ,218946 + ,41 + ,866 + ,112 + ,1 + ,76 + ,29 + ,108 + ,96 + ,80670 + ,25239 + ,102996 + ,112 + ,109 + ,2144 + ,244052 + ,68 + ,574 + ,66 + ,5 + ,101 + ,44 + ,166 + ,164 + ,143558 + ,29801 + ,160604 + ,137 + ,133 + ,4691 + ,341570 + ,168 + ,1276 + ,190 + ,1 + ,94 + ,21 + ,80 + ,78 + ,117105 + ,18450 + ,158051 + ,135 + ,123 + ,1112 + ,103597 + ,43 + ,379 + ,66 + ,1 + ,27 + ,16 + ,60 + ,49 + ,23789 + ,7132 + ,44547 + ,26 + ,26 + ,2694 + ,233328 + ,132 + ,825 + ,165 + ,5 + ,92 + ,28 + ,107 + ,102 + ,120733 + ,34861 + ,162647 + ,230 + ,230 + ,1973 + ,256462 + ,105 + ,798 + ,56 + ,0 + ,123 + ,35 + ,127 + ,124 + ,105195 + ,35940 + ,174141 + ,181 + ,166 + ,1769 + ,206161 + ,71 + ,663 + ,61 + ,12 + ,75 + ,28 + ,107 + ,99 + ,73107 + ,16688 + ,60622 + ,71 + ,68 + ,3148 + ,311473 + ,112 + ,1069 + ,53 + ,8 + ,128 + ,38 + ,146 + ,129 + ,132068 + ,24683 + ,179566 + ,147 + ,147 + ,2474 + ,235800 + ,94 + ,921 + ,127 + ,8 + ,105 + ,23 + ,84 + ,62 + ,149193 + ,46230 + ,184301 + ,190 + ,179 + ,2084 + ,177939 + ,82 + ,858 + ,63 + ,8 + ,55 + ,36 + ,141 + ,73 + ,46821 + ,10387 + ,75661 + ,64 + ,61 + ,1954 + ,207176 + ,70 + ,711 + ,38 + ,8 + ,56 + ,32 + ,123 + ,114 + ,87011 + ,21436 + ,96144 + ,105 + ,101 + ,1226 + ,196553 + ,57 + ,503 + ,50 + ,2 + ,41 + ,29 + ,111 + ,99 + ,95260 + ,30546 + ,129847 + ,107 + ,108 + ,1389 + ,174184 + ,53 + ,382 + ,52 + ,0 + ,72 + ,25 + ,98 + ,70 + ,55183 + ,19746 + ,117286 + ,94 + ,90 + ,1496 + ,143246 + ,103 + ,464 + ,42 + ,5 + ,67 + ,27 + ,105 + ,104 + ,106671 + ,15977 + ,71180 + ,116 + ,114 + ,2269 + ,187559 + ,121 + ,717 + ,76 + ,8 + ,75 + ,36 + ,135 + ,116 + ,73511 + ,22583 + ,109377 + ,106 + ,103 + ,1833 + ,187681 + ,62 + ,690 + ,67 + ,2 + ,114 + ,28 + ,107 + ,91 + ,92945 + ,17274 + ,85298 + ,143 + ,142 + ,1268 + ,119016 + ,52 + ,462 + ,50 + ,5 + ,118 + ,23 + ,85 + ,74 + ,78664 + ,16469 + ,73631 + ,81 + ,79 + ,1943 + ,182192 + ,52 + ,657 + ,53 + ,12 + ,77 + ,40 + ,155 + ,138 + ,70054 + ,14251 + ,86767 + ,89 + ,88 + ,893 + ,73566 + ,32 + ,385 + ,39 + ,6 + ,22 + ,23 + ,88 + ,67 + ,22618 + ,3007 + ,23824 + ,26 + ,25 + ,1762 + ,194979 + ,62 + ,577 + ,50 + ,7 + ,66 + ,40 + ,155 + ,151 + ,74011 + ,16851 + ,93487 + ,84 + ,83 + ,1403 + ,167488 + ,45 + ,619 + ,77 + ,2 + ,69 + ,28 + ,104 + ,72 + ,83737 + ,21113 + ,82981 + ,113 + ,113 + ,1425 + ,143756 + ,46 + ,479 + ,57 + ,0 + ,105 + ,34 + ,132 + ,120 + ,69094 + ,17401 + ,73815 + ,120 + ,118 + ,1857 + ,275541 + ,63 + ,817 + ,73 + ,4 + ,116 + ,33 + ,127 + ,115 + ,93133 + ,23958 + ,94552 + ,110 + ,110 + ,1840 + ,243199 + ,75 + ,752 + ,34 + ,3 + ,88 + ,28 + ,108 + ,105 + ,95536 + ,23567 + ,132190 + ,134 + ,129 + ,1502 + ,182999 + ,88 + ,430 + ,39 + ,6 + ,73 + ,34 + ,129 + ,104 + ,225920 + ,13065 + ,128754 + ,54 + ,51 + ,1441 + ,135649 + ,46 + ,451 + ,46 + ,2 + ,99 + ,30 + ,116 + ,108 + ,62133 + ,15358 + ,66363 + ,96 + ,93 + ,1420 + ,152299 + ,53 + ,537 + ,63 + ,0 + ,62 + ,33 + ,122 + ,98 + ,61370 + ,14587 + ,67808 + ,78 + ,76 + ,1416 + ,120221 + ,37 + ,519 + ,35 + ,1 + ,53 + ,22 + ,85 + ,69 + ,43836 + ,12770 + ,61724 + ,51 + ,49 + ,2970 + ,346485 + ,90 + ,1000 + ,106 + ,0 + ,118 + ,38 + ,147 + ,111 + ,106117 + ,24021 + ,131722 + ,121 + ,118 + ,1317 + ,145790 + ,63 + ,637 + ,43 + ,5 + ,30 + ,26 + ,99 + ,99 + ,38692 + ,9648 + ,68580 + ,38 + ,38 + ,1644 + ,193339 + ,78 + ,465 + ,47 + ,2 + ,100 + ,35 + ,87 + ,71 + ,84651 + ,20537 + ,106175 + ,145 + ,141 + ,870 + ,80953 + ,25 + ,437 + ,31 + ,0 + ,49 + ,8 + ,28 + ,27 + ,56622 + ,7905 + ,55792 + ,59 + ,58 + ,1654 + ,122774 + ,45 + ,711 + ,162 + ,0 + ,24 + ,24 + ,90 + ,69 + ,15986 + ,4527 + ,25157 + ,27 + ,27 + ,1054 + ,130585 + ,46 + ,299 + ,57 + ,5 + ,67 + ,29 + ,109 + ,107 + ,95364 + ,30495 + ,76669 + ,91 + ,91 + ,937 + ,112611 + ,41 + ,248 + ,36 + ,0 + ,46 + ,20 + ,78 + ,73 + ,26706 + ,7117 + ,57283 + ,48 + ,48 + ,3004 + ,286468 + ,144 + ,1162 + ,263 + ,1 + ,57 + ,29 + ,111 + ,107 + ,89691 + ,17719 + ,105805 + ,68 + ,63 + ,2008 + ,241066 + ,82 + ,714 + ,78 + ,0 + ,75 + ,45 + ,158 + ,93 + ,67267 + ,27056 + ,129484 + ,58 + ,56 + ,2547 + ,148446 + ,91 + ,905 + ,63 + ,1 + ,135 + ,37 + ,141 + ,129 + ,126846 + ,33473 + ,72413 + ,150 + ,144 + ,1885 + ,204713 + ,71 + ,649 + ,54 + ,1 + ,68 + ,33 + ,122 + ,69 + ,41140 + ,9758 + ,87831 + ,74 + ,73 + ,1626 + ,182079 + ,63 + ,512 + ,63 + ,2 + ,124 + ,33 + ,124 + ,118 + ,102860 + ,21115 + ,96971 + ,181 + ,168 + ,1468 + ,140344 + ,53 + ,472 + ,77 + ,6 + ,33 + ,25 + ,93 + ,73 + ,51715 + ,7236 + ,71299 + ,65 + ,64 + ,2445 + ,220516 + ,62 + ,905 + ,79 + ,1 + ,98 + ,32 + ,124 + ,119 + ,55801 + ,13790 + ,77494 + ,97 + ,97 + ,1964 + ,243060 + ,63 + ,786 + ,110 + ,4 + ,58 + ,29 + ,112 + ,104 + ,111813 + ,32902 + ,120336 + ,121 + ,117 + ,1381 + ,162765 + ,32 + ,489 + ,56 + ,2 + ,68 + ,28 + ,108 + ,107 + ,120293 + ,25131 + ,93913 + ,99 + ,100 + ,1369 + ,182613 + ,39 + ,479 + ,56 + ,3 + ,81 + ,28 + ,99 + ,99 + ,138599 + ,30910 + ,136048 + ,152 + ,149 + ,1659 + ,232138 + ,62 + ,617 + ,43 + ,0 + ,131 + ,31 + ,117 + ,90 + ,161647 + ,35947 + ,181248 + ,188 + ,187 + ,2888 + ,265318 + ,117 + ,925 + ,111 + ,10 + ,110 + ,52 + ,199 + ,197 + ,115929 + ,29848 + ,146123 + ,138 + ,127 + ,1290 + ,85574 + ,34 + ,351 + ,71 + ,0 + ,37 + ,21 + ,78 + ,36 + ,24266 + ,6943 + ,32036 + ,40 + ,37 + ,2845 + ,310839 + ,92 + ,1144 + ,62 + ,9 + ,130 + ,24 + ,91 + ,85 + ,162901 + ,42705 + ,186646 + ,254 + ,245 + ,1982 + ,225060 + ,93 + ,669 + ,56 + ,7 + ,93 + ,41 + ,158 + ,139 + ,109825 + ,31808 + ,102255 + ,87 + ,87 + ,1904 + ,232317 + ,54 + ,707 + ,74 + ,0 + ,118 + ,33 + ,126 + ,106 + ,129838 + ,26675 + ,168237 + ,178 + ,177 + ,1391 + ,144966 + ,144 + ,458 + ,60 + ,0 + ,39 + ,32 + ,122 + ,50 + ,37510 + ,8435 + ,64219 + ,51 + ,49 + ,602 + ,43287 + ,14 + ,214 + ,43 + ,4 + ,13 + ,19 + ,71 + ,64 + ,43750 + ,7409 + ,19630 + ,49 + ,49 + ,1743 + ,155754 + ,61 + ,599 + ,68 + ,4 + ,74 + ,20 + ,75 + ,31 + ,40652 + ,14993 + ,76825 + ,73 + ,73 + ,1559 + ,164709 + ,109 + ,572 + ,53 + ,0 + ,81 + ,31 + ,115 + ,63 + ,87771 + ,36867 + ,115338 + ,176 + ,177 + ,2014 + ,201940 + ,38 + ,897 + ,87 + ,0 + ,109 + ,31 + ,119 + ,92 + ,85872 + ,33835 + ,109427 + ,94 + ,94 + ,2143 + ,235454 + ,73 + ,819 + ,46 + ,0 + ,151 + ,32 + ,124 + ,106 + ,89275 + ,24164 + ,118168 + ,120 + ,117 + ,2146 + ,220801 + ,75 + ,720 + ,105 + ,1 + ,51 + ,18 + ,72 + ,63 + ,44418 + ,12607 + ,84845 + ,66 + ,60 + ,874 + ,99466 + ,50 + ,273 + ,32 + ,0 + ,28 + ,23 + ,91 + ,69 + ,192565 + ,22609 + ,153197 + ,56 + ,55 + ,1590 + ,92661 + ,61 + ,508 + ,133 + ,1 + ,40 + ,17 + ,45 + ,41 + ,35232 + ,5892 + ,29877 + ,39 + ,39 + ,1590 + ,133328 + ,55 + ,506 + ,79 + ,0 + ,56 + ,20 + ,78 + ,56 + ,40909 + ,17014 + ,63506 + ,66 + ,64 + ,1210 + ,61361 + ,77 + ,451 + ,51 + ,0 + ,27 + ,12 + ,39 + ,25 + ,13294 + ,5394 + ,22445 + ,27 + ,26 + ,2072 + ,125930 + ,75 + ,699 + ,207 + ,4 + ,37 + ,17 + ,68 + ,65 + ,32387 + ,9178 + ,47695 + ,65 + ,64 + ,1281 + ,100750 + ,72 + ,407 + ,67 + ,0 + ,83 + ,30 + ,119 + ,93 + ,140867 + ,6440 + ,68370 + ,58 + ,58 + ,1401 + ,224549 + ,50 + ,465 + ,47 + ,4 + ,54 + ,31 + ,117 + ,114 + ,120662 + ,21916 + ,146304 + ,98 + ,95 + ,834 + ,82316 + ,32 + ,245 + ,34 + ,4 + ,27 + ,10 + ,39 + ,38 + ,21233 + ,4011 + ,38233 + ,25 + ,25 + ,1105 + ,102010 + ,53 + ,370 + ,66 + ,3 + ,28 + ,13 + ,50 + ,44 + ,44332 + ,5818 + ,42071 + ,26 + ,26 + ,1272 + ,101523 + ,42 + ,316 + ,76 + ,0 + ,59 + ,22 + ,88 + ,87 + ,61056 + ,18647 + ,50517 + ,77 + ,76 + ,1944 + ,243511 + ,71 + ,603 + ,65 + ,0 + ,133 + ,42 + ,155 + ,110 + ,101338 + ,20556 + ,103950 + ,130 + ,129 + ,391 + ,22938 + ,10 + ,154 + ,9 + ,0 + ,12 + ,1 + ,0 + ,0 + ,1168 + ,238 + ,5841 + ,11 + ,11 + ,761 + ,41566 + ,35 + ,229 + ,42 + ,5 + ,0 + ,9 + ,36 + ,27 + ,13497 + ,70 + ,2341 + ,2 + ,2 + ,1605 + ,152474 + ,65 + ,577 + ,45 + ,0 + ,106 + ,32 + ,123 + ,83 + ,65567 + ,22392 + ,84396 + ,101 + ,101 + ,530 + ,61857 + ,25 + ,192 + ,25 + ,4 + ,23 + ,11 + ,32 + ,30 + ,25162 + ,3913 + ,24610 + ,31 + ,28 + ,1988 + ,99923 + ,66 + ,617 + ,115 + ,0 + ,44 + ,25 + ,99 + ,80 + ,32334 + ,12237 + ,35753 + ,36 + ,36 + ,1386 + ,132487 + ,41 + ,411 + ,97 + ,0 + ,71 + ,36 + ,136 + ,98 + ,40735 + ,8388 + ,55515 + ,120 + ,89 + ,2395 + ,317394 + ,86 + ,975 + ,53 + ,1 + ,116 + ,31 + ,117 + ,82 + ,91413 + ,22120 + ,209056 + ,195 + ,193 + ,387 + ,21054 + ,16 + ,146 + ,2 + ,0 + ,4 + ,0 + ,0 + ,0 + ,855 + ,338 + ,6622 + ,4 + ,4 + ,1742 + ,209641 + ,42 + ,705 + ,52 + ,5 + ,62 + ,24 + ,88 + ,60 + ,97068 + ,11727 + ,115814 + ,89 + ,84 + ,620 + ,22648 + ,19 + ,184 + ,44 + ,0 + ,12 + ,13 + ,39 + ,28 + ,44339 + ,3704 + ,11609 + ,24 + ,23 + ,449 + ,31414 + ,19 + ,200 + ,22 + ,0 + ,18 + ,8 + ,25 + ,9 + ,14116 + ,3988 + ,13155 + ,39 + ,39 + ,800 + ,46698 + ,45 + ,274 + ,35 + ,0 + ,14 + ,13 + ,52 + ,33 + ,10288 + ,3030 + ,18274 + ,14 + ,14 + ,1684 + ,131698 + ,65 + ,502 + ,74 + ,0 + ,60 + ,19 + ,75 + ,59 + ,65622 + ,13520 + ,72875 + ,78 + ,78 + ,1050 + ,91735 + ,35 + ,382 + ,103 + ,0 + ,7 + ,18 + ,71 + ,49 + ,16563 + ,1421 + ,10112 + ,15 + ,14 + ,2699 + ,244749 + ,95 + ,964 + ,144 + ,2 + ,98 + ,33 + ,124 + ,115 + ,76643 + ,20923 + ,142775 + ,106 + ,101 + ,1606 + ,184510 + ,49 + ,537 + ,60 + ,7 + ,64 + ,40 + ,151 + ,140 + ,110681 + ,20237 + ,68847 + ,83 + ,82 + ,1502 + ,79863 + ,37 + ,438 + ,134 + ,1 + ,29 + ,22 + ,71 + ,49 + ,29011 + ,3219 + ,17659 + ,24 + ,24 + ,1204 + ,128423 + ,64 + ,369 + ,89 + ,8 + ,32 + ,38 + ,145 + ,120 + ,92696 + ,3769 + ,20112 + ,37 + ,36 + ,1138 + ,97839 + ,38 + ,417 + ,42 + ,2 + ,25 + ,24 + ,87 + ,66 + ,94785 + ,12252 + ,61023 + ,77 + ,75 + ,568 + ,38214 + ,34 + ,276 + ,52 + ,0 + ,16 + ,8 + ,27 + ,21 + ,8773 + ,1888 + ,13983 + ,16 + ,16 + ,1459 + ,151101 + ,32 + ,514 + ,98 + ,2 + ,48 + ,35 + ,131 + ,124 + ,83209 + ,14497 + ,65176 + ,56 + ,55 + ,2158 + ,272458 + ,65 + ,822 + ,99 + ,0 + ,100 + ,43 + ,162 + ,152 + ,93815 + ,28864 + ,132432 + ,132 + ,131 + ,1111 + ,172494 + ,52 + ,389 + ,52 + ,0 + ,46 + ,43 + ,165 + ,139 + ,86687 + ,21721 + ,112494 + ,144 + ,131 + ,1421 + ,108043 + ,62 + ,466 + ,29 + ,1 + ,45 + ,14 + ,54 + ,38 + ,34553 + ,4821 + ,45109 + ,40 + ,39 + ,2833 + ,328107 + ,65 + ,1255 + ,125 + ,3 + ,129 + ,41 + ,159 + ,144 + ,105547 + ,33644 + ,170875 + ,153 + ,144 + ,1955 + ,250579 + ,83 + ,694 + ,106 + ,0 + ,130 + ,38 + ,147 + ,120 + ,103487 + ,15923 + ,180759 + ,143 + ,139 + ,2922 + ,351067 + ,95 + ,1024 + ,95 + ,3 + ,136 + ,45 + ,170 + ,160 + ,213688 + ,42935 + ,214921 + ,220 + ,211 + ,1002 + ,158015 + ,29 + ,400 + ,40 + ,0 + ,59 + ,31 + ,119 + ,114 + ,71220 + ,18864 + ,100226 + ,79 + ,78 + ,1060 + ,98866 + ,18 + ,397 + ,140 + ,0 + ,25 + ,13 + ,49 + ,39 + ,23517 + ,4977 + ,32043 + ,50 + ,50 + ,956 + ,85439 + ,33 + ,350 + ,43 + ,0 + ,32 + ,28 + ,104 + ,78 + ,56926 + ,7785 + ,54454 + ,39 + ,39 + ,2186 + ,229242 + ,247 + ,719 + ,128 + ,4 + ,63 + ,31 + ,120 + ,119 + ,91721 + ,17939 + ,78876 + ,95 + ,90 + ,3604 + ,351619 + ,139 + ,1277 + ,142 + ,4 + ,95 + ,40 + ,150 + ,141 + ,115168 + ,23436 + ,170745 + ,169 + ,166 + ,1035 + ,84207 + ,29 + ,356 + ,73 + ,11 + ,14 + ,30 + ,112 + ,101 + ,111194 + ,325 + ,6940 + ,12 + ,12 + ,1417 + ,120445 + ,118 + ,457 + ,72 + ,0 + ,36 + ,16 + ,59 + ,56 + ,51009 + ,13539 + ,49025 + ,63 + ,57 + ,3261 + ,324598 + ,110 + ,1402 + ,128 + ,0 + ,113 + ,37 + ,136 + ,133 + ,135777 + ,34538 + ,122037 + ,134 + ,133 + ,1587 + ,131069 + ,67 + ,600 + ,61 + ,4 + ,47 + ,30 + ,107 + ,83 + ,51513 + ,12198 + ,53782 + ,69 + ,69 + ,1424 + ,204271 + ,42 + ,480 + ,73 + ,0 + ,92 + ,35 + ,130 + ,116 + ,74163 + ,26924 + ,127748 + ,119 + ,119 + ,1701 + ,165543 + ,65 + ,595 + ,148 + ,1 + ,70 + ,32 + ,115 + ,90 + ,51633 + ,12716 + ,86839 + ,119 + ,119 + ,1249 + ,141722 + ,94 + ,436 + ,64 + ,0 + ,19 + ,27 + ,107 + ,36 + ,75345 + ,8172 + ,44830 + ,75 + ,65 + ,946 + ,116048 + ,64 + ,230 + ,45 + ,0 + ,50 + ,20 + ,75 + ,50 + ,33416 + ,10855 + ,77395 + ,63 + ,61 + ,1926 + ,250047 + ,81 + ,651 + ,58 + ,0 + ,41 + ,18 + ,71 + ,61 + ,83305 + ,11932 + ,89324 + ,55 + ,49 + ,3352 + ,299775 + ,95 + ,1367 + ,97 + ,9 + ,91 + ,31 + ,120 + ,97 + ,98952 + ,14300 + ,103300 + ,103 + ,101 + ,1641 + ,195838 + ,67 + ,564 + ,50 + ,1 + ,111 + ,31 + ,116 + ,98 + ,102372 + ,25515 + ,112283 + ,197 + ,196 + ,2035 + ,173260 + ,63 + ,716 + ,37 + ,3 + ,41 + ,21 + ,79 + ,78 + ,37238 + ,2805 + ,10901 + ,16 + ,15 + ,2312 + ,254488 + ,83 + ,747 + ,50 + ,10 + ,120 + ,39 + ,150 + ,117 + ,103772 + ,29402 + ,120691 + ,140 + ,136 + ,1369 + ,104389 + ,45 + ,467 + ,105 + ,5 + ,135 + ,41 + ,156 + ,148 + ,123969 + ,16440 + ,58106 + ,89 + ,89 + ,1577 + ,136084 + ,30 + ,671 + ,69 + ,0 + ,27 + ,13 + ,51 + ,41 + ,27142 + ,11221 + ,57140 + ,40 + ,40 + ,2201 + ,199476 + ,70 + ,861 + ,46 + ,2 + ,87 + ,32 + ,118 + ,105 + ,135400 + ,28732 + ,122422 + ,125 + ,123 + ,961 + ,92499 + ,32 + ,319 + ,57 + ,0 + ,25 + ,18 + ,71 + ,55 + ,21399 + ,5250 + ,25899 + ,21 + ,21 + ,1900 + ,224330 + ,83 + ,612 + ,52 + ,1 + ,131 + ,39 + ,144 + ,132 + ,130115 + ,28608 + ,139296 + ,167 + ,163 + ,1254 + ,135781 + ,31 + ,433 + ,98 + ,2 + ,45 + ,14 + ,47 + ,44 + ,24874 + ,8092 + ,52678 + ,32 + ,29 + ,1335 + ,74408 + ,67 + ,434 + ,61 + ,4 + ,29 + ,7 + ,28 + ,21 + ,34988 + ,4473 + ,23853 + ,36 + ,35 + ,1597 + ,81240 + ,66 + ,503 + ,89 + ,0 + ,58 + ,17 + ,68 + ,50 + ,45549 + ,1572 + ,17306 + ,13 + ,13 + ,207 + ,14688 + ,10 + ,85 + ,0 + ,0 + ,4 + ,0 + ,0 + ,0 + ,6023 + ,2065 + ,7953 + ,5 + ,5 + ,1645 + ,181633 + ,70 + ,564 + ,48 + ,2 + ,47 + ,30 + ,110 + ,73 + ,64466 + ,14817 + ,89455 + ,96 + ,96 + ,2429 + ,271856 + ,103 + ,824 + ,91 + ,1 + ,109 + ,37 + ,147 + ,86 + ,54990 + ,16714 + ,147866 + ,151 + ,151 + ,151 + ,7199 + ,5 + ,74 + ,0 + ,0 + ,7 + ,0 + ,0 + ,0 + ,1644 + ,556 + ,4245 + ,6 + ,6 + ,474 + ,46660 + ,20 + ,259 + ,7 + ,0 + ,12 + ,5 + ,15 + ,13 + ,6179 + ,2089 + ,21509 + ,13 + ,13 + ,141 + ,17547 + ,5 + ,69 + ,3 + ,0 + ,0 + ,1 + ,4 + ,4 + ,3926 + ,2658 + ,7670 + ,3 + ,3 + ,1639 + ,133368 + ,36 + ,535 + ,54 + ,1 + ,37 + ,16 + ,64 + ,57 + ,32755 + ,10695 + ,66675 + ,57 + ,56 + ,872 + ,95227 + ,34 + ,239 + ,70 + ,0 + ,37 + ,32 + ,111 + ,48 + ,34777 + ,1669 + ,14336 + ,23 + ,23 + ,1318 + ,152601 + ,48 + ,438 + ,36 + ,2 + ,46 + ,24 + ,85 + ,46 + ,73224 + ,16267 + ,53608 + ,61 + ,57 + ,1018 + ,98146 + ,40 + ,459 + ,37 + ,0 + ,15 + ,17 + ,68 + ,48 + ,27114 + ,7768 + ,30059 + ,21 + ,14 + ,1383 + ,79619 + ,43 + ,426 + ,123 + ,3 + ,42 + ,11 + ,40 + ,32 + ,20760 + ,7252 + ,29668 + ,43 + ,43 + ,1314 + ,59194 + ,31 + ,288 + ,247 + ,6 + ,7 + ,24 + ,80 + ,68 + ,37636 + ,6387 + ,22097 + ,20 + ,20 + ,1335 + ,139942 + ,42 + ,498 + ,46 + ,0 + ,54 + ,22 + ,88 + ,87 + ,65461 + ,18715 + ,96841 + ,82 + ,72 + ,1403 + ,118612 + ,46 + ,454 + ,72 + ,2 + ,54 + ,12 + ,48 + ,43 + ,30080 + ,7936 + ,41907 + ,90 + ,87 + ,910 + ,72880 + ,33 + ,376 + ,41 + ,0 + ,14 + ,19 + ,76 + ,67 + ,24094 + ,8643 + ,27080 + ,25 + ,21 + ,616 + ,65475 + ,18 + ,225 + ,24 + ,2 + ,16 + ,13 + ,51 + ,46 + ,69008 + ,7294 + ,35885 + ,60 + ,56 + ,1407 + ,99643 + ,55 + ,555 + ,45 + ,1 + ,33 + ,17 + ,67 + ,46 + ,54968 + ,4570 + ,41247 + ,61 + ,59 + ,771 + ,71965 + ,35 + ,252 + ,33 + ,1 + ,32 + ,15 + ,59 + ,56 + ,46090 + ,7185 + ,28313 + ,85 + ,82 + ,766 + ,77272 + ,59 + ,208 + ,27 + ,2 + ,21 + ,16 + ,61 + ,48 + ,27507 + ,10058 + ,36845 + ,43 + ,43 + ,473 + ,49289 + ,19 + ,130 + ,36 + ,1 + ,15 + ,24 + ,76 + ,44 + ,10672 + ,2342 + ,16548 + ,25 + ,25 + ,1376 + ,135131 + ,66 + ,481 + ,87 + ,0 + ,38 + ,15 + ,60 + ,60 + ,34029 + ,8509 + ,36134 + ,41 + ,38 + ,1232 + ,108446 + ,60 + ,389 + ,90 + ,1 + ,22 + ,17 + ,68 + ,65 + ,46300 + ,13275 + ,55764 + ,26 + ,25 + ,1521 + ,89746 + ,36 + ,565 + ,114 + ,3 + ,28 + ,18 + ,71 + ,55 + ,24760 + ,6816 + ,28910 + ,38 + ,38 + ,572 + ,44296 + ,25 + ,173 + ,31 + ,0 + ,10 + ,20 + ,76 + ,38 + ,18779 + ,1930 + ,13339 + ,12 + ,12 + ,1059 + ,77648 + ,47 + ,278 + ,45 + ,0 + ,31 + ,16 + ,62 + ,52 + ,21280 + ,8086 + ,25319 + ,29 + ,29 + ,1544 + ,181528 + ,54 + ,609 + ,69 + ,0 + ,32 + ,16 + ,61 + ,60 + ,40662 + ,10737 + ,66956 + ,49 + ,47 + ,1230 + ,134019 + ,53 + ,422 + ,51 + ,0 + ,32 + ,18 + ,67 + ,54 + ,28987 + ,8033 + ,47487 + ,46 + ,45 + ,1206 + ,124064 + ,40 + ,445 + ,34 + ,1 + ,43 + ,22 + ,88 + ,86 + ,22827 + ,7058 + ,52785 + ,41 + ,40 + ,1205 + ,92630 + ,40 + ,387 + ,60 + ,4 + ,27 + ,8 + ,30 + ,24 + ,18513 + ,6782 + ,44683 + ,31 + ,30 + ,1255 + ,121848 + ,39 + ,339 + ,45 + ,0 + ,37 + ,17 + ,64 + ,52 + ,30594 + ,5401 + ,35619 + ,41 + ,41 + ,613 + ,52915 + ,14 + ,181 + ,54 + ,0 + ,20 + ,18 + ,68 + ,49 + ,24006 + ,6521 + ,21920 + ,26 + ,25 + ,721 + ,81872 + ,45 + ,245 + ,25 + ,0 + ,32 + ,16 + ,64 + ,61 + ,27913 + ,10856 + ,45608 + ,23 + ,23 + ,1109 + ,58981 + ,36 + ,384 + ,38 + ,7 + ,0 + ,23 + ,91 + ,61 + ,42744 + ,2154 + ,7721 + ,14 + ,14 + ,740 + ,53515 + ,28 + ,212 + ,52 + ,2 + ,5 + ,22 + ,88 + ,81 + ,12934 + ,6117 + ,20634 + ,16 + ,16 + ,1126 + ,60812 + ,44 + ,399 + ,67 + ,0 + ,26 + ,13 + ,52 + ,43 + ,22574 + ,5238 + ,29788 + ,25 + ,26 + ,728 + ,56375 + ,30 + ,229 + ,74 + ,7 + ,10 + ,13 + ,49 + ,40 + ,41385 + ,4820 + ,31931 + ,21 + ,21 + ,689 + ,65490 + ,22 + ,224 + ,38 + ,3 + ,27 + ,16 + ,62 + ,40 + ,18653 + ,5615 + ,37754 + ,32 + ,27 + ,592 + ,80949 + ,17 + ,203 + ,30 + ,0 + ,11 + ,16 + ,61 + ,56 + ,18472 + ,4272 + ,32505 + ,9 + ,9 + ,995 + ,76302 + ,31 + ,333 + ,26 + ,0 + ,29 + ,20 + ,76 + ,68 + ,30976 + ,8702 + ,40557 + ,35 + ,33 + ,1613 + ,104011 + ,55 + ,384 + ,67 + ,6 + ,25 + ,22 + ,88 + ,79 + ,63339 + ,15340 + ,94238 + ,42 + ,42 + ,2048 + ,98104 + ,54 + ,636 + ,132 + ,2 + ,55 + ,17 + ,66 + ,47 + ,25568 + ,8030 + ,44197 + ,68 + ,68 + ,705 + ,67989 + ,21 + ,185 + ,42 + ,0 + ,23 + ,18 + ,71 + ,57 + ,33747 + ,9526 + ,43228 + ,32 + ,32 + ,301 + ,30989 + ,14 + ,93 + ,35 + ,0 + ,5 + ,17 + ,68 + ,41 + ,4154 + ,1278 + ,4103 + ,6 + ,6 + ,1803 + ,135458 + ,81 + ,581 + ,118 + ,3 + ,43 + ,12 + ,48 + ,29 + ,19474 + ,4236 + ,44144 + ,68 + ,67 + ,799 + ,73504 + ,35 + ,248 + ,68 + ,0 + ,23 + ,7 + ,25 + ,3 + ,35130 + ,3023 + ,32868 + ,33 + ,33 + ,861 + ,63123 + ,43 + ,304 + ,43 + ,1 + ,34 + ,17 + ,68 + ,60 + ,39067 + ,7196 + ,27640 + ,84 + ,77 + ,1186 + ,61254 + ,46 + ,344 + ,76 + ,1 + ,36 + ,14 + ,41 + ,30 + ,13310 + ,3394 + ,14063 + ,46 + ,46 + ,1451 + ,74914 + ,30 + ,407 + ,64 + ,0 + ,35 + ,23 + ,90 + ,79 + ,65892 + ,6371 + ,28990 + ,30 + ,30 + ,628 + ,31774 + ,23 + ,170 + ,48 + ,1 + ,0 + ,17 + ,66 + ,47 + ,4143 + ,1574 + ,4694 + ,0 + ,0 + ,1161 + ,81437 + ,38 + ,312 + ,64 + ,0 + ,37 + ,14 + ,54 + ,40 + ,28579 + ,9620 + ,42648 + ,36 + ,36 + ,1463 + ,87186 + ,54 + ,507 + ,56 + ,0 + ,28 + ,15 + ,59 + ,48 + ,51776 + ,6978 + ,64329 + ,47 + ,46 + ,742 + ,50090 + ,20 + ,224 + ,71 + ,0 + ,16 + ,17 + ,60 + ,36 + ,21152 + ,4911 + ,21928 + ,20 + ,18 + ,979 + ,65745 + ,53 + ,340 + ,75 + ,0 + ,26 + ,21 + ,77 + ,42 + ,38084 + ,8645 + ,25836 + ,50 + ,48 + ,675 + ,56653 + ,45 + ,168 + ,39 + ,0 + ,38 + ,18 + ,68 + ,49 + ,27717 + ,8987 + ,22779 + ,30 + ,29 + ,1241 + ,158399 + ,39 + ,443 + ,42 + ,0 + ,23 + ,18 + ,72 + ,57 + ,32928 + ,5544 + ,40820 + ,30 + ,28 + ,676 + ,46455 + ,20 + ,204 + ,39 + ,0 + ,22 + ,17 + ,67 + ,12 + ,11342 + ,3083 + ,27530 + ,34 + ,34 + ,1049 + ,73624 + ,24 + ,367 + ,93 + ,0 + ,30 + ,17 + ,64 + ,40 + ,19499 + ,6909 + ,32378 + ,33 + ,33 + ,620 + ,38395 + ,31 + ,210 + ,38 + ,0 + ,16 + ,16 + ,63 + ,43 + ,16380 + ,3189 + ,10824 + ,34 + ,34 + ,1081 + ,91899 + ,35 + ,335 + ,60 + ,0 + ,18 + ,15 + ,59 + ,33 + ,36874 + ,6745 + ,39613 + ,37 + ,33 + ,1688 + ,139526 + ,151 + ,364 + ,71 + ,0 + ,28 + ,21 + ,84 + ,77 + ,48259 + ,16724 + ,60865 + ,83 + ,80 + ,736 + ,52164 + ,52 + ,178 + ,52 + ,0 + ,32 + ,16 + ,64 + ,43 + ,16734 + ,4850 + ,19787 + ,32 + ,32 + ,617 + ,51567 + ,30 + ,206 + ,27 + ,2 + ,21 + ,14 + ,56 + ,45 + ,28207 + ,7025 + ,20107 + ,30 + ,30 + ,812 + ,70551 + ,31 + ,279 + ,59 + ,0 + ,23 + ,15 + ,54 + ,47 + ,30143 + ,6047 + ,36605 + ,43 + ,41 + ,1051 + ,84856 + ,29 + ,387 + ,40 + ,1 + ,29 + ,17 + ,67 + ,43 + ,41369 + ,7377 + ,40961 + ,41 + ,41 + ,1656 + ,102538 + ,57 + ,490 + ,79 + ,1 + ,50 + ,15 + ,58 + ,45 + ,45833 + ,9078 + ,48231 + ,51 + ,51 + ,705 + ,86678 + ,40 + ,238 + ,44 + ,0 + ,12 + ,15 + ,59 + ,50 + ,29156 + ,4605 + ,39725 + ,19 + ,18 + ,945 + ,85709 + ,44 + ,343 + ,65 + ,0 + ,21 + ,10 + ,40 + ,35 + ,35944 + ,3238 + ,21455 + ,37 + ,34 + ,554 + ,34662 + ,25 + ,232 + ,10 + ,0 + ,18 + ,6 + ,22 + ,7 + ,36278 + ,8100 + ,23430 + ,33 + ,31 + ,1597 + ,150580 + ,77 + ,530 + ,124 + ,0 + ,27 + ,22 + ,83 + ,71 + ,45588 + ,9653 + ,62991 + ,41 + ,39 + ,982 + ,99611 + ,35 + ,291 + ,81 + ,0 + ,41 + ,21 + ,81 + ,67 + ,45097 + ,8914 + ,49363 + ,54 + ,54 + ,222 + ,19349 + ,11 + ,67 + ,15 + ,0 + ,13 + ,1 + ,2 + ,0 + ,3895 + ,786 + ,9604 + ,14 + ,14 + ,1212 + ,99373 + ,63 + ,397 + ,92 + ,1 + ,12 + ,18 + ,72 + ,62 + ,28394 + ,6700 + ,24552 + ,25 + ,24 + ,1143 + ,86230 + ,44 + ,467 + ,42 + ,0 + ,21 + ,17 + ,61 + ,54 + ,18632 + ,5788 + ,31493 + ,25 + ,24 + ,435 + ,30837 + ,19 + ,178 + ,10 + ,0 + ,8 + ,4 + ,15 + ,4 + ,2325 + ,593 + ,3439 + ,8 + ,8 + ,532 + ,31706 + ,13 + ,175 + ,24 + ,0 + ,26 + ,10 + ,32 + ,25 + ,25139 + ,4506 + ,19555 + ,26 + ,26 + ,882 + ,89806 + ,42 + ,299 + ,64 + ,0 + ,27 + ,16 + ,62 + ,40 + ,27975 + ,6382 + ,21228 + ,20 + ,19 + ,608 + ,62088 + ,38 + ,154 + ,45 + ,1 + ,13 + ,16 + ,58 + ,38 + ,14483 + ,5621 + ,23177 + ,11 + ,11 + ,459 + ,40151 + ,29 + ,106 + ,22 + ,0 + ,16 + ,9 + ,36 + ,19 + ,13127 + ,3997 + ,22094 + ,14 + ,14 + ,578 + ,27634 + ,20 + ,189 + ,56 + ,0 + ,2 + ,16 + ,59 + ,17 + ,5839 + ,520 + ,2342 + ,3 + ,1 + ,826 + ,76990 + ,27 + ,194 + ,94 + ,0 + ,42 + ,17 + ,68 + ,67 + ,24069 + ,8891 + ,38798 + ,40 + ,39 + ,509 + ,37460 + ,20 + ,135 + ,19 + ,0 + ,5 + ,7 + ,21 + ,14 + ,3738 + ,999 + ,3255 + ,5 + ,5 + ,717 + ,54157 + ,19 + ,201 + ,35 + ,0 + ,37 + ,15 + ,55 + ,30 + ,18625 + ,7067 + ,24261 + ,38 + ,37 + ,637 + ,49862 + ,37 + ,207 + ,32 + ,0 + ,17 + ,14 + ,54 + ,54 + ,36341 + ,4639 + ,18511 + ,32 + ,32 + ,857 + ,84337 + ,26 + ,280 + ,35 + ,0 + ,38 + ,14 + ,55 + ,35 + ,24548 + ,5654 + ,40798 + ,41 + ,38 + ,830 + ,64175 + ,42 + ,260 + ,48 + ,0 + ,37 + ,18 + ,72 + ,59 + ,21792 + ,6928 + ,28893 + ,46 + ,47 + ,652 + ,59382 + ,49 + ,227 + ,49 + ,0 + ,29 + ,12 + ,41 + ,24 + ,26263 + ,1514 + ,21425 + ,47 + ,47 + ,707 + ,119308 + ,30 + ,239 + ,48 + ,0 + ,32 + ,16 + ,61 + ,58 + ,23686 + ,9238 + ,50276 + ,37 + ,37 + ,954 + ,76702 + ,49 + ,333 + ,62 + ,0 + ,35 + ,21 + ,67 + ,42 + ,49303 + ,8204 + ,37643 + ,51 + ,51 + ,1461 + ,103425 + ,67 + ,428 + ,96 + ,1 + ,17 + ,19 + ,76 + ,46 + ,25659 + ,5926 + ,30377 + ,49 + ,45 + ,672 + ,70344 + ,28 + ,230 + ,45 + ,0 + ,20 + ,16 + ,64 + ,61 + ,28904 + ,5785 + ,27126 + ,21 + ,21 + ,778 + ,43410 + ,19 + ,292 + ,63 + ,0 + ,7 + ,1 + ,3 + ,3 + ,2781 + ,4 + ,13 + ,1 + ,1 + ,1141 + ,104838 + ,49 + ,350 + ,71 + ,1 + ,46 + ,16 + ,63 + ,52 + ,29236 + ,5930 + ,42097 + ,44 + ,42 + ,680 + ,62215 + ,27 + ,186 + ,26 + ,0 + ,24 + ,10 + ,40 + ,25 + ,19546 + ,3710 + ,24451 + ,26 + ,26 + ,1090 + ,69304 + ,30 + ,326 + ,48 + ,6 + ,40 + ,19 + ,69 + ,40 + ,22818 + ,705 + ,14335 + ,21 + ,21 + ,616 + ,53117 + ,22 + ,155 + ,29 + ,3 + ,3 + ,12 + ,48 + ,32 + ,32689 + ,443 + ,5084 + ,4 + ,4 + ,285 + ,19764 + ,12 + ,75 + ,19 + ,1 + ,10 + ,2 + ,8 + ,4 + ,5752 + ,2416 + ,9927 + ,10 + ,10 + ,1145 + ,86680 + ,31 + ,361 + ,45 + ,2 + ,37 + ,14 + ,52 + ,49 + ,22197 + ,7747 + ,43527 + ,43 + ,43 + ,733 + ,84105 + ,20 + ,261 + ,45 + ,0 + ,17 + ,17 + ,66 + ,63 + ,20055 + ,5432 + ,27184 + ,34 + ,34 + ,888 + ,77945 + ,20 + ,299 + ,67 + ,0 + ,28 + ,19 + ,76 + ,67 + ,25272 + ,4913 + ,21610 + ,32 + ,31 + ,849 + ,89113 + ,39 + ,300 + ,30 + ,0 + ,19 + ,14 + ,43 + ,32 + ,82206 + ,2650 + ,20484 + ,20 + ,19 + ,1182 + ,91005 + ,29 + ,450 + ,36 + ,3 + ,29 + ,11 + ,39 + ,23 + ,32073 + ,2370 + ,20156 + ,34 + ,34 + ,528 + ,40248 + ,16 + ,183 + ,34 + ,1 + ,8 + ,4 + ,14 + ,7 + ,5444 + ,775 + ,6012 + ,6 + ,6 + ,642 + ,64187 + ,27 + ,238 + ,36 + ,0 + ,10 + ,16 + ,61 + ,54 + ,20154 + ,5576 + ,18475 + ,12 + ,11 + ,947 + ,50857 + ,21 + ,165 + ,34 + ,0 + ,15 + ,20 + ,71 + ,37 + ,36944 + ,1352 + ,12645 + ,24 + ,24 + ,819 + ,56613 + ,19 + ,234 + ,37 + ,1 + ,15 + ,12 + ,44 + ,35 + ,8019 + ,3080 + ,11017 + ,16 + ,16 + ,757 + ,62792 + ,35 + ,176 + ,46 + ,0 + ,28 + ,15 + ,60 + ,51 + ,30884 + ,10205 + ,37623 + ,72 + ,72 + ,894 + ,72535 + ,14 + ,329 + ,44 + ,0 + ,17 + ,16 + ,64 + ,39 + ,19540 + ,6095 + ,35873 + ,27 + ,21) + ,dim=c(15 + ,289) + ,dimnames=list(c('pageviews' + ,'time_in_rfc' + ,'logins' + ,'compendium_views_info' + ,'compendium_views_pr' + ,'shared_compendiums' + ,'blogged_computations' + ,'compendiums_reviewed' + ,'feedback_messages_p1' + ,'feedback_messages_p120' + ,'totsize' + ,'totrevisions' + ,'totseconds' + ,'tothyperlinks' + ,'totblogs ') + ,1:289)) > y <- array(NA,dim=c(15,289),dimnames=list(c('pageviews','time_in_rfc','logins','compendium_views_info','compendium_views_pr','shared_compendiums','blogged_computations','compendiums_reviewed','feedback_messages_p1','feedback_messages_p120','totsize','totrevisions','totseconds','tothyperlinks','totblogs '),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 = '2' > par3 <- 'Linear Trend' > par2 <- 'Do not include Seasonal Dummies' > par1 <- '2' > #'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 time_in_rfc pageviews logins compendium_views_info compendium_views_pr 1 210907 1418 56 396 81 2 120982 869 56 297 55 3 176508 1530 54 559 50 4 179321 2172 89 967 125 5 123185 901 40 270 40 6 52746 463 25 143 37 7 385534 3201 92 1562 63 8 33170 371 18 109 44 9 101645 1192 63 371 88 10 149061 1583 44 656 66 11 165446 1439 33 511 57 12 237213 1764 84 655 74 13 173326 1495 88 465 49 14 133131 1373 55 525 52 15 258873 2187 60 885 88 16 180083 1491 66 497 36 17 324799 4041 154 1436 108 18 230964 1706 53 612 43 19 236785 2152 119 865 75 20 135473 1036 41 385 32 21 202925 1882 61 567 44 22 215147 1929 58 639 85 23 344297 2242 75 963 86 24 153935 1220 33 398 56 25 132943 1289 40 410 50 26 174724 2515 92 966 135 27 174415 2147 100 801 63 28 225548 2352 112 892 81 29 223632 1638 73 513 52 30 124817 1222 40 469 44 31 221698 1812 45 683 113 32 210767 1677 60 643 39 33 170266 1579 62 535 73 34 260561 1731 75 625 48 35 84853 807 31 264 33 36 294424 2452 77 992 59 37 101011 829 34 238 41 38 215641 1940 46 818 69 39 325107 2662 99 937 64 40 7176 186 17 70 1 41 167542 1499 66 507 59 42 106408 865 30 260 32 43 96560 1793 76 503 129 44 265769 2527 146 927 37 45 269651 2747 67 1269 31 46 149112 1324 56 537 65 47 175824 2702 107 910 107 48 152871 1383 58 532 74 49 111665 1179 34 345 54 50 116408 2099 61 918 76 51 362301 4308 119 1635 715 52 78800 918 42 330 57 53 183167 1831 66 557 66 54 277965 3373 89 1178 106 55 150629 1713 44 740 54 56 168809 1438 66 452 32 57 24188 496 24 218 20 58 329267 2253 259 764 71 59 65029 744 17 255 21 60 101097 1161 64 454 70 61 218946 2352 41 866 112 62 244052 2144 68 574 66 63 341570 4691 168 1276 190 64 103597 1112 43 379 66 65 233328 2694 132 825 165 66 256462 1973 105 798 56 67 206161 1769 71 663 61 68 311473 3148 112 1069 53 69 235800 2474 94 921 127 70 177939 2084 82 858 63 71 207176 1954 70 711 38 72 196553 1226 57 503 50 73 174184 1389 53 382 52 74 143246 1496 103 464 42 75 187559 2269 121 717 76 76 187681 1833 62 690 67 77 119016 1268 52 462 50 78 182192 1943 52 657 53 79 73566 893 32 385 39 80 194979 1762 62 577 50 81 167488 1403 45 619 77 82 143756 1425 46 479 57 83 275541 1857 63 817 73 84 243199 1840 75 752 34 85 182999 1502 88 430 39 86 135649 1441 46 451 46 87 152299 1420 53 537 63 88 120221 1416 37 519 35 89 346485 2970 90 1000 106 90 145790 1317 63 637 43 91 193339 1644 78 465 47 92 80953 870 25 437 31 93 122774 1654 45 711 162 94 130585 1054 46 299 57 95 112611 937 41 248 36 96 286468 3004 144 1162 263 97 241066 2008 82 714 78 98 148446 2547 91 905 63 99 204713 1885 71 649 54 100 182079 1626 63 512 63 101 140344 1468 53 472 77 102 220516 2445 62 905 79 103 243060 1964 63 786 110 104 162765 1381 32 489 56 105 182613 1369 39 479 56 106 232138 1659 62 617 43 107 265318 2888 117 925 111 108 85574 1290 34 351 71 109 310839 2845 92 1144 62 110 225060 1982 93 669 56 111 232317 1904 54 707 74 112 144966 1391 144 458 60 113 43287 602 14 214 43 114 155754 1743 61 599 68 115 164709 1559 109 572 53 116 201940 2014 38 897 87 117 235454 2143 73 819 46 118 220801 2146 75 720 105 119 99466 874 50 273 32 120 92661 1590 61 508 133 121 133328 1590 55 506 79 122 61361 1210 77 451 51 123 125930 2072 75 699 207 124 100750 1281 72 407 67 125 224549 1401 50 465 47 126 82316 834 32 245 34 127 102010 1105 53 370 66 128 101523 1272 42 316 76 129 243511 1944 71 603 65 130 22938 391 10 154 9 131 41566 761 35 229 42 132 152474 1605 65 577 45 133 61857 530 25 192 25 134 99923 1988 66 617 115 135 132487 1386 41 411 97 136 317394 2395 86 975 53 137 21054 387 16 146 2 138 209641 1742 42 705 52 139 22648 620 19 184 44 140 31414 449 19 200 22 141 46698 800 45 274 35 142 131698 1684 65 502 74 143 91735 1050 35 382 103 144 244749 2699 95 964 144 145 184510 1606 49 537 60 146 79863 1502 37 438 134 147 128423 1204 64 369 89 148 97839 1138 38 417 42 149 38214 568 34 276 52 150 151101 1459 32 514 98 151 272458 2158 65 822 99 152 172494 1111 52 389 52 153 108043 1421 62 466 29 154 328107 2833 65 1255 125 155 250579 1955 83 694 106 156 351067 2922 95 1024 95 157 158015 1002 29 400 40 158 98866 1060 18 397 140 159 85439 956 33 350 43 160 229242 2186 247 719 128 161 351619 3604 139 1277 142 162 84207 1035 29 356 73 163 120445 1417 118 457 72 164 324598 3261 110 1402 128 165 131069 1587 67 600 61 166 204271 1424 42 480 73 167 165543 1701 65 595 148 168 141722 1249 94 436 64 169 116048 946 64 230 45 170 250047 1926 81 651 58 171 299775 3352 95 1367 97 172 195838 1641 67 564 50 173 173260 2035 63 716 37 174 254488 2312 83 747 50 175 104389 1369 45 467 105 176 136084 1577 30 671 69 177 199476 2201 70 861 46 178 92499 961 32 319 57 179 224330 1900 83 612 52 180 135781 1254 31 433 98 181 74408 1335 67 434 61 182 81240 1597 66 503 89 183 14688 207 10 85 0 184 181633 1645 70 564 48 185 271856 2429 103 824 91 186 7199 151 5 74 0 187 46660 474 20 259 7 188 17547 141 5 69 3 189 133368 1639 36 535 54 190 95227 872 34 239 70 191 152601 1318 48 438 36 192 98146 1018 40 459 37 193 79619 1383 43 426 123 194 59194 1314 31 288 247 195 139942 1335 42 498 46 196 118612 1403 46 454 72 197 72880 910 33 376 41 198 65475 616 18 225 24 199 99643 1407 55 555 45 200 71965 771 35 252 33 201 77272 766 59 208 27 202 49289 473 19 130 36 203 135131 1376 66 481 87 204 108446 1232 60 389 90 205 89746 1521 36 565 114 206 44296 572 25 173 31 207 77648 1059 47 278 45 208 181528 1544 54 609 69 209 134019 1230 53 422 51 210 124064 1206 40 445 34 211 92630 1205 40 387 60 212 121848 1255 39 339 45 213 52915 613 14 181 54 214 81872 721 45 245 25 215 58981 1109 36 384 38 216 53515 740 28 212 52 217 60812 1126 44 399 67 218 56375 728 30 229 74 219 65490 689 22 224 38 220 80949 592 17 203 30 221 76302 995 31 333 26 222 104011 1613 55 384 67 223 98104 2048 54 636 132 224 67989 705 21 185 42 225 30989 301 14 93 35 226 135458 1803 81 581 118 227 73504 799 35 248 68 228 63123 861 43 304 43 229 61254 1186 46 344 76 230 74914 1451 30 407 64 231 31774 628 23 170 48 232 81437 1161 38 312 64 233 87186 1463 54 507 56 234 50090 742 20 224 71 235 65745 979 53 340 75 236 56653 675 45 168 39 237 158399 1241 39 443 42 238 46455 676 20 204 39 239 73624 1049 24 367 93 240 38395 620 31 210 38 241 91899 1081 35 335 60 242 139526 1688 151 364 71 243 52164 736 52 178 52 244 51567 617 30 206 27 245 70551 812 31 279 59 246 84856 1051 29 387 40 247 102538 1656 57 490 79 248 86678 705 40 238 44 249 85709 945 44 343 65 250 34662 554 25 232 10 251 150580 1597 77 530 124 252 99611 982 35 291 81 253 19349 222 11 67 15 254 99373 1212 63 397 92 255 86230 1143 44 467 42 256 30837 435 19 178 10 257 31706 532 13 175 24 258 89806 882 42 299 64 259 62088 608 38 154 45 260 40151 459 29 106 22 261 27634 578 20 189 56 262 76990 826 27 194 94 263 37460 509 20 135 19 264 54157 717 19 201 35 265 49862 637 37 207 32 266 84337 857 26 280 35 267 64175 830 42 260 48 268 59382 652 49 227 49 269 119308 707 30 239 48 270 76702 954 49 333 62 271 103425 1461 67 428 96 272 70344 672 28 230 45 273 43410 778 19 292 63 274 104838 1141 49 350 71 275 62215 680 27 186 26 276 69304 1090 30 326 48 277 53117 616 22 155 29 278 19764 285 12 75 19 279 86680 1145 31 361 45 280 84105 733 20 261 45 281 77945 888 20 299 67 282 89113 849 39 300 30 283 91005 1182 29 450 36 284 40248 528 16 183 34 285 64187 642 27 238 36 286 50857 947 21 165 34 287 56613 819 19 234 37 288 62792 757 35 176 46 289 72535 894 14 329 44 shared_compendiums blogged_computations compendiums_reviewed 1 3 79 30 2 4 58 28 3 12 60 38 4 2 108 30 5 1 49 22 6 3 0 26 7 0 121 25 8 0 1 18 9 0 20 11 10 5 43 26 11 0 69 25 12 0 78 38 13 7 86 44 14 7 44 30 15 3 104 40 16 9 63 34 17 0 158 47 18 4 102 30 19 3 77 31 20 0 82 23 21 7 115 36 22 0 101 36 23 1 80 30 24 5 50 25 25 7 83 39 26 0 123 34 27 0 73 31 28 5 81 31 29 0 105 33 30 0 47 25 31 0 105 33 32 3 94 35 33 4 44 42 34 1 114 43 35 4 38 30 36 2 107 33 37 0 30 13 38 0 71 32 39 0 84 36 40 0 0 0 41 2 59 28 42 1 33 14 43 0 42 17 44 2 96 32 45 10 106 30 46 6 56 35 47 0 57 20 48 5 59 28 49 4 39 28 50 1 34 39 51 2 76 34 52 2 20 26 53 0 91 39 54 8 115 39 55 3 85 33 56 0 76 28 57 0 8 4 58 8 79 39 59 5 21 18 60 3 30 14 61 1 76 29 62 5 101 44 63 1 94 21 64 1 27 16 65 5 92 28 66 0 123 35 67 12 75 28 68 8 128 38 69 8 105 23 70 8 55 36 71 8 56 32 72 2 41 29 73 0 72 25 74 5 67 27 75 8 75 36 76 2 114 28 77 5 118 23 78 12 77 40 79 6 22 23 80 7 66 40 81 2 69 28 82 0 105 34 83 4 116 33 84 3 88 28 85 6 73 34 86 2 99 30 87 0 62 33 88 1 53 22 89 0 118 38 90 5 30 26 91 2 100 35 92 0 49 8 93 0 24 24 94 5 67 29 95 0 46 20 96 1 57 29 97 0 75 45 98 1 135 37 99 1 68 33 100 2 124 33 101 6 33 25 102 1 98 32 103 4 58 29 104 2 68 28 105 3 81 28 106 0 131 31 107 10 110 52 108 0 37 21 109 9 130 24 110 7 93 41 111 0 118 33 112 0 39 32 113 4 13 19 114 4 74 20 115 0 81 31 116 0 109 31 117 0 151 32 118 1 51 18 119 0 28 23 120 1 40 17 121 0 56 20 122 0 27 12 123 4 37 17 124 0 83 30 125 4 54 31 126 4 27 10 127 3 28 13 128 0 59 22 129 0 133 42 130 0 12 1 131 5 0 9 132 0 106 32 133 4 23 11 134 0 44 25 135 0 71 36 136 1 116 31 137 0 4 0 138 5 62 24 139 0 12 13 140 0 18 8 141 0 14 13 142 0 60 19 143 0 7 18 144 2 98 33 145 7 64 40 146 1 29 22 147 8 32 38 148 2 25 24 149 0 16 8 150 2 48 35 151 0 100 43 152 0 46 43 153 1 45 14 154 3 129 41 155 0 130 38 156 3 136 45 157 0 59 31 158 0 25 13 159 0 32 28 160 4 63 31 161 4 95 40 162 11 14 30 163 0 36 16 164 0 113 37 165 4 47 30 166 0 92 35 167 1 70 32 168 0 19 27 169 0 50 20 170 0 41 18 171 9 91 31 172 1 111 31 173 3 41 21 174 10 120 39 175 5 135 41 176 0 27 13 177 2 87 32 178 0 25 18 179 1 131 39 180 2 45 14 181 4 29 7 182 0 58 17 183 0 4 0 184 2 47 30 185 1 109 37 186 0 7 0 187 0 12 5 188 0 0 1 189 1 37 16 190 0 37 32 191 2 46 24 192 0 15 17 193 3 42 11 194 6 7 24 195 0 54 22 196 2 54 12 197 0 14 19 198 2 16 13 199 1 33 17 200 1 32 15 201 2 21 16 202 1 15 24 203 0 38 15 204 1 22 17 205 3 28 18 206 0 10 20 207 0 31 16 208 0 32 16 209 0 32 18 210 1 43 22 211 4 27 8 212 0 37 17 213 0 20 18 214 0 32 16 215 7 0 23 216 2 5 22 217 0 26 13 218 7 10 13 219 3 27 16 220 0 11 16 221 0 29 20 222 6 25 22 223 2 55 17 224 0 23 18 225 0 5 17 226 3 43 12 227 0 23 7 228 1 34 17 229 1 36 14 230 0 35 23 231 1 0 17 232 0 37 14 233 0 28 15 234 0 16 17 235 0 26 21 236 0 38 18 237 0 23 18 238 0 22 17 239 0 30 17 240 0 16 16 241 0 18 15 242 0 28 21 243 0 32 16 244 2 21 14 245 0 23 15 246 1 29 17 247 1 50 15 248 0 12 15 249 0 21 10 250 0 18 6 251 0 27 22 252 0 41 21 253 0 13 1 254 1 12 18 255 0 21 17 256 0 8 4 257 0 26 10 258 0 27 16 259 1 13 16 260 0 16 9 261 0 2 16 262 0 42 17 263 0 5 7 264 0 37 15 265 0 17 14 266 0 38 14 267 0 37 18 268 0 29 12 269 0 32 16 270 0 35 21 271 1 17 19 272 0 20 16 273 0 7 1 274 1 46 16 275 0 24 10 276 6 40 19 277 3 3 12 278 1 10 2 279 2 37 14 280 0 17 17 281 0 28 19 282 0 19 14 283 3 29 11 284 1 8 4 285 0 10 16 286 0 15 20 287 1 15 12 288 0 28 15 289 0 17 16 feedback_messages_p1 feedback_messages_p120 totsize totrevisions totseconds 1 115 94 112285 24188 146283 2 109 103 84786 18273 98364 3 146 93 83123 14130 86146 4 116 103 101193 32287 96933 5 68 51 38361 8654 79234 6 101 70 68504 9245 42551 7 96 91 119182 33251 195663 8 67 22 22807 1271 6853 9 44 38 17140 5279 21529 10 100 93 116174 27101 95757 11 93 60 57635 16373 85584 12 140 123 66198 19716 143983 13 166 148 71701 17753 75851 14 99 90 57793 9028 59238 15 139 124 80444 18653 93163 16 130 70 53855 8828 96037 17 181 168 97668 29498 151511 18 116 115 133824 27563 136368 19 116 71 101481 18293 112642 20 88 66 99645 22530 94728 21 139 134 114789 15977 105499 22 135 117 99052 35082 121527 23 108 108 67654 16116 127766 24 89 84 65553 15849 98958 25 156 156 97500 16026 77900 26 129 120 69112 26569 85646 27 118 114 82753 24785 98579 28 118 94 85323 17569 130767 29 125 120 72654 23825 131741 30 95 81 30727 7869 53907 31 126 110 77873 14975 178812 32 135 133 117478 37791 146761 33 154 122 74007 9605 82036 34 165 158 90183 27295 163253 35 113 109 61542 2746 27032 36 127 124 101494 34461 171975 37 52 39 27570 8098 65990 38 121 92 55813 4787 86572 39 136 126 79215 24919 159676 40 0 0 1423 603 1929 41 108 70 55461 16329 85371 42 46 37 31081 12558 58391 43 54 38 22996 7784 31580 44 124 120 83122 28522 136815 45 115 93 70106 22265 120642 46 128 95 60578 14459 69107 47 80 77 39992 14526 50495 48 97 90 79892 22240 108016 49 104 80 49810 11802 46341 50 59 31 71570 7623 78348 51 125 110 100708 11912 79336 52 82 66 33032 7935 56968 53 149 138 82875 18220 93176 54 149 133 139077 19199 161632 55 122 113 71595 19918 87850 56 118 100 72260 21884 127969 57 12 7 5950 2694 15049 58 144 140 115762 15808 155135 59 67 61 32551 3597 25109 60 52 41 31701 5296 45824 61 108 96 80670 25239 102996 62 166 164 143558 29801 160604 63 80 78 117105 18450 158051 64 60 49 23789 7132 44547 65 107 102 120733 34861 162647 66 127 124 105195 35940 174141 67 107 99 73107 16688 60622 68 146 129 132068 24683 179566 69 84 62 149193 46230 184301 70 141 73 46821 10387 75661 71 123 114 87011 21436 96144 72 111 99 95260 30546 129847 73 98 70 55183 19746 117286 74 105 104 106671 15977 71180 75 135 116 73511 22583 109377 76 107 91 92945 17274 85298 77 85 74 78664 16469 73631 78 155 138 70054 14251 86767 79 88 67 22618 3007 23824 80 155 151 74011 16851 93487 81 104 72 83737 21113 82981 82 132 120 69094 17401 73815 83 127 115 93133 23958 94552 84 108 105 95536 23567 132190 85 129 104 225920 13065 128754 86 116 108 62133 15358 66363 87 122 98 61370 14587 67808 88 85 69 43836 12770 61724 89 147 111 106117 24021 131722 90 99 99 38692 9648 68580 91 87 71 84651 20537 106175 92 28 27 56622 7905 55792 93 90 69 15986 4527 25157 94 109 107 95364 30495 76669 95 78 73 26706 7117 57283 96 111 107 89691 17719 105805 97 158 93 67267 27056 129484 98 141 129 126846 33473 72413 99 122 69 41140 9758 87831 100 124 118 102860 21115 96971 101 93 73 51715 7236 71299 102 124 119 55801 13790 77494 103 112 104 111813 32902 120336 104 108 107 120293 25131 93913 105 99 99 138599 30910 136048 106 117 90 161647 35947 181248 107 199 197 115929 29848 146123 108 78 36 24266 6943 32036 109 91 85 162901 42705 186646 110 158 139 109825 31808 102255 111 126 106 129838 26675 168237 112 122 50 37510 8435 64219 113 71 64 43750 7409 19630 114 75 31 40652 14993 76825 115 115 63 87771 36867 115338 116 119 92 85872 33835 109427 117 124 106 89275 24164 118168 118 72 63 44418 12607 84845 119 91 69 192565 22609 153197 120 45 41 35232 5892 29877 121 78 56 40909 17014 63506 122 39 25 13294 5394 22445 123 68 65 32387 9178 47695 124 119 93 140867 6440 68370 125 117 114 120662 21916 146304 126 39 38 21233 4011 38233 127 50 44 44332 5818 42071 128 88 87 61056 18647 50517 129 155 110 101338 20556 103950 130 0 0 1168 238 5841 131 36 27 13497 70 2341 132 123 83 65567 22392 84396 133 32 30 25162 3913 24610 134 99 80 32334 12237 35753 135 136 98 40735 8388 55515 136 117 82 91413 22120 209056 137 0 0 855 338 6622 138 88 60 97068 11727 115814 139 39 28 44339 3704 11609 140 25 9 14116 3988 13155 141 52 33 10288 3030 18274 142 75 59 65622 13520 72875 143 71 49 16563 1421 10112 144 124 115 76643 20923 142775 145 151 140 110681 20237 68847 146 71 49 29011 3219 17659 147 145 120 92696 3769 20112 148 87 66 94785 12252 61023 149 27 21 8773 1888 13983 150 131 124 83209 14497 65176 151 162 152 93815 28864 132432 152 165 139 86687 21721 112494 153 54 38 34553 4821 45109 154 159 144 105547 33644 170875 155 147 120 103487 15923 180759 156 170 160 213688 42935 214921 157 119 114 71220 18864 100226 158 49 39 23517 4977 32043 159 104 78 56926 7785 54454 160 120 119 91721 17939 78876 161 150 141 115168 23436 170745 162 112 101 111194 325 6940 163 59 56 51009 13539 49025 164 136 133 135777 34538 122037 165 107 83 51513 12198 53782 166 130 116 74163 26924 127748 167 115 90 51633 12716 86839 168 107 36 75345 8172 44830 169 75 50 33416 10855 77395 170 71 61 83305 11932 89324 171 120 97 98952 14300 103300 172 116 98 102372 25515 112283 173 79 78 37238 2805 10901 174 150 117 103772 29402 120691 175 156 148 123969 16440 58106 176 51 41 27142 11221 57140 177 118 105 135400 28732 122422 178 71 55 21399 5250 25899 179 144 132 130115 28608 139296 180 47 44 24874 8092 52678 181 28 21 34988 4473 23853 182 68 50 45549 1572 17306 183 0 0 6023 2065 7953 184 110 73 64466 14817 89455 185 147 86 54990 16714 147866 186 0 0 1644 556 4245 187 15 13 6179 2089 21509 188 4 4 3926 2658 7670 189 64 57 32755 10695 66675 190 111 48 34777 1669 14336 191 85 46 73224 16267 53608 192 68 48 27114 7768 30059 193 40 32 20760 7252 29668 194 80 68 37636 6387 22097 195 88 87 65461 18715 96841 196 48 43 30080 7936 41907 197 76 67 24094 8643 27080 198 51 46 69008 7294 35885 199 67 46 54968 4570 41247 200 59 56 46090 7185 28313 201 61 48 27507 10058 36845 202 76 44 10672 2342 16548 203 60 60 34029 8509 36134 204 68 65 46300 13275 55764 205 71 55 24760 6816 28910 206 76 38 18779 1930 13339 207 62 52 21280 8086 25319 208 61 60 40662 10737 66956 209 67 54 28987 8033 47487 210 88 86 22827 7058 52785 211 30 24 18513 6782 44683 212 64 52 30594 5401 35619 213 68 49 24006 6521 21920 214 64 61 27913 10856 45608 215 91 61 42744 2154 7721 216 88 81 12934 6117 20634 217 52 43 22574 5238 29788 218 49 40 41385 4820 31931 219 62 40 18653 5615 37754 220 61 56 18472 4272 32505 221 76 68 30976 8702 40557 222 88 79 63339 15340 94238 223 66 47 25568 8030 44197 224 71 57 33747 9526 43228 225 68 41 4154 1278 4103 226 48 29 19474 4236 44144 227 25 3 35130 3023 32868 228 68 60 39067 7196 27640 229 41 30 13310 3394 14063 230 90 79 65892 6371 28990 231 66 47 4143 1574 4694 232 54 40 28579 9620 42648 233 59 48 51776 6978 64329 234 60 36 21152 4911 21928 235 77 42 38084 8645 25836 236 68 49 27717 8987 22779 237 72 57 32928 5544 40820 238 67 12 11342 3083 27530 239 64 40 19499 6909 32378 240 63 43 16380 3189 10824 241 59 33 36874 6745 39613 242 84 77 48259 16724 60865 243 64 43 16734 4850 19787 244 56 45 28207 7025 20107 245 54 47 30143 6047 36605 246 67 43 41369 7377 40961 247 58 45 45833 9078 48231 248 59 50 29156 4605 39725 249 40 35 35944 3238 21455 250 22 7 36278 8100 23430 251 83 71 45588 9653 62991 252 81 67 45097 8914 49363 253 2 0 3895 786 9604 254 72 62 28394 6700 24552 255 61 54 18632 5788 31493 256 15 4 2325 593 3439 257 32 25 25139 4506 19555 258 62 40 27975 6382 21228 259 58 38 14483 5621 23177 260 36 19 13127 3997 22094 261 59 17 5839 520 2342 262 68 67 24069 8891 38798 263 21 14 3738 999 3255 264 55 30 18625 7067 24261 265 54 54 36341 4639 18511 266 55 35 24548 5654 40798 267 72 59 21792 6928 28893 268 41 24 26263 1514 21425 269 61 58 23686 9238 50276 270 67 42 49303 8204 37643 271 76 46 25659 5926 30377 272 64 61 28904 5785 27126 273 3 3 2781 4 13 274 63 52 29236 5930 42097 275 40 25 19546 3710 24451 276 69 40 22818 705 14335 277 48 32 32689 443 5084 278 8 4 5752 2416 9927 279 52 49 22197 7747 43527 280 66 63 20055 5432 27184 281 76 67 25272 4913 21610 282 43 32 82206 2650 20484 283 39 23 32073 2370 20156 284 14 7 5444 775 6012 285 61 54 20154 5576 18475 286 71 37 36944 1352 12645 287 44 35 8019 3080 11017 288 60 51 30884 10205 37623 289 64 39 19540 6095 35873 tothyperlinks totblogs\r t 1 144 145 1 2 103 101 2 3 98 98 3 4 135 132 4 5 61 60 5 6 39 38 6 7 150 144 7 8 5 5 8 9 28 28 9 10 84 84 10 11 80 79 11 12 130 127 12 13 82 78 13 14 60 60 14 15 131 131 15 16 84 84 16 17 140 133 17 18 151 150 18 19 91 91 19 20 138 132 20 21 150 136 21 22 124 124 22 23 119 118 23 24 73 70 24 25 110 107 25 26 123 119 26 27 90 89 27 28 116 112 28 29 113 108 29 30 56 52 30 31 115 112 31 32 119 116 32 33 129 123 33 34 127 125 34 35 27 27 35 36 175 162 36 37 35 32 37 38 64 64 38 39 96 92 39 40 0 0 40 41 84 83 41 42 41 41 42 43 47 47 43 44 126 120 44 45 105 105 45 46 80 79 46 47 70 65 47 48 73 70 48 49 57 55 49 50 40 39 50 51 68 67 51 52 21 21 52 53 127 127 53 54 154 152 54 55 116 113 55 56 102 99 56 57 7 7 57 58 148 141 58 59 21 21 59 60 35 35 60 61 112 109 61 62 137 133 62 63 135 123 63 64 26 26 64 65 230 230 65 66 181 166 66 67 71 68 67 68 147 147 68 69 190 179 69 70 64 61 70 71 105 101 71 72 107 108 72 73 94 90 73 74 116 114 74 75 106 103 75 76 143 142 76 77 81 79 77 78 89 88 78 79 26 25 79 80 84 83 80 81 113 113 81 82 120 118 82 83 110 110 83 84 134 129 84 85 54 51 85 86 96 93 86 87 78 76 87 88 51 49 88 89 121 118 89 90 38 38 90 91 145 141 91 92 59 58 92 93 27 27 93 94 91 91 94 95 48 48 95 96 68 63 96 97 58 56 97 98 150 144 98 99 74 73 99 100 181 168 100 101 65 64 101 102 97 97 102 103 121 117 103 104 99 100 104 105 152 149 105 106 188 187 106 107 138 127 107 108 40 37 108 109 254 245 109 110 87 87 110 111 178 177 111 112 51 49 112 113 49 49 113 114 73 73 114 115 176 177 115 116 94 94 116 117 120 117 117 118 66 60 118 119 56 55 119 120 39 39 120 121 66 64 121 122 27 26 122 123 65 64 123 124 58 58 124 125 98 95 125 126 25 25 126 127 26 26 127 128 77 76 128 129 130 129 129 130 11 11 130 131 2 2 131 132 101 101 132 133 31 28 133 134 36 36 134 135 120 89 135 136 195 193 136 137 4 4 137 138 89 84 138 139 24 23 139 140 39 39 140 141 14 14 141 142 78 78 142 143 15 14 143 144 106 101 144 145 83 82 145 146 24 24 146 147 37 36 147 148 77 75 148 149 16 16 149 150 56 55 150 151 132 131 151 152 144 131 152 153 40 39 153 154 153 144 154 155 143 139 155 156 220 211 156 157 79 78 157 158 50 50 158 159 39 39 159 160 95 90 160 161 169 166 161 162 12 12 162 163 63 57 163 164 134 133 164 165 69 69 165 166 119 119 166 167 119 119 167 168 75 65 168 169 63 61 169 170 55 49 170 171 103 101 171 172 197 196 172 173 16 15 173 174 140 136 174 175 89 89 175 176 40 40 176 177 125 123 177 178 21 21 178 179 167 163 179 180 32 29 180 181 36 35 181 182 13 13 182 183 5 5 183 184 96 96 184 185 151 151 185 186 6 6 186 187 13 13 187 188 3 3 188 189 57 56 189 190 23 23 190 191 61 57 191 192 21 14 192 193 43 43 193 194 20 20 194 195 82 72 195 196 90 87 196 197 25 21 197 198 60 56 198 199 61 59 199 200 85 82 200 201 43 43 201 202 25 25 202 203 41 38 203 204 26 25 204 205 38 38 205 206 12 12 206 207 29 29 207 208 49 47 208 209 46 45 209 210 41 40 210 211 31 30 211 212 41 41 212 213 26 25 213 214 23 23 214 215 14 14 215 216 16 16 216 217 25 26 217 218 21 21 218 219 32 27 219 220 9 9 220 221 35 33 221 222 42 42 222 223 68 68 223 224 32 32 224 225 6 6 225 226 68 67 226 227 33 33 227 228 84 77 228 229 46 46 229 230 30 30 230 231 0 0 231 232 36 36 232 233 47 46 233 234 20 18 234 235 50 48 235 236 30 29 236 237 30 28 237 238 34 34 238 239 33 33 239 240 34 34 240 241 37 33 241 242 83 80 242 243 32 32 243 244 30 30 244 245 43 41 245 246 41 41 246 247 51 51 247 248 19 18 248 249 37 34 249 250 33 31 250 251 41 39 251 252 54 54 252 253 14 14 253 254 25 24 254 255 25 24 255 256 8 8 256 257 26 26 257 258 20 19 258 259 11 11 259 260 14 14 260 261 3 1 261 262 40 39 262 263 5 5 263 264 38 37 264 265 32 32 265 266 41 38 266 267 46 47 267 268 47 47 268 269 37 37 269 270 51 51 270 271 49 45 271 272 21 21 272 273 1 1 273 274 44 42 274 275 26 26 275 276 21 21 276 277 4 4 277 278 10 10 278 279 43 43 279 280 34 34 280 281 32 31 281 282 20 19 282 283 34 34 283 284 6 6 284 285 12 11 285 286 24 24 286 287 16 16 287 288 72 72 288 289 27 21 289 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) pageviews logins -9437.5600 5.0350 176.1018 compendium_views_info compendium_views_pr shared_compendiums 113.7962 16.0210 -114.9774 blogged_computations compendiums_reviewed feedback_messages_p1 43.1976 -815.9969 353.6426 feedback_messages_p120 totsize totrevisions 191.0420 -0.1022 -0.8045 totseconds tothyperlinks `totblogs\\r` 0.8692 96.0070 -88.3276 t 17.5342 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -61576 -11517 470 9777 88013 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -9.438e+03 5.726e+03 -1.648 0.10045 pageviews 5.035e+00 8.183e+00 0.615 0.53886 logins 1.761e+02 6.106e+01 2.884 0.00424 ** compendium_views_info 1.138e+02 1.749e+01 6.506 3.68e-10 *** compendium_views_pr 1.602e+01 3.315e+01 0.483 0.62923 shared_compendiums -1.150e+02 5.609e+02 -0.205 0.83773 blogged_computations 4.320e+01 9.196e+01 0.470 0.63890 compendiums_reviewed -8.160e+02 7.012e+02 -1.164 0.24557 feedback_messages_p1 3.536e+02 2.107e+02 1.678 0.09446 . feedback_messages_p120 1.910e+02 1.041e+02 1.836 0.06744 . totsize -1.022e-01 6.485e-02 -1.576 0.11626 totrevisions -8.045e-01 3.188e-01 -2.523 0.01219 * totseconds 8.692e-01 7.197e-02 12.077 < 2e-16 *** tothyperlinks 9.601e+01 4.369e+02 0.220 0.82624 `totblogs\\r` -8.833e+01 4.528e+02 -0.195 0.84548 t 1.753e+01 1.978e+01 0.886 0.37617 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 21540 on 273 degrees of freedom Multiple R-squared: 0.9351, Adjusted R-squared: 0.9315 F-statistic: 262.3 on 15 and 273 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.5945930 8.108140e-01 4.054070e-01 [2,] 0.6928307 6.143387e-01 3.071693e-01 [3,] 0.6723134 6.553732e-01 3.276866e-01 [4,] 0.6293667 7.412665e-01 3.706333e-01 [5,] 0.6761007 6.477985e-01 3.238993e-01 [6,] 0.7259328 5.481344e-01 2.740672e-01 [7,] 0.9247819 1.504363e-01 7.521813e-02 [8,] 0.9648675 7.026493e-02 3.513246e-02 [9,] 0.9719675 5.606501e-02 2.803251e-02 [10,] 0.9748995 5.020102e-02 2.510051e-02 [11,] 0.9701531 5.969379e-02 2.984690e-02 [12,] 0.9604659 7.906811e-02 3.953405e-02 [13,] 0.9841232 3.175354e-02 1.587677e-02 [14,] 0.9776001 4.479979e-02 2.239990e-02 [15,] 0.9674973 6.500534e-02 3.250267e-02 [16,] 0.9551825 8.963497e-02 4.481748e-02 [17,] 0.9377356 1.245288e-01 6.226441e-02 [18,] 0.9226967 1.546066e-01 7.730332e-02 [19,] 0.8982521 2.034958e-01 1.017479e-01 [20,] 0.8694283 2.611433e-01 1.305717e-01 [21,] 0.9152384 1.695232e-01 8.476158e-02 [22,] 0.9272772 1.454456e-01 7.272281e-02 [23,] 0.9106299 1.787403e-01 8.937015e-02 [24,] 0.8916148 2.167704e-01 1.083852e-01 [25,] 0.8635880 2.728240e-01 1.364120e-01 [26,] 0.8327008 3.345984e-01 1.672992e-01 [27,] 0.8493373 3.013255e-01 1.506627e-01 [28,] 0.8277481 3.445039e-01 1.722519e-01 [29,] 0.7950757 4.098486e-01 2.049243e-01 [30,] 0.7591701 4.816598e-01 2.408299e-01 [31,] 0.7291583 5.416833e-01 2.708417e-01 [32,] 0.8866641 2.266718e-01 1.133359e-01 [33,] 0.9896963 2.060742e-02 1.030371e-02 [34,] 0.9891674 2.166515e-02 1.083258e-02 [35,] 0.9863964 2.720730e-02 1.360365e-02 [36,] 0.9953420 9.316068e-03 4.658034e-03 [37,] 0.9969305 6.138906e-03 3.069453e-03 [38,] 0.9963216 7.356716e-03 3.678358e-03 [39,] 0.9951891 9.621887e-03 4.810943e-03 [40,] 0.9971203 5.759332e-03 2.879666e-03 [41,] 0.9961207 7.758563e-03 3.879282e-03 [42,] 0.9948744 1.025120e-02 5.125598e-03 [43,] 0.9943097 1.138052e-02 5.690262e-03 [44,] 0.9933171 1.336590e-02 6.682948e-03 [45,] 0.9910023 1.799532e-02 8.997662e-03 [46,] 0.9883458 2.330833e-02 1.165417e-02 [47,] 0.9941381 1.172385e-02 5.861923e-03 [48,] 0.9928888 1.422232e-02 7.111162e-03 [49,] 0.9985463 2.907456e-03 1.453728e-03 [50,] 0.9980806 3.838757e-03 1.919379e-03 [51,] 0.9978790 4.241936e-03 2.120968e-03 [52,] 0.9976260 4.747961e-03 2.373980e-03 [53,] 0.9972470 5.505968e-03 2.752984e-03 [54,] 0.9970060 5.988091e-03 2.994046e-03 [55,] 0.9965250 6.950026e-03 3.475013e-03 [56,] 0.9953595 9.280971e-03 4.640486e-03 [57,] 0.9965338 6.932498e-03 3.466249e-03 [58,] 0.9962155 7.568927e-03 3.784464e-03 [59,] 0.9952036 9.592836e-03 4.796418e-03 [60,] 0.9940881 1.182384e-02 5.911921e-03 [61,] 0.9929688 1.406231e-02 7.031155e-03 [62,] 0.9909139 1.817218e-02 9.086088e-03 [63,] 0.9905850 1.883006e-02 9.415032e-03 [64,] 0.9878834 2.423325e-02 1.211662e-02 [65,] 0.9994115 1.177000e-03 5.884999e-04 [66,] 0.9992882 1.423509e-03 7.117547e-04 [67,] 0.9990162 1.967598e-03 9.837992e-04 [68,] 0.9986532 2.693697e-03 1.346848e-03 [69,] 0.9981997 3.600645e-03 1.800323e-03 [70,] 0.9977800 4.440002e-03 2.220001e-03 [71,] 0.9999540 9.209414e-05 4.604707e-05 [72,] 0.9999560 8.806449e-05 4.403225e-05 [73,] 0.9999753 4.932813e-05 2.466406e-05 [74,] 0.9999742 5.167618e-05 2.583809e-05 [75,] 0.9999650 6.994792e-05 3.497396e-05 [76,] 0.9999614 7.712401e-05 3.856200e-05 [77,] 0.9999478 1.043973e-04 5.219867e-05 [78,] 0.9999393 1.214381e-04 6.071907e-05 [79,] 0.9999205 1.589247e-04 7.946235e-05 [80,] 0.9999866 2.685454e-05 1.342727e-05 [81,] 0.9999832 3.355333e-05 1.677667e-05 [82,] 0.9999767 4.656121e-05 2.328060e-05 [83,] 0.9999662 6.751914e-05 3.375957e-05 [84,] 0.9999510 9.806400e-05 4.903200e-05 [85,] 0.9999736 5.272040e-05 2.636020e-05 [86,] 0.9999667 6.650877e-05 3.325439e-05 [87,] 0.9999523 9.541491e-05 4.770745e-05 [88,] 0.9999389 1.221891e-04 6.109456e-05 [89,] 0.9999709 5.819394e-05 2.909697e-05 [90,] 0.9999576 8.482397e-05 4.241198e-05 [91,] 0.9999433 1.134182e-04 5.670909e-05 [92,] 0.9999414 1.172839e-04 5.864196e-05 [93,] 0.9999402 1.195110e-04 5.975551e-05 [94,] 0.9999189 1.622885e-04 8.114425e-05 [95,] 0.9998901 2.197294e-04 1.098647e-04 [96,] 0.9998467 3.066992e-04 1.533496e-04 [97,] 0.9998167 3.666012e-04 1.833006e-04 [98,] 0.9997584 4.831591e-04 2.415795e-04 [99,] 0.9996651 6.697366e-04 3.348683e-04 [100,] 0.9998257 3.486022e-04 1.743011e-04 [101,] 0.9999797 4.063184e-05 2.031592e-05 [102,] 0.9999737 5.250391e-05 2.625195e-05 [103,] 0.9999622 7.568146e-05 3.784073e-05 [104,] 0.9999757 4.855460e-05 2.427730e-05 [105,] 0.9999776 4.479320e-05 2.239660e-05 [106,] 0.9999872 2.569959e-05 1.284980e-05 [107,] 0.9999876 2.479780e-05 1.239890e-05 [108,] 0.9999832 3.363904e-05 1.681952e-05 [109,] 0.9999767 4.656882e-05 2.328441e-05 [110,] 0.9999672 6.567311e-05 3.283655e-05 [111,] 0.9999918 1.631611e-05 8.158053e-06 [112,] 0.9999880 2.402893e-05 1.201446e-05 [113,] 0.9999824 3.525202e-05 1.762601e-05 [114,] 0.9999779 4.412535e-05 2.206268e-05 [115,] 0.9999712 5.765323e-05 2.882661e-05 [116,] 0.9999839 3.229934e-05 1.614967e-05 [117,] 0.9999790 4.206251e-05 2.103126e-05 [118,] 0.9999724 5.526001e-05 2.763000e-05 [119,] 0.9999599 8.022989e-05 4.011494e-05 [120,] 0.9999482 1.035866e-04 5.179332e-05 [121,] 0.9999390 1.219381e-04 6.096906e-05 [122,] 0.9999140 1.720550e-04 8.602749e-05 [123,] 0.9999033 1.934339e-04 9.671697e-05 [124,] 0.9998725 2.549487e-04 1.274743e-04 [125,] 0.9998706 2.587771e-04 1.293885e-04 [126,] 0.9999307 1.386037e-04 6.930186e-05 [127,] 0.9999532 9.350015e-05 4.675008e-05 [128,] 0.9999394 1.211609e-04 6.058047e-05 [129,] 0.9999614 7.714605e-05 3.857302e-05 [130,] 0.9999563 8.738028e-05 4.369014e-05 [131,] 0.9999476 1.047621e-04 5.238105e-05 [132,] 0.9999274 1.451658e-04 7.258289e-05 [133,] 0.9999398 1.203454e-04 6.017269e-05 [134,] 0.9999188 1.623144e-04 8.115721e-05 [135,] 0.9998974 2.052795e-04 1.026398e-04 [136,] 0.9998552 2.896632e-04 1.448316e-04 [137,] 0.9999115 1.770012e-04 8.850060e-05 [138,] 0.9998892 2.215504e-04 1.107752e-04 [139,] 0.9998452 3.096480e-04 1.548240e-04 [140,] 0.9998633 2.734449e-04 1.367224e-04 [141,] 0.9998941 2.118414e-04 1.059207e-04 [142,] 0.9998618 2.763073e-04 1.381536e-04 [143,] 0.9998126 3.748796e-04 1.874398e-04 [144,] 0.9997969 4.061178e-04 2.030589e-04 [145,] 0.9997859 4.281296e-04 2.140648e-04 [146,] 0.9997915 4.169595e-04 2.084797e-04 [147,] 0.9997774 4.452176e-04 2.226088e-04 [148,] 0.9997746 4.507338e-04 2.253669e-04 [149,] 0.9997075 5.850409e-04 2.925204e-04 [150,] 0.9996995 6.009013e-04 3.004506e-04 [151,] 0.9996091 7.818442e-04 3.909221e-04 [152,] 0.9999881 2.387762e-05 1.193881e-05 [153,] 0.9999839 3.228610e-05 1.614305e-05 [154,] 0.9999831 3.382452e-05 1.691226e-05 [155,] 0.9999951 9.750658e-06 4.875329e-06 [156,] 0.9999989 2.203865e-06 1.101933e-06 [157,] 0.9999992 1.673482e-06 8.367408e-07 [158,] 0.9999988 2.442779e-06 1.221389e-06 [159,] 0.9999988 2.492523e-06 1.246262e-06 [160,] 0.9999986 2.888761e-06 1.444381e-06 [161,] 0.9999978 4.309722e-06 2.154861e-06 [162,] 0.9999983 3.359902e-06 1.679951e-06 [163,] 0.9999976 4.898303e-06 2.449152e-06 [164,] 0.9999984 3.270713e-06 1.635357e-06 [165,] 0.9999974 5.226990e-06 2.613495e-06 [166,] 0.9999971 5.818967e-06 2.909484e-06 [167,] 0.9999966 6.859039e-06 3.429520e-06 [168,] 0.9999949 1.025190e-05 5.125950e-06 [169,] 0.9999932 1.363827e-05 6.819136e-06 [170,] 0.9999895 2.105661e-05 1.052830e-05 [171,] 0.9999856 2.876579e-05 1.438290e-05 [172,] 0.9999887 2.259638e-05 1.129819e-05 [173,] 0.9999997 6.838541e-07 3.419270e-07 [174,] 0.9999994 1.144207e-06 5.721033e-07 [175,] 0.9999991 1.818696e-06 9.093479e-07 [176,] 0.9999987 2.679696e-06 1.339848e-06 [177,] 0.9999986 2.756822e-06 1.378411e-06 [178,] 0.9999988 2.428762e-06 1.214381e-06 [179,] 0.9999985 3.080173e-06 1.540087e-06 [180,] 0.9999977 4.558693e-06 2.279346e-06 [181,] 0.9999977 4.574836e-06 2.287418e-06 [182,] 0.9999969 6.147345e-06 3.073672e-06 [183,] 0.9999959 8.275675e-06 4.137837e-06 [184,] 0.9999937 1.253449e-05 6.267246e-06 [185,] 0.9999937 1.251406e-05 6.257032e-06 [186,] 0.9999899 2.012585e-05 1.006292e-05 [187,] 0.9999858 2.842073e-05 1.421037e-05 [188,] 0.9999774 4.525136e-05 2.262568e-05 [189,] 0.9999669 6.621095e-05 3.310548e-05 [190,] 0.9999936 1.281516e-05 6.407580e-06 [191,] 0.9999974 5.127375e-06 2.563687e-06 [192,] 0.9999957 8.587271e-06 4.293635e-06 [193,] 0.9999944 1.118507e-05 5.592534e-06 [194,] 0.9999995 9.538656e-07 4.769328e-07 [195,] 0.9999994 1.260679e-06 6.303393e-07 [196,] 0.9999989 2.232061e-06 1.116031e-06 [197,] 0.9999982 3.546227e-06 1.773114e-06 [198,] 0.9999970 5.902187e-06 2.951093e-06 [199,] 0.9999980 3.946971e-06 1.973485e-06 [200,] 0.9999966 6.741647e-06 3.370823e-06 [201,] 0.9999942 1.167817e-05 5.839086e-06 [202,] 0.9999943 1.146948e-05 5.734738e-06 [203,] 0.9999906 1.881856e-05 9.409282e-06 [204,] 0.9999958 8.322586e-06 4.161293e-06 [205,] 0.9999960 7.904181e-06 3.952091e-06 [206,] 0.9999932 1.365265e-05 6.826325e-06 [207,] 0.9999886 2.287598e-05 1.143799e-05 [208,] 0.9999856 2.887745e-05 1.443872e-05 [209,] 0.9999845 3.100145e-05 1.550073e-05 [210,] 0.9999742 5.169788e-05 2.584894e-05 [211,] 0.9999629 7.422507e-05 3.711254e-05 [212,] 0.9999662 6.762887e-05 3.381444e-05 [213,] 0.9999529 9.420772e-05 4.710386e-05 [214,] 0.9999187 1.626947e-04 8.134733e-05 [215,] 0.9999995 1.090066e-06 5.450332e-07 [216,] 0.9999989 2.213571e-06 1.106786e-06 [217,] 0.9999979 4.285060e-06 2.142530e-06 [218,] 0.9999960 7.960374e-06 3.980187e-06 [219,] 1.0000000 5.375424e-09 2.687712e-09 [220,] 1.0000000 1.315458e-08 6.577290e-09 [221,] 1.0000000 3.297627e-08 1.648814e-08 [222,] 1.0000000 8.140718e-08 4.070359e-08 [223,] 0.9999999 1.251924e-07 6.259622e-08 [224,] 0.9999999 2.497451e-07 1.248725e-07 [225,] 0.9999997 5.811973e-07 2.905987e-07 [226,] 0.9999993 1.333100e-06 6.665498e-07 [227,] 0.9999985 3.095040e-06 1.547520e-06 [228,] 0.9999964 7.111829e-06 3.555914e-06 [229,] 0.9999962 7.649377e-06 3.824689e-06 [230,] 0.9999916 1.687519e-05 8.437597e-06 [231,] 0.9999905 1.906537e-05 9.532685e-06 [232,] 0.9999814 3.717656e-05 1.858828e-05 [233,] 0.9999679 6.416286e-05 3.208143e-05 [234,] 0.9999306 1.387833e-04 6.939164e-05 [235,] 0.9998510 2.980473e-04 1.490237e-04 [236,] 0.9996782 6.436175e-04 3.218087e-04 [237,] 0.9996264 7.471962e-04 3.735981e-04 [238,] 0.9992255 1.548912e-03 7.744560e-04 [239,] 0.9988516 2.296738e-03 1.148369e-03 [240,] 0.9992037 1.592514e-03 7.962570e-04 [241,] 0.9983364 3.327151e-03 1.663575e-03 [242,] 0.9971044 5.791246e-03 2.895623e-03 [243,] 0.9950211 9.957867e-03 4.978933e-03 [244,] 0.9900493 1.990130e-02 9.950651e-03 [245,] 0.9855996 2.880089e-02 1.440045e-02 [246,] 0.9906743 1.865144e-02 9.325720e-03 [247,] 0.9803535 3.929307e-02 1.964653e-02 [248,] 0.9773676 4.526487e-02 2.263244e-02 [249,] 0.9501348 9.973034e-02 4.986517e-02 [250,] 0.9066215 1.867569e-01 9.337846e-02 [251,] 0.9984126 3.174726e-03 1.587363e-03 [252,] 0.9979280 4.143998e-03 2.071999e-03 > postscript(file="/var/wessaorg/rcomp/tmp/1njjr1353319398.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/2ng2p1353319398.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/3l9311353319398.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/4dxte1353319398.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/5fyv61353319398.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 6 22524.08798 -19049.18590 8899.09502 -39854.42987 13206.85550 -11984.14960 7 8 9 10 11 12 14796.92312 8426.92362 22564.05449 -16498.80835 19612.63513 -2105.98905 13 14 15 16 17 18 6526.37742 -2450.82089 41245.29140 8176.92208 -42914.47444 25449.52181 19 20 21 22 23 24 8675.89699 3681.18388 7506.35581 16882.60258 88013.27010 8717.10707 25 26 27 28 29 30 -20937.19821 -50662.53788 -36630.66115 -29238.26452 18289.10880 -2410.62399 31 32 33 34 35 36 -45728.67244 -7796.15791 -4263.74252 5461.37370 892.30655 3558.82619 37 38 39 40 41 42 6813.03300 8849.75097 37220.84746 2954.13911 11284.01392 24362.15295 43 44 45 46 47 48 -9969.24708 -3633.05932 -10205.67033 -463.89236 -11516.66994 -16709.95685 49 50 51 52 53 54 11448.89544 -53387.11028 40465.16067 -24100.01077 -3604.42757 -44460.52171 55 56 57 58 59 60 -37646.98069 -22098.81437 -12246.53486 32272.34594 581.13945 -8764.51219 61 62 63 64 65 66 9776.54070 2986.78711 5790.93305 5230.83005 -32334.14867 -15150.07400 67 68 69 70 71 72 48220.85006 -12076.16438 -28399.57656 -26696.91091 12434.71482 15151.26792 73 74 75 76 77 78 9481.52709 -4226.95675 -33171.37253 7046.03962 -15271.52337 -12338.44537 79 80 81 82 83 84 -14626.71725 2826.14338 10958.35261 -9533.76448 70199.14420 16069.52283 85 86 87 88 89 90 -868.22506 -5035.06975 2288.85480 -12077.50330 75272.30587 -18864.55971 91 92 93 94 95 96 37062.02088 -17767.48760 -12435.69098 19327.77499 10747.71909 9410.47726 97 98 99 100 101 102 17067.16970 -51439.72103 18539.93300 7293.15395 -2281.39230 4790.18276 103 104 105 106 107 108 32667.77200 12966.66728 6165.21891 -3014.39711 -31660.45522 280.40158 109 110 111 112 113 114 9280.39748 24267.87964 -16856.81887 -7037.17649 -8733.96120 4441.86213 115 116 117 118 119 120 -14642.64764 -9864.17043 5416.15513 35544.62465 -61575.74854 -8107.01265 121 122 123 124 125 126 1713.20407 -27180.19204 -28369.41916 -37753.88418 24319.40113 9661.20863 127 128 129 130 131 132 7262.04483 2054.82236 46890.34154 4144.84539 1394.20751 -13547.25460 133 134 135 136 137 138 14862.93086 -36672.92369 -4474.01954 -7437.53279 1077.15167 11897.99816 139 140 141 142 143 144 -10366.22564 -2213.88238 -17222.27944 -10829.37749 15853.72729 -34245.23200 145 146 147 148 149 150 30181.29242 -6961.25585 25527.46630 -14771.35931 -13633.47699 4565.51747 151 152 153 154 155 156 22851.61465 -6084.96544 -5722.86663 -3827.74821 -32525.76376 15876.61081 157 158 159 160 161 162 1920.89741 12528.04862 -24816.40522 9442.05596 -2313.83565 10541.36685 163 164 165 166 167 168 -1629.98308 20058.40247 -14414.80862 14184.50708 -13582.05310 25075.34978 169 170 171 172 173 174 1898.53014 73209.76794 7790.25124 9883.88179 45499.37568 29133.51239 175 176 177 178 179 180 -39419.91658 -2846.75659 -26108.53500 12858.43162 -285.91820 26884.50823 181 182 183 184 185 186 -11131.01660 -22298.22538 3594.12524 12138.16918 -973.38802 -109.21924 187 188 189 190 191 192 -3384.31312 8685.13950 -5321.80677 30985.20165 45314.73302 1343.14002 193 194 195 196 197 198 -10223.97705 -15261.43738 -21623.73455 11209.40855 -14417.92816 3420.04840 199 200 201 202 203 204 -22178.49163 3018.41689 5021.38160 6255.13766 22582.99680 -4944.21080 205 206 207 208 209 210 -23727.96049 -3807.03199 3798.80589 32311.76639 22431.44483 -4092.21601 211 212 213 214 215 216 -931.25398 30988.15242 378.48706 -3121.89250 -16561.05587 -14892.50688 217 218 219 220 221 222 -30043.89900 -8000.08089 -6765.81908 14287.10699 -17245.01097 -45089.25971 223 224 225 226 227 228 -41135.41710 -4482.78207 972.62061 1768.17170 11969.21276 -17124.96723 229 230 231 232 233 234 -6491.74951 -22842.94523 -10888.11295 -5857.53479 -47390.29679 -6653.12974 235 236 237 238 239 240 -14503.60331 469.87415 48761.84692 -12660.04755 -12433.07969 -13021.36899 241 242 243 244 245 246 4841.84047 3163.70469 -7617.15344 -2364.36904 -7755.60678 -10027.14952 247 248 249 250 251 252 -16994.86809 6612.27787 11702.37811 -9606.36323 2770.59968 1490.71903 253 254 255 256 257 258 5553.50598 4158.74174 -16138.56857 4347.51710 -8687.18215 19455.14547 259 260 261 262 263 264 9597.37430 735.18240 -8705.10687 43.32399 14247.43703 157.94071 265 266 267 268 269 270 -6750.79997 2456.46238 -14221.28108 -34.11583 32111.67627 -8479.04069 271 272 273 274 275 276 -1809.43715 1578.99116 5707.70370 3756.30154 8801.71997 -2034.53180 277 278 279 280 281 282 18785.13605 3488.07624 -9155.92588 11291.44095 1400.23558 29613.59491 283 284 285 286 287 288 9923.49655 10513.78497 4022.78901 4675.30911 6938.83049 -5378.62402 289 -10983.73711 > postscript(file="/var/wessaorg/rcomp/tmp/6rqgk1353319398.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 22524.08798 NA 1 -19049.18590 22524.08798 2 8899.09502 -19049.18590 3 -39854.42987 8899.09502 4 13206.85550 -39854.42987 5 -11984.14960 13206.85550 6 14796.92312 -11984.14960 7 8426.92362 14796.92312 8 22564.05449 8426.92362 9 -16498.80835 22564.05449 10 19612.63513 -16498.80835 11 -2105.98905 19612.63513 12 6526.37742 -2105.98905 13 -2450.82089 6526.37742 14 41245.29140 -2450.82089 15 8176.92208 41245.29140 16 -42914.47444 8176.92208 17 25449.52181 -42914.47444 18 8675.89699 25449.52181 19 3681.18388 8675.89699 20 7506.35581 3681.18388 21 16882.60258 7506.35581 22 88013.27010 16882.60258 23 8717.10707 88013.27010 24 -20937.19821 8717.10707 25 -50662.53788 -20937.19821 26 -36630.66115 -50662.53788 27 -29238.26452 -36630.66115 28 18289.10880 -29238.26452 29 -2410.62399 18289.10880 30 -45728.67244 -2410.62399 31 -7796.15791 -45728.67244 32 -4263.74252 -7796.15791 33 5461.37370 -4263.74252 34 892.30655 5461.37370 35 3558.82619 892.30655 36 6813.03300 3558.82619 37 8849.75097 6813.03300 38 37220.84746 8849.75097 39 2954.13911 37220.84746 40 11284.01392 2954.13911 41 24362.15295 11284.01392 42 -9969.24708 24362.15295 43 -3633.05932 -9969.24708 44 -10205.67033 -3633.05932 45 -463.89236 -10205.67033 46 -11516.66994 -463.89236 47 -16709.95685 -11516.66994 48 11448.89544 -16709.95685 49 -53387.11028 11448.89544 50 40465.16067 -53387.11028 51 -24100.01077 40465.16067 52 -3604.42757 -24100.01077 53 -44460.52171 -3604.42757 54 -37646.98069 -44460.52171 55 -22098.81437 -37646.98069 56 -12246.53486 -22098.81437 57 32272.34594 -12246.53486 58 581.13945 32272.34594 59 -8764.51219 581.13945 60 9776.54070 -8764.51219 61 2986.78711 9776.54070 62 5790.93305 2986.78711 63 5230.83005 5790.93305 64 -32334.14867 5230.83005 65 -15150.07400 -32334.14867 66 48220.85006 -15150.07400 67 -12076.16438 48220.85006 68 -28399.57656 -12076.16438 69 -26696.91091 -28399.57656 70 12434.71482 -26696.91091 71 15151.26792 12434.71482 72 9481.52709 15151.26792 73 -4226.95675 9481.52709 74 -33171.37253 -4226.95675 75 7046.03962 -33171.37253 76 -15271.52337 7046.03962 77 -12338.44537 -15271.52337 78 -14626.71725 -12338.44537 79 2826.14338 -14626.71725 80 10958.35261 2826.14338 81 -9533.76448 10958.35261 82 70199.14420 -9533.76448 83 16069.52283 70199.14420 84 -868.22506 16069.52283 85 -5035.06975 -868.22506 86 2288.85480 -5035.06975 87 -12077.50330 2288.85480 88 75272.30587 -12077.50330 89 -18864.55971 75272.30587 90 37062.02088 -18864.55971 91 -17767.48760 37062.02088 92 -12435.69098 -17767.48760 93 19327.77499 -12435.69098 94 10747.71909 19327.77499 95 9410.47726 10747.71909 96 17067.16970 9410.47726 97 -51439.72103 17067.16970 98 18539.93300 -51439.72103 99 7293.15395 18539.93300 100 -2281.39230 7293.15395 101 4790.18276 -2281.39230 102 32667.77200 4790.18276 103 12966.66728 32667.77200 104 6165.21891 12966.66728 105 -3014.39711 6165.21891 106 -31660.45522 -3014.39711 107 280.40158 -31660.45522 108 9280.39748 280.40158 109 24267.87964 9280.39748 110 -16856.81887 24267.87964 111 -7037.17649 -16856.81887 112 -8733.96120 -7037.17649 113 4441.86213 -8733.96120 114 -14642.64764 4441.86213 115 -9864.17043 -14642.64764 116 5416.15513 -9864.17043 117 35544.62465 5416.15513 118 -61575.74854 35544.62465 119 -8107.01265 -61575.74854 120 1713.20407 -8107.01265 121 -27180.19204 1713.20407 122 -28369.41916 -27180.19204 123 -37753.88418 -28369.41916 124 24319.40113 -37753.88418 125 9661.20863 24319.40113 126 7262.04483 9661.20863 127 2054.82236 7262.04483 128 46890.34154 2054.82236 129 4144.84539 46890.34154 130 1394.20751 4144.84539 131 -13547.25460 1394.20751 132 14862.93086 -13547.25460 133 -36672.92369 14862.93086 134 -4474.01954 -36672.92369 135 -7437.53279 -4474.01954 136 1077.15167 -7437.53279 137 11897.99816 1077.15167 138 -10366.22564 11897.99816 139 -2213.88238 -10366.22564 140 -17222.27944 -2213.88238 141 -10829.37749 -17222.27944 142 15853.72729 -10829.37749 143 -34245.23200 15853.72729 144 30181.29242 -34245.23200 145 -6961.25585 30181.29242 146 25527.46630 -6961.25585 147 -14771.35931 25527.46630 148 -13633.47699 -14771.35931 149 4565.51747 -13633.47699 150 22851.61465 4565.51747 151 -6084.96544 22851.61465 152 -5722.86663 -6084.96544 153 -3827.74821 -5722.86663 154 -32525.76376 -3827.74821 155 15876.61081 -32525.76376 156 1920.89741 15876.61081 157 12528.04862 1920.89741 158 -24816.40522 12528.04862 159 9442.05596 -24816.40522 160 -2313.83565 9442.05596 161 10541.36685 -2313.83565 162 -1629.98308 10541.36685 163 20058.40247 -1629.98308 164 -14414.80862 20058.40247 165 14184.50708 -14414.80862 166 -13582.05310 14184.50708 167 25075.34978 -13582.05310 168 1898.53014 25075.34978 169 73209.76794 1898.53014 170 7790.25124 73209.76794 171 9883.88179 7790.25124 172 45499.37568 9883.88179 173 29133.51239 45499.37568 174 -39419.91658 29133.51239 175 -2846.75659 -39419.91658 176 -26108.53500 -2846.75659 177 12858.43162 -26108.53500 178 -285.91820 12858.43162 179 26884.50823 -285.91820 180 -11131.01660 26884.50823 181 -22298.22538 -11131.01660 182 3594.12524 -22298.22538 183 12138.16918 3594.12524 184 -973.38802 12138.16918 185 -109.21924 -973.38802 186 -3384.31312 -109.21924 187 8685.13950 -3384.31312 188 -5321.80677 8685.13950 189 30985.20165 -5321.80677 190 45314.73302 30985.20165 191 1343.14002 45314.73302 192 -10223.97705 1343.14002 193 -15261.43738 -10223.97705 194 -21623.73455 -15261.43738 195 11209.40855 -21623.73455 196 -14417.92816 11209.40855 197 3420.04840 -14417.92816 198 -22178.49163 3420.04840 199 3018.41689 -22178.49163 200 5021.38160 3018.41689 201 6255.13766 5021.38160 202 22582.99680 6255.13766 203 -4944.21080 22582.99680 204 -23727.96049 -4944.21080 205 -3807.03199 -23727.96049 206 3798.80589 -3807.03199 207 32311.76639 3798.80589 208 22431.44483 32311.76639 209 -4092.21601 22431.44483 210 -931.25398 -4092.21601 211 30988.15242 -931.25398 212 378.48706 30988.15242 213 -3121.89250 378.48706 214 -16561.05587 -3121.89250 215 -14892.50688 -16561.05587 216 -30043.89900 -14892.50688 217 -8000.08089 -30043.89900 218 -6765.81908 -8000.08089 219 14287.10699 -6765.81908 220 -17245.01097 14287.10699 221 -45089.25971 -17245.01097 222 -41135.41710 -45089.25971 223 -4482.78207 -41135.41710 224 972.62061 -4482.78207 225 1768.17170 972.62061 226 11969.21276 1768.17170 227 -17124.96723 11969.21276 228 -6491.74951 -17124.96723 229 -22842.94523 -6491.74951 230 -10888.11295 -22842.94523 231 -5857.53479 -10888.11295 232 -47390.29679 -5857.53479 233 -6653.12974 -47390.29679 234 -14503.60331 -6653.12974 235 469.87415 -14503.60331 236 48761.84692 469.87415 237 -12660.04755 48761.84692 238 -12433.07969 -12660.04755 239 -13021.36899 -12433.07969 240 4841.84047 -13021.36899 241 3163.70469 4841.84047 242 -7617.15344 3163.70469 243 -2364.36904 -7617.15344 244 -7755.60678 -2364.36904 245 -10027.14952 -7755.60678 246 -16994.86809 -10027.14952 247 6612.27787 -16994.86809 248 11702.37811 6612.27787 249 -9606.36323 11702.37811 250 2770.59968 -9606.36323 251 1490.71903 2770.59968 252 5553.50598 1490.71903 253 4158.74174 5553.50598 254 -16138.56857 4158.74174 255 4347.51710 -16138.56857 256 -8687.18215 4347.51710 257 19455.14547 -8687.18215 258 9597.37430 19455.14547 259 735.18240 9597.37430 260 -8705.10687 735.18240 261 43.32399 -8705.10687 262 14247.43703 43.32399 263 157.94071 14247.43703 264 -6750.79997 157.94071 265 2456.46238 -6750.79997 266 -14221.28108 2456.46238 267 -34.11583 -14221.28108 268 32111.67627 -34.11583 269 -8479.04069 32111.67627 270 -1809.43715 -8479.04069 271 1578.99116 -1809.43715 272 5707.70370 1578.99116 273 3756.30154 5707.70370 274 8801.71997 3756.30154 275 -2034.53180 8801.71997 276 18785.13605 -2034.53180 277 3488.07624 18785.13605 278 -9155.92588 3488.07624 279 11291.44095 -9155.92588 280 1400.23558 11291.44095 281 29613.59491 1400.23558 282 9923.49655 29613.59491 283 10513.78497 9923.49655 284 4022.78901 10513.78497 285 4675.30911 4022.78901 286 6938.83049 4675.30911 287 -5378.62402 6938.83049 288 -10983.73711 -5378.62402 289 NA -10983.73711 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -19049.18590 22524.08798 [2,] 8899.09502 -19049.18590 [3,] -39854.42987 8899.09502 [4,] 13206.85550 -39854.42987 [5,] -11984.14960 13206.85550 [6,] 14796.92312 -11984.14960 [7,] 8426.92362 14796.92312 [8,] 22564.05449 8426.92362 [9,] -16498.80835 22564.05449 [10,] 19612.63513 -16498.80835 [11,] -2105.98905 19612.63513 [12,] 6526.37742 -2105.98905 [13,] -2450.82089 6526.37742 [14,] 41245.29140 -2450.82089 [15,] 8176.92208 41245.29140 [16,] -42914.47444 8176.92208 [17,] 25449.52181 -42914.47444 [18,] 8675.89699 25449.52181 [19,] 3681.18388 8675.89699 [20,] 7506.35581 3681.18388 [21,] 16882.60258 7506.35581 [22,] 88013.27010 16882.60258 [23,] 8717.10707 88013.27010 [24,] -20937.19821 8717.10707 [25,] -50662.53788 -20937.19821 [26,] -36630.66115 -50662.53788 [27,] -29238.26452 -36630.66115 [28,] 18289.10880 -29238.26452 [29,] -2410.62399 18289.10880 [30,] -45728.67244 -2410.62399 [31,] -7796.15791 -45728.67244 [32,] -4263.74252 -7796.15791 [33,] 5461.37370 -4263.74252 [34,] 892.30655 5461.37370 [35,] 3558.82619 892.30655 [36,] 6813.03300 3558.82619 [37,] 8849.75097 6813.03300 [38,] 37220.84746 8849.75097 [39,] 2954.13911 37220.84746 [40,] 11284.01392 2954.13911 [41,] 24362.15295 11284.01392 [42,] -9969.24708 24362.15295 [43,] -3633.05932 -9969.24708 [44,] -10205.67033 -3633.05932 [45,] -463.89236 -10205.67033 [46,] -11516.66994 -463.89236 [47,] -16709.95685 -11516.66994 [48,] 11448.89544 -16709.95685 [49,] -53387.11028 11448.89544 [50,] 40465.16067 -53387.11028 [51,] -24100.01077 40465.16067 [52,] -3604.42757 -24100.01077 [53,] -44460.52171 -3604.42757 [54,] -37646.98069 -44460.52171 [55,] -22098.81437 -37646.98069 [56,] -12246.53486 -22098.81437 [57,] 32272.34594 -12246.53486 [58,] 581.13945 32272.34594 [59,] -8764.51219 581.13945 [60,] 9776.54070 -8764.51219 [61,] 2986.78711 9776.54070 [62,] 5790.93305 2986.78711 [63,] 5230.83005 5790.93305 [64,] -32334.14867 5230.83005 [65,] -15150.07400 -32334.14867 [66,] 48220.85006 -15150.07400 [67,] -12076.16438 48220.85006 [68,] -28399.57656 -12076.16438 [69,] -26696.91091 -28399.57656 [70,] 12434.71482 -26696.91091 [71,] 15151.26792 12434.71482 [72,] 9481.52709 15151.26792 [73,] -4226.95675 9481.52709 [74,] -33171.37253 -4226.95675 [75,] 7046.03962 -33171.37253 [76,] -15271.52337 7046.03962 [77,] -12338.44537 -15271.52337 [78,] -14626.71725 -12338.44537 [79,] 2826.14338 -14626.71725 [80,] 10958.35261 2826.14338 [81,] -9533.76448 10958.35261 [82,] 70199.14420 -9533.76448 [83,] 16069.52283 70199.14420 [84,] -868.22506 16069.52283 [85,] -5035.06975 -868.22506 [86,] 2288.85480 -5035.06975 [87,] -12077.50330 2288.85480 [88,] 75272.30587 -12077.50330 [89,] -18864.55971 75272.30587 [90,] 37062.02088 -18864.55971 [91,] -17767.48760 37062.02088 [92,] -12435.69098 -17767.48760 [93,] 19327.77499 -12435.69098 [94,] 10747.71909 19327.77499 [95,] 9410.47726 10747.71909 [96,] 17067.16970 9410.47726 [97,] -51439.72103 17067.16970 [98,] 18539.93300 -51439.72103 [99,] 7293.15395 18539.93300 [100,] -2281.39230 7293.15395 [101,] 4790.18276 -2281.39230 [102,] 32667.77200 4790.18276 [103,] 12966.66728 32667.77200 [104,] 6165.21891 12966.66728 [105,] -3014.39711 6165.21891 [106,] -31660.45522 -3014.39711 [107,] 280.40158 -31660.45522 [108,] 9280.39748 280.40158 [109,] 24267.87964 9280.39748 [110,] -16856.81887 24267.87964 [111,] -7037.17649 -16856.81887 [112,] -8733.96120 -7037.17649 [113,] 4441.86213 -8733.96120 [114,] -14642.64764 4441.86213 [115,] -9864.17043 -14642.64764 [116,] 5416.15513 -9864.17043 [117,] 35544.62465 5416.15513 [118,] -61575.74854 35544.62465 [119,] -8107.01265 -61575.74854 [120,] 1713.20407 -8107.01265 [121,] -27180.19204 1713.20407 [122,] -28369.41916 -27180.19204 [123,] -37753.88418 -28369.41916 [124,] 24319.40113 -37753.88418 [125,] 9661.20863 24319.40113 [126,] 7262.04483 9661.20863 [127,] 2054.82236 7262.04483 [128,] 46890.34154 2054.82236 [129,] 4144.84539 46890.34154 [130,] 1394.20751 4144.84539 [131,] -13547.25460 1394.20751 [132,] 14862.93086 -13547.25460 [133,] -36672.92369 14862.93086 [134,] -4474.01954 -36672.92369 [135,] -7437.53279 -4474.01954 [136,] 1077.15167 -7437.53279 [137,] 11897.99816 1077.15167 [138,] -10366.22564 11897.99816 [139,] -2213.88238 -10366.22564 [140,] -17222.27944 -2213.88238 [141,] -10829.37749 -17222.27944 [142,] 15853.72729 -10829.37749 [143,] -34245.23200 15853.72729 [144,] 30181.29242 -34245.23200 [145,] -6961.25585 30181.29242 [146,] 25527.46630 -6961.25585 [147,] -14771.35931 25527.46630 [148,] -13633.47699 -14771.35931 [149,] 4565.51747 -13633.47699 [150,] 22851.61465 4565.51747 [151,] -6084.96544 22851.61465 [152,] -5722.86663 -6084.96544 [153,] -3827.74821 -5722.86663 [154,] -32525.76376 -3827.74821 [155,] 15876.61081 -32525.76376 [156,] 1920.89741 15876.61081 [157,] 12528.04862 1920.89741 [158,] -24816.40522 12528.04862 [159,] 9442.05596 -24816.40522 [160,] -2313.83565 9442.05596 [161,] 10541.36685 -2313.83565 [162,] -1629.98308 10541.36685 [163,] 20058.40247 -1629.98308 [164,] -14414.80862 20058.40247 [165,] 14184.50708 -14414.80862 [166,] -13582.05310 14184.50708 [167,] 25075.34978 -13582.05310 [168,] 1898.53014 25075.34978 [169,] 73209.76794 1898.53014 [170,] 7790.25124 73209.76794 [171,] 9883.88179 7790.25124 [172,] 45499.37568 9883.88179 [173,] 29133.51239 45499.37568 [174,] -39419.91658 29133.51239 [175,] -2846.75659 -39419.91658 [176,] -26108.53500 -2846.75659 [177,] 12858.43162 -26108.53500 [178,] -285.91820 12858.43162 [179,] 26884.50823 -285.91820 [180,] -11131.01660 26884.50823 [181,] -22298.22538 -11131.01660 [182,] 3594.12524 -22298.22538 [183,] 12138.16918 3594.12524 [184,] -973.38802 12138.16918 [185,] -109.21924 -973.38802 [186,] -3384.31312 -109.21924 [187,] 8685.13950 -3384.31312 [188,] -5321.80677 8685.13950 [189,] 30985.20165 -5321.80677 [190,] 45314.73302 30985.20165 [191,] 1343.14002 45314.73302 [192,] -10223.97705 1343.14002 [193,] -15261.43738 -10223.97705 [194,] -21623.73455 -15261.43738 [195,] 11209.40855 -21623.73455 [196,] -14417.92816 11209.40855 [197,] 3420.04840 -14417.92816 [198,] -22178.49163 3420.04840 [199,] 3018.41689 -22178.49163 [200,] 5021.38160 3018.41689 [201,] 6255.13766 5021.38160 [202,] 22582.99680 6255.13766 [203,] -4944.21080 22582.99680 [204,] -23727.96049 -4944.21080 [205,] -3807.03199 -23727.96049 [206,] 3798.80589 -3807.03199 [207,] 32311.76639 3798.80589 [208,] 22431.44483 32311.76639 [209,] -4092.21601 22431.44483 [210,] -931.25398 -4092.21601 [211,] 30988.15242 -931.25398 [212,] 378.48706 30988.15242 [213,] -3121.89250 378.48706 [214,] -16561.05587 -3121.89250 [215,] -14892.50688 -16561.05587 [216,] -30043.89900 -14892.50688 [217,] -8000.08089 -30043.89900 [218,] -6765.81908 -8000.08089 [219,] 14287.10699 -6765.81908 [220,] -17245.01097 14287.10699 [221,] -45089.25971 -17245.01097 [222,] -41135.41710 -45089.25971 [223,] -4482.78207 -41135.41710 [224,] 972.62061 -4482.78207 [225,] 1768.17170 972.62061 [226,] 11969.21276 1768.17170 [227,] -17124.96723 11969.21276 [228,] -6491.74951 -17124.96723 [229,] -22842.94523 -6491.74951 [230,] -10888.11295 -22842.94523 [231,] -5857.53479 -10888.11295 [232,] -47390.29679 -5857.53479 [233,] -6653.12974 -47390.29679 [234,] -14503.60331 -6653.12974 [235,] 469.87415 -14503.60331 [236,] 48761.84692 469.87415 [237,] -12660.04755 48761.84692 [238,] -12433.07969 -12660.04755 [239,] -13021.36899 -12433.07969 [240,] 4841.84047 -13021.36899 [241,] 3163.70469 4841.84047 [242,] -7617.15344 3163.70469 [243,] -2364.36904 -7617.15344 [244,] -7755.60678 -2364.36904 [245,] -10027.14952 -7755.60678 [246,] -16994.86809 -10027.14952 [247,] 6612.27787 -16994.86809 [248,] 11702.37811 6612.27787 [249,] -9606.36323 11702.37811 [250,] 2770.59968 -9606.36323 [251,] 1490.71903 2770.59968 [252,] 5553.50598 1490.71903 [253,] 4158.74174 5553.50598 [254,] -16138.56857 4158.74174 [255,] 4347.51710 -16138.56857 [256,] -8687.18215 4347.51710 [257,] 19455.14547 -8687.18215 [258,] 9597.37430 19455.14547 [259,] 735.18240 9597.37430 [260,] -8705.10687 735.18240 [261,] 43.32399 -8705.10687 [262,] 14247.43703 43.32399 [263,] 157.94071 14247.43703 [264,] -6750.79997 157.94071 [265,] 2456.46238 -6750.79997 [266,] -14221.28108 2456.46238 [267,] -34.11583 -14221.28108 [268,] 32111.67627 -34.11583 [269,] -8479.04069 32111.67627 [270,] -1809.43715 -8479.04069 [271,] 1578.99116 -1809.43715 [272,] 5707.70370 1578.99116 [273,] 3756.30154 5707.70370 [274,] 8801.71997 3756.30154 [275,] -2034.53180 8801.71997 [276,] 18785.13605 -2034.53180 [277,] 3488.07624 18785.13605 [278,] -9155.92588 3488.07624 [279,] 11291.44095 -9155.92588 [280,] 1400.23558 11291.44095 [281,] 29613.59491 1400.23558 [282,] 9923.49655 29613.59491 [283,] 10513.78497 9923.49655 [284,] 4022.78901 10513.78497 [285,] 4675.30911 4022.78901 [286,] 6938.83049 4675.30911 [287,] -5378.62402 6938.83049 [288,] -10983.73711 -5378.62402 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -19049.18590 22524.08798 2 8899.09502 -19049.18590 3 -39854.42987 8899.09502 4 13206.85550 -39854.42987 5 -11984.14960 13206.85550 6 14796.92312 -11984.14960 7 8426.92362 14796.92312 8 22564.05449 8426.92362 9 -16498.80835 22564.05449 10 19612.63513 -16498.80835 11 -2105.98905 19612.63513 12 6526.37742 -2105.98905 13 -2450.82089 6526.37742 14 41245.29140 -2450.82089 15 8176.92208 41245.29140 16 -42914.47444 8176.92208 17 25449.52181 -42914.47444 18 8675.89699 25449.52181 19 3681.18388 8675.89699 20 7506.35581 3681.18388 21 16882.60258 7506.35581 22 88013.27010 16882.60258 23 8717.10707 88013.27010 24 -20937.19821 8717.10707 25 -50662.53788 -20937.19821 26 -36630.66115 -50662.53788 27 -29238.26452 -36630.66115 28 18289.10880 -29238.26452 29 -2410.62399 18289.10880 30 -45728.67244 -2410.62399 31 -7796.15791 -45728.67244 32 -4263.74252 -7796.15791 33 5461.37370 -4263.74252 34 892.30655 5461.37370 35 3558.82619 892.30655 36 6813.03300 3558.82619 37 8849.75097 6813.03300 38 37220.84746 8849.75097 39 2954.13911 37220.84746 40 11284.01392 2954.13911 41 24362.15295 11284.01392 42 -9969.24708 24362.15295 43 -3633.05932 -9969.24708 44 -10205.67033 -3633.05932 45 -463.89236 -10205.67033 46 -11516.66994 -463.89236 47 -16709.95685 -11516.66994 48 11448.89544 -16709.95685 49 -53387.11028 11448.89544 50 40465.16067 -53387.11028 51 -24100.01077 40465.16067 52 -3604.42757 -24100.01077 53 -44460.52171 -3604.42757 54 -37646.98069 -44460.52171 55 -22098.81437 -37646.98069 56 -12246.53486 -22098.81437 57 32272.34594 -12246.53486 58 581.13945 32272.34594 59 -8764.51219 581.13945 60 9776.54070 -8764.51219 61 2986.78711 9776.54070 62 5790.93305 2986.78711 63 5230.83005 5790.93305 64 -32334.14867 5230.83005 65 -15150.07400 -32334.14867 66 48220.85006 -15150.07400 67 -12076.16438 48220.85006 68 -28399.57656 -12076.16438 69 -26696.91091 -28399.57656 70 12434.71482 -26696.91091 71 15151.26792 12434.71482 72 9481.52709 15151.26792 73 -4226.95675 9481.52709 74 -33171.37253 -4226.95675 75 7046.03962 -33171.37253 76 -15271.52337 7046.03962 77 -12338.44537 -15271.52337 78 -14626.71725 -12338.44537 79 2826.14338 -14626.71725 80 10958.35261 2826.14338 81 -9533.76448 10958.35261 82 70199.14420 -9533.76448 83 16069.52283 70199.14420 84 -868.22506 16069.52283 85 -5035.06975 -868.22506 86 2288.85480 -5035.06975 87 -12077.50330 2288.85480 88 75272.30587 -12077.50330 89 -18864.55971 75272.30587 90 37062.02088 -18864.55971 91 -17767.48760 37062.02088 92 -12435.69098 -17767.48760 93 19327.77499 -12435.69098 94 10747.71909 19327.77499 95 9410.47726 10747.71909 96 17067.16970 9410.47726 97 -51439.72103 17067.16970 98 18539.93300 -51439.72103 99 7293.15395 18539.93300 100 -2281.39230 7293.15395 101 4790.18276 -2281.39230 102 32667.77200 4790.18276 103 12966.66728 32667.77200 104 6165.21891 12966.66728 105 -3014.39711 6165.21891 106 -31660.45522 -3014.39711 107 280.40158 -31660.45522 108 9280.39748 280.40158 109 24267.87964 9280.39748 110 -16856.81887 24267.87964 111 -7037.17649 -16856.81887 112 -8733.96120 -7037.17649 113 4441.86213 -8733.96120 114 -14642.64764 4441.86213 115 -9864.17043 -14642.64764 116 5416.15513 -9864.17043 117 35544.62465 5416.15513 118 -61575.74854 35544.62465 119 -8107.01265 -61575.74854 120 1713.20407 -8107.01265 121 -27180.19204 1713.20407 122 -28369.41916 -27180.19204 123 -37753.88418 -28369.41916 124 24319.40113 -37753.88418 125 9661.20863 24319.40113 126 7262.04483 9661.20863 127 2054.82236 7262.04483 128 46890.34154 2054.82236 129 4144.84539 46890.34154 130 1394.20751 4144.84539 131 -13547.25460 1394.20751 132 14862.93086 -13547.25460 133 -36672.92369 14862.93086 134 -4474.01954 -36672.92369 135 -7437.53279 -4474.01954 136 1077.15167 -7437.53279 137 11897.99816 1077.15167 138 -10366.22564 11897.99816 139 -2213.88238 -10366.22564 140 -17222.27944 -2213.88238 141 -10829.37749 -17222.27944 142 15853.72729 -10829.37749 143 -34245.23200 15853.72729 144 30181.29242 -34245.23200 145 -6961.25585 30181.29242 146 25527.46630 -6961.25585 147 -14771.35931 25527.46630 148 -13633.47699 -14771.35931 149 4565.51747 -13633.47699 150 22851.61465 4565.51747 151 -6084.96544 22851.61465 152 -5722.86663 -6084.96544 153 -3827.74821 -5722.86663 154 -32525.76376 -3827.74821 155 15876.61081 -32525.76376 156 1920.89741 15876.61081 157 12528.04862 1920.89741 158 -24816.40522 12528.04862 159 9442.05596 -24816.40522 160 -2313.83565 9442.05596 161 10541.36685 -2313.83565 162 -1629.98308 10541.36685 163 20058.40247 -1629.98308 164 -14414.80862 20058.40247 165 14184.50708 -14414.80862 166 -13582.05310 14184.50708 167 25075.34978 -13582.05310 168 1898.53014 25075.34978 169 73209.76794 1898.53014 170 7790.25124 73209.76794 171 9883.88179 7790.25124 172 45499.37568 9883.88179 173 29133.51239 45499.37568 174 -39419.91658 29133.51239 175 -2846.75659 -39419.91658 176 -26108.53500 -2846.75659 177 12858.43162 -26108.53500 178 -285.91820 12858.43162 179 26884.50823 -285.91820 180 -11131.01660 26884.50823 181 -22298.22538 -11131.01660 182 3594.12524 -22298.22538 183 12138.16918 3594.12524 184 -973.38802 12138.16918 185 -109.21924 -973.38802 186 -3384.31312 -109.21924 187 8685.13950 -3384.31312 188 -5321.80677 8685.13950 189 30985.20165 -5321.80677 190 45314.73302 30985.20165 191 1343.14002 45314.73302 192 -10223.97705 1343.14002 193 -15261.43738 -10223.97705 194 -21623.73455 -15261.43738 195 11209.40855 -21623.73455 196 -14417.92816 11209.40855 197 3420.04840 -14417.92816 198 -22178.49163 3420.04840 199 3018.41689 -22178.49163 200 5021.38160 3018.41689 201 6255.13766 5021.38160 202 22582.99680 6255.13766 203 -4944.21080 22582.99680 204 -23727.96049 -4944.21080 205 -3807.03199 -23727.96049 206 3798.80589 -3807.03199 207 32311.76639 3798.80589 208 22431.44483 32311.76639 209 -4092.21601 22431.44483 210 -931.25398 -4092.21601 211 30988.15242 -931.25398 212 378.48706 30988.15242 213 -3121.89250 378.48706 214 -16561.05587 -3121.89250 215 -14892.50688 -16561.05587 216 -30043.89900 -14892.50688 217 -8000.08089 -30043.89900 218 -6765.81908 -8000.08089 219 14287.10699 -6765.81908 220 -17245.01097 14287.10699 221 -45089.25971 -17245.01097 222 -41135.41710 -45089.25971 223 -4482.78207 -41135.41710 224 972.62061 -4482.78207 225 1768.17170 972.62061 226 11969.21276 1768.17170 227 -17124.96723 11969.21276 228 -6491.74951 -17124.96723 229 -22842.94523 -6491.74951 230 -10888.11295 -22842.94523 231 -5857.53479 -10888.11295 232 -47390.29679 -5857.53479 233 -6653.12974 -47390.29679 234 -14503.60331 -6653.12974 235 469.87415 -14503.60331 236 48761.84692 469.87415 237 -12660.04755 48761.84692 238 -12433.07969 -12660.04755 239 -13021.36899 -12433.07969 240 4841.84047 -13021.36899 241 3163.70469 4841.84047 242 -7617.15344 3163.70469 243 -2364.36904 -7617.15344 244 -7755.60678 -2364.36904 245 -10027.14952 -7755.60678 246 -16994.86809 -10027.14952 247 6612.27787 -16994.86809 248 11702.37811 6612.27787 249 -9606.36323 11702.37811 250 2770.59968 -9606.36323 251 1490.71903 2770.59968 252 5553.50598 1490.71903 253 4158.74174 5553.50598 254 -16138.56857 4158.74174 255 4347.51710 -16138.56857 256 -8687.18215 4347.51710 257 19455.14547 -8687.18215 258 9597.37430 19455.14547 259 735.18240 9597.37430 260 -8705.10687 735.18240 261 43.32399 -8705.10687 262 14247.43703 43.32399 263 157.94071 14247.43703 264 -6750.79997 157.94071 265 2456.46238 -6750.79997 266 -14221.28108 2456.46238 267 -34.11583 -14221.28108 268 32111.67627 -34.11583 269 -8479.04069 32111.67627 270 -1809.43715 -8479.04069 271 1578.99116 -1809.43715 272 5707.70370 1578.99116 273 3756.30154 5707.70370 274 8801.71997 3756.30154 275 -2034.53180 8801.71997 276 18785.13605 -2034.53180 277 3488.07624 18785.13605 278 -9155.92588 3488.07624 279 11291.44095 -9155.92588 280 1400.23558 11291.44095 281 29613.59491 1400.23558 282 9923.49655 29613.59491 283 10513.78497 9923.49655 284 4022.78901 10513.78497 285 4675.30911 4022.78901 286 6938.83049 4675.30911 287 -5378.62402 6938.83049 288 -10983.73711 -5378.62402 > 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/703hq1353319398.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/8jlfh1353319398.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/9ss5u1353319398.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/10327n1353319398.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/11kien1353319398.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/12aove1353319398.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/13al901353319398.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/14sbml1353319398.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/15iv2f1353319398.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/16fnfe1353319398.tab") + } > > try(system("convert tmp/1njjr1353319398.ps tmp/1njjr1353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/2ng2p1353319398.ps tmp/2ng2p1353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/3l9311353319398.ps tmp/3l9311353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/4dxte1353319398.ps tmp/4dxte1353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/5fyv61353319398.ps tmp/5fyv61353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/6rqgk1353319398.ps tmp/6rqgk1353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/703hq1353319398.ps tmp/703hq1353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/8jlfh1353319398.ps tmp/8jlfh1353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/9ss5u1353319398.ps tmp/9ss5u1353319398.png",intern=TRUE)) character(0) > try(system("convert tmp/10327n1353319398.ps tmp/10327n1353319398.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 22.296 1.308 23.580