R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(162556
+ ,213118
+ ,34643609608
+ ,6282154
+ ,29790
+ ,81767
+ ,2435838930
+ ,4321023
+ ,87550
+ ,153198
+ ,13412484900
+ ,4111912
+ ,84738
+ ,-26007
+ ,-2203781166
+ ,223193
+ ,54660
+ ,126942
+ ,6938649720
+ ,1491348
+ ,42634
+ ,157214
+ ,6702661676
+ ,1629616
+ ,40949
+ ,129352
+ ,5296835048
+ ,1398893
+ ,45187
+ ,234817
+ ,10610675779
+ ,1926517
+ ,37704
+ ,60448
+ ,2279131392
+ ,983660
+ ,16275
+ ,47818
+ ,778237950
+ ,1443586
+ ,25830
+ ,245546
+ ,6342453180
+ ,1073089
+ ,12679
+ ,48020
+ ,608845580
+ ,984885
+ ,18014
+ ,-1710
+ ,-30803940
+ ,1405225
+ ,43556
+ ,32648
+ ,1422016288
+ ,227132
+ ,24811
+ ,95350
+ ,2365728850
+ ,929118
+ ,6575
+ ,151352
+ ,995139400
+ ,1071292
+ ,7123
+ ,288170
+ ,2052634910
+ ,638830
+ ,21950
+ ,114337
+ ,2509697150
+ ,856956
+ ,37597
+ ,37884
+ ,1424324748
+ ,992426
+ ,17821
+ ,122844
+ ,2189202924
+ ,444477
+ ,12988
+ ,82340
+ ,1069431920
+ ,857217
+ ,22330
+ ,79801
+ ,1781956330
+ ,711969
+ ,13326
+ ,165548
+ ,2206092648
+ ,702380
+ ,16189
+ ,116384
+ ,1884140576
+ ,358589
+ ,7146
+ ,134028
+ ,957764088
+ ,297978
+ ,15824
+ ,63838
+ ,1010172512
+ ,585715
+ ,27664
+ ,74996
+ ,2074689344
+ ,657954
+ ,11920
+ ,31080
+ ,370473600
+ ,209458
+ ,8568
+ ,32168
+ ,275615424
+ ,786690
+ ,14416
+ ,49857
+ ,718738512
+ ,439798
+ ,3369
+ ,87161
+ ,293645409
+ ,688779
+ ,11819
+ ,106113
+ ,1254149547
+ ,574339
+ ,6984
+ ,80570
+ ,562700880
+ ,741409
+ ,4519
+ ,102129
+ ,461520951
+ ,597793
+ ,2220
+ ,301670
+ ,669707400
+ ,644190
+ ,18562
+ ,102313
+ ,1899133906
+ ,377934
+ ,10327
+ ,88577
+ ,914734679
+ ,640273
+ ,5336
+ ,112477
+ ,600177272
+ ,697458
+ ,2365
+ ,191778
+ ,453554970
+ ,550608
+ ,4069
+ ,79804
+ ,324722476
+ ,207393
+ ,8636
+ ,128294
+ ,1107946984
+ ,301607
+ ,13718
+ ,96448
+ ,1323073664
+ ,345783
+ ,4525
+ ,93811
+ ,424494775
+ ,501749
+ ,6869
+ ,117520
+ ,807244880
+ ,379983
+ ,4628
+ ,69159
+ ,320067852
+ ,387475
+ ,3689
+ ,101792
+ ,375510688
+ ,377305
+ ,4891
+ ,210568
+ ,1029888088
+ ,370837
+ ,7489
+ ,136996
+ ,1025963044
+ ,430866
+ ,4901
+ ,121920
+ ,597529920
+ ,469107
+ ,2284
+ ,76403
+ ,174504452
+ ,194493
+ ,3160
+ ,108094
+ ,341577040
+ ,530670
+ ,4150
+ ,134759
+ ,559249850
+ ,518365
+ ,7285
+ ,188873
+ ,1375939805
+ ,491303
+ ,1134
+ ,146216
+ ,165808944
+ ,527021
+ ,4658
+ ,156608
+ ,729480064
+ ,233773
+ ,2384
+ ,61348
+ ,146253632
+ ,405972
+ ,3748
+ ,50350
+ ,188711800
+ ,652925
+ ,5371
+ ,87720
+ ,471144120
+ ,446211
+ ,1285
+ ,99489
+ ,127843365
+ ,341340
+ ,9327
+ ,87419
+ ,815357013
+ ,387699
+ ,5565
+ ,94355
+ ,525085575
+ ,493408
+ ,1528
+ ,60326
+ ,92178128
+ ,146494
+ ,3122
+ ,94670
+ ,295559740
+ ,414462
+ ,7561
+ ,82425
+ ,623215425
+ ,364304
+ ,2675
+ ,59017
+ ,157870475
+ ,355178
+ ,13253
+ ,90829
+ ,1203756737
+ ,357760
+ ,880
+ ,80791
+ ,71096080
+ ,261216
+ ,2053
+ ,100423
+ ,206168419
+ ,397144
+ ,1424
+ ,131116
+ ,186709184
+ ,374943
+ ,4036
+ ,100269
+ ,404685684
+ ,424898
+ ,3045
+ ,27330
+ ,83219850
+ ,202055
+ ,5119
+ ,39039
+ ,199840641
+ ,378525
+ ,1431
+ ,106885
+ ,152952435
+ ,310768
+ ,554
+ ,79285
+ ,43923890
+ ,325738
+ ,1975
+ ,118881
+ ,234789975
+ ,394510
+ ,1765
+ ,77623
+ ,137004595
+ ,247060
+ ,1012
+ ,114768
+ ,116145216
+ ,368078
+ ,810
+ ,74015
+ ,59952150
+ ,236761
+ ,1280
+ ,69465
+ ,88915200
+ ,312378
+ ,666
+ ,117869
+ ,78500754
+ ,339836
+ ,1380
+ ,60982
+ ,84155160
+ ,347385
+ ,4677
+ ,90131
+ ,421542687
+ ,426280
+ ,876
+ ,138971
+ ,121738596
+ ,352850
+ ,814
+ ,39625
+ ,32254750
+ ,301881
+ ,514
+ ,102725
+ ,52800650
+ ,377516
+ ,5692
+ ,64239
+ ,365648388
+ ,357312
+ ,3642
+ ,90262
+ ,328734204
+ ,458343
+ ,540
+ ,103960
+ ,56138400
+ ,354228
+ ,2099
+ ,106611
+ ,223776489
+ ,308636
+ ,567
+ ,103345
+ ,58596615
+ ,386212
+ ,2001
+ ,95551
+ ,191197551
+ ,393343
+ ,2949
+ ,82903
+ ,244480947
+ ,378509
+ ,2253
+ ,63593
+ ,143275029
+ ,452469
+ ,6533
+ ,126910
+ ,829103030
+ ,364839
+ ,1889
+ ,37527
+ ,70888503
+ ,358649
+ ,3055
+ ,60247
+ ,184054585
+ ,376641
+ ,272
+ ,112995
+ ,30734640
+ ,429112
+ ,1414
+ ,70184
+ ,99240176
+ ,330546
+ ,2564
+ ,130140
+ ,333678960
+ ,403560
+ ,1383
+ ,73221
+ ,101264643
+ ,317892
+ ,1261
+ ,76114
+ ,95979754
+ ,307528
+ ,975
+ ,90534
+ ,88270650
+ ,235133
+ ,3366
+ ,108479
+ ,365140314
+ ,299243
+ ,576
+ ,113761
+ ,65526336
+ ,314073
+ ,1686
+ ,68696
+ ,115821456
+ ,368186
+ ,746
+ ,71561
+ ,53384506
+ ,269661
+ ,3192
+ ,59831
+ ,190980552
+ ,125390
+ ,2045
+ ,97890
+ ,200185050
+ ,510834
+ ,5702
+ ,101481
+ ,578644662
+ ,321896
+ ,1932
+ ,72954
+ ,140947128
+ ,249898
+ ,936
+ ,67939
+ ,63590904
+ ,408881
+ ,3437
+ ,48022
+ ,165051614
+ ,158492
+ ,5131
+ ,86111
+ ,441835541
+ ,292154
+ ,2397
+ ,74020
+ ,177425940
+ ,289513
+ ,1389
+ ,57530
+ ,79909170
+ ,378049
+ ,1503
+ ,56364
+ ,84715092
+ ,343466
+ ,402
+ ,84990
+ ,34165980
+ ,332743
+ ,2239
+ ,88590
+ ,198353010
+ ,442882
+ ,2234
+ ,77200
+ ,172464800
+ ,214215
+ ,837
+ ,61262
+ ,51276294
+ ,315688
+ ,10579
+ ,110309
+ ,1166958911
+ ,375195
+ ,875
+ ,67000
+ ,58625000
+ ,334280
+ ,1585
+ ,93099
+ ,147561915
+ ,355864
+ ,1659
+ ,107577
+ ,178470243
+ ,480382
+ ,2647
+ ,62920
+ ,166549240
+ ,353058
+ ,3294
+ ,75832
+ ,249790608
+ ,217193
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,94
+ ,60793
+ ,5714542
+ ,314533
+ ,422
+ ,57935
+ ,24448570
+ ,318056
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,34
+ ,60630
+ ,2061420
+ ,314353
+ ,1558
+ ,55637
+ ,86682446
+ ,369448
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,43
+ ,60887
+ ,2618141
+ ,312846
+ ,645
+ ,60720
+ ,39164400
+ ,312075
+ ,316
+ ,60505
+ ,19119580
+ ,315009
+ ,115
+ ,60945
+ ,7008675
+ ,318903
+ ,5
+ ,60720
+ ,303600
+ ,314887
+ ,897
+ ,60720
+ ,54465840
+ ,314913
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,389
+ ,58990
+ ,22947110
+ ,325506
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1002
+ ,56750
+ ,56863500
+ ,298568
+ ,36
+ ,60894
+ ,2192184
+ ,315834
+ ,460
+ ,63346
+ ,29139160
+ ,329784
+ ,309
+ ,56535
+ ,17469315
+ ,312878
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,9
+ ,60835
+ ,547515
+ ,314987
+ ,271
+ ,60720
+ ,16455120
+ ,325249
+ ,14
+ ,61016
+ ,854224
+ ,315877
+ ,520
+ ,58650
+ ,30498000
+ ,291650
+ ,1766
+ ,60438
+ ,106733508
+ ,305959
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,458
+ ,58625
+ ,26850250
+ ,297765
+ ,20
+ ,60938
+ ,1218760
+ ,315245
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,98
+ ,61490
+ ,6026020
+ ,315236
+ ,405
+ ,60845
+ ,24642225
+ ,336425
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,483
+ ,60830
+ ,29380890
+ ,306268
+ ,454
+ ,63261
+ ,28720494
+ ,302187
+ ,47
+ ,60720
+ ,2853840
+ ,314882
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,757
+ ,45689
+ ,34586573
+ ,382712
+ ,4655
+ ,60720
+ ,282651600
+ ,341570
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,36
+ ,61564
+ ,2216304
+ ,312412
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,203
+ ,61938
+ ,12573414
+ ,309596
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,126
+ ,60951
+ ,7679826
+ ,315547
+ ,400
+ ,60720
+ ,24288000
+ ,313267
+ ,71
+ ,60745
+ ,4312895
+ ,316176
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,972
+ ,71642
+ ,69636024
+ ,359335
+ ,531
+ ,71641
+ ,38041371
+ ,330068
+ ,2461
+ ,55792
+ ,137304112
+ ,314289
+ ,378
+ ,71873
+ ,27167994
+ ,297413
+ ,23
+ ,62555
+ ,1438765
+ ,314806
+ ,638
+ ,60370
+ ,38516060
+ ,333210
+ ,2300
+ ,64873
+ ,149207900
+ ,352108
+ ,149
+ ,62041
+ ,9244109
+ ,313332
+ ,226
+ ,65745
+ ,14858370
+ ,291787
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,275
+ ,59500
+ ,16362500
+ ,318745
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,141
+ ,61630
+ ,8689830
+ ,315366
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,28
+ ,60890
+ ,1704920
+ ,315688
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,4980
+ ,113521
+ ,565334580
+ ,409642
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,472
+ ,80045
+ ,37781240
+ ,269587
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,203
+ ,50804
+ ,10313212
+ ,300962
+ ,496
+ ,87390
+ ,43345440
+ ,325479
+ ,10
+ ,61656
+ ,616560
+ ,316155
+ ,63
+ ,65688
+ ,4138344
+ ,318574
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1136
+ ,48522
+ ,55120992
+ ,343613
+ ,265
+ ,60720
+ ,16090800
+ ,306948
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,267
+ ,57640
+ ,15389880
+ ,330059
+ ,474
+ ,61977
+ ,29377098
+ ,288985
+ ,534
+ ,62620
+ ,33439080
+ ,304485
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,15
+ ,60831
+ ,912465
+ ,315688
+ ,397
+ ,60646
+ ,24076462
+ ,317736
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1866
+ ,56225
+ ,104915850
+ ,322331
+ ,288
+ ,60510
+ ,17426880
+ ,296656
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,3
+ ,60698
+ ,182094
+ ,315354
+ ,468
+ ,60720
+ ,28416960
+ ,312161
+ ,20
+ ,60805
+ ,1216100
+ ,315576
+ ,278
+ ,61404
+ ,17070312
+ ,314922
+ ,61
+ ,60720
+ ,3703920
+ ,314551
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,192
+ ,65276
+ ,12532992
+ ,312339
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,317
+ ,63915
+ ,20261055
+ ,298700
+ ,738
+ ,60720
+ ,44811360
+ ,321376
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,368
+ ,61686
+ ,22700448
+ ,303230
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,2
+ ,60743
+ ,121486
+ ,315487
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,53
+ ,60349
+ ,3198497
+ ,315793
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,94
+ ,61360
+ ,5767840
+ ,312887
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,24
+ ,59818
+ ,1435632
+ ,315637
+ ,2332
+ ,72680
+ ,169489760
+ ,324385
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,131
+ ,61808
+ ,8096848
+ ,308989
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,206
+ ,53110
+ ,10940660
+ ,296702
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,167
+ ,64245
+ ,10728915
+ ,307322
+ ,622
+ ,73007
+ ,45410354
+ ,304376
+ ,2328
+ ,82732
+ ,192600096
+ ,253588
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,365
+ ,54820
+ ,20009300
+ ,309560
+ ,364
+ ,47705
+ ,17364620
+ ,298466
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,226
+ ,72835
+ ,16460710
+ ,343929
+ ,307
+ ,58856
+ ,18068792
+ ,331955
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,188
+ ,77655
+ ,14599140
+ ,381180
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,138
+ ,69817
+ ,9634746
+ ,331420
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,125
+ ,60798
+ ,7599750
+ ,310201
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,282
+ ,62452
+ ,17611464
+ ,320016
+ ,335
+ ,64175
+ ,21498625
+ ,320398
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1324
+ ,67440
+ ,89290560
+ ,291841
+ ,176
+ ,68136
+ ,11991936
+ ,310670
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,249
+ ,56726
+ ,14124774
+ ,313491
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,333
+ ,70811
+ ,23580063
+ ,331323
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,601
+ ,60720
+ ,36492720
+ ,319210
+ ,30
+ ,62045
+ ,1861350
+ ,318098
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,249
+ ,54323
+ ,13526427
+ ,292754
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,165
+ ,62841
+ ,10368765
+ ,325176
+ ,453
+ ,81125
+ ,36749625
+ ,365959
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,53
+ ,59506
+ ,3153818
+ ,302409
+ ,382
+ ,59365
+ ,22677430
+ ,340968
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,30
+ ,60798
+ ,1823940
+ ,313164
+ ,290
+ ,58790
+ ,17049100
+ ,301164
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,366
+ ,61808
+ ,22621728
+ ,344425
+ ,2
+ ,60735
+ ,121470
+ ,315394
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,209
+ ,64016
+ ,13379344
+ ,316647
+ ,384
+ ,54683
+ ,20998272
+ ,309836
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,365
+ ,87192
+ ,31825080
+ ,346611
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,49
+ ,64107
+ ,3141243
+ ,322031
+ ,3
+ ,60761
+ ,182283
+ ,315656
+ ,133
+ ,65990
+ ,8776670
+ ,339445
+ ,32
+ ,59988
+ ,1919616
+ ,314964
+ ,368
+ ,61167
+ ,22509456
+ ,297141
+ ,1
+ ,60719
+ ,60719
+ ,315372
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,22
+ ,60722
+ ,1335884
+ ,312502
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,96
+ ,60379
+ ,5796384
+ ,313729
+ ,1
+ ,60727
+ ,60727
+ ,315388
+ ,314
+ ,60720
+ ,19066080
+ ,315371
+ ,844
+ ,60925
+ ,51420700
+ ,296139
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,26
+ ,60896
+ ,1583296
+ ,313880
+ ,125
+ ,59734
+ ,7466750
+ ,317698
+ ,304
+ ,62969
+ ,19142576
+ ,295580
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,621
+ ,60720
+ ,37707120
+ ,308256
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,119
+ ,59118
+ ,7035042
+ ,303677
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1595
+ ,60720
+ ,96848400
+ ,319369
+ ,312
+ ,58598
+ ,18282576
+ ,318690
+ ,60
+ ,61124
+ ,3667440
+ ,314049
+ ,587
+ ,59595
+ ,34982265
+ ,325699
+ ,135
+ ,62065
+ ,8378775
+ ,314210
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,514
+ ,78780
+ ,40492920
+ ,322378
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,60722
+ ,60722
+ ,315398
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1763
+ ,61600
+ ,108600800
+ ,308336
+ ,180
+ ,59635
+ ,10734300
+ ,316386
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,218
+ ,60720
+ ,13236960
+ ,315553
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,448
+ ,59781
+ ,26781888
+ ,323361
+ ,227
+ ,76644
+ ,17398188
+ ,336639
+ ,174
+ ,64820
+ ,11278680
+ ,307424
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,121
+ ,56178
+ ,6797538
+ ,295370
+ ,607
+ ,60436
+ ,36684652
+ ,322340
+ ,2212
+ ,60720
+ ,134312640
+ ,319864
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,530
+ ,73433
+ ,38919490
+ ,317291
+ ,571
+ ,41477
+ ,23683367
+ ,280398
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,78
+ ,62700
+ ,4890600
+ ,317330
+ ,2489
+ ,67804
+ ,168764156
+ ,238125
+ ,131
+ ,59661
+ ,7815591
+ ,327071
+ ,923
+ ,58620
+ ,54106260
+ ,309038
+ ,72
+ ,60398
+ ,4348656
+ ,314210
+ ,572
+ ,58580
+ ,33507760
+ ,307930
+ ,397
+ ,62710
+ ,24895870
+ ,322327
+ ,450
+ ,59325
+ ,26696250
+ ,292136
+ ,622
+ ,60950
+ ,37910900
+ ,263276
+ ,694
+ ,68060
+ ,47233640
+ ,367655
+ ,3425
+ ,83620
+ ,286398500
+ ,283910
+ ,562
+ ,58456
+ ,32852272
+ ,283587
+ ,4917
+ ,52811
+ ,259671687
+ ,243650
+ ,1442
+ ,121173
+ ,174731466
+ ,438493
+ ,529
+ ,63870
+ ,33787230
+ ,296261
+ ,2126
+ ,21001
+ ,44648126
+ ,230621
+ ,1061
+ ,70415
+ ,74710315
+ ,304252
+ ,776
+ ,64230
+ ,49842480
+ ,333505
+ ,611
+ ,59190
+ ,36165090
+ ,296919
+ ,1526
+ ,69351
+ ,105829626
+ ,278990
+ ,592
+ ,64270
+ ,38047840
+ ,276898
+ ,1182
+ ,70694
+ ,83560308
+ ,327007
+ ,621
+ ,68005
+ ,42231105
+ ,317046
+ ,989
+ ,58930
+ ,58281770
+ ,304555
+ ,438
+ ,58320
+ ,25544160
+ ,298096
+ ,726
+ ,69980
+ ,50805480
+ ,231861
+ ,1303
+ ,69863
+ ,91031489
+ ,309422
+ ,7419
+ ,63255
+ ,469288845
+ ,286963
+ ,1164
+ ,57320
+ ,66720480
+ ,269753
+ ,3310
+ ,75230
+ ,249011300
+ ,448243
+ ,1920
+ ,79420
+ ,152486400
+ ,165404
+ ,965
+ ,73490
+ ,70917850
+ ,204325
+ ,3256
+ ,35250
+ ,114774000
+ ,407159
+ ,1135
+ ,62285
+ ,70693475
+ ,290476
+ ,1270
+ ,69206
+ ,87891620
+ ,275311
+ ,661
+ ,65920
+ ,43573120
+ ,246541
+ ,1013
+ ,69770
+ ,70677010
+ ,253468
+ ,2844
+ ,72683
+ ,206710452
+ ,240897
+ ,11528
+ ,-14545
+ ,-167674760
+ ,-83265
+ ,6526
+ ,55830
+ ,364346580
+ ,-42143
+ ,2264
+ ,55174
+ ,124913936
+ ,272713
+ ,5109
+ ,67038
+ ,342497142
+ ,215362
+ ,3999
+ ,51252
+ ,204956748
+ ,42754
+ ,35624
+ ,157278
+ ,5602871472
+ ,306275
+ ,9252
+ ,79510
+ ,735626520
+ ,253537
+ ,15236
+ ,77440
+ ,1179875840
+ ,372631
+ ,18073
+ ,27284
+ ,493103732
+ ,-7170)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Costs'
+ ,'Dividends'
+ ,'C_D'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Costs','Dividends','C_D','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 = '4'
> #'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 Costs Dividends C_D
1 6282154 162556 213118 34643609608
2 4321023 29790 81767 2435838930
3 4111912 87550 153198 13412484900
4 223193 84738 -26007 -2203781166
5 1491348 54660 126942 6938649720
6 1629616 42634 157214 6702661676
7 1398893 40949 129352 5296835048
8 1926517 45187 234817 10610675779
9 983660 37704 60448 2279131392
10 1443586 16275 47818 778237950
11 1073089 25830 245546 6342453180
12 984885 12679 48020 608845580
13 1405225 18014 -1710 -30803940
14 227132 43556 32648 1422016288
15 929118 24811 95350 2365728850
16 1071292 6575 151352 995139400
17 638830 7123 288170 2052634910
18 856956 21950 114337 2509697150
19 992426 37597 37884 1424324748
20 444477 17821 122844 2189202924
21 857217 12988 82340 1069431920
22 711969 22330 79801 1781956330
23 702380 13326 165548 2206092648
24 358589 16189 116384 1884140576
25 297978 7146 134028 957764088
26 585715 15824 63838 1010172512
27 657954 27664 74996 2074689344
28 209458 11920 31080 370473600
29 786690 8568 32168 275615424
30 439798 14416 49857 718738512
31 688779 3369 87161 293645409
32 574339 11819 106113 1254149547
33 741409 6984 80570 562700880
34 597793 4519 102129 461520951
35 644190 2220 301670 669707400
36 377934 18562 102313 1899133906
37 640273 10327 88577 914734679
38 697458 5336 112477 600177272
39 550608 2365 191778 453554970
40 207393 4069 79804 324722476
41 301607 8636 128294 1107946984
42 345783 13718 96448 1323073664
43 501749 4525 93811 424494775
44 379983 6869 117520 807244880
45 387475 4628 69159 320067852
46 377305 3689 101792 375510688
47 370837 4891 210568 1029888088
48 430866 7489 136996 1025963044
49 469107 4901 121920 597529920
50 194493 2284 76403 174504452
51 530670 3160 108094 341577040
52 518365 4150 134759 559249850
53 491303 7285 188873 1375939805
54 527021 1134 146216 165808944
55 233773 4658 156608 729480064
56 405972 2384 61348 146253632
57 652925 3748 50350 188711800
58 446211 5371 87720 471144120
59 341340 1285 99489 127843365
60 387699 9327 87419 815357013
61 493408 5565 94355 525085575
62 146494 1528 60326 92178128
63 414462 3122 94670 295559740
64 364304 7561 82425 623215425
65 355178 2675 59017 157870475
66 357760 13253 90829 1203756737
67 261216 880 80791 71096080
68 397144 2053 100423 206168419
69 374943 1424 131116 186709184
70 424898 4036 100269 404685684
71 202055 3045 27330 83219850
72 378525 5119 39039 199840641
73 310768 1431 106885 152952435
74 325738 554 79285 43923890
75 394510 1975 118881 234789975
76 247060 1765 77623 137004595
77 368078 1012 114768 116145216
78 236761 810 74015 59952150
79 312378 1280 69465 88915200
80 339836 666 117869 78500754
81 347385 1380 60982 84155160
82 426280 4677 90131 421542687
83 352850 876 138971 121738596
84 301881 814 39625 32254750
85 377516 514 102725 52800650
86 357312 5692 64239 365648388
87 458343 3642 90262 328734204
88 354228 540 103960 56138400
89 308636 2099 106611 223776489
90 386212 567 103345 58596615
91 393343 2001 95551 191197551
92 378509 2949 82903 244480947
93 452469 2253 63593 143275029
94 364839 6533 126910 829103030
95 358649 1889 37527 70888503
96 376641 3055 60247 184054585
97 429112 272 112995 30734640
98 330546 1414 70184 99240176
99 403560 2564 130140 333678960
100 317892 1383 73221 101264643
101 307528 1261 76114 95979754
102 235133 975 90534 88270650
103 299243 3366 108479 365140314
104 314073 576 113761 65526336
105 368186 1686 68696 115821456
106 269661 746 71561 53384506
107 125390 3192 59831 190980552
108 510834 2045 97890 200185050
109 321896 5702 101481 578644662
110 249898 1932 72954 140947128
111 408881 936 67939 63590904
112 158492 3437 48022 165051614
113 292154 5131 86111 441835541
114 289513 2397 74020 177425940
115 378049 1389 57530 79909170
116 343466 1503 56364 84715092
117 332743 402 84990 34165980
118 442882 2239 88590 198353010
119 214215 2234 77200 172464800
120 315688 837 61262 51276294
121 375195 10579 110309 1166958911
122 334280 875 67000 58625000
123 355864 1585 93099 147561915
124 480382 1659 107577 178470243
125 353058 2647 62920 166549240
126 217193 3294 75832 249790608
127 315380 0 60720 0
128 314533 94 60793 5714542
129 318056 422 57935 24448570
130 315380 0 60720 0
131 314353 34 60630 2061420
132 369448 1558 55637 86682446
133 315380 0 60720 0
134 312846 43 60887 2618141
135 312075 645 60720 39164400
136 315009 316 60505 19119580
137 318903 115 60945 7008675
138 314887 5 60720 303600
139 314913 897 60720 54465840
140 315380 0 60720 0
141 325506 389 58990 22947110
142 315380 0 60720 0
143 298568 1002 56750 56863500
144 315834 36 60894 2192184
145 329784 460 63346 29139160
146 312878 309 56535 17469315
147 315380 0 60720 0
148 314987 9 60835 547515
149 325249 271 60720 16455120
150 315877 14 61016 854224
151 291650 520 58650 30498000
152 305959 1766 60438 106733508
153 315380 0 60720 0
154 297765 458 58625 26850250
155 315245 20 60938 1218760
156 315380 0 60720 0
157 315380 0 60720 0
158 315236 98 61490 6026020
159 336425 405 60845 24642225
160 315380 0 60720 0
161 315380 0 60720 0
162 315380 0 60720 0
163 315380 0 60720 0
164 306268 483 60830 29380890
165 302187 454 63261 28720494
166 314882 47 60720 2853840
167 315380 0 60720 0
168 382712 757 45689 34586573
169 341570 4655 60720 282651600
170 315380 0 60720 0
171 315380 0 60720 0
172 312412 36 61564 2216304
173 315380 0 60720 0
174 309596 203 61938 12573414
175 315380 0 60720 0
176 315547 126 60951 7679826
177 313267 400 60720 24288000
178 316176 71 60745 4312895
179 315380 0 60720 0
180 315380 0 60720 0
181 359335 972 71642 69636024
182 330068 531 71641 38041371
183 314289 2461 55792 137304112
184 297413 378 71873 27167994
185 314806 23 62555 1438765
186 333210 638 60370 38516060
187 352108 2300 64873 149207900
188 313332 149 62041 9244109
189 291787 226 65745 14858370
190 315380 0 60720 0
191 318745 275 59500 16362500
192 315380 0 60720 0
193 315366 141 61630 8689830
194 315380 0 60720 0
195 315688 28 60890 1704920
196 315380 0 60720 0
197 409642 4980 113521 565334580
198 315380 0 60720 0
199 315380 0 60720 0
200 269587 472 80045 37781240
201 315380 0 60720 0
202 315380 0 60720 0
203 315380 0 60720 0
204 300962 203 50804 10313212
205 325479 496 87390 43345440
206 316155 10 61656 616560
207 318574 63 65688 4138344
208 315380 0 60720 0
209 343613 1136 48522 55120992
210 306948 265 60720 16090800
211 315380 0 60720 0
212 315380 0 60720 0
213 330059 267 57640 15389880
214 288985 474 61977 29377098
215 304485 534 62620 33439080
216 315380 0 60720 0
217 315688 15 60831 912465
218 317736 397 60646 24076462
219 315380 0 60720 0
220 322331 1866 56225 104915850
221 296656 288 60510 17426880
222 315380 0 60720 0
223 315354 3 60698 182094
224 312161 468 60720 28416960
225 315576 20 60805 1216100
226 314922 278 61404 17070312
227 314551 61 60720 3703920
228 315380 0 60720 0
229 312339 192 65276 12532992
230 315380 0 60720 0
231 298700 317 63915 20261055
232 321376 738 60720 44811360
233 315380 0 60720 0
234 303230 368 61686 22700448
235 315380 0 60720 0
236 315487 2 60743 121486
237 315380 0 60720 0
238 315793 53 60349 3198497
239 315380 0 60720 0
240 315380 0 60720 0
241 315380 0 60720 0
242 312887 94 61360 5767840
243 315380 0 60720 0
244 315637 24 59818 1435632
245 324385 2332 72680 169489760
246 315380 0 60720 0
247 315380 0 60720 0
248 308989 131 61808 8096848
249 315380 0 60720 0
250 315380 0 60720 0
251 296702 206 53110 10940660
252 315380 0 60720 0
253 307322 167 64245 10728915
254 304376 622 73007 45410354
255 253588 2328 82732 192600096
256 315380 0 60720 0
257 309560 365 54820 20009300
258 298466 364 47705 17364620
259 315380 0 60720 0
260 315380 0 60720 0
261 315380 0 60720 0
262 315380 0 60720 0
263 343929 226 72835 16460710
264 331955 307 58856 18068792
265 315380 0 60720 0
266 315380 0 60720 0
267 315380 0 60720 0
268 381180 188 77655 14599140
269 315380 0 60720 0
270 331420 138 69817 9634746
271 315380 0 60720 0
272 315380 0 60720 0
273 315380 0 60720 0
274 310201 125 60798 7599750
275 315380 0 60720 0
276 320016 282 62452 17611464
277 320398 335 64175 21498625
278 315380 0 60720 0
279 291841 1324 67440 89290560
280 310670 176 68136 11991936
281 315380 0 60720 0
282 315380 0 60720 0
283 313491 249 56726 14124774
284 315380 0 60720 0
285 331323 333 70811 23580063
286 315380 0 60720 0
287 319210 601 60720 36492720
288 318098 30 62045 1861350
289 315380 0 60720 0
290 292754 249 54323 13526427
291 315380 0 60720 0
292 325176 165 62841 10368765
293 365959 453 81125 36749625
294 315380 0 60720 0
295 302409 53 59506 3153818
296 340968 382 59365 22677430
297 315380 0 60720 0
298 315380 0 60720 0
299 315380 0 60720 0
300 315380 0 60720 0
301 313164 30 60798 1823940
302 301164 290 58790 17049100
303 315380 0 60720 0
304 315380 0 60720 0
305 344425 366 61808 22621728
306 315394 2 60735 121470
307 315380 0 60720 0
308 316647 209 64016 13379344
309 309836 384 54683 20998272
310 315380 0 60720 0
311 315380 0 60720 0
312 346611 365 87192 31825080
313 315380 0 60720 0
314 322031 49 64107 3141243
315 315656 3 60761 182283
316 339445 133 65990 8776670
317 314964 32 59988 1919616
318 297141 368 61167 22509456
319 315372 1 60719 60719
320 315380 0 60720 0
321 315380 0 60720 0
322 315380 0 60720 0
323 315380 0 60720 0
324 315380 0 60720 0
325 315380 0 60720 0
326 312502 22 60722 1335884
327 315380 0 60720 0
328 315380 0 60720 0
329 315380 0 60720 0
330 315380 0 60720 0
331 315380 0 60720 0
332 315380 0 60720 0
333 315380 0 60720 0
334 313729 96 60379 5796384
335 315388 1 60727 60727
336 315371 314 60720 19066080
337 296139 844 60925 51420700
338 315380 0 60720 0
339 313880 26 60896 1583296
340 317698 125 59734 7466750
341 295580 304 62969 19142576
342 315380 0 60720 0
343 315380 0 60720 0
344 315380 0 60720 0
345 308256 621 60720 37707120
346 315380 0 60720 0
347 303677 119 59118 7035042
348 315380 0 60720 0
349 315380 0 60720 0
350 319369 1595 60720 96848400
351 318690 312 58598 18282576
352 314049 60 61124 3667440
353 325699 587 59595 34982265
354 314210 135 62065 8378775
355 315380 0 60720 0
356 315380 0 60720 0
357 322378 514 78780 40492920
358 315380 0 60720 0
359 315380 0 60720 0
360 315380 0 60720 0
361 315398 1 60722 60722
362 315380 0 60720 0
363 315380 0 60720 0
364 308336 1763 61600 108600800
365 316386 180 59635 10734300
366 315380 0 60720 0
367 315380 0 60720 0
368 315380 0 60720 0
369 315380 0 60720 0
370 315553 218 60720 13236960
371 315380 0 60720 0
372 323361 448 59781 26781888
373 336639 227 76644 17398188
374 307424 174 64820 11278680
375 315380 0 60720 0
376 315380 0 60720 0
377 295370 121 56178 6797538
378 322340 607 60436 36684652
379 319864 2212 60720 134312640
380 315380 0 60720 0
381 315380 0 60720 0
382 317291 530 73433 38919490
383 280398 571 41477 23683367
384 315380 0 60720 0
385 317330 78 62700 4890600
386 238125 2489 67804 168764156
387 327071 131 59661 7815591
388 309038 923 58620 54106260
389 314210 72 60398 4348656
390 307930 572 58580 33507760
391 322327 397 62710 24895870
392 292136 450 59325 26696250
393 263276 622 60950 37910900
394 367655 694 68060 47233640
395 283910 3425 83620 286398500
396 283587 562 58456 32852272
397 243650 4917 52811 259671687
398 438493 1442 121173 174731466
399 296261 529 63870 33787230
400 230621 2126 21001 44648126
401 304252 1061 70415 74710315
402 333505 776 64230 49842480
403 296919 611 59190 36165090
404 278990 1526 69351 105829626
405 276898 592 64270 38047840
406 327007 1182 70694 83560308
407 317046 621 68005 42231105
408 304555 989 58930 58281770
409 298096 438 58320 25544160
410 231861 726 69980 50805480
411 309422 1303 69863 91031489
412 286963 7419 63255 469288845
413 269753 1164 57320 66720480
414 448243 3310 75230 249011300
415 165404 1920 79420 152486400
416 204325 965 73490 70917850
417 407159 3256 35250 114774000
418 290476 1135 62285 70693475
419 275311 1270 69206 87891620
420 246541 661 65920 43573120
421 253468 1013 69770 70677010
422 240897 2844 72683 206710452
423 -83265 11528 -14545 -167674760
424 -42143 6526 55830 364346580
425 272713 2264 55174 124913936
426 215362 5109 67038 342497142
427 42754 3999 51252 204956748
428 306275 35624 157278 5602871472
429 253537 9252 79510 735626520
430 372631 15236 77440 1179875840
431 -7170 18073 27284 493103732
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Dividends C_D
3.109e+05 8.950e+00 -5.554e-02 1.382e-04
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-1089227 -15812 7886 7886 3411356
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.109e+05 2.933e+04 10.598 < 2e-16 ***
Costs 8.950e+00 1.791e+00 4.997 8.49e-07 ***
Dividends -5.554e-02 3.826e-01 -0.145 0.885
C_D 1.382e-04 1.138e-05 12.144 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 220700 on 427 degrees of freedom
Multiple R-squared: 0.739, Adjusted R-squared: 0.7372
F-statistic: 403.1 on 3 and 427 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 7.558739e-55 3.779370e-55
[2,] 1.0000000 1.634946e-59 8.174728e-60
[3,] 1.0000000 3.024277e-79 1.512138e-79
[4,] 1.0000000 4.762220e-93 2.381110e-93
[5,] 1.0000000 3.753788e-94 1.876894e-94
[6,] 1.0000000 3.630726e-102 1.815363e-102
[7,] 1.0000000 1.268983e-117 6.344917e-118
[8,] 1.0000000 2.504837e-129 1.252418e-129
[9,] 1.0000000 2.276825e-132 1.138412e-132
[10,] 1.0000000 8.948970e-140 4.474485e-140
[11,] 1.0000000 6.240380e-139 3.120190e-139
[12,] 1.0000000 1.918912e-142 9.594560e-143
[13,] 1.0000000 8.871789e-146 4.435895e-146
[14,] 1.0000000 6.379859e-148 3.189929e-148
[15,] 1.0000000 1.526249e-153 7.631243e-154
[16,] 1.0000000 1.205595e-156 6.027975e-157
[17,] 1.0000000 1.164456e-157 5.822281e-158
[18,] 1.0000000 2.697554e-159 1.348777e-159
[19,] 1.0000000 4.650409e-160 2.325205e-160
[20,] 1.0000000 8.462954e-162 4.231477e-162
[21,] 1.0000000 5.562963e-165 2.781481e-165
[22,] 1.0000000 8.265555e-167 4.132777e-167
[23,] 1.0000000 1.674151e-176 8.370754e-177
[24,] 1.0000000 8.793400e-178 4.396700e-178
[25,] 1.0000000 5.965908e-182 2.982954e-182
[26,] 1.0000000 2.487706e-183 1.243853e-183
[27,] 1.0000000 2.447110e-191 1.223555e-191
[28,] 1.0000000 2.205342e-193 1.102671e-193
[29,] 1.0000000 3.531335e-193 1.765667e-193
[30,] 1.0000000 8.341821e-194 4.170910e-194
[31,] 1.0000000 2.350140e-199 1.175070e-199
[32,] 1.0000000 2.427829e-205 1.213914e-205
[33,] 1.0000000 1.700204e-204 8.501019e-205
[34,] 1.0000000 3.727426e-205 1.863713e-205
[35,] 1.0000000 3.345291e-205 1.672646e-205
[36,] 1.0000000 3.585482e-205 1.792741e-205
[37,] 1.0000000 5.042446e-206 2.521223e-206
[38,] 1.0000000 3.112342e-205 1.556171e-205
[39,] 1.0000000 5.741230e-205 2.870615e-205
[40,] 1.0000000 4.480651e-204 2.240325e-204
[41,] 1.0000000 6.860258e-205 3.430129e-205
[42,] 1.0000000 5.184691e-204 2.592346e-204
[43,] 1.0000000 1.713328e-203 8.566640e-204
[44,] 1.0000000 3.876786e-204 1.938393e-204
[45,] 1.0000000 4.609710e-205 2.304855e-205
[46,] 1.0000000 5.106033e-205 2.553016e-205
[47,] 1.0000000 5.489092e-204 2.744546e-204
[48,] 1.0000000 7.389679e-204 3.694839e-204
[49,] 1.0000000 3.111298e-206 1.555649e-206
[50,] 1.0000000 3.465386e-206 1.732693e-206
[51,] 1.0000000 2.930881e-218 1.465440e-218
[52,] 1.0000000 7.686953e-219 3.843476e-219
[53,] 1.0000000 7.280669e-218 3.640335e-218
[54,] 1.0000000 4.973890e-218 2.486945e-218
[55,] 1.0000000 7.673854e-220 3.836927e-220
[56,] 1.0000000 8.951492e-222 4.475746e-222
[57,] 1.0000000 2.686801e-221 1.343401e-221
[58,] 1.0000000 4.330265e-221 2.165133e-221
[59,] 1.0000000 1.342294e-220 6.711471e-221
[60,] 1.0000000 3.855560e-221 1.927780e-221
[61,] 1.0000000 8.985430e-221 4.492715e-221
[62,] 1.0000000 6.594840e-220 3.297420e-220
[63,] 1.0000000 6.627296e-219 3.313648e-219
[64,] 1.0000000 1.161522e-218 5.807608e-219
[65,] 1.0000000 4.618133e-218 2.309067e-218
[66,] 1.0000000 8.512287e-220 4.256143e-220
[67,] 1.0000000 4.009951e-219 2.004976e-219
[68,] 1.0000000 4.352955e-218 2.176477e-218
[69,] 1.0000000 4.935109e-217 2.467554e-217
[70,] 1.0000000 9.628619e-217 4.814310e-217
[71,] 1.0000000 1.054554e-215 5.272768e-216
[72,] 1.0000000 1.095268e-215 5.476338e-216
[73,] 1.0000000 1.138856e-214 5.694279e-215
[74,] 1.0000000 6.205893e-214 3.102946e-214
[75,] 1.0000000 3.822271e-213 1.911136e-213
[76,] 1.0000000 1.772649e-213 8.863243e-214
[77,] 1.0000000 4.945300e-213 2.472650e-213
[78,] 1.0000000 3.941939e-212 1.970970e-212
[79,] 1.0000000 4.419260e-211 2.209630e-211
[80,] 1.0000000 2.809602e-211 1.404801e-211
[81,] 1.0000000 2.430974e-212 1.215487e-212
[82,] 1.0000000 2.441341e-211 1.220671e-211
[83,] 1.0000000 9.372766e-211 4.686383e-211
[84,] 1.0000000 1.004517e-209 5.022585e-210
[85,] 1.0000000 6.529960e-209 3.264980e-209
[86,] 1.0000000 2.577418e-208 1.288709e-208
[87,] 1.0000000 1.853775e-210 9.268873e-211
[88,] 1.0000000 1.660109e-209 8.300545e-210
[89,] 1.0000000 6.771394e-210 3.385697e-210
[90,] 1.0000000 3.565175e-210 1.782588e-210
[91,] 1.0000000 2.394584e-209 1.197292e-209
[92,] 1.0000000 2.287818e-208 1.143909e-208
[93,] 1.0000000 2.505376e-207 1.252688e-207
[94,] 1.0000000 2.574386e-206 1.287193e-206
[95,] 1.0000000 2.410556e-205 1.205278e-205
[96,] 1.0000000 3.631700e-206 1.815850e-206
[97,] 1.0000000 8.829328e-206 4.414664e-206
[98,] 1.0000000 8.089014e-206 4.044507e-206
[99,] 1.0000000 2.968634e-205 1.484317e-205
[100,] 1.0000000 1.052902e-204 5.264510e-205
[101,] 1.0000000 3.714804e-207 1.857402e-207
[102,] 1.0000000 1.745658e-209 8.728292e-210
[103,] 1.0000000 1.520965e-208 7.604824e-209
[104,] 1.0000000 3.336526e-208 1.668263e-208
[105,] 1.0000000 1.800098e-208 9.000488e-209
[106,] 1.0000000 5.547082e-209 2.773541e-209
[107,] 1.0000000 4.611185e-208 2.305593e-208
[108,] 1.0000000 3.936575e-207 1.968288e-207
[109,] 1.0000000 3.205608e-207 1.602804e-207
[110,] 1.0000000 1.209478e-206 6.047392e-207
[111,] 1.0000000 1.180234e-205 5.901172e-206
[112,] 1.0000000 2.266630e-206 1.133315e-206
[113,] 1.0000000 3.342081e-207 1.671040e-207
[114,] 1.0000000 3.534126e-206 1.767063e-206
[115,] 1.0000000 6.052355e-206 3.026177e-206
[116,] 1.0000000 5.754796e-205 2.877398e-205
[117,] 1.0000000 6.262828e-204 3.131414e-204
[118,] 1.0000000 4.026910e-205 2.013455e-205
[119,] 1.0000000 7.201418e-205 3.600709e-205
[120,] 1.0000000 4.234174e-205 2.117087e-205
[121,] 1.0000000 4.703443e-204 2.351722e-204
[122,] 1.0000000 5.189406e-203 2.594703e-203
[123,] 1.0000000 5.366136e-202 2.683068e-202
[124,] 1.0000000 5.865806e-201 2.932903e-201
[125,] 1.0000000 6.372923e-200 3.186462e-200
[126,] 1.0000000 4.897763e-200 2.448882e-200
[127,] 1.0000000 5.314889e-199 2.657445e-199
[128,] 1.0000000 5.724130e-198 2.862065e-198
[129,] 1.0000000 6.028740e-197 3.014370e-197
[130,] 1.0000000 6.369703e-196 3.184852e-196
[131,] 1.0000000 6.694837e-195 3.347419e-195
[132,] 1.0000000 7.074921e-194 3.537461e-194
[133,] 1.0000000 7.010041e-193 3.505021e-193
[134,] 1.0000000 7.332266e-192 3.666133e-192
[135,] 1.0000000 6.713162e-191 3.356581e-191
[136,] 1.0000000 6.950890e-190 3.475445e-190
[137,] 1.0000000 6.992126e-189 3.496063e-189
[138,] 1.0000000 7.161138e-188 3.580569e-188
[139,] 1.0000000 6.508263e-187 3.254131e-187
[140,] 1.0000000 6.421796e-186 3.210898e-186
[141,] 1.0000000 6.482457e-185 3.241229e-185
[142,] 1.0000000 6.510466e-184 3.255233e-184
[143,] 1.0000000 6.030276e-183 3.015138e-183
[144,] 1.0000000 5.994008e-182 2.997004e-182
[145,] 1.0000000 5.433609e-181 2.716805e-181
[146,] 1.0000000 4.862022e-180 2.431011e-180
[147,] 1.0000000 4.764204e-179 2.382102e-179
[148,] 1.0000000 4.477571e-178 2.238786e-178
[149,] 1.0000000 4.344059e-177 2.172030e-177
[150,] 1.0000000 4.193122e-176 2.096561e-176
[151,] 1.0000000 4.026944e-175 2.013472e-175
[152,] 1.0000000 3.845954e-174 1.922977e-174
[153,] 1.0000000 2.803278e-173 1.401639e-173
[154,] 1.0000000 2.654106e-172 1.327053e-172
[155,] 1.0000000 2.500278e-171 1.250139e-171
[156,] 1.0000000 2.343587e-170 1.171793e-170
[157,] 1.0000000 2.185755e-169 1.092878e-169
[158,] 1.0000000 2.012778e-168 1.006389e-168
[159,] 1.0000000 1.779679e-167 8.898394e-168
[160,] 1.0000000 1.635131e-166 8.175654e-167
[161,] 1.0000000 1.495471e-165 7.477357e-166
[162,] 1.0000000 3.646600e-166 1.823300e-166
[163,] 1.0000000 2.174807e-166 1.087404e-166
[164,] 1.0000000 2.008492e-165 1.004246e-165
[165,] 1.0000000 1.845789e-164 9.228946e-165
[166,] 1.0000000 1.686453e-163 8.432263e-164
[167,] 1.0000000 1.534749e-162 7.673743e-163
[168,] 1.0000000 1.383008e-161 6.915038e-162
[169,] 1.0000000 1.246432e-160 6.232162e-161
[170,] 1.0000000 1.114494e-159 5.572469e-160
[171,] 1.0000000 9.855525e-159 4.927763e-159
[172,] 1.0000000 8.725100e-158 4.362550e-158
[173,] 1.0000000 7.712300e-157 3.856150e-157
[174,] 1.0000000 6.784104e-156 3.392052e-156
[175,] 1.0000000 3.198533e-155 1.599267e-155
[176,] 1.0000000 2.743457e-154 1.371728e-154
[177,] 1.0000000 1.327097e-153 6.635487e-154
[178,] 1.0000000 8.779941e-153 4.389970e-153
[179,] 1.0000000 7.594495e-152 3.797248e-152
[180,] 1.0000000 4.910273e-151 2.455136e-151
[181,] 1.0000000 1.117115e-150 5.585576e-151
[182,] 1.0000000 9.585494e-150 4.792747e-150
[183,] 1.0000000 6.433622e-149 3.216811e-149
[184,] 1.0000000 5.468699e-148 2.734350e-148
[185,] 1.0000000 4.443757e-147 2.221879e-147
[186,] 1.0000000 3.741778e-146 1.870889e-146
[187,] 1.0000000 3.131839e-145 1.565920e-145
[188,] 1.0000000 2.612093e-144 1.306047e-144
[189,] 1.0000000 2.166570e-143 1.083285e-143
[190,] 1.0000000 1.789901e-142 8.949503e-143
[191,] 1.0000000 3.679373e-142 1.839686e-142
[192,] 1.0000000 3.030549e-141 1.515274e-141
[193,] 1.0000000 2.484324e-140 1.242162e-140
[194,] 1.0000000 5.030474e-140 2.515237e-140
[195,] 1.0000000 4.111890e-139 2.055945e-139
[196,] 1.0000000 3.345163e-138 1.672582e-138
[197,] 1.0000000 2.708555e-137 1.354278e-137
[198,] 1.0000000 2.166090e-136 1.083045e-136
[199,] 1.0000000 1.555402e-135 7.777011e-136
[200,] 1.0000000 1.243118e-134 6.215588e-135
[201,] 1.0000000 9.900281e-134 4.950140e-134
[202,] 1.0000000 7.833742e-133 3.916871e-133
[203,] 1.0000000 1.617030e-132 8.085151e-133
[204,] 1.0000000 1.271440e-131 6.357201e-132
[205,] 1.0000000 9.976882e-131 4.988441e-131
[206,] 1.0000000 7.791925e-130 3.895963e-130
[207,] 1.0000000 4.926640e-129 2.463320e-129
[208,] 1.0000000 3.331382e-128 1.665691e-128
[209,] 1.0000000 2.547574e-127 1.273787e-127
[210,] 1.0000000 1.955889e-126 9.779443e-127
[211,] 1.0000000 1.493649e-125 7.468244e-126
[212,] 1.0000000 1.100506e-124 5.502532e-125
[213,] 1.0000000 8.331147e-124 4.165574e-124
[214,] 1.0000000 3.669983e-123 1.834992e-123
[215,] 1.0000000 2.627360e-122 1.313680e-122
[216,] 1.0000000 1.966481e-121 9.832403e-122
[217,] 1.0000000 1.464775e-120 7.323875e-121
[218,] 1.0000000 1.076157e-119 5.380784e-120
[219,] 1.0000000 7.936421e-119 3.968210e-119
[220,] 1.0000000 5.795578e-118 2.897789e-118
[221,] 1.0000000 4.236248e-117 2.118124e-117
[222,] 1.0000000 3.082147e-116 1.541074e-116
[223,] 1.0000000 2.223500e-115 1.111750e-115
[224,] 1.0000000 1.602636e-114 8.013182e-115
[225,] 1.0000000 1.076545e-113 5.382726e-114
[226,] 1.0000000 6.957862e-113 3.478931e-113
[227,] 1.0000000 4.948663e-112 2.474332e-112
[228,] 1.0000000 3.463428e-111 1.731714e-111
[229,] 1.0000000 2.440390e-110 1.220195e-110
[230,] 1.0000000 1.711081e-109 8.555403e-110
[231,] 1.0000000 1.194230e-108 5.971151e-109
[232,] 1.0000000 8.264723e-108 4.132361e-108
[233,] 1.0000000 5.713525e-107 2.856763e-107
[234,] 1.0000000 3.931021e-106 1.965510e-106
[235,] 1.0000000 2.691718e-105 1.345859e-105
[236,] 1.0000000 1.838261e-104 9.191307e-105
[237,] 1.0000000 1.246743e-103 6.233717e-104
[238,] 1.0000000 8.378779e-103 4.189390e-103
[239,] 1.0000000 4.883270e-102 2.441635e-102
[240,] 1.0000000 3.266965e-101 1.633482e-101
[241,] 1.0000000 2.175133e-100 1.087566e-100
[242,] 1.0000000 1.440290e-99 7.201452e-100
[243,] 1.0000000 9.497991e-99 4.748995e-99
[244,] 1.0000000 6.233213e-98 3.116607e-98
[245,] 1.0000000 4.062438e-97 2.031219e-97
[246,] 1.0000000 2.640528e-96 1.320264e-96
[247,] 1.0000000 1.682969e-95 8.414843e-96
[248,] 1.0000000 9.893016e-95 4.946508e-95
[249,] 1.0000000 1.772827e-94 8.864137e-95
[250,] 1.0000000 1.140099e-93 5.700497e-94
[251,] 1.0000000 7.060370e-93 3.530185e-93
[252,] 1.0000000 4.330226e-92 2.165113e-92
[253,] 1.0000000 2.744630e-91 1.372315e-91
[254,] 1.0000000 1.731084e-90 8.655422e-91
[255,] 1.0000000 1.086449e-89 5.432246e-90
[256,] 1.0000000 6.785047e-89 3.392523e-89
[257,] 1.0000000 3.966436e-88 1.983218e-88
[258,] 1.0000000 2.041567e-87 1.020783e-87
[259,] 1.0000000 1.257588e-86 6.287938e-87
[260,] 1.0000000 7.708079e-86 3.854040e-86
[261,] 1.0000000 4.700912e-85 2.350456e-85
[262,] 1.0000000 1.670528e-84 8.352639e-85
[263,] 1.0000000 1.011834e-83 5.059172e-84
[264,] 1.0000000 6.027321e-83 3.013661e-83
[265,] 1.0000000 3.614248e-82 1.807124e-82
[266,] 1.0000000 2.156316e-81 1.078158e-81
[267,] 1.0000000 1.279976e-80 6.399882e-81
[268,] 1.0000000 7.584774e-80 3.792387e-80
[269,] 1.0000000 4.456761e-79 2.228380e-79
[270,] 1.0000000 2.558711e-78 1.279355e-78
[271,] 1.0000000 1.471547e-77 7.357733e-78
[272,] 1.0000000 8.515274e-77 4.257637e-77
[273,] 1.0000000 4.651648e-76 2.325824e-76
[274,] 1.0000000 2.614948e-75 1.307474e-75
[275,] 1.0000000 1.491190e-74 7.455949e-75
[276,] 1.0000000 8.459391e-74 4.229696e-74
[277,] 1.0000000 4.645322e-73 2.322661e-73
[278,] 1.0000000 2.607978e-72 1.303989e-72
[279,] 1.0000000 1.438040e-71 7.190198e-72
[280,] 1.0000000 7.989222e-71 3.994611e-71
[281,] 1.0000000 4.194530e-70 2.097265e-70
[282,] 1.0000000 2.300912e-69 1.150456e-69
[283,] 1.0000000 1.258340e-68 6.291698e-69
[284,] 1.0000000 6.811215e-68 3.405608e-68
[285,] 1.0000000 3.685979e-67 1.842990e-67
[286,] 1.0000000 1.920092e-66 9.600461e-67
[287,] 1.0000000 8.540716e-66 4.270358e-66
[288,] 1.0000000 4.553659e-65 2.276830e-65
[289,] 1.0000000 2.406991e-64 1.203495e-64
[290,] 1.0000000 9.521619e-64 4.760809e-64
[291,] 1.0000000 5.003373e-63 2.501687e-63
[292,] 1.0000000 2.614645e-62 1.307323e-62
[293,] 1.0000000 1.358783e-61 6.793915e-62
[294,] 1.0000000 7.022058e-61 3.511029e-61
[295,] 1.0000000 3.620643e-60 1.810322e-60
[296,] 1.0000000 1.855116e-59 9.275579e-60
[297,] 1.0000000 9.427912e-59 4.713956e-59
[298,] 1.0000000 4.764232e-58 2.382116e-58
[299,] 1.0000000 1.827981e-57 9.139905e-58
[300,] 1.0000000 9.147311e-57 4.573655e-57
[301,] 1.0000000 4.551219e-56 2.275610e-56
[302,] 1.0000000 2.256528e-55 1.128264e-55
[303,] 1.0000000 1.072925e-54 5.364624e-55
[304,] 1.0000000 5.246372e-54 2.623186e-54
[305,] 1.0000000 2.550315e-53 1.275157e-53
[306,] 1.0000000 1.239058e-52 6.195288e-53
[307,] 1.0000000 5.952688e-52 2.976344e-52
[308,] 1.0000000 2.826933e-51 1.413467e-51
[309,] 1.0000000 1.341179e-50 6.705894e-51
[310,] 1.0000000 5.740491e-50 2.870246e-50
[311,] 1.0000000 2.686186e-49 1.343093e-49
[312,] 1.0000000 1.246388e-48 6.231939e-49
[313,] 1.0000000 5.779544e-48 2.889772e-48
[314,] 1.0000000 2.663502e-47 1.331751e-47
[315,] 1.0000000 1.219858e-46 6.099288e-47
[316,] 1.0000000 5.551942e-46 2.775971e-46
[317,] 1.0000000 2.510981e-45 1.255490e-45
[318,] 1.0000000 1.128462e-44 5.642310e-45
[319,] 1.0000000 5.039154e-44 2.519577e-44
[320,] 1.0000000 2.246771e-43 1.123386e-43
[321,] 1.0000000 9.904629e-43 4.952315e-43
[322,] 1.0000000 4.337989e-42 2.168995e-42
[323,] 1.0000000 1.887513e-41 9.437563e-42
[324,] 1.0000000 8.158725e-41 4.079363e-41
[325,] 1.0000000 3.503202e-40 1.751601e-40
[326,] 1.0000000 1.494161e-39 7.470806e-40
[327,] 1.0000000 6.329908e-39 3.164954e-39
[328,] 1.0000000 2.662968e-38 1.331484e-38
[329,] 1.0000000 1.112823e-37 5.564113e-38
[330,] 1.0000000 4.566980e-37 2.283490e-37
[331,] 1.0000000 1.893955e-36 9.469773e-37
[332,] 1.0000000 7.753152e-36 3.876576e-36
[333,] 1.0000000 3.161261e-35 1.580630e-35
[334,] 1.0000000 1.254796e-34 6.273982e-35
[335,] 1.0000000 4.985541e-34 2.492770e-34
[336,] 1.0000000 1.984781e-33 9.923907e-34
[337,] 1.0000000 7.843958e-33 3.921979e-33
[338,] 1.0000000 3.077171e-32 1.538586e-32
[339,] 1.0000000 1.196176e-31 5.980881e-32
[340,] 1.0000000 4.623043e-31 2.311521e-31
[341,] 1.0000000 1.790784e-30 8.953922e-31
[342,] 1.0000000 6.817450e-30 3.408725e-30
[343,] 1.0000000 2.575384e-29 1.287692e-29
[344,] 1.0000000 8.562889e-29 4.281444e-29
[345,] 1.0000000 3.067162e-28 1.533581e-28
[346,] 1.0000000 1.136078e-27 5.680389e-28
[347,] 1.0000000 3.801943e-27 1.900971e-27
[348,] 1.0000000 1.387979e-26 6.939897e-27
[349,] 1.0000000 5.005389e-26 2.502694e-26
[350,] 1.0000000 1.790094e-25 8.950468e-26
[351,] 1.0000000 6.420425e-25 3.210213e-25
[352,] 1.0000000 2.258778e-24 1.129389e-24
[353,] 1.0000000 7.878521e-24 3.939260e-24
[354,] 1.0000000 2.724158e-23 1.362079e-23
[355,] 1.0000000 9.336000e-23 4.668000e-23
[356,] 1.0000000 3.171325e-22 1.585663e-22
[357,] 1.0000000 1.067563e-21 5.337816e-22
[358,] 1.0000000 3.398774e-21 1.699387e-21
[359,] 1.0000000 1.106767e-20 5.533836e-21
[360,] 1.0000000 3.624978e-20 1.812489e-20
[361,] 1.0000000 1.176011e-19 5.880053e-20
[362,] 1.0000000 3.778459e-19 1.889229e-19
[363,] 1.0000000 1.202138e-18 6.010691e-19
[364,] 1.0000000 3.757288e-18 1.878644e-18
[365,] 1.0000000 1.171640e-17 5.858201e-18
[366,] 1.0000000 3.408895e-17 1.704447e-17
[367,] 1.0000000 1.050848e-16 5.254238e-17
[368,] 1.0000000 3.248350e-16 1.624175e-16
[369,] 1.0000000 9.722489e-16 4.861245e-16
[370,] 1.0000000 2.878293e-15 1.439146e-15
[371,] 1.0000000 8.598528e-15 4.299264e-15
[372,] 1.0000000 2.357267e-14 1.178634e-14
[373,] 1.0000000 5.835260e-14 2.917630e-14
[374,] 1.0000000 1.654031e-13 8.270153e-14
[375,] 1.0000000 4.632287e-13 2.316144e-13
[376,] 1.0000000 1.319653e-12 6.598263e-13
[377,] 1.0000000 3.516607e-12 1.758304e-12
[378,] 1.0000000 9.485563e-12 4.742781e-12
[379,] 1.0000000 2.532512e-11 1.266256e-11
[380,] 1.0000000 6.061433e-11 3.030716e-11
[381,] 1.0000000 1.492052e-10 7.460262e-11
[382,] 1.0000000 3.743283e-10 1.871641e-10
[383,] 1.0000000 9.495643e-10 4.747822e-10
[384,] 1.0000000 2.353804e-09 1.176902e-09
[385,] 1.0000000 5.656970e-09 2.828485e-09
[386,] 1.0000000 1.424074e-08 7.120370e-09
[387,] 1.0000000 3.504551e-08 1.752275e-08
[388,] 1.0000000 6.521758e-08 3.260879e-08
[389,] 0.9999999 1.564404e-07 7.822018e-08
[390,] 0.9999998 3.764172e-07 1.882086e-07
[391,] 0.9999996 8.549540e-07 4.274770e-07
[392,] 0.9999991 1.786796e-06 8.933981e-07
[393,] 0.9999979 4.106415e-06 2.053208e-06
[394,] 0.9999958 8.340300e-06 4.170150e-06
[395,] 0.9999907 1.855249e-05 9.276247e-06
[396,] 0.9999820 3.604372e-05 1.802186e-05
[397,] 0.9999624 7.529082e-05 3.764541e-05
[398,] 0.9999197 1.605082e-04 8.025408e-05
[399,] 0.9998323 3.353303e-04 1.676652e-04
[400,] 0.9996773 6.454466e-04 3.227233e-04
[401,] 0.9993797 1.240696e-03 6.203481e-04
[402,] 0.9988706 2.258706e-03 1.129353e-03
[403,] 0.9979484 4.103177e-03 2.051589e-03
[404,] 0.9963258 7.348326e-03 3.674163e-03
[405,] 0.9934303 1.313937e-02 6.569685e-03
[406,] 0.9897730 2.045392e-02 1.022696e-02
[407,] 0.9823029 3.539422e-02 1.769711e-02
[408,] 0.9861369 2.772617e-02 1.386308e-02
[409,] 0.9803969 3.920630e-02 1.960315e-02
[410,] 0.9684618 6.307648e-02 3.153824e-02
[411,] 0.9945327 1.093459e-02 5.467297e-03
[412,] 0.9907649 1.847019e-02 9.235095e-03
[413,] 0.9808068 3.838647e-02 1.919323e-02
[414,] 0.9610995 7.780094e-02 3.890047e-02
[415,] 0.9235743 1.528513e-01 7.642566e-02
[416,] 0.8551765 2.896469e-01 1.448235e-01
[417,] 0.8224317 3.551366e-01 1.775683e-01
[418,] 0.8561185 2.877631e-01 1.438815e-01
> postscript(file="/var/www/html/rcomp/tmp/12hau1292955824.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/22hau1292955824.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/3v89f1292955824.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/4v89f1292955824.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/5v89f1292955824.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
-2.608799e+05 3.411356e+06 1.171853e+06 -5.428272e+05 -2.608628e+05
6 7 8 9 10
1.933686e+04 -3.509149e+03 -2.425525e+05 2.364602e+04 8.821366e+05
11 12 13 14 15
-3.321008e+05 4.790465e+05 9.373054e+05 -6.683091e+05 7.445716e+04
16 17 18 19 20
5.724179e+05 -3.540487e+03 9.051730e+03 1.502868e+05 -3.216965e+05
21 22 23 24 25
2.868466e+05 -4.064967e+04 -2.352861e+04 -3.511651e+05 -2.018015e+05
26 27 28 29 30
-2.871585e+03 -1.831363e+05 -2.575756e+05 3.628291e+05 -9.667566e+04
31 32 33 34 35
3.120082e+05 -9.784852e+03 2.947249e+05 1.883540e+05 2.376275e+05
36 37 38 39 40
-3.559115e+05 1.154494e+05 2.621138e+05 1.665263e+05 -1.803472e+05
41 42 43 44 45
-2.325873e+05 -2.654010e+05 9.691298e+04 -9.742618e+04 -5.215692e+03
46 47 48 49 50
-1.283432e+04 -1.144812e+05 -8.124632e+04 3.854628e+04 -1.566947e+05
51 52 53 54 55
1.503060e+05 1.005303e+05 -6.448411e+04 1.912044e+05 -2.109276e+05
56 57 58 59 60
5.695869e+04 2.852243e+05 2.701651e+04 6.825481e+03 -1.145012e+05
61 62 63 64 65
6.538881e+04 -1.874398e+05 4.005412e+04 -9.580666e+04 1.824990e+03
66 67 68 69 70
-2.330799e+05 -6.286751e+04 4.498030e+04 3.280319e+04 2.753540e+04
71 72 73 74 75
-1.460492e+05 -3.612186e+03 -2.811354e+04 8.244743e+03 4.011274e+04
76 77 78 79 80
-9.423110e+04 3.847241e+04 -8.553179e+04 -1.837766e+04 1.870321e+04
81 82 83 84 85
1.592131e+04 2.028738e+04 2.503243e+04 -1.852851e+04 6.045534e+04
86 87 88 89 90
-5.147537e+04 7.445046e+04 3.654182e+04 -4.602990e+04 6.791020e+04
91 92 93 94 95
4.344371e+04 1.205713e+04 1.051645e+05 -1.120634e+05 2.316142e+04
96 97 98 99 100
1.633573e+04 1.178379e+05 -2.796317e+03 3.084612e+04 -1.528409e+04
101 102 103 104 105
-2.366499e+04 -9.163387e+04 -8.620064e+04 -4.688865e+03 3.003454e+04
106 107 108 109 110
-5.128738e+04 -2.371219e+05 1.594284e+05 -1.143577e+05 -9.369199e+04
111 112 113 114 115
8.462011e+04 -2.032839e+05 -1.209303e+05 -6.322223e+04 4.690003e+04
116 117 118 119 120
1.056765e+04 1.827586e+04 8.947695e+04 -1.361990e+05 -6.355304e+03
121 122 123 124 125
-1.855450e+05 1.119938e+04 1.558382e+04 1.359708e+05 -1.027438e+03
126 127 128 129 130
-1.534732e+05 7.885890e+03 5.411699e+03 3.250702e+03 7.885890e+03
131 132 133 134 135
6.264633e+03 3.574508e+04 7.885890e+03 4.614398e+03 -6.605722e+03
136 137 138 139 140
2.031769e+03 9.423296e+03 7.306172e+03 -8.138305e+03 7.885890e+03
141 142 143 144 145
1.126219e+04 7.885890e+03 -2.597494e+04 7.724318e+03 1.429069e+04
146 147 148 149 150
-2.892442e+01 7.885890e+03 7.343041e+03 1.305478e+04 8.155946e+03
151 152 153 154 155
-2.482893e+04 -3.211069e+04 7.885890e+03 -1.765617e+04 7.415522e+03
156 157 158 159 160
7.885890e+03 7.885890e+03 6.074550e+03 2.190668e+04 7.885890e+03
161 162 163 164 165
7.885890e+03 7.885890e+03 7.885890e+03 -9.604297e+03 -1.319946e+04
166 167 168 169 170
6.572742e+03 7.885890e+03 6.282702e+04 -4.665849e+04 7.885890e+03
171 172 173 174 175
7.885890e+03 4.336193e+03 7.885890e+03 -1.385396e+03 7.885890e+03
176 177 178 179 180
5.876404e+03 -1.164544e+03 7.451639e+03 7.885890e+03 7.885890e+03
181 182 183 184 185
3.412188e+04 1.316928e+04 -3.448478e+04 -1.660040e+04 7.009061e+03
186 187 188 189 190
1.466212e+04 3.633812e+03 3.299849e+03 -1.950469e+04 7.885890e+03
191 192 193 194 195
6.460031e+03 7.885890e+03 5.459244e+03 7.885890e+03 7.717053e+03
196 197 198 199 200
7.885890e+03 -1.764143e+04 7.885890e+03 7.885890e+03 -4.628102e+04
201 202 203 204 205
7.885890e+03 7.885890e+03 7.885890e+03 -1.032528e+04 9.034898e+03
206 207 208 209 210
8.538142e+03 1.021988e+04 7.885890e+03 1.765476e+04 -5.142160e+03
211 212 213 214 215
7.885890e+03 7.885890e+03 1.787679e+04 -2.674253e+04 -1.230533e+04
216 217 218 219 220
7.885890e+03 7.939671e+03 3.356439e+03 7.885890e+03 -1.661634e+04
221 222 223 224 225
-1.583636e+04 7.885890e+03 7.806647e+03 -3.449907e+03 7.739504e+03
226 227 228 229 230
2.618073e+03 5.998931e+03 7.885890e+03 1.647015e+03 7.885890e+03
231 232 233 234 235
-1.425460e+04 1.082325e+03 7.885890e+03 -1.064204e+04 7.885890e+03
236 237 238 239 240
7.959474e+03 7.885890e+03 7.361795e+03 7.885890e+03 7.885890e+03
241 242 243 244 245
7.885890e+03 3.789819e+03 7.885890e+03 7.679544e+03 -2.674581e+04
246 247 248 249 250
7.885890e+03 7.885890e+03 -7.364000e+02 7.885890e+03 7.885890e+03
251 252 253 254 255
-1.457080e+04 7.885890e+03 -2.954105e+03 -1.427997e+04 -1.001436e+05
256 257 258 259 260
7.885890e+03 -4.294474e+03 -1.540906e+04 7.885890e+03 7.885890e+03
261 262 263 264 265
7.885890e+03 7.885890e+03 3.280955e+04 1.911200e+04 7.885890e+03
266 267 268 269 270
7.885890e+03 7.885890e+03 7.092566e+04 7.885890e+03 2.186414e+04
271 272 273 274 275
7.885890e+03 7.885890e+03 7.885890e+03 5.419260e+02 7.885890e+03
276 277 278 279 280
7.659666e+03 7.125660e+03 7.885890e+03 -3.947279e+04 3.548361e+02
281 282 283 284 285
7.885890e+03 7.885890e+03 1.594011e+03 7.885890e+03 1.814935e+04
286 287 288 289 290
7.885890e+03 1.292396e+03 1.015167e+04 7.885890e+03 -1.919372e+04
291 292 293 294 295
7.885890e+03 1.488961e+04 5.046362e+04 7.885890e+03 -6.062845e+03
296 297 298 299 300
2.684495e+04 7.885890e+03 7.885890e+03 7.885890e+03 7.885890e+03
301 302 303 304 305
5.153591e+03 -1.138956e+04 7.885890e+03 7.885890e+03 3.058851e+04
306 307 308 309 310
7.866032e+03 7.885890e+03 5.615897e+03 -4.332841e+03 7.885890e+03
311 312 313 314 315
7.885890e+03 3.292089e+04 7.885890e+03 1.385221e+04 8.112119e+03
316 317 318 319 320
2.983997e+04 6.877481e+03 -1.673346e+04 7.860491e+03 7.885890e+03
321 322 323 324 325
7.885890e+03 7.885890e+03 7.885890e+03 7.885890e+03 7.885890e+03
326 327 328 329 330
4.626436e+03 7.885890e+03 7.885890e+03 7.885890e+03 7.885890e+03
331 332 333 334 335
7.885890e+03 7.885890e+03 7.885890e+03 4.555494e+03 7.876934e+03
336 337 338 339 340
2.431005e+03 -2.600563e+04 7.885890e+03 5.944099e+03 7.998222e+03
341 342 343 344 345
-1.715618e+04 7.885890e+03 7.885890e+03 7.885890e+03 -1.000848e+04
346 347 348 349 350
7.885890e+03 -5.943610e+03 7.885890e+03 7.885890e+03 -1.578813e+04
351 352 353 354 355
5.758370e+03 5.533360e+03 8.053020e+03 4.424100e+03 7.885890e+03
356 357 358 359 360
7.885890e+03 5.688983e+03 7.885890e+03 7.885890e+03 7.885890e+03
361 362 363 364 365
7.886657e+03 7.885890e+03 7.885890e+03 -2.990045e+04 5.736787e+03
366 367 368 369 370
7.885890e+03 7.885890e+03 7.885890e+03 7.885890e+03 4.277989e+03
371 372 373 374 375
7.885890e+03 8.102971e+03 2.559254e+04 -2.958819e+03 7.885890e+03
376 377 378 379 380
7.885890e+03 -1.439895e+04 4.326394e+03 -2.599412e+04 7.885890e+03
381 382 383 384 385
7.885890e+03 3.793556e+02 -3.654900e+04 7.885890e+03 8.571701e+03
386 387 388 389 390
-1.145814e+05 1.726525e+04 -1.431291e+04 5.452475e+03 -9.434267e+03
391 392 393 394 395
7.948788e+03 -2.315341e+04 -5.501282e+04 4.782787e+04 -9.255679e+04
396 397 398 399 400
-3.360404e+04 -1.441857e+05 9.729573e+04 -2.046329e+04 -1.042779e+05
401 402 403 404 405
-2.252725e+04 1.237066e+04 -2.112778e+04 -5.631186e+04 -4.095689e+04
406 407 408 409 410
-2.063086e+03 -1.439303e+03 -1.994659e+04 -1.698257e+04 -8.863966e+04
411 412 413 414 415
-2.180996e+04 -1.516623e+05 -5.757076e+04 7.750794e+04 -1.793147e+05
416 417 418 419 420
-1.209000e+05 5.324414e+04 -3.686172e+04 -5.522805e+04 -7.260360e+04
421 422 423 424 425
-7.235992e+04 -1.199613e+05 -4.749294e+05 -4.586813e+05 -7.261920e+04
426 427 428 429 430
-1.848519e+05 -3.293887e+05 -1.089227e+06 -2.374091e+05 -2.333980e+05
431
-5.464335e+05
> postscript(file="/var/www/html/rcomp/tmp/6o0801292955824.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 -2.608799e+05 NA
1 3.411356e+06 -2.608799e+05
2 1.171853e+06 3.411356e+06
3 -5.428272e+05 1.171853e+06
4 -2.608628e+05 -5.428272e+05
5 1.933686e+04 -2.608628e+05
6 -3.509149e+03 1.933686e+04
7 -2.425525e+05 -3.509149e+03
8 2.364602e+04 -2.425525e+05
9 8.821366e+05 2.364602e+04
10 -3.321008e+05 8.821366e+05
11 4.790465e+05 -3.321008e+05
12 9.373054e+05 4.790465e+05
13 -6.683091e+05 9.373054e+05
14 7.445716e+04 -6.683091e+05
15 5.724179e+05 7.445716e+04
16 -3.540487e+03 5.724179e+05
17 9.051730e+03 -3.540487e+03
18 1.502868e+05 9.051730e+03
19 -3.216965e+05 1.502868e+05
20 2.868466e+05 -3.216965e+05
21 -4.064967e+04 2.868466e+05
22 -2.352861e+04 -4.064967e+04
23 -3.511651e+05 -2.352861e+04
24 -2.018015e+05 -3.511651e+05
25 -2.871585e+03 -2.018015e+05
26 -1.831363e+05 -2.871585e+03
27 -2.575756e+05 -1.831363e+05
28 3.628291e+05 -2.575756e+05
29 -9.667566e+04 3.628291e+05
30 3.120082e+05 -9.667566e+04
31 -9.784852e+03 3.120082e+05
32 2.947249e+05 -9.784852e+03
33 1.883540e+05 2.947249e+05
34 2.376275e+05 1.883540e+05
35 -3.559115e+05 2.376275e+05
36 1.154494e+05 -3.559115e+05
37 2.621138e+05 1.154494e+05
38 1.665263e+05 2.621138e+05
39 -1.803472e+05 1.665263e+05
40 -2.325873e+05 -1.803472e+05
41 -2.654010e+05 -2.325873e+05
42 9.691298e+04 -2.654010e+05
43 -9.742618e+04 9.691298e+04
44 -5.215692e+03 -9.742618e+04
45 -1.283432e+04 -5.215692e+03
46 -1.144812e+05 -1.283432e+04
47 -8.124632e+04 -1.144812e+05
48 3.854628e+04 -8.124632e+04
49 -1.566947e+05 3.854628e+04
50 1.503060e+05 -1.566947e+05
51 1.005303e+05 1.503060e+05
52 -6.448411e+04 1.005303e+05
53 1.912044e+05 -6.448411e+04
54 -2.109276e+05 1.912044e+05
55 5.695869e+04 -2.109276e+05
56 2.852243e+05 5.695869e+04
57 2.701651e+04 2.852243e+05
58 6.825481e+03 2.701651e+04
59 -1.145012e+05 6.825481e+03
60 6.538881e+04 -1.145012e+05
61 -1.874398e+05 6.538881e+04
62 4.005412e+04 -1.874398e+05
63 -9.580666e+04 4.005412e+04
64 1.824990e+03 -9.580666e+04
65 -2.330799e+05 1.824990e+03
66 -6.286751e+04 -2.330799e+05
67 4.498030e+04 -6.286751e+04
68 3.280319e+04 4.498030e+04
69 2.753540e+04 3.280319e+04
70 -1.460492e+05 2.753540e+04
71 -3.612186e+03 -1.460492e+05
72 -2.811354e+04 -3.612186e+03
73 8.244743e+03 -2.811354e+04
74 4.011274e+04 8.244743e+03
75 -9.423110e+04 4.011274e+04
76 3.847241e+04 -9.423110e+04
77 -8.553179e+04 3.847241e+04
78 -1.837766e+04 -8.553179e+04
79 1.870321e+04 -1.837766e+04
80 1.592131e+04 1.870321e+04
81 2.028738e+04 1.592131e+04
82 2.503243e+04 2.028738e+04
83 -1.852851e+04 2.503243e+04
84 6.045534e+04 -1.852851e+04
85 -5.147537e+04 6.045534e+04
86 7.445046e+04 -5.147537e+04
87 3.654182e+04 7.445046e+04
88 -4.602990e+04 3.654182e+04
89 6.791020e+04 -4.602990e+04
90 4.344371e+04 6.791020e+04
91 1.205713e+04 4.344371e+04
92 1.051645e+05 1.205713e+04
93 -1.120634e+05 1.051645e+05
94 2.316142e+04 -1.120634e+05
95 1.633573e+04 2.316142e+04
96 1.178379e+05 1.633573e+04
97 -2.796317e+03 1.178379e+05
98 3.084612e+04 -2.796317e+03
99 -1.528409e+04 3.084612e+04
100 -2.366499e+04 -1.528409e+04
101 -9.163387e+04 -2.366499e+04
102 -8.620064e+04 -9.163387e+04
103 -4.688865e+03 -8.620064e+04
104 3.003454e+04 -4.688865e+03
105 -5.128738e+04 3.003454e+04
106 -2.371219e+05 -5.128738e+04
107 1.594284e+05 -2.371219e+05
108 -1.143577e+05 1.594284e+05
109 -9.369199e+04 -1.143577e+05
110 8.462011e+04 -9.369199e+04
111 -2.032839e+05 8.462011e+04
112 -1.209303e+05 -2.032839e+05
113 -6.322223e+04 -1.209303e+05
114 4.690003e+04 -6.322223e+04
115 1.056765e+04 4.690003e+04
116 1.827586e+04 1.056765e+04
117 8.947695e+04 1.827586e+04
118 -1.361990e+05 8.947695e+04
119 -6.355304e+03 -1.361990e+05
120 -1.855450e+05 -6.355304e+03
121 1.119938e+04 -1.855450e+05
122 1.558382e+04 1.119938e+04
123 1.359708e+05 1.558382e+04
124 -1.027438e+03 1.359708e+05
125 -1.534732e+05 -1.027438e+03
126 7.885890e+03 -1.534732e+05
127 5.411699e+03 7.885890e+03
128 3.250702e+03 5.411699e+03
129 7.885890e+03 3.250702e+03
130 6.264633e+03 7.885890e+03
131 3.574508e+04 6.264633e+03
132 7.885890e+03 3.574508e+04
133 4.614398e+03 7.885890e+03
134 -6.605722e+03 4.614398e+03
135 2.031769e+03 -6.605722e+03
136 9.423296e+03 2.031769e+03
137 7.306172e+03 9.423296e+03
138 -8.138305e+03 7.306172e+03
139 7.885890e+03 -8.138305e+03
140 1.126219e+04 7.885890e+03
141 7.885890e+03 1.126219e+04
142 -2.597494e+04 7.885890e+03
143 7.724318e+03 -2.597494e+04
144 1.429069e+04 7.724318e+03
145 -2.892442e+01 1.429069e+04
146 7.885890e+03 -2.892442e+01
147 7.343041e+03 7.885890e+03
148 1.305478e+04 7.343041e+03
149 8.155946e+03 1.305478e+04
150 -2.482893e+04 8.155946e+03
151 -3.211069e+04 -2.482893e+04
152 7.885890e+03 -3.211069e+04
153 -1.765617e+04 7.885890e+03
154 7.415522e+03 -1.765617e+04
155 7.885890e+03 7.415522e+03
156 7.885890e+03 7.885890e+03
157 6.074550e+03 7.885890e+03
158 2.190668e+04 6.074550e+03
159 7.885890e+03 2.190668e+04
160 7.885890e+03 7.885890e+03
161 7.885890e+03 7.885890e+03
162 7.885890e+03 7.885890e+03
163 -9.604297e+03 7.885890e+03
164 -1.319946e+04 -9.604297e+03
165 6.572742e+03 -1.319946e+04
166 7.885890e+03 6.572742e+03
167 6.282702e+04 7.885890e+03
168 -4.665849e+04 6.282702e+04
169 7.885890e+03 -4.665849e+04
170 7.885890e+03 7.885890e+03
171 4.336193e+03 7.885890e+03
172 7.885890e+03 4.336193e+03
173 -1.385396e+03 7.885890e+03
174 7.885890e+03 -1.385396e+03
175 5.876404e+03 7.885890e+03
176 -1.164544e+03 5.876404e+03
177 7.451639e+03 -1.164544e+03
178 7.885890e+03 7.451639e+03
179 7.885890e+03 7.885890e+03
180 3.412188e+04 7.885890e+03
181 1.316928e+04 3.412188e+04
182 -3.448478e+04 1.316928e+04
183 -1.660040e+04 -3.448478e+04
184 7.009061e+03 -1.660040e+04
185 1.466212e+04 7.009061e+03
186 3.633812e+03 1.466212e+04
187 3.299849e+03 3.633812e+03
188 -1.950469e+04 3.299849e+03
189 7.885890e+03 -1.950469e+04
190 6.460031e+03 7.885890e+03
191 7.885890e+03 6.460031e+03
192 5.459244e+03 7.885890e+03
193 7.885890e+03 5.459244e+03
194 7.717053e+03 7.885890e+03
195 7.885890e+03 7.717053e+03
196 -1.764143e+04 7.885890e+03
197 7.885890e+03 -1.764143e+04
198 7.885890e+03 7.885890e+03
199 -4.628102e+04 7.885890e+03
200 7.885890e+03 -4.628102e+04
201 7.885890e+03 7.885890e+03
202 7.885890e+03 7.885890e+03
203 -1.032528e+04 7.885890e+03
204 9.034898e+03 -1.032528e+04
205 8.538142e+03 9.034898e+03
206 1.021988e+04 8.538142e+03
207 7.885890e+03 1.021988e+04
208 1.765476e+04 7.885890e+03
209 -5.142160e+03 1.765476e+04
210 7.885890e+03 -5.142160e+03
211 7.885890e+03 7.885890e+03
212 1.787679e+04 7.885890e+03
213 -2.674253e+04 1.787679e+04
214 -1.230533e+04 -2.674253e+04
215 7.885890e+03 -1.230533e+04
216 7.939671e+03 7.885890e+03
217 3.356439e+03 7.939671e+03
218 7.885890e+03 3.356439e+03
219 -1.661634e+04 7.885890e+03
220 -1.583636e+04 -1.661634e+04
221 7.885890e+03 -1.583636e+04
222 7.806647e+03 7.885890e+03
223 -3.449907e+03 7.806647e+03
224 7.739504e+03 -3.449907e+03
225 2.618073e+03 7.739504e+03
226 5.998931e+03 2.618073e+03
227 7.885890e+03 5.998931e+03
228 1.647015e+03 7.885890e+03
229 7.885890e+03 1.647015e+03
230 -1.425460e+04 7.885890e+03
231 1.082325e+03 -1.425460e+04
232 7.885890e+03 1.082325e+03
233 -1.064204e+04 7.885890e+03
234 7.885890e+03 -1.064204e+04
235 7.959474e+03 7.885890e+03
236 7.885890e+03 7.959474e+03
237 7.361795e+03 7.885890e+03
238 7.885890e+03 7.361795e+03
239 7.885890e+03 7.885890e+03
240 7.885890e+03 7.885890e+03
241 3.789819e+03 7.885890e+03
242 7.885890e+03 3.789819e+03
243 7.679544e+03 7.885890e+03
244 -2.674581e+04 7.679544e+03
245 7.885890e+03 -2.674581e+04
246 7.885890e+03 7.885890e+03
247 -7.364000e+02 7.885890e+03
248 7.885890e+03 -7.364000e+02
249 7.885890e+03 7.885890e+03
250 -1.457080e+04 7.885890e+03
251 7.885890e+03 -1.457080e+04
252 -2.954105e+03 7.885890e+03
253 -1.427997e+04 -2.954105e+03
254 -1.001436e+05 -1.427997e+04
255 7.885890e+03 -1.001436e+05
256 -4.294474e+03 7.885890e+03
257 -1.540906e+04 -4.294474e+03
258 7.885890e+03 -1.540906e+04
259 7.885890e+03 7.885890e+03
260 7.885890e+03 7.885890e+03
261 7.885890e+03 7.885890e+03
262 3.280955e+04 7.885890e+03
263 1.911200e+04 3.280955e+04
264 7.885890e+03 1.911200e+04
265 7.885890e+03 7.885890e+03
266 7.885890e+03 7.885890e+03
267 7.092566e+04 7.885890e+03
268 7.885890e+03 7.092566e+04
269 2.186414e+04 7.885890e+03
270 7.885890e+03 2.186414e+04
271 7.885890e+03 7.885890e+03
272 7.885890e+03 7.885890e+03
273 5.419260e+02 7.885890e+03
274 7.885890e+03 5.419260e+02
275 7.659666e+03 7.885890e+03
276 7.125660e+03 7.659666e+03
277 7.885890e+03 7.125660e+03
278 -3.947279e+04 7.885890e+03
279 3.548361e+02 -3.947279e+04
280 7.885890e+03 3.548361e+02
281 7.885890e+03 7.885890e+03
282 1.594011e+03 7.885890e+03
283 7.885890e+03 1.594011e+03
284 1.814935e+04 7.885890e+03
285 7.885890e+03 1.814935e+04
286 1.292396e+03 7.885890e+03
287 1.015167e+04 1.292396e+03
288 7.885890e+03 1.015167e+04
289 -1.919372e+04 7.885890e+03
290 7.885890e+03 -1.919372e+04
291 1.488961e+04 7.885890e+03
292 5.046362e+04 1.488961e+04
293 7.885890e+03 5.046362e+04
294 -6.062845e+03 7.885890e+03
295 2.684495e+04 -6.062845e+03
296 7.885890e+03 2.684495e+04
297 7.885890e+03 7.885890e+03
298 7.885890e+03 7.885890e+03
299 7.885890e+03 7.885890e+03
300 5.153591e+03 7.885890e+03
301 -1.138956e+04 5.153591e+03
302 7.885890e+03 -1.138956e+04
303 7.885890e+03 7.885890e+03
304 3.058851e+04 7.885890e+03
305 7.866032e+03 3.058851e+04
306 7.885890e+03 7.866032e+03
307 5.615897e+03 7.885890e+03
308 -4.332841e+03 5.615897e+03
309 7.885890e+03 -4.332841e+03
310 7.885890e+03 7.885890e+03
311 3.292089e+04 7.885890e+03
312 7.885890e+03 3.292089e+04
313 1.385221e+04 7.885890e+03
314 8.112119e+03 1.385221e+04
315 2.983997e+04 8.112119e+03
316 6.877481e+03 2.983997e+04
317 -1.673346e+04 6.877481e+03
318 7.860491e+03 -1.673346e+04
319 7.885890e+03 7.860491e+03
320 7.885890e+03 7.885890e+03
321 7.885890e+03 7.885890e+03
322 7.885890e+03 7.885890e+03
323 7.885890e+03 7.885890e+03
324 7.885890e+03 7.885890e+03
325 4.626436e+03 7.885890e+03
326 7.885890e+03 4.626436e+03
327 7.885890e+03 7.885890e+03
328 7.885890e+03 7.885890e+03
329 7.885890e+03 7.885890e+03
330 7.885890e+03 7.885890e+03
331 7.885890e+03 7.885890e+03
332 7.885890e+03 7.885890e+03
333 4.555494e+03 7.885890e+03
334 7.876934e+03 4.555494e+03
335 2.431005e+03 7.876934e+03
336 -2.600563e+04 2.431005e+03
337 7.885890e+03 -2.600563e+04
338 5.944099e+03 7.885890e+03
339 7.998222e+03 5.944099e+03
340 -1.715618e+04 7.998222e+03
341 7.885890e+03 -1.715618e+04
342 7.885890e+03 7.885890e+03
343 7.885890e+03 7.885890e+03
344 -1.000848e+04 7.885890e+03
345 7.885890e+03 -1.000848e+04
346 -5.943610e+03 7.885890e+03
347 7.885890e+03 -5.943610e+03
348 7.885890e+03 7.885890e+03
349 -1.578813e+04 7.885890e+03
350 5.758370e+03 -1.578813e+04
351 5.533360e+03 5.758370e+03
352 8.053020e+03 5.533360e+03
353 4.424100e+03 8.053020e+03
354 7.885890e+03 4.424100e+03
355 7.885890e+03 7.885890e+03
356 5.688983e+03 7.885890e+03
357 7.885890e+03 5.688983e+03
358 7.885890e+03 7.885890e+03
359 7.885890e+03 7.885890e+03
360 7.886657e+03 7.885890e+03
361 7.885890e+03 7.886657e+03
362 7.885890e+03 7.885890e+03
363 -2.990045e+04 7.885890e+03
364 5.736787e+03 -2.990045e+04
365 7.885890e+03 5.736787e+03
366 7.885890e+03 7.885890e+03
367 7.885890e+03 7.885890e+03
368 7.885890e+03 7.885890e+03
369 4.277989e+03 7.885890e+03
370 7.885890e+03 4.277989e+03
371 8.102971e+03 7.885890e+03
372 2.559254e+04 8.102971e+03
373 -2.958819e+03 2.559254e+04
374 7.885890e+03 -2.958819e+03
375 7.885890e+03 7.885890e+03
376 -1.439895e+04 7.885890e+03
377 4.326394e+03 -1.439895e+04
378 -2.599412e+04 4.326394e+03
379 7.885890e+03 -2.599412e+04
380 7.885890e+03 7.885890e+03
381 3.793556e+02 7.885890e+03
382 -3.654900e+04 3.793556e+02
383 7.885890e+03 -3.654900e+04
384 8.571701e+03 7.885890e+03
385 -1.145814e+05 8.571701e+03
386 1.726525e+04 -1.145814e+05
387 -1.431291e+04 1.726525e+04
388 5.452475e+03 -1.431291e+04
389 -9.434267e+03 5.452475e+03
390 7.948788e+03 -9.434267e+03
391 -2.315341e+04 7.948788e+03
392 -5.501282e+04 -2.315341e+04
393 4.782787e+04 -5.501282e+04
394 -9.255679e+04 4.782787e+04
395 -3.360404e+04 -9.255679e+04
396 -1.441857e+05 -3.360404e+04
397 9.729573e+04 -1.441857e+05
398 -2.046329e+04 9.729573e+04
399 -1.042779e+05 -2.046329e+04
400 -2.252725e+04 -1.042779e+05
401 1.237066e+04 -2.252725e+04
402 -2.112778e+04 1.237066e+04
403 -5.631186e+04 -2.112778e+04
404 -4.095689e+04 -5.631186e+04
405 -2.063086e+03 -4.095689e+04
406 -1.439303e+03 -2.063086e+03
407 -1.994659e+04 -1.439303e+03
408 -1.698257e+04 -1.994659e+04
409 -8.863966e+04 -1.698257e+04
410 -2.180996e+04 -8.863966e+04
411 -1.516623e+05 -2.180996e+04
412 -5.757076e+04 -1.516623e+05
413 7.750794e+04 -5.757076e+04
414 -1.793147e+05 7.750794e+04
415 -1.209000e+05 -1.793147e+05
416 5.324414e+04 -1.209000e+05
417 -3.686172e+04 5.324414e+04
418 -5.522805e+04 -3.686172e+04
419 -7.260360e+04 -5.522805e+04
420 -7.235992e+04 -7.260360e+04
421 -1.199613e+05 -7.235992e+04
422 -4.749294e+05 -1.199613e+05
423 -4.586813e+05 -4.749294e+05
424 -7.261920e+04 -4.586813e+05
425 -1.848519e+05 -7.261920e+04
426 -3.293887e+05 -1.848519e+05
427 -1.089227e+06 -3.293887e+05
428 -2.374091e+05 -1.089227e+06
429 -2.333980e+05 -2.374091e+05
430 -5.464335e+05 -2.333980e+05
431 NA -5.464335e+05
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3.411356e+06 -2.608799e+05
[2,] 1.171853e+06 3.411356e+06
[3,] -5.428272e+05 1.171853e+06
[4,] -2.608628e+05 -5.428272e+05
[5,] 1.933686e+04 -2.608628e+05
[6,] -3.509149e+03 1.933686e+04
[7,] -2.425525e+05 -3.509149e+03
[8,] 2.364602e+04 -2.425525e+05
[9,] 8.821366e+05 2.364602e+04
[10,] -3.321008e+05 8.821366e+05
[11,] 4.790465e+05 -3.321008e+05
[12,] 9.373054e+05 4.790465e+05
[13,] -6.683091e+05 9.373054e+05
[14,] 7.445716e+04 -6.683091e+05
[15,] 5.724179e+05 7.445716e+04
[16,] -3.540487e+03 5.724179e+05
[17,] 9.051730e+03 -3.540487e+03
[18,] 1.502868e+05 9.051730e+03
[19,] -3.216965e+05 1.502868e+05
[20,] 2.868466e+05 -3.216965e+05
[21,] -4.064967e+04 2.868466e+05
[22,] -2.352861e+04 -4.064967e+04
[23,] -3.511651e+05 -2.352861e+04
[24,] -2.018015e+05 -3.511651e+05
[25,] -2.871585e+03 -2.018015e+05
[26,] -1.831363e+05 -2.871585e+03
[27,] -2.575756e+05 -1.831363e+05
[28,] 3.628291e+05 -2.575756e+05
[29,] -9.667566e+04 3.628291e+05
[30,] 3.120082e+05 -9.667566e+04
[31,] -9.784852e+03 3.120082e+05
[32,] 2.947249e+05 -9.784852e+03
[33,] 1.883540e+05 2.947249e+05
[34,] 2.376275e+05 1.883540e+05
[35,] -3.559115e+05 2.376275e+05
[36,] 1.154494e+05 -3.559115e+05
[37,] 2.621138e+05 1.154494e+05
[38,] 1.665263e+05 2.621138e+05
[39,] -1.803472e+05 1.665263e+05
[40,] -2.325873e+05 -1.803472e+05
[41,] -2.654010e+05 -2.325873e+05
[42,] 9.691298e+04 -2.654010e+05
[43,] -9.742618e+04 9.691298e+04
[44,] -5.215692e+03 -9.742618e+04
[45,] -1.283432e+04 -5.215692e+03
[46,] -1.144812e+05 -1.283432e+04
[47,] -8.124632e+04 -1.144812e+05
[48,] 3.854628e+04 -8.124632e+04
[49,] -1.566947e+05 3.854628e+04
[50,] 1.503060e+05 -1.566947e+05
[51,] 1.005303e+05 1.503060e+05
[52,] -6.448411e+04 1.005303e+05
[53,] 1.912044e+05 -6.448411e+04
[54,] -2.109276e+05 1.912044e+05
[55,] 5.695869e+04 -2.109276e+05
[56,] 2.852243e+05 5.695869e+04
[57,] 2.701651e+04 2.852243e+05
[58,] 6.825481e+03 2.701651e+04
[59,] -1.145012e+05 6.825481e+03
[60,] 6.538881e+04 -1.145012e+05
[61,] -1.874398e+05 6.538881e+04
[62,] 4.005412e+04 -1.874398e+05
[63,] -9.580666e+04 4.005412e+04
[64,] 1.824990e+03 -9.580666e+04
[65,] -2.330799e+05 1.824990e+03
[66,] -6.286751e+04 -2.330799e+05
[67,] 4.498030e+04 -6.286751e+04
[68,] 3.280319e+04 4.498030e+04
[69,] 2.753540e+04 3.280319e+04
[70,] -1.460492e+05 2.753540e+04
[71,] -3.612186e+03 -1.460492e+05
[72,] -2.811354e+04 -3.612186e+03
[73,] 8.244743e+03 -2.811354e+04
[74,] 4.011274e+04 8.244743e+03
[75,] -9.423110e+04 4.011274e+04
[76,] 3.847241e+04 -9.423110e+04
[77,] -8.553179e+04 3.847241e+04
[78,] -1.837766e+04 -8.553179e+04
[79,] 1.870321e+04 -1.837766e+04
[80,] 1.592131e+04 1.870321e+04
[81,] 2.028738e+04 1.592131e+04
[82,] 2.503243e+04 2.028738e+04
[83,] -1.852851e+04 2.503243e+04
[84,] 6.045534e+04 -1.852851e+04
[85,] -5.147537e+04 6.045534e+04
[86,] 7.445046e+04 -5.147537e+04
[87,] 3.654182e+04 7.445046e+04
[88,] -4.602990e+04 3.654182e+04
[89,] 6.791020e+04 -4.602990e+04
[90,] 4.344371e+04 6.791020e+04
[91,] 1.205713e+04 4.344371e+04
[92,] 1.051645e+05 1.205713e+04
[93,] -1.120634e+05 1.051645e+05
[94,] 2.316142e+04 -1.120634e+05
[95,] 1.633573e+04 2.316142e+04
[96,] 1.178379e+05 1.633573e+04
[97,] -2.796317e+03 1.178379e+05
[98,] 3.084612e+04 -2.796317e+03
[99,] -1.528409e+04 3.084612e+04
[100,] -2.366499e+04 -1.528409e+04
[101,] -9.163387e+04 -2.366499e+04
[102,] -8.620064e+04 -9.163387e+04
[103,] -4.688865e+03 -8.620064e+04
[104,] 3.003454e+04 -4.688865e+03
[105,] -5.128738e+04 3.003454e+04
[106,] -2.371219e+05 -5.128738e+04
[107,] 1.594284e+05 -2.371219e+05
[108,] -1.143577e+05 1.594284e+05
[109,] -9.369199e+04 -1.143577e+05
[110,] 8.462011e+04 -9.369199e+04
[111,] -2.032839e+05 8.462011e+04
[112,] -1.209303e+05 -2.032839e+05
[113,] -6.322223e+04 -1.209303e+05
[114,] 4.690003e+04 -6.322223e+04
[115,] 1.056765e+04 4.690003e+04
[116,] 1.827586e+04 1.056765e+04
[117,] 8.947695e+04 1.827586e+04
[118,] -1.361990e+05 8.947695e+04
[119,] -6.355304e+03 -1.361990e+05
[120,] -1.855450e+05 -6.355304e+03
[121,] 1.119938e+04 -1.855450e+05
[122,] 1.558382e+04 1.119938e+04
[123,] 1.359708e+05 1.558382e+04
[124,] -1.027438e+03 1.359708e+05
[125,] -1.534732e+05 -1.027438e+03
[126,] 7.885890e+03 -1.534732e+05
[127,] 5.411699e+03 7.885890e+03
[128,] 3.250702e+03 5.411699e+03
[129,] 7.885890e+03 3.250702e+03
[130,] 6.264633e+03 7.885890e+03
[131,] 3.574508e+04 6.264633e+03
[132,] 7.885890e+03 3.574508e+04
[133,] 4.614398e+03 7.885890e+03
[134,] -6.605722e+03 4.614398e+03
[135,] 2.031769e+03 -6.605722e+03
[136,] 9.423296e+03 2.031769e+03
[137,] 7.306172e+03 9.423296e+03
[138,] -8.138305e+03 7.306172e+03
[139,] 7.885890e+03 -8.138305e+03
[140,] 1.126219e+04 7.885890e+03
[141,] 7.885890e+03 1.126219e+04
[142,] -2.597494e+04 7.885890e+03
[143,] 7.724318e+03 -2.597494e+04
[144,] 1.429069e+04 7.724318e+03
[145,] -2.892442e+01 1.429069e+04
[146,] 7.885890e+03 -2.892442e+01
[147,] 7.343041e+03 7.885890e+03
[148,] 1.305478e+04 7.343041e+03
[149,] 8.155946e+03 1.305478e+04
[150,] -2.482893e+04 8.155946e+03
[151,] -3.211069e+04 -2.482893e+04
[152,] 7.885890e+03 -3.211069e+04
[153,] -1.765617e+04 7.885890e+03
[154,] 7.415522e+03 -1.765617e+04
[155,] 7.885890e+03 7.415522e+03
[156,] 7.885890e+03 7.885890e+03
[157,] 6.074550e+03 7.885890e+03
[158,] 2.190668e+04 6.074550e+03
[159,] 7.885890e+03 2.190668e+04
[160,] 7.885890e+03 7.885890e+03
[161,] 7.885890e+03 7.885890e+03
[162,] 7.885890e+03 7.885890e+03
[163,] -9.604297e+03 7.885890e+03
[164,] -1.319946e+04 -9.604297e+03
[165,] 6.572742e+03 -1.319946e+04
[166,] 7.885890e+03 6.572742e+03
[167,] 6.282702e+04 7.885890e+03
[168,] -4.665849e+04 6.282702e+04
[169,] 7.885890e+03 -4.665849e+04
[170,] 7.885890e+03 7.885890e+03
[171,] 4.336193e+03 7.885890e+03
[172,] 7.885890e+03 4.336193e+03
[173,] -1.385396e+03 7.885890e+03
[174,] 7.885890e+03 -1.385396e+03
[175,] 5.876404e+03 7.885890e+03
[176,] -1.164544e+03 5.876404e+03
[177,] 7.451639e+03 -1.164544e+03
[178,] 7.885890e+03 7.451639e+03
[179,] 7.885890e+03 7.885890e+03
[180,] 3.412188e+04 7.885890e+03
[181,] 1.316928e+04 3.412188e+04
[182,] -3.448478e+04 1.316928e+04
[183,] -1.660040e+04 -3.448478e+04
[184,] 7.009061e+03 -1.660040e+04
[185,] 1.466212e+04 7.009061e+03
[186,] 3.633812e+03 1.466212e+04
[187,] 3.299849e+03 3.633812e+03
[188,] -1.950469e+04 3.299849e+03
[189,] 7.885890e+03 -1.950469e+04
[190,] 6.460031e+03 7.885890e+03
[191,] 7.885890e+03 6.460031e+03
[192,] 5.459244e+03 7.885890e+03
[193,] 7.885890e+03 5.459244e+03
[194,] 7.717053e+03 7.885890e+03
[195,] 7.885890e+03 7.717053e+03
[196,] -1.764143e+04 7.885890e+03
[197,] 7.885890e+03 -1.764143e+04
[198,] 7.885890e+03 7.885890e+03
[199,] -4.628102e+04 7.885890e+03
[200,] 7.885890e+03 -4.628102e+04
[201,] 7.885890e+03 7.885890e+03
[202,] 7.885890e+03 7.885890e+03
[203,] -1.032528e+04 7.885890e+03
[204,] 9.034898e+03 -1.032528e+04
[205,] 8.538142e+03 9.034898e+03
[206,] 1.021988e+04 8.538142e+03
[207,] 7.885890e+03 1.021988e+04
[208,] 1.765476e+04 7.885890e+03
[209,] -5.142160e+03 1.765476e+04
[210,] 7.885890e+03 -5.142160e+03
[211,] 7.885890e+03 7.885890e+03
[212,] 1.787679e+04 7.885890e+03
[213,] -2.674253e+04 1.787679e+04
[214,] -1.230533e+04 -2.674253e+04
[215,] 7.885890e+03 -1.230533e+04
[216,] 7.939671e+03 7.885890e+03
[217,] 3.356439e+03 7.939671e+03
[218,] 7.885890e+03 3.356439e+03
[219,] -1.661634e+04 7.885890e+03
[220,] -1.583636e+04 -1.661634e+04
[221,] 7.885890e+03 -1.583636e+04
[222,] 7.806647e+03 7.885890e+03
[223,] -3.449907e+03 7.806647e+03
[224,] 7.739504e+03 -3.449907e+03
[225,] 2.618073e+03 7.739504e+03
[226,] 5.998931e+03 2.618073e+03
[227,] 7.885890e+03 5.998931e+03
[228,] 1.647015e+03 7.885890e+03
[229,] 7.885890e+03 1.647015e+03
[230,] -1.425460e+04 7.885890e+03
[231,] 1.082325e+03 -1.425460e+04
[232,] 7.885890e+03 1.082325e+03
[233,] -1.064204e+04 7.885890e+03
[234,] 7.885890e+03 -1.064204e+04
[235,] 7.959474e+03 7.885890e+03
[236,] 7.885890e+03 7.959474e+03
[237,] 7.361795e+03 7.885890e+03
[238,] 7.885890e+03 7.361795e+03
[239,] 7.885890e+03 7.885890e+03
[240,] 7.885890e+03 7.885890e+03
[241,] 3.789819e+03 7.885890e+03
[242,] 7.885890e+03 3.789819e+03
[243,] 7.679544e+03 7.885890e+03
[244,] -2.674581e+04 7.679544e+03
[245,] 7.885890e+03 -2.674581e+04
[246,] 7.885890e+03 7.885890e+03
[247,] -7.364000e+02 7.885890e+03
[248,] 7.885890e+03 -7.364000e+02
[249,] 7.885890e+03 7.885890e+03
[250,] -1.457080e+04 7.885890e+03
[251,] 7.885890e+03 -1.457080e+04
[252,] -2.954105e+03 7.885890e+03
[253,] -1.427997e+04 -2.954105e+03
[254,] -1.001436e+05 -1.427997e+04
[255,] 7.885890e+03 -1.001436e+05
[256,] -4.294474e+03 7.885890e+03
[257,] -1.540906e+04 -4.294474e+03
[258,] 7.885890e+03 -1.540906e+04
[259,] 7.885890e+03 7.885890e+03
[260,] 7.885890e+03 7.885890e+03
[261,] 7.885890e+03 7.885890e+03
[262,] 3.280955e+04 7.885890e+03
[263,] 1.911200e+04 3.280955e+04
[264,] 7.885890e+03 1.911200e+04
[265,] 7.885890e+03 7.885890e+03
[266,] 7.885890e+03 7.885890e+03
[267,] 7.092566e+04 7.885890e+03
[268,] 7.885890e+03 7.092566e+04
[269,] 2.186414e+04 7.885890e+03
[270,] 7.885890e+03 2.186414e+04
[271,] 7.885890e+03 7.885890e+03
[272,] 7.885890e+03 7.885890e+03
[273,] 5.419260e+02 7.885890e+03
[274,] 7.885890e+03 5.419260e+02
[275,] 7.659666e+03 7.885890e+03
[276,] 7.125660e+03 7.659666e+03
[277,] 7.885890e+03 7.125660e+03
[278,] -3.947279e+04 7.885890e+03
[279,] 3.548361e+02 -3.947279e+04
[280,] 7.885890e+03 3.548361e+02
[281,] 7.885890e+03 7.885890e+03
[282,] 1.594011e+03 7.885890e+03
[283,] 7.885890e+03 1.594011e+03
[284,] 1.814935e+04 7.885890e+03
[285,] 7.885890e+03 1.814935e+04
[286,] 1.292396e+03 7.885890e+03
[287,] 1.015167e+04 1.292396e+03
[288,] 7.885890e+03 1.015167e+04
[289,] -1.919372e+04 7.885890e+03
[290,] 7.885890e+03 -1.919372e+04
[291,] 1.488961e+04 7.885890e+03
[292,] 5.046362e+04 1.488961e+04
[293,] 7.885890e+03 5.046362e+04
[294,] -6.062845e+03 7.885890e+03
[295,] 2.684495e+04 -6.062845e+03
[296,] 7.885890e+03 2.684495e+04
[297,] 7.885890e+03 7.885890e+03
[298,] 7.885890e+03 7.885890e+03
[299,] 7.885890e+03 7.885890e+03
[300,] 5.153591e+03 7.885890e+03
[301,] -1.138956e+04 5.153591e+03
[302,] 7.885890e+03 -1.138956e+04
[303,] 7.885890e+03 7.885890e+03
[304,] 3.058851e+04 7.885890e+03
[305,] 7.866032e+03 3.058851e+04
[306,] 7.885890e+03 7.866032e+03
[307,] 5.615897e+03 7.885890e+03
[308,] -4.332841e+03 5.615897e+03
[309,] 7.885890e+03 -4.332841e+03
[310,] 7.885890e+03 7.885890e+03
[311,] 3.292089e+04 7.885890e+03
[312,] 7.885890e+03 3.292089e+04
[313,] 1.385221e+04 7.885890e+03
[314,] 8.112119e+03 1.385221e+04
[315,] 2.983997e+04 8.112119e+03
[316,] 6.877481e+03 2.983997e+04
[317,] -1.673346e+04 6.877481e+03
[318,] 7.860491e+03 -1.673346e+04
[319,] 7.885890e+03 7.860491e+03
[320,] 7.885890e+03 7.885890e+03
[321,] 7.885890e+03 7.885890e+03
[322,] 7.885890e+03 7.885890e+03
[323,] 7.885890e+03 7.885890e+03
[324,] 7.885890e+03 7.885890e+03
[325,] 4.626436e+03 7.885890e+03
[326,] 7.885890e+03 4.626436e+03
[327,] 7.885890e+03 7.885890e+03
[328,] 7.885890e+03 7.885890e+03
[329,] 7.885890e+03 7.885890e+03
[330,] 7.885890e+03 7.885890e+03
[331,] 7.885890e+03 7.885890e+03
[332,] 7.885890e+03 7.885890e+03
[333,] 4.555494e+03 7.885890e+03
[334,] 7.876934e+03 4.555494e+03
[335,] 2.431005e+03 7.876934e+03
[336,] -2.600563e+04 2.431005e+03
[337,] 7.885890e+03 -2.600563e+04
[338,] 5.944099e+03 7.885890e+03
[339,] 7.998222e+03 5.944099e+03
[340,] -1.715618e+04 7.998222e+03
[341,] 7.885890e+03 -1.715618e+04
[342,] 7.885890e+03 7.885890e+03
[343,] 7.885890e+03 7.885890e+03
[344,] -1.000848e+04 7.885890e+03
[345,] 7.885890e+03 -1.000848e+04
[346,] -5.943610e+03 7.885890e+03
[347,] 7.885890e+03 -5.943610e+03
[348,] 7.885890e+03 7.885890e+03
[349,] -1.578813e+04 7.885890e+03
[350,] 5.758370e+03 -1.578813e+04
[351,] 5.533360e+03 5.758370e+03
[352,] 8.053020e+03 5.533360e+03
[353,] 4.424100e+03 8.053020e+03
[354,] 7.885890e+03 4.424100e+03
[355,] 7.885890e+03 7.885890e+03
[356,] 5.688983e+03 7.885890e+03
[357,] 7.885890e+03 5.688983e+03
[358,] 7.885890e+03 7.885890e+03
[359,] 7.885890e+03 7.885890e+03
[360,] 7.886657e+03 7.885890e+03
[361,] 7.885890e+03 7.886657e+03
[362,] 7.885890e+03 7.885890e+03
[363,] -2.990045e+04 7.885890e+03
[364,] 5.736787e+03 -2.990045e+04
[365,] 7.885890e+03 5.736787e+03
[366,] 7.885890e+03 7.885890e+03
[367,] 7.885890e+03 7.885890e+03
[368,] 7.885890e+03 7.885890e+03
[369,] 4.277989e+03 7.885890e+03
[370,] 7.885890e+03 4.277989e+03
[371,] 8.102971e+03 7.885890e+03
[372,] 2.559254e+04 8.102971e+03
[373,] -2.958819e+03 2.559254e+04
[374,] 7.885890e+03 -2.958819e+03
[375,] 7.885890e+03 7.885890e+03
[376,] -1.439895e+04 7.885890e+03
[377,] 4.326394e+03 -1.439895e+04
[378,] -2.599412e+04 4.326394e+03
[379,] 7.885890e+03 -2.599412e+04
[380,] 7.885890e+03 7.885890e+03
[381,] 3.793556e+02 7.885890e+03
[382,] -3.654900e+04 3.793556e+02
[383,] 7.885890e+03 -3.654900e+04
[384,] 8.571701e+03 7.885890e+03
[385,] -1.145814e+05 8.571701e+03
[386,] 1.726525e+04 -1.145814e+05
[387,] -1.431291e+04 1.726525e+04
[388,] 5.452475e+03 -1.431291e+04
[389,] -9.434267e+03 5.452475e+03
[390,] 7.948788e+03 -9.434267e+03
[391,] -2.315341e+04 7.948788e+03
[392,] -5.501282e+04 -2.315341e+04
[393,] 4.782787e+04 -5.501282e+04
[394,] -9.255679e+04 4.782787e+04
[395,] -3.360404e+04 -9.255679e+04
[396,] -1.441857e+05 -3.360404e+04
[397,] 9.729573e+04 -1.441857e+05
[398,] -2.046329e+04 9.729573e+04
[399,] -1.042779e+05 -2.046329e+04
[400,] -2.252725e+04 -1.042779e+05
[401,] 1.237066e+04 -2.252725e+04
[402,] -2.112778e+04 1.237066e+04
[403,] -5.631186e+04 -2.112778e+04
[404,] -4.095689e+04 -5.631186e+04
[405,] -2.063086e+03 -4.095689e+04
[406,] -1.439303e+03 -2.063086e+03
[407,] -1.994659e+04 -1.439303e+03
[408,] -1.698257e+04 -1.994659e+04
[409,] -8.863966e+04 -1.698257e+04
[410,] -2.180996e+04 -8.863966e+04
[411,] -1.516623e+05 -2.180996e+04
[412,] -5.757076e+04 -1.516623e+05
[413,] 7.750794e+04 -5.757076e+04
[414,] -1.793147e+05 7.750794e+04
[415,] -1.209000e+05 -1.793147e+05
[416,] 5.324414e+04 -1.209000e+05
[417,] -3.686172e+04 5.324414e+04
[418,] -5.522805e+04 -3.686172e+04
[419,] -7.260360e+04 -5.522805e+04
[420,] -7.235992e+04 -7.260360e+04
[421,] -1.199613e+05 -7.235992e+04
[422,] -4.749294e+05 -1.199613e+05
[423,] -4.586813e+05 -4.749294e+05
[424,] -7.261920e+04 -4.586813e+05
[425,] -1.848519e+05 -7.261920e+04
[426,] -3.293887e+05 -1.848519e+05
[427,] -1.089227e+06 -3.293887e+05
[428,] -2.374091e+05 -1.089227e+06
[429,] -2.333980e+05 -2.374091e+05
[430,] -5.464335e+05 -2.333980e+05
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3.411356e+06 -2.608799e+05
2 1.171853e+06 3.411356e+06
3 -5.428272e+05 1.171853e+06
4 -2.608628e+05 -5.428272e+05
5 1.933686e+04 -2.608628e+05
6 -3.509149e+03 1.933686e+04
7 -2.425525e+05 -3.509149e+03
8 2.364602e+04 -2.425525e+05
9 8.821366e+05 2.364602e+04
10 -3.321008e+05 8.821366e+05
11 4.790465e+05 -3.321008e+05
12 9.373054e+05 4.790465e+05
13 -6.683091e+05 9.373054e+05
14 7.445716e+04 -6.683091e+05
15 5.724179e+05 7.445716e+04
16 -3.540487e+03 5.724179e+05
17 9.051730e+03 -3.540487e+03
18 1.502868e+05 9.051730e+03
19 -3.216965e+05 1.502868e+05
20 2.868466e+05 -3.216965e+05
21 -4.064967e+04 2.868466e+05
22 -2.352861e+04 -4.064967e+04
23 -3.511651e+05 -2.352861e+04
24 -2.018015e+05 -3.511651e+05
25 -2.871585e+03 -2.018015e+05
26 -1.831363e+05 -2.871585e+03
27 -2.575756e+05 -1.831363e+05
28 3.628291e+05 -2.575756e+05
29 -9.667566e+04 3.628291e+05
30 3.120082e+05 -9.667566e+04
31 -9.784852e+03 3.120082e+05
32 2.947249e+05 -9.784852e+03
33 1.883540e+05 2.947249e+05
34 2.376275e+05 1.883540e+05
35 -3.559115e+05 2.376275e+05
36 1.154494e+05 -3.559115e+05
37 2.621138e+05 1.154494e+05
38 1.665263e+05 2.621138e+05
39 -1.803472e+05 1.665263e+05
40 -2.325873e+05 -1.803472e+05
41 -2.654010e+05 -2.325873e+05
42 9.691298e+04 -2.654010e+05
43 -9.742618e+04 9.691298e+04
44 -5.215692e+03 -9.742618e+04
45 -1.283432e+04 -5.215692e+03
46 -1.144812e+05 -1.283432e+04
47 -8.124632e+04 -1.144812e+05
48 3.854628e+04 -8.124632e+04
49 -1.566947e+05 3.854628e+04
50 1.503060e+05 -1.566947e+05
51 1.005303e+05 1.503060e+05
52 -6.448411e+04 1.005303e+05
53 1.912044e+05 -6.448411e+04
54 -2.109276e+05 1.912044e+05
55 5.695869e+04 -2.109276e+05
56 2.852243e+05 5.695869e+04
57 2.701651e+04 2.852243e+05
58 6.825481e+03 2.701651e+04
59 -1.145012e+05 6.825481e+03
60 6.538881e+04 -1.145012e+05
61 -1.874398e+05 6.538881e+04
62 4.005412e+04 -1.874398e+05
63 -9.580666e+04 4.005412e+04
64 1.824990e+03 -9.580666e+04
65 -2.330799e+05 1.824990e+03
66 -6.286751e+04 -2.330799e+05
67 4.498030e+04 -6.286751e+04
68 3.280319e+04 4.498030e+04
69 2.753540e+04 3.280319e+04
70 -1.460492e+05 2.753540e+04
71 -3.612186e+03 -1.460492e+05
72 -2.811354e+04 -3.612186e+03
73 8.244743e+03 -2.811354e+04
74 4.011274e+04 8.244743e+03
75 -9.423110e+04 4.011274e+04
76 3.847241e+04 -9.423110e+04
77 -8.553179e+04 3.847241e+04
78 -1.837766e+04 -8.553179e+04
79 1.870321e+04 -1.837766e+04
80 1.592131e+04 1.870321e+04
81 2.028738e+04 1.592131e+04
82 2.503243e+04 2.028738e+04
83 -1.852851e+04 2.503243e+04
84 6.045534e+04 -1.852851e+04
85 -5.147537e+04 6.045534e+04
86 7.445046e+04 -5.147537e+04
87 3.654182e+04 7.445046e+04
88 -4.602990e+04 3.654182e+04
89 6.791020e+04 -4.602990e+04
90 4.344371e+04 6.791020e+04
91 1.205713e+04 4.344371e+04
92 1.051645e+05 1.205713e+04
93 -1.120634e+05 1.051645e+05
94 2.316142e+04 -1.120634e+05
95 1.633573e+04 2.316142e+04
96 1.178379e+05 1.633573e+04
97 -2.796317e+03 1.178379e+05
98 3.084612e+04 -2.796317e+03
99 -1.528409e+04 3.084612e+04
100 -2.366499e+04 -1.528409e+04
101 -9.163387e+04 -2.366499e+04
102 -8.620064e+04 -9.163387e+04
103 -4.688865e+03 -8.620064e+04
104 3.003454e+04 -4.688865e+03
105 -5.128738e+04 3.003454e+04
106 -2.371219e+05 -5.128738e+04
107 1.594284e+05 -2.371219e+05
108 -1.143577e+05 1.594284e+05
109 -9.369199e+04 -1.143577e+05
110 8.462011e+04 -9.369199e+04
111 -2.032839e+05 8.462011e+04
112 -1.209303e+05 -2.032839e+05
113 -6.322223e+04 -1.209303e+05
114 4.690003e+04 -6.322223e+04
115 1.056765e+04 4.690003e+04
116 1.827586e+04 1.056765e+04
117 8.947695e+04 1.827586e+04
118 -1.361990e+05 8.947695e+04
119 -6.355304e+03 -1.361990e+05
120 -1.855450e+05 -6.355304e+03
121 1.119938e+04 -1.855450e+05
122 1.558382e+04 1.119938e+04
123 1.359708e+05 1.558382e+04
124 -1.027438e+03 1.359708e+05
125 -1.534732e+05 -1.027438e+03
126 7.885890e+03 -1.534732e+05
127 5.411699e+03 7.885890e+03
128 3.250702e+03 5.411699e+03
129 7.885890e+03 3.250702e+03
130 6.264633e+03 7.885890e+03
131 3.574508e+04 6.264633e+03
132 7.885890e+03 3.574508e+04
133 4.614398e+03 7.885890e+03
134 -6.605722e+03 4.614398e+03
135 2.031769e+03 -6.605722e+03
136 9.423296e+03 2.031769e+03
137 7.306172e+03 9.423296e+03
138 -8.138305e+03 7.306172e+03
139 7.885890e+03 -8.138305e+03
140 1.126219e+04 7.885890e+03
141 7.885890e+03 1.126219e+04
142 -2.597494e+04 7.885890e+03
143 7.724318e+03 -2.597494e+04
144 1.429069e+04 7.724318e+03
145 -2.892442e+01 1.429069e+04
146 7.885890e+03 -2.892442e+01
147 7.343041e+03 7.885890e+03
148 1.305478e+04 7.343041e+03
149 8.155946e+03 1.305478e+04
150 -2.482893e+04 8.155946e+03
151 -3.211069e+04 -2.482893e+04
152 7.885890e+03 -3.211069e+04
153 -1.765617e+04 7.885890e+03
154 7.415522e+03 -1.765617e+04
155 7.885890e+03 7.415522e+03
156 7.885890e+03 7.885890e+03
157 6.074550e+03 7.885890e+03
158 2.190668e+04 6.074550e+03
159 7.885890e+03 2.190668e+04
160 7.885890e+03 7.885890e+03
161 7.885890e+03 7.885890e+03
162 7.885890e+03 7.885890e+03
163 -9.604297e+03 7.885890e+03
164 -1.319946e+04 -9.604297e+03
165 6.572742e+03 -1.319946e+04
166 7.885890e+03 6.572742e+03
167 6.282702e+04 7.885890e+03
168 -4.665849e+04 6.282702e+04
169 7.885890e+03 -4.665849e+04
170 7.885890e+03 7.885890e+03
171 4.336193e+03 7.885890e+03
172 7.885890e+03 4.336193e+03
173 -1.385396e+03 7.885890e+03
174 7.885890e+03 -1.385396e+03
175 5.876404e+03 7.885890e+03
176 -1.164544e+03 5.876404e+03
177 7.451639e+03 -1.164544e+03
178 7.885890e+03 7.451639e+03
179 7.885890e+03 7.885890e+03
180 3.412188e+04 7.885890e+03
181 1.316928e+04 3.412188e+04
182 -3.448478e+04 1.316928e+04
183 -1.660040e+04 -3.448478e+04
184 7.009061e+03 -1.660040e+04
185 1.466212e+04 7.009061e+03
186 3.633812e+03 1.466212e+04
187 3.299849e+03 3.633812e+03
188 -1.950469e+04 3.299849e+03
189 7.885890e+03 -1.950469e+04
190 6.460031e+03 7.885890e+03
191 7.885890e+03 6.460031e+03
192 5.459244e+03 7.885890e+03
193 7.885890e+03 5.459244e+03
194 7.717053e+03 7.885890e+03
195 7.885890e+03 7.717053e+03
196 -1.764143e+04 7.885890e+03
197 7.885890e+03 -1.764143e+04
198 7.885890e+03 7.885890e+03
199 -4.628102e+04 7.885890e+03
200 7.885890e+03 -4.628102e+04
201 7.885890e+03 7.885890e+03
202 7.885890e+03 7.885890e+03
203 -1.032528e+04 7.885890e+03
204 9.034898e+03 -1.032528e+04
205 8.538142e+03 9.034898e+03
206 1.021988e+04 8.538142e+03
207 7.885890e+03 1.021988e+04
208 1.765476e+04 7.885890e+03
209 -5.142160e+03 1.765476e+04
210 7.885890e+03 -5.142160e+03
211 7.885890e+03 7.885890e+03
212 1.787679e+04 7.885890e+03
213 -2.674253e+04 1.787679e+04
214 -1.230533e+04 -2.674253e+04
215 7.885890e+03 -1.230533e+04
216 7.939671e+03 7.885890e+03
217 3.356439e+03 7.939671e+03
218 7.885890e+03 3.356439e+03
219 -1.661634e+04 7.885890e+03
220 -1.583636e+04 -1.661634e+04
221 7.885890e+03 -1.583636e+04
222 7.806647e+03 7.885890e+03
223 -3.449907e+03 7.806647e+03
224 7.739504e+03 -3.449907e+03
225 2.618073e+03 7.739504e+03
226 5.998931e+03 2.618073e+03
227 7.885890e+03 5.998931e+03
228 1.647015e+03 7.885890e+03
229 7.885890e+03 1.647015e+03
230 -1.425460e+04 7.885890e+03
231 1.082325e+03 -1.425460e+04
232 7.885890e+03 1.082325e+03
233 -1.064204e+04 7.885890e+03
234 7.885890e+03 -1.064204e+04
235 7.959474e+03 7.885890e+03
236 7.885890e+03 7.959474e+03
237 7.361795e+03 7.885890e+03
238 7.885890e+03 7.361795e+03
239 7.885890e+03 7.885890e+03
240 7.885890e+03 7.885890e+03
241 3.789819e+03 7.885890e+03
242 7.885890e+03 3.789819e+03
243 7.679544e+03 7.885890e+03
244 -2.674581e+04 7.679544e+03
245 7.885890e+03 -2.674581e+04
246 7.885890e+03 7.885890e+03
247 -7.364000e+02 7.885890e+03
248 7.885890e+03 -7.364000e+02
249 7.885890e+03 7.885890e+03
250 -1.457080e+04 7.885890e+03
251 7.885890e+03 -1.457080e+04
252 -2.954105e+03 7.885890e+03
253 -1.427997e+04 -2.954105e+03
254 -1.001436e+05 -1.427997e+04
255 7.885890e+03 -1.001436e+05
256 -4.294474e+03 7.885890e+03
257 -1.540906e+04 -4.294474e+03
258 7.885890e+03 -1.540906e+04
259 7.885890e+03 7.885890e+03
260 7.885890e+03 7.885890e+03
261 7.885890e+03 7.885890e+03
262 3.280955e+04 7.885890e+03
263 1.911200e+04 3.280955e+04
264 7.885890e+03 1.911200e+04
265 7.885890e+03 7.885890e+03
266 7.885890e+03 7.885890e+03
267 7.092566e+04 7.885890e+03
268 7.885890e+03 7.092566e+04
269 2.186414e+04 7.885890e+03
270 7.885890e+03 2.186414e+04
271 7.885890e+03 7.885890e+03
272 7.885890e+03 7.885890e+03
273 5.419260e+02 7.885890e+03
274 7.885890e+03 5.419260e+02
275 7.659666e+03 7.885890e+03
276 7.125660e+03 7.659666e+03
277 7.885890e+03 7.125660e+03
278 -3.947279e+04 7.885890e+03
279 3.548361e+02 -3.947279e+04
280 7.885890e+03 3.548361e+02
281 7.885890e+03 7.885890e+03
282 1.594011e+03 7.885890e+03
283 7.885890e+03 1.594011e+03
284 1.814935e+04 7.885890e+03
285 7.885890e+03 1.814935e+04
286 1.292396e+03 7.885890e+03
287 1.015167e+04 1.292396e+03
288 7.885890e+03 1.015167e+04
289 -1.919372e+04 7.885890e+03
290 7.885890e+03 -1.919372e+04
291 1.488961e+04 7.885890e+03
292 5.046362e+04 1.488961e+04
293 7.885890e+03 5.046362e+04
294 -6.062845e+03 7.885890e+03
295 2.684495e+04 -6.062845e+03
296 7.885890e+03 2.684495e+04
297 7.885890e+03 7.885890e+03
298 7.885890e+03 7.885890e+03
299 7.885890e+03 7.885890e+03
300 5.153591e+03 7.885890e+03
301 -1.138956e+04 5.153591e+03
302 7.885890e+03 -1.138956e+04
303 7.885890e+03 7.885890e+03
304 3.058851e+04 7.885890e+03
305 7.866032e+03 3.058851e+04
306 7.885890e+03 7.866032e+03
307 5.615897e+03 7.885890e+03
308 -4.332841e+03 5.615897e+03
309 7.885890e+03 -4.332841e+03
310 7.885890e+03 7.885890e+03
311 3.292089e+04 7.885890e+03
312 7.885890e+03 3.292089e+04
313 1.385221e+04 7.885890e+03
314 8.112119e+03 1.385221e+04
315 2.983997e+04 8.112119e+03
316 6.877481e+03 2.983997e+04
317 -1.673346e+04 6.877481e+03
318 7.860491e+03 -1.673346e+04
319 7.885890e+03 7.860491e+03
320 7.885890e+03 7.885890e+03
321 7.885890e+03 7.885890e+03
322 7.885890e+03 7.885890e+03
323 7.885890e+03 7.885890e+03
324 7.885890e+03 7.885890e+03
325 4.626436e+03 7.885890e+03
326 7.885890e+03 4.626436e+03
327 7.885890e+03 7.885890e+03
328 7.885890e+03 7.885890e+03
329 7.885890e+03 7.885890e+03
330 7.885890e+03 7.885890e+03
331 7.885890e+03 7.885890e+03
332 7.885890e+03 7.885890e+03
333 4.555494e+03 7.885890e+03
334 7.876934e+03 4.555494e+03
335 2.431005e+03 7.876934e+03
336 -2.600563e+04 2.431005e+03
337 7.885890e+03 -2.600563e+04
338 5.944099e+03 7.885890e+03
339 7.998222e+03 5.944099e+03
340 -1.715618e+04 7.998222e+03
341 7.885890e+03 -1.715618e+04
342 7.885890e+03 7.885890e+03
343 7.885890e+03 7.885890e+03
344 -1.000848e+04 7.885890e+03
345 7.885890e+03 -1.000848e+04
346 -5.943610e+03 7.885890e+03
347 7.885890e+03 -5.943610e+03
348 7.885890e+03 7.885890e+03
349 -1.578813e+04 7.885890e+03
350 5.758370e+03 -1.578813e+04
351 5.533360e+03 5.758370e+03
352 8.053020e+03 5.533360e+03
353 4.424100e+03 8.053020e+03
354 7.885890e+03 4.424100e+03
355 7.885890e+03 7.885890e+03
356 5.688983e+03 7.885890e+03
357 7.885890e+03 5.688983e+03
358 7.885890e+03 7.885890e+03
359 7.885890e+03 7.885890e+03
360 7.886657e+03 7.885890e+03
361 7.885890e+03 7.886657e+03
362 7.885890e+03 7.885890e+03
363 -2.990045e+04 7.885890e+03
364 5.736787e+03 -2.990045e+04
365 7.885890e+03 5.736787e+03
366 7.885890e+03 7.885890e+03
367 7.885890e+03 7.885890e+03
368 7.885890e+03 7.885890e+03
369 4.277989e+03 7.885890e+03
370 7.885890e+03 4.277989e+03
371 8.102971e+03 7.885890e+03
372 2.559254e+04 8.102971e+03
373 -2.958819e+03 2.559254e+04
374 7.885890e+03 -2.958819e+03
375 7.885890e+03 7.885890e+03
376 -1.439895e+04 7.885890e+03
377 4.326394e+03 -1.439895e+04
378 -2.599412e+04 4.326394e+03
379 7.885890e+03 -2.599412e+04
380 7.885890e+03 7.885890e+03
381 3.793556e+02 7.885890e+03
382 -3.654900e+04 3.793556e+02
383 7.885890e+03 -3.654900e+04
384 8.571701e+03 7.885890e+03
385 -1.145814e+05 8.571701e+03
386 1.726525e+04 -1.145814e+05
387 -1.431291e+04 1.726525e+04
388 5.452475e+03 -1.431291e+04
389 -9.434267e+03 5.452475e+03
390 7.948788e+03 -9.434267e+03
391 -2.315341e+04 7.948788e+03
392 -5.501282e+04 -2.315341e+04
393 4.782787e+04 -5.501282e+04
394 -9.255679e+04 4.782787e+04
395 -3.360404e+04 -9.255679e+04
396 -1.441857e+05 -3.360404e+04
397 9.729573e+04 -1.441857e+05
398 -2.046329e+04 9.729573e+04
399 -1.042779e+05 -2.046329e+04
400 -2.252725e+04 -1.042779e+05
401 1.237066e+04 -2.252725e+04
402 -2.112778e+04 1.237066e+04
403 -5.631186e+04 -2.112778e+04
404 -4.095689e+04 -5.631186e+04
405 -2.063086e+03 -4.095689e+04
406 -1.439303e+03 -2.063086e+03
407 -1.994659e+04 -1.439303e+03
408 -1.698257e+04 -1.994659e+04
409 -8.863966e+04 -1.698257e+04
410 -2.180996e+04 -8.863966e+04
411 -1.516623e+05 -2.180996e+04
412 -5.757076e+04 -1.516623e+05
413 7.750794e+04 -5.757076e+04
414 -1.793147e+05 7.750794e+04
415 -1.209000e+05 -1.793147e+05
416 5.324414e+04 -1.209000e+05
417 -3.686172e+04 5.324414e+04
418 -5.522805e+04 -3.686172e+04
419 -7.260360e+04 -5.522805e+04
420 -7.235992e+04 -7.260360e+04
421 -1.199613e+05 -7.235992e+04
422 -4.749294e+05 -1.199613e+05
423 -4.586813e+05 -4.749294e+05
424 -7.261920e+04 -4.586813e+05
425 -1.848519e+05 -7.261920e+04
426 -3.293887e+05 -1.848519e+05
427 -1.089227e+06 -3.293887e+05
428 -2.374091e+05 -1.089227e+06
429 -2.333980e+05 -2.374091e+05
430 -5.464335e+05 -2.333980e+05
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7g9ql1292955824.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/8g9ql1292955824.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/9g9ql1292955824.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/109ipo1292955824.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/11u05u1292955824.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/12yj401292955824.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/13ct281292955824.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/14421b1292955824.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/1583zh1292955824.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/164cx81292955824.tab")
+ }
>
> try(system("convert tmp/12hau1292955824.ps tmp/12hau1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/22hau1292955824.ps tmp/22hau1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/3v89f1292955824.ps tmp/3v89f1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/4v89f1292955824.ps tmp/4v89f1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/5v89f1292955824.ps tmp/5v89f1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/6o0801292955824.ps tmp/6o0801292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/7g9ql1292955824.ps tmp/7g9ql1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/8g9ql1292955824.ps tmp/8g9ql1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/9g9ql1292955824.ps tmp/9g9ql1292955824.png",intern=TRUE))
character(0)
> try(system("convert tmp/109ipo1292955824.ps tmp/109ipo1292955824.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.682 2.237 25.068