R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(162556
+ ,807
+ ,213118
+ ,6282154
+ ,29790
+ ,444
+ ,81767
+ ,4321023
+ ,87550
+ ,412
+ ,153198
+ ,4111912
+ ,84738
+ ,428
+ ,-26007
+ ,223193
+ ,54660
+ ,315
+ ,126942
+ ,1491348
+ ,42634
+ ,168
+ ,157214
+ ,1629616
+ ,40949
+ ,263
+ ,129352
+ ,1398893
+ ,45187
+ ,267
+ ,234817
+ ,1926517
+ ,37704
+ ,228
+ ,60448
+ ,983660
+ ,16275
+ ,129
+ ,47818
+ ,1443586
+ ,25830
+ ,104
+ ,245546
+ ,1073089
+ ,12679
+ ,122
+ ,48020
+ ,984885
+ ,18014
+ ,393
+ ,-1710
+ ,1405225
+ ,43556
+ ,190
+ ,32648
+ ,227132
+ ,24811
+ ,280
+ ,95350
+ ,929118
+ ,6575
+ ,63
+ ,151352
+ ,1071292
+ ,7123
+ ,102
+ ,288170
+ ,638830
+ ,21950
+ ,265
+ ,114337
+ ,856956
+ ,37597
+ ,234
+ ,37884
+ ,992426
+ ,17821
+ ,277
+ ,122844
+ ,444477
+ ,12988
+ ,73
+ ,82340
+ ,857217
+ ,22330
+ ,67
+ ,79801
+ ,711969
+ ,13326
+ ,103
+ ,165548
+ ,702380
+ ,16189
+ ,290
+ ,116384
+ ,358589
+ ,7146
+ ,83
+ ,134028
+ ,297978
+ ,15824
+ ,56
+ ,63838
+ ,585715
+ ,27664
+ ,236
+ ,74996
+ ,657954
+ ,11920
+ ,73
+ ,31080
+ ,209458
+ ,8568
+ ,34
+ ,32168
+ ,786690
+ ,14416
+ ,139
+ ,49857
+ ,439798
+ ,3369
+ ,26
+ ,87161
+ ,688779
+ ,11819
+ ,70
+ ,106113
+ ,574339
+ ,6984
+ ,40
+ ,80570
+ ,741409
+ ,4519
+ ,42
+ ,102129
+ ,597793
+ ,2220
+ ,12
+ ,301670
+ ,644190
+ ,18562
+ ,211
+ ,102313
+ ,377934
+ ,10327
+ ,74
+ ,88577
+ ,640273
+ ,5336
+ ,80
+ ,112477
+ ,697458
+ ,2365
+ ,83
+ ,191778
+ ,550608
+ ,4069
+ ,131
+ ,79804
+ ,207393
+ ,8636
+ ,203
+ ,128294
+ ,301607
+ ,13718
+ ,56
+ ,96448
+ ,345783
+ ,4525
+ ,89
+ ,93811
+ ,501749
+ ,6869
+ ,88
+ ,117520
+ ,379983
+ ,4628
+ ,39
+ ,69159
+ ,387475
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,4891
+ ,49
+ ,210568
+ ,370837
+ ,7489
+ ,149
+ ,136996
+ ,430866
+ ,4901
+ ,58
+ ,121920
+ ,469107
+ ,2284
+ ,41
+ ,76403
+ ,194493
+ ,3160
+ ,90
+ ,108094
+ ,530670
+ ,4150
+ ,136
+ ,134759
+ ,518365
+ ,7285
+ ,97
+ ,188873
+ ,491303
+ ,1134
+ ,63
+ ,146216
+ ,527021
+ ,4658
+ ,114
+ ,156608
+ ,233773
+ ,2384
+ ,77
+ ,61348
+ ,405972
+ ,3748
+ ,6
+ ,50350
+ ,652925
+ ,5371
+ ,47
+ ,87720
+ ,446211
+ ,1285
+ ,51
+ ,99489
+ ,341340
+ ,9327
+ ,85
+ ,87419
+ ,387699
+ ,5565
+ ,43
+ ,94355
+ ,493408
+ ,1528
+ ,32
+ ,60326
+ ,146494
+ ,3122
+ ,25
+ ,94670
+ ,414462
+ ,7561
+ ,77
+ ,82425
+ ,364304
+ ,2675
+ ,54
+ ,59017
+ ,355178
+ ,13253
+ ,251
+ ,90829
+ ,357760
+ ,880
+ ,15
+ ,80791
+ ,261216
+ ,2053
+ ,44
+ ,100423
+ ,397144
+ ,1424
+ ,73
+ ,131116
+ ,374943
+ ,4036
+ ,85
+ ,100269
+ ,424898
+ ,3045
+ ,49
+ ,27330
+ ,202055
+ ,5119
+ ,38
+ ,39039
+ ,378525
+ ,1431
+ ,35
+ ,106885
+ ,310768
+ ,554
+ ,9
+ ,79285
+ ,325738
+ ,1975
+ ,34
+ ,118881
+ ,394510
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,1012
+ ,29
+ ,114768
+ ,368078
+ ,810
+ ,11
+ ,74015
+ ,236761
+ ,1280
+ ,52
+ ,69465
+ ,312378
+ ,666
+ ,13
+ ,117869
+ ,339836
+ ,1380
+ ,29
+ ,60982
+ ,347385
+ ,4677
+ ,66
+ ,90131
+ ,426280
+ ,876
+ ,33
+ ,138971
+ ,352850
+ ,814
+ ,15
+ ,39625
+ ,301881
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,5692
+ ,68
+ ,64239
+ ,357312
+ ,3642
+ ,100
+ ,90262
+ ,458343
+ ,540
+ ,13
+ ,103960
+ ,354228
+ ,2099
+ ,45
+ ,106611
+ ,308636
+ ,567
+ ,14
+ ,103345
+ ,386212
+ ,2001
+ ,36
+ ,95551
+ ,393343
+ ,2949
+ ,40
+ ,82903
+ ,378509
+ ,2253
+ ,68
+ ,63593
+ ,452469
+ ,6533
+ ,29
+ ,126910
+ ,364839
+ ,1889
+ ,43
+ ,37527
+ ,358649
+ ,3055
+ ,30
+ ,60247
+ ,376641
+ ,272
+ ,9
+ ,112995
+ ,429112
+ ,1414
+ ,22
+ ,70184
+ ,330546
+ ,2564
+ ,19
+ ,130140
+ ,403560
+ ,1383
+ ,9
+ ,73221
+ ,317892
+ ,1261
+ ,31
+ ,76114
+ ,307528
+ ,975
+ ,19
+ ,90534
+ ,235133
+ ,3366
+ ,55
+ ,108479
+ ,299243
+ ,576
+ ,8
+ ,113761
+ ,314073
+ ,1686
+ ,28
+ ,68696
+ ,368186
+ ,746
+ ,29
+ ,71561
+ ,269661
+ ,3192
+ ,48
+ ,59831
+ ,125390
+ ,2045
+ ,16
+ ,97890
+ ,510834
+ ,5702
+ ,47
+ ,101481
+ ,321896
+ ,1932
+ ,20
+ ,72954
+ ,249898
+ ,936
+ ,22
+ ,67939
+ ,408881
+ ,3437
+ ,33
+ ,48022
+ ,158492
+ ,5131
+ ,44
+ ,86111
+ ,292154
+ ,2397
+ ,13
+ ,74020
+ ,289513
+ ,1389
+ ,6
+ ,57530
+ ,378049
+ ,1503
+ ,35
+ ,56364
+ ,343466
+ ,402
+ ,8
+ ,84990
+ ,332743
+ ,2239
+ ,17
+ ,88590
+ ,442882
+ ,2234
+ ,11
+ ,77200
+ ,214215
+ ,837
+ ,21
+ ,61262
+ ,315688
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,875
+ ,12
+ ,67000
+ ,334280
+ ,1585
+ ,112
+ ,93099
+ ,355864
+ ,1659
+ ,25
+ ,107577
+ ,480382
+ ,2647
+ ,17
+ ,62920
+ ,353058
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,94
+ ,10
+ ,60793
+ ,314533
+ ,422
+ ,23
+ ,57935
+ ,318056
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,34
+ ,7
+ ,60630
+ ,314353
+ ,1558
+ ,25
+ ,55637
+ ,369448
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,43
+ ,20
+ ,60887
+ ,312846
+ ,645
+ ,4
+ ,60720
+ ,312075
+ ,316
+ ,4
+ ,60505
+ ,315009
+ ,115
+ ,10
+ ,60945
+ ,318903
+ ,5
+ ,1
+ ,60720
+ ,314887
+ ,897
+ ,4
+ ,60720
+ ,314913
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,389
+ ,8
+ ,58990
+ ,325506
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1002
+ ,11
+ ,56750
+ ,298568
+ ,36
+ ,4
+ ,60894
+ ,315834
+ ,460
+ ,15
+ ,63346
+ ,329784
+ ,309
+ ,9
+ ,56535
+ ,312878
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,9
+ ,7
+ ,60835
+ ,314987
+ ,271
+ ,2
+ ,60720
+ ,325249
+ ,14
+ ,0
+ ,61016
+ ,315877
+ ,520
+ ,7
+ ,58650
+ ,291650
+ ,1766
+ ,46
+ ,60438
+ ,305959
+ ,0
+ ,5
+ ,60720
+ ,315380
+ ,458
+ ,7
+ ,58625
+ ,297765
+ ,20
+ ,2
+ ,60938
+ ,315245
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,98
+ ,2
+ ,61490
+ ,315236
+ ,405
+ ,5
+ ,60845
+ ,336425
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,483
+ ,7
+ ,60830
+ ,306268
+ ,454
+ ,24
+ ,63261
+ ,302187
+ ,47
+ ,1
+ ,60720
+ ,314882
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,757
+ ,18
+ ,45689
+ ,382712
+ ,4655
+ ,55
+ ,60720
+ ,341570
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,36
+ ,3
+ ,61564
+ ,312412
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,9
+ ,61938
+ ,309596
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,126
+ ,8
+ ,60951
+ ,315547
+ ,400
+ ,113
+ ,60720
+ ,313267
+ ,71
+ ,0
+ ,60745
+ ,316176
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,972
+ ,19
+ ,71642
+ ,359335
+ ,531
+ ,11
+ ,71641
+ ,330068
+ ,2461
+ ,25
+ ,55792
+ ,314289
+ ,378
+ ,16
+ ,71873
+ ,297413
+ ,23
+ ,5
+ ,62555
+ ,314806
+ ,638
+ ,11
+ ,60370
+ ,333210
+ ,2300
+ ,23
+ ,64873
+ ,352108
+ ,149
+ ,6
+ ,62041
+ ,313332
+ ,226
+ ,5
+ ,65745
+ ,291787
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,275
+ ,7
+ ,59500
+ ,318745
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,141
+ ,7
+ ,61630
+ ,315366
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,28
+ ,3
+ ,60890
+ ,315688
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,4980
+ ,89
+ ,113521
+ ,409642
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,472
+ ,19
+ ,80045
+ ,269587
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,203
+ ,12
+ ,50804
+ ,300962
+ ,496
+ ,12
+ ,87390
+ ,325479
+ ,10
+ ,5
+ ,61656
+ ,316155
+ ,63
+ ,2
+ ,65688
+ ,318574
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1136
+ ,26
+ ,48522
+ ,343613
+ ,265
+ ,3
+ ,60720
+ ,306948
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,267
+ ,11
+ ,57640
+ ,330059
+ ,474
+ ,10
+ ,61977
+ ,288985
+ ,534
+ ,5
+ ,62620
+ ,304485
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,15
+ ,6
+ ,60831
+ ,315688
+ ,397
+ ,7
+ ,60646
+ ,317736
+ ,0
+ ,2
+ ,60720
+ ,315380
+ ,1866
+ ,28
+ ,56225
+ ,322331
+ ,288
+ ,3
+ ,60510
+ ,296656
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,3
+ ,1
+ ,60698
+ ,315354
+ ,468
+ ,20
+ ,60720
+ ,312161
+ ,20
+ ,1
+ ,60805
+ ,315576
+ ,278
+ ,22
+ ,61404
+ ,314922
+ ,61
+ ,9
+ ,60720
+ ,314551
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,192
+ ,2
+ ,65276
+ ,312339
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,317
+ ,7
+ ,63915
+ ,298700
+ ,738
+ ,9
+ ,60720
+ ,321376
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,368
+ ,13
+ ,61686
+ ,303230
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,2
+ ,0
+ ,60743
+ ,315487
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,53
+ ,6
+ ,60349
+ ,315793
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,94
+ ,3
+ ,61360
+ ,312887
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,24
+ ,7
+ ,59818
+ ,315637
+ ,2332
+ ,2
+ ,72680
+ ,324385
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,131
+ ,15
+ ,61808
+ ,308989
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,206
+ ,9
+ ,53110
+ ,296702
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,167
+ ,1
+ ,64245
+ ,307322
+ ,622
+ ,38
+ ,73007
+ ,304376
+ ,2328
+ ,57
+ ,82732
+ ,253588
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,365
+ ,7
+ ,54820
+ ,309560
+ ,364
+ ,26
+ ,47705
+ ,298466
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,226
+ ,13
+ ,72835
+ ,343929
+ ,307
+ ,10
+ ,58856
+ ,331955
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,188
+ ,9
+ ,77655
+ ,381180
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,138
+ ,26
+ ,69817
+ ,331420
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,125
+ ,19
+ ,60798
+ ,310201
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,282
+ ,12
+ ,62452
+ ,320016
+ ,335
+ ,23
+ ,64175
+ ,320398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1324
+ ,29
+ ,67440
+ ,291841
+ ,176
+ ,8
+ ,68136
+ ,310670
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,26
+ ,56726
+ ,313491
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,333
+ ,9
+ ,70811
+ ,331323
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,601
+ ,5
+ ,60720
+ ,319210
+ ,30
+ ,3
+ ,62045
+ ,318098
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,249
+ ,13
+ ,54323
+ ,292754
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,165
+ ,12
+ ,62841
+ ,325176
+ ,453
+ ,19
+ ,81125
+ ,365959
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,53
+ ,10
+ ,59506
+ ,302409
+ ,382
+ ,9
+ ,59365
+ ,340968
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,9
+ ,60720
+ ,315380
+ ,30
+ ,4
+ ,60798
+ ,313164
+ ,290
+ ,1
+ ,58790
+ ,301164
+ ,0
+ ,1
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,366
+ ,14
+ ,61808
+ ,344425
+ ,2
+ ,12
+ ,60735
+ ,315394
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,209
+ ,19
+ ,64016
+ ,316647
+ ,384
+ ,17
+ ,54683
+ ,309836
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,365
+ ,32
+ ,87192
+ ,346611
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,49
+ ,14
+ ,64107
+ ,322031
+ ,3
+ ,8
+ ,60761
+ ,315656
+ ,133
+ ,4
+ ,65990
+ ,339445
+ ,32
+ ,0
+ ,59988
+ ,314964
+ ,368
+ ,20
+ ,61167
+ ,297141
+ ,1
+ ,5
+ ,60719
+ ,315372
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,22
+ ,1
+ ,60722
+ ,312502
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,96
+ ,4
+ ,60379
+ ,313729
+ ,1
+ ,1
+ ,60727
+ ,315388
+ ,314
+ ,4
+ ,60720
+ ,315371
+ ,844
+ ,20
+ ,60925
+ ,296139
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,26
+ ,1
+ ,60896
+ ,313880
+ ,125
+ ,10
+ ,59734
+ ,317698
+ ,304
+ ,12
+ ,62969
+ ,295580
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,621
+ ,13
+ ,60720
+ ,308256
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,119
+ ,3
+ ,59118
+ ,303677
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1595
+ ,10
+ ,60720
+ ,319369
+ ,312
+ ,3
+ ,58598
+ ,318690
+ ,60
+ ,7
+ ,61124
+ ,314049
+ ,587
+ ,10
+ ,59595
+ ,325699
+ ,135
+ ,1
+ ,62065
+ ,314210
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,514
+ ,15
+ ,78780
+ ,322378
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1
+ ,4
+ ,60722
+ ,315398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,1763
+ ,28
+ ,61600
+ ,308336
+ ,180
+ ,9
+ ,59635
+ ,316386
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,218
+ ,7
+ ,60720
+ ,315553
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,448
+ ,7
+ ,59781
+ ,323361
+ ,227
+ ,7
+ ,76644
+ ,336639
+ ,174
+ ,3
+ ,64820
+ ,307424
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,121
+ ,11
+ ,56178
+ ,295370
+ ,607
+ ,7
+ ,60436
+ ,322340
+ ,2212
+ ,10
+ ,60720
+ ,319864
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,530
+ ,18
+ ,73433
+ ,317291
+ ,571
+ ,14
+ ,41477
+ ,280398
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,78
+ ,12
+ ,62700
+ ,317330
+ ,2489
+ ,29
+ ,67804
+ ,238125
+ ,131
+ ,3
+ ,59661
+ ,327071
+ ,923
+ ,6
+ ,58620
+ ,309038
+ ,72
+ ,3
+ ,60398
+ ,314210
+ ,572
+ ,8
+ ,58580
+ ,307930
+ ,397
+ ,10
+ ,62710
+ ,322327
+ ,450
+ ,6
+ ,59325
+ ,292136
+ ,622
+ ,8
+ ,60950
+ ,263276
+ ,694
+ ,6
+ ,68060
+ ,367655
+ ,3425
+ ,9
+ ,83620
+ ,283910
+ ,562
+ ,8
+ ,58456
+ ,283587
+ ,4917
+ ,26
+ ,52811
+ ,243650
+ ,1442
+ ,239
+ ,121173
+ ,438493
+ ,529
+ ,7
+ ,63870
+ ,296261
+ ,2126
+ ,41
+ ,21001
+ ,230621
+ ,1061
+ ,3
+ ,70415
+ ,304252
+ ,776
+ ,8
+ ,64230
+ ,333505
+ ,611
+ ,6
+ ,59190
+ ,296919
+ ,1526
+ ,21
+ ,69351
+ ,278990
+ ,592
+ ,7
+ ,64270
+ ,276898
+ ,1182
+ ,11
+ ,70694
+ ,327007
+ ,621
+ ,11
+ ,68005
+ ,317046
+ ,989
+ ,12
+ ,58930
+ ,304555
+ ,438
+ ,9
+ ,58320
+ ,298096
+ ,726
+ ,3
+ ,69980
+ ,231861
+ ,1303
+ ,57
+ ,69863
+ ,309422
+ ,7419
+ ,21
+ ,63255
+ ,286963
+ ,1164
+ ,15
+ ,57320
+ ,269753
+ ,3310
+ ,32
+ ,75230
+ ,448243
+ ,1920
+ ,11
+ ,79420
+ ,165404
+ ,965
+ ,2
+ ,73490
+ ,204325
+ ,3256
+ ,23
+ ,35250
+ ,407159
+ ,1135
+ ,20
+ ,62285
+ ,290476
+ ,1270
+ ,24
+ ,69206
+ ,275311
+ ,661
+ ,1
+ ,65920
+ ,246541
+ ,1013
+ ,1
+ ,69770
+ ,253468
+ ,2844
+ ,74
+ ,72683
+ ,240897
+ ,11528
+ ,68
+ ,-14545
+ ,-83265
+ ,6526
+ ,20
+ ,55830
+ ,-42143
+ ,2264
+ ,20
+ ,55174
+ ,272713
+ ,5109
+ ,82
+ ,67038
+ ,215362
+ ,3999
+ ,21
+ ,51252
+ ,42754
+ ,35624
+ ,244
+ ,157278
+ ,306275
+ ,9252
+ ,32
+ ,79510
+ ,253537
+ ,15236
+ ,86
+ ,77440
+ ,372631
+ ,18073
+ ,69
+ ,27284
+ ,-7170)
+ ,dim=c(4
+ ,431)
+ ,dimnames=list(c('1'
+ ,'2'
+ ,'3'
+ ,'4')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('1','2','3','4'),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 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
1 2 3 4 t
1 162556 807 213118 6282154 1
2 29790 444 81767 4321023 2
3 87550 412 153198 4111912 3
4 84738 428 -26007 223193 4
5 54660 315 126942 1491348 5
6 42634 168 157214 1629616 6
7 40949 263 129352 1398893 7
8 45187 267 234817 1926517 8
9 37704 228 60448 983660 9
10 16275 129 47818 1443586 10
11 25830 104 245546 1073089 11
12 12679 122 48020 984885 12
13 18014 393 -1710 1405225 13
14 43556 190 32648 227132 14
15 24811 280 95350 929118 15
16 6575 63 151352 1071292 16
17 7123 102 288170 638830 17
18 21950 265 114337 856956 18
19 37597 234 37884 992426 19
20 17821 277 122844 444477 20
21 12988 73 82340 857217 21
22 22330 67 79801 711969 22
23 13326 103 165548 702380 23
24 16189 290 116384 358589 24
25 7146 83 134028 297978 25
26 15824 56 63838 585715 26
27 27664 236 74996 657954 27
28 11920 73 31080 209458 28
29 8568 34 32168 786690 29
30 14416 139 49857 439798 30
31 3369 26 87161 688779 31
32 11819 70 106113 574339 32
33 6984 40 80570 741409 33
34 4519 42 102129 597793 34
35 2220 12 301670 644190 35
36 18562 211 102313 377934 36
37 10327 74 88577 640273 37
38 5336 80 112477 697458 38
39 2365 83 191778 550608 39
40 4069 131 79804 207393 40
41 8636 203 128294 301607 41
42 13718 56 96448 345783 42
43 4525 89 93811 501749 43
44 6869 88 117520 379983 44
45 4628 39 69159 387475 45
46 3689 25 101792 377305 46
47 4891 49 210568 370837 47
48 7489 149 136996 430866 48
49 4901 58 121920 469107 49
50 2284 41 76403 194493 50
51 3160 90 108094 530670 51
52 4150 136 134759 518365 52
53 7285 97 188873 491303 53
54 1134 63 146216 527021 54
55 4658 114 156608 233773 55
56 2384 77 61348 405972 56
57 3748 6 50350 652925 57
58 5371 47 87720 446211 58
59 1285 51 99489 341340 59
60 9327 85 87419 387699 60
61 5565 43 94355 493408 61
62 1528 32 60326 146494 62
63 3122 25 94670 414462 63
64 7561 77 82425 364304 64
65 2675 54 59017 355178 65
66 13253 251 90829 357760 66
67 880 15 80791 261216 67
68 2053 44 100423 397144 68
69 1424 73 131116 374943 69
70 4036 85 100269 424898 70
71 3045 49 27330 202055 71
72 5119 38 39039 378525 72
73 1431 35 106885 310768 73
74 554 9 79285 325738 74
75 1975 34 118881 394510 75
76 1765 20 77623 247060 76
77 1012 29 114768 368078 77
78 810 11 74015 236761 78
79 1280 52 69465 312378 79
80 666 13 117869 339836 80
81 1380 29 60982 347385 81
82 4677 66 90131 426280 82
83 876 33 138971 352850 83
84 814 15 39625 301881 84
85 514 15 102725 377516 85
86 5692 68 64239 357312 86
87 3642 100 90262 458343 87
88 540 13 103960 354228 88
89 2099 45 106611 308636 89
90 567 14 103345 386212 90
91 2001 36 95551 393343 91
92 2949 40 82903 378509 92
93 2253 68 63593 452469 93
94 6533 29 126910 364839 94
95 1889 43 37527 358649 95
96 3055 30 60247 376641 96
97 272 9 112995 429112 97
98 1414 22 70184 330546 98
99 2564 19 130140 403560 99
100 1383 9 73221 317892 100
101 1261 31 76114 307528 101
102 975 19 90534 235133 102
103 3366 55 108479 299243 103
104 576 8 113761 314073 104
105 1686 28 68696 368186 105
106 746 29 71561 269661 106
107 3192 48 59831 125390 107
108 2045 16 97890 510834 108
109 5702 47 101481 321896 109
110 1932 20 72954 249898 110
111 936 22 67939 408881 111
112 3437 33 48022 158492 112
113 5131 44 86111 292154 113
114 2397 13 74020 289513 114
115 1389 6 57530 378049 115
116 1503 35 56364 343466 116
117 402 8 84990 332743 117
118 2239 17 88590 442882 118
119 2234 11 77200 214215 119
120 837 21 61262 315688 120
121 10579 92 110309 375195 121
122 875 12 67000 334280 122
123 1585 112 93099 355864 123
124 1659 25 107577 480382 124
125 2647 17 62920 353058 125
126 3294 23 75832 217193 126
127 0 0 60720 315380 127
128 94 10 60793 314533 128
129 422 23 57935 318056 129
130 0 0 60720 315380 130
131 34 7 60630 314353 131
132 1558 25 55637 369448 132
133 0 0 60720 315380 133
134 43 20 60887 312846 134
135 645 4 60720 312075 135
136 316 4 60505 315009 136
137 115 10 60945 318903 137
138 5 1 60720 314887 138
139 897 4 60720 314913 139
140 0 0 60720 315380 140
141 389 8 58990 325506 141
142 0 0 60720 315380 142
143 1002 11 56750 298568 143
144 36 4 60894 315834 144
145 460 15 63346 329784 145
146 309 9 56535 312878 146
147 0 0 60720 315380 147
148 9 7 60835 314987 148
149 271 2 60720 325249 149
150 14 0 61016 315877 150
151 520 7 58650 291650 151
152 1766 46 60438 305959 152
153 0 5 60720 315380 153
154 458 7 58625 297765 154
155 20 2 60938 315245 155
156 0 0 60720 315380 156
157 0 0 60720 315380 157
158 98 2 61490 315236 158
159 405 5 60845 336425 159
160 0 0 60720 315380 160
161 0 0 60720 315380 161
162 0 0 60720 315380 162
163 0 0 60720 315380 163
164 483 7 60830 306268 164
165 454 24 63261 302187 165
166 47 1 60720 314882 166
167 0 0 60720 315380 167
168 757 18 45689 382712 168
169 4655 55 60720 341570 169
170 0 0 60720 315380 170
171 0 0 60720 315380 171
172 36 3 61564 312412 172
173 0 0 60720 315380 173
174 203 9 61938 309596 174
175 0 0 60720 315380 175
176 126 8 60951 315547 176
177 400 113 60720 313267 177
178 71 0 60745 316176 178
179 0 0 60720 315380 179
180 0 0 60720 315380 180
181 972 19 71642 359335 181
182 531 11 71641 330068 182
183 2461 25 55792 314289 183
184 378 16 71873 297413 184
185 23 5 62555 314806 185
186 638 11 60370 333210 186
187 2300 23 64873 352108 187
188 149 6 62041 313332 188
189 226 5 65745 291787 189
190 0 0 60720 315380 190
191 275 7 59500 318745 191
192 0 0 60720 315380 192
193 141 7 61630 315366 193
194 0 0 60720 315380 194
195 28 3 60890 315688 195
196 0 0 60720 315380 196
197 4980 89 113521 409642 197
198 0 0 60720 315380 198
199 0 0 60720 315380 199
200 472 19 80045 269587 200
201 0 0 60720 315380 201
202 0 0 60720 315380 202
203 0 0 60720 315380 203
204 203 12 50804 300962 204
205 496 12 87390 325479 205
206 10 5 61656 316155 206
207 63 2 65688 318574 207
208 0 0 60720 315380 208
209 1136 26 48522 343613 209
210 265 3 60720 306948 210
211 0 0 60720 315380 211
212 0 0 60720 315380 212
213 267 11 57640 330059 213
214 474 10 61977 288985 214
215 534 5 62620 304485 215
216 0 2 60720 315380 216
217 15 6 60831 315688 217
218 397 7 60646 317736 218
219 0 2 60720 315380 219
220 1866 28 56225 322331 220
221 288 3 60510 296656 221
222 0 0 60720 315380 222
223 3 1 60698 315354 223
224 468 20 60720 312161 224
225 20 1 60805 315576 225
226 278 22 61404 314922 226
227 61 9 60720 314551 227
228 0 0 60720 315380 228
229 192 2 65276 312339 229
230 0 0 60720 315380 230
231 317 7 63915 298700 231
232 738 9 60720 321376 232
233 0 0 60720 315380 233
234 368 13 61686 303230 234
235 0 0 60720 315380 235
236 2 0 60743 315487 236
237 0 0 60720 315380 237
238 53 6 60349 315793 238
239 0 0 60720 315380 239
240 0 0 60720 315380 240
241 0 0 60720 315380 241
242 94 3 61360 312887 242
243 0 0 60720 315380 243
244 24 7 59818 315637 244
245 2332 2 72680 324385 245
246 0 0 60720 315380 246
247 0 0 60720 315380 247
248 131 15 61808 308989 248
249 0 0 60720 315380 249
250 0 0 60720 315380 250
251 206 9 53110 296702 251
252 0 0 60720 315380 252
253 167 1 64245 307322 253
254 622 38 73007 304376 254
255 2328 57 82732 253588 255
256 0 0 60720 315380 256
257 365 7 54820 309560 257
258 364 26 47705 298466 258
259 0 0 60720 315380 259
260 0 0 60720 315380 260
261 0 0 60720 315380 261
262 0 0 60720 315380 262
263 226 13 72835 343929 263
264 307 10 58856 331955 264
265 0 0 60720 315380 265
266 0 0 60720 315380 266
267 0 0 60720 315380 267
268 188 9 77655 381180 268
269 0 0 60720 315380 269
270 138 26 69817 331420 270
271 0 0 60720 315380 271
272 0 0 60720 315380 272
273 0 0 60720 315380 273
274 125 19 60798 310201 274
275 0 0 60720 315380 275
276 282 12 62452 320016 276
277 335 23 64175 320398 277
278 0 0 60720 315380 278
279 1324 29 67440 291841 279
280 176 8 68136 310670 280
281 0 0 60720 315380 281
282 0 0 60720 315380 282
283 249 26 56726 313491 283
284 0 0 60720 315380 284
285 333 9 70811 331323 285
286 0 0 60720 315380 286
287 601 5 60720 319210 287
288 30 3 62045 318098 288
289 0 0 60720 315380 289
290 249 13 54323 292754 290
291 0 0 60720 315380 291
292 165 12 62841 325176 292
293 453 19 81125 365959 293
294 0 0 60720 315380 294
295 53 10 59506 302409 295
296 382 9 59365 340968 296
297 0 0 60720 315380 297
298 0 0 60720 315380 298
299 0 0 60720 315380 299
300 0 9 60720 315380 300
301 30 4 60798 313164 301
302 290 1 58790 301164 302
303 0 1 60720 315380 303
304 0 0 60720 315380 304
305 366 14 61808 344425 305
306 2 12 60735 315394 306
307 0 0 60720 315380 307
308 209 19 64016 316647 308
309 384 17 54683 309836 309
310 0 0 60720 315380 310
311 0 0 60720 315380 311
312 365 32 87192 346611 312
313 0 0 60720 315380 313
314 49 14 64107 322031 314
315 3 8 60761 315656 315
316 133 4 65990 339445 316
317 32 0 59988 314964 317
318 368 20 61167 297141 318
319 1 5 60719 315372 319
320 0 0 60720 315380 320
321 0 0 60720 315380 321
322 0 0 60720 315380 322
323 0 0 60720 315380 323
324 0 0 60720 315380 324
325 0 0 60720 315380 325
326 22 1 60722 312502 326
327 0 0 60720 315380 327
328 0 0 60720 315380 328
329 0 0 60720 315380 329
330 0 0 60720 315380 330
331 0 0 60720 315380 331
332 0 0 60720 315380 332
333 0 0 60720 315380 333
334 96 4 60379 313729 334
335 1 1 60727 315388 335
336 314 4 60720 315371 336
337 844 20 60925 296139 337
338 0 0 60720 315380 338
339 26 1 60896 313880 339
340 125 10 59734 317698 340
341 304 12 62969 295580 341
342 0 0 60720 315380 342
343 0 0 60720 315380 343
344 0 0 60720 315380 344
345 621 13 60720 308256 345
346 0 0 60720 315380 346
347 119 3 59118 303677 347
348 0 0 60720 315380 348
349 0 0 60720 315380 349
350 1595 10 60720 319369 350
351 312 3 58598 318690 351
352 60 7 61124 314049 352
353 587 10 59595 325699 353
354 135 1 62065 314210 354
355 0 0 60720 315380 355
356 0 0 60720 315380 356
357 514 15 78780 322378 357
358 0 0 60720 315380 358
359 0 0 60720 315380 359
360 0 0 60720 315380 360
361 1 4 60722 315398 361
362 0 0 60720 315380 362
363 0 0 60720 315380 363
364 1763 28 61600 308336 364
365 180 9 59635 316386 365
366 0 0 60720 315380 366
367 0 0 60720 315380 367
368 0 0 60720 315380 368
369 0 0 60720 315380 369
370 218 7 60720 315553 370
371 0 0 60720 315380 371
372 448 7 59781 323361 372
373 227 7 76644 336639 373
374 174 3 64820 307424 374
375 0 0 60720 315380 375
376 0 0 60720 315380 376
377 121 11 56178 295370 377
378 607 7 60436 322340 378
379 2212 10 60720 319864 379
380 0 0 60720 315380 380
381 0 0 60720 315380 381
382 530 18 73433 317291 382
383 571 14 41477 280398 383
384 0 0 60720 315380 384
385 78 12 62700 317330 385
386 2489 29 67804 238125 386
387 131 3 59661 327071 387
388 923 6 58620 309038 388
389 72 3 60398 314210 389
390 572 8 58580 307930 390
391 397 10 62710 322327 391
392 450 6 59325 292136 392
393 622 8 60950 263276 393
394 694 6 68060 367655 394
395 3425 9 83620 283910 395
396 562 8 58456 283587 396
397 4917 26 52811 243650 397
398 1442 239 121173 438493 398
399 529 7 63870 296261 399
400 2126 41 21001 230621 400
401 1061 3 70415 304252 401
402 776 8 64230 333505 402
403 611 6 59190 296919 403
404 1526 21 69351 278990 404
405 592 7 64270 276898 405
406 1182 11 70694 327007 406
407 621 11 68005 317046 407
408 989 12 58930 304555 408
409 438 9 58320 298096 409
410 726 3 69980 231861 410
411 1303 57 69863 309422 411
412 7419 21 63255 286963 412
413 1164 15 57320 269753 413
414 3310 32 75230 448243 414
415 1920 11 79420 165404 415
416 965 2 73490 204325 416
417 3256 23 35250 407159 417
418 1135 20 62285 290476 418
419 1270 24 69206 275311 419
420 661 1 65920 246541 420
421 1013 1 69770 253468 421
422 2844 74 72683 240897 422
423 11528 68 -14545 -83265 423
424 6526 20 55830 -42143 424
425 2264 20 55174 272713 425
426 5109 82 67038 215362 426
427 3999 21 51252 42754 427
428 35624 244 157278 306275 428
429 9252 32 79510 253537 429
430 15236 86 77440 372631 430
431 18073 69 27284 -7170 431
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) `2` `3` `4` t
-3.391e+03 1.059e+02 -1.292e-02 8.837e-03 3.862e+00
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-50962.3 -641.5 114.4 639.6 40491.7
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -3.391e+03 9.542e+02 -3.554 0.000422 ***
`2` 1.059e+02 5.251e+00 20.163 < 2e-16 ***
`3` -1.292e-02 8.774e-03 -1.473 0.141520
`4` 8.837e-03 8.730e-04 10.122 < 2e-16 ***
t 3.862e+00 2.302e+00 1.678 0.094117 .
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 5093 on 426 degrees of freedom
Multiple R-squared: 0.816, Adjusted R-squared: 0.8143
F-statistic: 472.4 on 4 and 426 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 1.710754e-25 8.553768e-26
[2,] 1.0000000 6.654162e-36 3.327081e-36
[3,] 1.0000000 8.692591e-38 4.346295e-38
[4,] 1.0000000 2.085906e-39 1.042953e-39
[5,] 1.0000000 1.362213e-38 6.811064e-39
[6,] 1.0000000 1.397681e-55 6.988406e-56
[7,] 1.0000000 2.295201e-78 1.147600e-78
[8,] 1.0000000 3.317768e-79 1.658884e-79
[9,] 1.0000000 2.276276e-79 1.138138e-79
[10,] 1.0000000 2.036840e-81 1.018420e-81
[11,] 1.0000000 1.029413e-80 5.147066e-81
[12,] 1.0000000 4.014202e-92 2.007101e-92
[13,] 1.0000000 4.597675e-93 2.298838e-93
[14,] 1.0000000 2.937780e-98 1.468890e-98
[15,] 1.0000000 6.775108e-112 3.387554e-112
[16,] 1.0000000 2.754994e-112 1.377497e-112
[17,] 1.0000000 8.046763e-114 4.023382e-114
[18,] 1.0000000 1.676688e-113 8.383440e-114
[19,] 1.0000000 1.101032e-120 5.505159e-121
[20,] 1.0000000 3.547343e-124 1.773672e-124
[21,] 1.0000000 1.362573e-126 6.812867e-127
[22,] 1.0000000 1.510882e-128 7.554411e-129
[23,] 1.0000000 5.401048e-129 2.700524e-129
[24,] 1.0000000 1.506448e-128 7.532241e-129
[25,] 1.0000000 9.519440e-131 4.759720e-131
[26,] 1.0000000 8.060314e-132 4.030157e-132
[27,] 1.0000000 2.024811e-131 1.012405e-131
[28,] 1.0000000 1.330001e-130 6.650007e-131
[29,] 1.0000000 6.628154e-131 3.314077e-131
[30,] 1.0000000 5.140435e-133 2.570218e-133
[31,] 1.0000000 1.510334e-132 7.551670e-133
[32,] 1.0000000 4.870180e-132 2.435090e-132
[33,] 1.0000000 1.850058e-132 9.250288e-133
[34,] 1.0000000 2.272215e-133 1.136108e-133
[35,] 1.0000000 1.717043e-139 8.585217e-140
[36,] 1.0000000 1.353526e-138 6.767629e-139
[37,] 1.0000000 6.827691e-138 3.413845e-138
[38,] 1.0000000 7.039825e-138 3.519912e-138
[39,] 1.0000000 1.011027e-137 5.055135e-138
[40,] 1.0000000 5.010244e-137 2.505122e-137
[41,] 1.0000000 1.978649e-136 9.893243e-137
[42,] 1.0000000 6.235696e-136 3.117848e-136
[43,] 1.0000000 2.277196e-135 1.138598e-135
[44,] 1.0000000 1.418098e-134 7.090489e-135
[45,] 1.0000000 1.298094e-134 6.490471e-135
[46,] 1.0000000 6.980535e-134 3.490267e-134
[47,] 1.0000000 3.263186e-133 1.631593e-133
[48,] 1.0000000 1.524396e-133 7.621978e-134
[49,] 1.0000000 7.387874e-133 3.693937e-133
[50,] 1.0000000 8.626677e-136 4.313338e-136
[51,] 1.0000000 4.448729e-136 2.224365e-136
[52,] 1.0000000 1.418381e-135 7.091904e-136
[53,] 1.0000000 3.250002e-136 1.625001e-136
[54,] 1.0000000 5.029296e-137 2.514648e-137
[55,] 1.0000000 8.907800e-137 4.453900e-137
[56,] 1.0000000 1.344697e-136 6.723486e-137
[57,] 1.0000000 1.513833e-136 7.569165e-137
[58,] 1.0000000 9.374310e-136 4.687155e-136
[59,] 1.0000000 6.039303e-136 3.019651e-136
[60,] 1.0000000 1.949325e-135 9.746625e-136
[61,] 1.0000000 1.107673e-134 5.538366e-135
[62,] 1.0000000 9.438040e-135 4.719020e-135
[63,] 1.0000000 5.808994e-134 2.904497e-134
[64,] 1.0000000 2.428833e-133 1.214416e-133
[65,] 1.0000000 4.721881e-134 2.360941e-134
[66,] 1.0000000 1.519807e-133 7.599036e-134
[67,] 1.0000000 6.116898e-133 3.058449e-133
[68,] 1.0000000 3.082550e-132 1.541275e-132
[69,] 1.0000000 1.095411e-131 5.477053e-132
[70,] 1.0000000 4.732582e-131 2.366291e-131
[71,] 1.0000000 1.566722e-130 7.833608e-131
[72,] 1.0000000 5.313058e-130 2.656529e-130
[73,] 1.0000000 1.974743e-129 9.873715e-130
[74,] 1.0000000 1.055690e-128 5.278449e-129
[75,] 1.0000000 4.543856e-128 2.271928e-128
[76,] 1.0000000 7.975849e-128 3.987925e-128
[77,] 1.0000000 3.550344e-127 1.775172e-127
[78,] 1.0000000 1.712311e-126 8.561556e-127
[79,] 1.0000000 4.691571e-126 2.345785e-126
[80,] 1.0000000 1.914981e-125 9.574905e-126
[81,] 1.0000000 8.063458e-125 4.031729e-125
[82,] 1.0000000 2.121884e-124 1.060942e-124
[83,] 1.0000000 1.007424e-123 5.037120e-124
[84,] 1.0000000 5.417755e-123 2.708877e-123
[85,] 1.0000000 2.505281e-122 1.252641e-122
[86,] 1.0000000 1.377794e-121 6.888969e-122
[87,] 1.0000000 5.908783e-122 2.954392e-122
[88,] 1.0000000 2.935640e-121 1.467820e-121
[89,] 1.0000000 5.651767e-121 2.825884e-121
[90,] 1.0000000 2.819552e-120 1.409776e-120
[91,] 1.0000000 1.362380e-119 6.811902e-120
[92,] 1.0000000 5.389623e-119 2.694812e-119
[93,] 1.0000000 2.062973e-118 1.031486e-118
[94,] 1.0000000 9.547514e-118 4.773757e-118
[95,] 1.0000000 2.200585e-117 1.100293e-117
[96,] 1.0000000 7.042896e-117 3.521448e-117
[97,] 1.0000000 2.024851e-116 1.012425e-116
[98,] 1.0000000 9.567836e-116 4.783918e-116
[99,] 1.0000000 2.879035e-115 1.439517e-115
[100,] 1.0000000 4.174271e-115 2.087136e-115
[101,] 1.0000000 4.322663e-115 2.161331e-115
[102,] 1.0000000 1.381055e-114 6.905274e-115
[103,] 1.0000000 5.700019e-114 2.850010e-114
[104,] 1.0000000 2.343260e-113 1.171630e-113
[105,] 1.0000000 8.902223e-113 4.451112e-113
[106,] 1.0000000 3.283000e-112 1.641500e-112
[107,] 1.0000000 1.236095e-111 6.180475e-112
[108,] 1.0000000 2.096969e-111 1.048484e-111
[109,] 1.0000000 1.075105e-110 5.375525e-111
[110,] 1.0000000 5.101800e-110 2.550900e-110
[111,] 1.0000000 7.903806e-110 3.951903e-110
[112,] 1.0000000 2.840990e-109 1.420495e-109
[113,] 1.0000000 1.439173e-108 7.195863e-109
[114,] 1.0000000 4.721834e-109 2.360917e-109
[115,] 1.0000000 2.178959e-108 1.089479e-108
[116,] 1.0000000 4.471762e-110 2.235881e-110
[117,] 1.0000000 1.348201e-109 6.741004e-110
[118,] 1.0000000 2.306906e-109 1.153453e-109
[119,] 1.0000000 9.255384e-109 4.627692e-109
[120,] 1.0000000 4.318849e-108 2.159424e-108
[121,] 1.0000000 2.146713e-107 1.073356e-107
[122,] 1.0000000 1.042782e-106 5.213911e-107
[123,] 1.0000000 4.818900e-106 2.409450e-106
[124,] 1.0000000 2.351100e-105 1.175550e-105
[125,] 1.0000000 8.072534e-105 4.036267e-105
[126,] 1.0000000 3.690608e-104 1.845304e-104
[127,] 1.0000000 1.617037e-103 8.085186e-104
[128,] 1.0000000 6.890352e-103 3.445176e-103
[129,] 1.0000000 3.119316e-102 1.559658e-102
[130,] 1.0000000 1.497617e-101 7.488087e-102
[131,] 1.0000000 6.824097e-101 3.412049e-101
[132,] 1.0000000 2.625529e-100 1.312765e-100
[133,] 1.0000000 1.174618e-99 5.873091e-100
[134,] 1.0000000 5.248782e-99 2.624391e-99
[135,] 1.0000000 2.330251e-98 1.165125e-98
[136,] 1.0000000 1.022586e-97 5.112931e-98
[137,] 1.0000000 4.668558e-97 2.334279e-97
[138,] 1.0000000 2.167809e-96 1.083904e-96
[139,] 1.0000000 9.843627e-96 4.921814e-96
[140,] 1.0000000 4.281870e-95 2.140935e-95
[141,] 1.0000000 1.949910e-94 9.749549e-95
[142,] 1.0000000 7.949215e-94 3.974607e-94
[143,] 1.0000000 3.408587e-93 1.704294e-93
[144,] 1.0000000 1.502944e-92 7.514722e-93
[145,] 1.0000000 5.721286e-92 2.860643e-92
[146,] 1.0000000 2.530528e-91 1.265264e-91
[147,] 1.0000000 1.101835e-90 5.509175e-91
[148,] 1.0000000 4.729651e-90 2.364826e-90
[149,] 1.0000000 1.984243e-89 9.921213e-90
[150,] 1.0000000 8.289693e-89 4.144847e-89
[151,] 1.0000000 3.485405e-88 1.742702e-88
[152,] 1.0000000 1.325829e-87 6.629146e-88
[153,] 1.0000000 5.463802e-87 2.731901e-87
[154,] 1.0000000 2.241756e-86 1.120878e-86
[155,] 1.0000000 9.156854e-86 4.578427e-86
[156,] 1.0000000 3.723475e-85 1.861737e-85
[157,] 1.0000000 1.541536e-84 7.707679e-85
[158,] 1.0000000 5.672045e-84 2.836022e-84
[159,] 1.0000000 2.291072e-83 1.145536e-83
[160,] 1.0000000 9.167074e-83 4.583537e-83
[161,] 1.0000000 2.262462e-82 1.131231e-82
[162,] 1.0000000 5.245471e-82 2.622735e-82
[163,] 1.0000000 2.064060e-81 1.032030e-81
[164,] 1.0000000 8.084259e-81 4.042130e-81
[165,] 1.0000000 3.254799e-80 1.627400e-80
[166,] 1.0000000 1.263225e-79 6.316126e-80
[167,] 1.0000000 5.122530e-79 2.561265e-79
[168,] 1.0000000 1.970578e-78 9.852891e-79
[169,] 1.0000000 7.888558e-78 3.944279e-78
[170,] 1.0000000 2.741897e-80 1.370949e-80
[171,] 1.0000000 1.065276e-79 5.326378e-80
[172,] 1.0000000 4.174439e-79 2.087219e-79
[173,] 1.0000000 1.628537e-78 8.142683e-79
[174,] 1.0000000 6.279471e-78 3.139735e-78
[175,] 1.0000000 2.517428e-77 1.258714e-77
[176,] 1.0000000 8.099773e-77 4.049887e-77
[177,] 1.0000000 2.872214e-76 1.436107e-76
[178,] 1.0000000 1.138852e-75 5.694262e-76
[179,] 1.0000000 4.231372e-75 2.115686e-75
[180,] 1.0000000 1.124359e-74 5.621796e-75
[181,] 1.0000000 4.394526e-74 2.197263e-74
[182,] 1.0000000 1.690443e-73 8.452215e-74
[183,] 1.0000000 6.311345e-73 3.155672e-73
[184,] 1.0000000 2.393711e-72 1.196855e-72
[185,] 1.0000000 8.853066e-72 4.426533e-72
[186,] 1.0000000 3.390136e-71 1.695068e-71
[187,] 1.0000000 1.242557e-70 6.212784e-71
[188,] 1.0000000 4.643715e-70 2.321857e-70
[189,] 1.0000000 1.686197e-69 8.430983e-70
[190,] 1.0000000 5.220293e-69 2.610147e-69
[191,] 1.0000000 1.879528e-68 9.397641e-69
[192,] 1.0000000 6.735223e-68 3.367612e-68
[193,] 1.0000000 1.600054e-67 8.000270e-68
[194,] 1.0000000 5.723029e-67 2.861514e-67
[195,] 1.0000000 2.037355e-66 1.018678e-66
[196,] 1.0000000 7.218482e-66 3.609241e-66
[197,] 1.0000000 2.643477e-65 1.321739e-65
[198,] 1.0000000 9.223163e-65 4.611582e-65
[199,] 1.0000000 3.337883e-64 1.668942e-64
[200,] 1.0000000 1.187738e-63 5.938689e-64
[201,] 1.0000000 4.124349e-63 2.062175e-63
[202,] 1.0000000 1.341353e-62 6.706763e-63
[203,] 1.0000000 4.665119e-62 2.332559e-62
[204,] 1.0000000 1.595646e-61 7.978229e-62
[205,] 1.0000000 5.431371e-61 2.715686e-61
[206,] 1.0000000 1.866014e-60 9.330068e-61
[207,] 1.0000000 6.464081e-60 3.232040e-60
[208,] 1.0000000 2.193570e-59 1.096785e-59
[209,] 1.0000000 7.459248e-59 3.729624e-59
[210,] 1.0000000 2.569799e-58 1.284900e-58
[211,] 1.0000000 8.586124e-58 4.293062e-58
[212,] 1.0000000 2.878696e-57 1.439348e-57
[213,] 1.0000000 9.164902e-57 4.582451e-57
[214,] 1.0000000 3.064730e-56 1.532365e-56
[215,] 1.0000000 9.957318e-56 4.978659e-56
[216,] 1.0000000 3.247031e-55 1.623516e-55
[217,] 1.0000000 1.057548e-54 5.287742e-55
[218,] 1.0000000 3.413693e-54 1.706847e-54
[219,] 1.0000000 1.053533e-53 5.267665e-54
[220,] 1.0000000 3.473906e-53 1.736953e-53
[221,] 1.0000000 1.101714e-52 5.508569e-53
[222,] 1.0000000 3.535316e-52 1.767658e-52
[223,] 1.0000000 1.110443e-51 5.552215e-52
[224,] 1.0000000 3.589113e-51 1.794556e-51
[225,] 1.0000000 1.085226e-50 5.426130e-51
[226,] 1.0000000 3.360505e-50 1.680252e-50
[227,] 1.0000000 1.066151e-49 5.330754e-50
[228,] 1.0000000 3.273547e-49 1.636774e-49
[229,] 1.0000000 9.996456e-49 4.998228e-49
[230,] 1.0000000 3.038266e-48 1.519133e-48
[231,] 1.0000000 9.472478e-48 4.736239e-48
[232,] 1.0000000 2.851242e-47 1.425621e-47
[233,] 1.0000000 8.538227e-47 4.269114e-47
[234,] 1.0000000 2.543631e-46 1.271816e-46
[235,] 1.0000000 7.691918e-46 3.845959e-46
[236,] 1.0000000 2.268512e-45 1.134256e-45
[237,] 1.0000000 6.890957e-45 3.445478e-45
[238,] 1.0000000 1.072591e-44 5.362953e-45
[239,] 1.0000000 3.112327e-44 1.556163e-44
[240,] 1.0000000 8.982371e-44 4.491186e-44
[241,] 1.0000000 2.645354e-43 1.322677e-43
[242,] 1.0000000 7.573813e-43 3.786906e-43
[243,] 1.0000000 2.156693e-42 1.078347e-42
[244,] 1.0000000 6.375871e-42 3.187935e-42
[245,] 1.0000000 1.797981e-41 8.989906e-42
[246,] 1.0000000 5.143568e-41 2.571784e-41
[247,] 1.0000000 9.942264e-41 4.971132e-41
[248,] 1.0000000 1.125724e-40 5.628620e-41
[249,] 1.0000000 3.194452e-40 1.597226e-40
[250,] 1.0000000 9.047241e-40 4.523621e-40
[251,] 1.0000000 2.433076e-39 1.216538e-39
[252,] 1.0000000 6.817395e-39 3.408698e-39
[253,] 1.0000000 1.900086e-38 9.500428e-39
[254,] 1.0000000 5.267541e-38 2.633771e-38
[255,] 1.0000000 1.452482e-37 7.262408e-38
[256,] 1.0000000 4.093246e-37 2.046623e-37
[257,] 1.0000000 1.116986e-36 5.584932e-37
[258,] 1.0000000 3.033793e-36 1.516896e-36
[259,] 1.0000000 8.195223e-36 4.097611e-36
[260,] 1.0000000 2.201707e-35 1.100853e-35
[261,] 1.0000000 5.698574e-35 2.849287e-35
[262,] 1.0000000 1.513667e-34 7.568335e-35
[263,] 1.0000000 3.643535e-34 1.821768e-34
[264,] 1.0000000 9.620643e-34 4.810322e-34
[265,] 1.0000000 2.526287e-33 1.263144e-33
[266,] 1.0000000 6.596980e-33 3.298490e-33
[267,] 1.0000000 1.685226e-32 8.426128e-33
[268,] 1.0000000 4.367518e-32 2.183759e-32
[269,] 1.0000000 1.160001e-31 5.800007e-32
[270,] 1.0000000 2.903955e-31 1.451978e-31
[271,] 1.0000000 7.435475e-31 3.717738e-31
[272,] 1.0000000 1.823400e-30 9.116998e-31
[273,] 1.0000000 4.745341e-30 2.372670e-30
[274,] 1.0000000 1.201885e-29 6.009427e-30
[275,] 1.0000000 3.027155e-29 1.513578e-29
[276,] 1.0000000 7.127631e-29 3.563816e-29
[277,] 1.0000000 1.782279e-28 8.911397e-29
[278,] 1.0000000 4.528513e-28 2.264256e-28
[279,] 1.0000000 1.120443e-27 5.602216e-28
[280,] 1.0000000 2.665753e-27 1.332877e-27
[281,] 1.0000000 6.609466e-27 3.304733e-27
[282,] 1.0000000 1.607357e-26 8.036783e-27
[283,] 1.0000000 3.959537e-26 1.979768e-26
[284,] 1.0000000 9.538666e-26 4.769333e-26
[285,] 1.0000000 2.337968e-25 1.168984e-25
[286,] 1.0000000 5.662636e-25 2.831318e-25
[287,] 1.0000000 1.343999e-24 6.719993e-25
[288,] 1.0000000 3.228810e-24 1.614405e-24
[289,] 1.0000000 7.373812e-24 3.686906e-24
[290,] 1.0000000 1.722396e-23 8.611978e-24
[291,] 1.0000000 3.999568e-23 1.999784e-23
[292,] 1.0000000 9.232427e-23 4.616213e-23
[293,] 1.0000000 2.172455e-22 1.086228e-22
[294,] 1.0000000 5.057659e-22 2.528830e-22
[295,] 1.0000000 1.146722e-21 5.733611e-22
[296,] 1.0000000 2.604117e-21 1.302058e-21
[297,] 1.0000000 5.844932e-21 2.922466e-21
[298,] 1.0000000 1.309131e-20 6.545655e-21
[299,] 1.0000000 2.968757e-20 1.484379e-20
[300,] 1.0000000 6.556118e-20 3.278059e-20
[301,] 1.0000000 1.436613e-19 7.183065e-20
[302,] 1.0000000 3.206352e-19 1.603176e-19
[303,] 1.0000000 6.984208e-19 3.492104e-19
[304,] 1.0000000 1.511904e-18 7.559522e-19
[305,] 1.0000000 2.762060e-18 1.381030e-18
[306,] 1.0000000 5.950144e-18 2.975072e-18
[307,] 1.0000000 1.284300e-17 6.421498e-18
[308,] 1.0000000 2.786895e-17 1.393448e-17
[309,] 1.0000000 5.867596e-17 2.933798e-17
[310,] 1.0000000 1.232109e-16 6.160544e-17
[311,] 1.0000000 2.510663e-16 1.255332e-16
[312,] 1.0000000 5.317880e-16 2.658940e-16
[313,] 1.0000000 1.103953e-15 5.519763e-16
[314,] 1.0000000 2.277122e-15 1.138561e-15
[315,] 1.0000000 4.666874e-15 2.333437e-15
[316,] 1.0000000 9.502766e-15 4.751383e-15
[317,] 1.0000000 1.922372e-14 9.611860e-15
[318,] 1.0000000 3.863381e-14 1.931690e-14
[319,] 1.0000000 7.758306e-14 3.879153e-14
[320,] 1.0000000 1.538832e-13 7.694161e-14
[321,] 1.0000000 3.031764e-13 1.515882e-13
[322,] 1.0000000 5.932764e-13 2.966382e-13
[323,] 1.0000000 1.153064e-12 5.765318e-13
[324,] 1.0000000 2.225671e-12 1.112835e-12
[325,] 1.0000000 4.266349e-12 2.133174e-12
[326,] 1.0000000 8.121112e-12 4.060556e-12
[327,] 1.0000000 1.552788e-11 7.763938e-12
[328,] 1.0000000 2.928059e-11 1.464029e-11
[329,] 1.0000000 5.450244e-11 2.725122e-11
[330,] 1.0000000 1.029919e-10 5.149593e-11
[331,] 1.0000000 1.896899e-10 9.484496e-11
[332,] 1.0000000 3.484734e-10 1.742367e-10
[333,] 1.0000000 6.451233e-10 3.225617e-10
[334,] 1.0000000 1.186619e-09 5.933094e-10
[335,] 1.0000000 2.129895e-09 1.064947e-09
[336,] 1.0000000 3.794605e-09 1.897302e-09
[337,] 1.0000000 6.709788e-09 3.354894e-09
[338,] 1.0000000 1.198716e-08 5.993581e-09
[339,] 1.0000000 2.089391e-08 1.044696e-08
[340,] 1.0000000 3.658100e-08 1.829050e-08
[341,] 1.0000000 6.281278e-08 3.140639e-08
[342,] 0.9999999 1.070146e-07 5.350730e-08
[343,] 0.9999999 1.658959e-07 8.294796e-08
[344,] 0.9999999 2.732675e-07 1.366338e-07
[345,] 0.9999998 4.636839e-07 2.318420e-07
[346,] 0.9999996 7.585449e-07 3.792725e-07
[347,] 0.9999994 1.240475e-06 6.202377e-07
[348,] 0.9999990 2.011056e-06 1.005528e-06
[349,] 0.9999984 3.233106e-06 1.616553e-06
[350,] 0.9999973 5.324713e-06 2.662357e-06
[351,] 0.9999958 8.442116e-06 4.221058e-06
[352,] 0.9999934 1.327065e-05 6.635327e-06
[353,] 0.9999897 2.068172e-05 1.034086e-05
[354,] 0.9999838 3.242124e-05 1.621062e-05
[355,] 0.9999752 4.968367e-05 2.484183e-05
[356,] 0.9999623 7.546769e-05 3.773385e-05
[357,] 0.9999423 1.154014e-04 5.770070e-05
[358,] 0.9999122 1.755309e-04 8.776544e-05
[359,] 0.9998700 2.599641e-04 1.299820e-04
[360,] 0.9998092 3.815281e-04 1.907640e-04
[361,] 0.9997226 5.548335e-04 2.774168e-04
[362,] 0.9996003 7.994517e-04 3.997258e-04
[363,] 0.9994216 1.156854e-03 5.784270e-04
[364,] 0.9991814 1.637119e-03 8.185593e-04
[365,] 0.9988620 2.275975e-03 1.137987e-03
[366,] 0.9983834 3.233216e-03 1.616608e-03
[367,] 0.9977395 4.521025e-03 2.260513e-03
[368,] 0.9969152 6.169679e-03 3.084839e-03
[369,] 0.9958304 8.339131e-03 4.169565e-03
[370,] 0.9942328 1.153432e-02 5.767162e-03
[371,] 0.9924475 1.510494e-02 7.552469e-03
[372,] 0.9912932 1.741369e-02 8.706843e-03
[373,] 0.9886896 2.262072e-02 1.131036e-02
[374,] 0.9854522 2.909557e-02 1.454778e-02
[375,] 0.9806191 3.876187e-02 1.938093e-02
[376,] 0.9751826 4.963481e-02 2.481741e-02
[377,] 0.9688839 6.223226e-02 3.111613e-02
[378,] 0.9600740 7.985192e-02 3.992596e-02
[379,] 0.9491476 1.017048e-01 5.085242e-02
[380,] 0.9382556 1.234888e-01 6.174442e-02
[381,] 0.9273810 1.452379e-01 7.261897e-02
[382,] 0.9124950 1.750101e-01 8.750504e-02
[383,] 0.8963113 2.073773e-01 1.036887e-01
[384,] 0.8768588 2.462825e-01 1.231412e-01
[385,] 0.8552740 2.894520e-01 1.447260e-01
[386,] 0.8289075 3.421850e-01 1.710925e-01
[387,] 0.8113656 3.772689e-01 1.886344e-01
[388,] 0.8109618 3.780763e-01 1.890382e-01
[389,] 0.7863540 4.272920e-01 2.136460e-01
[390,] 0.8113666 3.772667e-01 1.886334e-01
[391,] 0.9994130 1.173947e-03 5.869737e-04
[392,] 0.9990277 1.944681e-03 9.723407e-04
[393,] 0.9984607 3.078524e-03 1.539262e-03
[394,] 0.9977372 4.525584e-03 2.262792e-03
[395,] 0.9964909 7.018215e-03 3.509107e-03
[396,] 0.9946008 1.079846e-02 5.399228e-03
[397,] 0.9912203 1.755939e-02 8.779695e-03
[398,] 0.9865361 2.692778e-02 1.346389e-02
[399,] 0.9804174 3.916510e-02 1.958255e-02
[400,] 0.9707785 5.844304e-02 2.922152e-02
[401,] 0.9578760 8.424805e-02 4.212403e-02
[402,] 0.9396744 1.206512e-01 6.032562e-02
[403,] 0.9182080 1.635840e-01 8.179199e-02
[404,] 0.9167186 1.665627e-01 8.328136e-02
[405,] 0.9485738 1.028524e-01 5.142621e-02
[406,] 0.9237498 1.525005e-01 7.625023e-02
[407,] 0.9001299 1.997401e-01 9.987007e-02
[408,] 0.8641954 2.716091e-01 1.358046e-01
[409,] 0.8337720 3.324559e-01 1.662280e-01
[410,] 0.8383154 3.233693e-01 1.616846e-01
[411,] 0.7887647 4.224706e-01 2.112353e-01
[412,] 0.7131482 5.737036e-01 2.868518e-01
[413,] 0.6877046 6.245909e-01 3.122954e-01
[414,] 0.7737841 4.524318e-01 2.262159e-01
[415,] 0.6682213 6.635574e-01 3.317787e-01
[416,] 0.6291731 7.416538e-01 3.708269e-01
> postscript(file="/var/www/html/rcomp/tmp/1ey931291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/2ey931291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/3pp861291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4pp861291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5pp861291213708.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
27742.767049 -50962.285449 12952.860263 40491.729428 13143.453408
6 7 8 9 10
15846.235052 5778.285078 6289.315915 9010.028560 -6169.038316
11 12 13 14 15
11858.409601 -4975.490298 -32693.119223 25192.057311 -8478.469344
16 17 18 19 20
-4276.584071 -2271.594400 -8879.882398 7860.094024 -10532.048973
21 22 23 24 25
2058.266456 13282.387520 1656.023871 -12880.335853 752.106490
26 27 28 29 30
8834.946546 1119.818840 6025.020315 1711.242347 -267.153836
31 32 33 34 35
-1072.538372 3971.440740 502.245915 -630.598784 2411.540941
36 37 38 39 40
-2542.537020 1227.347464 -4599.216983 -5569.108487 -7364.993437
41 42 43 44 45
-10630.587780 9208.885620 -4894.132081 -1065.661377 1186.024864
46 47 48 49 50
2236.996572 2357.147916 -7117.296068 -607.515537 409.976400
51 52 53 54 55
-6466.885385 -9897.539900 -1698.858363 -5120.979530 -4274.606629
56 57 58 59 60
-5388.027116 1164.608914 752.664720 -2681.840140 1190.958877
61 62 63 64 65
1027.238158 776.879414 1183.942165 398.693482 -2277.969651
66 67 68 69 70
-12172.444652 1160.090672 -1688.559331 -4798.869066 -4301.318376
71 72 73 74 75
-458.143035 1368.450188 -530.194060 852.647411 -473.041342
76 77 78 79 80
1565.123474 -733.977048 1599.637498 -3002.039193 892.052746
81 82 83 84 85
-893.680778 -1838.314630 -869.274515 137.057585 -19.699991
86 87 88 89 90
-775.659513 -6773.945940 428.217089 -967.413640 51.028841
91 92 93 94 95
-1011.778298 -523.503202 -5090.946949 4906.913402 -2323.635098
96 97 98 99 100
349.481789 3.965096 83.506320 1676.896306 1572.201226
101 102 103 104 105
-753.885336 1052.843331 -706.060809 1413.306144 -658.620191
106 107 108 109 110
-800.660274 753.225915 75.980592 2163.123179 1515.395401
111 112 113 114 115
-1165.961738 2121.866454 1958.485205 2369.754650 1103.493217
116 117 118 119 120
-1566.136330 652.282041 605.790104 3105.698267 -456.590262
121 122 123 124 125
1872.593323 436.399257 -9298.175625 -930.367574 1448.778830
126 127 128 129 130
2824.195183 898.420356 -61.740547 -1182.017152 886.834418
131 132 133 134 135
183.775794 -753.199877 875.248480 -1178.516846 1118.239960
136 137 138 139 140
756.671529 -112.151741 759.422504 1329.712566 848.214625
141 142 143 144 145
274.528869 840.490666 771.289147 443.512513 -392.536825
146 147 148 149 150
149.212224 821.180769 90.168854 785.498746 823.028279
151 152 153 154 155
767.574032 -2222.665716 268.645232 639.626590 602.549860
156 157 158 159 160
786.422955 782.560976 676.177402 466.113640 770.975038
161 162 163 164 165
767.113058 763.251079 759.389100 579.362587 -1185.854145
166 167 168 169 170
693.331268 743.941182 -1197.902430 -663.224739 732.355244
171 172 173 174 175
728.493265 480.149080 720.769306 33.907297 713.045348
176 177 178 179 180
-10.288872 -10839.624726 765.748233 697.597430 693.735451
181 182 183 184 185
-596.985960 63.754208 442.283834 -334.762033 196.849548
186 187 188 189 190
-18.123973 260.735104 211.773811 629.048084 655.115657
191 192 193 194 195
139.640907 647.391699 55.304984 639.667740 345.662809
196 197 198 199 200
631.943781 -3965.196616 624.219823 620.357843 -268.659414
201 202 203 204 205
612.633885 608.771905 604.909926 -467.165174 78.146239
206 207 208 209 210
79.208324 476.696604 585.600029 -1442.093427 599.771669
211 212 213 214 215
574.014091 570.152112 -500.833967 227.199344 684.037151
216 217 218 219 220
342.958730 -70.681439 181.094716 331.372792 -678.698942
221 222 223 224 225
668.526557 531.532318 424.743046 -1097.199902 433.440116
226 227 228 229 230
-1522.228462 -372.306275 508.360442 570.507298 500.636484
231 232 233 234 235
261.358549 225.071132 489.050546 -403.302692 481.326587
236 237 238 239 240
478.816294 473.602628 -120.940172 465.878670 462.016690
241 242 243 244 245
458.154711 260.976481 450.430752 -284.468751 2637.952759
246 247 248 249 250
438.844814 434.982835 -955.431590 427.258876 423.396897
251 252 253 254 255
-260.611959 415.672938 589.703426 -2737.177295 -2472.122247
256 257 258 259 260
400.225021 -4.565114 -2014.924002 388.639083 384.777103
261 262 263 264 265
380.915124 377.053145 -872.870551 -552.961821 365.467207
266 267 268 269 270
361.605227 357.743248 -773.584747 350.019289 -2292.711657
271 272 273 274 275
342.295331 338.433351 334.571372 -1509.097510 326.847413
276 277 278 279 280
-684.071752 -1780.641795 315.261475 -1140.046963 -225.978934
281 282 283 284 285
303.675537 299.813558 -2242.664042 292.089599 -342.101433
286 287 288 289 290
284.365640 318.294194 -17.871493 272.779702 -741.154919
291 292 293 294 295
265.055744 -903.435107 -1484.507552 253.469806 -657.183957
296 297 298 299 300
-568.742312 241.883868 238.021888 234.159909 -722.556661
301 302 303 304 305
-146.464074 507.385422 112.839259 214.850012 -1147.840501
306 307 308 309 310
-1061.276595 203.264074 -1571.779432 -1249.324944 191.678136
311 312 313 314 315
187.816157 -2772.843126 180.092198 -1271.990122 -673.522761
316 317 318 319 320
-266.539167 190.860250 -1421.716898 -371.385567 153.058343
321 322 323 324 325
149.196363 145.334384 141.472405 137.610425 133.748446
326 327 328 329 330
71.472540 126.024487 122.162508 118.300529 114.438549
331 332 333 334 335
110.576570 106.714591 102.852611 -218.317390 -9.724309
336 337 338 339 340
-18.144723 -1013.367370 83.542715 15.338132 -891.135674
341 342 343 344 345
-690.477252 68.094797 64.232818 60.370838 -635.881694
346 347 348 349 350
52.646880 -67.117545 44.922921 41.060942 538.220747
351 352 353 354 355
-28.956105 -634.650820 -551.842853 78.600165 17.889066
356 357 358 359 360
14.027086 -892.363260 6.303128 2.441148 -1.420831
361 362 363 364 365
-427.906958 -9.144790 -13.006769 -1144.684290 -816.497706
366 367 368 369 370
-24.592707 -28.454686 -32.316666 -36.178645 -564.678556
371 372 373 374 375
-43.902604 -423.537470 -547.803157 -75.811721 -59.350521
376 377 378 379 380
-63.212501 -992.545636 -270.221644 1038.849018 -78.660418
381 382 383 384 385
-82.522397 -1314.680639 -941.020790 -94.108335 -1282.086145
386 387 388 389 390
91.115155 -409.312452 207.111753 -352.858578 -354.083199
391 392 393 394 395
-818.541933 -122.861476 109.568825 -441.059198 2909.611869
396 397 398 399 400
-173.737927 2551.660241 -24316.428277 -154.482048 -2134.985671
401 402 403 404 405
807.254631 -349.414307 -48.355455 -435.550501 61.626919
406 407 408 409 410
-135.517885 -647.106399 -395.741461 -583.790427 1071.596039
411 412 413 414 415
-4759.314072 5277.312480 -270.930752 -1274.484240 2108.587284
416 417 418 419 420
1681.995813 -540.849417 -967.567234 -1036.460926 997.523489
421 422 423 424 425
1334.204129 -4418.629989 6634.053519 7256.200101 199.469731
426 427 428 429 430
-2863.361271 3802.339283 10855.376843 6385.521807 5569.343320
431
12910.421832
> postscript(file="/var/www/html/rcomp/tmp/6hyqr1291213708.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 27742.767049 NA
1 -50962.285449 27742.767049
2 12952.860263 -50962.285449
3 40491.729428 12952.860263
4 13143.453408 40491.729428
5 15846.235052 13143.453408
6 5778.285078 15846.235052
7 6289.315915 5778.285078
8 9010.028560 6289.315915
9 -6169.038316 9010.028560
10 11858.409601 -6169.038316
11 -4975.490298 11858.409601
12 -32693.119223 -4975.490298
13 25192.057311 -32693.119223
14 -8478.469344 25192.057311
15 -4276.584071 -8478.469344
16 -2271.594400 -4276.584071
17 -8879.882398 -2271.594400
18 7860.094024 -8879.882398
19 -10532.048973 7860.094024
20 2058.266456 -10532.048973
21 13282.387520 2058.266456
22 1656.023871 13282.387520
23 -12880.335853 1656.023871
24 752.106490 -12880.335853
25 8834.946546 752.106490
26 1119.818840 8834.946546
27 6025.020315 1119.818840
28 1711.242347 6025.020315
29 -267.153836 1711.242347
30 -1072.538372 -267.153836
31 3971.440740 -1072.538372
32 502.245915 3971.440740
33 -630.598784 502.245915
34 2411.540941 -630.598784
35 -2542.537020 2411.540941
36 1227.347464 -2542.537020
37 -4599.216983 1227.347464
38 -5569.108487 -4599.216983
39 -7364.993437 -5569.108487
40 -10630.587780 -7364.993437
41 9208.885620 -10630.587780
42 -4894.132081 9208.885620
43 -1065.661377 -4894.132081
44 1186.024864 -1065.661377
45 2236.996572 1186.024864
46 2357.147916 2236.996572
47 -7117.296068 2357.147916
48 -607.515537 -7117.296068
49 409.976400 -607.515537
50 -6466.885385 409.976400
51 -9897.539900 -6466.885385
52 -1698.858363 -9897.539900
53 -5120.979530 -1698.858363
54 -4274.606629 -5120.979530
55 -5388.027116 -4274.606629
56 1164.608914 -5388.027116
57 752.664720 1164.608914
58 -2681.840140 752.664720
59 1190.958877 -2681.840140
60 1027.238158 1190.958877
61 776.879414 1027.238158
62 1183.942165 776.879414
63 398.693482 1183.942165
64 -2277.969651 398.693482
65 -12172.444652 -2277.969651
66 1160.090672 -12172.444652
67 -1688.559331 1160.090672
68 -4798.869066 -1688.559331
69 -4301.318376 -4798.869066
70 -458.143035 -4301.318376
71 1368.450188 -458.143035
72 -530.194060 1368.450188
73 852.647411 -530.194060
74 -473.041342 852.647411
75 1565.123474 -473.041342
76 -733.977048 1565.123474
77 1599.637498 -733.977048
78 -3002.039193 1599.637498
79 892.052746 -3002.039193
80 -893.680778 892.052746
81 -1838.314630 -893.680778
82 -869.274515 -1838.314630
83 137.057585 -869.274515
84 -19.699991 137.057585
85 -775.659513 -19.699991
86 -6773.945940 -775.659513
87 428.217089 -6773.945940
88 -967.413640 428.217089
89 51.028841 -967.413640
90 -1011.778298 51.028841
91 -523.503202 -1011.778298
92 -5090.946949 -523.503202
93 4906.913402 -5090.946949
94 -2323.635098 4906.913402
95 349.481789 -2323.635098
96 3.965096 349.481789
97 83.506320 3.965096
98 1676.896306 83.506320
99 1572.201226 1676.896306
100 -753.885336 1572.201226
101 1052.843331 -753.885336
102 -706.060809 1052.843331
103 1413.306144 -706.060809
104 -658.620191 1413.306144
105 -800.660274 -658.620191
106 753.225915 -800.660274
107 75.980592 753.225915
108 2163.123179 75.980592
109 1515.395401 2163.123179
110 -1165.961738 1515.395401
111 2121.866454 -1165.961738
112 1958.485205 2121.866454
113 2369.754650 1958.485205
114 1103.493217 2369.754650
115 -1566.136330 1103.493217
116 652.282041 -1566.136330
117 605.790104 652.282041
118 3105.698267 605.790104
119 -456.590262 3105.698267
120 1872.593323 -456.590262
121 436.399257 1872.593323
122 -9298.175625 436.399257
123 -930.367574 -9298.175625
124 1448.778830 -930.367574
125 2824.195183 1448.778830
126 898.420356 2824.195183
127 -61.740547 898.420356
128 -1182.017152 -61.740547
129 886.834418 -1182.017152
130 183.775794 886.834418
131 -753.199877 183.775794
132 875.248480 -753.199877
133 -1178.516846 875.248480
134 1118.239960 -1178.516846
135 756.671529 1118.239960
136 -112.151741 756.671529
137 759.422504 -112.151741
138 1329.712566 759.422504
139 848.214625 1329.712566
140 274.528869 848.214625
141 840.490666 274.528869
142 771.289147 840.490666
143 443.512513 771.289147
144 -392.536825 443.512513
145 149.212224 -392.536825
146 821.180769 149.212224
147 90.168854 821.180769
148 785.498746 90.168854
149 823.028279 785.498746
150 767.574032 823.028279
151 -2222.665716 767.574032
152 268.645232 -2222.665716
153 639.626590 268.645232
154 602.549860 639.626590
155 786.422955 602.549860
156 782.560976 786.422955
157 676.177402 782.560976
158 466.113640 676.177402
159 770.975038 466.113640
160 767.113058 770.975038
161 763.251079 767.113058
162 759.389100 763.251079
163 579.362587 759.389100
164 -1185.854145 579.362587
165 693.331268 -1185.854145
166 743.941182 693.331268
167 -1197.902430 743.941182
168 -663.224739 -1197.902430
169 732.355244 -663.224739
170 728.493265 732.355244
171 480.149080 728.493265
172 720.769306 480.149080
173 33.907297 720.769306
174 713.045348 33.907297
175 -10.288872 713.045348
176 -10839.624726 -10.288872
177 765.748233 -10839.624726
178 697.597430 765.748233
179 693.735451 697.597430
180 -596.985960 693.735451
181 63.754208 -596.985960
182 442.283834 63.754208
183 -334.762033 442.283834
184 196.849548 -334.762033
185 -18.123973 196.849548
186 260.735104 -18.123973
187 211.773811 260.735104
188 629.048084 211.773811
189 655.115657 629.048084
190 139.640907 655.115657
191 647.391699 139.640907
192 55.304984 647.391699
193 639.667740 55.304984
194 345.662809 639.667740
195 631.943781 345.662809
196 -3965.196616 631.943781
197 624.219823 -3965.196616
198 620.357843 624.219823
199 -268.659414 620.357843
200 612.633885 -268.659414
201 608.771905 612.633885
202 604.909926 608.771905
203 -467.165174 604.909926
204 78.146239 -467.165174
205 79.208324 78.146239
206 476.696604 79.208324
207 585.600029 476.696604
208 -1442.093427 585.600029
209 599.771669 -1442.093427
210 574.014091 599.771669
211 570.152112 574.014091
212 -500.833967 570.152112
213 227.199344 -500.833967
214 684.037151 227.199344
215 342.958730 684.037151
216 -70.681439 342.958730
217 181.094716 -70.681439
218 331.372792 181.094716
219 -678.698942 331.372792
220 668.526557 -678.698942
221 531.532318 668.526557
222 424.743046 531.532318
223 -1097.199902 424.743046
224 433.440116 -1097.199902
225 -1522.228462 433.440116
226 -372.306275 -1522.228462
227 508.360442 -372.306275
228 570.507298 508.360442
229 500.636484 570.507298
230 261.358549 500.636484
231 225.071132 261.358549
232 489.050546 225.071132
233 -403.302692 489.050546
234 481.326587 -403.302692
235 478.816294 481.326587
236 473.602628 478.816294
237 -120.940172 473.602628
238 465.878670 -120.940172
239 462.016690 465.878670
240 458.154711 462.016690
241 260.976481 458.154711
242 450.430752 260.976481
243 -284.468751 450.430752
244 2637.952759 -284.468751
245 438.844814 2637.952759
246 434.982835 438.844814
247 -955.431590 434.982835
248 427.258876 -955.431590
249 423.396897 427.258876
250 -260.611959 423.396897
251 415.672938 -260.611959
252 589.703426 415.672938
253 -2737.177295 589.703426
254 -2472.122247 -2737.177295
255 400.225021 -2472.122247
256 -4.565114 400.225021
257 -2014.924002 -4.565114
258 388.639083 -2014.924002
259 384.777103 388.639083
260 380.915124 384.777103
261 377.053145 380.915124
262 -872.870551 377.053145
263 -552.961821 -872.870551
264 365.467207 -552.961821
265 361.605227 365.467207
266 357.743248 361.605227
267 -773.584747 357.743248
268 350.019289 -773.584747
269 -2292.711657 350.019289
270 342.295331 -2292.711657
271 338.433351 342.295331
272 334.571372 338.433351
273 -1509.097510 334.571372
274 326.847413 -1509.097510
275 -684.071752 326.847413
276 -1780.641795 -684.071752
277 315.261475 -1780.641795
278 -1140.046963 315.261475
279 -225.978934 -1140.046963
280 303.675537 -225.978934
281 299.813558 303.675537
282 -2242.664042 299.813558
283 292.089599 -2242.664042
284 -342.101433 292.089599
285 284.365640 -342.101433
286 318.294194 284.365640
287 -17.871493 318.294194
288 272.779702 -17.871493
289 -741.154919 272.779702
290 265.055744 -741.154919
291 -903.435107 265.055744
292 -1484.507552 -903.435107
293 253.469806 -1484.507552
294 -657.183957 253.469806
295 -568.742312 -657.183957
296 241.883868 -568.742312
297 238.021888 241.883868
298 234.159909 238.021888
299 -722.556661 234.159909
300 -146.464074 -722.556661
301 507.385422 -146.464074
302 112.839259 507.385422
303 214.850012 112.839259
304 -1147.840501 214.850012
305 -1061.276595 -1147.840501
306 203.264074 -1061.276595
307 -1571.779432 203.264074
308 -1249.324944 -1571.779432
309 191.678136 -1249.324944
310 187.816157 191.678136
311 -2772.843126 187.816157
312 180.092198 -2772.843126
313 -1271.990122 180.092198
314 -673.522761 -1271.990122
315 -266.539167 -673.522761
316 190.860250 -266.539167
317 -1421.716898 190.860250
318 -371.385567 -1421.716898
319 153.058343 -371.385567
320 149.196363 153.058343
321 145.334384 149.196363
322 141.472405 145.334384
323 137.610425 141.472405
324 133.748446 137.610425
325 71.472540 133.748446
326 126.024487 71.472540
327 122.162508 126.024487
328 118.300529 122.162508
329 114.438549 118.300529
330 110.576570 114.438549
331 106.714591 110.576570
332 102.852611 106.714591
333 -218.317390 102.852611
334 -9.724309 -218.317390
335 -18.144723 -9.724309
336 -1013.367370 -18.144723
337 83.542715 -1013.367370
338 15.338132 83.542715
339 -891.135674 15.338132
340 -690.477252 -891.135674
341 68.094797 -690.477252
342 64.232818 68.094797
343 60.370838 64.232818
344 -635.881694 60.370838
345 52.646880 -635.881694
346 -67.117545 52.646880
347 44.922921 -67.117545
348 41.060942 44.922921
349 538.220747 41.060942
350 -28.956105 538.220747
351 -634.650820 -28.956105
352 -551.842853 -634.650820
353 78.600165 -551.842853
354 17.889066 78.600165
355 14.027086 17.889066
356 -892.363260 14.027086
357 6.303128 -892.363260
358 2.441148 6.303128
359 -1.420831 2.441148
360 -427.906958 -1.420831
361 -9.144790 -427.906958
362 -13.006769 -9.144790
363 -1144.684290 -13.006769
364 -816.497706 -1144.684290
365 -24.592707 -816.497706
366 -28.454686 -24.592707
367 -32.316666 -28.454686
368 -36.178645 -32.316666
369 -564.678556 -36.178645
370 -43.902604 -564.678556
371 -423.537470 -43.902604
372 -547.803157 -423.537470
373 -75.811721 -547.803157
374 -59.350521 -75.811721
375 -63.212501 -59.350521
376 -992.545636 -63.212501
377 -270.221644 -992.545636
378 1038.849018 -270.221644
379 -78.660418 1038.849018
380 -82.522397 -78.660418
381 -1314.680639 -82.522397
382 -941.020790 -1314.680639
383 -94.108335 -941.020790
384 -1282.086145 -94.108335
385 91.115155 -1282.086145
386 -409.312452 91.115155
387 207.111753 -409.312452
388 -352.858578 207.111753
389 -354.083199 -352.858578
390 -818.541933 -354.083199
391 -122.861476 -818.541933
392 109.568825 -122.861476
393 -441.059198 109.568825
394 2909.611869 -441.059198
395 -173.737927 2909.611869
396 2551.660241 -173.737927
397 -24316.428277 2551.660241
398 -154.482048 -24316.428277
399 -2134.985671 -154.482048
400 807.254631 -2134.985671
401 -349.414307 807.254631
402 -48.355455 -349.414307
403 -435.550501 -48.355455
404 61.626919 -435.550501
405 -135.517885 61.626919
406 -647.106399 -135.517885
407 -395.741461 -647.106399
408 -583.790427 -395.741461
409 1071.596039 -583.790427
410 -4759.314072 1071.596039
411 5277.312480 -4759.314072
412 -270.930752 5277.312480
413 -1274.484240 -270.930752
414 2108.587284 -1274.484240
415 1681.995813 2108.587284
416 -540.849417 1681.995813
417 -967.567234 -540.849417
418 -1036.460926 -967.567234
419 997.523489 -1036.460926
420 1334.204129 997.523489
421 -4418.629989 1334.204129
422 6634.053519 -4418.629989
423 7256.200101 6634.053519
424 199.469731 7256.200101
425 -2863.361271 199.469731
426 3802.339283 -2863.361271
427 10855.376843 3802.339283
428 6385.521807 10855.376843
429 5569.343320 6385.521807
430 12910.421832 5569.343320
431 NA 12910.421832
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -50962.285449 27742.767049
[2,] 12952.860263 -50962.285449
[3,] 40491.729428 12952.860263
[4,] 13143.453408 40491.729428
[5,] 15846.235052 13143.453408
[6,] 5778.285078 15846.235052
[7,] 6289.315915 5778.285078
[8,] 9010.028560 6289.315915
[9,] -6169.038316 9010.028560
[10,] 11858.409601 -6169.038316
[11,] -4975.490298 11858.409601
[12,] -32693.119223 -4975.490298
[13,] 25192.057311 -32693.119223
[14,] -8478.469344 25192.057311
[15,] -4276.584071 -8478.469344
[16,] -2271.594400 -4276.584071
[17,] -8879.882398 -2271.594400
[18,] 7860.094024 -8879.882398
[19,] -10532.048973 7860.094024
[20,] 2058.266456 -10532.048973
[21,] 13282.387520 2058.266456
[22,] 1656.023871 13282.387520
[23,] -12880.335853 1656.023871
[24,] 752.106490 -12880.335853
[25,] 8834.946546 752.106490
[26,] 1119.818840 8834.946546
[27,] 6025.020315 1119.818840
[28,] 1711.242347 6025.020315
[29,] -267.153836 1711.242347
[30,] -1072.538372 -267.153836
[31,] 3971.440740 -1072.538372
[32,] 502.245915 3971.440740
[33,] -630.598784 502.245915
[34,] 2411.540941 -630.598784
[35,] -2542.537020 2411.540941
[36,] 1227.347464 -2542.537020
[37,] -4599.216983 1227.347464
[38,] -5569.108487 -4599.216983
[39,] -7364.993437 -5569.108487
[40,] -10630.587780 -7364.993437
[41,] 9208.885620 -10630.587780
[42,] -4894.132081 9208.885620
[43,] -1065.661377 -4894.132081
[44,] 1186.024864 -1065.661377
[45,] 2236.996572 1186.024864
[46,] 2357.147916 2236.996572
[47,] -7117.296068 2357.147916
[48,] -607.515537 -7117.296068
[49,] 409.976400 -607.515537
[50,] -6466.885385 409.976400
[51,] -9897.539900 -6466.885385
[52,] -1698.858363 -9897.539900
[53,] -5120.979530 -1698.858363
[54,] -4274.606629 -5120.979530
[55,] -5388.027116 -4274.606629
[56,] 1164.608914 -5388.027116
[57,] 752.664720 1164.608914
[58,] -2681.840140 752.664720
[59,] 1190.958877 -2681.840140
[60,] 1027.238158 1190.958877
[61,] 776.879414 1027.238158
[62,] 1183.942165 776.879414
[63,] 398.693482 1183.942165
[64,] -2277.969651 398.693482
[65,] -12172.444652 -2277.969651
[66,] 1160.090672 -12172.444652
[67,] -1688.559331 1160.090672
[68,] -4798.869066 -1688.559331
[69,] -4301.318376 -4798.869066
[70,] -458.143035 -4301.318376
[71,] 1368.450188 -458.143035
[72,] -530.194060 1368.450188
[73,] 852.647411 -530.194060
[74,] -473.041342 852.647411
[75,] 1565.123474 -473.041342
[76,] -733.977048 1565.123474
[77,] 1599.637498 -733.977048
[78,] -3002.039193 1599.637498
[79,] 892.052746 -3002.039193
[80,] -893.680778 892.052746
[81,] -1838.314630 -893.680778
[82,] -869.274515 -1838.314630
[83,] 137.057585 -869.274515
[84,] -19.699991 137.057585
[85,] -775.659513 -19.699991
[86,] -6773.945940 -775.659513
[87,] 428.217089 -6773.945940
[88,] -967.413640 428.217089
[89,] 51.028841 -967.413640
[90,] -1011.778298 51.028841
[91,] -523.503202 -1011.778298
[92,] -5090.946949 -523.503202
[93,] 4906.913402 -5090.946949
[94,] -2323.635098 4906.913402
[95,] 349.481789 -2323.635098
[96,] 3.965096 349.481789
[97,] 83.506320 3.965096
[98,] 1676.896306 83.506320
[99,] 1572.201226 1676.896306
[100,] -753.885336 1572.201226
[101,] 1052.843331 -753.885336
[102,] -706.060809 1052.843331
[103,] 1413.306144 -706.060809
[104,] -658.620191 1413.306144
[105,] -800.660274 -658.620191
[106,] 753.225915 -800.660274
[107,] 75.980592 753.225915
[108,] 2163.123179 75.980592
[109,] 1515.395401 2163.123179
[110,] -1165.961738 1515.395401
[111,] 2121.866454 -1165.961738
[112,] 1958.485205 2121.866454
[113,] 2369.754650 1958.485205
[114,] 1103.493217 2369.754650
[115,] -1566.136330 1103.493217
[116,] 652.282041 -1566.136330
[117,] 605.790104 652.282041
[118,] 3105.698267 605.790104
[119,] -456.590262 3105.698267
[120,] 1872.593323 -456.590262
[121,] 436.399257 1872.593323
[122,] -9298.175625 436.399257
[123,] -930.367574 -9298.175625
[124,] 1448.778830 -930.367574
[125,] 2824.195183 1448.778830
[126,] 898.420356 2824.195183
[127,] -61.740547 898.420356
[128,] -1182.017152 -61.740547
[129,] 886.834418 -1182.017152
[130,] 183.775794 886.834418
[131,] -753.199877 183.775794
[132,] 875.248480 -753.199877
[133,] -1178.516846 875.248480
[134,] 1118.239960 -1178.516846
[135,] 756.671529 1118.239960
[136,] -112.151741 756.671529
[137,] 759.422504 -112.151741
[138,] 1329.712566 759.422504
[139,] 848.214625 1329.712566
[140,] 274.528869 848.214625
[141,] 840.490666 274.528869
[142,] 771.289147 840.490666
[143,] 443.512513 771.289147
[144,] -392.536825 443.512513
[145,] 149.212224 -392.536825
[146,] 821.180769 149.212224
[147,] 90.168854 821.180769
[148,] 785.498746 90.168854
[149,] 823.028279 785.498746
[150,] 767.574032 823.028279
[151,] -2222.665716 767.574032
[152,] 268.645232 -2222.665716
[153,] 639.626590 268.645232
[154,] 602.549860 639.626590
[155,] 786.422955 602.549860
[156,] 782.560976 786.422955
[157,] 676.177402 782.560976
[158,] 466.113640 676.177402
[159,] 770.975038 466.113640
[160,] 767.113058 770.975038
[161,] 763.251079 767.113058
[162,] 759.389100 763.251079
[163,] 579.362587 759.389100
[164,] -1185.854145 579.362587
[165,] 693.331268 -1185.854145
[166,] 743.941182 693.331268
[167,] -1197.902430 743.941182
[168,] -663.224739 -1197.902430
[169,] 732.355244 -663.224739
[170,] 728.493265 732.355244
[171,] 480.149080 728.493265
[172,] 720.769306 480.149080
[173,] 33.907297 720.769306
[174,] 713.045348 33.907297
[175,] -10.288872 713.045348
[176,] -10839.624726 -10.288872
[177,] 765.748233 -10839.624726
[178,] 697.597430 765.748233
[179,] 693.735451 697.597430
[180,] -596.985960 693.735451
[181,] 63.754208 -596.985960
[182,] 442.283834 63.754208
[183,] -334.762033 442.283834
[184,] 196.849548 -334.762033
[185,] -18.123973 196.849548
[186,] 260.735104 -18.123973
[187,] 211.773811 260.735104
[188,] 629.048084 211.773811
[189,] 655.115657 629.048084
[190,] 139.640907 655.115657
[191,] 647.391699 139.640907
[192,] 55.304984 647.391699
[193,] 639.667740 55.304984
[194,] 345.662809 639.667740
[195,] 631.943781 345.662809
[196,] -3965.196616 631.943781
[197,] 624.219823 -3965.196616
[198,] 620.357843 624.219823
[199,] -268.659414 620.357843
[200,] 612.633885 -268.659414
[201,] 608.771905 612.633885
[202,] 604.909926 608.771905
[203,] -467.165174 604.909926
[204,] 78.146239 -467.165174
[205,] 79.208324 78.146239
[206,] 476.696604 79.208324
[207,] 585.600029 476.696604
[208,] -1442.093427 585.600029
[209,] 599.771669 -1442.093427
[210,] 574.014091 599.771669
[211,] 570.152112 574.014091
[212,] -500.833967 570.152112
[213,] 227.199344 -500.833967
[214,] 684.037151 227.199344
[215,] 342.958730 684.037151
[216,] -70.681439 342.958730
[217,] 181.094716 -70.681439
[218,] 331.372792 181.094716
[219,] -678.698942 331.372792
[220,] 668.526557 -678.698942
[221,] 531.532318 668.526557
[222,] 424.743046 531.532318
[223,] -1097.199902 424.743046
[224,] 433.440116 -1097.199902
[225,] -1522.228462 433.440116
[226,] -372.306275 -1522.228462
[227,] 508.360442 -372.306275
[228,] 570.507298 508.360442
[229,] 500.636484 570.507298
[230,] 261.358549 500.636484
[231,] 225.071132 261.358549
[232,] 489.050546 225.071132
[233,] -403.302692 489.050546
[234,] 481.326587 -403.302692
[235,] 478.816294 481.326587
[236,] 473.602628 478.816294
[237,] -120.940172 473.602628
[238,] 465.878670 -120.940172
[239,] 462.016690 465.878670
[240,] 458.154711 462.016690
[241,] 260.976481 458.154711
[242,] 450.430752 260.976481
[243,] -284.468751 450.430752
[244,] 2637.952759 -284.468751
[245,] 438.844814 2637.952759
[246,] 434.982835 438.844814
[247,] -955.431590 434.982835
[248,] 427.258876 -955.431590
[249,] 423.396897 427.258876
[250,] -260.611959 423.396897
[251,] 415.672938 -260.611959
[252,] 589.703426 415.672938
[253,] -2737.177295 589.703426
[254,] -2472.122247 -2737.177295
[255,] 400.225021 -2472.122247
[256,] -4.565114 400.225021
[257,] -2014.924002 -4.565114
[258,] 388.639083 -2014.924002
[259,] 384.777103 388.639083
[260,] 380.915124 384.777103
[261,] 377.053145 380.915124
[262,] -872.870551 377.053145
[263,] -552.961821 -872.870551
[264,] 365.467207 -552.961821
[265,] 361.605227 365.467207
[266,] 357.743248 361.605227
[267,] -773.584747 357.743248
[268,] 350.019289 -773.584747
[269,] -2292.711657 350.019289
[270,] 342.295331 -2292.711657
[271,] 338.433351 342.295331
[272,] 334.571372 338.433351
[273,] -1509.097510 334.571372
[274,] 326.847413 -1509.097510
[275,] -684.071752 326.847413
[276,] -1780.641795 -684.071752
[277,] 315.261475 -1780.641795
[278,] -1140.046963 315.261475
[279,] -225.978934 -1140.046963
[280,] 303.675537 -225.978934
[281,] 299.813558 303.675537
[282,] -2242.664042 299.813558
[283,] 292.089599 -2242.664042
[284,] -342.101433 292.089599
[285,] 284.365640 -342.101433
[286,] 318.294194 284.365640
[287,] -17.871493 318.294194
[288,] 272.779702 -17.871493
[289,] -741.154919 272.779702
[290,] 265.055744 -741.154919
[291,] -903.435107 265.055744
[292,] -1484.507552 -903.435107
[293,] 253.469806 -1484.507552
[294,] -657.183957 253.469806
[295,] -568.742312 -657.183957
[296,] 241.883868 -568.742312
[297,] 238.021888 241.883868
[298,] 234.159909 238.021888
[299,] -722.556661 234.159909
[300,] -146.464074 -722.556661
[301,] 507.385422 -146.464074
[302,] 112.839259 507.385422
[303,] 214.850012 112.839259
[304,] -1147.840501 214.850012
[305,] -1061.276595 -1147.840501
[306,] 203.264074 -1061.276595
[307,] -1571.779432 203.264074
[308,] -1249.324944 -1571.779432
[309,] 191.678136 -1249.324944
[310,] 187.816157 191.678136
[311,] -2772.843126 187.816157
[312,] 180.092198 -2772.843126
[313,] -1271.990122 180.092198
[314,] -673.522761 -1271.990122
[315,] -266.539167 -673.522761
[316,] 190.860250 -266.539167
[317,] -1421.716898 190.860250
[318,] -371.385567 -1421.716898
[319,] 153.058343 -371.385567
[320,] 149.196363 153.058343
[321,] 145.334384 149.196363
[322,] 141.472405 145.334384
[323,] 137.610425 141.472405
[324,] 133.748446 137.610425
[325,] 71.472540 133.748446
[326,] 126.024487 71.472540
[327,] 122.162508 126.024487
[328,] 118.300529 122.162508
[329,] 114.438549 118.300529
[330,] 110.576570 114.438549
[331,] 106.714591 110.576570
[332,] 102.852611 106.714591
[333,] -218.317390 102.852611
[334,] -9.724309 -218.317390
[335,] -18.144723 -9.724309
[336,] -1013.367370 -18.144723
[337,] 83.542715 -1013.367370
[338,] 15.338132 83.542715
[339,] -891.135674 15.338132
[340,] -690.477252 -891.135674
[341,] 68.094797 -690.477252
[342,] 64.232818 68.094797
[343,] 60.370838 64.232818
[344,] -635.881694 60.370838
[345,] 52.646880 -635.881694
[346,] -67.117545 52.646880
[347,] 44.922921 -67.117545
[348,] 41.060942 44.922921
[349,] 538.220747 41.060942
[350,] -28.956105 538.220747
[351,] -634.650820 -28.956105
[352,] -551.842853 -634.650820
[353,] 78.600165 -551.842853
[354,] 17.889066 78.600165
[355,] 14.027086 17.889066
[356,] -892.363260 14.027086
[357,] 6.303128 -892.363260
[358,] 2.441148 6.303128
[359,] -1.420831 2.441148
[360,] -427.906958 -1.420831
[361,] -9.144790 -427.906958
[362,] -13.006769 -9.144790
[363,] -1144.684290 -13.006769
[364,] -816.497706 -1144.684290
[365,] -24.592707 -816.497706
[366,] -28.454686 -24.592707
[367,] -32.316666 -28.454686
[368,] -36.178645 -32.316666
[369,] -564.678556 -36.178645
[370,] -43.902604 -564.678556
[371,] -423.537470 -43.902604
[372,] -547.803157 -423.537470
[373,] -75.811721 -547.803157
[374,] -59.350521 -75.811721
[375,] -63.212501 -59.350521
[376,] -992.545636 -63.212501
[377,] -270.221644 -992.545636
[378,] 1038.849018 -270.221644
[379,] -78.660418 1038.849018
[380,] -82.522397 -78.660418
[381,] -1314.680639 -82.522397
[382,] -941.020790 -1314.680639
[383,] -94.108335 -941.020790
[384,] -1282.086145 -94.108335
[385,] 91.115155 -1282.086145
[386,] -409.312452 91.115155
[387,] 207.111753 -409.312452
[388,] -352.858578 207.111753
[389,] -354.083199 -352.858578
[390,] -818.541933 -354.083199
[391,] -122.861476 -818.541933
[392,] 109.568825 -122.861476
[393,] -441.059198 109.568825
[394,] 2909.611869 -441.059198
[395,] -173.737927 2909.611869
[396,] 2551.660241 -173.737927
[397,] -24316.428277 2551.660241
[398,] -154.482048 -24316.428277
[399,] -2134.985671 -154.482048
[400,] 807.254631 -2134.985671
[401,] -349.414307 807.254631
[402,] -48.355455 -349.414307
[403,] -435.550501 -48.355455
[404,] 61.626919 -435.550501
[405,] -135.517885 61.626919
[406,] -647.106399 -135.517885
[407,] -395.741461 -647.106399
[408,] -583.790427 -395.741461
[409,] 1071.596039 -583.790427
[410,] -4759.314072 1071.596039
[411,] 5277.312480 -4759.314072
[412,] -270.930752 5277.312480
[413,] -1274.484240 -270.930752
[414,] 2108.587284 -1274.484240
[415,] 1681.995813 2108.587284
[416,] -540.849417 1681.995813
[417,] -967.567234 -540.849417
[418,] -1036.460926 -967.567234
[419,] 997.523489 -1036.460926
[420,] 1334.204129 997.523489
[421,] -4418.629989 1334.204129
[422,] 6634.053519 -4418.629989
[423,] 7256.200101 6634.053519
[424,] 199.469731 7256.200101
[425,] -2863.361271 199.469731
[426,] 3802.339283 -2863.361271
[427,] 10855.376843 3802.339283
[428,] 6385.521807 10855.376843
[429,] 5569.343320 6385.521807
[430,] 12910.421832 5569.343320
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -50962.285449 27742.767049
2 12952.860263 -50962.285449
3 40491.729428 12952.860263
4 13143.453408 40491.729428
5 15846.235052 13143.453408
6 5778.285078 15846.235052
7 6289.315915 5778.285078
8 9010.028560 6289.315915
9 -6169.038316 9010.028560
10 11858.409601 -6169.038316
11 -4975.490298 11858.409601
12 -32693.119223 -4975.490298
13 25192.057311 -32693.119223
14 -8478.469344 25192.057311
15 -4276.584071 -8478.469344
16 -2271.594400 -4276.584071
17 -8879.882398 -2271.594400
18 7860.094024 -8879.882398
19 -10532.048973 7860.094024
20 2058.266456 -10532.048973
21 13282.387520 2058.266456
22 1656.023871 13282.387520
23 -12880.335853 1656.023871
24 752.106490 -12880.335853
25 8834.946546 752.106490
26 1119.818840 8834.946546
27 6025.020315 1119.818840
28 1711.242347 6025.020315
29 -267.153836 1711.242347
30 -1072.538372 -267.153836
31 3971.440740 -1072.538372
32 502.245915 3971.440740
33 -630.598784 502.245915
34 2411.540941 -630.598784
35 -2542.537020 2411.540941
36 1227.347464 -2542.537020
37 -4599.216983 1227.347464
38 -5569.108487 -4599.216983
39 -7364.993437 -5569.108487
40 -10630.587780 -7364.993437
41 9208.885620 -10630.587780
42 -4894.132081 9208.885620
43 -1065.661377 -4894.132081
44 1186.024864 -1065.661377
45 2236.996572 1186.024864
46 2357.147916 2236.996572
47 -7117.296068 2357.147916
48 -607.515537 -7117.296068
49 409.976400 -607.515537
50 -6466.885385 409.976400
51 -9897.539900 -6466.885385
52 -1698.858363 -9897.539900
53 -5120.979530 -1698.858363
54 -4274.606629 -5120.979530
55 -5388.027116 -4274.606629
56 1164.608914 -5388.027116
57 752.664720 1164.608914
58 -2681.840140 752.664720
59 1190.958877 -2681.840140
60 1027.238158 1190.958877
61 776.879414 1027.238158
62 1183.942165 776.879414
63 398.693482 1183.942165
64 -2277.969651 398.693482
65 -12172.444652 -2277.969651
66 1160.090672 -12172.444652
67 -1688.559331 1160.090672
68 -4798.869066 -1688.559331
69 -4301.318376 -4798.869066
70 -458.143035 -4301.318376
71 1368.450188 -458.143035
72 -530.194060 1368.450188
73 852.647411 -530.194060
74 -473.041342 852.647411
75 1565.123474 -473.041342
76 -733.977048 1565.123474
77 1599.637498 -733.977048
78 -3002.039193 1599.637498
79 892.052746 -3002.039193
80 -893.680778 892.052746
81 -1838.314630 -893.680778
82 -869.274515 -1838.314630
83 137.057585 -869.274515
84 -19.699991 137.057585
85 -775.659513 -19.699991
86 -6773.945940 -775.659513
87 428.217089 -6773.945940
88 -967.413640 428.217089
89 51.028841 -967.413640
90 -1011.778298 51.028841
91 -523.503202 -1011.778298
92 -5090.946949 -523.503202
93 4906.913402 -5090.946949
94 -2323.635098 4906.913402
95 349.481789 -2323.635098
96 3.965096 349.481789
97 83.506320 3.965096
98 1676.896306 83.506320
99 1572.201226 1676.896306
100 -753.885336 1572.201226
101 1052.843331 -753.885336
102 -706.060809 1052.843331
103 1413.306144 -706.060809
104 -658.620191 1413.306144
105 -800.660274 -658.620191
106 753.225915 -800.660274
107 75.980592 753.225915
108 2163.123179 75.980592
109 1515.395401 2163.123179
110 -1165.961738 1515.395401
111 2121.866454 -1165.961738
112 1958.485205 2121.866454
113 2369.754650 1958.485205
114 1103.493217 2369.754650
115 -1566.136330 1103.493217
116 652.282041 -1566.136330
117 605.790104 652.282041
118 3105.698267 605.790104
119 -456.590262 3105.698267
120 1872.593323 -456.590262
121 436.399257 1872.593323
122 -9298.175625 436.399257
123 -930.367574 -9298.175625
124 1448.778830 -930.367574
125 2824.195183 1448.778830
126 898.420356 2824.195183
127 -61.740547 898.420356
128 -1182.017152 -61.740547
129 886.834418 -1182.017152
130 183.775794 886.834418
131 -753.199877 183.775794
132 875.248480 -753.199877
133 -1178.516846 875.248480
134 1118.239960 -1178.516846
135 756.671529 1118.239960
136 -112.151741 756.671529
137 759.422504 -112.151741
138 1329.712566 759.422504
139 848.214625 1329.712566
140 274.528869 848.214625
141 840.490666 274.528869
142 771.289147 840.490666
143 443.512513 771.289147
144 -392.536825 443.512513
145 149.212224 -392.536825
146 821.180769 149.212224
147 90.168854 821.180769
148 785.498746 90.168854
149 823.028279 785.498746
150 767.574032 823.028279
151 -2222.665716 767.574032
152 268.645232 -2222.665716
153 639.626590 268.645232
154 602.549860 639.626590
155 786.422955 602.549860
156 782.560976 786.422955
157 676.177402 782.560976
158 466.113640 676.177402
159 770.975038 466.113640
160 767.113058 770.975038
161 763.251079 767.113058
162 759.389100 763.251079
163 579.362587 759.389100
164 -1185.854145 579.362587
165 693.331268 -1185.854145
166 743.941182 693.331268
167 -1197.902430 743.941182
168 -663.224739 -1197.902430
169 732.355244 -663.224739
170 728.493265 732.355244
171 480.149080 728.493265
172 720.769306 480.149080
173 33.907297 720.769306
174 713.045348 33.907297
175 -10.288872 713.045348
176 -10839.624726 -10.288872
177 765.748233 -10839.624726
178 697.597430 765.748233
179 693.735451 697.597430
180 -596.985960 693.735451
181 63.754208 -596.985960
182 442.283834 63.754208
183 -334.762033 442.283834
184 196.849548 -334.762033
185 -18.123973 196.849548
186 260.735104 -18.123973
187 211.773811 260.735104
188 629.048084 211.773811
189 655.115657 629.048084
190 139.640907 655.115657
191 647.391699 139.640907
192 55.304984 647.391699
193 639.667740 55.304984
194 345.662809 639.667740
195 631.943781 345.662809
196 -3965.196616 631.943781
197 624.219823 -3965.196616
198 620.357843 624.219823
199 -268.659414 620.357843
200 612.633885 -268.659414
201 608.771905 612.633885
202 604.909926 608.771905
203 -467.165174 604.909926
204 78.146239 -467.165174
205 79.208324 78.146239
206 476.696604 79.208324
207 585.600029 476.696604
208 -1442.093427 585.600029
209 599.771669 -1442.093427
210 574.014091 599.771669
211 570.152112 574.014091
212 -500.833967 570.152112
213 227.199344 -500.833967
214 684.037151 227.199344
215 342.958730 684.037151
216 -70.681439 342.958730
217 181.094716 -70.681439
218 331.372792 181.094716
219 -678.698942 331.372792
220 668.526557 -678.698942
221 531.532318 668.526557
222 424.743046 531.532318
223 -1097.199902 424.743046
224 433.440116 -1097.199902
225 -1522.228462 433.440116
226 -372.306275 -1522.228462
227 508.360442 -372.306275
228 570.507298 508.360442
229 500.636484 570.507298
230 261.358549 500.636484
231 225.071132 261.358549
232 489.050546 225.071132
233 -403.302692 489.050546
234 481.326587 -403.302692
235 478.816294 481.326587
236 473.602628 478.816294
237 -120.940172 473.602628
238 465.878670 -120.940172
239 462.016690 465.878670
240 458.154711 462.016690
241 260.976481 458.154711
242 450.430752 260.976481
243 -284.468751 450.430752
244 2637.952759 -284.468751
245 438.844814 2637.952759
246 434.982835 438.844814
247 -955.431590 434.982835
248 427.258876 -955.431590
249 423.396897 427.258876
250 -260.611959 423.396897
251 415.672938 -260.611959
252 589.703426 415.672938
253 -2737.177295 589.703426
254 -2472.122247 -2737.177295
255 400.225021 -2472.122247
256 -4.565114 400.225021
257 -2014.924002 -4.565114
258 388.639083 -2014.924002
259 384.777103 388.639083
260 380.915124 384.777103
261 377.053145 380.915124
262 -872.870551 377.053145
263 -552.961821 -872.870551
264 365.467207 -552.961821
265 361.605227 365.467207
266 357.743248 361.605227
267 -773.584747 357.743248
268 350.019289 -773.584747
269 -2292.711657 350.019289
270 342.295331 -2292.711657
271 338.433351 342.295331
272 334.571372 338.433351
273 -1509.097510 334.571372
274 326.847413 -1509.097510
275 -684.071752 326.847413
276 -1780.641795 -684.071752
277 315.261475 -1780.641795
278 -1140.046963 315.261475
279 -225.978934 -1140.046963
280 303.675537 -225.978934
281 299.813558 303.675537
282 -2242.664042 299.813558
283 292.089599 -2242.664042
284 -342.101433 292.089599
285 284.365640 -342.101433
286 318.294194 284.365640
287 -17.871493 318.294194
288 272.779702 -17.871493
289 -741.154919 272.779702
290 265.055744 -741.154919
291 -903.435107 265.055744
292 -1484.507552 -903.435107
293 253.469806 -1484.507552
294 -657.183957 253.469806
295 -568.742312 -657.183957
296 241.883868 -568.742312
297 238.021888 241.883868
298 234.159909 238.021888
299 -722.556661 234.159909
300 -146.464074 -722.556661
301 507.385422 -146.464074
302 112.839259 507.385422
303 214.850012 112.839259
304 -1147.840501 214.850012
305 -1061.276595 -1147.840501
306 203.264074 -1061.276595
307 -1571.779432 203.264074
308 -1249.324944 -1571.779432
309 191.678136 -1249.324944
310 187.816157 191.678136
311 -2772.843126 187.816157
312 180.092198 -2772.843126
313 -1271.990122 180.092198
314 -673.522761 -1271.990122
315 -266.539167 -673.522761
316 190.860250 -266.539167
317 -1421.716898 190.860250
318 -371.385567 -1421.716898
319 153.058343 -371.385567
320 149.196363 153.058343
321 145.334384 149.196363
322 141.472405 145.334384
323 137.610425 141.472405
324 133.748446 137.610425
325 71.472540 133.748446
326 126.024487 71.472540
327 122.162508 126.024487
328 118.300529 122.162508
329 114.438549 118.300529
330 110.576570 114.438549
331 106.714591 110.576570
332 102.852611 106.714591
333 -218.317390 102.852611
334 -9.724309 -218.317390
335 -18.144723 -9.724309
336 -1013.367370 -18.144723
337 83.542715 -1013.367370
338 15.338132 83.542715
339 -891.135674 15.338132
340 -690.477252 -891.135674
341 68.094797 -690.477252
342 64.232818 68.094797
343 60.370838 64.232818
344 -635.881694 60.370838
345 52.646880 -635.881694
346 -67.117545 52.646880
347 44.922921 -67.117545
348 41.060942 44.922921
349 538.220747 41.060942
350 -28.956105 538.220747
351 -634.650820 -28.956105
352 -551.842853 -634.650820
353 78.600165 -551.842853
354 17.889066 78.600165
355 14.027086 17.889066
356 -892.363260 14.027086
357 6.303128 -892.363260
358 2.441148 6.303128
359 -1.420831 2.441148
360 -427.906958 -1.420831
361 -9.144790 -427.906958
362 -13.006769 -9.144790
363 -1144.684290 -13.006769
364 -816.497706 -1144.684290
365 -24.592707 -816.497706
366 -28.454686 -24.592707
367 -32.316666 -28.454686
368 -36.178645 -32.316666
369 -564.678556 -36.178645
370 -43.902604 -564.678556
371 -423.537470 -43.902604
372 -547.803157 -423.537470
373 -75.811721 -547.803157
374 -59.350521 -75.811721
375 -63.212501 -59.350521
376 -992.545636 -63.212501
377 -270.221644 -992.545636
378 1038.849018 -270.221644
379 -78.660418 1038.849018
380 -82.522397 -78.660418
381 -1314.680639 -82.522397
382 -941.020790 -1314.680639
383 -94.108335 -941.020790
384 -1282.086145 -94.108335
385 91.115155 -1282.086145
386 -409.312452 91.115155
387 207.111753 -409.312452
388 -352.858578 207.111753
389 -354.083199 -352.858578
390 -818.541933 -354.083199
391 -122.861476 -818.541933
392 109.568825 -122.861476
393 -441.059198 109.568825
394 2909.611869 -441.059198
395 -173.737927 2909.611869
396 2551.660241 -173.737927
397 -24316.428277 2551.660241
398 -154.482048 -24316.428277
399 -2134.985671 -154.482048
400 807.254631 -2134.985671
401 -349.414307 807.254631
402 -48.355455 -349.414307
403 -435.550501 -48.355455
404 61.626919 -435.550501
405 -135.517885 61.626919
406 -647.106399 -135.517885
407 -395.741461 -647.106399
408 -583.790427 -395.741461
409 1071.596039 -583.790427
410 -4759.314072 1071.596039
411 5277.312480 -4759.314072
412 -270.930752 5277.312480
413 -1274.484240 -270.930752
414 2108.587284 -1274.484240
415 1681.995813 2108.587284
416 -540.849417 1681.995813
417 -967.567234 -540.849417
418 -1036.460926 -967.567234
419 997.523489 -1036.460926
420 1334.204129 997.523489
421 -4418.629989 1334.204129
422 6634.053519 -4418.629989
423 7256.200101 6634.053519
424 199.469731 7256.200101
425 -2863.361271 199.469731
426 3802.339283 -2863.361271
427 10855.376843 3802.339283
428 6385.521807 10855.376843
429 5569.343320 6385.521807
430 12910.421832 5569.343320
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7s7pu1291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8s7pu1291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/9s7pu1291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/103h6x1291213708.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11oh5l1291213708.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/12s0l81291213708.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/136a1h1291213708.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/149az51291213708.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/1553jo1291213709.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/168lzc1291213709.tab")
+ }
>
> try(system("convert tmp/1ey931291213708.ps tmp/1ey931291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ey931291213708.ps tmp/2ey931291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/3pp861291213708.ps tmp/3pp861291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/4pp861291213708.ps tmp/4pp861291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/5pp861291213708.ps tmp/5pp861291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/6hyqr1291213708.ps tmp/6hyqr1291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/7s7pu1291213708.ps tmp/7s7pu1291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/8s7pu1291213708.ps tmp/8s7pu1291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/9s7pu1291213708.ps tmp/9s7pu1291213708.png",intern=TRUE))
character(0)
> try(system("convert tmp/103h6x1291213708.ps tmp/103h6x1291213708.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.097 2.254 26.171