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(1
+ ,1
+ ,1
+ ,162556
+ ,162556
+ ,1081
+ ,1081
+ ,213118
+ ,213118
+ ,6282154
+ ,1
+ ,2
+ ,2
+ ,29790
+ ,29790
+ ,309
+ ,309
+ ,81767
+ ,81767
+ ,4321023
+ ,1
+ ,3
+ ,3
+ ,87550
+ ,87550
+ ,458
+ ,458
+ ,153198
+ ,153198
+ ,4111912
+ ,1
+ ,4
+ ,4
+ ,84738
+ ,84738
+ ,588
+ ,588
+ ,-26007
+ ,-26007
+ ,223193
+ ,1
+ ,5
+ ,5
+ ,54660
+ ,54660
+ ,302
+ ,302
+ ,126942
+ ,126942
+ ,1491348
+ ,0
+ ,6
+ ,0
+ ,42634
+ ,0
+ ,156
+ ,0
+ ,157214
+ ,0
+ ,1629616
+ ,1
+ ,7
+ ,7
+ ,40949
+ ,40949
+ ,481
+ ,481
+ ,129352
+ ,129352
+ ,1398893
+ ,0
+ ,8
+ ,0
+ ,45187
+ ,0
+ ,353
+ ,0
+ ,234817
+ ,0
+ ,1926517
+ ,1
+ ,9
+ ,9
+ ,37704
+ ,37704
+ ,452
+ ,452
+ ,60448
+ ,60448
+ ,983660
+ ,1
+ ,10
+ ,10
+ ,16275
+ ,16275
+ ,109
+ ,109
+ ,47818
+ ,47818
+ ,1443586
+ ,1
+ ,11
+ ,11
+ ,25830
+ ,25830
+ ,115
+ ,115
+ ,245546
+ ,245546
+ ,1073089
+ ,1
+ ,12
+ ,12
+ ,12679
+ ,12679
+ ,110
+ ,110
+ ,48020
+ ,48020
+ ,984885
+ ,1
+ ,13
+ ,13
+ ,18014
+ ,18014
+ ,239
+ ,239
+ ,-1710
+ ,-1710
+ ,1405225
+ ,1
+ ,14
+ ,14
+ ,43556
+ ,43556
+ ,247
+ ,247
+ ,32648
+ ,32648
+ ,227132
+ ,1
+ ,15
+ ,15
+ ,24811
+ ,24811
+ ,505
+ ,505
+ ,95350
+ ,95350
+ ,929118
+ ,0
+ ,16
+ ,0
+ ,6575
+ ,0
+ ,159
+ ,0
+ ,151352
+ ,0
+ ,1071292
+ ,0
+ ,17
+ ,0
+ ,7123
+ ,0
+ ,109
+ ,0
+ ,288170
+ ,0
+ ,638830
+ ,1
+ ,18
+ ,18
+ ,21950
+ ,21950
+ ,519
+ ,519
+ ,114337
+ ,114337
+ ,856956
+ ,1
+ ,19
+ ,19
+ ,37597
+ ,37597
+ ,248
+ ,248
+ ,37884
+ ,37884
+ ,992426
+ ,1
+ ,20
+ ,20
+ ,17821
+ ,17821
+ ,373
+ ,373
+ ,122844
+ ,122844
+ ,444477
+ ,1
+ ,21
+ ,21
+ ,12988
+ ,12988
+ ,119
+ ,119
+ ,82340
+ ,82340
+ ,857217
+ ,0
+ ,22
+ ,0
+ ,22330
+ ,0
+ ,84
+ ,0
+ ,79801
+ ,0
+ ,711969
+ ,0
+ ,23
+ ,0
+ ,13326
+ ,0
+ ,102
+ ,0
+ ,165548
+ ,0
+ ,702380
+ ,1
+ ,24
+ ,24
+ ,16189
+ ,16189
+ ,295
+ ,295
+ ,116384
+ ,116384
+ ,358589
+ ,1
+ ,25
+ ,25
+ ,7146
+ ,7146
+ ,105
+ ,105
+ ,134028
+ ,134028
+ ,297978
+ ,1
+ ,26
+ ,26
+ ,15824
+ ,15824
+ ,64
+ ,64
+ ,63838
+ ,63838
+ ,585715
+ ,1
+ ,27
+ ,27
+ ,27664
+ ,27664
+ ,282
+ ,282
+ ,74996
+ ,74996
+ ,657954
+ ,1
+ ,28
+ ,28
+ ,11920
+ ,11920
+ ,182
+ ,182
+ ,31080
+ ,31080
+ ,209458
+ ,1
+ ,29
+ ,29
+ ,8568
+ ,8568
+ ,37
+ ,37
+ ,32168
+ ,32168
+ ,786690
+ ,1
+ ,30
+ ,30
+ ,14416
+ ,14416
+ ,361
+ ,361
+ ,49857
+ ,49857
+ ,439798
+ ,1
+ ,31
+ ,31
+ ,3369
+ ,3369
+ ,28
+ ,28
+ ,87161
+ ,87161
+ ,688779
+ ,1
+ ,32
+ ,32
+ ,11819
+ ,11819
+ ,85
+ ,85
+ ,106113
+ ,106113
+ ,574339
+ ,0
+ ,33
+ ,0
+ ,6984
+ ,0
+ ,45
+ ,0
+ ,80570
+ ,0
+ ,741409
+ ,1
+ ,34
+ ,34
+ ,4519
+ ,4519
+ ,49
+ ,49
+ ,102129
+ ,102129
+ ,597793
+ ,1
+ ,35
+ ,35
+ ,2220
+ ,2220
+ ,22
+ ,22
+ ,301670
+ ,301670
+ ,644190
+ ,1
+ ,36
+ ,36
+ ,18562
+ ,18562
+ ,155
+ ,155
+ ,102313
+ ,102313
+ ,377934
+ ,1
+ ,37
+ ,37
+ ,10327
+ ,10327
+ ,91
+ ,91
+ ,88577
+ ,88577
+ ,640273
+ ,1
+ ,38
+ ,38
+ ,5336
+ ,5336
+ ,81
+ ,81
+ ,112477
+ ,112477
+ ,697458
+ ,1
+ ,39
+ ,39
+ ,2365
+ ,2365
+ ,79
+ ,79
+ ,191778
+ ,191778
+ ,550608
+ ,1
+ ,40
+ ,40
+ ,4069
+ ,4069
+ ,145
+ ,145
+ ,79804
+ ,79804
+ ,207393
+ ,0
+ ,41
+ ,0
+ ,8636
+ ,0
+ ,855
+ ,0
+ ,128294
+ ,0
+ ,301607
+ ,1
+ ,42
+ ,42
+ ,13718
+ ,13718
+ ,61
+ ,61
+ ,96448
+ ,96448
+ ,345783
+ ,0
+ ,43
+ ,0
+ ,4525
+ ,0
+ ,226
+ ,0
+ ,93811
+ ,0
+ ,501749
+ ,0
+ ,44
+ ,0
+ ,6869
+ ,0
+ ,105
+ ,0
+ ,117520
+ ,0
+ ,379983
+ ,0
+ ,45
+ ,0
+ ,4628
+ ,0
+ ,62
+ ,0
+ ,69159
+ ,0
+ ,387475
+ ,1
+ ,46
+ ,46
+ ,3689
+ ,3689
+ ,25
+ ,25
+ ,101792
+ ,101792
+ ,377305
+ ,1
+ ,47
+ ,47
+ ,4891
+ ,4891
+ ,217
+ ,217
+ ,210568
+ ,210568
+ ,370837
+ ,1
+ ,48
+ ,48
+ ,7489
+ ,7489
+ ,322
+ ,322
+ ,136996
+ ,136996
+ ,430866
+ ,1
+ ,49
+ ,49
+ ,4901
+ ,4901
+ ,84
+ ,84
+ ,121920
+ ,121920
+ ,469107
+ ,1
+ ,50
+ ,50
+ ,2284
+ ,2284
+ ,33
+ ,33
+ ,76403
+ ,76403
+ ,194493
+ ,1
+ ,51
+ ,51
+ ,3160
+ ,3160
+ ,108
+ ,108
+ ,108094
+ ,108094
+ ,530670
+ ,1
+ ,52
+ ,52
+ ,4150
+ ,4150
+ ,150
+ ,150
+ ,134759
+ ,134759
+ ,518365
+ ,1
+ ,53
+ ,53
+ ,7285
+ ,7285
+ ,115
+ ,115
+ ,188873
+ ,188873
+ ,491303
+ ,0
+ ,54
+ ,0
+ ,1134
+ ,0
+ ,162
+ ,0
+ ,146216
+ ,0
+ ,527021
+ ,1
+ ,55
+ ,55
+ ,4658
+ ,4658
+ ,158
+ ,158
+ ,156608
+ ,156608
+ ,233773
+ ,0
+ ,56
+ ,0
+ ,2384
+ ,0
+ ,97
+ ,0
+ ,61348
+ ,0
+ ,405972
+ ,1
+ ,57
+ ,57
+ ,3748
+ ,3748
+ ,9
+ ,9
+ ,50350
+ ,50350
+ ,652925
+ ,0
+ ,58
+ ,0
+ ,5371
+ ,0
+ ,66
+ ,0
+ ,87720
+ ,0
+ ,446211
+ ,1
+ ,59
+ ,59
+ ,1285
+ ,1285
+ ,107
+ ,107
+ ,99489
+ ,99489
+ ,341340
+ ,1
+ ,60
+ ,60
+ ,9327
+ ,9327
+ ,101
+ ,101
+ ,87419
+ ,87419
+ ,387699
+ ,0
+ ,61
+ ,0
+ ,5565
+ ,0
+ ,47
+ ,0
+ ,94355
+ ,0
+ ,493408
+ ,0
+ ,62
+ ,0
+ ,1528
+ ,0
+ ,38
+ ,0
+ ,60326
+ ,0
+ ,146494
+ ,1
+ ,63
+ ,63
+ ,3122
+ ,3122
+ ,34
+ ,34
+ ,94670
+ ,94670
+ ,414462
+ ,0
+ ,64
+ ,0
+ ,7561
+ ,0
+ ,87
+ ,0
+ ,82425
+ ,0
+ ,364304
+ ,0
+ ,65
+ ,0
+ ,2675
+ ,0
+ ,79
+ ,0
+ ,59017
+ ,0
+ ,355178
+ ,1
+ ,66
+ ,66
+ ,13253
+ ,13253
+ ,947
+ ,947
+ ,90829
+ ,90829
+ ,357760
+ ,1
+ ,67
+ ,67
+ ,880
+ ,880
+ ,74
+ ,74
+ ,80791
+ ,80791
+ ,261216
+ ,1
+ ,68
+ ,68
+ ,2053
+ ,2053
+ ,53
+ ,53
+ ,100423
+ ,100423
+ ,397144
+ ,1
+ ,69
+ ,69
+ ,1424
+ ,1424
+ ,94
+ ,94
+ ,131116
+ ,131116
+ ,374943
+ ,1
+ ,70
+ ,70
+ ,4036
+ ,4036
+ ,63
+ ,63
+ ,100269
+ ,100269
+ ,424898
+ ,1
+ ,71
+ ,71
+ ,3045
+ ,3045
+ ,58
+ ,58
+ ,27330
+ ,27330
+ ,202055
+ ,0
+ ,72
+ ,0
+ ,5119
+ ,0
+ ,49
+ ,0
+ ,39039
+ ,0
+ ,378525
+ ,0
+ ,73
+ ,0
+ ,1431
+ ,0
+ ,34
+ ,0
+ ,106885
+ ,0
+ ,310768
+ ,1
+ ,74
+ ,74
+ ,554
+ ,554
+ ,11
+ ,11
+ ,79285
+ ,79285
+ ,325738
+ ,1
+ ,75
+ ,75
+ ,1975
+ ,1975
+ ,35
+ ,35
+ ,118881
+ ,118881
+ ,394510
+ ,1
+ ,76
+ ,76
+ ,1765
+ ,1765
+ ,20
+ ,20
+ ,77623
+ ,77623
+ ,247060
+ ,1
+ ,77
+ ,77
+ ,1012
+ ,1012
+ ,47
+ ,47
+ ,114768
+ ,114768
+ ,368078
+ ,1
+ ,78
+ ,78
+ ,810
+ ,810
+ ,43
+ ,43
+ ,74015
+ ,74015
+ ,236761
+ ,1
+ ,79
+ ,79
+ ,1280
+ ,1280
+ ,117
+ ,117
+ ,69465
+ ,69465
+ ,312378
+ ,1
+ ,80
+ ,80
+ ,666
+ ,666
+ ,171
+ ,171
+ ,117869
+ ,117869
+ ,339836
+ ,1
+ ,81
+ ,81
+ ,1380
+ ,1380
+ ,26
+ ,26
+ ,60982
+ ,60982
+ ,347385
+ ,0
+ ,82
+ ,0
+ ,4677
+ ,0
+ ,75
+ ,0
+ ,90131
+ ,0
+ ,426280
+ ,1
+ ,83
+ ,83
+ ,876
+ ,876
+ ,59
+ ,59
+ ,138971
+ ,138971
+ ,352850
+ ,0
+ ,84
+ ,0
+ ,814
+ ,0
+ ,18
+ ,0
+ ,39625
+ ,0
+ ,301881
+ ,1
+ ,85
+ ,85
+ ,514
+ ,514
+ ,15
+ ,15
+ ,102725
+ ,102725
+ ,377516
+ ,1
+ ,86
+ ,86
+ ,5692
+ ,5692
+ ,72
+ ,72
+ ,64239
+ ,64239
+ ,357312
+ ,1
+ ,87
+ ,87
+ ,3642
+ ,3642
+ ,86
+ ,86
+ ,90262
+ ,90262
+ ,458343
+ ,0
+ ,88
+ ,0
+ ,540
+ ,0
+ ,14
+ ,0
+ ,103960
+ ,0
+ ,354228
+ ,0
+ ,89
+ ,0
+ ,2099
+ ,0
+ ,64
+ ,0
+ ,106611
+ ,0
+ ,308636
+ ,0
+ ,90
+ ,0
+ ,567
+ ,0
+ ,11
+ ,0
+ ,103345
+ ,0
+ ,386212
+ ,0
+ ,91
+ ,0
+ ,2001
+ ,0
+ ,52
+ ,0
+ ,95551
+ ,0
+ ,393343
+ ,1
+ ,92
+ ,92
+ ,2949
+ ,2949
+ ,41
+ ,41
+ ,82903
+ ,82903
+ ,378509
+ ,1
+ ,93
+ ,93
+ ,2253
+ ,2253
+ ,99
+ ,99
+ ,63593
+ ,63593
+ ,452469
+ ,1
+ ,94
+ ,94
+ ,6533
+ ,6533
+ ,75
+ ,75
+ ,126910
+ ,126910
+ ,364839
+ ,1
+ ,95
+ ,95
+ ,1889
+ ,1889
+ ,45
+ ,45
+ ,37527
+ ,37527
+ ,358649
+ ,1
+ ,96
+ ,96
+ ,3055
+ ,3055
+ ,43
+ ,43
+ ,60247
+ ,60247
+ ,376641
+ ,1
+ ,97
+ ,97
+ ,272
+ ,272
+ ,8
+ ,8
+ ,112995
+ ,112995
+ ,429112
+ ,0
+ ,98
+ ,0
+ ,1414
+ ,0
+ ,198
+ ,0
+ ,70184
+ ,0
+ ,330546
+ ,0
+ ,99
+ ,0
+ ,2564
+ ,0
+ ,22
+ ,0
+ ,130140
+ ,0
+ ,403560
+ ,1
+ ,100
+ ,100
+ ,1383
+ ,1383
+ ,11
+ ,11
+ ,73221
+ ,73221
+ ,317892
+ ,1
+ ,101
+ ,101
+ ,1261
+ ,1261
+ ,33
+ ,33
+ ,76114
+ ,76114
+ ,307528
+ ,1
+ ,102
+ ,102
+ ,975
+ ,975
+ ,23
+ ,23
+ ,90534
+ ,90534
+ ,235133
+ ,1
+ ,103
+ ,103
+ ,3366
+ ,3366
+ ,80
+ ,80
+ ,108479
+ ,108479
+ ,299243
+ ,1
+ ,104
+ ,104
+ ,576
+ ,576
+ ,18
+ ,18
+ ,113761
+ ,113761
+ ,314073
+ ,1
+ ,105
+ ,105
+ ,1686
+ ,1686
+ ,40
+ ,40
+ ,68696
+ ,68696
+ ,368186
+ ,1
+ ,106
+ ,106
+ ,746
+ ,746
+ ,23
+ ,23
+ ,71561
+ ,71561
+ ,269661
+ ,1
+ ,107
+ ,107
+ ,3192
+ ,3192
+ ,60
+ ,60
+ ,59831
+ ,59831
+ ,125390
+ ,1
+ ,108
+ ,108
+ ,2045
+ ,2045
+ ,20
+ ,20
+ ,97890
+ ,97890
+ ,510834
+ ,0
+ ,109
+ ,0
+ ,5702
+ ,0
+ ,61
+ ,0
+ ,101481
+ ,0
+ ,321896
+ ,1
+ ,110
+ ,110
+ ,1932
+ ,1932
+ ,36
+ ,36
+ ,72954
+ ,72954
+ ,249898
+ ,1
+ ,111
+ ,111
+ ,936
+ ,936
+ ,30
+ ,30
+ ,67939
+ ,67939
+ ,408881
+ ,1
+ ,112
+ ,112
+ ,3437
+ ,3437
+ ,47
+ ,47
+ ,48022
+ ,48022
+ ,158492
+ ,0
+ ,113
+ ,0
+ ,5131
+ ,0
+ ,71
+ ,0
+ ,86111
+ ,0
+ ,292154
+ ,1
+ ,114
+ ,114
+ ,2397
+ ,2397
+ ,14
+ ,14
+ ,74020
+ ,74020
+ ,289513
+ ,1
+ ,115
+ ,115
+ ,1389
+ ,1389
+ ,9
+ ,9
+ ,57530
+ ,57530
+ ,378049
+ ,0
+ ,116
+ ,0
+ ,1503
+ ,0
+ ,39
+ ,0
+ ,56364
+ ,0
+ ,343466
+ ,0
+ ,117
+ ,0
+ ,402
+ ,0
+ ,26
+ ,0
+ ,84990
+ ,0
+ ,332743
+ ,0
+ ,118
+ ,0
+ ,2239
+ ,0
+ ,21
+ ,0
+ ,88590
+ ,0
+ ,442882
+ ,1
+ ,119
+ ,119
+ ,2234
+ ,2234
+ ,16
+ ,16
+ ,77200
+ ,77200
+ ,214215
+ ,1
+ ,120
+ ,120
+ ,837
+ ,837
+ ,69
+ ,69
+ ,61262
+ ,61262
+ ,315688
+ ,1
+ ,121
+ ,121
+ ,10579
+ ,10579
+ ,92
+ ,92
+ ,110309
+ ,110309
+ ,375195
+ ,0
+ ,122
+ ,0
+ ,875
+ ,0
+ ,14
+ ,0
+ ,67000
+ ,0
+ ,334280
+ ,1
+ ,123
+ ,123
+ ,1585
+ ,1585
+ ,107
+ ,107
+ ,93099
+ ,93099
+ ,355864
+ ,1
+ ,124
+ ,124
+ ,1659
+ ,1659
+ ,29
+ ,29
+ ,107577
+ ,107577
+ ,480382
+ ,1
+ ,125
+ ,125
+ ,2647
+ ,2647
+ ,37
+ ,37
+ ,62920
+ ,62920
+ ,353058
+ ,1
+ ,126
+ ,126
+ ,3294
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,75832
+ ,217193
+ ,0
+ ,127
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,128
+ ,0
+ ,94
+ ,0
+ ,7
+ ,0
+ ,60793
+ ,0
+ ,314533
+ ,0
+ ,129
+ ,0
+ ,422
+ ,0
+ ,28
+ ,0
+ ,57935
+ ,0
+ ,318056
+ ,1
+ ,130
+ ,130
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,131
+ ,0
+ ,34
+ ,0
+ ,8
+ ,0
+ ,60630
+ ,0
+ ,314353
+ ,1
+ ,132
+ ,132
+ ,1558
+ ,1558
+ ,63
+ ,63
+ ,55637
+ ,55637
+ ,369448
+ ,1
+ ,133
+ ,133
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,134
+ ,0
+ ,43
+ ,0
+ ,3
+ ,0
+ ,60887
+ ,0
+ ,312846
+ ,1
+ ,135
+ ,135
+ ,645
+ ,645
+ ,5
+ ,5
+ ,60720
+ ,60720
+ ,312075
+ ,0
+ ,136
+ ,0
+ ,316
+ ,0
+ ,9
+ ,0
+ ,60505
+ ,0
+ ,315009
+ ,1
+ ,137
+ ,137
+ ,115
+ ,115
+ ,13
+ ,13
+ ,60945
+ ,60945
+ ,318903
+ ,0
+ ,138
+ ,0
+ ,5
+ ,0
+ ,2
+ ,0
+ ,60720
+ ,0
+ ,314887
+ ,0
+ ,139
+ ,0
+ ,897
+ ,0
+ ,5
+ ,0
+ ,60720
+ ,0
+ ,314913
+ ,0
+ ,140
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,141
+ ,0
+ ,389
+ ,0
+ ,14
+ ,0
+ ,58990
+ ,0
+ ,325506
+ ,1
+ ,142
+ ,142
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,143
+ ,0
+ ,1002
+ ,0
+ ,15
+ ,0
+ ,56750
+ ,0
+ ,298568
+ ,0
+ ,144
+ ,0
+ ,36
+ ,0
+ ,3
+ ,0
+ ,60894
+ ,0
+ ,315834
+ ,1
+ ,145
+ ,145
+ ,460
+ ,460
+ ,15
+ ,15
+ ,63346
+ ,63346
+ ,329784
+ ,0
+ ,146
+ ,0
+ ,309
+ ,0
+ ,11
+ ,0
+ ,56535
+ ,0
+ ,312878
+ ,1
+ ,147
+ ,147
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,148
+ ,0
+ ,9
+ ,0
+ ,6
+ ,0
+ ,60835
+ ,0
+ ,314987
+ ,1
+ ,149
+ ,149
+ ,271
+ ,271
+ ,2
+ ,2
+ ,60720
+ ,60720
+ ,325249
+ ,1
+ ,150
+ ,150
+ ,14
+ ,14
+ ,1
+ ,1
+ ,61016
+ ,61016
+ ,315877
+ ,0
+ ,151
+ ,0
+ ,520
+ ,0
+ ,10
+ ,0
+ ,58650
+ ,0
+ ,291650
+ ,0
+ ,152
+ ,0
+ ,1766
+ ,0
+ ,73
+ ,0
+ ,60438
+ ,0
+ ,305959
+ ,1
+ ,153
+ ,153
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,154
+ ,154
+ ,458
+ ,458
+ ,11
+ ,11
+ ,58625
+ ,58625
+ ,297765
+ ,0
+ ,155
+ ,0
+ ,20
+ ,0
+ ,3
+ ,0
+ ,60938
+ ,0
+ ,315245
+ ,0
+ ,156
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,157
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,158
+ ,0
+ ,98
+ ,0
+ ,2
+ ,0
+ ,61490
+ ,0
+ ,315236
+ ,0
+ ,159
+ ,0
+ ,405
+ ,0
+ ,7
+ ,0
+ ,60845
+ ,0
+ ,336425
+ ,0
+ ,160
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,161
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,162
+ ,162
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,163
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,164
+ ,164
+ ,483
+ ,483
+ ,27
+ ,27
+ ,60830
+ ,60830
+ ,306268
+ ,1
+ ,165
+ ,165
+ ,454
+ ,454
+ ,51
+ ,51
+ ,63261
+ ,63261
+ ,302187
+ ,0
+ ,166
+ ,0
+ ,47
+ ,0
+ ,3
+ ,0
+ ,60720
+ ,0
+ ,314882
+ ,1
+ ,167
+ ,167
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,168
+ ,168
+ ,757
+ ,757
+ ,19
+ ,19
+ ,45689
+ ,45689
+ ,382712
+ ,1
+ ,169
+ ,169
+ ,4655
+ ,4655
+ ,393
+ ,393
+ ,60720
+ ,60720
+ ,341570
+ ,0
+ ,170
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,171
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,172
+ ,172
+ ,36
+ ,36
+ ,4
+ ,4
+ ,61564
+ ,61564
+ ,312412
+ ,1
+ ,173
+ ,173
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,174
+ ,174
+ ,203
+ ,203
+ ,9
+ ,9
+ ,61938
+ ,61938
+ ,309596
+ ,1
+ ,175
+ ,175
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,176
+ ,0
+ ,126
+ ,0
+ ,10
+ ,0
+ ,60951
+ ,0
+ ,315547
+ ,1
+ ,177
+ ,177
+ ,400
+ ,400
+ ,152
+ ,152
+ ,60720
+ ,60720
+ ,313267
+ ,1
+ ,178
+ ,178
+ ,71
+ ,71
+ ,1
+ ,1
+ ,60745
+ ,60745
+ ,316176
+ ,0
+ ,179
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,180
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,181
+ ,181
+ ,972
+ ,972
+ ,34
+ ,34
+ ,71642
+ ,71642
+ ,359335
+ ,0
+ ,182
+ ,0
+ ,531
+ ,0
+ ,10
+ ,0
+ ,71641
+ ,0
+ ,330068
+ ,0
+ ,183
+ ,0
+ ,2461
+ ,0
+ ,57
+ ,0
+ ,55792
+ ,0
+ ,314289
+ ,0
+ ,184
+ ,0
+ ,378
+ ,0
+ ,52
+ ,0
+ ,71873
+ ,0
+ ,297413
+ ,1
+ ,185
+ ,185
+ ,23
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,62555
+ ,314806
+ ,1
+ ,186
+ ,186
+ ,638
+ ,638
+ ,14
+ ,14
+ ,60370
+ ,60370
+ ,333210
+ ,0
+ ,187
+ ,0
+ ,2300
+ ,0
+ ,29
+ ,0
+ ,64873
+ ,0
+ ,352108
+ ,0
+ ,188
+ ,0
+ ,149
+ ,0
+ ,5
+ ,0
+ ,62041
+ ,0
+ ,313332
+ ,1
+ ,189
+ ,189
+ ,226
+ ,226
+ ,5
+ ,5
+ ,65745
+ ,65745
+ ,291787
+ ,0
+ ,190
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,191
+ ,191
+ ,275
+ ,275
+ ,4
+ ,4
+ ,59500
+ ,59500
+ ,318745
+ ,0
+ ,192
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,193
+ ,0
+ ,141
+ ,0
+ ,6
+ ,0
+ ,61630
+ ,0
+ ,315366
+ ,1
+ ,194
+ ,194
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,195
+ ,0
+ ,28
+ ,0
+ ,2
+ ,0
+ ,60890
+ ,0
+ ,315688
+ ,1
+ ,196
+ ,196
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,197
+ ,197
+ ,4980
+ ,4980
+ ,91
+ ,91
+ ,113521
+ ,113521
+ ,409642
+ ,1
+ ,198
+ ,198
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,199
+ ,199
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,200
+ ,200
+ ,472
+ ,472
+ ,20
+ ,20
+ ,80045
+ ,80045
+ ,269587
+ ,1
+ ,201
+ ,201
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,202
+ ,202
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,203
+ ,203
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,204
+ ,204
+ ,203
+ ,203
+ ,27
+ ,27
+ ,50804
+ ,50804
+ ,300962
+ ,1
+ ,205
+ ,205
+ ,496
+ ,496
+ ,17
+ ,17
+ ,87390
+ ,87390
+ ,325479
+ ,1
+ ,206
+ ,206
+ ,10
+ ,10
+ ,2
+ ,2
+ ,61656
+ ,61656
+ ,316155
+ ,1
+ ,207
+ ,207
+ ,63
+ ,63
+ ,4
+ ,4
+ ,65688
+ ,65688
+ ,318574
+ ,1
+ ,208
+ ,208
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,209
+ ,0
+ ,1136
+ ,0
+ ,32
+ ,0
+ ,48522
+ ,0
+ ,343613
+ ,0
+ ,210
+ ,0
+ ,265
+ ,0
+ ,31
+ ,0
+ ,60720
+ ,0
+ ,306948
+ ,0
+ ,211
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,212
+ ,212
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,213
+ ,213
+ ,267
+ ,267
+ ,32
+ ,32
+ ,57640
+ ,57640
+ ,330059
+ ,1
+ ,214
+ ,214
+ ,474
+ ,474
+ ,20
+ ,20
+ ,61977
+ ,61977
+ ,288985
+ ,1
+ ,215
+ ,215
+ ,534
+ ,534
+ ,7
+ ,7
+ ,62620
+ ,62620
+ ,304485
+ ,1
+ ,216
+ ,216
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,217
+ ,217
+ ,15
+ ,15
+ ,8
+ ,8
+ ,60831
+ ,60831
+ ,315688
+ ,0
+ ,218
+ ,0
+ ,397
+ ,0
+ ,28
+ ,0
+ ,60646
+ ,0
+ ,317736
+ ,0
+ ,219
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,220
+ ,0
+ ,1866
+ ,0
+ ,29
+ ,0
+ ,56225
+ ,0
+ ,322331
+ ,1
+ ,221
+ ,221
+ ,288
+ ,288
+ ,4
+ ,4
+ ,60510
+ ,60510
+ ,296656
+ ,1
+ ,222
+ ,222
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,223
+ ,0
+ ,3
+ ,0
+ ,2
+ ,0
+ ,60698
+ ,0
+ ,315354
+ ,1
+ ,224
+ ,224
+ ,468
+ ,468
+ ,21
+ ,21
+ ,60720
+ ,60720
+ ,312161
+ ,0
+ ,225
+ ,0
+ ,20
+ ,0
+ ,2
+ ,0
+ ,60805
+ ,0
+ ,315576
+ ,0
+ ,226
+ ,0
+ ,278
+ ,0
+ ,26
+ ,0
+ ,61404
+ ,0
+ ,314922
+ ,0
+ ,227
+ ,0
+ ,61
+ ,0
+ ,14
+ ,0
+ ,60720
+ ,0
+ ,314551
+ ,1
+ ,228
+ ,228
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,229
+ ,229
+ ,192
+ ,192
+ ,4
+ ,4
+ ,65276
+ ,65276
+ ,312339
+ ,0
+ ,230
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,231
+ ,231
+ ,317
+ ,317
+ ,9
+ ,9
+ ,63915
+ ,63915
+ ,298700
+ ,1
+ ,232
+ ,232
+ ,738
+ ,738
+ ,10
+ ,10
+ ,60720
+ ,60720
+ ,321376
+ ,1
+ ,233
+ ,233
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,234
+ ,0
+ ,368
+ ,0
+ ,17
+ ,0
+ ,61686
+ ,0
+ ,303230
+ ,0
+ ,235
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,236
+ ,0
+ ,2
+ ,0
+ ,1
+ ,0
+ ,60743
+ ,0
+ ,315487
+ ,0
+ ,237
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,238
+ ,0
+ ,53
+ ,0
+ ,6
+ ,0
+ ,60349
+ ,0
+ ,315793
+ ,1
+ ,239
+ ,239
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,240
+ ,240
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,241
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,242
+ ,242
+ ,94
+ ,94
+ ,3
+ ,3
+ ,61360
+ ,61360
+ ,312887
+ ,1
+ ,243
+ ,243
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,244
+ ,244
+ ,24
+ ,24
+ ,8
+ ,8
+ ,59818
+ ,59818
+ ,315637
+ ,1
+ ,245
+ ,245
+ ,2332
+ ,2332
+ ,4
+ ,4
+ ,72680
+ ,72680
+ ,324385
+ ,1
+ ,246
+ ,246
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,247
+ ,247
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,248
+ ,0
+ ,131
+ ,0
+ ,11
+ ,0
+ ,61808
+ ,0
+ ,308989
+ ,1
+ ,249
+ ,249
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,250
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,251
+ ,251
+ ,206
+ ,206
+ ,9
+ ,9
+ ,53110
+ ,53110
+ ,296702
+ ,1
+ ,252
+ ,252
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,253
+ ,253
+ ,167
+ ,167
+ ,2
+ ,2
+ ,64245
+ ,64245
+ ,307322
+ ,0
+ ,254
+ ,0
+ ,622
+ ,0
+ ,73
+ ,0
+ ,73007
+ ,0
+ ,304376
+ ,0
+ ,255
+ ,0
+ ,2328
+ ,0
+ ,94
+ ,0
+ ,82732
+ ,0
+ ,253588
+ ,0
+ ,256
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,257
+ ,0
+ ,365
+ ,0
+ ,8
+ ,0
+ ,54820
+ ,0
+ ,309560
+ ,0
+ ,258
+ ,0
+ ,364
+ ,0
+ ,35
+ ,0
+ ,47705
+ ,0
+ ,298466
+ ,1
+ ,259
+ ,259
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,260
+ ,260
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,261
+ ,261
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,262
+ ,262
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,263
+ ,0
+ ,226
+ ,0
+ ,12
+ ,0
+ ,72835
+ ,0
+ ,343929
+ ,1
+ ,264
+ ,264
+ ,307
+ ,307
+ ,15
+ ,15
+ ,58856
+ ,58856
+ ,331955
+ ,1
+ ,265
+ ,265
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,266
+ ,266
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,267
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,268
+ ,268
+ ,188
+ ,188
+ ,11
+ ,11
+ ,77655
+ ,77655
+ ,381180
+ ,1
+ ,269
+ ,269
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,270
+ ,270
+ ,138
+ ,138
+ ,6
+ ,6
+ ,69817
+ ,69817
+ ,331420
+ ,1
+ ,271
+ ,271
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,272
+ ,272
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,273
+ ,273
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,274
+ ,0
+ ,125
+ ,0
+ ,12
+ ,0
+ ,60798
+ ,0
+ ,310201
+ ,0
+ ,275
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,276
+ ,276
+ ,282
+ ,282
+ ,30
+ ,30
+ ,62452
+ ,62452
+ ,320016
+ ,0
+ ,277
+ ,0
+ ,335
+ ,0
+ ,33
+ ,0
+ ,64175
+ ,0
+ ,320398
+ ,1
+ ,278
+ ,278
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,279
+ ,279
+ ,1324
+ ,1324
+ ,117
+ ,117
+ ,67440
+ ,67440
+ ,291841
+ ,1
+ ,280
+ ,280
+ ,176
+ ,176
+ ,28
+ ,28
+ ,68136
+ ,68136
+ ,310670
+ ,1
+ ,281
+ ,281
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,282
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,283
+ ,283
+ ,249
+ ,249
+ ,72
+ ,72
+ ,56726
+ ,56726
+ ,313491
+ ,0
+ ,284
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,285
+ ,0
+ ,333
+ ,0
+ ,13
+ ,0
+ ,70811
+ ,0
+ ,331323
+ ,0
+ ,286
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,287
+ ,287
+ ,601
+ ,601
+ ,6
+ ,6
+ ,60720
+ ,60720
+ ,319210
+ ,0
+ ,288
+ ,0
+ ,30
+ ,0
+ ,4
+ ,0
+ ,62045
+ ,0
+ ,318098
+ ,0
+ ,289
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,290
+ ,0
+ ,249
+ ,0
+ ,62
+ ,0
+ ,54323
+ ,0
+ ,292754
+ ,1
+ ,291
+ ,291
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,292
+ ,0
+ ,165
+ ,0
+ ,24
+ ,0
+ ,62841
+ ,0
+ ,325176
+ ,1
+ ,293
+ ,293
+ ,453
+ ,453
+ ,21
+ ,21
+ ,81125
+ ,81125
+ ,365959
+ ,0
+ ,294
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,295
+ ,295
+ ,53
+ ,53
+ ,14
+ ,14
+ ,59506
+ ,59506
+ ,302409
+ ,1
+ ,296
+ ,296
+ ,382
+ ,382
+ ,21
+ ,21
+ ,59365
+ ,59365
+ ,340968
+ ,0
+ ,297
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,298
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,299
+ ,299
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,300
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,301
+ ,301
+ ,30
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,60798
+ ,313164
+ ,1
+ ,302
+ ,302
+ ,290
+ ,290
+ ,2
+ ,2
+ ,58790
+ ,58790
+ ,301164
+ ,1
+ ,303
+ ,303
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,304
+ ,304
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,305
+ ,0
+ ,366
+ ,0
+ ,53
+ ,0
+ ,61808
+ ,0
+ ,344425
+ ,0
+ ,306
+ ,0
+ ,2
+ ,0
+ ,9
+ ,0
+ ,60735
+ ,0
+ ,315394
+ ,1
+ ,307
+ ,307
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,308
+ ,308
+ ,209
+ ,209
+ ,13
+ ,13
+ ,64016
+ ,64016
+ ,316647
+ ,0
+ ,309
+ ,0
+ ,384
+ ,0
+ ,22
+ ,0
+ ,54683
+ ,0
+ ,309836
+ ,1
+ ,310
+ ,310
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,311
+ ,311
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,312
+ ,312
+ ,365
+ ,365
+ ,83
+ ,83
+ ,87192
+ ,87192
+ ,346611
+ ,1
+ ,313
+ ,313
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,314
+ ,314
+ ,49
+ ,49
+ ,8
+ ,8
+ ,64107
+ ,64107
+ ,322031
+ ,0
+ ,315
+ ,0
+ ,3
+ ,0
+ ,4
+ ,0
+ ,60761
+ ,0
+ ,315656
+ ,0
+ ,316
+ ,0
+ ,133
+ ,0
+ ,14
+ ,0
+ ,65990
+ ,0
+ ,339445
+ ,0
+ ,317
+ ,0
+ ,32
+ ,0
+ ,1
+ ,0
+ ,59988
+ ,0
+ ,314964
+ ,0
+ ,318
+ ,0
+ ,368
+ ,0
+ ,17
+ ,0
+ ,61167
+ ,0
+ ,297141
+ ,0
+ ,319
+ ,0
+ ,1
+ ,0
+ ,6
+ ,0
+ ,60719
+ ,0
+ ,315372
+ ,0
+ ,320
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,321
+ ,321
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,322
+ ,322
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,323
+ ,323
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,324
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,325
+ ,325
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,326
+ ,326
+ ,22
+ ,22
+ ,2
+ ,2
+ ,60722
+ ,60722
+ ,312502
+ ,0
+ ,327
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,328
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,329
+ ,329
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,330
+ ,330
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,331
+ ,331
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,332
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,333
+ ,333
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,334
+ ,334
+ ,96
+ ,96
+ ,5
+ ,5
+ ,60379
+ ,60379
+ ,313729
+ ,0
+ ,335
+ ,0
+ ,1
+ ,0
+ ,2
+ ,0
+ ,60727
+ ,0
+ ,315388
+ ,0
+ ,336
+ ,0
+ ,314
+ ,0
+ ,5
+ ,0
+ ,60720
+ ,0
+ ,315371
+ ,1
+ ,337
+ ,337
+ ,844
+ ,844
+ ,78
+ ,78
+ ,60925
+ ,60925
+ ,296139
+ ,0
+ ,338
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,339
+ ,0
+ ,26
+ ,0
+ ,1
+ ,0
+ ,60896
+ ,0
+ ,313880
+ ,1
+ ,340
+ ,340
+ ,125
+ ,125
+ ,13
+ ,13
+ ,59734
+ ,59734
+ ,317698
+ ,0
+ ,341
+ ,0
+ ,304
+ ,0
+ ,15
+ ,0
+ ,62969
+ ,0
+ ,295580
+ ,0
+ ,342
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,343
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,344
+ ,344
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,345
+ ,0
+ ,621
+ ,0
+ ,48
+ ,0
+ ,60720
+ ,0
+ ,308256
+ ,0
+ ,346
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,347
+ ,0
+ ,119
+ ,0
+ ,6
+ ,0
+ ,59118
+ ,0
+ ,303677
+ ,0
+ ,348
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,349
+ ,349
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,350
+ ,0
+ ,1595
+ ,0
+ ,17
+ ,0
+ ,60720
+ ,0
+ ,319369
+ ,1
+ ,351
+ ,351
+ ,312
+ ,312
+ ,14
+ ,14
+ ,58598
+ ,58598
+ ,318690
+ ,0
+ ,352
+ ,0
+ ,60
+ ,0
+ ,10
+ ,0
+ ,61124
+ ,0
+ ,314049
+ ,1
+ ,353
+ ,353
+ ,587
+ ,587
+ ,12
+ ,12
+ ,59595
+ ,59595
+ ,325699
+ ,0
+ ,354
+ ,0
+ ,135
+ ,0
+ ,2
+ ,0
+ ,62065
+ ,0
+ ,314210
+ ,0
+ ,355
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,356
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,357
+ ,0
+ ,514
+ ,0
+ ,52
+ ,0
+ ,78780
+ ,0
+ ,322378
+ ,0
+ ,358
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,359
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,360
+ ,360
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,361
+ ,0
+ ,1
+ ,0
+ ,4
+ ,0
+ ,60722
+ ,0
+ ,315398
+ ,1
+ ,362
+ ,362
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,363
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,364
+ ,364
+ ,1763
+ ,1763
+ ,24
+ ,24
+ ,61600
+ ,61600
+ ,308336
+ ,1
+ ,365
+ ,365
+ ,180
+ ,180
+ ,11
+ ,11
+ ,59635
+ ,59635
+ ,316386
+ ,1
+ ,366
+ ,366
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,367
+ ,367
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,368
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,369
+ ,369
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,370
+ ,0
+ ,218
+ ,0
+ ,21
+ ,0
+ ,60720
+ ,0
+ ,315553
+ ,1
+ ,371
+ ,371
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,372
+ ,372
+ ,448
+ ,448
+ ,40
+ ,40
+ ,59781
+ ,59781
+ ,323361
+ ,0
+ ,373
+ ,0
+ ,227
+ ,0
+ ,9
+ ,0
+ ,76644
+ ,0
+ ,336639
+ ,1
+ ,374
+ ,374
+ ,174
+ ,174
+ ,1
+ ,1
+ ,64820
+ ,64820
+ ,307424
+ ,1
+ ,375
+ ,375
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,376
+ ,376
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,377
+ ,377
+ ,121
+ ,121
+ ,24
+ ,24
+ ,56178
+ ,56178
+ ,295370
+ ,0
+ ,378
+ ,0
+ ,607
+ ,0
+ ,11
+ ,0
+ ,60436
+ ,0
+ ,322340
+ ,0
+ ,379
+ ,0
+ ,2212
+ ,0
+ ,14
+ ,0
+ ,60720
+ ,0
+ ,319864
+ ,1
+ ,380
+ ,380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,1
+ ,381
+ ,381
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,382
+ ,0
+ ,530
+ ,0
+ ,60
+ ,0
+ ,73433
+ ,0
+ ,317291
+ ,0
+ ,383
+ ,0
+ ,571
+ ,0
+ ,80
+ ,0
+ ,41477
+ ,0
+ ,280398
+ ,1
+ ,384
+ ,384
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,385
+ ,0
+ ,78
+ ,0
+ ,16
+ ,0
+ ,62700
+ ,0
+ ,317330
+ ,1
+ ,386
+ ,386
+ ,2489
+ ,2489
+ ,40
+ ,40
+ ,67804
+ ,67804
+ ,238125
+ ,0
+ ,387
+ ,0
+ ,131
+ ,0
+ ,6
+ ,0
+ ,59661
+ ,0
+ ,327071
+ ,1
+ ,388
+ ,388
+ ,923
+ ,923
+ ,8
+ ,8
+ ,58620
+ ,58620
+ ,309038
+ ,0
+ ,389
+ ,0
+ ,72
+ ,0
+ ,3
+ ,0
+ ,60398
+ ,0
+ ,314210
+ ,1
+ ,390
+ ,390
+ ,572
+ ,572
+ ,16
+ ,16
+ ,58580
+ ,58580
+ ,307930
+ ,1
+ ,391
+ ,391
+ ,397
+ ,397
+ ,10
+ ,10
+ ,62710
+ ,62710
+ ,322327
+ ,0
+ ,392
+ ,0
+ ,450
+ ,0
+ ,8
+ ,0
+ ,59325
+ ,0
+ ,292136
+ ,1
+ ,393
+ ,393
+ ,622
+ ,622
+ ,7
+ ,7
+ ,60950
+ ,60950
+ ,263276
+ ,0
+ ,394
+ ,0
+ ,694
+ ,0
+ ,8
+ ,0
+ ,68060
+ ,0
+ ,367655
+ ,0
+ ,395
+ ,0
+ ,3425
+ ,0
+ ,12
+ ,0
+ ,83620
+ ,0
+ ,283910
+ ,1
+ ,396
+ ,396
+ ,562
+ ,562
+ ,13
+ ,13
+ ,58456
+ ,58456
+ ,283587
+ ,1
+ ,397
+ ,397
+ ,4917
+ ,4917
+ ,42
+ ,42
+ ,52811
+ ,52811
+ ,243650
+ ,1
+ ,398
+ ,398
+ ,1442
+ ,1442
+ ,118
+ ,118
+ ,121173
+ ,121173
+ ,438493
+ ,1
+ ,399
+ ,399
+ ,529
+ ,529
+ ,9
+ ,9
+ ,63870
+ ,63870
+ ,296261
+ ,1
+ ,400
+ ,400
+ ,2126
+ ,2126
+ ,138
+ ,138
+ ,21001
+ ,21001
+ ,230621
+ ,0
+ ,401
+ ,0
+ ,1061
+ ,0
+ ,5
+ ,0
+ ,70415
+ ,0
+ ,304252
+ ,0
+ ,402
+ ,0
+ ,776
+ ,0
+ ,9
+ ,0
+ ,64230
+ ,0
+ ,333505
+ ,0
+ ,403
+ ,0
+ ,611
+ ,0
+ ,8
+ ,0
+ ,59190
+ ,0
+ ,296919
+ ,0
+ ,404
+ ,0
+ ,1526
+ ,0
+ ,25
+ ,0
+ ,69351
+ ,0
+ ,278990
+ ,0
+ ,405
+ ,0
+ ,592
+ ,0
+ ,7
+ ,0
+ ,64270
+ ,0
+ ,276898
+ ,0
+ ,406
+ ,0
+ ,1182
+ ,0
+ ,13
+ ,0
+ ,70694
+ ,0
+ ,327007
+ ,1
+ ,407
+ ,407
+ ,621
+ ,621
+ ,16
+ ,16
+ ,68005
+ ,68005
+ ,317046
+ ,1
+ ,408
+ ,408
+ ,989
+ ,989
+ ,11
+ ,11
+ ,58930
+ ,58930
+ ,304555
+ ,0
+ ,409
+ ,0
+ ,438
+ ,0
+ ,11
+ ,0
+ ,58320
+ ,0
+ ,298096
+ ,1
+ ,410
+ ,410
+ ,726
+ ,726
+ ,3
+ ,3
+ ,69980
+ ,69980
+ ,231861
+ ,1
+ ,411
+ ,411
+ ,1303
+ ,1303
+ ,61
+ ,61
+ ,69863
+ ,69863
+ ,309422
+ ,0
+ ,412
+ ,0
+ ,7419
+ ,0
+ ,29
+ ,0
+ ,63255
+ ,0
+ ,286963
+ ,1
+ ,413
+ ,413
+ ,1164
+ ,1164
+ ,17
+ ,17
+ ,57320
+ ,57320
+ ,269753
+ ,1
+ ,414
+ ,414
+ ,3310
+ ,3310
+ ,33
+ ,33
+ ,75230
+ ,75230
+ ,448243
+ ,1
+ ,415
+ ,415
+ ,1920
+ ,1920
+ ,15
+ ,15
+ ,79420
+ ,79420
+ ,165404
+ ,0
+ ,416
+ ,0
+ ,965
+ ,0
+ ,3
+ ,0
+ ,73490
+ ,0
+ ,204325
+ ,1
+ ,417
+ ,417
+ ,3256
+ ,3256
+ ,66
+ ,66
+ ,35250
+ ,35250
+ ,407159
+ ,0
+ ,418
+ ,0
+ ,1135
+ ,0
+ ,17
+ ,0
+ ,62285
+ ,0
+ ,290476
+ ,0
+ ,419
+ ,0
+ ,1270
+ ,0
+ ,26
+ ,0
+ ,69206
+ ,0
+ ,275311
+ ,0
+ ,420
+ ,0
+ ,661
+ ,0
+ ,3
+ ,0
+ ,65920
+ ,0
+ ,246541
+ ,0
+ ,421
+ ,0
+ ,1013
+ ,0
+ ,2
+ ,0
+ ,69770
+ ,0
+ ,253468
+ ,1
+ ,422
+ ,422
+ ,2844
+ ,2844
+ ,67
+ ,67
+ ,72683
+ ,72683
+ ,240897
+ ,1
+ ,423
+ ,423
+ ,11528
+ ,11528
+ ,70
+ ,70
+ ,-14545
+ ,-14545
+ ,-83265
+ ,1
+ ,424
+ ,424
+ ,6526
+ ,6526
+ ,26
+ ,26
+ ,55830
+ ,55830
+ ,-42143
+ ,1
+ ,425
+ ,425
+ ,2264
+ ,2264
+ ,24
+ ,24
+ ,55174
+ ,55174
+ ,272713
+ ,1
+ ,426
+ ,426
+ ,5109
+ ,5109
+ ,97
+ ,97
+ ,67038
+ ,67038
+ ,215362
+ ,1
+ ,427
+ ,427
+ ,3999
+ ,3999
+ ,30
+ ,30
+ ,51252
+ ,51252
+ ,42754
+ ,1
+ ,428
+ ,428
+ ,35624
+ ,35624
+ ,223
+ ,223
+ ,157278
+ ,157278
+ ,306275
+ ,1
+ ,429
+ ,429
+ ,9252
+ ,9252
+ ,48
+ ,48
+ ,79510
+ ,79510
+ ,253537
+ ,1
+ ,430
+ ,430
+ ,15236
+ ,15236
+ ,90
+ ,90
+ ,77440
+ ,77440
+ ,372631
+ ,1
+ ,431
+ ,431
+ ,18073
+ ,18073
+ ,180
+ ,180
+ ,27284
+ ,27284
+ ,-7170)
+ ,dim=c(10
+ ,431)
+ ,dimnames=list(c('pop'
+ ,'t'
+ ,'pop_t'
+ ,'Costs'
+ ,'Costs_p'
+ ,'Trades'
+ ,'Trades_p'
+ ,'Dividends'
+ ,'Dividends_p'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(10,431),dimnames=list(c('pop','t','pop_t','Costs','Costs_p','Trades','Trades_p','Dividends','Dividends_p','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 = '10'
> #'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
Wealth\r pop t pop_t Costs Costs_p Trades Trades_p Dividends Dividends_p
1 6282154 1 1 1 162556 162556 1081 1081 213118 213118
2 4321023 1 2 2 29790 29790 309 309 81767 81767
3 4111912 1 3 3 87550 87550 458 458 153198 153198
4 223193 1 4 4 84738 84738 588 588 -26007 -26007
5 1491348 1 5 5 54660 54660 302 302 126942 126942
6 1629616 0 6 0 42634 0 156 0 157214 0
7 1398893 1 7 7 40949 40949 481 481 129352 129352
8 1926517 0 8 0 45187 0 353 0 234817 0
9 983660 1 9 9 37704 37704 452 452 60448 60448
10 1443586 1 10 10 16275 16275 109 109 47818 47818
11 1073089 1 11 11 25830 25830 115 115 245546 245546
12 984885 1 12 12 12679 12679 110 110 48020 48020
13 1405225 1 13 13 18014 18014 239 239 -1710 -1710
14 227132 1 14 14 43556 43556 247 247 32648 32648
15 929118 1 15 15 24811 24811 505 505 95350 95350
16 1071292 0 16 0 6575 0 159 0 151352 0
17 638830 0 17 0 7123 0 109 0 288170 0
18 856956 1 18 18 21950 21950 519 519 114337 114337
19 992426 1 19 19 37597 37597 248 248 37884 37884
20 444477 1 20 20 17821 17821 373 373 122844 122844
21 857217 1 21 21 12988 12988 119 119 82340 82340
22 711969 0 22 0 22330 0 84 0 79801 0
23 702380 0 23 0 13326 0 102 0 165548 0
24 358589 1 24 24 16189 16189 295 295 116384 116384
25 297978 1 25 25 7146 7146 105 105 134028 134028
26 585715 1 26 26 15824 15824 64 64 63838 63838
27 657954 1 27 27 27664 27664 282 282 74996 74996
28 209458 1 28 28 11920 11920 182 182 31080 31080
29 786690 1 29 29 8568 8568 37 37 32168 32168
30 439798 1 30 30 14416 14416 361 361 49857 49857
31 688779 1 31 31 3369 3369 28 28 87161 87161
32 574339 1 32 32 11819 11819 85 85 106113 106113
33 741409 0 33 0 6984 0 45 0 80570 0
34 597793 1 34 34 4519 4519 49 49 102129 102129
35 644190 1 35 35 2220 2220 22 22 301670 301670
36 377934 1 36 36 18562 18562 155 155 102313 102313
37 640273 1 37 37 10327 10327 91 91 88577 88577
38 697458 1 38 38 5336 5336 81 81 112477 112477
39 550608 1 39 39 2365 2365 79 79 191778 191778
40 207393 1 40 40 4069 4069 145 145 79804 79804
41 301607 0 41 0 8636 0 855 0 128294 0
42 345783 1 42 42 13718 13718 61 61 96448 96448
43 501749 0 43 0 4525 0 226 0 93811 0
44 379983 0 44 0 6869 0 105 0 117520 0
45 387475 0 45 0 4628 0 62 0 69159 0
46 377305 1 46 46 3689 3689 25 25 101792 101792
47 370837 1 47 47 4891 4891 217 217 210568 210568
48 430866 1 48 48 7489 7489 322 322 136996 136996
49 469107 1 49 49 4901 4901 84 84 121920 121920
50 194493 1 50 50 2284 2284 33 33 76403 76403
51 530670 1 51 51 3160 3160 108 108 108094 108094
52 518365 1 52 52 4150 4150 150 150 134759 134759
53 491303 1 53 53 7285 7285 115 115 188873 188873
54 527021 0 54 0 1134 0 162 0 146216 0
55 233773 1 55 55 4658 4658 158 158 156608 156608
56 405972 0 56 0 2384 0 97 0 61348 0
57 652925 1 57 57 3748 3748 9 9 50350 50350
58 446211 0 58 0 5371 0 66 0 87720 0
59 341340 1 59 59 1285 1285 107 107 99489 99489
60 387699 1 60 60 9327 9327 101 101 87419 87419
61 493408 0 61 0 5565 0 47 0 94355 0
62 146494 0 62 0 1528 0 38 0 60326 0
63 414462 1 63 63 3122 3122 34 34 94670 94670
64 364304 0 64 0 7561 0 87 0 82425 0
65 355178 0 65 0 2675 0 79 0 59017 0
66 357760 1 66 66 13253 13253 947 947 90829 90829
67 261216 1 67 67 880 880 74 74 80791 80791
68 397144 1 68 68 2053 2053 53 53 100423 100423
69 374943 1 69 69 1424 1424 94 94 131116 131116
70 424898 1 70 70 4036 4036 63 63 100269 100269
71 202055 1 71 71 3045 3045 58 58 27330 27330
72 378525 0 72 0 5119 0 49 0 39039 0
73 310768 0 73 0 1431 0 34 0 106885 0
74 325738 1 74 74 554 554 11 11 79285 79285
75 394510 1 75 75 1975 1975 35 35 118881 118881
76 247060 1 76 76 1765 1765 20 20 77623 77623
77 368078 1 77 77 1012 1012 47 47 114768 114768
78 236761 1 78 78 810 810 43 43 74015 74015
79 312378 1 79 79 1280 1280 117 117 69465 69465
80 339836 1 80 80 666 666 171 171 117869 117869
81 347385 1 81 81 1380 1380 26 26 60982 60982
82 426280 0 82 0 4677 0 75 0 90131 0
83 352850 1 83 83 876 876 59 59 138971 138971
84 301881 0 84 0 814 0 18 0 39625 0
85 377516 1 85 85 514 514 15 15 102725 102725
86 357312 1 86 86 5692 5692 72 72 64239 64239
87 458343 1 87 87 3642 3642 86 86 90262 90262
88 354228 0 88 0 540 0 14 0 103960 0
89 308636 0 89 0 2099 0 64 0 106611 0
90 386212 0 90 0 567 0 11 0 103345 0
91 393343 0 91 0 2001 0 52 0 95551 0
92 378509 1 92 92 2949 2949 41 41 82903 82903
93 452469 1 93 93 2253 2253 99 99 63593 63593
94 364839 1 94 94 6533 6533 75 75 126910 126910
95 358649 1 95 95 1889 1889 45 45 37527 37527
96 376641 1 96 96 3055 3055 43 43 60247 60247
97 429112 1 97 97 272 272 8 8 112995 112995
98 330546 0 98 0 1414 0 198 0 70184 0
99 403560 0 99 0 2564 0 22 0 130140 0
100 317892 1 100 100 1383 1383 11 11 73221 73221
101 307528 1 101 101 1261 1261 33 33 76114 76114
102 235133 1 102 102 975 975 23 23 90534 90534
103 299243 1 103 103 3366 3366 80 80 108479 108479
104 314073 1 104 104 576 576 18 18 113761 113761
105 368186 1 105 105 1686 1686 40 40 68696 68696
106 269661 1 106 106 746 746 23 23 71561 71561
107 125390 1 107 107 3192 3192 60 60 59831 59831
108 510834 1 108 108 2045 2045 20 20 97890 97890
109 321896 0 109 0 5702 0 61 0 101481 0
110 249898 1 110 110 1932 1932 36 36 72954 72954
111 408881 1 111 111 936 936 30 30 67939 67939
112 158492 1 112 112 3437 3437 47 47 48022 48022
113 292154 0 113 0 5131 0 71 0 86111 0
114 289513 1 114 114 2397 2397 14 14 74020 74020
115 378049 1 115 115 1389 1389 9 9 57530 57530
116 343466 0 116 0 1503 0 39 0 56364 0
117 332743 0 117 0 402 0 26 0 84990 0
118 442882 0 118 0 2239 0 21 0 88590 0
119 214215 1 119 119 2234 2234 16 16 77200 77200
120 315688 1 120 120 837 837 69 69 61262 61262
121 375195 1 121 121 10579 10579 92 92 110309 110309
122 334280 0 122 0 875 0 14 0 67000 0
123 355864 1 123 123 1585 1585 107 107 93099 93099
124 480382 1 124 124 1659 1659 29 29 107577 107577
125 353058 1 125 125 2647 2647 37 37 62920 62920
126 217193 1 126 126 3294 3294 23 23 75832 75832
127 315380 0 127 0 0 0 0 0 60720 0
128 314533 0 128 0 94 0 7 0 60793 0
129 318056 0 129 0 422 0 28 0 57935 0
130 315380 1 130 130 0 0 0 0 60720 60720
131 314353 0 131 0 34 0 8 0 60630 0
132 369448 1 132 132 1558 1558 63 63 55637 55637
133 315380 1 133 133 0 0 0 0 60720 60720
134 312846 0 134 0 43 0 3 0 60887 0
135 312075 1 135 135 645 645 5 5 60720 60720
136 315009 0 136 0 316 0 9 0 60505 0
137 318903 1 137 137 115 115 13 13 60945 60945
138 314887 0 138 0 5 0 2 0 60720 0
139 314913 0 139 0 897 0 5 0 60720 0
140 315380 0 140 0 0 0 0 0 60720 0
141 325506 0 141 0 389 0 14 0 58990 0
142 315380 1 142 142 0 0 0 0 60720 60720
143 298568 0 143 0 1002 0 15 0 56750 0
144 315834 0 144 0 36 0 3 0 60894 0
145 329784 1 145 145 460 460 15 15 63346 63346
146 312878 0 146 0 309 0 11 0 56535 0
147 315380 1 147 147 0 0 0 0 60720 60720
148 314987 0 148 0 9 0 6 0 60835 0
149 325249 1 149 149 271 271 2 2 60720 60720
150 315877 1 150 150 14 14 1 1 61016 61016
151 291650 0 151 0 520 0 10 0 58650 0
152 305959 0 152 0 1766 0 73 0 60438 0
153 315380 1 153 153 0 0 0 0 60720 60720
154 297765 1 154 154 458 458 11 11 58625 58625
155 315245 0 155 0 20 0 3 0 60938 0
156 315380 0 156 0 0 0 0 0 60720 0
157 315380 0 157 0 0 0 0 0 60720 0
158 315236 0 158 0 98 0 2 0 61490 0
159 336425 0 159 0 405 0 7 0 60845 0
160 315380 0 160 0 0 0 0 0 60720 0
161 315380 0 161 0 0 0 0 0 60720 0
162 315380 1 162 162 0 0 0 0 60720 60720
163 315380 0 163 0 0 0 0 0 60720 0
164 306268 1 164 164 483 483 27 27 60830 60830
165 302187 1 165 165 454 454 51 51 63261 63261
166 314882 0 166 0 47 0 3 0 60720 0
167 315380 1 167 167 0 0 0 0 60720 60720
168 382712 1 168 168 757 757 19 19 45689 45689
169 341570 1 169 169 4655 4655 393 393 60720 60720
170 315380 0 170 0 0 0 0 0 60720 0
171 315380 0 171 0 0 0 0 0 60720 0
172 312412 1 172 172 36 36 4 4 61564 61564
173 315380 1 173 173 0 0 0 0 60720 60720
174 309596 1 174 174 203 203 9 9 61938 61938
175 315380 1 175 175 0 0 0 0 60720 60720
176 315547 0 176 0 126 0 10 0 60951 0
177 313267 1 177 177 400 400 152 152 60720 60720
178 316176 1 178 178 71 71 1 1 60745 60745
179 315380 0 179 0 0 0 0 0 60720 0
180 315380 0 180 0 0 0 0 0 60720 0
181 359335 1 181 181 972 972 34 34 71642 71642
182 330068 0 182 0 531 0 10 0 71641 0
183 314289 0 183 0 2461 0 57 0 55792 0
184 297413 0 184 0 378 0 52 0 71873 0
185 314806 1 185 185 23 23 5 5 62555 62555
186 333210 1 186 186 638 638 14 14 60370 60370
187 352108 0 187 0 2300 0 29 0 64873 0
188 313332 0 188 0 149 0 5 0 62041 0
189 291787 1 189 189 226 226 5 5 65745 65745
190 315380 0 190 0 0 0 0 0 60720 0
191 318745 1 191 191 275 275 4 4 59500 59500
192 315380 0 192 0 0 0 0 0 60720 0
193 315366 0 193 0 141 0 6 0 61630 0
194 315380 1 194 194 0 0 0 0 60720 60720
195 315688 0 195 0 28 0 2 0 60890 0
196 315380 1 196 196 0 0 0 0 60720 60720
197 409642 1 197 197 4980 4980 91 91 113521 113521
198 315380 1 198 198 0 0 0 0 60720 60720
199 315380 1 199 199 0 0 0 0 60720 60720
200 269587 1 200 200 472 472 20 20 80045 80045
201 315380 1 201 201 0 0 0 0 60720 60720
202 315380 1 202 202 0 0 0 0 60720 60720
203 315380 1 203 203 0 0 0 0 60720 60720
204 300962 1 204 204 203 203 27 27 50804 50804
205 325479 1 205 205 496 496 17 17 87390 87390
206 316155 1 206 206 10 10 2 2 61656 61656
207 318574 1 207 207 63 63 4 4 65688 65688
208 315380 1 208 208 0 0 0 0 60720 60720
209 343613 0 209 0 1136 0 32 0 48522 0
210 306948 0 210 0 265 0 31 0 60720 0
211 315380 0 211 0 0 0 0 0 60720 0
212 315380 1 212 212 0 0 0 0 60720 60720
213 330059 1 213 213 267 267 32 32 57640 57640
214 288985 1 214 214 474 474 20 20 61977 61977
215 304485 1 215 215 534 534 7 7 62620 62620
216 315380 1 216 216 0 0 0 0 60720 60720
217 315688 1 217 217 15 15 8 8 60831 60831
218 317736 0 218 0 397 0 28 0 60646 0
219 315380 0 219 0 0 0 0 0 60720 0
220 322331 0 220 0 1866 0 29 0 56225 0
221 296656 1 221 221 288 288 4 4 60510 60510
222 315380 1 222 222 0 0 0 0 60720 60720
223 315354 0 223 0 3 0 2 0 60698 0
224 312161 1 224 224 468 468 21 21 60720 60720
225 315576 0 225 0 20 0 2 0 60805 0
226 314922 0 226 0 278 0 26 0 61404 0
227 314551 0 227 0 61 0 14 0 60720 0
228 315380 1 228 228 0 0 0 0 60720 60720
229 312339 1 229 229 192 192 4 4 65276 65276
230 315380 0 230 0 0 0 0 0 60720 0
231 298700 1 231 231 317 317 9 9 63915 63915
232 321376 1 232 232 738 738 10 10 60720 60720
233 315380 1 233 233 0 0 0 0 60720 60720
234 303230 0 234 0 368 0 17 0 61686 0
235 315380 0 235 0 0 0 0 0 60720 0
236 315487 0 236 0 2 0 1 0 60743 0
237 315380 0 237 0 0 0 0 0 60720 0
238 315793 0 238 0 53 0 6 0 60349 0
239 315380 1 239 239 0 0 0 0 60720 60720
240 315380 1 240 240 0 0 0 0 60720 60720
241 315380 0 241 0 0 0 0 0 60720 0
242 312887 1 242 242 94 94 3 3 61360 61360
243 315380 1 243 243 0 0 0 0 60720 60720
244 315637 1 244 244 24 24 8 8 59818 59818
245 324385 1 245 245 2332 2332 4 4 72680 72680
246 315380 1 246 246 0 0 0 0 60720 60720
247 315380 1 247 247 0 0 0 0 60720 60720
248 308989 0 248 0 131 0 11 0 61808 0
249 315380 1 249 249 0 0 0 0 60720 60720
250 315380 0 250 0 0 0 0 0 60720 0
251 296702 1 251 251 206 206 9 9 53110 53110
252 315380 1 252 252 0 0 0 0 60720 60720
253 307322 1 253 253 167 167 2 2 64245 64245
254 304376 0 254 0 622 0 73 0 73007 0
255 253588 0 255 0 2328 0 94 0 82732 0
256 315380 0 256 0 0 0 0 0 60720 0
257 309560 0 257 0 365 0 8 0 54820 0
258 298466 0 258 0 364 0 35 0 47705 0
259 315380 1 259 259 0 0 0 0 60720 60720
260 315380 1 260 260 0 0 0 0 60720 60720
261 315380 1 261 261 0 0 0 0 60720 60720
262 315380 1 262 262 0 0 0 0 60720 60720
263 343929 0 263 0 226 0 12 0 72835 0
264 331955 1 264 264 307 307 15 15 58856 58856
265 315380 1 265 265 0 0 0 0 60720 60720
266 315380 1 266 266 0 0 0 0 60720 60720
267 315380 0 267 0 0 0 0 0 60720 0
268 381180 1 268 268 188 188 11 11 77655 77655
269 315380 1 269 269 0 0 0 0 60720 60720
270 331420 1 270 270 138 138 6 6 69817 69817
271 315380 1 271 271 0 0 0 0 60720 60720
272 315380 1 272 272 0 0 0 0 60720 60720
273 315380 1 273 273 0 0 0 0 60720 60720
274 310201 0 274 0 125 0 12 0 60798 0
275 315380 0 275 0 0 0 0 0 60720 0
276 320016 1 276 276 282 282 30 30 62452 62452
277 320398 0 277 0 335 0 33 0 64175 0
278 315380 1 278 278 0 0 0 0 60720 60720
279 291841 1 279 279 1324 1324 117 117 67440 67440
280 310670 1 280 280 176 176 28 28 68136 68136
281 315380 1 281 281 0 0 0 0 60720 60720
282 315380 0 282 0 0 0 0 0 60720 0
283 313491 1 283 283 249 249 72 72 56726 56726
284 315380 0 284 0 0 0 0 0 60720 0
285 331323 0 285 0 333 0 13 0 70811 0
286 315380 0 286 0 0 0 0 0 60720 0
287 319210 1 287 287 601 601 6 6 60720 60720
288 318098 0 288 0 30 0 4 0 62045 0
289 315380 0 289 0 0 0 0 0 60720 0
290 292754 0 290 0 249 0 62 0 54323 0
291 315380 1 291 291 0 0 0 0 60720 60720
292 325176 0 292 0 165 0 24 0 62841 0
293 365959 1 293 293 453 453 21 21 81125 81125
294 315380 0 294 0 0 0 0 0 60720 0
295 302409 1 295 295 53 53 14 14 59506 59506
296 340968 1 296 296 382 382 21 21 59365 59365
297 315380 0 297 0 0 0 0 0 60720 0
298 315380 0 298 0 0 0 0 0 60720 0
299 315380 1 299 299 0 0 0 0 60720 60720
300 315380 0 300 0 0 0 0 0 60720 0
301 313164 1 301 301 30 30 4 4 60798 60798
302 301164 1 302 302 290 290 2 2 58790 58790
303 315380 1 303 303 0 0 0 0 60720 60720
304 315380 1 304 304 0 0 0 0 60720 60720
305 344425 0 305 0 366 0 53 0 61808 0
306 315394 0 306 0 2 0 9 0 60735 0
307 315380 1 307 307 0 0 0 0 60720 60720
308 316647 1 308 308 209 209 13 13 64016 64016
309 309836 0 309 0 384 0 22 0 54683 0
310 315380 1 310 310 0 0 0 0 60720 60720
311 315380 1 311 311 0 0 0 0 60720 60720
312 346611 1 312 312 365 365 83 83 87192 87192
313 315380 1 313 313 0 0 0 0 60720 60720
314 322031 1 314 314 49 49 8 8 64107 64107
315 315656 0 315 0 3 0 4 0 60761 0
316 339445 0 316 0 133 0 14 0 65990 0
317 314964 0 317 0 32 0 1 0 59988 0
318 297141 0 318 0 368 0 17 0 61167 0
319 315372 0 319 0 1 0 6 0 60719 0
320 315380 0 320 0 0 0 0 0 60720 0
321 315380 1 321 321 0 0 0 0 60720 60720
322 315380 1 322 322 0 0 0 0 60720 60720
323 315380 1 323 323 0 0 0 0 60720 60720
324 315380 0 324 0 0 0 0 0 60720 0
325 315380 1 325 325 0 0 0 0 60720 60720
326 312502 1 326 326 22 22 2 2 60722 60722
327 315380 0 327 0 0 0 0 0 60720 0
328 315380 0 328 0 0 0 0 0 60720 0
329 315380 1 329 329 0 0 0 0 60720 60720
330 315380 1 330 330 0 0 0 0 60720 60720
331 315380 1 331 331 0 0 0 0 60720 60720
332 315380 0 332 0 0 0 0 0 60720 0
333 315380 1 333 333 0 0 0 0 60720 60720
334 313729 1 334 334 96 96 5 5 60379 60379
335 315388 0 335 0 1 0 2 0 60727 0
336 315371 0 336 0 314 0 5 0 60720 0
337 296139 1 337 337 844 844 78 78 60925 60925
338 315380 0 338 0 0 0 0 0 60720 0
339 313880 0 339 0 26 0 1 0 60896 0
340 317698 1 340 340 125 125 13 13 59734 59734
341 295580 0 341 0 304 0 15 0 62969 0
342 315380 0 342 0 0 0 0 0 60720 0
343 315380 0 343 0 0 0 0 0 60720 0
344 315380 1 344 344 0 0 0 0 60720 60720
345 308256 0 345 0 621 0 48 0 60720 0
346 315380 0 346 0 0 0 0 0 60720 0
347 303677 0 347 0 119 0 6 0 59118 0
348 315380 0 348 0 0 0 0 0 60720 0
349 315380 1 349 349 0 0 0 0 60720 60720
350 319369 0 350 0 1595 0 17 0 60720 0
351 318690 1 351 351 312 312 14 14 58598 58598
352 314049 0 352 0 60 0 10 0 61124 0
353 325699 1 353 353 587 587 12 12 59595 59595
354 314210 0 354 0 135 0 2 0 62065 0
355 315380 0 355 0 0 0 0 0 60720 0
356 315380 0 356 0 0 0 0 0 60720 0
357 322378 0 357 0 514 0 52 0 78780 0
358 315380 0 358 0 0 0 0 0 60720 0
359 315380 0 359 0 0 0 0 0 60720 0
360 315380 1 360 360 0 0 0 0 60720 60720
361 315398 0 361 0 1 0 4 0 60722 0
362 315380 1 362 362 0 0 0 0 60720 60720
363 315380 0 363 0 0 0 0 0 60720 0
364 308336 1 364 364 1763 1763 24 24 61600 61600
365 316386 1 365 365 180 180 11 11 59635 59635
366 315380 1 366 366 0 0 0 0 60720 60720
367 315380 1 367 367 0 0 0 0 60720 60720
368 315380 0 368 0 0 0 0 0 60720 0
369 315380 1 369 369 0 0 0 0 60720 60720
370 315553 0 370 0 218 0 21 0 60720 0
371 315380 1 371 371 0 0 0 0 60720 60720
372 323361 1 372 372 448 448 40 40 59781 59781
373 336639 0 373 0 227 0 9 0 76644 0
374 307424 1 374 374 174 174 1 1 64820 64820
375 315380 1 375 375 0 0 0 0 60720 60720
376 315380 1 376 376 0 0 0 0 60720 60720
377 295370 1 377 377 121 121 24 24 56178 56178
378 322340 0 378 0 607 0 11 0 60436 0
379 319864 0 379 0 2212 0 14 0 60720 0
380 315380 1 380 380 0 0 0 0 60720 60720
381 315380 1 381 381 0 0 0 0 60720 60720
382 317291 0 382 0 530 0 60 0 73433 0
383 280398 0 383 0 571 0 80 0 41477 0
384 315380 1 384 384 0 0 0 0 60720 60720
385 317330 0 385 0 78 0 16 0 62700 0
386 238125 1 386 386 2489 2489 40 40 67804 67804
387 327071 0 387 0 131 0 6 0 59661 0
388 309038 1 388 388 923 923 8 8 58620 58620
389 314210 0 389 0 72 0 3 0 60398 0
390 307930 1 390 390 572 572 16 16 58580 58580
391 322327 1 391 391 397 397 10 10 62710 62710
392 292136 0 392 0 450 0 8 0 59325 0
393 263276 1 393 393 622 622 7 7 60950 60950
394 367655 0 394 0 694 0 8 0 68060 0
395 283910 0 395 0 3425 0 12 0 83620 0
396 283587 1 396 396 562 562 13 13 58456 58456
397 243650 1 397 397 4917 4917 42 42 52811 52811
398 438493 1 398 398 1442 1442 118 118 121173 121173
399 296261 1 399 399 529 529 9 9 63870 63870
400 230621 1 400 400 2126 2126 138 138 21001 21001
401 304252 0 401 0 1061 0 5 0 70415 0
402 333505 0 402 0 776 0 9 0 64230 0
403 296919 0 403 0 611 0 8 0 59190 0
404 278990 0 404 0 1526 0 25 0 69351 0
405 276898 0 405 0 592 0 7 0 64270 0
406 327007 0 406 0 1182 0 13 0 70694 0
407 317046 1 407 407 621 621 16 16 68005 68005
408 304555 1 408 408 989 989 11 11 58930 58930
409 298096 0 409 0 438 0 11 0 58320 0
410 231861 1 410 410 726 726 3 3 69980 69980
411 309422 1 411 411 1303 1303 61 61 69863 69863
412 286963 0 412 0 7419 0 29 0 63255 0
413 269753 1 413 413 1164 1164 17 17 57320 57320
414 448243 1 414 414 3310 3310 33 33 75230 75230
415 165404 1 415 415 1920 1920 15 15 79420 79420
416 204325 0 416 0 965 0 3 0 73490 0
417 407159 1 417 417 3256 3256 66 66 35250 35250
418 290476 0 418 0 1135 0 17 0 62285 0
419 275311 0 419 0 1270 0 26 0 69206 0
420 246541 0 420 0 661 0 3 0 65920 0
421 253468 0 421 0 1013 0 2 0 69770 0
422 240897 1 422 422 2844 2844 67 67 72683 72683
423 -83265 1 423 423 11528 11528 70 70 -14545 -14545
424 -42143 1 424 424 6526 6526 26 26 55830 55830
425 272713 1 425 425 2264 2264 24 24 55174 55174
426 215362 1 426 426 5109 5109 97 97 67038 67038
427 42754 1 427 427 3999 3999 30 30 51252 51252
428 306275 1 428 428 35624 35624 223 223 157278 157278
429 253537 1 429 429 9252 9252 48 48 79510 79510
430 372631 1 430 430 15236 15236 90 90 77440 77440
431 -7170 1 431 431 18073 18073 180 180 27284 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) pop t pop_t Costs Costs_p
232191.937 -55988.497 -54.709 -107.741 28.039 2.577
Trades Trades_p Dividends Dividends_p
-293.801 -203.666 1.345 0.794
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2198518 -15868 17840 48514 3211947
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.322e+05 8.774e+04 2.646 0.00844 **
pop -5.599e+04 1.051e+05 -0.533 0.59447
t -5.471e+01 1.917e+02 -0.285 0.77545
pop_t -1.077e+02 2.394e+02 -0.450 0.65288
Costs 2.804e+01 4.991e+00 5.618 3.52e-08 ***
Costs_p 2.577e+00 5.307e+00 0.486 0.62756
Trades -2.938e+02 3.044e+02 -0.965 0.33507
Trades_p -2.037e+02 3.695e+02 -0.551 0.58175
Dividends 1.345e+00 9.529e-01 1.411 0.15892
Dividends_p 7.941e-01 1.081e+00 0.734 0.46308
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 254800 on 421 degrees of freedom
Multiple R-squared: 0.6569, Adjusted R-squared: 0.6496
F-statistic: 89.57 on 9 and 421 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 4.400530e-148 2.200265e-148
[2,] 1.0000000 1.554867e-150 7.774335e-151
[3,] 1.0000000 5.653359e-153 2.826679e-153
[4,] 1.0000000 3.682659e-160 1.841329e-160
[5,] 1.0000000 7.314090e-159 3.657045e-159
[6,] 1.0000000 3.987893e-162 1.993947e-162
[7,] 1.0000000 2.316144e-175 1.158072e-175
[8,] 1.0000000 4.204986e-174 2.102493e-174
[9,] 1.0000000 1.018454e-183 5.092269e-184
[10,] 1.0000000 1.414147e-182 7.070737e-183
[11,] 1.0000000 6.261065e-182 3.130533e-182
[12,] 1.0000000 4.609470e-181 2.304735e-181
[13,] 1.0000000 5.460259e-181 2.730130e-181
[14,] 1.0000000 5.537477e-184 2.768738e-184
[15,] 1.0000000 4.201303e-187 2.100652e-187
[16,] 1.0000000 9.633093e-188 4.816546e-188
[17,] 1.0000000 8.595219e-201 4.297609e-201
[18,] 1.0000000 1.325597e-201 6.627984e-202
[19,] 1.0000000 1.091894e-206 5.459470e-207
[20,] 1.0000000 9.257938e-208 4.628969e-208
[21,] 1.0000000 1.651817e-213 8.259084e-214
[22,] 1.0000000 1.266692e-215 6.333458e-216
[23,] 1.0000000 1.690792e-215 8.453960e-216
[24,] 1.0000000 1.835519e-214 9.177597e-215
[25,] 1.0000000 2.700805e-220 1.350402e-220
[26,] 1.0000000 4.275343e-227 2.137672e-227
[27,] 1.0000000 4.304887e-226 2.152444e-226
[28,] 1.0000000 1.804841e-226 9.024203e-227
[29,] 1.0000000 6.069717e-227 3.034858e-227
[30,] 1.0000000 6.948441e-226 3.474220e-226
[31,] 1.0000000 1.280941e-225 6.404706e-226
[32,] 1.0000000 1.644216e-224 8.221078e-225
[33,] 1.0000000 1.972109e-223 9.860544e-224
[34,] 1.0000000 7.881304e-223 3.940652e-223
[35,] 1.0000000 1.926405e-223 9.632025e-224
[36,] 1.0000000 9.108583e-223 4.554292e-223
[37,] 1.0000000 1.623387e-222 8.116934e-223
[38,] 1.0000000 2.092253e-223 1.046127e-223
[39,] 1.0000000 8.933483e-225 4.466742e-225
[40,] 1.0000000 7.344624e-225 3.672312e-225
[41,] 1.0000000 9.033955e-224 4.516978e-224
[42,] 1.0000000 4.869981e-224 2.434991e-224
[43,] 1.0000000 1.197466e-226 5.987329e-227
[44,] 1.0000000 5.379843e-226 2.689922e-226
[45,] 1.0000000 1.620947e-239 8.104737e-240
[46,] 1.0000000 5.148309e-239 2.574155e-239
[47,] 1.0000000 1.787517e-238 8.937585e-239
[48,] 1.0000000 4.016810e-238 2.008405e-238
[49,] 1.0000000 9.104649e-239 4.552324e-239
[50,] 1.0000000 3.118493e-241 1.559247e-241
[51,] 1.0000000 8.985215e-241 4.492607e-241
[52,] 1.0000000 1.243017e-239 6.215085e-240
[53,] 1.0000000 1.464565e-238 7.322825e-239
[54,] 1.0000000 7.056586e-238 3.528293e-238
[55,] 1.0000000 7.087356e-238 3.543678e-238
[56,] 1.0000000 5.078236e-237 2.539118e-237
[57,] 1.0000000 1.856625e-236 9.283126e-237
[58,] 1.0000000 7.048237e-236 3.524119e-236
[59,] 1.0000000 4.359799e-235 2.179900e-235
[60,] 1.0000000 2.045147e-234 1.022574e-234
[61,] 1.0000000 1.444045e-233 7.220226e-234
[62,] 1.0000000 1.090777e-232 5.453883e-233
[63,] 1.0000000 1.198003e-231 5.990016e-232
[64,] 1.0000000 1.768849e-231 8.844247e-232
[65,] 1.0000000 1.144539e-230 5.722695e-231
[66,] 1.0000000 7.054315e-231 3.527157e-231
[67,] 1.0000000 5.233201e-230 2.616601e-230
[68,] 1.0000000 3.401701e-230 1.700851e-230
[69,] 1.0000000 2.295402e-229 1.147701e-229
[70,] 1.0000000 3.386206e-229 1.693103e-229
[71,] 1.0000000 4.863370e-230 2.431685e-230
[72,] 1.0000000 5.958570e-229 2.979285e-229
[73,] 1.0000000 4.984891e-228 2.492446e-228
[74,] 1.0000000 1.664518e-227 8.322590e-228
[75,] 1.0000000 1.451721e-227 7.258606e-228
[76,] 1.0000000 1.784513e-226 8.922565e-227
[77,] 1.0000000 1.181634e-225 5.908168e-226
[78,] 1.0000000 1.290064e-224 6.450322e-225
[79,] 1.0000000 8.982328e-224 4.491164e-224
[80,] 1.0000000 7.429534e-223 3.714767e-223
[81,] 1.0000000 2.720259e-224 1.360130e-224
[82,] 1.0000000 1.711687e-223 8.558436e-224
[83,] 1.0000000 4.491791e-224 2.245895e-224
[84,] 1.0000000 4.670870e-224 2.335435e-224
[85,] 1.0000000 5.354134e-223 2.677067e-223
[86,] 1.0000000 5.639386e-222 2.819693e-222
[87,] 1.0000000 4.347163e-221 2.173582e-221
[88,] 1.0000000 4.987834e-220 2.493917e-220
[89,] 1.0000000 4.418621e-219 2.209311e-219
[90,] 1.0000000 1.305923e-220 6.529613e-221
[91,] 1.0000000 3.482411e-221 1.741206e-221
[92,] 1.0000000 8.992411e-223 4.496206e-223
[93,] 1.0000000 5.696748e-222 2.848374e-222
[94,] 1.0000000 1.301779e-221 6.508894e-222
[95,] 1.0000000 1.464940e-224 7.324698e-225
[96,] 1.0000000 2.842890e-225 1.421445e-225
[97,] 1.0000000 3.545371e-224 1.772685e-224
[98,] 1.0000000 2.911824e-224 1.455912e-224
[99,] 1.0000000 4.146804e-224 2.073402e-224
[100,] 1.0000000 1.161914e-224 5.809569e-225
[101,] 1.0000000 1.112810e-223 5.564050e-224
[102,] 1.0000000 6.474692e-223 3.237346e-223
[103,] 1.0000000 8.623666e-223 4.311833e-223
[104,] 1.0000000 7.845539e-222 3.922769e-222
[105,] 1.0000000 9.648979e-221 4.824490e-221
[106,] 1.0000000 5.078233e-222 2.539117e-222
[107,] 1.0000000 7.047939e-224 3.523970e-224
[108,] 1.0000000 7.812463e-223 3.906232e-223
[109,] 1.0000000 7.282684e-222 3.641342e-222
[110,] 1.0000000 8.076629e-221 4.038314e-221
[111,] 1.0000000 3.198762e-220 1.599381e-220
[112,] 1.0000000 1.733986e-219 8.669932e-220
[113,] 1.0000000 9.933919e-219 4.966960e-219
[114,] 1.0000000 3.993644e-220 1.996822e-220
[115,] 1.0000000 4.965811e-219 2.482906e-219
[116,] 1.0000000 6.197022e-218 3.098511e-218
[117,] 1.0000000 7.809642e-217 3.904821e-217
[118,] 1.0000000 8.564758e-216 4.282379e-216
[119,] 1.0000000 1.058332e-214 5.291658e-215
[120,] 1.0000000 2.064065e-214 1.032032e-214
[121,] 1.0000000 2.326691e-213 1.163346e-213
[122,] 1.0000000 2.829655e-212 1.414827e-212
[123,] 1.0000000 3.312155e-211 1.656077e-211
[124,] 1.0000000 4.072867e-210 2.036434e-210
[125,] 1.0000000 4.594916e-209 2.297458e-209
[126,] 1.0000000 5.535146e-208 2.767573e-208
[127,] 1.0000000 6.793192e-207 3.396596e-207
[128,] 1.0000000 8.120784e-206 4.060392e-206
[129,] 1.0000000 9.649778e-205 4.824889e-205
[130,] 1.0000000 1.062695e-203 5.313474e-204
[131,] 1.0000000 1.111056e-202 5.555278e-203
[132,] 1.0000000 1.307590e-201 6.537949e-202
[133,] 1.0000000 1.464795e-200 7.323973e-201
[134,] 1.0000000 1.686464e-199 8.432319e-200
[135,] 1.0000000 1.839419e-198 9.197094e-199
[136,] 1.0000000 2.113677e-197 1.056838e-197
[137,] 1.0000000 2.325656e-196 1.162828e-196
[138,] 1.0000000 2.526576e-195 1.263288e-195
[139,] 1.0000000 2.125288e-194 1.062644e-194
[140,] 1.0000000 2.432797e-193 1.216398e-193
[141,] 1.0000000 2.623405e-192 1.311702e-192
[142,] 1.0000000 2.637486e-191 1.318743e-191
[143,] 1.0000000 2.915979e-190 1.457989e-190
[144,] 1.0000000 3.200168e-189 1.600084e-189
[145,] 1.0000000 3.491286e-188 1.745643e-188
[146,] 1.0000000 3.802355e-187 1.901178e-187
[147,] 1.0000000 3.878080e-186 1.939040e-186
[148,] 1.0000000 4.177593e-185 2.088796e-185
[149,] 1.0000000 4.471484e-184 2.235742e-184
[150,] 1.0000000 4.657277e-183 2.328638e-183
[151,] 1.0000000 4.930019e-182 2.465009e-182
[152,] 1.0000000 4.814813e-181 2.407407e-181
[153,] 1.0000000 3.725081e-180 1.862541e-180
[154,] 1.0000000 3.887906e-179 1.943953e-179
[155,] 1.0000000 3.977637e-178 1.988819e-178
[156,] 1.0000000 6.042146e-179 3.021073e-179
[157,] 1.0000000 4.834982e-178 2.417491e-178
[158,] 1.0000000 5.031321e-177 2.515660e-177
[159,] 1.0000000 5.199840e-176 2.599920e-176
[160,] 1.0000000 5.293732e-175 2.646866e-175
[161,] 1.0000000 5.498719e-174 2.749360e-174
[162,] 1.0000000 5.446150e-173 2.723075e-173
[163,] 1.0000000 5.609030e-172 2.804515e-172
[164,] 1.0000000 5.724687e-171 2.862344e-171
[165,] 1.0000000 5.359759e-170 2.679880e-170
[166,] 1.0000000 5.477694e-169 2.738847e-169
[167,] 1.0000000 5.446859e-168 2.723430e-168
[168,] 1.0000000 5.378303e-167 2.689151e-167
[169,] 1.0000000 4.976371e-166 2.488186e-166
[170,] 1.0000000 4.838148e-165 2.419074e-165
[171,] 1.0000000 4.832162e-164 2.416081e-164
[172,] 1.0000000 4.370994e-163 2.185497e-163
[173,] 1.0000000 4.215330e-162 2.107665e-162
[174,] 1.0000000 3.624587e-161 1.812293e-161
[175,] 1.0000000 2.004847e-160 1.002424e-160
[176,] 1.0000000 1.932483e-159 9.662415e-160
[177,] 1.0000000 1.116701e-158 5.583506e-159
[178,] 1.0000000 1.069366e-157 5.346828e-158
[179,] 1.0000000 1.029583e-156 5.147915e-157
[180,] 1.0000000 9.752616e-156 4.876308e-156
[181,] 1.0000000 9.255062e-155 4.627531e-155
[182,] 1.0000000 8.836963e-154 4.418482e-154
[183,] 1.0000000 8.259820e-153 4.129910e-153
[184,] 1.0000000 7.815987e-152 3.907994e-152
[185,] 1.0000000 6.087461e-151 3.043731e-151
[186,] 1.0000000 5.701433e-150 2.850716e-150
[187,] 1.0000000 5.315938e-149 2.657969e-149
[188,] 1.0000000 1.445937e-149 7.229685e-150
[189,] 1.0000000 1.345736e-148 6.728682e-149
[190,] 1.0000000 1.246677e-147 6.233384e-148
[191,] 1.0000000 1.149511e-146 5.747555e-147
[192,] 1.0000000 1.047247e-145 5.236235e-146
[193,] 1.0000000 1.821450e-145 9.107252e-146
[194,] 1.0000000 1.637536e-144 8.187679e-145
[195,] 1.0000000 1.341044e-143 6.705221e-144
[196,] 1.0000000 1.205937e-142 6.029684e-143
[197,] 1.0000000 8.394403e-142 4.197202e-142
[198,] 1.0000000 7.422797e-141 3.711398e-141
[199,] 1.0000000 6.629118e-140 3.314559e-140
[200,] 1.0000000 5.871784e-139 2.935892e-139
[201,] 1.0000000 4.873052e-138 2.436526e-138
[202,] 1.0000000 2.873095e-137 1.436548e-137
[203,] 1.0000000 2.224464e-136 1.112232e-136
[204,] 1.0000000 1.928743e-135 9.643715e-136
[205,] 1.0000000 1.654778e-134 8.273890e-135
[206,] 1.0000000 1.458119e-133 7.290593e-134
[207,] 1.0000000 1.259735e-132 6.298676e-133
[208,] 1.0000000 1.038748e-131 5.193742e-132
[209,] 1.0000000 7.404112e-131 3.702056e-131
[210,] 1.0000000 6.233967e-130 3.116983e-130
[211,] 1.0000000 5.309976e-129 2.654988e-129
[212,] 1.0000000 4.391629e-128 2.195815e-128
[213,] 1.0000000 3.708265e-127 1.854132e-127
[214,] 1.0000000 3.147358e-126 1.573679e-126
[215,] 1.0000000 2.639053e-125 1.319526e-125
[216,] 1.0000000 2.166219e-124 1.083110e-124
[217,] 1.0000000 1.513178e-123 7.565891e-124
[218,] 1.0000000 1.244700e-122 6.223502e-123
[219,] 1.0000000 7.236334e-122 3.618167e-122
[220,] 1.0000000 5.985766e-121 2.992883e-121
[221,] 1.0000000 4.773860e-120 2.386930e-120
[222,] 1.0000000 3.664252e-119 1.832126e-119
[223,] 1.0000000 2.940787e-118 1.470393e-118
[224,] 1.0000000 2.348296e-117 1.174148e-117
[225,] 1.0000000 1.860665e-116 9.303323e-117
[226,] 1.0000000 1.475580e-115 7.377902e-116
[227,] 1.0000000 1.150838e-114 5.754188e-115
[228,] 1.0000000 8.927817e-114 4.463909e-114
[229,] 1.0000000 6.923057e-113 3.461528e-113
[230,] 1.0000000 5.187130e-112 2.593565e-112
[231,] 1.0000000 3.964653e-111 1.982326e-111
[232,] 1.0000000 3.046834e-110 1.523417e-110
[233,] 1.0000000 2.143365e-109 1.071683e-109
[234,] 1.0000000 1.610412e-108 8.052059e-109
[235,] 1.0000000 1.203078e-107 6.015388e-108
[236,] 1.0000000 8.833769e-107 4.416884e-107
[237,] 1.0000000 6.540705e-106 3.270352e-106
[238,] 1.0000000 4.858114e-105 2.429057e-105
[239,] 1.0000000 3.632745e-104 1.816372e-104
[240,] 1.0000000 2.652348e-103 1.326174e-103
[241,] 1.0000000 1.602253e-102 8.011267e-103
[242,] 1.0000000 1.187704e-101 5.938519e-102
[243,] 1.0000000 3.454700e-101 1.727350e-101
[244,] 1.0000000 2.484293e-100 1.242147e-100
[245,] 1.0000000 1.749354e-99 8.746769e-100
[246,] 1.0000000 1.106949e-98 5.534747e-99
[247,] 1.0000000 7.882784e-98 3.941392e-98
[248,] 1.0000000 5.581643e-97 2.790822e-97
[249,] 1.0000000 3.929493e-96 1.964746e-96
[250,] 1.0000000 2.750170e-95 1.375085e-95
[251,] 1.0000000 1.873074e-94 9.365370e-95
[252,] 1.0000000 1.284709e-93 6.423544e-94
[253,] 1.0000000 8.908050e-93 4.454025e-93
[254,] 1.0000000 6.140790e-92 3.070395e-92
[255,] 1.0000000 4.202071e-91 2.101035e-91
[256,] 1.0000000 2.896791e-90 1.448396e-90
[257,] 1.0000000 1.974691e-89 9.873456e-90
[258,] 1.0000000 1.248036e-88 6.240179e-89
[259,] 1.0000000 8.393463e-88 4.196731e-88
[260,] 1.0000000 5.609987e-87 2.804993e-87
[261,] 1.0000000 3.725768e-86 1.862884e-86
[262,] 1.0000000 2.398899e-85 1.199449e-85
[263,] 1.0000000 1.580199e-84 7.900994e-85
[264,] 1.0000000 1.033068e-83 5.165342e-84
[265,] 1.0000000 6.864920e-83 3.432460e-83
[266,] 1.0000000 4.457063e-82 2.228532e-82
[267,] 1.0000000 1.302027e-81 6.510135e-82
[268,] 1.0000000 5.528190e-81 2.764095e-81
[269,] 1.0000000 3.477925e-80 1.738963e-80
[270,] 1.0000000 2.237616e-79 1.118808e-79
[271,] 1.0000000 1.348296e-78 6.741482e-79
[272,] 1.0000000 8.583511e-78 4.291756e-78
[273,] 1.0000000 5.549383e-77 2.774692e-77
[274,] 1.0000000 3.497419e-76 1.748710e-76
[275,] 1.0000000 2.221971e-75 1.110986e-75
[276,] 1.0000000 1.399194e-74 6.995971e-75
[277,] 1.0000000 8.658805e-74 4.329403e-74
[278,] 1.0000000 4.201124e-73 2.100562e-73
[279,] 1.0000000 2.547711e-72 1.273856e-72
[280,] 1.0000000 1.588609e-71 7.943046e-72
[281,] 1.0000000 8.943066e-71 4.471533e-71
[282,] 1.0000000 5.379007e-70 2.689504e-70
[283,] 1.0000000 2.871047e-69 1.435524e-69
[284,] 1.0000000 1.685823e-68 8.429113e-69
[285,] 1.0000000 9.991909e-68 4.995955e-68
[286,] 1.0000000 5.872551e-67 2.936276e-67
[287,] 1.0000000 3.438723e-66 1.719362e-66
[288,] 1.0000000 1.995985e-65 9.979923e-66
[289,] 1.0000000 1.137124e-64 5.685620e-65
[290,] 1.0000000 6.253698e-64 3.126849e-64
[291,] 1.0000000 3.567994e-63 1.783997e-63
[292,] 1.0000000 2.020804e-62 1.010402e-62
[293,] 1.0000000 1.121866e-61 5.609330e-62
[294,] 1.0000000 6.357849e-61 3.178924e-61
[295,] 1.0000000 3.554129e-60 1.777065e-60
[296,] 1.0000000 1.843685e-59 9.218424e-60
[297,] 1.0000000 1.006452e-58 5.032260e-59
[298,] 1.0000000 5.513159e-58 2.756579e-58
[299,] 1.0000000 2.996146e-57 1.498073e-57
[300,] 1.0000000 4.539927e-57 2.269963e-57
[301,] 1.0000000 2.380789e-56 1.190394e-56
[302,] 1.0000000 1.167864e-55 5.839319e-56
[303,] 1.0000000 6.417653e-55 3.208827e-55
[304,] 1.0000000 3.443491e-54 1.721745e-54
[305,] 1.0000000 1.873301e-53 9.366504e-54
[306,] 1.0000000 8.896568e-53 4.448284e-53
[307,] 1.0000000 4.751388e-52 2.375694e-52
[308,] 1.0000000 2.520105e-51 1.260053e-51
[309,] 1.0000000 1.283463e-50 6.417315e-51
[310,] 1.0000000 6.472546e-50 3.236273e-50
[311,] 1.0000000 3.229464e-49 1.614732e-49
[312,] 1.0000000 1.678057e-48 8.390286e-49
[313,] 1.0000000 8.272975e-48 4.136488e-48
[314,] 1.0000000 3.921405e-47 1.960702e-47
[315,] 1.0000000 2.004913e-46 1.002456e-46
[316,] 1.0000000 1.016978e-45 5.084890e-46
[317,] 1.0000000 4.880059e-45 2.440029e-45
[318,] 1.0000000 2.311775e-44 1.155888e-44
[319,] 1.0000000 1.079161e-43 5.395804e-44
[320,] 1.0000000 5.359048e-43 2.679524e-43
[321,] 1.0000000 2.462904e-42 1.231452e-42
[322,] 1.0000000 1.095383e-41 5.476916e-42
[323,] 1.0000000 5.349603e-41 2.674801e-41
[324,] 1.0000000 2.594422e-40 1.297211e-40
[325,] 1.0000000 5.578702e-40 2.789351e-40
[326,] 1.0000000 2.694067e-39 1.347033e-39
[327,] 1.0000000 1.284442e-38 6.422212e-39
[328,] 1.0000000 5.392528e-38 2.696264e-38
[329,] 1.0000000 2.262606e-37 1.131303e-37
[330,] 1.0000000 1.058563e-36 5.292817e-37
[331,] 1.0000000 4.908415e-36 2.454207e-36
[332,] 1.0000000 2.048210e-35 1.024105e-35
[333,] 1.0000000 9.099905e-35 4.549952e-35
[334,] 1.0000000 4.130316e-34 2.065158e-34
[335,] 1.0000000 1.745749e-33 8.728747e-34
[336,] 1.0000000 7.739839e-33 3.869919e-33
[337,] 1.0000000 3.216749e-32 1.608374e-32
[338,] 1.0000000 1.424072e-31 7.120359e-32
[339,] 1.0000000 6.023781e-31 3.011891e-31
[340,] 1.0000000 2.582872e-30 1.291436e-30
[341,] 1.0000000 1.114017e-29 5.570084e-30
[342,] 1.0000000 4.698128e-29 2.349064e-29
[343,] 1.0000000 1.966209e-28 9.831043e-29
[344,] 1.0000000 8.112987e-28 4.056493e-28
[345,] 1.0000000 3.403691e-27 1.701845e-27
[346,] 1.0000000 1.373102e-26 6.865509e-27
[347,] 1.0000000 5.428195e-26 2.714098e-26
[348,] 1.0000000 2.184595e-25 1.092298e-25
[349,] 1.0000000 8.425168e-25 4.212584e-25
[350,] 1.0000000 3.351458e-24 1.675729e-24
[351,] 1.0000000 1.244699e-23 6.223496e-24
[352,] 1.0000000 4.751814e-23 2.375907e-23
[353,] 1.0000000 1.849582e-22 9.247909e-23
[354,] 1.0000000 7.136954e-22 3.568477e-22
[355,] 1.0000000 2.729400e-21 1.364700e-21
[356,] 1.0000000 9.703850e-21 4.851925e-21
[357,] 1.0000000 3.657971e-20 1.828985e-20
[358,] 1.0000000 1.303662e-19 6.518309e-20
[359,] 1.0000000 4.835018e-19 2.417509e-19
[360,] 1.0000000 1.757195e-18 8.785974e-19
[361,] 1.0000000 6.404980e-18 3.202490e-18
[362,] 1.0000000 2.193347e-17 1.096674e-17
[363,] 1.0000000 7.803584e-17 3.901792e-17
[364,] 1.0000000 2.746502e-16 1.373251e-16
[365,] 1.0000000 9.162375e-16 4.581188e-16
[366,] 1.0000000 3.114774e-15 1.557387e-15
[367,] 1.0000000 1.029861e-14 5.149303e-15
[368,] 1.0000000 3.461331e-14 1.730666e-14
[369,] 1.0000000 1.148074e-13 5.740370e-14
[370,] 1.0000000 3.701822e-13 1.850911e-13
[371,] 1.0000000 9.939450e-13 4.969725e-13
[372,] 1.0000000 3.169042e-12 1.584521e-12
[373,] 1.0000000 9.496433e-12 4.748217e-12
[374,] 1.0000000 1.666220e-11 8.331101e-12
[375,] 1.0000000 5.177557e-11 2.588779e-11
[376,] 1.0000000 1.590055e-10 7.950274e-11
[377,] 1.0000000 4.725275e-10 2.362638e-10
[378,] 1.0000000 1.409216e-09 7.046079e-10
[379,] 1.0000000 4.105441e-09 2.052721e-09
[380,] 1.0000000 1.015422e-08 5.077110e-09
[381,] 1.0000000 2.678623e-08 1.339311e-08
[382,] 1.0000000 6.914849e-08 3.457424e-08
[383,] 0.9999999 1.892890e-07 9.464452e-08
[384,] 0.9999997 5.101021e-07 2.550511e-07
[385,] 0.9999993 1.309857e-06 6.549287e-07
[386,] 0.9999989 2.295345e-06 1.147673e-06
[387,] 0.9999971 5.823078e-06 2.911539e-06
[388,] 0.9999937 1.262323e-05 6.311614e-06
[389,] 0.9999842 3.152539e-05 1.576269e-05
[390,] 0.9999629 7.419085e-05 3.709542e-05
[391,] 0.9999109 1.782377e-04 8.911887e-05
[392,] 0.9997980 4.040955e-04 2.020477e-04
[393,] 0.9995464 9.072898e-04 4.536449e-04
[394,] 0.9990635 1.873061e-03 9.365306e-04
[395,] 0.9979730 4.053961e-03 2.026980e-03
[396,] 0.9958077 8.384612e-03 4.192306e-03
[397,] 0.9915409 1.691811e-02 8.459057e-03
[398,] 0.9851035 2.979295e-02 1.489647e-02
[399,] 0.9735390 5.292200e-02 2.646100e-02
[400,] 0.9516090 9.678210e-02 4.839105e-02
[401,] 0.9148740 1.702520e-01 8.512599e-02
[402,] 0.8902781 2.194439e-01 1.097219e-01
[403,] 0.8595830 2.808340e-01 1.404170e-01
[404,] 0.7644545 4.710909e-01 2.355455e-01
[405,] 0.8256509 3.486982e-01 1.743491e-01
[406,] 0.6805720 6.388561e-01 3.194280e-01
> postscript(file="/var/www/html/rcomp/tmp/16rmo1291239858.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/26rmo1291239858.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/3yjmr1291239858.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/4yjmr1291239858.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/5yjmr1291239858.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
1.211341e+06 3.211947e+06 1.155998e+06 -2.198518e+06 -4.787482e+05
6 7 8 9 10
3.677997e+04 -6.721782e+04 2.157326e+05 -2.498352e+05 7.226928e+05
11 12 13 14 15
-3.600831e+05 3.744753e+05 8.021790e+05 -1.227235e+06 4.304057e+04
16 17 18 19 20
4.988105e+05 -1.476334e+05 2.531232e+04 -2.893929e+05 -3.512555e+05
21 22 23 24 25
1.698836e+05 -2.277514e+05 -9.484326e+04 -4.115152e+05 -3.273638e+05
26 27 28 29 30
-1.754206e+05 -3.809225e+05 -3.030669e+05 3.024907e+05 -9.993122e+04
31 32 33 34 35
2.419794e+05 -1.431770e+05 2.200782e+05 9.470666e+04 -2.285539e+05
36 37 38 39 40
-5.024220e+05 9.737859e+03 1.637956e+05 -6.253546e+04 -1.854367e+05
41 42 43 44 45
-9.180229e+04 -4.195152e+05 8.528477e+04 -1.695808e+05 -4.680120e+04
46 47 48 49 50
-1.096396e+05 -2.898785e+05 -9.963793e+04 -6.815482e+04 -1.905060e+05
51 52 53 54 55
8.854448e+04 9.955848e+03 -2.460721e+05 1.169643e+05 -3.324518e+05
56 57 58 59 60
5.600280e+04 2.680244e+05 -3.197073e+04 -2.417528e+04 -2.010330e+05
61 62 63 64 65
-4.553505e+03 -1.951059e+05 -3.265241e+04 -1.616639e+05 -4.611929e+03
66 67 68 69 70
6.337081e+04 -6.702604e+04 -1.928274e+04 -6.731344e+04 -4.661021e+04
71 72 73 74 75
-8.543808e+04 -3.135755e+04 -9.129420e+04 -1.950571e+04 -6.682351e+04
76 77 78 79 80
-1.269023e+05 -4.868155e+04 -8.848015e+04 1.945406e+04 -1.078963e+04
81 82 83 84 85
2.459810e+04 -3.172784e+04 -1.045663e+05 3.465232e+03 -1.285915e+04
86 87 88 89 90
-8.075890e+04 3.450340e+04 -2.397391e+04 -1.020981e+05 7.308063e+03
91 92 93 94 95
-3.186995e+03 -2.994864e+04 1.356350e+05 -2.302266e+05 8.217000e+04
96 97 98 99 100
1.503895e+04 2.264769e+04 2.786399e+04 -6.364448e+04 -3.553867e+04
101 102 103 104 105
-3.724832e+04 -1.365407e+05 -1.554944e+05 -9.722487e+04 3.039562e+04
106 107 108 109 110
-5.377295e+04 -2.292727e+05 9.015128e+04 -1.827500e+05 -1.057083e+05
111 112 113 114 115
9.167128e+04 -1.840696e+05 -1.726568e+05 -9.290390e+04 5.943598e+04
116 117 118 119 120
1.114294e+04 -1.096814e+04 4.140871e+04 -1.682057e+05 3.665309e+04
121 122 123 124 125
-2.953918e+05 -1.754003e+03 5.227681e+03 5.787494e+04 -4.364669e+01
126 127 128 129 130
-1.901348e+05 8.484910e+03 7.015435e+03 1.140951e+04 3.042868e+04
131 132 133 134 135
9.194867e+03 7.933457e+04 3.091603e+04 5.785050e+03 1.067631e+04
136 137 138 139 140
2.679422e+03 3.755391e+04 9.041119e+03 -1.500719e+04 9.196128e+03
141 142 143 144 145
1.490940e+04 3.237808e+04 -2.580087e+04 9.506997e+03 3.503192e+04
146 147 148 149 150
7.217910e+03 3.319033e+04 1.059662e+04 3.608238e+04 3.361046e+04
151 152 153 154 155
-2.279056e+04 -2.725782e+04 3.416503e+04 1.264349e+04 9.909246e+03
156 157 158 159 160
1.007147e+04 1.012618e+04 6.841279e+03 2.181349e+04 1.029031e+04
161 162 163 164 165
1.034502e+04 3.562709e+04 1.045443e+04 2.524909e+04 2.895824e+04
166 167 168 169 170
9.684151e+03 3.643934e+04 1.223577e+05 1.159441e+05 1.083740e+04
171 172 173 174 175
1.089211e+04 3.336618e+04 3.741404e+04 2.744973e+04 3.773894e+04
176 177 178 179 180
1.042717e+04 9.931973e+04 3.729259e+04 1.132978e+04 1.138449e+04
181 182 183 184 185
4.646463e+04 -4.542231e+02 -3.517192e+04 -1.668223e+04 3.664797e+04
186 187 188 189 190
4.553636e+04 -1.305767e+04 5.289045e+03 1.241147e+03 1.193158e+04
191 192 193 194 195
3.988307e+04 1.204100e+04 8.667379e+03 4.082549e+04 1.208704e+04
196 197 198 199 200
4.115040e+04 -8.454995e+04 4.147530e+04 4.163775e+04 -4.982577e+04
201 202 203 204 205
4.196265e+04 4.212510e+04 4.228755e+04 5.645678e+04 -1.105799e+04
206 207 208 209 210
4.223679e+04 3.556701e+04 4.309980e+04 3.515670e+04 6.271365e+03
211 212 213 214 215
1.308047e+04 4.374960e+04 7.292312e+04 1.042869e+04 1.641191e+04
216 217 218 219 220
4.439940e+04 4.815296e+04 1.301405e+04 1.351814e+04 -1.723142e+04
221 222 223 224 225
2.010942e+04 4.537410e+04 1.424405e+04 3.859880e+04 1.395492e+04
226 227 228 229 230
1.236742e+04 1.552967e+04 4.634880e+04 2.983770e+04 1.411994e+04
231 232 233 234 235
1.809488e+04 3.537509e+04 4.716106e+04 -4.433804e+03 1.439348e+04
236 237 238 239 240
1.476199e+04 1.450290e+04 1.574626e+04 4.813576e+04 4.829821e+04
241 242 243 244 245
1.472174e+04 4.337584e+04 4.878556e+04 5.437915e+04 -3.686962e+04
246 247 248 249 250
4.927291e+04 4.943536e+04 6.809405e+03 4.976026e+04 1.521412e+04
251 252 253 254 255
4.585367e+04 5.024761e+04 3.069504e+04 -8.086104e+03 -1.135608e+05
256 257 258 259 260
1.554237e+04 9.827235e+03 1.631627e+04 5.138476e+04 5.154721e+04
261 262 263 264 265
5.170966e+04 5.187211e+04 2.537199e+04 7.082176e+04 5.235946e+04
266 267 268 269 270
5.252191e+04 1.614417e+04 8.214310e+04 5.300927e+04 4.851515e+04
271 272 273 274 275
5.333417e+04 5.349662e+04 5.365907e+04 1.126404e+04 1.658184e+04
276 277 278 279 280
6.136853e+04 1.736577e+04 5.447132e+04 3.439106e+04 4.276584e+04
281 282 283 284 285
5.495867e+04 1.696481e+04 9.013125e+04 1.707422e+04 1.398496e+04
286 287 288 289 290
1.718364e+04 4.434830e+04 1.856336e+04 1.734777e+04 1.461269e+04
291 292 293 294 295
5.658317e+04 2.688061e+04 6.042339e+04 1.762131e+04 5.220037e+04
296 297 298 299 300
8.463318e+04 1.778544e+04 1.784015e+04 5.788277e+04 1.794957e+04
301 302 303 304 305
5.689626e+04 4.039842e+04 5.853258e+04 5.869503e+04 5.111439e+04
306 307 308 309 310
2.085978e+04 5.918238e+04 5.363088e+04 1.671281e+04 5.966973e+04
311 312 313 314 315
5.983218e+04 6.472314e+04 6.015708e+04 6.220608e+04 2.008216e+04
316 317 318 319 320
3.618734e+04 1.884452e+04 -5.229339e+03 2.071715e+04 1.904375e+04
321 322 323 324 325
6.145668e+04 6.161913e+04 6.178158e+04 1.926259e+04 6.210648e+04
326 327 328 329 330
5.970805e+04 1.942671e+04 1.948142e+04 6.275628e+04 6.291873e+04
331 332 333 334 335
6.308118e+04 1.970026e+04 6.340608e+04 6.219511e+04 2.042253e+04
336 337 338 339 340
1.257498e+04 5.733947e+04 2.002851e+04 1.791135e+04 7.161022e+04
341 342 343 344 345
-6.748344e+03 2.024735e+04 2.030206e+04 6.519304e+04 9.977962e+03
346 347 348 349 350
2.046618e+04 9.398343e+03 2.057560e+04 6.600529e+04 -1.505291e+04
351 352 353 354 355
7.159120e+04 2.017587e+04 6.737856e+04 1.472760e+04 2.095856e+04
356 357 358 359 360
2.101327e+04 4.646220e+03 2.112269e+04 2.117740e+04 6.779224e+04
361 362 363 364 365
2.244929e+04 6.811714e+04 2.139624e+04 1.748012e+04 7.189243e+04
366 367 368 369 370
6.876694e+04 6.892939e+04 2.166978e+04 6.925429e+04 2.200961e+04
371 372 373 374 375
6.957919e+04 8.591394e+04 1.806851e+04 4.851195e+04 7.022900e+04
376 377 378 379 380
7.039145e+04 6.849296e+04 1.577116e+04 -3.115255e+04 7.104125e+04
381 382 383 384 385
7.120370e+04 1.001894e+04 2.087885e+04 7.169105e+04 2.440110e+04
386 387 388 389 390
-7.669326e+04 3.391406e+04 4.621197e+04 2.094429e+04 6.024018e+04
391 392 393 394 395
6.833938e+04 -8.652277e+03 4.996653e+03 4.838863e+04 -1.316236e+05
396 397 398 399 400
3.595084e+04 -1.106540e+05 8.233669e+04 3.655330e+04 7.804335e+04
401 402 403 404 405
-2.896978e+04 1.782120e+04 -7.600154e+03 -5.979880e+04 -3.410396e+04
406 407 408 409 410
-7.358674e+03 5.045972e+04 4.378674e+04 8.070822e+02 -4.814371e+04
411 412 413 414 415
4.101799e+04 -2.072469e+05 1.086752e+04 9.347326e+04 -1.645637e+05
416 417 418 419 420
-1.301070e+05 1.564545e+05 -2.943243e+04 -5.499050e+04 -6.896898e+04
421 422 423 424 425
-7.732781e+04 -7.594499e+04 -4.777555e+05 -4.557375e+05 -9.827975e+03
426 427 428 429 430
-1.431767e+05 -2.812068e+05 -1.116491e+06 -2.824049e+05 -3.210303e+05
431
-6.354804e+05
> postscript(file="/var/www/html/rcomp/tmp/6rs3u1291239858.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 1.211341e+06 NA
1 3.211947e+06 1.211341e+06
2 1.155998e+06 3.211947e+06
3 -2.198518e+06 1.155998e+06
4 -4.787482e+05 -2.198518e+06
5 3.677997e+04 -4.787482e+05
6 -6.721782e+04 3.677997e+04
7 2.157326e+05 -6.721782e+04
8 -2.498352e+05 2.157326e+05
9 7.226928e+05 -2.498352e+05
10 -3.600831e+05 7.226928e+05
11 3.744753e+05 -3.600831e+05
12 8.021790e+05 3.744753e+05
13 -1.227235e+06 8.021790e+05
14 4.304057e+04 -1.227235e+06
15 4.988105e+05 4.304057e+04
16 -1.476334e+05 4.988105e+05
17 2.531232e+04 -1.476334e+05
18 -2.893929e+05 2.531232e+04
19 -3.512555e+05 -2.893929e+05
20 1.698836e+05 -3.512555e+05
21 -2.277514e+05 1.698836e+05
22 -9.484326e+04 -2.277514e+05
23 -4.115152e+05 -9.484326e+04
24 -3.273638e+05 -4.115152e+05
25 -1.754206e+05 -3.273638e+05
26 -3.809225e+05 -1.754206e+05
27 -3.030669e+05 -3.809225e+05
28 3.024907e+05 -3.030669e+05
29 -9.993122e+04 3.024907e+05
30 2.419794e+05 -9.993122e+04
31 -1.431770e+05 2.419794e+05
32 2.200782e+05 -1.431770e+05
33 9.470666e+04 2.200782e+05
34 -2.285539e+05 9.470666e+04
35 -5.024220e+05 -2.285539e+05
36 9.737859e+03 -5.024220e+05
37 1.637956e+05 9.737859e+03
38 -6.253546e+04 1.637956e+05
39 -1.854367e+05 -6.253546e+04
40 -9.180229e+04 -1.854367e+05
41 -4.195152e+05 -9.180229e+04
42 8.528477e+04 -4.195152e+05
43 -1.695808e+05 8.528477e+04
44 -4.680120e+04 -1.695808e+05
45 -1.096396e+05 -4.680120e+04
46 -2.898785e+05 -1.096396e+05
47 -9.963793e+04 -2.898785e+05
48 -6.815482e+04 -9.963793e+04
49 -1.905060e+05 -6.815482e+04
50 8.854448e+04 -1.905060e+05
51 9.955848e+03 8.854448e+04
52 -2.460721e+05 9.955848e+03
53 1.169643e+05 -2.460721e+05
54 -3.324518e+05 1.169643e+05
55 5.600280e+04 -3.324518e+05
56 2.680244e+05 5.600280e+04
57 -3.197073e+04 2.680244e+05
58 -2.417528e+04 -3.197073e+04
59 -2.010330e+05 -2.417528e+04
60 -4.553505e+03 -2.010330e+05
61 -1.951059e+05 -4.553505e+03
62 -3.265241e+04 -1.951059e+05
63 -1.616639e+05 -3.265241e+04
64 -4.611929e+03 -1.616639e+05
65 6.337081e+04 -4.611929e+03
66 -6.702604e+04 6.337081e+04
67 -1.928274e+04 -6.702604e+04
68 -6.731344e+04 -1.928274e+04
69 -4.661021e+04 -6.731344e+04
70 -8.543808e+04 -4.661021e+04
71 -3.135755e+04 -8.543808e+04
72 -9.129420e+04 -3.135755e+04
73 -1.950571e+04 -9.129420e+04
74 -6.682351e+04 -1.950571e+04
75 -1.269023e+05 -6.682351e+04
76 -4.868155e+04 -1.269023e+05
77 -8.848015e+04 -4.868155e+04
78 1.945406e+04 -8.848015e+04
79 -1.078963e+04 1.945406e+04
80 2.459810e+04 -1.078963e+04
81 -3.172784e+04 2.459810e+04
82 -1.045663e+05 -3.172784e+04
83 3.465232e+03 -1.045663e+05
84 -1.285915e+04 3.465232e+03
85 -8.075890e+04 -1.285915e+04
86 3.450340e+04 -8.075890e+04
87 -2.397391e+04 3.450340e+04
88 -1.020981e+05 -2.397391e+04
89 7.308063e+03 -1.020981e+05
90 -3.186995e+03 7.308063e+03
91 -2.994864e+04 -3.186995e+03
92 1.356350e+05 -2.994864e+04
93 -2.302266e+05 1.356350e+05
94 8.217000e+04 -2.302266e+05
95 1.503895e+04 8.217000e+04
96 2.264769e+04 1.503895e+04
97 2.786399e+04 2.264769e+04
98 -6.364448e+04 2.786399e+04
99 -3.553867e+04 -6.364448e+04
100 -3.724832e+04 -3.553867e+04
101 -1.365407e+05 -3.724832e+04
102 -1.554944e+05 -1.365407e+05
103 -9.722487e+04 -1.554944e+05
104 3.039562e+04 -9.722487e+04
105 -5.377295e+04 3.039562e+04
106 -2.292727e+05 -5.377295e+04
107 9.015128e+04 -2.292727e+05
108 -1.827500e+05 9.015128e+04
109 -1.057083e+05 -1.827500e+05
110 9.167128e+04 -1.057083e+05
111 -1.840696e+05 9.167128e+04
112 -1.726568e+05 -1.840696e+05
113 -9.290390e+04 -1.726568e+05
114 5.943598e+04 -9.290390e+04
115 1.114294e+04 5.943598e+04
116 -1.096814e+04 1.114294e+04
117 4.140871e+04 -1.096814e+04
118 -1.682057e+05 4.140871e+04
119 3.665309e+04 -1.682057e+05
120 -2.953918e+05 3.665309e+04
121 -1.754003e+03 -2.953918e+05
122 5.227681e+03 -1.754003e+03
123 5.787494e+04 5.227681e+03
124 -4.364669e+01 5.787494e+04
125 -1.901348e+05 -4.364669e+01
126 8.484910e+03 -1.901348e+05
127 7.015435e+03 8.484910e+03
128 1.140951e+04 7.015435e+03
129 3.042868e+04 1.140951e+04
130 9.194867e+03 3.042868e+04
131 7.933457e+04 9.194867e+03
132 3.091603e+04 7.933457e+04
133 5.785050e+03 3.091603e+04
134 1.067631e+04 5.785050e+03
135 2.679422e+03 1.067631e+04
136 3.755391e+04 2.679422e+03
137 9.041119e+03 3.755391e+04
138 -1.500719e+04 9.041119e+03
139 9.196128e+03 -1.500719e+04
140 1.490940e+04 9.196128e+03
141 3.237808e+04 1.490940e+04
142 -2.580087e+04 3.237808e+04
143 9.506997e+03 -2.580087e+04
144 3.503192e+04 9.506997e+03
145 7.217910e+03 3.503192e+04
146 3.319033e+04 7.217910e+03
147 1.059662e+04 3.319033e+04
148 3.608238e+04 1.059662e+04
149 3.361046e+04 3.608238e+04
150 -2.279056e+04 3.361046e+04
151 -2.725782e+04 -2.279056e+04
152 3.416503e+04 -2.725782e+04
153 1.264349e+04 3.416503e+04
154 9.909246e+03 1.264349e+04
155 1.007147e+04 9.909246e+03
156 1.012618e+04 1.007147e+04
157 6.841279e+03 1.012618e+04
158 2.181349e+04 6.841279e+03
159 1.029031e+04 2.181349e+04
160 1.034502e+04 1.029031e+04
161 3.562709e+04 1.034502e+04
162 1.045443e+04 3.562709e+04
163 2.524909e+04 1.045443e+04
164 2.895824e+04 2.524909e+04
165 9.684151e+03 2.895824e+04
166 3.643934e+04 9.684151e+03
167 1.223577e+05 3.643934e+04
168 1.159441e+05 1.223577e+05
169 1.083740e+04 1.159441e+05
170 1.089211e+04 1.083740e+04
171 3.336618e+04 1.089211e+04
172 3.741404e+04 3.336618e+04
173 2.744973e+04 3.741404e+04
174 3.773894e+04 2.744973e+04
175 1.042717e+04 3.773894e+04
176 9.931973e+04 1.042717e+04
177 3.729259e+04 9.931973e+04
178 1.132978e+04 3.729259e+04
179 1.138449e+04 1.132978e+04
180 4.646463e+04 1.138449e+04
181 -4.542231e+02 4.646463e+04
182 -3.517192e+04 -4.542231e+02
183 -1.668223e+04 -3.517192e+04
184 3.664797e+04 -1.668223e+04
185 4.553636e+04 3.664797e+04
186 -1.305767e+04 4.553636e+04
187 5.289045e+03 -1.305767e+04
188 1.241147e+03 5.289045e+03
189 1.193158e+04 1.241147e+03
190 3.988307e+04 1.193158e+04
191 1.204100e+04 3.988307e+04
192 8.667379e+03 1.204100e+04
193 4.082549e+04 8.667379e+03
194 1.208704e+04 4.082549e+04
195 4.115040e+04 1.208704e+04
196 -8.454995e+04 4.115040e+04
197 4.147530e+04 -8.454995e+04
198 4.163775e+04 4.147530e+04
199 -4.982577e+04 4.163775e+04
200 4.196265e+04 -4.982577e+04
201 4.212510e+04 4.196265e+04
202 4.228755e+04 4.212510e+04
203 5.645678e+04 4.228755e+04
204 -1.105799e+04 5.645678e+04
205 4.223679e+04 -1.105799e+04
206 3.556701e+04 4.223679e+04
207 4.309980e+04 3.556701e+04
208 3.515670e+04 4.309980e+04
209 6.271365e+03 3.515670e+04
210 1.308047e+04 6.271365e+03
211 4.374960e+04 1.308047e+04
212 7.292312e+04 4.374960e+04
213 1.042869e+04 7.292312e+04
214 1.641191e+04 1.042869e+04
215 4.439940e+04 1.641191e+04
216 4.815296e+04 4.439940e+04
217 1.301405e+04 4.815296e+04
218 1.351814e+04 1.301405e+04
219 -1.723142e+04 1.351814e+04
220 2.010942e+04 -1.723142e+04
221 4.537410e+04 2.010942e+04
222 1.424405e+04 4.537410e+04
223 3.859880e+04 1.424405e+04
224 1.395492e+04 3.859880e+04
225 1.236742e+04 1.395492e+04
226 1.552967e+04 1.236742e+04
227 4.634880e+04 1.552967e+04
228 2.983770e+04 4.634880e+04
229 1.411994e+04 2.983770e+04
230 1.809488e+04 1.411994e+04
231 3.537509e+04 1.809488e+04
232 4.716106e+04 3.537509e+04
233 -4.433804e+03 4.716106e+04
234 1.439348e+04 -4.433804e+03
235 1.476199e+04 1.439348e+04
236 1.450290e+04 1.476199e+04
237 1.574626e+04 1.450290e+04
238 4.813576e+04 1.574626e+04
239 4.829821e+04 4.813576e+04
240 1.472174e+04 4.829821e+04
241 4.337584e+04 1.472174e+04
242 4.878556e+04 4.337584e+04
243 5.437915e+04 4.878556e+04
244 -3.686962e+04 5.437915e+04
245 4.927291e+04 -3.686962e+04
246 4.943536e+04 4.927291e+04
247 6.809405e+03 4.943536e+04
248 4.976026e+04 6.809405e+03
249 1.521412e+04 4.976026e+04
250 4.585367e+04 1.521412e+04
251 5.024761e+04 4.585367e+04
252 3.069504e+04 5.024761e+04
253 -8.086104e+03 3.069504e+04
254 -1.135608e+05 -8.086104e+03
255 1.554237e+04 -1.135608e+05
256 9.827235e+03 1.554237e+04
257 1.631627e+04 9.827235e+03
258 5.138476e+04 1.631627e+04
259 5.154721e+04 5.138476e+04
260 5.170966e+04 5.154721e+04
261 5.187211e+04 5.170966e+04
262 2.537199e+04 5.187211e+04
263 7.082176e+04 2.537199e+04
264 5.235946e+04 7.082176e+04
265 5.252191e+04 5.235946e+04
266 1.614417e+04 5.252191e+04
267 8.214310e+04 1.614417e+04
268 5.300927e+04 8.214310e+04
269 4.851515e+04 5.300927e+04
270 5.333417e+04 4.851515e+04
271 5.349662e+04 5.333417e+04
272 5.365907e+04 5.349662e+04
273 1.126404e+04 5.365907e+04
274 1.658184e+04 1.126404e+04
275 6.136853e+04 1.658184e+04
276 1.736577e+04 6.136853e+04
277 5.447132e+04 1.736577e+04
278 3.439106e+04 5.447132e+04
279 4.276584e+04 3.439106e+04
280 5.495867e+04 4.276584e+04
281 1.696481e+04 5.495867e+04
282 9.013125e+04 1.696481e+04
283 1.707422e+04 9.013125e+04
284 1.398496e+04 1.707422e+04
285 1.718364e+04 1.398496e+04
286 4.434830e+04 1.718364e+04
287 1.856336e+04 4.434830e+04
288 1.734777e+04 1.856336e+04
289 1.461269e+04 1.734777e+04
290 5.658317e+04 1.461269e+04
291 2.688061e+04 5.658317e+04
292 6.042339e+04 2.688061e+04
293 1.762131e+04 6.042339e+04
294 5.220037e+04 1.762131e+04
295 8.463318e+04 5.220037e+04
296 1.778544e+04 8.463318e+04
297 1.784015e+04 1.778544e+04
298 5.788277e+04 1.784015e+04
299 1.794957e+04 5.788277e+04
300 5.689626e+04 1.794957e+04
301 4.039842e+04 5.689626e+04
302 5.853258e+04 4.039842e+04
303 5.869503e+04 5.853258e+04
304 5.111439e+04 5.869503e+04
305 2.085978e+04 5.111439e+04
306 5.918238e+04 2.085978e+04
307 5.363088e+04 5.918238e+04
308 1.671281e+04 5.363088e+04
309 5.966973e+04 1.671281e+04
310 5.983218e+04 5.966973e+04
311 6.472314e+04 5.983218e+04
312 6.015708e+04 6.472314e+04
313 6.220608e+04 6.015708e+04
314 2.008216e+04 6.220608e+04
315 3.618734e+04 2.008216e+04
316 1.884452e+04 3.618734e+04
317 -5.229339e+03 1.884452e+04
318 2.071715e+04 -5.229339e+03
319 1.904375e+04 2.071715e+04
320 6.145668e+04 1.904375e+04
321 6.161913e+04 6.145668e+04
322 6.178158e+04 6.161913e+04
323 1.926259e+04 6.178158e+04
324 6.210648e+04 1.926259e+04
325 5.970805e+04 6.210648e+04
326 1.942671e+04 5.970805e+04
327 1.948142e+04 1.942671e+04
328 6.275628e+04 1.948142e+04
329 6.291873e+04 6.275628e+04
330 6.308118e+04 6.291873e+04
331 1.970026e+04 6.308118e+04
332 6.340608e+04 1.970026e+04
333 6.219511e+04 6.340608e+04
334 2.042253e+04 6.219511e+04
335 1.257498e+04 2.042253e+04
336 5.733947e+04 1.257498e+04
337 2.002851e+04 5.733947e+04
338 1.791135e+04 2.002851e+04
339 7.161022e+04 1.791135e+04
340 -6.748344e+03 7.161022e+04
341 2.024735e+04 -6.748344e+03
342 2.030206e+04 2.024735e+04
343 6.519304e+04 2.030206e+04
344 9.977962e+03 6.519304e+04
345 2.046618e+04 9.977962e+03
346 9.398343e+03 2.046618e+04
347 2.057560e+04 9.398343e+03
348 6.600529e+04 2.057560e+04
349 -1.505291e+04 6.600529e+04
350 7.159120e+04 -1.505291e+04
351 2.017587e+04 7.159120e+04
352 6.737856e+04 2.017587e+04
353 1.472760e+04 6.737856e+04
354 2.095856e+04 1.472760e+04
355 2.101327e+04 2.095856e+04
356 4.646220e+03 2.101327e+04
357 2.112269e+04 4.646220e+03
358 2.117740e+04 2.112269e+04
359 6.779224e+04 2.117740e+04
360 2.244929e+04 6.779224e+04
361 6.811714e+04 2.244929e+04
362 2.139624e+04 6.811714e+04
363 1.748012e+04 2.139624e+04
364 7.189243e+04 1.748012e+04
365 6.876694e+04 7.189243e+04
366 6.892939e+04 6.876694e+04
367 2.166978e+04 6.892939e+04
368 6.925429e+04 2.166978e+04
369 2.200961e+04 6.925429e+04
370 6.957919e+04 2.200961e+04
371 8.591394e+04 6.957919e+04
372 1.806851e+04 8.591394e+04
373 4.851195e+04 1.806851e+04
374 7.022900e+04 4.851195e+04
375 7.039145e+04 7.022900e+04
376 6.849296e+04 7.039145e+04
377 1.577116e+04 6.849296e+04
378 -3.115255e+04 1.577116e+04
379 7.104125e+04 -3.115255e+04
380 7.120370e+04 7.104125e+04
381 1.001894e+04 7.120370e+04
382 2.087885e+04 1.001894e+04
383 7.169105e+04 2.087885e+04
384 2.440110e+04 7.169105e+04
385 -7.669326e+04 2.440110e+04
386 3.391406e+04 -7.669326e+04
387 4.621197e+04 3.391406e+04
388 2.094429e+04 4.621197e+04
389 6.024018e+04 2.094429e+04
390 6.833938e+04 6.024018e+04
391 -8.652277e+03 6.833938e+04
392 4.996653e+03 -8.652277e+03
393 4.838863e+04 4.996653e+03
394 -1.316236e+05 4.838863e+04
395 3.595084e+04 -1.316236e+05
396 -1.106540e+05 3.595084e+04
397 8.233669e+04 -1.106540e+05
398 3.655330e+04 8.233669e+04
399 7.804335e+04 3.655330e+04
400 -2.896978e+04 7.804335e+04
401 1.782120e+04 -2.896978e+04
402 -7.600154e+03 1.782120e+04
403 -5.979880e+04 -7.600154e+03
404 -3.410396e+04 -5.979880e+04
405 -7.358674e+03 -3.410396e+04
406 5.045972e+04 -7.358674e+03
407 4.378674e+04 5.045972e+04
408 8.070822e+02 4.378674e+04
409 -4.814371e+04 8.070822e+02
410 4.101799e+04 -4.814371e+04
411 -2.072469e+05 4.101799e+04
412 1.086752e+04 -2.072469e+05
413 9.347326e+04 1.086752e+04
414 -1.645637e+05 9.347326e+04
415 -1.301070e+05 -1.645637e+05
416 1.564545e+05 -1.301070e+05
417 -2.943243e+04 1.564545e+05
418 -5.499050e+04 -2.943243e+04
419 -6.896898e+04 -5.499050e+04
420 -7.732781e+04 -6.896898e+04
421 -7.594499e+04 -7.732781e+04
422 -4.777555e+05 -7.594499e+04
423 -4.557375e+05 -4.777555e+05
424 -9.827975e+03 -4.557375e+05
425 -1.431767e+05 -9.827975e+03
426 -2.812068e+05 -1.431767e+05
427 -1.116491e+06 -2.812068e+05
428 -2.824049e+05 -1.116491e+06
429 -3.210303e+05 -2.824049e+05
430 -6.354804e+05 -3.210303e+05
431 NA -6.354804e+05
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3.211947e+06 1.211341e+06
[2,] 1.155998e+06 3.211947e+06
[3,] -2.198518e+06 1.155998e+06
[4,] -4.787482e+05 -2.198518e+06
[5,] 3.677997e+04 -4.787482e+05
[6,] -6.721782e+04 3.677997e+04
[7,] 2.157326e+05 -6.721782e+04
[8,] -2.498352e+05 2.157326e+05
[9,] 7.226928e+05 -2.498352e+05
[10,] -3.600831e+05 7.226928e+05
[11,] 3.744753e+05 -3.600831e+05
[12,] 8.021790e+05 3.744753e+05
[13,] -1.227235e+06 8.021790e+05
[14,] 4.304057e+04 -1.227235e+06
[15,] 4.988105e+05 4.304057e+04
[16,] -1.476334e+05 4.988105e+05
[17,] 2.531232e+04 -1.476334e+05
[18,] -2.893929e+05 2.531232e+04
[19,] -3.512555e+05 -2.893929e+05
[20,] 1.698836e+05 -3.512555e+05
[21,] -2.277514e+05 1.698836e+05
[22,] -9.484326e+04 -2.277514e+05
[23,] -4.115152e+05 -9.484326e+04
[24,] -3.273638e+05 -4.115152e+05
[25,] -1.754206e+05 -3.273638e+05
[26,] -3.809225e+05 -1.754206e+05
[27,] -3.030669e+05 -3.809225e+05
[28,] 3.024907e+05 -3.030669e+05
[29,] -9.993122e+04 3.024907e+05
[30,] 2.419794e+05 -9.993122e+04
[31,] -1.431770e+05 2.419794e+05
[32,] 2.200782e+05 -1.431770e+05
[33,] 9.470666e+04 2.200782e+05
[34,] -2.285539e+05 9.470666e+04
[35,] -5.024220e+05 -2.285539e+05
[36,] 9.737859e+03 -5.024220e+05
[37,] 1.637956e+05 9.737859e+03
[38,] -6.253546e+04 1.637956e+05
[39,] -1.854367e+05 -6.253546e+04
[40,] -9.180229e+04 -1.854367e+05
[41,] -4.195152e+05 -9.180229e+04
[42,] 8.528477e+04 -4.195152e+05
[43,] -1.695808e+05 8.528477e+04
[44,] -4.680120e+04 -1.695808e+05
[45,] -1.096396e+05 -4.680120e+04
[46,] -2.898785e+05 -1.096396e+05
[47,] -9.963793e+04 -2.898785e+05
[48,] -6.815482e+04 -9.963793e+04
[49,] -1.905060e+05 -6.815482e+04
[50,] 8.854448e+04 -1.905060e+05
[51,] 9.955848e+03 8.854448e+04
[52,] -2.460721e+05 9.955848e+03
[53,] 1.169643e+05 -2.460721e+05
[54,] -3.324518e+05 1.169643e+05
[55,] 5.600280e+04 -3.324518e+05
[56,] 2.680244e+05 5.600280e+04
[57,] -3.197073e+04 2.680244e+05
[58,] -2.417528e+04 -3.197073e+04
[59,] -2.010330e+05 -2.417528e+04
[60,] -4.553505e+03 -2.010330e+05
[61,] -1.951059e+05 -4.553505e+03
[62,] -3.265241e+04 -1.951059e+05
[63,] -1.616639e+05 -3.265241e+04
[64,] -4.611929e+03 -1.616639e+05
[65,] 6.337081e+04 -4.611929e+03
[66,] -6.702604e+04 6.337081e+04
[67,] -1.928274e+04 -6.702604e+04
[68,] -6.731344e+04 -1.928274e+04
[69,] -4.661021e+04 -6.731344e+04
[70,] -8.543808e+04 -4.661021e+04
[71,] -3.135755e+04 -8.543808e+04
[72,] -9.129420e+04 -3.135755e+04
[73,] -1.950571e+04 -9.129420e+04
[74,] -6.682351e+04 -1.950571e+04
[75,] -1.269023e+05 -6.682351e+04
[76,] -4.868155e+04 -1.269023e+05
[77,] -8.848015e+04 -4.868155e+04
[78,] 1.945406e+04 -8.848015e+04
[79,] -1.078963e+04 1.945406e+04
[80,] 2.459810e+04 -1.078963e+04
[81,] -3.172784e+04 2.459810e+04
[82,] -1.045663e+05 -3.172784e+04
[83,] 3.465232e+03 -1.045663e+05
[84,] -1.285915e+04 3.465232e+03
[85,] -8.075890e+04 -1.285915e+04
[86,] 3.450340e+04 -8.075890e+04
[87,] -2.397391e+04 3.450340e+04
[88,] -1.020981e+05 -2.397391e+04
[89,] 7.308063e+03 -1.020981e+05
[90,] -3.186995e+03 7.308063e+03
[91,] -2.994864e+04 -3.186995e+03
[92,] 1.356350e+05 -2.994864e+04
[93,] -2.302266e+05 1.356350e+05
[94,] 8.217000e+04 -2.302266e+05
[95,] 1.503895e+04 8.217000e+04
[96,] 2.264769e+04 1.503895e+04
[97,] 2.786399e+04 2.264769e+04
[98,] -6.364448e+04 2.786399e+04
[99,] -3.553867e+04 -6.364448e+04
[100,] -3.724832e+04 -3.553867e+04
[101,] -1.365407e+05 -3.724832e+04
[102,] -1.554944e+05 -1.365407e+05
[103,] -9.722487e+04 -1.554944e+05
[104,] 3.039562e+04 -9.722487e+04
[105,] -5.377295e+04 3.039562e+04
[106,] -2.292727e+05 -5.377295e+04
[107,] 9.015128e+04 -2.292727e+05
[108,] -1.827500e+05 9.015128e+04
[109,] -1.057083e+05 -1.827500e+05
[110,] 9.167128e+04 -1.057083e+05
[111,] -1.840696e+05 9.167128e+04
[112,] -1.726568e+05 -1.840696e+05
[113,] -9.290390e+04 -1.726568e+05
[114,] 5.943598e+04 -9.290390e+04
[115,] 1.114294e+04 5.943598e+04
[116,] -1.096814e+04 1.114294e+04
[117,] 4.140871e+04 -1.096814e+04
[118,] -1.682057e+05 4.140871e+04
[119,] 3.665309e+04 -1.682057e+05
[120,] -2.953918e+05 3.665309e+04
[121,] -1.754003e+03 -2.953918e+05
[122,] 5.227681e+03 -1.754003e+03
[123,] 5.787494e+04 5.227681e+03
[124,] -4.364669e+01 5.787494e+04
[125,] -1.901348e+05 -4.364669e+01
[126,] 8.484910e+03 -1.901348e+05
[127,] 7.015435e+03 8.484910e+03
[128,] 1.140951e+04 7.015435e+03
[129,] 3.042868e+04 1.140951e+04
[130,] 9.194867e+03 3.042868e+04
[131,] 7.933457e+04 9.194867e+03
[132,] 3.091603e+04 7.933457e+04
[133,] 5.785050e+03 3.091603e+04
[134,] 1.067631e+04 5.785050e+03
[135,] 2.679422e+03 1.067631e+04
[136,] 3.755391e+04 2.679422e+03
[137,] 9.041119e+03 3.755391e+04
[138,] -1.500719e+04 9.041119e+03
[139,] 9.196128e+03 -1.500719e+04
[140,] 1.490940e+04 9.196128e+03
[141,] 3.237808e+04 1.490940e+04
[142,] -2.580087e+04 3.237808e+04
[143,] 9.506997e+03 -2.580087e+04
[144,] 3.503192e+04 9.506997e+03
[145,] 7.217910e+03 3.503192e+04
[146,] 3.319033e+04 7.217910e+03
[147,] 1.059662e+04 3.319033e+04
[148,] 3.608238e+04 1.059662e+04
[149,] 3.361046e+04 3.608238e+04
[150,] -2.279056e+04 3.361046e+04
[151,] -2.725782e+04 -2.279056e+04
[152,] 3.416503e+04 -2.725782e+04
[153,] 1.264349e+04 3.416503e+04
[154,] 9.909246e+03 1.264349e+04
[155,] 1.007147e+04 9.909246e+03
[156,] 1.012618e+04 1.007147e+04
[157,] 6.841279e+03 1.012618e+04
[158,] 2.181349e+04 6.841279e+03
[159,] 1.029031e+04 2.181349e+04
[160,] 1.034502e+04 1.029031e+04
[161,] 3.562709e+04 1.034502e+04
[162,] 1.045443e+04 3.562709e+04
[163,] 2.524909e+04 1.045443e+04
[164,] 2.895824e+04 2.524909e+04
[165,] 9.684151e+03 2.895824e+04
[166,] 3.643934e+04 9.684151e+03
[167,] 1.223577e+05 3.643934e+04
[168,] 1.159441e+05 1.223577e+05
[169,] 1.083740e+04 1.159441e+05
[170,] 1.089211e+04 1.083740e+04
[171,] 3.336618e+04 1.089211e+04
[172,] 3.741404e+04 3.336618e+04
[173,] 2.744973e+04 3.741404e+04
[174,] 3.773894e+04 2.744973e+04
[175,] 1.042717e+04 3.773894e+04
[176,] 9.931973e+04 1.042717e+04
[177,] 3.729259e+04 9.931973e+04
[178,] 1.132978e+04 3.729259e+04
[179,] 1.138449e+04 1.132978e+04
[180,] 4.646463e+04 1.138449e+04
[181,] -4.542231e+02 4.646463e+04
[182,] -3.517192e+04 -4.542231e+02
[183,] -1.668223e+04 -3.517192e+04
[184,] 3.664797e+04 -1.668223e+04
[185,] 4.553636e+04 3.664797e+04
[186,] -1.305767e+04 4.553636e+04
[187,] 5.289045e+03 -1.305767e+04
[188,] 1.241147e+03 5.289045e+03
[189,] 1.193158e+04 1.241147e+03
[190,] 3.988307e+04 1.193158e+04
[191,] 1.204100e+04 3.988307e+04
[192,] 8.667379e+03 1.204100e+04
[193,] 4.082549e+04 8.667379e+03
[194,] 1.208704e+04 4.082549e+04
[195,] 4.115040e+04 1.208704e+04
[196,] -8.454995e+04 4.115040e+04
[197,] 4.147530e+04 -8.454995e+04
[198,] 4.163775e+04 4.147530e+04
[199,] -4.982577e+04 4.163775e+04
[200,] 4.196265e+04 -4.982577e+04
[201,] 4.212510e+04 4.196265e+04
[202,] 4.228755e+04 4.212510e+04
[203,] 5.645678e+04 4.228755e+04
[204,] -1.105799e+04 5.645678e+04
[205,] 4.223679e+04 -1.105799e+04
[206,] 3.556701e+04 4.223679e+04
[207,] 4.309980e+04 3.556701e+04
[208,] 3.515670e+04 4.309980e+04
[209,] 6.271365e+03 3.515670e+04
[210,] 1.308047e+04 6.271365e+03
[211,] 4.374960e+04 1.308047e+04
[212,] 7.292312e+04 4.374960e+04
[213,] 1.042869e+04 7.292312e+04
[214,] 1.641191e+04 1.042869e+04
[215,] 4.439940e+04 1.641191e+04
[216,] 4.815296e+04 4.439940e+04
[217,] 1.301405e+04 4.815296e+04
[218,] 1.351814e+04 1.301405e+04
[219,] -1.723142e+04 1.351814e+04
[220,] 2.010942e+04 -1.723142e+04
[221,] 4.537410e+04 2.010942e+04
[222,] 1.424405e+04 4.537410e+04
[223,] 3.859880e+04 1.424405e+04
[224,] 1.395492e+04 3.859880e+04
[225,] 1.236742e+04 1.395492e+04
[226,] 1.552967e+04 1.236742e+04
[227,] 4.634880e+04 1.552967e+04
[228,] 2.983770e+04 4.634880e+04
[229,] 1.411994e+04 2.983770e+04
[230,] 1.809488e+04 1.411994e+04
[231,] 3.537509e+04 1.809488e+04
[232,] 4.716106e+04 3.537509e+04
[233,] -4.433804e+03 4.716106e+04
[234,] 1.439348e+04 -4.433804e+03
[235,] 1.476199e+04 1.439348e+04
[236,] 1.450290e+04 1.476199e+04
[237,] 1.574626e+04 1.450290e+04
[238,] 4.813576e+04 1.574626e+04
[239,] 4.829821e+04 4.813576e+04
[240,] 1.472174e+04 4.829821e+04
[241,] 4.337584e+04 1.472174e+04
[242,] 4.878556e+04 4.337584e+04
[243,] 5.437915e+04 4.878556e+04
[244,] -3.686962e+04 5.437915e+04
[245,] 4.927291e+04 -3.686962e+04
[246,] 4.943536e+04 4.927291e+04
[247,] 6.809405e+03 4.943536e+04
[248,] 4.976026e+04 6.809405e+03
[249,] 1.521412e+04 4.976026e+04
[250,] 4.585367e+04 1.521412e+04
[251,] 5.024761e+04 4.585367e+04
[252,] 3.069504e+04 5.024761e+04
[253,] -8.086104e+03 3.069504e+04
[254,] -1.135608e+05 -8.086104e+03
[255,] 1.554237e+04 -1.135608e+05
[256,] 9.827235e+03 1.554237e+04
[257,] 1.631627e+04 9.827235e+03
[258,] 5.138476e+04 1.631627e+04
[259,] 5.154721e+04 5.138476e+04
[260,] 5.170966e+04 5.154721e+04
[261,] 5.187211e+04 5.170966e+04
[262,] 2.537199e+04 5.187211e+04
[263,] 7.082176e+04 2.537199e+04
[264,] 5.235946e+04 7.082176e+04
[265,] 5.252191e+04 5.235946e+04
[266,] 1.614417e+04 5.252191e+04
[267,] 8.214310e+04 1.614417e+04
[268,] 5.300927e+04 8.214310e+04
[269,] 4.851515e+04 5.300927e+04
[270,] 5.333417e+04 4.851515e+04
[271,] 5.349662e+04 5.333417e+04
[272,] 5.365907e+04 5.349662e+04
[273,] 1.126404e+04 5.365907e+04
[274,] 1.658184e+04 1.126404e+04
[275,] 6.136853e+04 1.658184e+04
[276,] 1.736577e+04 6.136853e+04
[277,] 5.447132e+04 1.736577e+04
[278,] 3.439106e+04 5.447132e+04
[279,] 4.276584e+04 3.439106e+04
[280,] 5.495867e+04 4.276584e+04
[281,] 1.696481e+04 5.495867e+04
[282,] 9.013125e+04 1.696481e+04
[283,] 1.707422e+04 9.013125e+04
[284,] 1.398496e+04 1.707422e+04
[285,] 1.718364e+04 1.398496e+04
[286,] 4.434830e+04 1.718364e+04
[287,] 1.856336e+04 4.434830e+04
[288,] 1.734777e+04 1.856336e+04
[289,] 1.461269e+04 1.734777e+04
[290,] 5.658317e+04 1.461269e+04
[291,] 2.688061e+04 5.658317e+04
[292,] 6.042339e+04 2.688061e+04
[293,] 1.762131e+04 6.042339e+04
[294,] 5.220037e+04 1.762131e+04
[295,] 8.463318e+04 5.220037e+04
[296,] 1.778544e+04 8.463318e+04
[297,] 1.784015e+04 1.778544e+04
[298,] 5.788277e+04 1.784015e+04
[299,] 1.794957e+04 5.788277e+04
[300,] 5.689626e+04 1.794957e+04
[301,] 4.039842e+04 5.689626e+04
[302,] 5.853258e+04 4.039842e+04
[303,] 5.869503e+04 5.853258e+04
[304,] 5.111439e+04 5.869503e+04
[305,] 2.085978e+04 5.111439e+04
[306,] 5.918238e+04 2.085978e+04
[307,] 5.363088e+04 5.918238e+04
[308,] 1.671281e+04 5.363088e+04
[309,] 5.966973e+04 1.671281e+04
[310,] 5.983218e+04 5.966973e+04
[311,] 6.472314e+04 5.983218e+04
[312,] 6.015708e+04 6.472314e+04
[313,] 6.220608e+04 6.015708e+04
[314,] 2.008216e+04 6.220608e+04
[315,] 3.618734e+04 2.008216e+04
[316,] 1.884452e+04 3.618734e+04
[317,] -5.229339e+03 1.884452e+04
[318,] 2.071715e+04 -5.229339e+03
[319,] 1.904375e+04 2.071715e+04
[320,] 6.145668e+04 1.904375e+04
[321,] 6.161913e+04 6.145668e+04
[322,] 6.178158e+04 6.161913e+04
[323,] 1.926259e+04 6.178158e+04
[324,] 6.210648e+04 1.926259e+04
[325,] 5.970805e+04 6.210648e+04
[326,] 1.942671e+04 5.970805e+04
[327,] 1.948142e+04 1.942671e+04
[328,] 6.275628e+04 1.948142e+04
[329,] 6.291873e+04 6.275628e+04
[330,] 6.308118e+04 6.291873e+04
[331,] 1.970026e+04 6.308118e+04
[332,] 6.340608e+04 1.970026e+04
[333,] 6.219511e+04 6.340608e+04
[334,] 2.042253e+04 6.219511e+04
[335,] 1.257498e+04 2.042253e+04
[336,] 5.733947e+04 1.257498e+04
[337,] 2.002851e+04 5.733947e+04
[338,] 1.791135e+04 2.002851e+04
[339,] 7.161022e+04 1.791135e+04
[340,] -6.748344e+03 7.161022e+04
[341,] 2.024735e+04 -6.748344e+03
[342,] 2.030206e+04 2.024735e+04
[343,] 6.519304e+04 2.030206e+04
[344,] 9.977962e+03 6.519304e+04
[345,] 2.046618e+04 9.977962e+03
[346,] 9.398343e+03 2.046618e+04
[347,] 2.057560e+04 9.398343e+03
[348,] 6.600529e+04 2.057560e+04
[349,] -1.505291e+04 6.600529e+04
[350,] 7.159120e+04 -1.505291e+04
[351,] 2.017587e+04 7.159120e+04
[352,] 6.737856e+04 2.017587e+04
[353,] 1.472760e+04 6.737856e+04
[354,] 2.095856e+04 1.472760e+04
[355,] 2.101327e+04 2.095856e+04
[356,] 4.646220e+03 2.101327e+04
[357,] 2.112269e+04 4.646220e+03
[358,] 2.117740e+04 2.112269e+04
[359,] 6.779224e+04 2.117740e+04
[360,] 2.244929e+04 6.779224e+04
[361,] 6.811714e+04 2.244929e+04
[362,] 2.139624e+04 6.811714e+04
[363,] 1.748012e+04 2.139624e+04
[364,] 7.189243e+04 1.748012e+04
[365,] 6.876694e+04 7.189243e+04
[366,] 6.892939e+04 6.876694e+04
[367,] 2.166978e+04 6.892939e+04
[368,] 6.925429e+04 2.166978e+04
[369,] 2.200961e+04 6.925429e+04
[370,] 6.957919e+04 2.200961e+04
[371,] 8.591394e+04 6.957919e+04
[372,] 1.806851e+04 8.591394e+04
[373,] 4.851195e+04 1.806851e+04
[374,] 7.022900e+04 4.851195e+04
[375,] 7.039145e+04 7.022900e+04
[376,] 6.849296e+04 7.039145e+04
[377,] 1.577116e+04 6.849296e+04
[378,] -3.115255e+04 1.577116e+04
[379,] 7.104125e+04 -3.115255e+04
[380,] 7.120370e+04 7.104125e+04
[381,] 1.001894e+04 7.120370e+04
[382,] 2.087885e+04 1.001894e+04
[383,] 7.169105e+04 2.087885e+04
[384,] 2.440110e+04 7.169105e+04
[385,] -7.669326e+04 2.440110e+04
[386,] 3.391406e+04 -7.669326e+04
[387,] 4.621197e+04 3.391406e+04
[388,] 2.094429e+04 4.621197e+04
[389,] 6.024018e+04 2.094429e+04
[390,] 6.833938e+04 6.024018e+04
[391,] -8.652277e+03 6.833938e+04
[392,] 4.996653e+03 -8.652277e+03
[393,] 4.838863e+04 4.996653e+03
[394,] -1.316236e+05 4.838863e+04
[395,] 3.595084e+04 -1.316236e+05
[396,] -1.106540e+05 3.595084e+04
[397,] 8.233669e+04 -1.106540e+05
[398,] 3.655330e+04 8.233669e+04
[399,] 7.804335e+04 3.655330e+04
[400,] -2.896978e+04 7.804335e+04
[401,] 1.782120e+04 -2.896978e+04
[402,] -7.600154e+03 1.782120e+04
[403,] -5.979880e+04 -7.600154e+03
[404,] -3.410396e+04 -5.979880e+04
[405,] -7.358674e+03 -3.410396e+04
[406,] 5.045972e+04 -7.358674e+03
[407,] 4.378674e+04 5.045972e+04
[408,] 8.070822e+02 4.378674e+04
[409,] -4.814371e+04 8.070822e+02
[410,] 4.101799e+04 -4.814371e+04
[411,] -2.072469e+05 4.101799e+04
[412,] 1.086752e+04 -2.072469e+05
[413,] 9.347326e+04 1.086752e+04
[414,] -1.645637e+05 9.347326e+04
[415,] -1.301070e+05 -1.645637e+05
[416,] 1.564545e+05 -1.301070e+05
[417,] -2.943243e+04 1.564545e+05
[418,] -5.499050e+04 -2.943243e+04
[419,] -6.896898e+04 -5.499050e+04
[420,] -7.732781e+04 -6.896898e+04
[421,] -7.594499e+04 -7.732781e+04
[422,] -4.777555e+05 -7.594499e+04
[423,] -4.557375e+05 -4.777555e+05
[424,] -9.827975e+03 -4.557375e+05
[425,] -1.431767e+05 -9.827975e+03
[426,] -2.812068e+05 -1.431767e+05
[427,] -1.116491e+06 -2.812068e+05
[428,] -2.824049e+05 -1.116491e+06
[429,] -3.210303e+05 -2.824049e+05
[430,] -6.354804e+05 -3.210303e+05
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3.211947e+06 1.211341e+06
2 1.155998e+06 3.211947e+06
3 -2.198518e+06 1.155998e+06
4 -4.787482e+05 -2.198518e+06
5 3.677997e+04 -4.787482e+05
6 -6.721782e+04 3.677997e+04
7 2.157326e+05 -6.721782e+04
8 -2.498352e+05 2.157326e+05
9 7.226928e+05 -2.498352e+05
10 -3.600831e+05 7.226928e+05
11 3.744753e+05 -3.600831e+05
12 8.021790e+05 3.744753e+05
13 -1.227235e+06 8.021790e+05
14 4.304057e+04 -1.227235e+06
15 4.988105e+05 4.304057e+04
16 -1.476334e+05 4.988105e+05
17 2.531232e+04 -1.476334e+05
18 -2.893929e+05 2.531232e+04
19 -3.512555e+05 -2.893929e+05
20 1.698836e+05 -3.512555e+05
21 -2.277514e+05 1.698836e+05
22 -9.484326e+04 -2.277514e+05
23 -4.115152e+05 -9.484326e+04
24 -3.273638e+05 -4.115152e+05
25 -1.754206e+05 -3.273638e+05
26 -3.809225e+05 -1.754206e+05
27 -3.030669e+05 -3.809225e+05
28 3.024907e+05 -3.030669e+05
29 -9.993122e+04 3.024907e+05
30 2.419794e+05 -9.993122e+04
31 -1.431770e+05 2.419794e+05
32 2.200782e+05 -1.431770e+05
33 9.470666e+04 2.200782e+05
34 -2.285539e+05 9.470666e+04
35 -5.024220e+05 -2.285539e+05
36 9.737859e+03 -5.024220e+05
37 1.637956e+05 9.737859e+03
38 -6.253546e+04 1.637956e+05
39 -1.854367e+05 -6.253546e+04
40 -9.180229e+04 -1.854367e+05
41 -4.195152e+05 -9.180229e+04
42 8.528477e+04 -4.195152e+05
43 -1.695808e+05 8.528477e+04
44 -4.680120e+04 -1.695808e+05
45 -1.096396e+05 -4.680120e+04
46 -2.898785e+05 -1.096396e+05
47 -9.963793e+04 -2.898785e+05
48 -6.815482e+04 -9.963793e+04
49 -1.905060e+05 -6.815482e+04
50 8.854448e+04 -1.905060e+05
51 9.955848e+03 8.854448e+04
52 -2.460721e+05 9.955848e+03
53 1.169643e+05 -2.460721e+05
54 -3.324518e+05 1.169643e+05
55 5.600280e+04 -3.324518e+05
56 2.680244e+05 5.600280e+04
57 -3.197073e+04 2.680244e+05
58 -2.417528e+04 -3.197073e+04
59 -2.010330e+05 -2.417528e+04
60 -4.553505e+03 -2.010330e+05
61 -1.951059e+05 -4.553505e+03
62 -3.265241e+04 -1.951059e+05
63 -1.616639e+05 -3.265241e+04
64 -4.611929e+03 -1.616639e+05
65 6.337081e+04 -4.611929e+03
66 -6.702604e+04 6.337081e+04
67 -1.928274e+04 -6.702604e+04
68 -6.731344e+04 -1.928274e+04
69 -4.661021e+04 -6.731344e+04
70 -8.543808e+04 -4.661021e+04
71 -3.135755e+04 -8.543808e+04
72 -9.129420e+04 -3.135755e+04
73 -1.950571e+04 -9.129420e+04
74 -6.682351e+04 -1.950571e+04
75 -1.269023e+05 -6.682351e+04
76 -4.868155e+04 -1.269023e+05
77 -8.848015e+04 -4.868155e+04
78 1.945406e+04 -8.848015e+04
79 -1.078963e+04 1.945406e+04
80 2.459810e+04 -1.078963e+04
81 -3.172784e+04 2.459810e+04
82 -1.045663e+05 -3.172784e+04
83 3.465232e+03 -1.045663e+05
84 -1.285915e+04 3.465232e+03
85 -8.075890e+04 -1.285915e+04
86 3.450340e+04 -8.075890e+04
87 -2.397391e+04 3.450340e+04
88 -1.020981e+05 -2.397391e+04
89 7.308063e+03 -1.020981e+05
90 -3.186995e+03 7.308063e+03
91 -2.994864e+04 -3.186995e+03
92 1.356350e+05 -2.994864e+04
93 -2.302266e+05 1.356350e+05
94 8.217000e+04 -2.302266e+05
95 1.503895e+04 8.217000e+04
96 2.264769e+04 1.503895e+04
97 2.786399e+04 2.264769e+04
98 -6.364448e+04 2.786399e+04
99 -3.553867e+04 -6.364448e+04
100 -3.724832e+04 -3.553867e+04
101 -1.365407e+05 -3.724832e+04
102 -1.554944e+05 -1.365407e+05
103 -9.722487e+04 -1.554944e+05
104 3.039562e+04 -9.722487e+04
105 -5.377295e+04 3.039562e+04
106 -2.292727e+05 -5.377295e+04
107 9.015128e+04 -2.292727e+05
108 -1.827500e+05 9.015128e+04
109 -1.057083e+05 -1.827500e+05
110 9.167128e+04 -1.057083e+05
111 -1.840696e+05 9.167128e+04
112 -1.726568e+05 -1.840696e+05
113 -9.290390e+04 -1.726568e+05
114 5.943598e+04 -9.290390e+04
115 1.114294e+04 5.943598e+04
116 -1.096814e+04 1.114294e+04
117 4.140871e+04 -1.096814e+04
118 -1.682057e+05 4.140871e+04
119 3.665309e+04 -1.682057e+05
120 -2.953918e+05 3.665309e+04
121 -1.754003e+03 -2.953918e+05
122 5.227681e+03 -1.754003e+03
123 5.787494e+04 5.227681e+03
124 -4.364669e+01 5.787494e+04
125 -1.901348e+05 -4.364669e+01
126 8.484910e+03 -1.901348e+05
127 7.015435e+03 8.484910e+03
128 1.140951e+04 7.015435e+03
129 3.042868e+04 1.140951e+04
130 9.194867e+03 3.042868e+04
131 7.933457e+04 9.194867e+03
132 3.091603e+04 7.933457e+04
133 5.785050e+03 3.091603e+04
134 1.067631e+04 5.785050e+03
135 2.679422e+03 1.067631e+04
136 3.755391e+04 2.679422e+03
137 9.041119e+03 3.755391e+04
138 -1.500719e+04 9.041119e+03
139 9.196128e+03 -1.500719e+04
140 1.490940e+04 9.196128e+03
141 3.237808e+04 1.490940e+04
142 -2.580087e+04 3.237808e+04
143 9.506997e+03 -2.580087e+04
144 3.503192e+04 9.506997e+03
145 7.217910e+03 3.503192e+04
146 3.319033e+04 7.217910e+03
147 1.059662e+04 3.319033e+04
148 3.608238e+04 1.059662e+04
149 3.361046e+04 3.608238e+04
150 -2.279056e+04 3.361046e+04
151 -2.725782e+04 -2.279056e+04
152 3.416503e+04 -2.725782e+04
153 1.264349e+04 3.416503e+04
154 9.909246e+03 1.264349e+04
155 1.007147e+04 9.909246e+03
156 1.012618e+04 1.007147e+04
157 6.841279e+03 1.012618e+04
158 2.181349e+04 6.841279e+03
159 1.029031e+04 2.181349e+04
160 1.034502e+04 1.029031e+04
161 3.562709e+04 1.034502e+04
162 1.045443e+04 3.562709e+04
163 2.524909e+04 1.045443e+04
164 2.895824e+04 2.524909e+04
165 9.684151e+03 2.895824e+04
166 3.643934e+04 9.684151e+03
167 1.223577e+05 3.643934e+04
168 1.159441e+05 1.223577e+05
169 1.083740e+04 1.159441e+05
170 1.089211e+04 1.083740e+04
171 3.336618e+04 1.089211e+04
172 3.741404e+04 3.336618e+04
173 2.744973e+04 3.741404e+04
174 3.773894e+04 2.744973e+04
175 1.042717e+04 3.773894e+04
176 9.931973e+04 1.042717e+04
177 3.729259e+04 9.931973e+04
178 1.132978e+04 3.729259e+04
179 1.138449e+04 1.132978e+04
180 4.646463e+04 1.138449e+04
181 -4.542231e+02 4.646463e+04
182 -3.517192e+04 -4.542231e+02
183 -1.668223e+04 -3.517192e+04
184 3.664797e+04 -1.668223e+04
185 4.553636e+04 3.664797e+04
186 -1.305767e+04 4.553636e+04
187 5.289045e+03 -1.305767e+04
188 1.241147e+03 5.289045e+03
189 1.193158e+04 1.241147e+03
190 3.988307e+04 1.193158e+04
191 1.204100e+04 3.988307e+04
192 8.667379e+03 1.204100e+04
193 4.082549e+04 8.667379e+03
194 1.208704e+04 4.082549e+04
195 4.115040e+04 1.208704e+04
196 -8.454995e+04 4.115040e+04
197 4.147530e+04 -8.454995e+04
198 4.163775e+04 4.147530e+04
199 -4.982577e+04 4.163775e+04
200 4.196265e+04 -4.982577e+04
201 4.212510e+04 4.196265e+04
202 4.228755e+04 4.212510e+04
203 5.645678e+04 4.228755e+04
204 -1.105799e+04 5.645678e+04
205 4.223679e+04 -1.105799e+04
206 3.556701e+04 4.223679e+04
207 4.309980e+04 3.556701e+04
208 3.515670e+04 4.309980e+04
209 6.271365e+03 3.515670e+04
210 1.308047e+04 6.271365e+03
211 4.374960e+04 1.308047e+04
212 7.292312e+04 4.374960e+04
213 1.042869e+04 7.292312e+04
214 1.641191e+04 1.042869e+04
215 4.439940e+04 1.641191e+04
216 4.815296e+04 4.439940e+04
217 1.301405e+04 4.815296e+04
218 1.351814e+04 1.301405e+04
219 -1.723142e+04 1.351814e+04
220 2.010942e+04 -1.723142e+04
221 4.537410e+04 2.010942e+04
222 1.424405e+04 4.537410e+04
223 3.859880e+04 1.424405e+04
224 1.395492e+04 3.859880e+04
225 1.236742e+04 1.395492e+04
226 1.552967e+04 1.236742e+04
227 4.634880e+04 1.552967e+04
228 2.983770e+04 4.634880e+04
229 1.411994e+04 2.983770e+04
230 1.809488e+04 1.411994e+04
231 3.537509e+04 1.809488e+04
232 4.716106e+04 3.537509e+04
233 -4.433804e+03 4.716106e+04
234 1.439348e+04 -4.433804e+03
235 1.476199e+04 1.439348e+04
236 1.450290e+04 1.476199e+04
237 1.574626e+04 1.450290e+04
238 4.813576e+04 1.574626e+04
239 4.829821e+04 4.813576e+04
240 1.472174e+04 4.829821e+04
241 4.337584e+04 1.472174e+04
242 4.878556e+04 4.337584e+04
243 5.437915e+04 4.878556e+04
244 -3.686962e+04 5.437915e+04
245 4.927291e+04 -3.686962e+04
246 4.943536e+04 4.927291e+04
247 6.809405e+03 4.943536e+04
248 4.976026e+04 6.809405e+03
249 1.521412e+04 4.976026e+04
250 4.585367e+04 1.521412e+04
251 5.024761e+04 4.585367e+04
252 3.069504e+04 5.024761e+04
253 -8.086104e+03 3.069504e+04
254 -1.135608e+05 -8.086104e+03
255 1.554237e+04 -1.135608e+05
256 9.827235e+03 1.554237e+04
257 1.631627e+04 9.827235e+03
258 5.138476e+04 1.631627e+04
259 5.154721e+04 5.138476e+04
260 5.170966e+04 5.154721e+04
261 5.187211e+04 5.170966e+04
262 2.537199e+04 5.187211e+04
263 7.082176e+04 2.537199e+04
264 5.235946e+04 7.082176e+04
265 5.252191e+04 5.235946e+04
266 1.614417e+04 5.252191e+04
267 8.214310e+04 1.614417e+04
268 5.300927e+04 8.214310e+04
269 4.851515e+04 5.300927e+04
270 5.333417e+04 4.851515e+04
271 5.349662e+04 5.333417e+04
272 5.365907e+04 5.349662e+04
273 1.126404e+04 5.365907e+04
274 1.658184e+04 1.126404e+04
275 6.136853e+04 1.658184e+04
276 1.736577e+04 6.136853e+04
277 5.447132e+04 1.736577e+04
278 3.439106e+04 5.447132e+04
279 4.276584e+04 3.439106e+04
280 5.495867e+04 4.276584e+04
281 1.696481e+04 5.495867e+04
282 9.013125e+04 1.696481e+04
283 1.707422e+04 9.013125e+04
284 1.398496e+04 1.707422e+04
285 1.718364e+04 1.398496e+04
286 4.434830e+04 1.718364e+04
287 1.856336e+04 4.434830e+04
288 1.734777e+04 1.856336e+04
289 1.461269e+04 1.734777e+04
290 5.658317e+04 1.461269e+04
291 2.688061e+04 5.658317e+04
292 6.042339e+04 2.688061e+04
293 1.762131e+04 6.042339e+04
294 5.220037e+04 1.762131e+04
295 8.463318e+04 5.220037e+04
296 1.778544e+04 8.463318e+04
297 1.784015e+04 1.778544e+04
298 5.788277e+04 1.784015e+04
299 1.794957e+04 5.788277e+04
300 5.689626e+04 1.794957e+04
301 4.039842e+04 5.689626e+04
302 5.853258e+04 4.039842e+04
303 5.869503e+04 5.853258e+04
304 5.111439e+04 5.869503e+04
305 2.085978e+04 5.111439e+04
306 5.918238e+04 2.085978e+04
307 5.363088e+04 5.918238e+04
308 1.671281e+04 5.363088e+04
309 5.966973e+04 1.671281e+04
310 5.983218e+04 5.966973e+04
311 6.472314e+04 5.983218e+04
312 6.015708e+04 6.472314e+04
313 6.220608e+04 6.015708e+04
314 2.008216e+04 6.220608e+04
315 3.618734e+04 2.008216e+04
316 1.884452e+04 3.618734e+04
317 -5.229339e+03 1.884452e+04
318 2.071715e+04 -5.229339e+03
319 1.904375e+04 2.071715e+04
320 6.145668e+04 1.904375e+04
321 6.161913e+04 6.145668e+04
322 6.178158e+04 6.161913e+04
323 1.926259e+04 6.178158e+04
324 6.210648e+04 1.926259e+04
325 5.970805e+04 6.210648e+04
326 1.942671e+04 5.970805e+04
327 1.948142e+04 1.942671e+04
328 6.275628e+04 1.948142e+04
329 6.291873e+04 6.275628e+04
330 6.308118e+04 6.291873e+04
331 1.970026e+04 6.308118e+04
332 6.340608e+04 1.970026e+04
333 6.219511e+04 6.340608e+04
334 2.042253e+04 6.219511e+04
335 1.257498e+04 2.042253e+04
336 5.733947e+04 1.257498e+04
337 2.002851e+04 5.733947e+04
338 1.791135e+04 2.002851e+04
339 7.161022e+04 1.791135e+04
340 -6.748344e+03 7.161022e+04
341 2.024735e+04 -6.748344e+03
342 2.030206e+04 2.024735e+04
343 6.519304e+04 2.030206e+04
344 9.977962e+03 6.519304e+04
345 2.046618e+04 9.977962e+03
346 9.398343e+03 2.046618e+04
347 2.057560e+04 9.398343e+03
348 6.600529e+04 2.057560e+04
349 -1.505291e+04 6.600529e+04
350 7.159120e+04 -1.505291e+04
351 2.017587e+04 7.159120e+04
352 6.737856e+04 2.017587e+04
353 1.472760e+04 6.737856e+04
354 2.095856e+04 1.472760e+04
355 2.101327e+04 2.095856e+04
356 4.646220e+03 2.101327e+04
357 2.112269e+04 4.646220e+03
358 2.117740e+04 2.112269e+04
359 6.779224e+04 2.117740e+04
360 2.244929e+04 6.779224e+04
361 6.811714e+04 2.244929e+04
362 2.139624e+04 6.811714e+04
363 1.748012e+04 2.139624e+04
364 7.189243e+04 1.748012e+04
365 6.876694e+04 7.189243e+04
366 6.892939e+04 6.876694e+04
367 2.166978e+04 6.892939e+04
368 6.925429e+04 2.166978e+04
369 2.200961e+04 6.925429e+04
370 6.957919e+04 2.200961e+04
371 8.591394e+04 6.957919e+04
372 1.806851e+04 8.591394e+04
373 4.851195e+04 1.806851e+04
374 7.022900e+04 4.851195e+04
375 7.039145e+04 7.022900e+04
376 6.849296e+04 7.039145e+04
377 1.577116e+04 6.849296e+04
378 -3.115255e+04 1.577116e+04
379 7.104125e+04 -3.115255e+04
380 7.120370e+04 7.104125e+04
381 1.001894e+04 7.120370e+04
382 2.087885e+04 1.001894e+04
383 7.169105e+04 2.087885e+04
384 2.440110e+04 7.169105e+04
385 -7.669326e+04 2.440110e+04
386 3.391406e+04 -7.669326e+04
387 4.621197e+04 3.391406e+04
388 2.094429e+04 4.621197e+04
389 6.024018e+04 2.094429e+04
390 6.833938e+04 6.024018e+04
391 -8.652277e+03 6.833938e+04
392 4.996653e+03 -8.652277e+03
393 4.838863e+04 4.996653e+03
394 -1.316236e+05 4.838863e+04
395 3.595084e+04 -1.316236e+05
396 -1.106540e+05 3.595084e+04
397 8.233669e+04 -1.106540e+05
398 3.655330e+04 8.233669e+04
399 7.804335e+04 3.655330e+04
400 -2.896978e+04 7.804335e+04
401 1.782120e+04 -2.896978e+04
402 -7.600154e+03 1.782120e+04
403 -5.979880e+04 -7.600154e+03
404 -3.410396e+04 -5.979880e+04
405 -7.358674e+03 -3.410396e+04
406 5.045972e+04 -7.358674e+03
407 4.378674e+04 5.045972e+04
408 8.070822e+02 4.378674e+04
409 -4.814371e+04 8.070822e+02
410 4.101799e+04 -4.814371e+04
411 -2.072469e+05 4.101799e+04
412 1.086752e+04 -2.072469e+05
413 9.347326e+04 1.086752e+04
414 -1.645637e+05 9.347326e+04
415 -1.301070e+05 -1.645637e+05
416 1.564545e+05 -1.301070e+05
417 -2.943243e+04 1.564545e+05
418 -5.499050e+04 -2.943243e+04
419 -6.896898e+04 -5.499050e+04
420 -7.732781e+04 -6.896898e+04
421 -7.594499e+04 -7.732781e+04
422 -4.777555e+05 -7.594499e+04
423 -4.557375e+05 -4.777555e+05
424 -9.827975e+03 -4.557375e+05
425 -1.431767e+05 -9.827975e+03
426 -2.812068e+05 -1.431767e+05
427 -1.116491e+06 -2.812068e+05
428 -2.824049e+05 -1.116491e+06
429 -3.210303e+05 -2.824049e+05
430 -6.354804e+05 -3.210303e+05
> 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/7kjkf1291239858.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/8kjkf1291239858.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/9kjkf1291239858.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/10vs101291239858.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/11gti51291239858.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/12jbgb1291239858.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/138ud51291239858.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/14jmv81291239858.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/154mbe1291239858.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/160w951291239858.tab")
+ }
>
> try(system("convert tmp/16rmo1291239858.ps tmp/16rmo1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/26rmo1291239858.ps tmp/26rmo1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/3yjmr1291239858.ps tmp/3yjmr1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/4yjmr1291239858.ps tmp/4yjmr1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/5yjmr1291239858.ps tmp/5yjmr1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/6rs3u1291239858.ps tmp/6rs3u1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/7kjkf1291239858.ps tmp/7kjkf1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/8kjkf1291239858.ps tmp/8kjkf1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/9kjkf1291239858.ps tmp/9kjkf1291239858.png",intern=TRUE))
character(0)
> try(system("convert tmp/10vs101291239858.ps tmp/10vs101291239858.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
15.640 2.377 32.886