R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(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 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '11' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x totsize pageviews time_in_rfc logins compendium_views_info 1 112285 1418 210907 56 396 2 84786 869 120982 56 297 3 83123 1530 176508 54 559 4 101193 2172 179321 89 967 5 38361 901 123185 40 270 6 68504 463 52746 25 143 7 119182 3201 385534 92 1562 8 22807 371 33170 18 109 9 17140 1192 101645 63 371 10 116174 1583 149061 44 656 11 57635 1439 165446 33 511 12 66198 1764 237213 84 655 13 71701 1495 173326 88 465 14 57793 1373 133131 55 525 15 80444 2187 258873 60 885 16 53855 1491 180083 66 497 17 97668 4041 324799 154 1436 18 133824 1706 230964 53 612 19 101481 2152 236785 119 865 20 99645 1036 135473 41 385 21 114789 1882 202925 61 567 22 99052 1929 215147 58 639 23 67654 2242 344297 75 963 24 65553 1220 153935 33 398 25 97500 1289 132943 40 410 26 69112 2515 174724 92 966 27 82753 2147 174415 100 801 28 85323 2352 225548 112 892 29 72654 1638 223632 73 513 30 30727 1222 124817 40 469 31 77873 1812 221698 45 683 32 117478 1677 210767 60 643 33 74007 1579 170266 62 535 34 90183 1731 260561 75 625 35 61542 807 84853 31 264 36 101494 2452 294424 77 992 37 27570 829 101011 34 238 38 55813 1940 215641 46 818 39 79215 2662 325107 99 937 40 1423 186 7176 17 70 41 55461 1499 167542 66 507 42 31081 865 106408 30 260 43 22996 1793 96560 76 503 44 83122 2527 265769 146 927 45 70106 2747 269651 67 1269 46 60578 1324 149112 56 537 47 39992 2702 175824 107 910 48 79892 1383 152871 58 532 49 49810 1179 111665 34 345 50 71570 2099 116408 61 918 51 100708 4308 362301 119 1635 52 33032 918 78800 42 330 53 82875 1831 183167 66 557 54 139077 3373 277965 89 1178 55 71595 1713 150629 44 740 56 72260 1438 168809 66 452 57 5950 496 24188 24 218 58 115762 2253 329267 259 764 59 32551 744 65029 17 255 60 31701 1161 101097 64 454 61 80670 2352 218946 41 866 62 143558 2144 244052 68 574 63 117105 4691 341570 168 1276 64 23789 1112 103597 43 379 65 120733 2694 233328 132 825 66 105195 1973 256462 105 798 67 73107 1769 206161 71 663 68 132068 3148 311473 112 1069 69 149193 2474 235800 94 921 70 46821 2084 177939 82 858 71 87011 1954 207176 70 711 72 95260 1226 196553 57 503 73 55183 1389 174184 53 382 74 106671 1496 143246 103 464 75 73511 2269 187559 121 717 76 92945 1833 187681 62 690 77 78664 1268 119016 52 462 78 70054 1943 182192 52 657 79 22618 893 73566 32 385 80 74011 1762 194979 62 577 81 83737 1403 167488 45 619 82 69094 1425 143756 46 479 83 93133 1857 275541 63 817 84 95536 1840 243199 75 752 85 225920 1502 182999 88 430 86 62133 1441 135649 46 451 87 61370 1420 152299 53 537 88 43836 1416 120221 37 519 89 106117 2970 346485 90 1000 90 38692 1317 145790 63 637 91 84651 1644 193339 78 465 92 56622 870 80953 25 437 93 15986 1654 122774 45 711 94 95364 1054 130585 46 299 95 26706 937 112611 41 248 96 89691 3004 286468 144 1162 97 67267 2008 241066 82 714 98 126846 2547 148446 91 905 99 41140 1885 204713 71 649 100 102860 1626 182079 63 512 101 51715 1468 140344 53 472 102 55801 2445 220516 62 905 103 111813 1964 243060 63 786 104 120293 1381 162765 32 489 105 138599 1369 182613 39 479 106 161647 1659 232138 62 617 107 115929 2888 265318 117 925 108 24266 1290 85574 34 351 109 162901 2845 310839 92 1144 110 109825 1982 225060 93 669 111 129838 1904 232317 54 707 112 37510 1391 144966 144 458 113 43750 602 43287 14 214 114 40652 1743 155754 61 599 115 87771 1559 164709 109 572 116 85872 2014 201940 38 897 117 89275 2143 235454 73 819 118 44418 2146 220801 75 720 119 192565 874 99466 50 273 120 35232 1590 92661 61 508 121 40909 1590 133328 55 506 122 13294 1210 61361 77 451 123 32387 2072 125930 75 699 124 140867 1281 100750 72 407 125 120662 1401 224549 50 465 126 21233 834 82316 32 245 127 44332 1105 102010 53 370 128 61056 1272 101523 42 316 129 101338 1944 243511 71 603 130 1168 391 22938 10 154 131 13497 761 41566 35 229 132 65567 1605 152474 65 577 133 25162 530 61857 25 192 134 32334 1988 99923 66 617 135 40735 1386 132487 41 411 136 91413 2395 317394 86 975 137 855 387 21054 16 146 138 97068 1742 209641 42 705 139 44339 620 22648 19 184 140 14116 449 31414 19 200 141 10288 800 46698 45 274 142 65622 1684 131698 65 502 143 16563 1050 91735 35 382 144 76643 2699 244749 95 964 145 110681 1606 184510 49 537 146 29011 1502 79863 37 438 147 92696 1204 128423 64 369 148 94785 1138 97839 38 417 149 8773 568 38214 34 276 150 83209 1459 151101 32 514 151 93815 2158 272458 65 822 152 86687 1111 172494 52 389 153 34553 1421 108043 62 466 154 105547 2833 328107 65 1255 155 103487 1955 250579 83 694 156 213688 2922 351067 95 1024 157 71220 1002 158015 29 400 158 23517 1060 98866 18 397 159 56926 956 85439 33 350 160 91721 2186 229242 247 719 161 115168 3604 351619 139 1277 162 111194 1035 84207 29 356 163 51009 1417 120445 118 457 164 135777 3261 324598 110 1402 165 51513 1587 131069 67 600 166 74163 1424 204271 42 480 167 51633 1701 165543 65 595 168 75345 1249 141722 94 436 169 33416 946 116048 64 230 170 83305 1926 250047 81 651 171 98952 3352 299775 95 1367 172 102372 1641 195838 67 564 173 37238 2035 173260 63 716 174 103772 2312 254488 83 747 175 123969 1369 104389 45 467 176 27142 1577 136084 30 671 177 135400 2201 199476 70 861 178 21399 961 92499 32 319 179 130115 1900 224330 83 612 180 24874 1254 135781 31 433 181 34988 1335 74408 67 434 182 45549 1597 81240 66 503 183 6023 207 14688 10 85 184 64466 1645 181633 70 564 185 54990 2429 271856 103 824 186 1644 151 7199 5 74 187 6179 474 46660 20 259 188 3926 141 17547 5 69 189 32755 1639 133368 36 535 190 34777 872 95227 34 239 191 73224 1318 152601 48 438 192 27114 1018 98146 40 459 193 20760 1383 79619 43 426 194 37636 1314 59194 31 288 195 65461 1335 139942 42 498 196 30080 1403 118612 46 454 197 24094 910 72880 33 376 198 69008 616 65475 18 225 199 54968 1407 99643 55 555 200 46090 771 71965 35 252 201 27507 766 77272 59 208 202 10672 473 49289 19 130 203 34029 1376 135131 66 481 204 46300 1232 108446 60 389 205 24760 1521 89746 36 565 206 18779 572 44296 25 173 207 21280 1059 77648 47 278 208 40662 1544 181528 54 609 209 28987 1230 134019 53 422 210 22827 1206 124064 40 445 211 18513 1205 92630 40 387 212 30594 1255 121848 39 339 213 24006 613 52915 14 181 214 27913 721 81872 45 245 215 42744 1109 58981 36 384 216 12934 740 53515 28 212 217 22574 1126 60812 44 399 218 41385 728 56375 30 229 219 18653 689 65490 22 224 220 18472 592 80949 17 203 221 30976 995 76302 31 333 222 63339 1613 104011 55 384 223 25568 2048 98104 54 636 224 33747 705 67989 21 185 225 4154 301 30989 14 93 226 19474 1803 135458 81 581 227 35130 799 73504 35 248 228 39067 861 63123 43 304 229 13310 1186 61254 46 344 230 65892 1451 74914 30 407 231 4143 628 31774 23 170 232 28579 1161 81437 38 312 233 51776 1463 87186 54 507 234 21152 742 50090 20 224 235 38084 979 65745 53 340 236 27717 675 56653 45 168 237 32928 1241 158399 39 443 238 11342 676 46455 20 204 239 19499 1049 73624 24 367 240 16380 620 38395 31 210 241 36874 1081 91899 35 335 242 48259 1688 139526 151 364 243 16734 736 52164 52 178 244 28207 617 51567 30 206 245 30143 812 70551 31 279 246 41369 1051 84856 29 387 247 45833 1656 102538 57 490 248 29156 705 86678 40 238 249 35944 945 85709 44 343 250 36278 554 34662 25 232 251 45588 1597 150580 77 530 252 45097 982 99611 35 291 253 3895 222 19349 11 67 254 28394 1212 99373 63 397 255 18632 1143 86230 44 467 256 2325 435 30837 19 178 257 25139 532 31706 13 175 258 27975 882 89806 42 299 259 14483 608 62088 38 154 260 13127 459 40151 29 106 261 5839 578 27634 20 189 262 24069 826 76990 27 194 263 3738 509 37460 20 135 264 18625 717 54157 19 201 265 36341 637 49862 37 207 266 24548 857 84337 26 280 267 21792 830 64175 42 260 268 26263 652 59382 49 227 269 23686 707 119308 30 239 270 49303 954 76702 49 333 271 25659 1461 103425 67 428 272 28904 672 70344 28 230 273 2781 778 43410 19 292 274 29236 1141 104838 49 350 275 19546 680 62215 27 186 276 22818 1090 69304 30 326 277 32689 616 53117 22 155 278 5752 285 19764 12 75 279 22197 1145 86680 31 361 280 20055 733 84105 20 261 281 25272 888 77945 20 299 282 82206 849 89113 39 300 283 32073 1182 91005 29 450 284 5444 528 40248 16 183 285 20154 642 64187 27 238 286 36944 947 50857 21 165 287 8019 819 56613 19 234 288 30884 757 62792 35 176 289 19540 894 72535 14 329 compendium_views_pr shared_compendiums blogged_computations 1 81 3 79 2 55 4 58 3 50 12 60 4 125 2 108 5 40 1 49 6 37 3 0 7 63 0 121 8 44 0 1 9 88 0 20 10 66 5 43 11 57 0 69 12 74 0 78 13 49 7 86 14 52 7 44 15 88 3 104 16 36 9 63 17 108 0 158 18 43 4 102 19 75 3 77 20 32 0 82 21 44 7 115 22 85 0 101 23 86 1 80 24 56 5 50 25 50 7 83 26 135 0 123 27 63 0 73 28 81 5 81 29 52 0 105 30 44 0 47 31 113 0 105 32 39 3 94 33 73 4 44 34 48 1 114 35 33 4 38 36 59 2 107 37 41 0 30 38 69 0 71 39 64 0 84 40 1 0 0 41 59 2 59 42 32 1 33 43 129 0 42 44 37 2 96 45 31 10 106 46 65 6 56 47 107 0 57 48 74 5 59 49 54 4 39 50 76 1 34 51 715 2 76 52 57 2 20 53 66 0 91 54 106 8 115 55 54 3 85 56 32 0 76 57 20 0 8 58 71 8 79 59 21 5 21 60 70 3 30 61 112 1 76 62 66 5 101 63 190 1 94 64 66 1 27 65 165 5 92 66 56 0 123 67 61 12 75 68 53 8 128 69 127 8 105 70 63 8 55 71 38 8 56 72 50 2 41 73 52 0 72 74 42 5 67 75 76 8 75 76 67 2 114 77 50 5 118 78 53 12 77 79 39 6 22 80 50 7 66 81 77 2 69 82 57 0 105 83 73 4 116 84 34 3 88 85 39 6 73 86 46 2 99 87 63 0 62 88 35 1 53 89 106 0 118 90 43 5 30 91 47 2 100 92 31 0 49 93 162 0 24 94 57 5 67 95 36 0 46 96 263 1 57 97 78 0 75 98 63 1 135 99 54 1 68 100 63 2 124 101 77 6 33 102 79 1 98 103 110 4 58 104 56 2 68 105 56 3 81 106 43 0 131 107 111 10 110 108 71 0 37 109 62 9 130 110 56 7 93 111 74 0 118 112 60 0 39 113 43 4 13 114 68 4 74 115 53 0 81 116 87 0 109 117 46 0 151 118 105 1 51 119 32 0 28 120 133 1 40 121 79 0 56 122 51 0 27 123 207 4 37 124 67 0 83 125 47 4 54 126 34 4 27 127 66 3 28 128 76 0 59 129 65 0 133 130 9 0 12 131 42 5 0 132 45 0 106 133 25 4 23 134 115 0 44 135 97 0 71 136 53 1 116 137 2 0 4 138 52 5 62 139 44 0 12 140 22 0 18 141 35 0 14 142 74 0 60 143 103 0 7 144 144 2 98 145 60 7 64 146 134 1 29 147 89 8 32 148 42 2 25 149 52 0 16 150 98 2 48 151 99 0 100 152 52 0 46 153 29 1 45 154 125 3 129 155 106 0 130 156 95 3 136 157 40 0 59 158 140 0 25 159 43 0 32 160 128 4 63 161 142 4 95 162 73 11 14 163 72 0 36 164 128 0 113 165 61 4 47 166 73 0 92 167 148 1 70 168 64 0 19 169 45 0 50 170 58 0 41 171 97 9 91 172 50 1 111 173 37 3 41 174 50 10 120 175 105 5 135 176 69 0 27 177 46 2 87 178 57 0 25 179 52 1 131 180 98 2 45 181 61 4 29 182 89 0 58 183 0 0 4 184 48 2 47 185 91 1 109 186 0 0 7 187 7 0 12 188 3 0 0 189 54 1 37 190 70 0 37 191 36 2 46 192 37 0 15 193 123 3 42 194 247 6 7 195 46 0 54 196 72 2 54 197 41 0 14 198 24 2 16 199 45 1 33 200 33 1 32 201 27 2 21 202 36 1 15 203 87 0 38 204 90 1 22 205 114 3 28 206 31 0 10 207 45 0 31 208 69 0 32 209 51 0 32 210 34 1 43 211 60 4 27 212 45 0 37 213 54 0 20 214 25 0 32 215 38 7 0 216 52 2 5 217 67 0 26 218 74 7 10 219 38 3 27 220 30 0 11 221 26 0 29 222 67 6 25 223 132 2 55 224 42 0 23 225 35 0 5 226 118 3 43 227 68 0 23 228 43 1 34 229 76 1 36 230 64 0 35 231 48 1 0 232 64 0 37 233 56 0 28 234 71 0 16 235 75 0 26 236 39 0 38 237 42 0 23 238 39 0 22 239 93 0 30 240 38 0 16 241 60 0 18 242 71 0 28 243 52 0 32 244 27 2 21 245 59 0 23 246 40 1 29 247 79 1 50 248 44 0 12 249 65 0 21 250 10 0 18 251 124 0 27 252 81 0 41 253 15 0 13 254 92 1 12 255 42 0 21 256 10 0 8 257 24 0 26 258 64 0 27 259 45 1 13 260 22 0 16 261 56 0 2 262 94 0 42 263 19 0 5 264 35 0 37 265 32 0 17 266 35 0 38 267 48 0 37 268 49 0 29 269 48 0 32 270 62 0 35 271 96 1 17 272 45 0 20 273 63 0 7 274 71 1 46 275 26 0 24 276 48 6 40 277 29 3 3 278 19 1 10 279 45 2 37 280 45 0 17 281 67 0 28 282 30 0 19 283 36 3 29 284 34 1 8 285 36 0 10 286 34 0 15 287 37 1 15 288 46 0 28 289 44 0 17 compendiums_reviewed feedback_messages_p1 feedback_messages_p120 1 30 115 94 2 28 109 103 3 38 146 93 4 30 116 103 5 22 68 51 6 26 101 70 7 25 96 91 8 18 67 22 9 11 44 38 10 26 100 93 11 25 93 60 12 38 140 123 13 44 166 148 14 30 99 90 15 40 139 124 16 34 130 70 17 47 181 168 18 30 116 115 19 31 116 71 20 23 88 66 21 36 139 134 22 36 135 117 23 30 108 108 24 25 89 84 25 39 156 156 26 34 129 120 27 31 118 114 28 31 118 94 29 33 125 120 30 25 95 81 31 33 126 110 32 35 135 133 33 42 154 122 34 43 165 158 35 30 113 109 36 33 127 124 37 13 52 39 38 32 121 92 39 36 136 126 40 0 0 0 41 28 108 70 42 14 46 37 43 17 54 38 44 32 124 120 45 30 115 93 46 35 128 95 47 20 80 77 48 28 97 90 49 28 104 80 50 39 59 31 51 34 125 110 52 26 82 66 53 39 149 138 54 39 149 133 55 33 122 113 56 28 118 100 57 4 12 7 58 39 144 140 59 18 67 61 60 14 52 41 61 29 108 96 62 44 166 164 63 21 80 78 64 16 60 49 65 28 107 102 66 35 127 124 67 28 107 99 68 38 146 129 69 23 84 62 70 36 141 73 71 32 123 114 72 29 111 99 73 25 98 70 74 27 105 104 75 36 135 116 76 28 107 91 77 23 85 74 78 40 155 138 79 23 88 67 80 40 155 151 81 28 104 72 82 34 132 120 83 33 127 115 84 28 108 105 85 34 129 104 86 30 116 108 87 33 122 98 88 22 85 69 89 38 147 111 90 26 99 99 91 35 87 71 92 8 28 27 93 24 90 69 94 29 109 107 95 20 78 73 96 29 111 107 97 45 158 93 98 37 141 129 99 33 122 69 100 33 124 118 101 25 93 73 102 32 124 119 103 29 112 104 104 28 108 107 105 28 99 99 106 31 117 90 107 52 199 197 108 21 78 36 109 24 91 85 110 41 158 139 111 33 126 106 112 32 122 50 113 19 71 64 114 20 75 31 115 31 115 63 116 31 119 92 117 32 124 106 118 18 72 63 119 23 91 69 120 17 45 41 121 20 78 56 122 12 39 25 123 17 68 65 124 30 119 93 125 31 117 114 126 10 39 38 127 13 50 44 128 22 88 87 129 42 155 110 130 1 0 0 131 9 36 27 132 32 123 83 133 11 32 30 134 25 99 80 135 36 136 98 136 31 117 82 137 0 0 0 138 24 88 60 139 13 39 28 140 8 25 9 141 13 52 33 142 19 75 59 143 18 71 49 144 33 124 115 145 40 151 140 146 22 71 49 147 38 145 120 148 24 87 66 149 8 27 21 150 35 131 124 151 43 162 152 152 43 165 139 153 14 54 38 154 41 159 144 155 38 147 120 156 45 170 160 157 31 119 114 158 13 49 39 159 28 104 78 160 31 120 119 161 40 150 141 162 30 112 101 163 16 59 56 164 37 136 133 165 30 107 83 166 35 130 116 167 32 115 90 168 27 107 36 169 20 75 50 170 18 71 61 171 31 120 97 172 31 116 98 173 21 79 78 174 39 150 117 175 41 156 148 176 13 51 41 177 32 118 105 178 18 71 55 179 39 144 132 180 14 47 44 181 7 28 21 182 17 68 50 183 0 0 0 184 30 110 73 185 37 147 86 186 0 0 0 187 5 15 13 188 1 4 4 189 16 64 57 190 32 111 48 191 24 85 46 192 17 68 48 193 11 40 32 194 24 80 68 195 22 88 87 196 12 48 43 197 19 76 67 198 13 51 46 199 17 67 46 200 15 59 56 201 16 61 48 202 24 76 44 203 15 60 60 204 17 68 65 205 18 71 55 206 20 76 38 207 16 62 52 208 16 61 60 209 18 67 54 210 22 88 86 211 8 30 24 212 17 64 52 213 18 68 49 214 16 64 61 215 23 91 61 216 22 88 81 217 13 52 43 218 13 49 40 219 16 62 40 220 16 61 56 221 20 76 68 222 22 88 79 223 17 66 47 224 18 71 57 225 17 68 41 226 12 48 29 227 7 25 3 228 17 68 60 229 14 41 30 230 23 90 79 231 17 66 47 232 14 54 40 233 15 59 48 234 17 60 36 235 21 77 42 236 18 68 49 237 18 72 57 238 17 67 12 239 17 64 40 240 16 63 43 241 15 59 33 242 21 84 77 243 16 64 43 244 14 56 45 245 15 54 47 246 17 67 43 247 15 58 45 248 15 59 50 249 10 40 35 250 6 22 7 251 22 83 71 252 21 81 67 253 1 2 0 254 18 72 62 255 17 61 54 256 4 15 4 257 10 32 25 258 16 62 40 259 16 58 38 260 9 36 19 261 16 59 17 262 17 68 67 263 7 21 14 264 15 55 30 265 14 54 54 266 14 55 35 267 18 72 59 268 12 41 24 269 16 61 58 270 21 67 42 271 19 76 46 272 16 64 61 273 1 3 3 274 16 63 52 275 10 40 25 276 19 69 40 277 12 48 32 278 2 8 4 279 14 52 49 280 17 66 63 281 19 76 67 282 14 43 32 283 11 39 23 284 4 14 7 285 16 61 54 286 20 71 37 287 12 44 35 288 15 60 51 289 16 64 39 totrevisions totseconds tothyperlinks totblogs 1 24188 146283 144 145 2 18273 98364 103 101 3 14130 86146 98 98 4 32287 96933 135 132 5 8654 79234 61 60 6 9245 42551 39 38 7 33251 195663 150 144 8 1271 6853 5 5 9 5279 21529 28 28 10 27101 95757 84 84 11 16373 85584 80 79 12 19716 143983 130 127 13 17753 75851 82 78 14 9028 59238 60 60 15 18653 93163 131 131 16 8828 96037 84 84 17 29498 151511 140 133 18 27563 136368 151 150 19 18293 112642 91 91 20 22530 94728 138 132 21 15977 105499 150 136 22 35082 121527 124 124 23 16116 127766 119 118 24 15849 98958 73 70 25 16026 77900 110 107 26 26569 85646 123 119 27 24785 98579 90 89 28 17569 130767 116 112 29 23825 131741 113 108 30 7869 53907 56 52 31 14975 178812 115 112 32 37791 146761 119 116 33 9605 82036 129 123 34 27295 163253 127 125 35 2746 27032 27 27 36 34461 171975 175 162 37 8098 65990 35 32 38 4787 86572 64 64 39 24919 159676 96 92 40 603 1929 0 0 41 16329 85371 84 83 42 12558 58391 41 41 43 7784 31580 47 47 44 28522 136815 126 120 45 22265 120642 105 105 46 14459 69107 80 79 47 14526 50495 70 65 48 22240 108016 73 70 49 11802 46341 57 55 50 7623 78348 40 39 51 11912 79336 68 67 52 7935 56968 21 21 53 18220 93176 127 127 54 19199 161632 154 152 55 19918 87850 116 113 56 21884 127969 102 99 57 2694 15049 7 7 58 15808 155135 148 141 59 3597 25109 21 21 60 5296 45824 35 35 61 25239 102996 112 109 62 29801 160604 137 133 63 18450 158051 135 123 64 7132 44547 26 26 65 34861 162647 230 230 66 35940 174141 181 166 67 16688 60622 71 68 68 24683 179566 147 147 69 46230 184301 190 179 70 10387 75661 64 61 71 21436 96144 105 101 72 30546 129847 107 108 73 19746 117286 94 90 74 15977 71180 116 114 75 22583 109377 106 103 76 17274 85298 143 142 77 16469 73631 81 79 78 14251 86767 89 88 79 3007 23824 26 25 80 16851 93487 84 83 81 21113 82981 113 113 82 17401 73815 120 118 83 23958 94552 110 110 84 23567 132190 134 129 85 13065 128754 54 51 86 15358 66363 96 93 87 14587 67808 78 76 88 12770 61724 51 49 89 24021 131722 121 118 90 9648 68580 38 38 91 20537 106175 145 141 92 7905 55792 59 58 93 4527 25157 27 27 94 30495 76669 91 91 95 7117 57283 48 48 96 17719 105805 68 63 97 27056 129484 58 56 98 33473 72413 150 144 99 9758 87831 74 73 100 21115 96971 181 168 101 7236 71299 65 64 102 13790 77494 97 97 103 32902 120336 121 117 104 25131 93913 99 100 105 30910 136048 152 149 106 35947 181248 188 187 107 29848 146123 138 127 108 6943 32036 40 37 109 42705 186646 254 245 110 31808 102255 87 87 111 26675 168237 178 177 112 8435 64219 51 49 113 7409 19630 49 49 114 14993 76825 73 73 115 36867 115338 176 177 116 33835 109427 94 94 117 24164 118168 120 117 118 12607 84845 66 60 119 22609 153197 56 55 120 5892 29877 39 39 121 17014 63506 66 64 122 5394 22445 27 26 123 9178 47695 65 64 124 6440 68370 58 58 125 21916 146304 98 95 126 4011 38233 25 25 127 5818 42071 26 26 128 18647 50517 77 76 129 20556 103950 130 129 130 238 5841 11 11 131 70 2341 2 2 132 22392 84396 101 101 133 3913 24610 31 28 134 12237 35753 36 36 135 8388 55515 120 89 136 22120 209056 195 193 137 338 6622 4 4 138 11727 115814 89 84 139 3704 11609 24 23 140 3988 13155 39 39 141 3030 18274 14 14 142 13520 72875 78 78 143 1421 10112 15 14 144 20923 142775 106 101 145 20237 68847 83 82 146 3219 17659 24 24 147 3769 20112 37 36 148 12252 61023 77 75 149 1888 13983 16 16 150 14497 65176 56 55 151 28864 132432 132 131 152 21721 112494 144 131 153 4821 45109 40 39 154 33644 170875 153 144 155 15923 180759 143 139 156 42935 214921 220 211 157 18864 100226 79 78 158 4977 32043 50 50 159 7785 54454 39 39 160 17939 78876 95 90 161 23436 170745 169 166 162 325 6940 12 12 163 13539 49025 63 57 164 34538 122037 134 133 165 12198 53782 69 69 166 26924 127748 119 119 167 12716 86839 119 119 168 8172 44830 75 65 169 10855 77395 63 61 170 11932 89324 55 49 171 14300 103300 103 101 172 25515 112283 197 196 173 2805 10901 16 15 174 29402 120691 140 136 175 16440 58106 89 89 176 11221 57140 40 40 177 28732 122422 125 123 178 5250 25899 21 21 179 28608 139296 167 163 180 8092 52678 32 29 181 4473 23853 36 35 182 1572 17306 13 13 183 2065 7953 5 5 184 14817 89455 96 96 185 16714 147866 151 151 186 556 4245 6 6 187 2089 21509 13 13 188 2658 7670 3 3 189 10695 66675 57 56 190 1669 14336 23 23 191 16267 53608 61 57 192 7768 30059 21 14 193 7252 29668 43 43 194 6387 22097 20 20 195 18715 96841 82 72 196 7936 41907 90 87 197 8643 27080 25 21 198 7294 35885 60 56 199 4570 41247 61 59 200 7185 28313 85 82 201 10058 36845 43 43 202 2342 16548 25 25 203 8509 36134 41 38 204 13275 55764 26 25 205 6816 28910 38 38 206 1930 13339 12 12 207 8086 25319 29 29 208 10737 66956 49 47 209 8033 47487 46 45 210 7058 52785 41 40 211 6782 44683 31 30 212 5401 35619 41 41 213 6521 21920 26 25 214 10856 45608 23 23 215 2154 7721 14 14 216 6117 20634 16 16 217 5238 29788 25 26 218 4820 31931 21 21 219 5615 37754 32 27 220 4272 32505 9 9 221 8702 40557 35 33 222 15340 94238 42 42 223 8030 44197 68 68 224 9526 43228 32 32 225 1278 4103 6 6 226 4236 44144 68 67 227 3023 32868 33 33 228 7196 27640 84 77 229 3394 14063 46 46 230 6371 28990 30 30 231 1574 4694 0 0 232 9620 42648 36 36 233 6978 64329 47 46 234 4911 21928 20 18 235 8645 25836 50 48 236 8987 22779 30 29 237 5544 40820 30 28 238 3083 27530 34 34 239 6909 32378 33 33 240 3189 10824 34 34 241 6745 39613 37 33 242 16724 60865 83 80 243 4850 19787 32 32 244 7025 20107 30 30 245 6047 36605 43 41 246 7377 40961 41 41 247 9078 48231 51 51 248 4605 39725 19 18 249 3238 21455 37 34 250 8100 23430 33 31 251 9653 62991 41 39 252 8914 49363 54 54 253 786 9604 14 14 254 6700 24552 25 24 255 5788 31493 25 24 256 593 3439 8 8 257 4506 19555 26 26 258 6382 21228 20 19 259 5621 23177 11 11 260 3997 22094 14 14 261 520 2342 3 1 262 8891 38798 40 39 263 999 3255 5 5 264 7067 24261 38 37 265 4639 18511 32 32 266 5654 40798 41 38 267 6928 28893 46 47 268 1514 21425 47 47 269 9238 50276 37 37 270 8204 37643 51 51 271 5926 30377 49 45 272 5785 27126 21 21 273 4 13 1 1 274 5930 42097 44 42 275 3710 24451 26 26 276 705 14335 21 21 277 443 5084 4 4 278 2416 9927 10 10 279 7747 43527 43 43 280 5432 27184 34 34 281 4913 21610 32 31 282 2650 20484 20 19 283 2370 20156 34 34 284 775 6012 6 6 285 5576 18475 12 11 286 1352 12645 24 24 287 3080 11017 16 16 288 10205 37623 72 72 289 6095 35873 27 21 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) pageviews time_in_rfc 1046.85044 0.72345 -0.08485 logins compendium_views_info compendium_views_pr 39.75784 -7.71169 40.06514 shared_compendiums blogged_computations compendiums_reviewed 2166.01565 47.61852 886.23668 feedback_messages_p1 feedback_messages_p120 totrevisions -198.87838 187.16938 1.00593 totseconds tothyperlinks totblogs 0.36978 18.18514 55.13571 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -35704 -11474 -3683 8165 132608 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1046.85044 3275.40858 0.320 0.749508 pageviews 0.72345 7.52564 0.096 0.923486 time_in_rfc -0.08485 0.05598 -1.516 0.130720 logins 39.75784 57.40904 0.693 0.489187 compendium_views_info -7.71169 17.26933 -0.447 0.655550 compendium_views_pr 40.06514 30.70711 1.305 0.193072 shared_compendiums 2166.01565 493.97916 4.385 1.66e-05 *** blogged_computations 47.61852 85.17113 0.559 0.576555 compendiums_reviewed 886.23668 641.42927 1.382 0.168202 feedback_messages_p1 -198.87838 194.96261 -1.020 0.308588 feedback_messages_p120 187.16938 96.56789 1.938 0.053624 . totrevisions 1.00593 0.29369 3.425 0.000709 *** totseconds 0.36978 0.07852 4.709 3.95e-06 *** tothyperlinks 18.18514 406.38324 0.045 0.964340 totblogs 55.13571 421.16916 0.131 0.895942 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 20040 on 274 degrees of freedom Multiple R-squared: 0.7663, Adjusted R-squared: 0.7543 F-statistic: 64.16 on 14 and 274 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.2183073928 4.366148e-01 7.816926e-01 [2,] 0.1644663165 3.289326e-01 8.355337e-01 [3,] 0.0826251108 1.652502e-01 9.173749e-01 [4,] 0.4170847189 8.341694e-01 5.829153e-01 [5,] 0.3775630532 7.551261e-01 6.224369e-01 [6,] 0.2762473631 5.524947e-01 7.237526e-01 [7,] 0.1944009260 3.888019e-01 8.055991e-01 [8,] 0.1500114855 3.000230e-01 8.499885e-01 [9,] 0.1018819174 2.037638e-01 8.981181e-01 [10,] 0.0885333567 1.770667e-01 9.114666e-01 [11,] 0.0574494270 1.148989e-01 9.425506e-01 [12,] 0.0396550182 7.931004e-02 9.603450e-01 [13,] 0.0304402287 6.088046e-02 9.695598e-01 [14,] 0.0510687256 1.021375e-01 9.489313e-01 [15,] 0.0364578683 7.291574e-02 9.635421e-01 [16,] 0.0240755432 4.815109e-02 9.759245e-01 [17,] 0.0166690026 3.333801e-02 9.833310e-01 [18,] 0.0278742737 5.574855e-02 9.721257e-01 [19,] 0.0327870344 6.557407e-02 9.672130e-01 [20,] 0.0272719799 5.454396e-02 9.727280e-01 [21,] 0.0198628583 3.972572e-02 9.801371e-01 [22,] 0.0145250395 2.905008e-02 9.854750e-01 [23,] 0.0214847246 4.296945e-02 9.785153e-01 [24,] 0.0170808166 3.416163e-02 9.829192e-01 [25,] 0.0130454458 2.609089e-02 9.869546e-01 [26,] 0.0089486713 1.789734e-02 9.910513e-01 [27,] 0.0071207659 1.424153e-02 9.928792e-01 [28,] 0.0244922556 4.898451e-02 9.755077e-01 [29,] 0.0189430927 3.788619e-02 9.810569e-01 [30,] 0.0140520407 2.810408e-02 9.859480e-01 [31,] 0.0106823648 2.136473e-02 9.893176e-01 [32,] 0.0074311047 1.486221e-02 9.925689e-01 [33,] 0.0105117171 2.102343e-02 9.894883e-01 [34,] 0.0157048042 3.140961e-02 9.842952e-01 [35,] 0.0134525239 2.690505e-02 9.865475e-01 [36,] 0.0095083772 1.901675e-02 9.904916e-01 [37,] 0.0109707902 2.194158e-02 9.890292e-01 [38,] 0.0114017481 2.280350e-02 9.885983e-01 [39,] 0.0086146242 1.722925e-02 9.913854e-01 [40,] 0.0062016490 1.240330e-02 9.937984e-01 [41,] 0.0061255703 1.225114e-02 9.938744e-01 [42,] 0.0042628531 8.525706e-03 9.957371e-01 [43,] 0.0029941266 5.988253e-03 9.970059e-01 [44,] 0.0021139120 4.227824e-03 9.978861e-01 [45,] 0.0024635930 4.927186e-03 9.975364e-01 [46,] 0.0025111718 5.022344e-03 9.974888e-01 [47,] 0.0018193373 3.638675e-03 9.981807e-01 [48,] 0.0053764528 1.075291e-02 9.946235e-01 [49,] 0.0051284853 1.025697e-02 9.948715e-01 [50,] 0.0037687376 7.537475e-03 9.962313e-01 [51,] 0.0030058854 6.011771e-03 9.969941e-01 [52,] 0.0020940129 4.188026e-03 9.979060e-01 [53,] 0.0020210891 4.042178e-03 9.979789e-01 [54,] 0.0015304914 3.060983e-03 9.984695e-01 [55,] 0.0010701357 2.140271e-03 9.989299e-01 [56,] 0.0009566743 1.913349e-03 9.990433e-01 [57,] 0.0017563771 3.512754e-03 9.982436e-01 [58,] 0.0034093712 6.818742e-03 9.965906e-01 [59,] 0.0032649013 6.529803e-03 9.967351e-01 [60,] 0.0027752955 5.550591e-03 9.972247e-01 [61,] 0.0043004798 8.600960e-03 9.956995e-01 [62,] 0.0038030061 7.606012e-03 9.961970e-01 [63,] 0.0037039999 7.408000e-03 9.962960e-01 [64,] 0.0029830964 5.966193e-03 9.970169e-01 [65,] 0.0021974616 4.394923e-03 9.978025e-01 [66,] 0.0017531672 3.506334e-03 9.982468e-01 [67,] 0.0012919306 2.583861e-03 9.987081e-01 [68,] 0.9739033183 5.219336e-02 2.609668e-02 [69,] 0.9687170662 6.256587e-02 3.128293e-02 [70,] 0.9624627881 7.507442e-02 3.753721e-02 [71,] 0.9550932981 8.981340e-02 4.490670e-02 [72,] 0.9559903496 8.801930e-02 4.400965e-02 [73,] 0.9610090496 7.798190e-02 3.899095e-02 [74,] 0.9573675593 8.526488e-02 4.263244e-02 [75,] 0.9568939709 8.621206e-02 4.310603e-02 [76,] 0.9491225975 1.017548e-01 5.087740e-02 [77,] 0.9392810371 1.214379e-01 6.071896e-02 [78,] 0.9358260933 1.283478e-01 6.417391e-02 [79,] 0.9350918491 1.298163e-01 6.490815e-02 [80,] 0.9406525598 1.186949e-01 5.934744e-02 [81,] 0.9544489688 9.110206e-02 4.555103e-02 [82,] 0.9511360217 9.772796e-02 4.886398e-02 [83,] 0.9427612795 1.144774e-01 5.723872e-02 [84,] 0.9355651564 1.288697e-01 6.443484e-02 [85,] 0.9264716393 1.470567e-01 7.352836e-02 [86,] 0.9194736911 1.610526e-01 8.052631e-02 [87,] 0.9448296907 1.103406e-01 5.517031e-02 [88,] 0.9490004738 1.019991e-01 5.099953e-02 [89,] 0.9638366583 7.232668e-02 3.616334e-02 [90,] 0.9777808911 4.443822e-02 2.221911e-02 [91,] 0.9726742095 5.465158e-02 2.732579e-02 [92,] 0.9689915355 6.201693e-02 3.100846e-02 [93,] 0.9640675928 7.186481e-02 3.593241e-02 [94,] 0.9611437979 7.771240e-02 3.885620e-02 [95,] 0.9581187052 8.376259e-02 4.188129e-02 [96,] 0.9500643889 9.987122e-02 4.993561e-02 [97,] 0.9514137749 9.717245e-02 4.858623e-02 [98,] 0.9471391377 1.057217e-01 5.286086e-02 [99,] 0.9374318846 1.251362e-01 6.256812e-02 [100,] 0.9260441465 1.479117e-01 7.395585e-02 [101,] 0.9137118223 1.725764e-01 8.628818e-02 [102,] 0.9997117298 5.765404e-04 2.882702e-04 [103,] 0.9996062704 7.874593e-04 3.937296e-04 [104,] 0.9994970956 1.005809e-03 5.029044e-04 [105,] 0.9994248432 1.150314e-03 5.751568e-04 [106,] 0.9994000388 1.199922e-03 5.999612e-04 [107,] 0.9999994219 1.156290e-06 5.781451e-07 [108,] 0.9999994384 1.123264e-06 5.616318e-07 [109,] 0.9999992843 1.431326e-06 7.156629e-07 [110,] 0.9999990116 1.976811e-06 9.884053e-07 [111,] 0.9999985401 2.919832e-06 1.459916e-06 [112,] 0.9999982261 3.547804e-06 1.773902e-06 [113,] 0.9999973369 5.326269e-06 2.663134e-06 [114,] 0.9999961623 7.675407e-06 3.837703e-06 [115,] 0.9999949653 1.006939e-05 5.034697e-06 [116,] 0.9999939400 1.211992e-05 6.059960e-06 [117,] 0.9999917815 1.643697e-05 8.218485e-06 [118,] 0.9999967927 6.414591e-06 3.207295e-06 [119,] 0.9999972204 5.559129e-06 2.779565e-06 [120,] 0.9999957888 8.422430e-06 4.211215e-06 [121,] 0.9999963473 7.305304e-06 3.652652e-06 [122,] 0.9999966128 6.774430e-06 3.387215e-06 [123,] 0.9999949497 1.010054e-05 5.050272e-06 [124,] 0.9999927197 1.456052e-05 7.280259e-06 [125,] 0.9999923904 1.521922e-05 7.609612e-06 [126,] 0.9999891301 2.173986e-05 1.086993e-05 [127,] 0.9999917181 1.656370e-05 8.281852e-06 [128,] 0.9999920240 1.595209e-05 7.976046e-06 [129,] 0.9999884247 2.315061e-05 1.157530e-05 [130,] 0.9999950206 9.958715e-06 4.979358e-06 [131,] 0.9999985608 2.878460e-06 1.439230e-06 [132,] 0.9999978642 4.271511e-06 2.135755e-06 [133,] 0.9999974312 5.137522e-06 2.568761e-06 [134,] 0.9999967315 6.536910e-06 3.268455e-06 [135,] 0.9999971447 5.710600e-06 2.855300e-06 [136,] 0.9999956771 8.645839e-06 4.322919e-06 [137,] 0.9999981943 3.611488e-06 1.805744e-06 [138,] 0.9999977055 4.589083e-06 2.294541e-06 [139,] 0.9999999102 1.796446e-07 8.982231e-08 [140,] 0.9999998581 2.837707e-07 1.418853e-07 [141,] 0.9999997938 4.123120e-07 2.061560e-07 [142,] 0.9999997217 5.566213e-07 2.783106e-07 [143,] 0.9999996165 7.670428e-07 3.835214e-07 [144,] 0.9999994072 1.185696e-06 5.928482e-07 [145,] 0.9999999963 7.393781e-09 3.696890e-09 [146,] 0.9999999948 1.041328e-08 5.206641e-09 [147,] 0.9999999967 6.623118e-09 3.311559e-09 [148,] 0.9999999966 6.716676e-09 3.358338e-09 [149,] 0.9999999964 7.100222e-09 3.550111e-09 [150,] 0.9999999963 7.455314e-09 3.727657e-09 [151,] 0.9999999993 1.369205e-09 6.846027e-10 [152,] 0.9999999992 1.636854e-09 8.184271e-10 [153,] 0.9999999999 1.477511e-10 7.387553e-11 [154,] 0.9999999999 2.190841e-10 1.095421e-10 [155,] 0.9999999998 3.990233e-10 1.995116e-10 [156,] 0.9999999997 5.919290e-10 2.959645e-10 [157,] 0.9999999998 4.858556e-10 2.429278e-10 [158,] 0.9999999999 1.211039e-10 6.055195e-11 [159,] 0.9999999999 1.748533e-10 8.742667e-11 [160,] 1.0000000000 4.639241e-11 2.319621e-11 [161,] 1.0000000000 8.895010e-11 4.447505e-11 [162,] 1.0000000000 3.249298e-11 1.624649e-11 [163,] 1.0000000000 4.559026e-11 2.279513e-11 [164,] 1.0000000000 8.240439e-11 4.120219e-11 [165,] 1.0000000000 1.696289e-11 8.481446e-12 [166,] 1.0000000000 3.358294e-11 1.679147e-11 [167,] 1.0000000000 6.322986e-11 3.161493e-11 [168,] 1.0000000000 3.753636e-11 1.876818e-11 [169,] 1.0000000000 7.398518e-11 3.699259e-11 [170,] 0.9999999999 1.059731e-10 5.298654e-11 [171,] 0.9999999999 2.009244e-10 1.004622e-10 [172,] 0.9999999999 2.085018e-10 1.042509e-10 [173,] 0.9999999998 3.260076e-10 1.630038e-10 [174,] 0.9999999999 2.420759e-10 1.210380e-10 [175,] 0.9999999998 4.614632e-10 2.307316e-10 [176,] 0.9999999996 7.868076e-10 3.934038e-10 [177,] 0.9999999995 1.047215e-09 5.236074e-10 [178,] 0.9999999990 1.995796e-09 9.978981e-10 [179,] 0.9999999985 2.945125e-09 1.472563e-09 [180,] 0.9999999977 4.631151e-09 2.315575e-09 [181,] 0.9999999996 8.527825e-10 4.263912e-10 [182,] 0.9999999995 9.033448e-10 4.516724e-10 [183,] 0.9999999994 1.247250e-09 6.236249e-10 [184,] 0.9999999990 2.030926e-09 1.015463e-09 [185,] 0.9999999996 7.908456e-10 3.954228e-10 [186,] 0.9999999993 1.420805e-09 7.104026e-10 [187,] 0.9999999990 1.939784e-09 9.698920e-10 [188,] 0.9999999982 3.626105e-09 1.813053e-09 [189,] 0.9999999965 6.983076e-09 3.491538e-09 [190,] 0.9999999938 1.245756e-08 6.228778e-09 [191,] 0.9999999886 2.289443e-08 1.144721e-08 [192,] 0.9999999848 3.048092e-08 1.524046e-08 [193,] 0.9999999895 2.092449e-08 1.046225e-08 [194,] 0.9999999876 2.482064e-08 1.241032e-08 [195,] 0.9999999782 4.367617e-08 2.183809e-08 [196,] 0.9999999586 8.280455e-08 4.140227e-08 [197,] 0.9999999261 1.478283e-07 7.391413e-08 [198,] 0.9999998982 2.036075e-07 1.018038e-07 [199,] 0.9999998872 2.255288e-07 1.127644e-07 [200,] 0.9999997879 4.242415e-07 2.121207e-07 [201,] 0.9999998530 2.940816e-07 1.470408e-07 [202,] 0.9999997749 4.502644e-07 2.251322e-07 [203,] 0.9999996917 6.165479e-07 3.082740e-07 [204,] 0.9999996531 6.937193e-07 3.468597e-07 [205,] 0.9999994862 1.027580e-06 5.137898e-07 [206,] 0.9999991944 1.611115e-06 8.055574e-07 [207,] 0.9999985552 2.889608e-06 1.444804e-06 [208,] 0.9999975161 4.967751e-06 2.483875e-06 [209,] 0.9999967364 6.527276e-06 3.263638e-06 [210,] 0.9999976454 4.709195e-06 2.354598e-06 [211,] 0.9999956539 8.692290e-06 4.346145e-06 [212,] 0.9999983534 3.293132e-06 1.646566e-06 [213,] 0.9999996679 6.641801e-07 3.320900e-07 [214,] 0.9999994190 1.161991e-06 5.809953e-07 [215,] 0.9999988606 2.278838e-06 1.139419e-06 [216,] 0.9999987827 2.434576e-06 1.217288e-06 [217,] 0.9999976380 4.724038e-06 2.362019e-06 [218,] 0.9999954763 9.047344e-06 4.523672e-06 [219,] 0.9999912950 1.741003e-05 8.705016e-06 [220,] 0.9999883162 2.336767e-05 1.168383e-05 [221,] 0.9999793595 4.128093e-05 2.064047e-05 [222,] 0.9999622328 7.553442e-05 3.776721e-05 [223,] 0.9999312984 1.374032e-04 6.870161e-05 [224,] 0.9998905940 2.188120e-04 1.094060e-04 [225,] 0.9998881088 2.237825e-04 1.118912e-04 [226,] 0.9998099798 3.800405e-04 1.900202e-04 [227,] 0.9996887524 6.224953e-04 3.112476e-04 [228,] 0.9994462111 1.107578e-03 5.537889e-04 [229,] 0.9995259149 9.481701e-04 4.740851e-04 [230,] 0.9995632912 8.734175e-04 4.367088e-04 [231,] 0.9992530710 1.493858e-03 7.469290e-04 [232,] 0.9987552961 2.489408e-03 1.244704e-03 [233,] 0.9995025956 9.948088e-04 4.974044e-04 [234,] 0.9991183487 1.763303e-03 8.816513e-04 [235,] 0.9992584733 1.483053e-03 7.415267e-04 [236,] 0.9986696774 2.660645e-03 1.330323e-03 [237,] 0.9976964309 4.607138e-03 2.303569e-03 [238,] 0.9972996938 5.400612e-03 2.700306e-03 [239,] 0.9957248043 8.550391e-03 4.275196e-03 [240,] 0.9925019877 1.499602e-02 7.498012e-03 [241,] 0.9877541496 2.449170e-02 1.224585e-02 [242,] 0.9824797102 3.504058e-02 1.752029e-02 [243,] 0.9712932025 5.741359e-02 2.870680e-02 [244,] 0.9528017815 9.439644e-02 4.719822e-02 [245,] 0.9458958041 1.082084e-01 5.410420e-02 [246,] 0.9838117255 3.237655e-02 1.618827e-02 [247,] 0.9706074719 5.878506e-02 2.939253e-02 [248,] 0.9469227554 1.061545e-01 5.307724e-02 [249,] 0.9081833001 1.836334e-01 9.181670e-02 [250,] 0.8489517667 3.020965e-01 1.510482e-01 [251,] 0.9003320311 1.993359e-01 9.966797e-02 [252,] 0.8496615965 3.006768e-01 1.503384e-01 [253,] 0.7983448653 4.033103e-01 2.016551e-01 [254,] 0.7134599860 5.730800e-01 2.865400e-01 > postscript(file="/var/www/rcomp/tmp/10nv31323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/2tamq1323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/3t7m61323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/4dqw21323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5sr731323768015.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 5083.48536 -4800.76920 -8043.89878 2375.70234 -16436.07273 19830.28558 7 8 9 10 11 12 12071.46232 11696.37084 -3700.86884 25027.12200 -2173.23653 -32261.30065 13 14 15 16 17 18 -22233.30128 -9558.52325 -4588.48851 -23649.19944 -18319.78623 24298.90482 19 20 21 22 23 24 23320.22095 22781.29675 14495.16514 -7002.41992 -8226.26911 -14912.55878 25 26 27 28 29 30 1691.68102 -20534.23561 -2117.59936 -9705.85260 -24168.82126 -12745.27790 31 32 33 34 35 36 -25949.20119 -6768.62149 -4311.96334 -29578.29108 16865.86144 -21116.14047 37 38 39 40 41 42 -11520.66697 6958.93408 -18816.00091 -645.52210 -11209.11462 -13184.37475 43 44 45 46 47 48 -11858.53404 -20592.52130 -26704.37502 -14378.89819 -6695.63307 -13861.11403 49 50 51 52 53 54 -6172.42326 8557.31534 26455.37045 -18360.29159 -2114.48728 17584.88095 55 56 57 58 59 60 -13655.42774 -16253.83483 -5115.16399 -3177.14353 -3384.65006 -5959.95249 61 62 63 64 65 66 -2435.51326 9241.28877 27013.89593 -11126.28380 -17872.68381 -23599.53510 67 68 69 70 71 72 -7752.38017 2624.43538 -947.73671 -20229.29690 -6515.43292 -5467.94012 73 74 75 76 77 78 -22434.17430 27177.26939 -34024.72121 16993.83501 2889.87194 -29065.35284 79 80 81 82 83 84 -15850.28126 -20599.75939 10353.12905 -6032.36082 8164.73331 -564.89400 85 86 87 88 89 90 132607.86690 -8279.90796 94.91248 -6401.20169 19793.66281 -22244.59881 91 92 93 94 95 96 -9007.88347 21187.35377 -11900.64469 -940.98459 -16549.97234 13285.71527 97 98 99 100 101 102 -24331.71565 28785.71284 -14847.49852 8446.32055 -10283.15755 -9398.40806 103 104 105 106 107 108 8554.95004 34059.32689 25565.84324 35469.09788 -30507.82064 -5138.79986 109 110 111 112 113 114 13983.00827 -1048.15915 15781.24918 -8012.45384 3027.04961 -18611.55913 115 116 117 118 119 120 -15419.46178 -5184.17275 861.84933 -7665.38091 98184.92438 -634.52131 121 122 123 124 125 126 -12111.92312 -9542.01725 -21072.21950 84741.84472 17426.78698 -12351.45079 127 128 129 130 131 132 8283.61894 1308.14094 16859.73758 -2449.42331 -3664.60148 -9258.87716 133 134 135 136 137 138 -5437.35369 -11851.29816 -15055.46292 -23443.56577 -1548.13237 24580.09645 139 140 141 142 143 144 24365.76431 -1079.80948 -7307.99819 10764.44052 2418.21319 -24147.68301 145 146 147 148 149 150 22239.65095 569.21328 38372.42558 38970.99939 -5404.52501 15166.67681 151 152 153 154 155 156 -12616.36309 -8467.28978 4115.81231 -19681.42939 -6121.88619 53340.09731 157 158 159 160 161 162 -6542.74309 -3958.40085 11020.67386 7693.62901 -4153.20707 63728.41988 163 164 165 166 167 168 4427.40501 35566.92673 -9002.38745 -25777.99910 -21984.03050 42444.28788 169 170 171 172 173 174 -19637.54427 37621.72139 18179.81442 4664.78334 17034.97488 -15150.37859 175 176 177 178 179 180 33575.13231 -7968.80423 35454.85787 -3304.71929 14719.39321 -14948.81996 181 182 183 184 185 186 7013.26562 25332.82318 755.45414 -2681.83425 -35703.83763 -1431.62282 187 188 189 190 191 192 -5790.59168 -2090.15439 -13215.69904 13169.75089 24263.31846 2801.79890 193 194 195 196 197 198 -15583.24768 -15724.97028 -5297.27155 -9083.02604 -6733.31607 33261.69580 199 200 201 202 203 204 23348.52976 10326.41207 -12620.01279 -15308.65792 1890.17829 -2174.88336 205 206 207 208 209 210 -11704.15866 2222.73603 -8153.06201 -157.95627 -5690.30853 -19961.74832 211 212 213 214 215 216 -17683.72874 2796.23867 -3387.29176 -11474.07479 10993.74415 -22352.22803 217 218 219 220 221 222 -4122.09586 -893.35584 -16804.38726 -6477.14820 -6947.55255 -16285.33372 223 224 225 226 227 228 -18261.59819 -4641.12939 -8419.87186 -15428.43762 17035.82560 1875.51234 229 230 231 232 233 234 -11354.07007 31268.10109 -12410.76397 -7070.00647 10513.00605 -3131.84042 235 236 237 238 239 240 5215.27167 -4084.41497 8449.33829 -7845.58923 -10505.01416 -3156.93538 241 242 243 244 245 246 9112.12501 -11111.42584 -8762.66804 -1322.70843 -2057.06271 8081.08372 247 248 249 250 251 252 5201.95444 1286.82579 17581.19517 15935.37064 -1243.23221 1225.17846 253 254 255 256 257 258 -2668.30154 -1798.23473 -9042.88249 -349.49334 4183.09042 5679.50345 259 260 261 262 263 264 -11435.07675 -4413.67007 -2084.38125 -16636.84538 -2474.68687 -7496.20363 265 266 267 268 269 270 11465.12509 -3683.28547 -11516.70696 6376.99395 -14045.98340 11484.69168 271 272 273 274 275 276 -4753.54481 986.45543 2560.43349 -6221.52354 1196.40415 -6791.87269 277 278 279 280 281 282 18524.46283 -4890.61380 -16946.62533 -7521.39460 -1798.78534 65266.10451 283 284 285 286 287 288 12261.48287 -1217.24202 -2488.29780 19563.18185 -8910.56503 -8956.28845 289 -6386.91361 > postscript(file="/var/www/rcomp/tmp/64sef1323768015.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 5083.48536 NA 1 -4800.76920 5083.48536 2 -8043.89878 -4800.76920 3 2375.70234 -8043.89878 4 -16436.07273 2375.70234 5 19830.28558 -16436.07273 6 12071.46232 19830.28558 7 11696.37084 12071.46232 8 -3700.86884 11696.37084 9 25027.12200 -3700.86884 10 -2173.23653 25027.12200 11 -32261.30065 -2173.23653 12 -22233.30128 -32261.30065 13 -9558.52325 -22233.30128 14 -4588.48851 -9558.52325 15 -23649.19944 -4588.48851 16 -18319.78623 -23649.19944 17 24298.90482 -18319.78623 18 23320.22095 24298.90482 19 22781.29675 23320.22095 20 14495.16514 22781.29675 21 -7002.41992 14495.16514 22 -8226.26911 -7002.41992 23 -14912.55878 -8226.26911 24 1691.68102 -14912.55878 25 -20534.23561 1691.68102 26 -2117.59936 -20534.23561 27 -9705.85260 -2117.59936 28 -24168.82126 -9705.85260 29 -12745.27790 -24168.82126 30 -25949.20119 -12745.27790 31 -6768.62149 -25949.20119 32 -4311.96334 -6768.62149 33 -29578.29108 -4311.96334 34 16865.86144 -29578.29108 35 -21116.14047 16865.86144 36 -11520.66697 -21116.14047 37 6958.93408 -11520.66697 38 -18816.00091 6958.93408 39 -645.52210 -18816.00091 40 -11209.11462 -645.52210 41 -13184.37475 -11209.11462 42 -11858.53404 -13184.37475 43 -20592.52130 -11858.53404 44 -26704.37502 -20592.52130 45 -14378.89819 -26704.37502 46 -6695.63307 -14378.89819 47 -13861.11403 -6695.63307 48 -6172.42326 -13861.11403 49 8557.31534 -6172.42326 50 26455.37045 8557.31534 51 -18360.29159 26455.37045 52 -2114.48728 -18360.29159 53 17584.88095 -2114.48728 54 -13655.42774 17584.88095 55 -16253.83483 -13655.42774 56 -5115.16399 -16253.83483 57 -3177.14353 -5115.16399 58 -3384.65006 -3177.14353 59 -5959.95249 -3384.65006 60 -2435.51326 -5959.95249 61 9241.28877 -2435.51326 62 27013.89593 9241.28877 63 -11126.28380 27013.89593 64 -17872.68381 -11126.28380 65 -23599.53510 -17872.68381 66 -7752.38017 -23599.53510 67 2624.43538 -7752.38017 68 -947.73671 2624.43538 69 -20229.29690 -947.73671 70 -6515.43292 -20229.29690 71 -5467.94012 -6515.43292 72 -22434.17430 -5467.94012 73 27177.26939 -22434.17430 74 -34024.72121 27177.26939 75 16993.83501 -34024.72121 76 2889.87194 16993.83501 77 -29065.35284 2889.87194 78 -15850.28126 -29065.35284 79 -20599.75939 -15850.28126 80 10353.12905 -20599.75939 81 -6032.36082 10353.12905 82 8164.73331 -6032.36082 83 -564.89400 8164.73331 84 132607.86690 -564.89400 85 -8279.90796 132607.86690 86 94.91248 -8279.90796 87 -6401.20169 94.91248 88 19793.66281 -6401.20169 89 -22244.59881 19793.66281 90 -9007.88347 -22244.59881 91 21187.35377 -9007.88347 92 -11900.64469 21187.35377 93 -940.98459 -11900.64469 94 -16549.97234 -940.98459 95 13285.71527 -16549.97234 96 -24331.71565 13285.71527 97 28785.71284 -24331.71565 98 -14847.49852 28785.71284 99 8446.32055 -14847.49852 100 -10283.15755 8446.32055 101 -9398.40806 -10283.15755 102 8554.95004 -9398.40806 103 34059.32689 8554.95004 104 25565.84324 34059.32689 105 35469.09788 25565.84324 106 -30507.82064 35469.09788 107 -5138.79986 -30507.82064 108 13983.00827 -5138.79986 109 -1048.15915 13983.00827 110 15781.24918 -1048.15915 111 -8012.45384 15781.24918 112 3027.04961 -8012.45384 113 -18611.55913 3027.04961 114 -15419.46178 -18611.55913 115 -5184.17275 -15419.46178 116 861.84933 -5184.17275 117 -7665.38091 861.84933 118 98184.92438 -7665.38091 119 -634.52131 98184.92438 120 -12111.92312 -634.52131 121 -9542.01725 -12111.92312 122 -21072.21950 -9542.01725 123 84741.84472 -21072.21950 124 17426.78698 84741.84472 125 -12351.45079 17426.78698 126 8283.61894 -12351.45079 127 1308.14094 8283.61894 128 16859.73758 1308.14094 129 -2449.42331 16859.73758 130 -3664.60148 -2449.42331 131 -9258.87716 -3664.60148 132 -5437.35369 -9258.87716 133 -11851.29816 -5437.35369 134 -15055.46292 -11851.29816 135 -23443.56577 -15055.46292 136 -1548.13237 -23443.56577 137 24580.09645 -1548.13237 138 24365.76431 24580.09645 139 -1079.80948 24365.76431 140 -7307.99819 -1079.80948 141 10764.44052 -7307.99819 142 2418.21319 10764.44052 143 -24147.68301 2418.21319 144 22239.65095 -24147.68301 145 569.21328 22239.65095 146 38372.42558 569.21328 147 38970.99939 38372.42558 148 -5404.52501 38970.99939 149 15166.67681 -5404.52501 150 -12616.36309 15166.67681 151 -8467.28978 -12616.36309 152 4115.81231 -8467.28978 153 -19681.42939 4115.81231 154 -6121.88619 -19681.42939 155 53340.09731 -6121.88619 156 -6542.74309 53340.09731 157 -3958.40085 -6542.74309 158 11020.67386 -3958.40085 159 7693.62901 11020.67386 160 -4153.20707 7693.62901 161 63728.41988 -4153.20707 162 4427.40501 63728.41988 163 35566.92673 4427.40501 164 -9002.38745 35566.92673 165 -25777.99910 -9002.38745 166 -21984.03050 -25777.99910 167 42444.28788 -21984.03050 168 -19637.54427 42444.28788 169 37621.72139 -19637.54427 170 18179.81442 37621.72139 171 4664.78334 18179.81442 172 17034.97488 4664.78334 173 -15150.37859 17034.97488 174 33575.13231 -15150.37859 175 -7968.80423 33575.13231 176 35454.85787 -7968.80423 177 -3304.71929 35454.85787 178 14719.39321 -3304.71929 179 -14948.81996 14719.39321 180 7013.26562 -14948.81996 181 25332.82318 7013.26562 182 755.45414 25332.82318 183 -2681.83425 755.45414 184 -35703.83763 -2681.83425 185 -1431.62282 -35703.83763 186 -5790.59168 -1431.62282 187 -2090.15439 -5790.59168 188 -13215.69904 -2090.15439 189 13169.75089 -13215.69904 190 24263.31846 13169.75089 191 2801.79890 24263.31846 192 -15583.24768 2801.79890 193 -15724.97028 -15583.24768 194 -5297.27155 -15724.97028 195 -9083.02604 -5297.27155 196 -6733.31607 -9083.02604 197 33261.69580 -6733.31607 198 23348.52976 33261.69580 199 10326.41207 23348.52976 200 -12620.01279 10326.41207 201 -15308.65792 -12620.01279 202 1890.17829 -15308.65792 203 -2174.88336 1890.17829 204 -11704.15866 -2174.88336 205 2222.73603 -11704.15866 206 -8153.06201 2222.73603 207 -157.95627 -8153.06201 208 -5690.30853 -157.95627 209 -19961.74832 -5690.30853 210 -17683.72874 -19961.74832 211 2796.23867 -17683.72874 212 -3387.29176 2796.23867 213 -11474.07479 -3387.29176 214 10993.74415 -11474.07479 215 -22352.22803 10993.74415 216 -4122.09586 -22352.22803 217 -893.35584 -4122.09586 218 -16804.38726 -893.35584 219 -6477.14820 -16804.38726 220 -6947.55255 -6477.14820 221 -16285.33372 -6947.55255 222 -18261.59819 -16285.33372 223 -4641.12939 -18261.59819 224 -8419.87186 -4641.12939 225 -15428.43762 -8419.87186 226 17035.82560 -15428.43762 227 1875.51234 17035.82560 228 -11354.07007 1875.51234 229 31268.10109 -11354.07007 230 -12410.76397 31268.10109 231 -7070.00647 -12410.76397 232 10513.00605 -7070.00647 233 -3131.84042 10513.00605 234 5215.27167 -3131.84042 235 -4084.41497 5215.27167 236 8449.33829 -4084.41497 237 -7845.58923 8449.33829 238 -10505.01416 -7845.58923 239 -3156.93538 -10505.01416 240 9112.12501 -3156.93538 241 -11111.42584 9112.12501 242 -8762.66804 -11111.42584 243 -1322.70843 -8762.66804 244 -2057.06271 -1322.70843 245 8081.08372 -2057.06271 246 5201.95444 8081.08372 247 1286.82579 5201.95444 248 17581.19517 1286.82579 249 15935.37064 17581.19517 250 -1243.23221 15935.37064 251 1225.17846 -1243.23221 252 -2668.30154 1225.17846 253 -1798.23473 -2668.30154 254 -9042.88249 -1798.23473 255 -349.49334 -9042.88249 256 4183.09042 -349.49334 257 5679.50345 4183.09042 258 -11435.07675 5679.50345 259 -4413.67007 -11435.07675 260 -2084.38125 -4413.67007 261 -16636.84538 -2084.38125 262 -2474.68687 -16636.84538 263 -7496.20363 -2474.68687 264 11465.12509 -7496.20363 265 -3683.28547 11465.12509 266 -11516.70696 -3683.28547 267 6376.99395 -11516.70696 268 -14045.98340 6376.99395 269 11484.69168 -14045.98340 270 -4753.54481 11484.69168 271 986.45543 -4753.54481 272 2560.43349 986.45543 273 -6221.52354 2560.43349 274 1196.40415 -6221.52354 275 -6791.87269 1196.40415 276 18524.46283 -6791.87269 277 -4890.61380 18524.46283 278 -16946.62533 -4890.61380 279 -7521.39460 -16946.62533 280 -1798.78534 -7521.39460 281 65266.10451 -1798.78534 282 12261.48287 65266.10451 283 -1217.24202 12261.48287 284 -2488.29780 -1217.24202 285 19563.18185 -2488.29780 286 -8910.56503 19563.18185 287 -8956.28845 -8910.56503 288 -6386.91361 -8956.28845 289 NA -6386.91361 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -4800.76920 5083.48536 [2,] -8043.89878 -4800.76920 [3,] 2375.70234 -8043.89878 [4,] -16436.07273 2375.70234 [5,] 19830.28558 -16436.07273 [6,] 12071.46232 19830.28558 [7,] 11696.37084 12071.46232 [8,] -3700.86884 11696.37084 [9,] 25027.12200 -3700.86884 [10,] -2173.23653 25027.12200 [11,] -32261.30065 -2173.23653 [12,] -22233.30128 -32261.30065 [13,] -9558.52325 -22233.30128 [14,] -4588.48851 -9558.52325 [15,] -23649.19944 -4588.48851 [16,] -18319.78623 -23649.19944 [17,] 24298.90482 -18319.78623 [18,] 23320.22095 24298.90482 [19,] 22781.29675 23320.22095 [20,] 14495.16514 22781.29675 [21,] -7002.41992 14495.16514 [22,] -8226.26911 -7002.41992 [23,] -14912.55878 -8226.26911 [24,] 1691.68102 -14912.55878 [25,] -20534.23561 1691.68102 [26,] -2117.59936 -20534.23561 [27,] -9705.85260 -2117.59936 [28,] -24168.82126 -9705.85260 [29,] -12745.27790 -24168.82126 [30,] -25949.20119 -12745.27790 [31,] -6768.62149 -25949.20119 [32,] -4311.96334 -6768.62149 [33,] -29578.29108 -4311.96334 [34,] 16865.86144 -29578.29108 [35,] -21116.14047 16865.86144 [36,] -11520.66697 -21116.14047 [37,] 6958.93408 -11520.66697 [38,] -18816.00091 6958.93408 [39,] -645.52210 -18816.00091 [40,] -11209.11462 -645.52210 [41,] -13184.37475 -11209.11462 [42,] -11858.53404 -13184.37475 [43,] -20592.52130 -11858.53404 [44,] -26704.37502 -20592.52130 [45,] -14378.89819 -26704.37502 [46,] -6695.63307 -14378.89819 [47,] -13861.11403 -6695.63307 [48,] -6172.42326 -13861.11403 [49,] 8557.31534 -6172.42326 [50,] 26455.37045 8557.31534 [51,] -18360.29159 26455.37045 [52,] -2114.48728 -18360.29159 [53,] 17584.88095 -2114.48728 [54,] -13655.42774 17584.88095 [55,] -16253.83483 -13655.42774 [56,] -5115.16399 -16253.83483 [57,] -3177.14353 -5115.16399 [58,] -3384.65006 -3177.14353 [59,] -5959.95249 -3384.65006 [60,] -2435.51326 -5959.95249 [61,] 9241.28877 -2435.51326 [62,] 27013.89593 9241.28877 [63,] -11126.28380 27013.89593 [64,] -17872.68381 -11126.28380 [65,] -23599.53510 -17872.68381 [66,] -7752.38017 -23599.53510 [67,] 2624.43538 -7752.38017 [68,] -947.73671 2624.43538 [69,] -20229.29690 -947.73671 [70,] -6515.43292 -20229.29690 [71,] -5467.94012 -6515.43292 [72,] -22434.17430 -5467.94012 [73,] 27177.26939 -22434.17430 [74,] -34024.72121 27177.26939 [75,] 16993.83501 -34024.72121 [76,] 2889.87194 16993.83501 [77,] -29065.35284 2889.87194 [78,] -15850.28126 -29065.35284 [79,] -20599.75939 -15850.28126 [80,] 10353.12905 -20599.75939 [81,] -6032.36082 10353.12905 [82,] 8164.73331 -6032.36082 [83,] -564.89400 8164.73331 [84,] 132607.86690 -564.89400 [85,] -8279.90796 132607.86690 [86,] 94.91248 -8279.90796 [87,] -6401.20169 94.91248 [88,] 19793.66281 -6401.20169 [89,] -22244.59881 19793.66281 [90,] -9007.88347 -22244.59881 [91,] 21187.35377 -9007.88347 [92,] -11900.64469 21187.35377 [93,] -940.98459 -11900.64469 [94,] -16549.97234 -940.98459 [95,] 13285.71527 -16549.97234 [96,] -24331.71565 13285.71527 [97,] 28785.71284 -24331.71565 [98,] -14847.49852 28785.71284 [99,] 8446.32055 -14847.49852 [100,] -10283.15755 8446.32055 [101,] -9398.40806 -10283.15755 [102,] 8554.95004 -9398.40806 [103,] 34059.32689 8554.95004 [104,] 25565.84324 34059.32689 [105,] 35469.09788 25565.84324 [106,] -30507.82064 35469.09788 [107,] -5138.79986 -30507.82064 [108,] 13983.00827 -5138.79986 [109,] -1048.15915 13983.00827 [110,] 15781.24918 -1048.15915 [111,] -8012.45384 15781.24918 [112,] 3027.04961 -8012.45384 [113,] -18611.55913 3027.04961 [114,] -15419.46178 -18611.55913 [115,] -5184.17275 -15419.46178 [116,] 861.84933 -5184.17275 [117,] -7665.38091 861.84933 [118,] 98184.92438 -7665.38091 [119,] -634.52131 98184.92438 [120,] -12111.92312 -634.52131 [121,] -9542.01725 -12111.92312 [122,] -21072.21950 -9542.01725 [123,] 84741.84472 -21072.21950 [124,] 17426.78698 84741.84472 [125,] -12351.45079 17426.78698 [126,] 8283.61894 -12351.45079 [127,] 1308.14094 8283.61894 [128,] 16859.73758 1308.14094 [129,] -2449.42331 16859.73758 [130,] -3664.60148 -2449.42331 [131,] -9258.87716 -3664.60148 [132,] -5437.35369 -9258.87716 [133,] -11851.29816 -5437.35369 [134,] -15055.46292 -11851.29816 [135,] -23443.56577 -15055.46292 [136,] -1548.13237 -23443.56577 [137,] 24580.09645 -1548.13237 [138,] 24365.76431 24580.09645 [139,] -1079.80948 24365.76431 [140,] -7307.99819 -1079.80948 [141,] 10764.44052 -7307.99819 [142,] 2418.21319 10764.44052 [143,] -24147.68301 2418.21319 [144,] 22239.65095 -24147.68301 [145,] 569.21328 22239.65095 [146,] 38372.42558 569.21328 [147,] 38970.99939 38372.42558 [148,] -5404.52501 38970.99939 [149,] 15166.67681 -5404.52501 [150,] -12616.36309 15166.67681 [151,] -8467.28978 -12616.36309 [152,] 4115.81231 -8467.28978 [153,] -19681.42939 4115.81231 [154,] -6121.88619 -19681.42939 [155,] 53340.09731 -6121.88619 [156,] -6542.74309 53340.09731 [157,] -3958.40085 -6542.74309 [158,] 11020.67386 -3958.40085 [159,] 7693.62901 11020.67386 [160,] -4153.20707 7693.62901 [161,] 63728.41988 -4153.20707 [162,] 4427.40501 63728.41988 [163,] 35566.92673 4427.40501 [164,] -9002.38745 35566.92673 [165,] -25777.99910 -9002.38745 [166,] -21984.03050 -25777.99910 [167,] 42444.28788 -21984.03050 [168,] -19637.54427 42444.28788 [169,] 37621.72139 -19637.54427 [170,] 18179.81442 37621.72139 [171,] 4664.78334 18179.81442 [172,] 17034.97488 4664.78334 [173,] -15150.37859 17034.97488 [174,] 33575.13231 -15150.37859 [175,] -7968.80423 33575.13231 [176,] 35454.85787 -7968.80423 [177,] -3304.71929 35454.85787 [178,] 14719.39321 -3304.71929 [179,] -14948.81996 14719.39321 [180,] 7013.26562 -14948.81996 [181,] 25332.82318 7013.26562 [182,] 755.45414 25332.82318 [183,] -2681.83425 755.45414 [184,] -35703.83763 -2681.83425 [185,] -1431.62282 -35703.83763 [186,] -5790.59168 -1431.62282 [187,] -2090.15439 -5790.59168 [188,] -13215.69904 -2090.15439 [189,] 13169.75089 -13215.69904 [190,] 24263.31846 13169.75089 [191,] 2801.79890 24263.31846 [192,] -15583.24768 2801.79890 [193,] -15724.97028 -15583.24768 [194,] -5297.27155 -15724.97028 [195,] -9083.02604 -5297.27155 [196,] -6733.31607 -9083.02604 [197,] 33261.69580 -6733.31607 [198,] 23348.52976 33261.69580 [199,] 10326.41207 23348.52976 [200,] -12620.01279 10326.41207 [201,] -15308.65792 -12620.01279 [202,] 1890.17829 -15308.65792 [203,] -2174.88336 1890.17829 [204,] -11704.15866 -2174.88336 [205,] 2222.73603 -11704.15866 [206,] -8153.06201 2222.73603 [207,] -157.95627 -8153.06201 [208,] -5690.30853 -157.95627 [209,] -19961.74832 -5690.30853 [210,] -17683.72874 -19961.74832 [211,] 2796.23867 -17683.72874 [212,] -3387.29176 2796.23867 [213,] -11474.07479 -3387.29176 [214,] 10993.74415 -11474.07479 [215,] -22352.22803 10993.74415 [216,] -4122.09586 -22352.22803 [217,] -893.35584 -4122.09586 [218,] -16804.38726 -893.35584 [219,] -6477.14820 -16804.38726 [220,] -6947.55255 -6477.14820 [221,] -16285.33372 -6947.55255 [222,] -18261.59819 -16285.33372 [223,] -4641.12939 -18261.59819 [224,] -8419.87186 -4641.12939 [225,] -15428.43762 -8419.87186 [226,] 17035.82560 -15428.43762 [227,] 1875.51234 17035.82560 [228,] -11354.07007 1875.51234 [229,] 31268.10109 -11354.07007 [230,] -12410.76397 31268.10109 [231,] -7070.00647 -12410.76397 [232,] 10513.00605 -7070.00647 [233,] -3131.84042 10513.00605 [234,] 5215.27167 -3131.84042 [235,] -4084.41497 5215.27167 [236,] 8449.33829 -4084.41497 [237,] -7845.58923 8449.33829 [238,] -10505.01416 -7845.58923 [239,] -3156.93538 -10505.01416 [240,] 9112.12501 -3156.93538 [241,] -11111.42584 9112.12501 [242,] -8762.66804 -11111.42584 [243,] -1322.70843 -8762.66804 [244,] -2057.06271 -1322.70843 [245,] 8081.08372 -2057.06271 [246,] 5201.95444 8081.08372 [247,] 1286.82579 5201.95444 [248,] 17581.19517 1286.82579 [249,] 15935.37064 17581.19517 [250,] -1243.23221 15935.37064 [251,] 1225.17846 -1243.23221 [252,] -2668.30154 1225.17846 [253,] -1798.23473 -2668.30154 [254,] -9042.88249 -1798.23473 [255,] -349.49334 -9042.88249 [256,] 4183.09042 -349.49334 [257,] 5679.50345 4183.09042 [258,] -11435.07675 5679.50345 [259,] -4413.67007 -11435.07675 [260,] -2084.38125 -4413.67007 [261,] -16636.84538 -2084.38125 [262,] -2474.68687 -16636.84538 [263,] -7496.20363 -2474.68687 [264,] 11465.12509 -7496.20363 [265,] -3683.28547 11465.12509 [266,] -11516.70696 -3683.28547 [267,] 6376.99395 -11516.70696 [268,] -14045.98340 6376.99395 [269,] 11484.69168 -14045.98340 [270,] -4753.54481 11484.69168 [271,] 986.45543 -4753.54481 [272,] 2560.43349 986.45543 [273,] -6221.52354 2560.43349 [274,] 1196.40415 -6221.52354 [275,] -6791.87269 1196.40415 [276,] 18524.46283 -6791.87269 [277,] -4890.61380 18524.46283 [278,] -16946.62533 -4890.61380 [279,] -7521.39460 -16946.62533 [280,] -1798.78534 -7521.39460 [281,] 65266.10451 -1798.78534 [282,] 12261.48287 65266.10451 [283,] -1217.24202 12261.48287 [284,] -2488.29780 -1217.24202 [285,] 19563.18185 -2488.29780 [286,] -8910.56503 19563.18185 [287,] -8956.28845 -8910.56503 [288,] -6386.91361 -8956.28845 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -4800.76920 5083.48536 2 -8043.89878 -4800.76920 3 2375.70234 -8043.89878 4 -16436.07273 2375.70234 5 19830.28558 -16436.07273 6 12071.46232 19830.28558 7 11696.37084 12071.46232 8 -3700.86884 11696.37084 9 25027.12200 -3700.86884 10 -2173.23653 25027.12200 11 -32261.30065 -2173.23653 12 -22233.30128 -32261.30065 13 -9558.52325 -22233.30128 14 -4588.48851 -9558.52325 15 -23649.19944 -4588.48851 16 -18319.78623 -23649.19944 17 24298.90482 -18319.78623 18 23320.22095 24298.90482 19 22781.29675 23320.22095 20 14495.16514 22781.29675 21 -7002.41992 14495.16514 22 -8226.26911 -7002.41992 23 -14912.55878 -8226.26911 24 1691.68102 -14912.55878 25 -20534.23561 1691.68102 26 -2117.59936 -20534.23561 27 -9705.85260 -2117.59936 28 -24168.82126 -9705.85260 29 -12745.27790 -24168.82126 30 -25949.20119 -12745.27790 31 -6768.62149 -25949.20119 32 -4311.96334 -6768.62149 33 -29578.29108 -4311.96334 34 16865.86144 -29578.29108 35 -21116.14047 16865.86144 36 -11520.66697 -21116.14047 37 6958.93408 -11520.66697 38 -18816.00091 6958.93408 39 -645.52210 -18816.00091 40 -11209.11462 -645.52210 41 -13184.37475 -11209.11462 42 -11858.53404 -13184.37475 43 -20592.52130 -11858.53404 44 -26704.37502 -20592.52130 45 -14378.89819 -26704.37502 46 -6695.63307 -14378.89819 47 -13861.11403 -6695.63307 48 -6172.42326 -13861.11403 49 8557.31534 -6172.42326 50 26455.37045 8557.31534 51 -18360.29159 26455.37045 52 -2114.48728 -18360.29159 53 17584.88095 -2114.48728 54 -13655.42774 17584.88095 55 -16253.83483 -13655.42774 56 -5115.16399 -16253.83483 57 -3177.14353 -5115.16399 58 -3384.65006 -3177.14353 59 -5959.95249 -3384.65006 60 -2435.51326 -5959.95249 61 9241.28877 -2435.51326 62 27013.89593 9241.28877 63 -11126.28380 27013.89593 64 -17872.68381 -11126.28380 65 -23599.53510 -17872.68381 66 -7752.38017 -23599.53510 67 2624.43538 -7752.38017 68 -947.73671 2624.43538 69 -20229.29690 -947.73671 70 -6515.43292 -20229.29690 71 -5467.94012 -6515.43292 72 -22434.17430 -5467.94012 73 27177.26939 -22434.17430 74 -34024.72121 27177.26939 75 16993.83501 -34024.72121 76 2889.87194 16993.83501 77 -29065.35284 2889.87194 78 -15850.28126 -29065.35284 79 -20599.75939 -15850.28126 80 10353.12905 -20599.75939 81 -6032.36082 10353.12905 82 8164.73331 -6032.36082 83 -564.89400 8164.73331 84 132607.86690 -564.89400 85 -8279.90796 132607.86690 86 94.91248 -8279.90796 87 -6401.20169 94.91248 88 19793.66281 -6401.20169 89 -22244.59881 19793.66281 90 -9007.88347 -22244.59881 91 21187.35377 -9007.88347 92 -11900.64469 21187.35377 93 -940.98459 -11900.64469 94 -16549.97234 -940.98459 95 13285.71527 -16549.97234 96 -24331.71565 13285.71527 97 28785.71284 -24331.71565 98 -14847.49852 28785.71284 99 8446.32055 -14847.49852 100 -10283.15755 8446.32055 101 -9398.40806 -10283.15755 102 8554.95004 -9398.40806 103 34059.32689 8554.95004 104 25565.84324 34059.32689 105 35469.09788 25565.84324 106 -30507.82064 35469.09788 107 -5138.79986 -30507.82064 108 13983.00827 -5138.79986 109 -1048.15915 13983.00827 110 15781.24918 -1048.15915 111 -8012.45384 15781.24918 112 3027.04961 -8012.45384 113 -18611.55913 3027.04961 114 -15419.46178 -18611.55913 115 -5184.17275 -15419.46178 116 861.84933 -5184.17275 117 -7665.38091 861.84933 118 98184.92438 -7665.38091 119 -634.52131 98184.92438 120 -12111.92312 -634.52131 121 -9542.01725 -12111.92312 122 -21072.21950 -9542.01725 123 84741.84472 -21072.21950 124 17426.78698 84741.84472 125 -12351.45079 17426.78698 126 8283.61894 -12351.45079 127 1308.14094 8283.61894 128 16859.73758 1308.14094 129 -2449.42331 16859.73758 130 -3664.60148 -2449.42331 131 -9258.87716 -3664.60148 132 -5437.35369 -9258.87716 133 -11851.29816 -5437.35369 134 -15055.46292 -11851.29816 135 -23443.56577 -15055.46292 136 -1548.13237 -23443.56577 137 24580.09645 -1548.13237 138 24365.76431 24580.09645 139 -1079.80948 24365.76431 140 -7307.99819 -1079.80948 141 10764.44052 -7307.99819 142 2418.21319 10764.44052 143 -24147.68301 2418.21319 144 22239.65095 -24147.68301 145 569.21328 22239.65095 146 38372.42558 569.21328 147 38970.99939 38372.42558 148 -5404.52501 38970.99939 149 15166.67681 -5404.52501 150 -12616.36309 15166.67681 151 -8467.28978 -12616.36309 152 4115.81231 -8467.28978 153 -19681.42939 4115.81231 154 -6121.88619 -19681.42939 155 53340.09731 -6121.88619 156 -6542.74309 53340.09731 157 -3958.40085 -6542.74309 158 11020.67386 -3958.40085 159 7693.62901 11020.67386 160 -4153.20707 7693.62901 161 63728.41988 -4153.20707 162 4427.40501 63728.41988 163 35566.92673 4427.40501 164 -9002.38745 35566.92673 165 -25777.99910 -9002.38745 166 -21984.03050 -25777.99910 167 42444.28788 -21984.03050 168 -19637.54427 42444.28788 169 37621.72139 -19637.54427 170 18179.81442 37621.72139 171 4664.78334 18179.81442 172 17034.97488 4664.78334 173 -15150.37859 17034.97488 174 33575.13231 -15150.37859 175 -7968.80423 33575.13231 176 35454.85787 -7968.80423 177 -3304.71929 35454.85787 178 14719.39321 -3304.71929 179 -14948.81996 14719.39321 180 7013.26562 -14948.81996 181 25332.82318 7013.26562 182 755.45414 25332.82318 183 -2681.83425 755.45414 184 -35703.83763 -2681.83425 185 -1431.62282 -35703.83763 186 -5790.59168 -1431.62282 187 -2090.15439 -5790.59168 188 -13215.69904 -2090.15439 189 13169.75089 -13215.69904 190 24263.31846 13169.75089 191 2801.79890 24263.31846 192 -15583.24768 2801.79890 193 -15724.97028 -15583.24768 194 -5297.27155 -15724.97028 195 -9083.02604 -5297.27155 196 -6733.31607 -9083.02604 197 33261.69580 -6733.31607 198 23348.52976 33261.69580 199 10326.41207 23348.52976 200 -12620.01279 10326.41207 201 -15308.65792 -12620.01279 202 1890.17829 -15308.65792 203 -2174.88336 1890.17829 204 -11704.15866 -2174.88336 205 2222.73603 -11704.15866 206 -8153.06201 2222.73603 207 -157.95627 -8153.06201 208 -5690.30853 -157.95627 209 -19961.74832 -5690.30853 210 -17683.72874 -19961.74832 211 2796.23867 -17683.72874 212 -3387.29176 2796.23867 213 -11474.07479 -3387.29176 214 10993.74415 -11474.07479 215 -22352.22803 10993.74415 216 -4122.09586 -22352.22803 217 -893.35584 -4122.09586 218 -16804.38726 -893.35584 219 -6477.14820 -16804.38726 220 -6947.55255 -6477.14820 221 -16285.33372 -6947.55255 222 -18261.59819 -16285.33372 223 -4641.12939 -18261.59819 224 -8419.87186 -4641.12939 225 -15428.43762 -8419.87186 226 17035.82560 -15428.43762 227 1875.51234 17035.82560 228 -11354.07007 1875.51234 229 31268.10109 -11354.07007 230 -12410.76397 31268.10109 231 -7070.00647 -12410.76397 232 10513.00605 -7070.00647 233 -3131.84042 10513.00605 234 5215.27167 -3131.84042 235 -4084.41497 5215.27167 236 8449.33829 -4084.41497 237 -7845.58923 8449.33829 238 -10505.01416 -7845.58923 239 -3156.93538 -10505.01416 240 9112.12501 -3156.93538 241 -11111.42584 9112.12501 242 -8762.66804 -11111.42584 243 -1322.70843 -8762.66804 244 -2057.06271 -1322.70843 245 8081.08372 -2057.06271 246 5201.95444 8081.08372 247 1286.82579 5201.95444 248 17581.19517 1286.82579 249 15935.37064 17581.19517 250 -1243.23221 15935.37064 251 1225.17846 -1243.23221 252 -2668.30154 1225.17846 253 -1798.23473 -2668.30154 254 -9042.88249 -1798.23473 255 -349.49334 -9042.88249 256 4183.09042 -349.49334 257 5679.50345 4183.09042 258 -11435.07675 5679.50345 259 -4413.67007 -11435.07675 260 -2084.38125 -4413.67007 261 -16636.84538 -2084.38125 262 -2474.68687 -16636.84538 263 -7496.20363 -2474.68687 264 11465.12509 -7496.20363 265 -3683.28547 11465.12509 266 -11516.70696 -3683.28547 267 6376.99395 -11516.70696 268 -14045.98340 6376.99395 269 11484.69168 -14045.98340 270 -4753.54481 11484.69168 271 986.45543 -4753.54481 272 2560.43349 986.45543 273 -6221.52354 2560.43349 274 1196.40415 -6221.52354 275 -6791.87269 1196.40415 276 18524.46283 -6791.87269 277 -4890.61380 18524.46283 278 -16946.62533 -4890.61380 279 -7521.39460 -16946.62533 280 -1798.78534 -7521.39460 281 65266.10451 -1798.78534 282 12261.48287 65266.10451 283 -1217.24202 12261.48287 284 -2488.29780 -1217.24202 285 19563.18185 -2488.29780 286 -8910.56503 19563.18185 287 -8956.28845 -8910.56503 288 -6386.91361 -8956.28845 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/713u81323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/80m911323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/9a4431323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/rcomp/tmp/10fqu91323768015.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/11xgfb1323768015.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/127cez1323768015.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/133m271323768015.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/14ehng1323768015.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/rcomp/tmp/15khok1323768015.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/rcomp/tmp/164yle1323768015.tab") + } > > try(system("convert tmp/10nv31323768015.ps tmp/10nv31323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/2tamq1323768015.ps tmp/2tamq1323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/3t7m61323768015.ps tmp/3t7m61323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/4dqw21323768015.ps tmp/4dqw21323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/5sr731323768015.ps tmp/5sr731323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/64sef1323768015.ps tmp/64sef1323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/713u81323768015.ps tmp/713u81323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/80m911323768015.ps tmp/80m911323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/9a4431323768015.ps tmp/9a4431323768015.png",intern=TRUE)) character(0) > try(system("convert tmp/10fqu91323768015.ps tmp/10fqu91323768015.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 16.348 0.680 17.586