R version 2.10.1 (2009-12-14)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
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(349
+ ,336
+ ,331
+ ,327
+ ,323
+ ,322
+ ,385
+ ,405
+ ,412
+ ,411
+ ,410
+ ,415
+ ,414
+ ,411
+ ,408
+ ,410
+ ,411
+ ,416
+ ,479
+ ,498
+ ,502
+ ,498
+ ,499
+ ,506
+ ,510
+ ,509
+ ,502
+ ,495
+ ,490
+ ,490
+ ,553
+ ,570
+ ,573
+ ,572
+ ,575
+ ,580
+ ,580
+ ,574
+ ,563
+ ,556
+ ,546
+ ,545
+ ,605
+ ,628
+ ,631
+ ,626
+ ,614
+ ,606
+ ,602
+ ,589
+ ,574
+ ,558
+ ,552
+ ,546
+ ,607
+ ,636
+ ,631
+ ,623
+ ,618
+ ,605
+ ,619
+ ,596
+ ,570
+ ,546
+ ,528
+ ,506
+ ,555
+ ,568
+ ,564
+ ,553
+ ,541
+ ,542
+ ,540
+ ,521
+ ,505
+ ,491
+ ,482
+ ,478
+ ,523
+ ,531
+ ,532
+ ,540
+ ,525
+ ,533
+ ,531
+ ,508
+ ,495
+ ,482
+ ,470
+ ,466
+ ,515
+ ,518
+ ,516
+ ,511
+ ,500
+ ,498
+ ,494
+ ,476
+ ,458
+ ,443
+ ,430
+ ,424
+ ,476
+ ,481
+ ,470
+ ,460
+ ,451
+ ,450
+ ,444
+ ,429
+ ,421
+ ,400
+ ,389
+ ,384
+ ,432
+ ,446
+ ,431
+ ,423
+ ,416
+ ,416
+ ,413
+ ,399
+ ,386
+ ,374
+ ,365
+ ,365
+ ,418
+ ,428
+ ,424
+ ,421
+ ,417
+ ,423
+ ,423
+ ,419
+ ,406
+ ,398
+ ,390
+ ,391
+ ,444
+ ,460
+ ,455
+ ,456
+ ,452
+ ,459
+ ,461
+ ,451
+ ,443
+ ,439
+ ,430
+ ,436
+ ,488
+ ,506
+ ,502
+ ,501
+ ,501
+ ,515
+ ,521
+ ,520
+ ,512
+ ,509
+ ,505
+ ,511
+ ,570
+ ,592
+ ,594
+ ,586
+ ,586
+ ,592
+ ,594
+ ,594
+ ,586
+ ,586
+ ,572
+ ,572
+ ,563
+ ,563
+ ,555
+ ,555
+ ,554
+ ,554
+ ,601
+ ,601
+ ,622
+ ,622
+ ,617
+ ,617
+ ,606
+ ,606
+ ,595
+ ,595
+ ,599
+ ,599
+ ,600
+ ,600
+ ,592
+ ,592
+ ,575
+ ,575
+ ,567
+ ,567
+ ,555
+ ,555
+ ,555
+ ,555
+ ,608
+ ,608
+ ,631
+ ,631
+ ,629
+ ,629
+ ,624
+ ,624
+ ,610
+ ,610
+ ,616
+ ,616
+ ,621
+ ,621
+ ,604
+ ,604
+ ,584
+ ,584
+ ,574
+ ,574
+ ,555
+ ,555
+ ,545
+ ,545
+ ,599
+ ,599
+ ,620
+ ,620
+ ,608
+ ,608
+ ,590
+ ,590
+ ,579
+ ,579
+ ,580
+ ,580
+ ,579
+ ,579
+ ,572
+ ,572
+ ,560
+ ,560
+ ,551
+ ,551
+ ,537
+ ,537
+ ,541
+ ,541
+ ,588
+ ,588
+ ,607
+ ,607
+ ,599
+ ,599
+ ,578
+ ,578
+ ,563
+ ,563
+ ,566
+ ,566
+ ,561
+ ,561
+ ,554
+ ,554
+ ,540
+ ,540
+ ,526
+ ,526
+ ,512
+ ,512
+ ,505
+ ,505
+ ,554
+ ,554
+ ,584
+ ,584
+ ,569
+ ,569
+ ,540
+ ,540
+ ,522
+ ,522
+ ,526
+ ,526
+ ,527
+ ,527
+ ,516
+ ,516
+ ,503
+ ,503
+ ,489
+ ,489
+ ,479
+ ,479
+ ,475
+ ,475
+ ,524
+ ,524
+ ,552
+ ,552
+ ,532
+ ,532
+ ,511
+ ,511
+ ,492
+ ,492
+ ,492
+ ,492
+ ,493
+ ,493
+ ,481
+ ,481
+ ,462
+ ,462
+ ,457
+ ,457
+ ,442
+ ,442
+ ,439
+ ,439
+ ,488
+ ,488
+ ,521
+ ,521
+ ,501
+ ,501
+ ,485
+ ,485
+ ,464
+ ,464
+ ,460
+ ,460
+ ,467
+ ,467
+ ,460
+ ,460
+ ,448
+ ,448
+ ,443
+ ,443
+ ,436
+ ,436
+ ,431
+ ,431
+ ,484
+ ,484
+ ,510
+ ,510
+ ,513
+ ,513
+ ,503
+ ,503
+ ,471
+ ,471
+ ,471
+ ,471
+ ,476
+ ,476
+ ,475
+ ,475
+ ,470
+ ,470
+ ,461
+ ,461
+ ,455
+ ,455
+ ,456
+ ,456
+ ,517
+ ,517
+ ,525
+ ,525
+ ,523
+ ,523
+ ,519
+ ,519
+ ,509
+ ,509
+ ,512
+ ,512
+ ,519
+ ,519
+ ,517
+ ,517
+ ,510
+ ,510
+ ,509
+ ,509
+ ,501
+ ,501
+ ,507
+ ,507
+ ,569
+ ,569
+ ,580
+ ,580
+ ,578
+ ,578
+ ,565
+ ,565
+ ,547
+ ,547
+ ,555
+ ,555
+ ,562
+ ,561
+ ,555
+ ,544
+ ,537
+ ,543
+ ,594
+ ,611
+ ,613
+ ,611
+ ,594
+ ,595
+ ,591
+ ,589
+ ,584
+ ,573
+ ,567
+ ,569
+ ,621
+ ,629
+ ,628
+ ,612
+ ,595
+ ,597
+ ,593
+ ,590
+ ,580
+ ,574
+ ,573
+ ,573
+ ,620
+ ,626
+ ,620
+ ,588
+ ,566
+ ,557
+ ,561
+ ,549
+ ,532
+ ,526
+ ,511
+ ,499
+ ,555
+ ,565
+ ,542
+ ,527
+ ,510
+ ,514
+ ,517
+ ,508
+ ,493
+ ,490
+ ,469
+ ,478
+ ,528
+ ,534
+ ,518
+ ,506
+ ,502
+ ,516
+ ,528
+ ,533
+ ,536
+ ,537
+ ,524
+ ,536
+ ,587
+ ,597
+ ,581
+ ,564
+ ,558
+ ,575
+ ,580
+ ,575
+ ,563
+ ,552
+ ,537
+ ,545
+ ,601
+ ,604
+ ,586
+ ,564
+ ,549
+ ,551)
+ ,dim=c(1
+ ,492)
+ ,dimnames=list(c('Werkloosheid')
+ ,1:492))
> y <- array(NA,dim=c(1,492),dimnames=list(c('Werkloosheid'),1:492))
> 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
Werkloosheid t
1 349 1
2 336 2
3 331 3
4 327 4
5 323 5
6 322 6
7 385 7
8 405 8
9 412 9
10 411 10
11 410 11
12 415 12
13 414 13
14 411 14
15 408 15
16 410 16
17 411 17
18 416 18
19 479 19
20 498 20
21 502 21
22 498 22
23 499 23
24 506 24
25 510 25
26 509 26
27 502 27
28 495 28
29 490 29
30 490 30
31 553 31
32 570 32
33 573 33
34 572 34
35 575 35
36 580 36
37 580 37
38 574 38
39 563 39
40 556 40
41 546 41
42 545 42
43 605 43
44 628 44
45 631 45
46 626 46
47 614 47
48 606 48
49 602 49
50 589 50
51 574 51
52 558 52
53 552 53
54 546 54
55 607 55
56 636 56
57 631 57
58 623 58
59 618 59
60 605 60
61 619 61
62 596 62
63 570 63
64 546 64
65 528 65
66 506 66
67 555 67
68 568 68
69 564 69
70 553 70
71 541 71
72 542 72
73 540 73
74 521 74
75 505 75
76 491 76
77 482 77
78 478 78
79 523 79
80 531 80
81 532 81
82 540 82
83 525 83
84 533 84
85 531 85
86 508 86
87 495 87
88 482 88
89 470 89
90 466 90
91 515 91
92 518 92
93 516 93
94 511 94
95 500 95
96 498 96
97 494 97
98 476 98
99 458 99
100 443 100
101 430 101
102 424 102
103 476 103
104 481 104
105 470 105
106 460 106
107 451 107
108 450 108
109 444 109
110 429 110
111 421 111
112 400 112
113 389 113
114 384 114
115 432 115
116 446 116
117 431 117
118 423 118
119 416 119
120 416 120
121 413 121
122 399 122
123 386 123
124 374 124
125 365 125
126 365 126
127 418 127
128 428 128
129 424 129
130 421 130
131 417 131
132 423 132
133 423 133
134 419 134
135 406 135
136 398 136
137 390 137
138 391 138
139 444 139
140 460 140
141 455 141
142 456 142
143 452 143
144 459 144
145 461 145
146 451 146
147 443 147
148 439 148
149 430 149
150 436 150
151 488 151
152 506 152
153 502 153
154 501 154
155 501 155
156 515 156
157 521 157
158 520 158
159 512 159
160 509 160
161 505 161
162 511 162
163 570 163
164 592 164
165 594 165
166 586 166
167 586 167
168 592 168
169 594 169
170 594 170
171 586 171
172 586 172
173 572 173
174 572 174
175 563 175
176 563 176
177 555 177
178 555 178
179 554 179
180 554 180
181 601 181
182 601 182
183 622 183
184 622 184
185 617 185
186 617 186
187 606 187
188 606 188
189 595 189
190 595 190
191 599 191
192 599 192
193 600 193
194 600 194
195 592 195
196 592 196
197 575 197
198 575 198
199 567 199
200 567 200
201 555 201
202 555 202
203 555 203
204 555 204
205 608 205
206 608 206
207 631 207
208 631 208
209 629 209
210 629 210
211 624 211
212 624 212
213 610 213
214 610 214
215 616 215
216 616 216
217 621 217
218 621 218
219 604 219
220 604 220
221 584 221
222 584 222
223 574 223
224 574 224
225 555 225
226 555 226
227 545 227
228 545 228
229 599 229
230 599 230
231 620 231
232 620 232
233 608 233
234 608 234
235 590 235
236 590 236
237 579 237
238 579 238
239 580 239
240 580 240
241 579 241
242 579 242
243 572 243
244 572 244
245 560 245
246 560 246
247 551 247
248 551 248
249 537 249
250 537 250
251 541 251
252 541 252
253 588 253
254 588 254
255 607 255
256 607 256
257 599 257
258 599 258
259 578 259
260 578 260
261 563 261
262 563 262
263 566 263
264 566 264
265 561 265
266 561 266
267 554 267
268 554 268
269 540 269
270 540 270
271 526 271
272 526 272
273 512 273
274 512 274
275 505 275
276 505 276
277 554 277
278 554 278
279 584 279
280 584 280
281 569 281
282 569 282
283 540 283
284 540 284
285 522 285
286 522 286
287 526 287
288 526 288
289 527 289
290 527 290
291 516 291
292 516 292
293 503 293
294 503 294
295 489 295
296 489 296
297 479 297
298 479 298
299 475 299
300 475 300
301 524 301
302 524 302
303 552 303
304 552 304
305 532 305
306 532 306
307 511 307
308 511 308
309 492 309
310 492 310
311 492 311
312 492 312
313 493 313
314 493 314
315 481 315
316 481 316
317 462 317
318 462 318
319 457 319
320 457 320
321 442 321
322 442 322
323 439 323
324 439 324
325 488 325
326 488 326
327 521 327
328 521 328
329 501 329
330 501 330
331 485 331
332 485 332
333 464 333
334 464 334
335 460 335
336 460 336
337 467 337
338 467 338
339 460 339
340 460 340
341 448 341
342 448 342
343 443 343
344 443 344
345 436 345
346 436 346
347 431 347
348 431 348
349 484 349
350 484 350
351 510 351
352 510 352
353 513 353
354 513 354
355 503 355
356 503 356
357 471 357
358 471 358
359 471 359
360 471 360
361 476 361
362 476 362
363 475 363
364 475 364
365 470 365
366 470 366
367 461 367
368 461 368
369 455 369
370 455 370
371 456 371
372 456 372
373 517 373
374 517 374
375 525 375
376 525 376
377 523 377
378 523 378
379 519 379
380 519 380
381 509 381
382 509 382
383 512 383
384 512 384
385 519 385
386 519 386
387 517 387
388 517 388
389 510 389
390 510 390
391 509 391
392 509 392
393 501 393
394 501 394
395 507 395
396 507 396
397 569 397
398 569 398
399 580 399
400 580 400
401 578 401
402 578 402
403 565 403
404 565 404
405 547 405
406 547 406
407 555 407
408 555 408
409 562 409
410 561 410
411 555 411
412 544 412
413 537 413
414 543 414
415 594 415
416 611 416
417 613 417
418 611 418
419 594 419
420 595 420
421 591 421
422 589 422
423 584 423
424 573 424
425 567 425
426 569 426
427 621 427
428 629 428
429 628 429
430 612 430
431 595 431
432 597 432
433 593 433
434 590 434
435 580 435
436 574 436
437 573 437
438 573 438
439 620 439
440 626 440
441 620 441
442 588 442
443 566 443
444 557 444
445 561 445
446 549 446
447 532 447
448 526 448
449 511 449
450 499 450
451 555 451
452 565 452
453 542 453
454 527 454
455 510 455
456 514 456
457 517 457
458 508 458
459 493 459
460 490 460
461 469 461
462 478 462
463 528 463
464 534 464
465 518 465
466 506 466
467 502 467
468 516 468
469 528 469
470 533 470
471 536 471
472 537 472
473 524 473
474 536 474
475 587 475
476 597 476
477 581 477
478 564 478
479 558 479
480 575 480
481 580 481
482 575 482
483 563 483
484 552 484
485 537 485
486 545 486
487 601 487
488 604 488
489 586 489
490 564 490
491 549 491
492 551 492
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) t
499.3033 0.1022
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-177.9165 -44.7236 0.2510 49.1475 130.9734
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 499.30327 5.70448 87.528 < 2e-16 ***
t 0.10220 0.02005 5.097 4.94e-07 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 63.17 on 490 degrees of freedom
Multiple R-squared: 0.05035, Adjusted R-squared: 0.04841
F-statistic: 25.98 on 1 and 490 DF, p-value: 4.939e-07
> 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,] 2.177875e-04 4.355750e-04 9.997822e-01
[2,] 2.981224e-05 5.962449e-05 9.999702e-01
[3,] 2.096144e-02 4.192288e-02 9.790386e-01
[4,] 2.680237e-02 5.360475e-02 9.731976e-01
[5,] 1.624438e-02 3.248877e-02 9.837556e-01
[6,] 6.987586e-03 1.397517e-02 9.930124e-01
[7,] 2.800421e-03 5.600842e-03 9.971996e-01
[8,] 1.101694e-03 2.203387e-03 9.988983e-01
[9,] 4.759062e-04 9.518125e-04 9.995241e-01
[10,] 2.499589e-04 4.999177e-04 9.997500e-01
[11,] 1.598831e-04 3.197663e-04 9.998401e-01
[12,] 9.847317e-05 1.969463e-04 9.999015e-01
[13,] 6.085964e-05 1.217193e-04 9.999391e-01
[14,] 3.278856e-05 6.557713e-05 9.999672e-01
[15,] 3.396103e-05 6.792206e-05 9.999660e-01
[16,] 4.238886e-05 8.477772e-05 9.999576e-01
[17,] 3.116620e-05 6.233240e-05 9.999688e-01
[18,] 1.436333e-05 2.872666e-05 9.999856e-01
[19,] 5.808088e-06 1.161618e-05 9.999942e-01
[20,] 2.271187e-06 4.542373e-06 9.999977e-01
[21,] 8.599782e-07 1.719956e-06 9.999991e-01
[22,] 3.401638e-07 6.803276e-07 9.999997e-01
[23,] 1.842236e-07 3.684472e-07 9.999998e-01
[24,] 1.613642e-07 3.227284e-07 9.999998e-01
[25,] 2.079138e-07 4.158276e-07 9.999998e-01
[26,] 2.764919e-07 5.529838e-07 9.999997e-01
[27,] 1.357055e-07 2.714110e-07 9.999999e-01
[28,] 8.345517e-08 1.669103e-07 9.999999e-01
[29,] 4.209355e-08 8.418710e-08 1.000000e+00
[30,] 1.761001e-08 3.522002e-08 1.000000e+00
[31,] 7.014263e-09 1.402853e-08 1.000000e+00
[32,] 2.753155e-09 5.506310e-09 1.000000e+00
[33,] 1.103320e-09 2.206640e-09 1.000000e+00
[34,] 5.565713e-10 1.113143e-09 1.000000e+00
[35,] 5.368860e-10 1.073772e-09 1.000000e+00
[36,] 9.131361e-10 1.826272e-09 1.000000e+00
[37,] 3.163662e-09 6.327323e-09 1.000000e+00
[38,] 1.032917e-08 2.065834e-08 1.000000e+00
[39,] 4.797893e-09 9.595786e-09 1.000000e+00
[40,] 2.648705e-09 5.297409e-09 1.000000e+00
[41,] 1.379157e-09 2.758314e-09 1.000000e+00
[42,] 6.677324e-10 1.335465e-09 1.000000e+00
[43,] 3.986230e-10 7.972460e-10 1.000000e+00
[44,] 3.503651e-10 7.007302e-10 1.000000e+00
[45,] 4.218567e-10 8.437134e-10 1.000000e+00
[46,] 1.112891e-09 2.225782e-09 1.000000e+00
[47,] 7.004771e-09 1.400954e-08 1.000000e+00
[48,] 8.871040e-08 1.774208e-07 9.999999e-01
[49,] 8.999910e-07 1.799982e-06 9.999991e-01
[50,] 7.240340e-06 1.448068e-05 9.999928e-01
[51,] 6.476292e-06 1.295258e-05 9.999935e-01
[52,] 4.647870e-06 9.295741e-06 9.999954e-01
[53,] 3.554816e-06 7.109631e-06 9.999964e-01
[54,] 3.125463e-06 6.250926e-06 9.999969e-01
[55,] 3.156945e-06 6.313890e-06 9.999968e-01
[56,] 4.431174e-06 8.862347e-06 9.999956e-01
[57,] 4.796542e-06 9.593084e-06 9.999952e-01
[58,] 9.228629e-06 1.845726e-05 9.999908e-01
[59,] 3.790947e-05 7.581894e-05 9.999621e-01
[60,] 2.674753e-04 5.349506e-04 9.997325e-01
[61,] 1.987366e-03 3.974732e-03 9.980126e-01
[62,] 1.412393e-02 2.824786e-02 9.858761e-01
[63,] 2.365502e-02 4.731004e-02 9.763450e-01
[64,] 3.159802e-02 6.319604e-02 9.684020e-01
[65,] 4.248512e-02 8.497025e-02 9.575149e-01
[66,] 6.111873e-02 1.222375e-01 9.388813e-01
[67,] 9.283176e-02 1.856635e-01 9.071682e-01
[68,] 1.291102e-01 2.582204e-01 8.708898e-01
[69,] 1.715469e-01 3.430938e-01 8.284531e-01
[70,] 2.457280e-01 4.914559e-01 7.542720e-01
[71,] 3.556170e-01 7.112340e-01 6.443830e-01
[72,] 4.950980e-01 9.901959e-01 5.049020e-01
[73,] 6.366715e-01 7.266571e-01 3.633285e-01
[74,] 7.531615e-01 4.936771e-01 2.468385e-01
[75,] 7.811847e-01 4.376306e-01 2.188153e-01
[76,] 7.973619e-01 4.052762e-01 2.026381e-01
[77,] 8.101389e-01 3.797221e-01 1.898611e-01
[78,] 8.155484e-01 3.689033e-01 1.844516e-01
[79,] 8.298109e-01 3.403782e-01 1.701891e-01
[80,] 8.365013e-01 3.269974e-01 1.634987e-01
[81,] 8.430636e-01 3.138728e-01 1.569364e-01
[82,] 8.630974e-01 2.738052e-01 1.369026e-01
[83,] 8.877983e-01 2.244033e-01 1.122017e-01
[84,] 9.145740e-01 1.708521e-01 8.542603e-02
[85,] 9.398780e-01 1.202441e-01 6.012204e-02
[86,] 9.583858e-01 8.322848e-02 4.161424e-02
[87,] 9.589658e-01 8.206843e-02 4.103422e-02
[88,] 9.586455e-01 8.270903e-02 4.135451e-02
[89,] 9.583703e-01 8.325932e-02 4.162966e-02
[90,] 9.586797e-01 8.264053e-02 4.132026e-02
[91,] 9.607659e-01 7.846827e-02 3.923413e-02
[92,] 9.626816e-01 7.463683e-02 3.731841e-02
[93,] 9.648266e-01 7.034675e-02 3.517337e-02
[94,] 9.699425e-01 6.011494e-02 3.005747e-02
[95,] 9.772228e-01 4.555449e-02 2.277725e-02
[96,] 9.846447e-01 3.071054e-02 1.535527e-02
[97,] 9.907893e-01 1.842146e-02 9.210732e-03
[98,] 9.947371e-01 1.052573e-02 5.262863e-03
[99,] 9.949332e-01 1.013363e-02 5.066813e-03
[100,] 9.948929e-01 1.021420e-02 5.107099e-03
[101,] 9.951513e-01 9.697378e-03 4.848689e-03
[102,] 9.956603e-01 8.679467e-03 4.339734e-03
[103,] 9.963290e-01 7.342008e-03 3.671004e-03
[104,] 9.968608e-01 6.278475e-03 3.139238e-03
[105,] 9.974067e-01 5.186615e-03 2.593308e-03
[106,] 9.981185e-01 3.763093e-03 1.881546e-03
[107,] 9.987256e-01 2.548799e-03 1.274400e-03
[108,] 9.993298e-01 1.340454e-03 6.702271e-04
[109,] 9.996976e-01 6.048779e-04 3.024389e-04
[110,] 9.998734e-01 2.532525e-04 1.266263e-04
[111,] 9.998946e-01 2.108599e-04 1.054299e-04
[112,] 9.998977e-01 2.046125e-04 1.023063e-04
[113,] 9.999135e-01 1.730314e-04 8.651568e-05
[114,] 9.999323e-01 1.353690e-04 6.768452e-05
[115,] 9.999508e-01 9.837741e-05 4.918870e-05
[116,] 9.999639e-01 7.227805e-05 3.613902e-05
[117,] 9.999742e-01 5.154892e-05 2.577446e-05
[118,] 9.999848e-01 3.031156e-05 1.515578e-05
[119,] 9.999928e-01 1.439378e-05 7.196890e-06
[120,] 9.999973e-01 5.392777e-06 2.696388e-06
[121,] 9.999992e-01 1.623499e-06 8.117493e-07
[122,] 9.999998e-01 4.673970e-07 2.336985e-07
[123,] 9.999998e-01 3.535290e-07 1.767645e-07
[124,] 9.999998e-01 3.043645e-07 1.521822e-07
[125,] 9.999999e-01 2.472807e-07 1.236403e-07
[126,] 9.999999e-01 1.912641e-07 9.563203e-08
[127,] 9.999999e-01 1.381375e-07 6.906877e-08
[128,] 9.999999e-01 1.074398e-07 5.371990e-08
[129,] 1.000000e+00 8.247646e-08 4.123823e-08
[130,] 1.000000e+00 5.879610e-08 2.939805e-08
[131,] 1.000000e+00 3.308685e-08 1.654343e-08
[132,] 1.000000e+00 1.545121e-08 7.725603e-09
[133,] 1.000000e+00 5.791377e-09 2.895688e-09
[134,] 1.000000e+00 2.081298e-09 1.040649e-09
[135,] 1.000000e+00 1.822372e-09 9.111860e-10
[136,] 1.000000e+00 1.853445e-09 9.267226e-10
[137,] 1.000000e+00 1.773724e-09 8.868618e-10
[138,] 1.000000e+00 1.687575e-09 8.437873e-10
[139,] 1.000000e+00 1.513900e-09 7.569502e-10
[140,] 1.000000e+00 1.434489e-09 7.172445e-10
[141,] 1.000000e+00 1.363307e-09 6.816536e-10
[142,] 1.000000e+00 1.143835e-09 5.719174e-10
[143,] 1.000000e+00 8.438224e-10 4.219112e-10
[144,] 1.000000e+00 5.686065e-10 2.843033e-10
[145,] 1.000000e+00 3.164614e-10 1.582307e-10
[146,] 1.000000e+00 1.867738e-10 9.338691e-11
[147,] 1.000000e+00 1.897514e-10 9.487570e-11
[148,] 1.000000e+00 2.009613e-10 1.004806e-10
[149,] 1.000000e+00 2.121983e-10 1.060991e-10
[150,] 1.000000e+00 2.234407e-10 1.117203e-10
[151,] 1.000000e+00 2.348940e-10 1.174470e-10
[152,] 1.000000e+00 2.496069e-10 1.248034e-10
[153,] 1.000000e+00 2.644764e-10 1.322382e-10
[154,] 1.000000e+00 2.829150e-10 1.414575e-10
[155,] 1.000000e+00 3.054787e-10 1.527394e-10
[156,] 1.000000e+00 3.292407e-10 1.646204e-10
[157,] 1.000000e+00 3.515498e-10 1.757749e-10
[158,] 1.000000e+00 3.797320e-10 1.898660e-10
[159,] 1.000000e+00 2.971307e-10 1.485654e-10
[160,] 1.000000e+00 1.726037e-10 8.630183e-11
[161,] 1.000000e+00 1.015258e-10 5.076291e-11
[162,] 1.000000e+00 7.159185e-11 3.579592e-11
[163,] 1.000000e+00 5.251985e-11 2.625993e-11
[164,] 1.000000e+00 3.634439e-11 1.817220e-11
[165,] 1.000000e+00 2.528821e-11 1.264410e-11
[166,] 1.000000e+00 1.827300e-11 9.136498e-12
[167,] 1.000000e+00 1.547981e-11 7.739906e-12
[168,] 1.000000e+00 1.350404e-11 6.752022e-12
[169,] 1.000000e+00 1.435136e-11 7.175678e-12
[170,] 1.000000e+00 1.554150e-11 7.770750e-12
[171,] 1.000000e+00 1.853687e-11 9.268435e-12
[172,] 1.000000e+00 2.238533e-11 1.119267e-11
[173,] 1.000000e+00 2.874964e-11 1.437482e-11
[174,] 1.000000e+00 3.720398e-11 1.860199e-11
[175,] 1.000000e+00 4.871496e-11 2.435748e-11
[176,] 1.000000e+00 6.417684e-11 3.208842e-11
[177,] 1.000000e+00 5.326695e-11 2.663347e-11
[178,] 1.000000e+00 4.524425e-11 2.262212e-11
[179,] 1.000000e+00 2.697984e-11 1.348992e-11
[180,] 1.000000e+00 1.648261e-11 8.241306e-12
[181,] 1.000000e+00 1.136749e-11 5.683747e-12
[182,] 1.000000e+00 7.997880e-12 3.998940e-12
[183,] 1.000000e+00 6.962621e-12 3.481310e-12
[184,] 1.000000e+00 6.153970e-12 3.076985e-12
[185,] 1.000000e+00 6.486164e-12 3.243082e-12
[186,] 1.000000e+00 6.907915e-12 3.453958e-12
[187,] 1.000000e+00 7.033369e-12 3.516684e-12
[188,] 1.000000e+00 7.228701e-12 3.614350e-12
[189,] 1.000000e+00 7.387873e-12 3.693937e-12
[190,] 1.000000e+00 7.608867e-12 3.804434e-12
[191,] 1.000000e+00 8.757877e-12 4.378938e-12
[192,] 1.000000e+00 1.012915e-11 5.064574e-12
[193,] 1.000000e+00 1.389323e-11 6.946617e-12
[194,] 1.000000e+00 1.907087e-11 9.535433e-12
[195,] 1.000000e+00 2.759881e-11 1.379941e-11
[196,] 1.000000e+00 3.989711e-11 1.994856e-11
[197,] 1.000000e+00 6.043812e-11 3.021906e-11
[198,] 1.000000e+00 9.128856e-11 4.564428e-11
[199,] 1.000000e+00 1.374743e-10 6.873717e-11
[200,] 1.000000e+00 2.063920e-10 1.031960e-10
[201,] 1.000000e+00 1.976011e-10 9.880054e-11
[202,] 1.000000e+00 1.894275e-10 9.471375e-11
[203,] 1.000000e+00 1.218433e-10 6.092166e-11
[204,] 1.000000e+00 7.776122e-11 3.888061e-11
[205,] 1.000000e+00 5.119402e-11 2.559701e-11
[206,] 1.000000e+00 3.334325e-11 1.667163e-11
[207,] 1.000000e+00 2.372297e-11 1.186148e-11
[208,] 1.000000e+00 1.668256e-11 8.341281e-12
[209,] 1.000000e+00 1.496797e-11 7.483983e-12
[210,] 1.000000e+00 1.331937e-11 6.659686e-12
[211,] 1.000000e+00 1.056101e-11 5.280505e-12
[212,] 1.000000e+00 8.262513e-12 4.131256e-12
[213,] 1.000000e+00 5.771270e-12 2.885635e-12
[214,] 1.000000e+00 3.950103e-12 1.975052e-12
[215,] 1.000000e+00 3.640283e-12 1.820142e-12
[216,] 1.000000e+00 3.311270e-12 1.655635e-12
[217,] 1.000000e+00 3.942842e-12 1.971421e-12
[218,] 1.000000e+00 4.656605e-12 2.328302e-12
[219,] 1.000000e+00 6.041520e-12 3.020760e-12
[220,] 1.000000e+00 7.783674e-12 3.891837e-12
[221,] 1.000000e+00 1.125958e-11 5.629791e-12
[222,] 1.000000e+00 1.620618e-11 8.103090e-12
[223,] 1.000000e+00 2.384459e-11 1.192229e-11
[224,] 1.000000e+00 3.494349e-11 1.747175e-11
[225,] 1.000000e+00 3.234871e-11 1.617436e-11
[226,] 1.000000e+00 2.948146e-11 1.474073e-11
[227,] 1.000000e+00 1.776343e-11 8.881715e-12
[228,] 1.000000e+00 1.030180e-11 5.150901e-12
[229,] 1.000000e+00 7.411705e-12 3.705852e-12
[230,] 1.000000e+00 5.175146e-12 2.587573e-12
[231,] 1.000000e+00 4.847556e-12 2.423778e-12
[232,] 1.000000e+00 4.458077e-12 2.229039e-12
[233,] 1.000000e+00 4.715393e-12 2.357696e-12
[234,] 1.000000e+00 4.918918e-12 2.459459e-12
[235,] 1.000000e+00 4.990138e-12 2.495069e-12
[236,] 1.000000e+00 4.984044e-12 2.492022e-12
[237,] 1.000000e+00 4.966304e-12 2.483152e-12
[238,] 1.000000e+00 4.867896e-12 2.433948e-12
[239,] 1.000000e+00 5.148104e-12 2.574052e-12
[240,] 1.000000e+00 5.369665e-12 2.684832e-12
[241,] 1.000000e+00 6.294409e-12 3.147205e-12
[242,] 1.000000e+00 7.310748e-12 3.655374e-12
[243,] 1.000000e+00 9.048317e-12 4.524158e-12
[244,] 1.000000e+00 1.112732e-11 5.563662e-12
[245,] 1.000000e+00 1.457316e-11 7.286580e-12
[246,] 1.000000e+00 1.903623e-11 9.518113e-12
[247,] 1.000000e+00 2.443863e-11 1.221932e-11
[248,] 1.000000e+00 3.125611e-11 1.562805e-11
[249,] 1.000000e+00 2.341502e-11 1.170751e-11
[250,] 1.000000e+00 1.698239e-11 8.491195e-12
[251,] 1.000000e+00 7.756413e-12 3.878207e-12
[252,] 1.000000e+00 3.287654e-12 1.643827e-12
[253,] 1.000000e+00 1.589456e-12 7.947278e-13
[254,] 1.000000e+00 7.166239e-13 3.583120e-13
[255,] 1.000000e+00 4.961779e-13 2.480889e-13
[256,] 1.000000e+00 3.304717e-13 1.652359e-13
[257,] 1.000000e+00 2.810807e-13 1.405404e-13
[258,] 1.000000e+00 2.334466e-13 1.167233e-13
[259,] 1.000000e+00 1.791666e-13 8.958332e-14
[260,] 1.000000e+00 1.333733e-13 6.668665e-14
[261,] 1.000000e+00 1.055687e-13 5.278434e-14
[262,] 1.000000e+00 8.127281e-14 4.063641e-14
[263,] 1.000000e+00 6.874826e-14 3.437413e-14
[264,] 1.000000e+00 5.686515e-14 2.843258e-14
[265,] 1.000000e+00 5.636223e-14 2.818112e-14
[266,] 1.000000e+00 5.525355e-14 2.762678e-14
[267,] 1.000000e+00 6.145106e-14 3.072553e-14
[268,] 1.000000e+00 6.817226e-14 3.408613e-14
[269,] 1.000000e+00 8.072640e-14 4.036320e-14
[270,] 1.000000e+00 9.594485e-14 4.797242e-14
[271,] 1.000000e+00 1.153274e-13 5.766369e-14
[272,] 1.000000e+00 1.394462e-13 6.972310e-14
[273,] 1.000000e+00 1.076883e-13 5.384417e-14
[274,] 1.000000e+00 8.079812e-14 4.039906e-14
[275,] 1.000000e+00 2.677651e-14 1.338825e-14
[276,] 1.000000e+00 7.911671e-15 3.955836e-15
[277,] 1.000000e+00 3.387341e-15 1.693671e-15
[278,] 1.000000e+00 1.336132e-15 6.680661e-16
[279,] 1.000000e+00 1.028801e-15 5.144003e-16
[280,] 1.000000e+00 7.725255e-16 3.862627e-16
[281,] 1.000000e+00 7.599782e-16 3.799891e-16
[282,] 1.000000e+00 7.429346e-16 3.714673e-16
[283,] 1.000000e+00 6.823526e-16 3.411763e-16
[284,] 1.000000e+00 6.197791e-16 3.098896e-16
[285,] 1.000000e+00 5.475541e-16 2.737770e-16
[286,] 1.000000e+00 4.772013e-16 2.386006e-16
[287,] 1.000000e+00 4.808160e-16 2.404080e-16
[288,] 1.000000e+00 4.825635e-16 2.412818e-16
[289,] 1.000000e+00 5.447529e-16 2.723764e-16
[290,] 1.000000e+00 6.176464e-16 3.088232e-16
[291,] 1.000000e+00 7.364716e-16 3.682358e-16
[292,] 1.000000e+00 8.870290e-16 4.435145e-16
[293,] 1.000000e+00 1.055468e-15 5.277342e-16
[294,] 1.000000e+00 1.271377e-15 6.356885e-16
[295,] 1.000000e+00 1.516409e-15 7.582046e-16
[296,] 1.000000e+00 1.830623e-15 9.153115e-16
[297,] 1.000000e+00 1.768115e-15 8.840577e-16
[298,] 1.000000e+00 1.687592e-15 8.437962e-16
[299,] 1.000000e+00 8.203039e-16 4.101520e-16
[300,] 1.000000e+00 3.717812e-16 1.858906e-16
[301,] 1.000000e+00 2.753761e-16 1.376880e-16
[302,] 1.000000e+00 1.980494e-16 9.902469e-17
[303,] 1.000000e+00 2.072706e-16 1.036353e-16
[304,] 1.000000e+00 2.159840e-16 1.079920e-16
[305,] 1.000000e+00 2.693574e-16 1.346787e-16
[306,] 1.000000e+00 3.381629e-16 1.690815e-16
[307,] 1.000000e+00 4.272927e-16 2.136464e-16
[308,] 1.000000e+00 5.433038e-16 2.716519e-16
[309,] 1.000000e+00 6.913295e-16 3.456648e-16
[310,] 1.000000e+00 8.845556e-16 4.422778e-16
[311,] 1.000000e+00 1.175408e-15 5.877038e-16
[312,] 1.000000e+00 1.575115e-15 7.875575e-16
[313,] 1.000000e+00 1.937986e-15 9.689932e-16
[314,] 1.000000e+00 2.406406e-15 1.203203e-15
[315,] 1.000000e+00 2.851732e-15 1.425866e-15
[316,] 1.000000e+00 3.403688e-15 1.701844e-15
[317,] 1.000000e+00 3.206055e-15 1.603028e-15
[318,] 1.000000e+00 3.020939e-15 1.510470e-15
[319,] 1.000000e+00 2.666265e-15 1.333132e-15
[320,] 1.000000e+00 2.340726e-15 1.170363e-15
[321,] 1.000000e+00 3.555727e-15 1.777863e-15
[322,] 1.000000e+00 5.418677e-15 2.709338e-15
[323,] 1.000000e+00 6.401479e-15 3.200740e-15
[324,] 1.000000e+00 7.500214e-15 3.750107e-15
[325,] 1.000000e+00 1.093255e-14 5.466274e-15
[326,] 1.000000e+00 1.595082e-14 7.975411e-15
[327,] 1.000000e+00 2.439060e-14 1.219530e-14
[328,] 1.000000e+00 3.739656e-14 1.869828e-14
[329,] 1.000000e+00 5.085319e-14 2.542660e-14
[330,] 1.000000e+00 6.914508e-14 3.457254e-14
[331,] 1.000000e+00 8.953191e-14 4.476595e-14
[332,] 1.000000e+00 1.156002e-13 5.780012e-14
[333,] 1.000000e+00 1.616957e-13 8.084784e-14
[334,] 1.000000e+00 2.256429e-13 1.128214e-13
[335,] 1.000000e+00 2.876595e-13 1.438297e-13
[336,] 1.000000e+00 3.643715e-13 1.821857e-13
[337,] 1.000000e+00 3.718563e-13 1.859281e-13
[338,] 1.000000e+00 3.730190e-13 1.865095e-13
[339,] 1.000000e+00 3.281189e-13 1.640595e-13
[340,] 1.000000e+00 2.806169e-13 1.403085e-13
[341,] 1.000000e+00 1.927065e-13 9.635325e-14
[342,] 1.000000e+00 1.259918e-13 6.299591e-14
[343,] 1.000000e+00 6.630954e-14 3.315477e-14
[344,] 1.000000e+00 3.227109e-14 1.613555e-14
[345,] 1.000000e+00 5.062201e-14 2.531100e-14
[346,] 1.000000e+00 7.868019e-14 3.934010e-14
[347,] 1.000000e+00 1.393196e-13 6.965978e-14
[348,] 1.000000e+00 2.458186e-13 1.229093e-13
[349,] 1.000000e+00 4.307812e-13 2.153906e-13
[350,] 1.000000e+00 7.522785e-13 3.761393e-13
[351,] 1.000000e+00 1.299291e-12 6.496456e-13
[352,] 1.000000e+00 2.233030e-12 1.116515e-12
[353,] 1.000000e+00 2.740370e-12 1.370185e-12
[354,] 1.000000e+00 3.291482e-12 1.645741e-12
[355,] 1.000000e+00 3.862388e-12 1.931194e-12
[356,] 1.000000e+00 4.419133e-12 2.209567e-12
[357,] 1.000000e+00 5.405490e-12 2.702745e-12
[358,] 1.000000e+00 6.455436e-12 3.227718e-12
[359,] 1.000000e+00 7.368956e-12 3.684478e-12
[360,] 1.000000e+00 8.168255e-12 4.084128e-12
[361,] 1.000000e+00 7.850809e-12 3.925405e-12
[362,] 1.000000e+00 7.224176e-12 3.612088e-12
[363,] 1.000000e+00 4.940973e-12 2.470486e-12
[364,] 1.000000e+00 3.123036e-12 1.561518e-12
[365,] 1.000000e+00 1.461943e-12 7.309713e-13
[366,] 1.000000e+00 5.997138e-13 2.998569e-13
[367,] 1.000000e+00 2.197936e-13 1.098968e-13
[368,] 1.000000e+00 6.753142e-14 3.376571e-14
[369,] 1.000000e+00 1.097611e-13 5.488055e-14
[370,] 1.000000e+00 1.757143e-13 8.785717e-14
[371,] 1.000000e+00 3.037236e-13 1.518618e-13
[372,] 1.000000e+00 5.192255e-13 2.596128e-13
[373,] 1.000000e+00 8.592840e-13 4.296420e-13
[374,] 1.000000e+00 1.403633e-12 7.018165e-13
[375,] 1.000000e+00 2.148558e-12 1.074279e-12
[376,] 1.000000e+00 3.232388e-12 1.616194e-12
[377,] 1.000000e+00 4.006558e-12 2.003279e-12
[378,] 1.000000e+00 4.809561e-12 2.404781e-12
[379,] 1.000000e+00 5.941918e-12 2.970959e-12
[380,] 1.000000e+00 7.107951e-12 3.553976e-12
[381,] 1.000000e+00 9.452716e-12 4.726358e-12
[382,] 1.000000e+00 1.224578e-11 6.122891e-12
[383,] 1.000000e+00 1.481487e-11 7.407436e-12
[384,] 1.000000e+00 1.733055e-11 8.665273e-12
[385,] 1.000000e+00 1.646378e-11 8.231889e-12
[386,] 1.000000e+00 1.476760e-11 7.383798e-12
[387,] 1.000000e+00 1.205648e-11 6.028239e-12
[388,] 1.000000e+00 9.120912e-12 4.560456e-12
[389,] 1.000000e+00 4.715843e-12 2.357921e-12
[390,] 1.000000e+00 2.104592e-12 1.052296e-12
[391,] 1.000000e+00 1.029359e-12 5.146795e-13
[392,] 1.000000e+00 4.292784e-13 2.146392e-13
[393,] 1.000000e+00 7.951440e-13 3.975720e-13
[394,] 1.000000e+00 1.467222e-12 7.336112e-13
[395,] 1.000000e+00 2.743494e-12 1.371747e-12
[396,] 1.000000e+00 5.110990e-12 2.555495e-12
[397,] 1.000000e+00 9.518783e-12 4.759391e-12
[398,] 1.000000e+00 1.765894e-11 8.829472e-12
[399,] 1.000000e+00 3.159063e-11 1.579531e-11
[400,] 1.000000e+00 5.621350e-11 2.810675e-11
[401,] 1.000000e+00 7.998973e-11 3.999486e-11
[402,] 1.000000e+00 1.115655e-10 5.578277e-11
[403,] 1.000000e+00 1.736859e-10 8.684296e-11
[404,] 1.000000e+00 2.670296e-10 1.335148e-10
[405,] 1.000000e+00 4.422407e-10 2.211203e-10
[406,] 1.000000e+00 7.188240e-10 3.594120e-10
[407,] 1.000000e+00 1.070730e-09 5.353651e-10
[408,] 1.000000e+00 1.281563e-09 6.407813e-10
[409,] 1.000000e+00 1.243873e-09 6.219365e-10
[410,] 1.000000e+00 1.338455e-09 6.692276e-10
[411,] 1.000000e+00 2.359903e-09 1.179952e-09
[412,] 1.000000e+00 3.472234e-09 1.736117e-09
[413,] 1.000000e+00 4.821081e-09 2.410541e-09
[414,] 1.000000e+00 6.726093e-09 3.363046e-09
[415,] 1.000000e+00 1.141154e-08 5.705772e-09
[416,] 1.000000e+00 1.894118e-08 9.470591e-09
[417,] 1.000000e+00 3.211542e-08 1.605771e-08
[418,] 1.000000e+00 5.461797e-08 2.730898e-08
[419,] 1.000000e+00 9.476821e-08 4.738411e-08
[420,] 9.999999e-01 1.675822e-07 8.379109e-08
[421,] 9.999999e-01 2.915640e-07 1.457820e-07
[422,] 9.999997e-01 5.064041e-07 2.532021e-07
[423,] 9.999998e-01 4.906624e-07 2.453312e-07
[424,] 9.999998e-01 3.465254e-07 1.732627e-07
[425,] 9.999999e-01 2.179999e-07 1.090000e-07
[426,] 9.999999e-01 1.955073e-07 9.775367e-08
[427,] 9.999999e-01 2.449139e-07 1.224569e-07
[428,] 9.999999e-01 2.778229e-07 1.389115e-07
[429,] 9.999998e-01 3.241172e-07 1.620586e-07
[430,] 9.999998e-01 3.806014e-07 1.903007e-07
[431,] 9.999997e-01 5.207489e-07 2.603744e-07
[432,] 9.999996e-01 7.615264e-07 3.807632e-07
[433,] 9.999995e-01 1.095616e-06 5.478078e-07
[434,] 9.999992e-01 1.522158e-06 7.610790e-07
[435,] 9.999998e-01 4.058126e-07 2.029063e-07
[436,] 1.000000e+00 4.067531e-08 2.033766e-08
[437,] 1.000000e+00 2.143614e-09 1.071807e-09
[438,] 1.000000e+00 4.826014e-10 2.413007e-10
[439,] 1.000000e+00 2.830644e-10 1.415322e-10
[440,] 1.000000e+00 2.174966e-10 1.087483e-10
[441,] 1.000000e+00 1.075853e-10 5.379264e-11
[442,] 1.000000e+00 8.525459e-11 4.262729e-11
[443,] 1.000000e+00 1.322733e-10 6.613664e-11
[444,] 1.000000e+00 2.429422e-10 1.214711e-10
[445,] 1.000000e+00 5.905729e-10 2.952865e-10
[446,] 1.000000e+00 1.430138e-09 7.150692e-10
[447,] 1.000000e+00 7.512381e-10 3.756191e-10
[448,] 1.000000e+00 1.101343e-10 5.506717e-11
[449,] 1.000000e+00 6.342947e-11 3.171473e-11
[450,] 1.000000e+00 8.035128e-11 4.017564e-11
[451,] 1.000000e+00 1.946080e-10 9.730401e-11
[452,] 1.000000e+00 4.311444e-10 2.155722e-10
[453,] 1.000000e+00 8.813286e-10 4.406643e-10
[454,] 1.000000e+00 2.285371e-09 1.142685e-09
[455,] 1.000000e+00 6.055295e-09 3.027647e-09
[456,] 1.000000e+00 1.414386e-08 7.071932e-09
[457,] 1.000000e+00 1.006123e-08 5.030615e-09
[458,] 1.000000e+00 7.402192e-09 3.701096e-09
[459,] 1.000000e+00 2.116993e-08 1.058497e-08
[460,] 1.000000e+00 5.553417e-08 2.776709e-08
[461,] 9.999999e-01 1.523806e-07 7.619031e-08
[462,] 9.999999e-01 2.790894e-07 1.395447e-07
[463,] 9.999998e-01 3.167850e-07 1.583925e-07
[464,] 9.999997e-01 5.220112e-07 2.610056e-07
[465,] 9.999994e-01 1.121355e-06 5.606774e-07
[466,] 9.999988e-01 2.450959e-06 1.225480e-06
[467,] 9.999974e-01 5.118808e-06 2.559404e-06
[468,] 9.999953e-01 9.356194e-06 4.678097e-06
[469,] 9.999974e-01 5.180996e-06 2.590498e-06
[470,] 9.999985e-01 3.058491e-06 1.529246e-06
[471,] 9.999953e-01 9.310892e-06 4.655446e-06
[472,] 9.999916e-01 1.680188e-05 8.400941e-06
[473,] 9.999764e-01 4.718509e-05 2.359254e-05
[474,] 9.999258e-01 1.484970e-04 7.424848e-05
[475,] 9.997994e-01 4.011644e-04 2.005822e-04
[476,] 9.994092e-01 1.181691e-03 5.908453e-04
[477,] 9.984956e-01 3.008757e-03 1.504378e-03
[478,] 9.962076e-01 7.584743e-03 3.792371e-03
[479,] 9.899139e-01 2.017230e-02 1.008615e-02
[480,] 9.766000e-01 4.679998e-02 2.339999e-02
[481,] 9.764616e-01 4.707675e-02 2.353838e-02
[482,] 9.994118e-01 1.176329e-03 5.881647e-04
[483,] 9.970489e-01 5.902139e-03 2.951069e-03
> postscript(file="/var/yougetitorg/rcomp/tmp/1tl341296725581.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/yougetitorg/rcomp/tmp/2kemt1296725581.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/yougetitorg/rcomp/tmp/3nrt31296725581.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/yougetitorg/rcomp/tmp/4hba91296725581.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/yougetitorg/rcomp/tmp/5mzeo1296725581.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 = 492
Frequency = 1
1 2 3 4 5
-150.40547338 -163.50767485 -168.60987632 -172.71207779 -176.81427927
6 7 8 9 10
-177.91648074 -115.01868221 -95.12088369 -88.22308516 -89.32528663
11 12 13 14 15
-90.42748810 -85.52968958 -86.63189105 -89.73409252 -92.83629400
16 17 18 19 20
-90.93849547 -90.04069694 -85.14289841 -22.24509989 -3.34730136
21 22 23 24 25
0.55049717 -3.55170431 -2.65390578 4.24389275 8.14169128
26 27 28 29 30
7.03948980 -0.06271167 -7.16491314 -12.26711461 -12.36931609
31 32 33 34 35
50.52848244 67.42628097 70.32407949 69.22187802 72.11967655
36 37 38 39 40
77.01747508 76.91527360 70.81307213 59.71087066 52.60866918
41 42 43 44 45
42.50646771 41.40426624 101.30206477 124.19986329 127.09766182
46 47 48 49 50
121.99546035 109.89325887 101.79105740 97.68885593 84.58665446
51 52 53 54 55
69.48445298 53.38225151 47.28005004 41.17784856 102.07564709
56 57 58 59 60
130.97344562 125.87124415 117.76904267 112.66684120 99.56463973
61 62 63 64 65
113.46243825 90.36023678 64.25803531 40.15583384 22.05363236
66 67 68 69 70
-0.04856911 48.84922942 61.74702794 57.64482647 46.54262500
71 72 73 74 75
34.44042353 35.33822205 33.23602058 14.13381911 -1.96838237
76 77 78 79 80
-16.07058384 -25.17278531 -29.27498678 15.62281174 23.52061027
81 82 83 84 85
24.41840880 32.31620732 17.21400585 25.11180438 23.00960291
86 87 88 89 90
-0.09259857 -13.19480004 -26.29700151 -38.39920299 -42.50140446
91 92 93 94 95
6.39639407 9.29419260 7.19199112 2.08978965 -9.01241182
96 97 98 99 100
-11.11461330 -15.21681477 -33.31901624 -51.42121771 -66.52341919
101 102 103 104 105
-79.62562066 -85.72782213 -33.83002361 -28.93222508 -40.03442655
106 107 108 109 110
-50.13662802 -59.23882950 -60.34103097 -66.44323244 -81.54543392
111 112 113 114 115
-89.64763539 -110.74983686 -121.85203833 -126.95423981 -79.05644128
116 117 118 119 120
-65.15864275 -80.26084423 -88.36304570 -95.46524717 -95.56744864
121 122 123 124 125
-98.66965012 -112.77185159 -125.87405306 -137.97625454 -147.07845601
126 127 128 129 130
-147.18065748 -94.28285895 -84.38506043 -88.48726190 -91.58946337
131 132 133 134 135
-95.69166485 -89.79386632 -89.89606779 -93.99826926 -107.10047074
136 137 138 139 140
-115.20267221 -123.30487368 -122.40707516 -69.50927663 -53.61147810
141 142 143 144 145
-58.71367957 -57.81588105 -61.91808252 -55.02028399 -53.12248547
146 147 148 149 150
-63.22468694 -71.32688841 -75.42908988 -84.53129136 -78.63349283
151 152 153 154 155
-26.73569430 -8.83789578 -12.94009725 -14.04229872 -14.14450019
156 157 158 159 160
-0.24670167 5.65109686 4.54889539 -3.55330609 -6.65550756
161 162 163 164 165
-10.75770903 -4.85991050 54.03788802 75.93568655 77.83348508
166 167 168 169 170
69.73128360 69.62908213 75.52688066 77.42467919 77.32247771
171 172 173 174 175
69.22027624 69.11807477 55.01587329 54.91367182 45.81147035
176 177 178 179 180
45.70926888 37.60706740 37.50486593 36.40266446 36.30046299
181 182 183 184 185
83.19826151 83.09606004 103.99385857 103.89165709 98.78945562
186 187 188 189 190
98.68725415 87.58505268 87.48285120 76.38064973 76.27844826
191 192 193 194 195
80.17624678 80.07404531 80.97184384 80.86964237 72.76744089
196 197 198 199 200
72.66523942 55.56303795 55.46083647 47.35863500 47.25643353
201 202 203 204 205
35.15423206 35.05203058 34.94982911 34.84762764 87.74542616
206 207 208 209 210
87.64322469 110.54102322 110.43882175 108.33662027 108.23441880
211 212 213 214 215
103.13221733 103.03001585 88.92781438 88.82561291 94.72341144
216 217 218 219 220
94.62120996 99.51900849 99.41680702 82.31460554 82.21240407
221 222 223 224 225
62.11020260 62.00800113 51.90579965 51.80359818 32.70139671
226 227 228 229 230
32.59919523 22.49699376 22.39479229 76.29259082 76.19038934
231 232 233 234 235
97.08818787 96.98598640 84.88378492 84.78158345 66.67938198
236 237 238 239 240
66.57718051 55.47497903 55.37277756 56.27057609 56.16837461
241 242 243 244 245
55.06617314 54.96397167 47.86177020 47.75956872 35.65736725
246 247 248 249 250
35.55516578 26.45296430 26.35076283 12.24856136 12.14635989
251 252 253 254 255
16.04415841 15.94195694 62.83975547 62.73755399 81.63535252
256 257 258 259 260
81.53315105 73.43094958 73.32874810 52.22654663 52.12434516
261 262 263 264 265
37.02214368 36.91994221 39.81774074 39.71553927 34.61333779
266 267 268 269 270
34.51113632 27.40893485 27.30673337 13.20453190 13.10233043
271 272 273 274 275
-0.99987104 -1.10207252 -15.20427399 -15.30647546 -22.40867694
276 277 278 279 280
-22.51087841 26.38692012 26.28471865 56.18251717 56.08031570
281 282 283 284 285
40.97811423 40.87591275 11.77371128 11.67150981 -6.43069166
286 287 288 289 290
-6.53289314 -2.63509461 -2.73729608 -1.83949756 -1.94169903
291 292 293 294 295
-13.04390050 -13.14610197 -26.24830345 -26.35050492 -40.45270639
296 297 298 299 300
-40.55490787 -50.65710934 -50.75931081 -54.86151228 -54.96371376
301 302 303 304 305
-6.06591523 -6.16811670 21.72968182 21.62748035 1.52527888
306 307 308 309 310
1.42307741 -19.67912407 -19.78132554 -38.88352701 -38.98572849
311 312 313 314 315
-39.08792996 -39.19013143 -38.29233290 -38.39453438 -50.49673585
316 317 318 319 320
-50.59893732 -69.70113880 -69.80334027 -74.90554174 -75.00774321
321 322 323 324 325
-90.10994469 -90.21214616 -93.31434763 -93.41654911 -44.51875058
326 327 328 329 330
-44.62095205 -11.72315352 -11.82535500 -31.92755647 -32.02975794
331 332 333 334 335
-48.13195941 -48.23416089 -69.33636236 -69.43856383 -73.54076531
336 337 338 339 340
-73.64296678 -66.74516825 -66.84736972 -73.94957120 -74.05177267
341 342 343 344 345
-86.15397414 -86.25617562 -91.35837709 -91.46057856 -98.56278003
346 347 348 349 350
-98.66498151 -103.76718298 -103.86938445 -50.97158593 -51.07378740
351 352 353 354 355
-25.17598887 -25.27819034 -22.38039182 -22.48259329 -32.58479476
356 357 358 359 360
-32.68699624 -64.78919771 -64.89139918 -64.99360065 -65.09580213
361 362 363 364 365
-60.19800360 -60.30020507 -61.40240655 -61.50460802 -66.60680949
366 367 368 369 370
-66.70901096 -75.81121244 -75.91341391 -82.01561538 -82.11781686
371 372 373 374 375
-81.22001833 -81.32221980 -20.42442127 -20.52662275 -12.62882422
376 377 378 379 380
-12.73102569 -14.83322717 -14.93542864 -19.03763011 -19.13983158
381 382 383 384 385
-29.24203306 -29.34423453 -26.44643600 -26.54863748 -19.65083895
386 387 388 389 390
-19.75304042 -21.85524189 -21.95744337 -29.05964484 -29.16184631
391 392 393 394 395
-30.26404779 -30.36624926 -38.46845073 -38.57065220 -32.67285368
396 397 398 399 400
-32.77505515 29.12274338 29.02054190 39.91834043 39.81613896
401 402 403 404 405
37.71393749 37.61173601 24.50953454 24.40733307 6.30513159
406 407 408 409 410
6.20293012 14.10072865 13.99852718 20.89632570 19.79412423
411 412 413 414 415
13.69192276 2.58972128 -4.51248019 1.38531834 52.28311687
416 417 418 419 420
69.18091539 71.07871392 68.97651245 51.87431097 52.77210950
421 422 423 424 425
48.66990803 46.56770656 41.46550508 30.36330361 24.26110214
426 427 428 429 430
26.15890066 78.05669919 85.95449772 84.85229625 68.75009477
431 432 433 434 435
51.64789330 53.54569183 49.44349035 46.34128888 36.23908741
436 437 438 439 440
30.13688594 29.03468446 28.93248299 75.83028152 81.72808004
441 442 443 444 445
75.62587857 43.52367710 21.42147563 12.31927415 16.21707268
446 447 448 449 450
4.11487121 -12.98733027 -19.08953174 -34.19173321 -46.29393468
451 452 453 454 455
9.60386384 19.50166237 -3.60053910 -18.70274058 -35.80494205
456 457 458 459 460
-31.90714352 -29.00934499 -38.11154647 -53.21374794 -56.31594941
461 462 463 464 465
-77.41815089 -68.52035236 -18.62255383 -12.72475530 -28.82695678
466 467 468 469 470
-40.92915825 -45.03135972 -31.13356120 -19.23576267 -14.33796414
471 472 473 474 475
-11.44016561 -10.54236709 -23.64456856 -11.74677003 39.15102849
476 477 478 479 480
49.04882702 32.94662555 15.84442408 9.74222260 26.64002113
481 482 483 484 485
31.53781966 26.43561818 14.33341671 3.23121524 -11.87098623
486 487 488 489 490
-3.97318771 51.92461082 54.82240935 36.72020788 14.61800640
491 492
-0.48419507 1.41360346
> postscript(file="/var/yougetitorg/rcomp/tmp/6q5jy1296725581.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 = 492
Frequency = 1
lag(myerror, k = 1) myerror
0 -150.40547338 NA
1 -163.50767485 -150.40547338
2 -168.60987632 -163.50767485
3 -172.71207779 -168.60987632
4 -176.81427927 -172.71207779
5 -177.91648074 -176.81427927
6 -115.01868221 -177.91648074
7 -95.12088369 -115.01868221
8 -88.22308516 -95.12088369
9 -89.32528663 -88.22308516
10 -90.42748810 -89.32528663
11 -85.52968958 -90.42748810
12 -86.63189105 -85.52968958
13 -89.73409252 -86.63189105
14 -92.83629400 -89.73409252
15 -90.93849547 -92.83629400
16 -90.04069694 -90.93849547
17 -85.14289841 -90.04069694
18 -22.24509989 -85.14289841
19 -3.34730136 -22.24509989
20 0.55049717 -3.34730136
21 -3.55170431 0.55049717
22 -2.65390578 -3.55170431
23 4.24389275 -2.65390578
24 8.14169128 4.24389275
25 7.03948980 8.14169128
26 -0.06271167 7.03948980
27 -7.16491314 -0.06271167
28 -12.26711461 -7.16491314
29 -12.36931609 -12.26711461
30 50.52848244 -12.36931609
31 67.42628097 50.52848244
32 70.32407949 67.42628097
33 69.22187802 70.32407949
34 72.11967655 69.22187802
35 77.01747508 72.11967655
36 76.91527360 77.01747508
37 70.81307213 76.91527360
38 59.71087066 70.81307213
39 52.60866918 59.71087066
40 42.50646771 52.60866918
41 41.40426624 42.50646771
42 101.30206477 41.40426624
43 124.19986329 101.30206477
44 127.09766182 124.19986329
45 121.99546035 127.09766182
46 109.89325887 121.99546035
47 101.79105740 109.89325887
48 97.68885593 101.79105740
49 84.58665446 97.68885593
50 69.48445298 84.58665446
51 53.38225151 69.48445298
52 47.28005004 53.38225151
53 41.17784856 47.28005004
54 102.07564709 41.17784856
55 130.97344562 102.07564709
56 125.87124415 130.97344562
57 117.76904267 125.87124415
58 112.66684120 117.76904267
59 99.56463973 112.66684120
60 113.46243825 99.56463973
61 90.36023678 113.46243825
62 64.25803531 90.36023678
63 40.15583384 64.25803531
64 22.05363236 40.15583384
65 -0.04856911 22.05363236
66 48.84922942 -0.04856911
67 61.74702794 48.84922942
68 57.64482647 61.74702794
69 46.54262500 57.64482647
70 34.44042353 46.54262500
71 35.33822205 34.44042353
72 33.23602058 35.33822205
73 14.13381911 33.23602058
74 -1.96838237 14.13381911
75 -16.07058384 -1.96838237
76 -25.17278531 -16.07058384
77 -29.27498678 -25.17278531
78 15.62281174 -29.27498678
79 23.52061027 15.62281174
80 24.41840880 23.52061027
81 32.31620732 24.41840880
82 17.21400585 32.31620732
83 25.11180438 17.21400585
84 23.00960291 25.11180438
85 -0.09259857 23.00960291
86 -13.19480004 -0.09259857
87 -26.29700151 -13.19480004
88 -38.39920299 -26.29700151
89 -42.50140446 -38.39920299
90 6.39639407 -42.50140446
91 9.29419260 6.39639407
92 7.19199112 9.29419260
93 2.08978965 7.19199112
94 -9.01241182 2.08978965
95 -11.11461330 -9.01241182
96 -15.21681477 -11.11461330
97 -33.31901624 -15.21681477
98 -51.42121771 -33.31901624
99 -66.52341919 -51.42121771
100 -79.62562066 -66.52341919
101 -85.72782213 -79.62562066
102 -33.83002361 -85.72782213
103 -28.93222508 -33.83002361
104 -40.03442655 -28.93222508
105 -50.13662802 -40.03442655
106 -59.23882950 -50.13662802
107 -60.34103097 -59.23882950
108 -66.44323244 -60.34103097
109 -81.54543392 -66.44323244
110 -89.64763539 -81.54543392
111 -110.74983686 -89.64763539
112 -121.85203833 -110.74983686
113 -126.95423981 -121.85203833
114 -79.05644128 -126.95423981
115 -65.15864275 -79.05644128
116 -80.26084423 -65.15864275
117 -88.36304570 -80.26084423
118 -95.46524717 -88.36304570
119 -95.56744864 -95.46524717
120 -98.66965012 -95.56744864
121 -112.77185159 -98.66965012
122 -125.87405306 -112.77185159
123 -137.97625454 -125.87405306
124 -147.07845601 -137.97625454
125 -147.18065748 -147.07845601
126 -94.28285895 -147.18065748
127 -84.38506043 -94.28285895
128 -88.48726190 -84.38506043
129 -91.58946337 -88.48726190
130 -95.69166485 -91.58946337
131 -89.79386632 -95.69166485
132 -89.89606779 -89.79386632
133 -93.99826926 -89.89606779
134 -107.10047074 -93.99826926
135 -115.20267221 -107.10047074
136 -123.30487368 -115.20267221
137 -122.40707516 -123.30487368
138 -69.50927663 -122.40707516
139 -53.61147810 -69.50927663
140 -58.71367957 -53.61147810
141 -57.81588105 -58.71367957
142 -61.91808252 -57.81588105
143 -55.02028399 -61.91808252
144 -53.12248547 -55.02028399
145 -63.22468694 -53.12248547
146 -71.32688841 -63.22468694
147 -75.42908988 -71.32688841
148 -84.53129136 -75.42908988
149 -78.63349283 -84.53129136
150 -26.73569430 -78.63349283
151 -8.83789578 -26.73569430
152 -12.94009725 -8.83789578
153 -14.04229872 -12.94009725
154 -14.14450019 -14.04229872
155 -0.24670167 -14.14450019
156 5.65109686 -0.24670167
157 4.54889539 5.65109686
158 -3.55330609 4.54889539
159 -6.65550756 -3.55330609
160 -10.75770903 -6.65550756
161 -4.85991050 -10.75770903
162 54.03788802 -4.85991050
163 75.93568655 54.03788802
164 77.83348508 75.93568655
165 69.73128360 77.83348508
166 69.62908213 69.73128360
167 75.52688066 69.62908213
168 77.42467919 75.52688066
169 77.32247771 77.42467919
170 69.22027624 77.32247771
171 69.11807477 69.22027624
172 55.01587329 69.11807477
173 54.91367182 55.01587329
174 45.81147035 54.91367182
175 45.70926888 45.81147035
176 37.60706740 45.70926888
177 37.50486593 37.60706740
178 36.40266446 37.50486593
179 36.30046299 36.40266446
180 83.19826151 36.30046299
181 83.09606004 83.19826151
182 103.99385857 83.09606004
183 103.89165709 103.99385857
184 98.78945562 103.89165709
185 98.68725415 98.78945562
186 87.58505268 98.68725415
187 87.48285120 87.58505268
188 76.38064973 87.48285120
189 76.27844826 76.38064973
190 80.17624678 76.27844826
191 80.07404531 80.17624678
192 80.97184384 80.07404531
193 80.86964237 80.97184384
194 72.76744089 80.86964237
195 72.66523942 72.76744089
196 55.56303795 72.66523942
197 55.46083647 55.56303795
198 47.35863500 55.46083647
199 47.25643353 47.35863500
200 35.15423206 47.25643353
201 35.05203058 35.15423206
202 34.94982911 35.05203058
203 34.84762764 34.94982911
204 87.74542616 34.84762764
205 87.64322469 87.74542616
206 110.54102322 87.64322469
207 110.43882175 110.54102322
208 108.33662027 110.43882175
209 108.23441880 108.33662027
210 103.13221733 108.23441880
211 103.03001585 103.13221733
212 88.92781438 103.03001585
213 88.82561291 88.92781438
214 94.72341144 88.82561291
215 94.62120996 94.72341144
216 99.51900849 94.62120996
217 99.41680702 99.51900849
218 82.31460554 99.41680702
219 82.21240407 82.31460554
220 62.11020260 82.21240407
221 62.00800113 62.11020260
222 51.90579965 62.00800113
223 51.80359818 51.90579965
224 32.70139671 51.80359818
225 32.59919523 32.70139671
226 22.49699376 32.59919523
227 22.39479229 22.49699376
228 76.29259082 22.39479229
229 76.19038934 76.29259082
230 97.08818787 76.19038934
231 96.98598640 97.08818787
232 84.88378492 96.98598640
233 84.78158345 84.88378492
234 66.67938198 84.78158345
235 66.57718051 66.67938198
236 55.47497903 66.57718051
237 55.37277756 55.47497903
238 56.27057609 55.37277756
239 56.16837461 56.27057609
240 55.06617314 56.16837461
241 54.96397167 55.06617314
242 47.86177020 54.96397167
243 47.75956872 47.86177020
244 35.65736725 47.75956872
245 35.55516578 35.65736725
246 26.45296430 35.55516578
247 26.35076283 26.45296430
248 12.24856136 26.35076283
249 12.14635989 12.24856136
250 16.04415841 12.14635989
251 15.94195694 16.04415841
252 62.83975547 15.94195694
253 62.73755399 62.83975547
254 81.63535252 62.73755399
255 81.53315105 81.63535252
256 73.43094958 81.53315105
257 73.32874810 73.43094958
258 52.22654663 73.32874810
259 52.12434516 52.22654663
260 37.02214368 52.12434516
261 36.91994221 37.02214368
262 39.81774074 36.91994221
263 39.71553927 39.81774074
264 34.61333779 39.71553927
265 34.51113632 34.61333779
266 27.40893485 34.51113632
267 27.30673337 27.40893485
268 13.20453190 27.30673337
269 13.10233043 13.20453190
270 -0.99987104 13.10233043
271 -1.10207252 -0.99987104
272 -15.20427399 -1.10207252
273 -15.30647546 -15.20427399
274 -22.40867694 -15.30647546
275 -22.51087841 -22.40867694
276 26.38692012 -22.51087841
277 26.28471865 26.38692012
278 56.18251717 26.28471865
279 56.08031570 56.18251717
280 40.97811423 56.08031570
281 40.87591275 40.97811423
282 11.77371128 40.87591275
283 11.67150981 11.77371128
284 -6.43069166 11.67150981
285 -6.53289314 -6.43069166
286 -2.63509461 -6.53289314
287 -2.73729608 -2.63509461
288 -1.83949756 -2.73729608
289 -1.94169903 -1.83949756
290 -13.04390050 -1.94169903
291 -13.14610197 -13.04390050
292 -26.24830345 -13.14610197
293 -26.35050492 -26.24830345
294 -40.45270639 -26.35050492
295 -40.55490787 -40.45270639
296 -50.65710934 -40.55490787
297 -50.75931081 -50.65710934
298 -54.86151228 -50.75931081
299 -54.96371376 -54.86151228
300 -6.06591523 -54.96371376
301 -6.16811670 -6.06591523
302 21.72968182 -6.16811670
303 21.62748035 21.72968182
304 1.52527888 21.62748035
305 1.42307741 1.52527888
306 -19.67912407 1.42307741
307 -19.78132554 -19.67912407
308 -38.88352701 -19.78132554
309 -38.98572849 -38.88352701
310 -39.08792996 -38.98572849
311 -39.19013143 -39.08792996
312 -38.29233290 -39.19013143
313 -38.39453438 -38.29233290
314 -50.49673585 -38.39453438
315 -50.59893732 -50.49673585
316 -69.70113880 -50.59893732
317 -69.80334027 -69.70113880
318 -74.90554174 -69.80334027
319 -75.00774321 -74.90554174
320 -90.10994469 -75.00774321
321 -90.21214616 -90.10994469
322 -93.31434763 -90.21214616
323 -93.41654911 -93.31434763
324 -44.51875058 -93.41654911
325 -44.62095205 -44.51875058
326 -11.72315352 -44.62095205
327 -11.82535500 -11.72315352
328 -31.92755647 -11.82535500
329 -32.02975794 -31.92755647
330 -48.13195941 -32.02975794
331 -48.23416089 -48.13195941
332 -69.33636236 -48.23416089
333 -69.43856383 -69.33636236
334 -73.54076531 -69.43856383
335 -73.64296678 -73.54076531
336 -66.74516825 -73.64296678
337 -66.84736972 -66.74516825
338 -73.94957120 -66.84736972
339 -74.05177267 -73.94957120
340 -86.15397414 -74.05177267
341 -86.25617562 -86.15397414
342 -91.35837709 -86.25617562
343 -91.46057856 -91.35837709
344 -98.56278003 -91.46057856
345 -98.66498151 -98.56278003
346 -103.76718298 -98.66498151
347 -103.86938445 -103.76718298
348 -50.97158593 -103.86938445
349 -51.07378740 -50.97158593
350 -25.17598887 -51.07378740
351 -25.27819034 -25.17598887
352 -22.38039182 -25.27819034
353 -22.48259329 -22.38039182
354 -32.58479476 -22.48259329
355 -32.68699624 -32.58479476
356 -64.78919771 -32.68699624
357 -64.89139918 -64.78919771
358 -64.99360065 -64.89139918
359 -65.09580213 -64.99360065
360 -60.19800360 -65.09580213
361 -60.30020507 -60.19800360
362 -61.40240655 -60.30020507
363 -61.50460802 -61.40240655
364 -66.60680949 -61.50460802
365 -66.70901096 -66.60680949
366 -75.81121244 -66.70901096
367 -75.91341391 -75.81121244
368 -82.01561538 -75.91341391
369 -82.11781686 -82.01561538
370 -81.22001833 -82.11781686
371 -81.32221980 -81.22001833
372 -20.42442127 -81.32221980
373 -20.52662275 -20.42442127
374 -12.62882422 -20.52662275
375 -12.73102569 -12.62882422
376 -14.83322717 -12.73102569
377 -14.93542864 -14.83322717
378 -19.03763011 -14.93542864
379 -19.13983158 -19.03763011
380 -29.24203306 -19.13983158
381 -29.34423453 -29.24203306
382 -26.44643600 -29.34423453
383 -26.54863748 -26.44643600
384 -19.65083895 -26.54863748
385 -19.75304042 -19.65083895
386 -21.85524189 -19.75304042
387 -21.95744337 -21.85524189
388 -29.05964484 -21.95744337
389 -29.16184631 -29.05964484
390 -30.26404779 -29.16184631
391 -30.36624926 -30.26404779
392 -38.46845073 -30.36624926
393 -38.57065220 -38.46845073
394 -32.67285368 -38.57065220
395 -32.77505515 -32.67285368
396 29.12274338 -32.77505515
397 29.02054190 29.12274338
398 39.91834043 29.02054190
399 39.81613896 39.91834043
400 37.71393749 39.81613896
401 37.61173601 37.71393749
402 24.50953454 37.61173601
403 24.40733307 24.50953454
404 6.30513159 24.40733307
405 6.20293012 6.30513159
406 14.10072865 6.20293012
407 13.99852718 14.10072865
408 20.89632570 13.99852718
409 19.79412423 20.89632570
410 13.69192276 19.79412423
411 2.58972128 13.69192276
412 -4.51248019 2.58972128
413 1.38531834 -4.51248019
414 52.28311687 1.38531834
415 69.18091539 52.28311687
416 71.07871392 69.18091539
417 68.97651245 71.07871392
418 51.87431097 68.97651245
419 52.77210950 51.87431097
420 48.66990803 52.77210950
421 46.56770656 48.66990803
422 41.46550508 46.56770656
423 30.36330361 41.46550508
424 24.26110214 30.36330361
425 26.15890066 24.26110214
426 78.05669919 26.15890066
427 85.95449772 78.05669919
428 84.85229625 85.95449772
429 68.75009477 84.85229625
430 51.64789330 68.75009477
431 53.54569183 51.64789330
432 49.44349035 53.54569183
433 46.34128888 49.44349035
434 36.23908741 46.34128888
435 30.13688594 36.23908741
436 29.03468446 30.13688594
437 28.93248299 29.03468446
438 75.83028152 28.93248299
439 81.72808004 75.83028152
440 75.62587857 81.72808004
441 43.52367710 75.62587857
442 21.42147563 43.52367710
443 12.31927415 21.42147563
444 16.21707268 12.31927415
445 4.11487121 16.21707268
446 -12.98733027 4.11487121
447 -19.08953174 -12.98733027
448 -34.19173321 -19.08953174
449 -46.29393468 -34.19173321
450 9.60386384 -46.29393468
451 19.50166237 9.60386384
452 -3.60053910 19.50166237
453 -18.70274058 -3.60053910
454 -35.80494205 -18.70274058
455 -31.90714352 -35.80494205
456 -29.00934499 -31.90714352
457 -38.11154647 -29.00934499
458 -53.21374794 -38.11154647
459 -56.31594941 -53.21374794
460 -77.41815089 -56.31594941
461 -68.52035236 -77.41815089
462 -18.62255383 -68.52035236
463 -12.72475530 -18.62255383
464 -28.82695678 -12.72475530
465 -40.92915825 -28.82695678
466 -45.03135972 -40.92915825
467 -31.13356120 -45.03135972
468 -19.23576267 -31.13356120
469 -14.33796414 -19.23576267
470 -11.44016561 -14.33796414
471 -10.54236709 -11.44016561
472 -23.64456856 -10.54236709
473 -11.74677003 -23.64456856
474 39.15102849 -11.74677003
475 49.04882702 39.15102849
476 32.94662555 49.04882702
477 15.84442408 32.94662555
478 9.74222260 15.84442408
479 26.64002113 9.74222260
480 31.53781966 26.64002113
481 26.43561818 31.53781966
482 14.33341671 26.43561818
483 3.23121524 14.33341671
484 -11.87098623 3.23121524
485 -3.97318771 -11.87098623
486 51.92461082 -3.97318771
487 54.82240935 51.92461082
488 36.72020788 54.82240935
489 14.61800640 36.72020788
490 -0.48419507 14.61800640
491 1.41360346 -0.48419507
492 NA 1.41360346
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -163.50767485 -150.40547338
[2,] -168.60987632 -163.50767485
[3,] -172.71207779 -168.60987632
[4,] -176.81427927 -172.71207779
[5,] -177.91648074 -176.81427927
[6,] -115.01868221 -177.91648074
[7,] -95.12088369 -115.01868221
[8,] -88.22308516 -95.12088369
[9,] -89.32528663 -88.22308516
[10,] -90.42748810 -89.32528663
[11,] -85.52968958 -90.42748810
[12,] -86.63189105 -85.52968958
[13,] -89.73409252 -86.63189105
[14,] -92.83629400 -89.73409252
[15,] -90.93849547 -92.83629400
[16,] -90.04069694 -90.93849547
[17,] -85.14289841 -90.04069694
[18,] -22.24509989 -85.14289841
[19,] -3.34730136 -22.24509989
[20,] 0.55049717 -3.34730136
[21,] -3.55170431 0.55049717
[22,] -2.65390578 -3.55170431
[23,] 4.24389275 -2.65390578
[24,] 8.14169128 4.24389275
[25,] 7.03948980 8.14169128
[26,] -0.06271167 7.03948980
[27,] -7.16491314 -0.06271167
[28,] -12.26711461 -7.16491314
[29,] -12.36931609 -12.26711461
[30,] 50.52848244 -12.36931609
[31,] 67.42628097 50.52848244
[32,] 70.32407949 67.42628097
[33,] 69.22187802 70.32407949
[34,] 72.11967655 69.22187802
[35,] 77.01747508 72.11967655
[36,] 76.91527360 77.01747508
[37,] 70.81307213 76.91527360
[38,] 59.71087066 70.81307213
[39,] 52.60866918 59.71087066
[40,] 42.50646771 52.60866918
[41,] 41.40426624 42.50646771
[42,] 101.30206477 41.40426624
[43,] 124.19986329 101.30206477
[44,] 127.09766182 124.19986329
[45,] 121.99546035 127.09766182
[46,] 109.89325887 121.99546035
[47,] 101.79105740 109.89325887
[48,] 97.68885593 101.79105740
[49,] 84.58665446 97.68885593
[50,] 69.48445298 84.58665446
[51,] 53.38225151 69.48445298
[52,] 47.28005004 53.38225151
[53,] 41.17784856 47.28005004
[54,] 102.07564709 41.17784856
[55,] 130.97344562 102.07564709
[56,] 125.87124415 130.97344562
[57,] 117.76904267 125.87124415
[58,] 112.66684120 117.76904267
[59,] 99.56463973 112.66684120
[60,] 113.46243825 99.56463973
[61,] 90.36023678 113.46243825
[62,] 64.25803531 90.36023678
[63,] 40.15583384 64.25803531
[64,] 22.05363236 40.15583384
[65,] -0.04856911 22.05363236
[66,] 48.84922942 -0.04856911
[67,] 61.74702794 48.84922942
[68,] 57.64482647 61.74702794
[69,] 46.54262500 57.64482647
[70,] 34.44042353 46.54262500
[71,] 35.33822205 34.44042353
[72,] 33.23602058 35.33822205
[73,] 14.13381911 33.23602058
[74,] -1.96838237 14.13381911
[75,] -16.07058384 -1.96838237
[76,] -25.17278531 -16.07058384
[77,] -29.27498678 -25.17278531
[78,] 15.62281174 -29.27498678
[79,] 23.52061027 15.62281174
[80,] 24.41840880 23.52061027
[81,] 32.31620732 24.41840880
[82,] 17.21400585 32.31620732
[83,] 25.11180438 17.21400585
[84,] 23.00960291 25.11180438
[85,] -0.09259857 23.00960291
[86,] -13.19480004 -0.09259857
[87,] -26.29700151 -13.19480004
[88,] -38.39920299 -26.29700151
[89,] -42.50140446 -38.39920299
[90,] 6.39639407 -42.50140446
[91,] 9.29419260 6.39639407
[92,] 7.19199112 9.29419260
[93,] 2.08978965 7.19199112
[94,] -9.01241182 2.08978965
[95,] -11.11461330 -9.01241182
[96,] -15.21681477 -11.11461330
[97,] -33.31901624 -15.21681477
[98,] -51.42121771 -33.31901624
[99,] -66.52341919 -51.42121771
[100,] -79.62562066 -66.52341919
[101,] -85.72782213 -79.62562066
[102,] -33.83002361 -85.72782213
[103,] -28.93222508 -33.83002361
[104,] -40.03442655 -28.93222508
[105,] -50.13662802 -40.03442655
[106,] -59.23882950 -50.13662802
[107,] -60.34103097 -59.23882950
[108,] -66.44323244 -60.34103097
[109,] -81.54543392 -66.44323244
[110,] -89.64763539 -81.54543392
[111,] -110.74983686 -89.64763539
[112,] -121.85203833 -110.74983686
[113,] -126.95423981 -121.85203833
[114,] -79.05644128 -126.95423981
[115,] -65.15864275 -79.05644128
[116,] -80.26084423 -65.15864275
[117,] -88.36304570 -80.26084423
[118,] -95.46524717 -88.36304570
[119,] -95.56744864 -95.46524717
[120,] -98.66965012 -95.56744864
[121,] -112.77185159 -98.66965012
[122,] -125.87405306 -112.77185159
[123,] -137.97625454 -125.87405306
[124,] -147.07845601 -137.97625454
[125,] -147.18065748 -147.07845601
[126,] -94.28285895 -147.18065748
[127,] -84.38506043 -94.28285895
[128,] -88.48726190 -84.38506043
[129,] -91.58946337 -88.48726190
[130,] -95.69166485 -91.58946337
[131,] -89.79386632 -95.69166485
[132,] -89.89606779 -89.79386632
[133,] -93.99826926 -89.89606779
[134,] -107.10047074 -93.99826926
[135,] -115.20267221 -107.10047074
[136,] -123.30487368 -115.20267221
[137,] -122.40707516 -123.30487368
[138,] -69.50927663 -122.40707516
[139,] -53.61147810 -69.50927663
[140,] -58.71367957 -53.61147810
[141,] -57.81588105 -58.71367957
[142,] -61.91808252 -57.81588105
[143,] -55.02028399 -61.91808252
[144,] -53.12248547 -55.02028399
[145,] -63.22468694 -53.12248547
[146,] -71.32688841 -63.22468694
[147,] -75.42908988 -71.32688841
[148,] -84.53129136 -75.42908988
[149,] -78.63349283 -84.53129136
[150,] -26.73569430 -78.63349283
[151,] -8.83789578 -26.73569430
[152,] -12.94009725 -8.83789578
[153,] -14.04229872 -12.94009725
[154,] -14.14450019 -14.04229872
[155,] -0.24670167 -14.14450019
[156,] 5.65109686 -0.24670167
[157,] 4.54889539 5.65109686
[158,] -3.55330609 4.54889539
[159,] -6.65550756 -3.55330609
[160,] -10.75770903 -6.65550756
[161,] -4.85991050 -10.75770903
[162,] 54.03788802 -4.85991050
[163,] 75.93568655 54.03788802
[164,] 77.83348508 75.93568655
[165,] 69.73128360 77.83348508
[166,] 69.62908213 69.73128360
[167,] 75.52688066 69.62908213
[168,] 77.42467919 75.52688066
[169,] 77.32247771 77.42467919
[170,] 69.22027624 77.32247771
[171,] 69.11807477 69.22027624
[172,] 55.01587329 69.11807477
[173,] 54.91367182 55.01587329
[174,] 45.81147035 54.91367182
[175,] 45.70926888 45.81147035
[176,] 37.60706740 45.70926888
[177,] 37.50486593 37.60706740
[178,] 36.40266446 37.50486593
[179,] 36.30046299 36.40266446
[180,] 83.19826151 36.30046299
[181,] 83.09606004 83.19826151
[182,] 103.99385857 83.09606004
[183,] 103.89165709 103.99385857
[184,] 98.78945562 103.89165709
[185,] 98.68725415 98.78945562
[186,] 87.58505268 98.68725415
[187,] 87.48285120 87.58505268
[188,] 76.38064973 87.48285120
[189,] 76.27844826 76.38064973
[190,] 80.17624678 76.27844826
[191,] 80.07404531 80.17624678
[192,] 80.97184384 80.07404531
[193,] 80.86964237 80.97184384
[194,] 72.76744089 80.86964237
[195,] 72.66523942 72.76744089
[196,] 55.56303795 72.66523942
[197,] 55.46083647 55.56303795
[198,] 47.35863500 55.46083647
[199,] 47.25643353 47.35863500
[200,] 35.15423206 47.25643353
[201,] 35.05203058 35.15423206
[202,] 34.94982911 35.05203058
[203,] 34.84762764 34.94982911
[204,] 87.74542616 34.84762764
[205,] 87.64322469 87.74542616
[206,] 110.54102322 87.64322469
[207,] 110.43882175 110.54102322
[208,] 108.33662027 110.43882175
[209,] 108.23441880 108.33662027
[210,] 103.13221733 108.23441880
[211,] 103.03001585 103.13221733
[212,] 88.92781438 103.03001585
[213,] 88.82561291 88.92781438
[214,] 94.72341144 88.82561291
[215,] 94.62120996 94.72341144
[216,] 99.51900849 94.62120996
[217,] 99.41680702 99.51900849
[218,] 82.31460554 99.41680702
[219,] 82.21240407 82.31460554
[220,] 62.11020260 82.21240407
[221,] 62.00800113 62.11020260
[222,] 51.90579965 62.00800113
[223,] 51.80359818 51.90579965
[224,] 32.70139671 51.80359818
[225,] 32.59919523 32.70139671
[226,] 22.49699376 32.59919523
[227,] 22.39479229 22.49699376
[228,] 76.29259082 22.39479229
[229,] 76.19038934 76.29259082
[230,] 97.08818787 76.19038934
[231,] 96.98598640 97.08818787
[232,] 84.88378492 96.98598640
[233,] 84.78158345 84.88378492
[234,] 66.67938198 84.78158345
[235,] 66.57718051 66.67938198
[236,] 55.47497903 66.57718051
[237,] 55.37277756 55.47497903
[238,] 56.27057609 55.37277756
[239,] 56.16837461 56.27057609
[240,] 55.06617314 56.16837461
[241,] 54.96397167 55.06617314
[242,] 47.86177020 54.96397167
[243,] 47.75956872 47.86177020
[244,] 35.65736725 47.75956872
[245,] 35.55516578 35.65736725
[246,] 26.45296430 35.55516578
[247,] 26.35076283 26.45296430
[248,] 12.24856136 26.35076283
[249,] 12.14635989 12.24856136
[250,] 16.04415841 12.14635989
[251,] 15.94195694 16.04415841
[252,] 62.83975547 15.94195694
[253,] 62.73755399 62.83975547
[254,] 81.63535252 62.73755399
[255,] 81.53315105 81.63535252
[256,] 73.43094958 81.53315105
[257,] 73.32874810 73.43094958
[258,] 52.22654663 73.32874810
[259,] 52.12434516 52.22654663
[260,] 37.02214368 52.12434516
[261,] 36.91994221 37.02214368
[262,] 39.81774074 36.91994221
[263,] 39.71553927 39.81774074
[264,] 34.61333779 39.71553927
[265,] 34.51113632 34.61333779
[266,] 27.40893485 34.51113632
[267,] 27.30673337 27.40893485
[268,] 13.20453190 27.30673337
[269,] 13.10233043 13.20453190
[270,] -0.99987104 13.10233043
[271,] -1.10207252 -0.99987104
[272,] -15.20427399 -1.10207252
[273,] -15.30647546 -15.20427399
[274,] -22.40867694 -15.30647546
[275,] -22.51087841 -22.40867694
[276,] 26.38692012 -22.51087841
[277,] 26.28471865 26.38692012
[278,] 56.18251717 26.28471865
[279,] 56.08031570 56.18251717
[280,] 40.97811423 56.08031570
[281,] 40.87591275 40.97811423
[282,] 11.77371128 40.87591275
[283,] 11.67150981 11.77371128
[284,] -6.43069166 11.67150981
[285,] -6.53289314 -6.43069166
[286,] -2.63509461 -6.53289314
[287,] -2.73729608 -2.63509461
[288,] -1.83949756 -2.73729608
[289,] -1.94169903 -1.83949756
[290,] -13.04390050 -1.94169903
[291,] -13.14610197 -13.04390050
[292,] -26.24830345 -13.14610197
[293,] -26.35050492 -26.24830345
[294,] -40.45270639 -26.35050492
[295,] -40.55490787 -40.45270639
[296,] -50.65710934 -40.55490787
[297,] -50.75931081 -50.65710934
[298,] -54.86151228 -50.75931081
[299,] -54.96371376 -54.86151228
[300,] -6.06591523 -54.96371376
[301,] -6.16811670 -6.06591523
[302,] 21.72968182 -6.16811670
[303,] 21.62748035 21.72968182
[304,] 1.52527888 21.62748035
[305,] 1.42307741 1.52527888
[306,] -19.67912407 1.42307741
[307,] -19.78132554 -19.67912407
[308,] -38.88352701 -19.78132554
[309,] -38.98572849 -38.88352701
[310,] -39.08792996 -38.98572849
[311,] -39.19013143 -39.08792996
[312,] -38.29233290 -39.19013143
[313,] -38.39453438 -38.29233290
[314,] -50.49673585 -38.39453438
[315,] -50.59893732 -50.49673585
[316,] -69.70113880 -50.59893732
[317,] -69.80334027 -69.70113880
[318,] -74.90554174 -69.80334027
[319,] -75.00774321 -74.90554174
[320,] -90.10994469 -75.00774321
[321,] -90.21214616 -90.10994469
[322,] -93.31434763 -90.21214616
[323,] -93.41654911 -93.31434763
[324,] -44.51875058 -93.41654911
[325,] -44.62095205 -44.51875058
[326,] -11.72315352 -44.62095205
[327,] -11.82535500 -11.72315352
[328,] -31.92755647 -11.82535500
[329,] -32.02975794 -31.92755647
[330,] -48.13195941 -32.02975794
[331,] -48.23416089 -48.13195941
[332,] -69.33636236 -48.23416089
[333,] -69.43856383 -69.33636236
[334,] -73.54076531 -69.43856383
[335,] -73.64296678 -73.54076531
[336,] -66.74516825 -73.64296678
[337,] -66.84736972 -66.74516825
[338,] -73.94957120 -66.84736972
[339,] -74.05177267 -73.94957120
[340,] -86.15397414 -74.05177267
[341,] -86.25617562 -86.15397414
[342,] -91.35837709 -86.25617562
[343,] -91.46057856 -91.35837709
[344,] -98.56278003 -91.46057856
[345,] -98.66498151 -98.56278003
[346,] -103.76718298 -98.66498151
[347,] -103.86938445 -103.76718298
[348,] -50.97158593 -103.86938445
[349,] -51.07378740 -50.97158593
[350,] -25.17598887 -51.07378740
[351,] -25.27819034 -25.17598887
[352,] -22.38039182 -25.27819034
[353,] -22.48259329 -22.38039182
[354,] -32.58479476 -22.48259329
[355,] -32.68699624 -32.58479476
[356,] -64.78919771 -32.68699624
[357,] -64.89139918 -64.78919771
[358,] -64.99360065 -64.89139918
[359,] -65.09580213 -64.99360065
[360,] -60.19800360 -65.09580213
[361,] -60.30020507 -60.19800360
[362,] -61.40240655 -60.30020507
[363,] -61.50460802 -61.40240655
[364,] -66.60680949 -61.50460802
[365,] -66.70901096 -66.60680949
[366,] -75.81121244 -66.70901096
[367,] -75.91341391 -75.81121244
[368,] -82.01561538 -75.91341391
[369,] -82.11781686 -82.01561538
[370,] -81.22001833 -82.11781686
[371,] -81.32221980 -81.22001833
[372,] -20.42442127 -81.32221980
[373,] -20.52662275 -20.42442127
[374,] -12.62882422 -20.52662275
[375,] -12.73102569 -12.62882422
[376,] -14.83322717 -12.73102569
[377,] -14.93542864 -14.83322717
[378,] -19.03763011 -14.93542864
[379,] -19.13983158 -19.03763011
[380,] -29.24203306 -19.13983158
[381,] -29.34423453 -29.24203306
[382,] -26.44643600 -29.34423453
[383,] -26.54863748 -26.44643600
[384,] -19.65083895 -26.54863748
[385,] -19.75304042 -19.65083895
[386,] -21.85524189 -19.75304042
[387,] -21.95744337 -21.85524189
[388,] -29.05964484 -21.95744337
[389,] -29.16184631 -29.05964484
[390,] -30.26404779 -29.16184631
[391,] -30.36624926 -30.26404779
[392,] -38.46845073 -30.36624926
[393,] -38.57065220 -38.46845073
[394,] -32.67285368 -38.57065220
[395,] -32.77505515 -32.67285368
[396,] 29.12274338 -32.77505515
[397,] 29.02054190 29.12274338
[398,] 39.91834043 29.02054190
[399,] 39.81613896 39.91834043
[400,] 37.71393749 39.81613896
[401,] 37.61173601 37.71393749
[402,] 24.50953454 37.61173601
[403,] 24.40733307 24.50953454
[404,] 6.30513159 24.40733307
[405,] 6.20293012 6.30513159
[406,] 14.10072865 6.20293012
[407,] 13.99852718 14.10072865
[408,] 20.89632570 13.99852718
[409,] 19.79412423 20.89632570
[410,] 13.69192276 19.79412423
[411,] 2.58972128 13.69192276
[412,] -4.51248019 2.58972128
[413,] 1.38531834 -4.51248019
[414,] 52.28311687 1.38531834
[415,] 69.18091539 52.28311687
[416,] 71.07871392 69.18091539
[417,] 68.97651245 71.07871392
[418,] 51.87431097 68.97651245
[419,] 52.77210950 51.87431097
[420,] 48.66990803 52.77210950
[421,] 46.56770656 48.66990803
[422,] 41.46550508 46.56770656
[423,] 30.36330361 41.46550508
[424,] 24.26110214 30.36330361
[425,] 26.15890066 24.26110214
[426,] 78.05669919 26.15890066
[427,] 85.95449772 78.05669919
[428,] 84.85229625 85.95449772
[429,] 68.75009477 84.85229625
[430,] 51.64789330 68.75009477
[431,] 53.54569183 51.64789330
[432,] 49.44349035 53.54569183
[433,] 46.34128888 49.44349035
[434,] 36.23908741 46.34128888
[435,] 30.13688594 36.23908741
[436,] 29.03468446 30.13688594
[437,] 28.93248299 29.03468446
[438,] 75.83028152 28.93248299
[439,] 81.72808004 75.83028152
[440,] 75.62587857 81.72808004
[441,] 43.52367710 75.62587857
[442,] 21.42147563 43.52367710
[443,] 12.31927415 21.42147563
[444,] 16.21707268 12.31927415
[445,] 4.11487121 16.21707268
[446,] -12.98733027 4.11487121
[447,] -19.08953174 -12.98733027
[448,] -34.19173321 -19.08953174
[449,] -46.29393468 -34.19173321
[450,] 9.60386384 -46.29393468
[451,] 19.50166237 9.60386384
[452,] -3.60053910 19.50166237
[453,] -18.70274058 -3.60053910
[454,] -35.80494205 -18.70274058
[455,] -31.90714352 -35.80494205
[456,] -29.00934499 -31.90714352
[457,] -38.11154647 -29.00934499
[458,] -53.21374794 -38.11154647
[459,] -56.31594941 -53.21374794
[460,] -77.41815089 -56.31594941
[461,] -68.52035236 -77.41815089
[462,] -18.62255383 -68.52035236
[463,] -12.72475530 -18.62255383
[464,] -28.82695678 -12.72475530
[465,] -40.92915825 -28.82695678
[466,] -45.03135972 -40.92915825
[467,] -31.13356120 -45.03135972
[468,] -19.23576267 -31.13356120
[469,] -14.33796414 -19.23576267
[470,] -11.44016561 -14.33796414
[471,] -10.54236709 -11.44016561
[472,] -23.64456856 -10.54236709
[473,] -11.74677003 -23.64456856
[474,] 39.15102849 -11.74677003
[475,] 49.04882702 39.15102849
[476,] 32.94662555 49.04882702
[477,] 15.84442408 32.94662555
[478,] 9.74222260 15.84442408
[479,] 26.64002113 9.74222260
[480,] 31.53781966 26.64002113
[481,] 26.43561818 31.53781966
[482,] 14.33341671 26.43561818
[483,] 3.23121524 14.33341671
[484,] -11.87098623 3.23121524
[485,] -3.97318771 -11.87098623
[486,] 51.92461082 -3.97318771
[487,] 54.82240935 51.92461082
[488,] 36.72020788 54.82240935
[489,] 14.61800640 36.72020788
[490,] -0.48419507 14.61800640
[491,] 1.41360346 -0.48419507
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -163.50767485 -150.40547338
2 -168.60987632 -163.50767485
3 -172.71207779 -168.60987632
4 -176.81427927 -172.71207779
5 -177.91648074 -176.81427927
6 -115.01868221 -177.91648074
7 -95.12088369 -115.01868221
8 -88.22308516 -95.12088369
9 -89.32528663 -88.22308516
10 -90.42748810 -89.32528663
11 -85.52968958 -90.42748810
12 -86.63189105 -85.52968958
13 -89.73409252 -86.63189105
14 -92.83629400 -89.73409252
15 -90.93849547 -92.83629400
16 -90.04069694 -90.93849547
17 -85.14289841 -90.04069694
18 -22.24509989 -85.14289841
19 -3.34730136 -22.24509989
20 0.55049717 -3.34730136
21 -3.55170431 0.55049717
22 -2.65390578 -3.55170431
23 4.24389275 -2.65390578
24 8.14169128 4.24389275
25 7.03948980 8.14169128
26 -0.06271167 7.03948980
27 -7.16491314 -0.06271167
28 -12.26711461 -7.16491314
29 -12.36931609 -12.26711461
30 50.52848244 -12.36931609
31 67.42628097 50.52848244
32 70.32407949 67.42628097
33 69.22187802 70.32407949
34 72.11967655 69.22187802
35 77.01747508 72.11967655
36 76.91527360 77.01747508
37 70.81307213 76.91527360
38 59.71087066 70.81307213
39 52.60866918 59.71087066
40 42.50646771 52.60866918
41 41.40426624 42.50646771
42 101.30206477 41.40426624
43 124.19986329 101.30206477
44 127.09766182 124.19986329
45 121.99546035 127.09766182
46 109.89325887 121.99546035
47 101.79105740 109.89325887
48 97.68885593 101.79105740
49 84.58665446 97.68885593
50 69.48445298 84.58665446
51 53.38225151 69.48445298
52 47.28005004 53.38225151
53 41.17784856 47.28005004
54 102.07564709 41.17784856
55 130.97344562 102.07564709
56 125.87124415 130.97344562
57 117.76904267 125.87124415
58 112.66684120 117.76904267
59 99.56463973 112.66684120
60 113.46243825 99.56463973
61 90.36023678 113.46243825
62 64.25803531 90.36023678
63 40.15583384 64.25803531
64 22.05363236 40.15583384
65 -0.04856911 22.05363236
66 48.84922942 -0.04856911
67 61.74702794 48.84922942
68 57.64482647 61.74702794
69 46.54262500 57.64482647
70 34.44042353 46.54262500
71 35.33822205 34.44042353
72 33.23602058 35.33822205
73 14.13381911 33.23602058
74 -1.96838237 14.13381911
75 -16.07058384 -1.96838237
76 -25.17278531 -16.07058384
77 -29.27498678 -25.17278531
78 15.62281174 -29.27498678
79 23.52061027 15.62281174
80 24.41840880 23.52061027
81 32.31620732 24.41840880
82 17.21400585 32.31620732
83 25.11180438 17.21400585
84 23.00960291 25.11180438
85 -0.09259857 23.00960291
86 -13.19480004 -0.09259857
87 -26.29700151 -13.19480004
88 -38.39920299 -26.29700151
89 -42.50140446 -38.39920299
90 6.39639407 -42.50140446
91 9.29419260 6.39639407
92 7.19199112 9.29419260
93 2.08978965 7.19199112
94 -9.01241182 2.08978965
95 -11.11461330 -9.01241182
96 -15.21681477 -11.11461330
97 -33.31901624 -15.21681477
98 -51.42121771 -33.31901624
99 -66.52341919 -51.42121771
100 -79.62562066 -66.52341919
101 -85.72782213 -79.62562066
102 -33.83002361 -85.72782213
103 -28.93222508 -33.83002361
104 -40.03442655 -28.93222508
105 -50.13662802 -40.03442655
106 -59.23882950 -50.13662802
107 -60.34103097 -59.23882950
108 -66.44323244 -60.34103097
109 -81.54543392 -66.44323244
110 -89.64763539 -81.54543392
111 -110.74983686 -89.64763539
112 -121.85203833 -110.74983686
113 -126.95423981 -121.85203833
114 -79.05644128 -126.95423981
115 -65.15864275 -79.05644128
116 -80.26084423 -65.15864275
117 -88.36304570 -80.26084423
118 -95.46524717 -88.36304570
119 -95.56744864 -95.46524717
120 -98.66965012 -95.56744864
121 -112.77185159 -98.66965012
122 -125.87405306 -112.77185159
123 -137.97625454 -125.87405306
124 -147.07845601 -137.97625454
125 -147.18065748 -147.07845601
126 -94.28285895 -147.18065748
127 -84.38506043 -94.28285895
128 -88.48726190 -84.38506043
129 -91.58946337 -88.48726190
130 -95.69166485 -91.58946337
131 -89.79386632 -95.69166485
132 -89.89606779 -89.79386632
133 -93.99826926 -89.89606779
134 -107.10047074 -93.99826926
135 -115.20267221 -107.10047074
136 -123.30487368 -115.20267221
137 -122.40707516 -123.30487368
138 -69.50927663 -122.40707516
139 -53.61147810 -69.50927663
140 -58.71367957 -53.61147810
141 -57.81588105 -58.71367957
142 -61.91808252 -57.81588105
143 -55.02028399 -61.91808252
144 -53.12248547 -55.02028399
145 -63.22468694 -53.12248547
146 -71.32688841 -63.22468694
147 -75.42908988 -71.32688841
148 -84.53129136 -75.42908988
149 -78.63349283 -84.53129136
150 -26.73569430 -78.63349283
151 -8.83789578 -26.73569430
152 -12.94009725 -8.83789578
153 -14.04229872 -12.94009725
154 -14.14450019 -14.04229872
155 -0.24670167 -14.14450019
156 5.65109686 -0.24670167
157 4.54889539 5.65109686
158 -3.55330609 4.54889539
159 -6.65550756 -3.55330609
160 -10.75770903 -6.65550756
161 -4.85991050 -10.75770903
162 54.03788802 -4.85991050
163 75.93568655 54.03788802
164 77.83348508 75.93568655
165 69.73128360 77.83348508
166 69.62908213 69.73128360
167 75.52688066 69.62908213
168 77.42467919 75.52688066
169 77.32247771 77.42467919
170 69.22027624 77.32247771
171 69.11807477 69.22027624
172 55.01587329 69.11807477
173 54.91367182 55.01587329
174 45.81147035 54.91367182
175 45.70926888 45.81147035
176 37.60706740 45.70926888
177 37.50486593 37.60706740
178 36.40266446 37.50486593
179 36.30046299 36.40266446
180 83.19826151 36.30046299
181 83.09606004 83.19826151
182 103.99385857 83.09606004
183 103.89165709 103.99385857
184 98.78945562 103.89165709
185 98.68725415 98.78945562
186 87.58505268 98.68725415
187 87.48285120 87.58505268
188 76.38064973 87.48285120
189 76.27844826 76.38064973
190 80.17624678 76.27844826
191 80.07404531 80.17624678
192 80.97184384 80.07404531
193 80.86964237 80.97184384
194 72.76744089 80.86964237
195 72.66523942 72.76744089
196 55.56303795 72.66523942
197 55.46083647 55.56303795
198 47.35863500 55.46083647
199 47.25643353 47.35863500
200 35.15423206 47.25643353
201 35.05203058 35.15423206
202 34.94982911 35.05203058
203 34.84762764 34.94982911
204 87.74542616 34.84762764
205 87.64322469 87.74542616
206 110.54102322 87.64322469
207 110.43882175 110.54102322
208 108.33662027 110.43882175
209 108.23441880 108.33662027
210 103.13221733 108.23441880
211 103.03001585 103.13221733
212 88.92781438 103.03001585
213 88.82561291 88.92781438
214 94.72341144 88.82561291
215 94.62120996 94.72341144
216 99.51900849 94.62120996
217 99.41680702 99.51900849
218 82.31460554 99.41680702
219 82.21240407 82.31460554
220 62.11020260 82.21240407
221 62.00800113 62.11020260
222 51.90579965 62.00800113
223 51.80359818 51.90579965
224 32.70139671 51.80359818
225 32.59919523 32.70139671
226 22.49699376 32.59919523
227 22.39479229 22.49699376
228 76.29259082 22.39479229
229 76.19038934 76.29259082
230 97.08818787 76.19038934
231 96.98598640 97.08818787
232 84.88378492 96.98598640
233 84.78158345 84.88378492
234 66.67938198 84.78158345
235 66.57718051 66.67938198
236 55.47497903 66.57718051
237 55.37277756 55.47497903
238 56.27057609 55.37277756
239 56.16837461 56.27057609
240 55.06617314 56.16837461
241 54.96397167 55.06617314
242 47.86177020 54.96397167
243 47.75956872 47.86177020
244 35.65736725 47.75956872
245 35.55516578 35.65736725
246 26.45296430 35.55516578
247 26.35076283 26.45296430
248 12.24856136 26.35076283
249 12.14635989 12.24856136
250 16.04415841 12.14635989
251 15.94195694 16.04415841
252 62.83975547 15.94195694
253 62.73755399 62.83975547
254 81.63535252 62.73755399
255 81.53315105 81.63535252
256 73.43094958 81.53315105
257 73.32874810 73.43094958
258 52.22654663 73.32874810
259 52.12434516 52.22654663
260 37.02214368 52.12434516
261 36.91994221 37.02214368
262 39.81774074 36.91994221
263 39.71553927 39.81774074
264 34.61333779 39.71553927
265 34.51113632 34.61333779
266 27.40893485 34.51113632
267 27.30673337 27.40893485
268 13.20453190 27.30673337
269 13.10233043 13.20453190
270 -0.99987104 13.10233043
271 -1.10207252 -0.99987104
272 -15.20427399 -1.10207252
273 -15.30647546 -15.20427399
274 -22.40867694 -15.30647546
275 -22.51087841 -22.40867694
276 26.38692012 -22.51087841
277 26.28471865 26.38692012
278 56.18251717 26.28471865
279 56.08031570 56.18251717
280 40.97811423 56.08031570
281 40.87591275 40.97811423
282 11.77371128 40.87591275
283 11.67150981 11.77371128
284 -6.43069166 11.67150981
285 -6.53289314 -6.43069166
286 -2.63509461 -6.53289314
287 -2.73729608 -2.63509461
288 -1.83949756 -2.73729608
289 -1.94169903 -1.83949756
290 -13.04390050 -1.94169903
291 -13.14610197 -13.04390050
292 -26.24830345 -13.14610197
293 -26.35050492 -26.24830345
294 -40.45270639 -26.35050492
295 -40.55490787 -40.45270639
296 -50.65710934 -40.55490787
297 -50.75931081 -50.65710934
298 -54.86151228 -50.75931081
299 -54.96371376 -54.86151228
300 -6.06591523 -54.96371376
301 -6.16811670 -6.06591523
302 21.72968182 -6.16811670
303 21.62748035 21.72968182
304 1.52527888 21.62748035
305 1.42307741 1.52527888
306 -19.67912407 1.42307741
307 -19.78132554 -19.67912407
308 -38.88352701 -19.78132554
309 -38.98572849 -38.88352701
310 -39.08792996 -38.98572849
311 -39.19013143 -39.08792996
312 -38.29233290 -39.19013143
313 -38.39453438 -38.29233290
314 -50.49673585 -38.39453438
315 -50.59893732 -50.49673585
316 -69.70113880 -50.59893732
317 -69.80334027 -69.70113880
318 -74.90554174 -69.80334027
319 -75.00774321 -74.90554174
320 -90.10994469 -75.00774321
321 -90.21214616 -90.10994469
322 -93.31434763 -90.21214616
323 -93.41654911 -93.31434763
324 -44.51875058 -93.41654911
325 -44.62095205 -44.51875058
326 -11.72315352 -44.62095205
327 -11.82535500 -11.72315352
328 -31.92755647 -11.82535500
329 -32.02975794 -31.92755647
330 -48.13195941 -32.02975794
331 -48.23416089 -48.13195941
332 -69.33636236 -48.23416089
333 -69.43856383 -69.33636236
334 -73.54076531 -69.43856383
335 -73.64296678 -73.54076531
336 -66.74516825 -73.64296678
337 -66.84736972 -66.74516825
338 -73.94957120 -66.84736972
339 -74.05177267 -73.94957120
340 -86.15397414 -74.05177267
341 -86.25617562 -86.15397414
342 -91.35837709 -86.25617562
343 -91.46057856 -91.35837709
344 -98.56278003 -91.46057856
345 -98.66498151 -98.56278003
346 -103.76718298 -98.66498151
347 -103.86938445 -103.76718298
348 -50.97158593 -103.86938445
349 -51.07378740 -50.97158593
350 -25.17598887 -51.07378740
351 -25.27819034 -25.17598887
352 -22.38039182 -25.27819034
353 -22.48259329 -22.38039182
354 -32.58479476 -22.48259329
355 -32.68699624 -32.58479476
356 -64.78919771 -32.68699624
357 -64.89139918 -64.78919771
358 -64.99360065 -64.89139918
359 -65.09580213 -64.99360065
360 -60.19800360 -65.09580213
361 -60.30020507 -60.19800360
362 -61.40240655 -60.30020507
363 -61.50460802 -61.40240655
364 -66.60680949 -61.50460802
365 -66.70901096 -66.60680949
366 -75.81121244 -66.70901096
367 -75.91341391 -75.81121244
368 -82.01561538 -75.91341391
369 -82.11781686 -82.01561538
370 -81.22001833 -82.11781686
371 -81.32221980 -81.22001833
372 -20.42442127 -81.32221980
373 -20.52662275 -20.42442127
374 -12.62882422 -20.52662275
375 -12.73102569 -12.62882422
376 -14.83322717 -12.73102569
377 -14.93542864 -14.83322717
378 -19.03763011 -14.93542864
379 -19.13983158 -19.03763011
380 -29.24203306 -19.13983158
381 -29.34423453 -29.24203306
382 -26.44643600 -29.34423453
383 -26.54863748 -26.44643600
384 -19.65083895 -26.54863748
385 -19.75304042 -19.65083895
386 -21.85524189 -19.75304042
387 -21.95744337 -21.85524189
388 -29.05964484 -21.95744337
389 -29.16184631 -29.05964484
390 -30.26404779 -29.16184631
391 -30.36624926 -30.26404779
392 -38.46845073 -30.36624926
393 -38.57065220 -38.46845073
394 -32.67285368 -38.57065220
395 -32.77505515 -32.67285368
396 29.12274338 -32.77505515
397 29.02054190 29.12274338
398 39.91834043 29.02054190
399 39.81613896 39.91834043
400 37.71393749 39.81613896
401 37.61173601 37.71393749
402 24.50953454 37.61173601
403 24.40733307 24.50953454
404 6.30513159 24.40733307
405 6.20293012 6.30513159
406 14.10072865 6.20293012
407 13.99852718 14.10072865
408 20.89632570 13.99852718
409 19.79412423 20.89632570
410 13.69192276 19.79412423
411 2.58972128 13.69192276
412 -4.51248019 2.58972128
413 1.38531834 -4.51248019
414 52.28311687 1.38531834
415 69.18091539 52.28311687
416 71.07871392 69.18091539
417 68.97651245 71.07871392
418 51.87431097 68.97651245
419 52.77210950 51.87431097
420 48.66990803 52.77210950
421 46.56770656 48.66990803
422 41.46550508 46.56770656
423 30.36330361 41.46550508
424 24.26110214 30.36330361
425 26.15890066 24.26110214
426 78.05669919 26.15890066
427 85.95449772 78.05669919
428 84.85229625 85.95449772
429 68.75009477 84.85229625
430 51.64789330 68.75009477
431 53.54569183 51.64789330
432 49.44349035 53.54569183
433 46.34128888 49.44349035
434 36.23908741 46.34128888
435 30.13688594 36.23908741
436 29.03468446 30.13688594
437 28.93248299 29.03468446
438 75.83028152 28.93248299
439 81.72808004 75.83028152
440 75.62587857 81.72808004
441 43.52367710 75.62587857
442 21.42147563 43.52367710
443 12.31927415 21.42147563
444 16.21707268 12.31927415
445 4.11487121 16.21707268
446 -12.98733027 4.11487121
447 -19.08953174 -12.98733027
448 -34.19173321 -19.08953174
449 -46.29393468 -34.19173321
450 9.60386384 -46.29393468
451 19.50166237 9.60386384
452 -3.60053910 19.50166237
453 -18.70274058 -3.60053910
454 -35.80494205 -18.70274058
455 -31.90714352 -35.80494205
456 -29.00934499 -31.90714352
457 -38.11154647 -29.00934499
458 -53.21374794 -38.11154647
459 -56.31594941 -53.21374794
460 -77.41815089 -56.31594941
461 -68.52035236 -77.41815089
462 -18.62255383 -68.52035236
463 -12.72475530 -18.62255383
464 -28.82695678 -12.72475530
465 -40.92915825 -28.82695678
466 -45.03135972 -40.92915825
467 -31.13356120 -45.03135972
468 -19.23576267 -31.13356120
469 -14.33796414 -19.23576267
470 -11.44016561 -14.33796414
471 -10.54236709 -11.44016561
472 -23.64456856 -10.54236709
473 -11.74677003 -23.64456856
474 39.15102849 -11.74677003
475 49.04882702 39.15102849
476 32.94662555 49.04882702
477 15.84442408 32.94662555
478 9.74222260 15.84442408
479 26.64002113 9.74222260
480 31.53781966 26.64002113
481 26.43561818 31.53781966
482 14.33341671 26.43561818
483 3.23121524 14.33341671
484 -11.87098623 3.23121524
485 -3.97318771 -11.87098623
486 51.92461082 -3.97318771
487 54.82240935 51.92461082
488 36.72020788 54.82240935
489 14.61800640 36.72020788
490 -0.48419507 14.61800640
491 1.41360346 -0.48419507
> 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/yougetitorg/rcomp/tmp/7h96i1296725581.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/yougetitorg/rcomp/tmp/8uuxf1296725581.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/yougetitorg/rcomp/tmp/9hygv1296725581.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/yougetitorg/rcomp/tmp/109eam1296725581.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/yougetitorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/yougetitorg/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/yougetitorg/rcomp/tmp/11ugg21296725581.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/yougetitorg/rcomp/tmp/12uqur1296725581.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/yougetitorg/rcomp/tmp/13tsxt1296725581.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/yougetitorg/rcomp/tmp/147dbd1296725581.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/yougetitorg/rcomp/tmp/153hvh1296725581.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/yougetitorg/rcomp/tmp/1659b81296725582.tab")
+ }
>
> try(system("convert tmp/1tl341296725581.ps tmp/1tl341296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/2kemt1296725581.ps tmp/2kemt1296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/3nrt31296725581.ps tmp/3nrt31296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/4hba91296725581.ps tmp/4hba91296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/5mzeo1296725581.ps tmp/5mzeo1296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/6q5jy1296725581.ps tmp/6q5jy1296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/7h96i1296725581.ps tmp/7h96i1296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/8uuxf1296725581.ps tmp/8uuxf1296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/9hygv1296725581.ps tmp/9hygv1296725581.png",intern=TRUE))
character(0)
> try(system("convert tmp/109eam1296725581.ps tmp/109eam1296725581.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
17.780 2.720 18.825