R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(-45.6
+ ,16.1
+ ,23.9
+ ,39.3
+ ,-39.4
+ ,-0.3
+ ,17.3
+ ,17.7
+ ,31.4
+ ,-28.6
+ ,-17.2
+ ,-79
+ ,-47.9
+ ,9.1
+ ,10.6
+ ,-23.9
+ ,-45
+ ,-42.2
+ ,43.2
+ ,32.1
+ ,-15.3
+ ,21.8
+ ,-12
+ ,-95.8
+ ,-14.3
+ ,47.8
+ ,64.8
+ ,40.2
+ ,-28.8
+ ,23.5
+ ,70.3
+ ,12.3
+ ,43.5
+ ,-30.1
+ ,-5.3
+ ,-24
+ ,11.1
+ ,21.5
+ ,38.5
+ ,16.8
+ ,-36.2
+ ,6
+ ,26.6
+ ,-8
+ ,13.2
+ ,-23.6
+ ,19.4
+ ,-46.2
+ ,-8.2
+ ,33.8
+ ,16.6
+ ,5.4
+ ,-25
+ ,-5.3
+ ,16.7
+ ,19
+ ,24.8
+ ,-11.4
+ ,4.9
+ ,-58.7
+ ,16.8
+ ,13.6
+ ,6.4
+ ,22.8
+ ,-19.6
+ ,2.2
+ ,19.8
+ ,-10.7
+ ,4.7
+ ,-44.5
+ ,-34.7
+ ,-119.7
+ ,-42.2
+ ,-5.4
+ ,19.1
+ ,18.8
+ ,-2.3
+ ,0.2
+ ,20.9
+ ,3.7
+ ,50.4
+ ,-18.6
+ ,10.6
+ ,-66
+ ,10
+ ,27.2
+ ,13.5
+ ,47.2
+ ,-20.3
+ ,23.1
+ ,12.6
+ ,19.8
+ ,5.4
+ ,-25.2
+ ,-6.5
+ ,-46.5
+ ,-2.6
+ ,-0.3
+ ,38.5
+ ,-8.9
+ ,-38
+ ,19.5
+ ,51.7
+ ,19.4
+ ,18.2
+ ,-50.8
+ ,-6.1
+ ,-54.6
+ ,12.1
+ ,26.3
+ ,19.5
+ ,-0.8
+ ,-49.6
+ ,28.8
+ ,31.7
+ ,2.3
+ ,3.8
+ ,-66.2
+ ,-20.5
+ ,-113.2
+ ,-65.2
+ ,-3.9
+ ,9.1
+ ,23.2
+ ,-39.1
+ ,12.5
+ ,49.1
+ ,54.9
+ ,30.8
+ ,-3.5
+ ,-28.3
+ ,-61
+ ,-2
+ ,40
+ ,74
+ ,23.1
+ ,-45.3
+ ,17.5
+ ,25.8
+ ,15.2
+ ,-3.6
+ ,-40.5
+ ,11.5
+ ,-59.8
+ ,23.3
+ ,-27.8
+ ,55.7
+ ,22.7
+ ,-79.2
+ ,28.8
+ ,17.3
+ ,39.6
+ ,-22.2
+ ,-43
+ ,-50.3
+ ,-86.5
+ ,-31.9
+ ,23.1
+ ,53.6
+ ,21.6
+ ,-64.2
+ ,35.2
+ ,52.1
+ ,40.6
+ ,17.1
+ ,-7.8
+ ,-10
+ ,-58
+ ,14
+ ,15.8
+ ,46
+ ,-8.9
+ ,-26.7
+ ,39
+ ,-1.3
+ ,38.7
+ ,22.1
+ ,-49.2
+ ,-3.4
+ ,-86.7
+ ,-24.3
+ ,42.8
+ ,44.9
+ ,4.4
+ ,-60.5
+ ,41.4
+ ,38.5
+ ,28.5
+ ,7.6
+ ,-46.4
+ ,7
+ ,-73
+ ,5.7
+ ,23.6
+ ,39.4
+ ,30.3
+ ,-92.5
+ ,77.8
+ ,12.4
+ ,28.9
+ ,6.4
+ ,-12
+ ,-9.1
+ ,-53.2
+ ,-23.1
+ ,47.3
+ ,20.7
+ ,27.8
+ ,-84.3
+ ,62.8
+ ,26.4
+ ,32.3
+ ,13.3
+ ,-17.9
+ ,10
+ ,-45.6
+ ,13.5
+ ,11.9
+ ,26
+ ,-6.3
+ ,-79.9
+ ,54.2
+ ,22.9
+ ,31.8
+ ,3.8
+ ,-11.4
+ ,-8.6
+ ,-49.4
+ ,-2.5
+ ,23
+ ,29
+ ,20.6
+ ,-117
+ ,37.9
+ ,30.7
+ ,4.7
+ ,-5.7
+ ,4.9
+ ,18.3
+ ,-35.4
+ ,-21.3
+ ,35.8
+ ,43.8
+ ,18.7
+ ,-131.1
+ ,39.8
+ ,44.5
+ ,16.5
+ ,9.7
+ ,-6.6
+ ,15.8
+ ,-45.7
+ ,-4.8
+ ,17.6
+ ,20.5
+ ,24.2
+ ,-109
+ ,20.8
+ ,31.2
+ ,-8.8
+ ,11.8
+ ,13
+ ,8.3
+ ,-77.9
+ ,-38.8
+ ,6.1
+ ,18.1
+ ,16.8
+ ,-128.5
+ ,15.9
+ ,29
+ ,-7.2
+ ,3.3
+ ,-34.8
+ ,-2.9
+ ,-77.8
+ ,-2.8
+ ,26.7
+ ,48.1
+ ,30
+ ,-109.6
+ ,16
+ ,26.9
+ ,22.1
+ ,27
+ ,-24.5
+ ,12
+ ,-75.2
+ ,3.5
+ ,19.7
+ ,51.8
+ ,35.3
+ ,-108.2
+ ,25.3
+ ,31.6
+ ,19.9
+ ,18.8
+ ,20.4
+ ,15
+ ,-55.9
+ ,-17
+ ,33.3
+ ,33.8
+ ,37.5
+ ,-104.8
+ ,29.7
+ ,34.2
+ ,4.3
+ ,40.2
+ ,-29.3
+ ,-0.2
+ ,-95
+ ,-13.2
+ ,38.5
+ ,45.4
+ ,15.7
+ ,-123.6
+ ,12
+ ,37.5
+ ,-31.7
+ ,15.8
+ ,-64.1
+ ,-42.1
+ ,-207.4
+ ,-12.9
+ ,-5
+ ,53.9
+ ,19.7
+ ,-94.6
+ ,36
+ ,51.3
+ ,17.4
+ ,27.8
+ ,1.3
+ ,3.6
+ ,-97.9
+ ,14.1
+ ,50.8
+ ,63.5
+ ,58.6
+ ,-135.1
+ ,7.8
+ ,25.5
+ ,29.6
+ ,19.3
+ ,-26.2
+ ,7.3
+ ,-82.6
+ ,-26.1
+ ,55.3
+ ,98.8
+ ,41.7
+ ,-130.2
+ ,51.2
+ ,18.4
+ ,32
+ ,21.6
+ ,-12.5
+ ,46.6
+ ,-101.7
+ ,15.8
+ ,26
+ ,79.1
+ ,23.1
+ ,-86.9
+ ,-11.2
+ ,50.7
+ ,13.4
+ ,33.7
+ ,-16.9
+ ,-9.6)
+ ,dim=c(1
+ ,371)
+ ,dimnames=list(c('y')
+ ,1:371))
> y <- array(NA,dim=c(1,371),dimnames=list(c('y'),1:371))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '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
y t
1 -45.6 1
2 16.1 2
3 23.9 3
4 39.3 4
5 -39.4 5
6 -0.3 6
7 17.3 7
8 17.7 8
9 31.4 9
10 -28.6 10
11 -17.2 11
12 -79.0 12
13 -47.9 13
14 9.1 14
15 10.6 15
16 -23.9 16
17 -45.0 17
18 -42.2 18
19 43.2 19
20 32.1 20
21 -15.3 21
22 21.8 22
23 -12.0 23
24 -95.8 24
25 -14.3 25
26 47.8 26
27 64.8 27
28 40.2 28
29 -28.8 29
30 23.5 30
31 70.3 31
32 12.3 32
33 43.5 33
34 -30.1 34
35 -5.3 35
36 -24.0 36
37 11.1 37
38 21.5 38
39 38.5 39
40 16.8 40
41 -36.2 41
42 6.0 42
43 26.6 43
44 -8.0 44
45 13.2 45
46 -23.6 46
47 19.4 47
48 -46.2 48
49 -8.2 49
50 33.8 50
51 16.6 51
52 5.4 52
53 -25.0 53
54 -5.3 54
55 16.7 55
56 19.0 56
57 24.8 57
58 -11.4 58
59 4.9 59
60 -58.7 60
61 16.8 61
62 13.6 62
63 6.4 63
64 22.8 64
65 -19.6 65
66 2.2 66
67 19.8 67
68 -10.7 68
69 4.7 69
70 -44.5 70
71 -34.7 71
72 -119.7 72
73 -42.2 73
74 -5.4 74
75 19.1 75
76 18.8 76
77 -2.3 77
78 0.2 78
79 20.9 79
80 3.7 80
81 50.4 81
82 -18.6 82
83 10.6 83
84 -66.0 84
85 10.0 85
86 27.2 86
87 13.5 87
88 47.2 88
89 -20.3 89
90 23.1 90
91 12.6 91
92 19.8 92
93 5.4 93
94 -25.2 94
95 -6.5 95
96 -46.5 96
97 -2.6 97
98 -0.3 98
99 38.5 99
100 -8.9 100
101 -38.0 101
102 19.5 102
103 51.7 103
104 19.4 104
105 18.2 105
106 -50.8 106
107 -6.1 107
108 -54.6 108
109 12.1 109
110 26.3 110
111 19.5 111
112 -0.8 112
113 -49.6 113
114 28.8 114
115 31.7 115
116 2.3 116
117 3.8 117
118 -66.2 118
119 -20.5 119
120 -113.2 120
121 -65.2 121
122 -3.9 122
123 9.1 123
124 23.2 124
125 -39.1 125
126 12.5 126
127 49.1 127
128 54.9 128
129 30.8 129
130 -3.5 130
131 -28.3 131
132 -61.0 132
133 -2.0 133
134 40.0 134
135 74.0 135
136 23.1 136
137 -45.3 137
138 17.5 138
139 25.8 139
140 15.2 140
141 -3.6 141
142 -40.5 142
143 11.5 143
144 -59.8 144
145 23.3 145
146 -27.8 146
147 55.7 147
148 22.7 148
149 -79.2 149
150 28.8 150
151 17.3 151
152 39.6 152
153 -22.2 153
154 -43.0 154
155 -50.3 155
156 -86.5 156
157 -31.9 157
158 23.1 158
159 53.6 159
160 21.6 160
161 -64.2 161
162 35.2 162
163 52.1 163
164 40.6 164
165 17.1 165
166 -7.8 166
167 -10.0 167
168 -58.0 168
169 14.0 169
170 15.8 170
171 46.0 171
172 -8.9 172
173 -26.7 173
174 39.0 174
175 -1.3 175
176 38.7 176
177 22.1 177
178 -49.2 178
179 -3.4 179
180 -86.7 180
181 -24.3 181
182 42.8 182
183 44.9 183
184 4.4 184
185 -60.5 185
186 41.4 186
187 38.5 187
188 28.5 188
189 7.6 189
190 -46.4 190
191 7.0 191
192 -73.0 192
193 5.7 193
194 23.6 194
195 39.4 195
196 30.3 196
197 -92.5 197
198 77.8 198
199 12.4 199
200 28.9 200
201 6.4 201
202 -12.0 202
203 -9.1 203
204 -53.2 204
205 -23.1 205
206 47.3 206
207 20.7 207
208 27.8 208
209 -84.3 209
210 62.8 210
211 26.4 211
212 32.3 212
213 13.3 213
214 -17.9 214
215 10.0 215
216 -45.6 216
217 13.5 217
218 11.9 218
219 26.0 219
220 -6.3 220
221 -79.9 221
222 54.2 222
223 22.9 223
224 31.8 224
225 3.8 225
226 -11.4 226
227 -8.6 227
228 -49.4 228
229 -2.5 229
230 23.0 230
231 29.0 231
232 20.6 232
233 -117.0 233
234 37.9 234
235 30.7 235
236 4.7 236
237 -5.7 237
238 4.9 238
239 18.3 239
240 -35.4 240
241 -21.3 241
242 35.8 242
243 43.8 243
244 18.7 244
245 -131.1 245
246 39.8 246
247 44.5 247
248 16.5 248
249 9.7 249
250 -6.6 250
251 15.8 251
252 -45.7 252
253 -4.8 253
254 17.6 254
255 20.5 255
256 24.2 256
257 -109.0 257
258 20.8 258
259 31.2 259
260 -8.8 260
261 11.8 261
262 13.0 262
263 8.3 263
264 -77.9 264
265 -38.8 265
266 6.1 266
267 18.1 267
268 16.8 268
269 -128.5 269
270 15.9 270
271 29.0 271
272 -7.2 272
273 3.3 273
274 -34.8 274
275 -2.9 275
276 -77.8 276
277 -2.8 277
278 26.7 278
279 48.1 279
280 30.0 280
281 -109.6 281
282 16.0 282
283 26.9 283
284 22.1 284
285 27.0 285
286 -24.5 286
287 12.0 287
288 -75.2 288
289 3.5 289
290 19.7 290
291 51.8 291
292 35.3 292
293 -108.2 293
294 25.3 294
295 31.6 295
296 19.9 296
297 18.8 297
298 20.4 298
299 15.0 299
300 -55.9 300
301 -17.0 301
302 33.3 302
303 33.8 303
304 37.5 304
305 -104.8 305
306 29.7 306
307 34.2 307
308 4.3 308
309 40.2 309
310 -29.3 310
311 -0.2 311
312 -95.0 312
313 -13.2 313
314 38.5 314
315 45.4 315
316 15.7 316
317 -123.6 317
318 12.0 318
319 37.5 319
320 -31.7 320
321 15.8 321
322 -64.1 322
323 -42.1 323
324 -207.4 324
325 -12.9 325
326 -5.0 326
327 53.9 327
328 19.7 328
329 -94.6 329
330 36.0 330
331 51.3 331
332 17.4 332
333 27.8 333
334 1.3 334
335 3.6 335
336 -97.9 336
337 14.1 337
338 50.8 338
339 63.5 339
340 58.6 340
341 -135.1 341
342 7.8 342
343 25.5 343
344 29.6 344
345 19.3 345
346 -26.2 346
347 7.3 347
348 -82.6 348
349 -26.1 349
350 55.3 350
351 98.8 351
352 41.7 352
353 -130.2 353
354 51.2 354
355 18.4 355
356 32.0 356
357 21.6 357
358 -12.5 358
359 46.6 359
360 -101.7 360
361 15.8 361
362 26.0 362
363 79.1 363
364 23.1 364
365 -86.9 365
366 -11.2 366
367 50.7 367
368 13.4 368
369 33.7 369
370 -16.9 370
371 -9.6 371
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) t
-0.7345174 -0.0009404
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-206.36 -20.00 12.78 27.23 99.86
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.7345174 4.4875842 -0.164 0.870
t -0.0009404 0.0209085 -0.045 0.964
Residual standard error: 43.13 on 369 degrees of freedom
Multiple R-squared: 5.482e-06, Adjusted R-squared: -0.002705
F-statistic: 0.002023 on 1 and 369 DF, p-value: 0.9641
> 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,] 6.283900e-01 7.432200e-01 0.37161001
[2,] 4.674244e-01 9.348488e-01 0.53257561
[3,] 3.281897e-01 6.563794e-01 0.67181030
[4,] 2.144046e-01 4.288092e-01 0.78559542
[5,] 1.369879e-01 2.739758e-01 0.86301208
[6,] 1.758634e-01 3.517268e-01 0.82413660
[7,] 1.331785e-01 2.663571e-01 0.86682147
[8,] 2.793008e-01 5.586017e-01 0.72069915
[9,] 2.220783e-01 4.441565e-01 0.77792175
[10,] 2.168177e-01 4.336354e-01 0.78318231
[11,] 1.950636e-01 3.901273e-01 0.80493636
[12,] 1.423927e-01 2.847854e-01 0.85760732
[13,] 1.123222e-01 2.246444e-01 0.88767780
[14,] 8.190608e-02 1.638122e-01 0.91809392
[15,] 1.664812e-01 3.329624e-01 0.83351880
[16,] 1.850100e-01 3.700199e-01 0.81499004
[17,] 1.411490e-01 2.822979e-01 0.85885103
[18,] 1.261850e-01 2.523701e-01 0.87381496
[19,] 9.419107e-02 1.883821e-01 0.90580893
[20,] 2.186883e-01 4.373765e-01 0.78131174
[21,] 1.754396e-01 3.508793e-01 0.82456036
[22,] 2.417773e-01 4.835546e-01 0.75822272
[23,] 3.424569e-01 6.849138e-01 0.65754312
[24,] 3.338486e-01 6.676971e-01 0.66615145
[25,] 3.102315e-01 6.204631e-01 0.68976846
[26,] 2.733809e-01 5.467617e-01 0.72661915
[27,] 3.344498e-01 6.688995e-01 0.66555024
[28,] 2.854168e-01 5.708337e-01 0.71458317
[29,] 2.620576e-01 5.241153e-01 0.73794237
[30,] 2.687102e-01 5.374204e-01 0.73128980
[31,] 2.332699e-01 4.665397e-01 0.76673015
[32,] 2.195399e-01 4.390798e-01 0.78046011
[33,] 1.829143e-01 3.658287e-01 0.81708565
[34,] 1.531744e-01 3.063488e-01 0.84682558
[35,] 1.366037e-01 2.732074e-01 0.86339628
[36,] 1.109711e-01 2.219422e-01 0.88902891
[37,] 1.214849e-01 2.429698e-01 0.87851512
[38,] 9.843444e-02 1.968689e-01 0.90156556
[39,] 8.151482e-02 1.630296e-01 0.91848518
[40,] 6.827681e-02 1.365536e-01 0.93172319
[41,] 5.377262e-02 1.075452e-01 0.94622738
[42,] 5.010178e-02 1.002036e-01 0.94989822
[43,] 3.961012e-02 7.922024e-02 0.96038988
[44,] 4.899042e-02 9.798083e-02 0.95100958
[45,] 3.935907e-02 7.871814e-02 0.96064093
[46,] 3.454275e-02 6.908550e-02 0.96545725
[47,] 2.698485e-02 5.396971e-02 0.97301515
[48,] 2.063005e-02 4.126010e-02 0.97936995
[49,] 1.882016e-02 3.764032e-02 0.98117984
[50,] 1.452520e-02 2.905040e-02 0.98547480
[51,] 1.111417e-02 2.222835e-02 0.98888583
[52,] 8.496656e-03 1.699331e-02 0.99150334
[53,] 6.639844e-03 1.327969e-02 0.99336016
[54,] 5.257788e-03 1.051558e-02 0.99474221
[55,] 3.823597e-03 7.647195e-03 0.99617640
[56,] 6.647200e-03 1.329440e-02 0.99335280
[57,] 5.056956e-03 1.011391e-02 0.99494304
[58,] 3.754375e-03 7.508749e-03 0.99624563
[59,] 2.719895e-03 5.439789e-03 0.99728011
[60,] 2.087206e-03 4.174411e-03 0.99791279
[61,] 1.726209e-03 3.452418e-03 0.99827379
[62,] 1.228697e-03 2.457394e-03 0.99877130
[63,] 9.114174e-04 1.822835e-03 0.99908858
[64,] 6.804827e-04 1.360965e-03 0.99931952
[65,] 4.731352e-04 9.462704e-04 0.99952686
[66,] 5.935361e-04 1.187072e-03 0.99940646
[67,] 5.756615e-04 1.151323e-03 0.99942434
[68,] 7.520202e-03 1.504040e-02 0.99247980
[69,] 7.220754e-03 1.444151e-02 0.99277925
[70,] 5.477486e-03 1.095497e-02 0.99452251
[71,] 4.629475e-03 9.258949e-03 0.99537053
[72,] 3.855656e-03 7.711313e-03 0.99614434
[73,] 2.878630e-03 5.757261e-03 0.99712137
[74,] 2.139762e-03 4.279524e-03 0.99786024
[75,] 1.777480e-03 3.554961e-03 0.99822252
[76,] 1.313196e-03 2.626392e-03 0.99868680
[77,] 1.674062e-03 3.348123e-03 0.99832594
[78,] 1.305491e-03 2.610982e-03 0.99869451
[79,] 9.798918e-04 1.959784e-03 0.99902011
[80,] 1.628791e-03 3.257582e-03 0.99837121
[81,] 1.242272e-03 2.484544e-03 0.99875773
[82,] 1.087350e-03 2.174700e-03 0.99891265
[83,] 8.320902e-04 1.664180e-03 0.99916791
[84,] 9.654048e-04 1.930810e-03 0.99903460
[85,] 7.669924e-04 1.533985e-03 0.99923301
[86,] 6.223785e-04 1.244757e-03 0.99937762
[87,] 4.640438e-04 9.280876e-04 0.99953596
[88,] 3.596618e-04 7.193237e-04 0.99964034
[89,] 2.582217e-04 5.164433e-04 0.99974178
[90,] 2.150437e-04 4.300873e-04 0.99978496
[91,] 1.542108e-04 3.084216e-04 0.99984579
[92,] 1.777530e-04 3.555061e-04 0.99982225
[93,] 1.258122e-04 2.516244e-04 0.99987419
[94,] 8.841690e-05 1.768338e-04 0.99991158
[95,] 8.840341e-05 1.768068e-04 0.99991160
[96,] 6.301677e-05 1.260335e-04 0.99993698
[97,] 6.132443e-05 1.226489e-04 0.99993868
[98,] 4.717355e-05 9.434709e-05 0.99995283
[99,] 6.187598e-05 1.237520e-04 0.99993812
[100,] 4.672904e-05 9.345808e-05 0.99995327
[101,] 3.467502e-05 6.935005e-05 0.99996532
[102,] 4.611436e-05 9.222872e-05 0.99995389
[103,] 3.242707e-05 6.485413e-05 0.99996757
[104,] 4.520869e-05 9.041738e-05 0.99995479
[105,] 3.286365e-05 6.572730e-05 0.99996714
[106,] 2.711666e-05 5.423332e-05 0.99997288
[107,] 2.056613e-05 4.113226e-05 0.99997943
[108,] 1.417645e-05 2.835290e-05 0.99998582
[109,] 1.769648e-05 3.539296e-05 0.99998230
[110,] 1.506726e-05 3.013452e-05 0.99998493
[111,] 1.323074e-05 2.646147e-05 0.99998677
[112,] 9.090883e-06 1.818177e-05 0.99999091
[113,] 6.220147e-06 1.244029e-05 0.99999378
[114,] 1.255193e-05 2.510385e-05 0.99998745
[115,] 9.461478e-06 1.892296e-05 0.99999054
[116,] 1.026523e-04 2.053046e-04 0.99989735
[117,] 1.557618e-04 3.115236e-04 0.99984424
[118,] 1.134931e-04 2.269862e-04 0.99988651
[119,] 8.568869e-05 1.713774e-04 0.99991431
[120,] 7.254790e-05 1.450958e-04 0.99992745
[121,] 6.672211e-05 1.334442e-04 0.99993328
[122,] 5.107743e-05 1.021549e-04 0.99994892
[123,] 6.609692e-05 1.321938e-04 0.99993390
[124,] 9.381472e-05 1.876294e-04 0.99990619
[125,] 8.418048e-05 1.683610e-04 0.99991582
[126,] 6.082766e-05 1.216553e-04 0.99993917
[127,] 5.042845e-05 1.008569e-04 0.99994957
[128,] 7.227760e-05 1.445552e-04 0.99992772
[129,] 5.223266e-05 1.044653e-04 0.99994777
[130,] 5.457425e-05 1.091485e-04 0.99994543
[131,] 1.213647e-04 2.427295e-04 0.99987864
[132,] 9.890643e-05 1.978129e-04 0.99990109
[133,] 1.051541e-04 2.103081e-04 0.99989485
[134,] 8.204650e-05 1.640930e-04 0.99991795
[135,] 6.840350e-05 1.368070e-04 0.99993160
[136,] 5.199163e-05 1.039833e-04 0.99994801
[137,] 3.757352e-05 7.514704e-05 0.99996243
[138,] 3.716269e-05 7.432539e-05 0.99996284
[139,] 2.751136e-05 5.502272e-05 0.99997249
[140,] 3.934163e-05 7.868327e-05 0.99996066
[141,] 3.203261e-05 6.406522e-05 0.99996797
[142,] 2.641621e-05 5.283243e-05 0.99997358
[143,] 3.630595e-05 7.261191e-05 0.99996369
[144,] 2.912769e-05 5.825538e-05 0.99997087
[145,] 6.862629e-05 1.372526e-04 0.99993137
[146,] 5.947187e-05 1.189437e-04 0.99994053
[147,] 4.613283e-05 9.226565e-05 0.99995387
[148,] 4.562880e-05 9.125759e-05 0.99995437
[149,] 3.620693e-05 7.241385e-05 0.99996379
[150,] 3.684119e-05 7.368237e-05 0.99996316
[151,] 4.207051e-05 8.414102e-05 0.99995793
[152,] 1.128878e-04 2.257756e-04 0.99988711
[153,] 9.698229e-05 1.939646e-04 0.99990302
[154,] 8.100449e-05 1.620090e-04 0.99991900
[155,] 1.047626e-04 2.095252e-04 0.99989524
[156,] 8.544175e-05 1.708835e-04 0.99991456
[157,] 1.245038e-04 2.490076e-04 0.99987550
[158,] 1.179759e-04 2.359518e-04 0.99988202
[159,] 1.441516e-04 2.883032e-04 0.99985585
[160,] 1.444110e-04 2.888220e-04 0.99985559
[161,] 1.135172e-04 2.270344e-04 0.99988648
[162,] 8.484427e-05 1.696885e-04 0.99991516
[163,] 6.349410e-05 1.269882e-04 0.99993651
[164,] 8.280976e-05 1.656195e-04 0.99991719
[165,] 6.366864e-05 1.273373e-04 0.99993633
[166,] 4.921689e-05 9.843379e-05 0.99995078
[167,] 5.341002e-05 1.068200e-04 0.99994659
[168,] 3.959292e-05 7.918584e-05 0.99996041
[169,] 3.263758e-05 6.527515e-05 0.99996736
[170,] 3.169648e-05 6.339295e-05 0.99996830
[171,] 2.301614e-05 4.603228e-05 0.99997698
[172,] 2.212718e-05 4.425435e-05 0.99997787
[173,] 1.753374e-05 3.506747e-05 0.99998247
[174,] 1.983944e-05 3.967888e-05 0.99998016
[175,] 1.432883e-05 2.865766e-05 0.99998567
[176,] 4.026566e-05 8.053133e-05 0.99995973
[177,] 3.239894e-05 6.479789e-05 0.99996760
[178,] 3.338434e-05 6.676868e-05 0.99996662
[179,] 3.533598e-05 7.067196e-05 0.99996466
[180,] 2.587317e-05 5.174633e-05 0.99997413
[181,] 3.574937e-05 7.149874e-05 0.99996425
[182,] 3.578409e-05 7.156818e-05 0.99996422
[183,] 3.421337e-05 6.842674e-05 0.99996579
[184,] 2.892910e-05 5.785821e-05 0.99997107
[185,] 2.127074e-05 4.254147e-05 0.99997873
[186,] 2.264149e-05 4.528297e-05 0.99997736
[187,] 1.656559e-05 3.313117e-05 0.99998343
[188,] 3.079678e-05 6.159356e-05 0.99996920
[189,] 2.259979e-05 4.519958e-05 0.99997740
[190,] 1.821692e-05 3.643383e-05 0.99998178
[191,] 1.757430e-05 3.514859e-05 0.99998243
[192,] 1.505749e-05 3.011498e-05 0.99998494
[193,] 4.854859e-05 9.709717e-05 0.99995145
[194,] 9.915473e-05 1.983095e-04 0.99990085
[195,] 7.573080e-05 1.514616e-04 0.99992427
[196,] 6.460445e-05 1.292089e-04 0.99993540
[197,] 4.805722e-05 9.611444e-05 0.99995194
[198,] 3.613280e-05 7.226561e-05 0.99996387
[199,] 2.676155e-05 5.352310e-05 0.99997324
[200,] 3.159415e-05 6.318830e-05 0.99996841
[201,] 2.516875e-05 5.033751e-05 0.99997483
[202,] 2.725084e-05 5.450168e-05 0.99997275
[203,] 2.143746e-05 4.287491e-05 0.99997856
[204,] 1.788957e-05 3.577913e-05 0.99998211
[205,] 4.359581e-05 8.719163e-05 0.99995640
[206,] 6.239476e-05 1.247895e-04 0.99993761
[207,] 5.191534e-05 1.038307e-04 0.99994808
[208,] 4.575174e-05 9.150348e-05 0.99995425
[209,] 3.471357e-05 6.942714e-05 0.99996529
[210,] 2.669397e-05 5.338795e-05 0.99997331
[211,] 1.981040e-05 3.962080e-05 0.99998019
[212,] 2.038186e-05 4.076371e-05 0.99997962
[213,] 1.528781e-05 3.057563e-05 0.99998471
[214,] 1.133550e-05 2.267100e-05 0.99998866
[215,] 9.261551e-06 1.852310e-05 0.99999074
[216,] 6.656595e-06 1.331319e-05 0.99999334
[217,] 1.456310e-05 2.912620e-05 0.99998544
[218,] 1.787312e-05 3.574624e-05 0.99998213
[219,] 1.427607e-05 2.855215e-05 0.99998572
[220,] 1.246260e-05 2.492520e-05 0.99998754
[221,] 8.995305e-06 1.799061e-05 0.99999100
[222,] 6.545893e-06 1.309179e-05 0.99999345
[223,] 4.696511e-06 9.393021e-06 0.99999530
[224,] 5.084632e-06 1.016926e-05 0.99999492
[225,] 3.593676e-06 7.187353e-06 0.99999641
[226,] 2.826888e-06 5.653776e-06 0.99999717
[227,] 2.366357e-06 4.732715e-06 0.99999763
[228,] 1.821016e-06 3.642032e-06 0.99999818
[229,] 1.456265e-05 2.912530e-05 0.99998544
[230,] 1.371441e-05 2.742882e-05 0.99998629
[231,] 1.183545e-05 2.367089e-05 0.99998816
[232,] 8.531996e-06 1.706399e-05 0.99999147
[233,] 6.083628e-06 1.216726e-05 0.99999392
[234,] 4.341977e-06 8.683955e-06 0.99999566
[235,] 3.304075e-06 6.608150e-06 0.99999670
[236,] 2.842740e-06 5.685480e-06 0.99999716
[237,] 2.126585e-06 4.253169e-06 0.99999787
[238,] 1.946191e-06 3.892382e-06 0.99999805
[239,] 2.036970e-06 4.073940e-06 0.99999796
[240,] 1.550140e-06 3.100279e-06 0.99999845
[241,] 2.118610e-05 4.237219e-05 0.99997881
[242,] 2.063613e-05 4.127226e-05 0.99997936
[243,] 2.175006e-05 4.350012e-05 0.99997825
[244,] 1.674795e-05 3.349591e-05 0.99998325
[245,] 1.238035e-05 2.476070e-05 0.99998762
[246,] 8.845831e-06 1.769166e-05 0.99999115
[247,] 6.722441e-06 1.344488e-05 0.99999328
[248,] 6.505216e-06 1.301043e-05 0.99999349
[249,] 4.582657e-06 9.165313e-06 0.99999542
[250,] 3.491509e-06 6.983018e-06 0.99999651
[251,] 2.726952e-06 5.453904e-06 0.99999727
[252,] 2.214432e-06 4.428864e-06 0.99999779
[253,] 1.123877e-05 2.247753e-05 0.99998876
[254,] 8.854978e-06 1.770996e-05 0.99999115
[255,] 7.815536e-06 1.563107e-05 0.99999218
[256,] 5.519191e-06 1.103838e-05 0.99999448
[257,] 4.060084e-06 8.120168e-06 0.99999594
[258,] 3.003608e-06 6.007217e-06 0.99999700
[259,] 2.156831e-06 4.313662e-06 0.99999784
[260,] 3.978494e-06 7.956988e-06 0.99999602
[261,] 3.424172e-06 6.848343e-06 0.99999658
[262,] 2.410501e-06 4.821002e-06 0.99999759
[263,] 1.821583e-06 3.643167e-06 0.99999818
[264,] 1.361773e-06 2.723545e-06 0.99999864
[265,] 1.527511e-05 3.055021e-05 0.99998472
[266,] 1.146555e-05 2.293111e-05 0.99998853
[267,] 9.704211e-06 1.940842e-05 0.99999030
[268,] 6.776992e-06 1.355398e-05 0.99999322
[269,] 4.742670e-06 9.485340e-06 0.99999526
[270,] 3.874248e-06 7.748496e-06 0.99999613
[271,] 2.652807e-06 5.305614e-06 0.99999735
[272,] 5.016374e-06 1.003275e-05 0.99999498
[273,] 3.438651e-06 6.877302e-06 0.99999656
[274,] 2.747369e-06 5.494739e-06 0.99999725
[275,] 3.053963e-06 6.107926e-06 0.99999695
[276,] 2.583237e-06 5.166475e-06 0.99999742
[277,] 1.325346e-05 2.650692e-05 0.99998675
[278,] 9.766046e-06 1.953209e-05 0.99999023
[279,] 7.905556e-06 1.581111e-05 0.99999209
[280,] 6.119390e-06 1.223878e-05 0.99999388
[281,] 5.002786e-06 1.000557e-05 0.99999500
[282,] 3.651684e-06 7.303368e-06 0.99999635
[283,] 2.594757e-06 5.189514e-06 0.99999741
[284,] 4.454612e-06 8.909224e-06 0.99999555
[285,] 3.028739e-06 6.057477e-06 0.99999697
[286,] 2.243510e-06 4.487019e-06 0.99999776
[287,] 2.718325e-06 5.436650e-06 0.99999728
[288,] 2.498855e-06 4.997711e-06 0.99999750
[289,] 1.167089e-05 2.334178e-05 0.99998833
[290,] 9.266010e-06 1.853202e-05 0.99999073
[291,] 7.992399e-06 1.598480e-05 0.99999201
[292,] 6.093800e-06 1.218760e-05 0.99999391
[293,] 4.614928e-06 9.229857e-06 0.99999539
[294,] 3.567999e-06 7.135999e-06 0.99999643
[295,] 2.631460e-06 5.262921e-06 0.99999737
[296,] 2.690305e-06 5.380610e-06 0.99999731
[297,] 1.803896e-06 3.607792e-06 0.99999820
[298,] 1.626531e-06 3.253061e-06 0.99999837
[299,] 1.514916e-06 3.029832e-06 0.99999849
[300,] 1.553850e-06 3.107699e-06 0.99999845
[301,] 5.672203e-06 1.134441e-05 0.99999433
[302,] 4.952883e-06 9.905767e-06 0.99999505
[303,] 4.759354e-06 9.518709e-06 0.99999524
[304,] 3.293873e-06 6.587746e-06 0.99999671
[305,] 3.721432e-06 7.442864e-06 0.99999628
[306,] 2.563087e-06 5.126175e-06 0.99999744
[307,] 1.728744e-06 3.457488e-06 0.99999827
[308,] 4.228031e-06 8.456062e-06 0.99999577
[309,] 2.741013e-06 5.482027e-06 0.99999726
[310,] 2.853357e-06 5.706714e-06 0.99999715
[311,] 3.637704e-06 7.275409e-06 0.99999636
[312,] 2.866183e-06 5.732365e-06 0.99999713
[313,] 1.817191e-05 3.634382e-05 0.99998183
[314,] 1.324696e-05 2.649392e-05 0.99998675
[315,] 1.410299e-05 2.820598e-05 0.99998590
[316,] 9.733441e-06 1.946688e-05 0.99999027
[317,] 7.496132e-06 1.499226e-05 0.99999250
[318,] 7.810839e-06 1.562168e-05 0.99999219
[319,] 5.959579e-06 1.191916e-05 0.99999404
[320,] 4.093203e-03 8.186407e-03 0.99590680
[321,] 3.097915e-03 6.195831e-03 0.99690208
[322,] 2.263169e-03 4.526338e-03 0.99773683
[323,] 2.374489e-03 4.748978e-03 0.99762551
[324,] 1.779282e-03 3.558563e-03 0.99822072
[325,] 5.295343e-03 1.059069e-02 0.99470466
[326,] 4.345420e-03 8.690840e-03 0.99565458
[327,] 4.329322e-03 8.658645e-03 0.99567068
[328,] 3.204399e-03 6.408798e-03 0.99679560
[329,] 2.537191e-03 5.074382e-03 0.99746281
[330,] 1.764935e-03 3.529871e-03 0.99823506
[331,] 1.212287e-03 2.424574e-03 0.99878771
[332,] 4.200351e-03 8.400702e-03 0.99579965
[333,] 2.957668e-03 5.915336e-03 0.99704233
[334,] 2.805486e-03 5.610972e-03 0.99719451
[335,] 3.589205e-03 7.178411e-03 0.99641079
[336,] 4.714741e-03 9.429482e-03 0.99528526
[337,] 4.203893e-02 8.407786e-02 0.95796107
[338,] 3.148873e-02 6.297746e-02 0.96851127
[339,] 2.448127e-02 4.896254e-02 0.97551873
[340,] 1.963241e-02 3.926482e-02 0.98036759
[341,] 1.475370e-02 2.950741e-02 0.98524630
[342,] 1.146942e-02 2.293884e-02 0.98853058
[343,] 7.877612e-03 1.575522e-02 0.99212239
[344,] 1.998586e-02 3.997173e-02 0.98001414
[345,] 1.926792e-02 3.853584e-02 0.98073208
[346,] 1.624355e-02 3.248710e-02 0.98375645
[347,] 3.879394e-02 7.758787e-02 0.96120606
[348,] 3.679664e-02 7.359329e-02 0.96320336
[349,] 3.010967e-01 6.021934e-01 0.69890332
[350,] 2.760326e-01 5.520653e-01 0.72396736
[351,] 2.187870e-01 4.375740e-01 0.78121299
[352,] 1.805713e-01 3.611426e-01 0.81942868
[353,] 1.411455e-01 2.822911e-01 0.85885446
[354,] 1.035268e-01 2.070537e-01 0.89647316
[355,] 1.058854e-01 2.117709e-01 0.89411457
[356,] 4.039021e-01 8.078043e-01 0.59609785
[357,] 3.167089e-01 6.334178e-01 0.68329108
[358,] 2.309314e-01 4.618627e-01 0.76906864
[359,] 3.689731e-01 7.379461e-01 0.63102693
[360,] 3.536148e-01 7.072296e-01 0.64638522
[361,] 7.752343e-01 4.495315e-01 0.22476573
[362,] 9.149374e-01 1.701252e-01 0.08506258
> postscript(file="/var/www/rcomp/tmp/12zvf1322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/2po2l1322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/3qzl51322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/4aqff1322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/5am5x1322152297.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 = 371
Frequency = 1
1 2 3 4 5
-44.86454221 16.83639820 24.63733861 40.03827902 -38.66078057
6 7 8 9 10
0.44015985 18.04110026 18.44204067 32.14298108 -27.85607851
11 12 13 14 15
-16.45513810 -78.25419768 -47.15325727 9.84768314 11.34862355
16 17 18 19 20
-23.15043604 -44.24949563 -41.44855521 43.95238520 32.85332561
21 22 23 24 25
-14.54573398 22.55520643 -11.24385315 -95.04291274 -13.54197233
26 27 28 29 30
48.55896808 65.55990849 40.96084890 -28.03821068 24.26272973
31 32 33 34 35
71.06367014 13.06461055 44.26555096 -29.33350863 -4.53256821
36 37 38 39 40
-23.23162780 11.86931261 22.27025302 39.27119343 17.57213385
41 42 43 44 45
-35.42692574 6.77401467 27.37495508 -7.22410451 13.97683590
46 47 48 49 50
-22.82222368 20.17871673 -45.42034286 -7.41940245 34.58153796
51 52 53 54 55
17.38247837 6.18341879 -24.21564080 -4.51470039 17.48624002
56 57 58 59 60
19.78718043 25.58812084 -10.61093874 5.69000167 -57.90905792
61 62 63 64 65
17.59188249 14.39282290 7.19376332 23.59470373 -18.80435586
66 67 68 69 70
2.99658455 20.59752496 -9.90153463 5.49940579 -43.69965380
71 72 73 74 75
-33.89871339 -118.89777298 -41.39683257 -4.59589216 19.90504826
76 77 78 79 80
19.60598867 -1.49307092 1.00786949 21.70880990 4.50975032
81 82 83 84 85
51.21069073 -17.78836886 11.41257155 -65.18648804 10.81445237
86 87 88 89 90
28.01539279 14.31633320 48.01727361 -19.48178598 23.91915443
91 92 93 94 95
13.42009484 20.62103526 6.22197567 -24.37708392 -5.67614351
96 97 98 99 100
-45.67520310 -1.77426268 0.52667773 39.32761814 -8.07144145
101 102 103 104 105
-37.17050104 20.33043937 52.53137979 20.23232020 19.03326061
106 107 108 109 110
-49.96579898 -5.26485857 -53.76391816 12.93702226 27.13796267
111 112 113 114 115
20.33890308 0.03984349 -48.75921610 29.64172432 32.54266473
116 117 118 119 120
3.14360514 4.64454555 -65.35451404 -19.65357363 -112.35263321
121 122 123 124 125
-64.35169280 -3.05075239 9.95018802 24.05112843 -38.24793116
126 127 128 129 130
13.35300926 49.95394967 55.75489008 31.65583049 -2.64322910
131 132 133 134 135
-27.44228868 -60.14134827 -1.14040786 40.86053255 74.86147296
136 137 138 139 140
23.96241337 -44.43664621 18.36429420 26.66523461 16.06617502
141 142 143 144 145
-2.73288457 -39.63194416 12.36899626 -58.93006333 24.17087708
146 147 148 149 150
-26.92818251 56.57275790 23.57369832 -78.32536127 29.67557914
151 152 153 154 155
18.17651955 40.47745996 -21.32159963 -42.12065921 -49.41971880
156 157 158 159 160
-85.61877839 -31.01783798 23.98310243 54.48404284 22.48498326
161 162 163 164 165
-63.31407633 36.08686408 52.98780449 41.48874490 17.98968532
166 167 168 169 170
-6.90937427 -9.10843386 -57.10749345 14.89344696 16.69438737
171 172 173 174 175
46.89532779 -8.00373180 -25.80279139 39.89814902 -0.40091057
176 177 178 179 180
39.60002984 23.00097026 -48.29808933 -2.49714892 -85.79620851
181 182 183 184 185
-23.39526810 43.70567232 45.80661273 5.30755314 -59.59150645
186 187 188 189 190
42.30943396 39.41037437 29.41131479 8.51225520 -45.48680439
191 192 193 194 195
7.91413602 -72.08492357 6.61601684 24.51695726 40.31789767
196 197 198 199 200
31.21883808 -91.58022151 78.72071890 13.32165932 29.82259973
201 202 203 204 205
7.32354014 -11.07551945 -8.17457904 -52.27363863 -22.17269821
206 207 208 209 210
48.22824220 21.62918261 28.73012302 -83.36893657 63.73200384
211 212 213 214 215
27.33294426 33.23388467 14.23482508 -16.96423451 10.93670590
216 217 218 219 220
-44.66235368 14.43858673 12.83952714 26.94046755 -5.35859204
221 222 223 224 225
-78.95765163 55.14328879 23.84422920 32.74516961 4.74611002
226 227 228 229 230
-10.45294957 -7.65200916 -48.45106874 -1.55012833 23.95081208
231 232 233 234 235
29.95175249 21.55269290 -116.04636668 38.85457373 31.65551414
236 237 238 239 240
5.65645455 -4.74260504 5.85833537 19.25927579 -34.43978380
241 242 243 244 245
-20.33884339 36.76209702 44.76303743 19.66397784 -130.13508174
246 247 248 249 250
40.76585867 45.46679908 17.46773949 10.66867990 -5.63037968
251 252 253 254 255
16.77056073 -44.72849886 -3.82755845 18.57338196 21.47432237
256 257 258 259 260
25.17526279 -108.02379680 21.77714361 32.17808402 -7.82097557
261 262 263 264 265
12.77996484 13.98090526 9.28184567 -76.91721392 -37.81627351
266 267 268 269 270
7.08466690 19.08560731 17.78654773 -127.51251186 16.88842855
271 272 273 274 275
29.98936896 -6.20969063 4.29124979 -33.80780980 -1.90686939
276 277 278 279 280
-76.80592898 -1.80498857 27.69595184 49.09689226 30.99783267
281 282 283 284 285
-108.60122692 16.99971349 27.90065390 23.10159431 28.00253473
286 287 288 289 290
-23.49652486 13.00441555 -74.19464404 4.50629637 20.70723679
291 292 293 294 295
52.80817720 36.30911761 -107.18994198 26.31099843 32.61193884
296 297 298 299 300
20.91287926 19.81381967 21.41476008 16.01570049 -54.88335910
301 302 303 304 305
-15.98241869 34.31852173 34.81946214 38.52040255 -103.77865704
306 307 308 309 310
30.72228337 35.22322379 5.32416420 41.22510461 -28.27395498
311 312 313 314 315
0.82698543 -93.97207416 -12.17113374 39.52980667 46.43074708
316 317 318 319 320
16.73168749 -122.56737210 13.03356831 38.53450873 -30.66455086
321 322 323 324 325
16.83638955 -63.06267004 -41.06172963 -206.36078921 -11.85984880
326 327 328 329 330
-3.95890839 54.94203202 20.74297243 -93.55608716 37.04485326
331 332 333 334 335
52.34579367 18.44673408 28.84767449 2.34861490 4.64955531
336 337 338 339 340
-96.84950427 15.15143614 51.85237655 64.55331696 59.65425737
341 342 343 344 345
-134.04480221 8.85613820 26.55707861 30.65801902 20.35895943
346 347 348 349 350
-25.14010016 8.36084026 -81.53821933 -25.03727892 56.36366149
351 352 353 354 355
99.86460190 42.76554231 -129.13351727 52.26742314 19.46836355
356 357 358 359 360
33.06930396 22.67024437 -11.42881521 47.67212520 -100.62693439
361 362 363 364 365
16.87400602 27.07494643 80.17588684 24.17682726 -85.82223233
366 367 368 369 370
-10.12129192 51.77964849 14.48058890 34.78152931 -15.81753027
371
-8.51658986
> postscript(file="/var/www/rcomp/tmp/627ut1322152297.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 = 371
Frequency = 1
lag(myerror, k = 1) myerror
0 -44.86454221 NA
1 16.83639820 -44.86454221
2 24.63733861 16.83639820
3 40.03827902 24.63733861
4 -38.66078057 40.03827902
5 0.44015985 -38.66078057
6 18.04110026 0.44015985
7 18.44204067 18.04110026
8 32.14298108 18.44204067
9 -27.85607851 32.14298108
10 -16.45513810 -27.85607851
11 -78.25419768 -16.45513810
12 -47.15325727 -78.25419768
13 9.84768314 -47.15325727
14 11.34862355 9.84768314
15 -23.15043604 11.34862355
16 -44.24949563 -23.15043604
17 -41.44855521 -44.24949563
18 43.95238520 -41.44855521
19 32.85332561 43.95238520
20 -14.54573398 32.85332561
21 22.55520643 -14.54573398
22 -11.24385315 22.55520643
23 -95.04291274 -11.24385315
24 -13.54197233 -95.04291274
25 48.55896808 -13.54197233
26 65.55990849 48.55896808
27 40.96084890 65.55990849
28 -28.03821068 40.96084890
29 24.26272973 -28.03821068
30 71.06367014 24.26272973
31 13.06461055 71.06367014
32 44.26555096 13.06461055
33 -29.33350863 44.26555096
34 -4.53256821 -29.33350863
35 -23.23162780 -4.53256821
36 11.86931261 -23.23162780
37 22.27025302 11.86931261
38 39.27119343 22.27025302
39 17.57213385 39.27119343
40 -35.42692574 17.57213385
41 6.77401467 -35.42692574
42 27.37495508 6.77401467
43 -7.22410451 27.37495508
44 13.97683590 -7.22410451
45 -22.82222368 13.97683590
46 20.17871673 -22.82222368
47 -45.42034286 20.17871673
48 -7.41940245 -45.42034286
49 34.58153796 -7.41940245
50 17.38247837 34.58153796
51 6.18341879 17.38247837
52 -24.21564080 6.18341879
53 -4.51470039 -24.21564080
54 17.48624002 -4.51470039
55 19.78718043 17.48624002
56 25.58812084 19.78718043
57 -10.61093874 25.58812084
58 5.69000167 -10.61093874
59 -57.90905792 5.69000167
60 17.59188249 -57.90905792
61 14.39282290 17.59188249
62 7.19376332 14.39282290
63 23.59470373 7.19376332
64 -18.80435586 23.59470373
65 2.99658455 -18.80435586
66 20.59752496 2.99658455
67 -9.90153463 20.59752496
68 5.49940579 -9.90153463
69 -43.69965380 5.49940579
70 -33.89871339 -43.69965380
71 -118.89777298 -33.89871339
72 -41.39683257 -118.89777298
73 -4.59589216 -41.39683257
74 19.90504826 -4.59589216
75 19.60598867 19.90504826
76 -1.49307092 19.60598867
77 1.00786949 -1.49307092
78 21.70880990 1.00786949
79 4.50975032 21.70880990
80 51.21069073 4.50975032
81 -17.78836886 51.21069073
82 11.41257155 -17.78836886
83 -65.18648804 11.41257155
84 10.81445237 -65.18648804
85 28.01539279 10.81445237
86 14.31633320 28.01539279
87 48.01727361 14.31633320
88 -19.48178598 48.01727361
89 23.91915443 -19.48178598
90 13.42009484 23.91915443
91 20.62103526 13.42009484
92 6.22197567 20.62103526
93 -24.37708392 6.22197567
94 -5.67614351 -24.37708392
95 -45.67520310 -5.67614351
96 -1.77426268 -45.67520310
97 0.52667773 -1.77426268
98 39.32761814 0.52667773
99 -8.07144145 39.32761814
100 -37.17050104 -8.07144145
101 20.33043937 -37.17050104
102 52.53137979 20.33043937
103 20.23232020 52.53137979
104 19.03326061 20.23232020
105 -49.96579898 19.03326061
106 -5.26485857 -49.96579898
107 -53.76391816 -5.26485857
108 12.93702226 -53.76391816
109 27.13796267 12.93702226
110 20.33890308 27.13796267
111 0.03984349 20.33890308
112 -48.75921610 0.03984349
113 29.64172432 -48.75921610
114 32.54266473 29.64172432
115 3.14360514 32.54266473
116 4.64454555 3.14360514
117 -65.35451404 4.64454555
118 -19.65357363 -65.35451404
119 -112.35263321 -19.65357363
120 -64.35169280 -112.35263321
121 -3.05075239 -64.35169280
122 9.95018802 -3.05075239
123 24.05112843 9.95018802
124 -38.24793116 24.05112843
125 13.35300926 -38.24793116
126 49.95394967 13.35300926
127 55.75489008 49.95394967
128 31.65583049 55.75489008
129 -2.64322910 31.65583049
130 -27.44228868 -2.64322910
131 -60.14134827 -27.44228868
132 -1.14040786 -60.14134827
133 40.86053255 -1.14040786
134 74.86147296 40.86053255
135 23.96241337 74.86147296
136 -44.43664621 23.96241337
137 18.36429420 -44.43664621
138 26.66523461 18.36429420
139 16.06617502 26.66523461
140 -2.73288457 16.06617502
141 -39.63194416 -2.73288457
142 12.36899626 -39.63194416
143 -58.93006333 12.36899626
144 24.17087708 -58.93006333
145 -26.92818251 24.17087708
146 56.57275790 -26.92818251
147 23.57369832 56.57275790
148 -78.32536127 23.57369832
149 29.67557914 -78.32536127
150 18.17651955 29.67557914
151 40.47745996 18.17651955
152 -21.32159963 40.47745996
153 -42.12065921 -21.32159963
154 -49.41971880 -42.12065921
155 -85.61877839 -49.41971880
156 -31.01783798 -85.61877839
157 23.98310243 -31.01783798
158 54.48404284 23.98310243
159 22.48498326 54.48404284
160 -63.31407633 22.48498326
161 36.08686408 -63.31407633
162 52.98780449 36.08686408
163 41.48874490 52.98780449
164 17.98968532 41.48874490
165 -6.90937427 17.98968532
166 -9.10843386 -6.90937427
167 -57.10749345 -9.10843386
168 14.89344696 -57.10749345
169 16.69438737 14.89344696
170 46.89532779 16.69438737
171 -8.00373180 46.89532779
172 -25.80279139 -8.00373180
173 39.89814902 -25.80279139
174 -0.40091057 39.89814902
175 39.60002984 -0.40091057
176 23.00097026 39.60002984
177 -48.29808933 23.00097026
178 -2.49714892 -48.29808933
179 -85.79620851 -2.49714892
180 -23.39526810 -85.79620851
181 43.70567232 -23.39526810
182 45.80661273 43.70567232
183 5.30755314 45.80661273
184 -59.59150645 5.30755314
185 42.30943396 -59.59150645
186 39.41037437 42.30943396
187 29.41131479 39.41037437
188 8.51225520 29.41131479
189 -45.48680439 8.51225520
190 7.91413602 -45.48680439
191 -72.08492357 7.91413602
192 6.61601684 -72.08492357
193 24.51695726 6.61601684
194 40.31789767 24.51695726
195 31.21883808 40.31789767
196 -91.58022151 31.21883808
197 78.72071890 -91.58022151
198 13.32165932 78.72071890
199 29.82259973 13.32165932
200 7.32354014 29.82259973
201 -11.07551945 7.32354014
202 -8.17457904 -11.07551945
203 -52.27363863 -8.17457904
204 -22.17269821 -52.27363863
205 48.22824220 -22.17269821
206 21.62918261 48.22824220
207 28.73012302 21.62918261
208 -83.36893657 28.73012302
209 63.73200384 -83.36893657
210 27.33294426 63.73200384
211 33.23388467 27.33294426
212 14.23482508 33.23388467
213 -16.96423451 14.23482508
214 10.93670590 -16.96423451
215 -44.66235368 10.93670590
216 14.43858673 -44.66235368
217 12.83952714 14.43858673
218 26.94046755 12.83952714
219 -5.35859204 26.94046755
220 -78.95765163 -5.35859204
221 55.14328879 -78.95765163
222 23.84422920 55.14328879
223 32.74516961 23.84422920
224 4.74611002 32.74516961
225 -10.45294957 4.74611002
226 -7.65200916 -10.45294957
227 -48.45106874 -7.65200916
228 -1.55012833 -48.45106874
229 23.95081208 -1.55012833
230 29.95175249 23.95081208
231 21.55269290 29.95175249
232 -116.04636668 21.55269290
233 38.85457373 -116.04636668
234 31.65551414 38.85457373
235 5.65645455 31.65551414
236 -4.74260504 5.65645455
237 5.85833537 -4.74260504
238 19.25927579 5.85833537
239 -34.43978380 19.25927579
240 -20.33884339 -34.43978380
241 36.76209702 -20.33884339
242 44.76303743 36.76209702
243 19.66397784 44.76303743
244 -130.13508174 19.66397784
245 40.76585867 -130.13508174
246 45.46679908 40.76585867
247 17.46773949 45.46679908
248 10.66867990 17.46773949
249 -5.63037968 10.66867990
250 16.77056073 -5.63037968
251 -44.72849886 16.77056073
252 -3.82755845 -44.72849886
253 18.57338196 -3.82755845
254 21.47432237 18.57338196
255 25.17526279 21.47432237
256 -108.02379680 25.17526279
257 21.77714361 -108.02379680
258 32.17808402 21.77714361
259 -7.82097557 32.17808402
260 12.77996484 -7.82097557
261 13.98090526 12.77996484
262 9.28184567 13.98090526
263 -76.91721392 9.28184567
264 -37.81627351 -76.91721392
265 7.08466690 -37.81627351
266 19.08560731 7.08466690
267 17.78654773 19.08560731
268 -127.51251186 17.78654773
269 16.88842855 -127.51251186
270 29.98936896 16.88842855
271 -6.20969063 29.98936896
272 4.29124979 -6.20969063
273 -33.80780980 4.29124979
274 -1.90686939 -33.80780980
275 -76.80592898 -1.90686939
276 -1.80498857 -76.80592898
277 27.69595184 -1.80498857
278 49.09689226 27.69595184
279 30.99783267 49.09689226
280 -108.60122692 30.99783267
281 16.99971349 -108.60122692
282 27.90065390 16.99971349
283 23.10159431 27.90065390
284 28.00253473 23.10159431
285 -23.49652486 28.00253473
286 13.00441555 -23.49652486
287 -74.19464404 13.00441555
288 4.50629637 -74.19464404
289 20.70723679 4.50629637
290 52.80817720 20.70723679
291 36.30911761 52.80817720
292 -107.18994198 36.30911761
293 26.31099843 -107.18994198
294 32.61193884 26.31099843
295 20.91287926 32.61193884
296 19.81381967 20.91287926
297 21.41476008 19.81381967
298 16.01570049 21.41476008
299 -54.88335910 16.01570049
300 -15.98241869 -54.88335910
301 34.31852173 -15.98241869
302 34.81946214 34.31852173
303 38.52040255 34.81946214
304 -103.77865704 38.52040255
305 30.72228337 -103.77865704
306 35.22322379 30.72228337
307 5.32416420 35.22322379
308 41.22510461 5.32416420
309 -28.27395498 41.22510461
310 0.82698543 -28.27395498
311 -93.97207416 0.82698543
312 -12.17113374 -93.97207416
313 39.52980667 -12.17113374
314 46.43074708 39.52980667
315 16.73168749 46.43074708
316 -122.56737210 16.73168749
317 13.03356831 -122.56737210
318 38.53450873 13.03356831
319 -30.66455086 38.53450873
320 16.83638955 -30.66455086
321 -63.06267004 16.83638955
322 -41.06172963 -63.06267004
323 -206.36078921 -41.06172963
324 -11.85984880 -206.36078921
325 -3.95890839 -11.85984880
326 54.94203202 -3.95890839
327 20.74297243 54.94203202
328 -93.55608716 20.74297243
329 37.04485326 -93.55608716
330 52.34579367 37.04485326
331 18.44673408 52.34579367
332 28.84767449 18.44673408
333 2.34861490 28.84767449
334 4.64955531 2.34861490
335 -96.84950427 4.64955531
336 15.15143614 -96.84950427
337 51.85237655 15.15143614
338 64.55331696 51.85237655
339 59.65425737 64.55331696
340 -134.04480221 59.65425737
341 8.85613820 -134.04480221
342 26.55707861 8.85613820
343 30.65801902 26.55707861
344 20.35895943 30.65801902
345 -25.14010016 20.35895943
346 8.36084026 -25.14010016
347 -81.53821933 8.36084026
348 -25.03727892 -81.53821933
349 56.36366149 -25.03727892
350 99.86460190 56.36366149
351 42.76554231 99.86460190
352 -129.13351727 42.76554231
353 52.26742314 -129.13351727
354 19.46836355 52.26742314
355 33.06930396 19.46836355
356 22.67024437 33.06930396
357 -11.42881521 22.67024437
358 47.67212520 -11.42881521
359 -100.62693439 47.67212520
360 16.87400602 -100.62693439
361 27.07494643 16.87400602
362 80.17588684 27.07494643
363 24.17682726 80.17588684
364 -85.82223233 24.17682726
365 -10.12129192 -85.82223233
366 51.77964849 -10.12129192
367 14.48058890 51.77964849
368 34.78152931 14.48058890
369 -15.81753027 34.78152931
370 -8.51658986 -15.81753027
371 NA -8.51658986
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 16.83639820 -44.86454221
[2,] 24.63733861 16.83639820
[3,] 40.03827902 24.63733861
[4,] -38.66078057 40.03827902
[5,] 0.44015985 -38.66078057
[6,] 18.04110026 0.44015985
[7,] 18.44204067 18.04110026
[8,] 32.14298108 18.44204067
[9,] -27.85607851 32.14298108
[10,] -16.45513810 -27.85607851
[11,] -78.25419768 -16.45513810
[12,] -47.15325727 -78.25419768
[13,] 9.84768314 -47.15325727
[14,] 11.34862355 9.84768314
[15,] -23.15043604 11.34862355
[16,] -44.24949563 -23.15043604
[17,] -41.44855521 -44.24949563
[18,] 43.95238520 -41.44855521
[19,] 32.85332561 43.95238520
[20,] -14.54573398 32.85332561
[21,] 22.55520643 -14.54573398
[22,] -11.24385315 22.55520643
[23,] -95.04291274 -11.24385315
[24,] -13.54197233 -95.04291274
[25,] 48.55896808 -13.54197233
[26,] 65.55990849 48.55896808
[27,] 40.96084890 65.55990849
[28,] -28.03821068 40.96084890
[29,] 24.26272973 -28.03821068
[30,] 71.06367014 24.26272973
[31,] 13.06461055 71.06367014
[32,] 44.26555096 13.06461055
[33,] -29.33350863 44.26555096
[34,] -4.53256821 -29.33350863
[35,] -23.23162780 -4.53256821
[36,] 11.86931261 -23.23162780
[37,] 22.27025302 11.86931261
[38,] 39.27119343 22.27025302
[39,] 17.57213385 39.27119343
[40,] -35.42692574 17.57213385
[41,] 6.77401467 -35.42692574
[42,] 27.37495508 6.77401467
[43,] -7.22410451 27.37495508
[44,] 13.97683590 -7.22410451
[45,] -22.82222368 13.97683590
[46,] 20.17871673 -22.82222368
[47,] -45.42034286 20.17871673
[48,] -7.41940245 -45.42034286
[49,] 34.58153796 -7.41940245
[50,] 17.38247837 34.58153796
[51,] 6.18341879 17.38247837
[52,] -24.21564080 6.18341879
[53,] -4.51470039 -24.21564080
[54,] 17.48624002 -4.51470039
[55,] 19.78718043 17.48624002
[56,] 25.58812084 19.78718043
[57,] -10.61093874 25.58812084
[58,] 5.69000167 -10.61093874
[59,] -57.90905792 5.69000167
[60,] 17.59188249 -57.90905792
[61,] 14.39282290 17.59188249
[62,] 7.19376332 14.39282290
[63,] 23.59470373 7.19376332
[64,] -18.80435586 23.59470373
[65,] 2.99658455 -18.80435586
[66,] 20.59752496 2.99658455
[67,] -9.90153463 20.59752496
[68,] 5.49940579 -9.90153463
[69,] -43.69965380 5.49940579
[70,] -33.89871339 -43.69965380
[71,] -118.89777298 -33.89871339
[72,] -41.39683257 -118.89777298
[73,] -4.59589216 -41.39683257
[74,] 19.90504826 -4.59589216
[75,] 19.60598867 19.90504826
[76,] -1.49307092 19.60598867
[77,] 1.00786949 -1.49307092
[78,] 21.70880990 1.00786949
[79,] 4.50975032 21.70880990
[80,] 51.21069073 4.50975032
[81,] -17.78836886 51.21069073
[82,] 11.41257155 -17.78836886
[83,] -65.18648804 11.41257155
[84,] 10.81445237 -65.18648804
[85,] 28.01539279 10.81445237
[86,] 14.31633320 28.01539279
[87,] 48.01727361 14.31633320
[88,] -19.48178598 48.01727361
[89,] 23.91915443 -19.48178598
[90,] 13.42009484 23.91915443
[91,] 20.62103526 13.42009484
[92,] 6.22197567 20.62103526
[93,] -24.37708392 6.22197567
[94,] -5.67614351 -24.37708392
[95,] -45.67520310 -5.67614351
[96,] -1.77426268 -45.67520310
[97,] 0.52667773 -1.77426268
[98,] 39.32761814 0.52667773
[99,] -8.07144145 39.32761814
[100,] -37.17050104 -8.07144145
[101,] 20.33043937 -37.17050104
[102,] 52.53137979 20.33043937
[103,] 20.23232020 52.53137979
[104,] 19.03326061 20.23232020
[105,] -49.96579898 19.03326061
[106,] -5.26485857 -49.96579898
[107,] -53.76391816 -5.26485857
[108,] 12.93702226 -53.76391816
[109,] 27.13796267 12.93702226
[110,] 20.33890308 27.13796267
[111,] 0.03984349 20.33890308
[112,] -48.75921610 0.03984349
[113,] 29.64172432 -48.75921610
[114,] 32.54266473 29.64172432
[115,] 3.14360514 32.54266473
[116,] 4.64454555 3.14360514
[117,] -65.35451404 4.64454555
[118,] -19.65357363 -65.35451404
[119,] -112.35263321 -19.65357363
[120,] -64.35169280 -112.35263321
[121,] -3.05075239 -64.35169280
[122,] 9.95018802 -3.05075239
[123,] 24.05112843 9.95018802
[124,] -38.24793116 24.05112843
[125,] 13.35300926 -38.24793116
[126,] 49.95394967 13.35300926
[127,] 55.75489008 49.95394967
[128,] 31.65583049 55.75489008
[129,] -2.64322910 31.65583049
[130,] -27.44228868 -2.64322910
[131,] -60.14134827 -27.44228868
[132,] -1.14040786 -60.14134827
[133,] 40.86053255 -1.14040786
[134,] 74.86147296 40.86053255
[135,] 23.96241337 74.86147296
[136,] -44.43664621 23.96241337
[137,] 18.36429420 -44.43664621
[138,] 26.66523461 18.36429420
[139,] 16.06617502 26.66523461
[140,] -2.73288457 16.06617502
[141,] -39.63194416 -2.73288457
[142,] 12.36899626 -39.63194416
[143,] -58.93006333 12.36899626
[144,] 24.17087708 -58.93006333
[145,] -26.92818251 24.17087708
[146,] 56.57275790 -26.92818251
[147,] 23.57369832 56.57275790
[148,] -78.32536127 23.57369832
[149,] 29.67557914 -78.32536127
[150,] 18.17651955 29.67557914
[151,] 40.47745996 18.17651955
[152,] -21.32159963 40.47745996
[153,] -42.12065921 -21.32159963
[154,] -49.41971880 -42.12065921
[155,] -85.61877839 -49.41971880
[156,] -31.01783798 -85.61877839
[157,] 23.98310243 -31.01783798
[158,] 54.48404284 23.98310243
[159,] 22.48498326 54.48404284
[160,] -63.31407633 22.48498326
[161,] 36.08686408 -63.31407633
[162,] 52.98780449 36.08686408
[163,] 41.48874490 52.98780449
[164,] 17.98968532 41.48874490
[165,] -6.90937427 17.98968532
[166,] -9.10843386 -6.90937427
[167,] -57.10749345 -9.10843386
[168,] 14.89344696 -57.10749345
[169,] 16.69438737 14.89344696
[170,] 46.89532779 16.69438737
[171,] -8.00373180 46.89532779
[172,] -25.80279139 -8.00373180
[173,] 39.89814902 -25.80279139
[174,] -0.40091057 39.89814902
[175,] 39.60002984 -0.40091057
[176,] 23.00097026 39.60002984
[177,] -48.29808933 23.00097026
[178,] -2.49714892 -48.29808933
[179,] -85.79620851 -2.49714892
[180,] -23.39526810 -85.79620851
[181,] 43.70567232 -23.39526810
[182,] 45.80661273 43.70567232
[183,] 5.30755314 45.80661273
[184,] -59.59150645 5.30755314
[185,] 42.30943396 -59.59150645
[186,] 39.41037437 42.30943396
[187,] 29.41131479 39.41037437
[188,] 8.51225520 29.41131479
[189,] -45.48680439 8.51225520
[190,] 7.91413602 -45.48680439
[191,] -72.08492357 7.91413602
[192,] 6.61601684 -72.08492357
[193,] 24.51695726 6.61601684
[194,] 40.31789767 24.51695726
[195,] 31.21883808 40.31789767
[196,] -91.58022151 31.21883808
[197,] 78.72071890 -91.58022151
[198,] 13.32165932 78.72071890
[199,] 29.82259973 13.32165932
[200,] 7.32354014 29.82259973
[201,] -11.07551945 7.32354014
[202,] -8.17457904 -11.07551945
[203,] -52.27363863 -8.17457904
[204,] -22.17269821 -52.27363863
[205,] 48.22824220 -22.17269821
[206,] 21.62918261 48.22824220
[207,] 28.73012302 21.62918261
[208,] -83.36893657 28.73012302
[209,] 63.73200384 -83.36893657
[210,] 27.33294426 63.73200384
[211,] 33.23388467 27.33294426
[212,] 14.23482508 33.23388467
[213,] -16.96423451 14.23482508
[214,] 10.93670590 -16.96423451
[215,] -44.66235368 10.93670590
[216,] 14.43858673 -44.66235368
[217,] 12.83952714 14.43858673
[218,] 26.94046755 12.83952714
[219,] -5.35859204 26.94046755
[220,] -78.95765163 -5.35859204
[221,] 55.14328879 -78.95765163
[222,] 23.84422920 55.14328879
[223,] 32.74516961 23.84422920
[224,] 4.74611002 32.74516961
[225,] -10.45294957 4.74611002
[226,] -7.65200916 -10.45294957
[227,] -48.45106874 -7.65200916
[228,] -1.55012833 -48.45106874
[229,] 23.95081208 -1.55012833
[230,] 29.95175249 23.95081208
[231,] 21.55269290 29.95175249
[232,] -116.04636668 21.55269290
[233,] 38.85457373 -116.04636668
[234,] 31.65551414 38.85457373
[235,] 5.65645455 31.65551414
[236,] -4.74260504 5.65645455
[237,] 5.85833537 -4.74260504
[238,] 19.25927579 5.85833537
[239,] -34.43978380 19.25927579
[240,] -20.33884339 -34.43978380
[241,] 36.76209702 -20.33884339
[242,] 44.76303743 36.76209702
[243,] 19.66397784 44.76303743
[244,] -130.13508174 19.66397784
[245,] 40.76585867 -130.13508174
[246,] 45.46679908 40.76585867
[247,] 17.46773949 45.46679908
[248,] 10.66867990 17.46773949
[249,] -5.63037968 10.66867990
[250,] 16.77056073 -5.63037968
[251,] -44.72849886 16.77056073
[252,] -3.82755845 -44.72849886
[253,] 18.57338196 -3.82755845
[254,] 21.47432237 18.57338196
[255,] 25.17526279 21.47432237
[256,] -108.02379680 25.17526279
[257,] 21.77714361 -108.02379680
[258,] 32.17808402 21.77714361
[259,] -7.82097557 32.17808402
[260,] 12.77996484 -7.82097557
[261,] 13.98090526 12.77996484
[262,] 9.28184567 13.98090526
[263,] -76.91721392 9.28184567
[264,] -37.81627351 -76.91721392
[265,] 7.08466690 -37.81627351
[266,] 19.08560731 7.08466690
[267,] 17.78654773 19.08560731
[268,] -127.51251186 17.78654773
[269,] 16.88842855 -127.51251186
[270,] 29.98936896 16.88842855
[271,] -6.20969063 29.98936896
[272,] 4.29124979 -6.20969063
[273,] -33.80780980 4.29124979
[274,] -1.90686939 -33.80780980
[275,] -76.80592898 -1.90686939
[276,] -1.80498857 -76.80592898
[277,] 27.69595184 -1.80498857
[278,] 49.09689226 27.69595184
[279,] 30.99783267 49.09689226
[280,] -108.60122692 30.99783267
[281,] 16.99971349 -108.60122692
[282,] 27.90065390 16.99971349
[283,] 23.10159431 27.90065390
[284,] 28.00253473 23.10159431
[285,] -23.49652486 28.00253473
[286,] 13.00441555 -23.49652486
[287,] -74.19464404 13.00441555
[288,] 4.50629637 -74.19464404
[289,] 20.70723679 4.50629637
[290,] 52.80817720 20.70723679
[291,] 36.30911761 52.80817720
[292,] -107.18994198 36.30911761
[293,] 26.31099843 -107.18994198
[294,] 32.61193884 26.31099843
[295,] 20.91287926 32.61193884
[296,] 19.81381967 20.91287926
[297,] 21.41476008 19.81381967
[298,] 16.01570049 21.41476008
[299,] -54.88335910 16.01570049
[300,] -15.98241869 -54.88335910
[301,] 34.31852173 -15.98241869
[302,] 34.81946214 34.31852173
[303,] 38.52040255 34.81946214
[304,] -103.77865704 38.52040255
[305,] 30.72228337 -103.77865704
[306,] 35.22322379 30.72228337
[307,] 5.32416420 35.22322379
[308,] 41.22510461 5.32416420
[309,] -28.27395498 41.22510461
[310,] 0.82698543 -28.27395498
[311,] -93.97207416 0.82698543
[312,] -12.17113374 -93.97207416
[313,] 39.52980667 -12.17113374
[314,] 46.43074708 39.52980667
[315,] 16.73168749 46.43074708
[316,] -122.56737210 16.73168749
[317,] 13.03356831 -122.56737210
[318,] 38.53450873 13.03356831
[319,] -30.66455086 38.53450873
[320,] 16.83638955 -30.66455086
[321,] -63.06267004 16.83638955
[322,] -41.06172963 -63.06267004
[323,] -206.36078921 -41.06172963
[324,] -11.85984880 -206.36078921
[325,] -3.95890839 -11.85984880
[326,] 54.94203202 -3.95890839
[327,] 20.74297243 54.94203202
[328,] -93.55608716 20.74297243
[329,] 37.04485326 -93.55608716
[330,] 52.34579367 37.04485326
[331,] 18.44673408 52.34579367
[332,] 28.84767449 18.44673408
[333,] 2.34861490 28.84767449
[334,] 4.64955531 2.34861490
[335,] -96.84950427 4.64955531
[336,] 15.15143614 -96.84950427
[337,] 51.85237655 15.15143614
[338,] 64.55331696 51.85237655
[339,] 59.65425737 64.55331696
[340,] -134.04480221 59.65425737
[341,] 8.85613820 -134.04480221
[342,] 26.55707861 8.85613820
[343,] 30.65801902 26.55707861
[344,] 20.35895943 30.65801902
[345,] -25.14010016 20.35895943
[346,] 8.36084026 -25.14010016
[347,] -81.53821933 8.36084026
[348,] -25.03727892 -81.53821933
[349,] 56.36366149 -25.03727892
[350,] 99.86460190 56.36366149
[351,] 42.76554231 99.86460190
[352,] -129.13351727 42.76554231
[353,] 52.26742314 -129.13351727
[354,] 19.46836355 52.26742314
[355,] 33.06930396 19.46836355
[356,] 22.67024437 33.06930396
[357,] -11.42881521 22.67024437
[358,] 47.67212520 -11.42881521
[359,] -100.62693439 47.67212520
[360,] 16.87400602 -100.62693439
[361,] 27.07494643 16.87400602
[362,] 80.17588684 27.07494643
[363,] 24.17682726 80.17588684
[364,] -85.82223233 24.17682726
[365,] -10.12129192 -85.82223233
[366,] 51.77964849 -10.12129192
[367,] 14.48058890 51.77964849
[368,] 34.78152931 14.48058890
[369,] -15.81753027 34.78152931
[370,] -8.51658986 -15.81753027
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 16.83639820 -44.86454221
2 24.63733861 16.83639820
3 40.03827902 24.63733861
4 -38.66078057 40.03827902
5 0.44015985 -38.66078057
6 18.04110026 0.44015985
7 18.44204067 18.04110026
8 32.14298108 18.44204067
9 -27.85607851 32.14298108
10 -16.45513810 -27.85607851
11 -78.25419768 -16.45513810
12 -47.15325727 -78.25419768
13 9.84768314 -47.15325727
14 11.34862355 9.84768314
15 -23.15043604 11.34862355
16 -44.24949563 -23.15043604
17 -41.44855521 -44.24949563
18 43.95238520 -41.44855521
19 32.85332561 43.95238520
20 -14.54573398 32.85332561
21 22.55520643 -14.54573398
22 -11.24385315 22.55520643
23 -95.04291274 -11.24385315
24 -13.54197233 -95.04291274
25 48.55896808 -13.54197233
26 65.55990849 48.55896808
27 40.96084890 65.55990849
28 -28.03821068 40.96084890
29 24.26272973 -28.03821068
30 71.06367014 24.26272973
31 13.06461055 71.06367014
32 44.26555096 13.06461055
33 -29.33350863 44.26555096
34 -4.53256821 -29.33350863
35 -23.23162780 -4.53256821
36 11.86931261 -23.23162780
37 22.27025302 11.86931261
38 39.27119343 22.27025302
39 17.57213385 39.27119343
40 -35.42692574 17.57213385
41 6.77401467 -35.42692574
42 27.37495508 6.77401467
43 -7.22410451 27.37495508
44 13.97683590 -7.22410451
45 -22.82222368 13.97683590
46 20.17871673 -22.82222368
47 -45.42034286 20.17871673
48 -7.41940245 -45.42034286
49 34.58153796 -7.41940245
50 17.38247837 34.58153796
51 6.18341879 17.38247837
52 -24.21564080 6.18341879
53 -4.51470039 -24.21564080
54 17.48624002 -4.51470039
55 19.78718043 17.48624002
56 25.58812084 19.78718043
57 -10.61093874 25.58812084
58 5.69000167 -10.61093874
59 -57.90905792 5.69000167
60 17.59188249 -57.90905792
61 14.39282290 17.59188249
62 7.19376332 14.39282290
63 23.59470373 7.19376332
64 -18.80435586 23.59470373
65 2.99658455 -18.80435586
66 20.59752496 2.99658455
67 -9.90153463 20.59752496
68 5.49940579 -9.90153463
69 -43.69965380 5.49940579
70 -33.89871339 -43.69965380
71 -118.89777298 -33.89871339
72 -41.39683257 -118.89777298
73 -4.59589216 -41.39683257
74 19.90504826 -4.59589216
75 19.60598867 19.90504826
76 -1.49307092 19.60598867
77 1.00786949 -1.49307092
78 21.70880990 1.00786949
79 4.50975032 21.70880990
80 51.21069073 4.50975032
81 -17.78836886 51.21069073
82 11.41257155 -17.78836886
83 -65.18648804 11.41257155
84 10.81445237 -65.18648804
85 28.01539279 10.81445237
86 14.31633320 28.01539279
87 48.01727361 14.31633320
88 -19.48178598 48.01727361
89 23.91915443 -19.48178598
90 13.42009484 23.91915443
91 20.62103526 13.42009484
92 6.22197567 20.62103526
93 -24.37708392 6.22197567
94 -5.67614351 -24.37708392
95 -45.67520310 -5.67614351
96 -1.77426268 -45.67520310
97 0.52667773 -1.77426268
98 39.32761814 0.52667773
99 -8.07144145 39.32761814
100 -37.17050104 -8.07144145
101 20.33043937 -37.17050104
102 52.53137979 20.33043937
103 20.23232020 52.53137979
104 19.03326061 20.23232020
105 -49.96579898 19.03326061
106 -5.26485857 -49.96579898
107 -53.76391816 -5.26485857
108 12.93702226 -53.76391816
109 27.13796267 12.93702226
110 20.33890308 27.13796267
111 0.03984349 20.33890308
112 -48.75921610 0.03984349
113 29.64172432 -48.75921610
114 32.54266473 29.64172432
115 3.14360514 32.54266473
116 4.64454555 3.14360514
117 -65.35451404 4.64454555
118 -19.65357363 -65.35451404
119 -112.35263321 -19.65357363
120 -64.35169280 -112.35263321
121 -3.05075239 -64.35169280
122 9.95018802 -3.05075239
123 24.05112843 9.95018802
124 -38.24793116 24.05112843
125 13.35300926 -38.24793116
126 49.95394967 13.35300926
127 55.75489008 49.95394967
128 31.65583049 55.75489008
129 -2.64322910 31.65583049
130 -27.44228868 -2.64322910
131 -60.14134827 -27.44228868
132 -1.14040786 -60.14134827
133 40.86053255 -1.14040786
134 74.86147296 40.86053255
135 23.96241337 74.86147296
136 -44.43664621 23.96241337
137 18.36429420 -44.43664621
138 26.66523461 18.36429420
139 16.06617502 26.66523461
140 -2.73288457 16.06617502
141 -39.63194416 -2.73288457
142 12.36899626 -39.63194416
143 -58.93006333 12.36899626
144 24.17087708 -58.93006333
145 -26.92818251 24.17087708
146 56.57275790 -26.92818251
147 23.57369832 56.57275790
148 -78.32536127 23.57369832
149 29.67557914 -78.32536127
150 18.17651955 29.67557914
151 40.47745996 18.17651955
152 -21.32159963 40.47745996
153 -42.12065921 -21.32159963
154 -49.41971880 -42.12065921
155 -85.61877839 -49.41971880
156 -31.01783798 -85.61877839
157 23.98310243 -31.01783798
158 54.48404284 23.98310243
159 22.48498326 54.48404284
160 -63.31407633 22.48498326
161 36.08686408 -63.31407633
162 52.98780449 36.08686408
163 41.48874490 52.98780449
164 17.98968532 41.48874490
165 -6.90937427 17.98968532
166 -9.10843386 -6.90937427
167 -57.10749345 -9.10843386
168 14.89344696 -57.10749345
169 16.69438737 14.89344696
170 46.89532779 16.69438737
171 -8.00373180 46.89532779
172 -25.80279139 -8.00373180
173 39.89814902 -25.80279139
174 -0.40091057 39.89814902
175 39.60002984 -0.40091057
176 23.00097026 39.60002984
177 -48.29808933 23.00097026
178 -2.49714892 -48.29808933
179 -85.79620851 -2.49714892
180 -23.39526810 -85.79620851
181 43.70567232 -23.39526810
182 45.80661273 43.70567232
183 5.30755314 45.80661273
184 -59.59150645 5.30755314
185 42.30943396 -59.59150645
186 39.41037437 42.30943396
187 29.41131479 39.41037437
188 8.51225520 29.41131479
189 -45.48680439 8.51225520
190 7.91413602 -45.48680439
191 -72.08492357 7.91413602
192 6.61601684 -72.08492357
193 24.51695726 6.61601684
194 40.31789767 24.51695726
195 31.21883808 40.31789767
196 -91.58022151 31.21883808
197 78.72071890 -91.58022151
198 13.32165932 78.72071890
199 29.82259973 13.32165932
200 7.32354014 29.82259973
201 -11.07551945 7.32354014
202 -8.17457904 -11.07551945
203 -52.27363863 -8.17457904
204 -22.17269821 -52.27363863
205 48.22824220 -22.17269821
206 21.62918261 48.22824220
207 28.73012302 21.62918261
208 -83.36893657 28.73012302
209 63.73200384 -83.36893657
210 27.33294426 63.73200384
211 33.23388467 27.33294426
212 14.23482508 33.23388467
213 -16.96423451 14.23482508
214 10.93670590 -16.96423451
215 -44.66235368 10.93670590
216 14.43858673 -44.66235368
217 12.83952714 14.43858673
218 26.94046755 12.83952714
219 -5.35859204 26.94046755
220 -78.95765163 -5.35859204
221 55.14328879 -78.95765163
222 23.84422920 55.14328879
223 32.74516961 23.84422920
224 4.74611002 32.74516961
225 -10.45294957 4.74611002
226 -7.65200916 -10.45294957
227 -48.45106874 -7.65200916
228 -1.55012833 -48.45106874
229 23.95081208 -1.55012833
230 29.95175249 23.95081208
231 21.55269290 29.95175249
232 -116.04636668 21.55269290
233 38.85457373 -116.04636668
234 31.65551414 38.85457373
235 5.65645455 31.65551414
236 -4.74260504 5.65645455
237 5.85833537 -4.74260504
238 19.25927579 5.85833537
239 -34.43978380 19.25927579
240 -20.33884339 -34.43978380
241 36.76209702 -20.33884339
242 44.76303743 36.76209702
243 19.66397784 44.76303743
244 -130.13508174 19.66397784
245 40.76585867 -130.13508174
246 45.46679908 40.76585867
247 17.46773949 45.46679908
248 10.66867990 17.46773949
249 -5.63037968 10.66867990
250 16.77056073 -5.63037968
251 -44.72849886 16.77056073
252 -3.82755845 -44.72849886
253 18.57338196 -3.82755845
254 21.47432237 18.57338196
255 25.17526279 21.47432237
256 -108.02379680 25.17526279
257 21.77714361 -108.02379680
258 32.17808402 21.77714361
259 -7.82097557 32.17808402
260 12.77996484 -7.82097557
261 13.98090526 12.77996484
262 9.28184567 13.98090526
263 -76.91721392 9.28184567
264 -37.81627351 -76.91721392
265 7.08466690 -37.81627351
266 19.08560731 7.08466690
267 17.78654773 19.08560731
268 -127.51251186 17.78654773
269 16.88842855 -127.51251186
270 29.98936896 16.88842855
271 -6.20969063 29.98936896
272 4.29124979 -6.20969063
273 -33.80780980 4.29124979
274 -1.90686939 -33.80780980
275 -76.80592898 -1.90686939
276 -1.80498857 -76.80592898
277 27.69595184 -1.80498857
278 49.09689226 27.69595184
279 30.99783267 49.09689226
280 -108.60122692 30.99783267
281 16.99971349 -108.60122692
282 27.90065390 16.99971349
283 23.10159431 27.90065390
284 28.00253473 23.10159431
285 -23.49652486 28.00253473
286 13.00441555 -23.49652486
287 -74.19464404 13.00441555
288 4.50629637 -74.19464404
289 20.70723679 4.50629637
290 52.80817720 20.70723679
291 36.30911761 52.80817720
292 -107.18994198 36.30911761
293 26.31099843 -107.18994198
294 32.61193884 26.31099843
295 20.91287926 32.61193884
296 19.81381967 20.91287926
297 21.41476008 19.81381967
298 16.01570049 21.41476008
299 -54.88335910 16.01570049
300 -15.98241869 -54.88335910
301 34.31852173 -15.98241869
302 34.81946214 34.31852173
303 38.52040255 34.81946214
304 -103.77865704 38.52040255
305 30.72228337 -103.77865704
306 35.22322379 30.72228337
307 5.32416420 35.22322379
308 41.22510461 5.32416420
309 -28.27395498 41.22510461
310 0.82698543 -28.27395498
311 -93.97207416 0.82698543
312 -12.17113374 -93.97207416
313 39.52980667 -12.17113374
314 46.43074708 39.52980667
315 16.73168749 46.43074708
316 -122.56737210 16.73168749
317 13.03356831 -122.56737210
318 38.53450873 13.03356831
319 -30.66455086 38.53450873
320 16.83638955 -30.66455086
321 -63.06267004 16.83638955
322 -41.06172963 -63.06267004
323 -206.36078921 -41.06172963
324 -11.85984880 -206.36078921
325 -3.95890839 -11.85984880
326 54.94203202 -3.95890839
327 20.74297243 54.94203202
328 -93.55608716 20.74297243
329 37.04485326 -93.55608716
330 52.34579367 37.04485326
331 18.44673408 52.34579367
332 28.84767449 18.44673408
333 2.34861490 28.84767449
334 4.64955531 2.34861490
335 -96.84950427 4.64955531
336 15.15143614 -96.84950427
337 51.85237655 15.15143614
338 64.55331696 51.85237655
339 59.65425737 64.55331696
340 -134.04480221 59.65425737
341 8.85613820 -134.04480221
342 26.55707861 8.85613820
343 30.65801902 26.55707861
344 20.35895943 30.65801902
345 -25.14010016 20.35895943
346 8.36084026 -25.14010016
347 -81.53821933 8.36084026
348 -25.03727892 -81.53821933
349 56.36366149 -25.03727892
350 99.86460190 56.36366149
351 42.76554231 99.86460190
352 -129.13351727 42.76554231
353 52.26742314 -129.13351727
354 19.46836355 52.26742314
355 33.06930396 19.46836355
356 22.67024437 33.06930396
357 -11.42881521 22.67024437
358 47.67212520 -11.42881521
359 -100.62693439 47.67212520
360 16.87400602 -100.62693439
361 27.07494643 16.87400602
362 80.17588684 27.07494643
363 24.17682726 80.17588684
364 -85.82223233 24.17682726
365 -10.12129192 -85.82223233
366 51.77964849 -10.12129192
367 14.48058890 51.77964849
368 34.78152931 14.48058890
369 -15.81753027 34.78152931
370 -8.51658986 -15.81753027
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/79oht1322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/87fba1322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/9f6q41322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/rcomp/tmp/10vsql1322152297.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/11rbr21322152297.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/1226hd1322152297.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/13tmaw1322152297.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/14t4p31322152297.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/15s7w01322152297.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/16m8b21322152297.tab")
+ }
>
> try(system("convert tmp/12zvf1322152297.ps tmp/12zvf1322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/2po2l1322152297.ps tmp/2po2l1322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/3qzl51322152297.ps tmp/3qzl51322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/4aqff1322152297.ps tmp/4aqff1322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/5am5x1322152297.ps tmp/5am5x1322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/627ut1322152297.ps tmp/627ut1322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/79oht1322152297.ps tmp/79oht1322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/87fba1322152297.ps tmp/87fba1322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/9f6q41322152297.ps tmp/9f6q41322152297.png",intern=TRUE))
character(0)
> try(system("convert tmp/10vsql1322152297.ps tmp/10vsql1322152297.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.152 0.756 14.401