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(1
+ ,162556
+ ,162556
+ ,213118
+ ,213118
+ ,6282929
+ ,1
+ ,29790
+ ,29790
+ ,81767
+ ,81767
+ ,4324047
+ ,1
+ ,87550
+ ,87550
+ ,153198
+ ,153198
+ ,4108272
+ ,0
+ ,84738
+ ,0
+ ,-26007
+ ,0
+ ,-1212617
+ ,1
+ ,54660
+ ,54660
+ ,126942
+ ,126942
+ ,1485329
+ ,1
+ ,42634
+ ,42634
+ ,157214
+ ,157214
+ ,1779876
+ ,0
+ ,40949
+ ,0
+ ,129352
+ ,0
+ ,1367203
+ ,1
+ ,42312
+ ,42312
+ ,234817
+ ,234817
+ ,2519076
+ ,1
+ ,37704
+ ,37704
+ ,60448
+ ,60448
+ ,912684
+ ,1
+ ,16275
+ ,16275
+ ,47818
+ ,47818
+ ,1443586
+ ,0
+ ,25830
+ ,0
+ ,245546
+ ,0
+ ,1220017
+ ,0
+ ,12679
+ ,0
+ ,48020
+ ,0
+ ,984885
+ ,1
+ ,18014
+ ,18014
+ ,-1710
+ ,-1710
+ ,1457425
+ ,0
+ ,43556
+ ,0
+ ,32648
+ ,0
+ ,-572920
+ ,1
+ ,24524
+ ,24524
+ ,95350
+ ,95350
+ ,929144
+ ,0
+ ,6532
+ ,0
+ ,151352
+ ,0
+ ,1151176
+ ,0
+ ,7123
+ ,0
+ ,288170
+ ,0
+ ,790090
+ ,1
+ ,20813
+ ,20813
+ ,114337
+ ,114337
+ ,774497
+ ,1
+ ,37597
+ ,37597
+ ,37884
+ ,37884
+ ,990576
+ ,0
+ ,17821
+ ,0
+ ,122844
+ ,0
+ ,454195
+ ,1
+ ,12988
+ ,12988
+ ,82340
+ ,82340
+ ,876607
+ ,1
+ ,22330
+ ,22330
+ ,79801
+ ,79801
+ ,711969
+ ,0
+ ,13326
+ ,0
+ ,165548
+ ,0
+ ,702380
+ ,0
+ ,16189
+ ,0
+ ,116384
+ ,0
+ ,264449
+ ,0
+ ,7146
+ ,0
+ ,134028
+ ,0
+ ,450033
+ ,0
+ ,15824
+ ,0
+ ,63838
+ ,0
+ ,541063
+ ,1
+ ,26088
+ ,26088
+ ,74996
+ ,74996
+ ,588864
+ ,0
+ ,11326
+ ,0
+ ,31080
+ ,0
+ ,-37216
+ ,0
+ ,8568
+ ,0
+ ,32168
+ ,0
+ ,783310
+ ,0
+ ,14416
+ ,0
+ ,49857
+ ,0
+ ,467359
+ ,1
+ ,3369
+ ,3369
+ ,87161
+ ,87161
+ ,688779
+ ,1
+ ,11819
+ ,11819
+ ,106113
+ ,106113
+ ,608419
+ ,1
+ ,6620
+ ,6620
+ ,80570
+ ,80570
+ ,696348
+ ,1
+ ,4519
+ ,4519
+ ,102129
+ ,102129
+ ,597793
+ ,0
+ ,2220
+ ,0
+ ,301670
+ ,0
+ ,821730
+ ,0
+ ,18562
+ ,0
+ ,102313
+ ,0
+ ,377934
+ ,0
+ ,10327
+ ,0
+ ,88577
+ ,0
+ ,651939
+ ,1
+ ,5336
+ ,5336
+ ,112477
+ ,112477
+ ,697458
+ ,1
+ ,2365
+ ,2365
+ ,191778
+ ,191778
+ ,700368
+ ,0
+ ,4069
+ ,0
+ ,79804
+ ,0
+ ,225986
+ ,0
+ ,7710
+ ,0
+ ,128294
+ ,0
+ ,348695
+ ,0
+ ,13718
+ ,0
+ ,96448
+ ,0
+ ,373683
+ ,0
+ ,4525
+ ,0
+ ,93811
+ ,0
+ ,501709
+ ,0
+ ,6869
+ ,0
+ ,117520
+ ,0
+ ,413743
+ ,0
+ ,4628
+ ,0
+ ,69159
+ ,0
+ ,379825
+ ,1
+ ,3653
+ ,3653
+ ,101792
+ ,101792
+ ,336260
+ ,1
+ ,1265
+ ,1265
+ ,210568
+ ,210568
+ ,636765
+ ,1
+ ,7489
+ ,7489
+ ,136996
+ ,136996
+ ,481231
+ ,0
+ ,4901
+ ,0
+ ,121920
+ ,0
+ ,469107
+ ,0
+ ,2284
+ ,0
+ ,76403
+ ,0
+ ,211928
+ ,1
+ ,3160
+ ,3160
+ ,108094
+ ,108094
+ ,563925
+ ,1
+ ,4150
+ ,4150
+ ,134759
+ ,134759
+ ,511939
+ ,1
+ ,7285
+ ,7285
+ ,188873
+ ,188873
+ ,521016
+ ,1
+ ,1134
+ ,1134
+ ,146216
+ ,146216
+ ,543856
+ ,1
+ ,4658
+ ,4658
+ ,156608
+ ,156608
+ ,329304
+ ,0
+ ,2384
+ ,0
+ ,61348
+ ,0
+ ,423262
+ ,0
+ ,3748
+ ,0
+ ,50350
+ ,0
+ ,509665
+ ,0
+ ,5371
+ ,0
+ ,87720
+ ,0
+ ,455881
+ ,0
+ ,1285
+ ,0
+ ,99489
+ ,0
+ ,367772
+ ,1
+ ,9327
+ ,9327
+ ,87419
+ ,87419
+ ,406339
+ ,1
+ ,5565
+ ,5565
+ ,94355
+ ,94355
+ ,493408
+ ,0
+ ,1528
+ ,0
+ ,60326
+ ,0
+ ,232942
+ ,1
+ ,3122
+ ,3122
+ ,94670
+ ,94670
+ ,416002
+ ,1
+ ,7317
+ ,7317
+ ,82425
+ ,82425
+ ,337430
+ ,0
+ ,2675
+ ,0
+ ,59017
+ ,0
+ ,361517
+ ,0
+ ,13253
+ ,0
+ ,90829
+ ,0
+ ,360962
+ ,0
+ ,880
+ ,0
+ ,80791
+ ,0
+ ,235561
+ ,1
+ ,2053
+ ,2053
+ ,100423
+ ,100423
+ ,408247
+ ,0
+ ,1424
+ ,0
+ ,131116
+ ,0
+ ,450296
+ ,1
+ ,4036
+ ,4036
+ ,100269
+ ,100269
+ ,418799
+ ,1
+ ,3045
+ ,3045
+ ,27330
+ ,27330
+ ,247405
+ ,0
+ ,5119
+ ,0
+ ,39039
+ ,0
+ ,378519
+ ,0
+ ,1431
+ ,0
+ ,106885
+ ,0
+ ,326638
+ ,0
+ ,554
+ ,0
+ ,79285
+ ,0
+ ,328233
+ ,0
+ ,1975
+ ,0
+ ,118881
+ ,0
+ ,386225
+ ,1
+ ,1286
+ ,1286
+ ,77623
+ ,77623
+ ,283662
+ ,0
+ ,1012
+ ,0
+ ,114768
+ ,0
+ ,370225
+ ,0
+ ,810
+ ,0
+ ,74015
+ ,0
+ ,269236
+ ,0
+ ,1280
+ ,0
+ ,69465
+ ,0
+ ,365732
+ ,1
+ ,666
+ ,666
+ ,117869
+ ,117869
+ ,420383
+ ,0
+ ,1380
+ ,0
+ ,60982
+ ,0
+ ,345811
+ ,1
+ ,4608
+ ,4608
+ ,90131
+ ,90131
+ ,431809
+ ,0
+ ,876
+ ,0
+ ,138971
+ ,0
+ ,418876
+ ,0
+ ,814
+ ,0
+ ,39625
+ ,0
+ ,297476
+ ,0
+ ,514
+ ,0
+ ,102725
+ ,0
+ ,416776
+ ,1
+ ,5692
+ ,5692
+ ,64239
+ ,64239
+ ,357257
+ ,0
+ ,3642
+ ,0
+ ,90262
+ ,0
+ ,458343
+ ,0
+ ,540
+ ,0
+ ,103960
+ ,0
+ ,388386
+ ,0
+ ,2099
+ ,0
+ ,106611
+ ,0
+ ,358934
+ ,0
+ ,567
+ ,0
+ ,103345
+ ,0
+ ,407560
+ ,0
+ ,2001
+ ,0
+ ,95551
+ ,0
+ ,392558
+ ,1
+ ,2949
+ ,2949
+ ,82903
+ ,82903
+ ,373177
+ ,0
+ ,2253
+ ,0
+ ,63593
+ ,0
+ ,428370
+ ,1
+ ,6533
+ ,6533
+ ,126910
+ ,126910
+ ,369419
+ ,0
+ ,1889
+ ,0
+ ,37527
+ ,0
+ ,358649
+ ,1
+ ,3055
+ ,3055
+ ,60247
+ ,60247
+ ,376641
+ ,0
+ ,272
+ ,0
+ ,112995
+ ,0
+ ,467427
+ ,1
+ ,1414
+ ,1414
+ ,70184
+ ,70184
+ ,364885
+ ,0
+ ,2564
+ ,0
+ ,130140
+ ,0
+ ,436230
+ ,1
+ ,1383
+ ,1383
+ ,73221
+ ,73221
+ ,329118
+ ,1
+ ,1261
+ ,1261
+ ,76114
+ ,76114
+ ,317365
+ ,0
+ ,975
+ ,0
+ ,90534
+ ,0
+ ,286849
+ ,0
+ ,3366
+ ,0
+ ,108479
+ ,0
+ ,376685
+ ,0
+ ,576
+ ,0
+ ,113761
+ ,0
+ ,407198
+ ,0
+ ,1306
+ ,0
+ ,68696
+ ,0
+ ,377772
+ ,0
+ ,746
+ ,0
+ ,71561
+ ,0
+ ,271483
+ ,1
+ ,3192
+ ,3192
+ ,59831
+ ,59831
+ ,153661
+ ,1
+ ,2045
+ ,2045
+ ,97890
+ ,97890
+ ,513294
+ ,0
+ ,5477
+ ,0
+ ,101481
+ ,0
+ ,324881
+ ,1
+ ,1932
+ ,1932
+ ,72954
+ ,72954
+ ,264512
+ ,0
+ ,936
+ ,0
+ ,67939
+ ,0
+ ,420968
+ ,1
+ ,3437
+ ,3437
+ ,48022
+ ,48022
+ ,129302
+ ,0
+ ,5131
+ ,0
+ ,86111
+ ,0
+ ,191521
+ ,1
+ ,2397
+ ,2397
+ ,74020
+ ,74020
+ ,268673
+ ,1
+ ,1389
+ ,1389
+ ,57530
+ ,57530
+ ,353179
+ ,0
+ ,1503
+ ,0
+ ,56364
+ ,0
+ ,354624
+ ,0
+ ,402
+ ,0
+ ,84990
+ ,0
+ ,363713
+ ,0
+ ,2239
+ ,0
+ ,88590
+ ,0
+ ,456657
+ ,1
+ ,2234
+ ,2234
+ ,77200
+ ,77200
+ ,211742
+ ,0
+ ,837
+ ,0
+ ,61262
+ ,0
+ ,338381
+ ,0
+ ,10579
+ ,0
+ ,110309
+ ,0
+ ,418530
+ ,0
+ ,875
+ ,0
+ ,67000
+ ,0
+ ,351483
+ ,0
+ ,1395
+ ,0
+ ,93099
+ ,0
+ ,372928
+ ,1
+ ,1659
+ ,1659
+ ,107577
+ ,107577
+ ,485538
+ ,1
+ ,2647
+ ,2647
+ ,62920
+ ,62920
+ ,279268
+ ,1
+ ,3294
+ ,3294
+ ,75832
+ ,75832
+ ,219060
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,0
+ ,60793
+ ,0
+ ,325314
+ ,0
+ ,422
+ ,0
+ ,57935
+ ,0
+ ,322046
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,34
+ ,0
+ ,60630
+ ,0
+ ,325599
+ ,0
+ ,1558
+ ,0
+ ,55637
+ ,0
+ ,377028
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,43
+ ,0
+ ,60887
+ ,0
+ ,323850
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,316
+ ,0
+ ,60505
+ ,0
+ ,331514
+ ,0
+ ,115
+ ,0
+ ,60945
+ ,0
+ ,325632
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,389
+ ,0
+ ,58990
+ ,0
+ ,322265
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,1002
+ ,0
+ ,56750
+ ,0
+ ,325906
+ ,0
+ ,36
+ ,0
+ ,60894
+ ,0
+ ,325985
+ ,0
+ ,460
+ ,0
+ ,63346
+ ,0
+ ,346145
+ ,0
+ ,309
+ ,0
+ ,56535
+ ,0
+ ,325898
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,9
+ ,0
+ ,60835
+ ,0
+ ,325356
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,14
+ ,0
+ ,61016
+ ,0
+ ,325930
+ ,0
+ ,520
+ ,0
+ ,58650
+ ,0
+ ,318020
+ ,0
+ ,1766
+ ,0
+ ,60438
+ ,0
+ ,326389
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,458
+ ,0
+ ,58625
+ ,0
+ ,302925
+ ,0
+ ,20
+ ,0
+ ,60938
+ ,0
+ ,325540
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,98
+ ,0
+ ,61490
+ ,0
+ ,326736
+ ,0
+ ,405
+ ,0
+ ,60845
+ ,0
+ ,340580
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,483
+ ,0
+ ,60830
+ ,0
+ ,331828
+ ,0
+ ,454
+ ,0
+ ,63261
+ ,0
+ ,323299
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,757
+ ,0
+ ,45689
+ ,0
+ ,387722
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,36
+ ,0
+ ,61564
+ ,0
+ ,324598
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,203
+ ,0
+ ,61938
+ ,0
+ ,328726
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,90
+ ,0
+ ,60951
+ ,0
+ ,325043
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,71
+ ,71
+ ,60745
+ ,60745
+ ,325806
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,972
+ ,0
+ ,71642
+ ,0
+ ,387732
+ ,1
+ ,531
+ ,531
+ ,71641
+ ,71641
+ ,349729
+ ,0
+ ,604
+ ,0
+ ,55792
+ ,0
+ ,332202
+ ,1
+ ,283
+ ,283
+ ,71873
+ ,71873
+ ,305442
+ ,1
+ ,23
+ ,23
+ ,62555
+ ,62555
+ ,329537
+ ,1
+ ,638
+ ,638
+ ,60370
+ ,60370
+ ,327055
+ ,1
+ ,699
+ ,699
+ ,64873
+ ,64873
+ ,356245
+ ,0
+ ,149
+ ,0
+ ,62041
+ ,0
+ ,328451
+ ,0
+ ,226
+ ,0
+ ,65745
+ ,0
+ ,307062
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,275
+ ,0
+ ,59500
+ ,0
+ ,331345
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,141
+ ,0
+ ,61630
+ ,0
+ ,331824
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,28
+ ,0
+ ,60890
+ ,0
+ ,325685
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,2566
+ ,2566
+ ,113521
+ ,113521
+ ,404480
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,472
+ ,472
+ ,80045
+ ,80045
+ ,318314
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,203
+ ,203
+ ,50804
+ ,50804
+ ,311807
+ ,1
+ ,496
+ ,496
+ ,87390
+ ,87390
+ ,337724
+ ,1
+ ,10
+ ,10
+ ,61656
+ ,61656
+ ,326431
+ ,1
+ ,63
+ ,63
+ ,65688
+ ,65688
+ ,327556
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,1136
+ ,1136
+ ,48522
+ ,48522
+ ,356850
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,267
+ ,0
+ ,57640
+ ,0
+ ,322741
+ ,0
+ ,474
+ ,0
+ ,61977
+ ,0
+ ,310902
+ ,0
+ ,534
+ ,0
+ ,62620
+ ,0
+ ,324295
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,15
+ ,0
+ ,60831
+ ,0
+ ,326156
+ ,0
+ ,397
+ ,0
+ ,60646
+ ,0
+ ,326960
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,1061
+ ,0
+ ,56225
+ ,0
+ ,333411
+ ,0
+ ,288
+ ,0
+ ,60510
+ ,0
+ ,297761
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,3
+ ,0
+ ,60698
+ ,0
+ ,325536
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,20
+ ,0
+ ,60805
+ ,0
+ ,325762
+ ,0
+ ,278
+ ,0
+ ,61404
+ ,0
+ ,327957
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,192
+ ,0
+ ,65276
+ ,0
+ ,318521
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,317
+ ,0
+ ,63915
+ ,0
+ ,319775
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,368
+ ,0
+ ,61686
+ ,0
+ ,332128
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,2
+ ,0
+ ,60743
+ ,0
+ ,325486
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,0
+ ,60349
+ ,0
+ ,325838
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,0
+ ,61360
+ ,0
+ ,331767
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,24
+ ,0
+ ,59818
+ ,0
+ ,324523
+ ,0
+ ,2332
+ ,0
+ ,72680
+ ,0
+ ,339995
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,131
+ ,0
+ ,61808
+ ,0
+ ,319582
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,206
+ ,0
+ ,53110
+ ,0
+ ,307245
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,167
+ ,0
+ ,64245
+ ,0
+ ,317967
+ ,0
+ ,622
+ ,0
+ ,73007
+ ,0
+ ,331488
+ ,0
+ ,885
+ ,0
+ ,82732
+ ,0
+ ,335452
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,365
+ ,0
+ ,54820
+ ,0
+ ,334184
+ ,0
+ ,364
+ ,0
+ ,47705
+ ,0
+ ,313213
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,226
+ ,0
+ ,72835
+ ,0
+ ,348678
+ ,0
+ ,307
+ ,0
+ ,58856
+ ,0
+ ,328727
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,188
+ ,0
+ ,77655
+ ,0
+ ,387978
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,138
+ ,0
+ ,69817
+ ,0
+ ,336704
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,125
+ ,0
+ ,60798
+ ,0
+ ,322076
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,282
+ ,0
+ ,62452
+ ,0
+ ,334272
+ ,0
+ ,335
+ ,0
+ ,64175
+ ,0
+ ,338197
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,813
+ ,813
+ ,67440
+ ,67440
+ ,321024
+ ,0
+ ,176
+ ,0
+ ,68136
+ ,0
+ ,322145
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,249
+ ,0
+ ,56726
+ ,0
+ ,323351
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,333
+ ,0
+ ,70811
+ ,0
+ ,327748
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,30
+ ,0
+ ,62045
+ ,0
+ ,328157
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,249
+ ,0
+ ,54323
+ ,0
+ ,311594
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,165
+ ,0
+ ,62841
+ ,0
+ ,335962
+ ,0
+ ,453
+ ,0
+ ,81125
+ ,0
+ ,372426
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,0
+ ,59506
+ ,0
+ ,319844
+ ,1
+ ,382
+ ,382
+ ,59365
+ ,59365
+ ,355822
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,30
+ ,30
+ ,60798
+ ,60798
+ ,324047
+ ,0
+ ,290
+ ,0
+ ,58790
+ ,0
+ ,311464
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,366
+ ,0
+ ,61808
+ ,0
+ ,353417
+ ,0
+ ,2
+ ,0
+ ,60735
+ ,0
+ ,325590
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,209
+ ,209
+ ,64016
+ ,64016
+ ,328576
+ ,0
+ ,384
+ ,0
+ ,54683
+ ,0
+ ,326126
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,365
+ ,0
+ ,87192
+ ,0
+ ,369376
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,49
+ ,49
+ ,64107
+ ,64107
+ ,332013
+ ,0
+ ,3
+ ,0
+ ,60761
+ ,0
+ ,325871
+ ,0
+ ,133
+ ,0
+ ,65990
+ ,0
+ ,342165
+ ,0
+ ,32
+ ,0
+ ,59988
+ ,0
+ ,324967
+ ,0
+ ,368
+ ,0
+ ,61167
+ ,0
+ ,314832
+ ,0
+ ,1
+ ,0
+ ,60719
+ ,0
+ ,325557
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,22
+ ,0
+ ,60722
+ ,0
+ ,322649
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,96
+ ,0
+ ,60379
+ ,0
+ ,324598
+ ,0
+ ,1
+ ,0
+ ,60727
+ ,0
+ ,325567
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,81
+ ,0
+ ,60925
+ ,0
+ ,324005
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,26
+ ,0
+ ,60896
+ ,0
+ ,325748
+ ,0
+ ,125
+ ,0
+ ,59734
+ ,0
+ ,323385
+ ,0
+ ,304
+ ,0
+ ,62969
+ ,0
+ ,315409
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,119
+ ,0
+ ,59118
+ ,0
+ ,312275
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,312
+ ,0
+ ,58598
+ ,0
+ ,320576
+ ,0
+ ,60
+ ,0
+ ,61124
+ ,0
+ ,325246
+ ,0
+ ,587
+ ,0
+ ,59595
+ ,0
+ ,332961
+ ,0
+ ,135
+ ,0
+ ,62065
+ ,0
+ ,323010
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,514
+ ,0
+ ,78780
+ ,0
+ ,345253
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,1
+ ,0
+ ,60722
+ ,0
+ ,325559
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,58
+ ,58
+ ,61600
+ ,61600
+ ,319634
+ ,0
+ ,180
+ ,0
+ ,59635
+ ,0
+ ,319951
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,448
+ ,0
+ ,59781
+ ,0
+ ,318519
+ ,0
+ ,227
+ ,0
+ ,76644
+ ,0
+ ,343222
+ ,0
+ ,174
+ ,0
+ ,64820
+ ,0
+ ,317234
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,121
+ ,0
+ ,56178
+ ,0
+ ,314025
+ ,0
+ ,607
+ ,0
+ ,60436
+ ,0
+ ,320249
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,530
+ ,0
+ ,73433
+ ,0
+ ,349365
+ ,0
+ ,571
+ ,0
+ ,41477
+ ,0
+ ,289197
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,78
+ ,0
+ ,62700
+ ,0
+ ,329245
+ ,0
+ ,2489
+ ,0
+ ,67804
+ ,0
+ ,240869
+ ,0
+ ,131
+ ,0
+ ,59661
+ ,0
+ ,327182
+ ,0
+ ,923
+ ,0
+ ,58620
+ ,0
+ ,322876
+ ,0
+ ,72
+ ,0
+ ,60398
+ ,0
+ ,323117
+ ,0
+ ,572
+ ,0
+ ,58580
+ ,0
+ ,306351
+ ,0
+ ,397
+ ,0
+ ,62710
+ ,0
+ ,335137
+ ,0
+ ,450
+ ,0
+ ,59325
+ ,0
+ ,308271
+ ,0
+ ,622
+ ,0
+ ,60950
+ ,0
+ ,301731
+ ,0
+ ,694
+ ,0
+ ,68060
+ ,0
+ ,382409
+ ,1
+ ,3425
+ ,3425
+ ,83620
+ ,83620
+ ,279230
+ ,0
+ ,562
+ ,0
+ ,58456
+ ,0
+ ,298731
+ ,0
+ ,4917
+ ,0
+ ,52811
+ ,0
+ ,243650
+ ,1
+ ,1442
+ ,1442
+ ,121173
+ ,121173
+ ,532682
+ ,0
+ ,529
+ ,0
+ ,63870
+ ,0
+ ,319771
+ ,1
+ ,2126
+ ,2126
+ ,21001
+ ,21001
+ ,171493
+ ,0
+ ,1061
+ ,0
+ ,70415
+ ,0
+ ,347262
+ ,0
+ ,776
+ ,0
+ ,64230
+ ,0
+ ,343945
+ ,0
+ ,611
+ ,0
+ ,59190
+ ,0
+ ,311874
+ ,1
+ ,1526
+ ,1526
+ ,69351
+ ,69351
+ ,302211
+ ,0
+ ,592
+ ,0
+ ,64270
+ ,0
+ ,316708
+ ,0
+ ,1182
+ ,0
+ ,70694
+ ,0
+ ,333463
+ ,0
+ ,621
+ ,0
+ ,68005
+ ,0
+ ,344282
+ ,0
+ ,989
+ ,0
+ ,58930
+ ,0
+ ,319635
+ ,0
+ ,438
+ ,0
+ ,58320
+ ,0
+ ,301186
+ ,0
+ ,726
+ ,0
+ ,69980
+ ,0
+ ,300381
+ ,0
+ ,1303
+ ,0
+ ,69863
+ ,0
+ ,318765
+ ,1
+ ,6341
+ ,6341
+ ,63255
+ ,63255
+ ,286146
+ ,1
+ ,1164
+ ,1164
+ ,57320
+ ,57320
+ ,306844
+ ,1
+ ,3310
+ ,3310
+ ,75230
+ ,75230
+ ,307705
+ ,0
+ ,1366
+ ,0
+ ,79420
+ ,0
+ ,312448
+ ,0
+ ,965
+ ,0
+ ,73490
+ ,0
+ ,299715
+ ,0
+ ,3256
+ ,0
+ ,35250
+ ,0
+ ,373399
+ ,1
+ ,1135
+ ,1135
+ ,62285
+ ,62285
+ ,299446
+ ,0
+ ,1270
+ ,0
+ ,69206
+ ,0
+ ,325586
+ ,0
+ ,661
+ ,0
+ ,65920
+ ,0
+ ,291221
+ ,0
+ ,1013
+ ,0
+ ,69770
+ ,0
+ ,261173
+ ,0
+ ,2844
+ ,0
+ ,72683
+ ,0
+ ,255027
+ ,1
+ ,11528
+ ,11528
+ ,-14545
+ ,-14545
+ ,-78375
+ ,0
+ ,6526
+ ,0
+ ,55830
+ ,0
+ ,-58143
+ ,0
+ ,2264
+ ,0
+ ,55174
+ ,0
+ ,227033
+ ,1
+ ,4461
+ ,4461
+ ,67038
+ ,67038
+ ,235098
+ ,0
+ ,3999
+ ,0
+ ,51252
+ ,0
+ ,21267
+ ,0
+ ,35624
+ ,0
+ ,157278
+ ,0
+ ,238675
+ ,0
+ ,9252
+ ,0
+ ,79510
+ ,0
+ ,197687
+ ,0
+ ,15236
+ ,0
+ ,77440
+ ,0
+ ,418341
+ ,0
+ ,18073
+ ,0
+ ,27284
+ ,0
+ ,-297706)
+ ,dim=c(6
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'GrCosts'
+ ,'Dividends'
+ ,'GrDiv'
+ ,'Wealth')
+ ,1:431))
> y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','GrCosts','Dividends','GrDiv','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 = '6'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Group Costs GrCosts Dividends GrDiv
1 6282929 1 162556 162556 213118 213118
2 4324047 1 29790 29790 81767 81767
3 4108272 1 87550 87550 153198 153198
4 -1212617 0 84738 0 -26007 0
5 1485329 1 54660 54660 126942 126942
6 1779876 1 42634 42634 157214 157214
7 1367203 0 40949 0 129352 0
8 2519076 1 42312 42312 234817 234817
9 912684 1 37704 37704 60448 60448
10 1443586 1 16275 16275 47818 47818
11 1220017 0 25830 0 245546 0
12 984885 0 12679 0 48020 0
13 1457425 1 18014 18014 -1710 -1710
14 -572920 0 43556 0 32648 0
15 929144 1 24524 24524 95350 95350
16 1151176 0 6532 0 151352 0
17 790090 0 7123 0 288170 0
18 774497 1 20813 20813 114337 114337
19 990576 1 37597 37597 37884 37884
20 454195 0 17821 0 122844 0
21 876607 1 12988 12988 82340 82340
22 711969 1 22330 22330 79801 79801
23 702380 0 13326 0 165548 0
24 264449 0 16189 0 116384 0
25 450033 0 7146 0 134028 0
26 541063 0 15824 0 63838 0
27 588864 1 26088 26088 74996 74996
28 -37216 0 11326 0 31080 0
29 783310 0 8568 0 32168 0
30 467359 0 14416 0 49857 0
31 688779 1 3369 3369 87161 87161
32 608419 1 11819 11819 106113 106113
33 696348 1 6620 6620 80570 80570
34 597793 1 4519 4519 102129 102129
35 821730 0 2220 0 301670 0
36 377934 0 18562 0 102313 0
37 651939 0 10327 0 88577 0
38 697458 1 5336 5336 112477 112477
39 700368 1 2365 2365 191778 191778
40 225986 0 4069 0 79804 0
41 348695 0 7710 0 128294 0
42 373683 0 13718 0 96448 0
43 501709 0 4525 0 93811 0
44 413743 0 6869 0 117520 0
45 379825 0 4628 0 69159 0
46 336260 1 3653 3653 101792 101792
47 636765 1 1265 1265 210568 210568
48 481231 1 7489 7489 136996 136996
49 469107 0 4901 0 121920 0
50 211928 0 2284 0 76403 0
51 563925 1 3160 3160 108094 108094
52 511939 1 4150 4150 134759 134759
53 521016 1 7285 7285 188873 188873
54 543856 1 1134 1134 146216 146216
55 329304 1 4658 4658 156608 156608
56 423262 0 2384 0 61348 0
57 509665 0 3748 0 50350 0
58 455881 0 5371 0 87720 0
59 367772 0 1285 0 99489 0
60 406339 1 9327 9327 87419 87419
61 493408 1 5565 5565 94355 94355
62 232942 0 1528 0 60326 0
63 416002 1 3122 3122 94670 94670
64 337430 1 7317 7317 82425 82425
65 361517 0 2675 0 59017 0
66 360962 0 13253 0 90829 0
67 235561 0 880 0 80791 0
68 408247 1 2053 2053 100423 100423
69 450296 0 1424 0 131116 0
70 418799 1 4036 4036 100269 100269
71 247405 1 3045 3045 27330 27330
72 378519 0 5119 0 39039 0
73 326638 0 1431 0 106885 0
74 328233 0 554 0 79285 0
75 386225 0 1975 0 118881 0
76 283662 1 1286 1286 77623 77623
77 370225 0 1012 0 114768 0
78 269236 0 810 0 74015 0
79 365732 0 1280 0 69465 0
80 420383 1 666 666 117869 117869
81 345811 0 1380 0 60982 0
82 431809 1 4608 4608 90131 90131
83 418876 0 876 0 138971 0
84 297476 0 814 0 39625 0
85 416776 0 514 0 102725 0
86 357257 1 5692 5692 64239 64239
87 458343 0 3642 0 90262 0
88 388386 0 540 0 103960 0
89 358934 0 2099 0 106611 0
90 407560 0 567 0 103345 0
91 392558 0 2001 0 95551 0
92 373177 1 2949 2949 82903 82903
93 428370 0 2253 0 63593 0
94 369419 1 6533 6533 126910 126910
95 358649 0 1889 0 37527 0
96 376641 1 3055 3055 60247 60247
97 467427 0 272 0 112995 0
98 364885 1 1414 1414 70184 70184
99 436230 0 2564 0 130140 0
100 329118 1 1383 1383 73221 73221
101 317365 1 1261 1261 76114 76114
102 286849 0 975 0 90534 0
103 376685 0 3366 0 108479 0
104 407198 0 576 0 113761 0
105 377772 0 1306 0 68696 0
106 271483 0 746 0 71561 0
107 153661 1 3192 3192 59831 59831
108 513294 1 2045 2045 97890 97890
109 324881 0 5477 0 101481 0
110 264512 1 1932 1932 72954 72954
111 420968 0 936 0 67939 0
112 129302 1 3437 3437 48022 48022
113 191521 0 5131 0 86111 0
114 268673 1 2397 2397 74020 74020
115 353179 1 1389 1389 57530 57530
116 354624 0 1503 0 56364 0
117 363713 0 402 0 84990 0
118 456657 0 2239 0 88590 0
119 211742 1 2234 2234 77200 77200
120 338381 0 837 0 61262 0
121 418530 0 10579 0 110309 0
122 351483 0 875 0 67000 0
123 372928 0 1395 0 93099 0
124 485538 1 1659 1659 107577 107577
125 279268 1 2647 2647 62920 62920
126 219060 1 3294 3294 75832 75832
127 325560 0 0 0 60720 0
128 325314 0 94 0 60793 0
129 322046 0 422 0 57935 0
130 325560 0 0 0 60720 0
131 325599 0 34 0 60630 0
132 377028 0 1558 0 55637 0
133 325560 0 0 0 60720 0
134 323850 0 43 0 60887 0
135 325560 0 0 0 60720 0
136 331514 0 316 0 60505 0
137 325632 0 115 0 60945 0
138 325560 0 0 0 60720 0
139 325560 0 0 0 60720 0
140 325560 0 0 0 60720 0
141 322265 0 389 0 58990 0
142 325560 0 0 0 60720 0
143 325906 0 1002 0 56750 0
144 325985 0 36 0 60894 0
145 346145 0 460 0 63346 0
146 325898 0 309 0 56535 0
147 325560 0 0 0 60720 0
148 325356 0 9 0 60835 0
149 325560 0 0 0 60720 0
150 325930 0 14 0 61016 0
151 318020 0 520 0 58650 0
152 326389 0 1766 0 60438 0
153 325560 0 0 0 60720 0
154 302925 0 458 0 58625 0
155 325540 0 20 0 60938 0
156 325560 0 0 0 60720 0
157 325560 0 0 0 60720 0
158 326736 0 98 0 61490 0
159 340580 0 405 0 60845 0
160 325560 0 0 0 60720 0
161 325560 0 0 0 60720 0
162 325560 0 0 0 60720 0
163 325560 0 0 0 60720 0
164 331828 0 483 0 60830 0
165 323299 0 454 0 63261 0
166 325560 0 0 0 60720 0
167 325560 0 0 0 60720 0
168 387722 0 757 0 45689 0
169 325560 0 0 0 60720 0
170 325560 0 0 0 60720 0
171 325560 0 0 0 60720 0
172 324598 0 36 0 61564 0
173 325560 0 0 0 60720 0
174 328726 0 203 0 61938 0
175 325560 0 0 0 60720 0
176 325043 0 90 0 60951 0
177 325560 0 0 0 60720 0
178 325806 1 71 71 60745 60745
179 325560 0 0 0 60720 0
180 325560 0 0 0 60720 0
181 387732 0 972 0 71642 0
182 349729 1 531 531 71641 71641
183 332202 0 604 0 55792 0
184 305442 1 283 283 71873 71873
185 329537 1 23 23 62555 62555
186 327055 1 638 638 60370 60370
187 356245 1 699 699 64873 64873
188 328451 0 149 0 62041 0
189 307062 0 226 0 65745 0
190 325560 0 0 0 60720 0
191 331345 0 275 0 59500 0
192 325560 0 0 0 60720 0
193 331824 0 141 0 61630 0
194 325560 0 0 0 60720 0
195 325685 0 28 0 60890 0
196 325560 1 0 0 60720 60720
197 404480 1 2566 2566 113521 113521
198 325560 1 0 0 60720 60720
199 325560 1 0 0 60720 60720
200 318314 1 472 472 80045 80045
201 325560 1 0 0 60720 60720
202 325560 1 0 0 60720 60720
203 325560 1 0 0 60720 60720
204 311807 1 203 203 50804 50804
205 337724 1 496 496 87390 87390
206 326431 1 10 10 61656 61656
207 327556 1 63 63 65688 65688
208 325560 1 0 0 60720 60720
209 356850 1 1136 1136 48522 48522
210 325560 1 0 0 60720 60720
211 325560 0 0 0 60720 0
212 325560 0 0 0 60720 0
213 322741 0 267 0 57640 0
214 310902 0 474 0 61977 0
215 324295 0 534 0 62620 0
216 325560 0 0 0 60720 0
217 326156 0 15 0 60831 0
218 326960 0 397 0 60646 0
219 325560 0 0 0 60720 0
220 333411 0 1061 0 56225 0
221 297761 0 288 0 60510 0
222 325560 0 0 0 60720 0
223 325536 0 3 0 60698 0
224 325560 0 0 0 60720 0
225 325762 0 20 0 60805 0
226 327957 0 278 0 61404 0
227 325560 0 0 0 60720 0
228 325560 0 0 0 60720 0
229 318521 0 192 0 65276 0
230 325560 0 0 0 60720 0
231 319775 0 317 0 63915 0
232 325560 0 0 0 60720 0
233 325560 0 0 0 60720 0
234 332128 0 368 0 61686 0
235 325560 0 0 0 60720 0
236 325486 0 2 0 60743 0
237 325560 0 0 0 60720 0
238 325838 0 53 0 60349 0
239 325560 0 0 0 60720 0
240 325560 0 0 0 60720 0
241 325560 0 0 0 60720 0
242 331767 0 94 0 61360 0
243 325560 0 0 0 60720 0
244 324523 0 24 0 59818 0
245 339995 0 2332 0 72680 0
246 325560 1 0 0 60720 60720
247 325560 1 0 0 60720 60720
248 319582 0 131 0 61808 0
249 325560 1 0 0 60720 60720
250 325560 1 0 0 60720 60720
251 307245 0 206 0 53110 0
252 325560 1 0 0 60720 60720
253 317967 0 167 0 64245 0
254 331488 0 622 0 73007 0
255 335452 0 885 0 82732 0
256 325560 0 0 0 60720 0
257 334184 0 365 0 54820 0
258 313213 0 364 0 47705 0
259 325560 1 0 0 60720 60720
260 325560 0 0 0 60720 0
261 325560 1 0 0 60720 60720
262 325560 1 0 0 60720 60720
263 348678 0 226 0 72835 0
264 328727 0 307 0 58856 0
265 325560 1 0 0 60720 60720
266 325560 0 0 0 60720 0
267 325560 1 0 0 60720 60720
268 387978 0 188 0 77655 0
269 325560 1 0 0 60720 60720
270 336704 0 138 0 69817 0
271 325560 1 0 0 60720 60720
272 325560 1 0 0 60720 60720
273 325560 1 0 0 60720 60720
274 322076 0 125 0 60798 0
275 325560 0 0 0 60720 0
276 334272 0 282 0 62452 0
277 338197 0 335 0 64175 0
278 325560 0 0 0 60720 0
279 321024 1 813 813 67440 67440
280 322145 0 176 0 68136 0
281 325560 1 0 0 60720 60720
282 325560 0 0 0 60720 0
283 323351 0 249 0 56726 0
284 325560 0 0 0 60720 0
285 327748 0 333 0 70811 0
286 325560 0 0 0 60720 0
287 325560 1 0 0 60720 60720
288 328157 0 30 0 62045 0
289 325560 1 0 0 60720 60720
290 311594 0 249 0 54323 0
291 325560 1 0 0 60720 60720
292 335962 0 165 0 62841 0
293 372426 0 453 0 81125 0
294 325560 0 0 0 60720 0
295 319844 0 53 0 59506 0
296 355822 1 382 382 59365 59365
297 325560 0 0 0 60720 0
298 325560 1 0 0 60720 60720
299 325560 0 0 0 60720 0
300 325560 0 0 0 60720 0
301 324047 1 30 30 60798 60798
302 311464 0 290 0 58790 0
303 325560 0 0 0 60720 0
304 325560 1 0 0 60720 60720
305 353417 0 366 0 61808 0
306 325590 0 2 0 60735 0
307 325560 1 0 0 60720 60720
308 328576 1 209 209 64016 64016
309 326126 0 384 0 54683 0
310 325560 1 0 0 60720 60720
311 325560 1 0 0 60720 60720
312 369376 0 365 0 87192 0
313 325560 1 0 0 60720 60720
314 332013 1 49 49 64107 64107
315 325871 0 3 0 60761 0
316 342165 0 133 0 65990 0
317 324967 0 32 0 59988 0
318 314832 0 368 0 61167 0
319 325557 0 1 0 60719 0
320 325560 0 0 0 60720 0
321 325560 0 0 0 60720 0
322 325560 0 0 0 60720 0
323 325560 0 0 0 60720 0
324 325560 1 0 0 60720 60720
325 325560 0 0 0 60720 0
326 322649 0 22 0 60722 0
327 325560 1 0 0 60720 60720
328 325560 1 0 0 60720 60720
329 325560 0 0 0 60720 0
330 325560 0 0 0 60720 0
331 325560 0 0 0 60720 0
332 325560 0 0 0 60720 0
333 325560 0 0 0 60720 0
334 324598 0 96 0 60379 0
335 325567 0 1 0 60727 0
336 325560 0 0 0 60720 0
337 324005 0 81 0 60925 0
338 325560 0 0 0 60720 0
339 325748 0 26 0 60896 0
340 323385 0 125 0 59734 0
341 315409 0 304 0 62969 0
342 325560 0 0 0 60720 0
343 325560 0 0 0 60720 0
344 325560 0 0 0 60720 0
345 325560 0 0 0 60720 0
346 325560 0 0 0 60720 0
347 312275 0 119 0 59118 0
348 325560 0 0 0 60720 0
349 325560 0 0 0 60720 0
350 325560 0 0 0 60720 0
351 320576 0 312 0 58598 0
352 325246 0 60 0 61124 0
353 332961 0 587 0 59595 0
354 323010 0 135 0 62065 0
355 325560 0 0 0 60720 0
356 325560 0 0 0 60720 0
357 345253 0 514 0 78780 0
358 325560 0 0 0 60720 0
359 325560 0 0 0 60720 0
360 325560 0 0 0 60720 0
361 325559 0 1 0 60722 0
362 325560 0 0 0 60720 0
363 325560 0 0 0 60720 0
364 319634 1 58 58 61600 61600
365 319951 0 180 0 59635 0
366 325560 0 0 0 60720 0
367 325560 0 0 0 60720 0
368 325560 0 0 0 60720 0
369 325560 0 0 0 60720 0
370 325560 0 0 0 60720 0
371 325560 0 0 0 60720 0
372 318519 0 448 0 59781 0
373 343222 0 227 0 76644 0
374 317234 0 174 0 64820 0
375 325560 0 0 0 60720 0
376 325560 0 0 0 60720 0
377 314025 0 121 0 56178 0
378 320249 0 607 0 60436 0
379 325560 0 0 0 60720 0
380 325560 0 0 0 60720 0
381 325560 0 0 0 60720 0
382 349365 0 530 0 73433 0
383 289197 0 571 0 41477 0
384 325560 0 0 0 60720 0
385 329245 0 78 0 62700 0
386 240869 0 2489 0 67804 0
387 327182 0 131 0 59661 0
388 322876 0 923 0 58620 0
389 323117 0 72 0 60398 0
390 306351 0 572 0 58580 0
391 335137 0 397 0 62710 0
392 308271 0 450 0 59325 0
393 301731 0 622 0 60950 0
394 382409 0 694 0 68060 0
395 279230 1 3425 3425 83620 83620
396 298731 0 562 0 58456 0
397 243650 0 4917 0 52811 0
398 532682 1 1442 1442 121173 121173
399 319771 0 529 0 63870 0
400 171493 1 2126 2126 21001 21001
401 347262 0 1061 0 70415 0
402 343945 0 776 0 64230 0
403 311874 0 611 0 59190 0
404 302211 1 1526 1526 69351 69351
405 316708 0 592 0 64270 0
406 333463 0 1182 0 70694 0
407 344282 0 621 0 68005 0
408 319635 0 989 0 58930 0
409 301186 0 438 0 58320 0
410 300381 0 726 0 69980 0
411 318765 0 1303 0 69863 0
412 286146 1 6341 6341 63255 63255
413 306844 1 1164 1164 57320 57320
414 307705 1 3310 3310 75230 75230
415 312448 0 1366 0 79420 0
416 299715 0 965 0 73490 0
417 373399 0 3256 0 35250 0
418 299446 1 1135 1135 62285 62285
419 325586 0 1270 0 69206 0
420 291221 0 661 0 65920 0
421 261173 0 1013 0 69770 0
422 255027 0 2844 0 72683 0
423 -78375 1 11528 11528 -14545 -14545
424 -58143 0 6526 0 55830 0
425 227033 0 2264 0 55174 0
426 235098 1 4461 4461 67038 67038
427 21267 0 3999 0 51252 0
428 238675 0 35624 0 157278 0
429 197687 0 9252 0 79510 0
430 418341 0 15236 0 77440 0
431 -297706 0 18073 0 27284 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs GrCosts Dividends GrDiv
117378.943 45851.725 -7.650 43.285 3.424 -1.613
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-855582.7 -18397.8 269.7 25822.7 2951183.2
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.174e+05 3.105e+04 3.780 0.000179 ***
Group 4.585e+04 5.561e+04 0.825 0.410105
Costs -7.650e+00 1.721e+00 -4.445 1.12e-05 ***
GrCosts 4.328e+01 2.065e+00 20.965 < 2e-16 ***
Dividends 3.424e+00 4.168e-01 8.215 2.59e-15 ***
GrDiv -1.613e+00 7.000e-01 -2.305 0.021671 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 210000 on 425 degrees of freedom
Multiple R-squared: 0.7822, Adjusted R-squared: 0.7796
F-statistic: 305.2 on 5 and 425 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 3.928338e-86 1.964169e-86
[2,] 1.0000000 7.940290e-92 3.970145e-92
[3,] 1.0000000 1.742086e-94 8.710430e-95
[4,] 1.0000000 1.382743e-105 6.913717e-106
[5,] 1.0000000 2.448957e-120 1.224478e-120
[6,] 1.0000000 3.577209e-129 1.788604e-129
[7,] 1.0000000 1.426528e-132 7.132642e-133
[8,] 1.0000000 2.813127e-143 1.406563e-143
[9,] 1.0000000 8.847252e-148 4.423626e-148
[10,] 1.0000000 1.496397e-149 7.481987e-150
[11,] 1.0000000 3.046286e-155 1.523143e-155
[12,] 1.0000000 6.495943e-155 3.247971e-155
[13,] 1.0000000 4.108295e-159 2.054147e-159
[14,] 1.0000000 1.390631e-161 6.953153e-162
[15,] 1.0000000 1.047230e-162 5.236151e-163
[16,] 1.0000000 6.338407e-163 3.169204e-163
[17,] 1.0000000 1.270486e-162 6.352429e-163
[18,] 1.0000000 1.959619e-164 9.798096e-165
[19,] 1.0000000 1.173600e-167 5.867998e-168
[20,] 1.0000000 9.010655e-170 4.505328e-170
[21,] 1.0000000 1.834766e-182 9.173829e-183
[22,] 1.0000000 1.405182e-184 7.025911e-185
[23,] 1.0000000 3.359450e-188 1.679725e-188
[24,] 1.0000000 1.034151e-189 5.170757e-190
[25,] 1.0000000 5.416736e-196 2.708368e-196
[26,] 1.0000000 2.345896e-197 1.172948e-197
[27,] 1.0000000 1.130549e-197 5.652746e-198
[28,] 1.0000000 2.421284e-197 1.210642e-197
[29,] 1.0000000 1.792957e-204 8.964787e-205
[30,] 1.0000000 2.213544e-209 1.106772e-209
[31,] 1.0000000 1.276771e-208 6.383857e-209
[32,] 1.0000000 4.138455e-209 2.069227e-209
[33,] 1.0000000 1.025320e-208 5.126599e-209
[34,] 1.0000000 1.991331e-208 9.956653e-209
[35,] 1.0000000 4.289176e-209 2.144588e-209
[36,] 1.0000000 2.473127e-208 1.236564e-208
[37,] 1.0000000 7.943502e-208 3.971751e-208
[38,] 1.0000000 3.095188e-207 1.547594e-207
[39,] 1.0000000 2.249467e-206 1.124734e-206
[40,] 1.0000000 9.188083e-206 4.594042e-206
[41,] 1.0000000 3.330634e-205 1.665317e-205
[42,] 1.0000000 9.456548e-206 4.728274e-206
[43,] 1.0000000 1.357956e-206 6.789779e-207
[44,] 1.0000000 8.643769e-206 4.321884e-206
[45,] 1.0000000 3.851816e-205 1.925908e-205
[46,] 1.0000000 4.092463e-204 2.046231e-204
[47,] 1.0000000 8.495437e-207 4.247718e-207
[48,] 1.0000000 6.144824e-207 3.072412e-207
[49,] 1.0000000 4.076781e-211 2.038390e-211
[50,] 1.0000000 1.229191e-211 6.145953e-212
[51,] 1.0000000 7.842217e-211 3.921108e-211
[52,] 1.0000000 2.662363e-211 1.331182e-211
[53,] 1.0000000 2.672890e-212 1.336445e-212
[54,] 1.0000000 5.689039e-212 2.844519e-212
[55,] 1.0000000 4.344623e-211 2.172312e-211
[56,] 1.0000000 1.043557e-210 5.217783e-211
[57,] 1.0000000 4.103002e-210 2.051501e-210
[58,] 1.0000000 1.922985e-210 9.614924e-211
[59,] 1.0000000 3.735838e-211 1.867919e-211
[60,] 1.0000000 4.247557e-210 2.123779e-210
[61,] 1.0000000 3.216746e-209 1.608373e-209
[62,] 1.0000000 1.922920e-208 9.614602e-209
[63,] 1.0000000 7.495868e-208 3.747934e-208
[64,] 1.0000000 1.919947e-209 9.599734e-210
[65,] 1.0000000 2.192575e-209 1.096288e-209
[66,] 1.0000000 1.566899e-208 7.834494e-209
[67,] 1.0000000 6.523269e-208 3.261634e-208
[68,] 1.0000000 2.368979e-207 1.184489e-207
[69,] 1.0000000 5.077727e-207 2.538863e-207
[70,] 1.0000000 9.815383e-207 4.907692e-207
[71,] 1.0000000 7.871039e-206 3.935519e-206
[72,] 1.0000000 4.154043e-205 2.077021e-205
[73,] 1.0000000 3.180501e-204 1.590250e-204
[74,] 1.0000000 2.828822e-204 1.414411e-204
[75,] 1.0000000 2.172224e-204 1.086112e-204
[76,] 1.0000000 1.601655e-203 8.008275e-204
[77,] 1.0000000 1.502963e-202 7.514817e-203
[78,] 1.0000000 2.546787e-203 1.273393e-203
[79,] 1.0000000 1.845281e-203 9.226405e-204
[80,] 1.0000000 1.051185e-202 5.255923e-203
[81,] 1.0000000 3.332448e-202 1.666224e-202
[82,] 1.0000000 2.673888e-201 1.336944e-201
[83,] 1.0000000 2.672146e-200 1.336073e-200
[84,] 1.0000000 1.675453e-199 8.377263e-200
[85,] 1.0000000 2.694082e-200 1.347041e-200
[86,] 1.0000000 1.583800e-199 7.919001e-200
[87,] 1.0000000 2.222162e-200 1.111081e-200
[88,] 1.0000000 7.114298e-201 3.557149e-201
[89,] 1.0000000 7.086767e-200 3.543384e-200
[90,] 1.0000000 4.437362e-199 2.218681e-199
[91,] 1.0000000 1.426962e-198 7.134811e-199
[92,] 1.0000000 1.470267e-197 7.351335e-198
[93,] 1.0000000 1.279797e-196 6.398987e-197
[94,] 1.0000000 2.264937e-197 1.132468e-197
[95,] 1.0000000 8.535786e-197 4.267893e-197
[96,] 1.0000000 7.159051e-197 3.579525e-197
[97,] 1.0000000 4.177779e-196 2.088889e-196
[98,] 1.0000000 5.938464e-196 2.969232e-196
[99,] 1.0000000 1.409756e-196 7.048780e-197
[100,] 1.0000000 8.485378e-198 4.242689e-198
[101,] 1.0000000 1.867575e-197 9.337874e-198
[102,] 1.0000000 7.653100e-197 3.826550e-197
[103,] 1.0000000 8.894456e-197 4.447228e-197
[104,] 1.0000000 2.892081e-197 1.446041e-197
[105,] 1.0000000 1.085784e-199 5.428922e-200
[106,] 1.0000000 5.269902e-199 2.634951e-199
[107,] 1.0000000 1.544519e-198 7.722594e-199
[108,] 1.0000000 3.206145e-198 1.603073e-198
[109,] 1.0000000 1.426600e-197 7.133002e-198
[110,] 1.0000000 5.069907e-197 2.534953e-197
[111,] 1.0000000 4.254588e-198 2.127294e-198
[112,] 1.0000000 3.629878e-197 1.814939e-197
[113,] 1.0000000 9.976852e-197 4.988426e-197
[114,] 1.0000000 9.432640e-196 4.716320e-196
[115,] 1.0000000 3.551342e-195 1.775671e-195
[116,] 1.0000000 6.585015e-195 3.292508e-195
[117,] 1.0000000 6.502460e-194 3.251230e-194
[118,] 1.0000000 3.583924e-194 1.791962e-194
[119,] 1.0000000 3.708065e-193 1.854033e-193
[120,] 1.0000000 3.810155e-192 1.905078e-192
[121,] 1.0000000 3.599949e-191 1.799975e-191
[122,] 1.0000000 3.663133e-190 1.831567e-190
[123,] 1.0000000 3.702144e-189 1.851072e-189
[124,] 1.0000000 1.329667e-189 6.648335e-190
[125,] 1.0000000 1.344765e-188 6.723827e-189
[126,] 1.0000000 1.353959e-187 6.769795e-188
[127,] 1.0000000 1.354422e-186 6.772108e-187
[128,] 1.0000000 1.264228e-185 6.321138e-186
[129,] 1.0000000 1.248822e-184 6.244110e-185
[130,] 1.0000000 1.229087e-183 6.145434e-184
[131,] 1.0000000 1.203124e-182 6.015622e-183
[132,] 1.0000000 1.171361e-181 5.856807e-182
[133,] 1.0000000 1.084620e-180 5.423102e-181
[134,] 1.0000000 1.044744e-179 5.223721e-180
[135,] 1.0000000 6.805478e-179 3.402739e-179
[136,] 1.0000000 6.489104e-178 3.244552e-178
[137,] 1.0000000 5.404755e-177 2.702378e-177
[138,] 1.0000000 4.191970e-176 2.095985e-176
[139,] 1.0000000 3.938757e-175 1.969378e-175
[140,] 1.0000000 3.684201e-174 1.842100e-174
[141,] 1.0000000 3.425382e-173 1.712691e-173
[142,] 1.0000000 3.170865e-172 1.585432e-172
[143,] 1.0000000 2.802378e-171 1.401189e-171
[144,] 1.0000000 1.838953e-170 9.194767e-171
[145,] 1.0000000 1.676088e-169 8.380441e-170
[146,] 1.0000000 1.508932e-168 7.544659e-169
[147,] 1.0000000 1.362555e-167 6.812776e-168
[148,] 1.0000000 1.222917e-166 6.114583e-167
[149,] 1.0000000 1.091946e-165 5.459731e-166
[150,] 1.0000000 9.705648e-165 4.852824e-165
[151,] 1.0000000 7.157104e-164 3.578552e-164
[152,] 1.0000000 6.295811e-163 3.147906e-163
[153,] 1.0000000 5.509979e-162 2.754989e-162
[154,] 1.0000000 4.797745e-161 2.398872e-161
[155,] 1.0000000 4.156410e-160 2.078205e-160
[156,] 1.0000000 3.298731e-159 1.649366e-159
[157,] 1.0000000 2.830485e-158 1.415242e-158
[158,] 1.0000000 2.415774e-157 1.207887e-157
[159,] 1.0000000 2.051471e-156 1.025735e-156
[160,] 1.0000000 3.197620e-158 1.598810e-158
[161,] 1.0000000 2.750689e-157 1.375344e-157
[162,] 1.0000000 2.354424e-156 1.177212e-156
[163,] 1.0000000 2.005213e-155 1.002606e-155
[164,] 1.0000000 1.707968e-154 8.539840e-155
[165,] 1.0000000 1.440299e-153 7.201497e-154
[166,] 1.0000000 1.203869e-152 6.019347e-153
[167,] 1.0000000 1.005179e-151 5.025894e-152
[168,] 1.0000000 8.352321e-151 4.176161e-151
[169,] 1.0000000 6.905243e-150 3.452622e-150
[170,] 1.0000000 5.522659e-149 2.761330e-149
[171,] 1.0000000 4.521798e-148 2.260899e-148
[172,] 1.0000000 3.684172e-147 1.842086e-147
[173,] 1.0000000 1.969256e-146 9.846281e-147
[174,] 1.0000000 1.554459e-145 7.772297e-146
[175,] 1.0000000 7.398860e-145 3.699430e-145
[176,] 1.0000000 3.163397e-144 1.581699e-144
[177,] 1.0000000 2.468818e-143 1.234409e-143
[178,] 1.0000000 1.858128e-142 9.290638e-143
[179,] 1.0000000 1.061147e-141 5.305733e-142
[180,] 1.0000000 8.384353e-141 4.192176e-141
[181,] 1.0000000 4.849314e-140 2.424657e-140
[182,] 1.0000000 3.799232e-139 1.899616e-139
[183,] 1.0000000 2.642958e-138 1.321479e-138
[184,] 1.0000000 2.051161e-137 1.025580e-137
[185,] 1.0000000 1.552772e-136 7.763860e-137
[186,] 1.0000000 1.193501e-135 5.967507e-136
[187,] 1.0000000 9.133105e-135 4.566552e-135
[188,] 1.0000000 6.828279e-134 3.414139e-134
[189,] 1.0000000 4.661675e-133 2.330837e-133
[190,] 1.0000000 3.458445e-132 1.729222e-132
[191,] 1.0000000 2.557236e-131 1.278618e-131
[192,] 1.0000000 8.332686e-131 4.166343e-131
[193,] 1.0000000 6.135147e-130 3.067573e-130
[194,] 1.0000000 4.501097e-129 2.250548e-129
[195,] 1.0000000 3.290254e-128 1.645127e-128
[196,] 1.0000000 2.211593e-127 1.105796e-127
[197,] 1.0000000 6.135759e-127 3.067879e-127
[198,] 1.0000000 4.446292e-126 2.223146e-126
[199,] 1.0000000 3.073861e-125 1.536930e-125
[200,] 1.0000000 2.212988e-124 1.106494e-124
[201,] 1.0000000 1.851559e-124 9.257793e-125
[202,] 1.0000000 1.336855e-123 6.684274e-124
[203,] 1.0000000 9.674995e-123 4.837498e-123
[204,] 1.0000000 6.968608e-122 3.484304e-122
[205,] 1.0000000 4.523724e-121 2.261862e-121
[206,] 1.0000000 3.189903e-120 1.594951e-120
[207,] 1.0000000 2.270371e-119 1.135186e-119
[208,] 1.0000000 1.605219e-118 8.026096e-119
[209,] 1.0000000 1.128320e-117 5.641599e-118
[210,] 1.0000000 7.588995e-117 3.794497e-117
[211,] 1.0000000 5.289653e-116 2.644827e-116
[212,] 1.0000000 1.899713e-115 9.498566e-116
[213,] 1.0000000 1.215582e-114 6.077910e-115
[214,] 1.0000000 8.379833e-114 4.189917e-114
[215,] 1.0000000 5.747219e-113 2.873609e-113
[216,] 1.0000000 3.924237e-112 1.962119e-112
[217,] 1.0000000 2.665412e-111 1.332706e-111
[218,] 1.0000000 1.774317e-110 8.871586e-111
[219,] 1.0000000 1.194242e-109 5.971211e-110
[220,] 1.0000000 7.999663e-109 3.999831e-109
[221,] 1.0000000 4.930728e-108 2.465364e-108
[222,] 1.0000000 3.274238e-107 1.637119e-107
[223,] 1.0000000 2.134144e-106 1.067072e-106
[224,] 1.0000000 1.404143e-105 7.020714e-106
[225,] 1.0000000 9.194109e-105 4.597054e-105
[226,] 1.0000000 5.741959e-104 2.870980e-104
[227,] 1.0000000 3.723910e-103 1.861955e-103
[228,] 1.0000000 2.404482e-102 1.202241e-102
[229,] 1.0000000 1.544430e-101 7.722151e-102
[230,] 1.0000000 9.775098e-101 4.887549e-101
[231,] 1.0000000 6.218180e-100 3.109090e-100
[232,] 1.0000000 3.936413e-99 1.968207e-99
[233,] 1.0000000 2.479880e-98 1.239940e-98
[234,] 1.0000000 1.521671e-97 7.608353e-98
[235,] 1.0000000 9.493580e-97 4.746790e-97
[236,] 1.0000000 5.823001e-96 2.911501e-96
[237,] 1.0000000 3.590292e-95 1.795146e-95
[238,] 1.0000000 2.203973e-94 1.101986e-94
[239,] 1.0000000 1.347161e-93 6.735807e-94
[240,] 1.0000000 8.268122e-93 4.134061e-93
[241,] 1.0000000 5.007893e-92 2.503947e-92
[242,] 1.0000000 3.019966e-91 1.509983e-91
[243,] 1.0000000 1.606507e-90 8.032535e-91
[244,] 1.0000000 9.606200e-90 4.803100e-90
[245,] 1.0000000 5.513644e-89 2.756822e-89
[246,] 1.0000000 2.641251e-88 1.320626e-88
[247,] 1.0000000 5.664043e-88 2.832021e-88
[248,] 1.0000000 3.353232e-87 1.676616e-87
[249,] 1.0000000 1.199015e-86 5.995076e-87
[250,] 1.0000000 3.241162e-86 1.620581e-86
[251,] 1.0000000 1.904104e-85 9.520519e-86
[252,] 1.0000000 1.111337e-84 5.556686e-85
[253,] 1.0000000 6.467634e-84 3.233817e-84
[254,] 1.0000000 3.746668e-83 1.873334e-83
[255,] 1.0000000 1.965936e-82 9.829681e-83
[256,] 1.0000000 1.006659e-81 5.033297e-82
[257,] 1.0000000 5.755550e-81 2.877775e-81
[258,] 1.0000000 3.265676e-80 1.632838e-80
[259,] 1.0000000 1.849259e-79 9.246295e-80
[260,] 1.0000000 1.046569e-78 5.232845e-79
[261,] 1.0000000 5.868955e-78 2.934478e-78
[262,] 1.0000000 2.986678e-77 1.493339e-77
[263,] 1.0000000 1.659591e-76 8.297957e-77
[264,] 1.0000000 9.177319e-76 4.588660e-76
[265,] 1.0000000 5.050318e-75 2.525159e-75
[266,] 1.0000000 2.768578e-74 1.384289e-74
[267,] 1.0000000 1.503327e-73 7.516636e-74
[268,] 1.0000000 7.947264e-73 3.973632e-73
[269,] 1.0000000 4.236644e-72 2.118322e-72
[270,] 1.0000000 2.264846e-71 1.132423e-71
[271,] 1.0000000 1.182934e-70 5.914672e-71
[272,] 1.0000000 5.398892e-70 2.699446e-70
[273,] 1.0000000 2.857594e-69 1.428797e-69
[274,] 1.0000000 1.499261e-68 7.496307e-69
[275,] 1.0000000 6.828905e-68 3.414452e-68
[276,] 1.0000000 3.546673e-67 1.773337e-67
[277,] 1.0000000 1.510713e-66 7.553564e-67
[278,] 1.0000000 7.780897e-66 3.890448e-66
[279,] 1.0000000 4.002323e-65 2.001161e-65
[280,] 1.0000000 2.052412e-64 1.026206e-64
[281,] 1.0000000 1.044705e-63 5.223523e-64
[282,] 1.0000000 4.649705e-63 2.324853e-63
[283,] 1.0000000 2.344256e-62 1.172128e-62
[284,] 1.0000000 1.156762e-61 5.783810e-62
[285,] 1.0000000 4.717028e-61 2.358514e-61
[286,] 1.0000000 2.333805e-60 1.166903e-60
[287,] 1.0000000 1.145721e-59 5.728604e-60
[288,] 1.0000000 4.001051e-59 2.000526e-59
[289,] 1.0000000 1.952411e-58 9.762054e-59
[290,] 1.0000000 9.501120e-58 4.750560e-58
[291,] 1.0000000 4.585141e-57 2.292571e-57
[292,] 1.0000000 2.200239e-56 1.100119e-56
[293,] 1.0000000 1.055220e-55 5.276101e-56
[294,] 1.0000000 5.011137e-55 2.505569e-55
[295,] 1.0000000 2.364138e-54 1.182069e-54
[296,] 1.0000000 1.112308e-53 5.561539e-54
[297,] 1.0000000 4.162377e-53 2.081189e-53
[298,] 1.0000000 1.932435e-52 9.662176e-53
[299,] 1.0000000 8.952636e-52 4.476318e-52
[300,] 1.0000000 4.144818e-51 2.072409e-51
[301,] 1.0000000 1.399651e-50 6.998256e-51
[302,] 1.0000000 6.387042e-50 3.193521e-50
[303,] 1.0000000 2.895279e-49 1.447640e-49
[304,] 1.0000000 6.045406e-49 3.022703e-49
[305,] 1.0000000 2.724850e-48 1.362425e-48
[306,] 1.0000000 1.226473e-47 6.132366e-48
[307,] 1.0000000 5.458510e-47 2.729255e-47
[308,] 1.0000000 2.438763e-46 1.219382e-46
[309,] 1.0000000 1.061359e-45 5.306793e-46
[310,] 1.0000000 4.681957e-45 2.340978e-45
[311,] 1.0000000 2.033932e-44 1.016966e-44
[312,] 1.0000000 8.780858e-44 4.390429e-44
[313,] 1.0000000 3.766972e-43 1.883486e-43
[314,] 1.0000000 1.605777e-42 8.028884e-43
[315,] 1.0000000 6.801385e-42 3.400693e-42
[316,] 1.0000000 2.860616e-41 1.430308e-41
[317,] 1.0000000 1.196369e-40 5.981843e-41
[318,] 1.0000000 4.997265e-40 2.498632e-40
[319,] 1.0000000 2.055771e-39 1.027886e-39
[320,] 1.0000000 8.355660e-39 4.177830e-39
[321,] 1.0000000 3.405775e-38 1.702888e-38
[322,] 1.0000000 1.378918e-37 6.894590e-38
[323,] 1.0000000 5.545323e-37 2.772661e-37
[324,] 1.0000000 2.214926e-36 1.107463e-36
[325,] 1.0000000 8.786449e-36 4.393224e-36
[326,] 1.0000000 3.440225e-35 1.720112e-35
[327,] 1.0000000 1.345946e-34 6.729730e-35
[328,] 1.0000000 5.228877e-34 2.614439e-34
[329,] 1.0000000 2.022646e-33 1.011323e-33
[330,] 1.0000000 7.746975e-33 3.873487e-33
[331,] 1.0000000 2.947632e-32 1.473816e-32
[332,] 1.0000000 1.097230e-31 5.486151e-32
[333,] 1.0000000 4.140778e-31 2.070389e-31
[334,] 1.0000000 1.541156e-30 7.705778e-31
[335,] 1.0000000 5.693466e-30 2.846733e-30
[336,] 1.0000000 2.087580e-29 1.043790e-29
[337,] 1.0000000 7.596533e-29 3.798266e-29
[338,] 1.0000000 2.743233e-28 1.371616e-28
[339,] 1.0000000 9.899384e-28 4.949692e-28
[340,] 1.0000000 3.519871e-27 1.759936e-27
[341,] 1.0000000 1.241716e-26 6.208580e-27
[342,] 1.0000000 4.345703e-26 2.172851e-26
[343,] 1.0000000 1.453869e-25 7.269346e-26
[344,] 1.0000000 5.022525e-25 2.511262e-25
[345,] 1.0000000 1.539108e-24 7.695542e-25
[346,] 1.0000000 5.284325e-24 2.642162e-24
[347,] 1.0000000 1.776639e-23 8.883194e-24
[348,] 1.0000000 5.922694e-23 2.961347e-23
[349,] 1.0000000 1.580674e-22 7.903370e-23
[350,] 1.0000000 5.210726e-22 2.605363e-22
[351,] 1.0000000 1.702699e-21 8.513497e-22
[352,] 1.0000000 5.514615e-21 2.757307e-21
[353,] 1.0000000 1.770012e-20 8.850058e-21
[354,] 1.0000000 5.629662e-20 2.814831e-20
[355,] 1.0000000 1.774095e-19 8.870477e-20
[356,] 1.0000000 5.614113e-19 2.807057e-19
[357,] 1.0000000 1.724782e-18 8.623908e-19
[358,] 1.0000000 5.285323e-18 2.642661e-18
[359,] 1.0000000 1.603918e-17 8.019591e-18
[360,] 1.0000000 4.819542e-17 2.409771e-17
[361,] 1.0000000 1.433770e-16 7.168849e-17
[362,] 1.0000000 4.222194e-16 2.111097e-16
[363,] 1.0000000 1.230593e-15 6.152963e-16
[364,] 1.0000000 3.501694e-15 1.750847e-15
[365,] 1.0000000 8.914665e-15 4.457332e-15
[366,] 1.0000000 2.551357e-14 1.275678e-14
[367,] 1.0000000 7.175477e-14 3.587739e-14
[368,] 1.0000000 1.995645e-13 9.978223e-14
[369,] 1.0000000 5.263768e-13 2.631884e-13
[370,] 1.0000000 1.408639e-12 7.043197e-13
[371,] 1.0000000 3.785000e-12 1.892500e-12
[372,] 1.0000000 1.004888e-11 5.024440e-12
[373,] 1.0000000 2.635440e-11 1.317720e-11
[374,] 1.0000000 6.900571e-11 3.450285e-11
[375,] 1.0000000 1.090952e-10 5.454762e-11
[376,] 1.0000000 2.781316e-10 1.390658e-10
[377,] 1.0000000 7.083336e-10 3.541668e-10
[378,] 1.0000000 1.463665e-09 7.318327e-10
[379,] 1.0000000 3.528912e-09 1.764456e-09
[380,] 1.0000000 7.904060e-09 3.952030e-09
[381,] 1.0000000 1.900859e-08 9.504296e-09
[382,] 1.0000000 4.492066e-08 2.246033e-08
[383,] 0.9999999 1.033325e-07 5.166626e-08
[384,] 0.9999999 2.391132e-07 1.195566e-07
[385,] 0.9999997 5.564545e-07 2.782272e-07
[386,] 0.9999994 1.103763e-06 5.518813e-07
[387,] 0.9999990 2.022141e-06 1.011070e-06
[388,] 0.9999978 4.460880e-06 2.230440e-06
[389,] 0.9999960 8.098457e-06 4.049229e-06
[390,] 0.9999916 1.682333e-05 8.411663e-06
[391,] 0.9999820 3.601403e-05 1.800701e-05
[392,] 0.9999622 7.554771e-05 3.777386e-05
[393,] 0.9999222 1.556781e-04 7.783906e-05
[394,] 0.9998525 2.949126e-04 1.474563e-04
[395,] 0.9997201 5.597361e-04 2.798680e-04
[396,] 0.9994492 1.101613e-03 5.508066e-04
[397,] 0.9989476 2.104760e-03 1.052380e-03
[398,] 0.9980163 3.967341e-03 1.983671e-03
[399,] 0.9964360 7.128077e-03 3.564038e-03
[400,] 0.9942112 1.157757e-02 5.788786e-03
[401,] 0.9902445 1.951092e-02 9.755458e-03
[402,] 0.9830383 3.392334e-02 1.696167e-02
[403,] 0.9715886 5.682286e-02 2.841143e-02
[404,] 0.9566296 8.674078e-02 4.337039e-02
[405,] 0.9306208 1.387585e-01 6.937924e-02
[406,] 0.8933875 2.132250e-01 1.066125e-01
[407,] 0.8425038 3.149924e-01 1.574962e-01
[408,] 0.7738484 4.523031e-01 2.261516e-01
[409,] 0.9108151 1.783697e-01 8.918486e-02
[410,] 0.8508713 2.982575e-01 1.491287e-01
[411,] 0.7709301 4.581398e-01 2.290699e-01
[412,] 0.6587239 6.825523e-01 3.412761e-01
[413,] 0.5139631 9.720739e-01 4.860369e-01
[414,] 0.3521063 7.042125e-01 6.478937e-01
> postscript(file="/var/www/html/rcomp/tmp/1jzi41291335824.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/2jzi41291335824.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/3jzi41291335824.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/4crz71291335824.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/5crz71291335824.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
-58916.68341 2951183.24003 547776.97946 -592723.57478 -855582.68318
6 7 8 9 10
-187303.82930 1120158.08279 422848.11045 -703591.95988 613805.40397
11 12 13 14 15
459456.04993 800071.51561 655359.62146 -468897.87437 -280661.11852
16 17 18 19 20
565520.53288 -259523.08584 -337447.59549 -581028.38399 52511.45770
21 22 23 24 25
101447.04410 -391496.71174 120088.18252 -127599.19151 -71606.42087
26 27 28 29 30
326145.45395 -639817.67921 -174375.28759 621327.36936 289543.00246
31 32 33 34 35
247663.51192 -168131.35366 151317.62041 88593.24434 -311615.01065
36 37 38 39 40
52219.10887 310262.12965 140406.31363 105590.98411 -133523.16907
41 42 43 44 45
-148996.18697 30995.22637 97726.72451 -53490.33926 61041.57451
46 47 48 49 50
-141469.50111 47161.91860 -196941.82015 -28247.02662 -149590.54896
51 52 53 54 55
92352.03082 -43197.40228 -243825.49661 75448.94087 -283498.95034
56 57 58 59 60
114058.25801 248553.73190 79226.58274 -80437.40445 -247557.75152
61 62 63 64 65
-38989.03249 -79310.46952 -29908.81807 -235797.07384 62520.90448
66 67 68 69 70
33957.12914 -151722.70630 -9987.33107 -105144.12475 -69820.92442
71 72 73 74 75
-73823.45845 166625.49732 -145779.19403 -56387.81872 -123106.25890
76 77 78 79 80
-65954.22615 -132389.61114 -95381.48058 20289.54469 19983.58345
81 82 83 84 85
30180.16223 -58836.44905 -167652.48064 50643.93943 -48411.72652
86 87 88 89 90
-125132.00610 59758.15692 -80831.59816 -107434.97899 -59345.23388
91 92 93 94 95
-36690.09821 -45261.40047 110477.03904 -256423.03336 127224.14916
96 97 98 99 100
-4549.51391 -34777.47183 24177.92868 -107147.52429 -15983.74517
101 102 103 104 105
-28627.86911 -133068.99569 -86388.00661 -95303.81873 35161.57075
106 107 108 109 110
-85221.31954 -231658.23841 99931.47560 -98081.59614 -99669.94893
111 112 113 114 115
78119.21937 -243364.25034 -181459.97158 -114009.58235 36276.51649
116 117 118 119 120
55746.57428 -41605.06999 53064.69920 -170890.36224 17637.61731
121 122 123 124 125
4368.30333 11382.81709 -52559.93239 68389.53214 -92223.61851
126 127 128 129 130
-198868.43790 269.66610 492.78022 9519.96729 269.66610
131 132 133 134 135
876.92556 81060.62993 269.66610 -1683.22084 269.66610
136 137 138 139 140
9377.15855 450.96100 269.66610 269.66610 269.66610
141 142 143 144 145
5874.10042 269.66610 21874.36284 374.26243 15381.84798
146 147 148 149 150
17301.29008 269.66610 -259.25801 269.66610 -266.77164
151 152 153 154 155
3795.40766 15573.66403 269.66610 -11688.27232 -343.79348
156 157 158 159 160
269.66610 269.66610 -441.21905 17959.78921 269.66610
161 162 163 164 165
269.66610 269.66610 269.66610 9855.82870 -7219.00178
166 167 168 169 170
269.66610 269.66610 119690.15922 269.66610 269.66610
171 172 173 174 175
269.66610 -3306.88500 269.66610 818.00422 269.66610
176 177 178 179 180
-349.82655 269.66610 50048.92944 269.66610 269.66610
181 182 183 184 185
32479.16452 37849.42693 28406.06360 1979.83537 52212.88945
186 187 188 189 190
31771.85850 50634.13806 -222.76290 -33705.60342 269.66610
191 192 193 194 195
12335.74127 269.66610 4496.34472 269.66610 26.76109
196 197 198 199 200
52378.29308 -55752.85014 52378.29308 52378.29308 -6680.96187
201 202 203 204 205
52378.29308 52378.29308 52378.29308 49347.13500 -1426.44177
206 207 208 209 210
51198.04438 43133.28603 52378.29308 65274.78044 52378.29308
211 212 213 214 215
269.66610 269.66610 10039.37072 -15066.46219 -3416.17584
216 217 218 219 220
269.66610 600.33670 4959.98750 269.66610 31628.34903
221 222 223 224 225
-24607.15405 269.66610 343.94547 269.66610 333.61190
226 227 228 229 230
2451.20280 269.66610 269.66610 -20900.79062 269.66610
231 232 233 234 235
-14030.37493 269.66610 269.66610 6345.08101 269.66610
236 237 238 239 240
132.21122 269.66610 2223.44248 269.66610 269.66610
241 242 243 244 245
269.66610 5004.31516 269.66610 2504.79811 -8408.43581
246 247 248 249 250
52378.29308 52378.29308 -8431.64237 52378.29308 52378.29308
251 252 253 254 255
9587.91386 52378.29308 -18115.78580 -31116.13407 -58439.63706
256 257 258 259 260
269.66610 31888.00688 35271.83324 52378.29308 269.66610
261 262 263 264 265
52378.29308 52378.29308 -16366.47699 12167.65305 52378.29308
266 267 268 269 270
269.66610 52378.29308 6138.66863 52378.29308 -18679.71641
271 272 273 274 275
52378.29308 52378.29308 52378.29308 -2525.19903 269.66610
276 277 278 279 280
5208.34419 3639.05382 269.66610 6702.44125 -27192.11395
281 282 283 284 285
52378.29308 269.66610 13641.30379 269.66610 -29547.57748
286 287 288 289 290
269.66610 52378.29308 -1440.77572 52378.29308 10112.41764
291 292 293 294 295
52378.29308 4671.35207 -19267.78540 269.66610 -884.04067
296 297 298 299 300
71481.29254 269.66610 52378.29308 269.66610 269.66610
301 302 303 304 305
49654.99785 -4999.40151 269.66610 52378.29308 27201.04130
306 307 308 309 310
263.60403 52378.29308 41978.19496 24444.45331 52378.29308
311 312 313 314 315
52378.29308 -43764.97916 52378.29308 50952.03341 463.22713
316 317 318 319 320
-152.93183 2427.89874 -9173.81075 277.73992 269.66610
321 322 323 324 325
269.66610 269.66610 269.66610 52378.29308 269.66610
326 327 328 329 330
-2479.88831 52378.29308 52378.29308 269.66610 269.66610
331 332 333 334 335
269.66610 269.66610 269.66610 1209.65733 260.34712
336 337 338 339 340
269.66610 -1367.64739 269.66610 53.91705 2427.04406
341 342 343 344 345
-15256.62205 269.66610 269.66610 269.66610 269.66610
346 347 348 349 350
269.66610 -6619.60825 269.66610 269.66610 269.66610
351 352 353 354 355
4938.31961 -968.68749 16013.16365 -5853.03737 269.66610
356 357 358 359 360
269.66610 -37944.63660 269.66610 269.66610 269.66610
361 362 363 364 365
269.46762 269.66610 269.66610 42791.96245 -247.23550
366 367 368 369 370
269.66610 269.66610 269.66610 269.66610 269.66610
371 372 373 374 375
269.66610 -129.02984 -34857.22662 -20764.09565 269.66610
376 377 378 379 380
269.66610 5212.54708 574.48936 269.66610 269.66610
381 382 383 384 385
269.66610 -15401.57498 34164.62363 269.66610 -2228.37519
386 387 388 389 390
-89637.51538 6519.90171 11836.96694 -519.99381 -7236.11997
391 392 393 394 395
6069.64377 -8800.34051 -19588.75261 37294.67149 -157469.05530
396 397 398 399 400
-14508.02866 -16945.45973 98646.88263 -12258.55024 -105526.30134
401 402 403 404 405
-3108.63911 12572.25382 -3503.48237 -40978.81913 -16209.25828
406 407 408 409 410
-16937.34734 -1202.43237 8039.37711 -12535.91598 -51062.81077
411 412 413 414 415
-27864.30386 -217588.38755 -1660.31473 -109703.50562 -66423.50149
416 417 418 419 420
-61919.12149 160227.99049 -17015.45784 -19046.11042 -46818.19381
421 422 423 424 425
-87356.28065 -89470.05262 -626069.48963 -316767.42372 -61948.31074
426 427 428 429 430
-208492.56134 -241012.72772 -144726.09880 -121166.99639 152350.80295
431
-370254.75590
> postscript(file="/var/www/html/rcomp/tmp/6crz71291335824.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 -58916.68341 NA
1 2951183.24003 -58916.68341
2 547776.97946 2951183.24003
3 -592723.57478 547776.97946
4 -855582.68318 -592723.57478
5 -187303.82930 -855582.68318
6 1120158.08279 -187303.82930
7 422848.11045 1120158.08279
8 -703591.95988 422848.11045
9 613805.40397 -703591.95988
10 459456.04993 613805.40397
11 800071.51561 459456.04993
12 655359.62146 800071.51561
13 -468897.87437 655359.62146
14 -280661.11852 -468897.87437
15 565520.53288 -280661.11852
16 -259523.08584 565520.53288
17 -337447.59549 -259523.08584
18 -581028.38399 -337447.59549
19 52511.45770 -581028.38399
20 101447.04410 52511.45770
21 -391496.71174 101447.04410
22 120088.18252 -391496.71174
23 -127599.19151 120088.18252
24 -71606.42087 -127599.19151
25 326145.45395 -71606.42087
26 -639817.67921 326145.45395
27 -174375.28759 -639817.67921
28 621327.36936 -174375.28759
29 289543.00246 621327.36936
30 247663.51192 289543.00246
31 -168131.35366 247663.51192
32 151317.62041 -168131.35366
33 88593.24434 151317.62041
34 -311615.01065 88593.24434
35 52219.10887 -311615.01065
36 310262.12965 52219.10887
37 140406.31363 310262.12965
38 105590.98411 140406.31363
39 -133523.16907 105590.98411
40 -148996.18697 -133523.16907
41 30995.22637 -148996.18697
42 97726.72451 30995.22637
43 -53490.33926 97726.72451
44 61041.57451 -53490.33926
45 -141469.50111 61041.57451
46 47161.91860 -141469.50111
47 -196941.82015 47161.91860
48 -28247.02662 -196941.82015
49 -149590.54896 -28247.02662
50 92352.03082 -149590.54896
51 -43197.40228 92352.03082
52 -243825.49661 -43197.40228
53 75448.94087 -243825.49661
54 -283498.95034 75448.94087
55 114058.25801 -283498.95034
56 248553.73190 114058.25801
57 79226.58274 248553.73190
58 -80437.40445 79226.58274
59 -247557.75152 -80437.40445
60 -38989.03249 -247557.75152
61 -79310.46952 -38989.03249
62 -29908.81807 -79310.46952
63 -235797.07384 -29908.81807
64 62520.90448 -235797.07384
65 33957.12914 62520.90448
66 -151722.70630 33957.12914
67 -9987.33107 -151722.70630
68 -105144.12475 -9987.33107
69 -69820.92442 -105144.12475
70 -73823.45845 -69820.92442
71 166625.49732 -73823.45845
72 -145779.19403 166625.49732
73 -56387.81872 -145779.19403
74 -123106.25890 -56387.81872
75 -65954.22615 -123106.25890
76 -132389.61114 -65954.22615
77 -95381.48058 -132389.61114
78 20289.54469 -95381.48058
79 19983.58345 20289.54469
80 30180.16223 19983.58345
81 -58836.44905 30180.16223
82 -167652.48064 -58836.44905
83 50643.93943 -167652.48064
84 -48411.72652 50643.93943
85 -125132.00610 -48411.72652
86 59758.15692 -125132.00610
87 -80831.59816 59758.15692
88 -107434.97899 -80831.59816
89 -59345.23388 -107434.97899
90 -36690.09821 -59345.23388
91 -45261.40047 -36690.09821
92 110477.03904 -45261.40047
93 -256423.03336 110477.03904
94 127224.14916 -256423.03336
95 -4549.51391 127224.14916
96 -34777.47183 -4549.51391
97 24177.92868 -34777.47183
98 -107147.52429 24177.92868
99 -15983.74517 -107147.52429
100 -28627.86911 -15983.74517
101 -133068.99569 -28627.86911
102 -86388.00661 -133068.99569
103 -95303.81873 -86388.00661
104 35161.57075 -95303.81873
105 -85221.31954 35161.57075
106 -231658.23841 -85221.31954
107 99931.47560 -231658.23841
108 -98081.59614 99931.47560
109 -99669.94893 -98081.59614
110 78119.21937 -99669.94893
111 -243364.25034 78119.21937
112 -181459.97158 -243364.25034
113 -114009.58235 -181459.97158
114 36276.51649 -114009.58235
115 55746.57428 36276.51649
116 -41605.06999 55746.57428
117 53064.69920 -41605.06999
118 -170890.36224 53064.69920
119 17637.61731 -170890.36224
120 4368.30333 17637.61731
121 11382.81709 4368.30333
122 -52559.93239 11382.81709
123 68389.53214 -52559.93239
124 -92223.61851 68389.53214
125 -198868.43790 -92223.61851
126 269.66610 -198868.43790
127 492.78022 269.66610
128 9519.96729 492.78022
129 269.66610 9519.96729
130 876.92556 269.66610
131 81060.62993 876.92556
132 269.66610 81060.62993
133 -1683.22084 269.66610
134 269.66610 -1683.22084
135 9377.15855 269.66610
136 450.96100 9377.15855
137 269.66610 450.96100
138 269.66610 269.66610
139 269.66610 269.66610
140 5874.10042 269.66610
141 269.66610 5874.10042
142 21874.36284 269.66610
143 374.26243 21874.36284
144 15381.84798 374.26243
145 17301.29008 15381.84798
146 269.66610 17301.29008
147 -259.25801 269.66610
148 269.66610 -259.25801
149 -266.77164 269.66610
150 3795.40766 -266.77164
151 15573.66403 3795.40766
152 269.66610 15573.66403
153 -11688.27232 269.66610
154 -343.79348 -11688.27232
155 269.66610 -343.79348
156 269.66610 269.66610
157 -441.21905 269.66610
158 17959.78921 -441.21905
159 269.66610 17959.78921
160 269.66610 269.66610
161 269.66610 269.66610
162 269.66610 269.66610
163 9855.82870 269.66610
164 -7219.00178 9855.82870
165 269.66610 -7219.00178
166 269.66610 269.66610
167 119690.15922 269.66610
168 269.66610 119690.15922
169 269.66610 269.66610
170 269.66610 269.66610
171 -3306.88500 269.66610
172 269.66610 -3306.88500
173 818.00422 269.66610
174 269.66610 818.00422
175 -349.82655 269.66610
176 269.66610 -349.82655
177 50048.92944 269.66610
178 269.66610 50048.92944
179 269.66610 269.66610
180 32479.16452 269.66610
181 37849.42693 32479.16452
182 28406.06360 37849.42693
183 1979.83537 28406.06360
184 52212.88945 1979.83537
185 31771.85850 52212.88945
186 50634.13806 31771.85850
187 -222.76290 50634.13806
188 -33705.60342 -222.76290
189 269.66610 -33705.60342
190 12335.74127 269.66610
191 269.66610 12335.74127
192 4496.34472 269.66610
193 269.66610 4496.34472
194 26.76109 269.66610
195 52378.29308 26.76109
196 -55752.85014 52378.29308
197 52378.29308 -55752.85014
198 52378.29308 52378.29308
199 -6680.96187 52378.29308
200 52378.29308 -6680.96187
201 52378.29308 52378.29308
202 52378.29308 52378.29308
203 49347.13500 52378.29308
204 -1426.44177 49347.13500
205 51198.04438 -1426.44177
206 43133.28603 51198.04438
207 52378.29308 43133.28603
208 65274.78044 52378.29308
209 52378.29308 65274.78044
210 269.66610 52378.29308
211 269.66610 269.66610
212 10039.37072 269.66610
213 -15066.46219 10039.37072
214 -3416.17584 -15066.46219
215 269.66610 -3416.17584
216 600.33670 269.66610
217 4959.98750 600.33670
218 269.66610 4959.98750
219 31628.34903 269.66610
220 -24607.15405 31628.34903
221 269.66610 -24607.15405
222 343.94547 269.66610
223 269.66610 343.94547
224 333.61190 269.66610
225 2451.20280 333.61190
226 269.66610 2451.20280
227 269.66610 269.66610
228 -20900.79062 269.66610
229 269.66610 -20900.79062
230 -14030.37493 269.66610
231 269.66610 -14030.37493
232 269.66610 269.66610
233 6345.08101 269.66610
234 269.66610 6345.08101
235 132.21122 269.66610
236 269.66610 132.21122
237 2223.44248 269.66610
238 269.66610 2223.44248
239 269.66610 269.66610
240 269.66610 269.66610
241 5004.31516 269.66610
242 269.66610 5004.31516
243 2504.79811 269.66610
244 -8408.43581 2504.79811
245 52378.29308 -8408.43581
246 52378.29308 52378.29308
247 -8431.64237 52378.29308
248 52378.29308 -8431.64237
249 52378.29308 52378.29308
250 9587.91386 52378.29308
251 52378.29308 9587.91386
252 -18115.78580 52378.29308
253 -31116.13407 -18115.78580
254 -58439.63706 -31116.13407
255 269.66610 -58439.63706
256 31888.00688 269.66610
257 35271.83324 31888.00688
258 52378.29308 35271.83324
259 269.66610 52378.29308
260 52378.29308 269.66610
261 52378.29308 52378.29308
262 -16366.47699 52378.29308
263 12167.65305 -16366.47699
264 52378.29308 12167.65305
265 269.66610 52378.29308
266 52378.29308 269.66610
267 6138.66863 52378.29308
268 52378.29308 6138.66863
269 -18679.71641 52378.29308
270 52378.29308 -18679.71641
271 52378.29308 52378.29308
272 52378.29308 52378.29308
273 -2525.19903 52378.29308
274 269.66610 -2525.19903
275 5208.34419 269.66610
276 3639.05382 5208.34419
277 269.66610 3639.05382
278 6702.44125 269.66610
279 -27192.11395 6702.44125
280 52378.29308 -27192.11395
281 269.66610 52378.29308
282 13641.30379 269.66610
283 269.66610 13641.30379
284 -29547.57748 269.66610
285 269.66610 -29547.57748
286 52378.29308 269.66610
287 -1440.77572 52378.29308
288 52378.29308 -1440.77572
289 10112.41764 52378.29308
290 52378.29308 10112.41764
291 4671.35207 52378.29308
292 -19267.78540 4671.35207
293 269.66610 -19267.78540
294 -884.04067 269.66610
295 71481.29254 -884.04067
296 269.66610 71481.29254
297 52378.29308 269.66610
298 269.66610 52378.29308
299 269.66610 269.66610
300 49654.99785 269.66610
301 -4999.40151 49654.99785
302 269.66610 -4999.40151
303 52378.29308 269.66610
304 27201.04130 52378.29308
305 263.60403 27201.04130
306 52378.29308 263.60403
307 41978.19496 52378.29308
308 24444.45331 41978.19496
309 52378.29308 24444.45331
310 52378.29308 52378.29308
311 -43764.97916 52378.29308
312 52378.29308 -43764.97916
313 50952.03341 52378.29308
314 463.22713 50952.03341
315 -152.93183 463.22713
316 2427.89874 -152.93183
317 -9173.81075 2427.89874
318 277.73992 -9173.81075
319 269.66610 277.73992
320 269.66610 269.66610
321 269.66610 269.66610
322 269.66610 269.66610
323 52378.29308 269.66610
324 269.66610 52378.29308
325 -2479.88831 269.66610
326 52378.29308 -2479.88831
327 52378.29308 52378.29308
328 269.66610 52378.29308
329 269.66610 269.66610
330 269.66610 269.66610
331 269.66610 269.66610
332 269.66610 269.66610
333 1209.65733 269.66610
334 260.34712 1209.65733
335 269.66610 260.34712
336 -1367.64739 269.66610
337 269.66610 -1367.64739
338 53.91705 269.66610
339 2427.04406 53.91705
340 -15256.62205 2427.04406
341 269.66610 -15256.62205
342 269.66610 269.66610
343 269.66610 269.66610
344 269.66610 269.66610
345 269.66610 269.66610
346 -6619.60825 269.66610
347 269.66610 -6619.60825
348 269.66610 269.66610
349 269.66610 269.66610
350 4938.31961 269.66610
351 -968.68749 4938.31961
352 16013.16365 -968.68749
353 -5853.03737 16013.16365
354 269.66610 -5853.03737
355 269.66610 269.66610
356 -37944.63660 269.66610
357 269.66610 -37944.63660
358 269.66610 269.66610
359 269.66610 269.66610
360 269.46762 269.66610
361 269.66610 269.46762
362 269.66610 269.66610
363 42791.96245 269.66610
364 -247.23550 42791.96245
365 269.66610 -247.23550
366 269.66610 269.66610
367 269.66610 269.66610
368 269.66610 269.66610
369 269.66610 269.66610
370 269.66610 269.66610
371 -129.02984 269.66610
372 -34857.22662 -129.02984
373 -20764.09565 -34857.22662
374 269.66610 -20764.09565
375 269.66610 269.66610
376 5212.54708 269.66610
377 574.48936 5212.54708
378 269.66610 574.48936
379 269.66610 269.66610
380 269.66610 269.66610
381 -15401.57498 269.66610
382 34164.62363 -15401.57498
383 269.66610 34164.62363
384 -2228.37519 269.66610
385 -89637.51538 -2228.37519
386 6519.90171 -89637.51538
387 11836.96694 6519.90171
388 -519.99381 11836.96694
389 -7236.11997 -519.99381
390 6069.64377 -7236.11997
391 -8800.34051 6069.64377
392 -19588.75261 -8800.34051
393 37294.67149 -19588.75261
394 -157469.05530 37294.67149
395 -14508.02866 -157469.05530
396 -16945.45973 -14508.02866
397 98646.88263 -16945.45973
398 -12258.55024 98646.88263
399 -105526.30134 -12258.55024
400 -3108.63911 -105526.30134
401 12572.25382 -3108.63911
402 -3503.48237 12572.25382
403 -40978.81913 -3503.48237
404 -16209.25828 -40978.81913
405 -16937.34734 -16209.25828
406 -1202.43237 -16937.34734
407 8039.37711 -1202.43237
408 -12535.91598 8039.37711
409 -51062.81077 -12535.91598
410 -27864.30386 -51062.81077
411 -217588.38755 -27864.30386
412 -1660.31473 -217588.38755
413 -109703.50562 -1660.31473
414 -66423.50149 -109703.50562
415 -61919.12149 -66423.50149
416 160227.99049 -61919.12149
417 -17015.45784 160227.99049
418 -19046.11042 -17015.45784
419 -46818.19381 -19046.11042
420 -87356.28065 -46818.19381
421 -89470.05262 -87356.28065
422 -626069.48963 -89470.05262
423 -316767.42372 -626069.48963
424 -61948.31074 -316767.42372
425 -208492.56134 -61948.31074
426 -241012.72772 -208492.56134
427 -144726.09880 -241012.72772
428 -121166.99639 -144726.09880
429 152350.80295 -121166.99639
430 -370254.75590 152350.80295
431 NA -370254.75590
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2951183.24003 -58916.68341
[2,] 547776.97946 2951183.24003
[3,] -592723.57478 547776.97946
[4,] -855582.68318 -592723.57478
[5,] -187303.82930 -855582.68318
[6,] 1120158.08279 -187303.82930
[7,] 422848.11045 1120158.08279
[8,] -703591.95988 422848.11045
[9,] 613805.40397 -703591.95988
[10,] 459456.04993 613805.40397
[11,] 800071.51561 459456.04993
[12,] 655359.62146 800071.51561
[13,] -468897.87437 655359.62146
[14,] -280661.11852 -468897.87437
[15,] 565520.53288 -280661.11852
[16,] -259523.08584 565520.53288
[17,] -337447.59549 -259523.08584
[18,] -581028.38399 -337447.59549
[19,] 52511.45770 -581028.38399
[20,] 101447.04410 52511.45770
[21,] -391496.71174 101447.04410
[22,] 120088.18252 -391496.71174
[23,] -127599.19151 120088.18252
[24,] -71606.42087 -127599.19151
[25,] 326145.45395 -71606.42087
[26,] -639817.67921 326145.45395
[27,] -174375.28759 -639817.67921
[28,] 621327.36936 -174375.28759
[29,] 289543.00246 621327.36936
[30,] 247663.51192 289543.00246
[31,] -168131.35366 247663.51192
[32,] 151317.62041 -168131.35366
[33,] 88593.24434 151317.62041
[34,] -311615.01065 88593.24434
[35,] 52219.10887 -311615.01065
[36,] 310262.12965 52219.10887
[37,] 140406.31363 310262.12965
[38,] 105590.98411 140406.31363
[39,] -133523.16907 105590.98411
[40,] -148996.18697 -133523.16907
[41,] 30995.22637 -148996.18697
[42,] 97726.72451 30995.22637
[43,] -53490.33926 97726.72451
[44,] 61041.57451 -53490.33926
[45,] -141469.50111 61041.57451
[46,] 47161.91860 -141469.50111
[47,] -196941.82015 47161.91860
[48,] -28247.02662 -196941.82015
[49,] -149590.54896 -28247.02662
[50,] 92352.03082 -149590.54896
[51,] -43197.40228 92352.03082
[52,] -243825.49661 -43197.40228
[53,] 75448.94087 -243825.49661
[54,] -283498.95034 75448.94087
[55,] 114058.25801 -283498.95034
[56,] 248553.73190 114058.25801
[57,] 79226.58274 248553.73190
[58,] -80437.40445 79226.58274
[59,] -247557.75152 -80437.40445
[60,] -38989.03249 -247557.75152
[61,] -79310.46952 -38989.03249
[62,] -29908.81807 -79310.46952
[63,] -235797.07384 -29908.81807
[64,] 62520.90448 -235797.07384
[65,] 33957.12914 62520.90448
[66,] -151722.70630 33957.12914
[67,] -9987.33107 -151722.70630
[68,] -105144.12475 -9987.33107
[69,] -69820.92442 -105144.12475
[70,] -73823.45845 -69820.92442
[71,] 166625.49732 -73823.45845
[72,] -145779.19403 166625.49732
[73,] -56387.81872 -145779.19403
[74,] -123106.25890 -56387.81872
[75,] -65954.22615 -123106.25890
[76,] -132389.61114 -65954.22615
[77,] -95381.48058 -132389.61114
[78,] 20289.54469 -95381.48058
[79,] 19983.58345 20289.54469
[80,] 30180.16223 19983.58345
[81,] -58836.44905 30180.16223
[82,] -167652.48064 -58836.44905
[83,] 50643.93943 -167652.48064
[84,] -48411.72652 50643.93943
[85,] -125132.00610 -48411.72652
[86,] 59758.15692 -125132.00610
[87,] -80831.59816 59758.15692
[88,] -107434.97899 -80831.59816
[89,] -59345.23388 -107434.97899
[90,] -36690.09821 -59345.23388
[91,] -45261.40047 -36690.09821
[92,] 110477.03904 -45261.40047
[93,] -256423.03336 110477.03904
[94,] 127224.14916 -256423.03336
[95,] -4549.51391 127224.14916
[96,] -34777.47183 -4549.51391
[97,] 24177.92868 -34777.47183
[98,] -107147.52429 24177.92868
[99,] -15983.74517 -107147.52429
[100,] -28627.86911 -15983.74517
[101,] -133068.99569 -28627.86911
[102,] -86388.00661 -133068.99569
[103,] -95303.81873 -86388.00661
[104,] 35161.57075 -95303.81873
[105,] -85221.31954 35161.57075
[106,] -231658.23841 -85221.31954
[107,] 99931.47560 -231658.23841
[108,] -98081.59614 99931.47560
[109,] -99669.94893 -98081.59614
[110,] 78119.21937 -99669.94893
[111,] -243364.25034 78119.21937
[112,] -181459.97158 -243364.25034
[113,] -114009.58235 -181459.97158
[114,] 36276.51649 -114009.58235
[115,] 55746.57428 36276.51649
[116,] -41605.06999 55746.57428
[117,] 53064.69920 -41605.06999
[118,] -170890.36224 53064.69920
[119,] 17637.61731 -170890.36224
[120,] 4368.30333 17637.61731
[121,] 11382.81709 4368.30333
[122,] -52559.93239 11382.81709
[123,] 68389.53214 -52559.93239
[124,] -92223.61851 68389.53214
[125,] -198868.43790 -92223.61851
[126,] 269.66610 -198868.43790
[127,] 492.78022 269.66610
[128,] 9519.96729 492.78022
[129,] 269.66610 9519.96729
[130,] 876.92556 269.66610
[131,] 81060.62993 876.92556
[132,] 269.66610 81060.62993
[133,] -1683.22084 269.66610
[134,] 269.66610 -1683.22084
[135,] 9377.15855 269.66610
[136,] 450.96100 9377.15855
[137,] 269.66610 450.96100
[138,] 269.66610 269.66610
[139,] 269.66610 269.66610
[140,] 5874.10042 269.66610
[141,] 269.66610 5874.10042
[142,] 21874.36284 269.66610
[143,] 374.26243 21874.36284
[144,] 15381.84798 374.26243
[145,] 17301.29008 15381.84798
[146,] 269.66610 17301.29008
[147,] -259.25801 269.66610
[148,] 269.66610 -259.25801
[149,] -266.77164 269.66610
[150,] 3795.40766 -266.77164
[151,] 15573.66403 3795.40766
[152,] 269.66610 15573.66403
[153,] -11688.27232 269.66610
[154,] -343.79348 -11688.27232
[155,] 269.66610 -343.79348
[156,] 269.66610 269.66610
[157,] -441.21905 269.66610
[158,] 17959.78921 -441.21905
[159,] 269.66610 17959.78921
[160,] 269.66610 269.66610
[161,] 269.66610 269.66610
[162,] 269.66610 269.66610
[163,] 9855.82870 269.66610
[164,] -7219.00178 9855.82870
[165,] 269.66610 -7219.00178
[166,] 269.66610 269.66610
[167,] 119690.15922 269.66610
[168,] 269.66610 119690.15922
[169,] 269.66610 269.66610
[170,] 269.66610 269.66610
[171,] -3306.88500 269.66610
[172,] 269.66610 -3306.88500
[173,] 818.00422 269.66610
[174,] 269.66610 818.00422
[175,] -349.82655 269.66610
[176,] 269.66610 -349.82655
[177,] 50048.92944 269.66610
[178,] 269.66610 50048.92944
[179,] 269.66610 269.66610
[180,] 32479.16452 269.66610
[181,] 37849.42693 32479.16452
[182,] 28406.06360 37849.42693
[183,] 1979.83537 28406.06360
[184,] 52212.88945 1979.83537
[185,] 31771.85850 52212.88945
[186,] 50634.13806 31771.85850
[187,] -222.76290 50634.13806
[188,] -33705.60342 -222.76290
[189,] 269.66610 -33705.60342
[190,] 12335.74127 269.66610
[191,] 269.66610 12335.74127
[192,] 4496.34472 269.66610
[193,] 269.66610 4496.34472
[194,] 26.76109 269.66610
[195,] 52378.29308 26.76109
[196,] -55752.85014 52378.29308
[197,] 52378.29308 -55752.85014
[198,] 52378.29308 52378.29308
[199,] -6680.96187 52378.29308
[200,] 52378.29308 -6680.96187
[201,] 52378.29308 52378.29308
[202,] 52378.29308 52378.29308
[203,] 49347.13500 52378.29308
[204,] -1426.44177 49347.13500
[205,] 51198.04438 -1426.44177
[206,] 43133.28603 51198.04438
[207,] 52378.29308 43133.28603
[208,] 65274.78044 52378.29308
[209,] 52378.29308 65274.78044
[210,] 269.66610 52378.29308
[211,] 269.66610 269.66610
[212,] 10039.37072 269.66610
[213,] -15066.46219 10039.37072
[214,] -3416.17584 -15066.46219
[215,] 269.66610 -3416.17584
[216,] 600.33670 269.66610
[217,] 4959.98750 600.33670
[218,] 269.66610 4959.98750
[219,] 31628.34903 269.66610
[220,] -24607.15405 31628.34903
[221,] 269.66610 -24607.15405
[222,] 343.94547 269.66610
[223,] 269.66610 343.94547
[224,] 333.61190 269.66610
[225,] 2451.20280 333.61190
[226,] 269.66610 2451.20280
[227,] 269.66610 269.66610
[228,] -20900.79062 269.66610
[229,] 269.66610 -20900.79062
[230,] -14030.37493 269.66610
[231,] 269.66610 -14030.37493
[232,] 269.66610 269.66610
[233,] 6345.08101 269.66610
[234,] 269.66610 6345.08101
[235,] 132.21122 269.66610
[236,] 269.66610 132.21122
[237,] 2223.44248 269.66610
[238,] 269.66610 2223.44248
[239,] 269.66610 269.66610
[240,] 269.66610 269.66610
[241,] 5004.31516 269.66610
[242,] 269.66610 5004.31516
[243,] 2504.79811 269.66610
[244,] -8408.43581 2504.79811
[245,] 52378.29308 -8408.43581
[246,] 52378.29308 52378.29308
[247,] -8431.64237 52378.29308
[248,] 52378.29308 -8431.64237
[249,] 52378.29308 52378.29308
[250,] 9587.91386 52378.29308
[251,] 52378.29308 9587.91386
[252,] -18115.78580 52378.29308
[253,] -31116.13407 -18115.78580
[254,] -58439.63706 -31116.13407
[255,] 269.66610 -58439.63706
[256,] 31888.00688 269.66610
[257,] 35271.83324 31888.00688
[258,] 52378.29308 35271.83324
[259,] 269.66610 52378.29308
[260,] 52378.29308 269.66610
[261,] 52378.29308 52378.29308
[262,] -16366.47699 52378.29308
[263,] 12167.65305 -16366.47699
[264,] 52378.29308 12167.65305
[265,] 269.66610 52378.29308
[266,] 52378.29308 269.66610
[267,] 6138.66863 52378.29308
[268,] 52378.29308 6138.66863
[269,] -18679.71641 52378.29308
[270,] 52378.29308 -18679.71641
[271,] 52378.29308 52378.29308
[272,] 52378.29308 52378.29308
[273,] -2525.19903 52378.29308
[274,] 269.66610 -2525.19903
[275,] 5208.34419 269.66610
[276,] 3639.05382 5208.34419
[277,] 269.66610 3639.05382
[278,] 6702.44125 269.66610
[279,] -27192.11395 6702.44125
[280,] 52378.29308 -27192.11395
[281,] 269.66610 52378.29308
[282,] 13641.30379 269.66610
[283,] 269.66610 13641.30379
[284,] -29547.57748 269.66610
[285,] 269.66610 -29547.57748
[286,] 52378.29308 269.66610
[287,] -1440.77572 52378.29308
[288,] 52378.29308 -1440.77572
[289,] 10112.41764 52378.29308
[290,] 52378.29308 10112.41764
[291,] 4671.35207 52378.29308
[292,] -19267.78540 4671.35207
[293,] 269.66610 -19267.78540
[294,] -884.04067 269.66610
[295,] 71481.29254 -884.04067
[296,] 269.66610 71481.29254
[297,] 52378.29308 269.66610
[298,] 269.66610 52378.29308
[299,] 269.66610 269.66610
[300,] 49654.99785 269.66610
[301,] -4999.40151 49654.99785
[302,] 269.66610 -4999.40151
[303,] 52378.29308 269.66610
[304,] 27201.04130 52378.29308
[305,] 263.60403 27201.04130
[306,] 52378.29308 263.60403
[307,] 41978.19496 52378.29308
[308,] 24444.45331 41978.19496
[309,] 52378.29308 24444.45331
[310,] 52378.29308 52378.29308
[311,] -43764.97916 52378.29308
[312,] 52378.29308 -43764.97916
[313,] 50952.03341 52378.29308
[314,] 463.22713 50952.03341
[315,] -152.93183 463.22713
[316,] 2427.89874 -152.93183
[317,] -9173.81075 2427.89874
[318,] 277.73992 -9173.81075
[319,] 269.66610 277.73992
[320,] 269.66610 269.66610
[321,] 269.66610 269.66610
[322,] 269.66610 269.66610
[323,] 52378.29308 269.66610
[324,] 269.66610 52378.29308
[325,] -2479.88831 269.66610
[326,] 52378.29308 -2479.88831
[327,] 52378.29308 52378.29308
[328,] 269.66610 52378.29308
[329,] 269.66610 269.66610
[330,] 269.66610 269.66610
[331,] 269.66610 269.66610
[332,] 269.66610 269.66610
[333,] 1209.65733 269.66610
[334,] 260.34712 1209.65733
[335,] 269.66610 260.34712
[336,] -1367.64739 269.66610
[337,] 269.66610 -1367.64739
[338,] 53.91705 269.66610
[339,] 2427.04406 53.91705
[340,] -15256.62205 2427.04406
[341,] 269.66610 -15256.62205
[342,] 269.66610 269.66610
[343,] 269.66610 269.66610
[344,] 269.66610 269.66610
[345,] 269.66610 269.66610
[346,] -6619.60825 269.66610
[347,] 269.66610 -6619.60825
[348,] 269.66610 269.66610
[349,] 269.66610 269.66610
[350,] 4938.31961 269.66610
[351,] -968.68749 4938.31961
[352,] 16013.16365 -968.68749
[353,] -5853.03737 16013.16365
[354,] 269.66610 -5853.03737
[355,] 269.66610 269.66610
[356,] -37944.63660 269.66610
[357,] 269.66610 -37944.63660
[358,] 269.66610 269.66610
[359,] 269.66610 269.66610
[360,] 269.46762 269.66610
[361,] 269.66610 269.46762
[362,] 269.66610 269.66610
[363,] 42791.96245 269.66610
[364,] -247.23550 42791.96245
[365,] 269.66610 -247.23550
[366,] 269.66610 269.66610
[367,] 269.66610 269.66610
[368,] 269.66610 269.66610
[369,] 269.66610 269.66610
[370,] 269.66610 269.66610
[371,] -129.02984 269.66610
[372,] -34857.22662 -129.02984
[373,] -20764.09565 -34857.22662
[374,] 269.66610 -20764.09565
[375,] 269.66610 269.66610
[376,] 5212.54708 269.66610
[377,] 574.48936 5212.54708
[378,] 269.66610 574.48936
[379,] 269.66610 269.66610
[380,] 269.66610 269.66610
[381,] -15401.57498 269.66610
[382,] 34164.62363 -15401.57498
[383,] 269.66610 34164.62363
[384,] -2228.37519 269.66610
[385,] -89637.51538 -2228.37519
[386,] 6519.90171 -89637.51538
[387,] 11836.96694 6519.90171
[388,] -519.99381 11836.96694
[389,] -7236.11997 -519.99381
[390,] 6069.64377 -7236.11997
[391,] -8800.34051 6069.64377
[392,] -19588.75261 -8800.34051
[393,] 37294.67149 -19588.75261
[394,] -157469.05530 37294.67149
[395,] -14508.02866 -157469.05530
[396,] -16945.45973 -14508.02866
[397,] 98646.88263 -16945.45973
[398,] -12258.55024 98646.88263
[399,] -105526.30134 -12258.55024
[400,] -3108.63911 -105526.30134
[401,] 12572.25382 -3108.63911
[402,] -3503.48237 12572.25382
[403,] -40978.81913 -3503.48237
[404,] -16209.25828 -40978.81913
[405,] -16937.34734 -16209.25828
[406,] -1202.43237 -16937.34734
[407,] 8039.37711 -1202.43237
[408,] -12535.91598 8039.37711
[409,] -51062.81077 -12535.91598
[410,] -27864.30386 -51062.81077
[411,] -217588.38755 -27864.30386
[412,] -1660.31473 -217588.38755
[413,] -109703.50562 -1660.31473
[414,] -66423.50149 -109703.50562
[415,] -61919.12149 -66423.50149
[416,] 160227.99049 -61919.12149
[417,] -17015.45784 160227.99049
[418,] -19046.11042 -17015.45784
[419,] -46818.19381 -19046.11042
[420,] -87356.28065 -46818.19381
[421,] -89470.05262 -87356.28065
[422,] -626069.48963 -89470.05262
[423,] -316767.42372 -626069.48963
[424,] -61948.31074 -316767.42372
[425,] -208492.56134 -61948.31074
[426,] -241012.72772 -208492.56134
[427,] -144726.09880 -241012.72772
[428,] -121166.99639 -144726.09880
[429,] 152350.80295 -121166.99639
[430,] -370254.75590 152350.80295
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2951183.24003 -58916.68341
2 547776.97946 2951183.24003
3 -592723.57478 547776.97946
4 -855582.68318 -592723.57478
5 -187303.82930 -855582.68318
6 1120158.08279 -187303.82930
7 422848.11045 1120158.08279
8 -703591.95988 422848.11045
9 613805.40397 -703591.95988
10 459456.04993 613805.40397
11 800071.51561 459456.04993
12 655359.62146 800071.51561
13 -468897.87437 655359.62146
14 -280661.11852 -468897.87437
15 565520.53288 -280661.11852
16 -259523.08584 565520.53288
17 -337447.59549 -259523.08584
18 -581028.38399 -337447.59549
19 52511.45770 -581028.38399
20 101447.04410 52511.45770
21 -391496.71174 101447.04410
22 120088.18252 -391496.71174
23 -127599.19151 120088.18252
24 -71606.42087 -127599.19151
25 326145.45395 -71606.42087
26 -639817.67921 326145.45395
27 -174375.28759 -639817.67921
28 621327.36936 -174375.28759
29 289543.00246 621327.36936
30 247663.51192 289543.00246
31 -168131.35366 247663.51192
32 151317.62041 -168131.35366
33 88593.24434 151317.62041
34 -311615.01065 88593.24434
35 52219.10887 -311615.01065
36 310262.12965 52219.10887
37 140406.31363 310262.12965
38 105590.98411 140406.31363
39 -133523.16907 105590.98411
40 -148996.18697 -133523.16907
41 30995.22637 -148996.18697
42 97726.72451 30995.22637
43 -53490.33926 97726.72451
44 61041.57451 -53490.33926
45 -141469.50111 61041.57451
46 47161.91860 -141469.50111
47 -196941.82015 47161.91860
48 -28247.02662 -196941.82015
49 -149590.54896 -28247.02662
50 92352.03082 -149590.54896
51 -43197.40228 92352.03082
52 -243825.49661 -43197.40228
53 75448.94087 -243825.49661
54 -283498.95034 75448.94087
55 114058.25801 -283498.95034
56 248553.73190 114058.25801
57 79226.58274 248553.73190
58 -80437.40445 79226.58274
59 -247557.75152 -80437.40445
60 -38989.03249 -247557.75152
61 -79310.46952 -38989.03249
62 -29908.81807 -79310.46952
63 -235797.07384 -29908.81807
64 62520.90448 -235797.07384
65 33957.12914 62520.90448
66 -151722.70630 33957.12914
67 -9987.33107 -151722.70630
68 -105144.12475 -9987.33107
69 -69820.92442 -105144.12475
70 -73823.45845 -69820.92442
71 166625.49732 -73823.45845
72 -145779.19403 166625.49732
73 -56387.81872 -145779.19403
74 -123106.25890 -56387.81872
75 -65954.22615 -123106.25890
76 -132389.61114 -65954.22615
77 -95381.48058 -132389.61114
78 20289.54469 -95381.48058
79 19983.58345 20289.54469
80 30180.16223 19983.58345
81 -58836.44905 30180.16223
82 -167652.48064 -58836.44905
83 50643.93943 -167652.48064
84 -48411.72652 50643.93943
85 -125132.00610 -48411.72652
86 59758.15692 -125132.00610
87 -80831.59816 59758.15692
88 -107434.97899 -80831.59816
89 -59345.23388 -107434.97899
90 -36690.09821 -59345.23388
91 -45261.40047 -36690.09821
92 110477.03904 -45261.40047
93 -256423.03336 110477.03904
94 127224.14916 -256423.03336
95 -4549.51391 127224.14916
96 -34777.47183 -4549.51391
97 24177.92868 -34777.47183
98 -107147.52429 24177.92868
99 -15983.74517 -107147.52429
100 -28627.86911 -15983.74517
101 -133068.99569 -28627.86911
102 -86388.00661 -133068.99569
103 -95303.81873 -86388.00661
104 35161.57075 -95303.81873
105 -85221.31954 35161.57075
106 -231658.23841 -85221.31954
107 99931.47560 -231658.23841
108 -98081.59614 99931.47560
109 -99669.94893 -98081.59614
110 78119.21937 -99669.94893
111 -243364.25034 78119.21937
112 -181459.97158 -243364.25034
113 -114009.58235 -181459.97158
114 36276.51649 -114009.58235
115 55746.57428 36276.51649
116 -41605.06999 55746.57428
117 53064.69920 -41605.06999
118 -170890.36224 53064.69920
119 17637.61731 -170890.36224
120 4368.30333 17637.61731
121 11382.81709 4368.30333
122 -52559.93239 11382.81709
123 68389.53214 -52559.93239
124 -92223.61851 68389.53214
125 -198868.43790 -92223.61851
126 269.66610 -198868.43790
127 492.78022 269.66610
128 9519.96729 492.78022
129 269.66610 9519.96729
130 876.92556 269.66610
131 81060.62993 876.92556
132 269.66610 81060.62993
133 -1683.22084 269.66610
134 269.66610 -1683.22084
135 9377.15855 269.66610
136 450.96100 9377.15855
137 269.66610 450.96100
138 269.66610 269.66610
139 269.66610 269.66610
140 5874.10042 269.66610
141 269.66610 5874.10042
142 21874.36284 269.66610
143 374.26243 21874.36284
144 15381.84798 374.26243
145 17301.29008 15381.84798
146 269.66610 17301.29008
147 -259.25801 269.66610
148 269.66610 -259.25801
149 -266.77164 269.66610
150 3795.40766 -266.77164
151 15573.66403 3795.40766
152 269.66610 15573.66403
153 -11688.27232 269.66610
154 -343.79348 -11688.27232
155 269.66610 -343.79348
156 269.66610 269.66610
157 -441.21905 269.66610
158 17959.78921 -441.21905
159 269.66610 17959.78921
160 269.66610 269.66610
161 269.66610 269.66610
162 269.66610 269.66610
163 9855.82870 269.66610
164 -7219.00178 9855.82870
165 269.66610 -7219.00178
166 269.66610 269.66610
167 119690.15922 269.66610
168 269.66610 119690.15922
169 269.66610 269.66610
170 269.66610 269.66610
171 -3306.88500 269.66610
172 269.66610 -3306.88500
173 818.00422 269.66610
174 269.66610 818.00422
175 -349.82655 269.66610
176 269.66610 -349.82655
177 50048.92944 269.66610
178 269.66610 50048.92944
179 269.66610 269.66610
180 32479.16452 269.66610
181 37849.42693 32479.16452
182 28406.06360 37849.42693
183 1979.83537 28406.06360
184 52212.88945 1979.83537
185 31771.85850 52212.88945
186 50634.13806 31771.85850
187 -222.76290 50634.13806
188 -33705.60342 -222.76290
189 269.66610 -33705.60342
190 12335.74127 269.66610
191 269.66610 12335.74127
192 4496.34472 269.66610
193 269.66610 4496.34472
194 26.76109 269.66610
195 52378.29308 26.76109
196 -55752.85014 52378.29308
197 52378.29308 -55752.85014
198 52378.29308 52378.29308
199 -6680.96187 52378.29308
200 52378.29308 -6680.96187
201 52378.29308 52378.29308
202 52378.29308 52378.29308
203 49347.13500 52378.29308
204 -1426.44177 49347.13500
205 51198.04438 -1426.44177
206 43133.28603 51198.04438
207 52378.29308 43133.28603
208 65274.78044 52378.29308
209 52378.29308 65274.78044
210 269.66610 52378.29308
211 269.66610 269.66610
212 10039.37072 269.66610
213 -15066.46219 10039.37072
214 -3416.17584 -15066.46219
215 269.66610 -3416.17584
216 600.33670 269.66610
217 4959.98750 600.33670
218 269.66610 4959.98750
219 31628.34903 269.66610
220 -24607.15405 31628.34903
221 269.66610 -24607.15405
222 343.94547 269.66610
223 269.66610 343.94547
224 333.61190 269.66610
225 2451.20280 333.61190
226 269.66610 2451.20280
227 269.66610 269.66610
228 -20900.79062 269.66610
229 269.66610 -20900.79062
230 -14030.37493 269.66610
231 269.66610 -14030.37493
232 269.66610 269.66610
233 6345.08101 269.66610
234 269.66610 6345.08101
235 132.21122 269.66610
236 269.66610 132.21122
237 2223.44248 269.66610
238 269.66610 2223.44248
239 269.66610 269.66610
240 269.66610 269.66610
241 5004.31516 269.66610
242 269.66610 5004.31516
243 2504.79811 269.66610
244 -8408.43581 2504.79811
245 52378.29308 -8408.43581
246 52378.29308 52378.29308
247 -8431.64237 52378.29308
248 52378.29308 -8431.64237
249 52378.29308 52378.29308
250 9587.91386 52378.29308
251 52378.29308 9587.91386
252 -18115.78580 52378.29308
253 -31116.13407 -18115.78580
254 -58439.63706 -31116.13407
255 269.66610 -58439.63706
256 31888.00688 269.66610
257 35271.83324 31888.00688
258 52378.29308 35271.83324
259 269.66610 52378.29308
260 52378.29308 269.66610
261 52378.29308 52378.29308
262 -16366.47699 52378.29308
263 12167.65305 -16366.47699
264 52378.29308 12167.65305
265 269.66610 52378.29308
266 52378.29308 269.66610
267 6138.66863 52378.29308
268 52378.29308 6138.66863
269 -18679.71641 52378.29308
270 52378.29308 -18679.71641
271 52378.29308 52378.29308
272 52378.29308 52378.29308
273 -2525.19903 52378.29308
274 269.66610 -2525.19903
275 5208.34419 269.66610
276 3639.05382 5208.34419
277 269.66610 3639.05382
278 6702.44125 269.66610
279 -27192.11395 6702.44125
280 52378.29308 -27192.11395
281 269.66610 52378.29308
282 13641.30379 269.66610
283 269.66610 13641.30379
284 -29547.57748 269.66610
285 269.66610 -29547.57748
286 52378.29308 269.66610
287 -1440.77572 52378.29308
288 52378.29308 -1440.77572
289 10112.41764 52378.29308
290 52378.29308 10112.41764
291 4671.35207 52378.29308
292 -19267.78540 4671.35207
293 269.66610 -19267.78540
294 -884.04067 269.66610
295 71481.29254 -884.04067
296 269.66610 71481.29254
297 52378.29308 269.66610
298 269.66610 52378.29308
299 269.66610 269.66610
300 49654.99785 269.66610
301 -4999.40151 49654.99785
302 269.66610 -4999.40151
303 52378.29308 269.66610
304 27201.04130 52378.29308
305 263.60403 27201.04130
306 52378.29308 263.60403
307 41978.19496 52378.29308
308 24444.45331 41978.19496
309 52378.29308 24444.45331
310 52378.29308 52378.29308
311 -43764.97916 52378.29308
312 52378.29308 -43764.97916
313 50952.03341 52378.29308
314 463.22713 50952.03341
315 -152.93183 463.22713
316 2427.89874 -152.93183
317 -9173.81075 2427.89874
318 277.73992 -9173.81075
319 269.66610 277.73992
320 269.66610 269.66610
321 269.66610 269.66610
322 269.66610 269.66610
323 52378.29308 269.66610
324 269.66610 52378.29308
325 -2479.88831 269.66610
326 52378.29308 -2479.88831
327 52378.29308 52378.29308
328 269.66610 52378.29308
329 269.66610 269.66610
330 269.66610 269.66610
331 269.66610 269.66610
332 269.66610 269.66610
333 1209.65733 269.66610
334 260.34712 1209.65733
335 269.66610 260.34712
336 -1367.64739 269.66610
337 269.66610 -1367.64739
338 53.91705 269.66610
339 2427.04406 53.91705
340 -15256.62205 2427.04406
341 269.66610 -15256.62205
342 269.66610 269.66610
343 269.66610 269.66610
344 269.66610 269.66610
345 269.66610 269.66610
346 -6619.60825 269.66610
347 269.66610 -6619.60825
348 269.66610 269.66610
349 269.66610 269.66610
350 4938.31961 269.66610
351 -968.68749 4938.31961
352 16013.16365 -968.68749
353 -5853.03737 16013.16365
354 269.66610 -5853.03737
355 269.66610 269.66610
356 -37944.63660 269.66610
357 269.66610 -37944.63660
358 269.66610 269.66610
359 269.66610 269.66610
360 269.46762 269.66610
361 269.66610 269.46762
362 269.66610 269.66610
363 42791.96245 269.66610
364 -247.23550 42791.96245
365 269.66610 -247.23550
366 269.66610 269.66610
367 269.66610 269.66610
368 269.66610 269.66610
369 269.66610 269.66610
370 269.66610 269.66610
371 -129.02984 269.66610
372 -34857.22662 -129.02984
373 -20764.09565 -34857.22662
374 269.66610 -20764.09565
375 269.66610 269.66610
376 5212.54708 269.66610
377 574.48936 5212.54708
378 269.66610 574.48936
379 269.66610 269.66610
380 269.66610 269.66610
381 -15401.57498 269.66610
382 34164.62363 -15401.57498
383 269.66610 34164.62363
384 -2228.37519 269.66610
385 -89637.51538 -2228.37519
386 6519.90171 -89637.51538
387 11836.96694 6519.90171
388 -519.99381 11836.96694
389 -7236.11997 -519.99381
390 6069.64377 -7236.11997
391 -8800.34051 6069.64377
392 -19588.75261 -8800.34051
393 37294.67149 -19588.75261
394 -157469.05530 37294.67149
395 -14508.02866 -157469.05530
396 -16945.45973 -14508.02866
397 98646.88263 -16945.45973
398 -12258.55024 98646.88263
399 -105526.30134 -12258.55024
400 -3108.63911 -105526.30134
401 12572.25382 -3108.63911
402 -3503.48237 12572.25382
403 -40978.81913 -3503.48237
404 -16209.25828 -40978.81913
405 -16937.34734 -16209.25828
406 -1202.43237 -16937.34734
407 8039.37711 -1202.43237
408 -12535.91598 8039.37711
409 -51062.81077 -12535.91598
410 -27864.30386 -51062.81077
411 -217588.38755 -27864.30386
412 -1660.31473 -217588.38755
413 -109703.50562 -1660.31473
414 -66423.50149 -109703.50562
415 -61919.12149 -66423.50149
416 160227.99049 -61919.12149
417 -17015.45784 160227.99049
418 -19046.11042 -17015.45784
419 -46818.19381 -19046.11042
420 -87356.28065 -46818.19381
421 -89470.05262 -87356.28065
422 -626069.48963 -89470.05262
423 -316767.42372 -626069.48963
424 -61948.31074 -316767.42372
425 -208492.56134 -61948.31074
426 -241012.72772 -208492.56134
427 -144726.09880 -241012.72772
428 -121166.99639 -144726.09880
429 152350.80295 -121166.99639
430 -370254.75590 152350.80295
> 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/75izs1291335824.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/8f9gv1291335824.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/9f9gv1291335824.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/10f9gv1291335824.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/11t1e41291335824.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/12eky11291335825.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/133lvd1291335825.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/14ecug1291335825.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/15zdbm1291335825.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/16dn8d1291335825.tab")
+ }
>
> try(system("convert tmp/1jzi41291335824.ps tmp/1jzi41291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/2jzi41291335824.ps tmp/2jzi41291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/3jzi41291335824.ps tmp/3jzi41291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/4crz71291335824.ps tmp/4crz71291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/5crz71291335824.ps tmp/5crz71291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/6crz71291335824.ps tmp/6crz71291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/75izs1291335824.ps tmp/75izs1291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/8f9gv1291335824.ps tmp/8f9gv1291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/9f9gv1291335824.ps tmp/9f9gv1291335824.png",intern=TRUE))
character(0)
> try(system("convert tmp/10f9gv1291335824.ps tmp/10f9gv1291335824.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
12.810 2.196 26.236