R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1) + ,dim=c(5 + ,154) + ,dimnames=list(c('Correctanalysis' + ,'T40Treatment' + ,'T40NoTreatment' + ,'T20Treatment' + ,'T20NoTreatment') + ,1:154)) > y <- array(NA,dim=c(5,154),dimnames=list(c('Correctanalysis','T40Treatment','T40NoTreatment','T20Treatment','T20NoTreatment'),1:154)) > 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' > 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 Attaching package: 'zoo' The following object(s) are masked from 'package:base': as.Date, as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Correctanalysis T40Treatment T40NoTreatment T20Treatment T20NoTreatment 1 0 1 0 0 0 2 0 0 1 0 0 3 0 0 1 0 0 4 0 0 1 0 0 5 0 0 1 0 0 6 0 0 1 0 0 7 0 0 1 0 0 8 0 1 0 0 0 9 0 0 1 0 0 10 0 0 1 0 0 11 0 1 0 0 0 12 0 0 1 0 0 13 0 0 1 0 0 14 0 1 0 0 0 15 0 0 1 0 0 16 0 1 0 0 0 17 1 1 0 0 0 18 0 1 0 0 0 19 0 0 1 0 0 20 1 1 0 0 0 21 0 0 1 0 0 22 0 0 1 0 0 23 0 0 1 0 0 24 0 0 1 0 0 25 0 1 0 0 0 26 0 0 1 0 0 27 0 0 1 0 0 28 0 0 1 0 0 29 0 0 1 0 0 30 0 0 1 0 0 31 0 0 1 0 0 32 0 0 1 0 0 33 0 0 1 0 0 34 0 1 0 0 0 35 0 0 1 0 0 36 0 0 1 0 0 37 0 1 0 0 0 38 0 0 1 0 0 39 0 0 1 0 0 40 0 1 0 0 0 41 1 0 1 0 0 42 0 0 1 0 0 43 0 0 1 0 0 44 0 1 0 0 0 45 0 0 1 0 0 46 0 0 1 0 0 47 0 0 1 0 0 48 0 0 1 0 0 49 0 0 1 0 0 50 0 0 1 0 0 51 0 1 0 0 0 52 1 1 0 0 0 53 0 0 1 0 0 54 1 0 1 0 0 55 0 0 1 0 0 56 0 1 0 0 0 57 0 0 1 0 0 58 0 0 1 0 0 59 0 0 1 0 0 60 1 1 0 0 0 61 0 1 0 0 0 62 0 0 1 0 0 63 0 0 1 0 0 64 0 1 0 0 0 65 0 0 1 0 0 66 0 0 1 0 0 67 1 1 0 0 0 68 0 0 1 0 0 69 0 0 1 0 0 70 0 0 1 0 0 71 0 0 1 0 0 72 0 0 1 0 0 73 0 0 1 0 0 74 0 0 1 0 0 75 0 0 1 0 0 76 0 1 0 0 0 77 0 0 1 0 0 78 0 0 1 0 0 79 1 1 0 0 0 80 0 1 0 0 0 81 0 0 1 0 0 82 0 0 1 0 0 83 0 0 1 0 0 84 1 0 1 0 0 85 0 0 1 0 0 86 0 0 1 0 0 87 0 0 0 0 1 88 0 0 0 1 0 89 0 0 0 0 1 90 0 0 0 0 1 91 0 0 0 0 1 92 0 0 0 1 0 93 0 0 0 0 1 94 0 0 0 0 1 95 0 0 0 1 0 96 0 0 0 0 1 97 0 0 0 1 0 98 0 0 0 0 1 99 0 0 0 0 1 100 0 0 0 0 1 101 0 0 0 0 1 102 0 0 0 0 1 103 0 0 0 0 1 104 0 0 0 0 1 105 0 0 0 1 0 106 0 0 0 0 1 107 0 0 0 0 1 108 0 0 0 1 0 109 0 0 0 0 1 110 0 0 0 0 1 111 0 0 0 1 0 112 0 0 0 1 0 113 0 0 0 0 1 114 0 0 0 1 0 115 0 0 0 0 1 116 0 0 0 0 1 117 0 0 0 0 1 118 0 0 0 0 1 119 0 0 0 0 1 120 0 0 0 0 1 121 0 0 0 0 1 122 0 0 0 0 1 123 0 0 0 1 0 124 0 0 0 0 1 125 0 0 0 0 1 126 0 0 0 1 0 127 0 0 0 0 1 128 0 0 0 0 1 129 0 0 0 0 1 130 0 0 0 0 1 131 0 0 0 0 1 132 0 0 0 0 1 133 0 0 0 0 1 134 0 0 0 0 1 135 0 0 0 0 1 136 0 0 0 0 1 137 0 0 0 0 1 138 0 0 0 1 0 139 0 0 0 1 0 140 0 0 0 0 1 141 1 0 0 0 1 142 0 0 0 1 0 143 0 0 0 0 1 144 0 0 0 0 1 145 0 0 0 0 1 146 0 0 0 1 0 147 0 0 0 1 0 148 0 0 0 1 0 149 0 0 0 0 1 150 0 0 0 0 1 151 0 0 0 0 1 152 1 0 0 0 1 153 1 0 0 0 1 154 0 0 0 0 1 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) T40Treatment T40NoTreatment T20Treatment T20NoTreatment 0.05882 0.20205 -0.01120 -0.05882 NA > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -0.26087 -0.05882 -0.04762 -0.04762 0.95238 Coefficients: (1 not defined because of singularities) Estimate Std. Error t value Pr(>|t|) (Intercept) 0.05882 0.03636 1.618 0.10784 T40Treatment 0.20205 0.06522 3.098 0.00233 ** T40NoTreatment -0.01120 0.04892 -0.229 0.81913 T20Treatment -0.05882 0.07273 -0.809 0.41989 T20NoTreatment NA NA NA NA --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.2597 on 150 degrees of freedom Multiple R-squared: 0.08581, Adjusted R-squared: 0.06753 F-statistic: 4.693 on 3 and 150 DF, p-value: 0.003668 > 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.000000e+00 0.0000000000 1.000000000 [2,] 0.000000e+00 0.0000000000 1.000000000 [3,] 0.000000e+00 0.0000000000 1.000000000 [4,] 0.000000e+00 0.0000000000 1.000000000 [5,] 0.000000e+00 0.0000000000 1.000000000 [6,] 0.000000e+00 0.0000000000 1.000000000 [7,] 0.000000e+00 0.0000000000 1.000000000 [8,] 0.000000e+00 0.0000000000 1.000000000 [9,] 0.000000e+00 0.0000000000 1.000000000 [10,] 5.460373e-01 0.9079254577 0.453962729 [11,] 4.963416e-01 0.9926831425 0.503658429 [12,] 4.139806e-01 0.8279612109 0.586019395 [13,] 8.870002e-01 0.2259995282 0.112999764 [14,] 8.479006e-01 0.3041988134 0.152099407 [15,] 8.016855e-01 0.3966290148 0.198314507 [16,] 7.488962e-01 0.5022076879 0.251103844 [17,] 6.905213e-01 0.6189574599 0.309478730 [18,] 6.774372e-01 0.6451256767 0.322562838 [19,] 6.152486e-01 0.7695027561 0.384751378 [20,] 5.508782e-01 0.8982436310 0.449121816 [21,] 4.860278e-01 0.9720556257 0.513972187 [22,] 4.223712e-01 0.8447423707 0.577628815 [23,] 3.614320e-01 0.7228639564 0.638568022 [24,] 3.044870e-01 0.6089739857 0.695513007 [25,] 2.525025e-01 0.5050049405 0.747497530 [26,] 2.061064e-01 0.4122127514 0.793893624 [27,] 1.974220e-01 0.3948439166 0.802578042 [28,] 1.585385e-01 0.3170770554 0.841461472 [29,] 1.253620e-01 0.2507239318 0.874638034 [30,] 1.181377e-01 0.2362754729 0.881862264 [31,] 9.189019e-02 0.1837803769 0.908109812 [32,] 7.041345e-02 0.1408268937 0.929586553 [33,] 6.582486e-02 0.1316497143 0.934175143 [34,] 6.404319e-01 0.7191361933 0.359568097 [35,] 5.898742e-01 0.8202515883 0.410125794 [36,] 5.381094e-01 0.9237811615 0.461890581 [37,] 5.294548e-01 0.9410904537 0.470545227 [38,] 4.777825e-01 0.9555650162 0.522217492 [39,] 4.267880e-01 0.8535759405 0.573212030 [40,] 3.772842e-01 0.7545683218 0.622715839 [41,] 3.299994e-01 0.6599988078 0.670000596 [42,] 2.855478e-01 0.5710956720 0.714452164 [43,] 2.444092e-01 0.4888184093 0.755590795 [44,] 2.444044e-01 0.4888087262 0.755595637 [45,] 6.081938e-01 0.7836123574 0.391806179 [46,] 5.609416e-01 0.8781167826 0.439058391 [47,] 9.458109e-01 0.1083782508 0.054189125 [48,] 9.314299e-01 0.1371402158 0.068570108 [49,] 9.357275e-01 0.1285450857 0.064272543 [50,] 9.195010e-01 0.1609980524 0.080499026 [51,] 9.004144e-01 0.1991712375 0.099585619 [52,] 8.782854e-01 0.2434292108 0.121714605 [53,] 9.717460e-01 0.0565079497 0.028253975 [54,] 9.739488e-01 0.0521024381 0.026051219 [55,] 9.660992e-01 0.0678016227 0.033900811 [56,] 9.564334e-01 0.0871332600 0.043566630 [57,] 9.643484e-01 0.0713032612 0.035651631 [58,] 9.543588e-01 0.0912823406 0.045641170 [59,] 9.422837e-01 0.1154326135 0.057716307 [60,] 9.870464e-01 0.0259072491 0.012953625 [61,] 9.827422e-01 0.0345156764 0.017257838 [62,] 9.773028e-01 0.0453944091 0.022697205 [63,] 9.705327e-01 0.0589345924 0.029467296 [64,] 9.622374e-01 0.0755251010 0.037762550 [65,] 9.522386e-01 0.0955227422 0.047761371 [66,] 9.403951e-01 0.1192097499 0.059604875 [67,] 9.266336e-01 0.1467327886 0.073366394 [68,] 9.109940e-01 0.1780120054 0.089006003 [69,] 9.244874e-01 0.1510251144 0.075512557 [70,] 9.092297e-01 0.1815405809 0.090770290 [71,] 8.929046e-01 0.2141907928 0.107095396 [72,] 9.794678e-01 0.0410643654 0.020532183 [73,] 9.755789e-01 0.0488422403 0.024421120 [74,] 9.705543e-01 0.0588914438 0.029445722 [75,] 9.661769e-01 0.0676462162 0.033823108 [76,] 9.644885e-01 0.0710230593 0.035511530 [77,] 9.989950e-01 0.0020099952 0.001004998 [78,] 9.985015e-01 0.0029970715 0.001498536 [79,] 9.977937e-01 0.0044126761 0.002206338 [80,] 9.968286e-01 0.0063428977 0.003171449 [81,] 9.954297e-01 0.0091406015 0.004570301 [82,] 9.935924e-01 0.0128152623 0.006407631 [83,] 9.911274e-01 0.0177452978 0.008872649 [84,] 9.878651e-01 0.0242697125 0.012134856 [85,] 9.833673e-01 0.0332653144 0.016632657 [86,] 9.778044e-01 0.0443912199 0.022195610 [87,] 9.707402e-01 0.0585196411 0.029259821 [88,] 9.613648e-01 0.0772704709 0.038635235 [89,] 9.502883e-01 0.0994233399 0.049711670 [90,] 9.359732e-01 0.1280536068 0.064026803 [91,] 9.195688e-01 0.1608624456 0.080431223 [92,] 9.001474e-01 0.1997052334 0.099852617 [93,] 8.774776e-01 0.2450447859 0.122522393 [94,] 8.513901e-01 0.2972198571 0.148609929 [95,] 8.217957e-01 0.3564085369 0.178204268 [96,] 7.887021e-01 0.4225957530 0.211297877 [97,] 7.522260e-01 0.4955480555 0.247774028 [98,] 7.094025e-01 0.5811950535 0.290597527 [99,] 6.665885e-01 0.6668230496 0.333411525 [100,] 6.214435e-01 0.7571130919 0.378556546 [101,] 5.704418e-01 0.8591164360 0.429558218 [102,] 5.221167e-01 0.9557666298 0.477883315 [103,] 4.734227e-01 0.9468453716 0.526577314 [104,] 4.204994e-01 0.8409988091 0.579500595 [105,] 3.686870e-01 0.7373739404 0.631313030 [106,] 3.230497e-01 0.6460993255 0.676950337 [107,] 2.757261e-01 0.5514521100 0.724273945 [108,] 2.356126e-01 0.4712251748 0.764387413 [109,] 1.989077e-01 0.3978154340 0.801092283 [110,] 1.658551e-01 0.3317101156 0.834144942 [111,] 1.365660e-01 0.2731319066 0.863434047 [112,] 1.110273e-01 0.2220546643 0.888972668 [113,] 8.911624e-02 0.1782324710 0.910883764 [114,] 7.061924e-02 0.1412384836 0.929380758 [115,] 5.525476e-02 0.1105095173 0.944745241 [116,] 4.070964e-02 0.0814192878 0.959290356 [117,] 3.085910e-02 0.0617182045 0.969140898 [118,] 2.309510e-02 0.0461901983 0.976904901 [119,] 1.598001e-02 0.0319600107 0.984019995 [120,] 1.156207e-02 0.0231241397 0.988437930 [121,] 8.263645e-03 0.0165272899 0.991736355 [122,] 5.840256e-03 0.0116805129 0.994159744 [123,] 4.087300e-03 0.0081746007 0.995912700 [124,] 2.838164e-03 0.0056763284 0.997161836 [125,] 1.960654e-03 0.0039213083 0.998039346 [126,] 1.352464e-03 0.0027049273 0.998647536 [127,] 9.363005e-04 0.0018726010 0.999063699 [128,] 6.551458e-04 0.0013102915 0.999344854 [129,] 4.679752e-04 0.0009359504 0.999532025 [130,] 3.461699e-04 0.0006923397 0.999653830 [131,] 1.720733e-04 0.0003441466 0.999827927 [132,] 8.116632e-05 0.0001623326 0.999918834 [133,] 5.892338e-05 0.0001178468 0.999941077 [134,] 2.335148e-03 0.0046702951 0.997664852 [135,] 1.121207e-03 0.0022424146 0.998878793 [136,] 6.923688e-04 0.0013847377 0.999307631 [137,] 4.421800e-04 0.0008843601 0.999557820 [138,] 3.060936e-04 0.0006121872 0.999693906 [139,] 9.797086e-05 0.0001959417 0.999902029 > postscript(file="/var/wessaorg/rcomp/tmp/1l9f51356032183.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/2at5r1356032183.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/3ib7o1356032183.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/42h531356032183.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/5aqdk1356032183.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 = 154 Frequency = 1 1 2 3 4 5 -2.608696e-01 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 6 7 8 9 10 -4.761905e-02 -4.761905e-02 -2.608696e-01 -4.761905e-02 -4.761905e-02 11 12 13 14 15 -2.608696e-01 -4.761905e-02 -4.761905e-02 -2.608696e-01 -4.761905e-02 16 17 18 19 20 -2.608696e-01 7.391304e-01 -2.608696e-01 -4.761905e-02 7.391304e-01 21 22 23 24 25 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 -2.608696e-01 26 27 28 29 30 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 31 32 33 34 35 -4.761905e-02 -4.761905e-02 -4.761905e-02 -2.608696e-01 -4.761905e-02 36 37 38 39 40 -4.761905e-02 -2.608696e-01 -4.761905e-02 -4.761905e-02 -2.608696e-01 41 42 43 44 45 9.523810e-01 -4.761905e-02 -4.761905e-02 -2.608696e-01 -4.761905e-02 46 47 48 49 50 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 51 52 53 54 55 -2.608696e-01 7.391304e-01 -4.761905e-02 9.523810e-01 -4.761905e-02 56 57 58 59 60 -2.608696e-01 -4.761905e-02 -4.761905e-02 -4.761905e-02 7.391304e-01 61 62 63 64 65 -2.608696e-01 -4.761905e-02 -4.761905e-02 -2.608696e-01 -4.761905e-02 66 67 68 69 70 -4.761905e-02 7.391304e-01 -4.761905e-02 -4.761905e-02 -4.761905e-02 71 72 73 74 75 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 -4.761905e-02 76 77 78 79 80 -2.608696e-01 -4.761905e-02 -4.761905e-02 7.391304e-01 -2.608696e-01 81 82 83 84 85 -4.761905e-02 -4.761905e-02 -4.761905e-02 9.523810e-01 -4.761905e-02 86 87 88 89 90 -4.761905e-02 -5.882353e-02 -2.297153e-18 -5.882353e-02 -5.882353e-02 91 92 93 94 95 -5.882353e-02 -2.297153e-18 -5.882353e-02 -5.882353e-02 -2.297153e-18 96 97 98 99 100 -5.882353e-02 -2.297153e-18 -5.882353e-02 -5.882353e-02 -5.882353e-02 101 102 103 104 105 -5.882353e-02 -5.882353e-02 -5.882353e-02 -5.882353e-02 -2.297153e-18 106 107 108 109 110 -5.882353e-02 -5.882353e-02 -2.297153e-18 -5.882353e-02 -5.882353e-02 111 112 113 114 115 -2.297153e-18 -2.297153e-18 -5.882353e-02 -2.297153e-18 -5.882353e-02 116 117 118 119 120 -5.882353e-02 -5.882353e-02 -5.882353e-02 -5.882353e-02 -5.882353e-02 121 122 123 124 125 -5.882353e-02 -5.882353e-02 -2.297153e-18 -5.882353e-02 -5.882353e-02 126 127 128 129 130 -2.297153e-18 -5.882353e-02 -5.882353e-02 -5.882353e-02 -5.882353e-02 131 132 133 134 135 -5.882353e-02 -5.882353e-02 -5.882353e-02 -5.882353e-02 -5.882353e-02 136 137 138 139 140 -5.882353e-02 -5.882353e-02 -2.297153e-18 -2.297153e-18 -5.882353e-02 141 142 143 144 145 9.411765e-01 -2.297153e-18 -5.882353e-02 -5.882353e-02 -5.882353e-02 146 147 148 149 150 -2.297153e-18 -2.297153e-18 -2.297153e-18 -5.882353e-02 -5.882353e-02 151 152 153 154 -5.882353e-02 9.411765e-01 9.411765e-01 -5.882353e-02 > postscript(file="/var/wessaorg/rcomp/tmp/6ifm91356032183.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 = 154 Frequency = 1 lag(myerror, k = 1) myerror 0 -2.608696e-01 NA 1 -4.761905e-02 -2.608696e-01 2 -4.761905e-02 -4.761905e-02 3 -4.761905e-02 -4.761905e-02 4 -4.761905e-02 -4.761905e-02 5 -4.761905e-02 -4.761905e-02 6 -4.761905e-02 -4.761905e-02 7 -2.608696e-01 -4.761905e-02 8 -4.761905e-02 -2.608696e-01 9 -4.761905e-02 -4.761905e-02 10 -2.608696e-01 -4.761905e-02 11 -4.761905e-02 -2.608696e-01 12 -4.761905e-02 -4.761905e-02 13 -2.608696e-01 -4.761905e-02 14 -4.761905e-02 -2.608696e-01 15 -2.608696e-01 -4.761905e-02 16 7.391304e-01 -2.608696e-01 17 -2.608696e-01 7.391304e-01 18 -4.761905e-02 -2.608696e-01 19 7.391304e-01 -4.761905e-02 20 -4.761905e-02 7.391304e-01 21 -4.761905e-02 -4.761905e-02 22 -4.761905e-02 -4.761905e-02 23 -4.761905e-02 -4.761905e-02 24 -2.608696e-01 -4.761905e-02 25 -4.761905e-02 -2.608696e-01 26 -4.761905e-02 -4.761905e-02 27 -4.761905e-02 -4.761905e-02 28 -4.761905e-02 -4.761905e-02 29 -4.761905e-02 -4.761905e-02 30 -4.761905e-02 -4.761905e-02 31 -4.761905e-02 -4.761905e-02 32 -4.761905e-02 -4.761905e-02 33 -2.608696e-01 -4.761905e-02 34 -4.761905e-02 -2.608696e-01 35 -4.761905e-02 -4.761905e-02 36 -2.608696e-01 -4.761905e-02 37 -4.761905e-02 -2.608696e-01 38 -4.761905e-02 -4.761905e-02 39 -2.608696e-01 -4.761905e-02 40 9.523810e-01 -2.608696e-01 41 -4.761905e-02 9.523810e-01 42 -4.761905e-02 -4.761905e-02 43 -2.608696e-01 -4.761905e-02 44 -4.761905e-02 -2.608696e-01 45 -4.761905e-02 -4.761905e-02 46 -4.761905e-02 -4.761905e-02 47 -4.761905e-02 -4.761905e-02 48 -4.761905e-02 -4.761905e-02 49 -4.761905e-02 -4.761905e-02 50 -2.608696e-01 -4.761905e-02 51 7.391304e-01 -2.608696e-01 52 -4.761905e-02 7.391304e-01 53 9.523810e-01 -4.761905e-02 54 -4.761905e-02 9.523810e-01 55 -2.608696e-01 -4.761905e-02 56 -4.761905e-02 -2.608696e-01 57 -4.761905e-02 -4.761905e-02 58 -4.761905e-02 -4.761905e-02 59 7.391304e-01 -4.761905e-02 60 -2.608696e-01 7.391304e-01 61 -4.761905e-02 -2.608696e-01 62 -4.761905e-02 -4.761905e-02 63 -2.608696e-01 -4.761905e-02 64 -4.761905e-02 -2.608696e-01 65 -4.761905e-02 -4.761905e-02 66 7.391304e-01 -4.761905e-02 67 -4.761905e-02 7.391304e-01 68 -4.761905e-02 -4.761905e-02 69 -4.761905e-02 -4.761905e-02 70 -4.761905e-02 -4.761905e-02 71 -4.761905e-02 -4.761905e-02 72 -4.761905e-02 -4.761905e-02 73 -4.761905e-02 -4.761905e-02 74 -4.761905e-02 -4.761905e-02 75 -2.608696e-01 -4.761905e-02 76 -4.761905e-02 -2.608696e-01 77 -4.761905e-02 -4.761905e-02 78 7.391304e-01 -4.761905e-02 79 -2.608696e-01 7.391304e-01 80 -4.761905e-02 -2.608696e-01 81 -4.761905e-02 -4.761905e-02 82 -4.761905e-02 -4.761905e-02 83 9.523810e-01 -4.761905e-02 84 -4.761905e-02 9.523810e-01 85 -4.761905e-02 -4.761905e-02 86 -5.882353e-02 -4.761905e-02 87 -2.297153e-18 -5.882353e-02 88 -5.882353e-02 -2.297153e-18 89 -5.882353e-02 -5.882353e-02 90 -5.882353e-02 -5.882353e-02 91 -2.297153e-18 -5.882353e-02 92 -5.882353e-02 -2.297153e-18 93 -5.882353e-02 -5.882353e-02 94 -2.297153e-18 -5.882353e-02 95 -5.882353e-02 -2.297153e-18 96 -2.297153e-18 -5.882353e-02 97 -5.882353e-02 -2.297153e-18 98 -5.882353e-02 -5.882353e-02 99 -5.882353e-02 -5.882353e-02 100 -5.882353e-02 -5.882353e-02 101 -5.882353e-02 -5.882353e-02 102 -5.882353e-02 -5.882353e-02 103 -5.882353e-02 -5.882353e-02 104 -2.297153e-18 -5.882353e-02 105 -5.882353e-02 -2.297153e-18 106 -5.882353e-02 -5.882353e-02 107 -2.297153e-18 -5.882353e-02 108 -5.882353e-02 -2.297153e-18 109 -5.882353e-02 -5.882353e-02 110 -2.297153e-18 -5.882353e-02 111 -2.297153e-18 -2.297153e-18 112 -5.882353e-02 -2.297153e-18 113 -2.297153e-18 -5.882353e-02 114 -5.882353e-02 -2.297153e-18 115 -5.882353e-02 -5.882353e-02 116 -5.882353e-02 -5.882353e-02 117 -5.882353e-02 -5.882353e-02 118 -5.882353e-02 -5.882353e-02 119 -5.882353e-02 -5.882353e-02 120 -5.882353e-02 -5.882353e-02 121 -5.882353e-02 -5.882353e-02 122 -2.297153e-18 -5.882353e-02 123 -5.882353e-02 -2.297153e-18 124 -5.882353e-02 -5.882353e-02 125 -2.297153e-18 -5.882353e-02 126 -5.882353e-02 -2.297153e-18 127 -5.882353e-02 -5.882353e-02 128 -5.882353e-02 -5.882353e-02 129 -5.882353e-02 -5.882353e-02 130 -5.882353e-02 -5.882353e-02 131 -5.882353e-02 -5.882353e-02 132 -5.882353e-02 -5.882353e-02 133 -5.882353e-02 -5.882353e-02 134 -5.882353e-02 -5.882353e-02 135 -5.882353e-02 -5.882353e-02 136 -5.882353e-02 -5.882353e-02 137 -2.297153e-18 -5.882353e-02 138 -2.297153e-18 -2.297153e-18 139 -5.882353e-02 -2.297153e-18 140 9.411765e-01 -5.882353e-02 141 -2.297153e-18 9.411765e-01 142 -5.882353e-02 -2.297153e-18 143 -5.882353e-02 -5.882353e-02 144 -5.882353e-02 -5.882353e-02 145 -2.297153e-18 -5.882353e-02 146 -2.297153e-18 -2.297153e-18 147 -2.297153e-18 -2.297153e-18 148 -5.882353e-02 -2.297153e-18 149 -5.882353e-02 -5.882353e-02 150 -5.882353e-02 -5.882353e-02 151 9.411765e-01 -5.882353e-02 152 9.411765e-01 9.411765e-01 153 -5.882353e-02 9.411765e-01 154 NA -5.882353e-02 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -4.761905e-02 -2.608696e-01 [2,] -4.761905e-02 -4.761905e-02 [3,] -4.761905e-02 -4.761905e-02 [4,] -4.761905e-02 -4.761905e-02 [5,] -4.761905e-02 -4.761905e-02 [6,] -4.761905e-02 -4.761905e-02 [7,] -2.608696e-01 -4.761905e-02 [8,] -4.761905e-02 -2.608696e-01 [9,] -4.761905e-02 -4.761905e-02 [10,] -2.608696e-01 -4.761905e-02 [11,] -4.761905e-02 -2.608696e-01 [12,] -4.761905e-02 -4.761905e-02 [13,] -2.608696e-01 -4.761905e-02 [14,] -4.761905e-02 -2.608696e-01 [15,] -2.608696e-01 -4.761905e-02 [16,] 7.391304e-01 -2.608696e-01 [17,] -2.608696e-01 7.391304e-01 [18,] -4.761905e-02 -2.608696e-01 [19,] 7.391304e-01 -4.761905e-02 [20,] -4.761905e-02 7.391304e-01 [21,] -4.761905e-02 -4.761905e-02 [22,] -4.761905e-02 -4.761905e-02 [23,] -4.761905e-02 -4.761905e-02 [24,] -2.608696e-01 -4.761905e-02 [25,] -4.761905e-02 -2.608696e-01 [26,] -4.761905e-02 -4.761905e-02 [27,] -4.761905e-02 -4.761905e-02 [28,] -4.761905e-02 -4.761905e-02 [29,] -4.761905e-02 -4.761905e-02 [30,] -4.761905e-02 -4.761905e-02 [31,] -4.761905e-02 -4.761905e-02 [32,] -4.761905e-02 -4.761905e-02 [33,] -2.608696e-01 -4.761905e-02 [34,] -4.761905e-02 -2.608696e-01 [35,] -4.761905e-02 -4.761905e-02 [36,] -2.608696e-01 -4.761905e-02 [37,] -4.761905e-02 -2.608696e-01 [38,] -4.761905e-02 -4.761905e-02 [39,] -2.608696e-01 -4.761905e-02 [40,] 9.523810e-01 -2.608696e-01 [41,] -4.761905e-02 9.523810e-01 [42,] -4.761905e-02 -4.761905e-02 [43,] -2.608696e-01 -4.761905e-02 [44,] -4.761905e-02 -2.608696e-01 [45,] -4.761905e-02 -4.761905e-02 [46,] -4.761905e-02 -4.761905e-02 [47,] -4.761905e-02 -4.761905e-02 [48,] -4.761905e-02 -4.761905e-02 [49,] -4.761905e-02 -4.761905e-02 [50,] -2.608696e-01 -4.761905e-02 [51,] 7.391304e-01 -2.608696e-01 [52,] -4.761905e-02 7.391304e-01 [53,] 9.523810e-01 -4.761905e-02 [54,] -4.761905e-02 9.523810e-01 [55,] -2.608696e-01 -4.761905e-02 [56,] -4.761905e-02 -2.608696e-01 [57,] -4.761905e-02 -4.761905e-02 [58,] -4.761905e-02 -4.761905e-02 [59,] 7.391304e-01 -4.761905e-02 [60,] -2.608696e-01 7.391304e-01 [61,] -4.761905e-02 -2.608696e-01 [62,] -4.761905e-02 -4.761905e-02 [63,] -2.608696e-01 -4.761905e-02 [64,] -4.761905e-02 -2.608696e-01 [65,] -4.761905e-02 -4.761905e-02 [66,] 7.391304e-01 -4.761905e-02 [67,] -4.761905e-02 7.391304e-01 [68,] -4.761905e-02 -4.761905e-02 [69,] -4.761905e-02 -4.761905e-02 [70,] -4.761905e-02 -4.761905e-02 [71,] -4.761905e-02 -4.761905e-02 [72,] -4.761905e-02 -4.761905e-02 [73,] -4.761905e-02 -4.761905e-02 [74,] -4.761905e-02 -4.761905e-02 [75,] -2.608696e-01 -4.761905e-02 [76,] -4.761905e-02 -2.608696e-01 [77,] -4.761905e-02 -4.761905e-02 [78,] 7.391304e-01 -4.761905e-02 [79,] -2.608696e-01 7.391304e-01 [80,] -4.761905e-02 -2.608696e-01 [81,] -4.761905e-02 -4.761905e-02 [82,] -4.761905e-02 -4.761905e-02 [83,] 9.523810e-01 -4.761905e-02 [84,] -4.761905e-02 9.523810e-01 [85,] -4.761905e-02 -4.761905e-02 [86,] -5.882353e-02 -4.761905e-02 [87,] -2.297153e-18 -5.882353e-02 [88,] -5.882353e-02 -2.297153e-18 [89,] -5.882353e-02 -5.882353e-02 [90,] -5.882353e-02 -5.882353e-02 [91,] -2.297153e-18 -5.882353e-02 [92,] -5.882353e-02 -2.297153e-18 [93,] -5.882353e-02 -5.882353e-02 [94,] -2.297153e-18 -5.882353e-02 [95,] -5.882353e-02 -2.297153e-18 [96,] -2.297153e-18 -5.882353e-02 [97,] -5.882353e-02 -2.297153e-18 [98,] -5.882353e-02 -5.882353e-02 [99,] -5.882353e-02 -5.882353e-02 [100,] -5.882353e-02 -5.882353e-02 [101,] -5.882353e-02 -5.882353e-02 [102,] -5.882353e-02 -5.882353e-02 [103,] -5.882353e-02 -5.882353e-02 [104,] -2.297153e-18 -5.882353e-02 [105,] -5.882353e-02 -2.297153e-18 [106,] -5.882353e-02 -5.882353e-02 [107,] -2.297153e-18 -5.882353e-02 [108,] -5.882353e-02 -2.297153e-18 [109,] -5.882353e-02 -5.882353e-02 [110,] -2.297153e-18 -5.882353e-02 [111,] -2.297153e-18 -2.297153e-18 [112,] -5.882353e-02 -2.297153e-18 [113,] -2.297153e-18 -5.882353e-02 [114,] -5.882353e-02 -2.297153e-18 [115,] -5.882353e-02 -5.882353e-02 [116,] -5.882353e-02 -5.882353e-02 [117,] -5.882353e-02 -5.882353e-02 [118,] -5.882353e-02 -5.882353e-02 [119,] -5.882353e-02 -5.882353e-02 [120,] -5.882353e-02 -5.882353e-02 [121,] -5.882353e-02 -5.882353e-02 [122,] -2.297153e-18 -5.882353e-02 [123,] -5.882353e-02 -2.297153e-18 [124,] -5.882353e-02 -5.882353e-02 [125,] -2.297153e-18 -5.882353e-02 [126,] -5.882353e-02 -2.297153e-18 [127,] -5.882353e-02 -5.882353e-02 [128,] -5.882353e-02 -5.882353e-02 [129,] -5.882353e-02 -5.882353e-02 [130,] -5.882353e-02 -5.882353e-02 [131,] -5.882353e-02 -5.882353e-02 [132,] -5.882353e-02 -5.882353e-02 [133,] -5.882353e-02 -5.882353e-02 [134,] -5.882353e-02 -5.882353e-02 [135,] -5.882353e-02 -5.882353e-02 [136,] -5.882353e-02 -5.882353e-02 [137,] -2.297153e-18 -5.882353e-02 [138,] -2.297153e-18 -2.297153e-18 [139,] -5.882353e-02 -2.297153e-18 [140,] 9.411765e-01 -5.882353e-02 [141,] -2.297153e-18 9.411765e-01 [142,] -5.882353e-02 -2.297153e-18 [143,] -5.882353e-02 -5.882353e-02 [144,] -5.882353e-02 -5.882353e-02 [145,] -2.297153e-18 -5.882353e-02 [146,] -2.297153e-18 -2.297153e-18 [147,] -2.297153e-18 -2.297153e-18 [148,] -5.882353e-02 -2.297153e-18 [149,] -5.882353e-02 -5.882353e-02 [150,] -5.882353e-02 -5.882353e-02 [151,] 9.411765e-01 -5.882353e-02 [152,] 9.411765e-01 9.411765e-01 [153,] -5.882353e-02 9.411765e-01 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -4.761905e-02 -2.608696e-01 2 -4.761905e-02 -4.761905e-02 3 -4.761905e-02 -4.761905e-02 4 -4.761905e-02 -4.761905e-02 5 -4.761905e-02 -4.761905e-02 6 -4.761905e-02 -4.761905e-02 7 -2.608696e-01 -4.761905e-02 8 -4.761905e-02 -2.608696e-01 9 -4.761905e-02 -4.761905e-02 10 -2.608696e-01 -4.761905e-02 11 -4.761905e-02 -2.608696e-01 12 -4.761905e-02 -4.761905e-02 13 -2.608696e-01 -4.761905e-02 14 -4.761905e-02 -2.608696e-01 15 -2.608696e-01 -4.761905e-02 16 7.391304e-01 -2.608696e-01 17 -2.608696e-01 7.391304e-01 18 -4.761905e-02 -2.608696e-01 19 7.391304e-01 -4.761905e-02 20 -4.761905e-02 7.391304e-01 21 -4.761905e-02 -4.761905e-02 22 -4.761905e-02 -4.761905e-02 23 -4.761905e-02 -4.761905e-02 24 -2.608696e-01 -4.761905e-02 25 -4.761905e-02 -2.608696e-01 26 -4.761905e-02 -4.761905e-02 27 -4.761905e-02 -4.761905e-02 28 -4.761905e-02 -4.761905e-02 29 -4.761905e-02 -4.761905e-02 30 -4.761905e-02 -4.761905e-02 31 -4.761905e-02 -4.761905e-02 32 -4.761905e-02 -4.761905e-02 33 -2.608696e-01 -4.761905e-02 34 -4.761905e-02 -2.608696e-01 35 -4.761905e-02 -4.761905e-02 36 -2.608696e-01 -4.761905e-02 37 -4.761905e-02 -2.608696e-01 38 -4.761905e-02 -4.761905e-02 39 -2.608696e-01 -4.761905e-02 40 9.523810e-01 -2.608696e-01 41 -4.761905e-02 9.523810e-01 42 -4.761905e-02 -4.761905e-02 43 -2.608696e-01 -4.761905e-02 44 -4.761905e-02 -2.608696e-01 45 -4.761905e-02 -4.761905e-02 46 -4.761905e-02 -4.761905e-02 47 -4.761905e-02 -4.761905e-02 48 -4.761905e-02 -4.761905e-02 49 -4.761905e-02 -4.761905e-02 50 -2.608696e-01 -4.761905e-02 51 7.391304e-01 -2.608696e-01 52 -4.761905e-02 7.391304e-01 53 9.523810e-01 -4.761905e-02 54 -4.761905e-02 9.523810e-01 55 -2.608696e-01 -4.761905e-02 56 -4.761905e-02 -2.608696e-01 57 -4.761905e-02 -4.761905e-02 58 -4.761905e-02 -4.761905e-02 59 7.391304e-01 -4.761905e-02 60 -2.608696e-01 7.391304e-01 61 -4.761905e-02 -2.608696e-01 62 -4.761905e-02 -4.761905e-02 63 -2.608696e-01 -4.761905e-02 64 -4.761905e-02 -2.608696e-01 65 -4.761905e-02 -4.761905e-02 66 7.391304e-01 -4.761905e-02 67 -4.761905e-02 7.391304e-01 68 -4.761905e-02 -4.761905e-02 69 -4.761905e-02 -4.761905e-02 70 -4.761905e-02 -4.761905e-02 71 -4.761905e-02 -4.761905e-02 72 -4.761905e-02 -4.761905e-02 73 -4.761905e-02 -4.761905e-02 74 -4.761905e-02 -4.761905e-02 75 -2.608696e-01 -4.761905e-02 76 -4.761905e-02 -2.608696e-01 77 -4.761905e-02 -4.761905e-02 78 7.391304e-01 -4.761905e-02 79 -2.608696e-01 7.391304e-01 80 -4.761905e-02 -2.608696e-01 81 -4.761905e-02 -4.761905e-02 82 -4.761905e-02 -4.761905e-02 83 9.523810e-01 -4.761905e-02 84 -4.761905e-02 9.523810e-01 85 -4.761905e-02 -4.761905e-02 86 -5.882353e-02 -4.761905e-02 87 -2.297153e-18 -5.882353e-02 88 -5.882353e-02 -2.297153e-18 89 -5.882353e-02 -5.882353e-02 90 -5.882353e-02 -5.882353e-02 91 -2.297153e-18 -5.882353e-02 92 -5.882353e-02 -2.297153e-18 93 -5.882353e-02 -5.882353e-02 94 -2.297153e-18 -5.882353e-02 95 -5.882353e-02 -2.297153e-18 96 -2.297153e-18 -5.882353e-02 97 -5.882353e-02 -2.297153e-18 98 -5.882353e-02 -5.882353e-02 99 -5.882353e-02 -5.882353e-02 100 -5.882353e-02 -5.882353e-02 101 -5.882353e-02 -5.882353e-02 102 -5.882353e-02 -5.882353e-02 103 -5.882353e-02 -5.882353e-02 104 -2.297153e-18 -5.882353e-02 105 -5.882353e-02 -2.297153e-18 106 -5.882353e-02 -5.882353e-02 107 -2.297153e-18 -5.882353e-02 108 -5.882353e-02 -2.297153e-18 109 -5.882353e-02 -5.882353e-02 110 -2.297153e-18 -5.882353e-02 111 -2.297153e-18 -2.297153e-18 112 -5.882353e-02 -2.297153e-18 113 -2.297153e-18 -5.882353e-02 114 -5.882353e-02 -2.297153e-18 115 -5.882353e-02 -5.882353e-02 116 -5.882353e-02 -5.882353e-02 117 -5.882353e-02 -5.882353e-02 118 -5.882353e-02 -5.882353e-02 119 -5.882353e-02 -5.882353e-02 120 -5.882353e-02 -5.882353e-02 121 -5.882353e-02 -5.882353e-02 122 -2.297153e-18 -5.882353e-02 123 -5.882353e-02 -2.297153e-18 124 -5.882353e-02 -5.882353e-02 125 -2.297153e-18 -5.882353e-02 126 -5.882353e-02 -2.297153e-18 127 -5.882353e-02 -5.882353e-02 128 -5.882353e-02 -5.882353e-02 129 -5.882353e-02 -5.882353e-02 130 -5.882353e-02 -5.882353e-02 131 -5.882353e-02 -5.882353e-02 132 -5.882353e-02 -5.882353e-02 133 -5.882353e-02 -5.882353e-02 134 -5.882353e-02 -5.882353e-02 135 -5.882353e-02 -5.882353e-02 136 -5.882353e-02 -5.882353e-02 137 -2.297153e-18 -5.882353e-02 138 -2.297153e-18 -2.297153e-18 139 -5.882353e-02 -2.297153e-18 140 9.411765e-01 -5.882353e-02 141 -2.297153e-18 9.411765e-01 142 -5.882353e-02 -2.297153e-18 143 -5.882353e-02 -5.882353e-02 144 -5.882353e-02 -5.882353e-02 145 -2.297153e-18 -5.882353e-02 146 -2.297153e-18 -2.297153e-18 147 -2.297153e-18 -2.297153e-18 148 -5.882353e-02 -2.297153e-18 149 -5.882353e-02 -5.882353e-02 150 -5.882353e-02 -5.882353e-02 151 9.411765e-01 -5.882353e-02 152 9.411765e-01 9.411765e-01 153 -5.882353e-02 9.411765e-01 > 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/7hbu81356032183.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/8uzya1356032183.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/9xj611356032184.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/1076yf1356032184.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='') + } + } Error: subscript out of bounds Execution halted