R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(210907 + ,30 + ,79 + ,112285 + ,120982 + ,28 + ,58 + ,84786 + ,176508 + ,38 + ,60 + ,83123 + ,179321 + ,30 + ,108 + ,101193 + ,123185 + ,22 + ,49 + ,38361 + ,52746 + ,26 + ,0 + ,68504 + ,385534 + ,25 + ,121 + ,119182 + ,33170 + ,18 + ,1 + ,22807 + ,101645 + ,11 + ,20 + ,17140 + ,149061 + ,26 + ,43 + ,116174 + ,165446 + ,25 + ,69 + ,57635 + ,237213 + ,38 + ,78 + ,66198 + ,173326 + ,44 + ,86 + ,71701 + ,133131 + ,30 + ,44 + ,57793 + ,258873 + ,40 + ,104 + ,80444 + ,180083 + ,34 + ,63 + ,53855 + ,324799 + ,47 + ,158 + ,97668 + ,230964 + ,30 + ,102 + ,133824 + ,236785 + ,31 + ,77 + ,101481 + ,135473 + ,23 + ,82 + ,99645 + ,202925 + ,36 + ,115 + ,114789 + ,215147 + ,36 + ,101 + ,99052 + ,344297 + ,30 + ,80 + ,67654 + ,153935 + ,25 + ,50 + ,65553 + ,132943 + ,39 + ,83 + ,97500 + ,174724 + ,34 + ,123 + ,69112 + ,174415 + ,31 + ,73 + ,82753 + ,225548 + ,31 + ,81 + ,85323 + ,223632 + ,33 + ,105 + ,72654 + ,124817 + ,25 + ,47 + ,30727 + ,221698 + ,33 + ,105 + ,77873 + ,210767 + ,35 + ,94 + ,117478 + ,170266 + ,42 + ,44 + ,74007 + ,260561 + ,43 + ,114 + ,90183 + ,84853 + ,30 + ,38 + ,61542 + ,294424 + ,33 + ,107 + ,101494 + ,101011 + ,13 + ,30 + ,27570 + ,215641 + ,32 + ,71 + ,55813 + ,325107 + ,36 + ,84 + ,79215 + ,7176 + ,0 + ,0 + ,1423 + ,167542 + ,28 + ,59 + ,55461 + ,106408 + ,14 + ,33 + ,31081 + ,96560 + ,17 + ,42 + ,22996 + ,265769 + ,32 + ,96 + ,83122 + ,269651 + ,30 + ,106 + ,70106 + ,149112 + ,35 + ,56 + ,60578 + ,175824 + ,20 + ,57 + ,39992 + ,152871 + ,28 + ,59 + ,79892 + ,111665 + ,28 + ,39 + ,49810 + ,116408 + ,39 + ,34 + ,71570 + ,362301 + ,34 + ,76 + ,100708 + ,78800 + ,26 + ,20 + ,33032 + ,183167 + ,39 + ,91 + ,82875 + ,277965 + ,39 + ,115 + ,139077 + ,150629 + ,33 + ,85 + ,71595 + ,168809 + ,28 + ,76 + ,72260 + ,24188 + ,4 + ,8 + ,5950 + ,329267 + ,39 + ,79 + ,115762 + ,65029 + ,18 + ,21 + ,32551 + ,101097 + ,14 + ,30 + ,31701 + ,218946 + ,29 + ,76 + ,80670 + ,244052 + ,44 + ,101 + ,143558 + ,341570 + ,21 + ,94 + ,117105 + ,103597 + ,16 + ,27 + ,23789 + ,233328 + ,28 + ,92 + ,120733 + ,256462 + ,35 + ,123 + ,105195 + ,206161 + ,28 + ,75 + ,73107 + ,311473 + ,38 + ,128 + ,132068 + ,235800 + ,23 + ,105 + ,149193 + ,177939 + ,36 + ,55 + ,46821 + ,207176 + ,32 + ,56 + ,87011 + ,196553 + ,29 + ,41 + ,95260 + ,174184 + ,25 + ,72 + ,55183 + ,143246 + ,27 + ,67 + ,106671 + ,187559 + ,36 + ,75 + ,73511 + ,187681 + ,28 + ,114 + ,92945 + ,119016 + ,23 + ,118 + ,78664 + ,182192 + ,40 + ,77 + ,70054 + ,73566 + ,23 + ,22 + ,22618 + ,194979 + ,40 + ,66 + ,74011 + ,167488 + ,28 + ,69 + ,83737 + ,143756 + ,34 + ,105 + ,69094 + ,275541 + ,33 + ,116 + ,93133 + ,243199 + ,28 + ,88 + ,95536 + ,182999 + ,34 + ,73 + ,225920 + ,135649 + ,30 + ,99 + ,62133 + ,152299 + ,33 + ,62 + ,61370 + ,120221 + ,22 + ,53 + ,43836 + ,346485 + ,38 + ,118 + ,106117 + ,145790 + ,26 + ,30 + ,38692 + ,193339 + ,35 + ,100 + ,84651 + ,80953 + ,8 + ,49 + ,56622 + ,122774 + ,24 + ,24 + ,15986 + ,130585 + ,29 + ,67 + ,95364 + ,112611 + ,20 + ,46 + ,26706 + ,286468 + ,29 + ,57 + ,89691 + ,241066 + ,45 + ,75 + ,67267 + ,148446 + ,37 + ,135 + ,126846 + ,204713 + ,33 + ,68 + ,41140 + ,182079 + ,33 + ,124 + ,102860 + ,140344 + ,25 + ,33 + ,51715 + ,220516 + ,32 + ,98 + ,55801 + ,243060 + ,29 + ,58 + ,111813 + ,162765 + ,28 + ,68 + ,120293 + ,182613 + ,28 + ,81 + ,138599 + ,232138 + ,31 + ,131 + ,161647 + ,265318 + ,52 + ,110 + ,115929 + ,85574 + ,21 + ,37 + ,24266 + ,310839 + ,24 + ,130 + ,162901 + ,225060 + ,41 + ,93 + ,109825 + ,232317 + ,33 + ,118 + ,129838 + ,144966 + ,32 + ,39 + ,37510 + ,43287 + ,19 + ,13 + ,43750 + ,155754 + ,20 + ,74 + ,40652 + ,164709 + ,31 + ,81 + ,87771 + ,201940 + ,31 + ,109 + ,85872 + ,235454 + ,32 + ,151 + ,89275 + ,220801 + ,18 + ,51 + ,44418 + ,99466 + ,23 + ,28 + ,192565 + ,92661 + ,17 + ,40 + ,35232 + ,133328 + ,20 + ,56 + ,40909 + ,61361 + ,12 + ,27 + ,13294 + ,125930 + ,17 + ,37 + ,32387 + ,100750 + ,30 + ,83 + ,140867 + ,224549 + ,31 + ,54 + ,120662 + ,82316 + ,10 + ,27 + ,21233 + ,102010 + ,13 + ,28 + ,44332 + ,101523 + ,22 + ,59 + ,61056 + ,243511 + ,42 + ,133 + ,101338 + ,22938 + ,1 + ,12 + ,1168 + ,41566 + ,9 + ,0 + ,13497 + ,152474 + ,32 + ,106 + ,65567 + ,61857 + ,11 + ,23 + ,25162 + ,99923 + ,25 + ,44 + ,32334 + ,132487 + ,36 + ,71 + ,40735 + ,317394 + ,31 + ,116 + ,91413 + ,21054 + ,0 + ,4 + ,855 + ,209641 + ,24 + ,62 + ,97068 + ,22648 + ,13 + ,12 + ,44339 + ,31414 + ,8 + ,18 + ,14116 + ,46698 + ,13 + ,14 + ,10288 + ,131698 + ,19 + ,60 + ,65622 + ,91735 + ,18 + ,7 + ,16563 + ,244749 + ,33 + ,98 + ,76643 + ,184510 + ,40 + ,64 + ,110681 + ,79863 + ,22 + ,29 + ,29011 + ,128423 + ,38 + ,32 + ,92696 + ,97839 + ,24 + ,25 + ,94785 + ,38214 + ,8 + ,16 + ,8773 + ,151101 + ,35 + ,48 + ,83209 + ,272458 + ,43 + ,100 + ,93815 + ,172494 + ,43 + ,46 + ,86687 + ,108043 + ,14 + ,45 + ,34553 + ,328107 + ,41 + ,129 + ,105547 + ,250579 + ,38 + ,130 + ,103487 + ,351067 + ,45 + ,136 + ,213688 + ,158015 + ,31 + ,59 + ,71220 + ,98866 + ,13 + ,25 + ,23517 + ,85439 + ,28 + ,32 + ,56926 + ,229242 + ,31 + ,63 + ,91721 + ,351619 + ,40 + ,95 + ,115168 + ,84207 + ,30 + ,14 + ,111194 + ,120445 + ,16 + ,36 + ,51009 + ,324598 + ,37 + ,113 + ,135777 + ,131069 + ,30 + ,47 + ,51513 + ,204271 + ,35 + ,92 + ,74163 + ,165543 + ,32 + ,70 + ,51633 + ,141722 + ,27 + ,19 + ,75345 + ,116048 + ,20 + ,50 + ,33416 + ,250047 + ,18 + ,41 + ,83305 + ,299775 + ,31 + ,91 + ,98952 + ,195838 + ,31 + ,111 + ,102372 + ,173260 + ,21 + ,41 + ,37238 + ,254488 + ,39 + ,120 + ,103772 + ,104389 + ,41 + ,135 + ,123969 + ,136084 + ,13 + ,27 + ,27142 + ,199476 + ,32 + ,87 + ,135400 + ,92499 + ,18 + ,25 + ,21399 + ,224330 + ,39 + ,131 + ,130115 + ,135781 + ,14 + ,45 + ,24874 + ,74408 + ,7 + ,29 + ,34988 + ,81240 + ,17 + ,58 + ,45549 + ,14688 + ,0 + ,4 + ,6023 + ,181633 + ,30 + ,47 + ,64466 + ,271856 + ,37 + ,109 + ,54990 + ,7199 + ,0 + ,7 + ,1644 + ,46660 + ,5 + ,12 + ,6179 + ,17547 + ,1 + ,0 + ,3926 + ,133368 + ,16 + ,37 + ,32755 + ,95227 + ,32 + ,37 + ,34777 + ,152601 + ,24 + ,46 + ,73224 + ,98146 + ,17 + ,15 + ,27114 + ,79619 + ,11 + ,42 + ,20760 + ,59194 + ,24 + ,7 + ,37636 + ,139942 + ,22 + ,54 + ,65461 + ,118612 + ,12 + ,54 + ,30080 + ,72880 + ,19 + ,14 + ,24094 + ,65475 + ,13 + ,16 + ,69008 + ,99643 + ,17 + ,33 + ,54968 + ,71965 + ,15 + ,32 + ,46090 + ,77272 + ,16 + ,21 + ,27507 + ,49289 + ,24 + ,15 + ,10672 + ,135131 + ,15 + ,38 + ,34029 + ,108446 + ,17 + ,22 + ,46300 + ,89746 + ,18 + ,28 + ,24760 + ,44296 + ,20 + ,10 + ,18779 + ,77648 + ,16 + ,31 + ,21280 + ,181528 + ,16 + ,32 + ,40662 + ,134019 + ,18 + ,32 + ,28987 + ,124064 + ,22 + ,43 + ,22827 + ,92630 + ,8 + ,27 + ,18513 + ,121848 + ,17 + ,37 + ,30594 + ,52915 + ,18 + ,20 + ,24006 + ,81872 + ,16 + ,32 + ,27913 + ,58981 + ,23 + ,0 + ,42744 + ,53515 + ,22 + ,5 + ,12934 + ,60812 + ,13 + ,26 + ,22574 + ,56375 + ,13 + ,10 + ,41385 + ,65490 + ,16 + ,27 + ,18653 + ,80949 + ,16 + ,11 + ,18472 + ,76302 + ,20 + ,29 + ,30976 + ,104011 + ,22 + ,25 + ,63339 + ,98104 + ,17 + ,55 + ,25568 + ,67989 + ,18 + ,23 + ,33747 + ,30989 + ,17 + ,5 + ,4154 + ,135458 + ,12 + ,43 + ,19474 + ,73504 + ,7 + ,23 + ,35130 + ,63123 + ,17 + ,34 + ,39067 + ,61254 + ,14 + ,36 + ,13310 + ,74914 + ,23 + ,35 + ,65892 + ,31774 + ,17 + ,0 + ,4143 + ,81437 + ,14 + ,37 + ,28579 + ,87186 + ,15 + ,28 + ,51776 + ,50090 + ,17 + ,16 + ,21152 + ,65745 + ,21 + ,26 + ,38084 + ,56653 + ,18 + ,38 + ,27717 + ,158399 + ,18 + ,23 + ,32928 + ,46455 + ,17 + ,22 + ,11342 + ,73624 + ,17 + ,30 + ,19499 + ,38395 + ,16 + ,16 + ,16380 + ,91899 + ,15 + ,18 + ,36874 + ,139526 + ,21 + ,28 + ,48259 + ,52164 + ,16 + ,32 + ,16734 + ,51567 + ,14 + ,21 + ,28207 + ,70551 + ,15 + ,23 + ,30143 + ,84856 + ,17 + ,29 + ,41369 + ,102538 + ,15 + ,50 + ,45833 + ,86678 + ,15 + ,12 + ,29156 + ,85709 + ,10 + ,21 + ,35944 + ,34662 + ,6 + ,18 + ,36278 + ,150580 + ,22 + ,27 + ,45588 + ,99611 + ,21 + ,41 + ,45097 + ,19349 + ,1 + ,13 + ,3895 + ,99373 + ,18 + ,12 + ,28394 + ,86230 + ,17 + ,21 + ,18632 + ,30837 + ,4 + ,8 + ,2325 + ,31706 + ,10 + ,26 + ,25139 + ,89806 + ,16 + ,27 + ,27975 + ,62088 + ,16 + ,13 + ,14483 + ,40151 + ,9 + ,16 + ,13127 + ,27634 + ,16 + ,2 + ,5839 + ,76990 + ,17 + ,42 + ,24069 + ,37460 + ,7 + ,5 + ,3738 + ,54157 + ,15 + ,37 + ,18625 + ,49862 + ,14 + ,17 + ,36341 + ,84337 + ,14 + ,38 + ,24548 + ,64175 + ,18 + ,37 + ,21792 + ,59382 + ,12 + ,29 + ,26263 + ,119308 + ,16 + ,32 + ,23686 + ,76702 + ,21 + ,35 + ,49303 + ,103425 + ,19 + ,17 + ,25659 + ,70344 + ,16 + ,20 + ,28904 + ,43410 + ,1 + ,7 + ,2781 + ,104838 + ,16 + ,46 + ,29236 + ,62215 + ,10 + ,24 + ,19546 + ,69304 + ,19 + ,40 + ,22818 + ,53117 + ,12 + ,3 + ,32689 + ,19764 + ,2 + ,10 + ,5752 + ,86680 + ,14 + ,37 + ,22197 + ,84105 + ,17 + ,17 + ,20055 + ,77945 + ,19 + ,28 + ,25272 + ,89113 + ,14 + ,19 + ,82206 + ,91005 + ,11 + ,29 + ,32073 + ,40248 + ,4 + ,8 + ,5444 + ,64187 + ,16 + ,10 + ,20154 + ,50857 + ,20 + ,15 + ,36944 + ,56613 + ,12 + ,15 + ,8019 + ,62792 + ,15 + ,28 + ,30884 + ,72535 + ,16 + ,17 + ,19540) + ,dim=c(4 + ,289) + ,dimnames=list(c('time' + ,'feedback' + ,'blogs' + ,'dgddg') + ,1:289)) > y <- array(NA,dim=c(4,289),dimnames=list(c('time','feedback','blogs','dgddg'),1:289)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '1' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x time feedback blogs dgddg 1 210907 30 79 112285 2 120982 28 58 84786 3 176508 38 60 83123 4 179321 30 108 101193 5 123185 22 49 38361 6 52746 26 0 68504 7 385534 25 121 119182 8 33170 18 1 22807 9 101645 11 20 17140 10 149061 26 43 116174 11 165446 25 69 57635 12 237213 38 78 66198 13 173326 44 86 71701 14 133131 30 44 57793 15 258873 40 104 80444 16 180083 34 63 53855 17 324799 47 158 97668 18 230964 30 102 133824 19 236785 31 77 101481 20 135473 23 82 99645 21 202925 36 115 114789 22 215147 36 101 99052 23 344297 30 80 67654 24 153935 25 50 65553 25 132943 39 83 97500 26 174724 34 123 69112 27 174415 31 73 82753 28 225548 31 81 85323 29 223632 33 105 72654 30 124817 25 47 30727 31 221698 33 105 77873 32 210767 35 94 117478 33 170266 42 44 74007 34 260561 43 114 90183 35 84853 30 38 61542 36 294424 33 107 101494 37 101011 13 30 27570 38 215641 32 71 55813 39 325107 36 84 79215 40 7176 0 0 1423 41 167542 28 59 55461 42 106408 14 33 31081 43 96560 17 42 22996 44 265769 32 96 83122 45 269651 30 106 70106 46 149112 35 56 60578 47 175824 20 57 39992 48 152871 28 59 79892 49 111665 28 39 49810 50 116408 39 34 71570 51 362301 34 76 100708 52 78800 26 20 33032 53 183167 39 91 82875 54 277965 39 115 139077 55 150629 33 85 71595 56 168809 28 76 72260 57 24188 4 8 5950 58 329267 39 79 115762 59 65029 18 21 32551 60 101097 14 30 31701 61 218946 29 76 80670 62 244052 44 101 143558 63 341570 21 94 117105 64 103597 16 27 23789 65 233328 28 92 120733 66 256462 35 123 105195 67 206161 28 75 73107 68 311473 38 128 132068 69 235800 23 105 149193 70 177939 36 55 46821 71 207176 32 56 87011 72 196553 29 41 95260 73 174184 25 72 55183 74 143246 27 67 106671 75 187559 36 75 73511 76 187681 28 114 92945 77 119016 23 118 78664 78 182192 40 77 70054 79 73566 23 22 22618 80 194979 40 66 74011 81 167488 28 69 83737 82 143756 34 105 69094 83 275541 33 116 93133 84 243199 28 88 95536 85 182999 34 73 225920 86 135649 30 99 62133 87 152299 33 62 61370 88 120221 22 53 43836 89 346485 38 118 106117 90 145790 26 30 38692 91 193339 35 100 84651 92 80953 8 49 56622 93 122774 24 24 15986 94 130585 29 67 95364 95 112611 20 46 26706 96 286468 29 57 89691 97 241066 45 75 67267 98 148446 37 135 126846 99 204713 33 68 41140 100 182079 33 124 102860 101 140344 25 33 51715 102 220516 32 98 55801 103 243060 29 58 111813 104 162765 28 68 120293 105 182613 28 81 138599 106 232138 31 131 161647 107 265318 52 110 115929 108 85574 21 37 24266 109 310839 24 130 162901 110 225060 41 93 109825 111 232317 33 118 129838 112 144966 32 39 37510 113 43287 19 13 43750 114 155754 20 74 40652 115 164709 31 81 87771 116 201940 31 109 85872 117 235454 32 151 89275 118 220801 18 51 44418 119 99466 23 28 192565 120 92661 17 40 35232 121 133328 20 56 40909 122 61361 12 27 13294 123 125930 17 37 32387 124 100750 30 83 140867 125 224549 31 54 120662 126 82316 10 27 21233 127 102010 13 28 44332 128 101523 22 59 61056 129 243511 42 133 101338 130 22938 1 12 1168 131 41566 9 0 13497 132 152474 32 106 65567 133 61857 11 23 25162 134 99923 25 44 32334 135 132487 36 71 40735 136 317394 31 116 91413 137 21054 0 4 855 138 209641 24 62 97068 139 22648 13 12 44339 140 31414 8 18 14116 141 46698 13 14 10288 142 131698 19 60 65622 143 91735 18 7 16563 144 244749 33 98 76643 145 184510 40 64 110681 146 79863 22 29 29011 147 128423 38 32 92696 148 97839 24 25 94785 149 38214 8 16 8773 150 151101 35 48 83209 151 272458 43 100 93815 152 172494 43 46 86687 153 108043 14 45 34553 154 328107 41 129 105547 155 250579 38 130 103487 156 351067 45 136 213688 157 158015 31 59 71220 158 98866 13 25 23517 159 85439 28 32 56926 160 229242 31 63 91721 161 351619 40 95 115168 162 84207 30 14 111194 163 120445 16 36 51009 164 324598 37 113 135777 165 131069 30 47 51513 166 204271 35 92 74163 167 165543 32 70 51633 168 141722 27 19 75345 169 116048 20 50 33416 170 250047 18 41 83305 171 299775 31 91 98952 172 195838 31 111 102372 173 173260 21 41 37238 174 254488 39 120 103772 175 104389 41 135 123969 176 136084 13 27 27142 177 199476 32 87 135400 178 92499 18 25 21399 179 224330 39 131 130115 180 135781 14 45 24874 181 74408 7 29 34988 182 81240 17 58 45549 183 14688 0 4 6023 184 181633 30 47 64466 185 271856 37 109 54990 186 7199 0 7 1644 187 46660 5 12 6179 188 17547 1 0 3926 189 133368 16 37 32755 190 95227 32 37 34777 191 152601 24 46 73224 192 98146 17 15 27114 193 79619 11 42 20760 194 59194 24 7 37636 195 139942 22 54 65461 196 118612 12 54 30080 197 72880 19 14 24094 198 65475 13 16 69008 199 99643 17 33 54968 200 71965 15 32 46090 201 77272 16 21 27507 202 49289 24 15 10672 203 135131 15 38 34029 204 108446 17 22 46300 205 89746 18 28 24760 206 44296 20 10 18779 207 77648 16 31 21280 208 181528 16 32 40662 209 134019 18 32 28987 210 124064 22 43 22827 211 92630 8 27 18513 212 121848 17 37 30594 213 52915 18 20 24006 214 81872 16 32 27913 215 58981 23 0 42744 216 53515 22 5 12934 217 60812 13 26 22574 218 56375 13 10 41385 219 65490 16 27 18653 220 80949 16 11 18472 221 76302 20 29 30976 222 104011 22 25 63339 223 98104 17 55 25568 224 67989 18 23 33747 225 30989 17 5 4154 226 135458 12 43 19474 227 73504 7 23 35130 228 63123 17 34 39067 229 61254 14 36 13310 230 74914 23 35 65892 231 31774 17 0 4143 232 81437 14 37 28579 233 87186 15 28 51776 234 50090 17 16 21152 235 65745 21 26 38084 236 56653 18 38 27717 237 158399 18 23 32928 238 46455 17 22 11342 239 73624 17 30 19499 240 38395 16 16 16380 241 91899 15 18 36874 242 139526 21 28 48259 243 52164 16 32 16734 244 51567 14 21 28207 245 70551 15 23 30143 246 84856 17 29 41369 247 102538 15 50 45833 248 86678 15 12 29156 249 85709 10 21 35944 250 34662 6 18 36278 251 150580 22 27 45588 252 99611 21 41 45097 253 19349 1 13 3895 254 99373 18 12 28394 255 86230 17 21 18632 256 30837 4 8 2325 257 31706 10 26 25139 258 89806 16 27 27975 259 62088 16 13 14483 260 40151 9 16 13127 261 27634 16 2 5839 262 76990 17 42 24069 263 37460 7 5 3738 264 54157 15 37 18625 265 49862 14 17 36341 266 84337 14 38 24548 267 64175 18 37 21792 268 59382 12 29 26263 269 119308 16 32 23686 270 76702 21 35 49303 271 103425 19 17 25659 272 70344 16 20 28904 273 43410 1 7 2781 274 104838 16 46 29236 275 62215 10 24 19546 276 69304 19 40 22818 277 53117 12 3 32689 278 19764 2 10 5752 279 86680 14 37 22197 280 84105 17 17 20055 281 77945 19 28 25272 282 89113 14 19 82206 283 91005 11 29 32073 284 40248 4 8 5444 285 64187 16 10 20154 286 50857 20 15 36944 287 56613 12 15 8019 288 62792 15 28 30884 289 72535 16 17 19540 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) feedback blogs dgddg 1.470e+04 1.774e+03 1.148e+03 3.816e-01 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -185310 -22442 -5975 18293 161620 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.470e+04 6.091e+03 2.413 0.016437 * feedback 1.774e+03 3.899e+02 4.549 7.97e-06 *** blogs 1.148e+03 1.170e+02 9.811 < 2e-16 *** dgddg 3.816e-01 1.019e-01 3.747 0.000217 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 42190 on 285 degrees of freedom Multiple R-squared: 0.7402, Adjusted R-squared: 0.7375 F-statistic: 270.7 on 3 and 285 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.9633273 7.334541e-02 3.667271e-02 [2,] 0.9257279 1.485443e-01 7.427213e-02 [3,] 0.8768804 2.462392e-01 1.231196e-01 [4,] 0.8200441 3.599118e-01 1.799559e-01 [5,] 0.7383030 5.233941e-01 2.616970e-01 [6,] 0.8112905 3.774190e-01 1.887095e-01 [7,] 0.7594717 4.810566e-01 2.405283e-01 [8,] 0.6868835 6.262331e-01 3.131165e-01 [9,] 0.6200817 7.598367e-01 3.799183e-01 [10,] 0.5574093 8.851814e-01 4.425907e-01 [11,] 0.4975217 9.950434e-01 5.024783e-01 [12,] 0.4600968 9.201937e-01 5.399032e-01 [13,] 0.4376141 8.752281e-01 5.623859e-01 [14,] 0.6120534 7.758932e-01 3.879466e-01 [15,] 0.6763612 6.472777e-01 3.236388e-01 [16,] 0.6224309 7.551381e-01 3.775691e-01 [17,] 0.9700739 5.985224e-02 2.992612e-02 [18,] 0.9581349 8.373013e-02 4.186507e-02 [19,] 0.9710654 5.786914e-02 2.893457e-02 [20,] 0.9904115 1.917696e-02 9.588478e-03 [21,] 0.9861807 2.763861e-02 1.381930e-02 [22,] 0.9827731 3.445381e-02 1.722691e-02 [23,] 0.9761684 4.766321e-02 2.383160e-02 [24,] 0.9675975 6.480504e-02 3.240252e-02 [25,] 0.9570187 8.596255e-02 4.298127e-02 [26,] 0.9441225 1.117549e-01 5.587747e-02 [27,] 0.9387014 1.225971e-01 6.129856e-02 [28,] 0.9223384 1.553233e-01 7.766164e-02 [29,] 0.9193413 1.613174e-01 8.065870e-02 [30,] 0.9272756 1.454487e-01 7.272436e-02 [31,] 0.9087922 1.824156e-01 9.120778e-02 [32,] 0.9046413 1.907174e-01 9.535870e-02 [33,] 0.9805305 3.893892e-02 1.946946e-02 [34,] 0.9775209 4.495829e-02 2.247914e-02 [35,] 0.9707514 5.849715e-02 2.924858e-02 [36,] 0.9622262 7.554764e-02 3.777382e-02 [37,] 0.9535637 9.287251e-02 4.643625e-02 [38,] 0.9528895 9.422097e-02 4.711048e-02 [39,] 0.9493447 1.013107e-01 5.065533e-02 [40,] 0.9369971 1.260059e-01 6.300293e-02 [41,] 0.9299416 1.401169e-01 7.005844e-02 [42,] 0.9142877 1.714246e-01 8.571230e-02 [43,] 0.8971648 2.056705e-01 1.028352e-01 [44,] 0.8797056 2.405888e-01 1.202944e-01 [45,] 0.9940341 1.193187e-02 5.965935e-03 [46,] 0.9921493 1.570142e-02 7.850710e-03 [47,] 0.9916068 1.678635e-02 8.393177e-03 [48,] 0.9889503 2.209941e-02 1.104971e-02 [49,] 0.9902236 1.955285e-02 9.776423e-03 [50,] 0.9877638 2.447239e-02 1.223619e-02 [51,] 0.9849806 3.003879e-02 1.501939e-02 [52,] 0.9964427 7.114664e-03 3.557332e-03 [53,] 0.9954115 9.177044e-03 4.588522e-03 [54,] 0.9939692 1.206169e-02 6.030843e-03 [55,] 0.9928776 1.424489e-02 7.122443e-03 [56,] 0.9914190 1.716195e-02 8.580977e-03 [57,] 0.9984559 3.088292e-03 1.544146e-03 [58,] 0.9979859 4.028297e-03 2.014148e-03 [59,] 0.9974259 5.148169e-03 2.574085e-03 [60,] 0.9968412 6.317583e-03 3.158792e-03 [61,] 0.9960796 7.840718e-03 3.920359e-03 [62,] 0.9952936 9.412861e-03 4.706431e-03 [63,] 0.9949850 1.003003e-02 5.015017e-03 [64,] 0.9939767 1.204652e-02 6.023260e-03 [65,] 0.9935483 1.290335e-02 6.451673e-03 [66,] 0.9936621 1.267576e-02 6.337879e-03 [67,] 0.9918799 1.624014e-02 8.120070e-03 [68,] 0.9922839 1.543223e-02 7.716117e-03 [69,] 0.9900850 1.982993e-02 9.914966e-03 [70,] 0.9920044 1.599126e-02 7.995630e-03 [71,] 0.9983734 3.253165e-03 1.626582e-03 [72,] 0.9979053 4.189477e-03 2.094739e-03 [73,] 0.9972958 5.408389e-03 2.704194e-03 [74,] 0.9964460 7.107905e-03 3.553952e-03 [75,] 0.9954686 9.062800e-03 4.531400e-03 [76,] 0.9975266 4.946870e-03 2.473435e-03 [77,] 0.9971925 5.615033e-03 2.807516e-03 [78,] 0.9969937 6.012553e-03 3.006277e-03 [79,] 0.9982258 3.548397e-03 1.774198e-03 [80,] 0.9989451 2.109741e-03 1.054870e-03 [81,] 0.9986283 2.743328e-03 1.371664e-03 [82,] 0.9982163 3.567464e-03 1.783732e-03 [83,] 0.9992817 1.436545e-03 7.182723e-04 [84,] 0.9992102 1.579503e-03 7.897516e-04 [85,] 0.9990926 1.814731e-03 9.073656e-04 [86,] 0.9989594 2.081250e-03 1.040625e-03 [87,] 0.9988095 2.380928e-03 1.190464e-03 [88,] 0.9989193 2.161383e-03 1.080692e-03 [89,] 0.9985568 2.886433e-03 1.443216e-03 [90,] 0.9998211 3.577045e-04 1.788522e-04 [91,] 0.9998047 3.906925e-04 1.953463e-04 [92,] 0.9999918 1.648890e-05 8.244452e-06 [93,] 0.9999911 1.785267e-05 8.926334e-06 [94,] 0.9999956 8.884210e-06 4.442105e-06 [95,] 0.9999941 1.184290e-05 5.921452e-06 [96,] 0.9999917 1.652321e-05 8.261606e-06 [97,] 0.9999954 9.242968e-06 4.621484e-06 [98,] 0.9999941 1.175476e-05 5.877380e-06 [99,] 0.9999926 1.470804e-05 7.354018e-06 [100,] 0.9999934 1.328948e-05 6.644742e-06 [101,] 0.9999904 1.923762e-05 9.618808e-06 [102,] 0.9999871 2.579673e-05 1.289837e-05 [103,] 0.9999872 2.555532e-05 1.277766e-05 [104,] 0.9999817 3.652449e-05 1.826224e-05 [105,] 0.9999765 4.703084e-05 2.351542e-05 [106,] 0.9999677 6.456489e-05 3.228244e-05 [107,] 0.9999660 6.803971e-05 3.401985e-05 [108,] 0.9999517 9.654824e-05 4.827412e-05 [109,] 0.9999429 1.141246e-04 5.706231e-05 [110,] 0.9999281 1.438702e-04 7.193509e-05 [111,] 0.9999277 1.446899e-04 7.234493e-05 [112,] 0.9999843 3.132663e-05 1.566331e-05 [113,] 0.9999907 1.850725e-05 9.253624e-06 [114,] 0.9999871 2.586945e-05 1.293473e-05 [115,] 0.9999813 3.732450e-05 1.866225e-05 [116,] 0.9999742 5.162130e-05 2.581065e-05 [117,] 0.9999675 6.494151e-05 3.247075e-05 [118,] 0.9999981 3.764178e-06 1.882089e-06 [119,] 0.9999982 3.671061e-06 1.835531e-06 [120,] 0.9999973 5.403397e-06 2.701699e-06 [121,] 0.9999961 7.820873e-06 3.910437e-06 [122,] 0.9999964 7.256887e-06 3.628443e-06 [123,] 0.9999959 8.128161e-06 4.064080e-06 [124,] 0.9999942 1.160523e-05 5.802613e-06 [125,] 0.9999915 1.691928e-05 8.459638e-06 [126,] 0.9999953 9.371661e-06 4.685831e-06 [127,] 0.9999933 1.345890e-05 6.729449e-06 [128,] 0.9999911 1.781234e-05 8.906169e-06 [129,] 0.9999910 1.797880e-05 8.989398e-06 [130,] 0.9999967 6.635661e-06 3.317831e-06 [131,] 0.9999951 9.828142e-06 4.914071e-06 [132,] 0.9999951 9.811005e-06 4.905502e-06 [133,] 0.9999958 8.416433e-06 4.208217e-06 [134,] 0.9999946 1.079116e-05 5.395581e-06 [135,] 0.9999923 1.546216e-05 7.731081e-06 [136,] 0.9999891 2.186850e-05 1.093425e-05 [137,] 0.9999874 2.514843e-05 1.257421e-05 [138,] 0.9999852 2.962176e-05 1.481088e-05 [139,] 0.9999799 4.024987e-05 2.012494e-05 [140,] 0.9999728 5.443902e-05 2.721951e-05 [141,] 0.9999663 6.730285e-05 3.365143e-05 [142,] 0.9999603 7.940902e-05 3.970451e-05 [143,] 0.9999453 1.093132e-04 5.465661e-05 [144,] 0.9999256 1.488948e-04 7.444742e-05 [145,] 0.9999165 1.670369e-04 8.351843e-05 [146,] 0.9998829 2.341095e-04 1.170547e-04 [147,] 0.9998372 3.256524e-04 1.628262e-04 [148,] 0.9998808 2.383041e-04 1.191520e-04 [149,] 0.9998417 3.165333e-04 1.582666e-04 [150,] 0.9997895 4.210215e-04 2.105107e-04 [151,] 0.9997111 5.777521e-04 2.888760e-04 [152,] 0.9996440 7.119719e-04 3.559859e-04 [153,] 0.9996270 7.459820e-04 3.729910e-04 [154,] 0.9996886 6.227039e-04 3.113519e-04 [155,] 0.9999771 4.576605e-05 2.288302e-05 [156,] 0.9999818 3.640790e-05 1.820395e-05 [157,] 0.9999749 5.024912e-05 2.512456e-05 [158,] 0.9999891 2.189695e-05 1.094848e-05 [159,] 0.9999840 3.193549e-05 1.596774e-05 [160,] 0.9999776 4.479505e-05 2.239752e-05 [161,] 0.9999684 6.321353e-05 3.160677e-05 [162,] 0.9999612 7.769866e-05 3.884933e-05 [163,] 0.9999445 1.109972e-04 5.549862e-05 [164,] 0.9999987 2.506967e-06 1.253484e-06 [165,] 1.0000000 9.460124e-08 4.730062e-08 [166,] 0.9999999 1.331074e-07 6.655368e-08 [167,] 1.0000000 4.346300e-08 2.173150e-08 [168,] 1.0000000 4.908674e-08 2.454337e-08 [169,] 1.0000000 8.367739e-13 4.183870e-13 [170,] 1.0000000 2.237062e-13 1.118531e-13 [171,] 1.0000000 3.804879e-13 1.902439e-13 [172,] 1.0000000 6.934349e-13 3.467174e-13 [173,] 1.0000000 1.239825e-13 6.199126e-14 [174,] 1.0000000 1.141113e-13 5.705567e-14 [175,] 1.0000000 2.330708e-13 1.165354e-13 [176,] 1.0000000 8.937278e-14 4.468639e-14 [177,] 1.0000000 1.812949e-13 9.064746e-14 [178,] 1.0000000 1.750911e-13 8.754557e-14 [179,] 1.0000000 7.559650e-14 3.779825e-14 [180,] 1.0000000 1.258091e-13 6.290454e-14 [181,] 1.0000000 2.529540e-13 1.264770e-13 [182,] 1.0000000 5.195727e-13 2.597864e-13 [183,] 1.0000000 3.959381e-13 1.979690e-13 [184,] 1.0000000 6.174988e-13 3.087494e-13 [185,] 1.0000000 1.062105e-12 5.310523e-13 [186,] 1.0000000 1.201499e-12 6.007495e-13 [187,] 1.0000000 2.355531e-12 1.177766e-12 [188,] 1.0000000 4.160587e-12 2.080294e-12 [189,] 1.0000000 8.378632e-12 4.189316e-12 [190,] 1.0000000 1.508320e-11 7.541601e-12 [191,] 1.0000000 2.950234e-11 1.475117e-11 [192,] 1.0000000 4.255835e-11 2.127918e-11 [193,] 1.0000000 8.318434e-11 4.159217e-11 [194,] 1.0000000 1.160628e-10 5.803142e-11 [195,] 1.0000000 2.248085e-10 1.124043e-10 [196,] 1.0000000 3.532603e-10 1.766301e-10 [197,] 1.0000000 2.474879e-10 1.237440e-10 [198,] 1.0000000 3.657206e-10 1.828603e-10 [199,] 1.0000000 6.794805e-10 3.397403e-10 [200,] 1.0000000 1.051239e-09 5.256193e-10 [201,] 1.0000000 1.997308e-09 9.986542e-10 [202,] 1.0000000 1.834332e-11 9.171661e-12 [203,] 1.0000000 6.875563e-12 3.437781e-12 [204,] 1.0000000 8.130955e-12 4.065478e-12 [205,] 1.0000000 8.190545e-12 4.095272e-12 [206,] 1.0000000 6.823110e-12 3.411555e-12 [207,] 1.0000000 1.097848e-11 5.489242e-12 [208,] 1.0000000 2.335513e-11 1.167756e-11 [209,] 1.0000000 4.166897e-11 2.083448e-11 [210,] 1.0000000 8.416908e-11 4.208454e-11 [211,] 1.0000000 1.675959e-10 8.379795e-11 [212,] 1.0000000 3.183636e-10 1.591818e-10 [213,] 1.0000000 6.344778e-10 3.172389e-10 [214,] 1.0000000 9.473124e-10 4.736562e-10 [215,] 1.0000000 1.801723e-09 9.008616e-10 [216,] 1.0000000 3.666680e-09 1.833340e-09 [217,] 1.0000000 7.174948e-09 3.587474e-09 [218,] 1.0000000 1.286365e-08 6.431824e-09 [219,] 1.0000000 1.999768e-08 9.998842e-09 [220,] 1.0000000 2.686172e-09 1.343086e-09 [221,] 1.0000000 5.010112e-09 2.505056e-09 [222,] 1.0000000 6.400808e-09 3.200404e-09 [223,] 1.0000000 1.242744e-08 6.213718e-09 [224,] 1.0000000 5.774308e-09 2.887154e-09 [225,] 1.0000000 9.135014e-09 4.567507e-09 [226,] 1.0000000 1.864724e-08 9.323622e-09 [227,] 1.0000000 3.791390e-08 1.895695e-08 [228,] 1.0000000 5.610118e-08 2.805059e-08 [229,] 1.0000000 5.719475e-08 2.859738e-08 [230,] 1.0000000 5.508622e-08 2.754311e-08 [231,] 1.0000000 8.407517e-10 4.203759e-10 [232,] 1.0000000 1.225375e-09 6.126876e-10 [233,] 1.0000000 2.817013e-09 1.408507e-09 [234,] 1.0000000 2.782753e-09 1.391377e-09 [235,] 1.0000000 5.178761e-09 2.589380e-09 [236,] 1.0000000 2.294107e-09 1.147053e-09 [237,] 1.0000000 2.974851e-09 1.487426e-09 [238,] 1.0000000 4.743892e-09 2.371946e-09 [239,] 1.0000000 1.115982e-08 5.579910e-09 [240,] 1.0000000 2.623947e-08 1.311974e-08 [241,] 1.0000000 5.701166e-08 2.850583e-08 [242,] 1.0000000 9.470399e-08 4.735200e-08 [243,] 0.9999999 1.275306e-07 6.376529e-08 [244,] 0.9999999 2.144321e-07 1.072160e-07 [245,] 1.0000000 1.441051e-08 7.205254e-09 [246,] 1.0000000 3.671126e-08 1.835563e-08 [247,] 1.0000000 7.811165e-08 3.905582e-08 [248,] 1.0000000 5.583792e-08 2.791896e-08 [249,] 1.0000000 9.224648e-08 4.612324e-08 [250,] 0.9999999 2.340378e-07 1.170189e-07 [251,] 0.9999999 1.249486e-07 6.247432e-08 [252,] 0.9999999 2.321923e-07 1.160961e-07 [253,] 0.9999997 5.987794e-07 2.993897e-07 [254,] 0.9999994 1.245337e-06 6.226685e-07 [255,] 0.9999992 1.678587e-06 8.392934e-07 [256,] 0.9999980 3.949302e-06 1.974651e-06 [257,] 0.9999951 9.808251e-06 4.904125e-06 [258,] 0.9999948 1.038637e-05 5.193186e-06 [259,] 0.9999918 1.645555e-05 8.227776e-06 [260,] 0.9999791 4.173944e-05 2.086972e-05 [261,] 0.9999733 5.331013e-05 2.665507e-05 [262,] 0.9999506 9.883743e-05 4.941872e-05 [263,] 0.9999869 2.616793e-05 1.308396e-05 [264,] 0.9999812 3.769588e-05 1.884794e-05 [265,] 0.9999944 1.118467e-05 5.592336e-06 [266,] 0.9999810 3.790612e-05 1.895306e-05 [267,] 0.9999483 1.034654e-04 5.173269e-05 [268,] 0.9998816 2.368044e-04 1.184022e-04 [269,] 0.9996250 7.499315e-04 3.749657e-04 [270,] 0.9995306 9.387744e-04 4.693872e-04 [271,] 0.9985531 2.893857e-03 1.446928e-03 [272,] 0.9983878 3.224401e-03 1.612200e-03 [273,] 0.9946113 1.077735e-02 5.388674e-03 [274,] 0.9932989 1.340220e-02 6.701101e-03 [275,] 0.9769704 4.605921e-02 2.302960e-02 [276,] 0.9370958 1.258085e-01 6.290424e-02 > postscript(file="/var/wessaorg/rcomp/tmp/17m4g1324127192.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/2kdvx1324127192.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/3650n1324127192.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/469sk1324127192.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/5i1h81324127192.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 9459.55117 -42318.42267 -6191.83910 -51181.44191 -1424.77691 6 7 8 9 10 -34217.21900 142113.49876 -23311.33777 37933.80275 -5452.25711 11 12 13 14 15 5202.08209 40310.65109 -45502.31467 -7345.85894 23141.98373 16 17 18 19 20 12204.28474 8093.56625 -5104.05047 39982.50566 -52177.46897 21 22 23 24 25 -51444.03012 -17146.35953 158733.95290 12478.74058 -83418.08288 26 27 28 29 30 -67848.80908 -10648.99555 30320.34322 2142.80883 94.76229 31 32 33 34 35 -1782.88705 -18749.30410 2315.48642 4313.23727 -50167.40189 36 37 38 39 40 59633.04958 18293.18762 41379.84449 119897.55241 -8067.82469 41 42 43 44 45 14284.84145 17132.88928 -5281.92931 52389.55026 53307.82657 46 47 48 49 50 -15071.18073 44956.57391 -9709.61554 -16478.38248 -33812.38589 51 52 53 54 55 161619.86137 -17583.46955 -36795.71870 9005.58453 -47498.83271 56 57 58 59 60 -10372.69665 -9061.65831 110528.14257 -18128.10104 15028.86516 61 62 63 64 65 34781.01161 -19416.58746 137029.71115 20444.19031 17282.06624 66 67 68 69 70 -1654.78295 27803.92871 32040.02889 2840.02156 18379.84862 71 72 73 74 75 38226.85155 46995.24938 11432.24145 -36963.17969 -5142.91568 76 77 78 79 80 -43013.30532 -101950.60841 -18581.69736 -15817.45803 5321.68653 81 82 83 84 85 -8038.56928 -78148.44311 33610.05983 41360.29097 -62022.46549 86 87 88 89 90 -69616.45843 -15525.92548 -11069.61226 88434.16372 35767.92906 91 92 93 94 95 -30536.89723 -25791.94407 31851.91723 -48856.82370 -559.69883 96 97 98 99 100 120669.73879 34782.44303 -135255.32333 37721.16333 -72746.88468 101 102 103 104 105 23682.28990 15267.17887 67671.59079 -25564.35686 -27624.55868 106 107 108 109 110 -49609.77672 -12128.11343 -18108.62133 42177.36240 -11030.87558 111 112 113 114 115 -25917.17252 14421.25911 -36734.85479 5121.06623 -31452.87244 116 117 118 119 120 -25637.27442 -43405.93363 98679.33854 -61660.48481 -11554.75924 121 122 123 124 125 3258.48559 -10691.33122 26243.54396 -116196.47094 46827.37628 126 127 128 129 130 10781.62257 15191.12727 -43226.34386 -37029.30306 -7756.67525 131 132 133 134 135 5749.94481 -65684.64897 -8359.16811 -21968.92476 -43115.36233 136 137 138 139 140 79667.12025 1435.49405 44157.36125 -45807.76922 -23525.94368 141 142 143 144 145 -11058.81079 -10620.19957 30749.35717 29772.53703 -16845.74526 146 147 148 149 150 -18221.37351 -25790.02002 -24302.53665 -12391.20456 -12535.82700 151 152 153 154 155 30894.23177 -4365.06184 3668.55944 52325.71859 -20242.49620 156 157 158 159 160 18886.15233 -6577.67244 23434.21451 -37384.99266 52234.20934 161 162 163 164 165 112967.21883 -42213.13655 16573.63549 62741.33280 -10454.84274 166 167 168 169 170 -6419.53683 -5975.10608 28564.89972 -4274.83988 124563.73369 171 172 173 174 175 87867.57984 -40331.79298 60035.50787 -6737.48329 -185309.72510 176 177 178 179 180 56973.10550 -23523.24308 9006.32656 -59575.07568 35100.29855 181 182 183 184 185 650.16261 -47574.47123 -6902.73900 35165.98104 45421.03807 186 187 188 189 190 -16164.19011 6957.67253 -425.86320 35314.93984 -31979.04053 191 192 193 194 195 14582.57981 25724.78799 -10726.61674 -20476.60749 -749.09312 196 197 198 199 200 9161.48249 -788.51354 -16986.49645 -4069.46484 -23663.87858 201 202 203 204 205 -412.52230 -29274.38424 37217.72342 20667.92260 1527.10537 206 207 208 209 210 -24526.57156 -9138.75869 86196.74319 39595.53682 12269.53655 211 212 213 214 215 25681.30704 22845.79586 -25833.26267 -8593.93189 -12830.08305 216 217 218 219 220 -10885.34145 -15407.77494 -8657.72993 -15702.78861 18191.06023 221 222 223 224 225 -18984.59768 -2582.31993 -19641.65844 -17920.24526 -21191.51468 226 227 228 229 230 42681.45731 6579.13757 -35669.12186 -24682.85265 -45906.04895 231 232 233 234 235 -14663.01217 -11474.73109 -6021.34917 -21203.83039 -30584.43213 236 237 238 239 240 -44172.96623 72802.30483 -27982.26411 -13109.05887 -29303.88693 241 242 243 244 245 15857.22161 37017.82138 -34035.75672 -22836.99206 -8661.37160 246 247 248 249 250 -9075.31576 -13654.29794 20468.76156 15447.71686 -25187.82865 251 252 253 254 255 48465.16659 -16612.67549 -13534.22494 28133.05847 10158.55754 256 257 258 259 260 -1029.27115 -40171.14001 5055.71207 -1443.36333 -13889.62917 261 262 263 264 265 -19972.13487 -25261.41191 3176.57626 -36729.87890 -23054.67831 266 267 268 269 270 -8184.26550 -33241.98285 -19915.33529 30455.19294 -34239.62063 271 272 273 274 275 25715.66205 -6725.79016 17839.06982 -2202.87352 -5231.99488 276 277 278 279 280 -33721.94544 1211.73279 -12158.15154 -3796.20644 12081.95025 281 282 283 284 285 -12243.11968 -3602.58689 11264.26255 7191.44346 1935.02952 286 287 288 289 -30637.07695 348.06676 -22442.45967 2482.32017 > postscript(file="/var/wessaorg/rcomp/tmp/6su7a1324127192.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 9459.55117 NA 1 -42318.42267 9459.55117 2 -6191.83910 -42318.42267 3 -51181.44191 -6191.83910 4 -1424.77691 -51181.44191 5 -34217.21900 -1424.77691 6 142113.49876 -34217.21900 7 -23311.33777 142113.49876 8 37933.80275 -23311.33777 9 -5452.25711 37933.80275 10 5202.08209 -5452.25711 11 40310.65109 5202.08209 12 -45502.31467 40310.65109 13 -7345.85894 -45502.31467 14 23141.98373 -7345.85894 15 12204.28474 23141.98373 16 8093.56625 12204.28474 17 -5104.05047 8093.56625 18 39982.50566 -5104.05047 19 -52177.46897 39982.50566 20 -51444.03012 -52177.46897 21 -17146.35953 -51444.03012 22 158733.95290 -17146.35953 23 12478.74058 158733.95290 24 -83418.08288 12478.74058 25 -67848.80908 -83418.08288 26 -10648.99555 -67848.80908 27 30320.34322 -10648.99555 28 2142.80883 30320.34322 29 94.76229 2142.80883 30 -1782.88705 94.76229 31 -18749.30410 -1782.88705 32 2315.48642 -18749.30410 33 4313.23727 2315.48642 34 -50167.40189 4313.23727 35 59633.04958 -50167.40189 36 18293.18762 59633.04958 37 41379.84449 18293.18762 38 119897.55241 41379.84449 39 -8067.82469 119897.55241 40 14284.84145 -8067.82469 41 17132.88928 14284.84145 42 -5281.92931 17132.88928 43 52389.55026 -5281.92931 44 53307.82657 52389.55026 45 -15071.18073 53307.82657 46 44956.57391 -15071.18073 47 -9709.61554 44956.57391 48 -16478.38248 -9709.61554 49 -33812.38589 -16478.38248 50 161619.86137 -33812.38589 51 -17583.46955 161619.86137 52 -36795.71870 -17583.46955 53 9005.58453 -36795.71870 54 -47498.83271 9005.58453 55 -10372.69665 -47498.83271 56 -9061.65831 -10372.69665 57 110528.14257 -9061.65831 58 -18128.10104 110528.14257 59 15028.86516 -18128.10104 60 34781.01161 15028.86516 61 -19416.58746 34781.01161 62 137029.71115 -19416.58746 63 20444.19031 137029.71115 64 17282.06624 20444.19031 65 -1654.78295 17282.06624 66 27803.92871 -1654.78295 67 32040.02889 27803.92871 68 2840.02156 32040.02889 69 18379.84862 2840.02156 70 38226.85155 18379.84862 71 46995.24938 38226.85155 72 11432.24145 46995.24938 73 -36963.17969 11432.24145 74 -5142.91568 -36963.17969 75 -43013.30532 -5142.91568 76 -101950.60841 -43013.30532 77 -18581.69736 -101950.60841 78 -15817.45803 -18581.69736 79 5321.68653 -15817.45803 80 -8038.56928 5321.68653 81 -78148.44311 -8038.56928 82 33610.05983 -78148.44311 83 41360.29097 33610.05983 84 -62022.46549 41360.29097 85 -69616.45843 -62022.46549 86 -15525.92548 -69616.45843 87 -11069.61226 -15525.92548 88 88434.16372 -11069.61226 89 35767.92906 88434.16372 90 -30536.89723 35767.92906 91 -25791.94407 -30536.89723 92 31851.91723 -25791.94407 93 -48856.82370 31851.91723 94 -559.69883 -48856.82370 95 120669.73879 -559.69883 96 34782.44303 120669.73879 97 -135255.32333 34782.44303 98 37721.16333 -135255.32333 99 -72746.88468 37721.16333 100 23682.28990 -72746.88468 101 15267.17887 23682.28990 102 67671.59079 15267.17887 103 -25564.35686 67671.59079 104 -27624.55868 -25564.35686 105 -49609.77672 -27624.55868 106 -12128.11343 -49609.77672 107 -18108.62133 -12128.11343 108 42177.36240 -18108.62133 109 -11030.87558 42177.36240 110 -25917.17252 -11030.87558 111 14421.25911 -25917.17252 112 -36734.85479 14421.25911 113 5121.06623 -36734.85479 114 -31452.87244 5121.06623 115 -25637.27442 -31452.87244 116 -43405.93363 -25637.27442 117 98679.33854 -43405.93363 118 -61660.48481 98679.33854 119 -11554.75924 -61660.48481 120 3258.48559 -11554.75924 121 -10691.33122 3258.48559 122 26243.54396 -10691.33122 123 -116196.47094 26243.54396 124 46827.37628 -116196.47094 125 10781.62257 46827.37628 126 15191.12727 10781.62257 127 -43226.34386 15191.12727 128 -37029.30306 -43226.34386 129 -7756.67525 -37029.30306 130 5749.94481 -7756.67525 131 -65684.64897 5749.94481 132 -8359.16811 -65684.64897 133 -21968.92476 -8359.16811 134 -43115.36233 -21968.92476 135 79667.12025 -43115.36233 136 1435.49405 79667.12025 137 44157.36125 1435.49405 138 -45807.76922 44157.36125 139 -23525.94368 -45807.76922 140 -11058.81079 -23525.94368 141 -10620.19957 -11058.81079 142 30749.35717 -10620.19957 143 29772.53703 30749.35717 144 -16845.74526 29772.53703 145 -18221.37351 -16845.74526 146 -25790.02002 -18221.37351 147 -24302.53665 -25790.02002 148 -12391.20456 -24302.53665 149 -12535.82700 -12391.20456 150 30894.23177 -12535.82700 151 -4365.06184 30894.23177 152 3668.55944 -4365.06184 153 52325.71859 3668.55944 154 -20242.49620 52325.71859 155 18886.15233 -20242.49620 156 -6577.67244 18886.15233 157 23434.21451 -6577.67244 158 -37384.99266 23434.21451 159 52234.20934 -37384.99266 160 112967.21883 52234.20934 161 -42213.13655 112967.21883 162 16573.63549 -42213.13655 163 62741.33280 16573.63549 164 -10454.84274 62741.33280 165 -6419.53683 -10454.84274 166 -5975.10608 -6419.53683 167 28564.89972 -5975.10608 168 -4274.83988 28564.89972 169 124563.73369 -4274.83988 170 87867.57984 124563.73369 171 -40331.79298 87867.57984 172 60035.50787 -40331.79298 173 -6737.48329 60035.50787 174 -185309.72510 -6737.48329 175 56973.10550 -185309.72510 176 -23523.24308 56973.10550 177 9006.32656 -23523.24308 178 -59575.07568 9006.32656 179 35100.29855 -59575.07568 180 650.16261 35100.29855 181 -47574.47123 650.16261 182 -6902.73900 -47574.47123 183 35165.98104 -6902.73900 184 45421.03807 35165.98104 185 -16164.19011 45421.03807 186 6957.67253 -16164.19011 187 -425.86320 6957.67253 188 35314.93984 -425.86320 189 -31979.04053 35314.93984 190 14582.57981 -31979.04053 191 25724.78799 14582.57981 192 -10726.61674 25724.78799 193 -20476.60749 -10726.61674 194 -749.09312 -20476.60749 195 9161.48249 -749.09312 196 -788.51354 9161.48249 197 -16986.49645 -788.51354 198 -4069.46484 -16986.49645 199 -23663.87858 -4069.46484 200 -412.52230 -23663.87858 201 -29274.38424 -412.52230 202 37217.72342 -29274.38424 203 20667.92260 37217.72342 204 1527.10537 20667.92260 205 -24526.57156 1527.10537 206 -9138.75869 -24526.57156 207 86196.74319 -9138.75869 208 39595.53682 86196.74319 209 12269.53655 39595.53682 210 25681.30704 12269.53655 211 22845.79586 25681.30704 212 -25833.26267 22845.79586 213 -8593.93189 -25833.26267 214 -12830.08305 -8593.93189 215 -10885.34145 -12830.08305 216 -15407.77494 -10885.34145 217 -8657.72993 -15407.77494 218 -15702.78861 -8657.72993 219 18191.06023 -15702.78861 220 -18984.59768 18191.06023 221 -2582.31993 -18984.59768 222 -19641.65844 -2582.31993 223 -17920.24526 -19641.65844 224 -21191.51468 -17920.24526 225 42681.45731 -21191.51468 226 6579.13757 42681.45731 227 -35669.12186 6579.13757 228 -24682.85265 -35669.12186 229 -45906.04895 -24682.85265 230 -14663.01217 -45906.04895 231 -11474.73109 -14663.01217 232 -6021.34917 -11474.73109 233 -21203.83039 -6021.34917 234 -30584.43213 -21203.83039 235 -44172.96623 -30584.43213 236 72802.30483 -44172.96623 237 -27982.26411 72802.30483 238 -13109.05887 -27982.26411 239 -29303.88693 -13109.05887 240 15857.22161 -29303.88693 241 37017.82138 15857.22161 242 -34035.75672 37017.82138 243 -22836.99206 -34035.75672 244 -8661.37160 -22836.99206 245 -9075.31576 -8661.37160 246 -13654.29794 -9075.31576 247 20468.76156 -13654.29794 248 15447.71686 20468.76156 249 -25187.82865 15447.71686 250 48465.16659 -25187.82865 251 -16612.67549 48465.16659 252 -13534.22494 -16612.67549 253 28133.05847 -13534.22494 254 10158.55754 28133.05847 255 -1029.27115 10158.55754 256 -40171.14001 -1029.27115 257 5055.71207 -40171.14001 258 -1443.36333 5055.71207 259 -13889.62917 -1443.36333 260 -19972.13487 -13889.62917 261 -25261.41191 -19972.13487 262 3176.57626 -25261.41191 263 -36729.87890 3176.57626 264 -23054.67831 -36729.87890 265 -8184.26550 -23054.67831 266 -33241.98285 -8184.26550 267 -19915.33529 -33241.98285 268 30455.19294 -19915.33529 269 -34239.62063 30455.19294 270 25715.66205 -34239.62063 271 -6725.79016 25715.66205 272 17839.06982 -6725.79016 273 -2202.87352 17839.06982 274 -5231.99488 -2202.87352 275 -33721.94544 -5231.99488 276 1211.73279 -33721.94544 277 -12158.15154 1211.73279 278 -3796.20644 -12158.15154 279 12081.95025 -3796.20644 280 -12243.11968 12081.95025 281 -3602.58689 -12243.11968 282 11264.26255 -3602.58689 283 7191.44346 11264.26255 284 1935.02952 7191.44346 285 -30637.07695 1935.02952 286 348.06676 -30637.07695 287 -22442.45967 348.06676 288 2482.32017 -22442.45967 289 NA 2482.32017 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -42318.42267 9459.55117 [2,] -6191.83910 -42318.42267 [3,] -51181.44191 -6191.83910 [4,] -1424.77691 -51181.44191 [5,] -34217.21900 -1424.77691 [6,] 142113.49876 -34217.21900 [7,] -23311.33777 142113.49876 [8,] 37933.80275 -23311.33777 [9,] -5452.25711 37933.80275 [10,] 5202.08209 -5452.25711 [11,] 40310.65109 5202.08209 [12,] -45502.31467 40310.65109 [13,] -7345.85894 -45502.31467 [14,] 23141.98373 -7345.85894 [15,] 12204.28474 23141.98373 [16,] 8093.56625 12204.28474 [17,] -5104.05047 8093.56625 [18,] 39982.50566 -5104.05047 [19,] -52177.46897 39982.50566 [20,] -51444.03012 -52177.46897 [21,] -17146.35953 -51444.03012 [22,] 158733.95290 -17146.35953 [23,] 12478.74058 158733.95290 [24,] -83418.08288 12478.74058 [25,] -67848.80908 -83418.08288 [26,] -10648.99555 -67848.80908 [27,] 30320.34322 -10648.99555 [28,] 2142.80883 30320.34322 [29,] 94.76229 2142.80883 [30,] -1782.88705 94.76229 [31,] -18749.30410 -1782.88705 [32,] 2315.48642 -18749.30410 [33,] 4313.23727 2315.48642 [34,] -50167.40189 4313.23727 [35,] 59633.04958 -50167.40189 [36,] 18293.18762 59633.04958 [37,] 41379.84449 18293.18762 [38,] 119897.55241 41379.84449 [39,] -8067.82469 119897.55241 [40,] 14284.84145 -8067.82469 [41,] 17132.88928 14284.84145 [42,] -5281.92931 17132.88928 [43,] 52389.55026 -5281.92931 [44,] 53307.82657 52389.55026 [45,] -15071.18073 53307.82657 [46,] 44956.57391 -15071.18073 [47,] -9709.61554 44956.57391 [48,] -16478.38248 -9709.61554 [49,] -33812.38589 -16478.38248 [50,] 161619.86137 -33812.38589 [51,] -17583.46955 161619.86137 [52,] -36795.71870 -17583.46955 [53,] 9005.58453 -36795.71870 [54,] -47498.83271 9005.58453 [55,] -10372.69665 -47498.83271 [56,] -9061.65831 -10372.69665 [57,] 110528.14257 -9061.65831 [58,] -18128.10104 110528.14257 [59,] 15028.86516 -18128.10104 [60,] 34781.01161 15028.86516 [61,] -19416.58746 34781.01161 [62,] 137029.71115 -19416.58746 [63,] 20444.19031 137029.71115 [64,] 17282.06624 20444.19031 [65,] -1654.78295 17282.06624 [66,] 27803.92871 -1654.78295 [67,] 32040.02889 27803.92871 [68,] 2840.02156 32040.02889 [69,] 18379.84862 2840.02156 [70,] 38226.85155 18379.84862 [71,] 46995.24938 38226.85155 [72,] 11432.24145 46995.24938 [73,] -36963.17969 11432.24145 [74,] -5142.91568 -36963.17969 [75,] -43013.30532 -5142.91568 [76,] -101950.60841 -43013.30532 [77,] -18581.69736 -101950.60841 [78,] -15817.45803 -18581.69736 [79,] 5321.68653 -15817.45803 [80,] -8038.56928 5321.68653 [81,] -78148.44311 -8038.56928 [82,] 33610.05983 -78148.44311 [83,] 41360.29097 33610.05983 [84,] -62022.46549 41360.29097 [85,] -69616.45843 -62022.46549 [86,] -15525.92548 -69616.45843 [87,] -11069.61226 -15525.92548 [88,] 88434.16372 -11069.61226 [89,] 35767.92906 88434.16372 [90,] -30536.89723 35767.92906 [91,] -25791.94407 -30536.89723 [92,] 31851.91723 -25791.94407 [93,] -48856.82370 31851.91723 [94,] -559.69883 -48856.82370 [95,] 120669.73879 -559.69883 [96,] 34782.44303 120669.73879 [97,] -135255.32333 34782.44303 [98,] 37721.16333 -135255.32333 [99,] -72746.88468 37721.16333 [100,] 23682.28990 -72746.88468 [101,] 15267.17887 23682.28990 [102,] 67671.59079 15267.17887 [103,] -25564.35686 67671.59079 [104,] -27624.55868 -25564.35686 [105,] -49609.77672 -27624.55868 [106,] -12128.11343 -49609.77672 [107,] -18108.62133 -12128.11343 [108,] 42177.36240 -18108.62133 [109,] -11030.87558 42177.36240 [110,] -25917.17252 -11030.87558 [111,] 14421.25911 -25917.17252 [112,] -36734.85479 14421.25911 [113,] 5121.06623 -36734.85479 [114,] -31452.87244 5121.06623 [115,] -25637.27442 -31452.87244 [116,] -43405.93363 -25637.27442 [117,] 98679.33854 -43405.93363 [118,] -61660.48481 98679.33854 [119,] -11554.75924 -61660.48481 [120,] 3258.48559 -11554.75924 [121,] -10691.33122 3258.48559 [122,] 26243.54396 -10691.33122 [123,] -116196.47094 26243.54396 [124,] 46827.37628 -116196.47094 [125,] 10781.62257 46827.37628 [126,] 15191.12727 10781.62257 [127,] -43226.34386 15191.12727 [128,] -37029.30306 -43226.34386 [129,] -7756.67525 -37029.30306 [130,] 5749.94481 -7756.67525 [131,] -65684.64897 5749.94481 [132,] -8359.16811 -65684.64897 [133,] -21968.92476 -8359.16811 [134,] -43115.36233 -21968.92476 [135,] 79667.12025 -43115.36233 [136,] 1435.49405 79667.12025 [137,] 44157.36125 1435.49405 [138,] -45807.76922 44157.36125 [139,] -23525.94368 -45807.76922 [140,] -11058.81079 -23525.94368 [141,] -10620.19957 -11058.81079 [142,] 30749.35717 -10620.19957 [143,] 29772.53703 30749.35717 [144,] -16845.74526 29772.53703 [145,] -18221.37351 -16845.74526 [146,] -25790.02002 -18221.37351 [147,] -24302.53665 -25790.02002 [148,] -12391.20456 -24302.53665 [149,] -12535.82700 -12391.20456 [150,] 30894.23177 -12535.82700 [151,] -4365.06184 30894.23177 [152,] 3668.55944 -4365.06184 [153,] 52325.71859 3668.55944 [154,] -20242.49620 52325.71859 [155,] 18886.15233 -20242.49620 [156,] -6577.67244 18886.15233 [157,] 23434.21451 -6577.67244 [158,] -37384.99266 23434.21451 [159,] 52234.20934 -37384.99266 [160,] 112967.21883 52234.20934 [161,] -42213.13655 112967.21883 [162,] 16573.63549 -42213.13655 [163,] 62741.33280 16573.63549 [164,] -10454.84274 62741.33280 [165,] -6419.53683 -10454.84274 [166,] -5975.10608 -6419.53683 [167,] 28564.89972 -5975.10608 [168,] -4274.83988 28564.89972 [169,] 124563.73369 -4274.83988 [170,] 87867.57984 124563.73369 [171,] -40331.79298 87867.57984 [172,] 60035.50787 -40331.79298 [173,] -6737.48329 60035.50787 [174,] -185309.72510 -6737.48329 [175,] 56973.10550 -185309.72510 [176,] -23523.24308 56973.10550 [177,] 9006.32656 -23523.24308 [178,] -59575.07568 9006.32656 [179,] 35100.29855 -59575.07568 [180,] 650.16261 35100.29855 [181,] -47574.47123 650.16261 [182,] -6902.73900 -47574.47123 [183,] 35165.98104 -6902.73900 [184,] 45421.03807 35165.98104 [185,] -16164.19011 45421.03807 [186,] 6957.67253 -16164.19011 [187,] -425.86320 6957.67253 [188,] 35314.93984 -425.86320 [189,] -31979.04053 35314.93984 [190,] 14582.57981 -31979.04053 [191,] 25724.78799 14582.57981 [192,] -10726.61674 25724.78799 [193,] -20476.60749 -10726.61674 [194,] -749.09312 -20476.60749 [195,] 9161.48249 -749.09312 [196,] -788.51354 9161.48249 [197,] -16986.49645 -788.51354 [198,] -4069.46484 -16986.49645 [199,] -23663.87858 -4069.46484 [200,] -412.52230 -23663.87858 [201,] -29274.38424 -412.52230 [202,] 37217.72342 -29274.38424 [203,] 20667.92260 37217.72342 [204,] 1527.10537 20667.92260 [205,] -24526.57156 1527.10537 [206,] -9138.75869 -24526.57156 [207,] 86196.74319 -9138.75869 [208,] 39595.53682 86196.74319 [209,] 12269.53655 39595.53682 [210,] 25681.30704 12269.53655 [211,] 22845.79586 25681.30704 [212,] -25833.26267 22845.79586 [213,] -8593.93189 -25833.26267 [214,] -12830.08305 -8593.93189 [215,] -10885.34145 -12830.08305 [216,] -15407.77494 -10885.34145 [217,] -8657.72993 -15407.77494 [218,] -15702.78861 -8657.72993 [219,] 18191.06023 -15702.78861 [220,] -18984.59768 18191.06023 [221,] -2582.31993 -18984.59768 [222,] -19641.65844 -2582.31993 [223,] -17920.24526 -19641.65844 [224,] -21191.51468 -17920.24526 [225,] 42681.45731 -21191.51468 [226,] 6579.13757 42681.45731 [227,] -35669.12186 6579.13757 [228,] -24682.85265 -35669.12186 [229,] -45906.04895 -24682.85265 [230,] -14663.01217 -45906.04895 [231,] -11474.73109 -14663.01217 [232,] -6021.34917 -11474.73109 [233,] -21203.83039 -6021.34917 [234,] -30584.43213 -21203.83039 [235,] -44172.96623 -30584.43213 [236,] 72802.30483 -44172.96623 [237,] -27982.26411 72802.30483 [238,] -13109.05887 -27982.26411 [239,] -29303.88693 -13109.05887 [240,] 15857.22161 -29303.88693 [241,] 37017.82138 15857.22161 [242,] -34035.75672 37017.82138 [243,] -22836.99206 -34035.75672 [244,] -8661.37160 -22836.99206 [245,] -9075.31576 -8661.37160 [246,] -13654.29794 -9075.31576 [247,] 20468.76156 -13654.29794 [248,] 15447.71686 20468.76156 [249,] -25187.82865 15447.71686 [250,] 48465.16659 -25187.82865 [251,] -16612.67549 48465.16659 [252,] -13534.22494 -16612.67549 [253,] 28133.05847 -13534.22494 [254,] 10158.55754 28133.05847 [255,] -1029.27115 10158.55754 [256,] -40171.14001 -1029.27115 [257,] 5055.71207 -40171.14001 [258,] -1443.36333 5055.71207 [259,] -13889.62917 -1443.36333 [260,] -19972.13487 -13889.62917 [261,] -25261.41191 -19972.13487 [262,] 3176.57626 -25261.41191 [263,] -36729.87890 3176.57626 [264,] -23054.67831 -36729.87890 [265,] -8184.26550 -23054.67831 [266,] -33241.98285 -8184.26550 [267,] -19915.33529 -33241.98285 [268,] 30455.19294 -19915.33529 [269,] -34239.62063 30455.19294 [270,] 25715.66205 -34239.62063 [271,] -6725.79016 25715.66205 [272,] 17839.06982 -6725.79016 [273,] -2202.87352 17839.06982 [274,] -5231.99488 -2202.87352 [275,] -33721.94544 -5231.99488 [276,] 1211.73279 -33721.94544 [277,] -12158.15154 1211.73279 [278,] -3796.20644 -12158.15154 [279,] 12081.95025 -3796.20644 [280,] -12243.11968 12081.95025 [281,] -3602.58689 -12243.11968 [282,] 11264.26255 -3602.58689 [283,] 7191.44346 11264.26255 [284,] 1935.02952 7191.44346 [285,] -30637.07695 1935.02952 [286,] 348.06676 -30637.07695 [287,] -22442.45967 348.06676 [288,] 2482.32017 -22442.45967 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -42318.42267 9459.55117 2 -6191.83910 -42318.42267 3 -51181.44191 -6191.83910 4 -1424.77691 -51181.44191 5 -34217.21900 -1424.77691 6 142113.49876 -34217.21900 7 -23311.33777 142113.49876 8 37933.80275 -23311.33777 9 -5452.25711 37933.80275 10 5202.08209 -5452.25711 11 40310.65109 5202.08209 12 -45502.31467 40310.65109 13 -7345.85894 -45502.31467 14 23141.98373 -7345.85894 15 12204.28474 23141.98373 16 8093.56625 12204.28474 17 -5104.05047 8093.56625 18 39982.50566 -5104.05047 19 -52177.46897 39982.50566 20 -51444.03012 -52177.46897 21 -17146.35953 -51444.03012 22 158733.95290 -17146.35953 23 12478.74058 158733.95290 24 -83418.08288 12478.74058 25 -67848.80908 -83418.08288 26 -10648.99555 -67848.80908 27 30320.34322 -10648.99555 28 2142.80883 30320.34322 29 94.76229 2142.80883 30 -1782.88705 94.76229 31 -18749.30410 -1782.88705 32 2315.48642 -18749.30410 33 4313.23727 2315.48642 34 -50167.40189 4313.23727 35 59633.04958 -50167.40189 36 18293.18762 59633.04958 37 41379.84449 18293.18762 38 119897.55241 41379.84449 39 -8067.82469 119897.55241 40 14284.84145 -8067.82469 41 17132.88928 14284.84145 42 -5281.92931 17132.88928 43 52389.55026 -5281.92931 44 53307.82657 52389.55026 45 -15071.18073 53307.82657 46 44956.57391 -15071.18073 47 -9709.61554 44956.57391 48 -16478.38248 -9709.61554 49 -33812.38589 -16478.38248 50 161619.86137 -33812.38589 51 -17583.46955 161619.86137 52 -36795.71870 -17583.46955 53 9005.58453 -36795.71870 54 -47498.83271 9005.58453 55 -10372.69665 -47498.83271 56 -9061.65831 -10372.69665 57 110528.14257 -9061.65831 58 -18128.10104 110528.14257 59 15028.86516 -18128.10104 60 34781.01161 15028.86516 61 -19416.58746 34781.01161 62 137029.71115 -19416.58746 63 20444.19031 137029.71115 64 17282.06624 20444.19031 65 -1654.78295 17282.06624 66 27803.92871 -1654.78295 67 32040.02889 27803.92871 68 2840.02156 32040.02889 69 18379.84862 2840.02156 70 38226.85155 18379.84862 71 46995.24938 38226.85155 72 11432.24145 46995.24938 73 -36963.17969 11432.24145 74 -5142.91568 -36963.17969 75 -43013.30532 -5142.91568 76 -101950.60841 -43013.30532 77 -18581.69736 -101950.60841 78 -15817.45803 -18581.69736 79 5321.68653 -15817.45803 80 -8038.56928 5321.68653 81 -78148.44311 -8038.56928 82 33610.05983 -78148.44311 83 41360.29097 33610.05983 84 -62022.46549 41360.29097 85 -69616.45843 -62022.46549 86 -15525.92548 -69616.45843 87 -11069.61226 -15525.92548 88 88434.16372 -11069.61226 89 35767.92906 88434.16372 90 -30536.89723 35767.92906 91 -25791.94407 -30536.89723 92 31851.91723 -25791.94407 93 -48856.82370 31851.91723 94 -559.69883 -48856.82370 95 120669.73879 -559.69883 96 34782.44303 120669.73879 97 -135255.32333 34782.44303 98 37721.16333 -135255.32333 99 -72746.88468 37721.16333 100 23682.28990 -72746.88468 101 15267.17887 23682.28990 102 67671.59079 15267.17887 103 -25564.35686 67671.59079 104 -27624.55868 -25564.35686 105 -49609.77672 -27624.55868 106 -12128.11343 -49609.77672 107 -18108.62133 -12128.11343 108 42177.36240 -18108.62133 109 -11030.87558 42177.36240 110 -25917.17252 -11030.87558 111 14421.25911 -25917.17252 112 -36734.85479 14421.25911 113 5121.06623 -36734.85479 114 -31452.87244 5121.06623 115 -25637.27442 -31452.87244 116 -43405.93363 -25637.27442 117 98679.33854 -43405.93363 118 -61660.48481 98679.33854 119 -11554.75924 -61660.48481 120 3258.48559 -11554.75924 121 -10691.33122 3258.48559 122 26243.54396 -10691.33122 123 -116196.47094 26243.54396 124 46827.37628 -116196.47094 125 10781.62257 46827.37628 126 15191.12727 10781.62257 127 -43226.34386 15191.12727 128 -37029.30306 -43226.34386 129 -7756.67525 -37029.30306 130 5749.94481 -7756.67525 131 -65684.64897 5749.94481 132 -8359.16811 -65684.64897 133 -21968.92476 -8359.16811 134 -43115.36233 -21968.92476 135 79667.12025 -43115.36233 136 1435.49405 79667.12025 137 44157.36125 1435.49405 138 -45807.76922 44157.36125 139 -23525.94368 -45807.76922 140 -11058.81079 -23525.94368 141 -10620.19957 -11058.81079 142 30749.35717 -10620.19957 143 29772.53703 30749.35717 144 -16845.74526 29772.53703 145 -18221.37351 -16845.74526 146 -25790.02002 -18221.37351 147 -24302.53665 -25790.02002 148 -12391.20456 -24302.53665 149 -12535.82700 -12391.20456 150 30894.23177 -12535.82700 151 -4365.06184 30894.23177 152 3668.55944 -4365.06184 153 52325.71859 3668.55944 154 -20242.49620 52325.71859 155 18886.15233 -20242.49620 156 -6577.67244 18886.15233 157 23434.21451 -6577.67244 158 -37384.99266 23434.21451 159 52234.20934 -37384.99266 160 112967.21883 52234.20934 161 -42213.13655 112967.21883 162 16573.63549 -42213.13655 163 62741.33280 16573.63549 164 -10454.84274 62741.33280 165 -6419.53683 -10454.84274 166 -5975.10608 -6419.53683 167 28564.89972 -5975.10608 168 -4274.83988 28564.89972 169 124563.73369 -4274.83988 170 87867.57984 124563.73369 171 -40331.79298 87867.57984 172 60035.50787 -40331.79298 173 -6737.48329 60035.50787 174 -185309.72510 -6737.48329 175 56973.10550 -185309.72510 176 -23523.24308 56973.10550 177 9006.32656 -23523.24308 178 -59575.07568 9006.32656 179 35100.29855 -59575.07568 180 650.16261 35100.29855 181 -47574.47123 650.16261 182 -6902.73900 -47574.47123 183 35165.98104 -6902.73900 184 45421.03807 35165.98104 185 -16164.19011 45421.03807 186 6957.67253 -16164.19011 187 -425.86320 6957.67253 188 35314.93984 -425.86320 189 -31979.04053 35314.93984 190 14582.57981 -31979.04053 191 25724.78799 14582.57981 192 -10726.61674 25724.78799 193 -20476.60749 -10726.61674 194 -749.09312 -20476.60749 195 9161.48249 -749.09312 196 -788.51354 9161.48249 197 -16986.49645 -788.51354 198 -4069.46484 -16986.49645 199 -23663.87858 -4069.46484 200 -412.52230 -23663.87858 201 -29274.38424 -412.52230 202 37217.72342 -29274.38424 203 20667.92260 37217.72342 204 1527.10537 20667.92260 205 -24526.57156 1527.10537 206 -9138.75869 -24526.57156 207 86196.74319 -9138.75869 208 39595.53682 86196.74319 209 12269.53655 39595.53682 210 25681.30704 12269.53655 211 22845.79586 25681.30704 212 -25833.26267 22845.79586 213 -8593.93189 -25833.26267 214 -12830.08305 -8593.93189 215 -10885.34145 -12830.08305 216 -15407.77494 -10885.34145 217 -8657.72993 -15407.77494 218 -15702.78861 -8657.72993 219 18191.06023 -15702.78861 220 -18984.59768 18191.06023 221 -2582.31993 -18984.59768 222 -19641.65844 -2582.31993 223 -17920.24526 -19641.65844 224 -21191.51468 -17920.24526 225 42681.45731 -21191.51468 226 6579.13757 42681.45731 227 -35669.12186 6579.13757 228 -24682.85265 -35669.12186 229 -45906.04895 -24682.85265 230 -14663.01217 -45906.04895 231 -11474.73109 -14663.01217 232 -6021.34917 -11474.73109 233 -21203.83039 -6021.34917 234 -30584.43213 -21203.83039 235 -44172.96623 -30584.43213 236 72802.30483 -44172.96623 237 -27982.26411 72802.30483 238 -13109.05887 -27982.26411 239 -29303.88693 -13109.05887 240 15857.22161 -29303.88693 241 37017.82138 15857.22161 242 -34035.75672 37017.82138 243 -22836.99206 -34035.75672 244 -8661.37160 -22836.99206 245 -9075.31576 -8661.37160 246 -13654.29794 -9075.31576 247 20468.76156 -13654.29794 248 15447.71686 20468.76156 249 -25187.82865 15447.71686 250 48465.16659 -25187.82865 251 -16612.67549 48465.16659 252 -13534.22494 -16612.67549 253 28133.05847 -13534.22494 254 10158.55754 28133.05847 255 -1029.27115 10158.55754 256 -40171.14001 -1029.27115 257 5055.71207 -40171.14001 258 -1443.36333 5055.71207 259 -13889.62917 -1443.36333 260 -19972.13487 -13889.62917 261 -25261.41191 -19972.13487 262 3176.57626 -25261.41191 263 -36729.87890 3176.57626 264 -23054.67831 -36729.87890 265 -8184.26550 -23054.67831 266 -33241.98285 -8184.26550 267 -19915.33529 -33241.98285 268 30455.19294 -19915.33529 269 -34239.62063 30455.19294 270 25715.66205 -34239.62063 271 -6725.79016 25715.66205 272 17839.06982 -6725.79016 273 -2202.87352 17839.06982 274 -5231.99488 -2202.87352 275 -33721.94544 -5231.99488 276 1211.73279 -33721.94544 277 -12158.15154 1211.73279 278 -3796.20644 -12158.15154 279 12081.95025 -3796.20644 280 -12243.11968 12081.95025 281 -3602.58689 -12243.11968 282 11264.26255 -3602.58689 283 7191.44346 11264.26255 284 1935.02952 7191.44346 285 -30637.07695 1935.02952 286 348.06676 -30637.07695 287 -22442.45967 348.06676 288 2482.32017 -22442.45967 > 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/7opul1324127192.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/81c781324127192.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/9ycpj1324127192.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/10pkhv1324127192.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/11phjv1324127192.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/12rbfa1324127192.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/13dwi51324127192.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/14rpfg1324127192.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/15ea2g1324127192.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/16m9pl1324127192.tab") + } > > try(system("convert tmp/17m4g1324127192.ps tmp/17m4g1324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/2kdvx1324127192.ps tmp/2kdvx1324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/3650n1324127192.ps tmp/3650n1324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/469sk1324127192.ps tmp/469sk1324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/5i1h81324127192.ps tmp/5i1h81324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/6su7a1324127192.ps tmp/6su7a1324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/7opul1324127192.ps tmp/7opul1324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/81c781324127192.ps tmp/81c781324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/9ycpj1324127192.ps tmp/9ycpj1324127192.png",intern=TRUE)) character(0) > try(system("convert tmp/10pkhv1324127192.ps tmp/10pkhv1324127192.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 7.526 0.637 8.255