R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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
+ ,115
+ ,112285
+ ,24188
+ ,146283
+ ,144
+ ,120982
+ ,109
+ ,84786
+ ,18273
+ ,98364
+ ,103
+ ,176508
+ ,146
+ ,83123
+ ,14130
+ ,86146
+ ,98
+ ,179321
+ ,116
+ ,101193
+ ,32287
+ ,96933
+ ,135
+ ,123185
+ ,68
+ ,38361
+ ,8654
+ ,79234
+ ,61
+ ,52746
+ ,101
+ ,68504
+ ,9245
+ ,42551
+ ,39
+ ,385534
+ ,96
+ ,119182
+ ,33251
+ ,195663
+ ,150
+ ,33170
+ ,67
+ ,22807
+ ,1271
+ ,6853
+ ,5
+ ,101645
+ ,44
+ ,17140
+ ,5279
+ ,21529
+ ,28
+ ,149061
+ ,100
+ ,116174
+ ,27101
+ ,95757
+ ,84
+ ,165446
+ ,93
+ ,57635
+ ,16373
+ ,85584
+ ,80
+ ,237213
+ ,140
+ ,66198
+ ,19716
+ ,143983
+ ,130
+ ,173326
+ ,166
+ ,71701
+ ,17753
+ ,75851
+ ,82
+ ,133131
+ ,99
+ ,57793
+ ,9028
+ ,59238
+ ,60
+ ,258873
+ ,139
+ ,80444
+ ,18653
+ ,93163
+ ,131
+ ,180083
+ ,130
+ ,53855
+ ,8828
+ ,96037
+ ,84
+ ,324799
+ ,181
+ ,97668
+ ,29498
+ ,151511
+ ,140
+ ,230964
+ ,116
+ ,133824
+ ,27563
+ ,136368
+ ,151
+ ,236785
+ ,116
+ ,101481
+ ,18293
+ ,112642
+ ,91
+ ,135473
+ ,88
+ ,99645
+ ,22530
+ ,94728
+ ,138
+ ,202925
+ ,139
+ ,114789
+ ,15977
+ ,105499
+ ,150
+ ,215147
+ ,135
+ ,99052
+ ,35082
+ ,121527
+ ,124
+ ,344297
+ ,108
+ ,67654
+ ,16116
+ ,127766
+ ,119
+ ,153935
+ ,89
+ ,65553
+ ,15849
+ ,98958
+ ,73
+ ,132943
+ ,156
+ ,97500
+ ,16026
+ ,77900
+ ,110
+ ,174724
+ ,129
+ ,69112
+ ,26569
+ ,85646
+ ,123
+ ,174415
+ ,118
+ ,82753
+ ,24785
+ ,98579
+ ,90
+ ,225548
+ ,118
+ ,85323
+ ,17569
+ ,130767
+ ,116
+ ,223632
+ ,125
+ ,72654
+ ,23825
+ ,131741
+ ,113
+ ,124817
+ ,95
+ ,30727
+ ,7869
+ ,53907
+ ,56
+ ,221698
+ ,126
+ ,77873
+ ,14975
+ ,178812
+ ,115
+ ,210767
+ ,135
+ ,117478
+ ,37791
+ ,146761
+ ,119
+ ,170266
+ ,154
+ ,74007
+ ,9605
+ ,82036
+ ,129
+ ,260561
+ ,165
+ ,90183
+ ,27295
+ ,163253
+ ,127
+ ,84853
+ ,113
+ ,61542
+ ,2746
+ ,27032
+ ,27
+ ,294424
+ ,127
+ ,101494
+ ,34461
+ ,171975
+ ,175
+ ,101011
+ ,52
+ ,27570
+ ,8098
+ ,65990
+ ,35
+ ,215641
+ ,121
+ ,55813
+ ,4787
+ ,86572
+ ,64
+ ,325107
+ ,136
+ ,79215
+ ,24919
+ ,159676
+ ,96
+ ,7176
+ ,0
+ ,1423
+ ,603
+ ,1929
+ ,0
+ ,167542
+ ,108
+ ,55461
+ ,16329
+ ,85371
+ ,84
+ ,106408
+ ,46
+ ,31081
+ ,12558
+ ,58391
+ ,41
+ ,96560
+ ,54
+ ,22996
+ ,7784
+ ,31580
+ ,47
+ ,265769
+ ,124
+ ,83122
+ ,28522
+ ,136815
+ ,126
+ ,269651
+ ,115
+ ,70106
+ ,22265
+ ,120642
+ ,105
+ ,149112
+ ,128
+ ,60578
+ ,14459
+ ,69107
+ ,80
+ ,175824
+ ,80
+ ,39992
+ ,14526
+ ,50495
+ ,70
+ ,152871
+ ,97
+ ,79892
+ ,22240
+ ,108016
+ ,73
+ ,111665
+ ,104
+ ,49810
+ ,11802
+ ,46341
+ ,57
+ ,116408
+ ,59
+ ,71570
+ ,7623
+ ,78348
+ ,40
+ ,362301
+ ,125
+ ,100708
+ ,11912
+ ,79336
+ ,68
+ ,78800
+ ,82
+ ,33032
+ ,7935
+ ,56968
+ ,21
+ ,183167
+ ,149
+ ,82875
+ ,18220
+ ,93176
+ ,127
+ ,277965
+ ,149
+ ,139077
+ ,19199
+ ,161632
+ ,154
+ ,150629
+ ,122
+ ,71595
+ ,19918
+ ,87850
+ ,116
+ ,168809
+ ,118
+ ,72260
+ ,21884
+ ,127969
+ ,102
+ ,24188
+ ,12
+ ,5950
+ ,2694
+ ,15049
+ ,7
+ ,329267
+ ,144
+ ,115762
+ ,15808
+ ,155135
+ ,148
+ ,65029
+ ,67
+ ,32551
+ ,3597
+ ,25109
+ ,21
+ ,101097
+ ,52
+ ,31701
+ ,5296
+ ,45824
+ ,35
+ ,218946
+ ,108
+ ,80670
+ ,25239
+ ,102996
+ ,112
+ ,244052
+ ,166
+ ,143558
+ ,29801
+ ,160604
+ ,137
+ ,341570
+ ,80
+ ,117105
+ ,18450
+ ,158051
+ ,135
+ ,103597
+ ,60
+ ,23789
+ ,7132
+ ,44547
+ ,26
+ ,233328
+ ,107
+ ,120733
+ ,34861
+ ,162647
+ ,230
+ ,256462
+ ,127
+ ,105195
+ ,35940
+ ,174141
+ ,181
+ ,206161
+ ,107
+ ,73107
+ ,16688
+ ,60622
+ ,71
+ ,311473
+ ,146
+ ,132068
+ ,24683
+ ,179566
+ ,147
+ ,235800
+ ,84
+ ,149193
+ ,46230
+ ,184301
+ ,190
+ ,177939
+ ,141
+ ,46821
+ ,10387
+ ,75661
+ ,64
+ ,207176
+ ,123
+ ,87011
+ ,21436
+ ,96144
+ ,105
+ ,196553
+ ,111
+ ,95260
+ ,30546
+ ,129847
+ ,107
+ ,174184
+ ,98
+ ,55183
+ ,19746
+ ,117286
+ ,94
+ ,143246
+ ,105
+ ,106671
+ ,15977
+ ,71180
+ ,116
+ ,187559
+ ,135
+ ,73511
+ ,22583
+ ,109377
+ ,106
+ ,187681
+ ,107
+ ,92945
+ ,17274
+ ,85298
+ ,143
+ ,119016
+ ,85
+ ,78664
+ ,16469
+ ,73631
+ ,81
+ ,182192
+ ,155
+ ,70054
+ ,14251
+ ,86767
+ ,89
+ ,73566
+ ,88
+ ,22618
+ ,3007
+ ,23824
+ ,26
+ ,194979
+ ,155
+ ,74011
+ ,16851
+ ,93487
+ ,84
+ ,167488
+ ,104
+ ,83737
+ ,21113
+ ,82981
+ ,113
+ ,143756
+ ,132
+ ,69094
+ ,17401
+ ,73815
+ ,120
+ ,275541
+ ,127
+ ,93133
+ ,23958
+ ,94552
+ ,110
+ ,243199
+ ,108
+ ,95536
+ ,23567
+ ,132190
+ ,134
+ ,182999
+ ,129
+ ,225920
+ ,13065
+ ,128754
+ ,54
+ ,135649
+ ,116
+ ,62133
+ ,15358
+ ,66363
+ ,96
+ ,152299
+ ,122
+ ,61370
+ ,14587
+ ,67808
+ ,78
+ ,120221
+ ,85
+ ,43836
+ ,12770
+ ,61724
+ ,51
+ ,346485
+ ,147
+ ,106117
+ ,24021
+ ,131722
+ ,121
+ ,145790
+ ,99
+ ,38692
+ ,9648
+ ,68580
+ ,38
+ ,193339
+ ,87
+ ,84651
+ ,20537
+ ,106175
+ ,145
+ ,80953
+ ,28
+ ,56622
+ ,7905
+ ,55792
+ ,59
+ ,122774
+ ,90
+ ,15986
+ ,4527
+ ,25157
+ ,27
+ ,130585
+ ,109
+ ,95364
+ ,30495
+ ,76669
+ ,91
+ ,112611
+ ,78
+ ,26706
+ ,7117
+ ,57283
+ ,48
+ ,286468
+ ,111
+ ,89691
+ ,17719
+ ,105805
+ ,68
+ ,241066
+ ,158
+ ,67267
+ ,27056
+ ,129484
+ ,58
+ ,148446
+ ,141
+ ,126846
+ ,33473
+ ,72413
+ ,150
+ ,204713
+ ,122
+ ,41140
+ ,9758
+ ,87831
+ ,74
+ ,182079
+ ,124
+ ,102860
+ ,21115
+ ,96971
+ ,181
+ ,140344
+ ,93
+ ,51715
+ ,7236
+ ,71299
+ ,65
+ ,220516
+ ,124
+ ,55801
+ ,13790
+ ,77494
+ ,97
+ ,243060
+ ,112
+ ,111813
+ ,32902
+ ,120336
+ ,121
+ ,162765
+ ,108
+ ,120293
+ ,25131
+ ,93913
+ ,99
+ ,182613
+ ,99
+ ,138599
+ ,30910
+ ,136048
+ ,152
+ ,232138
+ ,117
+ ,161647
+ ,35947
+ ,181248
+ ,188
+ ,265318
+ ,199
+ ,115929
+ ,29848
+ ,146123
+ ,138
+ ,85574
+ ,78
+ ,24266
+ ,6943
+ ,32036
+ ,40
+ ,310839
+ ,91
+ ,162901
+ ,42705
+ ,186646
+ ,254
+ ,225060
+ ,158
+ ,109825
+ ,31808
+ ,102255
+ ,87
+ ,232317
+ ,126
+ ,129838
+ ,26675
+ ,168237
+ ,178
+ ,144966
+ ,122
+ ,37510
+ ,8435
+ ,64219
+ ,51
+ ,43287
+ ,71
+ ,43750
+ ,7409
+ ,19630
+ ,49
+ ,155754
+ ,75
+ ,40652
+ ,14993
+ ,76825
+ ,73
+ ,164709
+ ,115
+ ,87771
+ ,36867
+ ,115338
+ ,176
+ ,201940
+ ,119
+ ,85872
+ ,33835
+ ,109427
+ ,94
+ ,235454
+ ,124
+ ,89275
+ ,24164
+ ,118168
+ ,120
+ ,220801
+ ,72
+ ,44418
+ ,12607
+ ,84845
+ ,66
+ ,99466
+ ,91
+ ,192565
+ ,22609
+ ,153197
+ ,56
+ ,92661
+ ,45
+ ,35232
+ ,5892
+ ,29877
+ ,39
+ ,133328
+ ,78
+ ,40909
+ ,17014
+ ,63506
+ ,66
+ ,61361
+ ,39
+ ,13294
+ ,5394
+ ,22445
+ ,27
+ ,125930
+ ,68
+ ,32387
+ ,9178
+ ,47695
+ ,65
+ ,100750
+ ,119
+ ,140867
+ ,6440
+ ,68370
+ ,58
+ ,224549
+ ,117
+ ,120662
+ ,21916
+ ,146304
+ ,98
+ ,82316
+ ,39
+ ,21233
+ ,4011
+ ,38233
+ ,25
+ ,102010
+ ,50
+ ,44332
+ ,5818
+ ,42071
+ ,26
+ ,101523
+ ,88
+ ,61056
+ ,18647
+ ,50517
+ ,77
+ ,243511
+ ,155
+ ,101338
+ ,20556
+ ,103950
+ ,130
+ ,22938
+ ,0
+ ,1168
+ ,238
+ ,5841
+ ,11
+ ,41566
+ ,36
+ ,13497
+ ,70
+ ,2341
+ ,2
+ ,152474
+ ,123
+ ,65567
+ ,22392
+ ,84396
+ ,101
+ ,61857
+ ,32
+ ,25162
+ ,3913
+ ,24610
+ ,31
+ ,99923
+ ,99
+ ,32334
+ ,12237
+ ,35753
+ ,36
+ ,132487
+ ,136
+ ,40735
+ ,8388
+ ,55515
+ ,120
+ ,317394
+ ,117
+ ,91413
+ ,22120
+ ,209056
+ ,195
+ ,21054
+ ,0
+ ,855
+ ,338
+ ,6622
+ ,4
+ ,209641
+ ,88
+ ,97068
+ ,11727
+ ,115814
+ ,89
+ ,22648
+ ,39
+ ,44339
+ ,3704
+ ,11609
+ ,24
+ ,31414
+ ,25
+ ,14116
+ ,3988
+ ,13155
+ ,39
+ ,46698
+ ,52
+ ,10288
+ ,3030
+ ,18274
+ ,14
+ ,131698
+ ,75
+ ,65622
+ ,13520
+ ,72875
+ ,78
+ ,91735
+ ,71
+ ,16563
+ ,1421
+ ,10112
+ ,15
+ ,244749
+ ,124
+ ,76643
+ ,20923
+ ,142775
+ ,106
+ ,184510
+ ,151
+ ,110681
+ ,20237
+ ,68847
+ ,83
+ ,79863
+ ,71
+ ,29011
+ ,3219
+ ,17659
+ ,24
+ ,128423
+ ,145
+ ,92696
+ ,3769
+ ,20112
+ ,37
+ ,97839
+ ,87
+ ,94785
+ ,12252
+ ,61023
+ ,77
+ ,38214
+ ,27
+ ,8773
+ ,1888
+ ,13983
+ ,16
+ ,151101
+ ,131
+ ,83209
+ ,14497
+ ,65176
+ ,56
+ ,272458
+ ,162
+ ,93815
+ ,28864
+ ,132432
+ ,132
+ ,172494
+ ,165
+ ,86687
+ ,21721
+ ,112494
+ ,144
+ ,108043
+ ,54
+ ,34553
+ ,4821
+ ,45109
+ ,40
+ ,328107
+ ,159
+ ,105547
+ ,33644
+ ,170875
+ ,153
+ ,250579
+ ,147
+ ,103487
+ ,15923
+ ,180759
+ ,143
+ ,351067
+ ,170
+ ,213688
+ ,42935
+ ,214921
+ ,220
+ ,158015
+ ,119
+ ,71220
+ ,18864
+ ,100226
+ ,79
+ ,98866
+ ,49
+ ,23517
+ ,4977
+ ,32043
+ ,50
+ ,85439
+ ,104
+ ,56926
+ ,7785
+ ,54454
+ ,39
+ ,229242
+ ,120
+ ,91721
+ ,17939
+ ,78876
+ ,95
+ ,351619
+ ,150
+ ,115168
+ ,23436
+ ,170745
+ ,169
+ ,84207
+ ,112
+ ,111194
+ ,325
+ ,6940
+ ,12
+ ,120445
+ ,59
+ ,51009
+ ,13539
+ ,49025
+ ,63
+ ,324598
+ ,136
+ ,135777
+ ,34538
+ ,122037
+ ,134
+ ,131069
+ ,107
+ ,51513
+ ,12198
+ ,53782
+ ,69
+ ,204271
+ ,130
+ ,74163
+ ,26924
+ ,127748
+ ,119
+ ,165543
+ ,115
+ ,51633
+ ,12716
+ ,86839
+ ,119
+ ,141722
+ ,107
+ ,75345
+ ,8172
+ ,44830
+ ,75
+ ,116048
+ ,75
+ ,33416
+ ,10855
+ ,77395
+ ,63
+ ,250047
+ ,71
+ ,83305
+ ,11932
+ ,89324
+ ,55
+ ,299775
+ ,120
+ ,98952
+ ,14300
+ ,103300
+ ,103
+ ,195838
+ ,116
+ ,102372
+ ,25515
+ ,112283
+ ,197
+ ,173260
+ ,79
+ ,37238
+ ,2805
+ ,10901
+ ,16
+ ,254488
+ ,150
+ ,103772
+ ,29402
+ ,120691
+ ,140
+ ,104389
+ ,156
+ ,123969
+ ,16440
+ ,58106
+ ,89
+ ,136084
+ ,51
+ ,27142
+ ,11221
+ ,57140
+ ,40
+ ,199476
+ ,118
+ ,135400
+ ,28732
+ ,122422
+ ,125
+ ,92499
+ ,71
+ ,21399
+ ,5250
+ ,25899
+ ,21
+ ,224330
+ ,144
+ ,130115
+ ,28608
+ ,139296
+ ,167
+ ,135781
+ ,47
+ ,24874
+ ,8092
+ ,52678
+ ,32
+ ,74408
+ ,28
+ ,34988
+ ,4473
+ ,23853
+ ,36
+ ,81240
+ ,68
+ ,45549
+ ,1572
+ ,17306
+ ,13
+ ,14688
+ ,0
+ ,6023
+ ,2065
+ ,7953
+ ,5
+ ,181633
+ ,110
+ ,64466
+ ,14817
+ ,89455
+ ,96
+ ,271856
+ ,147
+ ,54990
+ ,16714
+ ,147866
+ ,151
+ ,7199
+ ,0
+ ,1644
+ ,556
+ ,4245
+ ,6
+ ,46660
+ ,15
+ ,6179
+ ,2089
+ ,21509
+ ,13
+ ,17547
+ ,4
+ ,3926
+ ,2658
+ ,7670
+ ,3
+ ,133368
+ ,64
+ ,32755
+ ,10695
+ ,66675
+ ,57
+ ,95227
+ ,111
+ ,34777
+ ,1669
+ ,14336
+ ,23
+ ,152601
+ ,85
+ ,73224
+ ,16267
+ ,53608
+ ,61
+ ,98146
+ ,68
+ ,27114
+ ,7768
+ ,30059
+ ,21
+ ,79619
+ ,40
+ ,20760
+ ,7252
+ ,29668
+ ,43
+ ,59194
+ ,80
+ ,37636
+ ,6387
+ ,22097
+ ,20
+ ,139942
+ ,88
+ ,65461
+ ,18715
+ ,96841
+ ,82
+ ,118612
+ ,48
+ ,30080
+ ,7936
+ ,41907
+ ,90
+ ,72880
+ ,76
+ ,24094
+ ,8643
+ ,27080
+ ,25
+ ,65475
+ ,51
+ ,69008
+ ,7294
+ ,35885
+ ,60
+ ,99643
+ ,67
+ ,54968
+ ,4570
+ ,41247
+ ,61
+ ,71965
+ ,59
+ ,46090
+ ,7185
+ ,28313
+ ,85
+ ,77272
+ ,61
+ ,27507
+ ,10058
+ ,36845
+ ,43
+ ,49289
+ ,76
+ ,10672
+ ,2342
+ ,16548
+ ,25
+ ,135131
+ ,60
+ ,34029
+ ,8509
+ ,36134
+ ,41
+ ,108446
+ ,68
+ ,46300
+ ,13275
+ ,55764
+ ,26
+ ,89746
+ ,71
+ ,24760
+ ,6816
+ ,28910
+ ,38
+ ,44296
+ ,76
+ ,18779
+ ,1930
+ ,13339
+ ,12
+ ,77648
+ ,62
+ ,21280
+ ,8086
+ ,25319
+ ,29
+ ,181528
+ ,61
+ ,40662
+ ,10737
+ ,66956
+ ,49
+ ,134019
+ ,67
+ ,28987
+ ,8033
+ ,47487
+ ,46
+ ,124064
+ ,88
+ ,22827
+ ,7058
+ ,52785
+ ,41
+ ,92630
+ ,30
+ ,18513
+ ,6782
+ ,44683
+ ,31
+ ,121848
+ ,64
+ ,30594
+ ,5401
+ ,35619
+ ,41
+ ,52915
+ ,68
+ ,24006
+ ,6521
+ ,21920
+ ,26
+ ,81872
+ ,64
+ ,27913
+ ,10856
+ ,45608
+ ,23
+ ,58981
+ ,91
+ ,42744
+ ,2154
+ ,7721
+ ,14
+ ,53515
+ ,88
+ ,12934
+ ,6117
+ ,20634
+ ,16
+ ,60812
+ ,52
+ ,22574
+ ,5238
+ ,29788
+ ,25
+ ,56375
+ ,49
+ ,41385
+ ,4820
+ ,31931
+ ,21
+ ,65490
+ ,62
+ ,18653
+ ,5615
+ ,37754
+ ,32
+ ,80949
+ ,61
+ ,18472
+ ,4272
+ ,32505
+ ,9
+ ,76302
+ ,76
+ ,30976
+ ,8702
+ ,40557
+ ,35
+ ,104011
+ ,88
+ ,63339
+ ,15340
+ ,94238
+ ,42
+ ,98104
+ ,66
+ ,25568
+ ,8030
+ ,44197
+ ,68
+ ,67989
+ ,71
+ ,33747
+ ,9526
+ ,43228
+ ,32
+ ,30989
+ ,68
+ ,4154
+ ,1278
+ ,4103
+ ,6
+ ,135458
+ ,48
+ ,19474
+ ,4236
+ ,44144
+ ,68
+ ,73504
+ ,25
+ ,35130
+ ,3023
+ ,32868
+ ,33
+ ,63123
+ ,68
+ ,39067
+ ,7196
+ ,27640
+ ,84
+ ,61254
+ ,41
+ ,13310
+ ,3394
+ ,14063
+ ,46
+ ,74914
+ ,90
+ ,65892
+ ,6371
+ ,28990
+ ,30
+ ,31774
+ ,66
+ ,4143
+ ,1574
+ ,4694
+ ,0
+ ,81437
+ ,54
+ ,28579
+ ,9620
+ ,42648
+ ,36
+ ,87186
+ ,59
+ ,51776
+ ,6978
+ ,64329
+ ,47
+ ,50090
+ ,60
+ ,21152
+ ,4911
+ ,21928
+ ,20
+ ,65745
+ ,77
+ ,38084
+ ,8645
+ ,25836
+ ,50
+ ,56653
+ ,68
+ ,27717
+ ,8987
+ ,22779
+ ,30
+ ,158399
+ ,72
+ ,32928
+ ,5544
+ ,40820
+ ,30
+ ,46455
+ ,67
+ ,11342
+ ,3083
+ ,27530
+ ,34
+ ,73624
+ ,64
+ ,19499
+ ,6909
+ ,32378
+ ,33
+ ,38395
+ ,63
+ ,16380
+ ,3189
+ ,10824
+ ,34
+ ,91899
+ ,59
+ ,36874
+ ,6745
+ ,39613
+ ,37
+ ,139526
+ ,84
+ ,48259
+ ,16724
+ ,60865
+ ,83
+ ,52164
+ ,64
+ ,16734
+ ,4850
+ ,19787
+ ,32
+ ,51567
+ ,56
+ ,28207
+ ,7025
+ ,20107
+ ,30
+ ,70551
+ ,54
+ ,30143
+ ,6047
+ ,36605
+ ,43
+ ,84856
+ ,67
+ ,41369
+ ,7377
+ ,40961
+ ,41
+ ,102538
+ ,58
+ ,45833
+ ,9078
+ ,48231
+ ,51
+ ,86678
+ ,59
+ ,29156
+ ,4605
+ ,39725
+ ,19
+ ,85709
+ ,40
+ ,35944
+ ,3238
+ ,21455
+ ,37
+ ,34662
+ ,22
+ ,36278
+ ,8100
+ ,23430
+ ,33
+ ,150580
+ ,83
+ ,45588
+ ,9653
+ ,62991
+ ,41
+ ,99611
+ ,81
+ ,45097
+ ,8914
+ ,49363
+ ,54
+ ,19349
+ ,2
+ ,3895
+ ,786
+ ,9604
+ ,14
+ ,99373
+ ,72
+ ,28394
+ ,6700
+ ,24552
+ ,25
+ ,86230
+ ,61
+ ,18632
+ ,5788
+ ,31493
+ ,25
+ ,30837
+ ,15
+ ,2325
+ ,593
+ ,3439
+ ,8
+ ,31706
+ ,32
+ ,25139
+ ,4506
+ ,19555
+ ,26
+ ,89806
+ ,62
+ ,27975
+ ,6382
+ ,21228
+ ,20
+ ,62088
+ ,58
+ ,14483
+ ,5621
+ ,23177
+ ,11
+ ,40151
+ ,36
+ ,13127
+ ,3997
+ ,22094
+ ,14
+ ,27634
+ ,59
+ ,5839
+ ,520
+ ,2342
+ ,3
+ ,76990
+ ,68
+ ,24069
+ ,8891
+ ,38798
+ ,40
+ ,37460
+ ,21
+ ,3738
+ ,999
+ ,3255
+ ,5
+ ,54157
+ ,55
+ ,18625
+ ,7067
+ ,24261
+ ,38
+ ,49862
+ ,54
+ ,36341
+ ,4639
+ ,18511
+ ,32
+ ,84337
+ ,55
+ ,24548
+ ,5654
+ ,40798
+ ,41
+ ,64175
+ ,72
+ ,21792
+ ,6928
+ ,28893
+ ,46
+ ,59382
+ ,41
+ ,26263
+ ,1514
+ ,21425
+ ,47
+ ,119308
+ ,61
+ ,23686
+ ,9238
+ ,50276
+ ,37
+ ,76702
+ ,67
+ ,49303
+ ,8204
+ ,37643
+ ,51
+ ,103425
+ ,76
+ ,25659
+ ,5926
+ ,30377
+ ,49
+ ,70344
+ ,64
+ ,28904
+ ,5785
+ ,27126
+ ,21
+ ,43410
+ ,3
+ ,2781
+ ,4
+ ,13
+ ,1
+ ,104838
+ ,63
+ ,29236
+ ,5930
+ ,42097
+ ,44
+ ,62215
+ ,40
+ ,19546
+ ,3710
+ ,24451
+ ,26
+ ,69304
+ ,69
+ ,22818
+ ,705
+ ,14335
+ ,21
+ ,53117
+ ,48
+ ,32689
+ ,443
+ ,5084
+ ,4
+ ,19764
+ ,8
+ ,5752
+ ,2416
+ ,9927
+ ,10
+ ,86680
+ ,52
+ ,22197
+ ,7747
+ ,43527
+ ,43
+ ,84105
+ ,66
+ ,20055
+ ,5432
+ ,27184
+ ,34
+ ,77945
+ ,76
+ ,25272
+ ,4913
+ ,21610
+ ,32
+ ,89113
+ ,43
+ ,82206
+ ,2650
+ ,20484
+ ,20
+ ,91005
+ ,39
+ ,32073
+ ,2370
+ ,20156
+ ,34
+ ,40248
+ ,14
+ ,5444
+ ,775
+ ,6012
+ ,6
+ ,64187
+ ,61
+ ,20154
+ ,5576
+ ,18475
+ ,12
+ ,50857
+ ,71
+ ,36944
+ ,1352
+ ,12645
+ ,24
+ ,56613
+ ,44
+ ,8019
+ ,3080
+ ,11017
+ ,16
+ ,62792
+ ,60
+ ,30884
+ ,10205
+ ,37623
+ ,72
+ ,72535
+ ,64
+ ,19540
+ ,6095
+ ,35873
+ ,27)
+ ,dim=c(6
+ ,289)
+ ,dimnames=list(c('time_in_rfc'
+ ,'feedback_messages_p1'
+ ,'totsize'
+ ,'totrevisions'
+ ,'totseconds'
+ ,'tothyperlinks
')
+ ,1:289))
> y <- array(NA,dim=c(6,289),dimnames=list(c('time_in_rfc','feedback_messages_p1','totsize','totrevisions','totseconds','tothyperlinks
'),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'
> #'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
> 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 feedback_messages_p1 totsize totrevisions totseconds
1 210907 115 112285 24188 146283
2 120982 109 84786 18273 98364
3 176508 146 83123 14130 86146
4 179321 116 101193 32287 96933
5 123185 68 38361 8654 79234
6 52746 101 68504 9245 42551
7 385534 96 119182 33251 195663
8 33170 67 22807 1271 6853
9 101645 44 17140 5279 21529
10 149061 100 116174 27101 95757
11 165446 93 57635 16373 85584
12 237213 140 66198 19716 143983
13 173326 166 71701 17753 75851
14 133131 99 57793 9028 59238
15 258873 139 80444 18653 93163
16 180083 130 53855 8828 96037
17 324799 181 97668 29498 151511
18 230964 116 133824 27563 136368
19 236785 116 101481 18293 112642
20 135473 88 99645 22530 94728
21 202925 139 114789 15977 105499
22 215147 135 99052 35082 121527
23 344297 108 67654 16116 127766
24 153935 89 65553 15849 98958
25 132943 156 97500 16026 77900
26 174724 129 69112 26569 85646
27 174415 118 82753 24785 98579
28 225548 118 85323 17569 130767
29 223632 125 72654 23825 131741
30 124817 95 30727 7869 53907
31 221698 126 77873 14975 178812
32 210767 135 117478 37791 146761
33 170266 154 74007 9605 82036
34 260561 165 90183 27295 163253
35 84853 113 61542 2746 27032
36 294424 127 101494 34461 171975
37 101011 52 27570 8098 65990
38 215641 121 55813 4787 86572
39 325107 136 79215 24919 159676
40 7176 0 1423 603 1929
41 167542 108 55461 16329 85371
42 106408 46 31081 12558 58391
43 96560 54 22996 7784 31580
44 265769 124 83122 28522 136815
45 269651 115 70106 22265 120642
46 149112 128 60578 14459 69107
47 175824 80 39992 14526 50495
48 152871 97 79892 22240 108016
49 111665 104 49810 11802 46341
50 116408 59 71570 7623 78348
51 362301 125 100708 11912 79336
52 78800 82 33032 7935 56968
53 183167 149 82875 18220 93176
54 277965 149 139077 19199 161632
55 150629 122 71595 19918 87850
56 168809 118 72260 21884 127969
57 24188 12 5950 2694 15049
58 329267 144 115762 15808 155135
59 65029 67 32551 3597 25109
60 101097 52 31701 5296 45824
61 218946 108 80670 25239 102996
62 244052 166 143558 29801 160604
63 341570 80 117105 18450 158051
64 103597 60 23789 7132 44547
65 233328 107 120733 34861 162647
66 256462 127 105195 35940 174141
67 206161 107 73107 16688 60622
68 311473 146 132068 24683 179566
69 235800 84 149193 46230 184301
70 177939 141 46821 10387 75661
71 207176 123 87011 21436 96144
72 196553 111 95260 30546 129847
73 174184 98 55183 19746 117286
74 143246 105 106671 15977 71180
75 187559 135 73511 22583 109377
76 187681 107 92945 17274 85298
77 119016 85 78664 16469 73631
78 182192 155 70054 14251 86767
79 73566 88 22618 3007 23824
80 194979 155 74011 16851 93487
81 167488 104 83737 21113 82981
82 143756 132 69094 17401 73815
83 275541 127 93133 23958 94552
84 243199 108 95536 23567 132190
85 182999 129 225920 13065 128754
86 135649 116 62133 15358 66363
87 152299 122 61370 14587 67808
88 120221 85 43836 12770 61724
89 346485 147 106117 24021 131722
90 145790 99 38692 9648 68580
91 193339 87 84651 20537 106175
92 80953 28 56622 7905 55792
93 122774 90 15986 4527 25157
94 130585 109 95364 30495 76669
95 112611 78 26706 7117 57283
96 286468 111 89691 17719 105805
97 241066 158 67267 27056 129484
98 148446 141 126846 33473 72413
99 204713 122 41140 9758 87831
100 182079 124 102860 21115 96971
101 140344 93 51715 7236 71299
102 220516 124 55801 13790 77494
103 243060 112 111813 32902 120336
104 162765 108 120293 25131 93913
105 182613 99 138599 30910 136048
106 232138 117 161647 35947 181248
107 265318 199 115929 29848 146123
108 85574 78 24266 6943 32036
109 310839 91 162901 42705 186646
110 225060 158 109825 31808 102255
111 232317 126 129838 26675 168237
112 144966 122 37510 8435 64219
113 43287 71 43750 7409 19630
114 155754 75 40652 14993 76825
115 164709 115 87771 36867 115338
116 201940 119 85872 33835 109427
117 235454 124 89275 24164 118168
118 220801 72 44418 12607 84845
119 99466 91 192565 22609 153197
120 92661 45 35232 5892 29877
121 133328 78 40909 17014 63506
122 61361 39 13294 5394 22445
123 125930 68 32387 9178 47695
124 100750 119 140867 6440 68370
125 224549 117 120662 21916 146304
126 82316 39 21233 4011 38233
127 102010 50 44332 5818 42071
128 101523 88 61056 18647 50517
129 243511 155 101338 20556 103950
130 22938 0 1168 238 5841
131 41566 36 13497 70 2341
132 152474 123 65567 22392 84396
133 61857 32 25162 3913 24610
134 99923 99 32334 12237 35753
135 132487 136 40735 8388 55515
136 317394 117 91413 22120 209056
137 21054 0 855 338 6622
138 209641 88 97068 11727 115814
139 22648 39 44339 3704 11609
140 31414 25 14116 3988 13155
141 46698 52 10288 3030 18274
142 131698 75 65622 13520 72875
143 91735 71 16563 1421 10112
144 244749 124 76643 20923 142775
145 184510 151 110681 20237 68847
146 79863 71 29011 3219 17659
147 128423 145 92696 3769 20112
148 97839 87 94785 12252 61023
149 38214 27 8773 1888 13983
150 151101 131 83209 14497 65176
151 272458 162 93815 28864 132432
152 172494 165 86687 21721 112494
153 108043 54 34553 4821 45109
154 328107 159 105547 33644 170875
155 250579 147 103487 15923 180759
156 351067 170 213688 42935 214921
157 158015 119 71220 18864 100226
158 98866 49 23517 4977 32043
159 85439 104 56926 7785 54454
160 229242 120 91721 17939 78876
161 351619 150 115168 23436 170745
162 84207 112 111194 325 6940
163 120445 59 51009 13539 49025
164 324598 136 135777 34538 122037
165 131069 107 51513 12198 53782
166 204271 130 74163 26924 127748
167 165543 115 51633 12716 86839
168 141722 107 75345 8172 44830
169 116048 75 33416 10855 77395
170 250047 71 83305 11932 89324
171 299775 120 98952 14300 103300
172 195838 116 102372 25515 112283
173 173260 79 37238 2805 10901
174 254488 150 103772 29402 120691
175 104389 156 123969 16440 58106
176 136084 51 27142 11221 57140
177 199476 118 135400 28732 122422
178 92499 71 21399 5250 25899
179 224330 144 130115 28608 139296
180 135781 47 24874 8092 52678
181 74408 28 34988 4473 23853
182 81240 68 45549 1572 17306
183 14688 0 6023 2065 7953
184 181633 110 64466 14817 89455
185 271856 147 54990 16714 147866
186 7199 0 1644 556 4245
187 46660 15 6179 2089 21509
188 17547 4 3926 2658 7670
189 133368 64 32755 10695 66675
190 95227 111 34777 1669 14336
191 152601 85 73224 16267 53608
192 98146 68 27114 7768 30059
193 79619 40 20760 7252 29668
194 59194 80 37636 6387 22097
195 139942 88 65461 18715 96841
196 118612 48 30080 7936 41907
197 72880 76 24094 8643 27080
198 65475 51 69008 7294 35885
199 99643 67 54968 4570 41247
200 71965 59 46090 7185 28313
201 77272 61 27507 10058 36845
202 49289 76 10672 2342 16548
203 135131 60 34029 8509 36134
204 108446 68 46300 13275 55764
205 89746 71 24760 6816 28910
206 44296 76 18779 1930 13339
207 77648 62 21280 8086 25319
208 181528 61 40662 10737 66956
209 134019 67 28987 8033 47487
210 124064 88 22827 7058 52785
211 92630 30 18513 6782 44683
212 121848 64 30594 5401 35619
213 52915 68 24006 6521 21920
214 81872 64 27913 10856 45608
215 58981 91 42744 2154 7721
216 53515 88 12934 6117 20634
217 60812 52 22574 5238 29788
218 56375 49 41385 4820 31931
219 65490 62 18653 5615 37754
220 80949 61 18472 4272 32505
221 76302 76 30976 8702 40557
222 104011 88 63339 15340 94238
223 98104 66 25568 8030 44197
224 67989 71 33747 9526 43228
225 30989 68 4154 1278 4103
226 135458 48 19474 4236 44144
227 73504 25 35130 3023 32868
228 63123 68 39067 7196 27640
229 61254 41 13310 3394 14063
230 74914 90 65892 6371 28990
231 31774 66 4143 1574 4694
232 81437 54 28579 9620 42648
233 87186 59 51776 6978 64329
234 50090 60 21152 4911 21928
235 65745 77 38084 8645 25836
236 56653 68 27717 8987 22779
237 158399 72 32928 5544 40820
238 46455 67 11342 3083 27530
239 73624 64 19499 6909 32378
240 38395 63 16380 3189 10824
241 91899 59 36874 6745 39613
242 139526 84 48259 16724 60865
243 52164 64 16734 4850 19787
244 51567 56 28207 7025 20107
245 70551 54 30143 6047 36605
246 84856 67 41369 7377 40961
247 102538 58 45833 9078 48231
248 86678 59 29156 4605 39725
249 85709 40 35944 3238 21455
250 34662 22 36278 8100 23430
251 150580 83 45588 9653 62991
252 99611 81 45097 8914 49363
253 19349 2 3895 786 9604
254 99373 72 28394 6700 24552
255 86230 61 18632 5788 31493
256 30837 15 2325 593 3439
257 31706 32 25139 4506 19555
258 89806 62 27975 6382 21228
259 62088 58 14483 5621 23177
260 40151 36 13127 3997 22094
261 27634 59 5839 520 2342
262 76990 68 24069 8891 38798
263 37460 21 3738 999 3255
264 54157 55 18625 7067 24261
265 49862 54 36341 4639 18511
266 84337 55 24548 5654 40798
267 64175 72 21792 6928 28893
268 59382 41 26263 1514 21425
269 119308 61 23686 9238 50276
270 76702 67 49303 8204 37643
271 103425 76 25659 5926 30377
272 70344 64 28904 5785 27126
273 43410 3 2781 4 13
274 104838 63 29236 5930 42097
275 62215 40 19546 3710 24451
276 69304 69 22818 705 14335
277 53117 48 32689 443 5084
278 19764 8 5752 2416 9927
279 86680 52 22197 7747 43527
280 84105 66 20055 5432 27184
281 77945 76 25272 4913 21610
282 89113 43 82206 2650 20484
283 91005 39 32073 2370 20156
284 40248 14 5444 775 6012
285 64187 61 20154 5576 18475
286 50857 71 36944 1352 12645
287 56613 44 8019 3080 11017
288 62792 60 30884 10205 37623
289 72535 64 19540 6095 35873
tothyperlinks\r
1 144
2 103
3 98
4 135
5 61
6 39
7 150
8 5
9 28
10 84
11 80
12 130
13 82
14 60
15 131
16 84
17 140
18 151
19 91
20 138
21 150
22 124
23 119
24 73
25 110
26 123
27 90
28 116
29 113
30 56
31 115
32 119
33 129
34 127
35 27
36 175
37 35
38 64
39 96
40 0
41 84
42 41
43 47
44 126
45 105
46 80
47 70
48 73
49 57
50 40
51 68
52 21
53 127
54 154
55 116
56 102
57 7
58 148
59 21
60 35
61 112
62 137
63 135
64 26
65 230
66 181
67 71
68 147
69 190
70 64
71 105
72 107
73 94
74 116
75 106
76 143
77 81
78 89
79 26
80 84
81 113
82 120
83 110
84 134
85 54
86 96
87 78
88 51
89 121
90 38
91 145
92 59
93 27
94 91
95 48
96 68
97 58
98 150
99 74
100 181
101 65
102 97
103 121
104 99
105 152
106 188
107 138
108 40
109 254
110 87
111 178
112 51
113 49
114 73
115 176
116 94
117 120
118 66
119 56
120 39
121 66
122 27
123 65
124 58
125 98
126 25
127 26
128 77
129 130
130 11
131 2
132 101
133 31
134 36
135 120
136 195
137 4
138 89
139 24
140 39
141 14
142 78
143 15
144 106
145 83
146 24
147 37
148 77
149 16
150 56
151 132
152 144
153 40
154 153
155 143
156 220
157 79
158 50
159 39
160 95
161 169
162 12
163 63
164 134
165 69
166 119
167 119
168 75
169 63
170 55
171 103
172 197
173 16
174 140
175 89
176 40
177 125
178 21
179 167
180 32
181 36
182 13
183 5
184 96
185 151
186 6
187 13
188 3
189 57
190 23
191 61
192 21
193 43
194 20
195 82
196 90
197 25
198 60
199 61
200 85
201 43
202 25
203 41
204 26
205 38
206 12
207 29
208 49
209 46
210 41
211 31
212 41
213 26
214 23
215 14
216 16
217 25
218 21
219 32
220 9
221 35
222 42
223 68
224 32
225 6
226 68
227 33
228 84
229 46
230 30
231 0
232 36
233 47
234 20
235 50
236 30
237 30
238 34
239 33
240 34
241 37
242 83
243 32
244 30
245 43
246 41
247 51
248 19
249 37
250 33
251 41
252 54
253 14
254 25
255 25
256 8
257 26
258 20
259 11
260 14
261 3
262 40
263 5
264 38
265 32
266 41
267 46
268 47
269 37
270 51
271 49
272 21
273 1
274 44
275 26
276 21
277 4
278 10
279 43
280 34
281 32
282 20
283 34
284 6
285 12
286 24
287 16
288 72
289 27
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) feedback_messages_p1 totsize
14618.1413 505.4920 -0.1344
totrevisions totseconds `tothyperlinks\r`
-0.6224 1.2767 126.8120
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-123899 -18775 -6494 12528 195528
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.462e+04 4.914e+03 2.975 0.00318 **
feedback_messages_p1 5.055e+02 8.016e+01 6.306 1.1e-09 ***
totsize -1.344e-01 9.812e-02 -1.369 0.17195
totrevisions -6.224e-01 4.856e-01 -1.282 0.20103
totseconds 1.277e+00 1.050e-01 12.156 < 2e-16 ***
`tothyperlinks\r` 1.268e+02 9.935e+01 1.276 0.20285
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 34440 on 283 degrees of freedom
Multiple R-squared: 0.8281, Adjusted R-squared: 0.8251
F-statistic: 272.7 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.8551792 2.896416e-01 1.448208e-01
[2,] 0.7521185 4.957629e-01 2.478815e-01
[3,] 0.7021119 5.957761e-01 2.978881e-01
[4,] 0.7374555 5.250891e-01 2.625445e-01
[5,] 0.6674386 6.651227e-01 3.325614e-01
[6,] 0.6541736 6.916529e-01 3.458264e-01
[7,] 0.9148944 1.702112e-01 8.510562e-02
[8,] 0.8755567 2.488866e-01 1.244433e-01
[9,] 0.8503674 2.992653e-01 1.496326e-01
[10,] 0.8044548 3.910904e-01 1.955452e-01
[11,] 0.8510258 2.979484e-01 1.489742e-01
[12,] 0.8295912 3.408177e-01 1.704088e-01
[13,] 0.8055837 3.888327e-01 1.944163e-01
[14,] 0.7753430 4.493139e-01 2.246570e-01
[15,] 0.9590082 8.198358e-02 4.099179e-02
[16,] 0.9540909 9.181820e-02 4.590910e-02
[17,] 0.9445964 1.108072e-01 5.540361e-02
[18,] 0.9264197 1.471606e-01 7.358031e-02
[19,] 0.9038105 1.923789e-01 9.618947e-02
[20,] 0.8849823 2.300354e-01 1.150177e-01
[21,] 0.8792112 2.415775e-01 1.207888e-01
[22,] 0.8469384 3.061232e-01 1.530616e-01
[23,] 0.9596165 8.076698e-02 4.038349e-02
[24,] 0.9571106 8.577883e-02 4.288941e-02
[25,] 0.9466998 1.066004e-01 5.330020e-02
[26,] 0.9389004 1.221992e-01 6.109958e-02
[27,] 0.9327470 1.345061e-01 6.725304e-02
[28,] 0.9161117 1.677766e-01 8.388831e-02
[29,] 0.9015117 1.969765e-01 9.848827e-02
[30,] 0.9197491 1.605018e-01 8.025089e-02
[31,] 0.9397991 1.204018e-01 6.020090e-02
[32,] 0.9247329 1.505342e-01 7.526712e-02
[33,] 0.9055261 1.889479e-01 9.447393e-02
[34,] 0.8831393 2.337214e-01 1.168607e-01
[35,] 0.8624317 2.751367e-01 1.375683e-01
[36,] 0.8431438 3.137125e-01 1.568562e-01
[37,] 0.8640430 2.719140e-01 1.359570e-01
[38,] 0.8371803 3.256394e-01 1.628197e-01
[39,] 0.8766785 2.466430e-01 1.233215e-01
[40,] 0.8680788 2.638425e-01 1.319212e-01
[41,] 0.8421775 3.156450e-01 1.578225e-01
[42,] 0.8147155 3.705691e-01 1.852845e-01
[43,] 0.9999879 2.421718e-05 1.210859e-05
[44,] 0.9999887 2.264645e-05 1.132323e-05
[45,] 0.9999838 3.234917e-05 1.617458e-05
[46,] 0.9999756 4.871999e-05 2.435999e-05
[47,] 0.9999709 5.823486e-05 2.911743e-05
[48,] 0.9999842 3.167455e-05 1.583727e-05
[49,] 0.9999764 4.711673e-05 2.355836e-05
[50,] 0.9999810 3.791902e-05 1.895951e-05
[51,] 0.9999714 5.729901e-05 2.864951e-05
[52,] 0.9999569 8.618180e-05 4.309090e-05
[53,] 0.9999505 9.907577e-05 4.953788e-05
[54,] 0.9999542 9.163873e-05 4.581937e-05
[55,] 0.9999917 1.662155e-05 8.310775e-06
[56,] 0.9999875 2.502848e-05 1.251424e-05
[57,] 0.9999898 2.038649e-05 1.019324e-05
[58,] 0.9999885 2.304751e-05 1.152376e-05
[59,] 0.9999967 6.541838e-06 3.270919e-06
[60,] 0.9999950 9.924054e-06 4.962027e-06
[61,] 0.9999949 1.014255e-05 5.071277e-06
[62,] 0.9999922 1.560201e-05 7.801003e-06
[63,] 0.9999895 2.098749e-05 1.049374e-05
[64,] 0.9999861 2.788714e-05 1.394357e-05
[65,] 0.9999840 3.206404e-05 1.603202e-05
[66,] 0.9999766 4.688980e-05 2.344490e-05
[67,] 0.9999696 6.072933e-05 3.036466e-05
[68,] 0.9999579 8.425662e-05 4.212831e-05
[69,] 0.9999471 1.058222e-04 5.291109e-05
[70,] 0.9999261 1.477139e-04 7.385695e-05
[71,] 0.9998976 2.048141e-04 1.024071e-04
[72,] 0.9998548 2.904523e-04 1.452261e-04
[73,] 0.9997956 4.087317e-04 2.043659e-04
[74,] 0.9997516 4.968879e-04 2.484439e-04
[75,] 0.9999633 7.334563e-05 3.667281e-05
[76,] 0.9999499 1.002542e-04 5.012708e-05
[77,] 0.9999555 8.891983e-05 4.445992e-05
[78,] 0.9999401 1.198889e-04 5.994444e-05
[79,] 0.9999140 1.719990e-04 8.599952e-05
[80,] 0.9998794 2.412557e-04 1.206279e-04
[81,] 0.9999925 1.503413e-05 7.517063e-06
[82,] 0.9999888 2.244363e-05 1.122181e-05
[83,] 0.9999835 3.297592e-05 1.648796e-05
[84,] 0.9999772 4.557094e-05 2.278547e-05
[85,] 0.9999761 4.771024e-05 2.385512e-05
[86,] 0.9999679 6.410331e-05 3.205166e-05
[87,] 0.9999555 8.901881e-05 4.450941e-05
[88,] 0.9999961 7.819626e-06 3.909813e-06
[89,] 0.9999941 1.173912e-05 5.869559e-06
[90,] 0.9999918 1.637202e-05 8.186012e-06
[91,] 0.9999894 2.127427e-05 1.063714e-05
[92,] 0.9999855 2.901406e-05 1.450703e-05
[93,] 0.9999793 4.146710e-05 2.073355e-05
[94,] 0.9999855 2.907753e-05 1.453876e-05
[95,] 0.9999864 2.723202e-05 1.361601e-05
[96,] 0.9999804 3.915302e-05 1.957651e-05
[97,] 0.9999815 3.691777e-05 1.845888e-05
[98,] 0.9999888 2.247517e-05 1.123758e-05
[99,] 0.9999853 2.946813e-05 1.473406e-05
[100,] 0.9999788 4.239324e-05 2.119662e-05
[101,] 0.9999750 5.000026e-05 2.500013e-05
[102,] 0.9999695 6.103307e-05 3.051653e-05
[103,] 0.9999787 4.269608e-05 2.134804e-05
[104,] 0.9999698 6.037306e-05 3.018653e-05
[105,] 0.9999661 6.783439e-05 3.391720e-05
[106,] 0.9999532 9.362636e-05 4.681318e-05
[107,] 0.9999672 6.563868e-05 3.281934e-05
[108,] 0.9999536 9.278436e-05 4.639218e-05
[109,] 0.9999401 1.197353e-04 5.986764e-05
[110,] 0.9999792 4.163386e-05 2.081693e-05
[111,] 0.9999999 1.914901e-07 9.574504e-08
[112,] 0.9999999 2.612705e-07 1.306353e-07
[113,] 0.9999998 4.081050e-07 2.040525e-07
[114,] 0.9999997 6.368737e-07 3.184369e-07
[115,] 0.9999996 8.621036e-07 4.310518e-07
[116,] 0.9999998 3.881104e-07 1.940552e-07
[117,] 0.9999998 4.337595e-07 2.168798e-07
[118,] 0.9999997 6.765307e-07 3.382653e-07
[119,] 0.9999995 9.941277e-07 4.970638e-07
[120,] 0.9999993 1.414627e-06 7.073135e-07
[121,] 0.9999992 1.588699e-06 7.943495e-07
[122,] 0.9999988 2.423913e-06 1.211957e-06
[123,] 0.9999982 3.615174e-06 1.807587e-06
[124,] 0.9999977 4.627631e-06 2.313815e-06
[125,] 0.9999965 6.943789e-06 3.471895e-06
[126,] 0.9999948 1.030834e-05 5.154168e-06
[127,] 0.9999935 1.301446e-05 6.507230e-06
[128,] 0.9999919 1.615285e-05 8.076423e-06
[129,] 0.9999882 2.356949e-05 1.178474e-05
[130,] 0.9999833 3.333235e-05 1.666618e-05
[131,] 0.9999813 3.748928e-05 1.874464e-05
[132,] 0.9999745 5.097814e-05 2.548907e-05
[133,] 0.9999657 6.856355e-05 3.428177e-05
[134,] 0.9999529 9.413877e-05 4.706938e-05
[135,] 0.9999534 9.314462e-05 4.657231e-05
[136,] 0.9999345 1.309980e-04 6.549902e-05
[137,] 0.9999186 1.628234e-04 8.141170e-05
[138,] 0.9998900 2.199459e-04 1.099729e-04
[139,] 0.9998731 2.537082e-04 1.268541e-04
[140,] 0.9998871 2.257569e-04 1.128784e-04
[141,] 0.9998453 3.094906e-04 1.547453e-04
[142,] 0.9997855 4.290322e-04 2.145161e-04
[143,] 0.9997497 5.005117e-04 2.502559e-04
[144,] 0.9998685 2.629563e-04 1.314781e-04
[145,] 0.9998216 3.568872e-04 1.784436e-04
[146,] 0.9998147 3.705776e-04 1.852888e-04
[147,] 0.9999480 1.039646e-04 5.198232e-05
[148,] 0.9999559 8.828009e-05 4.414005e-05
[149,] 0.9999612 7.760678e-05 3.880339e-05
[150,] 0.9999508 9.834464e-05 4.917232e-05
[151,] 0.9999693 6.142669e-05 3.071335e-05
[152,] 0.9999872 2.559156e-05 1.279578e-05
[153,] 0.9999887 2.256948e-05 1.128474e-05
[154,] 0.9999854 2.910913e-05 1.455457e-05
[155,] 0.9999812 3.763909e-05 1.881955e-05
[156,] 0.9999998 4.563932e-07 2.281966e-07
[157,] 0.9999996 7.033055e-07 3.516527e-07
[158,] 0.9999995 9.466478e-07 4.733239e-07
[159,] 0.9999994 1.297623e-06 6.488113e-07
[160,] 0.9999991 1.799316e-06 8.996579e-07
[161,] 0.9999992 1.556912e-06 7.784561e-07
[162,] 0.9999999 1.179734e-07 5.898670e-08
[163,] 1.0000000 1.185999e-09 5.929993e-10
[164,] 1.0000000 1.961454e-09 9.807270e-10
[165,] 1.0000000 2.429080e-13 1.214540e-13
[166,] 1.0000000 1.137499e-13 5.687496e-14
[167,] 1.0000000 7.743489e-14 3.871745e-14
[168,] 1.0000000 6.328346e-14 3.164173e-14
[169,] 1.0000000 1.257159e-13 6.285794e-14
[170,] 1.0000000 1.874556e-13 9.372779e-14
[171,] 1.0000000 2.066695e-13 1.033348e-13
[172,] 1.0000000 1.256378e-13 6.281888e-14
[173,] 1.0000000 2.236812e-13 1.118406e-13
[174,] 1.0000000 3.830103e-13 1.915051e-13
[175,] 1.0000000 6.962177e-13 3.481089e-13
[176,] 1.0000000 1.372487e-12 6.862434e-13
[177,] 1.0000000 2.735081e-12 1.367541e-12
[178,] 1.0000000 4.245574e-12 2.122787e-12
[179,] 1.0000000 8.238048e-12 4.119024e-12
[180,] 1.0000000 1.494662e-11 7.473309e-12
[181,] 1.0000000 2.783045e-11 1.391522e-11
[182,] 1.0000000 3.696865e-11 1.848432e-11
[183,] 1.0000000 1.064470e-11 5.322349e-12
[184,] 1.0000000 1.113611e-11 5.568054e-12
[185,] 1.0000000 2.042101e-11 1.021050e-11
[186,] 1.0000000 3.759930e-11 1.879965e-11
[187,] 1.0000000 3.791327e-11 1.895664e-11
[188,] 1.0000000 5.197766e-11 2.598883e-11
[189,] 1.0000000 9.772241e-11 4.886121e-11
[190,] 1.0000000 1.241873e-10 6.209367e-11
[191,] 1.0000000 2.342487e-10 1.171243e-10
[192,] 1.0000000 3.965276e-10 1.982638e-10
[193,] 1.0000000 7.483104e-10 3.741552e-10
[194,] 1.0000000 1.117800e-09 5.589001e-10
[195,] 1.0000000 1.650484e-10 8.252418e-11
[196,] 1.0000000 2.857804e-10 1.428902e-10
[197,] 1.0000000 4.668404e-10 2.334202e-10
[198,] 1.0000000 7.045435e-10 3.522718e-10
[199,] 1.0000000 1.054688e-09 5.273441e-10
[200,] 1.0000000 5.625000e-11 2.812500e-11
[201,] 1.0000000 3.019983e-11 1.509992e-11
[202,] 1.0000000 5.408096e-11 2.704048e-11
[203,] 1.0000000 9.602840e-11 4.801420e-11
[204,] 1.0000000 4.988145e-11 2.494072e-11
[205,] 1.0000000 9.324951e-11 4.662476e-11
[206,] 1.0000000 1.857915e-10 9.289576e-11
[207,] 1.0000000 3.800967e-10 1.900484e-10
[208,] 1.0000000 6.584435e-10 3.292218e-10
[209,] 1.0000000 1.213086e-09 6.065429e-10
[210,] 1.0000000 1.693339e-09 8.466697e-10
[211,] 1.0000000 2.305586e-09 1.152793e-09
[212,] 1.0000000 4.520823e-09 2.260412e-09
[213,] 1.0000000 7.914610e-09 3.957305e-09
[214,] 1.0000000 5.588139e-10 2.794069e-10
[215,] 1.0000000 1.178389e-09 5.891944e-10
[216,] 1.0000000 1.144996e-09 5.724982e-10
[217,] 1.0000000 1.966504e-09 9.832518e-10
[218,] 1.0000000 4.547874e-10 2.273937e-10
[219,] 1.0000000 9.975623e-10 4.987812e-10
[220,] 1.0000000 2.017901e-09 1.008951e-09
[221,] 1.0000000 2.894259e-09 1.447129e-09
[222,] 1.0000000 3.209979e-09 1.604990e-09
[223,] 1.0000000 4.907911e-09 2.453956e-09
[224,] 1.0000000 1.003967e-08 5.019834e-09
[225,] 1.0000000 9.549148e-10 4.774574e-10
[226,] 1.0000000 1.394414e-09 6.972071e-10
[227,] 1.0000000 3.070003e-09 1.535002e-09
[228,] 1.0000000 6.361062e-09 3.180531e-09
[229,] 1.0000000 1.040864e-10 5.204319e-11
[230,] 1.0000000 6.360539e-11 3.180269e-11
[231,] 1.0000000 1.560285e-10 7.801424e-11
[232,] 1.0000000 3.102768e-10 1.551384e-10
[233,] 1.0000000 7.860069e-10 3.930034e-10
[234,] 1.0000000 2.635910e-10 1.317955e-10
[235,] 1.0000000 5.752014e-10 2.876007e-10
[236,] 1.0000000 1.419071e-09 7.095357e-10
[237,] 1.0000000 2.585012e-09 1.292506e-09
[238,] 1.0000000 4.883585e-09 2.441792e-09
[239,] 1.0000000 1.226037e-08 6.130183e-09
[240,] 1.0000000 1.771370e-08 8.856852e-09
[241,] 1.0000000 1.610060e-08 8.050300e-09
[242,] 1.0000000 3.256412e-08 1.628206e-08
[243,] 1.0000000 5.928379e-08 2.964189e-08
[244,] 0.9999999 1.211819e-07 6.059095e-08
[245,] 0.9999999 2.294297e-07 1.147148e-07
[246,] 0.9999999 1.329994e-07 6.649971e-08
[247,] 0.9999998 3.145292e-07 1.572646e-07
[248,] 0.9999996 8.033498e-07 4.016749e-07
[249,] 0.9999997 6.551614e-07 3.275807e-07
[250,] 0.9999999 2.949312e-07 1.474656e-07
[251,] 0.9999996 8.092950e-07 4.046475e-07
[252,] 0.9999997 6.706172e-07 3.353086e-07
[253,] 0.9999996 8.755408e-07 4.377704e-07
[254,] 0.9999989 2.287467e-06 1.143733e-06
[255,] 0.9999969 6.295950e-06 3.147975e-06
[256,] 0.9999917 1.656472e-05 8.282361e-06
[257,] 0.9999835 3.290371e-05 1.645186e-05
[258,] 0.9999674 6.519556e-05 3.259778e-05
[259,] 0.9999260 1.479515e-04 7.397574e-05
[260,] 0.9998743 2.513766e-04 1.256883e-04
[261,] 0.9998765 2.470288e-04 1.235144e-04
[262,] 0.9997446 5.107478e-04 2.553739e-04
[263,] 0.9997603 4.793187e-04 2.396594e-04
[264,] 0.9993439 1.312222e-03 6.561109e-04
[265,] 0.9989108 2.178481e-03 1.089240e-03
[266,] 0.9973980 5.203905e-03 2.601953e-03
[267,] 0.9935360 1.292794e-02 6.463971e-03
[268,] 0.9843067 3.138656e-02 1.569328e-02
[269,] 0.9640905 7.181894e-02 3.590947e-02
[270,] 0.9602075 7.958508e-02 3.979254e-02
[271,] 0.9052998 1.894004e-01 9.470018e-02
[272,] 0.8291404 3.417192e-01 1.708596e-01
> postscript(file="/var/wessaorg/rcomp/tmp/1162j1324143944.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/2ioyg1324143944.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/3aqj51324143944.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/4xeii1324143944.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/5h5n61324143944.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
-36726.19374 -64615.62528 -14361.63119 -1119.57500 -24162.07284
6 7 8 9 10
-57240.15403 90266.66998 -20844.06640 39336.30912 -16537.54058
11 12 13 14 15
2338.83722 -27321.44011 -11760.42687 -1386.25821 60853.23138
16 17 18 19 20
-20783.97495 38976.60406 -408.97338 33197.22406 -34659.53386
21 22 23 24 25
-10304.58846 -3451.02317 115993.05436 -22599.70131 -50863.58111
26 27 28 29 30
-4225.11484 -10578.27125 -7983.46935 -12109.73614 -4722.88606
31 32 33 34 35
-79706.52394 -35252.11326 -27372.49982 -32893.14372 -14844.00823
36 37 38 39 40
8935.68431 -19838.07024 31692.01867 51857.67063 -9338.45843
41 42 43 44 45
-3702.34462 780.41921 16300.44608 26735.90574 52835.83359
46 47 48 49 50
-11446.49393 61835.36613 -33367.99113 -7879.48745 -18774.90594
51 52 53 54 55
195527.99757 -43287.34856 -19359.88142 -7224.63759 -30513.66706
56 57 58 59 60
-58444.50224 -14121.07971 50417.38951 -11565.12679 4805.63521
61 62 63 64 65
30581.14392 -39062.23117 94822.25479 6112.95992 -34281.77825
66 67 68 69 70
-31134.80469 71262.94012 8261.97404 -31857.46280 87.63759
71 72 73 74 75
19349.73122 -21712.39901 -31931.07118 -5759.95982 -24455.03002
76 77 78 79 80
15178.17166 -22027.95419 -14559.32712 -14338.78254 -7568.03295
81 82 83 84 85
4416.23232 -26932.49939 89483.47079 15728.32359 -29572.14820
86 87 88 89 90
-16600.78278 -3128.36439 -8798.44385 103250.48627 -45.26150
91 92 93 94 95
4954.47517 -14004.16924 32084.42997 -16764.27137 -12639.57546
96 97 98 99 100
95111.93807 -213.68220 -11043.61619 18505.26173 -15016.53557
101 102 103 104 105
-9104.96996 48057.50563 38347.19072 -7097.98645 -37160.18704
106 107 108 109 110
-52775.68033 -19799.11663 -6864.63780 28181.22642 23542.71336
111 112 113 114 115
-49311.15894 -9490.07575 -28007.37667 10675.43362 -42876.57757
116 117 118 119 120
8135.58822 19103.40790 66908.26584 -123899.03729 20606.21775
121 122 123 124 125
5917.74806 91.89708 17865.86842 -45730.92117 -18577.08243
126 127 128 129 130
1349.50929 14684.55211 -12030.25242 27749.91501 -227.38513
131 132 133 134 135
7364.80607 -22132.36640 1527.89985 -2990.41814 -26279.25712
136 137 138 139 140
-21953.28665 -2200.64645 11731.38221 -21286.43397 -13203.74813
141 142 143 144 145
-16043.89921 -6533.17274 29524.35706 -4957.13114 22604.92743
146 147 148 149 150
9667.18999 24939.83015 -28070.11123 -7580.09059 152.83501
151 152 153 154 155
20700.87216 -62249.28626 11107.03861 30673.31311 -63445.71602
156 157 158 159 160
3654.64495 -33426.31991 18485.34960 -43724.89806 64703.31965
161 162 163 164 165
51811.28369 17734.50395 20702.20698 108171.26543 -538.56569
166 167 168 169 170
-27529.59598 -18315.43309 21479.39061 -32037.83494 97141.21982
171 172 173 174 175
101745.77266 -16119.00355 109510.65655 24445.08328 -47668.74259
176 177 178 179 180
28291.57846 -10866.26239 12404.61751 -26811.99110 34469.66254
181 182 183 184 185
18102.10445 15603.38577 -8623.54272 2910.70592 -7211.81642
186 187 188 189 190
-13032.79222 -2519.86545 -7084.27172 5101.59329 8990.97110
191 192 193 194 195
38800.56704 16591.93679 8753.14389 -17579.52992 -32754.32259
196 197 198 199 200
23794.15017 -9283.08061 -14535.45325 990.19278 -8739.53924
201 202 203 204 205
-10719.32915 -25152.60042 48718.82345 -555.09565 5077.82232
206 207 208 209 210
-23567.12976 3578.11226 56522.37831 27965.88550 -169.00990
211 212 213 214 215
8576.27889 31675.48864 -20075.53116 -15736.34004 -6185.96791
216 217 218 219 220
-28414.52055 -15000.07055 -17881.96471 -26727.31742 -2004.76387
221 222 223 224 225
-23374.25418 -62675.14987 -6494.31144 -31304.33838 -22648.34273
226 227 228 229 230
36846.10194 6701.94019 -22081.72525 6023.43656 -13196.36940
231 232 233 234 235
-20663.28969 -9665.62064 -34047.19671 -19491.43823 -16624.59090
236 237 238 239 240
-15908.08542 59339.82437 -38048.31915 -11948.39187 -22014.39972
241 242 243 244 245
1342.22851 11105.93924 -18859.26737 -12672.01540 -15737.61668
246 247 248 249 250
-10975.67090 2364.28163 -4108.06634 25631.80931 -15259.77291
251 252 253 254 255
20517.51812 -14215.69160 -9304.66092 21827.96483 3504.31029
256 257 258 259 260
3912.77846 -21169.19638 21939.57979 -7389.98596 -18396.83285
261 262 263 264 265
-19070.50909 -17841.05035 8560.72325 -17155.91168 -11974.04042
266 267 268 269 270
-8553.18336 -22320.55457 -4804.27121 13906.03012 -14580.83776
271 272 273 274 275
12528.35207 -6437.09671 27508.13540 6666.65700 -2201.91950
276 277 278 279 280
2346.68337 11905.11214 -10563.76331 -7444.79525 3181.60823
281 282 283 284 285
-285.22052 36764.94050 32411.85739 11330.22991 -197.10752
286 287 288 289
-13033.33015 6652.88026 -28819.43591 -17239.79177
> postscript(file="/var/wessaorg/rcomp/tmp/6d4x11324143944.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 -36726.19374 NA
1 -64615.62528 -36726.19374
2 -14361.63119 -64615.62528
3 -1119.57500 -14361.63119
4 -24162.07284 -1119.57500
5 -57240.15403 -24162.07284
6 90266.66998 -57240.15403
7 -20844.06640 90266.66998
8 39336.30912 -20844.06640
9 -16537.54058 39336.30912
10 2338.83722 -16537.54058
11 -27321.44011 2338.83722
12 -11760.42687 -27321.44011
13 -1386.25821 -11760.42687
14 60853.23138 -1386.25821
15 -20783.97495 60853.23138
16 38976.60406 -20783.97495
17 -408.97338 38976.60406
18 33197.22406 -408.97338
19 -34659.53386 33197.22406
20 -10304.58846 -34659.53386
21 -3451.02317 -10304.58846
22 115993.05436 -3451.02317
23 -22599.70131 115993.05436
24 -50863.58111 -22599.70131
25 -4225.11484 -50863.58111
26 -10578.27125 -4225.11484
27 -7983.46935 -10578.27125
28 -12109.73614 -7983.46935
29 -4722.88606 -12109.73614
30 -79706.52394 -4722.88606
31 -35252.11326 -79706.52394
32 -27372.49982 -35252.11326
33 -32893.14372 -27372.49982
34 -14844.00823 -32893.14372
35 8935.68431 -14844.00823
36 -19838.07024 8935.68431
37 31692.01867 -19838.07024
38 51857.67063 31692.01867
39 -9338.45843 51857.67063
40 -3702.34462 -9338.45843
41 780.41921 -3702.34462
42 16300.44608 780.41921
43 26735.90574 16300.44608
44 52835.83359 26735.90574
45 -11446.49393 52835.83359
46 61835.36613 -11446.49393
47 -33367.99113 61835.36613
48 -7879.48745 -33367.99113
49 -18774.90594 -7879.48745
50 195527.99757 -18774.90594
51 -43287.34856 195527.99757
52 -19359.88142 -43287.34856
53 -7224.63759 -19359.88142
54 -30513.66706 -7224.63759
55 -58444.50224 -30513.66706
56 -14121.07971 -58444.50224
57 50417.38951 -14121.07971
58 -11565.12679 50417.38951
59 4805.63521 -11565.12679
60 30581.14392 4805.63521
61 -39062.23117 30581.14392
62 94822.25479 -39062.23117
63 6112.95992 94822.25479
64 -34281.77825 6112.95992
65 -31134.80469 -34281.77825
66 71262.94012 -31134.80469
67 8261.97404 71262.94012
68 -31857.46280 8261.97404
69 87.63759 -31857.46280
70 19349.73122 87.63759
71 -21712.39901 19349.73122
72 -31931.07118 -21712.39901
73 -5759.95982 -31931.07118
74 -24455.03002 -5759.95982
75 15178.17166 -24455.03002
76 -22027.95419 15178.17166
77 -14559.32712 -22027.95419
78 -14338.78254 -14559.32712
79 -7568.03295 -14338.78254
80 4416.23232 -7568.03295
81 -26932.49939 4416.23232
82 89483.47079 -26932.49939
83 15728.32359 89483.47079
84 -29572.14820 15728.32359
85 -16600.78278 -29572.14820
86 -3128.36439 -16600.78278
87 -8798.44385 -3128.36439
88 103250.48627 -8798.44385
89 -45.26150 103250.48627
90 4954.47517 -45.26150
91 -14004.16924 4954.47517
92 32084.42997 -14004.16924
93 -16764.27137 32084.42997
94 -12639.57546 -16764.27137
95 95111.93807 -12639.57546
96 -213.68220 95111.93807
97 -11043.61619 -213.68220
98 18505.26173 -11043.61619
99 -15016.53557 18505.26173
100 -9104.96996 -15016.53557
101 48057.50563 -9104.96996
102 38347.19072 48057.50563
103 -7097.98645 38347.19072
104 -37160.18704 -7097.98645
105 -52775.68033 -37160.18704
106 -19799.11663 -52775.68033
107 -6864.63780 -19799.11663
108 28181.22642 -6864.63780
109 23542.71336 28181.22642
110 -49311.15894 23542.71336
111 -9490.07575 -49311.15894
112 -28007.37667 -9490.07575
113 10675.43362 -28007.37667
114 -42876.57757 10675.43362
115 8135.58822 -42876.57757
116 19103.40790 8135.58822
117 66908.26584 19103.40790
118 -123899.03729 66908.26584
119 20606.21775 -123899.03729
120 5917.74806 20606.21775
121 91.89708 5917.74806
122 17865.86842 91.89708
123 -45730.92117 17865.86842
124 -18577.08243 -45730.92117
125 1349.50929 -18577.08243
126 14684.55211 1349.50929
127 -12030.25242 14684.55211
128 27749.91501 -12030.25242
129 -227.38513 27749.91501
130 7364.80607 -227.38513
131 -22132.36640 7364.80607
132 1527.89985 -22132.36640
133 -2990.41814 1527.89985
134 -26279.25712 -2990.41814
135 -21953.28665 -26279.25712
136 -2200.64645 -21953.28665
137 11731.38221 -2200.64645
138 -21286.43397 11731.38221
139 -13203.74813 -21286.43397
140 -16043.89921 -13203.74813
141 -6533.17274 -16043.89921
142 29524.35706 -6533.17274
143 -4957.13114 29524.35706
144 22604.92743 -4957.13114
145 9667.18999 22604.92743
146 24939.83015 9667.18999
147 -28070.11123 24939.83015
148 -7580.09059 -28070.11123
149 152.83501 -7580.09059
150 20700.87216 152.83501
151 -62249.28626 20700.87216
152 11107.03861 -62249.28626
153 30673.31311 11107.03861
154 -63445.71602 30673.31311
155 3654.64495 -63445.71602
156 -33426.31991 3654.64495
157 18485.34960 -33426.31991
158 -43724.89806 18485.34960
159 64703.31965 -43724.89806
160 51811.28369 64703.31965
161 17734.50395 51811.28369
162 20702.20698 17734.50395
163 108171.26543 20702.20698
164 -538.56569 108171.26543
165 -27529.59598 -538.56569
166 -18315.43309 -27529.59598
167 21479.39061 -18315.43309
168 -32037.83494 21479.39061
169 97141.21982 -32037.83494
170 101745.77266 97141.21982
171 -16119.00355 101745.77266
172 109510.65655 -16119.00355
173 24445.08328 109510.65655
174 -47668.74259 24445.08328
175 28291.57846 -47668.74259
176 -10866.26239 28291.57846
177 12404.61751 -10866.26239
178 -26811.99110 12404.61751
179 34469.66254 -26811.99110
180 18102.10445 34469.66254
181 15603.38577 18102.10445
182 -8623.54272 15603.38577
183 2910.70592 -8623.54272
184 -7211.81642 2910.70592
185 -13032.79222 -7211.81642
186 -2519.86545 -13032.79222
187 -7084.27172 -2519.86545
188 5101.59329 -7084.27172
189 8990.97110 5101.59329
190 38800.56704 8990.97110
191 16591.93679 38800.56704
192 8753.14389 16591.93679
193 -17579.52992 8753.14389
194 -32754.32259 -17579.52992
195 23794.15017 -32754.32259
196 -9283.08061 23794.15017
197 -14535.45325 -9283.08061
198 990.19278 -14535.45325
199 -8739.53924 990.19278
200 -10719.32915 -8739.53924
201 -25152.60042 -10719.32915
202 48718.82345 -25152.60042
203 -555.09565 48718.82345
204 5077.82232 -555.09565
205 -23567.12976 5077.82232
206 3578.11226 -23567.12976
207 56522.37831 3578.11226
208 27965.88550 56522.37831
209 -169.00990 27965.88550
210 8576.27889 -169.00990
211 31675.48864 8576.27889
212 -20075.53116 31675.48864
213 -15736.34004 -20075.53116
214 -6185.96791 -15736.34004
215 -28414.52055 -6185.96791
216 -15000.07055 -28414.52055
217 -17881.96471 -15000.07055
218 -26727.31742 -17881.96471
219 -2004.76387 -26727.31742
220 -23374.25418 -2004.76387
221 -62675.14987 -23374.25418
222 -6494.31144 -62675.14987
223 -31304.33838 -6494.31144
224 -22648.34273 -31304.33838
225 36846.10194 -22648.34273
226 6701.94019 36846.10194
227 -22081.72525 6701.94019
228 6023.43656 -22081.72525
229 -13196.36940 6023.43656
230 -20663.28969 -13196.36940
231 -9665.62064 -20663.28969
232 -34047.19671 -9665.62064
233 -19491.43823 -34047.19671
234 -16624.59090 -19491.43823
235 -15908.08542 -16624.59090
236 59339.82437 -15908.08542
237 -38048.31915 59339.82437
238 -11948.39187 -38048.31915
239 -22014.39972 -11948.39187
240 1342.22851 -22014.39972
241 11105.93924 1342.22851
242 -18859.26737 11105.93924
243 -12672.01540 -18859.26737
244 -15737.61668 -12672.01540
245 -10975.67090 -15737.61668
246 2364.28163 -10975.67090
247 -4108.06634 2364.28163
248 25631.80931 -4108.06634
249 -15259.77291 25631.80931
250 20517.51812 -15259.77291
251 -14215.69160 20517.51812
252 -9304.66092 -14215.69160
253 21827.96483 -9304.66092
254 3504.31029 21827.96483
255 3912.77846 3504.31029
256 -21169.19638 3912.77846
257 21939.57979 -21169.19638
258 -7389.98596 21939.57979
259 -18396.83285 -7389.98596
260 -19070.50909 -18396.83285
261 -17841.05035 -19070.50909
262 8560.72325 -17841.05035
263 -17155.91168 8560.72325
264 -11974.04042 -17155.91168
265 -8553.18336 -11974.04042
266 -22320.55457 -8553.18336
267 -4804.27121 -22320.55457
268 13906.03012 -4804.27121
269 -14580.83776 13906.03012
270 12528.35207 -14580.83776
271 -6437.09671 12528.35207
272 27508.13540 -6437.09671
273 6666.65700 27508.13540
274 -2201.91950 6666.65700
275 2346.68337 -2201.91950
276 11905.11214 2346.68337
277 -10563.76331 11905.11214
278 -7444.79525 -10563.76331
279 3181.60823 -7444.79525
280 -285.22052 3181.60823
281 36764.94050 -285.22052
282 32411.85739 36764.94050
283 11330.22991 32411.85739
284 -197.10752 11330.22991
285 -13033.33015 -197.10752
286 6652.88026 -13033.33015
287 -28819.43591 6652.88026
288 -17239.79177 -28819.43591
289 NA -17239.79177
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -64615.62528 -36726.19374
[2,] -14361.63119 -64615.62528
[3,] -1119.57500 -14361.63119
[4,] -24162.07284 -1119.57500
[5,] -57240.15403 -24162.07284
[6,] 90266.66998 -57240.15403
[7,] -20844.06640 90266.66998
[8,] 39336.30912 -20844.06640
[9,] -16537.54058 39336.30912
[10,] 2338.83722 -16537.54058
[11,] -27321.44011 2338.83722
[12,] -11760.42687 -27321.44011
[13,] -1386.25821 -11760.42687
[14,] 60853.23138 -1386.25821
[15,] -20783.97495 60853.23138
[16,] 38976.60406 -20783.97495
[17,] -408.97338 38976.60406
[18,] 33197.22406 -408.97338
[19,] -34659.53386 33197.22406
[20,] -10304.58846 -34659.53386
[21,] -3451.02317 -10304.58846
[22,] 115993.05436 -3451.02317
[23,] -22599.70131 115993.05436
[24,] -50863.58111 -22599.70131
[25,] -4225.11484 -50863.58111
[26,] -10578.27125 -4225.11484
[27,] -7983.46935 -10578.27125
[28,] -12109.73614 -7983.46935
[29,] -4722.88606 -12109.73614
[30,] -79706.52394 -4722.88606
[31,] -35252.11326 -79706.52394
[32,] -27372.49982 -35252.11326
[33,] -32893.14372 -27372.49982
[34,] -14844.00823 -32893.14372
[35,] 8935.68431 -14844.00823
[36,] -19838.07024 8935.68431
[37,] 31692.01867 -19838.07024
[38,] 51857.67063 31692.01867
[39,] -9338.45843 51857.67063
[40,] -3702.34462 -9338.45843
[41,] 780.41921 -3702.34462
[42,] 16300.44608 780.41921
[43,] 26735.90574 16300.44608
[44,] 52835.83359 26735.90574
[45,] -11446.49393 52835.83359
[46,] 61835.36613 -11446.49393
[47,] -33367.99113 61835.36613
[48,] -7879.48745 -33367.99113
[49,] -18774.90594 -7879.48745
[50,] 195527.99757 -18774.90594
[51,] -43287.34856 195527.99757
[52,] -19359.88142 -43287.34856
[53,] -7224.63759 -19359.88142
[54,] -30513.66706 -7224.63759
[55,] -58444.50224 -30513.66706
[56,] -14121.07971 -58444.50224
[57,] 50417.38951 -14121.07971
[58,] -11565.12679 50417.38951
[59,] 4805.63521 -11565.12679
[60,] 30581.14392 4805.63521
[61,] -39062.23117 30581.14392
[62,] 94822.25479 -39062.23117
[63,] 6112.95992 94822.25479
[64,] -34281.77825 6112.95992
[65,] -31134.80469 -34281.77825
[66,] 71262.94012 -31134.80469
[67,] 8261.97404 71262.94012
[68,] -31857.46280 8261.97404
[69,] 87.63759 -31857.46280
[70,] 19349.73122 87.63759
[71,] -21712.39901 19349.73122
[72,] -31931.07118 -21712.39901
[73,] -5759.95982 -31931.07118
[74,] -24455.03002 -5759.95982
[75,] 15178.17166 -24455.03002
[76,] -22027.95419 15178.17166
[77,] -14559.32712 -22027.95419
[78,] -14338.78254 -14559.32712
[79,] -7568.03295 -14338.78254
[80,] 4416.23232 -7568.03295
[81,] -26932.49939 4416.23232
[82,] 89483.47079 -26932.49939
[83,] 15728.32359 89483.47079
[84,] -29572.14820 15728.32359
[85,] -16600.78278 -29572.14820
[86,] -3128.36439 -16600.78278
[87,] -8798.44385 -3128.36439
[88,] 103250.48627 -8798.44385
[89,] -45.26150 103250.48627
[90,] 4954.47517 -45.26150
[91,] -14004.16924 4954.47517
[92,] 32084.42997 -14004.16924
[93,] -16764.27137 32084.42997
[94,] -12639.57546 -16764.27137
[95,] 95111.93807 -12639.57546
[96,] -213.68220 95111.93807
[97,] -11043.61619 -213.68220
[98,] 18505.26173 -11043.61619
[99,] -15016.53557 18505.26173
[100,] -9104.96996 -15016.53557
[101,] 48057.50563 -9104.96996
[102,] 38347.19072 48057.50563
[103,] -7097.98645 38347.19072
[104,] -37160.18704 -7097.98645
[105,] -52775.68033 -37160.18704
[106,] -19799.11663 -52775.68033
[107,] -6864.63780 -19799.11663
[108,] 28181.22642 -6864.63780
[109,] 23542.71336 28181.22642
[110,] -49311.15894 23542.71336
[111,] -9490.07575 -49311.15894
[112,] -28007.37667 -9490.07575
[113,] 10675.43362 -28007.37667
[114,] -42876.57757 10675.43362
[115,] 8135.58822 -42876.57757
[116,] 19103.40790 8135.58822
[117,] 66908.26584 19103.40790
[118,] -123899.03729 66908.26584
[119,] 20606.21775 -123899.03729
[120,] 5917.74806 20606.21775
[121,] 91.89708 5917.74806
[122,] 17865.86842 91.89708
[123,] -45730.92117 17865.86842
[124,] -18577.08243 -45730.92117
[125,] 1349.50929 -18577.08243
[126,] 14684.55211 1349.50929
[127,] -12030.25242 14684.55211
[128,] 27749.91501 -12030.25242
[129,] -227.38513 27749.91501
[130,] 7364.80607 -227.38513
[131,] -22132.36640 7364.80607
[132,] 1527.89985 -22132.36640
[133,] -2990.41814 1527.89985
[134,] -26279.25712 -2990.41814
[135,] -21953.28665 -26279.25712
[136,] -2200.64645 -21953.28665
[137,] 11731.38221 -2200.64645
[138,] -21286.43397 11731.38221
[139,] -13203.74813 -21286.43397
[140,] -16043.89921 -13203.74813
[141,] -6533.17274 -16043.89921
[142,] 29524.35706 -6533.17274
[143,] -4957.13114 29524.35706
[144,] 22604.92743 -4957.13114
[145,] 9667.18999 22604.92743
[146,] 24939.83015 9667.18999
[147,] -28070.11123 24939.83015
[148,] -7580.09059 -28070.11123
[149,] 152.83501 -7580.09059
[150,] 20700.87216 152.83501
[151,] -62249.28626 20700.87216
[152,] 11107.03861 -62249.28626
[153,] 30673.31311 11107.03861
[154,] -63445.71602 30673.31311
[155,] 3654.64495 -63445.71602
[156,] -33426.31991 3654.64495
[157,] 18485.34960 -33426.31991
[158,] -43724.89806 18485.34960
[159,] 64703.31965 -43724.89806
[160,] 51811.28369 64703.31965
[161,] 17734.50395 51811.28369
[162,] 20702.20698 17734.50395
[163,] 108171.26543 20702.20698
[164,] -538.56569 108171.26543
[165,] -27529.59598 -538.56569
[166,] -18315.43309 -27529.59598
[167,] 21479.39061 -18315.43309
[168,] -32037.83494 21479.39061
[169,] 97141.21982 -32037.83494
[170,] 101745.77266 97141.21982
[171,] -16119.00355 101745.77266
[172,] 109510.65655 -16119.00355
[173,] 24445.08328 109510.65655
[174,] -47668.74259 24445.08328
[175,] 28291.57846 -47668.74259
[176,] -10866.26239 28291.57846
[177,] 12404.61751 -10866.26239
[178,] -26811.99110 12404.61751
[179,] 34469.66254 -26811.99110
[180,] 18102.10445 34469.66254
[181,] 15603.38577 18102.10445
[182,] -8623.54272 15603.38577
[183,] 2910.70592 -8623.54272
[184,] -7211.81642 2910.70592
[185,] -13032.79222 -7211.81642
[186,] -2519.86545 -13032.79222
[187,] -7084.27172 -2519.86545
[188,] 5101.59329 -7084.27172
[189,] 8990.97110 5101.59329
[190,] 38800.56704 8990.97110
[191,] 16591.93679 38800.56704
[192,] 8753.14389 16591.93679
[193,] -17579.52992 8753.14389
[194,] -32754.32259 -17579.52992
[195,] 23794.15017 -32754.32259
[196,] -9283.08061 23794.15017
[197,] -14535.45325 -9283.08061
[198,] 990.19278 -14535.45325
[199,] -8739.53924 990.19278
[200,] -10719.32915 -8739.53924
[201,] -25152.60042 -10719.32915
[202,] 48718.82345 -25152.60042
[203,] -555.09565 48718.82345
[204,] 5077.82232 -555.09565
[205,] -23567.12976 5077.82232
[206,] 3578.11226 -23567.12976
[207,] 56522.37831 3578.11226
[208,] 27965.88550 56522.37831
[209,] -169.00990 27965.88550
[210,] 8576.27889 -169.00990
[211,] 31675.48864 8576.27889
[212,] -20075.53116 31675.48864
[213,] -15736.34004 -20075.53116
[214,] -6185.96791 -15736.34004
[215,] -28414.52055 -6185.96791
[216,] -15000.07055 -28414.52055
[217,] -17881.96471 -15000.07055
[218,] -26727.31742 -17881.96471
[219,] -2004.76387 -26727.31742
[220,] -23374.25418 -2004.76387
[221,] -62675.14987 -23374.25418
[222,] -6494.31144 -62675.14987
[223,] -31304.33838 -6494.31144
[224,] -22648.34273 -31304.33838
[225,] 36846.10194 -22648.34273
[226,] 6701.94019 36846.10194
[227,] -22081.72525 6701.94019
[228,] 6023.43656 -22081.72525
[229,] -13196.36940 6023.43656
[230,] -20663.28969 -13196.36940
[231,] -9665.62064 -20663.28969
[232,] -34047.19671 -9665.62064
[233,] -19491.43823 -34047.19671
[234,] -16624.59090 -19491.43823
[235,] -15908.08542 -16624.59090
[236,] 59339.82437 -15908.08542
[237,] -38048.31915 59339.82437
[238,] -11948.39187 -38048.31915
[239,] -22014.39972 -11948.39187
[240,] 1342.22851 -22014.39972
[241,] 11105.93924 1342.22851
[242,] -18859.26737 11105.93924
[243,] -12672.01540 -18859.26737
[244,] -15737.61668 -12672.01540
[245,] -10975.67090 -15737.61668
[246,] 2364.28163 -10975.67090
[247,] -4108.06634 2364.28163
[248,] 25631.80931 -4108.06634
[249,] -15259.77291 25631.80931
[250,] 20517.51812 -15259.77291
[251,] -14215.69160 20517.51812
[252,] -9304.66092 -14215.69160
[253,] 21827.96483 -9304.66092
[254,] 3504.31029 21827.96483
[255,] 3912.77846 3504.31029
[256,] -21169.19638 3912.77846
[257,] 21939.57979 -21169.19638
[258,] -7389.98596 21939.57979
[259,] -18396.83285 -7389.98596
[260,] -19070.50909 -18396.83285
[261,] -17841.05035 -19070.50909
[262,] 8560.72325 -17841.05035
[263,] -17155.91168 8560.72325
[264,] -11974.04042 -17155.91168
[265,] -8553.18336 -11974.04042
[266,] -22320.55457 -8553.18336
[267,] -4804.27121 -22320.55457
[268,] 13906.03012 -4804.27121
[269,] -14580.83776 13906.03012
[270,] 12528.35207 -14580.83776
[271,] -6437.09671 12528.35207
[272,] 27508.13540 -6437.09671
[273,] 6666.65700 27508.13540
[274,] -2201.91950 6666.65700
[275,] 2346.68337 -2201.91950
[276,] 11905.11214 2346.68337
[277,] -10563.76331 11905.11214
[278,] -7444.79525 -10563.76331
[279,] 3181.60823 -7444.79525
[280,] -285.22052 3181.60823
[281,] 36764.94050 -285.22052
[282,] 32411.85739 36764.94050
[283,] 11330.22991 32411.85739
[284,] -197.10752 11330.22991
[285,] -13033.33015 -197.10752
[286,] 6652.88026 -13033.33015
[287,] -28819.43591 6652.88026
[288,] -17239.79177 -28819.43591
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -64615.62528 -36726.19374
2 -14361.63119 -64615.62528
3 -1119.57500 -14361.63119
4 -24162.07284 -1119.57500
5 -57240.15403 -24162.07284
6 90266.66998 -57240.15403
7 -20844.06640 90266.66998
8 39336.30912 -20844.06640
9 -16537.54058 39336.30912
10 2338.83722 -16537.54058
11 -27321.44011 2338.83722
12 -11760.42687 -27321.44011
13 -1386.25821 -11760.42687
14 60853.23138 -1386.25821
15 -20783.97495 60853.23138
16 38976.60406 -20783.97495
17 -408.97338 38976.60406
18 33197.22406 -408.97338
19 -34659.53386 33197.22406
20 -10304.58846 -34659.53386
21 -3451.02317 -10304.58846
22 115993.05436 -3451.02317
23 -22599.70131 115993.05436
24 -50863.58111 -22599.70131
25 -4225.11484 -50863.58111
26 -10578.27125 -4225.11484
27 -7983.46935 -10578.27125
28 -12109.73614 -7983.46935
29 -4722.88606 -12109.73614
30 -79706.52394 -4722.88606
31 -35252.11326 -79706.52394
32 -27372.49982 -35252.11326
33 -32893.14372 -27372.49982
34 -14844.00823 -32893.14372
35 8935.68431 -14844.00823
36 -19838.07024 8935.68431
37 31692.01867 -19838.07024
38 51857.67063 31692.01867
39 -9338.45843 51857.67063
40 -3702.34462 -9338.45843
41 780.41921 -3702.34462
42 16300.44608 780.41921
43 26735.90574 16300.44608
44 52835.83359 26735.90574
45 -11446.49393 52835.83359
46 61835.36613 -11446.49393
47 -33367.99113 61835.36613
48 -7879.48745 -33367.99113
49 -18774.90594 -7879.48745
50 195527.99757 -18774.90594
51 -43287.34856 195527.99757
52 -19359.88142 -43287.34856
53 -7224.63759 -19359.88142
54 -30513.66706 -7224.63759
55 -58444.50224 -30513.66706
56 -14121.07971 -58444.50224
57 50417.38951 -14121.07971
58 -11565.12679 50417.38951
59 4805.63521 -11565.12679
60 30581.14392 4805.63521
61 -39062.23117 30581.14392
62 94822.25479 -39062.23117
63 6112.95992 94822.25479
64 -34281.77825 6112.95992
65 -31134.80469 -34281.77825
66 71262.94012 -31134.80469
67 8261.97404 71262.94012
68 -31857.46280 8261.97404
69 87.63759 -31857.46280
70 19349.73122 87.63759
71 -21712.39901 19349.73122
72 -31931.07118 -21712.39901
73 -5759.95982 -31931.07118
74 -24455.03002 -5759.95982
75 15178.17166 -24455.03002
76 -22027.95419 15178.17166
77 -14559.32712 -22027.95419
78 -14338.78254 -14559.32712
79 -7568.03295 -14338.78254
80 4416.23232 -7568.03295
81 -26932.49939 4416.23232
82 89483.47079 -26932.49939
83 15728.32359 89483.47079
84 -29572.14820 15728.32359
85 -16600.78278 -29572.14820
86 -3128.36439 -16600.78278
87 -8798.44385 -3128.36439
88 103250.48627 -8798.44385
89 -45.26150 103250.48627
90 4954.47517 -45.26150
91 -14004.16924 4954.47517
92 32084.42997 -14004.16924
93 -16764.27137 32084.42997
94 -12639.57546 -16764.27137
95 95111.93807 -12639.57546
96 -213.68220 95111.93807
97 -11043.61619 -213.68220
98 18505.26173 -11043.61619
99 -15016.53557 18505.26173
100 -9104.96996 -15016.53557
101 48057.50563 -9104.96996
102 38347.19072 48057.50563
103 -7097.98645 38347.19072
104 -37160.18704 -7097.98645
105 -52775.68033 -37160.18704
106 -19799.11663 -52775.68033
107 -6864.63780 -19799.11663
108 28181.22642 -6864.63780
109 23542.71336 28181.22642
110 -49311.15894 23542.71336
111 -9490.07575 -49311.15894
112 -28007.37667 -9490.07575
113 10675.43362 -28007.37667
114 -42876.57757 10675.43362
115 8135.58822 -42876.57757
116 19103.40790 8135.58822
117 66908.26584 19103.40790
118 -123899.03729 66908.26584
119 20606.21775 -123899.03729
120 5917.74806 20606.21775
121 91.89708 5917.74806
122 17865.86842 91.89708
123 -45730.92117 17865.86842
124 -18577.08243 -45730.92117
125 1349.50929 -18577.08243
126 14684.55211 1349.50929
127 -12030.25242 14684.55211
128 27749.91501 -12030.25242
129 -227.38513 27749.91501
130 7364.80607 -227.38513
131 -22132.36640 7364.80607
132 1527.89985 -22132.36640
133 -2990.41814 1527.89985
134 -26279.25712 -2990.41814
135 -21953.28665 -26279.25712
136 -2200.64645 -21953.28665
137 11731.38221 -2200.64645
138 -21286.43397 11731.38221
139 -13203.74813 -21286.43397
140 -16043.89921 -13203.74813
141 -6533.17274 -16043.89921
142 29524.35706 -6533.17274
143 -4957.13114 29524.35706
144 22604.92743 -4957.13114
145 9667.18999 22604.92743
146 24939.83015 9667.18999
147 -28070.11123 24939.83015
148 -7580.09059 -28070.11123
149 152.83501 -7580.09059
150 20700.87216 152.83501
151 -62249.28626 20700.87216
152 11107.03861 -62249.28626
153 30673.31311 11107.03861
154 -63445.71602 30673.31311
155 3654.64495 -63445.71602
156 -33426.31991 3654.64495
157 18485.34960 -33426.31991
158 -43724.89806 18485.34960
159 64703.31965 -43724.89806
160 51811.28369 64703.31965
161 17734.50395 51811.28369
162 20702.20698 17734.50395
163 108171.26543 20702.20698
164 -538.56569 108171.26543
165 -27529.59598 -538.56569
166 -18315.43309 -27529.59598
167 21479.39061 -18315.43309
168 -32037.83494 21479.39061
169 97141.21982 -32037.83494
170 101745.77266 97141.21982
171 -16119.00355 101745.77266
172 109510.65655 -16119.00355
173 24445.08328 109510.65655
174 -47668.74259 24445.08328
175 28291.57846 -47668.74259
176 -10866.26239 28291.57846
177 12404.61751 -10866.26239
178 -26811.99110 12404.61751
179 34469.66254 -26811.99110
180 18102.10445 34469.66254
181 15603.38577 18102.10445
182 -8623.54272 15603.38577
183 2910.70592 -8623.54272
184 -7211.81642 2910.70592
185 -13032.79222 -7211.81642
186 -2519.86545 -13032.79222
187 -7084.27172 -2519.86545
188 5101.59329 -7084.27172
189 8990.97110 5101.59329
190 38800.56704 8990.97110
191 16591.93679 38800.56704
192 8753.14389 16591.93679
193 -17579.52992 8753.14389
194 -32754.32259 -17579.52992
195 23794.15017 -32754.32259
196 -9283.08061 23794.15017
197 -14535.45325 -9283.08061
198 990.19278 -14535.45325
199 -8739.53924 990.19278
200 -10719.32915 -8739.53924
201 -25152.60042 -10719.32915
202 48718.82345 -25152.60042
203 -555.09565 48718.82345
204 5077.82232 -555.09565
205 -23567.12976 5077.82232
206 3578.11226 -23567.12976
207 56522.37831 3578.11226
208 27965.88550 56522.37831
209 -169.00990 27965.88550
210 8576.27889 -169.00990
211 31675.48864 8576.27889
212 -20075.53116 31675.48864
213 -15736.34004 -20075.53116
214 -6185.96791 -15736.34004
215 -28414.52055 -6185.96791
216 -15000.07055 -28414.52055
217 -17881.96471 -15000.07055
218 -26727.31742 -17881.96471
219 -2004.76387 -26727.31742
220 -23374.25418 -2004.76387
221 -62675.14987 -23374.25418
222 -6494.31144 -62675.14987
223 -31304.33838 -6494.31144
224 -22648.34273 -31304.33838
225 36846.10194 -22648.34273
226 6701.94019 36846.10194
227 -22081.72525 6701.94019
228 6023.43656 -22081.72525
229 -13196.36940 6023.43656
230 -20663.28969 -13196.36940
231 -9665.62064 -20663.28969
232 -34047.19671 -9665.62064
233 -19491.43823 -34047.19671
234 -16624.59090 -19491.43823
235 -15908.08542 -16624.59090
236 59339.82437 -15908.08542
237 -38048.31915 59339.82437
238 -11948.39187 -38048.31915
239 -22014.39972 -11948.39187
240 1342.22851 -22014.39972
241 11105.93924 1342.22851
242 -18859.26737 11105.93924
243 -12672.01540 -18859.26737
244 -15737.61668 -12672.01540
245 -10975.67090 -15737.61668
246 2364.28163 -10975.67090
247 -4108.06634 2364.28163
248 25631.80931 -4108.06634
249 -15259.77291 25631.80931
250 20517.51812 -15259.77291
251 -14215.69160 20517.51812
252 -9304.66092 -14215.69160
253 21827.96483 -9304.66092
254 3504.31029 21827.96483
255 3912.77846 3504.31029
256 -21169.19638 3912.77846
257 21939.57979 -21169.19638
258 -7389.98596 21939.57979
259 -18396.83285 -7389.98596
260 -19070.50909 -18396.83285
261 -17841.05035 -19070.50909
262 8560.72325 -17841.05035
263 -17155.91168 8560.72325
264 -11974.04042 -17155.91168
265 -8553.18336 -11974.04042
266 -22320.55457 -8553.18336
267 -4804.27121 -22320.55457
268 13906.03012 -4804.27121
269 -14580.83776 13906.03012
270 12528.35207 -14580.83776
271 -6437.09671 12528.35207
272 27508.13540 -6437.09671
273 6666.65700 27508.13540
274 -2201.91950 6666.65700
275 2346.68337 -2201.91950
276 11905.11214 2346.68337
277 -10563.76331 11905.11214
278 -7444.79525 -10563.76331
279 3181.60823 -7444.79525
280 -285.22052 3181.60823
281 36764.94050 -285.22052
282 32411.85739 36764.94050
283 11330.22991 32411.85739
284 -197.10752 11330.22991
285 -13033.33015 -197.10752
286 6652.88026 -13033.33015
287 -28819.43591 6652.88026
288 -17239.79177 -28819.43591
> 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/7udhe1324143944.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/8vipx1324143944.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/9osg71324143944.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/10geor1324143944.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/1148w21324143944.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/12mek41324143944.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/13k8ay1324143944.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/14d8z81324143944.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/15lhmh1324143944.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/16zceo1324143944.tab")
+ }
>
> try(system("convert tmp/1162j1324143944.ps tmp/1162j1324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ioyg1324143944.ps tmp/2ioyg1324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/3aqj51324143944.ps tmp/3aqj51324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/4xeii1324143944.ps tmp/4xeii1324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/5h5n61324143944.ps tmp/5h5n61324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/6d4x11324143944.ps tmp/6d4x11324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/7udhe1324143944.ps tmp/7udhe1324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/8vipx1324143944.ps tmp/8vipx1324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/9osg71324143944.ps tmp/9osg71324143944.png",intern=TRUE))
character(0)
> try(system("convert tmp/10geor1324143944.ps tmp/10geor1324143944.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
7.933 0.612 8.617