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 = '2'
> #'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
Orders Costs Dividends Wealth
1 807 162556 213118 6282154
2 444 29790 81767 4321023
3 412 87550 153198 4111912
4 428 84738 -26007 223193
5 315 54660 126942 1491348
6 168 42634 157214 1629616
7 263 40949 129352 1398893
8 267 45187 234817 1926517
9 228 37704 60448 983660
10 129 16275 47818 1443586
11 104 25830 245546 1073089
12 122 12679 48020 984885
13 393 18014 -1710 1405225
14 190 43556 32648 227132
15 280 24811 95350 929118
16 63 6575 151352 1071292
17 102 7123 288170 638830
18 265 21950 114337 856956
19 234 37597 37884 992426
20 277 17821 122844 444477
21 73 12988 82340 857217
22 67 22330 79801 711969
23 103 13326 165548 702380
24 290 16189 116384 358589
25 83 7146 134028 297978
26 56 15824 63838 585715
27 236 27664 74996 657954
28 73 11920 31080 209458
29 34 8568 32168 786690
30 139 14416 49857 439798
31 26 3369 87161 688779
32 70 11819 106113 574339
33 40 6984 80570 741409
34 42 4519 102129 597793
35 12 2220 301670 644190
36 211 18562 102313 377934
37 74 10327 88577 640273
38 80 5336 112477 697458
39 83 2365 191778 550608
40 131 4069 79804 207393
41 203 8636 128294 301607
42 56 13718 96448 345783
43 89 4525 93811 501749
44 88 6869 117520 379983
45 39 4628 69159 387475
46 25 3689 101792 377305
47 49 4891 210568 370837
48 149 7489 136996 430866
49 58 4901 121920 469107
50 41 2284 76403 194493
51 90 3160 108094 530670
52 136 4150 134759 518365
53 97 7285 188873 491303
54 63 1134 146216 527021
55 114 4658 156608 233773
56 77 2384 61348 405972
57 6 3748 50350 652925
58 47 5371 87720 446211
59 51 1285 99489 341340
60 85 9327 87419 387699
61 43 5565 94355 493408
62 32 1528 60326 146494
63 25 3122 94670 414462
64 77 7561 82425 364304
65 54 2675 59017 355178
66 251 13253 90829 357760
67 15 880 80791 261216
68 44 2053 100423 397144
69 73 1424 131116 374943
70 85 4036 100269 424898
71 49 3045 27330 202055
72 38 5119 39039 378525
73 35 1431 106885 310768
74 9 554 79285 325738
75 34 1975 118881 394510
76 20 1765 77623 247060
77 29 1012 114768 368078
78 11 810 74015 236761
79 52 1280 69465 312378
80 13 666 117869 339836
81 29 1380 60982 347385
82 66 4677 90131 426280
83 33 876 138971 352850
84 15 814 39625 301881
85 15 514 102725 377516
86 68 5692 64239 357312
87 100 3642 90262 458343
88 13 540 103960 354228
89 45 2099 106611 308636
90 14 567 103345 386212
91 36 2001 95551 393343
92 40 2949 82903 378509
93 68 2253 63593 452469
94 29 6533 126910 364839
95 43 1889 37527 358649
96 30 3055 60247 376641
97 9 272 112995 429112
98 22 1414 70184 330546
99 19 2564 130140 403560
100 9 1383 73221 317892
101 31 1261 76114 307528
102 19 975 90534 235133
103 55 3366 108479 299243
104 8 576 113761 314073
105 28 1686 68696 368186
106 29 746 71561 269661
107 48 3192 59831 125390
108 16 2045 97890 510834
109 47 5702 101481 321896
110 20 1932 72954 249898
111 22 936 67939 408881
112 33 3437 48022 158492
113 44 5131 86111 292154
114 13 2397 74020 289513
115 6 1389 57530 378049
116 35 1503 56364 343466
117 8 402 84990 332743
118 17 2239 88590 442882
119 11 2234 77200 214215
120 21 837 61262 315688
121 92 10579 110309 375195
122 12 875 67000 334280
123 112 1585 93099 355864
124 25 1659 107577 480382
125 17 2647 62920 353058
126 23 3294 75832 217193
127 0 0 60720 315380
128 10 94 60793 314533
129 23 422 57935 318056
130 0 0 60720 315380
131 7 34 60630 314353
132 25 1558 55637 369448
133 1 0 60720 315380
134 20 43 60887 312846
135 4 645 60720 312075
136 4 316 60505 315009
137 10 115 60945 318903
138 1 5 60720 314887
139 4 897 60720 314913
140 0 0 60720 315380
141 8 389 58990 325506
142 0 0 60720 315380
143 11 1002 56750 298568
144 4 36 60894 315834
145 15 460 63346 329784
146 9 309 56535 312878
147 0 0 60720 315380
148 7 9 60835 314987
149 2 271 60720 325249
150 0 14 61016 315877
151 7 520 58650 291650
152 46 1766 60438 305959
153 5 0 60720 315380
154 7 458 58625 297765
155 2 20 60938 315245
156 0 0 60720 315380
157 0 0 60720 315380
158 2 98 61490 315236
159 5 405 60845 336425
160 0 0 60720 315380
161 0 0 60720 315380
162 0 0 60720 315380
163 0 0 60720 315380
164 7 483 60830 306268
165 24 454 63261 302187
166 1 47 60720 314882
167 0 0 60720 315380
168 18 757 45689 382712
169 55 4655 60720 341570
170 0 0 60720 315380
171 0 0 60720 315380
172 3 36 61564 312412
173 0 0 60720 315380
174 9 203 61938 309596
175 0 0 60720 315380
176 8 126 60951 315547
177 113 400 60720 313267
178 0 71 60745 316176
179 0 0 60720 315380
180 0 0 60720 315380
181 19 972 71642 359335
182 11 531 71641 330068
183 25 2461 55792 314289
184 16 378 71873 297413
185 5 23 62555 314806
186 11 638 60370 333210
187 23 2300 64873 352108
188 6 149 62041 313332
189 5 226 65745 291787
190 0 0 60720 315380
191 7 275 59500 318745
192 0 0 60720 315380
193 7 141 61630 315366
194 0 0 60720 315380
195 3 28 60890 315688
196 0 0 60720 315380
197 89 4980 113521 409642
198 0 0 60720 315380
199 0 0 60720 315380
200 19 472 80045 269587
201 0 0 60720 315380
202 0 0 60720 315380
203 0 0 60720 315380
204 12 203 50804 300962
205 12 496 87390 325479
206 5 10 61656 316155
207 2 63 65688 318574
208 0 0 60720 315380
209 26 1136 48522 343613
210 3 265 60720 306948
211 0 0 60720 315380
212 0 0 60720 315380
213 11 267 57640 330059
214 10 474 61977 288985
215 5 534 62620 304485
216 2 0 60720 315380
217 6 15 60831 315688
218 7 397 60646 317736
219 2 0 60720 315380
220 28 1866 56225 322331
221 3 288 60510 296656
222 0 0 60720 315380
223 1 3 60698 315354
224 20 468 60720 312161
225 1 20 60805 315576
226 22 278 61404 314922
227 9 61 60720 314551
228 0 0 60720 315380
229 2 192 65276 312339
230 0 0 60720 315380
231 7 317 63915 298700
232 9 738 60720 321376
233 0 0 60720 315380
234 13 368 61686 303230
235 0 0 60720 315380
236 0 2 60743 315487
237 0 0 60720 315380
238 6 53 60349 315793
239 0 0 60720 315380
240 0 0 60720 315380
241 0 0 60720 315380
242 3 94 61360 312887
243 0 0 60720 315380
244 7 24 59818 315637
245 2 2332 72680 324385
246 0 0 60720 315380
247 0 0 60720 315380
248 15 131 61808 308989
249 0 0 60720 315380
250 0 0 60720 315380
251 9 206 53110 296702
252 0 0 60720 315380
253 1 167 64245 307322
254 38 622 73007 304376
255 57 2328 82732 253588
256 0 0 60720 315380
257 7 365 54820 309560
258 26 364 47705 298466
259 0 0 60720 315380
260 0 0 60720 315380
261 0 0 60720 315380
262 0 0 60720 315380
263 13 226 72835 343929
264 10 307 58856 331955
265 0 0 60720 315380
266 0 0 60720 315380
267 0 0 60720 315380
268 9 188 77655 381180
269 0 0 60720 315380
270 26 138 69817 331420
271 0 0 60720 315380
272 0 0 60720 315380
273 0 0 60720 315380
274 19 125 60798 310201
275 0 0 60720 315380
276 12 282 62452 320016
277 23 335 64175 320398
278 0 0 60720 315380
279 29 1324 67440 291841
280 8 176 68136 310670
281 0 0 60720 315380
282 0 0 60720 315380
283 26 249 56726 313491
284 0 0 60720 315380
285 9 333 70811 331323
286 0 0 60720 315380
287 5 601 60720 319210
288 3 30 62045 318098
289 0 0 60720 315380
290 13 249 54323 292754
291 0 0 60720 315380
292 12 165 62841 325176
293 19 453 81125 365959
294 0 0 60720 315380
295 10 53 59506 302409
296 9 382 59365 340968
297 0 0 60720 315380
298 0 0 60720 315380
299 0 0 60720 315380
300 9 0 60720 315380
301 4 30 60798 313164
302 1 290 58790 301164
303 1 0 60720 315380
304 0 0 60720 315380
305 14 366 61808 344425
306 12 2 60735 315394
307 0 0 60720 315380
308 19 209 64016 316647
309 17 384 54683 309836
310 0 0 60720 315380
311 0 0 60720 315380
312 32 365 87192 346611
313 0 0 60720 315380
314 14 49 64107 322031
315 8 3 60761 315656
316 4 133 65990 339445
317 0 32 59988 314964
318 20 368 61167 297141
319 5 1 60719 315372
320 0 0 60720 315380
321 0 0 60720 315380
322 0 0 60720 315380
323 0 0 60720 315380
324 0 0 60720 315380
325 0 0 60720 315380
326 1 22 60722 312502
327 0 0 60720 315380
328 0 0 60720 315380
329 0 0 60720 315380
330 0 0 60720 315380
331 0 0 60720 315380
332 0 0 60720 315380
333 0 0 60720 315380
334 4 96 60379 313729
335 1 1 60727 315388
336 4 314 60720 315371
337 20 844 60925 296139
338 0 0 60720 315380
339 1 26 60896 313880
340 10 125 59734 317698
341 12 304 62969 295580
342 0 0 60720 315380
343 0 0 60720 315380
344 0 0 60720 315380
345 13 621 60720 308256
346 0 0 60720 315380
347 3 119 59118 303677
348 0 0 60720 315380
349 0 0 60720 315380
350 10 1595 60720 319369
351 3 312 58598 318690
352 7 60 61124 314049
353 10 587 59595 325699
354 1 135 62065 314210
355 0 0 60720 315380
356 0 0 60720 315380
357 15 514 78780 322378
358 0 0 60720 315380
359 0 0 60720 315380
360 0 0 60720 315380
361 4 1 60722 315398
362 0 0 60720 315380
363 0 0 60720 315380
364 28 1763 61600 308336
365 9 180 59635 316386
366 0 0 60720 315380
367 0 0 60720 315380
368 0 0 60720 315380
369 0 0 60720 315380
370 7 218 60720 315553
371 0 0 60720 315380
372 7 448 59781 323361
373 7 227 76644 336639
374 3 174 64820 307424
375 0 0 60720 315380
376 0 0 60720 315380
377 11 121 56178 295370
378 7 607 60436 322340
379 10 2212 60720 319864
380 0 0 60720 315380
381 0 0 60720 315380
382 18 530 73433 317291
383 14 571 41477 280398
384 0 0 60720 315380
385 12 78 62700 317330
386 29 2489 67804 238125
387 3 131 59661 327071
388 6 923 58620 309038
389 3 72 60398 314210
390 8 572 58580 307930
391 10 397 62710 322327
392 6 450 59325 292136
393 8 622 60950 263276
394 6 694 68060 367655
395 9 3425 83620 283910
396 8 562 58456 283587
397 26 4917 52811 243650
398 239 1442 121173 438493
399 7 529 63870 296261
400 41 2126 21001 230621
401 3 1061 70415 304252
402 8 776 64230 333505
403 6 611 59190 296919
404 21 1526 69351 278990
405 7 592 64270 276898
406 11 1182 70694 327007
407 11 621 68005 317046
408 12 989 58930 304555
409 9 438 58320 298096
410 3 726 69980 231861
411 57 1303 69863 309422
412 21 7419 63255 286963
413 15 1164 57320 269753
414 32 3310 75230 448243
415 11 1920 79420 165404
416 2 965 73490 204325
417 23 3256 35250 407159
418 20 1135 62285 290476
419 24 1270 69206 275311
420 1 661 65920 246541
421 1 1013 69770 253468
422 74 2844 72683 240897
423 68 11528 -14545 -83265
424 20 6526 55830 -42143
425 20 2264 55174 272713
426 82 5109 67038 215362
427 21 3999 51252 42754
428 244 35624 157278 306275
429 32 9252 79510 253537
430 86 15236 77440 372631
431 69 18073 27284 -7170
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Costs Dividends Wealth
-5.927e+00 4.792e-03 2.134e-04 1.896e-05
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-130.625 -13.010 -9.128 1.221 286.324
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -5.927e+00 4.329e+00 -1.369 0.171719
Costs 4.792e-03 2.340e-04 20.477 < 2e-16 ***
Dividends 2.134e-04 5.641e-05 3.783 0.000177 ***
Wealth 1.896e-05 6.569e-06 2.887 0.004092 **
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 34.75 on 427 degrees of freedom
Multiple R-squared: 0.7839, Adjusted R-squared: 0.7824
F-statistic: 516.3 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,] 0.9982348 3.530380e-03 1.765190e-03
[2,] 0.9998106 3.788738e-04 1.894369e-04
[3,] 0.9997402 5.195682e-04 2.597841e-04
[4,] 0.9999833 3.342067e-05 1.671033e-05
[5,] 0.9999848 3.035395e-05 1.517698e-05
[6,] 0.9999842 3.151213e-05 1.575607e-05
[7,] 1.0000000 1.588233e-11 7.941166e-12
[8,] 1.0000000 1.766135e-12 8.830675e-13
[9,] 1.0000000 1.601033e-14 8.005167e-15
[10,] 1.0000000 1.512065e-14 7.560323e-15
[11,] 1.0000000 2.447820e-16 1.223910e-16
[12,] 1.0000000 2.959264e-18 1.479632e-18
[13,] 1.0000000 4.413380e-18 2.206690e-18
[14,] 1.0000000 7.104719e-23 3.552359e-23
[15,] 1.0000000 5.737978e-25 2.868989e-25
[16,] 1.0000000 3.051290e-29 1.525645e-29
[17,] 1.0000000 5.699245e-29 2.849622e-29
[18,] 1.0000000 5.248810e-37 2.624405e-37
[19,] 1.0000000 1.519599e-36 7.597997e-37
[20,] 1.0000000 7.895915e-40 3.947958e-40
[21,] 1.0000000 1.375510e-39 6.877550e-40
[22,] 1.0000000 2.445891e-40 1.222946e-40
[23,] 1.0000000 3.075644e-43 1.537822e-43
[24,] 1.0000000 8.130189e-43 4.065095e-43
[25,] 1.0000000 8.139010e-44 4.069505e-44
[26,] 1.0000000 3.252055e-44 1.626028e-44
[27,] 1.0000000 3.212196e-45 1.606098e-45
[28,] 1.0000000 2.612578e-45 1.306289e-45
[29,] 1.0000000 2.005803e-47 1.002901e-47
[30,] 1.0000000 8.337355e-49 4.168678e-49
[31,] 1.0000000 7.471390e-49 3.735695e-49
[32,] 1.0000000 2.504666e-48 1.252333e-48
[33,] 1.0000000 6.909073e-48 3.454537e-48
[34,] 1.0000000 5.505645e-50 2.752823e-50
[35,] 1.0000000 1.777326e-55 8.886629e-56
[36,] 1.0000000 6.343819e-57 3.171909e-57
[37,] 1.0000000 1.399601e-56 6.998006e-57
[38,] 1.0000000 5.161566e-56 2.580783e-56
[39,] 1.0000000 3.925506e-56 1.962753e-56
[40,] 1.0000000 1.967231e-56 9.836156e-57
[41,] 1.0000000 8.463077e-57 4.231539e-57
[42,] 1.0000000 6.643086e-58 3.321543e-58
[43,] 1.0000000 1.537242e-57 7.686211e-58
[44,] 1.0000000 2.473981e-57 1.236991e-57
[45,] 1.0000000 4.066066e-57 2.033033e-57
[46,] 1.0000000 1.515954e-58 7.579770e-59
[47,] 1.0000000 3.948550e-58 1.974275e-58
[48,] 1.0000000 1.400942e-57 7.004708e-58
[49,] 1.0000000 4.211118e-58 2.105559e-58
[50,] 1.0000000 1.738842e-58 8.694210e-59
[51,] 1.0000000 2.188622e-60 1.094311e-60
[52,] 1.0000000 2.848182e-60 1.424091e-60
[53,] 1.0000000 5.834498e-60 2.917249e-60
[54,] 1.0000000 1.876161e-59 9.380807e-60
[55,] 1.0000000 1.472232e-59 7.361162e-60
[56,] 1.0000000 1.522912e-59 7.614560e-60
[57,] 1.0000000 1.149525e-59 5.747626e-60
[58,] 1.0000000 3.591374e-59 1.795687e-59
[59,] 1.0000000 5.351593e-59 2.675797e-59
[60,] 1.0000000 5.287666e-71 2.643833e-71
[61,] 1.0000000 6.685561e-71 3.342781e-71
[62,] 1.0000000 1.825038e-70 9.125191e-71
[63,] 1.0000000 2.055132e-70 1.027566e-70
[64,] 1.0000000 1.867579e-70 9.337895e-71
[65,] 1.0000000 7.429756e-71 3.714878e-71
[66,] 1.0000000 9.031350e-71 4.515675e-71
[67,] 1.0000000 2.420535e-70 1.210268e-70
[68,] 1.0000000 2.558882e-70 1.279441e-70
[69,] 1.0000000 6.313101e-70 3.156551e-70
[70,] 1.0000000 1.079854e-69 5.399269e-70
[71,] 1.0000000 2.785265e-69 1.392632e-69
[72,] 1.0000000 4.096927e-69 2.048464e-69
[73,] 1.0000000 3.331330e-69 1.665665e-69
[74,] 1.0000000 4.375153e-69 2.187576e-69
[75,] 1.0000000 8.081618e-69 4.040809e-69
[76,] 1.0000000 2.184961e-68 1.092481e-68
[77,] 1.0000000 6.634984e-68 3.317492e-68
[78,] 1.0000000 8.438484e-68 4.219242e-68
[79,] 1.0000000 1.392804e-67 6.964019e-68
[80,] 1.0000000 2.449459e-67 1.224729e-67
[81,] 1.0000000 1.497267e-68 7.486335e-69
[82,] 1.0000000 2.285478e-68 1.142739e-68
[83,] 1.0000000 6.951083e-68 3.475542e-68
[84,] 1.0000000 1.073633e-67 5.368165e-68
[85,] 1.0000000 3.054618e-67 1.527309e-67
[86,] 1.0000000 8.252023e-67 4.126011e-67
[87,] 1.0000000 2.693377e-67 1.346688e-67
[88,] 1.0000000 3.944519e-68 1.972259e-68
[89,] 1.0000000 2.088630e-68 1.044315e-68
[90,] 1.0000000 4.283686e-68 2.141843e-68
[91,] 1.0000000 3.363597e-68 1.681798e-68
[92,] 1.0000000 7.608381e-68 3.804191e-68
[93,] 1.0000000 2.037087e-68 1.018544e-68
[94,] 1.0000000 2.831503e-68 1.415752e-68
[95,] 1.0000000 7.457111e-68 3.728556e-68
[96,] 1.0000000 1.789345e-67 8.946726e-68
[97,] 1.0000000 6.237970e-67 3.118985e-67
[98,] 1.0000000 3.138240e-67 1.569120e-67
[99,] 1.0000000 7.578249e-67 3.789124e-67
[100,] 1.0000000 1.858245e-66 9.291226e-67
[101,] 1.0000000 2.382560e-66 1.191280e-66
[102,] 1.0000000 1.678983e-66 8.394917e-67
[103,] 1.0000000 3.634883e-66 1.817442e-66
[104,] 1.0000000 8.738220e-66 4.369110e-66
[105,] 1.0000000 2.093571e-65 1.046786e-65
[106,] 1.0000000 3.556209e-65 1.778104e-65
[107,] 1.0000000 1.017864e-64 5.089320e-65
[108,] 1.0000000 1.528917e-64 7.644584e-65
[109,] 1.0000000 2.329265e-64 1.164633e-64
[110,] 1.0000000 3.406101e-64 1.703050e-64
[111,] 1.0000000 5.203742e-64 2.601871e-64
[112,] 1.0000000 5.742088e-64 2.871044e-64
[113,] 1.0000000 8.516681e-64 4.258341e-64
[114,] 1.0000000 1.997791e-63 9.988954e-64
[115,] 1.0000000 6.168347e-63 3.084174e-63
[116,] 1.0000000 1.382126e-62 6.910629e-63
[117,] 1.0000000 3.821532e-66 1.910766e-66
[118,] 1.0000000 4.842284e-66 2.421142e-66
[119,] 1.0000000 1.006642e-65 5.033209e-66
[120,] 1.0000000 2.373082e-65 1.186541e-65
[121,] 1.0000000 4.433866e-65 2.216933e-65
[122,] 1.0000000 1.073637e-64 5.368186e-65
[123,] 1.0000000 2.109796e-64 1.054898e-64
[124,] 1.0000000 4.076799e-64 2.038400e-64
[125,] 1.0000000 9.721224e-64 4.860612e-64
[126,] 1.0000000 2.150096e-63 1.075048e-63
[127,] 1.0000000 4.450786e-63 2.225393e-63
[128,] 1.0000000 9.865985e-63 4.932993e-63
[129,] 1.0000000 2.076681e-62 1.038341e-62
[130,] 1.0000000 4.638756e-62 2.319378e-62
[131,] 1.0000000 1.177592e-61 5.887961e-62
[132,] 1.0000000 2.534487e-61 1.267243e-61
[133,] 1.0000000 5.246779e-61 2.623389e-61
[134,] 1.0000000 1.108694e-60 5.543470e-61
[135,] 1.0000000 2.758609e-60 1.379305e-60
[136,] 1.0000000 5.910029e-60 2.955014e-60
[137,] 1.0000000 1.496616e-59 7.483079e-60
[138,] 1.0000000 3.619837e-59 1.809918e-59
[139,] 1.0000000 9.720389e-59 4.860194e-59
[140,] 1.0000000 2.452340e-58 1.226170e-58
[141,] 1.0000000 5.377179e-58 2.688590e-58
[142,] 1.0000000 1.384565e-57 6.922824e-58
[143,] 1.0000000 3.107500e-57 1.553750e-57
[144,] 1.0000000 6.863367e-57 3.431683e-57
[145,] 1.0000000 1.743476e-56 8.717378e-57
[146,] 1.0000000 1.413780e-56 7.068898e-57
[147,] 1.0000000 3.575946e-56 1.787973e-56
[148,] 1.0000000 9.118855e-56 4.559428e-56
[149,] 1.0000000 2.165974e-55 1.082987e-55
[150,] 1.0000000 4.892129e-55 2.446065e-55
[151,] 1.0000000 1.108384e-54 5.541919e-55
[152,] 1.0000000 2.614059e-54 1.307030e-54
[153,] 1.0000000 6.368322e-54 3.184161e-54
[154,] 1.0000000 1.452670e-53 7.263350e-54
[155,] 1.0000000 3.320946e-53 1.660473e-53
[156,] 1.0000000 7.607087e-53 3.803544e-53
[157,] 1.0000000 1.745619e-52 8.728093e-53
[158,] 1.0000000 4.470249e-52 2.235124e-52
[159,] 1.0000000 9.923499e-52 4.961750e-52
[160,] 1.0000000 2.329486e-51 1.164743e-51
[161,] 1.0000000 5.356619e-51 2.678309e-51
[162,] 1.0000000 1.005883e-50 5.029415e-51
[163,] 1.0000000 1.368181e-50 6.840906e-51
[164,] 1.0000000 3.154729e-50 1.577364e-50
[165,] 1.0000000 7.279591e-50 3.639795e-50
[166,] 1.0000000 1.787881e-49 8.939406e-50
[167,] 1.0000000 4.127696e-49 2.063848e-49
[168,] 1.0000000 1.081393e-48 5.406964e-49
[169,] 1.0000000 2.494917e-48 1.247459e-48
[170,] 1.0000000 6.468277e-48 3.234139e-48
[171,] 1.0000000 5.120789e-55 2.560395e-55
[172,] 1.0000000 1.253448e-54 6.267241e-55
[173,] 1.0000000 3.103329e-54 1.551665e-54
[174,] 1.0000000 7.683765e-54 3.841883e-54
[175,] 1.0000000 2.241164e-53 1.120582e-53
[176,] 1.0000000 6.033518e-53 3.016759e-53
[177,] 1.0000000 1.502999e-52 7.514993e-53
[178,] 1.0000000 4.367340e-52 2.183670e-52
[179,] 1.0000000 1.173511e-51 5.867557e-52
[180,] 1.0000000 3.211067e-51 1.605534e-51
[181,] 1.0000000 9.027439e-51 4.513720e-51
[182,] 1.0000000 2.426183e-50 1.213092e-50
[183,] 1.0000000 6.232851e-50 3.116425e-50
[184,] 1.0000000 1.521008e-49 7.605039e-50
[185,] 1.0000000 4.076696e-49 2.038348e-49
[186,] 1.0000000 9.925991e-49 4.962995e-49
[187,] 1.0000000 2.666625e-48 1.333313e-48
[188,] 1.0000000 6.475185e-48 3.237593e-48
[189,] 1.0000000 1.666113e-47 8.330563e-48
[190,] 1.0000000 4.034962e-47 2.017481e-47
[191,] 1.0000000 5.044490e-47 2.522245e-47
[192,] 1.0000000 1.220001e-46 6.100006e-47
[193,] 1.0000000 2.946672e-46 1.473336e-46
[194,] 1.0000000 8.215387e-46 4.107694e-46
[195,] 1.0000000 1.974970e-45 9.874852e-46
[196,] 1.0000000 4.740755e-45 2.370378e-45
[197,] 1.0000000 1.136212e-44 5.681060e-45
[198,] 1.0000000 2.480295e-44 1.240148e-44
[199,] 1.0000000 4.843811e-44 2.421906e-44
[200,] 1.0000000 1.250937e-43 6.254684e-44
[201,] 1.0000000 2.914489e-43 1.457245e-43
[202,] 1.0000000 6.926240e-43 3.463120e-43
[203,] 1.0000000 8.825934e-43 4.412967e-43
[204,] 1.0000000 2.177320e-42 1.088660e-42
[205,] 1.0000000 5.165025e-42 2.582512e-42
[206,] 1.0000000 1.222805e-41 6.114024e-42
[207,] 1.0000000 3.022787e-41 1.511393e-41
[208,] 1.0000000 7.829506e-41 3.914753e-41
[209,] 1.0000000 1.900291e-40 9.501456e-41
[210,] 1.0000000 4.634540e-40 2.317270e-40
[211,] 1.0000000 1.174343e-39 5.871714e-40
[212,] 1.0000000 2.952201e-39 1.476101e-39
[213,] 1.0000000 7.138251e-39 3.569125e-39
[214,] 1.0000000 1.317443e-38 6.587213e-39
[215,] 1.0000000 3.164652e-38 1.582326e-38
[216,] 1.0000000 7.323336e-38 3.661668e-38
[217,] 1.0000000 1.724476e-37 8.622382e-38
[218,] 1.0000000 3.675468e-37 1.837734e-37
[219,] 1.0000000 8.592255e-37 4.296128e-37
[220,] 1.0000000 1.673135e-36 8.365675e-37
[221,] 1.0000000 4.110628e-36 2.055314e-36
[222,] 1.0000000 9.384328e-36 4.692164e-36
[223,] 1.0000000 2.065507e-35 1.032753e-35
[224,] 1.0000000 4.690526e-35 2.345263e-35
[225,] 1.0000000 1.129985e-34 5.649925e-35
[226,] 1.0000000 2.729701e-34 1.364850e-34
[227,] 1.0000000 6.143585e-34 3.071792e-34
[228,] 1.0000000 1.464666e-33 7.323331e-34
[229,] 1.0000000 3.277330e-33 1.638665e-33
[230,] 1.0000000 7.309171e-33 3.654585e-33
[231,] 1.0000000 1.626113e-32 8.130567e-33
[232,] 1.0000000 3.853634e-32 1.926817e-32
[233,] 1.0000000 8.521092e-32 4.260546e-32
[234,] 1.0000000 1.878538e-31 9.392691e-32
[235,] 1.0000000 4.128861e-31 2.064430e-31
[236,] 1.0000000 9.375418e-31 4.687709e-31
[237,] 1.0000000 2.047621e-30 1.023811e-30
[238,] 1.0000000 4.755825e-30 2.377913e-30
[239,] 1.0000000 5.007747e-30 2.503873e-30
[240,] 1.0000000 1.087612e-29 5.438061e-30
[241,] 1.0000000 2.354563e-29 1.177282e-29
[242,] 1.0000000 5.181283e-29 2.590641e-29
[243,] 1.0000000 1.114833e-28 5.574166e-29
[244,] 1.0000000 2.390853e-28 1.195427e-28
[245,] 1.0000000 5.140044e-28 2.570022e-28
[246,] 1.0000000 1.095793e-27 5.478963e-28
[247,] 1.0000000 2.239302e-27 1.119651e-27
[248,] 1.0000000 3.268411e-27 1.634206e-27
[249,] 1.0000000 3.782429e-27 1.891215e-27
[250,] 1.0000000 8.017740e-27 4.008870e-27
[251,] 1.0000000 1.774164e-26 8.870821e-27
[252,] 1.0000000 1.720013e-26 8.600067e-27
[253,] 1.0000000 3.636534e-26 1.818267e-26
[254,] 1.0000000 7.661660e-26 3.830830e-26
[255,] 1.0000000 1.608514e-25 8.042568e-26
[256,] 1.0000000 3.364968e-25 1.682484e-25
[257,] 1.0000000 7.448833e-25 3.724416e-25
[258,] 1.0000000 1.618162e-24 8.090809e-25
[259,] 1.0000000 3.347615e-24 1.673808e-24
[260,] 1.0000000 6.900169e-24 3.450085e-24
[261,] 1.0000000 1.417041e-23 7.085207e-24
[262,] 1.0000000 2.739323e-23 1.369661e-23
[263,] 1.0000000 5.583270e-23 2.791635e-23
[264,] 1.0000000 1.031040e-22 5.155202e-23
[265,] 1.0000000 2.089033e-22 1.044517e-22
[266,] 1.0000000 4.216589e-22 2.108295e-22
[267,] 1.0000000 8.478355e-22 4.239178e-22
[268,] 1.0000000 1.563417e-21 7.817083e-22
[269,] 1.0000000 3.124630e-21 1.562315e-21
[270,] 1.0000000 6.540469e-21 3.270234e-21
[271,] 1.0000000 1.179464e-20 5.897319e-21
[272,] 1.0000000 2.333991e-20 1.166995e-20
[273,] 1.0000000 4.308600e-20 2.154300e-20
[274,] 1.0000000 8.777971e-20 4.388986e-20
[275,] 1.0000000 1.718004e-19 8.590018e-20
[276,] 1.0000000 3.348860e-19 1.674430e-19
[277,] 1.0000000 4.430779e-19 2.215390e-19
[278,] 1.0000000 8.610882e-19 4.305441e-19
[279,] 1.0000000 1.693535e-18 8.467677e-19
[280,] 1.0000000 3.263284e-18 1.631642e-18
[281,] 1.0000000 6.407520e-18 3.203760e-18
[282,] 1.0000000 1.254245e-17 6.271223e-18
[283,] 1.0000000 2.385953e-17 1.192977e-17
[284,] 1.0000000 4.330813e-17 2.165407e-17
[285,] 1.0000000 8.180411e-17 4.090206e-17
[286,] 1.0000000 1.601745e-16 8.008725e-17
[287,] 1.0000000 3.118189e-16 1.559094e-16
[288,] 1.0000000 5.811872e-16 2.905936e-16
[289,] 1.0000000 1.112268e-15 5.561342e-16
[290,] 1.0000000 2.140147e-15 1.070074e-15
[291,] 1.0000000 3.938176e-15 1.969088e-15
[292,] 1.0000000 7.214428e-15 3.607214e-15
[293,] 1.0000000 1.315682e-14 6.578411e-15
[294,] 1.0000000 2.483252e-14 1.241626e-14
[295,] 1.0000000 4.639858e-14 2.319929e-14
[296,] 1.0000000 8.420475e-14 4.210238e-14
[297,] 1.0000000 1.524302e-13 7.621512e-14
[298,] 1.0000000 2.716597e-13 1.358299e-13
[299,] 1.0000000 4.952269e-13 2.476135e-13
[300,] 1.0000000 8.924292e-13 4.462146e-13
[301,] 1.0000000 1.569698e-12 7.848492e-13
[302,] 1.0000000 2.698681e-12 1.349341e-12
[303,] 1.0000000 4.319144e-12 2.159572e-12
[304,] 1.0000000 7.514713e-12 3.757356e-12
[305,] 1.0000000 1.301089e-11 6.505444e-12
[306,] 1.0000000 2.292172e-11 1.146086e-11
[307,] 1.0000000 3.930495e-11 1.965247e-11
[308,] 1.0000000 6.858879e-11 3.429440e-11
[309,] 1.0000000 1.204359e-10 6.021794e-11
[310,] 1.0000000 2.034082e-10 1.017041e-10
[311,] 1.0000000 3.428875e-10 1.714437e-10
[312,] 1.0000000 5.454366e-10 2.727183e-10
[313,] 1.0000000 9.367951e-10 4.683976e-10
[314,] 1.0000000 1.553456e-09 7.767279e-10
[315,] 1.0000000 2.562618e-09 1.281309e-09
[316,] 1.0000000 4.205170e-09 2.102585e-09
[317,] 1.0000000 6.864068e-09 3.432034e-09
[318,] 1.0000000 1.114453e-08 5.572264e-09
[319,] 1.0000000 1.799730e-08 8.998650e-09
[320,] 1.0000000 2.912784e-08 1.456392e-08
[321,] 1.0000000 4.652786e-08 2.326393e-08
[322,] 1.0000000 7.391512e-08 3.695756e-08
[323,] 0.9999999 1.167755e-07 5.838776e-08
[324,] 0.9999999 1.834643e-07 9.173215e-08
[325,] 0.9999999 2.866257e-07 1.433129e-07
[326,] 0.9999998 4.452718e-07 2.226359e-07
[327,] 0.9999997 6.878033e-07 3.439016e-07
[328,] 0.9999995 1.081409e-06 5.407047e-07
[329,] 0.9999992 1.663716e-06 8.318581e-07
[330,] 0.9999987 2.567589e-06 1.283794e-06
[331,] 0.9999981 3.866945e-06 1.933472e-06
[332,] 0.9999971 5.806983e-06 2.903492e-06
[333,] 0.9999956 8.718161e-06 4.359080e-06
[334,] 0.9999934 1.328884e-05 6.644420e-06
[335,] 0.9999899 2.017675e-05 1.008837e-05
[336,] 0.9999852 2.958701e-05 1.479350e-05
[337,] 0.9999784 4.312090e-05 2.156045e-05
[338,] 0.9999688 6.245861e-05 3.122931e-05
[339,] 0.9999538 9.231621e-05 4.615810e-05
[340,] 0.9999340 1.320790e-04 6.603949e-05
[341,] 0.9999044 1.912683e-04 9.563415e-05
[342,] 0.9998649 2.701715e-04 1.350858e-04
[343,] 0.9998104 3.791864e-04 1.895932e-04
[344,] 0.9997326 5.347336e-04 2.673668e-04
[345,] 0.9996243 7.513796e-04 3.756898e-04
[346,] 0.9994700 1.060030e-03 5.300149e-04
[347,] 0.9992571 1.485813e-03 7.429067e-04
[348,] 0.9989944 2.011191e-03 1.005595e-03
[349,] 0.9986437 2.712575e-03 1.356288e-03
[350,] 0.9981830 3.633996e-03 1.816998e-03
[351,] 0.9975754 4.849102e-03 2.424551e-03
[352,] 0.9967959 6.408232e-03 3.204116e-03
[353,] 0.9957946 8.410789e-03 4.205395e-03
[354,] 0.9945184 1.096325e-02 5.481627e-03
[355,] 0.9927849 1.443016e-02 7.215080e-03
[356,] 0.9907272 1.854568e-02 9.272838e-03
[357,] 0.9881658 2.366836e-02 1.183418e-02
[358,] 0.9849628 3.007431e-02 1.503716e-02
[359,] 0.9806262 3.874767e-02 1.937383e-02
[360,] 0.9757661 4.846790e-02 2.423395e-02
[361,] 0.9699038 6.019245e-02 3.009622e-02
[362,] 0.9628919 7.421618e-02 3.710809e-02
[363,] 0.9545760 9.084796e-02 4.542398e-02
[364,] 0.9438182 1.123636e-01 5.618178e-02
[365,] 0.9322370 1.355261e-01 6.776304e-02
[366,] 0.9176026 1.647947e-01 8.239736e-02
[367,] 0.9043554 1.912892e-01 9.564460e-02
[368,] 0.8872373 2.255255e-01 1.127627e-01
[369,] 0.8678234 2.643532e-01 1.321766e-01
[370,] 0.8461841 3.076319e-01 1.538159e-01
[371,] 0.8189067 3.621866e-01 1.810933e-01
[372,] 0.7900137 4.199726e-01 2.099863e-01
[373,] 0.7615561 4.768878e-01 2.384439e-01
[374,] 0.7303780 5.392439e-01 2.696220e-01
[375,] 0.6972588 6.054823e-01 3.027412e-01
[376,] 0.6577735 6.844529e-01 3.422265e-01
[377,] 0.6208907 7.582186e-01 3.791093e-01
[378,] 0.5824354 8.351292e-01 4.175646e-01
[379,] 0.5369887 9.260225e-01 4.630113e-01
[380,] 0.4906021 9.812042e-01 5.093979e-01
[381,] 0.4478182 8.956365e-01 5.521818e-01
[382,] 0.4052253 8.104506e-01 5.947747e-01
[383,] 0.3637518 7.275035e-01 6.362482e-01
[384,] 0.3212151 6.424303e-01 6.787849e-01
[385,] 0.2808053 5.616107e-01 7.191947e-01
[386,] 0.2433549 4.867098e-01 7.566451e-01
[387,] 0.2080829 4.161658e-01 7.919171e-01
[388,] 0.1831176 3.662352e-01 8.168824e-01
[389,] 0.1847397 3.694794e-01 8.152603e-01
[390,] 0.1542966 3.085932e-01 8.457034e-01
[391,] 0.1263188 2.526376e-01 8.736812e-01
[392,] 0.9998227 3.546202e-04 1.773101e-04
[393,] 0.9996823 6.353814e-04 3.176907e-04
[394,] 0.9997480 5.039037e-04 2.519518e-04
[395,] 0.9996026 7.948014e-04 3.974007e-04
[396,] 0.9992983 1.403340e-03 7.016698e-04
[397,] 0.9987702 2.459577e-03 1.229788e-03
[398,] 0.9978372 4.325645e-03 2.162822e-03
[399,] 0.9963361 7.327835e-03 3.663918e-03
[400,] 0.9940502 1.189964e-02 5.949820e-03
[401,] 0.9902170 1.956593e-02 9.782967e-03
[402,] 0.9841478 3.170447e-02 1.585224e-02
[403,] 0.9748771 5.024573e-02 2.512287e-02
[404,] 0.9636439 7.271222e-02 3.635611e-02
[405,] 0.9738500 5.229995e-02 2.614998e-02
[406,] 0.9745168 5.096648e-02 2.548324e-02
[407,] 0.9585277 8.294470e-02 4.147235e-02
[408,] 0.9351713 1.296575e-01 6.482874e-02
[409,] 0.9050946 1.898109e-01 9.490545e-02
[410,] 0.8721892 2.556215e-01 1.278108e-01
[411,] 0.8136384 3.727232e-01 1.863616e-01
[412,] 0.7373297 5.253405e-01 2.626703e-01
[413,] 0.6444105 7.111789e-01 3.555895e-01
[414,] 0.5614008 8.771984e-01 4.385992e-01
[415,] 0.4995467 9.990934e-01 5.004533e-01
[416,] 0.5616004 8.767993e-01 4.383996e-01
[417,] 0.6038777 7.922447e-01 3.961223e-01
[418,] 0.4811527 9.623054e-01 5.188473e-01
> postscript(file="/var/www/html/rcomp/tmp/1plp11292936115.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/2plp11292936115.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/3plp11292936115.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/4hd6m1292936115.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/5hd6m1292936115.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
-130.62456449 207.79334272 -112.26455656 29.18643365 3.63482457
6 7 8 9 10
-94.81885653 18.57565718 -30.24135586 21.70176209 19.36167303
11 12 13 14 15
-86.59101631 38.24812170 286.32422167 -24.06338721 129.07075690
16 17 18 19 20
-15.18946340 0.19059469 125.09730397 32.86298850 162.88926644
21 22 23 24 25
-17.13496710 -64.60510293 -3.57357844 186.71670527 20.43441848
26 27 28 29 30
-38.62869905 80.88433341 11.20338030 -22.91166378 56.86844426
31 32 33 34 35
-15.87633601 -14.24201358 -18.79068589 -6.85565114 -69.29690936
36 37 38 39 40
98.98114436 -0.60087043 23.13146189 26.23152884 96.46713205
41 42 43 44 45
134.44917314 -30.94569690 43.71166625 28.72929285 0.64519469
46 47 48 49 50
-15.62559042 -20.47347369 81.63767597 5.53095440 15.99101616
51 52 53 54 55
47.65652634 83.45607185 18.39970423 22.29965975 59.75602632
56 57 58 59 60
50.71422965 -29.15787504 0.01050758 23.06755342 20.22741122
61 62 63 64 65
-7.22984443 14.95436832 -12.09344580 22.19916826 27.78031977
66 67 68 69 70
167.25442132 -5.48263787 11.12991749 37.01572608 42.13414125
71 72 73 74 75
30.67218546 3.88908437 5.36947546 -10.82257940 -2.38493941
76 77 78 79 80
-3.77907105 -1.39149175 -7.23762732 31.04721357 -15.85966015
81 82 83 84 85
8.71432865 22.19956634 -1.61548123 2.84655442 -10.61434171
86 87 88 89 90
26.16837306 60.52327270 -12.56087242 12.26736841 -12.16550127
91 92 93 94 95
4.49075950 6.92813611 40.98126276 -30.37710520 25.06648983
96 97 98 99 100
1.28995575 -18.62446858 -0.09282251 -22.78143049 -13.35236550
101 102 103 104 105
8.81145997 -3.52225122 15.97552526 -19.06330757 4.20756189
106 107 108 109 110
10.96884534 23.48641421 -18.44702831 -2.15469877 -3.63686258
111 112 113 114 115
1.19137508 9.20452636 1.42510254 -13.84374069 -14.17365641
116 117 118 119 120
15.18462428 -12.44437004 -15.10372818 -14.31341775 3.85762648
121 122 123 124 125
16.58075841 -6.90138075 83.71807617 -9.08693641 -9.87813040
126 127 128 129 130
-7.15741302 -13.01004748 -3.46000365 8.51128591 -13.01004748
131 132 133 134 135
-6.13429458 4.58352942 -12.01004748 6.79631498 -12.03816439
136 137 138 139 140
-10.47138135 -3.67593098 -12.02465886 -13.29954097 -13.01004748
141 142 143 144 145
-6.69696279 -13.01004748 -5.64564122 -9.22829329 -1.04779182
146 147 148 149 150
-4.55031080 -13.01004748 -6.07026144 -12.49579126 -13.14971888
151 152 153 154 155
-7.61018199 24.76624166 -8.01004748 -7.42370042 -11.14984283
156 157 158 159 160
-13.01004748 -13.01004748 -11.64122750 -10.37649821 -13.01004748
161 162 163 164 165
-13.01004748 -13.01004748 -13.01004748 -8.17523433 8.52238728
166 167 168 169 170
-12.22582453 -13.01004748 3.29305043 19.18697244 -13.01004748
171 172 173 174 175
-13.01004748 -10.30637057 -13.01004748 -5.13302816 -13.01004748
176 177 178 179 180
-5.66628636 98.11325251 -13.37070166 -13.01004748 -13.01004748
181 182 183 184 185
-1.83179678 -7.16339185 1.26927222 -0.86053356 -8.50093012
186 187 188 189 190
-5.33069669 -2.61405333 -7.96708434 -9.71788902 -13.01004748
191 192 193 194 195
-7.13130717 -13.01004748 -6.87961816 -13.01004748 -10.18633601
196 197 198 199 200
-13.01004748 39.07214005 -13.01004748 -13.01004748 0.47291583
201 202 203 204 205
-13.01004748 -13.01004748 -13.01004748 0.40646372 -9.26918373
206 207 208 209 210
-8.27238595 -12.43257522 -13.01004748 9.61379766 -11.12001887
211 212 213 214 215
-13.01004748 -13.01004748 -2.91061907 -5.04913646 -10.76776363
216 217 218 219 220
-11.01004748 -7.11145167 -7.94132213 -11.01004748 6.87557740
221 222 223 224 225
-10.99026770 -13.01004748 -12.01923585 4.80837403 -12.12773964
226 227 228 229 230
7.52053233 -4.28663498 -13.01004748 -12.84459258 -13.01004748
231 232 233 234 235
-7.89455090 -7.66017679 -13.01004748 -1.74921133 -13.01004748
236 237 238 239 240
-13.02656797 -13.01004748 -7.19268622 -13.01004748 -13.01004748
241 242 243 244 245
-13.01004748 -10.54977902 -13.01004748 -5.93745771 -24.90757444
246 247 248 249 250
-13.01004748 -13.01004748 1.25123907 -13.01004748 -13.01004748
251 252 253 254 255
-3.01918929 -13.01004748 -13.40966807 19.59623438 29.30913776
256 257 258 259 260
-13.01004748 -6.38979561 14.34355929 -13.01004748 -13.01004748
261 262 263 264 265
-13.01004748 -13.01004748 -4.21946369 -4.39771772 -13.01004748
266 267 268 269 270
-13.01004748 -13.01004748 -9.77221151 -13.01004748 10.08340022
271 272 273 274 275
-13.01004748 -13.01004748 -13.01004748 5.47252280 -13.01004748
276 277 278 279 280
-2.81884947 7.55228145 -13.01004748 8.65788379 -7.34654199
281 282 283 284 285
-13.01004748 -13.01004748 12.68482492 -13.01004748 -8.06128403
286 287 288 289 290
-13.01004748 -10.96261291 -10.48807205 -13.01004748 0.59078953
291 292 293 294 295
-13.01004748 -2.43904337 -1.49388020 -13.01004748 -2.75902125
296 297 298 299 300
-6.03662514 -13.01004748 -13.01004748 -13.01004748 -4.01004748
301 302 303 304 305
-9.12842915 -12.71831761 -12.01004748 -13.01004748 -1.54679346
306 307 308 309 310
-1.02309748 -13.01004748 4.26111624 3.54315760 -13.01004748
311 312 313 314 315
-13.01004748 11.00010465 -13.01004748 -0.09368542 -5.03840529
316 317 318 319 320
-11.22820302 -12.99930635 5.47699173 -8.01447432 -13.01004748
321 322 323 324 325
-13.01004748 -13.01004748 -13.01004748 -13.01004748 -13.01004748
326 327 328 329 330
-12.06132417 -13.01004748 -13.01004748 -13.01004748 -13.01004748
331 332 333 334 335
-13.01004748 -13.01004748 -13.01004748 -9.36600274 -12.01648475
336 337 338 339 340
-10.51453845 3.26667748 -13.01004748 -12.14374943 -3.44259767
341 342 343 344 345
-2.57123733 -13.01004748 -13.01004748 -13.01004748 -2.85074308
346 347 348 349 350
-13.01004748 -10.01653981 -13.01004748 -13.01004748 -10.72879238
351 352 353 354 355
-11.11509622 -6.35852982 -5.77851683 -12.92176715 -13.01004748
356 357 358 359 360
-13.01004748 -4.45943392 -13.01004748 -13.01004748 -13.01004748
361 362 363 364 365
-9.01560747 -13.01004748 -13.01004748 6.48759728 -4.66015046
366 367 368 369 370
-13.01004748 -13.01004748 -13.01004748 -13.01004748 -7.05796557
371 372 373 374 375
-13.01004748 -8.10779640 -10.89878890 -11.56783921 -13.01004748
376 377 378 379 380
-13.01004748 -1.24126971 -8.99011503 -13.69479069 -13.01004748
381 382 383 384 385
-13.01004748 -0.29870166 3.02316739 -13.01004748 -1.84328535
386 387 388 389 390
4.01618617 -10.63350187 -10.86463146 -10.26417164 -7.15312399
391 392 393 394 395
-5.46879276 -8.42799497 -7.05170825 -12.89307885 -24.71203709
396 397 398 399 400
-6.61715727 -7.52414074 203.84608517 -8.85458702 27.88477530
401 402 403 404 405
-16.95198399 -9.82122245 -9.26138167 -0.47417404 -8.87467144
406 407 408 409 410
-10.02281589 -6.57189259 -5.16203979 -5.26905751 -13.88120415
411 412 413 414 415
35.90812542 -27.56335233 -1.99717252 -2.48676517 -12.35691311
416 417 418 419 420
-16.25330300 -1.91803219 1.68941382 3.85325725 -14.98176598
421 422 423 424 425
-17.62138299 46.22141168 23.36785417 -16.45936332 -1.86649331
426 427 428 429 430
45.05644409 -3.98314145 39.85181744 -28.18161075 -4.67298230
431
-17.36361924
> postscript(file="/var/www/html/rcomp/tmp/6sm571292936115.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 -130.62456449 NA
1 207.79334272 -130.62456449
2 -112.26455656 207.79334272
3 29.18643365 -112.26455656
4 3.63482457 29.18643365
5 -94.81885653 3.63482457
6 18.57565718 -94.81885653
7 -30.24135586 18.57565718
8 21.70176209 -30.24135586
9 19.36167303 21.70176209
10 -86.59101631 19.36167303
11 38.24812170 -86.59101631
12 286.32422167 38.24812170
13 -24.06338721 286.32422167
14 129.07075690 -24.06338721
15 -15.18946340 129.07075690
16 0.19059469 -15.18946340
17 125.09730397 0.19059469
18 32.86298850 125.09730397
19 162.88926644 32.86298850
20 -17.13496710 162.88926644
21 -64.60510293 -17.13496710
22 -3.57357844 -64.60510293
23 186.71670527 -3.57357844
24 20.43441848 186.71670527
25 -38.62869905 20.43441848
26 80.88433341 -38.62869905
27 11.20338030 80.88433341
28 -22.91166378 11.20338030
29 56.86844426 -22.91166378
30 -15.87633601 56.86844426
31 -14.24201358 -15.87633601
32 -18.79068589 -14.24201358
33 -6.85565114 -18.79068589
34 -69.29690936 -6.85565114
35 98.98114436 -69.29690936
36 -0.60087043 98.98114436
37 23.13146189 -0.60087043
38 26.23152884 23.13146189
39 96.46713205 26.23152884
40 134.44917314 96.46713205
41 -30.94569690 134.44917314
42 43.71166625 -30.94569690
43 28.72929285 43.71166625
44 0.64519469 28.72929285
45 -15.62559042 0.64519469
46 -20.47347369 -15.62559042
47 81.63767597 -20.47347369
48 5.53095440 81.63767597
49 15.99101616 5.53095440
50 47.65652634 15.99101616
51 83.45607185 47.65652634
52 18.39970423 83.45607185
53 22.29965975 18.39970423
54 59.75602632 22.29965975
55 50.71422965 59.75602632
56 -29.15787504 50.71422965
57 0.01050758 -29.15787504
58 23.06755342 0.01050758
59 20.22741122 23.06755342
60 -7.22984443 20.22741122
61 14.95436832 -7.22984443
62 -12.09344580 14.95436832
63 22.19916826 -12.09344580
64 27.78031977 22.19916826
65 167.25442132 27.78031977
66 -5.48263787 167.25442132
67 11.12991749 -5.48263787
68 37.01572608 11.12991749
69 42.13414125 37.01572608
70 30.67218546 42.13414125
71 3.88908437 30.67218546
72 5.36947546 3.88908437
73 -10.82257940 5.36947546
74 -2.38493941 -10.82257940
75 -3.77907105 -2.38493941
76 -1.39149175 -3.77907105
77 -7.23762732 -1.39149175
78 31.04721357 -7.23762732
79 -15.85966015 31.04721357
80 8.71432865 -15.85966015
81 22.19956634 8.71432865
82 -1.61548123 22.19956634
83 2.84655442 -1.61548123
84 -10.61434171 2.84655442
85 26.16837306 -10.61434171
86 60.52327270 26.16837306
87 -12.56087242 60.52327270
88 12.26736841 -12.56087242
89 -12.16550127 12.26736841
90 4.49075950 -12.16550127
91 6.92813611 4.49075950
92 40.98126276 6.92813611
93 -30.37710520 40.98126276
94 25.06648983 -30.37710520
95 1.28995575 25.06648983
96 -18.62446858 1.28995575
97 -0.09282251 -18.62446858
98 -22.78143049 -0.09282251
99 -13.35236550 -22.78143049
100 8.81145997 -13.35236550
101 -3.52225122 8.81145997
102 15.97552526 -3.52225122
103 -19.06330757 15.97552526
104 4.20756189 -19.06330757
105 10.96884534 4.20756189
106 23.48641421 10.96884534
107 -18.44702831 23.48641421
108 -2.15469877 -18.44702831
109 -3.63686258 -2.15469877
110 1.19137508 -3.63686258
111 9.20452636 1.19137508
112 1.42510254 9.20452636
113 -13.84374069 1.42510254
114 -14.17365641 -13.84374069
115 15.18462428 -14.17365641
116 -12.44437004 15.18462428
117 -15.10372818 -12.44437004
118 -14.31341775 -15.10372818
119 3.85762648 -14.31341775
120 16.58075841 3.85762648
121 -6.90138075 16.58075841
122 83.71807617 -6.90138075
123 -9.08693641 83.71807617
124 -9.87813040 -9.08693641
125 -7.15741302 -9.87813040
126 -13.01004748 -7.15741302
127 -3.46000365 -13.01004748
128 8.51128591 -3.46000365
129 -13.01004748 8.51128591
130 -6.13429458 -13.01004748
131 4.58352942 -6.13429458
132 -12.01004748 4.58352942
133 6.79631498 -12.01004748
134 -12.03816439 6.79631498
135 -10.47138135 -12.03816439
136 -3.67593098 -10.47138135
137 -12.02465886 -3.67593098
138 -13.29954097 -12.02465886
139 -13.01004748 -13.29954097
140 -6.69696279 -13.01004748
141 -13.01004748 -6.69696279
142 -5.64564122 -13.01004748
143 -9.22829329 -5.64564122
144 -1.04779182 -9.22829329
145 -4.55031080 -1.04779182
146 -13.01004748 -4.55031080
147 -6.07026144 -13.01004748
148 -12.49579126 -6.07026144
149 -13.14971888 -12.49579126
150 -7.61018199 -13.14971888
151 24.76624166 -7.61018199
152 -8.01004748 24.76624166
153 -7.42370042 -8.01004748
154 -11.14984283 -7.42370042
155 -13.01004748 -11.14984283
156 -13.01004748 -13.01004748
157 -11.64122750 -13.01004748
158 -10.37649821 -11.64122750
159 -13.01004748 -10.37649821
160 -13.01004748 -13.01004748
161 -13.01004748 -13.01004748
162 -13.01004748 -13.01004748
163 -8.17523433 -13.01004748
164 8.52238728 -8.17523433
165 -12.22582453 8.52238728
166 -13.01004748 -12.22582453
167 3.29305043 -13.01004748
168 19.18697244 3.29305043
169 -13.01004748 19.18697244
170 -13.01004748 -13.01004748
171 -10.30637057 -13.01004748
172 -13.01004748 -10.30637057
173 -5.13302816 -13.01004748
174 -13.01004748 -5.13302816
175 -5.66628636 -13.01004748
176 98.11325251 -5.66628636
177 -13.37070166 98.11325251
178 -13.01004748 -13.37070166
179 -13.01004748 -13.01004748
180 -1.83179678 -13.01004748
181 -7.16339185 -1.83179678
182 1.26927222 -7.16339185
183 -0.86053356 1.26927222
184 -8.50093012 -0.86053356
185 -5.33069669 -8.50093012
186 -2.61405333 -5.33069669
187 -7.96708434 -2.61405333
188 -9.71788902 -7.96708434
189 -13.01004748 -9.71788902
190 -7.13130717 -13.01004748
191 -13.01004748 -7.13130717
192 -6.87961816 -13.01004748
193 -13.01004748 -6.87961816
194 -10.18633601 -13.01004748
195 -13.01004748 -10.18633601
196 39.07214005 -13.01004748
197 -13.01004748 39.07214005
198 -13.01004748 -13.01004748
199 0.47291583 -13.01004748
200 -13.01004748 0.47291583
201 -13.01004748 -13.01004748
202 -13.01004748 -13.01004748
203 0.40646372 -13.01004748
204 -9.26918373 0.40646372
205 -8.27238595 -9.26918373
206 -12.43257522 -8.27238595
207 -13.01004748 -12.43257522
208 9.61379766 -13.01004748
209 -11.12001887 9.61379766
210 -13.01004748 -11.12001887
211 -13.01004748 -13.01004748
212 -2.91061907 -13.01004748
213 -5.04913646 -2.91061907
214 -10.76776363 -5.04913646
215 -11.01004748 -10.76776363
216 -7.11145167 -11.01004748
217 -7.94132213 -7.11145167
218 -11.01004748 -7.94132213
219 6.87557740 -11.01004748
220 -10.99026770 6.87557740
221 -13.01004748 -10.99026770
222 -12.01923585 -13.01004748
223 4.80837403 -12.01923585
224 -12.12773964 4.80837403
225 7.52053233 -12.12773964
226 -4.28663498 7.52053233
227 -13.01004748 -4.28663498
228 -12.84459258 -13.01004748
229 -13.01004748 -12.84459258
230 -7.89455090 -13.01004748
231 -7.66017679 -7.89455090
232 -13.01004748 -7.66017679
233 -1.74921133 -13.01004748
234 -13.01004748 -1.74921133
235 -13.02656797 -13.01004748
236 -13.01004748 -13.02656797
237 -7.19268622 -13.01004748
238 -13.01004748 -7.19268622
239 -13.01004748 -13.01004748
240 -13.01004748 -13.01004748
241 -10.54977902 -13.01004748
242 -13.01004748 -10.54977902
243 -5.93745771 -13.01004748
244 -24.90757444 -5.93745771
245 -13.01004748 -24.90757444
246 -13.01004748 -13.01004748
247 1.25123907 -13.01004748
248 -13.01004748 1.25123907
249 -13.01004748 -13.01004748
250 -3.01918929 -13.01004748
251 -13.01004748 -3.01918929
252 -13.40966807 -13.01004748
253 19.59623438 -13.40966807
254 29.30913776 19.59623438
255 -13.01004748 29.30913776
256 -6.38979561 -13.01004748
257 14.34355929 -6.38979561
258 -13.01004748 14.34355929
259 -13.01004748 -13.01004748
260 -13.01004748 -13.01004748
261 -13.01004748 -13.01004748
262 -4.21946369 -13.01004748
263 -4.39771772 -4.21946369
264 -13.01004748 -4.39771772
265 -13.01004748 -13.01004748
266 -13.01004748 -13.01004748
267 -9.77221151 -13.01004748
268 -13.01004748 -9.77221151
269 10.08340022 -13.01004748
270 -13.01004748 10.08340022
271 -13.01004748 -13.01004748
272 -13.01004748 -13.01004748
273 5.47252280 -13.01004748
274 -13.01004748 5.47252280
275 -2.81884947 -13.01004748
276 7.55228145 -2.81884947
277 -13.01004748 7.55228145
278 8.65788379 -13.01004748
279 -7.34654199 8.65788379
280 -13.01004748 -7.34654199
281 -13.01004748 -13.01004748
282 12.68482492 -13.01004748
283 -13.01004748 12.68482492
284 -8.06128403 -13.01004748
285 -13.01004748 -8.06128403
286 -10.96261291 -13.01004748
287 -10.48807205 -10.96261291
288 -13.01004748 -10.48807205
289 0.59078953 -13.01004748
290 -13.01004748 0.59078953
291 -2.43904337 -13.01004748
292 -1.49388020 -2.43904337
293 -13.01004748 -1.49388020
294 -2.75902125 -13.01004748
295 -6.03662514 -2.75902125
296 -13.01004748 -6.03662514
297 -13.01004748 -13.01004748
298 -13.01004748 -13.01004748
299 -4.01004748 -13.01004748
300 -9.12842915 -4.01004748
301 -12.71831761 -9.12842915
302 -12.01004748 -12.71831761
303 -13.01004748 -12.01004748
304 -1.54679346 -13.01004748
305 -1.02309748 -1.54679346
306 -13.01004748 -1.02309748
307 4.26111624 -13.01004748
308 3.54315760 4.26111624
309 -13.01004748 3.54315760
310 -13.01004748 -13.01004748
311 11.00010465 -13.01004748
312 -13.01004748 11.00010465
313 -0.09368542 -13.01004748
314 -5.03840529 -0.09368542
315 -11.22820302 -5.03840529
316 -12.99930635 -11.22820302
317 5.47699173 -12.99930635
318 -8.01447432 5.47699173
319 -13.01004748 -8.01447432
320 -13.01004748 -13.01004748
321 -13.01004748 -13.01004748
322 -13.01004748 -13.01004748
323 -13.01004748 -13.01004748
324 -13.01004748 -13.01004748
325 -12.06132417 -13.01004748
326 -13.01004748 -12.06132417
327 -13.01004748 -13.01004748
328 -13.01004748 -13.01004748
329 -13.01004748 -13.01004748
330 -13.01004748 -13.01004748
331 -13.01004748 -13.01004748
332 -13.01004748 -13.01004748
333 -9.36600274 -13.01004748
334 -12.01648475 -9.36600274
335 -10.51453845 -12.01648475
336 3.26667748 -10.51453845
337 -13.01004748 3.26667748
338 -12.14374943 -13.01004748
339 -3.44259767 -12.14374943
340 -2.57123733 -3.44259767
341 -13.01004748 -2.57123733
342 -13.01004748 -13.01004748
343 -13.01004748 -13.01004748
344 -2.85074308 -13.01004748
345 -13.01004748 -2.85074308
346 -10.01653981 -13.01004748
347 -13.01004748 -10.01653981
348 -13.01004748 -13.01004748
349 -10.72879238 -13.01004748
350 -11.11509622 -10.72879238
351 -6.35852982 -11.11509622
352 -5.77851683 -6.35852982
353 -12.92176715 -5.77851683
354 -13.01004748 -12.92176715
355 -13.01004748 -13.01004748
356 -4.45943392 -13.01004748
357 -13.01004748 -4.45943392
358 -13.01004748 -13.01004748
359 -13.01004748 -13.01004748
360 -9.01560747 -13.01004748
361 -13.01004748 -9.01560747
362 -13.01004748 -13.01004748
363 6.48759728 -13.01004748
364 -4.66015046 6.48759728
365 -13.01004748 -4.66015046
366 -13.01004748 -13.01004748
367 -13.01004748 -13.01004748
368 -13.01004748 -13.01004748
369 -7.05796557 -13.01004748
370 -13.01004748 -7.05796557
371 -8.10779640 -13.01004748
372 -10.89878890 -8.10779640
373 -11.56783921 -10.89878890
374 -13.01004748 -11.56783921
375 -13.01004748 -13.01004748
376 -1.24126971 -13.01004748
377 -8.99011503 -1.24126971
378 -13.69479069 -8.99011503
379 -13.01004748 -13.69479069
380 -13.01004748 -13.01004748
381 -0.29870166 -13.01004748
382 3.02316739 -0.29870166
383 -13.01004748 3.02316739
384 -1.84328535 -13.01004748
385 4.01618617 -1.84328535
386 -10.63350187 4.01618617
387 -10.86463146 -10.63350187
388 -10.26417164 -10.86463146
389 -7.15312399 -10.26417164
390 -5.46879276 -7.15312399
391 -8.42799497 -5.46879276
392 -7.05170825 -8.42799497
393 -12.89307885 -7.05170825
394 -24.71203709 -12.89307885
395 -6.61715727 -24.71203709
396 -7.52414074 -6.61715727
397 203.84608517 -7.52414074
398 -8.85458702 203.84608517
399 27.88477530 -8.85458702
400 -16.95198399 27.88477530
401 -9.82122245 -16.95198399
402 -9.26138167 -9.82122245
403 -0.47417404 -9.26138167
404 -8.87467144 -0.47417404
405 -10.02281589 -8.87467144
406 -6.57189259 -10.02281589
407 -5.16203979 -6.57189259
408 -5.26905751 -5.16203979
409 -13.88120415 -5.26905751
410 35.90812542 -13.88120415
411 -27.56335233 35.90812542
412 -1.99717252 -27.56335233
413 -2.48676517 -1.99717252
414 -12.35691311 -2.48676517
415 -16.25330300 -12.35691311
416 -1.91803219 -16.25330300
417 1.68941382 -1.91803219
418 3.85325725 1.68941382
419 -14.98176598 3.85325725
420 -17.62138299 -14.98176598
421 46.22141168 -17.62138299
422 23.36785417 46.22141168
423 -16.45936332 23.36785417
424 -1.86649331 -16.45936332
425 45.05644409 -1.86649331
426 -3.98314145 45.05644409
427 39.85181744 -3.98314145
428 -28.18161075 39.85181744
429 -4.67298230 -28.18161075
430 -17.36361924 -4.67298230
431 NA -17.36361924
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 207.79334272 -130.62456449
[2,] -112.26455656 207.79334272
[3,] 29.18643365 -112.26455656
[4,] 3.63482457 29.18643365
[5,] -94.81885653 3.63482457
[6,] 18.57565718 -94.81885653
[7,] -30.24135586 18.57565718
[8,] 21.70176209 -30.24135586
[9,] 19.36167303 21.70176209
[10,] -86.59101631 19.36167303
[11,] 38.24812170 -86.59101631
[12,] 286.32422167 38.24812170
[13,] -24.06338721 286.32422167
[14,] 129.07075690 -24.06338721
[15,] -15.18946340 129.07075690
[16,] 0.19059469 -15.18946340
[17,] 125.09730397 0.19059469
[18,] 32.86298850 125.09730397
[19,] 162.88926644 32.86298850
[20,] -17.13496710 162.88926644
[21,] -64.60510293 -17.13496710
[22,] -3.57357844 -64.60510293
[23,] 186.71670527 -3.57357844
[24,] 20.43441848 186.71670527
[25,] -38.62869905 20.43441848
[26,] 80.88433341 -38.62869905
[27,] 11.20338030 80.88433341
[28,] -22.91166378 11.20338030
[29,] 56.86844426 -22.91166378
[30,] -15.87633601 56.86844426
[31,] -14.24201358 -15.87633601
[32,] -18.79068589 -14.24201358
[33,] -6.85565114 -18.79068589
[34,] -69.29690936 -6.85565114
[35,] 98.98114436 -69.29690936
[36,] -0.60087043 98.98114436
[37,] 23.13146189 -0.60087043
[38,] 26.23152884 23.13146189
[39,] 96.46713205 26.23152884
[40,] 134.44917314 96.46713205
[41,] -30.94569690 134.44917314
[42,] 43.71166625 -30.94569690
[43,] 28.72929285 43.71166625
[44,] 0.64519469 28.72929285
[45,] -15.62559042 0.64519469
[46,] -20.47347369 -15.62559042
[47,] 81.63767597 -20.47347369
[48,] 5.53095440 81.63767597
[49,] 15.99101616 5.53095440
[50,] 47.65652634 15.99101616
[51,] 83.45607185 47.65652634
[52,] 18.39970423 83.45607185
[53,] 22.29965975 18.39970423
[54,] 59.75602632 22.29965975
[55,] 50.71422965 59.75602632
[56,] -29.15787504 50.71422965
[57,] 0.01050758 -29.15787504
[58,] 23.06755342 0.01050758
[59,] 20.22741122 23.06755342
[60,] -7.22984443 20.22741122
[61,] 14.95436832 -7.22984443
[62,] -12.09344580 14.95436832
[63,] 22.19916826 -12.09344580
[64,] 27.78031977 22.19916826
[65,] 167.25442132 27.78031977
[66,] -5.48263787 167.25442132
[67,] 11.12991749 -5.48263787
[68,] 37.01572608 11.12991749
[69,] 42.13414125 37.01572608
[70,] 30.67218546 42.13414125
[71,] 3.88908437 30.67218546
[72,] 5.36947546 3.88908437
[73,] -10.82257940 5.36947546
[74,] -2.38493941 -10.82257940
[75,] -3.77907105 -2.38493941
[76,] -1.39149175 -3.77907105
[77,] -7.23762732 -1.39149175
[78,] 31.04721357 -7.23762732
[79,] -15.85966015 31.04721357
[80,] 8.71432865 -15.85966015
[81,] 22.19956634 8.71432865
[82,] -1.61548123 22.19956634
[83,] 2.84655442 -1.61548123
[84,] -10.61434171 2.84655442
[85,] 26.16837306 -10.61434171
[86,] 60.52327270 26.16837306
[87,] -12.56087242 60.52327270
[88,] 12.26736841 -12.56087242
[89,] -12.16550127 12.26736841
[90,] 4.49075950 -12.16550127
[91,] 6.92813611 4.49075950
[92,] 40.98126276 6.92813611
[93,] -30.37710520 40.98126276
[94,] 25.06648983 -30.37710520
[95,] 1.28995575 25.06648983
[96,] -18.62446858 1.28995575
[97,] -0.09282251 -18.62446858
[98,] -22.78143049 -0.09282251
[99,] -13.35236550 -22.78143049
[100,] 8.81145997 -13.35236550
[101,] -3.52225122 8.81145997
[102,] 15.97552526 -3.52225122
[103,] -19.06330757 15.97552526
[104,] 4.20756189 -19.06330757
[105,] 10.96884534 4.20756189
[106,] 23.48641421 10.96884534
[107,] -18.44702831 23.48641421
[108,] -2.15469877 -18.44702831
[109,] -3.63686258 -2.15469877
[110,] 1.19137508 -3.63686258
[111,] 9.20452636 1.19137508
[112,] 1.42510254 9.20452636
[113,] -13.84374069 1.42510254
[114,] -14.17365641 -13.84374069
[115,] 15.18462428 -14.17365641
[116,] -12.44437004 15.18462428
[117,] -15.10372818 -12.44437004
[118,] -14.31341775 -15.10372818
[119,] 3.85762648 -14.31341775
[120,] 16.58075841 3.85762648
[121,] -6.90138075 16.58075841
[122,] 83.71807617 -6.90138075
[123,] -9.08693641 83.71807617
[124,] -9.87813040 -9.08693641
[125,] -7.15741302 -9.87813040
[126,] -13.01004748 -7.15741302
[127,] -3.46000365 -13.01004748
[128,] 8.51128591 -3.46000365
[129,] -13.01004748 8.51128591
[130,] -6.13429458 -13.01004748
[131,] 4.58352942 -6.13429458
[132,] -12.01004748 4.58352942
[133,] 6.79631498 -12.01004748
[134,] -12.03816439 6.79631498
[135,] -10.47138135 -12.03816439
[136,] -3.67593098 -10.47138135
[137,] -12.02465886 -3.67593098
[138,] -13.29954097 -12.02465886
[139,] -13.01004748 -13.29954097
[140,] -6.69696279 -13.01004748
[141,] -13.01004748 -6.69696279
[142,] -5.64564122 -13.01004748
[143,] -9.22829329 -5.64564122
[144,] -1.04779182 -9.22829329
[145,] -4.55031080 -1.04779182
[146,] -13.01004748 -4.55031080
[147,] -6.07026144 -13.01004748
[148,] -12.49579126 -6.07026144
[149,] -13.14971888 -12.49579126
[150,] -7.61018199 -13.14971888
[151,] 24.76624166 -7.61018199
[152,] -8.01004748 24.76624166
[153,] -7.42370042 -8.01004748
[154,] -11.14984283 -7.42370042
[155,] -13.01004748 -11.14984283
[156,] -13.01004748 -13.01004748
[157,] -11.64122750 -13.01004748
[158,] -10.37649821 -11.64122750
[159,] -13.01004748 -10.37649821
[160,] -13.01004748 -13.01004748
[161,] -13.01004748 -13.01004748
[162,] -13.01004748 -13.01004748
[163,] -8.17523433 -13.01004748
[164,] 8.52238728 -8.17523433
[165,] -12.22582453 8.52238728
[166,] -13.01004748 -12.22582453
[167,] 3.29305043 -13.01004748
[168,] 19.18697244 3.29305043
[169,] -13.01004748 19.18697244
[170,] -13.01004748 -13.01004748
[171,] -10.30637057 -13.01004748
[172,] -13.01004748 -10.30637057
[173,] -5.13302816 -13.01004748
[174,] -13.01004748 -5.13302816
[175,] -5.66628636 -13.01004748
[176,] 98.11325251 -5.66628636
[177,] -13.37070166 98.11325251
[178,] -13.01004748 -13.37070166
[179,] -13.01004748 -13.01004748
[180,] -1.83179678 -13.01004748
[181,] -7.16339185 -1.83179678
[182,] 1.26927222 -7.16339185
[183,] -0.86053356 1.26927222
[184,] -8.50093012 -0.86053356
[185,] -5.33069669 -8.50093012
[186,] -2.61405333 -5.33069669
[187,] -7.96708434 -2.61405333
[188,] -9.71788902 -7.96708434
[189,] -13.01004748 -9.71788902
[190,] -7.13130717 -13.01004748
[191,] -13.01004748 -7.13130717
[192,] -6.87961816 -13.01004748
[193,] -13.01004748 -6.87961816
[194,] -10.18633601 -13.01004748
[195,] -13.01004748 -10.18633601
[196,] 39.07214005 -13.01004748
[197,] -13.01004748 39.07214005
[198,] -13.01004748 -13.01004748
[199,] 0.47291583 -13.01004748
[200,] -13.01004748 0.47291583
[201,] -13.01004748 -13.01004748
[202,] -13.01004748 -13.01004748
[203,] 0.40646372 -13.01004748
[204,] -9.26918373 0.40646372
[205,] -8.27238595 -9.26918373
[206,] -12.43257522 -8.27238595
[207,] -13.01004748 -12.43257522
[208,] 9.61379766 -13.01004748
[209,] -11.12001887 9.61379766
[210,] -13.01004748 -11.12001887
[211,] -13.01004748 -13.01004748
[212,] -2.91061907 -13.01004748
[213,] -5.04913646 -2.91061907
[214,] -10.76776363 -5.04913646
[215,] -11.01004748 -10.76776363
[216,] -7.11145167 -11.01004748
[217,] -7.94132213 -7.11145167
[218,] -11.01004748 -7.94132213
[219,] 6.87557740 -11.01004748
[220,] -10.99026770 6.87557740
[221,] -13.01004748 -10.99026770
[222,] -12.01923585 -13.01004748
[223,] 4.80837403 -12.01923585
[224,] -12.12773964 4.80837403
[225,] 7.52053233 -12.12773964
[226,] -4.28663498 7.52053233
[227,] -13.01004748 -4.28663498
[228,] -12.84459258 -13.01004748
[229,] -13.01004748 -12.84459258
[230,] -7.89455090 -13.01004748
[231,] -7.66017679 -7.89455090
[232,] -13.01004748 -7.66017679
[233,] -1.74921133 -13.01004748
[234,] -13.01004748 -1.74921133
[235,] -13.02656797 -13.01004748
[236,] -13.01004748 -13.02656797
[237,] -7.19268622 -13.01004748
[238,] -13.01004748 -7.19268622
[239,] -13.01004748 -13.01004748
[240,] -13.01004748 -13.01004748
[241,] -10.54977902 -13.01004748
[242,] -13.01004748 -10.54977902
[243,] -5.93745771 -13.01004748
[244,] -24.90757444 -5.93745771
[245,] -13.01004748 -24.90757444
[246,] -13.01004748 -13.01004748
[247,] 1.25123907 -13.01004748
[248,] -13.01004748 1.25123907
[249,] -13.01004748 -13.01004748
[250,] -3.01918929 -13.01004748
[251,] -13.01004748 -3.01918929
[252,] -13.40966807 -13.01004748
[253,] 19.59623438 -13.40966807
[254,] 29.30913776 19.59623438
[255,] -13.01004748 29.30913776
[256,] -6.38979561 -13.01004748
[257,] 14.34355929 -6.38979561
[258,] -13.01004748 14.34355929
[259,] -13.01004748 -13.01004748
[260,] -13.01004748 -13.01004748
[261,] -13.01004748 -13.01004748
[262,] -4.21946369 -13.01004748
[263,] -4.39771772 -4.21946369
[264,] -13.01004748 -4.39771772
[265,] -13.01004748 -13.01004748
[266,] -13.01004748 -13.01004748
[267,] -9.77221151 -13.01004748
[268,] -13.01004748 -9.77221151
[269,] 10.08340022 -13.01004748
[270,] -13.01004748 10.08340022
[271,] -13.01004748 -13.01004748
[272,] -13.01004748 -13.01004748
[273,] 5.47252280 -13.01004748
[274,] -13.01004748 5.47252280
[275,] -2.81884947 -13.01004748
[276,] 7.55228145 -2.81884947
[277,] -13.01004748 7.55228145
[278,] 8.65788379 -13.01004748
[279,] -7.34654199 8.65788379
[280,] -13.01004748 -7.34654199
[281,] -13.01004748 -13.01004748
[282,] 12.68482492 -13.01004748
[283,] -13.01004748 12.68482492
[284,] -8.06128403 -13.01004748
[285,] -13.01004748 -8.06128403
[286,] -10.96261291 -13.01004748
[287,] -10.48807205 -10.96261291
[288,] -13.01004748 -10.48807205
[289,] 0.59078953 -13.01004748
[290,] -13.01004748 0.59078953
[291,] -2.43904337 -13.01004748
[292,] -1.49388020 -2.43904337
[293,] -13.01004748 -1.49388020
[294,] -2.75902125 -13.01004748
[295,] -6.03662514 -2.75902125
[296,] -13.01004748 -6.03662514
[297,] -13.01004748 -13.01004748
[298,] -13.01004748 -13.01004748
[299,] -4.01004748 -13.01004748
[300,] -9.12842915 -4.01004748
[301,] -12.71831761 -9.12842915
[302,] -12.01004748 -12.71831761
[303,] -13.01004748 -12.01004748
[304,] -1.54679346 -13.01004748
[305,] -1.02309748 -1.54679346
[306,] -13.01004748 -1.02309748
[307,] 4.26111624 -13.01004748
[308,] 3.54315760 4.26111624
[309,] -13.01004748 3.54315760
[310,] -13.01004748 -13.01004748
[311,] 11.00010465 -13.01004748
[312,] -13.01004748 11.00010465
[313,] -0.09368542 -13.01004748
[314,] -5.03840529 -0.09368542
[315,] -11.22820302 -5.03840529
[316,] -12.99930635 -11.22820302
[317,] 5.47699173 -12.99930635
[318,] -8.01447432 5.47699173
[319,] -13.01004748 -8.01447432
[320,] -13.01004748 -13.01004748
[321,] -13.01004748 -13.01004748
[322,] -13.01004748 -13.01004748
[323,] -13.01004748 -13.01004748
[324,] -13.01004748 -13.01004748
[325,] -12.06132417 -13.01004748
[326,] -13.01004748 -12.06132417
[327,] -13.01004748 -13.01004748
[328,] -13.01004748 -13.01004748
[329,] -13.01004748 -13.01004748
[330,] -13.01004748 -13.01004748
[331,] -13.01004748 -13.01004748
[332,] -13.01004748 -13.01004748
[333,] -9.36600274 -13.01004748
[334,] -12.01648475 -9.36600274
[335,] -10.51453845 -12.01648475
[336,] 3.26667748 -10.51453845
[337,] -13.01004748 3.26667748
[338,] -12.14374943 -13.01004748
[339,] -3.44259767 -12.14374943
[340,] -2.57123733 -3.44259767
[341,] -13.01004748 -2.57123733
[342,] -13.01004748 -13.01004748
[343,] -13.01004748 -13.01004748
[344,] -2.85074308 -13.01004748
[345,] -13.01004748 -2.85074308
[346,] -10.01653981 -13.01004748
[347,] -13.01004748 -10.01653981
[348,] -13.01004748 -13.01004748
[349,] -10.72879238 -13.01004748
[350,] -11.11509622 -10.72879238
[351,] -6.35852982 -11.11509622
[352,] -5.77851683 -6.35852982
[353,] -12.92176715 -5.77851683
[354,] -13.01004748 -12.92176715
[355,] -13.01004748 -13.01004748
[356,] -4.45943392 -13.01004748
[357,] -13.01004748 -4.45943392
[358,] -13.01004748 -13.01004748
[359,] -13.01004748 -13.01004748
[360,] -9.01560747 -13.01004748
[361,] -13.01004748 -9.01560747
[362,] -13.01004748 -13.01004748
[363,] 6.48759728 -13.01004748
[364,] -4.66015046 6.48759728
[365,] -13.01004748 -4.66015046
[366,] -13.01004748 -13.01004748
[367,] -13.01004748 -13.01004748
[368,] -13.01004748 -13.01004748
[369,] -7.05796557 -13.01004748
[370,] -13.01004748 -7.05796557
[371,] -8.10779640 -13.01004748
[372,] -10.89878890 -8.10779640
[373,] -11.56783921 -10.89878890
[374,] -13.01004748 -11.56783921
[375,] -13.01004748 -13.01004748
[376,] -1.24126971 -13.01004748
[377,] -8.99011503 -1.24126971
[378,] -13.69479069 -8.99011503
[379,] -13.01004748 -13.69479069
[380,] -13.01004748 -13.01004748
[381,] -0.29870166 -13.01004748
[382,] 3.02316739 -0.29870166
[383,] -13.01004748 3.02316739
[384,] -1.84328535 -13.01004748
[385,] 4.01618617 -1.84328535
[386,] -10.63350187 4.01618617
[387,] -10.86463146 -10.63350187
[388,] -10.26417164 -10.86463146
[389,] -7.15312399 -10.26417164
[390,] -5.46879276 -7.15312399
[391,] -8.42799497 -5.46879276
[392,] -7.05170825 -8.42799497
[393,] -12.89307885 -7.05170825
[394,] -24.71203709 -12.89307885
[395,] -6.61715727 -24.71203709
[396,] -7.52414074 -6.61715727
[397,] 203.84608517 -7.52414074
[398,] -8.85458702 203.84608517
[399,] 27.88477530 -8.85458702
[400,] -16.95198399 27.88477530
[401,] -9.82122245 -16.95198399
[402,] -9.26138167 -9.82122245
[403,] -0.47417404 -9.26138167
[404,] -8.87467144 -0.47417404
[405,] -10.02281589 -8.87467144
[406,] -6.57189259 -10.02281589
[407,] -5.16203979 -6.57189259
[408,] -5.26905751 -5.16203979
[409,] -13.88120415 -5.26905751
[410,] 35.90812542 -13.88120415
[411,] -27.56335233 35.90812542
[412,] -1.99717252 -27.56335233
[413,] -2.48676517 -1.99717252
[414,] -12.35691311 -2.48676517
[415,] -16.25330300 -12.35691311
[416,] -1.91803219 -16.25330300
[417,] 1.68941382 -1.91803219
[418,] 3.85325725 1.68941382
[419,] -14.98176598 3.85325725
[420,] -17.62138299 -14.98176598
[421,] 46.22141168 -17.62138299
[422,] 23.36785417 46.22141168
[423,] -16.45936332 23.36785417
[424,] -1.86649331 -16.45936332
[425,] 45.05644409 -1.86649331
[426,] -3.98314145 45.05644409
[427,] 39.85181744 -3.98314145
[428,] -28.18161075 39.85181744
[429,] -4.67298230 -28.18161075
[430,] -17.36361924 -4.67298230
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 207.79334272 -130.62456449
2 -112.26455656 207.79334272
3 29.18643365 -112.26455656
4 3.63482457 29.18643365
5 -94.81885653 3.63482457
6 18.57565718 -94.81885653
7 -30.24135586 18.57565718
8 21.70176209 -30.24135586
9 19.36167303 21.70176209
10 -86.59101631 19.36167303
11 38.24812170 -86.59101631
12 286.32422167 38.24812170
13 -24.06338721 286.32422167
14 129.07075690 -24.06338721
15 -15.18946340 129.07075690
16 0.19059469 -15.18946340
17 125.09730397 0.19059469
18 32.86298850 125.09730397
19 162.88926644 32.86298850
20 -17.13496710 162.88926644
21 -64.60510293 -17.13496710
22 -3.57357844 -64.60510293
23 186.71670527 -3.57357844
24 20.43441848 186.71670527
25 -38.62869905 20.43441848
26 80.88433341 -38.62869905
27 11.20338030 80.88433341
28 -22.91166378 11.20338030
29 56.86844426 -22.91166378
30 -15.87633601 56.86844426
31 -14.24201358 -15.87633601
32 -18.79068589 -14.24201358
33 -6.85565114 -18.79068589
34 -69.29690936 -6.85565114
35 98.98114436 -69.29690936
36 -0.60087043 98.98114436
37 23.13146189 -0.60087043
38 26.23152884 23.13146189
39 96.46713205 26.23152884
40 134.44917314 96.46713205
41 -30.94569690 134.44917314
42 43.71166625 -30.94569690
43 28.72929285 43.71166625
44 0.64519469 28.72929285
45 -15.62559042 0.64519469
46 -20.47347369 -15.62559042
47 81.63767597 -20.47347369
48 5.53095440 81.63767597
49 15.99101616 5.53095440
50 47.65652634 15.99101616
51 83.45607185 47.65652634
52 18.39970423 83.45607185
53 22.29965975 18.39970423
54 59.75602632 22.29965975
55 50.71422965 59.75602632
56 -29.15787504 50.71422965
57 0.01050758 -29.15787504
58 23.06755342 0.01050758
59 20.22741122 23.06755342
60 -7.22984443 20.22741122
61 14.95436832 -7.22984443
62 -12.09344580 14.95436832
63 22.19916826 -12.09344580
64 27.78031977 22.19916826
65 167.25442132 27.78031977
66 -5.48263787 167.25442132
67 11.12991749 -5.48263787
68 37.01572608 11.12991749
69 42.13414125 37.01572608
70 30.67218546 42.13414125
71 3.88908437 30.67218546
72 5.36947546 3.88908437
73 -10.82257940 5.36947546
74 -2.38493941 -10.82257940
75 -3.77907105 -2.38493941
76 -1.39149175 -3.77907105
77 -7.23762732 -1.39149175
78 31.04721357 -7.23762732
79 -15.85966015 31.04721357
80 8.71432865 -15.85966015
81 22.19956634 8.71432865
82 -1.61548123 22.19956634
83 2.84655442 -1.61548123
84 -10.61434171 2.84655442
85 26.16837306 -10.61434171
86 60.52327270 26.16837306
87 -12.56087242 60.52327270
88 12.26736841 -12.56087242
89 -12.16550127 12.26736841
90 4.49075950 -12.16550127
91 6.92813611 4.49075950
92 40.98126276 6.92813611
93 -30.37710520 40.98126276
94 25.06648983 -30.37710520
95 1.28995575 25.06648983
96 -18.62446858 1.28995575
97 -0.09282251 -18.62446858
98 -22.78143049 -0.09282251
99 -13.35236550 -22.78143049
100 8.81145997 -13.35236550
101 -3.52225122 8.81145997
102 15.97552526 -3.52225122
103 -19.06330757 15.97552526
104 4.20756189 -19.06330757
105 10.96884534 4.20756189
106 23.48641421 10.96884534
107 -18.44702831 23.48641421
108 -2.15469877 -18.44702831
109 -3.63686258 -2.15469877
110 1.19137508 -3.63686258
111 9.20452636 1.19137508
112 1.42510254 9.20452636
113 -13.84374069 1.42510254
114 -14.17365641 -13.84374069
115 15.18462428 -14.17365641
116 -12.44437004 15.18462428
117 -15.10372818 -12.44437004
118 -14.31341775 -15.10372818
119 3.85762648 -14.31341775
120 16.58075841 3.85762648
121 -6.90138075 16.58075841
122 83.71807617 -6.90138075
123 -9.08693641 83.71807617
124 -9.87813040 -9.08693641
125 -7.15741302 -9.87813040
126 -13.01004748 -7.15741302
127 -3.46000365 -13.01004748
128 8.51128591 -3.46000365
129 -13.01004748 8.51128591
130 -6.13429458 -13.01004748
131 4.58352942 -6.13429458
132 -12.01004748 4.58352942
133 6.79631498 -12.01004748
134 -12.03816439 6.79631498
135 -10.47138135 -12.03816439
136 -3.67593098 -10.47138135
137 -12.02465886 -3.67593098
138 -13.29954097 -12.02465886
139 -13.01004748 -13.29954097
140 -6.69696279 -13.01004748
141 -13.01004748 -6.69696279
142 -5.64564122 -13.01004748
143 -9.22829329 -5.64564122
144 -1.04779182 -9.22829329
145 -4.55031080 -1.04779182
146 -13.01004748 -4.55031080
147 -6.07026144 -13.01004748
148 -12.49579126 -6.07026144
149 -13.14971888 -12.49579126
150 -7.61018199 -13.14971888
151 24.76624166 -7.61018199
152 -8.01004748 24.76624166
153 -7.42370042 -8.01004748
154 -11.14984283 -7.42370042
155 -13.01004748 -11.14984283
156 -13.01004748 -13.01004748
157 -11.64122750 -13.01004748
158 -10.37649821 -11.64122750
159 -13.01004748 -10.37649821
160 -13.01004748 -13.01004748
161 -13.01004748 -13.01004748
162 -13.01004748 -13.01004748
163 -8.17523433 -13.01004748
164 8.52238728 -8.17523433
165 -12.22582453 8.52238728
166 -13.01004748 -12.22582453
167 3.29305043 -13.01004748
168 19.18697244 3.29305043
169 -13.01004748 19.18697244
170 -13.01004748 -13.01004748
171 -10.30637057 -13.01004748
172 -13.01004748 -10.30637057
173 -5.13302816 -13.01004748
174 -13.01004748 -5.13302816
175 -5.66628636 -13.01004748
176 98.11325251 -5.66628636
177 -13.37070166 98.11325251
178 -13.01004748 -13.37070166
179 -13.01004748 -13.01004748
180 -1.83179678 -13.01004748
181 -7.16339185 -1.83179678
182 1.26927222 -7.16339185
183 -0.86053356 1.26927222
184 -8.50093012 -0.86053356
185 -5.33069669 -8.50093012
186 -2.61405333 -5.33069669
187 -7.96708434 -2.61405333
188 -9.71788902 -7.96708434
189 -13.01004748 -9.71788902
190 -7.13130717 -13.01004748
191 -13.01004748 -7.13130717
192 -6.87961816 -13.01004748
193 -13.01004748 -6.87961816
194 -10.18633601 -13.01004748
195 -13.01004748 -10.18633601
196 39.07214005 -13.01004748
197 -13.01004748 39.07214005
198 -13.01004748 -13.01004748
199 0.47291583 -13.01004748
200 -13.01004748 0.47291583
201 -13.01004748 -13.01004748
202 -13.01004748 -13.01004748
203 0.40646372 -13.01004748
204 -9.26918373 0.40646372
205 -8.27238595 -9.26918373
206 -12.43257522 -8.27238595
207 -13.01004748 -12.43257522
208 9.61379766 -13.01004748
209 -11.12001887 9.61379766
210 -13.01004748 -11.12001887
211 -13.01004748 -13.01004748
212 -2.91061907 -13.01004748
213 -5.04913646 -2.91061907
214 -10.76776363 -5.04913646
215 -11.01004748 -10.76776363
216 -7.11145167 -11.01004748
217 -7.94132213 -7.11145167
218 -11.01004748 -7.94132213
219 6.87557740 -11.01004748
220 -10.99026770 6.87557740
221 -13.01004748 -10.99026770
222 -12.01923585 -13.01004748
223 4.80837403 -12.01923585
224 -12.12773964 4.80837403
225 7.52053233 -12.12773964
226 -4.28663498 7.52053233
227 -13.01004748 -4.28663498
228 -12.84459258 -13.01004748
229 -13.01004748 -12.84459258
230 -7.89455090 -13.01004748
231 -7.66017679 -7.89455090
232 -13.01004748 -7.66017679
233 -1.74921133 -13.01004748
234 -13.01004748 -1.74921133
235 -13.02656797 -13.01004748
236 -13.01004748 -13.02656797
237 -7.19268622 -13.01004748
238 -13.01004748 -7.19268622
239 -13.01004748 -13.01004748
240 -13.01004748 -13.01004748
241 -10.54977902 -13.01004748
242 -13.01004748 -10.54977902
243 -5.93745771 -13.01004748
244 -24.90757444 -5.93745771
245 -13.01004748 -24.90757444
246 -13.01004748 -13.01004748
247 1.25123907 -13.01004748
248 -13.01004748 1.25123907
249 -13.01004748 -13.01004748
250 -3.01918929 -13.01004748
251 -13.01004748 -3.01918929
252 -13.40966807 -13.01004748
253 19.59623438 -13.40966807
254 29.30913776 19.59623438
255 -13.01004748 29.30913776
256 -6.38979561 -13.01004748
257 14.34355929 -6.38979561
258 -13.01004748 14.34355929
259 -13.01004748 -13.01004748
260 -13.01004748 -13.01004748
261 -13.01004748 -13.01004748
262 -4.21946369 -13.01004748
263 -4.39771772 -4.21946369
264 -13.01004748 -4.39771772
265 -13.01004748 -13.01004748
266 -13.01004748 -13.01004748
267 -9.77221151 -13.01004748
268 -13.01004748 -9.77221151
269 10.08340022 -13.01004748
270 -13.01004748 10.08340022
271 -13.01004748 -13.01004748
272 -13.01004748 -13.01004748
273 5.47252280 -13.01004748
274 -13.01004748 5.47252280
275 -2.81884947 -13.01004748
276 7.55228145 -2.81884947
277 -13.01004748 7.55228145
278 8.65788379 -13.01004748
279 -7.34654199 8.65788379
280 -13.01004748 -7.34654199
281 -13.01004748 -13.01004748
282 12.68482492 -13.01004748
283 -13.01004748 12.68482492
284 -8.06128403 -13.01004748
285 -13.01004748 -8.06128403
286 -10.96261291 -13.01004748
287 -10.48807205 -10.96261291
288 -13.01004748 -10.48807205
289 0.59078953 -13.01004748
290 -13.01004748 0.59078953
291 -2.43904337 -13.01004748
292 -1.49388020 -2.43904337
293 -13.01004748 -1.49388020
294 -2.75902125 -13.01004748
295 -6.03662514 -2.75902125
296 -13.01004748 -6.03662514
297 -13.01004748 -13.01004748
298 -13.01004748 -13.01004748
299 -4.01004748 -13.01004748
300 -9.12842915 -4.01004748
301 -12.71831761 -9.12842915
302 -12.01004748 -12.71831761
303 -13.01004748 -12.01004748
304 -1.54679346 -13.01004748
305 -1.02309748 -1.54679346
306 -13.01004748 -1.02309748
307 4.26111624 -13.01004748
308 3.54315760 4.26111624
309 -13.01004748 3.54315760
310 -13.01004748 -13.01004748
311 11.00010465 -13.01004748
312 -13.01004748 11.00010465
313 -0.09368542 -13.01004748
314 -5.03840529 -0.09368542
315 -11.22820302 -5.03840529
316 -12.99930635 -11.22820302
317 5.47699173 -12.99930635
318 -8.01447432 5.47699173
319 -13.01004748 -8.01447432
320 -13.01004748 -13.01004748
321 -13.01004748 -13.01004748
322 -13.01004748 -13.01004748
323 -13.01004748 -13.01004748
324 -13.01004748 -13.01004748
325 -12.06132417 -13.01004748
326 -13.01004748 -12.06132417
327 -13.01004748 -13.01004748
328 -13.01004748 -13.01004748
329 -13.01004748 -13.01004748
330 -13.01004748 -13.01004748
331 -13.01004748 -13.01004748
332 -13.01004748 -13.01004748
333 -9.36600274 -13.01004748
334 -12.01648475 -9.36600274
335 -10.51453845 -12.01648475
336 3.26667748 -10.51453845
337 -13.01004748 3.26667748
338 -12.14374943 -13.01004748
339 -3.44259767 -12.14374943
340 -2.57123733 -3.44259767
341 -13.01004748 -2.57123733
342 -13.01004748 -13.01004748
343 -13.01004748 -13.01004748
344 -2.85074308 -13.01004748
345 -13.01004748 -2.85074308
346 -10.01653981 -13.01004748
347 -13.01004748 -10.01653981
348 -13.01004748 -13.01004748
349 -10.72879238 -13.01004748
350 -11.11509622 -10.72879238
351 -6.35852982 -11.11509622
352 -5.77851683 -6.35852982
353 -12.92176715 -5.77851683
354 -13.01004748 -12.92176715
355 -13.01004748 -13.01004748
356 -4.45943392 -13.01004748
357 -13.01004748 -4.45943392
358 -13.01004748 -13.01004748
359 -13.01004748 -13.01004748
360 -9.01560747 -13.01004748
361 -13.01004748 -9.01560747
362 -13.01004748 -13.01004748
363 6.48759728 -13.01004748
364 -4.66015046 6.48759728
365 -13.01004748 -4.66015046
366 -13.01004748 -13.01004748
367 -13.01004748 -13.01004748
368 -13.01004748 -13.01004748
369 -7.05796557 -13.01004748
370 -13.01004748 -7.05796557
371 -8.10779640 -13.01004748
372 -10.89878890 -8.10779640
373 -11.56783921 -10.89878890
374 -13.01004748 -11.56783921
375 -13.01004748 -13.01004748
376 -1.24126971 -13.01004748
377 -8.99011503 -1.24126971
378 -13.69479069 -8.99011503
379 -13.01004748 -13.69479069
380 -13.01004748 -13.01004748
381 -0.29870166 -13.01004748
382 3.02316739 -0.29870166
383 -13.01004748 3.02316739
384 -1.84328535 -13.01004748
385 4.01618617 -1.84328535
386 -10.63350187 4.01618617
387 -10.86463146 -10.63350187
388 -10.26417164 -10.86463146
389 -7.15312399 -10.26417164
390 -5.46879276 -7.15312399
391 -8.42799497 -5.46879276
392 -7.05170825 -8.42799497
393 -12.89307885 -7.05170825
394 -24.71203709 -12.89307885
395 -6.61715727 -24.71203709
396 -7.52414074 -6.61715727
397 203.84608517 -7.52414074
398 -8.85458702 203.84608517
399 27.88477530 -8.85458702
400 -16.95198399 27.88477530
401 -9.82122245 -16.95198399
402 -9.26138167 -9.82122245
403 -0.47417404 -9.26138167
404 -8.87467144 -0.47417404
405 -10.02281589 -8.87467144
406 -6.57189259 -10.02281589
407 -5.16203979 -6.57189259
408 -5.26905751 -5.16203979
409 -13.88120415 -5.26905751
410 35.90812542 -13.88120415
411 -27.56335233 35.90812542
412 -1.99717252 -27.56335233
413 -2.48676517 -1.99717252
414 -12.35691311 -2.48676517
415 -16.25330300 -12.35691311
416 -1.91803219 -16.25330300
417 1.68941382 -1.91803219
418 3.85325725 1.68941382
419 -14.98176598 3.85325725
420 -17.62138299 -14.98176598
421 46.22141168 -17.62138299
422 23.36785417 46.22141168
423 -16.45936332 23.36785417
424 -1.86649331 -16.45936332
425 45.05644409 -1.86649331
426 -3.98314145 45.05644409
427 39.85181744 -3.98314145
428 -28.18161075 39.85181744
429 -4.67298230 -28.18161075
430 -17.36361924 -4.67298230
> 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/7lv5a1292936115.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/8lv5a1292936115.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/9lv5a1292936115.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/10w54u1292936115.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/11h5li1292936115.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/12se231292936115.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/13gfhf1292936115.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/147j4r1292936115.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/155gwr1292936115.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/169hcx1292936115.tab")
+ }
>
> try(system("convert tmp/1plp11292936115.ps tmp/1plp11292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/2plp11292936115.ps tmp/2plp11292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/3plp11292936115.ps tmp/3plp11292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/4hd6m1292936115.ps tmp/4hd6m1292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/5hd6m1292936115.ps tmp/5hd6m1292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/6sm571292936115.ps tmp/6sm571292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/7lv5a1292936115.ps tmp/7lv5a1292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/8lv5a1292936115.ps tmp/8lv5a1292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/9lv5a1292936115.ps tmp/9lv5a1292936115.png",intern=TRUE))
character(0)
> try(system("convert tmp/10w54u1292936115.ps tmp/10w54u1292936115.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
11.597 2.085 26.507