R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(56 + ,3 + ,30 + ,112285 + ,24188 + ,146283 + ,144 + ,145 + ,1 + ,56 + ,4 + ,28 + ,84786 + ,18273 + ,98364 + ,103 + ,101 + ,1 + ,54 + ,12 + ,38 + ,83123 + ,14130 + ,86146 + ,98 + ,98 + ,1 + ,89 + ,2 + ,30 + ,101193 + ,32287 + ,96933 + ,135 + ,132 + ,1 + ,40 + ,1 + ,22 + ,38361 + ,8654 + ,79234 + ,61 + ,60 + ,1 + ,25 + ,3 + ,26 + ,68504 + ,9245 + ,42551 + ,39 + ,38 + ,1 + ,92 + ,0 + ,25 + ,119182 + ,33251 + ,195663 + ,150 + ,144 + ,1 + ,18 + ,0 + ,18 + ,22807 + ,1271 + ,6853 + ,5 + ,5 + ,1 + ,63 + ,0 + ,11 + ,17140 + ,5279 + ,21529 + ,28 + ,28 + ,0 + ,44 + ,5 + ,26 + ,116174 + ,27101 + ,95757 + ,84 + ,84 + ,1 + ,33 + ,0 + ,25 + ,57635 + ,16373 + ,85584 + ,80 + ,79 + ,1 + ,84 + ,0 + ,38 + ,66198 + ,19716 + ,143983 + ,130 + ,127 + ,1 + ,88 + ,7 + ,44 + ,71701 + ,17753 + ,75851 + ,82 + ,78 + ,1 + ,55 + ,7 + ,30 + ,57793 + ,9028 + ,59238 + ,60 + ,60 + ,1 + ,60 + ,3 + ,40 + ,80444 + ,18653 + ,93163 + ,131 + ,131 + ,1 + ,66 + ,9 + ,34 + ,53855 + ,8828 + ,96037 + ,84 + ,84 + ,1 + ,154 + ,0 + ,47 + ,97668 + ,29498 + ,151511 + ,140 + ,133 + ,1 + ,53 + ,4 + ,30 + ,133824 + ,27563 + ,136368 + ,151 + ,150 + ,1 + ,119 + ,3 + ,31 + ,101481 + ,18293 + ,112642 + ,91 + ,91 + ,1 + ,41 + ,0 + ,23 + ,99645 + ,22530 + ,94728 + ,138 + ,132 + ,1 + ,61 + ,7 + ,36 + ,114789 + ,15977 + ,105499 + ,150 + ,136 + ,1 + ,58 + ,0 + ,36 + ,99052 + ,35082 + ,121527 + ,124 + ,124 + ,1 + ,75 + ,1 + ,30 + ,67654 + ,16116 + ,127766 + ,119 + ,118 + ,1 + ,33 + ,5 + ,25 + ,65553 + ,15849 + ,98958 + ,73 + ,70 + ,1 + ,40 + ,7 + ,39 + ,97500 + ,16026 + ,77900 + ,110 + ,107 + ,1 + ,92 + ,0 + ,34 + ,69112 + ,26569 + ,85646 + ,123 + ,119 + ,1 + ,100 + ,0 + ,31 + ,82753 + ,24785 + ,98579 + ,90 + ,89 + ,1 + ,112 + ,5 + ,31 + ,85323 + ,17569 + ,130767 + ,116 + ,112 + ,1 + ,73 + ,0 + ,33 + ,72654 + ,23825 + ,131741 + ,113 + ,108 + ,1 + ,40 + ,0 + ,25 + ,30727 + ,7869 + ,53907 + ,56 + ,52 + ,1 + ,45 + ,0 + ,33 + ,77873 + ,14975 + ,178812 + ,115 + ,112 + ,1 + ,60 + ,3 + ,35 + ,117478 + ,37791 + ,146761 + ,119 + ,116 + ,1 + ,62 + ,4 + ,42 + ,74007 + ,9605 + ,82036 + ,129 + ,123 + ,1 + ,75 + ,1 + ,43 + ,90183 + ,27295 + ,163253 + ,127 + ,125 + ,1 + ,31 + ,4 + ,30 + ,61542 + ,2746 + ,27032 + ,27 + ,27 + ,1 + ,77 + ,2 + ,33 + ,101494 + ,34461 + ,171975 + ,175 + ,162 + ,1 + ,34 + ,0 + ,13 + ,27570 + ,8098 + ,65990 + ,35 + ,32 + ,0 + ,46 + ,0 + ,32 + ,55813 + ,4787 + ,86572 + ,64 + ,64 + ,1 + ,99 + ,0 + ,36 + ,79215 + ,24919 + ,159676 + ,96 + ,92 + ,1 + ,17 + ,0 + ,0 + ,1423 + ,603 + ,1929 + ,0 + ,0 + ,0 + ,66 + ,2 + ,28 + ,55461 + ,16329 + ,85371 + ,84 + ,83 + ,1 + ,30 + ,1 + ,14 + ,31081 + ,12558 + ,58391 + ,41 + ,41 + ,1 + ,76 + ,0 + ,17 + ,22996 + ,7784 + ,31580 + ,47 + ,47 + ,0 + ,146 + ,2 + ,32 + ,83122 + ,28522 + ,136815 + ,126 + ,120 + ,1 + ,67 + ,10 + ,30 + ,70106 + ,22265 + ,120642 + ,105 + ,105 + ,1 + ,56 + ,6 + ,35 + ,60578 + ,14459 + ,69107 + ,80 + ,79 + ,1 + ,107 + ,0 + ,20 + ,39992 + ,14526 + ,50495 + ,70 + ,65 + ,0 + ,58 + ,5 + ,28 + ,79892 + ,22240 + ,108016 + ,73 + ,70 + ,1 + ,34 + ,4 + ,28 + ,49810 + ,11802 + ,46341 + ,57 + ,55 + ,1 + ,61 + ,1 + ,39 + ,71570 + ,7623 + ,78348 + ,40 + ,39 + ,1 + ,119 + ,2 + ,34 + ,100708 + ,11912 + ,79336 + ,68 + ,67 + ,1 + ,42 + ,2 + ,26 + ,33032 + ,7935 + ,56968 + ,21 + ,21 + ,1 + ,66 + ,0 + ,39 + ,82875 + ,18220 + ,93176 + ,127 + ,127 + ,1 + ,89 + ,8 + ,39 + ,139077 + ,19199 + ,161632 + ,154 + ,152 + ,1 + ,44 + ,3 + ,33 + ,71595 + ,19918 + ,87850 + ,116 + ,113 + ,1 + ,66 + ,0 + ,28 + ,72260 + ,21884 + ,127969 + ,102 + ,99 + ,1 + ,24 + ,0 + ,4 + ,5950 + ,2694 + ,15049 + ,7 + ,7 + ,1 + ,259 + ,8 + ,39 + ,115762 + ,15808 + ,155135 + ,148 + ,141 + ,1 + ,17 + ,5 + ,18 + ,32551 + ,3597 + ,25109 + ,21 + ,21 + ,1 + ,64 + ,3 + ,14 + ,31701 + ,5296 + ,45824 + ,35 + ,35 + ,1 + ,41 + ,1 + ,29 + ,80670 + ,25239 + ,102996 + ,112 + ,109 + ,1 + ,68 + ,5 + ,44 + ,143558 + ,29801 + ,160604 + ,137 + ,133 + ,1 + ,168 + ,1 + ,21 + ,117105 + ,18450 + ,158051 + ,135 + ,123 + ,0 + ,43 + ,1 + ,16 + ,23789 + ,7132 + ,44547 + ,26 + ,26 + ,0 + ,132 + ,5 + ,28 + ,120733 + ,34861 + ,162647 + ,230 + ,230 + ,1 + ,105 + ,0 + ,35 + ,105195 + ,35940 + ,174141 + ,181 + ,166 + ,1 + ,71 + ,12 + ,28 + ,73107 + ,16688 + ,60622 + ,71 + ,68 + ,1 + ,112 + ,8 + ,38 + ,132068 + ,24683 + ,179566 + ,147 + ,147 + ,1 + ,94 + ,8 + ,23 + ,149193 + ,46230 + ,184301 + ,190 + ,179 + ,1 + ,82 + ,8 + ,36 + ,46821 + ,10387 + ,75661 + ,64 + ,61 + ,1 + ,70 + ,8 + ,32 + ,87011 + ,21436 + ,96144 + ,105 + ,101 + ,1 + ,57 + ,2 + ,29 + ,95260 + ,30546 + ,129847 + ,107 + ,108 + ,1 + ,53 + ,0 + ,25 + ,55183 + ,19746 + ,117286 + ,94 + ,90 + ,1 + ,103 + ,5 + ,27 + ,106671 + ,15977 + ,71180 + ,116 + ,114 + ,1 + ,121 + ,8 + ,36 + ,73511 + ,22583 + ,109377 + ,106 + ,103 + ,1 + ,62 + ,2 + ,28 + ,92945 + ,17274 + ,85298 + ,143 + ,142 + ,1 + ,52 + ,5 + ,23 + ,78664 + ,16469 + ,73631 + ,81 + ,79 + ,1 + ,52 + ,12 + ,40 + ,70054 + ,14251 + ,86767 + ,89 + ,88 + ,1 + ,32 + ,6 + ,23 + ,22618 + ,3007 + ,23824 + ,26 + ,25 + ,1 + ,62 + ,7 + ,40 + ,74011 + ,16851 + ,93487 + ,84 + ,83 + ,1 + ,45 + ,2 + ,28 + ,83737 + ,21113 + ,82981 + ,113 + ,113 + ,1 + ,46 + ,0 + ,34 + ,69094 + ,17401 + ,73815 + ,120 + ,118 + ,1 + ,63 + ,4 + ,33 + ,93133 + ,23958 + ,94552 + ,110 + ,110 + ,1 + ,75 + ,3 + ,28 + ,95536 + ,23567 + ,132190 + ,134 + ,129 + ,1 + ,88 + ,6 + ,34 + ,225920 + ,13065 + ,128754 + ,54 + ,51 + ,1 + ,46 + ,2 + ,30 + ,62133 + ,15358 + ,66363 + ,96 + ,93 + ,1 + ,53 + ,0 + ,33 + ,61370 + ,14587 + ,67808 + ,78 + ,76 + ,1 + ,37 + ,1 + ,22 + ,43836 + ,12770 + ,61724 + ,51 + ,49 + ,1 + ,90 + ,0 + ,38 + ,106117 + ,24021 + ,131722 + ,121 + ,118 + ,1 + ,63 + ,5 + ,26 + ,38692 + ,9648 + ,68580 + ,38 + ,38 + ,1 + ,78 + ,2 + ,35 + ,84651 + ,20537 + ,106175 + ,145 + ,141 + ,1 + ,25 + ,0 + ,8 + ,56622 + ,7905 + ,55792 + ,59 + ,58 + ,1 + ,45 + ,0 + ,24 + ,15986 + ,4527 + ,25157 + ,27 + ,27 + ,1 + ,46 + ,5 + ,29 + ,95364 + ,30495 + ,76669 + ,91 + ,91 + ,1 + ,41 + ,0 + ,20 + ,26706 + ,7117 + ,57283 + ,48 + ,48 + ,0 + ,144 + ,1 + ,29 + ,89691 + ,17719 + ,105805 + ,68 + ,63 + ,1 + ,82 + ,0 + ,45 + ,67267 + ,27056 + ,129484 + ,58 + ,56 + ,1 + ,91 + ,1 + ,37 + ,126846 + ,33473 + ,72413 + ,150 + ,144 + ,1 + ,71 + ,1 + ,33 + ,41140 + ,9758 + ,87831 + ,74 + ,73 + ,1 + ,63 + ,2 + ,33 + ,102860 + ,21115 + ,96971 + ,181 + ,168 + ,1 + ,53 + ,6 + ,25 + ,51715 + ,7236 + ,71299 + ,65 + ,64 + ,1 + ,62 + ,1 + ,32 + ,55801 + ,13790 + ,77494 + ,97 + ,97 + ,1 + ,63 + ,4 + ,29 + ,111813 + ,32902 + ,120336 + ,121 + ,117 + ,1 + ,32 + ,2 + ,28 + ,120293 + ,25131 + ,93913 + ,99 + ,100 + ,1 + ,39 + ,3 + ,28 + ,138599 + ,30910 + ,136048 + ,152 + ,149 + ,1 + ,62 + ,0 + ,31 + ,161647 + ,35947 + ,181248 + ,188 + ,187 + ,1 + ,117 + ,10 + ,52 + ,115929 + ,29848 + ,146123 + ,138 + ,127 + ,1 + ,34 + ,0 + ,21 + ,24266 + ,6943 + ,32036 + ,40 + ,37 + ,0 + ,92 + ,9 + ,24 + ,162901 + ,42705 + ,186646 + ,254 + ,245 + ,1 + ,93 + ,7 + ,41 + ,109825 + ,31808 + ,102255 + ,87 + ,87 + ,1 + ,54 + ,0 + ,33 + ,129838 + ,26675 + ,168237 + ,178 + ,177 + ,1 + ,144 + ,0 + ,32 + ,37510 + ,8435 + ,64219 + ,51 + ,49 + ,1 + ,14 + ,4 + ,19 + ,43750 + ,7409 + ,19630 + ,49 + ,49 + ,1 + ,61 + ,4 + ,20 + ,40652 + ,14993 + ,76825 + ,73 + ,73 + ,1 + ,109 + ,0 + ,31 + ,87771 + ,36867 + ,115338 + ,176 + ,177 + ,1 + ,38 + ,0 + ,31 + ,85872 + ,33835 + ,109427 + ,94 + ,94 + ,1 + ,73 + ,0 + ,32 + ,89275 + ,24164 + ,118168 + ,120 + ,117 + ,1 + ,75 + ,1 + ,18 + ,44418 + ,12607 + ,84845 + ,66 + ,60 + ,0 + ,50 + ,0 + ,23 + ,192565 + ,22609 + ,153197 + ,56 + ,55 + ,1 + ,61 + ,1 + ,17 + ,35232 + ,5892 + ,29877 + ,39 + ,39 + ,0 + ,55 + ,0 + ,20 + ,40909 + ,17014 + ,63506 + ,66 + ,64 + ,0 + ,77 + ,0 + ,12 + ,13294 + ,5394 + ,22445 + ,27 + ,26 + ,0 + ,75 + ,4 + ,17 + ,32387 + ,9178 + ,47695 + ,65 + ,64 + ,0 + ,72 + ,0 + ,30 + ,140867 + ,6440 + ,68370 + ,58 + ,58 + ,1 + ,50 + ,4 + ,31 + ,120662 + ,21916 + ,146304 + ,98 + ,95 + ,1 + ,32 + ,4 + ,10 + ,21233 + ,4011 + ,38233 + ,25 + ,25 + ,0 + ,53 + ,3 + ,13 + ,44332 + ,5818 + ,42071 + ,26 + ,26 + ,0 + ,42 + ,0 + ,22 + ,61056 + ,18647 + ,50517 + ,77 + ,76 + ,0 + ,71 + ,0 + ,42 + ,101338 + ,20556 + ,103950 + ,130 + ,129 + ,1 + ,10 + ,0 + ,1 + ,1168 + ,238 + ,5841 + ,11 + ,11 + ,1 + ,35 + ,5 + ,9 + ,13497 + ,70 + ,2341 + ,2 + ,2 + ,0 + ,65 + ,0 + ,32 + ,65567 + ,22392 + ,84396 + ,101 + ,101 + ,1 + ,25 + ,4 + ,11 + ,25162 + ,3913 + ,24610 + ,31 + ,28 + ,1 + ,66 + ,0 + ,25 + ,32334 + ,12237 + ,35753 + ,36 + ,36 + ,0 + ,41 + ,0 + ,36 + ,40735 + ,8388 + ,55515 + ,120 + ,89 + ,1 + ,86 + ,1 + ,31 + ,91413 + ,22120 + ,209056 + ,195 + ,193 + ,1 + ,16 + ,0 + ,0 + ,855 + ,338 + ,6622 + ,4 + ,4 + ,1 + ,42 + ,5 + ,24 + ,97068 + ,11727 + ,115814 + ,89 + ,84 + ,1 + ,19 + ,0 + ,13 + ,44339 + ,3704 + ,11609 + ,24 + ,23 + ,0 + ,19 + ,0 + ,8 + ,14116 + ,3988 + ,13155 + ,39 + ,39 + ,1 + ,45 + ,0 + ,13 + ,10288 + ,3030 + ,18274 + ,14 + ,14 + ,0 + ,65 + ,0 + ,19 + ,65622 + ,13520 + ,72875 + ,78 + ,78 + ,0 + ,35 + ,0 + ,18 + ,16563 + ,1421 + ,10112 + ,15 + ,14 + ,0 + ,95 + ,2 + ,33 + ,76643 + ,20923 + ,142775 + ,106 + ,101 + ,1 + ,49 + ,7 + ,40 + ,110681 + ,20237 + ,68847 + ,83 + ,82 + ,1 + ,37 + ,1 + ,22 + ,29011 + ,3219 + ,17659 + ,24 + ,24 + ,0 + ,64 + ,8 + ,38 + ,92696 + ,3769 + ,20112 + ,37 + ,36 + ,1 + ,38 + ,2 + ,24 + ,94785 + ,12252 + ,61023 + ,77 + ,75 + ,1 + ,34 + ,0 + ,8 + ,8773 + ,1888 + ,13983 + ,16 + ,16 + ,1 + ,32 + ,2 + ,35 + ,83209 + ,14497 + ,65176 + ,56 + ,55 + ,1 + ,65 + ,0 + ,43 + ,93815 + ,28864 + ,132432 + ,132 + ,131 + ,1 + ,52 + ,0 + ,43 + ,86687 + ,21721 + ,112494 + ,144 + ,131 + ,1 + ,62 + ,1 + ,14 + ,34553 + ,4821 + ,45109 + ,40 + ,39 + ,0 + ,65 + ,3 + ,41 + ,105547 + ,33644 + ,170875 + ,153 + ,144 + ,1 + ,83 + ,0 + ,38 + ,103487 + ,15923 + ,180759 + ,143 + ,139 + ,1 + ,95 + ,3 + ,45 + ,213688 + ,42935 + ,214921 + ,220 + ,211 + ,1 + ,29 + ,0 + ,31 + ,71220 + ,18864 + ,100226 + ,79 + ,78 + ,1 + ,18 + ,0 + ,13 + ,23517 + ,4977 + ,32043 + ,50 + ,50 + ,0 + ,33 + ,0 + ,28 + ,56926 + ,7785 + ,54454 + ,39 + ,39 + ,1 + ,247 + ,4 + ,31 + ,91721 + ,17939 + ,78876 + ,95 + ,90 + ,1 + ,139 + ,4 + ,40 + ,115168 + ,23436 + ,170745 + ,169 + ,166 + ,1 + ,29 + ,11 + ,30 + ,111194 + ,325 + ,6940 + ,12 + ,12 + ,1 + ,118 + ,0 + ,16 + ,51009 + ,13539 + ,49025 + ,63 + ,57 + ,0 + ,110 + ,0 + ,37 + ,135777 + ,34538 + ,122037 + ,134 + ,133 + ,1 + ,67 + ,4 + ,30 + ,51513 + ,12198 + ,53782 + ,69 + ,69 + ,1 + ,42 + ,0 + ,35 + ,74163 + ,26924 + ,127748 + ,119 + ,119 + ,1 + ,65 + ,1 + ,32 + ,51633 + ,12716 + ,86839 + ,119 + ,119 + ,1 + ,94 + ,0 + ,27 + ,75345 + ,8172 + ,44830 + ,75 + ,65 + ,1 + ,64 + ,0 + ,20 + ,33416 + ,10855 + ,77395 + ,63 + ,61 + ,0 + ,81 + ,0 + ,18 + ,83305 + ,11932 + ,89324 + ,55 + ,49 + ,0 + ,95 + ,9 + ,31 + ,98952 + ,14300 + ,103300 + ,103 + ,101 + ,1 + ,67 + ,1 + ,31 + ,102372 + ,25515 + ,112283 + ,197 + ,196 + ,1 + ,63 + ,3 + ,21 + ,37238 + ,2805 + ,10901 + ,16 + ,15 + ,1 + ,83 + ,10 + ,39 + ,103772 + ,29402 + ,120691 + ,140 + ,136 + ,1 + ,45 + ,5 + ,41 + ,123969 + ,16440 + ,58106 + ,89 + ,89 + ,1 + ,30 + ,0 + ,13 + ,27142 + ,11221 + ,57140 + ,40 + ,40 + ,0 + ,70 + ,2 + ,32 + ,135400 + ,28732 + ,122422 + ,125 + ,123 + ,1 + ,32 + ,0 + ,18 + ,21399 + ,5250 + ,25899 + ,21 + ,21 + ,0 + ,83 + ,1 + ,39 + ,130115 + ,28608 + ,139296 + ,167 + ,163 + ,1 + ,31 + ,2 + ,14 + ,24874 + ,8092 + ,52678 + ,32 + ,29 + ,0 + ,67 + ,4 + ,7 + ,34988 + ,4473 + ,23853 + ,36 + ,35 + ,0 + ,66 + ,0 + ,17 + ,45549 + ,1572 + ,17306 + ,13 + ,13 + ,0 + ,10 + ,0 + ,0 + ,6023 + ,2065 + ,7953 + ,5 + ,5 + ,1 + ,70 + ,2 + ,30 + ,64466 + ,14817 + ,89455 + ,96 + ,96 + ,1 + ,103 + ,1 + ,37 + ,54990 + ,16714 + ,147866 + ,151 + ,151 + ,1 + ,5 + ,0 + ,0 + ,1644 + ,556 + ,4245 + ,6 + ,6 + ,1 + ,20 + ,0 + ,5 + ,6179 + ,2089 + ,21509 + ,13 + ,13 + ,1 + ,5 + ,0 + ,1 + ,3926 + ,2658 + ,7670 + ,3 + ,3 + ,1 + ,36 + ,1 + ,16 + ,32755 + ,10695 + ,66675 + ,57 + ,56 + ,0 + ,34 + ,0 + ,32 + ,34777 + ,1669 + ,14336 + ,23 + ,23 + ,1 + ,48 + ,2 + ,24 + ,73224 + ,16267 + ,53608 + ,61 + ,57 + ,1 + ,40 + ,0 + ,17 + ,27114 + ,7768 + ,30059 + ,21 + ,14 + ,0 + ,43 + ,3 + ,11 + ,20760 + ,7252 + ,29668 + ,43 + ,43 + ,0 + ,31 + ,6 + ,24 + ,37636 + ,6387 + ,22097 + ,20 + ,20 + ,0 + ,42 + ,0 + ,22 + ,65461 + ,18715 + ,96841 + ,82 + ,72 + ,0 + ,46 + ,2 + ,12 + ,30080 + ,7936 + ,41907 + ,90 + ,87 + ,0 + ,33 + ,0 + ,19 + ,24094 + ,8643 + ,27080 + ,25 + ,21 + ,0 + ,18 + ,2 + ,13 + ,69008 + ,7294 + ,35885 + ,60 + ,56 + ,0 + ,55 + ,1 + ,17 + ,54968 + ,4570 + ,41247 + ,61 + ,59 + ,0 + ,35 + ,1 + ,15 + ,46090 + ,7185 + ,28313 + ,85 + ,82 + ,0 + ,59 + ,2 + ,16 + ,27507 + ,10058 + ,36845 + ,43 + ,43 + ,0 + ,19 + ,1 + ,24 + ,10672 + ,2342 + ,16548 + ,25 + ,25 + ,0 + ,66 + ,0 + ,15 + ,34029 + ,8509 + ,36134 + ,41 + ,38 + ,0 + ,60 + ,1 + ,17 + ,46300 + ,13275 + ,55764 + ,26 + ,25 + ,0 + ,36 + ,3 + ,18 + ,24760 + ,6816 + ,28910 + ,38 + ,38 + ,0 + ,25 + ,0 + ,20 + ,18779 + ,1930 + ,13339 + ,12 + ,12 + ,0 + ,47 + ,0 + ,16 + ,21280 + ,8086 + ,25319 + ,29 + ,29 + ,0 + ,54 + ,0 + ,16 + ,40662 + ,10737 + ,66956 + ,49 + ,47 + ,0 + ,53 + ,0 + ,18 + ,28987 + ,8033 + ,47487 + ,46 + ,45 + ,0 + ,40 + ,1 + ,22 + ,22827 + ,7058 + ,52785 + ,41 + ,40 + ,0 + ,40 + ,4 + ,8 + ,18513 + ,6782 + ,44683 + ,31 + ,30 + ,0 + ,39 + ,0 + ,17 + ,30594 + ,5401 + ,35619 + ,41 + ,41 + ,0 + ,14 + ,0 + ,18 + ,24006 + ,6521 + ,21920 + ,26 + ,25 + ,0 + ,45 + ,0 + ,16 + ,27913 + ,10856 + ,45608 + ,23 + ,23 + ,0 + ,36 + ,7 + ,23 + ,42744 + ,2154 + ,7721 + ,14 + ,14 + ,0 + ,28 + ,2 + ,22 + ,12934 + ,6117 + ,20634 + ,16 + ,16 + ,0 + ,44 + ,0 + ,13 + ,22574 + ,5238 + ,29788 + ,25 + ,26 + ,0 + ,30 + ,7 + ,13 + ,41385 + ,4820 + ,31931 + ,21 + ,21 + ,0 + ,22 + ,3 + ,16 + ,18653 + ,5615 + ,37754 + ,32 + ,27 + ,0 + ,17 + ,0 + ,16 + ,18472 + ,4272 + ,32505 + ,9 + ,9 + ,0 + ,31 + ,0 + ,20 + ,30976 + ,8702 + ,40557 + ,35 + ,33 + ,0 + ,55 + ,6 + ,22 + ,63339 + ,15340 + ,94238 + ,42 + ,42 + ,0 + ,54 + ,2 + ,17 + ,25568 + ,8030 + ,44197 + ,68 + ,68 + ,0 + ,21 + ,0 + ,18 + ,33747 + ,9526 + ,43228 + ,32 + ,32 + ,0 + ,14 + ,0 + ,17 + ,4154 + ,1278 + ,4103 + ,6 + ,6 + ,0 + ,81 + ,3 + ,12 + ,19474 + ,4236 + ,44144 + ,68 + ,67 + ,0 + ,35 + ,0 + ,7 + ,35130 + ,3023 + ,32868 + ,33 + ,33 + ,0 + ,43 + ,1 + ,17 + ,39067 + ,7196 + ,27640 + ,84 + ,77 + ,0 + ,46 + ,1 + ,14 + ,13310 + ,3394 + ,14063 + ,46 + ,46 + ,0 + ,30 + ,0 + ,23 + ,65892 + ,6371 + ,28990 + ,30 + ,30 + ,0 + ,23 + ,1 + ,17 + ,4143 + ,1574 + ,4694 + ,0 + ,0 + ,0 + ,38 + ,0 + ,14 + ,28579 + ,9620 + ,42648 + ,36 + ,36 + ,0 + ,54 + ,0 + ,15 + ,51776 + ,6978 + ,64329 + ,47 + ,46 + ,0 + ,20 + ,0 + ,17 + ,21152 + ,4911 + ,21928 + ,20 + ,18 + ,0 + ,53 + ,0 + ,21 + ,38084 + ,8645 + ,25836 + ,50 + ,48 + ,0 + ,45 + ,0 + ,18 + ,27717 + ,8987 + ,22779 + ,30 + ,29 + ,0 + ,39 + ,0 + ,18 + ,32928 + ,5544 + ,40820 + ,30 + ,28 + ,0 + ,20 + ,0 + ,17 + ,11342 + ,3083 + ,27530 + ,34 + ,34 + ,0 + ,24 + ,0 + ,17 + ,19499 + ,6909 + ,32378 + ,33 + ,33 + ,0 + ,31 + ,0 + ,16 + ,16380 + ,3189 + ,10824 + ,34 + ,34 + ,0 + ,35 + ,0 + ,15 + ,36874 + ,6745 + ,39613 + ,37 + ,33 + ,0 + ,151 + ,0 + ,21 + ,48259 + ,16724 + ,60865 + ,83 + ,80 + ,0 + ,52 + ,0 + ,16 + ,16734 + ,4850 + ,19787 + ,32 + ,32 + ,0 + ,30 + ,2 + ,14 + ,28207 + ,7025 + ,20107 + ,30 + ,30 + ,0 + ,31 + ,0 + ,15 + ,30143 + ,6047 + ,36605 + ,43 + ,41 + ,0 + ,29 + ,1 + ,17 + ,41369 + ,7377 + ,40961 + ,41 + ,41 + ,0 + ,57 + ,1 + ,15 + ,45833 + ,9078 + ,48231 + ,51 + ,51 + ,0 + ,40 + ,0 + ,15 + ,29156 + ,4605 + ,39725 + ,19 + ,18 + ,0 + ,44 + ,0 + ,10 + ,35944 + ,3238 + ,21455 + ,37 + ,34 + ,0 + ,25 + ,0 + ,6 + ,36278 + ,8100 + ,23430 + ,33 + ,31 + ,0 + ,77 + ,0 + ,22 + ,45588 + ,9653 + ,62991 + ,41 + ,39 + ,0 + ,35 + ,0 + ,21 + ,45097 + ,8914 + ,49363 + ,54 + ,54 + ,0 + ,11 + ,0 + ,1 + ,3895 + ,786 + ,9604 + ,14 + ,14 + ,0 + ,63 + ,1 + ,18 + ,28394 + ,6700 + ,24552 + ,25 + ,24 + ,0 + ,44 + ,0 + ,17 + ,18632 + ,5788 + ,31493 + ,25 + ,24 + ,0 + ,19 + ,0 + ,4 + ,2325 + ,593 + ,3439 + ,8 + ,8 + ,0 + ,13 + ,0 + ,10 + ,25139 + ,4506 + ,19555 + ,26 + ,26 + ,0 + ,42 + ,0 + ,16 + ,27975 + ,6382 + ,21228 + ,20 + ,19 + ,0 + ,38 + ,1 + ,16 + ,14483 + ,5621 + ,23177 + ,11 + ,11 + ,0 + ,29 + ,0 + ,9 + ,13127 + ,3997 + ,22094 + ,14 + ,14 + ,0 + ,20 + ,0 + ,16 + ,5839 + ,520 + ,2342 + ,3 + ,1 + ,0 + ,27 + ,0 + ,17 + ,24069 + ,8891 + ,38798 + ,40 + ,39 + ,0 + ,20 + ,0 + ,7 + ,3738 + ,999 + ,3255 + ,5 + ,5 + ,0 + ,19 + ,0 + ,15 + ,18625 + ,7067 + ,24261 + ,38 + ,37 + ,0 + ,37 + ,0 + ,14 + ,36341 + ,4639 + ,18511 + ,32 + ,32 + ,0 + ,26 + ,0 + ,14 + ,24548 + ,5654 + ,40798 + ,41 + ,38 + ,0 + ,42 + ,0 + ,18 + ,21792 + ,6928 + ,28893 + ,46 + ,47 + ,0 + ,49 + ,0 + ,12 + ,26263 + ,1514 + ,21425 + ,47 + ,47 + ,0 + ,30 + ,0 + ,16 + ,23686 + ,9238 + ,50276 + ,37 + ,37 + ,0 + ,49 + ,0 + ,21 + ,49303 + ,8204 + ,37643 + ,51 + ,51 + ,0 + ,67 + ,1 + ,19 + ,25659 + ,5926 + ,30377 + ,49 + ,45 + ,0 + ,28 + ,0 + ,16 + ,28904 + ,5785 + ,27126 + ,21 + ,21 + ,0 + ,19 + ,0 + ,1 + ,2781 + ,4 + ,13 + ,1 + ,1 + ,0 + ,49 + ,1 + ,16 + ,29236 + ,5930 + ,42097 + ,44 + ,42 + ,0 + ,27 + ,0 + ,10 + ,19546 + ,3710 + ,24451 + ,26 + ,26 + ,0 + ,30 + ,6 + ,19 + ,22818 + ,705 + ,14335 + ,21 + ,21 + ,0 + ,22 + ,3 + ,12 + ,32689 + ,443 + ,5084 + ,4 + ,4 + ,0 + ,12 + ,1 + ,2 + ,5752 + ,2416 + ,9927 + ,10 + ,10 + ,0 + ,31 + ,2 + ,14 + ,22197 + ,7747 + ,43527 + ,43 + ,43 + ,0 + ,20 + ,0 + ,17 + ,20055 + ,5432 + ,27184 + ,34 + ,34 + ,0 + ,20 + ,0 + ,19 + ,25272 + ,4913 + ,21610 + ,32 + ,31 + ,0 + ,39 + ,0 + ,14 + ,82206 + ,2650 + ,20484 + ,20 + ,19 + ,0 + ,29 + ,3 + ,11 + ,32073 + ,2370 + ,20156 + ,34 + ,34 + ,0 + ,16 + ,1 + ,4 + ,5444 + ,775 + ,6012 + ,6 + ,6 + ,0 + ,27 + ,0 + ,16 + ,20154 + ,5576 + ,18475 + ,12 + ,11 + ,0 + ,21 + ,0 + ,20 + ,36944 + ,1352 + ,12645 + ,24 + ,24 + ,0 + ,19 + ,1 + ,12 + ,8019 + ,3080 + ,11017 + ,16 + ,16 + ,0 + ,35 + ,0 + ,15 + ,30884 + ,10205 + ,37623 + ,72 + ,72 + ,0 + ,14 + ,0 + ,16 + ,19540 + ,6095 + ,35873 + ,27 + ,21 + ,0) + ,dim=c(9 + ,289) + ,dimnames=list(c('logins' + ,'shared_compendiums' + ,'compendiums_reviewed' + ,'totsize' + ,'totrevisions' + ,'totseconds' + ,'tothyperlinks' + ,'totblogs' + ,'course_cid') + ,1:289)) > y <- array(NA,dim=c(9,289),dimnames=list(c('logins','shared_compendiums','compendiums_reviewed','totsize','totrevisions','totseconds','tothyperlinks','totblogs','course_cid'),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 = '4' > 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 logins shared_compendiums compendiums_reviewed totrevisions 1 112285 56 3 30 24188 2 84786 56 4 28 18273 3 83123 54 12 38 14130 4 101193 89 2 30 32287 5 38361 40 1 22 8654 6 68504 25 3 26 9245 7 119182 92 0 25 33251 8 22807 18 0 18 1271 9 17140 63 0 11 5279 10 116174 44 5 26 27101 11 57635 33 0 25 16373 12 66198 84 0 38 19716 13 71701 88 7 44 17753 14 57793 55 7 30 9028 15 80444 60 3 40 18653 16 53855 66 9 34 8828 17 97668 154 0 47 29498 18 133824 53 4 30 27563 19 101481 119 3 31 18293 20 99645 41 0 23 22530 21 114789 61 7 36 15977 22 99052 58 0 36 35082 23 67654 75 1 30 16116 24 65553 33 5 25 15849 25 97500 40 7 39 16026 26 69112 92 0 34 26569 27 82753 100 0 31 24785 28 85323 112 5 31 17569 29 72654 73 0 33 23825 30 30727 40 0 25 7869 31 77873 45 0 33 14975 32 117478 60 3 35 37791 33 74007 62 4 42 9605 34 90183 75 1 43 27295 35 61542 31 4 30 2746 36 101494 77 2 33 34461 37 27570 34 0 13 8098 38 55813 46 0 32 4787 39 79215 99 0 36 24919 40 1423 17 0 0 603 41 55461 66 2 28 16329 42 31081 30 1 14 12558 43 22996 76 0 17 7784 44 83122 146 2 32 28522 45 70106 67 10 30 22265 46 60578 56 6 35 14459 47 39992 107 0 20 14526 48 79892 58 5 28 22240 49 49810 34 4 28 11802 50 71570 61 1 39 7623 51 100708 119 2 34 11912 52 33032 42 2 26 7935 53 82875 66 0 39 18220 54 139077 89 8 39 19199 55 71595 44 3 33 19918 56 72260 66 0 28 21884 57 5950 24 0 4 2694 58 115762 259 8 39 15808 59 32551 17 5 18 3597 60 31701 64 3 14 5296 61 80670 41 1 29 25239 62 143558 68 5 44 29801 63 117105 168 1 21 18450 64 23789 43 1 16 7132 65 120733 132 5 28 34861 66 105195 105 0 35 35940 67 73107 71 12 28 16688 68 132068 112 8 38 24683 69 149193 94 8 23 46230 70 46821 82 8 36 10387 71 87011 70 8 32 21436 72 95260 57 2 29 30546 73 55183 53 0 25 19746 74 106671 103 5 27 15977 75 73511 121 8 36 22583 76 92945 62 2 28 17274 77 78664 52 5 23 16469 78 70054 52 12 40 14251 79 22618 32 6 23 3007 80 74011 62 7 40 16851 81 83737 45 2 28 21113 82 69094 46 0 34 17401 83 93133 63 4 33 23958 84 95536 75 3 28 23567 85 225920 88 6 34 13065 86 62133 46 2 30 15358 87 61370 53 0 33 14587 88 43836 37 1 22 12770 89 106117 90 0 38 24021 90 38692 63 5 26 9648 91 84651 78 2 35 20537 92 56622 25 0 8 7905 93 15986 45 0 24 4527 94 95364 46 5 29 30495 95 26706 41 0 20 7117 96 89691 144 1 29 17719 97 67267 82 0 45 27056 98 126846 91 1 37 33473 99 41140 71 1 33 9758 100 102860 63 2 33 21115 101 51715 53 6 25 7236 102 55801 62 1 32 13790 103 111813 63 4 29 32902 104 120293 32 2 28 25131 105 138599 39 3 28 30910 106 161647 62 0 31 35947 107 115929 117 10 52 29848 108 24266 34 0 21 6943 109 162901 92 9 24 42705 110 109825 93 7 41 31808 111 129838 54 0 33 26675 112 37510 144 0 32 8435 113 43750 14 4 19 7409 114 40652 61 4 20 14993 115 87771 109 0 31 36867 116 85872 38 0 31 33835 117 89275 73 0 32 24164 118 44418 75 1 18 12607 119 192565 50 0 23 22609 120 35232 61 1 17 5892 121 40909 55 0 20 17014 122 13294 77 0 12 5394 123 32387 75 4 17 9178 124 140867 72 0 30 6440 125 120662 50 4 31 21916 126 21233 32 4 10 4011 127 44332 53 3 13 5818 128 61056 42 0 22 18647 129 101338 71 0 42 20556 130 1168 10 0 1 238 131 13497 35 5 9 70 132 65567 65 0 32 22392 133 25162 25 4 11 3913 134 32334 66 0 25 12237 135 40735 41 0 36 8388 136 91413 86 1 31 22120 137 855 16 0 0 338 138 97068 42 5 24 11727 139 44339 19 0 13 3704 140 14116 19 0 8 3988 141 10288 45 0 13 3030 142 65622 65 0 19 13520 143 16563 35 0 18 1421 144 76643 95 2 33 20923 145 110681 49 7 40 20237 146 29011 37 1 22 3219 147 92696 64 8 38 3769 148 94785 38 2 24 12252 149 8773 34 0 8 1888 150 83209 32 2 35 14497 151 93815 65 0 43 28864 152 86687 52 0 43 21721 153 34553 62 1 14 4821 154 105547 65 3 41 33644 155 103487 83 0 38 15923 156 213688 95 3 45 42935 157 71220 29 0 31 18864 158 23517 18 0 13 4977 159 56926 33 0 28 7785 160 91721 247 4 31 17939 161 115168 139 4 40 23436 162 111194 29 11 30 325 163 51009 118 0 16 13539 164 135777 110 0 37 34538 165 51513 67 4 30 12198 166 74163 42 0 35 26924 167 51633 65 1 32 12716 168 75345 94 0 27 8172 169 33416 64 0 20 10855 170 83305 81 0 18 11932 171 98952 95 9 31 14300 172 102372 67 1 31 25515 173 37238 63 3 21 2805 174 103772 83 10 39 29402 175 123969 45 5 41 16440 176 27142 30 0 13 11221 177 135400 70 2 32 28732 178 21399 32 0 18 5250 179 130115 83 1 39 28608 180 24874 31 2 14 8092 181 34988 67 4 7 4473 182 45549 66 0 17 1572 183 6023 10 0 0 2065 184 64466 70 2 30 14817 185 54990 103 1 37 16714 186 1644 5 0 0 556 187 6179 20 0 5 2089 188 3926 5 0 1 2658 189 32755 36 1 16 10695 190 34777 34 0 32 1669 191 73224 48 2 24 16267 192 27114 40 0 17 7768 193 20760 43 3 11 7252 194 37636 31 6 24 6387 195 65461 42 0 22 18715 196 30080 46 2 12 7936 197 24094 33 0 19 8643 198 69008 18 2 13 7294 199 54968 55 1 17 4570 200 46090 35 1 15 7185 201 27507 59 2 16 10058 202 10672 19 1 24 2342 203 34029 66 0 15 8509 204 46300 60 1 17 13275 205 24760 36 3 18 6816 206 18779 25 0 20 1930 207 21280 47 0 16 8086 208 40662 54 0 16 10737 209 28987 53 0 18 8033 210 22827 40 1 22 7058 211 18513 40 4 8 6782 212 30594 39 0 17 5401 213 24006 14 0 18 6521 214 27913 45 0 16 10856 215 42744 36 7 23 2154 216 12934 28 2 22 6117 217 22574 44 0 13 5238 218 41385 30 7 13 4820 219 18653 22 3 16 5615 220 18472 17 0 16 4272 221 30976 31 0 20 8702 222 63339 55 6 22 15340 223 25568 54 2 17 8030 224 33747 21 0 18 9526 225 4154 14 0 17 1278 226 19474 81 3 12 4236 227 35130 35 0 7 3023 228 39067 43 1 17 7196 229 13310 46 1 14 3394 230 65892 30 0 23 6371 231 4143 23 1 17 1574 232 28579 38 0 14 9620 233 51776 54 0 15 6978 234 21152 20 0 17 4911 235 38084 53 0 21 8645 236 27717 45 0 18 8987 237 32928 39 0 18 5544 238 11342 20 0 17 3083 239 19499 24 0 17 6909 240 16380 31 0 16 3189 241 36874 35 0 15 6745 242 48259 151 0 21 16724 243 16734 52 0 16 4850 244 28207 30 2 14 7025 245 30143 31 0 15 6047 246 41369 29 1 17 7377 247 45833 57 1 15 9078 248 29156 40 0 15 4605 249 35944 44 0 10 3238 250 36278 25 0 6 8100 251 45588 77 0 22 9653 252 45097 35 0 21 8914 253 3895 11 0 1 786 254 28394 63 1 18 6700 255 18632 44 0 17 5788 256 2325 19 0 4 593 257 25139 13 0 10 4506 258 27975 42 0 16 6382 259 14483 38 1 16 5621 260 13127 29 0 9 3997 261 5839 20 0 16 520 262 24069 27 0 17 8891 263 3738 20 0 7 999 264 18625 19 0 15 7067 265 36341 37 0 14 4639 266 24548 26 0 14 5654 267 21792 42 0 18 6928 268 26263 49 0 12 1514 269 23686 30 0 16 9238 270 49303 49 0 21 8204 271 25659 67 1 19 5926 272 28904 28 0 16 5785 273 2781 19 0 1 4 274 29236 49 1 16 5930 275 19546 27 0 10 3710 276 22818 30 6 19 705 277 32689 22 3 12 443 278 5752 12 1 2 2416 279 22197 31 2 14 7747 280 20055 20 0 17 5432 281 25272 20 0 19 4913 282 82206 39 0 14 2650 283 32073 29 3 11 2370 284 5444 16 1 4 775 285 20154 27 0 16 5576 286 36944 21 0 20 1352 287 8019 19 1 12 3080 288 30884 35 0 15 10205 289 19540 14 0 16 6095 totseconds tothyperlinks totblogs course_cid 1 146283 144 145 1 2 98364 103 101 1 3 86146 98 98 1 4 96933 135 132 1 5 79234 61 60 1 6 42551 39 38 1 7 195663 150 144 1 8 6853 5 5 1 9 21529 28 28 0 10 95757 84 84 1 11 85584 80 79 1 12 143983 130 127 1 13 75851 82 78 1 14 59238 60 60 1 15 93163 131 131 1 16 96037 84 84 1 17 151511 140 133 1 18 136368 151 150 1 19 112642 91 91 1 20 94728 138 132 1 21 105499 150 136 1 22 121527 124 124 1 23 127766 119 118 1 24 98958 73 70 1 25 77900 110 107 1 26 85646 123 119 1 27 98579 90 89 1 28 130767 116 112 1 29 131741 113 108 1 30 53907 56 52 1 31 178812 115 112 1 32 146761 119 116 1 33 82036 129 123 1 34 163253 127 125 1 35 27032 27 27 1 36 171975 175 162 1 37 65990 35 32 0 38 86572 64 64 1 39 159676 96 92 1 40 1929 0 0 0 41 85371 84 83 1 42 58391 41 41 1 43 31580 47 47 0 44 136815 126 120 1 45 120642 105 105 1 46 69107 80 79 1 47 50495 70 65 0 48 108016 73 70 1 49 46341 57 55 1 50 78348 40 39 1 51 79336 68 67 1 52 56968 21 21 1 53 93176 127 127 1 54 161632 154 152 1 55 87850 116 113 1 56 127969 102 99 1 57 15049 7 7 1 58 155135 148 141 1 59 25109 21 21 1 60 45824 35 35 1 61 102996 112 109 1 62 160604 137 133 1 63 158051 135 123 0 64 44547 26 26 0 65 162647 230 230 1 66 174141 181 166 1 67 60622 71 68 1 68 179566 147 147 1 69 184301 190 179 1 70 75661 64 61 1 71 96144 105 101 1 72 129847 107 108 1 73 117286 94 90 1 74 71180 116 114 1 75 109377 106 103 1 76 85298 143 142 1 77 73631 81 79 1 78 86767 89 88 1 79 23824 26 25 1 80 93487 84 83 1 81 82981 113 113 1 82 73815 120 118 1 83 94552 110 110 1 84 132190 134 129 1 85 128754 54 51 1 86 66363 96 93 1 87 67808 78 76 1 88 61724 51 49 1 89 131722 121 118 1 90 68580 38 38 1 91 106175 145 141 1 92 55792 59 58 1 93 25157 27 27 1 94 76669 91 91 1 95 57283 48 48 0 96 105805 68 63 1 97 129484 58 56 1 98 72413 150 144 1 99 87831 74 73 1 100 96971 181 168 1 101 71299 65 64 1 102 77494 97 97 1 103 120336 121 117 1 104 93913 99 100 1 105 136048 152 149 1 106 181248 188 187 1 107 146123 138 127 1 108 32036 40 37 0 109 186646 254 245 1 110 102255 87 87 1 111 168237 178 177 1 112 64219 51 49 1 113 19630 49 49 1 114 76825 73 73 1 115 115338 176 177 1 116 109427 94 94 1 117 118168 120 117 1 118 84845 66 60 0 119 153197 56 55 1 120 29877 39 39 0 121 63506 66 64 0 122 22445 27 26 0 123 47695 65 64 0 124 68370 58 58 1 125 146304 98 95 1 126 38233 25 25 0 127 42071 26 26 0 128 50517 77 76 0 129 103950 130 129 1 130 5841 11 11 1 131 2341 2 2 0 132 84396 101 101 1 133 24610 31 28 1 134 35753 36 36 0 135 55515 120 89 1 136 209056 195 193 1 137 6622 4 4 1 138 115814 89 84 1 139 11609 24 23 0 140 13155 39 39 1 141 18274 14 14 0 142 72875 78 78 0 143 10112 15 14 0 144 142775 106 101 1 145 68847 83 82 1 146 17659 24 24 0 147 20112 37 36 1 148 61023 77 75 1 149 13983 16 16 1 150 65176 56 55 1 151 132432 132 131 1 152 112494 144 131 1 153 45109 40 39 0 154 170875 153 144 1 155 180759 143 139 1 156 214921 220 211 1 157 100226 79 78 1 158 32043 50 50 0 159 54454 39 39 1 160 78876 95 90 1 161 170745 169 166 1 162 6940 12 12 1 163 49025 63 57 0 164 122037 134 133 1 165 53782 69 69 1 166 127748 119 119 1 167 86839 119 119 1 168 44830 75 65 1 169 77395 63 61 0 170 89324 55 49 0 171 103300 103 101 1 172 112283 197 196 1 173 10901 16 15 1 174 120691 140 136 1 175 58106 89 89 1 176 57140 40 40 0 177 122422 125 123 1 178 25899 21 21 0 179 139296 167 163 1 180 52678 32 29 0 181 23853 36 35 0 182 17306 13 13 0 183 7953 5 5 1 184 89455 96 96 1 185 147866 151 151 1 186 4245 6 6 1 187 21509 13 13 1 188 7670 3 3 1 189 66675 57 56 0 190 14336 23 23 1 191 53608 61 57 1 192 30059 21 14 0 193 29668 43 43 0 194 22097 20 20 0 195 96841 82 72 0 196 41907 90 87 0 197 27080 25 21 0 198 35885 60 56 0 199 41247 61 59 0 200 28313 85 82 0 201 36845 43 43 0 202 16548 25 25 0 203 36134 41 38 0 204 55764 26 25 0 205 28910 38 38 0 206 13339 12 12 0 207 25319 29 29 0 208 66956 49 47 0 209 47487 46 45 0 210 52785 41 40 0 211 44683 31 30 0 212 35619 41 41 0 213 21920 26 25 0 214 45608 23 23 0 215 7721 14 14 0 216 20634 16 16 0 217 29788 25 26 0 218 31931 21 21 0 219 37754 32 27 0 220 32505 9 9 0 221 40557 35 33 0 222 94238 42 42 0 223 44197 68 68 0 224 43228 32 32 0 225 4103 6 6 0 226 44144 68 67 0 227 32868 33 33 0 228 27640 84 77 0 229 14063 46 46 0 230 28990 30 30 0 231 4694 0 0 0 232 42648 36 36 0 233 64329 47 46 0 234 21928 20 18 0 235 25836 50 48 0 236 22779 30 29 0 237 40820 30 28 0 238 27530 34 34 0 239 32378 33 33 0 240 10824 34 34 0 241 39613 37 33 0 242 60865 83 80 0 243 19787 32 32 0 244 20107 30 30 0 245 36605 43 41 0 246 40961 41 41 0 247 48231 51 51 0 248 39725 19 18 0 249 21455 37 34 0 250 23430 33 31 0 251 62991 41 39 0 252 49363 54 54 0 253 9604 14 14 0 254 24552 25 24 0 255 31493 25 24 0 256 3439 8 8 0 257 19555 26 26 0 258 21228 20 19 0 259 23177 11 11 0 260 22094 14 14 0 261 2342 3 1 0 262 38798 40 39 0 263 3255 5 5 0 264 24261 38 37 0 265 18511 32 32 0 266 40798 41 38 0 267 28893 46 47 0 268 21425 47 47 0 269 50276 37 37 0 270 37643 51 51 0 271 30377 49 45 0 272 27126 21 21 0 273 13 1 1 0 274 42097 44 42 0 275 24451 26 26 0 276 14335 21 21 0 277 5084 4 4 0 278 9927 10 10 0 279 43527 43 43 0 280 27184 34 34 0 281 21610 32 31 0 282 20484 20 19 0 283 20156 34 34 0 284 6012 6 6 0 285 18475 12 11 0 286 12645 24 24 0 287 11017 16 16 0 288 37623 72 72 0 289 35873 27 21 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) logins shared_compendiums 528.2800 9.5830 2014.6175 compendiums_reviewed totrevisions totseconds 539.6194 1.1471 0.2528 tothyperlinks totblogs course_cid 72.9181 4.3617 5358.7324 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -42081 -10786 -3571 7324 137064 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 5.283e+02 3.045e+03 0.173 0.86241 logins 9.583e+00 4.593e+01 0.209 0.83487 shared_compendiums 2.015e+03 4.824e+02 4.177 3.96e-05 *** compendiums_reviewed 5.396e+02 1.919e+02 2.813 0.00526 ** totrevisions 1.147e+00 2.782e-01 4.123 4.92e-05 *** totseconds 2.528e-01 6.111e-02 4.136 4.68e-05 *** tothyperlinks 7.292e+01 4.051e+02 0.180 0.85730 totblogs 4.362e+00 4.171e+02 0.010 0.99166 course_cid 5.359e+03 3.454e+03 1.552 0.12190 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 20090 on 280 degrees of freedom Multiple R-squared: 0.7598, Adjusted R-squared: 0.7529 F-statistic: 110.7 on 8 and 280 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 0.307162738 6.143255e-01 6.928373e-01 [2,] 0.168059014 3.361180e-01 8.319410e-01 [3,] 0.083835179 1.676704e-01 9.161648e-01 [4,] 0.043414971 8.682994e-02 9.565850e-01 [5,] 0.022473571 4.494714e-02 9.775264e-01 [6,] 0.023393471 4.678694e-02 9.766065e-01 [7,] 0.037387459 7.477492e-02 9.626125e-01 [8,] 0.050988930 1.019779e-01 9.490111e-01 [9,] 0.067598301 1.351966e-01 9.324017e-01 [10,] 0.101840233 2.036805e-01 8.981598e-01 [11,] 0.071925249 1.438505e-01 9.280748e-01 [12,] 0.049293410 9.858682e-02 9.507066e-01 [13,] 0.041486756 8.297351e-02 9.585132e-01 [14,] 0.037423428 7.484686e-02 9.625766e-01 [15,] 0.056035286 1.120706e-01 9.439647e-01 [16,] 0.038238467 7.647693e-02 9.617615e-01 [17,] 0.025780286 5.156057e-02 9.742197e-01 [18,] 0.020041284 4.008257e-02 9.799587e-01 [19,] 0.014719414 2.943883e-02 9.852806e-01 [20,] 0.012013663 2.402733e-02 9.879863e-01 [21,] 0.007785767 1.557153e-02 9.922142e-01 [22,] 0.005370678 1.074136e-02 9.946293e-01 [23,] 0.003445684 6.891368e-03 9.965543e-01 [24,] 0.009151172 1.830234e-02 9.908488e-01 [25,] 0.013286346 2.657269e-02 9.867137e-01 [26,] 0.008937157 1.787431e-02 9.910628e-01 [27,] 0.008307689 1.661538e-02 9.916923e-01 [28,] 0.005828222 1.165644e-02 9.941718e-01 [29,] 0.004057219 8.114439e-03 9.959428e-01 [30,] 0.003798831 7.597662e-03 9.962012e-01 [31,] 0.005377918 1.075584e-02 9.946221e-01 [32,] 0.003538664 7.077328e-03 9.964613e-01 [33,] 0.002736755 5.473511e-03 9.972632e-01 [34,] 0.008298208 1.659642e-02 9.917018e-01 [35,] 0.006915767 1.383153e-02 9.930842e-01 [36,] 0.004756638 9.513277e-03 9.952434e-01 [37,] 0.003255439 6.510878e-03 9.967446e-01 [38,] 0.002411606 4.823212e-03 9.975884e-01 [39,] 0.004319437 8.638875e-03 9.956806e-01 [40,] 0.021544202 4.308840e-02 9.784558e-01 [41,] 0.018071360 3.614272e-02 9.819286e-01 [42,] 0.013181950 2.636390e-02 9.868180e-01 [43,] 0.022806789 4.561358e-02 9.771932e-01 [44,] 0.019582198 3.916440e-02 9.804178e-01 [45,] 0.015556283 3.111257e-02 9.844437e-01 [46,] 0.013203490 2.640698e-02 9.867965e-01 [47,] 0.010250464 2.050093e-02 9.897495e-01 [48,] 0.007550920 1.510184e-02 9.924491e-01 [49,] 0.005802150 1.160430e-02 9.941979e-01 [50,] 0.004195665 8.391331e-03 9.958043e-01 [51,] 0.006780204 1.356041e-02 9.932198e-01 [52,] 0.014691704 2.938341e-02 9.853083e-01 [53,] 0.011344779 2.268956e-02 9.886552e-01 [54,] 0.010333444 2.066689e-02 9.896666e-01 [55,] 0.009761510 1.952302e-02 9.902385e-01 [56,] 0.007870389 1.574078e-02 9.921296e-01 [57,] 0.006238688 1.247738e-02 9.937613e-01 [58,] 0.004707051 9.414102e-03 9.952929e-01 [59,] 0.006766971 1.353394e-02 9.932330e-01 [60,] 0.005193656 1.038731e-02 9.948063e-01 [61,] 0.003889173 7.778345e-03 9.961108e-01 [62,] 0.004024209 8.048417e-03 9.959758e-01 [63,] 0.006293945 1.258789e-02 9.937061e-01 [64,] 0.008771261 1.754252e-02 9.912287e-01 [65,] 0.007218922 1.443784e-02 9.927811e-01 [66,] 0.005826155 1.165231e-02 9.941738e-01 [67,] 0.006802754 1.360551e-02 9.931972e-01 [68,] 0.006349927 1.269985e-02 9.936501e-01 [69,] 0.005527758 1.105552e-02 9.944722e-01 [70,] 0.004173140 8.346281e-03 9.958269e-01 [71,] 0.003145073 6.290147e-03 9.968549e-01 [72,] 0.002369561 4.739122e-03 9.976304e-01 [73,] 0.001748532 3.497063e-03 9.982515e-01 [74,] 0.978200843 4.359831e-02 2.179916e-02 [75,] 0.973424623 5.315075e-02 2.657538e-02 [76,] 0.967597075 6.480585e-02 3.240292e-02 [77,] 0.962095257 7.580949e-02 3.790474e-02 [78,] 0.956352198 8.729560e-02 4.364780e-02 [79,] 0.958957792 8.208442e-02 4.104221e-02 [80,] 0.951166037 9.766793e-02 4.883396e-02 [81,] 0.948990328 1.020193e-01 5.100967e-02 [82,] 0.945218764 1.095625e-01 5.478124e-02 [83,] 0.937889090 1.242218e-01 6.211091e-02 [84,] 0.928587651 1.428247e-01 7.141235e-02 [85,] 0.920234463 1.595311e-01 7.976554e-02 [86,] 0.936153370 1.276933e-01 6.384663e-02 [87,] 0.952870858 9.425828e-02 4.712914e-02 [88,] 0.957187258 8.562548e-02 4.281274e-02 [89,] 0.951312852 9.737430e-02 4.868715e-02 [90,] 0.946227406 1.075452e-01 5.377259e-02 [91,] 0.940748635 1.185027e-01 5.925137e-02 [92,] 0.931064576 1.378708e-01 6.893542e-02 [93,] 0.952931418 9.413716e-02 4.706858e-02 [94,] 0.961263931 7.747214e-02 3.873607e-02 [95,] 0.973853941 5.229212e-02 2.614606e-02 [96,] 0.976115960 4.776808e-02 2.388404e-02 [97,] 0.971118131 5.776374e-02 2.888187e-02 [98,] 0.967491990 6.501602e-02 3.250801e-02 [99,] 0.962943469 7.411306e-02 3.705653e-02 [100,] 0.961433403 7.713319e-02 3.856660e-02 [101,] 0.963393057 7.321389e-02 3.660694e-02 [102,] 0.956060944 8.787811e-02 4.393906e-02 [103,] 0.964534338 7.093132e-02 3.546566e-02 [104,] 0.965324811 6.935038e-02 3.467519e-02 [105,] 0.961986192 7.602762e-02 3.801381e-02 [106,] 0.954783016 9.043397e-02 4.521698e-02 [107,] 0.947522131 1.049557e-01 5.247787e-02 [108,] 0.999928435 1.431310e-04 7.156549e-05 [109,] 0.999904582 1.908363e-04 9.541813e-05 [110,] 0.999875948 2.481044e-04 1.240522e-04 [111,] 0.999835319 3.293629e-04 1.646814e-04 [112,] 0.999804111 3.917781e-04 1.958890e-04 [113,] 0.999999911 1.777208e-07 8.886040e-08 [114,] 0.999999917 1.667079e-07 8.335397e-08 [115,] 0.999999877 2.460083e-07 1.230042e-07 [116,] 0.999999839 3.213822e-07 1.606911e-07 [117,] 0.999999775 4.504507e-07 2.252253e-07 [118,] 0.999999696 6.073353e-07 3.036676e-07 [119,] 0.999999561 8.787621e-07 4.393811e-07 [120,] 0.999999357 1.286556e-06 6.432782e-07 [121,] 0.999999272 1.455034e-06 7.275168e-07 [122,] 0.999999019 1.961487e-06 9.807434e-07 [123,] 0.999998713 2.574002e-06 1.287001e-06 [124,] 0.999999449 1.101180e-06 5.505899e-07 [125,] 0.999999590 8.205944e-07 4.102972e-07 [126,] 0.999999403 1.193574e-06 5.967870e-07 [127,] 0.999999374 1.252054e-06 6.260271e-07 [128,] 0.999999623 7.533485e-07 3.766743e-07 [129,] 0.999999467 1.066331e-06 5.331655e-07 [130,] 0.999999205 1.589921e-06 7.949605e-07 [131,] 0.999999167 1.665042e-06 8.325210e-07 [132,] 0.999998735 2.530984e-06 1.265492e-06 [133,] 0.999998812 2.376830e-06 1.188415e-06 [134,] 0.999998705 2.589800e-06 1.294900e-06 [135,] 0.999998062 3.876975e-06 1.938488e-06 [136,] 0.999998948 2.103321e-06 1.051660e-06 [137,] 0.999999583 8.330324e-07 4.165162e-07 [138,] 0.999999416 1.167809e-06 5.839044e-07 [139,] 0.999999281 1.437485e-06 7.187424e-07 [140,] 0.999999102 1.795980e-06 8.979899e-07 [141,] 0.999999276 1.447216e-06 7.236082e-07 [142,] 0.999998929 2.141327e-06 1.070663e-06 [143,] 0.999999573 8.548993e-07 4.274497e-07 [144,] 0.999999379 1.242287e-06 6.211437e-07 [145,] 0.999999974 5.213414e-08 2.606707e-08 [146,] 0.999999958 8.419119e-08 4.209560e-08 [147,] 0.999999932 1.365238e-07 6.826188e-08 [148,] 0.999999901 1.975431e-07 9.877155e-08 [149,] 0.999999875 2.490769e-07 1.245384e-07 [150,] 0.999999803 3.948862e-07 1.974431e-07 [151,] 0.999999998 3.981023e-09 1.990512e-09 [152,] 0.999999997 6.097076e-09 3.048538e-09 [153,] 0.999999998 3.970667e-09 1.985334e-09 [154,] 0.999999998 4.918144e-09 2.459072e-09 [155,] 0.999999998 4.681161e-09 2.340581e-09 [156,] 0.999999998 4.010240e-09 2.005120e-09 [157,] 0.999999998 4.594541e-09 2.297271e-09 [158,] 0.999999997 5.876421e-09 2.938210e-09 [159,] 0.999999999 1.244701e-09 6.223507e-10 [160,] 0.999999999 1.687825e-09 8.439124e-10 [161,] 0.999999999 2.861579e-09 1.430790e-09 [162,] 0.999999998 4.929332e-09 2.464666e-09 [163,] 0.999999999 2.749608e-09 1.374804e-09 [164,] 1.000000000 3.130401e-10 1.565201e-10 [165,] 1.000000000 4.964079e-10 2.482040e-10 [166,] 1.000000000 4.429616e-11 2.214808e-11 [167,] 1.000000000 8.235144e-11 4.117572e-11 [168,] 1.000000000 1.582494e-11 7.912472e-12 [169,] 1.000000000 2.408434e-11 1.204217e-11 [170,] 1.000000000 4.147902e-11 2.073951e-11 [171,] 1.000000000 1.774128e-11 8.870640e-12 [172,] 1.000000000 3.426284e-11 1.713142e-11 [173,] 1.000000000 5.840963e-11 2.920481e-11 [174,] 1.000000000 1.101317e-11 5.506584e-12 [175,] 1.000000000 2.025930e-11 1.012965e-11 [176,] 1.000000000 2.628185e-11 1.314092e-11 [177,] 1.000000000 3.360839e-11 1.680420e-11 [178,] 1.000000000 4.958218e-11 2.479109e-11 [179,] 1.000000000 6.742550e-11 3.371275e-11 [180,] 1.000000000 1.251071e-10 6.255356e-11 [181,] 1.000000000 2.372460e-10 1.186230e-10 [182,] 1.000000000 3.806129e-10 1.903064e-10 [183,] 1.000000000 6.937333e-10 3.468667e-10 [184,] 0.999999999 1.313558e-09 6.567789e-10 [185,] 0.999999999 1.754780e-09 8.773898e-10 [186,] 0.999999998 3.133475e-09 1.566738e-09 [187,] 1.000000000 3.122058e-10 1.561029e-10 [188,] 1.000000000 2.473805e-10 1.236903e-10 [189,] 1.000000000 2.872610e-10 1.436305e-10 [190,] 1.000000000 4.977831e-10 2.488915e-10 [191,] 1.000000000 5.251394e-10 2.625697e-10 [192,] 0.999999999 1.014451e-09 5.072257e-10 [193,] 0.999999999 1.629296e-09 8.146479e-10 [194,] 0.999999999 2.806081e-09 1.403041e-09 [195,] 0.999999997 5.217100e-09 2.608550e-09 [196,] 0.999999995 9.326431e-09 4.663215e-09 [197,] 0.999999991 1.760038e-08 8.800191e-09 [198,] 0.999999985 2.931322e-08 1.465661e-08 [199,] 0.999999987 2.576705e-08 1.288353e-08 [200,] 0.999999983 3.459434e-08 1.729717e-08 [201,] 0.999999968 6.447059e-08 3.223530e-08 [202,] 0.999999940 1.193767e-07 5.968833e-08 [203,] 0.999999895 2.098091e-07 1.049045e-07 [204,] 0.999999880 2.396521e-07 1.198260e-07 [205,] 0.999999866 2.683968e-07 1.341984e-07 [206,] 0.999999762 4.761564e-07 2.380782e-07 [207,] 0.999999726 5.488897e-07 2.744448e-07 [208,] 0.999999630 7.405597e-07 3.702799e-07 [209,] 0.999999422 1.155673e-06 5.778363e-07 [210,] 0.999999000 1.999631e-06 9.998155e-07 [211,] 0.999998467 3.065058e-06 1.532529e-06 [212,] 0.999997980 4.039400e-06 2.019700e-06 [213,] 0.999996456 7.087053e-06 3.543526e-06 [214,] 0.999995767 8.465410e-06 4.232705e-06 [215,] 0.999996677 6.645797e-06 3.322898e-06 [216,] 0.999995389 9.222882e-06 4.611441e-06 [217,] 0.999992260 1.547981e-05 7.739903e-06 [218,] 0.999990521 1.895881e-05 9.479404e-06 [219,] 0.999998999 2.001142e-06 1.000571e-06 [220,] 0.999998684 2.632614e-06 1.316307e-06 [221,] 0.999997553 4.894663e-06 2.447331e-06 [222,] 0.999996517 6.966669e-06 3.483334e-06 [223,] 0.999993652 1.269537e-05 6.347683e-06 [224,] 0.999989180 2.163951e-05 1.081975e-05 [225,] 0.999980912 3.817602e-05 1.908801e-05 [226,] 0.999966311 6.737810e-05 3.368905e-05 [227,] 0.999968837 6.232667e-05 3.116333e-05 [228,] 0.999952054 9.589244e-05 4.794622e-05 [229,] 0.999928695 1.426106e-04 7.130529e-05 [230,] 0.999890713 2.185732e-04 1.092866e-04 [231,] 0.999830922 3.381565e-04 1.690782e-04 [232,] 0.999808224 3.835522e-04 1.917761e-04 [233,] 0.999692817 6.143654e-04 3.071827e-04 [234,] 0.999482175 1.035651e-03 5.178255e-04 [235,] 0.999304606 1.390787e-03 6.953936e-04 [236,] 0.999013760 1.972481e-03 9.862404e-04 [237,] 0.998387188 3.225624e-03 1.612812e-03 [238,] 0.997811582 4.376837e-03 2.188418e-03 [239,] 0.999144049 1.711902e-03 8.559511e-04 [240,] 0.998568058 2.863885e-03 1.431942e-03 [241,] 0.997739714 4.520572e-03 2.260286e-03 [242,] 0.996358273 7.283454e-03 3.641727e-03 [243,] 0.994129414 1.174117e-02 5.870586e-03 [244,] 0.992856107 1.428779e-02 7.143893e-03 [245,] 0.989780084 2.043983e-02 1.021992e-02 [246,] 0.986473726 2.705255e-02 1.352627e-02 [247,] 0.979766522 4.046696e-02 2.023348e-02 [248,] 0.974534305 5.093139e-02 2.546569e-02 [249,] 0.965133332 6.973334e-02 3.486667e-02 [250,] 0.967235260 6.552948e-02 3.276474e-02 [251,] 0.950804342 9.839132e-02 4.919566e-02 [252,] 0.941933261 1.161335e-01 5.806674e-02 [253,] 0.915618749 1.687625e-01 8.438125e-02 [254,] 0.888690903 2.226182e-01 1.113091e-01 [255,] 0.844006963 3.119861e-01 1.559930e-01 [256,] 0.816470330 3.670593e-01 1.835297e-01 [257,] 0.810857298 3.782854e-01 1.891427e-01 [258,] 0.754453280 4.910934e-01 2.455467e-01 [259,] 0.682422743 6.351545e-01 3.175773e-01 [260,] 0.765388242 4.692235e-01 2.346118e-01 [261,] 0.676564394 6.468712e-01 3.234356e-01 [262,] 0.655566810 6.888664e-01 3.444332e-01 [263,] 0.914071325 1.718573e-01 8.592867e-02 [264,] 0.944295850 1.114083e-01 5.570415e-02 [265,] 0.890838547 2.183229e-01 1.091615e-01 [266,] 0.905851460 1.882971e-01 9.414854e-02 > postscript(file="/var/wessaorg/rcomp/tmp/1euqg1324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/2njje1324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/365ss1324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/4av0q1324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5y8n21324145294.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 7775.7772 1420.0847 -13518.6087 2278.8402 -16459.5432 17933.8491 7 8 9 10 11 12 -240.8192 3457.8234 -3588.8140 23979.8658 -8650.3435 -30042.2078 13 14 15 16 17 18 -18730.7895 -8877.7754 -8714.8850 -30035.6479 -17978.4392 25431.4449 19 20 21 22 23 24 15193.8497 20527.9996 18264.7705 -7359.1211 -17127.7298 -13035.1420 25 26 27 28 29 30 9521.1825 -17616.5121 -1118.7795 -10591.7092 -21079.2198 -15996.1021 31 32 33 34 35 36 -17501.3966 -3542.7005 -4892.9376 -24020.6381 19041.7336 -23433.2243 37 38 39 40 41 42 -8959.7482 -101.8572 -23392.6932 -447.4581 -16993.4841 -16995.3554 43 44 45 46 47 48 -7977.3759 -22471.0041 -36905.7478 -17051.3628 -7167.5393 -10175.0672 49 50 51 52 53 54 -9217.9742 10404.1422 32336.4386 -16441.1940 1044.6059 20405.2553 55 56 57 58 59 60 -12569.0572 -14686.7899 -9760.5008 1478.5744 -5380.7072 -8760.1904 61 62 63 64 65 66 -6900.4229 17854.2787 30126.6432 -9249.9531 -10475.1836 -19749.5174 67 68 69 70 71 72 -12684.3672 3424.0765 -372.0612 -31366.8962 -9918.9255 -6983.9013 73 74 75 76 77 78 -24245.0321 29879.8610 -30808.2930 14903.8073 6040.9775 -27243.5644 79 80 81 82 83 84 -19550.6279 -17603.8757 4354.7794 -3463.7447 894.6354 -3002.5614 85 86 87 88 89 90 137064.0312 -6209.3563 -2723.2346 -10474.1071 8675.7866 -23240.0762 91 92 93 94 95 96 -6482.1840 18453.4946 -16921.2513 1922.4606 -11359.7933 12458.5582 97 98 99 100 101 102 -31926.3117 29841.1763 -24357.4463 11870.4100 -11598.9544 -12864.5674 103 104 105 106 107 108 4123.7824 34740.7519 29607.5292 36866.5858 -21074.7208 -7059.8949 109 110 111 112 113 114 9296.9503 -2235.8083 18752.1998 -16865.1354 2170.4463 -26928.5223 115 116 117 118 119 120 -20936.8408 -10842.1287 -1426.3803 -9537.9737 104807.4306 5606.7126 121 122 123 124 125 126 -11598.8216 -8390.3988 -13694.2456 88950.5901 19829.2014 -9253.4243 127 128 129 130 131 132 10919.9404 8149.0512 12210.5116 -7953.9351 -3122.9209 -13033.5359 133 134 135 136 137 138 -8050.4768 -8173.1341 -17763.4138 -27316.8568 -7555.9711 18173.1399 139 140 141 142 143 144 27580.1378 -7183.5975 -6863.1242 14261.6272 645.3913 -20249.5024 145 146 147 148 149 150 21612.0802 4231.2019 37311.2998 36133.5522 -8693.3373 16672.8377 151 152 153 154 155 156 -12678.2859 -7323.5627 3842.4926 -22698.8930 1311.3232 56027.6916 157 158 159 160 161 162 -4745.7783 -1871.1133 9905.5222 10846.1812 -4782.4715 63625.3697 163 164 165 166 167 168 7951.5718 28054.6084 -12181.5878 -23383.3754 -19891.5878 27529.8044 169 170 171 172 173 174 -15392.0247 31798.3711 6830.0927 4231.6693 6166.3621 -19136.0587 175 176 177 178 179 180 45030.3670 -11094.2591 33992.2758 -3340.4312 19460.0563 -12592.3428 181 182 183 184 185 186 8044.1058 28032.5465 -4725.1808 -9335.6821 -42081.2029 -6465.3613 187 188 189 190 191 192 -11435.3391 -7768.0095 -12288.3663 3980.8218 12990.6559 -1071.7764 193 194 195 196 197 198 -11300.6262 -2685.2076 419.8854 -8031.6142 -5676.8469 35206.4092 199 200 201 202 203 204 22351.2908 13163.4787 -10423.2119 -13805.0176 2724.7412 2681.1442 205 206 207 208 209 210 -9932.7629 705.9302 -6248.6783 -2035.8463 -6530.2779 -16573.1385 211 212 213 214 215 216 -16239.0376 2151.4416 -1395.1613 -7438.5538 9852.8583 -17232.1273 217 218 219 220 221 222 -865.0432 4229.0820 -15198.6789 -4665.0070 -3570.9939 -6337.5012 223 224 225 226 227 228 -14317.9593 -1022.1312 -8648.7023 -15617.4573 16163.3042 5236.7946 229 230 231 232 233 234 -8231.0444 35710.9558 -10785.8056 -4464.9053 14743.8378 -1454.2300 235 236 237 238 239 240 5413.6810 -1336.1648 3325.9676 -11673.9977 -9092.1587 -2100.7122 241 242 243 244 245 246 7324.3709 -6017.6006 -5964.2079 348.4075 1720.2725 7390.7695 247 248 249 250 251 252 8104.1205 3362.8387 17614.3349 14517.3573 2295.9042 6025.9996 253 254 255 256 257 258 -1689.3718 -284.6652 -8018.6510 -2711.5422 6969.1586 4182.7887 259 260 261 262 263 264 -10214.0525 -3787.0978 -4926.4181 -8983.7503 -3114.3484 -7350.5978 265 266 267 268 269 270 15430.3350 -3737.3091 -7661.1836 8005.4396 -11927.6581 14106.4902 271 272 273 274 275 276 -6023.7573 4358.3133 1445.8700 -3244.7077 917.5384 -6393.1730 277 278 279 280 281 282 17328.2852 -4038.4484 -13423.7245 -5568.0219 732.8970 63990.7185 283 284 285 286 287 288 8846.3634 -2282.9797 -1255.8275 18820.3326 -8735.5727 -4853.7493 289 -7875.5501 > postscript(file="/var/wessaorg/rcomp/tmp/6sxlr1324145294.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 7775.7772 NA 1 1420.0847 7775.7772 2 -13518.6087 1420.0847 3 2278.8402 -13518.6087 4 -16459.5432 2278.8402 5 17933.8491 -16459.5432 6 -240.8192 17933.8491 7 3457.8234 -240.8192 8 -3588.8140 3457.8234 9 23979.8658 -3588.8140 10 -8650.3435 23979.8658 11 -30042.2078 -8650.3435 12 -18730.7895 -30042.2078 13 -8877.7754 -18730.7895 14 -8714.8850 -8877.7754 15 -30035.6479 -8714.8850 16 -17978.4392 -30035.6479 17 25431.4449 -17978.4392 18 15193.8497 25431.4449 19 20527.9996 15193.8497 20 18264.7705 20527.9996 21 -7359.1211 18264.7705 22 -17127.7298 -7359.1211 23 -13035.1420 -17127.7298 24 9521.1825 -13035.1420 25 -17616.5121 9521.1825 26 -1118.7795 -17616.5121 27 -10591.7092 -1118.7795 28 -21079.2198 -10591.7092 29 -15996.1021 -21079.2198 30 -17501.3966 -15996.1021 31 -3542.7005 -17501.3966 32 -4892.9376 -3542.7005 33 -24020.6381 -4892.9376 34 19041.7336 -24020.6381 35 -23433.2243 19041.7336 36 -8959.7482 -23433.2243 37 -101.8572 -8959.7482 38 -23392.6932 -101.8572 39 -447.4581 -23392.6932 40 -16993.4841 -447.4581 41 -16995.3554 -16993.4841 42 -7977.3759 -16995.3554 43 -22471.0041 -7977.3759 44 -36905.7478 -22471.0041 45 -17051.3628 -36905.7478 46 -7167.5393 -17051.3628 47 -10175.0672 -7167.5393 48 -9217.9742 -10175.0672 49 10404.1422 -9217.9742 50 32336.4386 10404.1422 51 -16441.1940 32336.4386 52 1044.6059 -16441.1940 53 20405.2553 1044.6059 54 -12569.0572 20405.2553 55 -14686.7899 -12569.0572 56 -9760.5008 -14686.7899 57 1478.5744 -9760.5008 58 -5380.7072 1478.5744 59 -8760.1904 -5380.7072 60 -6900.4229 -8760.1904 61 17854.2787 -6900.4229 62 30126.6432 17854.2787 63 -9249.9531 30126.6432 64 -10475.1836 -9249.9531 65 -19749.5174 -10475.1836 66 -12684.3672 -19749.5174 67 3424.0765 -12684.3672 68 -372.0612 3424.0765 69 -31366.8962 -372.0612 70 -9918.9255 -31366.8962 71 -6983.9013 -9918.9255 72 -24245.0321 -6983.9013 73 29879.8610 -24245.0321 74 -30808.2930 29879.8610 75 14903.8073 -30808.2930 76 6040.9775 14903.8073 77 -27243.5644 6040.9775 78 -19550.6279 -27243.5644 79 -17603.8757 -19550.6279 80 4354.7794 -17603.8757 81 -3463.7447 4354.7794 82 894.6354 -3463.7447 83 -3002.5614 894.6354 84 137064.0312 -3002.5614 85 -6209.3563 137064.0312 86 -2723.2346 -6209.3563 87 -10474.1071 -2723.2346 88 8675.7866 -10474.1071 89 -23240.0762 8675.7866 90 -6482.1840 -23240.0762 91 18453.4946 -6482.1840 92 -16921.2513 18453.4946 93 1922.4606 -16921.2513 94 -11359.7933 1922.4606 95 12458.5582 -11359.7933 96 -31926.3117 12458.5582 97 29841.1763 -31926.3117 98 -24357.4463 29841.1763 99 11870.4100 -24357.4463 100 -11598.9544 11870.4100 101 -12864.5674 -11598.9544 102 4123.7824 -12864.5674 103 34740.7519 4123.7824 104 29607.5292 34740.7519 105 36866.5858 29607.5292 106 -21074.7208 36866.5858 107 -7059.8949 -21074.7208 108 9296.9503 -7059.8949 109 -2235.8083 9296.9503 110 18752.1998 -2235.8083 111 -16865.1354 18752.1998 112 2170.4463 -16865.1354 113 -26928.5223 2170.4463 114 -20936.8408 -26928.5223 115 -10842.1287 -20936.8408 116 -1426.3803 -10842.1287 117 -9537.9737 -1426.3803 118 104807.4306 -9537.9737 119 5606.7126 104807.4306 120 -11598.8216 5606.7126 121 -8390.3988 -11598.8216 122 -13694.2456 -8390.3988 123 88950.5901 -13694.2456 124 19829.2014 88950.5901 125 -9253.4243 19829.2014 126 10919.9404 -9253.4243 127 8149.0512 10919.9404 128 12210.5116 8149.0512 129 -7953.9351 12210.5116 130 -3122.9209 -7953.9351 131 -13033.5359 -3122.9209 132 -8050.4768 -13033.5359 133 -8173.1341 -8050.4768 134 -17763.4138 -8173.1341 135 -27316.8568 -17763.4138 136 -7555.9711 -27316.8568 137 18173.1399 -7555.9711 138 27580.1378 18173.1399 139 -7183.5975 27580.1378 140 -6863.1242 -7183.5975 141 14261.6272 -6863.1242 142 645.3913 14261.6272 143 -20249.5024 645.3913 144 21612.0802 -20249.5024 145 4231.2019 21612.0802 146 37311.2998 4231.2019 147 36133.5522 37311.2998 148 -8693.3373 36133.5522 149 16672.8377 -8693.3373 150 -12678.2859 16672.8377 151 -7323.5627 -12678.2859 152 3842.4926 -7323.5627 153 -22698.8930 3842.4926 154 1311.3232 -22698.8930 155 56027.6916 1311.3232 156 -4745.7783 56027.6916 157 -1871.1133 -4745.7783 158 9905.5222 -1871.1133 159 10846.1812 9905.5222 160 -4782.4715 10846.1812 161 63625.3697 -4782.4715 162 7951.5718 63625.3697 163 28054.6084 7951.5718 164 -12181.5878 28054.6084 165 -23383.3754 -12181.5878 166 -19891.5878 -23383.3754 167 27529.8044 -19891.5878 168 -15392.0247 27529.8044 169 31798.3711 -15392.0247 170 6830.0927 31798.3711 171 4231.6693 6830.0927 172 6166.3621 4231.6693 173 -19136.0587 6166.3621 174 45030.3670 -19136.0587 175 -11094.2591 45030.3670 176 33992.2758 -11094.2591 177 -3340.4312 33992.2758 178 19460.0563 -3340.4312 179 -12592.3428 19460.0563 180 8044.1058 -12592.3428 181 28032.5465 8044.1058 182 -4725.1808 28032.5465 183 -9335.6821 -4725.1808 184 -42081.2029 -9335.6821 185 -6465.3613 -42081.2029 186 -11435.3391 -6465.3613 187 -7768.0095 -11435.3391 188 -12288.3663 -7768.0095 189 3980.8218 -12288.3663 190 12990.6559 3980.8218 191 -1071.7764 12990.6559 192 -11300.6262 -1071.7764 193 -2685.2076 -11300.6262 194 419.8854 -2685.2076 195 -8031.6142 419.8854 196 -5676.8469 -8031.6142 197 35206.4092 -5676.8469 198 22351.2908 35206.4092 199 13163.4787 22351.2908 200 -10423.2119 13163.4787 201 -13805.0176 -10423.2119 202 2724.7412 -13805.0176 203 2681.1442 2724.7412 204 -9932.7629 2681.1442 205 705.9302 -9932.7629 206 -6248.6783 705.9302 207 -2035.8463 -6248.6783 208 -6530.2779 -2035.8463 209 -16573.1385 -6530.2779 210 -16239.0376 -16573.1385 211 2151.4416 -16239.0376 212 -1395.1613 2151.4416 213 -7438.5538 -1395.1613 214 9852.8583 -7438.5538 215 -17232.1273 9852.8583 216 -865.0432 -17232.1273 217 4229.0820 -865.0432 218 -15198.6789 4229.0820 219 -4665.0070 -15198.6789 220 -3570.9939 -4665.0070 221 -6337.5012 -3570.9939 222 -14317.9593 -6337.5012 223 -1022.1312 -14317.9593 224 -8648.7023 -1022.1312 225 -15617.4573 -8648.7023 226 16163.3042 -15617.4573 227 5236.7946 16163.3042 228 -8231.0444 5236.7946 229 35710.9558 -8231.0444 230 -10785.8056 35710.9558 231 -4464.9053 -10785.8056 232 14743.8378 -4464.9053 233 -1454.2300 14743.8378 234 5413.6810 -1454.2300 235 -1336.1648 5413.6810 236 3325.9676 -1336.1648 237 -11673.9977 3325.9676 238 -9092.1587 -11673.9977 239 -2100.7122 -9092.1587 240 7324.3709 -2100.7122 241 -6017.6006 7324.3709 242 -5964.2079 -6017.6006 243 348.4075 -5964.2079 244 1720.2725 348.4075 245 7390.7695 1720.2725 246 8104.1205 7390.7695 247 3362.8387 8104.1205 248 17614.3349 3362.8387 249 14517.3573 17614.3349 250 2295.9042 14517.3573 251 6025.9996 2295.9042 252 -1689.3718 6025.9996 253 -284.6652 -1689.3718 254 -8018.6510 -284.6652 255 -2711.5422 -8018.6510 256 6969.1586 -2711.5422 257 4182.7887 6969.1586 258 -10214.0525 4182.7887 259 -3787.0978 -10214.0525 260 -4926.4181 -3787.0978 261 -8983.7503 -4926.4181 262 -3114.3484 -8983.7503 263 -7350.5978 -3114.3484 264 15430.3350 -7350.5978 265 -3737.3091 15430.3350 266 -7661.1836 -3737.3091 267 8005.4396 -7661.1836 268 -11927.6581 8005.4396 269 14106.4902 -11927.6581 270 -6023.7573 14106.4902 271 4358.3133 -6023.7573 272 1445.8700 4358.3133 273 -3244.7077 1445.8700 274 917.5384 -3244.7077 275 -6393.1730 917.5384 276 17328.2852 -6393.1730 277 -4038.4484 17328.2852 278 -13423.7245 -4038.4484 279 -5568.0219 -13423.7245 280 732.8970 -5568.0219 281 63990.7185 732.8970 282 8846.3634 63990.7185 283 -2282.9797 8846.3634 284 -1255.8275 -2282.9797 285 18820.3326 -1255.8275 286 -8735.5727 18820.3326 287 -4853.7493 -8735.5727 288 -7875.5501 -4853.7493 289 NA -7875.5501 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 1420.0847 7775.7772 [2,] -13518.6087 1420.0847 [3,] 2278.8402 -13518.6087 [4,] -16459.5432 2278.8402 [5,] 17933.8491 -16459.5432 [6,] -240.8192 17933.8491 [7,] 3457.8234 -240.8192 [8,] -3588.8140 3457.8234 [9,] 23979.8658 -3588.8140 [10,] -8650.3435 23979.8658 [11,] -30042.2078 -8650.3435 [12,] -18730.7895 -30042.2078 [13,] -8877.7754 -18730.7895 [14,] -8714.8850 -8877.7754 [15,] -30035.6479 -8714.8850 [16,] -17978.4392 -30035.6479 [17,] 25431.4449 -17978.4392 [18,] 15193.8497 25431.4449 [19,] 20527.9996 15193.8497 [20,] 18264.7705 20527.9996 [21,] -7359.1211 18264.7705 [22,] -17127.7298 -7359.1211 [23,] -13035.1420 -17127.7298 [24,] 9521.1825 -13035.1420 [25,] -17616.5121 9521.1825 [26,] -1118.7795 -17616.5121 [27,] -10591.7092 -1118.7795 [28,] -21079.2198 -10591.7092 [29,] -15996.1021 -21079.2198 [30,] -17501.3966 -15996.1021 [31,] -3542.7005 -17501.3966 [32,] -4892.9376 -3542.7005 [33,] -24020.6381 -4892.9376 [34,] 19041.7336 -24020.6381 [35,] -23433.2243 19041.7336 [36,] -8959.7482 -23433.2243 [37,] -101.8572 -8959.7482 [38,] -23392.6932 -101.8572 [39,] -447.4581 -23392.6932 [40,] -16993.4841 -447.4581 [41,] -16995.3554 -16993.4841 [42,] -7977.3759 -16995.3554 [43,] -22471.0041 -7977.3759 [44,] -36905.7478 -22471.0041 [45,] -17051.3628 -36905.7478 [46,] -7167.5393 -17051.3628 [47,] -10175.0672 -7167.5393 [48,] -9217.9742 -10175.0672 [49,] 10404.1422 -9217.9742 [50,] 32336.4386 10404.1422 [51,] -16441.1940 32336.4386 [52,] 1044.6059 -16441.1940 [53,] 20405.2553 1044.6059 [54,] -12569.0572 20405.2553 [55,] -14686.7899 -12569.0572 [56,] -9760.5008 -14686.7899 [57,] 1478.5744 -9760.5008 [58,] -5380.7072 1478.5744 [59,] -8760.1904 -5380.7072 [60,] -6900.4229 -8760.1904 [61,] 17854.2787 -6900.4229 [62,] 30126.6432 17854.2787 [63,] -9249.9531 30126.6432 [64,] -10475.1836 -9249.9531 [65,] -19749.5174 -10475.1836 [66,] -12684.3672 -19749.5174 [67,] 3424.0765 -12684.3672 [68,] -372.0612 3424.0765 [69,] -31366.8962 -372.0612 [70,] -9918.9255 -31366.8962 [71,] -6983.9013 -9918.9255 [72,] -24245.0321 -6983.9013 [73,] 29879.8610 -24245.0321 [74,] -30808.2930 29879.8610 [75,] 14903.8073 -30808.2930 [76,] 6040.9775 14903.8073 [77,] -27243.5644 6040.9775 [78,] -19550.6279 -27243.5644 [79,] -17603.8757 -19550.6279 [80,] 4354.7794 -17603.8757 [81,] -3463.7447 4354.7794 [82,] 894.6354 -3463.7447 [83,] -3002.5614 894.6354 [84,] 137064.0312 -3002.5614 [85,] -6209.3563 137064.0312 [86,] -2723.2346 -6209.3563 [87,] -10474.1071 -2723.2346 [88,] 8675.7866 -10474.1071 [89,] -23240.0762 8675.7866 [90,] -6482.1840 -23240.0762 [91,] 18453.4946 -6482.1840 [92,] -16921.2513 18453.4946 [93,] 1922.4606 -16921.2513 [94,] -11359.7933 1922.4606 [95,] 12458.5582 -11359.7933 [96,] -31926.3117 12458.5582 [97,] 29841.1763 -31926.3117 [98,] -24357.4463 29841.1763 [99,] 11870.4100 -24357.4463 [100,] -11598.9544 11870.4100 [101,] -12864.5674 -11598.9544 [102,] 4123.7824 -12864.5674 [103,] 34740.7519 4123.7824 [104,] 29607.5292 34740.7519 [105,] 36866.5858 29607.5292 [106,] -21074.7208 36866.5858 [107,] -7059.8949 -21074.7208 [108,] 9296.9503 -7059.8949 [109,] -2235.8083 9296.9503 [110,] 18752.1998 -2235.8083 [111,] -16865.1354 18752.1998 [112,] 2170.4463 -16865.1354 [113,] -26928.5223 2170.4463 [114,] -20936.8408 -26928.5223 [115,] -10842.1287 -20936.8408 [116,] -1426.3803 -10842.1287 [117,] -9537.9737 -1426.3803 [118,] 104807.4306 -9537.9737 [119,] 5606.7126 104807.4306 [120,] -11598.8216 5606.7126 [121,] -8390.3988 -11598.8216 [122,] -13694.2456 -8390.3988 [123,] 88950.5901 -13694.2456 [124,] 19829.2014 88950.5901 [125,] -9253.4243 19829.2014 [126,] 10919.9404 -9253.4243 [127,] 8149.0512 10919.9404 [128,] 12210.5116 8149.0512 [129,] -7953.9351 12210.5116 [130,] -3122.9209 -7953.9351 [131,] -13033.5359 -3122.9209 [132,] -8050.4768 -13033.5359 [133,] -8173.1341 -8050.4768 [134,] -17763.4138 -8173.1341 [135,] -27316.8568 -17763.4138 [136,] -7555.9711 -27316.8568 [137,] 18173.1399 -7555.9711 [138,] 27580.1378 18173.1399 [139,] -7183.5975 27580.1378 [140,] -6863.1242 -7183.5975 [141,] 14261.6272 -6863.1242 [142,] 645.3913 14261.6272 [143,] -20249.5024 645.3913 [144,] 21612.0802 -20249.5024 [145,] 4231.2019 21612.0802 [146,] 37311.2998 4231.2019 [147,] 36133.5522 37311.2998 [148,] -8693.3373 36133.5522 [149,] 16672.8377 -8693.3373 [150,] -12678.2859 16672.8377 [151,] -7323.5627 -12678.2859 [152,] 3842.4926 -7323.5627 [153,] -22698.8930 3842.4926 [154,] 1311.3232 -22698.8930 [155,] 56027.6916 1311.3232 [156,] -4745.7783 56027.6916 [157,] -1871.1133 -4745.7783 [158,] 9905.5222 -1871.1133 [159,] 10846.1812 9905.5222 [160,] -4782.4715 10846.1812 [161,] 63625.3697 -4782.4715 [162,] 7951.5718 63625.3697 [163,] 28054.6084 7951.5718 [164,] -12181.5878 28054.6084 [165,] -23383.3754 -12181.5878 [166,] -19891.5878 -23383.3754 [167,] 27529.8044 -19891.5878 [168,] -15392.0247 27529.8044 [169,] 31798.3711 -15392.0247 [170,] 6830.0927 31798.3711 [171,] 4231.6693 6830.0927 [172,] 6166.3621 4231.6693 [173,] -19136.0587 6166.3621 [174,] 45030.3670 -19136.0587 [175,] -11094.2591 45030.3670 [176,] 33992.2758 -11094.2591 [177,] -3340.4312 33992.2758 [178,] 19460.0563 -3340.4312 [179,] -12592.3428 19460.0563 [180,] 8044.1058 -12592.3428 [181,] 28032.5465 8044.1058 [182,] -4725.1808 28032.5465 [183,] -9335.6821 -4725.1808 [184,] -42081.2029 -9335.6821 [185,] -6465.3613 -42081.2029 [186,] -11435.3391 -6465.3613 [187,] -7768.0095 -11435.3391 [188,] -12288.3663 -7768.0095 [189,] 3980.8218 -12288.3663 [190,] 12990.6559 3980.8218 [191,] -1071.7764 12990.6559 [192,] -11300.6262 -1071.7764 [193,] -2685.2076 -11300.6262 [194,] 419.8854 -2685.2076 [195,] -8031.6142 419.8854 [196,] -5676.8469 -8031.6142 [197,] 35206.4092 -5676.8469 [198,] 22351.2908 35206.4092 [199,] 13163.4787 22351.2908 [200,] -10423.2119 13163.4787 [201,] -13805.0176 -10423.2119 [202,] 2724.7412 -13805.0176 [203,] 2681.1442 2724.7412 [204,] -9932.7629 2681.1442 [205,] 705.9302 -9932.7629 [206,] -6248.6783 705.9302 [207,] -2035.8463 -6248.6783 [208,] -6530.2779 -2035.8463 [209,] -16573.1385 -6530.2779 [210,] -16239.0376 -16573.1385 [211,] 2151.4416 -16239.0376 [212,] -1395.1613 2151.4416 [213,] -7438.5538 -1395.1613 [214,] 9852.8583 -7438.5538 [215,] -17232.1273 9852.8583 [216,] -865.0432 -17232.1273 [217,] 4229.0820 -865.0432 [218,] -15198.6789 4229.0820 [219,] -4665.0070 -15198.6789 [220,] -3570.9939 -4665.0070 [221,] -6337.5012 -3570.9939 [222,] -14317.9593 -6337.5012 [223,] -1022.1312 -14317.9593 [224,] -8648.7023 -1022.1312 [225,] -15617.4573 -8648.7023 [226,] 16163.3042 -15617.4573 [227,] 5236.7946 16163.3042 [228,] -8231.0444 5236.7946 [229,] 35710.9558 -8231.0444 [230,] -10785.8056 35710.9558 [231,] -4464.9053 -10785.8056 [232,] 14743.8378 -4464.9053 [233,] -1454.2300 14743.8378 [234,] 5413.6810 -1454.2300 [235,] -1336.1648 5413.6810 [236,] 3325.9676 -1336.1648 [237,] -11673.9977 3325.9676 [238,] -9092.1587 -11673.9977 [239,] -2100.7122 -9092.1587 [240,] 7324.3709 -2100.7122 [241,] -6017.6006 7324.3709 [242,] -5964.2079 -6017.6006 [243,] 348.4075 -5964.2079 [244,] 1720.2725 348.4075 [245,] 7390.7695 1720.2725 [246,] 8104.1205 7390.7695 [247,] 3362.8387 8104.1205 [248,] 17614.3349 3362.8387 [249,] 14517.3573 17614.3349 [250,] 2295.9042 14517.3573 [251,] 6025.9996 2295.9042 [252,] -1689.3718 6025.9996 [253,] -284.6652 -1689.3718 [254,] -8018.6510 -284.6652 [255,] -2711.5422 -8018.6510 [256,] 6969.1586 -2711.5422 [257,] 4182.7887 6969.1586 [258,] -10214.0525 4182.7887 [259,] -3787.0978 -10214.0525 [260,] -4926.4181 -3787.0978 [261,] -8983.7503 -4926.4181 [262,] -3114.3484 -8983.7503 [263,] -7350.5978 -3114.3484 [264,] 15430.3350 -7350.5978 [265,] -3737.3091 15430.3350 [266,] -7661.1836 -3737.3091 [267,] 8005.4396 -7661.1836 [268,] -11927.6581 8005.4396 [269,] 14106.4902 -11927.6581 [270,] -6023.7573 14106.4902 [271,] 4358.3133 -6023.7573 [272,] 1445.8700 4358.3133 [273,] -3244.7077 1445.8700 [274,] 917.5384 -3244.7077 [275,] -6393.1730 917.5384 [276,] 17328.2852 -6393.1730 [277,] -4038.4484 17328.2852 [278,] -13423.7245 -4038.4484 [279,] -5568.0219 -13423.7245 [280,] 732.8970 -5568.0219 [281,] 63990.7185 732.8970 [282,] 8846.3634 63990.7185 [283,] -2282.9797 8846.3634 [284,] -1255.8275 -2282.9797 [285,] 18820.3326 -1255.8275 [286,] -8735.5727 18820.3326 [287,] -4853.7493 -8735.5727 [288,] -7875.5501 -4853.7493 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 1420.0847 7775.7772 2 -13518.6087 1420.0847 3 2278.8402 -13518.6087 4 -16459.5432 2278.8402 5 17933.8491 -16459.5432 6 -240.8192 17933.8491 7 3457.8234 -240.8192 8 -3588.8140 3457.8234 9 23979.8658 -3588.8140 10 -8650.3435 23979.8658 11 -30042.2078 -8650.3435 12 -18730.7895 -30042.2078 13 -8877.7754 -18730.7895 14 -8714.8850 -8877.7754 15 -30035.6479 -8714.8850 16 -17978.4392 -30035.6479 17 25431.4449 -17978.4392 18 15193.8497 25431.4449 19 20527.9996 15193.8497 20 18264.7705 20527.9996 21 -7359.1211 18264.7705 22 -17127.7298 -7359.1211 23 -13035.1420 -17127.7298 24 9521.1825 -13035.1420 25 -17616.5121 9521.1825 26 -1118.7795 -17616.5121 27 -10591.7092 -1118.7795 28 -21079.2198 -10591.7092 29 -15996.1021 -21079.2198 30 -17501.3966 -15996.1021 31 -3542.7005 -17501.3966 32 -4892.9376 -3542.7005 33 -24020.6381 -4892.9376 34 19041.7336 -24020.6381 35 -23433.2243 19041.7336 36 -8959.7482 -23433.2243 37 -101.8572 -8959.7482 38 -23392.6932 -101.8572 39 -447.4581 -23392.6932 40 -16993.4841 -447.4581 41 -16995.3554 -16993.4841 42 -7977.3759 -16995.3554 43 -22471.0041 -7977.3759 44 -36905.7478 -22471.0041 45 -17051.3628 -36905.7478 46 -7167.5393 -17051.3628 47 -10175.0672 -7167.5393 48 -9217.9742 -10175.0672 49 10404.1422 -9217.9742 50 32336.4386 10404.1422 51 -16441.1940 32336.4386 52 1044.6059 -16441.1940 53 20405.2553 1044.6059 54 -12569.0572 20405.2553 55 -14686.7899 -12569.0572 56 -9760.5008 -14686.7899 57 1478.5744 -9760.5008 58 -5380.7072 1478.5744 59 -8760.1904 -5380.7072 60 -6900.4229 -8760.1904 61 17854.2787 -6900.4229 62 30126.6432 17854.2787 63 -9249.9531 30126.6432 64 -10475.1836 -9249.9531 65 -19749.5174 -10475.1836 66 -12684.3672 -19749.5174 67 3424.0765 -12684.3672 68 -372.0612 3424.0765 69 -31366.8962 -372.0612 70 -9918.9255 -31366.8962 71 -6983.9013 -9918.9255 72 -24245.0321 -6983.9013 73 29879.8610 -24245.0321 74 -30808.2930 29879.8610 75 14903.8073 -30808.2930 76 6040.9775 14903.8073 77 -27243.5644 6040.9775 78 -19550.6279 -27243.5644 79 -17603.8757 -19550.6279 80 4354.7794 -17603.8757 81 -3463.7447 4354.7794 82 894.6354 -3463.7447 83 -3002.5614 894.6354 84 137064.0312 -3002.5614 85 -6209.3563 137064.0312 86 -2723.2346 -6209.3563 87 -10474.1071 -2723.2346 88 8675.7866 -10474.1071 89 -23240.0762 8675.7866 90 -6482.1840 -23240.0762 91 18453.4946 -6482.1840 92 -16921.2513 18453.4946 93 1922.4606 -16921.2513 94 -11359.7933 1922.4606 95 12458.5582 -11359.7933 96 -31926.3117 12458.5582 97 29841.1763 -31926.3117 98 -24357.4463 29841.1763 99 11870.4100 -24357.4463 100 -11598.9544 11870.4100 101 -12864.5674 -11598.9544 102 4123.7824 -12864.5674 103 34740.7519 4123.7824 104 29607.5292 34740.7519 105 36866.5858 29607.5292 106 -21074.7208 36866.5858 107 -7059.8949 -21074.7208 108 9296.9503 -7059.8949 109 -2235.8083 9296.9503 110 18752.1998 -2235.8083 111 -16865.1354 18752.1998 112 2170.4463 -16865.1354 113 -26928.5223 2170.4463 114 -20936.8408 -26928.5223 115 -10842.1287 -20936.8408 116 -1426.3803 -10842.1287 117 -9537.9737 -1426.3803 118 104807.4306 -9537.9737 119 5606.7126 104807.4306 120 -11598.8216 5606.7126 121 -8390.3988 -11598.8216 122 -13694.2456 -8390.3988 123 88950.5901 -13694.2456 124 19829.2014 88950.5901 125 -9253.4243 19829.2014 126 10919.9404 -9253.4243 127 8149.0512 10919.9404 128 12210.5116 8149.0512 129 -7953.9351 12210.5116 130 -3122.9209 -7953.9351 131 -13033.5359 -3122.9209 132 -8050.4768 -13033.5359 133 -8173.1341 -8050.4768 134 -17763.4138 -8173.1341 135 -27316.8568 -17763.4138 136 -7555.9711 -27316.8568 137 18173.1399 -7555.9711 138 27580.1378 18173.1399 139 -7183.5975 27580.1378 140 -6863.1242 -7183.5975 141 14261.6272 -6863.1242 142 645.3913 14261.6272 143 -20249.5024 645.3913 144 21612.0802 -20249.5024 145 4231.2019 21612.0802 146 37311.2998 4231.2019 147 36133.5522 37311.2998 148 -8693.3373 36133.5522 149 16672.8377 -8693.3373 150 -12678.2859 16672.8377 151 -7323.5627 -12678.2859 152 3842.4926 -7323.5627 153 -22698.8930 3842.4926 154 1311.3232 -22698.8930 155 56027.6916 1311.3232 156 -4745.7783 56027.6916 157 -1871.1133 -4745.7783 158 9905.5222 -1871.1133 159 10846.1812 9905.5222 160 -4782.4715 10846.1812 161 63625.3697 -4782.4715 162 7951.5718 63625.3697 163 28054.6084 7951.5718 164 -12181.5878 28054.6084 165 -23383.3754 -12181.5878 166 -19891.5878 -23383.3754 167 27529.8044 -19891.5878 168 -15392.0247 27529.8044 169 31798.3711 -15392.0247 170 6830.0927 31798.3711 171 4231.6693 6830.0927 172 6166.3621 4231.6693 173 -19136.0587 6166.3621 174 45030.3670 -19136.0587 175 -11094.2591 45030.3670 176 33992.2758 -11094.2591 177 -3340.4312 33992.2758 178 19460.0563 -3340.4312 179 -12592.3428 19460.0563 180 8044.1058 -12592.3428 181 28032.5465 8044.1058 182 -4725.1808 28032.5465 183 -9335.6821 -4725.1808 184 -42081.2029 -9335.6821 185 -6465.3613 -42081.2029 186 -11435.3391 -6465.3613 187 -7768.0095 -11435.3391 188 -12288.3663 -7768.0095 189 3980.8218 -12288.3663 190 12990.6559 3980.8218 191 -1071.7764 12990.6559 192 -11300.6262 -1071.7764 193 -2685.2076 -11300.6262 194 419.8854 -2685.2076 195 -8031.6142 419.8854 196 -5676.8469 -8031.6142 197 35206.4092 -5676.8469 198 22351.2908 35206.4092 199 13163.4787 22351.2908 200 -10423.2119 13163.4787 201 -13805.0176 -10423.2119 202 2724.7412 -13805.0176 203 2681.1442 2724.7412 204 -9932.7629 2681.1442 205 705.9302 -9932.7629 206 -6248.6783 705.9302 207 -2035.8463 -6248.6783 208 -6530.2779 -2035.8463 209 -16573.1385 -6530.2779 210 -16239.0376 -16573.1385 211 2151.4416 -16239.0376 212 -1395.1613 2151.4416 213 -7438.5538 -1395.1613 214 9852.8583 -7438.5538 215 -17232.1273 9852.8583 216 -865.0432 -17232.1273 217 4229.0820 -865.0432 218 -15198.6789 4229.0820 219 -4665.0070 -15198.6789 220 -3570.9939 -4665.0070 221 -6337.5012 -3570.9939 222 -14317.9593 -6337.5012 223 -1022.1312 -14317.9593 224 -8648.7023 -1022.1312 225 -15617.4573 -8648.7023 226 16163.3042 -15617.4573 227 5236.7946 16163.3042 228 -8231.0444 5236.7946 229 35710.9558 -8231.0444 230 -10785.8056 35710.9558 231 -4464.9053 -10785.8056 232 14743.8378 -4464.9053 233 -1454.2300 14743.8378 234 5413.6810 -1454.2300 235 -1336.1648 5413.6810 236 3325.9676 -1336.1648 237 -11673.9977 3325.9676 238 -9092.1587 -11673.9977 239 -2100.7122 -9092.1587 240 7324.3709 -2100.7122 241 -6017.6006 7324.3709 242 -5964.2079 -6017.6006 243 348.4075 -5964.2079 244 1720.2725 348.4075 245 7390.7695 1720.2725 246 8104.1205 7390.7695 247 3362.8387 8104.1205 248 17614.3349 3362.8387 249 14517.3573 17614.3349 250 2295.9042 14517.3573 251 6025.9996 2295.9042 252 -1689.3718 6025.9996 253 -284.6652 -1689.3718 254 -8018.6510 -284.6652 255 -2711.5422 -8018.6510 256 6969.1586 -2711.5422 257 4182.7887 6969.1586 258 -10214.0525 4182.7887 259 -3787.0978 -10214.0525 260 -4926.4181 -3787.0978 261 -8983.7503 -4926.4181 262 -3114.3484 -8983.7503 263 -7350.5978 -3114.3484 264 15430.3350 -7350.5978 265 -3737.3091 15430.3350 266 -7661.1836 -3737.3091 267 8005.4396 -7661.1836 268 -11927.6581 8005.4396 269 14106.4902 -11927.6581 270 -6023.7573 14106.4902 271 4358.3133 -6023.7573 272 1445.8700 4358.3133 273 -3244.7077 1445.8700 274 917.5384 -3244.7077 275 -6393.1730 917.5384 276 17328.2852 -6393.1730 277 -4038.4484 17328.2852 278 -13423.7245 -4038.4484 279 -5568.0219 -13423.7245 280 732.8970 -5568.0219 281 63990.7185 732.8970 282 8846.3634 63990.7185 283 -2282.9797 8846.3634 284 -1255.8275 -2282.9797 285 18820.3326 -1255.8275 286 -8735.5727 18820.3326 287 -4853.7493 -8735.5727 288 -7875.5501 -4853.7493 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/72zji1324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/8nhn11324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/9beiz1324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/wessaorg/rcomp/tmp/10x1d31324145294.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/1166wh1324145294.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/121x801324145294.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/135v2y1324145294.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/14tnpv1324145294.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/15n9ea1324145294.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/wessaorg/rcomp/tmp/16q2c21324145294.tab") + } > > try(system("convert tmp/1euqg1324145294.ps tmp/1euqg1324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/2njje1324145294.ps tmp/2njje1324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/365ss1324145294.ps tmp/365ss1324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/4av0q1324145294.ps tmp/4av0q1324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/5y8n21324145294.ps tmp/5y8n21324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/6sxlr1324145294.ps tmp/6sxlr1324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/72zji1324145294.ps tmp/72zji1324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/8nhn11324145294.ps tmp/8nhn11324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/9beiz1324145294.ps tmp/9beiz1324145294.png",intern=TRUE)) character(0) > try(system("convert tmp/10x1d31324145294.ps tmp/10x1d31324145294.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 9.270 0.695 9.978