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
+ ,162556
+ ,807
+ ,213118
+ ,0
+ ,29790
+ ,444
+ ,81767
+ ,0
+ ,87550
+ ,412
+ ,153198
+ ,1
+ ,84738
+ ,428
+ ,-26007
+ ,0
+ ,54660
+ ,315
+ ,126942
+ ,0
+ ,42634
+ ,168
+ ,157214
+ ,1
+ ,40949
+ ,263
+ ,129352
+ ,0
+ ,45187
+ ,267
+ ,234817
+ ,0
+ ,37704
+ ,228
+ ,60448
+ ,0
+ ,16275
+ ,129
+ ,47818
+ ,1
+ ,25830
+ ,104
+ ,245546
+ ,1
+ ,12679
+ ,122
+ ,48020
+ ,0
+ ,18014
+ ,393
+ ,-1710
+ ,1
+ ,43556
+ ,190
+ ,32648
+ ,0
+ ,24811
+ ,280
+ ,95350
+ ,1
+ ,6575
+ ,63
+ ,151352
+ ,1
+ ,7123
+ ,102
+ ,288170
+ ,0
+ ,21950
+ ,265
+ ,114337
+ ,0
+ ,37597
+ ,234
+ ,37884
+ ,1
+ ,17821
+ ,277
+ ,122844
+ ,0
+ ,12988
+ ,73
+ ,82340
+ ,0
+ ,22330
+ ,67
+ ,79801
+ ,1
+ ,13326
+ ,103
+ ,165548
+ ,1
+ ,16189
+ ,290
+ ,116384
+ ,1
+ ,7146
+ ,83
+ ,134028
+ ,1
+ ,15824
+ ,56
+ ,63838
+ ,0
+ ,27664
+ ,236
+ ,74996
+ ,1
+ ,11920
+ ,73
+ ,31080
+ ,1
+ ,8568
+ ,34
+ ,32168
+ ,1
+ ,14416
+ ,139
+ ,49857
+ ,0
+ ,3369
+ ,26
+ ,87161
+ ,0
+ ,11819
+ ,70
+ ,106113
+ ,0
+ ,6984
+ ,40
+ ,80570
+ ,0
+ ,4519
+ ,42
+ ,102129
+ ,1
+ ,2220
+ ,12
+ ,301670
+ ,1
+ ,18562
+ ,211
+ ,102313
+ ,1
+ ,10327
+ ,74
+ ,88577
+ ,0
+ ,5336
+ ,80
+ ,112477
+ ,0
+ ,2365
+ ,83
+ ,191778
+ ,1
+ ,4069
+ ,131
+ ,79804
+ ,1
+ ,8636
+ ,203
+ ,128294
+ ,1
+ ,13718
+ ,56
+ ,96448
+ ,1
+ ,4525
+ ,89
+ ,93811
+ ,1
+ ,6869
+ ,88
+ ,117520
+ ,1
+ ,4628
+ ,39
+ ,69159
+ ,0
+ ,3689
+ ,25
+ ,101792
+ ,0
+ ,4891
+ ,49
+ ,210568
+ ,0
+ ,7489
+ ,149
+ ,136996
+ ,1
+ ,4901
+ ,58
+ ,121920
+ ,1
+ ,2284
+ ,41
+ ,76403
+ ,0
+ ,3160
+ ,90
+ ,108094
+ ,0
+ ,4150
+ ,136
+ ,134759
+ ,0
+ ,7285
+ ,97
+ ,188873
+ ,0
+ ,1134
+ ,63
+ ,146216
+ ,0
+ ,4658
+ ,114
+ ,156608
+ ,1
+ ,2384
+ ,77
+ ,61348
+ ,1
+ ,3748
+ ,6
+ ,50350
+ ,1
+ ,5371
+ ,47
+ ,87720
+ ,1
+ ,1285
+ ,51
+ ,99489
+ ,0
+ ,9327
+ ,85
+ ,87419
+ ,0
+ ,5565
+ ,43
+ ,94355
+ ,1
+ ,1528
+ ,32
+ ,60326
+ ,0
+ ,3122
+ ,25
+ ,94670
+ ,0
+ ,7561
+ ,77
+ ,82425
+ ,1
+ ,2675
+ ,54
+ ,59017
+ ,1
+ ,13253
+ ,251
+ ,90829
+ ,1
+ ,880
+ ,15
+ ,80791
+ ,0
+ ,2053
+ ,44
+ ,100423
+ ,1
+ ,1424
+ ,73
+ ,131116
+ ,0
+ ,4036
+ ,85
+ ,100269
+ ,0
+ ,3045
+ ,49
+ ,27330
+ ,1
+ ,5119
+ ,38
+ ,39039
+ ,1
+ ,1431
+ ,35
+ ,106885
+ ,1
+ ,554
+ ,9
+ ,79285
+ ,1
+ ,1975
+ ,34
+ ,118881
+ ,0
+ ,1765
+ ,20
+ ,77623
+ ,1
+ ,1012
+ ,29
+ ,114768
+ ,1
+ ,810
+ ,11
+ ,74015
+ ,1
+ ,1280
+ ,52
+ ,69465
+ ,0
+ ,666
+ ,13
+ ,117869
+ ,1
+ ,1380
+ ,29
+ ,60982
+ ,0
+ ,4677
+ ,66
+ ,90131
+ ,1
+ ,876
+ ,33
+ ,138971
+ ,1
+ ,814
+ ,15
+ ,39625
+ ,1
+ ,514
+ ,15
+ ,102725
+ ,0
+ ,5692
+ ,68
+ ,64239
+ ,1
+ ,3642
+ ,100
+ ,90262
+ ,1
+ ,540
+ ,13
+ ,103960
+ ,1
+ ,2099
+ ,45
+ ,106611
+ ,1
+ ,567
+ ,14
+ ,103345
+ ,1
+ ,2001
+ ,36
+ ,95551
+ ,0
+ ,2949
+ ,40
+ ,82903
+ ,1
+ ,2253
+ ,68
+ ,63593
+ ,0
+ ,6533
+ ,29
+ ,126910
+ ,1
+ ,1889
+ ,43
+ ,37527
+ ,0
+ ,3055
+ ,30
+ ,60247
+ ,1
+ ,272
+ ,9
+ ,112995
+ ,0
+ ,1414
+ ,22
+ ,70184
+ ,1
+ ,2564
+ ,19
+ ,130140
+ ,0
+ ,1383
+ ,9
+ ,73221
+ ,0
+ ,1261
+ ,31
+ ,76114
+ ,1
+ ,975
+ ,19
+ ,90534
+ ,1
+ ,3366
+ ,55
+ ,108479
+ ,1
+ ,576
+ ,8
+ ,113761
+ ,1
+ ,1686
+ ,28
+ ,68696
+ ,1
+ ,746
+ ,29
+ ,71561
+ ,0
+ ,3192
+ ,48
+ ,59831
+ ,0
+ ,2045
+ ,16
+ ,97890
+ ,1
+ ,5702
+ ,47
+ ,101481
+ ,0
+ ,1932
+ ,20
+ ,72954
+ ,1
+ ,936
+ ,22
+ ,67939
+ ,0
+ ,3437
+ ,33
+ ,48022
+ ,1
+ ,5131
+ ,44
+ ,86111
+ ,0
+ ,2397
+ ,13
+ ,74020
+ ,0
+ ,1389
+ ,6
+ ,57530
+ ,1
+ ,1503
+ ,35
+ ,56364
+ ,1
+ ,402
+ ,8
+ ,84990
+ ,1
+ ,2239
+ ,17
+ ,88590
+ ,0
+ ,2234
+ ,11
+ ,77200
+ ,1
+ ,837
+ ,21
+ ,61262
+ ,1
+ ,10579
+ ,92
+ ,110309
+ ,1
+ ,875
+ ,12
+ ,67000
+ ,1
+ ,1585
+ ,112
+ ,93099
+ ,0
+ ,1659
+ ,25
+ ,107577
+ ,0
+ ,2647
+ ,17
+ ,62920
+ ,0
+ ,3294
+ ,23
+ ,75832
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,94
+ ,10
+ ,60793
+ ,1
+ ,422
+ ,23
+ ,57935
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,34
+ ,7
+ ,60630
+ ,1
+ ,1558
+ ,25
+ ,55637
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,43
+ ,20
+ ,60887
+ ,1
+ ,645
+ ,4
+ ,60720
+ ,1
+ ,316
+ ,4
+ ,60505
+ ,1
+ ,115
+ ,10
+ ,60945
+ ,1
+ ,5
+ ,1
+ ,60720
+ ,1
+ ,897
+ ,4
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,389
+ ,8
+ ,58990
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1002
+ ,11
+ ,56750
+ ,1
+ ,36
+ ,4
+ ,60894
+ ,1
+ ,460
+ ,15
+ ,63346
+ ,1
+ ,309
+ ,9
+ ,56535
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,9
+ ,7
+ ,60835
+ ,1
+ ,271
+ ,2
+ ,60720
+ ,1
+ ,14
+ ,0
+ ,61016
+ ,1
+ ,520
+ ,7
+ ,58650
+ ,1
+ ,1766
+ ,46
+ ,60438
+ ,1
+ ,0
+ ,5
+ ,60720
+ ,1
+ ,458
+ ,7
+ ,58625
+ ,1
+ ,20
+ ,2
+ ,60938
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,98
+ ,2
+ ,61490
+ ,1
+ ,405
+ ,5
+ ,60845
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,483
+ ,7
+ ,60830
+ ,1
+ ,454
+ ,24
+ ,63261
+ ,1
+ ,47
+ ,1
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,757
+ ,18
+ ,45689
+ ,1
+ ,4655
+ ,55
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,36
+ ,3
+ ,61564
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,203
+ ,9
+ ,61938
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,126
+ ,8
+ ,60951
+ ,1
+ ,400
+ ,113
+ ,60720
+ ,0
+ ,71
+ ,0
+ ,60745
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,972
+ ,19
+ ,71642
+ ,0
+ ,531
+ ,11
+ ,71641
+ ,1
+ ,2461
+ ,25
+ ,55792
+ ,0
+ ,378
+ ,16
+ ,71873
+ ,0
+ ,23
+ ,5
+ ,62555
+ ,0
+ ,638
+ ,11
+ ,60370
+ ,0
+ ,2300
+ ,23
+ ,64873
+ ,1
+ ,149
+ ,6
+ ,62041
+ ,1
+ ,226
+ ,5
+ ,65745
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,275
+ ,7
+ ,59500
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,141
+ ,7
+ ,61630
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,28
+ ,3
+ ,60890
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,4980
+ ,89
+ ,113521
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,472
+ ,19
+ ,80045
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,203
+ ,12
+ ,50804
+ ,0
+ ,496
+ ,12
+ ,87390
+ ,0
+ ,10
+ ,5
+ ,61656
+ ,0
+ ,63
+ ,2
+ ,65688
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1136
+ ,26
+ ,48522
+ ,0
+ ,265
+ ,3
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,267
+ ,11
+ ,57640
+ ,1
+ ,474
+ ,10
+ ,61977
+ ,1
+ ,534
+ ,5
+ ,62620
+ ,1
+ ,0
+ ,2
+ ,60720
+ ,1
+ ,15
+ ,6
+ ,60831
+ ,1
+ ,397
+ ,7
+ ,60646
+ ,1
+ ,0
+ ,2
+ ,60720
+ ,1
+ ,1866
+ ,28
+ ,56225
+ ,1
+ ,288
+ ,3
+ ,60510
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,3
+ ,1
+ ,60698
+ ,1
+ ,468
+ ,20
+ ,60720
+ ,1
+ ,20
+ ,1
+ ,60805
+ ,1
+ ,278
+ ,22
+ ,61404
+ ,1
+ ,61
+ ,9
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,192
+ ,2
+ ,65276
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,317
+ ,7
+ ,63915
+ ,1
+ ,738
+ ,9
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,368
+ ,13
+ ,61686
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,2
+ ,0
+ ,60743
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,53
+ ,6
+ ,60349
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,94
+ ,3
+ ,61360
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,24
+ ,7
+ ,59818
+ ,1
+ ,2332
+ ,2
+ ,72680
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,131
+ ,15
+ ,61808
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,206
+ ,9
+ ,53110
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,167
+ ,1
+ ,64245
+ ,1
+ ,622
+ ,38
+ ,73007
+ ,1
+ ,2328
+ ,57
+ ,82732
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,365
+ ,7
+ ,54820
+ ,1
+ ,364
+ ,26
+ ,47705
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,226
+ ,13
+ ,72835
+ ,1
+ ,307
+ ,10
+ ,58856
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,188
+ ,9
+ ,77655
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,138
+ ,26
+ ,69817
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,125
+ ,19
+ ,60798
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,282
+ ,12
+ ,62452
+ ,1
+ ,335
+ ,23
+ ,64175
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1324
+ ,29
+ ,67440
+ ,1
+ ,176
+ ,8
+ ,68136
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,249
+ ,26
+ ,56726
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,333
+ ,9
+ ,70811
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,601
+ ,5
+ ,60720
+ ,1
+ ,30
+ ,3
+ ,62045
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,249
+ ,13
+ ,54323
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,165
+ ,12
+ ,62841
+ ,1
+ ,453
+ ,19
+ ,81125
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,53
+ ,10
+ ,59506
+ ,0
+ ,382
+ ,9
+ ,59365
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,9
+ ,60720
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,1
+ ,290
+ ,1
+ ,58790
+ ,1
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,366
+ ,14
+ ,61808
+ ,1
+ ,2
+ ,12
+ ,60735
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,209
+ ,19
+ ,64016
+ ,1
+ ,384
+ ,17
+ ,54683
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,365
+ ,32
+ ,87192
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,49
+ ,14
+ ,64107
+ ,1
+ ,3
+ ,8
+ ,60761
+ ,1
+ ,133
+ ,4
+ ,65990
+ ,1
+ ,32
+ ,0
+ ,59988
+ ,1
+ ,368
+ ,20
+ ,61167
+ ,1
+ ,1
+ ,5
+ ,60719
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,22
+ ,1
+ ,60722
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,96
+ ,4
+ ,60379
+ ,1
+ ,1
+ ,1
+ ,60727
+ ,1
+ ,314
+ ,4
+ ,60720
+ ,1
+ ,844
+ ,20
+ ,60925
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,26
+ ,1
+ ,60896
+ ,1
+ ,125
+ ,10
+ ,59734
+ ,1
+ ,304
+ ,12
+ ,62969
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,621
+ ,13
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,119
+ ,3
+ ,59118
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1595
+ ,10
+ ,60720
+ ,1
+ ,312
+ ,3
+ ,58598
+ ,1
+ ,60
+ ,7
+ ,61124
+ ,1
+ ,587
+ ,10
+ ,59595
+ ,1
+ ,135
+ ,1
+ ,62065
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,514
+ ,15
+ ,78780
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1
+ ,4
+ ,60722
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1763
+ ,28
+ ,61600
+ ,1
+ ,180
+ ,9
+ ,59635
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,218
+ ,7
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,448
+ ,7
+ ,59781
+ ,1
+ ,227
+ ,7
+ ,76644
+ ,1
+ ,174
+ ,3
+ ,64820
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,121
+ ,11
+ ,56178
+ ,1
+ ,607
+ ,7
+ ,60436
+ ,1
+ ,2212
+ ,10
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,530
+ ,18
+ ,73433
+ ,1
+ ,571
+ ,14
+ ,41477
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,78
+ ,12
+ ,62700
+ ,1
+ ,2489
+ ,29
+ ,67804
+ ,1
+ ,131
+ ,3
+ ,59661
+ ,1
+ ,923
+ ,6
+ ,58620
+ ,1
+ ,72
+ ,3
+ ,60398
+ ,1
+ ,572
+ ,8
+ ,58580
+ ,1
+ ,397
+ ,10
+ ,62710
+ ,1
+ ,450
+ ,6
+ ,59325
+ ,1
+ ,622
+ ,8
+ ,60950
+ ,1
+ ,694
+ ,6
+ ,68060
+ ,0
+ ,3425
+ ,9
+ ,83620
+ ,1
+ ,562
+ ,8
+ ,58456
+ ,1
+ ,4917
+ ,26
+ ,52811
+ ,0
+ ,1442
+ ,239
+ ,121173
+ ,1
+ ,529
+ ,7
+ ,63870
+ ,0
+ ,2126
+ ,41
+ ,21001
+ ,1
+ ,1061
+ ,3
+ ,70415
+ ,1
+ ,776
+ ,8
+ ,64230
+ ,1
+ ,611
+ ,6
+ ,59190
+ ,0
+ ,1526
+ ,21
+ ,69351
+ ,1
+ ,592
+ ,7
+ ,64270
+ ,1
+ ,1182
+ ,11
+ ,70694
+ ,1
+ ,621
+ ,11
+ ,68005
+ ,1
+ ,989
+ ,12
+ ,58930
+ ,1
+ ,438
+ ,9
+ ,58320
+ ,1
+ ,726
+ ,3
+ ,69980
+ ,1
+ ,1303
+ ,57
+ ,69863
+ ,0
+ ,7419
+ ,21
+ ,63255
+ ,0
+ ,1164
+ ,15
+ ,57320
+ ,0
+ ,3310
+ ,32
+ ,75230
+ ,1
+ ,1920
+ ,11
+ ,79420
+ ,1
+ ,965
+ ,2
+ ,73490
+ ,1
+ ,3256
+ ,23
+ ,35250
+ ,0
+ ,1135
+ ,20
+ ,62285
+ ,1
+ ,1270
+ ,24
+ ,69206
+ ,1
+ ,661
+ ,1
+ ,65920
+ ,1
+ ,1013
+ ,1
+ ,69770
+ ,1
+ ,2844
+ ,74
+ ,72683
+ ,0
+ ,11528
+ ,68
+ ,-14545
+ ,1
+ ,6526
+ ,20
+ ,55830
+ ,1
+ ,2264
+ ,20
+ ,55174
+ ,0
+ ,5109
+ ,82
+ ,67038
+ ,1
+ ,3999
+ ,21
+ ,51252
+ ,1
+ ,35624
+ ,244
+ ,157278
+ ,1
+ ,9252
+ ,32
+ ,79510
+ ,1
+ ,15236
+ ,86
+ ,77440
+ ,1
+ ,18073
+ ,69
+ ,27284)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Orders'
+ ,'Dividends
')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Group','Costs','Orders','Dividends
'),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
Orders Group Costs Dividends\r
1 807 0 162556 213118
2 444 0 29790 81767
3 412 0 87550 153198
4 428 1 84738 -26007
5 315 0 54660 126942
6 168 0 42634 157214
7 263 1 40949 129352
8 267 0 45187 234817
9 228 0 37704 60448
10 129 0 16275 47818
11 104 1 25830 245546
12 122 1 12679 48020
13 393 0 18014 -1710
14 190 1 43556 32648
15 280 0 24811 95350
16 63 1 6575 151352
17 102 1 7123 288170
18 265 0 21950 114337
19 234 0 37597 37884
20 277 1 17821 122844
21 73 0 12988 82340
22 67 0 22330 79801
23 103 1 13326 165548
24 290 1 16189 116384
25 83 1 7146 134028
26 56 1 15824 63838
27 236 0 27664 74996
28 73 1 11920 31080
29 34 1 8568 32168
30 139 1 14416 49857
31 26 0 3369 87161
32 70 0 11819 106113
33 40 0 6984 80570
34 42 0 4519 102129
35 12 1 2220 301670
36 211 1 18562 102313
37 74 1 10327 88577
38 80 0 5336 112477
39 83 0 2365 191778
40 131 1 4069 79804
41 203 1 8636 128294
42 56 1 13718 96448
43 89 1 4525 93811
44 88 1 6869 117520
45 39 1 4628 69159
46 25 0 3689 101792
47 49 0 4891 210568
48 149 0 7489 136996
49 58 1 4901 121920
50 41 1 2284 76403
51 90 0 3160 108094
52 136 0 4150 134759
53 97 0 7285 188873
54 63 0 1134 146216
55 114 0 4658 156608
56 77 1 2384 61348
57 6 1 3748 50350
58 47 1 5371 87720
59 51 1 1285 99489
60 85 0 9327 87419
61 43 0 5565 94355
62 32 1 1528 60326
63 25 0 3122 94670
64 77 0 7561 82425
65 54 1 2675 59017
66 251 1 13253 90829
67 15 1 880 80791
68 44 0 2053 100423
69 73 1 1424 131116
70 85 0 4036 100269
71 49 0 3045 27330
72 38 1 5119 39039
73 35 1 1431 106885
74 9 1 554 79285
75 34 1 1975 118881
76 20 0 1765 77623
77 29 1 1012 114768
78 11 1 810 74015
79 52 1 1280 69465
80 13 0 666 117869
81 29 1 1380 60982
82 66 0 4677 90131
83 33 1 876 138971
84 15 1 814 39625
85 15 1 514 102725
86 68 0 5692 64239
87 100 1 3642 90262
88 13 1 540 103960
89 45 1 2099 106611
90 14 1 567 103345
91 36 1 2001 95551
92 40 0 2949 82903
93 68 1 2253 63593
94 29 0 6533 126910
95 43 1 1889 37527
96 30 0 3055 60247
97 9 1 272 112995
98 22 0 1414 70184
99 19 1 2564 130140
100 9 0 1383 73221
101 31 0 1261 76114
102 19 1 975 90534
103 55 1 3366 108479
104 8 1 576 113761
105 28 1 1686 68696
106 29 1 746 71561
107 48 0 3192 59831
108 16 0 2045 97890
109 47 1 5702 101481
110 20 0 1932 72954
111 22 1 936 67939
112 33 0 3437 48022
113 44 1 5131 86111
114 13 0 2397 74020
115 6 0 1389 57530
116 35 1 1503 56364
117 8 1 402 84990
118 17 1 2239 88590
119 11 0 2234 77200
120 21 1 837 61262
121 92 1 10579 110309
122 12 1 875 67000
123 112 1 1585 93099
124 25 0 1659 107577
125 17 0 2647 62920
126 23 0 3294 75832
127 0 1 0 60720
128 10 1 94 60793
129 23 1 422 57935
130 0 1 0 60720
131 7 1 34 60630
132 25 1 1558 55637
133 0 1 0 60720
134 20 1 43 60887
135 4 1 645 60720
136 4 1 316 60505
137 10 1 115 60945
138 1 1 5 60720
139 4 1 897 60720
140 0 1 0 60720
141 8 1 389 58990
142 0 1 0 60720
143 11 1 1002 56750
144 4 1 36 60894
145 15 1 460 63346
146 9 1 309 56535
147 0 1 0 60720
148 7 1 9 60835
149 2 1 271 60720
150 0 1 14 61016
151 7 1 520 58650
152 46 1 1766 60438
153 5 1 0 60720
154 7 1 458 58625
155 2 1 20 60938
156 0 1 0 60720
157 0 1 0 60720
158 2 1 98 61490
159 5 1 405 60845
160 0 1 0 60720
161 0 1 0 60720
162 0 1 0 60720
163 0 1 0 60720
164 7 1 483 60830
165 24 1 454 63261
166 1 1 47 60720
167 0 1 0 60720
168 18 1 757 45689
169 55 1 4655 60720
170 0 1 0 60720
171 0 1 0 60720
172 3 1 36 61564
173 0 1 0 60720
174 9 1 203 61938
175 0 1 0 60720
176 8 1 126 60951
177 113 1 400 60720
178 0 0 71 60745
179 0 1 0 60720
180 0 1 0 60720
181 19 1 972 71642
182 11 0 531 71641
183 25 1 2461 55792
184 16 0 378 71873
185 5 0 23 62555
186 11 0 638 60370
187 23 0 2300 64873
188 6 1 149 62041
189 5 1 226 65745
190 0 1 0 60720
191 7 1 275 59500
192 0 1 0 60720
193 7 1 141 61630
194 0 1 0 60720
195 3 1 28 60890
196 0 0 0 60720
197 89 0 4980 113521
198 0 0 0 60720
199 0 0 0 60720
200 19 0 472 80045
201 0 0 0 60720
202 0 0 0 60720
203 0 0 0 60720
204 12 0 203 50804
205 12 0 496 87390
206 5 0 10 61656
207 2 0 63 65688
208 0 0 0 60720
209 26 0 1136 48522
210 3 0 265 60720
211 0 1 0 60720
212 0 1 0 60720
213 11 1 267 57640
214 10 1 474 61977
215 5 1 534 62620
216 2 1 0 60720
217 6 1 15 60831
218 7 1 397 60646
219 2 1 0 60720
220 28 1 1866 56225
221 3 1 288 60510
222 0 1 0 60720
223 1 1 3 60698
224 20 1 468 60720
225 1 1 20 60805
226 22 1 278 61404
227 9 1 61 60720
228 0 1 0 60720
229 2 1 192 65276
230 0 1 0 60720
231 7 1 317 63915
232 9 1 738 60720
233 0 1 0 60720
234 13 1 368 61686
235 0 1 0 60720
236 0 1 2 60743
237 0 1 0 60720
238 6 1 53 60349
239 0 1 0 60720
240 0 1 0 60720
241 0 1 0 60720
242 3 1 94 61360
243 0 1 0 60720
244 7 1 24 59818
245 2 1 2332 72680
246 0 0 0 60720
247 0 0 0 60720
248 15 1 131 61808
249 0 0 0 60720
250 0 0 0 60720
251 9 1 206 53110
252 0 0 0 60720
253 1 1 167 64245
254 38 1 622 73007
255 57 1 2328 82732
256 0 1 0 60720
257 7 1 365 54820
258 26 1 364 47705
259 0 0 0 60720
260 0 1 0 60720
261 0 0 0 60720
262 0 0 0 60720
263 13 1 226 72835
264 10 1 307 58856
265 0 0 0 60720
266 0 1 0 60720
267 0 0 0 60720
268 9 1 188 77655
269 0 0 0 60720
270 26 1 138 69817
271 0 0 0 60720
272 0 0 0 60720
273 0 0 0 60720
274 19 1 125 60798
275 0 1 0 60720
276 12 1 282 62452
277 23 1 335 64175
278 0 1 0 60720
279 29 0 1324 67440
280 8 1 176 68136
281 0 0 0 60720
282 0 1 0 60720
283 26 1 249 56726
284 0 1 0 60720
285 9 1 333 70811
286 0 1 0 60720
287 5 0 601 60720
288 3 1 30 62045
289 0 0 0 60720
290 13 1 249 54323
291 0 0 0 60720
292 12 1 165 62841
293 19 1 453 81125
294 0 1 0 60720
295 10 1 53 59506
296 9 0 382 59365
297 0 1 0 60720
298 0 0 0 60720
299 0 1 0 60720
300 9 1 0 60720
301 4 0 30 60798
302 1 1 290 58790
303 1 1 0 60720
304 0 0 0 60720
305 14 1 366 61808
306 12 1 2 60735
307 0 0 0 60720
308 19 0 209 64016
309 17 1 384 54683
310 0 0 0 60720
311 0 0 0 60720
312 32 1 365 87192
313 0 0 0 60720
314 14 0 49 64107
315 8 1 3 60761
316 4 1 133 65990
317 0 1 32 59988
318 20 1 368 61167
319 5 1 1 60719
320 0 1 0 60720
321 0 1 0 60720
322 0 1 0 60720
323 0 1 0 60720
324 0 0 0 60720
325 0 1 0 60720
326 1 1 22 60722
327 0 0 0 60720
328 0 0 0 60720
329 0 1 0 60720
330 0 1 0 60720
331 0 1 0 60720
332 0 1 0 60720
333 0 1 0 60720
334 4 1 96 60379
335 1 1 1 60727
336 4 1 314 60720
337 20 1 844 60925
338 0 1 0 60720
339 1 1 26 60896
340 10 1 125 59734
341 12 1 304 62969
342 0 1 0 60720
343 0 1 0 60720
344 0 1 0 60720
345 13 1 621 60720
346 0 1 0 60720
347 3 1 119 59118
348 0 1 0 60720
349 0 1 0 60720
350 10 1 1595 60720
351 3 1 312 58598
352 7 1 60 61124
353 10 1 587 59595
354 1 1 135 62065
355 0 1 0 60720
356 0 1 0 60720
357 15 1 514 78780
358 0 1 0 60720
359 0 1 0 60720
360 0 1 0 60720
361 4 1 1 60722
362 0 1 0 60720
363 0 1 0 60720
364 28 0 1763 61600
365 9 1 180 59635
366 0 1 0 60720
367 0 1 0 60720
368 0 1 0 60720
369 0 1 0 60720
370 7 1 218 60720
371 0 1 0 60720
372 7 1 448 59781
373 7 1 227 76644
374 3 1 174 64820
375 0 1 0 60720
376 0 1 0 60720
377 11 1 121 56178
378 7 1 607 60436
379 10 1 2212 60720
380 0 1 0 60720
381 0 1 0 60720
382 18 1 530 73433
383 14 1 571 41477
384 0 1 0 60720
385 12 1 78 62700
386 29 1 2489 67804
387 3 1 131 59661
388 6 1 923 58620
389 3 1 72 60398
390 8 1 572 58580
391 10 1 397 62710
392 6 1 450 59325
393 8 1 622 60950
394 6 1 694 68060
395 9 0 3425 83620
396 8 1 562 58456
397 26 1 4917 52811
398 239 0 1442 121173
399 7 1 529 63870
400 41 0 2126 21001
401 3 1 1061 70415
402 8 1 776 64230
403 6 1 611 59190
404 21 0 1526 69351
405 7 1 592 64270
406 11 1 1182 70694
407 11 1 621 68005
408 12 1 989 58930
409 9 1 438 58320
410 3 1 726 69980
411 57 1 1303 69863
412 21 0 7419 63255
413 15 0 1164 57320
414 32 0 3310 75230
415 11 1 1920 79420
416 2 1 965 73490
417 23 1 3256 35250
418 20 0 1135 62285
419 24 1 1270 69206
420 1 1 661 65920
421 1 1 1013 69770
422 74 1 2844 72683
423 68 0 11528 -14545
424 20 1 6526 55830
425 20 1 2264 55174
426 82 0 5109 67038
427 21 1 3999 51252
428 244 1 35624 157278
429 32 1 9252 79510
430 86 1 15236 77440
431 69 1 18073 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs `Dividends\r`
3.6872375 -7.7818932 0.0052650 0.0002373
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-103.121 -10.315 -7.962 1.543 294.875
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.687e+00 5.309e+00 0.694 0.4878
Group -7.782e+00 3.837e+00 -2.028 0.0432 *
Costs 5.265e-03 1.522e-04 34.598 < 2e-16 ***
`Dividends\r` 2.373e-04 5.564e-05 4.265 2.46e-05 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 34.92 on 427 degrees of freedom
Multiple R-squared: 0.7818, Adjusted R-squared: 0.7802
F-statistic: 509.9 on 3 and 427 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.992110e-12 3.496055e-12
[2,] 1.0000000 1.783766e-11 8.918829e-12
[3,] 1.0000000 4.509206e-12 2.254603e-12
[4,] 1.0000000 3.262580e-12 1.631290e-12
[5,] 1.0000000 1.645105e-12 8.225525e-13
[6,] 1.0000000 5.441343e-12 2.720671e-12
[7,] 1.0000000 9.191492e-21 4.595746e-21
[8,] 1.0000000 1.521241e-22 7.606205e-23
[9,] 1.0000000 8.805669e-23 4.402835e-23
[10,] 1.0000000 2.834911e-22 1.417455e-22
[11,] 1.0000000 1.772261e-23 8.861304e-24
[12,] 1.0000000 9.183352e-24 4.591676e-24
[13,] 1.0000000 1.073623e-24 5.368117e-25
[14,] 1.0000000 4.910575e-30 2.455287e-30
[15,] 1.0000000 6.839123e-33 3.419561e-33
[16,] 1.0000000 6.494326e-39 3.247163e-39
[17,] 1.0000000 1.351486e-38 6.757429e-39
[18,] 1.0000000 3.749504e-46 1.874752e-46
[19,] 1.0000000 1.098941e-45 5.494704e-46
[20,] 1.0000000 4.529107e-48 2.264554e-48
[21,] 1.0000000 1.602584e-47 8.012920e-48
[22,] 1.0000000 9.495726e-48 4.747863e-48
[23,] 1.0000000 1.471032e-48 7.355160e-49
[24,] 1.0000000 3.708728e-48 1.854364e-48
[25,] 1.0000000 1.842286e-49 9.211431e-50
[26,] 1.0000000 2.316757e-50 1.158378e-50
[27,] 1.0000000 3.267304e-51 1.633652e-51
[28,] 1.0000000 2.439680e-51 1.219840e-51
[29,] 1.0000000 1.039451e-53 5.197257e-54
[30,] 1.0000000 7.887344e-55 3.943672e-55
[31,] 1.0000000 1.676388e-54 8.381938e-55
[32,] 1.0000000 5.273186e-54 2.636593e-54
[33,] 1.0000000 2.030178e-53 1.015089e-53
[34,] 1.0000000 2.262117e-55 1.131059e-55
[35,] 1.0000000 2.683886e-60 1.341943e-60
[36,] 1.0000000 1.047815e-61 5.239076e-62
[37,] 1.0000000 1.499437e-61 7.497183e-62
[38,] 1.0000000 5.680464e-61 2.840232e-61
[39,] 1.0000000 8.189325e-61 4.094662e-61
[40,] 1.0000000 2.776572e-61 1.388286e-61
[41,] 1.0000000 3.634413e-62 1.817206e-62
[42,] 1.0000000 5.807807e-63 2.903903e-63
[43,] 1.0000000 1.799503e-62 8.997515e-63
[44,] 1.0000000 4.109309e-62 2.054654e-62
[45,] 1.0000000 4.860539e-62 2.430269e-62
[46,] 1.0000000 1.811006e-63 9.055031e-64
[47,] 1.0000000 5.780898e-63 2.890449e-63
[48,] 1.0000000 1.977694e-62 9.888469e-63
[49,] 1.0000000 2.245873e-62 1.122936e-62
[50,] 1.0000000 6.333041e-63 3.166521e-63
[51,] 1.0000000 1.627810e-63 8.139049e-64
[52,] 1.0000000 3.783842e-63 1.891921e-63
[53,] 1.0000000 9.721985e-63 4.860992e-63
[54,] 1.0000000 2.226015e-62 1.113007e-62
[55,] 1.0000000 1.868257e-62 9.341286e-63
[56,] 1.0000000 3.656508e-62 1.828254e-62
[57,] 1.0000000 2.421168e-62 1.210584e-62
[58,] 1.0000000 5.961729e-62 2.980864e-62
[59,] 1.0000000 8.574549e-62 4.287275e-62
[60,] 1.0000000 3.099402e-74 1.549701e-74
[61,] 1.0000000 5.386321e-74 2.693161e-74
[62,] 1.0000000 1.354820e-73 6.774098e-74
[63,] 1.0000000 1.736474e-73 8.682370e-74
[64,] 1.0000000 1.314838e-73 6.574192e-74
[65,] 1.0000000 6.193171e-74 3.096586e-74
[66,] 1.0000000 1.012983e-73 5.064913e-74
[67,] 1.0000000 3.249078e-73 1.624539e-73
[68,] 1.0000000 5.239363e-73 2.619682e-73
[69,] 1.0000000 1.560232e-72 7.801161e-73
[70,] 1.0000000 2.026410e-72 1.013205e-72
[71,] 1.0000000 6.174976e-72 3.087488e-72
[72,] 1.0000000 1.161566e-71 5.807830e-72
[73,] 1.0000000 1.022715e-71 5.113575e-72
[74,] 1.0000000 8.489627e-72 4.244814e-72
[75,] 1.0000000 1.875998e-71 9.379989e-72
[76,] 1.0000000 4.215616e-71 2.107808e-71
[77,] 1.0000000 1.259376e-70 6.296881e-71
[78,] 1.0000000 2.195866e-70 1.097933e-70
[79,] 1.0000000 4.708007e-70 2.354003e-70
[80,] 1.0000000 6.457538e-70 3.228769e-70
[81,] 1.0000000 1.562963e-71 7.814817e-72
[82,] 1.0000000 2.946854e-71 1.473427e-71
[83,] 1.0000000 1.002246e-70 5.011231e-71
[84,] 1.0000000 2.025366e-70 1.012683e-70
[85,] 1.0000000 6.643732e-70 3.321866e-70
[86,] 1.0000000 1.584116e-69 7.920578e-70
[87,] 1.0000000 2.769924e-70 1.384962e-70
[88,] 1.0000000 2.046323e-71 1.023162e-71
[89,] 1.0000000 9.532588e-72 4.766294e-72
[90,] 1.0000000 1.833001e-71 9.165006e-72
[91,] 1.0000000 1.936678e-71 9.683388e-72
[92,] 1.0000000 4.098143e-71 2.049072e-71
[93,] 1.0000000 1.228848e-71 6.144241e-72
[94,] 1.0000000 1.450964e-71 7.254819e-72
[95,] 1.0000000 3.822171e-71 1.911085e-71
[96,] 1.0000000 9.711186e-71 4.855593e-71
[97,] 1.0000000 3.508204e-70 1.754102e-70
[98,] 1.0000000 1.910363e-70 9.551816e-71
[99,] 1.0000000 5.375739e-70 2.687870e-70
[100,] 1.0000000 1.430172e-69 7.150858e-70
[101,] 1.0000000 1.995901e-69 9.979507e-70
[102,] 1.0000000 1.499535e-69 7.497675e-70
[103,] 1.0000000 3.308274e-69 1.654137e-69
[104,] 1.0000000 6.895805e-69 3.447903e-69
[105,] 1.0000000 1.950584e-68 9.752918e-69
[106,] 1.0000000 3.380696e-68 1.690348e-68
[107,] 1.0000000 9.982100e-68 4.991050e-68
[108,] 1.0000000 1.257560e-67 6.287802e-68
[109,] 1.0000000 2.063069e-67 1.031534e-67
[110,] 1.0000000 2.962667e-67 1.481334e-67
[111,] 1.0000000 5.234343e-67 2.617171e-67
[112,] 1.0000000 7.744554e-67 3.872277e-67
[113,] 1.0000000 8.320241e-67 4.160120e-67
[114,] 1.0000000 2.150073e-66 1.075036e-66
[115,] 1.0000000 7.086468e-66 3.543234e-66
[116,] 1.0000000 1.867668e-65 9.338342e-66
[117,] 1.0000000 3.870475e-69 1.935237e-69
[118,] 1.0000000 4.874480e-69 2.437240e-69
[119,] 1.0000000 1.014539e-68 5.072693e-69
[120,] 1.0000000 1.883724e-68 9.418622e-69
[121,] 1.0000000 4.263225e-68 2.131613e-68
[122,] 1.0000000 1.171018e-67 5.855090e-68
[123,] 1.0000000 2.366865e-67 1.183433e-67
[124,] 1.0000000 5.496431e-67 2.748216e-67
[125,] 1.0000000 1.502610e-66 7.513050e-67
[126,] 1.0000000 3.328682e-66 1.664341e-66
[127,] 1.0000000 7.905247e-66 3.952624e-66
[128,] 1.0000000 1.834922e-65 9.174610e-66
[129,] 1.0000000 4.502132e-65 2.251066e-65
[130,] 1.0000000 1.166524e-64 5.832620e-65
[131,] 1.0000000 3.294386e-64 1.647193e-64
[132,] 1.0000000 8.289443e-64 4.144722e-64
[133,] 1.0000000 1.994589e-63 9.972947e-64
[134,] 1.0000000 4.931471e-63 2.465735e-63
[135,] 1.0000000 1.376827e-62 6.884135e-63
[136,] 1.0000000 3.431935e-62 1.715968e-62
[137,] 1.0000000 9.510795e-62 4.755397e-62
[138,] 1.0000000 2.604245e-61 1.302123e-61
[139,] 1.0000000 7.497224e-61 3.748612e-61
[140,] 1.0000000 2.062800e-60 1.031400e-60
[141,] 1.0000000 5.199464e-60 2.599732e-60
[142,] 1.0000000 1.474569e-59 7.372844e-60
[143,] 1.0000000 3.789851e-59 1.894926e-59
[144,] 1.0000000 9.556606e-59 4.778303e-59
[145,] 1.0000000 2.650627e-58 1.325314e-58
[146,] 1.0000000 1.931332e-58 9.656662e-59
[147,] 1.0000000 5.393672e-58 2.696836e-58
[148,] 1.0000000 1.497786e-57 7.488930e-58
[149,] 1.0000000 3.981498e-57 1.990749e-57
[150,] 1.0000000 1.014852e-56 5.074259e-57
[151,] 1.0000000 2.588848e-56 1.294424e-56
[152,] 1.0000000 6.801085e-56 3.400542e-56
[153,] 1.0000000 1.842844e-55 9.214220e-56
[154,] 1.0000000 4.702612e-55 2.351306e-55
[155,] 1.0000000 1.200097e-54 6.000485e-55
[156,] 1.0000000 3.062293e-54 1.531146e-54
[157,] 1.0000000 7.811966e-54 3.905983e-54
[158,] 1.0000000 2.153295e-53 1.076648e-53
[159,] 1.0000000 4.774363e-53 2.387182e-53
[160,] 1.0000000 1.234882e-52 6.174408e-53
[161,] 1.0000000 3.136989e-52 1.568494e-52
[162,] 1.0000000 5.412707e-52 2.706353e-52
[163,] 1.0000000 6.528112e-52 3.264056e-52
[164,] 1.0000000 1.657910e-51 8.289548e-52
[165,] 1.0000000 4.206196e-51 2.103098e-51
[166,] 1.0000000 1.117860e-50 5.589298e-51
[167,] 1.0000000 2.828150e-50 1.414075e-50
[168,] 1.0000000 7.791565e-50 3.895783e-50
[169,] 1.0000000 1.964139e-49 9.820694e-50
[170,] 1.0000000 5.357586e-49 2.678793e-49
[171,] 1.0000000 1.999787e-56 9.998934e-57
[172,] 1.0000000 4.064127e-56 2.032064e-56
[173,] 1.0000000 1.104421e-55 5.522104e-56
[174,] 1.0000000 2.996943e-55 1.498471e-55
[175,] 1.0000000 9.096166e-55 4.548083e-55
[176,] 1.0000000 2.118078e-54 1.059039e-54
[177,] 1.0000000 5.195834e-54 2.597917e-54
[178,] 1.0000000 1.416881e-53 7.084405e-54
[179,] 1.0000000 3.423805e-53 1.711903e-53
[180,] 1.0000000 8.965101e-53 4.482550e-53
[181,] 1.0000000 2.431945e-52 1.215973e-52
[182,] 1.0000000 6.928039e-52 3.464019e-52
[183,] 1.0000000 1.897589e-51 9.487945e-52
[184,] 1.0000000 5.030297e-51 2.515148e-51
[185,] 1.0000000 1.414732e-50 7.073661e-51
[186,] 1.0000000 3.734060e-50 1.867030e-50
[187,] 1.0000000 1.052504e-49 5.262522e-50
[188,] 1.0000000 2.764914e-49 1.382457e-49
[189,] 1.0000000 7.579097e-49 3.789548e-49
[190,] 1.0000000 1.595157e-48 7.975784e-49
[191,] 1.0000000 2.734595e-48 1.367298e-48
[192,] 1.0000000 5.783236e-48 2.891618e-48
[193,] 1.0000000 1.227474e-47 6.137372e-48
[194,] 1.0000000 3.188136e-47 1.594068e-47
[195,] 1.0000000 6.774001e-47 3.387000e-47
[196,] 1.0000000 1.441271e-46 7.206354e-47
[197,] 1.0000000 3.068292e-46 1.534146e-46
[198,] 1.0000000 7.419131e-46 3.709565e-46
[199,] 1.0000000 1.124032e-45 5.620160e-46
[200,] 1.0000000 2.750786e-45 1.375393e-45
[201,] 1.0000000 5.550674e-45 2.775337e-45
[202,] 1.0000000 1.173340e-44 5.866699e-45
[203,] 1.0000000 1.945521e-44 9.727607e-45
[204,] 1.0000000 4.396808e-44 2.198404e-44
[205,] 1.0000000 1.115967e-43 5.579833e-44
[206,] 1.0000000 2.825116e-43 1.412558e-43
[207,] 1.0000000 7.091953e-43 3.545977e-43
[208,] 1.0000000 1.892853e-42 9.464265e-43
[209,] 1.0000000 4.856959e-42 2.428480e-42
[210,] 1.0000000 1.253702e-41 6.268508e-42
[211,] 1.0000000 3.303322e-41 1.651661e-41
[212,] 1.0000000 8.655479e-41 4.327739e-41
[213,] 1.0000000 2.211852e-40 1.105926e-40
[214,] 1.0000000 3.916902e-40 1.958451e-40
[215,] 1.0000000 9.916675e-40 4.958338e-40
[216,] 1.0000000 2.443009e-39 1.221505e-39
[217,] 1.0000000 6.095579e-39 3.047790e-39
[218,] 1.0000000 1.269484e-38 6.347421e-39
[219,] 1.0000000 3.144292e-38 1.572146e-38
[220,] 1.0000000 5.895837e-38 2.947919e-38
[221,] 1.0000000 1.480531e-37 7.402655e-38
[222,] 1.0000000 3.590873e-37 1.795437e-37
[223,] 1.0000000 8.423045e-37 4.211522e-37
[224,] 1.0000000 2.029844e-36 1.014922e-36
[225,] 1.0000000 5.094191e-36 2.547095e-36
[226,] 1.0000000 1.274297e-35 6.371483e-36
[227,] 1.0000000 3.039284e-35 1.519642e-35
[228,] 1.0000000 7.310774e-35 3.655387e-35
[229,] 1.0000000 1.732236e-34 8.661181e-35
[230,] 1.0000000 4.088997e-34 2.044498e-34
[231,] 1.0000000 9.622425e-34 4.811212e-34
[232,] 1.0000000 2.351929e-33 1.175964e-33
[233,] 1.0000000 5.495746e-33 2.747873e-33
[234,] 1.0000000 1.279661e-32 6.398306e-33
[235,] 1.0000000 2.969038e-32 1.484519e-32
[236,] 1.0000000 7.048784e-32 3.524392e-32
[237,] 1.0000000 1.623483e-31 8.117413e-32
[238,] 1.0000000 3.856525e-31 1.928263e-31
[239,] 1.0000000 4.411217e-31 2.205608e-31
[240,] 1.0000000 8.623265e-31 4.311633e-31
[241,] 1.0000000 1.677939e-30 8.389696e-31
[242,] 1.0000000 3.669133e-30 1.834567e-30
[243,] 1.0000000 7.084032e-30 3.542016e-30
[244,] 1.0000000 1.359775e-29 6.798877e-30
[245,] 1.0000000 2.910427e-29 1.455213e-29
[246,] 1.0000000 5.536734e-29 2.768367e-29
[247,] 1.0000000 1.201548e-28 6.007740e-29
[248,] 1.0000000 1.632488e-28 8.162442e-29
[249,] 1.0000000 1.773064e-28 8.865322e-29
[250,] 1.0000000 3.975273e-28 1.987636e-28
[251,] 1.0000000 8.922528e-28 4.461264e-28
[252,] 1.0000000 7.428240e-28 3.714120e-28
[253,] 1.0000000 1.394616e-27 6.973078e-28
[254,] 1.0000000 3.111812e-27 1.555906e-27
[255,] 1.0000000 5.780581e-27 2.890291e-27
[256,] 1.0000000 1.064468e-26 5.322338e-27
[257,] 1.0000000 2.454721e-26 1.227360e-26
[258,] 1.0000000 5.402800e-26 2.701400e-26
[259,] 1.0000000 9.799723e-26 4.899861e-26
[260,] 1.0000000 2.143780e-25 1.071890e-25
[261,] 1.0000000 3.838445e-25 1.919222e-25
[262,] 1.0000000 7.947439e-25 3.973719e-25
[263,] 1.0000000 1.403131e-24 7.015655e-25
[264,] 1.0000000 2.502649e-24 1.251324e-24
[265,] 1.0000000 4.357512e-24 2.178756e-24
[266,] 1.0000000 7.484381e-24 3.742191e-24
[267,] 1.0000000 1.266300e-23 6.331501e-24
[268,] 1.0000000 2.255312e-23 1.127656e-23
[269,] 1.0000000 4.816236e-23 2.408118e-23
[270,] 1.0000000 1.029760e-22 5.148800e-23
[271,] 1.0000000 1.787494e-22 8.937470e-23
[272,] 1.0000000 3.779232e-22 1.889616e-22
[273,] 1.0000000 8.227022e-22 4.113511e-22
[274,] 1.0000000 1.766401e-21 8.832007e-22
[275,] 1.0000000 2.904006e-21 1.452003e-21
[276,] 1.0000000 6.043175e-21 3.021588e-21
[277,] 1.0000000 7.137288e-21 3.568644e-21
[278,] 1.0000000 1.481721e-20 7.408605e-21
[279,] 1.0000000 3.098226e-20 1.549113e-20
[280,] 1.0000000 6.373586e-20 3.186793e-20
[281,] 1.0000000 1.087795e-19 5.438976e-20
[282,] 1.0000000 2.258337e-19 1.129169e-19
[283,] 1.0000000 3.542595e-19 1.771298e-19
[284,] 1.0000000 6.289664e-19 3.144832e-19
[285,] 1.0000000 9.638186e-19 4.819093e-19
[286,] 1.0000000 1.937321e-18 9.686603e-19
[287,] 1.0000000 4.000215e-18 2.000107e-18
[288,] 1.0000000 8.004473e-18 4.002237e-18
[289,] 1.0000000 1.562749e-17 7.813744e-18
[290,] 1.0000000 2.910029e-17 1.455015e-17
[291,] 1.0000000 5.748653e-17 2.874326e-17
[292,] 1.0000000 8.395044e-17 4.197522e-17
[293,] 1.0000000 1.647761e-16 8.238805e-17
[294,] 1.0000000 3.215049e-16 1.607524e-16
[295,] 1.0000000 5.167018e-16 2.583509e-16
[296,] 1.0000000 1.006685e-15 5.033427e-16
[297,] 1.0000000 1.957017e-15 9.785084e-16
[298,] 1.0000000 2.690196e-15 1.345098e-15
[299,] 1.0000000 5.025930e-15 2.512965e-15
[300,] 1.0000000 9.254793e-15 4.627397e-15
[301,] 1.0000000 1.224917e-14 6.124584e-15
[302,] 1.0000000 2.315581e-14 1.157791e-14
[303,] 1.0000000 3.578035e-14 1.789018e-14
[304,] 1.0000000 4.469414e-14 2.234707e-14
[305,] 1.0000000 5.284523e-14 2.642261e-14
[306,] 1.0000000 9.813493e-14 4.906747e-14
[307,] 1.0000000 1.087161e-13 5.435804e-14
[308,] 1.0000000 1.707627e-13 8.538135e-14
[309,] 1.0000000 3.209443e-13 1.604721e-13
[310,] 1.0000000 5.986409e-13 2.993205e-13
[311,] 1.0000000 1.112760e-12 5.563799e-13
[312,] 1.0000000 1.770611e-12 8.853057e-13
[313,] 1.0000000 3.300800e-12 1.650400e-12
[314,] 1.0000000 6.041190e-12 3.020595e-12
[315,] 1.0000000 1.099596e-11 5.497981e-12
[316,] 1.0000000 1.990366e-11 9.951828e-12
[317,] 1.0000000 3.582637e-11 1.791319e-11
[318,] 1.0000000 3.531693e-11 1.765846e-11
[319,] 1.0000000 6.338293e-11 3.169147e-11
[320,] 1.0000000 1.137002e-10 5.685012e-11
[321,] 1.0000000 1.009264e-10 5.046319e-11
[322,] 1.0000000 7.721486e-11 3.860743e-11
[323,] 1.0000000 1.384523e-10 6.922617e-11
[324,] 1.0000000 2.467961e-10 1.233980e-10
[325,] 1.0000000 4.373145e-10 2.186572e-10
[326,] 1.0000000 7.702757e-10 3.851378e-10
[327,] 1.0000000 1.348571e-09 6.742857e-10
[328,] 1.0000000 2.380483e-09 1.190241e-09
[329,] 1.0000000 4.140865e-09 2.070432e-09
[330,] 1.0000000 7.199882e-09 3.599941e-09
[331,] 1.0000000 1.135360e-08 5.676801e-09
[332,] 1.0000000 1.931184e-08 9.655919e-09
[333,] 1.0000000 3.277872e-08 1.638936e-08
[334,] 1.0000000 5.448625e-08 2.724312e-08
[335,] 1.0000000 9.094295e-08 4.547148e-08
[336,] 0.9999999 1.509330e-07 7.546651e-08
[337,] 0.9999999 2.488522e-07 1.244261e-07
[338,] 0.9999998 4.075822e-07 2.037911e-07
[339,] 0.9999997 6.572311e-07 3.286156e-07
[340,] 0.9999995 1.062815e-06 5.314073e-07
[341,] 0.9999991 1.724235e-06 8.621177e-07
[342,] 0.9999986 2.750335e-06 1.375168e-06
[343,] 0.9999978 4.356700e-06 2.178350e-06
[344,] 0.9999965 6.909053e-06 3.454526e-06
[345,] 0.9999946 1.088545e-05 5.442724e-06
[346,] 0.9999915 1.702572e-05 8.512858e-06
[347,] 0.9999869 2.626959e-05 1.313479e-05
[348,] 0.9999800 4.006691e-05 2.003345e-05
[349,] 0.9999696 6.082596e-05 3.041298e-05
[350,] 0.9999542 9.165769e-05 4.582884e-05
[351,] 0.9999317 1.366200e-04 6.831000e-05
[352,] 0.9998986 2.027828e-04 1.013914e-04
[353,] 0.9998507 2.986989e-04 1.493494e-04
[354,] 0.9997817 4.366062e-04 2.183031e-04
[355,] 0.9996804 6.392882e-04 3.196441e-04
[356,] 0.9995400 9.199195e-04 4.599597e-04
[357,] 0.9993434 1.313284e-03 6.566422e-04
[358,] 0.9991559 1.688181e-03 8.440904e-04
[359,] 0.9988074 2.385298e-03 1.192649e-03
[360,] 0.9983323 3.335351e-03 1.667675e-03
[361,] 0.9976873 4.625492e-03 2.312746e-03
[362,] 0.9968193 6.361484e-03 3.180742e-03
[363,] 0.9956621 8.675794e-03 4.337897e-03
[364,] 0.9940837 1.183253e-02 5.916267e-03
[365,] 0.9920673 1.586540e-02 7.932699e-03
[366,] 0.9893704 2.125924e-02 1.062962e-02
[367,] 0.9864092 2.718163e-02 1.359082e-02
[368,] 0.9822819 3.543628e-02 1.771814e-02
[369,] 0.9770407 4.591865e-02 2.295932e-02
[370,] 0.9705134 5.897327e-02 2.948663e-02
[371,] 0.9627919 7.441629e-02 3.720814e-02
[372,] 0.9526758 9.464843e-02 4.732421e-02
[373,] 0.9408165 1.183670e-01 5.918348e-02
[374,] 0.9266109 1.467782e-01 7.338908e-02
[375,] 0.9098331 1.803337e-01 9.016687e-02
[376,] 0.8894295 2.211409e-01 1.105705e-01
[377,] 0.8762230 2.475540e-01 1.237770e-01
[378,] 0.8512575 2.974850e-01 1.487425e-01
[379,] 0.8226052 3.547897e-01 1.773948e-01
[380,] 0.7911589 4.176822e-01 2.088411e-01
[381,] 0.7550630 4.898741e-01 2.449370e-01
[382,] 0.7155580 5.688840e-01 2.844420e-01
[383,] 0.6729005 6.541990e-01 3.270995e-01
[384,] 0.6274248 7.451503e-01 3.725752e-01
[385,] 0.5795208 8.409584e-01 4.204792e-01
[386,] 0.5299879 9.400242e-01 4.700121e-01
[387,] 0.4794956 9.589913e-01 5.205044e-01
[388,] 0.4316219 8.632437e-01 5.683781e-01
[389,] 0.5672728 8.654543e-01 4.327272e-01
[390,] 0.5149140 9.701720e-01 4.850860e-01
[391,] 0.4637006 9.274013e-01 5.362994e-01
[392,] 0.9999254 1.491608e-04 7.458038e-05
[393,] 0.9998579 2.841312e-04 1.420656e-04
[394,] 0.9998615 2.769138e-04 1.384569e-04
[395,] 0.9997647 4.706756e-04 2.353378e-04
[396,] 0.9995600 8.799155e-04 4.399578e-04
[397,] 0.9991909 1.618236e-03 8.091180e-04
[398,] 0.9985932 2.813665e-03 1.406833e-03
[399,] 0.9975167 4.966569e-03 2.483284e-03
[400,] 0.9957463 8.507382e-03 4.253691e-03
[401,] 0.9927585 1.448305e-02 7.241527e-03
[402,] 0.9879732 2.405367e-02 1.202684e-02
[403,] 0.9804978 3.900447e-02 1.950224e-02
[404,] 0.9707397 5.852053e-02 2.926026e-02
[405,] 0.9838116 3.237689e-02 1.618844e-02
[406,] 0.9915559 1.688811e-02 8.444057e-03
[407,] 0.9874102 2.517951e-02 1.258976e-02
[408,] 0.9855522 2.889565e-02 1.444782e-02
[409,] 0.9762763 4.744735e-02 2.372368e-02
[410,] 0.9635633 7.287342e-02 3.643671e-02
[411,] 0.9466721 1.066559e-01 5.332794e-02
[412,] 0.9600247 7.995055e-02 3.997527e-02
[413,] 0.9305068 1.389864e-01 6.949318e-02
[414,] 0.8876499 2.247002e-01 1.123501e-01
[415,] 0.8448250 3.103501e-01 1.551750e-01
[416,] 0.9783325 4.333496e-02 2.166748e-02
[417,] 0.9411137 1.177727e-01 5.888635e-02
[418,] 0.8662929 2.674141e-01 1.337071e-01
> postscript(file="/var/www/html/rcomp/tmp/1wgnn1291216940.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/2wgnn1291216940.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/3774q1291216940.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/4774q1291216940.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/5774q1291216940.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
-103.12144864 264.06371393 -88.99460634 -7.87939772 -6.59766571
6 7 8 9 10
-97.46474114 20.80076843 -30.32269121 11.45580057 28.27686821
11 12 13 14 15
-86.17227154 47.94377848 294.87478806 -42.97573934 123.05471294
16 17 18 19 20
-3.44092154 0.20482440 118.61198031 23.37394291 158.11423255
21 22 23 24 25
-18.60966894 -73.19279320 -2.35389504 181.23977737 17.66401481
26 27 28 29 30
-38.36851597 68.86397589 6.96004113 -14.64986439 55.36251778
31 32 33 34 35
-16.10969250 -21.09658048 -19.57853588 -9.71658559 -67.18461972
36 37 38 39 40
93.08518901 2.70226347 21.52616528 21.34913616 94.73262800
41 42 43 44 45
131.17992576 -35.01927573 47.00770795 28.04002519 2.31571418
46 47 48 49 50
-22.26665686 -30.40942241 73.37187276 7.35736507 14.93777033
51 52 53 54 55
44.02296804 78.48259622 10.13471750 18.64292493 48.62286782
56 57 58 59 60
53.98405474 -21.58741500 1.99900318 24.71884384 11.46018567
61 62 63 64 65
-12.37888920 13.73343469 -17.59123856 13.94333651 30.00512098
66 67 68 69 70
163.76242696 -4.71150432 5.67177477 38.48143101 36.26781822
71 72 73 74 75
22.79499240 5.87853997 6.19496782 -8.63771584 -0.51603113
76 77 78 79 80
-11.40111070 0.53024646 -6.73490421 32.87032938 -22.16587073
81 82 83 84 85
11.35697634 16.29885511 -0.49745998 5.40531575 -7.98979095
86 87 88 89 90
19.09944882 63.49893927 -10.41976570 12.74296952 -9.41597178
91 92 93 94 95
6.88364933 1.11209862 45.14099699 -39.20121921 28.24332426
96 97 98 99 100
-4.06937189 -15.15289012 -5.78770427 -21.28905374 -19.34521631
101 102 103 104 105
2.61056046 -3.52384460 15.62890398 -17.93523504 6.91523333
106 107 108 109 110
12.18442828 13.30804579 -21.68498505 -3.00941175 -11.17234032
111 112 113 114 115
5.04363424 -0.17942150 0.64444480 -20.87354555 -18.65308885
116 117 118 119 120
17.80530380 -10.19132013 -11.71746846 -22.77001316 6.14942522
121 122 123 124 125
14.21814203 -4.41236126 85.65578864 -12.95156869 -15.55559455
126 127 128 129 130
-16.02626991 -10.31514134 -0.82737576 11.12395059 -10.31514134
131 132 133 134 135
-3.47279309 7.68825694 -10.31514134 9.41883179 -9.71106899
136 137 138 139 140
-7.92785981 -0.97401280 -9.34146636 -11.03785002 -10.31514134
141 142 143 144 145
-3.95267213 -10.31514134 -3.64853293 -6.54597439 1.63976627
146 147 148 149 150
-1.94886225 -10.31514134 -3.38981766 -9.74195746 -10.45909679
151 152 153 154 155
-5.56170040 26.45378438 -5.31514134 -5.22933725 -8.47217620
156 157 158 159 160
-10.31514134 -10.31514134 -9.01384468 -7.47713244 -10.31514134
161 162 163 164 165
-10.31514134 -10.31514134 -10.31514134 -5.88424303 10.69152811
166 167 168 169 170
-9.56259654 -10.31514134 7.26633980 20.17626456 -10.31514134
171 172 173 174 175
-10.31514134 -7.70497577 -10.31514134 -2.67298745 -10.31514134
176 177 178 179 180
-3.03335174 100.57885702 -18.47678277 -10.31514134 -10.31514134
181 182 183 184 185
0.97531480 -12.48447433 2.89717434 -6.73398589 -13.65360363
186 187 188 189 190
-10.37304679 -8.19211529 -5.41312073 -7.69754265 -10.31514134
191 192 193 194 195
-4.47349256 -10.31514134 -4.27346403 -10.31514134 -7.50290509
196 197 198 199 200
-18.09703459 32.15274937 -18.09703459 -18.09703459 -6.16823853
201 202 203 204 205
-18.09703459 -18.09703459 -18.09703459 -4.81260995 -15.03768096
206 207 208 209 210
-13.37181194 -17.60771325 -18.09703459 4.81669519 -16.49226068
211 212 213 214 215
-10.31514134 -10.31514134 0.01003430 -3.10905887 -8.57755298
216 217 218 219 220
-8.31514134 -4.42045843 -5.38778663 -8.31514134 8.92709417
221 222 223 224 225
-8.78162627 -10.31514134 -9.32571542 7.22083674 -9.44061324
226 227 228 229 230
10.05886372 -1.63630659 -10.31514134 -10.40723154 -10.31514134
231 232 233 234 235
-5.74237066 -5.20071437 -10.31514134 0.51809038 -10.31514134
236 237 238 239 240
-10.33112961 -10.31514134 -4.50614251 -10.31514134 -10.31514134
241 242 243 244 245
-10.31514134 -7.96193365 -10.31514134 -3.22744286 -23.43142427
246 247 248 249 250
-18.09703459 -18.09703459 3.73694386 -18.09703459 -18.09703459
251 252 253 254 255
-0.59376125 -18.09703459 -11.03093416 21.49413058 29.20414036
256 257 258 259 260
-10.31514134 -3.83670640 16.85705837 -18.09703459 -10.31514134
261 262 263 264 265
-18.09703459 -18.09703459 -1.38010953 -1.48914151 -18.09703459
266 267 268 269 270
-10.31514134 -18.09703459 -6.32390008 -18.09703459 12.79942900
271 272 273 274 275
-18.09703459 -18.09703459 -18.09703459 8.00822253 -10.31514134
276 277 278 279 280
-0.21090294 10.10115724 -10.31514134 2.33733983 -5.00171380
281 282 283 284 285
-18.09703459 -10.31514134 15.32171073 -10.31514134 -5.46313840
286 287 288 289 290
-10.31514134 -16.26130206 -7.78753450 -18.09703459 2.89197988
291 292 293 294 295
-18.09703459 0.31278681 1.45738906 -10.31514134 -0.30608555
296 297 298 299 300
-10.78670366 -10.31514134 -18.09703459 -10.31514134 -1.31514134
301 302 303 304 305
-14.27349532 -10.38397362 -9.31514134 -18.09703459 1.49966789
306 307 308 309 310
1.67076892 -18.09703459 -0.97961233 6.09577074 -18.09703459
311 312 313 314 315
-18.09703459 13.48091629 -18.09703459 -5.15880739 -2.34066629
316 317 318 319 320
-8.26603956 -10.30990653 7.64125712 -5.32016903 -10.31514134
321 322 323 324 325
-10.31514134 -10.31514134 -10.31514134 -18.09703459 -10.31514134
326 327 328 329 330
-9.43144607 -18.09703459 -18.09703459 -10.31514134 -10.31514134
331 332 333 334 335
-10.31514134 -10.31514134 -10.31514134 -6.73965715 -9.32206756
336 337 338 339 340
-7.96835263 5.19254552 -10.31514134 -9.49379898 -0.73927378
341 342 343 344 345
-0.44942515 -10.31514134 -10.31514134 -10.31514134 -0.58470889
346 347 348 349 350
-10.31514134 -7.56149740 -10.31514134 -10.31514134 -8.71282289
351 352 353 354 355
-8.45423914 -3.72691705 -3.13871882 -10.34510624 -10.31514134
356 357 358 359 360
-10.31514134 -2.30727135 -10.31514134 -10.31514134 -10.31514134
361 362 363 364 365
-6.32088098 -10.31514134 -10.31514134 0.41192554 -2.00535477
366 367 368 369 370
-10.31514134 -10.31514134 -10.31514134 -10.31514134 -4.46291224
371 372 373 374 375
-10.31514134 -5.45102393 -8.28930927 -9.20424560 -10.31514134
376 377 378 379 380
-10.31514134 1.12568015 -6.44360123 -11.96133042 -10.31514134
381 382 383 384 385
-10.31514134 1.87741455 5.24520344 -10.31514134 0.80430365
386 387 388 389 390
3.89912046 -7.75353977 -8.67637758 -7.61780605 -4.81886852
391 392 393 394 395
-2.87760581 -6.35333807 -5.64455646 -9.71094995 -32.56419855
396 397 398 399 400
-4.73679136 -8.32623824 198.96439569 -6.84787233 21.13550096
401 402 403 404 405
-15.20208445 -7.23376193 -7.16896614 -7.17970092 -7.27449379
406 407 408 409 410
-7.90536097 -4.31355229 -3.09743566 -3.05165594 -13.33507583
411 412 413 414 415
37.65478271 -36.75969481 -8.41862667 -6.96764605 -13.86174903
416 417 418 419 420
-16.42638921 1.58643104 -4.44421300 4.98444413 -14.02934965
421 422 423 424 425
-16.79629576 45.87218169 7.06954909 -23.51408533 -0.91895888
426 427 428 429 430
35.50470013 -8.12298961 23.20964267 -31.48611750 -8.50065897
431
-28.53467948
> postscript(file="/var/www/html/rcomp/tmp/6zgls1291216940.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 -103.12144864 NA
1 264.06371393 -103.12144864
2 -88.99460634 264.06371393
3 -7.87939772 -88.99460634
4 -6.59766571 -7.87939772
5 -97.46474114 -6.59766571
6 20.80076843 -97.46474114
7 -30.32269121 20.80076843
8 11.45580057 -30.32269121
9 28.27686821 11.45580057
10 -86.17227154 28.27686821
11 47.94377848 -86.17227154
12 294.87478806 47.94377848
13 -42.97573934 294.87478806
14 123.05471294 -42.97573934
15 -3.44092154 123.05471294
16 0.20482440 -3.44092154
17 118.61198031 0.20482440
18 23.37394291 118.61198031
19 158.11423255 23.37394291
20 -18.60966894 158.11423255
21 -73.19279320 -18.60966894
22 -2.35389504 -73.19279320
23 181.23977737 -2.35389504
24 17.66401481 181.23977737
25 -38.36851597 17.66401481
26 68.86397589 -38.36851597
27 6.96004113 68.86397589
28 -14.64986439 6.96004113
29 55.36251778 -14.64986439
30 -16.10969250 55.36251778
31 -21.09658048 -16.10969250
32 -19.57853588 -21.09658048
33 -9.71658559 -19.57853588
34 -67.18461972 -9.71658559
35 93.08518901 -67.18461972
36 2.70226347 93.08518901
37 21.52616528 2.70226347
38 21.34913616 21.52616528
39 94.73262800 21.34913616
40 131.17992576 94.73262800
41 -35.01927573 131.17992576
42 47.00770795 -35.01927573
43 28.04002519 47.00770795
44 2.31571418 28.04002519
45 -22.26665686 2.31571418
46 -30.40942241 -22.26665686
47 73.37187276 -30.40942241
48 7.35736507 73.37187276
49 14.93777033 7.35736507
50 44.02296804 14.93777033
51 78.48259622 44.02296804
52 10.13471750 78.48259622
53 18.64292493 10.13471750
54 48.62286782 18.64292493
55 53.98405474 48.62286782
56 -21.58741500 53.98405474
57 1.99900318 -21.58741500
58 24.71884384 1.99900318
59 11.46018567 24.71884384
60 -12.37888920 11.46018567
61 13.73343469 -12.37888920
62 -17.59123856 13.73343469
63 13.94333651 -17.59123856
64 30.00512098 13.94333651
65 163.76242696 30.00512098
66 -4.71150432 163.76242696
67 5.67177477 -4.71150432
68 38.48143101 5.67177477
69 36.26781822 38.48143101
70 22.79499240 36.26781822
71 5.87853997 22.79499240
72 6.19496782 5.87853997
73 -8.63771584 6.19496782
74 -0.51603113 -8.63771584
75 -11.40111070 -0.51603113
76 0.53024646 -11.40111070
77 -6.73490421 0.53024646
78 32.87032938 -6.73490421
79 -22.16587073 32.87032938
80 11.35697634 -22.16587073
81 16.29885511 11.35697634
82 -0.49745998 16.29885511
83 5.40531575 -0.49745998
84 -7.98979095 5.40531575
85 19.09944882 -7.98979095
86 63.49893927 19.09944882
87 -10.41976570 63.49893927
88 12.74296952 -10.41976570
89 -9.41597178 12.74296952
90 6.88364933 -9.41597178
91 1.11209862 6.88364933
92 45.14099699 1.11209862
93 -39.20121921 45.14099699
94 28.24332426 -39.20121921
95 -4.06937189 28.24332426
96 -15.15289012 -4.06937189
97 -5.78770427 -15.15289012
98 -21.28905374 -5.78770427
99 -19.34521631 -21.28905374
100 2.61056046 -19.34521631
101 -3.52384460 2.61056046
102 15.62890398 -3.52384460
103 -17.93523504 15.62890398
104 6.91523333 -17.93523504
105 12.18442828 6.91523333
106 13.30804579 12.18442828
107 -21.68498505 13.30804579
108 -3.00941175 -21.68498505
109 -11.17234032 -3.00941175
110 5.04363424 -11.17234032
111 -0.17942150 5.04363424
112 0.64444480 -0.17942150
113 -20.87354555 0.64444480
114 -18.65308885 -20.87354555
115 17.80530380 -18.65308885
116 -10.19132013 17.80530380
117 -11.71746846 -10.19132013
118 -22.77001316 -11.71746846
119 6.14942522 -22.77001316
120 14.21814203 6.14942522
121 -4.41236126 14.21814203
122 85.65578864 -4.41236126
123 -12.95156869 85.65578864
124 -15.55559455 -12.95156869
125 -16.02626991 -15.55559455
126 -10.31514134 -16.02626991
127 -0.82737576 -10.31514134
128 11.12395059 -0.82737576
129 -10.31514134 11.12395059
130 -3.47279309 -10.31514134
131 7.68825694 -3.47279309
132 -10.31514134 7.68825694
133 9.41883179 -10.31514134
134 -9.71106899 9.41883179
135 -7.92785981 -9.71106899
136 -0.97401280 -7.92785981
137 -9.34146636 -0.97401280
138 -11.03785002 -9.34146636
139 -10.31514134 -11.03785002
140 -3.95267213 -10.31514134
141 -10.31514134 -3.95267213
142 -3.64853293 -10.31514134
143 -6.54597439 -3.64853293
144 1.63976627 -6.54597439
145 -1.94886225 1.63976627
146 -10.31514134 -1.94886225
147 -3.38981766 -10.31514134
148 -9.74195746 -3.38981766
149 -10.45909679 -9.74195746
150 -5.56170040 -10.45909679
151 26.45378438 -5.56170040
152 -5.31514134 26.45378438
153 -5.22933725 -5.31514134
154 -8.47217620 -5.22933725
155 -10.31514134 -8.47217620
156 -10.31514134 -10.31514134
157 -9.01384468 -10.31514134
158 -7.47713244 -9.01384468
159 -10.31514134 -7.47713244
160 -10.31514134 -10.31514134
161 -10.31514134 -10.31514134
162 -10.31514134 -10.31514134
163 -5.88424303 -10.31514134
164 10.69152811 -5.88424303
165 -9.56259654 10.69152811
166 -10.31514134 -9.56259654
167 7.26633980 -10.31514134
168 20.17626456 7.26633980
169 -10.31514134 20.17626456
170 -10.31514134 -10.31514134
171 -7.70497577 -10.31514134
172 -10.31514134 -7.70497577
173 -2.67298745 -10.31514134
174 -10.31514134 -2.67298745
175 -3.03335174 -10.31514134
176 100.57885702 -3.03335174
177 -18.47678277 100.57885702
178 -10.31514134 -18.47678277
179 -10.31514134 -10.31514134
180 0.97531480 -10.31514134
181 -12.48447433 0.97531480
182 2.89717434 -12.48447433
183 -6.73398589 2.89717434
184 -13.65360363 -6.73398589
185 -10.37304679 -13.65360363
186 -8.19211529 -10.37304679
187 -5.41312073 -8.19211529
188 -7.69754265 -5.41312073
189 -10.31514134 -7.69754265
190 -4.47349256 -10.31514134
191 -10.31514134 -4.47349256
192 -4.27346403 -10.31514134
193 -10.31514134 -4.27346403
194 -7.50290509 -10.31514134
195 -18.09703459 -7.50290509
196 32.15274937 -18.09703459
197 -18.09703459 32.15274937
198 -18.09703459 -18.09703459
199 -6.16823853 -18.09703459
200 -18.09703459 -6.16823853
201 -18.09703459 -18.09703459
202 -18.09703459 -18.09703459
203 -4.81260995 -18.09703459
204 -15.03768096 -4.81260995
205 -13.37181194 -15.03768096
206 -17.60771325 -13.37181194
207 -18.09703459 -17.60771325
208 4.81669519 -18.09703459
209 -16.49226068 4.81669519
210 -10.31514134 -16.49226068
211 -10.31514134 -10.31514134
212 0.01003430 -10.31514134
213 -3.10905887 0.01003430
214 -8.57755298 -3.10905887
215 -8.31514134 -8.57755298
216 -4.42045843 -8.31514134
217 -5.38778663 -4.42045843
218 -8.31514134 -5.38778663
219 8.92709417 -8.31514134
220 -8.78162627 8.92709417
221 -10.31514134 -8.78162627
222 -9.32571542 -10.31514134
223 7.22083674 -9.32571542
224 -9.44061324 7.22083674
225 10.05886372 -9.44061324
226 -1.63630659 10.05886372
227 -10.31514134 -1.63630659
228 -10.40723154 -10.31514134
229 -10.31514134 -10.40723154
230 -5.74237066 -10.31514134
231 -5.20071437 -5.74237066
232 -10.31514134 -5.20071437
233 0.51809038 -10.31514134
234 -10.31514134 0.51809038
235 -10.33112961 -10.31514134
236 -10.31514134 -10.33112961
237 -4.50614251 -10.31514134
238 -10.31514134 -4.50614251
239 -10.31514134 -10.31514134
240 -10.31514134 -10.31514134
241 -7.96193365 -10.31514134
242 -10.31514134 -7.96193365
243 -3.22744286 -10.31514134
244 -23.43142427 -3.22744286
245 -18.09703459 -23.43142427
246 -18.09703459 -18.09703459
247 3.73694386 -18.09703459
248 -18.09703459 3.73694386
249 -18.09703459 -18.09703459
250 -0.59376125 -18.09703459
251 -18.09703459 -0.59376125
252 -11.03093416 -18.09703459
253 21.49413058 -11.03093416
254 29.20414036 21.49413058
255 -10.31514134 29.20414036
256 -3.83670640 -10.31514134
257 16.85705837 -3.83670640
258 -18.09703459 16.85705837
259 -10.31514134 -18.09703459
260 -18.09703459 -10.31514134
261 -18.09703459 -18.09703459
262 -1.38010953 -18.09703459
263 -1.48914151 -1.38010953
264 -18.09703459 -1.48914151
265 -10.31514134 -18.09703459
266 -18.09703459 -10.31514134
267 -6.32390008 -18.09703459
268 -18.09703459 -6.32390008
269 12.79942900 -18.09703459
270 -18.09703459 12.79942900
271 -18.09703459 -18.09703459
272 -18.09703459 -18.09703459
273 8.00822253 -18.09703459
274 -10.31514134 8.00822253
275 -0.21090294 -10.31514134
276 10.10115724 -0.21090294
277 -10.31514134 10.10115724
278 2.33733983 -10.31514134
279 -5.00171380 2.33733983
280 -18.09703459 -5.00171380
281 -10.31514134 -18.09703459
282 15.32171073 -10.31514134
283 -10.31514134 15.32171073
284 -5.46313840 -10.31514134
285 -10.31514134 -5.46313840
286 -16.26130206 -10.31514134
287 -7.78753450 -16.26130206
288 -18.09703459 -7.78753450
289 2.89197988 -18.09703459
290 -18.09703459 2.89197988
291 0.31278681 -18.09703459
292 1.45738906 0.31278681
293 -10.31514134 1.45738906
294 -0.30608555 -10.31514134
295 -10.78670366 -0.30608555
296 -10.31514134 -10.78670366
297 -18.09703459 -10.31514134
298 -10.31514134 -18.09703459
299 -1.31514134 -10.31514134
300 -14.27349532 -1.31514134
301 -10.38397362 -14.27349532
302 -9.31514134 -10.38397362
303 -18.09703459 -9.31514134
304 1.49966789 -18.09703459
305 1.67076892 1.49966789
306 -18.09703459 1.67076892
307 -0.97961233 -18.09703459
308 6.09577074 -0.97961233
309 -18.09703459 6.09577074
310 -18.09703459 -18.09703459
311 13.48091629 -18.09703459
312 -18.09703459 13.48091629
313 -5.15880739 -18.09703459
314 -2.34066629 -5.15880739
315 -8.26603956 -2.34066629
316 -10.30990653 -8.26603956
317 7.64125712 -10.30990653
318 -5.32016903 7.64125712
319 -10.31514134 -5.32016903
320 -10.31514134 -10.31514134
321 -10.31514134 -10.31514134
322 -10.31514134 -10.31514134
323 -18.09703459 -10.31514134
324 -10.31514134 -18.09703459
325 -9.43144607 -10.31514134
326 -18.09703459 -9.43144607
327 -18.09703459 -18.09703459
328 -10.31514134 -18.09703459
329 -10.31514134 -10.31514134
330 -10.31514134 -10.31514134
331 -10.31514134 -10.31514134
332 -10.31514134 -10.31514134
333 -6.73965715 -10.31514134
334 -9.32206756 -6.73965715
335 -7.96835263 -9.32206756
336 5.19254552 -7.96835263
337 -10.31514134 5.19254552
338 -9.49379898 -10.31514134
339 -0.73927378 -9.49379898
340 -0.44942515 -0.73927378
341 -10.31514134 -0.44942515
342 -10.31514134 -10.31514134
343 -10.31514134 -10.31514134
344 -0.58470889 -10.31514134
345 -10.31514134 -0.58470889
346 -7.56149740 -10.31514134
347 -10.31514134 -7.56149740
348 -10.31514134 -10.31514134
349 -8.71282289 -10.31514134
350 -8.45423914 -8.71282289
351 -3.72691705 -8.45423914
352 -3.13871882 -3.72691705
353 -10.34510624 -3.13871882
354 -10.31514134 -10.34510624
355 -10.31514134 -10.31514134
356 -2.30727135 -10.31514134
357 -10.31514134 -2.30727135
358 -10.31514134 -10.31514134
359 -10.31514134 -10.31514134
360 -6.32088098 -10.31514134
361 -10.31514134 -6.32088098
362 -10.31514134 -10.31514134
363 0.41192554 -10.31514134
364 -2.00535477 0.41192554
365 -10.31514134 -2.00535477
366 -10.31514134 -10.31514134
367 -10.31514134 -10.31514134
368 -10.31514134 -10.31514134
369 -4.46291224 -10.31514134
370 -10.31514134 -4.46291224
371 -5.45102393 -10.31514134
372 -8.28930927 -5.45102393
373 -9.20424560 -8.28930927
374 -10.31514134 -9.20424560
375 -10.31514134 -10.31514134
376 1.12568015 -10.31514134
377 -6.44360123 1.12568015
378 -11.96133042 -6.44360123
379 -10.31514134 -11.96133042
380 -10.31514134 -10.31514134
381 1.87741455 -10.31514134
382 5.24520344 1.87741455
383 -10.31514134 5.24520344
384 0.80430365 -10.31514134
385 3.89912046 0.80430365
386 -7.75353977 3.89912046
387 -8.67637758 -7.75353977
388 -7.61780605 -8.67637758
389 -4.81886852 -7.61780605
390 -2.87760581 -4.81886852
391 -6.35333807 -2.87760581
392 -5.64455646 -6.35333807
393 -9.71094995 -5.64455646
394 -32.56419855 -9.71094995
395 -4.73679136 -32.56419855
396 -8.32623824 -4.73679136
397 198.96439569 -8.32623824
398 -6.84787233 198.96439569
399 21.13550096 -6.84787233
400 -15.20208445 21.13550096
401 -7.23376193 -15.20208445
402 -7.16896614 -7.23376193
403 -7.17970092 -7.16896614
404 -7.27449379 -7.17970092
405 -7.90536097 -7.27449379
406 -4.31355229 -7.90536097
407 -3.09743566 -4.31355229
408 -3.05165594 -3.09743566
409 -13.33507583 -3.05165594
410 37.65478271 -13.33507583
411 -36.75969481 37.65478271
412 -8.41862667 -36.75969481
413 -6.96764605 -8.41862667
414 -13.86174903 -6.96764605
415 -16.42638921 -13.86174903
416 1.58643104 -16.42638921
417 -4.44421300 1.58643104
418 4.98444413 -4.44421300
419 -14.02934965 4.98444413
420 -16.79629576 -14.02934965
421 45.87218169 -16.79629576
422 7.06954909 45.87218169
423 -23.51408533 7.06954909
424 -0.91895888 -23.51408533
425 35.50470013 -0.91895888
426 -8.12298961 35.50470013
427 23.20964267 -8.12298961
428 -31.48611750 23.20964267
429 -8.50065897 -31.48611750
430 -28.53467948 -8.50065897
431 NA -28.53467948
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 264.06371393 -103.12144864
[2,] -88.99460634 264.06371393
[3,] -7.87939772 -88.99460634
[4,] -6.59766571 -7.87939772
[5,] -97.46474114 -6.59766571
[6,] 20.80076843 -97.46474114
[7,] -30.32269121 20.80076843
[8,] 11.45580057 -30.32269121
[9,] 28.27686821 11.45580057
[10,] -86.17227154 28.27686821
[11,] 47.94377848 -86.17227154
[12,] 294.87478806 47.94377848
[13,] -42.97573934 294.87478806
[14,] 123.05471294 -42.97573934
[15,] -3.44092154 123.05471294
[16,] 0.20482440 -3.44092154
[17,] 118.61198031 0.20482440
[18,] 23.37394291 118.61198031
[19,] 158.11423255 23.37394291
[20,] -18.60966894 158.11423255
[21,] -73.19279320 -18.60966894
[22,] -2.35389504 -73.19279320
[23,] 181.23977737 -2.35389504
[24,] 17.66401481 181.23977737
[25,] -38.36851597 17.66401481
[26,] 68.86397589 -38.36851597
[27,] 6.96004113 68.86397589
[28,] -14.64986439 6.96004113
[29,] 55.36251778 -14.64986439
[30,] -16.10969250 55.36251778
[31,] -21.09658048 -16.10969250
[32,] -19.57853588 -21.09658048
[33,] -9.71658559 -19.57853588
[34,] -67.18461972 -9.71658559
[35,] 93.08518901 -67.18461972
[36,] 2.70226347 93.08518901
[37,] 21.52616528 2.70226347
[38,] 21.34913616 21.52616528
[39,] 94.73262800 21.34913616
[40,] 131.17992576 94.73262800
[41,] -35.01927573 131.17992576
[42,] 47.00770795 -35.01927573
[43,] 28.04002519 47.00770795
[44,] 2.31571418 28.04002519
[45,] -22.26665686 2.31571418
[46,] -30.40942241 -22.26665686
[47,] 73.37187276 -30.40942241
[48,] 7.35736507 73.37187276
[49,] 14.93777033 7.35736507
[50,] 44.02296804 14.93777033
[51,] 78.48259622 44.02296804
[52,] 10.13471750 78.48259622
[53,] 18.64292493 10.13471750
[54,] 48.62286782 18.64292493
[55,] 53.98405474 48.62286782
[56,] -21.58741500 53.98405474
[57,] 1.99900318 -21.58741500
[58,] 24.71884384 1.99900318
[59,] 11.46018567 24.71884384
[60,] -12.37888920 11.46018567
[61,] 13.73343469 -12.37888920
[62,] -17.59123856 13.73343469
[63,] 13.94333651 -17.59123856
[64,] 30.00512098 13.94333651
[65,] 163.76242696 30.00512098
[66,] -4.71150432 163.76242696
[67,] 5.67177477 -4.71150432
[68,] 38.48143101 5.67177477
[69,] 36.26781822 38.48143101
[70,] 22.79499240 36.26781822
[71,] 5.87853997 22.79499240
[72,] 6.19496782 5.87853997
[73,] -8.63771584 6.19496782
[74,] -0.51603113 -8.63771584
[75,] -11.40111070 -0.51603113
[76,] 0.53024646 -11.40111070
[77,] -6.73490421 0.53024646
[78,] 32.87032938 -6.73490421
[79,] -22.16587073 32.87032938
[80,] 11.35697634 -22.16587073
[81,] 16.29885511 11.35697634
[82,] -0.49745998 16.29885511
[83,] 5.40531575 -0.49745998
[84,] -7.98979095 5.40531575
[85,] 19.09944882 -7.98979095
[86,] 63.49893927 19.09944882
[87,] -10.41976570 63.49893927
[88,] 12.74296952 -10.41976570
[89,] -9.41597178 12.74296952
[90,] 6.88364933 -9.41597178
[91,] 1.11209862 6.88364933
[92,] 45.14099699 1.11209862
[93,] -39.20121921 45.14099699
[94,] 28.24332426 -39.20121921
[95,] -4.06937189 28.24332426
[96,] -15.15289012 -4.06937189
[97,] -5.78770427 -15.15289012
[98,] -21.28905374 -5.78770427
[99,] -19.34521631 -21.28905374
[100,] 2.61056046 -19.34521631
[101,] -3.52384460 2.61056046
[102,] 15.62890398 -3.52384460
[103,] -17.93523504 15.62890398
[104,] 6.91523333 -17.93523504
[105,] 12.18442828 6.91523333
[106,] 13.30804579 12.18442828
[107,] -21.68498505 13.30804579
[108,] -3.00941175 -21.68498505
[109,] -11.17234032 -3.00941175
[110,] 5.04363424 -11.17234032
[111,] -0.17942150 5.04363424
[112,] 0.64444480 -0.17942150
[113,] -20.87354555 0.64444480
[114,] -18.65308885 -20.87354555
[115,] 17.80530380 -18.65308885
[116,] -10.19132013 17.80530380
[117,] -11.71746846 -10.19132013
[118,] -22.77001316 -11.71746846
[119,] 6.14942522 -22.77001316
[120,] 14.21814203 6.14942522
[121,] -4.41236126 14.21814203
[122,] 85.65578864 -4.41236126
[123,] -12.95156869 85.65578864
[124,] -15.55559455 -12.95156869
[125,] -16.02626991 -15.55559455
[126,] -10.31514134 -16.02626991
[127,] -0.82737576 -10.31514134
[128,] 11.12395059 -0.82737576
[129,] -10.31514134 11.12395059
[130,] -3.47279309 -10.31514134
[131,] 7.68825694 -3.47279309
[132,] -10.31514134 7.68825694
[133,] 9.41883179 -10.31514134
[134,] -9.71106899 9.41883179
[135,] -7.92785981 -9.71106899
[136,] -0.97401280 -7.92785981
[137,] -9.34146636 -0.97401280
[138,] -11.03785002 -9.34146636
[139,] -10.31514134 -11.03785002
[140,] -3.95267213 -10.31514134
[141,] -10.31514134 -3.95267213
[142,] -3.64853293 -10.31514134
[143,] -6.54597439 -3.64853293
[144,] 1.63976627 -6.54597439
[145,] -1.94886225 1.63976627
[146,] -10.31514134 -1.94886225
[147,] -3.38981766 -10.31514134
[148,] -9.74195746 -3.38981766
[149,] -10.45909679 -9.74195746
[150,] -5.56170040 -10.45909679
[151,] 26.45378438 -5.56170040
[152,] -5.31514134 26.45378438
[153,] -5.22933725 -5.31514134
[154,] -8.47217620 -5.22933725
[155,] -10.31514134 -8.47217620
[156,] -10.31514134 -10.31514134
[157,] -9.01384468 -10.31514134
[158,] -7.47713244 -9.01384468
[159,] -10.31514134 -7.47713244
[160,] -10.31514134 -10.31514134
[161,] -10.31514134 -10.31514134
[162,] -10.31514134 -10.31514134
[163,] -5.88424303 -10.31514134
[164,] 10.69152811 -5.88424303
[165,] -9.56259654 10.69152811
[166,] -10.31514134 -9.56259654
[167,] 7.26633980 -10.31514134
[168,] 20.17626456 7.26633980
[169,] -10.31514134 20.17626456
[170,] -10.31514134 -10.31514134
[171,] -7.70497577 -10.31514134
[172,] -10.31514134 -7.70497577
[173,] -2.67298745 -10.31514134
[174,] -10.31514134 -2.67298745
[175,] -3.03335174 -10.31514134
[176,] 100.57885702 -3.03335174
[177,] -18.47678277 100.57885702
[178,] -10.31514134 -18.47678277
[179,] -10.31514134 -10.31514134
[180,] 0.97531480 -10.31514134
[181,] -12.48447433 0.97531480
[182,] 2.89717434 -12.48447433
[183,] -6.73398589 2.89717434
[184,] -13.65360363 -6.73398589
[185,] -10.37304679 -13.65360363
[186,] -8.19211529 -10.37304679
[187,] -5.41312073 -8.19211529
[188,] -7.69754265 -5.41312073
[189,] -10.31514134 -7.69754265
[190,] -4.47349256 -10.31514134
[191,] -10.31514134 -4.47349256
[192,] -4.27346403 -10.31514134
[193,] -10.31514134 -4.27346403
[194,] -7.50290509 -10.31514134
[195,] -18.09703459 -7.50290509
[196,] 32.15274937 -18.09703459
[197,] -18.09703459 32.15274937
[198,] -18.09703459 -18.09703459
[199,] -6.16823853 -18.09703459
[200,] -18.09703459 -6.16823853
[201,] -18.09703459 -18.09703459
[202,] -18.09703459 -18.09703459
[203,] -4.81260995 -18.09703459
[204,] -15.03768096 -4.81260995
[205,] -13.37181194 -15.03768096
[206,] -17.60771325 -13.37181194
[207,] -18.09703459 -17.60771325
[208,] 4.81669519 -18.09703459
[209,] -16.49226068 4.81669519
[210,] -10.31514134 -16.49226068
[211,] -10.31514134 -10.31514134
[212,] 0.01003430 -10.31514134
[213,] -3.10905887 0.01003430
[214,] -8.57755298 -3.10905887
[215,] -8.31514134 -8.57755298
[216,] -4.42045843 -8.31514134
[217,] -5.38778663 -4.42045843
[218,] -8.31514134 -5.38778663
[219,] 8.92709417 -8.31514134
[220,] -8.78162627 8.92709417
[221,] -10.31514134 -8.78162627
[222,] -9.32571542 -10.31514134
[223,] 7.22083674 -9.32571542
[224,] -9.44061324 7.22083674
[225,] 10.05886372 -9.44061324
[226,] -1.63630659 10.05886372
[227,] -10.31514134 -1.63630659
[228,] -10.40723154 -10.31514134
[229,] -10.31514134 -10.40723154
[230,] -5.74237066 -10.31514134
[231,] -5.20071437 -5.74237066
[232,] -10.31514134 -5.20071437
[233,] 0.51809038 -10.31514134
[234,] -10.31514134 0.51809038
[235,] -10.33112961 -10.31514134
[236,] -10.31514134 -10.33112961
[237,] -4.50614251 -10.31514134
[238,] -10.31514134 -4.50614251
[239,] -10.31514134 -10.31514134
[240,] -10.31514134 -10.31514134
[241,] -7.96193365 -10.31514134
[242,] -10.31514134 -7.96193365
[243,] -3.22744286 -10.31514134
[244,] -23.43142427 -3.22744286
[245,] -18.09703459 -23.43142427
[246,] -18.09703459 -18.09703459
[247,] 3.73694386 -18.09703459
[248,] -18.09703459 3.73694386
[249,] -18.09703459 -18.09703459
[250,] -0.59376125 -18.09703459
[251,] -18.09703459 -0.59376125
[252,] -11.03093416 -18.09703459
[253,] 21.49413058 -11.03093416
[254,] 29.20414036 21.49413058
[255,] -10.31514134 29.20414036
[256,] -3.83670640 -10.31514134
[257,] 16.85705837 -3.83670640
[258,] -18.09703459 16.85705837
[259,] -10.31514134 -18.09703459
[260,] -18.09703459 -10.31514134
[261,] -18.09703459 -18.09703459
[262,] -1.38010953 -18.09703459
[263,] -1.48914151 -1.38010953
[264,] -18.09703459 -1.48914151
[265,] -10.31514134 -18.09703459
[266,] -18.09703459 -10.31514134
[267,] -6.32390008 -18.09703459
[268,] -18.09703459 -6.32390008
[269,] 12.79942900 -18.09703459
[270,] -18.09703459 12.79942900
[271,] -18.09703459 -18.09703459
[272,] -18.09703459 -18.09703459
[273,] 8.00822253 -18.09703459
[274,] -10.31514134 8.00822253
[275,] -0.21090294 -10.31514134
[276,] 10.10115724 -0.21090294
[277,] -10.31514134 10.10115724
[278,] 2.33733983 -10.31514134
[279,] -5.00171380 2.33733983
[280,] -18.09703459 -5.00171380
[281,] -10.31514134 -18.09703459
[282,] 15.32171073 -10.31514134
[283,] -10.31514134 15.32171073
[284,] -5.46313840 -10.31514134
[285,] -10.31514134 -5.46313840
[286,] -16.26130206 -10.31514134
[287,] -7.78753450 -16.26130206
[288,] -18.09703459 -7.78753450
[289,] 2.89197988 -18.09703459
[290,] -18.09703459 2.89197988
[291,] 0.31278681 -18.09703459
[292,] 1.45738906 0.31278681
[293,] -10.31514134 1.45738906
[294,] -0.30608555 -10.31514134
[295,] -10.78670366 -0.30608555
[296,] -10.31514134 -10.78670366
[297,] -18.09703459 -10.31514134
[298,] -10.31514134 -18.09703459
[299,] -1.31514134 -10.31514134
[300,] -14.27349532 -1.31514134
[301,] -10.38397362 -14.27349532
[302,] -9.31514134 -10.38397362
[303,] -18.09703459 -9.31514134
[304,] 1.49966789 -18.09703459
[305,] 1.67076892 1.49966789
[306,] -18.09703459 1.67076892
[307,] -0.97961233 -18.09703459
[308,] 6.09577074 -0.97961233
[309,] -18.09703459 6.09577074
[310,] -18.09703459 -18.09703459
[311,] 13.48091629 -18.09703459
[312,] -18.09703459 13.48091629
[313,] -5.15880739 -18.09703459
[314,] -2.34066629 -5.15880739
[315,] -8.26603956 -2.34066629
[316,] -10.30990653 -8.26603956
[317,] 7.64125712 -10.30990653
[318,] -5.32016903 7.64125712
[319,] -10.31514134 -5.32016903
[320,] -10.31514134 -10.31514134
[321,] -10.31514134 -10.31514134
[322,] -10.31514134 -10.31514134
[323,] -18.09703459 -10.31514134
[324,] -10.31514134 -18.09703459
[325,] -9.43144607 -10.31514134
[326,] -18.09703459 -9.43144607
[327,] -18.09703459 -18.09703459
[328,] -10.31514134 -18.09703459
[329,] -10.31514134 -10.31514134
[330,] -10.31514134 -10.31514134
[331,] -10.31514134 -10.31514134
[332,] -10.31514134 -10.31514134
[333,] -6.73965715 -10.31514134
[334,] -9.32206756 -6.73965715
[335,] -7.96835263 -9.32206756
[336,] 5.19254552 -7.96835263
[337,] -10.31514134 5.19254552
[338,] -9.49379898 -10.31514134
[339,] -0.73927378 -9.49379898
[340,] -0.44942515 -0.73927378
[341,] -10.31514134 -0.44942515
[342,] -10.31514134 -10.31514134
[343,] -10.31514134 -10.31514134
[344,] -0.58470889 -10.31514134
[345,] -10.31514134 -0.58470889
[346,] -7.56149740 -10.31514134
[347,] -10.31514134 -7.56149740
[348,] -10.31514134 -10.31514134
[349,] -8.71282289 -10.31514134
[350,] -8.45423914 -8.71282289
[351,] -3.72691705 -8.45423914
[352,] -3.13871882 -3.72691705
[353,] -10.34510624 -3.13871882
[354,] -10.31514134 -10.34510624
[355,] -10.31514134 -10.31514134
[356,] -2.30727135 -10.31514134
[357,] -10.31514134 -2.30727135
[358,] -10.31514134 -10.31514134
[359,] -10.31514134 -10.31514134
[360,] -6.32088098 -10.31514134
[361,] -10.31514134 -6.32088098
[362,] -10.31514134 -10.31514134
[363,] 0.41192554 -10.31514134
[364,] -2.00535477 0.41192554
[365,] -10.31514134 -2.00535477
[366,] -10.31514134 -10.31514134
[367,] -10.31514134 -10.31514134
[368,] -10.31514134 -10.31514134
[369,] -4.46291224 -10.31514134
[370,] -10.31514134 -4.46291224
[371,] -5.45102393 -10.31514134
[372,] -8.28930927 -5.45102393
[373,] -9.20424560 -8.28930927
[374,] -10.31514134 -9.20424560
[375,] -10.31514134 -10.31514134
[376,] 1.12568015 -10.31514134
[377,] -6.44360123 1.12568015
[378,] -11.96133042 -6.44360123
[379,] -10.31514134 -11.96133042
[380,] -10.31514134 -10.31514134
[381,] 1.87741455 -10.31514134
[382,] 5.24520344 1.87741455
[383,] -10.31514134 5.24520344
[384,] 0.80430365 -10.31514134
[385,] 3.89912046 0.80430365
[386,] -7.75353977 3.89912046
[387,] -8.67637758 -7.75353977
[388,] -7.61780605 -8.67637758
[389,] -4.81886852 -7.61780605
[390,] -2.87760581 -4.81886852
[391,] -6.35333807 -2.87760581
[392,] -5.64455646 -6.35333807
[393,] -9.71094995 -5.64455646
[394,] -32.56419855 -9.71094995
[395,] -4.73679136 -32.56419855
[396,] -8.32623824 -4.73679136
[397,] 198.96439569 -8.32623824
[398,] -6.84787233 198.96439569
[399,] 21.13550096 -6.84787233
[400,] -15.20208445 21.13550096
[401,] -7.23376193 -15.20208445
[402,] -7.16896614 -7.23376193
[403,] -7.17970092 -7.16896614
[404,] -7.27449379 -7.17970092
[405,] -7.90536097 -7.27449379
[406,] -4.31355229 -7.90536097
[407,] -3.09743566 -4.31355229
[408,] -3.05165594 -3.09743566
[409,] -13.33507583 -3.05165594
[410,] 37.65478271 -13.33507583
[411,] -36.75969481 37.65478271
[412,] -8.41862667 -36.75969481
[413,] -6.96764605 -8.41862667
[414,] -13.86174903 -6.96764605
[415,] -16.42638921 -13.86174903
[416,] 1.58643104 -16.42638921
[417,] -4.44421300 1.58643104
[418,] 4.98444413 -4.44421300
[419,] -14.02934965 4.98444413
[420,] -16.79629576 -14.02934965
[421,] 45.87218169 -16.79629576
[422,] 7.06954909 45.87218169
[423,] -23.51408533 7.06954909
[424,] -0.91895888 -23.51408533
[425,] 35.50470013 -0.91895888
[426,] -8.12298961 35.50470013
[427,] 23.20964267 -8.12298961
[428,] -31.48611750 23.20964267
[429,] -8.50065897 -31.48611750
[430,] -28.53467948 -8.50065897
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 264.06371393 -103.12144864
2 -88.99460634 264.06371393
3 -7.87939772 -88.99460634
4 -6.59766571 -7.87939772
5 -97.46474114 -6.59766571
6 20.80076843 -97.46474114
7 -30.32269121 20.80076843
8 11.45580057 -30.32269121
9 28.27686821 11.45580057
10 -86.17227154 28.27686821
11 47.94377848 -86.17227154
12 294.87478806 47.94377848
13 -42.97573934 294.87478806
14 123.05471294 -42.97573934
15 -3.44092154 123.05471294
16 0.20482440 -3.44092154
17 118.61198031 0.20482440
18 23.37394291 118.61198031
19 158.11423255 23.37394291
20 -18.60966894 158.11423255
21 -73.19279320 -18.60966894
22 -2.35389504 -73.19279320
23 181.23977737 -2.35389504
24 17.66401481 181.23977737
25 -38.36851597 17.66401481
26 68.86397589 -38.36851597
27 6.96004113 68.86397589
28 -14.64986439 6.96004113
29 55.36251778 -14.64986439
30 -16.10969250 55.36251778
31 -21.09658048 -16.10969250
32 -19.57853588 -21.09658048
33 -9.71658559 -19.57853588
34 -67.18461972 -9.71658559
35 93.08518901 -67.18461972
36 2.70226347 93.08518901
37 21.52616528 2.70226347
38 21.34913616 21.52616528
39 94.73262800 21.34913616
40 131.17992576 94.73262800
41 -35.01927573 131.17992576
42 47.00770795 -35.01927573
43 28.04002519 47.00770795
44 2.31571418 28.04002519
45 -22.26665686 2.31571418
46 -30.40942241 -22.26665686
47 73.37187276 -30.40942241
48 7.35736507 73.37187276
49 14.93777033 7.35736507
50 44.02296804 14.93777033
51 78.48259622 44.02296804
52 10.13471750 78.48259622
53 18.64292493 10.13471750
54 48.62286782 18.64292493
55 53.98405474 48.62286782
56 -21.58741500 53.98405474
57 1.99900318 -21.58741500
58 24.71884384 1.99900318
59 11.46018567 24.71884384
60 -12.37888920 11.46018567
61 13.73343469 -12.37888920
62 -17.59123856 13.73343469
63 13.94333651 -17.59123856
64 30.00512098 13.94333651
65 163.76242696 30.00512098
66 -4.71150432 163.76242696
67 5.67177477 -4.71150432
68 38.48143101 5.67177477
69 36.26781822 38.48143101
70 22.79499240 36.26781822
71 5.87853997 22.79499240
72 6.19496782 5.87853997
73 -8.63771584 6.19496782
74 -0.51603113 -8.63771584
75 -11.40111070 -0.51603113
76 0.53024646 -11.40111070
77 -6.73490421 0.53024646
78 32.87032938 -6.73490421
79 -22.16587073 32.87032938
80 11.35697634 -22.16587073
81 16.29885511 11.35697634
82 -0.49745998 16.29885511
83 5.40531575 -0.49745998
84 -7.98979095 5.40531575
85 19.09944882 -7.98979095
86 63.49893927 19.09944882
87 -10.41976570 63.49893927
88 12.74296952 -10.41976570
89 -9.41597178 12.74296952
90 6.88364933 -9.41597178
91 1.11209862 6.88364933
92 45.14099699 1.11209862
93 -39.20121921 45.14099699
94 28.24332426 -39.20121921
95 -4.06937189 28.24332426
96 -15.15289012 -4.06937189
97 -5.78770427 -15.15289012
98 -21.28905374 -5.78770427
99 -19.34521631 -21.28905374
100 2.61056046 -19.34521631
101 -3.52384460 2.61056046
102 15.62890398 -3.52384460
103 -17.93523504 15.62890398
104 6.91523333 -17.93523504
105 12.18442828 6.91523333
106 13.30804579 12.18442828
107 -21.68498505 13.30804579
108 -3.00941175 -21.68498505
109 -11.17234032 -3.00941175
110 5.04363424 -11.17234032
111 -0.17942150 5.04363424
112 0.64444480 -0.17942150
113 -20.87354555 0.64444480
114 -18.65308885 -20.87354555
115 17.80530380 -18.65308885
116 -10.19132013 17.80530380
117 -11.71746846 -10.19132013
118 -22.77001316 -11.71746846
119 6.14942522 -22.77001316
120 14.21814203 6.14942522
121 -4.41236126 14.21814203
122 85.65578864 -4.41236126
123 -12.95156869 85.65578864
124 -15.55559455 -12.95156869
125 -16.02626991 -15.55559455
126 -10.31514134 -16.02626991
127 -0.82737576 -10.31514134
128 11.12395059 -0.82737576
129 -10.31514134 11.12395059
130 -3.47279309 -10.31514134
131 7.68825694 -3.47279309
132 -10.31514134 7.68825694
133 9.41883179 -10.31514134
134 -9.71106899 9.41883179
135 -7.92785981 -9.71106899
136 -0.97401280 -7.92785981
137 -9.34146636 -0.97401280
138 -11.03785002 -9.34146636
139 -10.31514134 -11.03785002
140 -3.95267213 -10.31514134
141 -10.31514134 -3.95267213
142 -3.64853293 -10.31514134
143 -6.54597439 -3.64853293
144 1.63976627 -6.54597439
145 -1.94886225 1.63976627
146 -10.31514134 -1.94886225
147 -3.38981766 -10.31514134
148 -9.74195746 -3.38981766
149 -10.45909679 -9.74195746
150 -5.56170040 -10.45909679
151 26.45378438 -5.56170040
152 -5.31514134 26.45378438
153 -5.22933725 -5.31514134
154 -8.47217620 -5.22933725
155 -10.31514134 -8.47217620
156 -10.31514134 -10.31514134
157 -9.01384468 -10.31514134
158 -7.47713244 -9.01384468
159 -10.31514134 -7.47713244
160 -10.31514134 -10.31514134
161 -10.31514134 -10.31514134
162 -10.31514134 -10.31514134
163 -5.88424303 -10.31514134
164 10.69152811 -5.88424303
165 -9.56259654 10.69152811
166 -10.31514134 -9.56259654
167 7.26633980 -10.31514134
168 20.17626456 7.26633980
169 -10.31514134 20.17626456
170 -10.31514134 -10.31514134
171 -7.70497577 -10.31514134
172 -10.31514134 -7.70497577
173 -2.67298745 -10.31514134
174 -10.31514134 -2.67298745
175 -3.03335174 -10.31514134
176 100.57885702 -3.03335174
177 -18.47678277 100.57885702
178 -10.31514134 -18.47678277
179 -10.31514134 -10.31514134
180 0.97531480 -10.31514134
181 -12.48447433 0.97531480
182 2.89717434 -12.48447433
183 -6.73398589 2.89717434
184 -13.65360363 -6.73398589
185 -10.37304679 -13.65360363
186 -8.19211529 -10.37304679
187 -5.41312073 -8.19211529
188 -7.69754265 -5.41312073
189 -10.31514134 -7.69754265
190 -4.47349256 -10.31514134
191 -10.31514134 -4.47349256
192 -4.27346403 -10.31514134
193 -10.31514134 -4.27346403
194 -7.50290509 -10.31514134
195 -18.09703459 -7.50290509
196 32.15274937 -18.09703459
197 -18.09703459 32.15274937
198 -18.09703459 -18.09703459
199 -6.16823853 -18.09703459
200 -18.09703459 -6.16823853
201 -18.09703459 -18.09703459
202 -18.09703459 -18.09703459
203 -4.81260995 -18.09703459
204 -15.03768096 -4.81260995
205 -13.37181194 -15.03768096
206 -17.60771325 -13.37181194
207 -18.09703459 -17.60771325
208 4.81669519 -18.09703459
209 -16.49226068 4.81669519
210 -10.31514134 -16.49226068
211 -10.31514134 -10.31514134
212 0.01003430 -10.31514134
213 -3.10905887 0.01003430
214 -8.57755298 -3.10905887
215 -8.31514134 -8.57755298
216 -4.42045843 -8.31514134
217 -5.38778663 -4.42045843
218 -8.31514134 -5.38778663
219 8.92709417 -8.31514134
220 -8.78162627 8.92709417
221 -10.31514134 -8.78162627
222 -9.32571542 -10.31514134
223 7.22083674 -9.32571542
224 -9.44061324 7.22083674
225 10.05886372 -9.44061324
226 -1.63630659 10.05886372
227 -10.31514134 -1.63630659
228 -10.40723154 -10.31514134
229 -10.31514134 -10.40723154
230 -5.74237066 -10.31514134
231 -5.20071437 -5.74237066
232 -10.31514134 -5.20071437
233 0.51809038 -10.31514134
234 -10.31514134 0.51809038
235 -10.33112961 -10.31514134
236 -10.31514134 -10.33112961
237 -4.50614251 -10.31514134
238 -10.31514134 -4.50614251
239 -10.31514134 -10.31514134
240 -10.31514134 -10.31514134
241 -7.96193365 -10.31514134
242 -10.31514134 -7.96193365
243 -3.22744286 -10.31514134
244 -23.43142427 -3.22744286
245 -18.09703459 -23.43142427
246 -18.09703459 -18.09703459
247 3.73694386 -18.09703459
248 -18.09703459 3.73694386
249 -18.09703459 -18.09703459
250 -0.59376125 -18.09703459
251 -18.09703459 -0.59376125
252 -11.03093416 -18.09703459
253 21.49413058 -11.03093416
254 29.20414036 21.49413058
255 -10.31514134 29.20414036
256 -3.83670640 -10.31514134
257 16.85705837 -3.83670640
258 -18.09703459 16.85705837
259 -10.31514134 -18.09703459
260 -18.09703459 -10.31514134
261 -18.09703459 -18.09703459
262 -1.38010953 -18.09703459
263 -1.48914151 -1.38010953
264 -18.09703459 -1.48914151
265 -10.31514134 -18.09703459
266 -18.09703459 -10.31514134
267 -6.32390008 -18.09703459
268 -18.09703459 -6.32390008
269 12.79942900 -18.09703459
270 -18.09703459 12.79942900
271 -18.09703459 -18.09703459
272 -18.09703459 -18.09703459
273 8.00822253 -18.09703459
274 -10.31514134 8.00822253
275 -0.21090294 -10.31514134
276 10.10115724 -0.21090294
277 -10.31514134 10.10115724
278 2.33733983 -10.31514134
279 -5.00171380 2.33733983
280 -18.09703459 -5.00171380
281 -10.31514134 -18.09703459
282 15.32171073 -10.31514134
283 -10.31514134 15.32171073
284 -5.46313840 -10.31514134
285 -10.31514134 -5.46313840
286 -16.26130206 -10.31514134
287 -7.78753450 -16.26130206
288 -18.09703459 -7.78753450
289 2.89197988 -18.09703459
290 -18.09703459 2.89197988
291 0.31278681 -18.09703459
292 1.45738906 0.31278681
293 -10.31514134 1.45738906
294 -0.30608555 -10.31514134
295 -10.78670366 -0.30608555
296 -10.31514134 -10.78670366
297 -18.09703459 -10.31514134
298 -10.31514134 -18.09703459
299 -1.31514134 -10.31514134
300 -14.27349532 -1.31514134
301 -10.38397362 -14.27349532
302 -9.31514134 -10.38397362
303 -18.09703459 -9.31514134
304 1.49966789 -18.09703459
305 1.67076892 1.49966789
306 -18.09703459 1.67076892
307 -0.97961233 -18.09703459
308 6.09577074 -0.97961233
309 -18.09703459 6.09577074
310 -18.09703459 -18.09703459
311 13.48091629 -18.09703459
312 -18.09703459 13.48091629
313 -5.15880739 -18.09703459
314 -2.34066629 -5.15880739
315 -8.26603956 -2.34066629
316 -10.30990653 -8.26603956
317 7.64125712 -10.30990653
318 -5.32016903 7.64125712
319 -10.31514134 -5.32016903
320 -10.31514134 -10.31514134
321 -10.31514134 -10.31514134
322 -10.31514134 -10.31514134
323 -18.09703459 -10.31514134
324 -10.31514134 -18.09703459
325 -9.43144607 -10.31514134
326 -18.09703459 -9.43144607
327 -18.09703459 -18.09703459
328 -10.31514134 -18.09703459
329 -10.31514134 -10.31514134
330 -10.31514134 -10.31514134
331 -10.31514134 -10.31514134
332 -10.31514134 -10.31514134
333 -6.73965715 -10.31514134
334 -9.32206756 -6.73965715
335 -7.96835263 -9.32206756
336 5.19254552 -7.96835263
337 -10.31514134 5.19254552
338 -9.49379898 -10.31514134
339 -0.73927378 -9.49379898
340 -0.44942515 -0.73927378
341 -10.31514134 -0.44942515
342 -10.31514134 -10.31514134
343 -10.31514134 -10.31514134
344 -0.58470889 -10.31514134
345 -10.31514134 -0.58470889
346 -7.56149740 -10.31514134
347 -10.31514134 -7.56149740
348 -10.31514134 -10.31514134
349 -8.71282289 -10.31514134
350 -8.45423914 -8.71282289
351 -3.72691705 -8.45423914
352 -3.13871882 -3.72691705
353 -10.34510624 -3.13871882
354 -10.31514134 -10.34510624
355 -10.31514134 -10.31514134
356 -2.30727135 -10.31514134
357 -10.31514134 -2.30727135
358 -10.31514134 -10.31514134
359 -10.31514134 -10.31514134
360 -6.32088098 -10.31514134
361 -10.31514134 -6.32088098
362 -10.31514134 -10.31514134
363 0.41192554 -10.31514134
364 -2.00535477 0.41192554
365 -10.31514134 -2.00535477
366 -10.31514134 -10.31514134
367 -10.31514134 -10.31514134
368 -10.31514134 -10.31514134
369 -4.46291224 -10.31514134
370 -10.31514134 -4.46291224
371 -5.45102393 -10.31514134
372 -8.28930927 -5.45102393
373 -9.20424560 -8.28930927
374 -10.31514134 -9.20424560
375 -10.31514134 -10.31514134
376 1.12568015 -10.31514134
377 -6.44360123 1.12568015
378 -11.96133042 -6.44360123
379 -10.31514134 -11.96133042
380 -10.31514134 -10.31514134
381 1.87741455 -10.31514134
382 5.24520344 1.87741455
383 -10.31514134 5.24520344
384 0.80430365 -10.31514134
385 3.89912046 0.80430365
386 -7.75353977 3.89912046
387 -8.67637758 -7.75353977
388 -7.61780605 -8.67637758
389 -4.81886852 -7.61780605
390 -2.87760581 -4.81886852
391 -6.35333807 -2.87760581
392 -5.64455646 -6.35333807
393 -9.71094995 -5.64455646
394 -32.56419855 -9.71094995
395 -4.73679136 -32.56419855
396 -8.32623824 -4.73679136
397 198.96439569 -8.32623824
398 -6.84787233 198.96439569
399 21.13550096 -6.84787233
400 -15.20208445 21.13550096
401 -7.23376193 -15.20208445
402 -7.16896614 -7.23376193
403 -7.17970092 -7.16896614
404 -7.27449379 -7.17970092
405 -7.90536097 -7.27449379
406 -4.31355229 -7.90536097
407 -3.09743566 -4.31355229
408 -3.05165594 -3.09743566
409 -13.33507583 -3.05165594
410 37.65478271 -13.33507583
411 -36.75969481 37.65478271
412 -8.41862667 -36.75969481
413 -6.96764605 -8.41862667
414 -13.86174903 -6.96764605
415 -16.42638921 -13.86174903
416 1.58643104 -16.42638921
417 -4.44421300 1.58643104
418 4.98444413 -4.44421300
419 -14.02934965 4.98444413
420 -16.79629576 -14.02934965
421 45.87218169 -16.79629576
422 7.06954909 45.87218169
423 -23.51408533 7.06954909
424 -0.91895888 -23.51408533
425 35.50470013 -0.91895888
426 -8.12298961 35.50470013
427 23.20964267 -8.12298961
428 -31.48611750 23.20964267
429 -8.50065897 -31.48611750
430 -28.53467948 -8.50065897
> 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/7s7kv1291216940.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/8s7kv1291216940.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/9s7kv1291216940.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/103hjy1291216940.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/116zim1291216940.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/12siza1291216940.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/13oswj1291216940.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/14radp1291216940.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/15usbu1291216940.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/16929l1291216940.tab")
+ }
>
> try(system("convert tmp/1wgnn1291216940.ps tmp/1wgnn1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/2wgnn1291216940.ps tmp/2wgnn1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/3774q1291216940.ps tmp/3774q1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/4774q1291216940.ps tmp/4774q1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/5774q1291216940.ps tmp/5774q1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/6zgls1291216940.ps tmp/6zgls1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/7s7kv1291216940.ps tmp/7s7kv1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/8s7kv1291216940.ps tmp/8s7kv1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/9s7kv1291216940.ps tmp/9s7kv1291216940.png",intern=TRUE))
character(0)
> try(system("convert tmp/103hjy1291216940.ps tmp/103hjy1291216940.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.760 2.169 26.446