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(162556
+ ,1081
+ ,213118
+ ,6282154
+ ,1
+ ,5626
+ ,37
+ ,56289
+ ,29790
+ ,309
+ ,81767
+ ,4321023
+ ,2
+ ,13337
+ ,138
+ ,28328
+ ,87550
+ ,458
+ ,153198
+ ,4111912
+ ,3
+ ,8541
+ ,45
+ ,17936
+ ,84738
+ ,588
+ ,-26007
+ ,223193
+ ,415
+ ,39
+ ,0
+ ,4145
+ ,54660
+ ,302
+ ,126942
+ ,1491348
+ ,6
+ ,4276
+ ,24
+ ,3040
+ ,42634
+ ,156
+ ,157214
+ ,1629616
+ ,5
+ ,9164
+ ,34
+ ,2964
+ ,40949
+ ,481
+ ,129352
+ ,1398893
+ ,9
+ ,2493
+ ,29
+ ,2865
+ ,45187
+ ,353
+ ,234817
+ ,1926517
+ ,4
+ ,4891
+ ,38
+ ,2854
+ ,37704
+ ,452
+ ,60448
+ ,983660
+ ,14
+ ,1734
+ ,21
+ ,2167
+ ,16275
+ ,109
+ ,47818
+ ,1443586
+ ,7
+ ,11409
+ ,76
+ ,1974
+ ,25830
+ ,115
+ ,245546
+ ,1073089
+ ,10
+ ,7592
+ ,34
+ ,1910
+ ,12679
+ ,110
+ ,48020
+ ,984885
+ ,13
+ ,7135
+ ,62
+ ,1871
+ ,18014
+ ,239
+ ,-1710
+ ,1405225
+ ,8
+ ,5043
+ ,67
+ ,943
+ ,43556
+ ,247
+ ,32648
+ ,227132
+ ,414
+ ,110
+ ,1
+ ,929
+ ,24811
+ ,505
+ ,95350
+ ,929118
+ ,15
+ ,1444
+ ,29
+ ,822
+ ,6575
+ ,159
+ ,151352
+ ,1071292
+ ,11
+ ,5480
+ ,133
+ ,819
+ ,7123
+ ,109
+ ,288170
+ ,638830
+ ,28
+ ,4026
+ ,62
+ ,769
+ ,21950
+ ,519
+ ,114337
+ ,856956
+ ,17
+ ,1266
+ ,30
+ ,745
+ ,37597
+ ,248
+ ,37884
+ ,992426
+ ,12
+ ,3195
+ ,21
+ ,652
+ ,17821
+ ,373
+ ,122844
+ ,444477
+ ,46
+ ,655
+ ,14
+ ,643
+ ,12988
+ ,119
+ ,82340
+ ,857217
+ ,16
+ ,5523
+ ,51
+ ,601
+ ,22330
+ ,84
+ ,79801
+ ,711969
+ ,20
+ ,6095
+ ,23
+ ,446
+ ,13326
+ ,102
+ ,165548
+ ,702380
+ ,21
+ ,4925
+ ,38
+ ,436
+ ,16189
+ ,295
+ ,116384
+ ,358589
+ ,89
+ ,538
+ ,10
+ ,379
+ ,7146
+ ,105
+ ,134028
+ ,297978
+ ,367
+ ,933
+ ,14
+ ,305
+ ,15824
+ ,64
+ ,63838
+ ,585715
+ ,30
+ ,6027
+ ,24
+ ,284
+ ,27664
+ ,282
+ ,74996
+ ,657954
+ ,24
+ ,1624
+ ,17
+ ,247
+ ,11920
+ ,182
+ ,31080
+ ,209458
+ ,419
+ ,52
+ ,1
+ ,238
+ ,8568
+ ,37
+ ,32168
+ ,786690
+ ,18
+ ,15856
+ ,68
+ ,223
+ ,14416
+ ,361
+ ,49857
+ ,439798
+ ,48
+ ,664
+ ,17
+ ,220
+ ,3369
+ ,28
+ ,87161
+ ,688779
+ ,23
+ ,17456
+ ,145
+ ,217
+ ,11819
+ ,85
+ ,106113
+ ,574339
+ ,31
+ ,4404
+ ,32
+ ,199
+ ,6984
+ ,45
+ ,80570
+ ,741409
+ ,19
+ ,12031
+ ,78
+ ,195
+ ,4519
+ ,49
+ ,102129
+ ,597793
+ ,29
+ ,8118
+ ,88
+ ,163
+ ,2220
+ ,22
+ ,301670
+ ,644190
+ ,26
+ ,20190
+ ,200
+ ,154
+ ,18562
+ ,155
+ ,102313
+ ,377934
+ ,72
+ ,1148
+ ,10
+ ,143
+ ,10327
+ ,91
+ ,88577
+ ,640273
+ ,27
+ ,4838
+ ,43
+ ,135
+ ,5336
+ ,81
+ ,112477
+ ,697458
+ ,22
+ ,6141
+ ,93
+ ,132
+ ,2365
+ ,79
+ ,191778
+ ,550608
+ ,32
+ ,4438
+ ,148
+ ,120
+ ,4069
+ ,145
+ ,79804
+ ,207393
+ ,420
+ ,51
+ ,2
+ ,119
+ ,8636
+ ,855
+ ,128294
+ ,301607
+ ,359
+ ,119
+ ,12
+ ,108
+ ,13718
+ ,61
+ ,96448
+ ,345783
+ ,100
+ ,2390
+ ,11
+ ,104
+ ,4525
+ ,226
+ ,93811
+ ,501749
+ ,37
+ ,1335
+ ,67
+ ,101
+ ,6869
+ ,105
+ ,117520
+ ,379983
+ ,68
+ ,1714
+ ,26
+ ,90
+ ,4628
+ ,62
+ ,69159
+ ,387475
+ ,64
+ ,3024
+ ,41
+ ,85
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,74
+ ,7092
+ ,48
+ ,74
+ ,4891
+ ,217
+ ,210568
+ ,370837
+ ,79
+ ,787
+ ,35
+ ,73
+ ,7489
+ ,322
+ ,136996
+ ,430866
+ ,50
+ ,717
+ ,31
+ ,70
+ ,4901
+ ,84
+ ,121920
+ ,469107
+ ,41
+ ,3204
+ ,55
+ ,67
+ ,2284
+ ,33
+ ,76403
+ ,194493
+ ,423
+ ,-167
+ ,-2
+ ,66
+ ,3160
+ ,108
+ ,108094
+ ,530670
+ ,33
+ ,3062
+ ,105
+ ,66
+ ,4150
+ ,150
+ ,134759
+ ,518365
+ ,35
+ ,2122
+ ,77
+ ,66
+ ,7285
+ ,115
+ ,188873
+ ,491303
+ ,39
+ ,2533
+ ,40
+ ,59
+ ,1134
+ ,162
+ ,146216
+ ,527021
+ ,34
+ ,2019
+ ,288
+ ,58
+ ,4658
+ ,158
+ ,156608
+ ,233773
+ ,411
+ ,214
+ ,7
+ ,58
+ ,2384
+ ,97
+ ,61348
+ ,405972
+ ,58
+ ,2123
+ ,86
+ ,54
+ ,3748
+ ,9
+ ,50350
+ ,652925
+ ,25
+ ,50325
+ ,121
+ ,53
+ ,5371
+ ,66
+ ,87720
+ ,446211
+ ,45
+ ,3730
+ ,46
+ ,49
+ ,1285
+ ,107
+ ,99489
+ ,341340
+ ,106
+ ,1321
+ ,110
+ ,49
+ ,9327
+ ,101
+ ,87419
+ ,387699
+ ,63
+ ,1858
+ ,20
+ ,44
+ ,5565
+ ,47
+ ,94355
+ ,493408
+ ,38
+ ,6243
+ ,53
+ ,39
+ ,1528
+ ,38
+ ,60326
+ ,146494
+ ,426
+ ,-1408
+ ,-35
+ ,38
+ ,3122
+ ,34
+ ,94670
+ ,414462
+ ,54
+ ,6308
+ ,69
+ ,38
+ ,7561
+ ,87
+ ,82425
+ ,364304
+ ,86
+ ,1889
+ ,22
+ ,37
+ ,2675
+ ,79
+ ,59017
+ ,355178
+ ,93
+ ,1964
+ ,58
+ ,36
+ ,13253
+ ,947
+ ,90829
+ ,357760
+ ,90
+ ,167
+ ,12
+ ,35
+ ,880
+ ,74
+ ,80791
+ ,261216
+ ,399
+ ,827
+ ,70
+ ,34
+ ,2053
+ ,53
+ ,100423
+ ,397144
+ ,60
+ ,3720
+ ,96
+ ,33
+ ,1424
+ ,94
+ ,131116
+ ,374943
+ ,77
+ ,1861
+ ,123
+ ,32
+ ,4036
+ ,63
+ ,100269
+ ,424898
+ ,53
+ ,3570
+ ,56
+ ,32
+ ,3045
+ ,58
+ ,27330
+ ,202055
+ ,422
+ ,35
+ ,1
+ ,31
+ ,5119
+ ,49
+ ,39039
+ ,378525
+ ,69
+ ,3643
+ ,35
+ ,31
+ ,1431
+ ,34
+ ,106885
+ ,310768
+ ,330
+ ,3258
+ ,77
+ ,30
+ ,554
+ ,11
+ ,79285
+ ,325738
+ ,125
+ ,11431
+ ,227
+ ,30
+ ,1975
+ ,35
+ ,118881
+ ,394510
+ ,61
+ ,5557
+ ,98
+ ,30
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,404
+ ,2353
+ ,27
+ ,28
+ ,1012
+ ,47
+ ,114768
+ ,368078
+ ,82
+ ,3576
+ ,166
+ ,26
+ ,810
+ ,43
+ ,74015
+ ,236761
+ ,409
+ ,855
+ ,45
+ ,26
+ ,1280
+ ,117
+ ,69465
+ ,312378
+ ,326
+ ,960
+ ,88
+ ,26
+ ,666
+ ,171
+ ,117869
+ ,339836
+ ,108
+ ,818
+ ,210
+ ,25
+ ,1380
+ ,26
+ ,60982
+ ,347385
+ ,98
+ ,5669
+ ,107
+ ,25
+ ,4677
+ ,75
+ ,90131
+ ,426280
+ ,52
+ ,3017
+ ,48
+ ,24
+ ,876
+ ,59
+ ,138971
+ ,352850
+ ,96
+ ,2591
+ ,174
+ ,24
+ ,814
+ ,18
+ ,39625
+ ,301881
+ ,358
+ ,5660
+ ,125
+ ,23
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,73
+ ,11834
+ ,346
+ ,23
+ ,5692
+ ,72
+ ,64239
+ ,357312
+ ,91
+ ,2185
+ ,28
+ ,22
+ ,3642
+ ,86
+ ,90262
+ ,458343
+ ,42
+ ,3004
+ ,71
+ ,22
+ ,540
+ ,14
+ ,103960
+ ,354228
+ ,94
+ ,11016
+ ,286
+ ,22
+ ,2099
+ ,64
+ ,106611
+ ,308636
+ ,339
+ ,1697
+ ,52
+ ,20
+ ,567
+ ,11
+ ,103345
+ ,386212
+ ,65
+ ,16928
+ ,329
+ ,20
+ ,2001
+ ,52
+ ,95551
+ ,393343
+ ,62
+ ,3718
+ ,97
+ ,19
+ ,2949
+ ,41
+ ,82903
+ ,378509
+ ,70
+ ,4354
+ ,61
+ ,18
+ ,2253
+ ,99
+ ,63593
+ ,452469
+ ,43
+ ,2550
+ ,112
+ ,17
+ ,6533
+ ,75
+ ,126910
+ ,364839
+ ,85
+ ,2198
+ ,25
+ ,16
+ ,1889
+ ,45
+ ,37527
+ ,358649
+ ,88
+ ,3526
+ ,84
+ ,14
+ ,3055
+ ,43
+ ,60247
+ ,376641
+ ,75
+ ,4108
+ ,58
+ ,14
+ ,272
+ ,8
+ ,112995
+ ,429112
+ ,51
+ ,28639
+ ,844
+ ,14
+ ,1414
+ ,198
+ ,70184
+ ,330546
+ ,119
+ ,659
+ ,92
+ ,13
+ ,2564
+ ,22
+ ,130140
+ ,403560
+ ,59
+ ,9253
+ ,79
+ ,13
+ ,1383
+ ,11
+ ,73221
+ ,317892
+ ,151
+ ,10717
+ ,85
+ ,13
+ ,1261
+ ,33
+ ,76114
+ ,307528
+ ,343
+ ,3258
+ ,85
+ ,13
+ ,975
+ ,23
+ ,90534
+ ,235133
+ ,410
+ ,1528
+ ,36
+ ,13
+ ,3366
+ ,80
+ ,108479
+ ,299243
+ ,362
+ ,1241
+ ,29
+ ,12
+ ,576
+ ,18
+ ,113761
+ ,314073
+ ,313
+ ,6337
+ ,198
+ ,12
+ ,1686
+ ,40
+ ,68696
+ ,368186
+ ,81
+ ,4205
+ ,100
+ ,11
+ ,746
+ ,23
+ ,71561
+ ,269661
+ ,396
+ ,3029
+ ,93
+ ,10
+ ,3192
+ ,60
+ ,59831
+ ,125390
+ ,427
+ ,-1243
+ ,-23
+ ,10
+ ,2045
+ ,20
+ ,97890
+ ,510834
+ ,36
+ ,15542
+ ,152
+ ,10
+ ,5702
+ ,61
+ ,101481
+ ,321896
+ ,138
+ ,1998
+ ,21
+ ,10
+ ,1932
+ ,36
+ ,72954
+ ,249898
+ ,403
+ ,1386
+ ,26
+ ,10
+ ,936
+ ,30
+ ,67939
+ ,408881
+ ,56
+ ,6963
+ ,223
+ ,9
+ ,3437
+ ,47
+ ,48022
+ ,158492
+ ,425
+ ,-883
+ ,-12
+ ,9
+ ,5131
+ ,71
+ ,86111
+ ,292154
+ ,379
+ ,1298
+ ,18
+ ,9
+ ,2397
+ ,14
+ ,74020
+ ,289513
+ ,385
+ ,6394
+ ,37
+ ,6
+ ,1389
+ ,9
+ ,57530
+ ,378049
+ ,71
+ ,19783
+ ,128
+ ,6
+ ,1503
+ ,39
+ ,56364
+ ,343466
+ ,104
+ ,3679
+ ,95
+ ,6
+ ,402
+ ,26
+ ,84990
+ ,332743
+ ,115
+ ,5106
+ ,330
+ ,5
+ ,2239
+ ,21
+ ,88590
+ ,442882
+ ,47
+ ,11566
+ ,108
+ ,5
+ ,2234
+ ,16
+ ,77200
+ ,214215
+ ,418
+ ,888
+ ,6
+ ,5
+ ,837
+ ,69
+ ,61262
+ ,315688
+ ,165
+ ,1677
+ ,138
+ ,4
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,76
+ ,1904
+ ,17
+ ,4
+ ,875
+ ,14
+ ,67000
+ ,334280
+ ,112
+ ,9591
+ ,153
+ ,3
+ ,1585
+ ,107
+ ,93099
+ ,355864
+ ,92
+ ,1457
+ ,98
+ ,2
+ ,1659
+ ,29
+ ,107577
+ ,480382
+ ,40
+ ,9668
+ ,169
+ ,2
+ ,2647
+ ,37
+ ,62920
+ ,353058
+ ,95
+ ,4137
+ ,58
+ ,1
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,416
+ ,748
+ ,5
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,191
+ ,0
+ ,94
+ ,7
+ ,60793
+ ,314533
+ ,308
+ ,16362
+ ,1212
+ ,0
+ ,422
+ ,28
+ ,57935
+ ,318056
+ ,150
+ ,4216
+ ,280
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,197
+ ,0
+ ,34
+ ,8
+ ,60630
+ ,314353
+ ,309
+ ,14294
+ ,3340
+ ,0
+ ,1558
+ ,63
+ ,55637
+ ,369448
+ ,80
+ ,2690
+ ,109
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,196
+ ,0
+ ,43
+ ,3
+ ,60887
+ ,312846
+ ,323
+ ,37615
+ ,2613
+ ,0
+ ,645
+ ,5
+ ,60720
+ ,312075
+ ,329
+ ,22415
+ ,174
+ ,0
+ ,316
+ ,9
+ ,60505
+ ,315009
+ ,299
+ ,12779
+ ,364
+ ,0
+ ,115
+ ,13
+ ,60945
+ ,318903
+ ,145
+ ,9146
+ ,1037
+ ,0
+ ,5
+ ,2
+ ,60720
+ ,314887
+ ,304
+ ,57443
+ ,24814
+ ,0
+ ,897
+ ,5
+ ,60720
+ ,314913
+ ,303
+ ,22983
+ ,128
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,188
+ ,0
+ ,389
+ ,14
+ ,58990
+ ,325506
+ ,127
+ ,8965
+ ,322
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,220
+ ,0
+ ,1002
+ ,15
+ ,56750
+ ,298568
+ ,364
+ ,6571
+ ,98
+ ,0
+ ,36
+ ,3
+ ,60894
+ ,315834
+ ,162
+ ,38611
+ ,3215
+ ,0
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,122
+ ,8652
+ ,282
+ ,0
+ ,309
+ ,11
+ ,56535
+ ,312878
+ ,322
+ ,10262
+ ,365
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,0
+ ,9
+ ,6
+ ,60835
+ ,314987
+ ,300
+ ,19165
+ ,12458
+ ,0
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,129
+ ,62625
+ ,463
+ ,0
+ ,14
+ ,1
+ ,61016
+ ,315877
+ ,161
+ ,115877
+ ,8041
+ ,0
+ ,520
+ ,10
+ ,58650
+ ,291650
+ ,383
+ ,9165
+ ,176
+ ,0
+ ,1766
+ ,73
+ ,60438
+ ,305959
+ ,349
+ ,1451
+ ,60
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,208
+ ,0
+ ,458
+ ,11
+ ,58625
+ ,297765
+ ,368
+ ,8888
+ ,213
+ ,0
+ ,20
+ ,3
+ ,60938
+ ,315245
+ ,297
+ ,38415
+ ,5674
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,217
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,219
+ ,0
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,298
+ ,57618
+ ,1179
+ ,0
+ ,405
+ ,7
+ ,60845
+ ,336425
+ ,111
+ ,19489
+ ,337
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,213
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,212
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,207
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,218
+ ,0
+ ,483
+ ,27
+ ,60830
+ ,306268
+ ,348
+ ,3936
+ ,220
+ ,0
+ ,454
+ ,51
+ ,63261
+ ,302187
+ ,357
+ ,2004
+ ,225
+ ,0
+ ,47
+ ,3
+ ,60720
+ ,314882
+ ,305
+ ,38294
+ ,2432
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,193
+ ,0
+ ,757
+ ,19
+ ,45689
+ ,382712
+ ,66
+ ,9616
+ ,242
+ ,0
+ ,4655
+ ,393
+ ,60720
+ ,341570
+ ,105
+ ,360
+ ,30
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,192
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,211
+ ,0
+ ,36
+ ,4
+ ,61564
+ ,312412
+ ,325
+ ,28103
+ ,3146
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,189
+ ,0
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,334
+ ,12177
+ ,541
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,205
+ ,0
+ ,126
+ ,10
+ ,60951
+ ,315547
+ ,171
+ ,11555
+ ,921
+ ,0
+ ,400
+ ,152
+ ,60720
+ ,313267
+ ,319
+ ,745
+ ,283
+ ,0
+ ,71
+ ,1
+ ,60745
+ ,316176
+ ,159
+ ,116176
+ ,1646
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,209
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,190
+ ,0
+ ,972
+ ,34
+ ,71642
+ ,359335
+ ,87
+ ,4686
+ ,164
+ ,0
+ ,531
+ ,10
+ ,71641
+ ,330068
+ ,120
+ ,13007
+ ,245
+ ,0
+ ,2461
+ ,57
+ ,55792
+ ,314289
+ ,310
+ ,2005
+ ,46
+ ,0
+ ,378
+ ,52
+ ,71873
+ ,297413
+ ,369
+ ,1873
+ ,258
+ ,0
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,306
+ ,22961
+ ,4974
+ ,0
+ ,638
+ ,14
+ ,60370
+ ,333210
+ ,114
+ ,9515
+ ,209
+ ,0
+ ,2300
+ ,29
+ ,64873
+ ,352108
+ ,97
+ ,5245
+ ,66
+ ,0
+ ,149
+ ,5
+ ,62041
+ ,313332
+ ,318
+ ,22666
+ ,761
+ ,0
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,382
+ ,18357
+ ,406
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,184
+ ,0
+ ,275
+ ,4
+ ,59500
+ ,318745
+ ,146
+ ,29686
+ ,431
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,179
+ ,0
+ ,141
+ ,6
+ ,61630
+ ,315366
+ ,295
+ ,19228
+ ,819
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,186
+ ,0
+ ,28
+ ,2
+ ,60890
+ ,315688
+ ,166
+ ,57844
+ ,4148
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,263
+ ,0
+ ,4980
+ ,91
+ ,113521
+ ,409642
+ ,55
+ ,2304
+ ,42
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,259
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,261
+ ,0
+ ,472
+ ,20
+ ,80045
+ ,269587
+ ,397
+ ,3479
+ ,148
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,243
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,250
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,247
+ ,0
+ ,203
+ ,27
+ ,50804
+ ,300962
+ ,361
+ ,3739
+ ,498
+ ,0
+ ,496
+ ,17
+ ,87390
+ ,325479
+ ,128
+ ,7381
+ ,253
+ ,0
+ ,10
+ ,2
+ ,61656
+ ,316155
+ ,160
+ ,58077
+ ,11639
+ ,0
+ ,63
+ ,4
+ ,65688
+ ,318574
+ ,148
+ ,29643
+ ,1891
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,228
+ ,0
+ ,1136
+ ,32
+ ,48522
+ ,343613
+ ,103
+ ,4488
+ ,126
+ ,0
+ ,265
+ ,31
+ ,60720
+ ,306948
+ ,346
+ ,3450
+ ,404
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,210
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,206
+ ,0
+ ,267
+ ,32
+ ,57640
+ ,330059
+ ,121
+ ,4064
+ ,487
+ ,0
+ ,474
+ ,20
+ ,61977
+ ,288985
+ ,386
+ ,4449
+ ,188
+ ,0
+ ,534
+ ,7
+ ,62620
+ ,304485
+ ,351
+ ,14926
+ ,196
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,201
+ ,0
+ ,15
+ ,8
+ ,60831
+ ,315688
+ ,164
+ ,14461
+ ,7687
+ ,0
+ ,397
+ ,28
+ ,60646
+ ,317736
+ ,152
+ ,4205
+ ,297
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,199
+ ,0
+ ,1866
+ ,29
+ ,56225
+ ,322331
+ ,135
+ ,4218
+ ,66
+ ,0
+ ,288
+ ,4
+ ,60510
+ ,296656
+ ,373
+ ,24164
+ ,335
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,198
+ ,0
+ ,3
+ ,2
+ ,60698
+ ,315354
+ ,296
+ ,57677
+ ,36505
+ ,0
+ ,468
+ ,21
+ ,60720
+ ,312161
+ ,328
+ ,5341
+ ,239
+ ,0
+ ,20
+ ,2
+ ,60805
+ ,315576
+ ,169
+ ,57788
+ ,5924
+ ,0
+ ,278
+ ,26
+ ,61404
+ ,314922
+ ,302
+ ,4420
+ ,413
+ ,0
+ ,61
+ ,14
+ ,60720
+ ,314551
+ ,307
+ ,8182
+ ,1874
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,215
+ ,0
+ ,192
+ ,4
+ ,65276
+ ,312339
+ ,327
+ ,28085
+ ,586
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,214
+ ,0
+ ,317
+ ,9
+ ,63915
+ ,298700
+ ,363
+ ,10967
+ ,311
+ ,0
+ ,738
+ ,10
+ ,60720
+ ,321376
+ ,139
+ ,12138
+ ,164
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,216
+ ,0
+ ,368
+ ,17
+ ,61686
+ ,303230
+ ,355
+ ,6072
+ ,281
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,204
+ ,0
+ ,2
+ ,1
+ ,60743
+ ,315487
+ ,172
+ ,115487
+ ,52257
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,202
+ ,0
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,163
+ ,19299
+ ,2197
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,195
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,194
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,200
+ ,0
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,321
+ ,37629
+ ,1200
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,237
+ ,0
+ ,24
+ ,8
+ ,59818
+ ,315637
+ ,168
+ ,14455
+ ,4772
+ ,0
+ ,2332
+ ,4
+ ,72680
+ ,324385
+ ,131
+ ,31096
+ ,53
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,236
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,233
+ ,0
+ ,131
+ ,11
+ ,61808
+ ,308989
+ ,338
+ ,9908
+ ,834
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,260
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,234
+ ,0
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,372
+ ,10745
+ ,469
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,240
+ ,0
+ ,167
+ ,2
+ ,64245
+ ,307322
+ ,345
+ ,53661
+ ,642
+ ,0
+ ,622
+ ,73
+ ,73007
+ ,304376
+ ,352
+ ,1430
+ ,168
+ ,0
+ ,2328
+ ,94
+ ,82732
+ ,253588
+ ,400
+ ,570
+ ,23
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,238
+ ,0
+ ,365
+ ,8
+ ,54820
+ ,309560
+ ,335
+ ,13695
+ ,300
+ ,0
+ ,364
+ ,35
+ ,47705
+ ,298466
+ ,365
+ ,2813
+ ,271
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,239
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,235
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,262
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,283
+ ,0
+ ,226
+ ,12
+ ,72835
+ ,343929
+ ,102
+ ,11994
+ ,636
+ ,0
+ ,307
+ ,15
+ ,58856
+ ,331955
+ ,116
+ ,8797
+ ,430
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,288
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,281
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,282
+ ,0
+ ,188
+ ,11
+ ,77655
+ ,381180
+ ,67
+ ,16471
+ ,964
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,221
+ ,0
+ ,138
+ ,6
+ ,69817
+ ,331420
+ ,117
+ ,21903
+ ,954
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,286
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,285
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,265
+ ,0
+ ,125
+ ,12
+ ,60798
+ ,310201
+ ,332
+ ,9183
+ ,883
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,284
+ ,0
+ ,282
+ ,30
+ ,62452
+ ,320016
+ ,141
+ ,4001
+ ,425
+ ,0
+ ,335
+ ,33
+ ,64175
+ ,320398
+ ,140
+ ,3648
+ ,359
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,232
+ ,0
+ ,1324
+ ,117
+ ,67440
+ ,291841
+ ,381
+ ,785
+ ,69
+ ,0
+ ,176
+ ,28
+ ,68136
+ ,310670
+ ,331
+ ,3952
+ ,629
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,256
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,255
+ ,0
+ ,249
+ ,72
+ ,56726
+ ,313491
+ ,317
+ ,1576
+ ,455
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,242
+ ,0
+ ,333
+ ,13
+ ,70811
+ ,331323
+ ,118
+ ,10102
+ ,394
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,241
+ ,0
+ ,601
+ ,6
+ ,60720
+ ,319210
+ ,144
+ ,19868
+ ,198
+ ,0
+ ,30
+ ,4
+ ,62045
+ ,318098
+ ,149
+ ,29525
+ ,3912
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,254
+ ,0
+ ,249
+ ,62
+ ,54323
+ ,292754
+ ,378
+ ,1496
+ ,373
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,258
+ ,0
+ ,165
+ ,24
+ ,62841
+ ,325176
+ ,130
+ ,5216
+ ,759
+ ,0
+ ,453
+ ,21
+ ,81125
+ ,365959
+ ,84
+ ,7903
+ ,366
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,257
+ ,0
+ ,53
+ ,14
+ ,59506
+ ,302409
+ ,356
+ ,7315
+ ,1932
+ ,0
+ ,382
+ ,21
+ ,59365
+ ,340968
+ ,107
+ ,6713
+ ,369
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,229
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,230
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,224
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,223
+ ,0
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,320
+ ,28291
+ ,3808
+ ,0
+ ,290
+ ,2
+ ,58790
+ ,301164
+ ,360
+ ,50582
+ ,348
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,225
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,231
+ ,0
+ ,366
+ ,53
+ ,61808
+ ,344425
+ ,101
+ ,2725
+ ,394
+ ,0
+ ,2
+ ,9
+ ,60735
+ ,315394
+ ,174
+ ,12822
+ ,63057
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,222
+ ,0
+ ,209
+ ,13
+ ,64016
+ ,316647
+ ,157
+ ,8973
+ ,558
+ ,0
+ ,384
+ ,22
+ ,54683
+ ,309836
+ ,333
+ ,4993
+ ,286
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,227
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,244
+ ,0
+ ,365
+ ,83
+ ,87192
+ ,346611
+ ,99
+ ,1766
+ ,402
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,226
+ ,0
+ ,49
+ ,8
+ ,64107
+ ,322031
+ ,137
+ ,15254
+ ,2500
+ ,0
+ ,3
+ ,4
+ ,60761
+ ,315656
+ ,167
+ ,28914
+ ,37920
+ ,0
+ ,133
+ ,14
+ ,65990
+ ,339445
+ ,109
+ ,9960
+ ,1052
+ ,0
+ ,32
+ ,1
+ ,59988
+ ,314964
+ ,301
+ ,114964
+ ,3603
+ ,0
+ ,368
+ ,17
+ ,61167
+ ,297141
+ ,370
+ ,5714
+ ,264
+ ,0
+ ,1
+ ,6
+ ,60719
+ ,315372
+ ,293
+ ,19229
+ ,84213
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,274
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,246
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,177
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,185
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,245
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,273
+ ,0
+ ,22
+ ,2
+ ,60722
+ ,312502
+ ,324
+ ,56251
+ ,5154
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,252
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,253
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,277
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,278
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,275
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,276
+ ,0
+ ,96
+ ,5
+ ,60379
+ ,313729
+ ,316
+ ,22746
+ ,1188
+ ,0
+ ,1
+ ,2
+ ,60727
+ ,315388
+ ,175
+ ,57694
+ ,117743
+ ,0
+ ,314
+ ,5
+ ,60720
+ ,315371
+ ,294
+ ,23074
+ ,367
+ ,0
+ ,844
+ ,78
+ ,60925
+ ,296139
+ ,375
+ ,1233
+ ,114
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,180
+ ,0
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,315
+ ,113880
+ ,4328
+ ,0
+ ,125
+ ,13
+ ,59734
+ ,317698
+ ,153
+ ,9054
+ ,944
+ ,0
+ ,304
+ ,15
+ ,62969
+ ,295580
+ ,376
+ ,6372
+ ,314
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,187
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,183
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,182
+ ,0
+ ,621
+ ,48
+ ,60720
+ ,308256
+ ,341
+ ,2255
+ ,174
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,181
+ ,0
+ ,119
+ ,6
+ ,59118
+ ,303677
+ ,354
+ ,17280
+ ,872
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,178
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,176
+ ,0
+ ,1595
+ ,17
+ ,60720
+ ,319369
+ ,143
+ ,7022
+ ,75
+ ,0
+ ,312
+ ,14
+ ,58598
+ ,318690
+ ,147
+ ,8478
+ ,381
+ ,0
+ ,60
+ ,10
+ ,61124
+ ,314049
+ ,314
+ ,11405
+ ,1916
+ ,0
+ ,587
+ ,12
+ ,59595
+ ,325699
+ ,126
+ ,10475
+ ,214
+ ,0
+ ,135
+ ,2
+ ,62065
+ ,314210
+ ,311
+ ,57105
+ ,847
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,272
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,269
+ ,0
+ ,514
+ ,52
+ ,78780
+ ,322378
+ ,133
+ ,2353
+ ,238
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,279
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,251
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,270
+ ,0
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,173
+ ,28849
+ ,172235
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,268
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,289
+ ,0
+ ,1763
+ ,24
+ ,61600
+ ,308336
+ ,340
+ ,4514
+ ,61
+ ,0
+ ,180
+ ,11
+ ,59635
+ ,316386
+ ,158
+ ,10581
+ ,647
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,287
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,248
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,291
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,292
+ ,0
+ ,218
+ ,21
+ ,60720
+ ,315553
+ ,170
+ ,5503
+ ,531
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,290
+ ,0
+ ,448
+ ,40
+ ,59781
+ ,323361
+ ,132
+ ,3084
+ ,275
+ ,0
+ ,227
+ ,9
+ ,76644
+ ,336639
+ ,110
+ ,15182
+ ,603
+ ,0
+ ,174
+ ,1
+ ,64820
+ ,307424
+ ,344
+ ,107424
+ ,618
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,266
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,280
+ ,0
+ ,121
+ ,24
+ ,56178
+ ,295370
+ ,377
+ ,3974
+ ,788
+ ,0
+ ,607
+ ,11
+ ,60436
+ ,322340
+ ,134
+ ,11122
+ ,201
+ ,0
+ ,2212
+ ,14
+ ,60720
+ ,319864
+ ,142
+ ,8562
+ ,54
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,271
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,264
+ ,0
+ ,530
+ ,60
+ ,73433
+ ,317291
+ ,155
+ ,1955
+ ,221
+ ,0
+ ,571
+ ,80
+ ,41477
+ ,280398
+ ,390
+ ,1005
+ ,141
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,267
+ ,0
+ ,78
+ ,16
+ ,62700
+ ,317330
+ ,154
+ ,7333
+ ,1511
+ ,0
+ ,2489
+ ,40
+ ,67804
+ ,238125
+ ,408
+ ,953
+ ,15
+ ,0
+ ,131
+ ,6
+ ,59661
+ ,327071
+ ,123
+ ,21179
+ ,969
+ ,0
+ ,923
+ ,8
+ ,58620
+ ,309038
+ ,337
+ ,13630
+ ,118
+ ,0
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,312
+ ,38070
+ ,1591
+ ,0
+ ,572
+ ,16
+ ,58580
+ ,307930
+ ,342
+ ,6746
+ ,189
+ ,-1
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,136
+ ,12233
+ ,308
+ ,-1
+ ,450
+ ,8
+ ,59325
+ ,292136
+ ,380
+ ,11517
+ ,205
+ ,-1
+ ,622
+ ,7
+ ,60950
+ ,263276
+ ,398
+ ,9039
+ ,102
+ ,-2
+ ,694
+ ,8
+ ,68060
+ ,367655
+ ,83
+ ,20957
+ ,241
+ ,-2
+ ,3425
+ ,12
+ ,83620
+ ,283910
+ ,388
+ ,6993
+ ,25
+ ,-3
+ ,562
+ ,13
+ ,58456
+ ,283587
+ ,389
+ ,6430
+ ,149
+ ,-3
+ ,4917
+ ,42
+ ,52811
+ ,243650
+ ,406
+ ,1039
+ ,9
+ ,-3
+ ,1442
+ ,118
+ ,121173
+ ,438493
+ ,49
+ ,2021
+ ,165
+ ,-3
+ ,529
+ ,9
+ ,63870
+ ,296261
+ ,374
+ ,10696
+ ,182
+ ,-3
+ ,2126
+ ,138
+ ,21001
+ ,230621
+ ,413
+ ,222
+ ,14
+ ,-3
+ ,1061
+ ,5
+ ,70415
+ ,304252
+ ,353
+ ,20850
+ ,98
+ ,-4
+ ,776
+ ,9
+ ,64230
+ ,333505
+ ,113
+ ,14834
+ ,172
+ ,-4
+ ,611
+ ,8
+ ,59190
+ ,296919
+ ,371
+ ,12115
+ ,159
+ ,-4
+ ,1526
+ ,25
+ ,69351
+ ,278990
+ ,391
+ ,3160
+ ,52
+ ,-4
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,392
+ ,10985
+ ,130
+ ,-4
+ ,1182
+ ,13
+ ,70694
+ ,327007
+ ,124
+ ,9770
+ ,107
+ ,-5
+ ,621
+ ,16
+ ,68005
+ ,317046
+ ,156
+ ,7315
+ ,188
+ ,-5
+ ,989
+ ,11
+ ,58930
+ ,304555
+ ,350
+ ,9505
+ ,106
+ ,-5
+ ,438
+ ,11
+ ,58320
+ ,298096
+ ,366
+ ,8918
+ ,224
+ ,-5
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,412
+ ,10620
+ ,44
+ ,-5
+ ,1303
+ ,61
+ ,69863
+ ,309422
+ ,336
+ ,1794
+ ,84
+ ,-5
+ ,7419
+ ,29
+ ,63255
+ ,286963
+ ,387
+ ,2999
+ ,12
+ ,-5
+ ,1164
+ ,17
+ ,57320
+ ,269753
+ ,395
+ ,4103
+ ,60
+ ,-6
+ ,3310
+ ,33
+ ,75230
+ ,448243
+ ,44
+ ,7523
+ ,75
+ ,-6
+ ,1920
+ ,15
+ ,79420
+ ,165404
+ ,424
+ ,-2306
+ ,-18
+ ,-6
+ ,965
+ ,3
+ ,73490
+ ,204325
+ ,421
+ ,1442
+ ,4
+ ,-8
+ ,3256
+ ,66
+ ,35250
+ ,407159
+ ,57
+ ,3139
+ ,64
+ ,-8
+ ,1135
+ ,17
+ ,62285
+ ,290476
+ ,384
+ ,5322
+ ,80
+ ,-9
+ ,1270
+ ,26
+ ,69206
+ ,275311
+ ,393
+ ,2897
+ ,59
+ ,-10
+ ,661
+ ,3
+ ,65920
+ ,246541
+ ,405
+ ,15514
+ ,70
+ ,-12
+ ,1013
+ ,2
+ ,69770
+ ,253468
+ ,402
+ ,26734
+ ,53
+ ,-14
+ ,2844
+ ,67
+ ,72683
+ ,240897
+ ,407
+ ,610
+ ,14
+ ,-15
+ ,11528
+ ,70
+ ,-14545
+ ,-83265
+ ,431
+ ,-4047
+ ,-25
+ ,-25
+ ,6526
+ ,26
+ ,55830
+ ,-42143
+ ,430
+ ,-9313
+ ,-37
+ ,-28
+ ,2264
+ ,24
+ ,55174
+ ,272713
+ ,394
+ ,3030
+ ,32
+ ,-31
+ ,5109
+ ,97
+ ,67038
+ ,215362
+ ,417
+ ,158
+ ,3
+ ,-31
+ ,3999
+ ,30
+ ,51252
+ ,42754
+ ,428
+ ,-5242
+ ,-39
+ ,-58
+ ,35624
+ ,223
+ ,157278
+ ,306275
+ ,347
+ ,477
+ ,3
+ ,-70
+ ,9252
+ ,48
+ ,79510
+ ,253537
+ ,401
+ ,1115
+ ,6
+ ,-129
+ ,15236
+ ,90
+ ,77440
+ ,372631
+ ,78
+ ,1918
+ ,11
+ ,-158
+ ,18073
+ ,180
+ ,27284
+ ,-7170
+ ,429
+ ,-1151
+ ,-11
+ ,-366)
+ ,dim=c(8
+ ,431)
+ ,dimnames=list(c('Costs'
+ ,'Trades'
+ ,'Dividends'
+ ,'Wealth'
+ ,'WRank'
+ ,'Profit/Trades'
+ ,'Profit/Cost'
+ ,'CCScore
')
+ ,1:431))
> y <- array(NA,dim=c(8,431),dimnames=list(c('Costs','Trades','Dividends','Wealth','WRank','Profit/Trades','Profit/Cost','CCScore
'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '8'
> ylab = ''
> xlab = ''
> main = ''
> #'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
CCScore\r Costs Trades Dividends Wealth WRank Profit/Trades Profit/Cost
1 56289 162556 1081 213118 6282154 1 5626 37
2 28328 29790 309 81767 4321023 2 13337 138
3 17936 87550 458 153198 4111912 3 8541 45
4 4145 84738 588 -26007 223193 415 39 0
5 3040 54660 302 126942 1491348 6 4276 24
6 2964 42634 156 157214 1629616 5 9164 34
7 2865 40949 481 129352 1398893 9 2493 29
8 2854 45187 353 234817 1926517 4 4891 38
9 2167 37704 452 60448 983660 14 1734 21
10 1974 16275 109 47818 1443586 7 11409 76
11 1910 25830 115 245546 1073089 10 7592 34
12 1871 12679 110 48020 984885 13 7135 62
13 943 18014 239 -1710 1405225 8 5043 67
14 929 43556 247 32648 227132 414 110 1
15 822 24811 505 95350 929118 15 1444 29
16 819 6575 159 151352 1071292 11 5480 133
17 769 7123 109 288170 638830 28 4026 62
18 745 21950 519 114337 856956 17 1266 30
19 652 37597 248 37884 992426 12 3195 21
20 643 17821 373 122844 444477 46 655 14
21 601 12988 119 82340 857217 16 5523 51
22 446 22330 84 79801 711969 20 6095 23
23 436 13326 102 165548 702380 21 4925 38
24 379 16189 295 116384 358589 89 538 10
25 305 7146 105 134028 297978 367 933 14
26 284 15824 64 63838 585715 30 6027 24
27 247 27664 282 74996 657954 24 1624 17
28 238 11920 182 31080 209458 419 52 1
29 223 8568 37 32168 786690 18 15856 68
30 220 14416 361 49857 439798 48 664 17
31 217 3369 28 87161 688779 23 17456 145
32 199 11819 85 106113 574339 31 4404 32
33 195 6984 45 80570 741409 19 12031 78
34 163 4519 49 102129 597793 29 8118 88
35 154 2220 22 301670 644190 26 20190 200
36 143 18562 155 102313 377934 72 1148 10
37 135 10327 91 88577 640273 27 4838 43
38 132 5336 81 112477 697458 22 6141 93
39 120 2365 79 191778 550608 32 4438 148
40 119 4069 145 79804 207393 420 51 2
41 108 8636 855 128294 301607 359 119 12
42 104 13718 61 96448 345783 100 2390 11
43 101 4525 226 93811 501749 37 1335 67
44 90 6869 105 117520 379983 68 1714 26
45 85 4628 62 69159 387475 64 3024 41
46 74 3689 25 101792 377305 74 7092 48
47 73 4891 217 210568 370837 79 787 35
48 70 7489 322 136996 430866 50 717 31
49 67 4901 84 121920 469107 41 3204 55
50 66 2284 33 76403 194493 423 -167 -2
51 66 3160 108 108094 530670 33 3062 105
52 66 4150 150 134759 518365 35 2122 77
53 59 7285 115 188873 491303 39 2533 40
54 58 1134 162 146216 527021 34 2019 288
55 58 4658 158 156608 233773 411 214 7
56 54 2384 97 61348 405972 58 2123 86
57 53 3748 9 50350 652925 25 50325 121
58 49 5371 66 87720 446211 45 3730 46
59 49 1285 107 99489 341340 106 1321 110
60 44 9327 101 87419 387699 63 1858 20
61 39 5565 47 94355 493408 38 6243 53
62 38 1528 38 60326 146494 426 -1408 -35
63 38 3122 34 94670 414462 54 6308 69
64 37 7561 87 82425 364304 86 1889 22
65 36 2675 79 59017 355178 93 1964 58
66 35 13253 947 90829 357760 90 167 12
67 34 880 74 80791 261216 399 827 70
68 33 2053 53 100423 397144 60 3720 96
69 32 1424 94 131116 374943 77 1861 123
70 32 4036 63 100269 424898 53 3570 56
71 31 3045 58 27330 202055 422 35 1
72 31 5119 49 39039 378525 69 3643 35
73 30 1431 34 106885 310768 330 3258 77
74 30 554 11 79285 325738 125 11431 227
75 30 1975 35 118881 394510 61 5557 98
76 28 1765 20 77623 247060 404 2353 27
77 26 1012 47 114768 368078 82 3576 166
78 26 810 43 74015 236761 409 855 45
79 26 1280 117 69465 312378 326 960 88
80 25 666 171 117869 339836 108 818 210
81 25 1380 26 60982 347385 98 5669 107
82 24 4677 75 90131 426280 52 3017 48
83 24 876 59 138971 352850 96 2591 174
84 23 814 18 39625 301881 358 5660 125
85 23 514 15 102725 377516 73 11834 346
86 22 5692 72 64239 357312 91 2185 28
87 22 3642 86 90262 458343 42 3004 71
88 22 540 14 103960 354228 94 11016 286
89 20 2099 64 106611 308636 339 1697 52
90 20 567 11 103345 386212 65 16928 329
91 19 2001 52 95551 393343 62 3718 97
92 18 2949 41 82903 378509 70 4354 61
93 17 2253 99 63593 452469 43 2550 112
94 16 6533 75 126910 364839 85 2198 25
95 14 1889 45 37527 358649 88 3526 84
96 14 3055 43 60247 376641 75 4108 58
97 14 272 8 112995 429112 51 28639 844
98 13 1414 198 70184 330546 119 659 92
99 13 2564 22 130140 403560 59 9253 79
100 13 1383 11 73221 317892 151 10717 85
101 13 1261 33 76114 307528 343 3258 85
102 13 975 23 90534 235133 410 1528 36
103 12 3366 80 108479 299243 362 1241 29
104 12 576 18 113761 314073 313 6337 198
105 11 1686 40 68696 368186 81 4205 100
106 10 746 23 71561 269661 396 3029 93
107 10 3192 60 59831 125390 427 -1243 -23
108 10 2045 20 97890 510834 36 15542 152
109 10 5702 61 101481 321896 138 1998 21
110 10 1932 36 72954 249898 403 1386 26
111 9 936 30 67939 408881 56 6963 223
112 9 3437 47 48022 158492 425 -883 -12
113 9 5131 71 86111 292154 379 1298 18
114 6 2397 14 74020 289513 385 6394 37
115 6 1389 9 57530 378049 71 19783 128
116 6 1503 39 56364 343466 104 3679 95
117 5 402 26 84990 332743 115 5106 330
118 5 2239 21 88590 442882 47 11566 108
119 5 2234 16 77200 214215 418 888 6
120 4 837 69 61262 315688 165 1677 138
121 4 10579 92 110309 375195 76 1904 17
122 3 875 14 67000 334280 112 9591 153
123 2 1585 107 93099 355864 92 1457 98
124 2 1659 29 107577 480382 40 9668 169
125 1 2647 37 62920 353058 95 4137 58
126 1 3294 23 75832 217193 416 748 5
127 7 0 0 60720 315380 191 0 94
128 28 60793 314533 308 16362 1212 0 422
129 0 57935 318056 150 4216 280 0 0
130 314353 60720 315380 197 0 34 8 60630
131 369448 309 14294 3340 0 1558 63 55637
132 315380 80 2690 109 0 0 0 60720
133 37615 196 0 43 3 60887 312846 323
134 22415 2613 0 645 5 60720 312075 329
135 12779 174 0 316 9 60505 315009 299
136 9146 364 0 115 13 60945 318903 145
137 57443 1037 0 5 2 60720 314887 304
138 22983 24814 0 897 5 60720 314913 303
139 0 128 0 0 0 60720 315380 188
140 0 389 14 58990 325506 127 8965 322
141 15 0 0 60720 315380 220 0 1002
142 3 56750 298568 364 6571 98 0 36
143 15 60894 315834 162 38611 3215 0 460
144 11 63346 329784 122 8652 282 0 309
145 0 56535 312878 322 10262 365 0 0
146 314987 60720 315380 203 0 9 6 60835
147 325249 300 19165 12458 0 271 2 60720
148 315877 129 62625 463 0 14 1 61016
149 291650 161 115877 8041 0 520 10 58650
150 305959 383 9165 176 0 1766 73 60438
151 315380 349 1451 60 0 0 0 60720
152 8888 208 0 458 11 58625 297765 368
153 38415 213 0 20 3 60938 315245 297
154 0 5674 0 0 0 60720 315380 217
155 2 0 0 60720 315380 219 0 98
156 7 61490 315236 298 57618 1179 0 405
157 0 60845 336425 111 19489 337 0 0
158 315380 60720 315380 213 0 0 0 60720
159 0 212 0 0 0 60720 315380 207
160 27 0 0 60720 315380 218 0 483
161 51 60830 306268 348 3936 220 0 454
162 3 63261 302187 357 2004 225 0 47
163 0 60720 314882 305 38294 2432 0 0
164 382712 60720 315380 193 0 757 19 45689
165 341570 66 9616 242 0 4655 393 60720
166 315380 105 360 30 0 0 0 60720
167 0 192 0 0 0 60720 315380 211
168 0 36 4 61564 312412 325 28103 3146
169 9 0 0 60720 315380 189 0 203
170 0 61938 309596 334 12177 541 0 0
171 315547 60720 315380 205 0 126 10 60951
172 313267 171 11555 921 0 400 152 60720
173 316176 319 745 283 0 71 1 60745
174 315380 159 116176 1646 0 0 0 60720
175 0 209 0 0 0 60720 315380 190
176 0 972 34 71642 359335 87 4686 164
177 0 531 10 71641 330068 120 13007 245
178 0 2461 57 55792 314289 310 2005 46
179 0 378 52 71873 297413 369 1873 258
180 0 23 5 62555 314806 306 22961 4974
181 0 638 14 60370 333210 114 9515 209
182 0 2300 29 64873 352108 97 5245 66
183 0 149 5 62041 313332 318 22666 761
184 0 226 5 65745 291787 382 18357 406
185 4 0 0 60720 315380 184 0 275
186 0 59500 318745 146 29686 431 0 0
187 315366 60720 315380 179 0 141 6 61630
188 315380 295 19228 819 0 0 0 60720
189 57844 186 0 28 2 60890 315688 166
190 0 4148 0 0 0 60720 315380 263
191 0 4980 91 113521 409642 55 2304 42
192 0 0 0 60720 315380 259 0 0
193 269587 60720 315380 261 0 472 20 80045
194 315380 397 3479 148 0 0 0 60720
195 0 243 0 0 0 60720 315380 250
196 27 0 0 60720 315380 247 0 203
197 17 50804 300962 361 3739 498 0 496
198 2 87390 325479 128 7381 253 0 10
199 4 61656 316155 160 58077 11639 0 63
200 0 65688 318574 148 29643 1891 0 0
201 343613 60720 315380 228 0 1136 32 48522
202 306948 103 4488 126 0 265 31 60720
203 210 346 3450 404 0 0 60720 315380
204 267 0 0 0 60720 315380 206 0
205 474 32 57640 330059 121 4064 487 0
206 534 20 61977 288985 386 4449 188 0
207 0 7 62620 304485 351 14926 196 0
208 60831 0 60720 315380 201 0 15 8
209 60646 315688 164 14461 7687 0 397 28
210 60720 317736 152 4205 297 0 0 0
211 135 315380 199 0 1866 29 56225 322331
212 373 4218 66 0 288 4 60510 296656
213 198 24164 335 0 0 0 60720 315380
214 36505 0 3 2 60698 315354 296 57677
215 239 0 468 21 60720 312161 328 5341
216 5924 0 20 2 60805 315576 169 57788
217 413 0 278 26 61404 314922 302 4420
218 1874 0 61 14 60720 314551 307 8182
219 192 0 0 0 60720 315380 215 0
220 0 4 65276 312339 327 28085 586 0
221 63915 0 60720 315380 214 0 317 9
222 60720 298700 363 10967 311 0 738 10
223 60720 321376 139 12138 164 0 0 0
224 355 315380 216 0 368 17 61686 303230
225 204 6072 281 0 0 0 60720 315380
226 52257 0 2 1 60743 315487 172 115487
227 53 0 0 0 60720 315380 202 0
228 0 6 60349 315793 163 19299 2197 0
229 60720 0 60720 315380 195 0 0 0
230 200 315380 194 0 0 0 60720 315380
231 1200 0 94 3 61360 312887 321 37629
232 24 0 0 0 60720 315380 237 0
233 2332 8 59818 315637 168 14455 4772 0
234 0 4 72680 324385 131 31096 53 0
235 60720 0 60720 315380 236 0 0 0
236 338 315380 233 0 131 11 61808 308989
237 260 9908 834 0 0 0 60720 315380
238 206 0 0 0 60720 315380 234 0
239 0 9 53110 296702 372 10745 469 0
240 64245 0 60720 315380 240 0 167 2
241 73007 307322 345 53661 642 0 622 73
242 82732 304376 352 1430 168 0 2328 94
243 60720 253588 400 570 23 0 0 0
244 335 315380 238 0 365 8 54820 309560
245 365 13695 300 0 364 35 47705 298466
246 239 2813 271 0 0 0 60720 315380
247 0 0 0 0 60720 315380 235 0
248 60720 0 60720 315380 262 0 0 0
249 102 315380 283 0 226 12 72835 343929
250 116 11994 636 0 307 15 58856 331955
251 288 8797 430 0 0 0 60720 315380
252 0 0 0 0 60720 315380 281 0
253 77655 0 60720 315380 282 0 188 11
254 60720 381180 67 16471 964 0 0 0
255 117 315380 221 0 138 6 69817 331420
256 286 21903 954 0 0 0 60720 315380
257 0 0 0 0 60720 315380 285 0
258 60798 0 60720 315380 265 0 125 12
259 60720 310201 332 9183 883 0 0 0
260 141 315380 284 0 282 30 62452 320016
261 140 4001 425 0 335 33 64175 320398
262 232 3648 359 0 0 0 60720 315380
263 69 0 1324 117 67440 291841 381 785
264 629 0 176 28 68136 310670 331 3952
265 0 0 0 0 60720 315380 256 0
266 56726 0 60720 315380 255 0 249 72
267 60720 313491 317 1576 455 0 0 0
268 118 315380 242 0 333 13 70811 331323
269 241 10102 394 0 0 0 60720 315380
270 198 0 601 6 60720 319210 144 19868
271 3912 0 30 4 62045 318098 149 29525
272 249 0 0 0 60720 315380 254 0
273 0 62 54323 292754 378 1496 373 0
274 62841 0 60720 315380 258 0 165 24
275 81125 325176 130 5216 759 0 453 21
276 60720 365959 84 7903 366 0 0 0
277 356 315380 257 0 53 14 59506 302409
278 107 7315 1932 0 382 21 59365 340968
279 229 6713 369 0 0 0 60720 315380
280 0 0 0 0 60720 315380 230 0
281 60720 0 60720 315380 224 0 0 0
282 320 315380 223 0 30 4 60798 313164
283 360 28291 3808 0 290 2 58790 301164
284 225 50582 348 0 0 0 60720 315380
285 366 0 0 0 60720 315380 231 0
286 2 53 61808 344425 101 2725 394 0
287 0 9 60735 315394 174 12822 63057 0
288 64016 0 60720 315380 222 0 209 13
289 54683 316647 157 8973 558 0 384 22
290 60720 309836 333 4993 286 0 0 0
291 244 315380 227 0 0 0 60720 315380
292 402 0 365 83 87192 346611 99 1766
293 49 0 0 0 60720 315380 226 0
294 3 8 64107 322031 137 15254 2500 0
295 133 4 60761 315656 167 28914 37920 0
296 32 14 65990 339445 109 9960 1052 0
297 368 1 59988 314964 301 114964 3603 0
298 1 17 61167 297141 370 5714 264 0
299 0 6 60719 315372 293 19229 84213 0
300 60720 0 60720 315380 274 0 0 0
301 177 315380 246 0 0 0 60720 315380
302 0 0 0 0 60720 315380 185 0
303 60720 0 60720 315380 245 0 0 0
304 324 315380 273 0 22 2 60722 312502
305 252 56251 5154 0 0 0 60720 315380
306 0 0 0 0 60720 315380 249 0
307 60720 0 60720 315380 253 0 0 0
308 278 315380 277 0 0 0 60720 315380
309 0 0 0 0 60720 315380 275 0
310 60379 0 60720 315380 276 0 96 5
311 60727 313729 316 22746 1188 0 1 2
312 60720 315388 175 57694 117743 0 314 5
313 60925 315371 294 23074 367 0 844 78
314 60720 296139 375 1233 114 0 0 0
315 315 315380 180 0 26 1 60896 313880
316 153 113880 4328 0 125 13 59734 317698
317 376 9054 944 0 304 15 62969 295580
318 187 6372 314 0 0 0 60720 315380
319 0 0 0 0 60720 315380 183 0
320 60720 0 60720 315380 182 0 621 48
321 60720 308256 341 2255 174 0 0 0
322 354 315380 181 0 119 6 59118 303677
323 178 17280 872 0 0 0 60720 315380
324 1595 0 0 0 60720 315380 176 0
325 312 17 60720 319369 143 7022 75 0
326 60 14 58598 318690 147 8478 381 0
327 587 10 61124 314049 314 11405 1916 0
328 135 12 59595 325699 126 10475 214 0
329 0 2 62065 314210 311 57105 847 0
330 60720 0 60720 315380 272 0 0 0
331 133 315380 269 0 514 52 78780 322378
332 279 2353 238 0 0 0 60720 315380
333 0 0 0 0 60720 315380 251 0
334 60722 0 60720 315380 270 0 1 4
335 60720 315398 173 28849 172235 0 0 0
336 289 315380 268 0 0 0 60720 315380
337 61 0 1763 24 61600 308336 340 4514
338 647 0 180 11 59635 316386 158 10581
339 0 0 0 0 60720 315380 287 0
340 60720 0 60720 315380 248 0 0 0
341 292 315380 291 0 0 0 60720 315380
342 531 0 218 21 60720 315553 170 5503
343 448 0 0 0 60720 315380 290 0
344 227 40 59781 323361 132 3084 275 0
345 174 9 76644 336639 110 15182 603 0
346 0 1 64820 307424 344 107424 618 0
347 60720 0 60720 315380 266 0 0 0
348 377 315380 280 0 121 24 56178 295370
349 134 3974 788 0 607 11 60436 322340
350 142 11122 201 0 2212 14 60720 319864
351 271 8562 54 0 0 0 60720 315380
352 530 0 0 0 60720 315380 264 0
353 571 60 73433 317291 155 1955 221 0
354 0 80 41477 280398 390 1005 141 0
355 62700 0 60720 315380 267 0 78 16
356 67804 317330 154 7333 1511 0 2489 40
357 59661 238125 408 953 15 0 131 6
358 58620 327071 123 21179 969 0 923 8
359 60398 309038 337 13630 118 0 72 3
360 58580 314210 312 38070 1591 0 572 16
361 62710 307930 342 6746 189 -1 397 10
362 59325 322327 136 12233 308 -1 450 8
363 60950 292136 380 11517 205 -1 622 7
364 68060 263276 398 9039 102 -2 694 8
365 83620 367655 83 20957 241 -2 3425 12
366 58456 283910 388 6993 25 -3 562 13
367 52811 283587 389 6430 149 -3 4917 42
368 121173 243650 406 1039 9 -3 1442 118
369 63870 438493 49 2021 165 -3 529 9
370 21001 296261 374 10696 182 -3 2126 138
371 70415 230621 413 222 14 -3 1061 5
372 64230 304252 353 20850 98 -4 776 9
373 59190 333505 113 14834 172 -4 611 8
374 69351 296919 371 12115 159 -4 1526 25
375 64270 278990 391 3160 52 -4 592 7
376 70694 276898 392 10985 130 -4 1182 13
377 68005 327007 124 9770 107 -5 621 16
378 58930 317046 156 7315 188 -5 989 11
379 58320 304555 350 9505 106 -5 438 11
380 69980 298096 366 8918 224 -5 726 3
381 69863 231861 412 10620 44 -5 1303 61
382 63255 309422 336 1794 84 -5 7419 29
383 57320 286963 387 2999 12 -5 1164 17
384 75230 269753 395 4103 60 -6 3310 33
385 79420 448243 44 7523 75 -6 1920 15
386 73490 165404 424 -2306 -18 -6 965 3
387 35250 204325 421 1442 4 -8 3256 66
388 62285 407159 57 3139 64 -8 1135 17
389 69206 290476 384 5322 80 -9 1270 26
390 65920 275311 393 2897 59 -10 661 3
391 69770 246541 405 15514 70 -12 1013 2
392 72683 253468 402 26734 53 -14 2844 67
393 -14545 240897 407 610 14 -15 11528 70
394 55830 -83265 431 -4047 -25 -25 6526 26
395 55174 -42143 430 -9313 -37 -28 2264 24
396 67038 272713 394 3030 32 -31 5109 97
397 51252 215362 417 158 3 -31 3999 30
398 157278 42754 428 -5242 -39 -58 35624 223
399 79510 306275 347 477 3 -70 9252 48
400 77440 253537 401 1115 6 -129 15236 90
401 27284 372631 78 1918 11 -158 18073 180
402 213118 -7170 429 -1151 -11 -366 162556 1081
403 81767 6282154 1 5626 37 56289 29790 309
404 153198 4321023 2 13337 138 28328 87550 458
405 -26007 4111912 3 8541 45 17936 84738 588
406 126942 223193 415 39 0 4145 54660 302
407 157214 1491348 6 4276 24 3040 42634 156
408 129352 1629616 5 9164 34 2964 40949 481
409 234817 1398893 9 2493 29 2865 45187 353
410 60448 1926517 4 4891 38 2854 37704 452
411 47818 983660 14 1734 21 2167 16275 109
412 245546 1443586 7 11409 76 1974 25830 115
413 48020 1073089 10 7592 34 1910 12679 110
414 -1710 984885 13 7135 62 1871 18014 239
415 32648 1405225 8 5043 67 943 43556 247
416 95350 227132 414 110 1 929 24811 505
417 151352 929118 15 1444 29 822 6575 159
418 288170 1071292 11 5480 133 819 7123 109
419 114337 638830 28 4026 62 769 21950 519
420 37884 856956 17 1266 30 745 37597 248
421 122844 992426 12 3195 21 652 17821 373
422 82340 444477 46 655 14 643 12988 119
423 79801 857217 16 5523 51 601 22330 84
424 165548 711969 20 6095 23 446 13326 102
425 116384 702380 21 4925 38 436 16189 295
426 134028 358589 89 538 10 379 7146 105
427 63838 297978 367 933 14 305 15824 64
428 74996 585715 30 6027 24 284 27664 282
429 31080 657954 24 1624 17 247 11920 182
430 32168 209458 419 52 1 238 8568 37
431 49857 786690 18 15856 68 223 14416 361
t
1 1
2 2
3 3
4 4
5 5
6 6
7 7
8 8
9 9
10 10
11 11
12 12
13 13
14 14
15 15
16 16
17 17
18 18
19 19
20 20
21 21
22 22
23 23
24 24
25 25
26 26
27 27
28 28
29 29
30 30
31 31
32 32
33 33
34 34
35 35
36 36
37 37
38 38
39 39
40 40
41 41
42 42
43 43
44 44
45 45
46 46
47 47
48 48
49 49
50 50
51 51
52 52
53 53
54 54
55 55
56 56
57 57
58 58
59 59
60 60
61 61
62 62
63 63
64 64
65 65
66 66
67 67
68 68
69 69
70 70
71 71
72 72
73 73
74 74
75 75
76 76
77 77
78 78
79 79
80 80
81 81
82 82
83 83
84 84
85 85
86 86
87 87
88 88
89 89
90 90
91 91
92 92
93 93
94 94
95 95
96 96
97 97
98 98
99 99
100 100
101 101
102 102
103 103
104 104
105 105
106 106
107 107
108 108
109 109
110 110
111 111
112 112
113 113
114 114
115 115
116 116
117 117
118 118
119 119
120 120
121 121
122 122
123 123
124 124
125 125
126 126
127 127
128 128
129 129
130 130
131 131
132 132
133 133
134 134
135 135
136 136
137 137
138 138
139 139
140 140
141 141
142 142
143 143
144 144
145 145
146 146
147 147
148 148
149 149
150 150
151 151
152 152
153 153
154 154
155 155
156 156
157 157
158 158
159 159
160 160
161 161
162 162
163 163
164 164
165 165
166 166
167 167
168 168
169 169
170 170
171 171
172 172
173 173
174 174
175 175
176 176
177 177
178 178
179 179
180 180
181 181
182 182
183 183
184 184
185 185
186 186
187 187
188 188
189 189
190 190
191 191
192 192
193 193
194 194
195 195
196 196
197 197
198 198
199 199
200 200
201 201
202 202
203 203
204 204
205 205
206 206
207 207
208 208
209 209
210 210
211 211
212 212
213 213
214 214
215 215
216 216
217 217
218 218
219 219
220 220
221 221
222 222
223 223
224 224
225 225
226 226
227 227
228 228
229 229
230 230
231 231
232 232
233 233
234 234
235 235
236 236
237 237
238 238
239 239
240 240
241 241
242 242
243 243
244 244
245 245
246 246
247 247
248 248
249 249
250 250
251 251
252 252
253 253
254 254
255 255
256 256
257 257
258 258
259 259
260 260
261 261
262 262
263 263
264 264
265 265
266 266
267 267
268 268
269 269
270 270
271 271
272 272
273 273
274 274
275 275
276 276
277 277
278 278
279 279
280 280
281 281
282 282
283 283
284 284
285 285
286 286
287 287
288 288
289 289
290 290
291 291
292 292
293 293
294 294
295 295
296 296
297 297
298 298
299 299
300 300
301 301
302 302
303 303
304 304
305 305
306 306
307 307
308 308
309 309
310 310
311 311
312 312
313 313
314 314
315 315
316 316
317 317
318 318
319 319
320 320
321 321
322 322
323 323
324 324
325 325
326 326
327 327
328 328
329 329
330 330
331 331
332 332
333 333
334 334
335 335
336 336
337 337
338 338
339 339
340 340
341 341
342 342
343 343
344 344
345 345
346 346
347 347
348 348
349 349
350 350
351 351
352 352
353 353
354 354
355 355
356 356
357 357
358 358
359 359
360 360
361 361
362 362
363 363
364 364
365 365
366 366
367 367
368 368
369 369
370 370
371 371
372 372
373 373
374 374
375 375
376 376
377 377
378 378
379 379
380 380
381 381
382 382
383 383
384 384
385 385
386 386
387 387
388 388
389 389
390 390
391 391
392 392
393 393
394 394
395 395
396 396
397 397
398 398
399 399
400 400
401 401
402 402
403 403
404 404
405 405
406 406
407 407
408 408
409 409
410 410
411 411
412 412
413 413
414 414
415 415
416 416
417 417
418 418
419 419
420 420
421 421
422 422
423 423
424 424
425 425
426 426
427 427
428 428
429 429
430 430
431 431
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Trades Dividends
3.427e+04 3.277e-03 2.111e-01 -2.005e-01
Wealth WRank `Profit/Trades` `Profit/Cost`
-7.774e-03 -2.020e-01 -1.109e-01 -1.258e-01
t
1.263e+02
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-128142.7 -26899.9 -15305.5 -357.2 323605.2
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.427e+04 9.895e+03 3.463 0.000589 ***
Costs 3.277e-03 8.093e-03 0.405 0.685695
Trades 2.111e-01 4.739e-02 4.455 1.08e-05 ***
Dividends -2.005e-01 3.731e-02 -5.375 1.27e-07 ***
Wealth -7.774e-03 8.408e-03 -0.925 0.355678
WRank -2.020e-01 4.140e-02 -4.878 1.52e-06 ***
`Profit/Trades` -1.109e-01 5.838e-02 -1.900 0.058075 .
`Profit/Cost` -1.258e-01 3.806e-02 -3.304 0.001034 **
t 1.263e+02 3.541e+01 3.567 0.000403 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 70400 on 422 degrees of freedom
Multiple R-squared: 0.1984, Adjusted R-squared: 0.1832
F-statistic: 13.06 on 8 and 422 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,] 5.904270e-04 1.180854e-03 9.994096e-01
[2,] 4.359540e-05 8.719080e-05 9.999564e-01
[3,] 4.923609e-06 9.847217e-06 9.999951e-01
[4,] 3.635685e-07 7.271370e-07 9.999996e-01
[5,] 3.342353e-08 6.684706e-08 1.000000e+00
[6,] 2.949171e-09 5.898341e-09 1.000000e+00
[7,] 2.250026e-10 4.500052e-10 1.000000e+00
[8,] 1.350280e-11 2.700560e-11 1.000000e+00
[9,] 8.972003e-13 1.794401e-12 1.000000e+00
[10,] 5.015618e-14 1.003124e-13 1.000000e+00
[11,] 2.739866e-15 5.479732e-15 1.000000e+00
[12,] 1.412347e-16 2.824694e-16 1.000000e+00
[13,] 7.052698e-18 1.410540e-17 1.000000e+00
[14,] 3.810817e-19 7.621634e-19 1.000000e+00
[15,] 1.757920e-20 3.515840e-20 1.000000e+00
[16,] 9.946025e-22 1.989205e-21 1.000000e+00
[17,] 5.979291e-23 1.195858e-22 1.000000e+00
[18,] 8.085465e-24 1.617093e-23 1.000000e+00
[19,] 3.725087e-25 7.450174e-25 1.000000e+00
[20,] 2.919857e-26 5.839715e-26 1.000000e+00
[21,] 1.306371e-27 2.612742e-27 1.000000e+00
[22,] 6.764805e-29 1.352961e-28 1.000000e+00
[23,] 3.091844e-30 6.183688e-30 1.000000e+00
[24,] 2.165589e-31 4.331177e-31 1.000000e+00
[25,] 9.174534e-33 1.834907e-32 1.000000e+00
[26,] 4.421312e-34 8.842625e-34 1.000000e+00
[27,] 2.531413e-35 5.062826e-35 1.000000e+00
[28,] 1.156928e-36 2.313856e-36 1.000000e+00
[29,] 9.675211e-38 1.935042e-37 1.000000e+00
[30,] 2.178948e-38 4.357896e-38 1.000000e+00
[31,] 9.115287e-40 1.823057e-39 1.000000e+00
[32,] 4.143779e-41 8.287558e-41 1.000000e+00
[33,] 1.674656e-42 3.349312e-42 1.000000e+00
[34,] 6.678943e-44 1.335789e-43 1.000000e+00
[35,] 2.629583e-45 5.259165e-45 1.000000e+00
[36,] 1.040641e-46 2.081283e-46 1.000000e+00
[37,] 4.410589e-48 8.821178e-48 1.000000e+00
[38,] 1.839064e-49 3.678128e-49 1.000000e+00
[39,] 1.027481e-50 2.054963e-50 1.000000e+00
[40,] 5.649760e-52 1.129952e-51 1.000000e+00
[41,] 2.728232e-53 5.456464e-53 1.000000e+00
[42,] 1.173487e-54 2.346973e-54 1.000000e+00
[43,] 7.717151e-56 1.543430e-55 1.000000e+00
[44,] 5.045944e-57 1.009189e-56 1.000000e+00
[45,] 1.907477e-58 3.814954e-58 1.000000e+00
[46,] 1.039713e-59 2.079427e-59 1.000000e+00
[47,] 4.052644e-61 8.105288e-61 1.000000e+00
[48,] 1.541825e-62 3.083650e-62 1.000000e+00
[49,] 5.930796e-64 1.186159e-63 1.000000e+00
[50,] 2.451766e-65 4.903533e-65 1.000000e+00
[51,] 1.338843e-66 2.677687e-66 1.000000e+00
[52,] 4.845142e-68 9.690285e-68 1.000000e+00
[53,] 1.795200e-69 3.590399e-69 1.000000e+00
[54,] 6.509314e-71 1.301863e-70 1.000000e+00
[55,] 4.204360e-72 8.408720e-72 1.000000e+00
[56,] 3.030240e-73 6.060479e-73 1.000000e+00
[57,] 1.094466e-74 2.188933e-74 1.000000e+00
[58,] 4.058736e-76 8.117471e-76 1.000000e+00
[59,] 1.489497e-77 2.978994e-77 1.000000e+00
[60,] 9.179564e-79 1.835913e-78 1.000000e+00
[61,] 3.163016e-80 6.326032e-80 1.000000e+00
[62,] 1.725395e-81 3.450791e-81 1.000000e+00
[63,] 6.309589e-83 1.261918e-82 1.000000e+00
[64,] 2.192305e-84 4.384611e-84 1.000000e+00
[65,] 1.177396e-85 2.354793e-85 1.000000e+00
[66,] 4.204348e-87 8.408696e-87 1.000000e+00
[67,] 2.240509e-88 4.481019e-88 1.000000e+00
[68,] 1.164902e-89 2.329804e-89 1.000000e+00
[69,] 4.320961e-91 8.641922e-91 1.000000e+00
[70,] 1.434755e-92 2.869510e-92 1.000000e+00
[71,] 5.203801e-94 1.040760e-93 1.000000e+00
[72,] 1.789184e-95 3.578369e-95 1.000000e+00
[73,] 8.833761e-97 1.766752e-96 1.000000e+00
[74,] 2.887147e-98 5.774293e-98 1.000000e+00
[75,] 9.947777e-100 1.989555e-99 1.000000e+00
[76,] 3.821275e-101 7.642549e-101 1.000000e+00
[77,] 1.219902e-102 2.439804e-102 1.000000e+00
[78,] 6.302632e-104 1.260526e-103 1.000000e+00
[79,] 1.992380e-105 3.984760e-105 1.000000e+00
[80,] 6.536001e-107 1.307200e-106 1.000000e+00
[81,] 2.094131e-108 4.188262e-108 1.000000e+00
[82,] 7.611727e-110 1.522345e-109 1.000000e+00
[83,] 2.597137e-111 5.194273e-111 1.000000e+00
[84,] 8.044313e-113 1.608863e-112 1.000000e+00
[85,] 2.498742e-114 4.997483e-114 1.000000e+00
[86,] 8.701460e-116 1.740292e-115 1.000000e+00
[87,] 2.932462e-117 5.864924e-117 1.000000e+00
[88,] 9.634789e-119 1.926958e-118 1.000000e+00
[89,] 2.961227e-120 5.922453e-120 1.000000e+00
[90,] 1.525570e-121 3.051140e-121 1.000000e+00
[91,] 7.873043e-123 1.574609e-122 1.000000e+00
[92,] 4.279264e-124 8.558527e-124 1.000000e+00
[93,] 1.782467e-125 3.564933e-125 1.000000e+00
[94,] 5.348421e-127 1.069684e-126 1.000000e+00
[95,] 2.484215e-128 4.968430e-128 1.000000e+00
[96,] 1.368629e-129 2.737258e-129 1.000000e+00
[97,] 5.785746e-131 1.157149e-130 1.000000e+00
[98,] 1.818864e-132 3.637728e-132 1.000000e+00
[99,] 8.254031e-134 1.650806e-133 1.000000e+00
[100,] 2.551692e-135 5.103385e-135 1.000000e+00
[101,] 1.350277e-136 2.700553e-136 1.000000e+00
[102,] 6.551344e-138 1.310269e-137 1.000000e+00
[103,] 2.616008e-139 5.232017e-139 1.000000e+00
[104,] 7.205877e-141 1.441175e-140 1.000000e+00
[105,] 2.085084e-142 4.170168e-142 1.000000e+00
[106,] 6.403349e-144 1.280670e-143 1.000000e+00
[107,] 2.036021e-145 4.072042e-145 1.000000e+00
[108,] 8.815361e-147 1.763072e-146 1.000000e+00
[109,] 2.744402e-148 5.488804e-148 1.000000e+00
[110,] 8.778013e-150 1.755603e-149 1.000000e+00
[111,] 2.410712e-151 4.821423e-151 1.000000e+00
[112,] 7.014048e-153 1.402810e-152 1.000000e+00
[113,] 2.750677e-154 5.501354e-154 1.000000e+00
[114,] 7.501974e-156 1.500395e-155 1.000000e+00
[115,] 3.260762e-157 6.521524e-157 1.000000e+00
[116,] 9.634362e-159 1.926872e-158 1.000000e+00
[117,] 1.416128e-159 2.832256e-159 1.000000e+00
[118,] 7.668229e-160 1.533646e-159 1.000000e+00
[119,] 8.590370e-160 1.718074e-159 1.000000e+00
[120,] 4.315460e-123 8.630920e-123 1.000000e+00
[121,] 4.905942e-116 9.811885e-116 1.000000e+00
[122,] 8.115143e-101 1.623029e-100 1.000000e+00
[123,] 4.644884e-101 9.289769e-101 1.000000e+00
[124,] 4.794334e-101 9.588668e-101 1.000000e+00
[125,] 2.646195e-101 5.292390e-101 1.000000e+00
[126,] 4.326771e-97 8.653542e-97 1.000000e+00
[127,] 6.056449e-98 1.211290e-97 1.000000e+00
[128,] 2.518146e-96 5.036292e-96 1.000000e+00
[129,] 2.740599e-97 5.481198e-97 1.000000e+00
[130,] 5.310676e-98 1.062135e-97 1.000000e+00
[131,] 2.651619e-98 5.303239e-98 1.000000e+00
[132,] 6.005093e-99 1.201019e-98 1.000000e+00
[133,] 2.392399e-99 4.784798e-99 1.000000e+00
[134,] 1.171499e-99 2.342998e-99 1.000000e+00
[135,] 1.959053e-98 3.918105e-98 1.000000e+00
[136,] 2.727530e-98 5.455060e-98 1.000000e+00
[137,] 1.851847e-97 3.703693e-97 1.000000e+00
[138,] 1.023800e-95 2.047600e-95 1.000000e+00
[139,] 1.931111e-94 3.862221e-94 1.000000e+00
[140,] 2.704535e-94 5.409070e-94 1.000000e+00
[141,] 1.777910e-94 3.555820e-94 1.000000e+00
[142,] 1.402033e-94 2.804066e-94 1.000000e+00
[143,] 7.190115e-94 1.438023e-93 1.000000e+00
[144,] 9.495889e-95 1.899178e-94 1.000000e+00
[145,] 3.010480e-95 6.020960e-95 1.000000e+00
[146,] 1.647604e-95 3.295208e-95 1.000000e+00
[147,] 2.321440e-95 4.642881e-95 1.000000e+00
[148,] 5.511402e-95 1.102280e-94 1.000000e+00
[149,] 8.170093e-96 1.634019e-95 1.000000e+00
[150,] 4.861590e-96 9.723179e-96 1.000000e+00
[151,] 4.134672e-96 8.269343e-96 1.000000e+00
[152,] 3.357971e-96 6.715941e-96 1.000000e+00
[153,] 1.876021e-62 3.752041e-62 1.000000e+00
[154,] 1.730938e-61 3.461875e-61 1.000000e+00
[155,] 1.667357e-60 3.334714e-60 1.000000e+00
[156,] 9.645912e-61 1.929182e-60 1.000000e+00
[157,] 3.151641e-61 6.303282e-61 1.000000e+00
[158,] 7.461043e-62 1.492209e-61 1.000000e+00
[159,] 1.480234e-61 2.960468e-61 1.000000e+00
[160,] 1.491146e-60 2.982292e-60 1.000000e+00
[161,] 3.370352e-59 6.740704e-59 1.000000e+00
[162,] 1.458755e-57 2.917511e-57 1.000000e+00
[163,] 8.232158e-56 1.646432e-55 1.000000e+00
[164,] 4.721167e-56 9.442335e-56 1.000000e+00
[165,] 1.251629e-56 2.503259e-56 1.000000e+00
[166,] 3.327049e-57 6.654097e-57 1.000000e+00
[167,] 8.682457e-58 1.736491e-57 1.000000e+00
[168,] 2.259557e-58 4.519114e-58 1.000000e+00
[169,] 1.529816e-58 3.059633e-58 1.000000e+00
[170,] 4.032464e-59 8.064928e-59 1.000000e+00
[171,] 1.050563e-59 2.101126e-59 1.000000e+00
[172,] 2.762677e-60 5.525355e-60 1.000000e+00
[173,] 7.333456e-61 1.466691e-60 1.000000e+00
[174,] 1.863433e-61 3.726866e-61 1.000000e+00
[175,] 4.026934e-61 8.053869e-61 1.000000e+00
[176,] 4.183634e-59 8.367267e-59 1.000000e+00
[177,] 3.600209e-56 7.200418e-56 1.000000e+00
[178,] 6.218418e-55 1.243684e-54 1.000000e+00
[179,] 5.056776e-55 1.011355e-54 1.000000e+00
[180,] 1.574265e-55 3.148530e-55 1.000000e+00
[181,] 4.411123e-56 8.822246e-56 1.000000e+00
[182,] 1.264823e-36 2.529646e-36 1.000000e+00
[183,] 6.894685e-33 1.378937e-32 1.000000e+00
[184,] 1.961556e-32 3.923113e-32 1.000000e+00
[185,] 8.060923e-33 1.612185e-32 1.000000e+00
[186,] 1.039689e-32 2.079379e-32 1.000000e+00
[187,] 3.205093e-32 6.410186e-32 1.000000e+00
[188,] 6.335968e-31 1.267194e-30 1.000000e+00
[189,] 6.210115e-26 1.242023e-25 1.000000e+00
[190,] 1.205103e-22 2.410206e-22 1.000000e+00
[191,] 2.869576e-17 5.739152e-17 1.000000e+00
[192,] 1.000000e+00 1.474361e-23 7.371803e-24
[193,] 1.000000e+00 4.503053e-28 2.251527e-28
[194,] 1.000000e+00 6.889273e-28 3.444636e-28
[195,] 1.000000e+00 1.186007e-27 5.930036e-28
[196,] 1.000000e+00 2.069678e-27 1.034839e-27
[197,] 1.000000e+00 2.132726e-27 1.066363e-27
[198,] 1.000000e+00 4.542678e-27 2.271339e-27
[199,] 1.000000e+00 9.846110e-27 4.923055e-27
[200,] 1.000000e+00 4.271587e-31 2.135794e-31
[201,] 1.000000e+00 1.151135e-33 5.755676e-34
[202,] 1.000000e+00 6.207046e-35 3.103523e-35
[203,] 1.000000e+00 6.729252e-35 3.364626e-35
[204,] 1.000000e+00 1.633606e-34 8.168029e-35
[205,] 1.000000e+00 3.296216e-34 1.648108e-34
[206,] 1.000000e+00 8.135102e-34 4.067551e-34
[207,] 1.000000e+00 1.992881e-33 9.964404e-34
[208,] 1.000000e+00 4.882634e-33 2.441317e-33
[209,] 1.000000e+00 1.003128e-32 5.015638e-33
[210,] 1.000000e+00 1.224183e-32 6.120913e-33
[211,] 1.000000e+00 2.495347e-32 1.247673e-32
[212,] 1.000000e+00 5.291688e-32 2.645844e-32
[213,] 1.000000e+00 5.350090e-32 2.675045e-32
[214,] 1.000000e+00 1.949754e-32 9.748771e-33
[215,] 1.000000e+00 4.790391e-33 2.395195e-33
[216,] 1.000000e+00 1.171887e-32 5.859437e-33
[217,] 1.000000e+00 2.413908e-32 1.206954e-32
[218,] 1.000000e+00 3.482572e-32 1.741286e-32
[219,] 1.000000e+00 5.958315e-32 2.979157e-32
[220,] 1.000000e+00 1.333746e-31 6.668730e-32
[221,] 1.000000e+00 3.209895e-31 1.604947e-31
[222,] 1.000000e+00 6.612984e-31 3.306492e-31
[223,] 1.000000e+00 1.417258e-30 7.086290e-31
[224,] 1.000000e+00 2.150580e-30 1.075290e-30
[225,] 1.000000e+00 4.086586e-30 2.043293e-30
[226,] 1.000000e+00 3.191033e-30 1.595516e-30
[227,] 1.000000e+00 7.552324e-30 3.776162e-30
[228,] 1.000000e+00 1.457540e-29 7.287699e-30
[229,] 1.000000e+00 2.059534e-29 1.029767e-29
[230,] 1.000000e+00 3.841734e-29 1.920867e-29
[231,] 1.000000e+00 7.585726e-29 3.792863e-29
[232,] 1.000000e+00 1.775449e-28 8.877244e-29
[233,] 1.000000e+00 3.507533e-28 1.753767e-28
[234,] 1.000000e+00 3.891220e-28 1.945610e-28
[235,] 1.000000e+00 4.916826e-28 2.458413e-28
[236,] 1.000000e+00 1.129364e-27 5.646819e-28
[237,] 1.000000e+00 1.700182e-27 8.500909e-28
[238,] 1.000000e+00 3.636779e-27 1.818389e-27
[239,] 1.000000e+00 5.123546e-27 2.561773e-27
[240,] 1.000000e+00 7.979282e-27 3.989641e-27
[241,] 1.000000e+00 1.804294e-26 9.021469e-27
[242,] 1.000000e+00 1.629284e-26 8.146419e-27
[243,] 1.000000e+00 3.650495e-26 1.825247e-26
[244,] 1.000000e+00 7.988970e-26 3.994485e-26
[245,] 1.000000e+00 1.336092e-25 6.680461e-26
[246,] 1.000000e+00 2.973075e-25 1.486537e-25
[247,] 1.000000e+00 4.205903e-25 2.102952e-25
[248,] 1.000000e+00 9.482376e-25 4.741188e-25
[249,] 1.000000e+00 2.048927e-24 1.024463e-24
[250,] 1.000000e+00 3.463721e-24 1.731860e-24
[251,] 1.000000e+00 6.013124e-24 3.006562e-24
[252,] 1.000000e+00 1.299269e-23 6.496343e-24
[253,] 1.000000e+00 2.814913e-23 1.407456e-23
[254,] 1.000000e+00 6.088515e-23 3.044257e-23
[255,] 1.000000e+00 9.183942e-23 4.591971e-23
[256,] 1.000000e+00 2.005668e-22 1.002834e-22
[257,] 1.000000e+00 4.260634e-22 2.130317e-22
[258,] 1.000000e+00 7.452462e-22 3.726231e-22
[259,] 1.000000e+00 1.538350e-21 7.691751e-22
[260,] 1.000000e+00 3.011075e-21 1.505537e-21
[261,] 1.000000e+00 6.292620e-21 3.146310e-21
[262,] 1.000000e+00 1.161086e-20 5.805430e-21
[263,] 1.000000e+00 1.510851e-20 7.554256e-21
[264,] 1.000000e+00 2.677517e-20 1.338759e-20
[265,] 1.000000e+00 5.560293e-20 2.780147e-20
[266,] 1.000000e+00 1.138944e-19 5.694719e-20
[267,] 1.000000e+00 1.900986e-19 9.504930e-20
[268,] 1.000000e+00 3.332276e-19 1.666138e-19
[269,] 1.000000e+00 6.734376e-19 3.367188e-19
[270,] 1.000000e+00 8.694694e-19 4.347347e-19
[271,] 1.000000e+00 1.745485e-18 8.727424e-19
[272,] 1.000000e+00 3.095562e-18 1.547781e-18
[273,] 1.000000e+00 5.602839e-18 2.801420e-18
[274,] 1.000000e+00 1.108989e-17 5.544947e-18
[275,] 1.000000e+00 2.138777e-17 1.069388e-17
[276,] 1.000000e+00 2.706295e-17 1.353148e-17
[277,] 1.000000e+00 3.360065e-17 1.680033e-17
[278,] 1.000000e+00 6.718357e-17 3.359179e-17
[279,] 1.000000e+00 1.305809e-16 6.529043e-17
[280,] 1.000000e+00 2.539504e-16 1.269752e-16
[281,] 1.000000e+00 4.828629e-16 2.414315e-16
[282,] 1.000000e+00 9.263845e-16 4.631923e-16
[283,] 1.000000e+00 1.726144e-15 8.630721e-16
[284,] 1.000000e+00 2.685963e-15 1.342981e-15
[285,] 1.000000e+00 4.930642e-15 2.465321e-15
[286,] 1.000000e+00 9.597994e-15 4.798997e-15
[287,] 1.000000e+00 1.653813e-14 8.269064e-15
[288,] 1.000000e+00 1.038782e-14 5.193910e-15
[289,] 1.000000e+00 1.575689e-14 7.878445e-15
[290,] 1.000000e+00 2.996851e-14 1.498425e-14
[291,] 1.000000e+00 5.664639e-14 2.832320e-14
[292,] 1.000000e+00 8.322961e-14 4.161480e-14
[293,] 1.000000e+00 1.564089e-13 7.820443e-14
[294,] 1.000000e+00 2.686678e-13 1.343339e-13
[295,] 1.000000e+00 4.988416e-13 2.494208e-13
[296,] 1.000000e+00 7.084620e-13 3.542310e-13
[297,] 1.000000e+00 1.305915e-12 6.529573e-13
[298,] 1.000000e+00 2.391656e-12 1.195828e-12
[299,] 1.000000e+00 3.297296e-12 1.648648e-12
[300,] 1.000000e+00 5.979418e-12 2.989709e-12
[301,] 1.000000e+00 1.097838e-11 5.489189e-12
[302,] 1.000000e+00 1.963611e-11 9.818056e-12
[303,] 1.000000e+00 3.457303e-11 1.728651e-11
[304,] 1.000000e+00 6.155306e-11 3.077653e-11
[305,] 1.000000e+00 1.028677e-10 5.143384e-11
[306,] 1.000000e+00 1.677771e-10 8.388854e-11
[307,] 1.000000e+00 2.739325e-10 1.369663e-10
[308,] 1.000000e+00 4.776508e-10 2.388254e-10
[309,] 1.000000e+00 6.325045e-10 3.162522e-10
[310,] 1.000000e+00 1.068288e-09 5.341439e-10
[311,] 1.000000e+00 1.845227e-09 9.226133e-10
[312,] 1.000000e+00 2.972095e-09 1.486047e-09
[313,] 1.000000e+00 5.052324e-09 2.526162e-09
[314,] 1.000000e+00 8.131178e-09 4.065589e-09
[315,] 1.000000e+00 1.278456e-08 6.392281e-09
[316,] 1.000000e+00 2.000727e-08 1.000364e-08
[317,] 1.000000e+00 3.042861e-08 1.521430e-08
[318,] 1.000000e+00 4.891314e-08 2.445657e-08
[319,] 1.000000e+00 6.953548e-08 3.476774e-08
[320,] 9.999999e-01 1.151794e-07 5.758968e-08
[321,] 9.999999e-01 1.798390e-07 8.991949e-08
[322,] 9.999999e-01 2.934950e-07 1.467475e-07
[323,] 9.999998e-01 4.021364e-07 2.010682e-07
[324,] 9.999997e-01 6.434325e-07 3.217163e-07
[325,] 9.999995e-01 1.034654e-06 5.173270e-07
[326,] 9.999992e-01 1.649441e-06 8.247206e-07
[327,] 9.999987e-01 2.609608e-06 1.304804e-06
[328,] 9.999979e-01 4.112455e-06 2.056227e-06
[329,] 9.999973e-01 5.470855e-06 2.735428e-06
[330,] 9.999957e-01 8.539513e-06 4.269756e-06
[331,] 9.999934e-01 1.322380e-05 6.611898e-06
[332,] 9.999898e-01 2.034374e-05 1.017187e-05
[333,] 9.999855e-01 2.891465e-05 1.445732e-05
[334,] 9.999786e-01 4.270397e-05 2.135199e-05
[335,] 9.999681e-01 6.376998e-05 3.188499e-05
[336,] 9.999579e-01 8.414040e-05 4.207020e-05
[337,] 9.999372e-01 1.255094e-04 6.275472e-05
[338,] 9.999101e-01 1.798633e-04 8.993166e-05
[339,] 9.998714e-01 2.572188e-04 1.286094e-04
[340,] 9.998179e-01 3.641578e-04 1.820789e-04
[341,] 9.997381e-01 5.238451e-04 2.619226e-04
[342,] 9.996572e-01 6.856531e-04 3.428266e-04
[343,] 9.995920e-01 8.160356e-04 4.080178e-04
[344,] 9.996210e-01 7.580385e-04 3.790193e-04
[345,] 9.994529e-01 1.094191e-03 5.470955e-04
[346,] 9.992229e-01 1.554169e-03 7.770846e-04
[347,] 9.989165e-01 2.166969e-03 1.083484e-03
[348,] 9.984738e-01 3.052345e-03 1.526172e-03
[349,] 9.979075e-01 4.185084e-03 2.092542e-03
[350,] 9.971000e-01 5.800033e-03 2.900017e-03
[351,] 9.961853e-01 7.629317e-03 3.814658e-03
[352,] 9.947999e-01 1.040010e-02 5.200052e-03
[353,] 9.930281e-01 1.394386e-02 6.971929e-03
[354,] 9.906759e-01 1.864825e-02 9.324127e-03
[355,] 9.876601e-01 2.467980e-02 1.233990e-02
[356,] 9.839275e-01 3.214494e-02 1.607247e-02
[357,] 9.854980e-01 2.900392e-02 1.450196e-02
[358,] 9.814153e-01 3.716946e-02 1.858473e-02
[359,] 9.791657e-01 4.166856e-02 2.083428e-02
[360,] 9.739030e-01 5.219406e-02 2.609703e-02
[361,] 9.664426e-01 6.711473e-02 3.355737e-02
[362,] 9.594525e-01 8.109507e-02 4.054753e-02
[363,] 9.485656e-01 1.028689e-01 5.143444e-02
[364,] 9.356918e-01 1.286165e-01 6.430824e-02
[365,] 9.198175e-01 1.603650e-01 8.018249e-02
[366,] 9.021974e-01 1.956052e-01 9.780259e-02
[367,] 8.940703e-01 2.118595e-01 1.059297e-01
[368,] 8.726429e-01 2.547142e-01 1.273571e-01
[369,] 8.713776e-01 2.572449e-01 1.286224e-01
[370,] 8.455135e-01 3.089730e-01 1.544865e-01
[371,] 8.264222e-01 3.471556e-01 1.735778e-01
[372,] 7.938341e-01 4.123318e-01 2.061659e-01
[373,] 7.571763e-01 4.856474e-01 2.428237e-01
[374,] 7.254314e-01 5.491373e-01 2.745686e-01
[375,] 6.921841e-01 6.156318e-01 3.078159e-01
[376,] 6.497313e-01 7.005374e-01 3.502687e-01
[377,] 6.392751e-01 7.214497e-01 3.607249e-01
[378,] 6.068760e-01 7.862480e-01 3.931240e-01
[379,] 5.722109e-01 8.555782e-01 4.277891e-01
[380,] 5.251709e-01 9.496583e-01 4.748291e-01
[381,] 4.863944e-01 9.727887e-01 5.136056e-01
[382,] 5.434968e-01 9.130064e-01 4.565032e-01
[383,] 4.967866e-01 9.935732e-01 5.032134e-01
[384,] 4.534750e-01 9.069500e-01 5.465250e-01
[385,] 4.156657e-01 8.313314e-01 5.843343e-01
[386,] 3.829232e-01 7.658465e-01 6.170768e-01
[387,] 3.730577e-01 7.461154e-01 6.269423e-01
[388,] 3.194651e-01 6.389303e-01 6.805349e-01
[389,] 2.725796e-01 5.451593e-01 7.274204e-01
[390,] 3.823551e-01 7.647102e-01 6.176449e-01
[391,] 4.055104e-01 8.110207e-01 5.944896e-01
[392,] 3.673070e-01 7.346139e-01 6.326930e-01
[393,] 3.398847e-01 6.797694e-01 6.601153e-01
[394,] 2.962868e-01 5.925737e-01 7.037132e-01
[395,] 2.403008e-01 4.806015e-01 7.596992e-01
[396,] 2.093850e-01 4.187699e-01 7.906150e-01
[397,] 1.726435e-01 3.452870e-01 8.273565e-01
[398,] 6.159476e-01 7.681047e-01 3.840524e-01
[399,] 6.614322e-01 6.771357e-01 3.385678e-01
[400,] 5.787666e-01 8.424669e-01 4.212334e-01
[401,] 9.868875e-01 2.622509e-02 1.311254e-02
[402,] 9.804284e-01 3.914315e-02 1.957157e-02
[403,] 9.776531e-01 4.469371e-02 2.234686e-02
[404,] 9.838742e-01 3.225166e-02 1.612583e-02
[405,] 9.701988e-01 5.960245e-02 2.980122e-02
[406,] 9.725344e-01 5.493115e-02 2.746558e-02
[407,] 9.775781e-01 4.484381e-02 2.242190e-02
[408,] 9.284908e-01 1.430184e-01 7.150920e-02
> postscript(file="/var/www/html/rcomp/tmp/1v60k1291403555.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/2oyio1291403555.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/3oyio1291403555.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/4oyio1291403555.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/5g7zq1291403555.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
113335.10800 45130.75225 46545.05013 -34420.67941 5424.34423
6 7 8 9 10
12980.44965 4571.76742 29962.50628 -13490.37303 -11544.92494
11 12 13 14 15
24570.17846 -15889.25753 -23924.29567 -26893.30496 -9018.85866
16 17 18 19 20
3779.00899 27515.99415 -6241.46070 -20519.37076 -8116.62598
21 22 23 24 25
-12589.46637 -14471.53417 2408.76670 -10831.80199 -7795.27155
26 27 28 29 30
-19300.84783 -17240.67763 -29692.25331 -23405.34088 -24458.12284
31 32 33 34 35
-13191.30356 -11925.65536 -15004.86879 -12380.50849 29221.02536
36 37 38 39 40
-15168.29959 -15572.33782 -10296.54961 4154.31758 -21539.42178
41 42 43 44 45
-11389.91736 -17211.70982 -16784.66783 -13053.68753 -22660.77010
46 47 48 49 50
-15868.74392 5019.74555 -9439.62113 -11958.64841 -23632.18095
51 52 53 54 55
-14516.32925 -9511.20176 1034.40043 -7384.99477 -7877.54491
56 57 58 59 60
-25598.95212 -20651.73288 -20088.45877 -18912.41934 -21088.97828
61 62 63 64 65
-18498.37566 -28912.29181 -19279.34657 -22767.30551 -27627.12892
66 67 68 69 70
-21779.41785 -24302.19435 -19179.95872 -13531.60941 -19280.35643
71 72 73 74 75
-36085.41103 -32163.61247 -19183.31351 -23835.32948 -16178.15683
76 77 78 79 80
-26016.70509 -17671.31133 -27239.45869 -27707.07779 -17969.23635
81 82 83 84 85
-28891.75382 -22893.12164 -13804.05613 -33851.37690 -20080.61241
86 87 88 89 90
-29217.73237 -23250.56212 -20487.84929 -21468.58773 -19961.07012
91 92 93 94 95
-23104.48714 -25816.29203 -29449.61222 -17612.92941 -35534.63835
96 97 98 99 100
-30910.18592 -17219.82343 -29926.51131 -16490.23932 -28507.42842
101 102 103 104 105
-28927.26054 -26906.54069 -22999.55706 -21352.47342 -30399.08372
106 107 108 109 110
-30778.94492 -34876.96197 -22557.95064 -24952.28897 -31354.39071
111 112 113 114 115
-30673.27702 -37577.05489 -28801.14202 -30785.34286 -32091.85803
116 117 118 119 120
-34511.73258 -28786.13027 -26664.21778 -31973.59307 -34461.31157
121 122 123 124 125
-24335.34401 -32539.10786 -28200.86751 -23530.47538 -34223.25258
126 127 128 129 130
-33133.81793 -35622.47982 -116526.42167 -117782.59050 204552.47730
131 132 133 134 135
323605.22394 271531.40850 33606.58577 18274.57571 8732.66803
136 137 138 139 140
5434.01687 53109.35305 18626.76887 -4544.16520 -36533.29349
141 142 143 144 145
-37262.56333 -115273.95080 -118155.84423 -122091.80543 -118606.80328
146 147 148 149 150
203187.47182 278557.86978 257464.28629 223192.63679 258811.75501
151 152 153 154 155
269382.69623 438.79131 32149.22038 -6453.11745 -39157.57744
156 157 158 159 160
-119924.46185 -125084.29672 202049.99327 -7067.94810 -39715.83725
161 162 163 164 165
-119211.15246 -118595.16054 -120684.95157 266884.96032 293101.83845
166 167 168 169 170
267713.40636 -8077.73832 -37133.48134 -40911.55992 -121036.09624
171 172 173 174 175
200629.16357 262755.06917 267611.67777 242574.06279 -9090.79373
176 177 178 179 180
-38788.89093 -38196.34003 -42846.76193 -39848.03568 -38776.32782
181 182 183 184 185
-41331.41486 -40911.67935 -39830.26732 -39891.42128 -42929.23286
186 187 188 189 190
-124904.29186 198510.20739 261108.80095 47056.27266 -10988.94590
191 192 193 194 195
-32206.29979 -43836.73229 154374.09869 263541.32315 -11609.25673
196 197 198 199 200
-44291.80666 -122576.12028 -128142.72863 -123507.23374 -126361.65486
201 202 203 204 205
223554.32665 253939.44249 -13945.58892 4427.44480 -4797.64726
206 207 208 209 210
-13968.71196 -9540.17068 50715.44407 1920.75171 -297.17549
211 212 213 214 215
-15062.26497 -16673.38236 -14721.91863 46660.10862 2950.66652
216 217 218 219 220
15868.46698 3357.44202 5128.96156 2459.02061 -7467.16835
221 222 223 224 225
52191.34466 -355.77664 -358.57927 -18297.94332 -16160.76077
226 227 228 229 230
68180.04213 1308.21947 -8341.74193 47949.53529 -17791.59277
231 232 233 234 235
6177.57183 651.62846 -7252.99262 -7835.72560 47192.08494
236 237 238 239 240
-19099.35521 -17749.62883 75.52651 -13948.26798 50104.42190
241 242 243 244 245
18065.37654 17387.64145 -5038.01339 -20816.06567 -22115.81690
246 247 248 249 250
-18765.16020 -1267.01618 45550.45403 -15369.41567 -17617.46693
251 252 253 254 255
-19400.81691 -1893.38674 61876.37714 -3579.48902 -18008.92973
256 257 258 259 260
-20187.87962 -2524.41719 44380.90666 -5496.23620 -20875.01508
261 262 263 264 265
-19771.34257 -20814.19431 -8061.80218 -3202.55109 -3537.99352
266 267 268 269 270
39319.77330 -8042.81988 -19553.18127 -21717.80034 -837.50405
271 272 273 274 275
3871.05902 -4173.27937 -21168.80318 44409.08172 12138.12916
276 277 278 279 280
-8034.30096 -25347.39844 -20229.30361 -22976.36334 -5435.30095
281 282 283 284 285
41382.42858 -24514.02999 -26146.59111 -23751.18007 -5700.66425
286 287 288 289 290
-13779.98950 -10510.31834 43819.17241 -15305.52162 -10255.17083
291 292 293 294 295
-25458.54604 111.78307 -7028.57779 -17000.41990 -10879.52977
296 297 298 299 300
-15360.00666 2700.35048 -24050.64828 -8384.68807 38983.21516
301 302 303 304 305
-26792.50610 -8218.78037 38604.10516 -27391.22282 -27409.66719
306 307 308 309 310
-8716.85899 38098.98796 -27582.11523 -9092.85885 37390.56217
311 312 313 314 315
-9342.49803 -1502.78960 -9235.35143 -14005.48023 -28577.40319
316 317 318 319 320
-28726.68074 -29992.71083 -27931.13243 -10366.01471 36531.53903
321 322 323 324 325
-14716.68987 -30901.31175 -28725.17001 -9403.26560 -22356.84771
326 327 328 329 330
-22095.20193 -22395.61037 -20693.27009 -14290.21928 35194.35413
331 332 333 334 335
-27731.91095 -29578.04217 -12126.59802 34691.77317 -9802.73191
336 337 338 339 340
-31105.47079 -13776.48290 -10631.80979 -12880.37293 33931.21906
341 342 343 344 345
-31738.80115 -12056.07443 -12937.21948 -24616.08131 -23213.63104
346 347 348 349 350
-8241.84614 33047.29505 -35550.08400 -31140.79258 -31425.35871
351 352 353 354 355
-31967.14482 -13994.75781 -29742.15094 -31288.82146 34027.61022
356 357 358 359 360
-10733.49685 -20353.61462 -17601.18599 -17550.87755 -14537.81175
361 362 363 364 365
-16831.91553 -19240.14729 -17819.41964 -11234.68222 6617.70766
366 367 368 369 370
-21581.84772 -26977.38552 39927.58782 -17981.80119 -58646.59526
371 372 373 374 375
-11388.33382 -13822.81897 -20258.44233 -10599.93479 -17654.98405
376 377 378 379 380
-9714.81934 -12943.62358 -22570.45093 -22929.42119 -11463.76019
381 382 383 384 385
-11088.50753 -19155.80349 -25608.70333 -7308.63846 -3226.43011
386 387 388 389 390
-10514.98416 -47994.81721 -21574.78127 -14012.95385 -17934.51060
391 392 393 394 395
-11550.59963 -6325.45106 -97914.37325 -28104.92934 -30551.49193
396 397 398 399 400
-17038.43413 -33475.29375 75432.31451 -5111.89402 -6361.50908
401 402 403 404 405
-56484.65217 145879.80998 -8145.41065 71913.35675 -110089.68199
406 407 408 409 410
47527.37604 72877.11387 45254.61363 150445.11606 -26117.47344
411 412 413 414 415
-38978.33012 160079.75443 -38597.23496 -87656.23635 -52573.07115
416 417 418 419 420
10738.23220 62576.78751 199667.35070 28516.23724 -47632.17632
421 422 423 424 425
34948.10466 -4972.85033 -6984.54617 78197.79206 29044.13301
426 427 428 429 430
45755.64161 -23398.74984 -10880.62655 -57808.15262 -56167.22682
431
-36555.74229
> postscript(file="/var/www/html/rcomp/tmp/6g7zq1291403555.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 113335.10800 NA
1 45130.75225 113335.10800
2 46545.05013 45130.75225
3 -34420.67941 46545.05013
4 5424.34423 -34420.67941
5 12980.44965 5424.34423
6 4571.76742 12980.44965
7 29962.50628 4571.76742
8 -13490.37303 29962.50628
9 -11544.92494 -13490.37303
10 24570.17846 -11544.92494
11 -15889.25753 24570.17846
12 -23924.29567 -15889.25753
13 -26893.30496 -23924.29567
14 -9018.85866 -26893.30496
15 3779.00899 -9018.85866
16 27515.99415 3779.00899
17 -6241.46070 27515.99415
18 -20519.37076 -6241.46070
19 -8116.62598 -20519.37076
20 -12589.46637 -8116.62598
21 -14471.53417 -12589.46637
22 2408.76670 -14471.53417
23 -10831.80199 2408.76670
24 -7795.27155 -10831.80199
25 -19300.84783 -7795.27155
26 -17240.67763 -19300.84783
27 -29692.25331 -17240.67763
28 -23405.34088 -29692.25331
29 -24458.12284 -23405.34088
30 -13191.30356 -24458.12284
31 -11925.65536 -13191.30356
32 -15004.86879 -11925.65536
33 -12380.50849 -15004.86879
34 29221.02536 -12380.50849
35 -15168.29959 29221.02536
36 -15572.33782 -15168.29959
37 -10296.54961 -15572.33782
38 4154.31758 -10296.54961
39 -21539.42178 4154.31758
40 -11389.91736 -21539.42178
41 -17211.70982 -11389.91736
42 -16784.66783 -17211.70982
43 -13053.68753 -16784.66783
44 -22660.77010 -13053.68753
45 -15868.74392 -22660.77010
46 5019.74555 -15868.74392
47 -9439.62113 5019.74555
48 -11958.64841 -9439.62113
49 -23632.18095 -11958.64841
50 -14516.32925 -23632.18095
51 -9511.20176 -14516.32925
52 1034.40043 -9511.20176
53 -7384.99477 1034.40043
54 -7877.54491 -7384.99477
55 -25598.95212 -7877.54491
56 -20651.73288 -25598.95212
57 -20088.45877 -20651.73288
58 -18912.41934 -20088.45877
59 -21088.97828 -18912.41934
60 -18498.37566 -21088.97828
61 -28912.29181 -18498.37566
62 -19279.34657 -28912.29181
63 -22767.30551 -19279.34657
64 -27627.12892 -22767.30551
65 -21779.41785 -27627.12892
66 -24302.19435 -21779.41785
67 -19179.95872 -24302.19435
68 -13531.60941 -19179.95872
69 -19280.35643 -13531.60941
70 -36085.41103 -19280.35643
71 -32163.61247 -36085.41103
72 -19183.31351 -32163.61247
73 -23835.32948 -19183.31351
74 -16178.15683 -23835.32948
75 -26016.70509 -16178.15683
76 -17671.31133 -26016.70509
77 -27239.45869 -17671.31133
78 -27707.07779 -27239.45869
79 -17969.23635 -27707.07779
80 -28891.75382 -17969.23635
81 -22893.12164 -28891.75382
82 -13804.05613 -22893.12164
83 -33851.37690 -13804.05613
84 -20080.61241 -33851.37690
85 -29217.73237 -20080.61241
86 -23250.56212 -29217.73237
87 -20487.84929 -23250.56212
88 -21468.58773 -20487.84929
89 -19961.07012 -21468.58773
90 -23104.48714 -19961.07012
91 -25816.29203 -23104.48714
92 -29449.61222 -25816.29203
93 -17612.92941 -29449.61222
94 -35534.63835 -17612.92941
95 -30910.18592 -35534.63835
96 -17219.82343 -30910.18592
97 -29926.51131 -17219.82343
98 -16490.23932 -29926.51131
99 -28507.42842 -16490.23932
100 -28927.26054 -28507.42842
101 -26906.54069 -28927.26054
102 -22999.55706 -26906.54069
103 -21352.47342 -22999.55706
104 -30399.08372 -21352.47342
105 -30778.94492 -30399.08372
106 -34876.96197 -30778.94492
107 -22557.95064 -34876.96197
108 -24952.28897 -22557.95064
109 -31354.39071 -24952.28897
110 -30673.27702 -31354.39071
111 -37577.05489 -30673.27702
112 -28801.14202 -37577.05489
113 -30785.34286 -28801.14202
114 -32091.85803 -30785.34286
115 -34511.73258 -32091.85803
116 -28786.13027 -34511.73258
117 -26664.21778 -28786.13027
118 -31973.59307 -26664.21778
119 -34461.31157 -31973.59307
120 -24335.34401 -34461.31157
121 -32539.10786 -24335.34401
122 -28200.86751 -32539.10786
123 -23530.47538 -28200.86751
124 -34223.25258 -23530.47538
125 -33133.81793 -34223.25258
126 -35622.47982 -33133.81793
127 -116526.42167 -35622.47982
128 -117782.59050 -116526.42167
129 204552.47730 -117782.59050
130 323605.22394 204552.47730
131 271531.40850 323605.22394
132 33606.58577 271531.40850
133 18274.57571 33606.58577
134 8732.66803 18274.57571
135 5434.01687 8732.66803
136 53109.35305 5434.01687
137 18626.76887 53109.35305
138 -4544.16520 18626.76887
139 -36533.29349 -4544.16520
140 -37262.56333 -36533.29349
141 -115273.95080 -37262.56333
142 -118155.84423 -115273.95080
143 -122091.80543 -118155.84423
144 -118606.80328 -122091.80543
145 203187.47182 -118606.80328
146 278557.86978 203187.47182
147 257464.28629 278557.86978
148 223192.63679 257464.28629
149 258811.75501 223192.63679
150 269382.69623 258811.75501
151 438.79131 269382.69623
152 32149.22038 438.79131
153 -6453.11745 32149.22038
154 -39157.57744 -6453.11745
155 -119924.46185 -39157.57744
156 -125084.29672 -119924.46185
157 202049.99327 -125084.29672
158 -7067.94810 202049.99327
159 -39715.83725 -7067.94810
160 -119211.15246 -39715.83725
161 -118595.16054 -119211.15246
162 -120684.95157 -118595.16054
163 266884.96032 -120684.95157
164 293101.83845 266884.96032
165 267713.40636 293101.83845
166 -8077.73832 267713.40636
167 -37133.48134 -8077.73832
168 -40911.55992 -37133.48134
169 -121036.09624 -40911.55992
170 200629.16357 -121036.09624
171 262755.06917 200629.16357
172 267611.67777 262755.06917
173 242574.06279 267611.67777
174 -9090.79373 242574.06279
175 -38788.89093 -9090.79373
176 -38196.34003 -38788.89093
177 -42846.76193 -38196.34003
178 -39848.03568 -42846.76193
179 -38776.32782 -39848.03568
180 -41331.41486 -38776.32782
181 -40911.67935 -41331.41486
182 -39830.26732 -40911.67935
183 -39891.42128 -39830.26732
184 -42929.23286 -39891.42128
185 -124904.29186 -42929.23286
186 198510.20739 -124904.29186
187 261108.80095 198510.20739
188 47056.27266 261108.80095
189 -10988.94590 47056.27266
190 -32206.29979 -10988.94590
191 -43836.73229 -32206.29979
192 154374.09869 -43836.73229
193 263541.32315 154374.09869
194 -11609.25673 263541.32315
195 -44291.80666 -11609.25673
196 -122576.12028 -44291.80666
197 -128142.72863 -122576.12028
198 -123507.23374 -128142.72863
199 -126361.65486 -123507.23374
200 223554.32665 -126361.65486
201 253939.44249 223554.32665
202 -13945.58892 253939.44249
203 4427.44480 -13945.58892
204 -4797.64726 4427.44480
205 -13968.71196 -4797.64726
206 -9540.17068 -13968.71196
207 50715.44407 -9540.17068
208 1920.75171 50715.44407
209 -297.17549 1920.75171
210 -15062.26497 -297.17549
211 -16673.38236 -15062.26497
212 -14721.91863 -16673.38236
213 46660.10862 -14721.91863
214 2950.66652 46660.10862
215 15868.46698 2950.66652
216 3357.44202 15868.46698
217 5128.96156 3357.44202
218 2459.02061 5128.96156
219 -7467.16835 2459.02061
220 52191.34466 -7467.16835
221 -355.77664 52191.34466
222 -358.57927 -355.77664
223 -18297.94332 -358.57927
224 -16160.76077 -18297.94332
225 68180.04213 -16160.76077
226 1308.21947 68180.04213
227 -8341.74193 1308.21947
228 47949.53529 -8341.74193
229 -17791.59277 47949.53529
230 6177.57183 -17791.59277
231 651.62846 6177.57183
232 -7252.99262 651.62846
233 -7835.72560 -7252.99262
234 47192.08494 -7835.72560
235 -19099.35521 47192.08494
236 -17749.62883 -19099.35521
237 75.52651 -17749.62883
238 -13948.26798 75.52651
239 50104.42190 -13948.26798
240 18065.37654 50104.42190
241 17387.64145 18065.37654
242 -5038.01339 17387.64145
243 -20816.06567 -5038.01339
244 -22115.81690 -20816.06567
245 -18765.16020 -22115.81690
246 -1267.01618 -18765.16020
247 45550.45403 -1267.01618
248 -15369.41567 45550.45403
249 -17617.46693 -15369.41567
250 -19400.81691 -17617.46693
251 -1893.38674 -19400.81691
252 61876.37714 -1893.38674
253 -3579.48902 61876.37714
254 -18008.92973 -3579.48902
255 -20187.87962 -18008.92973
256 -2524.41719 -20187.87962
257 44380.90666 -2524.41719
258 -5496.23620 44380.90666
259 -20875.01508 -5496.23620
260 -19771.34257 -20875.01508
261 -20814.19431 -19771.34257
262 -8061.80218 -20814.19431
263 -3202.55109 -8061.80218
264 -3537.99352 -3202.55109
265 39319.77330 -3537.99352
266 -8042.81988 39319.77330
267 -19553.18127 -8042.81988
268 -21717.80034 -19553.18127
269 -837.50405 -21717.80034
270 3871.05902 -837.50405
271 -4173.27937 3871.05902
272 -21168.80318 -4173.27937
273 44409.08172 -21168.80318
274 12138.12916 44409.08172
275 -8034.30096 12138.12916
276 -25347.39844 -8034.30096
277 -20229.30361 -25347.39844
278 -22976.36334 -20229.30361
279 -5435.30095 -22976.36334
280 41382.42858 -5435.30095
281 -24514.02999 41382.42858
282 -26146.59111 -24514.02999
283 -23751.18007 -26146.59111
284 -5700.66425 -23751.18007
285 -13779.98950 -5700.66425
286 -10510.31834 -13779.98950
287 43819.17241 -10510.31834
288 -15305.52162 43819.17241
289 -10255.17083 -15305.52162
290 -25458.54604 -10255.17083
291 111.78307 -25458.54604
292 -7028.57779 111.78307
293 -17000.41990 -7028.57779
294 -10879.52977 -17000.41990
295 -15360.00666 -10879.52977
296 2700.35048 -15360.00666
297 -24050.64828 2700.35048
298 -8384.68807 -24050.64828
299 38983.21516 -8384.68807
300 -26792.50610 38983.21516
301 -8218.78037 -26792.50610
302 38604.10516 -8218.78037
303 -27391.22282 38604.10516
304 -27409.66719 -27391.22282
305 -8716.85899 -27409.66719
306 38098.98796 -8716.85899
307 -27582.11523 38098.98796
308 -9092.85885 -27582.11523
309 37390.56217 -9092.85885
310 -9342.49803 37390.56217
311 -1502.78960 -9342.49803
312 -9235.35143 -1502.78960
313 -14005.48023 -9235.35143
314 -28577.40319 -14005.48023
315 -28726.68074 -28577.40319
316 -29992.71083 -28726.68074
317 -27931.13243 -29992.71083
318 -10366.01471 -27931.13243
319 36531.53903 -10366.01471
320 -14716.68987 36531.53903
321 -30901.31175 -14716.68987
322 -28725.17001 -30901.31175
323 -9403.26560 -28725.17001
324 -22356.84771 -9403.26560
325 -22095.20193 -22356.84771
326 -22395.61037 -22095.20193
327 -20693.27009 -22395.61037
328 -14290.21928 -20693.27009
329 35194.35413 -14290.21928
330 -27731.91095 35194.35413
331 -29578.04217 -27731.91095
332 -12126.59802 -29578.04217
333 34691.77317 -12126.59802
334 -9802.73191 34691.77317
335 -31105.47079 -9802.73191
336 -13776.48290 -31105.47079
337 -10631.80979 -13776.48290
338 -12880.37293 -10631.80979
339 33931.21906 -12880.37293
340 -31738.80115 33931.21906
341 -12056.07443 -31738.80115
342 -12937.21948 -12056.07443
343 -24616.08131 -12937.21948
344 -23213.63104 -24616.08131
345 -8241.84614 -23213.63104
346 33047.29505 -8241.84614
347 -35550.08400 33047.29505
348 -31140.79258 -35550.08400
349 -31425.35871 -31140.79258
350 -31967.14482 -31425.35871
351 -13994.75781 -31967.14482
352 -29742.15094 -13994.75781
353 -31288.82146 -29742.15094
354 34027.61022 -31288.82146
355 -10733.49685 34027.61022
356 -20353.61462 -10733.49685
357 -17601.18599 -20353.61462
358 -17550.87755 -17601.18599
359 -14537.81175 -17550.87755
360 -16831.91553 -14537.81175
361 -19240.14729 -16831.91553
362 -17819.41964 -19240.14729
363 -11234.68222 -17819.41964
364 6617.70766 -11234.68222
365 -21581.84772 6617.70766
366 -26977.38552 -21581.84772
367 39927.58782 -26977.38552
368 -17981.80119 39927.58782
369 -58646.59526 -17981.80119
370 -11388.33382 -58646.59526
371 -13822.81897 -11388.33382
372 -20258.44233 -13822.81897
373 -10599.93479 -20258.44233
374 -17654.98405 -10599.93479
375 -9714.81934 -17654.98405
376 -12943.62358 -9714.81934
377 -22570.45093 -12943.62358
378 -22929.42119 -22570.45093
379 -11463.76019 -22929.42119
380 -11088.50753 -11463.76019
381 -19155.80349 -11088.50753
382 -25608.70333 -19155.80349
383 -7308.63846 -25608.70333
384 -3226.43011 -7308.63846
385 -10514.98416 -3226.43011
386 -47994.81721 -10514.98416
387 -21574.78127 -47994.81721
388 -14012.95385 -21574.78127
389 -17934.51060 -14012.95385
390 -11550.59963 -17934.51060
391 -6325.45106 -11550.59963
392 -97914.37325 -6325.45106
393 -28104.92934 -97914.37325
394 -30551.49193 -28104.92934
395 -17038.43413 -30551.49193
396 -33475.29375 -17038.43413
397 75432.31451 -33475.29375
398 -5111.89402 75432.31451
399 -6361.50908 -5111.89402
400 -56484.65217 -6361.50908
401 145879.80998 -56484.65217
402 -8145.41065 145879.80998
403 71913.35675 -8145.41065
404 -110089.68199 71913.35675
405 47527.37604 -110089.68199
406 72877.11387 47527.37604
407 45254.61363 72877.11387
408 150445.11606 45254.61363
409 -26117.47344 150445.11606
410 -38978.33012 -26117.47344
411 160079.75443 -38978.33012
412 -38597.23496 160079.75443
413 -87656.23635 -38597.23496
414 -52573.07115 -87656.23635
415 10738.23220 -52573.07115
416 62576.78751 10738.23220
417 199667.35070 62576.78751
418 28516.23724 199667.35070
419 -47632.17632 28516.23724
420 34948.10466 -47632.17632
421 -4972.85033 34948.10466
422 -6984.54617 -4972.85033
423 78197.79206 -6984.54617
424 29044.13301 78197.79206
425 45755.64161 29044.13301
426 -23398.74984 45755.64161
427 -10880.62655 -23398.74984
428 -57808.15262 -10880.62655
429 -56167.22682 -57808.15262
430 -36555.74229 -56167.22682
431 NA -36555.74229
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 45130.75225 113335.10800
[2,] 46545.05013 45130.75225
[3,] -34420.67941 46545.05013
[4,] 5424.34423 -34420.67941
[5,] 12980.44965 5424.34423
[6,] 4571.76742 12980.44965
[7,] 29962.50628 4571.76742
[8,] -13490.37303 29962.50628
[9,] -11544.92494 -13490.37303
[10,] 24570.17846 -11544.92494
[11,] -15889.25753 24570.17846
[12,] -23924.29567 -15889.25753
[13,] -26893.30496 -23924.29567
[14,] -9018.85866 -26893.30496
[15,] 3779.00899 -9018.85866
[16,] 27515.99415 3779.00899
[17,] -6241.46070 27515.99415
[18,] -20519.37076 -6241.46070
[19,] -8116.62598 -20519.37076
[20,] -12589.46637 -8116.62598
[21,] -14471.53417 -12589.46637
[22,] 2408.76670 -14471.53417
[23,] -10831.80199 2408.76670
[24,] -7795.27155 -10831.80199
[25,] -19300.84783 -7795.27155
[26,] -17240.67763 -19300.84783
[27,] -29692.25331 -17240.67763
[28,] -23405.34088 -29692.25331
[29,] -24458.12284 -23405.34088
[30,] -13191.30356 -24458.12284
[31,] -11925.65536 -13191.30356
[32,] -15004.86879 -11925.65536
[33,] -12380.50849 -15004.86879
[34,] 29221.02536 -12380.50849
[35,] -15168.29959 29221.02536
[36,] -15572.33782 -15168.29959
[37,] -10296.54961 -15572.33782
[38,] 4154.31758 -10296.54961
[39,] -21539.42178 4154.31758
[40,] -11389.91736 -21539.42178
[41,] -17211.70982 -11389.91736
[42,] -16784.66783 -17211.70982
[43,] -13053.68753 -16784.66783
[44,] -22660.77010 -13053.68753
[45,] -15868.74392 -22660.77010
[46,] 5019.74555 -15868.74392
[47,] -9439.62113 5019.74555
[48,] -11958.64841 -9439.62113
[49,] -23632.18095 -11958.64841
[50,] -14516.32925 -23632.18095
[51,] -9511.20176 -14516.32925
[52,] 1034.40043 -9511.20176
[53,] -7384.99477 1034.40043
[54,] -7877.54491 -7384.99477
[55,] -25598.95212 -7877.54491
[56,] -20651.73288 -25598.95212
[57,] -20088.45877 -20651.73288
[58,] -18912.41934 -20088.45877
[59,] -21088.97828 -18912.41934
[60,] -18498.37566 -21088.97828
[61,] -28912.29181 -18498.37566
[62,] -19279.34657 -28912.29181
[63,] -22767.30551 -19279.34657
[64,] -27627.12892 -22767.30551
[65,] -21779.41785 -27627.12892
[66,] -24302.19435 -21779.41785
[67,] -19179.95872 -24302.19435
[68,] -13531.60941 -19179.95872
[69,] -19280.35643 -13531.60941
[70,] -36085.41103 -19280.35643
[71,] -32163.61247 -36085.41103
[72,] -19183.31351 -32163.61247
[73,] -23835.32948 -19183.31351
[74,] -16178.15683 -23835.32948
[75,] -26016.70509 -16178.15683
[76,] -17671.31133 -26016.70509
[77,] -27239.45869 -17671.31133
[78,] -27707.07779 -27239.45869
[79,] -17969.23635 -27707.07779
[80,] -28891.75382 -17969.23635
[81,] -22893.12164 -28891.75382
[82,] -13804.05613 -22893.12164
[83,] -33851.37690 -13804.05613
[84,] -20080.61241 -33851.37690
[85,] -29217.73237 -20080.61241
[86,] -23250.56212 -29217.73237
[87,] -20487.84929 -23250.56212
[88,] -21468.58773 -20487.84929
[89,] -19961.07012 -21468.58773
[90,] -23104.48714 -19961.07012
[91,] -25816.29203 -23104.48714
[92,] -29449.61222 -25816.29203
[93,] -17612.92941 -29449.61222
[94,] -35534.63835 -17612.92941
[95,] -30910.18592 -35534.63835
[96,] -17219.82343 -30910.18592
[97,] -29926.51131 -17219.82343
[98,] -16490.23932 -29926.51131
[99,] -28507.42842 -16490.23932
[100,] -28927.26054 -28507.42842
[101,] -26906.54069 -28927.26054
[102,] -22999.55706 -26906.54069
[103,] -21352.47342 -22999.55706
[104,] -30399.08372 -21352.47342
[105,] -30778.94492 -30399.08372
[106,] -34876.96197 -30778.94492
[107,] -22557.95064 -34876.96197
[108,] -24952.28897 -22557.95064
[109,] -31354.39071 -24952.28897
[110,] -30673.27702 -31354.39071
[111,] -37577.05489 -30673.27702
[112,] -28801.14202 -37577.05489
[113,] -30785.34286 -28801.14202
[114,] -32091.85803 -30785.34286
[115,] -34511.73258 -32091.85803
[116,] -28786.13027 -34511.73258
[117,] -26664.21778 -28786.13027
[118,] -31973.59307 -26664.21778
[119,] -34461.31157 -31973.59307
[120,] -24335.34401 -34461.31157
[121,] -32539.10786 -24335.34401
[122,] -28200.86751 -32539.10786
[123,] -23530.47538 -28200.86751
[124,] -34223.25258 -23530.47538
[125,] -33133.81793 -34223.25258
[126,] -35622.47982 -33133.81793
[127,] -116526.42167 -35622.47982
[128,] -117782.59050 -116526.42167
[129,] 204552.47730 -117782.59050
[130,] 323605.22394 204552.47730
[131,] 271531.40850 323605.22394
[132,] 33606.58577 271531.40850
[133,] 18274.57571 33606.58577
[134,] 8732.66803 18274.57571
[135,] 5434.01687 8732.66803
[136,] 53109.35305 5434.01687
[137,] 18626.76887 53109.35305
[138,] -4544.16520 18626.76887
[139,] -36533.29349 -4544.16520
[140,] -37262.56333 -36533.29349
[141,] -115273.95080 -37262.56333
[142,] -118155.84423 -115273.95080
[143,] -122091.80543 -118155.84423
[144,] -118606.80328 -122091.80543
[145,] 203187.47182 -118606.80328
[146,] 278557.86978 203187.47182
[147,] 257464.28629 278557.86978
[148,] 223192.63679 257464.28629
[149,] 258811.75501 223192.63679
[150,] 269382.69623 258811.75501
[151,] 438.79131 269382.69623
[152,] 32149.22038 438.79131
[153,] -6453.11745 32149.22038
[154,] -39157.57744 -6453.11745
[155,] -119924.46185 -39157.57744
[156,] -125084.29672 -119924.46185
[157,] 202049.99327 -125084.29672
[158,] -7067.94810 202049.99327
[159,] -39715.83725 -7067.94810
[160,] -119211.15246 -39715.83725
[161,] -118595.16054 -119211.15246
[162,] -120684.95157 -118595.16054
[163,] 266884.96032 -120684.95157
[164,] 293101.83845 266884.96032
[165,] 267713.40636 293101.83845
[166,] -8077.73832 267713.40636
[167,] -37133.48134 -8077.73832
[168,] -40911.55992 -37133.48134
[169,] -121036.09624 -40911.55992
[170,] 200629.16357 -121036.09624
[171,] 262755.06917 200629.16357
[172,] 267611.67777 262755.06917
[173,] 242574.06279 267611.67777
[174,] -9090.79373 242574.06279
[175,] -38788.89093 -9090.79373
[176,] -38196.34003 -38788.89093
[177,] -42846.76193 -38196.34003
[178,] -39848.03568 -42846.76193
[179,] -38776.32782 -39848.03568
[180,] -41331.41486 -38776.32782
[181,] -40911.67935 -41331.41486
[182,] -39830.26732 -40911.67935
[183,] -39891.42128 -39830.26732
[184,] -42929.23286 -39891.42128
[185,] -124904.29186 -42929.23286
[186,] 198510.20739 -124904.29186
[187,] 261108.80095 198510.20739
[188,] 47056.27266 261108.80095
[189,] -10988.94590 47056.27266
[190,] -32206.29979 -10988.94590
[191,] -43836.73229 -32206.29979
[192,] 154374.09869 -43836.73229
[193,] 263541.32315 154374.09869
[194,] -11609.25673 263541.32315
[195,] -44291.80666 -11609.25673
[196,] -122576.12028 -44291.80666
[197,] -128142.72863 -122576.12028
[198,] -123507.23374 -128142.72863
[199,] -126361.65486 -123507.23374
[200,] 223554.32665 -126361.65486
[201,] 253939.44249 223554.32665
[202,] -13945.58892 253939.44249
[203,] 4427.44480 -13945.58892
[204,] -4797.64726 4427.44480
[205,] -13968.71196 -4797.64726
[206,] -9540.17068 -13968.71196
[207,] 50715.44407 -9540.17068
[208,] 1920.75171 50715.44407
[209,] -297.17549 1920.75171
[210,] -15062.26497 -297.17549
[211,] -16673.38236 -15062.26497
[212,] -14721.91863 -16673.38236
[213,] 46660.10862 -14721.91863
[214,] 2950.66652 46660.10862
[215,] 15868.46698 2950.66652
[216,] 3357.44202 15868.46698
[217,] 5128.96156 3357.44202
[218,] 2459.02061 5128.96156
[219,] -7467.16835 2459.02061
[220,] 52191.34466 -7467.16835
[221,] -355.77664 52191.34466
[222,] -358.57927 -355.77664
[223,] -18297.94332 -358.57927
[224,] -16160.76077 -18297.94332
[225,] 68180.04213 -16160.76077
[226,] 1308.21947 68180.04213
[227,] -8341.74193 1308.21947
[228,] 47949.53529 -8341.74193
[229,] -17791.59277 47949.53529
[230,] 6177.57183 -17791.59277
[231,] 651.62846 6177.57183
[232,] -7252.99262 651.62846
[233,] -7835.72560 -7252.99262
[234,] 47192.08494 -7835.72560
[235,] -19099.35521 47192.08494
[236,] -17749.62883 -19099.35521
[237,] 75.52651 -17749.62883
[238,] -13948.26798 75.52651
[239,] 50104.42190 -13948.26798
[240,] 18065.37654 50104.42190
[241,] 17387.64145 18065.37654
[242,] -5038.01339 17387.64145
[243,] -20816.06567 -5038.01339
[244,] -22115.81690 -20816.06567
[245,] -18765.16020 -22115.81690
[246,] -1267.01618 -18765.16020
[247,] 45550.45403 -1267.01618
[248,] -15369.41567 45550.45403
[249,] -17617.46693 -15369.41567
[250,] -19400.81691 -17617.46693
[251,] -1893.38674 -19400.81691
[252,] 61876.37714 -1893.38674
[253,] -3579.48902 61876.37714
[254,] -18008.92973 -3579.48902
[255,] -20187.87962 -18008.92973
[256,] -2524.41719 -20187.87962
[257,] 44380.90666 -2524.41719
[258,] -5496.23620 44380.90666
[259,] -20875.01508 -5496.23620
[260,] -19771.34257 -20875.01508
[261,] -20814.19431 -19771.34257
[262,] -8061.80218 -20814.19431
[263,] -3202.55109 -8061.80218
[264,] -3537.99352 -3202.55109
[265,] 39319.77330 -3537.99352
[266,] -8042.81988 39319.77330
[267,] -19553.18127 -8042.81988
[268,] -21717.80034 -19553.18127
[269,] -837.50405 -21717.80034
[270,] 3871.05902 -837.50405
[271,] -4173.27937 3871.05902
[272,] -21168.80318 -4173.27937
[273,] 44409.08172 -21168.80318
[274,] 12138.12916 44409.08172
[275,] -8034.30096 12138.12916
[276,] -25347.39844 -8034.30096
[277,] -20229.30361 -25347.39844
[278,] -22976.36334 -20229.30361
[279,] -5435.30095 -22976.36334
[280,] 41382.42858 -5435.30095
[281,] -24514.02999 41382.42858
[282,] -26146.59111 -24514.02999
[283,] -23751.18007 -26146.59111
[284,] -5700.66425 -23751.18007
[285,] -13779.98950 -5700.66425
[286,] -10510.31834 -13779.98950
[287,] 43819.17241 -10510.31834
[288,] -15305.52162 43819.17241
[289,] -10255.17083 -15305.52162
[290,] -25458.54604 -10255.17083
[291,] 111.78307 -25458.54604
[292,] -7028.57779 111.78307
[293,] -17000.41990 -7028.57779
[294,] -10879.52977 -17000.41990
[295,] -15360.00666 -10879.52977
[296,] 2700.35048 -15360.00666
[297,] -24050.64828 2700.35048
[298,] -8384.68807 -24050.64828
[299,] 38983.21516 -8384.68807
[300,] -26792.50610 38983.21516
[301,] -8218.78037 -26792.50610
[302,] 38604.10516 -8218.78037
[303,] -27391.22282 38604.10516
[304,] -27409.66719 -27391.22282
[305,] -8716.85899 -27409.66719
[306,] 38098.98796 -8716.85899
[307,] -27582.11523 38098.98796
[308,] -9092.85885 -27582.11523
[309,] 37390.56217 -9092.85885
[310,] -9342.49803 37390.56217
[311,] -1502.78960 -9342.49803
[312,] -9235.35143 -1502.78960
[313,] -14005.48023 -9235.35143
[314,] -28577.40319 -14005.48023
[315,] -28726.68074 -28577.40319
[316,] -29992.71083 -28726.68074
[317,] -27931.13243 -29992.71083
[318,] -10366.01471 -27931.13243
[319,] 36531.53903 -10366.01471
[320,] -14716.68987 36531.53903
[321,] -30901.31175 -14716.68987
[322,] -28725.17001 -30901.31175
[323,] -9403.26560 -28725.17001
[324,] -22356.84771 -9403.26560
[325,] -22095.20193 -22356.84771
[326,] -22395.61037 -22095.20193
[327,] -20693.27009 -22395.61037
[328,] -14290.21928 -20693.27009
[329,] 35194.35413 -14290.21928
[330,] -27731.91095 35194.35413
[331,] -29578.04217 -27731.91095
[332,] -12126.59802 -29578.04217
[333,] 34691.77317 -12126.59802
[334,] -9802.73191 34691.77317
[335,] -31105.47079 -9802.73191
[336,] -13776.48290 -31105.47079
[337,] -10631.80979 -13776.48290
[338,] -12880.37293 -10631.80979
[339,] 33931.21906 -12880.37293
[340,] -31738.80115 33931.21906
[341,] -12056.07443 -31738.80115
[342,] -12937.21948 -12056.07443
[343,] -24616.08131 -12937.21948
[344,] -23213.63104 -24616.08131
[345,] -8241.84614 -23213.63104
[346,] 33047.29505 -8241.84614
[347,] -35550.08400 33047.29505
[348,] -31140.79258 -35550.08400
[349,] -31425.35871 -31140.79258
[350,] -31967.14482 -31425.35871
[351,] -13994.75781 -31967.14482
[352,] -29742.15094 -13994.75781
[353,] -31288.82146 -29742.15094
[354,] 34027.61022 -31288.82146
[355,] -10733.49685 34027.61022
[356,] -20353.61462 -10733.49685
[357,] -17601.18599 -20353.61462
[358,] -17550.87755 -17601.18599
[359,] -14537.81175 -17550.87755
[360,] -16831.91553 -14537.81175
[361,] -19240.14729 -16831.91553
[362,] -17819.41964 -19240.14729
[363,] -11234.68222 -17819.41964
[364,] 6617.70766 -11234.68222
[365,] -21581.84772 6617.70766
[366,] -26977.38552 -21581.84772
[367,] 39927.58782 -26977.38552
[368,] -17981.80119 39927.58782
[369,] -58646.59526 -17981.80119
[370,] -11388.33382 -58646.59526
[371,] -13822.81897 -11388.33382
[372,] -20258.44233 -13822.81897
[373,] -10599.93479 -20258.44233
[374,] -17654.98405 -10599.93479
[375,] -9714.81934 -17654.98405
[376,] -12943.62358 -9714.81934
[377,] -22570.45093 -12943.62358
[378,] -22929.42119 -22570.45093
[379,] -11463.76019 -22929.42119
[380,] -11088.50753 -11463.76019
[381,] -19155.80349 -11088.50753
[382,] -25608.70333 -19155.80349
[383,] -7308.63846 -25608.70333
[384,] -3226.43011 -7308.63846
[385,] -10514.98416 -3226.43011
[386,] -47994.81721 -10514.98416
[387,] -21574.78127 -47994.81721
[388,] -14012.95385 -21574.78127
[389,] -17934.51060 -14012.95385
[390,] -11550.59963 -17934.51060
[391,] -6325.45106 -11550.59963
[392,] -97914.37325 -6325.45106
[393,] -28104.92934 -97914.37325
[394,] -30551.49193 -28104.92934
[395,] -17038.43413 -30551.49193
[396,] -33475.29375 -17038.43413
[397,] 75432.31451 -33475.29375
[398,] -5111.89402 75432.31451
[399,] -6361.50908 -5111.89402
[400,] -56484.65217 -6361.50908
[401,] 145879.80998 -56484.65217
[402,] -8145.41065 145879.80998
[403,] 71913.35675 -8145.41065
[404,] -110089.68199 71913.35675
[405,] 47527.37604 -110089.68199
[406,] 72877.11387 47527.37604
[407,] 45254.61363 72877.11387
[408,] 150445.11606 45254.61363
[409,] -26117.47344 150445.11606
[410,] -38978.33012 -26117.47344
[411,] 160079.75443 -38978.33012
[412,] -38597.23496 160079.75443
[413,] -87656.23635 -38597.23496
[414,] -52573.07115 -87656.23635
[415,] 10738.23220 -52573.07115
[416,] 62576.78751 10738.23220
[417,] 199667.35070 62576.78751
[418,] 28516.23724 199667.35070
[419,] -47632.17632 28516.23724
[420,] 34948.10466 -47632.17632
[421,] -4972.85033 34948.10466
[422,] -6984.54617 -4972.85033
[423,] 78197.79206 -6984.54617
[424,] 29044.13301 78197.79206
[425,] 45755.64161 29044.13301
[426,] -23398.74984 45755.64161
[427,] -10880.62655 -23398.74984
[428,] -57808.15262 -10880.62655
[429,] -56167.22682 -57808.15262
[430,] -36555.74229 -56167.22682
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 45130.75225 113335.10800
2 46545.05013 45130.75225
3 -34420.67941 46545.05013
4 5424.34423 -34420.67941
5 12980.44965 5424.34423
6 4571.76742 12980.44965
7 29962.50628 4571.76742
8 -13490.37303 29962.50628
9 -11544.92494 -13490.37303
10 24570.17846 -11544.92494
11 -15889.25753 24570.17846
12 -23924.29567 -15889.25753
13 -26893.30496 -23924.29567
14 -9018.85866 -26893.30496
15 3779.00899 -9018.85866
16 27515.99415 3779.00899
17 -6241.46070 27515.99415
18 -20519.37076 -6241.46070
19 -8116.62598 -20519.37076
20 -12589.46637 -8116.62598
21 -14471.53417 -12589.46637
22 2408.76670 -14471.53417
23 -10831.80199 2408.76670
24 -7795.27155 -10831.80199
25 -19300.84783 -7795.27155
26 -17240.67763 -19300.84783
27 -29692.25331 -17240.67763
28 -23405.34088 -29692.25331
29 -24458.12284 -23405.34088
30 -13191.30356 -24458.12284
31 -11925.65536 -13191.30356
32 -15004.86879 -11925.65536
33 -12380.50849 -15004.86879
34 29221.02536 -12380.50849
35 -15168.29959 29221.02536
36 -15572.33782 -15168.29959
37 -10296.54961 -15572.33782
38 4154.31758 -10296.54961
39 -21539.42178 4154.31758
40 -11389.91736 -21539.42178
41 -17211.70982 -11389.91736
42 -16784.66783 -17211.70982
43 -13053.68753 -16784.66783
44 -22660.77010 -13053.68753
45 -15868.74392 -22660.77010
46 5019.74555 -15868.74392
47 -9439.62113 5019.74555
48 -11958.64841 -9439.62113
49 -23632.18095 -11958.64841
50 -14516.32925 -23632.18095
51 -9511.20176 -14516.32925
52 1034.40043 -9511.20176
53 -7384.99477 1034.40043
54 -7877.54491 -7384.99477
55 -25598.95212 -7877.54491
56 -20651.73288 -25598.95212
57 -20088.45877 -20651.73288
58 -18912.41934 -20088.45877
59 -21088.97828 -18912.41934
60 -18498.37566 -21088.97828
61 -28912.29181 -18498.37566
62 -19279.34657 -28912.29181
63 -22767.30551 -19279.34657
64 -27627.12892 -22767.30551
65 -21779.41785 -27627.12892
66 -24302.19435 -21779.41785
67 -19179.95872 -24302.19435
68 -13531.60941 -19179.95872
69 -19280.35643 -13531.60941
70 -36085.41103 -19280.35643
71 -32163.61247 -36085.41103
72 -19183.31351 -32163.61247
73 -23835.32948 -19183.31351
74 -16178.15683 -23835.32948
75 -26016.70509 -16178.15683
76 -17671.31133 -26016.70509
77 -27239.45869 -17671.31133
78 -27707.07779 -27239.45869
79 -17969.23635 -27707.07779
80 -28891.75382 -17969.23635
81 -22893.12164 -28891.75382
82 -13804.05613 -22893.12164
83 -33851.37690 -13804.05613
84 -20080.61241 -33851.37690
85 -29217.73237 -20080.61241
86 -23250.56212 -29217.73237
87 -20487.84929 -23250.56212
88 -21468.58773 -20487.84929
89 -19961.07012 -21468.58773
90 -23104.48714 -19961.07012
91 -25816.29203 -23104.48714
92 -29449.61222 -25816.29203
93 -17612.92941 -29449.61222
94 -35534.63835 -17612.92941
95 -30910.18592 -35534.63835
96 -17219.82343 -30910.18592
97 -29926.51131 -17219.82343
98 -16490.23932 -29926.51131
99 -28507.42842 -16490.23932
100 -28927.26054 -28507.42842
101 -26906.54069 -28927.26054
102 -22999.55706 -26906.54069
103 -21352.47342 -22999.55706
104 -30399.08372 -21352.47342
105 -30778.94492 -30399.08372
106 -34876.96197 -30778.94492
107 -22557.95064 -34876.96197
108 -24952.28897 -22557.95064
109 -31354.39071 -24952.28897
110 -30673.27702 -31354.39071
111 -37577.05489 -30673.27702
112 -28801.14202 -37577.05489
113 -30785.34286 -28801.14202
114 -32091.85803 -30785.34286
115 -34511.73258 -32091.85803
116 -28786.13027 -34511.73258
117 -26664.21778 -28786.13027
118 -31973.59307 -26664.21778
119 -34461.31157 -31973.59307
120 -24335.34401 -34461.31157
121 -32539.10786 -24335.34401
122 -28200.86751 -32539.10786
123 -23530.47538 -28200.86751
124 -34223.25258 -23530.47538
125 -33133.81793 -34223.25258
126 -35622.47982 -33133.81793
127 -116526.42167 -35622.47982
128 -117782.59050 -116526.42167
129 204552.47730 -117782.59050
130 323605.22394 204552.47730
131 271531.40850 323605.22394
132 33606.58577 271531.40850
133 18274.57571 33606.58577
134 8732.66803 18274.57571
135 5434.01687 8732.66803
136 53109.35305 5434.01687
137 18626.76887 53109.35305
138 -4544.16520 18626.76887
139 -36533.29349 -4544.16520
140 -37262.56333 -36533.29349
141 -115273.95080 -37262.56333
142 -118155.84423 -115273.95080
143 -122091.80543 -118155.84423
144 -118606.80328 -122091.80543
145 203187.47182 -118606.80328
146 278557.86978 203187.47182
147 257464.28629 278557.86978
148 223192.63679 257464.28629
149 258811.75501 223192.63679
150 269382.69623 258811.75501
151 438.79131 269382.69623
152 32149.22038 438.79131
153 -6453.11745 32149.22038
154 -39157.57744 -6453.11745
155 -119924.46185 -39157.57744
156 -125084.29672 -119924.46185
157 202049.99327 -125084.29672
158 -7067.94810 202049.99327
159 -39715.83725 -7067.94810
160 -119211.15246 -39715.83725
161 -118595.16054 -119211.15246
162 -120684.95157 -118595.16054
163 266884.96032 -120684.95157
164 293101.83845 266884.96032
165 267713.40636 293101.83845
166 -8077.73832 267713.40636
167 -37133.48134 -8077.73832
168 -40911.55992 -37133.48134
169 -121036.09624 -40911.55992
170 200629.16357 -121036.09624
171 262755.06917 200629.16357
172 267611.67777 262755.06917
173 242574.06279 267611.67777
174 -9090.79373 242574.06279
175 -38788.89093 -9090.79373
176 -38196.34003 -38788.89093
177 -42846.76193 -38196.34003
178 -39848.03568 -42846.76193
179 -38776.32782 -39848.03568
180 -41331.41486 -38776.32782
181 -40911.67935 -41331.41486
182 -39830.26732 -40911.67935
183 -39891.42128 -39830.26732
184 -42929.23286 -39891.42128
185 -124904.29186 -42929.23286
186 198510.20739 -124904.29186
187 261108.80095 198510.20739
188 47056.27266 261108.80095
189 -10988.94590 47056.27266
190 -32206.29979 -10988.94590
191 -43836.73229 -32206.29979
192 154374.09869 -43836.73229
193 263541.32315 154374.09869
194 -11609.25673 263541.32315
195 -44291.80666 -11609.25673
196 -122576.12028 -44291.80666
197 -128142.72863 -122576.12028
198 -123507.23374 -128142.72863
199 -126361.65486 -123507.23374
200 223554.32665 -126361.65486
201 253939.44249 223554.32665
202 -13945.58892 253939.44249
203 4427.44480 -13945.58892
204 -4797.64726 4427.44480
205 -13968.71196 -4797.64726
206 -9540.17068 -13968.71196
207 50715.44407 -9540.17068
208 1920.75171 50715.44407
209 -297.17549 1920.75171
210 -15062.26497 -297.17549
211 -16673.38236 -15062.26497
212 -14721.91863 -16673.38236
213 46660.10862 -14721.91863
214 2950.66652 46660.10862
215 15868.46698 2950.66652
216 3357.44202 15868.46698
217 5128.96156 3357.44202
218 2459.02061 5128.96156
219 -7467.16835 2459.02061
220 52191.34466 -7467.16835
221 -355.77664 52191.34466
222 -358.57927 -355.77664
223 -18297.94332 -358.57927
224 -16160.76077 -18297.94332
225 68180.04213 -16160.76077
226 1308.21947 68180.04213
227 -8341.74193 1308.21947
228 47949.53529 -8341.74193
229 -17791.59277 47949.53529
230 6177.57183 -17791.59277
231 651.62846 6177.57183
232 -7252.99262 651.62846
233 -7835.72560 -7252.99262
234 47192.08494 -7835.72560
235 -19099.35521 47192.08494
236 -17749.62883 -19099.35521
237 75.52651 -17749.62883
238 -13948.26798 75.52651
239 50104.42190 -13948.26798
240 18065.37654 50104.42190
241 17387.64145 18065.37654
242 -5038.01339 17387.64145
243 -20816.06567 -5038.01339
244 -22115.81690 -20816.06567
245 -18765.16020 -22115.81690
246 -1267.01618 -18765.16020
247 45550.45403 -1267.01618
248 -15369.41567 45550.45403
249 -17617.46693 -15369.41567
250 -19400.81691 -17617.46693
251 -1893.38674 -19400.81691
252 61876.37714 -1893.38674
253 -3579.48902 61876.37714
254 -18008.92973 -3579.48902
255 -20187.87962 -18008.92973
256 -2524.41719 -20187.87962
257 44380.90666 -2524.41719
258 -5496.23620 44380.90666
259 -20875.01508 -5496.23620
260 -19771.34257 -20875.01508
261 -20814.19431 -19771.34257
262 -8061.80218 -20814.19431
263 -3202.55109 -8061.80218
264 -3537.99352 -3202.55109
265 39319.77330 -3537.99352
266 -8042.81988 39319.77330
267 -19553.18127 -8042.81988
268 -21717.80034 -19553.18127
269 -837.50405 -21717.80034
270 3871.05902 -837.50405
271 -4173.27937 3871.05902
272 -21168.80318 -4173.27937
273 44409.08172 -21168.80318
274 12138.12916 44409.08172
275 -8034.30096 12138.12916
276 -25347.39844 -8034.30096
277 -20229.30361 -25347.39844
278 -22976.36334 -20229.30361
279 -5435.30095 -22976.36334
280 41382.42858 -5435.30095
281 -24514.02999 41382.42858
282 -26146.59111 -24514.02999
283 -23751.18007 -26146.59111
284 -5700.66425 -23751.18007
285 -13779.98950 -5700.66425
286 -10510.31834 -13779.98950
287 43819.17241 -10510.31834
288 -15305.52162 43819.17241
289 -10255.17083 -15305.52162
290 -25458.54604 -10255.17083
291 111.78307 -25458.54604
292 -7028.57779 111.78307
293 -17000.41990 -7028.57779
294 -10879.52977 -17000.41990
295 -15360.00666 -10879.52977
296 2700.35048 -15360.00666
297 -24050.64828 2700.35048
298 -8384.68807 -24050.64828
299 38983.21516 -8384.68807
300 -26792.50610 38983.21516
301 -8218.78037 -26792.50610
302 38604.10516 -8218.78037
303 -27391.22282 38604.10516
304 -27409.66719 -27391.22282
305 -8716.85899 -27409.66719
306 38098.98796 -8716.85899
307 -27582.11523 38098.98796
308 -9092.85885 -27582.11523
309 37390.56217 -9092.85885
310 -9342.49803 37390.56217
311 -1502.78960 -9342.49803
312 -9235.35143 -1502.78960
313 -14005.48023 -9235.35143
314 -28577.40319 -14005.48023
315 -28726.68074 -28577.40319
316 -29992.71083 -28726.68074
317 -27931.13243 -29992.71083
318 -10366.01471 -27931.13243
319 36531.53903 -10366.01471
320 -14716.68987 36531.53903
321 -30901.31175 -14716.68987
322 -28725.17001 -30901.31175
323 -9403.26560 -28725.17001
324 -22356.84771 -9403.26560
325 -22095.20193 -22356.84771
326 -22395.61037 -22095.20193
327 -20693.27009 -22395.61037
328 -14290.21928 -20693.27009
329 35194.35413 -14290.21928
330 -27731.91095 35194.35413
331 -29578.04217 -27731.91095
332 -12126.59802 -29578.04217
333 34691.77317 -12126.59802
334 -9802.73191 34691.77317
335 -31105.47079 -9802.73191
336 -13776.48290 -31105.47079
337 -10631.80979 -13776.48290
338 -12880.37293 -10631.80979
339 33931.21906 -12880.37293
340 -31738.80115 33931.21906
341 -12056.07443 -31738.80115
342 -12937.21948 -12056.07443
343 -24616.08131 -12937.21948
344 -23213.63104 -24616.08131
345 -8241.84614 -23213.63104
346 33047.29505 -8241.84614
347 -35550.08400 33047.29505
348 -31140.79258 -35550.08400
349 -31425.35871 -31140.79258
350 -31967.14482 -31425.35871
351 -13994.75781 -31967.14482
352 -29742.15094 -13994.75781
353 -31288.82146 -29742.15094
354 34027.61022 -31288.82146
355 -10733.49685 34027.61022
356 -20353.61462 -10733.49685
357 -17601.18599 -20353.61462
358 -17550.87755 -17601.18599
359 -14537.81175 -17550.87755
360 -16831.91553 -14537.81175
361 -19240.14729 -16831.91553
362 -17819.41964 -19240.14729
363 -11234.68222 -17819.41964
364 6617.70766 -11234.68222
365 -21581.84772 6617.70766
366 -26977.38552 -21581.84772
367 39927.58782 -26977.38552
368 -17981.80119 39927.58782
369 -58646.59526 -17981.80119
370 -11388.33382 -58646.59526
371 -13822.81897 -11388.33382
372 -20258.44233 -13822.81897
373 -10599.93479 -20258.44233
374 -17654.98405 -10599.93479
375 -9714.81934 -17654.98405
376 -12943.62358 -9714.81934
377 -22570.45093 -12943.62358
378 -22929.42119 -22570.45093
379 -11463.76019 -22929.42119
380 -11088.50753 -11463.76019
381 -19155.80349 -11088.50753
382 -25608.70333 -19155.80349
383 -7308.63846 -25608.70333
384 -3226.43011 -7308.63846
385 -10514.98416 -3226.43011
386 -47994.81721 -10514.98416
387 -21574.78127 -47994.81721
388 -14012.95385 -21574.78127
389 -17934.51060 -14012.95385
390 -11550.59963 -17934.51060
391 -6325.45106 -11550.59963
392 -97914.37325 -6325.45106
393 -28104.92934 -97914.37325
394 -30551.49193 -28104.92934
395 -17038.43413 -30551.49193
396 -33475.29375 -17038.43413
397 75432.31451 -33475.29375
398 -5111.89402 75432.31451
399 -6361.50908 -5111.89402
400 -56484.65217 -6361.50908
401 145879.80998 -56484.65217
402 -8145.41065 145879.80998
403 71913.35675 -8145.41065
404 -110089.68199 71913.35675
405 47527.37604 -110089.68199
406 72877.11387 47527.37604
407 45254.61363 72877.11387
408 150445.11606 45254.61363
409 -26117.47344 150445.11606
410 -38978.33012 -26117.47344
411 160079.75443 -38978.33012
412 -38597.23496 160079.75443
413 -87656.23635 -38597.23496
414 -52573.07115 -87656.23635
415 10738.23220 -52573.07115
416 62576.78751 10738.23220
417 199667.35070 62576.78751
418 28516.23724 199667.35070
419 -47632.17632 28516.23724
420 34948.10466 -47632.17632
421 -4972.85033 34948.10466
422 -6984.54617 -4972.85033
423 78197.79206 -6984.54617
424 29044.13301 78197.79206
425 45755.64161 29044.13301
426 -23398.74984 45755.64161
427 -10880.62655 -23398.74984
428 -57808.15262 -10880.62655
429 -56167.22682 -57808.15262
430 -36555.74229 -56167.22682
> 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/7rggb1291403555.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/8rggb1291403555.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/9rggb1291403555.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/10k7xe1291403555.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/11nqek1291403555.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/129qcq1291403555.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/1350sz1291403555.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/14qj951291403555.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/15u1pt1291403555.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/168bnj1291403555.tab")
+ }
>
> try(system("convert tmp/1v60k1291403555.ps tmp/1v60k1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/2oyio1291403555.ps tmp/2oyio1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/3oyio1291403555.ps tmp/3oyio1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/4oyio1291403555.ps tmp/4oyio1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/5g7zq1291403555.ps tmp/5g7zq1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/6g7zq1291403555.ps tmp/6g7zq1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/7rggb1291403555.ps tmp/7rggb1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/8rggb1291403555.ps tmp/8rggb1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/9rggb1291403555.ps tmp/9rggb1291403555.png",intern=TRUE))
character(0)
> try(system("convert tmp/10k7xe1291403555.ps tmp/10k7xe1291403555.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.385 2.389 60.280