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 = '3'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Trades Group Costs Dividends
1 1081 1 162556 213118
2 309 1 29790 81767
3 458 1 87550 153198
4 588 0 84738 -26007
5 299 1 54660 126942
6 156 1 42634 157214
7 481 0 40949 129352
8 323 1 42312 234817
9 452 1 37704 60448
10 109 1 16275 47818
11 115 0 25830 245546
12 110 0 12679 48020
13 239 1 18014 -1710
14 247 0 43556 32648
15 497 1 24524 95350
16 103 0 6532 151352
17 109 0 7123 288170
18 502 1 20813 114337
19 248 1 37597 37884
20 373 0 17821 122844
21 119 1 12988 82340
22 84 1 22330 79801
23 102 0 13326 165548
24 295 0 16189 116384
25 105 0 7146 134028
26 64 0 15824 63838
27 267 1 26088 74996
28 129 0 11326 31080
29 37 0 8568 32168
30 361 0 14416 49857
31 28 1 3369 87161
32 85 1 11819 106113
33 44 1 6620 80570
34 49 1 4519 102129
35 22 0 2220 301670
36 155 0 18562 102313
37 91 0 10327 88577
38 81 1 5336 112477
39 79 1 2365 191778
40 145 0 4069 79804
41 816 0 7710 128294
42 61 0 13718 96448
43 226 0 4525 93811
44 105 0 6869 117520
45 62 0 4628 69159
46 24 1 3653 101792
47 26 1 1265 210568
48 322 1 7489 136996
49 84 0 4901 121920
50 33 0 2284 76403
51 108 1 3160 108094
52 150 1 4150 134759
53 115 1 7285 188873
54 162 1 1134 146216
55 158 1 4658 156608
56 97 0 2384 61348
57 9 0 3748 50350
58 66 0 5371 87720
59 107 0 1285 99489
60 101 1 9327 87419
61 47 1 5565 94355
62 38 0 1528 60326
63 34 1 3122 94670
64 84 1 7317 82425
65 79 0 2675 59017
66 947 0 13253 90829
67 74 0 880 80791
68 53 1 2053 100423
69 94 0 1424 131116
70 63 1 4036 100269
71 58 1 3045 27330
72 49 0 5119 39039
73 34 0 1431 106885
74 11 0 554 79285
75 35 0 1975 118881
76 17 1 1286 77623
77 47 0 1012 114768
78 43 0 810 74015
79 117 0 1280 69465
80 171 1 666 117869
81 26 0 1380 60982
82 73 1 4608 90131
83 59 0 876 138971
84 18 0 814 39625
85 15 0 514 102725
86 72 1 5692 64239
87 86 0 3642 90262
88 14 0 540 103960
89 64 0 2099 106611
90 11 0 567 103345
91 52 0 2001 95551
92 41 1 2949 82903
93 99 0 2253 63593
94 75 1 6533 126910
95 45 0 1889 37527
96 43 1 3055 60247
97 8 0 272 112995
98 198 1 1414 70184
99 22 0 2564 130140
100 11 1 1383 73221
101 33 1 1261 76114
102 23 0 975 90534
103 80 0 3366 108479
104 18 0 576 113761
105 28 0 1306 68696
106 23 0 746 71561
107 60 1 3192 59831
108 20 1 2045 97890
109 59 0 5477 101481
110 36 1 1932 72954
111 30 0 936 67939
112 47 1 3437 48022
113 71 0 5131 86111
114 14 1 2397 74020
115 9 1 1389 57530
116 39 0 1503 56364
117 26 0 402 84990
118 21 0 2239 88590
119 16 1 2234 77200
120 69 0 837 61262
121 92 0 10579 110309
122 14 0 875 67000
123 103 0 1395 93099
124 29 1 1659 107577
125 37 1 2647 62920
126 23 1 3294 75832
127 0 0 0 60720
128 7 0 94 60793
129 28 0 422 57935
130 0 0 0 60720
131 8 0 34 60630
132 63 0 1558 55637
133 0 0 0 60720
134 3 0 43 60887
135 0 0 0 60720
136 9 0 316 60505
137 13 0 115 60945
138 0 0 0 60720
139 0 0 0 60720
140 0 0 0 60720
141 14 0 389 58990
142 0 0 0 60720
143 15 0 1002 56750
144 3 0 36 60894
145 15 0 460 63346
146 11 0 309 56535
147 0 0 0 60720
148 6 0 9 60835
149 0 0 0 60720
150 1 0 14 61016
151 10 0 520 58650
152 73 0 1766 60438
153 0 0 0 60720
154 11 0 458 58625
155 3 0 20 60938
156 0 0 0 60720
157 0 0 0 60720
158 2 0 98 61490
159 7 0 405 60845
160 0 0 0 60720
161 0 0 0 60720
162 0 0 0 60720
163 0 0 0 60720
164 27 0 483 60830
165 51 0 454 63261
166 0 0 0 60720
167 0 0 0 60720
168 19 0 757 45689
169 0 0 0 60720
170 0 0 0 60720
171 0 0 0 60720
172 4 0 36 61564
173 0 0 0 60720
174 9 0 203 61938
175 0 0 0 60720
176 8 0 90 60951
177 0 0 0 60720
178 1 1 71 60745
179 0 0 0 60720
180 0 0 0 60720
181 34 0 972 71642
182 10 1 531 71641
183 38 0 604 55792
184 10 1 283 71873
185 5 1 23 62555
186 14 1 638 60370
187 16 1 699 64873
188 5 0 149 62041
189 5 0 226 65745
190 0 0 0 60720
191 4 0 275 59500
192 0 0 0 60720
193 6 0 141 61630
194 0 0 0 60720
195 2 0 28 60890
196 0 1 0 60720
197 80 1 2566 113521
198 0 1 0 60720
199 0 1 0 60720
200 20 1 472 80045
201 0 1 0 60720
202 0 1 0 60720
203 0 1 0 60720
204 27 1 203 50804
205 17 1 496 87390
206 2 1 10 61656
207 4 1 63 65688
208 0 1 0 60720
209 32 1 1136 48522
210 0 1 0 60720
211 0 0 0 60720
212 0 0 0 60720
213 32 0 267 57640
214 20 0 474 61977
215 7 0 534 62620
216 0 0 0 60720
217 8 0 15 60831
218 28 0 397 60646
219 0 0 0 60720
220 20 0 1061 56225
221 4 0 288 60510
222 0 0 0 60720
223 2 0 3 60698
224 0 0 0 60720
225 2 0 20 60805
226 26 0 278 61404
227 0 0 0 60720
228 0 0 0 60720
229 4 0 192 65276
230 0 0 0 60720
231 9 0 317 63915
232 0 0 0 60720
233 0 0 0 60720
234 17 0 368 61686
235 0 0 0 60720
236 1 0 2 60743
237 0 0 0 60720
238 6 0 53 60349
239 0 0 0 60720
240 0 0 0 60720
241 0 0 0 60720
242 3 0 94 61360
243 0 0 0 60720
244 8 0 24 59818
245 4 0 2332 72680
246 0 1 0 60720
247 0 1 0 60720
248 11 0 131 61808
249 0 1 0 60720
250 0 1 0 60720
251 9 0 206 53110
252 0 1 0 60720
253 2 0 167 64245
254 73 0 622 73007
255 85 0 885 82732
256 0 0 0 60720
257 8 0 365 54820
258 35 0 364 47705
259 0 1 0 60720
260 0 0 0 60720
261 0 1 0 60720
262 0 1 0 60720
263 12 0 226 72835
264 15 0 307 58856
265 0 1 0 60720
266 0 0 0 60720
267 0 1 0 60720
268 11 0 188 77655
269 0 1 0 60720
270 6 0 138 69817
271 0 1 0 60720
272 0 1 0 60720
273 0 1 0 60720
274 12 0 125 60798
275 0 0 0 60720
276 30 0 282 62452
277 33 0 335 64175
278 0 0 0 60720
279 82 1 813 67440
280 28 0 176 68136
281 0 1 0 60720
282 0 0 0 60720
283 72 0 249 56726
284 0 0 0 60720
285 13 0 333 70811
286 0 0 0 60720
287 0 1 0 60720
288 4 0 30 62045
289 0 1 0 60720
290 62 0 249 54323
291 0 1 0 60720
292 24 0 165 62841
293 21 0 453 81125
294 0 0 0 60720
295 14 0 53 59506
296 21 1 382 59365
297 0 0 0 60720
298 0 1 0 60720
299 0 0 0 60720
300 0 0 0 60720
301 4 1 30 60798
302 2 0 290 58790
303 0 0 0 60720
304 0 1 0 60720
305 53 0 366 61808
306 9 0 2 60735
307 0 1 0 60720
308 13 1 209 64016
309 22 0 384 54683
310 0 1 0 60720
311 0 1 0 60720
312 83 0 365 87192
313 0 1 0 60720
314 8 1 49 64107
315 4 0 3 60761
316 14 0 133 65990
317 1 0 32 59988
318 17 0 368 61167
319 6 0 1 60719
320 0 0 0 60720
321 0 0 0 60720
322 0 0 0 60720
323 0 0 0 60720
324 0 1 0 60720
325 0 0 0 60720
326 2 0 22 60722
327 0 1 0 60720
328 0 1 0 60720
329 0 0 0 60720
330 0 0 0 60720
331 0 0 0 60720
332 0 0 0 60720
333 0 0 0 60720
334 5 0 96 60379
335 2 0 1 60727
336 0 0 0 60720
337 7 0 81 60925
338 0 0 0 60720
339 1 0 26 60896
340 13 0 125 59734
341 15 0 304 62969
342 0 0 0 60720
343 0 0 0 60720
344 0 0 0 60720
345 0 0 0 60720
346 0 0 0 60720
347 6 0 119 59118
348 0 0 0 60720
349 0 0 0 60720
350 0 0 0 60720
351 14 0 312 58598
352 10 0 60 61124
353 12 0 587 59595
354 2 0 135 62065
355 0 0 0 60720
356 0 0 0 60720
357 52 0 514 78780
358 0 0 0 60720
359 0 0 0 60720
360 0 0 0 60720
361 4 0 1 60722
362 0 0 0 60720
363 0 0 0 60720
364 3 1 58 61600
365 11 0 180 59635
366 0 0 0 60720
367 0 0 0 60720
368 0 0 0 60720
369 0 0 0 60720
370 0 0 0 60720
371 0 0 0 60720
372 40 0 448 59781
373 9 0 227 76644
374 1 0 174 64820
375 0 0 0 60720
376 0 0 0 60720
377 24 0 121 56178
378 11 0 607 60436
379 0 0 0 60720
380 0 0 0 60720
381 0 0 0 60720
382 60 0 530 73433
383 80 0 571 41477
384 0 0 0 60720
385 16 0 78 62700
386 40 0 2489 67804
387 6 0 131 59661
388 8 0 923 58620
389 3 0 72 60398
390 16 0 572 58580
391 10 0 397 62710
392 8 0 450 59325
393 7 0 622 60950
394 8 0 694 68060
395 12 1 3425 83620
396 13 0 562 58456
397 42 0 4917 52811
398 118 1 1442 121173
399 9 0 529 63870
400 138 1 2126 21001
401 5 0 1061 70415
402 9 0 776 64230
403 8 0 611 59190
404 25 1 1526 69351
405 7 0 592 64270
406 13 0 1182 70694
407 16 0 621 68005
408 11 0 989 58930
409 11 0 438 58320
410 3 0 726 69980
411 61 0 1303 69863
412 24 1 6341 63255
413 17 1 1164 57320
414 33 1 3310 75230
415 7 0 1366 79420
416 3 0 965 73490
417 66 0 3256 35250
418 17 1 1135 62285
419 26 0 1270 69206
420 3 0 661 65920
421 2 0 1013 69770
422 67 0 2844 72683
423 70 1 11528 -14545
424 26 0 6526 55830
425 24 0 2264 55174
426 94 1 4461 67038
427 30 0 3999 51252
428 223 0 35624 157278
429 48 0 9252 79510
430 90 0 15236 77440
431 180 0 18073 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Dividends
-6.395068 1.946248 0.006791 0.000398
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-193.048 -17.772 -15.941 -1.327 827.248
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -6.3950677 8.5726889 -0.746 0.456089
Group 1.9462480 7.7200912 0.252 0.801084
Costs 0.0067907 0.0003064 22.165 < 2e-16 ***
Dividends 0.0003980 0.0001118 3.559 0.000414 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 70.29 on 427 degrees of freedom
Multiple R-squared: 0.5965, Adjusted R-squared: 0.5937
F-statistic: 210.4 on 3 and 427 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 0.9891982 2.160369e-02 1.080184e-02
[2,] 0.9774971 4.500589e-02 2.250295e-02
[3,] 0.9990876 1.824775e-03 9.123873e-04
[4,] 0.9983066 3.386856e-03 1.693428e-03
[5,] 0.9993478 1.304490e-03 6.522451e-04
[6,] 0.9990296 1.940748e-03 9.703741e-04
[7,] 0.9984426 3.114777e-03 1.557389e-03
[8,] 0.9990658 1.868425e-03 9.342124e-04
[9,] 0.9999981 3.812962e-06 1.906481e-06
[10,] 0.9999959 8.172160e-06 4.086080e-06
[11,] 0.9999937 1.258336e-05 6.291678e-06
[12,] 1.0000000 2.387285e-08 1.193643e-08
[13,] 1.0000000 1.585370e-08 7.926851e-09
[14,] 1.0000000 1.309150e-09 6.545748e-10
[15,] 1.0000000 1.590912e-09 7.954561e-10
[16,] 1.0000000 2.375140e-10 1.187570e-10
[17,] 1.0000000 3.129417e-10 1.564709e-10
[18,] 1.0000000 2.416363e-10 1.208181e-10
[19,] 1.0000000 4.868240e-10 2.434120e-10
[20,] 1.0000000 1.981510e-10 9.907549e-11
[21,] 1.0000000 4.117460e-10 2.058730e-10
[22,] 1.0000000 8.181911e-10 4.090955e-10
[23,] 1.0000000 7.088641e-10 3.544320e-10
[24,] 1.0000000 3.464968e-11 1.732484e-11
[25,] 1.0000000 3.613769e-11 1.806884e-11
[26,] 1.0000000 4.301578e-11 2.150789e-11
[27,] 1.0000000 5.092505e-11 2.546253e-11
[28,] 1.0000000 8.195431e-11 4.097716e-11
[29,] 1.0000000 7.559975e-11 3.779988e-11
[30,] 1.0000000 1.122401e-10 5.612005e-11
[31,] 1.0000000 1.695180e-10 8.475901e-11
[32,] 1.0000000 3.227800e-10 1.613900e-10
[33,] 1.0000000 5.922317e-10 2.961158e-10
[34,] 1.0000000 8.910480e-10 4.455240e-10
[35,] 1.0000000 5.478572e-37 2.739286e-37
[36,] 1.0000000 1.051627e-37 5.258134e-38
[37,] 1.0000000 2.099185e-38 1.049592e-38
[38,] 1.0000000 5.455465e-38 2.727733e-38
[39,] 1.0000000 9.716315e-38 4.858158e-38
[40,] 1.0000000 1.520509e-37 7.602546e-38
[41,] 1.0000000 2.076130e-37 1.038065e-37
[42,] 1.0000000 4.025627e-40 2.012814e-40
[43,] 1.0000000 9.997115e-40 4.998557e-40
[44,] 1.0000000 1.559454e-39 7.797270e-40
[45,] 1.0000000 4.066155e-39 2.033078e-39
[46,] 1.0000000 7.221868e-39 3.610934e-39
[47,] 1.0000000 1.904544e-38 9.522720e-39
[48,] 1.0000000 1.317310e-38 6.586548e-39
[49,] 1.0000000 2.298676e-38 1.149338e-38
[50,] 1.0000000 4.376797e-38 2.188398e-38
[51,] 1.0000000 3.518336e-38 1.759168e-38
[52,] 1.0000000 7.244352e-38 3.622176e-38
[53,] 1.0000000 1.295794e-37 6.478972e-38
[54,] 1.0000000 3.163701e-37 1.581850e-37
[55,] 1.0000000 5.297626e-37 2.648813e-37
[56,] 1.0000000 1.029067e-36 5.145337e-37
[57,] 1.0000000 1.937379e-36 9.686894e-37
[58,] 1.0000000 4.744750e-36 2.372375e-36
[59,] 1.0000000 1.027904e-35 5.139521e-36
[60,] 1.0000000 1.250701e-141 6.253507e-142
[61,] 1.0000000 8.672096e-142 4.336048e-142
[62,] 1.0000000 4.897131e-141 2.448565e-141
[63,] 1.0000000 3.921327e-141 1.960663e-141
[64,] 1.0000000 2.268191e-140 1.134096e-140
[65,] 1.0000000 5.182763e-140 2.591381e-140
[66,] 1.0000000 1.433662e-139 7.168308e-140
[67,] 1.0000000 5.146311e-139 2.573156e-139
[68,] 1.0000000 1.359581e-138 6.797903e-139
[69,] 1.0000000 4.796382e-138 2.398191e-138
[70,] 1.0000000 2.076795e-137 1.038398e-137
[71,] 1.0000000 8.886387e-137 4.443194e-137
[72,] 1.0000000 2.685011e-136 1.342505e-136
[73,] 1.0000000 7.396435e-139 3.698217e-139
[74,] 1.0000000 6.404281e-146 3.202140e-146
[75,] 1.0000000 2.437103e-145 1.218552e-145
[76,] 1.0000000 1.011191e-144 5.055957e-145
[77,] 1.0000000 2.940968e-144 1.470484e-144
[78,] 1.0000000 1.099805e-143 5.499026e-144
[79,] 1.0000000 4.042514e-143 2.021257e-143
[80,] 1.0000000 1.771039e-142 8.855195e-143
[81,] 1.0000000 1.534188e-142 7.670942e-143
[82,] 1.0000000 5.638469e-142 2.819235e-142
[83,] 1.0000000 1.266967e-141 6.334834e-142
[84,] 1.0000000 4.428055e-141 2.214028e-141
[85,] 1.0000000 1.487023e-140 7.435116e-141
[86,] 1.0000000 8.202697e-140 4.101349e-140
[87,] 1.0000000 2.395847e-141 1.197924e-141
[88,] 1.0000000 1.399691e-140 6.998453e-141
[89,] 1.0000000 3.603428e-140 1.801714e-140
[90,] 1.0000000 1.885401e-139 9.427007e-140
[91,] 1.0000000 6.294604e-139 3.147302e-139
[92,] 1.0000000 4.329681e-154 2.164840e-154
[93,] 1.0000000 1.270951e-153 6.354756e-154
[94,] 1.0000000 5.755779e-153 2.877890e-153
[95,] 1.0000000 3.057046e-152 1.528523e-152
[96,] 1.0000000 1.574049e-151 7.870245e-152
[97,] 1.0000000 1.453396e-151 7.266980e-152
[98,] 1.0000000 7.222224e-151 3.611112e-151
[99,] 1.0000000 3.565327e-150 1.782663e-150
[100,] 1.0000000 1.804830e-149 9.024149e-150
[101,] 1.0000000 4.094214e-149 2.047107e-149
[102,] 1.0000000 2.000488e-148 1.000244e-148
[103,] 1.0000000 1.087923e-147 5.439616e-148
[104,] 1.0000000 5.721563e-147 2.860781e-147
[105,] 1.0000000 2.569612e-146 1.284806e-146
[106,] 1.0000000 1.130927e-145 5.654637e-146
[107,] 1.0000000 3.539989e-145 1.769994e-145
[108,] 1.0000000 1.465998e-144 7.329991e-145
[109,] 1.0000000 7.134642e-144 3.567321e-144
[110,] 1.0000000 2.415429e-143 1.207715e-143
[111,] 1.0000000 1.176578e-142 5.882888e-143
[112,] 1.0000000 5.372144e-142 2.686072e-142
[113,] 1.0000000 2.475437e-141 1.237718e-141
[114,] 1.0000000 3.331312e-142 1.665656e-142
[115,] 1.0000000 1.860192e-141 9.300962e-142
[116,] 1.0000000 9.543678e-141 4.771839e-141
[117,] 1.0000000 8.775016e-144 4.387508e-144
[118,] 1.0000000 5.257925e-143 2.628962e-143
[119,] 1.0000000 2.929025e-142 1.464513e-142
[120,] 1.0000000 1.421743e-141 7.108713e-142
[121,] 1.0000000 6.353445e-141 3.176722e-141
[122,] 1.0000000 3.301362e-140 1.650681e-140
[123,] 1.0000000 1.305971e-139 6.529856e-140
[124,] 1.0000000 5.957029e-139 2.978515e-139
[125,] 1.0000000 3.185499e-138 1.592750e-138
[126,] 1.0000000 1.286678e-138 6.433388e-139
[127,] 1.0000000 5.988342e-138 2.994171e-138
[128,] 1.0000000 3.007498e-137 1.503749e-137
[129,] 1.0000000 1.413154e-136 7.065768e-137
[130,] 1.0000000 7.656010e-136 3.828005e-136
[131,] 1.0000000 4.200473e-135 2.100236e-135
[132,] 1.0000000 1.992773e-134 9.963867e-135
[133,] 1.0000000 9.482542e-134 4.741271e-134
[134,] 1.0000000 4.524082e-133 2.262041e-133
[135,] 1.0000000 2.506426e-132 1.253213e-132
[136,] 1.0000000 1.199994e-131 5.999968e-132
[137,] 1.0000000 6.661202e-131 3.330601e-131
[138,] 1.0000000 3.420384e-130 1.710192e-130
[139,] 1.0000000 1.889940e-129 9.449702e-130
[140,] 1.0000000 1.048276e-128 5.241380e-129
[141,] 1.0000000 5.038018e-128 2.519009e-128
[142,] 1.0000000 2.721374e-127 1.360687e-127
[143,] 1.0000000 1.308294e-126 6.541472e-127
[144,] 1.0000000 6.435612e-126 3.217806e-126
[145,] 1.0000000 3.504333e-125 1.752167e-125
[146,] 1.0000000 5.071826e-126 2.535913e-126
[147,] 1.0000000 2.463283e-125 1.231642e-125
[148,] 1.0000000 1.363164e-124 6.815819e-125
[149,] 1.0000000 7.059965e-124 3.529982e-124
[150,] 1.0000000 3.420401e-123 1.710201e-123
[151,] 1.0000000 1.655464e-122 8.277319e-123
[152,] 1.0000000 8.256864e-122 4.128432e-122
[153,] 1.0000000 4.360474e-121 2.180237e-121
[154,] 1.0000000 2.101920e-120 1.050960e-120
[155,] 1.0000000 1.011656e-119 5.058281e-120
[156,] 1.0000000 4.861056e-119 2.430528e-119
[157,] 1.0000000 2.331623e-118 1.165811e-118
[158,] 1.0000000 9.760366e-118 4.880183e-118
[159,] 1.0000000 7.405988e-118 3.702994e-118
[160,] 1.0000000 3.566359e-117 1.783179e-117
[161,] 1.0000000 1.713885e-116 8.569426e-117
[162,] 1.0000000 8.900876e-116 4.450438e-116
[163,] 1.0000000 4.261205e-115 2.130602e-115
[164,] 1.0000000 2.035344e-114 1.017672e-114
[165,] 1.0000000 9.698649e-114 4.849325e-114
[166,] 1.0000000 4.971168e-113 2.485584e-113
[167,] 1.0000000 2.356926e-112 1.178463e-112
[168,] 1.0000000 1.251475e-111 6.257377e-112
[169,] 1.0000000 5.902652e-111 2.951326e-111
[170,] 1.0000000 3.114755e-110 1.557377e-110
[171,] 1.0000000 1.461110e-109 7.305549e-110
[172,] 1.0000000 6.662650e-109 3.331325e-109
[173,] 1.0000000 3.107434e-108 1.553717e-108
[174,] 1.0000000 1.445048e-107 7.225241e-108
[175,] 1.0000000 5.262908e-107 2.631454e-107
[176,] 1.0000000 2.602057e-106 1.301029e-106
[177,] 1.0000000 6.021674e-106 3.010837e-106
[178,] 1.0000000 3.030230e-105 1.515115e-105
[179,] 1.0000000 1.486250e-104 7.431250e-105
[180,] 1.0000000 7.652878e-104 3.826439e-104
[181,] 1.0000000 3.927600e-103 1.963800e-103
[182,] 1.0000000 1.926760e-102 9.633802e-103
[183,] 1.0000000 9.264119e-102 4.632059e-102
[184,] 1.0000000 4.200739e-101 2.100370e-101
[185,] 1.0000000 1.993971e-100 9.969854e-101
[186,] 1.0000000 8.983541e-100 4.491771e-100
[187,] 1.0000000 4.398009e-99 2.199004e-99
[188,] 1.0000000 1.968545e-98 9.842725e-99
[189,] 1.0000000 9.124439e-98 4.562220e-98
[190,] 1.0000000 3.919208e-97 1.959604e-97
[191,] 1.0000000 1.090923e-97 5.454616e-98
[192,] 1.0000000 4.779380e-97 2.389690e-97
[193,] 1.0000000 2.084532e-96 1.042266e-96
[194,] 1.0000000 9.819681e-96 4.909840e-96
[195,] 1.0000000 4.265338e-95 2.132669e-95
[196,] 1.0000000 1.842975e-94 9.214873e-95
[197,] 1.0000000 7.918455e-94 3.959227e-94
[198,] 1.0000000 2.928884e-93 1.464442e-93
[199,] 1.0000000 1.403227e-92 7.016135e-93
[200,] 1.0000000 6.285189e-92 3.142595e-92
[201,] 1.0000000 2.874887e-91 1.437444e-91
[202,] 1.0000000 1.219313e-90 6.096563e-91
[203,] 1.0000000 4.650031e-90 2.325016e-90
[204,] 1.0000000 1.967478e-89 9.837391e-90
[205,] 1.0000000 8.400587e-89 4.200293e-89
[206,] 1.0000000 3.574226e-88 1.787113e-88
[207,] 1.0000000 1.017031e-87 5.085157e-88
[208,] 1.0000000 4.515457e-87 2.257728e-87
[209,] 1.0000000 2.030873e-86 1.015437e-86
[210,] 1.0000000 8.566640e-86 4.283320e-86
[211,] 1.0000000 3.994428e-85 1.997214e-85
[212,] 1.0000000 1.388367e-84 6.941837e-85
[213,] 1.0000000 5.811781e-84 2.905891e-84
[214,] 1.0000000 2.655780e-83 1.327890e-83
[215,] 1.0000000 1.150192e-82 5.750958e-83
[216,] 1.0000000 4.764753e-82 2.382377e-82
[217,] 1.0000000 2.044037e-81 1.022019e-81
[218,] 1.0000000 8.403165e-81 4.201583e-81
[219,] 1.0000000 3.569533e-80 1.784766e-80
[220,] 1.0000000 1.274200e-79 6.371001e-80
[221,] 1.0000000 5.192956e-79 2.596478e-79
[222,] 1.0000000 2.108031e-78 1.054016e-78
[223,] 1.0000000 8.905835e-78 4.452917e-78
[224,] 1.0000000 3.586870e-77 1.793435e-77
[225,] 1.0000000 1.578234e-76 7.891171e-77
[226,] 1.0000000 6.307112e-76 3.153556e-76
[227,] 1.0000000 2.510220e-75 1.255110e-75
[228,] 1.0000000 1.077868e-74 5.389340e-75
[229,] 1.0000000 4.259715e-74 2.129858e-74
[230,] 1.0000000 1.710020e-73 8.550098e-74
[231,] 1.0000000 6.701773e-73 3.350887e-73
[232,] 1.0000000 2.850260e-72 1.425130e-72
[233,] 1.0000000 1.107886e-71 5.539429e-72
[234,] 1.0000000 4.287914e-71 2.143957e-71
[235,] 1.0000000 1.652436e-70 8.262181e-71
[236,] 1.0000000 6.622965e-70 3.311482e-70
[237,] 1.0000000 2.530378e-69 1.265189e-69
[238,] 1.0000000 1.066697e-68 5.333483e-69
[239,] 1.0000000 2.674216e-68 1.337108e-68
[240,] 1.0000000 1.000851e-67 5.004253e-68
[241,] 1.0000000 3.718517e-67 1.859258e-67
[242,] 1.0000000 1.554053e-66 7.770263e-67
[243,] 1.0000000 5.707571e-66 2.853785e-66
[244,] 1.0000000 2.079401e-65 1.039700e-65
[245,] 1.0000000 8.569470e-65 4.284735e-65
[246,] 1.0000000 3.083868e-64 1.541934e-64
[247,] 1.0000000 1.150111e-63 5.750555e-64
[248,] 1.0000000 1.555306e-64 7.776529e-65
[249,] 1.0000000 4.021326e-66 2.010663e-66
[250,] 1.0000000 1.551332e-65 7.756660e-66
[251,] 1.0000000 6.478633e-65 3.239317e-65
[252,] 1.0000000 1.562509e-64 7.812544e-65
[253,] 1.0000000 5.690763e-64 2.845381e-64
[254,] 1.0000000 2.173438e-63 1.086719e-63
[255,] 1.0000000 7.811068e-63 3.905534e-63
[256,] 1.0000000 2.777771e-62 1.388886e-62
[257,] 1.0000000 1.144608e-61 5.723042e-62
[258,] 1.0000000 4.675386e-61 2.337693e-61
[259,] 1.0000000 1.632270e-60 8.161349e-61
[260,] 1.0000000 6.099872e-60 3.049936e-60
[261,] 1.0000000 2.097547e-59 1.048774e-59
[262,] 1.0000000 8.453727e-59 4.226864e-59
[263,] 1.0000000 2.860325e-58 1.430162e-58
[264,] 1.0000000 1.127097e-57 5.635487e-58
[265,] 1.0000000 3.749235e-57 1.874617e-57
[266,] 1.0000000 1.229256e-56 6.146278e-57
[267,] 1.0000000 3.967805e-56 1.983902e-56
[268,] 1.0000000 1.572626e-55 7.863130e-56
[269,] 1.0000000 5.686270e-55 2.843135e-55
[270,] 1.0000000 1.588219e-54 7.941095e-55
[271,] 1.0000000 3.954957e-54 1.977479e-54
[272,] 1.0000000 1.425712e-53 7.128560e-54
[273,] 1.0000000 1.061983e-54 5.309915e-55
[274,] 1.0000000 3.048283e-54 1.524142e-54
[275,] 1.0000000 1.027284e-53 5.136420e-54
[276,] 1.0000000 3.748464e-53 1.874232e-53
[277,] 1.0000000 2.886949e-54 1.443474e-54
[278,] 1.0000000 1.079786e-53 5.398930e-54
[279,] 1.0000000 4.313983e-53 2.156992e-53
[280,] 1.0000000 1.599275e-52 7.996376e-53
[281,] 1.0000000 5.315801e-52 2.657901e-52
[282,] 1.0000000 2.063709e-51 1.031854e-51
[283,] 1.0000000 6.719849e-51 3.359924e-51
[284,] 1.0000000 1.569063e-51 7.845314e-52
[285,] 1.0000000 5.060422e-51 2.530211e-51
[286,] 1.0000000 1.607378e-50 8.036888e-51
[287,] 1.0000000 5.705470e-50 2.852735e-50
[288,] 1.0000000 2.113430e-49 1.056715e-49
[289,] 1.0000000 8.097832e-49 4.048916e-49
[290,] 1.0000000 3.159880e-48 1.579940e-48
[291,] 1.0000000 1.154857e-47 5.774287e-48
[292,] 1.0000000 3.600099e-47 1.800049e-47
[293,] 1.0000000 1.304462e-46 6.522310e-47
[294,] 1.0000000 4.701524e-46 2.350762e-46
[295,] 1.0000000 1.583687e-45 7.918434e-46
[296,] 1.0000000 5.663694e-45 2.831847e-45
[297,] 1.0000000 2.011319e-44 1.005659e-44
[298,] 1.0000000 5.922736e-44 2.961368e-44
[299,] 1.0000000 3.738945e-44 1.869472e-44
[300,] 1.0000000 1.426521e-43 7.132603e-44
[301,] 1.0000000 4.086166e-43 2.043083e-43
[302,] 1.0000000 1.503150e-42 7.515749e-43
[303,] 1.0000000 4.984722e-42 2.492361e-42
[304,] 1.0000000 1.356269e-41 6.781346e-42
[305,] 1.0000000 3.535271e-41 1.767636e-41
[306,] 1.0000000 4.560062e-43 2.280031e-43
[307,] 1.0000000 1.144672e-42 5.723358e-43
[308,] 1.0000000 3.686286e-42 1.843143e-42
[309,] 1.0000000 1.416751e-41 7.083757e-42
[310,] 1.0000000 5.244590e-41 2.622295e-41
[311,] 1.0000000 1.928356e-40 9.641780e-41
[312,] 1.0000000 6.972229e-40 3.486114e-40
[313,] 1.0000000 2.648462e-39 1.324231e-39
[314,] 1.0000000 9.487427e-39 4.743713e-39
[315,] 1.0000000 3.378465e-38 1.689233e-38
[316,] 1.0000000 1.195874e-37 5.979369e-38
[317,] 1.0000000 4.207502e-37 2.103751e-37
[318,] 1.0000000 9.293467e-37 4.646733e-37
[319,] 1.0000000 3.245657e-36 1.622828e-36
[320,] 1.0000000 1.154739e-35 5.773695e-36
[321,] 1.0000000 2.322325e-35 1.161163e-35
[322,] 1.0000000 4.147773e-35 2.073886e-35
[323,] 1.0000000 1.430524e-34 7.152620e-35
[324,] 1.0000000 4.902364e-34 2.451182e-34
[325,] 1.0000000 1.669249e-33 8.346243e-34
[326,] 1.0000000 5.646973e-33 2.823487e-33
[327,] 1.0000000 1.897856e-32 9.489280e-33
[328,] 1.0000000 6.660870e-32 3.330435e-32
[329,] 1.0000000 2.270410e-31 1.135205e-31
[330,] 1.0000000 7.481771e-31 3.740885e-31
[331,] 1.0000000 2.604751e-30 1.302375e-30
[332,] 1.0000000 8.471710e-30 4.235855e-30
[333,] 1.0000000 2.770070e-29 1.385035e-29
[334,] 1.0000000 9.399415e-29 4.699708e-29
[335,] 1.0000000 3.140567e-28 1.570284e-28
[336,] 1.0000000 9.968470e-28 4.984235e-28
[337,] 1.0000000 3.141401e-27 1.570700e-27
[338,] 1.0000000 9.827847e-27 4.913924e-27
[339,] 1.0000000 3.052100e-26 1.526050e-26
[340,] 1.0000000 9.408236e-26 4.704118e-26
[341,] 1.0000000 3.042023e-25 1.521012e-25
[342,] 1.0000000 9.240190e-25 4.620095e-25
[343,] 1.0000000 2.785207e-24 1.392604e-24
[344,] 1.0000000 8.330128e-24 4.165064e-24
[345,] 1.0000000 2.646414e-23 1.323207e-23
[346,] 1.0000000 8.370806e-23 4.185403e-23
[347,] 1.0000000 2.623470e-22 1.311735e-22
[348,] 1.0000000 7.754081e-22 3.877040e-22
[349,] 1.0000000 2.236804e-21 1.118402e-21
[350,] 1.0000000 6.398850e-21 3.199425e-21
[351,] 1.0000000 6.028194e-21 3.014097e-21
[352,] 1.0000000 1.741338e-20 8.706692e-21
[353,] 1.0000000 4.986935e-20 2.493467e-20
[354,] 1.0000000 1.415745e-19 7.078723e-20
[355,] 1.0000000 4.163113e-19 2.081557e-19
[356,] 1.0000000 1.161342e-18 5.806711e-19
[357,] 1.0000000 3.210207e-18 1.605104e-18
[358,] 1.0000000 5.750857e-18 2.875429e-18
[359,] 1.0000000 1.691906e-17 8.459529e-18
[360,] 1.0000000 4.589172e-17 2.294586e-17
[361,] 1.0000000 1.232807e-16 6.164036e-17
[362,] 1.0000000 3.279338e-16 1.639669e-16
[363,] 1.0000000 8.636370e-16 4.318185e-16
[364,] 1.0000000 2.251390e-15 1.125695e-15
[365,] 1.0000000 5.808450e-15 2.904225e-15
[366,] 1.0000000 1.101481e-14 5.507404e-15
[367,] 1.0000000 3.037083e-14 1.518542e-14
[368,] 1.0000000 7.756398e-14 3.878199e-14
[369,] 1.0000000 1.952448e-13 9.762238e-14
[370,] 1.0000000 4.859602e-13 2.429801e-13
[371,] 1.0000000 1.220539e-12 6.102694e-13
[372,] 1.0000000 3.191122e-12 1.595561e-12
[373,] 1.0000000 7.735629e-12 3.867815e-12
[374,] 1.0000000 1.852535e-11 9.262676e-12
[375,] 1.0000000 4.381506e-11 2.190753e-11
[376,] 1.0000000 3.955096e-11 1.977548e-11
[377,] 1.0000000 1.147273e-11 5.736365e-12
[378,] 1.0000000 2.934782e-11 1.467391e-11
[379,] 1.0000000 7.713928e-11 3.856964e-11
[380,] 1.0000000 1.801751e-10 9.008757e-11
[381,] 1.0000000 4.667238e-10 2.333619e-10
[382,] 1.0000000 1.170567e-09 5.852837e-10
[383,] 1.0000000 2.905017e-09 1.452509e-09
[384,] 1.0000000 7.299480e-09 3.649740e-09
[385,] 1.0000000 1.807066e-08 9.035328e-09
[386,] 1.0000000 4.370946e-08 2.185473e-08
[387,] 0.9999999 1.030032e-07 5.150162e-08
[388,] 0.9999999 2.403812e-07 1.201906e-07
[389,] 0.9999999 2.666555e-07 1.333277e-07
[390,] 0.9999997 6.320726e-07 3.160363e-07
[391,] 0.9999993 1.474829e-06 7.374144e-07
[392,] 1.0000000 9.363543e-08 4.681771e-08
[393,] 0.9999999 2.421164e-07 1.210582e-07
[394,] 1.0000000 1.320441e-09 6.602203e-10
[395,] 1.0000000 3.890420e-09 1.945210e-09
[396,] 1.0000000 1.205808e-08 6.029042e-09
[397,] 1.0000000 3.624668e-08 1.812334e-08
[398,] 0.9999999 1.094402e-07 5.472012e-08
[399,] 0.9999998 3.160931e-07 1.580465e-07
[400,] 0.9999995 9.210285e-07 4.605143e-07
[401,] 0.9999987 2.672100e-06 1.336050e-06
[402,] 0.9999964 7.242726e-06 3.621363e-06
[403,] 0.9999902 1.969726e-05 9.848629e-06
[404,] 0.9999762 4.753620e-05 2.376810e-05
[405,] 0.9999764 4.715663e-05 2.357831e-05
[406,] 0.9999622 7.554642e-05 3.777321e-05
[407,] 0.9998989 2.021139e-04 1.010570e-04
[408,] 0.9997311 5.377719e-04 2.688859e-04
[409,] 0.9993360 1.328068e-03 6.640339e-04
[410,] 0.9984744 3.051271e-03 1.525636e-03
[411,] 0.9975549 4.890260e-03 2.445130e-03
[412,] 0.9946803 1.063941e-02 5.319704e-03
[413,] 0.9876586 2.468280e-02 1.234140e-02
[414,] 0.9734041 5.319178e-02 2.659589e-02
[415,] 0.9482038 1.035924e-01 5.179621e-02
[416,] 0.9274583 1.450835e-01 7.254173e-02
[417,] 0.9959520 8.096072e-03 4.048036e-03
[418,] 0.9916997 1.660060e-02 8.300302e-03
> postscript(file="/var/www/html/freestat/rcomp/tmp/1b3iv1291317431.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/2b3iv1291317431.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/3lczy1291317431.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/4lczy1291317431.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/5lczy1291317431.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 6
-103.2385301 78.6110629 -193.0483656 29.3169664 -118.2528983 -191.6364156
7 8 9 10 11 12
157.8412239 -53.3358500 176.3545702 -16.1011360 -151.7357115 11.1840122
13 14 15 16 17 18
121.8020360 -55.3738445 296.9647588 4.8001512 -47.6667778 319.6081477
19 20 21 22 23 24
-17.9383429 209.4863290 2.4800653 -94.9479725 -47.9857649 145.1398071
25 26 27 28 29 30
9.5256328 -62.4682582 64.4450331 46.1139352 -27.5903855 249.6574688
31 32 33 34 35 36
-25.1191149 -33.0432926 -28.5724041 -17.8856722 -106.7450658 -5.3742185
37 38 39 40 41 42
-7.9860084 4.4478303 -8.9388896 92.0017448 718.9778225 -64.1458767
43 44 45 46 47 48
164.3304000 17.9768446 9.4424685 -36.8708145 -61.9475680 221.0689149
49 50 51 52 53 54
8.5897065 -6.5232862 47.9687929 72.6333328 -5.1928587 100.5541410
55 56 57 58 59 60
68.4877525 62.7895433 -30.0957386 1.0097033 65.0723661 7.3193113
61 62 63 64 65 66
-23.8946707 10.0091245 -20.4304021 5.9561988 43.7411938 827.2481558
67 68 69 70 71 72
42.2644065 3.5391411 38.5408987 0.1345086 30.8938355 5.0960188
73 74 75 76 77 78
-11.8626854 -17.9224420 -19.3312313 -18.1779932 0.8451727 14.4366059
79 80 81 82 83 84
87.0558872 124.0143014 -1.2469427 10.2851672 4.1358990 3.0966810
85 86 87 88 89 90
-22.9799469 12.2290963 31.7390770 -24.6480353 13.7101904 -27.5866134
91 92 93 94 95 96
6.7775631 -7.5723418 64.7856116 -15.4249588 23.6317018 2.7249456
97 98 99 100 101 102
-32.4240669 164.9135253 -40.8120314 -23.0846911 -1.4076433 -13.2584274
103 104 105 106 107 108
20.3629300 -24.7933030 -1.8146082 -4.1520971 18.9601902 -28.3983981
109 110 111 112 113 114
-12.1869943 -1.7065100 2.9992310 8.9964610 8.2798501 -27.2884462
115 116 117 118 119 120
-18.8804091 12.7557697 -4.1608512 -23.0681378 -25.4472065 45.3289581
121 122 123 124 125 126
-17.3466078 -12.2128149 62.8686143 -20.6326255 -1.5683111 -25.1008658
127 128 129 130 131 132
-17.7715239 -11.4389022 8.4712393 -17.7715239 -9.9665871 36.6716285
133 134 135 136 137 138
-17.7715239 -15.1299894 -17.7715239 -10.8318097 -5.6420026 -17.7715239
139 140 141 142 143 144
-17.7715239 -17.7715239 -5.7245588 -17.7715239 -7.9957263 -15.0852406
145 146 147 148 149 150
-6.9403875 -7.2042128 -17.7715239 -11.8784102 -17.7715239 -16.9844017
151 152 153 154 155 156
-10.4788181 43.3483659 -17.7715239 -9.0478457 -14.9941017 -17.7715239
157 158 159 160 161 162
-17.7715239 -16.7434713 -13.5715006 -17.7715239 -17.7715239 -17.7715239
163 164 165 166 167 168
-17.7715239 5.9047961 29.1341866 -17.7715239 -17.7715239 2.0702748
169 170 171 172 173 174
-17.7715239 -17.7715239 -17.7715239 -14.3519010 -17.7715239 -10.6347972
175 176 177 178 179 180
-17.7715239 -10.4746235 -17.7715239 -19.2098605 -17.7715239 -17.7715239
181 182 183 184 185 186
5.2809704 -17.6701883 18.0882501 -16.0784351 -15.6042886 -9.9109276
187 188 189 190 191 192
-10.1173556 -14.3090944 -16.3061709 -17.7715239 -15.1534012 -17.7715239
193 194 195 196 197 198
-13.0911907 -17.7715239 -16.0293232 -19.7177719 21.8425098 -19.7177719
199 200 201 202 203 204
-19.7177719 -10.6143344 -19.7177719 -19.7177719 -19.7177719 9.8502926
205 206 207 208 209 210
-16.7006246 -18.1582073 -18.1228516 -19.7177719 9.4228224 -19.7177719
211 212 213 214 215 216
-17.7715239 -17.7715239 13.6412053 -1.4905944 -15.1539497 -17.7715239
217 218 219 220 221 222
-9.9175622 7.5620269 -17.7715239 -3.1874263 -15.6436606 -17.7715239
223 224 225 226 227 228
-15.7831400 -17.7715239 -15.9411677 6.0684338 -17.7715239 -17.7715239
229 230 231 232 233 234
-16.8886255 -17.7715239 -12.1957821 -17.7715239 -17.7715239 -3.6549638
235 236 237 238 239 240
-17.7715239 -16.7942593 -17.7715239 -11.9837720 -17.7715239 -17.7715239
241 242 243 244 245 246
-17.7715239 -15.6645685 -17.7715239 -9.5755039 -34.3674831 -19.7177719
247 248 249 250 251 252
-19.7177719 -8.0941280 -19.7177719 -19.7177719 -7.1416207 -19.7177719
253 254 255 256 257 258
-18.3085199 46.1144388 52.4579341 -17.7715239 -9.9019202 19.9366442
259 260 261 262 263 264
-19.7177719 -17.7715239 -19.7177719 -19.7177719 -12.1279946 -4.1143907
265 266 267 268 269 270
-19.7177719 -17.7715239 -19.7177719 -14.7883112 -19.7177719 -16.3292489
271 272 273 274 275 276
-19.7177719 -19.7177719 -19.7177719 -6.6514034 -17.7715239 9.6241665
277 278 279 280 281 282
11.5785054 -17.7715239 54.0868392 6.0817440 -19.7177719 -17.7715239
283 284 285 286 287 288
54.1272101 -17.7715239 -11.0490447 -17.7715239 -19.7177719 -14.5025951
289 290 291 292 293 294
-19.7177719 45.0836053 -19.7177719 4.2638542 -7.9689040 -17.7715239
295 296 297 298 299 300
-3.6482575 -0.7725222 -17.7715239 -19.7177719 -17.7715239 -17.7715239
301 302 303 304 305 306
-15.9525365 -16.9726810 -17.7715239 -19.7177719 32.3100615 -8.7910753
307 308 309 310 311 312
-19.7177719 -9.4488344 4.0235829 -19.7177719 -19.7177719 52.2140070
313 314 315 316 317 318
-19.7177719 -13.3985432 -13.8082140 -6.7721475 -16.6974894 -3.4484015
319 320 321 322 323 324
-11.7779166 -17.7715239 -17.7715239 -17.7715239 -17.7715239 -19.7177719
325 326 327 328 329 330
-17.7715239 -15.9217150 -19.7177719 -19.7177719 -17.7715239 -17.7715239
331 332 333 334 335 336
-17.7715239 -17.7715239 -17.7715239 -13.2877113 -15.7811006 -17.7715239
337 338 339 340 341 342
-11.4031594 -17.7715239 -17.0181298 -5.2279308 -5.7309948 -17.7715239
343 344 345 346 347 348
-17.7715239 -17.7715239 -17.7715239 -17.7715239 -11.9420184 -17.7715239
349 350 351 352 353 354
-17.7715239 -17.7715239 -5.0456600 -8.3397571 -9.3099043 -17.2235769
355 356 357 358 359 360
-17.7715239 -17.7715239 23.5501756 -17.7715239 -17.7715239 -17.7715239
361 362 363 364 365 366
-13.7791106 -17.7715239 -17.7715239 -17.4618720 -7.5620163 -17.7715239
367 368 369 370 371 372
-17.7715239 -17.7715239 -17.7715239 -17.7715239 -17.7715239 19.5599725
373 374 375 376 377 378
-16.6507693 -19.5849049 -17.7715239 -17.7715239 7.2145218 -10.7804364
379 380 381 382 383 384
-17.7715239 -17.7715239 -17.7715239 33.5696334 66.0097200 -17.7715239
385 386 387 388 389 390
-3.0892383 2.5070302 -12.2396209 -15.2035233 -15.1322970 -4.8040735
391 392 393 394 395 396
-11.2594462 -12.2721206 -15.0868689 -17.4055818 -40.0900733 -7.6868146
397 398 399 400 401 402
-6.0135264 64.4297376 -13.6174969 119.6534186 -23.8350537 -15.4380758
403 404 405 406 407 408
-13.3116905 -8.5154968 -16.2045102 -16.7677685 -8.8879719 -12.7750886
409 410 411 412 413 414
-8.7906419 -23.3870447 30.7412973 -39.7864245 -9.2689253 -14.9699204
415 416 417 418 419 420
-27.4902067 -26.4069997 36.2550853 -11.0480680 -3.7731238 -21.3297682
421 422 423 424 425 426
-26.2523906 25.1544932 1.9547430 -34.1412990 -6.9383195 41.4744236
427 428 429 430 431
-11.1591966 -75.1129514 -40.0773534 -37.8889397 52.8080066
> postscript(file="/var/www/html/freestat/rcomp/tmp/6emzj1291317431.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 431
Frequency = 1
lag(myerror, k = 1) myerror
0 -103.2385301 NA
1 78.6110629 -103.2385301
2 -193.0483656 78.6110629
3 29.3169664 -193.0483656
4 -118.2528983 29.3169664
5 -191.6364156 -118.2528983
6 157.8412239 -191.6364156
7 -53.3358500 157.8412239
8 176.3545702 -53.3358500
9 -16.1011360 176.3545702
10 -151.7357115 -16.1011360
11 11.1840122 -151.7357115
12 121.8020360 11.1840122
13 -55.3738445 121.8020360
14 296.9647588 -55.3738445
15 4.8001512 296.9647588
16 -47.6667778 4.8001512
17 319.6081477 -47.6667778
18 -17.9383429 319.6081477
19 209.4863290 -17.9383429
20 2.4800653 209.4863290
21 -94.9479725 2.4800653
22 -47.9857649 -94.9479725
23 145.1398071 -47.9857649
24 9.5256328 145.1398071
25 -62.4682582 9.5256328
26 64.4450331 -62.4682582
27 46.1139352 64.4450331
28 -27.5903855 46.1139352
29 249.6574688 -27.5903855
30 -25.1191149 249.6574688
31 -33.0432926 -25.1191149
32 -28.5724041 -33.0432926
33 -17.8856722 -28.5724041
34 -106.7450658 -17.8856722
35 -5.3742185 -106.7450658
36 -7.9860084 -5.3742185
37 4.4478303 -7.9860084
38 -8.9388896 4.4478303
39 92.0017448 -8.9388896
40 718.9778225 92.0017448
41 -64.1458767 718.9778225
42 164.3304000 -64.1458767
43 17.9768446 164.3304000
44 9.4424685 17.9768446
45 -36.8708145 9.4424685
46 -61.9475680 -36.8708145
47 221.0689149 -61.9475680
48 8.5897065 221.0689149
49 -6.5232862 8.5897065
50 47.9687929 -6.5232862
51 72.6333328 47.9687929
52 -5.1928587 72.6333328
53 100.5541410 -5.1928587
54 68.4877525 100.5541410
55 62.7895433 68.4877525
56 -30.0957386 62.7895433
57 1.0097033 -30.0957386
58 65.0723661 1.0097033
59 7.3193113 65.0723661
60 -23.8946707 7.3193113
61 10.0091245 -23.8946707
62 -20.4304021 10.0091245
63 5.9561988 -20.4304021
64 43.7411938 5.9561988
65 827.2481558 43.7411938
66 42.2644065 827.2481558
67 3.5391411 42.2644065
68 38.5408987 3.5391411
69 0.1345086 38.5408987
70 30.8938355 0.1345086
71 5.0960188 30.8938355
72 -11.8626854 5.0960188
73 -17.9224420 -11.8626854
74 -19.3312313 -17.9224420
75 -18.1779932 -19.3312313
76 0.8451727 -18.1779932
77 14.4366059 0.8451727
78 87.0558872 14.4366059
79 124.0143014 87.0558872
80 -1.2469427 124.0143014
81 10.2851672 -1.2469427
82 4.1358990 10.2851672
83 3.0966810 4.1358990
84 -22.9799469 3.0966810
85 12.2290963 -22.9799469
86 31.7390770 12.2290963
87 -24.6480353 31.7390770
88 13.7101904 -24.6480353
89 -27.5866134 13.7101904
90 6.7775631 -27.5866134
91 -7.5723418 6.7775631
92 64.7856116 -7.5723418
93 -15.4249588 64.7856116
94 23.6317018 -15.4249588
95 2.7249456 23.6317018
96 -32.4240669 2.7249456
97 164.9135253 -32.4240669
98 -40.8120314 164.9135253
99 -23.0846911 -40.8120314
100 -1.4076433 -23.0846911
101 -13.2584274 -1.4076433
102 20.3629300 -13.2584274
103 -24.7933030 20.3629300
104 -1.8146082 -24.7933030
105 -4.1520971 -1.8146082
106 18.9601902 -4.1520971
107 -28.3983981 18.9601902
108 -12.1869943 -28.3983981
109 -1.7065100 -12.1869943
110 2.9992310 -1.7065100
111 8.9964610 2.9992310
112 8.2798501 8.9964610
113 -27.2884462 8.2798501
114 -18.8804091 -27.2884462
115 12.7557697 -18.8804091
116 -4.1608512 12.7557697
117 -23.0681378 -4.1608512
118 -25.4472065 -23.0681378
119 45.3289581 -25.4472065
120 -17.3466078 45.3289581
121 -12.2128149 -17.3466078
122 62.8686143 -12.2128149
123 -20.6326255 62.8686143
124 -1.5683111 -20.6326255
125 -25.1008658 -1.5683111
126 -17.7715239 -25.1008658
127 -11.4389022 -17.7715239
128 8.4712393 -11.4389022
129 -17.7715239 8.4712393
130 -9.9665871 -17.7715239
131 36.6716285 -9.9665871
132 -17.7715239 36.6716285
133 -15.1299894 -17.7715239
134 -17.7715239 -15.1299894
135 -10.8318097 -17.7715239
136 -5.6420026 -10.8318097
137 -17.7715239 -5.6420026
138 -17.7715239 -17.7715239
139 -17.7715239 -17.7715239
140 -5.7245588 -17.7715239
141 -17.7715239 -5.7245588
142 -7.9957263 -17.7715239
143 -15.0852406 -7.9957263
144 -6.9403875 -15.0852406
145 -7.2042128 -6.9403875
146 -17.7715239 -7.2042128
147 -11.8784102 -17.7715239
148 -17.7715239 -11.8784102
149 -16.9844017 -17.7715239
150 -10.4788181 -16.9844017
151 43.3483659 -10.4788181
152 -17.7715239 43.3483659
153 -9.0478457 -17.7715239
154 -14.9941017 -9.0478457
155 -17.7715239 -14.9941017
156 -17.7715239 -17.7715239
157 -16.7434713 -17.7715239
158 -13.5715006 -16.7434713
159 -17.7715239 -13.5715006
160 -17.7715239 -17.7715239
161 -17.7715239 -17.7715239
162 -17.7715239 -17.7715239
163 5.9047961 -17.7715239
164 29.1341866 5.9047961
165 -17.7715239 29.1341866
166 -17.7715239 -17.7715239
167 2.0702748 -17.7715239
168 -17.7715239 2.0702748
169 -17.7715239 -17.7715239
170 -17.7715239 -17.7715239
171 -14.3519010 -17.7715239
172 -17.7715239 -14.3519010
173 -10.6347972 -17.7715239
174 -17.7715239 -10.6347972
175 -10.4746235 -17.7715239
176 -17.7715239 -10.4746235
177 -19.2098605 -17.7715239
178 -17.7715239 -19.2098605
179 -17.7715239 -17.7715239
180 5.2809704 -17.7715239
181 -17.6701883 5.2809704
182 18.0882501 -17.6701883
183 -16.0784351 18.0882501
184 -15.6042886 -16.0784351
185 -9.9109276 -15.6042886
186 -10.1173556 -9.9109276
187 -14.3090944 -10.1173556
188 -16.3061709 -14.3090944
189 -17.7715239 -16.3061709
190 -15.1534012 -17.7715239
191 -17.7715239 -15.1534012
192 -13.0911907 -17.7715239
193 -17.7715239 -13.0911907
194 -16.0293232 -17.7715239
195 -19.7177719 -16.0293232
196 21.8425098 -19.7177719
197 -19.7177719 21.8425098
198 -19.7177719 -19.7177719
199 -10.6143344 -19.7177719
200 -19.7177719 -10.6143344
201 -19.7177719 -19.7177719
202 -19.7177719 -19.7177719
203 9.8502926 -19.7177719
204 -16.7006246 9.8502926
205 -18.1582073 -16.7006246
206 -18.1228516 -18.1582073
207 -19.7177719 -18.1228516
208 9.4228224 -19.7177719
209 -19.7177719 9.4228224
210 -17.7715239 -19.7177719
211 -17.7715239 -17.7715239
212 13.6412053 -17.7715239
213 -1.4905944 13.6412053
214 -15.1539497 -1.4905944
215 -17.7715239 -15.1539497
216 -9.9175622 -17.7715239
217 7.5620269 -9.9175622
218 -17.7715239 7.5620269
219 -3.1874263 -17.7715239
220 -15.6436606 -3.1874263
221 -17.7715239 -15.6436606
222 -15.7831400 -17.7715239
223 -17.7715239 -15.7831400
224 -15.9411677 -17.7715239
225 6.0684338 -15.9411677
226 -17.7715239 6.0684338
227 -17.7715239 -17.7715239
228 -16.8886255 -17.7715239
229 -17.7715239 -16.8886255
230 -12.1957821 -17.7715239
231 -17.7715239 -12.1957821
232 -17.7715239 -17.7715239
233 -3.6549638 -17.7715239
234 -17.7715239 -3.6549638
235 -16.7942593 -17.7715239
236 -17.7715239 -16.7942593
237 -11.9837720 -17.7715239
238 -17.7715239 -11.9837720
239 -17.7715239 -17.7715239
240 -17.7715239 -17.7715239
241 -15.6645685 -17.7715239
242 -17.7715239 -15.6645685
243 -9.5755039 -17.7715239
244 -34.3674831 -9.5755039
245 -19.7177719 -34.3674831
246 -19.7177719 -19.7177719
247 -8.0941280 -19.7177719
248 -19.7177719 -8.0941280
249 -19.7177719 -19.7177719
250 -7.1416207 -19.7177719
251 -19.7177719 -7.1416207
252 -18.3085199 -19.7177719
253 46.1144388 -18.3085199
254 52.4579341 46.1144388
255 -17.7715239 52.4579341
256 -9.9019202 -17.7715239
257 19.9366442 -9.9019202
258 -19.7177719 19.9366442
259 -17.7715239 -19.7177719
260 -19.7177719 -17.7715239
261 -19.7177719 -19.7177719
262 -12.1279946 -19.7177719
263 -4.1143907 -12.1279946
264 -19.7177719 -4.1143907
265 -17.7715239 -19.7177719
266 -19.7177719 -17.7715239
267 -14.7883112 -19.7177719
268 -19.7177719 -14.7883112
269 -16.3292489 -19.7177719
270 -19.7177719 -16.3292489
271 -19.7177719 -19.7177719
272 -19.7177719 -19.7177719
273 -6.6514034 -19.7177719
274 -17.7715239 -6.6514034
275 9.6241665 -17.7715239
276 11.5785054 9.6241665
277 -17.7715239 11.5785054
278 54.0868392 -17.7715239
279 6.0817440 54.0868392
280 -19.7177719 6.0817440
281 -17.7715239 -19.7177719
282 54.1272101 -17.7715239
283 -17.7715239 54.1272101
284 -11.0490447 -17.7715239
285 -17.7715239 -11.0490447
286 -19.7177719 -17.7715239
287 -14.5025951 -19.7177719
288 -19.7177719 -14.5025951
289 45.0836053 -19.7177719
290 -19.7177719 45.0836053
291 4.2638542 -19.7177719
292 -7.9689040 4.2638542
293 -17.7715239 -7.9689040
294 -3.6482575 -17.7715239
295 -0.7725222 -3.6482575
296 -17.7715239 -0.7725222
297 -19.7177719 -17.7715239
298 -17.7715239 -19.7177719
299 -17.7715239 -17.7715239
300 -15.9525365 -17.7715239
301 -16.9726810 -15.9525365
302 -17.7715239 -16.9726810
303 -19.7177719 -17.7715239
304 32.3100615 -19.7177719
305 -8.7910753 32.3100615
306 -19.7177719 -8.7910753
307 -9.4488344 -19.7177719
308 4.0235829 -9.4488344
309 -19.7177719 4.0235829
310 -19.7177719 -19.7177719
311 52.2140070 -19.7177719
312 -19.7177719 52.2140070
313 -13.3985432 -19.7177719
314 -13.8082140 -13.3985432
315 -6.7721475 -13.8082140
316 -16.6974894 -6.7721475
317 -3.4484015 -16.6974894
318 -11.7779166 -3.4484015
319 -17.7715239 -11.7779166
320 -17.7715239 -17.7715239
321 -17.7715239 -17.7715239
322 -17.7715239 -17.7715239
323 -19.7177719 -17.7715239
324 -17.7715239 -19.7177719
325 -15.9217150 -17.7715239
326 -19.7177719 -15.9217150
327 -19.7177719 -19.7177719
328 -17.7715239 -19.7177719
329 -17.7715239 -17.7715239
330 -17.7715239 -17.7715239
331 -17.7715239 -17.7715239
332 -17.7715239 -17.7715239
333 -13.2877113 -17.7715239
334 -15.7811006 -13.2877113
335 -17.7715239 -15.7811006
336 -11.4031594 -17.7715239
337 -17.7715239 -11.4031594
338 -17.0181298 -17.7715239
339 -5.2279308 -17.0181298
340 -5.7309948 -5.2279308
341 -17.7715239 -5.7309948
342 -17.7715239 -17.7715239
343 -17.7715239 -17.7715239
344 -17.7715239 -17.7715239
345 -17.7715239 -17.7715239
346 -11.9420184 -17.7715239
347 -17.7715239 -11.9420184
348 -17.7715239 -17.7715239
349 -17.7715239 -17.7715239
350 -5.0456600 -17.7715239
351 -8.3397571 -5.0456600
352 -9.3099043 -8.3397571
353 -17.2235769 -9.3099043
354 -17.7715239 -17.2235769
355 -17.7715239 -17.7715239
356 23.5501756 -17.7715239
357 -17.7715239 23.5501756
358 -17.7715239 -17.7715239
359 -17.7715239 -17.7715239
360 -13.7791106 -17.7715239
361 -17.7715239 -13.7791106
362 -17.7715239 -17.7715239
363 -17.4618720 -17.7715239
364 -7.5620163 -17.4618720
365 -17.7715239 -7.5620163
366 -17.7715239 -17.7715239
367 -17.7715239 -17.7715239
368 -17.7715239 -17.7715239
369 -17.7715239 -17.7715239
370 -17.7715239 -17.7715239
371 19.5599725 -17.7715239
372 -16.6507693 19.5599725
373 -19.5849049 -16.6507693
374 -17.7715239 -19.5849049
375 -17.7715239 -17.7715239
376 7.2145218 -17.7715239
377 -10.7804364 7.2145218
378 -17.7715239 -10.7804364
379 -17.7715239 -17.7715239
380 -17.7715239 -17.7715239
381 33.5696334 -17.7715239
382 66.0097200 33.5696334
383 -17.7715239 66.0097200
384 -3.0892383 -17.7715239
385 2.5070302 -3.0892383
386 -12.2396209 2.5070302
387 -15.2035233 -12.2396209
388 -15.1322970 -15.2035233
389 -4.8040735 -15.1322970
390 -11.2594462 -4.8040735
391 -12.2721206 -11.2594462
392 -15.0868689 -12.2721206
393 -17.4055818 -15.0868689
394 -40.0900733 -17.4055818
395 -7.6868146 -40.0900733
396 -6.0135264 -7.6868146
397 64.4297376 -6.0135264
398 -13.6174969 64.4297376
399 119.6534186 -13.6174969
400 -23.8350537 119.6534186
401 -15.4380758 -23.8350537
402 -13.3116905 -15.4380758
403 -8.5154968 -13.3116905
404 -16.2045102 -8.5154968
405 -16.7677685 -16.2045102
406 -8.8879719 -16.7677685
407 -12.7750886 -8.8879719
408 -8.7906419 -12.7750886
409 -23.3870447 -8.7906419
410 30.7412973 -23.3870447
411 -39.7864245 30.7412973
412 -9.2689253 -39.7864245
413 -14.9699204 -9.2689253
414 -27.4902067 -14.9699204
415 -26.4069997 -27.4902067
416 36.2550853 -26.4069997
417 -11.0480680 36.2550853
418 -3.7731238 -11.0480680
419 -21.3297682 -3.7731238
420 -26.2523906 -21.3297682
421 25.1544932 -26.2523906
422 1.9547430 25.1544932
423 -34.1412990 1.9547430
424 -6.9383195 -34.1412990
425 41.4744236 -6.9383195
426 -11.1591966 41.4744236
427 -75.1129514 -11.1591966
428 -40.0773534 -75.1129514
429 -37.8889397 -40.0773534
430 52.8080066 -37.8889397
431 NA 52.8080066
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 78.6110629 -103.2385301
[2,] -193.0483656 78.6110629
[3,] 29.3169664 -193.0483656
[4,] -118.2528983 29.3169664
[5,] -191.6364156 -118.2528983
[6,] 157.8412239 -191.6364156
[7,] -53.3358500 157.8412239
[8,] 176.3545702 -53.3358500
[9,] -16.1011360 176.3545702
[10,] -151.7357115 -16.1011360
[11,] 11.1840122 -151.7357115
[12,] 121.8020360 11.1840122
[13,] -55.3738445 121.8020360
[14,] 296.9647588 -55.3738445
[15,] 4.8001512 296.9647588
[16,] -47.6667778 4.8001512
[17,] 319.6081477 -47.6667778
[18,] -17.9383429 319.6081477
[19,] 209.4863290 -17.9383429
[20,] 2.4800653 209.4863290
[21,] -94.9479725 2.4800653
[22,] -47.9857649 -94.9479725
[23,] 145.1398071 -47.9857649
[24,] 9.5256328 145.1398071
[25,] -62.4682582 9.5256328
[26,] 64.4450331 -62.4682582
[27,] 46.1139352 64.4450331
[28,] -27.5903855 46.1139352
[29,] 249.6574688 -27.5903855
[30,] -25.1191149 249.6574688
[31,] -33.0432926 -25.1191149
[32,] -28.5724041 -33.0432926
[33,] -17.8856722 -28.5724041
[34,] -106.7450658 -17.8856722
[35,] -5.3742185 -106.7450658
[36,] -7.9860084 -5.3742185
[37,] 4.4478303 -7.9860084
[38,] -8.9388896 4.4478303
[39,] 92.0017448 -8.9388896
[40,] 718.9778225 92.0017448
[41,] -64.1458767 718.9778225
[42,] 164.3304000 -64.1458767
[43,] 17.9768446 164.3304000
[44,] 9.4424685 17.9768446
[45,] -36.8708145 9.4424685
[46,] -61.9475680 -36.8708145
[47,] 221.0689149 -61.9475680
[48,] 8.5897065 221.0689149
[49,] -6.5232862 8.5897065
[50,] 47.9687929 -6.5232862
[51,] 72.6333328 47.9687929
[52,] -5.1928587 72.6333328
[53,] 100.5541410 -5.1928587
[54,] 68.4877525 100.5541410
[55,] 62.7895433 68.4877525
[56,] -30.0957386 62.7895433
[57,] 1.0097033 -30.0957386
[58,] 65.0723661 1.0097033
[59,] 7.3193113 65.0723661
[60,] -23.8946707 7.3193113
[61,] 10.0091245 -23.8946707
[62,] -20.4304021 10.0091245
[63,] 5.9561988 -20.4304021
[64,] 43.7411938 5.9561988
[65,] 827.2481558 43.7411938
[66,] 42.2644065 827.2481558
[67,] 3.5391411 42.2644065
[68,] 38.5408987 3.5391411
[69,] 0.1345086 38.5408987
[70,] 30.8938355 0.1345086
[71,] 5.0960188 30.8938355
[72,] -11.8626854 5.0960188
[73,] -17.9224420 -11.8626854
[74,] -19.3312313 -17.9224420
[75,] -18.1779932 -19.3312313
[76,] 0.8451727 -18.1779932
[77,] 14.4366059 0.8451727
[78,] 87.0558872 14.4366059
[79,] 124.0143014 87.0558872
[80,] -1.2469427 124.0143014
[81,] 10.2851672 -1.2469427
[82,] 4.1358990 10.2851672
[83,] 3.0966810 4.1358990
[84,] -22.9799469 3.0966810
[85,] 12.2290963 -22.9799469
[86,] 31.7390770 12.2290963
[87,] -24.6480353 31.7390770
[88,] 13.7101904 -24.6480353
[89,] -27.5866134 13.7101904
[90,] 6.7775631 -27.5866134
[91,] -7.5723418 6.7775631
[92,] 64.7856116 -7.5723418
[93,] -15.4249588 64.7856116
[94,] 23.6317018 -15.4249588
[95,] 2.7249456 23.6317018
[96,] -32.4240669 2.7249456
[97,] 164.9135253 -32.4240669
[98,] -40.8120314 164.9135253
[99,] -23.0846911 -40.8120314
[100,] -1.4076433 -23.0846911
[101,] -13.2584274 -1.4076433
[102,] 20.3629300 -13.2584274
[103,] -24.7933030 20.3629300
[104,] -1.8146082 -24.7933030
[105,] -4.1520971 -1.8146082
[106,] 18.9601902 -4.1520971
[107,] -28.3983981 18.9601902
[108,] -12.1869943 -28.3983981
[109,] -1.7065100 -12.1869943
[110,] 2.9992310 -1.7065100
[111,] 8.9964610 2.9992310
[112,] 8.2798501 8.9964610
[113,] -27.2884462 8.2798501
[114,] -18.8804091 -27.2884462
[115,] 12.7557697 -18.8804091
[116,] -4.1608512 12.7557697
[117,] -23.0681378 -4.1608512
[118,] -25.4472065 -23.0681378
[119,] 45.3289581 -25.4472065
[120,] -17.3466078 45.3289581
[121,] -12.2128149 -17.3466078
[122,] 62.8686143 -12.2128149
[123,] -20.6326255 62.8686143
[124,] -1.5683111 -20.6326255
[125,] -25.1008658 -1.5683111
[126,] -17.7715239 -25.1008658
[127,] -11.4389022 -17.7715239
[128,] 8.4712393 -11.4389022
[129,] -17.7715239 8.4712393
[130,] -9.9665871 -17.7715239
[131,] 36.6716285 -9.9665871
[132,] -17.7715239 36.6716285
[133,] -15.1299894 -17.7715239
[134,] -17.7715239 -15.1299894
[135,] -10.8318097 -17.7715239
[136,] -5.6420026 -10.8318097
[137,] -17.7715239 -5.6420026
[138,] -17.7715239 -17.7715239
[139,] -17.7715239 -17.7715239
[140,] -5.7245588 -17.7715239
[141,] -17.7715239 -5.7245588
[142,] -7.9957263 -17.7715239
[143,] -15.0852406 -7.9957263
[144,] -6.9403875 -15.0852406
[145,] -7.2042128 -6.9403875
[146,] -17.7715239 -7.2042128
[147,] -11.8784102 -17.7715239
[148,] -17.7715239 -11.8784102
[149,] -16.9844017 -17.7715239
[150,] -10.4788181 -16.9844017
[151,] 43.3483659 -10.4788181
[152,] -17.7715239 43.3483659
[153,] -9.0478457 -17.7715239
[154,] -14.9941017 -9.0478457
[155,] -17.7715239 -14.9941017
[156,] -17.7715239 -17.7715239
[157,] -16.7434713 -17.7715239
[158,] -13.5715006 -16.7434713
[159,] -17.7715239 -13.5715006
[160,] -17.7715239 -17.7715239
[161,] -17.7715239 -17.7715239
[162,] -17.7715239 -17.7715239
[163,] 5.9047961 -17.7715239
[164,] 29.1341866 5.9047961
[165,] -17.7715239 29.1341866
[166,] -17.7715239 -17.7715239
[167,] 2.0702748 -17.7715239
[168,] -17.7715239 2.0702748
[169,] -17.7715239 -17.7715239
[170,] -17.7715239 -17.7715239
[171,] -14.3519010 -17.7715239
[172,] -17.7715239 -14.3519010
[173,] -10.6347972 -17.7715239
[174,] -17.7715239 -10.6347972
[175,] -10.4746235 -17.7715239
[176,] -17.7715239 -10.4746235
[177,] -19.2098605 -17.7715239
[178,] -17.7715239 -19.2098605
[179,] -17.7715239 -17.7715239
[180,] 5.2809704 -17.7715239
[181,] -17.6701883 5.2809704
[182,] 18.0882501 -17.6701883
[183,] -16.0784351 18.0882501
[184,] -15.6042886 -16.0784351
[185,] -9.9109276 -15.6042886
[186,] -10.1173556 -9.9109276
[187,] -14.3090944 -10.1173556
[188,] -16.3061709 -14.3090944
[189,] -17.7715239 -16.3061709
[190,] -15.1534012 -17.7715239
[191,] -17.7715239 -15.1534012
[192,] -13.0911907 -17.7715239
[193,] -17.7715239 -13.0911907
[194,] -16.0293232 -17.7715239
[195,] -19.7177719 -16.0293232
[196,] 21.8425098 -19.7177719
[197,] -19.7177719 21.8425098
[198,] -19.7177719 -19.7177719
[199,] -10.6143344 -19.7177719
[200,] -19.7177719 -10.6143344
[201,] -19.7177719 -19.7177719
[202,] -19.7177719 -19.7177719
[203,] 9.8502926 -19.7177719
[204,] -16.7006246 9.8502926
[205,] -18.1582073 -16.7006246
[206,] -18.1228516 -18.1582073
[207,] -19.7177719 -18.1228516
[208,] 9.4228224 -19.7177719
[209,] -19.7177719 9.4228224
[210,] -17.7715239 -19.7177719
[211,] -17.7715239 -17.7715239
[212,] 13.6412053 -17.7715239
[213,] -1.4905944 13.6412053
[214,] -15.1539497 -1.4905944
[215,] -17.7715239 -15.1539497
[216,] -9.9175622 -17.7715239
[217,] 7.5620269 -9.9175622
[218,] -17.7715239 7.5620269
[219,] -3.1874263 -17.7715239
[220,] -15.6436606 -3.1874263
[221,] -17.7715239 -15.6436606
[222,] -15.7831400 -17.7715239
[223,] -17.7715239 -15.7831400
[224,] -15.9411677 -17.7715239
[225,] 6.0684338 -15.9411677
[226,] -17.7715239 6.0684338
[227,] -17.7715239 -17.7715239
[228,] -16.8886255 -17.7715239
[229,] -17.7715239 -16.8886255
[230,] -12.1957821 -17.7715239
[231,] -17.7715239 -12.1957821
[232,] -17.7715239 -17.7715239
[233,] -3.6549638 -17.7715239
[234,] -17.7715239 -3.6549638
[235,] -16.7942593 -17.7715239
[236,] -17.7715239 -16.7942593
[237,] -11.9837720 -17.7715239
[238,] -17.7715239 -11.9837720
[239,] -17.7715239 -17.7715239
[240,] -17.7715239 -17.7715239
[241,] -15.6645685 -17.7715239
[242,] -17.7715239 -15.6645685
[243,] -9.5755039 -17.7715239
[244,] -34.3674831 -9.5755039
[245,] -19.7177719 -34.3674831
[246,] -19.7177719 -19.7177719
[247,] -8.0941280 -19.7177719
[248,] -19.7177719 -8.0941280
[249,] -19.7177719 -19.7177719
[250,] -7.1416207 -19.7177719
[251,] -19.7177719 -7.1416207
[252,] -18.3085199 -19.7177719
[253,] 46.1144388 -18.3085199
[254,] 52.4579341 46.1144388
[255,] -17.7715239 52.4579341
[256,] -9.9019202 -17.7715239
[257,] 19.9366442 -9.9019202
[258,] -19.7177719 19.9366442
[259,] -17.7715239 -19.7177719
[260,] -19.7177719 -17.7715239
[261,] -19.7177719 -19.7177719
[262,] -12.1279946 -19.7177719
[263,] -4.1143907 -12.1279946
[264,] -19.7177719 -4.1143907
[265,] -17.7715239 -19.7177719
[266,] -19.7177719 -17.7715239
[267,] -14.7883112 -19.7177719
[268,] -19.7177719 -14.7883112
[269,] -16.3292489 -19.7177719
[270,] -19.7177719 -16.3292489
[271,] -19.7177719 -19.7177719
[272,] -19.7177719 -19.7177719
[273,] -6.6514034 -19.7177719
[274,] -17.7715239 -6.6514034
[275,] 9.6241665 -17.7715239
[276,] 11.5785054 9.6241665
[277,] -17.7715239 11.5785054
[278,] 54.0868392 -17.7715239
[279,] 6.0817440 54.0868392
[280,] -19.7177719 6.0817440
[281,] -17.7715239 -19.7177719
[282,] 54.1272101 -17.7715239
[283,] -17.7715239 54.1272101
[284,] -11.0490447 -17.7715239
[285,] -17.7715239 -11.0490447
[286,] -19.7177719 -17.7715239
[287,] -14.5025951 -19.7177719
[288,] -19.7177719 -14.5025951
[289,] 45.0836053 -19.7177719
[290,] -19.7177719 45.0836053
[291,] 4.2638542 -19.7177719
[292,] -7.9689040 4.2638542
[293,] -17.7715239 -7.9689040
[294,] -3.6482575 -17.7715239
[295,] -0.7725222 -3.6482575
[296,] -17.7715239 -0.7725222
[297,] -19.7177719 -17.7715239
[298,] -17.7715239 -19.7177719
[299,] -17.7715239 -17.7715239
[300,] -15.9525365 -17.7715239
[301,] -16.9726810 -15.9525365
[302,] -17.7715239 -16.9726810
[303,] -19.7177719 -17.7715239
[304,] 32.3100615 -19.7177719
[305,] -8.7910753 32.3100615
[306,] -19.7177719 -8.7910753
[307,] -9.4488344 -19.7177719
[308,] 4.0235829 -9.4488344
[309,] -19.7177719 4.0235829
[310,] -19.7177719 -19.7177719
[311,] 52.2140070 -19.7177719
[312,] -19.7177719 52.2140070
[313,] -13.3985432 -19.7177719
[314,] -13.8082140 -13.3985432
[315,] -6.7721475 -13.8082140
[316,] -16.6974894 -6.7721475
[317,] -3.4484015 -16.6974894
[318,] -11.7779166 -3.4484015
[319,] -17.7715239 -11.7779166
[320,] -17.7715239 -17.7715239
[321,] -17.7715239 -17.7715239
[322,] -17.7715239 -17.7715239
[323,] -19.7177719 -17.7715239
[324,] -17.7715239 -19.7177719
[325,] -15.9217150 -17.7715239
[326,] -19.7177719 -15.9217150
[327,] -19.7177719 -19.7177719
[328,] -17.7715239 -19.7177719
[329,] -17.7715239 -17.7715239
[330,] -17.7715239 -17.7715239
[331,] -17.7715239 -17.7715239
[332,] -17.7715239 -17.7715239
[333,] -13.2877113 -17.7715239
[334,] -15.7811006 -13.2877113
[335,] -17.7715239 -15.7811006
[336,] -11.4031594 -17.7715239
[337,] -17.7715239 -11.4031594
[338,] -17.0181298 -17.7715239
[339,] -5.2279308 -17.0181298
[340,] -5.7309948 -5.2279308
[341,] -17.7715239 -5.7309948
[342,] -17.7715239 -17.7715239
[343,] -17.7715239 -17.7715239
[344,] -17.7715239 -17.7715239
[345,] -17.7715239 -17.7715239
[346,] -11.9420184 -17.7715239
[347,] -17.7715239 -11.9420184
[348,] -17.7715239 -17.7715239
[349,] -17.7715239 -17.7715239
[350,] -5.0456600 -17.7715239
[351,] -8.3397571 -5.0456600
[352,] -9.3099043 -8.3397571
[353,] -17.2235769 -9.3099043
[354,] -17.7715239 -17.2235769
[355,] -17.7715239 -17.7715239
[356,] 23.5501756 -17.7715239
[357,] -17.7715239 23.5501756
[358,] -17.7715239 -17.7715239
[359,] -17.7715239 -17.7715239
[360,] -13.7791106 -17.7715239
[361,] -17.7715239 -13.7791106
[362,] -17.7715239 -17.7715239
[363,] -17.4618720 -17.7715239
[364,] -7.5620163 -17.4618720
[365,] -17.7715239 -7.5620163
[366,] -17.7715239 -17.7715239
[367,] -17.7715239 -17.7715239
[368,] -17.7715239 -17.7715239
[369,] -17.7715239 -17.7715239
[370,] -17.7715239 -17.7715239
[371,] 19.5599725 -17.7715239
[372,] -16.6507693 19.5599725
[373,] -19.5849049 -16.6507693
[374,] -17.7715239 -19.5849049
[375,] -17.7715239 -17.7715239
[376,] 7.2145218 -17.7715239
[377,] -10.7804364 7.2145218
[378,] -17.7715239 -10.7804364
[379,] -17.7715239 -17.7715239
[380,] -17.7715239 -17.7715239
[381,] 33.5696334 -17.7715239
[382,] 66.0097200 33.5696334
[383,] -17.7715239 66.0097200
[384,] -3.0892383 -17.7715239
[385,] 2.5070302 -3.0892383
[386,] -12.2396209 2.5070302
[387,] -15.2035233 -12.2396209
[388,] -15.1322970 -15.2035233
[389,] -4.8040735 -15.1322970
[390,] -11.2594462 -4.8040735
[391,] -12.2721206 -11.2594462
[392,] -15.0868689 -12.2721206
[393,] -17.4055818 -15.0868689
[394,] -40.0900733 -17.4055818
[395,] -7.6868146 -40.0900733
[396,] -6.0135264 -7.6868146
[397,] 64.4297376 -6.0135264
[398,] -13.6174969 64.4297376
[399,] 119.6534186 -13.6174969
[400,] -23.8350537 119.6534186
[401,] -15.4380758 -23.8350537
[402,] -13.3116905 -15.4380758
[403,] -8.5154968 -13.3116905
[404,] -16.2045102 -8.5154968
[405,] -16.7677685 -16.2045102
[406,] -8.8879719 -16.7677685
[407,] -12.7750886 -8.8879719
[408,] -8.7906419 -12.7750886
[409,] -23.3870447 -8.7906419
[410,] 30.7412973 -23.3870447
[411,] -39.7864245 30.7412973
[412,] -9.2689253 -39.7864245
[413,] -14.9699204 -9.2689253
[414,] -27.4902067 -14.9699204
[415,] -26.4069997 -27.4902067
[416,] 36.2550853 -26.4069997
[417,] -11.0480680 36.2550853
[418,] -3.7731238 -11.0480680
[419,] -21.3297682 -3.7731238
[420,] -26.2523906 -21.3297682
[421,] 25.1544932 -26.2523906
[422,] 1.9547430 25.1544932
[423,] -34.1412990 1.9547430
[424,] -6.9383195 -34.1412990
[425,] 41.4744236 -6.9383195
[426,] -11.1591966 41.4744236
[427,] -75.1129514 -11.1591966
[428,] -40.0773534 -75.1129514
[429,] -37.8889397 -40.0773534
[430,] 52.8080066 -37.8889397
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 78.6110629 -103.2385301
2 -193.0483656 78.6110629
3 29.3169664 -193.0483656
4 -118.2528983 29.3169664
5 -191.6364156 -118.2528983
6 157.8412239 -191.6364156
7 -53.3358500 157.8412239
8 176.3545702 -53.3358500
9 -16.1011360 176.3545702
10 -151.7357115 -16.1011360
11 11.1840122 -151.7357115
12 121.8020360 11.1840122
13 -55.3738445 121.8020360
14 296.9647588 -55.3738445
15 4.8001512 296.9647588
16 -47.6667778 4.8001512
17 319.6081477 -47.6667778
18 -17.9383429 319.6081477
19 209.4863290 -17.9383429
20 2.4800653 209.4863290
21 -94.9479725 2.4800653
22 -47.9857649 -94.9479725
23 145.1398071 -47.9857649
24 9.5256328 145.1398071
25 -62.4682582 9.5256328
26 64.4450331 -62.4682582
27 46.1139352 64.4450331
28 -27.5903855 46.1139352
29 249.6574688 -27.5903855
30 -25.1191149 249.6574688
31 -33.0432926 -25.1191149
32 -28.5724041 -33.0432926
33 -17.8856722 -28.5724041
34 -106.7450658 -17.8856722
35 -5.3742185 -106.7450658
36 -7.9860084 -5.3742185
37 4.4478303 -7.9860084
38 -8.9388896 4.4478303
39 92.0017448 -8.9388896
40 718.9778225 92.0017448
41 -64.1458767 718.9778225
42 164.3304000 -64.1458767
43 17.9768446 164.3304000
44 9.4424685 17.9768446
45 -36.8708145 9.4424685
46 -61.9475680 -36.8708145
47 221.0689149 -61.9475680
48 8.5897065 221.0689149
49 -6.5232862 8.5897065
50 47.9687929 -6.5232862
51 72.6333328 47.9687929
52 -5.1928587 72.6333328
53 100.5541410 -5.1928587
54 68.4877525 100.5541410
55 62.7895433 68.4877525
56 -30.0957386 62.7895433
57 1.0097033 -30.0957386
58 65.0723661 1.0097033
59 7.3193113 65.0723661
60 -23.8946707 7.3193113
61 10.0091245 -23.8946707
62 -20.4304021 10.0091245
63 5.9561988 -20.4304021
64 43.7411938 5.9561988
65 827.2481558 43.7411938
66 42.2644065 827.2481558
67 3.5391411 42.2644065
68 38.5408987 3.5391411
69 0.1345086 38.5408987
70 30.8938355 0.1345086
71 5.0960188 30.8938355
72 -11.8626854 5.0960188
73 -17.9224420 -11.8626854
74 -19.3312313 -17.9224420
75 -18.1779932 -19.3312313
76 0.8451727 -18.1779932
77 14.4366059 0.8451727
78 87.0558872 14.4366059
79 124.0143014 87.0558872
80 -1.2469427 124.0143014
81 10.2851672 -1.2469427
82 4.1358990 10.2851672
83 3.0966810 4.1358990
84 -22.9799469 3.0966810
85 12.2290963 -22.9799469
86 31.7390770 12.2290963
87 -24.6480353 31.7390770
88 13.7101904 -24.6480353
89 -27.5866134 13.7101904
90 6.7775631 -27.5866134
91 -7.5723418 6.7775631
92 64.7856116 -7.5723418
93 -15.4249588 64.7856116
94 23.6317018 -15.4249588
95 2.7249456 23.6317018
96 -32.4240669 2.7249456
97 164.9135253 -32.4240669
98 -40.8120314 164.9135253
99 -23.0846911 -40.8120314
100 -1.4076433 -23.0846911
101 -13.2584274 -1.4076433
102 20.3629300 -13.2584274
103 -24.7933030 20.3629300
104 -1.8146082 -24.7933030
105 -4.1520971 -1.8146082
106 18.9601902 -4.1520971
107 -28.3983981 18.9601902
108 -12.1869943 -28.3983981
109 -1.7065100 -12.1869943
110 2.9992310 -1.7065100
111 8.9964610 2.9992310
112 8.2798501 8.9964610
113 -27.2884462 8.2798501
114 -18.8804091 -27.2884462
115 12.7557697 -18.8804091
116 -4.1608512 12.7557697
117 -23.0681378 -4.1608512
118 -25.4472065 -23.0681378
119 45.3289581 -25.4472065
120 -17.3466078 45.3289581
121 -12.2128149 -17.3466078
122 62.8686143 -12.2128149
123 -20.6326255 62.8686143
124 -1.5683111 -20.6326255
125 -25.1008658 -1.5683111
126 -17.7715239 -25.1008658
127 -11.4389022 -17.7715239
128 8.4712393 -11.4389022
129 -17.7715239 8.4712393
130 -9.9665871 -17.7715239
131 36.6716285 -9.9665871
132 -17.7715239 36.6716285
133 -15.1299894 -17.7715239
134 -17.7715239 -15.1299894
135 -10.8318097 -17.7715239
136 -5.6420026 -10.8318097
137 -17.7715239 -5.6420026
138 -17.7715239 -17.7715239
139 -17.7715239 -17.7715239
140 -5.7245588 -17.7715239
141 -17.7715239 -5.7245588
142 -7.9957263 -17.7715239
143 -15.0852406 -7.9957263
144 -6.9403875 -15.0852406
145 -7.2042128 -6.9403875
146 -17.7715239 -7.2042128
147 -11.8784102 -17.7715239
148 -17.7715239 -11.8784102
149 -16.9844017 -17.7715239
150 -10.4788181 -16.9844017
151 43.3483659 -10.4788181
152 -17.7715239 43.3483659
153 -9.0478457 -17.7715239
154 -14.9941017 -9.0478457
155 -17.7715239 -14.9941017
156 -17.7715239 -17.7715239
157 -16.7434713 -17.7715239
158 -13.5715006 -16.7434713
159 -17.7715239 -13.5715006
160 -17.7715239 -17.7715239
161 -17.7715239 -17.7715239
162 -17.7715239 -17.7715239
163 5.9047961 -17.7715239
164 29.1341866 5.9047961
165 -17.7715239 29.1341866
166 -17.7715239 -17.7715239
167 2.0702748 -17.7715239
168 -17.7715239 2.0702748
169 -17.7715239 -17.7715239
170 -17.7715239 -17.7715239
171 -14.3519010 -17.7715239
172 -17.7715239 -14.3519010
173 -10.6347972 -17.7715239
174 -17.7715239 -10.6347972
175 -10.4746235 -17.7715239
176 -17.7715239 -10.4746235
177 -19.2098605 -17.7715239
178 -17.7715239 -19.2098605
179 -17.7715239 -17.7715239
180 5.2809704 -17.7715239
181 -17.6701883 5.2809704
182 18.0882501 -17.6701883
183 -16.0784351 18.0882501
184 -15.6042886 -16.0784351
185 -9.9109276 -15.6042886
186 -10.1173556 -9.9109276
187 -14.3090944 -10.1173556
188 -16.3061709 -14.3090944
189 -17.7715239 -16.3061709
190 -15.1534012 -17.7715239
191 -17.7715239 -15.1534012
192 -13.0911907 -17.7715239
193 -17.7715239 -13.0911907
194 -16.0293232 -17.7715239
195 -19.7177719 -16.0293232
196 21.8425098 -19.7177719
197 -19.7177719 21.8425098
198 -19.7177719 -19.7177719
199 -10.6143344 -19.7177719
200 -19.7177719 -10.6143344
201 -19.7177719 -19.7177719
202 -19.7177719 -19.7177719
203 9.8502926 -19.7177719
204 -16.7006246 9.8502926
205 -18.1582073 -16.7006246
206 -18.1228516 -18.1582073
207 -19.7177719 -18.1228516
208 9.4228224 -19.7177719
209 -19.7177719 9.4228224
210 -17.7715239 -19.7177719
211 -17.7715239 -17.7715239
212 13.6412053 -17.7715239
213 -1.4905944 13.6412053
214 -15.1539497 -1.4905944
215 -17.7715239 -15.1539497
216 -9.9175622 -17.7715239
217 7.5620269 -9.9175622
218 -17.7715239 7.5620269
219 -3.1874263 -17.7715239
220 -15.6436606 -3.1874263
221 -17.7715239 -15.6436606
222 -15.7831400 -17.7715239
223 -17.7715239 -15.7831400
224 -15.9411677 -17.7715239
225 6.0684338 -15.9411677
226 -17.7715239 6.0684338
227 -17.7715239 -17.7715239
228 -16.8886255 -17.7715239
229 -17.7715239 -16.8886255
230 -12.1957821 -17.7715239
231 -17.7715239 -12.1957821
232 -17.7715239 -17.7715239
233 -3.6549638 -17.7715239
234 -17.7715239 -3.6549638
235 -16.7942593 -17.7715239
236 -17.7715239 -16.7942593
237 -11.9837720 -17.7715239
238 -17.7715239 -11.9837720
239 -17.7715239 -17.7715239
240 -17.7715239 -17.7715239
241 -15.6645685 -17.7715239
242 -17.7715239 -15.6645685
243 -9.5755039 -17.7715239
244 -34.3674831 -9.5755039
245 -19.7177719 -34.3674831
246 -19.7177719 -19.7177719
247 -8.0941280 -19.7177719
248 -19.7177719 -8.0941280
249 -19.7177719 -19.7177719
250 -7.1416207 -19.7177719
251 -19.7177719 -7.1416207
252 -18.3085199 -19.7177719
253 46.1144388 -18.3085199
254 52.4579341 46.1144388
255 -17.7715239 52.4579341
256 -9.9019202 -17.7715239
257 19.9366442 -9.9019202
258 -19.7177719 19.9366442
259 -17.7715239 -19.7177719
260 -19.7177719 -17.7715239
261 -19.7177719 -19.7177719
262 -12.1279946 -19.7177719
263 -4.1143907 -12.1279946
264 -19.7177719 -4.1143907
265 -17.7715239 -19.7177719
266 -19.7177719 -17.7715239
267 -14.7883112 -19.7177719
268 -19.7177719 -14.7883112
269 -16.3292489 -19.7177719
270 -19.7177719 -16.3292489
271 -19.7177719 -19.7177719
272 -19.7177719 -19.7177719
273 -6.6514034 -19.7177719
274 -17.7715239 -6.6514034
275 9.6241665 -17.7715239
276 11.5785054 9.6241665
277 -17.7715239 11.5785054
278 54.0868392 -17.7715239
279 6.0817440 54.0868392
280 -19.7177719 6.0817440
281 -17.7715239 -19.7177719
282 54.1272101 -17.7715239
283 -17.7715239 54.1272101
284 -11.0490447 -17.7715239
285 -17.7715239 -11.0490447
286 -19.7177719 -17.7715239
287 -14.5025951 -19.7177719
288 -19.7177719 -14.5025951
289 45.0836053 -19.7177719
290 -19.7177719 45.0836053
291 4.2638542 -19.7177719
292 -7.9689040 4.2638542
293 -17.7715239 -7.9689040
294 -3.6482575 -17.7715239
295 -0.7725222 -3.6482575
296 -17.7715239 -0.7725222
297 -19.7177719 -17.7715239
298 -17.7715239 -19.7177719
299 -17.7715239 -17.7715239
300 -15.9525365 -17.7715239
301 -16.9726810 -15.9525365
302 -17.7715239 -16.9726810
303 -19.7177719 -17.7715239
304 32.3100615 -19.7177719
305 -8.7910753 32.3100615
306 -19.7177719 -8.7910753
307 -9.4488344 -19.7177719
308 4.0235829 -9.4488344
309 -19.7177719 4.0235829
310 -19.7177719 -19.7177719
311 52.2140070 -19.7177719
312 -19.7177719 52.2140070
313 -13.3985432 -19.7177719
314 -13.8082140 -13.3985432
315 -6.7721475 -13.8082140
316 -16.6974894 -6.7721475
317 -3.4484015 -16.6974894
318 -11.7779166 -3.4484015
319 -17.7715239 -11.7779166
320 -17.7715239 -17.7715239
321 -17.7715239 -17.7715239
322 -17.7715239 -17.7715239
323 -19.7177719 -17.7715239
324 -17.7715239 -19.7177719
325 -15.9217150 -17.7715239
326 -19.7177719 -15.9217150
327 -19.7177719 -19.7177719
328 -17.7715239 -19.7177719
329 -17.7715239 -17.7715239
330 -17.7715239 -17.7715239
331 -17.7715239 -17.7715239
332 -17.7715239 -17.7715239
333 -13.2877113 -17.7715239
334 -15.7811006 -13.2877113
335 -17.7715239 -15.7811006
336 -11.4031594 -17.7715239
337 -17.7715239 -11.4031594
338 -17.0181298 -17.7715239
339 -5.2279308 -17.0181298
340 -5.7309948 -5.2279308
341 -17.7715239 -5.7309948
342 -17.7715239 -17.7715239
343 -17.7715239 -17.7715239
344 -17.7715239 -17.7715239
345 -17.7715239 -17.7715239
346 -11.9420184 -17.7715239
347 -17.7715239 -11.9420184
348 -17.7715239 -17.7715239
349 -17.7715239 -17.7715239
350 -5.0456600 -17.7715239
351 -8.3397571 -5.0456600
352 -9.3099043 -8.3397571
353 -17.2235769 -9.3099043
354 -17.7715239 -17.2235769
355 -17.7715239 -17.7715239
356 23.5501756 -17.7715239
357 -17.7715239 23.5501756
358 -17.7715239 -17.7715239
359 -17.7715239 -17.7715239
360 -13.7791106 -17.7715239
361 -17.7715239 -13.7791106
362 -17.7715239 -17.7715239
363 -17.4618720 -17.7715239
364 -7.5620163 -17.4618720
365 -17.7715239 -7.5620163
366 -17.7715239 -17.7715239
367 -17.7715239 -17.7715239
368 -17.7715239 -17.7715239
369 -17.7715239 -17.7715239
370 -17.7715239 -17.7715239
371 19.5599725 -17.7715239
372 -16.6507693 19.5599725
373 -19.5849049 -16.6507693
374 -17.7715239 -19.5849049
375 -17.7715239 -17.7715239
376 7.2145218 -17.7715239
377 -10.7804364 7.2145218
378 -17.7715239 -10.7804364
379 -17.7715239 -17.7715239
380 -17.7715239 -17.7715239
381 33.5696334 -17.7715239
382 66.0097200 33.5696334
383 -17.7715239 66.0097200
384 -3.0892383 -17.7715239
385 2.5070302 -3.0892383
386 -12.2396209 2.5070302
387 -15.2035233 -12.2396209
388 -15.1322970 -15.2035233
389 -4.8040735 -15.1322970
390 -11.2594462 -4.8040735
391 -12.2721206 -11.2594462
392 -15.0868689 -12.2721206
393 -17.4055818 -15.0868689
394 -40.0900733 -17.4055818
395 -7.6868146 -40.0900733
396 -6.0135264 -7.6868146
397 64.4297376 -6.0135264
398 -13.6174969 64.4297376
399 119.6534186 -13.6174969
400 -23.8350537 119.6534186
401 -15.4380758 -23.8350537
402 -13.3116905 -15.4380758
403 -8.5154968 -13.3116905
404 -16.2045102 -8.5154968
405 -16.7677685 -16.2045102
406 -8.8879719 -16.7677685
407 -12.7750886 -8.8879719
408 -8.7906419 -12.7750886
409 -23.3870447 -8.7906419
410 30.7412973 -23.3870447
411 -39.7864245 30.7412973
412 -9.2689253 -39.7864245
413 -14.9699204 -9.2689253
414 -27.4902067 -14.9699204
415 -26.4069997 -27.4902067
416 36.2550853 -26.4069997
417 -11.0480680 36.2550853
418 -3.7731238 -11.0480680
419 -21.3297682 -3.7731238
420 -26.2523906 -21.3297682
421 25.1544932 -26.2523906
422 1.9547430 25.1544932
423 -34.1412990 1.9547430
424 -6.9383195 -34.1412990
425 41.4744236 -6.9383195
426 -11.1591966 41.4744236
427 -75.1129514 -11.1591966
428 -40.0773534 -75.1129514
429 -37.8889397 -40.0773534
430 52.8080066 -37.8889397
> 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/7emzj1291317431.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/87vgm1291317431.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/97vgm1291317431.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/10h4x71291317431.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/113nwv1291317431.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/1265cj1291317431.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/132xsa1291317431.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/14nfqx1291317431.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/159gpl1291317431.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/165q5u1291317431.tab")
+ }
>
> try(system("convert tmp/1b3iv1291317431.ps tmp/1b3iv1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/2b3iv1291317431.ps tmp/2b3iv1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/3lczy1291317431.ps tmp/3lczy1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/4lczy1291317431.ps tmp/4lczy1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/5lczy1291317431.ps tmp/5lczy1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/6emzj1291317431.ps tmp/6emzj1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/7emzj1291317431.ps tmp/7emzj1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/87vgm1291317431.ps tmp/87vgm1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/97vgm1291317431.ps tmp/97vgm1291317431.png",intern=TRUE))
character(0)
> try(system("convert tmp/10h4x71291317431.ps tmp/10h4x71291317431.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.288 3.100 13.792