R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
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(807
+ ,213118
+ ,171986226
+ ,6282154
+ ,444
+ ,81767
+ ,36304548
+ ,4321023
+ ,412
+ ,153198
+ ,63117576
+ ,4111912
+ ,428
+ ,-26007
+ ,-11130996
+ ,223193
+ ,315
+ ,126942
+ ,39986730
+ ,1491348
+ ,168
+ ,157214
+ ,26411952
+ ,1629616
+ ,263
+ ,129352
+ ,34019576
+ ,1398893
+ ,267
+ ,234817
+ ,62696139
+ ,1926517
+ ,228
+ ,60448
+ ,13782144
+ ,983660
+ ,129
+ ,47818
+ ,6168522
+ ,1443586
+ ,104
+ ,245546
+ ,25536784
+ ,1073089
+ ,122
+ ,48020
+ ,5858440
+ ,984885
+ ,393
+ ,-1710
+ ,-672030
+ ,1405225
+ ,190
+ ,32648
+ ,6203120
+ ,227132
+ ,280
+ ,95350
+ ,26698000
+ ,929118
+ ,63
+ ,151352
+ ,9535176
+ ,1071292
+ ,102
+ ,288170
+ ,29393340
+ ,638830
+ ,265
+ ,114337
+ ,30299305
+ ,856956
+ ,234
+ ,37884
+ ,8864856
+ ,992426
+ ,277
+ ,122844
+ ,34027788
+ ,444477
+ ,73
+ ,82340
+ ,6010820
+ ,857217
+ ,67
+ ,79801
+ ,5346667
+ ,711969
+ ,103
+ ,165548
+ ,17051444
+ ,702380
+ ,290
+ ,116384
+ ,33751360
+ ,358589
+ ,83
+ ,134028
+ ,11124324
+ ,297978
+ ,56
+ ,63838
+ ,3574928
+ ,585715
+ ,236
+ ,74996
+ ,17699056
+ ,657954
+ ,73
+ ,31080
+ ,2268840
+ ,209458
+ ,34
+ ,32168
+ ,1093712
+ ,786690
+ ,139
+ ,49857
+ ,6930123
+ ,439798
+ ,26
+ ,87161
+ ,2266186
+ ,688779
+ ,70
+ ,106113
+ ,7427910
+ ,574339
+ ,40
+ ,80570
+ ,3222800
+ ,741409
+ ,42
+ ,102129
+ ,4289418
+ ,597793
+ ,12
+ ,301670
+ ,3620040
+ ,644190
+ ,211
+ ,102313
+ ,21588043
+ ,377934
+ ,74
+ ,88577
+ ,6554698
+ ,640273
+ ,80
+ ,112477
+ ,8998160
+ ,697458
+ ,83
+ ,191778
+ ,15917574
+ ,550608
+ ,131
+ ,79804
+ ,10454324
+ ,207393
+ ,203
+ ,128294
+ ,26043682
+ ,301607
+ ,56
+ ,96448
+ ,5401088
+ ,345783
+ ,89
+ ,93811
+ ,8349179
+ ,501749
+ ,88
+ ,117520
+ ,10341760
+ ,379983
+ ,39
+ ,69159
+ ,2697201
+ ,387475
+ ,25
+ ,101792
+ ,2544800
+ ,377305
+ ,49
+ ,210568
+ ,10317832
+ ,370837
+ ,149
+ ,136996
+ ,20412404
+ ,430866
+ ,58
+ ,121920
+ ,7071360
+ ,469107
+ ,41
+ ,76403
+ ,3132523
+ ,194493
+ ,90
+ ,108094
+ ,9728460
+ ,530670
+ ,136
+ ,134759
+ ,18327224
+ ,518365
+ ,97
+ ,188873
+ ,18320681
+ ,491303
+ ,63
+ ,146216
+ ,9211608
+ ,527021
+ ,114
+ ,156608
+ ,17853312
+ ,233773
+ ,77
+ ,61348
+ ,4723796
+ ,405972
+ ,6
+ ,50350
+ ,302100
+ ,652925
+ ,47
+ ,87720
+ ,4122840
+ ,446211
+ ,51
+ ,99489
+ ,5073939
+ ,341340
+ ,85
+ ,87419
+ ,7430615
+ ,387699
+ ,43
+ ,94355
+ ,4057265
+ ,493408
+ ,32
+ ,60326
+ ,1930432
+ ,146494
+ ,25
+ ,94670
+ ,2366750
+ ,414462
+ ,77
+ ,82425
+ ,6346725
+ ,364304
+ ,54
+ ,59017
+ ,3186918
+ ,355178
+ ,251
+ ,90829
+ ,22798079
+ ,357760
+ ,15
+ ,80791
+ ,1211865
+ ,261216
+ ,44
+ ,100423
+ ,4418612
+ ,397144
+ ,73
+ ,131116
+ ,9571468
+ ,374943
+ ,85
+ ,100269
+ ,8522865
+ ,424898
+ ,49
+ ,27330
+ ,1339170
+ ,202055
+ ,38
+ ,39039
+ ,1483482
+ ,378525
+ ,35
+ ,106885
+ ,3740975
+ ,310768
+ ,9
+ ,79285
+ ,713565
+ ,325738
+ ,34
+ ,118881
+ ,4041954
+ ,394510
+ ,20
+ ,77623
+ ,1552460
+ ,247060
+ ,29
+ ,114768
+ ,3328272
+ ,368078
+ ,11
+ ,74015
+ ,814165
+ ,236761
+ ,52
+ ,69465
+ ,3612180
+ ,312378
+ ,13
+ ,117869
+ ,1532297
+ ,339836
+ ,29
+ ,60982
+ ,1768478
+ ,347385
+ ,66
+ ,90131
+ ,5948646
+ ,426280
+ ,33
+ ,138971
+ ,4586043
+ ,352850
+ ,15
+ ,39625
+ ,594375
+ ,301881
+ ,15
+ ,102725
+ ,1540875
+ ,377516
+ ,68
+ ,64239
+ ,4368252
+ ,357312
+ ,100
+ ,90262
+ ,9026200
+ ,458343
+ ,13
+ ,103960
+ ,1351480
+ ,354228
+ ,45
+ ,106611
+ ,4797495
+ ,308636
+ ,14
+ ,103345
+ ,1446830
+ ,386212
+ ,36
+ ,95551
+ ,3439836
+ ,393343
+ ,40
+ ,82903
+ ,3316120
+ ,378509
+ ,68
+ ,63593
+ ,4324324
+ ,452469
+ ,29
+ ,126910
+ ,3680390
+ ,364839
+ ,43
+ ,37527
+ ,1613661
+ ,358649
+ ,30
+ ,60247
+ ,1807410
+ ,376641
+ ,9
+ ,112995
+ ,1016955
+ ,429112
+ ,22
+ ,70184
+ ,1544048
+ ,330546
+ ,19
+ ,130140
+ ,2472660
+ ,403560
+ ,9
+ ,73221
+ ,658989
+ ,317892
+ ,31
+ ,76114
+ ,2359534
+ ,307528
+ ,19
+ ,90534
+ ,1720146
+ ,235133
+ ,55
+ ,108479
+ ,5966345
+ ,299243
+ ,8
+ ,113761
+ ,910088
+ ,314073
+ ,28
+ ,68696
+ ,1923488
+ ,368186
+ ,29
+ ,71561
+ ,2075269
+ ,269661
+ ,48
+ ,59831
+ ,2871888
+ ,125390
+ ,16
+ ,97890
+ ,1566240
+ ,510834
+ ,47
+ ,101481
+ ,4769607
+ ,321896
+ ,20
+ ,72954
+ ,1459080
+ ,249898
+ ,22
+ ,67939
+ ,1494658
+ ,408881
+ ,33
+ ,48022
+ ,1584726
+ ,158492
+ ,44
+ ,86111
+ ,3788884
+ ,292154
+ ,13
+ ,74020
+ ,962260
+ ,289513
+ ,6
+ ,57530
+ ,345180
+ ,378049
+ ,35
+ ,56364
+ ,1972740
+ ,343466
+ ,8
+ ,84990
+ ,679920
+ ,332743
+ ,17
+ ,88590
+ ,1506030
+ ,442882
+ ,11
+ ,77200
+ ,849200
+ ,214215
+ ,21
+ ,61262
+ ,1286502
+ ,315688
+ ,92
+ ,110309
+ ,10148428
+ ,375195
+ ,12
+ ,67000
+ ,804000
+ ,334280
+ ,112
+ ,93099
+ ,10427088
+ ,355864
+ ,25
+ ,107577
+ ,2689425
+ ,480382
+ ,17
+ ,62920
+ ,1069640
+ ,353058
+ ,23
+ ,75832
+ ,1744136
+ ,217193
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,10
+ ,60793
+ ,607930
+ ,314533
+ ,23
+ ,57935
+ ,1332505
+ ,318056
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,60630
+ ,424410
+ ,314353
+ ,25
+ ,55637
+ ,1390925
+ ,369448
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,20
+ ,60887
+ ,1217740
+ ,312846
+ ,4
+ ,60720
+ ,242880
+ ,312075
+ ,4
+ ,60505
+ ,242020
+ ,315009
+ ,10
+ ,60945
+ ,609450
+ ,318903
+ ,1
+ ,60720
+ ,60720
+ ,314887
+ ,4
+ ,60720
+ ,242880
+ ,314913
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,8
+ ,58990
+ ,471920
+ ,325506
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,11
+ ,56750
+ ,624250
+ ,298568
+ ,4
+ ,60894
+ ,243576
+ ,315834
+ ,15
+ ,63346
+ ,950190
+ ,329784
+ ,9
+ ,56535
+ ,508815
+ ,312878
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,60835
+ ,425845
+ ,314987
+ ,2
+ ,60720
+ ,121440
+ ,325249
+ ,0
+ ,61016
+ ,0
+ ,315877
+ ,7
+ ,58650
+ ,410550
+ ,291650
+ ,46
+ ,60438
+ ,2780148
+ ,305959
+ ,5
+ ,60720
+ ,303600
+ ,315380
+ ,7
+ ,58625
+ ,410375
+ ,297765
+ ,2
+ ,60938
+ ,121876
+ ,315245
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,2
+ ,61490
+ ,122980
+ ,315236
+ ,5
+ ,60845
+ ,304225
+ ,336425
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,60830
+ ,425810
+ ,306268
+ ,24
+ ,63261
+ ,1518264
+ ,302187
+ ,1
+ ,60720
+ ,60720
+ ,314882
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,18
+ ,45689
+ ,822402
+ ,382712
+ ,55
+ ,60720
+ ,3339600
+ ,341570
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,3
+ ,61564
+ ,184692
+ ,312412
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,9
+ ,61938
+ ,557442
+ ,309596
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,8
+ ,60951
+ ,487608
+ ,315547
+ ,113
+ ,60720
+ ,6861360
+ ,313267
+ ,0
+ ,60745
+ ,0
+ ,316176
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,19
+ ,71642
+ ,1361198
+ ,359335
+ ,11
+ ,71641
+ ,788051
+ ,330068
+ ,25
+ ,55792
+ ,1394800
+ ,314289
+ ,16
+ ,71873
+ ,1149968
+ ,297413
+ ,5
+ ,62555
+ ,312775
+ ,314806
+ ,11
+ ,60370
+ ,664070
+ ,333210
+ ,23
+ ,64873
+ ,1492079
+ ,352108
+ ,6
+ ,62041
+ ,372246
+ ,313332
+ ,5
+ ,65745
+ ,328725
+ ,291787
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,59500
+ ,416500
+ ,318745
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,61630
+ ,431410
+ ,315366
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,3
+ ,60890
+ ,182670
+ ,315688
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,89
+ ,113521
+ ,10103369
+ ,409642
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,19
+ ,80045
+ ,1520855
+ ,269587
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,12
+ ,50804
+ ,609648
+ ,300962
+ ,12
+ ,87390
+ ,1048680
+ ,325479
+ ,5
+ ,61656
+ ,308280
+ ,316155
+ ,2
+ ,65688
+ ,131376
+ ,318574
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,26
+ ,48522
+ ,1261572
+ ,343613
+ ,3
+ ,60720
+ ,182160
+ ,306948
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,11
+ ,57640
+ ,634040
+ ,330059
+ ,10
+ ,61977
+ ,619770
+ ,288985
+ ,5
+ ,62620
+ ,313100
+ ,304485
+ ,2
+ ,60720
+ ,121440
+ ,315380
+ ,6
+ ,60831
+ ,364986
+ ,315688
+ ,7
+ ,60646
+ ,424522
+ ,317736
+ ,2
+ ,60720
+ ,121440
+ ,315380
+ ,28
+ ,56225
+ ,1574300
+ ,322331
+ ,3
+ ,60510
+ ,181530
+ ,296656
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,60698
+ ,60698
+ ,315354
+ ,20
+ ,60720
+ ,1214400
+ ,312161
+ ,1
+ ,60805
+ ,60805
+ ,315576
+ ,22
+ ,61404
+ ,1350888
+ ,314922
+ ,9
+ ,60720
+ ,546480
+ ,314551
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,2
+ ,65276
+ ,130552
+ ,312339
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,63915
+ ,447405
+ ,298700
+ ,9
+ ,60720
+ ,546480
+ ,321376
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,13
+ ,61686
+ ,801918
+ ,303230
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60743
+ ,0
+ ,315487
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,6
+ ,60349
+ ,362094
+ ,315793
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,3
+ ,61360
+ ,184080
+ ,312887
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,59818
+ ,418726
+ ,315637
+ ,2
+ ,72680
+ ,145360
+ ,324385
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,15
+ ,61808
+ ,927120
+ ,308989
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,9
+ ,53110
+ ,477990
+ ,296702
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,64245
+ ,64245
+ ,307322
+ ,38
+ ,73007
+ ,2774266
+ ,304376
+ ,57
+ ,82732
+ ,4715724
+ ,253588
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,54820
+ ,383740
+ ,309560
+ ,26
+ ,47705
+ ,1240330
+ ,298466
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,13
+ ,72835
+ ,946855
+ ,343929
+ ,10
+ ,58856
+ ,588560
+ ,331955
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,9
+ ,77655
+ ,698895
+ ,381180
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,26
+ ,69817
+ ,1815242
+ ,331420
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,19
+ ,60798
+ ,1155162
+ ,310201
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,12
+ ,62452
+ ,749424
+ ,320016
+ ,23
+ ,64175
+ ,1476025
+ ,320398
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,29
+ ,67440
+ ,1955760
+ ,291841
+ ,8
+ ,68136
+ ,545088
+ ,310670
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,26
+ ,56726
+ ,1474876
+ ,313491
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,9
+ ,70811
+ ,637299
+ ,331323
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,5
+ ,60720
+ ,303600
+ ,319210
+ ,3
+ ,62045
+ ,186135
+ ,318098
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,13
+ ,54323
+ ,706199
+ ,292754
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,12
+ ,62841
+ ,754092
+ ,325176
+ ,19
+ ,81125
+ ,1541375
+ ,365959
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,10
+ ,59506
+ ,595060
+ ,302409
+ ,9
+ ,59365
+ ,534285
+ ,340968
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,9
+ ,60720
+ ,546480
+ ,315380
+ ,4
+ ,60798
+ ,243192
+ ,313164
+ ,1
+ ,58790
+ ,58790
+ ,301164
+ ,1
+ ,60720
+ ,60720
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,14
+ ,61808
+ ,865312
+ ,344425
+ ,12
+ ,60735
+ ,728820
+ ,315394
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,19
+ ,64016
+ ,1216304
+ ,316647
+ ,17
+ ,54683
+ ,929611
+ ,309836
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,32
+ ,87192
+ ,2790144
+ ,346611
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,14
+ ,64107
+ ,897498
+ ,322031
+ ,8
+ ,60761
+ ,486088
+ ,315656
+ ,4
+ ,65990
+ ,263960
+ ,339445
+ ,0
+ ,59988
+ ,0
+ ,314964
+ ,20
+ ,61167
+ ,1223340
+ ,297141
+ ,5
+ ,60719
+ ,303595
+ ,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
+ ,1
+ ,60722
+ ,60722
+ ,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
+ ,4
+ ,60379
+ ,241516
+ ,313729
+ ,1
+ ,60727
+ ,60727
+ ,315388
+ ,4
+ ,60720
+ ,242880
+ ,315371
+ ,20
+ ,60925
+ ,1218500
+ ,296139
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,1
+ ,60896
+ ,60896
+ ,313880
+ ,10
+ ,59734
+ ,597340
+ ,317698
+ ,12
+ ,62969
+ ,755628
+ ,295580
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,13
+ ,60720
+ ,789360
+ ,308256
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,3
+ ,59118
+ ,177354
+ ,303677
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,10
+ ,60720
+ ,607200
+ ,319369
+ ,3
+ ,58598
+ ,175794
+ ,318690
+ ,7
+ ,61124
+ ,427868
+ ,314049
+ ,10
+ ,59595
+ ,595950
+ ,325699
+ ,1
+ ,62065
+ ,62065
+ ,314210
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,15
+ ,78780
+ ,1181700
+ ,322378
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,4
+ ,60722
+ ,242888
+ ,315398
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,28
+ ,61600
+ ,1724800
+ ,308336
+ ,9
+ ,59635
+ ,536715
+ ,316386
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,60720
+ ,425040
+ ,315553
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,7
+ ,59781
+ ,418467
+ ,323361
+ ,7
+ ,76644
+ ,536508
+ ,336639
+ ,3
+ ,64820
+ ,194460
+ ,307424
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,11
+ ,56178
+ ,617958
+ ,295370
+ ,7
+ ,60436
+ ,423052
+ ,322340
+ ,10
+ ,60720
+ ,607200
+ ,319864
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,18
+ ,73433
+ ,1321794
+ ,317291
+ ,14
+ ,41477
+ ,580678
+ ,280398
+ ,0
+ ,60720
+ ,0
+ ,315380
+ ,12
+ ,62700
+ ,752400
+ ,317330
+ ,29
+ ,67804
+ ,1966316
+ ,238125
+ ,3
+ ,59661
+ ,178983
+ ,327071
+ ,6
+ ,58620
+ ,351720
+ ,309038
+ ,3
+ ,60398
+ ,181194
+ ,314210
+ ,8
+ ,58580
+ ,468640
+ ,307930
+ ,10
+ ,62710
+ ,627100
+ ,322327
+ ,6
+ ,59325
+ ,355950
+ ,292136
+ ,8
+ ,60950
+ ,487600
+ ,263276
+ ,6
+ ,68060
+ ,408360
+ ,367655
+ ,9
+ ,83620
+ ,752580
+ ,283910
+ ,8
+ ,58456
+ ,467648
+ ,283587
+ ,26
+ ,52811
+ ,1373086
+ ,243650
+ ,239
+ ,121173
+ ,28960347
+ ,438493
+ ,7
+ ,63870
+ ,447090
+ ,296261
+ ,41
+ ,21001
+ ,861041
+ ,230621
+ ,3
+ ,70415
+ ,211245
+ ,304252
+ ,8
+ ,64230
+ ,513840
+ ,333505
+ ,6
+ ,59190
+ ,355140
+ ,296919
+ ,21
+ ,69351
+ ,1456371
+ ,278990
+ ,7
+ ,64270
+ ,449890
+ ,276898
+ ,11
+ ,70694
+ ,777634
+ ,327007
+ ,11
+ ,68005
+ ,748055
+ ,317046
+ ,12
+ ,58930
+ ,707160
+ ,304555
+ ,9
+ ,58320
+ ,524880
+ ,298096
+ ,3
+ ,69980
+ ,209940
+ ,231861
+ ,57
+ ,69863
+ ,3982191
+ ,309422
+ ,21
+ ,63255
+ ,1328355
+ ,286963
+ ,15
+ ,57320
+ ,859800
+ ,269753
+ ,32
+ ,75230
+ ,2407360
+ ,448243
+ ,11
+ ,79420
+ ,873620
+ ,165404
+ ,2
+ ,73490
+ ,146980
+ ,204325
+ ,23
+ ,35250
+ ,810750
+ ,407159
+ ,20
+ ,62285
+ ,1245700
+ ,290476
+ ,24
+ ,69206
+ ,1660944
+ ,275311
+ ,1
+ ,65920
+ ,65920
+ ,246541
+ ,1
+ ,69770
+ ,69770
+ ,253468
+ ,74
+ ,72683
+ ,5378542
+ ,240897
+ ,68
+ ,-14545
+ ,-989060
+ ,-83265
+ ,20
+ ,55830
+ ,1116600
+ ,-42143
+ ,20
+ ,55174
+ ,1103480
+ ,272713
+ ,82
+ ,67038
+ ,5497116
+ ,215362
+ ,21
+ ,51252
+ ,1076292
+ ,42754
+ ,244
+ ,157278
+ ,38375832
+ ,306275
+ ,32
+ ,79510
+ ,2544320
+ ,253537
+ ,86
+ ,77440
+ ,6659840
+ ,372631
+ ,69
+ ,27284
+ ,1882596
+ ,-7170)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('Orders'
+ ,'Dividends'
+ ,'O_D'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Orders','Dividends','O_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
> 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\r Orders Dividends O_D
1 6282154 807 213118 171986226
2 4321023 444 81767 36304548
3 4111912 412 153198 63117576
4 223193 428 -26007 -11130996
5 1491348 315 126942 39986730
6 1629616 168 157214 26411952
7 1398893 263 129352 34019576
8 1926517 267 234817 62696139
9 983660 228 60448 13782144
10 1443586 129 47818 6168522
11 1073089 104 245546 25536784
12 984885 122 48020 5858440
13 1405225 393 -1710 -672030
14 227132 190 32648 6203120
15 929118 280 95350 26698000
16 1071292 63 151352 9535176
17 638830 102 288170 29393340
18 856956 265 114337 30299305
19 992426 234 37884 8864856
20 444477 277 122844 34027788
21 857217 73 82340 6010820
22 711969 67 79801 5346667
23 702380 103 165548 17051444
24 358589 290 116384 33751360
25 297978 83 134028 11124324
26 585715 56 63838 3574928
27 657954 236 74996 17699056
28 209458 73 31080 2268840
29 786690 34 32168 1093712
30 439798 139 49857 6930123
31 688779 26 87161 2266186
32 574339 70 106113 7427910
33 741409 40 80570 3222800
34 597793 42 102129 4289418
35 644190 12 301670 3620040
36 377934 211 102313 21588043
37 640273 74 88577 6554698
38 697458 80 112477 8998160
39 550608 83 191778 15917574
40 207393 131 79804 10454324
41 301607 203 128294 26043682
42 345783 56 96448 5401088
43 501749 89 93811 8349179
44 379983 88 117520 10341760
45 387475 39 69159 2697201
46 377305 25 101792 2544800
47 370837 49 210568 10317832
48 430866 149 136996 20412404
49 469107 58 121920 7071360
50 194493 41 76403 3132523
51 530670 90 108094 9728460
52 518365 136 134759 18327224
53 491303 97 188873 18320681
54 527021 63 146216 9211608
55 233773 114 156608 17853312
56 405972 77 61348 4723796
57 652925 6 50350 302100
58 446211 47 87720 4122840
59 341340 51 99489 5073939
60 387699 85 87419 7430615
61 493408 43 94355 4057265
62 146494 32 60326 1930432
63 414462 25 94670 2366750
64 364304 77 82425 6346725
65 355178 54 59017 3186918
66 357760 251 90829 22798079
67 261216 15 80791 1211865
68 397144 44 100423 4418612
69 374943 73 131116 9571468
70 424898 85 100269 8522865
71 202055 49 27330 1339170
72 378525 38 39039 1483482
73 310768 35 106885 3740975
74 325738 9 79285 713565
75 394510 34 118881 4041954
76 247060 20 77623 1552460
77 368078 29 114768 3328272
78 236761 11 74015 814165
79 312378 52 69465 3612180
80 339836 13 117869 1532297
81 347385 29 60982 1768478
82 426280 66 90131 5948646
83 352850 33 138971 4586043
84 301881 15 39625 594375
85 377516 15 102725 1540875
86 357312 68 64239 4368252
87 458343 100 90262 9026200
88 354228 13 103960 1351480
89 308636 45 106611 4797495
90 386212 14 103345 1446830
91 393343 36 95551 3439836
92 378509 40 82903 3316120
93 452469 68 63593 4324324
94 364839 29 126910 3680390
95 358649 43 37527 1613661
96 376641 30 60247 1807410
97 429112 9 112995 1016955
98 330546 22 70184 1544048
99 403560 19 130140 2472660
100 317892 9 73221 658989
101 307528 31 76114 2359534
102 235133 19 90534 1720146
103 299243 55 108479 5966345
104 314073 8 113761 910088
105 368186 28 68696 1923488
106 269661 29 71561 2075269
107 125390 48 59831 2871888
108 510834 16 97890 1566240
109 321896 47 101481 4769607
110 249898 20 72954 1459080
111 408881 22 67939 1494658
112 158492 33 48022 1584726
113 292154 44 86111 3788884
114 289513 13 74020 962260
115 378049 6 57530 345180
116 343466 35 56364 1972740
117 332743 8 84990 679920
118 442882 17 88590 1506030
119 214215 11 77200 849200
120 315688 21 61262 1286502
121 375195 92 110309 10148428
122 334280 12 67000 804000
123 355864 112 93099 10427088
124 480382 25 107577 2689425
125 353058 17 62920 1069640
126 217193 23 75832 1744136
127 315380 0 60720 0
128 314533 10 60793 607930
129 318056 23 57935 1332505
130 315380 0 60720 0
131 314353 7 60630 424410
132 369448 25 55637 1390925
133 315380 0 60720 0
134 312846 20 60887 1217740
135 312075 4 60720 242880
136 315009 4 60505 242020
137 318903 10 60945 609450
138 314887 1 60720 60720
139 314913 4 60720 242880
140 315380 0 60720 0
141 325506 8 58990 471920
142 315380 0 60720 0
143 298568 11 56750 624250
144 315834 4 60894 243576
145 329784 15 63346 950190
146 312878 9 56535 508815
147 315380 0 60720 0
148 314987 7 60835 425845
149 325249 2 60720 121440
150 315877 0 61016 0
151 291650 7 58650 410550
152 305959 46 60438 2780148
153 315380 5 60720 303600
154 297765 7 58625 410375
155 315245 2 60938 121876
156 315380 0 60720 0
157 315380 0 60720 0
158 315236 2 61490 122980
159 336425 5 60845 304225
160 315380 0 60720 0
161 315380 0 60720 0
162 315380 0 60720 0
163 315380 0 60720 0
164 306268 7 60830 425810
165 302187 24 63261 1518264
166 314882 1 60720 60720
167 315380 0 60720 0
168 382712 18 45689 822402
169 341570 55 60720 3339600
170 315380 0 60720 0
171 315380 0 60720 0
172 312412 3 61564 184692
173 315380 0 60720 0
174 309596 9 61938 557442
175 315380 0 60720 0
176 315547 8 60951 487608
177 313267 113 60720 6861360
178 316176 0 60745 0
179 315380 0 60720 0
180 315380 0 60720 0
181 359335 19 71642 1361198
182 330068 11 71641 788051
183 314289 25 55792 1394800
184 297413 16 71873 1149968
185 314806 5 62555 312775
186 333210 11 60370 664070
187 352108 23 64873 1492079
188 313332 6 62041 372246
189 291787 5 65745 328725
190 315380 0 60720 0
191 318745 7 59500 416500
192 315380 0 60720 0
193 315366 7 61630 431410
194 315380 0 60720 0
195 315688 3 60890 182670
196 315380 0 60720 0
197 409642 89 113521 10103369
198 315380 0 60720 0
199 315380 0 60720 0
200 269587 19 80045 1520855
201 315380 0 60720 0
202 315380 0 60720 0
203 315380 0 60720 0
204 300962 12 50804 609648
205 325479 12 87390 1048680
206 316155 5 61656 308280
207 318574 2 65688 131376
208 315380 0 60720 0
209 343613 26 48522 1261572
210 306948 3 60720 182160
211 315380 0 60720 0
212 315380 0 60720 0
213 330059 11 57640 634040
214 288985 10 61977 619770
215 304485 5 62620 313100
216 315380 2 60720 121440
217 315688 6 60831 364986
218 317736 7 60646 424522
219 315380 2 60720 121440
220 322331 28 56225 1574300
221 296656 3 60510 181530
222 315380 0 60720 0
223 315354 1 60698 60698
224 312161 20 60720 1214400
225 315576 1 60805 60805
226 314922 22 61404 1350888
227 314551 9 60720 546480
228 315380 0 60720 0
229 312339 2 65276 130552
230 315380 0 60720 0
231 298700 7 63915 447405
232 321376 9 60720 546480
233 315380 0 60720 0
234 303230 13 61686 801918
235 315380 0 60720 0
236 315487 0 60743 0
237 315380 0 60720 0
238 315793 6 60349 362094
239 315380 0 60720 0
240 315380 0 60720 0
241 315380 0 60720 0
242 312887 3 61360 184080
243 315380 0 60720 0
244 315637 7 59818 418726
245 324385 2 72680 145360
246 315380 0 60720 0
247 315380 0 60720 0
248 308989 15 61808 927120
249 315380 0 60720 0
250 315380 0 60720 0
251 296702 9 53110 477990
252 315380 0 60720 0
253 307322 1 64245 64245
254 304376 38 73007 2774266
255 253588 57 82732 4715724
256 315380 0 60720 0
257 309560 7 54820 383740
258 298466 26 47705 1240330
259 315380 0 60720 0
260 315380 0 60720 0
261 315380 0 60720 0
262 315380 0 60720 0
263 343929 13 72835 946855
264 331955 10 58856 588560
265 315380 0 60720 0
266 315380 0 60720 0
267 315380 0 60720 0
268 381180 9 77655 698895
269 315380 0 60720 0
270 331420 26 69817 1815242
271 315380 0 60720 0
272 315380 0 60720 0
273 315380 0 60720 0
274 310201 19 60798 1155162
275 315380 0 60720 0
276 320016 12 62452 749424
277 320398 23 64175 1476025
278 315380 0 60720 0
279 291841 29 67440 1955760
280 310670 8 68136 545088
281 315380 0 60720 0
282 315380 0 60720 0
283 313491 26 56726 1474876
284 315380 0 60720 0
285 331323 9 70811 637299
286 315380 0 60720 0
287 319210 5 60720 303600
288 318098 3 62045 186135
289 315380 0 60720 0
290 292754 13 54323 706199
291 315380 0 60720 0
292 325176 12 62841 754092
293 365959 19 81125 1541375
294 315380 0 60720 0
295 302409 10 59506 595060
296 340968 9 59365 534285
297 315380 0 60720 0
298 315380 0 60720 0
299 315380 0 60720 0
300 315380 9 60720 546480
301 313164 4 60798 243192
302 301164 1 58790 58790
303 315380 1 60720 60720
304 315380 0 60720 0
305 344425 14 61808 865312
306 315394 12 60735 728820
307 315380 0 60720 0
308 316647 19 64016 1216304
309 309836 17 54683 929611
310 315380 0 60720 0
311 315380 0 60720 0
312 346611 32 87192 2790144
313 315380 0 60720 0
314 322031 14 64107 897498
315 315656 8 60761 486088
316 339445 4 65990 263960
317 314964 0 59988 0
318 297141 20 61167 1223340
319 315372 5 60719 303595
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 1 60722 60722
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 4 60379 241516
335 315388 1 60727 60727
336 315371 4 60720 242880
337 296139 20 60925 1218500
338 315380 0 60720 0
339 313880 1 60896 60896
340 317698 10 59734 597340
341 295580 12 62969 755628
342 315380 0 60720 0
343 315380 0 60720 0
344 315380 0 60720 0
345 308256 13 60720 789360
346 315380 0 60720 0
347 303677 3 59118 177354
348 315380 0 60720 0
349 315380 0 60720 0
350 319369 10 60720 607200
351 318690 3 58598 175794
352 314049 7 61124 427868
353 325699 10 59595 595950
354 314210 1 62065 62065
355 315380 0 60720 0
356 315380 0 60720 0
357 322378 15 78780 1181700
358 315380 0 60720 0
359 315380 0 60720 0
360 315380 0 60720 0
361 315398 4 60722 242888
362 315380 0 60720 0
363 315380 0 60720 0
364 308336 28 61600 1724800
365 316386 9 59635 536715
366 315380 0 60720 0
367 315380 0 60720 0
368 315380 0 60720 0
369 315380 0 60720 0
370 315553 7 60720 425040
371 315380 0 60720 0
372 323361 7 59781 418467
373 336639 7 76644 536508
374 307424 3 64820 194460
375 315380 0 60720 0
376 315380 0 60720 0
377 295370 11 56178 617958
378 322340 7 60436 423052
379 319864 10 60720 607200
380 315380 0 60720 0
381 315380 0 60720 0
382 317291 18 73433 1321794
383 280398 14 41477 580678
384 315380 0 60720 0
385 317330 12 62700 752400
386 238125 29 67804 1966316
387 327071 3 59661 178983
388 309038 6 58620 351720
389 314210 3 60398 181194
390 307930 8 58580 468640
391 322327 10 62710 627100
392 292136 6 59325 355950
393 263276 8 60950 487600
394 367655 6 68060 408360
395 283910 9 83620 752580
396 283587 8 58456 467648
397 243650 26 52811 1373086
398 438493 239 121173 28960347
399 296261 7 63870 447090
400 230621 41 21001 861041
401 304252 3 70415 211245
402 333505 8 64230 513840
403 296919 6 59190 355140
404 278990 21 69351 1456371
405 276898 7 64270 449890
406 327007 11 70694 777634
407 317046 11 68005 748055
408 304555 12 58930 707160
409 298096 9 58320 524880
410 231861 3 69980 209940
411 309422 57 69863 3982191
412 286963 21 63255 1328355
413 269753 15 57320 859800
414 448243 32 75230 2407360
415 165404 11 79420 873620
416 204325 2 73490 146980
417 407159 23 35250 810750
418 290476 20 62285 1245700
419 275311 24 69206 1660944
420 246541 1 65920 65920
421 253468 1 69770 69770
422 240897 74 72683 5378542
423 -83265 68 -14545 -989060
424 -42143 20 55830 1116600
425 272713 20 55174 1103480
426 215362 82 67038 5497116
427 42754 21 51252 1076292
428 306275 244 157278 38375832
429 253537 32 79510 2544320
430 372631 86 77440 6659840
431 -7170 69 27284 1882596
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Orders Dividends O_D
3.917e+05 8.432e+02 -1.823e+00 3.057e-02
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-1177595 -23719 23639 34415 2594168
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.917e+05 3.341e+04 11.725 < 2e-16 ***
Orders 8.432e+02 2.857e+02 2.951 0.00334 **
Dividends -1.823e+00 4.490e-01 -4.061 5.81e-05 ***
O_D 3.057e-02 2.202e-03 13.883 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 234800 on 427 degrees of freedom
Multiple R-squared: 0.7047, Adjusted R-squared: 0.7026
F-statistic: 339.7 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 2.412205e-13 1.206102e-13
[2,] 1.0000000 1.766274e-72 8.831368e-73
[3,] 1.0000000 4.413854e-73 2.206927e-73
[4,] 1.0000000 2.307902e-87 1.153951e-87
[5,] 1.0000000 8.967041e-91 4.483520e-91
[6,] 1.0000000 1.097964e-93 5.489818e-94
[7,] 1.0000000 6.203572e-95 3.101786e-95
[8,] 1.0000000 2.514953e-100 1.257476e-100
[9,] 1.0000000 4.047570e-106 2.023785e-106
[10,] 1.0000000 9.648272e-114 4.824136e-114
[11,] 1.0000000 3.830036e-115 1.915018e-115
[12,] 1.0000000 2.375063e-122 1.187532e-122
[13,] 1.0000000 3.431091e-127 1.715546e-127
[14,] 1.0000000 5.530630e-137 2.765315e-137
[15,] 1.0000000 2.300575e-143 1.150288e-143
[16,] 1.0000000 2.623895e-146 1.311948e-146
[17,] 1.0000000 2.403232e-147 1.201616e-147
[18,] 1.0000000 4.039002e-157 2.019501e-157
[19,] 1.0000000 4.866210e-157 2.433105e-157
[20,] 1.0000000 3.271720e-158 1.635860e-158
[21,] 1.0000000 2.096207e-161 1.048103e-161
[22,] 1.0000000 5.982313e-161 2.991157e-161
[23,] 1.0000000 5.428947e-171 2.714473e-171
[24,] 1.0000000 3.620383e-171 1.810192e-171
[25,] 1.0000000 8.639791e-176 4.319896e-176
[26,] 1.0000000 1.617980e-176 8.089899e-177
[27,] 1.0000000 5.849597e-184 2.924798e-184
[28,] 1.0000000 8.062963e-186 4.031482e-186
[29,] 1.0000000 1.128703e-186 5.643516e-187
[30,] 1.0000000 3.965101e-189 1.982550e-189
[31,] 1.0000000 5.263843e-193 2.631921e-193
[32,] 1.0000000 3.451767e-198 1.725883e-198
[33,] 1.0000000 1.813116e-197 9.065582e-198
[34,] 1.0000000 1.141767e-198 5.708837e-199
[35,] 1.0000000 2.206972e-202 1.103486e-202
[36,] 1.0000000 2.400054e-201 1.200027e-201
[37,] 1.0000000 7.448426e-202 3.724213e-202
[38,] 1.0000000 5.198009e-201 2.599005e-201
[39,] 1.0000000 2.803767e-200 1.401883e-200
[40,] 1.0000000 2.917671e-199 1.458836e-199
[41,] 1.0000000 9.883787e-200 4.941893e-200
[42,] 1.0000000 5.845799e-200 2.922900e-200
[43,] 1.0000000 2.838734e-199 1.419367e-199
[44,] 1.0000000 1.771905e-199 8.859527e-200
[45,] 1.0000000 3.040286e-200 1.520143e-200
[46,] 1.0000000 1.067648e-200 5.338238e-201
[47,] 1.0000000 6.883826e-200 3.441913e-200
[48,] 1.0000000 1.191128e-199 5.955641e-200
[49,] 1.0000000 8.534312e-203 4.267156e-203
[50,] 1.0000000 9.299268e-203 4.649634e-203
[51,] 1.0000000 1.796220e-212 8.981098e-213
[52,] 1.0000000 2.208819e-212 1.104409e-212
[53,] 1.0000000 2.421322e-211 1.210661e-211
[54,] 1.0000000 1.167255e-210 5.836275e-211
[55,] 1.0000000 1.829768e-211 9.148838e-212
[56,] 1.0000000 1.440850e-212 7.204252e-213
[57,] 1.0000000 6.281024e-212 3.140512e-212
[58,] 1.0000000 4.107223e-211 2.053611e-211
[59,] 1.0000000 1.826123e-210 9.130615e-211
[60,] 1.0000000 4.783515e-213 2.391758e-213
[61,] 1.0000000 1.836706e-212 9.183530e-213
[62,] 1.0000000 1.414260e-211 7.071301e-212
[63,] 1.0000000 1.244691e-210 6.223456e-211
[64,] 1.0000000 3.003511e-210 1.501756e-210
[65,] 1.0000000 2.108897e-209 1.054449e-209
[66,] 1.0000000 1.150879e-209 5.754396e-210
[67,] 1.0000000 6.600556e-209 3.300278e-209
[68,] 1.0000000 6.424888e-208 3.212444e-208
[69,] 1.0000000 6.369586e-207 3.184793e-207
[70,] 1.0000000 1.767582e-206 8.837912e-207
[71,] 1.0000000 1.688796e-205 8.443982e-206
[72,] 1.0000000 3.262399e-205 1.631200e-205
[73,] 1.0000000 3.246616e-204 1.623308e-204
[74,] 1.0000000 1.513042e-203 7.565208e-204
[75,] 1.0000000 9.343020e-203 4.671510e-203
[76,] 1.0000000 1.353297e-202 6.766486e-203
[77,] 1.0000000 4.143209e-202 2.071605e-202
[78,] 1.0000000 3.422896e-201 1.711448e-201
[79,] 1.0000000 2.891949e-200 1.445975e-200
[80,] 1.0000000 7.406198e-200 3.703099e-200
[81,] 1.0000000 5.486298e-201 2.743149e-201
[82,] 1.0000000 4.335011e-200 2.167505e-200
[83,] 1.0000000 2.029235e-199 1.014617e-199
[84,] 1.0000000 1.610871e-198 8.054354e-199
[85,] 1.0000000 1.070382e-197 5.351908e-198
[86,] 1.0000000 5.763223e-197 2.881612e-197
[87,] 1.0000000 2.397773e-199 1.198886e-199
[88,] 1.0000000 1.464572e-198 7.322861e-199
[89,] 1.0000000 3.359115e-199 1.679558e-199
[90,] 1.0000000 4.821630e-199 2.410815e-199
[91,] 1.0000000 2.432947e-198 1.216473e-198
[92,] 1.0000000 2.242750e-197 1.121375e-197
[93,] 1.0000000 1.623685e-196 8.118426e-197
[94,] 1.0000000 1.529186e-195 7.645930e-196
[95,] 1.0000000 1.501032e-194 7.505160e-195
[96,] 1.0000000 5.254379e-195 2.627190e-195
[97,] 1.0000000 1.278745e-194 6.393724e-195
[98,] 1.0000000 8.344600e-195 4.172300e-195
[99,] 1.0000000 3.305689e-194 1.652844e-194
[100,] 1.0000000 2.015257e-193 1.007629e-193
[101,] 1.0000000 5.848798e-195 2.924399e-195
[102,] 1.0000000 4.088605e-196 2.044303e-196
[103,] 1.0000000 2.387243e-195 1.193622e-195
[104,] 1.0000000 5.714812e-195 2.857406e-195
[105,] 1.0000000 4.196784e-195 2.098392e-195
[106,] 1.0000000 2.103089e-195 1.051545e-195
[107,] 1.0000000 1.220652e-194 6.103260e-195
[108,] 1.0000000 7.554101e-194 3.777050e-194
[109,] 1.0000000 1.408838e-193 7.044192e-194
[110,] 1.0000000 4.118940e-193 2.059470e-193
[111,] 1.0000000 3.141841e-192 1.570920e-192
[112,] 1.0000000 3.576496e-192 1.788248e-192
[113,] 1.0000000 4.391048e-193 2.195524e-193
[114,] 1.0000000 4.032869e-192 2.016435e-192
[115,] 1.0000000 2.771426e-191 1.385713e-191
[116,] 1.0000000 2.451994e-190 1.225997e-190
[117,] 1.0000000 6.462261e-190 3.231130e-190
[118,] 1.0000000 4.241149e-190 2.120575e-190
[119,] 1.0000000 2.001628e-189 1.000814e-189
[120,] 1.0000000 7.152407e-190 3.576203e-190
[121,] 1.0000000 6.761905e-189 3.380952e-189
[122,] 1.0000000 6.418398e-188 3.209199e-188
[123,] 1.0000000 4.910069e-187 2.455034e-187
[124,] 1.0000000 4.584165e-186 2.292083e-186
[125,] 1.0000000 4.307217e-185 2.153609e-185
[126,] 1.0000000 4.452139e-185 2.226069e-185
[127,] 1.0000000 4.136262e-184 2.068131e-184
[128,] 1.0000000 3.686498e-183 1.843249e-183
[129,] 1.0000000 3.444550e-182 1.722275e-182
[130,] 1.0000000 3.175257e-181 1.587628e-181
[131,] 1.0000000 2.815349e-180 1.407675e-180
[132,] 1.0000000 2.566842e-179 1.283421e-179
[133,] 1.0000000 2.337149e-178 1.168575e-178
[134,] 1.0000000 2.106308e-177 1.053154e-177
[135,] 1.0000000 1.667822e-176 8.339112e-177
[136,] 1.0000000 1.490812e-175 7.454062e-176
[137,] 1.0000000 1.365502e-174 6.827509e-175
[138,] 1.0000000 1.211944e-173 6.059719e-174
[139,] 1.0000000 9.261667e-173 4.630834e-173
[140,] 1.0000000 7.789612e-172 3.894806e-172
[141,] 1.0000000 6.806988e-171 3.403494e-171
[142,] 1.0000000 5.928339e-170 2.964170e-170
[143,] 1.0000000 4.866655e-169 2.433328e-169
[144,] 1.0000000 4.194526e-168 2.097263e-168
[145,] 1.0000000 3.561164e-167 1.780582e-167
[146,] 1.0000000 2.158148e-166 1.079074e-166
[147,] 1.0000000 1.837582e-165 9.187910e-166
[148,] 1.0000000 1.581189e-164 7.905943e-165
[149,] 1.0000000 1.336549e-163 6.682747e-164
[150,] 1.0000000 1.121648e-162 5.608239e-163
[151,] 1.0000000 9.372090e-162 4.686045e-162
[152,] 1.0000000 7.824196e-161 3.912098e-161
[153,] 1.0000000 5.288853e-160 2.644427e-160
[154,] 1.0000000 4.364169e-159 2.182085e-159
[155,] 1.0000000 3.585304e-158 1.792652e-158
[156,] 1.0000000 2.932445e-157 1.466223e-157
[157,] 1.0000000 2.387858e-156 1.193929e-156
[158,] 1.0000000 1.967429e-155 9.837145e-156
[159,] 1.0000000 1.598483e-154 7.992413e-155
[160,] 1.0000000 1.285524e-153 6.427620e-154
[161,] 1.0000000 1.027301e-152 5.136507e-153
[162,] 1.0000000 2.312843e-153 1.156422e-153
[163,] 1.0000000 3.353002e-153 1.676501e-153
[164,] 1.0000000 2.696874e-152 1.348437e-152
[165,] 1.0000000 2.159534e-151 1.079767e-151
[166,] 1.0000000 1.735157e-150 8.675784e-151
[167,] 1.0000000 1.376993e-149 6.884964e-150
[168,] 1.0000000 1.102464e-148 5.512318e-149
[169,] 1.0000000 8.669010e-148 4.334505e-148
[170,] 1.0000000 6.700901e-147 3.350450e-147
[171,] 1.0000000 8.465111e-148 4.232556e-148
[172,] 1.0000000 6.674165e-147 3.337083e-147
[173,] 1.0000000 5.245210e-146 2.622605e-146
[174,] 1.0000000 4.103956e-145 2.051978e-145
[175,] 1.0000000 1.898971e-144 9.494854e-145
[176,] 1.0000000 1.459144e-143 7.295721e-144
[177,] 1.0000000 7.552188e-143 3.776094e-143
[178,] 1.0000000 5.285778e-142 2.642889e-142
[179,] 1.0000000 4.076854e-141 2.038427e-141
[180,] 1.0000000 2.376149e-140 1.188075e-140
[181,] 1.0000000 8.011954e-140 4.005977e-140
[182,] 1.0000000 6.123967e-139 3.061983e-139
[183,] 1.0000000 3.955153e-138 1.977577e-138
[184,] 1.0000000 2.984162e-137 1.492081e-137
[185,] 1.0000000 2.125321e-136 1.062660e-136
[186,] 1.0000000 1.589491e-135 7.947455e-136
[187,] 1.0000000 1.174327e-134 5.871633e-135
[188,] 1.0000000 8.704013e-134 4.352006e-134
[189,] 1.0000000 6.408503e-133 3.204251e-133
[190,] 1.0000000 4.707580e-132 2.353790e-132
[191,] 1.0000000 1.263677e-131 6.318384e-132
[192,] 1.0000000 9.241575e-131 4.620787e-131
[193,] 1.0000000 6.728548e-130 3.364274e-130
[194,] 1.0000000 2.026907e-129 1.013453e-129
[195,] 1.0000000 1.468098e-128 7.340491e-129
[196,] 1.0000000 1.058620e-127 5.293101e-128
[197,] 1.0000000 7.599537e-127 3.799769e-127
[198,] 1.0000000 4.930846e-126 2.465423e-126
[199,] 1.0000000 2.918631e-125 1.459316e-125
[200,] 1.0000000 2.059901e-124 1.029950e-124
[201,] 1.0000000 1.448669e-123 7.243345e-124
[202,] 1.0000000 1.017938e-122 5.089692e-123
[203,] 1.0000000 1.100163e-122 5.500815e-123
[204,] 1.0000000 7.765123e-122 3.882561e-122
[205,] 1.0000000 5.428828e-121 2.714414e-121
[206,] 1.0000000 3.778464e-120 1.889232e-120
[207,] 1.0000000 1.918485e-119 9.592423e-120
[208,] 1.0000000 1.267244e-118 6.336221e-119
[209,] 1.0000000 8.688477e-118 4.344239e-118
[210,] 1.0000000 5.939891e-117 2.969945e-117
[211,] 1.0000000 3.992126e-116 1.996063e-116
[212,] 1.0000000 2.620472e-115 1.310236e-115
[213,] 1.0000000 1.767101e-114 8.835504e-115
[214,] 1.0000000 6.117434e-114 3.058717e-114
[215,] 1.0000000 3.994508e-113 1.997254e-113
[216,] 1.0000000 2.669638e-112 1.334819e-112
[217,] 1.0000000 1.775367e-111 8.876837e-112
[218,] 1.0000000 1.042669e-110 5.213345e-111
[219,] 1.0000000 6.872403e-110 3.436201e-110
[220,] 1.0000000 3.825780e-109 1.912890e-109
[221,] 1.0000000 2.433006e-108 1.216503e-108
[222,] 1.0000000 1.583965e-107 7.919827e-108
[223,] 1.0000000 1.015766e-106 5.078831e-107
[224,] 1.0000000 6.553756e-106 3.276878e-106
[225,] 1.0000000 4.103610e-105 2.051805e-105
[226,] 1.0000000 2.428433e-104 1.214216e-104
[227,] 1.0000000 1.546124e-103 7.730620e-104
[228,] 1.0000000 9.839761e-103 4.919881e-103
[229,] 1.0000000 6.206677e-102 3.103338e-102
[230,] 1.0000000 3.896610e-101 1.948305e-101
[231,] 1.0000000 2.435400e-100 1.217700e-100
[232,] 1.0000000 1.484500e-99 7.422501e-100
[233,] 1.0000000 9.193155e-99 4.596577e-99
[234,] 1.0000000 5.666847e-98 2.833424e-98
[235,] 1.0000000 3.477022e-97 1.738511e-97
[236,] 1.0000000 2.130140e-96 1.065070e-96
[237,] 1.0000000 1.294866e-95 6.474332e-96
[238,] 1.0000000 7.594533e-95 3.797267e-95
[239,] 1.0000000 4.437576e-94 2.218788e-94
[240,] 1.0000000 2.660942e-93 1.330471e-93
[241,] 1.0000000 1.588151e-92 7.940756e-93
[242,] 1.0000000 9.203820e-92 4.601910e-92
[243,] 1.0000000 5.441907e-91 2.720954e-91
[244,] 1.0000000 3.202502e-90 1.601251e-90
[245,] 1.0000000 1.837726e-89 9.188630e-90
[246,] 1.0000000 1.071250e-88 5.356250e-89
[247,] 1.0000000 6.083797e-88 3.041898e-88
[248,] 1.0000000 3.395655e-87 1.697828e-87
[249,] 1.0000000 1.087657e-86 5.438283e-87
[250,] 1.0000000 6.241738e-86 3.120869e-86
[251,] 1.0000000 3.355179e-85 1.677590e-85
[252,] 1.0000000 1.193107e-84 5.965533e-85
[253,] 1.0000000 6.767021e-84 3.383510e-84
[254,] 1.0000000 3.819713e-83 1.909856e-83
[255,] 1.0000000 2.145722e-82 1.072861e-82
[256,] 1.0000000 1.199555e-81 5.997776e-82
[257,] 1.0000000 6.196906e-81 3.098453e-81
[258,] 1.0000000 2.685745e-80 1.342872e-80
[259,] 1.0000000 1.482240e-79 7.411200e-80
[260,] 1.0000000 8.140541e-79 4.070270e-79
[261,] 1.0000000 4.449006e-78 2.224503e-78
[262,] 1.0000000 1.672812e-77 8.364060e-78
[263,] 1.0000000 9.074498e-77 4.537249e-77
[264,] 1.0000000 4.175415e-76 2.087708e-76
[265,] 1.0000000 2.245004e-75 1.122502e-75
[266,] 1.0000000 1.201096e-74 6.005480e-75
[267,] 1.0000000 6.394035e-74 3.197018e-74
[268,] 1.0000000 3.124242e-73 1.562121e-73
[269,] 1.0000000 1.647256e-72 8.236279e-73
[270,] 1.0000000 8.123251e-72 4.061625e-72
[271,] 1.0000000 3.658849e-71 1.829425e-71
[272,] 1.0000000 1.902360e-70 9.511799e-71
[273,] 1.0000000 9.761299e-70 4.880650e-70
[274,] 1.0000000 5.026295e-69 2.513148e-69
[275,] 1.0000000 2.574047e-68 1.287024e-68
[276,] 1.0000000 1.311473e-67 6.557367e-68
[277,] 1.0000000 4.802342e-67 2.401171e-67
[278,] 1.0000000 2.426741e-66 1.213370e-66
[279,] 1.0000000 1.202440e-65 6.012198e-66
[280,] 1.0000000 6.014112e-65 3.007056e-65
[281,] 1.0000000 2.900420e-64 1.450210e-64
[282,] 1.0000000 1.424875e-63 7.124374e-64
[283,] 1.0000000 7.016138e-63 3.508069e-63
[284,] 1.0000000 3.360883e-62 1.680441e-62
[285,] 1.0000000 1.637498e-61 8.187488e-62
[286,] 1.0000000 7.226363e-61 3.613181e-61
[287,] 1.0000000 2.900484e-60 1.450242e-60
[288,] 1.0000000 1.393494e-59 6.967468e-60
[289,] 1.0000000 6.651419e-59 3.325710e-59
[290,] 1.0000000 2.310122e-58 1.155061e-58
[291,] 1.0000000 1.094015e-57 5.470076e-58
[292,] 1.0000000 5.152731e-57 2.576365e-57
[293,] 1.0000000 2.413604e-56 1.206802e-56
[294,] 1.0000000 1.083765e-55 5.418825e-56
[295,] 1.0000000 5.001869e-55 2.500934e-55
[296,] 1.0000000 2.324253e-54 1.162127e-54
[297,] 1.0000000 1.062764e-53 5.313822e-54
[298,] 1.0000000 4.842389e-53 2.421194e-53
[299,] 1.0000000 1.527104e-52 7.635522e-53
[300,] 1.0000000 6.504914e-52 3.252457e-52
[301,] 1.0000000 2.922012e-51 1.461006e-51
[302,] 1.0000000 1.208992e-50 6.044962e-51
[303,] 1.0000000 4.504451e-50 2.252225e-50
[304,] 1.0000000 1.991805e-49 9.959023e-50
[305,] 1.0000000 8.756181e-49 4.378091e-49
[306,] 1.0000000 3.726025e-48 1.863013e-48
[307,] 1.0000000 1.619686e-47 8.098431e-48
[308,] 1.0000000 6.513086e-47 3.256543e-47
[309,] 1.0000000 2.706139e-46 1.353069e-46
[310,] 1.0000000 1.055300e-45 5.276500e-46
[311,] 1.0000000 4.473746e-45 2.236873e-45
[312,] 1.0000000 1.853870e-44 9.269348e-45
[313,] 1.0000000 7.647664e-44 3.823832e-44
[314,] 1.0000000 3.190879e-43 1.595440e-43
[315,] 1.0000000 1.323121e-42 6.615607e-43
[316,] 1.0000000 5.452286e-42 2.726143e-42
[317,] 1.0000000 2.232693e-41 1.116346e-41
[318,] 1.0000000 9.085143e-41 4.542571e-41
[319,] 1.0000000 3.673397e-40 1.836698e-40
[320,] 1.0000000 1.482302e-39 7.411512e-40
[321,] 1.0000000 5.916782e-39 2.958391e-39
[322,] 1.0000000 2.346429e-38 1.173214e-38
[323,] 1.0000000 9.244446e-38 4.622223e-38
[324,] 1.0000000 3.618133e-37 1.809066e-37
[325,] 1.0000000 1.406682e-36 7.033410e-37
[326,] 1.0000000 5.432406e-36 2.716203e-36
[327,] 1.0000000 2.083771e-35 1.041886e-35
[328,] 1.0000000 7.862928e-35 3.931464e-35
[329,] 1.0000000 2.966948e-34 1.483474e-34
[330,] 1.0000000 1.100472e-33 5.502358e-34
[331,] 1.0000000 4.050514e-33 2.025257e-33
[332,] 1.0000000 1.501189e-32 7.505943e-33
[333,] 1.0000000 5.536218e-32 2.768109e-32
[334,] 1.0000000 1.893921e-31 9.469604e-32
[335,] 1.0000000 6.988680e-31 3.494340e-31
[336,] 1.0000000 2.518117e-30 1.259058e-30
[337,] 1.0000000 9.006482e-30 4.503241e-30
[338,] 1.0000000 3.197436e-29 1.598718e-29
[339,] 1.0000000 1.099070e-28 5.495349e-29
[340,] 1.0000000 3.844277e-28 1.922139e-28
[341,] 1.0000000 1.350689e-27 6.753445e-28
[342,] 1.0000000 4.652514e-27 2.326257e-27
[343,] 1.0000000 1.590083e-26 7.950415e-27
[344,] 1.0000000 5.077499e-26 2.538749e-26
[345,] 1.0000000 1.639629e-25 8.198145e-26
[346,] 1.0000000 5.382420e-25 2.691210e-25
[347,] 1.0000000 1.588019e-24 7.940097e-25
[348,] 1.0000000 5.260168e-24 2.630084e-24
[349,] 1.0000000 1.714120e-23 8.570599e-24
[350,] 1.0000000 5.538428e-23 2.769214e-23
[351,] 1.0000000 1.822020e-22 9.110101e-23
[352,] 1.0000000 5.790287e-22 2.895143e-22
[353,] 1.0000000 1.823969e-21 9.119847e-22
[354,] 1.0000000 5.694449e-21 2.847225e-21
[355,] 1.0000000 1.739242e-20 8.696212e-21
[356,] 1.0000000 5.331895e-20 2.665948e-20
[357,] 1.0000000 1.619365e-19 8.096826e-20
[358,] 1.0000000 4.459122e-19 2.229561e-19
[359,] 1.0000000 1.267156e-18 6.335782e-19
[360,] 1.0000000 3.747410e-18 1.873705e-18
[361,] 1.0000000 1.097285e-17 5.486424e-18
[362,] 1.0000000 3.180664e-17 1.590332e-17
[363,] 1.0000000 9.125245e-17 4.562623e-17
[364,] 1.0000000 2.528501e-16 1.264250e-16
[365,] 1.0000000 7.103848e-16 3.551924e-16
[366,] 1.0000000 1.834875e-15 9.174377e-16
[367,] 1.0000000 5.180395e-15 2.590198e-15
[368,] 1.0000000 1.463573e-14 7.317863e-15
[369,] 1.0000000 3.942440e-14 1.971220e-14
[370,] 1.0000000 1.049551e-13 5.247753e-14
[371,] 1.0000000 2.755193e-13 1.377596e-13
[372,] 1.0000000 6.762735e-13 3.381368e-13
[373,] 1.0000000 1.645781e-12 8.228906e-13
[374,] 1.0000000 4.173207e-12 2.086603e-12
[375,] 1.0000000 1.044065e-11 5.220326e-12
[376,] 1.0000000 2.710694e-11 1.355347e-11
[377,] 1.0000000 5.748747e-11 2.874374e-11
[378,] 1.0000000 1.379712e-10 6.898559e-11
[379,] 1.0000000 3.207677e-10 1.603838e-10
[380,] 1.0000000 7.492963e-10 3.746482e-10
[381,] 1.0000000 1.601327e-09 8.006635e-10
[382,] 1.0000000 3.597136e-09 1.798568e-09
[383,] 1.0000000 8.000891e-09 4.000446e-09
[384,] 1.0000000 1.736538e-08 8.682691e-09
[385,] 1.0000000 3.658660e-08 1.829330e-08
[386,] 1.0000000 8.269409e-08 4.134704e-08
[387,] 0.9999999 1.948187e-07 9.740936e-08
[388,] 0.9999998 3.321823e-07 1.660911e-07
[389,] 0.9999996 7.145376e-07 3.572688e-07
[390,] 0.9999992 1.553908e-06 7.769542e-07
[391,] 0.9999983 3.415269e-06 1.707634e-06
[392,] 0.9999989 2.211253e-06 1.105626e-06
[393,] 0.9999975 4.958514e-06 2.479257e-06
[394,] 0.9999969 6.285862e-06 3.142931e-06
[395,] 0.9999929 1.423416e-05 7.117079e-06
[396,] 0.9999865 2.695334e-05 1.347667e-05
[397,] 0.9999730 5.406222e-05 2.703111e-05
[398,] 0.9999415 1.170785e-04 5.853926e-05
[399,] 0.9998766 2.467473e-04 1.233736e-04
[400,] 0.9997595 4.809721e-04 2.404861e-04
[401,] 0.9995423 9.153458e-04 4.576729e-04
[402,] 0.9992143 1.571397e-03 7.856986e-04
[403,] 0.9986712 2.657539e-03 1.328770e-03
[404,] 0.9974800 5.040079e-03 2.520040e-03
[405,] 0.9956407 8.718587e-03 4.359293e-03
[406,] 0.9923942 1.521159e-02 7.605797e-03
[407,] 0.9872296 2.554074e-02 1.277037e-02
[408,] 0.9889947 2.201058e-02 1.100529e-02
[409,] 0.9843607 3.127861e-02 1.563931e-02
[410,] 0.9739277 5.214466e-02 2.607233e-02
[411,] 0.9976557 4.688672e-03 2.344336e-03
[412,] 0.9963847 7.230510e-03 3.615255e-03
[413,] 0.9920132 1.597370e-02 7.986849e-03
[414,] 0.9844017 3.119656e-02 1.559828e-02
[415,] 0.9707689 5.846214e-02 2.923107e-02
[416,] 0.9348681 1.302639e-01 6.513194e-02
[417,] 0.8855443 2.289113e-01 1.144557e-01
[418,] 0.8880029 2.239942e-01 1.119971e-01
> postscript(file="/var/www/rcomp/tmp/1nb4g1292956049.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/rcomp/tmp/2nb4g1292956049.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/rcomp/tmp/3nb4g1292956049.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/rcomp/tmp/4xk311292956049.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/rcomp/tmp/5xk311292956049.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
340675.6874 2594167.7524 1722550.3597 -236477.4865 -156944.7777
6 7 8 9 10
575467.7008 -18734.1681 -178890.7705 88611.2273 841746.1429
11 12 13 14 15
260708.2843 398795.4080 699614.9668 -454857.4185 -341002.5886
16 17 18 19 20
610942.5762 -212053.1222 -475997.8305 201509.3228 -997071.1842
21 22 23 24 25
370357.3971 245843.6206 4400.2931 -1097247.5545 -259403.5402
26 27 28 29 30
153925.6423 -337062.9511 -256461.0148 391564.1640 -190039.5441
31 32 33 34 35
364818.4903 90029.4346 364381.7852 225777.9029 681742.2419
36 37 38 39 40
-665091.9651 147314.3445 168314.3345 -48020.6220 -468840.2425
41 42 43 44 45
-823521.3136 -82380.0097 -49175.0083 -187787.8009 6554.6514
46 47 48 49 50
72345.5298 6328.8551 -460706.9652 34635.5935 -188214.8335
51 52 53 54 55
-37223.1732 -302577.5598 -197893.8234 67199.5523 -514295.7847
56 57 58 59 60
-83186.7776 338758.8715 48801.7867 -67061.1856 -143424.2107
61 62 63 64 65
113473.3182 -221187.1722 101960.8337 -136042.5409 -71852.6127
66 67 68 69 70
-776923.7604 -32851.0573 16382.4429 -131842.4345 -116188.8282
71 72 73 74 75
-222043.2780 -19362.1218 -29906.7408 49218.1317 67348.3632
76 77 78 79 80
-67411.4994 59451.8497 -54129.1308 -106918.2436 105260.6267
81 82 83 84 85
-11619.3933 -38571.8810 46526.3361 -48363.1395 113381.1226
86 87 88 89 90
-108117.6610 -129024.2839 99821.2971 -73269.7527 126925.8198
91 92 93 94 95
40367.0213 2882.4075 -12795.5059 67585.5605 -50190.5068
96 97 98 99 100
14263.1426 204777.8549 1082.3575 157549.7872 31984.5842
101 102 103 104 105
-43643.3516 -60082.1371 -123422.5483 95245.7345 19350.1987
106 107 108 109 110
-79434.6543 -285466.4373 236265.2051 -70196.6536 -70231.3401
111 112 113 114 115
76834.1596 -221896.4443 -95449.6266 -7581.8911 75656.5911
116 117 118 119 120
-35261.7376 68496.3892 152354.7718 -71939.2435 -21325.6342
121 122 123 124 125
-203185.2554 30067.5382 -279276.5402 181548.4371 29069.8249
126 127 128 129 130
-108933.2935 34415.4883 6684.3415 -28116.2424 34415.4883
131 132 133 134 135
14347.2110 15613.6376 34415.4883 -21905.9925 20312.5184
136 137 138 139 140
22880.8169 11285.0025 31222.9958 23150.5184 34415.4883
141 142 143 144 145
20214.4836 34415.4883 -17994.0601 24367.4815 11910.6657
146 147 148 149 150
1139.3479 34415.4883 15311.1015 38885.5034 35452.1628
151 152 153 154 155
-11542.0486 -99299.7376 20918.0259 -5467.2792 29265.6370
156 157 158 159 160
34415.4883 34415.4883 30229.3058 42171.8216 34415.4883
161 162 163 164 165
34415.4883 34415.4883 34415.4883 6584.0554 -40796.9083
166 167 168 169 170
31217.9958 34415.4883 34023.1804 -87866.5980 34415.4883
171 172 173 174 175
34415.4883 24810.4047 34415.4883 6221.6155 34415.4883
176 177 178 179 180
13351.2167 -272740.1617 35257.0689 34415.4883 34415.4883
181 182 183 184 185
40649.1911 35647.8860 -39381.2283 -11864.4442 23409.1469
186 187 188 189 190
22030.6439 13706.8251 18336.7011 5718.6132 34415.4883
191 192 193 194 195
16920.7901 34415.4883 16969.4330 34415.4883 26919.3675
196 197 198 199 200
34415.4883 -158975.0055 34415.4883 34415.4883 -38659.2549
201 202 203 204 205
34415.4883 34415.4883 34415.4883 -26837.7169 50961.7796
206 207 208 209 210
23256.4882 40964.5220 34415.4883 -20082.2236 17885.0109
211 212 213 214 215
34415.4883 34415.4883 14820.3119 -17066.9307 13196.7206
216 217 218 219 220
29016.5034 18708.5506 17755.9585 29016.5034 -38566.9285
221 222 223 224 225
7229.3939 34415.4883 31650.5575 -22793.3612 32064.3713
226 227 228 229 230
-24644.3051 9291.0560 34415.4883 34003.5446 34415.4883
231 232 233 234 235
3980.5053 16116.0560 34415.4883 -11450.5996 34415.4883
236 237 238 239 240
34564.4225 34415.4883 18023.1699 34415.4883 34415.4883
241 242 243 244 245
34415.4883 24932.1768 34415.4883 14324.5227 59095.9927
246 247 248 249 250
34415.4883 34415.4883 -10983.1623 34415.4883 34415.4883
251 252 253 254 255
-20338.8207 34415.4883 29977.0971 -71041.4238 -179472.7167
256 257 258 259 260
34415.4883 204.6311 -66069.3918 34415.4883 34415.4883
261 262 263 264 265
34415.4883 34415.4883 45144.5464 21166.9312 34415.4883
266 267 268 269 270
34415.4883 34415.4883 102136.7639 34415.4883 -10376.3459
271 272 273 274 275
34415.4883 34415.4883 34415.4883 -21956.9647 34415.4883
276 277 278 279 280
9179.9993 -18784.9860 34415.4883 -61114.5573 19816.8185
281 282 283 284 285
34415.4883 34415.4883 -41767.5768 34415.4883 41684.7152
286 287 288 289 290
34415.4883 24748.0259 31329.2606 34415.4883 -32424.7001
291 292 293 294 295
34415.4883 14906.5245 59054.4943 34415.4883 -7392.6895
296 297 298 299 300
33610.4098 34415.4883 34415.4883 34415.4883 10120.0560
301 302 303 304 305
21534.1915 14040.1758 31715.9958 34415.4883 27185.5923
306 307 308 309 310
2057.4240 34415.4883 -11513.0490 -24889.1373 34415.4883
311 312 313 314 315
34415.4883 1629.6336 34415.4883 7999.2017 13160.2732
316 317 318 319 320
56646.4587 32664.8879 -37271.6918 20908.3556 34415.4883
321 322 323 324 325
34415.4883 34415.4883 34415.4883 34415.4883 34415.4883
326 327 328 329 330
28841.5811 34415.4883 34415.4883 34415.4883 34415.4883
331 332 333 334 335
34415.4883 34415.4883 34415.4883 21386.4987 31736.5444
336 337 338 339 340
23608.5184 -38566.9446 34415.4883 30531.5027 8242.3020
341 342 343 344 345
-14503.0610 34415.4883 34415.4883 34415.4883 -7801.9139
346 347 348 349 350
34415.4883 11840.1331 34415.4883 34415.4883 11409.5635
351 352 353 354 355
25952.7483 14838.1667 16032.3686 32957.1139 34415.4883
356 357 358 359 360
34415.4883 25566.6364 34415.4883 34415.4883 34415.4883
361 362 363 364 365
23638.9203 34415.4883 34415.4883 -47363.1526 9446.3909
366 367 368 369 370
34415.4883 34415.4883 34415.4883 34415.4883 15692.0410
371 372 373 374 375
34415.4883 21988.9815 62403.2924 25460.1988 34415.4883
376 377 378 379 380
34415.4883 -22042.5868 22022.0219 11904.5635 34415.4883
381 382 383 384 385
34415.4883 3918.4103 -65207.6130 34415.4883 6855.1773
386 387 388 389 390
-114489.6195 36174.3426 8432.9664 24589.4649 1991.2370
391 392 393 394 395
17387.4017 -7312.9789 -38921.3620 82529.6201 14101.0516
396 397 398 399 400
-22547.5156 -115634.5973 -819140.4195 1469.0903 -183653.8205
401 402 403 404 405
31975.9930 36485.6108 -2751.3511 -48468.7720 -17250.2209
406 407 408 409 410
31178.7587 17219.3911 -11410.3112 -10879.3322 -41168.2135
411 412 413 414 415
-124676.3668 -47696.4693 -56343.7006 93154.6182 -117449.2482
416 417 418 419 420
-59536.7217 35577.8986 -42581.9122 -61195.8229 -27801.2093
421 422 423 424 425
-13972.4958 -245082.7760 -528552.8930 -383023.0073 -68961.9405
426 427 428 429 430
-291280.2801 -306083.6139 -1177595.1299 -97935.0948 -153965.3588
431
-464829.1094
> postscript(file="/var/www/rcomp/tmp/68u241292956049.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 340675.6874 NA
1 2594167.7524 340675.6874
2 1722550.3597 2594167.7524
3 -236477.4865 1722550.3597
4 -156944.7777 -236477.4865
5 575467.7008 -156944.7777
6 -18734.1681 575467.7008
7 -178890.7705 -18734.1681
8 88611.2273 -178890.7705
9 841746.1429 88611.2273
10 260708.2843 841746.1429
11 398795.4080 260708.2843
12 699614.9668 398795.4080
13 -454857.4185 699614.9668
14 -341002.5886 -454857.4185
15 610942.5762 -341002.5886
16 -212053.1222 610942.5762
17 -475997.8305 -212053.1222
18 201509.3228 -475997.8305
19 -997071.1842 201509.3228
20 370357.3971 -997071.1842
21 245843.6206 370357.3971
22 4400.2931 245843.6206
23 -1097247.5545 4400.2931
24 -259403.5402 -1097247.5545
25 153925.6423 -259403.5402
26 -337062.9511 153925.6423
27 -256461.0148 -337062.9511
28 391564.1640 -256461.0148
29 -190039.5441 391564.1640
30 364818.4903 -190039.5441
31 90029.4346 364818.4903
32 364381.7852 90029.4346
33 225777.9029 364381.7852
34 681742.2419 225777.9029
35 -665091.9651 681742.2419
36 147314.3445 -665091.9651
37 168314.3345 147314.3445
38 -48020.6220 168314.3345
39 -468840.2425 -48020.6220
40 -823521.3136 -468840.2425
41 -82380.0097 -823521.3136
42 -49175.0083 -82380.0097
43 -187787.8009 -49175.0083
44 6554.6514 -187787.8009
45 72345.5298 6554.6514
46 6328.8551 72345.5298
47 -460706.9652 6328.8551
48 34635.5935 -460706.9652
49 -188214.8335 34635.5935
50 -37223.1732 -188214.8335
51 -302577.5598 -37223.1732
52 -197893.8234 -302577.5598
53 67199.5523 -197893.8234
54 -514295.7847 67199.5523
55 -83186.7776 -514295.7847
56 338758.8715 -83186.7776
57 48801.7867 338758.8715
58 -67061.1856 48801.7867
59 -143424.2107 -67061.1856
60 113473.3182 -143424.2107
61 -221187.1722 113473.3182
62 101960.8337 -221187.1722
63 -136042.5409 101960.8337
64 -71852.6127 -136042.5409
65 -776923.7604 -71852.6127
66 -32851.0573 -776923.7604
67 16382.4429 -32851.0573
68 -131842.4345 16382.4429
69 -116188.8282 -131842.4345
70 -222043.2780 -116188.8282
71 -19362.1218 -222043.2780
72 -29906.7408 -19362.1218
73 49218.1317 -29906.7408
74 67348.3632 49218.1317
75 -67411.4994 67348.3632
76 59451.8497 -67411.4994
77 -54129.1308 59451.8497
78 -106918.2436 -54129.1308
79 105260.6267 -106918.2436
80 -11619.3933 105260.6267
81 -38571.8810 -11619.3933
82 46526.3361 -38571.8810
83 -48363.1395 46526.3361
84 113381.1226 -48363.1395
85 -108117.6610 113381.1226
86 -129024.2839 -108117.6610
87 99821.2971 -129024.2839
88 -73269.7527 99821.2971
89 126925.8198 -73269.7527
90 40367.0213 126925.8198
91 2882.4075 40367.0213
92 -12795.5059 2882.4075
93 67585.5605 -12795.5059
94 -50190.5068 67585.5605
95 14263.1426 -50190.5068
96 204777.8549 14263.1426
97 1082.3575 204777.8549
98 157549.7872 1082.3575
99 31984.5842 157549.7872
100 -43643.3516 31984.5842
101 -60082.1371 -43643.3516
102 -123422.5483 -60082.1371
103 95245.7345 -123422.5483
104 19350.1987 95245.7345
105 -79434.6543 19350.1987
106 -285466.4373 -79434.6543
107 236265.2051 -285466.4373
108 -70196.6536 236265.2051
109 -70231.3401 -70196.6536
110 76834.1596 -70231.3401
111 -221896.4443 76834.1596
112 -95449.6266 -221896.4443
113 -7581.8911 -95449.6266
114 75656.5911 -7581.8911
115 -35261.7376 75656.5911
116 68496.3892 -35261.7376
117 152354.7718 68496.3892
118 -71939.2435 152354.7718
119 -21325.6342 -71939.2435
120 -203185.2554 -21325.6342
121 30067.5382 -203185.2554
122 -279276.5402 30067.5382
123 181548.4371 -279276.5402
124 29069.8249 181548.4371
125 -108933.2935 29069.8249
126 34415.4883 -108933.2935
127 6684.3415 34415.4883
128 -28116.2424 6684.3415
129 34415.4883 -28116.2424
130 14347.2110 34415.4883
131 15613.6376 14347.2110
132 34415.4883 15613.6376
133 -21905.9925 34415.4883
134 20312.5184 -21905.9925
135 22880.8169 20312.5184
136 11285.0025 22880.8169
137 31222.9958 11285.0025
138 23150.5184 31222.9958
139 34415.4883 23150.5184
140 20214.4836 34415.4883
141 34415.4883 20214.4836
142 -17994.0601 34415.4883
143 24367.4815 -17994.0601
144 11910.6657 24367.4815
145 1139.3479 11910.6657
146 34415.4883 1139.3479
147 15311.1015 34415.4883
148 38885.5034 15311.1015
149 35452.1628 38885.5034
150 -11542.0486 35452.1628
151 -99299.7376 -11542.0486
152 20918.0259 -99299.7376
153 -5467.2792 20918.0259
154 29265.6370 -5467.2792
155 34415.4883 29265.6370
156 34415.4883 34415.4883
157 30229.3058 34415.4883
158 42171.8216 30229.3058
159 34415.4883 42171.8216
160 34415.4883 34415.4883
161 34415.4883 34415.4883
162 34415.4883 34415.4883
163 6584.0554 34415.4883
164 -40796.9083 6584.0554
165 31217.9958 -40796.9083
166 34415.4883 31217.9958
167 34023.1804 34415.4883
168 -87866.5980 34023.1804
169 34415.4883 -87866.5980
170 34415.4883 34415.4883
171 24810.4047 34415.4883
172 34415.4883 24810.4047
173 6221.6155 34415.4883
174 34415.4883 6221.6155
175 13351.2167 34415.4883
176 -272740.1617 13351.2167
177 35257.0689 -272740.1617
178 34415.4883 35257.0689
179 34415.4883 34415.4883
180 40649.1911 34415.4883
181 35647.8860 40649.1911
182 -39381.2283 35647.8860
183 -11864.4442 -39381.2283
184 23409.1469 -11864.4442
185 22030.6439 23409.1469
186 13706.8251 22030.6439
187 18336.7011 13706.8251
188 5718.6132 18336.7011
189 34415.4883 5718.6132
190 16920.7901 34415.4883
191 34415.4883 16920.7901
192 16969.4330 34415.4883
193 34415.4883 16969.4330
194 26919.3675 34415.4883
195 34415.4883 26919.3675
196 -158975.0055 34415.4883
197 34415.4883 -158975.0055
198 34415.4883 34415.4883
199 -38659.2549 34415.4883
200 34415.4883 -38659.2549
201 34415.4883 34415.4883
202 34415.4883 34415.4883
203 -26837.7169 34415.4883
204 50961.7796 -26837.7169
205 23256.4882 50961.7796
206 40964.5220 23256.4882
207 34415.4883 40964.5220
208 -20082.2236 34415.4883
209 17885.0109 -20082.2236
210 34415.4883 17885.0109
211 34415.4883 34415.4883
212 14820.3119 34415.4883
213 -17066.9307 14820.3119
214 13196.7206 -17066.9307
215 29016.5034 13196.7206
216 18708.5506 29016.5034
217 17755.9585 18708.5506
218 29016.5034 17755.9585
219 -38566.9285 29016.5034
220 7229.3939 -38566.9285
221 34415.4883 7229.3939
222 31650.5575 34415.4883
223 -22793.3612 31650.5575
224 32064.3713 -22793.3612
225 -24644.3051 32064.3713
226 9291.0560 -24644.3051
227 34415.4883 9291.0560
228 34003.5446 34415.4883
229 34415.4883 34003.5446
230 3980.5053 34415.4883
231 16116.0560 3980.5053
232 34415.4883 16116.0560
233 -11450.5996 34415.4883
234 34415.4883 -11450.5996
235 34564.4225 34415.4883
236 34415.4883 34564.4225
237 18023.1699 34415.4883
238 34415.4883 18023.1699
239 34415.4883 34415.4883
240 34415.4883 34415.4883
241 24932.1768 34415.4883
242 34415.4883 24932.1768
243 14324.5227 34415.4883
244 59095.9927 14324.5227
245 34415.4883 59095.9927
246 34415.4883 34415.4883
247 -10983.1623 34415.4883
248 34415.4883 -10983.1623
249 34415.4883 34415.4883
250 -20338.8207 34415.4883
251 34415.4883 -20338.8207
252 29977.0971 34415.4883
253 -71041.4238 29977.0971
254 -179472.7167 -71041.4238
255 34415.4883 -179472.7167
256 204.6311 34415.4883
257 -66069.3918 204.6311
258 34415.4883 -66069.3918
259 34415.4883 34415.4883
260 34415.4883 34415.4883
261 34415.4883 34415.4883
262 45144.5464 34415.4883
263 21166.9312 45144.5464
264 34415.4883 21166.9312
265 34415.4883 34415.4883
266 34415.4883 34415.4883
267 102136.7639 34415.4883
268 34415.4883 102136.7639
269 -10376.3459 34415.4883
270 34415.4883 -10376.3459
271 34415.4883 34415.4883
272 34415.4883 34415.4883
273 -21956.9647 34415.4883
274 34415.4883 -21956.9647
275 9179.9993 34415.4883
276 -18784.9860 9179.9993
277 34415.4883 -18784.9860
278 -61114.5573 34415.4883
279 19816.8185 -61114.5573
280 34415.4883 19816.8185
281 34415.4883 34415.4883
282 -41767.5768 34415.4883
283 34415.4883 -41767.5768
284 41684.7152 34415.4883
285 34415.4883 41684.7152
286 24748.0259 34415.4883
287 31329.2606 24748.0259
288 34415.4883 31329.2606
289 -32424.7001 34415.4883
290 34415.4883 -32424.7001
291 14906.5245 34415.4883
292 59054.4943 14906.5245
293 34415.4883 59054.4943
294 -7392.6895 34415.4883
295 33610.4098 -7392.6895
296 34415.4883 33610.4098
297 34415.4883 34415.4883
298 34415.4883 34415.4883
299 10120.0560 34415.4883
300 21534.1915 10120.0560
301 14040.1758 21534.1915
302 31715.9958 14040.1758
303 34415.4883 31715.9958
304 27185.5923 34415.4883
305 2057.4240 27185.5923
306 34415.4883 2057.4240
307 -11513.0490 34415.4883
308 -24889.1373 -11513.0490
309 34415.4883 -24889.1373
310 34415.4883 34415.4883
311 1629.6336 34415.4883
312 34415.4883 1629.6336
313 7999.2017 34415.4883
314 13160.2732 7999.2017
315 56646.4587 13160.2732
316 32664.8879 56646.4587
317 -37271.6918 32664.8879
318 20908.3556 -37271.6918
319 34415.4883 20908.3556
320 34415.4883 34415.4883
321 34415.4883 34415.4883
322 34415.4883 34415.4883
323 34415.4883 34415.4883
324 34415.4883 34415.4883
325 28841.5811 34415.4883
326 34415.4883 28841.5811
327 34415.4883 34415.4883
328 34415.4883 34415.4883
329 34415.4883 34415.4883
330 34415.4883 34415.4883
331 34415.4883 34415.4883
332 34415.4883 34415.4883
333 21386.4987 34415.4883
334 31736.5444 21386.4987
335 23608.5184 31736.5444
336 -38566.9446 23608.5184
337 34415.4883 -38566.9446
338 30531.5027 34415.4883
339 8242.3020 30531.5027
340 -14503.0610 8242.3020
341 34415.4883 -14503.0610
342 34415.4883 34415.4883
343 34415.4883 34415.4883
344 -7801.9139 34415.4883
345 34415.4883 -7801.9139
346 11840.1331 34415.4883
347 34415.4883 11840.1331
348 34415.4883 34415.4883
349 11409.5635 34415.4883
350 25952.7483 11409.5635
351 14838.1667 25952.7483
352 16032.3686 14838.1667
353 32957.1139 16032.3686
354 34415.4883 32957.1139
355 34415.4883 34415.4883
356 25566.6364 34415.4883
357 34415.4883 25566.6364
358 34415.4883 34415.4883
359 34415.4883 34415.4883
360 23638.9203 34415.4883
361 34415.4883 23638.9203
362 34415.4883 34415.4883
363 -47363.1526 34415.4883
364 9446.3909 -47363.1526
365 34415.4883 9446.3909
366 34415.4883 34415.4883
367 34415.4883 34415.4883
368 34415.4883 34415.4883
369 15692.0410 34415.4883
370 34415.4883 15692.0410
371 21988.9815 34415.4883
372 62403.2924 21988.9815
373 25460.1988 62403.2924
374 34415.4883 25460.1988
375 34415.4883 34415.4883
376 -22042.5868 34415.4883
377 22022.0219 -22042.5868
378 11904.5635 22022.0219
379 34415.4883 11904.5635
380 34415.4883 34415.4883
381 3918.4103 34415.4883
382 -65207.6130 3918.4103
383 34415.4883 -65207.6130
384 6855.1773 34415.4883
385 -114489.6195 6855.1773
386 36174.3426 -114489.6195
387 8432.9664 36174.3426
388 24589.4649 8432.9664
389 1991.2370 24589.4649
390 17387.4017 1991.2370
391 -7312.9789 17387.4017
392 -38921.3620 -7312.9789
393 82529.6201 -38921.3620
394 14101.0516 82529.6201
395 -22547.5156 14101.0516
396 -115634.5973 -22547.5156
397 -819140.4195 -115634.5973
398 1469.0903 -819140.4195
399 -183653.8205 1469.0903
400 31975.9930 -183653.8205
401 36485.6108 31975.9930
402 -2751.3511 36485.6108
403 -48468.7720 -2751.3511
404 -17250.2209 -48468.7720
405 31178.7587 -17250.2209
406 17219.3911 31178.7587
407 -11410.3112 17219.3911
408 -10879.3322 -11410.3112
409 -41168.2135 -10879.3322
410 -124676.3668 -41168.2135
411 -47696.4693 -124676.3668
412 -56343.7006 -47696.4693
413 93154.6182 -56343.7006
414 -117449.2482 93154.6182
415 -59536.7217 -117449.2482
416 35577.8986 -59536.7217
417 -42581.9122 35577.8986
418 -61195.8229 -42581.9122
419 -27801.2093 -61195.8229
420 -13972.4958 -27801.2093
421 -245082.7760 -13972.4958
422 -528552.8930 -245082.7760
423 -383023.0073 -528552.8930
424 -68961.9405 -383023.0073
425 -291280.2801 -68961.9405
426 -306083.6139 -291280.2801
427 -1177595.1299 -306083.6139
428 -97935.0948 -1177595.1299
429 -153965.3588 -97935.0948
430 -464829.1094 -153965.3588
431 NA -464829.1094
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2594167.7524 340675.6874
[2,] 1722550.3597 2594167.7524
[3,] -236477.4865 1722550.3597
[4,] -156944.7777 -236477.4865
[5,] 575467.7008 -156944.7777
[6,] -18734.1681 575467.7008
[7,] -178890.7705 -18734.1681
[8,] 88611.2273 -178890.7705
[9,] 841746.1429 88611.2273
[10,] 260708.2843 841746.1429
[11,] 398795.4080 260708.2843
[12,] 699614.9668 398795.4080
[13,] -454857.4185 699614.9668
[14,] -341002.5886 -454857.4185
[15,] 610942.5762 -341002.5886
[16,] -212053.1222 610942.5762
[17,] -475997.8305 -212053.1222
[18,] 201509.3228 -475997.8305
[19,] -997071.1842 201509.3228
[20,] 370357.3971 -997071.1842
[21,] 245843.6206 370357.3971
[22,] 4400.2931 245843.6206
[23,] -1097247.5545 4400.2931
[24,] -259403.5402 -1097247.5545
[25,] 153925.6423 -259403.5402
[26,] -337062.9511 153925.6423
[27,] -256461.0148 -337062.9511
[28,] 391564.1640 -256461.0148
[29,] -190039.5441 391564.1640
[30,] 364818.4903 -190039.5441
[31,] 90029.4346 364818.4903
[32,] 364381.7852 90029.4346
[33,] 225777.9029 364381.7852
[34,] 681742.2419 225777.9029
[35,] -665091.9651 681742.2419
[36,] 147314.3445 -665091.9651
[37,] 168314.3345 147314.3445
[38,] -48020.6220 168314.3345
[39,] -468840.2425 -48020.6220
[40,] -823521.3136 -468840.2425
[41,] -82380.0097 -823521.3136
[42,] -49175.0083 -82380.0097
[43,] -187787.8009 -49175.0083
[44,] 6554.6514 -187787.8009
[45,] 72345.5298 6554.6514
[46,] 6328.8551 72345.5298
[47,] -460706.9652 6328.8551
[48,] 34635.5935 -460706.9652
[49,] -188214.8335 34635.5935
[50,] -37223.1732 -188214.8335
[51,] -302577.5598 -37223.1732
[52,] -197893.8234 -302577.5598
[53,] 67199.5523 -197893.8234
[54,] -514295.7847 67199.5523
[55,] -83186.7776 -514295.7847
[56,] 338758.8715 -83186.7776
[57,] 48801.7867 338758.8715
[58,] -67061.1856 48801.7867
[59,] -143424.2107 -67061.1856
[60,] 113473.3182 -143424.2107
[61,] -221187.1722 113473.3182
[62,] 101960.8337 -221187.1722
[63,] -136042.5409 101960.8337
[64,] -71852.6127 -136042.5409
[65,] -776923.7604 -71852.6127
[66,] -32851.0573 -776923.7604
[67,] 16382.4429 -32851.0573
[68,] -131842.4345 16382.4429
[69,] -116188.8282 -131842.4345
[70,] -222043.2780 -116188.8282
[71,] -19362.1218 -222043.2780
[72,] -29906.7408 -19362.1218
[73,] 49218.1317 -29906.7408
[74,] 67348.3632 49218.1317
[75,] -67411.4994 67348.3632
[76,] 59451.8497 -67411.4994
[77,] -54129.1308 59451.8497
[78,] -106918.2436 -54129.1308
[79,] 105260.6267 -106918.2436
[80,] -11619.3933 105260.6267
[81,] -38571.8810 -11619.3933
[82,] 46526.3361 -38571.8810
[83,] -48363.1395 46526.3361
[84,] 113381.1226 -48363.1395
[85,] -108117.6610 113381.1226
[86,] -129024.2839 -108117.6610
[87,] 99821.2971 -129024.2839
[88,] -73269.7527 99821.2971
[89,] 126925.8198 -73269.7527
[90,] 40367.0213 126925.8198
[91,] 2882.4075 40367.0213
[92,] -12795.5059 2882.4075
[93,] 67585.5605 -12795.5059
[94,] -50190.5068 67585.5605
[95,] 14263.1426 -50190.5068
[96,] 204777.8549 14263.1426
[97,] 1082.3575 204777.8549
[98,] 157549.7872 1082.3575
[99,] 31984.5842 157549.7872
[100,] -43643.3516 31984.5842
[101,] -60082.1371 -43643.3516
[102,] -123422.5483 -60082.1371
[103,] 95245.7345 -123422.5483
[104,] 19350.1987 95245.7345
[105,] -79434.6543 19350.1987
[106,] -285466.4373 -79434.6543
[107,] 236265.2051 -285466.4373
[108,] -70196.6536 236265.2051
[109,] -70231.3401 -70196.6536
[110,] 76834.1596 -70231.3401
[111,] -221896.4443 76834.1596
[112,] -95449.6266 -221896.4443
[113,] -7581.8911 -95449.6266
[114,] 75656.5911 -7581.8911
[115,] -35261.7376 75656.5911
[116,] 68496.3892 -35261.7376
[117,] 152354.7718 68496.3892
[118,] -71939.2435 152354.7718
[119,] -21325.6342 -71939.2435
[120,] -203185.2554 -21325.6342
[121,] 30067.5382 -203185.2554
[122,] -279276.5402 30067.5382
[123,] 181548.4371 -279276.5402
[124,] 29069.8249 181548.4371
[125,] -108933.2935 29069.8249
[126,] 34415.4883 -108933.2935
[127,] 6684.3415 34415.4883
[128,] -28116.2424 6684.3415
[129,] 34415.4883 -28116.2424
[130,] 14347.2110 34415.4883
[131,] 15613.6376 14347.2110
[132,] 34415.4883 15613.6376
[133,] -21905.9925 34415.4883
[134,] 20312.5184 -21905.9925
[135,] 22880.8169 20312.5184
[136,] 11285.0025 22880.8169
[137,] 31222.9958 11285.0025
[138,] 23150.5184 31222.9958
[139,] 34415.4883 23150.5184
[140,] 20214.4836 34415.4883
[141,] 34415.4883 20214.4836
[142,] -17994.0601 34415.4883
[143,] 24367.4815 -17994.0601
[144,] 11910.6657 24367.4815
[145,] 1139.3479 11910.6657
[146,] 34415.4883 1139.3479
[147,] 15311.1015 34415.4883
[148,] 38885.5034 15311.1015
[149,] 35452.1628 38885.5034
[150,] -11542.0486 35452.1628
[151,] -99299.7376 -11542.0486
[152,] 20918.0259 -99299.7376
[153,] -5467.2792 20918.0259
[154,] 29265.6370 -5467.2792
[155,] 34415.4883 29265.6370
[156,] 34415.4883 34415.4883
[157,] 30229.3058 34415.4883
[158,] 42171.8216 30229.3058
[159,] 34415.4883 42171.8216
[160,] 34415.4883 34415.4883
[161,] 34415.4883 34415.4883
[162,] 34415.4883 34415.4883
[163,] 6584.0554 34415.4883
[164,] -40796.9083 6584.0554
[165,] 31217.9958 -40796.9083
[166,] 34415.4883 31217.9958
[167,] 34023.1804 34415.4883
[168,] -87866.5980 34023.1804
[169,] 34415.4883 -87866.5980
[170,] 34415.4883 34415.4883
[171,] 24810.4047 34415.4883
[172,] 34415.4883 24810.4047
[173,] 6221.6155 34415.4883
[174,] 34415.4883 6221.6155
[175,] 13351.2167 34415.4883
[176,] -272740.1617 13351.2167
[177,] 35257.0689 -272740.1617
[178,] 34415.4883 35257.0689
[179,] 34415.4883 34415.4883
[180,] 40649.1911 34415.4883
[181,] 35647.8860 40649.1911
[182,] -39381.2283 35647.8860
[183,] -11864.4442 -39381.2283
[184,] 23409.1469 -11864.4442
[185,] 22030.6439 23409.1469
[186,] 13706.8251 22030.6439
[187,] 18336.7011 13706.8251
[188,] 5718.6132 18336.7011
[189,] 34415.4883 5718.6132
[190,] 16920.7901 34415.4883
[191,] 34415.4883 16920.7901
[192,] 16969.4330 34415.4883
[193,] 34415.4883 16969.4330
[194,] 26919.3675 34415.4883
[195,] 34415.4883 26919.3675
[196,] -158975.0055 34415.4883
[197,] 34415.4883 -158975.0055
[198,] 34415.4883 34415.4883
[199,] -38659.2549 34415.4883
[200,] 34415.4883 -38659.2549
[201,] 34415.4883 34415.4883
[202,] 34415.4883 34415.4883
[203,] -26837.7169 34415.4883
[204,] 50961.7796 -26837.7169
[205,] 23256.4882 50961.7796
[206,] 40964.5220 23256.4882
[207,] 34415.4883 40964.5220
[208,] -20082.2236 34415.4883
[209,] 17885.0109 -20082.2236
[210,] 34415.4883 17885.0109
[211,] 34415.4883 34415.4883
[212,] 14820.3119 34415.4883
[213,] -17066.9307 14820.3119
[214,] 13196.7206 -17066.9307
[215,] 29016.5034 13196.7206
[216,] 18708.5506 29016.5034
[217,] 17755.9585 18708.5506
[218,] 29016.5034 17755.9585
[219,] -38566.9285 29016.5034
[220,] 7229.3939 -38566.9285
[221,] 34415.4883 7229.3939
[222,] 31650.5575 34415.4883
[223,] -22793.3612 31650.5575
[224,] 32064.3713 -22793.3612
[225,] -24644.3051 32064.3713
[226,] 9291.0560 -24644.3051
[227,] 34415.4883 9291.0560
[228,] 34003.5446 34415.4883
[229,] 34415.4883 34003.5446
[230,] 3980.5053 34415.4883
[231,] 16116.0560 3980.5053
[232,] 34415.4883 16116.0560
[233,] -11450.5996 34415.4883
[234,] 34415.4883 -11450.5996
[235,] 34564.4225 34415.4883
[236,] 34415.4883 34564.4225
[237,] 18023.1699 34415.4883
[238,] 34415.4883 18023.1699
[239,] 34415.4883 34415.4883
[240,] 34415.4883 34415.4883
[241,] 24932.1768 34415.4883
[242,] 34415.4883 24932.1768
[243,] 14324.5227 34415.4883
[244,] 59095.9927 14324.5227
[245,] 34415.4883 59095.9927
[246,] 34415.4883 34415.4883
[247,] -10983.1623 34415.4883
[248,] 34415.4883 -10983.1623
[249,] 34415.4883 34415.4883
[250,] -20338.8207 34415.4883
[251,] 34415.4883 -20338.8207
[252,] 29977.0971 34415.4883
[253,] -71041.4238 29977.0971
[254,] -179472.7167 -71041.4238
[255,] 34415.4883 -179472.7167
[256,] 204.6311 34415.4883
[257,] -66069.3918 204.6311
[258,] 34415.4883 -66069.3918
[259,] 34415.4883 34415.4883
[260,] 34415.4883 34415.4883
[261,] 34415.4883 34415.4883
[262,] 45144.5464 34415.4883
[263,] 21166.9312 45144.5464
[264,] 34415.4883 21166.9312
[265,] 34415.4883 34415.4883
[266,] 34415.4883 34415.4883
[267,] 102136.7639 34415.4883
[268,] 34415.4883 102136.7639
[269,] -10376.3459 34415.4883
[270,] 34415.4883 -10376.3459
[271,] 34415.4883 34415.4883
[272,] 34415.4883 34415.4883
[273,] -21956.9647 34415.4883
[274,] 34415.4883 -21956.9647
[275,] 9179.9993 34415.4883
[276,] -18784.9860 9179.9993
[277,] 34415.4883 -18784.9860
[278,] -61114.5573 34415.4883
[279,] 19816.8185 -61114.5573
[280,] 34415.4883 19816.8185
[281,] 34415.4883 34415.4883
[282,] -41767.5768 34415.4883
[283,] 34415.4883 -41767.5768
[284,] 41684.7152 34415.4883
[285,] 34415.4883 41684.7152
[286,] 24748.0259 34415.4883
[287,] 31329.2606 24748.0259
[288,] 34415.4883 31329.2606
[289,] -32424.7001 34415.4883
[290,] 34415.4883 -32424.7001
[291,] 14906.5245 34415.4883
[292,] 59054.4943 14906.5245
[293,] 34415.4883 59054.4943
[294,] -7392.6895 34415.4883
[295,] 33610.4098 -7392.6895
[296,] 34415.4883 33610.4098
[297,] 34415.4883 34415.4883
[298,] 34415.4883 34415.4883
[299,] 10120.0560 34415.4883
[300,] 21534.1915 10120.0560
[301,] 14040.1758 21534.1915
[302,] 31715.9958 14040.1758
[303,] 34415.4883 31715.9958
[304,] 27185.5923 34415.4883
[305,] 2057.4240 27185.5923
[306,] 34415.4883 2057.4240
[307,] -11513.0490 34415.4883
[308,] -24889.1373 -11513.0490
[309,] 34415.4883 -24889.1373
[310,] 34415.4883 34415.4883
[311,] 1629.6336 34415.4883
[312,] 34415.4883 1629.6336
[313,] 7999.2017 34415.4883
[314,] 13160.2732 7999.2017
[315,] 56646.4587 13160.2732
[316,] 32664.8879 56646.4587
[317,] -37271.6918 32664.8879
[318,] 20908.3556 -37271.6918
[319,] 34415.4883 20908.3556
[320,] 34415.4883 34415.4883
[321,] 34415.4883 34415.4883
[322,] 34415.4883 34415.4883
[323,] 34415.4883 34415.4883
[324,] 34415.4883 34415.4883
[325,] 28841.5811 34415.4883
[326,] 34415.4883 28841.5811
[327,] 34415.4883 34415.4883
[328,] 34415.4883 34415.4883
[329,] 34415.4883 34415.4883
[330,] 34415.4883 34415.4883
[331,] 34415.4883 34415.4883
[332,] 34415.4883 34415.4883
[333,] 21386.4987 34415.4883
[334,] 31736.5444 21386.4987
[335,] 23608.5184 31736.5444
[336,] -38566.9446 23608.5184
[337,] 34415.4883 -38566.9446
[338,] 30531.5027 34415.4883
[339,] 8242.3020 30531.5027
[340,] -14503.0610 8242.3020
[341,] 34415.4883 -14503.0610
[342,] 34415.4883 34415.4883
[343,] 34415.4883 34415.4883
[344,] -7801.9139 34415.4883
[345,] 34415.4883 -7801.9139
[346,] 11840.1331 34415.4883
[347,] 34415.4883 11840.1331
[348,] 34415.4883 34415.4883
[349,] 11409.5635 34415.4883
[350,] 25952.7483 11409.5635
[351,] 14838.1667 25952.7483
[352,] 16032.3686 14838.1667
[353,] 32957.1139 16032.3686
[354,] 34415.4883 32957.1139
[355,] 34415.4883 34415.4883
[356,] 25566.6364 34415.4883
[357,] 34415.4883 25566.6364
[358,] 34415.4883 34415.4883
[359,] 34415.4883 34415.4883
[360,] 23638.9203 34415.4883
[361,] 34415.4883 23638.9203
[362,] 34415.4883 34415.4883
[363,] -47363.1526 34415.4883
[364,] 9446.3909 -47363.1526
[365,] 34415.4883 9446.3909
[366,] 34415.4883 34415.4883
[367,] 34415.4883 34415.4883
[368,] 34415.4883 34415.4883
[369,] 15692.0410 34415.4883
[370,] 34415.4883 15692.0410
[371,] 21988.9815 34415.4883
[372,] 62403.2924 21988.9815
[373,] 25460.1988 62403.2924
[374,] 34415.4883 25460.1988
[375,] 34415.4883 34415.4883
[376,] -22042.5868 34415.4883
[377,] 22022.0219 -22042.5868
[378,] 11904.5635 22022.0219
[379,] 34415.4883 11904.5635
[380,] 34415.4883 34415.4883
[381,] 3918.4103 34415.4883
[382,] -65207.6130 3918.4103
[383,] 34415.4883 -65207.6130
[384,] 6855.1773 34415.4883
[385,] -114489.6195 6855.1773
[386,] 36174.3426 -114489.6195
[387,] 8432.9664 36174.3426
[388,] 24589.4649 8432.9664
[389,] 1991.2370 24589.4649
[390,] 17387.4017 1991.2370
[391,] -7312.9789 17387.4017
[392,] -38921.3620 -7312.9789
[393,] 82529.6201 -38921.3620
[394,] 14101.0516 82529.6201
[395,] -22547.5156 14101.0516
[396,] -115634.5973 -22547.5156
[397,] -819140.4195 -115634.5973
[398,] 1469.0903 -819140.4195
[399,] -183653.8205 1469.0903
[400,] 31975.9930 -183653.8205
[401,] 36485.6108 31975.9930
[402,] -2751.3511 36485.6108
[403,] -48468.7720 -2751.3511
[404,] -17250.2209 -48468.7720
[405,] 31178.7587 -17250.2209
[406,] 17219.3911 31178.7587
[407,] -11410.3112 17219.3911
[408,] -10879.3322 -11410.3112
[409,] -41168.2135 -10879.3322
[410,] -124676.3668 -41168.2135
[411,] -47696.4693 -124676.3668
[412,] -56343.7006 -47696.4693
[413,] 93154.6182 -56343.7006
[414,] -117449.2482 93154.6182
[415,] -59536.7217 -117449.2482
[416,] 35577.8986 -59536.7217
[417,] -42581.9122 35577.8986
[418,] -61195.8229 -42581.9122
[419,] -27801.2093 -61195.8229
[420,] -13972.4958 -27801.2093
[421,] -245082.7760 -13972.4958
[422,] -528552.8930 -245082.7760
[423,] -383023.0073 -528552.8930
[424,] -68961.9405 -383023.0073
[425,] -291280.2801 -68961.9405
[426,] -306083.6139 -291280.2801
[427,] -1177595.1299 -306083.6139
[428,] -97935.0948 -1177595.1299
[429,] -153965.3588 -97935.0948
[430,] -464829.1094 -153965.3588
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2594167.7524 340675.6874
2 1722550.3597 2594167.7524
3 -236477.4865 1722550.3597
4 -156944.7777 -236477.4865
5 575467.7008 -156944.7777
6 -18734.1681 575467.7008
7 -178890.7705 -18734.1681
8 88611.2273 -178890.7705
9 841746.1429 88611.2273
10 260708.2843 841746.1429
11 398795.4080 260708.2843
12 699614.9668 398795.4080
13 -454857.4185 699614.9668
14 -341002.5886 -454857.4185
15 610942.5762 -341002.5886
16 -212053.1222 610942.5762
17 -475997.8305 -212053.1222
18 201509.3228 -475997.8305
19 -997071.1842 201509.3228
20 370357.3971 -997071.1842
21 245843.6206 370357.3971
22 4400.2931 245843.6206
23 -1097247.5545 4400.2931
24 -259403.5402 -1097247.5545
25 153925.6423 -259403.5402
26 -337062.9511 153925.6423
27 -256461.0148 -337062.9511
28 391564.1640 -256461.0148
29 -190039.5441 391564.1640
30 364818.4903 -190039.5441
31 90029.4346 364818.4903
32 364381.7852 90029.4346
33 225777.9029 364381.7852
34 681742.2419 225777.9029
35 -665091.9651 681742.2419
36 147314.3445 -665091.9651
37 168314.3345 147314.3445
38 -48020.6220 168314.3345
39 -468840.2425 -48020.6220
40 -823521.3136 -468840.2425
41 -82380.0097 -823521.3136
42 -49175.0083 -82380.0097
43 -187787.8009 -49175.0083
44 6554.6514 -187787.8009
45 72345.5298 6554.6514
46 6328.8551 72345.5298
47 -460706.9652 6328.8551
48 34635.5935 -460706.9652
49 -188214.8335 34635.5935
50 -37223.1732 -188214.8335
51 -302577.5598 -37223.1732
52 -197893.8234 -302577.5598
53 67199.5523 -197893.8234
54 -514295.7847 67199.5523
55 -83186.7776 -514295.7847
56 338758.8715 -83186.7776
57 48801.7867 338758.8715
58 -67061.1856 48801.7867
59 -143424.2107 -67061.1856
60 113473.3182 -143424.2107
61 -221187.1722 113473.3182
62 101960.8337 -221187.1722
63 -136042.5409 101960.8337
64 -71852.6127 -136042.5409
65 -776923.7604 -71852.6127
66 -32851.0573 -776923.7604
67 16382.4429 -32851.0573
68 -131842.4345 16382.4429
69 -116188.8282 -131842.4345
70 -222043.2780 -116188.8282
71 -19362.1218 -222043.2780
72 -29906.7408 -19362.1218
73 49218.1317 -29906.7408
74 67348.3632 49218.1317
75 -67411.4994 67348.3632
76 59451.8497 -67411.4994
77 -54129.1308 59451.8497
78 -106918.2436 -54129.1308
79 105260.6267 -106918.2436
80 -11619.3933 105260.6267
81 -38571.8810 -11619.3933
82 46526.3361 -38571.8810
83 -48363.1395 46526.3361
84 113381.1226 -48363.1395
85 -108117.6610 113381.1226
86 -129024.2839 -108117.6610
87 99821.2971 -129024.2839
88 -73269.7527 99821.2971
89 126925.8198 -73269.7527
90 40367.0213 126925.8198
91 2882.4075 40367.0213
92 -12795.5059 2882.4075
93 67585.5605 -12795.5059
94 -50190.5068 67585.5605
95 14263.1426 -50190.5068
96 204777.8549 14263.1426
97 1082.3575 204777.8549
98 157549.7872 1082.3575
99 31984.5842 157549.7872
100 -43643.3516 31984.5842
101 -60082.1371 -43643.3516
102 -123422.5483 -60082.1371
103 95245.7345 -123422.5483
104 19350.1987 95245.7345
105 -79434.6543 19350.1987
106 -285466.4373 -79434.6543
107 236265.2051 -285466.4373
108 -70196.6536 236265.2051
109 -70231.3401 -70196.6536
110 76834.1596 -70231.3401
111 -221896.4443 76834.1596
112 -95449.6266 -221896.4443
113 -7581.8911 -95449.6266
114 75656.5911 -7581.8911
115 -35261.7376 75656.5911
116 68496.3892 -35261.7376
117 152354.7718 68496.3892
118 -71939.2435 152354.7718
119 -21325.6342 -71939.2435
120 -203185.2554 -21325.6342
121 30067.5382 -203185.2554
122 -279276.5402 30067.5382
123 181548.4371 -279276.5402
124 29069.8249 181548.4371
125 -108933.2935 29069.8249
126 34415.4883 -108933.2935
127 6684.3415 34415.4883
128 -28116.2424 6684.3415
129 34415.4883 -28116.2424
130 14347.2110 34415.4883
131 15613.6376 14347.2110
132 34415.4883 15613.6376
133 -21905.9925 34415.4883
134 20312.5184 -21905.9925
135 22880.8169 20312.5184
136 11285.0025 22880.8169
137 31222.9958 11285.0025
138 23150.5184 31222.9958
139 34415.4883 23150.5184
140 20214.4836 34415.4883
141 34415.4883 20214.4836
142 -17994.0601 34415.4883
143 24367.4815 -17994.0601
144 11910.6657 24367.4815
145 1139.3479 11910.6657
146 34415.4883 1139.3479
147 15311.1015 34415.4883
148 38885.5034 15311.1015
149 35452.1628 38885.5034
150 -11542.0486 35452.1628
151 -99299.7376 -11542.0486
152 20918.0259 -99299.7376
153 -5467.2792 20918.0259
154 29265.6370 -5467.2792
155 34415.4883 29265.6370
156 34415.4883 34415.4883
157 30229.3058 34415.4883
158 42171.8216 30229.3058
159 34415.4883 42171.8216
160 34415.4883 34415.4883
161 34415.4883 34415.4883
162 34415.4883 34415.4883
163 6584.0554 34415.4883
164 -40796.9083 6584.0554
165 31217.9958 -40796.9083
166 34415.4883 31217.9958
167 34023.1804 34415.4883
168 -87866.5980 34023.1804
169 34415.4883 -87866.5980
170 34415.4883 34415.4883
171 24810.4047 34415.4883
172 34415.4883 24810.4047
173 6221.6155 34415.4883
174 34415.4883 6221.6155
175 13351.2167 34415.4883
176 -272740.1617 13351.2167
177 35257.0689 -272740.1617
178 34415.4883 35257.0689
179 34415.4883 34415.4883
180 40649.1911 34415.4883
181 35647.8860 40649.1911
182 -39381.2283 35647.8860
183 -11864.4442 -39381.2283
184 23409.1469 -11864.4442
185 22030.6439 23409.1469
186 13706.8251 22030.6439
187 18336.7011 13706.8251
188 5718.6132 18336.7011
189 34415.4883 5718.6132
190 16920.7901 34415.4883
191 34415.4883 16920.7901
192 16969.4330 34415.4883
193 34415.4883 16969.4330
194 26919.3675 34415.4883
195 34415.4883 26919.3675
196 -158975.0055 34415.4883
197 34415.4883 -158975.0055
198 34415.4883 34415.4883
199 -38659.2549 34415.4883
200 34415.4883 -38659.2549
201 34415.4883 34415.4883
202 34415.4883 34415.4883
203 -26837.7169 34415.4883
204 50961.7796 -26837.7169
205 23256.4882 50961.7796
206 40964.5220 23256.4882
207 34415.4883 40964.5220
208 -20082.2236 34415.4883
209 17885.0109 -20082.2236
210 34415.4883 17885.0109
211 34415.4883 34415.4883
212 14820.3119 34415.4883
213 -17066.9307 14820.3119
214 13196.7206 -17066.9307
215 29016.5034 13196.7206
216 18708.5506 29016.5034
217 17755.9585 18708.5506
218 29016.5034 17755.9585
219 -38566.9285 29016.5034
220 7229.3939 -38566.9285
221 34415.4883 7229.3939
222 31650.5575 34415.4883
223 -22793.3612 31650.5575
224 32064.3713 -22793.3612
225 -24644.3051 32064.3713
226 9291.0560 -24644.3051
227 34415.4883 9291.0560
228 34003.5446 34415.4883
229 34415.4883 34003.5446
230 3980.5053 34415.4883
231 16116.0560 3980.5053
232 34415.4883 16116.0560
233 -11450.5996 34415.4883
234 34415.4883 -11450.5996
235 34564.4225 34415.4883
236 34415.4883 34564.4225
237 18023.1699 34415.4883
238 34415.4883 18023.1699
239 34415.4883 34415.4883
240 34415.4883 34415.4883
241 24932.1768 34415.4883
242 34415.4883 24932.1768
243 14324.5227 34415.4883
244 59095.9927 14324.5227
245 34415.4883 59095.9927
246 34415.4883 34415.4883
247 -10983.1623 34415.4883
248 34415.4883 -10983.1623
249 34415.4883 34415.4883
250 -20338.8207 34415.4883
251 34415.4883 -20338.8207
252 29977.0971 34415.4883
253 -71041.4238 29977.0971
254 -179472.7167 -71041.4238
255 34415.4883 -179472.7167
256 204.6311 34415.4883
257 -66069.3918 204.6311
258 34415.4883 -66069.3918
259 34415.4883 34415.4883
260 34415.4883 34415.4883
261 34415.4883 34415.4883
262 45144.5464 34415.4883
263 21166.9312 45144.5464
264 34415.4883 21166.9312
265 34415.4883 34415.4883
266 34415.4883 34415.4883
267 102136.7639 34415.4883
268 34415.4883 102136.7639
269 -10376.3459 34415.4883
270 34415.4883 -10376.3459
271 34415.4883 34415.4883
272 34415.4883 34415.4883
273 -21956.9647 34415.4883
274 34415.4883 -21956.9647
275 9179.9993 34415.4883
276 -18784.9860 9179.9993
277 34415.4883 -18784.9860
278 -61114.5573 34415.4883
279 19816.8185 -61114.5573
280 34415.4883 19816.8185
281 34415.4883 34415.4883
282 -41767.5768 34415.4883
283 34415.4883 -41767.5768
284 41684.7152 34415.4883
285 34415.4883 41684.7152
286 24748.0259 34415.4883
287 31329.2606 24748.0259
288 34415.4883 31329.2606
289 -32424.7001 34415.4883
290 34415.4883 -32424.7001
291 14906.5245 34415.4883
292 59054.4943 14906.5245
293 34415.4883 59054.4943
294 -7392.6895 34415.4883
295 33610.4098 -7392.6895
296 34415.4883 33610.4098
297 34415.4883 34415.4883
298 34415.4883 34415.4883
299 10120.0560 34415.4883
300 21534.1915 10120.0560
301 14040.1758 21534.1915
302 31715.9958 14040.1758
303 34415.4883 31715.9958
304 27185.5923 34415.4883
305 2057.4240 27185.5923
306 34415.4883 2057.4240
307 -11513.0490 34415.4883
308 -24889.1373 -11513.0490
309 34415.4883 -24889.1373
310 34415.4883 34415.4883
311 1629.6336 34415.4883
312 34415.4883 1629.6336
313 7999.2017 34415.4883
314 13160.2732 7999.2017
315 56646.4587 13160.2732
316 32664.8879 56646.4587
317 -37271.6918 32664.8879
318 20908.3556 -37271.6918
319 34415.4883 20908.3556
320 34415.4883 34415.4883
321 34415.4883 34415.4883
322 34415.4883 34415.4883
323 34415.4883 34415.4883
324 34415.4883 34415.4883
325 28841.5811 34415.4883
326 34415.4883 28841.5811
327 34415.4883 34415.4883
328 34415.4883 34415.4883
329 34415.4883 34415.4883
330 34415.4883 34415.4883
331 34415.4883 34415.4883
332 34415.4883 34415.4883
333 21386.4987 34415.4883
334 31736.5444 21386.4987
335 23608.5184 31736.5444
336 -38566.9446 23608.5184
337 34415.4883 -38566.9446
338 30531.5027 34415.4883
339 8242.3020 30531.5027
340 -14503.0610 8242.3020
341 34415.4883 -14503.0610
342 34415.4883 34415.4883
343 34415.4883 34415.4883
344 -7801.9139 34415.4883
345 34415.4883 -7801.9139
346 11840.1331 34415.4883
347 34415.4883 11840.1331
348 34415.4883 34415.4883
349 11409.5635 34415.4883
350 25952.7483 11409.5635
351 14838.1667 25952.7483
352 16032.3686 14838.1667
353 32957.1139 16032.3686
354 34415.4883 32957.1139
355 34415.4883 34415.4883
356 25566.6364 34415.4883
357 34415.4883 25566.6364
358 34415.4883 34415.4883
359 34415.4883 34415.4883
360 23638.9203 34415.4883
361 34415.4883 23638.9203
362 34415.4883 34415.4883
363 -47363.1526 34415.4883
364 9446.3909 -47363.1526
365 34415.4883 9446.3909
366 34415.4883 34415.4883
367 34415.4883 34415.4883
368 34415.4883 34415.4883
369 15692.0410 34415.4883
370 34415.4883 15692.0410
371 21988.9815 34415.4883
372 62403.2924 21988.9815
373 25460.1988 62403.2924
374 34415.4883 25460.1988
375 34415.4883 34415.4883
376 -22042.5868 34415.4883
377 22022.0219 -22042.5868
378 11904.5635 22022.0219
379 34415.4883 11904.5635
380 34415.4883 34415.4883
381 3918.4103 34415.4883
382 -65207.6130 3918.4103
383 34415.4883 -65207.6130
384 6855.1773 34415.4883
385 -114489.6195 6855.1773
386 36174.3426 -114489.6195
387 8432.9664 36174.3426
388 24589.4649 8432.9664
389 1991.2370 24589.4649
390 17387.4017 1991.2370
391 -7312.9789 17387.4017
392 -38921.3620 -7312.9789
393 82529.6201 -38921.3620
394 14101.0516 82529.6201
395 -22547.5156 14101.0516
396 -115634.5973 -22547.5156
397 -819140.4195 -115634.5973
398 1469.0903 -819140.4195
399 -183653.8205 1469.0903
400 31975.9930 -183653.8205
401 36485.6108 31975.9930
402 -2751.3511 36485.6108
403 -48468.7720 -2751.3511
404 -17250.2209 -48468.7720
405 31178.7587 -17250.2209
406 17219.3911 31178.7587
407 -11410.3112 17219.3911
408 -10879.3322 -11410.3112
409 -41168.2135 -10879.3322
410 -124676.3668 -41168.2135
411 -47696.4693 -124676.3668
412 -56343.7006 -47696.4693
413 93154.6182 -56343.7006
414 -117449.2482 93154.6182
415 -59536.7217 -117449.2482
416 35577.8986 -59536.7217
417 -42581.9122 35577.8986
418 -61195.8229 -42581.9122
419 -27801.2093 -61195.8229
420 -13972.4958 -27801.2093
421 -245082.7760 -13972.4958
422 -528552.8930 -245082.7760
423 -383023.0073 -528552.8930
424 -68961.9405 -383023.0073
425 -291280.2801 -68961.9405
426 -306083.6139 -291280.2801
427 -1177595.1299 -306083.6139
428 -97935.0948 -1177595.1299
429 -153965.3588 -97935.0948
430 -464829.1094 -153965.3588
> 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/rcomp/tmp/7jljo1292956049.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/rcomp/tmp/8jljo1292956049.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/rcomp/tmp/9jljo1292956049.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/rcomp/tmp/10jljo1292956049.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/114l0c1292956049.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/rcomp/tmp/12bn0y1292956050.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/rcomp/tmp/13oxg71292956050.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/rcomp/tmp/14h6xs1292956050.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/rcomp/tmp/15dydi1292956050.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/rcomp/tmp/16hyb61292956050.tab")
+ }
>
> try(system("convert tmp/1nb4g1292956049.ps tmp/1nb4g1292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nb4g1292956049.ps tmp/2nb4g1292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/3nb4g1292956049.ps tmp/3nb4g1292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/4xk311292956049.ps tmp/4xk311292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/5xk311292956049.ps tmp/5xk311292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/68u241292956049.ps tmp/68u241292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/7jljo1292956049.ps tmp/7jljo1292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/8jljo1292956049.ps tmp/8jljo1292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/9jljo1292956049.ps tmp/9jljo1292956049.png",intern=TRUE))
character(0)
> try(system("convert tmp/10jljo1292956049.ps tmp/10jljo1292956049.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.360 1.530 12.892