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
+ ,807
+ ,213118
+ ,6282154
+ ,1
+ ,5626
+ ,37
+ ,56289
+ ,29790
+ ,309
+ ,444
+ ,81767
+ ,4321023
+ ,2
+ ,13337
+ ,138
+ ,28328
+ ,87550
+ ,458
+ ,412
+ ,153198
+ ,4111912
+ ,3
+ ,8541
+ ,45
+ ,17936
+ ,84738
+ ,588
+ ,428
+ ,-26007
+ ,223193
+ ,415
+ ,39
+ ,0
+ ,4145
+ ,54660
+ ,302
+ ,315
+ ,126942
+ ,1491348
+ ,6
+ ,4276
+ ,24
+ ,3040
+ ,42634
+ ,156
+ ,168
+ ,157214
+ ,1629616
+ ,5
+ ,9164
+ ,34
+ ,2964
+ ,40949
+ ,481
+ ,263
+ ,129352
+ ,1398893
+ ,9
+ ,2493
+ ,29
+ ,2865
+ ,45187
+ ,353
+ ,267
+ ,234817
+ ,1926517
+ ,4
+ ,4891
+ ,38
+ ,2854
+ ,37704
+ ,452
+ ,228
+ ,60448
+ ,983660
+ ,14
+ ,1734
+ ,21
+ ,2167
+ ,16275
+ ,109
+ ,129
+ ,47818
+ ,1443586
+ ,7
+ ,11409
+ ,76
+ ,1974
+ ,25830
+ ,115
+ ,104
+ ,245546
+ ,1073089
+ ,10
+ ,7592
+ ,34
+ ,1910
+ ,12679
+ ,110
+ ,122
+ ,48020
+ ,984885
+ ,13
+ ,7135
+ ,62
+ ,1871
+ ,18014
+ ,239
+ ,393
+ ,-1710
+ ,1405225
+ ,8
+ ,5043
+ ,67
+ ,943
+ ,43556
+ ,247
+ ,190
+ ,32648
+ ,227132
+ ,414
+ ,110
+ ,1
+ ,929
+ ,24811
+ ,505
+ ,280
+ ,95350
+ ,929118
+ ,15
+ ,1444
+ ,29
+ ,822
+ ,6575
+ ,159
+ ,63
+ ,151352
+ ,1071292
+ ,11
+ ,5480
+ ,133
+ ,819
+ ,7123
+ ,109
+ ,102
+ ,288170
+ ,638830
+ ,28
+ ,4026
+ ,62
+ ,769
+ ,21950
+ ,519
+ ,265
+ ,114337
+ ,856956
+ ,17
+ ,1266
+ ,30
+ ,745
+ ,37597
+ ,248
+ ,234
+ ,37884
+ ,992426
+ ,12
+ ,3195
+ ,21
+ ,652
+ ,17821
+ ,373
+ ,277
+ ,122844
+ ,444477
+ ,46
+ ,655
+ ,14
+ ,643
+ ,12988
+ ,119
+ ,73
+ ,82340
+ ,857217
+ ,16
+ ,5523
+ ,51
+ ,601
+ ,22330
+ ,84
+ ,67
+ ,79801
+ ,711969
+ ,20
+ ,6095
+ ,23
+ ,446
+ ,13326
+ ,102
+ ,103
+ ,165548
+ ,702380
+ ,21
+ ,4925
+ ,38
+ ,436
+ ,16189
+ ,295
+ ,290
+ ,116384
+ ,358589
+ ,89
+ ,538
+ ,10
+ ,379
+ ,7146
+ ,105
+ ,83
+ ,134028
+ ,297978
+ ,367
+ ,933
+ ,14
+ ,305
+ ,15824
+ ,64
+ ,56
+ ,63838
+ ,585715
+ ,30
+ ,6027
+ ,24
+ ,284
+ ,27664
+ ,282
+ ,236
+ ,74996
+ ,657954
+ ,24
+ ,1624
+ ,17
+ ,247
+ ,11920
+ ,182
+ ,73
+ ,31080
+ ,209458
+ ,419
+ ,52
+ ,1
+ ,238
+ ,8568
+ ,37
+ ,34
+ ,32168
+ ,786690
+ ,18
+ ,15856
+ ,68
+ ,223
+ ,14416
+ ,361
+ ,139
+ ,49857
+ ,439798
+ ,48
+ ,664
+ ,17
+ ,220
+ ,3369
+ ,28
+ ,26
+ ,87161
+ ,688779
+ ,23
+ ,17456
+ ,145
+ ,217
+ ,11819
+ ,85
+ ,70
+ ,106113
+ ,574339
+ ,31
+ ,4404
+ ,32
+ ,199
+ ,6984
+ ,45
+ ,40
+ ,80570
+ ,741409
+ ,19
+ ,12031
+ ,78
+ ,195
+ ,4519
+ ,49
+ ,42
+ ,102129
+ ,597793
+ ,29
+ ,8118
+ ,88
+ ,163
+ ,2220
+ ,22
+ ,12
+ ,301670
+ ,644190
+ ,26
+ ,20190
+ ,200
+ ,154
+ ,18562
+ ,155
+ ,211
+ ,102313
+ ,377934
+ ,72
+ ,1148
+ ,10
+ ,143
+ ,10327
+ ,91
+ ,74
+ ,88577
+ ,640273
+ ,27
+ ,4838
+ ,43
+ ,135
+ ,5336
+ ,81
+ ,80
+ ,112477
+ ,697458
+ ,22
+ ,6141
+ ,93
+ ,132
+ ,2365
+ ,79
+ ,83
+ ,191778
+ ,550608
+ ,32
+ ,4438
+ ,148
+ ,120
+ ,4069
+ ,145
+ ,131
+ ,79804
+ ,207393
+ ,420
+ ,51
+ ,2
+ ,119
+ ,8636
+ ,855
+ ,203
+ ,128294
+ ,301607
+ ,359
+ ,119
+ ,12
+ ,108
+ ,13718
+ ,61
+ ,56
+ ,96448
+ ,345783
+ ,100
+ ,2390
+ ,11
+ ,104
+ ,4525
+ ,226
+ ,89
+ ,93811
+ ,501749
+ ,37
+ ,1335
+ ,67
+ ,101
+ ,6869
+ ,105
+ ,88
+ ,117520
+ ,379983
+ ,68
+ ,1714
+ ,26
+ ,90
+ ,4628
+ ,62
+ ,39
+ ,69159
+ ,387475
+ ,64
+ ,3024
+ ,41
+ ,85
+ ,3689
+ ,25
+ ,25
+ ,101792
+ ,377305
+ ,74
+ ,7092
+ ,48
+ ,74
+ ,4891
+ ,217
+ ,49
+ ,210568
+ ,370837
+ ,79
+ ,787
+ ,35
+ ,73
+ ,7489
+ ,322
+ ,149
+ ,136996
+ ,430866
+ ,50
+ ,717
+ ,31
+ ,70
+ ,4901
+ ,84
+ ,58
+ ,121920
+ ,469107
+ ,41
+ ,3204
+ ,55
+ ,67
+ ,2284
+ ,33
+ ,41
+ ,76403
+ ,194493
+ ,423
+ ,-167
+ ,-2
+ ,66
+ ,3160
+ ,108
+ ,90
+ ,108094
+ ,530670
+ ,33
+ ,3062
+ ,105
+ ,66
+ ,4150
+ ,150
+ ,136
+ ,134759
+ ,518365
+ ,35
+ ,2122
+ ,77
+ ,66
+ ,7285
+ ,115
+ ,97
+ ,188873
+ ,491303
+ ,39
+ ,2533
+ ,40
+ ,59
+ ,1134
+ ,162
+ ,63
+ ,146216
+ ,527021
+ ,34
+ ,2019
+ ,288
+ ,58
+ ,4658
+ ,158
+ ,114
+ ,156608
+ ,233773
+ ,411
+ ,214
+ ,7
+ ,58
+ ,2384
+ ,97
+ ,77
+ ,61348
+ ,405972
+ ,58
+ ,2123
+ ,86
+ ,54
+ ,3748
+ ,9
+ ,6
+ ,50350
+ ,652925
+ ,25
+ ,50325
+ ,121
+ ,53
+ ,5371
+ ,66
+ ,47
+ ,87720
+ ,446211
+ ,45
+ ,3730
+ ,46
+ ,49
+ ,1285
+ ,107
+ ,51
+ ,99489
+ ,341340
+ ,106
+ ,1321
+ ,110
+ ,49
+ ,9327
+ ,101
+ ,85
+ ,87419
+ ,387699
+ ,63
+ ,1858
+ ,20
+ ,44
+ ,5565
+ ,47
+ ,43
+ ,94355
+ ,493408
+ ,38
+ ,6243
+ ,53
+ ,39
+ ,1528
+ ,38
+ ,32
+ ,60326
+ ,146494
+ ,426
+ ,-1408
+ ,-35
+ ,38
+ ,3122
+ ,34
+ ,25
+ ,94670
+ ,414462
+ ,54
+ ,6308
+ ,69
+ ,38
+ ,7561
+ ,87
+ ,77
+ ,82425
+ ,364304
+ ,86
+ ,1889
+ ,22
+ ,37
+ ,2675
+ ,79
+ ,54
+ ,59017
+ ,355178
+ ,93
+ ,1964
+ ,58
+ ,36
+ ,13253
+ ,947
+ ,251
+ ,90829
+ ,357760
+ ,90
+ ,167
+ ,12
+ ,35
+ ,880
+ ,74
+ ,15
+ ,80791
+ ,261216
+ ,399
+ ,827
+ ,70
+ ,34
+ ,2053
+ ,53
+ ,44
+ ,100423
+ ,397144
+ ,60
+ ,3720
+ ,96
+ ,33
+ ,1424
+ ,94
+ ,73
+ ,131116
+ ,374943
+ ,77
+ ,1861
+ ,123
+ ,32
+ ,4036
+ ,63
+ ,85
+ ,100269
+ ,424898
+ ,53
+ ,3570
+ ,56
+ ,32
+ ,3045
+ ,58
+ ,49
+ ,27330
+ ,202055
+ ,422
+ ,35
+ ,1
+ ,31
+ ,5119
+ ,49
+ ,38
+ ,39039
+ ,378525
+ ,69
+ ,3643
+ ,35
+ ,31
+ ,1431
+ ,34
+ ,35
+ ,106885
+ ,310768
+ ,330
+ ,3258
+ ,77
+ ,30
+ ,554
+ ,11
+ ,9
+ ,79285
+ ,325738
+ ,125
+ ,11431
+ ,227
+ ,30
+ ,1975
+ ,35
+ ,34
+ ,118881
+ ,394510
+ ,61
+ ,5557
+ ,98
+ ,30
+ ,1765
+ ,20
+ ,20
+ ,77623
+ ,247060
+ ,404
+ ,2353
+ ,27
+ ,28
+ ,1012
+ ,47
+ ,29
+ ,114768
+ ,368078
+ ,82
+ ,3576
+ ,166
+ ,26
+ ,810
+ ,43
+ ,11
+ ,74015
+ ,236761
+ ,409
+ ,855
+ ,45
+ ,26
+ ,1280
+ ,117
+ ,52
+ ,69465
+ ,312378
+ ,326
+ ,960
+ ,88
+ ,26
+ ,666
+ ,171
+ ,13
+ ,117869
+ ,339836
+ ,108
+ ,818
+ ,210
+ ,25
+ ,1380
+ ,26
+ ,29
+ ,60982
+ ,347385
+ ,98
+ ,5669
+ ,107
+ ,25
+ ,4677
+ ,75
+ ,66
+ ,90131
+ ,426280
+ ,52
+ ,3017
+ ,48
+ ,24
+ ,876
+ ,59
+ ,33
+ ,138971
+ ,352850
+ ,96
+ ,2591
+ ,174
+ ,24
+ ,814
+ ,18
+ ,15
+ ,39625
+ ,301881
+ ,358
+ ,5660
+ ,125
+ ,23
+ ,514
+ ,15
+ ,15
+ ,102725
+ ,377516
+ ,73
+ ,11834
+ ,346
+ ,23
+ ,5692
+ ,72
+ ,68
+ ,64239
+ ,357312
+ ,91
+ ,2185
+ ,28
+ ,22
+ ,3642
+ ,86
+ ,100
+ ,90262
+ ,458343
+ ,42
+ ,3004
+ ,71
+ ,22
+ ,540
+ ,14
+ ,13
+ ,103960
+ ,354228
+ ,94
+ ,11016
+ ,286
+ ,22
+ ,2099
+ ,64
+ ,45
+ ,106611
+ ,308636
+ ,339
+ ,1697
+ ,52
+ ,20
+ ,567
+ ,11
+ ,14
+ ,103345
+ ,386212
+ ,65
+ ,16928
+ ,329
+ ,20
+ ,2001
+ ,52
+ ,36
+ ,95551
+ ,393343
+ ,62
+ ,3718
+ ,97
+ ,19
+ ,2949
+ ,41
+ ,40
+ ,82903
+ ,378509
+ ,70
+ ,4354
+ ,61
+ ,18
+ ,2253
+ ,99
+ ,68
+ ,63593
+ ,452469
+ ,43
+ ,2550
+ ,112
+ ,17
+ ,6533
+ ,75
+ ,29
+ ,126910
+ ,364839
+ ,85
+ ,2198
+ ,25
+ ,16
+ ,1889
+ ,45
+ ,43
+ ,37527
+ ,358649
+ ,88
+ ,3526
+ ,84
+ ,14
+ ,3055
+ ,43
+ ,30
+ ,60247
+ ,376641
+ ,75
+ ,4108
+ ,58
+ ,14
+ ,272
+ ,8
+ ,9
+ ,112995
+ ,429112
+ ,51
+ ,28639
+ ,844
+ ,14
+ ,1414
+ ,198
+ ,22
+ ,70184
+ ,330546
+ ,119
+ ,659
+ ,92
+ ,13
+ ,2564
+ ,22
+ ,19
+ ,130140
+ ,403560
+ ,59
+ ,9253
+ ,79
+ ,13
+ ,1383
+ ,11
+ ,9
+ ,73221
+ ,317892
+ ,151
+ ,10717
+ ,85
+ ,13
+ ,1261
+ ,33
+ ,31
+ ,76114
+ ,307528
+ ,343
+ ,3258
+ ,85
+ ,13
+ ,975
+ ,23
+ ,19
+ ,90534
+ ,235133
+ ,410
+ ,1528
+ ,36
+ ,13
+ ,3366
+ ,80
+ ,55
+ ,108479
+ ,299243
+ ,362
+ ,1241
+ ,29
+ ,12
+ ,576
+ ,18
+ ,8
+ ,113761
+ ,314073
+ ,313
+ ,6337
+ ,198
+ ,12
+ ,1686
+ ,40
+ ,28
+ ,68696
+ ,368186
+ ,81
+ ,4205
+ ,100
+ ,11
+ ,746
+ ,23
+ ,29
+ ,71561
+ ,269661
+ ,396
+ ,3029
+ ,93
+ ,10
+ ,3192
+ ,60
+ ,48
+ ,59831
+ ,125390
+ ,427
+ ,-1243
+ ,-23
+ ,10
+ ,2045
+ ,20
+ ,16
+ ,97890
+ ,510834
+ ,36
+ ,15542
+ ,152
+ ,10
+ ,5702
+ ,61
+ ,47
+ ,101481
+ ,321896
+ ,138
+ ,1998
+ ,21
+ ,10
+ ,1932
+ ,36
+ ,20
+ ,72954
+ ,249898
+ ,403
+ ,1386
+ ,26
+ ,10
+ ,936
+ ,30
+ ,22
+ ,67939
+ ,408881
+ ,56
+ ,6963
+ ,223
+ ,9
+ ,3437
+ ,47
+ ,33
+ ,48022
+ ,158492
+ ,425
+ ,-883
+ ,-12
+ ,9
+ ,5131
+ ,71
+ ,44
+ ,86111
+ ,292154
+ ,379
+ ,1298
+ ,18
+ ,9
+ ,2397
+ ,14
+ ,13
+ ,74020
+ ,289513
+ ,385
+ ,6394
+ ,37
+ ,6
+ ,1389
+ ,9
+ ,6
+ ,57530
+ ,378049
+ ,71
+ ,19783
+ ,128
+ ,6
+ ,1503
+ ,39
+ ,35
+ ,56364
+ ,343466
+ ,104
+ ,3679
+ ,95
+ ,6
+ ,402
+ ,26
+ ,8
+ ,84990
+ ,332743
+ ,115
+ ,5106
+ ,330
+ ,5
+ ,2239
+ ,21
+ ,17
+ ,88590
+ ,442882
+ ,47
+ ,11566
+ ,108
+ ,5
+ ,2234
+ ,16
+ ,11
+ ,77200
+ ,214215
+ ,418
+ ,888
+ ,6
+ ,5
+ ,837
+ ,69
+ ,21
+ ,61262
+ ,315688
+ ,165
+ ,1677
+ ,138
+ ,4
+ ,10579
+ ,92
+ ,92
+ ,110309
+ ,375195
+ ,76
+ ,1904
+ ,17
+ ,4
+ ,875
+ ,14
+ ,12
+ ,67000
+ ,334280
+ ,112
+ ,9591
+ ,153
+ ,3
+ ,1585
+ ,107
+ ,112
+ ,93099
+ ,355864
+ ,92
+ ,1457
+ ,98
+ ,2
+ ,1659
+ ,29
+ ,25
+ ,107577
+ ,480382
+ ,40
+ ,9668
+ ,169
+ ,2
+ ,2647
+ ,37
+ ,17
+ ,62920
+ ,353058
+ ,95
+ ,4137
+ ,58
+ ,1
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,217193
+ ,416
+ ,748
+ ,5
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,191
+ ,0
+ ,94
+ ,7
+ ,10
+ ,60793
+ ,314533
+ ,308
+ ,16362
+ ,1212
+ ,0
+ ,422
+ ,28
+ ,23
+ ,57935
+ ,318056
+ ,150
+ ,4216
+ ,280
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,197
+ ,0
+ ,34
+ ,8
+ ,7
+ ,60630
+ ,314353
+ ,309
+ ,14294
+ ,3340
+ ,0
+ ,1558
+ ,63
+ ,25
+ ,55637
+ ,369448
+ ,80
+ ,2690
+ ,109
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,196
+ ,0
+ ,43
+ ,3
+ ,20
+ ,60887
+ ,312846
+ ,323
+ ,37615
+ ,2613
+ ,0
+ ,645
+ ,5
+ ,4
+ ,60720
+ ,312075
+ ,329
+ ,22415
+ ,174
+ ,0
+ ,316
+ ,9
+ ,4
+ ,60505
+ ,315009
+ ,299
+ ,12779
+ ,364
+ ,0
+ ,115
+ ,13
+ ,10
+ ,60945
+ ,318903
+ ,145
+ ,9146
+ ,1037
+ ,0
+ ,5
+ ,2
+ ,1
+ ,60720
+ ,314887
+ ,304
+ ,57443
+ ,24814
+ ,0
+ ,897
+ ,5
+ ,4
+ ,60720
+ ,314913
+ ,303
+ ,22983
+ ,128
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,188
+ ,0
+ ,389
+ ,14
+ ,8
+ ,58990
+ ,325506
+ ,127
+ ,8965
+ ,322
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,220
+ ,0
+ ,1002
+ ,15
+ ,11
+ ,56750
+ ,298568
+ ,364
+ ,6571
+ ,98
+ ,0
+ ,36
+ ,3
+ ,4
+ ,60894
+ ,315834
+ ,162
+ ,38611
+ ,3215
+ ,0
+ ,460
+ ,15
+ ,15
+ ,63346
+ ,329784
+ ,122
+ ,8652
+ ,282
+ ,0
+ ,309
+ ,11
+ ,9
+ ,56535
+ ,312878
+ ,322
+ ,10262
+ ,365
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,0
+ ,9
+ ,6
+ ,7
+ ,60835
+ ,314987
+ ,300
+ ,19165
+ ,12458
+ ,0
+ ,271
+ ,2
+ ,2
+ ,60720
+ ,325249
+ ,129
+ ,62625
+ ,463
+ ,0
+ ,14
+ ,1
+ ,0
+ ,61016
+ ,315877
+ ,161
+ ,115877
+ ,8041
+ ,0
+ ,520
+ ,10
+ ,7
+ ,58650
+ ,291650
+ ,383
+ ,9165
+ ,176
+ ,0
+ ,1766
+ ,73
+ ,46
+ ,60438
+ ,305959
+ ,349
+ ,1451
+ ,60
+ ,0
+ ,0
+ ,0
+ ,5
+ ,60720
+ ,315380
+ ,208
+ ,0
+ ,458
+ ,11
+ ,7
+ ,58625
+ ,297765
+ ,368
+ ,8888
+ ,213
+ ,0
+ ,20
+ ,3
+ ,2
+ ,60938
+ ,315245
+ ,297
+ ,38415
+ ,5674
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,217
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,219
+ ,0
+ ,98
+ ,2
+ ,2
+ ,61490
+ ,315236
+ ,298
+ ,57618
+ ,1179
+ ,0
+ ,405
+ ,7
+ ,5
+ ,60845
+ ,336425
+ ,111
+ ,19489
+ ,337
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,213
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,212
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,207
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,218
+ ,0
+ ,483
+ ,27
+ ,7
+ ,60830
+ ,306268
+ ,348
+ ,3936
+ ,220
+ ,0
+ ,454
+ ,51
+ ,24
+ ,63261
+ ,302187
+ ,357
+ ,2004
+ ,225
+ ,0
+ ,47
+ ,3
+ ,1
+ ,60720
+ ,314882
+ ,305
+ ,38294
+ ,2432
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,193
+ ,0
+ ,757
+ ,19
+ ,18
+ ,45689
+ ,382712
+ ,66
+ ,9616
+ ,242
+ ,0
+ ,4655
+ ,393
+ ,55
+ ,60720
+ ,341570
+ ,105
+ ,360
+ ,30
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,192
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,211
+ ,0
+ ,36
+ ,4
+ ,3
+ ,61564
+ ,312412
+ ,325
+ ,28103
+ ,3146
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,189
+ ,0
+ ,203
+ ,9
+ ,9
+ ,61938
+ ,309596
+ ,334
+ ,12177
+ ,541
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,205
+ ,0
+ ,126
+ ,10
+ ,8
+ ,60951
+ ,315547
+ ,171
+ ,11555
+ ,921
+ ,0
+ ,400
+ ,152
+ ,113
+ ,60720
+ ,313267
+ ,319
+ ,745
+ ,283
+ ,0
+ ,71
+ ,1
+ ,0
+ ,60745
+ ,316176
+ ,159
+ ,116176
+ ,1646
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,209
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,190
+ ,0
+ ,972
+ ,34
+ ,19
+ ,71642
+ ,359335
+ ,87
+ ,4686
+ ,164
+ ,0
+ ,531
+ ,10
+ ,11
+ ,71641
+ ,330068
+ ,120
+ ,13007
+ ,245
+ ,0
+ ,2461
+ ,57
+ ,25
+ ,55792
+ ,314289
+ ,310
+ ,2005
+ ,46
+ ,0
+ ,378
+ ,52
+ ,16
+ ,71873
+ ,297413
+ ,369
+ ,1873
+ ,258
+ ,0
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,314806
+ ,306
+ ,22961
+ ,4974
+ ,0
+ ,638
+ ,14
+ ,11
+ ,60370
+ ,333210
+ ,114
+ ,9515
+ ,209
+ ,0
+ ,2300
+ ,29
+ ,23
+ ,64873
+ ,352108
+ ,97
+ ,5245
+ ,66
+ ,0
+ ,149
+ ,5
+ ,6
+ ,62041
+ ,313332
+ ,318
+ ,22666
+ ,761
+ ,0
+ ,226
+ ,5
+ ,5
+ ,65745
+ ,291787
+ ,382
+ ,18357
+ ,406
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,184
+ ,0
+ ,275
+ ,4
+ ,7
+ ,59500
+ ,318745
+ ,146
+ ,29686
+ ,431
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,179
+ ,0
+ ,141
+ ,6
+ ,7
+ ,61630
+ ,315366
+ ,295
+ ,19228
+ ,819
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,186
+ ,0
+ ,28
+ ,2
+ ,3
+ ,60890
+ ,315688
+ ,166
+ ,57844
+ ,4148
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,263
+ ,0
+ ,4980
+ ,91
+ ,89
+ ,113521
+ ,409642
+ ,55
+ ,2304
+ ,42
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,259
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,261
+ ,0
+ ,472
+ ,20
+ ,19
+ ,80045
+ ,269587
+ ,397
+ ,3479
+ ,148
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,243
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,250
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,247
+ ,0
+ ,203
+ ,27
+ ,12
+ ,50804
+ ,300962
+ ,361
+ ,3739
+ ,498
+ ,0
+ ,496
+ ,17
+ ,12
+ ,87390
+ ,325479
+ ,128
+ ,7381
+ ,253
+ ,0
+ ,10
+ ,2
+ ,5
+ ,61656
+ ,316155
+ ,160
+ ,58077
+ ,11639
+ ,0
+ ,63
+ ,4
+ ,2
+ ,65688
+ ,318574
+ ,148
+ ,29643
+ ,1891
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,228
+ ,0
+ ,1136
+ ,32
+ ,26
+ ,48522
+ ,343613
+ ,103
+ ,4488
+ ,126
+ ,0
+ ,265
+ ,31
+ ,3
+ ,60720
+ ,306948
+ ,346
+ ,3450
+ ,404
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,210
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,206
+ ,0
+ ,267
+ ,32
+ ,11
+ ,57640
+ ,330059
+ ,121
+ ,4064
+ ,487
+ ,0
+ ,474
+ ,20
+ ,10
+ ,61977
+ ,288985
+ ,386
+ ,4449
+ ,188
+ ,0
+ ,534
+ ,7
+ ,5
+ ,62620
+ ,304485
+ ,351
+ ,14926
+ ,196
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,201
+ ,0
+ ,15
+ ,8
+ ,6
+ ,60831
+ ,315688
+ ,164
+ ,14461
+ ,7687
+ ,0
+ ,397
+ ,28
+ ,7
+ ,60646
+ ,317736
+ ,152
+ ,4205
+ ,297
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,199
+ ,0
+ ,1866
+ ,29
+ ,28
+ ,56225
+ ,322331
+ ,135
+ ,4218
+ ,66
+ ,0
+ ,288
+ ,4
+ ,3
+ ,60510
+ ,296656
+ ,373
+ ,24164
+ ,335
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,198
+ ,0
+ ,3
+ ,2
+ ,1
+ ,60698
+ ,315354
+ ,296
+ ,57677
+ ,36505
+ ,0
+ ,468
+ ,21
+ ,20
+ ,60720
+ ,312161
+ ,328
+ ,5341
+ ,239
+ ,0
+ ,20
+ ,2
+ ,1
+ ,60805
+ ,315576
+ ,169
+ ,57788
+ ,5924
+ ,0
+ ,278
+ ,26
+ ,22
+ ,61404
+ ,314922
+ ,302
+ ,4420
+ ,413
+ ,0
+ ,61
+ ,14
+ ,9
+ ,60720
+ ,314551
+ ,307
+ ,8182
+ ,1874
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,215
+ ,0
+ ,192
+ ,4
+ ,2
+ ,65276
+ ,312339
+ ,327
+ ,28085
+ ,586
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,214
+ ,0
+ ,317
+ ,9
+ ,7
+ ,63915
+ ,298700
+ ,363
+ ,10967
+ ,311
+ ,0
+ ,738
+ ,10
+ ,9
+ ,60720
+ ,321376
+ ,139
+ ,12138
+ ,164
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,216
+ ,0
+ ,368
+ ,17
+ ,13
+ ,61686
+ ,303230
+ ,355
+ ,6072
+ ,281
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,204
+ ,0
+ ,2
+ ,1
+ ,0
+ ,60743
+ ,315487
+ ,172
+ ,115487
+ ,52257
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,202
+ ,0
+ ,53
+ ,6
+ ,6
+ ,60349
+ ,315793
+ ,163
+ ,19299
+ ,2197
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,195
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,194
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,200
+ ,0
+ ,94
+ ,3
+ ,3
+ ,61360
+ ,312887
+ ,321
+ ,37629
+ ,1200
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,237
+ ,0
+ ,24
+ ,8
+ ,7
+ ,59818
+ ,315637
+ ,168
+ ,14455
+ ,4772
+ ,0
+ ,2332
+ ,4
+ ,2
+ ,72680
+ ,324385
+ ,131
+ ,31096
+ ,53
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,236
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,233
+ ,0
+ ,131
+ ,11
+ ,15
+ ,61808
+ ,308989
+ ,338
+ ,9908
+ ,834
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,260
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,234
+ ,0
+ ,206
+ ,9
+ ,9
+ ,53110
+ ,296702
+ ,372
+ ,10745
+ ,469
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,240
+ ,0
+ ,167
+ ,2
+ ,1
+ ,64245
+ ,307322
+ ,345
+ ,53661
+ ,642
+ ,0
+ ,622
+ ,73
+ ,38
+ ,73007
+ ,304376
+ ,352
+ ,1430
+ ,168
+ ,0
+ ,2328
+ ,94
+ ,57
+ ,82732
+ ,253588
+ ,400
+ ,570
+ ,23
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,238
+ ,0
+ ,365
+ ,8
+ ,7
+ ,54820
+ ,309560
+ ,335
+ ,13695
+ ,300
+ ,0
+ ,364
+ ,35
+ ,26
+ ,47705
+ ,298466
+ ,365
+ ,2813
+ ,271
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,239
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,235
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,262
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,283
+ ,0
+ ,226
+ ,12
+ ,13
+ ,72835
+ ,343929
+ ,102
+ ,11994
+ ,636
+ ,0
+ ,307
+ ,15
+ ,10
+ ,58856
+ ,331955
+ ,116
+ ,8797
+ ,430
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,288
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,281
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,282
+ ,0
+ ,188
+ ,11
+ ,9
+ ,77655
+ ,381180
+ ,67
+ ,16471
+ ,964
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,221
+ ,0
+ ,138
+ ,6
+ ,26
+ ,69817
+ ,331420
+ ,117
+ ,21903
+ ,954
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,286
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,285
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,265
+ ,0
+ ,125
+ ,12
+ ,19
+ ,60798
+ ,310201
+ ,332
+ ,9183
+ ,883
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,284
+ ,0
+ ,282
+ ,30
+ ,12
+ ,62452
+ ,320016
+ ,141
+ ,4001
+ ,425
+ ,0
+ ,335
+ ,33
+ ,23
+ ,64175
+ ,320398
+ ,140
+ ,3648
+ ,359
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,232
+ ,0
+ ,1324
+ ,117
+ ,29
+ ,67440
+ ,291841
+ ,381
+ ,785
+ ,69
+ ,0
+ ,176
+ ,28
+ ,8
+ ,68136
+ ,310670
+ ,331
+ ,3952
+ ,629
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,256
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,255
+ ,0
+ ,249
+ ,72
+ ,26
+ ,56726
+ ,313491
+ ,317
+ ,1576
+ ,455
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,242
+ ,0
+ ,333
+ ,13
+ ,9
+ ,70811
+ ,331323
+ ,118
+ ,10102
+ ,394
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,241
+ ,0
+ ,601
+ ,6
+ ,5
+ ,60720
+ ,319210
+ ,144
+ ,19868
+ ,198
+ ,0
+ ,30
+ ,4
+ ,3
+ ,62045
+ ,318098
+ ,149
+ ,29525
+ ,3912
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,254
+ ,0
+ ,249
+ ,62
+ ,13
+ ,54323
+ ,292754
+ ,378
+ ,1496
+ ,373
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,258
+ ,0
+ ,165
+ ,24
+ ,12
+ ,62841
+ ,325176
+ ,130
+ ,5216
+ ,759
+ ,0
+ ,453
+ ,21
+ ,19
+ ,81125
+ ,365959
+ ,84
+ ,7903
+ ,366
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,257
+ ,0
+ ,53
+ ,14
+ ,10
+ ,59506
+ ,302409
+ ,356
+ ,7315
+ ,1932
+ ,0
+ ,382
+ ,21
+ ,9
+ ,59365
+ ,340968
+ ,107
+ ,6713
+ ,369
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,229
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,230
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,224
+ ,0
+ ,0
+ ,0
+ ,9
+ ,60720
+ ,315380
+ ,223
+ ,0
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,313164
+ ,320
+ ,28291
+ ,3808
+ ,0
+ ,290
+ ,2
+ ,1
+ ,58790
+ ,301164
+ ,360
+ ,50582
+ ,348
+ ,0
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,315380
+ ,225
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,231
+ ,0
+ ,366
+ ,53
+ ,14
+ ,61808
+ ,344425
+ ,101
+ ,2725
+ ,394
+ ,0
+ ,2
+ ,9
+ ,12
+ ,60735
+ ,315394
+ ,174
+ ,12822
+ ,63057
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,222
+ ,0
+ ,209
+ ,13
+ ,19
+ ,64016
+ ,316647
+ ,157
+ ,8973
+ ,558
+ ,0
+ ,384
+ ,22
+ ,17
+ ,54683
+ ,309836
+ ,333
+ ,4993
+ ,286
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,227
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,244
+ ,0
+ ,365
+ ,83
+ ,32
+ ,87192
+ ,346611
+ ,99
+ ,1766
+ ,402
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,226
+ ,0
+ ,49
+ ,8
+ ,14
+ ,64107
+ ,322031
+ ,137
+ ,15254
+ ,2500
+ ,0
+ ,3
+ ,4
+ ,8
+ ,60761
+ ,315656
+ ,167
+ ,28914
+ ,37920
+ ,0
+ ,133
+ ,14
+ ,4
+ ,65990
+ ,339445
+ ,109
+ ,9960
+ ,1052
+ ,0
+ ,32
+ ,1
+ ,0
+ ,59988
+ ,314964
+ ,301
+ ,114964
+ ,3603
+ ,0
+ ,368
+ ,17
+ ,20
+ ,61167
+ ,297141
+ ,370
+ ,5714
+ ,264
+ ,0
+ ,1
+ ,6
+ ,5
+ ,60719
+ ,315372
+ ,293
+ ,19229
+ ,84213
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,274
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,246
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,177
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,185
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,245
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,273
+ ,0
+ ,22
+ ,2
+ ,1
+ ,60722
+ ,312502
+ ,324
+ ,56251
+ ,5154
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,252
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,253
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,277
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,278
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,275
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,276
+ ,0
+ ,96
+ ,5
+ ,4
+ ,60379
+ ,313729
+ ,316
+ ,22746
+ ,1188
+ ,0
+ ,1
+ ,2
+ ,1
+ ,60727
+ ,315388
+ ,175
+ ,57694
+ ,117743
+ ,0
+ ,314
+ ,5
+ ,4
+ ,60720
+ ,315371
+ ,294
+ ,23074
+ ,367
+ ,0
+ ,844
+ ,78
+ ,20
+ ,60925
+ ,296139
+ ,375
+ ,1233
+ ,114
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,180
+ ,0
+ ,26
+ ,1
+ ,1
+ ,60896
+ ,313880
+ ,315
+ ,113880
+ ,4328
+ ,0
+ ,125
+ ,13
+ ,10
+ ,59734
+ ,317698
+ ,153
+ ,9054
+ ,944
+ ,0
+ ,304
+ ,15
+ ,12
+ ,62969
+ ,295580
+ ,376
+ ,6372
+ ,314
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,187
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,183
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,182
+ ,0
+ ,621
+ ,48
+ ,13
+ ,60720
+ ,308256
+ ,341
+ ,2255
+ ,174
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,181
+ ,0
+ ,119
+ ,6
+ ,3
+ ,59118
+ ,303677
+ ,354
+ ,17280
+ ,872
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,178
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,176
+ ,0
+ ,1595
+ ,17
+ ,10
+ ,60720
+ ,319369
+ ,143
+ ,7022
+ ,75
+ ,0
+ ,312
+ ,14
+ ,3
+ ,58598
+ ,318690
+ ,147
+ ,8478
+ ,381
+ ,0
+ ,60
+ ,10
+ ,7
+ ,61124
+ ,314049
+ ,314
+ ,11405
+ ,1916
+ ,0
+ ,587
+ ,12
+ ,10
+ ,59595
+ ,325699
+ ,126
+ ,10475
+ ,214
+ ,0
+ ,135
+ ,2
+ ,1
+ ,62065
+ ,314210
+ ,311
+ ,57105
+ ,847
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,272
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,269
+ ,0
+ ,514
+ ,52
+ ,15
+ ,78780
+ ,322378
+ ,133
+ ,2353
+ ,238
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,279
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,251
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,270
+ ,0
+ ,1
+ ,4
+ ,4
+ ,60722
+ ,315398
+ ,173
+ ,28849
+ ,172235
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,268
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,289
+ ,0
+ ,1763
+ ,24
+ ,28
+ ,61600
+ ,308336
+ ,340
+ ,4514
+ ,61
+ ,0
+ ,180
+ ,11
+ ,9
+ ,59635
+ ,316386
+ ,158
+ ,10581
+ ,647
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,287
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,248
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,291
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,292
+ ,0
+ ,218
+ ,21
+ ,7
+ ,60720
+ ,315553
+ ,170
+ ,5503
+ ,531
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,290
+ ,0
+ ,448
+ ,40
+ ,7
+ ,59781
+ ,323361
+ ,132
+ ,3084
+ ,275
+ ,0
+ ,227
+ ,9
+ ,7
+ ,76644
+ ,336639
+ ,110
+ ,15182
+ ,603
+ ,0
+ ,174
+ ,1
+ ,3
+ ,64820
+ ,307424
+ ,344
+ ,107424
+ ,618
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,266
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,280
+ ,0
+ ,121
+ ,24
+ ,11
+ ,56178
+ ,295370
+ ,377
+ ,3974
+ ,788
+ ,0
+ ,607
+ ,11
+ ,7
+ ,60436
+ ,322340
+ ,134
+ ,11122
+ ,201
+ ,0
+ ,2212
+ ,14
+ ,10
+ ,60720
+ ,319864
+ ,142
+ ,8562
+ ,54
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,271
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,264
+ ,0
+ ,530
+ ,60
+ ,18
+ ,73433
+ ,317291
+ ,155
+ ,1955
+ ,221
+ ,0
+ ,571
+ ,80
+ ,14
+ ,41477
+ ,280398
+ ,390
+ ,1005
+ ,141
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,267
+ ,0
+ ,78
+ ,16
+ ,12
+ ,62700
+ ,317330
+ ,154
+ ,7333
+ ,1511
+ ,0
+ ,2489
+ ,40
+ ,29
+ ,67804
+ ,238125
+ ,408
+ ,953
+ ,15
+ ,0
+ ,131
+ ,6
+ ,3
+ ,59661
+ ,327071
+ ,123
+ ,21179
+ ,969
+ ,0
+ ,923
+ ,8
+ ,6
+ ,58620
+ ,309038
+ ,337
+ ,13630
+ ,118
+ ,0
+ ,72
+ ,3
+ ,3
+ ,60398
+ ,314210
+ ,312
+ ,38070
+ ,1591
+ ,0
+ ,572
+ ,16
+ ,8
+ ,58580
+ ,307930
+ ,342
+ ,6746
+ ,189
+ ,-1
+ ,397
+ ,10
+ ,10
+ ,62710
+ ,322327
+ ,136
+ ,12233
+ ,308
+ ,-1
+ ,450
+ ,8
+ ,6
+ ,59325
+ ,292136
+ ,380
+ ,11517
+ ,205
+ ,-1
+ ,622
+ ,7
+ ,8
+ ,60950
+ ,263276
+ ,398
+ ,9039
+ ,102
+ ,-2
+ ,694
+ ,8
+ ,6
+ ,68060
+ ,367655
+ ,83
+ ,20957
+ ,241
+ ,-2
+ ,3425
+ ,12
+ ,9
+ ,83620
+ ,283910
+ ,388
+ ,6993
+ ,25
+ ,-3
+ ,562
+ ,13
+ ,8
+ ,58456
+ ,283587
+ ,389
+ ,6430
+ ,149
+ ,-3
+ ,4917
+ ,42
+ ,26
+ ,52811
+ ,243650
+ ,406
+ ,1039
+ ,9
+ ,-3
+ ,1442
+ ,118
+ ,239
+ ,121173
+ ,438493
+ ,49
+ ,2021
+ ,165
+ ,-3
+ ,529
+ ,9
+ ,7
+ ,63870
+ ,296261
+ ,374
+ ,10696
+ ,182
+ ,-3
+ ,2126
+ ,138
+ ,41
+ ,21001
+ ,230621
+ ,413
+ ,222
+ ,14
+ ,-3
+ ,1061
+ ,5
+ ,3
+ ,70415
+ ,304252
+ ,353
+ ,20850
+ ,98
+ ,-4
+ ,776
+ ,9
+ ,8
+ ,64230
+ ,333505
+ ,113
+ ,14834
+ ,172
+ ,-4
+ ,611
+ ,8
+ ,6
+ ,59190
+ ,296919
+ ,371
+ ,12115
+ ,159
+ ,-4
+ ,1526
+ ,25
+ ,21
+ ,69351
+ ,278990
+ ,391
+ ,3160
+ ,52
+ ,-4
+ ,592
+ ,7
+ ,7
+ ,64270
+ ,276898
+ ,392
+ ,10985
+ ,130
+ ,-4
+ ,1182
+ ,13
+ ,11
+ ,70694
+ ,327007
+ ,124
+ ,9770
+ ,107
+ ,-5
+ ,621
+ ,16
+ ,11
+ ,68005
+ ,317046
+ ,156
+ ,7315
+ ,188
+ ,-5
+ ,989
+ ,11
+ ,12
+ ,58930
+ ,304555
+ ,350
+ ,9505
+ ,106
+ ,-5
+ ,438
+ ,11
+ ,9
+ ,58320
+ ,298096
+ ,366
+ ,8918
+ ,224
+ ,-5
+ ,726
+ ,3
+ ,3
+ ,69980
+ ,231861
+ ,412
+ ,10620
+ ,44
+ ,-5
+ ,1303
+ ,61
+ ,57
+ ,69863
+ ,309422
+ ,336
+ ,1794
+ ,84
+ ,-5
+ ,7419
+ ,29
+ ,21
+ ,63255
+ ,286963
+ ,387
+ ,2999
+ ,12
+ ,-5
+ ,1164
+ ,17
+ ,15
+ ,57320
+ ,269753
+ ,395
+ ,4103
+ ,60
+ ,-6
+ ,3310
+ ,33
+ ,32
+ ,75230
+ ,448243
+ ,44
+ ,7523
+ ,75
+ ,-6
+ ,1920
+ ,15
+ ,11
+ ,79420
+ ,165404
+ ,424
+ ,-2306
+ ,-18
+ ,-6
+ ,965
+ ,3
+ ,2
+ ,73490
+ ,204325
+ ,421
+ ,1442
+ ,4
+ ,-8
+ ,3256
+ ,66
+ ,23
+ ,35250
+ ,407159
+ ,57
+ ,3139
+ ,64
+ ,-8
+ ,1135
+ ,17
+ ,20
+ ,62285
+ ,290476
+ ,384
+ ,5322
+ ,80
+ ,-9
+ ,1270
+ ,26
+ ,24
+ ,69206
+ ,275311
+ ,393
+ ,2897
+ ,59
+ ,-10
+ ,661
+ ,3
+ ,1
+ ,65920
+ ,246541
+ ,405
+ ,15514
+ ,70
+ ,-12
+ ,1013
+ ,2
+ ,1
+ ,69770
+ ,253468
+ ,402
+ ,26734
+ ,53
+ ,-14
+ ,2844
+ ,67
+ ,74
+ ,72683
+ ,240897
+ ,407
+ ,610
+ ,14
+ ,-15
+ ,11528
+ ,70
+ ,68
+ ,-14545
+ ,-83265
+ ,431
+ ,-4047
+ ,-25
+ ,-25
+ ,6526
+ ,26
+ ,20
+ ,55830
+ ,-42143
+ ,430
+ ,-9313
+ ,-37
+ ,-28
+ ,2264
+ ,24
+ ,20
+ ,55174
+ ,272713
+ ,394
+ ,3030
+ ,32
+ ,-31
+ ,5109
+ ,97
+ ,82
+ ,67038
+ ,215362
+ ,417
+ ,158
+ ,3
+ ,-31
+ ,3999
+ ,30
+ ,21
+ ,51252
+ ,42754
+ ,428
+ ,-5242
+ ,-39
+ ,-58
+ ,35624
+ ,223
+ ,244
+ ,157278
+ ,306275
+ ,347
+ ,477
+ ,3
+ ,-70
+ ,9252
+ ,48
+ ,32
+ ,79510
+ ,253537
+ ,401
+ ,1115
+ ,6
+ ,-129
+ ,15236
+ ,90
+ ,86
+ ,77440
+ ,372631
+ ,78
+ ,1918
+ ,11
+ ,-158
+ ,18073
+ ,180
+ ,69
+ ,27284
+ ,-7170
+ ,429
+ ,-1151
+ ,-11
+ ,-366)
+ ,dim=c(9
+ ,431)
+ ,dimnames=list(c('Costs'
+ ,'Trades'
+ ,'Orders'
+ ,'Dividends'
+ ,'Wealth'
+ ,'WRank'
+ ,'Profit/Trades'
+ ,'Profit/Cost'
+ ,'CCScore
')
+ ,1:431))
> y <- array(NA,dim=c(9,431),dimnames=list(c('Costs','Trades','Orders','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 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '9'
> 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 Orders Dividends Wealth WRank Profit/Trades
1 56289 162556 1081 807 213118 6282154 1 5626
2 28328 29790 309 444 81767 4321023 2 13337
3 17936 87550 458 412 153198 4111912 3 8541
4 4145 84738 588 428 -26007 223193 415 39
5 3040 54660 302 315 126942 1491348 6 4276
6 2964 42634 156 168 157214 1629616 5 9164
7 2865 40949 481 263 129352 1398893 9 2493
8 2854 45187 353 267 234817 1926517 4 4891
9 2167 37704 452 228 60448 983660 14 1734
10 1974 16275 109 129 47818 1443586 7 11409
11 1910 25830 115 104 245546 1073089 10 7592
12 1871 12679 110 122 48020 984885 13 7135
13 943 18014 239 393 -1710 1405225 8 5043
14 929 43556 247 190 32648 227132 414 110
15 822 24811 505 280 95350 929118 15 1444
16 819 6575 159 63 151352 1071292 11 5480
17 769 7123 109 102 288170 638830 28 4026
18 745 21950 519 265 114337 856956 17 1266
19 652 37597 248 234 37884 992426 12 3195
20 643 17821 373 277 122844 444477 46 655
21 601 12988 119 73 82340 857217 16 5523
22 446 22330 84 67 79801 711969 20 6095
23 436 13326 102 103 165548 702380 21 4925
24 379 16189 295 290 116384 358589 89 538
25 305 7146 105 83 134028 297978 367 933
26 284 15824 64 56 63838 585715 30 6027
27 247 27664 282 236 74996 657954 24 1624
28 238 11920 182 73 31080 209458 419 52
29 223 8568 37 34 32168 786690 18 15856
30 220 14416 361 139 49857 439798 48 664
31 217 3369 28 26 87161 688779 23 17456
32 199 11819 85 70 106113 574339 31 4404
33 195 6984 45 40 80570 741409 19 12031
34 163 4519 49 42 102129 597793 29 8118
35 154 2220 22 12 301670 644190 26 20190
36 143 18562 155 211 102313 377934 72 1148
37 135 10327 91 74 88577 640273 27 4838
38 132 5336 81 80 112477 697458 22 6141
39 120 2365 79 83 191778 550608 32 4438
40 119 4069 145 131 79804 207393 420 51
41 108 8636 855 203 128294 301607 359 119
42 104 13718 61 56 96448 345783 100 2390
43 101 4525 226 89 93811 501749 37 1335
44 90 6869 105 88 117520 379983 68 1714
45 85 4628 62 39 69159 387475 64 3024
46 74 3689 25 25 101792 377305 74 7092
47 73 4891 217 49 210568 370837 79 787
48 70 7489 322 149 136996 430866 50 717
49 67 4901 84 58 121920 469107 41 3204
50 66 2284 33 41 76403 194493 423 -167
51 66 3160 108 90 108094 530670 33 3062
52 66 4150 150 136 134759 518365 35 2122
53 59 7285 115 97 188873 491303 39 2533
54 58 1134 162 63 146216 527021 34 2019
55 58 4658 158 114 156608 233773 411 214
56 54 2384 97 77 61348 405972 58 2123
57 53 3748 9 6 50350 652925 25 50325
58 49 5371 66 47 87720 446211 45 3730
59 49 1285 107 51 99489 341340 106 1321
60 44 9327 101 85 87419 387699 63 1858
61 39 5565 47 43 94355 493408 38 6243
62 38 1528 38 32 60326 146494 426 -1408
63 38 3122 34 25 94670 414462 54 6308
64 37 7561 87 77 82425 364304 86 1889
65 36 2675 79 54 59017 355178 93 1964
66 35 13253 947 251 90829 357760 90 167
67 34 880 74 15 80791 261216 399 827
68 33 2053 53 44 100423 397144 60 3720
69 32 1424 94 73 131116 374943 77 1861
70 32 4036 63 85 100269 424898 53 3570
71 31 3045 58 49 27330 202055 422 35
72 31 5119 49 38 39039 378525 69 3643
73 30 1431 34 35 106885 310768 330 3258
74 30 554 11 9 79285 325738 125 11431
75 30 1975 35 34 118881 394510 61 5557
76 28 1765 20 20 77623 247060 404 2353
77 26 1012 47 29 114768 368078 82 3576
78 26 810 43 11 74015 236761 409 855
79 26 1280 117 52 69465 312378 326 960
80 25 666 171 13 117869 339836 108 818
81 25 1380 26 29 60982 347385 98 5669
82 24 4677 75 66 90131 426280 52 3017
83 24 876 59 33 138971 352850 96 2591
84 23 814 18 15 39625 301881 358 5660
85 23 514 15 15 102725 377516 73 11834
86 22 5692 72 68 64239 357312 91 2185
87 22 3642 86 100 90262 458343 42 3004
88 22 540 14 13 103960 354228 94 11016
89 20 2099 64 45 106611 308636 339 1697
90 20 567 11 14 103345 386212 65 16928
91 19 2001 52 36 95551 393343 62 3718
92 18 2949 41 40 82903 378509 70 4354
93 17 2253 99 68 63593 452469 43 2550
94 16 6533 75 29 126910 364839 85 2198
95 14 1889 45 43 37527 358649 88 3526
96 14 3055 43 30 60247 376641 75 4108
97 14 272 8 9 112995 429112 51 28639
98 13 1414 198 22 70184 330546 119 659
99 13 2564 22 19 130140 403560 59 9253
100 13 1383 11 9 73221 317892 151 10717
101 13 1261 33 31 76114 307528 343 3258
102 13 975 23 19 90534 235133 410 1528
103 12 3366 80 55 108479 299243 362 1241
104 12 576 18 8 113761 314073 313 6337
105 11 1686 40 28 68696 368186 81 4205
106 10 746 23 29 71561 269661 396 3029
107 10 3192 60 48 59831 125390 427 -1243
108 10 2045 20 16 97890 510834 36 15542
109 10 5702 61 47 101481 321896 138 1998
110 10 1932 36 20 72954 249898 403 1386
111 9 936 30 22 67939 408881 56 6963
112 9 3437 47 33 48022 158492 425 -883
113 9 5131 71 44 86111 292154 379 1298
114 6 2397 14 13 74020 289513 385 6394
115 6 1389 9 6 57530 378049 71 19783
116 6 1503 39 35 56364 343466 104 3679
117 5 402 26 8 84990 332743 115 5106
118 5 2239 21 17 88590 442882 47 11566
119 5 2234 16 11 77200 214215 418 888
120 4 837 69 21 61262 315688 165 1677
121 4 10579 92 92 110309 375195 76 1904
122 3 875 14 12 67000 334280 112 9591
123 2 1585 107 112 93099 355864 92 1457
124 2 1659 29 25 107577 480382 40 9668
125 1 2647 37 17 62920 353058 95 4137
126 1 3294 23 23 75832 217193 416 748
127 7 0 0 0 60720 315380 191 0
128 28 10 60793 314533 308 16362 1212 0
129 0 23 57935 318056 150 4216 280 0
130 60630 0 60720 315380 197 0 34 8
131 55637 314353 309 14294 3340 0 1558 63
132 60720 369448 80 2690 109 0 0 0
133 323 315380 196 0 43 3 20 60887
134 329 37615 2613 0 645 5 4 60720
135 299 22415 174 0 316 9 4 60505
136 145 12779 364 0 115 13 10 60945
137 304 9146 1037 0 5 2 1 60720
138 303 57443 24814 0 897 5 4 60720
139 188 22983 128 0 0 0 0 60720
140 322 0 389 14 8 58990 325506 127
141 1002 0 0 0 0 60720 315380 220
142 36 15 11 56750 298568 364 6571 98
143 460 3 4 60894 315834 162 38611 3215
144 309 15 15 63346 329784 122 8652 282
145 0 11 9 56535 312878 322 10262 365
146 7 0 0 60720 315380 203 0 9
147 2 60835 314987 300 19165 12458 0 271
148 0 60720 325249 129 62625 463 0 14
149 7 61016 315877 161 115877 8041 0 520
150 46 58650 291650 383 9165 176 0 1766
151 5 60438 305959 349 1451 60 0 0
152 297765 60720 315380 208 0 458 11 7
153 315245 368 8888 213 0 20 3 2
154 315380 297 38415 5674 0 0 0 0
155 0 217 0 0 0 0 60720 315380
156 0 98 2 2 61490 315236 298 57618
157 0 405 7 5 60845 336425 111 19489
158 0 0 0 0 60720 315380 213 0
159 60720 0 60720 315380 212 0 0 0
160 218 315380 207 0 0 0 0 60720
161 220 0 483 27 7 60830 306268 348
162 225 0 454 51 24 63261 302187 357
163 2432 0 47 3 1 60720 314882 305
164 757 0 0 0 0 60720 315380 193
165 4655 19 18 45689 382712 66 9616 242
166 0 393 55 60720 341570 105 360 30
167 0 0 0 60720 315380 192 0 0
168 312412 60720 315380 211 0 36 4 3
169 315380 325 28103 3146 0 0 0 0
170 12177 189 0 203 9 9 61938 309596
171 0 541 0 0 0 0 60720 315380
172 0 126 10 8 60951 315547 171 11555
173 0 400 152 113 60720 313267 319 745
174 0 71 1 0 60745 316176 159 116176
175 0 0 0 0 60720 315380 209 0
176 71642 0 60720 315380 190 0 972 34
177 71641 359335 87 4686 164 0 531 10
178 55792 330068 120 13007 245 0 2461 57
179 71873 314289 310 2005 46 0 378 52
180 62555 297413 369 1873 258 0 23 5
181 60370 314806 306 22961 4974 0 638 14
182 64873 333210 114 9515 209 0 2300 29
183 62041 352108 97 5245 66 0 149 5
184 65745 313332 318 22666 761 0 226 5
185 60720 291787 382 18357 406 0 0 0
186 146 315380 184 0 275 4 7 59500
187 179 29686 431 0 0 0 0 60720
188 819 0 141 6 7 61630 315366 295
189 28 0 0 0 0 60720 315380 186
190 0 2 3 60890 315688 166 57844 4148
191 89 0 0 60720 315380 263 0 4980
192 0 113521 409642 55 2304 42 0 0
193 315380 60720 315380 259 0 0 0 0
194 3479 261 0 472 20 19 80045 269587
195 0 148 0 0 0 0 60720 315380
196 0 0 0 0 60720 315380 250 0
197 50804 0 60720 315380 247 0 203 27
198 87390 300962 361 3739 498 0 496 17
199 61656 325479 128 7381 253 0 10 2
200 65688 316155 160 58077 11639 0 63 4
201 60720 318574 148 29643 1891 0 0 0
202 103 315380 228 0 1136 32 26 48522
203 346 4488 126 0 265 31 3 60720
204 0 3450 404 0 0 0 60720 315380
205 32 0 0 0 60720 315380 206 0
206 20 11 57640 330059 121 4064 487 0
207 7 10 61977 288985 386 4449 188 0
208 0 5 62620 304485 351 14926 196 0
209 60831 2 60720 315380 201 0 15 8
210 60646 315688 164 14461 7687 0 397 28
211 60720 317736 152 4205 297 0 0 0
212 135 315380 199 0 1866 29 28 56225
213 373 4218 66 0 288 4 3 60510
214 198 24164 335 0 0 0 0 60720
215 36505 0 3 2 1 60698 315354 296
216 239 0 468 21 20 60720 312161 328
217 5924 0 20 2 1 60805 315576 169
218 413 0 278 26 22 61404 314922 302
219 1874 0 61 14 9 60720 314551 307
220 192 0 0 0 0 60720 315380 215
221 0 4 2 65276 312339 327 28085 586
222 7 0 0 60720 315380 214 0 317
223 9 63915 298700 363 10967 311 0 738
224 0 60720 321376 139 12138 164 0 0
225 303230 60720 315380 216 0 368 17 13
226 315380 355 6072 281 0 0 0 0
227 115487 204 0 2 1 0 60743 315487
228 0 52257 0 0 0 0 60720 315380
229 0 53 6 6 60349 315793 163 19299
230 0 0 0 0 60720 315380 195 0
231 60720 0 60720 315380 194 0 0 0
232 321 315380 200 0 94 3 3 61360
233 237 37629 1200 0 0 0 0 60720
234 4772 0 24 8 7 59818 315637 168
235 53 0 2332 4 2 72680 324385 131
236 0 0 0 0 0 60720 315380 236
237 15 0 0 60720 315380 233 0 131
238 0 61808 308989 338 9908 834 0 0
239 315380 60720 315380 260 0 0 0 0
240 10745 234 0 206 9 9 53110 296702
241 0 469 0 0 0 0 60720 315380
242 0 167 2 1 64245 307322 345 53661
243 0 622 73 38 73007 304376 352 1430
244 0 2328 94 57 82732 253588 400 570
245 8 0 0 0 60720 315380 238 0
246 35 7 54820 309560 335 13695 300 0
247 0 26 47705 298466 365 2813 271 0
248 60720 0 60720 315380 239 0 0 0
249 262 315380 235 0 0 0 0 60720
250 226 0 0 0 0 60720 315380 283
251 307 12 13 72835 343929 102 11994 636
252 0 15 10 58856 331955 116 8797 430
253 0 0 0 60720 315380 288 0 0
254 315380 60720 315380 281 0 0 0 0
255 16471 282 0 188 11 9 77655 381180
256 0 964 0 0 0 0 60720 315380
257 0 138 6 26 69817 331420 117 21903
258 0 0 0 0 60720 315380 286 0
259 60720 0 60720 315380 285 0 0 0
260 332 315380 265 0 125 12 19 60798
261 284 9183 883 0 0 0 0 60720
262 425 0 282 30 12 62452 320016 141
263 359 0 335 33 23 64175 320398 140
264 1324 0 0 0 0 60720 315380 232
265 176 117 29 67440 291841 381 785 69
266 0 28 8 68136 310670 331 3952 629
267 0 0 0 60720 315380 256 0 0
268 313491 60720 315380 255 0 249 72 26
269 315380 317 1576 455 0 0 0 0
270 10102 242 0 333 13 9 70811 331323
271 0 394 0 0 0 0 60720 315380
272 0 601 6 5 60720 319210 144 19868
273 0 30 4 3 62045 318098 149 29525
274 62 0 0 0 60720 315380 254 0
275 0 13 54323 292754 378 1496 373 0
276 62841 0 60720 315380 258 0 165 24
277 81125 325176 130 5216 759 0 453 21
278 60720 365959 84 7903 366 0 0 0
279 356 315380 257 0 53 14 10 59506
280 107 7315 1932 0 382 21 9 59365
281 229 6713 369 0 0 0 0 60720
282 0 0 0 0 0 60720 315380 230
283 9 0 0 60720 315380 224 0 0
284 313164 60720 315380 223 0 30 4 4
285 301164 320 28291 3808 0 290 2 1
286 315380 360 50582 348 0 0 0 1
287 0 225 0 0 0 0 60720 315380
288 0 366 53 14 61808 344425 101 2725
289 0 2 9 12 60735 315394 174 12822
290 13 0 0 0 60720 315380 222 0
291 22 19 64016 316647 157 8973 558 0
292 0 17 54683 309836 333 4993 286 0
293 60720 0 60720 315380 227 0 0 0
294 99 315380 244 0 365 83 32 87192
295 226 1766 402 0 0 0 0 60720
296 2500 0 49 8 14 64107 322031 137
297 37920 0 3 4 8 60761 315656 167
298 1052 0 133 14 4 65990 339445 109
299 3603 0 32 1 0 59988 314964 301
300 264 0 368 17 20 61167 297141 370
301 84213 0 1 6 5 60719 315372 293
302 0 0 0 0 0 60720 315380 274
303 0 0 0 60720 315380 246 0 0
304 315380 60720 315380 177 0 0 0 0
305 0 185 0 0 0 0 60720 315380
306 2 0 0 0 60720 315380 273 0
307 0 1 60722 312502 324 56251 5154 0
308 60720 0 60720 315380 252 0 0 0
309 253 315380 249 0 0 0 0 60720
310 0 0 0 0 0 60720 315380 277
311 0 0 0 60720 315380 278 0 0
312 315380 60720 315380 275 0 0 0 0
313 22746 276 0 96 5 4 60379 313729
314 57694 1188 0 1 2 1 60727 315388
315 23074 117743 0 314 5 4 60720 315371
316 1233 367 0 844 78 20 60925 296139
317 0 114 0 0 0 0 60720 315380
318 0 26 1 1 60896 313880 315 113880
319 0 125 13 10 59734 317698 153 9054
320 0 304 15 12 62969 295580 376 6372
321 0 0 0 0 60720 315380 187 0
322 60720 0 60720 315380 183 0 0 0
323 341 315380 182 0 621 48 13 60720
324 181 2255 174 0 0 0 0 60720
325 872 0 119 6 3 59118 303677 354
326 0 0 0 0 0 60720 315380 178
327 10 0 0 60720 315380 176 0 1595
328 3 60720 319369 143 7022 75 0 312
329 7 58598 318690 147 8478 381 0 60
330 10 61124 314049 314 11405 1916 0 587
331 1 59595 325699 126 10475 214 0 135
332 0 62065 314210 311 57105 847 0 0
333 315380 60720 315380 272 0 0 0 0
334 2353 269 0 514 52 15 78780 322378
335 0 238 0 0 0 0 60720 315380
336 0 0 0 0 60720 315380 251 0
337 60722 0 60720 315380 270 0 1 4
338 60720 315398 173 28849 172235 0 0 0
339 289 315380 268 0 0 0 0 60720
340 61 0 1763 24 28 61600 308336 340
341 647 0 180 11 9 59635 316386 158
342 0 0 0 0 0 60720 315380 287
343 0 0 0 60720 315380 248 0 0
344 315380 60720 315380 291 0 0 0 0
345 5503 292 0 218 21 7 60720 315553
346 0 531 0 0 0 0 60720 315380
347 0 448 40 7 59781 323361 132 3084
348 0 227 9 7 76644 336639 110 15182
349 0 174 1 3 64820 307424 344 107424
350 0 0 0 0 60720 315380 266 0
351 56178 0 60720 315380 280 0 121 24
352 60436 295370 377 3974 788 0 607 11
353 60720 322340 134 11122 201 0 2212 14
354 60720 319864 142 8562 54 0 0 0
355 264 315380 271 0 0 0 0 60720
356 221 0 530 60 18 73433 317291 155
357 141 0 571 80 14 41477 280398 390
358 78 0 0 0 0 60720 315380 267
359 2489 16 12 62700 317330 154 7333 1511
360 131 40 29 67804 238125 408 953 15
361 923 6 3 59661 327071 123 21179 969
362 72 8 6 58620 309038 337 13630 118
363 572 3 3 60398 314210 312 38070 1591
364 397 16 8 58580 307930 342 6746 189
365 450 10 10 62710 322327 136 12233 308
366 622 8 6 59325 292136 380 11517 205
367 694 7 8 60950 263276 398 9039 102
368 3425 8 6 68060 367655 83 20957 241
369 562 12 9 83620 283910 388 6993 25
370 4917 13 8 58456 283587 389 6430 149
371 1442 42 26 52811 243650 406 1039 9
372 529 118 239 121173 438493 49 2021 165
373 2126 9 7 63870 296261 374 10696 182
374 1061 138 41 21001 230621 413 222 14
375 776 5 3 70415 304252 353 20850 98
376 611 9 8 64230 333505 113 14834 172
377 1526 8 6 59190 296919 371 12115 159
378 592 25 21 69351 278990 391 3160 52
379 1182 7 7 64270 276898 392 10985 130
380 621 13 11 70694 327007 124 9770 107
381 989 16 11 68005 317046 156 7315 188
382 438 11 12 58930 304555 350 9505 106
383 726 11 9 58320 298096 366 8918 224
384 1303 3 3 69980 231861 412 10620 44
385 7419 61 57 69863 309422 336 1794 84
386 1164 29 21 63255 286963 387 2999 12
387 3310 17 15 57320 269753 395 4103 60
388 1920 33 32 75230 448243 44 7523 75
389 965 15 11 79420 165404 424 -2306 -18
390 3256 3 2 73490 204325 421 1442 4
391 1135 66 23 35250 407159 57 3139 64
392 1270 17 20 62285 290476 384 5322 80
393 661 26 24 69206 275311 393 2897 59
394 1013 3 1 65920 246541 405 15514 70
395 2844 2 1 69770 253468 402 26734 53
396 11528 67 74 72683 240897 407 610 14
397 6526 70 68 -14545 -83265 431 -4047 -25
398 2264 26 20 55830 -42143 430 -9313 -37
399 5109 24 20 55174 272713 394 3030 32
400 3999 97 82 67038 215362 417 158 3
401 35624 30 21 51252 42754 428 -5242 -39
402 9252 223 244 157278 306275 347 477 3
403 15236 48 32 79510 253537 401 1115 6
404 18073 90 86 77440 372631 78 1918 11
405 162556 180 69 27284 -7170 429 -1151 -11
406 29790 1081 807 213118 6282154 1 5626 37
407 87550 309 444 81767 4321023 2 13337 138
408 84738 458 412 153198 4111912 3 8541 45
409 54660 588 428 -26007 223193 415 39 0
410 42634 302 315 126942 1491348 6 4276 24
411 40949 156 168 157214 1629616 5 9164 34
412 45187 481 263 129352 1398893 9 2493 29
413 37704 353 267 234817 1926517 4 4891 38
414 16275 452 228 60448 983660 14 1734 21
415 25830 109 129 47818 1443586 7 11409 76
416 12679 115 104 245546 1073089 10 7592 34
417 18014 110 122 48020 984885 13 7135 62
418 43556 239 393 -1710 1405225 8 5043 67
419 24811 247 190 32648 227132 414 110 1
420 6575 505 280 95350 929118 15 1444 29
421 7123 159 63 151352 1071292 11 5480 133
422 21950 109 102 288170 638830 28 4026 62
423 37597 519 265 114337 856956 17 1266 30
424 17821 248 234 37884 992426 12 3195 21
425 12988 373 277 122844 444477 46 655 14
426 22330 119 73 82340 857217 16 5523 51
427 13326 84 67 79801 711969 20 6095 23
428 16189 102 103 165548 702380 21 4925 38
429 7146 295 290 116384 358589 89 538 10
430 15824 105 83 134028 297978 367 933 14
431 27664 64 56 63838 585715 30 6027 24
Profit/Cost
1 37
2 138
3 45
4 0
5 24
6 34
7 29
8 38
9 21
10 76
11 34
12 62
13 67
14 1
15 29
16 133
17 62
18 30
19 21
20 14
21 51
22 23
23 38
24 10
25 14
26 24
27 17
28 1
29 68
30 17
31 145
32 32
33 78
34 88
35 200
36 10
37 43
38 93
39 148
40 2
41 12
42 11
43 67
44 26
45 41
46 48
47 35
48 31
49 55
50 -2
51 105
52 77
53 40
54 288
55 7
56 86
57 121
58 46
59 110
60 20
61 53
62 -35
63 69
64 22
65 58
66 12
67 70
68 96
69 123
70 56
71 1
72 35
73 77
74 227
75 98
76 27
77 166
78 45
79 88
80 210
81 107
82 48
83 174
84 125
85 346
86 28
87 71
88 286
89 52
90 329
91 97
92 61
93 112
94 25
95 84
96 58
97 844
98 92
99 79
100 85
101 85
102 36
103 29
104 198
105 100
106 93
107 -23
108 152
109 21
110 26
111 223
112 -12
113 18
114 37
115 128
116 95
117 330
118 108
119 6
120 138
121 17
122 153
123 98
124 169
125 58
126 5
127 94
128 422
129 0
130 7
131 25
132 0
133 312846
134 312075
135 315009
136 318903
137 314887
138 314913
139 315380
140 8965
141 0
142 0
143 0
144 0
145 0
146 6
147 2
148 1
149 10
150 73
151 0
152 58625
153 60938
154 60720
155 219
156 1179
157 337
158 0
159 0
160 315380
161 3936
162 2004
163 38294
164 0
165 0
166 0
167 0
168 61564
169 60720
170 334
171 205
172 921
173 283
174 1646
175 0
176 19
177 11
178 25
179 16
180 5
181 11
182 23
183 6
184 5
185 0
186 318745
187 315380
188 19228
189 0
190 0
191 91
192 0
193 60720
194 397
195 243
196 0
197 12
198 12
199 5
200 2
201 0
202 343613
203 306948
204 210
205 267
206 474
207 534
208 0
209 6
210 7
211 2
212 322331
213 296656
214 315380
215 57677
216 5341
217 57788
218 4420
219 8182
220 0
221 0
222 9
223 10
224 0
225 61686
226 60720
227 172
228 202
229 2197
230 0
231 0
232 312887
233 315380
234 14455
235 31096
236 0
237 11
238 0
239 60720
240 372
241 240
242 642
243 168
244 23
245 365
246 364
247 0
248 0
249 315380
250 0
251 0
252 0
253 0
254 60720
255 67
256 221
257 954
258 0
259 0
260 310201
261 315380
262 4001
263 3648
264 0
265 0
266 0
267 0
268 56726
269 60720
270 118
271 241
272 198
273 3912
274 249
275 0
276 12
277 19
278 0
279 302409
280 340968
281 315380
282 0
283 0
284 60798
285 58790
286 60720
287 231
288 394
289 63057
290 209
291 384
292 0
293 0
294 346611
295 315380
296 15254
297 28914
298 9960
299 114964
300 5714
301 19229
302 0
303 0
304 60720
305 245
306 22
307 0
308 0
309 315380
310 0
311 0
312 60720
313 316
314 175
315 294
316 375
317 180
318 4328
319 944
320 314
321 0
322 0
323 308256
324 315380
325 17280
326 0
327 17
328 14
329 10
330 12
331 2
332 0
333 60720
334 133
335 279
336 0
337 4
338 0
339 315380
340 4514
341 10581
342 0
343 0
344 60720
345 170
346 290
347 275
348 603
349 618
350 0
351 11
352 7
353 10
354 0
355 315380
356 1955
357 1005
358 0
359 0
360 0
361 0
362 0
363 0
364 -1
365 -1
366 -1
367 -2
368 -2
369 -3
370 -3
371 -3
372 -3
373 -3
374 -3
375 -4
376 -4
377 -4
378 -4
379 -4
380 -5
381 -5
382 -5
383 -5
384 -5
385 -5
386 -5
387 -6
388 -6
389 -6
390 -8
391 -8
392 -9
393 -10
394 -12
395 -14
396 -15
397 -25
398 -28
399 -31
400 -31
401 -58
402 -70
403 -129
404 -158
405 -366
406 56289
407 28328
408 17936
409 4145
410 3040
411 2964
412 2865
413 2854
414 2167
415 1974
416 1910
417 1871
418 943
419 929
420 822
421 819
422 769
423 745
424 652
425 643
426 601
427 446
428 436
429 379
430 305
431 284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Trades Orders
1.264e+04 8.168e-02 4.097e-01 -1.466e-02
Dividends Wealth WRank `Profit/Trades`
6.958e-03 -6.208e-03 -2.200e-02 -2.514e-02
`Profit/Cost`
5.197e-03
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-189742 -12391 -10675 -3356 301756
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.264e+04 4.777e+03 2.647 0.00843 **
Costs 8.168e-02 3.405e-02 2.399 0.01689 *
Trades 4.097e-01 3.750e-02 10.923 < 2e-16 ***
Orders -1.466e-02 3.840e-02 -0.382 0.70286
Dividends 6.958e-03 6.158e-03 1.130 0.25911
Wealth -6.208e-03 6.098e-03 -1.018 0.30922
WRank -2.200e-02 3.292e-02 -0.668 0.50422
`Profit/Trades` -2.514e-02 3.867e-02 -0.650 0.51598
`Profit/Cost` 5.197e-03 3.762e-02 0.138 0.89020
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 57400 on 422 degrees of freedom
Multiple R-squared: 0.2564, Adjusted R-squared: 0.2423
F-statistic: 18.19 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,] 1.163617e-03 2.327234e-03 9.988364e-01
[2,] 9.082369e-05 1.816474e-04 9.999092e-01
[3,] 1.939320e-05 3.878640e-05 9.999806e-01
[4,] 1.589270e-06 3.178540e-06 9.999984e-01
[5,] 1.500579e-07 3.001157e-07 9.999998e-01
[6,] 1.711452e-08 3.422905e-08 1.000000e+00
[7,] 1.277687e-09 2.555373e-09 1.000000e+00
[8,] 1.219229e-10 2.438458e-10 1.000000e+00
[9,] 1.197016e-11 2.394032e-11 1.000000e+00
[10,] 1.277401e-12 2.554801e-12 1.000000e+00
[11,] 1.735983e-13 3.471966e-13 1.000000e+00
[12,] 1.646775e-14 3.293551e-14 1.000000e+00
[13,] 1.430237e-15 2.860474e-15 1.000000e+00
[14,] 9.622610e-17 1.924522e-16 1.000000e+00
[15,] 1.192645e-17 2.385290e-17 1.000000e+00
[16,] 8.245700e-19 1.649140e-18 1.000000e+00
[17,] 5.369125e-20 1.073825e-19 1.000000e+00
[18,] 4.377187e-21 8.754375e-21 1.000000e+00
[19,] 3.101932e-22 6.203863e-22 1.000000e+00
[20,] 2.031266e-23 4.062532e-23 1.000000e+00
[21,] 1.854242e-24 3.708484e-24 1.000000e+00
[22,] 1.197557e-25 2.395114e-25 1.000000e+00
[23,] 8.853133e-27 1.770627e-26 1.000000e+00
[24,] 5.163738e-28 1.032748e-27 1.000000e+00
[25,] 5.120832e-29 1.024166e-28 1.000000e+00
[26,] 3.508592e-30 7.017185e-30 1.000000e+00
[27,] 2.101781e-31 4.203562e-31 1.000000e+00
[28,] 1.207007e-32 2.414014e-32 1.000000e+00
[29,] 6.560432e-34 1.312086e-33 1.000000e+00
[30,] 5.379738e-35 1.075948e-34 1.000000e+00
[31,] 4.449088e-36 8.898175e-36 1.000000e+00
[32,] 2.570921e-37 5.141842e-37 1.000000e+00
[33,] 1.950831e-38 3.901662e-38 1.000000e+00
[34,] 1.398825e-39 2.797650e-39 1.000000e+00
[35,] 1.133503e-40 2.267005e-40 1.000000e+00
[36,] 6.752378e-42 1.350476e-41 1.000000e+00
[37,] 3.915163e-43 7.830325e-43 1.000000e+00
[38,] 2.239204e-44 4.478407e-44 1.000000e+00
[39,] 1.126096e-45 2.252193e-45 1.000000e+00
[40,] 5.587147e-47 1.117429e-46 1.000000e+00
[41,] 2.799945e-48 5.599890e-48 1.000000e+00
[42,] 1.513476e-49 3.026953e-49 1.000000e+00
[43,] 1.086207e-50 2.172414e-50 1.000000e+00
[44,] 5.343272e-52 1.068654e-51 1.000000e+00
[45,] 2.904291e-53 5.808581e-53 1.000000e+00
[46,] 2.639192e-54 5.278383e-54 1.000000e+00
[47,] 1.565807e-55 3.131614e-55 1.000000e+00
[48,] 8.092220e-57 1.618444e-56 1.000000e+00
[49,] 4.903822e-58 9.807644e-58 1.000000e+00
[50,] 2.633995e-59 5.267990e-59 1.000000e+00
[51,] 1.218633e-60 2.437265e-60 1.000000e+00
[52,] 6.801830e-62 1.360366e-61 1.000000e+00
[53,] 3.829645e-63 7.659289e-63 1.000000e+00
[54,] 2.022889e-64 4.045778e-64 1.000000e+00
[55,] 8.911845e-66 1.782369e-65 1.000000e+00
[56,] 4.297606e-67 8.595212e-67 1.000000e+00
[57,] 2.159063e-68 4.318125e-68 1.000000e+00
[58,] 1.008703e-69 2.017406e-69 1.000000e+00
[59,] 5.228103e-71 1.045621e-70 1.000000e+00
[60,] 2.244866e-72 4.489732e-72 1.000000e+00
[61,] 1.183822e-73 2.367645e-73 1.000000e+00
[62,] 5.068321e-75 1.013664e-74 1.000000e+00
[63,] 2.116373e-76 4.232745e-76 1.000000e+00
[64,] 1.006387e-77 2.012774e-77 1.000000e+00
[65,] 4.192071e-79 8.384142e-79 1.000000e+00
[66,] 1.731731e-80 3.463463e-80 1.000000e+00
[67,] 7.278347e-82 1.455669e-81 1.000000e+00
[68,] 3.083519e-83 6.167038e-83 1.000000e+00
[69,] 1.226812e-84 2.453625e-84 1.000000e+00
[70,] 5.450037e-86 1.090007e-85 1.000000e+00
[71,] 2.519812e-87 5.039625e-87 1.000000e+00
[72,] 9.860129e-89 1.972026e-88 1.000000e+00
[73,] 4.197690e-90 8.395381e-90 1.000000e+00
[74,] 1.637472e-91 3.274945e-91 1.000000e+00
[75,] 7.717221e-93 1.543444e-92 1.000000e+00
[76,] 3.278987e-94 6.557974e-94 1.000000e+00
[77,] 1.223971e-95 2.447942e-95 1.000000e+00
[78,] 4.616864e-97 9.233728e-97 1.000000e+00
[79,] 1.699178e-98 3.398357e-98 1.000000e+00
[80,] 7.003710e-100 1.400742e-99 1.000000e+00
[81,] 3.117965e-101 6.235930e-101 1.000000e+00
[82,] 1.170106e-102 2.340213e-102 1.000000e+00
[83,] 5.029846e-104 1.005969e-103 1.000000e+00
[84,] 2.020616e-105 4.041233e-105 1.000000e+00
[85,] 8.354111e-107 1.670822e-106 1.000000e+00
[86,] 5.033559e-108 1.006712e-107 1.000000e+00
[87,] 1.803197e-109 3.606394e-109 1.000000e+00
[88,] 7.796235e-111 1.559247e-110 1.000000e+00
[89,] 3.172654e-112 6.345309e-112 1.000000e+00
[90,] 1.147217e-113 2.294434e-113 1.000000e+00
[91,] 4.069072e-115 8.138144e-115 1.000000e+00
[92,] 1.430383e-116 2.860765e-116 1.000000e+00
[93,] 5.162269e-118 1.032454e-117 1.000000e+00
[94,] 1.928849e-119 3.857698e-119 1.000000e+00
[95,] 6.859690e-121 1.371938e-120 1.000000e+00
[96,] 2.365039e-122 4.730078e-122 1.000000e+00
[97,] 8.496377e-124 1.699275e-123 1.000000e+00
[98,] 3.184745e-125 6.369490e-125 1.000000e+00
[99,] 1.079878e-126 2.159755e-126 1.000000e+00
[100,] 3.595989e-128 7.191978e-128 1.000000e+00
[101,] 1.184619e-129 2.369238e-129 1.000000e+00
[102,] 3.984958e-131 7.969917e-131 1.000000e+00
[103,] 1.296101e-132 2.592202e-132 1.000000e+00
[104,] 5.662257e-134 1.132451e-133 1.000000e+00
[105,] 1.984246e-135 3.968492e-135 1.000000e+00
[106,] 6.275883e-137 1.255177e-136 1.000000e+00
[107,] 2.232396e-138 4.464791e-138 1.000000e+00
[108,] 7.076081e-140 1.415216e-139 1.000000e+00
[109,] 2.204419e-141 4.408838e-141 1.000000e+00
[110,] 8.007902e-143 1.601580e-142 1.000000e+00
[111,] 2.720692e-144 5.441383e-144 1.000000e+00
[112,] 9.706163e-146 1.941233e-145 1.000000e+00
[113,] 3.064557e-147 6.129114e-147 1.000000e+00
[114,] 1.031702e-148 2.063403e-148 1.000000e+00
[115,] 3.153502e-150 6.307003e-150 1.000000e+00
[116,] 9.446220e-152 1.889244e-151 1.000000e+00
[117,] 6.164047e-153 1.232809e-152 1.000000e+00
[118,] 1.041407e-153 2.082815e-153 1.000000e+00
[119,] 3.817660e-113 7.635320e-113 1.000000e+00
[120,] 4.729266e-107 9.458532e-107 1.000000e+00
[121,] 4.383530e-108 8.767059e-108 1.000000e+00
[122,] 3.638944e-109 7.277887e-109 1.000000e+00
[123,] 5.860192e-105 1.172038e-104 1.000000e+00
[124,] 1.052371e-102 2.104742e-102 1.000000e+00
[125,] 1.775515e-102 3.551029e-102 1.000000e+00
[126,] 6.429888e-103 1.285978e-102 1.000000e+00
[127,] 1.046522e-103 2.093044e-103 1.000000e+00
[128,] 1.767786e-104 3.535571e-104 1.000000e+00
[129,] 1.583471e-105 3.166943e-105 1.000000e+00
[130,] 1.410100e-106 2.820199e-106 1.000000e+00
[131,] 1.373457e-107 2.746914e-107 1.000000e+00
[132,] 1.302901e-108 2.605802e-108 1.000000e+00
[133,] 1.223322e-109 2.446643e-109 1.000000e+00
[134,] 1.140925e-110 2.281850e-110 1.000000e+00
[135,] 1.056930e-111 2.113859e-111 1.000000e+00
[136,] 3.099620e-112 6.199239e-112 1.000000e+00
[137,] 1.110345e-112 2.220689e-112 1.000000e+00
[138,] 4.352453e-113 8.704906e-113 1.000000e+00
[139,] 1.724157e-113 3.448315e-113 1.000000e+00
[140,] 9.328383e-114 1.865677e-113 1.000000e+00
[141,] 2.668244e-34 5.336488e-34 1.000000e+00
[142,] 1.045326e-15 2.090652e-15 1.000000e+00
[143,] 2.981441e-08 5.962882e-08 1.000000e+00
[144,] 4.914200e-07 9.828399e-07 9.999995e-01
[145,] 3.325959e-07 6.651918e-07 9.999997e-01
[146,] 2.247343e-07 4.494685e-07 9.999998e-01
[147,] 1.517831e-07 3.035663e-07 9.999998e-01
[148,] 1.165051e-07 2.330103e-07 9.999999e-01
[149,] 9.030550e-08 1.806110e-07 9.999999e-01
[150,] 6.009220e-08 1.201844e-07 9.999999e-01
[151,] 3.973670e-08 7.947341e-08 1.000000e+00
[152,] 2.607965e-08 5.215930e-08 1.000000e+00
[153,] 1.704286e-08 3.408572e-08 1.000000e+00
[154,] 1.265428e-08 2.530855e-08 1.000000e+00
[155,] 8.758592e-09 1.751718e-08 1.000000e+00
[156,] 5.961044e-09 1.192209e-08 1.000000e+00
[157,] 3.727863e-06 7.455726e-06 9.999963e-01
[158,] 6.111401e-03 1.222280e-02 9.938886e-01
[159,] 5.545514e-03 1.109103e-02 9.944545e-01
[160,] 4.582230e-03 9.164460e-03 9.954178e-01
[161,] 3.717059e-03 7.434119e-03 9.962829e-01
[162,] 3.005096e-03 6.010192e-03 9.969949e-01
[163,] 2.406270e-03 4.812541e-03 9.975937e-01
[164,] 1.927960e-03 3.855920e-03 9.980720e-01
[165,] 1.678964e-03 3.357927e-03 9.983210e-01
[166,] 1.501230e-03 3.002461e-03 9.984988e-01
[167,] 1.218754e-03 2.437507e-03 9.987812e-01
[168,] 1.041221e-03 2.082442e-03 9.989588e-01
[169,] 8.468727e-04 1.693745e-03 9.991531e-01
[170,] 6.776524e-04 1.355305e-03 9.993223e-01
[171,] 5.437649e-04 1.087530e-03 9.994562e-01
[172,] 4.289855e-04 8.579710e-04 9.995710e-01
[173,] 3.437309e-04 6.874618e-04 9.996563e-01
[174,] 2.711307e-04 5.422615e-04 9.997289e-01
[175,] 2.436245e-04 4.872490e-04 9.997564e-01
[176,] 1.920169e-04 3.840337e-04 9.998080e-01
[177,] 1.460871e-04 2.921743e-04 9.998539e-01
[178,] 1.106882e-04 2.213765e-04 9.998893e-01
[179,] 8.726897e-05 1.745379e-04 9.999127e-01
[180,] 6.822263e-05 1.364453e-04 9.999318e-01
[181,] 4.739073e-04 9.478146e-04 9.995261e-01
[182,] 6.791975e-03 1.358395e-02 9.932080e-01
[183,] 5.625644e-03 1.125129e-02 9.943744e-01
[184,] 4.627578e-03 9.255156e-03 9.953724e-01
[185,] 3.779370e-03 7.558740e-03 9.962206e-01
[186,] 3.109619e-03 6.219237e-03 9.968904e-01
[187,] 2.874818e-03 5.749636e-03 9.971252e-01
[188,] 2.365217e-03 4.730435e-03 9.976348e-01
[189,] 1.985509e-03 3.971017e-03 9.980145e-01
[190,] 1.635163e-03 3.270327e-03 9.983648e-01
[191,] 1.461758e-03 2.923517e-03 9.985382e-01
[192,] 1.187675e-03 2.375351e-03 9.988123e-01
[193,] 9.448220e-04 1.889644e-03 9.990552e-01
[194,] 7.482221e-04 1.496444e-03 9.992518e-01
[195,] 6.594012e-04 1.318802e-03 9.993406e-01
[196,] 5.722416e-04 1.144483e-03 9.994278e-01
[197,] 4.877084e-04 9.754168e-04 9.995123e-01
[198,] 4.064155e-04 8.128309e-04 9.995936e-01
[199,] 3.279971e-04 6.559943e-04 9.996720e-01
[200,] 2.647203e-04 5.294405e-04 9.997353e-01
[201,] 2.268120e-04 4.536240e-04 9.997732e-01
[202,] 1.798196e-04 3.596393e-04 9.998202e-01
[203,] 1.431079e-04 2.862157e-04 9.998569e-01
[204,] 1.189685e-04 2.379371e-04 9.998810e-01
[205,] 9.050916e-05 1.810183e-04 9.999095e-01
[206,] 6.830378e-05 1.366076e-04 9.999317e-01
[207,] 5.145599e-05 1.029120e-04 9.999485e-01
[208,] 3.852855e-05 7.705711e-05 9.999615e-01
[209,] 2.875453e-05 5.750906e-05 9.999712e-01
[210,] 2.221370e-05 4.442740e-05 9.999778e-01
[211,] 1.710546e-05 3.421092e-05 9.999829e-01
[212,] 6.312628e-05 1.262526e-04 9.999369e-01
[213,] 2.779532e-04 5.559064e-04 9.997220e-01
[214,] 2.357386e-03 4.714773e-03 9.976426e-01
[215,] 1.480488e-01 2.960976e-01 8.519512e-01
[216,] 2.027941e-01 4.055881e-01 7.972059e-01
[217,] 1.839977e-01 3.679954e-01 8.160023e-01
[218,] 1.665364e-01 3.330728e-01 8.334636e-01
[219,] 1.502401e-01 3.004802e-01 8.497599e-01
[220,] 1.382768e-01 2.765536e-01 8.617232e-01
[221,] 1.283752e-01 2.567503e-01 8.716248e-01
[222,] 1.162598e-01 2.325196e-01 8.837402e-01
[223,] 1.031607e-01 2.063214e-01 8.968393e-01
[224,] 9.126296e-02 1.825259e-01 9.087370e-01
[225,] 8.042251e-02 1.608450e-01 9.195775e-01
[226,] 7.179800e-02 1.435960e-01 9.282020e-01
[227,] 1.534022e-01 3.068044e-01 8.465978e-01
[228,] 3.222566e-01 6.445132e-01 6.777434e-01
[229,] 2.977010e-01 5.954019e-01 7.022990e-01
[230,] 2.740340e-01 5.480680e-01 7.259660e-01
[231,] 2.515851e-01 5.031703e-01 7.484149e-01
[232,] 2.302473e-01 4.604946e-01 7.697527e-01
[233,] 2.099881e-01 4.199762e-01 7.900119e-01
[234,] 1.908753e-01 3.817506e-01 8.091247e-01
[235,] 1.775741e-01 3.551482e-01 8.224259e-01
[236,] 1.640475e-01 3.280949e-01 8.359525e-01
[237,] 1.515819e-01 3.031637e-01 8.484181e-01
[238,] 1.407286e-01 2.814572e-01 8.592714e-01
[239,] 1.257219e-01 2.514439e-01 8.742781e-01
[240,] 1.138301e-01 2.276603e-01 8.861699e-01
[241,] 1.026630e-01 2.053259e-01 8.973370e-01
[242,] 9.213250e-02 1.842650e-01 9.078675e-01
[243,] 2.064413e-01 4.128826e-01 7.935587e-01
[244,] 1.883890e-01 3.767781e-01 8.116110e-01
[245,] 1.699956e-01 3.399913e-01 8.300044e-01
[246,] 1.529919e-01 3.059837e-01 8.470081e-01
[247,] 1.372492e-01 2.744984e-01 8.627508e-01
[248,] 1.263743e-01 2.527486e-01 8.736257e-01
[249,] 1.166655e-01 2.333310e-01 8.833345e-01
[250,] 1.052642e-01 2.105284e-01 8.947358e-01
[251,] 9.297236e-02 1.859447e-01 9.070276e-01
[252,] 8.177053e-02 1.635411e-01 9.182295e-01
[253,] 7.159630e-02 1.431926e-01 9.284037e-01
[254,] 6.337425e-02 1.267485e-01 9.366257e-01
[255,] 5.597594e-02 1.119519e-01 9.440241e-01
[256,] 4.934057e-02 9.868113e-02 9.506594e-01
[257,] 1.205638e-01 2.411276e-01 8.794362e-01
[258,] 7.266030e-01 5.467941e-01 2.733970e-01
[259,] 7.023280e-01 5.953440e-01 2.976720e-01
[260,] 6.768823e-01 6.462355e-01 3.231177e-01
[261,] 6.510922e-01 6.978156e-01 3.489078e-01
[262,] 6.244711e-01 7.510579e-01 3.755289e-01
[263,] 5.973840e-01 8.052320e-01 4.026160e-01
[264,] 5.772156e-01 8.455689e-01 4.227844e-01
[265,] 5.573575e-01 8.852850e-01 4.426425e-01
[266,] 5.467997e-01 9.064005e-01 4.532003e-01
[267,] 5.255815e-01 9.488370e-01 4.744185e-01
[268,] 5.043967e-01 9.912067e-01 4.956033e-01
[269,] 4.810860e-01 9.621720e-01 5.189140e-01
[270,] 4.593575e-01 9.187150e-01 5.406425e-01
[271,] 4.306237e-01 8.612474e-01 5.693763e-01
[272,] 4.051225e-01 8.102451e-01 5.948775e-01
[273,] 6.056028e-01 7.887944e-01 3.943972e-01
[274,] 9.700192e-01 5.996168e-02 2.998084e-02
[275,] 9.999330e-01 1.340416e-04 6.702082e-05
[276,] 9.999090e-01 1.819447e-04 9.097233e-05
[277,] 9.998778e-01 2.443634e-04 1.221817e-04
[278,] 9.998364e-01 3.272596e-04 1.636298e-04
[279,] 9.997824e-01 4.352535e-04 2.176268e-04
[280,] 9.997466e-01 5.068559e-04 2.534280e-04
[281,] 9.997028e-01 5.944521e-04 2.972261e-04
[282,] 9.996235e-01 7.530450e-04 3.765225e-04
[283,] 9.995398e-01 9.204918e-04 4.602459e-04
[284,] 9.994007e-01 1.198695e-03 5.993476e-04
[285,] 9.992150e-01 1.569994e-03 7.849969e-04
[286,] 9.990628e-01 1.874479e-03 9.372394e-04
[287,] 9.987826e-01 2.434855e-03 1.217428e-03
[288,] 9.984255e-01 3.149084e-03 1.574542e-03
[289,] 9.979805e-01 4.039015e-03 2.019507e-03
[290,] 9.985199e-01 2.960161e-03 1.480081e-03
[291,] 9.980906e-01 3.818767e-03 1.909383e-03
[292,] 9.975885e-01 4.822920e-03 2.411460e-03
[293,] 9.998853e-01 2.293741e-04 1.146871e-04
[294,] 9.998439e-01 3.121566e-04 1.560783e-04
[295,] 9.997892e-01 4.216406e-04 2.108203e-04
[296,] 9.997548e-01 4.904478e-04 2.452239e-04
[297,] 9.996844e-01 6.312531e-04 3.156266e-04
[298,] 9.996175e-01 7.649380e-04 3.824690e-04
[299,] 9.994875e-01 1.025031e-03 5.125155e-04
[300,] 9.993293e-01 1.341321e-03 6.706605e-04
[301,] 9.999951e-01 9.856069e-06 4.928034e-06
[302,] 9.999930e-01 1.398252e-05 6.991259e-06
[303,] 9.999934e-01 1.321884e-05 6.609419e-06
[304,] 9.999905e-01 1.900553e-05 9.502766e-06
[305,] 9.999861e-01 2.770604e-05 1.385302e-05
[306,] 9.999799e-01 4.016022e-05 2.008011e-05
[307,] 9.999711e-01 5.770211e-05 2.885105e-05
[308,] 9.999590e-01 8.204093e-05 4.102047e-05
[309,] 9.999420e-01 1.159979e-04 5.799894e-05
[310,] 9.999185e-01 1.629998e-04 8.149988e-05
[311,] 9.998959e-01 2.082833e-04 1.041417e-04
[312,] 9.998737e-01 2.525704e-04 1.262852e-04
[313,] 9.998257e-01 3.486545e-04 1.743273e-04
[314,] 9.997579e-01 4.841041e-04 2.420521e-04
[315,] 9.996660e-01 6.679619e-04 3.339810e-04
[316,] 9.995513e-01 8.973689e-04 4.486844e-04
[317,] 9.998666e-01 2.667359e-04 1.333680e-04
[318,] 9.999778e-01 4.430335e-05 2.215168e-05
[319,] 9.999990e-01 2.050720e-06 1.025360e-06
[320,] 1.000000e+00 7.966356e-10 3.983178e-10
[321,] 1.000000e+00 3.827289e-26 1.913644e-26
[322,] 1.000000e+00 4.410911e-26 2.205455e-26
[323,] 1.000000e+00 1.371157e-25 6.855783e-26
[324,] 1.000000e+00 4.205891e-25 2.102945e-25
[325,] 1.000000e+00 1.273465e-24 6.367325e-25
[326,] 1.000000e+00 3.584761e-24 1.792380e-24
[327,] 1.000000e+00 1.021440e-23 5.107202e-24
[328,] 1.000000e+00 2.968128e-23 1.484064e-23
[329,] 1.000000e+00 8.843404e-23 4.421702e-23
[330,] 1.000000e+00 2.565767e-22 1.282884e-22
[331,] 1.000000e+00 7.477848e-22 3.738924e-22
[332,] 1.000000e+00 1.890453e-21 9.452265e-22
[333,] 1.000000e+00 2.456364e-21 1.228182e-21
[334,] 1.000000e+00 7.085765e-21 3.542883e-21
[335,] 1.000000e+00 2.035940e-20 1.017970e-20
[336,] 1.000000e+00 5.789686e-20 2.894843e-20
[337,] 1.000000e+00 1.639129e-19 8.195645e-20
[338,] 1.000000e+00 3.152652e-19 1.576326e-19
[339,] 1.000000e+00 3.476704e-19 1.738352e-19
[340,] 1.000000e+00 1.837036e-19 9.185179e-20
[341,] 1.000000e+00 4.981693e-19 2.490846e-19
[342,] 1.000000e+00 1.423027e-18 7.115137e-19
[343,] 1.000000e+00 1.705259e-18 8.526294e-19
[344,] 1.000000e+00 3.278363e-18 1.639182e-18
[345,] 1.000000e+00 6.628741e-18 3.314370e-18
[346,] 1.000000e+00 1.810101e-17 9.050505e-18
[347,] 1.000000e+00 2.822671e-17 1.411336e-17
[348,] 1.000000e+00 7.907613e-17 3.953806e-17
[349,] 1.000000e+00 1.957106e-16 9.785529e-17
[350,] 1.000000e+00 5.307747e-16 2.653874e-16
[351,] 1.000000e+00 1.421296e-15 7.106478e-16
[352,] 1.000000e+00 2.339586e-15 1.169793e-15
[353,] 1.000000e+00 6.343920e-15 3.171960e-15
[354,] 1.000000e+00 1.555324e-14 7.776620e-15
[355,] 1.000000e+00 4.179485e-14 2.089743e-14
[356,] 1.000000e+00 1.081012e-13 5.405059e-14
[357,] 1.000000e+00 2.387429e-13 1.193714e-13
[358,] 1.000000e+00 5.750546e-13 2.875273e-13
[359,] 1.000000e+00 1.492127e-12 7.460636e-13
[360,] 1.000000e+00 3.207914e-12 1.603957e-12
[361,] 1.000000e+00 8.063717e-12 4.031859e-12
[362,] 1.000000e+00 2.033570e-11 1.016785e-11
[363,] 1.000000e+00 3.688611e-11 1.844306e-11
[364,] 1.000000e+00 9.144119e-11 4.572060e-11
[365,] 1.000000e+00 2.151128e-10 1.075564e-10
[366,] 1.000000e+00 5.246148e-10 2.623074e-10
[367,] 1.000000e+00 1.136414e-09 5.682071e-10
[368,] 1.000000e+00 2.691624e-09 1.345812e-09
[369,] 1.000000e+00 6.330281e-09 3.165140e-09
[370,] 1.000000e+00 1.417540e-08 7.087699e-09
[371,] 1.000000e+00 3.188909e-08 1.594454e-08
[372,] 1.000000e+00 7.139663e-08 3.569831e-08
[373,] 9.999999e-01 1.517410e-07 7.587052e-08
[374,] 9.999998e-01 3.339547e-07 1.669773e-07
[375,] 9.999997e-01 6.142487e-07 3.071244e-07
[376,] 9.999994e-01 1.242503e-06 6.212513e-07
[377,] 9.999987e-01 2.645882e-06 1.322941e-06
[378,] 9.999977e-01 4.591986e-06 2.295993e-06
[379,] 9.999959e-01 8.257482e-06 4.128741e-06
[380,] 9.999917e-01 1.654594e-05 8.272968e-06
[381,] 9.999848e-01 3.038097e-05 1.519049e-05
[382,] 9.999750e-01 4.998039e-05 2.499020e-05
[383,] 9.999564e-01 8.728440e-05 4.364220e-05
[384,] 9.999455e-01 1.089660e-04 5.448298e-05
[385,] 9.999089e-01 1.821976e-04 9.109878e-05
[386,] 9.998542e-01 2.916046e-04 1.458023e-04
[387,] 9.997586e-01 4.828271e-04 2.414135e-04
[388,] 9.997338e-01 5.324842e-04 2.662421e-04
[389,] 9.997661e-01 4.677472e-04 2.338736e-04
[390,] 9.995902e-01 8.195735e-04 4.097867e-04
[391,] 9.994171e-01 1.165773e-03 5.828866e-04
[392,] 9.996705e-01 6.589715e-04 3.294857e-04
[393,] 9.993387e-01 1.322594e-03 6.612970e-04
[394,] 1.000000e+00 2.022534e-08 1.011267e-08
[395,] 1.000000e+00 1.066081e-09 5.330403e-10
[396,] 1.000000e+00 5.246524e-09 2.623262e-09
[397,] 1.000000e+00 1.971987e-08 9.859936e-09
[398,] 1.000000e+00 6.253527e-08 3.126763e-08
[399,] 9.999998e-01 3.162083e-07 1.581042e-07
[400,] 9.999993e-01 1.498040e-06 7.490201e-07
[401,] 9.999989e-01 2.241416e-06 1.120708e-06
[402,] 9.999945e-01 1.094310e-05 5.471551e-06
[403,] 9.999738e-01 5.246902e-05 2.623451e-05
[404,] 9.998994e-01 2.011993e-04 1.005996e-04
[405,] 9.997359e-01 5.282838e-04 2.641419e-04
[406,] 9.987247e-01 2.550574e-03 1.275287e-03
[407,] 9.957717e-01 8.456674e-03 4.228337e-03
[408,] 9.938845e-01 1.223103e-02 6.115513e-03
> postscript(file="/var/www/html/rcomp/tmp/1sb5g1291301535.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/2sb5g1291301535.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/3sb5g1291301535.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/432m11291301535.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/532m11291301535.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 431
Frequency = 1
1 2 3 4 5 6
67597.3278 39723.7974 22636.8972 -14077.8071 -5704.0262 -3969.4508
7 8 9 10 11 12
-5468.7088 -3171.3922 -8007.7062 -3125.1589 -7744.1998 -5891.4039
13 14 15 16 17 18
-4401.0919 -14175.6398 -8909.2994 -6690.3245 -10436.7778 -9343.1248
19 20 21 22 23 24
-9182.1943 -11682.5550 -8262.8788 -10036.2001 -10003.2218 -12271.3588
25 26 27 28 29 30
-12014.9410 -10332.9962 -11163.7273 -12357.7540 -8075.8499 -11345.4600
31 32 33 34 35 36
-8604.2072 -10504.8608 -8691.7598 -9663.8672 -10272.0322 -12412.0269
37 38 39 40 41 42
-9907.2548 -9277.2744 -10552.1168 -12171.3138 -12597.3304 -12146.0802
43 44 45 46 47 48
-10506.8167 -11570.2580 -10959.5187 -11066.6252 -12199.4867 -11574.0036
49 50 51 52 53 54
-10864.9241 -12095.7642 -10258.6656 -10641.4424 -11424.6950 -10438.7596
55 56 57 58 59 60
-12652.7717 -10674.8294 -7931.6590 -10804.9200 -11280.4274 -11554.5445
61 62 63 64 65 66
-10513.2053 -12281.2829 -10800.0118 -11520.8955 -11011.7762 -12479.8449
67 68 69 70 71 72
-11622.4581 -10937.8147 -11301.7111 -10934.6249 -11809.5967 -10878.5941
73 74 75 76 77 78
-11469.1560 -10903.3164 -11025.8128 -11705.7003 -11141.3152 -11715.7240
79 80 81 82 83 84
-11282.1145 -11430.9917 -10864.7099 -10934.9377 -11424.5128 -10945.9050
85 86 87 88 89 90
-10741.8724 -11286.5622 -10658.9196 -10917.5817 -11596.1079 -10569.9313
91 92 93 94 95 96
-10936.9838 -10998.5909 -10418.7580 -11752.1999 -10745.7005 -10872.2244
97 98 99 100 101 102
-10060.0624 -11244.0322 -11015.0061 -11011.1878 -11277.6763 -11841.9552
103 104 105 106 107 108
-11796.1607 -11362.0832 -10871.1308 -11442.5885 -12577.4531 -9927.2119
109 110 111 112 113 114
-11777.8299 -11718.1146 -10481.7946 -12296.7286 -11826.3110 -11387.2132
115 116 117 118 119 120
-10309.3095 -10941.0249 -11078.0281 -10405.1009 -12002.8881 -11156.8685
121 122 123 124 125 126
-11928.4702 -10865.3644 -11213.2301 -10311.3701 -11013.4744 -12071.7149
127 128 129 130 131 132
-11097.0398 -32785.7374 -31684.9904 27735.0273 17411.7203 17904.8315
133 134 135 136 137 138
-38256.9643 -16557.0026 -14364.6072 -13816.9600 -13621.1866 -27313.8478
139 140 141 142 143 144
-14497.6276 -4995.3696 -4319.3008 -13709.3539 -12558.8073 -13509.6039
145 146 147 148 149 150
-13759.2662 -13939.2648 -146690.8873 -151273.1997 -147761.4556 -136876.5586
151 152 153 154 155 156
-142916.9026 150666.9555 298617.3353 286743.2309 -3397.6479 -9673.8952
157 158 159 160 161 162
-10523.2124 -11103.0672 27824.0100 -38383.9976 -5515.8671 -5563.1901
163 164 165 166 167 168
-3116.3550 -4564.9795 -9772.4107 -14175.2629 -13946.5282 165295.8507
169 170 171 172 173 174
290928.2368 8665.4050 -3424.0406 -10833.1352 -11189.8796 -8193.6683
175 176 177 178 179 180
-11103.1552 38768.3065 29689.7823 16382.8778 33469.1296 25493.1049
181 182 183 184 185 186
22203.2177 25154.5139 20676.3832 27709.4339 24352.3221 -38496.4683
187 188 189 190 191 192
-15178.2771 -4652.7249 -5294.1555 -12570.6859 -13732.3666 -189742.3901
193 194 195 196 197 198
168268.5537 -642.2560 -3392.1365 -11102.2531 17912.8496 50077.8428
199 200 201 202 203 204
22480.7138 27926.3424 22415.2256 -38968.1121 -12785.8360 -3827.1839
205 206 207 208 209 210
-11072.6088 -31365.8697 -33763.8745 -33738.4296 27935.4232 22316.9405
211 212 213 214 215 216
22120.2418 -38625.0343 -12664.3209 -14668.8941 30883.9432 -5369.7218
217 218 219 220 221 222
297.7586 -5048.6971 -3530.7233 -5129.4265 -13226.1934 -13931.4692
223 224 225 226 227 228
-140269.1905 -149337.3644 156115.8886 299908.8713 112093.8766 -7648.3683
229 230 231 232 233 234
-10631.9757 -11103.4633 27824.1353 -38249.6540 -16084.1128 -635.4399
235 236 237 238 239 240
-6115.0196 -5320.8986 -13928.0376 -144329.2619 168268.5683 6711.1829
241 242 243 244 245 246
-3418.3413 -9843.5010 -11298.9017 -11849.3718 -11096.4140 -30441.0379
247 248 249 250 251 252
-27792.0379 27823.8222 -38351.4679 -5093.7170 -13387.5614 -13890.6039
253 254 255 256 257 258
-13945.9322 168268.8762 15098.3623 -3458.6759 -10536.5767 -11101.4609
259 260 261 262 263 264
27823.5021 -38265.2582 -13583.6840 -4921.4852 -4988.3175 -3996.9991
265 266 267 268 269 270
-13509.4897 -13706.9800 -13946.1309 166404.0359 301756.3289 7330.5145
271 272 273 274 275 276
-3412.2202 -10633.8450 -10378.9590 -11041.4591 -30591.8150 29948.8615
277 278 279 280 281 282
41948.5409 18262.8119 -38229.6500 -14207.2020 -13226.3638 -5321.0494
283 284 285 286 287 288
-13937.3295 166051.9954 276657.2340 281675.7833 -3398.3637 -10917.7235
289 290 291 292 293 294
-11112.9927 -11090.9554 -34140.8152 -30469.1148 27823.9057 -38016.2974
295 296 297 298 299 300
-12838.7941 -2755.3435 32452.2003 -3815.9283 -2340.5326 -5632.4148
301 302 303 304 305 306
78793.0612 -5319.9433 -13946.1929 168267.3517 -3395.1691 -11099.8613
307 308 309 310 311 312
-32477.2195 27823.7317 -38366.2030 -5319.8678 -13945.9943 168268.7882
313 314 315 316 317 318
19295.4177 54217.6282 10080.3950 -2644.9526 -3389.0318 -8273.5090
319 320 321 322 323 324
-10875.8196 -11110.3020 -11103.6393 27824.2118 -38217.4694 -12830.3363
325 326 327 328 329 330
-4851.8719 -5322.3566 -13896.6187 -148469.3140 -148028.3140 -146325.5946
331 332 333 334 335 336
-151000.3286 -146817.7818 168268.7442 -467.8752 -3399.6751 -11102.2311
337 338 339 340 341 342
27825.7082 21467.4703 -38337.9865 -6152.3353 -4789.0893 -5319.6165
343 344 345 346 347 348
-13946.1805 168269.0228 2106.9224 -3423.6655 -11025.5434 -10727.7474
349 350 351 352 353 354
-8495.3349 -11101.9010 23284.7455 23577.8067 21902.5861 22016.0832
355 356 357 358 359 360
-38364.2154 -5207.4807 -6303.2821 -5242.1192 -11249.1952 -13166.5866
361 362 363 364 365 366
-12632.1763 -13560.5184 -12494.1668 -13379.5186 -13244.0418 -12926.6095
367 368 369 370 371 372
-12687.7083 -10314.2730 -12678.8124 -8699.3621 -12110.9923 -13447.8481
373 374 375 376 377 378
-11403.8518 -12899.4066 -12490.3652 -13083.9550 -12045.9084 -12913.3510
379 380 381 382 383 384
-12201.9231 -13048.5309 -12702.6642 -13252.4426 -12937.0619 -11691.9427
385 386 387 388 389 390
-6337.8772 -12491.0807 -10283.3339 -12587.6097 -11719.2953 -9698.3613
391 392 393 394 395 396
-13768.4306 -12369.5281 -12827.7549 -12034.4520 -9948.6925 -1745.4977
397 398 399 400 401 402
-5871.5347 -9481.0027 -8563.1971 -9195.3685 23310.0873 -3322.0735
403 404 405 406 407 408
2004.8449 3973.2002 150298.5544 -24029.0406 45980.9817 45618.1761
409 410 411 412 413 414
39841.0064 21399.4789 19376.4440 24599.6358 15053.0460 -2429.6048
415 416 417 418 419 420
4023.7643 -3725.8022 -688.5390 21036.9928 10968.0252 -11263.2979
421 422 423 424 425 426
-10675.1240 9121.2801 20540.6647 -1221.1366 -1079.6174 5009.0938
427 428 429 430 431
-3003.4174 1141.6817 -6418.6668 3050.9329 11984.6767
> postscript(file="/var/www/html/rcomp/tmp/632m11291301535.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 67597.3278 NA
1 39723.7974 67597.3278
2 22636.8972 39723.7974
3 -14077.8071 22636.8972
4 -5704.0262 -14077.8071
5 -3969.4508 -5704.0262
6 -5468.7088 -3969.4508
7 -3171.3922 -5468.7088
8 -8007.7062 -3171.3922
9 -3125.1589 -8007.7062
10 -7744.1998 -3125.1589
11 -5891.4039 -7744.1998
12 -4401.0919 -5891.4039
13 -14175.6398 -4401.0919
14 -8909.2994 -14175.6398
15 -6690.3245 -8909.2994
16 -10436.7778 -6690.3245
17 -9343.1248 -10436.7778
18 -9182.1943 -9343.1248
19 -11682.5550 -9182.1943
20 -8262.8788 -11682.5550
21 -10036.2001 -8262.8788
22 -10003.2218 -10036.2001
23 -12271.3588 -10003.2218
24 -12014.9410 -12271.3588
25 -10332.9962 -12014.9410
26 -11163.7273 -10332.9962
27 -12357.7540 -11163.7273
28 -8075.8499 -12357.7540
29 -11345.4600 -8075.8499
30 -8604.2072 -11345.4600
31 -10504.8608 -8604.2072
32 -8691.7598 -10504.8608
33 -9663.8672 -8691.7598
34 -10272.0322 -9663.8672
35 -12412.0269 -10272.0322
36 -9907.2548 -12412.0269
37 -9277.2744 -9907.2548
38 -10552.1168 -9277.2744
39 -12171.3138 -10552.1168
40 -12597.3304 -12171.3138
41 -12146.0802 -12597.3304
42 -10506.8167 -12146.0802
43 -11570.2580 -10506.8167
44 -10959.5187 -11570.2580
45 -11066.6252 -10959.5187
46 -12199.4867 -11066.6252
47 -11574.0036 -12199.4867
48 -10864.9241 -11574.0036
49 -12095.7642 -10864.9241
50 -10258.6656 -12095.7642
51 -10641.4424 -10258.6656
52 -11424.6950 -10641.4424
53 -10438.7596 -11424.6950
54 -12652.7717 -10438.7596
55 -10674.8294 -12652.7717
56 -7931.6590 -10674.8294
57 -10804.9200 -7931.6590
58 -11280.4274 -10804.9200
59 -11554.5445 -11280.4274
60 -10513.2053 -11554.5445
61 -12281.2829 -10513.2053
62 -10800.0118 -12281.2829
63 -11520.8955 -10800.0118
64 -11011.7762 -11520.8955
65 -12479.8449 -11011.7762
66 -11622.4581 -12479.8449
67 -10937.8147 -11622.4581
68 -11301.7111 -10937.8147
69 -10934.6249 -11301.7111
70 -11809.5967 -10934.6249
71 -10878.5941 -11809.5967
72 -11469.1560 -10878.5941
73 -10903.3164 -11469.1560
74 -11025.8128 -10903.3164
75 -11705.7003 -11025.8128
76 -11141.3152 -11705.7003
77 -11715.7240 -11141.3152
78 -11282.1145 -11715.7240
79 -11430.9917 -11282.1145
80 -10864.7099 -11430.9917
81 -10934.9377 -10864.7099
82 -11424.5128 -10934.9377
83 -10945.9050 -11424.5128
84 -10741.8724 -10945.9050
85 -11286.5622 -10741.8724
86 -10658.9196 -11286.5622
87 -10917.5817 -10658.9196
88 -11596.1079 -10917.5817
89 -10569.9313 -11596.1079
90 -10936.9838 -10569.9313
91 -10998.5909 -10936.9838
92 -10418.7580 -10998.5909
93 -11752.1999 -10418.7580
94 -10745.7005 -11752.1999
95 -10872.2244 -10745.7005
96 -10060.0624 -10872.2244
97 -11244.0322 -10060.0624
98 -11015.0061 -11244.0322
99 -11011.1878 -11015.0061
100 -11277.6763 -11011.1878
101 -11841.9552 -11277.6763
102 -11796.1607 -11841.9552
103 -11362.0832 -11796.1607
104 -10871.1308 -11362.0832
105 -11442.5885 -10871.1308
106 -12577.4531 -11442.5885
107 -9927.2119 -12577.4531
108 -11777.8299 -9927.2119
109 -11718.1146 -11777.8299
110 -10481.7946 -11718.1146
111 -12296.7286 -10481.7946
112 -11826.3110 -12296.7286
113 -11387.2132 -11826.3110
114 -10309.3095 -11387.2132
115 -10941.0249 -10309.3095
116 -11078.0281 -10941.0249
117 -10405.1009 -11078.0281
118 -12002.8881 -10405.1009
119 -11156.8685 -12002.8881
120 -11928.4702 -11156.8685
121 -10865.3644 -11928.4702
122 -11213.2301 -10865.3644
123 -10311.3701 -11213.2301
124 -11013.4744 -10311.3701
125 -12071.7149 -11013.4744
126 -11097.0398 -12071.7149
127 -32785.7374 -11097.0398
128 -31684.9904 -32785.7374
129 27735.0273 -31684.9904
130 17411.7203 27735.0273
131 17904.8315 17411.7203
132 -38256.9643 17904.8315
133 -16557.0026 -38256.9643
134 -14364.6072 -16557.0026
135 -13816.9600 -14364.6072
136 -13621.1866 -13816.9600
137 -27313.8478 -13621.1866
138 -14497.6276 -27313.8478
139 -4995.3696 -14497.6276
140 -4319.3008 -4995.3696
141 -13709.3539 -4319.3008
142 -12558.8073 -13709.3539
143 -13509.6039 -12558.8073
144 -13759.2662 -13509.6039
145 -13939.2648 -13759.2662
146 -146690.8873 -13939.2648
147 -151273.1997 -146690.8873
148 -147761.4556 -151273.1997
149 -136876.5586 -147761.4556
150 -142916.9026 -136876.5586
151 150666.9555 -142916.9026
152 298617.3353 150666.9555
153 286743.2309 298617.3353
154 -3397.6479 286743.2309
155 -9673.8952 -3397.6479
156 -10523.2124 -9673.8952
157 -11103.0672 -10523.2124
158 27824.0100 -11103.0672
159 -38383.9976 27824.0100
160 -5515.8671 -38383.9976
161 -5563.1901 -5515.8671
162 -3116.3550 -5563.1901
163 -4564.9795 -3116.3550
164 -9772.4107 -4564.9795
165 -14175.2629 -9772.4107
166 -13946.5282 -14175.2629
167 165295.8507 -13946.5282
168 290928.2368 165295.8507
169 8665.4050 290928.2368
170 -3424.0406 8665.4050
171 -10833.1352 -3424.0406
172 -11189.8796 -10833.1352
173 -8193.6683 -11189.8796
174 -11103.1552 -8193.6683
175 38768.3065 -11103.1552
176 29689.7823 38768.3065
177 16382.8778 29689.7823
178 33469.1296 16382.8778
179 25493.1049 33469.1296
180 22203.2177 25493.1049
181 25154.5139 22203.2177
182 20676.3832 25154.5139
183 27709.4339 20676.3832
184 24352.3221 27709.4339
185 -38496.4683 24352.3221
186 -15178.2771 -38496.4683
187 -4652.7249 -15178.2771
188 -5294.1555 -4652.7249
189 -12570.6859 -5294.1555
190 -13732.3666 -12570.6859
191 -189742.3901 -13732.3666
192 168268.5537 -189742.3901
193 -642.2560 168268.5537
194 -3392.1365 -642.2560
195 -11102.2531 -3392.1365
196 17912.8496 -11102.2531
197 50077.8428 17912.8496
198 22480.7138 50077.8428
199 27926.3424 22480.7138
200 22415.2256 27926.3424
201 -38968.1121 22415.2256
202 -12785.8360 -38968.1121
203 -3827.1839 -12785.8360
204 -11072.6088 -3827.1839
205 -31365.8697 -11072.6088
206 -33763.8745 -31365.8697
207 -33738.4296 -33763.8745
208 27935.4232 -33738.4296
209 22316.9405 27935.4232
210 22120.2418 22316.9405
211 -38625.0343 22120.2418
212 -12664.3209 -38625.0343
213 -14668.8941 -12664.3209
214 30883.9432 -14668.8941
215 -5369.7218 30883.9432
216 297.7586 -5369.7218
217 -5048.6971 297.7586
218 -3530.7233 -5048.6971
219 -5129.4265 -3530.7233
220 -13226.1934 -5129.4265
221 -13931.4692 -13226.1934
222 -140269.1905 -13931.4692
223 -149337.3644 -140269.1905
224 156115.8886 -149337.3644
225 299908.8713 156115.8886
226 112093.8766 299908.8713
227 -7648.3683 112093.8766
228 -10631.9757 -7648.3683
229 -11103.4633 -10631.9757
230 27824.1353 -11103.4633
231 -38249.6540 27824.1353
232 -16084.1128 -38249.6540
233 -635.4399 -16084.1128
234 -6115.0196 -635.4399
235 -5320.8986 -6115.0196
236 -13928.0376 -5320.8986
237 -144329.2619 -13928.0376
238 168268.5683 -144329.2619
239 6711.1829 168268.5683
240 -3418.3413 6711.1829
241 -9843.5010 -3418.3413
242 -11298.9017 -9843.5010
243 -11849.3718 -11298.9017
244 -11096.4140 -11849.3718
245 -30441.0379 -11096.4140
246 -27792.0379 -30441.0379
247 27823.8222 -27792.0379
248 -38351.4679 27823.8222
249 -5093.7170 -38351.4679
250 -13387.5614 -5093.7170
251 -13890.6039 -13387.5614
252 -13945.9322 -13890.6039
253 168268.8762 -13945.9322
254 15098.3623 168268.8762
255 -3458.6759 15098.3623
256 -10536.5767 -3458.6759
257 -11101.4609 -10536.5767
258 27823.5021 -11101.4609
259 -38265.2582 27823.5021
260 -13583.6840 -38265.2582
261 -4921.4852 -13583.6840
262 -4988.3175 -4921.4852
263 -3996.9991 -4988.3175
264 -13509.4897 -3996.9991
265 -13706.9800 -13509.4897
266 -13946.1309 -13706.9800
267 166404.0359 -13946.1309
268 301756.3289 166404.0359
269 7330.5145 301756.3289
270 -3412.2202 7330.5145
271 -10633.8450 -3412.2202
272 -10378.9590 -10633.8450
273 -11041.4591 -10378.9590
274 -30591.8150 -11041.4591
275 29948.8615 -30591.8150
276 41948.5409 29948.8615
277 18262.8119 41948.5409
278 -38229.6500 18262.8119
279 -14207.2020 -38229.6500
280 -13226.3638 -14207.2020
281 -5321.0494 -13226.3638
282 -13937.3295 -5321.0494
283 166051.9954 -13937.3295
284 276657.2340 166051.9954
285 281675.7833 276657.2340
286 -3398.3637 281675.7833
287 -10917.7235 -3398.3637
288 -11112.9927 -10917.7235
289 -11090.9554 -11112.9927
290 -34140.8152 -11090.9554
291 -30469.1148 -34140.8152
292 27823.9057 -30469.1148
293 -38016.2974 27823.9057
294 -12838.7941 -38016.2974
295 -2755.3435 -12838.7941
296 32452.2003 -2755.3435
297 -3815.9283 32452.2003
298 -2340.5326 -3815.9283
299 -5632.4148 -2340.5326
300 78793.0612 -5632.4148
301 -5319.9433 78793.0612
302 -13946.1929 -5319.9433
303 168267.3517 -13946.1929
304 -3395.1691 168267.3517
305 -11099.8613 -3395.1691
306 -32477.2195 -11099.8613
307 27823.7317 -32477.2195
308 -38366.2030 27823.7317
309 -5319.8678 -38366.2030
310 -13945.9943 -5319.8678
311 168268.7882 -13945.9943
312 19295.4177 168268.7882
313 54217.6282 19295.4177
314 10080.3950 54217.6282
315 -2644.9526 10080.3950
316 -3389.0318 -2644.9526
317 -8273.5090 -3389.0318
318 -10875.8196 -8273.5090
319 -11110.3020 -10875.8196
320 -11103.6393 -11110.3020
321 27824.2118 -11103.6393
322 -38217.4694 27824.2118
323 -12830.3363 -38217.4694
324 -4851.8719 -12830.3363
325 -5322.3566 -4851.8719
326 -13896.6187 -5322.3566
327 -148469.3140 -13896.6187
328 -148028.3140 -148469.3140
329 -146325.5946 -148028.3140
330 -151000.3286 -146325.5946
331 -146817.7818 -151000.3286
332 168268.7442 -146817.7818
333 -467.8752 168268.7442
334 -3399.6751 -467.8752
335 -11102.2311 -3399.6751
336 27825.7082 -11102.2311
337 21467.4703 27825.7082
338 -38337.9865 21467.4703
339 -6152.3353 -38337.9865
340 -4789.0893 -6152.3353
341 -5319.6165 -4789.0893
342 -13946.1805 -5319.6165
343 168269.0228 -13946.1805
344 2106.9224 168269.0228
345 -3423.6655 2106.9224
346 -11025.5434 -3423.6655
347 -10727.7474 -11025.5434
348 -8495.3349 -10727.7474
349 -11101.9010 -8495.3349
350 23284.7455 -11101.9010
351 23577.8067 23284.7455
352 21902.5861 23577.8067
353 22016.0832 21902.5861
354 -38364.2154 22016.0832
355 -5207.4807 -38364.2154
356 -6303.2821 -5207.4807
357 -5242.1192 -6303.2821
358 -11249.1952 -5242.1192
359 -13166.5866 -11249.1952
360 -12632.1763 -13166.5866
361 -13560.5184 -12632.1763
362 -12494.1668 -13560.5184
363 -13379.5186 -12494.1668
364 -13244.0418 -13379.5186
365 -12926.6095 -13244.0418
366 -12687.7083 -12926.6095
367 -10314.2730 -12687.7083
368 -12678.8124 -10314.2730
369 -8699.3621 -12678.8124
370 -12110.9923 -8699.3621
371 -13447.8481 -12110.9923
372 -11403.8518 -13447.8481
373 -12899.4066 -11403.8518
374 -12490.3652 -12899.4066
375 -13083.9550 -12490.3652
376 -12045.9084 -13083.9550
377 -12913.3510 -12045.9084
378 -12201.9231 -12913.3510
379 -13048.5309 -12201.9231
380 -12702.6642 -13048.5309
381 -13252.4426 -12702.6642
382 -12937.0619 -13252.4426
383 -11691.9427 -12937.0619
384 -6337.8772 -11691.9427
385 -12491.0807 -6337.8772
386 -10283.3339 -12491.0807
387 -12587.6097 -10283.3339
388 -11719.2953 -12587.6097
389 -9698.3613 -11719.2953
390 -13768.4306 -9698.3613
391 -12369.5281 -13768.4306
392 -12827.7549 -12369.5281
393 -12034.4520 -12827.7549
394 -9948.6925 -12034.4520
395 -1745.4977 -9948.6925
396 -5871.5347 -1745.4977
397 -9481.0027 -5871.5347
398 -8563.1971 -9481.0027
399 -9195.3685 -8563.1971
400 23310.0873 -9195.3685
401 -3322.0735 23310.0873
402 2004.8449 -3322.0735
403 3973.2002 2004.8449
404 150298.5544 3973.2002
405 -24029.0406 150298.5544
406 45980.9817 -24029.0406
407 45618.1761 45980.9817
408 39841.0064 45618.1761
409 21399.4789 39841.0064
410 19376.4440 21399.4789
411 24599.6358 19376.4440
412 15053.0460 24599.6358
413 -2429.6048 15053.0460
414 4023.7643 -2429.6048
415 -3725.8022 4023.7643
416 -688.5390 -3725.8022
417 21036.9928 -688.5390
418 10968.0252 21036.9928
419 -11263.2979 10968.0252
420 -10675.1240 -11263.2979
421 9121.2801 -10675.1240
422 20540.6647 9121.2801
423 -1221.1366 20540.6647
424 -1079.6174 -1221.1366
425 5009.0938 -1079.6174
426 -3003.4174 5009.0938
427 1141.6817 -3003.4174
428 -6418.6668 1141.6817
429 3050.9329 -6418.6668
430 11984.6767 3050.9329
431 NA 11984.6767
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 39723.7974 67597.3278
[2,] 22636.8972 39723.7974
[3,] -14077.8071 22636.8972
[4,] -5704.0262 -14077.8071
[5,] -3969.4508 -5704.0262
[6,] -5468.7088 -3969.4508
[7,] -3171.3922 -5468.7088
[8,] -8007.7062 -3171.3922
[9,] -3125.1589 -8007.7062
[10,] -7744.1998 -3125.1589
[11,] -5891.4039 -7744.1998
[12,] -4401.0919 -5891.4039
[13,] -14175.6398 -4401.0919
[14,] -8909.2994 -14175.6398
[15,] -6690.3245 -8909.2994
[16,] -10436.7778 -6690.3245
[17,] -9343.1248 -10436.7778
[18,] -9182.1943 -9343.1248
[19,] -11682.5550 -9182.1943
[20,] -8262.8788 -11682.5550
[21,] -10036.2001 -8262.8788
[22,] -10003.2218 -10036.2001
[23,] -12271.3588 -10003.2218
[24,] -12014.9410 -12271.3588
[25,] -10332.9962 -12014.9410
[26,] -11163.7273 -10332.9962
[27,] -12357.7540 -11163.7273
[28,] -8075.8499 -12357.7540
[29,] -11345.4600 -8075.8499
[30,] -8604.2072 -11345.4600
[31,] -10504.8608 -8604.2072
[32,] -8691.7598 -10504.8608
[33,] -9663.8672 -8691.7598
[34,] -10272.0322 -9663.8672
[35,] -12412.0269 -10272.0322
[36,] -9907.2548 -12412.0269
[37,] -9277.2744 -9907.2548
[38,] -10552.1168 -9277.2744
[39,] -12171.3138 -10552.1168
[40,] -12597.3304 -12171.3138
[41,] -12146.0802 -12597.3304
[42,] -10506.8167 -12146.0802
[43,] -11570.2580 -10506.8167
[44,] -10959.5187 -11570.2580
[45,] -11066.6252 -10959.5187
[46,] -12199.4867 -11066.6252
[47,] -11574.0036 -12199.4867
[48,] -10864.9241 -11574.0036
[49,] -12095.7642 -10864.9241
[50,] -10258.6656 -12095.7642
[51,] -10641.4424 -10258.6656
[52,] -11424.6950 -10641.4424
[53,] -10438.7596 -11424.6950
[54,] -12652.7717 -10438.7596
[55,] -10674.8294 -12652.7717
[56,] -7931.6590 -10674.8294
[57,] -10804.9200 -7931.6590
[58,] -11280.4274 -10804.9200
[59,] -11554.5445 -11280.4274
[60,] -10513.2053 -11554.5445
[61,] -12281.2829 -10513.2053
[62,] -10800.0118 -12281.2829
[63,] -11520.8955 -10800.0118
[64,] -11011.7762 -11520.8955
[65,] -12479.8449 -11011.7762
[66,] -11622.4581 -12479.8449
[67,] -10937.8147 -11622.4581
[68,] -11301.7111 -10937.8147
[69,] -10934.6249 -11301.7111
[70,] -11809.5967 -10934.6249
[71,] -10878.5941 -11809.5967
[72,] -11469.1560 -10878.5941
[73,] -10903.3164 -11469.1560
[74,] -11025.8128 -10903.3164
[75,] -11705.7003 -11025.8128
[76,] -11141.3152 -11705.7003
[77,] -11715.7240 -11141.3152
[78,] -11282.1145 -11715.7240
[79,] -11430.9917 -11282.1145
[80,] -10864.7099 -11430.9917
[81,] -10934.9377 -10864.7099
[82,] -11424.5128 -10934.9377
[83,] -10945.9050 -11424.5128
[84,] -10741.8724 -10945.9050
[85,] -11286.5622 -10741.8724
[86,] -10658.9196 -11286.5622
[87,] -10917.5817 -10658.9196
[88,] -11596.1079 -10917.5817
[89,] -10569.9313 -11596.1079
[90,] -10936.9838 -10569.9313
[91,] -10998.5909 -10936.9838
[92,] -10418.7580 -10998.5909
[93,] -11752.1999 -10418.7580
[94,] -10745.7005 -11752.1999
[95,] -10872.2244 -10745.7005
[96,] -10060.0624 -10872.2244
[97,] -11244.0322 -10060.0624
[98,] -11015.0061 -11244.0322
[99,] -11011.1878 -11015.0061
[100,] -11277.6763 -11011.1878
[101,] -11841.9552 -11277.6763
[102,] -11796.1607 -11841.9552
[103,] -11362.0832 -11796.1607
[104,] -10871.1308 -11362.0832
[105,] -11442.5885 -10871.1308
[106,] -12577.4531 -11442.5885
[107,] -9927.2119 -12577.4531
[108,] -11777.8299 -9927.2119
[109,] -11718.1146 -11777.8299
[110,] -10481.7946 -11718.1146
[111,] -12296.7286 -10481.7946
[112,] -11826.3110 -12296.7286
[113,] -11387.2132 -11826.3110
[114,] -10309.3095 -11387.2132
[115,] -10941.0249 -10309.3095
[116,] -11078.0281 -10941.0249
[117,] -10405.1009 -11078.0281
[118,] -12002.8881 -10405.1009
[119,] -11156.8685 -12002.8881
[120,] -11928.4702 -11156.8685
[121,] -10865.3644 -11928.4702
[122,] -11213.2301 -10865.3644
[123,] -10311.3701 -11213.2301
[124,] -11013.4744 -10311.3701
[125,] -12071.7149 -11013.4744
[126,] -11097.0398 -12071.7149
[127,] -32785.7374 -11097.0398
[128,] -31684.9904 -32785.7374
[129,] 27735.0273 -31684.9904
[130,] 17411.7203 27735.0273
[131,] 17904.8315 17411.7203
[132,] -38256.9643 17904.8315
[133,] -16557.0026 -38256.9643
[134,] -14364.6072 -16557.0026
[135,] -13816.9600 -14364.6072
[136,] -13621.1866 -13816.9600
[137,] -27313.8478 -13621.1866
[138,] -14497.6276 -27313.8478
[139,] -4995.3696 -14497.6276
[140,] -4319.3008 -4995.3696
[141,] -13709.3539 -4319.3008
[142,] -12558.8073 -13709.3539
[143,] -13509.6039 -12558.8073
[144,] -13759.2662 -13509.6039
[145,] -13939.2648 -13759.2662
[146,] -146690.8873 -13939.2648
[147,] -151273.1997 -146690.8873
[148,] -147761.4556 -151273.1997
[149,] -136876.5586 -147761.4556
[150,] -142916.9026 -136876.5586
[151,] 150666.9555 -142916.9026
[152,] 298617.3353 150666.9555
[153,] 286743.2309 298617.3353
[154,] -3397.6479 286743.2309
[155,] -9673.8952 -3397.6479
[156,] -10523.2124 -9673.8952
[157,] -11103.0672 -10523.2124
[158,] 27824.0100 -11103.0672
[159,] -38383.9976 27824.0100
[160,] -5515.8671 -38383.9976
[161,] -5563.1901 -5515.8671
[162,] -3116.3550 -5563.1901
[163,] -4564.9795 -3116.3550
[164,] -9772.4107 -4564.9795
[165,] -14175.2629 -9772.4107
[166,] -13946.5282 -14175.2629
[167,] 165295.8507 -13946.5282
[168,] 290928.2368 165295.8507
[169,] 8665.4050 290928.2368
[170,] -3424.0406 8665.4050
[171,] -10833.1352 -3424.0406
[172,] -11189.8796 -10833.1352
[173,] -8193.6683 -11189.8796
[174,] -11103.1552 -8193.6683
[175,] 38768.3065 -11103.1552
[176,] 29689.7823 38768.3065
[177,] 16382.8778 29689.7823
[178,] 33469.1296 16382.8778
[179,] 25493.1049 33469.1296
[180,] 22203.2177 25493.1049
[181,] 25154.5139 22203.2177
[182,] 20676.3832 25154.5139
[183,] 27709.4339 20676.3832
[184,] 24352.3221 27709.4339
[185,] -38496.4683 24352.3221
[186,] -15178.2771 -38496.4683
[187,] -4652.7249 -15178.2771
[188,] -5294.1555 -4652.7249
[189,] -12570.6859 -5294.1555
[190,] -13732.3666 -12570.6859
[191,] -189742.3901 -13732.3666
[192,] 168268.5537 -189742.3901
[193,] -642.2560 168268.5537
[194,] -3392.1365 -642.2560
[195,] -11102.2531 -3392.1365
[196,] 17912.8496 -11102.2531
[197,] 50077.8428 17912.8496
[198,] 22480.7138 50077.8428
[199,] 27926.3424 22480.7138
[200,] 22415.2256 27926.3424
[201,] -38968.1121 22415.2256
[202,] -12785.8360 -38968.1121
[203,] -3827.1839 -12785.8360
[204,] -11072.6088 -3827.1839
[205,] -31365.8697 -11072.6088
[206,] -33763.8745 -31365.8697
[207,] -33738.4296 -33763.8745
[208,] 27935.4232 -33738.4296
[209,] 22316.9405 27935.4232
[210,] 22120.2418 22316.9405
[211,] -38625.0343 22120.2418
[212,] -12664.3209 -38625.0343
[213,] -14668.8941 -12664.3209
[214,] 30883.9432 -14668.8941
[215,] -5369.7218 30883.9432
[216,] 297.7586 -5369.7218
[217,] -5048.6971 297.7586
[218,] -3530.7233 -5048.6971
[219,] -5129.4265 -3530.7233
[220,] -13226.1934 -5129.4265
[221,] -13931.4692 -13226.1934
[222,] -140269.1905 -13931.4692
[223,] -149337.3644 -140269.1905
[224,] 156115.8886 -149337.3644
[225,] 299908.8713 156115.8886
[226,] 112093.8766 299908.8713
[227,] -7648.3683 112093.8766
[228,] -10631.9757 -7648.3683
[229,] -11103.4633 -10631.9757
[230,] 27824.1353 -11103.4633
[231,] -38249.6540 27824.1353
[232,] -16084.1128 -38249.6540
[233,] -635.4399 -16084.1128
[234,] -6115.0196 -635.4399
[235,] -5320.8986 -6115.0196
[236,] -13928.0376 -5320.8986
[237,] -144329.2619 -13928.0376
[238,] 168268.5683 -144329.2619
[239,] 6711.1829 168268.5683
[240,] -3418.3413 6711.1829
[241,] -9843.5010 -3418.3413
[242,] -11298.9017 -9843.5010
[243,] -11849.3718 -11298.9017
[244,] -11096.4140 -11849.3718
[245,] -30441.0379 -11096.4140
[246,] -27792.0379 -30441.0379
[247,] 27823.8222 -27792.0379
[248,] -38351.4679 27823.8222
[249,] -5093.7170 -38351.4679
[250,] -13387.5614 -5093.7170
[251,] -13890.6039 -13387.5614
[252,] -13945.9322 -13890.6039
[253,] 168268.8762 -13945.9322
[254,] 15098.3623 168268.8762
[255,] -3458.6759 15098.3623
[256,] -10536.5767 -3458.6759
[257,] -11101.4609 -10536.5767
[258,] 27823.5021 -11101.4609
[259,] -38265.2582 27823.5021
[260,] -13583.6840 -38265.2582
[261,] -4921.4852 -13583.6840
[262,] -4988.3175 -4921.4852
[263,] -3996.9991 -4988.3175
[264,] -13509.4897 -3996.9991
[265,] -13706.9800 -13509.4897
[266,] -13946.1309 -13706.9800
[267,] 166404.0359 -13946.1309
[268,] 301756.3289 166404.0359
[269,] 7330.5145 301756.3289
[270,] -3412.2202 7330.5145
[271,] -10633.8450 -3412.2202
[272,] -10378.9590 -10633.8450
[273,] -11041.4591 -10378.9590
[274,] -30591.8150 -11041.4591
[275,] 29948.8615 -30591.8150
[276,] 41948.5409 29948.8615
[277,] 18262.8119 41948.5409
[278,] -38229.6500 18262.8119
[279,] -14207.2020 -38229.6500
[280,] -13226.3638 -14207.2020
[281,] -5321.0494 -13226.3638
[282,] -13937.3295 -5321.0494
[283,] 166051.9954 -13937.3295
[284,] 276657.2340 166051.9954
[285,] 281675.7833 276657.2340
[286,] -3398.3637 281675.7833
[287,] -10917.7235 -3398.3637
[288,] -11112.9927 -10917.7235
[289,] -11090.9554 -11112.9927
[290,] -34140.8152 -11090.9554
[291,] -30469.1148 -34140.8152
[292,] 27823.9057 -30469.1148
[293,] -38016.2974 27823.9057
[294,] -12838.7941 -38016.2974
[295,] -2755.3435 -12838.7941
[296,] 32452.2003 -2755.3435
[297,] -3815.9283 32452.2003
[298,] -2340.5326 -3815.9283
[299,] -5632.4148 -2340.5326
[300,] 78793.0612 -5632.4148
[301,] -5319.9433 78793.0612
[302,] -13946.1929 -5319.9433
[303,] 168267.3517 -13946.1929
[304,] -3395.1691 168267.3517
[305,] -11099.8613 -3395.1691
[306,] -32477.2195 -11099.8613
[307,] 27823.7317 -32477.2195
[308,] -38366.2030 27823.7317
[309,] -5319.8678 -38366.2030
[310,] -13945.9943 -5319.8678
[311,] 168268.7882 -13945.9943
[312,] 19295.4177 168268.7882
[313,] 54217.6282 19295.4177
[314,] 10080.3950 54217.6282
[315,] -2644.9526 10080.3950
[316,] -3389.0318 -2644.9526
[317,] -8273.5090 -3389.0318
[318,] -10875.8196 -8273.5090
[319,] -11110.3020 -10875.8196
[320,] -11103.6393 -11110.3020
[321,] 27824.2118 -11103.6393
[322,] -38217.4694 27824.2118
[323,] -12830.3363 -38217.4694
[324,] -4851.8719 -12830.3363
[325,] -5322.3566 -4851.8719
[326,] -13896.6187 -5322.3566
[327,] -148469.3140 -13896.6187
[328,] -148028.3140 -148469.3140
[329,] -146325.5946 -148028.3140
[330,] -151000.3286 -146325.5946
[331,] -146817.7818 -151000.3286
[332,] 168268.7442 -146817.7818
[333,] -467.8752 168268.7442
[334,] -3399.6751 -467.8752
[335,] -11102.2311 -3399.6751
[336,] 27825.7082 -11102.2311
[337,] 21467.4703 27825.7082
[338,] -38337.9865 21467.4703
[339,] -6152.3353 -38337.9865
[340,] -4789.0893 -6152.3353
[341,] -5319.6165 -4789.0893
[342,] -13946.1805 -5319.6165
[343,] 168269.0228 -13946.1805
[344,] 2106.9224 168269.0228
[345,] -3423.6655 2106.9224
[346,] -11025.5434 -3423.6655
[347,] -10727.7474 -11025.5434
[348,] -8495.3349 -10727.7474
[349,] -11101.9010 -8495.3349
[350,] 23284.7455 -11101.9010
[351,] 23577.8067 23284.7455
[352,] 21902.5861 23577.8067
[353,] 22016.0832 21902.5861
[354,] -38364.2154 22016.0832
[355,] -5207.4807 -38364.2154
[356,] -6303.2821 -5207.4807
[357,] -5242.1192 -6303.2821
[358,] -11249.1952 -5242.1192
[359,] -13166.5866 -11249.1952
[360,] -12632.1763 -13166.5866
[361,] -13560.5184 -12632.1763
[362,] -12494.1668 -13560.5184
[363,] -13379.5186 -12494.1668
[364,] -13244.0418 -13379.5186
[365,] -12926.6095 -13244.0418
[366,] -12687.7083 -12926.6095
[367,] -10314.2730 -12687.7083
[368,] -12678.8124 -10314.2730
[369,] -8699.3621 -12678.8124
[370,] -12110.9923 -8699.3621
[371,] -13447.8481 -12110.9923
[372,] -11403.8518 -13447.8481
[373,] -12899.4066 -11403.8518
[374,] -12490.3652 -12899.4066
[375,] -13083.9550 -12490.3652
[376,] -12045.9084 -13083.9550
[377,] -12913.3510 -12045.9084
[378,] -12201.9231 -12913.3510
[379,] -13048.5309 -12201.9231
[380,] -12702.6642 -13048.5309
[381,] -13252.4426 -12702.6642
[382,] -12937.0619 -13252.4426
[383,] -11691.9427 -12937.0619
[384,] -6337.8772 -11691.9427
[385,] -12491.0807 -6337.8772
[386,] -10283.3339 -12491.0807
[387,] -12587.6097 -10283.3339
[388,] -11719.2953 -12587.6097
[389,] -9698.3613 -11719.2953
[390,] -13768.4306 -9698.3613
[391,] -12369.5281 -13768.4306
[392,] -12827.7549 -12369.5281
[393,] -12034.4520 -12827.7549
[394,] -9948.6925 -12034.4520
[395,] -1745.4977 -9948.6925
[396,] -5871.5347 -1745.4977
[397,] -9481.0027 -5871.5347
[398,] -8563.1971 -9481.0027
[399,] -9195.3685 -8563.1971
[400,] 23310.0873 -9195.3685
[401,] -3322.0735 23310.0873
[402,] 2004.8449 -3322.0735
[403,] 3973.2002 2004.8449
[404,] 150298.5544 3973.2002
[405,] -24029.0406 150298.5544
[406,] 45980.9817 -24029.0406
[407,] 45618.1761 45980.9817
[408,] 39841.0064 45618.1761
[409,] 21399.4789 39841.0064
[410,] 19376.4440 21399.4789
[411,] 24599.6358 19376.4440
[412,] 15053.0460 24599.6358
[413,] -2429.6048 15053.0460
[414,] 4023.7643 -2429.6048
[415,] -3725.8022 4023.7643
[416,] -688.5390 -3725.8022
[417,] 21036.9928 -688.5390
[418,] 10968.0252 21036.9928
[419,] -11263.2979 10968.0252
[420,] -10675.1240 -11263.2979
[421,] 9121.2801 -10675.1240
[422,] 20540.6647 9121.2801
[423,] -1221.1366 20540.6647
[424,] -1079.6174 -1221.1366
[425,] 5009.0938 -1079.6174
[426,] -3003.4174 5009.0938
[427,] 1141.6817 -3003.4174
[428,] -6418.6668 1141.6817
[429,] 3050.9329 -6418.6668
[430,] 11984.6767 3050.9329
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 39723.7974 67597.3278
2 22636.8972 39723.7974
3 -14077.8071 22636.8972
4 -5704.0262 -14077.8071
5 -3969.4508 -5704.0262
6 -5468.7088 -3969.4508
7 -3171.3922 -5468.7088
8 -8007.7062 -3171.3922
9 -3125.1589 -8007.7062
10 -7744.1998 -3125.1589
11 -5891.4039 -7744.1998
12 -4401.0919 -5891.4039
13 -14175.6398 -4401.0919
14 -8909.2994 -14175.6398
15 -6690.3245 -8909.2994
16 -10436.7778 -6690.3245
17 -9343.1248 -10436.7778
18 -9182.1943 -9343.1248
19 -11682.5550 -9182.1943
20 -8262.8788 -11682.5550
21 -10036.2001 -8262.8788
22 -10003.2218 -10036.2001
23 -12271.3588 -10003.2218
24 -12014.9410 -12271.3588
25 -10332.9962 -12014.9410
26 -11163.7273 -10332.9962
27 -12357.7540 -11163.7273
28 -8075.8499 -12357.7540
29 -11345.4600 -8075.8499
30 -8604.2072 -11345.4600
31 -10504.8608 -8604.2072
32 -8691.7598 -10504.8608
33 -9663.8672 -8691.7598
34 -10272.0322 -9663.8672
35 -12412.0269 -10272.0322
36 -9907.2548 -12412.0269
37 -9277.2744 -9907.2548
38 -10552.1168 -9277.2744
39 -12171.3138 -10552.1168
40 -12597.3304 -12171.3138
41 -12146.0802 -12597.3304
42 -10506.8167 -12146.0802
43 -11570.2580 -10506.8167
44 -10959.5187 -11570.2580
45 -11066.6252 -10959.5187
46 -12199.4867 -11066.6252
47 -11574.0036 -12199.4867
48 -10864.9241 -11574.0036
49 -12095.7642 -10864.9241
50 -10258.6656 -12095.7642
51 -10641.4424 -10258.6656
52 -11424.6950 -10641.4424
53 -10438.7596 -11424.6950
54 -12652.7717 -10438.7596
55 -10674.8294 -12652.7717
56 -7931.6590 -10674.8294
57 -10804.9200 -7931.6590
58 -11280.4274 -10804.9200
59 -11554.5445 -11280.4274
60 -10513.2053 -11554.5445
61 -12281.2829 -10513.2053
62 -10800.0118 -12281.2829
63 -11520.8955 -10800.0118
64 -11011.7762 -11520.8955
65 -12479.8449 -11011.7762
66 -11622.4581 -12479.8449
67 -10937.8147 -11622.4581
68 -11301.7111 -10937.8147
69 -10934.6249 -11301.7111
70 -11809.5967 -10934.6249
71 -10878.5941 -11809.5967
72 -11469.1560 -10878.5941
73 -10903.3164 -11469.1560
74 -11025.8128 -10903.3164
75 -11705.7003 -11025.8128
76 -11141.3152 -11705.7003
77 -11715.7240 -11141.3152
78 -11282.1145 -11715.7240
79 -11430.9917 -11282.1145
80 -10864.7099 -11430.9917
81 -10934.9377 -10864.7099
82 -11424.5128 -10934.9377
83 -10945.9050 -11424.5128
84 -10741.8724 -10945.9050
85 -11286.5622 -10741.8724
86 -10658.9196 -11286.5622
87 -10917.5817 -10658.9196
88 -11596.1079 -10917.5817
89 -10569.9313 -11596.1079
90 -10936.9838 -10569.9313
91 -10998.5909 -10936.9838
92 -10418.7580 -10998.5909
93 -11752.1999 -10418.7580
94 -10745.7005 -11752.1999
95 -10872.2244 -10745.7005
96 -10060.0624 -10872.2244
97 -11244.0322 -10060.0624
98 -11015.0061 -11244.0322
99 -11011.1878 -11015.0061
100 -11277.6763 -11011.1878
101 -11841.9552 -11277.6763
102 -11796.1607 -11841.9552
103 -11362.0832 -11796.1607
104 -10871.1308 -11362.0832
105 -11442.5885 -10871.1308
106 -12577.4531 -11442.5885
107 -9927.2119 -12577.4531
108 -11777.8299 -9927.2119
109 -11718.1146 -11777.8299
110 -10481.7946 -11718.1146
111 -12296.7286 -10481.7946
112 -11826.3110 -12296.7286
113 -11387.2132 -11826.3110
114 -10309.3095 -11387.2132
115 -10941.0249 -10309.3095
116 -11078.0281 -10941.0249
117 -10405.1009 -11078.0281
118 -12002.8881 -10405.1009
119 -11156.8685 -12002.8881
120 -11928.4702 -11156.8685
121 -10865.3644 -11928.4702
122 -11213.2301 -10865.3644
123 -10311.3701 -11213.2301
124 -11013.4744 -10311.3701
125 -12071.7149 -11013.4744
126 -11097.0398 -12071.7149
127 -32785.7374 -11097.0398
128 -31684.9904 -32785.7374
129 27735.0273 -31684.9904
130 17411.7203 27735.0273
131 17904.8315 17411.7203
132 -38256.9643 17904.8315
133 -16557.0026 -38256.9643
134 -14364.6072 -16557.0026
135 -13816.9600 -14364.6072
136 -13621.1866 -13816.9600
137 -27313.8478 -13621.1866
138 -14497.6276 -27313.8478
139 -4995.3696 -14497.6276
140 -4319.3008 -4995.3696
141 -13709.3539 -4319.3008
142 -12558.8073 -13709.3539
143 -13509.6039 -12558.8073
144 -13759.2662 -13509.6039
145 -13939.2648 -13759.2662
146 -146690.8873 -13939.2648
147 -151273.1997 -146690.8873
148 -147761.4556 -151273.1997
149 -136876.5586 -147761.4556
150 -142916.9026 -136876.5586
151 150666.9555 -142916.9026
152 298617.3353 150666.9555
153 286743.2309 298617.3353
154 -3397.6479 286743.2309
155 -9673.8952 -3397.6479
156 -10523.2124 -9673.8952
157 -11103.0672 -10523.2124
158 27824.0100 -11103.0672
159 -38383.9976 27824.0100
160 -5515.8671 -38383.9976
161 -5563.1901 -5515.8671
162 -3116.3550 -5563.1901
163 -4564.9795 -3116.3550
164 -9772.4107 -4564.9795
165 -14175.2629 -9772.4107
166 -13946.5282 -14175.2629
167 165295.8507 -13946.5282
168 290928.2368 165295.8507
169 8665.4050 290928.2368
170 -3424.0406 8665.4050
171 -10833.1352 -3424.0406
172 -11189.8796 -10833.1352
173 -8193.6683 -11189.8796
174 -11103.1552 -8193.6683
175 38768.3065 -11103.1552
176 29689.7823 38768.3065
177 16382.8778 29689.7823
178 33469.1296 16382.8778
179 25493.1049 33469.1296
180 22203.2177 25493.1049
181 25154.5139 22203.2177
182 20676.3832 25154.5139
183 27709.4339 20676.3832
184 24352.3221 27709.4339
185 -38496.4683 24352.3221
186 -15178.2771 -38496.4683
187 -4652.7249 -15178.2771
188 -5294.1555 -4652.7249
189 -12570.6859 -5294.1555
190 -13732.3666 -12570.6859
191 -189742.3901 -13732.3666
192 168268.5537 -189742.3901
193 -642.2560 168268.5537
194 -3392.1365 -642.2560
195 -11102.2531 -3392.1365
196 17912.8496 -11102.2531
197 50077.8428 17912.8496
198 22480.7138 50077.8428
199 27926.3424 22480.7138
200 22415.2256 27926.3424
201 -38968.1121 22415.2256
202 -12785.8360 -38968.1121
203 -3827.1839 -12785.8360
204 -11072.6088 -3827.1839
205 -31365.8697 -11072.6088
206 -33763.8745 -31365.8697
207 -33738.4296 -33763.8745
208 27935.4232 -33738.4296
209 22316.9405 27935.4232
210 22120.2418 22316.9405
211 -38625.0343 22120.2418
212 -12664.3209 -38625.0343
213 -14668.8941 -12664.3209
214 30883.9432 -14668.8941
215 -5369.7218 30883.9432
216 297.7586 -5369.7218
217 -5048.6971 297.7586
218 -3530.7233 -5048.6971
219 -5129.4265 -3530.7233
220 -13226.1934 -5129.4265
221 -13931.4692 -13226.1934
222 -140269.1905 -13931.4692
223 -149337.3644 -140269.1905
224 156115.8886 -149337.3644
225 299908.8713 156115.8886
226 112093.8766 299908.8713
227 -7648.3683 112093.8766
228 -10631.9757 -7648.3683
229 -11103.4633 -10631.9757
230 27824.1353 -11103.4633
231 -38249.6540 27824.1353
232 -16084.1128 -38249.6540
233 -635.4399 -16084.1128
234 -6115.0196 -635.4399
235 -5320.8986 -6115.0196
236 -13928.0376 -5320.8986
237 -144329.2619 -13928.0376
238 168268.5683 -144329.2619
239 6711.1829 168268.5683
240 -3418.3413 6711.1829
241 -9843.5010 -3418.3413
242 -11298.9017 -9843.5010
243 -11849.3718 -11298.9017
244 -11096.4140 -11849.3718
245 -30441.0379 -11096.4140
246 -27792.0379 -30441.0379
247 27823.8222 -27792.0379
248 -38351.4679 27823.8222
249 -5093.7170 -38351.4679
250 -13387.5614 -5093.7170
251 -13890.6039 -13387.5614
252 -13945.9322 -13890.6039
253 168268.8762 -13945.9322
254 15098.3623 168268.8762
255 -3458.6759 15098.3623
256 -10536.5767 -3458.6759
257 -11101.4609 -10536.5767
258 27823.5021 -11101.4609
259 -38265.2582 27823.5021
260 -13583.6840 -38265.2582
261 -4921.4852 -13583.6840
262 -4988.3175 -4921.4852
263 -3996.9991 -4988.3175
264 -13509.4897 -3996.9991
265 -13706.9800 -13509.4897
266 -13946.1309 -13706.9800
267 166404.0359 -13946.1309
268 301756.3289 166404.0359
269 7330.5145 301756.3289
270 -3412.2202 7330.5145
271 -10633.8450 -3412.2202
272 -10378.9590 -10633.8450
273 -11041.4591 -10378.9590
274 -30591.8150 -11041.4591
275 29948.8615 -30591.8150
276 41948.5409 29948.8615
277 18262.8119 41948.5409
278 -38229.6500 18262.8119
279 -14207.2020 -38229.6500
280 -13226.3638 -14207.2020
281 -5321.0494 -13226.3638
282 -13937.3295 -5321.0494
283 166051.9954 -13937.3295
284 276657.2340 166051.9954
285 281675.7833 276657.2340
286 -3398.3637 281675.7833
287 -10917.7235 -3398.3637
288 -11112.9927 -10917.7235
289 -11090.9554 -11112.9927
290 -34140.8152 -11090.9554
291 -30469.1148 -34140.8152
292 27823.9057 -30469.1148
293 -38016.2974 27823.9057
294 -12838.7941 -38016.2974
295 -2755.3435 -12838.7941
296 32452.2003 -2755.3435
297 -3815.9283 32452.2003
298 -2340.5326 -3815.9283
299 -5632.4148 -2340.5326
300 78793.0612 -5632.4148
301 -5319.9433 78793.0612
302 -13946.1929 -5319.9433
303 168267.3517 -13946.1929
304 -3395.1691 168267.3517
305 -11099.8613 -3395.1691
306 -32477.2195 -11099.8613
307 27823.7317 -32477.2195
308 -38366.2030 27823.7317
309 -5319.8678 -38366.2030
310 -13945.9943 -5319.8678
311 168268.7882 -13945.9943
312 19295.4177 168268.7882
313 54217.6282 19295.4177
314 10080.3950 54217.6282
315 -2644.9526 10080.3950
316 -3389.0318 -2644.9526
317 -8273.5090 -3389.0318
318 -10875.8196 -8273.5090
319 -11110.3020 -10875.8196
320 -11103.6393 -11110.3020
321 27824.2118 -11103.6393
322 -38217.4694 27824.2118
323 -12830.3363 -38217.4694
324 -4851.8719 -12830.3363
325 -5322.3566 -4851.8719
326 -13896.6187 -5322.3566
327 -148469.3140 -13896.6187
328 -148028.3140 -148469.3140
329 -146325.5946 -148028.3140
330 -151000.3286 -146325.5946
331 -146817.7818 -151000.3286
332 168268.7442 -146817.7818
333 -467.8752 168268.7442
334 -3399.6751 -467.8752
335 -11102.2311 -3399.6751
336 27825.7082 -11102.2311
337 21467.4703 27825.7082
338 -38337.9865 21467.4703
339 -6152.3353 -38337.9865
340 -4789.0893 -6152.3353
341 -5319.6165 -4789.0893
342 -13946.1805 -5319.6165
343 168269.0228 -13946.1805
344 2106.9224 168269.0228
345 -3423.6655 2106.9224
346 -11025.5434 -3423.6655
347 -10727.7474 -11025.5434
348 -8495.3349 -10727.7474
349 -11101.9010 -8495.3349
350 23284.7455 -11101.9010
351 23577.8067 23284.7455
352 21902.5861 23577.8067
353 22016.0832 21902.5861
354 -38364.2154 22016.0832
355 -5207.4807 -38364.2154
356 -6303.2821 -5207.4807
357 -5242.1192 -6303.2821
358 -11249.1952 -5242.1192
359 -13166.5866 -11249.1952
360 -12632.1763 -13166.5866
361 -13560.5184 -12632.1763
362 -12494.1668 -13560.5184
363 -13379.5186 -12494.1668
364 -13244.0418 -13379.5186
365 -12926.6095 -13244.0418
366 -12687.7083 -12926.6095
367 -10314.2730 -12687.7083
368 -12678.8124 -10314.2730
369 -8699.3621 -12678.8124
370 -12110.9923 -8699.3621
371 -13447.8481 -12110.9923
372 -11403.8518 -13447.8481
373 -12899.4066 -11403.8518
374 -12490.3652 -12899.4066
375 -13083.9550 -12490.3652
376 -12045.9084 -13083.9550
377 -12913.3510 -12045.9084
378 -12201.9231 -12913.3510
379 -13048.5309 -12201.9231
380 -12702.6642 -13048.5309
381 -13252.4426 -12702.6642
382 -12937.0619 -13252.4426
383 -11691.9427 -12937.0619
384 -6337.8772 -11691.9427
385 -12491.0807 -6337.8772
386 -10283.3339 -12491.0807
387 -12587.6097 -10283.3339
388 -11719.2953 -12587.6097
389 -9698.3613 -11719.2953
390 -13768.4306 -9698.3613
391 -12369.5281 -13768.4306
392 -12827.7549 -12369.5281
393 -12034.4520 -12827.7549
394 -9948.6925 -12034.4520
395 -1745.4977 -9948.6925
396 -5871.5347 -1745.4977
397 -9481.0027 -5871.5347
398 -8563.1971 -9481.0027
399 -9195.3685 -8563.1971
400 23310.0873 -9195.3685
401 -3322.0735 23310.0873
402 2004.8449 -3322.0735
403 3973.2002 2004.8449
404 150298.5544 3973.2002
405 -24029.0406 150298.5544
406 45980.9817 -24029.0406
407 45618.1761 45980.9817
408 39841.0064 45618.1761
409 21399.4789 39841.0064
410 19376.4440 21399.4789
411 24599.6358 19376.4440
412 15053.0460 24599.6358
413 -2429.6048 15053.0460
414 4023.7643 -2429.6048
415 -3725.8022 4023.7643
416 -688.5390 -3725.8022
417 21036.9928 -688.5390
418 10968.0252 21036.9928
419 -11263.2979 10968.0252
420 -10675.1240 -11263.2979
421 9121.2801 -10675.1240
422 20540.6647 9121.2801
423 -1221.1366 20540.6647
424 -1079.6174 -1221.1366
425 5009.0938 -1079.6174
426 -3003.4174 5009.0938
427 1141.6817 -3003.4174
428 -6418.6668 1141.6817
429 3050.9329 -6418.6668
430 11984.6767 3050.9329
> 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/7eb3m1291301535.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/862k71291301535.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/962k71291301535.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/10zcks1291301535.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/113uig1291301535.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/126dz41291301535.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/1324ec1291301535.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/1465d01291301535.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/1595c61291301535.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/16u6sc1291301535.tab")
+ }
>
> try(system("convert tmp/1sb5g1291301535.ps tmp/1sb5g1291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/2sb5g1291301535.ps tmp/2sb5g1291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/3sb5g1291301535.ps tmp/3sb5g1291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/432m11291301535.ps tmp/432m11291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/532m11291301535.ps tmp/532m11291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/632m11291301535.ps tmp/632m11291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/7eb3m1291301535.ps tmp/7eb3m1291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/862k71291301535.ps tmp/862k71291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/962k71291301535.ps tmp/962k71291301535.png",intern=TRUE))
character(0)
> try(system("convert tmp/10zcks1291301535.ps tmp/10zcks1291301535.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.333 2.210 32.513