R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(48
+ ,3
+ ,0
+ ,38
+ ,4
+ ,0
+ ,31
+ ,4
+ ,5
+ ,31
+ ,4
+ ,4
+ ,106
+ ,5
+ ,8
+ ,81
+ ,5
+ ,12
+ ,21
+ ,6
+ ,15
+ ,35
+ ,6
+ ,12
+ ,24
+ ,6
+ ,2
+ ,121
+ ,7
+ ,13
+ ,143
+ ,7
+ ,6
+ ,56
+ ,7
+ ,8
+ ,52
+ ,8
+ ,5
+ ,52
+ ,8
+ ,15
+ ,73
+ ,9
+ ,11
+ ,109
+ ,10
+ ,5
+ ,56
+ ,10
+ ,17
+ ,44
+ ,10
+ ,11
+ ,43
+ ,10
+ ,12
+ ,23
+ ,10
+ ,16
+ ,285
+ ,11
+ ,12
+ ,162
+ ,11
+ ,2
+ ,68
+ ,11
+ ,8
+ ,92
+ ,11
+ ,12
+ ,62
+ ,11
+ ,11
+ ,42
+ ,11
+ ,9
+ ,183
+ ,12
+ ,12
+ ,115
+ ,12
+ ,8
+ ,114
+ ,12
+ ,6
+ ,73
+ ,12
+ ,20
+ ,502
+ ,13
+ ,9
+ ,151
+ ,13
+ ,13
+ ,113
+ ,13
+ ,10
+ ,101
+ ,13
+ ,15
+ ,67
+ ,13
+ ,22
+ ,49
+ ,13
+ ,14
+ ,89
+ ,14
+ ,5
+ ,72
+ ,14
+ ,22
+ ,235
+ ,15
+ ,13
+ ,228
+ ,15
+ ,25
+ ,254
+ ,15
+ ,19
+ ,85
+ ,15
+ ,30
+ ,38
+ ,15
+ ,4
+ ,401
+ ,16
+ ,8
+ ,174
+ ,16
+ ,10
+ ,70
+ ,16
+ ,17
+ ,181
+ ,16
+ ,13
+ ,173
+ ,16
+ ,15
+ ,145
+ ,16
+ ,19
+ ,123
+ ,16
+ ,18
+ ,372
+ ,17
+ ,29
+ ,150
+ ,17
+ ,16
+ ,121
+ ,17
+ ,17
+ ,114
+ ,17
+ ,6
+ ,94
+ ,17
+ ,7
+ ,176
+ ,18
+ ,61
+ ,214
+ ,18
+ ,11
+ ,201
+ ,18
+ ,21
+ ,91
+ ,18
+ ,18
+ ,47
+ ,18
+ ,6
+ ,212
+ ,19
+ ,17
+ ,536
+ ,19
+ ,18
+ ,112
+ ,19
+ ,21
+ ,110
+ ,19
+ ,21
+ ,183
+ ,20
+ ,24
+ ,233
+ ,20
+ ,19
+ ,157
+ ,20
+ ,22
+ ,144
+ ,20
+ ,12
+ ,37
+ ,20
+ ,35
+ ,109
+ ,21
+ ,18
+ ,107
+ ,21
+ ,18
+ ,94
+ ,21
+ ,22
+ ,256
+ ,22
+ ,14
+ ,156
+ ,22
+ ,8
+ ,149
+ ,22
+ ,19
+ ,129
+ ,22
+ ,12
+ ,74
+ ,23
+ ,11
+ ,153
+ ,23
+ ,24
+ ,139
+ ,25
+ ,29
+ ,137
+ ,25
+ ,24
+ ,122
+ ,25
+ ,24
+ ,236
+ ,26
+ ,31
+ ,232
+ ,26
+ ,22
+ ,161
+ ,27
+ ,41
+ ,195
+ ,27
+ ,20
+ ,145
+ ,27
+ ,28
+ ,191
+ ,28
+ ,28
+ ,131
+ ,28
+ ,53
+ ,482
+ ,29
+ ,22
+ ,186
+ ,29
+ ,34
+ ,165
+ ,29
+ ,11
+ ,147
+ ,29
+ ,13
+ ,145
+ ,29
+ ,35
+ ,118
+ ,30
+ ,33
+ ,169
+ ,30
+ ,22
+ ,109
+ ,30
+ ,24
+ ,98
+ ,30
+ ,28
+ ,215
+ ,31
+ ,30
+ ,170
+ ,32
+ ,17
+ ,171
+ ,32
+ ,20
+ ,132
+ ,32
+ ,27
+ ,66
+ ,32
+ ,16
+ ,58
+ ,32
+ ,41
+ ,198
+ ,33
+ ,9
+ ,188
+ ,33
+ ,21
+ ,146
+ ,33
+ ,20
+ ,99
+ ,33
+ ,20
+ ,255
+ ,34
+ ,29
+ ,328
+ ,34
+ ,36
+ ,266
+ ,34
+ ,69
+ ,201
+ ,34
+ ,22
+ ,183
+ ,34
+ ,25
+ ,169
+ ,34
+ ,15
+ ,278
+ ,35
+ ,21
+ ,35
+ ,33
+ ,97
+ ,35
+ ,49
+ ,335
+ ,36
+ ,40
+ ,224
+ ,36
+ ,113
+ ,178
+ ,36
+ ,36
+ ,153
+ ,36
+ ,31
+ ,121
+ ,37
+ ,34
+ ,429
+ ,37
+ ,35
+ ,103
+ ,37
+ ,27
+ ,179
+ ,38
+ ,65
+ ,187
+ ,38
+ ,19
+ ,162
+ ,38
+ ,22
+ ,225
+ ,39
+ ,69
+ ,31
+ ,39
+ ,48
+ ,282
+ ,41
+ ,42
+ ,279
+ ,41
+ ,58
+ ,221
+ ,41
+ ,33
+ ,217
+ ,41
+ ,31
+ ,199
+ ,41
+ ,34
+ ,415
+ ,42
+ ,56
+ ,114
+ ,42
+ ,25
+ ,157
+ ,42
+ ,37
+ ,355
+ ,43
+ ,35
+ ,265
+ ,43
+ ,51
+ ,304
+ ,43
+ ,28
+ ,165
+ ,43
+ ,34
+ ,239
+ ,44
+ ,72
+ ,159
+ ,45
+ ,27
+ ,93
+ ,45
+ ,18
+ ,170
+ ,46
+ ,17
+ ,144
+ ,46
+ ,100
+ ,369
+ ,47
+ ,31
+ ,148
+ ,47
+ ,53
+ ,209
+ ,47
+ ,36
+ ,179
+ ,48
+ ,49
+ ,125
+ ,48
+ ,64
+ ,98
+ ,48
+ ,109
+ ,140
+ ,48
+ ,43
+ ,124
+ ,48
+ ,56
+ ,185
+ ,49
+ ,45
+ ,125
+ ,49
+ ,40
+ ,344
+ ,50
+ ,54
+ ,222
+ ,50
+ ,41
+ ,150
+ ,50
+ ,46
+ ,146
+ ,50
+ ,45
+ ,215
+ ,50
+ ,37
+ ,99
+ ,50
+ ,72
+ ,448
+ ,51
+ ,47
+ ,175
+ ,51
+ ,140
+ ,175
+ ,51
+ ,63
+ ,345
+ ,52
+ ,60
+ ,449
+ ,52
+ ,62
+ ,156
+ ,53
+ ,34
+ ,144
+ ,53
+ ,39
+ ,474
+ ,54
+ ,60
+ ,189
+ ,54
+ ,43
+ ,136
+ ,54
+ ,40
+ ,287
+ ,54
+ ,43
+ ,177
+ ,54
+ ,129
+ ,156
+ ,55
+ ,46
+ ,129
+ ,55
+ ,42
+ ,472
+ ,58
+ ,53
+ ,232
+ ,58
+ ,40
+ ,141
+ ,58
+ ,96
+ ,306
+ ,59
+ ,36
+ ,144
+ ,59
+ ,58
+ ,336
+ ,59
+ ,59
+ ,181
+ ,59
+ ,48
+ ,121
+ ,59
+ ,74
+ ,173
+ ,60
+ ,47
+ ,168
+ ,61
+ ,94
+ ,398
+ ,62
+ ,65
+ ,392
+ ,62
+ ,62
+ ,434
+ ,62
+ ,32
+ ,286
+ ,62
+ ,31
+ ,153
+ ,62
+ ,48
+ ,48
+ ,62
+ ,64
+ ,41
+ ,62
+ ,16
+ ,194
+ ,63
+ ,51
+ ,324
+ ,65
+ ,64
+ ,455
+ ,65
+ ,100
+ ,295
+ ,66
+ ,86
+ ,188
+ ,66
+ ,113
+ ,95
+ ,66
+ ,38
+ ,1
+ ,182
+ ,68
+ ,62
+ ,264
+ ,68
+ ,48
+ ,492
+ ,69
+ ,46
+ ,423
+ ,70
+ ,43
+ ,201
+ ,70
+ ,68
+ ,186
+ ,70
+ ,72
+ ,356
+ ,71
+ ,73
+ ,326
+ ,72
+ ,82
+ ,192
+ ,73
+ ,82
+ ,511
+ ,73
+ ,81
+ ,458
+ ,74
+ ,62
+ ,359
+ ,74
+ ,54
+ ,378
+ ,74
+ ,77
+ ,232
+ ,75
+ ,63
+ ,483
+ ,75
+ ,39
+ ,330
+ ,75
+ ,201
+ ,324
+ ,75
+ ,58
+ ,204
+ ,75
+ ,96
+ ,163
+ ,77
+ ,66
+ ,236
+ ,78
+ ,113
+ ,542
+ ,79
+ ,85
+ ,326
+ ,79
+ ,58
+ ,347
+ ,80
+ ,51
+ ,456
+ ,80
+ ,152
+ ,380
+ ,81
+ ,68
+ ,150
+ ,81
+ ,159
+ ,185
+ ,83
+ ,49
+ ,434
+ ,83
+ ,65
+ ,284
+ ,83
+ ,45
+ ,574
+ ,84
+ ,55
+ ,406
+ ,84
+ ,63
+ ,352
+ ,85
+ ,80
+ ,305
+ ,85
+ ,113
+ ,208
+ ,85
+ ,97
+ ,343
+ ,86
+ ,202
+ ,204
+ ,86
+ ,87
+ ,145
+ ,88
+ ,45
+ ,265
+ ,89
+ ,61
+ ,539
+ ,89
+ ,83
+ ,291
+ ,89
+ ,45
+ ,343
+ ,90
+ ,71
+ ,163
+ ,90
+ ,73
+ ,100
+ ,90
+ ,48
+ ,883
+ ,91
+ ,80
+ ,452
+ ,91
+ ,117
+ ,280
+ ,91
+ ,168
+ ,205
+ ,92
+ ,127
+ ,433
+ ,93
+ ,91
+ ,372
+ ,94
+ ,59
+ ,202
+ ,94
+ ,184
+ ,321
+ ,96
+ ,49
+ ,240
+ ,96
+ ,114
+ ,243
+ ,97
+ ,119
+ ,436
+ ,99
+ ,82
+ ,239
+ ,99
+ ,85
+ ,589
+ ,100
+ ,71
+ ,238
+ ,102
+ ,79
+ ,956
+ ,103
+ ,136
+ ,657
+ ,105
+ ,94
+ ,372
+ ,105
+ ,235
+ ,222
+ ,105
+ ,185
+ ,397
+ ,106
+ ,114
+ ,1
+ ,003
+ ,107
+ ,121
+ ,509
+ ,108
+ ,77
+ ,388
+ ,108
+ ,171
+ ,490
+ ,109
+ ,67
+ ,259
+ ,109
+ ,77
+ ,355
+ ,110
+ ,147
+ ,253
+ ,111
+ ,52
+ ,488
+ ,113
+ ,99
+ ,160
+ ,115
+ ,166
+ ,307
+ ,115
+ ,288
+ ,554
+ ,116
+ ,143
+ ,505
+ ,118
+ ,180
+ ,1
+ ,083
+ ,120
+ ,76
+ ,426
+ ,120
+ ,111
+ ,1
+ ,000
+ ,123
+ ,85
+ ,323
+ ,126
+ ,109
+ ,1
+ ,061
+ ,127
+ ,95
+ ,259
+ ,131
+ ,98
+ ,737
+ ,132
+ ,167
+ ,559
+ ,132
+ ,101
+ ,596
+ ,138
+ ,76
+ ,537
+ ,140
+ ,95
+ ,719
+ ,141
+ ,111
+ ,489
+ ,141
+ ,95
+ ,478
+ ,141
+ ,106
+ ,337
+ ,142
+ ,74
+ ,856
+ ,145
+ ,145
+ ,449
+ ,146
+ ,77
+ ,359
+ ,148
+ ,158
+ ,469
+ ,148
+ ,140
+ ,434
+ ,150
+ ,151
+ ,357
+ ,150
+ ,118
+ ,768
+ ,153
+ ,111
+ ,373
+ ,155
+ ,475
+ ,581
+ ,156
+ ,236
+ ,587
+ ,157
+ ,159
+ ,526
+ ,158
+ ,195
+ ,195
+ ,158
+ ,58
+ ,614
+ ,160
+ ,99
+ ,553
+ ,160
+ ,135
+ ,558
+ ,160
+ ,93
+ ,417
+ ,160
+ ,311
+ ,921
+ ,162
+ ,158
+ ,476
+ ,163
+ ,176
+ ,787
+ ,164
+ ,131
+ ,908
+ ,166
+ ,100
+ ,552
+ ,168
+ ,125
+ ,850
+ ,170
+ ,150
+ ,582
+ ,173
+ ,192
+ ,514
+ ,174
+ ,147
+ ,386
+ ,174
+ ,152
+ ,530
+ ,176
+ ,171
+ ,546
+ ,180
+ ,201
+ ,983
+ ,181
+ ,148
+ ,495
+ ,181
+ ,203
+ ,1
+ ,057
+ ,182
+ ,172
+ ,548
+ ,182
+ ,274
+ ,522
+ ,190
+ ,419
+ ,791
+ ,194
+ ,176
+ ,994
+ ,196
+ ,174
+ ,569
+ ,196
+ ,219
+ ,764
+ ,197
+ ,119
+ ,602
+ ,198
+ ,192
+ ,831
+ ,199
+ ,205
+ ,568
+ ,200
+ ,111
+ ,948
+ ,206
+ ,175
+ ,351
+ ,206
+ ,165
+ ,367
+ ,211
+ ,266
+ ,1
+ ,001
+ ,215
+ ,200
+ ,708
+ ,215
+ ,303
+ ,644
+ ,224
+ ,412
+ ,467
+ ,225
+ ,191
+ ,991
+ ,229
+ ,127
+ ,586
+ ,232
+ ,288
+ ,992
+ ,236
+ ,141
+ ,717
+ ,240
+ ,189
+ ,1
+ ,139
+ ,251
+ ,222
+ ,563
+ ,251
+ ,263
+ ,806
+ ,255
+ ,300
+ ,549
+ ,256
+ ,419
+ ,1
+ ,157
+ ,262
+ ,346
+ ,528
+ ,262
+ ,369
+ ,1
+ ,533
+ ,267
+ ,235
+ ,756
+ ,269
+ ,256
+ ,1
+ ,039
+ ,271
+ ,430
+ ,1
+ ,070
+ ,275
+ ,172
+ ,1
+ ,053
+ ,278
+ ,351
+ ,784
+ ,278
+ ,311
+ ,756
+ ,283
+ ,208
+ ,1
+ ,003
+ ,285
+ ,243
+ ,1
+ ,155
+ ,298
+ ,292
+ ,2
+ ,548
+ ,301
+ ,277
+ ,1
+ ,057
+ ,303
+ ,259
+ ,850
+ ,309
+ ,245
+ ,632
+ ,314
+ ,357
+ ,501
+ ,325
+ ,415
+ ,1
+ ,528
+ ,338
+ ,250
+ ,1
+ ,776
+ ,338
+ ,343
+ ,921
+ ,346
+ ,505
+ ,414
+ ,356
+ ,342
+ ,1
+ ,219
+ ,360
+ ,107
+ ,1
+ ,738
+ ,368
+ ,273
+ ,1
+ ,084
+ ,375
+ ,306
+ ,767
+ ,386
+ ,725
+ ,2
+ ,372
+ ,388
+ ,222
+ ,2
+ ,036
+ ,391
+ ,311
+ ,1
+ ,729
+ ,393
+ ,259
+ ,1
+ ,393
+ ,402
+ ,250
+ ,847
+ ,404
+ ,552
+ ,222
+ ,407
+ ,99
+ ,1
+ ,076
+ ,408
+ ,565
+ ,1
+ ,062
+ ,409
+ ,338
+ ,416
+ ,410
+ ,160
+ ,343
+ ,417
+ ,109
+ ,1
+ ,813
+ ,424
+ ,312
+ ,1
+ ,920
+ ,435
+ ,318
+ ,2
+ ,392
+ ,436
+ ,377
+ ,156
+ ,443
+ ,221
+ ,2
+ ,230
+ ,449
+ ,547
+ ,564
+ ,451
+ ,566
+ ,1
+ ,549
+ ,456
+ ,245
+ ,5
+ ,091
+ ,457
+ ,300
+ ,1
+ ,370
+ ,462
+ ,403
+ ,622
+ ,469
+ ,195
+ ,501
+ ,479
+ ,220
+ ,1
+ ,187
+ ,480
+ ,753
+ ,1
+ ,217
+ ,481
+ ,868
+ ,1
+ ,278
+ ,487
+ ,388
+ ,2
+ ,251
+ ,497
+ ,399
+ ,1
+ ,156
+ ,520
+ ,478
+ ,1
+ ,357
+ ,525
+ ,441
+ ,2
+ ,116
+ ,534
+ ,434
+ ,2
+ ,220
+ ,536
+ ,451
+ ,509
+ ,537
+ ,694
+ ,1
+ ,080
+ ,548
+ ,991
+ ,2
+ ,748
+ ,609
+ ,489
+ ,2
+ ,099
+ ,626
+ ,612
+ ,2
+ ,477
+ ,630
+ ,776
+ ,1
+ ,456
+ ,635
+ ,303
+ ,3
+ ,974
+ ,649
+ ,578
+ ,3
+ ,826
+ ,666
+ ,736
+ ,2
+ ,572
+ ,670
+ ,637
+ ,2
+ ,205
+ ,673
+ ,724
+ ,782
+ ,681
+ ,1
+ ,120
+ ,2
+ ,848
+ ,695
+ ,1
+ ,559
+ ,2
+ ,649
+ ,696
+ ,773
+ ,2
+ ,064
+ ,740
+ ,859
+ ,1
+ ,752
+ ,744
+ ,710
+ ,2
+ ,967
+ ,754
+ ,671
+ ,2
+ ,625
+ ,756
+ ,629
+ ,1
+ ,256
+ ,756
+ ,530
+ ,2
+ ,442
+ ,798
+ ,1
+ ,651
+ ,1
+ ,684
+ ,804
+ ,625
+ ,5
+ ,446
+ ,834
+ ,778
+ ,2
+ ,193
+ ,837
+ ,853
+ ,2
+ ,102
+ ,845
+ ,730
+ ,2
+ ,346
+ ,864
+ ,1
+ ,068
+ ,2
+ ,096
+ ,871
+ ,558
+ ,3
+ ,126
+ ,874
+ ,687
+ ,3
+ ,193
+ ,877
+ ,778
+ ,2
+ ,225
+ ,919
+ ,668
+ ,1
+ ,685
+ ,938
+ ,1
+ ,274
+ ,4
+ ,693
+ ,952
+ ,976
+ ,821
+ ,965
+ ,1
+ ,715
+ ,4
+ ,653
+ ,1
+ ,032
+ ,798
+ ,1
+ ,292
+ ,1
+ ,069
+ ,1
+ ,507
+ ,5
+ ,160
+ ,1
+ ,100
+ ,1
+ ,038
+ ,3
+ ,104
+ ,1
+ ,133
+ ,1
+ ,270
+ ,5
+ ,057
+ ,1
+ ,201
+ ,988
+ ,5
+ ,644
+ ,1
+ ,223
+ ,886
+ ,5
+ ,695
+ ,1
+ ,249
+ ,1
+ ,122
+ ,4
+ ,586
+ ,1
+ ,279
+ ,869
+ ,3
+ ,819
+ ,1
+ ,313
+ ,1
+ ,690
+ ,3
+ ,556
+ ,1
+ ,353
+ ,1
+ ,201
+ ,4
+ ,398
+ ,1
+ ,370
+ ,1
+ ,260
+ ,5
+ ,799
+ ,1
+ ,380
+ ,1
+ ,166
+ ,4
+ ,735
+ ,1
+ ,393
+ ,1
+ ,336
+ ,3
+ ,372
+ ,1
+ ,500
+ ,813
+ ,4
+ ,880
+ ,1
+ ,508
+ ,1
+ ,522
+ ,6
+ ,216
+ ,1
+ ,578
+ ,1
+ ,199
+ ,8
+ ,041
+ ,1
+ ,609
+ ,1
+ ,416
+ ,6
+ ,884
+ ,1
+ ,617
+ ,1
+ ,909
+ ,8
+ ,556
+ ,1
+ ,779
+ ,1
+ ,522
+ ,3
+ ,858
+ ,1
+ ,900
+ ,2
+ ,359
+ ,9
+ ,309
+ ,2
+ ,001
+ ,1
+ ,689
+ ,7
+ ,639
+ ,2
+ ,191
+ ,1
+ ,594
+ ,5
+ ,039
+ ,2
+ ,192
+ ,2
+ ,008
+ ,6
+ ,321
+ ,2
+ ,198
+ ,1
+ ,262
+ ,5
+ ,363
+ ,2
+ ,283
+ ,1
+ ,890
+ ,6
+ ,642
+ ,2
+ ,554
+ ,2
+ ,599
+ ,5
+ ,338
+ ,3
+ ,028
+ ,3
+ ,107
+ ,5
+ ,865
+ ,3
+ ,446
+ ,4
+ ,404
+ ,5
+ ,395
+ ,3
+ ,754
+ ,4
+ ,436
+ ,11
+ ,218
+ ,4
+ ,293
+ ,4
+ ,275
+ ,10
+ ,796
+ ,4
+ ,656
+ ,4
+ ,944
+ ,9
+ ,594
+ ,5
+ ,145
+ ,5
+ ,441
+ ,22
+ ,162
+ ,6
+ ,225
+ ,8
+ ,019
+ ,17
+ ,773
+ ,6
+ ,618
+ ,6
+ ,634
+ ,15
+ ,153
+ ,7
+ ,590
+ ,7
+ ,508
+ ,20
+ ,202
+ ,7
+ ,663
+ ,6
+ ,509
+ ,17
+ ,094
+ ,9
+ ,945
+ ,11
+ ,166
+ ,28
+ ,435
+ ,10
+ ,277
+ ,7
+ ,551
+ ,41
+ ,169
+ ,14
+ ,738
+ ,13
+ ,704
+ ,44
+ ,525
+ ,15
+ ,100
+ ,15
+ ,136
+ ,45
+ ,104
+ ,18
+ ,103
+ ,14
+ ,819
+ ,31
+ ,118
+ ,22
+ ,259
+ ,19
+ ,638
+ ,72
+ ,772
+ ,26
+ ,073
+ ,22
+ ,274)
+ ,dim=c(3
+ ,467)
+ ,dimnames=list(c('weekdag'
+ ,'zaterdag'
+ ,'zondag')
+ ,1:467))
> y <- array(NA,dim=c(3,467),dimnames=list(c('weekdag','zaterdag','zondag'),1:467))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '3'
> #'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
zondag weekdag zaterdag
1 0 48 3
2 0 38 4
3 5 31 4
4 4 31 4
5 8 106 5
6 12 81 5
7 15 21 6
8 12 35 6
9 2 24 6
10 13 121 7
11 6 143 7
12 8 56 7
13 5 52 8
14 15 52 8
15 11 73 9
16 5 109 10
17 17 56 10
18 11 44 10
19 12 43 10
20 16 23 10
21 12 285 11
22 2 162 11
23 8 68 11
24 12 92 11
25 11 62 11
26 9 42 11
27 12 183 12
28 8 115 12
29 6 114 12
30 20 73 12
31 9 502 13
32 13 151 13
33 10 113 13
34 15 101 13
35 22 67 13
36 14 49 13
37 5 89 14
38 22 72 14
39 13 235 15
40 25 228 15
41 19 254 15
42 30 85 15
43 4 38 15
44 8 401 16
45 10 174 16
46 17 70 16
47 13 181 16
48 15 173 16
49 19 145 16
50 18 123 16
51 29 372 17
52 16 150 17
53 17 121 17
54 6 114 17
55 7 94 17
56 61 176 18
57 11 214 18
58 21 201 18
59 18 91 18
60 6 47 18
61 17 212 19
62 18 536 19
63 21 112 19
64 21 110 19
65 24 183 20
66 19 233 20
67 22 157 20
68 12 144 20
69 35 37 20
70 18 109 21
71 18 107 21
72 22 94 21
73 14 256 22
74 8 156 22
75 19 149 22
76 12 129 22
77 11 74 23
78 24 153 23
79 29 139 25
80 24 137 25
81 24 122 25
82 31 236 26
83 22 232 26
84 41 161 27
85 20 195 27
86 28 145 27
87 28 191 28
88 53 131 28
89 22 482 29
90 34 186 29
91 11 165 29
92 13 147 29
93 35 145 29
94 33 118 30
95 22 169 30
96 24 109 30
97 28 98 30
98 30 215 31
99 17 170 32
100 20 171 32
101 27 132 32
102 16 66 32
103 41 58 32
104 9 198 33
105 21 188 33
106 20 146 33
107 20 99 33
108 29 255 34
109 36 328 34
110 69 266 34
111 22 201 34
112 25 183 34
113 15 169 34
114 21 278 35
115 97 35 33
116 335 35 49
117 224 36 40
118 178 36 113
119 153 36 36
120 121 36 31
121 429 37 34
122 103 37 35
123 179 37 27
124 187 38 65
125 162 38 19
126 225 38 22
127 31 39 69
128 282 39 48
129 279 41 42
130 221 41 58
131 217 41 33
132 199 41 31
133 415 41 34
134 114 42 56
135 157 42 25
136 355 42 37
137 265 43 35
138 304 43 51
139 165 43 28
140 239 43 34
141 159 44 72
142 93 45 27
143 170 45 18
144 144 46 17
145 369 46 100
146 148 47 31
147 209 47 53
148 179 47 36
149 125 48 49
150 98 48 64
151 140 48 109
152 124 48 43
153 185 48 56
154 125 49 45
155 344 49 40
156 222 50 54
157 150 50 41
158 146 50 46
159 215 50 45
160 99 50 37
161 448 50 72
162 175 51 47
163 175 51 140
164 345 51 63
165 449 52 60
166 156 52 62
167 144 53 34
168 474 53 39
169 189 54 60
170 136 54 43
171 287 54 40
172 177 54 43
173 156 54 129
174 129 55 46
175 472 55 42
176 232 58 53
177 141 58 40
178 306 58 96
179 144 59 36
180 336 59 58
181 181 59 59
182 121 59 48
183 173 59 74
184 168 60 47
185 398 61 94
186 392 62 65
187 434 62 62
188 286 62 32
189 153 62 31
190 48 62 48
191 41 62 64
192 194 62 16
193 324 63 51
194 455 65 64
195 295 65 100
196 188 66 86
197 95 66 113
198 1 66 38
199 62 182 68
200 48 264 68
201 46 492 69
202 43 423 70
203 68 201 70
204 72 186 70
205 73 356 71
206 82 326 72
207 82 192 73
208 81 511 73
209 62 458 74
210 54 359 74
211 77 378 74
212 63 232 75
213 39 483 75
214 201 330 75
215 58 324 75
216 96 204 75
217 66 163 77
218 113 236 78
219 85 542 79
220 58 326 79
221 51 347 80
222 152 456 80
223 68 380 81
224 159 150 81
225 49 185 83
226 65 434 83
227 45 284 83
228 55 574 84
229 63 406 84
230 80 352 85
231 113 305 85
232 97 208 85
233 202 343 86
234 87 204 86
235 45 145 88
236 61 265 89
237 83 539 89
238 45 291 89
239 71 343 90
240 73 163 90
241 48 100 90
242 80 883 91
243 117 452 91
244 168 280 91
245 127 205 92
246 91 433 93
247 59 372 94
248 184 202 94
249 49 321 96
250 114 240 96
251 119 243 97
252 82 436 99
253 85 239 99
254 71 589 100
255 79 238 102
256 136 956 103
257 94 657 105
258 235 372 105
259 185 222 105
260 114 397 106
261 107 1 3
262 108 121 509
263 108 77 388
264 109 171 490
265 109 67 259
266 110 77 355
267 111 147 253
268 113 52 488
269 115 99 160
270 115 166 307
271 116 288 554
272 118 143 505
273 83 180 1
274 426 120 76
275 1 120 111
276 85 0 123
277 109 323 126
278 127 1 61
279 131 95 259
280 132 98 737
281 132 167 559
282 138 101 596
283 140 76 537
284 141 95 719
285 141 111 489
286 141 95 478
287 142 106 337
288 145 74 856
289 146 145 449
290 148 77 359
291 148 158 469
292 150 140 434
293 150 151 357
294 153 118 768
295 155 111 373
296 156 475 581
297 157 236 587
298 158 159 526
299 158 195 195
300 160 58 614
301 160 99 553
302 160 135 558
303 160 93 417
304 162 311 921
305 163 158 476
306 164 176 787
307 166 131 908
308 168 100 552
309 170 125 850
310 173 150 582
311 174 192 514
312 174 147 386
313 176 152 530
314 180 171 546
315 181 201 983
316 181 148 495
317 57 203 1
318 548 182 172
319 522 182 274
320 791 190 419
321 994 194 176
322 569 196 174
323 764 196 219
324 602 197 119
325 831 198 192
326 568 199 205
327 948 200 111
328 351 206 175
329 367 206 165
330 1 211 266
331 200 1 215
332 303 708 215
333 412 644 224
334 191 467 225
335 127 991 229
336 288 586 232
337 141 992 236
338 189 717 240
339 251 1 139
340 251 222 563
341 255 263 806
342 256 300 549
343 157 419 1
344 528 262 346
345 1 262 369
346 235 533 267
347 256 756 269
348 271 1 39
349 70 430 1
350 1 275 172
351 351 53 278
352 311 784 278
353 208 756 283
354 285 1 3
355 155 243 1
356 2 298 292
357 277 548 301
358 303 1 57
359 309 259 850
360 314 245 632
361 325 357 501
362 528 415 1
363 1 338 250
364 343 776 338
365 505 921 346
366 342 414 356
367 360 1 219
368 738 107 1
369 1 368 273
370 306 84 375
371 725 767 386
372 388 2 372
373 36 222 2
374 1 391 311
375 259 729 393
376 402 1 393
377 404 250 847
378 407 552 222
379 76 99 1
380 1 408 565
381 338 62 409
382 160 416 410
383 109 343 417
384 424 1 813
385 920 312 1
386 2 435 318
387 377 392 436
388 221 156 443
389 449 2 230
390 451 547 564
391 549 566 1
392 5 456 245
393 300 91 457
394 462 1 370
395 469 403 622
396 479 195 501
397 187 220 1
398 1 480 753
399 868 217 481
400 487 1 278
401 251 388 2
402 1 497 399
403 478 156 520
404 525 1 357
405 116 441 2
406 2 534 434
407 451 220 536
408 694 509 537
409 548 1 80
410 748 991 2
411 2 609 489
412 612 99 626
413 630 2 477
414 456 776 1
415 3 635 303
416 578 974 649
417 666 3 826
418 572 736 2
419 2 670 637
420 724 205 673
421 1 782 681
422 848 120 2
423 559 695 1
424 696 2 649
425 64 773 2
426 1 740 859
427 710 752 744
428 754 2 967
429 625 671 2
430 1 756 629
431 530 256 756
432 798 2 442
433 1 1 651
434 625 684 804
435 834 5 446
436 193 778 2
437 2 837 853
438 730 102 845
439 864 2 346
440 2 1 68
441 558 96 871
442 874 3 126
443 193 687 3
444 2 877 778
445 668 225 919
446 938 1 685
447 4 1 274
448 976 693 952
449 1 821 965
450 653 715 4
451 798 1 32
452 1 1 292
453 507 69 1
454 1 5 160
455 38 100 1
456 1 3 104
457 270 133 1
458 1 5 57
459 5 201 988
460 223 644 1
461 695 886 5
462 1 1 249
463 586 122 4
464 869 1 279
465 1 3 819
466 690 313 1
467 1 3 556
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) weekdag zaterdag
137.42780 0.01666 0.25350
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-394.73 -127.48 -84.38 69.12 808.72
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 137.42780 14.92374 9.209 < 2e-16 ***
weekdag 0.01666 0.04535 0.367 0.714
zaterdag 0.25350 0.04045 6.267 8.39e-10 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 210.6 on 464 degrees of freedom
Multiple R-squared: 0.08214, Adjusted R-squared: 0.07819
F-statistic: 20.76 on 2 and 464 DF, p-value: 2.311e-09
> 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,] 3.982276e-06 7.964552e-06 0.9999960
[2,] 4.778479e-08 9.556958e-08 1.0000000
[3,] 7.399959e-10 1.479992e-09 1.0000000
[4,] 4.253715e-10 8.507431e-10 1.0000000
[5,] 1.020401e-11 2.040801e-11 1.0000000
[6,] 6.310355e-13 1.262071e-12 1.0000000
[7,] 1.699596e-14 3.399192e-14 1.0000000
[8,] 9.438103e-16 1.887621e-15 1.0000000
[9,] 2.935511e-17 5.871022e-17 1.0000000
[10,] 6.774302e-19 1.354860e-18 1.0000000
[11,] 5.620200e-20 1.124040e-19 1.0000000
[12,] 1.965112e-21 3.930224e-21 1.0000000
[13,] 4.663949e-23 9.327898e-23 1.0000000
[14,] 1.002978e-24 2.005955e-24 1.0000000
[15,] 2.586975e-26 5.173949e-26 1.0000000
[16,] 5.291270e-28 1.058254e-27 1.0000000
[17,] 1.222578e-28 2.445157e-28 1.0000000
[18,] 4.235753e-30 8.471505e-30 1.0000000
[19,] 9.342275e-32 1.868455e-31 1.0000000
[20,] 2.085603e-33 4.171207e-33 1.0000000
[21,] 5.669721e-35 1.133944e-34 1.0000000
[22,] 1.208281e-36 2.416561e-36 1.0000000
[23,] 3.603686e-38 7.207372e-38 1.0000000
[24,] 1.530161e-39 3.060322e-39 1.0000000
[25,] 1.460769e-40 2.921537e-40 1.0000000
[26,] 3.312572e-42 6.625145e-42 1.0000000
[27,] 7.064492e-44 1.412898e-43 1.0000000
[28,] 1.767976e-45 3.535953e-45 1.0000000
[29,] 4.108193e-47 8.216386e-47 1.0000000
[30,] 4.550331e-48 9.100662e-48 1.0000000
[31,] 9.658946e-50 1.931789e-49 1.0000000
[32,] 1.605556e-50 3.211111e-50 1.0000000
[33,] 1.236476e-51 2.472952e-51 1.0000000
[34,] 2.767322e-53 5.534645e-53 1.0000000
[35,] 7.588692e-54 1.517738e-53 1.0000000
[36,] 2.592139e-55 5.184277e-55 1.0000000
[37,] 1.279305e-55 2.558609e-55 1.0000000
[38,] 7.096855e-56 1.419371e-55 1.0000000
[39,] 3.373626e-57 6.747253e-57 1.0000000
[40,] 1.576449e-58 3.152898e-58 1.0000000
[41,] 4.190748e-60 8.381496e-60 1.0000000
[42,] 1.230198e-61 2.460396e-61 1.0000000
[43,] 3.227431e-63 6.454862e-63 1.0000000
[44,] 9.616756e-65 1.923351e-64 1.0000000
[45,] 2.561540e-66 5.123079e-66 1.0000000
[46,] 1.276241e-66 2.552481e-66 1.0000000
[47,] 3.578522e-68 7.157044e-68 1.0000000
[48,] 9.586279e-70 1.917256e-69 1.0000000
[49,] 1.913690e-70 3.827379e-70 1.0000000
[50,] 2.119404e-71 4.238808e-71 1.0000000
[51,] 5.278694e-65 1.055739e-64 1.0000000
[52,] 3.456291e-66 6.912581e-66 1.0000000
[53,] 1.362418e-67 2.724836e-67 1.0000000
[54,] 5.329716e-69 1.065943e-68 1.0000000
[55,] 7.400382e-70 1.480076e-69 1.0000000
[56,] 2.997659e-71 5.995318e-71 1.0000000
[57,] 1.167120e-72 2.334240e-72 1.0000000
[58,] 4.514293e-74 9.028587e-74 1.0000000
[59,] 1.729503e-75 3.459007e-75 1.0000000
[60,] 7.072090e-77 1.414418e-76 1.0000000
[61,] 2.687407e-78 5.374815e-78 1.0000000
[62,] 1.005501e-79 2.011001e-79 1.0000000
[63,] 6.719980e-81 1.343996e-80 1.0000000
[64,] 1.221366e-81 2.442732e-81 1.0000000
[65,] 5.264456e-83 1.052891e-82 1.0000000
[66,] 2.223106e-84 4.446212e-84 1.0000000
[67,] 8.235825e-86 1.647165e-85 1.0000000
[68,] 5.255999e-87 1.051200e-86 1.0000000
[69,] 9.283622e-88 1.856724e-87 1.0000000
[70,] 3.648310e-89 7.296621e-89 1.0000000
[71,] 2.687256e-90 5.374511e-90 1.0000000
[72,] 2.402323e-91 4.804646e-91 1.0000000
[73,] 9.639550e-93 1.927910e-92 1.0000000
[74,] 4.820992e-94 9.641984e-94 1.0000000
[75,] 1.806918e-95 3.613837e-95 1.0000000
[76,] 6.732685e-97 1.346537e-96 1.0000000
[77,] 3.585972e-98 7.171944e-98 1.0000000
[78,] 1.390185e-99 2.780370e-99 1.0000000
[79,] 3.493286e-100 6.986573e-100 1.0000000
[80,] 1.761459e-101 3.522918e-101 1.0000000
[81,] 6.784160e-103 1.356832e-102 1.0000000
[82,] 2.536517e-104 5.073033e-104 1.0000000
[83,] 1.409978e-103 2.819957e-103 1.0000000
[84,] 7.068659e-105 1.413732e-104 1.0000000
[85,] 3.229139e-106 6.458277e-106 1.0000000
[86,] 1.446914e-106 2.893827e-106 1.0000000
[87,] 3.051793e-107 6.103586e-107 1.0000000
[88,] 1.744967e-108 3.489934e-108 1.0000000
[89,] 7.723433e-110 1.544687e-109 1.0000000
[90,] 4.418697e-111 8.837394e-111 1.0000000
[91,] 2.116472e-112 4.232945e-112 1.0000000
[92,] 8.508776e-114 1.701755e-113 1.0000000
[93,] 3.392209e-115 6.784417e-115 1.0000000
[94,] 4.698125e-116 9.396250e-116 1.0000000
[95,] 3.603223e-117 7.206446e-117 1.0000000
[96,] 1.500654e-118 3.001308e-118 1.0000000
[97,] 2.032175e-119 4.064350e-119 1.0000000
[98,] 2.383781e-120 4.767563e-120 1.0000000
[99,] 1.761120e-120 3.522240e-120 1.0000000
[100,] 1.106105e-121 2.212210e-121 1.0000000
[101,] 7.560369e-123 1.512074e-122 1.0000000
[102,] 5.005268e-124 1.001054e-123 1.0000000
[103,] 2.092806e-125 4.185612e-125 1.0000000
[104,] 1.251661e-126 2.503323e-126 1.0000000
[105,] 9.214105e-124 1.842821e-123 1.0000000
[106,] 6.276641e-125 1.255328e-124 1.0000000
[107,] 3.314541e-126 6.629082e-126 1.0000000
[108,] 5.267080e-127 1.053416e-126 1.0000000
[109,] 3.993400e-128 7.986801e-128 1.0000000
[110,] 4.130642e-120 8.261285e-120 1.0000000
[111,] 2.645363e-77 5.290726e-77 1.0000000
[112,] 1.322809e-72 2.645618e-72 1.0000000
[113,] 5.982041e-73 1.196408e-72 1.0000000
[114,] 3.789176e-72 7.578353e-72 1.0000000
[115,] 3.774723e-72 7.549447e-72 1.0000000
[116,] 3.060124e-55 6.120248e-55 1.0000000
[117,] 7.555561e-56 1.511112e-55 1.0000000
[118,] 3.376827e-55 6.753654e-55 1.0000000
[119,] 8.858302e-56 1.771660e-55 1.0000000
[120,] 3.505119e-55 7.010238e-55 1.0000000
[121,] 2.305083e-53 4.610167e-53 1.0000000
[122,] 1.303405e-52 2.606810e-52 1.0000000
[123,] 3.778927e-51 7.557853e-51 1.0000000
[124,] 1.030991e-49 2.061981e-49 1.0000000
[125,] 5.311077e-50 1.062215e-49 1.0000000
[126,] 1.478016e-49 2.956033e-49 1.0000000
[127,] 2.325478e-49 4.650957e-49 1.0000000
[128,] 2.047748e-44 4.095497e-44 1.0000000
[129,] 7.486846e-45 1.497369e-44 1.0000000
[130,] 4.515167e-45 9.030335e-45 1.0000000
[131,] 7.265594e-43 1.453119e-42 1.0000000
[132,] 2.352501e-42 4.705001e-42 1.0000000
[133,] 5.059467e-42 1.011893e-41 1.0000000
[134,] 2.504647e-42 5.009295e-42 1.0000000
[135,] 3.305375e-42 6.610750e-42 1.0000000
[136,] 1.565775e-42 3.131549e-42 1.0000000
[137,] 4.955379e-43 9.910758e-43 1.0000000
[138,] 4.285922e-43 8.571844e-43 1.0000000
[139,] 2.463213e-43 4.926426e-43 1.0000000
[140,] 1.017550e-43 2.035099e-43 1.0000000
[141,] 3.673115e-44 7.346229e-44 1.0000000
[142,] 1.216331e-44 2.432663e-44 1.0000000
[143,] 4.911573e-45 9.823146e-45 1.0000000
[144,] 1.700076e-45 3.400151e-45 1.0000000
[145,] 1.556621e-45 3.113242e-45 1.0000000
[146,] 9.900166e-45 1.980033e-44 1.0000000
[147,] 3.107996e-45 6.215992e-45 1.0000000
[148,] 9.588485e-46 1.917697e-45 1.0000000
[149,] 3.009395e-46 6.018791e-46 1.0000000
[150,] 5.906502e-45 1.181300e-44 1.0000000
[151,] 2.340035e-45 4.680070e-45 1.0000000
[152,] 7.499160e-46 1.499832e-45 1.0000000
[153,] 2.323312e-46 4.646624e-46 1.0000000
[154,] 1.097479e-46 2.194959e-46 1.0000000
[155,] 3.661364e-47 7.322727e-47 1.0000000
[156,] 8.763361e-46 1.752672e-45 1.0000000
[157,] 2.822278e-46 5.644556e-46 1.0000000
[158,] 5.751113e-45 1.150223e-44 1.0000000
[159,] 1.243529e-44 2.487057e-44 1.0000000
[160,] 7.342386e-43 1.468477e-42 1.0000000
[161,] 2.671592e-43 5.343184e-43 1.0000000
[162,] 9.524254e-44 1.904851e-43 1.0000000
[163,] 9.121532e-41 1.824306e-40 1.0000000
[164,] 3.130895e-41 6.261789e-41 1.0000000
[165,] 1.105166e-41 2.210331e-41 1.0000000
[166,] 1.595877e-41 3.191753e-41 1.0000000
[167,] 5.879796e-42 1.175959e-41 1.0000000
[168,] 3.212270e-41 6.424539e-41 1.0000000
[169,] 1.174804e-41 2.349609e-41 1.0000000
[170,] 2.882687e-39 5.765373e-39 1.0000000
[171,] 1.296953e-39 2.593905e-39 1.0000000
[172,] 4.790291e-40 9.580581e-40 1.0000000
[173,] 1.799462e-40 3.598924e-40 1.0000000
[174,] 6.718088e-41 1.343618e-40 1.0000000
[175,] 9.094078e-41 1.818816e-40 1.0000000
[176,] 3.248470e-41 6.496940e-41 1.0000000
[177,] 1.287924e-41 2.575849e-41 1.0000000
[178,] 5.274591e-42 1.054918e-41 1.0000000
[179,] 1.902405e-42 3.804811e-42 1.0000000
[180,] 1.641496e-42 3.282992e-42 1.0000000
[181,] 4.396947e-42 8.793894e-42 1.0000000
[182,] 3.401486e-41 6.802973e-41 1.0000000
[183,] 5.394803e-41 1.078961e-40 1.0000000
[184,] 2.146630e-41 4.293260e-41 1.0000000
[185,] 1.717925e-41 3.435850e-41 1.0000000
[186,] 2.802453e-41 5.604906e-41 1.0000000
[187,] 2.137743e-41 4.275486e-41 1.0000000
[188,] 2.632300e-41 5.264599e-41 1.0000000
[189,] 2.306838e-40 4.613677e-40 1.0000000
[190,] 8.667482e-41 1.733496e-40 1.0000000
[191,] 4.331874e-41 8.663747e-41 1.0000000
[192,] 2.481017e-40 4.962034e-40 1.0000000
[193,] 2.531208e-40 5.062415e-40 1.0000000
[194,] 1.673615e-40 3.347230e-40 1.0000000
[195,] 9.457078e-41 1.891416e-40 1.0000000
[196,] 4.083661e-41 8.167322e-41 1.0000000
[197,] 1.796109e-41 3.592218e-41 1.0000000
[198,] 1.057754e-41 2.115509e-41 1.0000000
[199,] 6.229303e-42 1.245861e-41 1.0000000
[200,] 2.603782e-42 5.207565e-42 1.0000000
[201,] 1.083872e-42 2.167743e-42 1.0000000
[202,] 5.940884e-43 1.188177e-42 1.0000000
[203,] 2.813539e-43 5.627077e-43 1.0000000
[204,] 1.185200e-43 2.370399e-43 1.0000000
[205,] 5.445832e-44 1.089166e-43 1.0000000
[206,] 2.243254e-44 4.486508e-44 1.0000000
[207,] 1.338579e-44 2.677157e-44 1.0000000
[208,] 5.869129e-45 1.173826e-44 1.0000000
[209,] 3.260264e-45 6.520528e-45 1.0000000
[210,] 1.580990e-45 3.161979e-45 1.0000000
[211,] 7.724555e-46 1.544911e-45 1.0000000
[212,] 6.164559e-46 1.232912e-45 1.0000000
[213,] 2.590198e-46 5.180396e-46 1.0000000
[214,] 1.265240e-46 2.530479e-46 1.0000000
[215,] 6.385374e-47 1.277075e-46 1.0000000
[216,] 3.252761e-47 6.505522e-47 1.0000000
[217,] 1.797304e-47 3.594609e-47 1.0000000
[218,] 7.936599e-48 1.587320e-47 1.0000000
[219,] 3.125412e-48 6.250825e-48 1.0000000
[220,] 3.332817e-48 6.665634e-48 1.0000000
[221,] 1.432390e-48 2.864780e-48 1.0000000
[222,] 9.900397e-49 1.980079e-48 1.0000000
[223,] 4.637763e-49 9.275526e-49 1.0000000
[224,] 2.080619e-49 4.161238e-49 1.0000000
[225,] 9.244681e-50 1.848936e-49 1.0000000
[226,] 3.735437e-50 7.470874e-50 1.0000000
[227,] 2.044702e-50 4.089404e-50 1.0000000
[228,] 1.027426e-50 2.054852e-50 1.0000000
[229,] 6.472889e-51 1.294578e-50 1.0000000
[230,] 1.071136e-50 2.142272e-50 1.0000000
[231,] 7.387830e-51 1.477566e-50 1.0000000
[232,] 3.451128e-51 6.902255e-51 1.0000000
[233,] 2.610218e-51 5.220436e-51 1.0000000
[234,] 1.313689e-51 2.627378e-51 1.0000000
[235,] 1.241035e-51 2.482070e-51 1.0000000
[236,] 2.550309e-51 5.100618e-51 1.0000000
[237,] 4.664030e-51 9.328061e-51 1.0000000
[238,] 2.015304e-51 4.030609e-51 1.0000000
[239,] 7.735322e-52 1.547064e-51 1.0000000
[240,] 3.639205e-52 7.278410e-52 1.0000000
[241,] 1.592741e-52 3.185482e-52 1.0000000
[242,] 9.191696e-53 1.838339e-52 1.0000000
[243,] 3.409371e-53 6.818742e-53 1.0000000
[244,] 2.623818e-53 5.247636e-53 1.0000000
[245,] 1.319334e-53 2.638667e-53 1.0000000
[246,] 6.398299e-54 1.279660e-53 1.0000000
[247,] 2.991017e-54 5.982034e-54 1.0000000
[248,] 2.098904e-54 4.197807e-54 1.0000000
[249,] 1.037587e-54 2.075174e-54 1.0000000
[250,] 8.363986e-55 1.672797e-54 1.0000000
[251,] 3.202987e-54 6.405973e-54 1.0000000
[252,] 1.654817e-54 3.309633e-54 1.0000000
[253,] 8.617070e-55 1.723414e-54 1.0000000
[254,] 3.291747e-55 6.583494e-55 1.0000000
[255,] 1.501439e-55 3.002878e-55 1.0000000
[256,] 7.627110e-56 1.525422e-55 1.0000000
[257,] 2.125163e-48 4.250325e-48 1.0000000
[258,] 3.178707e-47 6.357414e-47 1.0000000
[259,] 2.296769e-46 4.593538e-46 1.0000000
[260,] 2.002304e-46 4.004607e-46 1.0000000
[261,] 2.397074e-46 4.794149e-46 1.0000000
[262,] 1.544384e-46 3.088768e-46 1.0000000
[263,] 2.506706e-46 5.013412e-46 1.0000000
[264,] 1.357599e-46 2.715198e-46 1.0000000
[265,] 8.216384e-47 1.643277e-46 1.0000000
[266,] 6.440593e-47 1.288119e-46 1.0000000
[267,] 4.881168e-47 9.762336e-47 1.0000000
[268,] 2.952400e-47 5.904801e-47 1.0000000
[269,] 3.423234e-46 6.846469e-46 1.0000000
[270,] 4.843771e-46 9.687542e-46 1.0000000
[271,] 3.498292e-46 6.996584e-46 1.0000000
[272,] 1.951967e-46 3.903934e-46 1.0000000
[273,] 1.109858e-46 2.219716e-46 1.0000000
[274,] 6.057741e-47 1.211548e-46 1.0000000
[275,] 5.479594e-47 1.095919e-46 1.0000000
[276,] 3.095164e-47 6.190328e-47 1.0000000
[277,] 1.781366e-47 3.562731e-47 1.0000000
[278,] 9.751067e-48 1.950213e-47 1.0000000
[279,] 5.593750e-48 1.118750e-47 1.0000000
[280,] 2.844735e-48 5.689469e-48 1.0000000
[281,] 1.460369e-48 2.920738e-48 1.0000000
[282,] 7.474104e-49 1.494821e-48 1.0000000
[283,] 4.296606e-49 8.593213e-49 1.0000000
[284,] 2.111721e-49 4.223441e-49 1.0000000
[285,] 1.076593e-49 2.153186e-49 1.0000000
[286,] 5.253226e-50 1.050645e-49 1.0000000
[287,] 2.587785e-50 5.175569e-50 1.0000000
[288,] 1.301562e-50 2.603124e-50 1.0000000
[289,] 6.471480e-51 1.294296e-50 1.0000000
[290,] 3.256433e-51 6.512865e-51 1.0000000
[291,] 1.517278e-51 3.034556e-51 1.0000000
[292,] 7.009445e-52 1.401889e-51 1.0000000
[293,] 3.356955e-52 6.713910e-52 1.0000000
[294,] 1.837628e-52 3.675255e-52 1.0000000
[295,] 9.274973e-53 1.854995e-52 1.0000000
[296,] 4.585944e-53 9.171888e-53 1.0000000
[297,] 2.238796e-53 4.477592e-53 1.0000000
[298,] 1.145793e-53 2.291586e-53 1.0000000
[299,] 5.058749e-54 1.011750e-53 1.0000000
[300,] 2.502432e-54 5.004863e-54 1.0000000
[301,] 1.201520e-54 2.403040e-54 1.0000000
[302,] 6.118044e-55 1.223609e-54 1.0000000
[303,] 3.088754e-55 6.177508e-55 1.0000000
[304,] 1.562873e-55 3.125746e-55 1.0000000
[305,] 7.740894e-56 1.548179e-55 1.0000000
[306,] 3.863961e-56 7.727922e-56 1.0000000
[307,] 2.026731e-56 4.053462e-56 1.0000000
[308,] 1.025796e-56 2.051592e-56 1.0000000
[309,] 5.144014e-57 1.028803e-56 1.0000000
[310,] 2.566950e-57 5.133899e-57 1.0000000
[311,] 1.335847e-57 2.671694e-57 1.0000000
[312,] 1.220593e-57 2.441187e-57 1.0000000
[313,] 5.905658e-55 1.181132e-54 1.0000000
[314,] 5.330670e-53 1.066134e-52 1.0000000
[315,] 2.631926e-47 5.263852e-47 1.0000000
[316,] 3.854748e-38 7.709496e-38 1.0000000
[317,] 1.459304e-36 2.918608e-36 1.0000000
[318,] 4.261263e-33 8.522527e-33 1.0000000
[319,] 1.685535e-31 3.371071e-31 1.0000000
[320,] 7.618045e-28 1.523609e-27 1.0000000
[321,] 6.965265e-27 1.393053e-26 1.0000000
[322,] 1.877419e-22 3.754837e-22 1.0000000
[323,] 1.656676e-22 3.313352e-22 1.0000000
[324,] 1.586369e-22 3.172738e-22 1.0000000
[325,] 2.384253e-22 4.768506e-22 1.0000000
[326,] 1.554020e-22 3.108040e-22 1.0000000
[327,] 1.421653e-22 2.843307e-22 1.0000000
[328,] 2.360544e-22 4.721088e-22 1.0000000
[329,] 1.440775e-22 2.881549e-22 1.0000000
[330,] 8.522924e-23 1.704585e-22 1.0000000
[331,] 6.018362e-23 1.203672e-22 1.0000000
[332,] 3.485496e-23 6.970991e-23 1.0000000
[333,] 2.040560e-23 4.081120e-23 1.0000000
[334,] 1.300154e-23 2.600308e-23 1.0000000
[335,] 7.372601e-24 1.474520e-23 1.0000000
[336,] 4.077728e-24 8.155457e-24 1.0000000
[337,] 2.282177e-24 4.564354e-24 1.0000000
[338,] 1.458244e-24 2.916488e-24 1.0000000
[339,] 3.671343e-24 7.342687e-24 1.0000000
[340,] 6.388388e-24 1.277678e-23 1.0000000
[341,] 3.748113e-24 7.496226e-24 1.0000000
[342,] 2.253056e-24 4.506112e-24 1.0000000
[343,] 1.516826e-24 3.033651e-24 1.0000000
[344,] 1.273703e-24 2.547406e-24 1.0000000
[345,] 2.042522e-24 4.085043e-24 1.0000000
[346,] 1.445189e-24 2.890377e-24 1.0000000
[347,] 1.003510e-24 2.007019e-24 1.0000000
[348,] 5.587168e-25 1.117434e-24 1.0000000
[349,] 3.972736e-25 7.945471e-25 1.0000000
[350,] 2.837612e-25 5.675224e-25 1.0000000
[351,] 5.157458e-25 1.031492e-24 1.0000000
[352,] 3.055097e-25 6.110195e-25 1.0000000
[353,] 2.198226e-25 4.396452e-25 1.0000000
[354,] 1.114443e-25 2.228886e-25 1.0000000
[355,] 5.968394e-26 1.193679e-25 1.0000000
[356,] 3.401125e-26 6.802249e-26 1.0000000
[357,] 1.312708e-25 2.625415e-25 1.0000000
[358,] 2.373391e-25 4.746782e-25 1.0000000
[359,] 1.607996e-25 3.215993e-25 1.0000000
[360,] 3.434121e-25 6.868242e-25 1.0000000
[361,] 2.110903e-25 4.221806e-25 1.0000000
[362,] 1.540536e-25 3.081072e-25 1.0000000
[363,] 7.030922e-24 1.406184e-23 1.0000000
[364,] 1.236944e-23 2.473888e-23 1.0000000
[365,] 7.331745e-24 1.466349e-23 1.0000000
[366,] 1.691306e-22 3.382612e-22 1.0000000
[367,] 1.167515e-22 2.335030e-22 1.0000000
[368,] 1.614271e-22 3.228543e-22 1.0000000
[369,] 2.998714e-22 5.997427e-22 1.0000000
[370,] 1.542643e-22 3.085287e-22 1.0000000
[371,] 1.078491e-22 2.156982e-22 1.0000000
[372,] 5.806475e-23 1.161295e-22 1.0000000
[373,] 4.649975e-23 9.299950e-23 1.0000000
[374,] 6.095673e-23 1.219135e-22 1.0000000
[375,] 1.387774e-22 2.775548e-22 1.0000000
[376,] 8.360638e-23 1.672128e-22 1.0000000
[377,] 5.928595e-23 1.185719e-22 1.0000000
[378,] 5.861222e-23 1.172244e-22 1.0000000
[379,] 3.271638e-23 6.543277e-23 1.0000000
[380,] 2.014787e-20 4.029573e-20 1.0000000
[381,] 3.845382e-20 7.690763e-20 1.0000000
[382,] 2.299894e-20 4.599787e-20 1.0000000
[383,] 1.572570e-20 3.145140e-20 1.0000000
[384,] 1.337207e-20 2.674415e-20 1.0000000
[385,] 9.461613e-21 1.892323e-20 1.0000000
[386,] 2.137392e-20 4.274783e-20 1.0000000
[387,] 3.822508e-20 7.645015e-20 1.0000000
[388,] 2.281092e-20 4.562184e-20 1.0000000
[389,] 1.762084e-20 3.524167e-20 1.0000000
[390,] 1.203139e-20 2.406278e-20 1.0000000
[391,] 8.869180e-21 1.773836e-20 1.0000000
[392,] 6.702166e-21 1.340433e-20 1.0000000
[393,] 1.816282e-20 3.632565e-20 1.0000000
[394,] 6.010416e-19 1.202083e-18 1.0000000
[395,] 5.061987e-19 1.012397e-18 1.0000000
[396,] 3.112359e-19 6.224719e-19 1.0000000
[397,] 6.392993e-19 1.278599e-18 1.0000000
[398,] 4.336757e-19 8.673513e-19 1.0000000
[399,] 3.887862e-19 7.775724e-19 1.0000000
[400,] 3.519071e-19 7.038142e-19 1.0000000
[401,] 7.568782e-19 1.513756e-18 1.0000000
[402,] 4.545865e-19 9.091731e-19 1.0000000
[403,] 1.462722e-18 2.925445e-18 1.0000000
[404,] 1.771067e-18 3.542133e-18 1.0000000
[405,] 3.074962e-17 6.149924e-17 1.0000000
[406,] 5.734496e-17 1.146899e-16 1.0000000
[407,] 6.153234e-17 1.230647e-16 1.0000000
[408,] 7.973174e-17 1.594635e-16 1.0000000
[409,] 6.962728e-17 1.392546e-16 1.0000000
[410,] 1.140238e-16 2.280477e-16 1.0000000
[411,] 1.439790e-16 2.879580e-16 1.0000000
[412,] 1.720992e-16 3.441984e-16 1.0000000
[413,] 2.977940e-16 5.955879e-16 1.0000000
[414,] 5.404336e-16 1.080867e-15 1.0000000
[415,] 1.240752e-15 2.481504e-15 1.0000000
[416,] 2.258337e-15 4.516674e-15 1.0000000
[417,] 3.952639e-14 7.905278e-14 1.0000000
[418,] 5.544892e-14 1.108978e-13 1.0000000
[419,] 8.663569e-14 1.732714e-13 1.0000000
[420,] 7.480613e-14 1.496123e-13 1.0000000
[421,] 1.808237e-13 3.616473e-13 1.0000000
[422,] 3.542845e-13 7.085690e-13 1.0000000
[423,] 7.256124e-13 1.451225e-12 1.0000000
[424,] 1.344083e-12 2.688167e-12 1.0000000
[425,] 2.391010e-12 4.782020e-12 1.0000000
[426,] 1.522393e-12 3.044787e-12 1.0000000
[427,] 7.482763e-12 1.496553e-11 1.0000000
[428,] 1.572917e-11 3.145834e-11 1.0000000
[429,] 1.636907e-11 3.273814e-11 1.0000000
[430,] 1.076788e-10 2.153576e-10 1.0000000
[431,] 6.389257e-11 1.277851e-10 1.0000000
[432,] 1.447557e-10 2.895115e-10 1.0000000
[433,] 3.112861e-10 6.225722e-10 1.0000000
[434,] 4.060839e-09 8.121679e-09 1.0000000
[435,] 5.179924e-09 1.035985e-08 1.0000000
[436,] 4.927312e-09 9.854625e-09 1.0000000
[437,] 8.143015e-08 1.628603e-07 0.9999999
[438,] 6.956544e-08 1.391309e-07 0.9999999
[439,] 3.085827e-07 6.171654e-07 0.9999997
[440,] 6.068255e-07 1.213651e-06 0.9999994
[441,] 9.458210e-05 1.891642e-04 0.9999054
[442,] 7.520017e-05 1.504003e-04 0.9999248
[443,] 4.967267e-03 9.934534e-03 0.9950327
[444,] 3.984766e-03 7.969532e-03 0.9960152
[445,] 3.077845e-03 6.155691e-03 0.9969222
[446,] 1.814194e-02 3.628388e-02 0.9818581
[447,] 1.415379e-02 2.830758e-02 0.9858462
[448,] 1.380801e-02 2.761603e-02 0.9861920
[449,] 1.140729e-02 2.281458e-02 0.9885927
[450,] 1.239458e-02 2.478916e-02 0.9876054
[451,] 1.462012e-02 2.924024e-02 0.9853799
[452,] 9.591274e-03 1.918255e-02 0.9904087
[453,] 2.810221e-02 5.620443e-02 0.9718978
[454,] 1.829481e-02 3.658962e-02 0.9817052
[455,] 3.628732e-02 7.257464e-02 0.9637127
[456,] 1.808695e-02 3.617391e-02 0.9819130
> postscript(file="/var/wessaorg/rcomp/tmp/1a7rn1322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/23nv41322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/3rwn51322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/46qtc1322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/5ldtv1322129486.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 = 467
Frequency = 1
1 2 3 4 5 6
-138.987852 -139.074778 -133.958177 -134.958177 -132.460974 -128.044542
7 8 9 10 11 12
-124.298603 -127.531805 -137.348575 -128.217833 -135.584293 -132.135108
13 14 15 16 17 18
-135.321978 -125.321978 -129.925281 -136.778443 -123.895606 -129.695718
19 20 21 22 23 24
-128.679061 -124.345915 -132.963628 -140.914779 -133.348993 -129.748768
25 26 27 28 29 30
-130.249049 -131.915903 -131.518082 -134.385385 -136.368728 -121.685779
31 32 33 34 35 36
-140.085262 -130.238548 -132.605570 -127.405682 -119.839334 -127.539503
37 38 39 40 41 42
-137.459294 -120.176120 -132.144760 -120.028159 -126.461249 -112.646164
43 44 45 46 47 48
-137.863271 -140.163372 -134.382164 -125.649804 -131.498765 -129.365506
49 50 51 52 53 54
-124.899102 -125.532641 -118.933809 -128.235888 -126.752826 -137.636225
55 56 57 58 59 60
-136.303079 -83.922477 -134.555454 -124.338909 -125.506606 -136.773685
61 62 63 64 65 66
-128.775639 -133.172606 -123.109909 -123.076594 -121.546077 -127.378942
67 68 69 70 71 72
-123.112987 -132.896442 -108.114110 -126.566935 -126.533621 -122.317076
73 74 75 76 77 78
-133.269058 -137.603328 -126.486727 -133.153581 -133.490928 -121.806855
79 80 81 82 83 84
-117.080652 -122.047337 -121.797478 -116.949910 -125.883280 -105.954111
85 86 87 88 89 90
-127.520459 -118.687594 -119.707330 -93.707891 -130.808105 -113.877542
91 92 93 94 95 96
-136.527739 -134.227907 -112.194593 -113.998345 -125.847867 -122.848429
97 98 99 100 101 102
-118.665199 -118.867603 -131.371523 -128.388181 -120.738546 -130.639164
103 104 105 106 107 108
-105.505905 -140.091427 -127.924854 -128.225247 -127.442354 -121.294393
109 110 111 112 113 114
-115.510376 -81.477623 -127.394898 -124.095067 -133.861865 -129.931010
115 116 117 118 119 120
-49.376287 184.567724 75.832561 11.327111 5.846558 -24.885945
121 122 123 124 125 126
282.336900 -43.916600 34.111395 32.461764 19.122732 81.362234
127 128 129 130 131 132
-124.568891 131.754595 130.242276 68.186287 70.523770 53.030768
133 134 135 136 137 138
268.270270 -38.323372 12.535107 207.493115 117.983456 152.927467
139 140 141 142 143 144
19.757952 92.236956 2.587325 -52.021864 27.259630 1.496472
145 146 147 148 149 150
205.456029 1.930824 57.353840 31.663328 -25.648820 -56.451310
151 152 153 154 155 156
-25.858779 -25.127825 32.576684 -24.651481 195.616016 70.050368
157 158 159 160 161 162
1.345859 -3.921637 65.331862 -48.640143 291.487381 24.808206
163 164 165 166 167 168
1.232770 190.752217 295.496058 1.989059 -2.929617 325.802886
169 170 171 172 173 174
35.462743 -13.227768 138.532729 27.772232 -15.028709 -21.004924
175 176 177 178 179 180
323.009074 80.170609 -7.533900 143.270139 -3.536560 182.886455
181 182 183 184 185 186
27.632956 -29.578552 15.830466 17.658290 235.727166 237.061988
187 188 189 190 191 192
279.822486 139.427466 6.680965 -102.628523 -113.684512 51.483455
193 194 195 196 197 198
172.594321 300.265516 131.139540 27.671874 -72.172608 -147.160160
199 200 201 202 203 204
-95.697386 -111.063285 -117.114650 -119.218795 -90.520874 -86.271014
205 206 207 208 209 210
-88.356255 -79.110035 -77.131456 -83.445136 -101.815798 -108.166725
211 212 213 214 215 216
-85.483214 -97.304747 -125.485730 39.062838 -103.837219 -63.838342
217 218 219 220 221 222
-93.662391 -48.131874 -81.482508 -104.884531 -112.487833 -13.303479
223 224 225 226 227 228
-96.291024 -1.459844 -112.549848 -100.697517 -118.198921 -113.283039
229 230 231 232 233 234
-102.484611 -84.838616 -51.055723 -65.439964 37.057800 -75.626835
235 236 237 238 239 240
-117.151052 -103.403428 -85.967530 -119.836518 -94.956197 -89.957882
241 242 243 244 245 246
-113.908472 -95.204641 -51.025343 2.839714 -37.164488 -77.215852
247 248 249 250 251 252
-108.453256 19.378486 -118.110732 -51.761491 -47.064962 -87.786820
253 254 255 256 257 258
-81.505331 -101.588887 -88.249172 -43.462616 -80.989080 64.758251
259 260 261 262 263 264
17.256847 -56.911681 -31.204959 -160.474487 -129.068149 -155.490865
265 266 267 268 269 270
-95.200165 -118.702671 -93.011753 -149.001647 -64.636767 -103.017205
271 272 273 274 275 276
-166.663726 -149.826950 -57.679618 267.307372 -166.565104 -83.608219
277 278 279 280 281 282
-65.749026 -25.907919 -73.666569 -193.889212 -149.915688 -152.195781
283 284 285 286 287 288
-134.822889 -180.276252 -122.237927 -119.182918 -82.622746 -210.655854
289 290 291 292 293 294
-107.664303 -81.716669 -110.950835 -99.778527 -80.442311 -181.080836
295 296 297 298 299 300
-78.832007 -136.623123 -133.163023 -115.416953 -32.108344 -134.042504
301 302 303 304 305 306
-119.261996 -121.129155 -84.686146 -214.081091 -97.725330 -175.863447
307 308 309 310 311 312
-203.787285 -111.025154 -184.984381 -114.462998 -96.924652 -63.727161
313 314 315 316 317 318
-98.314349 -98.686826 -208.965744 -84.375244 -84.062736 363.938685
319 320 321 322 323 324
312.081756 544.191097 808.724801 384.198485 567.791016 431.124289
325 326 327 328 329 330
641.602182 375.290034 779.102312 165.778412 184.313405 -207.373312
331 332 333 334 335 336
8.053187 99.276473 207.061047 -11.244110 -84.986535 81.999176
337 338 339 340 341 342
-72.777687 -21.210926 78.319135 -32.845837 -91.129119 -25.596117
343 344 345 346 347 348
12.339286 298.497221 -234.333263 21.009537 37.787959 123.669066
349 350 351 352 353 354
-74.843944 -184.610444 142.216551 90.040061 -13.761031 146.795041
355 356 357 358 359 360
13.270972 -214.413479 54.140701 151.106078 -48.216460 12.279592
361 362 363 364 365 366
54.622384 383.405916 -205.432800 106.963361 264.520057 107.430318
367 368 369 370 371 372
167.039190 598.536365 -211.763004 72.110741 476.945310 156.237138
373 374 375 376 377 378
-105.632724 -221.779095 9.803792 164.930310 47.693954 204.100517
379 380 381 382 383 384
-63.330376 -286.451095 95.858225 -88.291960 -139.850472 80.460599
385 386 387 388 389 390
777.121618 -223.286512 122.516833 -31.326538 253.234040 161.487039
391 392 393 394 395 396
401.890663 -202.130866 45.207197 230.760794 167.182731 211.320867
397 398 399 400 401 402
45.654090 -335.308291 605.024393 279.082731 106.602164 -245.852709
403 404 405 406 407 408
206.154015 297.056285 -29.280674 -254.341505 174.031959 411.964498
409 410 411 412 413 414
390.275594 593.557809 -269.533265 314.232555 371.619711 305.392629
415 416 417 418 419 420
-221.815483 259.826929 319.131794 421.805421 -308.067259 412.552413
421 422 423 424 425 426
-322.086847 708.066321 409.741870 394.017829 -86.810899 -366.510117
427 428 429 430 431 432
371.442416 371.405048 475.888146 -308.471792 196.662447 548.492186
433 434 435 436 437 438
-301.472512 272.365154 583.428217 42.105815 -365.604880 376.666234
439 440 441 442 443 444
638.828120 -152.682414 198.175195 704.581311 43.368130 -347.258724
445 446 447 448 449 450
293.858436 626.908511 -202.903272 585.697340 -394.730286 502.648226
451 452 453 454 455 456
652.443561 -210.466260 368.169343 -177.070980 -101.347033 -162.841704
457 458 459 460 461 462
130.103276 -150.960551 -386.233241 74.591393 541.546328 -199.565789
463 464 465 466 467
445.526008 660.829231 -344.093711 547.104961 -277.423392
> postscript(file="/var/wessaorg/rcomp/tmp/6h6731322129486.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 = 467
Frequency = 1
lag(myerror, k = 1) myerror
0 -138.987852 NA
1 -139.074778 -138.987852
2 -133.958177 -139.074778
3 -134.958177 -133.958177
4 -132.460974 -134.958177
5 -128.044542 -132.460974
6 -124.298603 -128.044542
7 -127.531805 -124.298603
8 -137.348575 -127.531805
9 -128.217833 -137.348575
10 -135.584293 -128.217833
11 -132.135108 -135.584293
12 -135.321978 -132.135108
13 -125.321978 -135.321978
14 -129.925281 -125.321978
15 -136.778443 -129.925281
16 -123.895606 -136.778443
17 -129.695718 -123.895606
18 -128.679061 -129.695718
19 -124.345915 -128.679061
20 -132.963628 -124.345915
21 -140.914779 -132.963628
22 -133.348993 -140.914779
23 -129.748768 -133.348993
24 -130.249049 -129.748768
25 -131.915903 -130.249049
26 -131.518082 -131.915903
27 -134.385385 -131.518082
28 -136.368728 -134.385385
29 -121.685779 -136.368728
30 -140.085262 -121.685779
31 -130.238548 -140.085262
32 -132.605570 -130.238548
33 -127.405682 -132.605570
34 -119.839334 -127.405682
35 -127.539503 -119.839334
36 -137.459294 -127.539503
37 -120.176120 -137.459294
38 -132.144760 -120.176120
39 -120.028159 -132.144760
40 -126.461249 -120.028159
41 -112.646164 -126.461249
42 -137.863271 -112.646164
43 -140.163372 -137.863271
44 -134.382164 -140.163372
45 -125.649804 -134.382164
46 -131.498765 -125.649804
47 -129.365506 -131.498765
48 -124.899102 -129.365506
49 -125.532641 -124.899102
50 -118.933809 -125.532641
51 -128.235888 -118.933809
52 -126.752826 -128.235888
53 -137.636225 -126.752826
54 -136.303079 -137.636225
55 -83.922477 -136.303079
56 -134.555454 -83.922477
57 -124.338909 -134.555454
58 -125.506606 -124.338909
59 -136.773685 -125.506606
60 -128.775639 -136.773685
61 -133.172606 -128.775639
62 -123.109909 -133.172606
63 -123.076594 -123.109909
64 -121.546077 -123.076594
65 -127.378942 -121.546077
66 -123.112987 -127.378942
67 -132.896442 -123.112987
68 -108.114110 -132.896442
69 -126.566935 -108.114110
70 -126.533621 -126.566935
71 -122.317076 -126.533621
72 -133.269058 -122.317076
73 -137.603328 -133.269058
74 -126.486727 -137.603328
75 -133.153581 -126.486727
76 -133.490928 -133.153581
77 -121.806855 -133.490928
78 -117.080652 -121.806855
79 -122.047337 -117.080652
80 -121.797478 -122.047337
81 -116.949910 -121.797478
82 -125.883280 -116.949910
83 -105.954111 -125.883280
84 -127.520459 -105.954111
85 -118.687594 -127.520459
86 -119.707330 -118.687594
87 -93.707891 -119.707330
88 -130.808105 -93.707891
89 -113.877542 -130.808105
90 -136.527739 -113.877542
91 -134.227907 -136.527739
92 -112.194593 -134.227907
93 -113.998345 -112.194593
94 -125.847867 -113.998345
95 -122.848429 -125.847867
96 -118.665199 -122.848429
97 -118.867603 -118.665199
98 -131.371523 -118.867603
99 -128.388181 -131.371523
100 -120.738546 -128.388181
101 -130.639164 -120.738546
102 -105.505905 -130.639164
103 -140.091427 -105.505905
104 -127.924854 -140.091427
105 -128.225247 -127.924854
106 -127.442354 -128.225247
107 -121.294393 -127.442354
108 -115.510376 -121.294393
109 -81.477623 -115.510376
110 -127.394898 -81.477623
111 -124.095067 -127.394898
112 -133.861865 -124.095067
113 -129.931010 -133.861865
114 -49.376287 -129.931010
115 184.567724 -49.376287
116 75.832561 184.567724
117 11.327111 75.832561
118 5.846558 11.327111
119 -24.885945 5.846558
120 282.336900 -24.885945
121 -43.916600 282.336900
122 34.111395 -43.916600
123 32.461764 34.111395
124 19.122732 32.461764
125 81.362234 19.122732
126 -124.568891 81.362234
127 131.754595 -124.568891
128 130.242276 131.754595
129 68.186287 130.242276
130 70.523770 68.186287
131 53.030768 70.523770
132 268.270270 53.030768
133 -38.323372 268.270270
134 12.535107 -38.323372
135 207.493115 12.535107
136 117.983456 207.493115
137 152.927467 117.983456
138 19.757952 152.927467
139 92.236956 19.757952
140 2.587325 92.236956
141 -52.021864 2.587325
142 27.259630 -52.021864
143 1.496472 27.259630
144 205.456029 1.496472
145 1.930824 205.456029
146 57.353840 1.930824
147 31.663328 57.353840
148 -25.648820 31.663328
149 -56.451310 -25.648820
150 -25.858779 -56.451310
151 -25.127825 -25.858779
152 32.576684 -25.127825
153 -24.651481 32.576684
154 195.616016 -24.651481
155 70.050368 195.616016
156 1.345859 70.050368
157 -3.921637 1.345859
158 65.331862 -3.921637
159 -48.640143 65.331862
160 291.487381 -48.640143
161 24.808206 291.487381
162 1.232770 24.808206
163 190.752217 1.232770
164 295.496058 190.752217
165 1.989059 295.496058
166 -2.929617 1.989059
167 325.802886 -2.929617
168 35.462743 325.802886
169 -13.227768 35.462743
170 138.532729 -13.227768
171 27.772232 138.532729
172 -15.028709 27.772232
173 -21.004924 -15.028709
174 323.009074 -21.004924
175 80.170609 323.009074
176 -7.533900 80.170609
177 143.270139 -7.533900
178 -3.536560 143.270139
179 182.886455 -3.536560
180 27.632956 182.886455
181 -29.578552 27.632956
182 15.830466 -29.578552
183 17.658290 15.830466
184 235.727166 17.658290
185 237.061988 235.727166
186 279.822486 237.061988
187 139.427466 279.822486
188 6.680965 139.427466
189 -102.628523 6.680965
190 -113.684512 -102.628523
191 51.483455 -113.684512
192 172.594321 51.483455
193 300.265516 172.594321
194 131.139540 300.265516
195 27.671874 131.139540
196 -72.172608 27.671874
197 -147.160160 -72.172608
198 -95.697386 -147.160160
199 -111.063285 -95.697386
200 -117.114650 -111.063285
201 -119.218795 -117.114650
202 -90.520874 -119.218795
203 -86.271014 -90.520874
204 -88.356255 -86.271014
205 -79.110035 -88.356255
206 -77.131456 -79.110035
207 -83.445136 -77.131456
208 -101.815798 -83.445136
209 -108.166725 -101.815798
210 -85.483214 -108.166725
211 -97.304747 -85.483214
212 -125.485730 -97.304747
213 39.062838 -125.485730
214 -103.837219 39.062838
215 -63.838342 -103.837219
216 -93.662391 -63.838342
217 -48.131874 -93.662391
218 -81.482508 -48.131874
219 -104.884531 -81.482508
220 -112.487833 -104.884531
221 -13.303479 -112.487833
222 -96.291024 -13.303479
223 -1.459844 -96.291024
224 -112.549848 -1.459844
225 -100.697517 -112.549848
226 -118.198921 -100.697517
227 -113.283039 -118.198921
228 -102.484611 -113.283039
229 -84.838616 -102.484611
230 -51.055723 -84.838616
231 -65.439964 -51.055723
232 37.057800 -65.439964
233 -75.626835 37.057800
234 -117.151052 -75.626835
235 -103.403428 -117.151052
236 -85.967530 -103.403428
237 -119.836518 -85.967530
238 -94.956197 -119.836518
239 -89.957882 -94.956197
240 -113.908472 -89.957882
241 -95.204641 -113.908472
242 -51.025343 -95.204641
243 2.839714 -51.025343
244 -37.164488 2.839714
245 -77.215852 -37.164488
246 -108.453256 -77.215852
247 19.378486 -108.453256
248 -118.110732 19.378486
249 -51.761491 -118.110732
250 -47.064962 -51.761491
251 -87.786820 -47.064962
252 -81.505331 -87.786820
253 -101.588887 -81.505331
254 -88.249172 -101.588887
255 -43.462616 -88.249172
256 -80.989080 -43.462616
257 64.758251 -80.989080
258 17.256847 64.758251
259 -56.911681 17.256847
260 -31.204959 -56.911681
261 -160.474487 -31.204959
262 -129.068149 -160.474487
263 -155.490865 -129.068149
264 -95.200165 -155.490865
265 -118.702671 -95.200165
266 -93.011753 -118.702671
267 -149.001647 -93.011753
268 -64.636767 -149.001647
269 -103.017205 -64.636767
270 -166.663726 -103.017205
271 -149.826950 -166.663726
272 -57.679618 -149.826950
273 267.307372 -57.679618
274 -166.565104 267.307372
275 -83.608219 -166.565104
276 -65.749026 -83.608219
277 -25.907919 -65.749026
278 -73.666569 -25.907919
279 -193.889212 -73.666569
280 -149.915688 -193.889212
281 -152.195781 -149.915688
282 -134.822889 -152.195781
283 -180.276252 -134.822889
284 -122.237927 -180.276252
285 -119.182918 -122.237927
286 -82.622746 -119.182918
287 -210.655854 -82.622746
288 -107.664303 -210.655854
289 -81.716669 -107.664303
290 -110.950835 -81.716669
291 -99.778527 -110.950835
292 -80.442311 -99.778527
293 -181.080836 -80.442311
294 -78.832007 -181.080836
295 -136.623123 -78.832007
296 -133.163023 -136.623123
297 -115.416953 -133.163023
298 -32.108344 -115.416953
299 -134.042504 -32.108344
300 -119.261996 -134.042504
301 -121.129155 -119.261996
302 -84.686146 -121.129155
303 -214.081091 -84.686146
304 -97.725330 -214.081091
305 -175.863447 -97.725330
306 -203.787285 -175.863447
307 -111.025154 -203.787285
308 -184.984381 -111.025154
309 -114.462998 -184.984381
310 -96.924652 -114.462998
311 -63.727161 -96.924652
312 -98.314349 -63.727161
313 -98.686826 -98.314349
314 -208.965744 -98.686826
315 -84.375244 -208.965744
316 -84.062736 -84.375244
317 363.938685 -84.062736
318 312.081756 363.938685
319 544.191097 312.081756
320 808.724801 544.191097
321 384.198485 808.724801
322 567.791016 384.198485
323 431.124289 567.791016
324 641.602182 431.124289
325 375.290034 641.602182
326 779.102312 375.290034
327 165.778412 779.102312
328 184.313405 165.778412
329 -207.373312 184.313405
330 8.053187 -207.373312
331 99.276473 8.053187
332 207.061047 99.276473
333 -11.244110 207.061047
334 -84.986535 -11.244110
335 81.999176 -84.986535
336 -72.777687 81.999176
337 -21.210926 -72.777687
338 78.319135 -21.210926
339 -32.845837 78.319135
340 -91.129119 -32.845837
341 -25.596117 -91.129119
342 12.339286 -25.596117
343 298.497221 12.339286
344 -234.333263 298.497221
345 21.009537 -234.333263
346 37.787959 21.009537
347 123.669066 37.787959
348 -74.843944 123.669066
349 -184.610444 -74.843944
350 142.216551 -184.610444
351 90.040061 142.216551
352 -13.761031 90.040061
353 146.795041 -13.761031
354 13.270972 146.795041
355 -214.413479 13.270972
356 54.140701 -214.413479
357 151.106078 54.140701
358 -48.216460 151.106078
359 12.279592 -48.216460
360 54.622384 12.279592
361 383.405916 54.622384
362 -205.432800 383.405916
363 106.963361 -205.432800
364 264.520057 106.963361
365 107.430318 264.520057
366 167.039190 107.430318
367 598.536365 167.039190
368 -211.763004 598.536365
369 72.110741 -211.763004
370 476.945310 72.110741
371 156.237138 476.945310
372 -105.632724 156.237138
373 -221.779095 -105.632724
374 9.803792 -221.779095
375 164.930310 9.803792
376 47.693954 164.930310
377 204.100517 47.693954
378 -63.330376 204.100517
379 -286.451095 -63.330376
380 95.858225 -286.451095
381 -88.291960 95.858225
382 -139.850472 -88.291960
383 80.460599 -139.850472
384 777.121618 80.460599
385 -223.286512 777.121618
386 122.516833 -223.286512
387 -31.326538 122.516833
388 253.234040 -31.326538
389 161.487039 253.234040
390 401.890663 161.487039
391 -202.130866 401.890663
392 45.207197 -202.130866
393 230.760794 45.207197
394 167.182731 230.760794
395 211.320867 167.182731
396 45.654090 211.320867
397 -335.308291 45.654090
398 605.024393 -335.308291
399 279.082731 605.024393
400 106.602164 279.082731
401 -245.852709 106.602164
402 206.154015 -245.852709
403 297.056285 206.154015
404 -29.280674 297.056285
405 -254.341505 -29.280674
406 174.031959 -254.341505
407 411.964498 174.031959
408 390.275594 411.964498
409 593.557809 390.275594
410 -269.533265 593.557809
411 314.232555 -269.533265
412 371.619711 314.232555
413 305.392629 371.619711
414 -221.815483 305.392629
415 259.826929 -221.815483
416 319.131794 259.826929
417 421.805421 319.131794
418 -308.067259 421.805421
419 412.552413 -308.067259
420 -322.086847 412.552413
421 708.066321 -322.086847
422 409.741870 708.066321
423 394.017829 409.741870
424 -86.810899 394.017829
425 -366.510117 -86.810899
426 371.442416 -366.510117
427 371.405048 371.442416
428 475.888146 371.405048
429 -308.471792 475.888146
430 196.662447 -308.471792
431 548.492186 196.662447
432 -301.472512 548.492186
433 272.365154 -301.472512
434 583.428217 272.365154
435 42.105815 583.428217
436 -365.604880 42.105815
437 376.666234 -365.604880
438 638.828120 376.666234
439 -152.682414 638.828120
440 198.175195 -152.682414
441 704.581311 198.175195
442 43.368130 704.581311
443 -347.258724 43.368130
444 293.858436 -347.258724
445 626.908511 293.858436
446 -202.903272 626.908511
447 585.697340 -202.903272
448 -394.730286 585.697340
449 502.648226 -394.730286
450 652.443561 502.648226
451 -210.466260 652.443561
452 368.169343 -210.466260
453 -177.070980 368.169343
454 -101.347033 -177.070980
455 -162.841704 -101.347033
456 130.103276 -162.841704
457 -150.960551 130.103276
458 -386.233241 -150.960551
459 74.591393 -386.233241
460 541.546328 74.591393
461 -199.565789 541.546328
462 445.526008 -199.565789
463 660.829231 445.526008
464 -344.093711 660.829231
465 547.104961 -344.093711
466 -277.423392 547.104961
467 NA -277.423392
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -139.074778 -138.987852
[2,] -133.958177 -139.074778
[3,] -134.958177 -133.958177
[4,] -132.460974 -134.958177
[5,] -128.044542 -132.460974
[6,] -124.298603 -128.044542
[7,] -127.531805 -124.298603
[8,] -137.348575 -127.531805
[9,] -128.217833 -137.348575
[10,] -135.584293 -128.217833
[11,] -132.135108 -135.584293
[12,] -135.321978 -132.135108
[13,] -125.321978 -135.321978
[14,] -129.925281 -125.321978
[15,] -136.778443 -129.925281
[16,] -123.895606 -136.778443
[17,] -129.695718 -123.895606
[18,] -128.679061 -129.695718
[19,] -124.345915 -128.679061
[20,] -132.963628 -124.345915
[21,] -140.914779 -132.963628
[22,] -133.348993 -140.914779
[23,] -129.748768 -133.348993
[24,] -130.249049 -129.748768
[25,] -131.915903 -130.249049
[26,] -131.518082 -131.915903
[27,] -134.385385 -131.518082
[28,] -136.368728 -134.385385
[29,] -121.685779 -136.368728
[30,] -140.085262 -121.685779
[31,] -130.238548 -140.085262
[32,] -132.605570 -130.238548
[33,] -127.405682 -132.605570
[34,] -119.839334 -127.405682
[35,] -127.539503 -119.839334
[36,] -137.459294 -127.539503
[37,] -120.176120 -137.459294
[38,] -132.144760 -120.176120
[39,] -120.028159 -132.144760
[40,] -126.461249 -120.028159
[41,] -112.646164 -126.461249
[42,] -137.863271 -112.646164
[43,] -140.163372 -137.863271
[44,] -134.382164 -140.163372
[45,] -125.649804 -134.382164
[46,] -131.498765 -125.649804
[47,] -129.365506 -131.498765
[48,] -124.899102 -129.365506
[49,] -125.532641 -124.899102
[50,] -118.933809 -125.532641
[51,] -128.235888 -118.933809
[52,] -126.752826 -128.235888
[53,] -137.636225 -126.752826
[54,] -136.303079 -137.636225
[55,] -83.922477 -136.303079
[56,] -134.555454 -83.922477
[57,] -124.338909 -134.555454
[58,] -125.506606 -124.338909
[59,] -136.773685 -125.506606
[60,] -128.775639 -136.773685
[61,] -133.172606 -128.775639
[62,] -123.109909 -133.172606
[63,] -123.076594 -123.109909
[64,] -121.546077 -123.076594
[65,] -127.378942 -121.546077
[66,] -123.112987 -127.378942
[67,] -132.896442 -123.112987
[68,] -108.114110 -132.896442
[69,] -126.566935 -108.114110
[70,] -126.533621 -126.566935
[71,] -122.317076 -126.533621
[72,] -133.269058 -122.317076
[73,] -137.603328 -133.269058
[74,] -126.486727 -137.603328
[75,] -133.153581 -126.486727
[76,] -133.490928 -133.153581
[77,] -121.806855 -133.490928
[78,] -117.080652 -121.806855
[79,] -122.047337 -117.080652
[80,] -121.797478 -122.047337
[81,] -116.949910 -121.797478
[82,] -125.883280 -116.949910
[83,] -105.954111 -125.883280
[84,] -127.520459 -105.954111
[85,] -118.687594 -127.520459
[86,] -119.707330 -118.687594
[87,] -93.707891 -119.707330
[88,] -130.808105 -93.707891
[89,] -113.877542 -130.808105
[90,] -136.527739 -113.877542
[91,] -134.227907 -136.527739
[92,] -112.194593 -134.227907
[93,] -113.998345 -112.194593
[94,] -125.847867 -113.998345
[95,] -122.848429 -125.847867
[96,] -118.665199 -122.848429
[97,] -118.867603 -118.665199
[98,] -131.371523 -118.867603
[99,] -128.388181 -131.371523
[100,] -120.738546 -128.388181
[101,] -130.639164 -120.738546
[102,] -105.505905 -130.639164
[103,] -140.091427 -105.505905
[104,] -127.924854 -140.091427
[105,] -128.225247 -127.924854
[106,] -127.442354 -128.225247
[107,] -121.294393 -127.442354
[108,] -115.510376 -121.294393
[109,] -81.477623 -115.510376
[110,] -127.394898 -81.477623
[111,] -124.095067 -127.394898
[112,] -133.861865 -124.095067
[113,] -129.931010 -133.861865
[114,] -49.376287 -129.931010
[115,] 184.567724 -49.376287
[116,] 75.832561 184.567724
[117,] 11.327111 75.832561
[118,] 5.846558 11.327111
[119,] -24.885945 5.846558
[120,] 282.336900 -24.885945
[121,] -43.916600 282.336900
[122,] 34.111395 -43.916600
[123,] 32.461764 34.111395
[124,] 19.122732 32.461764
[125,] 81.362234 19.122732
[126,] -124.568891 81.362234
[127,] 131.754595 -124.568891
[128,] 130.242276 131.754595
[129,] 68.186287 130.242276
[130,] 70.523770 68.186287
[131,] 53.030768 70.523770
[132,] 268.270270 53.030768
[133,] -38.323372 268.270270
[134,] 12.535107 -38.323372
[135,] 207.493115 12.535107
[136,] 117.983456 207.493115
[137,] 152.927467 117.983456
[138,] 19.757952 152.927467
[139,] 92.236956 19.757952
[140,] 2.587325 92.236956
[141,] -52.021864 2.587325
[142,] 27.259630 -52.021864
[143,] 1.496472 27.259630
[144,] 205.456029 1.496472
[145,] 1.930824 205.456029
[146,] 57.353840 1.930824
[147,] 31.663328 57.353840
[148,] -25.648820 31.663328
[149,] -56.451310 -25.648820
[150,] -25.858779 -56.451310
[151,] -25.127825 -25.858779
[152,] 32.576684 -25.127825
[153,] -24.651481 32.576684
[154,] 195.616016 -24.651481
[155,] 70.050368 195.616016
[156,] 1.345859 70.050368
[157,] -3.921637 1.345859
[158,] 65.331862 -3.921637
[159,] -48.640143 65.331862
[160,] 291.487381 -48.640143
[161,] 24.808206 291.487381
[162,] 1.232770 24.808206
[163,] 190.752217 1.232770
[164,] 295.496058 190.752217
[165,] 1.989059 295.496058
[166,] -2.929617 1.989059
[167,] 325.802886 -2.929617
[168,] 35.462743 325.802886
[169,] -13.227768 35.462743
[170,] 138.532729 -13.227768
[171,] 27.772232 138.532729
[172,] -15.028709 27.772232
[173,] -21.004924 -15.028709
[174,] 323.009074 -21.004924
[175,] 80.170609 323.009074
[176,] -7.533900 80.170609
[177,] 143.270139 -7.533900
[178,] -3.536560 143.270139
[179,] 182.886455 -3.536560
[180,] 27.632956 182.886455
[181,] -29.578552 27.632956
[182,] 15.830466 -29.578552
[183,] 17.658290 15.830466
[184,] 235.727166 17.658290
[185,] 237.061988 235.727166
[186,] 279.822486 237.061988
[187,] 139.427466 279.822486
[188,] 6.680965 139.427466
[189,] -102.628523 6.680965
[190,] -113.684512 -102.628523
[191,] 51.483455 -113.684512
[192,] 172.594321 51.483455
[193,] 300.265516 172.594321
[194,] 131.139540 300.265516
[195,] 27.671874 131.139540
[196,] -72.172608 27.671874
[197,] -147.160160 -72.172608
[198,] -95.697386 -147.160160
[199,] -111.063285 -95.697386
[200,] -117.114650 -111.063285
[201,] -119.218795 -117.114650
[202,] -90.520874 -119.218795
[203,] -86.271014 -90.520874
[204,] -88.356255 -86.271014
[205,] -79.110035 -88.356255
[206,] -77.131456 -79.110035
[207,] -83.445136 -77.131456
[208,] -101.815798 -83.445136
[209,] -108.166725 -101.815798
[210,] -85.483214 -108.166725
[211,] -97.304747 -85.483214
[212,] -125.485730 -97.304747
[213,] 39.062838 -125.485730
[214,] -103.837219 39.062838
[215,] -63.838342 -103.837219
[216,] -93.662391 -63.838342
[217,] -48.131874 -93.662391
[218,] -81.482508 -48.131874
[219,] -104.884531 -81.482508
[220,] -112.487833 -104.884531
[221,] -13.303479 -112.487833
[222,] -96.291024 -13.303479
[223,] -1.459844 -96.291024
[224,] -112.549848 -1.459844
[225,] -100.697517 -112.549848
[226,] -118.198921 -100.697517
[227,] -113.283039 -118.198921
[228,] -102.484611 -113.283039
[229,] -84.838616 -102.484611
[230,] -51.055723 -84.838616
[231,] -65.439964 -51.055723
[232,] 37.057800 -65.439964
[233,] -75.626835 37.057800
[234,] -117.151052 -75.626835
[235,] -103.403428 -117.151052
[236,] -85.967530 -103.403428
[237,] -119.836518 -85.967530
[238,] -94.956197 -119.836518
[239,] -89.957882 -94.956197
[240,] -113.908472 -89.957882
[241,] -95.204641 -113.908472
[242,] -51.025343 -95.204641
[243,] 2.839714 -51.025343
[244,] -37.164488 2.839714
[245,] -77.215852 -37.164488
[246,] -108.453256 -77.215852
[247,] 19.378486 -108.453256
[248,] -118.110732 19.378486
[249,] -51.761491 -118.110732
[250,] -47.064962 -51.761491
[251,] -87.786820 -47.064962
[252,] -81.505331 -87.786820
[253,] -101.588887 -81.505331
[254,] -88.249172 -101.588887
[255,] -43.462616 -88.249172
[256,] -80.989080 -43.462616
[257,] 64.758251 -80.989080
[258,] 17.256847 64.758251
[259,] -56.911681 17.256847
[260,] -31.204959 -56.911681
[261,] -160.474487 -31.204959
[262,] -129.068149 -160.474487
[263,] -155.490865 -129.068149
[264,] -95.200165 -155.490865
[265,] -118.702671 -95.200165
[266,] -93.011753 -118.702671
[267,] -149.001647 -93.011753
[268,] -64.636767 -149.001647
[269,] -103.017205 -64.636767
[270,] -166.663726 -103.017205
[271,] -149.826950 -166.663726
[272,] -57.679618 -149.826950
[273,] 267.307372 -57.679618
[274,] -166.565104 267.307372
[275,] -83.608219 -166.565104
[276,] -65.749026 -83.608219
[277,] -25.907919 -65.749026
[278,] -73.666569 -25.907919
[279,] -193.889212 -73.666569
[280,] -149.915688 -193.889212
[281,] -152.195781 -149.915688
[282,] -134.822889 -152.195781
[283,] -180.276252 -134.822889
[284,] -122.237927 -180.276252
[285,] -119.182918 -122.237927
[286,] -82.622746 -119.182918
[287,] -210.655854 -82.622746
[288,] -107.664303 -210.655854
[289,] -81.716669 -107.664303
[290,] -110.950835 -81.716669
[291,] -99.778527 -110.950835
[292,] -80.442311 -99.778527
[293,] -181.080836 -80.442311
[294,] -78.832007 -181.080836
[295,] -136.623123 -78.832007
[296,] -133.163023 -136.623123
[297,] -115.416953 -133.163023
[298,] -32.108344 -115.416953
[299,] -134.042504 -32.108344
[300,] -119.261996 -134.042504
[301,] -121.129155 -119.261996
[302,] -84.686146 -121.129155
[303,] -214.081091 -84.686146
[304,] -97.725330 -214.081091
[305,] -175.863447 -97.725330
[306,] -203.787285 -175.863447
[307,] -111.025154 -203.787285
[308,] -184.984381 -111.025154
[309,] -114.462998 -184.984381
[310,] -96.924652 -114.462998
[311,] -63.727161 -96.924652
[312,] -98.314349 -63.727161
[313,] -98.686826 -98.314349
[314,] -208.965744 -98.686826
[315,] -84.375244 -208.965744
[316,] -84.062736 -84.375244
[317,] 363.938685 -84.062736
[318,] 312.081756 363.938685
[319,] 544.191097 312.081756
[320,] 808.724801 544.191097
[321,] 384.198485 808.724801
[322,] 567.791016 384.198485
[323,] 431.124289 567.791016
[324,] 641.602182 431.124289
[325,] 375.290034 641.602182
[326,] 779.102312 375.290034
[327,] 165.778412 779.102312
[328,] 184.313405 165.778412
[329,] -207.373312 184.313405
[330,] 8.053187 -207.373312
[331,] 99.276473 8.053187
[332,] 207.061047 99.276473
[333,] -11.244110 207.061047
[334,] -84.986535 -11.244110
[335,] 81.999176 -84.986535
[336,] -72.777687 81.999176
[337,] -21.210926 -72.777687
[338,] 78.319135 -21.210926
[339,] -32.845837 78.319135
[340,] -91.129119 -32.845837
[341,] -25.596117 -91.129119
[342,] 12.339286 -25.596117
[343,] 298.497221 12.339286
[344,] -234.333263 298.497221
[345,] 21.009537 -234.333263
[346,] 37.787959 21.009537
[347,] 123.669066 37.787959
[348,] -74.843944 123.669066
[349,] -184.610444 -74.843944
[350,] 142.216551 -184.610444
[351,] 90.040061 142.216551
[352,] -13.761031 90.040061
[353,] 146.795041 -13.761031
[354,] 13.270972 146.795041
[355,] -214.413479 13.270972
[356,] 54.140701 -214.413479
[357,] 151.106078 54.140701
[358,] -48.216460 151.106078
[359,] 12.279592 -48.216460
[360,] 54.622384 12.279592
[361,] 383.405916 54.622384
[362,] -205.432800 383.405916
[363,] 106.963361 -205.432800
[364,] 264.520057 106.963361
[365,] 107.430318 264.520057
[366,] 167.039190 107.430318
[367,] 598.536365 167.039190
[368,] -211.763004 598.536365
[369,] 72.110741 -211.763004
[370,] 476.945310 72.110741
[371,] 156.237138 476.945310
[372,] -105.632724 156.237138
[373,] -221.779095 -105.632724
[374,] 9.803792 -221.779095
[375,] 164.930310 9.803792
[376,] 47.693954 164.930310
[377,] 204.100517 47.693954
[378,] -63.330376 204.100517
[379,] -286.451095 -63.330376
[380,] 95.858225 -286.451095
[381,] -88.291960 95.858225
[382,] -139.850472 -88.291960
[383,] 80.460599 -139.850472
[384,] 777.121618 80.460599
[385,] -223.286512 777.121618
[386,] 122.516833 -223.286512
[387,] -31.326538 122.516833
[388,] 253.234040 -31.326538
[389,] 161.487039 253.234040
[390,] 401.890663 161.487039
[391,] -202.130866 401.890663
[392,] 45.207197 -202.130866
[393,] 230.760794 45.207197
[394,] 167.182731 230.760794
[395,] 211.320867 167.182731
[396,] 45.654090 211.320867
[397,] -335.308291 45.654090
[398,] 605.024393 -335.308291
[399,] 279.082731 605.024393
[400,] 106.602164 279.082731
[401,] -245.852709 106.602164
[402,] 206.154015 -245.852709
[403,] 297.056285 206.154015
[404,] -29.280674 297.056285
[405,] -254.341505 -29.280674
[406,] 174.031959 -254.341505
[407,] 411.964498 174.031959
[408,] 390.275594 411.964498
[409,] 593.557809 390.275594
[410,] -269.533265 593.557809
[411,] 314.232555 -269.533265
[412,] 371.619711 314.232555
[413,] 305.392629 371.619711
[414,] -221.815483 305.392629
[415,] 259.826929 -221.815483
[416,] 319.131794 259.826929
[417,] 421.805421 319.131794
[418,] -308.067259 421.805421
[419,] 412.552413 -308.067259
[420,] -322.086847 412.552413
[421,] 708.066321 -322.086847
[422,] 409.741870 708.066321
[423,] 394.017829 409.741870
[424,] -86.810899 394.017829
[425,] -366.510117 -86.810899
[426,] 371.442416 -366.510117
[427,] 371.405048 371.442416
[428,] 475.888146 371.405048
[429,] -308.471792 475.888146
[430,] 196.662447 -308.471792
[431,] 548.492186 196.662447
[432,] -301.472512 548.492186
[433,] 272.365154 -301.472512
[434,] 583.428217 272.365154
[435,] 42.105815 583.428217
[436,] -365.604880 42.105815
[437,] 376.666234 -365.604880
[438,] 638.828120 376.666234
[439,] -152.682414 638.828120
[440,] 198.175195 -152.682414
[441,] 704.581311 198.175195
[442,] 43.368130 704.581311
[443,] -347.258724 43.368130
[444,] 293.858436 -347.258724
[445,] 626.908511 293.858436
[446,] -202.903272 626.908511
[447,] 585.697340 -202.903272
[448,] -394.730286 585.697340
[449,] 502.648226 -394.730286
[450,] 652.443561 502.648226
[451,] -210.466260 652.443561
[452,] 368.169343 -210.466260
[453,] -177.070980 368.169343
[454,] -101.347033 -177.070980
[455,] -162.841704 -101.347033
[456,] 130.103276 -162.841704
[457,] -150.960551 130.103276
[458,] -386.233241 -150.960551
[459,] 74.591393 -386.233241
[460,] 541.546328 74.591393
[461,] -199.565789 541.546328
[462,] 445.526008 -199.565789
[463,] 660.829231 445.526008
[464,] -344.093711 660.829231
[465,] 547.104961 -344.093711
[466,] -277.423392 547.104961
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -139.074778 -138.987852
2 -133.958177 -139.074778
3 -134.958177 -133.958177
4 -132.460974 -134.958177
5 -128.044542 -132.460974
6 -124.298603 -128.044542
7 -127.531805 -124.298603
8 -137.348575 -127.531805
9 -128.217833 -137.348575
10 -135.584293 -128.217833
11 -132.135108 -135.584293
12 -135.321978 -132.135108
13 -125.321978 -135.321978
14 -129.925281 -125.321978
15 -136.778443 -129.925281
16 -123.895606 -136.778443
17 -129.695718 -123.895606
18 -128.679061 -129.695718
19 -124.345915 -128.679061
20 -132.963628 -124.345915
21 -140.914779 -132.963628
22 -133.348993 -140.914779
23 -129.748768 -133.348993
24 -130.249049 -129.748768
25 -131.915903 -130.249049
26 -131.518082 -131.915903
27 -134.385385 -131.518082
28 -136.368728 -134.385385
29 -121.685779 -136.368728
30 -140.085262 -121.685779
31 -130.238548 -140.085262
32 -132.605570 -130.238548
33 -127.405682 -132.605570
34 -119.839334 -127.405682
35 -127.539503 -119.839334
36 -137.459294 -127.539503
37 -120.176120 -137.459294
38 -132.144760 -120.176120
39 -120.028159 -132.144760
40 -126.461249 -120.028159
41 -112.646164 -126.461249
42 -137.863271 -112.646164
43 -140.163372 -137.863271
44 -134.382164 -140.163372
45 -125.649804 -134.382164
46 -131.498765 -125.649804
47 -129.365506 -131.498765
48 -124.899102 -129.365506
49 -125.532641 -124.899102
50 -118.933809 -125.532641
51 -128.235888 -118.933809
52 -126.752826 -128.235888
53 -137.636225 -126.752826
54 -136.303079 -137.636225
55 -83.922477 -136.303079
56 -134.555454 -83.922477
57 -124.338909 -134.555454
58 -125.506606 -124.338909
59 -136.773685 -125.506606
60 -128.775639 -136.773685
61 -133.172606 -128.775639
62 -123.109909 -133.172606
63 -123.076594 -123.109909
64 -121.546077 -123.076594
65 -127.378942 -121.546077
66 -123.112987 -127.378942
67 -132.896442 -123.112987
68 -108.114110 -132.896442
69 -126.566935 -108.114110
70 -126.533621 -126.566935
71 -122.317076 -126.533621
72 -133.269058 -122.317076
73 -137.603328 -133.269058
74 -126.486727 -137.603328
75 -133.153581 -126.486727
76 -133.490928 -133.153581
77 -121.806855 -133.490928
78 -117.080652 -121.806855
79 -122.047337 -117.080652
80 -121.797478 -122.047337
81 -116.949910 -121.797478
82 -125.883280 -116.949910
83 -105.954111 -125.883280
84 -127.520459 -105.954111
85 -118.687594 -127.520459
86 -119.707330 -118.687594
87 -93.707891 -119.707330
88 -130.808105 -93.707891
89 -113.877542 -130.808105
90 -136.527739 -113.877542
91 -134.227907 -136.527739
92 -112.194593 -134.227907
93 -113.998345 -112.194593
94 -125.847867 -113.998345
95 -122.848429 -125.847867
96 -118.665199 -122.848429
97 -118.867603 -118.665199
98 -131.371523 -118.867603
99 -128.388181 -131.371523
100 -120.738546 -128.388181
101 -130.639164 -120.738546
102 -105.505905 -130.639164
103 -140.091427 -105.505905
104 -127.924854 -140.091427
105 -128.225247 -127.924854
106 -127.442354 -128.225247
107 -121.294393 -127.442354
108 -115.510376 -121.294393
109 -81.477623 -115.510376
110 -127.394898 -81.477623
111 -124.095067 -127.394898
112 -133.861865 -124.095067
113 -129.931010 -133.861865
114 -49.376287 -129.931010
115 184.567724 -49.376287
116 75.832561 184.567724
117 11.327111 75.832561
118 5.846558 11.327111
119 -24.885945 5.846558
120 282.336900 -24.885945
121 -43.916600 282.336900
122 34.111395 -43.916600
123 32.461764 34.111395
124 19.122732 32.461764
125 81.362234 19.122732
126 -124.568891 81.362234
127 131.754595 -124.568891
128 130.242276 131.754595
129 68.186287 130.242276
130 70.523770 68.186287
131 53.030768 70.523770
132 268.270270 53.030768
133 -38.323372 268.270270
134 12.535107 -38.323372
135 207.493115 12.535107
136 117.983456 207.493115
137 152.927467 117.983456
138 19.757952 152.927467
139 92.236956 19.757952
140 2.587325 92.236956
141 -52.021864 2.587325
142 27.259630 -52.021864
143 1.496472 27.259630
144 205.456029 1.496472
145 1.930824 205.456029
146 57.353840 1.930824
147 31.663328 57.353840
148 -25.648820 31.663328
149 -56.451310 -25.648820
150 -25.858779 -56.451310
151 -25.127825 -25.858779
152 32.576684 -25.127825
153 -24.651481 32.576684
154 195.616016 -24.651481
155 70.050368 195.616016
156 1.345859 70.050368
157 -3.921637 1.345859
158 65.331862 -3.921637
159 -48.640143 65.331862
160 291.487381 -48.640143
161 24.808206 291.487381
162 1.232770 24.808206
163 190.752217 1.232770
164 295.496058 190.752217
165 1.989059 295.496058
166 -2.929617 1.989059
167 325.802886 -2.929617
168 35.462743 325.802886
169 -13.227768 35.462743
170 138.532729 -13.227768
171 27.772232 138.532729
172 -15.028709 27.772232
173 -21.004924 -15.028709
174 323.009074 -21.004924
175 80.170609 323.009074
176 -7.533900 80.170609
177 143.270139 -7.533900
178 -3.536560 143.270139
179 182.886455 -3.536560
180 27.632956 182.886455
181 -29.578552 27.632956
182 15.830466 -29.578552
183 17.658290 15.830466
184 235.727166 17.658290
185 237.061988 235.727166
186 279.822486 237.061988
187 139.427466 279.822486
188 6.680965 139.427466
189 -102.628523 6.680965
190 -113.684512 -102.628523
191 51.483455 -113.684512
192 172.594321 51.483455
193 300.265516 172.594321
194 131.139540 300.265516
195 27.671874 131.139540
196 -72.172608 27.671874
197 -147.160160 -72.172608
198 -95.697386 -147.160160
199 -111.063285 -95.697386
200 -117.114650 -111.063285
201 -119.218795 -117.114650
202 -90.520874 -119.218795
203 -86.271014 -90.520874
204 -88.356255 -86.271014
205 -79.110035 -88.356255
206 -77.131456 -79.110035
207 -83.445136 -77.131456
208 -101.815798 -83.445136
209 -108.166725 -101.815798
210 -85.483214 -108.166725
211 -97.304747 -85.483214
212 -125.485730 -97.304747
213 39.062838 -125.485730
214 -103.837219 39.062838
215 -63.838342 -103.837219
216 -93.662391 -63.838342
217 -48.131874 -93.662391
218 -81.482508 -48.131874
219 -104.884531 -81.482508
220 -112.487833 -104.884531
221 -13.303479 -112.487833
222 -96.291024 -13.303479
223 -1.459844 -96.291024
224 -112.549848 -1.459844
225 -100.697517 -112.549848
226 -118.198921 -100.697517
227 -113.283039 -118.198921
228 -102.484611 -113.283039
229 -84.838616 -102.484611
230 -51.055723 -84.838616
231 -65.439964 -51.055723
232 37.057800 -65.439964
233 -75.626835 37.057800
234 -117.151052 -75.626835
235 -103.403428 -117.151052
236 -85.967530 -103.403428
237 -119.836518 -85.967530
238 -94.956197 -119.836518
239 -89.957882 -94.956197
240 -113.908472 -89.957882
241 -95.204641 -113.908472
242 -51.025343 -95.204641
243 2.839714 -51.025343
244 -37.164488 2.839714
245 -77.215852 -37.164488
246 -108.453256 -77.215852
247 19.378486 -108.453256
248 -118.110732 19.378486
249 -51.761491 -118.110732
250 -47.064962 -51.761491
251 -87.786820 -47.064962
252 -81.505331 -87.786820
253 -101.588887 -81.505331
254 -88.249172 -101.588887
255 -43.462616 -88.249172
256 -80.989080 -43.462616
257 64.758251 -80.989080
258 17.256847 64.758251
259 -56.911681 17.256847
260 -31.204959 -56.911681
261 -160.474487 -31.204959
262 -129.068149 -160.474487
263 -155.490865 -129.068149
264 -95.200165 -155.490865
265 -118.702671 -95.200165
266 -93.011753 -118.702671
267 -149.001647 -93.011753
268 -64.636767 -149.001647
269 -103.017205 -64.636767
270 -166.663726 -103.017205
271 -149.826950 -166.663726
272 -57.679618 -149.826950
273 267.307372 -57.679618
274 -166.565104 267.307372
275 -83.608219 -166.565104
276 -65.749026 -83.608219
277 -25.907919 -65.749026
278 -73.666569 -25.907919
279 -193.889212 -73.666569
280 -149.915688 -193.889212
281 -152.195781 -149.915688
282 -134.822889 -152.195781
283 -180.276252 -134.822889
284 -122.237927 -180.276252
285 -119.182918 -122.237927
286 -82.622746 -119.182918
287 -210.655854 -82.622746
288 -107.664303 -210.655854
289 -81.716669 -107.664303
290 -110.950835 -81.716669
291 -99.778527 -110.950835
292 -80.442311 -99.778527
293 -181.080836 -80.442311
294 -78.832007 -181.080836
295 -136.623123 -78.832007
296 -133.163023 -136.623123
297 -115.416953 -133.163023
298 -32.108344 -115.416953
299 -134.042504 -32.108344
300 -119.261996 -134.042504
301 -121.129155 -119.261996
302 -84.686146 -121.129155
303 -214.081091 -84.686146
304 -97.725330 -214.081091
305 -175.863447 -97.725330
306 -203.787285 -175.863447
307 -111.025154 -203.787285
308 -184.984381 -111.025154
309 -114.462998 -184.984381
310 -96.924652 -114.462998
311 -63.727161 -96.924652
312 -98.314349 -63.727161
313 -98.686826 -98.314349
314 -208.965744 -98.686826
315 -84.375244 -208.965744
316 -84.062736 -84.375244
317 363.938685 -84.062736
318 312.081756 363.938685
319 544.191097 312.081756
320 808.724801 544.191097
321 384.198485 808.724801
322 567.791016 384.198485
323 431.124289 567.791016
324 641.602182 431.124289
325 375.290034 641.602182
326 779.102312 375.290034
327 165.778412 779.102312
328 184.313405 165.778412
329 -207.373312 184.313405
330 8.053187 -207.373312
331 99.276473 8.053187
332 207.061047 99.276473
333 -11.244110 207.061047
334 -84.986535 -11.244110
335 81.999176 -84.986535
336 -72.777687 81.999176
337 -21.210926 -72.777687
338 78.319135 -21.210926
339 -32.845837 78.319135
340 -91.129119 -32.845837
341 -25.596117 -91.129119
342 12.339286 -25.596117
343 298.497221 12.339286
344 -234.333263 298.497221
345 21.009537 -234.333263
346 37.787959 21.009537
347 123.669066 37.787959
348 -74.843944 123.669066
349 -184.610444 -74.843944
350 142.216551 -184.610444
351 90.040061 142.216551
352 -13.761031 90.040061
353 146.795041 -13.761031
354 13.270972 146.795041
355 -214.413479 13.270972
356 54.140701 -214.413479
357 151.106078 54.140701
358 -48.216460 151.106078
359 12.279592 -48.216460
360 54.622384 12.279592
361 383.405916 54.622384
362 -205.432800 383.405916
363 106.963361 -205.432800
364 264.520057 106.963361
365 107.430318 264.520057
366 167.039190 107.430318
367 598.536365 167.039190
368 -211.763004 598.536365
369 72.110741 -211.763004
370 476.945310 72.110741
371 156.237138 476.945310
372 -105.632724 156.237138
373 -221.779095 -105.632724
374 9.803792 -221.779095
375 164.930310 9.803792
376 47.693954 164.930310
377 204.100517 47.693954
378 -63.330376 204.100517
379 -286.451095 -63.330376
380 95.858225 -286.451095
381 -88.291960 95.858225
382 -139.850472 -88.291960
383 80.460599 -139.850472
384 777.121618 80.460599
385 -223.286512 777.121618
386 122.516833 -223.286512
387 -31.326538 122.516833
388 253.234040 -31.326538
389 161.487039 253.234040
390 401.890663 161.487039
391 -202.130866 401.890663
392 45.207197 -202.130866
393 230.760794 45.207197
394 167.182731 230.760794
395 211.320867 167.182731
396 45.654090 211.320867
397 -335.308291 45.654090
398 605.024393 -335.308291
399 279.082731 605.024393
400 106.602164 279.082731
401 -245.852709 106.602164
402 206.154015 -245.852709
403 297.056285 206.154015
404 -29.280674 297.056285
405 -254.341505 -29.280674
406 174.031959 -254.341505
407 411.964498 174.031959
408 390.275594 411.964498
409 593.557809 390.275594
410 -269.533265 593.557809
411 314.232555 -269.533265
412 371.619711 314.232555
413 305.392629 371.619711
414 -221.815483 305.392629
415 259.826929 -221.815483
416 319.131794 259.826929
417 421.805421 319.131794
418 -308.067259 421.805421
419 412.552413 -308.067259
420 -322.086847 412.552413
421 708.066321 -322.086847
422 409.741870 708.066321
423 394.017829 409.741870
424 -86.810899 394.017829
425 -366.510117 -86.810899
426 371.442416 -366.510117
427 371.405048 371.442416
428 475.888146 371.405048
429 -308.471792 475.888146
430 196.662447 -308.471792
431 548.492186 196.662447
432 -301.472512 548.492186
433 272.365154 -301.472512
434 583.428217 272.365154
435 42.105815 583.428217
436 -365.604880 42.105815
437 376.666234 -365.604880
438 638.828120 376.666234
439 -152.682414 638.828120
440 198.175195 -152.682414
441 704.581311 198.175195
442 43.368130 704.581311
443 -347.258724 43.368130
444 293.858436 -347.258724
445 626.908511 293.858436
446 -202.903272 626.908511
447 585.697340 -202.903272
448 -394.730286 585.697340
449 502.648226 -394.730286
450 652.443561 502.648226
451 -210.466260 652.443561
452 368.169343 -210.466260
453 -177.070980 368.169343
454 -101.347033 -177.070980
455 -162.841704 -101.347033
456 130.103276 -162.841704
457 -150.960551 130.103276
458 -386.233241 -150.960551
459 74.591393 -386.233241
460 541.546328 74.591393
461 -199.565789 541.546328
462 445.526008 -199.565789
463 660.829231 445.526008
464 -344.093711 660.829231
465 547.104961 -344.093711
466 -277.423392 547.104961
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/72djq1322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/8ogez1322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/9o6mr1322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/wessaorg/rcomp/tmp/10thwl1322129486.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/11djrk1322129486.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/129kxv1322129487.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/13yqqg1322129487.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/14wyjz1322129487.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/15b2p01322129487.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/16j4071322129487.tab")
+ }
>
> try(system("convert tmp/1a7rn1322129486.ps tmp/1a7rn1322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/23nv41322129486.ps tmp/23nv41322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/3rwn51322129486.ps tmp/3rwn51322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/46qtc1322129486.ps tmp/46qtc1322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ldtv1322129486.ps tmp/5ldtv1322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/6h6731322129486.ps tmp/6h6731322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/72djq1322129486.ps tmp/72djq1322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/8ogez1322129486.ps tmp/8ogez1322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/9o6mr1322129486.ps tmp/9o6mr1322129486.png",intern=TRUE))
character(0)
> try(system("convert tmp/10thwl1322129486.ps tmp/10thwl1322129486.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.821 0.572 13.426