R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
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(0
+ ,162556
+ ,1081
+ ,807
+ ,213118
+ ,6282154
+ ,0
+ ,29790
+ ,309
+ ,444
+ ,81767
+ ,4321023
+ ,0
+ ,87550
+ ,458
+ ,412
+ ,153198
+ ,4111912
+ ,1
+ ,84738
+ ,588
+ ,428
+ ,-26007
+ ,223193
+ ,0
+ ,54660
+ ,302
+ ,315
+ ,126942
+ ,1491348
+ ,0
+ ,42634
+ ,156
+ ,168
+ ,157214
+ ,1629616
+ ,1
+ ,40949
+ ,481
+ ,263
+ ,129352
+ ,1398893
+ ,0
+ ,45187
+ ,353
+ ,267
+ ,234817
+ ,1926517
+ ,0
+ ,37704
+ ,452
+ ,228
+ ,60448
+ ,983660
+ ,0
+ ,16275
+ ,109
+ ,129
+ ,47818
+ ,1443586
+ ,1
+ ,25830
+ ,115
+ ,104
+ ,245546
+ ,1073089
+ ,1
+ ,12679
+ ,110
+ ,122
+ ,48020
+ ,984885
+ ,0
+ ,18014
+ ,239
+ ,393
+ ,-1710
+ ,1405225
+ ,1
+ ,43556
+ ,247
+ ,190
+ ,32648
+ ,227132
+ ,0
+ ,24811
+ ,505
+ ,280
+ ,95350
+ ,929118
+ ,1
+ ,6575
+ ,159
+ ,63
+ ,151352
+ ,1071292
+ ,1
+ ,7123
+ ,109
+ ,102
+ ,288170
+ ,638830
+ ,0
+ ,21950
+ ,519
+ ,265
+ ,114337
+ ,856956
+ ,0
+ ,37597
+ ,248
+ ,234
+ ,37884
+ ,992426
+ ,1
+ ,17821
+ ,373
+ ,277
+ ,122844
+ ,444477
+ ,0
+ ,12988
+ ,119
+ ,73
+ ,82340
+ ,857217
+ ,0
+ ,22330
+ ,84
+ ,67
+ ,79801
+ ,711969
+ ,1
+ ,13326
+ ,102
+ ,103
+ ,165548
+ ,702380
+ ,1
+ ,16189
+ ,295
+ ,290
+ ,116384
+ ,358589
+ ,1
+ ,7146
+ ,105
+ ,83
+ ,134028
+ ,297978
+ ,1
+ ,15824
+ ,64
+ ,56
+ ,63838
+ ,585715
+ ,0
+ ,27664
+ ,282
+ ,236
+ ,74996
+ ,657954
+ ,1
+ ,11920
+ ,182
+ ,73
+ ,31080
+ ,209458
+ ,1
+ ,8568
+ ,37
+ ,34
+ ,32168
+ ,786690
+ ,1
+ ,14416
+ ,361
+ ,139
+ ,49857
+ ,439798
+ ,0
+ ,3369
+ ,28
+ ,26
+ ,87161
+ ,688779
+ ,0
+ ,11819
+ ,85
+ ,70
+ ,106113
+ ,574339
+ ,0
+ ,6984
+ ,45
+ ,40
+ ,80570
+ ,741409
+ ,0
+ ,4519
+ ,49
+ ,42
+ ,102129
+ ,597793
+ ,1
+ ,2220
+ ,22
+ ,12
+ ,301670
+ ,644190
+ ,1
+ ,18562
+ ,155
+ ,211
+ ,102313
+ ,377934
+ ,1
+ ,10327
+ ,91
+ ,74
+ ,88577
+ ,640273
+ ,0
+ ,5336
+ ,81
+ ,80
+ ,112477
+ ,697458
+ ,0
+ ,2365
+ ,79
+ ,83
+ ,191778
+ ,550608
+ ,1
+ ,4069
+ ,145
+ ,131
+ ,79804
+ ,207393
+ ,1
+ ,8636
+ ,855
+ ,203
+ ,128294
+ ,301607
+ ,1
+ ,13718
+ ,61
+ ,56
+ ,96448
+ ,345783
+ ,1
+ ,4525
+ ,226
+ ,89
+ ,93811
+ ,501749
+ ,1
+ ,6869
+ ,105
+ ,88
+ ,117520
+ ,379983
+ ,1
+ ,4628
+ ,62
+ ,39
+ ,69159
+ ,387475
+ ,0
+ ,3689
+ ,25
+ ,25
+ ,101792
+ ,377305
+ ,0
+ ,4891
+ ,217
+ ,49
+ ,210568
+ ,370837
+ ,0
+ ,7489
+ ,322
+ ,149
+ ,136996
+ ,430866
+ ,1
+ ,4901
+ ,84
+ ,58
+ ,121920
+ ,469107
+ ,1
+ ,2284
+ ,33
+ ,41
+ ,76403
+ ,194493
+ ,0
+ ,3160
+ ,108
+ ,90
+ ,108094
+ ,530670
+ ,0
+ ,4150
+ ,150
+ ,136
+ ,134759
+ ,518365
+ ,0
+ ,7285
+ ,115
+ ,97
+ ,188873
+ ,491303
+ ,0
+ ,1134
+ ,162
+ ,63
+ ,146216
+ ,527021
+ ,0
+ ,4658
+ ,158
+ ,114
+ ,156608
+ ,233773
+ ,1
+ ,2384
+ ,97
+ ,77
+ ,61348
+ ,405972
+ ,1
+ ,3748
+ ,9
+ ,6
+ ,50350
+ ,652925
+ ,1
+ ,5371
+ ,66
+ ,47
+ ,87720
+ ,446211
+ ,1
+ ,1285
+ ,107
+ ,51
+ ,99489
+ ,341340
+ ,0
+ ,9327
+ ,101
+ ,85
+ ,87419
+ ,387699
+ ,0
+ ,5565
+ ,47
+ ,43
+ ,94355
+ ,493408
+ ,1
+ ,1528
+ ,38
+ ,32
+ ,60326
+ ,146494
+ ,0
+ ,3122
+ ,34
+ ,25
+ ,94670
+ ,414462
+ ,0
+ ,7561
+ ,87
+ ,77
+ ,82425
+ ,364304
+ ,1
+ ,2675
+ ,79
+ ,54
+ ,59017
+ ,355178
+ ,1
+ ,13253
+ ,947
+ ,251
+ ,90829
+ ,357760
+ ,1
+ ,880
+ ,74
+ ,15
+ ,80791
+ ,261216
+ ,0
+ ,2053
+ ,53
+ ,44
+ ,100423
+ ,397144
+ ,1
+ ,1424
+ ,94
+ ,73
+ ,131116
+ ,374943
+ ,0
+ ,4036
+ ,63
+ ,85
+ ,100269
+ ,424898
+ ,0
+ ,3045
+ ,58
+ ,49
+ ,27330
+ ,202055
+ ,1
+ ,5119
+ ,49
+ ,38
+ ,39039
+ ,378525
+ ,1
+ ,1431
+ ,34
+ ,35
+ ,106885
+ ,310768
+ ,1
+ ,554
+ ,11
+ ,9
+ ,79285
+ ,325738
+ ,1
+ ,1975
+ ,35
+ ,34
+ ,118881
+ ,394510
+ ,0
+ ,1765
+ ,20
+ ,20
+ ,77623
+ ,247060
+ ,1
+ ,1012
+ ,47
+ ,29
+ ,114768
+ ,368078
+ ,1
+ ,810
+ ,43
+ ,11
+ ,74015
+ ,236761
+ ,1
+ ,1280
+ ,117
+ ,52
+ ,69465
+ ,312378
+ ,0
+ ,666
+ ,171
+ ,13
+ ,117869
+ ,339836
+ ,1
+ ,1380
+ ,26
+ ,29
+ ,60982
+ ,347385
+ ,0
+ ,4677
+ ,75
+ ,66
+ ,90131
+ ,426280
+ ,1
+ ,876
+ ,59
+ ,33
+ ,138971
+ ,352850
+ ,1
+ ,814
+ ,18
+ ,15
+ ,39625
+ ,301881
+ ,1
+ ,514
+ ,15
+ ,15
+ ,102725
+ ,377516
+ ,0
+ ,5692
+ ,72
+ ,68
+ ,64239
+ ,357312
+ ,1
+ ,3642
+ ,86
+ ,100
+ ,90262
+ ,458343
+ ,1
+ ,540
+ ,14
+ ,13
+ ,103960
+ ,354228
+ ,1
+ ,2099
+ ,64
+ ,45
+ ,106611
+ ,308636
+ ,1
+ ,567
+ ,11
+ ,14
+ ,103345
+ ,386212
+ ,1
+ ,2001
+ ,52
+ ,36
+ ,95551
+ ,393343
+ ,0
+ ,2949
+ ,41
+ ,40
+ ,82903
+ ,378509
+ ,1
+ ,2253
+ ,99
+ ,68
+ ,63593
+ ,452469
+ ,0
+ ,6533
+ ,75
+ ,29
+ ,126910
+ ,364839
+ ,1
+ ,1889
+ ,45
+ ,43
+ ,37527
+ ,358649
+ ,0
+ ,3055
+ ,43
+ ,30
+ ,60247
+ ,376641
+ ,1
+ ,272
+ ,8
+ ,9
+ ,112995
+ ,429112
+ ,0
+ ,1414
+ ,198
+ ,22
+ ,70184
+ ,330546
+ ,1
+ ,2564
+ ,22
+ ,19
+ ,130140
+ ,403560
+ ,0
+ ,1383
+ ,11
+ ,9
+ ,73221
+ ,317892
+ ,0
+ ,1261
+ ,33
+ ,31
+ ,76114
+ ,307528
+ ,1
+ ,975
+ ,23
+ ,19
+ ,90534
+ ,235133
+ ,1
+ ,3366
+ ,80
+ ,55
+ ,108479
+ ,299243
+ ,1
+ ,576
+ ,18
+ ,8
+ ,113761
+ ,314073
+ ,1
+ ,1686
+ ,40
+ ,28
+ ,68696
+ ,368186
+ ,1
+ ,746
+ ,23
+ ,29
+ ,71561
+ ,269661
+ ,0
+ ,3192
+ ,60
+ ,48
+ ,59831
+ ,125390
+ ,0
+ ,2045
+ ,20
+ ,16
+ ,97890
+ ,510834
+ ,1
+ ,5702
+ ,61
+ ,47
+ ,101481
+ ,321896
+ ,0
+ ,1932
+ ,36
+ ,20
+ ,72954
+ ,249898
+ ,1
+ ,936
+ ,30
+ ,22
+ ,67939
+ ,408881
+ ,0
+ ,3437
+ ,47
+ ,33
+ ,48022
+ ,158492
+ ,1
+ ,5131
+ ,71
+ ,44
+ ,86111
+ ,292154
+ ,0
+ ,2397
+ ,14
+ ,13
+ ,74020
+ ,289513
+ ,0
+ ,1389
+ ,9
+ ,6
+ ,57530
+ ,378049
+ ,1
+ ,1503
+ ,39
+ ,35
+ ,56364
+ ,343466
+ ,1
+ ,402
+ ,26
+ ,8
+ ,84990
+ ,332743
+ ,1
+ ,2239
+ ,21
+ ,17
+ ,88590
+ ,442882
+ ,0
+ ,2234
+ ,16
+ ,11
+ ,77200
+ ,214215
+ ,1
+ ,837
+ ,69
+ ,21
+ ,61262
+ ,315688
+ ,1
+ ,10579
+ ,92
+ ,92
+ ,110309
+ ,375195
+ ,1
+ ,875
+ ,14
+ ,12
+ ,67000
+ ,334280
+ ,1
+ ,1585
+ ,107
+ ,112
+ ,93099
+ ,355864
+ ,0
+ ,1659
+ ,29
+ ,25
+ ,107577
+ ,480382
+ ,0
+ ,2647
+ ,37
+ ,17
+ ,62920
+ ,353058
+ ,0
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,217193
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,94
+ ,7
+ ,10
+ ,60793
+ ,314533
+ ,1
+ ,422
+ ,28
+ ,23
+ ,57935
+ ,318056
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,34
+ ,8
+ ,7
+ ,60630
+ ,314353
+ ,1
+ ,1558
+ ,63
+ ,25
+ ,55637
+ ,369448
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,43
+ ,3
+ ,20
+ ,60887
+ ,312846
+ ,1
+ ,645
+ ,5
+ ,4
+ ,60720
+ ,312075
+ ,1
+ ,316
+ ,9
+ ,4
+ ,60505
+ ,315009
+ ,1
+ ,115
+ ,13
+ ,10
+ ,60945
+ ,318903
+ ,1
+ ,5
+ ,2
+ ,1
+ ,60720
+ ,314887
+ ,1
+ ,897
+ ,5
+ ,4
+ ,60720
+ ,314913
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,389
+ ,14
+ ,8
+ ,58990
+ ,325506
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1002
+ ,15
+ ,11
+ ,56750
+ ,298568
+ ,1
+ ,36
+ ,3
+ ,4
+ ,60894
+ ,315834
+ ,1
+ ,460
+ ,15
+ ,15
+ ,63346
+ ,329784
+ ,1
+ ,309
+ ,11
+ ,9
+ ,56535
+ ,312878
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,9
+ ,6
+ ,7
+ ,60835
+ ,314987
+ ,1
+ ,271
+ ,2
+ ,2
+ ,60720
+ ,325249
+ ,1
+ ,14
+ ,1
+ ,0
+ ,61016
+ ,315877
+ ,1
+ ,520
+ ,10
+ ,7
+ ,58650
+ ,291650
+ ,1
+ ,1766
+ ,73
+ ,46
+ ,60438
+ ,305959
+ ,1
+ ,0
+ ,0
+ ,5
+ ,60720
+ ,315380
+ ,1
+ ,458
+ ,11
+ ,7
+ ,58625
+ ,297765
+ ,1
+ ,20
+ ,3
+ ,2
+ ,60938
+ ,315245
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,98
+ ,2
+ ,2
+ ,61490
+ ,315236
+ ,1
+ ,405
+ ,7
+ ,5
+ ,60845
+ ,336425
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,483
+ ,27
+ ,7
+ ,60830
+ ,306268
+ ,1
+ ,454
+ ,51
+ ,24
+ ,63261
+ ,302187
+ ,1
+ ,47
+ ,3
+ ,1
+ ,60720
+ ,314882
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,757
+ ,19
+ ,18
+ ,45689
+ ,382712
+ ,1
+ ,4655
+ ,393
+ ,55
+ ,60720
+ ,341570
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,36
+ ,4
+ ,3
+ ,61564
+ ,312412
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,9
+ ,9
+ ,61938
+ ,309596
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,126
+ ,10
+ ,8
+ ,60951
+ ,315547
+ ,1
+ ,400
+ ,152
+ ,113
+ ,60720
+ ,313267
+ ,0
+ ,71
+ ,1
+ ,0
+ ,60745
+ ,316176
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,972
+ ,34
+ ,19
+ ,71642
+ ,359335
+ ,0
+ ,531
+ ,10
+ ,11
+ ,71641
+ ,330068
+ ,1
+ ,2461
+ ,57
+ ,25
+ ,55792
+ ,314289
+ ,0
+ ,378
+ ,52
+ ,16
+ ,71873
+ ,297413
+ ,0
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,314806
+ ,0
+ ,638
+ ,14
+ ,11
+ ,60370
+ ,333210
+ ,0
+ ,2300
+ ,29
+ ,23
+ ,64873
+ ,352108
+ ,1
+ ,149
+ ,5
+ ,6
+ ,62041
+ ,313332
+ ,1
+ ,226
+ ,5
+ ,5
+ ,65745
+ ,291787
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,275
+ ,4
+ ,7
+ ,59500
+ ,318745
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,141
+ ,6
+ ,7
+ ,61630
+ ,315366
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,28
+ ,2
+ ,3
+ ,60890
+ ,315688
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,4980
+ ,91
+ ,89
+ ,113521
+ ,409642
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,472
+ ,20
+ ,19
+ ,80045
+ ,269587
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,203
+ ,27
+ ,12
+ ,50804
+ ,300962
+ ,0
+ ,496
+ ,17
+ ,12
+ ,87390
+ ,325479
+ ,0
+ ,10
+ ,2
+ ,5
+ ,61656
+ ,316155
+ ,0
+ ,63
+ ,4
+ ,2
+ ,65688
+ ,318574
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1136
+ ,32
+ ,26
+ ,48522
+ ,343613
+ ,0
+ ,265
+ ,31
+ ,3
+ ,60720
+ ,306948
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,267
+ ,32
+ ,11
+ ,57640
+ ,330059
+ ,1
+ ,474
+ ,20
+ ,10
+ ,61977
+ ,288985
+ ,1
+ ,534
+ ,7
+ ,5
+ ,62620
+ ,304485
+ ,1
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,1
+ ,15
+ ,8
+ ,6
+ ,60831
+ ,315688
+ ,1
+ ,397
+ ,28
+ ,7
+ ,60646
+ ,317736
+ ,1
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,1
+ ,1866
+ ,29
+ ,28
+ ,56225
+ ,322331
+ ,1
+ ,288
+ ,4
+ ,3
+ ,60510
+ ,296656
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,3
+ ,2
+ ,1
+ ,60698
+ ,315354
+ ,1
+ ,468
+ ,21
+ ,20
+ ,60720
+ ,312161
+ ,1
+ ,20
+ ,2
+ ,1
+ ,60805
+ ,315576
+ ,1
+ ,278
+ ,26
+ ,22
+ ,61404
+ ,314922
+ ,1
+ ,61
+ ,14
+ ,9
+ ,60720
+ ,314551
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,192
+ ,4
+ ,2
+ ,65276
+ ,312339
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,317
+ ,9
+ ,7
+ ,63915
+ ,298700
+ ,1
+ ,738
+ ,10
+ ,9
+ ,60720
+ ,321376
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,368
+ ,17
+ ,13
+ ,61686
+ ,303230
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,2
+ ,1
+ ,0
+ ,60743
+ ,315487
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,53
+ ,6
+ ,6
+ ,60349
+ ,315793
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,94
+ ,3
+ ,3
+ ,61360
+ ,312887
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,24
+ ,8
+ ,7
+ ,59818
+ ,315637
+ ,1
+ ,2332
+ ,4
+ ,2
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,131
+ ,11
+ ,15
+ ,61808
+ ,308989
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,206
+ ,9
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,167
+ ,2
+ ,1
+ ,64245
+ ,307322
+ ,1
+ ,622
+ ,73
+ ,38
+ ,73007
+ ,304376
+ ,1
+ ,2328
+ ,94
+ ,57
+ ,82732
+ ,253588
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,8
+ ,7
+ ,54820
+ ,309560
+ ,1
+ ,364
+ ,35
+ ,26
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,226
+ ,12
+ ,13
+ ,72835
+ ,343929
+ ,1
+ ,307
+ ,15
+ ,10
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,188
+ ,11
+ ,9
+ ,77655
+ ,381180
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,138
+ ,6
+ ,26
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,125
+ ,12
+ ,19
+ ,60798
+ ,310201
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,282
+ ,30
+ ,12
+ ,62452
+ ,320016
+ ,1
+ ,335
+ ,33
+ ,23
+ ,64175
+ ,320398
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1324
+ ,117
+ ,29
+ ,67440
+ ,291841
+ ,1
+ ,176
+ ,28
+ ,8
+ ,68136
+ ,310670
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,249
+ ,72
+ ,26
+ ,56726
+ ,313491
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,333
+ ,13
+ ,9
+ ,70811
+ ,331323
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,601
+ ,6
+ ,5
+ ,60720
+ ,319210
+ ,1
+ ,30
+ ,4
+ ,3
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,249
+ ,62
+ ,13
+ ,54323
+ ,292754
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,165
+ ,24
+ ,12
+ ,62841
+ ,325176
+ ,1
+ ,453
+ ,21
+ ,19
+ ,81125
+ ,365959
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,53
+ ,14
+ ,10
+ ,59506
+ ,302409
+ ,0
+ ,382
+ ,21
+ ,9
+ ,59365
+ ,340968
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,9
+ ,60720
+ ,315380
+ ,0
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,313164
+ ,1
+ ,290
+ ,2
+ ,1
+ ,58790
+ ,301164
+ ,1
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,366
+ ,53
+ ,14
+ ,61808
+ ,344425
+ ,1
+ ,2
+ ,9
+ ,12
+ ,60735
+ ,315394
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,209
+ ,13
+ ,19
+ ,64016
+ ,316647
+ ,1
+ ,384
+ ,22
+ ,17
+ ,54683
+ ,309836
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,365
+ ,83
+ ,32
+ ,87192
+ ,346611
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,49
+ ,8
+ ,14
+ ,64107
+ ,322031
+ ,1
+ ,3
+ ,4
+ ,8
+ ,60761
+ ,315656
+ ,1
+ ,133
+ ,14
+ ,4
+ ,65990
+ ,339445
+ ,1
+ ,32
+ ,1
+ ,0
+ ,59988
+ ,314964
+ ,1
+ ,368
+ ,17
+ ,20
+ ,61167
+ ,297141
+ ,1
+ ,1
+ ,6
+ ,5
+ ,60719
+ ,315372
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,22
+ ,2
+ ,1
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,96
+ ,5
+ ,4
+ ,60379
+ ,313729
+ ,1
+ ,1
+ ,2
+ ,1
+ ,60727
+ ,315388
+ ,1
+ ,314
+ ,5
+ ,4
+ ,60720
+ ,315371
+ ,1
+ ,844
+ ,78
+ ,20
+ ,60925
+ ,296139
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,26
+ ,1
+ ,1
+ ,60896
+ ,313880
+ ,1
+ ,125
+ ,13
+ ,10
+ ,59734
+ ,317698
+ ,1
+ ,304
+ ,15
+ ,12
+ ,62969
+ ,295580
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,621
+ ,48
+ ,13
+ ,60720
+ ,308256
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,119
+ ,6
+ ,3
+ ,59118
+ ,303677
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1595
+ ,17
+ ,10
+ ,60720
+ ,319369
+ ,1
+ ,312
+ ,14
+ ,3
+ ,58598
+ ,318690
+ ,1
+ ,60
+ ,10
+ ,7
+ ,61124
+ ,314049
+ ,1
+ ,587
+ ,12
+ ,10
+ ,59595
+ ,325699
+ ,1
+ ,135
+ ,2
+ ,1
+ ,62065
+ ,314210
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,514
+ ,52
+ ,15
+ ,78780
+ ,322378
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,1
+ ,4
+ ,4
+ ,60722
+ ,315398
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,1763
+ ,24
+ ,28
+ ,61600
+ ,308336
+ ,1
+ ,180
+ ,11
+ ,9
+ ,59635
+ ,316386
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,218
+ ,21
+ ,7
+ ,60720
+ ,315553
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,448
+ ,40
+ ,7
+ ,59781
+ ,323361
+ ,1
+ ,227
+ ,9
+ ,7
+ ,76644
+ ,336639
+ ,1
+ ,174
+ ,1
+ ,3
+ ,64820
+ ,307424
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,121
+ ,24
+ ,11
+ ,56178
+ ,295370
+ ,1
+ ,607
+ ,11
+ ,7
+ ,60436
+ ,322340
+ ,1
+ ,2212
+ ,14
+ ,10
+ ,60720
+ ,319864
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,530
+ ,60
+ ,18
+ ,73433
+ ,317291
+ ,1
+ ,571
+ ,80
+ ,14
+ ,41477
+ ,280398
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,78
+ ,16
+ ,12
+ ,62700
+ ,317330
+ ,1
+ ,2489
+ ,40
+ ,29
+ ,67804
+ ,238125
+ ,1
+ ,131
+ ,6
+ ,3
+ ,59661
+ ,327071
+ ,1
+ ,923
+ ,8
+ ,6
+ ,58620
+ ,309038
+ ,1
+ ,72
+ ,3
+ ,3
+ ,60398
+ ,314210
+ ,1
+ ,572
+ ,16
+ ,8
+ ,58580
+ ,307930
+ ,1
+ ,397
+ ,10
+ ,10
+ ,62710
+ ,322327
+ ,1
+ ,450
+ ,8
+ ,6
+ ,59325
+ ,292136
+ ,1
+ ,622
+ ,7
+ ,8
+ ,60950
+ ,263276
+ ,1
+ ,694
+ ,8
+ ,6
+ ,68060
+ ,367655
+ ,0
+ ,3425
+ ,12
+ ,9
+ ,83620
+ ,283910
+ ,1
+ ,562
+ ,13
+ ,8
+ ,58456
+ ,283587
+ ,1
+ ,4917
+ ,42
+ ,26
+ ,52811
+ ,243650
+ ,0
+ ,1442
+ ,118
+ ,239
+ ,121173
+ ,438493
+ ,1
+ ,529
+ ,9
+ ,7
+ ,63870
+ ,296261
+ ,0
+ ,2126
+ ,138
+ ,41
+ ,21001
+ ,230621
+ ,1
+ ,1061
+ ,5
+ ,3
+ ,70415
+ ,304252
+ ,1
+ ,776
+ ,9
+ ,8
+ ,64230
+ ,333505
+ ,1
+ ,611
+ ,8
+ ,6
+ ,59190
+ ,296919
+ ,0
+ ,1526
+ ,25
+ ,21
+ ,69351
+ ,278990
+ ,1
+ ,592
+ ,7
+ ,7
+ ,64270
+ ,276898
+ ,1
+ ,1182
+ ,13
+ ,11
+ ,70694
+ ,327007
+ ,1
+ ,621
+ ,16
+ ,11
+ ,68005
+ ,317046
+ ,1
+ ,989
+ ,11
+ ,12
+ ,58930
+ ,304555
+ ,1
+ ,438
+ ,11
+ ,9
+ ,58320
+ ,298096
+ ,1
+ ,726
+ ,3
+ ,3
+ ,69980
+ ,231861
+ ,1
+ ,1303
+ ,61
+ ,57
+ ,69863
+ ,309422
+ ,0
+ ,7419
+ ,29
+ ,21
+ ,63255
+ ,286963
+ ,0
+ ,1164
+ ,17
+ ,15
+ ,57320
+ ,269753
+ ,0
+ ,3310
+ ,33
+ ,32
+ ,75230
+ ,448243
+ ,1
+ ,1920
+ ,15
+ ,11
+ ,79420
+ ,165404
+ ,1
+ ,965
+ ,3
+ ,2
+ ,73490
+ ,204325
+ ,1
+ ,3256
+ ,66
+ ,23
+ ,35250
+ ,407159
+ ,0
+ ,1135
+ ,17
+ ,20
+ ,62285
+ ,290476
+ ,1
+ ,1270
+ ,26
+ ,24
+ ,69206
+ ,275311
+ ,1
+ ,661
+ ,3
+ ,1
+ ,65920
+ ,246541
+ ,1
+ ,1013
+ ,2
+ ,1
+ ,69770
+ ,253468
+ ,1
+ ,2844
+ ,67
+ ,74
+ ,72683
+ ,240897
+ ,0
+ ,11528
+ ,70
+ ,68
+ ,-14545
+ ,-83265
+ ,1
+ ,6526
+ ,26
+ ,20
+ ,55830
+ ,-42143
+ ,1
+ ,2264
+ ,24
+ ,20
+ ,55174
+ ,272713
+ ,0
+ ,5109
+ ,97
+ ,82
+ ,67038
+ ,215362
+ ,1
+ ,3999
+ ,30
+ ,21
+ ,51252
+ ,42754
+ ,1
+ ,35624
+ ,223
+ ,244
+ ,157278
+ ,306275
+ ,1
+ ,9252
+ ,48
+ ,32
+ ,79510
+ ,253537
+ ,1
+ ,15236
+ ,90
+ ,86
+ ,77440
+ ,372631
+ ,1
+ ,18073
+ ,180
+ ,69
+ ,27284
+ ,-7170)
+ ,dim=c(6
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Orders'
+ ,'Dividends'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','Trades','Orders','Dividends','Wealth'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Group Costs Trades Orders Dividends Wealth
1 0 162556 1081 807 213118 6282154
2 0 29790 309 444 81767 4321023
3 0 87550 458 412 153198 4111912
4 1 84738 588 428 -26007 223193
5 0 54660 302 315 126942 1491348
6 0 42634 156 168 157214 1629616
7 1 40949 481 263 129352 1398893
8 0 45187 353 267 234817 1926517
9 0 37704 452 228 60448 983660
10 0 16275 109 129 47818 1443586
11 1 25830 115 104 245546 1073089
12 1 12679 110 122 48020 984885
13 0 18014 239 393 -1710 1405225
14 1 43556 247 190 32648 227132
15 0 24811 505 280 95350 929118
16 1 6575 159 63 151352 1071292
17 1 7123 109 102 288170 638830
18 0 21950 519 265 114337 856956
19 0 37597 248 234 37884 992426
20 1 17821 373 277 122844 444477
21 0 12988 119 73 82340 857217
22 0 22330 84 67 79801 711969
23 1 13326 102 103 165548 702380
24 1 16189 295 290 116384 358589
25 1 7146 105 83 134028 297978
26 1 15824 64 56 63838 585715
27 0 27664 282 236 74996 657954
28 1 11920 182 73 31080 209458
29 1 8568 37 34 32168 786690
30 1 14416 361 139 49857 439798
31 0 3369 28 26 87161 688779
32 0 11819 85 70 106113 574339
33 0 6984 45 40 80570 741409
34 0 4519 49 42 102129 597793
35 1 2220 22 12 301670 644190
36 1 18562 155 211 102313 377934
37 1 10327 91 74 88577 640273
38 0 5336 81 80 112477 697458
39 0 2365 79 83 191778 550608
40 1 4069 145 131 79804 207393
41 1 8636 855 203 128294 301607
42 1 13718 61 56 96448 345783
43 1 4525 226 89 93811 501749
44 1 6869 105 88 117520 379983
45 1 4628 62 39 69159 387475
46 0 3689 25 25 101792 377305
47 0 4891 217 49 210568 370837
48 0 7489 322 149 136996 430866
49 1 4901 84 58 121920 469107
50 1 2284 33 41 76403 194493
51 0 3160 108 90 108094 530670
52 0 4150 150 136 134759 518365
53 0 7285 115 97 188873 491303
54 0 1134 162 63 146216 527021
55 0 4658 158 114 156608 233773
56 1 2384 97 77 61348 405972
57 1 3748 9 6 50350 652925
58 1 5371 66 47 87720 446211
59 1 1285 107 51 99489 341340
60 0 9327 101 85 87419 387699
61 0 5565 47 43 94355 493408
62 1 1528 38 32 60326 146494
63 0 3122 34 25 94670 414462
64 0 7561 87 77 82425 364304
65 1 2675 79 54 59017 355178
66 1 13253 947 251 90829 357760
67 1 880 74 15 80791 261216
68 0 2053 53 44 100423 397144
69 1 1424 94 73 131116 374943
70 0 4036 63 85 100269 424898
71 0 3045 58 49 27330 202055
72 1 5119 49 38 39039 378525
73 1 1431 34 35 106885 310768
74 1 554 11 9 79285 325738
75 1 1975 35 34 118881 394510
76 0 1765 20 20 77623 247060
77 1 1012 47 29 114768 368078
78 1 810 43 11 74015 236761
79 1 1280 117 52 69465 312378
80 0 666 171 13 117869 339836
81 1 1380 26 29 60982 347385
82 0 4677 75 66 90131 426280
83 1 876 59 33 138971 352850
84 1 814 18 15 39625 301881
85 1 514 15 15 102725 377516
86 0 5692 72 68 64239 357312
87 1 3642 86 100 90262 458343
88 1 540 14 13 103960 354228
89 1 2099 64 45 106611 308636
90 1 567 11 14 103345 386212
91 1 2001 52 36 95551 393343
92 0 2949 41 40 82903 378509
93 1 2253 99 68 63593 452469
94 0 6533 75 29 126910 364839
95 1 1889 45 43 37527 358649
96 0 3055 43 30 60247 376641
97 1 272 8 9 112995 429112
98 0 1414 198 22 70184 330546
99 1 2564 22 19 130140 403560
100 0 1383 11 9 73221 317892
101 0 1261 33 31 76114 307528
102 1 975 23 19 90534 235133
103 1 3366 80 55 108479 299243
104 1 576 18 8 113761 314073
105 1 1686 40 28 68696 368186
106 1 746 23 29 71561 269661
107 0 3192 60 48 59831 125390
108 0 2045 20 16 97890 510834
109 1 5702 61 47 101481 321896
110 0 1932 36 20 72954 249898
111 1 936 30 22 67939 408881
112 0 3437 47 33 48022 158492
113 1 5131 71 44 86111 292154
114 0 2397 14 13 74020 289513
115 0 1389 9 6 57530 378049
116 1 1503 39 35 56364 343466
117 1 402 26 8 84990 332743
118 1 2239 21 17 88590 442882
119 0 2234 16 11 77200 214215
120 1 837 69 21 61262 315688
121 1 10579 92 92 110309 375195
122 1 875 14 12 67000 334280
123 1 1585 107 112 93099 355864
124 0 1659 29 25 107577 480382
125 0 2647 37 17 62920 353058
126 0 3294 23 23 75832 217193
127 1 0 0 0 60720 315380
128 1 94 7 10 60793 314533
129 1 422 28 23 57935 318056
130 1 0 0 0 60720 315380
131 1 34 8 7 60630 314353
132 1 1558 63 25 55637 369448
133 1 0 0 0 60720 315380
134 1 43 3 20 60887 312846
135 1 645 5 4 60720 312075
136 1 316 9 4 60505 315009
137 1 115 13 10 60945 318903
138 1 5 2 1 60720 314887
139 1 897 5 4 60720 314913
140 1 0 0 0 60720 315380
141 1 389 14 8 58990 325506
142 1 0 0 0 60720 315380
143 1 1002 15 11 56750 298568
144 1 36 3 4 60894 315834
145 1 460 15 15 63346 329784
146 1 309 11 9 56535 312878
147 1 0 0 0 60720 315380
148 1 9 6 7 60835 314987
149 1 271 2 2 60720 325249
150 1 14 1 0 61016 315877
151 1 520 10 7 58650 291650
152 1 1766 73 46 60438 305959
153 1 0 0 5 60720 315380
154 1 458 11 7 58625 297765
155 1 20 3 2 60938 315245
156 1 0 0 0 60720 315380
157 1 0 0 0 60720 315380
158 1 98 2 2 61490 315236
159 1 405 7 5 60845 336425
160 1 0 0 0 60720 315380
161 1 0 0 0 60720 315380
162 1 0 0 0 60720 315380
163 1 0 0 0 60720 315380
164 1 483 27 7 60830 306268
165 1 454 51 24 63261 302187
166 1 47 3 1 60720 314882
167 1 0 0 0 60720 315380
168 1 757 19 18 45689 382712
169 1 4655 393 55 60720 341570
170 1 0 0 0 60720 315380
171 1 0 0 0 60720 315380
172 1 36 4 3 61564 312412
173 1 0 0 0 60720 315380
174 1 203 9 9 61938 309596
175 1 0 0 0 60720 315380
176 1 126 10 8 60951 315547
177 1 400 152 113 60720 313267
178 0 71 1 0 60745 316176
179 1 0 0 0 60720 315380
180 1 0 0 0 60720 315380
181 1 972 34 19 71642 359335
182 0 531 10 11 71641 330068
183 1 2461 57 25 55792 314289
184 0 378 52 16 71873 297413
185 0 23 5 5 62555 314806
186 0 638 14 11 60370 333210
187 0 2300 29 23 64873 352108
188 1 149 5 6 62041 313332
189 1 226 5 5 65745 291787
190 1 0 0 0 60720 315380
191 1 275 4 7 59500 318745
192 1 0 0 0 60720 315380
193 1 141 6 7 61630 315366
194 1 0 0 0 60720 315380
195 1 28 2 3 60890 315688
196 0 0 0 0 60720 315380
197 0 4980 91 89 113521 409642
198 0 0 0 0 60720 315380
199 0 0 0 0 60720 315380
200 0 472 20 19 80045 269587
201 0 0 0 0 60720 315380
202 0 0 0 0 60720 315380
203 0 0 0 0 60720 315380
204 0 203 27 12 50804 300962
205 0 496 17 12 87390 325479
206 0 10 2 5 61656 316155
207 0 63 4 2 65688 318574
208 0 0 0 0 60720 315380
209 0 1136 32 26 48522 343613
210 0 265 31 3 60720 306948
211 1 0 0 0 60720 315380
212 1 0 0 0 60720 315380
213 1 267 32 11 57640 330059
214 1 474 20 10 61977 288985
215 1 534 7 5 62620 304485
216 1 0 0 2 60720 315380
217 1 15 8 6 60831 315688
218 1 397 28 7 60646 317736
219 1 0 0 2 60720 315380
220 1 1866 29 28 56225 322331
221 1 288 4 3 60510 296656
222 1 0 0 0 60720 315380
223 1 3 2 1 60698 315354
224 1 468 21 20 60720 312161
225 1 20 2 1 60805 315576
226 1 278 26 22 61404 314922
227 1 61 14 9 60720 314551
228 1 0 0 0 60720 315380
229 1 192 4 2 65276 312339
230 1 0 0 0 60720 315380
231 1 317 9 7 63915 298700
232 1 738 10 9 60720 321376
233 1 0 0 0 60720 315380
234 1 368 17 13 61686 303230
235 1 0 0 0 60720 315380
236 1 2 1 0 60743 315487
237 1 0 0 0 60720 315380
238 1 53 6 6 60349 315793
239 1 0 0 0 60720 315380
240 1 0 0 0 60720 315380
241 1 0 0 0 60720 315380
242 1 94 3 3 61360 312887
243 1 0 0 0 60720 315380
244 1 24 8 7 59818 315637
245 1 2332 4 2 72680 324385
246 0 0 0 0 60720 315380
247 0 0 0 0 60720 315380
248 1 131 11 15 61808 308989
249 0 0 0 0 60720 315380
250 0 0 0 0 60720 315380
251 1 206 9 9 53110 296702
252 0 0 0 0 60720 315380
253 1 167 2 1 64245 307322
254 1 622 73 38 73007 304376
255 1 2328 94 57 82732 253588
256 1 0 0 0 60720 315380
257 1 365 8 7 54820 309560
258 1 364 35 26 47705 298466
259 0 0 0 0 60720 315380
260 1 0 0 0 60720 315380
261 0 0 0 0 60720 315380
262 0 0 0 0 60720 315380
263 1 226 12 13 72835 343929
264 1 307 15 10 58856 331955
265 0 0 0 0 60720 315380
266 1 0 0 0 60720 315380
267 0 0 0 0 60720 315380
268 1 188 11 9 77655 381180
269 0 0 0 0 60720 315380
270 1 138 6 26 69817 331420
271 0 0 0 0 60720 315380
272 0 0 0 0 60720 315380
273 0 0 0 0 60720 315380
274 1 125 12 19 60798 310201
275 1 0 0 0 60720 315380
276 1 282 30 12 62452 320016
277 1 335 33 23 64175 320398
278 1 0 0 0 60720 315380
279 0 1324 117 29 67440 291841
280 1 176 28 8 68136 310670
281 0 0 0 0 60720 315380
282 1 0 0 0 60720 315380
283 1 249 72 26 56726 313491
284 1 0 0 0 60720 315380
285 1 333 13 9 70811 331323
286 1 0 0 0 60720 315380
287 0 601 6 5 60720 319210
288 1 30 4 3 62045 318098
289 0 0 0 0 60720 315380
290 1 249 62 13 54323 292754
291 0 0 0 0 60720 315380
292 1 165 24 12 62841 325176
293 1 453 21 19 81125 365959
294 1 0 0 0 60720 315380
295 1 53 14 10 59506 302409
296 0 382 21 9 59365 340968
297 1 0 0 0 60720 315380
298 0 0 0 0 60720 315380
299 1 0 0 0 60720 315380
300 1 0 0 9 60720 315380
301 0 30 4 4 60798 313164
302 1 290 2 1 58790 301164
303 1 0 0 1 60720 315380
304 0 0 0 0 60720 315380
305 1 366 53 14 61808 344425
306 1 2 9 12 60735 315394
307 0 0 0 0 60720 315380
308 0 209 13 19 64016 316647
309 1 384 22 17 54683 309836
310 0 0 0 0 60720 315380
311 0 0 0 0 60720 315380
312 1 365 83 32 87192 346611
313 0 0 0 0 60720 315380
314 0 49 8 14 64107 322031
315 1 3 4 8 60761 315656
316 1 133 14 4 65990 339445
317 1 32 1 0 59988 314964
318 1 368 17 20 61167 297141
319 1 1 6 5 60719 315372
320 1 0 0 0 60720 315380
321 1 0 0 0 60720 315380
322 1 0 0 0 60720 315380
323 1 0 0 0 60720 315380
324 0 0 0 0 60720 315380
325 1 0 0 0 60720 315380
326 1 22 2 1 60722 312502
327 0 0 0 0 60720 315380
328 0 0 0 0 60720 315380
329 1 0 0 0 60720 315380
330 1 0 0 0 60720 315380
331 1 0 0 0 60720 315380
332 1 0 0 0 60720 315380
333 1 0 0 0 60720 315380
334 1 96 5 4 60379 313729
335 1 1 2 1 60727 315388
336 1 314 5 4 60720 315371
337 1 844 78 20 60925 296139
338 1 0 0 0 60720 315380
339 1 26 1 1 60896 313880
340 1 125 13 10 59734 317698
341 1 304 15 12 62969 295580
342 1 0 0 0 60720 315380
343 1 0 0 0 60720 315380
344 1 0 0 0 60720 315380
345 1 621 48 13 60720 308256
346 1 0 0 0 60720 315380
347 1 119 6 3 59118 303677
348 1 0 0 0 60720 315380
349 1 0 0 0 60720 315380
350 1 1595 17 10 60720 319369
351 1 312 14 3 58598 318690
352 1 60 10 7 61124 314049
353 1 587 12 10 59595 325699
354 1 135 2 1 62065 314210
355 1 0 0 0 60720 315380
356 1 0 0 0 60720 315380
357 1 514 52 15 78780 322378
358 1 0 0 0 60720 315380
359 1 0 0 0 60720 315380
360 1 0 0 0 60720 315380
361 1 1 4 4 60722 315398
362 1 0 0 0 60720 315380
363 1 0 0 0 60720 315380
364 0 1763 24 28 61600 308336
365 1 180 11 9 59635 316386
366 1 0 0 0 60720 315380
367 1 0 0 0 60720 315380
368 1 0 0 0 60720 315380
369 1 0 0 0 60720 315380
370 1 218 21 7 60720 315553
371 1 0 0 0 60720 315380
372 1 448 40 7 59781 323361
373 1 227 9 7 76644 336639
374 1 174 1 3 64820 307424
375 1 0 0 0 60720 315380
376 1 0 0 0 60720 315380
377 1 121 24 11 56178 295370
378 1 607 11 7 60436 322340
379 1 2212 14 10 60720 319864
380 1 0 0 0 60720 315380
381 1 0 0 0 60720 315380
382 1 530 60 18 73433 317291
383 1 571 80 14 41477 280398
384 1 0 0 0 60720 315380
385 1 78 16 12 62700 317330
386 1 2489 40 29 67804 238125
387 1 131 6 3 59661 327071
388 1 923 8 6 58620 309038
389 1 72 3 3 60398 314210
390 1 572 16 8 58580 307930
391 1 397 10 10 62710 322327
392 1 450 8 6 59325 292136
393 1 622 7 8 60950 263276
394 1 694 8 6 68060 367655
395 0 3425 12 9 83620 283910
396 1 562 13 8 58456 283587
397 1 4917 42 26 52811 243650
398 0 1442 118 239 121173 438493
399 1 529 9 7 63870 296261
400 0 2126 138 41 21001 230621
401 1 1061 5 3 70415 304252
402 1 776 9 8 64230 333505
403 1 611 8 6 59190 296919
404 0 1526 25 21 69351 278990
405 1 592 7 7 64270 276898
406 1 1182 13 11 70694 327007
407 1 621 16 11 68005 317046
408 1 989 11 12 58930 304555
409 1 438 11 9 58320 298096
410 1 726 3 3 69980 231861
411 1 1303 61 57 69863 309422
412 0 7419 29 21 63255 286963
413 0 1164 17 15 57320 269753
414 0 3310 33 32 75230 448243
415 1 1920 15 11 79420 165404
416 1 965 3 2 73490 204325
417 1 3256 66 23 35250 407159
418 0 1135 17 20 62285 290476
419 1 1270 26 24 69206 275311
420 1 661 3 1 65920 246541
421 1 1013 2 1 69770 253468
422 1 2844 67 74 72683 240897
423 0 11528 70 68 -14545 -83265
424 1 6526 26 20 55830 -42143
425 1 2264 24 20 55174 272713
426 0 5109 97 82 67038 215362
427 1 3999 30 21 51252 42754
428 1 35624 223 244 157278 306275
429 1 9252 48 32 79510 253537
430 1 15236 90 86 77440 372631
431 1 18073 180 69 27284 -7170
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Trades Orders Dividends Wealth
8.582e-01 3.305e-06 3.963e-04 -1.632e-03 -9.164e-07 -1.184e-07
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-0.8362 -0.5269 0.2348 0.2474 0.5940
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 8.582e-01 5.467e-02 15.698 <2e-16 ***
Costs 3.305e-06 4.166e-06 0.793 0.4281
Trades 3.963e-04 3.840e-04 1.032 0.3026
Orders -1.632e-03 8.091e-04 -2.018 0.0443 *
Dividends -9.164e-07 7.264e-07 -1.262 0.2078
Wealth -1.184e-07 8.658e-08 -1.368 0.1721
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 0.4373 on 425 degrees of freedom
Multiple R-squared: 0.06138, Adjusted R-squared: 0.05033
F-statistic: 5.558 on 5 and 425 DF, p-value: 5.668e-05
> 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.7056360 5.887279e-01 2.943640e-01
[2,] 0.5570343 8.859315e-01 4.429657e-01
[3,] 0.7773859 4.452283e-01 2.226141e-01
[4,] 0.8569903 2.860194e-01 1.430097e-01
[5,] 0.8146467 3.707066e-01 1.853533e-01
[6,] 0.7857493 4.285014e-01 2.142507e-01
[7,] 0.7756332 4.487336e-01 2.243668e-01
[8,] 0.7935156 4.129689e-01 2.064844e-01
[9,] 0.7798065 4.403870e-01 2.201935e-01
[10,] 0.7625889 4.748221e-01 2.374111e-01
[11,] 0.7618083 4.763834e-01 2.381917e-01
[12,] 0.7757835 4.484329e-01 2.242165e-01
[13,] 0.7706985 4.586031e-01 2.293015e-01
[14,] 0.7650809 4.698381e-01 2.349191e-01
[15,] 0.7568272 4.863457e-01 2.431728e-01
[16,] 0.7269607 5.460787e-01 2.730393e-01
[17,] 0.7051935 5.896130e-01 2.948065e-01
[18,] 0.7282006 5.435989e-01 2.717994e-01
[19,] 0.7479358 5.041284e-01 2.520642e-01
[20,] 0.7570631 4.858738e-01 2.429369e-01
[21,] 0.7732055 4.535890e-01 2.267945e-01
[22,] 0.7665527 4.668946e-01 2.334473e-01
[23,] 0.7926680 4.146640e-01 2.073320e-01
[24,] 0.8179267 3.641466e-01 1.820733e-01
[25,] 0.8221741 3.556517e-01 1.778259e-01
[26,] 0.8305459 3.389082e-01 1.694541e-01
[27,] 0.8120179 3.759643e-01 1.879821e-01
[28,] 0.8024680 3.950641e-01 1.975320e-01
[29,] 0.8082406 3.835188e-01 1.917594e-01
[30,] 0.8183567 3.632866e-01 1.816433e-01
[31,] 0.8458854 3.082292e-01 1.541146e-01
[32,] 0.8402686 3.194627e-01 1.597314e-01
[33,] 0.8150962 3.698077e-01 1.849038e-01
[34,] 0.8098223 3.803554e-01 1.901777e-01
[35,] 0.8031964 3.936072e-01 1.968036e-01
[36,] 0.7939970 4.120060e-01 2.060030e-01
[37,] 0.7890946 4.218107e-01 2.109054e-01
[38,] 0.8136760 3.726480e-01 1.863240e-01
[39,] 0.8562821 2.874359e-01 1.437179e-01
[40,] 0.8778355 2.443290e-01 1.221645e-01
[41,] 0.8752199 2.495603e-01 1.247801e-01
[42,] 0.8684346 2.631307e-01 1.315654e-01
[43,] 0.8786289 2.427423e-01 1.213711e-01
[44,] 0.8875958 2.248083e-01 1.124042e-01
[45,] 0.8976334 2.047333e-01 1.023666e-01
[46,] 0.9057327 1.885346e-01 9.426730e-02
[47,] 0.9159305 1.681391e-01 8.406953e-02
[48,] 0.9160847 1.678306e-01 8.391529e-02
[49,] 0.9166360 1.667281e-01 8.336404e-02
[50,] 0.9140689 1.718622e-01 8.593111e-02
[51,] 0.9096541 1.806917e-01 9.034587e-02
[52,] 0.9194115 1.611770e-01 8.058849e-02
[53,] 0.9269588 1.460824e-01 7.304120e-02
[54,] 0.9218430 1.563140e-01 7.815701e-02
[55,] 0.9309053 1.381894e-01 6.909468e-02
[56,] 0.9379902 1.240197e-01 6.200983e-02
[57,] 0.9359926 1.280147e-01 6.400737e-02
[58,] 0.9268462 1.463075e-01 7.315376e-02
[59,] 0.9221140 1.557719e-01 7.788596e-02
[60,] 0.9302615 1.394770e-01 6.973851e-02
[61,] 0.9295750 1.408500e-01 7.042502e-02
[62,] 0.9349175 1.301650e-01 6.508252e-02
[63,] 0.9435288 1.129425e-01 5.647123e-02
[64,] 0.9417001 1.165998e-01 5.829990e-02
[65,] 0.9400332 1.199336e-01 5.996678e-02
[66,] 0.9372235 1.255531e-01 6.277653e-02
[67,] 0.9353104 1.293791e-01 6.468955e-02
[68,] 0.9449893 1.100214e-01 5.501068e-02
[69,] 0.9428550 1.142901e-01 5.714504e-02
[70,] 0.9386629 1.226741e-01 6.133707e-02
[71,] 0.9348331 1.303338e-01 6.516692e-02
[72,] 0.9473498 1.053004e-01 5.265021e-02
[73,] 0.9445708 1.108584e-01 5.542921e-02
[74,] 0.9511335 9.773299e-02 4.886649e-02
[75,] 0.9488627 1.022747e-01 5.113733e-02
[76,] 0.9452036 1.095929e-01 5.479643e-02
[77,] 0.9421035 1.157931e-01 5.789654e-02
[78,] 0.9501882 9.962365e-02 4.981182e-02
[79,] 0.9495057 1.009885e-01 5.049427e-02
[80,] 0.9463901 1.072198e-01 5.360990e-02
[81,] 0.9432187 1.135626e-01 5.678132e-02
[82,] 0.9397085 1.205830e-01 6.029149e-02
[83,] 0.9359874 1.280253e-01 6.401263e-02
[84,] 0.9458296 1.083407e-01 5.417036e-02
[85,] 0.9428504 1.142992e-01 5.714960e-02
[86,] 0.9528291 9.434180e-02 4.717090e-02
[87,] 0.9495013 1.009974e-01 5.049870e-02
[88,] 0.9584596 8.308081e-02 4.154040e-02
[89,] 0.9558095 8.838102e-02 4.419051e-02
[90,] 0.9668554 6.628914e-02 3.314457e-02
[91,] 0.9646211 7.075781e-02 3.537890e-02
[92,] 0.9713155 5.736908e-02 2.868454e-02
[93,] 0.9763222 4.735552e-02 2.367776e-02
[94,] 0.9744470 5.110599e-02 2.555300e-02
[95,] 0.9726583 5.468346e-02 2.734173e-02
[96,] 0.9705642 5.887161e-02 2.943580e-02
[97,] 0.9683735 6.325303e-02 3.162652e-02
[98,] 0.9660287 6.794266e-02 3.397133e-02
[99,] 0.9725127 5.497453e-02 2.748727e-02
[100,] 0.9775028 4.499439e-02 2.249719e-02
[101,] 0.9757839 4.843219e-02 2.421610e-02
[102,] 0.9806601 3.867984e-02 1.933992e-02
[103,] 0.9792965 4.140700e-02 2.070350e-02
[104,] 0.9836417 3.271661e-02 1.635831e-02
[105,] 0.9822635 3.547294e-02 1.773647e-02
[106,] 0.9859131 2.817389e-02 1.408695e-02
[107,] 0.9887382 2.252354e-02 1.126177e-02
[108,] 0.9880046 2.399074e-02 1.199537e-02
[109,] 0.9870156 2.596884e-02 1.298442e-02
[110,] 0.9860997 2.780069e-02 1.390035e-02
[111,] 0.9892762 2.144759e-02 1.072379e-02
[112,] 0.9883002 2.339957e-02 1.169978e-02
[113,] 0.9875898 2.482045e-02 1.241022e-02
[114,] 0.9865752 2.684959e-02 1.342480e-02
[115,] 0.9861530 2.769401e-02 1.384701e-02
[116,] 0.9888309 2.233811e-02 1.116905e-02
[117,] 0.9913585 1.728291e-02 8.641453e-03
[118,] 0.9934247 1.315058e-02 6.575292e-03
[119,] 0.9928675 1.426498e-02 7.132490e-03
[120,] 0.9922768 1.544647e-02 7.723234e-03
[121,] 0.9916538 1.669236e-02 8.346178e-03
[122,] 0.9908722 1.825553e-02 9.127766e-03
[123,] 0.9900241 1.995172e-02 9.975859e-03
[124,] 0.9891036 2.179286e-02 1.089643e-02
[125,] 0.9880163 2.396735e-02 1.198367e-02
[126,] 0.9870066 2.598682e-02 1.299341e-02
[127,] 0.9856802 2.863952e-02 1.431976e-02
[128,] 0.9842024 3.159529e-02 1.579765e-02
[129,] 0.9826364 3.472719e-02 1.736359e-02
[130,] 0.9808241 3.835175e-02 1.917588e-02
[131,] 0.9788184 4.236315e-02 2.118157e-02
[132,] 0.9765940 4.681205e-02 2.340602e-02
[133,] 0.9742173 5.156541e-02 2.578270e-02
[134,] 0.9715162 5.696755e-02 2.848378e-02
[135,] 0.9686201 6.275976e-02 3.137988e-02
[136,] 0.9654401 6.911982e-02 3.455991e-02
[137,] 0.9621755 7.564907e-02 3.782454e-02
[138,] 0.9584715 8.305703e-02 4.152852e-02
[139,] 0.9542510 9.149794e-02 4.574897e-02
[140,] 0.9498352 1.003296e-01 5.016482e-02
[141,] 0.9449117 1.101767e-01 5.508835e-02
[142,] 0.9394781 1.210438e-01 6.052188e-02
[143,] 0.9336520 1.326959e-01 6.634795e-02
[144,] 0.9283468 1.433064e-01 7.165322e-02
[145,] 0.9218494 1.563012e-01 7.815062e-02
[146,] 0.9146696 1.706608e-01 8.533040e-02
[147,] 0.9069227 1.861546e-01 9.307729e-02
[148,] 0.8985545 2.028911e-01 1.014455e-01
[149,] 0.8896031 2.207938e-01 1.103969e-01
[150,] 0.8801303 2.397394e-01 1.198697e-01
[151,] 0.8702814 2.594372e-01 1.297186e-01
[152,] 0.8595085 2.809830e-01 1.404915e-01
[153,] 0.8480976 3.038048e-01 1.519024e-01
[154,] 0.8360418 3.279165e-01 1.639582e-01
[155,] 0.8233373 3.533253e-01 1.766627e-01
[156,] 0.8096827 3.806346e-01 1.903173e-01
[157,] 0.7963394 4.073213e-01 2.036606e-01
[158,] 0.7817054 4.365892e-01 2.182946e-01
[159,] 0.7664471 4.671058e-01 2.335529e-01
[160,] 0.7535768 4.928465e-01 2.464232e-01
[161,] 0.7330994 5.338012e-01 2.669006e-01
[162,] 0.7160969 5.678061e-01 2.839031e-01
[163,] 0.6985608 6.028784e-01 3.014392e-01
[164,] 0.6806758 6.386484e-01 3.193242e-01
[165,] 0.6621755 6.756491e-01 3.378245e-01
[166,] 0.6439955 7.120090e-01 3.560045e-01
[167,] 0.6246993 7.506014e-01 3.753007e-01
[168,] 0.6057609 7.884781e-01 3.942391e-01
[169,] 0.6005112 7.989776e-01 3.994888e-01
[170,] 0.6557533 6.884935e-01 3.442467e-01
[171,] 0.6367193 7.265613e-01 3.632807e-01
[172,] 0.6173215 7.653570e-01 3.826785e-01
[173,] 0.5989703 8.020594e-01 4.010297e-01
[174,] 0.6512192 6.975616e-01 3.487808e-01
[175,] 0.6333421 7.333159e-01 3.666579e-01
[176,] 0.6897730 6.204539e-01 3.102270e-01
[177,] 0.7367337 5.265326e-01 2.632663e-01
[178,] 0.7770063 4.459873e-01 2.229937e-01
[179,] 0.8118354 3.763293e-01 1.881646e-01
[180,] 0.7986999 4.026003e-01 2.013001e-01
[181,] 0.7843672 4.312655e-01 2.156328e-01
[182,] 0.7696028 4.607944e-01 2.303972e-01
[183,] 0.7552404 4.895192e-01 2.447596e-01
[184,] 0.7393904 5.212191e-01 2.606096e-01
[185,] 0.7237353 5.525295e-01 2.762647e-01
[186,] 0.7068700 5.862601e-01 2.931300e-01
[187,] 0.6899230 6.201541e-01 3.100770e-01
[188,] 0.7376297 5.247406e-01 2.623703e-01
[189,] 0.7709097 4.581806e-01 2.290903e-01
[190,] 0.8109040 3.781920e-01 1.890960e-01
[191,] 0.8461403 3.077193e-01 1.538597e-01
[192,] 0.8799962 2.400076e-01 1.200038e-01
[193,] 0.9055043 1.889914e-01 9.449568e-02
[194,] 0.9267838 1.464323e-01 7.321615e-02
[195,] 0.9442236 1.115528e-01 5.577638e-02
[196,] 0.9567797 8.644066e-02 4.322033e-02
[197,] 0.9711563 5.768744e-02 2.884372e-02
[198,] 0.9789801 4.203987e-02 2.101993e-02
[199,] 0.9855775 2.884495e-02 1.442247e-02
[200,] 0.9900882 1.982368e-02 9.911840e-03
[201,] 0.9923015 1.539703e-02 7.698517e-03
[202,] 0.9952107 9.578533e-03 4.789267e-03
[203,] 0.9945227 1.095458e-02 5.477292e-03
[204,] 0.9937452 1.250950e-02 6.254751e-03
[205,] 0.9928846 1.423083e-02 7.115416e-03
[206,] 0.9918733 1.625333e-02 8.126667e-03
[207,] 0.9907478 1.850449e-02 9.252247e-03
[208,] 0.9895315 2.093698e-02 1.046849e-02
[209,] 0.9881845 2.363093e-02 1.181547e-02
[210,] 0.9865499 2.690024e-02 1.345012e-02
[211,] 0.9848561 3.028783e-02 1.514391e-02
[212,] 0.9834955 3.300899e-02 1.650449e-02
[213,] 0.9814097 3.718056e-02 1.859028e-02
[214,] 0.9791090 4.178205e-02 2.089102e-02
[215,] 0.9765720 4.685603e-02 2.342802e-02
[216,] 0.9741837 5.163261e-02 2.581630e-02
[217,] 0.9711429 5.771412e-02 2.885706e-02
[218,] 0.9682820 6.343604e-02 3.171802e-02
[219,] 0.9648127 7.037457e-02 3.518729e-02
[220,] 0.9608698 7.826035e-02 3.913017e-02
[221,] 0.9563286 8.734270e-02 4.367135e-02
[222,] 0.9516239 9.675223e-02 4.837612e-02
[223,] 0.9464043 1.071913e-01 5.359566e-02
[224,] 0.9412176 1.175648e-01 5.878238e-02
[225,] 0.9352570 1.294861e-01 6.474304e-02
[226,] 0.9291269 1.417461e-01 7.087305e-02
[227,] 0.9222368 1.555265e-01 7.776324e-02
[228,] 0.9148010 1.703981e-01 8.519903e-02
[229,] 0.9068881 1.862239e-01 9.311193e-02
[230,] 0.8988915 2.022170e-01 1.011085e-01
[231,] 0.8899279 2.201442e-01 1.100721e-01
[232,] 0.8804072 2.391857e-01 1.195928e-01
[233,] 0.8703201 2.593597e-01 1.296799e-01
[234,] 0.8597607 2.804786e-01 1.402393e-01
[235,] 0.8485277 3.029447e-01 1.514723e-01
[236,] 0.8374700 3.250600e-01 1.625300e-01
[237,] 0.8232790 3.534419e-01 1.767210e-01
[238,] 0.8585539 2.828922e-01 1.414461e-01
[239,] 0.8890873 2.218254e-01 1.109127e-01
[240,] 0.8806207 2.387585e-01 1.193793e-01
[241,] 0.9078983 1.842033e-01 9.210165e-02
[242,] 0.9305964 1.388071e-01 6.940355e-02
[243,] 0.9249311 1.501378e-01 7.506888e-02
[244,] 0.9445165 1.109671e-01 5.548353e-02
[245,] 0.9382846 1.234308e-01 6.171540e-02
[246,] 0.9315174 1.369652e-01 6.848260e-02
[247,] 0.9240259 1.519483e-01 7.597414e-02
[248,] 0.9164570 1.670860e-01 8.354299e-02
[249,] 0.9094831 1.810338e-01 9.051688e-02
[250,] 0.9048766 1.902468e-01 9.512342e-02
[251,] 0.9288095 1.423810e-01 7.119052e-02
[252,] 0.9215512 1.568977e-01 7.844885e-02
[253,] 0.9424524 1.150952e-01 5.754761e-02
[254,] 0.9590349 8.193027e-02 4.096514e-02
[255,] 0.9541315 9.173698e-02 4.586849e-02
[256,] 0.9494871 1.010258e-01 5.051290e-02
[257,] 0.9646320 7.073592e-02 3.536796e-02
[258,] 0.9603074 7.938519e-02 3.969260e-02
[259,] 0.9729209 5.415829e-02 2.707914e-02
[260,] 0.9691546 6.169085e-02 3.084542e-02
[261,] 0.9795438 4.091247e-02 2.045623e-02
[262,] 0.9774016 4.519685e-02 2.259842e-02
[263,] 0.9855036 2.899282e-02 1.449641e-02
[264,] 0.9911146 1.777088e-02 8.885440e-03
[265,] 0.9948224 1.035528e-02 5.177640e-03
[266,] 0.9941376 1.172489e-02 5.862445e-03
[267,] 0.9931359 1.372828e-02 6.864138e-03
[268,] 0.9920229 1.595418e-02 7.977090e-03
[269,] 0.9909245 1.815104e-02 9.075519e-03
[270,] 0.9894510 2.109803e-02 1.054902e-02
[271,] 0.9951806 9.638775e-03 4.819388e-03
[272,] 0.9942177 1.156454e-02 5.782268e-03
[273,] 0.9967782 6.443577e-03 3.221788e-03
[274,] 0.9961847 7.630584e-03 3.815292e-03
[275,] 0.9954830 9.034100e-03 4.517050e-03
[276,] 0.9946779 1.064414e-02 5.322068e-03
[277,] 0.9936685 1.266296e-02 6.331479e-03
[278,] 0.9925820 1.483601e-02 7.418005e-03
[279,] 0.9957550 8.490098e-03 4.245049e-03
[280,] 0.9949797 1.004069e-02 5.020343e-03
[281,] 0.9973349 5.330287e-03 2.665143e-03
[282,] 0.9967665 6.466923e-03 3.233462e-03
[283,] 0.9983837 3.232637e-03 1.616318e-03
[284,] 0.9980506 3.898885e-03 1.949443e-03
[285,] 0.9976179 4.764188e-03 2.382094e-03
[286,] 0.9971322 5.735694e-03 2.867847e-03
[287,] 0.9966098 6.780357e-03 3.390179e-03
[288,] 0.9982545 3.491045e-03 1.745523e-03
[289,] 0.9978808 4.238466e-03 2.119233e-03
[290,] 0.9990263 1.947488e-03 9.737441e-04
[291,] 0.9988013 2.397429e-03 1.198714e-03
[292,] 0.9985720 2.855994e-03 1.427997e-03
[293,] 0.9993598 1.280371e-03 6.401854e-04
[294,] 0.9992060 1.588008e-03 7.940039e-04
[295,] 0.9990179 1.964204e-03 9.821022e-04
[296,] 0.9996040 7.920943e-04 3.960471e-04
[297,] 0.9994951 1.009828e-03 5.049138e-04
[298,] 0.9993849 1.230220e-03 6.151101e-04
[299,] 0.9997725 4.549800e-04 2.274900e-04
[300,] 0.9999081 1.838831e-04 9.194155e-05
[301,] 0.9998878 2.243433e-04 1.121716e-04
[302,] 0.9999679 6.415857e-05 3.207929e-05
[303,] 0.9999924 1.511891e-05 7.559454e-06
[304,] 0.9999895 2.107204e-05 1.053602e-05
[305,] 0.9999980 4.051899e-06 2.025949e-06
[306,] 0.9999996 7.813511e-07 3.906756e-07
[307,] 0.9999994 1.105512e-06 5.527560e-07
[308,] 0.9999992 1.617642e-06 8.088208e-07
[309,] 0.9999988 2.331371e-06 1.165685e-06
[310,] 0.9999984 3.123875e-06 1.561938e-06
[311,] 0.9999978 4.416212e-06 2.208106e-06
[312,] 0.9999969 6.297446e-06 3.148723e-06
[313,] 0.9999955 8.940832e-06 4.470416e-06
[314,] 0.9999937 1.263775e-05 6.318875e-06
[315,] 0.9999911 1.778359e-05 8.891796e-06
[316,] 0.9999986 2.700422e-06 1.350211e-06
[317,] 0.9999980 3.931241e-06 1.965620e-06
[318,] 0.9999972 5.687950e-06 2.843975e-06
[319,] 0.9999997 6.470672e-07 3.235336e-07
[320,] 1.0000000 4.596317e-08 2.298158e-08
[321,] 1.0000000 7.282579e-08 3.641289e-08
[322,] 0.9999999 1.148189e-07 5.740945e-08
[323,] 0.9999999 1.801224e-07 9.006119e-08
[324,] 0.9999999 2.811391e-07 1.405696e-07
[325,] 0.9999998 4.365632e-07 2.182816e-07
[326,] 0.9999997 6.669655e-07 3.334827e-07
[327,] 0.9999995 1.022999e-06 5.114995e-07
[328,] 0.9999992 1.551186e-06 7.755931e-07
[329,] 0.9999988 2.308018e-06 1.154009e-06
[330,] 0.9999983 3.494050e-06 1.747025e-06
[331,] 0.9999974 5.251293e-06 2.625647e-06
[332,] 0.9999962 7.598105e-06 3.799052e-06
[333,] 0.9999945 1.104293e-05 5.521466e-06
[334,] 0.9999918 1.636080e-05 8.180401e-06
[335,] 0.9999879 2.410050e-05 1.205025e-05
[336,] 0.9999824 3.529559e-05 1.764779e-05
[337,] 0.9999747 5.061577e-05 2.530789e-05
[338,] 0.9999633 7.330384e-05 3.665192e-05
[339,] 0.9999478 1.044632e-04 5.223158e-05
[340,] 0.9999253 1.494888e-04 7.474440e-05
[341,] 0.9998937 2.126111e-04 1.063055e-04
[342,] 0.9998515 2.969205e-04 1.484602e-04
[343,] 0.9997929 4.142819e-04 2.071409e-04
[344,] 0.9997152 5.696573e-04 2.848286e-04
[345,] 0.9996156 7.687655e-04 3.843828e-04
[346,] 0.9994685 1.062944e-03 5.314720e-04
[347,] 0.9992717 1.456505e-03 7.282524e-04
[348,] 0.9990087 1.982649e-03 9.913247e-04
[349,] 0.9986377 2.724583e-03 1.362291e-03
[350,] 0.9981697 3.660601e-03 1.830301e-03
[351,] 0.9975576 4.884887e-03 2.442444e-03
[352,] 0.9967630 6.474076e-03 3.237038e-03
[353,] 0.9957771 8.445877e-03 4.222938e-03
[354,] 0.9944798 1.104049e-02 5.520244e-03
[355,] 0.9928346 1.433080e-02 7.165402e-03
[356,] 0.9971114 5.777179e-03 2.888590e-03
[357,] 0.9962287 7.542623e-03 3.771311e-03
[358,] 0.9949957 1.000851e-02 5.004255e-03
[359,] 0.9934099 1.318011e-02 6.590057e-03
[360,] 0.9913879 1.722419e-02 8.612094e-03
[361,] 0.9888323 2.233550e-02 1.116775e-02
[362,] 0.9857224 2.855529e-02 1.427765e-02
[363,] 0.9817657 3.646851e-02 1.823426e-02
[364,] 0.9768955 4.620907e-02 2.310453e-02
[365,] 0.9704872 5.902560e-02 2.951280e-02
[366,] 0.9630149 7.397010e-02 3.698505e-02
[367,] 0.9542157 9.156869e-02 4.578435e-02
[368,] 0.9437836 1.124329e-01 5.621644e-02
[369,] 0.9332733 1.334535e-01 6.672673e-02
[370,] 0.9200214 1.599573e-01 7.997864e-02
[371,] 0.9045822 1.908355e-01 9.541777e-02
[372,] 0.8864679 2.270641e-01 1.135321e-01
[373,] 0.8660338 2.679325e-01 1.339662e-01
[374,] 0.8408458 3.183083e-01 1.591542e-01
[375,] 0.8248154 3.503692e-01 1.751846e-01
[376,] 0.7973036 4.053927e-01 2.026964e-01
[377,] 0.7706138 4.587724e-01 2.293862e-01
[378,] 0.7408099 5.183801e-01 2.591901e-01
[379,] 0.7085939 5.828121e-01 2.914061e-01
[380,] 0.6752046 6.495909e-01 3.247954e-01
[381,] 0.6396707 7.206587e-01 3.603293e-01
[382,] 0.6057325 7.885351e-01 3.942675e-01
[383,] 0.5707726 8.584548e-01 4.292274e-01
[384,] 0.5345039 9.309922e-01 4.654961e-01
[385,] 0.4963204 9.926407e-01 5.036796e-01
[386,] 0.4574908 9.149817e-01 5.425092e-01
[387,] 0.6799367 6.401266e-01 3.200633e-01
[388,] 0.6450563 7.098875e-01 3.549437e-01
[389,] 0.6157480 7.685041e-01 3.842520e-01
[390,] 0.5665791 8.668419e-01 4.334209e-01
[391,] 0.5230604 9.538793e-01 4.769396e-01
[392,] 0.5608803 8.782393e-01 4.391197e-01
[393,] 0.5083019 9.833963e-01 4.916981e-01
[394,] 0.4690639 9.381278e-01 5.309361e-01
[395,] 0.4340420 8.680841e-01 5.659580e-01
[396,] 0.5860884 8.278231e-01 4.139116e-01
[397,] 0.5373235 9.253530e-01 4.626765e-01
[398,] 0.4847728 9.695457e-01 5.152272e-01
[399,] 0.4331903 8.663806e-01 5.668097e-01
[400,] 0.4096764 8.193528e-01 5.903236e-01
[401,] 0.3860032 7.720064e-01 6.139968e-01
[402,] 0.3293756 6.587512e-01 6.706244e-01
[403,] 0.3119451 6.238901e-01 6.880549e-01
[404,] 0.4577870 9.155741e-01 5.422130e-01
[405,] 0.5692825 8.614349e-01 4.307175e-01
[406,] 0.7821115 4.357771e-01 2.178885e-01
[407,] 0.7060240 5.879520e-01 2.939760e-01
[408,] 0.6156869 7.686263e-01 3.843131e-01
[409,] 0.5786610 8.426780e-01 4.213390e-01
[410,] 0.8528958 2.942084e-01 1.471042e-01
[411,] 0.7727166 4.545667e-01 2.272834e-01
[412,] 0.6577749 6.844503e-01 3.422251e-01
[413,] 0.5233202 9.533595e-01 4.766798e-01
[414,] 0.7330536 5.338927e-01 2.669464e-01
> postscript(file="/var/www/html/rcomp/tmp/13tdy1291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/23tdy1291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/33tdy1291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4wkcj1291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5wkcj1291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 431
Frequency = 1
1 2 3 4 5 6
0.43293581 0.23239332 -0.02906483 0.33008759 -0.35131265 -0.44957037
7 8 9 10 11 12
0.52943665 -0.26818786 -0.61779874 -0.52978349 0.53276565 0.41613882
13 14 15 16 17 18
-0.20601138 0.26699770 -0.48578614 0.42550073 0.58132657 -0.49751323
19 20 21 22 23 24
-0.54645068 0.55252200 -0.65210778 -0.69843164 0.46039586 0.59395424
25 26 27 28 29 30
0.37020477 0.28345647 -0.52943794 0.20277225 0.27699963 0.27581767
31 32 33 34 35 36
-0.67651550 -0.65138491 -0.67215069 -0.65957760 0.49809064 0.50202496
37 38 39 40 41 42
0.34943671 -0.59164014 -0.52085581 0.38246263 0.25914808 0.29307248
43 44 45 46 47 48
0.32799116 0.37386669 0.27489395 -0.70149628 -0.64346182 -0.59072091
49 50 51 52 53 54
0.35430531 0.28118013 -0.60259299 -0.52443854 -0.53821088 -0.62687094
55 56 57 58 59 60
-0.57888273 0.32550517 0.25913457 0.30787744 0.31002748 -0.66423740
61 62 63 64 65 66
-0.68009488 0.24658854 -0.70531491 -0.67326022 0.28597908 0.25811008
67 68 69 70 71 72
0.23905231 -0.67507442 0.38359411 -0.61551461 -0.76225453 0.24812997
73 74 75 76 77 78
0.31551042 0.26156090 0.33259396 -0.73889110 0.31595963 0.23593282
79 80 81 82 83 84
0.27677127 -0.75865595 0.27132753 -0.66253050 0.33855867 0.22855453
85 86 87 88 89 90
0.29751384 -0.69332482 0.39595052 0.29293311 0.31723574 0.29888932
91 92 93 94 95 96
0.30751962 -0.69807119 0.31801737 -0.70263724 0.26481137 -0.73652044
97 98 99 100 101 102
0.30681413 -0.80193158 0.32270161 -0.74766448 -0.71864211 0.27131681
103 104 105 106 107 108
0.32363232 0.28729268 0.27266835 0.27510089 -0.74446203 -0.69653640
109 110 111 112 113 114
0.30665270 -0.74972558 0.27344052 -0.77150768 0.28207300 -0.74830300
115 116 117 118 119 120
-0.75904337 0.27086835 0.26054414 0.28748925 -0.75782506 0.23952578
121 122 123 124 125 126
0.36611330 0.25396284 0.40447888 -0.67886500 -0.75435964 -0.74541345
127 128 129 130 131 132
0.23481989 0.24802628 0.25764067 0.23481989 0.24276031 0.24726269
133 134 135 136 137 138
0.23481989 0.26599053 0.23684546 0.23649803 0.24623610 0.23558488
139 140 141 142 143 144
0.23634882 0.23481989 0.24066004 0.23481989 0.23789258 0.24025505
145 146 147 148 149 150
0.25595443 0.24000048 0.23481989 0.24389839 0.23756545 0.23470745
151 152 153 154 155 156
0.23585876 0.27377432 0.24298203 0.23636865 0.23701363 0.23481989
157 158 159 160 161 162
0.23481989 0.23765691 0.24147668 0.23481989 0.23481989 0.23481989
163 164 165 166 167 168
0.23481989 0.23297331 0.25305441 0.23504923 0.23481989 0.24837324
169 170 171 172 173 174
0.15658986 0.23481989 0.23481989 0.23843506 0.23481989 0.24570568
175 176 177 178 179 180
0.23481989 0.24373175 0.35747989 -0.76569382 0.23481989 0.23481989
181 182 183 184 185 186
0.26436488 -0.74119379 0.24026604 -0.75282388 -0.75746177 -0.75308857
187 188 189 190 191 192
-0.73857113 0.24310873 0.24206452 0.23481989 0.24303364 0.23481989
193 194 195 196 197 198
0.24423558 0.23481989 0.23902437 -0.76518011 -0.61286300 -0.76518011
199 200 201 202 203 204
-0.76518011 -0.73136359 -0.76518011 -0.76518011 -0.76518011 -0.76775502
205 206 207 208 209 210
-0.72833137 -0.75689406 -0.75877780 -0.76518011 -0.74700553 -0.77444136
211 212 213 214 215 216
0.23481989 0.23481989 0.23812980 0.23967849 0.23889437 0.23808475
217 218 219 220 221 222
0.24153295 0.23405074 0.23808475 0.25957415 0.23477055 0.23481989
223 224 225 226 227 228
0.23562663 0.25721916 0.23569479 0.26008430 0.24366427 0.23481989
229 230 231 232 233 234
0.23967999 0.23481989 0.24258536 0.24382045 0.23481989 0.24753519
235 236 237 238 239 240
0.23481989 0.23445076 0.23481989 0.24177067 0.23481989 0.23481989
241 242 243 244 245 246
0.23481989 0.23850898 0.23481989 0.24220134 0.24081959 -0.76518011
247 248 249 250 251 252
-0.76518011 0.25475464 -0.76518011 -0.76518011 0.23607919 -0.76518011
253 254 255 256 257 258
0.23738379 0.27582543 0.29577938 0.23481989 0.23577488 0.24826147
259 260 261 262 263 264
-0.76518011 0.23481989 -0.76518011 -0.76518011 0.26502205 0.24444054
265 266 267 268 269 270
-0.76518011 0.23481989 -0.76518011 0.26784252 -0.76518011 0.28466508
271 272 273 274 275 276
-0.76518011 -0.76518011 -0.76518011 0.26012593 0.23481989 0.24372533
277 278 279 280 281 282
0.26194223 0.23481989 -0.76520776 0.24244016 -0.76518011 0.23481989
283 284 285 286 287 288
0.24402542 0.23481989 0.25439490 0.23481989 -0.76092800 0.23956903
289 290 291 292 293 294
-0.76518011 0.22210868 -0.76518011 0.24745710 0.28070563 0.23481989
295 296 297 298 299 300
0.24277274 -0.75828352 0.23481989 -0.76518011 0.23481989 0.24951175
301 302 303 304 305 306
-0.76052556 0.23124935 0.23645232 -0.76518011 0.23989902 0.25085144
307 308 309 310 311 312
-0.76518011 -0.73683571 0.24639581 -0.76518011 -0.76518011 0.28091768
313 314 315 316 317 318
-0.76518011 -0.74176680 0.24635460 0.24304150 0.23359784 0.25776551
319 320 321 322 323 324
0.24059927 0.23481989 0.23481989 0.23481989 0.23481989 -0.76518011
325 326 327 328 329 330
0.23481989 0.23524808 -0.76518011 -0.76518011 0.23481989 0.23481989
331 332 333 334 335 336
0.23481989 0.23481989 0.23481989 0.23854304 0.23566384 0.23832959
337 338 339 340 341 342
0.23167988 0.23481989 0.23595377 0.24495064 0.24717650 0.23481989
343 344 345 346 347 348
0.23481989 0.23481989 0.23412490 0.23481989 0.23409239 0.23481989
349 350 351 352 353 354
0.23481989 0.23960937 0.23158593 0.24229854 0.24464039 0.23630761
355 356 357 358 359 360
0.23481989 0.23481989 0.25438006 0.23481989 0.23481989 0.23481989
361 362 363 364 365 366
0.23976520 0.23481989 0.23481989 -0.73483613 0.24368290 0.23481989
367 368 369 370 371 372
0.23481989 0.23481989 0.23481989 0.23722540 0.23481989 0.22900052
373 374 375 376 377 378
0.25904003 0.24156077 0.23481989 0.23481989 0.23633457 0.24044613
379 380 381 382 383 384
0.23881790 0.23481989 0.23481989 0.25055216 0.20230949 0.23481989
385 386 387 388 389 390
0.24985634 0.25542715 0.23732079 0.23571886 0.23785683 0.23680560
391 392 393 394 395 396
0.24851588 0.23592627 0.23709030 0.25206781 -0.74930377 0.23503095
397 398 399 400 401 402
0.22862939 -0.35657770 0.24155472 -0.80639508 0.24179596 0.24711152
403 404 405 406 407 408
0.23583698 -0.74224893 0.24021251 0.25423588 0.25125719 0.24385971
409 410 411 412 413 414
0.23945931 0.23472385 0.30706299 -0.76794935 -0.75979571 -0.70792623
415 416 417 418 419 420
0.23986259 0.23225704 0.22298226 -0.74453388 0.26252948 0.22969189
421 422 423 424 425 426
0.23327328 0.32181336 -0.83618760 0.18877903 0.24034166 -0.69269676
427 428 429 430 431
0.20303590 0.51444855 0.24735795 0.31129886 0.14756937
> postscript(file="/var/www/html/rcomp/tmp/6wkcj1291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 431
Frequency = 1
lag(myerror, k = 1) myerror
0 0.43293581 NA
1 0.23239332 0.43293581
2 -0.02906483 0.23239332
3 0.33008759 -0.02906483
4 -0.35131265 0.33008759
5 -0.44957037 -0.35131265
6 0.52943665 -0.44957037
7 -0.26818786 0.52943665
8 -0.61779874 -0.26818786
9 -0.52978349 -0.61779874
10 0.53276565 -0.52978349
11 0.41613882 0.53276565
12 -0.20601138 0.41613882
13 0.26699770 -0.20601138
14 -0.48578614 0.26699770
15 0.42550073 -0.48578614
16 0.58132657 0.42550073
17 -0.49751323 0.58132657
18 -0.54645068 -0.49751323
19 0.55252200 -0.54645068
20 -0.65210778 0.55252200
21 -0.69843164 -0.65210778
22 0.46039586 -0.69843164
23 0.59395424 0.46039586
24 0.37020477 0.59395424
25 0.28345647 0.37020477
26 -0.52943794 0.28345647
27 0.20277225 -0.52943794
28 0.27699963 0.20277225
29 0.27581767 0.27699963
30 -0.67651550 0.27581767
31 -0.65138491 -0.67651550
32 -0.67215069 -0.65138491
33 -0.65957760 -0.67215069
34 0.49809064 -0.65957760
35 0.50202496 0.49809064
36 0.34943671 0.50202496
37 -0.59164014 0.34943671
38 -0.52085581 -0.59164014
39 0.38246263 -0.52085581
40 0.25914808 0.38246263
41 0.29307248 0.25914808
42 0.32799116 0.29307248
43 0.37386669 0.32799116
44 0.27489395 0.37386669
45 -0.70149628 0.27489395
46 -0.64346182 -0.70149628
47 -0.59072091 -0.64346182
48 0.35430531 -0.59072091
49 0.28118013 0.35430531
50 -0.60259299 0.28118013
51 -0.52443854 -0.60259299
52 -0.53821088 -0.52443854
53 -0.62687094 -0.53821088
54 -0.57888273 -0.62687094
55 0.32550517 -0.57888273
56 0.25913457 0.32550517
57 0.30787744 0.25913457
58 0.31002748 0.30787744
59 -0.66423740 0.31002748
60 -0.68009488 -0.66423740
61 0.24658854 -0.68009488
62 -0.70531491 0.24658854
63 -0.67326022 -0.70531491
64 0.28597908 -0.67326022
65 0.25811008 0.28597908
66 0.23905231 0.25811008
67 -0.67507442 0.23905231
68 0.38359411 -0.67507442
69 -0.61551461 0.38359411
70 -0.76225453 -0.61551461
71 0.24812997 -0.76225453
72 0.31551042 0.24812997
73 0.26156090 0.31551042
74 0.33259396 0.26156090
75 -0.73889110 0.33259396
76 0.31595963 -0.73889110
77 0.23593282 0.31595963
78 0.27677127 0.23593282
79 -0.75865595 0.27677127
80 0.27132753 -0.75865595
81 -0.66253050 0.27132753
82 0.33855867 -0.66253050
83 0.22855453 0.33855867
84 0.29751384 0.22855453
85 -0.69332482 0.29751384
86 0.39595052 -0.69332482
87 0.29293311 0.39595052
88 0.31723574 0.29293311
89 0.29888932 0.31723574
90 0.30751962 0.29888932
91 -0.69807119 0.30751962
92 0.31801737 -0.69807119
93 -0.70263724 0.31801737
94 0.26481137 -0.70263724
95 -0.73652044 0.26481137
96 0.30681413 -0.73652044
97 -0.80193158 0.30681413
98 0.32270161 -0.80193158
99 -0.74766448 0.32270161
100 -0.71864211 -0.74766448
101 0.27131681 -0.71864211
102 0.32363232 0.27131681
103 0.28729268 0.32363232
104 0.27266835 0.28729268
105 0.27510089 0.27266835
106 -0.74446203 0.27510089
107 -0.69653640 -0.74446203
108 0.30665270 -0.69653640
109 -0.74972558 0.30665270
110 0.27344052 -0.74972558
111 -0.77150768 0.27344052
112 0.28207300 -0.77150768
113 -0.74830300 0.28207300
114 -0.75904337 -0.74830300
115 0.27086835 -0.75904337
116 0.26054414 0.27086835
117 0.28748925 0.26054414
118 -0.75782506 0.28748925
119 0.23952578 -0.75782506
120 0.36611330 0.23952578
121 0.25396284 0.36611330
122 0.40447888 0.25396284
123 -0.67886500 0.40447888
124 -0.75435964 -0.67886500
125 -0.74541345 -0.75435964
126 0.23481989 -0.74541345
127 0.24802628 0.23481989
128 0.25764067 0.24802628
129 0.23481989 0.25764067
130 0.24276031 0.23481989
131 0.24726269 0.24276031
132 0.23481989 0.24726269
133 0.26599053 0.23481989
134 0.23684546 0.26599053
135 0.23649803 0.23684546
136 0.24623610 0.23649803
137 0.23558488 0.24623610
138 0.23634882 0.23558488
139 0.23481989 0.23634882
140 0.24066004 0.23481989
141 0.23481989 0.24066004
142 0.23789258 0.23481989
143 0.24025505 0.23789258
144 0.25595443 0.24025505
145 0.24000048 0.25595443
146 0.23481989 0.24000048
147 0.24389839 0.23481989
148 0.23756545 0.24389839
149 0.23470745 0.23756545
150 0.23585876 0.23470745
151 0.27377432 0.23585876
152 0.24298203 0.27377432
153 0.23636865 0.24298203
154 0.23701363 0.23636865
155 0.23481989 0.23701363
156 0.23481989 0.23481989
157 0.23765691 0.23481989
158 0.24147668 0.23765691
159 0.23481989 0.24147668
160 0.23481989 0.23481989
161 0.23481989 0.23481989
162 0.23481989 0.23481989
163 0.23297331 0.23481989
164 0.25305441 0.23297331
165 0.23504923 0.25305441
166 0.23481989 0.23504923
167 0.24837324 0.23481989
168 0.15658986 0.24837324
169 0.23481989 0.15658986
170 0.23481989 0.23481989
171 0.23843506 0.23481989
172 0.23481989 0.23843506
173 0.24570568 0.23481989
174 0.23481989 0.24570568
175 0.24373175 0.23481989
176 0.35747989 0.24373175
177 -0.76569382 0.35747989
178 0.23481989 -0.76569382
179 0.23481989 0.23481989
180 0.26436488 0.23481989
181 -0.74119379 0.26436488
182 0.24026604 -0.74119379
183 -0.75282388 0.24026604
184 -0.75746177 -0.75282388
185 -0.75308857 -0.75746177
186 -0.73857113 -0.75308857
187 0.24310873 -0.73857113
188 0.24206452 0.24310873
189 0.23481989 0.24206452
190 0.24303364 0.23481989
191 0.23481989 0.24303364
192 0.24423558 0.23481989
193 0.23481989 0.24423558
194 0.23902437 0.23481989
195 -0.76518011 0.23902437
196 -0.61286300 -0.76518011
197 -0.76518011 -0.61286300
198 -0.76518011 -0.76518011
199 -0.73136359 -0.76518011
200 -0.76518011 -0.73136359
201 -0.76518011 -0.76518011
202 -0.76518011 -0.76518011
203 -0.76775502 -0.76518011
204 -0.72833137 -0.76775502
205 -0.75689406 -0.72833137
206 -0.75877780 -0.75689406
207 -0.76518011 -0.75877780
208 -0.74700553 -0.76518011
209 -0.77444136 -0.74700553
210 0.23481989 -0.77444136
211 0.23481989 0.23481989
212 0.23812980 0.23481989
213 0.23967849 0.23812980
214 0.23889437 0.23967849
215 0.23808475 0.23889437
216 0.24153295 0.23808475
217 0.23405074 0.24153295
218 0.23808475 0.23405074
219 0.25957415 0.23808475
220 0.23477055 0.25957415
221 0.23481989 0.23477055
222 0.23562663 0.23481989
223 0.25721916 0.23562663
224 0.23569479 0.25721916
225 0.26008430 0.23569479
226 0.24366427 0.26008430
227 0.23481989 0.24366427
228 0.23967999 0.23481989
229 0.23481989 0.23967999
230 0.24258536 0.23481989
231 0.24382045 0.24258536
232 0.23481989 0.24382045
233 0.24753519 0.23481989
234 0.23481989 0.24753519
235 0.23445076 0.23481989
236 0.23481989 0.23445076
237 0.24177067 0.23481989
238 0.23481989 0.24177067
239 0.23481989 0.23481989
240 0.23481989 0.23481989
241 0.23850898 0.23481989
242 0.23481989 0.23850898
243 0.24220134 0.23481989
244 0.24081959 0.24220134
245 -0.76518011 0.24081959
246 -0.76518011 -0.76518011
247 0.25475464 -0.76518011
248 -0.76518011 0.25475464
249 -0.76518011 -0.76518011
250 0.23607919 -0.76518011
251 -0.76518011 0.23607919
252 0.23738379 -0.76518011
253 0.27582543 0.23738379
254 0.29577938 0.27582543
255 0.23481989 0.29577938
256 0.23577488 0.23481989
257 0.24826147 0.23577488
258 -0.76518011 0.24826147
259 0.23481989 -0.76518011
260 -0.76518011 0.23481989
261 -0.76518011 -0.76518011
262 0.26502205 -0.76518011
263 0.24444054 0.26502205
264 -0.76518011 0.24444054
265 0.23481989 -0.76518011
266 -0.76518011 0.23481989
267 0.26784252 -0.76518011
268 -0.76518011 0.26784252
269 0.28466508 -0.76518011
270 -0.76518011 0.28466508
271 -0.76518011 -0.76518011
272 -0.76518011 -0.76518011
273 0.26012593 -0.76518011
274 0.23481989 0.26012593
275 0.24372533 0.23481989
276 0.26194223 0.24372533
277 0.23481989 0.26194223
278 -0.76520776 0.23481989
279 0.24244016 -0.76520776
280 -0.76518011 0.24244016
281 0.23481989 -0.76518011
282 0.24402542 0.23481989
283 0.23481989 0.24402542
284 0.25439490 0.23481989
285 0.23481989 0.25439490
286 -0.76092800 0.23481989
287 0.23956903 -0.76092800
288 -0.76518011 0.23956903
289 0.22210868 -0.76518011
290 -0.76518011 0.22210868
291 0.24745710 -0.76518011
292 0.28070563 0.24745710
293 0.23481989 0.28070563
294 0.24277274 0.23481989
295 -0.75828352 0.24277274
296 0.23481989 -0.75828352
297 -0.76518011 0.23481989
298 0.23481989 -0.76518011
299 0.24951175 0.23481989
300 -0.76052556 0.24951175
301 0.23124935 -0.76052556
302 0.23645232 0.23124935
303 -0.76518011 0.23645232
304 0.23989902 -0.76518011
305 0.25085144 0.23989902
306 -0.76518011 0.25085144
307 -0.73683571 -0.76518011
308 0.24639581 -0.73683571
309 -0.76518011 0.24639581
310 -0.76518011 -0.76518011
311 0.28091768 -0.76518011
312 -0.76518011 0.28091768
313 -0.74176680 -0.76518011
314 0.24635460 -0.74176680
315 0.24304150 0.24635460
316 0.23359784 0.24304150
317 0.25776551 0.23359784
318 0.24059927 0.25776551
319 0.23481989 0.24059927
320 0.23481989 0.23481989
321 0.23481989 0.23481989
322 0.23481989 0.23481989
323 -0.76518011 0.23481989
324 0.23481989 -0.76518011
325 0.23524808 0.23481989
326 -0.76518011 0.23524808
327 -0.76518011 -0.76518011
328 0.23481989 -0.76518011
329 0.23481989 0.23481989
330 0.23481989 0.23481989
331 0.23481989 0.23481989
332 0.23481989 0.23481989
333 0.23854304 0.23481989
334 0.23566384 0.23854304
335 0.23832959 0.23566384
336 0.23167988 0.23832959
337 0.23481989 0.23167988
338 0.23595377 0.23481989
339 0.24495064 0.23595377
340 0.24717650 0.24495064
341 0.23481989 0.24717650
342 0.23481989 0.23481989
343 0.23481989 0.23481989
344 0.23412490 0.23481989
345 0.23481989 0.23412490
346 0.23409239 0.23481989
347 0.23481989 0.23409239
348 0.23481989 0.23481989
349 0.23960937 0.23481989
350 0.23158593 0.23960937
351 0.24229854 0.23158593
352 0.24464039 0.24229854
353 0.23630761 0.24464039
354 0.23481989 0.23630761
355 0.23481989 0.23481989
356 0.25438006 0.23481989
357 0.23481989 0.25438006
358 0.23481989 0.23481989
359 0.23481989 0.23481989
360 0.23976520 0.23481989
361 0.23481989 0.23976520
362 0.23481989 0.23481989
363 -0.73483613 0.23481989
364 0.24368290 -0.73483613
365 0.23481989 0.24368290
366 0.23481989 0.23481989
367 0.23481989 0.23481989
368 0.23481989 0.23481989
369 0.23722540 0.23481989
370 0.23481989 0.23722540
371 0.22900052 0.23481989
372 0.25904003 0.22900052
373 0.24156077 0.25904003
374 0.23481989 0.24156077
375 0.23481989 0.23481989
376 0.23633457 0.23481989
377 0.24044613 0.23633457
378 0.23881790 0.24044613
379 0.23481989 0.23881790
380 0.23481989 0.23481989
381 0.25055216 0.23481989
382 0.20230949 0.25055216
383 0.23481989 0.20230949
384 0.24985634 0.23481989
385 0.25542715 0.24985634
386 0.23732079 0.25542715
387 0.23571886 0.23732079
388 0.23785683 0.23571886
389 0.23680560 0.23785683
390 0.24851588 0.23680560
391 0.23592627 0.24851588
392 0.23709030 0.23592627
393 0.25206781 0.23709030
394 -0.74930377 0.25206781
395 0.23503095 -0.74930377
396 0.22862939 0.23503095
397 -0.35657770 0.22862939
398 0.24155472 -0.35657770
399 -0.80639508 0.24155472
400 0.24179596 -0.80639508
401 0.24711152 0.24179596
402 0.23583698 0.24711152
403 -0.74224893 0.23583698
404 0.24021251 -0.74224893
405 0.25423588 0.24021251
406 0.25125719 0.25423588
407 0.24385971 0.25125719
408 0.23945931 0.24385971
409 0.23472385 0.23945931
410 0.30706299 0.23472385
411 -0.76794935 0.30706299
412 -0.75979571 -0.76794935
413 -0.70792623 -0.75979571
414 0.23986259 -0.70792623
415 0.23225704 0.23986259
416 0.22298226 0.23225704
417 -0.74453388 0.22298226
418 0.26252948 -0.74453388
419 0.22969189 0.26252948
420 0.23327328 0.22969189
421 0.32181336 0.23327328
422 -0.83618760 0.32181336
423 0.18877903 -0.83618760
424 0.24034166 0.18877903
425 -0.69269676 0.24034166
426 0.20303590 -0.69269676
427 0.51444855 0.20303590
428 0.24735795 0.51444855
429 0.31129886 0.24735795
430 0.14756937 0.31129886
431 NA 0.14756937
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 0.23239332 0.43293581
[2,] -0.02906483 0.23239332
[3,] 0.33008759 -0.02906483
[4,] -0.35131265 0.33008759
[5,] -0.44957037 -0.35131265
[6,] 0.52943665 -0.44957037
[7,] -0.26818786 0.52943665
[8,] -0.61779874 -0.26818786
[9,] -0.52978349 -0.61779874
[10,] 0.53276565 -0.52978349
[11,] 0.41613882 0.53276565
[12,] -0.20601138 0.41613882
[13,] 0.26699770 -0.20601138
[14,] -0.48578614 0.26699770
[15,] 0.42550073 -0.48578614
[16,] 0.58132657 0.42550073
[17,] -0.49751323 0.58132657
[18,] -0.54645068 -0.49751323
[19,] 0.55252200 -0.54645068
[20,] -0.65210778 0.55252200
[21,] -0.69843164 -0.65210778
[22,] 0.46039586 -0.69843164
[23,] 0.59395424 0.46039586
[24,] 0.37020477 0.59395424
[25,] 0.28345647 0.37020477
[26,] -0.52943794 0.28345647
[27,] 0.20277225 -0.52943794
[28,] 0.27699963 0.20277225
[29,] 0.27581767 0.27699963
[30,] -0.67651550 0.27581767
[31,] -0.65138491 -0.67651550
[32,] -0.67215069 -0.65138491
[33,] -0.65957760 -0.67215069
[34,] 0.49809064 -0.65957760
[35,] 0.50202496 0.49809064
[36,] 0.34943671 0.50202496
[37,] -0.59164014 0.34943671
[38,] -0.52085581 -0.59164014
[39,] 0.38246263 -0.52085581
[40,] 0.25914808 0.38246263
[41,] 0.29307248 0.25914808
[42,] 0.32799116 0.29307248
[43,] 0.37386669 0.32799116
[44,] 0.27489395 0.37386669
[45,] -0.70149628 0.27489395
[46,] -0.64346182 -0.70149628
[47,] -0.59072091 -0.64346182
[48,] 0.35430531 -0.59072091
[49,] 0.28118013 0.35430531
[50,] -0.60259299 0.28118013
[51,] -0.52443854 -0.60259299
[52,] -0.53821088 -0.52443854
[53,] -0.62687094 -0.53821088
[54,] -0.57888273 -0.62687094
[55,] 0.32550517 -0.57888273
[56,] 0.25913457 0.32550517
[57,] 0.30787744 0.25913457
[58,] 0.31002748 0.30787744
[59,] -0.66423740 0.31002748
[60,] -0.68009488 -0.66423740
[61,] 0.24658854 -0.68009488
[62,] -0.70531491 0.24658854
[63,] -0.67326022 -0.70531491
[64,] 0.28597908 -0.67326022
[65,] 0.25811008 0.28597908
[66,] 0.23905231 0.25811008
[67,] -0.67507442 0.23905231
[68,] 0.38359411 -0.67507442
[69,] -0.61551461 0.38359411
[70,] -0.76225453 -0.61551461
[71,] 0.24812997 -0.76225453
[72,] 0.31551042 0.24812997
[73,] 0.26156090 0.31551042
[74,] 0.33259396 0.26156090
[75,] -0.73889110 0.33259396
[76,] 0.31595963 -0.73889110
[77,] 0.23593282 0.31595963
[78,] 0.27677127 0.23593282
[79,] -0.75865595 0.27677127
[80,] 0.27132753 -0.75865595
[81,] -0.66253050 0.27132753
[82,] 0.33855867 -0.66253050
[83,] 0.22855453 0.33855867
[84,] 0.29751384 0.22855453
[85,] -0.69332482 0.29751384
[86,] 0.39595052 -0.69332482
[87,] 0.29293311 0.39595052
[88,] 0.31723574 0.29293311
[89,] 0.29888932 0.31723574
[90,] 0.30751962 0.29888932
[91,] -0.69807119 0.30751962
[92,] 0.31801737 -0.69807119
[93,] -0.70263724 0.31801737
[94,] 0.26481137 -0.70263724
[95,] -0.73652044 0.26481137
[96,] 0.30681413 -0.73652044
[97,] -0.80193158 0.30681413
[98,] 0.32270161 -0.80193158
[99,] -0.74766448 0.32270161
[100,] -0.71864211 -0.74766448
[101,] 0.27131681 -0.71864211
[102,] 0.32363232 0.27131681
[103,] 0.28729268 0.32363232
[104,] 0.27266835 0.28729268
[105,] 0.27510089 0.27266835
[106,] -0.74446203 0.27510089
[107,] -0.69653640 -0.74446203
[108,] 0.30665270 -0.69653640
[109,] -0.74972558 0.30665270
[110,] 0.27344052 -0.74972558
[111,] -0.77150768 0.27344052
[112,] 0.28207300 -0.77150768
[113,] -0.74830300 0.28207300
[114,] -0.75904337 -0.74830300
[115,] 0.27086835 -0.75904337
[116,] 0.26054414 0.27086835
[117,] 0.28748925 0.26054414
[118,] -0.75782506 0.28748925
[119,] 0.23952578 -0.75782506
[120,] 0.36611330 0.23952578
[121,] 0.25396284 0.36611330
[122,] 0.40447888 0.25396284
[123,] -0.67886500 0.40447888
[124,] -0.75435964 -0.67886500
[125,] -0.74541345 -0.75435964
[126,] 0.23481989 -0.74541345
[127,] 0.24802628 0.23481989
[128,] 0.25764067 0.24802628
[129,] 0.23481989 0.25764067
[130,] 0.24276031 0.23481989
[131,] 0.24726269 0.24276031
[132,] 0.23481989 0.24726269
[133,] 0.26599053 0.23481989
[134,] 0.23684546 0.26599053
[135,] 0.23649803 0.23684546
[136,] 0.24623610 0.23649803
[137,] 0.23558488 0.24623610
[138,] 0.23634882 0.23558488
[139,] 0.23481989 0.23634882
[140,] 0.24066004 0.23481989
[141,] 0.23481989 0.24066004
[142,] 0.23789258 0.23481989
[143,] 0.24025505 0.23789258
[144,] 0.25595443 0.24025505
[145,] 0.24000048 0.25595443
[146,] 0.23481989 0.24000048
[147,] 0.24389839 0.23481989
[148,] 0.23756545 0.24389839
[149,] 0.23470745 0.23756545
[150,] 0.23585876 0.23470745
[151,] 0.27377432 0.23585876
[152,] 0.24298203 0.27377432
[153,] 0.23636865 0.24298203
[154,] 0.23701363 0.23636865
[155,] 0.23481989 0.23701363
[156,] 0.23481989 0.23481989
[157,] 0.23765691 0.23481989
[158,] 0.24147668 0.23765691
[159,] 0.23481989 0.24147668
[160,] 0.23481989 0.23481989
[161,] 0.23481989 0.23481989
[162,] 0.23481989 0.23481989
[163,] 0.23297331 0.23481989
[164,] 0.25305441 0.23297331
[165,] 0.23504923 0.25305441
[166,] 0.23481989 0.23504923
[167,] 0.24837324 0.23481989
[168,] 0.15658986 0.24837324
[169,] 0.23481989 0.15658986
[170,] 0.23481989 0.23481989
[171,] 0.23843506 0.23481989
[172,] 0.23481989 0.23843506
[173,] 0.24570568 0.23481989
[174,] 0.23481989 0.24570568
[175,] 0.24373175 0.23481989
[176,] 0.35747989 0.24373175
[177,] -0.76569382 0.35747989
[178,] 0.23481989 -0.76569382
[179,] 0.23481989 0.23481989
[180,] 0.26436488 0.23481989
[181,] -0.74119379 0.26436488
[182,] 0.24026604 -0.74119379
[183,] -0.75282388 0.24026604
[184,] -0.75746177 -0.75282388
[185,] -0.75308857 -0.75746177
[186,] -0.73857113 -0.75308857
[187,] 0.24310873 -0.73857113
[188,] 0.24206452 0.24310873
[189,] 0.23481989 0.24206452
[190,] 0.24303364 0.23481989
[191,] 0.23481989 0.24303364
[192,] 0.24423558 0.23481989
[193,] 0.23481989 0.24423558
[194,] 0.23902437 0.23481989
[195,] -0.76518011 0.23902437
[196,] -0.61286300 -0.76518011
[197,] -0.76518011 -0.61286300
[198,] -0.76518011 -0.76518011
[199,] -0.73136359 -0.76518011
[200,] -0.76518011 -0.73136359
[201,] -0.76518011 -0.76518011
[202,] -0.76518011 -0.76518011
[203,] -0.76775502 -0.76518011
[204,] -0.72833137 -0.76775502
[205,] -0.75689406 -0.72833137
[206,] -0.75877780 -0.75689406
[207,] -0.76518011 -0.75877780
[208,] -0.74700553 -0.76518011
[209,] -0.77444136 -0.74700553
[210,] 0.23481989 -0.77444136
[211,] 0.23481989 0.23481989
[212,] 0.23812980 0.23481989
[213,] 0.23967849 0.23812980
[214,] 0.23889437 0.23967849
[215,] 0.23808475 0.23889437
[216,] 0.24153295 0.23808475
[217,] 0.23405074 0.24153295
[218,] 0.23808475 0.23405074
[219,] 0.25957415 0.23808475
[220,] 0.23477055 0.25957415
[221,] 0.23481989 0.23477055
[222,] 0.23562663 0.23481989
[223,] 0.25721916 0.23562663
[224,] 0.23569479 0.25721916
[225,] 0.26008430 0.23569479
[226,] 0.24366427 0.26008430
[227,] 0.23481989 0.24366427
[228,] 0.23967999 0.23481989
[229,] 0.23481989 0.23967999
[230,] 0.24258536 0.23481989
[231,] 0.24382045 0.24258536
[232,] 0.23481989 0.24382045
[233,] 0.24753519 0.23481989
[234,] 0.23481989 0.24753519
[235,] 0.23445076 0.23481989
[236,] 0.23481989 0.23445076
[237,] 0.24177067 0.23481989
[238,] 0.23481989 0.24177067
[239,] 0.23481989 0.23481989
[240,] 0.23481989 0.23481989
[241,] 0.23850898 0.23481989
[242,] 0.23481989 0.23850898
[243,] 0.24220134 0.23481989
[244,] 0.24081959 0.24220134
[245,] -0.76518011 0.24081959
[246,] -0.76518011 -0.76518011
[247,] 0.25475464 -0.76518011
[248,] -0.76518011 0.25475464
[249,] -0.76518011 -0.76518011
[250,] 0.23607919 -0.76518011
[251,] -0.76518011 0.23607919
[252,] 0.23738379 -0.76518011
[253,] 0.27582543 0.23738379
[254,] 0.29577938 0.27582543
[255,] 0.23481989 0.29577938
[256,] 0.23577488 0.23481989
[257,] 0.24826147 0.23577488
[258,] -0.76518011 0.24826147
[259,] 0.23481989 -0.76518011
[260,] -0.76518011 0.23481989
[261,] -0.76518011 -0.76518011
[262,] 0.26502205 -0.76518011
[263,] 0.24444054 0.26502205
[264,] -0.76518011 0.24444054
[265,] 0.23481989 -0.76518011
[266,] -0.76518011 0.23481989
[267,] 0.26784252 -0.76518011
[268,] -0.76518011 0.26784252
[269,] 0.28466508 -0.76518011
[270,] -0.76518011 0.28466508
[271,] -0.76518011 -0.76518011
[272,] -0.76518011 -0.76518011
[273,] 0.26012593 -0.76518011
[274,] 0.23481989 0.26012593
[275,] 0.24372533 0.23481989
[276,] 0.26194223 0.24372533
[277,] 0.23481989 0.26194223
[278,] -0.76520776 0.23481989
[279,] 0.24244016 -0.76520776
[280,] -0.76518011 0.24244016
[281,] 0.23481989 -0.76518011
[282,] 0.24402542 0.23481989
[283,] 0.23481989 0.24402542
[284,] 0.25439490 0.23481989
[285,] 0.23481989 0.25439490
[286,] -0.76092800 0.23481989
[287,] 0.23956903 -0.76092800
[288,] -0.76518011 0.23956903
[289,] 0.22210868 -0.76518011
[290,] -0.76518011 0.22210868
[291,] 0.24745710 -0.76518011
[292,] 0.28070563 0.24745710
[293,] 0.23481989 0.28070563
[294,] 0.24277274 0.23481989
[295,] -0.75828352 0.24277274
[296,] 0.23481989 -0.75828352
[297,] -0.76518011 0.23481989
[298,] 0.23481989 -0.76518011
[299,] 0.24951175 0.23481989
[300,] -0.76052556 0.24951175
[301,] 0.23124935 -0.76052556
[302,] 0.23645232 0.23124935
[303,] -0.76518011 0.23645232
[304,] 0.23989902 -0.76518011
[305,] 0.25085144 0.23989902
[306,] -0.76518011 0.25085144
[307,] -0.73683571 -0.76518011
[308,] 0.24639581 -0.73683571
[309,] -0.76518011 0.24639581
[310,] -0.76518011 -0.76518011
[311,] 0.28091768 -0.76518011
[312,] -0.76518011 0.28091768
[313,] -0.74176680 -0.76518011
[314,] 0.24635460 -0.74176680
[315,] 0.24304150 0.24635460
[316,] 0.23359784 0.24304150
[317,] 0.25776551 0.23359784
[318,] 0.24059927 0.25776551
[319,] 0.23481989 0.24059927
[320,] 0.23481989 0.23481989
[321,] 0.23481989 0.23481989
[322,] 0.23481989 0.23481989
[323,] -0.76518011 0.23481989
[324,] 0.23481989 -0.76518011
[325,] 0.23524808 0.23481989
[326,] -0.76518011 0.23524808
[327,] -0.76518011 -0.76518011
[328,] 0.23481989 -0.76518011
[329,] 0.23481989 0.23481989
[330,] 0.23481989 0.23481989
[331,] 0.23481989 0.23481989
[332,] 0.23481989 0.23481989
[333,] 0.23854304 0.23481989
[334,] 0.23566384 0.23854304
[335,] 0.23832959 0.23566384
[336,] 0.23167988 0.23832959
[337,] 0.23481989 0.23167988
[338,] 0.23595377 0.23481989
[339,] 0.24495064 0.23595377
[340,] 0.24717650 0.24495064
[341,] 0.23481989 0.24717650
[342,] 0.23481989 0.23481989
[343,] 0.23481989 0.23481989
[344,] 0.23412490 0.23481989
[345,] 0.23481989 0.23412490
[346,] 0.23409239 0.23481989
[347,] 0.23481989 0.23409239
[348,] 0.23481989 0.23481989
[349,] 0.23960937 0.23481989
[350,] 0.23158593 0.23960937
[351,] 0.24229854 0.23158593
[352,] 0.24464039 0.24229854
[353,] 0.23630761 0.24464039
[354,] 0.23481989 0.23630761
[355,] 0.23481989 0.23481989
[356,] 0.25438006 0.23481989
[357,] 0.23481989 0.25438006
[358,] 0.23481989 0.23481989
[359,] 0.23481989 0.23481989
[360,] 0.23976520 0.23481989
[361,] 0.23481989 0.23976520
[362,] 0.23481989 0.23481989
[363,] -0.73483613 0.23481989
[364,] 0.24368290 -0.73483613
[365,] 0.23481989 0.24368290
[366,] 0.23481989 0.23481989
[367,] 0.23481989 0.23481989
[368,] 0.23481989 0.23481989
[369,] 0.23722540 0.23481989
[370,] 0.23481989 0.23722540
[371,] 0.22900052 0.23481989
[372,] 0.25904003 0.22900052
[373,] 0.24156077 0.25904003
[374,] 0.23481989 0.24156077
[375,] 0.23481989 0.23481989
[376,] 0.23633457 0.23481989
[377,] 0.24044613 0.23633457
[378,] 0.23881790 0.24044613
[379,] 0.23481989 0.23881790
[380,] 0.23481989 0.23481989
[381,] 0.25055216 0.23481989
[382,] 0.20230949 0.25055216
[383,] 0.23481989 0.20230949
[384,] 0.24985634 0.23481989
[385,] 0.25542715 0.24985634
[386,] 0.23732079 0.25542715
[387,] 0.23571886 0.23732079
[388,] 0.23785683 0.23571886
[389,] 0.23680560 0.23785683
[390,] 0.24851588 0.23680560
[391,] 0.23592627 0.24851588
[392,] 0.23709030 0.23592627
[393,] 0.25206781 0.23709030
[394,] -0.74930377 0.25206781
[395,] 0.23503095 -0.74930377
[396,] 0.22862939 0.23503095
[397,] -0.35657770 0.22862939
[398,] 0.24155472 -0.35657770
[399,] -0.80639508 0.24155472
[400,] 0.24179596 -0.80639508
[401,] 0.24711152 0.24179596
[402,] 0.23583698 0.24711152
[403,] -0.74224893 0.23583698
[404,] 0.24021251 -0.74224893
[405,] 0.25423588 0.24021251
[406,] 0.25125719 0.25423588
[407,] 0.24385971 0.25125719
[408,] 0.23945931 0.24385971
[409,] 0.23472385 0.23945931
[410,] 0.30706299 0.23472385
[411,] -0.76794935 0.30706299
[412,] -0.75979571 -0.76794935
[413,] -0.70792623 -0.75979571
[414,] 0.23986259 -0.70792623
[415,] 0.23225704 0.23986259
[416,] 0.22298226 0.23225704
[417,] -0.74453388 0.22298226
[418,] 0.26252948 -0.74453388
[419,] 0.22969189 0.26252948
[420,] 0.23327328 0.22969189
[421,] 0.32181336 0.23327328
[422,] -0.83618760 0.32181336
[423,] 0.18877903 -0.83618760
[424,] 0.24034166 0.18877903
[425,] -0.69269676 0.24034166
[426,] 0.20303590 -0.69269676
[427,] 0.51444855 0.20303590
[428,] 0.24735795 0.51444855
[429,] 0.31129886 0.24735795
[430,] 0.14756937 0.31129886
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 0.23239332 0.43293581
2 -0.02906483 0.23239332
3 0.33008759 -0.02906483
4 -0.35131265 0.33008759
5 -0.44957037 -0.35131265
6 0.52943665 -0.44957037
7 -0.26818786 0.52943665
8 -0.61779874 -0.26818786
9 -0.52978349 -0.61779874
10 0.53276565 -0.52978349
11 0.41613882 0.53276565
12 -0.20601138 0.41613882
13 0.26699770 -0.20601138
14 -0.48578614 0.26699770
15 0.42550073 -0.48578614
16 0.58132657 0.42550073
17 -0.49751323 0.58132657
18 -0.54645068 -0.49751323
19 0.55252200 -0.54645068
20 -0.65210778 0.55252200
21 -0.69843164 -0.65210778
22 0.46039586 -0.69843164
23 0.59395424 0.46039586
24 0.37020477 0.59395424
25 0.28345647 0.37020477
26 -0.52943794 0.28345647
27 0.20277225 -0.52943794
28 0.27699963 0.20277225
29 0.27581767 0.27699963
30 -0.67651550 0.27581767
31 -0.65138491 -0.67651550
32 -0.67215069 -0.65138491
33 -0.65957760 -0.67215069
34 0.49809064 -0.65957760
35 0.50202496 0.49809064
36 0.34943671 0.50202496
37 -0.59164014 0.34943671
38 -0.52085581 -0.59164014
39 0.38246263 -0.52085581
40 0.25914808 0.38246263
41 0.29307248 0.25914808
42 0.32799116 0.29307248
43 0.37386669 0.32799116
44 0.27489395 0.37386669
45 -0.70149628 0.27489395
46 -0.64346182 -0.70149628
47 -0.59072091 -0.64346182
48 0.35430531 -0.59072091
49 0.28118013 0.35430531
50 -0.60259299 0.28118013
51 -0.52443854 -0.60259299
52 -0.53821088 -0.52443854
53 -0.62687094 -0.53821088
54 -0.57888273 -0.62687094
55 0.32550517 -0.57888273
56 0.25913457 0.32550517
57 0.30787744 0.25913457
58 0.31002748 0.30787744
59 -0.66423740 0.31002748
60 -0.68009488 -0.66423740
61 0.24658854 -0.68009488
62 -0.70531491 0.24658854
63 -0.67326022 -0.70531491
64 0.28597908 -0.67326022
65 0.25811008 0.28597908
66 0.23905231 0.25811008
67 -0.67507442 0.23905231
68 0.38359411 -0.67507442
69 -0.61551461 0.38359411
70 -0.76225453 -0.61551461
71 0.24812997 -0.76225453
72 0.31551042 0.24812997
73 0.26156090 0.31551042
74 0.33259396 0.26156090
75 -0.73889110 0.33259396
76 0.31595963 -0.73889110
77 0.23593282 0.31595963
78 0.27677127 0.23593282
79 -0.75865595 0.27677127
80 0.27132753 -0.75865595
81 -0.66253050 0.27132753
82 0.33855867 -0.66253050
83 0.22855453 0.33855867
84 0.29751384 0.22855453
85 -0.69332482 0.29751384
86 0.39595052 -0.69332482
87 0.29293311 0.39595052
88 0.31723574 0.29293311
89 0.29888932 0.31723574
90 0.30751962 0.29888932
91 -0.69807119 0.30751962
92 0.31801737 -0.69807119
93 -0.70263724 0.31801737
94 0.26481137 -0.70263724
95 -0.73652044 0.26481137
96 0.30681413 -0.73652044
97 -0.80193158 0.30681413
98 0.32270161 -0.80193158
99 -0.74766448 0.32270161
100 -0.71864211 -0.74766448
101 0.27131681 -0.71864211
102 0.32363232 0.27131681
103 0.28729268 0.32363232
104 0.27266835 0.28729268
105 0.27510089 0.27266835
106 -0.74446203 0.27510089
107 -0.69653640 -0.74446203
108 0.30665270 -0.69653640
109 -0.74972558 0.30665270
110 0.27344052 -0.74972558
111 -0.77150768 0.27344052
112 0.28207300 -0.77150768
113 -0.74830300 0.28207300
114 -0.75904337 -0.74830300
115 0.27086835 -0.75904337
116 0.26054414 0.27086835
117 0.28748925 0.26054414
118 -0.75782506 0.28748925
119 0.23952578 -0.75782506
120 0.36611330 0.23952578
121 0.25396284 0.36611330
122 0.40447888 0.25396284
123 -0.67886500 0.40447888
124 -0.75435964 -0.67886500
125 -0.74541345 -0.75435964
126 0.23481989 -0.74541345
127 0.24802628 0.23481989
128 0.25764067 0.24802628
129 0.23481989 0.25764067
130 0.24276031 0.23481989
131 0.24726269 0.24276031
132 0.23481989 0.24726269
133 0.26599053 0.23481989
134 0.23684546 0.26599053
135 0.23649803 0.23684546
136 0.24623610 0.23649803
137 0.23558488 0.24623610
138 0.23634882 0.23558488
139 0.23481989 0.23634882
140 0.24066004 0.23481989
141 0.23481989 0.24066004
142 0.23789258 0.23481989
143 0.24025505 0.23789258
144 0.25595443 0.24025505
145 0.24000048 0.25595443
146 0.23481989 0.24000048
147 0.24389839 0.23481989
148 0.23756545 0.24389839
149 0.23470745 0.23756545
150 0.23585876 0.23470745
151 0.27377432 0.23585876
152 0.24298203 0.27377432
153 0.23636865 0.24298203
154 0.23701363 0.23636865
155 0.23481989 0.23701363
156 0.23481989 0.23481989
157 0.23765691 0.23481989
158 0.24147668 0.23765691
159 0.23481989 0.24147668
160 0.23481989 0.23481989
161 0.23481989 0.23481989
162 0.23481989 0.23481989
163 0.23297331 0.23481989
164 0.25305441 0.23297331
165 0.23504923 0.25305441
166 0.23481989 0.23504923
167 0.24837324 0.23481989
168 0.15658986 0.24837324
169 0.23481989 0.15658986
170 0.23481989 0.23481989
171 0.23843506 0.23481989
172 0.23481989 0.23843506
173 0.24570568 0.23481989
174 0.23481989 0.24570568
175 0.24373175 0.23481989
176 0.35747989 0.24373175
177 -0.76569382 0.35747989
178 0.23481989 -0.76569382
179 0.23481989 0.23481989
180 0.26436488 0.23481989
181 -0.74119379 0.26436488
182 0.24026604 -0.74119379
183 -0.75282388 0.24026604
184 -0.75746177 -0.75282388
185 -0.75308857 -0.75746177
186 -0.73857113 -0.75308857
187 0.24310873 -0.73857113
188 0.24206452 0.24310873
189 0.23481989 0.24206452
190 0.24303364 0.23481989
191 0.23481989 0.24303364
192 0.24423558 0.23481989
193 0.23481989 0.24423558
194 0.23902437 0.23481989
195 -0.76518011 0.23902437
196 -0.61286300 -0.76518011
197 -0.76518011 -0.61286300
198 -0.76518011 -0.76518011
199 -0.73136359 -0.76518011
200 -0.76518011 -0.73136359
201 -0.76518011 -0.76518011
202 -0.76518011 -0.76518011
203 -0.76775502 -0.76518011
204 -0.72833137 -0.76775502
205 -0.75689406 -0.72833137
206 -0.75877780 -0.75689406
207 -0.76518011 -0.75877780
208 -0.74700553 -0.76518011
209 -0.77444136 -0.74700553
210 0.23481989 -0.77444136
211 0.23481989 0.23481989
212 0.23812980 0.23481989
213 0.23967849 0.23812980
214 0.23889437 0.23967849
215 0.23808475 0.23889437
216 0.24153295 0.23808475
217 0.23405074 0.24153295
218 0.23808475 0.23405074
219 0.25957415 0.23808475
220 0.23477055 0.25957415
221 0.23481989 0.23477055
222 0.23562663 0.23481989
223 0.25721916 0.23562663
224 0.23569479 0.25721916
225 0.26008430 0.23569479
226 0.24366427 0.26008430
227 0.23481989 0.24366427
228 0.23967999 0.23481989
229 0.23481989 0.23967999
230 0.24258536 0.23481989
231 0.24382045 0.24258536
232 0.23481989 0.24382045
233 0.24753519 0.23481989
234 0.23481989 0.24753519
235 0.23445076 0.23481989
236 0.23481989 0.23445076
237 0.24177067 0.23481989
238 0.23481989 0.24177067
239 0.23481989 0.23481989
240 0.23481989 0.23481989
241 0.23850898 0.23481989
242 0.23481989 0.23850898
243 0.24220134 0.23481989
244 0.24081959 0.24220134
245 -0.76518011 0.24081959
246 -0.76518011 -0.76518011
247 0.25475464 -0.76518011
248 -0.76518011 0.25475464
249 -0.76518011 -0.76518011
250 0.23607919 -0.76518011
251 -0.76518011 0.23607919
252 0.23738379 -0.76518011
253 0.27582543 0.23738379
254 0.29577938 0.27582543
255 0.23481989 0.29577938
256 0.23577488 0.23481989
257 0.24826147 0.23577488
258 -0.76518011 0.24826147
259 0.23481989 -0.76518011
260 -0.76518011 0.23481989
261 -0.76518011 -0.76518011
262 0.26502205 -0.76518011
263 0.24444054 0.26502205
264 -0.76518011 0.24444054
265 0.23481989 -0.76518011
266 -0.76518011 0.23481989
267 0.26784252 -0.76518011
268 -0.76518011 0.26784252
269 0.28466508 -0.76518011
270 -0.76518011 0.28466508
271 -0.76518011 -0.76518011
272 -0.76518011 -0.76518011
273 0.26012593 -0.76518011
274 0.23481989 0.26012593
275 0.24372533 0.23481989
276 0.26194223 0.24372533
277 0.23481989 0.26194223
278 -0.76520776 0.23481989
279 0.24244016 -0.76520776
280 -0.76518011 0.24244016
281 0.23481989 -0.76518011
282 0.24402542 0.23481989
283 0.23481989 0.24402542
284 0.25439490 0.23481989
285 0.23481989 0.25439490
286 -0.76092800 0.23481989
287 0.23956903 -0.76092800
288 -0.76518011 0.23956903
289 0.22210868 -0.76518011
290 -0.76518011 0.22210868
291 0.24745710 -0.76518011
292 0.28070563 0.24745710
293 0.23481989 0.28070563
294 0.24277274 0.23481989
295 -0.75828352 0.24277274
296 0.23481989 -0.75828352
297 -0.76518011 0.23481989
298 0.23481989 -0.76518011
299 0.24951175 0.23481989
300 -0.76052556 0.24951175
301 0.23124935 -0.76052556
302 0.23645232 0.23124935
303 -0.76518011 0.23645232
304 0.23989902 -0.76518011
305 0.25085144 0.23989902
306 -0.76518011 0.25085144
307 -0.73683571 -0.76518011
308 0.24639581 -0.73683571
309 -0.76518011 0.24639581
310 -0.76518011 -0.76518011
311 0.28091768 -0.76518011
312 -0.76518011 0.28091768
313 -0.74176680 -0.76518011
314 0.24635460 -0.74176680
315 0.24304150 0.24635460
316 0.23359784 0.24304150
317 0.25776551 0.23359784
318 0.24059927 0.25776551
319 0.23481989 0.24059927
320 0.23481989 0.23481989
321 0.23481989 0.23481989
322 0.23481989 0.23481989
323 -0.76518011 0.23481989
324 0.23481989 -0.76518011
325 0.23524808 0.23481989
326 -0.76518011 0.23524808
327 -0.76518011 -0.76518011
328 0.23481989 -0.76518011
329 0.23481989 0.23481989
330 0.23481989 0.23481989
331 0.23481989 0.23481989
332 0.23481989 0.23481989
333 0.23854304 0.23481989
334 0.23566384 0.23854304
335 0.23832959 0.23566384
336 0.23167988 0.23832959
337 0.23481989 0.23167988
338 0.23595377 0.23481989
339 0.24495064 0.23595377
340 0.24717650 0.24495064
341 0.23481989 0.24717650
342 0.23481989 0.23481989
343 0.23481989 0.23481989
344 0.23412490 0.23481989
345 0.23481989 0.23412490
346 0.23409239 0.23481989
347 0.23481989 0.23409239
348 0.23481989 0.23481989
349 0.23960937 0.23481989
350 0.23158593 0.23960937
351 0.24229854 0.23158593
352 0.24464039 0.24229854
353 0.23630761 0.24464039
354 0.23481989 0.23630761
355 0.23481989 0.23481989
356 0.25438006 0.23481989
357 0.23481989 0.25438006
358 0.23481989 0.23481989
359 0.23481989 0.23481989
360 0.23976520 0.23481989
361 0.23481989 0.23976520
362 0.23481989 0.23481989
363 -0.73483613 0.23481989
364 0.24368290 -0.73483613
365 0.23481989 0.24368290
366 0.23481989 0.23481989
367 0.23481989 0.23481989
368 0.23481989 0.23481989
369 0.23722540 0.23481989
370 0.23481989 0.23722540
371 0.22900052 0.23481989
372 0.25904003 0.22900052
373 0.24156077 0.25904003
374 0.23481989 0.24156077
375 0.23481989 0.23481989
376 0.23633457 0.23481989
377 0.24044613 0.23633457
378 0.23881790 0.24044613
379 0.23481989 0.23881790
380 0.23481989 0.23481989
381 0.25055216 0.23481989
382 0.20230949 0.25055216
383 0.23481989 0.20230949
384 0.24985634 0.23481989
385 0.25542715 0.24985634
386 0.23732079 0.25542715
387 0.23571886 0.23732079
388 0.23785683 0.23571886
389 0.23680560 0.23785683
390 0.24851588 0.23680560
391 0.23592627 0.24851588
392 0.23709030 0.23592627
393 0.25206781 0.23709030
394 -0.74930377 0.25206781
395 0.23503095 -0.74930377
396 0.22862939 0.23503095
397 -0.35657770 0.22862939
398 0.24155472 -0.35657770
399 -0.80639508 0.24155472
400 0.24179596 -0.80639508
401 0.24711152 0.24179596
402 0.23583698 0.24711152
403 -0.74224893 0.23583698
404 0.24021251 -0.74224893
405 0.25423588 0.24021251
406 0.25125719 0.25423588
407 0.24385971 0.25125719
408 0.23945931 0.24385971
409 0.23472385 0.23945931
410 0.30706299 0.23472385
411 -0.76794935 0.30706299
412 -0.75979571 -0.76794935
413 -0.70792623 -0.75979571
414 0.23986259 -0.70792623
415 0.23225704 0.23986259
416 0.22298226 0.23225704
417 -0.74453388 0.22298226
418 0.26252948 -0.74453388
419 0.22969189 0.26252948
420 0.23327328 0.22969189
421 0.32181336 0.23327328
422 -0.83618760 0.32181336
423 0.18877903 -0.83618760
424 0.24034166 0.18877903
425 -0.69269676 0.24034166
426 0.20303590 -0.69269676
427 0.51444855 0.20303590
428 0.24735795 0.51444855
429 0.31129886 0.24735795
430 0.14756937 0.31129886
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7pbt41291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8z2t71291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/9z2t71291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/10z2t71291208775.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11ec8y1291208775.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/1263811291208775.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13d5nc1291208775.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/146vmx1291208775.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/15rwkl1291208775.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/16n60c1291208775.tab")
+ }
>
> try(system("convert tmp/13tdy1291208775.ps tmp/13tdy1291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/23tdy1291208775.ps tmp/23tdy1291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/33tdy1291208775.ps tmp/33tdy1291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/4wkcj1291208775.ps tmp/4wkcj1291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/5wkcj1291208775.ps tmp/5wkcj1291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/6wkcj1291208775.ps tmp/6wkcj1291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/7pbt41291208775.ps tmp/7pbt41291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/8z2t71291208775.ps tmp/8z2t71291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/9z2t71291208775.ps tmp/9z2t71291208775.png",intern=TRUE))
character(0)
> try(system("convert tmp/10z2t71291208775.ps tmp/10z2t71291208775.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.106 2.255 27.747