R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(162556
+ ,807
+ ,213118
+ ,6282154
+ ,29790
+ ,444
+ ,81767
+ ,4321023
+ ,87550
+ ,412
+ ,153198
+ ,4111912
+ ,84738
+ ,428
+ ,-26007
+ ,223193
+ ,54660
+ ,315
+ ,126942
+ ,1491348
+ ,42634
+ ,168
+ ,157214
+ ,1629616
+ ,40949
+ ,263
+ ,129352
+ ,1398893
+ ,45187
+ ,267
+ ,234817
+ ,1926517
+ ,37704
+ ,228
+ ,60448
+ ,983660
+ ,16275
+ ,129
+ ,47818
+ ,1443586
+ ,25830
+ ,104
+ ,245546
+ ,1073089
+ ,12679
+ ,122
+ ,48020
+ ,984885
+ ,18014
+ ,393
+ ,-1710
+ ,1405225
+ ,43556
+ ,190
+ ,32648
+ ,227132
+ ,24811
+ ,280
+ ,95350
+ ,929118
+ ,6575
+ ,63
+ ,151352
+ ,1071292
+ ,7123
+ ,102
+ ,288170
+ ,638830
+ ,21950
+ ,265
+ ,114337
+ ,856956
+ ,37597
+ ,234
+ ,37884
+ ,992426
+ ,17821
+ ,277
+ ,122844
+ ,444477
+ ,12988
+ ,73
+ ,82340
+ ,857217
+ ,22330
+ ,67
+ ,79801
+ ,711969
+ ,13326
+ ,103
+ ,165548
+ ,702380
+ ,16189
+ ,290
+ ,116384
+ ,358589
+ ,7146
+ ,83
+ ,134028
+ ,297978
+ ,15824
+ ,56
+ ,63838
+ ,585715
+ ,27664
+ ,236
+ ,74996
+ ,657954
+ ,11920
+ ,73
+ ,31080
+ ,209458
+ ,8568
+ ,34
+ ,32168
+ ,786690
+ ,14416
+ ,139
+ ,49857
+ ,439798
+ ,3369
+ ,26
+ ,87161
+ ,688779
+ ,11819
+ ,70
+ ,106113
+ ,574339
+ ,6984
+ ,40
+ ,80570
+ ,741409
+ ,4519
+ ,42
+ ,102129
+ ,597793
+ ,2220
+ ,12
+ ,301670
+ ,644190
+ ,18562
+ ,211
+ ,102313
+ ,377934
+ ,10327
+ ,74
+ ,88577
+ ,640273
+ ,5336
+ ,80
+ ,112477
+ ,697458
+ ,2365
+ ,83
+ ,191778
+ ,550608
+ ,4069
+ ,131
+ ,79804
+ ,207393
+ ,8636
+ ,203
+ ,128294
+ ,301607
+ ,13718
+ ,56
+ ,96448
+ ,345783
+ ,4525
+ ,89
+ ,93811
+ ,501749
+ ,6869
+ ,88
+ ,117520
+ ,379983
+ ,4628
+ ,39
+ ,69159
+ ,387475
+ ,3689
+ ,25
+ ,101792
+ ,377305
+ ,4891
+ ,49
+ ,210568
+ ,370837
+ ,7489
+ ,149
+ ,136996
+ ,430866
+ ,4901
+ ,58
+ ,121920
+ ,469107
+ ,2284
+ ,41
+ ,76403
+ ,194493
+ ,3160
+ ,90
+ ,108094
+ ,530670
+ ,4150
+ ,136
+ ,134759
+ ,518365
+ ,7285
+ ,97
+ ,188873
+ ,491303
+ ,1134
+ ,63
+ ,146216
+ ,527021
+ ,4658
+ ,114
+ ,156608
+ ,233773
+ ,2384
+ ,77
+ ,61348
+ ,405972
+ ,3748
+ ,6
+ ,50350
+ ,652925
+ ,5371
+ ,47
+ ,87720
+ ,446211
+ ,1285
+ ,51
+ ,99489
+ ,341340
+ ,9327
+ ,85
+ ,87419
+ ,387699
+ ,5565
+ ,43
+ ,94355
+ ,493408
+ ,1528
+ ,32
+ ,60326
+ ,146494
+ ,3122
+ ,25
+ ,94670
+ ,414462
+ ,7561
+ ,77
+ ,82425
+ ,364304
+ ,2675
+ ,54
+ ,59017
+ ,355178
+ ,13253
+ ,251
+ ,90829
+ ,357760
+ ,880
+ ,15
+ ,80791
+ ,261216
+ ,2053
+ ,44
+ ,100423
+ ,397144
+ ,1424
+ ,73
+ ,131116
+ ,374943
+ ,4036
+ ,85
+ ,100269
+ ,424898
+ ,3045
+ ,49
+ ,27330
+ ,202055
+ ,5119
+ ,38
+ ,39039
+ ,378525
+ ,1431
+ ,35
+ ,106885
+ ,310768
+ ,554
+ ,9
+ ,79285
+ ,325738
+ ,1975
+ ,34
+ ,118881
+ ,394510
+ ,1765
+ ,20
+ ,77623
+ ,247060
+ ,1012
+ ,29
+ ,114768
+ ,368078
+ ,810
+ ,11
+ ,74015
+ ,236761
+ ,1280
+ ,52
+ ,69465
+ ,312378
+ ,666
+ ,13
+ ,117869
+ ,339836
+ ,1380
+ ,29
+ ,60982
+ ,347385
+ ,4677
+ ,66
+ ,90131
+ ,426280
+ ,876
+ ,33
+ ,138971
+ ,352850
+ ,814
+ ,15
+ ,39625
+ ,301881
+ ,514
+ ,15
+ ,102725
+ ,377516
+ ,5692
+ ,68
+ ,64239
+ ,357312
+ ,3642
+ ,100
+ ,90262
+ ,458343
+ ,540
+ ,13
+ ,103960
+ ,354228
+ ,2099
+ ,45
+ ,106611
+ ,308636
+ ,567
+ ,14
+ ,103345
+ ,386212
+ ,2001
+ ,36
+ ,95551
+ ,393343
+ ,2949
+ ,40
+ ,82903
+ ,378509
+ ,2253
+ ,68
+ ,63593
+ ,452469
+ ,6533
+ ,29
+ ,126910
+ ,364839
+ ,1889
+ ,43
+ ,37527
+ ,358649
+ ,3055
+ ,30
+ ,60247
+ ,376641
+ ,272
+ ,9
+ ,112995
+ ,429112
+ ,1414
+ ,22
+ ,70184
+ ,330546
+ ,2564
+ ,19
+ ,130140
+ ,403560
+ ,1383
+ ,9
+ ,73221
+ ,317892
+ ,1261
+ ,31
+ ,76114
+ ,307528
+ ,975
+ ,19
+ ,90534
+ ,235133
+ ,3366
+ ,55
+ ,108479
+ ,299243
+ ,576
+ ,8
+ ,113761
+ ,314073
+ ,1686
+ ,28
+ ,68696
+ ,368186
+ ,746
+ ,29
+ ,71561
+ ,269661
+ ,3192
+ ,48
+ ,59831
+ ,125390
+ ,2045
+ ,16
+ ,97890
+ ,510834
+ ,5702
+ ,47
+ ,101481
+ ,321896
+ ,1932
+ ,20
+ ,72954
+ ,249898
+ ,936
+ ,22
+ ,67939
+ ,408881
+ ,3437
+ ,33
+ ,48022
+ ,158492
+ ,5131
+ ,44
+ ,86111
+ ,292154
+ ,2397
+ ,13
+ ,74020
+ ,289513
+ ,1389
+ ,6
+ ,57530
+ ,378049
+ ,1503
+ ,35
+ ,56364
+ ,343466
+ ,402
+ ,8
+ ,84990
+ ,332743
+ ,2239
+ ,17
+ ,88590
+ ,442882
+ ,2234
+ ,11
+ ,77200
+ ,214215
+ ,837
+ ,21
+ ,61262
+ ,315688
+ ,10579
+ ,92
+ ,110309
+ ,375195
+ ,875
+ ,12
+ ,67000
+ ,334280
+ ,1585
+ ,112
+ ,93099
+ ,355864
+ ,1659
+ ,25
+ ,107577
+ ,480382
+ ,2647
+ ,17
+ ,62920
+ ,353058
+ ,3294
+ ,23
+ ,75832
+ ,217193
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,94
+ ,10
+ ,60793
+ ,314533
+ ,422
+ ,23
+ ,57935
+ ,318056
+ ,0
+ ,0
+ ,60720
+ ,315380
+ ,34
+ ,7
+ ,60630
+ ,314353
+ ,1558
+ ,25
+ ,55637
+ ,369448
+ ,0
+ ,1
+ ,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('Costs'
+ ,'Orders'
+ ,'Dividends'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(4,431),dimnames=list(c('Costs','Orders','Dividends','Wealth'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '4'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Costs Orders Dividends
1 6282154 162556 807 213118
2 4321023 29790 444 81767
3 4111912 87550 412 153198
4 223193 84738 428 -26007
5 1491348 54660 315 126942
6 1629616 42634 168 157214
7 1398893 40949 263 129352
8 1926517 45187 267 234817
9 983660 37704 228 60448
10 1443586 16275 129 47818
11 1073089 25830 104 245546
12 984885 12679 122 48020
13 1405225 18014 393 -1710
14 227132 43556 190 32648
15 929118 24811 280 95350
16 1071292 6575 63 151352
17 638830 7123 102 288170
18 856956 21950 265 114337
19 992426 37597 234 37884
20 444477 17821 277 122844
21 857217 12988 73 82340
22 711969 22330 67 79801
23 702380 13326 103 165548
24 358589 16189 290 116384
25 297978 7146 83 134028
26 585715 15824 56 63838
27 657954 27664 236 74996
28 209458 11920 73 31080
29 786690 8568 34 32168
30 439798 14416 139 49857
31 688779 3369 26 87161
32 574339 11819 70 106113
33 741409 6984 40 80570
34 597793 4519 42 102129
35 644190 2220 12 301670
36 377934 18562 211 102313
37 640273 10327 74 88577
38 697458 5336 80 112477
39 550608 2365 83 191778
40 207393 4069 131 79804
41 301607 8636 203 128294
42 345783 13718 56 96448
43 501749 4525 89 93811
44 379983 6869 88 117520
45 387475 4628 39 69159
46 377305 3689 25 101792
47 370837 4891 49 210568
48 430866 7489 149 136996
49 469107 4901 58 121920
50 194493 2284 41 76403
51 530670 3160 90 108094
52 518365 4150 136 134759
53 491303 7285 97 188873
54 527021 1134 63 146216
55 233773 4658 114 156608
56 405972 2384 77 61348
57 652925 3748 6 50350
58 446211 5371 47 87720
59 341340 1285 51 99489
60 387699 9327 85 87419
61 493408 5565 43 94355
62 146494 1528 32 60326
63 414462 3122 25 94670
64 364304 7561 77 82425
65 355178 2675 54 59017
66 357760 13253 251 90829
67 261216 880 15 80791
68 397144 2053 44 100423
69 374943 1424 73 131116
70 424898 4036 85 100269
71 202055 3045 49 27330
72 378525 5119 38 39039
73 310768 1431 35 106885
74 325738 554 9 79285
75 394510 1975 34 118881
76 247060 1765 20 77623
77 368078 1012 29 114768
78 236761 810 11 74015
79 312378 1280 52 69465
80 339836 666 13 117869
81 347385 1380 29 60982
82 426280 4677 66 90131
83 352850 876 33 138971
84 301881 814 15 39625
85 377516 514 15 102725
86 357312 5692 68 64239
87 458343 3642 100 90262
88 354228 540 13 103960
89 308636 2099 45 106611
90 386212 567 14 103345
91 393343 2001 36 95551
92 378509 2949 40 82903
93 452469 2253 68 63593
94 364839 6533 29 126910
95 358649 1889 43 37527
96 376641 3055 30 60247
97 429112 272 9 112995
98 330546 1414 22 70184
99 403560 2564 19 130140
100 317892 1383 9 73221
101 307528 1261 31 76114
102 235133 975 19 90534
103 299243 3366 55 108479
104 314073 576 8 113761
105 368186 1686 28 68696
106 269661 746 29 71561
107 125390 3192 48 59831
108 510834 2045 16 97890
109 321896 5702 47 101481
110 249898 1932 20 72954
111 408881 936 22 67939
112 158492 3437 33 48022
113 292154 5131 44 86111
114 289513 2397 13 74020
115 378049 1389 6 57530
116 343466 1503 35 56364
117 332743 402 8 84990
118 442882 2239 17 88590
119 214215 2234 11 77200
120 315688 837 21 61262
121 375195 10579 92 110309
122 334280 875 12 67000
123 355864 1585 112 93099
124 480382 1659 25 107577
125 353058 2647 17 62920
126 217193 3294 23 75832
127 315380 0 0 60720
128 314533 94 10 60793
129 318056 422 23 57935
130 315380 0 0 60720
131 314353 34 7 60630
132 369448 1558 25 55637
133 315380 0 1 60720
134 312846 43 20 60887
135 312075 645 4 60720
136 315009 316 4 60505
137 318903 115 10 60945
138 314887 5 1 60720
139 314913 897 4 60720
140 315380 0 0 60720
141 325506 389 8 58990
142 315380 0 0 60720
143 298568 1002 11 56750
144 315834 36 4 60894
145 329784 460 15 63346
146 312878 309 9 56535
147 315380 0 0 60720
148 314987 9 7 60835
149 325249 271 2 60720
150 315877 14 0 61016
151 291650 520 7 58650
152 305959 1766 46 60438
153 315380 0 5 60720
154 297765 458 7 58625
155 315245 20 2 60938
156 315380 0 0 60720
157 315380 0 0 60720
158 315236 98 2 61490
159 336425 405 5 60845
160 315380 0 0 60720
161 315380 0 0 60720
162 315380 0 0 60720
163 315380 0 0 60720
164 306268 483 7 60830
165 302187 454 24 63261
166 314882 47 1 60720
167 315380 0 0 60720
168 382712 757 18 45689
169 341570 4655 55 60720
170 315380 0 0 60720
171 315380 0 0 60720
172 312412 36 3 61564
173 315380 0 0 60720
174 309596 203 9 61938
175 315380 0 0 60720
176 315547 126 8 60951
177 313267 400 113 60720
178 316176 71 0 60745
179 315380 0 0 60720
180 315380 0 0 60720
181 359335 972 19 71642
182 330068 531 11 71641
183 314289 2461 25 55792
184 297413 378 16 71873
185 314806 23 5 62555
186 333210 638 11 60370
187 352108 2300 23 64873
188 313332 149 6 62041
189 291787 226 5 65745
190 315380 0 0 60720
191 318745 275 7 59500
192 315380 0 0 60720
193 315366 141 7 61630
194 315380 0 0 60720
195 315688 28 3 60890
196 315380 0 0 60720
197 409642 4980 89 113521
198 315380 0 0 60720
199 315380 0 0 60720
200 269587 472 19 80045
201 315380 0 0 60720
202 315380 0 0 60720
203 315380 0 0 60720
204 300962 203 12 50804
205 325479 496 12 87390
206 316155 10 5 61656
207 318574 63 2 65688
208 315380 0 0 60720
209 343613 1136 26 48522
210 306948 265 3 60720
211 315380 0 0 60720
212 315380 0 0 60720
213 330059 267 11 57640
214 288985 474 10 61977
215 304485 534 5 62620
216 315380 0 2 60720
217 315688 15 6 60831
218 317736 397 7 60646
219 315380 0 2 60720
220 322331 1866 28 56225
221 296656 288 3 60510
222 315380 0 0 60720
223 315354 3 1 60698
224 312161 468 20 60720
225 315576 20 1 60805
226 314922 278 22 61404
227 314551 61 9 60720
228 315380 0 0 60720
229 312339 192 2 65276
230 315380 0 0 60720
231 298700 317 7 63915
232 321376 738 9 60720
233 315380 0 0 60720
234 303230 368 13 61686
235 315380 0 0 60720
236 315487 2 0 60743
237 315380 0 0 60720
238 315793 53 6 60349
239 315380 0 0 60720
240 315380 0 0 60720
241 315380 0 0 60720
242 312887 94 3 61360
243 315380 0 0 60720
244 315637 24 7 59818
245 324385 2332 2 72680
246 315380 0 0 60720
247 315380 0 0 60720
248 308989 131 15 61808
249 315380 0 0 60720
250 315380 0 0 60720
251 296702 206 9 53110
252 315380 0 0 60720
253 307322 167 1 64245
254 304376 622 38 73007
255 253588 2328 57 82732
256 315380 0 0 60720
257 309560 365 7 54820
258 298466 364 26 47705
259 315380 0 0 60720
260 315380 0 0 60720
261 315380 0 0 60720
262 315380 0 0 60720
263 343929 226 13 72835
264 331955 307 10 58856
265 315380 0 0 60720
266 315380 0 0 60720
267 315380 0 0 60720
268 381180 188 9 77655
269 315380 0 0 60720
270 331420 138 26 69817
271 315380 0 0 60720
272 315380 0 0 60720
273 315380 0 0 60720
274 310201 125 19 60798
275 315380 0 0 60720
276 320016 282 12 62452
277 320398 335 23 64175
278 315380 0 0 60720
279 291841 1324 29 67440
280 310670 176 8 68136
281 315380 0 0 60720
282 315380 0 0 60720
283 313491 249 26 56726
284 315380 0 0 60720
285 331323 333 9 70811
286 315380 0 0 60720
287 319210 601 5 60720
288 318098 30 3 62045
289 315380 0 0 60720
290 292754 249 13 54323
291 315380 0 0 60720
292 325176 165 12 62841
293 365959 453 19 81125
294 315380 0 0 60720
295 302409 53 10 59506
296 340968 382 9 59365
297 315380 0 0 60720
298 315380 0 0 60720
299 315380 0 0 60720
300 315380 0 9 60720
301 313164 30 4 60798
302 301164 290 1 58790
303 315380 0 1 60720
304 315380 0 0 60720
305 344425 366 14 61808
306 315394 2 12 60735
307 315380 0 0 60720
308 316647 209 19 64016
309 309836 384 17 54683
310 315380 0 0 60720
311 315380 0 0 60720
312 346611 365 32 87192
313 315380 0 0 60720
314 322031 49 14 64107
315 315656 3 8 60761
316 339445 133 4 65990
317 314964 32 0 59988
318 297141 368 20 61167
319 315372 1 5 60719
320 315380 0 0 60720
321 315380 0 0 60720
322 315380 0 0 60720
323 315380 0 0 60720
324 315380 0 0 60720
325 315380 0 0 60720
326 312502 22 1 60722
327 315380 0 0 60720
328 315380 0 0 60720
329 315380 0 0 60720
330 315380 0 0 60720
331 315380 0 0 60720
332 315380 0 0 60720
333 315380 0 0 60720
334 313729 96 4 60379
335 315388 1 1 60727
336 315371 314 4 60720
337 296139 844 20 60925
338 315380 0 0 60720
339 313880 26 1 60896
340 317698 125 10 59734
341 295580 304 12 62969
342 315380 0 0 60720
343 315380 0 0 60720
344 315380 0 0 60720
345 308256 621 13 60720
346 315380 0 0 60720
347 303677 119 3 59118
348 315380 0 0 60720
349 315380 0 0 60720
350 319369 1595 10 60720
351 318690 312 3 58598
352 314049 60 7 61124
353 325699 587 10 59595
354 314210 135 1 62065
355 315380 0 0 60720
356 315380 0 0 60720
357 322378 514 15 78780
358 315380 0 0 60720
359 315380 0 0 60720
360 315380 0 0 60720
361 315398 1 4 60722
362 315380 0 0 60720
363 315380 0 0 60720
364 308336 1763 28 61600
365 316386 180 9 59635
366 315380 0 0 60720
367 315380 0 0 60720
368 315380 0 0 60720
369 315380 0 0 60720
370 315553 218 7 60720
371 315380 0 0 60720
372 323361 448 7 59781
373 336639 227 7 76644
374 307424 174 3 64820
375 315380 0 0 60720
376 315380 0 0 60720
377 295370 121 11 56178
378 322340 607 7 60436
379 319864 2212 10 60720
380 315380 0 0 60720
381 315380 0 0 60720
382 317291 530 18 73433
383 280398 571 14 41477
384 315380 0 0 60720
385 317330 78 12 62700
386 238125 2489 29 67804
387 327071 131 3 59661
388 309038 923 6 58620
389 314210 72 3 60398
390 307930 572 8 58580
391 322327 397 10 62710
392 292136 450 6 59325
393 263276 622 8 60950
394 367655 694 6 68060
395 283910 3425 9 83620
396 283587 562 8 58456
397 243650 4917 26 52811
398 438493 1442 239 121173
399 296261 529 7 63870
400 230621 2126 41 21001
401 304252 1061 3 70415
402 333505 776 8 64230
403 296919 611 6 59190
404 278990 1526 21 69351
405 276898 592 7 64270
406 327007 1182 11 70694
407 317046 621 11 68005
408 304555 989 12 58930
409 298096 438 9 58320
410 231861 726 3 69980
411 309422 1303 57 69863
412 286963 7419 21 63255
413 269753 1164 15 57320
414 448243 3310 32 75230
415 165404 1920 11 79420
416 204325 965 2 73490
417 407159 3256 23 35250
418 290476 1135 20 62285
419 275311 1270 24 69206
420 246541 661 1 65920
421 253468 1013 1 69770
422 240897 2844 74 72683
423 -83265 11528 68 -14545
424 -42143 6526 20 55830
425 272713 2264 20 55174
426 215362 5109 82 67038
427 42754 3999 21 51252
428 306275 35624 244 157278
429 253537 9252 32 79510
430 372631 15236 86 77440
431 -7170 18073 69 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Orders Dividends
1.650e+05 2.201e+01 1.009e+03 1.568e+00
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-2198537 -21925 41823 55157 2923782
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.650e+05 3.063e+04 5.386 1.19e-07 ***
Costs 2.201e+01 2.155e+00 10.217 < 2e-16 ***
Orders 1.009e+03 3.497e+02 2.887 0.004092 **
Dividends 1.568e+00 4.115e-01 3.812 0.000158 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 253500 on 427 degrees of freedom
Multiple R-squared: 0.6556, Adjusted R-squared: 0.6532
F-statistic: 271 on 3 and 427 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 2.430225e-32 1.215113e-32
[2,] 1.0000000 8.737333e-53 4.368667e-53
[3,] 1.0000000 2.605460e-53 1.302730e-53
[4,] 1.0000000 3.343472e-75 1.671736e-75
[5,] 1.0000000 3.787983e-78 1.893992e-78
[6,] 1.0000000 8.070447e-84 4.035223e-84
[7,] 1.0000000 3.228478e-105 1.614239e-105
[8,] 1.0000000 2.615742e-107 1.307871e-107
[9,] 1.0000000 8.429364e-116 4.214682e-116
[10,] 1.0000000 2.058424e-123 1.029212e-123
[11,] 1.0000000 1.769716e-133 8.848580e-134
[12,] 1.0000000 1.349267e-139 6.746335e-140
[13,] 1.0000000 1.674554e-144 8.372772e-145
[14,] 1.0000000 9.980505e-153 4.990253e-153
[15,] 1.0000000 4.714535e-161 2.357268e-161
[16,] 1.0000000 2.522695e-164 1.261347e-164
[17,] 1.0000000 7.254357e-165 3.627178e-165
[18,] 1.0000000 4.523802e-171 2.261901e-171
[19,] 1.0000000 1.148182e-170 5.740911e-171
[20,] 1.0000000 7.161049e-173 3.580524e-173
[21,] 1.0000000 6.187980e-176 3.093990e-176
[22,] 1.0000000 3.089039e-175 1.544519e-175
[23,] 1.0000000 4.705676e-189 2.352838e-189
[24,] 1.0000000 4.328027e-189 2.164014e-189
[25,] 1.0000000 6.441233e-195 3.220616e-195
[26,] 1.0000000 2.904189e-196 1.452095e-196
[27,] 1.0000000 1.152662e-205 5.763311e-206
[28,] 1.0000000 4.250207e-208 2.125104e-208
[29,] 1.0000000 2.165470e-207 1.082735e-207
[30,] 1.0000000 1.983385e-208 9.916925e-209
[31,] 1.0000000 4.715109e-214 2.357554e-214
[32,] 1.0000000 6.242833e-220 3.121416e-220
[33,] 1.0000000 4.577117e-219 2.288558e-219
[34,] 1.0000000 1.143453e-219 5.717266e-220
[35,] 1.0000000 3.043080e-221 1.521540e-221
[36,] 1.0000000 1.874232e-220 9.371158e-221
[37,] 1.0000000 8.048933e-221 4.024467e-221
[38,] 1.0000000 9.401952e-220 4.700976e-220
[39,] 1.0000000 2.087604e-219 1.043802e-219
[40,] 1.0000000 1.807373e-218 9.036866e-219
[41,] 1.0000000 4.582867e-219 2.291434e-219
[42,] 1.0000000 2.172893e-218 1.086446e-218
[43,] 1.0000000 8.204378e-218 4.102189e-218
[44,] 1.0000000 4.341007e-218 2.170504e-218
[45,] 1.0000000 1.573934e-218 7.869672e-219
[46,] 1.0000000 4.482458e-218 2.241229e-218
[47,] 1.0000000 2.729742e-217 1.364871e-217
[48,] 1.0000000 8.420524e-217 4.210262e-217
[49,] 1.0000000 5.048895e-220 2.524448e-220
[50,] 1.0000000 1.097779e-219 5.488896e-220
[51,] 1.0000000 3.514781e-233 1.757391e-233
[52,] 1.0000000 1.415033e-233 7.075164e-234
[53,] 1.0000000 1.300345e-232 6.501724e-233
[54,] 1.0000000 3.363462e-232 1.681731e-232
[55,] 1.0000000 6.710222e-234 3.355111e-234
[56,] 1.0000000 1.970304e-235 9.851518e-236
[57,] 1.0000000 5.446042e-235 2.723021e-235
[58,] 1.0000000 3.000963e-234 1.500481e-234
[59,] 1.0000000 1.669619e-233 8.348097e-234
[60,] 1.0000000 8.471872e-234 4.235936e-234
[61,] 1.0000000 2.271781e-233 1.135891e-233
[62,] 1.0000000 2.140684e-232 1.070342e-232
[63,] 1.0000000 1.213507e-231 6.067534e-232
[64,] 1.0000000 6.701009e-231 3.350504e-231
[65,] 1.0000000 4.718952e-230 2.359476e-230
[66,] 1.0000000 3.813748e-231 1.906874e-231
[67,] 1.0000000 1.309411e-230 6.547056e-231
[68,] 1.0000000 1.123061e-229 5.615305e-230
[69,] 1.0000000 1.348557e-228 6.742784e-229
[70,] 1.0000000 3.161448e-228 1.580724e-228
[71,] 1.0000000 2.852140e-227 1.426070e-227
[72,] 1.0000000 3.834064e-227 1.917032e-227
[73,] 1.0000000 4.032626e-226 2.016313e-226
[74,] 1.0000000 1.524307e-225 7.621534e-226
[75,] 1.0000000 9.068393e-225 4.534197e-225
[76,] 1.0000000 1.496334e-224 7.481670e-225
[77,] 1.0000000 1.458437e-224 7.292187e-225
[78,] 1.0000000 8.173524e-224 4.086762e-224
[79,] 1.0000000 7.939794e-223 3.969897e-223
[80,] 1.0000000 2.122205e-222 1.061102e-222
[81,] 1.0000000 1.693093e-222 8.465467e-223
[82,] 1.0000000 1.266378e-221 6.331888e-222
[83,] 1.0000000 2.736905e-221 1.368453e-221
[84,] 1.0000000 2.662629e-220 1.331314e-220
[85,] 1.0000000 2.420510e-219 1.210255e-219
[86,] 1.0000000 1.576943e-218 7.884716e-219
[87,] 1.0000000 5.307071e-220 2.653536e-220
[88,] 1.0000000 5.400431e-219 2.700215e-219
[89,] 1.0000000 1.923374e-219 9.616871e-220
[90,] 1.0000000 1.707259e-219 8.536295e-220
[91,] 1.0000000 1.532807e-218 7.664036e-219
[92,] 1.0000000 1.505760e-217 7.528801e-218
[93,] 1.0000000 1.317677e-216 6.588387e-217
[94,] 1.0000000 1.295529e-215 6.477644e-216
[95,] 1.0000000 1.169988e-214 5.849940e-215
[96,] 1.0000000 9.190388e-216 4.595194e-216
[97,] 1.0000000 8.431388e-216 4.215694e-216
[98,] 1.0000000 1.403601e-216 7.018006e-217
[99,] 1.0000000 6.191384e-216 3.095692e-216
[100,] 1.0000000 2.187477e-215 1.093738e-215
[101,] 1.0000000 1.345981e-217 6.729903e-218
[102,] 1.0000000 5.648113e-219 2.824057e-219
[103,] 1.0000000 4.298096e-218 2.149048e-218
[104,] 1.0000000 8.618840e-218 4.309420e-218
[105,] 1.0000000 6.133438e-218 3.066719e-218
[106,] 1.0000000 3.810103e-218 1.905052e-218
[107,] 1.0000000 2.946051e-217 1.473026e-217
[108,] 1.0000000 2.387174e-216 1.193587e-216
[109,] 1.0000000 1.584190e-216 7.920952e-217
[110,] 1.0000000 6.364207e-216 3.182103e-216
[111,] 1.0000000 4.584798e-215 2.292399e-215
[112,] 1.0000000 3.685365e-215 1.842682e-215
[113,] 1.0000000 3.977775e-216 1.988887e-216
[114,] 1.0000000 4.087402e-215 2.043701e-215
[115,] 1.0000000 1.975149e-214 9.875747e-215
[116,] 1.0000000 1.832532e-213 9.162659e-214
[117,] 1.0000000 1.589165e-212 7.945826e-213
[118,] 1.0000000 2.197853e-212 1.098926e-212
[119,] 1.0000000 6.315757e-212 3.157878e-212
[120,] 1.0000000 2.062314e-212 1.031157e-212
[121,] 1.0000000 1.925751e-211 9.628753e-212
[122,] 1.0000000 1.901750e-210 9.508748e-211
[123,] 1.0000000 1.818254e-209 9.091270e-210
[124,] 1.0000000 1.709494e-208 8.547469e-209
[125,] 1.0000000 1.668385e-207 8.341924e-208
[126,] 1.0000000 1.434972e-207 7.174861e-208
[127,] 1.0000000 1.372971e-206 6.864856e-207
[128,] 1.0000000 1.422863e-205 7.114317e-206
[129,] 1.0000000 1.407273e-204 7.036365e-205
[130,] 1.0000000 1.363319e-203 6.816593e-204
[131,] 1.0000000 1.326767e-202 6.633836e-203
[132,] 1.0000000 1.274203e-201 6.371017e-202
[133,] 1.0000000 1.218434e-200 6.092172e-201
[134,] 1.0000000 1.162635e-199 5.813177e-200
[135,] 1.0000000 9.621111e-199 4.810556e-199
[136,] 1.0000000 9.180154e-198 4.590077e-198
[137,] 1.0000000 9.398313e-197 4.699157e-197
[138,] 1.0000000 9.064287e-196 4.532144e-196
[139,] 1.0000000 8.295840e-195 4.147920e-195
[140,] 1.0000000 7.606577e-194 3.803289e-194
[141,] 1.0000000 7.212280e-193 3.606140e-193
[142,] 1.0000000 6.995395e-192 3.497697e-192
[143,] 1.0000000 6.022484e-191 3.011242e-191
[144,] 1.0000000 5.692807e-190 2.846403e-190
[145,] 1.0000000 5.542452e-189 2.771226e-189
[146,] 1.0000000 5.721135e-188 2.860568e-188
[147,] 1.0000000 5.433775e-187 2.716888e-187
[148,] 1.0000000 5.331204e-186 2.665602e-186
[149,] 1.0000000 5.012125e-185 2.506062e-185
[150,] 1.0000000 4.663512e-184 2.331756e-184
[151,] 1.0000000 4.330508e-183 2.165254e-183
[152,] 1.0000000 4.066533e-182 2.033267e-182
[153,] 1.0000000 2.844088e-181 1.422044e-181
[154,] 1.0000000 2.625283e-180 1.312641e-180
[155,] 1.0000000 2.417622e-179 1.208811e-179
[156,] 1.0000000 2.220982e-178 1.110491e-178
[157,] 1.0000000 2.035229e-177 1.017615e-177
[158,] 1.0000000 1.943694e-176 9.718470e-177
[159,] 1.0000000 1.813072e-175 9.065358e-176
[160,] 1.0000000 1.653068e-174 8.265341e-175
[161,] 1.0000000 1.496385e-173 7.481924e-174
[162,] 1.0000000 1.854320e-174 9.271601e-175
[163,] 1.0000000 3.310271e-174 1.655135e-174
[164,] 1.0000000 3.032966e-173 1.516483e-173
[165,] 1.0000000 2.770809e-172 1.385405e-172
[166,] 1.0000000 2.578881e-171 1.289440e-171
[167,] 1.0000000 2.341563e-170 1.170782e-170
[168,] 1.0000000 2.197854e-169 1.098927e-169
[169,] 1.0000000 1.982565e-168 9.912826e-169
[170,] 1.0000000 1.808351e-167 9.041757e-168
[171,] 1.0000000 1.698804e-166 8.494022e-167
[172,] 1.0000000 1.507778e-165 7.538891e-166
[173,] 1.0000000 1.340523e-164 6.702614e-165
[174,] 1.0000000 1.187954e-163 5.939772e-164
[175,] 1.0000000 7.923520e-163 3.961760e-163
[176,] 1.0000000 7.261153e-162 3.630577e-162
[177,] 1.0000000 4.232342e-161 2.116171e-161
[178,] 1.0000000 2.572383e-160 1.286192e-160
[179,] 1.0000000 2.296461e-159 1.148231e-159
[180,] 1.0000000 1.533100e-158 7.665501e-159
[181,] 1.0000000 5.734398e-158 2.867199e-158
[182,] 1.0000000 5.100126e-157 2.550063e-157
[183,] 1.0000000 3.605273e-156 1.802636e-156
[184,] 1.0000000 3.105905e-155 1.552952e-155
[185,] 1.0000000 2.543144e-154 1.271572e-154
[186,] 1.0000000 2.175442e-153 1.087721e-153
[187,] 1.0000000 1.885944e-152 9.429722e-153
[188,] 1.0000000 1.601479e-151 8.007393e-152
[189,] 1.0000000 1.360697e-150 6.803483e-151
[190,] 1.0000000 1.147050e-149 5.735252e-150
[191,] 1.0000000 8.798543e-149 4.399272e-149
[192,] 1.0000000 7.360103e-148 3.680051e-148
[193,] 1.0000000 6.134258e-147 3.067129e-147
[194,] 1.0000000 6.352474e-147 3.176237e-147
[195,] 1.0000000 5.307375e-146 2.653687e-146
[196,] 1.0000000 4.417763e-145 2.208881e-145
[197,] 1.0000000 3.663540e-144 1.831770e-144
[198,] 1.0000000 2.811964e-143 1.405982e-143
[199,] 1.0000000 1.236683e-142 6.183415e-143
[200,] 1.0000000 1.028928e-141 5.144639e-142
[201,] 1.0000000 8.559826e-141 4.279913e-141
[202,] 1.0000000 6.988784e-140 3.494392e-140
[203,] 1.0000000 1.006209e-139 5.031044e-140
[204,] 1.0000000 8.412725e-139 4.206363e-139
[205,] 1.0000000 6.838441e-138 3.419221e-138
[206,] 1.0000000 5.537162e-137 2.768581e-137
[207,] 1.0000000 3.356023e-136 1.678011e-136
[208,] 1.0000000 2.523185e-135 1.261592e-135
[209,] 1.0000000 2.072258e-134 1.036129e-134
[210,] 1.0000000 1.656891e-133 8.284456e-134
[211,] 1.0000000 1.323388e-132 6.616938e-133
[212,] 1.0000000 1.022267e-131 5.111336e-132
[213,] 1.0000000 8.074498e-131 4.037249e-131
[214,] 1.0000000 3.870213e-130 1.935107e-130
[215,] 1.0000000 3.048905e-129 1.524452e-129
[216,] 1.0000000 2.379499e-128 1.189749e-128
[217,] 1.0000000 1.850683e-127 9.253415e-128
[218,] 1.0000000 1.452543e-126 7.262714e-127
[219,] 1.0000000 1.120077e-125 5.600387e-126
[220,] 1.0000000 8.730633e-125 4.365317e-125
[221,] 1.0000000 6.729103e-124 3.364551e-124
[222,] 1.0000000 5.121344e-123 2.560672e-123
[223,] 1.0000000 3.921146e-122 1.960573e-122
[224,] 1.0000000 2.960148e-121 1.480074e-121
[225,] 1.0000000 2.148883e-120 1.074442e-120
[226,] 1.0000000 1.482857e-119 7.414286e-120
[227,] 1.0000000 1.106131e-118 5.530656e-119
[228,] 1.0000000 8.426124e-118 4.213062e-118
[229,] 1.0000000 6.232351e-117 3.116175e-117
[230,] 1.0000000 4.589173e-116 2.294586e-116
[231,] 1.0000000 3.365783e-115 1.682891e-115
[232,] 1.0000000 2.447622e-114 1.223811e-114
[233,] 1.0000000 1.779829e-113 8.899144e-114
[234,] 1.0000000 1.288705e-112 6.443523e-113
[235,] 1.0000000 9.291018e-112 4.645509e-112
[236,] 1.0000000 6.764258e-111 3.382129e-111
[237,] 1.0000000 4.834640e-110 2.417320e-110
[238,] 1.0000000 3.407681e-109 1.703841e-109
[239,] 1.0000000 2.473491e-108 1.236746e-108
[240,] 1.0000000 1.744328e-107 8.721641e-108
[241,] 1.0000000 1.224741e-106 6.123706e-107
[242,] 1.0000000 8.784312e-106 4.392156e-106
[243,] 1.0000000 6.113228e-105 3.056614e-105
[244,] 1.0000000 4.235583e-104 2.117791e-104
[245,] 1.0000000 2.901977e-103 1.450988e-103
[246,] 1.0000000 1.992722e-102 9.963610e-103
[247,] 1.0000000 1.368284e-101 6.841422e-102
[248,] 1.0000000 7.701955e-101 3.850978e-101
[249,] 1.0000000 6.137527e-101 3.068764e-101
[250,] 1.0000000 4.201373e-100 2.100687e-100
[251,] 1.0000000 2.627559e-99 1.313780e-99
[252,] 1.0000000 1.495040e-98 7.475198e-99
[253,] 1.0000000 1.010580e-97 5.052898e-98
[254,] 1.0000000 6.800028e-97 3.400014e-97
[255,] 1.0000000 4.554762e-96 2.277381e-96
[256,] 1.0000000 3.036884e-95 1.518442e-95
[257,] 1.0000000 2.041187e-94 1.020594e-94
[258,] 1.0000000 1.071082e-93 5.355410e-94
[259,] 1.0000000 7.049114e-93 3.524557e-93
[260,] 1.0000000 4.617682e-92 2.308841e-92
[261,] 1.0000000 3.010811e-91 1.505406e-91
[262,] 1.0000000 1.605316e-90 8.026580e-91
[263,] 1.0000000 1.037604e-89 5.188020e-90
[264,] 1.0000000 6.849416e-89 3.424708e-89
[265,] 1.0000000 4.385247e-88 2.192624e-88
[266,] 1.0000000 2.794251e-87 1.397125e-87
[267,] 1.0000000 1.771981e-86 8.859904e-87
[268,] 1.0000000 1.143321e-85 5.716604e-86
[269,] 1.0000000 7.180531e-85 3.590266e-85
[270,] 1.0000000 4.472021e-84 2.236011e-84
[271,] 1.0000000 2.825713e-83 1.412856e-83
[272,] 1.0000000 1.748768e-82 8.743841e-83
[273,] 1.0000000 9.991851e-82 4.995925e-82
[274,] 1.0000000 5.989676e-81 2.994838e-81
[275,] 1.0000000 3.659053e-80 1.829526e-80
[276,] 1.0000000 2.224255e-79 1.112127e-79
[277,] 1.0000000 1.274548e-78 6.372738e-79
[278,] 1.0000000 7.672426e-78 3.836213e-78
[279,] 1.0000000 4.695320e-77 2.347660e-77
[280,] 1.0000000 2.798374e-76 1.399187e-76
[281,] 1.0000000 1.583245e-75 7.916224e-76
[282,] 1.0000000 9.375018e-75 4.687509e-75
[283,] 1.0000000 5.503239e-74 2.751619e-74
[284,] 1.0000000 3.271996e-73 1.635998e-73
[285,] 1.0000000 1.900963e-72 9.504813e-73
[286,] 1.0000000 1.075367e-71 5.376835e-72
[287,] 1.0000000 6.232573e-71 3.116286e-71
[288,] 1.0000000 3.564920e-70 1.782460e-70
[289,] 1.0000000 2.076835e-69 1.038418e-69
[290,] 1.0000000 8.469871e-69 4.234935e-69
[291,] 1.0000000 4.776692e-68 2.388346e-68
[292,] 1.0000000 2.679562e-67 1.339781e-67
[293,] 1.0000000 1.495114e-66 7.475570e-67
[294,] 1.0000000 8.328321e-66 4.164161e-66
[295,] 1.0000000 4.641899e-65 2.320950e-65
[296,] 1.0000000 2.602864e-64 1.301432e-64
[297,] 1.0000000 1.421785e-63 7.108925e-64
[298,] 1.0000000 7.720436e-63 3.860218e-63
[299,] 1.0000000 3.201798e-62 1.600899e-62
[300,] 1.0000000 1.729495e-61 8.647474e-62
[301,] 1.0000000 9.248263e-61 4.624132e-61
[302,] 1.0000000 5.033930e-60 2.516965e-60
[303,] 1.0000000 2.447381e-59 1.223691e-59
[304,] 1.0000000 1.287200e-58 6.435999e-59
[305,] 1.0000000 6.731408e-58 3.365704e-58
[306,] 1.0000000 3.172027e-57 1.586013e-57
[307,] 1.0000000 1.643278e-56 8.216392e-57
[308,] 1.0000000 8.566117e-56 4.283058e-56
[309,] 1.0000000 4.395012e-55 2.197506e-55
[310,] 1.0000000 2.102942e-54 1.051471e-54
[311,] 1.0000000 1.055943e-53 5.279714e-54
[312,] 1.0000000 5.451191e-53 2.725596e-53
[313,] 1.0000000 2.731798e-52 1.365899e-52
[314,] 1.0000000 1.358346e-51 6.791732e-52
[315,] 1.0000000 6.713189e-51 3.356595e-51
[316,] 1.0000000 3.297504e-50 1.648752e-50
[317,] 1.0000000 1.609760e-49 8.048801e-50
[318,] 1.0000000 7.809751e-49 3.904875e-49
[319,] 1.0000000 3.765249e-48 1.882624e-48
[320,] 1.0000000 1.822590e-47 9.112948e-48
[321,] 1.0000000 8.676706e-47 4.338353e-47
[322,] 1.0000000 4.104312e-46 2.052156e-46
[323,] 1.0000000 1.928963e-45 9.644817e-46
[324,] 1.0000000 9.007075e-45 4.503538e-45
[325,] 1.0000000 4.178279e-44 2.089140e-44
[326,] 1.0000000 1.925494e-43 9.627469e-44
[327,] 1.0000000 8.814431e-43 4.407215e-43
[328,] 1.0000000 4.014123e-42 2.007061e-42
[329,] 1.0000000 1.813715e-41 9.068577e-42
[330,] 1.0000000 8.059688e-41 4.029844e-41
[331,] 1.0000000 3.711068e-40 1.855534e-40
[332,] 1.0000000 1.642302e-39 8.211508e-40
[333,] 1.0000000 7.269284e-39 3.634642e-39
[334,] 1.0000000 3.095741e-38 1.547870e-38
[335,] 1.0000000 1.374910e-37 6.874549e-38
[336,] 1.0000000 5.920336e-37 2.960168e-37
[337,] 1.0000000 2.530858e-36 1.265429e-36
[338,] 1.0000000 1.074009e-35 5.370047e-36
[339,] 1.0000000 4.578785e-35 2.289392e-35
[340,] 1.0000000 1.914534e-34 9.572672e-35
[341,] 1.0000000 8.122871e-34 4.061436e-34
[342,] 1.0000000 3.345712e-33 1.672856e-33
[343,] 1.0000000 1.367492e-32 6.837462e-33
[344,] 1.0000000 5.012340e-32 2.506170e-32
[345,] 1.0000000 1.897526e-31 9.487632e-32
[346,] 1.0000000 7.659029e-31 3.829515e-31
[347,] 1.0000000 2.728585e-30 1.364293e-30
[348,] 1.0000000 1.088799e-29 5.443995e-30
[349,] 1.0000000 4.246838e-29 2.123419e-29
[350,] 1.0000000 1.642728e-28 8.213639e-29
[351,] 1.0000000 6.318902e-28 3.159451e-28
[352,] 1.0000000 2.409522e-27 1.204761e-27
[353,] 1.0000000 9.109048e-27 4.554524e-27
[354,] 1.0000000 3.413653e-26 1.706826e-26
[355,] 1.0000000 1.270354e-25 6.351768e-26
[356,] 1.0000000 4.676789e-25 2.338395e-25
[357,] 1.0000000 1.706157e-24 8.530787e-25
[358,] 1.0000000 6.082174e-24 3.041087e-24
[359,] 1.0000000 2.135228e-23 1.067614e-23
[360,] 1.0000000 7.578590e-23 3.789295e-23
[361,] 1.0000000 2.664120e-22 1.332060e-22
[362,] 1.0000000 9.274171e-22 4.637086e-22
[363,] 1.0000000 3.196572e-21 1.598286e-21
[364,] 1.0000000 1.085744e-20 5.428722e-21
[365,] 1.0000000 3.667283e-20 1.833642e-20
[366,] 1.0000000 1.139737e-19 5.698684e-20
[367,] 1.0000000 3.983713e-19 1.991857e-19
[368,] 1.0000000 1.374268e-18 6.871339e-19
[369,] 1.0000000 4.456894e-18 2.228447e-18
[370,] 1.0000000 1.429328e-17 7.146640e-18
[371,] 1.0000000 4.676969e-17 2.338484e-17
[372,] 1.0000000 1.379875e-16 6.899375e-17
[373,] 1.0000000 3.754796e-16 1.877398e-16
[374,] 1.0000000 1.150579e-15 5.752895e-16
[375,] 1.0000000 3.482185e-15 1.741092e-15
[376,] 1.0000000 1.120536e-14 5.602678e-15
[377,] 1.0000000 2.962742e-14 1.481371e-14
[378,] 1.0000000 8.659144e-14 4.329572e-14
[379,] 1.0000000 2.549300e-13 1.274650e-13
[380,] 1.0000000 6.813281e-13 3.406640e-13
[381,] 1.0000000 1.781602e-12 8.908010e-13
[382,] 1.0000000 4.844327e-12 2.422163e-12
[383,] 1.0000000 1.336749e-11 6.683747e-12
[384,] 1.0000000 3.605059e-11 1.802530e-11
[385,] 1.0000000 9.509091e-11 4.754545e-11
[386,] 1.0000000 2.668014e-10 1.334007e-10
[387,] 1.0000000 7.603558e-10 3.801779e-10
[388,] 1.0000000 1.526761e-09 7.633805e-10
[389,] 1.0000000 4.142377e-09 2.071188e-09
[390,] 1.0000000 1.119048e-08 5.595240e-09
[391,] 1.0000000 2.889729e-08 1.444865e-08
[392,] 1.0000000 5.852557e-08 2.926279e-08
[393,] 0.9999999 1.525159e-07 7.625794e-08
[394,] 0.9999998 3.216941e-07 1.608470e-07
[395,] 0.9999996 8.187020e-07 4.093510e-07
[396,] 0.9999991 1.759423e-06 8.797116e-07
[397,] 0.9999980 4.097108e-06 2.048554e-06
[398,] 0.9999949 1.011585e-05 5.057927e-06
[399,] 0.9999878 2.439913e-05 1.219956e-05
[400,] 0.9999732 5.355787e-05 2.677894e-05
[401,] 0.9999412 1.176265e-04 5.881327e-05
[402,] 0.9998805 2.390873e-04 1.195437e-04
[403,] 0.9997562 4.875105e-04 2.437553e-04
[404,] 0.9994845 1.030921e-03 5.154604e-04
[405,] 0.9989024 2.195244e-03 1.097622e-03
[406,] 0.9982032 3.593599e-03 1.796799e-03
[407,] 0.9964341 7.131734e-03 3.565867e-03
[408,] 0.9968652 6.269549e-03 3.134775e-03
[409,] 0.9949695 1.006091e-02 5.030456e-03
[410,] 0.9905243 1.895148e-02 9.475738e-03
[411,] 0.9977197 4.560506e-03 2.280253e-03
[412,] 0.9952420 9.516004e-03 4.758002e-03
[413,] 0.9890590 2.188210e-02 1.094105e-02
[414,] 0.9754770 4.904606e-02 2.452303e-02
[415,] 0.9482403 1.035194e-01 5.175969e-02
[416,] 0.8942298 2.115405e-01 1.057702e-01
[417,] 0.8019672 3.960655e-01 1.980328e-01
[418,] 0.8245660 3.508681e-01 1.754340e-01
> postscript(file="/var/www/html/rcomp/tmp/1ulcq1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/2mvcb1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/3mvcb1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4mvcb1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/5f4bw1292936583.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
1389642.9555 2923781.5148 1363354.5309 -2198536.5879 -394043.7941
6 7 8 9 10
109880.6241 -135942.2407 128929.0956 -336335.1695 715090.7022
11 12 13 14 15
-150648.2836 342304.4900 449640.5419 -1139733.8466 -214270.3673
16 17 18 19 20
460576.1545 -237908.5126 -238086.3678 -295879.7711 -585121.5240
21 22 23 24 25
203466.5526 -137406.0789 -119601.0233 -638071.4974 -318325.0883
26 27 28 29 30
-84290.3138 -471904.3808 -340382.1733 348304.7876 -261063.8100
31 32 33 34 35
286671.4176 -87934.3094 255923.1787 130741.0876 -54935.4616
36 37 38 39 40
-569154.1520 34312.4748 157831.6796 -51019.2550 -304575.1497
41 42 43 44 45
-459635.3342 -329005.7740 168.5371 -209377.2166 -27236.7467
46 47 48 49 50
-53785.8121 -281550.6842 -264265.0929 -53546.3755 -181996.3972
51 52 53 54 55
35728.4327 -86627.1882 -228213.2479 44143.8787 -394464.4083
56 57 58 59 60
14554.8677 320398.0866 -22045.0506 -59459.7426 -205533.9695
61 62 63 64 65
14512.1235 -179051.4887 7024.0939 -174142.4494 -15772.7669
66 67 68 69 70
-494815.7464 -65001.8342 -14962.1865 -100729.1806 -72008.4810
71 72 73 74 75
-122294.9668 1254.8171 -88695.3911 15115.6276 -34735.1194
76 77 78 79 80
-98719.3209 -28468.6427 -73250.3450 -42230.1873 -37806.4900
81 82 83 84 85
27097.0836 -49658.1948 -82701.2953 41682.6321 24953.4226
86 87 88 89 90
-102380.0855 -29335.5087 1174.8496 -115197.8123 32519.6180
91 92 93 94 95
-1901.5977 -21805.7745 69499.4547 -172296.0652 49811.4065
96 97 98 99 100
19621.4688 71825.6583 2142.7165 -41169.6913 -1469.7132
101 102 103 104 105
-35892.6448 -112495.2638 -165507.7762 -50097.8850 30071.9758
106 107 108 109 110
-53262.0014 -252162.6722 131140.5689 -175230.3204 -92234.7579
111 112 113 114 115
94522.0028 -190791.2910 -165266.5841 -57462.7139 86193.8813
116 117 118 119 120
21656.7901 17528.3553 72494.7632 -132141.0857 14990.3877
121 122 123 124 125
-288568.2348 32830.8847 -103092.5235 84908.0487 13950.5787
126 127 128 129 130
-162466.3229 55156.8048 42031.7011 29693.9470 55156.8048
131 132 133 134 135
46456.5140 57662.4310 54147.3844 31225.8242 33614.4964
136 137 138 139 140
44128.6218 45700.9845 53544.3097 30904.7354 55156.8048
141 142 143 144 145
51357.0400 55156.8048 11408.9397 50507.6768 40173.9040
146 147 148 149 150
43331.3385 55156.8048 47319.3565 57040.9193 54881.3370
151 152 153 154 155
16159.7755 -39133.5926 50109.7026 23678.9119 52220.7451
156 157 158 159 160
55156.8048 55156.8048 49628.8010 62042.6027 55156.8048
161 162 163 164 165
55156.8048 55156.8048 55156.8048 28173.1243 3757.5271
166 167 168 169 170
52614.6829 55156.8048 111229.1895 -76650.7953 55156.8048
171 172 173 174 175
55156.8048 47044.2411 55156.8048 33908.6287 55156.8048
176 177 178 179 180
44112.2506 -69826.6760 54350.5341 55156.8048 55156.8048
181 182 183 184 185
41403.7863 29922.3002 -17619.1547 -4775.3971 46151.2683
186 187 188 189 190
48386.6125 11520.0678 41700.1465 13659.9087 55156.8048
191 192 193 194 195
47315.2565 55156.8048 43545.4751 55156.8048 51553.4904
196 197 198 199 200
55156.8048 -132869.2446 55156.8048 55156.8048 -50516.3690
201 202 203 204 205
55156.8048 55156.8048 55156.8048 39709.4000 393.0099
206 207 208 209 210
49196.4916 47153.0040 55156.8048 51267.7733 37862.5884
211 212 213 214 215
55156.8048 55156.8048 57684.9961 6260.9944 24478.6930
216 217 218 219 220
53137.9639 48903.9612 41823.0012 53137.9639 -185.6711
221 222 223 224 225
27393.6179 55156.8048 54089.8453 21446.4109 53769.7683
226 227 228 229 230
25298.5915 43900.1103 55156.8048 38724.2768 55156.8048
231 232 233 234 235
19420.9584 35821.0062 55156.8048 20267.7317 55156.8048
236 237 238 239 240
55183.7008 55156.8048 48928.3830 55156.8048 55156.8048
241 242 243 244 245
55156.8048 46562.3377 55156.8048 49234.2381 -7954.4067
246 247 248 249 250
55156.8048 55156.8048 29034.0780 55156.8048 55156.8048
251 252 253 254 255
34794.7895 55156.8048 36884.1340 -27169.9018 -149947.4502
256 257 258 259 260
55156.8048 43489.2216 24397.7124 55156.8048 55156.8048
261 262 263 264 265
55156.8048 55156.8048 46606.2917 57802.5942 55156.8048
266 267 268 269 270
55156.8048 55156.8048 81171.6505 55156.8048 27645.6970
271 272 273 274 275
55156.8048 55156.8048 55156.8048 27924.6122 55156.8048
276 277 278 279 280
38755.0091 24164.1677 55156.8048 -37343.0787 26865.2489
281 282 283 284 285
55156.8048 55156.8048 27805.5139 55156.8048 38856.9030
286 287 288 289 290
55156.8048 40708.7327 52107.9101 55156.8048 23959.9458
291 292 293 294 295
55156.8048 45880.6314 44579.9973 55156.8048 32828.8978
296 297 298 299 300
65375.5583 55156.8048 55156.8048 55156.8048 46072.0207
301 302 303 304 305
48120.3368 36574.1491 54147.3844 55156.8048 60305.9911
306 307 308 309 310
42990.2029 55156.8048 27474.1123 33467.6100 55156.8048
311 312 313 314 315
55156.8048 4531.1075 55156.8048 41284.8743 47227.0903
316 317 318 319 320
63990.4489 55184.4267 7926.8099 50081.2561 55156.8048
321 322 323 324 325
55156.8048 55156.8048 55156.8048 55156.8048 55156.8048
326 327 328 329 330
50781.9191 55156.8048 55156.8048 55156.8048 55156.8048
331 332 333 334 335
55156.8048 55156.8048 55156.8048 47889.5289 54122.3903
336 337 338 339 340
44197.4366 -3174.7314 55156.8048 51798.9505 46175.2185
341 342 343 344 345
13023.7963 55156.8048 55156.8048 55156.8048 21239.0706
346 347 348 349 350
55156.8048 40318.4114 55156.8048 55156.8048 13937.7951
351 352 353 354 355
51898.1207 44805.3156 44223.3366 47895.8151 55156.8048
356 357 358 359 360
55156.8048 7371.7650 55156.8048 55156.8048 55156.8048
361 362 363 364 365
51111.9712 55156.8048 55156.8048 -20343.5093 44817.0938
366 367 368 369 370
55156.8048 55156.8048 55156.8048 55156.8048 43464.6080
371 372 373 374 375
55156.8048 47681.9423 39376.6041 33911.3345 55156.8048
376 377 378 379 380
55156.8048 28503.2373 42133.2398 849.5865 55156.8048
381 382 383 384 385
55156.8048 7290.7240 23653.9248 55156.8048 40171.0802
386 387 388 389 390
-117277.3791 62596.5683 25732.2345 49878.5072 30395.3699
391 392 393 394 395
40148.4756 18137.5424 -19076.5837 74584.5603 -96716.4176
396 397 398 399 400
6467.0060 -138660.7052 -189544.2257 12385.3742 -55495.2182
401 402 403 404 405
2436.6634 42617.6282 19587.8792 -49563.0225 -8991.9428
406 407 408 409 410
14014.8986 20621.8119 13253.5098 22909.7443 -61897.0645
411 412 413 414 415
-51363.8727 -161762.7527 -25904.1719 60090.8563 -177521.3404
416 417 418 419 420
-99190.4423 91986.7555 -17377.1556 -50407.0390 -37399.3791
421 422 423 424 425
-44260.1342 -175397.0259 -547868.0565 -458554.3229 -48841.8151
426 427 428 429 430
-250017.3396 -311854.7009 -1135952.0870 -272140.7589 -336046.0989
431
-682476.5187
> postscript(file="/var/www/html/rcomp/tmp/6f4bw1292936583.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 1389642.9555 NA
1 2923781.5148 1389642.9555
2 1363354.5309 2923781.5148
3 -2198536.5879 1363354.5309
4 -394043.7941 -2198536.5879
5 109880.6241 -394043.7941
6 -135942.2407 109880.6241
7 128929.0956 -135942.2407
8 -336335.1695 128929.0956
9 715090.7022 -336335.1695
10 -150648.2836 715090.7022
11 342304.4900 -150648.2836
12 449640.5419 342304.4900
13 -1139733.8466 449640.5419
14 -214270.3673 -1139733.8466
15 460576.1545 -214270.3673
16 -237908.5126 460576.1545
17 -238086.3678 -237908.5126
18 -295879.7711 -238086.3678
19 -585121.5240 -295879.7711
20 203466.5526 -585121.5240
21 -137406.0789 203466.5526
22 -119601.0233 -137406.0789
23 -638071.4974 -119601.0233
24 -318325.0883 -638071.4974
25 -84290.3138 -318325.0883
26 -471904.3808 -84290.3138
27 -340382.1733 -471904.3808
28 348304.7876 -340382.1733
29 -261063.8100 348304.7876
30 286671.4176 -261063.8100
31 -87934.3094 286671.4176
32 255923.1787 -87934.3094
33 130741.0876 255923.1787
34 -54935.4616 130741.0876
35 -569154.1520 -54935.4616
36 34312.4748 -569154.1520
37 157831.6796 34312.4748
38 -51019.2550 157831.6796
39 -304575.1497 -51019.2550
40 -459635.3342 -304575.1497
41 -329005.7740 -459635.3342
42 168.5371 -329005.7740
43 -209377.2166 168.5371
44 -27236.7467 -209377.2166
45 -53785.8121 -27236.7467
46 -281550.6842 -53785.8121
47 -264265.0929 -281550.6842
48 -53546.3755 -264265.0929
49 -181996.3972 -53546.3755
50 35728.4327 -181996.3972
51 -86627.1882 35728.4327
52 -228213.2479 -86627.1882
53 44143.8787 -228213.2479
54 -394464.4083 44143.8787
55 14554.8677 -394464.4083
56 320398.0866 14554.8677
57 -22045.0506 320398.0866
58 -59459.7426 -22045.0506
59 -205533.9695 -59459.7426
60 14512.1235 -205533.9695
61 -179051.4887 14512.1235
62 7024.0939 -179051.4887
63 -174142.4494 7024.0939
64 -15772.7669 -174142.4494
65 -494815.7464 -15772.7669
66 -65001.8342 -494815.7464
67 -14962.1865 -65001.8342
68 -100729.1806 -14962.1865
69 -72008.4810 -100729.1806
70 -122294.9668 -72008.4810
71 1254.8171 -122294.9668
72 -88695.3911 1254.8171
73 15115.6276 -88695.3911
74 -34735.1194 15115.6276
75 -98719.3209 -34735.1194
76 -28468.6427 -98719.3209
77 -73250.3450 -28468.6427
78 -42230.1873 -73250.3450
79 -37806.4900 -42230.1873
80 27097.0836 -37806.4900
81 -49658.1948 27097.0836
82 -82701.2953 -49658.1948
83 41682.6321 -82701.2953
84 24953.4226 41682.6321
85 -102380.0855 24953.4226
86 -29335.5087 -102380.0855
87 1174.8496 -29335.5087
88 -115197.8123 1174.8496
89 32519.6180 -115197.8123
90 -1901.5977 32519.6180
91 -21805.7745 -1901.5977
92 69499.4547 -21805.7745
93 -172296.0652 69499.4547
94 49811.4065 -172296.0652
95 19621.4688 49811.4065
96 71825.6583 19621.4688
97 2142.7165 71825.6583
98 -41169.6913 2142.7165
99 -1469.7132 -41169.6913
100 -35892.6448 -1469.7132
101 -112495.2638 -35892.6448
102 -165507.7762 -112495.2638
103 -50097.8850 -165507.7762
104 30071.9758 -50097.8850
105 -53262.0014 30071.9758
106 -252162.6722 -53262.0014
107 131140.5689 -252162.6722
108 -175230.3204 131140.5689
109 -92234.7579 -175230.3204
110 94522.0028 -92234.7579
111 -190791.2910 94522.0028
112 -165266.5841 -190791.2910
113 -57462.7139 -165266.5841
114 86193.8813 -57462.7139
115 21656.7901 86193.8813
116 17528.3553 21656.7901
117 72494.7632 17528.3553
118 -132141.0857 72494.7632
119 14990.3877 -132141.0857
120 -288568.2348 14990.3877
121 32830.8847 -288568.2348
122 -103092.5235 32830.8847
123 84908.0487 -103092.5235
124 13950.5787 84908.0487
125 -162466.3229 13950.5787
126 55156.8048 -162466.3229
127 42031.7011 55156.8048
128 29693.9470 42031.7011
129 55156.8048 29693.9470
130 46456.5140 55156.8048
131 57662.4310 46456.5140
132 54147.3844 57662.4310
133 31225.8242 54147.3844
134 33614.4964 31225.8242
135 44128.6218 33614.4964
136 45700.9845 44128.6218
137 53544.3097 45700.9845
138 30904.7354 53544.3097
139 55156.8048 30904.7354
140 51357.0400 55156.8048
141 55156.8048 51357.0400
142 11408.9397 55156.8048
143 50507.6768 11408.9397
144 40173.9040 50507.6768
145 43331.3385 40173.9040
146 55156.8048 43331.3385
147 47319.3565 55156.8048
148 57040.9193 47319.3565
149 54881.3370 57040.9193
150 16159.7755 54881.3370
151 -39133.5926 16159.7755
152 50109.7026 -39133.5926
153 23678.9119 50109.7026
154 52220.7451 23678.9119
155 55156.8048 52220.7451
156 55156.8048 55156.8048
157 49628.8010 55156.8048
158 62042.6027 49628.8010
159 55156.8048 62042.6027
160 55156.8048 55156.8048
161 55156.8048 55156.8048
162 55156.8048 55156.8048
163 28173.1243 55156.8048
164 3757.5271 28173.1243
165 52614.6829 3757.5271
166 55156.8048 52614.6829
167 111229.1895 55156.8048
168 -76650.7953 111229.1895
169 55156.8048 -76650.7953
170 55156.8048 55156.8048
171 47044.2411 55156.8048
172 55156.8048 47044.2411
173 33908.6287 55156.8048
174 55156.8048 33908.6287
175 44112.2506 55156.8048
176 -69826.6760 44112.2506
177 54350.5341 -69826.6760
178 55156.8048 54350.5341
179 55156.8048 55156.8048
180 41403.7863 55156.8048
181 29922.3002 41403.7863
182 -17619.1547 29922.3002
183 -4775.3971 -17619.1547
184 46151.2683 -4775.3971
185 48386.6125 46151.2683
186 11520.0678 48386.6125
187 41700.1465 11520.0678
188 13659.9087 41700.1465
189 55156.8048 13659.9087
190 47315.2565 55156.8048
191 55156.8048 47315.2565
192 43545.4751 55156.8048
193 55156.8048 43545.4751
194 51553.4904 55156.8048
195 55156.8048 51553.4904
196 -132869.2446 55156.8048
197 55156.8048 -132869.2446
198 55156.8048 55156.8048
199 -50516.3690 55156.8048
200 55156.8048 -50516.3690
201 55156.8048 55156.8048
202 55156.8048 55156.8048
203 39709.4000 55156.8048
204 393.0099 39709.4000
205 49196.4916 393.0099
206 47153.0040 49196.4916
207 55156.8048 47153.0040
208 51267.7733 55156.8048
209 37862.5884 51267.7733
210 55156.8048 37862.5884
211 55156.8048 55156.8048
212 57684.9961 55156.8048
213 6260.9944 57684.9961
214 24478.6930 6260.9944
215 53137.9639 24478.6930
216 48903.9612 53137.9639
217 41823.0012 48903.9612
218 53137.9639 41823.0012
219 -185.6711 53137.9639
220 27393.6179 -185.6711
221 55156.8048 27393.6179
222 54089.8453 55156.8048
223 21446.4109 54089.8453
224 53769.7683 21446.4109
225 25298.5915 53769.7683
226 43900.1103 25298.5915
227 55156.8048 43900.1103
228 38724.2768 55156.8048
229 55156.8048 38724.2768
230 19420.9584 55156.8048
231 35821.0062 19420.9584
232 55156.8048 35821.0062
233 20267.7317 55156.8048
234 55156.8048 20267.7317
235 55183.7008 55156.8048
236 55156.8048 55183.7008
237 48928.3830 55156.8048
238 55156.8048 48928.3830
239 55156.8048 55156.8048
240 55156.8048 55156.8048
241 46562.3377 55156.8048
242 55156.8048 46562.3377
243 49234.2381 55156.8048
244 -7954.4067 49234.2381
245 55156.8048 -7954.4067
246 55156.8048 55156.8048
247 29034.0780 55156.8048
248 55156.8048 29034.0780
249 55156.8048 55156.8048
250 34794.7895 55156.8048
251 55156.8048 34794.7895
252 36884.1340 55156.8048
253 -27169.9018 36884.1340
254 -149947.4502 -27169.9018
255 55156.8048 -149947.4502
256 43489.2216 55156.8048
257 24397.7124 43489.2216
258 55156.8048 24397.7124
259 55156.8048 55156.8048
260 55156.8048 55156.8048
261 55156.8048 55156.8048
262 46606.2917 55156.8048
263 57802.5942 46606.2917
264 55156.8048 57802.5942
265 55156.8048 55156.8048
266 55156.8048 55156.8048
267 81171.6505 55156.8048
268 55156.8048 81171.6505
269 27645.6970 55156.8048
270 55156.8048 27645.6970
271 55156.8048 55156.8048
272 55156.8048 55156.8048
273 27924.6122 55156.8048
274 55156.8048 27924.6122
275 38755.0091 55156.8048
276 24164.1677 38755.0091
277 55156.8048 24164.1677
278 -37343.0787 55156.8048
279 26865.2489 -37343.0787
280 55156.8048 26865.2489
281 55156.8048 55156.8048
282 27805.5139 55156.8048
283 55156.8048 27805.5139
284 38856.9030 55156.8048
285 55156.8048 38856.9030
286 40708.7327 55156.8048
287 52107.9101 40708.7327
288 55156.8048 52107.9101
289 23959.9458 55156.8048
290 55156.8048 23959.9458
291 45880.6314 55156.8048
292 44579.9973 45880.6314
293 55156.8048 44579.9973
294 32828.8978 55156.8048
295 65375.5583 32828.8978
296 55156.8048 65375.5583
297 55156.8048 55156.8048
298 55156.8048 55156.8048
299 46072.0207 55156.8048
300 48120.3368 46072.0207
301 36574.1491 48120.3368
302 54147.3844 36574.1491
303 55156.8048 54147.3844
304 60305.9911 55156.8048
305 42990.2029 60305.9911
306 55156.8048 42990.2029
307 27474.1123 55156.8048
308 33467.6100 27474.1123
309 55156.8048 33467.6100
310 55156.8048 55156.8048
311 4531.1075 55156.8048
312 55156.8048 4531.1075
313 41284.8743 55156.8048
314 47227.0903 41284.8743
315 63990.4489 47227.0903
316 55184.4267 63990.4489
317 7926.8099 55184.4267
318 50081.2561 7926.8099
319 55156.8048 50081.2561
320 55156.8048 55156.8048
321 55156.8048 55156.8048
322 55156.8048 55156.8048
323 55156.8048 55156.8048
324 55156.8048 55156.8048
325 50781.9191 55156.8048
326 55156.8048 50781.9191
327 55156.8048 55156.8048
328 55156.8048 55156.8048
329 55156.8048 55156.8048
330 55156.8048 55156.8048
331 55156.8048 55156.8048
332 55156.8048 55156.8048
333 47889.5289 55156.8048
334 54122.3903 47889.5289
335 44197.4366 54122.3903
336 -3174.7314 44197.4366
337 55156.8048 -3174.7314
338 51798.9505 55156.8048
339 46175.2185 51798.9505
340 13023.7963 46175.2185
341 55156.8048 13023.7963
342 55156.8048 55156.8048
343 55156.8048 55156.8048
344 21239.0706 55156.8048
345 55156.8048 21239.0706
346 40318.4114 55156.8048
347 55156.8048 40318.4114
348 55156.8048 55156.8048
349 13937.7951 55156.8048
350 51898.1207 13937.7951
351 44805.3156 51898.1207
352 44223.3366 44805.3156
353 47895.8151 44223.3366
354 55156.8048 47895.8151
355 55156.8048 55156.8048
356 7371.7650 55156.8048
357 55156.8048 7371.7650
358 55156.8048 55156.8048
359 55156.8048 55156.8048
360 51111.9712 55156.8048
361 55156.8048 51111.9712
362 55156.8048 55156.8048
363 -20343.5093 55156.8048
364 44817.0938 -20343.5093
365 55156.8048 44817.0938
366 55156.8048 55156.8048
367 55156.8048 55156.8048
368 55156.8048 55156.8048
369 43464.6080 55156.8048
370 55156.8048 43464.6080
371 47681.9423 55156.8048
372 39376.6041 47681.9423
373 33911.3345 39376.6041
374 55156.8048 33911.3345
375 55156.8048 55156.8048
376 28503.2373 55156.8048
377 42133.2398 28503.2373
378 849.5865 42133.2398
379 55156.8048 849.5865
380 55156.8048 55156.8048
381 7290.7240 55156.8048
382 23653.9248 7290.7240
383 55156.8048 23653.9248
384 40171.0802 55156.8048
385 -117277.3791 40171.0802
386 62596.5683 -117277.3791
387 25732.2345 62596.5683
388 49878.5072 25732.2345
389 30395.3699 49878.5072
390 40148.4756 30395.3699
391 18137.5424 40148.4756
392 -19076.5837 18137.5424
393 74584.5603 -19076.5837
394 -96716.4176 74584.5603
395 6467.0060 -96716.4176
396 -138660.7052 6467.0060
397 -189544.2257 -138660.7052
398 12385.3742 -189544.2257
399 -55495.2182 12385.3742
400 2436.6634 -55495.2182
401 42617.6282 2436.6634
402 19587.8792 42617.6282
403 -49563.0225 19587.8792
404 -8991.9428 -49563.0225
405 14014.8986 -8991.9428
406 20621.8119 14014.8986
407 13253.5098 20621.8119
408 22909.7443 13253.5098
409 -61897.0645 22909.7443
410 -51363.8727 -61897.0645
411 -161762.7527 -51363.8727
412 -25904.1719 -161762.7527
413 60090.8563 -25904.1719
414 -177521.3404 60090.8563
415 -99190.4423 -177521.3404
416 91986.7555 -99190.4423
417 -17377.1556 91986.7555
418 -50407.0390 -17377.1556
419 -37399.3791 -50407.0390
420 -44260.1342 -37399.3791
421 -175397.0259 -44260.1342
422 -547868.0565 -175397.0259
423 -458554.3229 -547868.0565
424 -48841.8151 -458554.3229
425 -250017.3396 -48841.8151
426 -311854.7009 -250017.3396
427 -1135952.0870 -311854.7009
428 -272140.7589 -1135952.0870
429 -336046.0989 -272140.7589
430 -682476.5187 -336046.0989
431 NA -682476.5187
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2923781.5148 1389642.9555
[2,] 1363354.5309 2923781.5148
[3,] -2198536.5879 1363354.5309
[4,] -394043.7941 -2198536.5879
[5,] 109880.6241 -394043.7941
[6,] -135942.2407 109880.6241
[7,] 128929.0956 -135942.2407
[8,] -336335.1695 128929.0956
[9,] 715090.7022 -336335.1695
[10,] -150648.2836 715090.7022
[11,] 342304.4900 -150648.2836
[12,] 449640.5419 342304.4900
[13,] -1139733.8466 449640.5419
[14,] -214270.3673 -1139733.8466
[15,] 460576.1545 -214270.3673
[16,] -237908.5126 460576.1545
[17,] -238086.3678 -237908.5126
[18,] -295879.7711 -238086.3678
[19,] -585121.5240 -295879.7711
[20,] 203466.5526 -585121.5240
[21,] -137406.0789 203466.5526
[22,] -119601.0233 -137406.0789
[23,] -638071.4974 -119601.0233
[24,] -318325.0883 -638071.4974
[25,] -84290.3138 -318325.0883
[26,] -471904.3808 -84290.3138
[27,] -340382.1733 -471904.3808
[28,] 348304.7876 -340382.1733
[29,] -261063.8100 348304.7876
[30,] 286671.4176 -261063.8100
[31,] -87934.3094 286671.4176
[32,] 255923.1787 -87934.3094
[33,] 130741.0876 255923.1787
[34,] -54935.4616 130741.0876
[35,] -569154.1520 -54935.4616
[36,] 34312.4748 -569154.1520
[37,] 157831.6796 34312.4748
[38,] -51019.2550 157831.6796
[39,] -304575.1497 -51019.2550
[40,] -459635.3342 -304575.1497
[41,] -329005.7740 -459635.3342
[42,] 168.5371 -329005.7740
[43,] -209377.2166 168.5371
[44,] -27236.7467 -209377.2166
[45,] -53785.8121 -27236.7467
[46,] -281550.6842 -53785.8121
[47,] -264265.0929 -281550.6842
[48,] -53546.3755 -264265.0929
[49,] -181996.3972 -53546.3755
[50,] 35728.4327 -181996.3972
[51,] -86627.1882 35728.4327
[52,] -228213.2479 -86627.1882
[53,] 44143.8787 -228213.2479
[54,] -394464.4083 44143.8787
[55,] 14554.8677 -394464.4083
[56,] 320398.0866 14554.8677
[57,] -22045.0506 320398.0866
[58,] -59459.7426 -22045.0506
[59,] -205533.9695 -59459.7426
[60,] 14512.1235 -205533.9695
[61,] -179051.4887 14512.1235
[62,] 7024.0939 -179051.4887
[63,] -174142.4494 7024.0939
[64,] -15772.7669 -174142.4494
[65,] -494815.7464 -15772.7669
[66,] -65001.8342 -494815.7464
[67,] -14962.1865 -65001.8342
[68,] -100729.1806 -14962.1865
[69,] -72008.4810 -100729.1806
[70,] -122294.9668 -72008.4810
[71,] 1254.8171 -122294.9668
[72,] -88695.3911 1254.8171
[73,] 15115.6276 -88695.3911
[74,] -34735.1194 15115.6276
[75,] -98719.3209 -34735.1194
[76,] -28468.6427 -98719.3209
[77,] -73250.3450 -28468.6427
[78,] -42230.1873 -73250.3450
[79,] -37806.4900 -42230.1873
[80,] 27097.0836 -37806.4900
[81,] -49658.1948 27097.0836
[82,] -82701.2953 -49658.1948
[83,] 41682.6321 -82701.2953
[84,] 24953.4226 41682.6321
[85,] -102380.0855 24953.4226
[86,] -29335.5087 -102380.0855
[87,] 1174.8496 -29335.5087
[88,] -115197.8123 1174.8496
[89,] 32519.6180 -115197.8123
[90,] -1901.5977 32519.6180
[91,] -21805.7745 -1901.5977
[92,] 69499.4547 -21805.7745
[93,] -172296.0652 69499.4547
[94,] 49811.4065 -172296.0652
[95,] 19621.4688 49811.4065
[96,] 71825.6583 19621.4688
[97,] 2142.7165 71825.6583
[98,] -41169.6913 2142.7165
[99,] -1469.7132 -41169.6913
[100,] -35892.6448 -1469.7132
[101,] -112495.2638 -35892.6448
[102,] -165507.7762 -112495.2638
[103,] -50097.8850 -165507.7762
[104,] 30071.9758 -50097.8850
[105,] -53262.0014 30071.9758
[106,] -252162.6722 -53262.0014
[107,] 131140.5689 -252162.6722
[108,] -175230.3204 131140.5689
[109,] -92234.7579 -175230.3204
[110,] 94522.0028 -92234.7579
[111,] -190791.2910 94522.0028
[112,] -165266.5841 -190791.2910
[113,] -57462.7139 -165266.5841
[114,] 86193.8813 -57462.7139
[115,] 21656.7901 86193.8813
[116,] 17528.3553 21656.7901
[117,] 72494.7632 17528.3553
[118,] -132141.0857 72494.7632
[119,] 14990.3877 -132141.0857
[120,] -288568.2348 14990.3877
[121,] 32830.8847 -288568.2348
[122,] -103092.5235 32830.8847
[123,] 84908.0487 -103092.5235
[124,] 13950.5787 84908.0487
[125,] -162466.3229 13950.5787
[126,] 55156.8048 -162466.3229
[127,] 42031.7011 55156.8048
[128,] 29693.9470 42031.7011
[129,] 55156.8048 29693.9470
[130,] 46456.5140 55156.8048
[131,] 57662.4310 46456.5140
[132,] 54147.3844 57662.4310
[133,] 31225.8242 54147.3844
[134,] 33614.4964 31225.8242
[135,] 44128.6218 33614.4964
[136,] 45700.9845 44128.6218
[137,] 53544.3097 45700.9845
[138,] 30904.7354 53544.3097
[139,] 55156.8048 30904.7354
[140,] 51357.0400 55156.8048
[141,] 55156.8048 51357.0400
[142,] 11408.9397 55156.8048
[143,] 50507.6768 11408.9397
[144,] 40173.9040 50507.6768
[145,] 43331.3385 40173.9040
[146,] 55156.8048 43331.3385
[147,] 47319.3565 55156.8048
[148,] 57040.9193 47319.3565
[149,] 54881.3370 57040.9193
[150,] 16159.7755 54881.3370
[151,] -39133.5926 16159.7755
[152,] 50109.7026 -39133.5926
[153,] 23678.9119 50109.7026
[154,] 52220.7451 23678.9119
[155,] 55156.8048 52220.7451
[156,] 55156.8048 55156.8048
[157,] 49628.8010 55156.8048
[158,] 62042.6027 49628.8010
[159,] 55156.8048 62042.6027
[160,] 55156.8048 55156.8048
[161,] 55156.8048 55156.8048
[162,] 55156.8048 55156.8048
[163,] 28173.1243 55156.8048
[164,] 3757.5271 28173.1243
[165,] 52614.6829 3757.5271
[166,] 55156.8048 52614.6829
[167,] 111229.1895 55156.8048
[168,] -76650.7953 111229.1895
[169,] 55156.8048 -76650.7953
[170,] 55156.8048 55156.8048
[171,] 47044.2411 55156.8048
[172,] 55156.8048 47044.2411
[173,] 33908.6287 55156.8048
[174,] 55156.8048 33908.6287
[175,] 44112.2506 55156.8048
[176,] -69826.6760 44112.2506
[177,] 54350.5341 -69826.6760
[178,] 55156.8048 54350.5341
[179,] 55156.8048 55156.8048
[180,] 41403.7863 55156.8048
[181,] 29922.3002 41403.7863
[182,] -17619.1547 29922.3002
[183,] -4775.3971 -17619.1547
[184,] 46151.2683 -4775.3971
[185,] 48386.6125 46151.2683
[186,] 11520.0678 48386.6125
[187,] 41700.1465 11520.0678
[188,] 13659.9087 41700.1465
[189,] 55156.8048 13659.9087
[190,] 47315.2565 55156.8048
[191,] 55156.8048 47315.2565
[192,] 43545.4751 55156.8048
[193,] 55156.8048 43545.4751
[194,] 51553.4904 55156.8048
[195,] 55156.8048 51553.4904
[196,] -132869.2446 55156.8048
[197,] 55156.8048 -132869.2446
[198,] 55156.8048 55156.8048
[199,] -50516.3690 55156.8048
[200,] 55156.8048 -50516.3690
[201,] 55156.8048 55156.8048
[202,] 55156.8048 55156.8048
[203,] 39709.4000 55156.8048
[204,] 393.0099 39709.4000
[205,] 49196.4916 393.0099
[206,] 47153.0040 49196.4916
[207,] 55156.8048 47153.0040
[208,] 51267.7733 55156.8048
[209,] 37862.5884 51267.7733
[210,] 55156.8048 37862.5884
[211,] 55156.8048 55156.8048
[212,] 57684.9961 55156.8048
[213,] 6260.9944 57684.9961
[214,] 24478.6930 6260.9944
[215,] 53137.9639 24478.6930
[216,] 48903.9612 53137.9639
[217,] 41823.0012 48903.9612
[218,] 53137.9639 41823.0012
[219,] -185.6711 53137.9639
[220,] 27393.6179 -185.6711
[221,] 55156.8048 27393.6179
[222,] 54089.8453 55156.8048
[223,] 21446.4109 54089.8453
[224,] 53769.7683 21446.4109
[225,] 25298.5915 53769.7683
[226,] 43900.1103 25298.5915
[227,] 55156.8048 43900.1103
[228,] 38724.2768 55156.8048
[229,] 55156.8048 38724.2768
[230,] 19420.9584 55156.8048
[231,] 35821.0062 19420.9584
[232,] 55156.8048 35821.0062
[233,] 20267.7317 55156.8048
[234,] 55156.8048 20267.7317
[235,] 55183.7008 55156.8048
[236,] 55156.8048 55183.7008
[237,] 48928.3830 55156.8048
[238,] 55156.8048 48928.3830
[239,] 55156.8048 55156.8048
[240,] 55156.8048 55156.8048
[241,] 46562.3377 55156.8048
[242,] 55156.8048 46562.3377
[243,] 49234.2381 55156.8048
[244,] -7954.4067 49234.2381
[245,] 55156.8048 -7954.4067
[246,] 55156.8048 55156.8048
[247,] 29034.0780 55156.8048
[248,] 55156.8048 29034.0780
[249,] 55156.8048 55156.8048
[250,] 34794.7895 55156.8048
[251,] 55156.8048 34794.7895
[252,] 36884.1340 55156.8048
[253,] -27169.9018 36884.1340
[254,] -149947.4502 -27169.9018
[255,] 55156.8048 -149947.4502
[256,] 43489.2216 55156.8048
[257,] 24397.7124 43489.2216
[258,] 55156.8048 24397.7124
[259,] 55156.8048 55156.8048
[260,] 55156.8048 55156.8048
[261,] 55156.8048 55156.8048
[262,] 46606.2917 55156.8048
[263,] 57802.5942 46606.2917
[264,] 55156.8048 57802.5942
[265,] 55156.8048 55156.8048
[266,] 55156.8048 55156.8048
[267,] 81171.6505 55156.8048
[268,] 55156.8048 81171.6505
[269,] 27645.6970 55156.8048
[270,] 55156.8048 27645.6970
[271,] 55156.8048 55156.8048
[272,] 55156.8048 55156.8048
[273,] 27924.6122 55156.8048
[274,] 55156.8048 27924.6122
[275,] 38755.0091 55156.8048
[276,] 24164.1677 38755.0091
[277,] 55156.8048 24164.1677
[278,] -37343.0787 55156.8048
[279,] 26865.2489 -37343.0787
[280,] 55156.8048 26865.2489
[281,] 55156.8048 55156.8048
[282,] 27805.5139 55156.8048
[283,] 55156.8048 27805.5139
[284,] 38856.9030 55156.8048
[285,] 55156.8048 38856.9030
[286,] 40708.7327 55156.8048
[287,] 52107.9101 40708.7327
[288,] 55156.8048 52107.9101
[289,] 23959.9458 55156.8048
[290,] 55156.8048 23959.9458
[291,] 45880.6314 55156.8048
[292,] 44579.9973 45880.6314
[293,] 55156.8048 44579.9973
[294,] 32828.8978 55156.8048
[295,] 65375.5583 32828.8978
[296,] 55156.8048 65375.5583
[297,] 55156.8048 55156.8048
[298,] 55156.8048 55156.8048
[299,] 46072.0207 55156.8048
[300,] 48120.3368 46072.0207
[301,] 36574.1491 48120.3368
[302,] 54147.3844 36574.1491
[303,] 55156.8048 54147.3844
[304,] 60305.9911 55156.8048
[305,] 42990.2029 60305.9911
[306,] 55156.8048 42990.2029
[307,] 27474.1123 55156.8048
[308,] 33467.6100 27474.1123
[309,] 55156.8048 33467.6100
[310,] 55156.8048 55156.8048
[311,] 4531.1075 55156.8048
[312,] 55156.8048 4531.1075
[313,] 41284.8743 55156.8048
[314,] 47227.0903 41284.8743
[315,] 63990.4489 47227.0903
[316,] 55184.4267 63990.4489
[317,] 7926.8099 55184.4267
[318,] 50081.2561 7926.8099
[319,] 55156.8048 50081.2561
[320,] 55156.8048 55156.8048
[321,] 55156.8048 55156.8048
[322,] 55156.8048 55156.8048
[323,] 55156.8048 55156.8048
[324,] 55156.8048 55156.8048
[325,] 50781.9191 55156.8048
[326,] 55156.8048 50781.9191
[327,] 55156.8048 55156.8048
[328,] 55156.8048 55156.8048
[329,] 55156.8048 55156.8048
[330,] 55156.8048 55156.8048
[331,] 55156.8048 55156.8048
[332,] 55156.8048 55156.8048
[333,] 47889.5289 55156.8048
[334,] 54122.3903 47889.5289
[335,] 44197.4366 54122.3903
[336,] -3174.7314 44197.4366
[337,] 55156.8048 -3174.7314
[338,] 51798.9505 55156.8048
[339,] 46175.2185 51798.9505
[340,] 13023.7963 46175.2185
[341,] 55156.8048 13023.7963
[342,] 55156.8048 55156.8048
[343,] 55156.8048 55156.8048
[344,] 21239.0706 55156.8048
[345,] 55156.8048 21239.0706
[346,] 40318.4114 55156.8048
[347,] 55156.8048 40318.4114
[348,] 55156.8048 55156.8048
[349,] 13937.7951 55156.8048
[350,] 51898.1207 13937.7951
[351,] 44805.3156 51898.1207
[352,] 44223.3366 44805.3156
[353,] 47895.8151 44223.3366
[354,] 55156.8048 47895.8151
[355,] 55156.8048 55156.8048
[356,] 7371.7650 55156.8048
[357,] 55156.8048 7371.7650
[358,] 55156.8048 55156.8048
[359,] 55156.8048 55156.8048
[360,] 51111.9712 55156.8048
[361,] 55156.8048 51111.9712
[362,] 55156.8048 55156.8048
[363,] -20343.5093 55156.8048
[364,] 44817.0938 -20343.5093
[365,] 55156.8048 44817.0938
[366,] 55156.8048 55156.8048
[367,] 55156.8048 55156.8048
[368,] 55156.8048 55156.8048
[369,] 43464.6080 55156.8048
[370,] 55156.8048 43464.6080
[371,] 47681.9423 55156.8048
[372,] 39376.6041 47681.9423
[373,] 33911.3345 39376.6041
[374,] 55156.8048 33911.3345
[375,] 55156.8048 55156.8048
[376,] 28503.2373 55156.8048
[377,] 42133.2398 28503.2373
[378,] 849.5865 42133.2398
[379,] 55156.8048 849.5865
[380,] 55156.8048 55156.8048
[381,] 7290.7240 55156.8048
[382,] 23653.9248 7290.7240
[383,] 55156.8048 23653.9248
[384,] 40171.0802 55156.8048
[385,] -117277.3791 40171.0802
[386,] 62596.5683 -117277.3791
[387,] 25732.2345 62596.5683
[388,] 49878.5072 25732.2345
[389,] 30395.3699 49878.5072
[390,] 40148.4756 30395.3699
[391,] 18137.5424 40148.4756
[392,] -19076.5837 18137.5424
[393,] 74584.5603 -19076.5837
[394,] -96716.4176 74584.5603
[395,] 6467.0060 -96716.4176
[396,] -138660.7052 6467.0060
[397,] -189544.2257 -138660.7052
[398,] 12385.3742 -189544.2257
[399,] -55495.2182 12385.3742
[400,] 2436.6634 -55495.2182
[401,] 42617.6282 2436.6634
[402,] 19587.8792 42617.6282
[403,] -49563.0225 19587.8792
[404,] -8991.9428 -49563.0225
[405,] 14014.8986 -8991.9428
[406,] 20621.8119 14014.8986
[407,] 13253.5098 20621.8119
[408,] 22909.7443 13253.5098
[409,] -61897.0645 22909.7443
[410,] -51363.8727 -61897.0645
[411,] -161762.7527 -51363.8727
[412,] -25904.1719 -161762.7527
[413,] 60090.8563 -25904.1719
[414,] -177521.3404 60090.8563
[415,] -99190.4423 -177521.3404
[416,] 91986.7555 -99190.4423
[417,] -17377.1556 91986.7555
[418,] -50407.0390 -17377.1556
[419,] -37399.3791 -50407.0390
[420,] -44260.1342 -37399.3791
[421,] -175397.0259 -44260.1342
[422,] -547868.0565 -175397.0259
[423,] -458554.3229 -547868.0565
[424,] -48841.8151 -458554.3229
[425,] -250017.3396 -48841.8151
[426,] -311854.7009 -250017.3396
[427,] -1135952.0870 -311854.7009
[428,] -272140.7589 -1135952.0870
[429,] -336046.0989 -272140.7589
[430,] -682476.5187 -336046.0989
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2923781.5148 1389642.9555
2 1363354.5309 2923781.5148
3 -2198536.5879 1363354.5309
4 -394043.7941 -2198536.5879
5 109880.6241 -394043.7941
6 -135942.2407 109880.6241
7 128929.0956 -135942.2407
8 -336335.1695 128929.0956
9 715090.7022 -336335.1695
10 -150648.2836 715090.7022
11 342304.4900 -150648.2836
12 449640.5419 342304.4900
13 -1139733.8466 449640.5419
14 -214270.3673 -1139733.8466
15 460576.1545 -214270.3673
16 -237908.5126 460576.1545
17 -238086.3678 -237908.5126
18 -295879.7711 -238086.3678
19 -585121.5240 -295879.7711
20 203466.5526 -585121.5240
21 -137406.0789 203466.5526
22 -119601.0233 -137406.0789
23 -638071.4974 -119601.0233
24 -318325.0883 -638071.4974
25 -84290.3138 -318325.0883
26 -471904.3808 -84290.3138
27 -340382.1733 -471904.3808
28 348304.7876 -340382.1733
29 -261063.8100 348304.7876
30 286671.4176 -261063.8100
31 -87934.3094 286671.4176
32 255923.1787 -87934.3094
33 130741.0876 255923.1787
34 -54935.4616 130741.0876
35 -569154.1520 -54935.4616
36 34312.4748 -569154.1520
37 157831.6796 34312.4748
38 -51019.2550 157831.6796
39 -304575.1497 -51019.2550
40 -459635.3342 -304575.1497
41 -329005.7740 -459635.3342
42 168.5371 -329005.7740
43 -209377.2166 168.5371
44 -27236.7467 -209377.2166
45 -53785.8121 -27236.7467
46 -281550.6842 -53785.8121
47 -264265.0929 -281550.6842
48 -53546.3755 -264265.0929
49 -181996.3972 -53546.3755
50 35728.4327 -181996.3972
51 -86627.1882 35728.4327
52 -228213.2479 -86627.1882
53 44143.8787 -228213.2479
54 -394464.4083 44143.8787
55 14554.8677 -394464.4083
56 320398.0866 14554.8677
57 -22045.0506 320398.0866
58 -59459.7426 -22045.0506
59 -205533.9695 -59459.7426
60 14512.1235 -205533.9695
61 -179051.4887 14512.1235
62 7024.0939 -179051.4887
63 -174142.4494 7024.0939
64 -15772.7669 -174142.4494
65 -494815.7464 -15772.7669
66 -65001.8342 -494815.7464
67 -14962.1865 -65001.8342
68 -100729.1806 -14962.1865
69 -72008.4810 -100729.1806
70 -122294.9668 -72008.4810
71 1254.8171 -122294.9668
72 -88695.3911 1254.8171
73 15115.6276 -88695.3911
74 -34735.1194 15115.6276
75 -98719.3209 -34735.1194
76 -28468.6427 -98719.3209
77 -73250.3450 -28468.6427
78 -42230.1873 -73250.3450
79 -37806.4900 -42230.1873
80 27097.0836 -37806.4900
81 -49658.1948 27097.0836
82 -82701.2953 -49658.1948
83 41682.6321 -82701.2953
84 24953.4226 41682.6321
85 -102380.0855 24953.4226
86 -29335.5087 -102380.0855
87 1174.8496 -29335.5087
88 -115197.8123 1174.8496
89 32519.6180 -115197.8123
90 -1901.5977 32519.6180
91 -21805.7745 -1901.5977
92 69499.4547 -21805.7745
93 -172296.0652 69499.4547
94 49811.4065 -172296.0652
95 19621.4688 49811.4065
96 71825.6583 19621.4688
97 2142.7165 71825.6583
98 -41169.6913 2142.7165
99 -1469.7132 -41169.6913
100 -35892.6448 -1469.7132
101 -112495.2638 -35892.6448
102 -165507.7762 -112495.2638
103 -50097.8850 -165507.7762
104 30071.9758 -50097.8850
105 -53262.0014 30071.9758
106 -252162.6722 -53262.0014
107 131140.5689 -252162.6722
108 -175230.3204 131140.5689
109 -92234.7579 -175230.3204
110 94522.0028 -92234.7579
111 -190791.2910 94522.0028
112 -165266.5841 -190791.2910
113 -57462.7139 -165266.5841
114 86193.8813 -57462.7139
115 21656.7901 86193.8813
116 17528.3553 21656.7901
117 72494.7632 17528.3553
118 -132141.0857 72494.7632
119 14990.3877 -132141.0857
120 -288568.2348 14990.3877
121 32830.8847 -288568.2348
122 -103092.5235 32830.8847
123 84908.0487 -103092.5235
124 13950.5787 84908.0487
125 -162466.3229 13950.5787
126 55156.8048 -162466.3229
127 42031.7011 55156.8048
128 29693.9470 42031.7011
129 55156.8048 29693.9470
130 46456.5140 55156.8048
131 57662.4310 46456.5140
132 54147.3844 57662.4310
133 31225.8242 54147.3844
134 33614.4964 31225.8242
135 44128.6218 33614.4964
136 45700.9845 44128.6218
137 53544.3097 45700.9845
138 30904.7354 53544.3097
139 55156.8048 30904.7354
140 51357.0400 55156.8048
141 55156.8048 51357.0400
142 11408.9397 55156.8048
143 50507.6768 11408.9397
144 40173.9040 50507.6768
145 43331.3385 40173.9040
146 55156.8048 43331.3385
147 47319.3565 55156.8048
148 57040.9193 47319.3565
149 54881.3370 57040.9193
150 16159.7755 54881.3370
151 -39133.5926 16159.7755
152 50109.7026 -39133.5926
153 23678.9119 50109.7026
154 52220.7451 23678.9119
155 55156.8048 52220.7451
156 55156.8048 55156.8048
157 49628.8010 55156.8048
158 62042.6027 49628.8010
159 55156.8048 62042.6027
160 55156.8048 55156.8048
161 55156.8048 55156.8048
162 55156.8048 55156.8048
163 28173.1243 55156.8048
164 3757.5271 28173.1243
165 52614.6829 3757.5271
166 55156.8048 52614.6829
167 111229.1895 55156.8048
168 -76650.7953 111229.1895
169 55156.8048 -76650.7953
170 55156.8048 55156.8048
171 47044.2411 55156.8048
172 55156.8048 47044.2411
173 33908.6287 55156.8048
174 55156.8048 33908.6287
175 44112.2506 55156.8048
176 -69826.6760 44112.2506
177 54350.5341 -69826.6760
178 55156.8048 54350.5341
179 55156.8048 55156.8048
180 41403.7863 55156.8048
181 29922.3002 41403.7863
182 -17619.1547 29922.3002
183 -4775.3971 -17619.1547
184 46151.2683 -4775.3971
185 48386.6125 46151.2683
186 11520.0678 48386.6125
187 41700.1465 11520.0678
188 13659.9087 41700.1465
189 55156.8048 13659.9087
190 47315.2565 55156.8048
191 55156.8048 47315.2565
192 43545.4751 55156.8048
193 55156.8048 43545.4751
194 51553.4904 55156.8048
195 55156.8048 51553.4904
196 -132869.2446 55156.8048
197 55156.8048 -132869.2446
198 55156.8048 55156.8048
199 -50516.3690 55156.8048
200 55156.8048 -50516.3690
201 55156.8048 55156.8048
202 55156.8048 55156.8048
203 39709.4000 55156.8048
204 393.0099 39709.4000
205 49196.4916 393.0099
206 47153.0040 49196.4916
207 55156.8048 47153.0040
208 51267.7733 55156.8048
209 37862.5884 51267.7733
210 55156.8048 37862.5884
211 55156.8048 55156.8048
212 57684.9961 55156.8048
213 6260.9944 57684.9961
214 24478.6930 6260.9944
215 53137.9639 24478.6930
216 48903.9612 53137.9639
217 41823.0012 48903.9612
218 53137.9639 41823.0012
219 -185.6711 53137.9639
220 27393.6179 -185.6711
221 55156.8048 27393.6179
222 54089.8453 55156.8048
223 21446.4109 54089.8453
224 53769.7683 21446.4109
225 25298.5915 53769.7683
226 43900.1103 25298.5915
227 55156.8048 43900.1103
228 38724.2768 55156.8048
229 55156.8048 38724.2768
230 19420.9584 55156.8048
231 35821.0062 19420.9584
232 55156.8048 35821.0062
233 20267.7317 55156.8048
234 55156.8048 20267.7317
235 55183.7008 55156.8048
236 55156.8048 55183.7008
237 48928.3830 55156.8048
238 55156.8048 48928.3830
239 55156.8048 55156.8048
240 55156.8048 55156.8048
241 46562.3377 55156.8048
242 55156.8048 46562.3377
243 49234.2381 55156.8048
244 -7954.4067 49234.2381
245 55156.8048 -7954.4067
246 55156.8048 55156.8048
247 29034.0780 55156.8048
248 55156.8048 29034.0780
249 55156.8048 55156.8048
250 34794.7895 55156.8048
251 55156.8048 34794.7895
252 36884.1340 55156.8048
253 -27169.9018 36884.1340
254 -149947.4502 -27169.9018
255 55156.8048 -149947.4502
256 43489.2216 55156.8048
257 24397.7124 43489.2216
258 55156.8048 24397.7124
259 55156.8048 55156.8048
260 55156.8048 55156.8048
261 55156.8048 55156.8048
262 46606.2917 55156.8048
263 57802.5942 46606.2917
264 55156.8048 57802.5942
265 55156.8048 55156.8048
266 55156.8048 55156.8048
267 81171.6505 55156.8048
268 55156.8048 81171.6505
269 27645.6970 55156.8048
270 55156.8048 27645.6970
271 55156.8048 55156.8048
272 55156.8048 55156.8048
273 27924.6122 55156.8048
274 55156.8048 27924.6122
275 38755.0091 55156.8048
276 24164.1677 38755.0091
277 55156.8048 24164.1677
278 -37343.0787 55156.8048
279 26865.2489 -37343.0787
280 55156.8048 26865.2489
281 55156.8048 55156.8048
282 27805.5139 55156.8048
283 55156.8048 27805.5139
284 38856.9030 55156.8048
285 55156.8048 38856.9030
286 40708.7327 55156.8048
287 52107.9101 40708.7327
288 55156.8048 52107.9101
289 23959.9458 55156.8048
290 55156.8048 23959.9458
291 45880.6314 55156.8048
292 44579.9973 45880.6314
293 55156.8048 44579.9973
294 32828.8978 55156.8048
295 65375.5583 32828.8978
296 55156.8048 65375.5583
297 55156.8048 55156.8048
298 55156.8048 55156.8048
299 46072.0207 55156.8048
300 48120.3368 46072.0207
301 36574.1491 48120.3368
302 54147.3844 36574.1491
303 55156.8048 54147.3844
304 60305.9911 55156.8048
305 42990.2029 60305.9911
306 55156.8048 42990.2029
307 27474.1123 55156.8048
308 33467.6100 27474.1123
309 55156.8048 33467.6100
310 55156.8048 55156.8048
311 4531.1075 55156.8048
312 55156.8048 4531.1075
313 41284.8743 55156.8048
314 47227.0903 41284.8743
315 63990.4489 47227.0903
316 55184.4267 63990.4489
317 7926.8099 55184.4267
318 50081.2561 7926.8099
319 55156.8048 50081.2561
320 55156.8048 55156.8048
321 55156.8048 55156.8048
322 55156.8048 55156.8048
323 55156.8048 55156.8048
324 55156.8048 55156.8048
325 50781.9191 55156.8048
326 55156.8048 50781.9191
327 55156.8048 55156.8048
328 55156.8048 55156.8048
329 55156.8048 55156.8048
330 55156.8048 55156.8048
331 55156.8048 55156.8048
332 55156.8048 55156.8048
333 47889.5289 55156.8048
334 54122.3903 47889.5289
335 44197.4366 54122.3903
336 -3174.7314 44197.4366
337 55156.8048 -3174.7314
338 51798.9505 55156.8048
339 46175.2185 51798.9505
340 13023.7963 46175.2185
341 55156.8048 13023.7963
342 55156.8048 55156.8048
343 55156.8048 55156.8048
344 21239.0706 55156.8048
345 55156.8048 21239.0706
346 40318.4114 55156.8048
347 55156.8048 40318.4114
348 55156.8048 55156.8048
349 13937.7951 55156.8048
350 51898.1207 13937.7951
351 44805.3156 51898.1207
352 44223.3366 44805.3156
353 47895.8151 44223.3366
354 55156.8048 47895.8151
355 55156.8048 55156.8048
356 7371.7650 55156.8048
357 55156.8048 7371.7650
358 55156.8048 55156.8048
359 55156.8048 55156.8048
360 51111.9712 55156.8048
361 55156.8048 51111.9712
362 55156.8048 55156.8048
363 -20343.5093 55156.8048
364 44817.0938 -20343.5093
365 55156.8048 44817.0938
366 55156.8048 55156.8048
367 55156.8048 55156.8048
368 55156.8048 55156.8048
369 43464.6080 55156.8048
370 55156.8048 43464.6080
371 47681.9423 55156.8048
372 39376.6041 47681.9423
373 33911.3345 39376.6041
374 55156.8048 33911.3345
375 55156.8048 55156.8048
376 28503.2373 55156.8048
377 42133.2398 28503.2373
378 849.5865 42133.2398
379 55156.8048 849.5865
380 55156.8048 55156.8048
381 7290.7240 55156.8048
382 23653.9248 7290.7240
383 55156.8048 23653.9248
384 40171.0802 55156.8048
385 -117277.3791 40171.0802
386 62596.5683 -117277.3791
387 25732.2345 62596.5683
388 49878.5072 25732.2345
389 30395.3699 49878.5072
390 40148.4756 30395.3699
391 18137.5424 40148.4756
392 -19076.5837 18137.5424
393 74584.5603 -19076.5837
394 -96716.4176 74584.5603
395 6467.0060 -96716.4176
396 -138660.7052 6467.0060
397 -189544.2257 -138660.7052
398 12385.3742 -189544.2257
399 -55495.2182 12385.3742
400 2436.6634 -55495.2182
401 42617.6282 2436.6634
402 19587.8792 42617.6282
403 -49563.0225 19587.8792
404 -8991.9428 -49563.0225
405 14014.8986 -8991.9428
406 20621.8119 14014.8986
407 13253.5098 20621.8119
408 22909.7443 13253.5098
409 -61897.0645 22909.7443
410 -51363.8727 -61897.0645
411 -161762.7527 -51363.8727
412 -25904.1719 -161762.7527
413 60090.8563 -25904.1719
414 -177521.3404 60090.8563
415 -99190.4423 -177521.3404
416 91986.7555 -99190.4423
417 -17377.1556 91986.7555
418 -50407.0390 -17377.1556
419 -37399.3791 -50407.0390
420 -44260.1342 -37399.3791
421 -175397.0259 -44260.1342
422 -547868.0565 -175397.0259
423 -458554.3229 -547868.0565
424 -48841.8151 -458554.3229
425 -250017.3396 -48841.8151
426 -311854.7009 -250017.3396
427 -1135952.0870 -311854.7009
428 -272140.7589 -1135952.0870
429 -336046.0989 -272140.7589
430 -682476.5187 -336046.0989
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7qdsz1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8qdsz1292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/904921292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/rcomp/tmp/1004921292936583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11m5qq1292936583.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/12p5ow1292936583.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13w73p1292936583.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/14hpkv1292936583.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/15sg1y1292936583.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/rcomp/tmp/16oqh71292936583.tab")
+ }
>
> try(system("convert tmp/1ulcq1292936583.ps tmp/1ulcq1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/2mvcb1292936583.ps tmp/2mvcb1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/3mvcb1292936583.ps tmp/3mvcb1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/4mvcb1292936583.ps tmp/4mvcb1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/5f4bw1292936583.ps tmp/5f4bw1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/6f4bw1292936583.ps tmp/6f4bw1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/7qdsz1292936583.ps tmp/7qdsz1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/8qdsz1292936583.ps tmp/8qdsz1292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/904921292936583.ps tmp/904921292936583.png",intern=TRUE))
character(0)
> try(system("convert tmp/1004921292936583.ps tmp/1004921292936583.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.661 2.210 23.947