R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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.
Natural language support but running in an English locale
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(1081
+ ,1
+ ,213118
+ ,162556
+ ,309
+ ,1
+ ,81767
+ ,29790
+ ,458
+ ,1
+ ,153198
+ ,87550
+ ,588
+ ,0
+ ,-26007
+ ,84738
+ ,302
+ ,1
+ ,126942
+ ,54660
+ ,156
+ ,1
+ ,157214
+ ,42634
+ ,481
+ ,0
+ ,129352
+ ,40949
+ ,353
+ ,1
+ ,234817
+ ,45187
+ ,452
+ ,1
+ ,60448
+ ,37704
+ ,109
+ ,1
+ ,47818
+ ,16275
+ ,115
+ ,0
+ ,245546
+ ,25830
+ ,110
+ ,0
+ ,48020
+ ,12679
+ ,239
+ ,1
+ ,-1710
+ ,18014
+ ,247
+ ,0
+ ,32648
+ ,43556
+ ,505
+ ,1
+ ,95350
+ ,24811
+ ,159
+ ,0
+ ,151352
+ ,6575
+ ,109
+ ,0
+ ,288170
+ ,7123
+ ,519
+ ,1
+ ,114337
+ ,21950
+ ,248
+ ,1
+ ,37884
+ ,37597
+ ,373
+ ,0
+ ,122844
+ ,17821
+ ,119
+ ,1
+ ,82340
+ ,12988
+ ,84
+ ,1
+ ,79801
+ ,22330
+ ,102
+ ,0
+ ,165548
+ ,13326
+ ,295
+ ,0
+ ,116384
+ ,16189
+ ,105
+ ,0
+ ,134028
+ ,7146
+ ,64
+ ,0
+ ,63838
+ ,15824
+ ,282
+ ,1
+ ,74996
+ ,27664
+ ,182
+ ,0
+ ,31080
+ ,11920
+ ,37
+ ,0
+ ,32168
+ ,8568
+ ,361
+ ,0
+ ,49857
+ ,14416
+ ,28
+ ,1
+ ,87161
+ ,3369
+ ,85
+ ,1
+ ,106113
+ ,11819
+ ,45
+ ,1
+ ,80570
+ ,6984
+ ,49
+ ,1
+ ,102129
+ ,4519
+ ,22
+ ,0
+ ,301670
+ ,2220
+ ,155
+ ,0
+ ,102313
+ ,18562
+ ,91
+ ,0
+ ,88577
+ ,10327
+ ,81
+ ,1
+ ,112477
+ ,5336
+ ,79
+ ,1
+ ,191778
+ ,2365
+ ,145
+ ,0
+ ,79804
+ ,4069
+ ,855
+ ,0
+ ,128294
+ ,8636
+ ,61
+ ,0
+ ,96448
+ ,13718
+ ,226
+ ,0
+ ,93811
+ ,4525
+ ,105
+ ,0
+ ,117520
+ ,6869
+ ,62
+ ,0
+ ,69159
+ ,4628
+ ,25
+ ,1
+ ,101792
+ ,3689
+ ,217
+ ,1
+ ,210568
+ ,4891
+ ,322
+ ,1
+ ,136996
+ ,7489
+ ,84
+ ,0
+ ,121920
+ ,4901
+ ,33
+ ,0
+ ,76403
+ ,2284
+ ,108
+ ,1
+ ,108094
+ ,3160
+ ,150
+ ,1
+ ,134759
+ ,4150
+ ,115
+ ,1
+ ,188873
+ ,7285
+ ,162
+ ,1
+ ,146216
+ ,1134
+ ,158
+ ,1
+ ,156608
+ ,4658
+ ,97
+ ,0
+ ,61348
+ ,2384
+ ,9
+ ,0
+ ,50350
+ ,3748
+ ,66
+ ,0
+ ,87720
+ ,5371
+ ,107
+ ,0
+ ,99489
+ ,1285
+ ,101
+ ,1
+ ,87419
+ ,9327
+ ,47
+ ,1
+ ,94355
+ ,5565
+ ,38
+ ,0
+ ,60326
+ ,1528
+ ,34
+ ,1
+ ,94670
+ ,3122
+ ,87
+ ,1
+ ,82425
+ ,7561
+ ,79
+ ,0
+ ,59017
+ ,2675
+ ,947
+ ,0
+ ,90829
+ ,13253
+ ,74
+ ,0
+ ,80791
+ ,880
+ ,53
+ ,1
+ ,100423
+ ,2053
+ ,94
+ ,0
+ ,131116
+ ,1424
+ ,63
+ ,1
+ ,100269
+ ,4036
+ ,58
+ ,1
+ ,27330
+ ,3045
+ ,49
+ ,0
+ ,39039
+ ,5119
+ ,34
+ ,0
+ ,106885
+ ,1431
+ ,11
+ ,0
+ ,79285
+ ,554
+ ,35
+ ,0
+ ,118881
+ ,1975
+ ,20
+ ,1
+ ,77623
+ ,1765
+ ,47
+ ,0
+ ,114768
+ ,1012
+ ,43
+ ,0
+ ,74015
+ ,810
+ ,117
+ ,0
+ ,69465
+ ,1280
+ ,171
+ ,1
+ ,117869
+ ,666
+ ,26
+ ,0
+ ,60982
+ ,1380
+ ,75
+ ,1
+ ,90131
+ ,4677
+ ,59
+ ,0
+ ,138971
+ ,876
+ ,18
+ ,0
+ ,39625
+ ,814
+ ,15
+ ,0
+ ,102725
+ ,514
+ ,72
+ ,1
+ ,64239
+ ,5692
+ ,86
+ ,0
+ ,90262
+ ,3642
+ ,14
+ ,0
+ ,103960
+ ,540
+ ,64
+ ,0
+ ,106611
+ ,2099
+ ,11
+ ,0
+ ,103345
+ ,567
+ ,52
+ ,0
+ ,95551
+ ,2001
+ ,41
+ ,1
+ ,82903
+ ,2949
+ ,99
+ ,0
+ ,63593
+ ,2253
+ ,75
+ ,1
+ ,126910
+ ,6533
+ ,45
+ ,0
+ ,37527
+ ,1889
+ ,43
+ ,1
+ ,60247
+ ,3055
+ ,8
+ ,0
+ ,112995
+ ,272
+ ,198
+ ,1
+ ,70184
+ ,1414
+ ,22
+ ,0
+ ,130140
+ ,2564
+ ,11
+ ,1
+ ,73221
+ ,1383
+ ,33
+ ,1
+ ,76114
+ ,1261
+ ,23
+ ,0
+ ,90534
+ ,975
+ ,80
+ ,0
+ ,108479
+ ,3366
+ ,18
+ ,0
+ ,113761
+ ,576
+ ,40
+ ,0
+ ,68696
+ ,1686
+ ,23
+ ,0
+ ,71561
+ ,746
+ ,60
+ ,1
+ ,59831
+ ,3192
+ ,20
+ ,1
+ ,97890
+ ,2045
+ ,61
+ ,0
+ ,101481
+ ,5702
+ ,36
+ ,1
+ ,72954
+ ,1932
+ ,30
+ ,0
+ ,67939
+ ,936
+ ,47
+ ,1
+ ,48022
+ ,3437
+ ,71
+ ,0
+ ,86111
+ ,5131
+ ,14
+ ,1
+ ,74020
+ ,2397
+ ,9
+ ,1
+ ,57530
+ ,1389
+ ,39
+ ,0
+ ,56364
+ ,1503
+ ,26
+ ,0
+ ,84990
+ ,402
+ ,21
+ ,0
+ ,88590
+ ,2239
+ ,16
+ ,1
+ ,77200
+ ,2234
+ ,69
+ ,0
+ ,61262
+ ,837
+ ,92
+ ,0
+ ,110309
+ ,10579
+ ,14
+ ,0
+ ,67000
+ ,875
+ ,107
+ ,0
+ ,93099
+ ,1585
+ ,29
+ ,1
+ ,107577
+ ,1659
+ ,37
+ ,1
+ ,62920
+ ,2647
+ ,23
+ ,1
+ ,75832
+ ,3294
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,7
+ ,0
+ ,60793
+ ,94
+ ,28
+ ,0
+ ,57935
+ ,422
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,8
+ ,0
+ ,60630
+ ,34
+ ,63
+ ,0
+ ,55637
+ ,1558
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,3
+ ,0
+ ,60887
+ ,43
+ ,5
+ ,0
+ ,60720
+ ,645
+ ,9
+ ,0
+ ,60505
+ ,316
+ ,13
+ ,0
+ ,60945
+ ,115
+ ,2
+ ,0
+ ,60720
+ ,5
+ ,5
+ ,0
+ ,60720
+ ,897
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,14
+ ,0
+ ,58990
+ ,389
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,15
+ ,0
+ ,56750
+ ,1002
+ ,3
+ ,0
+ ,60894
+ ,36
+ ,15
+ ,0
+ ,63346
+ ,460
+ ,11
+ ,0
+ ,56535
+ ,309
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,6
+ ,0
+ ,60835
+ ,9
+ ,2
+ ,0
+ ,60720
+ ,271
+ ,1
+ ,0
+ ,61016
+ ,14
+ ,10
+ ,0
+ ,58650
+ ,520
+ ,73
+ ,0
+ ,60438
+ ,1766
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,11
+ ,0
+ ,58625
+ ,458
+ ,3
+ ,0
+ ,60938
+ ,20
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,2
+ ,0
+ ,61490
+ ,98
+ ,7
+ ,0
+ ,60845
+ ,405
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,27
+ ,0
+ ,60830
+ ,483
+ ,51
+ ,0
+ ,63261
+ ,454
+ ,3
+ ,0
+ ,60720
+ ,47
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,19
+ ,0
+ ,45689
+ ,757
+ ,393
+ ,0
+ ,60720
+ ,4655
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,4
+ ,0
+ ,61564
+ ,36
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,9
+ ,0
+ ,61938
+ ,203
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,10
+ ,0
+ ,60951
+ ,126
+ ,152
+ ,0
+ ,60720
+ ,400
+ ,1
+ ,1
+ ,60745
+ ,71
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,34
+ ,0
+ ,71642
+ ,972
+ ,10
+ ,1
+ ,71641
+ ,531
+ ,57
+ ,0
+ ,55792
+ ,2461
+ ,52
+ ,1
+ ,71873
+ ,378
+ ,5
+ ,1
+ ,62555
+ ,23
+ ,14
+ ,1
+ ,60370
+ ,638
+ ,29
+ ,1
+ ,64873
+ ,2300
+ ,5
+ ,0
+ ,62041
+ ,149
+ ,5
+ ,0
+ ,65745
+ ,226
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,4
+ ,0
+ ,59500
+ ,275
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,6
+ ,0
+ ,61630
+ ,141
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,2
+ ,0
+ ,60890
+ ,28
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,91
+ ,1
+ ,113521
+ ,4980
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,20
+ ,1
+ ,80045
+ ,472
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,27
+ ,1
+ ,50804
+ ,203
+ ,17
+ ,1
+ ,87390
+ ,496
+ ,2
+ ,1
+ ,61656
+ ,10
+ ,4
+ ,1
+ ,65688
+ ,63
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,32
+ ,1
+ ,48522
+ ,1136
+ ,31
+ ,1
+ ,60720
+ ,265
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,32
+ ,0
+ ,57640
+ ,267
+ ,20
+ ,0
+ ,61977
+ ,474
+ ,7
+ ,0
+ ,62620
+ ,534
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,8
+ ,0
+ ,60831
+ ,15
+ ,28
+ ,0
+ ,60646
+ ,397
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,29
+ ,0
+ ,56225
+ ,1866
+ ,4
+ ,0
+ ,60510
+ ,288
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,2
+ ,0
+ ,60698
+ ,3
+ ,21
+ ,0
+ ,60720
+ ,468
+ ,2
+ ,0
+ ,60805
+ ,20
+ ,26
+ ,0
+ ,61404
+ ,278
+ ,14
+ ,0
+ ,60720
+ ,61
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,4
+ ,0
+ ,65276
+ ,192
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,9
+ ,0
+ ,63915
+ ,317
+ ,10
+ ,0
+ ,60720
+ ,738
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,17
+ ,0
+ ,61686
+ ,368
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,0
+ ,60743
+ ,2
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,6
+ ,0
+ ,60349
+ ,53
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,3
+ ,0
+ ,61360
+ ,94
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,8
+ ,0
+ ,59818
+ ,24
+ ,4
+ ,0
+ ,72680
+ ,2332
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,11
+ ,0
+ ,61808
+ ,131
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,9
+ ,0
+ ,53110
+ ,206
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,2
+ ,0
+ ,64245
+ ,167
+ ,73
+ ,0
+ ,73007
+ ,622
+ ,94
+ ,0
+ ,82732
+ ,2328
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,8
+ ,0
+ ,54820
+ ,365
+ ,35
+ ,0
+ ,47705
+ ,364
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,12
+ ,0
+ ,72835
+ ,226
+ ,15
+ ,0
+ ,58856
+ ,307
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,11
+ ,0
+ ,77655
+ ,188
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,6
+ ,0
+ ,69817
+ ,138
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,12
+ ,0
+ ,60798
+ ,125
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,30
+ ,0
+ ,62452
+ ,282
+ ,33
+ ,0
+ ,64175
+ ,335
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,117
+ ,1
+ ,67440
+ ,1324
+ ,28
+ ,0
+ ,68136
+ ,176
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,72
+ ,0
+ ,56726
+ ,249
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,13
+ ,0
+ ,70811
+ ,333
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,6
+ ,1
+ ,60720
+ ,601
+ ,4
+ ,0
+ ,62045
+ ,30
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,62
+ ,0
+ ,54323
+ ,249
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,24
+ ,0
+ ,62841
+ ,165
+ ,21
+ ,0
+ ,81125
+ ,453
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,14
+ ,0
+ ,59506
+ ,53
+ ,21
+ ,1
+ ,59365
+ ,382
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,4
+ ,1
+ ,60798
+ ,30
+ ,2
+ ,0
+ ,58790
+ ,290
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,53
+ ,0
+ ,61808
+ ,366
+ ,9
+ ,0
+ ,60735
+ ,2
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,13
+ ,1
+ ,64016
+ ,209
+ ,22
+ ,0
+ ,54683
+ ,384
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,83
+ ,0
+ ,87192
+ ,365
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,8
+ ,1
+ ,64107
+ ,49
+ ,4
+ ,0
+ ,60761
+ ,3
+ ,14
+ ,0
+ ,65990
+ ,133
+ ,1
+ ,0
+ ,59988
+ ,32
+ ,17
+ ,0
+ ,61167
+ ,368
+ ,6
+ ,0
+ ,60719
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,2
+ ,0
+ ,60722
+ ,22
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,5
+ ,0
+ ,60379
+ ,96
+ ,2
+ ,0
+ ,60727
+ ,1
+ ,5
+ ,0
+ ,60720
+ ,314
+ ,78
+ ,0
+ ,60925
+ ,844
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,0
+ ,60896
+ ,26
+ ,13
+ ,0
+ ,59734
+ ,125
+ ,15
+ ,0
+ ,62969
+ ,304
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,48
+ ,0
+ ,60720
+ ,621
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,6
+ ,0
+ ,59118
+ ,119
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,17
+ ,0
+ ,60720
+ ,1595
+ ,14
+ ,0
+ ,58598
+ ,312
+ ,10
+ ,0
+ ,61124
+ ,60
+ ,12
+ ,0
+ ,59595
+ ,587
+ ,2
+ ,0
+ ,62065
+ ,135
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,52
+ ,0
+ ,78780
+ ,514
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,4
+ ,0
+ ,60722
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,24
+ ,1
+ ,61600
+ ,1763
+ ,11
+ ,0
+ ,59635
+ ,180
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,21
+ ,0
+ ,60720
+ ,218
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,40
+ ,0
+ ,59781
+ ,448
+ ,9
+ ,0
+ ,76644
+ ,227
+ ,1
+ ,0
+ ,64820
+ ,174
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,24
+ ,0
+ ,56178
+ ,121
+ ,11
+ ,0
+ ,60436
+ ,607
+ ,14
+ ,0
+ ,60720
+ ,2212
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,60
+ ,0
+ ,73433
+ ,530
+ ,80
+ ,0
+ ,41477
+ ,571
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,16
+ ,0
+ ,62700
+ ,78
+ ,40
+ ,0
+ ,67804
+ ,2489
+ ,6
+ ,0
+ ,59661
+ ,131
+ ,8
+ ,0
+ ,58620
+ ,923
+ ,3
+ ,0
+ ,60398
+ ,72
+ ,16
+ ,0
+ ,58580
+ ,572
+ ,10
+ ,0
+ ,62710
+ ,397
+ ,8
+ ,0
+ ,59325
+ ,450
+ ,7
+ ,0
+ ,60950
+ ,622
+ ,8
+ ,0
+ ,68060
+ ,694
+ ,12
+ ,1
+ ,83620
+ ,3425
+ ,13
+ ,0
+ ,58456
+ ,562
+ ,42
+ ,0
+ ,52811
+ ,4917
+ ,118
+ ,1
+ ,121173
+ ,1442
+ ,9
+ ,0
+ ,63870
+ ,529
+ ,138
+ ,1
+ ,21001
+ ,2126
+ ,5
+ ,0
+ ,70415
+ ,1061
+ ,9
+ ,0
+ ,64230
+ ,776
+ ,8
+ ,0
+ ,59190
+ ,611
+ ,25
+ ,1
+ ,69351
+ ,1526
+ ,7
+ ,0
+ ,64270
+ ,592
+ ,13
+ ,0
+ ,70694
+ ,1182
+ ,16
+ ,0
+ ,68005
+ ,621
+ ,11
+ ,0
+ ,58930
+ ,989
+ ,11
+ ,0
+ ,58320
+ ,438
+ ,3
+ ,0
+ ,69980
+ ,726
+ ,61
+ ,0
+ ,69863
+ ,1303
+ ,29
+ ,1
+ ,63255
+ ,7419
+ ,17
+ ,1
+ ,57320
+ ,1164
+ ,33
+ ,1
+ ,75230
+ ,3310
+ ,15
+ ,0
+ ,79420
+ ,1920
+ ,3
+ ,0
+ ,73490
+ ,965
+ ,66
+ ,0
+ ,35250
+ ,3256
+ ,17
+ ,1
+ ,62285
+ ,1135
+ ,26
+ ,0
+ ,69206
+ ,1270
+ ,3
+ ,0
+ ,65920
+ ,661
+ ,2
+ ,0
+ ,69770
+ ,1013
+ ,67
+ ,0
+ ,72683
+ ,2844
+ ,70
+ ,1
+ ,-14545
+ ,11528
+ ,26
+ ,0
+ ,55830
+ ,6526
+ ,24
+ ,0
+ ,55174
+ ,2264
+ ,97
+ ,1
+ ,67038
+ ,5109
+ ,30
+ ,0
+ ,51252
+ ,3999
+ ,223
+ ,0
+ ,157278
+ ,35624
+ ,48
+ ,0
+ ,79510
+ ,9252
+ ,90
+ ,0
+ ,77440
+ ,15236
+ ,180
+ ,0
+ ,27284
+ ,18073)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Trades'
+ ,'Group'
+ ,'Dividends'
+ ,'Costs
')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Trades','Group','Dividends','Costs
'),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 = '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
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
Trades Group Dividends Costs\r
1 1081 1 213118 162556
2 309 1 81767 29790
3 458 1 153198 87550
4 588 0 -26007 84738
5 302 1 126942 54660
6 156 1 157214 42634
7 481 0 129352 40949
8 353 1 234817 45187
9 452 1 60448 37704
10 109 1 47818 16275
11 115 0 245546 25830
12 110 0 48020 12679
13 239 1 -1710 18014
14 247 0 32648 43556
15 505 1 95350 24811
16 159 0 151352 6575
17 109 0 288170 7123
18 519 1 114337 21950
19 248 1 37884 37597
20 373 0 122844 17821
21 119 1 82340 12988
22 84 1 79801 22330
23 102 0 165548 13326
24 295 0 116384 16189
25 105 0 134028 7146
26 64 0 63838 15824
27 282 1 74996 27664
28 182 0 31080 11920
29 37 0 32168 8568
30 361 0 49857 14416
31 28 1 87161 3369
32 85 1 106113 11819
33 45 1 80570 6984
34 49 1 102129 4519
35 22 0 301670 2220
36 155 0 102313 18562
37 91 0 88577 10327
38 81 1 112477 5336
39 79 1 191778 2365
40 145 0 79804 4069
41 855 0 128294 8636
42 61 0 96448 13718
43 226 0 93811 4525
44 105 0 117520 6869
45 62 0 69159 4628
46 25 1 101792 3689
47 217 1 210568 4891
48 322 1 136996 7489
49 84 0 121920 4901
50 33 0 76403 2284
51 108 1 108094 3160
52 150 1 134759 4150
53 115 1 188873 7285
54 162 1 146216 1134
55 158 1 156608 4658
56 97 0 61348 2384
57 9 0 50350 3748
58 66 0 87720 5371
59 107 0 99489 1285
60 101 1 87419 9327
61 47 1 94355 5565
62 38 0 60326 1528
63 34 1 94670 3122
64 87 1 82425 7561
65 79 0 59017 2675
66 947 0 90829 13253
67 74 0 80791 880
68 53 1 100423 2053
69 94 0 131116 1424
70 63 1 100269 4036
71 58 1 27330 3045
72 49 0 39039 5119
73 34 0 106885 1431
74 11 0 79285 554
75 35 0 118881 1975
76 20 1 77623 1765
77 47 0 114768 1012
78 43 0 74015 810
79 117 0 69465 1280
80 171 1 117869 666
81 26 0 60982 1380
82 75 1 90131 4677
83 59 0 138971 876
84 18 0 39625 814
85 15 0 102725 514
86 72 1 64239 5692
87 86 0 90262 3642
88 14 0 103960 540
89 64 0 106611 2099
90 11 0 103345 567
91 52 0 95551 2001
92 41 1 82903 2949
93 99 0 63593 2253
94 75 1 126910 6533
95 45 0 37527 1889
96 43 1 60247 3055
97 8 0 112995 272
98 198 1 70184 1414
99 22 0 130140 2564
100 11 1 73221 1383
101 33 1 76114 1261
102 23 0 90534 975
103 80 0 108479 3366
104 18 0 113761 576
105 40 0 68696 1686
106 23 0 71561 746
107 60 1 59831 3192
108 20 1 97890 2045
109 61 0 101481 5702
110 36 1 72954 1932
111 30 0 67939 936
112 47 1 48022 3437
113 71 0 86111 5131
114 14 1 74020 2397
115 9 1 57530 1389
116 39 0 56364 1503
117 26 0 84990 402
118 21 0 88590 2239
119 16 1 77200 2234
120 69 0 61262 837
121 92 0 110309 10579
122 14 0 67000 875
123 107 0 93099 1585
124 29 1 107577 1659
125 37 1 62920 2647
126 23 1 75832 3294
127 0 0 60720 0
128 7 0 60793 94
129 28 0 57935 422
130 0 0 60720 0
131 8 0 60630 34
132 63 0 55637 1558
133 0 0 60720 0
134 3 0 60887 43
135 5 0 60720 645
136 9 0 60505 316
137 13 0 60945 115
138 2 0 60720 5
139 5 0 60720 897
140 0 0 60720 0
141 14 0 58990 389
142 0 0 60720 0
143 15 0 56750 1002
144 3 0 60894 36
145 15 0 63346 460
146 11 0 56535 309
147 0 0 60720 0
148 6 0 60835 9
149 2 0 60720 271
150 1 0 61016 14
151 10 0 58650 520
152 73 0 60438 1766
153 0 0 60720 0
154 11 0 58625 458
155 3 0 60938 20
156 0 0 60720 0
157 0 0 60720 0
158 2 0 61490 98
159 7 0 60845 405
160 0 0 60720 0
161 0 0 60720 0
162 0 0 60720 0
163 0 0 60720 0
164 27 0 60830 483
165 51 0 63261 454
166 3 0 60720 47
167 0 0 60720 0
168 19 0 45689 757
169 393 0 60720 4655
170 0 0 60720 0
171 0 0 60720 0
172 4 0 61564 36
173 0 0 60720 0
174 9 0 61938 203
175 0 0 60720 0
176 10 0 60951 126
177 152 0 60720 400
178 1 1 60745 71
179 0 0 60720 0
180 0 0 60720 0
181 34 0 71642 972
182 10 1 71641 531
183 57 0 55792 2461
184 52 1 71873 378
185 5 1 62555 23
186 14 1 60370 638
187 29 1 64873 2300
188 5 0 62041 149
189 5 0 65745 226
190 0 0 60720 0
191 4 0 59500 275
192 0 0 60720 0
193 6 0 61630 141
194 0 0 60720 0
195 2 0 60890 28
196 0 1 60720 0
197 91 1 113521 4980
198 0 1 60720 0
199 0 1 60720 0
200 20 1 80045 472
201 0 1 60720 0
202 0 1 60720 0
203 0 1 60720 0
204 27 1 50804 203
205 17 1 87390 496
206 2 1 61656 10
207 4 1 65688 63
208 0 1 60720 0
209 32 1 48522 1136
210 31 1 60720 265
211 0 0 60720 0
212 0 0 60720 0
213 32 0 57640 267
214 20 0 61977 474
215 7 0 62620 534
216 0 0 60720 0
217 8 0 60831 15
218 28 0 60646 397
219 0 0 60720 0
220 29 0 56225 1866
221 4 0 60510 288
222 0 0 60720 0
223 2 0 60698 3
224 21 0 60720 468
225 2 0 60805 20
226 26 0 61404 278
227 14 0 60720 61
228 0 0 60720 0
229 4 0 65276 192
230 0 0 60720 0
231 9 0 63915 317
232 10 0 60720 738
233 0 0 60720 0
234 17 0 61686 368
235 0 0 60720 0
236 1 0 60743 2
237 0 0 60720 0
238 6 0 60349 53
239 0 0 60720 0
240 0 0 60720 0
241 0 0 60720 0
242 3 0 61360 94
243 0 0 60720 0
244 8 0 59818 24
245 4 0 72680 2332
246 0 1 60720 0
247 0 1 60720 0
248 11 0 61808 131
249 0 1 60720 0
250 0 1 60720 0
251 9 0 53110 206
252 0 1 60720 0
253 2 0 64245 167
254 73 0 73007 622
255 94 0 82732 2328
256 0 0 60720 0
257 8 0 54820 365
258 35 0 47705 364
259 0 1 60720 0
260 0 0 60720 0
261 0 1 60720 0
262 0 1 60720 0
263 12 0 72835 226
264 15 0 58856 307
265 0 1 60720 0
266 0 0 60720 0
267 0 1 60720 0
268 11 0 77655 188
269 0 1 60720 0
270 6 0 69817 138
271 0 1 60720 0
272 0 1 60720 0
273 0 1 60720 0
274 12 0 60798 125
275 0 0 60720 0
276 30 0 62452 282
277 33 0 64175 335
278 0 0 60720 0
279 117 1 67440 1324
280 28 0 68136 176
281 0 1 60720 0
282 0 0 60720 0
283 72 0 56726 249
284 0 0 60720 0
285 13 0 70811 333
286 0 0 60720 0
287 6 1 60720 601
288 4 0 62045 30
289 0 1 60720 0
290 62 0 54323 249
291 0 1 60720 0
292 24 0 62841 165
293 21 0 81125 453
294 0 0 60720 0
295 14 0 59506 53
296 21 1 59365 382
297 0 0 60720 0
298 0 1 60720 0
299 0 0 60720 0
300 0 0 60720 0
301 4 1 60798 30
302 2 0 58790 290
303 0 0 60720 0
304 0 1 60720 0
305 53 0 61808 366
306 9 0 60735 2
307 0 1 60720 0
308 13 1 64016 209
309 22 0 54683 384
310 0 1 60720 0
311 0 1 60720 0
312 83 0 87192 365
313 0 1 60720 0
314 8 1 64107 49
315 4 0 60761 3
316 14 0 65990 133
317 1 0 59988 32
318 17 0 61167 368
319 6 0 60719 1
320 0 0 60720 0
321 0 0 60720 0
322 0 0 60720 0
323 0 0 60720 0
324 0 1 60720 0
325 0 0 60720 0
326 2 0 60722 22
327 0 1 60720 0
328 0 1 60720 0
329 0 0 60720 0
330 0 0 60720 0
331 0 0 60720 0
332 0 0 60720 0
333 0 0 60720 0
334 5 0 60379 96
335 2 0 60727 1
336 5 0 60720 314
337 78 0 60925 844
338 0 0 60720 0
339 1 0 60896 26
340 13 0 59734 125
341 15 0 62969 304
342 0 0 60720 0
343 0 0 60720 0
344 0 0 60720 0
345 48 0 60720 621
346 0 0 60720 0
347 6 0 59118 119
348 0 0 60720 0
349 0 0 60720 0
350 17 0 60720 1595
351 14 0 58598 312
352 10 0 61124 60
353 12 0 59595 587
354 2 0 62065 135
355 0 0 60720 0
356 0 0 60720 0
357 52 0 78780 514
358 0 0 60720 0
359 0 0 60720 0
360 0 0 60720 0
361 4 0 60722 1
362 0 0 60720 0
363 0 0 60720 0
364 24 1 61600 1763
365 11 0 59635 180
366 0 0 60720 0
367 0 0 60720 0
368 0 0 60720 0
369 0 0 60720 0
370 21 0 60720 218
371 0 0 60720 0
372 40 0 59781 448
373 9 0 76644 227
374 1 0 64820 174
375 0 0 60720 0
376 0 0 60720 0
377 24 0 56178 121
378 11 0 60436 607
379 14 0 60720 2212
380 0 0 60720 0
381 0 0 60720 0
382 60 0 73433 530
383 80 0 41477 571
384 0 0 60720 0
385 16 0 62700 78
386 40 0 67804 2489
387 6 0 59661 131
388 8 0 58620 923
389 3 0 60398 72
390 16 0 58580 572
391 10 0 62710 397
392 8 0 59325 450
393 7 0 60950 622
394 8 0 68060 694
395 12 1 83620 3425
396 13 0 58456 562
397 42 0 52811 4917
398 118 1 121173 1442
399 9 0 63870 529
400 138 1 21001 2126
401 5 0 70415 1061
402 9 0 64230 776
403 8 0 59190 611
404 25 1 69351 1526
405 7 0 64270 592
406 13 0 70694 1182
407 16 0 68005 621
408 11 0 58930 989
409 11 0 58320 438
410 3 0 69980 726
411 61 0 69863 1303
412 29 1 63255 7419
413 17 1 57320 1164
414 33 1 75230 3310
415 15 0 79420 1920
416 3 0 73490 965
417 66 0 35250 3256
418 17 1 62285 1135
419 26 0 69206 1270
420 3 0 65920 661
421 2 0 69770 1013
422 67 0 72683 2844
423 70 1 -14545 11528
424 26 0 55830 6526
425 24 0 55174 2264
426 97 1 67038 5109
427 30 0 51252 3999
428 223 0 157278 35624
429 48 0 79510 9252
430 90 0 77440 15236
431 180 0 27284 18073
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Dividends `Costs\r`
-7.0545694 1.2994476 0.0004474 0.0067798
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-198.3573 -20.1125 -17.4311 -0.1176 823.5639
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -7.0545694 8.9953847 -0.784 0.433332
Group 1.2994476 8.1004977 0.160 0.872629
Dividends 0.0004474 0.0001175 3.809 0.000160 ***
`Costs\r` 0.0067798 0.0003213 21.102 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 73.72 on 427 degrees of freedom
Multiple R-squared: 0.5779, Adjusted R-squared: 0.5749
F-statistic: 194.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,] 0.9822166 3.556675e-02 1.778338e-02
[2,] 0.9653323 6.933545e-02 3.466773e-02
[3,] 0.9979209 4.158152e-03 2.079076e-03
[4,] 0.9962346 7.530704e-03 3.765352e-03
[5,] 0.9984838 3.032313e-03 1.516157e-03
[6,] 0.9977332 4.533650e-03 2.266825e-03
[7,] 0.9964119 7.176230e-03 3.588115e-03
[8,] 0.9975970 4.805997e-03 2.402998e-03
[9,] 0.9999913 1.737335e-05 8.686674e-06
[10,] 0.9999868 2.638867e-05 1.319433e-05
[11,] 0.9999775 4.497556e-05 2.248778e-05
[12,] 0.9999999 1.454464e-07 7.272319e-08
[13,] 0.9999999 1.026389e-07 5.131944e-08
[14,] 1.0000000 1.565329e-08 7.826646e-09
[15,] 1.0000000 1.835985e-08 9.179925e-09
[16,] 1.0000000 3.268794e-09 1.634397e-09
[17,] 1.0000000 3.943350e-09 1.971675e-09
[18,] 1.0000000 3.731298e-09 1.865649e-09
[19,] 1.0000000 6.985717e-09 3.492859e-09
[20,] 1.0000000 3.042229e-09 1.521115e-09
[21,] 1.0000000 5.915186e-09 2.957593e-09
[22,] 1.0000000 1.172922e-08 5.864608e-09
[23,] 1.0000000 9.475281e-09 4.737640e-09
[24,] 1.0000000 9.037342e-10 4.518671e-10
[25,] 1.0000000 9.169261e-10 4.584630e-10
[26,] 1.0000000 1.060106e-09 5.300532e-10
[27,] 1.0000000 1.203669e-09 6.018347e-10
[28,] 1.0000000 1.854943e-09 9.274715e-10
[29,] 1.0000000 1.574582e-09 7.872909e-10
[30,] 1.0000000 2.188165e-09 1.094083e-09
[31,] 1.0000000 3.088649e-09 1.544324e-09
[32,] 1.0000000 5.604123e-09 2.802062e-09
[33,] 1.0000000 9.745248e-09 4.872624e-09
[34,] 1.0000000 1.524801e-08 7.624007e-09
[35,] 1.0000000 4.824515e-34 2.412258e-34
[36,] 1.0000000 9.792969e-35 4.896485e-35
[37,] 1.0000000 3.530437e-35 1.765218e-35
[38,] 1.0000000 8.403726e-35 4.201863e-35
[39,] 1.0000000 1.401715e-34 7.008575e-35
[40,] 1.0000000 2.186528e-34 1.093264e-34
[41,] 1.0000000 1.575171e-34 7.875853e-35
[42,] 1.0000000 1.277826e-36 6.389131e-37
[43,] 1.0000000 2.895761e-36 1.447880e-36
[44,] 1.0000000 4.347480e-36 2.173740e-36
[45,] 1.0000000 1.121558e-35 5.607792e-36
[46,] 1.0000000 2.212964e-35 1.106482e-35
[47,] 1.0000000 5.658321e-35 2.829161e-35
[48,] 1.0000000 5.044833e-35 2.522417e-35
[49,] 1.0000000 9.567070e-35 4.783535e-35
[50,] 1.0000000 1.920101e-34 9.600507e-35
[51,] 1.0000000 1.599445e-34 7.997225e-35
[52,] 1.0000000 3.102932e-34 1.551466e-34
[53,] 1.0000000 5.825096e-34 2.912548e-34
[54,] 1.0000000 1.315646e-33 6.578232e-34
[55,] 1.0000000 2.074044e-33 1.037022e-33
[56,] 1.0000000 3.936864e-33 1.968432e-33
[57,] 1.0000000 6.934487e-33 3.467244e-33
[58,] 1.0000000 1.603353e-32 8.016765e-33
[59,] 1.0000000 3.504505e-32 1.752252e-32
[60,] 1.0000000 5.115063e-106 2.557532e-106
[61,] 1.0000000 8.053715e-106 4.026857e-106
[62,] 1.0000000 3.508555e-105 1.754277e-105
[63,] 1.0000000 5.490621e-105 2.745311e-105
[64,] 1.0000000 2.431742e-104 1.215871e-104
[65,] 1.0000000 7.894660e-104 3.947330e-104
[66,] 1.0000000 2.022281e-103 1.011140e-103
[67,] 1.0000000 5.531728e-103 2.765864e-103
[68,] 1.0000000 1.194606e-102 5.973030e-103
[69,] 1.0000000 3.273720e-102 1.636860e-102
[70,] 1.0000000 1.112280e-101 5.561402e-102
[71,] 1.0000000 3.737587e-101 1.868793e-101
[72,] 1.0000000 1.131059e-100 5.655294e-101
[73,] 1.0000000 1.412410e-101 7.062052e-102
[74,] 1.0000000 6.390831e-105 3.195415e-105
[75,] 1.0000000 1.951933e-104 9.759667e-105
[76,] 1.0000000 7.527625e-104 3.763813e-104
[77,] 1.0000000 2.090931e-103 1.045466e-103
[78,] 1.0000000 6.444974e-103 3.222487e-103
[79,] 1.0000000 1.801577e-102 9.007883e-103
[80,] 1.0000000 7.429728e-102 3.714864e-102
[81,] 1.0000000 1.501027e-101 7.505133e-102
[82,] 1.0000000 4.212265e-101 2.106133e-101
[83,] 1.0000000 1.157339e-100 5.786693e-101
[84,] 1.0000000 3.164350e-100 1.582175e-100
[85,] 1.0000000 1.032007e-99 5.160034e-100
[86,] 1.0000000 4.152333e-99 2.076167e-99
[87,] 1.0000000 2.144352e-99 1.072176e-99
[88,] 1.0000000 8.947888e-99 4.473944e-99
[89,] 1.0000000 2.911255e-98 1.455628e-98
[90,] 1.0000000 1.187516e-97 5.937582e-98
[91,] 1.0000000 3.192685e-97 1.596342e-97
[92,] 1.0000000 2.022532e-103 1.011266e-103
[93,] 1.0000000 5.237092e-103 2.618546e-103
[94,] 1.0000000 1.736089e-102 8.680445e-103
[95,] 1.0000000 6.937118e-102 3.468559e-102
[96,] 1.0000000 2.445301e-101 1.222651e-101
[97,] 1.0000000 5.297226e-101 2.648613e-101
[98,] 1.0000000 1.790735e-100 8.953675e-101
[99,] 1.0000000 6.419752e-100 3.209876e-100
[100,] 1.0000000 2.350926e-99 1.175463e-99
[101,] 1.0000000 7.633737e-99 3.816868e-99
[102,] 1.0000000 2.667929e-98 1.333964e-98
[103,] 1.0000000 1.010420e-97 5.052099e-98
[104,] 1.0000000 3.996146e-97 1.998073e-97
[105,] 1.0000000 1.457149e-96 7.285746e-97
[106,] 1.0000000 5.756694e-96 2.878347e-96
[107,] 1.0000000 1.998433e-95 9.992167e-96
[108,] 1.0000000 6.250929e-95 3.125465e-95
[109,] 1.0000000 2.136534e-94 1.068267e-94
[110,] 1.0000000 7.454174e-94 3.727087e-94
[111,] 1.0000000 2.678614e-93 1.339307e-93
[112,] 1.0000000 8.515539e-93 4.257769e-93
[113,] 1.0000000 2.837636e-92 1.418818e-92
[114,] 1.0000000 3.367368e-92 1.683684e-92
[115,] 1.0000000 1.245806e-91 6.229032e-92
[116,] 1.0000000 4.294911e-91 2.147456e-91
[117,] 1.0000000 6.236790e-92 3.118395e-92
[118,] 1.0000000 2.393117e-91 1.196559e-91
[119,] 1.0000000 9.533747e-91 4.766874e-91
[120,] 1.0000000 3.201540e-90 1.600770e-90
[121,] 1.0000000 9.927288e-90 4.963644e-90
[122,] 1.0000000 3.410390e-89 1.705195e-89
[123,] 1.0000000 1.205413e-88 6.027065e-89
[124,] 1.0000000 3.791266e-88 1.895633e-88
[125,] 1.0000000 1.333829e-87 6.669146e-88
[126,] 1.0000000 2.484379e-87 1.242190e-87
[127,] 1.0000000 7.886741e-87 3.943371e-87
[128,] 1.0000000 2.624233e-86 1.312116e-86
[129,] 1.0000000 8.453496e-86 4.226748e-86
[130,] 1.0000000 2.974239e-85 1.487119e-85
[131,] 1.0000000 1.085686e-84 5.428431e-85
[132,] 1.0000000 3.601922e-84 1.800961e-84
[133,] 1.0000000 1.132276e-83 5.661378e-84
[134,] 1.0000000 3.654535e-83 1.827267e-83
[135,] 1.0000000 1.334438e-82 6.672188e-83
[136,] 1.0000000 4.313341e-82 2.156670e-82
[137,] 1.0000000 1.547223e-81 7.736113e-82
[138,] 1.0000000 5.209453e-81 2.604727e-81
[139,] 1.0000000 1.889244e-80 9.446218e-81
[140,] 1.0000000 6.788448e-80 3.394224e-80
[141,] 1.0000000 2.192874e-79 1.096437e-79
[142,] 1.0000000 7.631805e-79 3.815902e-79
[143,] 1.0000000 2.465311e-78 1.232655e-78
[144,] 1.0000000 8.055981e-78 4.027990e-78
[145,] 1.0000000 2.813335e-77 1.406668e-77
[146,] 1.0000000 3.592931e-77 1.796465e-77
[147,] 1.0000000 1.157258e-76 5.786288e-77
[148,] 1.0000000 4.073192e-76 2.036596e-76
[149,] 1.0000000 1.359615e-75 6.798073e-76
[150,] 1.0000000 4.361302e-75 2.180651e-75
[151,] 1.0000000 1.397073e-74 6.985367e-75
[152,] 1.0000000 4.545476e-74 2.272738e-74
[153,] 1.0000000 1.528797e-73 7.643983e-74
[154,] 1.0000000 4.871839e-73 2.435919e-73
[155,] 1.0000000 1.549670e-72 7.748351e-73
[156,] 1.0000000 4.919756e-72 2.459878e-72
[157,] 1.0000000 1.558700e-71 7.793500e-72
[158,] 1.0000000 5.266419e-71 2.633209e-71
[159,] 1.0000000 1.072454e-70 5.362268e-71
[160,] 1.0000000 3.496896e-70 1.748448e-70
[161,] 1.0000000 1.099435e-69 5.497175e-70
[162,] 1.0000000 3.859697e-69 1.929849e-69
[163,] 1.0000000 1.582703e-104 7.913516e-105
[164,] 1.0000000 6.727806e-104 3.363903e-104
[165,] 1.0000000 2.853717e-103 1.426859e-103
[166,] 1.0000000 1.296964e-102 6.484818e-103
[167,] 1.0000000 5.475819e-102 2.737910e-102
[168,] 1.0000000 2.593142e-101 1.296571e-101
[169,] 1.0000000 1.089470e-100 5.447349e-101
[170,] 1.0000000 5.189489e-100 2.594745e-100
[171,] 1.0000000 3.774377e-108 1.887188e-108
[172,] 1.0000000 1.652433e-107 8.262164e-108
[173,] 1.0000000 7.435809e-107 3.717904e-107
[174,] 1.0000000 3.337385e-106 1.668693e-106
[175,] 1.0000000 1.290135e-105 6.450674e-106
[176,] 1.0000000 6.166814e-105 3.083407e-105
[177,] 1.0000000 1.097839e-104 5.489193e-105
[178,] 1.0000000 1.258733e-104 6.293665e-105
[179,] 1.0000000 6.039581e-104 3.019791e-104
[180,] 1.0000000 3.070130e-103 1.535065e-103
[181,] 1.0000000 1.562085e-102 7.810424e-103
[182,] 1.0000000 7.466915e-102 3.733457e-102
[183,] 1.0000000 3.496371e-101 1.748185e-101
[184,] 1.0000000 1.541338e-100 7.706690e-101
[185,] 1.0000000 7.121632e-100 3.560816e-100
[186,] 1.0000000 3.121106e-99 1.560553e-99
[187,] 1.0000000 1.493217e-98 7.466087e-99
[188,] 1.0000000 6.504827e-98 3.252414e-98
[189,] 1.0000000 2.937831e-97 1.468915e-97
[190,] 1.0000000 1.233500e-96 6.167502e-97
[191,] 1.0000000 8.304744e-97 4.152372e-97
[192,] 1.0000000 3.536095e-96 1.768048e-96
[193,] 1.0000000 1.498509e-95 7.492547e-96
[194,] 1.0000000 7.112259e-95 3.556130e-95
[195,] 1.0000000 2.999977e-94 1.499989e-94
[196,] 1.0000000 1.258239e-93 6.291197e-94
[197,] 1.0000000 5.245278e-93 2.622639e-93
[198,] 1.0000000 2.078940e-92 1.039470e-92
[199,] 1.0000000 9.897178e-92 4.948589e-92
[200,] 1.0000000 4.301604e-91 2.150802e-91
[201,] 1.0000000 1.911531e-90 9.557653e-91
[202,] 1.0000000 7.845037e-90 3.922518e-90
[203,] 1.0000000 3.189477e-89 1.594738e-89
[204,] 1.0000000 1.104296e-88 5.521478e-89
[205,] 1.0000000 4.583244e-88 2.291622e-88
[206,] 1.0000000 1.895950e-87 9.479752e-88
[207,] 1.0000000 5.908707e-87 2.954353e-87
[208,] 1.0000000 2.648153e-86 1.324076e-86
[209,] 1.0000000 1.159940e-85 5.799698e-86
[210,] 1.0000000 4.757605e-85 2.378802e-85
[211,] 1.0000000 2.177397e-84 1.088698e-84
[212,] 1.0000000 8.003549e-84 4.001774e-84
[213,] 1.0000000 3.258149e-83 1.629075e-83
[214,] 1.0000000 1.437805e-82 7.189025e-83
[215,] 1.0000000 6.065549e-82 3.032774e-82
[216,] 1.0000000 2.445528e-81 1.222764e-81
[217,] 1.0000000 1.022166e-80 5.110831e-81
[218,] 1.0000000 4.373295e-80 2.186648e-80
[219,] 1.0000000 1.812255e-79 9.061274e-80
[220,] 1.0000000 6.757457e-79 3.378728e-79
[221,] 1.0000000 2.985115e-78 1.492557e-78
[222,] 1.0000000 1.182815e-77 5.914076e-78
[223,] 1.0000000 4.877397e-77 2.438698e-77
[224,] 1.0000000 1.918047e-76 9.590237e-77
[225,] 1.0000000 8.273964e-76 4.136982e-76
[226,] 1.0000000 3.506552e-75 1.753276e-75
[227,] 1.0000000 1.363487e-74 6.817436e-75
[228,] 1.0000000 5.842338e-74 2.921169e-74
[229,] 1.0000000 2.256088e-73 1.128044e-73
[230,] 1.0000000 8.850164e-73 4.425082e-73
[231,] 1.0000000 3.390040e-72 1.695020e-72
[232,] 1.0000000 1.412953e-71 7.064767e-72
[233,] 1.0000000 5.369059e-71 2.684529e-71
[234,] 1.0000000 2.031690e-70 1.015845e-70
[235,] 1.0000000 7.655800e-70 3.827900e-70
[236,] 1.0000000 3.000974e-69 1.500487e-69
[237,] 1.0000000 1.121340e-68 5.606698e-69
[238,] 1.0000000 4.648777e-68 2.324389e-68
[239,] 1.0000000 1.157931e-67 5.789657e-68
[240,] 1.0000000 4.256635e-67 2.128317e-67
[241,] 1.0000000 1.553094e-66 7.765470e-67
[242,] 1.0000000 6.407647e-66 3.203824e-66
[243,] 1.0000000 2.310790e-65 1.155395e-65
[244,] 1.0000000 8.264507e-65 4.132253e-65
[245,] 1.0000000 3.353939e-64 1.676969e-64
[246,] 1.0000000 1.184655e-63 5.923277e-64
[247,] 1.0000000 4.322681e-63 2.161340e-63
[248,] 1.0000000 9.347389e-64 4.673695e-64
[249,] 1.0000000 5.349341e-65 2.674670e-65
[250,] 1.0000000 2.003147e-64 1.001574e-64
[251,] 1.0000000 8.150694e-64 4.075347e-64
[252,] 1.0000000 2.166130e-63 1.083065e-63
[253,] 1.0000000 7.681461e-63 3.840730e-63
[254,] 1.0000000 2.846280e-62 1.423140e-62
[255,] 1.0000000 9.957621e-62 4.978810e-62
[256,] 1.0000000 3.445817e-61 1.722909e-61
[257,] 1.0000000 1.394259e-60 6.971296e-61
[258,] 1.0000000 5.636073e-60 2.818037e-60
[259,] 1.0000000 1.914144e-59 9.570719e-60
[260,] 1.0000000 6.942381e-59 3.471191e-59
[261,] 1.0000000 2.321476e-58 1.160738e-58
[262,] 1.0000000 9.170172e-58 4.585086e-58
[263,] 1.0000000 3.015860e-57 1.507930e-57
[264,] 1.0000000 1.155634e-56 5.778172e-57
[265,] 1.0000000 3.734674e-56 1.867337e-56
[266,] 1.0000000 1.188783e-55 5.943916e-56
[267,] 1.0000000 3.722359e-55 1.861180e-55
[268,] 1.0000000 1.453570e-54 7.267850e-55
[269,] 1.0000000 5.104171e-54 2.552086e-54
[270,] 1.0000000 1.518786e-53 7.593928e-54
[271,] 1.0000000 4.110401e-53 2.055200e-53
[272,] 1.0000000 1.437278e-52 7.186389e-53
[273,] 1.0000000 1.929817e-55 9.649085e-56
[274,] 1.0000000 5.966435e-55 2.983217e-55
[275,] 1.0000000 2.033688e-54 1.016844e-54
[276,] 1.0000000 7.364492e-54 3.682246e-54
[277,] 1.0000000 8.749608e-55 4.374804e-55
[278,] 1.0000000 3.237819e-54 1.618909e-54
[279,] 1.0000000 1.309161e-53 6.545804e-54
[280,] 1.0000000 4.801790e-53 2.400895e-53
[281,] 1.0000000 1.698034e-52 8.490169e-53
[282,] 1.0000000 6.559953e-52 3.279977e-52
[283,] 1.0000000 2.161664e-51 1.080832e-51
[284,] 1.0000000 7.184591e-52 3.592295e-52
[285,] 1.0000000 2.342678e-51 1.171339e-51
[286,] 1.0000000 7.904527e-51 3.952263e-51
[287,] 1.0000000 2.909843e-50 1.454922e-50
[288,] 1.0000000 1.063227e-49 5.316133e-50
[289,] 1.0000000 4.157662e-49 2.078831e-49
[290,] 1.0000000 1.636942e-48 8.184711e-49
[291,] 1.0000000 5.901677e-48 2.950838e-48
[292,] 1.0000000 1.858560e-47 9.292799e-48
[293,] 1.0000000 6.643123e-47 3.321562e-47
[294,] 1.0000000 2.361821e-46 1.180911e-46
[295,] 1.0000000 8.001385e-46 4.000693e-46
[296,] 1.0000000 2.824435e-45 1.412217e-45
[297,] 1.0000000 9.894432e-45 4.947216e-45
[298,] 1.0000000 2.945175e-44 1.472587e-44
[299,] 1.0000000 2.409635e-44 1.204817e-44
[300,] 1.0000000 9.230202e-44 4.615101e-44
[301,] 1.0000000 2.671648e-43 1.335824e-43
[302,] 1.0000000 9.829985e-43 4.914992e-43
[303,] 1.0000000 3.402339e-42 1.701170e-42
[304,] 1.0000000 9.354923e-42 4.677461e-42
[305,] 1.0000000 2.465916e-41 1.232958e-41
[306,] 1.0000000 6.530214e-43 3.265107e-43
[307,] 1.0000000 1.660780e-42 8.303900e-43
[308,] 1.0000000 5.326875e-42 2.663438e-42
[309,] 1.0000000 2.016461e-41 1.008231e-41
[310,] 1.0000000 7.588159e-41 3.794080e-41
[311,] 1.0000000 2.730298e-40 1.365149e-40
[312,] 1.0000000 1.008034e-39 5.040168e-40
[313,] 1.0000000 3.790874e-39 1.895437e-39
[314,] 1.0000000 1.326960e-38 6.634799e-39
[315,] 1.0000000 4.617184e-38 2.308592e-38
[316,] 1.0000000 1.596887e-37 7.984437e-38
[317,] 1.0000000 5.489403e-37 2.744702e-37
[318,] 1.0000000 1.228226e-36 6.141130e-37
[319,] 1.0000000 4.189347e-36 2.094673e-36
[320,] 1.0000000 1.460597e-35 7.302987e-36
[321,] 1.0000000 2.979037e-35 1.489518e-35
[322,] 1.0000000 5.407684e-35 2.703842e-35
[323,] 1.0000000 1.819514e-34 9.097570e-35
[324,] 1.0000000 6.082782e-34 3.041391e-34
[325,] 1.0000000 2.020338e-33 1.010169e-33
[326,] 1.0000000 6.666415e-33 3.333207e-33
[327,] 1.0000000 2.185134e-32 1.092567e-32
[328,] 1.0000000 7.549967e-32 3.774984e-32
[329,] 1.0000000 2.519307e-31 1.259653e-31
[330,] 1.0000000 8.461340e-31 4.230670e-31
[331,] 1.0000000 1.030248e-31 5.151239e-32
[332,] 1.0000000 3.425763e-31 1.712881e-31
[333,] 1.0000000 1.146926e-30 5.734630e-31
[334,] 1.0000000 4.071309e-30 2.035654e-30
[335,] 1.0000000 1.425964e-29 7.129821e-30
[336,] 1.0000000 4.627089e-29 2.313545e-29
[337,] 1.0000000 1.490542e-28 7.452711e-29
[338,] 1.0000000 4.766289e-28 2.383145e-28
[339,] 1.0000000 6.262642e-28 3.131321e-28
[340,] 1.0000000 2.012179e-27 1.006089e-27
[341,] 1.0000000 6.824653e-27 3.412326e-27
[342,] 1.0000000 2.160611e-26 1.080306e-26
[343,] 1.0000000 6.787395e-26 3.393697e-26
[344,] 1.0000000 2.276864e-25 1.138432e-25
[345,] 1.0000000 7.627347e-25 3.813674e-25
[346,] 1.0000000 2.539802e-24 1.269901e-24
[347,] 1.0000000 8.371450e-24 4.185725e-24
[348,] 1.0000000 2.584429e-23 1.292215e-23
[349,] 1.0000000 7.776458e-23 3.888229e-23
[350,] 1.0000000 2.320296e-22 1.160148e-22
[351,] 1.0000000 2.333854e-22 1.166927e-22
[352,] 1.0000000 7.031727e-22 3.515864e-22
[353,] 1.0000000 2.100247e-21 1.050124e-21
[354,] 1.0000000 6.217862e-21 3.108931e-21
[355,] 1.0000000 1.915061e-20 9.575304e-21
[356,] 1.0000000 5.570308e-20 2.785154e-20
[357,] 1.0000000 1.605300e-19 8.026500e-20
[358,] 1.0000000 4.044419e-19 2.022209e-19
[359,] 1.0000000 1.246704e-18 6.233520e-19
[360,] 1.0000000 3.512841e-18 1.756421e-18
[361,] 1.0000000 9.801791e-18 4.900895e-18
[362,] 1.0000000 2.707860e-17 1.353930e-17
[363,] 1.0000000 7.405232e-17 3.702616e-17
[364,] 1.0000000 2.116225e-16 1.058112e-16
[365,] 1.0000000 5.699292e-16 2.849646e-16
[366,] 1.0000000 1.134272e-15 5.671359e-16
[367,] 1.0000000 3.273603e-15 1.636802e-15
[368,] 1.0000000 8.725925e-15 4.362962e-15
[369,] 1.0000000 2.292806e-14 1.146403e-14
[370,] 1.0000000 5.956328e-14 2.978164e-14
[371,] 1.0000000 1.571483e-13 7.857414e-14
[372,] 1.0000000 4.302030e-13 2.151015e-13
[373,] 1.0000000 1.093966e-12 5.469831e-13
[374,] 1.0000000 2.735702e-12 1.367851e-12
[375,] 1.0000000 6.755723e-12 3.377862e-12
[376,] 1.0000000 6.368768e-12 3.184384e-12
[377,] 1.0000000 1.893586e-12 9.467932e-13
[378,] 1.0000000 5.058938e-12 2.529469e-12
[379,] 1.0000000 1.395227e-11 6.976137e-12
[380,] 1.0000000 3.419246e-11 1.709623e-11
[381,] 1.0000000 9.269245e-11 4.634622e-11
[382,] 1.0000000 2.430964e-10 1.215482e-10
[383,] 1.0000000 6.307088e-10 3.153544e-10
[384,] 1.0000000 1.661666e-09 8.308328e-10
[385,] 1.0000000 4.307703e-09 2.153852e-09
[386,] 1.0000000 1.090335e-08 5.451674e-09
[387,] 1.0000000 2.686629e-08 1.343314e-08
[388,] 1.0000000 6.555743e-08 3.277872e-08
[389,] 1.0000000 7.543605e-08 3.771802e-08
[390,] 0.9999999 1.873081e-07 9.365403e-08
[391,] 0.9999998 4.580422e-07 2.290211e-07
[392,] 1.0000000 2.881856e-08 1.440928e-08
[393,] 1.0000000 7.794760e-08 3.897380e-08
[394,] 1.0000000 3.844653e-10 1.922326e-10
[395,] 1.0000000 1.182595e-09 5.912976e-10
[396,] 1.0000000 3.839402e-09 1.919701e-09
[397,] 1.0000000 1.208336e-08 6.041678e-09
[398,] 1.0000000 3.820630e-08 1.910315e-08
[399,] 0.9999999 1.155342e-07 5.776709e-08
[400,] 0.9999998 3.530822e-07 1.765411e-07
[401,] 0.9999995 1.077779e-06 5.388894e-07
[402,] 0.9999985 3.058725e-06 1.529363e-06
[403,] 0.9999956 8.725485e-06 4.362742e-06
[404,] 0.9999890 2.193737e-05 1.096869e-05
[405,] 0.9999886 2.289329e-05 1.144664e-05
[406,] 0.9999822 3.569126e-05 1.784563e-05
[407,] 0.9999495 1.010669e-04 5.053346e-05
[408,] 0.9998577 2.845138e-04 1.422569e-04
[409,] 0.9996159 7.681873e-04 3.840936e-04
[410,] 0.9990701 1.859868e-03 9.299341e-04
[411,] 0.9984481 3.103792e-03 1.551896e-03
[412,] 0.9964145 7.171099e-03 3.585549e-03
[413,] 0.9912788 1.744233e-02 8.721167e-03
[414,] 0.9802293 3.954138e-02 1.977069e-02
[415,] 0.9595258 8.094848e-02 4.047424e-02
[416,] 0.9419372 1.161257e-01 5.806285e-02
[417,] 0.9966356 6.728711e-03 3.364355e-03
[418,] 0.9927857 1.442867e-02 7.214334e-03
> postscript(file="/var/www/html/freestat/rcomp/tmp/1u6us1291220076.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/freestat/rcomp/tmp/2u6us1291220076.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/freestat/rcomp/tmp/34xtd1291220076.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/freestat/rcomp/tmp/44xtd1291220076.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/freestat/rcomp/tmp/54xtd1291220076.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
-110.69016704 76.20183813 -198.35726405 32.18593712 -119.62313933
6 7 8 9 10
-197.63371809 152.55553020 -52.66323880 175.08514330 -16.98020945
11 12 13 14 15
-162.92800712 9.60893054 123.38933829 -55.85251084 299.88109729
16 17 18 19 20
53.76037664 -61.16938199 324.78296210 -18.09394600 204.26993947
21 22 23 24 25
-0.14075448 -97.34142370 -55.36138811 140.22483364 3.64014479
26 27 28 29 30
-64.79067361 66.64508715 94.33399176 -28.42699075 248.01056012
31 32 33 34 35
-28.08308894 -36.85159740 -32.64306179 -21.57673949 -120.96825026
36 37 38 39 40
-9.56799380 -11.59085419 0.25433282 -17.08342276 88.76214939
41 42 43 44 45
746.10375864 -68.10267497 160.40362848 12.90407123 6.73498924
46 47 48 49 50
-39.79874751 95.38393596 215.68730631 3.27804151 -9.61429376
51 52 53 54 55
43.96814327 67.32584201 -13.14017387 94.64760777 62.10614343
56 57 58 59 60
60.44356341 -31.88337792 -2.60685551 60.82967578 4.40757968
61 62 63 64 65
-27.19017879 7.70430889 -23.76812471 4.61505310 41.51357359
66 67 68 69 70
823.56394588 38.94125251 -0.09452463 32.73688905 -3.46991614
71 72 73 74 75
30.88285221 3.88226326 -16.46925357 -21.17473373 -24.52464286
76 77 78 79 80
-20.94088473 -4.15550133 11.44752161 84.29676625 119.50341882
81 82 83 84 85
-3.58578876 8.72014271 -2.06224013 1.80700853 -27.39095318
86 87 88 89 90
10.42314360 27.97804589 -29.11978499 9.12444860 -32.02767862
91 92 93 94 95
2.73727777 -10.33048993 62.32726696 -20.31860007 22.45742754
96 97 98 99 100
1.08749092 -37.34520097 162.76713732 -46.55537707 -25.38148939
101 102 103 104 105
-3.84872881 -17.06199162 15.69870248 -29.74897249 4.88823991
106 107 108 109 110
-7.02061610 17.34478646 -31.90698393 -16.00784045 -3.98412592
111 112 113 114 115
0.31176426 7.96726680 4.74018138 -29.61366572 -20.40177737
116 117 118 119 120
10.64646197 -7.69671122 -26.76185198 -27.93134251 42.97035273
121 122 123 124 125
-22.02258285 -14.85454668 61.65472597 -24.62410114 -3.34230158
126 127 128 129 130
-27.50583992 -20.11247676 -13.78243690 6.27250938 -20.11247676
131 132 133 134 135
-12.30272175 34.59884515 -20.11247676 -17.47872541 -19.48543151
136 137 138 139 140
-13.15869135 -7.99281927 -18.14637564 -21.19393476 -20.11247676
141 142 143 144 145
-7.97578100 -20.11247676 -10.12957309 -17.43439889 -9.40608528
146 147 148 149 150
-9.33499490 -20.11247676 -14.22494748 -19.94979574 -19.33982849
151 152 153 154 155
-12.71181036 41.04061200 -20.11247676 -11.28027895 -17.34560876
156 157 158 159 160
-20.11247676 -20.11247676 -19.12140434 -15.91421245 -20.11247676
161 162 163 164 165
-20.11247676 -20.11247676 -20.11247676 3.56367634 26.67262365
166 167 168 169 170
-17.43112618 -20.11247676 0.48033039 341.32767152 -20.11247676
171 172 173 174 175
-20.11247676 -16.73416702 -20.11247676 -13.03372267 -20.11247676
176 177 178 179 180
-11.07008128 129.17561331 -20.90447380 -20.11247676 -20.11247676
181 182 183 184 185
2.41091417 -19.89820536 22.40735916 23.03529987 -17.38886597
186 187 188 189 190
-11.58082544 -9.86352145 -16.71369858 -18.89296685 -20.11247676
191 192 193 194 195
-17.43106839 -20.11247676 -15.47557277 -20.11247676 -18.37837103
196 197 198 199 200
-21.41192441 12.20083126 -21.41192441 -21.41192441 -13.25827533
201 202 203 204 205
-21.41192441 -21.41192441 -21.41192441 8.64834965 -19.70725399
206 207 208 209 210
-19.89850271 -20.06180854 -21.41192441 8.34382106 7.79143527
211 212 213 214 215
-20.11247676 -20.11247676 11.45536194 -3.88849083 -17.58296524
216 217 218 219 220
-20.11247676 -12.26383646 5.22906136 -20.11247676 -1.75240558
221 222 223 224 225
-17.97109473 -20.11247676 -18.12297296 -2.28541137 -18.28610255
226 227 228 229 230
3.69671390 -6.52604303 -20.11247676 -19.45261682 -20.11247676
231 232 233 234 235
-14.69115669 -15.11595057 -20.11247676 -6.03963690 -20.11247676
236 237 238 239 240
-19.13632686 -20.11247676 -14.30581382 -20.11247676 -20.11247676
241 242 243 244 245
-20.11247676 -18.03612127 -20.11247676 -11.87162292 -37.27399646
246 247 248 249 250
-21.41192441 -21.41192441 -10.48741492 -21.41192441 -21.41192441
251 252 253 254 255
-9.10428130 -21.41192441 -20.82183746 43.17311371 48.25570585
256 257 258 259 260
-20.11247676 -11.94734535 18.24279301 -21.41192441 -20.11247676
261 262 263 264 265
-21.41192441 -21.41192441 -15.06514006 -6.35988584 -21.41192441
266 267 268 269 270
-20.11247676 -21.41192441 -17.96404950 -21.41192441 -19.11822104
271 272 273 274 275
-21.41192441 -21.41192441 -21.41192441 -8.99484699 -20.11247676
276 277 278 279 280
7.20070375 9.07047943 -20.11247676 83.60502414 3.37625233
281 282 283 284 285
-21.41192441 -20.11247676 51.98633530 -20.11247676 -13.88500776
286 287 288 289 290
-20.11247676 -19.48656906 -16.90869504 -21.41192441 43.06147384
291 292 293 294 295
-21.41192441 1.81989292 -11.31322030 -20.11247676 -5.92864287
296 297 298 299 300
-2.39555089 -20.11247676 -21.41192441 -20.11247676 -20.11247676
301 302 303 304 305
-17.65021603 -19.21510027 -20.11247676 -21.41192441 29.91933801
306 307 308 309 310
-11.13274754 -21.41192441 -11.30357758 1.98513480 -21.41192441
311 312 313 314 315
-21.41192441 48.56893235 -21.41192441 -15.25952839 -16.15116011
316 317 318 319 320
-9.37206450 -19.00192169 -5.80742845 -14.11880912 -20.11247676
321 322 323 324 325
-20.11247676 -20.11247676 -20.11247676 -21.41192441 -20.11247676
326 327 328 329 330
-18.26252664 -21.41192441 -21.41192441 -20.11247676 -20.11247676
331 332 333 334 335
-20.11247676 -20.11247676 -20.11247676 -15.61076659 -18.12238844
336 337 338 339 340
-17.24132605 52.07367320 -20.11247676 -19.36749597 -7.51879731
341 342 343 344 345
-8.17976491 -20.11247676 -20.11247676 -20.11247676 23.67728308
346 347 348 349 350
-20.11247676 -14.20251094 -20.11247676 -20.11247676 -13.92621758
351 352 353 354 355
-7.27835162 -10.70001896 -11.58886256 -19.62951970 -20.11247676
356 357 358 359 360
-20.11247676 20.32240190 -20.11247676 -20.11247676 -20.11247676
361 362 363 364 365
-16.12015137 -20.11247676 -20.11247676 -9.75839269 -9.84739082
366 367 368 369 370
-20.11247676 -20.11247676 -20.11247676 -20.11247676 -0.59046767
371 372 373 374 375
-20.11247676 17.27030692 -19.77612403 -22.12655958 -20.11247676
376 377 378 379 380
-20.11247676 5.09932996 -13.10073418 -21.10933863 -20.11247676
381 382 383 384 385
-20.11247676 30.60625416 64.62588100 -20.11247676 -5.52718114
386 387 388 389 390
-0.15682497 -14.52681465 -17.43063715 -17.45655288 -7.03303959
391 392 393 394 395
-13.69440345 -14.53923133 -17.43240217 -20.10166747 -42.87845938
396 397 398 399 400
-9.90976237 -7.91002328 59.76405401 -16.10833527 119.94515556
401 402 403 404 405
-26.64350743 -17.94400909 -15.57037403 -10.61950067 -18.71442713
406 407 408 409 410
-19.58868900 -11.58213607 -15.01680098 -11.00822184 -26.17765729
411 412 413 414 415
27.96276021 -43.84527101 -10.78237087 -17.34497241 -26.49630715
416 417 418 419 420
-29.36845055 35.20823961 -12.80717347 -6.51955548 -23.92046654
421 422 423 424 425
-29.02949549 22.25341660 4.10553075 -36.16942719 -8.98052266
426 427 428 429 430
38.12343739 -12.98866984 -81.83668364 -43.24588337 -40.88990650
431
52.31642520
> postscript(file="/var/www/html/freestat/rcomp/tmp/6xpby1291220076.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 -110.69016704 NA
1 76.20183813 -110.69016704
2 -198.35726405 76.20183813
3 32.18593712 -198.35726405
4 -119.62313933 32.18593712
5 -197.63371809 -119.62313933
6 152.55553020 -197.63371809
7 -52.66323880 152.55553020
8 175.08514330 -52.66323880
9 -16.98020945 175.08514330
10 -162.92800712 -16.98020945
11 9.60893054 -162.92800712
12 123.38933829 9.60893054
13 -55.85251084 123.38933829
14 299.88109729 -55.85251084
15 53.76037664 299.88109729
16 -61.16938199 53.76037664
17 324.78296210 -61.16938199
18 -18.09394600 324.78296210
19 204.26993947 -18.09394600
20 -0.14075448 204.26993947
21 -97.34142370 -0.14075448
22 -55.36138811 -97.34142370
23 140.22483364 -55.36138811
24 3.64014479 140.22483364
25 -64.79067361 3.64014479
26 66.64508715 -64.79067361
27 94.33399176 66.64508715
28 -28.42699075 94.33399176
29 248.01056012 -28.42699075
30 -28.08308894 248.01056012
31 -36.85159740 -28.08308894
32 -32.64306179 -36.85159740
33 -21.57673949 -32.64306179
34 -120.96825026 -21.57673949
35 -9.56799380 -120.96825026
36 -11.59085419 -9.56799380
37 0.25433282 -11.59085419
38 -17.08342276 0.25433282
39 88.76214939 -17.08342276
40 746.10375864 88.76214939
41 -68.10267497 746.10375864
42 160.40362848 -68.10267497
43 12.90407123 160.40362848
44 6.73498924 12.90407123
45 -39.79874751 6.73498924
46 95.38393596 -39.79874751
47 215.68730631 95.38393596
48 3.27804151 215.68730631
49 -9.61429376 3.27804151
50 43.96814327 -9.61429376
51 67.32584201 43.96814327
52 -13.14017387 67.32584201
53 94.64760777 -13.14017387
54 62.10614343 94.64760777
55 60.44356341 62.10614343
56 -31.88337792 60.44356341
57 -2.60685551 -31.88337792
58 60.82967578 -2.60685551
59 4.40757968 60.82967578
60 -27.19017879 4.40757968
61 7.70430889 -27.19017879
62 -23.76812471 7.70430889
63 4.61505310 -23.76812471
64 41.51357359 4.61505310
65 823.56394588 41.51357359
66 38.94125251 823.56394588
67 -0.09452463 38.94125251
68 32.73688905 -0.09452463
69 -3.46991614 32.73688905
70 30.88285221 -3.46991614
71 3.88226326 30.88285221
72 -16.46925357 3.88226326
73 -21.17473373 -16.46925357
74 -24.52464286 -21.17473373
75 -20.94088473 -24.52464286
76 -4.15550133 -20.94088473
77 11.44752161 -4.15550133
78 84.29676625 11.44752161
79 119.50341882 84.29676625
80 -3.58578876 119.50341882
81 8.72014271 -3.58578876
82 -2.06224013 8.72014271
83 1.80700853 -2.06224013
84 -27.39095318 1.80700853
85 10.42314360 -27.39095318
86 27.97804589 10.42314360
87 -29.11978499 27.97804589
88 9.12444860 -29.11978499
89 -32.02767862 9.12444860
90 2.73727777 -32.02767862
91 -10.33048993 2.73727777
92 62.32726696 -10.33048993
93 -20.31860007 62.32726696
94 22.45742754 -20.31860007
95 1.08749092 22.45742754
96 -37.34520097 1.08749092
97 162.76713732 -37.34520097
98 -46.55537707 162.76713732
99 -25.38148939 -46.55537707
100 -3.84872881 -25.38148939
101 -17.06199162 -3.84872881
102 15.69870248 -17.06199162
103 -29.74897249 15.69870248
104 4.88823991 -29.74897249
105 -7.02061610 4.88823991
106 17.34478646 -7.02061610
107 -31.90698393 17.34478646
108 -16.00784045 -31.90698393
109 -3.98412592 -16.00784045
110 0.31176426 -3.98412592
111 7.96726680 0.31176426
112 4.74018138 7.96726680
113 -29.61366572 4.74018138
114 -20.40177737 -29.61366572
115 10.64646197 -20.40177737
116 -7.69671122 10.64646197
117 -26.76185198 -7.69671122
118 -27.93134251 -26.76185198
119 42.97035273 -27.93134251
120 -22.02258285 42.97035273
121 -14.85454668 -22.02258285
122 61.65472597 -14.85454668
123 -24.62410114 61.65472597
124 -3.34230158 -24.62410114
125 -27.50583992 -3.34230158
126 -20.11247676 -27.50583992
127 -13.78243690 -20.11247676
128 6.27250938 -13.78243690
129 -20.11247676 6.27250938
130 -12.30272175 -20.11247676
131 34.59884515 -12.30272175
132 -20.11247676 34.59884515
133 -17.47872541 -20.11247676
134 -19.48543151 -17.47872541
135 -13.15869135 -19.48543151
136 -7.99281927 -13.15869135
137 -18.14637564 -7.99281927
138 -21.19393476 -18.14637564
139 -20.11247676 -21.19393476
140 -7.97578100 -20.11247676
141 -20.11247676 -7.97578100
142 -10.12957309 -20.11247676
143 -17.43439889 -10.12957309
144 -9.40608528 -17.43439889
145 -9.33499490 -9.40608528
146 -20.11247676 -9.33499490
147 -14.22494748 -20.11247676
148 -19.94979574 -14.22494748
149 -19.33982849 -19.94979574
150 -12.71181036 -19.33982849
151 41.04061200 -12.71181036
152 -20.11247676 41.04061200
153 -11.28027895 -20.11247676
154 -17.34560876 -11.28027895
155 -20.11247676 -17.34560876
156 -20.11247676 -20.11247676
157 -19.12140434 -20.11247676
158 -15.91421245 -19.12140434
159 -20.11247676 -15.91421245
160 -20.11247676 -20.11247676
161 -20.11247676 -20.11247676
162 -20.11247676 -20.11247676
163 3.56367634 -20.11247676
164 26.67262365 3.56367634
165 -17.43112618 26.67262365
166 -20.11247676 -17.43112618
167 0.48033039 -20.11247676
168 341.32767152 0.48033039
169 -20.11247676 341.32767152
170 -20.11247676 -20.11247676
171 -16.73416702 -20.11247676
172 -20.11247676 -16.73416702
173 -13.03372267 -20.11247676
174 -20.11247676 -13.03372267
175 -11.07008128 -20.11247676
176 129.17561331 -11.07008128
177 -20.90447380 129.17561331
178 -20.11247676 -20.90447380
179 -20.11247676 -20.11247676
180 2.41091417 -20.11247676
181 -19.89820536 2.41091417
182 22.40735916 -19.89820536
183 23.03529987 22.40735916
184 -17.38886597 23.03529987
185 -11.58082544 -17.38886597
186 -9.86352145 -11.58082544
187 -16.71369858 -9.86352145
188 -18.89296685 -16.71369858
189 -20.11247676 -18.89296685
190 -17.43106839 -20.11247676
191 -20.11247676 -17.43106839
192 -15.47557277 -20.11247676
193 -20.11247676 -15.47557277
194 -18.37837103 -20.11247676
195 -21.41192441 -18.37837103
196 12.20083126 -21.41192441
197 -21.41192441 12.20083126
198 -21.41192441 -21.41192441
199 -13.25827533 -21.41192441
200 -21.41192441 -13.25827533
201 -21.41192441 -21.41192441
202 -21.41192441 -21.41192441
203 8.64834965 -21.41192441
204 -19.70725399 8.64834965
205 -19.89850271 -19.70725399
206 -20.06180854 -19.89850271
207 -21.41192441 -20.06180854
208 8.34382106 -21.41192441
209 7.79143527 8.34382106
210 -20.11247676 7.79143527
211 -20.11247676 -20.11247676
212 11.45536194 -20.11247676
213 -3.88849083 11.45536194
214 -17.58296524 -3.88849083
215 -20.11247676 -17.58296524
216 -12.26383646 -20.11247676
217 5.22906136 -12.26383646
218 -20.11247676 5.22906136
219 -1.75240558 -20.11247676
220 -17.97109473 -1.75240558
221 -20.11247676 -17.97109473
222 -18.12297296 -20.11247676
223 -2.28541137 -18.12297296
224 -18.28610255 -2.28541137
225 3.69671390 -18.28610255
226 -6.52604303 3.69671390
227 -20.11247676 -6.52604303
228 -19.45261682 -20.11247676
229 -20.11247676 -19.45261682
230 -14.69115669 -20.11247676
231 -15.11595057 -14.69115669
232 -20.11247676 -15.11595057
233 -6.03963690 -20.11247676
234 -20.11247676 -6.03963690
235 -19.13632686 -20.11247676
236 -20.11247676 -19.13632686
237 -14.30581382 -20.11247676
238 -20.11247676 -14.30581382
239 -20.11247676 -20.11247676
240 -20.11247676 -20.11247676
241 -18.03612127 -20.11247676
242 -20.11247676 -18.03612127
243 -11.87162292 -20.11247676
244 -37.27399646 -11.87162292
245 -21.41192441 -37.27399646
246 -21.41192441 -21.41192441
247 -10.48741492 -21.41192441
248 -21.41192441 -10.48741492
249 -21.41192441 -21.41192441
250 -9.10428130 -21.41192441
251 -21.41192441 -9.10428130
252 -20.82183746 -21.41192441
253 43.17311371 -20.82183746
254 48.25570585 43.17311371
255 -20.11247676 48.25570585
256 -11.94734535 -20.11247676
257 18.24279301 -11.94734535
258 -21.41192441 18.24279301
259 -20.11247676 -21.41192441
260 -21.41192441 -20.11247676
261 -21.41192441 -21.41192441
262 -15.06514006 -21.41192441
263 -6.35988584 -15.06514006
264 -21.41192441 -6.35988584
265 -20.11247676 -21.41192441
266 -21.41192441 -20.11247676
267 -17.96404950 -21.41192441
268 -21.41192441 -17.96404950
269 -19.11822104 -21.41192441
270 -21.41192441 -19.11822104
271 -21.41192441 -21.41192441
272 -21.41192441 -21.41192441
273 -8.99484699 -21.41192441
274 -20.11247676 -8.99484699
275 7.20070375 -20.11247676
276 9.07047943 7.20070375
277 -20.11247676 9.07047943
278 83.60502414 -20.11247676
279 3.37625233 83.60502414
280 -21.41192441 3.37625233
281 -20.11247676 -21.41192441
282 51.98633530 -20.11247676
283 -20.11247676 51.98633530
284 -13.88500776 -20.11247676
285 -20.11247676 -13.88500776
286 -19.48656906 -20.11247676
287 -16.90869504 -19.48656906
288 -21.41192441 -16.90869504
289 43.06147384 -21.41192441
290 -21.41192441 43.06147384
291 1.81989292 -21.41192441
292 -11.31322030 1.81989292
293 -20.11247676 -11.31322030
294 -5.92864287 -20.11247676
295 -2.39555089 -5.92864287
296 -20.11247676 -2.39555089
297 -21.41192441 -20.11247676
298 -20.11247676 -21.41192441
299 -20.11247676 -20.11247676
300 -17.65021603 -20.11247676
301 -19.21510027 -17.65021603
302 -20.11247676 -19.21510027
303 -21.41192441 -20.11247676
304 29.91933801 -21.41192441
305 -11.13274754 29.91933801
306 -21.41192441 -11.13274754
307 -11.30357758 -21.41192441
308 1.98513480 -11.30357758
309 -21.41192441 1.98513480
310 -21.41192441 -21.41192441
311 48.56893235 -21.41192441
312 -21.41192441 48.56893235
313 -15.25952839 -21.41192441
314 -16.15116011 -15.25952839
315 -9.37206450 -16.15116011
316 -19.00192169 -9.37206450
317 -5.80742845 -19.00192169
318 -14.11880912 -5.80742845
319 -20.11247676 -14.11880912
320 -20.11247676 -20.11247676
321 -20.11247676 -20.11247676
322 -20.11247676 -20.11247676
323 -21.41192441 -20.11247676
324 -20.11247676 -21.41192441
325 -18.26252664 -20.11247676
326 -21.41192441 -18.26252664
327 -21.41192441 -21.41192441
328 -20.11247676 -21.41192441
329 -20.11247676 -20.11247676
330 -20.11247676 -20.11247676
331 -20.11247676 -20.11247676
332 -20.11247676 -20.11247676
333 -15.61076659 -20.11247676
334 -18.12238844 -15.61076659
335 -17.24132605 -18.12238844
336 52.07367320 -17.24132605
337 -20.11247676 52.07367320
338 -19.36749597 -20.11247676
339 -7.51879731 -19.36749597
340 -8.17976491 -7.51879731
341 -20.11247676 -8.17976491
342 -20.11247676 -20.11247676
343 -20.11247676 -20.11247676
344 23.67728308 -20.11247676
345 -20.11247676 23.67728308
346 -14.20251094 -20.11247676
347 -20.11247676 -14.20251094
348 -20.11247676 -20.11247676
349 -13.92621758 -20.11247676
350 -7.27835162 -13.92621758
351 -10.70001896 -7.27835162
352 -11.58886256 -10.70001896
353 -19.62951970 -11.58886256
354 -20.11247676 -19.62951970
355 -20.11247676 -20.11247676
356 20.32240190 -20.11247676
357 -20.11247676 20.32240190
358 -20.11247676 -20.11247676
359 -20.11247676 -20.11247676
360 -16.12015137 -20.11247676
361 -20.11247676 -16.12015137
362 -20.11247676 -20.11247676
363 -9.75839269 -20.11247676
364 -9.84739082 -9.75839269
365 -20.11247676 -9.84739082
366 -20.11247676 -20.11247676
367 -20.11247676 -20.11247676
368 -20.11247676 -20.11247676
369 -0.59046767 -20.11247676
370 -20.11247676 -0.59046767
371 17.27030692 -20.11247676
372 -19.77612403 17.27030692
373 -22.12655958 -19.77612403
374 -20.11247676 -22.12655958
375 -20.11247676 -20.11247676
376 5.09932996 -20.11247676
377 -13.10073418 5.09932996
378 -21.10933863 -13.10073418
379 -20.11247676 -21.10933863
380 -20.11247676 -20.11247676
381 30.60625416 -20.11247676
382 64.62588100 30.60625416
383 -20.11247676 64.62588100
384 -5.52718114 -20.11247676
385 -0.15682497 -5.52718114
386 -14.52681465 -0.15682497
387 -17.43063715 -14.52681465
388 -17.45655288 -17.43063715
389 -7.03303959 -17.45655288
390 -13.69440345 -7.03303959
391 -14.53923133 -13.69440345
392 -17.43240217 -14.53923133
393 -20.10166747 -17.43240217
394 -42.87845938 -20.10166747
395 -9.90976237 -42.87845938
396 -7.91002328 -9.90976237
397 59.76405401 -7.91002328
398 -16.10833527 59.76405401
399 119.94515556 -16.10833527
400 -26.64350743 119.94515556
401 -17.94400909 -26.64350743
402 -15.57037403 -17.94400909
403 -10.61950067 -15.57037403
404 -18.71442713 -10.61950067
405 -19.58868900 -18.71442713
406 -11.58213607 -19.58868900
407 -15.01680098 -11.58213607
408 -11.00822184 -15.01680098
409 -26.17765729 -11.00822184
410 27.96276021 -26.17765729
411 -43.84527101 27.96276021
412 -10.78237087 -43.84527101
413 -17.34497241 -10.78237087
414 -26.49630715 -17.34497241
415 -29.36845055 -26.49630715
416 35.20823961 -29.36845055
417 -12.80717347 35.20823961
418 -6.51955548 -12.80717347
419 -23.92046654 -6.51955548
420 -29.02949549 -23.92046654
421 22.25341660 -29.02949549
422 4.10553075 22.25341660
423 -36.16942719 4.10553075
424 -8.98052266 -36.16942719
425 38.12343739 -8.98052266
426 -12.98866984 38.12343739
427 -81.83668364 -12.98866984
428 -43.24588337 -81.83668364
429 -40.88990650 -43.24588337
430 52.31642520 -40.88990650
431 NA 52.31642520
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 76.20183813 -110.69016704
[2,] -198.35726405 76.20183813
[3,] 32.18593712 -198.35726405
[4,] -119.62313933 32.18593712
[5,] -197.63371809 -119.62313933
[6,] 152.55553020 -197.63371809
[7,] -52.66323880 152.55553020
[8,] 175.08514330 -52.66323880
[9,] -16.98020945 175.08514330
[10,] -162.92800712 -16.98020945
[11,] 9.60893054 -162.92800712
[12,] 123.38933829 9.60893054
[13,] -55.85251084 123.38933829
[14,] 299.88109729 -55.85251084
[15,] 53.76037664 299.88109729
[16,] -61.16938199 53.76037664
[17,] 324.78296210 -61.16938199
[18,] -18.09394600 324.78296210
[19,] 204.26993947 -18.09394600
[20,] -0.14075448 204.26993947
[21,] -97.34142370 -0.14075448
[22,] -55.36138811 -97.34142370
[23,] 140.22483364 -55.36138811
[24,] 3.64014479 140.22483364
[25,] -64.79067361 3.64014479
[26,] 66.64508715 -64.79067361
[27,] 94.33399176 66.64508715
[28,] -28.42699075 94.33399176
[29,] 248.01056012 -28.42699075
[30,] -28.08308894 248.01056012
[31,] -36.85159740 -28.08308894
[32,] -32.64306179 -36.85159740
[33,] -21.57673949 -32.64306179
[34,] -120.96825026 -21.57673949
[35,] -9.56799380 -120.96825026
[36,] -11.59085419 -9.56799380
[37,] 0.25433282 -11.59085419
[38,] -17.08342276 0.25433282
[39,] 88.76214939 -17.08342276
[40,] 746.10375864 88.76214939
[41,] -68.10267497 746.10375864
[42,] 160.40362848 -68.10267497
[43,] 12.90407123 160.40362848
[44,] 6.73498924 12.90407123
[45,] -39.79874751 6.73498924
[46,] 95.38393596 -39.79874751
[47,] 215.68730631 95.38393596
[48,] 3.27804151 215.68730631
[49,] -9.61429376 3.27804151
[50,] 43.96814327 -9.61429376
[51,] 67.32584201 43.96814327
[52,] -13.14017387 67.32584201
[53,] 94.64760777 -13.14017387
[54,] 62.10614343 94.64760777
[55,] 60.44356341 62.10614343
[56,] -31.88337792 60.44356341
[57,] -2.60685551 -31.88337792
[58,] 60.82967578 -2.60685551
[59,] 4.40757968 60.82967578
[60,] -27.19017879 4.40757968
[61,] 7.70430889 -27.19017879
[62,] -23.76812471 7.70430889
[63,] 4.61505310 -23.76812471
[64,] 41.51357359 4.61505310
[65,] 823.56394588 41.51357359
[66,] 38.94125251 823.56394588
[67,] -0.09452463 38.94125251
[68,] 32.73688905 -0.09452463
[69,] -3.46991614 32.73688905
[70,] 30.88285221 -3.46991614
[71,] 3.88226326 30.88285221
[72,] -16.46925357 3.88226326
[73,] -21.17473373 -16.46925357
[74,] -24.52464286 -21.17473373
[75,] -20.94088473 -24.52464286
[76,] -4.15550133 -20.94088473
[77,] 11.44752161 -4.15550133
[78,] 84.29676625 11.44752161
[79,] 119.50341882 84.29676625
[80,] -3.58578876 119.50341882
[81,] 8.72014271 -3.58578876
[82,] -2.06224013 8.72014271
[83,] 1.80700853 -2.06224013
[84,] -27.39095318 1.80700853
[85,] 10.42314360 -27.39095318
[86,] 27.97804589 10.42314360
[87,] -29.11978499 27.97804589
[88,] 9.12444860 -29.11978499
[89,] -32.02767862 9.12444860
[90,] 2.73727777 -32.02767862
[91,] -10.33048993 2.73727777
[92,] 62.32726696 -10.33048993
[93,] -20.31860007 62.32726696
[94,] 22.45742754 -20.31860007
[95,] 1.08749092 22.45742754
[96,] -37.34520097 1.08749092
[97,] 162.76713732 -37.34520097
[98,] -46.55537707 162.76713732
[99,] -25.38148939 -46.55537707
[100,] -3.84872881 -25.38148939
[101,] -17.06199162 -3.84872881
[102,] 15.69870248 -17.06199162
[103,] -29.74897249 15.69870248
[104,] 4.88823991 -29.74897249
[105,] -7.02061610 4.88823991
[106,] 17.34478646 -7.02061610
[107,] -31.90698393 17.34478646
[108,] -16.00784045 -31.90698393
[109,] -3.98412592 -16.00784045
[110,] 0.31176426 -3.98412592
[111,] 7.96726680 0.31176426
[112,] 4.74018138 7.96726680
[113,] -29.61366572 4.74018138
[114,] -20.40177737 -29.61366572
[115,] 10.64646197 -20.40177737
[116,] -7.69671122 10.64646197
[117,] -26.76185198 -7.69671122
[118,] -27.93134251 -26.76185198
[119,] 42.97035273 -27.93134251
[120,] -22.02258285 42.97035273
[121,] -14.85454668 -22.02258285
[122,] 61.65472597 -14.85454668
[123,] -24.62410114 61.65472597
[124,] -3.34230158 -24.62410114
[125,] -27.50583992 -3.34230158
[126,] -20.11247676 -27.50583992
[127,] -13.78243690 -20.11247676
[128,] 6.27250938 -13.78243690
[129,] -20.11247676 6.27250938
[130,] -12.30272175 -20.11247676
[131,] 34.59884515 -12.30272175
[132,] -20.11247676 34.59884515
[133,] -17.47872541 -20.11247676
[134,] -19.48543151 -17.47872541
[135,] -13.15869135 -19.48543151
[136,] -7.99281927 -13.15869135
[137,] -18.14637564 -7.99281927
[138,] -21.19393476 -18.14637564
[139,] -20.11247676 -21.19393476
[140,] -7.97578100 -20.11247676
[141,] -20.11247676 -7.97578100
[142,] -10.12957309 -20.11247676
[143,] -17.43439889 -10.12957309
[144,] -9.40608528 -17.43439889
[145,] -9.33499490 -9.40608528
[146,] -20.11247676 -9.33499490
[147,] -14.22494748 -20.11247676
[148,] -19.94979574 -14.22494748
[149,] -19.33982849 -19.94979574
[150,] -12.71181036 -19.33982849
[151,] 41.04061200 -12.71181036
[152,] -20.11247676 41.04061200
[153,] -11.28027895 -20.11247676
[154,] -17.34560876 -11.28027895
[155,] -20.11247676 -17.34560876
[156,] -20.11247676 -20.11247676
[157,] -19.12140434 -20.11247676
[158,] -15.91421245 -19.12140434
[159,] -20.11247676 -15.91421245
[160,] -20.11247676 -20.11247676
[161,] -20.11247676 -20.11247676
[162,] -20.11247676 -20.11247676
[163,] 3.56367634 -20.11247676
[164,] 26.67262365 3.56367634
[165,] -17.43112618 26.67262365
[166,] -20.11247676 -17.43112618
[167,] 0.48033039 -20.11247676
[168,] 341.32767152 0.48033039
[169,] -20.11247676 341.32767152
[170,] -20.11247676 -20.11247676
[171,] -16.73416702 -20.11247676
[172,] -20.11247676 -16.73416702
[173,] -13.03372267 -20.11247676
[174,] -20.11247676 -13.03372267
[175,] -11.07008128 -20.11247676
[176,] 129.17561331 -11.07008128
[177,] -20.90447380 129.17561331
[178,] -20.11247676 -20.90447380
[179,] -20.11247676 -20.11247676
[180,] 2.41091417 -20.11247676
[181,] -19.89820536 2.41091417
[182,] 22.40735916 -19.89820536
[183,] 23.03529987 22.40735916
[184,] -17.38886597 23.03529987
[185,] -11.58082544 -17.38886597
[186,] -9.86352145 -11.58082544
[187,] -16.71369858 -9.86352145
[188,] -18.89296685 -16.71369858
[189,] -20.11247676 -18.89296685
[190,] -17.43106839 -20.11247676
[191,] -20.11247676 -17.43106839
[192,] -15.47557277 -20.11247676
[193,] -20.11247676 -15.47557277
[194,] -18.37837103 -20.11247676
[195,] -21.41192441 -18.37837103
[196,] 12.20083126 -21.41192441
[197,] -21.41192441 12.20083126
[198,] -21.41192441 -21.41192441
[199,] -13.25827533 -21.41192441
[200,] -21.41192441 -13.25827533
[201,] -21.41192441 -21.41192441
[202,] -21.41192441 -21.41192441
[203,] 8.64834965 -21.41192441
[204,] -19.70725399 8.64834965
[205,] -19.89850271 -19.70725399
[206,] -20.06180854 -19.89850271
[207,] -21.41192441 -20.06180854
[208,] 8.34382106 -21.41192441
[209,] 7.79143527 8.34382106
[210,] -20.11247676 7.79143527
[211,] -20.11247676 -20.11247676
[212,] 11.45536194 -20.11247676
[213,] -3.88849083 11.45536194
[214,] -17.58296524 -3.88849083
[215,] -20.11247676 -17.58296524
[216,] -12.26383646 -20.11247676
[217,] 5.22906136 -12.26383646
[218,] -20.11247676 5.22906136
[219,] -1.75240558 -20.11247676
[220,] -17.97109473 -1.75240558
[221,] -20.11247676 -17.97109473
[222,] -18.12297296 -20.11247676
[223,] -2.28541137 -18.12297296
[224,] -18.28610255 -2.28541137
[225,] 3.69671390 -18.28610255
[226,] -6.52604303 3.69671390
[227,] -20.11247676 -6.52604303
[228,] -19.45261682 -20.11247676
[229,] -20.11247676 -19.45261682
[230,] -14.69115669 -20.11247676
[231,] -15.11595057 -14.69115669
[232,] -20.11247676 -15.11595057
[233,] -6.03963690 -20.11247676
[234,] -20.11247676 -6.03963690
[235,] -19.13632686 -20.11247676
[236,] -20.11247676 -19.13632686
[237,] -14.30581382 -20.11247676
[238,] -20.11247676 -14.30581382
[239,] -20.11247676 -20.11247676
[240,] -20.11247676 -20.11247676
[241,] -18.03612127 -20.11247676
[242,] -20.11247676 -18.03612127
[243,] -11.87162292 -20.11247676
[244,] -37.27399646 -11.87162292
[245,] -21.41192441 -37.27399646
[246,] -21.41192441 -21.41192441
[247,] -10.48741492 -21.41192441
[248,] -21.41192441 -10.48741492
[249,] -21.41192441 -21.41192441
[250,] -9.10428130 -21.41192441
[251,] -21.41192441 -9.10428130
[252,] -20.82183746 -21.41192441
[253,] 43.17311371 -20.82183746
[254,] 48.25570585 43.17311371
[255,] -20.11247676 48.25570585
[256,] -11.94734535 -20.11247676
[257,] 18.24279301 -11.94734535
[258,] -21.41192441 18.24279301
[259,] -20.11247676 -21.41192441
[260,] -21.41192441 -20.11247676
[261,] -21.41192441 -21.41192441
[262,] -15.06514006 -21.41192441
[263,] -6.35988584 -15.06514006
[264,] -21.41192441 -6.35988584
[265,] -20.11247676 -21.41192441
[266,] -21.41192441 -20.11247676
[267,] -17.96404950 -21.41192441
[268,] -21.41192441 -17.96404950
[269,] -19.11822104 -21.41192441
[270,] -21.41192441 -19.11822104
[271,] -21.41192441 -21.41192441
[272,] -21.41192441 -21.41192441
[273,] -8.99484699 -21.41192441
[274,] -20.11247676 -8.99484699
[275,] 7.20070375 -20.11247676
[276,] 9.07047943 7.20070375
[277,] -20.11247676 9.07047943
[278,] 83.60502414 -20.11247676
[279,] 3.37625233 83.60502414
[280,] -21.41192441 3.37625233
[281,] -20.11247676 -21.41192441
[282,] 51.98633530 -20.11247676
[283,] -20.11247676 51.98633530
[284,] -13.88500776 -20.11247676
[285,] -20.11247676 -13.88500776
[286,] -19.48656906 -20.11247676
[287,] -16.90869504 -19.48656906
[288,] -21.41192441 -16.90869504
[289,] 43.06147384 -21.41192441
[290,] -21.41192441 43.06147384
[291,] 1.81989292 -21.41192441
[292,] -11.31322030 1.81989292
[293,] -20.11247676 -11.31322030
[294,] -5.92864287 -20.11247676
[295,] -2.39555089 -5.92864287
[296,] -20.11247676 -2.39555089
[297,] -21.41192441 -20.11247676
[298,] -20.11247676 -21.41192441
[299,] -20.11247676 -20.11247676
[300,] -17.65021603 -20.11247676
[301,] -19.21510027 -17.65021603
[302,] -20.11247676 -19.21510027
[303,] -21.41192441 -20.11247676
[304,] 29.91933801 -21.41192441
[305,] -11.13274754 29.91933801
[306,] -21.41192441 -11.13274754
[307,] -11.30357758 -21.41192441
[308,] 1.98513480 -11.30357758
[309,] -21.41192441 1.98513480
[310,] -21.41192441 -21.41192441
[311,] 48.56893235 -21.41192441
[312,] -21.41192441 48.56893235
[313,] -15.25952839 -21.41192441
[314,] -16.15116011 -15.25952839
[315,] -9.37206450 -16.15116011
[316,] -19.00192169 -9.37206450
[317,] -5.80742845 -19.00192169
[318,] -14.11880912 -5.80742845
[319,] -20.11247676 -14.11880912
[320,] -20.11247676 -20.11247676
[321,] -20.11247676 -20.11247676
[322,] -20.11247676 -20.11247676
[323,] -21.41192441 -20.11247676
[324,] -20.11247676 -21.41192441
[325,] -18.26252664 -20.11247676
[326,] -21.41192441 -18.26252664
[327,] -21.41192441 -21.41192441
[328,] -20.11247676 -21.41192441
[329,] -20.11247676 -20.11247676
[330,] -20.11247676 -20.11247676
[331,] -20.11247676 -20.11247676
[332,] -20.11247676 -20.11247676
[333,] -15.61076659 -20.11247676
[334,] -18.12238844 -15.61076659
[335,] -17.24132605 -18.12238844
[336,] 52.07367320 -17.24132605
[337,] -20.11247676 52.07367320
[338,] -19.36749597 -20.11247676
[339,] -7.51879731 -19.36749597
[340,] -8.17976491 -7.51879731
[341,] -20.11247676 -8.17976491
[342,] -20.11247676 -20.11247676
[343,] -20.11247676 -20.11247676
[344,] 23.67728308 -20.11247676
[345,] -20.11247676 23.67728308
[346,] -14.20251094 -20.11247676
[347,] -20.11247676 -14.20251094
[348,] -20.11247676 -20.11247676
[349,] -13.92621758 -20.11247676
[350,] -7.27835162 -13.92621758
[351,] -10.70001896 -7.27835162
[352,] -11.58886256 -10.70001896
[353,] -19.62951970 -11.58886256
[354,] -20.11247676 -19.62951970
[355,] -20.11247676 -20.11247676
[356,] 20.32240190 -20.11247676
[357,] -20.11247676 20.32240190
[358,] -20.11247676 -20.11247676
[359,] -20.11247676 -20.11247676
[360,] -16.12015137 -20.11247676
[361,] -20.11247676 -16.12015137
[362,] -20.11247676 -20.11247676
[363,] -9.75839269 -20.11247676
[364,] -9.84739082 -9.75839269
[365,] -20.11247676 -9.84739082
[366,] -20.11247676 -20.11247676
[367,] -20.11247676 -20.11247676
[368,] -20.11247676 -20.11247676
[369,] -0.59046767 -20.11247676
[370,] -20.11247676 -0.59046767
[371,] 17.27030692 -20.11247676
[372,] -19.77612403 17.27030692
[373,] -22.12655958 -19.77612403
[374,] -20.11247676 -22.12655958
[375,] -20.11247676 -20.11247676
[376,] 5.09932996 -20.11247676
[377,] -13.10073418 5.09932996
[378,] -21.10933863 -13.10073418
[379,] -20.11247676 -21.10933863
[380,] -20.11247676 -20.11247676
[381,] 30.60625416 -20.11247676
[382,] 64.62588100 30.60625416
[383,] -20.11247676 64.62588100
[384,] -5.52718114 -20.11247676
[385,] -0.15682497 -5.52718114
[386,] -14.52681465 -0.15682497
[387,] -17.43063715 -14.52681465
[388,] -17.45655288 -17.43063715
[389,] -7.03303959 -17.45655288
[390,] -13.69440345 -7.03303959
[391,] -14.53923133 -13.69440345
[392,] -17.43240217 -14.53923133
[393,] -20.10166747 -17.43240217
[394,] -42.87845938 -20.10166747
[395,] -9.90976237 -42.87845938
[396,] -7.91002328 -9.90976237
[397,] 59.76405401 -7.91002328
[398,] -16.10833527 59.76405401
[399,] 119.94515556 -16.10833527
[400,] -26.64350743 119.94515556
[401,] -17.94400909 -26.64350743
[402,] -15.57037403 -17.94400909
[403,] -10.61950067 -15.57037403
[404,] -18.71442713 -10.61950067
[405,] -19.58868900 -18.71442713
[406,] -11.58213607 -19.58868900
[407,] -15.01680098 -11.58213607
[408,] -11.00822184 -15.01680098
[409,] -26.17765729 -11.00822184
[410,] 27.96276021 -26.17765729
[411,] -43.84527101 27.96276021
[412,] -10.78237087 -43.84527101
[413,] -17.34497241 -10.78237087
[414,] -26.49630715 -17.34497241
[415,] -29.36845055 -26.49630715
[416,] 35.20823961 -29.36845055
[417,] -12.80717347 35.20823961
[418,] -6.51955548 -12.80717347
[419,] -23.92046654 -6.51955548
[420,] -29.02949549 -23.92046654
[421,] 22.25341660 -29.02949549
[422,] 4.10553075 22.25341660
[423,] -36.16942719 4.10553075
[424,] -8.98052266 -36.16942719
[425,] 38.12343739 -8.98052266
[426,] -12.98866984 38.12343739
[427,] -81.83668364 -12.98866984
[428,] -43.24588337 -81.83668364
[429,] -40.88990650 -43.24588337
[430,] 52.31642520 -40.88990650
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 76.20183813 -110.69016704
2 -198.35726405 76.20183813
3 32.18593712 -198.35726405
4 -119.62313933 32.18593712
5 -197.63371809 -119.62313933
6 152.55553020 -197.63371809
7 -52.66323880 152.55553020
8 175.08514330 -52.66323880
9 -16.98020945 175.08514330
10 -162.92800712 -16.98020945
11 9.60893054 -162.92800712
12 123.38933829 9.60893054
13 -55.85251084 123.38933829
14 299.88109729 -55.85251084
15 53.76037664 299.88109729
16 -61.16938199 53.76037664
17 324.78296210 -61.16938199
18 -18.09394600 324.78296210
19 204.26993947 -18.09394600
20 -0.14075448 204.26993947
21 -97.34142370 -0.14075448
22 -55.36138811 -97.34142370
23 140.22483364 -55.36138811
24 3.64014479 140.22483364
25 -64.79067361 3.64014479
26 66.64508715 -64.79067361
27 94.33399176 66.64508715
28 -28.42699075 94.33399176
29 248.01056012 -28.42699075
30 -28.08308894 248.01056012
31 -36.85159740 -28.08308894
32 -32.64306179 -36.85159740
33 -21.57673949 -32.64306179
34 -120.96825026 -21.57673949
35 -9.56799380 -120.96825026
36 -11.59085419 -9.56799380
37 0.25433282 -11.59085419
38 -17.08342276 0.25433282
39 88.76214939 -17.08342276
40 746.10375864 88.76214939
41 -68.10267497 746.10375864
42 160.40362848 -68.10267497
43 12.90407123 160.40362848
44 6.73498924 12.90407123
45 -39.79874751 6.73498924
46 95.38393596 -39.79874751
47 215.68730631 95.38393596
48 3.27804151 215.68730631
49 -9.61429376 3.27804151
50 43.96814327 -9.61429376
51 67.32584201 43.96814327
52 -13.14017387 67.32584201
53 94.64760777 -13.14017387
54 62.10614343 94.64760777
55 60.44356341 62.10614343
56 -31.88337792 60.44356341
57 -2.60685551 -31.88337792
58 60.82967578 -2.60685551
59 4.40757968 60.82967578
60 -27.19017879 4.40757968
61 7.70430889 -27.19017879
62 -23.76812471 7.70430889
63 4.61505310 -23.76812471
64 41.51357359 4.61505310
65 823.56394588 41.51357359
66 38.94125251 823.56394588
67 -0.09452463 38.94125251
68 32.73688905 -0.09452463
69 -3.46991614 32.73688905
70 30.88285221 -3.46991614
71 3.88226326 30.88285221
72 -16.46925357 3.88226326
73 -21.17473373 -16.46925357
74 -24.52464286 -21.17473373
75 -20.94088473 -24.52464286
76 -4.15550133 -20.94088473
77 11.44752161 -4.15550133
78 84.29676625 11.44752161
79 119.50341882 84.29676625
80 -3.58578876 119.50341882
81 8.72014271 -3.58578876
82 -2.06224013 8.72014271
83 1.80700853 -2.06224013
84 -27.39095318 1.80700853
85 10.42314360 -27.39095318
86 27.97804589 10.42314360
87 -29.11978499 27.97804589
88 9.12444860 -29.11978499
89 -32.02767862 9.12444860
90 2.73727777 -32.02767862
91 -10.33048993 2.73727777
92 62.32726696 -10.33048993
93 -20.31860007 62.32726696
94 22.45742754 -20.31860007
95 1.08749092 22.45742754
96 -37.34520097 1.08749092
97 162.76713732 -37.34520097
98 -46.55537707 162.76713732
99 -25.38148939 -46.55537707
100 -3.84872881 -25.38148939
101 -17.06199162 -3.84872881
102 15.69870248 -17.06199162
103 -29.74897249 15.69870248
104 4.88823991 -29.74897249
105 -7.02061610 4.88823991
106 17.34478646 -7.02061610
107 -31.90698393 17.34478646
108 -16.00784045 -31.90698393
109 -3.98412592 -16.00784045
110 0.31176426 -3.98412592
111 7.96726680 0.31176426
112 4.74018138 7.96726680
113 -29.61366572 4.74018138
114 -20.40177737 -29.61366572
115 10.64646197 -20.40177737
116 -7.69671122 10.64646197
117 -26.76185198 -7.69671122
118 -27.93134251 -26.76185198
119 42.97035273 -27.93134251
120 -22.02258285 42.97035273
121 -14.85454668 -22.02258285
122 61.65472597 -14.85454668
123 -24.62410114 61.65472597
124 -3.34230158 -24.62410114
125 -27.50583992 -3.34230158
126 -20.11247676 -27.50583992
127 -13.78243690 -20.11247676
128 6.27250938 -13.78243690
129 -20.11247676 6.27250938
130 -12.30272175 -20.11247676
131 34.59884515 -12.30272175
132 -20.11247676 34.59884515
133 -17.47872541 -20.11247676
134 -19.48543151 -17.47872541
135 -13.15869135 -19.48543151
136 -7.99281927 -13.15869135
137 -18.14637564 -7.99281927
138 -21.19393476 -18.14637564
139 -20.11247676 -21.19393476
140 -7.97578100 -20.11247676
141 -20.11247676 -7.97578100
142 -10.12957309 -20.11247676
143 -17.43439889 -10.12957309
144 -9.40608528 -17.43439889
145 -9.33499490 -9.40608528
146 -20.11247676 -9.33499490
147 -14.22494748 -20.11247676
148 -19.94979574 -14.22494748
149 -19.33982849 -19.94979574
150 -12.71181036 -19.33982849
151 41.04061200 -12.71181036
152 -20.11247676 41.04061200
153 -11.28027895 -20.11247676
154 -17.34560876 -11.28027895
155 -20.11247676 -17.34560876
156 -20.11247676 -20.11247676
157 -19.12140434 -20.11247676
158 -15.91421245 -19.12140434
159 -20.11247676 -15.91421245
160 -20.11247676 -20.11247676
161 -20.11247676 -20.11247676
162 -20.11247676 -20.11247676
163 3.56367634 -20.11247676
164 26.67262365 3.56367634
165 -17.43112618 26.67262365
166 -20.11247676 -17.43112618
167 0.48033039 -20.11247676
168 341.32767152 0.48033039
169 -20.11247676 341.32767152
170 -20.11247676 -20.11247676
171 -16.73416702 -20.11247676
172 -20.11247676 -16.73416702
173 -13.03372267 -20.11247676
174 -20.11247676 -13.03372267
175 -11.07008128 -20.11247676
176 129.17561331 -11.07008128
177 -20.90447380 129.17561331
178 -20.11247676 -20.90447380
179 -20.11247676 -20.11247676
180 2.41091417 -20.11247676
181 -19.89820536 2.41091417
182 22.40735916 -19.89820536
183 23.03529987 22.40735916
184 -17.38886597 23.03529987
185 -11.58082544 -17.38886597
186 -9.86352145 -11.58082544
187 -16.71369858 -9.86352145
188 -18.89296685 -16.71369858
189 -20.11247676 -18.89296685
190 -17.43106839 -20.11247676
191 -20.11247676 -17.43106839
192 -15.47557277 -20.11247676
193 -20.11247676 -15.47557277
194 -18.37837103 -20.11247676
195 -21.41192441 -18.37837103
196 12.20083126 -21.41192441
197 -21.41192441 12.20083126
198 -21.41192441 -21.41192441
199 -13.25827533 -21.41192441
200 -21.41192441 -13.25827533
201 -21.41192441 -21.41192441
202 -21.41192441 -21.41192441
203 8.64834965 -21.41192441
204 -19.70725399 8.64834965
205 -19.89850271 -19.70725399
206 -20.06180854 -19.89850271
207 -21.41192441 -20.06180854
208 8.34382106 -21.41192441
209 7.79143527 8.34382106
210 -20.11247676 7.79143527
211 -20.11247676 -20.11247676
212 11.45536194 -20.11247676
213 -3.88849083 11.45536194
214 -17.58296524 -3.88849083
215 -20.11247676 -17.58296524
216 -12.26383646 -20.11247676
217 5.22906136 -12.26383646
218 -20.11247676 5.22906136
219 -1.75240558 -20.11247676
220 -17.97109473 -1.75240558
221 -20.11247676 -17.97109473
222 -18.12297296 -20.11247676
223 -2.28541137 -18.12297296
224 -18.28610255 -2.28541137
225 3.69671390 -18.28610255
226 -6.52604303 3.69671390
227 -20.11247676 -6.52604303
228 -19.45261682 -20.11247676
229 -20.11247676 -19.45261682
230 -14.69115669 -20.11247676
231 -15.11595057 -14.69115669
232 -20.11247676 -15.11595057
233 -6.03963690 -20.11247676
234 -20.11247676 -6.03963690
235 -19.13632686 -20.11247676
236 -20.11247676 -19.13632686
237 -14.30581382 -20.11247676
238 -20.11247676 -14.30581382
239 -20.11247676 -20.11247676
240 -20.11247676 -20.11247676
241 -18.03612127 -20.11247676
242 -20.11247676 -18.03612127
243 -11.87162292 -20.11247676
244 -37.27399646 -11.87162292
245 -21.41192441 -37.27399646
246 -21.41192441 -21.41192441
247 -10.48741492 -21.41192441
248 -21.41192441 -10.48741492
249 -21.41192441 -21.41192441
250 -9.10428130 -21.41192441
251 -21.41192441 -9.10428130
252 -20.82183746 -21.41192441
253 43.17311371 -20.82183746
254 48.25570585 43.17311371
255 -20.11247676 48.25570585
256 -11.94734535 -20.11247676
257 18.24279301 -11.94734535
258 -21.41192441 18.24279301
259 -20.11247676 -21.41192441
260 -21.41192441 -20.11247676
261 -21.41192441 -21.41192441
262 -15.06514006 -21.41192441
263 -6.35988584 -15.06514006
264 -21.41192441 -6.35988584
265 -20.11247676 -21.41192441
266 -21.41192441 -20.11247676
267 -17.96404950 -21.41192441
268 -21.41192441 -17.96404950
269 -19.11822104 -21.41192441
270 -21.41192441 -19.11822104
271 -21.41192441 -21.41192441
272 -21.41192441 -21.41192441
273 -8.99484699 -21.41192441
274 -20.11247676 -8.99484699
275 7.20070375 -20.11247676
276 9.07047943 7.20070375
277 -20.11247676 9.07047943
278 83.60502414 -20.11247676
279 3.37625233 83.60502414
280 -21.41192441 3.37625233
281 -20.11247676 -21.41192441
282 51.98633530 -20.11247676
283 -20.11247676 51.98633530
284 -13.88500776 -20.11247676
285 -20.11247676 -13.88500776
286 -19.48656906 -20.11247676
287 -16.90869504 -19.48656906
288 -21.41192441 -16.90869504
289 43.06147384 -21.41192441
290 -21.41192441 43.06147384
291 1.81989292 -21.41192441
292 -11.31322030 1.81989292
293 -20.11247676 -11.31322030
294 -5.92864287 -20.11247676
295 -2.39555089 -5.92864287
296 -20.11247676 -2.39555089
297 -21.41192441 -20.11247676
298 -20.11247676 -21.41192441
299 -20.11247676 -20.11247676
300 -17.65021603 -20.11247676
301 -19.21510027 -17.65021603
302 -20.11247676 -19.21510027
303 -21.41192441 -20.11247676
304 29.91933801 -21.41192441
305 -11.13274754 29.91933801
306 -21.41192441 -11.13274754
307 -11.30357758 -21.41192441
308 1.98513480 -11.30357758
309 -21.41192441 1.98513480
310 -21.41192441 -21.41192441
311 48.56893235 -21.41192441
312 -21.41192441 48.56893235
313 -15.25952839 -21.41192441
314 -16.15116011 -15.25952839
315 -9.37206450 -16.15116011
316 -19.00192169 -9.37206450
317 -5.80742845 -19.00192169
318 -14.11880912 -5.80742845
319 -20.11247676 -14.11880912
320 -20.11247676 -20.11247676
321 -20.11247676 -20.11247676
322 -20.11247676 -20.11247676
323 -21.41192441 -20.11247676
324 -20.11247676 -21.41192441
325 -18.26252664 -20.11247676
326 -21.41192441 -18.26252664
327 -21.41192441 -21.41192441
328 -20.11247676 -21.41192441
329 -20.11247676 -20.11247676
330 -20.11247676 -20.11247676
331 -20.11247676 -20.11247676
332 -20.11247676 -20.11247676
333 -15.61076659 -20.11247676
334 -18.12238844 -15.61076659
335 -17.24132605 -18.12238844
336 52.07367320 -17.24132605
337 -20.11247676 52.07367320
338 -19.36749597 -20.11247676
339 -7.51879731 -19.36749597
340 -8.17976491 -7.51879731
341 -20.11247676 -8.17976491
342 -20.11247676 -20.11247676
343 -20.11247676 -20.11247676
344 23.67728308 -20.11247676
345 -20.11247676 23.67728308
346 -14.20251094 -20.11247676
347 -20.11247676 -14.20251094
348 -20.11247676 -20.11247676
349 -13.92621758 -20.11247676
350 -7.27835162 -13.92621758
351 -10.70001896 -7.27835162
352 -11.58886256 -10.70001896
353 -19.62951970 -11.58886256
354 -20.11247676 -19.62951970
355 -20.11247676 -20.11247676
356 20.32240190 -20.11247676
357 -20.11247676 20.32240190
358 -20.11247676 -20.11247676
359 -20.11247676 -20.11247676
360 -16.12015137 -20.11247676
361 -20.11247676 -16.12015137
362 -20.11247676 -20.11247676
363 -9.75839269 -20.11247676
364 -9.84739082 -9.75839269
365 -20.11247676 -9.84739082
366 -20.11247676 -20.11247676
367 -20.11247676 -20.11247676
368 -20.11247676 -20.11247676
369 -0.59046767 -20.11247676
370 -20.11247676 -0.59046767
371 17.27030692 -20.11247676
372 -19.77612403 17.27030692
373 -22.12655958 -19.77612403
374 -20.11247676 -22.12655958
375 -20.11247676 -20.11247676
376 5.09932996 -20.11247676
377 -13.10073418 5.09932996
378 -21.10933863 -13.10073418
379 -20.11247676 -21.10933863
380 -20.11247676 -20.11247676
381 30.60625416 -20.11247676
382 64.62588100 30.60625416
383 -20.11247676 64.62588100
384 -5.52718114 -20.11247676
385 -0.15682497 -5.52718114
386 -14.52681465 -0.15682497
387 -17.43063715 -14.52681465
388 -17.45655288 -17.43063715
389 -7.03303959 -17.45655288
390 -13.69440345 -7.03303959
391 -14.53923133 -13.69440345
392 -17.43240217 -14.53923133
393 -20.10166747 -17.43240217
394 -42.87845938 -20.10166747
395 -9.90976237 -42.87845938
396 -7.91002328 -9.90976237
397 59.76405401 -7.91002328
398 -16.10833527 59.76405401
399 119.94515556 -16.10833527
400 -26.64350743 119.94515556
401 -17.94400909 -26.64350743
402 -15.57037403 -17.94400909
403 -10.61950067 -15.57037403
404 -18.71442713 -10.61950067
405 -19.58868900 -18.71442713
406 -11.58213607 -19.58868900
407 -15.01680098 -11.58213607
408 -11.00822184 -15.01680098
409 -26.17765729 -11.00822184
410 27.96276021 -26.17765729
411 -43.84527101 27.96276021
412 -10.78237087 -43.84527101
413 -17.34497241 -10.78237087
414 -26.49630715 -17.34497241
415 -29.36845055 -26.49630715
416 35.20823961 -29.36845055
417 -12.80717347 35.20823961
418 -6.51955548 -12.80717347
419 -23.92046654 -6.51955548
420 -29.02949549 -23.92046654
421 22.25341660 -29.02949549
422 4.10553075 22.25341660
423 -36.16942719 4.10553075
424 -8.98052266 -36.16942719
425 38.12343739 -8.98052266
426 -12.98866984 38.12343739
427 -81.83668364 -12.98866984
428 -43.24588337 -81.83668364
429 -40.88990650 -43.24588337
430 52.31642520 -40.88990650
> 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/freestat/rcomp/tmp/78ga11291220076.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/freestat/rcomp/tmp/88ga11291220076.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/freestat/rcomp/tmp/98ga11291220076.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/freestat/rcomp/tmp/101p941291220076.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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/11mq8s1291220076.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/freestat/rcomp/tmp/12pqoy1291220076.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/freestat/rcomp/tmp/13l0461291220076.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/freestat/rcomp/tmp/14p12c1291220076.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/freestat/rcomp/tmp/15sjji1291220076.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/freestat/rcomp/tmp/16e2zo1291220076.tab")
+ }
>
> try(system("convert tmp/1u6us1291220076.ps tmp/1u6us1291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/2u6us1291220076.ps tmp/2u6us1291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/34xtd1291220076.ps tmp/34xtd1291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/44xtd1291220076.ps tmp/44xtd1291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/54xtd1291220076.ps tmp/54xtd1291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/6xpby1291220076.ps tmp/6xpby1291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/78ga11291220076.ps tmp/78ga11291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/88ga11291220076.ps tmp/88ga11291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/98ga11291220076.ps tmp/98ga11291220076.png",intern=TRUE))
character(0)
> try(system("convert tmp/101p941291220076.ps tmp/101p941291220076.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.272 3.137 13.735