R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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.
Natural language support but running in an English locale
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 = '4'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
4 1 2 3 t
1 6282154 162556 807 213118 1
2 4321023 29790 444 81767 2
3 4111912 87550 412 153198 3
4 223193 84738 428 -26007 4
5 1491348 54660 315 126942 5
6 1629616 42634 168 157214 6
7 1398893 40949 263 129352 7
8 1926517 45187 267 234817 8
9 983660 37704 228 60448 9
10 1443586 16275 129 47818 10
11 1073089 25830 104 245546 11
12 984885 12679 122 48020 12
13 1405225 18014 393 -1710 13
14 227132 43556 190 32648 14
15 929118 24811 280 95350 15
16 1071292 6575 63 151352 16
17 638830 7123 102 288170 17
18 856956 21950 265 114337 18
19 992426 37597 234 37884 19
20 444477 17821 277 122844 20
21 857217 12988 73 82340 21
22 711969 22330 67 79801 22
23 702380 13326 103 165548 23
24 358589 16189 290 116384 24
25 297978 7146 83 134028 25
26 585715 15824 56 63838 26
27 657954 27664 236 74996 27
28 209458 11920 73 31080 28
29 786690 8568 34 32168 29
30 439798 14416 139 49857 30
31 688779 3369 26 87161 31
32 574339 11819 70 106113 32
33 741409 6984 40 80570 33
34 597793 4519 42 102129 34
35 644190 2220 12 301670 35
36 377934 18562 211 102313 36
37 640273 10327 74 88577 37
38 697458 5336 80 112477 38
39 550608 2365 83 191778 39
40 207393 4069 131 79804 40
41 301607 8636 203 128294 41
42 345783 13718 56 96448 42
43 501749 4525 89 93811 43
44 379983 6869 88 117520 44
45 387475 4628 39 69159 45
46 377305 3689 25 101792 46
47 370837 4891 49 210568 47
48 430866 7489 149 136996 48
49 469107 4901 58 121920 49
50 194493 2284 41 76403 50
51 530670 3160 90 108094 51
52 518365 4150 136 134759 52
53 491303 7285 97 188873 53
54 527021 1134 63 146216 54
55 233773 4658 114 156608 55
56 405972 2384 77 61348 56
57 652925 3748 6 50350 57
58 446211 5371 47 87720 58
59 341340 1285 51 99489 59
60 387699 9327 85 87419 60
61 493408 5565 43 94355 61
62 146494 1528 32 60326 62
63 414462 3122 25 94670 63
64 364304 7561 77 82425 64
65 355178 2675 54 59017 65
66 357760 13253 251 90829 66
67 261216 880 15 80791 67
68 397144 2053 44 100423 68
69 374943 1424 73 131116 69
70 424898 4036 85 100269 70
71 202055 3045 49 27330 71
72 378525 5119 38 39039 72
73 310768 1431 35 106885 73
74 325738 554 9 79285 74
75 394510 1975 34 118881 75
76 247060 1765 20 77623 76
77 368078 1012 29 114768 77
78 236761 810 11 74015 78
79 312378 1280 52 69465 79
80 339836 666 13 117869 80
81 347385 1380 29 60982 81
82 426280 4677 66 90131 82
83 352850 876 33 138971 83
84 301881 814 15 39625 84
85 377516 514 15 102725 85
86 357312 5692 68 64239 86
87 458343 3642 100 90262 87
88 354228 540 13 103960 88
89 308636 2099 45 106611 89
90 386212 567 14 103345 90
91 393343 2001 36 95551 91
92 378509 2949 40 82903 92
93 452469 2253 68 63593 93
94 364839 6533 29 126910 94
95 358649 1889 43 37527 95
96 376641 3055 30 60247 96
97 429112 272 9 112995 97
98 330546 1414 22 70184 98
99 403560 2564 19 130140 99
100 317892 1383 9 73221 100
101 307528 1261 31 76114 101
102 235133 975 19 90534 102
103 299243 3366 55 108479 103
104 314073 576 8 113761 104
105 368186 1686 28 68696 105
106 269661 746 29 71561 106
107 125390 3192 48 59831 107
108 510834 2045 16 97890 108
109 321896 5702 47 101481 109
110 249898 1932 20 72954 110
111 408881 936 22 67939 111
112 158492 3437 33 48022 112
113 292154 5131 44 86111 113
114 289513 2397 13 74020 114
115 378049 1389 6 57530 115
116 343466 1503 35 56364 116
117 332743 402 8 84990 117
118 442882 2239 17 88590 118
119 214215 2234 11 77200 119
120 315688 837 21 61262 120
121 375195 10579 92 110309 121
122 334280 875 12 67000 122
123 355864 1585 112 93099 123
124 480382 1659 25 107577 124
125 353058 2647 17 62920 125
126 217193 3294 23 75832 126
127 315380 0 0 60720 127
128 314533 94 10 60793 128
129 318056 422 23 57935 129
130 315380 0 0 60720 130
131 314353 34 7 60630 131
132 369448 1558 25 55637 132
133 315380 0 0 60720 133
134 312846 43 20 60887 134
135 312075 645 4 60720 135
136 315009 316 4 60505 136
137 318903 115 10 60945 137
138 314887 5 1 60720 138
139 314913 897 4 60720 139
140 315380 0 0 60720 140
141 325506 389 8 58990 141
142 315380 0 0 60720 142
143 298568 1002 11 56750 143
144 315834 36 4 60894 144
145 329784 460 15 63346 145
146 312878 309 9 56535 146
147 315380 0 0 60720 147
148 314987 9 7 60835 148
149 325249 271 2 60720 149
150 315877 14 0 61016 150
151 291650 520 7 58650 151
152 305959 1766 46 60438 152
153 315380 0 5 60720 153
154 297765 458 7 58625 154
155 315245 20 2 60938 155
156 315380 0 0 60720 156
157 315380 0 0 60720 157
158 315236 98 2 61490 158
159 336425 405 5 60845 159
160 315380 0 0 60720 160
161 315380 0 0 60720 161
162 315380 0 0 60720 162
163 315380 0 0 60720 163
164 306268 483 7 60830 164
165 302187 454 24 63261 165
166 314882 47 1 60720 166
167 315380 0 0 60720 167
168 382712 757 18 45689 168
169 341570 4655 55 60720 169
170 315380 0 0 60720 170
171 315380 0 0 60720 171
172 312412 36 3 61564 172
173 315380 0 0 60720 173
174 309596 203 9 61938 174
175 315380 0 0 60720 175
176 315547 126 8 60951 176
177 313267 400 113 60720 177
178 316176 71 0 60745 178
179 315380 0 0 60720 179
180 315380 0 0 60720 180
181 359335 972 19 71642 181
182 330068 531 11 71641 182
183 314289 2461 25 55792 183
184 297413 378 16 71873 184
185 314806 23 5 62555 185
186 333210 638 11 60370 186
187 352108 2300 23 64873 187
188 313332 149 6 62041 188
189 291787 226 5 65745 189
190 315380 0 0 60720 190
191 318745 275 7 59500 191
192 315380 0 0 60720 192
193 315366 141 7 61630 193
194 315380 0 0 60720 194
195 315688 28 3 60890 195
196 315380 0 0 60720 196
197 409642 4980 89 113521 197
198 315380 0 0 60720 198
199 315380 0 0 60720 199
200 269587 472 19 80045 200
201 315380 0 0 60720 201
202 315380 0 0 60720 202
203 315380 0 0 60720 203
204 300962 203 12 50804 204
205 325479 496 12 87390 205
206 316155 10 5 61656 206
207 318574 63 2 65688 207
208 315380 0 0 60720 208
209 343613 1136 26 48522 209
210 306948 265 3 60720 210
211 315380 0 0 60720 211
212 315380 0 0 60720 212
213 330059 267 11 57640 213
214 288985 474 10 61977 214
215 304485 534 5 62620 215
216 315380 0 2 60720 216
217 315688 15 6 60831 217
218 317736 397 7 60646 218
219 315380 0 2 60720 219
220 322331 1866 28 56225 220
221 296656 288 3 60510 221
222 315380 0 0 60720 222
223 315354 3 1 60698 223
224 312161 468 20 60720 224
225 315576 20 1 60805 225
226 314922 278 22 61404 226
227 314551 61 9 60720 227
228 315380 0 0 60720 228
229 312339 192 2 65276 229
230 315380 0 0 60720 230
231 298700 317 7 63915 231
232 321376 738 9 60720 232
233 315380 0 0 60720 233
234 303230 368 13 61686 234
235 315380 0 0 60720 235
236 315487 2 0 60743 236
237 315380 0 0 60720 237
238 315793 53 6 60349 238
239 315380 0 0 60720 239
240 315380 0 0 60720 240
241 315380 0 0 60720 241
242 312887 94 3 61360 242
243 315380 0 0 60720 243
244 315637 24 7 59818 244
245 324385 2332 2 72680 245
246 315380 0 0 60720 246
247 315380 0 0 60720 247
248 308989 131 15 61808 248
249 315380 0 0 60720 249
250 315380 0 0 60720 250
251 296702 206 9 53110 251
252 315380 0 0 60720 252
253 307322 167 1 64245 253
254 304376 622 38 73007 254
255 253588 2328 57 82732 255
256 315380 0 0 60720 256
257 309560 365 7 54820 257
258 298466 364 26 47705 258
259 315380 0 0 60720 259
260 315380 0 0 60720 260
261 315380 0 0 60720 261
262 315380 0 0 60720 262
263 343929 226 13 72835 263
264 331955 307 10 58856 264
265 315380 0 0 60720 265
266 315380 0 0 60720 266
267 315380 0 0 60720 267
268 381180 188 9 77655 268
269 315380 0 0 60720 269
270 331420 138 26 69817 270
271 315380 0 0 60720 271
272 315380 0 0 60720 272
273 315380 0 0 60720 273
274 310201 125 19 60798 274
275 315380 0 0 60720 275
276 320016 282 12 62452 276
277 320398 335 23 64175 277
278 315380 0 0 60720 278
279 291841 1324 29 67440 279
280 310670 176 8 68136 280
281 315380 0 0 60720 281
282 315380 0 0 60720 282
283 313491 249 26 56726 283
284 315380 0 0 60720 284
285 331323 333 9 70811 285
286 315380 0 0 60720 286
287 319210 601 5 60720 287
288 318098 30 3 62045 288
289 315380 0 0 60720 289
290 292754 249 13 54323 290
291 315380 0 0 60720 291
292 325176 165 12 62841 292
293 365959 453 19 81125 293
294 315380 0 0 60720 294
295 302409 53 10 59506 295
296 340968 382 9 59365 296
297 315380 0 0 60720 297
298 315380 0 0 60720 298
299 315380 0 0 60720 299
300 315380 0 9 60720 300
301 313164 30 4 60798 301
302 301164 290 1 58790 302
303 315380 0 1 60720 303
304 315380 0 0 60720 304
305 344425 366 14 61808 305
306 315394 2 12 60735 306
307 315380 0 0 60720 307
308 316647 209 19 64016 308
309 309836 384 17 54683 309
310 315380 0 0 60720 310
311 315380 0 0 60720 311
312 346611 365 32 87192 312
313 315380 0 0 60720 313
314 322031 49 14 64107 314
315 315656 3 8 60761 315
316 339445 133 4 65990 316
317 314964 32 0 59988 317
318 297141 368 20 61167 318
319 315372 1 5 60719 319
320 315380 0 0 60720 320
321 315380 0 0 60720 321
322 315380 0 0 60720 322
323 315380 0 0 60720 323
324 315380 0 0 60720 324
325 315380 0 0 60720 325
326 312502 22 1 60722 326
327 315380 0 0 60720 327
328 315380 0 0 60720 328
329 315380 0 0 60720 329
330 315380 0 0 60720 330
331 315380 0 0 60720 331
332 315380 0 0 60720 332
333 315380 0 0 60720 333
334 313729 96 4 60379 334
335 315388 1 1 60727 335
336 315371 314 4 60720 336
337 296139 844 20 60925 337
338 315380 0 0 60720 338
339 313880 26 1 60896 339
340 317698 125 10 59734 340
341 295580 304 12 62969 341
342 315380 0 0 60720 342
343 315380 0 0 60720 343
344 315380 0 0 60720 344
345 308256 621 13 60720 345
346 315380 0 0 60720 346
347 303677 119 3 59118 347
348 315380 0 0 60720 348
349 315380 0 0 60720 349
350 319369 1595 10 60720 350
351 318690 312 3 58598 351
352 314049 60 7 61124 352
353 325699 587 10 59595 353
354 314210 135 1 62065 354
355 315380 0 0 60720 355
356 315380 0 0 60720 356
357 322378 514 15 78780 357
358 315380 0 0 60720 358
359 315380 0 0 60720 359
360 315380 0 0 60720 360
361 315398 1 4 60722 361
362 315380 0 0 60720 362
363 315380 0 0 60720 363
364 308336 1763 28 61600 364
365 316386 180 9 59635 365
366 315380 0 0 60720 366
367 315380 0 0 60720 367
368 315380 0 0 60720 368
369 315380 0 0 60720 369
370 315553 218 7 60720 370
371 315380 0 0 60720 371
372 323361 448 7 59781 372
373 336639 227 7 76644 373
374 307424 174 3 64820 374
375 315380 0 0 60720 375
376 315380 0 0 60720 376
377 295370 121 11 56178 377
378 322340 607 7 60436 378
379 319864 2212 10 60720 379
380 315380 0 0 60720 380
381 315380 0 0 60720 381
382 317291 530 18 73433 382
383 280398 571 14 41477 383
384 315380 0 0 60720 384
385 317330 78 12 62700 385
386 238125 2489 29 67804 386
387 327071 131 3 59661 387
388 309038 923 6 58620 388
389 314210 72 3 60398 389
390 307930 572 8 58580 390
391 322327 397 10 62710 391
392 292136 450 6 59325 392
393 263276 622 8 60950 393
394 367655 694 6 68060 394
395 283910 3425 9 83620 395
396 283587 562 8 58456 396
397 243650 4917 26 52811 397
398 438493 1442 239 121173 398
399 296261 529 7 63870 399
400 230621 2126 41 21001 400
401 304252 1061 3 70415 401
402 333505 776 8 64230 402
403 296919 611 6 59190 403
404 278990 1526 21 69351 404
405 276898 592 7 64270 405
406 327007 1182 11 70694 406
407 317046 621 11 68005 407
408 304555 989 12 58930 408
409 298096 438 9 58320 409
410 231861 726 3 69980 410
411 309422 1303 57 69863 411
412 286963 7419 21 63255 412
413 269753 1164 15 57320 413
414 448243 3310 32 75230 414
415 165404 1920 11 79420 415
416 204325 965 2 73490 416
417 407159 3256 23 35250 417
418 290476 1135 20 62285 418
419 275311 1270 24 69206 419
420 246541 661 1 65920 420
421 253468 1013 1 69770 421
422 240897 2844 74 72683 422
423 -83265 11528 68 -14545 423
424 -42143 6526 20 55830 424
425 272713 2264 20 55174 425
426 215362 5109 82 67038 426
427 42754 3999 21 51252 427
428 306275 35624 244 157278 428
429 253537 9252 32 79510 429
430 372631 15236 86 77440 430
431 -7170 18073 69 27284 431
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) `1` `2` `3` t
1.521e+05 2.194e+01 1.042e+03 1.614e+00 4.063e+01
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2192330 -24812 40310 54069 2920559
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.521e+05 4.768e+04 3.190 0.001528 **
`1` 2.194e+01 2.167e+00 10.122 < 2e-16 ***
`2` 1.042e+03 3.623e+02 2.877 0.004216 **
`3` 1.614e+00 4.313e-01 3.741 0.000208 ***
t 4.063e+01 1.151e+02 0.353 0.724132
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 253800 on 426 degrees of freedom
Multiple R-squared: 0.6557, Adjusted R-squared: 0.6525
F-statistic: 202.9 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 2.049545e-12 1.024773e-12
[2,] 1.0000000 3.330072e-43 1.665036e-43
[3,] 1.0000000 1.626109e-77 8.130544e-78
[4,] 1.0000000 1.016329e-79 5.081643e-80
[5,] 1.0000000 2.684890e-88 1.342445e-88
[6,] 1.0000000 4.986150e-103 2.493075e-103
[7,] 1.0000000 3.954754e-107 1.977377e-107
[8,] 1.0000000 4.642465e-110 2.321233e-110
[9,] 1.0000000 7.044137e-122 3.522069e-122
[10,] 1.0000000 2.804502e-126 1.402251e-126
[11,] 1.0000000 2.382863e-128 1.191432e-128
[12,] 1.0000000 1.032210e-139 5.161048e-140
[13,] 1.0000000 3.557656e-141 1.778828e-141
[14,] 1.0000000 1.784468e-157 8.922340e-158
[15,] 1.0000000 7.160686e-165 3.580343e-165
[16,] 1.0000000 7.548287e-166 3.774144e-166
[17,] 1.0000000 3.907770e-167 1.953885e-167
[18,] 1.0000000 2.415530e-167 1.207765e-167
[19,] 1.0000000 1.099518e-172 5.497589e-173
[20,] 1.0000000 8.340245e-175 4.170122e-175
[21,] 1.0000000 5.046338e-176 2.523169e-176
[22,] 1.0000000 1.552772e-191 7.763858e-192
[23,] 1.0000000 1.760173e-191 8.800864e-192
[24,] 1.0000000 3.673260e-198 1.836630e-198
[25,] 1.0000000 1.422199e-199 7.110994e-200
[26,] 1.0000000 3.881845e-209 1.940923e-209
[27,] 1.0000000 1.286468e-211 6.432339e-212
[28,] 1.0000000 1.327234e-210 6.636171e-211
[29,] 1.0000000 2.135419e-210 1.067710e-210
[30,] 1.0000000 5.494673e-216 2.747336e-216
[31,] 1.0000000 3.961285e-222 1.980643e-222
[32,] 1.0000000 3.526740e-221 1.763370e-221
[33,] 1.0000000 8.043914e-222 4.021957e-222
[34,] 1.0000000 2.137789e-222 1.068895e-222
[35,] 1.0000000 1.699589e-221 8.497944e-222
[36,] 1.0000000 4.932411e-222 2.466206e-222
[37,] 1.0000000 6.345028e-221 3.172514e-221
[38,] 1.0000000 8.041461e-221 4.020731e-221
[39,] 1.0000000 3.975388e-220 1.987694e-220
[40,] 1.0000000 2.555460e-220 1.277730e-220
[41,] 1.0000000 2.889349e-219 1.444674e-219
[42,] 1.0000000 7.725807e-219 3.862904e-219
[43,] 1.0000000 9.978300e-220 4.989150e-220
[44,] 1.0000000 2.655385e-220 1.327692e-220
[45,] 1.0000000 8.671152e-220 4.335576e-220
[46,] 1.0000000 8.761674e-219 4.380837e-219
[47,] 1.0000000 1.575460e-218 7.877298e-219
[48,] 1.0000000 2.328749e-221 1.164374e-221
[49,] 1.0000000 3.402870e-221 1.701435e-221
[50,] 1.0000000 2.513238e-234 1.256619e-234
[51,] 1.0000000 1.171262e-234 5.856308e-235
[52,] 1.0000000 7.689298e-234 3.844649e-234
[53,] 1.0000000 3.111290e-233 1.555645e-233
[54,] 1.0000000 8.405666e-235 4.202833e-235
[55,] 1.0000000 6.666983e-237 3.333491e-237
[56,] 1.0000000 1.970578e-236 9.852890e-237
[57,] 1.0000000 1.501413e-235 7.507064e-236
[58,] 1.0000000 9.260166e-235 4.630083e-235
[59,] 1.0000000 1.332603e-234 6.663017e-235
[60,] 1.0000000 1.829533e-234 9.147664e-235
[61,] 1.0000000 1.564282e-233 7.821408e-234
[62,] 1.0000000 8.407455e-233 4.203728e-233
[63,] 1.0000000 5.312013e-232 2.656007e-232
[64,] 1.0000000 1.922456e-231 9.612281e-232
[65,] 1.0000000 3.788180e-232 1.894090e-232
[66,] 1.0000000 9.843170e-232 4.921585e-232
[67,] 1.0000000 6.800964e-231 3.400482e-231
[68,] 1.0000000 7.712337e-230 3.856169e-230
[69,] 1.0000000 1.039340e-229 5.196698e-230
[70,] 1.0000000 8.276236e-229 4.138118e-229
[71,] 1.0000000 5.609776e-229 2.804888e-229
[72,] 1.0000000 4.928620e-228 2.464310e-228
[73,] 1.0000000 1.491341e-227 7.456706e-228
[74,] 1.0000000 1.116697e-226 5.583485e-227
[75,] 1.0000000 2.867995e-226 1.433997e-226
[76,] 1.0000000 2.458527e-226 1.229264e-226
[77,] 1.0000000 1.757355e-225 8.786776e-226
[78,] 1.0000000 1.679467e-224 8.397333e-225
[79,] 1.0000000 7.687452e-224 3.843726e-224
[80,] 1.0000000 9.529757e-224 4.764878e-224
[81,] 1.0000000 6.501092e-223 3.250546e-223
[82,] 1.0000000 1.010119e-222 5.050594e-223
[83,] 1.0000000 1.003948e-221 5.019741e-222
[84,] 1.0000000 1.054173e-220 5.270864e-221
[85,] 1.0000000 8.916359e-220 4.458179e-220
[86,] 1.0000000 7.195645e-221 3.597822e-221
[87,] 1.0000000 6.541842e-220 3.270921e-220
[88,] 1.0000000 5.170201e-220 2.585100e-220
[89,] 1.0000000 8.063827e-220 4.031914e-220
[90,] 1.0000000 8.165304e-219 4.082652e-219
[91,] 1.0000000 9.001268e-218 4.500634e-218
[92,] 1.0000000 7.597372e-217 3.798686e-217
[93,] 1.0000000 7.753580e-216 3.876790e-216
[94,] 1.0000000 6.430493e-215 3.215246e-215
[95,] 1.0000000 2.752225e-216 1.376112e-216
[96,] 1.0000000 1.486490e-216 7.432450e-217
[97,] 1.0000000 1.517389e-217 7.586946e-218
[98,] 1.0000000 9.892778e-217 4.946389e-217
[99,] 1.0000000 2.342522e-216 1.171261e-216
[100,] 1.0000000 2.522600e-219 1.261300e-219
[101,] 1.0000000 2.341966e-220 1.170983e-220
[102,] 1.0000000 1.178566e-219 5.892829e-220
[103,] 1.0000000 1.144424e-219 5.722122e-220
[104,] 1.0000000 1.736820e-219 8.684101e-220
[105,] 1.0000000 2.925303e-220 1.462651e-220
[106,] 1.0000000 1.350341e-219 6.751704e-220
[107,] 1.0000000 7.573311e-219 3.786655e-219
[108,] 1.0000000 1.228767e-218 6.143834e-219
[109,] 1.0000000 8.561574e-218 4.280787e-218
[110,] 1.0000000 4.899369e-217 2.449684e-217
[111,] 1.0000000 8.307913e-217 4.153957e-217
[112,] 1.0000000 2.081188e-218 1.040594e-218
[113,] 1.0000000 2.424848e-217 1.212424e-217
[114,] 1.0000000 1.269171e-216 6.345855e-217
[115,] 1.0000000 1.429915e-215 7.149574e-216
[116,] 1.0000000 7.834940e-215 3.917470e-215
[117,] 1.0000000 2.189057e-214 1.094528e-214
[118,] 1.0000000 1.185418e-213 5.927091e-214
[119,] 1.0000000 8.834189e-215 4.417094e-215
[120,] 1.0000000 9.101181e-214 4.550591e-214
[121,] 1.0000000 9.682872e-213 4.841436e-213
[122,] 1.0000000 1.114808e-211 5.574041e-212
[123,] 1.0000000 1.161532e-210 5.807660e-211
[124,] 1.0000000 1.229455e-209 6.147273e-210
[125,] 1.0000000 2.558451e-209 1.279226e-209
[126,] 1.0000000 2.718409e-208 1.359205e-208
[127,] 1.0000000 2.915887e-207 1.457944e-207
[128,] 1.0000000 3.222851e-206 1.611426e-206
[129,] 1.0000000 3.534030e-205 1.767015e-205
[130,] 1.0000000 3.904286e-204 1.952143e-204
[131,] 1.0000000 4.145103e-203 2.072551e-203
[132,] 1.0000000 4.642596e-202 2.321298e-202
[133,] 1.0000000 4.925795e-201 2.462897e-201
[134,] 1.0000000 5.321324e-200 2.660662e-200
[135,] 1.0000000 5.639162e-199 2.819581e-199
[136,] 1.0000000 6.118033e-198 3.059016e-198
[137,] 1.0000000 6.494678e-197 3.247339e-197
[138,] 1.0000000 7.165212e-196 3.582606e-196
[139,] 1.0000000 7.825957e-195 3.912979e-195
[140,] 1.0000000 8.195410e-194 4.097705e-194
[141,] 1.0000000 8.558754e-193 4.279377e-193
[142,] 1.0000000 9.128375e-192 4.564188e-192
[143,] 1.0000000 9.480667e-191 4.740333e-191
[144,] 1.0000000 8.289324e-190 4.144662e-190
[145,] 1.0000000 8.239061e-189 4.119531e-189
[146,] 1.0000000 8.453780e-188 4.226890e-188
[147,] 1.0000000 7.916688e-187 3.958344e-187
[148,] 1.0000000 8.024859e-186 4.012429e-186
[149,] 1.0000000 8.120972e-185 4.060486e-185
[150,] 1.0000000 8.186254e-184 4.093127e-184
[151,] 1.0000000 8.167932e-183 4.083966e-183
[152,] 1.0000000 7.820314e-182 3.910157e-182
[153,] 1.0000000 7.804298e-181 3.902149e-181
[154,] 1.0000000 7.754825e-180 3.877412e-180
[155,] 1.0000000 7.671860e-179 3.835930e-179
[156,] 1.0000000 7.555880e-178 3.777940e-178
[157,] 1.0000000 7.044219e-177 3.522110e-177
[158,] 1.0000000 5.401004e-176 2.700502e-176
[159,] 1.0000000 5.233808e-175 2.616904e-175
[160,] 1.0000000 5.052368e-174 2.526184e-174
[161,] 1.0000000 1.885060e-174 9.425301e-175
[162,] 1.0000000 5.776604e-174 2.888302e-174
[163,] 1.0000000 5.673396e-173 2.836698e-173
[164,] 1.0000000 5.545939e-172 2.772969e-172
[165,] 1.0000000 5.227816e-171 2.613908e-171
[166,] 1.0000000 5.060671e-170 2.530336e-170
[167,] 1.0000000 4.589327e-169 2.294663e-169
[168,] 1.0000000 4.398487e-168 2.199244e-168
[169,] 1.0000000 4.190331e-167 2.095166e-167
[170,] 1.0000000 2.764189e-166 1.382094e-166
[171,] 1.0000000 2.630388e-165 1.315194e-165
[172,] 1.0000000 2.476354e-164 1.238177e-164
[173,] 1.0000000 2.320174e-163 1.160087e-163
[174,] 1.0000000 1.957354e-162 9.786769e-163
[175,] 1.0000000 1.710260e-161 8.551298e-162
[176,] 1.0000000 1.296522e-160 6.482609e-161
[177,] 1.0000000 5.275363e-160 2.637682e-160
[178,] 1.0000000 4.687470e-159 2.343735e-159
[179,] 1.0000000 3.997488e-158 1.998744e-158
[180,] 1.0000000 2.228521e-157 1.114261e-157
[181,] 1.0000000 1.982023e-156 9.910117e-157
[182,] 1.0000000 1.054603e-155 5.273013e-156
[183,] 1.0000000 9.538764e-155 4.769382e-155
[184,] 1.0000000 8.691757e-154 4.345879e-154
[185,] 1.0000000 7.790837e-153 3.895418e-153
[186,] 1.0000000 6.876207e-152 3.438104e-152
[187,] 1.0000000 6.104203e-151 3.052101e-151
[188,] 1.0000000 5.385718e-150 2.692859e-150
[189,] 1.0000000 4.735138e-149 2.367569e-149
[190,] 1.0000000 2.541282e-148 1.270641e-148
[191,] 1.0000000 2.220579e-147 1.110290e-147
[192,] 1.0000000 1.931129e-146 9.655643e-147
[193,] 1.0000000 7.854352e-147 3.927176e-147
[194,] 1.0000000 6.810993e-146 3.405496e-146
[195,] 1.0000000 5.878131e-145 2.939065e-145
[196,] 1.0000000 5.048869e-144 2.524434e-144
[197,] 1.0000000 4.380650e-143 2.190325e-143
[198,] 1.0000000 1.139059e-142 5.695293e-143
[199,] 1.0000000 9.536621e-142 4.768310e-142
[200,] 1.0000000 7.448065e-141 3.724033e-141
[201,] 1.0000000 6.251751e-140 3.125876e-140
[202,] 1.0000000 1.787956e-139 8.939778e-140
[203,] 1.0000000 1.432400e-138 7.161999e-139
[204,] 1.0000000 1.198620e-137 5.993102e-138
[205,] 1.0000000 9.982143e-137 4.991072e-137
[206,] 1.0000000 7.738719e-136 3.869359e-136
[207,] 1.0000000 4.470903e-135 2.235451e-135
[208,] 1.0000000 3.293740e-134 1.646870e-134
[209,] 1.0000000 2.692748e-133 1.346374e-133
[210,] 1.0000000 2.187244e-132 1.093622e-132
[211,] 1.0000000 1.798440e-131 8.992202e-132
[212,] 1.0000000 1.450250e-130 7.251249e-131
[213,] 1.0000000 9.423269e-130 4.711634e-130
[214,] 1.0000000 6.453251e-129 3.226626e-129
[215,] 1.0000000 5.155729e-128 2.577864e-128
[216,] 1.0000000 4.098718e-127 2.049359e-127
[217,] 1.0000000 3.181582e-126 1.590791e-126
[218,] 1.0000000 2.507079e-125 1.253540e-125
[219,] 1.0000000 1.924120e-124 9.620602e-125
[220,] 1.0000000 1.490570e-123 7.452851e-124
[221,] 1.0000000 1.158274e-122 5.791372e-123
[222,] 1.0000000 7.950795e-122 3.975397e-122
[223,] 1.0000000 6.116264e-121 3.058132e-121
[224,] 1.0000000 3.567815e-120 1.783907e-120
[225,] 1.0000000 2.746564e-119 1.373282e-119
[226,] 1.0000000 2.084118e-118 1.042059e-118
[227,] 1.0000000 1.399589e-117 6.997946e-118
[228,] 1.0000000 1.051524e-116 5.257619e-117
[229,] 1.0000000 7.865248e-116 3.932624e-116
[230,] 1.0000000 5.853286e-115 2.926643e-115
[231,] 1.0000000 4.355257e-114 2.177628e-114
[232,] 1.0000000 3.211141e-113 1.605570e-113
[233,] 1.0000000 2.356340e-112 1.178170e-112
[234,] 1.0000000 1.720876e-111 8.604378e-112
[235,] 1.0000000 1.222055e-110 6.110276e-111
[236,] 1.0000000 8.838614e-110 4.419307e-110
[237,] 1.0000000 6.411972e-109 3.205986e-109
[238,] 1.0000000 4.367455e-108 2.183728e-108
[239,] 1.0000000 3.115592e-107 1.557796e-107
[240,] 1.0000000 2.211972e-106 1.105986e-106
[241,] 1.0000000 1.450680e-105 7.253401e-106
[242,] 1.0000000 1.019867e-104 5.099335e-105
[243,] 1.0000000 7.135728e-104 3.567864e-104
[244,] 1.0000000 4.981615e-103 2.490807e-103
[245,] 1.0000000 3.453495e-102 1.726747e-102
[246,] 1.0000000 2.055096e-101 1.027548e-101
[247,] 1.0000000 7.298345e-101 3.649172e-101
[248,] 1.0000000 1.232723e-101 6.163615e-102
[249,] 1.0000000 8.441044e-101 4.220522e-101
[250,] 1.0000000 5.993491e-100 2.996746e-100
[251,] 1.0000000 4.159178e-99 2.079589e-99
[252,] 1.0000000 2.819596e-98 1.409798e-98
[253,] 1.0000000 1.902144e-97 9.510722e-98
[254,] 1.0000000 1.276952e-96 6.384760e-97
[255,] 1.0000000 8.530561e-96 4.265281e-96
[256,] 1.0000000 5.541413e-95 2.770707e-95
[257,] 1.0000000 3.707030e-94 1.853515e-94
[258,] 1.0000000 2.449153e-93 1.224576e-93
[259,] 1.0000000 1.610200e-92 8.051000e-93
[260,] 1.0000000 1.053453e-91 5.267264e-92
[261,] 1.0000000 6.940826e-91 3.470413e-91
[262,] 1.0000000 4.511618e-90 2.255809e-90
[263,] 1.0000000 2.670130e-89 1.335065e-89
[264,] 1.0000000 1.718038e-88 8.590190e-89
[265,] 1.0000000 1.100023e-87 5.500114e-88
[266,] 1.0000000 7.008676e-87 3.504338e-87
[267,] 1.0000000 4.192747e-86 2.096373e-86
[268,] 1.0000000 2.644247e-85 1.322123e-85
[269,] 1.0000000 1.659495e-84 8.297473e-85
[270,] 1.0000000 9.959005e-84 4.979502e-84
[271,] 1.0000000 6.190178e-83 3.095089e-83
[272,] 1.0000000 2.273109e-82 1.136554e-82
[273,] 1.0000000 1.045523e-81 5.227616e-82
[274,] 1.0000000 6.376613e-81 3.188306e-81
[275,] 1.0000000 3.869417e-80 1.934709e-80
[276,] 1.0000000 2.391071e-79 1.195535e-79
[277,] 1.0000000 1.438320e-78 7.191599e-79
[278,] 1.0000000 7.949586e-78 3.974793e-78
[279,] 1.0000000 4.729177e-77 2.364589e-77
[280,] 1.0000000 2.896300e-76 1.448150e-76
[281,] 1.0000000 1.696594e-75 8.482968e-76
[282,] 1.0000000 9.951494e-75 4.975747e-75
[283,] 1.0000000 5.505402e-74 2.752701e-74
[284,] 1.0000000 3.194072e-73 1.597036e-73
[285,] 1.0000000 1.879396e-72 9.396979e-73
[286,] 1.0000000 1.064516e-71 5.322580e-72
[287,] 1.0000000 6.086496e-71 3.043248e-71
[288,] 1.0000000 3.134627e-70 1.567313e-70
[289,] 1.0000000 1.765606e-69 8.828028e-70
[290,] 1.0000000 9.986554e-69 4.993277e-69
[291,] 1.0000000 5.619210e-68 2.809605e-68
[292,] 1.0000000 3.145333e-67 1.572666e-67
[293,] 1.0000000 1.734036e-66 8.670179e-67
[294,] 1.0000000 9.419122e-66 4.709561e-66
[295,] 1.0000000 4.898885e-65 2.449442e-65
[296,] 1.0000000 2.678501e-64 1.339251e-64
[297,] 1.0000000 1.458233e-63 7.291166e-64
[298,] 1.0000000 8.009358e-63 4.004679e-63
[299,] 1.0000000 4.280786e-62 2.140393e-62
[300,] 1.0000000 2.305605e-61 1.152802e-61
[301,] 1.0000000 1.153463e-60 5.767314e-61
[302,] 1.0000000 6.339271e-60 3.169635e-60
[303,] 1.0000000 3.363436e-59 1.681718e-59
[304,] 1.0000000 1.774965e-58 8.874826e-59
[305,] 1.0000000 5.502149e-58 2.751074e-58
[306,] 1.0000000 2.864765e-57 1.432383e-57
[307,] 1.0000000 1.439110e-56 7.195548e-57
[308,] 1.0000000 7.345919e-56 3.672960e-56
[309,] 1.0000000 3.925280e-55 1.962640e-55
[310,] 1.0000000 2.021098e-54 1.010549e-54
[311,] 1.0000000 8.396060e-54 4.198030e-54
[312,] 1.0000000 4.199148e-53 2.099574e-53
[313,] 1.0000000 2.101321e-52 1.050661e-52
[314,] 1.0000000 1.045648e-51 5.228238e-52
[315,] 1.0000000 5.174065e-51 2.587032e-51
[316,] 1.0000000 2.545798e-50 1.272899e-50
[317,] 1.0000000 1.245528e-49 6.227639e-50
[318,] 1.0000000 6.059150e-49 3.029575e-49
[319,] 1.0000000 2.876549e-48 1.438275e-48
[320,] 1.0000000 1.382779e-47 6.913896e-48
[321,] 1.0000000 6.608976e-47 3.304488e-47
[322,] 1.0000000 3.140556e-46 1.570278e-46
[323,] 1.0000000 1.483750e-45 7.418749e-46
[324,] 1.0000000 6.969266e-45 3.484633e-45
[325,] 1.0000000 3.254437e-44 1.627218e-44
[326,] 1.0000000 1.510835e-43 7.554175e-44
[327,] 1.0000000 6.934608e-43 3.467304e-43
[328,] 1.0000000 3.177623e-42 1.588811e-42
[329,] 1.0000000 1.457985e-41 7.289924e-42
[330,] 1.0000000 5.678135e-41 2.839067e-41
[331,] 1.0000000 2.553022e-40 1.276511e-40
[332,] 1.0000000 1.127867e-39 5.639333e-40
[333,] 1.0000000 5.092999e-39 2.546500e-39
[334,] 1.0000000 1.766471e-38 8.832354e-39
[335,] 1.0000000 7.726309e-38 3.863154e-38
[336,] 1.0000000 3.358570e-37 1.679285e-37
[337,] 1.0000000 1.450917e-36 7.254584e-37
[338,] 1.0000000 5.988279e-36 2.994140e-36
[339,] 1.0000000 2.552443e-35 1.276221e-35
[340,] 1.0000000 1.028623e-34 5.143117e-35
[341,] 1.0000000 4.323169e-34 2.161584e-34
[342,] 1.0000000 1.805438e-33 9.027189e-34
[343,] 1.0000000 7.780755e-33 3.890377e-33
[344,] 1.0000000 3.339604e-32 1.669802e-32
[345,] 1.0000000 1.349280e-31 6.746398e-32
[346,] 1.0000000 5.746229e-31 2.873115e-31
[347,] 1.0000000 2.295931e-30 1.147966e-30
[348,] 1.0000000 9.290337e-30 4.645169e-30
[349,] 1.0000000 3.734444e-29 1.867222e-29
[350,] 1.0000000 1.049958e-28 5.249791e-29
[351,] 1.0000000 4.146631e-28 2.073316e-28
[352,] 1.0000000 1.626632e-27 8.133160e-28
[353,] 1.0000000 6.337733e-27 3.168866e-27
[354,] 1.0000000 2.442982e-26 1.221491e-26
[355,] 1.0000000 9.387696e-26 4.693848e-26
[356,] 1.0000000 3.582463e-25 1.791232e-25
[357,] 1.0000000 1.314851e-24 6.574254e-25
[358,] 1.0000000 4.997395e-24 2.498698e-24
[359,] 1.0000000 1.867743e-23 9.338717e-24
[360,] 1.0000000 6.930392e-23 3.465196e-23
[361,] 1.0000000 2.552871e-22 1.276436e-22
[362,] 1.0000000 9.334479e-22 4.667240e-22
[363,] 1.0000000 3.385332e-21 1.692666e-21
[364,] 1.0000000 1.219398e-20 6.096991e-21
[365,] 1.0000000 4.454049e-20 2.227025e-20
[366,] 1.0000000 1.450759e-19 7.253795e-20
[367,] 1.0000000 4.695025e-19 2.347512e-19
[368,] 1.0000000 1.634997e-18 8.174985e-19
[369,] 1.0000000 5.647176e-18 2.823588e-18
[370,] 1.0000000 1.864363e-17 9.321815e-18
[371,] 1.0000000 6.422569e-17 3.211285e-17
[372,] 1.0000000 2.177250e-16 1.088625e-16
[373,] 1.0000000 7.284053e-16 3.642026e-16
[374,] 1.0000000 2.413466e-15 1.206733e-15
[375,] 1.0000000 7.012941e-15 3.506470e-15
[376,] 1.0000000 2.282268e-14 1.141134e-14
[377,] 1.0000000 7.340541e-14 3.670270e-14
[378,] 1.0000000 2.318502e-13 1.159251e-13
[379,] 1.0000000 3.771027e-13 1.885514e-13
[380,] 1.0000000 1.192594e-12 5.962968e-13
[381,] 1.0000000 3.706694e-12 1.853347e-12
[382,] 1.0000000 1.137015e-11 5.685077e-12
[383,] 1.0000000 3.439899e-11 1.719949e-11
[384,] 1.0000000 1.028068e-10 5.140339e-11
[385,] 1.0000000 2.956548e-10 1.478274e-10
[386,] 1.0000000 7.250624e-10 3.625312e-10
[387,] 1.0000000 1.941714e-09 9.708568e-10
[388,] 1.0000000 3.970079e-09 1.985039e-09
[389,] 1.0000000 1.069071e-08 5.345357e-09
[390,] 1.0000000 2.672203e-08 1.336102e-08
[391,] 1.0000000 4.588624e-08 2.294312e-08
[392,] 0.9999999 1.200968e-07 6.004840e-08
[393,] 0.9999998 3.221606e-07 1.610803e-07
[394,] 0.9999996 8.220961e-07 4.110480e-07
[395,] 0.9999989 2.122023e-06 1.061011e-06
[396,] 0.9999973 5.407133e-06 2.703566e-06
[397,] 0.9999940 1.205979e-05 6.029895e-06
[398,] 0.9999861 2.774546e-05 1.387273e-05
[399,] 0.9999666 6.678666e-05 3.339333e-05
[400,] 0.9999216 1.568585e-04 7.842925e-05
[401,] 0.9998208 3.584182e-04 1.792091e-04
[402,] 0.9996009 7.981064e-04 3.990532e-04
[403,] 0.9992772 1.445526e-03 7.227630e-04
[404,] 0.9985135 2.972919e-03 1.486460e-03
[405,] 0.9969994 6.001128e-03 3.000564e-03
[406,] 0.9939295 1.214096e-02 6.070480e-03
[407,] 0.9931498 1.370039e-02 6.850194e-03
[408,] 0.9904063 1.918734e-02 9.593669e-03
[409,] 0.9850379 2.992422e-02 1.496211e-02
[410,] 0.9942056 1.158886e-02 5.794429e-03
[411,] 0.9892813 2.143749e-02 1.071875e-02
[412,] 0.9783294 4.334119e-02 2.167060e-02
[413,] 0.9591700 8.166006e-02 4.083003e-02
[414,] 0.9437181 1.125637e-01 5.628187e-02
[415,] 0.8785829 2.428342e-01 1.214171e-01
[416,] 0.8364419 3.271163e-01 1.635581e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/11lt11291213940.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/freestat/rcomp/tmp/21lt11291213940.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/freestat/rcomp/tmp/31lt11291213940.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/freestat/rcomp/tmp/4cual1291213940.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/freestat/rcomp/tmp/5cual1291213940.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
1378556.5473 2920559.1318 1362238.5911 -2192329.7703 -393334.6417
6 7 8 9 10
113103.7594 -134743.8218 125506.6345 -331198.5858 722401.0795
11 12 13 14 15
-150778.8045 349484.8648 450535.3685 -1131858.9834 -213626.6275
16 17 18 19 20
464402.8104 -241544.4240 -238144.6973 -290334.6560 -586348.5220
21 22 23 24 25
210361.0732 -129541.0024 -117507.2879 -639717.8407 -314694.6713
26 27 28 29 30
-75992.6191 -469171.0226 -331535.4341 358090.3559 -255125.5876
31 32 33 34 35
293763.5421 -82550.2598 263043.2554 136595.9837 -57327.0272
36 37 38 39 40
-567882.3723 40042.8987 161870.7193 -50925.6784 -300909.0262
41 42 43 44 45
-460221.7187 -322989.6990 4491.9591 -205957.5863 -20232.1883
46 47 48 49 50
-47907.6524 -281326.6209 -263843.2642 -49691.7965 -175763.5452
51 52 53 54 55
38946.2512 -86090.2880 -228648.1861 46251.1114 -394276.9548
56 57 58 59 60
20050.5834 328781.9677 -16614.4277 -55038.8645 -201121.6032
61 62 63 64 65
19666.7573 -172340.9734 12490.8868 -169536.6759 -9760.9545
66 67 68 69 70
-495954.0026 -58910.1095 -10662.0331 -98855.4293 -68979.4098
71 72 73 74 75
-114903.1380 8593.0754 -84641.7246 21163.4329 -31231.3166
76 77 78 79 80
-92947.9122 -24768.2825 -67173.3567 -37298.2761 -33869.1040
81 82 83 84 85
33093.2582 -45987.1372 -80479.9940 48939.0844 29295.0480
86 87 88 89 90
-97691.0212 -27066.6374 5406.3288 -112059.9978 36666.8338
91 92 93 94 95
1942.7568 -17490.6982 73675.7476 -169421.7860 55877.4815
96 97 98 99 100
25133.9722 75394.2829 7264.4673 -38613.6608 3857.7440
101 102 103 104 105
-31467.5557 -108390.1822 -163256.4947 -46792.7400 34800.0371
106 107 108 109 110
-48807.2575 -247658.6891 134847.9625 -172469.1537 -87621.3880
111 112 113 114 115
99181.1773 -185446.2654 -161917.6497 -52795.5546 91719.8485
116 117 118 119 120
26251.8232 21592.4992 76197.8356 -127767.3956 19611.0134
121 122 123 124 125
-287804.9497 37409.0695 -102962.4649 87203.2981 18559.8907
126 127 128 129 130
-158629.6053 60143.8504 46653.7379 34002.5266 60021.9485
131 132 133 134 135
51057.9569 61972.4103 59900.0467 35267.8752 38193.5252
136 137 138 139 140
48652.0961 49952.0189 58051.9448 35340.0993 59615.6090
141 142 143 144 145
55620.0338 59534.3411 15639.3696 54667.5314 43853.1344
146 147 148 149 150
47463.3110 59331.1713 51218.8854 61088.6887 58921.4727
151 152 153 154 155
20074.4171 -36526.7328 53876.2080 27468.1387 55996.0625
156 157 158 159 160
58965.4658 58924.8318 53263.1106 65589.9836 58802.9300
161 162 163 164 165
58762.2960 58721.6621 58681.0281 31458.2314 6332.1673
166 167 168 169 170
55987.7124 58518.4923 114695.6319 -74826.4184 58396.5905
171 172 173 174 175
58355.9565 50068.8768 58274.6886 36650.7308 58193.4207
176 177 178 179 180
46844.7364 -70549.0715 57269.4351 58030.8849 57990.2509
181 182 183 184 185
43152.3619 31859.7547 -15320.8909 -3105.2102 48536.2790
186 187 188 189 190
50678.9229 13298.9575 43963.1090 15753.4171 57583.9114
191 192 193 194 195
49547.7785 57502.6435 45590.4319 57421.3756 53673.4067
196 197 198 199 200
57340.1077 -135660.1478 57258.8398 57218.2059 -49957.0442
201 202 203 204 205
57136.9380 57096.3040 57055.6701 41637.2367 648.6942
206 207 208 209 210
50767.8451 48603.9056 56852.5003 52706.0792 39398.4260
211 212 213 214 215
56730.5985 56689.9645 58975.7867 7363.4582 25680.0137
216 217 218 219 220
54442.9648 50033.1890 42915.7504 54321.0630 446.5766
221 222 223 224 225
28493.6951 56283.6250 55144.4390 21870.7797 54739.5314
226 227 228 229 230
25530.9287 44533.0254 56039.8213 39309.5114 55958.5534
231 232 233 234 235
19831.7371 36301.4485 55836.6515 20464.2956 55755.3836
236 237 238 239 240
55740.7560 55674.1157 49228.9269 55592.8478 55552.2139
241 242 243 244 245
55511.5799 46756.1595 55430.3120 49279.9930 -8193.6765
246 247 248 249 250
55308.4102 55267.7762 28572.8824 55186.5083 55145.8744
251 252 253 254 255
34807.2727 55064.6065 36571.6879 -29098.9220 -152852.2784
256 257 258 259 260
54902.0706 43258.1493 23824.0670 54780.1688 54739.5348
261 262 263 264 265
54698.9009 54658.2669 45109.9720 57001.9044 54536.3651
266 267 268 269 270
54495.7311 54455.0972 79382.7272 54373.8293 25568.1954
271 272 273 274 275
54292.5614 54251.9274 54211.2935 26320.8842 54130.0256
276 277 278 279 280
37236.7001 22170.3984 54008.1237 -39688.4814 25050.8307
281 282 283 284 285
53886.2219 53845.5879 25799.7149 53764.3200 36697.3672
286 287 288 289 290
53683.0521 39075.2944 50396.8199 53561.1503 22204.8590
291 292 293 294 295
53479.8824 43685.8372 41310.1326 53357.9805 30720.1589
296 297 298 299 300
63290.0063 53236.0786 53195.4447 53154.8107 43734.0894
301 302 303 304 305
45904.5503 34526.3823 51950.0430 52951.6410 57579.0697
306 307 308 309 310
40309.5053 52829.7391 24349.6867 30803.0719 52707.8373
311 312 313 314 315
52667.2033 -218.1321 52585.9354 38064.6083 44310.8329
316 317 318 319 320
60938.2202 52486.4997 4503.8971 47102.6457 52301.4978
321 322 323 324 325
52260.8638 52220.2299 52179.5959 52138.9620 52098.3280
326 327 328 329 330
47651.5540 52017.0601 51976.4262 51935.7922 51895.1583
331 332 333 334 335
51854.5243 51813.8903 51773.2564 44356.7007 50624.5211
336 337 338 339 340
40584.2542 -7323.1072 51570.0866 48132.7798 42233.0392
341 342 343 344 345
8842.5627 51407.5508 51366.9169 51326.2829 16987.8690
346 347 348 349 350
51245.0150 36348.8650 51163.7471 51123.1132 9654.7960
351 352 353 354 355
47803.9920 40406.2766 39793.7947 43575.4660 50879.3095
356 357 358 359 360
50838.6755 1743.1207 50757.4076 50716.7737 50676.1397
361 362 363 364 365
46459.4107 50594.8718 50554.2379 -25813.1803 39900.4701
366 367 368 369 370
50432.3360 50391.7021 50351.0681 50310.4342 38364.2479
371 372 373 374 375
50229.1663 42559.9748 33435.3504 28591.0951 50066.6304
376 377 378 379 380
50025.9965 23185.1999 36749.7733 -4565.5933 49863.4607
381 382 383 384 385
49822.8267 790.6458 18691.7412 49700.9249 34197.1859
386 387 388 389 390
-123899.8243 56978.0222 20080.9755 44140.9672 24572.7433
391 392 393 394 395
34019.8369 12256.4655 -25124.4773 68245.7197 -103693.0314
396 397 398 399 400
405.4308 -144772.2788 -206035.0405 5987.5583 -60992.2552
401 402 403 404 405
-4167.1290 36067.3258 13277.9864 -56796.4158 -15646.9222
406 407 408 409 410
6941.8966 13587.6916 6583.6352 16283.9769 -68871.9896
411 412 413 414 415
-60102.7556 -168604.4202 -32789.7838 51958.0975 -185299.1439
416 417 418 419 420
-106517.0924 85830.2016 -24856.5535 -58360.9935 -44536.4016
421 422 423 424 425
-51585.4272 -185152.7264 -552875.0594 -465582.0946 -56199.7327
426 427 428 429 430
-259773.3347 -319019.5228 -1152897.9866 -280631.4909 -345807.6229
431
-689241.6941
> postscript(file="/var/www/html/freestat/rcomp/tmp/653ro1291213940.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 1378556.5473 NA
1 2920559.1318 1378556.5473
2 1362238.5911 2920559.1318
3 -2192329.7703 1362238.5911
4 -393334.6417 -2192329.7703
5 113103.7594 -393334.6417
6 -134743.8218 113103.7594
7 125506.6345 -134743.8218
8 -331198.5858 125506.6345
9 722401.0795 -331198.5858
10 -150778.8045 722401.0795
11 349484.8648 -150778.8045
12 450535.3685 349484.8648
13 -1131858.9834 450535.3685
14 -213626.6275 -1131858.9834
15 464402.8104 -213626.6275
16 -241544.4240 464402.8104
17 -238144.6973 -241544.4240
18 -290334.6560 -238144.6973
19 -586348.5220 -290334.6560
20 210361.0732 -586348.5220
21 -129541.0024 210361.0732
22 -117507.2879 -129541.0024
23 -639717.8407 -117507.2879
24 -314694.6713 -639717.8407
25 -75992.6191 -314694.6713
26 -469171.0226 -75992.6191
27 -331535.4341 -469171.0226
28 358090.3559 -331535.4341
29 -255125.5876 358090.3559
30 293763.5421 -255125.5876
31 -82550.2598 293763.5421
32 263043.2554 -82550.2598
33 136595.9837 263043.2554
34 -57327.0272 136595.9837
35 -567882.3723 -57327.0272
36 40042.8987 -567882.3723
37 161870.7193 40042.8987
38 -50925.6784 161870.7193
39 -300909.0262 -50925.6784
40 -460221.7187 -300909.0262
41 -322989.6990 -460221.7187
42 4491.9591 -322989.6990
43 -205957.5863 4491.9591
44 -20232.1883 -205957.5863
45 -47907.6524 -20232.1883
46 -281326.6209 -47907.6524
47 -263843.2642 -281326.6209
48 -49691.7965 -263843.2642
49 -175763.5452 -49691.7965
50 38946.2512 -175763.5452
51 -86090.2880 38946.2512
52 -228648.1861 -86090.2880
53 46251.1114 -228648.1861
54 -394276.9548 46251.1114
55 20050.5834 -394276.9548
56 328781.9677 20050.5834
57 -16614.4277 328781.9677
58 -55038.8645 -16614.4277
59 -201121.6032 -55038.8645
60 19666.7573 -201121.6032
61 -172340.9734 19666.7573
62 12490.8868 -172340.9734
63 -169536.6759 12490.8868
64 -9760.9545 -169536.6759
65 -495954.0026 -9760.9545
66 -58910.1095 -495954.0026
67 -10662.0331 -58910.1095
68 -98855.4293 -10662.0331
69 -68979.4098 -98855.4293
70 -114903.1380 -68979.4098
71 8593.0754 -114903.1380
72 -84641.7246 8593.0754
73 21163.4329 -84641.7246
74 -31231.3166 21163.4329
75 -92947.9122 -31231.3166
76 -24768.2825 -92947.9122
77 -67173.3567 -24768.2825
78 -37298.2761 -67173.3567
79 -33869.1040 -37298.2761
80 33093.2582 -33869.1040
81 -45987.1372 33093.2582
82 -80479.9940 -45987.1372
83 48939.0844 -80479.9940
84 29295.0480 48939.0844
85 -97691.0212 29295.0480
86 -27066.6374 -97691.0212
87 5406.3288 -27066.6374
88 -112059.9978 5406.3288
89 36666.8338 -112059.9978
90 1942.7568 36666.8338
91 -17490.6982 1942.7568
92 73675.7476 -17490.6982
93 -169421.7860 73675.7476
94 55877.4815 -169421.7860
95 25133.9722 55877.4815
96 75394.2829 25133.9722
97 7264.4673 75394.2829
98 -38613.6608 7264.4673
99 3857.7440 -38613.6608
100 -31467.5557 3857.7440
101 -108390.1822 -31467.5557
102 -163256.4947 -108390.1822
103 -46792.7400 -163256.4947
104 34800.0371 -46792.7400
105 -48807.2575 34800.0371
106 -247658.6891 -48807.2575
107 134847.9625 -247658.6891
108 -172469.1537 134847.9625
109 -87621.3880 -172469.1537
110 99181.1773 -87621.3880
111 -185446.2654 99181.1773
112 -161917.6497 -185446.2654
113 -52795.5546 -161917.6497
114 91719.8485 -52795.5546
115 26251.8232 91719.8485
116 21592.4992 26251.8232
117 76197.8356 21592.4992
118 -127767.3956 76197.8356
119 19611.0134 -127767.3956
120 -287804.9497 19611.0134
121 37409.0695 -287804.9497
122 -102962.4649 37409.0695
123 87203.2981 -102962.4649
124 18559.8907 87203.2981
125 -158629.6053 18559.8907
126 60143.8504 -158629.6053
127 46653.7379 60143.8504
128 34002.5266 46653.7379
129 60021.9485 34002.5266
130 51057.9569 60021.9485
131 61972.4103 51057.9569
132 59900.0467 61972.4103
133 35267.8752 59900.0467
134 38193.5252 35267.8752
135 48652.0961 38193.5252
136 49952.0189 48652.0961
137 58051.9448 49952.0189
138 35340.0993 58051.9448
139 59615.6090 35340.0993
140 55620.0338 59615.6090
141 59534.3411 55620.0338
142 15639.3696 59534.3411
143 54667.5314 15639.3696
144 43853.1344 54667.5314
145 47463.3110 43853.1344
146 59331.1713 47463.3110
147 51218.8854 59331.1713
148 61088.6887 51218.8854
149 58921.4727 61088.6887
150 20074.4171 58921.4727
151 -36526.7328 20074.4171
152 53876.2080 -36526.7328
153 27468.1387 53876.2080
154 55996.0625 27468.1387
155 58965.4658 55996.0625
156 58924.8318 58965.4658
157 53263.1106 58924.8318
158 65589.9836 53263.1106
159 58802.9300 65589.9836
160 58762.2960 58802.9300
161 58721.6621 58762.2960
162 58681.0281 58721.6621
163 31458.2314 58681.0281
164 6332.1673 31458.2314
165 55987.7124 6332.1673
166 58518.4923 55987.7124
167 114695.6319 58518.4923
168 -74826.4184 114695.6319
169 58396.5905 -74826.4184
170 58355.9565 58396.5905
171 50068.8768 58355.9565
172 58274.6886 50068.8768
173 36650.7308 58274.6886
174 58193.4207 36650.7308
175 46844.7364 58193.4207
176 -70549.0715 46844.7364
177 57269.4351 -70549.0715
178 58030.8849 57269.4351
179 57990.2509 58030.8849
180 43152.3619 57990.2509
181 31859.7547 43152.3619
182 -15320.8909 31859.7547
183 -3105.2102 -15320.8909
184 48536.2790 -3105.2102
185 50678.9229 48536.2790
186 13298.9575 50678.9229
187 43963.1090 13298.9575
188 15753.4171 43963.1090
189 57583.9114 15753.4171
190 49547.7785 57583.9114
191 57502.6435 49547.7785
192 45590.4319 57502.6435
193 57421.3756 45590.4319
194 53673.4067 57421.3756
195 57340.1077 53673.4067
196 -135660.1478 57340.1077
197 57258.8398 -135660.1478
198 57218.2059 57258.8398
199 -49957.0442 57218.2059
200 57136.9380 -49957.0442
201 57096.3040 57136.9380
202 57055.6701 57096.3040
203 41637.2367 57055.6701
204 648.6942 41637.2367
205 50767.8451 648.6942
206 48603.9056 50767.8451
207 56852.5003 48603.9056
208 52706.0792 56852.5003
209 39398.4260 52706.0792
210 56730.5985 39398.4260
211 56689.9645 56730.5985
212 58975.7867 56689.9645
213 7363.4582 58975.7867
214 25680.0137 7363.4582
215 54442.9648 25680.0137
216 50033.1890 54442.9648
217 42915.7504 50033.1890
218 54321.0630 42915.7504
219 446.5766 54321.0630
220 28493.6951 446.5766
221 56283.6250 28493.6951
222 55144.4390 56283.6250
223 21870.7797 55144.4390
224 54739.5314 21870.7797
225 25530.9287 54739.5314
226 44533.0254 25530.9287
227 56039.8213 44533.0254
228 39309.5114 56039.8213
229 55958.5534 39309.5114
230 19831.7371 55958.5534
231 36301.4485 19831.7371
232 55836.6515 36301.4485
233 20464.2956 55836.6515
234 55755.3836 20464.2956
235 55740.7560 55755.3836
236 55674.1157 55740.7560
237 49228.9269 55674.1157
238 55592.8478 49228.9269
239 55552.2139 55592.8478
240 55511.5799 55552.2139
241 46756.1595 55511.5799
242 55430.3120 46756.1595
243 49279.9930 55430.3120
244 -8193.6765 49279.9930
245 55308.4102 -8193.6765
246 55267.7762 55308.4102
247 28572.8824 55267.7762
248 55186.5083 28572.8824
249 55145.8744 55186.5083
250 34807.2727 55145.8744
251 55064.6065 34807.2727
252 36571.6879 55064.6065
253 -29098.9220 36571.6879
254 -152852.2784 -29098.9220
255 54902.0706 -152852.2784
256 43258.1493 54902.0706
257 23824.0670 43258.1493
258 54780.1688 23824.0670
259 54739.5348 54780.1688
260 54698.9009 54739.5348
261 54658.2669 54698.9009
262 45109.9720 54658.2669
263 57001.9044 45109.9720
264 54536.3651 57001.9044
265 54495.7311 54536.3651
266 54455.0972 54495.7311
267 79382.7272 54455.0972
268 54373.8293 79382.7272
269 25568.1954 54373.8293
270 54292.5614 25568.1954
271 54251.9274 54292.5614
272 54211.2935 54251.9274
273 26320.8842 54211.2935
274 54130.0256 26320.8842
275 37236.7001 54130.0256
276 22170.3984 37236.7001
277 54008.1237 22170.3984
278 -39688.4814 54008.1237
279 25050.8307 -39688.4814
280 53886.2219 25050.8307
281 53845.5879 53886.2219
282 25799.7149 53845.5879
283 53764.3200 25799.7149
284 36697.3672 53764.3200
285 53683.0521 36697.3672
286 39075.2944 53683.0521
287 50396.8199 39075.2944
288 53561.1503 50396.8199
289 22204.8590 53561.1503
290 53479.8824 22204.8590
291 43685.8372 53479.8824
292 41310.1326 43685.8372
293 53357.9805 41310.1326
294 30720.1589 53357.9805
295 63290.0063 30720.1589
296 53236.0786 63290.0063
297 53195.4447 53236.0786
298 53154.8107 53195.4447
299 43734.0894 53154.8107
300 45904.5503 43734.0894
301 34526.3823 45904.5503
302 51950.0430 34526.3823
303 52951.6410 51950.0430
304 57579.0697 52951.6410
305 40309.5053 57579.0697
306 52829.7391 40309.5053
307 24349.6867 52829.7391
308 30803.0719 24349.6867
309 52707.8373 30803.0719
310 52667.2033 52707.8373
311 -218.1321 52667.2033
312 52585.9354 -218.1321
313 38064.6083 52585.9354
314 44310.8329 38064.6083
315 60938.2202 44310.8329
316 52486.4997 60938.2202
317 4503.8971 52486.4997
318 47102.6457 4503.8971
319 52301.4978 47102.6457
320 52260.8638 52301.4978
321 52220.2299 52260.8638
322 52179.5959 52220.2299
323 52138.9620 52179.5959
324 52098.3280 52138.9620
325 47651.5540 52098.3280
326 52017.0601 47651.5540
327 51976.4262 52017.0601
328 51935.7922 51976.4262
329 51895.1583 51935.7922
330 51854.5243 51895.1583
331 51813.8903 51854.5243
332 51773.2564 51813.8903
333 44356.7007 51773.2564
334 50624.5211 44356.7007
335 40584.2542 50624.5211
336 -7323.1072 40584.2542
337 51570.0866 -7323.1072
338 48132.7798 51570.0866
339 42233.0392 48132.7798
340 8842.5627 42233.0392
341 51407.5508 8842.5627
342 51366.9169 51407.5508
343 51326.2829 51366.9169
344 16987.8690 51326.2829
345 51245.0150 16987.8690
346 36348.8650 51245.0150
347 51163.7471 36348.8650
348 51123.1132 51163.7471
349 9654.7960 51123.1132
350 47803.9920 9654.7960
351 40406.2766 47803.9920
352 39793.7947 40406.2766
353 43575.4660 39793.7947
354 50879.3095 43575.4660
355 50838.6755 50879.3095
356 1743.1207 50838.6755
357 50757.4076 1743.1207
358 50716.7737 50757.4076
359 50676.1397 50716.7737
360 46459.4107 50676.1397
361 50594.8718 46459.4107
362 50554.2379 50594.8718
363 -25813.1803 50554.2379
364 39900.4701 -25813.1803
365 50432.3360 39900.4701
366 50391.7021 50432.3360
367 50351.0681 50391.7021
368 50310.4342 50351.0681
369 38364.2479 50310.4342
370 50229.1663 38364.2479
371 42559.9748 50229.1663
372 33435.3504 42559.9748
373 28591.0951 33435.3504
374 50066.6304 28591.0951
375 50025.9965 50066.6304
376 23185.1999 50025.9965
377 36749.7733 23185.1999
378 -4565.5933 36749.7733
379 49863.4607 -4565.5933
380 49822.8267 49863.4607
381 790.6458 49822.8267
382 18691.7412 790.6458
383 49700.9249 18691.7412
384 34197.1859 49700.9249
385 -123899.8243 34197.1859
386 56978.0222 -123899.8243
387 20080.9755 56978.0222
388 44140.9672 20080.9755
389 24572.7433 44140.9672
390 34019.8369 24572.7433
391 12256.4655 34019.8369
392 -25124.4773 12256.4655
393 68245.7197 -25124.4773
394 -103693.0314 68245.7197
395 405.4308 -103693.0314
396 -144772.2788 405.4308
397 -206035.0405 -144772.2788
398 5987.5583 -206035.0405
399 -60992.2552 5987.5583
400 -4167.1290 -60992.2552
401 36067.3258 -4167.1290
402 13277.9864 36067.3258
403 -56796.4158 13277.9864
404 -15646.9222 -56796.4158
405 6941.8966 -15646.9222
406 13587.6916 6941.8966
407 6583.6352 13587.6916
408 16283.9769 6583.6352
409 -68871.9896 16283.9769
410 -60102.7556 -68871.9896
411 -168604.4202 -60102.7556
412 -32789.7838 -168604.4202
413 51958.0975 -32789.7838
414 -185299.1439 51958.0975
415 -106517.0924 -185299.1439
416 85830.2016 -106517.0924
417 -24856.5535 85830.2016
418 -58360.9935 -24856.5535
419 -44536.4016 -58360.9935
420 -51585.4272 -44536.4016
421 -185152.7264 -51585.4272
422 -552875.0594 -185152.7264
423 -465582.0946 -552875.0594
424 -56199.7327 -465582.0946
425 -259773.3347 -56199.7327
426 -319019.5228 -259773.3347
427 -1152897.9866 -319019.5228
428 -280631.4909 -1152897.9866
429 -345807.6229 -280631.4909
430 -689241.6941 -345807.6229
431 NA -689241.6941
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2920559.1318 1378556.5473
[2,] 1362238.5911 2920559.1318
[3,] -2192329.7703 1362238.5911
[4,] -393334.6417 -2192329.7703
[5,] 113103.7594 -393334.6417
[6,] -134743.8218 113103.7594
[7,] 125506.6345 -134743.8218
[8,] -331198.5858 125506.6345
[9,] 722401.0795 -331198.5858
[10,] -150778.8045 722401.0795
[11,] 349484.8648 -150778.8045
[12,] 450535.3685 349484.8648
[13,] -1131858.9834 450535.3685
[14,] -213626.6275 -1131858.9834
[15,] 464402.8104 -213626.6275
[16,] -241544.4240 464402.8104
[17,] -238144.6973 -241544.4240
[18,] -290334.6560 -238144.6973
[19,] -586348.5220 -290334.6560
[20,] 210361.0732 -586348.5220
[21,] -129541.0024 210361.0732
[22,] -117507.2879 -129541.0024
[23,] -639717.8407 -117507.2879
[24,] -314694.6713 -639717.8407
[25,] -75992.6191 -314694.6713
[26,] -469171.0226 -75992.6191
[27,] -331535.4341 -469171.0226
[28,] 358090.3559 -331535.4341
[29,] -255125.5876 358090.3559
[30,] 293763.5421 -255125.5876
[31,] -82550.2598 293763.5421
[32,] 263043.2554 -82550.2598
[33,] 136595.9837 263043.2554
[34,] -57327.0272 136595.9837
[35,] -567882.3723 -57327.0272
[36,] 40042.8987 -567882.3723
[37,] 161870.7193 40042.8987
[38,] -50925.6784 161870.7193
[39,] -300909.0262 -50925.6784
[40,] -460221.7187 -300909.0262
[41,] -322989.6990 -460221.7187
[42,] 4491.9591 -322989.6990
[43,] -205957.5863 4491.9591
[44,] -20232.1883 -205957.5863
[45,] -47907.6524 -20232.1883
[46,] -281326.6209 -47907.6524
[47,] -263843.2642 -281326.6209
[48,] -49691.7965 -263843.2642
[49,] -175763.5452 -49691.7965
[50,] 38946.2512 -175763.5452
[51,] -86090.2880 38946.2512
[52,] -228648.1861 -86090.2880
[53,] 46251.1114 -228648.1861
[54,] -394276.9548 46251.1114
[55,] 20050.5834 -394276.9548
[56,] 328781.9677 20050.5834
[57,] -16614.4277 328781.9677
[58,] -55038.8645 -16614.4277
[59,] -201121.6032 -55038.8645
[60,] 19666.7573 -201121.6032
[61,] -172340.9734 19666.7573
[62,] 12490.8868 -172340.9734
[63,] -169536.6759 12490.8868
[64,] -9760.9545 -169536.6759
[65,] -495954.0026 -9760.9545
[66,] -58910.1095 -495954.0026
[67,] -10662.0331 -58910.1095
[68,] -98855.4293 -10662.0331
[69,] -68979.4098 -98855.4293
[70,] -114903.1380 -68979.4098
[71,] 8593.0754 -114903.1380
[72,] -84641.7246 8593.0754
[73,] 21163.4329 -84641.7246
[74,] -31231.3166 21163.4329
[75,] -92947.9122 -31231.3166
[76,] -24768.2825 -92947.9122
[77,] -67173.3567 -24768.2825
[78,] -37298.2761 -67173.3567
[79,] -33869.1040 -37298.2761
[80,] 33093.2582 -33869.1040
[81,] -45987.1372 33093.2582
[82,] -80479.9940 -45987.1372
[83,] 48939.0844 -80479.9940
[84,] 29295.0480 48939.0844
[85,] -97691.0212 29295.0480
[86,] -27066.6374 -97691.0212
[87,] 5406.3288 -27066.6374
[88,] -112059.9978 5406.3288
[89,] 36666.8338 -112059.9978
[90,] 1942.7568 36666.8338
[91,] -17490.6982 1942.7568
[92,] 73675.7476 -17490.6982
[93,] -169421.7860 73675.7476
[94,] 55877.4815 -169421.7860
[95,] 25133.9722 55877.4815
[96,] 75394.2829 25133.9722
[97,] 7264.4673 75394.2829
[98,] -38613.6608 7264.4673
[99,] 3857.7440 -38613.6608
[100,] -31467.5557 3857.7440
[101,] -108390.1822 -31467.5557
[102,] -163256.4947 -108390.1822
[103,] -46792.7400 -163256.4947
[104,] 34800.0371 -46792.7400
[105,] -48807.2575 34800.0371
[106,] -247658.6891 -48807.2575
[107,] 134847.9625 -247658.6891
[108,] -172469.1537 134847.9625
[109,] -87621.3880 -172469.1537
[110,] 99181.1773 -87621.3880
[111,] -185446.2654 99181.1773
[112,] -161917.6497 -185446.2654
[113,] -52795.5546 -161917.6497
[114,] 91719.8485 -52795.5546
[115,] 26251.8232 91719.8485
[116,] 21592.4992 26251.8232
[117,] 76197.8356 21592.4992
[118,] -127767.3956 76197.8356
[119,] 19611.0134 -127767.3956
[120,] -287804.9497 19611.0134
[121,] 37409.0695 -287804.9497
[122,] -102962.4649 37409.0695
[123,] 87203.2981 -102962.4649
[124,] 18559.8907 87203.2981
[125,] -158629.6053 18559.8907
[126,] 60143.8504 -158629.6053
[127,] 46653.7379 60143.8504
[128,] 34002.5266 46653.7379
[129,] 60021.9485 34002.5266
[130,] 51057.9569 60021.9485
[131,] 61972.4103 51057.9569
[132,] 59900.0467 61972.4103
[133,] 35267.8752 59900.0467
[134,] 38193.5252 35267.8752
[135,] 48652.0961 38193.5252
[136,] 49952.0189 48652.0961
[137,] 58051.9448 49952.0189
[138,] 35340.0993 58051.9448
[139,] 59615.6090 35340.0993
[140,] 55620.0338 59615.6090
[141,] 59534.3411 55620.0338
[142,] 15639.3696 59534.3411
[143,] 54667.5314 15639.3696
[144,] 43853.1344 54667.5314
[145,] 47463.3110 43853.1344
[146,] 59331.1713 47463.3110
[147,] 51218.8854 59331.1713
[148,] 61088.6887 51218.8854
[149,] 58921.4727 61088.6887
[150,] 20074.4171 58921.4727
[151,] -36526.7328 20074.4171
[152,] 53876.2080 -36526.7328
[153,] 27468.1387 53876.2080
[154,] 55996.0625 27468.1387
[155,] 58965.4658 55996.0625
[156,] 58924.8318 58965.4658
[157,] 53263.1106 58924.8318
[158,] 65589.9836 53263.1106
[159,] 58802.9300 65589.9836
[160,] 58762.2960 58802.9300
[161,] 58721.6621 58762.2960
[162,] 58681.0281 58721.6621
[163,] 31458.2314 58681.0281
[164,] 6332.1673 31458.2314
[165,] 55987.7124 6332.1673
[166,] 58518.4923 55987.7124
[167,] 114695.6319 58518.4923
[168,] -74826.4184 114695.6319
[169,] 58396.5905 -74826.4184
[170,] 58355.9565 58396.5905
[171,] 50068.8768 58355.9565
[172,] 58274.6886 50068.8768
[173,] 36650.7308 58274.6886
[174,] 58193.4207 36650.7308
[175,] 46844.7364 58193.4207
[176,] -70549.0715 46844.7364
[177,] 57269.4351 -70549.0715
[178,] 58030.8849 57269.4351
[179,] 57990.2509 58030.8849
[180,] 43152.3619 57990.2509
[181,] 31859.7547 43152.3619
[182,] -15320.8909 31859.7547
[183,] -3105.2102 -15320.8909
[184,] 48536.2790 -3105.2102
[185,] 50678.9229 48536.2790
[186,] 13298.9575 50678.9229
[187,] 43963.1090 13298.9575
[188,] 15753.4171 43963.1090
[189,] 57583.9114 15753.4171
[190,] 49547.7785 57583.9114
[191,] 57502.6435 49547.7785
[192,] 45590.4319 57502.6435
[193,] 57421.3756 45590.4319
[194,] 53673.4067 57421.3756
[195,] 57340.1077 53673.4067
[196,] -135660.1478 57340.1077
[197,] 57258.8398 -135660.1478
[198,] 57218.2059 57258.8398
[199,] -49957.0442 57218.2059
[200,] 57136.9380 -49957.0442
[201,] 57096.3040 57136.9380
[202,] 57055.6701 57096.3040
[203,] 41637.2367 57055.6701
[204,] 648.6942 41637.2367
[205,] 50767.8451 648.6942
[206,] 48603.9056 50767.8451
[207,] 56852.5003 48603.9056
[208,] 52706.0792 56852.5003
[209,] 39398.4260 52706.0792
[210,] 56730.5985 39398.4260
[211,] 56689.9645 56730.5985
[212,] 58975.7867 56689.9645
[213,] 7363.4582 58975.7867
[214,] 25680.0137 7363.4582
[215,] 54442.9648 25680.0137
[216,] 50033.1890 54442.9648
[217,] 42915.7504 50033.1890
[218,] 54321.0630 42915.7504
[219,] 446.5766 54321.0630
[220,] 28493.6951 446.5766
[221,] 56283.6250 28493.6951
[222,] 55144.4390 56283.6250
[223,] 21870.7797 55144.4390
[224,] 54739.5314 21870.7797
[225,] 25530.9287 54739.5314
[226,] 44533.0254 25530.9287
[227,] 56039.8213 44533.0254
[228,] 39309.5114 56039.8213
[229,] 55958.5534 39309.5114
[230,] 19831.7371 55958.5534
[231,] 36301.4485 19831.7371
[232,] 55836.6515 36301.4485
[233,] 20464.2956 55836.6515
[234,] 55755.3836 20464.2956
[235,] 55740.7560 55755.3836
[236,] 55674.1157 55740.7560
[237,] 49228.9269 55674.1157
[238,] 55592.8478 49228.9269
[239,] 55552.2139 55592.8478
[240,] 55511.5799 55552.2139
[241,] 46756.1595 55511.5799
[242,] 55430.3120 46756.1595
[243,] 49279.9930 55430.3120
[244,] -8193.6765 49279.9930
[245,] 55308.4102 -8193.6765
[246,] 55267.7762 55308.4102
[247,] 28572.8824 55267.7762
[248,] 55186.5083 28572.8824
[249,] 55145.8744 55186.5083
[250,] 34807.2727 55145.8744
[251,] 55064.6065 34807.2727
[252,] 36571.6879 55064.6065
[253,] -29098.9220 36571.6879
[254,] -152852.2784 -29098.9220
[255,] 54902.0706 -152852.2784
[256,] 43258.1493 54902.0706
[257,] 23824.0670 43258.1493
[258,] 54780.1688 23824.0670
[259,] 54739.5348 54780.1688
[260,] 54698.9009 54739.5348
[261,] 54658.2669 54698.9009
[262,] 45109.9720 54658.2669
[263,] 57001.9044 45109.9720
[264,] 54536.3651 57001.9044
[265,] 54495.7311 54536.3651
[266,] 54455.0972 54495.7311
[267,] 79382.7272 54455.0972
[268,] 54373.8293 79382.7272
[269,] 25568.1954 54373.8293
[270,] 54292.5614 25568.1954
[271,] 54251.9274 54292.5614
[272,] 54211.2935 54251.9274
[273,] 26320.8842 54211.2935
[274,] 54130.0256 26320.8842
[275,] 37236.7001 54130.0256
[276,] 22170.3984 37236.7001
[277,] 54008.1237 22170.3984
[278,] -39688.4814 54008.1237
[279,] 25050.8307 -39688.4814
[280,] 53886.2219 25050.8307
[281,] 53845.5879 53886.2219
[282,] 25799.7149 53845.5879
[283,] 53764.3200 25799.7149
[284,] 36697.3672 53764.3200
[285,] 53683.0521 36697.3672
[286,] 39075.2944 53683.0521
[287,] 50396.8199 39075.2944
[288,] 53561.1503 50396.8199
[289,] 22204.8590 53561.1503
[290,] 53479.8824 22204.8590
[291,] 43685.8372 53479.8824
[292,] 41310.1326 43685.8372
[293,] 53357.9805 41310.1326
[294,] 30720.1589 53357.9805
[295,] 63290.0063 30720.1589
[296,] 53236.0786 63290.0063
[297,] 53195.4447 53236.0786
[298,] 53154.8107 53195.4447
[299,] 43734.0894 53154.8107
[300,] 45904.5503 43734.0894
[301,] 34526.3823 45904.5503
[302,] 51950.0430 34526.3823
[303,] 52951.6410 51950.0430
[304,] 57579.0697 52951.6410
[305,] 40309.5053 57579.0697
[306,] 52829.7391 40309.5053
[307,] 24349.6867 52829.7391
[308,] 30803.0719 24349.6867
[309,] 52707.8373 30803.0719
[310,] 52667.2033 52707.8373
[311,] -218.1321 52667.2033
[312,] 52585.9354 -218.1321
[313,] 38064.6083 52585.9354
[314,] 44310.8329 38064.6083
[315,] 60938.2202 44310.8329
[316,] 52486.4997 60938.2202
[317,] 4503.8971 52486.4997
[318,] 47102.6457 4503.8971
[319,] 52301.4978 47102.6457
[320,] 52260.8638 52301.4978
[321,] 52220.2299 52260.8638
[322,] 52179.5959 52220.2299
[323,] 52138.9620 52179.5959
[324,] 52098.3280 52138.9620
[325,] 47651.5540 52098.3280
[326,] 52017.0601 47651.5540
[327,] 51976.4262 52017.0601
[328,] 51935.7922 51976.4262
[329,] 51895.1583 51935.7922
[330,] 51854.5243 51895.1583
[331,] 51813.8903 51854.5243
[332,] 51773.2564 51813.8903
[333,] 44356.7007 51773.2564
[334,] 50624.5211 44356.7007
[335,] 40584.2542 50624.5211
[336,] -7323.1072 40584.2542
[337,] 51570.0866 -7323.1072
[338,] 48132.7798 51570.0866
[339,] 42233.0392 48132.7798
[340,] 8842.5627 42233.0392
[341,] 51407.5508 8842.5627
[342,] 51366.9169 51407.5508
[343,] 51326.2829 51366.9169
[344,] 16987.8690 51326.2829
[345,] 51245.0150 16987.8690
[346,] 36348.8650 51245.0150
[347,] 51163.7471 36348.8650
[348,] 51123.1132 51163.7471
[349,] 9654.7960 51123.1132
[350,] 47803.9920 9654.7960
[351,] 40406.2766 47803.9920
[352,] 39793.7947 40406.2766
[353,] 43575.4660 39793.7947
[354,] 50879.3095 43575.4660
[355,] 50838.6755 50879.3095
[356,] 1743.1207 50838.6755
[357,] 50757.4076 1743.1207
[358,] 50716.7737 50757.4076
[359,] 50676.1397 50716.7737
[360,] 46459.4107 50676.1397
[361,] 50594.8718 46459.4107
[362,] 50554.2379 50594.8718
[363,] -25813.1803 50554.2379
[364,] 39900.4701 -25813.1803
[365,] 50432.3360 39900.4701
[366,] 50391.7021 50432.3360
[367,] 50351.0681 50391.7021
[368,] 50310.4342 50351.0681
[369,] 38364.2479 50310.4342
[370,] 50229.1663 38364.2479
[371,] 42559.9748 50229.1663
[372,] 33435.3504 42559.9748
[373,] 28591.0951 33435.3504
[374,] 50066.6304 28591.0951
[375,] 50025.9965 50066.6304
[376,] 23185.1999 50025.9965
[377,] 36749.7733 23185.1999
[378,] -4565.5933 36749.7733
[379,] 49863.4607 -4565.5933
[380,] 49822.8267 49863.4607
[381,] 790.6458 49822.8267
[382,] 18691.7412 790.6458
[383,] 49700.9249 18691.7412
[384,] 34197.1859 49700.9249
[385,] -123899.8243 34197.1859
[386,] 56978.0222 -123899.8243
[387,] 20080.9755 56978.0222
[388,] 44140.9672 20080.9755
[389,] 24572.7433 44140.9672
[390,] 34019.8369 24572.7433
[391,] 12256.4655 34019.8369
[392,] -25124.4773 12256.4655
[393,] 68245.7197 -25124.4773
[394,] -103693.0314 68245.7197
[395,] 405.4308 -103693.0314
[396,] -144772.2788 405.4308
[397,] -206035.0405 -144772.2788
[398,] 5987.5583 -206035.0405
[399,] -60992.2552 5987.5583
[400,] -4167.1290 -60992.2552
[401,] 36067.3258 -4167.1290
[402,] 13277.9864 36067.3258
[403,] -56796.4158 13277.9864
[404,] -15646.9222 -56796.4158
[405,] 6941.8966 -15646.9222
[406,] 13587.6916 6941.8966
[407,] 6583.6352 13587.6916
[408,] 16283.9769 6583.6352
[409,] -68871.9896 16283.9769
[410,] -60102.7556 -68871.9896
[411,] -168604.4202 -60102.7556
[412,] -32789.7838 -168604.4202
[413,] 51958.0975 -32789.7838
[414,] -185299.1439 51958.0975
[415,] -106517.0924 -185299.1439
[416,] 85830.2016 -106517.0924
[417,] -24856.5535 85830.2016
[418,] -58360.9935 -24856.5535
[419,] -44536.4016 -58360.9935
[420,] -51585.4272 -44536.4016
[421,] -185152.7264 -51585.4272
[422,] -552875.0594 -185152.7264
[423,] -465582.0946 -552875.0594
[424,] -56199.7327 -465582.0946
[425,] -259773.3347 -56199.7327
[426,] -319019.5228 -259773.3347
[427,] -1152897.9866 -319019.5228
[428,] -280631.4909 -1152897.9866
[429,] -345807.6229 -280631.4909
[430,] -689241.6941 -345807.6229
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2920559.1318 1378556.5473
2 1362238.5911 2920559.1318
3 -2192329.7703 1362238.5911
4 -393334.6417 -2192329.7703
5 113103.7594 -393334.6417
6 -134743.8218 113103.7594
7 125506.6345 -134743.8218
8 -331198.5858 125506.6345
9 722401.0795 -331198.5858
10 -150778.8045 722401.0795
11 349484.8648 -150778.8045
12 450535.3685 349484.8648
13 -1131858.9834 450535.3685
14 -213626.6275 -1131858.9834
15 464402.8104 -213626.6275
16 -241544.4240 464402.8104
17 -238144.6973 -241544.4240
18 -290334.6560 -238144.6973
19 -586348.5220 -290334.6560
20 210361.0732 -586348.5220
21 -129541.0024 210361.0732
22 -117507.2879 -129541.0024
23 -639717.8407 -117507.2879
24 -314694.6713 -639717.8407
25 -75992.6191 -314694.6713
26 -469171.0226 -75992.6191
27 -331535.4341 -469171.0226
28 358090.3559 -331535.4341
29 -255125.5876 358090.3559
30 293763.5421 -255125.5876
31 -82550.2598 293763.5421
32 263043.2554 -82550.2598
33 136595.9837 263043.2554
34 -57327.0272 136595.9837
35 -567882.3723 -57327.0272
36 40042.8987 -567882.3723
37 161870.7193 40042.8987
38 -50925.6784 161870.7193
39 -300909.0262 -50925.6784
40 -460221.7187 -300909.0262
41 -322989.6990 -460221.7187
42 4491.9591 -322989.6990
43 -205957.5863 4491.9591
44 -20232.1883 -205957.5863
45 -47907.6524 -20232.1883
46 -281326.6209 -47907.6524
47 -263843.2642 -281326.6209
48 -49691.7965 -263843.2642
49 -175763.5452 -49691.7965
50 38946.2512 -175763.5452
51 -86090.2880 38946.2512
52 -228648.1861 -86090.2880
53 46251.1114 -228648.1861
54 -394276.9548 46251.1114
55 20050.5834 -394276.9548
56 328781.9677 20050.5834
57 -16614.4277 328781.9677
58 -55038.8645 -16614.4277
59 -201121.6032 -55038.8645
60 19666.7573 -201121.6032
61 -172340.9734 19666.7573
62 12490.8868 -172340.9734
63 -169536.6759 12490.8868
64 -9760.9545 -169536.6759
65 -495954.0026 -9760.9545
66 -58910.1095 -495954.0026
67 -10662.0331 -58910.1095
68 -98855.4293 -10662.0331
69 -68979.4098 -98855.4293
70 -114903.1380 -68979.4098
71 8593.0754 -114903.1380
72 -84641.7246 8593.0754
73 21163.4329 -84641.7246
74 -31231.3166 21163.4329
75 -92947.9122 -31231.3166
76 -24768.2825 -92947.9122
77 -67173.3567 -24768.2825
78 -37298.2761 -67173.3567
79 -33869.1040 -37298.2761
80 33093.2582 -33869.1040
81 -45987.1372 33093.2582
82 -80479.9940 -45987.1372
83 48939.0844 -80479.9940
84 29295.0480 48939.0844
85 -97691.0212 29295.0480
86 -27066.6374 -97691.0212
87 5406.3288 -27066.6374
88 -112059.9978 5406.3288
89 36666.8338 -112059.9978
90 1942.7568 36666.8338
91 -17490.6982 1942.7568
92 73675.7476 -17490.6982
93 -169421.7860 73675.7476
94 55877.4815 -169421.7860
95 25133.9722 55877.4815
96 75394.2829 25133.9722
97 7264.4673 75394.2829
98 -38613.6608 7264.4673
99 3857.7440 -38613.6608
100 -31467.5557 3857.7440
101 -108390.1822 -31467.5557
102 -163256.4947 -108390.1822
103 -46792.7400 -163256.4947
104 34800.0371 -46792.7400
105 -48807.2575 34800.0371
106 -247658.6891 -48807.2575
107 134847.9625 -247658.6891
108 -172469.1537 134847.9625
109 -87621.3880 -172469.1537
110 99181.1773 -87621.3880
111 -185446.2654 99181.1773
112 -161917.6497 -185446.2654
113 -52795.5546 -161917.6497
114 91719.8485 -52795.5546
115 26251.8232 91719.8485
116 21592.4992 26251.8232
117 76197.8356 21592.4992
118 -127767.3956 76197.8356
119 19611.0134 -127767.3956
120 -287804.9497 19611.0134
121 37409.0695 -287804.9497
122 -102962.4649 37409.0695
123 87203.2981 -102962.4649
124 18559.8907 87203.2981
125 -158629.6053 18559.8907
126 60143.8504 -158629.6053
127 46653.7379 60143.8504
128 34002.5266 46653.7379
129 60021.9485 34002.5266
130 51057.9569 60021.9485
131 61972.4103 51057.9569
132 59900.0467 61972.4103
133 35267.8752 59900.0467
134 38193.5252 35267.8752
135 48652.0961 38193.5252
136 49952.0189 48652.0961
137 58051.9448 49952.0189
138 35340.0993 58051.9448
139 59615.6090 35340.0993
140 55620.0338 59615.6090
141 59534.3411 55620.0338
142 15639.3696 59534.3411
143 54667.5314 15639.3696
144 43853.1344 54667.5314
145 47463.3110 43853.1344
146 59331.1713 47463.3110
147 51218.8854 59331.1713
148 61088.6887 51218.8854
149 58921.4727 61088.6887
150 20074.4171 58921.4727
151 -36526.7328 20074.4171
152 53876.2080 -36526.7328
153 27468.1387 53876.2080
154 55996.0625 27468.1387
155 58965.4658 55996.0625
156 58924.8318 58965.4658
157 53263.1106 58924.8318
158 65589.9836 53263.1106
159 58802.9300 65589.9836
160 58762.2960 58802.9300
161 58721.6621 58762.2960
162 58681.0281 58721.6621
163 31458.2314 58681.0281
164 6332.1673 31458.2314
165 55987.7124 6332.1673
166 58518.4923 55987.7124
167 114695.6319 58518.4923
168 -74826.4184 114695.6319
169 58396.5905 -74826.4184
170 58355.9565 58396.5905
171 50068.8768 58355.9565
172 58274.6886 50068.8768
173 36650.7308 58274.6886
174 58193.4207 36650.7308
175 46844.7364 58193.4207
176 -70549.0715 46844.7364
177 57269.4351 -70549.0715
178 58030.8849 57269.4351
179 57990.2509 58030.8849
180 43152.3619 57990.2509
181 31859.7547 43152.3619
182 -15320.8909 31859.7547
183 -3105.2102 -15320.8909
184 48536.2790 -3105.2102
185 50678.9229 48536.2790
186 13298.9575 50678.9229
187 43963.1090 13298.9575
188 15753.4171 43963.1090
189 57583.9114 15753.4171
190 49547.7785 57583.9114
191 57502.6435 49547.7785
192 45590.4319 57502.6435
193 57421.3756 45590.4319
194 53673.4067 57421.3756
195 57340.1077 53673.4067
196 -135660.1478 57340.1077
197 57258.8398 -135660.1478
198 57218.2059 57258.8398
199 -49957.0442 57218.2059
200 57136.9380 -49957.0442
201 57096.3040 57136.9380
202 57055.6701 57096.3040
203 41637.2367 57055.6701
204 648.6942 41637.2367
205 50767.8451 648.6942
206 48603.9056 50767.8451
207 56852.5003 48603.9056
208 52706.0792 56852.5003
209 39398.4260 52706.0792
210 56730.5985 39398.4260
211 56689.9645 56730.5985
212 58975.7867 56689.9645
213 7363.4582 58975.7867
214 25680.0137 7363.4582
215 54442.9648 25680.0137
216 50033.1890 54442.9648
217 42915.7504 50033.1890
218 54321.0630 42915.7504
219 446.5766 54321.0630
220 28493.6951 446.5766
221 56283.6250 28493.6951
222 55144.4390 56283.6250
223 21870.7797 55144.4390
224 54739.5314 21870.7797
225 25530.9287 54739.5314
226 44533.0254 25530.9287
227 56039.8213 44533.0254
228 39309.5114 56039.8213
229 55958.5534 39309.5114
230 19831.7371 55958.5534
231 36301.4485 19831.7371
232 55836.6515 36301.4485
233 20464.2956 55836.6515
234 55755.3836 20464.2956
235 55740.7560 55755.3836
236 55674.1157 55740.7560
237 49228.9269 55674.1157
238 55592.8478 49228.9269
239 55552.2139 55592.8478
240 55511.5799 55552.2139
241 46756.1595 55511.5799
242 55430.3120 46756.1595
243 49279.9930 55430.3120
244 -8193.6765 49279.9930
245 55308.4102 -8193.6765
246 55267.7762 55308.4102
247 28572.8824 55267.7762
248 55186.5083 28572.8824
249 55145.8744 55186.5083
250 34807.2727 55145.8744
251 55064.6065 34807.2727
252 36571.6879 55064.6065
253 -29098.9220 36571.6879
254 -152852.2784 -29098.9220
255 54902.0706 -152852.2784
256 43258.1493 54902.0706
257 23824.0670 43258.1493
258 54780.1688 23824.0670
259 54739.5348 54780.1688
260 54698.9009 54739.5348
261 54658.2669 54698.9009
262 45109.9720 54658.2669
263 57001.9044 45109.9720
264 54536.3651 57001.9044
265 54495.7311 54536.3651
266 54455.0972 54495.7311
267 79382.7272 54455.0972
268 54373.8293 79382.7272
269 25568.1954 54373.8293
270 54292.5614 25568.1954
271 54251.9274 54292.5614
272 54211.2935 54251.9274
273 26320.8842 54211.2935
274 54130.0256 26320.8842
275 37236.7001 54130.0256
276 22170.3984 37236.7001
277 54008.1237 22170.3984
278 -39688.4814 54008.1237
279 25050.8307 -39688.4814
280 53886.2219 25050.8307
281 53845.5879 53886.2219
282 25799.7149 53845.5879
283 53764.3200 25799.7149
284 36697.3672 53764.3200
285 53683.0521 36697.3672
286 39075.2944 53683.0521
287 50396.8199 39075.2944
288 53561.1503 50396.8199
289 22204.8590 53561.1503
290 53479.8824 22204.8590
291 43685.8372 53479.8824
292 41310.1326 43685.8372
293 53357.9805 41310.1326
294 30720.1589 53357.9805
295 63290.0063 30720.1589
296 53236.0786 63290.0063
297 53195.4447 53236.0786
298 53154.8107 53195.4447
299 43734.0894 53154.8107
300 45904.5503 43734.0894
301 34526.3823 45904.5503
302 51950.0430 34526.3823
303 52951.6410 51950.0430
304 57579.0697 52951.6410
305 40309.5053 57579.0697
306 52829.7391 40309.5053
307 24349.6867 52829.7391
308 30803.0719 24349.6867
309 52707.8373 30803.0719
310 52667.2033 52707.8373
311 -218.1321 52667.2033
312 52585.9354 -218.1321
313 38064.6083 52585.9354
314 44310.8329 38064.6083
315 60938.2202 44310.8329
316 52486.4997 60938.2202
317 4503.8971 52486.4997
318 47102.6457 4503.8971
319 52301.4978 47102.6457
320 52260.8638 52301.4978
321 52220.2299 52260.8638
322 52179.5959 52220.2299
323 52138.9620 52179.5959
324 52098.3280 52138.9620
325 47651.5540 52098.3280
326 52017.0601 47651.5540
327 51976.4262 52017.0601
328 51935.7922 51976.4262
329 51895.1583 51935.7922
330 51854.5243 51895.1583
331 51813.8903 51854.5243
332 51773.2564 51813.8903
333 44356.7007 51773.2564
334 50624.5211 44356.7007
335 40584.2542 50624.5211
336 -7323.1072 40584.2542
337 51570.0866 -7323.1072
338 48132.7798 51570.0866
339 42233.0392 48132.7798
340 8842.5627 42233.0392
341 51407.5508 8842.5627
342 51366.9169 51407.5508
343 51326.2829 51366.9169
344 16987.8690 51326.2829
345 51245.0150 16987.8690
346 36348.8650 51245.0150
347 51163.7471 36348.8650
348 51123.1132 51163.7471
349 9654.7960 51123.1132
350 47803.9920 9654.7960
351 40406.2766 47803.9920
352 39793.7947 40406.2766
353 43575.4660 39793.7947
354 50879.3095 43575.4660
355 50838.6755 50879.3095
356 1743.1207 50838.6755
357 50757.4076 1743.1207
358 50716.7737 50757.4076
359 50676.1397 50716.7737
360 46459.4107 50676.1397
361 50594.8718 46459.4107
362 50554.2379 50594.8718
363 -25813.1803 50554.2379
364 39900.4701 -25813.1803
365 50432.3360 39900.4701
366 50391.7021 50432.3360
367 50351.0681 50391.7021
368 50310.4342 50351.0681
369 38364.2479 50310.4342
370 50229.1663 38364.2479
371 42559.9748 50229.1663
372 33435.3504 42559.9748
373 28591.0951 33435.3504
374 50066.6304 28591.0951
375 50025.9965 50066.6304
376 23185.1999 50025.9965
377 36749.7733 23185.1999
378 -4565.5933 36749.7733
379 49863.4607 -4565.5933
380 49822.8267 49863.4607
381 790.6458 49822.8267
382 18691.7412 790.6458
383 49700.9249 18691.7412
384 34197.1859 49700.9249
385 -123899.8243 34197.1859
386 56978.0222 -123899.8243
387 20080.9755 56978.0222
388 44140.9672 20080.9755
389 24572.7433 44140.9672
390 34019.8369 24572.7433
391 12256.4655 34019.8369
392 -25124.4773 12256.4655
393 68245.7197 -25124.4773
394 -103693.0314 68245.7197
395 405.4308 -103693.0314
396 -144772.2788 405.4308
397 -206035.0405 -144772.2788
398 5987.5583 -206035.0405
399 -60992.2552 5987.5583
400 -4167.1290 -60992.2552
401 36067.3258 -4167.1290
402 13277.9864 36067.3258
403 -56796.4158 13277.9864
404 -15646.9222 -56796.4158
405 6941.8966 -15646.9222
406 13587.6916 6941.8966
407 6583.6352 13587.6916
408 16283.9769 6583.6352
409 -68871.9896 16283.9769
410 -60102.7556 -68871.9896
411 -168604.4202 -60102.7556
412 -32789.7838 -168604.4202
413 51958.0975 -32789.7838
414 -185299.1439 51958.0975
415 -106517.0924 -185299.1439
416 85830.2016 -106517.0924
417 -24856.5535 85830.2016
418 -58360.9935 -24856.5535
419 -44536.4016 -58360.9935
420 -51585.4272 -44536.4016
421 -185152.7264 -51585.4272
422 -552875.0594 -185152.7264
423 -465582.0946 -552875.0594
424 -56199.7327 -465582.0946
425 -259773.3347 -56199.7327
426 -319019.5228 -259773.3347
427 -1152897.9866 -319019.5228
428 -280631.4909 -1152897.9866
429 -345807.6229 -280631.4909
430 -689241.6941 -345807.6229
> 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/freestat/rcomp/tmp/753ro1291213940.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/freestat/rcomp/tmp/8fd8r1291213940.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/freestat/rcomp/tmp/9fd8r1291213940.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/freestat/rcomp/tmp/1084pu1291213940.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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/11tmo01291213940.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/freestat/rcomp/tmp/12f5m61291213940.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/freestat/rcomp/tmp/13tfkf1291213940.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/freestat/rcomp/tmp/14ef1k1291213940.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/freestat/rcomp/tmp/1576051291213940.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/freestat/rcomp/tmp/163yye1291213940.tab")
+ }
>
> try(system("convert tmp/11lt11291213940.ps tmp/11lt11291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/21lt11291213940.ps tmp/21lt11291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/31lt11291213940.ps tmp/31lt11291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/4cual1291213940.ps tmp/4cual1291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/5cual1291213940.ps tmp/5cual1291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/653ro1291213940.ps tmp/653ro1291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/753ro1291213940.ps tmp/753ro1291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/8fd8r1291213940.ps tmp/8fd8r1291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/9fd8r1291213940.ps tmp/9fd8r1291213940.png",intern=TRUE))
character(0)
> try(system("convert tmp/1084pu1291213940.ps tmp/1084pu1291213940.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
13.656 3.204 14.140