R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(162556
+ ,807
+ ,6282154
+ ,29790
+ ,444
+ ,4321023
+ ,87550
+ ,412
+ ,4111912
+ ,84738
+ ,428
+ ,223193
+ ,54660
+ ,315
+ ,1491348
+ ,42634
+ ,168
+ ,1629616
+ ,40949
+ ,263
+ ,1398893
+ ,45187
+ ,267
+ ,1926517
+ ,37704
+ ,228
+ ,983660
+ ,16275
+ ,129
+ ,1443586
+ ,25830
+ ,104
+ ,1073089
+ ,12679
+ ,122
+ ,984885
+ ,18014
+ ,393
+ ,1405225
+ ,43556
+ ,190
+ ,227132
+ ,24811
+ ,280
+ ,929118
+ ,6575
+ ,63
+ ,1071292
+ ,7123
+ ,102
+ ,638830
+ ,21950
+ ,265
+ ,856956
+ ,37597
+ ,234
+ ,992426
+ ,17821
+ ,277
+ ,444477
+ ,12988
+ ,73
+ ,857217
+ ,22330
+ ,67
+ ,711969
+ ,13326
+ ,103
+ ,702380
+ ,16189
+ ,290
+ ,358589
+ ,7146
+ ,83
+ ,297978
+ ,15824
+ ,56
+ ,585715
+ ,27664
+ ,236
+ ,657954
+ ,11920
+ ,73
+ ,209458
+ ,8568
+ ,34
+ ,786690
+ ,14416
+ ,139
+ ,439798
+ ,3369
+ ,26
+ ,688779
+ ,11819
+ ,70
+ ,574339
+ ,6984
+ ,40
+ ,741409
+ ,4519
+ ,42
+ ,597793
+ ,2220
+ ,12
+ ,644190
+ ,18562
+ ,211
+ ,377934
+ ,10327
+ ,74
+ ,640273
+ ,5336
+ ,80
+ ,697458
+ ,2365
+ ,83
+ ,550608
+ ,4069
+ ,131
+ ,207393
+ ,8636
+ ,203
+ ,301607
+ ,13718
+ ,56
+ ,345783
+ ,4525
+ ,89
+ ,501749
+ ,6869
+ ,88
+ ,379983
+ ,4628
+ ,39
+ ,387475
+ ,3689
+ ,25
+ ,377305
+ ,4891
+ ,49
+ ,370837
+ ,7489
+ ,149
+ ,430866
+ ,4901
+ ,58
+ ,469107
+ ,2284
+ ,41
+ ,194493
+ ,3160
+ ,90
+ ,530670
+ ,4150
+ ,136
+ ,518365
+ ,7285
+ ,97
+ ,491303
+ ,1134
+ ,63
+ ,527021
+ ,4658
+ ,114
+ ,233773
+ ,2384
+ ,77
+ ,405972
+ ,3748
+ ,6
+ ,652925
+ ,5371
+ ,47
+ ,446211
+ ,1285
+ ,51
+ ,341340
+ ,9327
+ ,85
+ ,387699
+ ,5565
+ ,43
+ ,493408
+ ,1528
+ ,32
+ ,146494
+ ,3122
+ ,25
+ ,414462
+ ,7561
+ ,77
+ ,364304
+ ,2675
+ ,54
+ ,355178
+ ,13253
+ ,251
+ ,357760
+ ,880
+ ,15
+ ,261216
+ ,2053
+ ,44
+ ,397144
+ ,1424
+ ,73
+ ,374943
+ ,4036
+ ,85
+ ,424898
+ ,3045
+ ,49
+ ,202055
+ ,5119
+ ,38
+ ,378525
+ ,1431
+ ,35
+ ,310768
+ ,554
+ ,9
+ ,325738
+ ,1975
+ ,34
+ ,394510
+ ,1765
+ ,20
+ ,247060
+ ,1012
+ ,29
+ ,368078
+ ,810
+ ,11
+ ,236761
+ ,1280
+ ,52
+ ,312378
+ ,666
+ ,13
+ ,339836
+ ,1380
+ ,29
+ ,347385
+ ,4677
+ ,66
+ ,426280
+ ,876
+ ,33
+ ,352850
+ ,814
+ ,15
+ ,301881
+ ,514
+ ,15
+ ,377516
+ ,5692
+ ,68
+ ,357312
+ ,3642
+ ,100
+ ,458343
+ ,540
+ ,13
+ ,354228
+ ,2099
+ ,45
+ ,308636
+ ,567
+ ,14
+ ,386212
+ ,2001
+ ,36
+ ,393343
+ ,2949
+ ,40
+ ,378509
+ ,2253
+ ,68
+ ,452469
+ ,6533
+ ,29
+ ,364839
+ ,1889
+ ,43
+ ,358649
+ ,3055
+ ,30
+ ,376641
+ ,272
+ ,9
+ ,429112
+ ,1414
+ ,22
+ ,330546
+ ,2564
+ ,19
+ ,403560
+ ,1383
+ ,9
+ ,317892
+ ,1261
+ ,31
+ ,307528
+ ,975
+ ,19
+ ,235133
+ ,3366
+ ,55
+ ,299243
+ ,576
+ ,8
+ ,314073
+ ,1686
+ ,28
+ ,368186
+ ,746
+ ,29
+ ,269661
+ ,3192
+ ,48
+ ,125390
+ ,2045
+ ,16
+ ,510834
+ ,5702
+ ,47
+ ,321896
+ ,1932
+ ,20
+ ,249898
+ ,936
+ ,22
+ ,408881
+ ,3437
+ ,33
+ ,158492
+ ,5131
+ ,44
+ ,292154
+ ,2397
+ ,13
+ ,289513
+ ,1389
+ ,6
+ ,378049
+ ,1503
+ ,35
+ ,343466
+ ,402
+ ,8
+ ,332743
+ ,2239
+ ,17
+ ,442882
+ ,2234
+ ,11
+ ,214215
+ ,837
+ ,21
+ ,315688
+ ,10579
+ ,92
+ ,375195
+ ,875
+ ,12
+ ,334280
+ ,1585
+ ,112
+ ,355864
+ ,1659
+ ,25
+ ,480382
+ ,2647
+ ,17
+ ,353058
+ ,3294
+ ,23
+ ,217193
+ ,0
+ ,0
+ ,315380
+ ,94
+ ,10
+ ,314533
+ ,422
+ ,23
+ ,318056
+ ,0
+ ,0
+ ,315380
+ ,34
+ ,7
+ ,314353
+ ,1558
+ ,25
+ ,369448
+ ,0
+ ,0
+ ,315380
+ ,43
+ ,20
+ ,312846
+ ,645
+ ,4
+ ,312075
+ ,316
+ ,4
+ ,315009
+ ,115
+ ,10
+ ,318903
+ ,5
+ ,1
+ ,314887
+ ,897
+ ,4
+ ,314913
+ ,0
+ ,0
+ ,315380
+ ,389
+ ,8
+ ,325506
+ ,0
+ ,0
+ ,315380
+ ,1002
+ ,11
+ ,298568
+ ,36
+ ,4
+ ,315834
+ ,460
+ ,15
+ ,329784
+ ,309
+ ,9
+ ,312878
+ ,0
+ ,0
+ ,315380
+ ,9
+ ,7
+ ,314987
+ ,271
+ ,2
+ ,325249
+ ,14
+ ,0
+ ,315877
+ ,520
+ ,7
+ ,291650
+ ,1766
+ ,46
+ ,305959
+ ,0
+ ,5
+ ,315380
+ ,458
+ ,7
+ ,297765
+ ,20
+ ,2
+ ,315245
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,98
+ ,2
+ ,315236
+ ,405
+ ,5
+ ,336425
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,483
+ ,7
+ ,306268
+ ,454
+ ,24
+ ,302187
+ ,47
+ ,1
+ ,314882
+ ,0
+ ,0
+ ,315380
+ ,757
+ ,18
+ ,382712
+ ,4655
+ ,55
+ ,341570
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,36
+ ,3
+ ,312412
+ ,0
+ ,0
+ ,315380
+ ,203
+ ,9
+ ,309596
+ ,0
+ ,0
+ ,315380
+ ,126
+ ,8
+ ,315547
+ ,400
+ ,113
+ ,313267
+ ,71
+ ,0
+ ,316176
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,972
+ ,19
+ ,359335
+ ,531
+ ,11
+ ,330068
+ ,2461
+ ,25
+ ,314289
+ ,378
+ ,16
+ ,297413
+ ,23
+ ,5
+ ,314806
+ ,638
+ ,11
+ ,333210
+ ,2300
+ ,23
+ ,352108
+ ,149
+ ,6
+ ,313332
+ ,226
+ ,5
+ ,291787
+ ,0
+ ,0
+ ,315380
+ ,275
+ ,7
+ ,318745
+ ,0
+ ,0
+ ,315380
+ ,141
+ ,7
+ ,315366
+ ,0
+ ,0
+ ,315380
+ ,28
+ ,3
+ ,315688
+ ,0
+ ,0
+ ,315380
+ ,4980
+ ,89
+ ,409642
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,472
+ ,19
+ ,269587
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,203
+ ,12
+ ,300962
+ ,496
+ ,12
+ ,325479
+ ,10
+ ,5
+ ,316155
+ ,63
+ ,2
+ ,318574
+ ,0
+ ,0
+ ,315380
+ ,1136
+ ,26
+ ,343613
+ ,265
+ ,3
+ ,306948
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,267
+ ,11
+ ,330059
+ ,474
+ ,10
+ ,288985
+ ,534
+ ,5
+ ,304485
+ ,0
+ ,2
+ ,315380
+ ,15
+ ,6
+ ,315688
+ ,397
+ ,7
+ ,317736
+ ,0
+ ,2
+ ,315380
+ ,1866
+ ,28
+ ,322331
+ ,288
+ ,3
+ ,296656
+ ,0
+ ,0
+ ,315380
+ ,3
+ ,1
+ ,315354
+ ,468
+ ,20
+ ,312161
+ ,20
+ ,1
+ ,315576
+ ,278
+ ,22
+ ,314922
+ ,61
+ ,9
+ ,314551
+ ,0
+ ,0
+ ,315380
+ ,192
+ ,2
+ ,312339
+ ,0
+ ,0
+ ,315380
+ ,317
+ ,7
+ ,298700
+ ,738
+ ,9
+ ,321376
+ ,0
+ ,0
+ ,315380
+ ,368
+ ,13
+ ,303230
+ ,0
+ ,0
+ ,315380
+ ,2
+ ,0
+ ,315487
+ ,0
+ ,0
+ ,315380
+ ,53
+ ,6
+ ,315793
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,94
+ ,3
+ ,312887
+ ,0
+ ,0
+ ,315380
+ ,24
+ ,7
+ ,315637
+ ,2332
+ ,2
+ ,324385
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,131
+ ,15
+ ,308989
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,206
+ ,9
+ ,296702
+ ,0
+ ,0
+ ,315380
+ ,167
+ ,1
+ ,307322
+ ,622
+ ,38
+ ,304376
+ ,2328
+ ,57
+ ,253588
+ ,0
+ ,0
+ ,315380
+ ,365
+ ,7
+ ,309560
+ ,364
+ ,26
+ ,298466
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,226
+ ,13
+ ,343929
+ ,307
+ ,10
+ ,331955
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,188
+ ,9
+ ,381180
+ ,0
+ ,0
+ ,315380
+ ,138
+ ,26
+ ,331420
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,125
+ ,19
+ ,310201
+ ,0
+ ,0
+ ,315380
+ ,282
+ ,12
+ ,320016
+ ,335
+ ,23
+ ,320398
+ ,0
+ ,0
+ ,315380
+ ,1324
+ ,29
+ ,291841
+ ,176
+ ,8
+ ,310670
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,249
+ ,26
+ ,313491
+ ,0
+ ,0
+ ,315380
+ ,333
+ ,9
+ ,331323
+ ,0
+ ,0
+ ,315380
+ ,601
+ ,5
+ ,319210
+ ,30
+ ,3
+ ,318098
+ ,0
+ ,0
+ ,315380
+ ,249
+ ,13
+ ,292754
+ ,0
+ ,0
+ ,315380
+ ,165
+ ,12
+ ,325176
+ ,453
+ ,19
+ ,365959
+ ,0
+ ,0
+ ,315380
+ ,53
+ ,10
+ ,302409
+ ,382
+ ,9
+ ,340968
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,9
+ ,315380
+ ,30
+ ,4
+ ,313164
+ ,290
+ ,1
+ ,301164
+ ,0
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,366
+ ,14
+ ,344425
+ ,2
+ ,12
+ ,315394
+ ,0
+ ,0
+ ,315380
+ ,209
+ ,19
+ ,316647
+ ,384
+ ,17
+ ,309836
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,365
+ ,32
+ ,346611
+ ,0
+ ,0
+ ,315380
+ ,49
+ ,14
+ ,322031
+ ,3
+ ,8
+ ,315656
+ ,133
+ ,4
+ ,339445
+ ,32
+ ,0
+ ,314964
+ ,368
+ ,20
+ ,297141
+ ,1
+ ,5
+ ,315372
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,22
+ ,1
+ ,312502
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,96
+ ,4
+ ,313729
+ ,1
+ ,1
+ ,315388
+ ,314
+ ,4
+ ,315371
+ ,844
+ ,20
+ ,296139
+ ,0
+ ,0
+ ,315380
+ ,26
+ ,1
+ ,313880
+ ,125
+ ,10
+ ,317698
+ ,304
+ ,12
+ ,295580
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,621
+ ,13
+ ,308256
+ ,0
+ ,0
+ ,315380
+ ,119
+ ,3
+ ,303677
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1595
+ ,10
+ ,319369
+ ,312
+ ,3
+ ,318690
+ ,60
+ ,7
+ ,314049
+ ,587
+ ,10
+ ,325699
+ ,135
+ ,1
+ ,314210
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,514
+ ,15
+ ,322378
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,4
+ ,315398
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1763
+ ,28
+ ,308336
+ ,180
+ ,9
+ ,316386
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,218
+ ,7
+ ,315553
+ ,0
+ ,0
+ ,315380
+ ,448
+ ,7
+ ,323361
+ ,227
+ ,7
+ ,336639
+ ,174
+ ,3
+ ,307424
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,121
+ ,11
+ ,295370
+ ,607
+ ,7
+ ,322340
+ ,2212
+ ,10
+ ,319864
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,530
+ ,18
+ ,317291
+ ,571
+ ,14
+ ,280398
+ ,0
+ ,0
+ ,315380
+ ,78
+ ,12
+ ,317330
+ ,2489
+ ,29
+ ,238125
+ ,131
+ ,3
+ ,327071
+ ,923
+ ,6
+ ,309038
+ ,72
+ ,3
+ ,314210
+ ,572
+ ,8
+ ,307930
+ ,397
+ ,10
+ ,322327
+ ,450
+ ,6
+ ,292136
+ ,622
+ ,8
+ ,263276
+ ,694
+ ,6
+ ,367655
+ ,3425
+ ,9
+ ,283910
+ ,562
+ ,8
+ ,283587
+ ,4917
+ ,26
+ ,243650
+ ,1442
+ ,239
+ ,438493
+ ,529
+ ,7
+ ,296261
+ ,2126
+ ,41
+ ,230621
+ ,1061
+ ,3
+ ,304252
+ ,776
+ ,8
+ ,333505
+ ,611
+ ,6
+ ,296919
+ ,1526
+ ,21
+ ,278990
+ ,592
+ ,7
+ ,276898
+ ,1182
+ ,11
+ ,327007
+ ,621
+ ,11
+ ,317046
+ ,989
+ ,12
+ ,304555
+ ,438
+ ,9
+ ,298096
+ ,726
+ ,3
+ ,231861
+ ,1303
+ ,57
+ ,309422
+ ,7419
+ ,21
+ ,286963
+ ,1164
+ ,15
+ ,269753
+ ,3310
+ ,32
+ ,448243
+ ,1920
+ ,11
+ ,165404
+ ,965
+ ,2
+ ,204325
+ ,3256
+ ,23
+ ,407159
+ ,1135
+ ,20
+ ,290476
+ ,1270
+ ,24
+ ,275311
+ ,661
+ ,1
+ ,246541
+ ,1013
+ ,1
+ ,253468
+ ,2844
+ ,74
+ ,240897
+ ,11528
+ ,68
+ ,-83265
+ ,6526
+ ,20
+ ,-42143
+ ,2264
+ ,20
+ ,272713
+ ,5109
+ ,82
+ ,215362
+ ,3999
+ ,21
+ ,42754
+ ,35624
+ ,244
+ ,306275
+ ,9252
+ ,32
+ ,253537
+ ,15236
+ ,86
+ ,372631
+ ,18073
+ ,69
+ ,-7170)
+ ,dim=c(3
+ ,431)
+ ,dimnames=list(c('kosten'
+ ,'orders'
+ ,'rijkdom')
+ ,1:431))
> y <- array(NA,dim=c(3,431),dimnames=list(c('kosten','orders','rijkdom'),1:431))
> 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
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
rijkdom kosten orders
1 6282154 162556 807
2 4321023 29790 444
3 4111912 87550 412
4 223193 84738 428
5 1491348 54660 315
6 1629616 42634 168
7 1398893 40949 263
8 1926517 45187 267
9 983660 37704 228
10 1443586 16275 129
11 1073089 25830 104
12 984885 12679 122
13 1405225 18014 393
14 227132 43556 190
15 929118 24811 280
16 1071292 6575 63
17 638830 7123 102
18 856956 21950 265
19 992426 37597 234
20 444477 17821 277
21 857217 12988 73
22 711969 22330 67
23 702380 13326 103
24 358589 16189 290
25 297978 7146 83
26 585715 15824 56
27 657954 27664 236
28 209458 11920 73
29 786690 8568 34
30 439798 14416 139
31 688779 3369 26
32 574339 11819 70
33 741409 6984 40
34 597793 4519 42
35 644190 2220 12
36 377934 18562 211
37 640273 10327 74
38 697458 5336 80
39 550608 2365 83
40 207393 4069 131
41 301607 8636 203
42 345783 13718 56
43 501749 4525 89
44 379983 6869 88
45 387475 4628 39
46 377305 3689 25
47 370837 4891 49
48 430866 7489 149
49 469107 4901 58
50 194493 2284 41
51 530670 3160 90
52 518365 4150 136
53 491303 7285 97
54 527021 1134 63
55 233773 4658 114
56 405972 2384 77
57 652925 3748 6
58 446211 5371 47
59 341340 1285 51
60 387699 9327 85
61 493408 5565 43
62 146494 1528 32
63 414462 3122 25
64 364304 7561 77
65 355178 2675 54
66 357760 13253 251
67 261216 880 15
68 397144 2053 44
69 374943 1424 73
70 424898 4036 85
71 202055 3045 49
72 378525 5119 38
73 310768 1431 35
74 325738 554 9
75 394510 1975 34
76 247060 1765 20
77 368078 1012 29
78 236761 810 11
79 312378 1280 52
80 339836 666 13
81 347385 1380 29
82 426280 4677 66
83 352850 876 33
84 301881 814 15
85 377516 514 15
86 357312 5692 68
87 458343 3642 100
88 354228 540 13
89 308636 2099 45
90 386212 567 14
91 393343 2001 36
92 378509 2949 40
93 452469 2253 68
94 364839 6533 29
95 358649 1889 43
96 376641 3055 30
97 429112 272 9
98 330546 1414 22
99 403560 2564 19
100 317892 1383 9
101 307528 1261 31
102 235133 975 19
103 299243 3366 55
104 314073 576 8
105 368186 1686 28
106 269661 746 29
107 125390 3192 48
108 510834 2045 16
109 321896 5702 47
110 249898 1932 20
111 408881 936 22
112 158492 3437 33
113 292154 5131 44
114 289513 2397 13
115 378049 1389 6
116 343466 1503 35
117 332743 402 8
118 442882 2239 17
119 214215 2234 11
120 315688 837 21
121 375195 10579 92
122 334280 875 12
123 355864 1585 112
124 480382 1659 25
125 353058 2647 17
126 217193 3294 23
127 315380 0 0
128 314533 94 10
129 318056 422 23
130 315380 0 0
131 314353 34 7
132 369448 1558 25
133 315380 0 0
134 312846 43 20
135 312075 645 4
136 315009 316 4
137 318903 115 10
138 314887 5 1
139 314913 897 4
140 315380 0 0
141 325506 389 8
142 315380 0 0
143 298568 1002 11
144 315834 36 4
145 329784 460 15
146 312878 309 9
147 315380 0 0
148 314987 9 7
149 325249 271 2
150 315877 14 0
151 291650 520 7
152 305959 1766 46
153 315380 0 5
154 297765 458 7
155 315245 20 2
156 315380 0 0
157 315380 0 0
158 315236 98 2
159 336425 405 5
160 315380 0 0
161 315380 0 0
162 315380 0 0
163 315380 0 0
164 306268 483 7
165 302187 454 24
166 314882 47 1
167 315380 0 0
168 382712 757 18
169 341570 4655 55
170 315380 0 0
171 315380 0 0
172 312412 36 3
173 315380 0 0
174 309596 203 9
175 315380 0 0
176 315547 126 8
177 313267 400 113
178 316176 71 0
179 315380 0 0
180 315380 0 0
181 359335 972 19
182 330068 531 11
183 314289 2461 25
184 297413 378 16
185 314806 23 5
186 333210 638 11
187 352108 2300 23
188 313332 149 6
189 291787 226 5
190 315380 0 0
191 318745 275 7
192 315380 0 0
193 315366 141 7
194 315380 0 0
195 315688 28 3
196 315380 0 0
197 409642 4980 89
198 315380 0 0
199 315380 0 0
200 269587 472 19
201 315380 0 0
202 315380 0 0
203 315380 0 0
204 300962 203 12
205 325479 496 12
206 316155 10 5
207 318574 63 2
208 315380 0 0
209 343613 1136 26
210 306948 265 3
211 315380 0 0
212 315380 0 0
213 330059 267 11
214 288985 474 10
215 304485 534 5
216 315380 0 2
217 315688 15 6
218 317736 397 7
219 315380 0 2
220 322331 1866 28
221 296656 288 3
222 315380 0 0
223 315354 3 1
224 312161 468 20
225 315576 20 1
226 314922 278 22
227 314551 61 9
228 315380 0 0
229 312339 192 2
230 315380 0 0
231 298700 317 7
232 321376 738 9
233 315380 0 0
234 303230 368 13
235 315380 0 0
236 315487 2 0
237 315380 0 0
238 315793 53 6
239 315380 0 0
240 315380 0 0
241 315380 0 0
242 312887 94 3
243 315380 0 0
244 315637 24 7
245 324385 2332 2
246 315380 0 0
247 315380 0 0
248 308989 131 15
249 315380 0 0
250 315380 0 0
251 296702 206 9
252 315380 0 0
253 307322 167 1
254 304376 622 38
255 253588 2328 57
256 315380 0 0
257 309560 365 7
258 298466 364 26
259 315380 0 0
260 315380 0 0
261 315380 0 0
262 315380 0 0
263 343929 226 13
264 331955 307 10
265 315380 0 0
266 315380 0 0
267 315380 0 0
268 381180 188 9
269 315380 0 0
270 331420 138 26
271 315380 0 0
272 315380 0 0
273 315380 0 0
274 310201 125 19
275 315380 0 0
276 320016 282 12
277 320398 335 23
278 315380 0 0
279 291841 1324 29
280 310670 176 8
281 315380 0 0
282 315380 0 0
283 313491 249 26
284 315380 0 0
285 331323 333 9
286 315380 0 0
287 319210 601 5
288 318098 30 3
289 315380 0 0
290 292754 249 13
291 315380 0 0
292 325176 165 12
293 365959 453 19
294 315380 0 0
295 302409 53 10
296 340968 382 9
297 315380 0 0
298 315380 0 0
299 315380 0 0
300 315380 0 9
301 313164 30 4
302 301164 290 1
303 315380 0 1
304 315380 0 0
305 344425 366 14
306 315394 2 12
307 315380 0 0
308 316647 209 19
309 309836 384 17
310 315380 0 0
311 315380 0 0
312 346611 365 32
313 315380 0 0
314 322031 49 14
315 315656 3 8
316 339445 133 4
317 314964 32 0
318 297141 368 20
319 315372 1 5
320 315380 0 0
321 315380 0 0
322 315380 0 0
323 315380 0 0
324 315380 0 0
325 315380 0 0
326 312502 22 1
327 315380 0 0
328 315380 0 0
329 315380 0 0
330 315380 0 0
331 315380 0 0
332 315380 0 0
333 315380 0 0
334 313729 96 4
335 315388 1 1
336 315371 314 4
337 296139 844 20
338 315380 0 0
339 313880 26 1
340 317698 125 10
341 295580 304 12
342 315380 0 0
343 315380 0 0
344 315380 0 0
345 308256 621 13
346 315380 0 0
347 303677 119 3
348 315380 0 0
349 315380 0 0
350 319369 1595 10
351 318690 312 3
352 314049 60 7
353 325699 587 10
354 314210 135 1
355 315380 0 0
356 315380 0 0
357 322378 514 15
358 315380 0 0
359 315380 0 0
360 315380 0 0
361 315398 1 4
362 315380 0 0
363 315380 0 0
364 308336 1763 28
365 316386 180 9
366 315380 0 0
367 315380 0 0
368 315380 0 0
369 315380 0 0
370 315553 218 7
371 315380 0 0
372 323361 448 7
373 336639 227 7
374 307424 174 3
375 315380 0 0
376 315380 0 0
377 295370 121 11
378 322340 607 7
379 319864 2212 10
380 315380 0 0
381 315380 0 0
382 317291 530 18
383 280398 571 14
384 315380 0 0
385 317330 78 12
386 238125 2489 29
387 327071 131 3
388 309038 923 6
389 314210 72 3
390 307930 572 8
391 322327 397 10
392 292136 450 6
393 263276 622 8
394 367655 694 6
395 283910 3425 9
396 283587 562 8
397 243650 4917 26
398 438493 1442 239
399 296261 529 7
400 230621 2126 41
401 304252 1061 3
402 333505 776 8
403 296919 611 6
404 278990 1526 21
405 276898 592 7
406 327007 1182 11
407 317046 621 11
408 304555 989 12
409 298096 438 9
410 231861 726 3
411 309422 1303 57
412 286963 7419 21
413 269753 1164 15
414 448243 3310 32
415 165404 1920 11
416 204325 965 2
417 407159 3256 23
418 290476 1135 20
419 275311 1270 24
420 246541 661 1
421 253468 1013 1
422 240897 2844 74
423 -83265 11528 68
424 -42143 6526 20
425 272713 2264 20
426 215362 5109 82
427 42754 3999 21
428 306275 35624 244
429 253537 9252 32
430 372631 15236 86
431 -7170 18073 69
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) kosten orders
269503.75 21.83 1290.12
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2448703 -8351 34265 45876 2828255
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.695e+05 1.388e+04 19.423 < 2e-16 ***
kosten 2.184e+01 2.188e+00 9.979 < 2e-16 ***
orders 1.290e+03 3.472e+02 3.716 0.000229 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 257500 on 428 degrees of freedom
Multiple R-squared: 0.6439, Adjusted R-squared: 0.6422
F-statistic: 387 on 2 and 428 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 3.031121e-83 1.515560e-83
[2,] 1.0000000 7.131567e-84 3.565783e-84
[3,] 1.0000000 2.374157e-91 1.187079e-91
[4,] 1.0000000 2.005009e-91 1.002504e-91
[5,] 1.0000000 8.553319e-103 4.276660e-103
[6,] 1.0000000 9.376522e-108 4.688261e-108
[7,] 1.0000000 2.438239e-110 1.219119e-110
[8,] 1.0000000 1.251548e-128 6.257742e-129
[9,] 1.0000000 4.637866e-136 2.318933e-136
[10,] 1.0000000 1.224253e-139 6.121264e-140
[11,] 1.0000000 2.282291e-151 1.141145e-151
[12,] 1.0000000 1.302023e-151 6.510116e-152
[13,] 1.0000000 4.153631e-155 2.076816e-155
[14,] 1.0000000 3.358898e-158 1.679449e-158
[15,] 1.0000000 8.154764e-162 4.077382e-162
[16,] 1.0000000 3.041757e-167 1.520879e-167
[17,] 1.0000000 6.439966e-169 3.219983e-169
[18,] 1.0000000 8.009691e-171 4.004845e-171
[19,] 1.0000000 1.523329e-174 7.616643e-175
[20,] 1.0000000 1.307145e-173 6.535727e-174
[21,] 1.0000000 2.520764e-174 1.260382e-174
[22,] 1.0000000 3.349988e-176 1.674994e-176
[23,] 1.0000000 8.574509e-176 4.287255e-176
[24,] 1.0000000 4.261454e-183 2.130727e-183
[25,] 1.0000000 1.961060e-182 9.805302e-183
[26,] 1.0000000 2.904256e-187 1.452128e-187
[27,] 1.0000000 1.191407e-188 5.957033e-189
[28,] 1.0000000 1.056506e-195 5.282530e-196
[29,] 1.0000000 3.178940e-198 1.589470e-198
[30,] 1.0000000 5.465698e-203 2.732849e-203
[31,] 1.0000000 1.492737e-203 7.463685e-204
[32,] 1.0000000 7.124162e-208 3.562081e-208
[33,] 1.0000000 1.312931e-213 6.564656e-214
[34,] 1.0000000 1.606901e-214 8.034503e-215
[35,] 1.0000000 9.667649e-216 4.833824e-216
[36,] 1.0000000 1.631686e-216 8.158431e-217
[37,] 1.0000000 5.859008e-216 2.929504e-216
[38,] 1.0000000 4.555031e-216 2.277516e-216
[39,] 1.0000000 4.292153e-215 2.146076e-215
[40,] 1.0000000 2.224175e-214 1.112087e-214
[41,] 1.0000000 1.162643e-213 5.813215e-214
[42,] 1.0000000 9.458867e-213 4.729434e-213
[43,] 1.0000000 5.942340e-212 2.971170e-212
[44,] 1.0000000 3.438009e-212 1.719005e-212
[45,] 1.0000000 3.270241e-212 1.635121e-212
[46,] 1.0000000 7.167938e-213 3.583969e-213
[47,] 1.0000000 7.790233e-213 3.895117e-213
[48,] 1.0000000 1.981194e-213 9.905972e-214
[49,] 1.0000000 1.464028e-214 7.320138e-215
[50,] 1.0000000 9.801019e-215 4.900510e-215
[51,] 1.0000000 7.250524e-214 3.625262e-214
[52,] 1.0000000 6.733777e-223 3.366889e-223
[53,] 1.0000000 2.130061e-223 1.065031e-223
[54,] 1.0000000 2.616428e-222 1.308214e-222
[55,] 1.0000000 6.340337e-222 3.170169e-222
[56,] 1.0000000 5.267335e-224 2.633668e-224
[57,] 1.0000000 1.847421e-225 9.237107e-226
[58,] 1.0000000 1.739899e-225 8.699497e-226
[59,] 1.0000000 9.231271e-225 4.615635e-225
[60,] 1.0000000 1.001221e-223 5.006107e-224
[61,] 1.0000000 5.188583e-224 2.594291e-224
[62,] 1.0000000 3.527322e-223 1.763661e-223
[63,] 1.0000000 1.567839e-222 7.839196e-223
[64,] 1.0000000 1.781478e-221 8.907390e-222
[65,] 1.0000000 5.434636e-221 2.717318e-221
[66,] 1.0000000 3.807592e-221 1.903796e-221
[67,] 1.0000000 1.051930e-220 5.259649e-221
[68,] 1.0000000 1.173747e-219 5.868737e-220
[69,] 1.0000000 1.157295e-218 5.786474e-219
[70,] 1.0000000 4.012713e-218 2.006356e-218
[71,] 1.0000000 2.139418e-217 1.069709e-217
[72,] 1.0000000 1.529183e-216 7.645913e-217
[73,] 1.0000000 5.646611e-216 2.823305e-216
[74,] 1.0000000 5.536921e-215 2.768461e-215
[75,] 1.0000000 4.917277e-214 2.458639e-214
[76,] 1.0000000 4.578538e-213 2.289269e-213
[77,] 1.0000000 5.020193e-213 2.510097e-213
[78,] 1.0000000 4.679761e-212 2.339881e-212
[79,] 1.0000000 4.780762e-211 2.390381e-211
[80,] 1.0000000 2.180355e-210 1.090177e-210
[81,] 1.0000000 1.426621e-209 7.133106e-210
[82,] 1.0000000 1.430879e-209 7.154395e-210
[83,] 1.0000000 1.007994e-208 5.039971e-209
[84,] 1.0000000 1.053723e-207 5.268613e-208
[85,] 1.0000000 3.484210e-207 1.742105e-207
[86,] 1.0000000 1.095539e-206 5.477695e-207
[87,] 1.0000000 4.543387e-206 2.271694e-206
[88,] 1.0000000 2.669290e-206 1.334645e-206
[89,] 1.0000000 4.322523e-206 2.161261e-206
[90,] 1.0000000 3.352014e-205 1.676007e-205
[91,] 1.0000000 1.054615e-204 5.273073e-205
[92,] 1.0000000 5.049718e-205 2.524859e-205
[93,] 1.0000000 4.789167e-204 2.394584e-204
[94,] 1.0000000 4.358908e-204 2.179454e-204
[95,] 1.0000000 4.201895e-203 2.100947e-203
[96,] 1.0000000 4.297084e-202 2.148542e-202
[97,] 1.0000000 1.243226e-201 6.216129e-202
[98,] 1.0000000 1.201700e-200 6.008500e-201
[99,] 1.0000000 1.167627e-199 5.838133e-200
[100,] 1.0000000 5.869277e-199 2.934638e-199
[101,] 1.0000000 3.405227e-198 1.702614e-198
[102,] 1.0000000 1.425923e-200 7.129617e-201
[103,] 1.0000000 4.431270e-204 2.215635e-204
[104,] 1.0000000 3.331433e-203 1.665716e-203
[105,] 1.0000000 1.738070e-202 8.690348e-203
[106,] 1.0000000 2.015645e-202 1.007823e-202
[107,] 1.0000000 1.503619e-203 7.518093e-204
[108,] 1.0000000 1.518862e-202 7.594310e-203
[109,] 1.0000000 1.591812e-201 7.959060e-202
[110,] 1.0000000 3.404485e-201 1.702242e-201
[111,] 1.0000000 3.055910e-200 1.527955e-200
[112,] 1.0000000 2.646932e-199 1.323466e-199
[113,] 1.0000000 1.391162e-200 6.955808e-201
[114,] 1.0000000 2.738460e-200 1.369230e-200
[115,] 1.0000000 2.827946e-199 1.413973e-199
[116,] 1.0000000 1.056223e-199 5.281113e-200
[117,] 1.0000000 8.913536e-199 4.456768e-199
[118,] 1.0000000 8.572013e-198 4.286006e-198
[119,] 1.0000000 3.105122e-200 1.552561e-200
[120,] 1.0000000 1.233619e-199 6.168094e-200
[121,] 1.0000000 2.829858e-199 1.414929e-199
[122,] 1.0000000 2.751382e-198 1.375691e-198
[123,] 1.0000000 2.766446e-197 1.383223e-197
[124,] 1.0000000 2.838145e-196 1.419072e-196
[125,] 1.0000000 2.729189e-195 1.364594e-195
[126,] 1.0000000 2.689797e-194 1.344898e-194
[127,] 1.0000000 1.021587e-193 5.107936e-194
[128,] 1.0000000 9.745014e-193 4.872507e-193
[129,] 1.0000000 9.557103e-192 4.778552e-192
[130,] 1.0000000 9.268271e-191 4.634135e-191
[131,] 1.0000000 8.875615e-190 4.437807e-190
[132,] 1.0000000 8.540957e-189 4.270479e-189
[133,] 1.0000000 8.042221e-188 4.021110e-188
[134,] 1.0000000 7.477459e-187 3.738729e-187
[135,] 1.0000000 6.944076e-186 3.472038e-186
[136,] 1.0000000 6.189332e-185 3.094666e-185
[137,] 1.0000000 5.704780e-184 2.852390e-184
[138,] 1.0000000 5.458611e-183 2.729306e-183
[139,] 1.0000000 5.055158e-182 2.527579e-182
[140,] 1.0000000 4.453346e-181 2.226673e-181
[141,] 1.0000000 4.176469e-180 2.088234e-180
[142,] 1.0000000 3.771655e-179 1.885828e-179
[143,] 1.0000000 3.471104e-178 1.735552e-178
[144,] 1.0000000 2.895904e-177 1.447952e-177
[145,] 1.0000000 2.578447e-176 1.289223e-176
[146,] 1.0000000 2.242139e-175 1.121070e-175
[147,] 1.0000000 1.967866e-174 9.839331e-175
[148,] 1.0000000 1.762457e-173 8.812283e-174
[149,] 1.0000000 1.563083e-172 7.815413e-173
[150,] 1.0000000 1.376205e-171 6.881023e-172
[151,] 1.0000000 1.197838e-170 5.989192e-171
[152,] 1.0000000 1.038451e-169 5.192256e-170
[153,] 1.0000000 9.025131e-169 4.512565e-169
[154,] 1.0000000 6.395145e-168 3.197573e-168
[155,] 1.0000000 5.479831e-167 2.739915e-167
[156,] 1.0000000 4.676536e-166 2.338268e-166
[157,] 1.0000000 3.974795e-165 1.987397e-165
[158,] 1.0000000 3.364585e-164 1.682292e-164
[159,] 1.0000000 2.928573e-163 1.464287e-163
[160,] 1.0000000 2.385561e-162 1.192781e-162
[161,] 1.0000000 2.003256e-161 1.001628e-161
[162,] 1.0000000 1.667338e-160 8.336690e-161
[163,] 1.0000000 4.228547e-160 2.114274e-160
[164,] 1.0000000 2.193356e-159 1.096678e-159
[165,] 1.0000000 1.811629e-158 9.058143e-159
[166,] 1.0000000 1.490160e-157 7.450801e-158
[167,] 1.0000000 1.241878e-156 6.209390e-157
[168,] 1.0000000 1.012994e-155 5.064972e-156
[169,] 1.0000000 8.456779e-155 4.228390e-155
[170,] 1.0000000 6.839812e-154 3.419906e-154
[171,] 1.0000000 5.620298e-153 2.810149e-153
[172,] 1.0000000 1.473969e-152 7.369844e-153
[173,] 1.0000000 1.174411e-151 5.872056e-152
[174,] 1.0000000 9.369974e-151 4.684987e-151
[175,] 1.0000000 7.444461e-150 3.722230e-150
[176,] 1.0000000 3.510510e-149 1.755255e-149
[177,] 1.0000000 2.580623e-148 1.290312e-148
[178,] 1.0000000 2.002368e-147 1.001184e-147
[179,] 1.0000000 1.513582e-146 7.567909e-147
[180,] 1.0000000 1.198769e-145 5.993847e-146
[181,] 1.0000000 8.291007e-145 4.145504e-145
[182,] 1.0000000 3.549216e-144 1.774608e-144
[183,] 1.0000000 2.794352e-143 1.397176e-143
[184,] 1.0000000 2.081553e-142 1.040776e-142
[185,] 1.0000000 1.590676e-141 7.953378e-142
[186,] 1.0000000 1.214343e-140 6.071716e-141
[187,] 1.0000000 9.200386e-140 4.600193e-140
[188,] 1.0000000 7.070252e-139 3.535126e-139
[189,] 1.0000000 5.310867e-138 2.655434e-138
[190,] 1.0000000 4.005299e-137 2.002650e-137
[191,] 1.0000000 2.982943e-136 1.491471e-136
[192,] 1.0000000 3.311066e-136 1.655533e-136
[193,] 1.0000000 2.465500e-135 1.232750e-135
[194,] 1.0000000 1.828035e-134 9.140175e-135
[195,] 1.0000000 8.884029e-134 4.442015e-134
[196,] 1.0000000 6.545488e-133 3.272744e-133
[197,] 1.0000000 4.801884e-132 2.400942e-132
[198,] 1.0000000 3.507646e-131 1.753823e-131
[199,] 1.0000000 2.547654e-130 1.273827e-130
[200,] 1.0000000 1.794754e-129 8.973768e-130
[201,] 1.0000000 1.309192e-128 6.545959e-129
[202,] 1.0000000 9.310178e-128 4.655089e-128
[203,] 1.0000000 6.656406e-127 3.328203e-127
[204,] 1.0000000 3.972948e-126 1.986474e-126
[205,] 1.0000000 2.885767e-125 1.442883e-125
[206,] 1.0000000 2.037376e-124 1.018688e-124
[207,] 1.0000000 1.432171e-123 7.160856e-124
[208,] 1.0000000 9.495477e-123 4.747739e-123
[209,] 1.0000000 6.313217e-122 3.156609e-122
[210,] 1.0000000 4.503926e-121 2.251963e-121
[211,] 1.0000000 3.131531e-120 1.565766e-120
[212,] 1.0000000 2.186161e-119 1.093081e-119
[213,] 1.0000000 1.497999e-118 7.489996e-119
[214,] 1.0000000 1.027789e-117 5.138946e-118
[215,] 1.0000000 6.805576e-117 3.402788e-117
[216,] 1.0000000 4.651250e-116 2.325625e-116
[217,] 1.0000000 3.131166e-115 1.565583e-115
[218,] 1.0000000 2.104950e-114 1.052475e-114
[219,] 1.0000000 1.443321e-113 7.216607e-114
[220,] 1.0000000 9.604622e-113 4.802311e-113
[221,] 1.0000000 6.517694e-112 3.258847e-112
[222,] 1.0000000 4.380233e-111 2.190117e-111
[223,] 1.0000000 2.871270e-110 1.435635e-110
[224,] 1.0000000 1.895349e-109 9.476744e-110
[225,] 1.0000000 1.231377e-108 6.156884e-109
[226,] 1.0000000 8.052550e-108 4.026275e-108
[227,] 1.0000000 5.036219e-107 2.518110e-107
[228,] 1.0000000 3.228984e-106 1.614492e-106
[229,] 1.0000000 2.092817e-105 1.046409e-105
[230,] 1.0000000 1.330031e-104 6.650154e-105
[231,] 1.0000000 8.410791e-104 4.205396e-104
[232,] 1.0000000 5.297394e-103 2.648697e-103
[233,] 1.0000000 3.363310e-102 1.681655e-102
[234,] 1.0000000 2.099244e-101 1.049622e-101
[235,] 1.0000000 1.304346e-100 6.521729e-101
[236,] 1.0000000 8.067742e-100 4.033871e-100
[237,] 1.0000000 5.033656e-99 2.516828e-99
[238,] 1.0000000 3.085202e-98 1.542601e-98
[239,] 1.0000000 1.910990e-97 9.554949e-98
[240,] 1.0000000 8.739258e-97 4.369629e-97
[241,] 1.0000000 5.287413e-96 2.643706e-96
[242,] 1.0000000 3.184315e-95 1.592158e-95
[243,] 1.0000000 1.949304e-94 9.746520e-95
[244,] 1.0000000 1.163345e-93 5.816726e-94
[245,] 1.0000000 6.910807e-93 3.455403e-93
[246,] 1.0000000 4.082010e-92 2.041005e-92
[247,] 1.0000000 2.403511e-91 1.201756e-91
[248,] 1.0000000 1.437097e-90 7.185487e-91
[249,] 1.0000000 7.976304e-90 3.988152e-90
[250,] 1.0000000 2.028640e-89 1.014320e-89
[251,] 1.0000000 1.179327e-88 5.896637e-89
[252,] 1.0000000 6.974840e-88 3.487420e-88
[253,] 1.0000000 3.792069e-87 1.896034e-87
[254,] 1.0000000 2.175337e-86 1.087668e-86
[255,] 1.0000000 1.242010e-85 6.210051e-86
[256,] 1.0000000 7.057753e-85 3.528877e-85
[257,] 1.0000000 3.991577e-84 1.995788e-84
[258,] 1.0000000 1.940624e-83 9.703122e-84
[259,] 1.0000000 1.018520e-82 5.092598e-83
[260,] 1.0000000 5.682157e-82 2.841079e-82
[261,] 1.0000000 3.154782e-81 1.577391e-81
[262,] 1.0000000 1.743136e-80 8.715682e-81
[263,] 1.0000000 4.189782e-80 2.094891e-80
[264,] 1.0000000 2.301701e-79 1.150850e-79
[265,] 1.0000000 1.273005e-78 6.365025e-79
[266,] 1.0000000 6.925541e-78 3.462770e-78
[267,] 1.0000000 3.749287e-77 1.874643e-77
[268,] 1.0000000 2.019798e-76 1.009899e-76
[269,] 1.0000000 1.105614e-75 5.528068e-76
[270,] 1.0000000 5.899069e-75 2.949534e-75
[271,] 1.0000000 3.159306e-74 1.579653e-74
[272,] 1.0000000 1.712204e-73 8.561021e-74
[273,] 1.0000000 9.000298e-73 4.500149e-73
[274,] 1.0000000 4.548468e-72 2.274234e-72
[275,] 1.0000000 2.425400e-71 1.212700e-71
[276,] 1.0000000 1.256916e-70 6.284579e-71
[277,] 1.0000000 6.480823e-70 3.240412e-70
[278,] 1.0000000 3.391832e-69 1.695916e-69
[279,] 1.0000000 1.731629e-68 8.658144e-69
[280,] 1.0000000 8.235681e-68 4.117841e-68
[281,] 1.0000000 4.162352e-67 2.081176e-67
[282,] 1.0000000 2.041252e-66 1.020626e-66
[283,] 1.0000000 1.018799e-65 5.093996e-66
[284,] 1.0000000 5.069690e-65 2.534845e-65
[285,] 1.0000000 2.476666e-64 1.238333e-64
[286,] 1.0000000 1.220625e-63 6.103125e-64
[287,] 1.0000000 5.943787e-63 2.971893e-63
[288,] 1.0000000 1.995089e-62 9.975446e-63
[289,] 1.0000000 9.696541e-62 4.848271e-62
[290,] 1.0000000 4.791553e-61 2.395776e-61
[291,] 1.0000000 1.981954e-60 9.909772e-61
[292,] 1.0000000 9.488062e-60 4.744031e-60
[293,] 1.0000000 4.517650e-59 2.258825e-59
[294,] 1.0000000 2.139378e-58 1.069689e-58
[295,] 1.0000000 1.028776e-57 5.143880e-58
[296,] 1.0000000 4.897574e-57 2.448787e-57
[297,] 1.0000000 2.346128e-56 1.173064e-56
[298,] 1.0000000 1.090289e-55 5.451445e-56
[299,] 1.0000000 5.024048e-55 2.512024e-55
[300,] 1.0000000 1.999610e-54 9.998050e-55
[301,] 1.0000000 9.349052e-54 4.674526e-54
[302,] 1.0000000 4.238831e-53 2.119415e-53
[303,] 1.0000000 1.966225e-52 9.831123e-53
[304,] 1.0000000 9.086901e-52 4.543450e-52
[305,] 1.0000000 4.051376e-51 2.025688e-51
[306,] 1.0000000 1.795901e-50 8.979505e-51
[307,] 1.0000000 7.537517e-50 3.768758e-50
[308,] 1.0000000 3.303520e-49 1.651760e-49
[309,] 1.0000000 1.457987e-48 7.289935e-49
[310,] 1.0000000 6.431942e-48 3.215971e-48
[311,] 1.0000000 2.407674e-47 1.203837e-47
[312,] 1.0000000 1.031654e-46 5.158269e-47
[313,] 1.0000000 4.446723e-46 2.223362e-46
[314,] 1.0000000 1.907777e-45 9.538883e-46
[315,] 1.0000000 8.030875e-45 4.015437e-45
[316,] 1.0000000 3.359996e-44 1.679998e-44
[317,] 1.0000000 1.397132e-43 6.985662e-44
[318,] 1.0000000 5.773528e-43 2.886764e-43
[319,] 1.0000000 2.371000e-42 1.185500e-42
[320,] 1.0000000 9.675861e-42 4.837930e-42
[321,] 1.0000000 3.969942e-41 1.984971e-41
[322,] 1.0000000 1.599673e-40 7.998363e-41
[323,] 1.0000000 6.404527e-40 3.202263e-40
[324,] 1.0000000 2.547601e-39 1.273800e-39
[325,] 1.0000000 1.006799e-38 5.033997e-39
[326,] 1.0000000 3.952759e-38 1.976379e-38
[327,] 1.0000000 1.541634e-37 7.708169e-38
[328,] 1.0000000 5.972596e-37 2.986298e-37
[329,] 1.0000000 2.329532e-36 1.164766e-36
[330,] 1.0000000 8.930177e-36 4.465088e-36
[331,] 1.0000000 3.393832e-35 1.696916e-35
[332,] 1.0000000 1.322450e-34 6.612248e-35
[333,] 1.0000000 4.955462e-34 2.477731e-34
[334,] 1.0000000 1.857433e-33 9.287165e-34
[335,] 1.0000000 6.967922e-33 3.483961e-33
[336,] 1.0000000 2.648756e-32 1.324378e-32
[337,] 1.0000000 9.656592e-32 4.828296e-32
[338,] 1.0000000 3.495095e-31 1.747548e-31
[339,] 1.0000000 1.255792e-30 6.278962e-31
[340,] 1.0000000 4.620686e-30 2.310343e-30
[341,] 1.0000000 1.636147e-29 8.180736e-30
[342,] 1.0000000 5.945144e-29 2.972572e-29
[343,] 1.0000000 2.074084e-28 1.037042e-28
[344,] 1.0000000 7.180569e-28 3.590284e-28
[345,] 1.0000000 2.339259e-27 1.169629e-27
[346,] 1.0000000 7.855837e-27 3.927919e-27
[347,] 1.0000000 2.712840e-26 1.356420e-26
[348,] 1.0000000 8.808039e-26 4.404019e-26
[349,] 1.0000000 2.940289e-25 1.470144e-25
[350,] 1.0000000 9.709432e-25 4.854716e-25
[351,] 1.0000000 3.179831e-24 1.589915e-24
[352,] 1.0000000 1.034631e-23 5.173156e-24
[353,] 1.0000000 3.332345e-23 1.666173e-23
[354,] 1.0000000 1.064121e-22 5.320605e-23
[355,] 1.0000000 3.368691e-22 1.684345e-22
[356,] 1.0000000 1.069586e-21 5.347929e-22
[357,] 1.0000000 3.326869e-21 1.663435e-21
[358,] 1.0000000 1.025496e-20 5.127480e-21
[359,] 1.0000000 3.245582e-20 1.622791e-20
[360,] 1.0000000 9.994554e-20 4.997277e-20
[361,] 1.0000000 2.998708e-19 1.499354e-19
[362,] 1.0000000 8.911874e-19 4.455937e-19
[363,] 1.0000000 2.623042e-18 1.311521e-18
[364,] 1.0000000 7.644992e-18 3.822496e-18
[365,] 1.0000000 2.236332e-17 1.118166e-17
[366,] 1.0000000 6.389299e-17 3.194649e-17
[367,] 1.0000000 1.761165e-16 8.805823e-17
[368,] 1.0000000 4.532034e-16 2.266017e-16
[369,] 1.0000000 1.294156e-15 6.470781e-16
[370,] 1.0000000 3.546293e-15 1.773146e-15
[371,] 1.0000000 9.610652e-15 4.805326e-15
[372,] 1.0000000 2.760638e-14 1.380319e-14
[373,] 1.0000000 7.128828e-14 3.564414e-14
[374,] 1.0000000 1.726107e-13 8.630536e-14
[375,] 1.0000000 4.475471e-13 2.237735e-13
[376,] 1.0000000 1.146299e-12 5.731497e-13
[377,] 1.0000000 3.001393e-12 1.500697e-12
[378,] 1.0000000 8.171533e-12 4.085766e-12
[379,] 1.0000000 2.022368e-11 1.011184e-11
[380,] 1.0000000 5.091401e-11 2.545701e-11
[381,] 1.0000000 1.254559e-10 6.272795e-11
[382,] 1.0000000 2.880587e-10 1.440293e-10
[383,] 1.0000000 6.897764e-10 3.448882e-10
[384,] 1.0000000 1.628101e-09 8.140504e-10
[385,] 1.0000000 3.871869e-09 1.935935e-09
[386,] 1.0000000 8.730066e-09 4.365033e-09
[387,] 1.0000000 2.097880e-08 1.048940e-08
[388,] 1.0000000 5.106215e-08 2.553108e-08
[389,] 1.0000000 8.174199e-08 4.087100e-08
[390,] 0.9999999 1.802571e-07 9.012855e-08
[391,] 0.9999998 4.186767e-07 2.093384e-07
[392,] 0.9999995 9.596505e-07 4.798252e-07
[393,] 0.9999992 1.654665e-06 8.273326e-07
[394,] 0.9999982 3.603466e-06 1.801733e-06
[395,] 0.9999965 7.082044e-06 3.541022e-06
[396,] 0.9999930 1.407388e-05 7.036939e-06
[397,] 0.9999872 2.563926e-05 1.281963e-05
[398,] 0.9999739 5.210188e-05 2.605094e-05
[399,] 0.9999449 1.102990e-04 5.514949e-05
[400,] 0.9998871 2.258447e-04 1.129223e-04
[401,] 0.9998032 3.935247e-04 1.967623e-04
[402,] 0.9996429 7.141098e-04 3.570549e-04
[403,] 0.9993487 1.302538e-03 6.512691e-04
[404,] 0.9988154 2.369239e-03 1.184619e-03
[405,] 0.9977313 4.537322e-03 2.268661e-03
[406,] 0.9957575 8.484940e-03 4.242470e-03
[407,] 0.9942121 1.157583e-02 5.787914e-03
[408,] 0.9897312 2.053756e-02 1.026878e-02
[409,] 0.9918348 1.633031e-02 8.165155e-03
[410,] 0.9851772 2.964564e-02 1.482282e-02
[411,] 0.9735342 5.293160e-02 2.646580e-02
[412,] 0.9773280 4.534402e-02 2.267201e-02
[413,] 0.9638779 7.224420e-02 3.612210e-02
[414,] 0.9416571 1.166859e-01 5.834294e-02
[415,] 0.9137935 1.724131e-01 8.620654e-02
[416,] 0.8869231 2.261538e-01 1.130769e-01
[417,] 0.8118077 3.763847e-01 1.881923e-01
[418,] 0.8517170 2.965659e-01 1.482830e-01
[419,] 0.8304385 3.391229e-01 1.695615e-01
[420,] 0.7214623 5.570754e-01 2.785377e-01
> postscript(file="/var/www/html/rcomp/tmp/17eoc1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/20o6x1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/30o6x1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/40o6x1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5axnh1292875566.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 = 431
Frequency = 1
1 2 3 4 5
1422177.7855 2828254.5535 1399259.2669 -2448702.6109 -378022.7100
6 7 8 9 10
212475.4498 -104017.1613 325911.2340 -403242.7649 652298.8444
11 12 13 14 15
105424.8794 281145.9654 235376.9863 -1238522.6120 -243356.9332
16 17 18 19 20
576948.3241 82206.4518 -233698.3394 -399881.1519 -571503.5540
21 22 23 24 25
209946.7167 -131539.6641 9026.1537 -638528.9441 -234635.5535
26 27 28 29 30
-101546.3212 -520050.0461 -414492.9031 286243.2595 -323799.7397
31 32 33 34 35
312171.4007 -43536.2599 267807.6146 175433.7361 310732.0000
36 37 38 39 40
-569078.4006 49814.5402 208235.4649 122385.7917 -319960.9457
41 42 43 44 45
-418353.9839 -295494.5996 18623.3084 -153032.9388 -33393.8921
46 47 48 49 50
-4999.5642 -68675.5518 -194384.4411 17765.0640 -177775.7603
51 52 53 54 55
76058.4571 -17208.1270 -62407.1809 151479.5222 -284509.5798
56 57 58 59 60
-14924.3743 293844.3840 -1201.9415 -22017.1224 -195116.1026
61 62 63 64 65
46919.6024 -197656.7454 44537.6686 -169630.2312 -42399.5980
66 67 68 69 70
-524936.9858 -46853.9476 26048.6923 -19831.6724 -42390.0984
71 72 73 74 75
-197150.8329 -51774.5780 -35135.1327 32526.8323 38018.9463
76 77 78 79 80
-86784.1679 39064.2674 -64620.0625 -52160.0646 39018.8929
81 82 83 84 85
10336.1251 -30491.9024 21645.3152 -4747.8623 77437.5254
86 87 88 89 90
-124202.2780 -19693.9857 56162.0558 -64753.8158 86266.4056
91 92 93 94 95
33704.0154 -6989.6710 46043.9995 -84723.7008 -7575.3138
96 97 98 99 100
1729.0119 142058.1967 1785.5551 53560.0816 6579.9276
101 102 103 104 105
-29502.7850 -80171.6983 -114712.4397 21671.5858 25745.8450
106 107 108 109 110
-53544.7222 -275735.4076 176036.5982 -132744.2026 -87592.5503
111 112 113 114 115
90557.5062 -228631.1611 -146148.2854 -49099.8440 70476.2658
116 117 118 119 120
-4009.2257 44140.8107 102558.5655 -118258.5694 816.2494
121 122 123 124 125
-243987.8612 30189.5715 -92740.5446 142401.7259 3826.0383
126 127 128 129 130
-153906.6564 45876.2528 30075.6448 9665.3884 45876.2528
131 132 133 134 135
35076.0683 33673.0231 45876.2528 16601.0575 23327.4580
136 137 138 139 140
33445.0498 33987.1176 43983.9643 20663.1323 45876.2528
141 142 143 144 145
37187.6608 45876.2528 -7005.3106 40383.7450 30884.5952
146 147 148 149 150
25016.3156 45876.2528 36255.9339 47247.8386 46067.5680
151 152 153 154 155
1761.4402 -61450.0008 39425.6762 9230.1870 42724.3296
156 157 158 159 160
45876.2528 45876.2528 41012.2288 51627.6528 45876.2528
161 162 163 164 165
45876.2528 45876.2528 45876.2528 17187.3214 -8192.4349
166 167 168 169 170
43061.9100 45876.2528 73457.3654 -100530.2721 45876.2528
171 172 173 174 175
45876.2528 38251.8603 45876.2528 24048.7859 45876.2528
176 177 178 179 180
32971.1674 -110753.6285 45121.9943 45876.2528 45876.2528
181 182 183 184 185
44095.8056 34778.7980 -41202.6439 -986.0808 38349.4798
186 187 188 189 190
35584.4931 2711.9615 32834.2016 10898.0508 45876.2528
191 192 193 194 195
34205.9235 45876.2528 33752.7633 45876.2528 41702.5373
196 197 198 199 200
45876.2528 -83418.4463 45876.2528 45876.2528 -34734.8816
201 202 203 204 205
45876.2528 45876.2528 45876.2528 11544.4399 29663.8946
206 207 208 209 210
39982.3299 45114.4407 45876.2528 15762.1198 27787.7310
211 212 213 214 215
45876.2528 45876.2528 40534.1392 -3769.5130 16870.9861
216 217 218 219 220
43296.0221 38116.0415 30533.0992 43296.0221 -24039.3876
221 222 223 224 225
16993.5346 45876.2528 44494.6336 6636.3416 44345.4449
226 227 228 229 230
10965.6898 32104.3027 45876.2528 36062.7740 45876.2528
231 232 233 234 235
13243.8692 24147.2612 45876.2528 8919.6114 45876.2528
236 237 238 239 240
45939.5835 45876.2528 37391.3257 45876.2528 45876.2528
241 242 243 244 245
45876.2528 37460.4520 45876.2528 36578.4145 1382.6751
246 247 248 249 250
45876.2528 45876.2528 17273.1870 45876.2528 45876.2528
251 252 253 254 255
11089.2820 45876.2528 32881.7550 -27733.2665 -140283.3289
256 257 258 259 260
45876.2528 23055.8072 -12528.5493 45876.2528 45876.2528
261 262 263 264 265
45876.2528 45876.2528 52719.1282 42846.8695 45876.2528
266 267 268 269 270
45876.2528 45876.2528 95960.3053 45876.2528 25360.0761
271 272 273 274 275
45876.2528 45876.2528 45876.2528 13455.7335 45876.2528
276 277 278 279 280
28873.5045 13907.0008 45876.2528 -43985.1358 27002.4361
281 282 283 284 285
45876.2528 45876.2528 5007.4327 45876.2528 42937.2846
286 287 288 289 290
45876.2528 30133.0662 44068.8680 45876.2528 1041.9318
291 292 293 294 295
45876.2528 36588.1557 62051.9763 45876.2528 18846.8644
296 297 298 299 300
51512.3879 45876.2528 45876.2528 45876.2528 34265.2149
301 302 303 304 305
37844.7527 24038.0960 44586.1374 45876.2528 48868.1653
306 307 308 309 310
30365.1997 45876.2528 18067.6250 10015.7961 45876.2528
311 312 313 314 315
45876.2528 27853.9242 45876.2528 33395.7416 35765.8263
316 317 318 319 320
61876.7863 44761.5447 -6200.1958 39395.8415 45876.2528
321 322 323 324 325
45876.2528 45876.2528 45876.2528 45876.2528 45876.2528
326 327 328 329 330
41227.7757 45876.2528 45876.2528 45876.2528 45876.2528
331 332 333 334 335
45876.2528 45876.2528 45876.2528 36968.6674 44572.3028
336 337 338 339 340
33850.7190 -17595.4776 45876.2528 42518.4372 32563.7714
341 342 343 344 345
3957.1427 45876.2528 45876.2528 45876.2528 8421.4511
346 347 348 349 350
45876.2528 27704.5864 45876.2528 45876.2528 2137.8717
351 352 353 354 355
38503.5036 34204.3680 30477.1743 40468.4630 45876.2528
356 357 358 359 360
45876.2528 22299.5254 45876.2528 45876.2528 45876.2528
361 362 363 364 365
40711.9569 45876.2528 45876.2528 -35785.4212 31340.9823
366 367 368 369 370
45876.2528 45876.2528 45876.2528 45876.2528 32258.4971
371 372 373 374 375
45876.2528 35044.5333 53147.9855 30250.6819 45876.2528
376 377 378 379 380
45876.2528 9032.9946 30551.8278 -10839.0923 45876.2528
381 382 383 384 385
45876.2528 12992.8254 -19634.9329 45876.2528 30641.7681
386 387 388 389 390
-123138.4747 50836.5709 11640.2014 39263.8138 15615.9244
391 392 393 394 395
31253.7532 5065.9793 -30129.8069 75257.3307 -71988.3779
396 397 398 399 400
-8508.7294 -166757.5998 -170833.8386 6175.9286 -138197.8894
401 402 403 404 405
7711.3690 36736.6607 6333.6046 -50925.8076 -14562.6528
406 407 408 409 410
17503.4568 19791.6817 -2024.5758 7417.6489 -57365.0314
411 412 413 414 415
-62068.8376 -171624.2565 -44517.9813 65182.9517 -160213.4970
416 417 418 419 420
-88829.3916 36889.0594 -29612.3537 -52885.4894 -38685.5501
421 422 423 424 425
-39444.3383 -186172.9561 -692206.1533 -479941.8205 -72026.6461
426 427 428 429 430
-271484.3057 -341158.8368 -1055853.5887 -259264.3938 -340495.0209
431
-760308.8934
> postscript(file="/var/www/html/rcomp/tmp/6axnh1292875566.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 = 431
Frequency = 1
lag(myerror, k = 1) myerror
0 1422177.7855 NA
1 2828254.5535 1422177.7855
2 1399259.2669 2828254.5535
3 -2448702.6109 1399259.2669
4 -378022.7100 -2448702.6109
5 212475.4498 -378022.7100
6 -104017.1613 212475.4498
7 325911.2340 -104017.1613
8 -403242.7649 325911.2340
9 652298.8444 -403242.7649
10 105424.8794 652298.8444
11 281145.9654 105424.8794
12 235376.9863 281145.9654
13 -1238522.6120 235376.9863
14 -243356.9332 -1238522.6120
15 576948.3241 -243356.9332
16 82206.4518 576948.3241
17 -233698.3394 82206.4518
18 -399881.1519 -233698.3394
19 -571503.5540 -399881.1519
20 209946.7167 -571503.5540
21 -131539.6641 209946.7167
22 9026.1537 -131539.6641
23 -638528.9441 9026.1537
24 -234635.5535 -638528.9441
25 -101546.3212 -234635.5535
26 -520050.0461 -101546.3212
27 -414492.9031 -520050.0461
28 286243.2595 -414492.9031
29 -323799.7397 286243.2595
30 312171.4007 -323799.7397
31 -43536.2599 312171.4007
32 267807.6146 -43536.2599
33 175433.7361 267807.6146
34 310732.0000 175433.7361
35 -569078.4006 310732.0000
36 49814.5402 -569078.4006
37 208235.4649 49814.5402
38 122385.7917 208235.4649
39 -319960.9457 122385.7917
40 -418353.9839 -319960.9457
41 -295494.5996 -418353.9839
42 18623.3084 -295494.5996
43 -153032.9388 18623.3084
44 -33393.8921 -153032.9388
45 -4999.5642 -33393.8921
46 -68675.5518 -4999.5642
47 -194384.4411 -68675.5518
48 17765.0640 -194384.4411
49 -177775.7603 17765.0640
50 76058.4571 -177775.7603
51 -17208.1270 76058.4571
52 -62407.1809 -17208.1270
53 151479.5222 -62407.1809
54 -284509.5798 151479.5222
55 -14924.3743 -284509.5798
56 293844.3840 -14924.3743
57 -1201.9415 293844.3840
58 -22017.1224 -1201.9415
59 -195116.1026 -22017.1224
60 46919.6024 -195116.1026
61 -197656.7454 46919.6024
62 44537.6686 -197656.7454
63 -169630.2312 44537.6686
64 -42399.5980 -169630.2312
65 -524936.9858 -42399.5980
66 -46853.9476 -524936.9858
67 26048.6923 -46853.9476
68 -19831.6724 26048.6923
69 -42390.0984 -19831.6724
70 -197150.8329 -42390.0984
71 -51774.5780 -197150.8329
72 -35135.1327 -51774.5780
73 32526.8323 -35135.1327
74 38018.9463 32526.8323
75 -86784.1679 38018.9463
76 39064.2674 -86784.1679
77 -64620.0625 39064.2674
78 -52160.0646 -64620.0625
79 39018.8929 -52160.0646
80 10336.1251 39018.8929
81 -30491.9024 10336.1251
82 21645.3152 -30491.9024
83 -4747.8623 21645.3152
84 77437.5254 -4747.8623
85 -124202.2780 77437.5254
86 -19693.9857 -124202.2780
87 56162.0558 -19693.9857
88 -64753.8158 56162.0558
89 86266.4056 -64753.8158
90 33704.0154 86266.4056
91 -6989.6710 33704.0154
92 46043.9995 -6989.6710
93 -84723.7008 46043.9995
94 -7575.3138 -84723.7008
95 1729.0119 -7575.3138
96 142058.1967 1729.0119
97 1785.5551 142058.1967
98 53560.0816 1785.5551
99 6579.9276 53560.0816
100 -29502.7850 6579.9276
101 -80171.6983 -29502.7850
102 -114712.4397 -80171.6983
103 21671.5858 -114712.4397
104 25745.8450 21671.5858
105 -53544.7222 25745.8450
106 -275735.4076 -53544.7222
107 176036.5982 -275735.4076
108 -132744.2026 176036.5982
109 -87592.5503 -132744.2026
110 90557.5062 -87592.5503
111 -228631.1611 90557.5062
112 -146148.2854 -228631.1611
113 -49099.8440 -146148.2854
114 70476.2658 -49099.8440
115 -4009.2257 70476.2658
116 44140.8107 -4009.2257
117 102558.5655 44140.8107
118 -118258.5694 102558.5655
119 816.2494 -118258.5694
120 -243987.8612 816.2494
121 30189.5715 -243987.8612
122 -92740.5446 30189.5715
123 142401.7259 -92740.5446
124 3826.0383 142401.7259
125 -153906.6564 3826.0383
126 45876.2528 -153906.6564
127 30075.6448 45876.2528
128 9665.3884 30075.6448
129 45876.2528 9665.3884
130 35076.0683 45876.2528
131 33673.0231 35076.0683
132 45876.2528 33673.0231
133 16601.0575 45876.2528
134 23327.4580 16601.0575
135 33445.0498 23327.4580
136 33987.1176 33445.0498
137 43983.9643 33987.1176
138 20663.1323 43983.9643
139 45876.2528 20663.1323
140 37187.6608 45876.2528
141 45876.2528 37187.6608
142 -7005.3106 45876.2528
143 40383.7450 -7005.3106
144 30884.5952 40383.7450
145 25016.3156 30884.5952
146 45876.2528 25016.3156
147 36255.9339 45876.2528
148 47247.8386 36255.9339
149 46067.5680 47247.8386
150 1761.4402 46067.5680
151 -61450.0008 1761.4402
152 39425.6762 -61450.0008
153 9230.1870 39425.6762
154 42724.3296 9230.1870
155 45876.2528 42724.3296
156 45876.2528 45876.2528
157 41012.2288 45876.2528
158 51627.6528 41012.2288
159 45876.2528 51627.6528
160 45876.2528 45876.2528
161 45876.2528 45876.2528
162 45876.2528 45876.2528
163 17187.3214 45876.2528
164 -8192.4349 17187.3214
165 43061.9100 -8192.4349
166 45876.2528 43061.9100
167 73457.3654 45876.2528
168 -100530.2721 73457.3654
169 45876.2528 -100530.2721
170 45876.2528 45876.2528
171 38251.8603 45876.2528
172 45876.2528 38251.8603
173 24048.7859 45876.2528
174 45876.2528 24048.7859
175 32971.1674 45876.2528
176 -110753.6285 32971.1674
177 45121.9943 -110753.6285
178 45876.2528 45121.9943
179 45876.2528 45876.2528
180 44095.8056 45876.2528
181 34778.7980 44095.8056
182 -41202.6439 34778.7980
183 -986.0808 -41202.6439
184 38349.4798 -986.0808
185 35584.4931 38349.4798
186 2711.9615 35584.4931
187 32834.2016 2711.9615
188 10898.0508 32834.2016
189 45876.2528 10898.0508
190 34205.9235 45876.2528
191 45876.2528 34205.9235
192 33752.7633 45876.2528
193 45876.2528 33752.7633
194 41702.5373 45876.2528
195 45876.2528 41702.5373
196 -83418.4463 45876.2528
197 45876.2528 -83418.4463
198 45876.2528 45876.2528
199 -34734.8816 45876.2528
200 45876.2528 -34734.8816
201 45876.2528 45876.2528
202 45876.2528 45876.2528
203 11544.4399 45876.2528
204 29663.8946 11544.4399
205 39982.3299 29663.8946
206 45114.4407 39982.3299
207 45876.2528 45114.4407
208 15762.1198 45876.2528
209 27787.7310 15762.1198
210 45876.2528 27787.7310
211 45876.2528 45876.2528
212 40534.1392 45876.2528
213 -3769.5130 40534.1392
214 16870.9861 -3769.5130
215 43296.0221 16870.9861
216 38116.0415 43296.0221
217 30533.0992 38116.0415
218 43296.0221 30533.0992
219 -24039.3876 43296.0221
220 16993.5346 -24039.3876
221 45876.2528 16993.5346
222 44494.6336 45876.2528
223 6636.3416 44494.6336
224 44345.4449 6636.3416
225 10965.6898 44345.4449
226 32104.3027 10965.6898
227 45876.2528 32104.3027
228 36062.7740 45876.2528
229 45876.2528 36062.7740
230 13243.8692 45876.2528
231 24147.2612 13243.8692
232 45876.2528 24147.2612
233 8919.6114 45876.2528
234 45876.2528 8919.6114
235 45939.5835 45876.2528
236 45876.2528 45939.5835
237 37391.3257 45876.2528
238 45876.2528 37391.3257
239 45876.2528 45876.2528
240 45876.2528 45876.2528
241 37460.4520 45876.2528
242 45876.2528 37460.4520
243 36578.4145 45876.2528
244 1382.6751 36578.4145
245 45876.2528 1382.6751
246 45876.2528 45876.2528
247 17273.1870 45876.2528
248 45876.2528 17273.1870
249 45876.2528 45876.2528
250 11089.2820 45876.2528
251 45876.2528 11089.2820
252 32881.7550 45876.2528
253 -27733.2665 32881.7550
254 -140283.3289 -27733.2665
255 45876.2528 -140283.3289
256 23055.8072 45876.2528
257 -12528.5493 23055.8072
258 45876.2528 -12528.5493
259 45876.2528 45876.2528
260 45876.2528 45876.2528
261 45876.2528 45876.2528
262 52719.1282 45876.2528
263 42846.8695 52719.1282
264 45876.2528 42846.8695
265 45876.2528 45876.2528
266 45876.2528 45876.2528
267 95960.3053 45876.2528
268 45876.2528 95960.3053
269 25360.0761 45876.2528
270 45876.2528 25360.0761
271 45876.2528 45876.2528
272 45876.2528 45876.2528
273 13455.7335 45876.2528
274 45876.2528 13455.7335
275 28873.5045 45876.2528
276 13907.0008 28873.5045
277 45876.2528 13907.0008
278 -43985.1358 45876.2528
279 27002.4361 -43985.1358
280 45876.2528 27002.4361
281 45876.2528 45876.2528
282 5007.4327 45876.2528
283 45876.2528 5007.4327
284 42937.2846 45876.2528
285 45876.2528 42937.2846
286 30133.0662 45876.2528
287 44068.8680 30133.0662
288 45876.2528 44068.8680
289 1041.9318 45876.2528
290 45876.2528 1041.9318
291 36588.1557 45876.2528
292 62051.9763 36588.1557
293 45876.2528 62051.9763
294 18846.8644 45876.2528
295 51512.3879 18846.8644
296 45876.2528 51512.3879
297 45876.2528 45876.2528
298 45876.2528 45876.2528
299 34265.2149 45876.2528
300 37844.7527 34265.2149
301 24038.0960 37844.7527
302 44586.1374 24038.0960
303 45876.2528 44586.1374
304 48868.1653 45876.2528
305 30365.1997 48868.1653
306 45876.2528 30365.1997
307 18067.6250 45876.2528
308 10015.7961 18067.6250
309 45876.2528 10015.7961
310 45876.2528 45876.2528
311 27853.9242 45876.2528
312 45876.2528 27853.9242
313 33395.7416 45876.2528
314 35765.8263 33395.7416
315 61876.7863 35765.8263
316 44761.5447 61876.7863
317 -6200.1958 44761.5447
318 39395.8415 -6200.1958
319 45876.2528 39395.8415
320 45876.2528 45876.2528
321 45876.2528 45876.2528
322 45876.2528 45876.2528
323 45876.2528 45876.2528
324 45876.2528 45876.2528
325 41227.7757 45876.2528
326 45876.2528 41227.7757
327 45876.2528 45876.2528
328 45876.2528 45876.2528
329 45876.2528 45876.2528
330 45876.2528 45876.2528
331 45876.2528 45876.2528
332 45876.2528 45876.2528
333 36968.6674 45876.2528
334 44572.3028 36968.6674
335 33850.7190 44572.3028
336 -17595.4776 33850.7190
337 45876.2528 -17595.4776
338 42518.4372 45876.2528
339 32563.7714 42518.4372
340 3957.1427 32563.7714
341 45876.2528 3957.1427
342 45876.2528 45876.2528
343 45876.2528 45876.2528
344 8421.4511 45876.2528
345 45876.2528 8421.4511
346 27704.5864 45876.2528
347 45876.2528 27704.5864
348 45876.2528 45876.2528
349 2137.8717 45876.2528
350 38503.5036 2137.8717
351 34204.3680 38503.5036
352 30477.1743 34204.3680
353 40468.4630 30477.1743
354 45876.2528 40468.4630
355 45876.2528 45876.2528
356 22299.5254 45876.2528
357 45876.2528 22299.5254
358 45876.2528 45876.2528
359 45876.2528 45876.2528
360 40711.9569 45876.2528
361 45876.2528 40711.9569
362 45876.2528 45876.2528
363 -35785.4212 45876.2528
364 31340.9823 -35785.4212
365 45876.2528 31340.9823
366 45876.2528 45876.2528
367 45876.2528 45876.2528
368 45876.2528 45876.2528
369 32258.4971 45876.2528
370 45876.2528 32258.4971
371 35044.5333 45876.2528
372 53147.9855 35044.5333
373 30250.6819 53147.9855
374 45876.2528 30250.6819
375 45876.2528 45876.2528
376 9032.9946 45876.2528
377 30551.8278 9032.9946
378 -10839.0923 30551.8278
379 45876.2528 -10839.0923
380 45876.2528 45876.2528
381 12992.8254 45876.2528
382 -19634.9329 12992.8254
383 45876.2528 -19634.9329
384 30641.7681 45876.2528
385 -123138.4747 30641.7681
386 50836.5709 -123138.4747
387 11640.2014 50836.5709
388 39263.8138 11640.2014
389 15615.9244 39263.8138
390 31253.7532 15615.9244
391 5065.9793 31253.7532
392 -30129.8069 5065.9793
393 75257.3307 -30129.8069
394 -71988.3779 75257.3307
395 -8508.7294 -71988.3779
396 -166757.5998 -8508.7294
397 -170833.8386 -166757.5998
398 6175.9286 -170833.8386
399 -138197.8894 6175.9286
400 7711.3690 -138197.8894
401 36736.6607 7711.3690
402 6333.6046 36736.6607
403 -50925.8076 6333.6046
404 -14562.6528 -50925.8076
405 17503.4568 -14562.6528
406 19791.6817 17503.4568
407 -2024.5758 19791.6817
408 7417.6489 -2024.5758
409 -57365.0314 7417.6489
410 -62068.8376 -57365.0314
411 -171624.2565 -62068.8376
412 -44517.9813 -171624.2565
413 65182.9517 -44517.9813
414 -160213.4970 65182.9517
415 -88829.3916 -160213.4970
416 36889.0594 -88829.3916
417 -29612.3537 36889.0594
418 -52885.4894 -29612.3537
419 -38685.5501 -52885.4894
420 -39444.3383 -38685.5501
421 -186172.9561 -39444.3383
422 -692206.1533 -186172.9561
423 -479941.8205 -692206.1533
424 -72026.6461 -479941.8205
425 -271484.3057 -72026.6461
426 -341158.8368 -271484.3057
427 -1055853.5887 -341158.8368
428 -259264.3938 -1055853.5887
429 -340495.0209 -259264.3938
430 -760308.8934 -340495.0209
431 NA -760308.8934
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2828254.5535 1422177.7855
[2,] 1399259.2669 2828254.5535
[3,] -2448702.6109 1399259.2669
[4,] -378022.7100 -2448702.6109
[5,] 212475.4498 -378022.7100
[6,] -104017.1613 212475.4498
[7,] 325911.2340 -104017.1613
[8,] -403242.7649 325911.2340
[9,] 652298.8444 -403242.7649
[10,] 105424.8794 652298.8444
[11,] 281145.9654 105424.8794
[12,] 235376.9863 281145.9654
[13,] -1238522.6120 235376.9863
[14,] -243356.9332 -1238522.6120
[15,] 576948.3241 -243356.9332
[16,] 82206.4518 576948.3241
[17,] -233698.3394 82206.4518
[18,] -399881.1519 -233698.3394
[19,] -571503.5540 -399881.1519
[20,] 209946.7167 -571503.5540
[21,] -131539.6641 209946.7167
[22,] 9026.1537 -131539.6641
[23,] -638528.9441 9026.1537
[24,] -234635.5535 -638528.9441
[25,] -101546.3212 -234635.5535
[26,] -520050.0461 -101546.3212
[27,] -414492.9031 -520050.0461
[28,] 286243.2595 -414492.9031
[29,] -323799.7397 286243.2595
[30,] 312171.4007 -323799.7397
[31,] -43536.2599 312171.4007
[32,] 267807.6146 -43536.2599
[33,] 175433.7361 267807.6146
[34,] 310732.0000 175433.7361
[35,] -569078.4006 310732.0000
[36,] 49814.5402 -569078.4006
[37,] 208235.4649 49814.5402
[38,] 122385.7917 208235.4649
[39,] -319960.9457 122385.7917
[40,] -418353.9839 -319960.9457
[41,] -295494.5996 -418353.9839
[42,] 18623.3084 -295494.5996
[43,] -153032.9388 18623.3084
[44,] -33393.8921 -153032.9388
[45,] -4999.5642 -33393.8921
[46,] -68675.5518 -4999.5642
[47,] -194384.4411 -68675.5518
[48,] 17765.0640 -194384.4411
[49,] -177775.7603 17765.0640
[50,] 76058.4571 -177775.7603
[51,] -17208.1270 76058.4571
[52,] -62407.1809 -17208.1270
[53,] 151479.5222 -62407.1809
[54,] -284509.5798 151479.5222
[55,] -14924.3743 -284509.5798
[56,] 293844.3840 -14924.3743
[57,] -1201.9415 293844.3840
[58,] -22017.1224 -1201.9415
[59,] -195116.1026 -22017.1224
[60,] 46919.6024 -195116.1026
[61,] -197656.7454 46919.6024
[62,] 44537.6686 -197656.7454
[63,] -169630.2312 44537.6686
[64,] -42399.5980 -169630.2312
[65,] -524936.9858 -42399.5980
[66,] -46853.9476 -524936.9858
[67,] 26048.6923 -46853.9476
[68,] -19831.6724 26048.6923
[69,] -42390.0984 -19831.6724
[70,] -197150.8329 -42390.0984
[71,] -51774.5780 -197150.8329
[72,] -35135.1327 -51774.5780
[73,] 32526.8323 -35135.1327
[74,] 38018.9463 32526.8323
[75,] -86784.1679 38018.9463
[76,] 39064.2674 -86784.1679
[77,] -64620.0625 39064.2674
[78,] -52160.0646 -64620.0625
[79,] 39018.8929 -52160.0646
[80,] 10336.1251 39018.8929
[81,] -30491.9024 10336.1251
[82,] 21645.3152 -30491.9024
[83,] -4747.8623 21645.3152
[84,] 77437.5254 -4747.8623
[85,] -124202.2780 77437.5254
[86,] -19693.9857 -124202.2780
[87,] 56162.0558 -19693.9857
[88,] -64753.8158 56162.0558
[89,] 86266.4056 -64753.8158
[90,] 33704.0154 86266.4056
[91,] -6989.6710 33704.0154
[92,] 46043.9995 -6989.6710
[93,] -84723.7008 46043.9995
[94,] -7575.3138 -84723.7008
[95,] 1729.0119 -7575.3138
[96,] 142058.1967 1729.0119
[97,] 1785.5551 142058.1967
[98,] 53560.0816 1785.5551
[99,] 6579.9276 53560.0816
[100,] -29502.7850 6579.9276
[101,] -80171.6983 -29502.7850
[102,] -114712.4397 -80171.6983
[103,] 21671.5858 -114712.4397
[104,] 25745.8450 21671.5858
[105,] -53544.7222 25745.8450
[106,] -275735.4076 -53544.7222
[107,] 176036.5982 -275735.4076
[108,] -132744.2026 176036.5982
[109,] -87592.5503 -132744.2026
[110,] 90557.5062 -87592.5503
[111,] -228631.1611 90557.5062
[112,] -146148.2854 -228631.1611
[113,] -49099.8440 -146148.2854
[114,] 70476.2658 -49099.8440
[115,] -4009.2257 70476.2658
[116,] 44140.8107 -4009.2257
[117,] 102558.5655 44140.8107
[118,] -118258.5694 102558.5655
[119,] 816.2494 -118258.5694
[120,] -243987.8612 816.2494
[121,] 30189.5715 -243987.8612
[122,] -92740.5446 30189.5715
[123,] 142401.7259 -92740.5446
[124,] 3826.0383 142401.7259
[125,] -153906.6564 3826.0383
[126,] 45876.2528 -153906.6564
[127,] 30075.6448 45876.2528
[128,] 9665.3884 30075.6448
[129,] 45876.2528 9665.3884
[130,] 35076.0683 45876.2528
[131,] 33673.0231 35076.0683
[132,] 45876.2528 33673.0231
[133,] 16601.0575 45876.2528
[134,] 23327.4580 16601.0575
[135,] 33445.0498 23327.4580
[136,] 33987.1176 33445.0498
[137,] 43983.9643 33987.1176
[138,] 20663.1323 43983.9643
[139,] 45876.2528 20663.1323
[140,] 37187.6608 45876.2528
[141,] 45876.2528 37187.6608
[142,] -7005.3106 45876.2528
[143,] 40383.7450 -7005.3106
[144,] 30884.5952 40383.7450
[145,] 25016.3156 30884.5952
[146,] 45876.2528 25016.3156
[147,] 36255.9339 45876.2528
[148,] 47247.8386 36255.9339
[149,] 46067.5680 47247.8386
[150,] 1761.4402 46067.5680
[151,] -61450.0008 1761.4402
[152,] 39425.6762 -61450.0008
[153,] 9230.1870 39425.6762
[154,] 42724.3296 9230.1870
[155,] 45876.2528 42724.3296
[156,] 45876.2528 45876.2528
[157,] 41012.2288 45876.2528
[158,] 51627.6528 41012.2288
[159,] 45876.2528 51627.6528
[160,] 45876.2528 45876.2528
[161,] 45876.2528 45876.2528
[162,] 45876.2528 45876.2528
[163,] 17187.3214 45876.2528
[164,] -8192.4349 17187.3214
[165,] 43061.9100 -8192.4349
[166,] 45876.2528 43061.9100
[167,] 73457.3654 45876.2528
[168,] -100530.2721 73457.3654
[169,] 45876.2528 -100530.2721
[170,] 45876.2528 45876.2528
[171,] 38251.8603 45876.2528
[172,] 45876.2528 38251.8603
[173,] 24048.7859 45876.2528
[174,] 45876.2528 24048.7859
[175,] 32971.1674 45876.2528
[176,] -110753.6285 32971.1674
[177,] 45121.9943 -110753.6285
[178,] 45876.2528 45121.9943
[179,] 45876.2528 45876.2528
[180,] 44095.8056 45876.2528
[181,] 34778.7980 44095.8056
[182,] -41202.6439 34778.7980
[183,] -986.0808 -41202.6439
[184,] 38349.4798 -986.0808
[185,] 35584.4931 38349.4798
[186,] 2711.9615 35584.4931
[187,] 32834.2016 2711.9615
[188,] 10898.0508 32834.2016
[189,] 45876.2528 10898.0508
[190,] 34205.9235 45876.2528
[191,] 45876.2528 34205.9235
[192,] 33752.7633 45876.2528
[193,] 45876.2528 33752.7633
[194,] 41702.5373 45876.2528
[195,] 45876.2528 41702.5373
[196,] -83418.4463 45876.2528
[197,] 45876.2528 -83418.4463
[198,] 45876.2528 45876.2528
[199,] -34734.8816 45876.2528
[200,] 45876.2528 -34734.8816
[201,] 45876.2528 45876.2528
[202,] 45876.2528 45876.2528
[203,] 11544.4399 45876.2528
[204,] 29663.8946 11544.4399
[205,] 39982.3299 29663.8946
[206,] 45114.4407 39982.3299
[207,] 45876.2528 45114.4407
[208,] 15762.1198 45876.2528
[209,] 27787.7310 15762.1198
[210,] 45876.2528 27787.7310
[211,] 45876.2528 45876.2528
[212,] 40534.1392 45876.2528
[213,] -3769.5130 40534.1392
[214,] 16870.9861 -3769.5130
[215,] 43296.0221 16870.9861
[216,] 38116.0415 43296.0221
[217,] 30533.0992 38116.0415
[218,] 43296.0221 30533.0992
[219,] -24039.3876 43296.0221
[220,] 16993.5346 -24039.3876
[221,] 45876.2528 16993.5346
[222,] 44494.6336 45876.2528
[223,] 6636.3416 44494.6336
[224,] 44345.4449 6636.3416
[225,] 10965.6898 44345.4449
[226,] 32104.3027 10965.6898
[227,] 45876.2528 32104.3027
[228,] 36062.7740 45876.2528
[229,] 45876.2528 36062.7740
[230,] 13243.8692 45876.2528
[231,] 24147.2612 13243.8692
[232,] 45876.2528 24147.2612
[233,] 8919.6114 45876.2528
[234,] 45876.2528 8919.6114
[235,] 45939.5835 45876.2528
[236,] 45876.2528 45939.5835
[237,] 37391.3257 45876.2528
[238,] 45876.2528 37391.3257
[239,] 45876.2528 45876.2528
[240,] 45876.2528 45876.2528
[241,] 37460.4520 45876.2528
[242,] 45876.2528 37460.4520
[243,] 36578.4145 45876.2528
[244,] 1382.6751 36578.4145
[245,] 45876.2528 1382.6751
[246,] 45876.2528 45876.2528
[247,] 17273.1870 45876.2528
[248,] 45876.2528 17273.1870
[249,] 45876.2528 45876.2528
[250,] 11089.2820 45876.2528
[251,] 45876.2528 11089.2820
[252,] 32881.7550 45876.2528
[253,] -27733.2665 32881.7550
[254,] -140283.3289 -27733.2665
[255,] 45876.2528 -140283.3289
[256,] 23055.8072 45876.2528
[257,] -12528.5493 23055.8072
[258,] 45876.2528 -12528.5493
[259,] 45876.2528 45876.2528
[260,] 45876.2528 45876.2528
[261,] 45876.2528 45876.2528
[262,] 52719.1282 45876.2528
[263,] 42846.8695 52719.1282
[264,] 45876.2528 42846.8695
[265,] 45876.2528 45876.2528
[266,] 45876.2528 45876.2528
[267,] 95960.3053 45876.2528
[268,] 45876.2528 95960.3053
[269,] 25360.0761 45876.2528
[270,] 45876.2528 25360.0761
[271,] 45876.2528 45876.2528
[272,] 45876.2528 45876.2528
[273,] 13455.7335 45876.2528
[274,] 45876.2528 13455.7335
[275,] 28873.5045 45876.2528
[276,] 13907.0008 28873.5045
[277,] 45876.2528 13907.0008
[278,] -43985.1358 45876.2528
[279,] 27002.4361 -43985.1358
[280,] 45876.2528 27002.4361
[281,] 45876.2528 45876.2528
[282,] 5007.4327 45876.2528
[283,] 45876.2528 5007.4327
[284,] 42937.2846 45876.2528
[285,] 45876.2528 42937.2846
[286,] 30133.0662 45876.2528
[287,] 44068.8680 30133.0662
[288,] 45876.2528 44068.8680
[289,] 1041.9318 45876.2528
[290,] 45876.2528 1041.9318
[291,] 36588.1557 45876.2528
[292,] 62051.9763 36588.1557
[293,] 45876.2528 62051.9763
[294,] 18846.8644 45876.2528
[295,] 51512.3879 18846.8644
[296,] 45876.2528 51512.3879
[297,] 45876.2528 45876.2528
[298,] 45876.2528 45876.2528
[299,] 34265.2149 45876.2528
[300,] 37844.7527 34265.2149
[301,] 24038.0960 37844.7527
[302,] 44586.1374 24038.0960
[303,] 45876.2528 44586.1374
[304,] 48868.1653 45876.2528
[305,] 30365.1997 48868.1653
[306,] 45876.2528 30365.1997
[307,] 18067.6250 45876.2528
[308,] 10015.7961 18067.6250
[309,] 45876.2528 10015.7961
[310,] 45876.2528 45876.2528
[311,] 27853.9242 45876.2528
[312,] 45876.2528 27853.9242
[313,] 33395.7416 45876.2528
[314,] 35765.8263 33395.7416
[315,] 61876.7863 35765.8263
[316,] 44761.5447 61876.7863
[317,] -6200.1958 44761.5447
[318,] 39395.8415 -6200.1958
[319,] 45876.2528 39395.8415
[320,] 45876.2528 45876.2528
[321,] 45876.2528 45876.2528
[322,] 45876.2528 45876.2528
[323,] 45876.2528 45876.2528
[324,] 45876.2528 45876.2528
[325,] 41227.7757 45876.2528
[326,] 45876.2528 41227.7757
[327,] 45876.2528 45876.2528
[328,] 45876.2528 45876.2528
[329,] 45876.2528 45876.2528
[330,] 45876.2528 45876.2528
[331,] 45876.2528 45876.2528
[332,] 45876.2528 45876.2528
[333,] 36968.6674 45876.2528
[334,] 44572.3028 36968.6674
[335,] 33850.7190 44572.3028
[336,] -17595.4776 33850.7190
[337,] 45876.2528 -17595.4776
[338,] 42518.4372 45876.2528
[339,] 32563.7714 42518.4372
[340,] 3957.1427 32563.7714
[341,] 45876.2528 3957.1427
[342,] 45876.2528 45876.2528
[343,] 45876.2528 45876.2528
[344,] 8421.4511 45876.2528
[345,] 45876.2528 8421.4511
[346,] 27704.5864 45876.2528
[347,] 45876.2528 27704.5864
[348,] 45876.2528 45876.2528
[349,] 2137.8717 45876.2528
[350,] 38503.5036 2137.8717
[351,] 34204.3680 38503.5036
[352,] 30477.1743 34204.3680
[353,] 40468.4630 30477.1743
[354,] 45876.2528 40468.4630
[355,] 45876.2528 45876.2528
[356,] 22299.5254 45876.2528
[357,] 45876.2528 22299.5254
[358,] 45876.2528 45876.2528
[359,] 45876.2528 45876.2528
[360,] 40711.9569 45876.2528
[361,] 45876.2528 40711.9569
[362,] 45876.2528 45876.2528
[363,] -35785.4212 45876.2528
[364,] 31340.9823 -35785.4212
[365,] 45876.2528 31340.9823
[366,] 45876.2528 45876.2528
[367,] 45876.2528 45876.2528
[368,] 45876.2528 45876.2528
[369,] 32258.4971 45876.2528
[370,] 45876.2528 32258.4971
[371,] 35044.5333 45876.2528
[372,] 53147.9855 35044.5333
[373,] 30250.6819 53147.9855
[374,] 45876.2528 30250.6819
[375,] 45876.2528 45876.2528
[376,] 9032.9946 45876.2528
[377,] 30551.8278 9032.9946
[378,] -10839.0923 30551.8278
[379,] 45876.2528 -10839.0923
[380,] 45876.2528 45876.2528
[381,] 12992.8254 45876.2528
[382,] -19634.9329 12992.8254
[383,] 45876.2528 -19634.9329
[384,] 30641.7681 45876.2528
[385,] -123138.4747 30641.7681
[386,] 50836.5709 -123138.4747
[387,] 11640.2014 50836.5709
[388,] 39263.8138 11640.2014
[389,] 15615.9244 39263.8138
[390,] 31253.7532 15615.9244
[391,] 5065.9793 31253.7532
[392,] -30129.8069 5065.9793
[393,] 75257.3307 -30129.8069
[394,] -71988.3779 75257.3307
[395,] -8508.7294 -71988.3779
[396,] -166757.5998 -8508.7294
[397,] -170833.8386 -166757.5998
[398,] 6175.9286 -170833.8386
[399,] -138197.8894 6175.9286
[400,] 7711.3690 -138197.8894
[401,] 36736.6607 7711.3690
[402,] 6333.6046 36736.6607
[403,] -50925.8076 6333.6046
[404,] -14562.6528 -50925.8076
[405,] 17503.4568 -14562.6528
[406,] 19791.6817 17503.4568
[407,] -2024.5758 19791.6817
[408,] 7417.6489 -2024.5758
[409,] -57365.0314 7417.6489
[410,] -62068.8376 -57365.0314
[411,] -171624.2565 -62068.8376
[412,] -44517.9813 -171624.2565
[413,] 65182.9517 -44517.9813
[414,] -160213.4970 65182.9517
[415,] -88829.3916 -160213.4970
[416,] 36889.0594 -88829.3916
[417,] -29612.3537 36889.0594
[418,] -52885.4894 -29612.3537
[419,] -38685.5501 -52885.4894
[420,] -39444.3383 -38685.5501
[421,] -186172.9561 -39444.3383
[422,] -692206.1533 -186172.9561
[423,] -479941.8205 -692206.1533
[424,] -72026.6461 -479941.8205
[425,] -271484.3057 -72026.6461
[426,] -341158.8368 -271484.3057
[427,] -1055853.5887 -341158.8368
[428,] -259264.3938 -1055853.5887
[429,] -340495.0209 -259264.3938
[430,] -760308.8934 -340495.0209
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2828254.5535 1422177.7855
2 1399259.2669 2828254.5535
3 -2448702.6109 1399259.2669
4 -378022.7100 -2448702.6109
5 212475.4498 -378022.7100
6 -104017.1613 212475.4498
7 325911.2340 -104017.1613
8 -403242.7649 325911.2340
9 652298.8444 -403242.7649
10 105424.8794 652298.8444
11 281145.9654 105424.8794
12 235376.9863 281145.9654
13 -1238522.6120 235376.9863
14 -243356.9332 -1238522.6120
15 576948.3241 -243356.9332
16 82206.4518 576948.3241
17 -233698.3394 82206.4518
18 -399881.1519 -233698.3394
19 -571503.5540 -399881.1519
20 209946.7167 -571503.5540
21 -131539.6641 209946.7167
22 9026.1537 -131539.6641
23 -638528.9441 9026.1537
24 -234635.5535 -638528.9441
25 -101546.3212 -234635.5535
26 -520050.0461 -101546.3212
27 -414492.9031 -520050.0461
28 286243.2595 -414492.9031
29 -323799.7397 286243.2595
30 312171.4007 -323799.7397
31 -43536.2599 312171.4007
32 267807.6146 -43536.2599
33 175433.7361 267807.6146
34 310732.0000 175433.7361
35 -569078.4006 310732.0000
36 49814.5402 -569078.4006
37 208235.4649 49814.5402
38 122385.7917 208235.4649
39 -319960.9457 122385.7917
40 -418353.9839 -319960.9457
41 -295494.5996 -418353.9839
42 18623.3084 -295494.5996
43 -153032.9388 18623.3084
44 -33393.8921 -153032.9388
45 -4999.5642 -33393.8921
46 -68675.5518 -4999.5642
47 -194384.4411 -68675.5518
48 17765.0640 -194384.4411
49 -177775.7603 17765.0640
50 76058.4571 -177775.7603
51 -17208.1270 76058.4571
52 -62407.1809 -17208.1270
53 151479.5222 -62407.1809
54 -284509.5798 151479.5222
55 -14924.3743 -284509.5798
56 293844.3840 -14924.3743
57 -1201.9415 293844.3840
58 -22017.1224 -1201.9415
59 -195116.1026 -22017.1224
60 46919.6024 -195116.1026
61 -197656.7454 46919.6024
62 44537.6686 -197656.7454
63 -169630.2312 44537.6686
64 -42399.5980 -169630.2312
65 -524936.9858 -42399.5980
66 -46853.9476 -524936.9858
67 26048.6923 -46853.9476
68 -19831.6724 26048.6923
69 -42390.0984 -19831.6724
70 -197150.8329 -42390.0984
71 -51774.5780 -197150.8329
72 -35135.1327 -51774.5780
73 32526.8323 -35135.1327
74 38018.9463 32526.8323
75 -86784.1679 38018.9463
76 39064.2674 -86784.1679
77 -64620.0625 39064.2674
78 -52160.0646 -64620.0625
79 39018.8929 -52160.0646
80 10336.1251 39018.8929
81 -30491.9024 10336.1251
82 21645.3152 -30491.9024
83 -4747.8623 21645.3152
84 77437.5254 -4747.8623
85 -124202.2780 77437.5254
86 -19693.9857 -124202.2780
87 56162.0558 -19693.9857
88 -64753.8158 56162.0558
89 86266.4056 -64753.8158
90 33704.0154 86266.4056
91 -6989.6710 33704.0154
92 46043.9995 -6989.6710
93 -84723.7008 46043.9995
94 -7575.3138 -84723.7008
95 1729.0119 -7575.3138
96 142058.1967 1729.0119
97 1785.5551 142058.1967
98 53560.0816 1785.5551
99 6579.9276 53560.0816
100 -29502.7850 6579.9276
101 -80171.6983 -29502.7850
102 -114712.4397 -80171.6983
103 21671.5858 -114712.4397
104 25745.8450 21671.5858
105 -53544.7222 25745.8450
106 -275735.4076 -53544.7222
107 176036.5982 -275735.4076
108 -132744.2026 176036.5982
109 -87592.5503 -132744.2026
110 90557.5062 -87592.5503
111 -228631.1611 90557.5062
112 -146148.2854 -228631.1611
113 -49099.8440 -146148.2854
114 70476.2658 -49099.8440
115 -4009.2257 70476.2658
116 44140.8107 -4009.2257
117 102558.5655 44140.8107
118 -118258.5694 102558.5655
119 816.2494 -118258.5694
120 -243987.8612 816.2494
121 30189.5715 -243987.8612
122 -92740.5446 30189.5715
123 142401.7259 -92740.5446
124 3826.0383 142401.7259
125 -153906.6564 3826.0383
126 45876.2528 -153906.6564
127 30075.6448 45876.2528
128 9665.3884 30075.6448
129 45876.2528 9665.3884
130 35076.0683 45876.2528
131 33673.0231 35076.0683
132 45876.2528 33673.0231
133 16601.0575 45876.2528
134 23327.4580 16601.0575
135 33445.0498 23327.4580
136 33987.1176 33445.0498
137 43983.9643 33987.1176
138 20663.1323 43983.9643
139 45876.2528 20663.1323
140 37187.6608 45876.2528
141 45876.2528 37187.6608
142 -7005.3106 45876.2528
143 40383.7450 -7005.3106
144 30884.5952 40383.7450
145 25016.3156 30884.5952
146 45876.2528 25016.3156
147 36255.9339 45876.2528
148 47247.8386 36255.9339
149 46067.5680 47247.8386
150 1761.4402 46067.5680
151 -61450.0008 1761.4402
152 39425.6762 -61450.0008
153 9230.1870 39425.6762
154 42724.3296 9230.1870
155 45876.2528 42724.3296
156 45876.2528 45876.2528
157 41012.2288 45876.2528
158 51627.6528 41012.2288
159 45876.2528 51627.6528
160 45876.2528 45876.2528
161 45876.2528 45876.2528
162 45876.2528 45876.2528
163 17187.3214 45876.2528
164 -8192.4349 17187.3214
165 43061.9100 -8192.4349
166 45876.2528 43061.9100
167 73457.3654 45876.2528
168 -100530.2721 73457.3654
169 45876.2528 -100530.2721
170 45876.2528 45876.2528
171 38251.8603 45876.2528
172 45876.2528 38251.8603
173 24048.7859 45876.2528
174 45876.2528 24048.7859
175 32971.1674 45876.2528
176 -110753.6285 32971.1674
177 45121.9943 -110753.6285
178 45876.2528 45121.9943
179 45876.2528 45876.2528
180 44095.8056 45876.2528
181 34778.7980 44095.8056
182 -41202.6439 34778.7980
183 -986.0808 -41202.6439
184 38349.4798 -986.0808
185 35584.4931 38349.4798
186 2711.9615 35584.4931
187 32834.2016 2711.9615
188 10898.0508 32834.2016
189 45876.2528 10898.0508
190 34205.9235 45876.2528
191 45876.2528 34205.9235
192 33752.7633 45876.2528
193 45876.2528 33752.7633
194 41702.5373 45876.2528
195 45876.2528 41702.5373
196 -83418.4463 45876.2528
197 45876.2528 -83418.4463
198 45876.2528 45876.2528
199 -34734.8816 45876.2528
200 45876.2528 -34734.8816
201 45876.2528 45876.2528
202 45876.2528 45876.2528
203 11544.4399 45876.2528
204 29663.8946 11544.4399
205 39982.3299 29663.8946
206 45114.4407 39982.3299
207 45876.2528 45114.4407
208 15762.1198 45876.2528
209 27787.7310 15762.1198
210 45876.2528 27787.7310
211 45876.2528 45876.2528
212 40534.1392 45876.2528
213 -3769.5130 40534.1392
214 16870.9861 -3769.5130
215 43296.0221 16870.9861
216 38116.0415 43296.0221
217 30533.0992 38116.0415
218 43296.0221 30533.0992
219 -24039.3876 43296.0221
220 16993.5346 -24039.3876
221 45876.2528 16993.5346
222 44494.6336 45876.2528
223 6636.3416 44494.6336
224 44345.4449 6636.3416
225 10965.6898 44345.4449
226 32104.3027 10965.6898
227 45876.2528 32104.3027
228 36062.7740 45876.2528
229 45876.2528 36062.7740
230 13243.8692 45876.2528
231 24147.2612 13243.8692
232 45876.2528 24147.2612
233 8919.6114 45876.2528
234 45876.2528 8919.6114
235 45939.5835 45876.2528
236 45876.2528 45939.5835
237 37391.3257 45876.2528
238 45876.2528 37391.3257
239 45876.2528 45876.2528
240 45876.2528 45876.2528
241 37460.4520 45876.2528
242 45876.2528 37460.4520
243 36578.4145 45876.2528
244 1382.6751 36578.4145
245 45876.2528 1382.6751
246 45876.2528 45876.2528
247 17273.1870 45876.2528
248 45876.2528 17273.1870
249 45876.2528 45876.2528
250 11089.2820 45876.2528
251 45876.2528 11089.2820
252 32881.7550 45876.2528
253 -27733.2665 32881.7550
254 -140283.3289 -27733.2665
255 45876.2528 -140283.3289
256 23055.8072 45876.2528
257 -12528.5493 23055.8072
258 45876.2528 -12528.5493
259 45876.2528 45876.2528
260 45876.2528 45876.2528
261 45876.2528 45876.2528
262 52719.1282 45876.2528
263 42846.8695 52719.1282
264 45876.2528 42846.8695
265 45876.2528 45876.2528
266 45876.2528 45876.2528
267 95960.3053 45876.2528
268 45876.2528 95960.3053
269 25360.0761 45876.2528
270 45876.2528 25360.0761
271 45876.2528 45876.2528
272 45876.2528 45876.2528
273 13455.7335 45876.2528
274 45876.2528 13455.7335
275 28873.5045 45876.2528
276 13907.0008 28873.5045
277 45876.2528 13907.0008
278 -43985.1358 45876.2528
279 27002.4361 -43985.1358
280 45876.2528 27002.4361
281 45876.2528 45876.2528
282 5007.4327 45876.2528
283 45876.2528 5007.4327
284 42937.2846 45876.2528
285 45876.2528 42937.2846
286 30133.0662 45876.2528
287 44068.8680 30133.0662
288 45876.2528 44068.8680
289 1041.9318 45876.2528
290 45876.2528 1041.9318
291 36588.1557 45876.2528
292 62051.9763 36588.1557
293 45876.2528 62051.9763
294 18846.8644 45876.2528
295 51512.3879 18846.8644
296 45876.2528 51512.3879
297 45876.2528 45876.2528
298 45876.2528 45876.2528
299 34265.2149 45876.2528
300 37844.7527 34265.2149
301 24038.0960 37844.7527
302 44586.1374 24038.0960
303 45876.2528 44586.1374
304 48868.1653 45876.2528
305 30365.1997 48868.1653
306 45876.2528 30365.1997
307 18067.6250 45876.2528
308 10015.7961 18067.6250
309 45876.2528 10015.7961
310 45876.2528 45876.2528
311 27853.9242 45876.2528
312 45876.2528 27853.9242
313 33395.7416 45876.2528
314 35765.8263 33395.7416
315 61876.7863 35765.8263
316 44761.5447 61876.7863
317 -6200.1958 44761.5447
318 39395.8415 -6200.1958
319 45876.2528 39395.8415
320 45876.2528 45876.2528
321 45876.2528 45876.2528
322 45876.2528 45876.2528
323 45876.2528 45876.2528
324 45876.2528 45876.2528
325 41227.7757 45876.2528
326 45876.2528 41227.7757
327 45876.2528 45876.2528
328 45876.2528 45876.2528
329 45876.2528 45876.2528
330 45876.2528 45876.2528
331 45876.2528 45876.2528
332 45876.2528 45876.2528
333 36968.6674 45876.2528
334 44572.3028 36968.6674
335 33850.7190 44572.3028
336 -17595.4776 33850.7190
337 45876.2528 -17595.4776
338 42518.4372 45876.2528
339 32563.7714 42518.4372
340 3957.1427 32563.7714
341 45876.2528 3957.1427
342 45876.2528 45876.2528
343 45876.2528 45876.2528
344 8421.4511 45876.2528
345 45876.2528 8421.4511
346 27704.5864 45876.2528
347 45876.2528 27704.5864
348 45876.2528 45876.2528
349 2137.8717 45876.2528
350 38503.5036 2137.8717
351 34204.3680 38503.5036
352 30477.1743 34204.3680
353 40468.4630 30477.1743
354 45876.2528 40468.4630
355 45876.2528 45876.2528
356 22299.5254 45876.2528
357 45876.2528 22299.5254
358 45876.2528 45876.2528
359 45876.2528 45876.2528
360 40711.9569 45876.2528
361 45876.2528 40711.9569
362 45876.2528 45876.2528
363 -35785.4212 45876.2528
364 31340.9823 -35785.4212
365 45876.2528 31340.9823
366 45876.2528 45876.2528
367 45876.2528 45876.2528
368 45876.2528 45876.2528
369 32258.4971 45876.2528
370 45876.2528 32258.4971
371 35044.5333 45876.2528
372 53147.9855 35044.5333
373 30250.6819 53147.9855
374 45876.2528 30250.6819
375 45876.2528 45876.2528
376 9032.9946 45876.2528
377 30551.8278 9032.9946
378 -10839.0923 30551.8278
379 45876.2528 -10839.0923
380 45876.2528 45876.2528
381 12992.8254 45876.2528
382 -19634.9329 12992.8254
383 45876.2528 -19634.9329
384 30641.7681 45876.2528
385 -123138.4747 30641.7681
386 50836.5709 -123138.4747
387 11640.2014 50836.5709
388 39263.8138 11640.2014
389 15615.9244 39263.8138
390 31253.7532 15615.9244
391 5065.9793 31253.7532
392 -30129.8069 5065.9793
393 75257.3307 -30129.8069
394 -71988.3779 75257.3307
395 -8508.7294 -71988.3779
396 -166757.5998 -8508.7294
397 -170833.8386 -166757.5998
398 6175.9286 -170833.8386
399 -138197.8894 6175.9286
400 7711.3690 -138197.8894
401 36736.6607 7711.3690
402 6333.6046 36736.6607
403 -50925.8076 6333.6046
404 -14562.6528 -50925.8076
405 17503.4568 -14562.6528
406 19791.6817 17503.4568
407 -2024.5758 19791.6817
408 7417.6489 -2024.5758
409 -57365.0314 7417.6489
410 -62068.8376 -57365.0314
411 -171624.2565 -62068.8376
412 -44517.9813 -171624.2565
413 65182.9517 -44517.9813
414 -160213.4970 65182.9517
415 -88829.3916 -160213.4970
416 36889.0594 -88829.3916
417 -29612.3537 36889.0594
418 -52885.4894 -29612.3537
419 -38685.5501 -52885.4894
420 -39444.3383 -38685.5501
421 -186172.9561 -39444.3383
422 -692206.1533 -186172.9561
423 -479941.8205 -692206.1533
424 -72026.6461 -479941.8205
425 -271484.3057 -72026.6461
426 -341158.8368 -271484.3057
427 -1055853.5887 -341158.8368
428 -259264.3938 -1055853.5887
429 -340495.0209 -259264.3938
430 -760308.8934 -340495.0209
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7l64k1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8l64k1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/9eg3n1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/10eg3n1292875566.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11hykb1292875566.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/122yih1292875566.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13mtnw1292875566.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/14k9fe1292875566.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/15nadk1292875566.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/16j1ba1292875566.tab")
+ }
>
> try(system("convert tmp/17eoc1292875566.ps tmp/17eoc1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/20o6x1292875566.ps tmp/20o6x1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/30o6x1292875566.ps tmp/30o6x1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/40o6x1292875566.ps tmp/40o6x1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/5axnh1292875566.ps tmp/5axnh1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/6axnh1292875566.ps tmp/6axnh1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/7l64k1292875566.ps tmp/7l64k1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/8l64k1292875566.ps tmp/8l64k1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/9eg3n1292875566.ps tmp/9eg3n1292875566.png",intern=TRUE))
character(0)
> try(system("convert tmp/10eg3n1292875566.ps tmp/10eg3n1292875566.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.253 2.312 23.541