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(0
+ ,807
+ ,6282154
+ ,0
+ ,444
+ ,4321023
+ ,0
+ ,412
+ ,4111912
+ ,1
+ ,428
+ ,223193
+ ,0
+ ,315
+ ,1491348
+ ,0
+ ,168
+ ,1629616
+ ,1
+ ,263
+ ,1398893
+ ,0
+ ,267
+ ,1926517
+ ,0
+ ,228
+ ,983660
+ ,0
+ ,129
+ ,1443586
+ ,1
+ ,104
+ ,1073089
+ ,1
+ ,122
+ ,984885
+ ,0
+ ,393
+ ,1405225
+ ,1
+ ,190
+ ,227132
+ ,0
+ ,280
+ ,929118
+ ,1
+ ,63
+ ,1071292
+ ,1
+ ,102
+ ,638830
+ ,0
+ ,265
+ ,856956
+ ,0
+ ,234
+ ,992426
+ ,1
+ ,277
+ ,444477
+ ,0
+ ,73
+ ,857217
+ ,0
+ ,67
+ ,711969
+ ,1
+ ,103
+ ,702380
+ ,1
+ ,290
+ ,358589
+ ,1
+ ,83
+ ,297978
+ ,1
+ ,56
+ ,585715
+ ,0
+ ,236
+ ,657954
+ ,1
+ ,73
+ ,209458
+ ,1
+ ,34
+ ,786690
+ ,1
+ ,139
+ ,439798
+ ,0
+ ,26
+ ,688779
+ ,0
+ ,70
+ ,574339
+ ,0
+ ,40
+ ,741409
+ ,0
+ ,42
+ ,597793
+ ,1
+ ,12
+ ,644190
+ ,1
+ ,211
+ ,377934
+ ,1
+ ,74
+ ,640273
+ ,0
+ ,80
+ ,697458
+ ,0
+ ,83
+ ,550608
+ ,1
+ ,131
+ ,207393
+ ,1
+ ,203
+ ,301607
+ ,1
+ ,56
+ ,345783
+ ,1
+ ,89
+ ,501749
+ ,1
+ ,88
+ ,379983
+ ,1
+ ,39
+ ,387475
+ ,0
+ ,25
+ ,377305
+ ,0
+ ,49
+ ,370837
+ ,0
+ ,149
+ ,430866
+ ,0
+ ,58
+ ,469107
+ ,1
+ ,41
+ ,194493
+ ,0
+ ,90
+ ,530670
+ ,0
+ ,136
+ ,518365
+ ,0
+ ,97
+ ,491303
+ ,0
+ ,63
+ ,527021
+ ,0
+ ,114
+ ,233773
+ ,1
+ ,77
+ ,405972
+ ,1
+ ,6
+ ,652925
+ ,1
+ ,47
+ ,446211
+ ,1
+ ,51
+ ,341340
+ ,0
+ ,85
+ ,387699
+ ,0
+ ,43
+ ,493408
+ ,1
+ ,32
+ ,146494
+ ,0
+ ,25
+ ,414462
+ ,0
+ ,77
+ ,364304
+ ,1
+ ,54
+ ,355178
+ ,1
+ ,251
+ ,357760
+ ,1
+ ,15
+ ,261216
+ ,0
+ ,44
+ ,397144
+ ,1
+ ,73
+ ,374943
+ ,0
+ ,85
+ ,424898
+ ,0
+ ,49
+ ,202055
+ ,1
+ ,38
+ ,378525
+ ,1
+ ,35
+ ,310768
+ ,1
+ ,9
+ ,325738
+ ,1
+ ,34
+ ,394510
+ ,0
+ ,20
+ ,247060
+ ,1
+ ,29
+ ,368078
+ ,1
+ ,11
+ ,236761
+ ,1
+ ,52
+ ,312378
+ ,0
+ ,13
+ ,339836
+ ,1
+ ,29
+ ,347385
+ ,0
+ ,66
+ ,426280
+ ,1
+ ,33
+ ,352850
+ ,1
+ ,15
+ ,301881
+ ,1
+ ,15
+ ,377516
+ ,0
+ ,68
+ ,357312
+ ,1
+ ,100
+ ,458343
+ ,1
+ ,13
+ ,354228
+ ,1
+ ,45
+ ,308636
+ ,1
+ ,14
+ ,386212
+ ,1
+ ,36
+ ,393343
+ ,0
+ ,40
+ ,378509
+ ,1
+ ,68
+ ,452469
+ ,0
+ ,29
+ ,364839
+ ,1
+ ,43
+ ,358649
+ ,0
+ ,30
+ ,376641
+ ,1
+ ,9
+ ,429112
+ ,0
+ ,22
+ ,330546
+ ,1
+ ,19
+ ,403560
+ ,0
+ ,9
+ ,317892
+ ,0
+ ,31
+ ,307528
+ ,1
+ ,19
+ ,235133
+ ,1
+ ,55
+ ,299243
+ ,1
+ ,8
+ ,314073
+ ,1
+ ,28
+ ,368186
+ ,1
+ ,29
+ ,269661
+ ,0
+ ,48
+ ,125390
+ ,0
+ ,16
+ ,510834
+ ,1
+ ,47
+ ,321896
+ ,0
+ ,20
+ ,249898
+ ,1
+ ,22
+ ,408881
+ ,0
+ ,33
+ ,158492
+ ,1
+ ,44
+ ,292154
+ ,0
+ ,13
+ ,289513
+ ,0
+ ,6
+ ,378049
+ ,1
+ ,35
+ ,343466
+ ,1
+ ,8
+ ,332743
+ ,1
+ ,17
+ ,442882
+ ,0
+ ,11
+ ,214215
+ ,1
+ ,21
+ ,315688
+ ,1
+ ,92
+ ,375195
+ ,1
+ ,12
+ ,334280
+ ,1
+ ,112
+ ,355864
+ ,0
+ ,25
+ ,480382
+ ,0
+ ,17
+ ,353058
+ ,0
+ ,23
+ ,217193
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,10
+ ,314533
+ ,1
+ ,23
+ ,318056
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,314353
+ ,1
+ ,25
+ ,369448
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,20
+ ,312846
+ ,1
+ ,4
+ ,312075
+ ,1
+ ,4
+ ,315009
+ ,1
+ ,10
+ ,318903
+ ,1
+ ,1
+ ,314887
+ ,1
+ ,4
+ ,314913
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,8
+ ,325506
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,11
+ ,298568
+ ,1
+ ,4
+ ,315834
+ ,1
+ ,15
+ ,329784
+ ,1
+ ,9
+ ,312878
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,314987
+ ,1
+ ,2
+ ,325249
+ ,1
+ ,0
+ ,315877
+ ,1
+ ,7
+ ,291650
+ ,1
+ ,46
+ ,305959
+ ,1
+ ,5
+ ,315380
+ ,1
+ ,7
+ ,297765
+ ,1
+ ,2
+ ,315245
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,2
+ ,315236
+ ,1
+ ,5
+ ,336425
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,306268
+ ,1
+ ,24
+ ,302187
+ ,1
+ ,1
+ ,314882
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,18
+ ,382712
+ ,1
+ ,55
+ ,341570
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,3
+ ,312412
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,9
+ ,309596
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,8
+ ,315547
+ ,1
+ ,113
+ ,313267
+ ,0
+ ,0
+ ,316176
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,19
+ ,359335
+ ,0
+ ,11
+ ,330068
+ ,1
+ ,25
+ ,314289
+ ,0
+ ,16
+ ,297413
+ ,0
+ ,5
+ ,314806
+ ,0
+ ,11
+ ,333210
+ ,0
+ ,23
+ ,352108
+ ,1
+ ,6
+ ,313332
+ ,1
+ ,5
+ ,291787
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,318745
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,315366
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,3
+ ,315688
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,89
+ ,409642
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,19
+ ,269587
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,12
+ ,300962
+ ,0
+ ,12
+ ,325479
+ ,0
+ ,5
+ ,316155
+ ,0
+ ,2
+ ,318574
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,26
+ ,343613
+ ,0
+ ,3
+ ,306948
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,11
+ ,330059
+ ,1
+ ,10
+ ,288985
+ ,1
+ ,5
+ ,304485
+ ,1
+ ,2
+ ,315380
+ ,1
+ ,6
+ ,315688
+ ,1
+ ,7
+ ,317736
+ ,1
+ ,2
+ ,315380
+ ,1
+ ,28
+ ,322331
+ ,1
+ ,3
+ ,296656
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,1
+ ,315354
+ ,1
+ ,20
+ ,312161
+ ,1
+ ,1
+ ,315576
+ ,1
+ ,22
+ ,314922
+ ,1
+ ,9
+ ,314551
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,2
+ ,312339
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,298700
+ ,1
+ ,9
+ ,321376
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,13
+ ,303230
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315487
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,6
+ ,315793
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,3
+ ,312887
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,315637
+ ,1
+ ,2
+ ,324385
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,15
+ ,308989
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,9
+ ,296702
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,1
+ ,307322
+ ,1
+ ,38
+ ,304376
+ ,1
+ ,57
+ ,253588
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,309560
+ ,1
+ ,26
+ ,298466
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,13
+ ,343929
+ ,1
+ ,10
+ ,331955
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,9
+ ,381180
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,26
+ ,331420
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,19
+ ,310201
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,12
+ ,320016
+ ,1
+ ,23
+ ,320398
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,29
+ ,291841
+ ,1
+ ,8
+ ,310670
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,26
+ ,313491
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,9
+ ,331323
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,5
+ ,319210
+ ,1
+ ,3
+ ,318098
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,13
+ ,292754
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,12
+ ,325176
+ ,1
+ ,19
+ ,365959
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,10
+ ,302409
+ ,0
+ ,9
+ ,340968
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,9
+ ,315380
+ ,0
+ ,4
+ ,313164
+ ,1
+ ,1
+ ,301164
+ ,1
+ ,1
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,14
+ ,344425
+ ,1
+ ,12
+ ,315394
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,19
+ ,316647
+ ,1
+ ,17
+ ,309836
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,32
+ ,346611
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,14
+ ,322031
+ ,1
+ ,8
+ ,315656
+ ,1
+ ,4
+ ,339445
+ ,1
+ ,0
+ ,314964
+ ,1
+ ,20
+ ,297141
+ ,1
+ ,5
+ ,315372
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,1
+ ,312502
+ ,0
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,4
+ ,313729
+ ,1
+ ,1
+ ,315388
+ ,1
+ ,4
+ ,315371
+ ,1
+ ,20
+ ,296139
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,1
+ ,313880
+ ,1
+ ,10
+ ,317698
+ ,1
+ ,12
+ ,295580
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,13
+ ,308256
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,3
+ ,303677
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,10
+ ,319369
+ ,1
+ ,3
+ ,318690
+ ,1
+ ,7
+ ,314049
+ ,1
+ ,10
+ ,325699
+ ,1
+ ,1
+ ,314210
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,15
+ ,322378
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,4
+ ,315398
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,0
+ ,28
+ ,308336
+ ,1
+ ,9
+ ,316386
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,315553
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,7
+ ,323361
+ ,1
+ ,7
+ ,336639
+ ,1
+ ,3
+ ,307424
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,11
+ ,295370
+ ,1
+ ,7
+ ,322340
+ ,1
+ ,10
+ ,319864
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,18
+ ,317291
+ ,1
+ ,14
+ ,280398
+ ,1
+ ,0
+ ,315380
+ ,1
+ ,12
+ ,317330
+ ,1
+ ,29
+ ,238125
+ ,1
+ ,3
+ ,327071
+ ,1
+ ,6
+ ,309038
+ ,1
+ ,3
+ ,314210
+ ,1
+ ,8
+ ,307930
+ ,1
+ ,10
+ ,322327
+ ,1
+ ,6
+ ,292136
+ ,1
+ ,8
+ ,263276
+ ,1
+ ,6
+ ,367655
+ ,0
+ ,9
+ ,283910
+ ,1
+ ,8
+ ,283587
+ ,0
+ ,26
+ ,243650
+ ,0
+ ,239
+ ,438493
+ ,1
+ ,7
+ ,296261
+ ,0
+ ,41
+ ,230621
+ ,1
+ ,3
+ ,304252
+ ,1
+ ,8
+ ,333505
+ ,1
+ ,6
+ ,296919
+ ,0
+ ,21
+ ,278990
+ ,1
+ ,7
+ ,276898
+ ,1
+ ,11
+ ,327007
+ ,1
+ ,11
+ ,317046
+ ,1
+ ,12
+ ,304555
+ ,1
+ ,9
+ ,298096
+ ,1
+ ,3
+ ,231861
+ ,1
+ ,57
+ ,309422
+ ,0
+ ,21
+ ,286963
+ ,0
+ ,15
+ ,269753
+ ,0
+ ,32
+ ,448243
+ ,1
+ ,11
+ ,165404
+ ,1
+ ,2
+ ,204325
+ ,1
+ ,23
+ ,407159
+ ,0
+ ,20
+ ,290476
+ ,1
+ ,24
+ ,275311
+ ,1
+ ,1
+ ,246541
+ ,1
+ ,1
+ ,253468
+ ,1
+ ,74
+ ,240897
+ ,0
+ ,68
+ ,-83265
+ ,1
+ ,20
+ ,-42143
+ ,1
+ ,20
+ ,272713
+ ,0
+ ,82
+ ,215362
+ ,1
+ ,21
+ ,42754
+ ,1
+ ,244
+ ,306275
+ ,1
+ ,32
+ ,253537
+ ,1
+ ,86
+ ,372631
+ ,1
+ ,69
+ ,-7170)
+ ,dim=c(3
+ ,431)
+ ,dimnames=list(c('group'
+ ,'orders'
+ ,'wealth
')
+ ,1:431))
> y <- array(NA,dim=c(3,431),dimnames=list(c('group','orders','wealth
'),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
wealth\r group orders
1 6282154 0 807
2 4321023 0 444
3 4111912 0 412
4 223193 1 428
5 1491348 0 315
6 1629616 0 168
7 1398893 1 263
8 1926517 0 267
9 983660 0 228
10 1443586 0 129
11 1073089 1 104
12 984885 1 122
13 1405225 0 393
14 227132 1 190
15 929118 0 280
16 1071292 1 63
17 638830 1 102
18 856956 0 265
19 992426 0 234
20 444477 1 277
21 857217 0 73
22 711969 0 67
23 702380 1 103
24 358589 1 290
25 297978 1 83
26 585715 1 56
27 657954 0 236
28 209458 1 73
29 786690 1 34
30 439798 1 139
31 688779 0 26
32 574339 0 70
33 741409 0 40
34 597793 0 42
35 644190 1 12
36 377934 1 211
37 640273 1 74
38 697458 0 80
39 550608 0 83
40 207393 1 131
41 301607 1 203
42 345783 1 56
43 501749 1 89
44 379983 1 88
45 387475 1 39
46 377305 0 25
47 370837 0 49
48 430866 0 149
49 469107 0 58
50 194493 1 41
51 530670 0 90
52 518365 0 136
53 491303 0 97
54 527021 0 63
55 233773 0 114
56 405972 1 77
57 652925 1 6
58 446211 1 47
59 341340 1 51
60 387699 0 85
61 493408 0 43
62 146494 1 32
63 414462 0 25
64 364304 0 77
65 355178 1 54
66 357760 1 251
67 261216 1 15
68 397144 0 44
69 374943 1 73
70 424898 0 85
71 202055 0 49
72 378525 1 38
73 310768 1 35
74 325738 1 9
75 394510 1 34
76 247060 0 20
77 368078 1 29
78 236761 1 11
79 312378 1 52
80 339836 0 13
81 347385 1 29
82 426280 0 66
83 352850 1 33
84 301881 1 15
85 377516 1 15
86 357312 0 68
87 458343 1 100
88 354228 1 13
89 308636 1 45
90 386212 1 14
91 393343 1 36
92 378509 0 40
93 452469 1 68
94 364839 0 29
95 358649 1 43
96 376641 0 30
97 429112 1 9
98 330546 0 22
99 403560 1 19
100 317892 0 9
101 307528 0 31
102 235133 1 19
103 299243 1 55
104 314073 1 8
105 368186 1 28
106 269661 1 29
107 125390 0 48
108 510834 0 16
109 321896 1 47
110 249898 0 20
111 408881 1 22
112 158492 0 33
113 292154 1 44
114 289513 0 13
115 378049 0 6
116 343466 1 35
117 332743 1 8
118 442882 1 17
119 214215 0 11
120 315688 1 21
121 375195 1 92
122 334280 1 12
123 355864 1 112
124 480382 0 25
125 353058 0 17
126 217193 0 23
127 315380 1 0
128 314533 1 10
129 318056 1 23
130 315380 1 0
131 314353 1 7
132 369448 1 25
133 315380 1 0
134 312846 1 20
135 312075 1 4
136 315009 1 4
137 318903 1 10
138 314887 1 1
139 314913 1 4
140 315380 1 0
141 325506 1 8
142 315380 1 0
143 298568 1 11
144 315834 1 4
145 329784 1 15
146 312878 1 9
147 315380 1 0
148 314987 1 7
149 325249 1 2
150 315877 1 0
151 291650 1 7
152 305959 1 46
153 315380 1 5
154 297765 1 7
155 315245 1 2
156 315380 1 0
157 315380 1 0
158 315236 1 2
159 336425 1 5
160 315380 1 0
161 315380 1 0
162 315380 1 0
163 315380 1 0
164 306268 1 7
165 302187 1 24
166 314882 1 1
167 315380 1 0
168 382712 1 18
169 341570 1 55
170 315380 1 0
171 315380 1 0
172 312412 1 3
173 315380 1 0
174 309596 1 9
175 315380 1 0
176 315547 1 8
177 313267 1 113
178 316176 0 0
179 315380 1 0
180 315380 1 0
181 359335 1 19
182 330068 0 11
183 314289 1 25
184 297413 0 16
185 314806 0 5
186 333210 0 11
187 352108 0 23
188 313332 1 6
189 291787 1 5
190 315380 1 0
191 318745 1 7
192 315380 1 0
193 315366 1 7
194 315380 1 0
195 315688 1 3
196 315380 0 0
197 409642 0 89
198 315380 0 0
199 315380 0 0
200 269587 0 19
201 315380 0 0
202 315380 0 0
203 315380 0 0
204 300962 0 12
205 325479 0 12
206 316155 0 5
207 318574 0 2
208 315380 0 0
209 343613 0 26
210 306948 0 3
211 315380 1 0
212 315380 1 0
213 330059 1 11
214 288985 1 10
215 304485 1 5
216 315380 1 2
217 315688 1 6
218 317736 1 7
219 315380 1 2
220 322331 1 28
221 296656 1 3
222 315380 1 0
223 315354 1 1
224 312161 1 20
225 315576 1 1
226 314922 1 22
227 314551 1 9
228 315380 1 0
229 312339 1 2
230 315380 1 0
231 298700 1 7
232 321376 1 9
233 315380 1 0
234 303230 1 13
235 315380 1 0
236 315487 1 0
237 315380 1 0
238 315793 1 6
239 315380 1 0
240 315380 1 0
241 315380 1 0
242 312887 1 3
243 315380 1 0
244 315637 1 7
245 324385 1 2
246 315380 0 0
247 315380 0 0
248 308989 1 15
249 315380 0 0
250 315380 0 0
251 296702 1 9
252 315380 0 0
253 307322 1 1
254 304376 1 38
255 253588 1 57
256 315380 1 0
257 309560 1 7
258 298466 1 26
259 315380 0 0
260 315380 1 0
261 315380 0 0
262 315380 0 0
263 343929 1 13
264 331955 1 10
265 315380 0 0
266 315380 1 0
267 315380 0 0
268 381180 1 9
269 315380 0 0
270 331420 1 26
271 315380 0 0
272 315380 0 0
273 315380 0 0
274 310201 1 19
275 315380 1 0
276 320016 1 12
277 320398 1 23
278 315380 1 0
279 291841 0 29
280 310670 1 8
281 315380 0 0
282 315380 1 0
283 313491 1 26
284 315380 1 0
285 331323 1 9
286 315380 1 0
287 319210 0 5
288 318098 1 3
289 315380 0 0
290 292754 1 13
291 315380 0 0
292 325176 1 12
293 365959 1 19
294 315380 1 0
295 302409 1 10
296 340968 0 9
297 315380 1 0
298 315380 0 0
299 315380 1 0
300 315380 1 9
301 313164 0 4
302 301164 1 1
303 315380 1 1
304 315380 0 0
305 344425 1 14
306 315394 1 12
307 315380 0 0
308 316647 0 19
309 309836 1 17
310 315380 0 0
311 315380 0 0
312 346611 1 32
313 315380 0 0
314 322031 0 14
315 315656 1 8
316 339445 1 4
317 314964 1 0
318 297141 1 20
319 315372 1 5
320 315380 1 0
321 315380 1 0
322 315380 1 0
323 315380 1 0
324 315380 0 0
325 315380 1 0
326 312502 1 1
327 315380 0 0
328 315380 0 0
329 315380 1 0
330 315380 1 0
331 315380 1 0
332 315380 1 0
333 315380 1 0
334 313729 1 4
335 315388 1 1
336 315371 1 4
337 296139 1 20
338 315380 1 0
339 313880 1 1
340 317698 1 10
341 295580 1 12
342 315380 1 0
343 315380 1 0
344 315380 1 0
345 308256 1 13
346 315380 1 0
347 303677 1 3
348 315380 1 0
349 315380 1 0
350 319369 1 10
351 318690 1 3
352 314049 1 7
353 325699 1 10
354 314210 1 1
355 315380 1 0
356 315380 1 0
357 322378 1 15
358 315380 1 0
359 315380 1 0
360 315380 1 0
361 315398 1 4
362 315380 1 0
363 315380 1 0
364 308336 0 28
365 316386 1 9
366 315380 1 0
367 315380 1 0
368 315380 1 0
369 315380 1 0
370 315553 1 7
371 315380 1 0
372 323361 1 7
373 336639 1 7
374 307424 1 3
375 315380 1 0
376 315380 1 0
377 295370 1 11
378 322340 1 7
379 319864 1 10
380 315380 1 0
381 315380 1 0
382 317291 1 18
383 280398 1 14
384 315380 1 0
385 317330 1 12
386 238125 1 29
387 327071 1 3
388 309038 1 6
389 314210 1 3
390 307930 1 8
391 322327 1 10
392 292136 1 6
393 263276 1 8
394 367655 1 6
395 283910 0 9
396 283587 1 8
397 243650 0 26
398 438493 0 239
399 296261 1 7
400 230621 0 41
401 304252 1 3
402 333505 1 8
403 296919 1 6
404 278990 0 21
405 276898 1 7
406 327007 1 11
407 317046 1 11
408 304555 1 12
409 298096 1 9
410 231861 1 3
411 309422 1 57
412 286963 0 21
413 269753 0 15
414 448243 0 32
415 165404 1 11
416 204325 1 2
417 407159 1 23
418 290476 0 20
419 275311 1 24
420 246541 1 1
421 253468 1 1
422 240897 1 74
423 -83265 0 68
424 -42143 1 20
425 272713 1 20
426 215362 0 82
427 42754 1 21
428 306275 1 244
429 253537 1 32
430 372631 1 86
431 -7170 1 69
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) group orders
282194 -50104 4263
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-1833510 -30828 34083 83290 2559626
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 282193.7 28180.1 10.014 <2e-16 ***
group -50104.1 31230.8 -1.604 0.109
orders 4263.1 189.1 22.544 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 285000 on 428 degrees of freedom
Multiple R-squared: 0.5637, Adjusted R-squared: 0.5616
F-statistic: 276.5 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 6.726328e-23 3.363164e-23
[2,] 1.0000000 3.082179e-46 1.541090e-46
[3,] 1.0000000 1.377808e-54 6.889039e-55
[4,] 1.0000000 1.009301e-60 5.046506e-61
[5,] 1.0000000 3.857151e-68 1.928575e-68
[6,] 1.0000000 5.449678e-90 2.724839e-90
[7,] 1.0000000 1.344880e-97 6.724398e-98
[8,] 1.0000000 3.642576e-117 1.821288e-117
[9,] 1.0000000 4.090349e-118 2.045175e-118
[10,] 1.0000000 4.188647e-126 2.094324e-126
[11,] 1.0000000 2.376011e-141 1.188005e-141
[12,] 1.0000000 4.164841e-142 2.082421e-142
[13,] 1.0000000 2.348058e-148 1.174029e-148
[14,] 1.0000000 6.209040e-155 3.104520e-155
[15,] 1.0000000 2.594458e-157 1.297229e-157
[16,] 1.0000000 2.757674e-162 1.378837e-162
[17,] 1.0000000 2.027162e-164 1.013581e-164
[18,] 1.0000000 1.741313e-167 8.706567e-168
[19,] 1.0000000 1.289633e-170 6.448166e-171
[20,] 1.0000000 7.914985e-170 3.957492e-170
[21,] 1.0000000 4.878305e-172 2.439153e-172
[22,] 1.0000000 4.309772e-176 2.154886e-176
[23,] 1.0000000 1.187706e-175 5.938531e-176
[24,] 1.0000000 2.671289e-184 1.335644e-184
[25,] 1.0000000 2.097498e-183 1.048749e-183
[26,] 1.0000000 9.444548e-187 4.722274e-187
[27,] 1.0000000 1.803520e-187 9.017600e-188
[28,] 1.0000000 1.049531e-192 5.247657e-193
[29,] 1.0000000 2.317962e-194 1.158981e-194
[30,] 1.0000000 4.402765e-200 2.201383e-200
[31,] 1.0000000 3.347369e-200 1.673685e-200
[32,] 1.0000000 6.966684e-204 3.483342e-204
[33,] 1.0000000 8.854423e-209 4.427212e-209
[34,] 1.0000000 4.370387e-210 2.185194e-210
[35,] 1.0000000 1.656524e-210 8.282618e-211
[36,] 1.0000000 5.495587e-211 2.747794e-211
[37,] 1.0000000 4.317012e-210 2.158506e-210
[38,] 1.0000000 1.729027e-210 8.645136e-211
[39,] 1.0000000 1.693906e-209 8.469529e-210
[40,] 1.0000000 5.109832e-209 2.554916e-209
[41,] 1.0000000 3.853302e-208 1.926651e-208
[42,] 1.0000000 2.332572e-207 1.166286e-207
[43,] 1.0000000 4.025577e-208 2.012789e-208
[44,] 1.0000000 4.224760e-208 2.112380e-208
[45,] 1.0000000 5.132186e-208 2.566093e-208
[46,] 1.0000000 3.952471e-209 1.976236e-209
[47,] 1.0000000 1.776731e-210 8.883655e-211
[48,] 1.0000000 4.066384e-211 2.033192e-211
[49,] 1.0000000 1.630917e-212 8.154586e-213
[50,] 1.0000000 2.673706e-213 1.336853e-213
[51,] 1.0000000 1.115930e-212 5.579648e-213
[52,] 1.0000000 6.899460e-221 3.449730e-221
[53,] 1.0000000 3.472489e-221 1.736245e-221
[54,] 1.0000000 3.208432e-220 1.604216e-220
[55,] 1.0000000 1.140296e-219 5.701479e-220
[56,] 1.0000000 1.046703e-220 5.233516e-221
[57,] 1.0000000 1.032100e-221 5.160502e-222
[58,] 1.0000000 2.361317e-221 1.180658e-221
[59,] 1.0000000 1.229232e-220 6.146158e-221
[60,] 1.0000000 9.967937e-220 4.983968e-220
[61,] 1.0000000 8.192262e-221 4.096131e-221
[62,] 1.0000000 3.526720e-220 1.763360e-220
[63,] 1.0000000 1.321497e-219 6.607486e-220
[64,] 1.0000000 9.152588e-219 4.576294e-219
[65,] 1.0000000 9.860747e-219 4.930374e-219
[66,] 1.0000000 8.071744e-219 4.035872e-219
[67,] 1.0000000 2.956260e-218 1.478130e-218
[68,] 1.0000000 2.541529e-217 1.270765e-217
[69,] 1.0000000 1.262817e-216 6.314086e-217
[70,] 1.0000000 2.962382e-216 1.481191e-216
[71,] 1.0000000 1.562874e-215 7.814370e-216
[72,] 1.0000000 6.717016e-215 3.358508e-215
[73,] 1.0000000 2.002139e-214 1.001070e-214
[74,] 1.0000000 2.151002e-213 1.075501e-213
[75,] 1.0000000 2.108402e-212 1.054201e-212
[76,] 1.0000000 1.370947e-211 6.854733e-212
[77,] 1.0000000 1.543117e-211 7.715586e-212
[78,] 1.0000000 9.842808e-211 4.921404e-211
[79,] 1.0000000 7.137662e-210 3.568831e-210
[80,] 1.0000000 2.037030e-209 1.018515e-209
[81,] 1.0000000 1.142408e-208 5.712042e-209
[82,] 1.0000000 3.221133e-209 1.610567e-209
[83,] 1.0000000 1.486053e-208 7.430263e-209
[84,] 1.0000000 1.536360e-207 7.681802e-208
[85,] 1.0000000 3.434063e-207 1.717031e-207
[86,] 1.0000000 8.682390e-207 4.341195e-207
[87,] 1.0000000 3.838682e-206 1.919341e-206
[88,] 1.0000000 9.262027e-207 4.631014e-207
[89,] 1.0000000 5.494880e-206 2.747440e-206
[90,] 1.0000000 3.339792e-205 1.669896e-205
[91,] 1.0000000 1.459655e-204 7.298277e-205
[92,] 1.0000000 5.429811e-205 2.714906e-205
[93,] 1.0000000 5.179814e-204 2.589907e-204
[94,] 1.0000000 6.676329e-204 3.338164e-204
[95,] 1.0000000 6.885725e-203 3.442863e-203
[96,] 1.0000000 6.857358e-202 3.428679e-202
[97,] 1.0000000 2.398643e-201 1.199322e-201
[98,] 1.0000000 2.466999e-200 1.233499e-200
[99,] 1.0000000 1.808342e-199 9.041711e-200
[100,] 1.0000000 7.830099e-199 3.915050e-199
[101,] 1.0000000 5.810222e-198 2.905111e-198
[102,] 1.0000000 6.730912e-200 3.365456e-200
[103,] 1.0000000 1.887291e-202 9.436453e-203
[104,] 1.0000000 1.885190e-201 9.425951e-202
[105,] 1.0000000 9.856601e-201 4.928300e-201
[106,] 1.0000000 9.718754e-201 4.859377e-201
[107,] 1.0000000 7.827697e-202 3.913848e-202
[108,] 1.0000000 7.796084e-201 3.898042e-201
[109,] 1.0000000 7.491519e-200 3.745760e-200
[110,] 1.0000000 2.804295e-199 1.402147e-199
[111,] 1.0000000 2.081273e-198 1.040636e-198
[112,] 1.0000000 1.397426e-197 6.987130e-198
[113,] 1.0000000 2.093055e-198 1.046527e-198
[114,] 1.0000000 3.015376e-198 1.507688e-198
[115,] 1.0000000 2.725664e-197 1.362832e-197
[116,] 1.0000000 8.266059e-197 4.133030e-197
[117,] 1.0000000 5.780207e-196 2.890103e-196
[118,] 1.0000000 2.018862e-195 1.009431e-195
[119,] 1.0000000 2.959550e-197 1.479775e-197
[120,] 1.0000000 1.923554e-196 9.617768e-197
[121,] 1.0000000 3.018466e-196 1.509233e-196
[122,] 1.0000000 2.212589e-195 1.106294e-195
[123,] 1.0000000 1.835511e-194 9.177556e-195
[124,] 1.0000000 1.659612e-193 8.298060e-194
[125,] 1.0000000 1.235918e-192 6.179591e-193
[126,] 1.0000000 1.004367e-191 5.021837e-192
[127,] 1.0000000 3.623713e-191 1.811856e-191
[128,] 1.0000000 2.744110e-190 1.372055e-190
[129,] 1.0000000 2.504589e-189 1.252295e-189
[130,] 1.0000000 2.015884e-188 1.007942e-188
[131,] 1.0000000 1.607680e-187 8.038400e-188
[132,] 1.0000000 1.324437e-186 6.622186e-187
[133,] 1.0000000 1.030988e-185 5.154942e-186
[134,] 1.0000000 8.278354e-185 4.139177e-185
[135,] 1.0000000 6.402712e-184 3.201356e-184
[136,] 1.0000000 4.954088e-183 2.477044e-183
[137,] 1.0000000 3.848546e-182 1.924273e-182
[138,] 1.0000000 3.328470e-181 1.664235e-181
[139,] 1.0000000 2.674965e-180 1.337482e-180
[140,] 1.0000000 2.070347e-179 1.035173e-179
[141,] 1.0000000 1.745965e-178 8.729823e-179
[142,] 1.0000000 1.359607e-177 6.798034e-178
[143,] 1.0000000 1.119825e-176 5.599125e-177
[144,] 1.0000000 8.307593e-176 4.153796e-176
[145,] 1.0000000 6.464785e-175 3.232392e-175
[146,] 1.0000000 5.239891e-174 2.619945e-174
[147,] 1.0000000 4.797956e-173 2.398978e-173
[148,] 1.0000000 3.858196e-172 1.929098e-172
[149,] 1.0000000 3.179787e-171 1.589893e-171
[150,] 1.0000000 2.499233e-170 1.249616e-170
[151,] 1.0000000 1.933679e-169 9.668393e-170
[152,] 1.0000000 1.495589e-168 7.477945e-169
[153,] 1.0000000 1.173045e-167 5.865223e-168
[154,] 1.0000000 7.802621e-167 3.901311e-167
[155,] 1.0000000 6.022576e-166 3.011288e-166
[156,] 1.0000000 4.644313e-165 2.322156e-165
[157,] 1.0000000 3.577677e-164 1.788838e-164
[158,] 1.0000000 2.752752e-163 1.376376e-163
[159,] 1.0000000 2.250704e-162 1.125352e-162
[160,] 1.0000000 1.943666e-161 9.718332e-162
[161,] 1.0000000 1.496661e-160 7.483306e-161
[162,] 1.0000000 1.141159e-159 5.705795e-160
[163,] 1.0000000 2.817966e-159 1.408983e-159
[164,] 1.0000000 1.692629e-158 8.463147e-159
[165,] 1.0000000 1.284836e-157 6.424178e-158
[166,] 1.0000000 9.735294e-157 4.867647e-157
[167,] 1.0000000 7.566151e-156 3.783076e-156
[168,] 1.0000000 5.708948e-155 2.854474e-155
[169,] 1.0000000 4.568908e-154 2.284454e-154
[170,] 1.0000000 3.430162e-153 1.715081e-153
[171,] 1.0000000 2.665102e-152 1.332551e-152
[172,] 1.0000000 1.121869e-151 5.609345e-152
[173,] 1.0000000 9.169723e-151 4.584862e-151
[174,] 1.0000000 6.791610e-150 3.395805e-150
[175,] 1.0000000 5.018933e-149 2.509466e-149
[176,] 1.0000000 2.223902e-148 1.111951e-148
[177,] 1.0000000 1.664524e-147 8.322621e-148
[178,] 1.0000000 1.314104e-146 6.570518e-147
[179,] 1.0000000 1.027910e-145 5.139552e-146
[180,] 1.0000000 8.210924e-145 4.105462e-145
[181,] 1.0000000 5.866128e-144 2.933064e-144
[182,] 1.0000000 3.143643e-143 1.571821e-143
[183,] 1.0000000 2.338737e-142 1.169368e-142
[184,] 1.0000000 1.705457e-141 8.527283e-142
[185,] 1.0000000 1.219372e-140 6.096860e-141
[186,] 1.0000000 8.803511e-140 4.401756e-140
[187,] 1.0000000 6.259172e-139 3.129586e-139
[188,] 1.0000000 4.563136e-138 2.281568e-138
[189,] 1.0000000 3.225440e-137 1.612720e-137
[190,] 1.0000000 2.299282e-136 1.149641e-136
[191,] 1.0000000 1.740652e-135 8.703262e-136
[192,] 1.0000000 8.219029e-136 4.109515e-136
[193,] 1.0000000 6.213304e-135 3.106652e-135
[194,] 1.0000000 4.674899e-134 2.337450e-134
[195,] 1.0000000 2.794588e-133 1.397294e-133
[196,] 1.0000000 2.087543e-132 1.043772e-132
[197,] 1.0000000 1.552039e-131 7.760193e-132
[198,] 1.0000000 1.148432e-130 5.742158e-131
[199,] 1.0000000 8.534374e-130 4.267187e-130
[200,] 1.0000000 5.910929e-129 2.955464e-129
[201,] 1.0000000 4.304301e-128 2.152151e-128
[202,] 1.0000000 3.087786e-127 1.543893e-127
[203,] 1.0000000 2.228433e-126 1.114216e-126
[204,] 1.0000000 1.177517e-125 5.887585e-126
[205,] 1.0000000 8.568856e-125 4.284428e-125
[206,] 1.0000000 5.743945e-124 2.871972e-124
[207,] 1.0000000 3.839176e-123 1.919588e-123
[208,] 1.0000000 2.396006e-122 1.198003e-122
[209,] 1.0000000 1.617302e-121 8.086509e-122
[210,] 1.0000000 1.111701e-120 5.558507e-121
[211,] 1.0000000 7.386867e-120 3.693433e-120
[212,] 1.0000000 4.944522e-119 2.472261e-119
[213,] 1.0000000 3.276635e-118 1.638317e-118
[214,] 1.0000000 2.154521e-117 1.077260e-117
[215,] 1.0000000 1.395665e-116 6.978324e-117
[216,] 1.0000000 9.231025e-116 4.615512e-116
[217,] 1.0000000 5.964796e-115 2.982398e-115
[218,] 1.0000000 3.854450e-114 1.927225e-114
[219,] 1.0000000 2.593910e-113 1.296955e-113
[220,] 1.0000000 1.662066e-112 8.310331e-113
[221,] 1.0000000 1.097416e-111 5.487082e-112
[222,] 1.0000000 7.140217e-111 3.570109e-111
[223,] 1.0000000 4.511684e-110 2.255842e-110
[224,] 1.0000000 2.884669e-109 1.442334e-109
[225,] 1.0000000 1.809976e-108 9.049880e-109
[226,] 1.0000000 1.172815e-107 5.864074e-108
[227,] 1.0000000 7.247830e-107 3.623915e-107
[228,] 1.0000000 4.497404e-106 2.248702e-106
[229,] 1.0000000 2.923441e-105 1.461721e-105
[230,] 1.0000000 1.800147e-104 9.000733e-105
[231,] 1.0000000 1.104032e-103 5.520159e-104
[232,] 1.0000000 6.748767e-103 3.374384e-103
[233,] 1.0000000 4.160306e-102 2.080153e-102
[234,] 1.0000000 2.523734e-101 1.261867e-101
[235,] 1.0000000 1.525225e-100 7.626124e-101
[236,] 1.0000000 9.182982e-100 4.591491e-100
[237,] 1.0000000 5.590293e-99 2.795147e-99
[238,] 1.0000000 3.339761e-98 1.669881e-98
[239,] 1.0000000 2.013517e-97 1.006759e-97
[240,] 1.0000000 1.147796e-96 5.738981e-97
[241,] 1.0000000 7.010722e-96 3.505361e-96
[242,] 1.0000000 4.259977e-95 2.129988e-95
[243,] 1.0000000 2.594147e-94 1.297073e-94
[244,] 1.0000000 1.560807e-93 7.804037e-94
[245,] 1.0000000 9.341138e-93 4.670569e-93
[246,] 1.0000000 5.580283e-92 2.790142e-92
[247,] 1.0000000 3.306510e-91 1.653255e-91
[248,] 1.0000000 1.931361e-90 9.656807e-91
[249,] 1.0000000 1.133052e-89 5.665262e-90
[250,] 1.0000000 4.781360e-89 2.390680e-89
[251,] 1.0000000 2.706072e-88 1.353036e-88
[252,] 1.0000000 1.568709e-87 7.843546e-88
[253,] 1.0000000 9.206815e-87 4.603408e-87
[254,] 1.0000000 5.284544e-86 2.642272e-86
[255,] 1.0000000 2.939047e-85 1.469524e-85
[256,] 1.0000000 1.670017e-84 8.350085e-85
[257,] 1.0000000 9.435649e-84 4.717824e-84
[258,] 1.0000000 4.280041e-83 2.140021e-83
[259,] 1.0000000 2.162712e-82 1.081356e-82
[260,] 1.0000000 1.206026e-81 6.030128e-82
[261,] 1.0000000 6.537767e-81 3.268884e-81
[262,] 1.0000000 3.607863e-80 1.803931e-80
[263,] 1.0000000 8.882770e-80 4.441385e-80
[264,] 1.0000000 4.874070e-79 2.437035e-79
[265,] 1.0000000 2.395643e-78 1.197822e-78
[266,] 1.0000000 1.301515e-77 6.507575e-78
[267,] 1.0000000 7.026577e-77 3.513289e-77
[268,] 1.0000000 3.769073e-76 1.884537e-76
[269,] 1.0000000 2.043665e-75 1.021832e-75
[270,] 1.0000000 1.067336e-74 5.336678e-75
[271,] 1.0000000 5.521608e-74 2.760804e-74
[272,] 1.0000000 2.833298e-73 1.416649e-73
[273,] 1.0000000 1.459212e-72 7.296059e-73
[274,] 1.0000000 7.516732e-72 3.758366e-72
[275,] 1.0000000 3.931302e-71 1.965651e-71
[276,] 1.0000000 2.028213e-70 1.014107e-70
[277,] 1.0000000 1.025376e-69 5.126879e-70
[278,] 1.0000000 5.248852e-69 2.624426e-69
[279,] 1.0000000 2.628153e-68 1.314076e-68
[280,] 1.0000000 1.214814e-67 6.074072e-68
[281,] 1.0000000 6.026698e-67 3.013349e-67
[282,] 1.0000000 2.963398e-66 1.481699e-66
[283,] 1.0000000 1.447834e-65 7.239170e-66
[284,] 1.0000000 7.157983e-65 3.578991e-65
[285,] 1.0000000 3.606523e-64 1.803261e-64
[286,] 1.0000000 1.760832e-63 8.804158e-64
[287,] 1.0000000 8.210000e-63 4.105000e-63
[288,] 1.0000000 2.489541e-62 1.244770e-62
[289,] 1.0000000 1.190921e-61 5.954603e-62
[290,] 1.0000000 5.886735e-61 2.943368e-61
[291,] 1.0000000 2.331152e-60 1.165576e-60
[292,] 1.0000000 1.100140e-59 5.500702e-60
[293,] 1.0000000 5.159553e-59 2.579777e-59
[294,] 1.0000000 2.410856e-58 1.205428e-58
[295,] 1.0000000 1.130542e-57 5.652709e-58
[296,] 1.0000000 5.244498e-57 2.622249e-57
[297,] 1.0000000 2.482107e-56 1.241053e-56
[298,] 1.0000000 1.137868e-55 5.689340e-56
[299,] 1.0000000 5.135263e-55 2.567632e-55
[300,] 1.0000000 1.946900e-54 9.734502e-55
[301,] 1.0000000 8.862554e-54 4.431277e-54
[302,] 1.0000000 3.925449e-53 1.962724e-53
[303,] 1.0000000 1.675474e-52 8.377371e-53
[304,] 1.0000000 7.640250e-52 3.820125e-52
[305,] 1.0000000 3.287417e-51 1.643708e-51
[306,] 1.0000000 1.395196e-50 6.975980e-51
[307,] 1.0000000 4.961187e-50 2.480594e-50
[308,] 1.0000000 2.069994e-49 1.034997e-49
[309,] 1.0000000 8.001725e-49 4.000862e-49
[310,] 1.0000000 3.471617e-48 1.735808e-48
[311,] 1.0000000 1.302932e-47 6.514661e-48
[312,] 1.0000000 5.564422e-47 2.782211e-47
[313,] 1.0000000 2.461785e-46 1.230893e-46
[314,] 1.0000000 1.042985e-45 5.214923e-46
[315,] 1.0000000 4.372211e-45 2.186106e-45
[316,] 1.0000000 1.822279e-44 9.111393e-45
[317,] 1.0000000 7.550888e-44 3.775444e-44
[318,] 1.0000000 3.110496e-43 1.555248e-43
[319,] 1.0000000 1.191362e-42 5.956812e-43
[320,] 1.0000000 4.853245e-42 2.426623e-42
[321,] 1.0000000 1.985223e-41 9.926113e-42
[322,] 1.0000000 7.342479e-41 3.671240e-41
[323,] 1.0000000 2.634627e-40 1.317313e-40
[324,] 1.0000000 1.049485e-39 5.247427e-40
[325,] 1.0000000 4.154705e-39 2.077353e-39
[326,] 1.0000000 1.634505e-38 8.172524e-39
[327,] 1.0000000 6.389832e-38 3.194916e-38
[328,] 1.0000000 2.482131e-37 1.241066e-37
[329,] 1.0000000 9.667222e-37 4.833611e-37
[330,] 1.0000000 3.710400e-36 1.855200e-36
[331,] 1.0000000 1.418496e-35 7.092482e-36
[332,] 1.0000000 5.608775e-35 2.804387e-35
[333,] 1.0000000 2.108687e-34 1.054344e-34
[334,] 1.0000000 7.921023e-34 3.960511e-34
[335,] 1.0000000 2.934168e-33 1.467084e-33
[336,] 1.0000000 1.129957e-32 5.649785e-33
[337,] 1.0000000 4.134485e-32 2.067243e-32
[338,] 1.0000000 1.502245e-31 7.511226e-32
[339,] 1.0000000 5.419857e-31 2.709929e-31
[340,] 1.0000000 1.995358e-30 9.976789e-31
[341,] 1.0000000 7.096139e-30 3.548070e-30
[342,] 1.0000000 2.587570e-29 1.293785e-29
[343,] 1.0000000 9.069535e-29 4.534767e-29
[344,] 1.0000000 3.155273e-28 1.577637e-28
[345,] 1.0000000 1.080671e-27 5.403355e-28
[346,] 1.0000000 3.667457e-27 1.833728e-27
[347,] 1.0000000 1.258212e-26 6.291059e-27
[348,] 1.0000000 4.097393e-26 2.048696e-26
[349,] 1.0000000 1.377795e-25 6.888977e-26
[350,] 1.0000000 4.574410e-25 2.287205e-25
[351,] 1.0000000 1.506466e-24 7.532332e-25
[352,] 1.0000000 4.827474e-24 2.413737e-24
[353,] 1.0000000 1.563724e-23 7.818622e-24
[354,] 1.0000000 5.022637e-23 2.511318e-23
[355,] 1.0000000 1.599494e-22 7.997471e-23
[356,] 1.0000000 5.063385e-22 2.531692e-22
[357,] 1.0000000 1.584486e-21 7.922432e-22
[358,] 1.0000000 4.914154e-21 2.457077e-21
[359,] 1.0000000 1.436833e-20 7.184164e-21
[360,] 1.0000000 4.395458e-20 2.197729e-20
[361,] 1.0000000 1.329370e-19 6.646848e-20
[362,] 1.0000000 3.982698e-19 1.991349e-19
[363,] 1.0000000 1.181761e-18 5.908806e-19
[364,] 1.0000000 3.472403e-18 1.736201e-18
[365,] 1.0000000 1.014466e-17 5.072328e-18
[366,] 1.0000000 2.921684e-17 1.460842e-17
[367,] 1.0000000 8.109839e-17 4.054920e-17
[368,] 1.0000000 2.083822e-16 1.041911e-16
[369,] 1.0000000 5.986175e-16 2.993087e-16
[370,] 1.0000000 1.652587e-15 8.262936e-16
[371,] 1.0000000 4.511336e-15 2.255668e-15
[372,] 1.0000000 1.296956e-14 6.484782e-15
[373,] 1.0000000 3.386516e-14 1.693258e-14
[374,] 1.0000000 8.847481e-14 4.423741e-14
[375,] 1.0000000 2.306212e-13 1.153106e-13
[376,] 1.0000000 5.936228e-13 2.968114e-13
[377,] 1.0000000 1.517267e-12 7.586334e-13
[378,] 1.0000000 4.179049e-12 2.089525e-12
[379,] 1.0000000 1.037201e-11 5.186005e-12
[380,] 1.0000000 2.545341e-11 1.272671e-11
[381,] 1.0000000 6.573680e-11 3.286840e-11
[382,] 1.0000000 1.505304e-10 7.526522e-11
[383,] 1.0000000 3.660334e-10 1.830167e-10
[384,] 1.0000000 8.585822e-10 4.292911e-10
[385,] 1.0000000 2.037454e-09 1.018727e-09
[386,] 1.0000000 4.511531e-09 2.255765e-09
[387,] 1.0000000 1.085894e-08 5.429469e-09
[388,] 1.0000000 2.695118e-08 1.347559e-08
[389,] 1.0000000 4.361734e-08 2.180867e-08
[390,] 0.9999999 1.037101e-07 5.185506e-08
[391,] 0.9999999 2.411496e-07 1.205748e-07
[392,] 0.9999997 5.663044e-07 2.831522e-07
[393,] 0.9999998 4.405580e-07 2.202790e-07
[394,] 0.9999995 1.010151e-06 5.050757e-07
[395,] 0.9999988 2.341153e-06 1.170577e-06
[396,] 0.9999975 5.089456e-06 2.544728e-06
[397,] 0.9999951 9.744057e-06 4.872029e-06
[398,] 0.9999896 2.070270e-05 1.035135e-05
[399,] 0.9999776 4.471630e-05 2.235815e-05
[400,] 0.9999523 9.537199e-05 4.768599e-05
[401,] 0.9999134 1.732003e-04 8.660015e-05
[402,] 0.9998414 3.172405e-04 1.586202e-04
[403,] 0.9997054 5.891871e-04 2.945936e-04
[404,] 0.9994584 1.083114e-03 5.415568e-04
[405,] 0.9989082 2.183586e-03 1.091793e-03
[406,] 0.9981264 3.747132e-03 1.873566e-03
[407,] 0.9965397 6.920645e-03 3.460322e-03
[408,] 0.9936916 1.261690e-02 6.308449e-03
[409,] 0.9962002 7.599505e-03 3.799753e-03
[410,] 0.9926660 1.466805e-02 7.334027e-03
[411,] 0.9858898 2.822038e-02 1.411019e-02
[412,] 0.9858736 2.825282e-02 1.412641e-02
[413,] 0.9861077 2.778468e-02 1.389234e-02
[414,] 0.9761110 4.777808e-02 2.388904e-02
[415,] 0.9594301 8.113985e-02 4.056993e-02
[416,] 0.9381104 1.237791e-01 6.188955e-02
[417,] 0.8907581 2.184838e-01 1.092419e-01
[418,] 0.8826826 2.346348e-01 1.173174e-01
[419,] 0.8624195 2.751610e-01 1.375805e-01
[420,] 0.7740979 4.518042e-01 2.259021e-01
> postscript(file="/var/www/html/rcomp/tmp/1hxin1291232858.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/2hxin1291232858.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/3r6hq1291232858.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/4r6hq1291232858.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/5r6hq1291232858.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
2559625.7921 2146005.8721 2073314.5816 -1833510.1393 -133727.1739
6 7 8 9 10
631218.8668 45603.9879 506071.3904 -270524.0886 611450.3878
11 12 13 14 15
397635.4196 232695.3330 -552373.2158 -814949.5497 -546748.1166
16 17 18 19 20
570626.1725 -28097.3485 -554963.3777 -287336.7841 -968495.6350
21 22 23 24 25
263815.8794 144146.5750 31189.5356 -1109804.1420 -287950.1460
26 27 28 29 30
114890.9839 -630335.0160 -333838.9868 409654.5342 -384864.6377
31 32 33 34 35
295744.3278 -6272.7728 288690.7049 136548.4730 360943.0845
36 37 38 39 40
-753672.9841 92712.8973 74215.0680 -85424.2798 -583164.7103
41 42 43 44 45
-795895.0567 -125041.0161 -109757.8415 -227260.7256 -10876.0454
46 47 48 49 50
-11466.5563 -120249.3384 -486531.9307 -60347.3817 -212384.2772
51 52 53 54 55
-135204.0912 -343612.4237 -204412.9027 -23748.9613 -534415.8734
56 57 58 59 60
-154377.4505 395256.7801 13755.0272 -108168.4365 -256859.5116
61 62 63 64 65
27900.3571 -222015.2339 25690.4437 -246149.5843 -107119.7842
66 67 68 69 70
-944371.6210 -34820.2632 -72626.7588 -168353.9868 -219660.5116
71 72 73 74 75
-289031.3384 -15562.9295 -70530.5817 55280.4323 17474.5342
76 77 78 79 80
-120395.9767 12358.1138 -42222.7996 -141393.5524 2221.8348
81 82 83 84 85
-8334.8862 -137279.3091 -19922.3499 5844.7368 81479.7368
86 87 88 89 90
-214773.5409 -200058.1167 66717.9686 -115293.7409 94438.8527
91 92 93 94 95
7781.3024 -74209.2951 -69512.4071 -40985.0200 -56754.5091
96 97 98 99 100
-33446.1359 158654.4323 -45436.2085 90471.2731 -2669.7015
101 102 103 104 105
-106822.2518 -77955.7269 -167317.9002 47878.5482 16729.2298
106 107 108 109 110
-86058.8862 -361433.2225 160430.4870 -110559.9728 -117557.9767
111 112 113 114 115
83002.9253 -264384.4836 -127512.6250 -48101.1652 70276.6463
116 117 118 119 120
-37832.5817 66548.5482 138319.5049 -114872.9333 -5926.9588
121 122 123 124 125
-249101.1893 51033.0845 -353694.5078 91610.4437 -1608.6289
126 127 128 129 130
-163052.3244 83290.4756 39812.3164 -12085.1906 83290.4756
131 132 133 134 135
52421.6641 30780.5775 83290.4756 -4505.8429 62933.0119
136 137 138 139 140
65867.0119 44182.3164 78534.3597 65771.0119 83290.4756
141 142 143 144 145
59311.5482 83290.4756 19584.2004 66692.0119 33747.7368
146 147 148 149 150
42420.4323 83290.4756 53055.6641 84633.2438 83787.4756
151 152 153 154 155
29718.6641 -122233.8568 61974.8960 35833.6641 74629.2438
156 157 158 159 160
83290.4756 83290.4756 74620.2438 83019.8960 83290.4756
161 162 163 164 165
83290.4756 83290.4756 83290.4756 44336.6641 -32217.3065
166 167 168 169 170
78529.3597 83290.4756 73886.3890 -124990.9002 83290.4756
171 172 173 174 175
83290.4756 67533.1278 83290.4756 39138.4323 83290.4756
176 177 178 179 180
49352.5482 -400554.6237 33982.3418 83290.4756 83290.4756
181 182 183 184 185
46246.2731 980.0667 -24378.4225 -52990.5130 11296.7622
186 187 188 189 190
4122.0667 -28137.3244 55663.7801 38381.8960 83290.4756
191 192 193 194 195
56813.6641 83290.4756 53434.6641 83290.4756 70809.1278
196 197 198 199 200
33186.3418 -251968.9753 33186.3418 33186.3418 -93605.8607
201 202 203 204 205
33186.3418 33186.3418 33186.3418 -32389.0493 -7872.0493
206 207 208 209 210
12645.7622 27854.1100 33186.3418 -49421.6722 11964.9940
211 212 213 214 215
83290.4756 83290.4756 51075.2004 14264.3164 51079.8960
216 217 218 219 220
74764.2438 58019.7801 55804.6641 74764.2438 -29125.7702
221 222 223 224 225
51777.1278 83290.4756 79001.3597 -5190.8429 79223.3597
226 227 228 229 230
-10956.0747 44093.4323 83290.4756 71723.2438 83290.4756
231 232 233 234 235
36768.6641 50918.4323 83290.4756 15719.9686 83290.4756
236 237 238 239 240
83397.4756 83290.4756 58124.7801 83290.4756 83290.4756
241 242 243 244 245
83290.4756 68008.1278 83290.4756 53705.6641 83769.2438
246 247 248 249 250
33186.3418 33186.3418 12952.7368 33186.3418 33186.3418
251 252 253 254 255
26244.4323 33186.3418 70969.3597 -89711.9295 -221499.1320
256 257 258 259 260
83290.4756 47628.6641 -44464.5384 33186.3418 83290.4756
261 262 263 264 265
33186.3418 33186.3418 56418.9686 57234.3164 33186.3418
266 267 268 269 270
83290.4756 33186.3418 110722.4323 33186.3418 -11510.5384
271 272 273 274 275
33186.3418 33186.3418 33186.3418 -2887.7269 83290.4756
276 277 278 279 280
36769.0845 -9743.1906 83290.4756 -113983.0200 44475.5482
281 282 283 284 285
33186.3418 83290.4756 -29439.5384 83290.4756 60865.4323
286 287 288 289 290
83290.4756 15700.7622 73219.1278 33186.3418 5243.9686
291 292 293 294 295
33186.3418 41929.0845 52870.2731 83290.4756 27688.3164
296 297 298 299 300
20406.2985 83290.4756 33186.3418 83290.4756 44922.4323
301 302 303 304 305
13917.8781 64811.3597 79027.3597 33186.3418 52651.8527
306 307 308 309 310
32147.0845 33186.3418 -46545.8607 5273.5049 33186.3418
311 312 313 314 315
33186.3418 -21898.2339 33186.3418 -19846.2811 49461.5482
316 317 318 319 320
90303.0119 82874.4756 -20210.8429 61966.8960 83290.4756
321 322 323 324 325
83290.4756 83290.4756 83290.4756 33186.3418 83290.4756
326 327 328 329 330
76149.3597 33186.3418 33186.3418 83290.4756 83290.4756
331 332 333 334 335
83290.4756 83290.4756 83290.4756 64587.0119 79035.3597
336 337 338 339 340
66229.0119 -21212.8429 83290.4756 77527.3597 42977.3164
341 342 343 344 345
12333.0845 83290.4756 83290.4756 83290.4756 20745.9686
346 347 348 349 350
83290.4756 58798.1278 83290.4756 83290.4756 44648.3164
351 352 353 354 355
73811.1278 52117.6641 50978.3164 77857.3597 83290.4756
356 357 358 359 360
83290.4756 26341.7368 83290.4756 83290.4756 83290.4756
361 362 363 364 365
66256.0119 83290.4756 83290.4756 -93224.9040 45928.4323
366 367 368 369 370
83290.4756 83290.4756 83290.4756 83290.4756 53621.6641
371 372 373 374 375
83290.4756 61429.6641 74707.6641 62545.1278 83290.4756
376 377 378 379 380
83290.4756 16386.2004 60408.6641 45143.3164 83290.4756
381 382 383 384 385
83290.4756 8465.3890 -11375.1473 83290.4756 34083.0845
386 387 388 389 390
-117594.8862 82192.1278 51369.7801 69331.1278 41735.5482
391 392 393 394 395
47606.3164 34467.7801 -2918.4518 109986.7801 -36651.7015
396 397 398 399 400
17392.5482 -149384.6722 -862585.3637 34329.6641 -226360.4110
401 402 403 404 405
59373.1278 67310.5482 39250.7801 -92729.0926 14966.6641
406 407 408 409 410
48023.2004 38062.2004 21308.0845 27638.4323 -13017.8722
411 412 413 414 415
-165665.1320 -84756.0926 -76387.3970 29629.6323 -113579.7996
416 417 418 419 420
-36290.7562 77017.8094 -76979.9767 -59093.3065 10188.3597
421 422 423 424 425
17115.3597 -306663.1027 -655350.5409 -359494.8429 -44638.8429
426 427 428 429 430
-416407.1639 -278860.9588 -966014.8096 -114972.2339 -226086.4938
431
-533414.5231
> postscript(file="/var/www/html/rcomp/tmp/62fzt1291232858.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 2559625.7921 NA
1 2146005.8721 2559625.7921
2 2073314.5816 2146005.8721
3 -1833510.1393 2073314.5816
4 -133727.1739 -1833510.1393
5 631218.8668 -133727.1739
6 45603.9879 631218.8668
7 506071.3904 45603.9879
8 -270524.0886 506071.3904
9 611450.3878 -270524.0886
10 397635.4196 611450.3878
11 232695.3330 397635.4196
12 -552373.2158 232695.3330
13 -814949.5497 -552373.2158
14 -546748.1166 -814949.5497
15 570626.1725 -546748.1166
16 -28097.3485 570626.1725
17 -554963.3777 -28097.3485
18 -287336.7841 -554963.3777
19 -968495.6350 -287336.7841
20 263815.8794 -968495.6350
21 144146.5750 263815.8794
22 31189.5356 144146.5750
23 -1109804.1420 31189.5356
24 -287950.1460 -1109804.1420
25 114890.9839 -287950.1460
26 -630335.0160 114890.9839
27 -333838.9868 -630335.0160
28 409654.5342 -333838.9868
29 -384864.6377 409654.5342
30 295744.3278 -384864.6377
31 -6272.7728 295744.3278
32 288690.7049 -6272.7728
33 136548.4730 288690.7049
34 360943.0845 136548.4730
35 -753672.9841 360943.0845
36 92712.8973 -753672.9841
37 74215.0680 92712.8973
38 -85424.2798 74215.0680
39 -583164.7103 -85424.2798
40 -795895.0567 -583164.7103
41 -125041.0161 -795895.0567
42 -109757.8415 -125041.0161
43 -227260.7256 -109757.8415
44 -10876.0454 -227260.7256
45 -11466.5563 -10876.0454
46 -120249.3384 -11466.5563
47 -486531.9307 -120249.3384
48 -60347.3817 -486531.9307
49 -212384.2772 -60347.3817
50 -135204.0912 -212384.2772
51 -343612.4237 -135204.0912
52 -204412.9027 -343612.4237
53 -23748.9613 -204412.9027
54 -534415.8734 -23748.9613
55 -154377.4505 -534415.8734
56 395256.7801 -154377.4505
57 13755.0272 395256.7801
58 -108168.4365 13755.0272
59 -256859.5116 -108168.4365
60 27900.3571 -256859.5116
61 -222015.2339 27900.3571
62 25690.4437 -222015.2339
63 -246149.5843 25690.4437
64 -107119.7842 -246149.5843
65 -944371.6210 -107119.7842
66 -34820.2632 -944371.6210
67 -72626.7588 -34820.2632
68 -168353.9868 -72626.7588
69 -219660.5116 -168353.9868
70 -289031.3384 -219660.5116
71 -15562.9295 -289031.3384
72 -70530.5817 -15562.9295
73 55280.4323 -70530.5817
74 17474.5342 55280.4323
75 -120395.9767 17474.5342
76 12358.1138 -120395.9767
77 -42222.7996 12358.1138
78 -141393.5524 -42222.7996
79 2221.8348 -141393.5524
80 -8334.8862 2221.8348
81 -137279.3091 -8334.8862
82 -19922.3499 -137279.3091
83 5844.7368 -19922.3499
84 81479.7368 5844.7368
85 -214773.5409 81479.7368
86 -200058.1167 -214773.5409
87 66717.9686 -200058.1167
88 -115293.7409 66717.9686
89 94438.8527 -115293.7409
90 7781.3024 94438.8527
91 -74209.2951 7781.3024
92 -69512.4071 -74209.2951
93 -40985.0200 -69512.4071
94 -56754.5091 -40985.0200
95 -33446.1359 -56754.5091
96 158654.4323 -33446.1359
97 -45436.2085 158654.4323
98 90471.2731 -45436.2085
99 -2669.7015 90471.2731
100 -106822.2518 -2669.7015
101 -77955.7269 -106822.2518
102 -167317.9002 -77955.7269
103 47878.5482 -167317.9002
104 16729.2298 47878.5482
105 -86058.8862 16729.2298
106 -361433.2225 -86058.8862
107 160430.4870 -361433.2225
108 -110559.9728 160430.4870
109 -117557.9767 -110559.9728
110 83002.9253 -117557.9767
111 -264384.4836 83002.9253
112 -127512.6250 -264384.4836
113 -48101.1652 -127512.6250
114 70276.6463 -48101.1652
115 -37832.5817 70276.6463
116 66548.5482 -37832.5817
117 138319.5049 66548.5482
118 -114872.9333 138319.5049
119 -5926.9588 -114872.9333
120 -249101.1893 -5926.9588
121 51033.0845 -249101.1893
122 -353694.5078 51033.0845
123 91610.4437 -353694.5078
124 -1608.6289 91610.4437
125 -163052.3244 -1608.6289
126 83290.4756 -163052.3244
127 39812.3164 83290.4756
128 -12085.1906 39812.3164
129 83290.4756 -12085.1906
130 52421.6641 83290.4756
131 30780.5775 52421.6641
132 83290.4756 30780.5775
133 -4505.8429 83290.4756
134 62933.0119 -4505.8429
135 65867.0119 62933.0119
136 44182.3164 65867.0119
137 78534.3597 44182.3164
138 65771.0119 78534.3597
139 83290.4756 65771.0119
140 59311.5482 83290.4756
141 83290.4756 59311.5482
142 19584.2004 83290.4756
143 66692.0119 19584.2004
144 33747.7368 66692.0119
145 42420.4323 33747.7368
146 83290.4756 42420.4323
147 53055.6641 83290.4756
148 84633.2438 53055.6641
149 83787.4756 84633.2438
150 29718.6641 83787.4756
151 -122233.8568 29718.6641
152 61974.8960 -122233.8568
153 35833.6641 61974.8960
154 74629.2438 35833.6641
155 83290.4756 74629.2438
156 83290.4756 83290.4756
157 74620.2438 83290.4756
158 83019.8960 74620.2438
159 83290.4756 83019.8960
160 83290.4756 83290.4756
161 83290.4756 83290.4756
162 83290.4756 83290.4756
163 44336.6641 83290.4756
164 -32217.3065 44336.6641
165 78529.3597 -32217.3065
166 83290.4756 78529.3597
167 73886.3890 83290.4756
168 -124990.9002 73886.3890
169 83290.4756 -124990.9002
170 83290.4756 83290.4756
171 67533.1278 83290.4756
172 83290.4756 67533.1278
173 39138.4323 83290.4756
174 83290.4756 39138.4323
175 49352.5482 83290.4756
176 -400554.6237 49352.5482
177 33982.3418 -400554.6237
178 83290.4756 33982.3418
179 83290.4756 83290.4756
180 46246.2731 83290.4756
181 980.0667 46246.2731
182 -24378.4225 980.0667
183 -52990.5130 -24378.4225
184 11296.7622 -52990.5130
185 4122.0667 11296.7622
186 -28137.3244 4122.0667
187 55663.7801 -28137.3244
188 38381.8960 55663.7801
189 83290.4756 38381.8960
190 56813.6641 83290.4756
191 83290.4756 56813.6641
192 53434.6641 83290.4756
193 83290.4756 53434.6641
194 70809.1278 83290.4756
195 33186.3418 70809.1278
196 -251968.9753 33186.3418
197 33186.3418 -251968.9753
198 33186.3418 33186.3418
199 -93605.8607 33186.3418
200 33186.3418 -93605.8607
201 33186.3418 33186.3418
202 33186.3418 33186.3418
203 -32389.0493 33186.3418
204 -7872.0493 -32389.0493
205 12645.7622 -7872.0493
206 27854.1100 12645.7622
207 33186.3418 27854.1100
208 -49421.6722 33186.3418
209 11964.9940 -49421.6722
210 83290.4756 11964.9940
211 83290.4756 83290.4756
212 51075.2004 83290.4756
213 14264.3164 51075.2004
214 51079.8960 14264.3164
215 74764.2438 51079.8960
216 58019.7801 74764.2438
217 55804.6641 58019.7801
218 74764.2438 55804.6641
219 -29125.7702 74764.2438
220 51777.1278 -29125.7702
221 83290.4756 51777.1278
222 79001.3597 83290.4756
223 -5190.8429 79001.3597
224 79223.3597 -5190.8429
225 -10956.0747 79223.3597
226 44093.4323 -10956.0747
227 83290.4756 44093.4323
228 71723.2438 83290.4756
229 83290.4756 71723.2438
230 36768.6641 83290.4756
231 50918.4323 36768.6641
232 83290.4756 50918.4323
233 15719.9686 83290.4756
234 83290.4756 15719.9686
235 83397.4756 83290.4756
236 83290.4756 83397.4756
237 58124.7801 83290.4756
238 83290.4756 58124.7801
239 83290.4756 83290.4756
240 83290.4756 83290.4756
241 68008.1278 83290.4756
242 83290.4756 68008.1278
243 53705.6641 83290.4756
244 83769.2438 53705.6641
245 33186.3418 83769.2438
246 33186.3418 33186.3418
247 12952.7368 33186.3418
248 33186.3418 12952.7368
249 33186.3418 33186.3418
250 26244.4323 33186.3418
251 33186.3418 26244.4323
252 70969.3597 33186.3418
253 -89711.9295 70969.3597
254 -221499.1320 -89711.9295
255 83290.4756 -221499.1320
256 47628.6641 83290.4756
257 -44464.5384 47628.6641
258 33186.3418 -44464.5384
259 83290.4756 33186.3418
260 33186.3418 83290.4756
261 33186.3418 33186.3418
262 56418.9686 33186.3418
263 57234.3164 56418.9686
264 33186.3418 57234.3164
265 83290.4756 33186.3418
266 33186.3418 83290.4756
267 110722.4323 33186.3418
268 33186.3418 110722.4323
269 -11510.5384 33186.3418
270 33186.3418 -11510.5384
271 33186.3418 33186.3418
272 33186.3418 33186.3418
273 -2887.7269 33186.3418
274 83290.4756 -2887.7269
275 36769.0845 83290.4756
276 -9743.1906 36769.0845
277 83290.4756 -9743.1906
278 -113983.0200 83290.4756
279 44475.5482 -113983.0200
280 33186.3418 44475.5482
281 83290.4756 33186.3418
282 -29439.5384 83290.4756
283 83290.4756 -29439.5384
284 60865.4323 83290.4756
285 83290.4756 60865.4323
286 15700.7622 83290.4756
287 73219.1278 15700.7622
288 33186.3418 73219.1278
289 5243.9686 33186.3418
290 33186.3418 5243.9686
291 41929.0845 33186.3418
292 52870.2731 41929.0845
293 83290.4756 52870.2731
294 27688.3164 83290.4756
295 20406.2985 27688.3164
296 83290.4756 20406.2985
297 33186.3418 83290.4756
298 83290.4756 33186.3418
299 44922.4323 83290.4756
300 13917.8781 44922.4323
301 64811.3597 13917.8781
302 79027.3597 64811.3597
303 33186.3418 79027.3597
304 52651.8527 33186.3418
305 32147.0845 52651.8527
306 33186.3418 32147.0845
307 -46545.8607 33186.3418
308 5273.5049 -46545.8607
309 33186.3418 5273.5049
310 33186.3418 33186.3418
311 -21898.2339 33186.3418
312 33186.3418 -21898.2339
313 -19846.2811 33186.3418
314 49461.5482 -19846.2811
315 90303.0119 49461.5482
316 82874.4756 90303.0119
317 -20210.8429 82874.4756
318 61966.8960 -20210.8429
319 83290.4756 61966.8960
320 83290.4756 83290.4756
321 83290.4756 83290.4756
322 83290.4756 83290.4756
323 33186.3418 83290.4756
324 83290.4756 33186.3418
325 76149.3597 83290.4756
326 33186.3418 76149.3597
327 33186.3418 33186.3418
328 83290.4756 33186.3418
329 83290.4756 83290.4756
330 83290.4756 83290.4756
331 83290.4756 83290.4756
332 83290.4756 83290.4756
333 64587.0119 83290.4756
334 79035.3597 64587.0119
335 66229.0119 79035.3597
336 -21212.8429 66229.0119
337 83290.4756 -21212.8429
338 77527.3597 83290.4756
339 42977.3164 77527.3597
340 12333.0845 42977.3164
341 83290.4756 12333.0845
342 83290.4756 83290.4756
343 83290.4756 83290.4756
344 20745.9686 83290.4756
345 83290.4756 20745.9686
346 58798.1278 83290.4756
347 83290.4756 58798.1278
348 83290.4756 83290.4756
349 44648.3164 83290.4756
350 73811.1278 44648.3164
351 52117.6641 73811.1278
352 50978.3164 52117.6641
353 77857.3597 50978.3164
354 83290.4756 77857.3597
355 83290.4756 83290.4756
356 26341.7368 83290.4756
357 83290.4756 26341.7368
358 83290.4756 83290.4756
359 83290.4756 83290.4756
360 66256.0119 83290.4756
361 83290.4756 66256.0119
362 83290.4756 83290.4756
363 -93224.9040 83290.4756
364 45928.4323 -93224.9040
365 83290.4756 45928.4323
366 83290.4756 83290.4756
367 83290.4756 83290.4756
368 83290.4756 83290.4756
369 53621.6641 83290.4756
370 83290.4756 53621.6641
371 61429.6641 83290.4756
372 74707.6641 61429.6641
373 62545.1278 74707.6641
374 83290.4756 62545.1278
375 83290.4756 83290.4756
376 16386.2004 83290.4756
377 60408.6641 16386.2004
378 45143.3164 60408.6641
379 83290.4756 45143.3164
380 83290.4756 83290.4756
381 8465.3890 83290.4756
382 -11375.1473 8465.3890
383 83290.4756 -11375.1473
384 34083.0845 83290.4756
385 -117594.8862 34083.0845
386 82192.1278 -117594.8862
387 51369.7801 82192.1278
388 69331.1278 51369.7801
389 41735.5482 69331.1278
390 47606.3164 41735.5482
391 34467.7801 47606.3164
392 -2918.4518 34467.7801
393 109986.7801 -2918.4518
394 -36651.7015 109986.7801
395 17392.5482 -36651.7015
396 -149384.6722 17392.5482
397 -862585.3637 -149384.6722
398 34329.6641 -862585.3637
399 -226360.4110 34329.6641
400 59373.1278 -226360.4110
401 67310.5482 59373.1278
402 39250.7801 67310.5482
403 -92729.0926 39250.7801
404 14966.6641 -92729.0926
405 48023.2004 14966.6641
406 38062.2004 48023.2004
407 21308.0845 38062.2004
408 27638.4323 21308.0845
409 -13017.8722 27638.4323
410 -165665.1320 -13017.8722
411 -84756.0926 -165665.1320
412 -76387.3970 -84756.0926
413 29629.6323 -76387.3970
414 -113579.7996 29629.6323
415 -36290.7562 -113579.7996
416 77017.8094 -36290.7562
417 -76979.9767 77017.8094
418 -59093.3065 -76979.9767
419 10188.3597 -59093.3065
420 17115.3597 10188.3597
421 -306663.1027 17115.3597
422 -655350.5409 -306663.1027
423 -359494.8429 -655350.5409
424 -44638.8429 -359494.8429
425 -416407.1639 -44638.8429
426 -278860.9588 -416407.1639
427 -966014.8096 -278860.9588
428 -114972.2339 -966014.8096
429 -226086.4938 -114972.2339
430 -533414.5231 -226086.4938
431 NA -533414.5231
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2146005.8721 2559625.7921
[2,] 2073314.5816 2146005.8721
[3,] -1833510.1393 2073314.5816
[4,] -133727.1739 -1833510.1393
[5,] 631218.8668 -133727.1739
[6,] 45603.9879 631218.8668
[7,] 506071.3904 45603.9879
[8,] -270524.0886 506071.3904
[9,] 611450.3878 -270524.0886
[10,] 397635.4196 611450.3878
[11,] 232695.3330 397635.4196
[12,] -552373.2158 232695.3330
[13,] -814949.5497 -552373.2158
[14,] -546748.1166 -814949.5497
[15,] 570626.1725 -546748.1166
[16,] -28097.3485 570626.1725
[17,] -554963.3777 -28097.3485
[18,] -287336.7841 -554963.3777
[19,] -968495.6350 -287336.7841
[20,] 263815.8794 -968495.6350
[21,] 144146.5750 263815.8794
[22,] 31189.5356 144146.5750
[23,] -1109804.1420 31189.5356
[24,] -287950.1460 -1109804.1420
[25,] 114890.9839 -287950.1460
[26,] -630335.0160 114890.9839
[27,] -333838.9868 -630335.0160
[28,] 409654.5342 -333838.9868
[29,] -384864.6377 409654.5342
[30,] 295744.3278 -384864.6377
[31,] -6272.7728 295744.3278
[32,] 288690.7049 -6272.7728
[33,] 136548.4730 288690.7049
[34,] 360943.0845 136548.4730
[35,] -753672.9841 360943.0845
[36,] 92712.8973 -753672.9841
[37,] 74215.0680 92712.8973
[38,] -85424.2798 74215.0680
[39,] -583164.7103 -85424.2798
[40,] -795895.0567 -583164.7103
[41,] -125041.0161 -795895.0567
[42,] -109757.8415 -125041.0161
[43,] -227260.7256 -109757.8415
[44,] -10876.0454 -227260.7256
[45,] -11466.5563 -10876.0454
[46,] -120249.3384 -11466.5563
[47,] -486531.9307 -120249.3384
[48,] -60347.3817 -486531.9307
[49,] -212384.2772 -60347.3817
[50,] -135204.0912 -212384.2772
[51,] -343612.4237 -135204.0912
[52,] -204412.9027 -343612.4237
[53,] -23748.9613 -204412.9027
[54,] -534415.8734 -23748.9613
[55,] -154377.4505 -534415.8734
[56,] 395256.7801 -154377.4505
[57,] 13755.0272 395256.7801
[58,] -108168.4365 13755.0272
[59,] -256859.5116 -108168.4365
[60,] 27900.3571 -256859.5116
[61,] -222015.2339 27900.3571
[62,] 25690.4437 -222015.2339
[63,] -246149.5843 25690.4437
[64,] -107119.7842 -246149.5843
[65,] -944371.6210 -107119.7842
[66,] -34820.2632 -944371.6210
[67,] -72626.7588 -34820.2632
[68,] -168353.9868 -72626.7588
[69,] -219660.5116 -168353.9868
[70,] -289031.3384 -219660.5116
[71,] -15562.9295 -289031.3384
[72,] -70530.5817 -15562.9295
[73,] 55280.4323 -70530.5817
[74,] 17474.5342 55280.4323
[75,] -120395.9767 17474.5342
[76,] 12358.1138 -120395.9767
[77,] -42222.7996 12358.1138
[78,] -141393.5524 -42222.7996
[79,] 2221.8348 -141393.5524
[80,] -8334.8862 2221.8348
[81,] -137279.3091 -8334.8862
[82,] -19922.3499 -137279.3091
[83,] 5844.7368 -19922.3499
[84,] 81479.7368 5844.7368
[85,] -214773.5409 81479.7368
[86,] -200058.1167 -214773.5409
[87,] 66717.9686 -200058.1167
[88,] -115293.7409 66717.9686
[89,] 94438.8527 -115293.7409
[90,] 7781.3024 94438.8527
[91,] -74209.2951 7781.3024
[92,] -69512.4071 -74209.2951
[93,] -40985.0200 -69512.4071
[94,] -56754.5091 -40985.0200
[95,] -33446.1359 -56754.5091
[96,] 158654.4323 -33446.1359
[97,] -45436.2085 158654.4323
[98,] 90471.2731 -45436.2085
[99,] -2669.7015 90471.2731
[100,] -106822.2518 -2669.7015
[101,] -77955.7269 -106822.2518
[102,] -167317.9002 -77955.7269
[103,] 47878.5482 -167317.9002
[104,] 16729.2298 47878.5482
[105,] -86058.8862 16729.2298
[106,] -361433.2225 -86058.8862
[107,] 160430.4870 -361433.2225
[108,] -110559.9728 160430.4870
[109,] -117557.9767 -110559.9728
[110,] 83002.9253 -117557.9767
[111,] -264384.4836 83002.9253
[112,] -127512.6250 -264384.4836
[113,] -48101.1652 -127512.6250
[114,] 70276.6463 -48101.1652
[115,] -37832.5817 70276.6463
[116,] 66548.5482 -37832.5817
[117,] 138319.5049 66548.5482
[118,] -114872.9333 138319.5049
[119,] -5926.9588 -114872.9333
[120,] -249101.1893 -5926.9588
[121,] 51033.0845 -249101.1893
[122,] -353694.5078 51033.0845
[123,] 91610.4437 -353694.5078
[124,] -1608.6289 91610.4437
[125,] -163052.3244 -1608.6289
[126,] 83290.4756 -163052.3244
[127,] 39812.3164 83290.4756
[128,] -12085.1906 39812.3164
[129,] 83290.4756 -12085.1906
[130,] 52421.6641 83290.4756
[131,] 30780.5775 52421.6641
[132,] 83290.4756 30780.5775
[133,] -4505.8429 83290.4756
[134,] 62933.0119 -4505.8429
[135,] 65867.0119 62933.0119
[136,] 44182.3164 65867.0119
[137,] 78534.3597 44182.3164
[138,] 65771.0119 78534.3597
[139,] 83290.4756 65771.0119
[140,] 59311.5482 83290.4756
[141,] 83290.4756 59311.5482
[142,] 19584.2004 83290.4756
[143,] 66692.0119 19584.2004
[144,] 33747.7368 66692.0119
[145,] 42420.4323 33747.7368
[146,] 83290.4756 42420.4323
[147,] 53055.6641 83290.4756
[148,] 84633.2438 53055.6641
[149,] 83787.4756 84633.2438
[150,] 29718.6641 83787.4756
[151,] -122233.8568 29718.6641
[152,] 61974.8960 -122233.8568
[153,] 35833.6641 61974.8960
[154,] 74629.2438 35833.6641
[155,] 83290.4756 74629.2438
[156,] 83290.4756 83290.4756
[157,] 74620.2438 83290.4756
[158,] 83019.8960 74620.2438
[159,] 83290.4756 83019.8960
[160,] 83290.4756 83290.4756
[161,] 83290.4756 83290.4756
[162,] 83290.4756 83290.4756
[163,] 44336.6641 83290.4756
[164,] -32217.3065 44336.6641
[165,] 78529.3597 -32217.3065
[166,] 83290.4756 78529.3597
[167,] 73886.3890 83290.4756
[168,] -124990.9002 73886.3890
[169,] 83290.4756 -124990.9002
[170,] 83290.4756 83290.4756
[171,] 67533.1278 83290.4756
[172,] 83290.4756 67533.1278
[173,] 39138.4323 83290.4756
[174,] 83290.4756 39138.4323
[175,] 49352.5482 83290.4756
[176,] -400554.6237 49352.5482
[177,] 33982.3418 -400554.6237
[178,] 83290.4756 33982.3418
[179,] 83290.4756 83290.4756
[180,] 46246.2731 83290.4756
[181,] 980.0667 46246.2731
[182,] -24378.4225 980.0667
[183,] -52990.5130 -24378.4225
[184,] 11296.7622 -52990.5130
[185,] 4122.0667 11296.7622
[186,] -28137.3244 4122.0667
[187,] 55663.7801 -28137.3244
[188,] 38381.8960 55663.7801
[189,] 83290.4756 38381.8960
[190,] 56813.6641 83290.4756
[191,] 83290.4756 56813.6641
[192,] 53434.6641 83290.4756
[193,] 83290.4756 53434.6641
[194,] 70809.1278 83290.4756
[195,] 33186.3418 70809.1278
[196,] -251968.9753 33186.3418
[197,] 33186.3418 -251968.9753
[198,] 33186.3418 33186.3418
[199,] -93605.8607 33186.3418
[200,] 33186.3418 -93605.8607
[201,] 33186.3418 33186.3418
[202,] 33186.3418 33186.3418
[203,] -32389.0493 33186.3418
[204,] -7872.0493 -32389.0493
[205,] 12645.7622 -7872.0493
[206,] 27854.1100 12645.7622
[207,] 33186.3418 27854.1100
[208,] -49421.6722 33186.3418
[209,] 11964.9940 -49421.6722
[210,] 83290.4756 11964.9940
[211,] 83290.4756 83290.4756
[212,] 51075.2004 83290.4756
[213,] 14264.3164 51075.2004
[214,] 51079.8960 14264.3164
[215,] 74764.2438 51079.8960
[216,] 58019.7801 74764.2438
[217,] 55804.6641 58019.7801
[218,] 74764.2438 55804.6641
[219,] -29125.7702 74764.2438
[220,] 51777.1278 -29125.7702
[221,] 83290.4756 51777.1278
[222,] 79001.3597 83290.4756
[223,] -5190.8429 79001.3597
[224,] 79223.3597 -5190.8429
[225,] -10956.0747 79223.3597
[226,] 44093.4323 -10956.0747
[227,] 83290.4756 44093.4323
[228,] 71723.2438 83290.4756
[229,] 83290.4756 71723.2438
[230,] 36768.6641 83290.4756
[231,] 50918.4323 36768.6641
[232,] 83290.4756 50918.4323
[233,] 15719.9686 83290.4756
[234,] 83290.4756 15719.9686
[235,] 83397.4756 83290.4756
[236,] 83290.4756 83397.4756
[237,] 58124.7801 83290.4756
[238,] 83290.4756 58124.7801
[239,] 83290.4756 83290.4756
[240,] 83290.4756 83290.4756
[241,] 68008.1278 83290.4756
[242,] 83290.4756 68008.1278
[243,] 53705.6641 83290.4756
[244,] 83769.2438 53705.6641
[245,] 33186.3418 83769.2438
[246,] 33186.3418 33186.3418
[247,] 12952.7368 33186.3418
[248,] 33186.3418 12952.7368
[249,] 33186.3418 33186.3418
[250,] 26244.4323 33186.3418
[251,] 33186.3418 26244.4323
[252,] 70969.3597 33186.3418
[253,] -89711.9295 70969.3597
[254,] -221499.1320 -89711.9295
[255,] 83290.4756 -221499.1320
[256,] 47628.6641 83290.4756
[257,] -44464.5384 47628.6641
[258,] 33186.3418 -44464.5384
[259,] 83290.4756 33186.3418
[260,] 33186.3418 83290.4756
[261,] 33186.3418 33186.3418
[262,] 56418.9686 33186.3418
[263,] 57234.3164 56418.9686
[264,] 33186.3418 57234.3164
[265,] 83290.4756 33186.3418
[266,] 33186.3418 83290.4756
[267,] 110722.4323 33186.3418
[268,] 33186.3418 110722.4323
[269,] -11510.5384 33186.3418
[270,] 33186.3418 -11510.5384
[271,] 33186.3418 33186.3418
[272,] 33186.3418 33186.3418
[273,] -2887.7269 33186.3418
[274,] 83290.4756 -2887.7269
[275,] 36769.0845 83290.4756
[276,] -9743.1906 36769.0845
[277,] 83290.4756 -9743.1906
[278,] -113983.0200 83290.4756
[279,] 44475.5482 -113983.0200
[280,] 33186.3418 44475.5482
[281,] 83290.4756 33186.3418
[282,] -29439.5384 83290.4756
[283,] 83290.4756 -29439.5384
[284,] 60865.4323 83290.4756
[285,] 83290.4756 60865.4323
[286,] 15700.7622 83290.4756
[287,] 73219.1278 15700.7622
[288,] 33186.3418 73219.1278
[289,] 5243.9686 33186.3418
[290,] 33186.3418 5243.9686
[291,] 41929.0845 33186.3418
[292,] 52870.2731 41929.0845
[293,] 83290.4756 52870.2731
[294,] 27688.3164 83290.4756
[295,] 20406.2985 27688.3164
[296,] 83290.4756 20406.2985
[297,] 33186.3418 83290.4756
[298,] 83290.4756 33186.3418
[299,] 44922.4323 83290.4756
[300,] 13917.8781 44922.4323
[301,] 64811.3597 13917.8781
[302,] 79027.3597 64811.3597
[303,] 33186.3418 79027.3597
[304,] 52651.8527 33186.3418
[305,] 32147.0845 52651.8527
[306,] 33186.3418 32147.0845
[307,] -46545.8607 33186.3418
[308,] 5273.5049 -46545.8607
[309,] 33186.3418 5273.5049
[310,] 33186.3418 33186.3418
[311,] -21898.2339 33186.3418
[312,] 33186.3418 -21898.2339
[313,] -19846.2811 33186.3418
[314,] 49461.5482 -19846.2811
[315,] 90303.0119 49461.5482
[316,] 82874.4756 90303.0119
[317,] -20210.8429 82874.4756
[318,] 61966.8960 -20210.8429
[319,] 83290.4756 61966.8960
[320,] 83290.4756 83290.4756
[321,] 83290.4756 83290.4756
[322,] 83290.4756 83290.4756
[323,] 33186.3418 83290.4756
[324,] 83290.4756 33186.3418
[325,] 76149.3597 83290.4756
[326,] 33186.3418 76149.3597
[327,] 33186.3418 33186.3418
[328,] 83290.4756 33186.3418
[329,] 83290.4756 83290.4756
[330,] 83290.4756 83290.4756
[331,] 83290.4756 83290.4756
[332,] 83290.4756 83290.4756
[333,] 64587.0119 83290.4756
[334,] 79035.3597 64587.0119
[335,] 66229.0119 79035.3597
[336,] -21212.8429 66229.0119
[337,] 83290.4756 -21212.8429
[338,] 77527.3597 83290.4756
[339,] 42977.3164 77527.3597
[340,] 12333.0845 42977.3164
[341,] 83290.4756 12333.0845
[342,] 83290.4756 83290.4756
[343,] 83290.4756 83290.4756
[344,] 20745.9686 83290.4756
[345,] 83290.4756 20745.9686
[346,] 58798.1278 83290.4756
[347,] 83290.4756 58798.1278
[348,] 83290.4756 83290.4756
[349,] 44648.3164 83290.4756
[350,] 73811.1278 44648.3164
[351,] 52117.6641 73811.1278
[352,] 50978.3164 52117.6641
[353,] 77857.3597 50978.3164
[354,] 83290.4756 77857.3597
[355,] 83290.4756 83290.4756
[356,] 26341.7368 83290.4756
[357,] 83290.4756 26341.7368
[358,] 83290.4756 83290.4756
[359,] 83290.4756 83290.4756
[360,] 66256.0119 83290.4756
[361,] 83290.4756 66256.0119
[362,] 83290.4756 83290.4756
[363,] -93224.9040 83290.4756
[364,] 45928.4323 -93224.9040
[365,] 83290.4756 45928.4323
[366,] 83290.4756 83290.4756
[367,] 83290.4756 83290.4756
[368,] 83290.4756 83290.4756
[369,] 53621.6641 83290.4756
[370,] 83290.4756 53621.6641
[371,] 61429.6641 83290.4756
[372,] 74707.6641 61429.6641
[373,] 62545.1278 74707.6641
[374,] 83290.4756 62545.1278
[375,] 83290.4756 83290.4756
[376,] 16386.2004 83290.4756
[377,] 60408.6641 16386.2004
[378,] 45143.3164 60408.6641
[379,] 83290.4756 45143.3164
[380,] 83290.4756 83290.4756
[381,] 8465.3890 83290.4756
[382,] -11375.1473 8465.3890
[383,] 83290.4756 -11375.1473
[384,] 34083.0845 83290.4756
[385,] -117594.8862 34083.0845
[386,] 82192.1278 -117594.8862
[387,] 51369.7801 82192.1278
[388,] 69331.1278 51369.7801
[389,] 41735.5482 69331.1278
[390,] 47606.3164 41735.5482
[391,] 34467.7801 47606.3164
[392,] -2918.4518 34467.7801
[393,] 109986.7801 -2918.4518
[394,] -36651.7015 109986.7801
[395,] 17392.5482 -36651.7015
[396,] -149384.6722 17392.5482
[397,] -862585.3637 -149384.6722
[398,] 34329.6641 -862585.3637
[399,] -226360.4110 34329.6641
[400,] 59373.1278 -226360.4110
[401,] 67310.5482 59373.1278
[402,] 39250.7801 67310.5482
[403,] -92729.0926 39250.7801
[404,] 14966.6641 -92729.0926
[405,] 48023.2004 14966.6641
[406,] 38062.2004 48023.2004
[407,] 21308.0845 38062.2004
[408,] 27638.4323 21308.0845
[409,] -13017.8722 27638.4323
[410,] -165665.1320 -13017.8722
[411,] -84756.0926 -165665.1320
[412,] -76387.3970 -84756.0926
[413,] 29629.6323 -76387.3970
[414,] -113579.7996 29629.6323
[415,] -36290.7562 -113579.7996
[416,] 77017.8094 -36290.7562
[417,] -76979.9767 77017.8094
[418,] -59093.3065 -76979.9767
[419,] 10188.3597 -59093.3065
[420,] 17115.3597 10188.3597
[421,] -306663.1027 17115.3597
[422,] -655350.5409 -306663.1027
[423,] -359494.8429 -655350.5409
[424,] -44638.8429 -359494.8429
[425,] -416407.1639 -44638.8429
[426,] -278860.9588 -416407.1639
[427,] -966014.8096 -278860.9588
[428,] -114972.2339 -966014.8096
[429,] -226086.4938 -114972.2339
[430,] -533414.5231 -226086.4938
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2146005.8721 2559625.7921
2 2073314.5816 2146005.8721
3 -1833510.1393 2073314.5816
4 -133727.1739 -1833510.1393
5 631218.8668 -133727.1739
6 45603.9879 631218.8668
7 506071.3904 45603.9879
8 -270524.0886 506071.3904
9 611450.3878 -270524.0886
10 397635.4196 611450.3878
11 232695.3330 397635.4196
12 -552373.2158 232695.3330
13 -814949.5497 -552373.2158
14 -546748.1166 -814949.5497
15 570626.1725 -546748.1166
16 -28097.3485 570626.1725
17 -554963.3777 -28097.3485
18 -287336.7841 -554963.3777
19 -968495.6350 -287336.7841
20 263815.8794 -968495.6350
21 144146.5750 263815.8794
22 31189.5356 144146.5750
23 -1109804.1420 31189.5356
24 -287950.1460 -1109804.1420
25 114890.9839 -287950.1460
26 -630335.0160 114890.9839
27 -333838.9868 -630335.0160
28 409654.5342 -333838.9868
29 -384864.6377 409654.5342
30 295744.3278 -384864.6377
31 -6272.7728 295744.3278
32 288690.7049 -6272.7728
33 136548.4730 288690.7049
34 360943.0845 136548.4730
35 -753672.9841 360943.0845
36 92712.8973 -753672.9841
37 74215.0680 92712.8973
38 -85424.2798 74215.0680
39 -583164.7103 -85424.2798
40 -795895.0567 -583164.7103
41 -125041.0161 -795895.0567
42 -109757.8415 -125041.0161
43 -227260.7256 -109757.8415
44 -10876.0454 -227260.7256
45 -11466.5563 -10876.0454
46 -120249.3384 -11466.5563
47 -486531.9307 -120249.3384
48 -60347.3817 -486531.9307
49 -212384.2772 -60347.3817
50 -135204.0912 -212384.2772
51 -343612.4237 -135204.0912
52 -204412.9027 -343612.4237
53 -23748.9613 -204412.9027
54 -534415.8734 -23748.9613
55 -154377.4505 -534415.8734
56 395256.7801 -154377.4505
57 13755.0272 395256.7801
58 -108168.4365 13755.0272
59 -256859.5116 -108168.4365
60 27900.3571 -256859.5116
61 -222015.2339 27900.3571
62 25690.4437 -222015.2339
63 -246149.5843 25690.4437
64 -107119.7842 -246149.5843
65 -944371.6210 -107119.7842
66 -34820.2632 -944371.6210
67 -72626.7588 -34820.2632
68 -168353.9868 -72626.7588
69 -219660.5116 -168353.9868
70 -289031.3384 -219660.5116
71 -15562.9295 -289031.3384
72 -70530.5817 -15562.9295
73 55280.4323 -70530.5817
74 17474.5342 55280.4323
75 -120395.9767 17474.5342
76 12358.1138 -120395.9767
77 -42222.7996 12358.1138
78 -141393.5524 -42222.7996
79 2221.8348 -141393.5524
80 -8334.8862 2221.8348
81 -137279.3091 -8334.8862
82 -19922.3499 -137279.3091
83 5844.7368 -19922.3499
84 81479.7368 5844.7368
85 -214773.5409 81479.7368
86 -200058.1167 -214773.5409
87 66717.9686 -200058.1167
88 -115293.7409 66717.9686
89 94438.8527 -115293.7409
90 7781.3024 94438.8527
91 -74209.2951 7781.3024
92 -69512.4071 -74209.2951
93 -40985.0200 -69512.4071
94 -56754.5091 -40985.0200
95 -33446.1359 -56754.5091
96 158654.4323 -33446.1359
97 -45436.2085 158654.4323
98 90471.2731 -45436.2085
99 -2669.7015 90471.2731
100 -106822.2518 -2669.7015
101 -77955.7269 -106822.2518
102 -167317.9002 -77955.7269
103 47878.5482 -167317.9002
104 16729.2298 47878.5482
105 -86058.8862 16729.2298
106 -361433.2225 -86058.8862
107 160430.4870 -361433.2225
108 -110559.9728 160430.4870
109 -117557.9767 -110559.9728
110 83002.9253 -117557.9767
111 -264384.4836 83002.9253
112 -127512.6250 -264384.4836
113 -48101.1652 -127512.6250
114 70276.6463 -48101.1652
115 -37832.5817 70276.6463
116 66548.5482 -37832.5817
117 138319.5049 66548.5482
118 -114872.9333 138319.5049
119 -5926.9588 -114872.9333
120 -249101.1893 -5926.9588
121 51033.0845 -249101.1893
122 -353694.5078 51033.0845
123 91610.4437 -353694.5078
124 -1608.6289 91610.4437
125 -163052.3244 -1608.6289
126 83290.4756 -163052.3244
127 39812.3164 83290.4756
128 -12085.1906 39812.3164
129 83290.4756 -12085.1906
130 52421.6641 83290.4756
131 30780.5775 52421.6641
132 83290.4756 30780.5775
133 -4505.8429 83290.4756
134 62933.0119 -4505.8429
135 65867.0119 62933.0119
136 44182.3164 65867.0119
137 78534.3597 44182.3164
138 65771.0119 78534.3597
139 83290.4756 65771.0119
140 59311.5482 83290.4756
141 83290.4756 59311.5482
142 19584.2004 83290.4756
143 66692.0119 19584.2004
144 33747.7368 66692.0119
145 42420.4323 33747.7368
146 83290.4756 42420.4323
147 53055.6641 83290.4756
148 84633.2438 53055.6641
149 83787.4756 84633.2438
150 29718.6641 83787.4756
151 -122233.8568 29718.6641
152 61974.8960 -122233.8568
153 35833.6641 61974.8960
154 74629.2438 35833.6641
155 83290.4756 74629.2438
156 83290.4756 83290.4756
157 74620.2438 83290.4756
158 83019.8960 74620.2438
159 83290.4756 83019.8960
160 83290.4756 83290.4756
161 83290.4756 83290.4756
162 83290.4756 83290.4756
163 44336.6641 83290.4756
164 -32217.3065 44336.6641
165 78529.3597 -32217.3065
166 83290.4756 78529.3597
167 73886.3890 83290.4756
168 -124990.9002 73886.3890
169 83290.4756 -124990.9002
170 83290.4756 83290.4756
171 67533.1278 83290.4756
172 83290.4756 67533.1278
173 39138.4323 83290.4756
174 83290.4756 39138.4323
175 49352.5482 83290.4756
176 -400554.6237 49352.5482
177 33982.3418 -400554.6237
178 83290.4756 33982.3418
179 83290.4756 83290.4756
180 46246.2731 83290.4756
181 980.0667 46246.2731
182 -24378.4225 980.0667
183 -52990.5130 -24378.4225
184 11296.7622 -52990.5130
185 4122.0667 11296.7622
186 -28137.3244 4122.0667
187 55663.7801 -28137.3244
188 38381.8960 55663.7801
189 83290.4756 38381.8960
190 56813.6641 83290.4756
191 83290.4756 56813.6641
192 53434.6641 83290.4756
193 83290.4756 53434.6641
194 70809.1278 83290.4756
195 33186.3418 70809.1278
196 -251968.9753 33186.3418
197 33186.3418 -251968.9753
198 33186.3418 33186.3418
199 -93605.8607 33186.3418
200 33186.3418 -93605.8607
201 33186.3418 33186.3418
202 33186.3418 33186.3418
203 -32389.0493 33186.3418
204 -7872.0493 -32389.0493
205 12645.7622 -7872.0493
206 27854.1100 12645.7622
207 33186.3418 27854.1100
208 -49421.6722 33186.3418
209 11964.9940 -49421.6722
210 83290.4756 11964.9940
211 83290.4756 83290.4756
212 51075.2004 83290.4756
213 14264.3164 51075.2004
214 51079.8960 14264.3164
215 74764.2438 51079.8960
216 58019.7801 74764.2438
217 55804.6641 58019.7801
218 74764.2438 55804.6641
219 -29125.7702 74764.2438
220 51777.1278 -29125.7702
221 83290.4756 51777.1278
222 79001.3597 83290.4756
223 -5190.8429 79001.3597
224 79223.3597 -5190.8429
225 -10956.0747 79223.3597
226 44093.4323 -10956.0747
227 83290.4756 44093.4323
228 71723.2438 83290.4756
229 83290.4756 71723.2438
230 36768.6641 83290.4756
231 50918.4323 36768.6641
232 83290.4756 50918.4323
233 15719.9686 83290.4756
234 83290.4756 15719.9686
235 83397.4756 83290.4756
236 83290.4756 83397.4756
237 58124.7801 83290.4756
238 83290.4756 58124.7801
239 83290.4756 83290.4756
240 83290.4756 83290.4756
241 68008.1278 83290.4756
242 83290.4756 68008.1278
243 53705.6641 83290.4756
244 83769.2438 53705.6641
245 33186.3418 83769.2438
246 33186.3418 33186.3418
247 12952.7368 33186.3418
248 33186.3418 12952.7368
249 33186.3418 33186.3418
250 26244.4323 33186.3418
251 33186.3418 26244.4323
252 70969.3597 33186.3418
253 -89711.9295 70969.3597
254 -221499.1320 -89711.9295
255 83290.4756 -221499.1320
256 47628.6641 83290.4756
257 -44464.5384 47628.6641
258 33186.3418 -44464.5384
259 83290.4756 33186.3418
260 33186.3418 83290.4756
261 33186.3418 33186.3418
262 56418.9686 33186.3418
263 57234.3164 56418.9686
264 33186.3418 57234.3164
265 83290.4756 33186.3418
266 33186.3418 83290.4756
267 110722.4323 33186.3418
268 33186.3418 110722.4323
269 -11510.5384 33186.3418
270 33186.3418 -11510.5384
271 33186.3418 33186.3418
272 33186.3418 33186.3418
273 -2887.7269 33186.3418
274 83290.4756 -2887.7269
275 36769.0845 83290.4756
276 -9743.1906 36769.0845
277 83290.4756 -9743.1906
278 -113983.0200 83290.4756
279 44475.5482 -113983.0200
280 33186.3418 44475.5482
281 83290.4756 33186.3418
282 -29439.5384 83290.4756
283 83290.4756 -29439.5384
284 60865.4323 83290.4756
285 83290.4756 60865.4323
286 15700.7622 83290.4756
287 73219.1278 15700.7622
288 33186.3418 73219.1278
289 5243.9686 33186.3418
290 33186.3418 5243.9686
291 41929.0845 33186.3418
292 52870.2731 41929.0845
293 83290.4756 52870.2731
294 27688.3164 83290.4756
295 20406.2985 27688.3164
296 83290.4756 20406.2985
297 33186.3418 83290.4756
298 83290.4756 33186.3418
299 44922.4323 83290.4756
300 13917.8781 44922.4323
301 64811.3597 13917.8781
302 79027.3597 64811.3597
303 33186.3418 79027.3597
304 52651.8527 33186.3418
305 32147.0845 52651.8527
306 33186.3418 32147.0845
307 -46545.8607 33186.3418
308 5273.5049 -46545.8607
309 33186.3418 5273.5049
310 33186.3418 33186.3418
311 -21898.2339 33186.3418
312 33186.3418 -21898.2339
313 -19846.2811 33186.3418
314 49461.5482 -19846.2811
315 90303.0119 49461.5482
316 82874.4756 90303.0119
317 -20210.8429 82874.4756
318 61966.8960 -20210.8429
319 83290.4756 61966.8960
320 83290.4756 83290.4756
321 83290.4756 83290.4756
322 83290.4756 83290.4756
323 33186.3418 83290.4756
324 83290.4756 33186.3418
325 76149.3597 83290.4756
326 33186.3418 76149.3597
327 33186.3418 33186.3418
328 83290.4756 33186.3418
329 83290.4756 83290.4756
330 83290.4756 83290.4756
331 83290.4756 83290.4756
332 83290.4756 83290.4756
333 64587.0119 83290.4756
334 79035.3597 64587.0119
335 66229.0119 79035.3597
336 -21212.8429 66229.0119
337 83290.4756 -21212.8429
338 77527.3597 83290.4756
339 42977.3164 77527.3597
340 12333.0845 42977.3164
341 83290.4756 12333.0845
342 83290.4756 83290.4756
343 83290.4756 83290.4756
344 20745.9686 83290.4756
345 83290.4756 20745.9686
346 58798.1278 83290.4756
347 83290.4756 58798.1278
348 83290.4756 83290.4756
349 44648.3164 83290.4756
350 73811.1278 44648.3164
351 52117.6641 73811.1278
352 50978.3164 52117.6641
353 77857.3597 50978.3164
354 83290.4756 77857.3597
355 83290.4756 83290.4756
356 26341.7368 83290.4756
357 83290.4756 26341.7368
358 83290.4756 83290.4756
359 83290.4756 83290.4756
360 66256.0119 83290.4756
361 83290.4756 66256.0119
362 83290.4756 83290.4756
363 -93224.9040 83290.4756
364 45928.4323 -93224.9040
365 83290.4756 45928.4323
366 83290.4756 83290.4756
367 83290.4756 83290.4756
368 83290.4756 83290.4756
369 53621.6641 83290.4756
370 83290.4756 53621.6641
371 61429.6641 83290.4756
372 74707.6641 61429.6641
373 62545.1278 74707.6641
374 83290.4756 62545.1278
375 83290.4756 83290.4756
376 16386.2004 83290.4756
377 60408.6641 16386.2004
378 45143.3164 60408.6641
379 83290.4756 45143.3164
380 83290.4756 83290.4756
381 8465.3890 83290.4756
382 -11375.1473 8465.3890
383 83290.4756 -11375.1473
384 34083.0845 83290.4756
385 -117594.8862 34083.0845
386 82192.1278 -117594.8862
387 51369.7801 82192.1278
388 69331.1278 51369.7801
389 41735.5482 69331.1278
390 47606.3164 41735.5482
391 34467.7801 47606.3164
392 -2918.4518 34467.7801
393 109986.7801 -2918.4518
394 -36651.7015 109986.7801
395 17392.5482 -36651.7015
396 -149384.6722 17392.5482
397 -862585.3637 -149384.6722
398 34329.6641 -862585.3637
399 -226360.4110 34329.6641
400 59373.1278 -226360.4110
401 67310.5482 59373.1278
402 39250.7801 67310.5482
403 -92729.0926 39250.7801
404 14966.6641 -92729.0926
405 48023.2004 14966.6641
406 38062.2004 48023.2004
407 21308.0845 38062.2004
408 27638.4323 21308.0845
409 -13017.8722 27638.4323
410 -165665.1320 -13017.8722
411 -84756.0926 -165665.1320
412 -76387.3970 -84756.0926
413 29629.6323 -76387.3970
414 -113579.7996 29629.6323
415 -36290.7562 -113579.7996
416 77017.8094 -36290.7562
417 -76979.9767 77017.8094
418 -59093.3065 -76979.9767
419 10188.3597 -59093.3065
420 17115.3597 10188.3597
421 -306663.1027 17115.3597
422 -655350.5409 -306663.1027
423 -359494.8429 -655350.5409
424 -44638.8429 -359494.8429
425 -416407.1639 -44638.8429
426 -278860.9588 -416407.1639
427 -966014.8096 -278860.9588
428 -114972.2339 -966014.8096
429 -226086.4938 -114972.2339
430 -533414.5231 -226086.4938
> 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/7d6gw1291232858.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/8d6gw1291232858.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/9d6gw1291232858.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/106yfh1291232858.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/11rgd51291232858.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/12chub1291232858.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/138qa11291232858.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/14c9871291232858.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/15fapv1291232858.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/161a5j1291232858.tab")
+ }
>
> try(system("convert tmp/1hxin1291232858.ps tmp/1hxin1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/2hxin1291232858.ps tmp/2hxin1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/3r6hq1291232858.ps tmp/3r6hq1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/4r6hq1291232858.ps tmp/4r6hq1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/5r6hq1291232858.ps tmp/5r6hq1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/62fzt1291232858.ps tmp/62fzt1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/7d6gw1291232858.ps tmp/7d6gw1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/8d6gw1291232858.ps tmp/8d6gw1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/9d6gw1291232858.ps tmp/9d6gw1291232858.png",intern=TRUE))
character(0)
> try(system("convert tmp/106yfh1291232858.ps tmp/106yfh1291232858.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.085 2.220 23.406