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(1
+ ,162556
+ ,1081
+ ,213118
+ ,1
+ ,29790
+ ,309
+ ,81767
+ ,1
+ ,87550
+ ,458
+ ,153198
+ ,0
+ ,84738
+ ,588
+ ,-26007
+ ,1
+ ,54660
+ ,299
+ ,126942
+ ,1
+ ,42634
+ ,156
+ ,157214
+ ,0
+ ,40949
+ ,481
+ ,129352
+ ,1
+ ,42312
+ ,323
+ ,234817
+ ,1
+ ,37704
+ ,452
+ ,60448
+ ,1
+ ,16275
+ ,109
+ ,47818
+ ,0
+ ,25830
+ ,115
+ ,245546
+ ,0
+ ,12679
+ ,110
+ ,48020
+ ,1
+ ,18014
+ ,239
+ ,-1710
+ ,0
+ ,43556
+ ,247
+ ,32648
+ ,1
+ ,24524
+ ,497
+ ,95350
+ ,0
+ ,6532
+ ,103
+ ,151352
+ ,0
+ ,7123
+ ,109
+ ,288170
+ ,1
+ ,20813
+ ,502
+ ,114337
+ ,1
+ ,37597
+ ,248
+ ,37884
+ ,0
+ ,17821
+ ,373
+ ,122844
+ ,1
+ ,12988
+ ,119
+ ,82340
+ ,1
+ ,22330
+ ,84
+ ,79801
+ ,0
+ ,13326
+ ,102
+ ,165548
+ ,0
+ ,16189
+ ,295
+ ,116384
+ ,0
+ ,7146
+ ,105
+ ,134028
+ ,0
+ ,15824
+ ,64
+ ,63838
+ ,1
+ ,26088
+ ,267
+ ,74996
+ ,0
+ ,11326
+ ,129
+ ,31080
+ ,0
+ ,8568
+ ,37
+ ,32168
+ ,0
+ ,14416
+ ,361
+ ,49857
+ ,1
+ ,3369
+ ,28
+ ,87161
+ ,1
+ ,11819
+ ,85
+ ,106113
+ ,1
+ ,6620
+ ,44
+ ,80570
+ ,1
+ ,4519
+ ,49
+ ,102129
+ ,0
+ ,2220
+ ,22
+ ,301670
+ ,0
+ ,18562
+ ,155
+ ,102313
+ ,0
+ ,10327
+ ,91
+ ,88577
+ ,1
+ ,5336
+ ,81
+ ,112477
+ ,1
+ ,2365
+ ,79
+ ,191778
+ ,0
+ ,4069
+ ,145
+ ,79804
+ ,0
+ ,7710
+ ,816
+ ,128294
+ ,0
+ ,13718
+ ,61
+ ,96448
+ ,0
+ ,4525
+ ,226
+ ,93811
+ ,0
+ ,6869
+ ,105
+ ,117520
+ ,0
+ ,4628
+ ,62
+ ,69159
+ ,1
+ ,3653
+ ,24
+ ,101792
+ ,1
+ ,1265
+ ,26
+ ,210568
+ ,1
+ ,7489
+ ,322
+ ,136996
+ ,0
+ ,4901
+ ,84
+ ,121920
+ ,0
+ ,2284
+ ,33
+ ,76403
+ ,1
+ ,3160
+ ,108
+ ,108094
+ ,1
+ ,4150
+ ,150
+ ,134759
+ ,1
+ ,7285
+ ,115
+ ,188873
+ ,1
+ ,1134
+ ,162
+ ,146216
+ ,1
+ ,4658
+ ,158
+ ,156608
+ ,0
+ ,2384
+ ,97
+ ,61348
+ ,0
+ ,3748
+ ,9
+ ,50350
+ ,0
+ ,5371
+ ,66
+ ,87720
+ ,0
+ ,1285
+ ,107
+ ,99489
+ ,1
+ ,9327
+ ,101
+ ,87419
+ ,1
+ ,5565
+ ,47
+ ,94355
+ ,0
+ ,1528
+ ,38
+ ,60326
+ ,1
+ ,3122
+ ,34
+ ,94670
+ ,1
+ ,7317
+ ,84
+ ,82425
+ ,0
+ ,2675
+ ,79
+ ,59017
+ ,0
+ ,13253
+ ,947
+ ,90829
+ ,0
+ ,880
+ ,74
+ ,80791
+ ,1
+ ,2053
+ ,53
+ ,100423
+ ,0
+ ,1424
+ ,94
+ ,131116
+ ,1
+ ,4036
+ ,63
+ ,100269
+ ,1
+ ,3045
+ ,58
+ ,27330
+ ,0
+ ,5119
+ ,49
+ ,39039
+ ,0
+ ,1431
+ ,34
+ ,106885
+ ,0
+ ,554
+ ,11
+ ,79285
+ ,0
+ ,1975
+ ,35
+ ,118881
+ ,1
+ ,1286
+ ,17
+ ,77623
+ ,0
+ ,1012
+ ,47
+ ,114768
+ ,0
+ ,810
+ ,43
+ ,74015
+ ,0
+ ,1280
+ ,117
+ ,69465
+ ,1
+ ,666
+ ,171
+ ,117869
+ ,0
+ ,1380
+ ,26
+ ,60982
+ ,1
+ ,4608
+ ,73
+ ,90131
+ ,0
+ ,876
+ ,59
+ ,138971
+ ,0
+ ,814
+ ,18
+ ,39625
+ ,0
+ ,514
+ ,15
+ ,102725
+ ,1
+ ,5692
+ ,72
+ ,64239
+ ,0
+ ,3642
+ ,86
+ ,90262
+ ,0
+ ,540
+ ,14
+ ,103960
+ ,0
+ ,2099
+ ,64
+ ,106611
+ ,0
+ ,567
+ ,11
+ ,103345
+ ,0
+ ,2001
+ ,52
+ ,95551
+ ,1
+ ,2949
+ ,41
+ ,82903
+ ,0
+ ,2253
+ ,99
+ ,63593
+ ,1
+ ,6533
+ ,75
+ ,126910
+ ,0
+ ,1889
+ ,45
+ ,37527
+ ,1
+ ,3055
+ ,43
+ ,60247
+ ,0
+ ,272
+ ,8
+ ,112995
+ ,1
+ ,1414
+ ,198
+ ,70184
+ ,0
+ ,2564
+ ,22
+ ,130140
+ ,1
+ ,1383
+ ,11
+ ,73221
+ ,1
+ ,1261
+ ,33
+ ,76114
+ ,0
+ ,975
+ ,23
+ ,90534
+ ,0
+ ,3366
+ ,80
+ ,108479
+ ,0
+ ,576
+ ,18
+ ,113761
+ ,0
+ ,1306
+ ,28
+ ,68696
+ ,0
+ ,746
+ ,23
+ ,71561
+ ,1
+ ,3192
+ ,60
+ ,59831
+ ,1
+ ,2045
+ ,20
+ ,97890
+ ,0
+ ,5477
+ ,59
+ ,101481
+ ,1
+ ,1932
+ ,36
+ ,72954
+ ,0
+ ,936
+ ,30
+ ,67939
+ ,1
+ ,3437
+ ,47
+ ,48022
+ ,0
+ ,5131
+ ,71
+ ,86111
+ ,1
+ ,2397
+ ,14
+ ,74020
+ ,1
+ ,1389
+ ,9
+ ,57530
+ ,0
+ ,1503
+ ,39
+ ,56364
+ ,0
+ ,402
+ ,26
+ ,84990
+ ,0
+ ,2239
+ ,21
+ ,88590
+ ,1
+ ,2234
+ ,16
+ ,77200
+ ,0
+ ,837
+ ,69
+ ,61262
+ ,0
+ ,10579
+ ,92
+ ,110309
+ ,0
+ ,875
+ ,14
+ ,67000
+ ,0
+ ,1395
+ ,103
+ ,93099
+ ,1
+ ,1659
+ ,29
+ ,107577
+ ,1
+ ,2647
+ ,37
+ ,62920
+ ,1
+ ,3294
+ ,23
+ ,75832
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,94
+ ,7
+ ,60793
+ ,0
+ ,422
+ ,28
+ ,57935
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,34
+ ,8
+ ,60630
+ ,0
+ ,1558
+ ,63
+ ,55637
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,43
+ ,3
+ ,60887
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,316
+ ,9
+ ,60505
+ ,0
+ ,115
+ ,13
+ ,60945
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,389
+ ,14
+ ,58990
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1002
+ ,15
+ ,56750
+ ,0
+ ,36
+ ,3
+ ,60894
+ ,0
+ ,460
+ ,15
+ ,63346
+ ,0
+ ,309
+ ,11
+ ,56535
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,9
+ ,6
+ ,60835
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,14
+ ,1
+ ,61016
+ ,0
+ ,520
+ ,10
+ ,58650
+ ,0
+ ,1766
+ ,73
+ ,60438
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,458
+ ,11
+ ,58625
+ ,0
+ ,20
+ ,3
+ ,60938
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,98
+ ,2
+ ,61490
+ ,0
+ ,405
+ ,7
+ ,60845
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,483
+ ,27
+ ,60830
+ ,0
+ ,454
+ ,51
+ ,63261
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,757
+ ,19
+ ,45689
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,36
+ ,4
+ ,61564
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,203
+ ,9
+ ,61938
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,90
+ ,8
+ ,60951
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,71
+ ,1
+ ,60745
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,972
+ ,34
+ ,71642
+ ,1
+ ,531
+ ,10
+ ,71641
+ ,0
+ ,604
+ ,38
+ ,55792
+ ,1
+ ,283
+ ,10
+ ,71873
+ ,1
+ ,23
+ ,5
+ ,62555
+ ,1
+ ,638
+ ,14
+ ,60370
+ ,1
+ ,699
+ ,16
+ ,64873
+ ,0
+ ,149
+ ,5
+ ,62041
+ ,0
+ ,226
+ ,5
+ ,65745
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,275
+ ,4
+ ,59500
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,141
+ ,6
+ ,61630
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,28
+ ,2
+ ,60890
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,2566
+ ,80
+ ,113521
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,472
+ ,20
+ ,80045
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,203
+ ,27
+ ,50804
+ ,1
+ ,496
+ ,17
+ ,87390
+ ,1
+ ,10
+ ,2
+ ,61656
+ ,1
+ ,63
+ ,4
+ ,65688
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,1136
+ ,32
+ ,48522
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,267
+ ,32
+ ,57640
+ ,0
+ ,474
+ ,20
+ ,61977
+ ,0
+ ,534
+ ,7
+ ,62620
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,15
+ ,8
+ ,60831
+ ,0
+ ,397
+ ,28
+ ,60646
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1061
+ ,20
+ ,56225
+ ,0
+ ,288
+ ,4
+ ,60510
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,3
+ ,2
+ ,60698
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,20
+ ,2
+ ,60805
+ ,0
+ ,278
+ ,26
+ ,61404
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,192
+ ,4
+ ,65276
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,317
+ ,9
+ ,63915
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,368
+ ,17
+ ,61686
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,2
+ ,1
+ ,60743
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,53
+ ,6
+ ,60349
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,94
+ ,3
+ ,61360
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,24
+ ,8
+ ,59818
+ ,0
+ ,2332
+ ,4
+ ,72680
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,131
+ ,11
+ ,61808
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,206
+ ,9
+ ,53110
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,167
+ ,2
+ ,64245
+ ,0
+ ,622
+ ,73
+ ,73007
+ ,0
+ ,885
+ ,85
+ ,82732
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,365
+ ,8
+ ,54820
+ ,0
+ ,364
+ ,35
+ ,47705
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,226
+ ,12
+ ,72835
+ ,0
+ ,307
+ ,15
+ ,58856
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,188
+ ,11
+ ,77655
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,138
+ ,6
+ ,69817
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,125
+ ,12
+ ,60798
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,282
+ ,30
+ ,62452
+ ,0
+ ,335
+ ,33
+ ,64175
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,813
+ ,82
+ ,67440
+ ,0
+ ,176
+ ,28
+ ,68136
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,249
+ ,72
+ ,56726
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,333
+ ,13
+ ,70811
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,30
+ ,4
+ ,62045
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,249
+ ,62
+ ,54323
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,165
+ ,24
+ ,62841
+ ,0
+ ,453
+ ,21
+ ,81125
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,53
+ ,14
+ ,59506
+ ,1
+ ,382
+ ,21
+ ,59365
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,30
+ ,4
+ ,60798
+ ,0
+ ,290
+ ,2
+ ,58790
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,366
+ ,53
+ ,61808
+ ,0
+ ,2
+ ,9
+ ,60735
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,209
+ ,13
+ ,64016
+ ,0
+ ,384
+ ,22
+ ,54683
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,365
+ ,83
+ ,87192
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,49
+ ,8
+ ,64107
+ ,0
+ ,3
+ ,4
+ ,60761
+ ,0
+ ,133
+ ,14
+ ,65990
+ ,0
+ ,32
+ ,1
+ ,59988
+ ,0
+ ,368
+ ,17
+ ,61167
+ ,0
+ ,1
+ ,6
+ ,60719
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,22
+ ,2
+ ,60722
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,96
+ ,5
+ ,60379
+ ,0
+ ,1
+ ,2
+ ,60727
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,81
+ ,7
+ ,60925
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,26
+ ,1
+ ,60896
+ ,0
+ ,125
+ ,13
+ ,59734
+ ,0
+ ,304
+ ,15
+ ,62969
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,119
+ ,6
+ ,59118
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,312
+ ,14
+ ,58598
+ ,0
+ ,60
+ ,10
+ ,61124
+ ,0
+ ,587
+ ,12
+ ,59595
+ ,0
+ ,135
+ ,2
+ ,62065
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,514
+ ,52
+ ,78780
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,1
+ ,4
+ ,60722
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,1
+ ,58
+ ,3
+ ,61600
+ ,0
+ ,180
+ ,11
+ ,59635
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,448
+ ,40
+ ,59781
+ ,0
+ ,227
+ ,9
+ ,76644
+ ,0
+ ,174
+ ,1
+ ,64820
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,121
+ ,24
+ ,56178
+ ,0
+ ,607
+ ,11
+ ,60436
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,530
+ ,60
+ ,73433
+ ,0
+ ,571
+ ,80
+ ,41477
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,78
+ ,16
+ ,62700
+ ,0
+ ,2489
+ ,40
+ ,67804
+ ,0
+ ,131
+ ,6
+ ,59661
+ ,0
+ ,923
+ ,8
+ ,58620
+ ,0
+ ,72
+ ,3
+ ,60398
+ ,0
+ ,572
+ ,16
+ ,58580
+ ,0
+ ,397
+ ,10
+ ,62710
+ ,0
+ ,450
+ ,8
+ ,59325
+ ,0
+ ,622
+ ,7
+ ,60950
+ ,0
+ ,694
+ ,8
+ ,68060
+ ,1
+ ,3425
+ ,12
+ ,83620
+ ,0
+ ,562
+ ,13
+ ,58456
+ ,0
+ ,4917
+ ,42
+ ,52811
+ ,1
+ ,1442
+ ,118
+ ,121173
+ ,0
+ ,529
+ ,9
+ ,63870
+ ,1
+ ,2126
+ ,138
+ ,21001
+ ,0
+ ,1061
+ ,5
+ ,70415
+ ,0
+ ,776
+ ,9
+ ,64230
+ ,0
+ ,611
+ ,8
+ ,59190
+ ,1
+ ,1526
+ ,25
+ ,69351
+ ,0
+ ,592
+ ,7
+ ,64270
+ ,0
+ ,1182
+ ,13
+ ,70694
+ ,0
+ ,621
+ ,16
+ ,68005
+ ,0
+ ,989
+ ,11
+ ,58930
+ ,0
+ ,438
+ ,11
+ ,58320
+ ,0
+ ,726
+ ,3
+ ,69980
+ ,0
+ ,1303
+ ,61
+ ,69863
+ ,1
+ ,6341
+ ,24
+ ,63255
+ ,1
+ ,1164
+ ,17
+ ,57320
+ ,1
+ ,3310
+ ,33
+ ,75230
+ ,0
+ ,1366
+ ,7
+ ,79420
+ ,0
+ ,965
+ ,3
+ ,73490
+ ,0
+ ,3256
+ ,66
+ ,35250
+ ,1
+ ,1135
+ ,17
+ ,62285
+ ,0
+ ,1270
+ ,26
+ ,69206
+ ,0
+ ,661
+ ,3
+ ,65920
+ ,0
+ ,1013
+ ,2
+ ,69770
+ ,0
+ ,2844
+ ,67
+ ,72683
+ ,1
+ ,11528
+ ,70
+ ,-14545
+ ,0
+ ,6526
+ ,26
+ ,55830
+ ,0
+ ,2264
+ ,24
+ ,55174
+ ,1
+ ,4461
+ ,94
+ ,67038
+ ,0
+ ,3999
+ ,30
+ ,51252
+ ,0
+ ,35624
+ ,223
+ ,157278
+ ,0
+ ,9252
+ ,48
+ ,79510
+ ,0
+ ,15236
+ ,90
+ ,77440
+ ,0
+ ,18073
+ ,180
+ ,27284)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Dividends')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Group','Costs','Trades','Dividends'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '4'
> #'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
Dividends Group Costs Trades
1 213118 1 162556 1081
2 81767 1 29790 309
3 153198 1 87550 458
4 -26007 0 84738 588
5 126942 1 54660 299
6 157214 1 42634 156
7 129352 0 40949 481
8 234817 1 42312 323
9 60448 1 37704 452
10 47818 1 16275 109
11 245546 0 25830 115
12 48020 0 12679 110
13 -1710 1 18014 239
14 32648 0 43556 247
15 95350 1 24524 497
16 151352 0 6532 103
17 288170 0 7123 109
18 114337 1 20813 502
19 37884 1 37597 248
20 122844 0 17821 373
21 82340 1 12988 119
22 79801 1 22330 84
23 165548 0 13326 102
24 116384 0 16189 295
25 134028 0 7146 105
26 63838 0 15824 64
27 74996 1 26088 267
28 31080 0 11326 129
29 32168 0 8568 37
30 49857 0 14416 361
31 87161 1 3369 28
32 106113 1 11819 85
33 80570 1 6620 44
34 102129 1 4519 49
35 301670 0 2220 22
36 102313 0 18562 155
37 88577 0 10327 91
38 112477 1 5336 81
39 191778 1 2365 79
40 79804 0 4069 145
41 128294 0 7710 816
42 96448 0 13718 61
43 93811 0 4525 226
44 117520 0 6869 105
45 69159 0 4628 62
46 101792 1 3653 24
47 210568 1 1265 26
48 136996 1 7489 322
49 121920 0 4901 84
50 76403 0 2284 33
51 108094 1 3160 108
52 134759 1 4150 150
53 188873 1 7285 115
54 146216 1 1134 162
55 156608 1 4658 158
56 61348 0 2384 97
57 50350 0 3748 9
58 87720 0 5371 66
59 99489 0 1285 107
60 87419 1 9327 101
61 94355 1 5565 47
62 60326 0 1528 38
63 94670 1 3122 34
64 82425 1 7317 84
65 59017 0 2675 79
66 90829 0 13253 947
67 80791 0 880 74
68 100423 1 2053 53
69 131116 0 1424 94
70 100269 1 4036 63
71 27330 1 3045 58
72 39039 0 5119 49
73 106885 0 1431 34
74 79285 0 554 11
75 118881 0 1975 35
76 77623 1 1286 17
77 114768 0 1012 47
78 74015 0 810 43
79 69465 0 1280 117
80 117869 1 666 171
81 60982 0 1380 26
82 90131 1 4608 73
83 138971 0 876 59
84 39625 0 814 18
85 102725 0 514 15
86 64239 1 5692 72
87 90262 0 3642 86
88 103960 0 540 14
89 106611 0 2099 64
90 103345 0 567 11
91 95551 0 2001 52
92 82903 1 2949 41
93 63593 0 2253 99
94 126910 1 6533 75
95 37527 0 1889 45
96 60247 1 3055 43
97 112995 0 272 8
98 70184 1 1414 198
99 130140 0 2564 22
100 73221 1 1383 11
101 76114 1 1261 33
102 90534 0 975 23
103 108479 0 3366 80
104 113761 0 576 18
105 68696 0 1306 28
106 71561 0 746 23
107 59831 1 3192 60
108 97890 1 2045 20
109 101481 0 5477 59
110 72954 1 1932 36
111 67939 0 936 30
112 48022 1 3437 47
113 86111 0 5131 71
114 74020 1 2397 14
115 57530 1 1389 9
116 56364 0 1503 39
117 84990 0 402 26
118 88590 0 2239 21
119 77200 1 2234 16
120 61262 0 837 69
121 110309 0 10579 92
122 67000 0 875 14
123 93099 0 1395 103
124 107577 1 1659 29
125 62920 1 2647 37
126 75832 1 3294 23
127 60720 0 0 0
128 60793 0 94 7
129 57935 0 422 28
130 60720 0 0 0
131 60630 0 34 8
132 55637 0 1558 63
133 60720 0 0 0
134 60887 0 43 3
135 60720 0 0 0
136 60505 0 316 9
137 60945 0 115 13
138 60720 0 0 0
139 60720 0 0 0
140 60720 0 0 0
141 58990 0 389 14
142 60720 0 0 0
143 56750 0 1002 15
144 60894 0 36 3
145 63346 0 460 15
146 56535 0 309 11
147 60720 0 0 0
148 60835 0 9 6
149 60720 0 0 0
150 61016 0 14 1
151 58650 0 520 10
152 60438 0 1766 73
153 60720 0 0 0
154 58625 0 458 11
155 60938 0 20 3
156 60720 0 0 0
157 60720 0 0 0
158 61490 0 98 2
159 60845 0 405 7
160 60720 0 0 0
161 60720 0 0 0
162 60720 0 0 0
163 60720 0 0 0
164 60830 0 483 27
165 63261 0 454 51
166 60720 0 0 0
167 60720 0 0 0
168 45689 0 757 19
169 60720 0 0 0
170 60720 0 0 0
171 60720 0 0 0
172 61564 0 36 4
173 60720 0 0 0
174 61938 0 203 9
175 60720 0 0 0
176 60951 0 90 8
177 60720 0 0 0
178 60745 1 71 1
179 60720 0 0 0
180 60720 0 0 0
181 71642 0 972 34
182 71641 1 531 10
183 55792 0 604 38
184 71873 1 283 10
185 62555 1 23 5
186 60370 1 638 14
187 64873 1 699 16
188 62041 0 149 5
189 65745 0 226 5
190 60720 0 0 0
191 59500 0 275 4
192 60720 0 0 0
193 61630 0 141 6
194 60720 0 0 0
195 60890 0 28 2
196 60720 1 0 0
197 113521 1 2566 80
198 60720 1 0 0
199 60720 1 0 0
200 80045 1 472 20
201 60720 1 0 0
202 60720 1 0 0
203 60720 1 0 0
204 50804 1 203 27
205 87390 1 496 17
206 61656 1 10 2
207 65688 1 63 4
208 60720 1 0 0
209 48522 1 1136 32
210 60720 1 0 0
211 60720 0 0 0
212 60720 0 0 0
213 57640 0 267 32
214 61977 0 474 20
215 62620 0 534 7
216 60720 0 0 0
217 60831 0 15 8
218 60646 0 397 28
219 60720 0 0 0
220 56225 0 1061 20
221 60510 0 288 4
222 60720 0 0 0
223 60698 0 3 2
224 60720 0 0 0
225 60805 0 20 2
226 61404 0 278 26
227 60720 0 0 0
228 60720 0 0 0
229 65276 0 192 4
230 60720 0 0 0
231 63915 0 317 9
232 60720 0 0 0
233 60720 0 0 0
234 61686 0 368 17
235 60720 0 0 0
236 60743 0 2 1
237 60720 0 0 0
238 60349 0 53 6
239 60720 0 0 0
240 60720 0 0 0
241 60720 0 0 0
242 61360 0 94 3
243 60720 0 0 0
244 59818 0 24 8
245 72680 0 2332 4
246 60720 1 0 0
247 60720 1 0 0
248 61808 0 131 11
249 60720 1 0 0
250 60720 1 0 0
251 53110 0 206 9
252 60720 1 0 0
253 64245 0 167 2
254 73007 0 622 73
255 82732 0 885 85
256 60720 0 0 0
257 54820 0 365 8
258 47705 0 364 35
259 60720 1 0 0
260 60720 0 0 0
261 60720 1 0 0
262 60720 1 0 0
263 72835 0 226 12
264 58856 0 307 15
265 60720 1 0 0
266 60720 0 0 0
267 60720 1 0 0
268 77655 0 188 11
269 60720 1 0 0
270 69817 0 138 6
271 60720 1 0 0
272 60720 1 0 0
273 60720 1 0 0
274 60798 0 125 12
275 60720 0 0 0
276 62452 0 282 30
277 64175 0 335 33
278 60720 0 0 0
279 67440 1 813 82
280 68136 0 176 28
281 60720 1 0 0
282 60720 0 0 0
283 56726 0 249 72
284 60720 0 0 0
285 70811 0 333 13
286 60720 0 0 0
287 60720 1 0 0
288 62045 0 30 4
289 60720 1 0 0
290 54323 0 249 62
291 60720 1 0 0
292 62841 0 165 24
293 81125 0 453 21
294 60720 0 0 0
295 59506 0 53 14
296 59365 1 382 21
297 60720 0 0 0
298 60720 1 0 0
299 60720 0 0 0
300 60720 0 0 0
301 60798 1 30 4
302 58790 0 290 2
303 60720 0 0 0
304 60720 1 0 0
305 61808 0 366 53
306 60735 0 2 9
307 60720 1 0 0
308 64016 1 209 13
309 54683 0 384 22
310 60720 1 0 0
311 60720 1 0 0
312 87192 0 365 83
313 60720 1 0 0
314 64107 1 49 8
315 60761 0 3 4
316 65990 0 133 14
317 59988 0 32 1
318 61167 0 368 17
319 60719 0 1 6
320 60720 0 0 0
321 60720 0 0 0
322 60720 0 0 0
323 60720 0 0 0
324 60720 1 0 0
325 60720 0 0 0
326 60722 0 22 2
327 60720 1 0 0
328 60720 1 0 0
329 60720 0 0 0
330 60720 0 0 0
331 60720 0 0 0
332 60720 0 0 0
333 60720 0 0 0
334 60379 0 96 5
335 60727 0 1 2
336 60720 0 0 0
337 60925 0 81 7
338 60720 0 0 0
339 60896 0 26 1
340 59734 0 125 13
341 62969 0 304 15
342 60720 0 0 0
343 60720 0 0 0
344 60720 0 0 0
345 60720 0 0 0
346 60720 0 0 0
347 59118 0 119 6
348 60720 0 0 0
349 60720 0 0 0
350 60720 0 0 0
351 58598 0 312 14
352 61124 0 60 10
353 59595 0 587 12
354 62065 0 135 2
355 60720 0 0 0
356 60720 0 0 0
357 78780 0 514 52
358 60720 0 0 0
359 60720 0 0 0
360 60720 0 0 0
361 60722 0 1 4
362 60720 0 0 0
363 60720 0 0 0
364 61600 1 58 3
365 59635 0 180 11
366 60720 0 0 0
367 60720 0 0 0
368 60720 0 0 0
369 60720 0 0 0
370 60720 0 0 0
371 60720 0 0 0
372 59781 0 448 40
373 76644 0 227 9
374 64820 0 174 1
375 60720 0 0 0
376 60720 0 0 0
377 56178 0 121 24
378 60436 0 607 11
379 60720 0 0 0
380 60720 0 0 0
381 60720 0 0 0
382 73433 0 530 60
383 41477 0 571 80
384 60720 0 0 0
385 62700 0 78 16
386 67804 0 2489 40
387 59661 0 131 6
388 58620 0 923 8
389 60398 0 72 3
390 58580 0 572 16
391 62710 0 397 10
392 59325 0 450 8
393 60950 0 622 7
394 68060 0 694 8
395 83620 1 3425 12
396 58456 0 562 13
397 52811 0 4917 42
398 121173 1 1442 118
399 63870 0 529 9
400 21001 1 2126 138
401 70415 0 1061 5
402 64230 0 776 9
403 59190 0 611 8
404 69351 1 1526 25
405 64270 0 592 7
406 70694 0 1182 13
407 68005 0 621 16
408 58930 0 989 11
409 58320 0 438 11
410 69980 0 726 3
411 69863 0 1303 61
412 63255 1 6341 24
413 57320 1 1164 17
414 75230 1 3310 33
415 79420 0 1366 7
416 73490 0 965 3
417 35250 0 3256 66
418 62285 1 1135 17
419 69206 0 1270 26
420 65920 0 661 3
421 69770 0 1013 2
422 72683 0 2844 67
423 -14545 1 11528 70
424 55830 0 6526 26
425 55174 0 2264 24
426 67038 1 4461 94
427 51252 0 3999 30
428 157278 0 35624 223
429 79510 0 9252 48
430 77440 0 15236 90
431 27284 0 18073 180
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades
6.615e+04 6.265e+03 3.015e-01 7.238e+01
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-160275.4 -9222.0 -5434.6 770.5 233253.5
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.615e+04 1.770e+03 37.369 < 2e-16 ***
Group 6.265e+03 3.279e+03 1.911 0.056699 .
Costs 3.015e-01 1.910e-01 1.578 0.115219
Trades 7.238e+01 2.034e+01 3.559 0.000414 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 29980 on 427 degrees of freedom
Multiple R-squared: 0.1343, Adjusted R-squared: 0.1282
F-statistic: 22.08 on 3 and 427 DF, p-value: 2.598e-13
> 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.99986646 2.670792e-04 1.335396e-04
[2,] 0.99999060 1.879397e-05 9.396986e-06
[3,] 0.99999990 2.012276e-07 1.006138e-07
[4,] 0.99999998 4.972163e-08 2.486082e-08
[5,] 1.00000000 6.342488e-14 3.171244e-14
[6,] 1.00000000 3.883773e-14 1.941886e-14
[7,] 1.00000000 6.981416e-16 3.490708e-16
[8,] 1.00000000 4.144501e-17 2.072250e-17
[9,] 1.00000000 4.886100e-17 2.443050e-17
[10,] 1.00000000 3.686607e-18 1.843304e-18
[11,] 1.00000000 1.195852e-29 5.979260e-30
[12,] 1.00000000 2.154357e-29 1.077179e-29
[13,] 1.00000000 3.414838e-31 1.707419e-31
[14,] 1.00000000 1.147118e-30 5.735588e-31
[15,] 1.00000000 3.317939e-30 1.658970e-30
[16,] 1.00000000 6.797723e-30 3.398862e-30
[17,] 1.00000000 1.602595e-30 8.012974e-31
[18,] 1.00000000 5.475743e-30 2.737871e-30
[19,] 1.00000000 7.044084e-30 3.522042e-30
[20,] 1.00000000 2.462673e-30 1.231336e-30
[21,] 1.00000000 3.223955e-30 1.611978e-30
[22,] 1.00000000 9.706013e-32 4.853007e-32
[23,] 1.00000000 7.505708e-33 3.752854e-33
[24,] 1.00000000 1.870122e-33 9.350609e-34
[25,] 1.00000000 6.153462e-33 3.076731e-33
[26,] 1.00000000 1.856202e-32 9.281011e-33
[27,] 1.00000000 5.627372e-32 2.813686e-32
[28,] 1.00000000 1.457809e-31 7.289047e-32
[29,] 1.00000000 5.617756e-51 2.808878e-51
[30,] 1.00000000 1.926111e-50 9.630557e-51
[31,] 1.00000000 5.007204e-50 2.503602e-50
[32,] 1.00000000 1.001823e-49 5.009117e-50
[33,] 1.00000000 4.355349e-55 2.177674e-55
[34,] 1.00000000 1.149742e-54 5.748708e-55
[35,] 1.00000000 1.312610e-54 6.563048e-55
[36,] 1.00000000 4.013505e-54 2.006753e-54
[37,] 1.00000000 1.363902e-53 6.819509e-54
[38,] 1.00000000 1.860859e-53 9.304293e-54
[39,] 1.00000000 2.835582e-53 1.417791e-53
[40,] 1.00000000 6.408225e-53 3.204113e-53
[41,] 1.00000000 4.429224e-62 2.214612e-62
[42,] 1.00000000 3.152802e-62 1.576401e-62
[43,] 1.00000000 1.446069e-62 7.230343e-63
[44,] 1.00000000 2.946602e-62 1.473301e-62
[45,] 1.00000000 5.089443e-62 2.544721e-62
[46,] 1.00000000 8.829031e-63 4.414515e-63
[47,] 1.00000000 5.449854e-69 2.724927e-69
[48,] 1.00000000 2.872466e-71 1.436233e-71
[49,] 1.00000000 1.372050e-74 6.860248e-75
[50,] 1.00000000 1.528825e-74 7.644125e-75
[51,] 1.00000000 6.895122e-75 3.447561e-75
[52,] 1.00000000 1.891708e-74 9.458541e-75
[53,] 1.00000000 2.275301e-74 1.137650e-74
[54,] 1.00000000 5.833345e-74 2.916672e-74
[55,] 1.00000000 1.066701e-73 5.333506e-74
[56,] 1.00000000 1.568127e-73 7.840637e-74
[57,] 1.00000000 2.223599e-73 1.111800e-73
[58,] 1.00000000 5.377617e-73 2.688809e-73
[59,] 1.00000000 7.176763e-73 3.588381e-73
[60,] 1.00000000 9.638981e-75 4.819490e-75
[61,] 1.00000000 3.097137e-74 1.548569e-74
[62,] 1.00000000 3.361515e-74 1.680758e-74
[63,] 1.00000000 8.677872e-76 4.338936e-76
[64,] 1.00000000 1.049719e-75 5.248595e-76
[65,] 1.00000000 4.845647e-78 2.422824e-78
[66,] 1.00000000 4.525804e-79 2.262902e-79
[67,] 1.00000000 1.375547e-79 6.877734e-80
[68,] 1.00000000 3.270083e-79 1.635041e-79
[69,] 1.00000000 1.179666e-80 5.898328e-81
[70,] 1.00000000 2.484444e-80 1.242222e-80
[71,] 1.00000000 1.906570e-81 9.532851e-82
[72,] 1.00000000 5.603492e-81 2.801746e-81
[73,] 1.00000000 1.363052e-80 6.815259e-81
[74,] 1.00000000 5.293813e-81 2.646906e-81
[75,] 1.00000000 1.105846e-80 5.529230e-81
[76,] 1.00000000 2.436138e-80 1.218069e-80
[77,] 1.00000000 2.038087e-84 1.019043e-84
[78,] 1.00000000 4.662500e-85 2.331250e-85
[79,] 1.00000000 7.912703e-86 3.956352e-86
[80,] 1.00000000 1.129035e-85 5.645174e-86
[81,] 1.00000000 2.488185e-85 1.244092e-85
[82,] 1.00000000 2.899299e-86 1.449649e-86
[83,] 1.00000000 4.723106e-87 2.361553e-87
[84,] 1.00000000 4.681275e-88 2.340638e-88
[85,] 1.00000000 3.114237e-88 1.557118e-88
[86,] 1.00000000 6.665577e-88 3.332788e-88
[87,] 1.00000000 1.502590e-87 7.512951e-88
[88,] 1.00000000 5.558345e-90 2.779173e-90
[89,] 1.00000000 5.252860e-91 2.626430e-91
[90,] 1.00000000 7.807902e-91 3.903951e-91
[91,] 1.00000000 3.276427e-93 1.638213e-93
[92,] 1.00000000 7.403505e-93 3.701753e-93
[93,] 1.00000000 1.090648e-97 5.453240e-98
[94,] 1.00000000 2.689903e-97 1.344951e-97
[95,] 1.00000000 6.666311e-97 3.333155e-97
[96,] 1.00000000 3.561008e-97 1.780504e-97
[97,] 1.00000000 1.458147e-98 7.290736e-99
[98,] 1.00000000 1.051616e-101 5.258082e-102
[99,] 1.00000000 3.341925e-101 1.670963e-101
[100,] 1.00000000 9.764243e-101 4.882122e-101
[101,] 1.00000000 1.562500e-100 7.812501e-101
[102,] 1.00000000 1.826434e-101 9.132172e-102
[103,] 1.00000000 2.465954e-102 1.232977e-102
[104,] 1.00000000 7.078745e-102 3.539372e-102
[105,] 1.00000000 2.295805e-101 1.147902e-101
[106,] 1.00000000 7.651697e-102 3.825849e-102
[107,] 1.00000000 1.373542e-101 6.867709e-102
[108,] 1.00000000 3.833615e-101 1.916808e-101
[109,] 1.00000000 7.908236e-101 3.954118e-101
[110,] 1.00000000 1.553641e-100 7.768206e-101
[111,] 1.00000000 1.176273e-100 5.881367e-101
[112,] 1.00000000 6.070421e-101 3.035211e-101
[113,] 1.00000000 1.415729e-100 7.078644e-101
[114,] 1.00000000 3.989477e-100 1.994739e-100
[115,] 1.00000000 1.524206e-101 7.621028e-102
[116,] 1.00000000 5.077913e-101 2.538957e-101
[117,] 1.00000000 1.505274e-101 7.526371e-102
[118,] 1.00000000 3.200580e-104 1.600290e-104
[119,] 1.00000000 9.484600e-104 4.742300e-104
[120,] 1.00000000 2.479970e-103 1.239985e-103
[121,] 1.00000000 8.077570e-103 4.038785e-103
[122,] 1.00000000 2.655593e-102 1.327796e-102
[123,] 1.00000000 7.476191e-102 3.738096e-102
[124,] 1.00000000 2.521722e-101 1.260861e-101
[125,] 1.00000000 8.519248e-101 4.259624e-101
[126,] 1.00000000 1.764508e-100 8.822540e-101
[127,] 1.00000000 6.128811e-100 3.064405e-100
[128,] 1.00000000 2.148707e-99 1.074353e-99
[129,] 1.00000000 7.567824e-99 3.783912e-99
[130,] 1.00000000 2.639269e-98 1.319635e-98
[131,] 1.00000000 9.367933e-98 4.683966e-98
[132,] 1.00000000 3.355114e-97 1.677557e-97
[133,] 1.00000000 1.207202e-96 6.036010e-97
[134,] 1.00000000 4.361838e-96 2.180919e-96
[135,] 1.00000000 1.462902e-95 7.314510e-96
[136,] 1.00000000 5.327677e-95 2.663839e-95
[137,] 1.00000000 1.551209e-94 7.756043e-95
[138,] 1.00000000 5.700767e-94 2.850383e-94
[139,] 1.00000000 2.157078e-93 1.078539e-93
[140,] 1.00000000 6.528597e-93 3.264298e-93
[141,] 1.00000000 2.412518e-92 1.206259e-92
[142,] 1.00000000 8.912435e-92 4.456218e-92
[143,] 1.00000000 3.301925e-91 1.650962e-91
[144,] 1.00000000 1.230234e-90 6.151170e-91
[145,] 1.00000000 4.196679e-90 2.098340e-90
[146,] 1.00000000 1.348997e-89 6.744986e-90
[147,] 1.00000000 5.016551e-89 2.508276e-89
[148,] 1.00000000 1.716367e-88 8.581833e-89
[149,] 1.00000000 6.399264e-88 3.199632e-88
[150,] 1.00000000 2.380768e-87 1.190384e-87
[151,] 1.00000000 8.853584e-87 4.426792e-87
[152,] 1.00000000 3.323272e-86 1.661636e-86
[153,] 1.00000000 1.226976e-85 6.134882e-86
[154,] 1.00000000 4.552299e-85 2.276150e-85
[155,] 1.00000000 1.687296e-84 8.436481e-85
[156,] 1.00000000 6.246868e-84 3.123434e-84
[157,] 1.00000000 2.309888e-83 1.154944e-83
[158,] 1.00000000 8.321781e-83 4.160890e-83
[159,] 1.00000000 3.086395e-82 1.543198e-82
[160,] 1.00000000 1.135972e-81 5.679862e-82
[161,] 1.00000000 4.174088e-81 2.087044e-81
[162,] 1.00000000 4.080640e-81 2.040320e-81
[163,] 1.00000000 1.505317e-80 7.526584e-81
[164,] 1.00000000 5.541906e-80 2.770953e-80
[165,] 1.00000000 2.036040e-79 1.018020e-79
[166,] 1.00000000 7.526309e-79 3.763155e-79
[167,] 1.00000000 2.752638e-78 1.376319e-78
[168,] 1.00000000 1.013796e-77 5.068981e-78
[169,] 1.00000000 3.689966e-77 1.844983e-77
[170,] 1.00000000 1.336965e-76 6.684827e-77
[171,] 1.00000000 4.842226e-76 2.421113e-76
[172,] 1.00000000 1.421906e-75 7.109528e-76
[173,] 1.00000000 5.124765e-75 2.562382e-75
[174,] 1.00000000 1.842165e-74 9.210823e-75
[175,] 1.00000000 5.551149e-74 2.775574e-74
[176,] 1.00000000 1.413230e-73 7.066150e-74
[177,] 1.00000000 3.852212e-73 1.926106e-73
[178,] 1.00000000 9.452189e-73 4.726094e-73
[179,] 1.00000000 2.885003e-72 1.442501e-72
[180,] 1.00000000 8.399322e-72 4.199661e-72
[181,] 1.00000000 2.620670e-71 1.310335e-71
[182,] 1.00000000 9.330843e-71 4.665421e-71
[183,] 1.00000000 3.215407e-70 1.607703e-70
[184,] 1.00000000 1.126926e-69 5.634632e-70
[185,] 1.00000000 3.828459e-69 1.914229e-69
[186,] 1.00000000 1.334031e-68 6.670154e-69
[187,] 1.00000000 4.661385e-68 2.330692e-68
[188,] 1.00000000 1.614128e-67 8.070642e-68
[189,] 1.00000000 5.576040e-67 2.788020e-67
[190,] 1.00000000 1.658521e-66 8.292605e-67
[191,] 1.00000000 1.823530e-69 9.117649e-70
[192,] 1.00000000 5.617712e-69 2.808856e-69
[193,] 1.00000000 1.739551e-68 8.697757e-69
[194,] 1.00000000 1.786471e-68 8.932357e-69
[195,] 1.00000000 5.576437e-68 2.788218e-68
[196,] 1.00000000 1.747298e-67 8.736492e-68
[197,] 1.00000000 5.493170e-67 2.746585e-67
[198,] 1.00000000 1.012522e-66 5.062609e-67
[199,] 1.00000000 3.016324e-67 1.508162e-67
[200,] 1.00000000 9.661776e-67 4.830888e-67
[201,] 1.00000000 2.899638e-66 1.449819e-66
[202,] 1.00000000 9.285998e-66 4.642999e-66
[203,] 1.00000000 1.329761e-65 6.648804e-66
[204,] 1.00000000 4.305569e-65 2.152785e-65
[205,] 1.00000000 1.502536e-64 7.512678e-65
[206,] 1.00000000 5.225809e-64 2.612904e-64
[207,] 1.00000000 1.619760e-63 8.098799e-64
[208,] 1.00000000 5.611928e-63 2.805964e-63
[209,] 1.00000000 1.954042e-62 9.770210e-63
[210,] 1.00000000 6.707140e-62 3.353570e-62
[211,] 1.00000000 2.288682e-61 1.144341e-61
[212,] 1.00000000 7.649624e-61 3.824812e-61
[213,] 1.00000000 2.598224e-60 1.299112e-60
[214,] 1.00000000 7.228172e-60 3.614086e-60
[215,] 1.00000000 2.421554e-59 1.210777e-59
[216,] 1.00000000 8.145345e-59 4.072673e-59
[217,] 1.00000000 2.726518e-58 1.363259e-58
[218,] 1.00000000 9.103215e-58 4.551608e-58
[219,] 1.00000000 3.028019e-57 1.514010e-57
[220,] 1.00000000 9.983758e-57 4.991879e-57
[221,] 1.00000000 3.295912e-56 1.647956e-56
[222,] 1.00000000 1.083910e-55 5.419550e-56
[223,] 1.00000000 3.501305e-55 1.750652e-55
[224,] 1.00000000 1.142686e-54 5.713432e-55
[225,] 1.00000000 3.732220e-54 1.866110e-54
[226,] 1.00000000 1.208543e-53 6.042713e-54
[227,] 1.00000000 3.898061e-53 1.949031e-53
[228,] 1.00000000 1.253667e-52 6.268333e-53
[229,] 1.00000000 4.011695e-52 2.005848e-52
[230,] 1.00000000 1.278431e-51 6.392157e-52
[231,] 1.00000000 4.058308e-51 2.029154e-51
[232,] 1.00000000 1.273463e-50 6.367315e-51
[233,] 1.00000000 4.010154e-50 2.005077e-50
[234,] 1.00000000 1.257665e-49 6.288327e-50
[235,] 1.00000000 3.928169e-49 1.964084e-49
[236,] 1.00000000 1.227112e-48 6.135558e-49
[237,] 1.00000000 3.801175e-48 1.900587e-48
[238,] 1.00000000 1.153987e-47 5.769934e-48
[239,] 1.00000000 2.937289e-47 1.468644e-47
[240,] 1.00000000 8.387692e-47 4.193846e-47
[241,] 1.00000000 2.393476e-46 1.196738e-46
[242,] 1.00000000 7.306181e-46 3.653090e-46
[243,] 1.00000000 2.074667e-45 1.037333e-45
[244,] 1.00000000 5.884510e-45 2.942255e-45
[245,] 1.00000000 1.347796e-44 6.738981e-45
[246,] 1.00000000 3.809794e-44 1.904897e-44
[247,] 1.00000000 1.135976e-43 5.679881e-44
[248,] 1.00000000 2.631330e-43 1.315665e-43
[249,] 1.00000000 2.535938e-43 1.267969e-43
[250,] 1.00000000 7.535435e-43 3.767718e-43
[251,] 1.00000000 1.864271e-42 9.321357e-43
[252,] 1.00000000 2.828366e-42 1.414183e-42
[253,] 1.00000000 7.916666e-42 3.958333e-42
[254,] 1.00000000 2.325688e-41 1.162844e-41
[255,] 1.00000000 6.471059e-41 3.235529e-41
[256,] 1.00000000 1.796958e-40 8.984789e-41
[257,] 1.00000000 3.875936e-40 1.937968e-40
[258,] 1.00000000 1.081704e-39 5.408520e-40
[259,] 1.00000000 2.978969e-39 1.489484e-39
[260,] 1.00000000 8.549852e-39 4.274926e-39
[261,] 1.00000000 2.339141e-38 1.169570e-38
[262,] 1.00000000 3.504517e-38 1.752258e-38
[263,] 1.00000000 9.567191e-38 4.783595e-38
[264,] 1.00000000 2.334132e-37 1.167066e-37
[265,] 1.00000000 6.334146e-37 3.167073e-37
[266,] 1.00000000 1.714139e-36 8.570695e-37
[267,] 1.00000000 4.625374e-36 2.312687e-36
[268,] 1.00000000 1.286354e-35 6.431768e-36
[269,] 1.00000000 3.570180e-35 1.785090e-35
[270,] 1.00000000 9.859885e-35 4.929942e-35
[271,] 1.00000000 2.687714e-34 1.343857e-34
[272,] 1.00000000 7.356634e-34 3.678317e-34
[273,] 1.00000000 1.691452e-33 8.457262e-34
[274,] 1.00000000 4.177464e-33 2.088732e-33
[275,] 1.00000000 1.091540e-32 5.457702e-33
[276,] 1.00000000 2.937489e-32 1.468745e-32
[277,] 1.00000000 7.267325e-32 3.633662e-32
[278,] 1.00000000 1.938521e-31 9.692606e-32
[279,] 1.00000000 4.331804e-31 2.165902e-31
[280,] 1.00000000 1.146199e-30 5.730993e-31
[281,] 1.00000000 2.921419e-30 1.460709e-30
[282,] 1.00000000 7.686479e-30 3.843240e-30
[283,] 1.00000000 1.942613e-29 9.713064e-30
[284,] 1.00000000 4.321355e-29 2.160677e-29
[285,] 1.00000000 1.083715e-28 5.418574e-29
[286,] 1.00000000 2.788453e-28 1.394227e-28
[287,] 1.00000000 3.147563e-28 1.573781e-28
[288,] 1.00000000 8.083807e-28 4.041904e-28
[289,] 1.00000000 2.039406e-27 1.019703e-27
[290,] 1.00000000 4.982567e-27 2.491283e-27
[291,] 1.00000000 1.260702e-26 6.303510e-27
[292,] 1.00000000 3.086054e-26 1.543027e-26
[293,] 1.00000000 7.729595e-26 3.864797e-26
[294,] 1.00000000 1.925928e-25 9.629640e-26
[295,] 1.00000000 4.640717e-25 2.320359e-25
[296,] 1.00000000 1.118423e-24 5.592117e-25
[297,] 1.00000000 2.744041e-24 1.372021e-24
[298,] 1.00000000 6.526114e-24 3.263057e-24
[299,] 1.00000000 1.575235e-23 7.876176e-24
[300,] 1.00000000 3.800217e-23 1.900108e-23
[301,] 1.00000000 8.902900e-23 4.451450e-23
[302,] 1.00000000 2.006962e-22 1.003481e-22
[303,] 1.00000000 4.236754e-22 2.118377e-22
[304,] 1.00000000 9.773480e-22 4.886740e-22
[305,] 1.00000000 2.242040e-21 1.121020e-21
[306,] 1.00000000 1.296945e-21 6.484726e-22
[307,] 1.00000000 3.003689e-21 1.501844e-21
[308,] 1.00000000 6.627586e-21 3.313793e-21
[309,] 1.00000000 1.551758e-20 7.758789e-21
[310,] 1.00000000 3.481705e-20 1.740853e-20
[311,] 1.00000000 8.025737e-20 4.012868e-20
[312,] 1.00000000 1.846776e-19 9.233879e-20
[313,] 1.00000000 4.226749e-19 2.113374e-19
[314,] 1.00000000 9.619535e-19 4.809767e-19
[315,] 1.00000000 2.176383e-18 1.088191e-18
[316,] 1.00000000 4.894774e-18 2.447387e-18
[317,] 1.00000000 1.094280e-17 5.471399e-18
[318,] 1.00000000 2.379299e-17 1.189649e-17
[319,] 1.00000000 5.256691e-17 2.628345e-17
[320,] 1.00000000 1.154187e-16 5.770933e-17
[321,] 1.00000000 2.463939e-16 1.231969e-16
[322,] 1.00000000 5.221930e-16 2.610965e-16
[323,] 1.00000000 1.126269e-15 5.631344e-16
[324,] 1.00000000 2.413907e-15 1.206953e-15
[325,] 1.00000000 5.140979e-15 2.570489e-15
[326,] 1.00000000 1.087920e-14 5.439601e-15
[327,] 1.00000000 2.287457e-14 1.143728e-14
[328,] 1.00000000 4.769525e-14 2.384762e-14
[329,] 1.00000000 9.898950e-14 4.949475e-14
[330,] 1.00000000 2.040980e-13 1.020490e-13
[331,] 1.00000000 4.181467e-13 2.090734e-13
[332,] 1.00000000 8.507260e-13 4.253630e-13
[333,] 1.00000000 1.720063e-12 8.600316e-13
[334,] 1.00000000 3.435457e-12 1.717729e-12
[335,] 1.00000000 6.844508e-12 3.422254e-12
[336,] 1.00000000 1.355160e-11 6.775802e-12
[337,] 1.00000000 2.664437e-11 1.332218e-11
[338,] 1.00000000 5.201874e-11 2.600937e-11
[339,] 1.00000000 1.008386e-10 5.041928e-11
[340,] 1.00000000 1.940788e-10 9.703941e-11
[341,] 1.00000000 3.675346e-10 1.837673e-10
[342,] 1.00000000 6.972410e-10 3.486205e-10
[343,] 1.00000000 1.313001e-09 6.565006e-10
[344,] 1.00000000 2.454226e-09 1.227113e-09
[345,] 1.00000000 4.494098e-09 2.247049e-09
[346,] 1.00000000 8.283335e-09 4.141668e-09
[347,] 0.99999999 1.501543e-08 7.507716e-09
[348,] 0.99999999 2.729819e-08 1.364909e-08
[349,] 0.99999998 4.912312e-08 2.456156e-08
[350,] 0.99999996 8.770178e-08 4.385089e-08
[351,] 0.99999995 1.097023e-07 5.485115e-08
[352,] 0.99999990 1.942598e-07 9.712989e-08
[353,] 0.99999983 3.412004e-07 1.706002e-07
[354,] 0.99999970 5.943707e-07 2.971854e-07
[355,] 0.99999949 1.027026e-06 5.135129e-07
[356,] 0.99999912 1.759369e-06 8.796845e-07
[357,] 0.99999851 2.988365e-06 1.494182e-06
[358,] 0.99999752 4.960608e-06 2.480304e-06
[359,] 0.99999587 8.257117e-06 4.128559e-06
[360,] 0.99999316 1.367059e-05 6.835294e-06
[361,] 0.99998878 2.243227e-05 1.121614e-05
[362,] 0.99998176 3.647880e-05 1.823940e-05
[363,] 0.99997061 5.878166e-05 2.939083e-05
[364,] 0.99995308 9.384884e-05 4.692442e-05
[365,] 0.99992578 1.484405e-04 7.422024e-05
[366,] 0.99988392 2.321507e-04 1.160754e-04
[367,] 0.99984466 3.106884e-04 1.553442e-04
[368,] 0.99976078 4.784434e-04 2.392217e-04
[369,] 0.99963389 7.322254e-04 3.661127e-04
[370,] 0.99944525 1.109494e-03 5.547469e-04
[371,] 0.99918039 1.639229e-03 8.196143e-04
[372,] 0.99878442 2.431168e-03 1.215584e-03
[373,] 0.99821294 3.574115e-03 1.787058e-03
[374,] 0.99740029 5.199426e-03 2.599713e-03
[375,] 0.99625816 7.483678e-03 3.741839e-03
[376,] 0.99519885 9.602299e-03 4.801150e-03
[377,] 0.99418712 1.162576e-02 5.812878e-03
[378,] 0.99181087 1.637826e-02 8.189130e-03
[379,] 0.98859103 2.281794e-02 1.140897e-02
[380,] 0.98438212 3.123576e-02 1.561788e-02
[381,] 0.97876602 4.246797e-02 2.123398e-02
[382,] 0.97162922 5.674157e-02 2.837078e-02
[383,] 0.96226446 7.547109e-02 3.773554e-02
[384,] 0.95063051 9.873898e-02 4.936949e-02
[385,] 0.93579413 1.284117e-01 6.420587e-02
[386,] 0.91784583 1.643083e-01 8.215417e-02
[387,] 0.89589799 2.082040e-01 1.041020e-01
[388,] 0.87037653 2.592469e-01 1.296235e-01
[389,] 0.85990310 2.801938e-01 1.400969e-01
[390,] 0.82781286 3.443743e-01 1.721871e-01
[391,] 0.80113201 3.977360e-01 1.988680e-01
[392,] 0.96858913 6.282174e-02 3.141087e-02
[393,] 0.95643867 8.712266e-02 4.356133e-02
[394,] 0.94535819 1.092836e-01 5.464181e-02
[395,] 0.92682078 1.463584e-01 7.317922e-02
[396,] 0.90264690 1.947062e-01 9.735310e-02
[397,] 0.87364129 2.527174e-01 1.263587e-01
[398,] 0.85696837 2.860633e-01 1.430316e-01
[399,] 0.81690432 3.661914e-01 1.830957e-01
[400,] 0.77241274 4.551745e-01 2.275873e-01
[401,] 0.72122502 5.575500e-01 2.787750e-01
[402,] 0.66262930 6.747414e-01 3.373707e-01
[403,] 0.59831538 8.033692e-01 4.016846e-01
[404,] 0.53054181 9.389164e-01 4.694582e-01
[405,] 0.48107349 9.621470e-01 5.189265e-01
[406,] 0.40863956 8.172791e-01 5.913604e-01
[407,] 0.34627983 6.925597e-01 6.537202e-01
[408,] 0.35127520 7.025504e-01 6.487248e-01
[409,] 0.29981203 5.996241e-01 7.001880e-01
[410,] 0.24275265 4.855053e-01 7.572473e-01
[411,] 0.20942236 4.188447e-01 7.905776e-01
[412,] 0.22985232 4.597046e-01 7.701477e-01
[413,] 0.17116875 3.423375e-01 8.288313e-01
[414,] 0.12023811 2.404762e-01 8.797619e-01
[415,] 0.08520058 1.704012e-01 9.147994e-01
[416,] 0.06016729 1.203346e-01 9.398327e-01
[417,] 0.90334825 1.933035e-01 9.665175e-02
[418,] 0.88792378 2.241524e-01 1.120762e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1ibba1291317005.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/2ibba1291317005.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/3tlav1291317005.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/4tlav1291317005.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/5tlav1291317005.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
13434.20093 -22001.95596 21226.85096 -160275.42977 16397.58555
6 7 8 9 10
60646.84247 16032.81249 126258.78154 -56058.28033 -37398.80811
11 12 13 14 15
163278.40317 -29920.12320 -96861.12864 -64519.40595 -20439.27068
16 17 18 19 20
75772.13249 211977.61730 -695.17654 -63823.64434 24316.32064
21 22 23 24 25
-2609.48780 -5432.02860 87991.85379 23994.39941 58118.21857
26 27 28 29 30
-11720.79462 -24616.51957 -47827.43884 -39248.44760 -46775.32554
31 32 33 34 35
11698.98532 23977.07434 2969.53197 24800.14668 233253.49102
36 37 38 39 40
19341.62911 12721.39819 32585.49496 112927.13870 1926.69030
41 42 43 44 45
749.00328 21741.40118 9933.07007 41693.74516 -2878.98424
46 47 48 49 50
26533.88465 135885.19386 39010.69371 48207.24353 7170.96656
51 52 53 54 55
26904.27338 50230.61351 105932.73160 61728.44813 71347.35780
56 57 58 59 60
-12546.77449 -17586.26788 15168.43483 25201.77626 4876.36530
61 62 63 64 65
16855.50417 -9039.99005 18848.16067 1718.99167 -13662.60753
66 67 68 69 70
-47869.75964 9014.57767 23548.20767 57727.85645 22072.41209
71 72 73 74 75
-50205.84089 -32206.04615 37837.79602 12167.08259 49597.37408
76 77 78 79 80
3585.31908 44906.14704 4503.59482 -5544.55776 32871.11131
81 82 83 84 85
-7470.75206 11038.08948 68281.54642 -28078.00703 35329.60750
86 87 88 89 90
-15108.39573 16784.11407 36629.15163 35190.84211 36223.16257
91 92 93 94 95
25029.00311 6626.63789 -10407.04109 47091.85660 -32454.53520
96 97 98 99 100
-16206.09370 46179.26939 -16993.81329 61619.76125 -411.62527
101 102 103 104 105
925.71082 22420.52417 35518.64415 46129.75950 120.79356
106 107 108 109 110
3516.57683 -17893.93562 23406.29769 29404.16080 -2653.77505
111 112 113 114 115
-669.40498 -28835.81868 13269.88357 -135.53936 -15959.66616
116 117 118 119 120
-13066.83572 16832.15408 20240.14593 2948.84331 -10139.53525
121 122 123 124 125
34305.02594 -431.86428 19068.14355 32558.23458 -12975.76034
126 127 128 129 130
754.52169 -5434.63835 -5896.67232 -10373.64507 -5434.63835
131 132 133 134 135
-6113.96409 -15547.64057 -5434.63835 -5497.75708 -5434.63835
136 137 138 139 140
-6396.38254 -6185.30969 -5434.63835 -5434.63835 -5434.63835
141 142 143 144 145
-8295.31583 -5434.63835 -10792.54403 -5488.64630 -4033.10934
146 147 148 149 150
-10509.04011 -5434.63835 -5756.65724 -5434.63835 -5215.24408
151 152 153 154 155
-8385.28088 -11533.20262 -5434.63835 -8463.96957 -5439.82166
156 157 158 159 160
-5434.63835 -5434.63835 -4838.95761 -5938.45126 -5434.63835
161 162 163 164 165
-5434.63835 -5434.63835 -5434.63835 -7424.65483 -6722.13032
166 167 168 169 170
-5434.63835 -5434.63835 -22069.20342 -5434.63835 -5434.63835
171 172 173 174 175
-5434.63835 -4891.03048 -5434.63835 -4929.30852 -5434.63835
176 177 178 179 180
-5809.85033 -5434.63835 -11768.16304 -5434.63835 -5434.63835
181 182 183 184 185
2733.20289 -1662.32900 -13295.36707 -1355.54708 -10233.22581
186 187 188 189 190
-13255.13047 -8915.29276 -4520.48867 -839.70725 -5434.63835
191 192 193 194 195
-7027.09854 -5434.63835 -5001.46052 -5434.63835 -5417.84981
196 197 198 199 200
-11699.36953 34537.14498 -11699.36953 -11699.36953 6035.62013
201 202 203 204 205
-11699.36953 -11699.36953 -11699.36953 -23630.95481 13590.53569
206 207 208 209 210
-10911.15327 -7039.90324 -11699.36953 -26556.21251 -11699.36953
211 212 213 214 215
-5434.63835 -5434.63835 -10911.44305 -5768.25177 -4202.34992
216 217 218 219 220
-5434.63835 -5907.23483 -7655.10657 -5434.63835 -11697.25576
221 222 223 224 225
-6021.01856 -5434.63835 -5602.31131 -5434.63835 -5500.43749
226 227 228 229 230
-6716.45496 -5434.63835 -5434.63835 -1226.07072 -5434.63835
231 232 233 234 235
-2986.68408 -5434.63835 -5434.63835 -5810.13601 -5434.63835
236 237 238 239 240
-5484.62560 -5434.63835 -6255.92500 -5434.63835 -5434.63835
241 242 243 244 245
-5434.63835 -5040.13563 -5434.63835 -6922.94869 5532.63364
246 247 248 249 250
-11699.36953 -11699.36953 -5182.36599 -11699.36953 -11699.36953
251 252 253 254 255
-13758.21314 -11699.36953 -2104.76387 1380.75916 10157.84407
256 257 258 259 260
-5434.63835 -12023.77383 -21092.84495 -11699.36953 -5434.63835
261 262 263 264 265
-11699.36953 -11699.36953 5743.60354 -8476.97372 -11699.36953
266 267 268 269 270
-5434.63835 -11699.36953 10647.44623 -11699.36953 3186.44410
271 272 273 274 275
-11699.36953 -11699.36953 -11699.36953 -6262.94092 -5434.63835
276 277 278 279 280
-5959.19781 -4469.33195 -5434.63835 -11160.02371 -98.46622
281 282 283 284 285
-11699.36953 -5434.63835 -14715.38224 -5434.63835 3614.95458
286 287 288 289 290
-5434.63835 -11699.36953 -4408.22124 -11699.36953 -16394.54051
291 292 293 294 295
-11699.36953 -5100.61259 13313.69640 -5434.63835 -7677.99838
296 297 298 299 300
-14689.62544 -5434.63835 -11699.36953 -5434.63835 -5434.63835
301 302 303 304 305
-11919.95242 -7596.85330 -5434.63835 -11699.36953 -8293.36314
306 307 308 309 310
-6071.69898 -11699.36953 -9407.38564 -13179.88151 -11699.36953
311 312 313 314 315
-11699.36953 14919.41322 -11699.36953 -8906.21837 -5684.07966
316 317 318 319 320
-1218.12158 -6248.67180 -6329.13601 -5870.24492 -5434.63835
321 322 323 324 325
-5434.63835 -5434.63835 -5434.63835 -11699.36953 -5434.63835
326 327 328 329 330
-5584.04057 -11699.36953 -11699.36953 -5434.63835 -5434.63835
331 332 333 334 335
-5434.63835 -5434.63835 -5434.63835 -6166.50705 -5572.70823
336 337 338 339 340
-5434.63835 -5760.75230 -5434.63835 -5338.86256 -7399.32509
341 342 343 344 345
-4363.06910 -5434.63835 -5434.63835 -5434.63835 -5434.63835
346 347 348 349 350
-5434.63835 -7506.82664 -5434.63835 -5434.63835 -5434.63835
351 352 353 354 355
-8664.09725 -5772.57247 -7605.25241 -4275.11459 -5434.63835
356 357 358 359 360
-5434.63835 8706.39311 -5434.63835 -5434.63835 -5434.63835
361 362 363 364 365
-5722.47658 -5434.63835 -5434.63835 -11054.01137 -7370.14145
366 367 368 369 370
-5434.63835 -5434.63835 -5434.63835 -5434.63835 -5434.63835
371 372 373 374 375
-5434.63835 -9404.09518 9769.45452 -1459.49048 -5434.63835
376 377 378 379 380
-5434.63835 -11750.34483 -6697.89903 -5434.63835 -5434.63835
381 382 383 384 385
-5434.63835 2775.49509 -30640.55150 -5434.63835 -4636.30523
386 387 388 389 390
-1996.53835 -6967.44512 -8392.03316 -5995.50174 -8905.26600
391 392 393 394 395
-4288.19146 -7544.40473 -5898.88544 1117.01950 9299.24584
396 397 398 399 400
-8809.09808 -17866.44586 39777.47740 -3095.61057 -62048.45942
401 402 403 404 405
3578.50683 -2810.09095 -7727.95268 -5338.12391 -2569.83924
406 407 408 409 410
3241.94711 504.95854 -8319.08732 -8762.93877 3389.29108
411 412 413 414 415
-1099.97952 -12813.65492 -16680.89304 -576.14468 12346.76878
416 417 418 419 420
6827.22302 -36663.80803 -11707.14838 786.41735 -651.10881
421 422 423 424 425
3165.13327 821.04228 -95507.41492 -14174.47698 -13400.54509
426 427 428 429 430
-13530.65177 -18280.02205 64239.62961 7091.07313 176.52243
431
-57349.52214
> postscript(file="/var/www/html/freestat/rcomp/tmp/6mu9y1291317005.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 13434.20093 NA
1 -22001.95596 13434.20093
2 21226.85096 -22001.95596
3 -160275.42977 21226.85096
4 16397.58555 -160275.42977
5 60646.84247 16397.58555
6 16032.81249 60646.84247
7 126258.78154 16032.81249
8 -56058.28033 126258.78154
9 -37398.80811 -56058.28033
10 163278.40317 -37398.80811
11 -29920.12320 163278.40317
12 -96861.12864 -29920.12320
13 -64519.40595 -96861.12864
14 -20439.27068 -64519.40595
15 75772.13249 -20439.27068
16 211977.61730 75772.13249
17 -695.17654 211977.61730
18 -63823.64434 -695.17654
19 24316.32064 -63823.64434
20 -2609.48780 24316.32064
21 -5432.02860 -2609.48780
22 87991.85379 -5432.02860
23 23994.39941 87991.85379
24 58118.21857 23994.39941
25 -11720.79462 58118.21857
26 -24616.51957 -11720.79462
27 -47827.43884 -24616.51957
28 -39248.44760 -47827.43884
29 -46775.32554 -39248.44760
30 11698.98532 -46775.32554
31 23977.07434 11698.98532
32 2969.53197 23977.07434
33 24800.14668 2969.53197
34 233253.49102 24800.14668
35 19341.62911 233253.49102
36 12721.39819 19341.62911
37 32585.49496 12721.39819
38 112927.13870 32585.49496
39 1926.69030 112927.13870
40 749.00328 1926.69030
41 21741.40118 749.00328
42 9933.07007 21741.40118
43 41693.74516 9933.07007
44 -2878.98424 41693.74516
45 26533.88465 -2878.98424
46 135885.19386 26533.88465
47 39010.69371 135885.19386
48 48207.24353 39010.69371
49 7170.96656 48207.24353
50 26904.27338 7170.96656
51 50230.61351 26904.27338
52 105932.73160 50230.61351
53 61728.44813 105932.73160
54 71347.35780 61728.44813
55 -12546.77449 71347.35780
56 -17586.26788 -12546.77449
57 15168.43483 -17586.26788
58 25201.77626 15168.43483
59 4876.36530 25201.77626
60 16855.50417 4876.36530
61 -9039.99005 16855.50417
62 18848.16067 -9039.99005
63 1718.99167 18848.16067
64 -13662.60753 1718.99167
65 -47869.75964 -13662.60753
66 9014.57767 -47869.75964
67 23548.20767 9014.57767
68 57727.85645 23548.20767
69 22072.41209 57727.85645
70 -50205.84089 22072.41209
71 -32206.04615 -50205.84089
72 37837.79602 -32206.04615
73 12167.08259 37837.79602
74 49597.37408 12167.08259
75 3585.31908 49597.37408
76 44906.14704 3585.31908
77 4503.59482 44906.14704
78 -5544.55776 4503.59482
79 32871.11131 -5544.55776
80 -7470.75206 32871.11131
81 11038.08948 -7470.75206
82 68281.54642 11038.08948
83 -28078.00703 68281.54642
84 35329.60750 -28078.00703
85 -15108.39573 35329.60750
86 16784.11407 -15108.39573
87 36629.15163 16784.11407
88 35190.84211 36629.15163
89 36223.16257 35190.84211
90 25029.00311 36223.16257
91 6626.63789 25029.00311
92 -10407.04109 6626.63789
93 47091.85660 -10407.04109
94 -32454.53520 47091.85660
95 -16206.09370 -32454.53520
96 46179.26939 -16206.09370
97 -16993.81329 46179.26939
98 61619.76125 -16993.81329
99 -411.62527 61619.76125
100 925.71082 -411.62527
101 22420.52417 925.71082
102 35518.64415 22420.52417
103 46129.75950 35518.64415
104 120.79356 46129.75950
105 3516.57683 120.79356
106 -17893.93562 3516.57683
107 23406.29769 -17893.93562
108 29404.16080 23406.29769
109 -2653.77505 29404.16080
110 -669.40498 -2653.77505
111 -28835.81868 -669.40498
112 13269.88357 -28835.81868
113 -135.53936 13269.88357
114 -15959.66616 -135.53936
115 -13066.83572 -15959.66616
116 16832.15408 -13066.83572
117 20240.14593 16832.15408
118 2948.84331 20240.14593
119 -10139.53525 2948.84331
120 34305.02594 -10139.53525
121 -431.86428 34305.02594
122 19068.14355 -431.86428
123 32558.23458 19068.14355
124 -12975.76034 32558.23458
125 754.52169 -12975.76034
126 -5434.63835 754.52169
127 -5896.67232 -5434.63835
128 -10373.64507 -5896.67232
129 -5434.63835 -10373.64507
130 -6113.96409 -5434.63835
131 -15547.64057 -6113.96409
132 -5434.63835 -15547.64057
133 -5497.75708 -5434.63835
134 -5434.63835 -5497.75708
135 -6396.38254 -5434.63835
136 -6185.30969 -6396.38254
137 -5434.63835 -6185.30969
138 -5434.63835 -5434.63835
139 -5434.63835 -5434.63835
140 -8295.31583 -5434.63835
141 -5434.63835 -8295.31583
142 -10792.54403 -5434.63835
143 -5488.64630 -10792.54403
144 -4033.10934 -5488.64630
145 -10509.04011 -4033.10934
146 -5434.63835 -10509.04011
147 -5756.65724 -5434.63835
148 -5434.63835 -5756.65724
149 -5215.24408 -5434.63835
150 -8385.28088 -5215.24408
151 -11533.20262 -8385.28088
152 -5434.63835 -11533.20262
153 -8463.96957 -5434.63835
154 -5439.82166 -8463.96957
155 -5434.63835 -5439.82166
156 -5434.63835 -5434.63835
157 -4838.95761 -5434.63835
158 -5938.45126 -4838.95761
159 -5434.63835 -5938.45126
160 -5434.63835 -5434.63835
161 -5434.63835 -5434.63835
162 -5434.63835 -5434.63835
163 -7424.65483 -5434.63835
164 -6722.13032 -7424.65483
165 -5434.63835 -6722.13032
166 -5434.63835 -5434.63835
167 -22069.20342 -5434.63835
168 -5434.63835 -22069.20342
169 -5434.63835 -5434.63835
170 -5434.63835 -5434.63835
171 -4891.03048 -5434.63835
172 -5434.63835 -4891.03048
173 -4929.30852 -5434.63835
174 -5434.63835 -4929.30852
175 -5809.85033 -5434.63835
176 -5434.63835 -5809.85033
177 -11768.16304 -5434.63835
178 -5434.63835 -11768.16304
179 -5434.63835 -5434.63835
180 2733.20289 -5434.63835
181 -1662.32900 2733.20289
182 -13295.36707 -1662.32900
183 -1355.54708 -13295.36707
184 -10233.22581 -1355.54708
185 -13255.13047 -10233.22581
186 -8915.29276 -13255.13047
187 -4520.48867 -8915.29276
188 -839.70725 -4520.48867
189 -5434.63835 -839.70725
190 -7027.09854 -5434.63835
191 -5434.63835 -7027.09854
192 -5001.46052 -5434.63835
193 -5434.63835 -5001.46052
194 -5417.84981 -5434.63835
195 -11699.36953 -5417.84981
196 34537.14498 -11699.36953
197 -11699.36953 34537.14498
198 -11699.36953 -11699.36953
199 6035.62013 -11699.36953
200 -11699.36953 6035.62013
201 -11699.36953 -11699.36953
202 -11699.36953 -11699.36953
203 -23630.95481 -11699.36953
204 13590.53569 -23630.95481
205 -10911.15327 13590.53569
206 -7039.90324 -10911.15327
207 -11699.36953 -7039.90324
208 -26556.21251 -11699.36953
209 -11699.36953 -26556.21251
210 -5434.63835 -11699.36953
211 -5434.63835 -5434.63835
212 -10911.44305 -5434.63835
213 -5768.25177 -10911.44305
214 -4202.34992 -5768.25177
215 -5434.63835 -4202.34992
216 -5907.23483 -5434.63835
217 -7655.10657 -5907.23483
218 -5434.63835 -7655.10657
219 -11697.25576 -5434.63835
220 -6021.01856 -11697.25576
221 -5434.63835 -6021.01856
222 -5602.31131 -5434.63835
223 -5434.63835 -5602.31131
224 -5500.43749 -5434.63835
225 -6716.45496 -5500.43749
226 -5434.63835 -6716.45496
227 -5434.63835 -5434.63835
228 -1226.07072 -5434.63835
229 -5434.63835 -1226.07072
230 -2986.68408 -5434.63835
231 -5434.63835 -2986.68408
232 -5434.63835 -5434.63835
233 -5810.13601 -5434.63835
234 -5434.63835 -5810.13601
235 -5484.62560 -5434.63835
236 -5434.63835 -5484.62560
237 -6255.92500 -5434.63835
238 -5434.63835 -6255.92500
239 -5434.63835 -5434.63835
240 -5434.63835 -5434.63835
241 -5040.13563 -5434.63835
242 -5434.63835 -5040.13563
243 -6922.94869 -5434.63835
244 5532.63364 -6922.94869
245 -11699.36953 5532.63364
246 -11699.36953 -11699.36953
247 -5182.36599 -11699.36953
248 -11699.36953 -5182.36599
249 -11699.36953 -11699.36953
250 -13758.21314 -11699.36953
251 -11699.36953 -13758.21314
252 -2104.76387 -11699.36953
253 1380.75916 -2104.76387
254 10157.84407 1380.75916
255 -5434.63835 10157.84407
256 -12023.77383 -5434.63835
257 -21092.84495 -12023.77383
258 -11699.36953 -21092.84495
259 -5434.63835 -11699.36953
260 -11699.36953 -5434.63835
261 -11699.36953 -11699.36953
262 5743.60354 -11699.36953
263 -8476.97372 5743.60354
264 -11699.36953 -8476.97372
265 -5434.63835 -11699.36953
266 -11699.36953 -5434.63835
267 10647.44623 -11699.36953
268 -11699.36953 10647.44623
269 3186.44410 -11699.36953
270 -11699.36953 3186.44410
271 -11699.36953 -11699.36953
272 -11699.36953 -11699.36953
273 -6262.94092 -11699.36953
274 -5434.63835 -6262.94092
275 -5959.19781 -5434.63835
276 -4469.33195 -5959.19781
277 -5434.63835 -4469.33195
278 -11160.02371 -5434.63835
279 -98.46622 -11160.02371
280 -11699.36953 -98.46622
281 -5434.63835 -11699.36953
282 -14715.38224 -5434.63835
283 -5434.63835 -14715.38224
284 3614.95458 -5434.63835
285 -5434.63835 3614.95458
286 -11699.36953 -5434.63835
287 -4408.22124 -11699.36953
288 -11699.36953 -4408.22124
289 -16394.54051 -11699.36953
290 -11699.36953 -16394.54051
291 -5100.61259 -11699.36953
292 13313.69640 -5100.61259
293 -5434.63835 13313.69640
294 -7677.99838 -5434.63835
295 -14689.62544 -7677.99838
296 -5434.63835 -14689.62544
297 -11699.36953 -5434.63835
298 -5434.63835 -11699.36953
299 -5434.63835 -5434.63835
300 -11919.95242 -5434.63835
301 -7596.85330 -11919.95242
302 -5434.63835 -7596.85330
303 -11699.36953 -5434.63835
304 -8293.36314 -11699.36953
305 -6071.69898 -8293.36314
306 -11699.36953 -6071.69898
307 -9407.38564 -11699.36953
308 -13179.88151 -9407.38564
309 -11699.36953 -13179.88151
310 -11699.36953 -11699.36953
311 14919.41322 -11699.36953
312 -11699.36953 14919.41322
313 -8906.21837 -11699.36953
314 -5684.07966 -8906.21837
315 -1218.12158 -5684.07966
316 -6248.67180 -1218.12158
317 -6329.13601 -6248.67180
318 -5870.24492 -6329.13601
319 -5434.63835 -5870.24492
320 -5434.63835 -5434.63835
321 -5434.63835 -5434.63835
322 -5434.63835 -5434.63835
323 -11699.36953 -5434.63835
324 -5434.63835 -11699.36953
325 -5584.04057 -5434.63835
326 -11699.36953 -5584.04057
327 -11699.36953 -11699.36953
328 -5434.63835 -11699.36953
329 -5434.63835 -5434.63835
330 -5434.63835 -5434.63835
331 -5434.63835 -5434.63835
332 -5434.63835 -5434.63835
333 -6166.50705 -5434.63835
334 -5572.70823 -6166.50705
335 -5434.63835 -5572.70823
336 -5760.75230 -5434.63835
337 -5434.63835 -5760.75230
338 -5338.86256 -5434.63835
339 -7399.32509 -5338.86256
340 -4363.06910 -7399.32509
341 -5434.63835 -4363.06910
342 -5434.63835 -5434.63835
343 -5434.63835 -5434.63835
344 -5434.63835 -5434.63835
345 -5434.63835 -5434.63835
346 -7506.82664 -5434.63835
347 -5434.63835 -7506.82664
348 -5434.63835 -5434.63835
349 -5434.63835 -5434.63835
350 -8664.09725 -5434.63835
351 -5772.57247 -8664.09725
352 -7605.25241 -5772.57247
353 -4275.11459 -7605.25241
354 -5434.63835 -4275.11459
355 -5434.63835 -5434.63835
356 8706.39311 -5434.63835
357 -5434.63835 8706.39311
358 -5434.63835 -5434.63835
359 -5434.63835 -5434.63835
360 -5722.47658 -5434.63835
361 -5434.63835 -5722.47658
362 -5434.63835 -5434.63835
363 -11054.01137 -5434.63835
364 -7370.14145 -11054.01137
365 -5434.63835 -7370.14145
366 -5434.63835 -5434.63835
367 -5434.63835 -5434.63835
368 -5434.63835 -5434.63835
369 -5434.63835 -5434.63835
370 -5434.63835 -5434.63835
371 -9404.09518 -5434.63835
372 9769.45452 -9404.09518
373 -1459.49048 9769.45452
374 -5434.63835 -1459.49048
375 -5434.63835 -5434.63835
376 -11750.34483 -5434.63835
377 -6697.89903 -11750.34483
378 -5434.63835 -6697.89903
379 -5434.63835 -5434.63835
380 -5434.63835 -5434.63835
381 2775.49509 -5434.63835
382 -30640.55150 2775.49509
383 -5434.63835 -30640.55150
384 -4636.30523 -5434.63835
385 -1996.53835 -4636.30523
386 -6967.44512 -1996.53835
387 -8392.03316 -6967.44512
388 -5995.50174 -8392.03316
389 -8905.26600 -5995.50174
390 -4288.19146 -8905.26600
391 -7544.40473 -4288.19146
392 -5898.88544 -7544.40473
393 1117.01950 -5898.88544
394 9299.24584 1117.01950
395 -8809.09808 9299.24584
396 -17866.44586 -8809.09808
397 39777.47740 -17866.44586
398 -3095.61057 39777.47740
399 -62048.45942 -3095.61057
400 3578.50683 -62048.45942
401 -2810.09095 3578.50683
402 -7727.95268 -2810.09095
403 -5338.12391 -7727.95268
404 -2569.83924 -5338.12391
405 3241.94711 -2569.83924
406 504.95854 3241.94711
407 -8319.08732 504.95854
408 -8762.93877 -8319.08732
409 3389.29108 -8762.93877
410 -1099.97952 3389.29108
411 -12813.65492 -1099.97952
412 -16680.89304 -12813.65492
413 -576.14468 -16680.89304
414 12346.76878 -576.14468
415 6827.22302 12346.76878
416 -36663.80803 6827.22302
417 -11707.14838 -36663.80803
418 786.41735 -11707.14838
419 -651.10881 786.41735
420 3165.13327 -651.10881
421 821.04228 3165.13327
422 -95507.41492 821.04228
423 -14174.47698 -95507.41492
424 -13400.54509 -14174.47698
425 -13530.65177 -13400.54509
426 -18280.02205 -13530.65177
427 64239.62961 -18280.02205
428 7091.07313 64239.62961
429 176.52243 7091.07313
430 -57349.52214 176.52243
431 NA -57349.52214
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -22001.95596 13434.20093
[2,] 21226.85096 -22001.95596
[3,] -160275.42977 21226.85096
[4,] 16397.58555 -160275.42977
[5,] 60646.84247 16397.58555
[6,] 16032.81249 60646.84247
[7,] 126258.78154 16032.81249
[8,] -56058.28033 126258.78154
[9,] -37398.80811 -56058.28033
[10,] 163278.40317 -37398.80811
[11,] -29920.12320 163278.40317
[12,] -96861.12864 -29920.12320
[13,] -64519.40595 -96861.12864
[14,] -20439.27068 -64519.40595
[15,] 75772.13249 -20439.27068
[16,] 211977.61730 75772.13249
[17,] -695.17654 211977.61730
[18,] -63823.64434 -695.17654
[19,] 24316.32064 -63823.64434
[20,] -2609.48780 24316.32064
[21,] -5432.02860 -2609.48780
[22,] 87991.85379 -5432.02860
[23,] 23994.39941 87991.85379
[24,] 58118.21857 23994.39941
[25,] -11720.79462 58118.21857
[26,] -24616.51957 -11720.79462
[27,] -47827.43884 -24616.51957
[28,] -39248.44760 -47827.43884
[29,] -46775.32554 -39248.44760
[30,] 11698.98532 -46775.32554
[31,] 23977.07434 11698.98532
[32,] 2969.53197 23977.07434
[33,] 24800.14668 2969.53197
[34,] 233253.49102 24800.14668
[35,] 19341.62911 233253.49102
[36,] 12721.39819 19341.62911
[37,] 32585.49496 12721.39819
[38,] 112927.13870 32585.49496
[39,] 1926.69030 112927.13870
[40,] 749.00328 1926.69030
[41,] 21741.40118 749.00328
[42,] 9933.07007 21741.40118
[43,] 41693.74516 9933.07007
[44,] -2878.98424 41693.74516
[45,] 26533.88465 -2878.98424
[46,] 135885.19386 26533.88465
[47,] 39010.69371 135885.19386
[48,] 48207.24353 39010.69371
[49,] 7170.96656 48207.24353
[50,] 26904.27338 7170.96656
[51,] 50230.61351 26904.27338
[52,] 105932.73160 50230.61351
[53,] 61728.44813 105932.73160
[54,] 71347.35780 61728.44813
[55,] -12546.77449 71347.35780
[56,] -17586.26788 -12546.77449
[57,] 15168.43483 -17586.26788
[58,] 25201.77626 15168.43483
[59,] 4876.36530 25201.77626
[60,] 16855.50417 4876.36530
[61,] -9039.99005 16855.50417
[62,] 18848.16067 -9039.99005
[63,] 1718.99167 18848.16067
[64,] -13662.60753 1718.99167
[65,] -47869.75964 -13662.60753
[66,] 9014.57767 -47869.75964
[67,] 23548.20767 9014.57767
[68,] 57727.85645 23548.20767
[69,] 22072.41209 57727.85645
[70,] -50205.84089 22072.41209
[71,] -32206.04615 -50205.84089
[72,] 37837.79602 -32206.04615
[73,] 12167.08259 37837.79602
[74,] 49597.37408 12167.08259
[75,] 3585.31908 49597.37408
[76,] 44906.14704 3585.31908
[77,] 4503.59482 44906.14704
[78,] -5544.55776 4503.59482
[79,] 32871.11131 -5544.55776
[80,] -7470.75206 32871.11131
[81,] 11038.08948 -7470.75206
[82,] 68281.54642 11038.08948
[83,] -28078.00703 68281.54642
[84,] 35329.60750 -28078.00703
[85,] -15108.39573 35329.60750
[86,] 16784.11407 -15108.39573
[87,] 36629.15163 16784.11407
[88,] 35190.84211 36629.15163
[89,] 36223.16257 35190.84211
[90,] 25029.00311 36223.16257
[91,] 6626.63789 25029.00311
[92,] -10407.04109 6626.63789
[93,] 47091.85660 -10407.04109
[94,] -32454.53520 47091.85660
[95,] -16206.09370 -32454.53520
[96,] 46179.26939 -16206.09370
[97,] -16993.81329 46179.26939
[98,] 61619.76125 -16993.81329
[99,] -411.62527 61619.76125
[100,] 925.71082 -411.62527
[101,] 22420.52417 925.71082
[102,] 35518.64415 22420.52417
[103,] 46129.75950 35518.64415
[104,] 120.79356 46129.75950
[105,] 3516.57683 120.79356
[106,] -17893.93562 3516.57683
[107,] 23406.29769 -17893.93562
[108,] 29404.16080 23406.29769
[109,] -2653.77505 29404.16080
[110,] -669.40498 -2653.77505
[111,] -28835.81868 -669.40498
[112,] 13269.88357 -28835.81868
[113,] -135.53936 13269.88357
[114,] -15959.66616 -135.53936
[115,] -13066.83572 -15959.66616
[116,] 16832.15408 -13066.83572
[117,] 20240.14593 16832.15408
[118,] 2948.84331 20240.14593
[119,] -10139.53525 2948.84331
[120,] 34305.02594 -10139.53525
[121,] -431.86428 34305.02594
[122,] 19068.14355 -431.86428
[123,] 32558.23458 19068.14355
[124,] -12975.76034 32558.23458
[125,] 754.52169 -12975.76034
[126,] -5434.63835 754.52169
[127,] -5896.67232 -5434.63835
[128,] -10373.64507 -5896.67232
[129,] -5434.63835 -10373.64507
[130,] -6113.96409 -5434.63835
[131,] -15547.64057 -6113.96409
[132,] -5434.63835 -15547.64057
[133,] -5497.75708 -5434.63835
[134,] -5434.63835 -5497.75708
[135,] -6396.38254 -5434.63835
[136,] -6185.30969 -6396.38254
[137,] -5434.63835 -6185.30969
[138,] -5434.63835 -5434.63835
[139,] -5434.63835 -5434.63835
[140,] -8295.31583 -5434.63835
[141,] -5434.63835 -8295.31583
[142,] -10792.54403 -5434.63835
[143,] -5488.64630 -10792.54403
[144,] -4033.10934 -5488.64630
[145,] -10509.04011 -4033.10934
[146,] -5434.63835 -10509.04011
[147,] -5756.65724 -5434.63835
[148,] -5434.63835 -5756.65724
[149,] -5215.24408 -5434.63835
[150,] -8385.28088 -5215.24408
[151,] -11533.20262 -8385.28088
[152,] -5434.63835 -11533.20262
[153,] -8463.96957 -5434.63835
[154,] -5439.82166 -8463.96957
[155,] -5434.63835 -5439.82166
[156,] -5434.63835 -5434.63835
[157,] -4838.95761 -5434.63835
[158,] -5938.45126 -4838.95761
[159,] -5434.63835 -5938.45126
[160,] -5434.63835 -5434.63835
[161,] -5434.63835 -5434.63835
[162,] -5434.63835 -5434.63835
[163,] -7424.65483 -5434.63835
[164,] -6722.13032 -7424.65483
[165,] -5434.63835 -6722.13032
[166,] -5434.63835 -5434.63835
[167,] -22069.20342 -5434.63835
[168,] -5434.63835 -22069.20342
[169,] -5434.63835 -5434.63835
[170,] -5434.63835 -5434.63835
[171,] -4891.03048 -5434.63835
[172,] -5434.63835 -4891.03048
[173,] -4929.30852 -5434.63835
[174,] -5434.63835 -4929.30852
[175,] -5809.85033 -5434.63835
[176,] -5434.63835 -5809.85033
[177,] -11768.16304 -5434.63835
[178,] -5434.63835 -11768.16304
[179,] -5434.63835 -5434.63835
[180,] 2733.20289 -5434.63835
[181,] -1662.32900 2733.20289
[182,] -13295.36707 -1662.32900
[183,] -1355.54708 -13295.36707
[184,] -10233.22581 -1355.54708
[185,] -13255.13047 -10233.22581
[186,] -8915.29276 -13255.13047
[187,] -4520.48867 -8915.29276
[188,] -839.70725 -4520.48867
[189,] -5434.63835 -839.70725
[190,] -7027.09854 -5434.63835
[191,] -5434.63835 -7027.09854
[192,] -5001.46052 -5434.63835
[193,] -5434.63835 -5001.46052
[194,] -5417.84981 -5434.63835
[195,] -11699.36953 -5417.84981
[196,] 34537.14498 -11699.36953
[197,] -11699.36953 34537.14498
[198,] -11699.36953 -11699.36953
[199,] 6035.62013 -11699.36953
[200,] -11699.36953 6035.62013
[201,] -11699.36953 -11699.36953
[202,] -11699.36953 -11699.36953
[203,] -23630.95481 -11699.36953
[204,] 13590.53569 -23630.95481
[205,] -10911.15327 13590.53569
[206,] -7039.90324 -10911.15327
[207,] -11699.36953 -7039.90324
[208,] -26556.21251 -11699.36953
[209,] -11699.36953 -26556.21251
[210,] -5434.63835 -11699.36953
[211,] -5434.63835 -5434.63835
[212,] -10911.44305 -5434.63835
[213,] -5768.25177 -10911.44305
[214,] -4202.34992 -5768.25177
[215,] -5434.63835 -4202.34992
[216,] -5907.23483 -5434.63835
[217,] -7655.10657 -5907.23483
[218,] -5434.63835 -7655.10657
[219,] -11697.25576 -5434.63835
[220,] -6021.01856 -11697.25576
[221,] -5434.63835 -6021.01856
[222,] -5602.31131 -5434.63835
[223,] -5434.63835 -5602.31131
[224,] -5500.43749 -5434.63835
[225,] -6716.45496 -5500.43749
[226,] -5434.63835 -6716.45496
[227,] -5434.63835 -5434.63835
[228,] -1226.07072 -5434.63835
[229,] -5434.63835 -1226.07072
[230,] -2986.68408 -5434.63835
[231,] -5434.63835 -2986.68408
[232,] -5434.63835 -5434.63835
[233,] -5810.13601 -5434.63835
[234,] -5434.63835 -5810.13601
[235,] -5484.62560 -5434.63835
[236,] -5434.63835 -5484.62560
[237,] -6255.92500 -5434.63835
[238,] -5434.63835 -6255.92500
[239,] -5434.63835 -5434.63835
[240,] -5434.63835 -5434.63835
[241,] -5040.13563 -5434.63835
[242,] -5434.63835 -5040.13563
[243,] -6922.94869 -5434.63835
[244,] 5532.63364 -6922.94869
[245,] -11699.36953 5532.63364
[246,] -11699.36953 -11699.36953
[247,] -5182.36599 -11699.36953
[248,] -11699.36953 -5182.36599
[249,] -11699.36953 -11699.36953
[250,] -13758.21314 -11699.36953
[251,] -11699.36953 -13758.21314
[252,] -2104.76387 -11699.36953
[253,] 1380.75916 -2104.76387
[254,] 10157.84407 1380.75916
[255,] -5434.63835 10157.84407
[256,] -12023.77383 -5434.63835
[257,] -21092.84495 -12023.77383
[258,] -11699.36953 -21092.84495
[259,] -5434.63835 -11699.36953
[260,] -11699.36953 -5434.63835
[261,] -11699.36953 -11699.36953
[262,] 5743.60354 -11699.36953
[263,] -8476.97372 5743.60354
[264,] -11699.36953 -8476.97372
[265,] -5434.63835 -11699.36953
[266,] -11699.36953 -5434.63835
[267,] 10647.44623 -11699.36953
[268,] -11699.36953 10647.44623
[269,] 3186.44410 -11699.36953
[270,] -11699.36953 3186.44410
[271,] -11699.36953 -11699.36953
[272,] -11699.36953 -11699.36953
[273,] -6262.94092 -11699.36953
[274,] -5434.63835 -6262.94092
[275,] -5959.19781 -5434.63835
[276,] -4469.33195 -5959.19781
[277,] -5434.63835 -4469.33195
[278,] -11160.02371 -5434.63835
[279,] -98.46622 -11160.02371
[280,] -11699.36953 -98.46622
[281,] -5434.63835 -11699.36953
[282,] -14715.38224 -5434.63835
[283,] -5434.63835 -14715.38224
[284,] 3614.95458 -5434.63835
[285,] -5434.63835 3614.95458
[286,] -11699.36953 -5434.63835
[287,] -4408.22124 -11699.36953
[288,] -11699.36953 -4408.22124
[289,] -16394.54051 -11699.36953
[290,] -11699.36953 -16394.54051
[291,] -5100.61259 -11699.36953
[292,] 13313.69640 -5100.61259
[293,] -5434.63835 13313.69640
[294,] -7677.99838 -5434.63835
[295,] -14689.62544 -7677.99838
[296,] -5434.63835 -14689.62544
[297,] -11699.36953 -5434.63835
[298,] -5434.63835 -11699.36953
[299,] -5434.63835 -5434.63835
[300,] -11919.95242 -5434.63835
[301,] -7596.85330 -11919.95242
[302,] -5434.63835 -7596.85330
[303,] -11699.36953 -5434.63835
[304,] -8293.36314 -11699.36953
[305,] -6071.69898 -8293.36314
[306,] -11699.36953 -6071.69898
[307,] -9407.38564 -11699.36953
[308,] -13179.88151 -9407.38564
[309,] -11699.36953 -13179.88151
[310,] -11699.36953 -11699.36953
[311,] 14919.41322 -11699.36953
[312,] -11699.36953 14919.41322
[313,] -8906.21837 -11699.36953
[314,] -5684.07966 -8906.21837
[315,] -1218.12158 -5684.07966
[316,] -6248.67180 -1218.12158
[317,] -6329.13601 -6248.67180
[318,] -5870.24492 -6329.13601
[319,] -5434.63835 -5870.24492
[320,] -5434.63835 -5434.63835
[321,] -5434.63835 -5434.63835
[322,] -5434.63835 -5434.63835
[323,] -11699.36953 -5434.63835
[324,] -5434.63835 -11699.36953
[325,] -5584.04057 -5434.63835
[326,] -11699.36953 -5584.04057
[327,] -11699.36953 -11699.36953
[328,] -5434.63835 -11699.36953
[329,] -5434.63835 -5434.63835
[330,] -5434.63835 -5434.63835
[331,] -5434.63835 -5434.63835
[332,] -5434.63835 -5434.63835
[333,] -6166.50705 -5434.63835
[334,] -5572.70823 -6166.50705
[335,] -5434.63835 -5572.70823
[336,] -5760.75230 -5434.63835
[337,] -5434.63835 -5760.75230
[338,] -5338.86256 -5434.63835
[339,] -7399.32509 -5338.86256
[340,] -4363.06910 -7399.32509
[341,] -5434.63835 -4363.06910
[342,] -5434.63835 -5434.63835
[343,] -5434.63835 -5434.63835
[344,] -5434.63835 -5434.63835
[345,] -5434.63835 -5434.63835
[346,] -7506.82664 -5434.63835
[347,] -5434.63835 -7506.82664
[348,] -5434.63835 -5434.63835
[349,] -5434.63835 -5434.63835
[350,] -8664.09725 -5434.63835
[351,] -5772.57247 -8664.09725
[352,] -7605.25241 -5772.57247
[353,] -4275.11459 -7605.25241
[354,] -5434.63835 -4275.11459
[355,] -5434.63835 -5434.63835
[356,] 8706.39311 -5434.63835
[357,] -5434.63835 8706.39311
[358,] -5434.63835 -5434.63835
[359,] -5434.63835 -5434.63835
[360,] -5722.47658 -5434.63835
[361,] -5434.63835 -5722.47658
[362,] -5434.63835 -5434.63835
[363,] -11054.01137 -5434.63835
[364,] -7370.14145 -11054.01137
[365,] -5434.63835 -7370.14145
[366,] -5434.63835 -5434.63835
[367,] -5434.63835 -5434.63835
[368,] -5434.63835 -5434.63835
[369,] -5434.63835 -5434.63835
[370,] -5434.63835 -5434.63835
[371,] -9404.09518 -5434.63835
[372,] 9769.45452 -9404.09518
[373,] -1459.49048 9769.45452
[374,] -5434.63835 -1459.49048
[375,] -5434.63835 -5434.63835
[376,] -11750.34483 -5434.63835
[377,] -6697.89903 -11750.34483
[378,] -5434.63835 -6697.89903
[379,] -5434.63835 -5434.63835
[380,] -5434.63835 -5434.63835
[381,] 2775.49509 -5434.63835
[382,] -30640.55150 2775.49509
[383,] -5434.63835 -30640.55150
[384,] -4636.30523 -5434.63835
[385,] -1996.53835 -4636.30523
[386,] -6967.44512 -1996.53835
[387,] -8392.03316 -6967.44512
[388,] -5995.50174 -8392.03316
[389,] -8905.26600 -5995.50174
[390,] -4288.19146 -8905.26600
[391,] -7544.40473 -4288.19146
[392,] -5898.88544 -7544.40473
[393,] 1117.01950 -5898.88544
[394,] 9299.24584 1117.01950
[395,] -8809.09808 9299.24584
[396,] -17866.44586 -8809.09808
[397,] 39777.47740 -17866.44586
[398,] -3095.61057 39777.47740
[399,] -62048.45942 -3095.61057
[400,] 3578.50683 -62048.45942
[401,] -2810.09095 3578.50683
[402,] -7727.95268 -2810.09095
[403,] -5338.12391 -7727.95268
[404,] -2569.83924 -5338.12391
[405,] 3241.94711 -2569.83924
[406,] 504.95854 3241.94711
[407,] -8319.08732 504.95854
[408,] -8762.93877 -8319.08732
[409,] 3389.29108 -8762.93877
[410,] -1099.97952 3389.29108
[411,] -12813.65492 -1099.97952
[412,] -16680.89304 -12813.65492
[413,] -576.14468 -16680.89304
[414,] 12346.76878 -576.14468
[415,] 6827.22302 12346.76878
[416,] -36663.80803 6827.22302
[417,] -11707.14838 -36663.80803
[418,] 786.41735 -11707.14838
[419,] -651.10881 786.41735
[420,] 3165.13327 -651.10881
[421,] 821.04228 3165.13327
[422,] -95507.41492 821.04228
[423,] -14174.47698 -95507.41492
[424,] -13400.54509 -14174.47698
[425,] -13530.65177 -13400.54509
[426,] -18280.02205 -13530.65177
[427,] 64239.62961 -18280.02205
[428,] 7091.07313 64239.62961
[429,] 176.52243 7091.07313
[430,] -57349.52214 176.52243
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -22001.95596 13434.20093
2 21226.85096 -22001.95596
3 -160275.42977 21226.85096
4 16397.58555 -160275.42977
5 60646.84247 16397.58555
6 16032.81249 60646.84247
7 126258.78154 16032.81249
8 -56058.28033 126258.78154
9 -37398.80811 -56058.28033
10 163278.40317 -37398.80811
11 -29920.12320 163278.40317
12 -96861.12864 -29920.12320
13 -64519.40595 -96861.12864
14 -20439.27068 -64519.40595
15 75772.13249 -20439.27068
16 211977.61730 75772.13249
17 -695.17654 211977.61730
18 -63823.64434 -695.17654
19 24316.32064 -63823.64434
20 -2609.48780 24316.32064
21 -5432.02860 -2609.48780
22 87991.85379 -5432.02860
23 23994.39941 87991.85379
24 58118.21857 23994.39941
25 -11720.79462 58118.21857
26 -24616.51957 -11720.79462
27 -47827.43884 -24616.51957
28 -39248.44760 -47827.43884
29 -46775.32554 -39248.44760
30 11698.98532 -46775.32554
31 23977.07434 11698.98532
32 2969.53197 23977.07434
33 24800.14668 2969.53197
34 233253.49102 24800.14668
35 19341.62911 233253.49102
36 12721.39819 19341.62911
37 32585.49496 12721.39819
38 112927.13870 32585.49496
39 1926.69030 112927.13870
40 749.00328 1926.69030
41 21741.40118 749.00328
42 9933.07007 21741.40118
43 41693.74516 9933.07007
44 -2878.98424 41693.74516
45 26533.88465 -2878.98424
46 135885.19386 26533.88465
47 39010.69371 135885.19386
48 48207.24353 39010.69371
49 7170.96656 48207.24353
50 26904.27338 7170.96656
51 50230.61351 26904.27338
52 105932.73160 50230.61351
53 61728.44813 105932.73160
54 71347.35780 61728.44813
55 -12546.77449 71347.35780
56 -17586.26788 -12546.77449
57 15168.43483 -17586.26788
58 25201.77626 15168.43483
59 4876.36530 25201.77626
60 16855.50417 4876.36530
61 -9039.99005 16855.50417
62 18848.16067 -9039.99005
63 1718.99167 18848.16067
64 -13662.60753 1718.99167
65 -47869.75964 -13662.60753
66 9014.57767 -47869.75964
67 23548.20767 9014.57767
68 57727.85645 23548.20767
69 22072.41209 57727.85645
70 -50205.84089 22072.41209
71 -32206.04615 -50205.84089
72 37837.79602 -32206.04615
73 12167.08259 37837.79602
74 49597.37408 12167.08259
75 3585.31908 49597.37408
76 44906.14704 3585.31908
77 4503.59482 44906.14704
78 -5544.55776 4503.59482
79 32871.11131 -5544.55776
80 -7470.75206 32871.11131
81 11038.08948 -7470.75206
82 68281.54642 11038.08948
83 -28078.00703 68281.54642
84 35329.60750 -28078.00703
85 -15108.39573 35329.60750
86 16784.11407 -15108.39573
87 36629.15163 16784.11407
88 35190.84211 36629.15163
89 36223.16257 35190.84211
90 25029.00311 36223.16257
91 6626.63789 25029.00311
92 -10407.04109 6626.63789
93 47091.85660 -10407.04109
94 -32454.53520 47091.85660
95 -16206.09370 -32454.53520
96 46179.26939 -16206.09370
97 -16993.81329 46179.26939
98 61619.76125 -16993.81329
99 -411.62527 61619.76125
100 925.71082 -411.62527
101 22420.52417 925.71082
102 35518.64415 22420.52417
103 46129.75950 35518.64415
104 120.79356 46129.75950
105 3516.57683 120.79356
106 -17893.93562 3516.57683
107 23406.29769 -17893.93562
108 29404.16080 23406.29769
109 -2653.77505 29404.16080
110 -669.40498 -2653.77505
111 -28835.81868 -669.40498
112 13269.88357 -28835.81868
113 -135.53936 13269.88357
114 -15959.66616 -135.53936
115 -13066.83572 -15959.66616
116 16832.15408 -13066.83572
117 20240.14593 16832.15408
118 2948.84331 20240.14593
119 -10139.53525 2948.84331
120 34305.02594 -10139.53525
121 -431.86428 34305.02594
122 19068.14355 -431.86428
123 32558.23458 19068.14355
124 -12975.76034 32558.23458
125 754.52169 -12975.76034
126 -5434.63835 754.52169
127 -5896.67232 -5434.63835
128 -10373.64507 -5896.67232
129 -5434.63835 -10373.64507
130 -6113.96409 -5434.63835
131 -15547.64057 -6113.96409
132 -5434.63835 -15547.64057
133 -5497.75708 -5434.63835
134 -5434.63835 -5497.75708
135 -6396.38254 -5434.63835
136 -6185.30969 -6396.38254
137 -5434.63835 -6185.30969
138 -5434.63835 -5434.63835
139 -5434.63835 -5434.63835
140 -8295.31583 -5434.63835
141 -5434.63835 -8295.31583
142 -10792.54403 -5434.63835
143 -5488.64630 -10792.54403
144 -4033.10934 -5488.64630
145 -10509.04011 -4033.10934
146 -5434.63835 -10509.04011
147 -5756.65724 -5434.63835
148 -5434.63835 -5756.65724
149 -5215.24408 -5434.63835
150 -8385.28088 -5215.24408
151 -11533.20262 -8385.28088
152 -5434.63835 -11533.20262
153 -8463.96957 -5434.63835
154 -5439.82166 -8463.96957
155 -5434.63835 -5439.82166
156 -5434.63835 -5434.63835
157 -4838.95761 -5434.63835
158 -5938.45126 -4838.95761
159 -5434.63835 -5938.45126
160 -5434.63835 -5434.63835
161 -5434.63835 -5434.63835
162 -5434.63835 -5434.63835
163 -7424.65483 -5434.63835
164 -6722.13032 -7424.65483
165 -5434.63835 -6722.13032
166 -5434.63835 -5434.63835
167 -22069.20342 -5434.63835
168 -5434.63835 -22069.20342
169 -5434.63835 -5434.63835
170 -5434.63835 -5434.63835
171 -4891.03048 -5434.63835
172 -5434.63835 -4891.03048
173 -4929.30852 -5434.63835
174 -5434.63835 -4929.30852
175 -5809.85033 -5434.63835
176 -5434.63835 -5809.85033
177 -11768.16304 -5434.63835
178 -5434.63835 -11768.16304
179 -5434.63835 -5434.63835
180 2733.20289 -5434.63835
181 -1662.32900 2733.20289
182 -13295.36707 -1662.32900
183 -1355.54708 -13295.36707
184 -10233.22581 -1355.54708
185 -13255.13047 -10233.22581
186 -8915.29276 -13255.13047
187 -4520.48867 -8915.29276
188 -839.70725 -4520.48867
189 -5434.63835 -839.70725
190 -7027.09854 -5434.63835
191 -5434.63835 -7027.09854
192 -5001.46052 -5434.63835
193 -5434.63835 -5001.46052
194 -5417.84981 -5434.63835
195 -11699.36953 -5417.84981
196 34537.14498 -11699.36953
197 -11699.36953 34537.14498
198 -11699.36953 -11699.36953
199 6035.62013 -11699.36953
200 -11699.36953 6035.62013
201 -11699.36953 -11699.36953
202 -11699.36953 -11699.36953
203 -23630.95481 -11699.36953
204 13590.53569 -23630.95481
205 -10911.15327 13590.53569
206 -7039.90324 -10911.15327
207 -11699.36953 -7039.90324
208 -26556.21251 -11699.36953
209 -11699.36953 -26556.21251
210 -5434.63835 -11699.36953
211 -5434.63835 -5434.63835
212 -10911.44305 -5434.63835
213 -5768.25177 -10911.44305
214 -4202.34992 -5768.25177
215 -5434.63835 -4202.34992
216 -5907.23483 -5434.63835
217 -7655.10657 -5907.23483
218 -5434.63835 -7655.10657
219 -11697.25576 -5434.63835
220 -6021.01856 -11697.25576
221 -5434.63835 -6021.01856
222 -5602.31131 -5434.63835
223 -5434.63835 -5602.31131
224 -5500.43749 -5434.63835
225 -6716.45496 -5500.43749
226 -5434.63835 -6716.45496
227 -5434.63835 -5434.63835
228 -1226.07072 -5434.63835
229 -5434.63835 -1226.07072
230 -2986.68408 -5434.63835
231 -5434.63835 -2986.68408
232 -5434.63835 -5434.63835
233 -5810.13601 -5434.63835
234 -5434.63835 -5810.13601
235 -5484.62560 -5434.63835
236 -5434.63835 -5484.62560
237 -6255.92500 -5434.63835
238 -5434.63835 -6255.92500
239 -5434.63835 -5434.63835
240 -5434.63835 -5434.63835
241 -5040.13563 -5434.63835
242 -5434.63835 -5040.13563
243 -6922.94869 -5434.63835
244 5532.63364 -6922.94869
245 -11699.36953 5532.63364
246 -11699.36953 -11699.36953
247 -5182.36599 -11699.36953
248 -11699.36953 -5182.36599
249 -11699.36953 -11699.36953
250 -13758.21314 -11699.36953
251 -11699.36953 -13758.21314
252 -2104.76387 -11699.36953
253 1380.75916 -2104.76387
254 10157.84407 1380.75916
255 -5434.63835 10157.84407
256 -12023.77383 -5434.63835
257 -21092.84495 -12023.77383
258 -11699.36953 -21092.84495
259 -5434.63835 -11699.36953
260 -11699.36953 -5434.63835
261 -11699.36953 -11699.36953
262 5743.60354 -11699.36953
263 -8476.97372 5743.60354
264 -11699.36953 -8476.97372
265 -5434.63835 -11699.36953
266 -11699.36953 -5434.63835
267 10647.44623 -11699.36953
268 -11699.36953 10647.44623
269 3186.44410 -11699.36953
270 -11699.36953 3186.44410
271 -11699.36953 -11699.36953
272 -11699.36953 -11699.36953
273 -6262.94092 -11699.36953
274 -5434.63835 -6262.94092
275 -5959.19781 -5434.63835
276 -4469.33195 -5959.19781
277 -5434.63835 -4469.33195
278 -11160.02371 -5434.63835
279 -98.46622 -11160.02371
280 -11699.36953 -98.46622
281 -5434.63835 -11699.36953
282 -14715.38224 -5434.63835
283 -5434.63835 -14715.38224
284 3614.95458 -5434.63835
285 -5434.63835 3614.95458
286 -11699.36953 -5434.63835
287 -4408.22124 -11699.36953
288 -11699.36953 -4408.22124
289 -16394.54051 -11699.36953
290 -11699.36953 -16394.54051
291 -5100.61259 -11699.36953
292 13313.69640 -5100.61259
293 -5434.63835 13313.69640
294 -7677.99838 -5434.63835
295 -14689.62544 -7677.99838
296 -5434.63835 -14689.62544
297 -11699.36953 -5434.63835
298 -5434.63835 -11699.36953
299 -5434.63835 -5434.63835
300 -11919.95242 -5434.63835
301 -7596.85330 -11919.95242
302 -5434.63835 -7596.85330
303 -11699.36953 -5434.63835
304 -8293.36314 -11699.36953
305 -6071.69898 -8293.36314
306 -11699.36953 -6071.69898
307 -9407.38564 -11699.36953
308 -13179.88151 -9407.38564
309 -11699.36953 -13179.88151
310 -11699.36953 -11699.36953
311 14919.41322 -11699.36953
312 -11699.36953 14919.41322
313 -8906.21837 -11699.36953
314 -5684.07966 -8906.21837
315 -1218.12158 -5684.07966
316 -6248.67180 -1218.12158
317 -6329.13601 -6248.67180
318 -5870.24492 -6329.13601
319 -5434.63835 -5870.24492
320 -5434.63835 -5434.63835
321 -5434.63835 -5434.63835
322 -5434.63835 -5434.63835
323 -11699.36953 -5434.63835
324 -5434.63835 -11699.36953
325 -5584.04057 -5434.63835
326 -11699.36953 -5584.04057
327 -11699.36953 -11699.36953
328 -5434.63835 -11699.36953
329 -5434.63835 -5434.63835
330 -5434.63835 -5434.63835
331 -5434.63835 -5434.63835
332 -5434.63835 -5434.63835
333 -6166.50705 -5434.63835
334 -5572.70823 -6166.50705
335 -5434.63835 -5572.70823
336 -5760.75230 -5434.63835
337 -5434.63835 -5760.75230
338 -5338.86256 -5434.63835
339 -7399.32509 -5338.86256
340 -4363.06910 -7399.32509
341 -5434.63835 -4363.06910
342 -5434.63835 -5434.63835
343 -5434.63835 -5434.63835
344 -5434.63835 -5434.63835
345 -5434.63835 -5434.63835
346 -7506.82664 -5434.63835
347 -5434.63835 -7506.82664
348 -5434.63835 -5434.63835
349 -5434.63835 -5434.63835
350 -8664.09725 -5434.63835
351 -5772.57247 -8664.09725
352 -7605.25241 -5772.57247
353 -4275.11459 -7605.25241
354 -5434.63835 -4275.11459
355 -5434.63835 -5434.63835
356 8706.39311 -5434.63835
357 -5434.63835 8706.39311
358 -5434.63835 -5434.63835
359 -5434.63835 -5434.63835
360 -5722.47658 -5434.63835
361 -5434.63835 -5722.47658
362 -5434.63835 -5434.63835
363 -11054.01137 -5434.63835
364 -7370.14145 -11054.01137
365 -5434.63835 -7370.14145
366 -5434.63835 -5434.63835
367 -5434.63835 -5434.63835
368 -5434.63835 -5434.63835
369 -5434.63835 -5434.63835
370 -5434.63835 -5434.63835
371 -9404.09518 -5434.63835
372 9769.45452 -9404.09518
373 -1459.49048 9769.45452
374 -5434.63835 -1459.49048
375 -5434.63835 -5434.63835
376 -11750.34483 -5434.63835
377 -6697.89903 -11750.34483
378 -5434.63835 -6697.89903
379 -5434.63835 -5434.63835
380 -5434.63835 -5434.63835
381 2775.49509 -5434.63835
382 -30640.55150 2775.49509
383 -5434.63835 -30640.55150
384 -4636.30523 -5434.63835
385 -1996.53835 -4636.30523
386 -6967.44512 -1996.53835
387 -8392.03316 -6967.44512
388 -5995.50174 -8392.03316
389 -8905.26600 -5995.50174
390 -4288.19146 -8905.26600
391 -7544.40473 -4288.19146
392 -5898.88544 -7544.40473
393 1117.01950 -5898.88544
394 9299.24584 1117.01950
395 -8809.09808 9299.24584
396 -17866.44586 -8809.09808
397 39777.47740 -17866.44586
398 -3095.61057 39777.47740
399 -62048.45942 -3095.61057
400 3578.50683 -62048.45942
401 -2810.09095 3578.50683
402 -7727.95268 -2810.09095
403 -5338.12391 -7727.95268
404 -2569.83924 -5338.12391
405 3241.94711 -2569.83924
406 504.95854 3241.94711
407 -8319.08732 504.95854
408 -8762.93877 -8319.08732
409 3389.29108 -8762.93877
410 -1099.97952 3389.29108
411 -12813.65492 -1099.97952
412 -16680.89304 -12813.65492
413 -576.14468 -16680.89304
414 12346.76878 -576.14468
415 6827.22302 12346.76878
416 -36663.80803 6827.22302
417 -11707.14838 -36663.80803
418 786.41735 -11707.14838
419 -651.10881 786.41735
420 3165.13327 -651.10881
421 821.04228 3165.13327
422 -95507.41492 821.04228
423 -14174.47698 -95507.41492
424 -13400.54509 -14174.47698
425 -13530.65177 -13400.54509
426 -18280.02205 -13530.65177
427 64239.62961 -18280.02205
428 7091.07313 64239.62961
429 176.52243 7091.07313
430 -57349.52214 176.52243
> 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/7wl811291317005.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/8wl811291317005.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/9pc7m1291317005.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/10pc7m1291317005.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/11ad6s1291317005.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/12wwny1291317005.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/13kwj91291317005.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/14d61c1291317005.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/15h6zi1291317005.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/16dgxr1291317005.tab")
+ }
>
> try(system("convert tmp/1ibba1291317005.ps tmp/1ibba1291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ibba1291317005.ps tmp/2ibba1291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/3tlav1291317005.ps tmp/3tlav1291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/4tlav1291317005.ps tmp/4tlav1291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/5tlav1291317005.ps tmp/5tlav1291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/6mu9y1291317005.ps tmp/6mu9y1291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/7wl811291317005.ps tmp/7wl811291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/8wl811291317005.ps tmp/8wl811291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/9pc7m1291317005.ps tmp/9pc7m1291317005.png",intern=TRUE))
character(0)
> try(system("convert tmp/10pc7m1291317005.ps tmp/10pc7m1291317005.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.454 3.127 14.278