R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,162556
+ ,1081
+ ,807
+ ,213118
+ ,230380558
+ ,6282154
+ ,1
+ ,29790
+ ,309
+ ,444
+ ,81767
+ ,25266003
+ ,4321023
+ ,1
+ ,87550
+ ,458
+ ,412
+ ,153198
+ ,70164684
+ ,4111912
+ ,0
+ ,84738
+ ,588
+ ,428
+ ,-26007
+ ,-15292116
+ ,223193
+ ,1
+ ,54660
+ ,302
+ ,315
+ ,126942
+ ,37955658
+ ,1491348
+ ,1
+ ,42634
+ ,156
+ ,168
+ ,157214
+ ,24525384
+ ,1629616
+ ,0
+ ,40949
+ ,481
+ ,263
+ ,129352
+ ,62218312
+ ,1398893
+ ,1
+ ,45187
+ ,353
+ ,267
+ ,234817
+ ,75845891
+ ,1926517
+ ,1
+ ,37704
+ ,452
+ ,228
+ ,60448
+ ,27322496
+ ,983660
+ ,1
+ ,16275
+ ,109
+ ,129
+ ,47818
+ ,5212162
+ ,1443586
+ ,0
+ ,25830
+ ,115
+ ,104
+ ,245546
+ ,28237790
+ ,1073089
+ ,0
+ ,12679
+ ,110
+ ,122
+ ,48020
+ ,5282200
+ ,984885
+ ,1
+ ,18014
+ ,239
+ ,393
+ ,-1710
+ ,-408690
+ ,1405225
+ ,0
+ ,43556
+ ,247
+ ,190
+ ,32648
+ ,8064056
+ ,227132
+ ,1
+ ,24811
+ ,505
+ ,280
+ ,95350
+ ,47388950
+ ,929118
+ ,0
+ ,6575
+ ,159
+ ,63
+ ,151352
+ ,15589256
+ ,1071292
+ ,0
+ ,7123
+ ,109
+ ,102
+ ,288170
+ ,31410530
+ ,638830
+ ,1
+ ,21950
+ ,519
+ ,265
+ ,114337
+ ,57397174
+ ,856956
+ ,1
+ ,37597
+ ,248
+ ,234
+ ,37884
+ ,9395232
+ ,992426
+ ,0
+ ,17821
+ ,373
+ ,277
+ ,122844
+ ,45820812
+ ,444477
+ ,1
+ ,12988
+ ,119
+ ,73
+ ,82340
+ ,9798460
+ ,857217
+ ,1
+ ,22330
+ ,84
+ ,67
+ ,79801
+ ,6703284
+ ,711969
+ ,0
+ ,13326
+ ,102
+ ,103
+ ,165548
+ ,16885896
+ ,702380
+ ,0
+ ,16189
+ ,295
+ ,290
+ ,116384
+ ,34333280
+ ,358589
+ ,0
+ ,7146
+ ,105
+ ,83
+ ,134028
+ ,14072940
+ ,297978
+ ,0
+ ,15824
+ ,64
+ ,56
+ ,63838
+ ,4085632
+ ,585715
+ ,1
+ ,27664
+ ,282
+ ,236
+ ,74996
+ ,20023932
+ ,657954
+ ,0
+ ,11920
+ ,182
+ ,73
+ ,31080
+ ,4009320
+ ,209458
+ ,0
+ ,8568
+ ,37
+ ,34
+ ,32168
+ ,1190216
+ ,786690
+ ,0
+ ,14416
+ ,361
+ ,139
+ ,49857
+ ,17998377
+ ,439798
+ ,1
+ ,3369
+ ,28
+ ,26
+ ,87161
+ ,2440508
+ ,688779
+ ,1
+ ,11819
+ ,85
+ ,70
+ ,106113
+ ,9019605
+ ,574339
+ ,1
+ ,6984
+ ,45
+ ,40
+ ,80570
+ ,3545080
+ ,741409
+ ,1
+ ,4519
+ ,49
+ ,42
+ ,102129
+ ,5004321
+ ,597793
+ ,0
+ ,2220
+ ,22
+ ,12
+ ,301670
+ ,6636740
+ ,644190
+ ,0
+ ,18562
+ ,155
+ ,211
+ ,102313
+ ,15858515
+ ,377934
+ ,0
+ ,10327
+ ,91
+ ,74
+ ,88577
+ ,8060507
+ ,640273
+ ,1
+ ,5336
+ ,81
+ ,80
+ ,112477
+ ,9110637
+ ,697458
+ ,1
+ ,2365
+ ,79
+ ,83
+ ,191778
+ ,15150462
+ ,550608
+ ,0
+ ,4069
+ ,145
+ ,131
+ ,79804
+ ,11571580
+ ,207393
+ ,0
+ ,8636
+ ,855
+ ,203
+ ,128294
+ ,104687904
+ ,301607
+ ,0
+ ,13718
+ ,61
+ ,56
+ ,96448
+ ,5883328
+ ,345783
+ ,0
+ ,4525
+ ,226
+ ,89
+ ,93811
+ ,21201286
+ ,501749
+ ,0
+ ,6869
+ ,105
+ ,88
+ ,117520
+ ,12339600
+ ,379983
+ ,0
+ ,4628
+ ,62
+ ,39
+ ,69159
+ ,4287858
+ ,387475
+ ,1
+ ,3689
+ ,25
+ ,25
+ ,101792
+ ,2443008
+ ,377305
+ ,1
+ ,4891
+ ,217
+ ,49
+ ,210568
+ ,5474768
+ ,370837
+ ,1
+ ,7489
+ ,322
+ ,149
+ ,136996
+ ,44112712
+ ,430866
+ ,0
+ ,4901
+ ,84
+ ,58
+ ,121920
+ ,10241280
+ ,469107
+ ,0
+ ,2284
+ ,33
+ ,41
+ ,76403
+ ,2521299
+ ,194493
+ ,1
+ ,3160
+ ,108
+ ,90
+ ,108094
+ ,11674152
+ ,530670
+ ,1
+ ,4150
+ ,150
+ ,136
+ ,134759
+ ,20213850
+ ,518365
+ ,1
+ ,7285
+ ,115
+ ,97
+ ,188873
+ ,21720395
+ ,491303
+ ,1
+ ,1134
+ ,162
+ ,63
+ ,146216
+ ,23686992
+ ,527021
+ ,1
+ ,4658
+ ,158
+ ,114
+ ,156608
+ ,24744064
+ ,233773
+ ,0
+ ,2384
+ ,97
+ ,77
+ ,61348
+ ,5950756
+ ,405972
+ ,0
+ ,3748
+ ,9
+ ,6
+ ,50350
+ ,453150
+ ,652925
+ ,0
+ ,5371
+ ,66
+ ,47
+ ,87720
+ ,5789520
+ ,446211
+ ,0
+ ,1285
+ ,107
+ ,51
+ ,99489
+ ,10645323
+ ,341340
+ ,1
+ ,9327
+ ,101
+ ,85
+ ,87419
+ ,8829319
+ ,387699
+ ,1
+ ,5565
+ ,47
+ ,43
+ ,94355
+ ,4434685
+ ,493408
+ ,0
+ ,1528
+ ,38
+ ,32
+ ,60326
+ ,2292388
+ ,146494
+ ,1
+ ,3122
+ ,34
+ ,25
+ ,94670
+ ,3218780
+ ,414462
+ ,1
+ ,7561
+ ,87
+ ,77
+ ,82425
+ ,6923700
+ ,364304
+ ,0
+ ,2675
+ ,79
+ ,54
+ ,59017
+ ,4662343
+ ,355178
+ ,0
+ ,13253
+ ,947
+ ,251
+ ,90829
+ ,86015063
+ ,357760
+ ,0
+ ,880
+ ,74
+ ,15
+ ,80791
+ ,5978534
+ ,261216
+ ,1
+ ,2053
+ ,53
+ ,44
+ ,100423
+ ,5322419
+ ,397144
+ ,0
+ ,1424
+ ,94
+ ,73
+ ,131116
+ ,12324904
+ ,374943
+ ,1
+ ,4036
+ ,63
+ ,85
+ ,100269
+ ,6316947
+ ,424898
+ ,1
+ ,3045
+ ,58
+ ,49
+ ,27330
+ ,1585140
+ ,202055
+ ,0
+ ,5119
+ ,49
+ ,38
+ ,39039
+ ,1912911
+ ,378525
+ ,0
+ ,1431
+ ,34
+ ,35
+ ,106885
+ ,3634090
+ ,310768
+ ,0
+ ,554
+ ,11
+ ,9
+ ,79285
+ ,872135
+ ,325738
+ ,0
+ ,1975
+ ,35
+ ,34
+ ,118881
+ ,4160835
+ ,394510
+ ,1
+ ,1765
+ ,20
+ ,20
+ ,77623
+ ,1319591
+ ,247060
+ ,0
+ ,1012
+ ,47
+ ,29
+ ,114768
+ ,5394096
+ ,368078
+ ,0
+ ,810
+ ,43
+ ,11
+ ,74015
+ ,3182645
+ ,236761
+ ,0
+ ,1280
+ ,117
+ ,52
+ ,69465
+ ,8127405
+ ,312378
+ ,1
+ ,666
+ ,171
+ ,13
+ ,117869
+ ,20155599
+ ,339836
+ ,0
+ ,1380
+ ,26
+ ,29
+ ,60982
+ ,1585532
+ ,347385
+ ,1
+ ,4677
+ ,75
+ ,66
+ ,90131
+ ,6579563
+ ,426280
+ ,0
+ ,876
+ ,59
+ ,33
+ ,138971
+ ,8199289
+ ,352850
+ ,0
+ ,814
+ ,18
+ ,15
+ ,39625
+ ,713250
+ ,301881
+ ,0
+ ,514
+ ,15
+ ,15
+ ,102725
+ ,1540875
+ ,377516
+ ,1
+ ,5692
+ ,72
+ ,68
+ ,64239
+ ,4625208
+ ,357312
+ ,0
+ ,3642
+ ,86
+ ,100
+ ,90262
+ ,7762532
+ ,458343
+ ,0
+ ,540
+ ,14
+ ,13
+ ,103960
+ ,1455440
+ ,354228
+ ,0
+ ,2099
+ ,64
+ ,45
+ ,106611
+ ,6823104
+ ,308636
+ ,0
+ ,567
+ ,11
+ ,14
+ ,103345
+ ,1136795
+ ,386212
+ ,0
+ ,2001
+ ,52
+ ,36
+ ,95551
+ ,4968652
+ ,393343
+ ,1
+ ,2949
+ ,41
+ ,40
+ ,82903
+ ,3399023
+ ,378509
+ ,0
+ ,2253
+ ,99
+ ,68
+ ,63593
+ ,6295707
+ ,452469
+ ,1
+ ,6533
+ ,75
+ ,29
+ ,126910
+ ,9518250
+ ,364839
+ ,0
+ ,1889
+ ,45
+ ,43
+ ,37527
+ ,1688715
+ ,358649
+ ,1
+ ,3055
+ ,43
+ ,30
+ ,60247
+ ,2590621
+ ,376641
+ ,0
+ ,272
+ ,8
+ ,9
+ ,112995
+ ,903960
+ ,429112
+ ,1
+ ,1414
+ ,198
+ ,22
+ ,70184
+ ,13896432
+ ,330546
+ ,0
+ ,2564
+ ,22
+ ,19
+ ,130140
+ ,2863080
+ ,403560
+ ,1
+ ,1383
+ ,11
+ ,9
+ ,73221
+ ,805431
+ ,317892
+ ,1
+ ,1261
+ ,33
+ ,31
+ ,76114
+ ,2511762
+ ,307528
+ ,0
+ ,975
+ ,23
+ ,19
+ ,90534
+ ,2082282
+ ,235133
+ ,0
+ ,3366
+ ,80
+ ,55
+ ,108479
+ ,8678320
+ ,299243
+ ,0
+ ,576
+ ,18
+ ,8
+ ,113761
+ ,2047698
+ ,314073
+ ,0
+ ,1686
+ ,40
+ ,28
+ ,68696
+ ,1923488
+ ,368186
+ ,0
+ ,746
+ ,23
+ ,29
+ ,71561
+ ,1645903
+ ,269661
+ ,1
+ ,3192
+ ,60
+ ,48
+ ,59831
+ ,3589860
+ ,125390
+ ,1
+ ,2045
+ ,20
+ ,16
+ ,97890
+ ,1957800
+ ,510834
+ ,0
+ ,5702
+ ,61
+ ,47
+ ,101481
+ ,5987379
+ ,321896
+ ,1
+ ,1932
+ ,36
+ ,20
+ ,72954
+ ,2626344
+ ,249898
+ ,0
+ ,936
+ ,30
+ ,22
+ ,67939
+ ,2038170
+ ,408881
+ ,1
+ ,3437
+ ,47
+ ,33
+ ,48022
+ ,2257034
+ ,158492
+ ,0
+ ,5131
+ ,71
+ ,44
+ ,86111
+ ,6113881
+ ,292154
+ ,1
+ ,2397
+ ,14
+ ,13
+ ,74020
+ ,1036280
+ ,289513
+ ,1
+ ,1389
+ ,9
+ ,6
+ ,57530
+ ,517770
+ ,378049
+ ,0
+ ,1503
+ ,39
+ ,35
+ ,56364
+ ,2198196
+ ,343466
+ ,0
+ ,402
+ ,26
+ ,8
+ ,84990
+ ,2209740
+ ,332743
+ ,0
+ ,2239
+ ,21
+ ,17
+ ,88590
+ ,1860390
+ ,442882
+ ,1
+ ,2234
+ ,16
+ ,11
+ ,77200
+ ,1235200
+ ,214215
+ ,0
+ ,837
+ ,69
+ ,21
+ ,61262
+ ,4227078
+ ,315688
+ ,0
+ ,10579
+ ,92
+ ,92
+ ,110309
+ ,10148428
+ ,375195
+ ,0
+ ,875
+ ,14
+ ,12
+ ,67000
+ ,938000
+ ,334280
+ ,0
+ ,1585
+ ,107
+ ,112
+ ,93099
+ ,9589197
+ ,355864
+ ,1
+ ,1659
+ ,29
+ ,25
+ ,107577
+ ,3119733
+ ,480382
+ ,1
+ ,2647
+ ,37
+ ,17
+ ,62920
+ ,2328040
+ ,353058
+ ,1
+ ,3294
+ ,23
+ ,23
+ ,75832
+ ,1744136
+ ,217193
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,94
+ ,7
+ ,10
+ ,60793
+ ,425551
+ ,314533
+ ,0
+ ,422
+ ,28
+ ,23
+ ,57935
+ ,1622180
+ ,318056
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,34
+ ,8
+ ,7
+ ,60630
+ ,485040
+ ,314353
+ ,0
+ ,1558
+ ,63
+ ,25
+ ,55637
+ ,3505131
+ ,369448
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,43
+ ,3
+ ,20
+ ,60887
+ ,182661
+ ,312846
+ ,0
+ ,645
+ ,5
+ ,4
+ ,60720
+ ,0
+ ,312075
+ ,0
+ ,316
+ ,9
+ ,4
+ ,60505
+ ,544545
+ ,315009
+ ,0
+ ,115
+ ,13
+ ,10
+ ,60945
+ ,792285
+ ,318903
+ ,0
+ ,5
+ ,2
+ ,1
+ ,60720
+ ,0
+ ,314887
+ ,0
+ ,897
+ ,5
+ ,4
+ ,60720
+ ,0
+ ,314913
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,389
+ ,14
+ ,8
+ ,58990
+ ,825860
+ ,325506
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1002
+ ,15
+ ,11
+ ,56750
+ ,851250
+ ,298568
+ ,0
+ ,36
+ ,3
+ ,4
+ ,60894
+ ,182682
+ ,315834
+ ,0
+ ,460
+ ,15
+ ,15
+ ,63346
+ ,950190
+ ,329784
+ ,0
+ ,309
+ ,11
+ ,9
+ ,56535
+ ,621885
+ ,312878
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,9
+ ,6
+ ,7
+ ,60835
+ ,365010
+ ,314987
+ ,0
+ ,271
+ ,2
+ ,2
+ ,60720
+ ,0
+ ,325249
+ ,0
+ ,14
+ ,1
+ ,0
+ ,61016
+ ,61016
+ ,315877
+ ,0
+ ,520
+ ,10
+ ,7
+ ,58650
+ ,586500
+ ,291650
+ ,0
+ ,1766
+ ,73
+ ,46
+ ,60438
+ ,4411974
+ ,305959
+ ,0
+ ,0
+ ,0
+ ,5
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,458
+ ,11
+ ,7
+ ,58625
+ ,644875
+ ,297765
+ ,0
+ ,20
+ ,3
+ ,2
+ ,60938
+ ,182814
+ ,315245
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,98
+ ,2
+ ,2
+ ,61490
+ ,122980
+ ,315236
+ ,0
+ ,405
+ ,7
+ ,5
+ ,60845
+ ,425915
+ ,336425
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,483
+ ,27
+ ,7
+ ,60830
+ ,1642410
+ ,306268
+ ,0
+ ,454
+ ,51
+ ,24
+ ,63261
+ ,3226311
+ ,302187
+ ,0
+ ,47
+ ,3
+ ,1
+ ,60720
+ ,0
+ ,314882
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,757
+ ,19
+ ,18
+ ,45689
+ ,868091
+ ,382712
+ ,0
+ ,4655
+ ,393
+ ,55
+ ,60720
+ ,0
+ ,341570
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,36
+ ,4
+ ,3
+ ,61564
+ ,246256
+ ,312412
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,203
+ ,9
+ ,9
+ ,61938
+ ,557442
+ ,309596
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,126
+ ,10
+ ,8
+ ,60951
+ ,487608
+ ,315547
+ ,0
+ ,400
+ ,152
+ ,113
+ ,60720
+ ,0
+ ,313267
+ ,1
+ ,71
+ ,1
+ ,0
+ ,60745
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,972
+ ,34
+ ,19
+ ,71642
+ ,2435828
+ ,359335
+ ,1
+ ,531
+ ,10
+ ,11
+ ,71641
+ ,716410
+ ,330068
+ ,0
+ ,2461
+ ,57
+ ,25
+ ,55792
+ ,2120096
+ ,314289
+ ,1
+ ,378
+ ,52
+ ,16
+ ,71873
+ ,718730
+ ,297413
+ ,1
+ ,23
+ ,5
+ ,5
+ ,62555
+ ,312775
+ ,314806
+ ,1
+ ,638
+ ,14
+ ,11
+ ,60370
+ ,845180
+ ,333210
+ ,1
+ ,2300
+ ,29
+ ,23
+ ,64873
+ ,1037968
+ ,352108
+ ,0
+ ,149
+ ,5
+ ,6
+ ,62041
+ ,310205
+ ,313332
+ ,0
+ ,226
+ ,5
+ ,5
+ ,65745
+ ,328725
+ ,291787
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,275
+ ,4
+ ,7
+ ,59500
+ ,238000
+ ,318745
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,141
+ ,6
+ ,7
+ ,61630
+ ,369780
+ ,315366
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,28
+ ,2
+ ,3
+ ,60890
+ ,121780
+ ,315688
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,4980
+ ,91
+ ,89
+ ,113521
+ ,9081680
+ ,409642
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,472
+ ,20
+ ,19
+ ,80045
+ ,1600900
+ ,269587
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,203
+ ,27
+ ,12
+ ,50804
+ ,1371708
+ ,300962
+ ,1
+ ,496
+ ,17
+ ,12
+ ,87390
+ ,1485630
+ ,325479
+ ,1
+ ,10
+ ,2
+ ,5
+ ,61656
+ ,123312
+ ,316155
+ ,1
+ ,63
+ ,4
+ ,2
+ ,65688
+ ,262752
+ ,318574
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,1136
+ ,32
+ ,26
+ ,48522
+ ,1552704
+ ,343613
+ ,1
+ ,265
+ ,31
+ ,3
+ ,60720
+ ,0
+ ,306948
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,267
+ ,32
+ ,11
+ ,57640
+ ,1844480
+ ,330059
+ ,0
+ ,474
+ ,20
+ ,10
+ ,61977
+ ,1239540
+ ,288985
+ ,0
+ ,534
+ ,7
+ ,5
+ ,62620
+ ,438340
+ ,304485
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,15
+ ,8
+ ,6
+ ,60831
+ ,486648
+ ,315688
+ ,0
+ ,397
+ ,28
+ ,7
+ ,60646
+ ,1698088
+ ,317736
+ ,0
+ ,0
+ ,0
+ ,2
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1866
+ ,29
+ ,28
+ ,56225
+ ,1124500
+ ,322331
+ ,0
+ ,288
+ ,4
+ ,3
+ ,60510
+ ,242040
+ ,296656
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,3
+ ,2
+ ,1
+ ,60698
+ ,121396
+ ,315354
+ ,0
+ ,468
+ ,21
+ ,20
+ ,60720
+ ,0
+ ,312161
+ ,0
+ ,20
+ ,2
+ ,1
+ ,60805
+ ,121610
+ ,315576
+ ,0
+ ,278
+ ,26
+ ,22
+ ,61404
+ ,1596504
+ ,314922
+ ,0
+ ,61
+ ,14
+ ,9
+ ,60720
+ ,0
+ ,314551
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,192
+ ,4
+ ,2
+ ,65276
+ ,261104
+ ,312339
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,317
+ ,9
+ ,7
+ ,63915
+ ,575235
+ ,298700
+ ,0
+ ,738
+ ,10
+ ,9
+ ,60720
+ ,0
+ ,321376
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,368
+ ,17
+ ,13
+ ,61686
+ ,1048662
+ ,303230
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,2
+ ,1
+ ,0
+ ,60743
+ ,60743
+ ,315487
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,53
+ ,6
+ ,6
+ ,60349
+ ,362094
+ ,315793
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,94
+ ,3
+ ,3
+ ,61360
+ ,184080
+ ,312887
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,24
+ ,8
+ ,7
+ ,59818
+ ,478544
+ ,315637
+ ,0
+ ,2332
+ ,4
+ ,2
+ ,72680
+ ,290720
+ ,324385
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,131
+ ,11
+ ,15
+ ,61808
+ ,679888
+ ,308989
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,206
+ ,9
+ ,9
+ ,53110
+ ,477990
+ ,296702
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,167
+ ,2
+ ,1
+ ,64245
+ ,128490
+ ,307322
+ ,0
+ ,622
+ ,73
+ ,38
+ ,73007
+ ,5329511
+ ,304376
+ ,0
+ ,2328
+ ,94
+ ,57
+ ,82732
+ ,7032220
+ ,253588
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,365
+ ,8
+ ,7
+ ,54820
+ ,438560
+ ,309560
+ ,0
+ ,364
+ ,35
+ ,26
+ ,47705
+ ,1669675
+ ,298466
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,226
+ ,12
+ ,13
+ ,72835
+ ,874020
+ ,343929
+ ,0
+ ,307
+ ,15
+ ,10
+ ,58856
+ ,882840
+ ,331955
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,188
+ ,11
+ ,9
+ ,77655
+ ,854205
+ ,381180
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,138
+ ,6
+ ,26
+ ,69817
+ ,418902
+ ,331420
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,125
+ ,12
+ ,19
+ ,60798
+ ,729576
+ ,310201
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,282
+ ,30
+ ,12
+ ,62452
+ ,1873560
+ ,320016
+ ,0
+ ,335
+ ,33
+ ,23
+ ,64175
+ ,2117775
+ ,320398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,1324
+ ,117
+ ,29
+ ,67440
+ ,5530080
+ ,291841
+ ,0
+ ,176
+ ,28
+ ,8
+ ,68136
+ ,1907808
+ ,310670
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,249
+ ,72
+ ,26
+ ,56726
+ ,4084272
+ ,313491
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,333
+ ,13
+ ,9
+ ,70811
+ ,920543
+ ,331323
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,601
+ ,6
+ ,5
+ ,60720
+ ,0
+ ,319210
+ ,0
+ ,30
+ ,4
+ ,3
+ ,62045
+ ,248180
+ ,318098
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,249
+ ,62
+ ,13
+ ,54323
+ ,3368026
+ ,292754
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,165
+ ,24
+ ,12
+ ,62841
+ ,1508184
+ ,325176
+ ,0
+ ,453
+ ,21
+ ,19
+ ,81125
+ ,1703625
+ ,365959
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,53
+ ,14
+ ,10
+ ,59506
+ ,833084
+ ,302409
+ ,1
+ ,382
+ ,21
+ ,9
+ ,59365
+ ,1246665
+ ,340968
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,9
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,30
+ ,4
+ ,4
+ ,60798
+ ,243192
+ ,313164
+ ,0
+ ,290
+ ,2
+ ,1
+ ,58790
+ ,117580
+ ,301164
+ ,0
+ ,0
+ ,0
+ ,1
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,366
+ ,53
+ ,14
+ ,61808
+ ,3275824
+ ,344425
+ ,0
+ ,2
+ ,9
+ ,12
+ ,60735
+ ,546615
+ ,315394
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,209
+ ,13
+ ,19
+ ,64016
+ ,832208
+ ,316647
+ ,0
+ ,384
+ ,22
+ ,17
+ ,54683
+ ,1203026
+ ,309836
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,365
+ ,83
+ ,32
+ ,87192
+ ,7236936
+ ,346611
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,49
+ ,8
+ ,14
+ ,64107
+ ,512856
+ ,322031
+ ,0
+ ,3
+ ,4
+ ,8
+ ,60761
+ ,243044
+ ,315656
+ ,0
+ ,133
+ ,14
+ ,4
+ ,65990
+ ,923860
+ ,339445
+ ,0
+ ,32
+ ,1
+ ,0
+ ,59988
+ ,59988
+ ,314964
+ ,0
+ ,368
+ ,17
+ ,20
+ ,61167
+ ,1039839
+ ,297141
+ ,0
+ ,1
+ ,6
+ ,5
+ ,60719
+ ,364314
+ ,315372
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,22
+ ,2
+ ,1
+ ,60722
+ ,121444
+ ,312502
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,96
+ ,5
+ ,4
+ ,60379
+ ,301895
+ ,313729
+ ,0
+ ,1
+ ,2
+ ,1
+ ,60727
+ ,121454
+ ,315388
+ ,0
+ ,314
+ ,5
+ ,4
+ ,60720
+ ,0
+ ,315371
+ ,0
+ ,844
+ ,78
+ ,20
+ ,60925
+ ,426475
+ ,296139
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,26
+ ,1
+ ,1
+ ,60896
+ ,60896
+ ,313880
+ ,0
+ ,125
+ ,13
+ ,10
+ ,59734
+ ,776542
+ ,317698
+ ,0
+ ,304
+ ,15
+ ,12
+ ,62969
+ ,944535
+ ,295580
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,621
+ ,48
+ ,13
+ ,60720
+ ,0
+ ,308256
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,119
+ ,6
+ ,3
+ ,59118
+ ,354708
+ ,303677
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1595
+ ,17
+ ,10
+ ,60720
+ ,0
+ ,319369
+ ,0
+ ,312
+ ,14
+ ,3
+ ,58598
+ ,820372
+ ,318690
+ ,0
+ ,60
+ ,10
+ ,7
+ ,61124
+ ,611240
+ ,314049
+ ,0
+ ,587
+ ,12
+ ,10
+ ,59595
+ ,715140
+ ,325699
+ ,0
+ ,135
+ ,2
+ ,1
+ ,62065
+ ,124130
+ ,314210
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,514
+ ,52
+ ,15
+ ,78780
+ ,4096560
+ ,322378
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,1
+ ,4
+ ,4
+ ,60722
+ ,242888
+ ,315398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,1763
+ ,24
+ ,28
+ ,61600
+ ,184800
+ ,308336
+ ,0
+ ,180
+ ,11
+ ,9
+ ,59635
+ ,655985
+ ,316386
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,218
+ ,21
+ ,7
+ ,60720
+ ,0
+ ,315553
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,448
+ ,40
+ ,7
+ ,59781
+ ,2391240
+ ,323361
+ ,0
+ ,227
+ ,9
+ ,7
+ ,76644
+ ,689796
+ ,336639
+ ,0
+ ,174
+ ,1
+ ,3
+ ,64820
+ ,64820
+ ,307424
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,121
+ ,24
+ ,11
+ ,56178
+ ,1348272
+ ,295370
+ ,0
+ ,607
+ ,11
+ ,7
+ ,60436
+ ,664796
+ ,322340
+ ,0
+ ,2212
+ ,14
+ ,10
+ ,60720
+ ,0
+ ,319864
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,530
+ ,60
+ ,18
+ ,73433
+ ,4405980
+ ,317291
+ ,0
+ ,571
+ ,80
+ ,14
+ ,41477
+ ,3318160
+ ,280398
+ ,0
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,78
+ ,16
+ ,12
+ ,62700
+ ,1003200
+ ,317330
+ ,0
+ ,2489
+ ,40
+ ,29
+ ,67804
+ ,2712160
+ ,238125
+ ,0
+ ,131
+ ,6
+ ,3
+ ,59661
+ ,357966
+ ,327071
+ ,0
+ ,923
+ ,8
+ ,6
+ ,58620
+ ,468960
+ ,309038
+ ,0
+ ,72
+ ,3
+ ,3
+ ,60398
+ ,181194
+ ,314210
+ ,0
+ ,572
+ ,16
+ ,8
+ ,58580
+ ,937280
+ ,307930
+ ,0
+ ,397
+ ,10
+ ,10
+ ,62710
+ ,627100
+ ,322327
+ ,0
+ ,450
+ ,8
+ ,6
+ ,59325
+ ,474600
+ ,292136
+ ,0
+ ,622
+ ,7
+ ,8
+ ,60950
+ ,426650
+ ,263276
+ ,0
+ ,694
+ ,8
+ ,6
+ ,68060
+ ,544480
+ ,367655
+ ,1
+ ,3425
+ ,12
+ ,9
+ ,83620
+ ,1003440
+ ,283910
+ ,0
+ ,562
+ ,13
+ ,8
+ ,58456
+ ,759928
+ ,283587
+ ,0
+ ,4917
+ ,42
+ ,26
+ ,52811
+ ,2218062
+ ,243650
+ ,1
+ ,1442
+ ,118
+ ,239
+ ,121173
+ ,14298414
+ ,438493
+ ,0
+ ,529
+ ,9
+ ,7
+ ,63870
+ ,574830
+ ,296261
+ ,1
+ ,2126
+ ,138
+ ,41
+ ,21001
+ ,2898138
+ ,230621
+ ,0
+ ,1061
+ ,5
+ ,3
+ ,70415
+ ,352075
+ ,304252
+ ,0
+ ,776
+ ,9
+ ,8
+ ,64230
+ ,578070
+ ,333505
+ ,0
+ ,611
+ ,8
+ ,6
+ ,59190
+ ,473520
+ ,296919
+ ,1
+ ,1526
+ ,25
+ ,21
+ ,69351
+ ,1733775
+ ,278990
+ ,0
+ ,592
+ ,7
+ ,7
+ ,64270
+ ,449890
+ ,276898
+ ,0
+ ,1182
+ ,13
+ ,11
+ ,70694
+ ,919022
+ ,327007
+ ,0
+ ,621
+ ,16
+ ,11
+ ,68005
+ ,1088080
+ ,317046
+ ,0
+ ,989
+ ,11
+ ,12
+ ,58930
+ ,648230
+ ,304555
+ ,0
+ ,438
+ ,11
+ ,9
+ ,58320
+ ,641520
+ ,298096
+ ,0
+ ,726
+ ,3
+ ,3
+ ,69980
+ ,209940
+ ,231861
+ ,0
+ ,1303
+ ,61
+ ,57
+ ,69863
+ ,4261643
+ ,309422
+ ,1
+ ,7419
+ ,29
+ ,21
+ ,63255
+ ,1518120
+ ,286963
+ ,1
+ ,1164
+ ,17
+ ,15
+ ,57320
+ ,974440
+ ,269753
+ ,1
+ ,3310
+ ,33
+ ,32
+ ,75230
+ ,2482590
+ ,448243
+ ,0
+ ,1920
+ ,15
+ ,11
+ ,79420
+ ,555940
+ ,165404
+ ,0
+ ,965
+ ,3
+ ,2
+ ,73490
+ ,220470
+ ,204325
+ ,0
+ ,3256
+ ,66
+ ,23
+ ,35250
+ ,2326500
+ ,407159
+ ,1
+ ,1135
+ ,17
+ ,20
+ ,62285
+ ,1058845
+ ,290476
+ ,0
+ ,1270
+ ,26
+ ,24
+ ,69206
+ ,1799356
+ ,275311
+ ,0
+ ,661
+ ,3
+ ,1
+ ,65920
+ ,197760
+ ,246541
+ ,0
+ ,1013
+ ,2
+ ,1
+ ,69770
+ ,139540
+ ,253468
+ ,0
+ ,2844
+ ,67
+ ,74
+ ,72683
+ ,4869761
+ ,240897
+ ,1
+ ,11528
+ ,70
+ ,68
+ ,-14545
+ ,-1018150
+ ,-83265
+ ,0
+ ,6526
+ ,26
+ ,20
+ ,55830
+ ,1451580
+ ,-42143
+ ,0
+ ,2264
+ ,24
+ ,20
+ ,55174
+ ,1324176
+ ,272713
+ ,1
+ ,5109
+ ,97
+ ,82
+ ,67038
+ ,6301572
+ ,215362
+ ,0
+ ,3999
+ ,30
+ ,21
+ ,51252
+ ,1537560
+ ,42754
+ ,0
+ ,35624
+ ,223
+ ,244
+ ,157278
+ ,35072994
+ ,306275
+ ,0
+ ,9252
+ ,48
+ ,32
+ ,79510
+ ,3816480
+ ,253537
+ ,0
+ ,15236
+ ,90
+ ,86
+ ,77440
+ ,6969600
+ ,372631
+ ,0
+ ,18073
+ ,180
+ ,69
+ ,27284
+ ,4911120
+ ,-7170)
+ ,dim=c(7
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Orders'
+ ,'Dividends'
+ ,'TrDiv'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(7,431),dimnames=list(c('Group','Costs','Trades','Orders','Dividends','TrDiv','Wealth
'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '7'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth\r Group Costs Trades Orders Dividends TrDiv
1 6282154 1 162556 1081 807 213118 230380558
2 4321023 1 29790 309 444 81767 25266003
3 4111912 1 87550 458 412 153198 70164684
4 223193 0 84738 588 428 -26007 -15292116
5 1491348 1 54660 302 315 126942 37955658
6 1629616 1 42634 156 168 157214 24525384
7 1398893 0 40949 481 263 129352 62218312
8 1926517 1 45187 353 267 234817 75845891
9 983660 1 37704 452 228 60448 27322496
10 1443586 1 16275 109 129 47818 5212162
11 1073089 0 25830 115 104 245546 28237790
12 984885 0 12679 110 122 48020 5282200
13 1405225 1 18014 239 393 -1710 -408690
14 227132 0 43556 247 190 32648 8064056
15 929118 1 24811 505 280 95350 47388950
16 1071292 0 6575 159 63 151352 15589256
17 638830 0 7123 109 102 288170 31410530
18 856956 1 21950 519 265 114337 57397174
19 992426 1 37597 248 234 37884 9395232
20 444477 0 17821 373 277 122844 45820812
21 857217 1 12988 119 73 82340 9798460
22 711969 1 22330 84 67 79801 6703284
23 702380 0 13326 102 103 165548 16885896
24 358589 0 16189 295 290 116384 34333280
25 297978 0 7146 105 83 134028 14072940
26 585715 0 15824 64 56 63838 4085632
27 657954 1 27664 282 236 74996 20023932
28 209458 0 11920 182 73 31080 4009320
29 786690 0 8568 37 34 32168 1190216
30 439798 0 14416 361 139 49857 17998377
31 688779 1 3369 28 26 87161 2440508
32 574339 1 11819 85 70 106113 9019605
33 741409 1 6984 45 40 80570 3545080
34 597793 1 4519 49 42 102129 5004321
35 644190 0 2220 22 12 301670 6636740
36 377934 0 18562 155 211 102313 15858515
37 640273 0 10327 91 74 88577 8060507
38 697458 1 5336 81 80 112477 9110637
39 550608 1 2365 79 83 191778 15150462
40 207393 0 4069 145 131 79804 11571580
41 301607 0 8636 855 203 128294 104687904
42 345783 0 13718 61 56 96448 5883328
43 501749 0 4525 226 89 93811 21201286
44 379983 0 6869 105 88 117520 12339600
45 387475 0 4628 62 39 69159 4287858
46 377305 1 3689 25 25 101792 2443008
47 370837 1 4891 217 49 210568 5474768
48 430866 1 7489 322 149 136996 44112712
49 469107 0 4901 84 58 121920 10241280
50 194493 0 2284 33 41 76403 2521299
51 530670 1 3160 108 90 108094 11674152
52 518365 1 4150 150 136 134759 20213850
53 491303 1 7285 115 97 188873 21720395
54 527021 1 1134 162 63 146216 23686992
55 233773 1 4658 158 114 156608 24744064
56 405972 0 2384 97 77 61348 5950756
57 652925 0 3748 9 6 50350 453150
58 446211 0 5371 66 47 87720 5789520
59 341340 0 1285 107 51 99489 10645323
60 387699 1 9327 101 85 87419 8829319
61 493408 1 5565 47 43 94355 4434685
62 146494 0 1528 38 32 60326 2292388
63 414462 1 3122 34 25 94670 3218780
64 364304 1 7561 87 77 82425 6923700
65 355178 0 2675 79 54 59017 4662343
66 357760 0 13253 947 251 90829 86015063
67 261216 0 880 74 15 80791 5978534
68 397144 1 2053 53 44 100423 5322419
69 374943 0 1424 94 73 131116 12324904
70 424898 1 4036 63 85 100269 6316947
71 202055 1 3045 58 49 27330 1585140
72 378525 0 5119 49 38 39039 1912911
73 310768 0 1431 34 35 106885 3634090
74 325738 0 554 11 9 79285 872135
75 394510 0 1975 35 34 118881 4160835
76 247060 1 1765 20 20 77623 1319591
77 368078 0 1012 47 29 114768 5394096
78 236761 0 810 43 11 74015 3182645
79 312378 0 1280 117 52 69465 8127405
80 339836 1 666 171 13 117869 20155599
81 347385 0 1380 26 29 60982 1585532
82 426280 1 4677 75 66 90131 6579563
83 352850 0 876 59 33 138971 8199289
84 301881 0 814 18 15 39625 713250
85 377516 0 514 15 15 102725 1540875
86 357312 1 5692 72 68 64239 4625208
87 458343 0 3642 86 100 90262 7762532
88 354228 0 540 14 13 103960 1455440
89 308636 0 2099 64 45 106611 6823104
90 386212 0 567 11 14 103345 1136795
91 393343 0 2001 52 36 95551 4968652
92 378509 1 2949 41 40 82903 3399023
93 452469 0 2253 99 68 63593 6295707
94 364839 1 6533 75 29 126910 9518250
95 358649 0 1889 45 43 37527 1688715
96 376641 1 3055 43 30 60247 2590621
97 429112 0 272 8 9 112995 903960
98 330546 1 1414 198 22 70184 13896432
99 403560 0 2564 22 19 130140 2863080
100 317892 1 1383 11 9 73221 805431
101 307528 1 1261 33 31 76114 2511762
102 235133 0 975 23 19 90534 2082282
103 299243 0 3366 80 55 108479 8678320
104 314073 0 576 18 8 113761 2047698
105 368186 0 1686 40 28 68696 1923488
106 269661 0 746 23 29 71561 1645903
107 125390 1 3192 60 48 59831 3589860
108 510834 1 2045 20 16 97890 1957800
109 321896 0 5702 61 47 101481 5987379
110 249898 1 1932 36 20 72954 2626344
111 408881 0 936 30 22 67939 2038170
112 158492 1 3437 47 33 48022 2257034
113 292154 0 5131 71 44 86111 6113881
114 289513 1 2397 14 13 74020 1036280
115 378049 1 1389 9 6 57530 517770
116 343466 0 1503 39 35 56364 2198196
117 332743 0 402 26 8 84990 2209740
118 442882 0 2239 21 17 88590 1860390
119 214215 1 2234 16 11 77200 1235200
120 315688 0 837 69 21 61262 4227078
121 375195 0 10579 92 92 110309 10148428
122 334280 0 875 14 12 67000 938000
123 355864 0 1585 107 112 93099 9589197
124 480382 1 1659 29 25 107577 3119733
125 353058 1 2647 37 17 62920 2328040
126 217193 1 3294 23 23 75832 1744136
127 315380 0 0 0 0 60720 0
128 314533 0 94 7 10 60793 425551
129 318056 0 422 28 23 57935 1622180
130 315380 0 0 0 0 60720 0
131 314353 0 34 8 7 60630 485040
132 369448 0 1558 63 25 55637 3505131
133 315380 0 0 0 0 60720 0
134 312846 0 43 3 20 60887 182661
135 312075 0 645 5 4 60720 0
136 315009 0 316 9 4 60505 544545
137 318903 0 115 13 10 60945 792285
138 314887 0 5 2 1 60720 0
139 314913 0 897 5 4 60720 0
140 315380 0 0 0 0 60720 0
141 325506 0 389 14 8 58990 825860
142 315380 0 0 0 0 60720 0
143 298568 0 1002 15 11 56750 851250
144 315834 0 36 3 4 60894 182682
145 329784 0 460 15 15 63346 950190
146 312878 0 309 11 9 56535 621885
147 315380 0 0 0 0 60720 0
148 314987 0 9 6 7 60835 365010
149 325249 0 271 2 2 60720 0
150 315877 0 14 1 0 61016 61016
151 291650 0 520 10 7 58650 586500
152 305959 0 1766 73 46 60438 4411974
153 315380 0 0 0 5 60720 0
154 297765 0 458 11 7 58625 644875
155 315245 0 20 3 2 60938 182814
156 315380 0 0 0 0 60720 0
157 315380 0 0 0 0 60720 0
158 315236 0 98 2 2 61490 122980
159 336425 0 405 7 5 60845 425915
160 315380 0 0 0 0 60720 0
161 315380 0 0 0 0 60720 0
162 315380 0 0 0 0 60720 0
163 315380 0 0 0 0 60720 0
164 306268 0 483 27 7 60830 1642410
165 302187 0 454 51 24 63261 3226311
166 314882 0 47 3 1 60720 0
167 315380 0 0 0 0 60720 0
168 382712 0 757 19 18 45689 868091
169 341570 0 4655 393 55 60720 0
170 315380 0 0 0 0 60720 0
171 315380 0 0 0 0 60720 0
172 312412 0 36 4 3 61564 246256
173 315380 0 0 0 0 60720 0
174 309596 0 203 9 9 61938 557442
175 315380 0 0 0 0 60720 0
176 315547 0 126 10 8 60951 487608
177 313267 0 400 152 113 60720 0
178 316176 1 71 1 0 60745 60745
179 315380 0 0 0 0 60720 0
180 315380 0 0 0 0 60720 0
181 359335 0 972 34 19 71642 2435828
182 330068 1 531 10 11 71641 716410
183 314289 0 2461 57 25 55792 2120096
184 297413 1 378 52 16 71873 718730
185 314806 1 23 5 5 62555 312775
186 333210 1 638 14 11 60370 845180
187 352108 1 2300 29 23 64873 1037968
188 313332 0 149 5 6 62041 310205
189 291787 0 226 5 5 65745 328725
190 315380 0 0 0 0 60720 0
191 318745 0 275 4 7 59500 238000
192 315380 0 0 0 0 60720 0
193 315366 0 141 6 7 61630 369780
194 315380 0 0 0 0 60720 0
195 315688 0 28 2 3 60890 121780
196 315380 1 0 0 0 60720 0
197 409642 1 4980 91 89 113521 9081680
198 315380 1 0 0 0 60720 0
199 315380 1 0 0 0 60720 0
200 269587 1 472 20 19 80045 1600900
201 315380 1 0 0 0 60720 0
202 315380 1 0 0 0 60720 0
203 315380 1 0 0 0 60720 0
204 300962 1 203 27 12 50804 1371708
205 325479 1 496 17 12 87390 1485630
206 316155 1 10 2 5 61656 123312
207 318574 1 63 4 2 65688 262752
208 315380 1 0 0 0 60720 0
209 343613 1 1136 32 26 48522 1552704
210 306948 1 265 31 3 60720 0
211 315380 0 0 0 0 60720 0
212 315380 0 0 0 0 60720 0
213 330059 0 267 32 11 57640 1844480
214 288985 0 474 20 10 61977 1239540
215 304485 0 534 7 5 62620 438340
216 315380 0 0 0 2 60720 0
217 315688 0 15 8 6 60831 486648
218 317736 0 397 28 7 60646 1698088
219 315380 0 0 0 2 60720 0
220 322331 0 1866 29 28 56225 1124500
221 296656 0 288 4 3 60510 242040
222 315380 0 0 0 0 60720 0
223 315354 0 3 2 1 60698 121396
224 312161 0 468 21 20 60720 0
225 315576 0 20 2 1 60805 121610
226 314922 0 278 26 22 61404 1596504
227 314551 0 61 14 9 60720 0
228 315380 0 0 0 0 60720 0
229 312339 0 192 4 2 65276 261104
230 315380 0 0 0 0 60720 0
231 298700 0 317 9 7 63915 575235
232 321376 0 738 10 9 60720 0
233 315380 0 0 0 0 60720 0
234 303230 0 368 17 13 61686 1048662
235 315380 0 0 0 0 60720 0
236 315487 0 2 1 0 60743 60743
237 315380 0 0 0 0 60720 0
238 315793 0 53 6 6 60349 362094
239 315380 0 0 0 0 60720 0
240 315380 0 0 0 0 60720 0
241 315380 0 0 0 0 60720 0
242 312887 0 94 3 3 61360 184080
243 315380 0 0 0 0 60720 0
244 315637 0 24 8 7 59818 478544
245 324385 0 2332 4 2 72680 290720
246 315380 1 0 0 0 60720 0
247 315380 1 0 0 0 60720 0
248 308989 0 131 11 15 61808 679888
249 315380 1 0 0 0 60720 0
250 315380 1 0 0 0 60720 0
251 296702 0 206 9 9 53110 477990
252 315380 1 0 0 0 60720 0
253 307322 0 167 2 1 64245 128490
254 304376 0 622 73 38 73007 5329511
255 253588 0 2328 94 57 82732 7032220
256 315380 0 0 0 0 60720 0
257 309560 0 365 8 7 54820 438560
258 298466 0 364 35 26 47705 1669675
259 315380 1 0 0 0 60720 0
260 315380 0 0 0 0 60720 0
261 315380 1 0 0 0 60720 0
262 315380 1 0 0 0 60720 0
263 343929 0 226 12 13 72835 874020
264 331955 0 307 15 10 58856 882840
265 315380 1 0 0 0 60720 0
266 315380 0 0 0 0 60720 0
267 315380 1 0 0 0 60720 0
268 381180 0 188 11 9 77655 854205
269 315380 1 0 0 0 60720 0
270 331420 0 138 6 26 69817 418902
271 315380 1 0 0 0 60720 0
272 315380 1 0 0 0 60720 0
273 315380 1 0 0 0 60720 0
274 310201 0 125 12 19 60798 729576
275 315380 0 0 0 0 60720 0
276 320016 0 282 30 12 62452 1873560
277 320398 0 335 33 23 64175 2117775
278 315380 0 0 0 0 60720 0
279 291841 1 1324 117 29 67440 5530080
280 310670 0 176 28 8 68136 1907808
281 315380 1 0 0 0 60720 0
282 315380 0 0 0 0 60720 0
283 313491 0 249 72 26 56726 4084272
284 315380 0 0 0 0 60720 0
285 331323 0 333 13 9 70811 920543
286 315380 0 0 0 0 60720 0
287 319210 1 601 6 5 60720 0
288 318098 0 30 4 3 62045 248180
289 315380 1 0 0 0 60720 0
290 292754 0 249 62 13 54323 3368026
291 315380 1 0 0 0 60720 0
292 325176 0 165 24 12 62841 1508184
293 365959 0 453 21 19 81125 1703625
294 315380 0 0 0 0 60720 0
295 302409 0 53 14 10 59506 833084
296 340968 1 382 21 9 59365 1246665
297 315380 0 0 0 0 60720 0
298 315380 1 0 0 0 60720 0
299 315380 0 0 0 0 60720 0
300 315380 0 0 0 9 60720 0
301 313164 1 30 4 4 60798 243192
302 301164 0 290 2 1 58790 117580
303 315380 0 0 0 1 60720 0
304 315380 1 0 0 0 60720 0
305 344425 0 366 53 14 61808 3275824
306 315394 0 2 9 12 60735 546615
307 315380 1 0 0 0 60720 0
308 316647 1 209 13 19 64016 832208
309 309836 0 384 22 17 54683 1203026
310 315380 1 0 0 0 60720 0
311 315380 1 0 0 0 60720 0
312 346611 0 365 83 32 87192 7236936
313 315380 1 0 0 0 60720 0
314 322031 1 49 8 14 64107 512856
315 315656 0 3 4 8 60761 243044
316 339445 0 133 14 4 65990 923860
317 314964 0 32 1 0 59988 59988
318 297141 0 368 17 20 61167 1039839
319 315372 0 1 6 5 60719 364314
320 315380 0 0 0 0 60720 0
321 315380 0 0 0 0 60720 0
322 315380 0 0 0 0 60720 0
323 315380 0 0 0 0 60720 0
324 315380 1 0 0 0 60720 0
325 315380 0 0 0 0 60720 0
326 312502 0 22 2 1 60722 121444
327 315380 1 0 0 0 60720 0
328 315380 1 0 0 0 60720 0
329 315380 0 0 0 0 60720 0
330 315380 0 0 0 0 60720 0
331 315380 0 0 0 0 60720 0
332 315380 0 0 0 0 60720 0
333 315380 0 0 0 0 60720 0
334 313729 0 96 5 4 60379 301895
335 315388 0 1 2 1 60727 121454
336 315371 0 314 5 4 60720 0
337 296139 0 844 78 20 60925 426475
338 315380 0 0 0 0 60720 0
339 313880 0 26 1 1 60896 60896
340 317698 0 125 13 10 59734 776542
341 295580 0 304 15 12 62969 944535
342 315380 0 0 0 0 60720 0
343 315380 0 0 0 0 60720 0
344 315380 0 0 0 0 60720 0
345 308256 0 621 48 13 60720 0
346 315380 0 0 0 0 60720 0
347 303677 0 119 6 3 59118 354708
348 315380 0 0 0 0 60720 0
349 315380 0 0 0 0 60720 0
350 319369 0 1595 17 10 60720 0
351 318690 0 312 14 3 58598 820372
352 314049 0 60 10 7 61124 611240
353 325699 0 587 12 10 59595 715140
354 314210 0 135 2 1 62065 124130
355 315380 0 0 0 0 60720 0
356 315380 0 0 0 0 60720 0
357 322378 0 514 52 15 78780 4096560
358 315380 0 0 0 0 60720 0
359 315380 0 0 0 0 60720 0
360 315380 0 0 0 0 60720 0
361 315398 0 1 4 4 60722 242888
362 315380 0 0 0 0 60720 0
363 315380 0 0 0 0 60720 0
364 308336 1 1763 24 28 61600 184800
365 316386 0 180 11 9 59635 655985
366 315380 0 0 0 0 60720 0
367 315380 0 0 0 0 60720 0
368 315380 0 0 0 0 60720 0
369 315380 0 0 0 0 60720 0
370 315553 0 218 21 7 60720 0
371 315380 0 0 0 0 60720 0
372 323361 0 448 40 7 59781 2391240
373 336639 0 227 9 7 76644 689796
374 307424 0 174 1 3 64820 64820
375 315380 0 0 0 0 60720 0
376 315380 0 0 0 0 60720 0
377 295370 0 121 24 11 56178 1348272
378 322340 0 607 11 7 60436 664796
379 319864 0 2212 14 10 60720 0
380 315380 0 0 0 0 60720 0
381 315380 0 0 0 0 60720 0
382 317291 0 530 60 18 73433 4405980
383 280398 0 571 80 14 41477 3318160
384 315380 0 0 0 0 60720 0
385 317330 0 78 16 12 62700 1003200
386 238125 0 2489 40 29 67804 2712160
387 327071 0 131 6 3 59661 357966
388 309038 0 923 8 6 58620 468960
389 314210 0 72 3 3 60398 181194
390 307930 0 572 16 8 58580 937280
391 322327 0 397 10 10 62710 627100
392 292136 0 450 8 6 59325 474600
393 263276 0 622 7 8 60950 426650
394 367655 0 694 8 6 68060 544480
395 283910 1 3425 12 9 83620 1003440
396 283587 0 562 13 8 58456 759928
397 243650 0 4917 42 26 52811 2218062
398 438493 1 1442 118 239 121173 14298414
399 296261 0 529 9 7 63870 574830
400 230621 1 2126 138 41 21001 2898138
401 304252 0 1061 5 3 70415 352075
402 333505 0 776 9 8 64230 578070
403 296919 0 611 8 6 59190 473520
404 278990 1 1526 25 21 69351 1733775
405 276898 0 592 7 7 64270 449890
406 327007 0 1182 13 11 70694 919022
407 317046 0 621 16 11 68005 1088080
408 304555 0 989 11 12 58930 648230
409 298096 0 438 11 9 58320 641520
410 231861 0 726 3 3 69980 209940
411 309422 0 1303 61 57 69863 4261643
412 286963 1 7419 29 21 63255 1518120
413 269753 1 1164 17 15 57320 974440
414 448243 1 3310 33 32 75230 2482590
415 165404 0 1920 15 11 79420 555940
416 204325 0 965 3 2 73490 220470
417 407159 0 3256 66 23 35250 2326500
418 290476 1 1135 17 20 62285 1058845
419 275311 0 1270 26 24 69206 1799356
420 246541 0 661 3 1 65920 197760
421 253468 0 1013 2 1 69770 139540
422 240897 0 2844 67 74 72683 4869761
423 -83265 1 11528 70 68 -14545 -1018150
424 -42143 0 6526 26 20 55830 1451580
425 272713 0 2264 24 20 55174 1324176
426 215362 1 5109 97 82 67038 6301572
427 42754 0 3999 30 21 51252 1537560
428 306275 0 35624 223 244 157278 35072994
429 253537 0 9252 48 32 79510 3816480
430 372631 0 15236 90 86 77440 6969600
431 -7170 0 18073 180 69 27284 4911120
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Orders Dividends
3.217e+05 4.216e+04 1.408e+01 -2.725e+03 3.065e+03 -5.856e-01
TrDiv
1.778e-02
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-1188448 -15430 19435 29268 2617776
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.217e+05 2.993e+04 10.748 < 2e-16 ***
Group 4.216e+04 2.370e+04 1.779 0.0759 .
Costs 1.408e+01 1.940e+00 7.259 1.88e-12 ***
Trades -2.725e+03 2.277e+02 -11.968 < 2e-16 ***
Orders 3.065e+03 3.867e+02 7.926 2.02e-14 ***
Dividends -5.856e-01 4.011e-01 -1.460 0.1450
TrDiv 1.778e-02 1.557e-03 11.416 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 214100 on 424 degrees of freedom
Multiple R-squared: 0.7562, Adjusted R-squared: 0.7528
F-statistic: 219.2 on 6 and 424 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 3.908873e-42 1.954437e-42
[2,] 1.0000000 1.044946e-43 5.224728e-44
[3,] 1.0000000 2.145338e-45 1.072669e-45
[4,] 1.0000000 1.390369e-77 6.951846e-78
[5,] 1.0000000 9.215240e-81 4.607620e-81
[6,] 1.0000000 6.074299e-83 3.037149e-83
[7,] 1.0000000 2.638137e-94 1.319068e-94
[8,] 1.0000000 3.102184e-97 1.551092e-97
[9,] 1.0000000 1.055356e-99 5.276782e-100
[10,] 1.0000000 1.122485e-104 5.612427e-105
[11,] 1.0000000 1.189060e-110 5.945298e-111
[12,] 1.0000000 4.807750e-115 2.403875e-115
[13,] 1.0000000 8.390606e-117 4.195303e-117
[14,] 1.0000000 2.660184e-117 1.330092e-117
[15,] 1.0000000 1.186422e-124 5.932112e-125
[16,] 1.0000000 3.104100e-124 1.552050e-124
[17,] 1.0000000 1.138221e-125 5.691105e-126
[18,] 1.0000000 9.545324e-131 4.772662e-131
[19,] 1.0000000 6.637554e-131 3.318777e-131
[20,] 1.0000000 2.294628e-141 1.147314e-141
[21,] 1.0000000 3.684278e-145 1.842139e-145
[22,] 1.0000000 6.822887e-149 3.411443e-149
[23,] 1.0000000 8.966636e-151 4.483318e-151
[24,] 1.0000000 3.247504e-158 1.623752e-158
[25,] 1.0000000 3.243316e-160 1.621658e-160
[26,] 1.0000000 1.026129e-160 5.130646e-161
[27,] 1.0000000 3.089878e-164 1.544939e-164
[28,] 1.0000000 1.301769e-169 6.508847e-170
[29,] 1.0000000 4.695807e-175 2.347903e-175
[30,] 1.0000000 3.984533e-175 1.992267e-175
[31,] 1.0000000 3.919484e-176 1.959742e-176
[32,] 1.0000000 9.203015e-178 4.601508e-178
[33,] 1.0000000 3.881099e-177 1.940550e-177
[34,] 1.0000000 4.157070e-178 2.078535e-178
[35,] 1.0000000 3.561892e-177 1.780946e-177
[36,] 1.0000000 9.179370e-177 4.589685e-177
[37,] 1.0000000 2.859478e-176 1.429739e-176
[38,] 1.0000000 9.835408e-177 4.917704e-177
[39,] 1.0000000 2.389990e-177 1.194995e-177
[40,] 1.0000000 6.622025e-177 3.311012e-177
[41,] 1.0000000 2.786805e-177 1.393402e-177
[42,] 1.0000000 4.739549e-178 2.369775e-178
[43,] 1.0000000 7.221390e-179 3.610695e-179
[44,] 1.0000000 3.636777e-179 1.818389e-179
[45,] 1.0000000 1.009147e-178 5.045737e-179
[46,] 1.0000000 1.072704e-183 5.363522e-184
[47,] 1.0000000 2.443819e-183 1.221909e-183
[48,] 1.0000000 8.313775e-196 4.156888e-196
[49,] 1.0000000 2.982411e-196 1.491206e-196
[50,] 1.0000000 2.097453e-195 1.048727e-195
[51,] 1.0000000 1.061922e-195 5.309611e-196
[52,] 1.0000000 2.564566e-197 1.282283e-197
[53,] 1.0000000 6.161342e-199 3.080671e-199
[54,] 1.0000000 1.766825e-198 8.834127e-199
[55,] 1.0000000 2.406412e-198 1.203206e-198
[56,] 1.0000000 1.271599e-197 6.357996e-198
[57,] 1.0000000 1.195076e-198 5.975380e-199
[58,] 1.0000000 2.632198e-198 1.316099e-198
[59,] 1.0000000 1.915897e-197 9.579484e-198
[60,] 1.0000000 8.405656e-197 4.202828e-197
[61,] 1.0000000 1.539951e-196 7.699756e-197
[62,] 1.0000000 3.866107e-196 1.933054e-196
[63,] 1.0000000 4.167390e-197 2.083695e-197
[64,] 1.0000000 1.307778e-196 6.538890e-197
[65,] 1.0000000 1.117219e-195 5.586095e-196
[66,] 1.0000000 1.033454e-194 5.167271e-195
[67,] 1.0000000 1.564942e-194 7.824709e-195
[68,] 1.0000000 1.079778e-193 5.398891e-194
[69,] 1.0000000 1.303994e-193 6.519971e-194
[70,] 1.0000000 1.017860e-192 5.089298e-193
[71,] 1.0000000 8.121289e-193 4.060645e-193
[72,] 1.0000000 5.123164e-192 2.561582e-192
[73,] 1.0000000 5.731890e-192 2.865945e-192
[74,] 1.0000000 3.135882e-192 1.567941e-192
[75,] 1.0000000 2.381538e-191 1.190769e-191
[76,] 1.0000000 1.892175e-190 9.460876e-191
[77,] 1.0000000 2.074694e-190 1.037347e-190
[78,] 1.0000000 1.230158e-190 6.150791e-191
[79,] 1.0000000 8.278442e-190 4.139221e-190
[80,] 1.0000000 1.177325e-189 5.886625e-190
[81,] 1.0000000 9.189858e-189 4.594929e-189
[82,] 1.0000000 6.997138e-188 3.498569e-188
[83,] 1.0000000 3.132277e-187 1.566138e-187
[84,] 1.0000000 1.744358e-188 8.721792e-189
[85,] 1.0000000 1.202912e-187 6.014562e-188
[86,] 1.0000000 4.958419e-188 2.479210e-188
[87,] 1.0000000 4.372690e-188 2.186345e-188
[88,] 1.0000000 2.439241e-187 1.219621e-187
[89,] 1.0000000 6.159480e-187 3.079740e-187
[90,] 1.0000000 4.680896e-186 2.340448e-186
[91,] 1.0000000 4.144283e-185 2.072141e-185
[92,] 1.0000000 2.724700e-184 1.362350e-184
[93,] 1.0000000 2.513352e-185 1.256676e-185
[94,] 1.0000000 9.711445e-186 4.855722e-186
[95,] 1.0000000 1.863421e-186 9.317103e-187
[96,] 1.0000000 6.935786e-186 3.467893e-186
[97,] 1.0000000 2.343780e-185 1.171890e-185
[98,] 1.0000000 3.870440e-188 1.935220e-188
[99,] 1.0000000 8.919708e-190 4.459854e-190
[100,] 1.0000000 5.935222e-189 2.967611e-189
[101,] 1.0000000 7.442567e-189 3.721284e-189
[102,] 1.0000000 5.301103e-189 2.650551e-189
[103,] 1.0000000 1.180028e-189 5.900140e-190
[104,] 1.0000000 6.496618e-189 3.248309e-189
[105,] 1.0000000 4.546480e-188 2.273240e-188
[106,] 1.0000000 3.638798e-188 1.819399e-188
[107,] 1.0000000 1.374292e-187 6.871458e-188
[108,] 1.0000000 7.320454e-187 3.660227e-187
[109,] 1.0000000 2.781333e-187 1.390667e-187
[110,] 1.0000000 2.654477e-188 1.327239e-188
[111,] 1.0000000 2.067562e-187 1.033781e-187
[112,] 1.0000000 5.752587e-187 2.876293e-187
[113,] 1.0000000 4.714060e-186 2.357030e-186
[114,] 1.0000000 2.202281e-185 1.101140e-185
[115,] 1.0000000 1.740570e-185 8.702851e-186
[116,] 1.0000000 5.710215e-185 2.855107e-185
[117,] 1.0000000 1.815285e-185 9.076425e-186
[118,] 1.0000000 1.622450e-184 8.112252e-185
[119,] 1.0000000 1.509894e-183 7.549469e-184
[120,] 1.0000000 1.354318e-182 6.771590e-183
[121,] 1.0000000 1.200635e-181 6.003176e-182
[122,] 1.0000000 1.090444e-180 5.452220e-181
[123,] 1.0000000 1.599948e-180 7.999739e-181
[124,] 1.0000000 1.418869e-179 7.094347e-180
[125,] 1.0000000 1.317322e-178 6.586611e-179
[126,] 1.0000000 1.145655e-177 5.728273e-178
[127,] 1.0000000 1.002756e-176 5.013780e-177
[128,] 1.0000000 8.831197e-176 4.415599e-176
[129,] 1.0000000 7.678882e-175 3.839441e-175
[130,] 1.0000000 6.237804e-174 3.118902e-174
[131,] 1.0000000 5.392747e-173 2.696373e-173
[132,] 1.0000000 4.076088e-172 2.038044e-172
[133,] 1.0000000 3.499756e-171 1.749878e-171
[134,] 1.0000000 3.110761e-170 1.555380e-170
[135,] 1.0000000 2.673970e-169 1.336985e-169
[136,] 1.0000000 2.098609e-168 1.049304e-168
[137,] 1.0000000 1.724587e-167 8.622933e-168
[138,] 1.0000000 1.449684e-166 7.248418e-167
[139,] 1.0000000 1.238711e-165 6.193557e-166
[140,] 1.0000000 9.131562e-165 4.565781e-165
[141,] 1.0000000 7.573076e-164 3.786538e-164
[142,] 1.0000000 6.302734e-163 3.151367e-163
[143,] 1.0000000 5.412046e-162 2.706023e-162
[144,] 1.0000000 4.483365e-161 2.241683e-161
[145,] 1.0000000 3.755662e-160 1.877831e-160
[146,] 1.0000000 3.067906e-159 1.533953e-159
[147,] 1.0000000 2.480127e-158 1.240064e-158
[148,] 1.0000000 1.996510e-157 9.982548e-158
[149,] 1.0000000 1.612870e-156 8.064349e-157
[150,] 1.0000000 9.553520e-156 4.776760e-156
[151,] 1.0000000 7.597415e-155 3.798707e-155
[152,] 1.0000000 6.015653e-154 3.007827e-154
[153,] 1.0000000 4.742465e-153 2.371232e-153
[154,] 1.0000000 3.722376e-152 1.861188e-152
[155,] 1.0000000 2.857934e-151 1.428967e-151
[156,] 1.0000000 1.893190e-150 9.465952e-151
[157,] 1.0000000 1.462255e-149 7.311275e-150
[158,] 1.0000000 1.129286e-148 5.646429e-149
[159,] 1.0000000 1.521012e-149 7.605062e-150
[160,] 1.0000000 4.066609e-154 2.033305e-154
[161,] 1.0000000 3.332437e-153 1.666218e-153
[162,] 1.0000000 2.718575e-152 1.359288e-152
[163,] 1.0000000 2.238306e-151 1.119153e-151
[164,] 1.0000000 1.809575e-150 9.047876e-151
[165,] 1.0000000 1.487498e-149 7.437492e-150
[166,] 1.0000000 1.191570e-148 5.957848e-149
[167,] 1.0000000 9.534674e-148 4.767337e-148
[168,] 1.0000000 5.968665e-148 2.984332e-148
[169,] 1.0000000 4.844685e-147 2.422342e-147
[170,] 1.0000000 3.860345e-146 1.930173e-146
[171,] 1.0000000 3.062135e-145 1.531067e-145
[172,] 1.0000000 1.552456e-144 7.762282e-145
[173,] 1.0000000 1.224564e-143 6.122821e-144
[174,] 1.0000000 4.691638e-143 2.345819e-143
[175,] 1.0000000 3.631547e-142 1.815774e-142
[176,] 1.0000000 2.824773e-141 1.412386e-141
[177,] 1.0000000 1.859185e-140 9.295926e-141
[178,] 1.0000000 3.590788e-140 1.795394e-140
[179,] 1.0000000 2.809971e-139 1.404986e-139
[180,] 1.0000000 1.845239e-138 9.226197e-139
[181,] 1.0000000 1.414216e-137 7.071081e-138
[182,] 1.0000000 1.030555e-136 5.152775e-137
[183,] 1.0000000 7.826531e-136 3.913266e-136
[184,] 1.0000000 5.949083e-135 2.974542e-135
[185,] 1.0000000 4.476732e-134 2.238366e-134
[186,] 1.0000000 3.358439e-133 1.679220e-133
[187,] 1.0000000 2.533453e-132 1.266726e-132
[188,] 1.0000000 9.281480e-132 4.640740e-132
[189,] 1.0000000 6.968238e-131 3.484119e-131
[190,] 1.0000000 5.207917e-130 2.603959e-130
[191,] 1.0000000 7.281771e-130 3.640886e-130
[192,] 1.0000000 5.432854e-129 2.716427e-129
[193,] 1.0000000 4.034670e-128 2.017335e-128
[194,] 1.0000000 2.982435e-127 1.491218e-127
[195,] 1.0000000 2.184677e-126 1.092338e-126
[196,] 1.0000000 1.231170e-125 6.155851e-126
[197,] 1.0000000 8.941063e-125 4.470531e-125
[198,] 1.0000000 6.432923e-124 3.216462e-124
[199,] 1.0000000 4.647641e-123 2.323821e-123
[200,] 1.0000000 1.060691e-122 5.303455e-123
[201,] 1.0000000 7.377447e-122 3.688723e-122
[202,] 1.0000000 5.217086e-121 2.608543e-121
[203,] 1.0000000 3.672792e-120 1.836396e-120
[204,] 1.0000000 2.162860e-119 1.081430e-119
[205,] 1.0000000 1.346674e-118 6.733372e-119
[206,] 1.0000000 9.478394e-118 4.739197e-118
[207,] 1.0000000 6.572324e-117 3.286162e-117
[208,] 1.0000000 4.544835e-116 2.272418e-116
[209,] 1.0000000 3.058214e-115 1.529107e-115
[210,] 1.0000000 2.091731e-114 1.045866e-114
[211,] 1.0000000 7.033129e-114 3.516565e-114
[212,] 1.0000000 4.702885e-113 2.351442e-113
[213,] 1.0000000 3.178368e-112 1.589184e-112
[214,] 1.0000000 2.139505e-111 1.069752e-111
[215,] 1.0000000 1.220844e-110 6.104218e-111
[216,] 1.0000000 8.140874e-110 4.070437e-110
[217,] 1.0000000 5.419023e-109 2.709511e-109
[218,] 1.0000000 3.366314e-108 1.683157e-108
[219,] 1.0000000 2.214464e-107 1.107232e-107
[220,] 1.0000000 1.461147e-106 7.305736e-107
[221,] 1.0000000 9.524529e-106 4.762264e-106
[222,] 1.0000000 5.988260e-105 2.994130e-105
[223,] 1.0000000 3.130107e-104 1.565053e-104
[224,] 1.0000000 2.014392e-103 1.007196e-103
[225,] 1.0000000 1.303920e-102 6.519602e-103
[226,] 1.0000000 8.313912e-102 4.156956e-102
[227,] 1.0000000 5.275470e-101 2.637735e-101
[228,] 1.0000000 3.332327e-100 1.666163e-100
[229,] 1.0000000 2.087629e-99 1.043814e-99
[230,] 1.0000000 1.306239e-98 6.531196e-99
[231,] 1.0000000 8.134663e-98 4.067331e-98
[232,] 1.0000000 5.041957e-97 2.520979e-97
[233,] 1.0000000 3.136886e-96 1.568443e-96
[234,] 1.0000000 1.925868e-95 9.629338e-96
[235,] 1.0000000 1.171682e-94 5.858409e-95
[236,] 1.0000000 6.229775e-94 3.114888e-94
[237,] 1.0000000 3.818448e-93 1.909224e-93
[238,] 1.0000000 2.329318e-92 1.164659e-92
[239,] 1.0000000 1.417569e-91 7.087847e-92
[240,] 1.0000000 8.564235e-91 4.282117e-91
[241,] 1.0000000 5.149227e-90 2.574613e-90
[242,] 1.0000000 3.086889e-89 1.543444e-89
[243,] 1.0000000 1.838023e-88 9.190116e-89
[244,] 1.0000000 1.081284e-87 5.406419e-88
[245,] 1.0000000 4.510494e-87 2.255247e-87
[246,] 1.0000000 2.888843e-87 1.444422e-87
[247,] 1.0000000 1.691541e-86 8.457707e-87
[248,] 1.0000000 9.506920e-86 4.753460e-86
[249,] 1.0000000 5.245745e-85 2.622872e-85
[250,] 1.0000000 3.070208e-84 1.535104e-84
[251,] 1.0000000 1.762620e-83 8.813099e-84
[252,] 1.0000000 1.021562e-82 5.107811e-83
[253,] 1.0000000 5.891581e-82 2.945790e-82
[254,] 1.0000000 3.125666e-81 1.562833e-81
[255,] 1.0000000 1.463538e-80 7.317690e-81
[256,] 1.0000000 8.329938e-80 4.164969e-80
[257,] 1.0000000 4.646510e-79 2.323255e-79
[258,] 1.0000000 2.618481e-78 1.309240e-78
[259,] 1.0000000 8.576041e-78 4.288020e-78
[260,] 1.0000000 4.801570e-77 2.400785e-77
[261,] 1.0000000 2.438888e-76 1.219444e-76
[262,] 1.0000000 1.352480e-75 6.762402e-76
[263,] 1.0000000 7.462130e-75 3.731065e-75
[264,] 1.0000000 4.096175e-74 2.048087e-74
[265,] 1.0000000 2.214539e-73 1.107270e-73
[266,] 1.0000000 1.184138e-72 5.920690e-73
[267,] 1.0000000 6.333426e-72 3.166713e-72
[268,] 1.0000000 3.390794e-71 1.695397e-71
[269,] 1.0000000 1.785530e-70 8.927650e-71
[270,] 1.0000000 6.201449e-70 3.100724e-70
[271,] 1.0000000 2.995341e-69 1.497670e-69
[272,] 1.0000000 1.581789e-68 7.908944e-69
[273,] 1.0000000 8.174506e-68 4.087253e-68
[274,] 1.0000000 4.158217e-67 2.079109e-67
[275,] 1.0000000 2.127296e-66 1.063648e-66
[276,] 1.0000000 1.080600e-65 5.403002e-66
[277,] 1.0000000 5.470083e-65 2.735042e-65
[278,] 1.0000000 2.698814e-64 1.349407e-64
[279,] 1.0000000 1.350834e-63 6.754172e-64
[280,] 1.0000000 6.848623e-63 3.424312e-63
[281,] 1.0000000 2.968749e-62 1.484375e-62
[282,] 1.0000000 1.490526e-61 7.452632e-62
[283,] 1.0000000 7.306993e-61 3.653496e-61
[284,] 1.0000000 3.128628e-60 1.564314e-60
[285,] 1.0000000 1.521587e-59 7.607937e-60
[286,] 1.0000000 7.474342e-59 3.737171e-59
[287,] 1.0000000 3.293209e-58 1.646604e-58
[288,] 1.0000000 1.577055e-57 7.885275e-58
[289,] 1.0000000 7.642456e-57 3.821228e-57
[290,] 1.0000000 3.619388e-56 1.809694e-56
[291,] 1.0000000 1.672293e-55 8.361466e-56
[292,] 1.0000000 7.941475e-55 3.970737e-55
[293,] 1.0000000 3.768594e-54 1.884297e-54
[294,] 1.0000000 1.742584e-53 8.712922e-54
[295,] 1.0000000 8.166366e-53 4.083183e-53
[296,] 1.0000000 3.502691e-52 1.751345e-52
[297,] 1.0000000 1.592927e-51 7.964636e-52
[298,] 1.0000000 7.342702e-51 3.671351e-51
[299,] 1.0000000 3.323503e-50 1.661752e-50
[300,] 1.0000000 1.448429e-49 7.242144e-50
[301,] 1.0000000 6.563139e-49 3.281570e-49
[302,] 1.0000000 2.956061e-48 1.478030e-48
[303,] 1.0000000 6.519018e-48 3.259509e-48
[304,] 1.0000000 2.919433e-47 1.459717e-47
[305,] 1.0000000 1.292070e-46 6.460349e-47
[306,] 1.0000000 5.576111e-46 2.788056e-46
[307,] 1.0000000 2.288733e-45 1.144366e-45
[308,] 1.0000000 9.816873e-45 4.908437e-45
[309,] 1.0000000 4.252258e-44 2.126129e-44
[310,] 1.0000000 1.813617e-43 9.068084e-44
[311,] 1.0000000 7.648296e-43 3.824148e-43
[312,] 1.0000000 3.204841e-42 1.602421e-42
[313,] 1.0000000 1.334292e-41 6.671458e-42
[314,] 1.0000000 5.519206e-41 2.759603e-41
[315,] 1.0000000 2.313664e-40 1.156832e-40
[316,] 1.0000000 9.448908e-40 4.724454e-40
[317,] 1.0000000 3.874082e-39 1.937041e-39
[318,] 1.0000000 1.592963e-38 7.964815e-39
[319,] 1.0000000 6.505751e-38 3.252875e-38
[320,] 1.0000000 2.588829e-37 1.294415e-37
[321,] 1.0000000 1.023169e-36 5.115845e-37
[322,] 1.0000000 4.016102e-36 2.008051e-36
[323,] 1.0000000 1.565488e-35 7.827442e-36
[324,] 1.0000000 6.059778e-35 3.029889e-35
[325,] 1.0000000 2.339080e-34 1.169540e-34
[326,] 1.0000000 8.952955e-34 4.476477e-34
[327,] 1.0000000 3.284060e-33 1.642030e-33
[328,] 1.0000000 1.128166e-32 5.640828e-33
[329,] 1.0000000 4.225197e-32 2.112598e-32
[330,] 1.0000000 1.579129e-31 7.895643e-32
[331,] 1.0000000 5.779044e-31 2.889522e-31
[332,] 1.0000000 2.138975e-30 1.069488e-30
[333,] 1.0000000 7.783977e-30 3.891988e-30
[334,] 1.0000000 2.811085e-29 1.405542e-29
[335,] 1.0000000 1.007367e-28 5.036837e-29
[336,] 1.0000000 3.170151e-28 1.585076e-28
[337,] 1.0000000 1.121482e-27 5.607408e-28
[338,] 1.0000000 4.026490e-27 2.013245e-27
[339,] 1.0000000 1.402266e-26 7.011330e-27
[340,] 1.0000000 4.843943e-26 2.421971e-26
[341,] 1.0000000 1.277061e-25 6.385307e-26
[342,] 1.0000000 4.305731e-25 2.152866e-25
[343,] 1.0000000 1.474220e-24 7.371099e-25
[344,] 1.0000000 4.487394e-24 2.243697e-24
[345,] 1.0000000 1.503227e-23 7.516136e-24
[346,] 1.0000000 4.951910e-23 2.475955e-23
[347,] 1.0000000 1.616816e-22 8.084078e-23
[348,] 1.0000000 4.367584e-22 2.183792e-22
[349,] 1.0000000 1.407272e-21 7.036359e-22
[350,] 1.0000000 4.492721e-21 2.246361e-21
[351,] 1.0000000 1.420939e-20 7.104693e-21
[352,] 1.0000000 4.461536e-20 2.230768e-20
[353,] 1.0000000 1.384384e-19 6.921920e-20
[354,] 1.0000000 4.253691e-19 2.126845e-19
[355,] 1.0000000 1.183508e-18 5.917541e-19
[356,] 1.0000000 3.546460e-18 1.773230e-18
[357,] 1.0000000 1.060546e-17 5.302731e-18
[358,] 1.0000000 3.138536e-17 1.569268e-17
[359,] 1.0000000 9.189677e-17 4.594838e-17
[360,] 1.0000000 2.661682e-16 1.330841e-16
[361,] 1.0000000 7.332853e-16 3.666426e-16
[362,] 1.0000000 2.080129e-15 1.040065e-15
[363,] 1.0000000 6.049051e-15 3.024526e-15
[364,] 1.0000000 1.724903e-14 8.624513e-15
[365,] 1.0000000 4.908203e-14 2.454102e-14
[366,] 1.0000000 1.333673e-13 6.668367e-14
[367,] 1.0000000 3.579477e-13 1.789738e-13
[368,] 1.0000000 1.000170e-12 5.000849e-13
[369,] 1.0000000 2.498746e-12 1.249373e-12
[370,] 1.0000000 4.887412e-12 2.443706e-12
[371,] 1.0000000 1.252272e-11 6.261358e-12
[372,] 1.0000000 3.163682e-11 1.581841e-11
[373,] 1.0000000 7.443512e-11 3.721756e-11
[374,] 1.0000000 1.913755e-10 9.568776e-11
[375,] 1.0000000 4.700726e-10 2.350363e-10
[376,] 1.0000000 1.175939e-09 5.879697e-10
[377,] 1.0000000 2.553986e-09 1.276993e-09
[378,] 1.0000000 5.768896e-09 2.884448e-09
[379,] 1.0000000 1.310138e-08 6.550692e-09
[380,] 1.0000000 3.016517e-08 1.508258e-08
[381,] 1.0000000 7.001368e-08 3.500684e-08
[382,] 0.9999999 1.521280e-07 7.606399e-08
[383,] 0.9999998 3.528164e-07 1.764082e-07
[384,] 0.9999996 8.262271e-07 4.131135e-07
[385,] 0.9999993 1.347806e-06 6.739028e-07
[386,] 0.9999985 2.914259e-06 1.457130e-06
[387,] 0.9999968 6.495370e-06 3.247685e-06
[388,] 0.9999932 1.356252e-05 6.781259e-06
[389,] 0.9999881 2.378738e-05 1.189369e-05
[390,] 0.9999747 5.050905e-05 2.525452e-05
[391,] 0.9999476 1.047618e-04 5.238090e-05
[392,] 0.9998931 2.138492e-04 1.069246e-04
[393,] 0.9998151 3.698316e-04 1.849158e-04
[394,] 0.9996518 6.963006e-04 3.481503e-04
[395,] 0.9993570 1.285937e-03 6.429687e-04
[396,] 0.9987532 2.493687e-03 1.246843e-03
[397,] 0.9978238 4.352468e-03 2.176234e-03
[398,] 0.9961317 7.736677e-03 3.868338e-03
[399,] 0.9939039 1.219229e-02 6.096145e-03
[400,] 0.9902459 1.950828e-02 9.754138e-03
[401,] 0.9827845 3.443101e-02 1.721550e-02
[402,] 0.9709392 5.812166e-02 2.906083e-02
[403,] 0.9537908 9.241842e-02 4.620921e-02
[404,] 0.9250910 1.498180e-01 7.490902e-02
[405,] 0.9238618 1.522764e-01 7.613821e-02
[406,] 0.8936620 2.126759e-01 1.063380e-01
[407,] 0.8401800 3.196400e-01 1.598200e-01
[408,] 0.9696029 6.079416e-02 3.039708e-02
[409,] 0.9434311 1.131379e-01 5.656895e-02
[410,] 0.8923107 2.153785e-01 1.076893e-01
[411,] 0.8148940 3.702120e-01 1.851060e-01
[412,] 0.6874003 6.251993e-01 3.125997e-01
> postscript(file="/var/www/html/rcomp/tmp/1lsvn1291394141.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/2lsvn1291394141.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/3w1u81291394141.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/4w1u81291394141.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/5w1u81291394141.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
131691.2321 2617776.3344 1343261.9381 -744280.2782 -384882.4697
6 7 8 9 10
231829.6052 -24799.1285 -140602.1496 171701.7832 687631.7910
11 12 13 14 15
24217.3337 344779.5263 240841.5619 -741217.9665 -52513.3362
16 17 18 19 20
708793.0810 -188311.0806 -167054.1605 -86916.7805 -703145.6635
21 22 23 24 25
285123.3410 -15121.8885 -47851.8834 -818053.2764 -264218.2790
26 27 28 29 30
8784.4069 -362208.6164 -60845.5727 338692.8381 182219.7224
31 32 33 34 35
281788.1807 -36995.6012 263445.7065 145992.6598 373117.5424
36 37 38 39 40
-651360.7695 102982.2074 137960.6316 -42632.1754 -336872.4361
41 42 43 44 45
-219519.1761 -222531.9421 137566.0210 -172492.4321 14357.7386
46 47 48 49 50
-30779.2236 405337.5305 -321471.9194 18934.8649 -195140.8815
51 52 53 54 55
-3391.6077 -192350.0646 -234497.9226 60187.4370 -462598.5773
56 57 58 59 60
9233.5964 306052.0936 33188.9153 5900.7235 -198482.0380
61 62 63 64 65
23938.1649 -196629.2342 20928.2372 -179696.4179 -2681.0677
66 67 68 69 70
185234.2138 23891.0866 -21817.5531 -76648.7399 -138163.7547
71 72 73 74 75
-208938.8859 -9289.8180 -47672.1817 29587.7757 31861.7253
76 77 78 79 80
-126419.0457 42686.2541 -26071.1896 28375.2746 103544.1770
81 82 83 84 85
-4214.0980 -65473.3759 14124.1325 -17642.2383 76278.7647
86 87 88 89 90
-143459.0879 -71856.1787 58271.2104 -64948.5081 83938.2662
91 92 93 94 95
42507.8170 -49580.7887 85793.1811 -70344.4744 -6814.0289
96 97 98 99 100
-15737.3701 147930.4943 213063.0816 72826.2490 -34459.7025
101 102 103 104 105
-79216.1575 -79815.2978 -111110.0897 39046.6756 52007.6985
106 107 108 109 110
-76065.7715 -295762.5098 146195.4603 -104872.2276 -108290.3137
111 112 113 114 115
91917.1948 -238787.2513 -101375.4457 -84835.3560 25282.2672
116 117 118 119 120
-6420.5109 62241.0108 113622.0525 -147931.6276 66650.5828
121 122 123 124 125
-242477.7777 24226.5214 -155730.5335 103139.2098 -3850.1514
126 127 128 129 130
-187428.7167 29268.1262 8003.2527 1350.2192 29268.1262
131 132 133 134 135
19435.3611 91186.7223 29268.1262 -30144.0315 18248.3820
136 137 138 139 140
26910.1165 21999.4275 31090.4241 17538.2510 29268.1262
141 142 143 144 145
31858.1740 29268.1262 -11943.4238 21985.8429 16747.6812
146 147 148 149 150
11303.9014 29268.1262 17224.4191 34642.2003 31382.0373
151 152 153 154 155
-7622.8240 -25650.9295 13943.2545 1038.1341 27775.4888
156 157 158 159 160
29268.1262 29268.1262 25329.7892 40865.2598 29268.1262
161 162 163 164 165
29268.1262 29268.1262 29268.1262 36353.2418 19251.7140
166 167 168 169 170
33219.4047 29268.1262 58320.1911 892399.6433 29268.1262
171 172 173 174 175
29268.1262 23616.3821 29268.1262 8373.1483 29268.1262
176 177 178 179 180
21862.0309 89432.1141 -11439.4325 29268.1262 29268.1262
181 182 183 184 185
57060.1061 -18485.6793 31673.1233 50247.3889 -19977.3842
186 187 188 189 190
-14838.0182 -16030.6275 15618.3139 -2106.0886 29268.1262
191 192 193 194 195
13262.5146 29268.1262 16125.6181 29268.1262 23372.3899
196 197 198 199 200
-12895.9151 -144048.2459 -12895.9151 -12895.9151 -86204.1344
201 202 203 204 205
-12895.9151 -12895.9151 -12895.9151 -23558.2587 -11021.0744
206 207 208 209 210
-23779.8406 -7579.1442 -12895.9151 -27880.4863 50231.4022
211 212 213 214 215
29268.1262 29268.1262 59092.2205 -1242.1168 7927.4901
216 217 218 219 220
23138.1775 24191.9708 50659.9525 23138.1775 -19460.1087
221 222 223 224 225
3769.9938 29268.1262 29414.7303 15392.3051 29456.2253
226 227 228 229 230
345.8293 38150.1854 29268.1262 26321.6280 29268.1262
231 232 233 234 235
2843.3897 24541.7601 29268.1262 347.1197 29268.1262
236 237 238 239 240
31005.9853 29268.1262 20243.1219 29268.1262 29268.1262
241 242 243 244 245
29268.1262 21535.2133 29268.1262 20500.1432 12045.8516
246 247 248 249 250
-12895.9151 -12895.9151 -6412.0358 -12895.9151 -12895.9151
251 252 253 254 255
-8320.2368 -12895.9151 21024.5792 4443.0129 -95940.8109
256 257 258 259 260
29268.1262 7405.9368 -14375.4965 -12895.9151 29268.1262
261 262 263 264 265
-12895.9151 -12895.9151 39051.9120 34965.7507 -12895.9151
266 267 268 269 270
29268.1262 -12895.9151 89547.2459 -12895.9151 -22091.7116
271 272 273 274 275
-12895.9151 -12895.9151 -12895.9151 -16124.8231 29268.1262
276 277 278 279 280
42623.2566 13388.0309 29268.1262 80534.1352 44298.5890
281 282 283 284 285
-12895.9151 29268.1262 65466.0576 29268.1262 37912.3754
286 287 288 289 290
29268.1262 -16500.7829 29634.3242 -12895.9151 68645.4134
291 292 293 294 295
-12895.9151 39801.4098 54131.1182 29268.1262 7535.8105
296 297 298 299 300
14006.3083 29268.1262 -12895.9151 29268.1262 1683.3572
301 302 303 304 305
-21170.2426 10134.3494 26203.1519 -12895.9151 97098.1523
306 307 308 309 310
7294.2924 -12895.9151 -50240.2563 1249.7813 -12895.9151
311 312 313 314 315
-12895.9151 70339.2670 -12895.9151 -35175.0997 11587.0227
316 317 318 319 320
64018.5109 29631.8963 -27343.7769 23796.4518 29268.1262
321 322 323 324 325
29268.1262 29268.1262 29268.1262 -12895.9151 29268.1262
326 327 328 329 330
26308.4131 -12895.9151 -12895.9151 29268.1262 29268.1262
331 332 333 334 335
29268.1262 29268.1262 29268.1262 22065.9821 29492.8409
336 337 338 339 340
26204.8240 141959.3090 29268.1262 26082.9678 20224.3425
341 342 343 344 345
-6185.1728 29268.1262 29268.1262 29268.1262 104371.9698
346 347 348 349 350
29268.1262 15803.2182 29268.1262 29268.1262 26480.6749
351 352 353 354 355
41319.2054 22261.8646 20005.2039 27164.0740 29268.1262
356 357 358 359 360
29268.1262 62525.8537 29268.1262 29268.1262 29268.1262
361 362 363 364 365
23597.0153 29268.1262 29268.1262 -67943.7223 17837.3355
366 367 368 369 370
29268.1262 29268.1262 29268.1262 29268.1262 62148.9428
371 372 373 374 375
29268.1262 75442.7081 47466.9676 13641.2611 29268.1262
376 377 378 379 380
29268.1262 12620.8166 24221.5966 10112.3784 29268.1262
381 382 383 384 385
29268.1262 61189.8980 91110.8221 29268.1262 20271.8496
386 387 388 389 390
-106966.4255 39288.3148 3777.1431 22655.9444 14935.5040
391 392 393 394 395
17246.8195 -6152.4952 -44485.5778 69803.8539 -91897.9029
396 397 398 399 400
-14362.6556 -120978.1743 -539798.0202 -2599.6976 48066.8103
401 402 403 404 405
7051.7420 28254.8104 -3696.2117 -92768.8288 -25845.1997
406 407 408 409 410
15471.1292 17005.1009 -14854.1938 -4598.1795 -63801.2624
411 412 413 414 415
-73896.2632 -156603.4022 -93868.7126 29585.1004 -139508.1245
416 417 418 419 420
-89769.1764 128309.4869 -86655.2669 -58399.6203 -44237.0641
421 422 423 424 425
-41702.0976 -209029.6603 -617473.3856 -439248.1597 -67953.7316
426 427 428 429 430
-280137.2562 -315144.1159 -1188447.8307 -186944.8306 -260412.3078
431
-375551.4045
> postscript(file="/var/www/html/rcomp/tmp/66tbb1291394141.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 131691.2321 NA
1 2617776.3344 131691.2321
2 1343261.9381 2617776.3344
3 -744280.2782 1343261.9381
4 -384882.4697 -744280.2782
5 231829.6052 -384882.4697
6 -24799.1285 231829.6052
7 -140602.1496 -24799.1285
8 171701.7832 -140602.1496
9 687631.7910 171701.7832
10 24217.3337 687631.7910
11 344779.5263 24217.3337
12 240841.5619 344779.5263
13 -741217.9665 240841.5619
14 -52513.3362 -741217.9665
15 708793.0810 -52513.3362
16 -188311.0806 708793.0810
17 -167054.1605 -188311.0806
18 -86916.7805 -167054.1605
19 -703145.6635 -86916.7805
20 285123.3410 -703145.6635
21 -15121.8885 285123.3410
22 -47851.8834 -15121.8885
23 -818053.2764 -47851.8834
24 -264218.2790 -818053.2764
25 8784.4069 -264218.2790
26 -362208.6164 8784.4069
27 -60845.5727 -362208.6164
28 338692.8381 -60845.5727
29 182219.7224 338692.8381
30 281788.1807 182219.7224
31 -36995.6012 281788.1807
32 263445.7065 -36995.6012
33 145992.6598 263445.7065
34 373117.5424 145992.6598
35 -651360.7695 373117.5424
36 102982.2074 -651360.7695
37 137960.6316 102982.2074
38 -42632.1754 137960.6316
39 -336872.4361 -42632.1754
40 -219519.1761 -336872.4361
41 -222531.9421 -219519.1761
42 137566.0210 -222531.9421
43 -172492.4321 137566.0210
44 14357.7386 -172492.4321
45 -30779.2236 14357.7386
46 405337.5305 -30779.2236
47 -321471.9194 405337.5305
48 18934.8649 -321471.9194
49 -195140.8815 18934.8649
50 -3391.6077 -195140.8815
51 -192350.0646 -3391.6077
52 -234497.9226 -192350.0646
53 60187.4370 -234497.9226
54 -462598.5773 60187.4370
55 9233.5964 -462598.5773
56 306052.0936 9233.5964
57 33188.9153 306052.0936
58 5900.7235 33188.9153
59 -198482.0380 5900.7235
60 23938.1649 -198482.0380
61 -196629.2342 23938.1649
62 20928.2372 -196629.2342
63 -179696.4179 20928.2372
64 -2681.0677 -179696.4179
65 185234.2138 -2681.0677
66 23891.0866 185234.2138
67 -21817.5531 23891.0866
68 -76648.7399 -21817.5531
69 -138163.7547 -76648.7399
70 -208938.8859 -138163.7547
71 -9289.8180 -208938.8859
72 -47672.1817 -9289.8180
73 29587.7757 -47672.1817
74 31861.7253 29587.7757
75 -126419.0457 31861.7253
76 42686.2541 -126419.0457
77 -26071.1896 42686.2541
78 28375.2746 -26071.1896
79 103544.1770 28375.2746
80 -4214.0980 103544.1770
81 -65473.3759 -4214.0980
82 14124.1325 -65473.3759
83 -17642.2383 14124.1325
84 76278.7647 -17642.2383
85 -143459.0879 76278.7647
86 -71856.1787 -143459.0879
87 58271.2104 -71856.1787
88 -64948.5081 58271.2104
89 83938.2662 -64948.5081
90 42507.8170 83938.2662
91 -49580.7887 42507.8170
92 85793.1811 -49580.7887
93 -70344.4744 85793.1811
94 -6814.0289 -70344.4744
95 -15737.3701 -6814.0289
96 147930.4943 -15737.3701
97 213063.0816 147930.4943
98 72826.2490 213063.0816
99 -34459.7025 72826.2490
100 -79216.1575 -34459.7025
101 -79815.2978 -79216.1575
102 -111110.0897 -79815.2978
103 39046.6756 -111110.0897
104 52007.6985 39046.6756
105 -76065.7715 52007.6985
106 -295762.5098 -76065.7715
107 146195.4603 -295762.5098
108 -104872.2276 146195.4603
109 -108290.3137 -104872.2276
110 91917.1948 -108290.3137
111 -238787.2513 91917.1948
112 -101375.4457 -238787.2513
113 -84835.3560 -101375.4457
114 25282.2672 -84835.3560
115 -6420.5109 25282.2672
116 62241.0108 -6420.5109
117 113622.0525 62241.0108
118 -147931.6276 113622.0525
119 66650.5828 -147931.6276
120 -242477.7777 66650.5828
121 24226.5214 -242477.7777
122 -155730.5335 24226.5214
123 103139.2098 -155730.5335
124 -3850.1514 103139.2098
125 -187428.7167 -3850.1514
126 29268.1262 -187428.7167
127 8003.2527 29268.1262
128 1350.2192 8003.2527
129 29268.1262 1350.2192
130 19435.3611 29268.1262
131 91186.7223 19435.3611
132 29268.1262 91186.7223
133 -30144.0315 29268.1262
134 18248.3820 -30144.0315
135 26910.1165 18248.3820
136 21999.4275 26910.1165
137 31090.4241 21999.4275
138 17538.2510 31090.4241
139 29268.1262 17538.2510
140 31858.1740 29268.1262
141 29268.1262 31858.1740
142 -11943.4238 29268.1262
143 21985.8429 -11943.4238
144 16747.6812 21985.8429
145 11303.9014 16747.6812
146 29268.1262 11303.9014
147 17224.4191 29268.1262
148 34642.2003 17224.4191
149 31382.0373 34642.2003
150 -7622.8240 31382.0373
151 -25650.9295 -7622.8240
152 13943.2545 -25650.9295
153 1038.1341 13943.2545
154 27775.4888 1038.1341
155 29268.1262 27775.4888
156 29268.1262 29268.1262
157 25329.7892 29268.1262
158 40865.2598 25329.7892
159 29268.1262 40865.2598
160 29268.1262 29268.1262
161 29268.1262 29268.1262
162 29268.1262 29268.1262
163 36353.2418 29268.1262
164 19251.7140 36353.2418
165 33219.4047 19251.7140
166 29268.1262 33219.4047
167 58320.1911 29268.1262
168 892399.6433 58320.1911
169 29268.1262 892399.6433
170 29268.1262 29268.1262
171 23616.3821 29268.1262
172 29268.1262 23616.3821
173 8373.1483 29268.1262
174 29268.1262 8373.1483
175 21862.0309 29268.1262
176 89432.1141 21862.0309
177 -11439.4325 89432.1141
178 29268.1262 -11439.4325
179 29268.1262 29268.1262
180 57060.1061 29268.1262
181 -18485.6793 57060.1061
182 31673.1233 -18485.6793
183 50247.3889 31673.1233
184 -19977.3842 50247.3889
185 -14838.0182 -19977.3842
186 -16030.6275 -14838.0182
187 15618.3139 -16030.6275
188 -2106.0886 15618.3139
189 29268.1262 -2106.0886
190 13262.5146 29268.1262
191 29268.1262 13262.5146
192 16125.6181 29268.1262
193 29268.1262 16125.6181
194 23372.3899 29268.1262
195 -12895.9151 23372.3899
196 -144048.2459 -12895.9151
197 -12895.9151 -144048.2459
198 -12895.9151 -12895.9151
199 -86204.1344 -12895.9151
200 -12895.9151 -86204.1344
201 -12895.9151 -12895.9151
202 -12895.9151 -12895.9151
203 -23558.2587 -12895.9151
204 -11021.0744 -23558.2587
205 -23779.8406 -11021.0744
206 -7579.1442 -23779.8406
207 -12895.9151 -7579.1442
208 -27880.4863 -12895.9151
209 50231.4022 -27880.4863
210 29268.1262 50231.4022
211 29268.1262 29268.1262
212 59092.2205 29268.1262
213 -1242.1168 59092.2205
214 7927.4901 -1242.1168
215 23138.1775 7927.4901
216 24191.9708 23138.1775
217 50659.9525 24191.9708
218 23138.1775 50659.9525
219 -19460.1087 23138.1775
220 3769.9938 -19460.1087
221 29268.1262 3769.9938
222 29414.7303 29268.1262
223 15392.3051 29414.7303
224 29456.2253 15392.3051
225 345.8293 29456.2253
226 38150.1854 345.8293
227 29268.1262 38150.1854
228 26321.6280 29268.1262
229 29268.1262 26321.6280
230 2843.3897 29268.1262
231 24541.7601 2843.3897
232 29268.1262 24541.7601
233 347.1197 29268.1262
234 29268.1262 347.1197
235 31005.9853 29268.1262
236 29268.1262 31005.9853
237 20243.1219 29268.1262
238 29268.1262 20243.1219
239 29268.1262 29268.1262
240 29268.1262 29268.1262
241 21535.2133 29268.1262
242 29268.1262 21535.2133
243 20500.1432 29268.1262
244 12045.8516 20500.1432
245 -12895.9151 12045.8516
246 -12895.9151 -12895.9151
247 -6412.0358 -12895.9151
248 -12895.9151 -6412.0358
249 -12895.9151 -12895.9151
250 -8320.2368 -12895.9151
251 -12895.9151 -8320.2368
252 21024.5792 -12895.9151
253 4443.0129 21024.5792
254 -95940.8109 4443.0129
255 29268.1262 -95940.8109
256 7405.9368 29268.1262
257 -14375.4965 7405.9368
258 -12895.9151 -14375.4965
259 29268.1262 -12895.9151
260 -12895.9151 29268.1262
261 -12895.9151 -12895.9151
262 39051.9120 -12895.9151
263 34965.7507 39051.9120
264 -12895.9151 34965.7507
265 29268.1262 -12895.9151
266 -12895.9151 29268.1262
267 89547.2459 -12895.9151
268 -12895.9151 89547.2459
269 -22091.7116 -12895.9151
270 -12895.9151 -22091.7116
271 -12895.9151 -12895.9151
272 -12895.9151 -12895.9151
273 -16124.8231 -12895.9151
274 29268.1262 -16124.8231
275 42623.2566 29268.1262
276 13388.0309 42623.2566
277 29268.1262 13388.0309
278 80534.1352 29268.1262
279 44298.5890 80534.1352
280 -12895.9151 44298.5890
281 29268.1262 -12895.9151
282 65466.0576 29268.1262
283 29268.1262 65466.0576
284 37912.3754 29268.1262
285 29268.1262 37912.3754
286 -16500.7829 29268.1262
287 29634.3242 -16500.7829
288 -12895.9151 29634.3242
289 68645.4134 -12895.9151
290 -12895.9151 68645.4134
291 39801.4098 -12895.9151
292 54131.1182 39801.4098
293 29268.1262 54131.1182
294 7535.8105 29268.1262
295 14006.3083 7535.8105
296 29268.1262 14006.3083
297 -12895.9151 29268.1262
298 29268.1262 -12895.9151
299 1683.3572 29268.1262
300 -21170.2426 1683.3572
301 10134.3494 -21170.2426
302 26203.1519 10134.3494
303 -12895.9151 26203.1519
304 97098.1523 -12895.9151
305 7294.2924 97098.1523
306 -12895.9151 7294.2924
307 -50240.2563 -12895.9151
308 1249.7813 -50240.2563
309 -12895.9151 1249.7813
310 -12895.9151 -12895.9151
311 70339.2670 -12895.9151
312 -12895.9151 70339.2670
313 -35175.0997 -12895.9151
314 11587.0227 -35175.0997
315 64018.5109 11587.0227
316 29631.8963 64018.5109
317 -27343.7769 29631.8963
318 23796.4518 -27343.7769
319 29268.1262 23796.4518
320 29268.1262 29268.1262
321 29268.1262 29268.1262
322 29268.1262 29268.1262
323 -12895.9151 29268.1262
324 29268.1262 -12895.9151
325 26308.4131 29268.1262
326 -12895.9151 26308.4131
327 -12895.9151 -12895.9151
328 29268.1262 -12895.9151
329 29268.1262 29268.1262
330 29268.1262 29268.1262
331 29268.1262 29268.1262
332 29268.1262 29268.1262
333 22065.9821 29268.1262
334 29492.8409 22065.9821
335 26204.8240 29492.8409
336 141959.3090 26204.8240
337 29268.1262 141959.3090
338 26082.9678 29268.1262
339 20224.3425 26082.9678
340 -6185.1728 20224.3425
341 29268.1262 -6185.1728
342 29268.1262 29268.1262
343 29268.1262 29268.1262
344 104371.9698 29268.1262
345 29268.1262 104371.9698
346 15803.2182 29268.1262
347 29268.1262 15803.2182
348 29268.1262 29268.1262
349 26480.6749 29268.1262
350 41319.2054 26480.6749
351 22261.8646 41319.2054
352 20005.2039 22261.8646
353 27164.0740 20005.2039
354 29268.1262 27164.0740
355 29268.1262 29268.1262
356 62525.8537 29268.1262
357 29268.1262 62525.8537
358 29268.1262 29268.1262
359 29268.1262 29268.1262
360 23597.0153 29268.1262
361 29268.1262 23597.0153
362 29268.1262 29268.1262
363 -67943.7223 29268.1262
364 17837.3355 -67943.7223
365 29268.1262 17837.3355
366 29268.1262 29268.1262
367 29268.1262 29268.1262
368 29268.1262 29268.1262
369 62148.9428 29268.1262
370 29268.1262 62148.9428
371 75442.7081 29268.1262
372 47466.9676 75442.7081
373 13641.2611 47466.9676
374 29268.1262 13641.2611
375 29268.1262 29268.1262
376 12620.8166 29268.1262
377 24221.5966 12620.8166
378 10112.3784 24221.5966
379 29268.1262 10112.3784
380 29268.1262 29268.1262
381 61189.8980 29268.1262
382 91110.8221 61189.8980
383 29268.1262 91110.8221
384 20271.8496 29268.1262
385 -106966.4255 20271.8496
386 39288.3148 -106966.4255
387 3777.1431 39288.3148
388 22655.9444 3777.1431
389 14935.5040 22655.9444
390 17246.8195 14935.5040
391 -6152.4952 17246.8195
392 -44485.5778 -6152.4952
393 69803.8539 -44485.5778
394 -91897.9029 69803.8539
395 -14362.6556 -91897.9029
396 -120978.1743 -14362.6556
397 -539798.0202 -120978.1743
398 -2599.6976 -539798.0202
399 48066.8103 -2599.6976
400 7051.7420 48066.8103
401 28254.8104 7051.7420
402 -3696.2117 28254.8104
403 -92768.8288 -3696.2117
404 -25845.1997 -92768.8288
405 15471.1292 -25845.1997
406 17005.1009 15471.1292
407 -14854.1938 17005.1009
408 -4598.1795 -14854.1938
409 -63801.2624 -4598.1795
410 -73896.2632 -63801.2624
411 -156603.4022 -73896.2632
412 -93868.7126 -156603.4022
413 29585.1004 -93868.7126
414 -139508.1245 29585.1004
415 -89769.1764 -139508.1245
416 128309.4869 -89769.1764
417 -86655.2669 128309.4869
418 -58399.6203 -86655.2669
419 -44237.0641 -58399.6203
420 -41702.0976 -44237.0641
421 -209029.6603 -41702.0976
422 -617473.3856 -209029.6603
423 -439248.1597 -617473.3856
424 -67953.7316 -439248.1597
425 -280137.2562 -67953.7316
426 -315144.1159 -280137.2562
427 -1188447.8307 -315144.1159
428 -186944.8306 -1188447.8307
429 -260412.3078 -186944.8306
430 -375551.4045 -260412.3078
431 NA -375551.4045
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2617776.3344 131691.2321
[2,] 1343261.9381 2617776.3344
[3,] -744280.2782 1343261.9381
[4,] -384882.4697 -744280.2782
[5,] 231829.6052 -384882.4697
[6,] -24799.1285 231829.6052
[7,] -140602.1496 -24799.1285
[8,] 171701.7832 -140602.1496
[9,] 687631.7910 171701.7832
[10,] 24217.3337 687631.7910
[11,] 344779.5263 24217.3337
[12,] 240841.5619 344779.5263
[13,] -741217.9665 240841.5619
[14,] -52513.3362 -741217.9665
[15,] 708793.0810 -52513.3362
[16,] -188311.0806 708793.0810
[17,] -167054.1605 -188311.0806
[18,] -86916.7805 -167054.1605
[19,] -703145.6635 -86916.7805
[20,] 285123.3410 -703145.6635
[21,] -15121.8885 285123.3410
[22,] -47851.8834 -15121.8885
[23,] -818053.2764 -47851.8834
[24,] -264218.2790 -818053.2764
[25,] 8784.4069 -264218.2790
[26,] -362208.6164 8784.4069
[27,] -60845.5727 -362208.6164
[28,] 338692.8381 -60845.5727
[29,] 182219.7224 338692.8381
[30,] 281788.1807 182219.7224
[31,] -36995.6012 281788.1807
[32,] 263445.7065 -36995.6012
[33,] 145992.6598 263445.7065
[34,] 373117.5424 145992.6598
[35,] -651360.7695 373117.5424
[36,] 102982.2074 -651360.7695
[37,] 137960.6316 102982.2074
[38,] -42632.1754 137960.6316
[39,] -336872.4361 -42632.1754
[40,] -219519.1761 -336872.4361
[41,] -222531.9421 -219519.1761
[42,] 137566.0210 -222531.9421
[43,] -172492.4321 137566.0210
[44,] 14357.7386 -172492.4321
[45,] -30779.2236 14357.7386
[46,] 405337.5305 -30779.2236
[47,] -321471.9194 405337.5305
[48,] 18934.8649 -321471.9194
[49,] -195140.8815 18934.8649
[50,] -3391.6077 -195140.8815
[51,] -192350.0646 -3391.6077
[52,] -234497.9226 -192350.0646
[53,] 60187.4370 -234497.9226
[54,] -462598.5773 60187.4370
[55,] 9233.5964 -462598.5773
[56,] 306052.0936 9233.5964
[57,] 33188.9153 306052.0936
[58,] 5900.7235 33188.9153
[59,] -198482.0380 5900.7235
[60,] 23938.1649 -198482.0380
[61,] -196629.2342 23938.1649
[62,] 20928.2372 -196629.2342
[63,] -179696.4179 20928.2372
[64,] -2681.0677 -179696.4179
[65,] 185234.2138 -2681.0677
[66,] 23891.0866 185234.2138
[67,] -21817.5531 23891.0866
[68,] -76648.7399 -21817.5531
[69,] -138163.7547 -76648.7399
[70,] -208938.8859 -138163.7547
[71,] -9289.8180 -208938.8859
[72,] -47672.1817 -9289.8180
[73,] 29587.7757 -47672.1817
[74,] 31861.7253 29587.7757
[75,] -126419.0457 31861.7253
[76,] 42686.2541 -126419.0457
[77,] -26071.1896 42686.2541
[78,] 28375.2746 -26071.1896
[79,] 103544.1770 28375.2746
[80,] -4214.0980 103544.1770
[81,] -65473.3759 -4214.0980
[82,] 14124.1325 -65473.3759
[83,] -17642.2383 14124.1325
[84,] 76278.7647 -17642.2383
[85,] -143459.0879 76278.7647
[86,] -71856.1787 -143459.0879
[87,] 58271.2104 -71856.1787
[88,] -64948.5081 58271.2104
[89,] 83938.2662 -64948.5081
[90,] 42507.8170 83938.2662
[91,] -49580.7887 42507.8170
[92,] 85793.1811 -49580.7887
[93,] -70344.4744 85793.1811
[94,] -6814.0289 -70344.4744
[95,] -15737.3701 -6814.0289
[96,] 147930.4943 -15737.3701
[97,] 213063.0816 147930.4943
[98,] 72826.2490 213063.0816
[99,] -34459.7025 72826.2490
[100,] -79216.1575 -34459.7025
[101,] -79815.2978 -79216.1575
[102,] -111110.0897 -79815.2978
[103,] 39046.6756 -111110.0897
[104,] 52007.6985 39046.6756
[105,] -76065.7715 52007.6985
[106,] -295762.5098 -76065.7715
[107,] 146195.4603 -295762.5098
[108,] -104872.2276 146195.4603
[109,] -108290.3137 -104872.2276
[110,] 91917.1948 -108290.3137
[111,] -238787.2513 91917.1948
[112,] -101375.4457 -238787.2513
[113,] -84835.3560 -101375.4457
[114,] 25282.2672 -84835.3560
[115,] -6420.5109 25282.2672
[116,] 62241.0108 -6420.5109
[117,] 113622.0525 62241.0108
[118,] -147931.6276 113622.0525
[119,] 66650.5828 -147931.6276
[120,] -242477.7777 66650.5828
[121,] 24226.5214 -242477.7777
[122,] -155730.5335 24226.5214
[123,] 103139.2098 -155730.5335
[124,] -3850.1514 103139.2098
[125,] -187428.7167 -3850.1514
[126,] 29268.1262 -187428.7167
[127,] 8003.2527 29268.1262
[128,] 1350.2192 8003.2527
[129,] 29268.1262 1350.2192
[130,] 19435.3611 29268.1262
[131,] 91186.7223 19435.3611
[132,] 29268.1262 91186.7223
[133,] -30144.0315 29268.1262
[134,] 18248.3820 -30144.0315
[135,] 26910.1165 18248.3820
[136,] 21999.4275 26910.1165
[137,] 31090.4241 21999.4275
[138,] 17538.2510 31090.4241
[139,] 29268.1262 17538.2510
[140,] 31858.1740 29268.1262
[141,] 29268.1262 31858.1740
[142,] -11943.4238 29268.1262
[143,] 21985.8429 -11943.4238
[144,] 16747.6812 21985.8429
[145,] 11303.9014 16747.6812
[146,] 29268.1262 11303.9014
[147,] 17224.4191 29268.1262
[148,] 34642.2003 17224.4191
[149,] 31382.0373 34642.2003
[150,] -7622.8240 31382.0373
[151,] -25650.9295 -7622.8240
[152,] 13943.2545 -25650.9295
[153,] 1038.1341 13943.2545
[154,] 27775.4888 1038.1341
[155,] 29268.1262 27775.4888
[156,] 29268.1262 29268.1262
[157,] 25329.7892 29268.1262
[158,] 40865.2598 25329.7892
[159,] 29268.1262 40865.2598
[160,] 29268.1262 29268.1262
[161,] 29268.1262 29268.1262
[162,] 29268.1262 29268.1262
[163,] 36353.2418 29268.1262
[164,] 19251.7140 36353.2418
[165,] 33219.4047 19251.7140
[166,] 29268.1262 33219.4047
[167,] 58320.1911 29268.1262
[168,] 892399.6433 58320.1911
[169,] 29268.1262 892399.6433
[170,] 29268.1262 29268.1262
[171,] 23616.3821 29268.1262
[172,] 29268.1262 23616.3821
[173,] 8373.1483 29268.1262
[174,] 29268.1262 8373.1483
[175,] 21862.0309 29268.1262
[176,] 89432.1141 21862.0309
[177,] -11439.4325 89432.1141
[178,] 29268.1262 -11439.4325
[179,] 29268.1262 29268.1262
[180,] 57060.1061 29268.1262
[181,] -18485.6793 57060.1061
[182,] 31673.1233 -18485.6793
[183,] 50247.3889 31673.1233
[184,] -19977.3842 50247.3889
[185,] -14838.0182 -19977.3842
[186,] -16030.6275 -14838.0182
[187,] 15618.3139 -16030.6275
[188,] -2106.0886 15618.3139
[189,] 29268.1262 -2106.0886
[190,] 13262.5146 29268.1262
[191,] 29268.1262 13262.5146
[192,] 16125.6181 29268.1262
[193,] 29268.1262 16125.6181
[194,] 23372.3899 29268.1262
[195,] -12895.9151 23372.3899
[196,] -144048.2459 -12895.9151
[197,] -12895.9151 -144048.2459
[198,] -12895.9151 -12895.9151
[199,] -86204.1344 -12895.9151
[200,] -12895.9151 -86204.1344
[201,] -12895.9151 -12895.9151
[202,] -12895.9151 -12895.9151
[203,] -23558.2587 -12895.9151
[204,] -11021.0744 -23558.2587
[205,] -23779.8406 -11021.0744
[206,] -7579.1442 -23779.8406
[207,] -12895.9151 -7579.1442
[208,] -27880.4863 -12895.9151
[209,] 50231.4022 -27880.4863
[210,] 29268.1262 50231.4022
[211,] 29268.1262 29268.1262
[212,] 59092.2205 29268.1262
[213,] -1242.1168 59092.2205
[214,] 7927.4901 -1242.1168
[215,] 23138.1775 7927.4901
[216,] 24191.9708 23138.1775
[217,] 50659.9525 24191.9708
[218,] 23138.1775 50659.9525
[219,] -19460.1087 23138.1775
[220,] 3769.9938 -19460.1087
[221,] 29268.1262 3769.9938
[222,] 29414.7303 29268.1262
[223,] 15392.3051 29414.7303
[224,] 29456.2253 15392.3051
[225,] 345.8293 29456.2253
[226,] 38150.1854 345.8293
[227,] 29268.1262 38150.1854
[228,] 26321.6280 29268.1262
[229,] 29268.1262 26321.6280
[230,] 2843.3897 29268.1262
[231,] 24541.7601 2843.3897
[232,] 29268.1262 24541.7601
[233,] 347.1197 29268.1262
[234,] 29268.1262 347.1197
[235,] 31005.9853 29268.1262
[236,] 29268.1262 31005.9853
[237,] 20243.1219 29268.1262
[238,] 29268.1262 20243.1219
[239,] 29268.1262 29268.1262
[240,] 29268.1262 29268.1262
[241,] 21535.2133 29268.1262
[242,] 29268.1262 21535.2133
[243,] 20500.1432 29268.1262
[244,] 12045.8516 20500.1432
[245,] -12895.9151 12045.8516
[246,] -12895.9151 -12895.9151
[247,] -6412.0358 -12895.9151
[248,] -12895.9151 -6412.0358
[249,] -12895.9151 -12895.9151
[250,] -8320.2368 -12895.9151
[251,] -12895.9151 -8320.2368
[252,] 21024.5792 -12895.9151
[253,] 4443.0129 21024.5792
[254,] -95940.8109 4443.0129
[255,] 29268.1262 -95940.8109
[256,] 7405.9368 29268.1262
[257,] -14375.4965 7405.9368
[258,] -12895.9151 -14375.4965
[259,] 29268.1262 -12895.9151
[260,] -12895.9151 29268.1262
[261,] -12895.9151 -12895.9151
[262,] 39051.9120 -12895.9151
[263,] 34965.7507 39051.9120
[264,] -12895.9151 34965.7507
[265,] 29268.1262 -12895.9151
[266,] -12895.9151 29268.1262
[267,] 89547.2459 -12895.9151
[268,] -12895.9151 89547.2459
[269,] -22091.7116 -12895.9151
[270,] -12895.9151 -22091.7116
[271,] -12895.9151 -12895.9151
[272,] -12895.9151 -12895.9151
[273,] -16124.8231 -12895.9151
[274,] 29268.1262 -16124.8231
[275,] 42623.2566 29268.1262
[276,] 13388.0309 42623.2566
[277,] 29268.1262 13388.0309
[278,] 80534.1352 29268.1262
[279,] 44298.5890 80534.1352
[280,] -12895.9151 44298.5890
[281,] 29268.1262 -12895.9151
[282,] 65466.0576 29268.1262
[283,] 29268.1262 65466.0576
[284,] 37912.3754 29268.1262
[285,] 29268.1262 37912.3754
[286,] -16500.7829 29268.1262
[287,] 29634.3242 -16500.7829
[288,] -12895.9151 29634.3242
[289,] 68645.4134 -12895.9151
[290,] -12895.9151 68645.4134
[291,] 39801.4098 -12895.9151
[292,] 54131.1182 39801.4098
[293,] 29268.1262 54131.1182
[294,] 7535.8105 29268.1262
[295,] 14006.3083 7535.8105
[296,] 29268.1262 14006.3083
[297,] -12895.9151 29268.1262
[298,] 29268.1262 -12895.9151
[299,] 1683.3572 29268.1262
[300,] -21170.2426 1683.3572
[301,] 10134.3494 -21170.2426
[302,] 26203.1519 10134.3494
[303,] -12895.9151 26203.1519
[304,] 97098.1523 -12895.9151
[305,] 7294.2924 97098.1523
[306,] -12895.9151 7294.2924
[307,] -50240.2563 -12895.9151
[308,] 1249.7813 -50240.2563
[309,] -12895.9151 1249.7813
[310,] -12895.9151 -12895.9151
[311,] 70339.2670 -12895.9151
[312,] -12895.9151 70339.2670
[313,] -35175.0997 -12895.9151
[314,] 11587.0227 -35175.0997
[315,] 64018.5109 11587.0227
[316,] 29631.8963 64018.5109
[317,] -27343.7769 29631.8963
[318,] 23796.4518 -27343.7769
[319,] 29268.1262 23796.4518
[320,] 29268.1262 29268.1262
[321,] 29268.1262 29268.1262
[322,] 29268.1262 29268.1262
[323,] -12895.9151 29268.1262
[324,] 29268.1262 -12895.9151
[325,] 26308.4131 29268.1262
[326,] -12895.9151 26308.4131
[327,] -12895.9151 -12895.9151
[328,] 29268.1262 -12895.9151
[329,] 29268.1262 29268.1262
[330,] 29268.1262 29268.1262
[331,] 29268.1262 29268.1262
[332,] 29268.1262 29268.1262
[333,] 22065.9821 29268.1262
[334,] 29492.8409 22065.9821
[335,] 26204.8240 29492.8409
[336,] 141959.3090 26204.8240
[337,] 29268.1262 141959.3090
[338,] 26082.9678 29268.1262
[339,] 20224.3425 26082.9678
[340,] -6185.1728 20224.3425
[341,] 29268.1262 -6185.1728
[342,] 29268.1262 29268.1262
[343,] 29268.1262 29268.1262
[344,] 104371.9698 29268.1262
[345,] 29268.1262 104371.9698
[346,] 15803.2182 29268.1262
[347,] 29268.1262 15803.2182
[348,] 29268.1262 29268.1262
[349,] 26480.6749 29268.1262
[350,] 41319.2054 26480.6749
[351,] 22261.8646 41319.2054
[352,] 20005.2039 22261.8646
[353,] 27164.0740 20005.2039
[354,] 29268.1262 27164.0740
[355,] 29268.1262 29268.1262
[356,] 62525.8537 29268.1262
[357,] 29268.1262 62525.8537
[358,] 29268.1262 29268.1262
[359,] 29268.1262 29268.1262
[360,] 23597.0153 29268.1262
[361,] 29268.1262 23597.0153
[362,] 29268.1262 29268.1262
[363,] -67943.7223 29268.1262
[364,] 17837.3355 -67943.7223
[365,] 29268.1262 17837.3355
[366,] 29268.1262 29268.1262
[367,] 29268.1262 29268.1262
[368,] 29268.1262 29268.1262
[369,] 62148.9428 29268.1262
[370,] 29268.1262 62148.9428
[371,] 75442.7081 29268.1262
[372,] 47466.9676 75442.7081
[373,] 13641.2611 47466.9676
[374,] 29268.1262 13641.2611
[375,] 29268.1262 29268.1262
[376,] 12620.8166 29268.1262
[377,] 24221.5966 12620.8166
[378,] 10112.3784 24221.5966
[379,] 29268.1262 10112.3784
[380,] 29268.1262 29268.1262
[381,] 61189.8980 29268.1262
[382,] 91110.8221 61189.8980
[383,] 29268.1262 91110.8221
[384,] 20271.8496 29268.1262
[385,] -106966.4255 20271.8496
[386,] 39288.3148 -106966.4255
[387,] 3777.1431 39288.3148
[388,] 22655.9444 3777.1431
[389,] 14935.5040 22655.9444
[390,] 17246.8195 14935.5040
[391,] -6152.4952 17246.8195
[392,] -44485.5778 -6152.4952
[393,] 69803.8539 -44485.5778
[394,] -91897.9029 69803.8539
[395,] -14362.6556 -91897.9029
[396,] -120978.1743 -14362.6556
[397,] -539798.0202 -120978.1743
[398,] -2599.6976 -539798.0202
[399,] 48066.8103 -2599.6976
[400,] 7051.7420 48066.8103
[401,] 28254.8104 7051.7420
[402,] -3696.2117 28254.8104
[403,] -92768.8288 -3696.2117
[404,] -25845.1997 -92768.8288
[405,] 15471.1292 -25845.1997
[406,] 17005.1009 15471.1292
[407,] -14854.1938 17005.1009
[408,] -4598.1795 -14854.1938
[409,] -63801.2624 -4598.1795
[410,] -73896.2632 -63801.2624
[411,] -156603.4022 -73896.2632
[412,] -93868.7126 -156603.4022
[413,] 29585.1004 -93868.7126
[414,] -139508.1245 29585.1004
[415,] -89769.1764 -139508.1245
[416,] 128309.4869 -89769.1764
[417,] -86655.2669 128309.4869
[418,] -58399.6203 -86655.2669
[419,] -44237.0641 -58399.6203
[420,] -41702.0976 -44237.0641
[421,] -209029.6603 -41702.0976
[422,] -617473.3856 -209029.6603
[423,] -439248.1597 -617473.3856
[424,] -67953.7316 -439248.1597
[425,] -280137.2562 -67953.7316
[426,] -315144.1159 -280137.2562
[427,] -1188447.8307 -315144.1159
[428,] -186944.8306 -1188447.8307
[429,] -260412.3078 -186944.8306
[430,] -375551.4045 -260412.3078
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2617776.3344 131691.2321
2 1343261.9381 2617776.3344
3 -744280.2782 1343261.9381
4 -384882.4697 -744280.2782
5 231829.6052 -384882.4697
6 -24799.1285 231829.6052
7 -140602.1496 -24799.1285
8 171701.7832 -140602.1496
9 687631.7910 171701.7832
10 24217.3337 687631.7910
11 344779.5263 24217.3337
12 240841.5619 344779.5263
13 -741217.9665 240841.5619
14 -52513.3362 -741217.9665
15 708793.0810 -52513.3362
16 -188311.0806 708793.0810
17 -167054.1605 -188311.0806
18 -86916.7805 -167054.1605
19 -703145.6635 -86916.7805
20 285123.3410 -703145.6635
21 -15121.8885 285123.3410
22 -47851.8834 -15121.8885
23 -818053.2764 -47851.8834
24 -264218.2790 -818053.2764
25 8784.4069 -264218.2790
26 -362208.6164 8784.4069
27 -60845.5727 -362208.6164
28 338692.8381 -60845.5727
29 182219.7224 338692.8381
30 281788.1807 182219.7224
31 -36995.6012 281788.1807
32 263445.7065 -36995.6012
33 145992.6598 263445.7065
34 373117.5424 145992.6598
35 -651360.7695 373117.5424
36 102982.2074 -651360.7695
37 137960.6316 102982.2074
38 -42632.1754 137960.6316
39 -336872.4361 -42632.1754
40 -219519.1761 -336872.4361
41 -222531.9421 -219519.1761
42 137566.0210 -222531.9421
43 -172492.4321 137566.0210
44 14357.7386 -172492.4321
45 -30779.2236 14357.7386
46 405337.5305 -30779.2236
47 -321471.9194 405337.5305
48 18934.8649 -321471.9194
49 -195140.8815 18934.8649
50 -3391.6077 -195140.8815
51 -192350.0646 -3391.6077
52 -234497.9226 -192350.0646
53 60187.4370 -234497.9226
54 -462598.5773 60187.4370
55 9233.5964 -462598.5773
56 306052.0936 9233.5964
57 33188.9153 306052.0936
58 5900.7235 33188.9153
59 -198482.0380 5900.7235
60 23938.1649 -198482.0380
61 -196629.2342 23938.1649
62 20928.2372 -196629.2342
63 -179696.4179 20928.2372
64 -2681.0677 -179696.4179
65 185234.2138 -2681.0677
66 23891.0866 185234.2138
67 -21817.5531 23891.0866
68 -76648.7399 -21817.5531
69 -138163.7547 -76648.7399
70 -208938.8859 -138163.7547
71 -9289.8180 -208938.8859
72 -47672.1817 -9289.8180
73 29587.7757 -47672.1817
74 31861.7253 29587.7757
75 -126419.0457 31861.7253
76 42686.2541 -126419.0457
77 -26071.1896 42686.2541
78 28375.2746 -26071.1896
79 103544.1770 28375.2746
80 -4214.0980 103544.1770
81 -65473.3759 -4214.0980
82 14124.1325 -65473.3759
83 -17642.2383 14124.1325
84 76278.7647 -17642.2383
85 -143459.0879 76278.7647
86 -71856.1787 -143459.0879
87 58271.2104 -71856.1787
88 -64948.5081 58271.2104
89 83938.2662 -64948.5081
90 42507.8170 83938.2662
91 -49580.7887 42507.8170
92 85793.1811 -49580.7887
93 -70344.4744 85793.1811
94 -6814.0289 -70344.4744
95 -15737.3701 -6814.0289
96 147930.4943 -15737.3701
97 213063.0816 147930.4943
98 72826.2490 213063.0816
99 -34459.7025 72826.2490
100 -79216.1575 -34459.7025
101 -79815.2978 -79216.1575
102 -111110.0897 -79815.2978
103 39046.6756 -111110.0897
104 52007.6985 39046.6756
105 -76065.7715 52007.6985
106 -295762.5098 -76065.7715
107 146195.4603 -295762.5098
108 -104872.2276 146195.4603
109 -108290.3137 -104872.2276
110 91917.1948 -108290.3137
111 -238787.2513 91917.1948
112 -101375.4457 -238787.2513
113 -84835.3560 -101375.4457
114 25282.2672 -84835.3560
115 -6420.5109 25282.2672
116 62241.0108 -6420.5109
117 113622.0525 62241.0108
118 -147931.6276 113622.0525
119 66650.5828 -147931.6276
120 -242477.7777 66650.5828
121 24226.5214 -242477.7777
122 -155730.5335 24226.5214
123 103139.2098 -155730.5335
124 -3850.1514 103139.2098
125 -187428.7167 -3850.1514
126 29268.1262 -187428.7167
127 8003.2527 29268.1262
128 1350.2192 8003.2527
129 29268.1262 1350.2192
130 19435.3611 29268.1262
131 91186.7223 19435.3611
132 29268.1262 91186.7223
133 -30144.0315 29268.1262
134 18248.3820 -30144.0315
135 26910.1165 18248.3820
136 21999.4275 26910.1165
137 31090.4241 21999.4275
138 17538.2510 31090.4241
139 29268.1262 17538.2510
140 31858.1740 29268.1262
141 29268.1262 31858.1740
142 -11943.4238 29268.1262
143 21985.8429 -11943.4238
144 16747.6812 21985.8429
145 11303.9014 16747.6812
146 29268.1262 11303.9014
147 17224.4191 29268.1262
148 34642.2003 17224.4191
149 31382.0373 34642.2003
150 -7622.8240 31382.0373
151 -25650.9295 -7622.8240
152 13943.2545 -25650.9295
153 1038.1341 13943.2545
154 27775.4888 1038.1341
155 29268.1262 27775.4888
156 29268.1262 29268.1262
157 25329.7892 29268.1262
158 40865.2598 25329.7892
159 29268.1262 40865.2598
160 29268.1262 29268.1262
161 29268.1262 29268.1262
162 29268.1262 29268.1262
163 36353.2418 29268.1262
164 19251.7140 36353.2418
165 33219.4047 19251.7140
166 29268.1262 33219.4047
167 58320.1911 29268.1262
168 892399.6433 58320.1911
169 29268.1262 892399.6433
170 29268.1262 29268.1262
171 23616.3821 29268.1262
172 29268.1262 23616.3821
173 8373.1483 29268.1262
174 29268.1262 8373.1483
175 21862.0309 29268.1262
176 89432.1141 21862.0309
177 -11439.4325 89432.1141
178 29268.1262 -11439.4325
179 29268.1262 29268.1262
180 57060.1061 29268.1262
181 -18485.6793 57060.1061
182 31673.1233 -18485.6793
183 50247.3889 31673.1233
184 -19977.3842 50247.3889
185 -14838.0182 -19977.3842
186 -16030.6275 -14838.0182
187 15618.3139 -16030.6275
188 -2106.0886 15618.3139
189 29268.1262 -2106.0886
190 13262.5146 29268.1262
191 29268.1262 13262.5146
192 16125.6181 29268.1262
193 29268.1262 16125.6181
194 23372.3899 29268.1262
195 -12895.9151 23372.3899
196 -144048.2459 -12895.9151
197 -12895.9151 -144048.2459
198 -12895.9151 -12895.9151
199 -86204.1344 -12895.9151
200 -12895.9151 -86204.1344
201 -12895.9151 -12895.9151
202 -12895.9151 -12895.9151
203 -23558.2587 -12895.9151
204 -11021.0744 -23558.2587
205 -23779.8406 -11021.0744
206 -7579.1442 -23779.8406
207 -12895.9151 -7579.1442
208 -27880.4863 -12895.9151
209 50231.4022 -27880.4863
210 29268.1262 50231.4022
211 29268.1262 29268.1262
212 59092.2205 29268.1262
213 -1242.1168 59092.2205
214 7927.4901 -1242.1168
215 23138.1775 7927.4901
216 24191.9708 23138.1775
217 50659.9525 24191.9708
218 23138.1775 50659.9525
219 -19460.1087 23138.1775
220 3769.9938 -19460.1087
221 29268.1262 3769.9938
222 29414.7303 29268.1262
223 15392.3051 29414.7303
224 29456.2253 15392.3051
225 345.8293 29456.2253
226 38150.1854 345.8293
227 29268.1262 38150.1854
228 26321.6280 29268.1262
229 29268.1262 26321.6280
230 2843.3897 29268.1262
231 24541.7601 2843.3897
232 29268.1262 24541.7601
233 347.1197 29268.1262
234 29268.1262 347.1197
235 31005.9853 29268.1262
236 29268.1262 31005.9853
237 20243.1219 29268.1262
238 29268.1262 20243.1219
239 29268.1262 29268.1262
240 29268.1262 29268.1262
241 21535.2133 29268.1262
242 29268.1262 21535.2133
243 20500.1432 29268.1262
244 12045.8516 20500.1432
245 -12895.9151 12045.8516
246 -12895.9151 -12895.9151
247 -6412.0358 -12895.9151
248 -12895.9151 -6412.0358
249 -12895.9151 -12895.9151
250 -8320.2368 -12895.9151
251 -12895.9151 -8320.2368
252 21024.5792 -12895.9151
253 4443.0129 21024.5792
254 -95940.8109 4443.0129
255 29268.1262 -95940.8109
256 7405.9368 29268.1262
257 -14375.4965 7405.9368
258 -12895.9151 -14375.4965
259 29268.1262 -12895.9151
260 -12895.9151 29268.1262
261 -12895.9151 -12895.9151
262 39051.9120 -12895.9151
263 34965.7507 39051.9120
264 -12895.9151 34965.7507
265 29268.1262 -12895.9151
266 -12895.9151 29268.1262
267 89547.2459 -12895.9151
268 -12895.9151 89547.2459
269 -22091.7116 -12895.9151
270 -12895.9151 -22091.7116
271 -12895.9151 -12895.9151
272 -12895.9151 -12895.9151
273 -16124.8231 -12895.9151
274 29268.1262 -16124.8231
275 42623.2566 29268.1262
276 13388.0309 42623.2566
277 29268.1262 13388.0309
278 80534.1352 29268.1262
279 44298.5890 80534.1352
280 -12895.9151 44298.5890
281 29268.1262 -12895.9151
282 65466.0576 29268.1262
283 29268.1262 65466.0576
284 37912.3754 29268.1262
285 29268.1262 37912.3754
286 -16500.7829 29268.1262
287 29634.3242 -16500.7829
288 -12895.9151 29634.3242
289 68645.4134 -12895.9151
290 -12895.9151 68645.4134
291 39801.4098 -12895.9151
292 54131.1182 39801.4098
293 29268.1262 54131.1182
294 7535.8105 29268.1262
295 14006.3083 7535.8105
296 29268.1262 14006.3083
297 -12895.9151 29268.1262
298 29268.1262 -12895.9151
299 1683.3572 29268.1262
300 -21170.2426 1683.3572
301 10134.3494 -21170.2426
302 26203.1519 10134.3494
303 -12895.9151 26203.1519
304 97098.1523 -12895.9151
305 7294.2924 97098.1523
306 -12895.9151 7294.2924
307 -50240.2563 -12895.9151
308 1249.7813 -50240.2563
309 -12895.9151 1249.7813
310 -12895.9151 -12895.9151
311 70339.2670 -12895.9151
312 -12895.9151 70339.2670
313 -35175.0997 -12895.9151
314 11587.0227 -35175.0997
315 64018.5109 11587.0227
316 29631.8963 64018.5109
317 -27343.7769 29631.8963
318 23796.4518 -27343.7769
319 29268.1262 23796.4518
320 29268.1262 29268.1262
321 29268.1262 29268.1262
322 29268.1262 29268.1262
323 -12895.9151 29268.1262
324 29268.1262 -12895.9151
325 26308.4131 29268.1262
326 -12895.9151 26308.4131
327 -12895.9151 -12895.9151
328 29268.1262 -12895.9151
329 29268.1262 29268.1262
330 29268.1262 29268.1262
331 29268.1262 29268.1262
332 29268.1262 29268.1262
333 22065.9821 29268.1262
334 29492.8409 22065.9821
335 26204.8240 29492.8409
336 141959.3090 26204.8240
337 29268.1262 141959.3090
338 26082.9678 29268.1262
339 20224.3425 26082.9678
340 -6185.1728 20224.3425
341 29268.1262 -6185.1728
342 29268.1262 29268.1262
343 29268.1262 29268.1262
344 104371.9698 29268.1262
345 29268.1262 104371.9698
346 15803.2182 29268.1262
347 29268.1262 15803.2182
348 29268.1262 29268.1262
349 26480.6749 29268.1262
350 41319.2054 26480.6749
351 22261.8646 41319.2054
352 20005.2039 22261.8646
353 27164.0740 20005.2039
354 29268.1262 27164.0740
355 29268.1262 29268.1262
356 62525.8537 29268.1262
357 29268.1262 62525.8537
358 29268.1262 29268.1262
359 29268.1262 29268.1262
360 23597.0153 29268.1262
361 29268.1262 23597.0153
362 29268.1262 29268.1262
363 -67943.7223 29268.1262
364 17837.3355 -67943.7223
365 29268.1262 17837.3355
366 29268.1262 29268.1262
367 29268.1262 29268.1262
368 29268.1262 29268.1262
369 62148.9428 29268.1262
370 29268.1262 62148.9428
371 75442.7081 29268.1262
372 47466.9676 75442.7081
373 13641.2611 47466.9676
374 29268.1262 13641.2611
375 29268.1262 29268.1262
376 12620.8166 29268.1262
377 24221.5966 12620.8166
378 10112.3784 24221.5966
379 29268.1262 10112.3784
380 29268.1262 29268.1262
381 61189.8980 29268.1262
382 91110.8221 61189.8980
383 29268.1262 91110.8221
384 20271.8496 29268.1262
385 -106966.4255 20271.8496
386 39288.3148 -106966.4255
387 3777.1431 39288.3148
388 22655.9444 3777.1431
389 14935.5040 22655.9444
390 17246.8195 14935.5040
391 -6152.4952 17246.8195
392 -44485.5778 -6152.4952
393 69803.8539 -44485.5778
394 -91897.9029 69803.8539
395 -14362.6556 -91897.9029
396 -120978.1743 -14362.6556
397 -539798.0202 -120978.1743
398 -2599.6976 -539798.0202
399 48066.8103 -2599.6976
400 7051.7420 48066.8103
401 28254.8104 7051.7420
402 -3696.2117 28254.8104
403 -92768.8288 -3696.2117
404 -25845.1997 -92768.8288
405 15471.1292 -25845.1997
406 17005.1009 15471.1292
407 -14854.1938 17005.1009
408 -4598.1795 -14854.1938
409 -63801.2624 -4598.1795
410 -73896.2632 -63801.2624
411 -156603.4022 -73896.2632
412 -93868.7126 -156603.4022
413 29585.1004 -93868.7126
414 -139508.1245 29585.1004
415 -89769.1764 -139508.1245
416 128309.4869 -89769.1764
417 -86655.2669 128309.4869
418 -58399.6203 -86655.2669
419 -44237.0641 -58399.6203
420 -41702.0976 -44237.0641
421 -209029.6603 -41702.0976
422 -617473.3856 -209029.6603
423 -439248.1597 -617473.3856
424 -67953.7316 -439248.1597
425 -280137.2562 -67953.7316
426 -315144.1159 -280137.2562
427 -1188447.8307 -315144.1159
428 -186944.8306 -1188447.8307
429 -260412.3078 -186944.8306
430 -375551.4045 -260412.3078
> 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/7h2ae1291394141.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/8h2ae1291394141.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/9h2ae1291394141.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/10absh1291394141.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/11duqn1291394141.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/12zu7t1291394141.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/13vm521291394141.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/14y4381291394141.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/15252d1291394141.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/16ffz41291394141.tab")
+ }
>
> try(system("convert tmp/1lsvn1291394141.ps tmp/1lsvn1291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/2lsvn1291394141.ps tmp/2lsvn1291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/3w1u81291394141.ps tmp/3w1u81291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/4w1u81291394141.ps tmp/4w1u81291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/5w1u81291394141.ps tmp/5w1u81291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/66tbb1291394141.ps tmp/66tbb1291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/7h2ae1291394141.ps tmp/7h2ae1291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/8h2ae1291394141.ps tmp/8h2ae1291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/9h2ae1291394141.ps tmp/9h2ae1291394141.png",intern=TRUE))
character(0)
> try(system("convert tmp/10absh1291394141.ps tmp/10absh1291394141.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.199 2.131 27.349