R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(210907
+ ,56
+ ,79
+ ,30
+ ,112285
+ ,24188
+ ,120982
+ ,56
+ ,58
+ ,28
+ ,84786
+ ,18273
+ ,176508
+ ,54
+ ,60
+ ,38
+ ,83123
+ ,14130
+ ,179321
+ ,89
+ ,108
+ ,30
+ ,101193
+ ,32287
+ ,123185
+ ,40
+ ,49
+ ,22
+ ,38361
+ ,8654
+ ,52746
+ ,25
+ ,0
+ ,26
+ ,68504
+ ,9245
+ ,385534
+ ,92
+ ,121
+ ,25
+ ,119182
+ ,33251
+ ,33170
+ ,18
+ ,1
+ ,18
+ ,22807
+ ,1271
+ ,101645
+ ,63
+ ,20
+ ,11
+ ,17140
+ ,5279
+ ,149061
+ ,44
+ ,43
+ ,26
+ ,116174
+ ,27101
+ ,165446
+ ,33
+ ,69
+ ,25
+ ,57635
+ ,16373
+ ,237213
+ ,84
+ ,78
+ ,38
+ ,66198
+ ,19716
+ ,173326
+ ,88
+ ,86
+ ,44
+ ,71701
+ ,17753
+ ,133131
+ ,55
+ ,44
+ ,30
+ ,57793
+ ,9028
+ ,258873
+ ,60
+ ,104
+ ,40
+ ,80444
+ ,18653
+ ,180083
+ ,66
+ ,63
+ ,34
+ ,53855
+ ,8828
+ ,324799
+ ,154
+ ,158
+ ,47
+ ,97668
+ ,29498
+ ,230964
+ ,53
+ ,102
+ ,30
+ ,133824
+ ,27563
+ ,236785
+ ,119
+ ,77
+ ,31
+ ,101481
+ ,18293
+ ,135473
+ ,41
+ ,82
+ ,23
+ ,99645
+ ,22530
+ ,202925
+ ,61
+ ,115
+ ,36
+ ,114789
+ ,15977
+ ,215147
+ ,58
+ ,101
+ ,36
+ ,99052
+ ,35082
+ ,344297
+ ,75
+ ,80
+ ,30
+ ,67654
+ ,16116
+ ,153935
+ ,33
+ ,50
+ ,25
+ ,65553
+ ,15849
+ ,132943
+ ,40
+ ,83
+ ,39
+ ,97500
+ ,16026
+ ,174724
+ ,92
+ ,123
+ ,34
+ ,69112
+ ,26569
+ ,174415
+ ,100
+ ,73
+ ,31
+ ,82753
+ ,24785
+ ,225548
+ ,112
+ ,81
+ ,31
+ ,85323
+ ,17569
+ ,223632
+ ,73
+ ,105
+ ,33
+ ,72654
+ ,23825
+ ,124817
+ ,40
+ ,47
+ ,25
+ ,30727
+ ,7869
+ ,221698
+ ,45
+ ,105
+ ,33
+ ,77873
+ ,14975
+ ,210767
+ ,60
+ ,94
+ ,35
+ ,117478
+ ,37791
+ ,170266
+ ,62
+ ,44
+ ,42
+ ,74007
+ ,9605
+ ,260561
+ ,75
+ ,114
+ ,43
+ ,90183
+ ,27295
+ ,84853
+ ,31
+ ,38
+ ,30
+ ,61542
+ ,2746
+ ,294424
+ ,77
+ ,107
+ ,33
+ ,101494
+ ,34461
+ ,101011
+ ,34
+ ,30
+ ,13
+ ,27570
+ ,8098
+ ,215641
+ ,46
+ ,71
+ ,32
+ ,55813
+ ,4787
+ ,325107
+ ,99
+ ,84
+ ,36
+ ,79215
+ ,24919
+ ,7176
+ ,17
+ ,0
+ ,0
+ ,1423
+ ,603
+ ,167542
+ ,66
+ ,59
+ ,28
+ ,55461
+ ,16329
+ ,106408
+ ,30
+ ,33
+ ,14
+ ,31081
+ ,12558
+ ,96560
+ ,76
+ ,42
+ ,17
+ ,22996
+ ,7784
+ ,265769
+ ,146
+ ,96
+ ,32
+ ,83122
+ ,28522
+ ,269651
+ ,67
+ ,106
+ ,30
+ ,70106
+ ,22265
+ ,149112
+ ,56
+ ,56
+ ,35
+ ,60578
+ ,14459
+ ,175824
+ ,107
+ ,57
+ ,20
+ ,39992
+ ,14526
+ ,152871
+ ,58
+ ,59
+ ,28
+ ,79892
+ ,22240
+ ,111665
+ ,34
+ ,39
+ ,28
+ ,49810
+ ,11802
+ ,116408
+ ,61
+ ,34
+ ,39
+ ,71570
+ ,7623
+ ,362301
+ ,119
+ ,76
+ ,34
+ ,100708
+ ,11912
+ ,78800
+ ,42
+ ,20
+ ,26
+ ,33032
+ ,7935
+ ,183167
+ ,66
+ ,91
+ ,39
+ ,82875
+ ,18220
+ ,277965
+ ,89
+ ,115
+ ,39
+ ,139077
+ ,19199
+ ,150629
+ ,44
+ ,85
+ ,33
+ ,71595
+ ,19918
+ ,168809
+ ,66
+ ,76
+ ,28
+ ,72260
+ ,21884
+ ,24188
+ ,24
+ ,8
+ ,4
+ ,5950
+ ,2694
+ ,329267
+ ,259
+ ,79
+ ,39
+ ,115762
+ ,15808
+ ,65029
+ ,17
+ ,21
+ ,18
+ ,32551
+ ,3597
+ ,101097
+ ,64
+ ,30
+ ,14
+ ,31701
+ ,5296
+ ,218946
+ ,41
+ ,76
+ ,29
+ ,80670
+ ,25239
+ ,244052
+ ,68
+ ,101
+ ,44
+ ,143558
+ ,29801
+ ,341570
+ ,168
+ ,94
+ ,21
+ ,117105
+ ,18450
+ ,103597
+ ,43
+ ,27
+ ,16
+ ,23789
+ ,7132
+ ,233328
+ ,132
+ ,92
+ ,28
+ ,120733
+ ,34861
+ ,256462
+ ,105
+ ,123
+ ,35
+ ,105195
+ ,35940
+ ,206161
+ ,71
+ ,75
+ ,28
+ ,73107
+ ,16688
+ ,311473
+ ,112
+ ,128
+ ,38
+ ,132068
+ ,24683
+ ,235800
+ ,94
+ ,105
+ ,23
+ ,149193
+ ,46230
+ ,177939
+ ,82
+ ,55
+ ,36
+ ,46821
+ ,10387
+ ,207176
+ ,70
+ ,56
+ ,32
+ ,87011
+ ,21436
+ ,196553
+ ,57
+ ,41
+ ,29
+ ,95260
+ ,30546
+ ,174184
+ ,53
+ ,72
+ ,25
+ ,55183
+ ,19746
+ ,143246
+ ,103
+ ,67
+ ,27
+ ,106671
+ ,15977
+ ,187559
+ ,121
+ ,75
+ ,36
+ ,73511
+ ,22583
+ ,187681
+ ,62
+ ,114
+ ,28
+ ,92945
+ ,17274
+ ,119016
+ ,52
+ ,118
+ ,23
+ ,78664
+ ,16469
+ ,182192
+ ,52
+ ,77
+ ,40
+ ,70054
+ ,14251
+ ,73566
+ ,32
+ ,22
+ ,23
+ ,22618
+ ,3007
+ ,194979
+ ,62
+ ,66
+ ,40
+ ,74011
+ ,16851
+ ,167488
+ ,45
+ ,69
+ ,28
+ ,83737
+ ,21113
+ ,143756
+ ,46
+ ,105
+ ,34
+ ,69094
+ ,17401
+ ,275541
+ ,63
+ ,116
+ ,33
+ ,93133
+ ,23958
+ ,243199
+ ,75
+ ,88
+ ,28
+ ,95536
+ ,23567
+ ,182999
+ ,88
+ ,73
+ ,34
+ ,225920
+ ,13065
+ ,135649
+ ,46
+ ,99
+ ,30
+ ,62133
+ ,15358
+ ,152299
+ ,53
+ ,62
+ ,33
+ ,61370
+ ,14587
+ ,120221
+ ,37
+ ,53
+ ,22
+ ,43836
+ ,12770
+ ,346485
+ ,90
+ ,118
+ ,38
+ ,106117
+ ,24021
+ ,145790
+ ,63
+ ,30
+ ,26
+ ,38692
+ ,9648
+ ,193339
+ ,78
+ ,100
+ ,35
+ ,84651
+ ,20537
+ ,80953
+ ,25
+ ,49
+ ,8
+ ,56622
+ ,7905
+ ,122774
+ ,45
+ ,24
+ ,24
+ ,15986
+ ,4527
+ ,130585
+ ,46
+ ,67
+ ,29
+ ,95364
+ ,30495
+ ,112611
+ ,41
+ ,46
+ ,20
+ ,26706
+ ,7117
+ ,286468
+ ,144
+ ,57
+ ,29
+ ,89691
+ ,17719
+ ,241066
+ ,82
+ ,75
+ ,45
+ ,67267
+ ,27056
+ ,148446
+ ,91
+ ,135
+ ,37
+ ,126846
+ ,33473
+ ,204713
+ ,71
+ ,68
+ ,33
+ ,41140
+ ,9758
+ ,182079
+ ,63
+ ,124
+ ,33
+ ,102860
+ ,21115
+ ,140344
+ ,53
+ ,33
+ ,25
+ ,51715
+ ,7236
+ ,220516
+ ,62
+ ,98
+ ,32
+ ,55801
+ ,13790
+ ,243060
+ ,63
+ ,58
+ ,29
+ ,111813
+ ,32902
+ ,162765
+ ,32
+ ,68
+ ,28
+ ,120293
+ ,25131
+ ,182613
+ ,39
+ ,81
+ ,28
+ ,138599
+ ,30910
+ ,232138
+ ,62
+ ,131
+ ,31
+ ,161647
+ ,35947
+ ,265318
+ ,117
+ ,110
+ ,52
+ ,115929
+ ,29848
+ ,85574
+ ,34
+ ,37
+ ,21
+ ,24266
+ ,6943
+ ,310839
+ ,92
+ ,130
+ ,24
+ ,162901
+ ,42705
+ ,225060
+ ,93
+ ,93
+ ,41
+ ,109825
+ ,31808
+ ,232317
+ ,54
+ ,118
+ ,33
+ ,129838
+ ,26675
+ ,144966
+ ,144
+ ,39
+ ,32
+ ,37510
+ ,8435
+ ,43287
+ ,14
+ ,13
+ ,19
+ ,43750
+ ,7409
+ ,155754
+ ,61
+ ,74
+ ,20
+ ,40652
+ ,14993
+ ,164709
+ ,109
+ ,81
+ ,31
+ ,87771
+ ,36867
+ ,201940
+ ,38
+ ,109
+ ,31
+ ,85872
+ ,33835
+ ,235454
+ ,73
+ ,151
+ ,32
+ ,89275
+ ,24164
+ ,220801
+ ,75
+ ,51
+ ,18
+ ,44418
+ ,12607
+ ,99466
+ ,50
+ ,28
+ ,23
+ ,192565
+ ,22609
+ ,92661
+ ,61
+ ,40
+ ,17
+ ,35232
+ ,5892
+ ,133328
+ ,55
+ ,56
+ ,20
+ ,40909
+ ,17014
+ ,61361
+ ,77
+ ,27
+ ,12
+ ,13294
+ ,5394
+ ,125930
+ ,75
+ ,37
+ ,17
+ ,32387
+ ,9178
+ ,100750
+ ,72
+ ,83
+ ,30
+ ,140867
+ ,6440
+ ,224549
+ ,50
+ ,54
+ ,31
+ ,120662
+ ,21916
+ ,82316
+ ,32
+ ,27
+ ,10
+ ,21233
+ ,4011
+ ,102010
+ ,53
+ ,28
+ ,13
+ ,44332
+ ,5818
+ ,101523
+ ,42
+ ,59
+ ,22
+ ,61056
+ ,18647
+ ,243511
+ ,71
+ ,133
+ ,42
+ ,101338
+ ,20556
+ ,22938
+ ,10
+ ,12
+ ,1
+ ,1168
+ ,238
+ ,41566
+ ,35
+ ,0
+ ,9
+ ,13497
+ ,70
+ ,152474
+ ,65
+ ,106
+ ,32
+ ,65567
+ ,22392
+ ,61857
+ ,25
+ ,23
+ ,11
+ ,25162
+ ,3913
+ ,99923
+ ,66
+ ,44
+ ,25
+ ,32334
+ ,12237
+ ,132487
+ ,41
+ ,71
+ ,36
+ ,40735
+ ,8388
+ ,317394
+ ,86
+ ,116
+ ,31
+ ,91413
+ ,22120
+ ,21054
+ ,16
+ ,4
+ ,0
+ ,855
+ ,338
+ ,209641
+ ,42
+ ,62
+ ,24
+ ,97068
+ ,11727
+ ,22648
+ ,19
+ ,12
+ ,13
+ ,44339
+ ,3704
+ ,31414
+ ,19
+ ,18
+ ,8
+ ,14116
+ ,3988
+ ,46698
+ ,45
+ ,14
+ ,13
+ ,10288
+ ,3030
+ ,131698
+ ,65
+ ,60
+ ,19
+ ,65622
+ ,13520
+ ,91735
+ ,35
+ ,7
+ ,18
+ ,16563
+ ,1421
+ ,244749
+ ,95
+ ,98
+ ,33
+ ,76643
+ ,20923
+ ,184510
+ ,49
+ ,64
+ ,40
+ ,110681
+ ,20237
+ ,79863
+ ,37
+ ,29
+ ,22
+ ,29011
+ ,3219
+ ,128423
+ ,64
+ ,32
+ ,38
+ ,92696
+ ,3769
+ ,97839
+ ,38
+ ,25
+ ,24
+ ,94785
+ ,12252
+ ,38214
+ ,34
+ ,16
+ ,8
+ ,8773
+ ,1888
+ ,151101
+ ,32
+ ,48
+ ,35
+ ,83209
+ ,14497
+ ,272458
+ ,65
+ ,100
+ ,43
+ ,93815
+ ,28864
+ ,172494
+ ,52
+ ,46
+ ,43
+ ,86687
+ ,21721
+ ,108043
+ ,62
+ ,45
+ ,14
+ ,34553
+ ,4821
+ ,328107
+ ,65
+ ,129
+ ,41
+ ,105547
+ ,33644
+ ,250579
+ ,83
+ ,130
+ ,38
+ ,103487
+ ,15923
+ ,351067
+ ,95
+ ,136
+ ,45
+ ,213688
+ ,42935
+ ,158015
+ ,29
+ ,59
+ ,31
+ ,71220
+ ,18864
+ ,98866
+ ,18
+ ,25
+ ,13
+ ,23517
+ ,4977
+ ,85439
+ ,33
+ ,32
+ ,28
+ ,56926
+ ,7785
+ ,229242
+ ,247
+ ,63
+ ,31
+ ,91721
+ ,17939
+ ,351619
+ ,139
+ ,95
+ ,40
+ ,115168
+ ,23436
+ ,84207
+ ,29
+ ,14
+ ,30
+ ,111194
+ ,325
+ ,120445
+ ,118
+ ,36
+ ,16
+ ,51009
+ ,13539
+ ,324598
+ ,110
+ ,113
+ ,37
+ ,135777
+ ,34538
+ ,131069
+ ,67
+ ,47
+ ,30
+ ,51513
+ ,12198
+ ,204271
+ ,42
+ ,92
+ ,35
+ ,74163
+ ,26924
+ ,165543
+ ,65
+ ,70
+ ,32
+ ,51633
+ ,12716
+ ,141722
+ ,94
+ ,19
+ ,27
+ ,75345
+ ,8172
+ ,116048
+ ,64
+ ,50
+ ,20
+ ,33416
+ ,10855
+ ,250047
+ ,81
+ ,41
+ ,18
+ ,83305
+ ,11932
+ ,299775
+ ,95
+ ,91
+ ,31
+ ,98952
+ ,14300
+ ,195838
+ ,67
+ ,111
+ ,31
+ ,102372
+ ,25515
+ ,173260
+ ,63
+ ,41
+ ,21
+ ,37238
+ ,2805
+ ,254488
+ ,83
+ ,120
+ ,39
+ ,103772
+ ,29402
+ ,104389
+ ,45
+ ,135
+ ,41
+ ,123969
+ ,16440
+ ,136084
+ ,30
+ ,27
+ ,13
+ ,27142
+ ,11221
+ ,199476
+ ,70
+ ,87
+ ,32
+ ,135400
+ ,28732
+ ,92499
+ ,32
+ ,25
+ ,18
+ ,21399
+ ,5250
+ ,224330
+ ,83
+ ,131
+ ,39
+ ,130115
+ ,28608
+ ,135781
+ ,31
+ ,45
+ ,14
+ ,24874
+ ,8092
+ ,74408
+ ,67
+ ,29
+ ,7
+ ,34988
+ ,4473
+ ,81240
+ ,66
+ ,58
+ ,17
+ ,45549
+ ,1572
+ ,14688
+ ,10
+ ,4
+ ,0
+ ,6023
+ ,2065
+ ,181633
+ ,70
+ ,47
+ ,30
+ ,64466
+ ,14817
+ ,271856
+ ,103
+ ,109
+ ,37
+ ,54990
+ ,16714
+ ,7199
+ ,5
+ ,7
+ ,0
+ ,1644
+ ,556
+ ,46660
+ ,20
+ ,12
+ ,5
+ ,6179
+ ,2089
+ ,17547
+ ,5
+ ,0
+ ,1
+ ,3926
+ ,2658
+ ,133368
+ ,36
+ ,37
+ ,16
+ ,32755
+ ,10695
+ ,95227
+ ,34
+ ,37
+ ,32
+ ,34777
+ ,1669
+ ,152601
+ ,48
+ ,46
+ ,24
+ ,73224
+ ,16267
+ ,98146
+ ,40
+ ,15
+ ,17
+ ,27114
+ ,7768
+ ,79619
+ ,43
+ ,42
+ ,11
+ ,20760
+ ,7252
+ ,59194
+ ,31
+ ,7
+ ,24
+ ,37636
+ ,6387
+ ,139942
+ ,42
+ ,54
+ ,22
+ ,65461
+ ,18715
+ ,118612
+ ,46
+ ,54
+ ,12
+ ,30080
+ ,7936
+ ,72880
+ ,33
+ ,14
+ ,19
+ ,24094
+ ,8643
+ ,65475
+ ,18
+ ,16
+ ,13
+ ,69008
+ ,7294
+ ,99643
+ ,55
+ ,33
+ ,17
+ ,54968
+ ,4570
+ ,71965
+ ,35
+ ,32
+ ,15
+ ,46090
+ ,7185
+ ,77272
+ ,59
+ ,21
+ ,16
+ ,27507
+ ,10058
+ ,49289
+ ,19
+ ,15
+ ,24
+ ,10672
+ ,2342
+ ,135131
+ ,66
+ ,38
+ ,15
+ ,34029
+ ,8509
+ ,108446
+ ,60
+ ,22
+ ,17
+ ,46300
+ ,13275
+ ,89746
+ ,36
+ ,28
+ ,18
+ ,24760
+ ,6816
+ ,44296
+ ,25
+ ,10
+ ,20
+ ,18779
+ ,1930
+ ,77648
+ ,47
+ ,31
+ ,16
+ ,21280
+ ,8086
+ ,181528
+ ,54
+ ,32
+ ,16
+ ,40662
+ ,10737
+ ,134019
+ ,53
+ ,32
+ ,18
+ ,28987
+ ,8033
+ ,124064
+ ,40
+ ,43
+ ,22
+ ,22827
+ ,7058
+ ,92630
+ ,40
+ ,27
+ ,8
+ ,18513
+ ,6782
+ ,121848
+ ,39
+ ,37
+ ,17
+ ,30594
+ ,5401
+ ,52915
+ ,14
+ ,20
+ ,18
+ ,24006
+ ,6521
+ ,81872
+ ,45
+ ,32
+ ,16
+ ,27913
+ ,10856
+ ,58981
+ ,36
+ ,0
+ ,23
+ ,42744
+ ,2154
+ ,53515
+ ,28
+ ,5
+ ,22
+ ,12934
+ ,6117
+ ,60812
+ ,44
+ ,26
+ ,13
+ ,22574
+ ,5238
+ ,56375
+ ,30
+ ,10
+ ,13
+ ,41385
+ ,4820
+ ,65490
+ ,22
+ ,27
+ ,16
+ ,18653
+ ,5615
+ ,80949
+ ,17
+ ,11
+ ,16
+ ,18472
+ ,4272
+ ,76302
+ ,31
+ ,29
+ ,20
+ ,30976
+ ,8702
+ ,104011
+ ,55
+ ,25
+ ,22
+ ,63339
+ ,15340
+ ,98104
+ ,54
+ ,55
+ ,17
+ ,25568
+ ,8030
+ ,67989
+ ,21
+ ,23
+ ,18
+ ,33747
+ ,9526
+ ,30989
+ ,14
+ ,5
+ ,17
+ ,4154
+ ,1278
+ ,135458
+ ,81
+ ,43
+ ,12
+ ,19474
+ ,4236
+ ,73504
+ ,35
+ ,23
+ ,7
+ ,35130
+ ,3023
+ ,63123
+ ,43
+ ,34
+ ,17
+ ,39067
+ ,7196
+ ,61254
+ ,46
+ ,36
+ ,14
+ ,13310
+ ,3394
+ ,74914
+ ,30
+ ,35
+ ,23
+ ,65892
+ ,6371
+ ,31774
+ ,23
+ ,0
+ ,17
+ ,4143
+ ,1574
+ ,81437
+ ,38
+ ,37
+ ,14
+ ,28579
+ ,9620
+ ,87186
+ ,54
+ ,28
+ ,15
+ ,51776
+ ,6978
+ ,50090
+ ,20
+ ,16
+ ,17
+ ,21152
+ ,4911
+ ,65745
+ ,53
+ ,26
+ ,21
+ ,38084
+ ,8645
+ ,56653
+ ,45
+ ,38
+ ,18
+ ,27717
+ ,8987
+ ,158399
+ ,39
+ ,23
+ ,18
+ ,32928
+ ,5544
+ ,46455
+ ,20
+ ,22
+ ,17
+ ,11342
+ ,3083
+ ,73624
+ ,24
+ ,30
+ ,17
+ ,19499
+ ,6909
+ ,38395
+ ,31
+ ,16
+ ,16
+ ,16380
+ ,3189
+ ,91899
+ ,35
+ ,18
+ ,15
+ ,36874
+ ,6745
+ ,139526
+ ,151
+ ,28
+ ,21
+ ,48259
+ ,16724
+ ,52164
+ ,52
+ ,32
+ ,16
+ ,16734
+ ,4850
+ ,51567
+ ,30
+ ,21
+ ,14
+ ,28207
+ ,7025
+ ,70551
+ ,31
+ ,23
+ ,15
+ ,30143
+ ,6047
+ ,84856
+ ,29
+ ,29
+ ,17
+ ,41369
+ ,7377
+ ,102538
+ ,57
+ ,50
+ ,15
+ ,45833
+ ,9078
+ ,86678
+ ,40
+ ,12
+ ,15
+ ,29156
+ ,4605
+ ,85709
+ ,44
+ ,21
+ ,10
+ ,35944
+ ,3238
+ ,34662
+ ,25
+ ,18
+ ,6
+ ,36278
+ ,8100
+ ,150580
+ ,77
+ ,27
+ ,22
+ ,45588
+ ,9653
+ ,99611
+ ,35
+ ,41
+ ,21
+ ,45097
+ ,8914
+ ,19349
+ ,11
+ ,13
+ ,1
+ ,3895
+ ,786
+ ,99373
+ ,63
+ ,12
+ ,18
+ ,28394
+ ,6700
+ ,86230
+ ,44
+ ,21
+ ,17
+ ,18632
+ ,5788
+ ,30837
+ ,19
+ ,8
+ ,4
+ ,2325
+ ,593
+ ,31706
+ ,13
+ ,26
+ ,10
+ ,25139
+ ,4506
+ ,89806
+ ,42
+ ,27
+ ,16
+ ,27975
+ ,6382
+ ,62088
+ ,38
+ ,13
+ ,16
+ ,14483
+ ,5621
+ ,40151
+ ,29
+ ,16
+ ,9
+ ,13127
+ ,3997
+ ,27634
+ ,20
+ ,2
+ ,16
+ ,5839
+ ,520
+ ,76990
+ ,27
+ ,42
+ ,17
+ ,24069
+ ,8891
+ ,37460
+ ,20
+ ,5
+ ,7
+ ,3738
+ ,999
+ ,54157
+ ,19
+ ,37
+ ,15
+ ,18625
+ ,7067
+ ,49862
+ ,37
+ ,17
+ ,14
+ ,36341
+ ,4639
+ ,84337
+ ,26
+ ,38
+ ,14
+ ,24548
+ ,5654
+ ,64175
+ ,42
+ ,37
+ ,18
+ ,21792
+ ,6928
+ ,59382
+ ,49
+ ,29
+ ,12
+ ,26263
+ ,1514
+ ,119308
+ ,30
+ ,32
+ ,16
+ ,23686
+ ,9238
+ ,76702
+ ,49
+ ,35
+ ,21
+ ,49303
+ ,8204
+ ,103425
+ ,67
+ ,17
+ ,19
+ ,25659
+ ,5926
+ ,70344
+ ,28
+ ,20
+ ,16
+ ,28904
+ ,5785
+ ,43410
+ ,19
+ ,7
+ ,1
+ ,2781
+ ,4
+ ,104838
+ ,49
+ ,46
+ ,16
+ ,29236
+ ,5930
+ ,62215
+ ,27
+ ,24
+ ,10
+ ,19546
+ ,3710
+ ,69304
+ ,30
+ ,40
+ ,19
+ ,22818
+ ,705
+ ,53117
+ ,22
+ ,3
+ ,12
+ ,32689
+ ,443
+ ,19764
+ ,12
+ ,10
+ ,2
+ ,5752
+ ,2416
+ ,86680
+ ,31
+ ,37
+ ,14
+ ,22197
+ ,7747
+ ,84105
+ ,20
+ ,17
+ ,17
+ ,20055
+ ,5432
+ ,77945
+ ,20
+ ,28
+ ,19
+ ,25272
+ ,4913
+ ,89113
+ ,39
+ ,19
+ ,14
+ ,82206
+ ,2650
+ ,91005
+ ,29
+ ,29
+ ,11
+ ,32073
+ ,2370
+ ,40248
+ ,16
+ ,8
+ ,4
+ ,5444
+ ,775
+ ,64187
+ ,27
+ ,10
+ ,16
+ ,20154
+ ,5576
+ ,50857
+ ,21
+ ,15
+ ,20
+ ,36944
+ ,1352
+ ,56613
+ ,19
+ ,15
+ ,12
+ ,8019
+ ,3080
+ ,62792
+ ,35
+ ,28
+ ,15
+ ,30884
+ ,10205
+ ,72535
+ ,14
+ ,17
+ ,16
+ ,19540
+ ,6095)
+ ,dim=c(6
+ ,289)
+ ,dimnames=list(c('Time_in_RFC'
+ ,'Logins'
+ ,'blogged_computations'
+ ,'compendiums_reviewed'
+ ,'totsize'
+ ,'totrevisions')
+ ,1:289))
> y <- array(NA,dim=c(6,289),dimnames=list(c('Time_in_RFC','Logins','blogged_computations','compendiums_reviewed','totsize','totrevisions'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from 'package:base':
as.Date, 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
Time_in_RFC Logins blogged_computations compendiums_reviewed totsize
1 210907 56 79 30 112285
2 120982 56 58 28 84786
3 176508 54 60 38 83123
4 179321 89 108 30 101193
5 123185 40 49 22 38361
6 52746 25 0 26 68504
7 385534 92 121 25 119182
8 33170 18 1 18 22807
9 101645 63 20 11 17140
10 149061 44 43 26 116174
11 165446 33 69 25 57635
12 237213 84 78 38 66198
13 173326 88 86 44 71701
14 133131 55 44 30 57793
15 258873 60 104 40 80444
16 180083 66 63 34 53855
17 324799 154 158 47 97668
18 230964 53 102 30 133824
19 236785 119 77 31 101481
20 135473 41 82 23 99645
21 202925 61 115 36 114789
22 215147 58 101 36 99052
23 344297 75 80 30 67654
24 153935 33 50 25 65553
25 132943 40 83 39 97500
26 174724 92 123 34 69112
27 174415 100 73 31 82753
28 225548 112 81 31 85323
29 223632 73 105 33 72654
30 124817 40 47 25 30727
31 221698 45 105 33 77873
32 210767 60 94 35 117478
33 170266 62 44 42 74007
34 260561 75 114 43 90183
35 84853 31 38 30 61542
36 294424 77 107 33 101494
37 101011 34 30 13 27570
38 215641 46 71 32 55813
39 325107 99 84 36 79215
40 7176 17 0 0 1423
41 167542 66 59 28 55461
42 106408 30 33 14 31081
43 96560 76 42 17 22996
44 265769 146 96 32 83122
45 269651 67 106 30 70106
46 149112 56 56 35 60578
47 175824 107 57 20 39992
48 152871 58 59 28 79892
49 111665 34 39 28 49810
50 116408 61 34 39 71570
51 362301 119 76 34 100708
52 78800 42 20 26 33032
53 183167 66 91 39 82875
54 277965 89 115 39 139077
55 150629 44 85 33 71595
56 168809 66 76 28 72260
57 24188 24 8 4 5950
58 329267 259 79 39 115762
59 65029 17 21 18 32551
60 101097 64 30 14 31701
61 218946 41 76 29 80670
62 244052 68 101 44 143558
63 341570 168 94 21 117105
64 103597 43 27 16 23789
65 233328 132 92 28 120733
66 256462 105 123 35 105195
67 206161 71 75 28 73107
68 311473 112 128 38 132068
69 235800 94 105 23 149193
70 177939 82 55 36 46821
71 207176 70 56 32 87011
72 196553 57 41 29 95260
73 174184 53 72 25 55183
74 143246 103 67 27 106671
75 187559 121 75 36 73511
76 187681 62 114 28 92945
77 119016 52 118 23 78664
78 182192 52 77 40 70054
79 73566 32 22 23 22618
80 194979 62 66 40 74011
81 167488 45 69 28 83737
82 143756 46 105 34 69094
83 275541 63 116 33 93133
84 243199 75 88 28 95536
85 182999 88 73 34 225920
86 135649 46 99 30 62133
87 152299 53 62 33 61370
88 120221 37 53 22 43836
89 346485 90 118 38 106117
90 145790 63 30 26 38692
91 193339 78 100 35 84651
92 80953 25 49 8 56622
93 122774 45 24 24 15986
94 130585 46 67 29 95364
95 112611 41 46 20 26706
96 286468 144 57 29 89691
97 241066 82 75 45 67267
98 148446 91 135 37 126846
99 204713 71 68 33 41140
100 182079 63 124 33 102860
101 140344 53 33 25 51715
102 220516 62 98 32 55801
103 243060 63 58 29 111813
104 162765 32 68 28 120293
105 182613 39 81 28 138599
106 232138 62 131 31 161647
107 265318 117 110 52 115929
108 85574 34 37 21 24266
109 310839 92 130 24 162901
110 225060 93 93 41 109825
111 232317 54 118 33 129838
112 144966 144 39 32 37510
113 43287 14 13 19 43750
114 155754 61 74 20 40652
115 164709 109 81 31 87771
116 201940 38 109 31 85872
117 235454 73 151 32 89275
118 220801 75 51 18 44418
119 99466 50 28 23 192565
120 92661 61 40 17 35232
121 133328 55 56 20 40909
122 61361 77 27 12 13294
123 125930 75 37 17 32387
124 100750 72 83 30 140867
125 224549 50 54 31 120662
126 82316 32 27 10 21233
127 102010 53 28 13 44332
128 101523 42 59 22 61056
129 243511 71 133 42 101338
130 22938 10 12 1 1168
131 41566 35 0 9 13497
132 152474 65 106 32 65567
133 61857 25 23 11 25162
134 99923 66 44 25 32334
135 132487 41 71 36 40735
136 317394 86 116 31 91413
137 21054 16 4 0 855
138 209641 42 62 24 97068
139 22648 19 12 13 44339
140 31414 19 18 8 14116
141 46698 45 14 13 10288
142 131698 65 60 19 65622
143 91735 35 7 18 16563
144 244749 95 98 33 76643
145 184510 49 64 40 110681
146 79863 37 29 22 29011
147 128423 64 32 38 92696
148 97839 38 25 24 94785
149 38214 34 16 8 8773
150 151101 32 48 35 83209
151 272458 65 100 43 93815
152 172494 52 46 43 86687
153 108043 62 45 14 34553
154 328107 65 129 41 105547
155 250579 83 130 38 103487
156 351067 95 136 45 213688
157 158015 29 59 31 71220
158 98866 18 25 13 23517
159 85439 33 32 28 56926
160 229242 247 63 31 91721
161 351619 139 95 40 115168
162 84207 29 14 30 111194
163 120445 118 36 16 51009
164 324598 110 113 37 135777
165 131069 67 47 30 51513
166 204271 42 92 35 74163
167 165543 65 70 32 51633
168 141722 94 19 27 75345
169 116048 64 50 20 33416
170 250047 81 41 18 83305
171 299775 95 91 31 98952
172 195838 67 111 31 102372
173 173260 63 41 21 37238
174 254488 83 120 39 103772
175 104389 45 135 41 123969
176 136084 30 27 13 27142
177 199476 70 87 32 135400
178 92499 32 25 18 21399
179 224330 83 131 39 130115
180 135781 31 45 14 24874
181 74408 67 29 7 34988
182 81240 66 58 17 45549
183 14688 10 4 0 6023
184 181633 70 47 30 64466
185 271856 103 109 37 54990
186 7199 5 7 0 1644
187 46660 20 12 5 6179
188 17547 5 0 1 3926
189 133368 36 37 16 32755
190 95227 34 37 32 34777
191 152601 48 46 24 73224
192 98146 40 15 17 27114
193 79619 43 42 11 20760
194 59194 31 7 24 37636
195 139942 42 54 22 65461
196 118612 46 54 12 30080
197 72880 33 14 19 24094
198 65475 18 16 13 69008
199 99643 55 33 17 54968
200 71965 35 32 15 46090
201 77272 59 21 16 27507
202 49289 19 15 24 10672
203 135131 66 38 15 34029
204 108446 60 22 17 46300
205 89746 36 28 18 24760
206 44296 25 10 20 18779
207 77648 47 31 16 21280
208 181528 54 32 16 40662
209 134019 53 32 18 28987
210 124064 40 43 22 22827
211 92630 40 27 8 18513
212 121848 39 37 17 30594
213 52915 14 20 18 24006
214 81872 45 32 16 27913
215 58981 36 0 23 42744
216 53515 28 5 22 12934
217 60812 44 26 13 22574
218 56375 30 10 13 41385
219 65490 22 27 16 18653
220 80949 17 11 16 18472
221 76302 31 29 20 30976
222 104011 55 25 22 63339
223 98104 54 55 17 25568
224 67989 21 23 18 33747
225 30989 14 5 17 4154
226 135458 81 43 12 19474
227 73504 35 23 7 35130
228 63123 43 34 17 39067
229 61254 46 36 14 13310
230 74914 30 35 23 65892
231 31774 23 0 17 4143
232 81437 38 37 14 28579
233 87186 54 28 15 51776
234 50090 20 16 17 21152
235 65745 53 26 21 38084
236 56653 45 38 18 27717
237 158399 39 23 18 32928
238 46455 20 22 17 11342
239 73624 24 30 17 19499
240 38395 31 16 16 16380
241 91899 35 18 15 36874
242 139526 151 28 21 48259
243 52164 52 32 16 16734
244 51567 30 21 14 28207
245 70551 31 23 15 30143
246 84856 29 29 17 41369
247 102538 57 50 15 45833
248 86678 40 12 15 29156
249 85709 44 21 10 35944
250 34662 25 18 6 36278
251 150580 77 27 22 45588
252 99611 35 41 21 45097
253 19349 11 13 1 3895
254 99373 63 12 18 28394
255 86230 44 21 17 18632
256 30837 19 8 4 2325
257 31706 13 26 10 25139
258 89806 42 27 16 27975
259 62088 38 13 16 14483
260 40151 29 16 9 13127
261 27634 20 2 16 5839
262 76990 27 42 17 24069
263 37460 20 5 7 3738
264 54157 19 37 15 18625
265 49862 37 17 14 36341
266 84337 26 38 14 24548
267 64175 42 37 18 21792
268 59382 49 29 12 26263
269 119308 30 32 16 23686
270 76702 49 35 21 49303
271 103425 67 17 19 25659
272 70344 28 20 16 28904
273 43410 19 7 1 2781
274 104838 49 46 16 29236
275 62215 27 24 10 19546
276 69304 30 40 19 22818
277 53117 22 3 12 32689
278 19764 12 10 2 5752
279 86680 31 37 14 22197
280 84105 20 17 17 20055
281 77945 20 28 19 25272
282 89113 39 19 14 82206
283 91005 29 29 11 32073
284 40248 16 8 4 5444
285 64187 27 10 16 20154
286 50857 21 15 20 36944
287 56613 19 15 12 8019
288 62792 35 28 15 30884
289 72535 14 17 16 19540
totrevisions
1 24188
2 18273
3 14130
4 32287
5 8654
6 9245
7 33251
8 1271
9 5279
10 27101
11 16373
12 19716
13 17753
14 9028
15 18653
16 8828
17 29498
18 27563
19 18293
20 22530
21 15977
22 35082
23 16116
24 15849
25 16026
26 26569
27 24785
28 17569
29 23825
30 7869
31 14975
32 37791
33 9605
34 27295
35 2746
36 34461
37 8098
38 4787
39 24919
40 603
41 16329
42 12558
43 7784
44 28522
45 22265
46 14459
47 14526
48 22240
49 11802
50 7623
51 11912
52 7935
53 18220
54 19199
55 19918
56 21884
57 2694
58 15808
59 3597
60 5296
61 25239
62 29801
63 18450
64 7132
65 34861
66 35940
67 16688
68 24683
69 46230
70 10387
71 21436
72 30546
73 19746
74 15977
75 22583
76 17274
77 16469
78 14251
79 3007
80 16851
81 21113
82 17401
83 23958
84 23567
85 13065
86 15358
87 14587
88 12770
89 24021
90 9648
91 20537
92 7905
93 4527
94 30495
95 7117
96 17719
97 27056
98 33473
99 9758
100 21115
101 7236
102 13790
103 32902
104 25131
105 30910
106 35947
107 29848
108 6943
109 42705
110 31808
111 26675
112 8435
113 7409
114 14993
115 36867
116 33835
117 24164
118 12607
119 22609
120 5892
121 17014
122 5394
123 9178
124 6440
125 21916
126 4011
127 5818
128 18647
129 20556
130 238
131 70
132 22392
133 3913
134 12237
135 8388
136 22120
137 338
138 11727
139 3704
140 3988
141 3030
142 13520
143 1421
144 20923
145 20237
146 3219
147 3769
148 12252
149 1888
150 14497
151 28864
152 21721
153 4821
154 33644
155 15923
156 42935
157 18864
158 4977
159 7785
160 17939
161 23436
162 325
163 13539
164 34538
165 12198
166 26924
167 12716
168 8172
169 10855
170 11932
171 14300
172 25515
173 2805
174 29402
175 16440
176 11221
177 28732
178 5250
179 28608
180 8092
181 4473
182 1572
183 2065
184 14817
185 16714
186 556
187 2089
188 2658
189 10695
190 1669
191 16267
192 7768
193 7252
194 6387
195 18715
196 7936
197 8643
198 7294
199 4570
200 7185
201 10058
202 2342
203 8509
204 13275
205 6816
206 1930
207 8086
208 10737
209 8033
210 7058
211 6782
212 5401
213 6521
214 10856
215 2154
216 6117
217 5238
218 4820
219 5615
220 4272
221 8702
222 15340
223 8030
224 9526
225 1278
226 4236
227 3023
228 7196
229 3394
230 6371
231 1574
232 9620
233 6978
234 4911
235 8645
236 8987
237 5544
238 3083
239 6909
240 3189
241 6745
242 16724
243 4850
244 7025
245 6047
246 7377
247 9078
248 4605
249 3238
250 8100
251 9653
252 8914
253 786
254 6700
255 5788
256 593
257 4506
258 6382
259 5621
260 3997
261 520
262 8891
263 999
264 7067
265 4639
266 5654
267 6928
268 1514
269 9238
270 8204
271 5926
272 5785
273 4
274 5930
275 3710
276 705
277 443
278 2416
279 7747
280 5432
281 4913
282 2650
283 2370
284 775
285 5576
286 1352
287 3080
288 10205
289 6095
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Logins blogged_computations
1965.6875 747.2744 773.8105
compendiums_reviewed totsize totrevisions
1248.1845 0.1274 1.4121
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-135598 -18952 -2334 16248 155558
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.966e+03 5.418e+03 0.363 0.716997
Logins 7.473e+02 7.922e+01 9.433 < 2e-16 ***
blogged_computations 7.738e+02 1.176e+02 6.580 2.27e-10 ***
compendiums_reviewed 1.248e+03 3.396e+02 3.676 0.000284 ***
totsize 1.274e-01 9.914e-02 1.285 0.199808
totrevisions 1.412e+00 4.401e-01 3.209 0.001487 **
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 36260 on 283 degrees of freedom
Multiple R-squared: 0.8094, Adjusted R-squared: 0.8061
F-statistic: 240.4 on 5 and 283 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.7968652 4.062695e-01 2.031348e-01
[2,] 0.7499568 5.000863e-01 2.500432e-01
[3,] 0.7210613 5.578774e-01 2.789387e-01
[4,] 0.9476719 1.046562e-01 5.232808e-02
[5,] 0.9382944 1.234112e-01 6.170559e-02
[6,] 0.9041579 1.916843e-01 9.584214e-02
[7,] 0.8717375 2.565249e-01 1.282625e-01
[8,] 0.8188219 3.623561e-01 1.811781e-01
[9,] 0.8205687 3.588626e-01 1.794313e-01
[10,] 0.8429308 3.141384e-01 1.570692e-01
[11,] 0.7903536 4.192928e-01 2.096464e-01
[12,] 0.8887244 2.225513e-01 1.112756e-01
[13,] 0.8992042 2.015916e-01 1.007958e-01
[14,] 0.8665257 2.669486e-01 1.334743e-01
[15,] 0.9984846 3.030828e-03 1.515414e-03
[16,] 0.9977646 4.470735e-03 2.235367e-03
[17,] 0.9976787 4.642608e-03 2.321304e-03
[18,] 0.9996084 7.832273e-04 3.916136e-04
[19,] 0.9995641 8.718441e-04 4.359221e-04
[20,] 0.9992921 1.415845e-03 7.079225e-04
[21,] 0.9988633 2.273318e-03 1.136659e-03
[22,] 0.9982249 3.550104e-03 1.775052e-03
[23,] 0.9974861 5.027899e-03 2.513950e-03
[24,] 0.9963190 7.362060e-03 3.681030e-03
[25,] 0.9956475 8.704912e-03 4.352456e-03
[26,] 0.9944573 1.108532e-02 5.542661e-03
[27,] 0.9933658 1.326845e-02 6.634225e-03
[28,] 0.9948372 1.032551e-02 5.162756e-03
[29,] 0.9927461 1.450783e-02 7.253917e-03
[30,] 0.9949573 1.008541e-02 5.042705e-03
[31,] 0.9992369 1.526145e-03 7.630727e-04
[32,] 0.9991119 1.776180e-03 8.880898e-04
[33,] 0.9986816 2.636710e-03 1.318355e-03
[34,] 0.9981201 3.759805e-03 1.879903e-03
[35,] 0.9981544 3.691299e-03 1.845649e-03
[36,] 0.9974257 5.148588e-03 2.574294e-03
[37,] 0.9977638 4.472372e-03 2.236186e-03
[38,] 0.9968638 6.272351e-03 3.136176e-03
[39,] 0.9956451 8.709712e-03 4.354856e-03
[40,] 0.9941993 1.160138e-02 5.800690e-03
[41,] 0.9920891 1.582177e-02 7.910885e-03
[42,] 0.9900161 1.996776e-02 9.983878e-03
[43,] 0.9995082 9.836308e-04 4.918154e-04
[44,] 0.9993051 1.389736e-03 6.948682e-04
[45,] 0.9991669 1.666179e-03 8.330893e-04
[46,] 0.9989627 2.074611e-03 1.037305e-03
[47,] 0.9987522 2.495603e-03 1.247802e-03
[48,] 0.9983612 3.277597e-03 1.638798e-03
[49,] 0.9979043 4.191418e-03 2.095709e-03
[50,] 0.9980086 3.982706e-03 1.991353e-03
[51,] 0.9972600 5.480006e-03 2.740003e-03
[52,] 0.9963654 7.269113e-03 3.634557e-03
[53,] 0.9968499 6.300129e-03 3.150064e-03
[54,] 0.9957794 8.441178e-03 4.220589e-03
[55,] 0.9968773 6.245311e-03 3.122655e-03
[56,] 0.9959417 8.116533e-03 4.058267e-03
[57,] 0.9960711 7.857713e-03 3.928856e-03
[58,] 0.9954420 9.116006e-03 4.558003e-03
[59,] 0.9943897 1.122050e-02 5.610252e-03
[60,] 0.9933603 1.327932e-02 6.639660e-03
[61,] 0.9925937 1.481257e-02 7.406285e-03
[62,] 0.9905472 1.890559e-02 9.452796e-03
[63,] 0.9902171 1.956575e-02 9.782877e-03
[64,] 0.9915646 1.687086e-02 8.435432e-03
[65,] 0.9892640 2.147208e-02 1.073604e-02
[66,] 0.9942192 1.156167e-02 5.780835e-03
[67,] 0.9946442 1.071153e-02 5.355763e-03
[68,] 0.9949564 1.008713e-02 5.043565e-03
[69,] 0.9986157 2.768639e-03 1.384319e-03
[70,] 0.9981312 3.737557e-03 1.868778e-03
[71,] 0.9975157 4.968669e-03 2.484335e-03
[72,] 0.9968117 6.376528e-03 3.188264e-03
[73,] 0.9957937 8.412652e-03 4.206326e-03
[74,] 0.9966685 6.662994e-03 3.331497e-03
[75,] 0.9972040 5.592066e-03 2.796033e-03
[76,] 0.9970725 5.855035e-03 2.927518e-03
[77,] 0.9973330 5.334078e-03 2.667039e-03
[78,] 0.9977296 4.540881e-03 2.270441e-03
[79,] 0.9970212 5.957507e-03 2.978753e-03
[80,] 0.9960924 7.815197e-03 3.907599e-03
[81,] 0.9989982 2.003567e-03 1.001784e-03
[82,] 0.9987835 2.433076e-03 1.216538e-03
[83,] 0.9986328 2.734397e-03 1.367198e-03
[84,] 0.9982123 3.575350e-03 1.787675e-03
[85,] 0.9980015 3.997092e-03 1.998546e-03
[86,] 0.9983315 3.336953e-03 1.668477e-03
[87,] 0.9978006 4.398857e-03 2.199429e-03
[88,] 0.9986123 2.775334e-03 1.387667e-03
[89,] 0.9982792 3.441588e-03 1.720794e-03
[90,] 0.9999666 6.686594e-05 3.343297e-05
[91,] 0.9999668 6.644342e-05 3.322171e-05
[92,] 0.9999727 5.462121e-05 2.731061e-05
[93,] 0.9999666 6.688767e-05 3.344384e-05
[94,] 0.9999627 7.463440e-05 3.731720e-05
[95,] 0.9999744 5.121011e-05 2.560506e-05
[96,] 0.9999634 7.329385e-05 3.664692e-05
[97,] 0.9999492 1.015232e-04 5.076159e-05
[98,] 0.9999419 1.162312e-04 5.811559e-05
[99,] 0.9999326 1.347802e-04 6.739012e-05
[100,] 0.9999084 1.832975e-04 9.164876e-05
[101,] 0.9998941 2.118472e-04 1.059236e-04
[102,] 0.9998773 2.453025e-04 1.226512e-04
[103,] 0.9998301 3.398000e-04 1.699000e-04
[104,] 0.9998791 2.417546e-04 1.208773e-04
[105,] 0.9998500 3.000975e-04 1.500487e-04
[106,] 0.9997921 4.157188e-04 2.078594e-04
[107,] 0.9999589 8.210900e-05 4.105450e-05
[108,] 0.9999493 1.014657e-04 5.073283e-05
[109,] 0.9999388 1.223802e-04 6.119011e-05
[110,] 0.9999814 3.714201e-05 1.857101e-05
[111,] 0.9999894 2.112209e-05 1.056105e-05
[112,] 0.9999866 2.682358e-05 1.341179e-05
[113,] 0.9999813 3.737580e-05 1.868790e-05
[114,] 0.9999846 3.073556e-05 1.536778e-05
[115,] 0.9999779 4.428027e-05 2.214014e-05
[116,] 0.9999957 8.693504e-06 4.346752e-06
[117,] 0.9999972 5.513012e-06 2.756506e-06
[118,] 0.9999961 7.785256e-06 3.892628e-06
[119,] 0.9999943 1.138015e-05 5.690077e-06
[120,] 0.9999954 9.113585e-06 4.556793e-06
[121,] 0.9999934 1.329137e-05 6.645684e-06
[122,] 0.9999902 1.953681e-05 9.768406e-06
[123,] 0.9999858 2.835722e-05 1.417861e-05
[124,] 0.9999938 1.246823e-05 6.234113e-06
[125,] 0.9999908 1.836808e-05 9.184040e-06
[126,] 0.9999914 1.729891e-05 8.649455e-06
[127,] 0.9999881 2.388525e-05 1.194263e-05
[128,] 0.9999974 5.291095e-06 2.645548e-06
[129,] 0.9999960 7.906311e-06 3.953156e-06
[130,] 0.9999987 2.540453e-06 1.270226e-06
[131,] 0.9999986 2.748928e-06 1.374464e-06
[132,] 0.9999981 3.738056e-06 1.869028e-06
[133,] 0.9999976 4.893225e-06 2.446612e-06
[134,] 0.9999967 6.612117e-06 3.306059e-06
[135,] 0.9999967 6.685715e-06 3.342857e-06
[136,] 0.9999954 9.200817e-06 4.600409e-06
[137,] 0.9999932 1.358633e-05 6.793163e-06
[138,] 0.9999901 1.980831e-05 9.904156e-06
[139,] 0.9999858 2.846772e-05 1.423386e-05
[140,] 0.9999809 3.815206e-05 1.907603e-05
[141,] 0.9999741 5.181808e-05 2.590904e-05
[142,] 0.9999641 7.179055e-05 3.589528e-05
[143,] 0.9999636 7.280746e-05 3.640373e-05
[144,] 0.9999482 1.035904e-04 5.179522e-05
[145,] 0.9999268 1.464876e-04 7.324378e-05
[146,] 0.9999664 6.722854e-05 3.361427e-05
[147,] 0.9999552 8.960723e-05 4.480361e-05
[148,] 0.9999441 1.118076e-04 5.590378e-05
[149,] 0.9999242 1.515954e-04 7.579770e-05
[150,] 0.9999284 1.431138e-04 7.155690e-05
[151,] 0.9999072 1.856619e-04 9.283096e-05
[152,] 0.9999856 2.885398e-05 1.442699e-05
[153,] 0.9999967 6.664677e-06 3.332339e-06
[154,] 0.9999950 1.002916e-05 5.014580e-06
[155,] 0.9999973 5.400000e-06 2.700000e-06
[156,] 0.9999978 4.428757e-06 2.214378e-06
[157,] 0.9999969 6.180246e-06 3.090123e-06
[158,] 0.9999961 7.801397e-06 3.900699e-06
[159,] 0.9999943 1.148599e-05 5.742996e-06
[160,] 0.9999916 1.674810e-05 8.374050e-06
[161,] 0.9999885 2.299288e-05 1.149644e-05
[162,] 0.9999997 5.332544e-07 2.666272e-07
[163,] 1.0000000 5.629114e-09 2.814557e-09
[164,] 1.0000000 8.477043e-09 4.238522e-09
[165,] 1.0000000 9.683100e-10 4.841550e-10
[166,] 1.0000000 1.477484e-09 7.387420e-10
[167,] 1.0000000 5.763139e-12 2.881569e-12
[168,] 1.0000000 1.376309e-12 6.881544e-13
[169,] 1.0000000 2.359827e-12 1.179913e-12
[170,] 1.0000000 3.791402e-12 1.895701e-12
[171,] 1.0000000 1.339504e-12 6.697522e-13
[172,] 1.0000000 6.507633e-13 3.253816e-13
[173,] 1.0000000 1.075577e-12 5.377887e-13
[174,] 1.0000000 5.256248e-13 2.628124e-13
[175,] 1.0000000 1.063153e-12 5.315767e-13
[176,] 1.0000000 1.156677e-12 5.783385e-13
[177,] 1.0000000 4.824021e-13 2.412011e-13
[178,] 1.0000000 8.784898e-13 4.392449e-13
[179,] 1.0000000 1.708711e-12 8.543553e-13
[180,] 1.0000000 3.453975e-12 1.726988e-12
[181,] 1.0000000 1.919904e-12 9.599518e-13
[182,] 1.0000000 3.800355e-12 1.900177e-12
[183,] 1.0000000 5.146239e-12 2.573120e-12
[184,] 1.0000000 6.661523e-12 3.330762e-12
[185,] 1.0000000 1.211079e-11 6.055397e-12
[186,] 1.0000000 2.205832e-11 1.102916e-11
[187,] 1.0000000 3.987826e-11 1.993913e-11
[188,] 1.0000000 6.511287e-11 3.255644e-11
[189,] 1.0000000 1.268905e-10 6.344525e-11
[190,] 1.0000000 2.449442e-10 1.224721e-10
[191,] 1.0000000 4.720088e-10 2.360044e-10
[192,] 1.0000000 7.655417e-10 3.827709e-10
[193,] 1.0000000 1.170058e-09 5.850292e-10
[194,] 1.0000000 2.125750e-09 1.062875e-09
[195,] 1.0000000 2.385524e-09 1.192762e-09
[196,] 1.0000000 4.415989e-09 2.207995e-09
[197,] 1.0000000 7.723862e-09 3.861931e-09
[198,] 1.0000000 1.278767e-08 6.393837e-09
[199,] 1.0000000 2.163240e-08 1.081620e-08
[200,] 1.0000000 1.749602e-10 8.748012e-11
[201,] 1.0000000 7.864515e-11 3.932258e-11
[202,] 1.0000000 6.269342e-11 3.134671e-11
[203,] 1.0000000 7.522485e-11 3.761243e-11
[204,] 1.0000000 3.711891e-11 1.855946e-11
[205,] 1.0000000 7.551356e-11 3.775678e-11
[206,] 1.0000000 1.527192e-10 7.635958e-11
[207,] 1.0000000 2.764360e-10 1.382180e-10
[208,] 1.0000000 5.375101e-10 2.687551e-10
[209,] 1.0000000 8.984009e-10 4.492004e-10
[210,] 1.0000000 1.720386e-09 8.601929e-10
[211,] 1.0000000 3.468116e-09 1.734058e-09
[212,] 1.0000000 3.206711e-09 1.603356e-09
[213,] 1.0000000 6.384085e-09 3.192043e-09
[214,] 1.0000000 1.250626e-08 6.253132e-09
[215,] 1.0000000 2.378564e-08 1.189282e-08
[216,] 1.0000000 4.653280e-08 2.326640e-08
[217,] 1.0000000 8.121466e-08 4.060733e-08
[218,] 1.0000000 4.881457e-08 2.440728e-08
[219,] 1.0000000 8.639401e-08 4.319701e-08
[220,] 1.0000000 9.891456e-08 4.945728e-08
[221,] 0.9999999 1.661763e-07 8.308813e-08
[222,] 0.9999999 1.902209e-07 9.511043e-08
[223,] 0.9999999 2.898647e-07 1.449324e-07
[224,] 0.9999997 5.557691e-07 2.778845e-07
[225,] 0.9999995 1.016345e-06 5.081727e-07
[226,] 0.9999992 1.672186e-06 8.360930e-07
[227,] 0.9999994 1.246148e-06 6.230741e-07
[228,] 0.9999996 8.245109e-07 4.122554e-07
[229,] 1.0000000 2.348694e-09 1.174347e-09
[230,] 1.0000000 4.307594e-09 2.153797e-09
[231,] 1.0000000 9.893922e-09 4.946961e-09
[232,] 1.0000000 9.136113e-09 4.568056e-09
[233,] 1.0000000 1.453797e-08 7.268983e-09
[234,] 1.0000000 1.106407e-08 5.532034e-09
[235,] 1.0000000 3.589153e-09 1.794577e-09
[236,] 1.0000000 4.943164e-09 2.471582e-09
[237,] 1.0000000 1.210848e-08 6.054239e-09
[238,] 1.0000000 2.814048e-08 1.407024e-08
[239,] 1.0000000 6.213686e-08 3.106843e-08
[240,] 0.9999999 1.155058e-07 5.775292e-08
[241,] 0.9999999 2.127522e-07 1.063761e-07
[242,] 0.9999999 1.769975e-07 8.849875e-08
[243,] 0.9999999 1.344593e-07 6.722966e-08
[244,] 0.9999998 3.245374e-07 1.622687e-07
[245,] 0.9999997 6.815226e-07 3.407613e-07
[246,] 0.9999993 1.400111e-06 7.000553e-07
[247,] 0.9999986 2.860627e-06 1.430313e-06
[248,] 0.9999967 6.526649e-06 3.263325e-06
[249,] 0.9999974 5.163969e-06 2.581984e-06
[250,] 0.9999948 1.035525e-05 5.177626e-06
[251,] 0.9999881 2.381667e-05 1.190833e-05
[252,] 0.9999816 3.689166e-05 1.844583e-05
[253,] 0.9999720 5.603311e-05 2.801656e-05
[254,] 0.9999417 1.166220e-04 5.831101e-05
[255,] 0.9998733 2.533848e-04 1.266924e-04
[256,] 0.9998573 2.854024e-04 1.427012e-04
[257,] 0.9998347 3.306835e-04 1.653418e-04
[258,] 0.9996434 7.131659e-04 3.565830e-04
[259,] 0.9996650 6.700225e-04 3.350112e-04
[260,] 0.9995086 9.828192e-04 4.914096e-04
[261,] 0.9998837 2.325148e-04 1.162574e-04
[262,] 0.9999260 1.479302e-04 7.396511e-05
[263,] 0.9997986 4.028539e-04 2.014269e-04
[264,] 0.9994362 1.127501e-03 5.637503e-04
[265,] 0.9986924 2.615117e-03 1.307559e-03
[266,] 0.9971478 5.704367e-03 2.852183e-03
[267,] 0.9926642 1.467156e-02 7.335782e-03
[268,] 0.9895876 2.082473e-02 1.041236e-02
[269,] 0.9742433 5.151349e-02 2.575675e-02
[270,] 0.9695100 6.097993e-02 3.048997e-02
[271,] 0.9230117 1.539767e-01 7.698833e-02
[272,] 0.8934338 2.131325e-01 1.065662e-01
> postscript(file="/var/wessaorg/rcomp/tmp/1jgzx1355344485.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/wessaorg/rcomp/tmp/2go0p1355344485.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/wessaorg/rcomp/tmp/3c60h1355344485.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/wessaorg/rcomp/tmp/405bj1355344485.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/wessaorg/rcomp/tmp/50ele1355344485.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 = 289
Frequency = 1
1 2 3 4 5
20054.92650 -39267.37959 9786.07784 -68655.39762 8843.53987
6 7 8 9 10
-22137.28872 127843.98692 -10188.29732 13756.38536 -4582.90784
11 12 13 14 15
23758.78903 28412.35229 -50072.40810 -1539.85864 45077.60109
16 17 18 19 20
18281.13180 -27272.33110 17045.97188 8855.12632 -33802.37833
21 22 23 24 25
-15733.65770 -15410.48865 155557.88519 26681.36411 -46872.01729
26 27 28 29 30
-79932.23768 -43003.01463 2834.79234 1774.39902 10359.74197
31 32 33 34 35
32596.60991 -20793.22019 12505.46924 10629.31176 -18846.97710
36 37 38 39 40
49335.80770 19249.21558 70547.51005 93944.95077 -8526.16917
41 42 43 44 45
5527.39732 17320.23035 -29839.61265 -10393.64593 57775.35474
46 47 48 49 50
-9856.94210 -778.84802 -14625.18748 -3847.88783 -26013.24960
51 52 53 54 55
140509.74081 -17893.97584 -23502.58651 26993.79371 -28429.27699
56 57 58 59 60
-16345.08969 -11457.85825 -13124.67022 2415.69969 -900.68655
61 62 63 64 65
45416.38129 -2176.44131 74138.47965 15532.46808 -38028.18836
66 67 68 69 70
-26987.29162 25273.92420 27651.92408 -30658.93521 6569.48965
71 72 73 74 75
28269.58468 28797.29729 10779.11246 -57387.34914 -49053.31498
76 77 78 79 80
-20014.17536 -75104.51786 2807.65547 -5172.47956 12458.12129
81 82 83 84 85
3069.88846 -49648.22246 49847.30173 36691.62088 -30886.05635
86 87 88 89 90
-44347.75590 -6856.27117 -1483.79093 91083.16917 22525.05152
91 92 93 94 95
-27767.63674 -5973.59560 30223.63795 -49010.98847 5995.42331
96 97 98 99 100
60141.51087 16842.72809 -135598.22114 36860.54496 -47029.72501
101 102 103 104 105
25225.44827 29861.25678 52229.94578 -1496.10635 -7431.58973
106 107 108 109 110
-27578.34681 -31022.97585 -9537.96730 28512.50807 -28451.70469
111 112 113 114 115
3287.99066 -51418.05450 -18952.06033 -372.56244 -83325.77429
116 117 118 119 120
-10181.43502 -23347.00057 77396.17203 -46699.03837 -19869.02250
121 122 123 124 125
-7273.00787 -43326.71348 981.73627 -83732.48715 58418.77084
126 127 128 129 130
14693.53160 8681.75671 -39054.08098 -8790.57889 2480.76056
131 132 133 134 135
393.62456 -60004.51985 950.32838 -38014.93214 -17026.97673
136 137 138 139 140
79824.01143 3450.44625 69430.05665 -29907.57182 -16094.02315
141 142 143 144 145
-21544.30255 -16437.32710 31613.86609 15457.85536 3797.88625
146 147 148 149 150
-7894.20940 -10693.41551 -11202.36786 -15309.30010 13320.18080
151 152 153 154 155
38154.02052 685.41720 -3759.88493 65614.14323 2892.95766
156 157 158 159 160
28848.67729 14317.45410 37853.32040 -19143.97638 -81762.28252
161 162 163 164 165
74574.96492 -2334.05559 -43144.93606 40737.61938 -18566.97229
166 167 168 169 170
8573.59275 -3639.21181 -30.14844 -16983.59576 105895.45902
171 172 173 174 175
84907.24986 -29855.16574 57572.58778 -5778.68567 -125853.79245
176 177 178 179 180
55277.15603 -19886.48620 14667.88417 -46683.56981 43757.65126
181 182 183 184 185
-19576.99277 -44168.96173 -1529.10063 24406.53971 31784.39311
186 187 188 189 190
-4914.33595 10484.98267 6343.09546 36622.51786 -7506.50900
191 192 193 194 195
16914.08240 19039.09175 -13595.31643 -15124.62639 2576.73039
196 197 198 199 200
10468.63810 -3569.40924 2358.93874 -3634.30126 -15658.28958
201 202 203 204 205
-22711.67328 -13105.37921 19366.28898 -1244.10864 3964.73629
206 207 208 209 210
-14171.29845 -17528.39711 74133.93679 30181.71646 18598.27738
211 212 213 214 215
17959.20140 29363.70716 -9723.09919 -17340.36256 -7082.32644
216 217 218 219 220
-10989.39647 -20652.05036 -4052.56661 -4085.20378 29410.70205
221 222 223 224 225
-12468.31457 -15592.00306 -22590.18346 -7685.98041 -8860.67173
226 227 228 229 230
16248.12282 10104.16354 -33643.24839 -26906.62737 -22653.20035
231 232 233 234 235
-11348.68093 -12256.56991 -11972.36757 -10051.15966 -39217.20067
236 237 238 239 240
-47034.33114 75000.57199 -14497.79299 -2950.46343 -25678.31287
241 242 243 244 245
16904.55674 -52921.73748 -42373.71045 -20055.50577 -3480.16592
246 247 248 249 250
1871.76317 -18094.35535 16595.33607 12979.44426 -23463.55757
251 252 253 254 255
23281.73988 -4780.77125 -3750.60904 5497.11791 3367.81366
256 257 258 259 260
2356.26284 -22141.09574 3014.53448 -8087.44399 -14417.03685
261 262 263 264 265
-12273.97483 -14493.11392 6055.51487 -21713.16311 -21563.12080
266 267 268 269 270
4951.03102 -32834.22855 -22102.84879 34128.13524 -33041.99109
271 272 273 274 275
2884.23135 155.73900 20221.28186 -1409.15435 1290.32037
276 277 278 279 280
-13650.51311 12621.33607 -5548.01648 1675.38507 22594.06659
281 282 283 284 285
5494.01572 11611.06491 23764.81162 13354.70250 3894.04267
286 287 288 289
-9988.33311 8492.67806 -24063.40623 15885.26075
> postscript(file="/var/wessaorg/rcomp/tmp/6bhzv1355344485.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 = 289
Frequency = 1
lag(myerror, k = 1) myerror
0 20054.92650 NA
1 -39267.37959 20054.92650
2 9786.07784 -39267.37959
3 -68655.39762 9786.07784
4 8843.53987 -68655.39762
5 -22137.28872 8843.53987
6 127843.98692 -22137.28872
7 -10188.29732 127843.98692
8 13756.38536 -10188.29732
9 -4582.90784 13756.38536
10 23758.78903 -4582.90784
11 28412.35229 23758.78903
12 -50072.40810 28412.35229
13 -1539.85864 -50072.40810
14 45077.60109 -1539.85864
15 18281.13180 45077.60109
16 -27272.33110 18281.13180
17 17045.97188 -27272.33110
18 8855.12632 17045.97188
19 -33802.37833 8855.12632
20 -15733.65770 -33802.37833
21 -15410.48865 -15733.65770
22 155557.88519 -15410.48865
23 26681.36411 155557.88519
24 -46872.01729 26681.36411
25 -79932.23768 -46872.01729
26 -43003.01463 -79932.23768
27 2834.79234 -43003.01463
28 1774.39902 2834.79234
29 10359.74197 1774.39902
30 32596.60991 10359.74197
31 -20793.22019 32596.60991
32 12505.46924 -20793.22019
33 10629.31176 12505.46924
34 -18846.97710 10629.31176
35 49335.80770 -18846.97710
36 19249.21558 49335.80770
37 70547.51005 19249.21558
38 93944.95077 70547.51005
39 -8526.16917 93944.95077
40 5527.39732 -8526.16917
41 17320.23035 5527.39732
42 -29839.61265 17320.23035
43 -10393.64593 -29839.61265
44 57775.35474 -10393.64593
45 -9856.94210 57775.35474
46 -778.84802 -9856.94210
47 -14625.18748 -778.84802
48 -3847.88783 -14625.18748
49 -26013.24960 -3847.88783
50 140509.74081 -26013.24960
51 -17893.97584 140509.74081
52 -23502.58651 -17893.97584
53 26993.79371 -23502.58651
54 -28429.27699 26993.79371
55 -16345.08969 -28429.27699
56 -11457.85825 -16345.08969
57 -13124.67022 -11457.85825
58 2415.69969 -13124.67022
59 -900.68655 2415.69969
60 45416.38129 -900.68655
61 -2176.44131 45416.38129
62 74138.47965 -2176.44131
63 15532.46808 74138.47965
64 -38028.18836 15532.46808
65 -26987.29162 -38028.18836
66 25273.92420 -26987.29162
67 27651.92408 25273.92420
68 -30658.93521 27651.92408
69 6569.48965 -30658.93521
70 28269.58468 6569.48965
71 28797.29729 28269.58468
72 10779.11246 28797.29729
73 -57387.34914 10779.11246
74 -49053.31498 -57387.34914
75 -20014.17536 -49053.31498
76 -75104.51786 -20014.17536
77 2807.65547 -75104.51786
78 -5172.47956 2807.65547
79 12458.12129 -5172.47956
80 3069.88846 12458.12129
81 -49648.22246 3069.88846
82 49847.30173 -49648.22246
83 36691.62088 49847.30173
84 -30886.05635 36691.62088
85 -44347.75590 -30886.05635
86 -6856.27117 -44347.75590
87 -1483.79093 -6856.27117
88 91083.16917 -1483.79093
89 22525.05152 91083.16917
90 -27767.63674 22525.05152
91 -5973.59560 -27767.63674
92 30223.63795 -5973.59560
93 -49010.98847 30223.63795
94 5995.42331 -49010.98847
95 60141.51087 5995.42331
96 16842.72809 60141.51087
97 -135598.22114 16842.72809
98 36860.54496 -135598.22114
99 -47029.72501 36860.54496
100 25225.44827 -47029.72501
101 29861.25678 25225.44827
102 52229.94578 29861.25678
103 -1496.10635 52229.94578
104 -7431.58973 -1496.10635
105 -27578.34681 -7431.58973
106 -31022.97585 -27578.34681
107 -9537.96730 -31022.97585
108 28512.50807 -9537.96730
109 -28451.70469 28512.50807
110 3287.99066 -28451.70469
111 -51418.05450 3287.99066
112 -18952.06033 -51418.05450
113 -372.56244 -18952.06033
114 -83325.77429 -372.56244
115 -10181.43502 -83325.77429
116 -23347.00057 -10181.43502
117 77396.17203 -23347.00057
118 -46699.03837 77396.17203
119 -19869.02250 -46699.03837
120 -7273.00787 -19869.02250
121 -43326.71348 -7273.00787
122 981.73627 -43326.71348
123 -83732.48715 981.73627
124 58418.77084 -83732.48715
125 14693.53160 58418.77084
126 8681.75671 14693.53160
127 -39054.08098 8681.75671
128 -8790.57889 -39054.08098
129 2480.76056 -8790.57889
130 393.62456 2480.76056
131 -60004.51985 393.62456
132 950.32838 -60004.51985
133 -38014.93214 950.32838
134 -17026.97673 -38014.93214
135 79824.01143 -17026.97673
136 3450.44625 79824.01143
137 69430.05665 3450.44625
138 -29907.57182 69430.05665
139 -16094.02315 -29907.57182
140 -21544.30255 -16094.02315
141 -16437.32710 -21544.30255
142 31613.86609 -16437.32710
143 15457.85536 31613.86609
144 3797.88625 15457.85536
145 -7894.20940 3797.88625
146 -10693.41551 -7894.20940
147 -11202.36786 -10693.41551
148 -15309.30010 -11202.36786
149 13320.18080 -15309.30010
150 38154.02052 13320.18080
151 685.41720 38154.02052
152 -3759.88493 685.41720
153 65614.14323 -3759.88493
154 2892.95766 65614.14323
155 28848.67729 2892.95766
156 14317.45410 28848.67729
157 37853.32040 14317.45410
158 -19143.97638 37853.32040
159 -81762.28252 -19143.97638
160 74574.96492 -81762.28252
161 -2334.05559 74574.96492
162 -43144.93606 -2334.05559
163 40737.61938 -43144.93606
164 -18566.97229 40737.61938
165 8573.59275 -18566.97229
166 -3639.21181 8573.59275
167 -30.14844 -3639.21181
168 -16983.59576 -30.14844
169 105895.45902 -16983.59576
170 84907.24986 105895.45902
171 -29855.16574 84907.24986
172 57572.58778 -29855.16574
173 -5778.68567 57572.58778
174 -125853.79245 -5778.68567
175 55277.15603 -125853.79245
176 -19886.48620 55277.15603
177 14667.88417 -19886.48620
178 -46683.56981 14667.88417
179 43757.65126 -46683.56981
180 -19576.99277 43757.65126
181 -44168.96173 -19576.99277
182 -1529.10063 -44168.96173
183 24406.53971 -1529.10063
184 31784.39311 24406.53971
185 -4914.33595 31784.39311
186 10484.98267 -4914.33595
187 6343.09546 10484.98267
188 36622.51786 6343.09546
189 -7506.50900 36622.51786
190 16914.08240 -7506.50900
191 19039.09175 16914.08240
192 -13595.31643 19039.09175
193 -15124.62639 -13595.31643
194 2576.73039 -15124.62639
195 10468.63810 2576.73039
196 -3569.40924 10468.63810
197 2358.93874 -3569.40924
198 -3634.30126 2358.93874
199 -15658.28958 -3634.30126
200 -22711.67328 -15658.28958
201 -13105.37921 -22711.67328
202 19366.28898 -13105.37921
203 -1244.10864 19366.28898
204 3964.73629 -1244.10864
205 -14171.29845 3964.73629
206 -17528.39711 -14171.29845
207 74133.93679 -17528.39711
208 30181.71646 74133.93679
209 18598.27738 30181.71646
210 17959.20140 18598.27738
211 29363.70716 17959.20140
212 -9723.09919 29363.70716
213 -17340.36256 -9723.09919
214 -7082.32644 -17340.36256
215 -10989.39647 -7082.32644
216 -20652.05036 -10989.39647
217 -4052.56661 -20652.05036
218 -4085.20378 -4052.56661
219 29410.70205 -4085.20378
220 -12468.31457 29410.70205
221 -15592.00306 -12468.31457
222 -22590.18346 -15592.00306
223 -7685.98041 -22590.18346
224 -8860.67173 -7685.98041
225 16248.12282 -8860.67173
226 10104.16354 16248.12282
227 -33643.24839 10104.16354
228 -26906.62737 -33643.24839
229 -22653.20035 -26906.62737
230 -11348.68093 -22653.20035
231 -12256.56991 -11348.68093
232 -11972.36757 -12256.56991
233 -10051.15966 -11972.36757
234 -39217.20067 -10051.15966
235 -47034.33114 -39217.20067
236 75000.57199 -47034.33114
237 -14497.79299 75000.57199
238 -2950.46343 -14497.79299
239 -25678.31287 -2950.46343
240 16904.55674 -25678.31287
241 -52921.73748 16904.55674
242 -42373.71045 -52921.73748
243 -20055.50577 -42373.71045
244 -3480.16592 -20055.50577
245 1871.76317 -3480.16592
246 -18094.35535 1871.76317
247 16595.33607 -18094.35535
248 12979.44426 16595.33607
249 -23463.55757 12979.44426
250 23281.73988 -23463.55757
251 -4780.77125 23281.73988
252 -3750.60904 -4780.77125
253 5497.11791 -3750.60904
254 3367.81366 5497.11791
255 2356.26284 3367.81366
256 -22141.09574 2356.26284
257 3014.53448 -22141.09574
258 -8087.44399 3014.53448
259 -14417.03685 -8087.44399
260 -12273.97483 -14417.03685
261 -14493.11392 -12273.97483
262 6055.51487 -14493.11392
263 -21713.16311 6055.51487
264 -21563.12080 -21713.16311
265 4951.03102 -21563.12080
266 -32834.22855 4951.03102
267 -22102.84879 -32834.22855
268 34128.13524 -22102.84879
269 -33041.99109 34128.13524
270 2884.23135 -33041.99109
271 155.73900 2884.23135
272 20221.28186 155.73900
273 -1409.15435 20221.28186
274 1290.32037 -1409.15435
275 -13650.51311 1290.32037
276 12621.33607 -13650.51311
277 -5548.01648 12621.33607
278 1675.38507 -5548.01648
279 22594.06659 1675.38507
280 5494.01572 22594.06659
281 11611.06491 5494.01572
282 23764.81162 11611.06491
283 13354.70250 23764.81162
284 3894.04267 13354.70250
285 -9988.33311 3894.04267
286 8492.67806 -9988.33311
287 -24063.40623 8492.67806
288 15885.26075 -24063.40623
289 NA 15885.26075
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -39267.37959 20054.92650
[2,] 9786.07784 -39267.37959
[3,] -68655.39762 9786.07784
[4,] 8843.53987 -68655.39762
[5,] -22137.28872 8843.53987
[6,] 127843.98692 -22137.28872
[7,] -10188.29732 127843.98692
[8,] 13756.38536 -10188.29732
[9,] -4582.90784 13756.38536
[10,] 23758.78903 -4582.90784
[11,] 28412.35229 23758.78903
[12,] -50072.40810 28412.35229
[13,] -1539.85864 -50072.40810
[14,] 45077.60109 -1539.85864
[15,] 18281.13180 45077.60109
[16,] -27272.33110 18281.13180
[17,] 17045.97188 -27272.33110
[18,] 8855.12632 17045.97188
[19,] -33802.37833 8855.12632
[20,] -15733.65770 -33802.37833
[21,] -15410.48865 -15733.65770
[22,] 155557.88519 -15410.48865
[23,] 26681.36411 155557.88519
[24,] -46872.01729 26681.36411
[25,] -79932.23768 -46872.01729
[26,] -43003.01463 -79932.23768
[27,] 2834.79234 -43003.01463
[28,] 1774.39902 2834.79234
[29,] 10359.74197 1774.39902
[30,] 32596.60991 10359.74197
[31,] -20793.22019 32596.60991
[32,] 12505.46924 -20793.22019
[33,] 10629.31176 12505.46924
[34,] -18846.97710 10629.31176
[35,] 49335.80770 -18846.97710
[36,] 19249.21558 49335.80770
[37,] 70547.51005 19249.21558
[38,] 93944.95077 70547.51005
[39,] -8526.16917 93944.95077
[40,] 5527.39732 -8526.16917
[41,] 17320.23035 5527.39732
[42,] -29839.61265 17320.23035
[43,] -10393.64593 -29839.61265
[44,] 57775.35474 -10393.64593
[45,] -9856.94210 57775.35474
[46,] -778.84802 -9856.94210
[47,] -14625.18748 -778.84802
[48,] -3847.88783 -14625.18748
[49,] -26013.24960 -3847.88783
[50,] 140509.74081 -26013.24960
[51,] -17893.97584 140509.74081
[52,] -23502.58651 -17893.97584
[53,] 26993.79371 -23502.58651
[54,] -28429.27699 26993.79371
[55,] -16345.08969 -28429.27699
[56,] -11457.85825 -16345.08969
[57,] -13124.67022 -11457.85825
[58,] 2415.69969 -13124.67022
[59,] -900.68655 2415.69969
[60,] 45416.38129 -900.68655
[61,] -2176.44131 45416.38129
[62,] 74138.47965 -2176.44131
[63,] 15532.46808 74138.47965
[64,] -38028.18836 15532.46808
[65,] -26987.29162 -38028.18836
[66,] 25273.92420 -26987.29162
[67,] 27651.92408 25273.92420
[68,] -30658.93521 27651.92408
[69,] 6569.48965 -30658.93521
[70,] 28269.58468 6569.48965
[71,] 28797.29729 28269.58468
[72,] 10779.11246 28797.29729
[73,] -57387.34914 10779.11246
[74,] -49053.31498 -57387.34914
[75,] -20014.17536 -49053.31498
[76,] -75104.51786 -20014.17536
[77,] 2807.65547 -75104.51786
[78,] -5172.47956 2807.65547
[79,] 12458.12129 -5172.47956
[80,] 3069.88846 12458.12129
[81,] -49648.22246 3069.88846
[82,] 49847.30173 -49648.22246
[83,] 36691.62088 49847.30173
[84,] -30886.05635 36691.62088
[85,] -44347.75590 -30886.05635
[86,] -6856.27117 -44347.75590
[87,] -1483.79093 -6856.27117
[88,] 91083.16917 -1483.79093
[89,] 22525.05152 91083.16917
[90,] -27767.63674 22525.05152
[91,] -5973.59560 -27767.63674
[92,] 30223.63795 -5973.59560
[93,] -49010.98847 30223.63795
[94,] 5995.42331 -49010.98847
[95,] 60141.51087 5995.42331
[96,] 16842.72809 60141.51087
[97,] -135598.22114 16842.72809
[98,] 36860.54496 -135598.22114
[99,] -47029.72501 36860.54496
[100,] 25225.44827 -47029.72501
[101,] 29861.25678 25225.44827
[102,] 52229.94578 29861.25678
[103,] -1496.10635 52229.94578
[104,] -7431.58973 -1496.10635
[105,] -27578.34681 -7431.58973
[106,] -31022.97585 -27578.34681
[107,] -9537.96730 -31022.97585
[108,] 28512.50807 -9537.96730
[109,] -28451.70469 28512.50807
[110,] 3287.99066 -28451.70469
[111,] -51418.05450 3287.99066
[112,] -18952.06033 -51418.05450
[113,] -372.56244 -18952.06033
[114,] -83325.77429 -372.56244
[115,] -10181.43502 -83325.77429
[116,] -23347.00057 -10181.43502
[117,] 77396.17203 -23347.00057
[118,] -46699.03837 77396.17203
[119,] -19869.02250 -46699.03837
[120,] -7273.00787 -19869.02250
[121,] -43326.71348 -7273.00787
[122,] 981.73627 -43326.71348
[123,] -83732.48715 981.73627
[124,] 58418.77084 -83732.48715
[125,] 14693.53160 58418.77084
[126,] 8681.75671 14693.53160
[127,] -39054.08098 8681.75671
[128,] -8790.57889 -39054.08098
[129,] 2480.76056 -8790.57889
[130,] 393.62456 2480.76056
[131,] -60004.51985 393.62456
[132,] 950.32838 -60004.51985
[133,] -38014.93214 950.32838
[134,] -17026.97673 -38014.93214
[135,] 79824.01143 -17026.97673
[136,] 3450.44625 79824.01143
[137,] 69430.05665 3450.44625
[138,] -29907.57182 69430.05665
[139,] -16094.02315 -29907.57182
[140,] -21544.30255 -16094.02315
[141,] -16437.32710 -21544.30255
[142,] 31613.86609 -16437.32710
[143,] 15457.85536 31613.86609
[144,] 3797.88625 15457.85536
[145,] -7894.20940 3797.88625
[146,] -10693.41551 -7894.20940
[147,] -11202.36786 -10693.41551
[148,] -15309.30010 -11202.36786
[149,] 13320.18080 -15309.30010
[150,] 38154.02052 13320.18080
[151,] 685.41720 38154.02052
[152,] -3759.88493 685.41720
[153,] 65614.14323 -3759.88493
[154,] 2892.95766 65614.14323
[155,] 28848.67729 2892.95766
[156,] 14317.45410 28848.67729
[157,] 37853.32040 14317.45410
[158,] -19143.97638 37853.32040
[159,] -81762.28252 -19143.97638
[160,] 74574.96492 -81762.28252
[161,] -2334.05559 74574.96492
[162,] -43144.93606 -2334.05559
[163,] 40737.61938 -43144.93606
[164,] -18566.97229 40737.61938
[165,] 8573.59275 -18566.97229
[166,] -3639.21181 8573.59275
[167,] -30.14844 -3639.21181
[168,] -16983.59576 -30.14844
[169,] 105895.45902 -16983.59576
[170,] 84907.24986 105895.45902
[171,] -29855.16574 84907.24986
[172,] 57572.58778 -29855.16574
[173,] -5778.68567 57572.58778
[174,] -125853.79245 -5778.68567
[175,] 55277.15603 -125853.79245
[176,] -19886.48620 55277.15603
[177,] 14667.88417 -19886.48620
[178,] -46683.56981 14667.88417
[179,] 43757.65126 -46683.56981
[180,] -19576.99277 43757.65126
[181,] -44168.96173 -19576.99277
[182,] -1529.10063 -44168.96173
[183,] 24406.53971 -1529.10063
[184,] 31784.39311 24406.53971
[185,] -4914.33595 31784.39311
[186,] 10484.98267 -4914.33595
[187,] 6343.09546 10484.98267
[188,] 36622.51786 6343.09546
[189,] -7506.50900 36622.51786
[190,] 16914.08240 -7506.50900
[191,] 19039.09175 16914.08240
[192,] -13595.31643 19039.09175
[193,] -15124.62639 -13595.31643
[194,] 2576.73039 -15124.62639
[195,] 10468.63810 2576.73039
[196,] -3569.40924 10468.63810
[197,] 2358.93874 -3569.40924
[198,] -3634.30126 2358.93874
[199,] -15658.28958 -3634.30126
[200,] -22711.67328 -15658.28958
[201,] -13105.37921 -22711.67328
[202,] 19366.28898 -13105.37921
[203,] -1244.10864 19366.28898
[204,] 3964.73629 -1244.10864
[205,] -14171.29845 3964.73629
[206,] -17528.39711 -14171.29845
[207,] 74133.93679 -17528.39711
[208,] 30181.71646 74133.93679
[209,] 18598.27738 30181.71646
[210,] 17959.20140 18598.27738
[211,] 29363.70716 17959.20140
[212,] -9723.09919 29363.70716
[213,] -17340.36256 -9723.09919
[214,] -7082.32644 -17340.36256
[215,] -10989.39647 -7082.32644
[216,] -20652.05036 -10989.39647
[217,] -4052.56661 -20652.05036
[218,] -4085.20378 -4052.56661
[219,] 29410.70205 -4085.20378
[220,] -12468.31457 29410.70205
[221,] -15592.00306 -12468.31457
[222,] -22590.18346 -15592.00306
[223,] -7685.98041 -22590.18346
[224,] -8860.67173 -7685.98041
[225,] 16248.12282 -8860.67173
[226,] 10104.16354 16248.12282
[227,] -33643.24839 10104.16354
[228,] -26906.62737 -33643.24839
[229,] -22653.20035 -26906.62737
[230,] -11348.68093 -22653.20035
[231,] -12256.56991 -11348.68093
[232,] -11972.36757 -12256.56991
[233,] -10051.15966 -11972.36757
[234,] -39217.20067 -10051.15966
[235,] -47034.33114 -39217.20067
[236,] 75000.57199 -47034.33114
[237,] -14497.79299 75000.57199
[238,] -2950.46343 -14497.79299
[239,] -25678.31287 -2950.46343
[240,] 16904.55674 -25678.31287
[241,] -52921.73748 16904.55674
[242,] -42373.71045 -52921.73748
[243,] -20055.50577 -42373.71045
[244,] -3480.16592 -20055.50577
[245,] 1871.76317 -3480.16592
[246,] -18094.35535 1871.76317
[247,] 16595.33607 -18094.35535
[248,] 12979.44426 16595.33607
[249,] -23463.55757 12979.44426
[250,] 23281.73988 -23463.55757
[251,] -4780.77125 23281.73988
[252,] -3750.60904 -4780.77125
[253,] 5497.11791 -3750.60904
[254,] 3367.81366 5497.11791
[255,] 2356.26284 3367.81366
[256,] -22141.09574 2356.26284
[257,] 3014.53448 -22141.09574
[258,] -8087.44399 3014.53448
[259,] -14417.03685 -8087.44399
[260,] -12273.97483 -14417.03685
[261,] -14493.11392 -12273.97483
[262,] 6055.51487 -14493.11392
[263,] -21713.16311 6055.51487
[264,] -21563.12080 -21713.16311
[265,] 4951.03102 -21563.12080
[266,] -32834.22855 4951.03102
[267,] -22102.84879 -32834.22855
[268,] 34128.13524 -22102.84879
[269,] -33041.99109 34128.13524
[270,] 2884.23135 -33041.99109
[271,] 155.73900 2884.23135
[272,] 20221.28186 155.73900
[273,] -1409.15435 20221.28186
[274,] 1290.32037 -1409.15435
[275,] -13650.51311 1290.32037
[276,] 12621.33607 -13650.51311
[277,] -5548.01648 12621.33607
[278,] 1675.38507 -5548.01648
[279,] 22594.06659 1675.38507
[280,] 5494.01572 22594.06659
[281,] 11611.06491 5494.01572
[282,] 23764.81162 11611.06491
[283,] 13354.70250 23764.81162
[284,] 3894.04267 13354.70250
[285,] -9988.33311 3894.04267
[286,] 8492.67806 -9988.33311
[287,] -24063.40623 8492.67806
[288,] 15885.26075 -24063.40623
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -39267.37959 20054.92650
2 9786.07784 -39267.37959
3 -68655.39762 9786.07784
4 8843.53987 -68655.39762
5 -22137.28872 8843.53987
6 127843.98692 -22137.28872
7 -10188.29732 127843.98692
8 13756.38536 -10188.29732
9 -4582.90784 13756.38536
10 23758.78903 -4582.90784
11 28412.35229 23758.78903
12 -50072.40810 28412.35229
13 -1539.85864 -50072.40810
14 45077.60109 -1539.85864
15 18281.13180 45077.60109
16 -27272.33110 18281.13180
17 17045.97188 -27272.33110
18 8855.12632 17045.97188
19 -33802.37833 8855.12632
20 -15733.65770 -33802.37833
21 -15410.48865 -15733.65770
22 155557.88519 -15410.48865
23 26681.36411 155557.88519
24 -46872.01729 26681.36411
25 -79932.23768 -46872.01729
26 -43003.01463 -79932.23768
27 2834.79234 -43003.01463
28 1774.39902 2834.79234
29 10359.74197 1774.39902
30 32596.60991 10359.74197
31 -20793.22019 32596.60991
32 12505.46924 -20793.22019
33 10629.31176 12505.46924
34 -18846.97710 10629.31176
35 49335.80770 -18846.97710
36 19249.21558 49335.80770
37 70547.51005 19249.21558
38 93944.95077 70547.51005
39 -8526.16917 93944.95077
40 5527.39732 -8526.16917
41 17320.23035 5527.39732
42 -29839.61265 17320.23035
43 -10393.64593 -29839.61265
44 57775.35474 -10393.64593
45 -9856.94210 57775.35474
46 -778.84802 -9856.94210
47 -14625.18748 -778.84802
48 -3847.88783 -14625.18748
49 -26013.24960 -3847.88783
50 140509.74081 -26013.24960
51 -17893.97584 140509.74081
52 -23502.58651 -17893.97584
53 26993.79371 -23502.58651
54 -28429.27699 26993.79371
55 -16345.08969 -28429.27699
56 -11457.85825 -16345.08969
57 -13124.67022 -11457.85825
58 2415.69969 -13124.67022
59 -900.68655 2415.69969
60 45416.38129 -900.68655
61 -2176.44131 45416.38129
62 74138.47965 -2176.44131
63 15532.46808 74138.47965
64 -38028.18836 15532.46808
65 -26987.29162 -38028.18836
66 25273.92420 -26987.29162
67 27651.92408 25273.92420
68 -30658.93521 27651.92408
69 6569.48965 -30658.93521
70 28269.58468 6569.48965
71 28797.29729 28269.58468
72 10779.11246 28797.29729
73 -57387.34914 10779.11246
74 -49053.31498 -57387.34914
75 -20014.17536 -49053.31498
76 -75104.51786 -20014.17536
77 2807.65547 -75104.51786
78 -5172.47956 2807.65547
79 12458.12129 -5172.47956
80 3069.88846 12458.12129
81 -49648.22246 3069.88846
82 49847.30173 -49648.22246
83 36691.62088 49847.30173
84 -30886.05635 36691.62088
85 -44347.75590 -30886.05635
86 -6856.27117 -44347.75590
87 -1483.79093 -6856.27117
88 91083.16917 -1483.79093
89 22525.05152 91083.16917
90 -27767.63674 22525.05152
91 -5973.59560 -27767.63674
92 30223.63795 -5973.59560
93 -49010.98847 30223.63795
94 5995.42331 -49010.98847
95 60141.51087 5995.42331
96 16842.72809 60141.51087
97 -135598.22114 16842.72809
98 36860.54496 -135598.22114
99 -47029.72501 36860.54496
100 25225.44827 -47029.72501
101 29861.25678 25225.44827
102 52229.94578 29861.25678
103 -1496.10635 52229.94578
104 -7431.58973 -1496.10635
105 -27578.34681 -7431.58973
106 -31022.97585 -27578.34681
107 -9537.96730 -31022.97585
108 28512.50807 -9537.96730
109 -28451.70469 28512.50807
110 3287.99066 -28451.70469
111 -51418.05450 3287.99066
112 -18952.06033 -51418.05450
113 -372.56244 -18952.06033
114 -83325.77429 -372.56244
115 -10181.43502 -83325.77429
116 -23347.00057 -10181.43502
117 77396.17203 -23347.00057
118 -46699.03837 77396.17203
119 -19869.02250 -46699.03837
120 -7273.00787 -19869.02250
121 -43326.71348 -7273.00787
122 981.73627 -43326.71348
123 -83732.48715 981.73627
124 58418.77084 -83732.48715
125 14693.53160 58418.77084
126 8681.75671 14693.53160
127 -39054.08098 8681.75671
128 -8790.57889 -39054.08098
129 2480.76056 -8790.57889
130 393.62456 2480.76056
131 -60004.51985 393.62456
132 950.32838 -60004.51985
133 -38014.93214 950.32838
134 -17026.97673 -38014.93214
135 79824.01143 -17026.97673
136 3450.44625 79824.01143
137 69430.05665 3450.44625
138 -29907.57182 69430.05665
139 -16094.02315 -29907.57182
140 -21544.30255 -16094.02315
141 -16437.32710 -21544.30255
142 31613.86609 -16437.32710
143 15457.85536 31613.86609
144 3797.88625 15457.85536
145 -7894.20940 3797.88625
146 -10693.41551 -7894.20940
147 -11202.36786 -10693.41551
148 -15309.30010 -11202.36786
149 13320.18080 -15309.30010
150 38154.02052 13320.18080
151 685.41720 38154.02052
152 -3759.88493 685.41720
153 65614.14323 -3759.88493
154 2892.95766 65614.14323
155 28848.67729 2892.95766
156 14317.45410 28848.67729
157 37853.32040 14317.45410
158 -19143.97638 37853.32040
159 -81762.28252 -19143.97638
160 74574.96492 -81762.28252
161 -2334.05559 74574.96492
162 -43144.93606 -2334.05559
163 40737.61938 -43144.93606
164 -18566.97229 40737.61938
165 8573.59275 -18566.97229
166 -3639.21181 8573.59275
167 -30.14844 -3639.21181
168 -16983.59576 -30.14844
169 105895.45902 -16983.59576
170 84907.24986 105895.45902
171 -29855.16574 84907.24986
172 57572.58778 -29855.16574
173 -5778.68567 57572.58778
174 -125853.79245 -5778.68567
175 55277.15603 -125853.79245
176 -19886.48620 55277.15603
177 14667.88417 -19886.48620
178 -46683.56981 14667.88417
179 43757.65126 -46683.56981
180 -19576.99277 43757.65126
181 -44168.96173 -19576.99277
182 -1529.10063 -44168.96173
183 24406.53971 -1529.10063
184 31784.39311 24406.53971
185 -4914.33595 31784.39311
186 10484.98267 -4914.33595
187 6343.09546 10484.98267
188 36622.51786 6343.09546
189 -7506.50900 36622.51786
190 16914.08240 -7506.50900
191 19039.09175 16914.08240
192 -13595.31643 19039.09175
193 -15124.62639 -13595.31643
194 2576.73039 -15124.62639
195 10468.63810 2576.73039
196 -3569.40924 10468.63810
197 2358.93874 -3569.40924
198 -3634.30126 2358.93874
199 -15658.28958 -3634.30126
200 -22711.67328 -15658.28958
201 -13105.37921 -22711.67328
202 19366.28898 -13105.37921
203 -1244.10864 19366.28898
204 3964.73629 -1244.10864
205 -14171.29845 3964.73629
206 -17528.39711 -14171.29845
207 74133.93679 -17528.39711
208 30181.71646 74133.93679
209 18598.27738 30181.71646
210 17959.20140 18598.27738
211 29363.70716 17959.20140
212 -9723.09919 29363.70716
213 -17340.36256 -9723.09919
214 -7082.32644 -17340.36256
215 -10989.39647 -7082.32644
216 -20652.05036 -10989.39647
217 -4052.56661 -20652.05036
218 -4085.20378 -4052.56661
219 29410.70205 -4085.20378
220 -12468.31457 29410.70205
221 -15592.00306 -12468.31457
222 -22590.18346 -15592.00306
223 -7685.98041 -22590.18346
224 -8860.67173 -7685.98041
225 16248.12282 -8860.67173
226 10104.16354 16248.12282
227 -33643.24839 10104.16354
228 -26906.62737 -33643.24839
229 -22653.20035 -26906.62737
230 -11348.68093 -22653.20035
231 -12256.56991 -11348.68093
232 -11972.36757 -12256.56991
233 -10051.15966 -11972.36757
234 -39217.20067 -10051.15966
235 -47034.33114 -39217.20067
236 75000.57199 -47034.33114
237 -14497.79299 75000.57199
238 -2950.46343 -14497.79299
239 -25678.31287 -2950.46343
240 16904.55674 -25678.31287
241 -52921.73748 16904.55674
242 -42373.71045 -52921.73748
243 -20055.50577 -42373.71045
244 -3480.16592 -20055.50577
245 1871.76317 -3480.16592
246 -18094.35535 1871.76317
247 16595.33607 -18094.35535
248 12979.44426 16595.33607
249 -23463.55757 12979.44426
250 23281.73988 -23463.55757
251 -4780.77125 23281.73988
252 -3750.60904 -4780.77125
253 5497.11791 -3750.60904
254 3367.81366 5497.11791
255 2356.26284 3367.81366
256 -22141.09574 2356.26284
257 3014.53448 -22141.09574
258 -8087.44399 3014.53448
259 -14417.03685 -8087.44399
260 -12273.97483 -14417.03685
261 -14493.11392 -12273.97483
262 6055.51487 -14493.11392
263 -21713.16311 6055.51487
264 -21563.12080 -21713.16311
265 4951.03102 -21563.12080
266 -32834.22855 4951.03102
267 -22102.84879 -32834.22855
268 34128.13524 -22102.84879
269 -33041.99109 34128.13524
270 2884.23135 -33041.99109
271 155.73900 2884.23135
272 20221.28186 155.73900
273 -1409.15435 20221.28186
274 1290.32037 -1409.15435
275 -13650.51311 1290.32037
276 12621.33607 -13650.51311
277 -5548.01648 12621.33607
278 1675.38507 -5548.01648
279 22594.06659 1675.38507
280 5494.01572 22594.06659
281 11611.06491 5494.01572
282 23764.81162 11611.06491
283 13354.70250 23764.81162
284 3894.04267 13354.70250
285 -9988.33311 3894.04267
286 8492.67806 -9988.33311
287 -24063.40623 8492.67806
288 15885.26075 -24063.40623
> 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/wessaorg/rcomp/tmp/7jbbg1355344485.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/wessaorg/rcomp/tmp/83zqy1355344485.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/wessaorg/rcomp/tmp/9nv121355344485.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/wessaorg/rcomp/tmp/10q1571355344485.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/11k7oq1355344485.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/wessaorg/rcomp/tmp/127dg71355344485.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/wessaorg/rcomp/tmp/13qj0k1355344485.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/wessaorg/rcomp/tmp/143nhd1355344485.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/wessaorg/rcomp/tmp/15h9e81355344485.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/wessaorg/rcomp/tmp/16bv941355344485.tab")
+ }
>
> try(system("convert tmp/1jgzx1355344485.ps tmp/1jgzx1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/2go0p1355344485.ps tmp/2go0p1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/3c60h1355344485.ps tmp/3c60h1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/405bj1355344485.ps tmp/405bj1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/50ele1355344485.ps tmp/50ele1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/6bhzv1355344485.ps tmp/6bhzv1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/7jbbg1355344485.ps tmp/7jbbg1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/83zqy1355344485.ps tmp/83zqy1355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/9nv121355344485.ps tmp/9nv121355344485.png",intern=TRUE))
character(0)
> try(system("convert tmp/10q1571355344485.ps tmp/10q1571355344485.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
16.466 1.392 18.059