R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1
+ ,162556
+ ,162556
+ ,213118
+ ,213118
+ ,230380558
+ ,6282929
+ ,1
+ ,29790
+ ,29790
+ ,81767
+ ,81767
+ ,25266003
+ ,4324047
+ ,1
+ ,87550
+ ,87550
+ ,153198
+ ,153198
+ ,70164684
+ ,4108272
+ ,0
+ ,84738
+ ,0
+ ,-26007
+ ,0
+ ,-15292116
+ ,-1212617
+ ,1
+ ,54660
+ ,54660
+ ,126942
+ ,126942
+ ,37955658
+ ,1485329
+ ,1
+ ,42634
+ ,42634
+ ,157214
+ ,157214
+ ,24525384
+ ,1779876
+ ,0
+ ,40949
+ ,0
+ ,129352
+ ,0
+ ,62218312
+ ,1367203
+ ,1
+ ,42312
+ ,42312
+ ,234817
+ ,234817
+ ,75845891
+ ,2519076
+ ,1
+ ,37704
+ ,37704
+ ,60448
+ ,60448
+ ,27322496
+ ,912684
+ ,1
+ ,16275
+ ,16275
+ ,47818
+ ,47818
+ ,5212162
+ ,1443586
+ ,0
+ ,25830
+ ,0
+ ,245546
+ ,0
+ ,28237790
+ ,1220017
+ ,0
+ ,12679
+ ,0
+ ,48020
+ ,0
+ ,5282200
+ ,984885
+ ,1
+ ,18014
+ ,18014
+ ,-1710
+ ,-1710
+ ,-408690
+ ,1457425
+ ,0
+ ,43556
+ ,0
+ ,32648
+ ,0
+ ,8064056
+ ,-572920
+ ,1
+ ,24524
+ ,24524
+ ,95350
+ ,95350
+ ,47388950
+ ,929144
+ ,0
+ ,6532
+ ,0
+ ,151352
+ ,0
+ ,15589256
+ ,1151176
+ ,0
+ ,7123
+ ,0
+ ,288170
+ ,0
+ ,31410530
+ ,790090
+ ,1
+ ,20813
+ ,20813
+ ,114337
+ ,114337
+ ,57397174
+ ,774497
+ ,1
+ ,37597
+ ,37597
+ ,37884
+ ,37884
+ ,9395232
+ ,990576
+ ,0
+ ,17821
+ ,0
+ ,122844
+ ,0
+ ,45820812
+ ,454195
+ ,1
+ ,12988
+ ,12988
+ ,82340
+ ,82340
+ ,9798460
+ ,876607
+ ,1
+ ,22330
+ ,22330
+ ,79801
+ ,79801
+ ,6703284
+ ,711969
+ ,0
+ ,13326
+ ,0
+ ,165548
+ ,0
+ ,16885896
+ ,702380
+ ,0
+ ,16189
+ ,0
+ ,116384
+ ,0
+ ,34333280
+ ,264449
+ ,0
+ ,7146
+ ,0
+ ,134028
+ ,0
+ ,14072940
+ ,450033
+ ,0
+ ,15824
+ ,0
+ ,63838
+ ,0
+ ,4085632
+ ,541063
+ ,1
+ ,26088
+ ,26088
+ ,74996
+ ,74996
+ ,20023932
+ ,588864
+ ,0
+ ,11326
+ ,0
+ ,31080
+ ,0
+ ,4009320
+ ,-37216
+ ,0
+ ,8568
+ ,0
+ ,32168
+ ,0
+ ,1190216
+ ,783310
+ ,0
+ ,14416
+ ,0
+ ,49857
+ ,0
+ ,17998377
+ ,467359
+ ,1
+ ,3369
+ ,3369
+ ,87161
+ ,87161
+ ,2440508
+ ,688779
+ ,1
+ ,11819
+ ,11819
+ ,106113
+ ,106113
+ ,9019605
+ ,608419
+ ,1
+ ,6620
+ ,6620
+ ,80570
+ ,80570
+ ,3545080
+ ,696348
+ ,1
+ ,4519
+ ,4519
+ ,102129
+ ,102129
+ ,5004321
+ ,597793
+ ,0
+ ,2220
+ ,0
+ ,301670
+ ,0
+ ,6636740
+ ,821730
+ ,0
+ ,18562
+ ,0
+ ,102313
+ ,0
+ ,15858515
+ ,377934
+ ,0
+ ,10327
+ ,0
+ ,88577
+ ,0
+ ,8060507
+ ,651939
+ ,1
+ ,5336
+ ,5336
+ ,112477
+ ,112477
+ ,9110637
+ ,697458
+ ,1
+ ,2365
+ ,2365
+ ,191778
+ ,191778
+ ,15150462
+ ,700368
+ ,0
+ ,4069
+ ,0
+ ,79804
+ ,0
+ ,11571580
+ ,225986
+ ,0
+ ,7710
+ ,0
+ ,128294
+ ,0
+ ,104687904
+ ,348695
+ ,0
+ ,13718
+ ,0
+ ,96448
+ ,0
+ ,5883328
+ ,373683
+ ,0
+ ,4525
+ ,0
+ ,93811
+ ,0
+ ,21201286
+ ,501709
+ ,0
+ ,6869
+ ,0
+ ,117520
+ ,0
+ ,12339600
+ ,413743
+ ,0
+ ,4628
+ ,0
+ ,69159
+ ,0
+ ,4287858
+ ,379825
+ ,1
+ ,3653
+ ,3653
+ ,101792
+ ,101792
+ ,2443008
+ ,336260
+ ,1
+ ,1265
+ ,1265
+ ,210568
+ ,210568
+ ,5474768
+ ,636765
+ ,1
+ ,7489
+ ,7489
+ ,136996
+ ,136996
+ ,44112712
+ ,481231
+ ,0
+ ,4901
+ ,0
+ ,121920
+ ,0
+ ,10241280
+ ,469107
+ ,0
+ ,2284
+ ,0
+ ,76403
+ ,0
+ ,2521299
+ ,211928
+ ,1
+ ,3160
+ ,3160
+ ,108094
+ ,108094
+ ,11674152
+ ,563925
+ ,1
+ ,4150
+ ,4150
+ ,134759
+ ,134759
+ ,20213850
+ ,511939
+ ,1
+ ,7285
+ ,7285
+ ,188873
+ ,188873
+ ,21720395
+ ,521016
+ ,1
+ ,1134
+ ,1134
+ ,146216
+ ,146216
+ ,23686992
+ ,543856
+ ,1
+ ,4658
+ ,4658
+ ,156608
+ ,156608
+ ,24744064
+ ,329304
+ ,0
+ ,2384
+ ,0
+ ,61348
+ ,0
+ ,5950756
+ ,423262
+ ,0
+ ,3748
+ ,0
+ ,50350
+ ,0
+ ,453150
+ ,509665
+ ,0
+ ,5371
+ ,0
+ ,87720
+ ,0
+ ,5789520
+ ,455881
+ ,0
+ ,1285
+ ,0
+ ,99489
+ ,0
+ ,10645323
+ ,367772
+ ,1
+ ,9327
+ ,9327
+ ,87419
+ ,87419
+ ,8829319
+ ,406339
+ ,1
+ ,5565
+ ,5565
+ ,94355
+ ,94355
+ ,4434685
+ ,493408
+ ,0
+ ,1528
+ ,0
+ ,60326
+ ,0
+ ,2292388
+ ,232942
+ ,1
+ ,3122
+ ,3122
+ ,94670
+ ,94670
+ ,3218780
+ ,416002
+ ,1
+ ,7317
+ ,7317
+ ,82425
+ ,82425
+ ,6923700
+ ,337430
+ ,0
+ ,2675
+ ,0
+ ,59017
+ ,0
+ ,4662343
+ ,361517
+ ,0
+ ,13253
+ ,0
+ ,90829
+ ,0
+ ,86015063
+ ,360962
+ ,0
+ ,880
+ ,0
+ ,80791
+ ,0
+ ,5978534
+ ,235561
+ ,1
+ ,2053
+ ,2053
+ ,100423
+ ,100423
+ ,5322419
+ ,408247
+ ,0
+ ,1424
+ ,0
+ ,131116
+ ,0
+ ,12324904
+ ,450296
+ ,1
+ ,4036
+ ,4036
+ ,100269
+ ,100269
+ ,6316947
+ ,418799
+ ,1
+ ,3045
+ ,3045
+ ,27330
+ ,27330
+ ,1585140
+ ,247405
+ ,0
+ ,5119
+ ,0
+ ,39039
+ ,0
+ ,1912911
+ ,378519
+ ,0
+ ,1431
+ ,0
+ ,106885
+ ,0
+ ,3634090
+ ,326638
+ ,0
+ ,554
+ ,0
+ ,79285
+ ,0
+ ,872135
+ ,328233
+ ,0
+ ,1975
+ ,0
+ ,118881
+ ,0
+ ,4160835
+ ,386225
+ ,1
+ ,1286
+ ,1286
+ ,77623
+ ,77623
+ ,1319591
+ ,283662
+ ,0
+ ,1012
+ ,0
+ ,114768
+ ,0
+ ,5394096
+ ,370225
+ ,0
+ ,810
+ ,0
+ ,74015
+ ,0
+ ,3182645
+ ,269236
+ ,0
+ ,1280
+ ,0
+ ,69465
+ ,0
+ ,8127405
+ ,365732
+ ,1
+ ,666
+ ,666
+ ,117869
+ ,117869
+ ,20155599
+ ,420383
+ ,0
+ ,1380
+ ,0
+ ,60982
+ ,0
+ ,1585532
+ ,345811
+ ,1
+ ,4608
+ ,4608
+ ,90131
+ ,90131
+ ,6579563
+ ,431809
+ ,0
+ ,876
+ ,0
+ ,138971
+ ,0
+ ,8199289
+ ,418876
+ ,0
+ ,814
+ ,0
+ ,39625
+ ,0
+ ,713250
+ ,297476
+ ,0
+ ,514
+ ,0
+ ,102725
+ ,0
+ ,1540875
+ ,416776
+ ,1
+ ,5692
+ ,5692
+ ,64239
+ ,64239
+ ,4625208
+ ,357257
+ ,0
+ ,3642
+ ,0
+ ,90262
+ ,0
+ ,7762532
+ ,458343
+ ,0
+ ,540
+ ,0
+ ,103960
+ ,0
+ ,1455440
+ ,388386
+ ,0
+ ,2099
+ ,0
+ ,106611
+ ,0
+ ,6823104
+ ,358934
+ ,0
+ ,567
+ ,0
+ ,103345
+ ,0
+ ,1136795
+ ,407560
+ ,0
+ ,2001
+ ,0
+ ,95551
+ ,0
+ ,4968652
+ ,392558
+ ,1
+ ,2949
+ ,2949
+ ,82903
+ ,82903
+ ,3399023
+ ,373177
+ ,0
+ ,2253
+ ,0
+ ,63593
+ ,0
+ ,6295707
+ ,428370
+ ,1
+ ,6533
+ ,6533
+ ,126910
+ ,126910
+ ,9518250
+ ,369419
+ ,0
+ ,1889
+ ,0
+ ,37527
+ ,0
+ ,1688715
+ ,358649
+ ,1
+ ,3055
+ ,3055
+ ,60247
+ ,60247
+ ,2590621
+ ,376641
+ ,0
+ ,272
+ ,0
+ ,112995
+ ,0
+ ,903960
+ ,467427
+ ,1
+ ,1414
+ ,1414
+ ,70184
+ ,70184
+ ,13896432
+ ,364885
+ ,0
+ ,2564
+ ,0
+ ,130140
+ ,0
+ ,2863080
+ ,436230
+ ,1
+ ,1383
+ ,1383
+ ,73221
+ ,73221
+ ,805431
+ ,329118
+ ,1
+ ,1261
+ ,1261
+ ,76114
+ ,76114
+ ,2511762
+ ,317365
+ ,0
+ ,975
+ ,0
+ ,90534
+ ,0
+ ,2082282
+ ,286849
+ ,0
+ ,3366
+ ,0
+ ,108479
+ ,0
+ ,8678320
+ ,376685
+ ,0
+ ,576
+ ,0
+ ,113761
+ ,0
+ ,2047698
+ ,407198
+ ,0
+ ,1306
+ ,0
+ ,68696
+ ,0
+ ,1923488
+ ,377772
+ ,0
+ ,746
+ ,0
+ ,71561
+ ,0
+ ,1645903
+ ,271483
+ ,1
+ ,3192
+ ,3192
+ ,59831
+ ,59831
+ ,3589860
+ ,153661
+ ,1
+ ,2045
+ ,2045
+ ,97890
+ ,97890
+ ,1957800
+ ,513294
+ ,0
+ ,5477
+ ,0
+ ,101481
+ ,0
+ ,5987379
+ ,324881
+ ,1
+ ,1932
+ ,1932
+ ,72954
+ ,72954
+ ,2626344
+ ,264512
+ ,0
+ ,936
+ ,0
+ ,67939
+ ,0
+ ,2038170
+ ,420968
+ ,1
+ ,3437
+ ,3437
+ ,48022
+ ,48022
+ ,2257034
+ ,129302
+ ,0
+ ,5131
+ ,0
+ ,86111
+ ,0
+ ,6113881
+ ,191521
+ ,1
+ ,2397
+ ,2397
+ ,74020
+ ,74020
+ ,1036280
+ ,268673
+ ,1
+ ,1389
+ ,1389
+ ,57530
+ ,57530
+ ,517770
+ ,353179
+ ,0
+ ,1503
+ ,0
+ ,56364
+ ,0
+ ,2198196
+ ,354624
+ ,0
+ ,402
+ ,0
+ ,84990
+ ,0
+ ,2209740
+ ,363713
+ ,0
+ ,2239
+ ,0
+ ,88590
+ ,0
+ ,1860390
+ ,456657
+ ,1
+ ,2234
+ ,2234
+ ,77200
+ ,77200
+ ,1235200
+ ,211742
+ ,0
+ ,837
+ ,0
+ ,61262
+ ,0
+ ,4227078
+ ,338381
+ ,0
+ ,10579
+ ,0
+ ,110309
+ ,0
+ ,10148428
+ ,418530
+ ,0
+ ,875
+ ,0
+ ,67000
+ ,0
+ ,938000
+ ,351483
+ ,0
+ ,1395
+ ,0
+ ,93099
+ ,0
+ ,9589197
+ ,372928
+ ,1
+ ,1659
+ ,1659
+ ,107577
+ ,107577
+ ,3119733
+ ,485538
+ ,1
+ ,2647
+ ,2647
+ ,62920
+ ,62920
+ ,2328040
+ ,279268
+ ,1
+ ,3294
+ ,3294
+ ,75832
+ ,75832
+ ,1744136
+ ,219060
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,0
+ ,60793
+ ,0
+ ,425551
+ ,325314
+ ,0
+ ,422
+ ,0
+ ,57935
+ ,0
+ ,1622180
+ ,322046
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,34
+ ,0
+ ,60630
+ ,0
+ ,485040
+ ,325599
+ ,0
+ ,1558
+ ,0
+ ,55637
+ ,0
+ ,3505131
+ ,377028
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,43
+ ,0
+ ,60887
+ ,0
+ ,182661
+ ,323850
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,316
+ ,0
+ ,60505
+ ,0
+ ,544545
+ ,331514
+ ,0
+ ,115
+ ,0
+ ,60945
+ ,0
+ ,792285
+ ,325632
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,389
+ ,0
+ ,58990
+ ,0
+ ,825860
+ ,322265
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,1002
+ ,0
+ ,56750
+ ,0
+ ,851250
+ ,325906
+ ,0
+ ,36
+ ,0
+ ,60894
+ ,0
+ ,182682
+ ,325985
+ ,0
+ ,460
+ ,0
+ ,63346
+ ,0
+ ,950190
+ ,346145
+ ,0
+ ,309
+ ,0
+ ,56535
+ ,0
+ ,621885
+ ,325898
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,9
+ ,0
+ ,60835
+ ,0
+ ,365010
+ ,325356
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,14
+ ,0
+ ,61016
+ ,0
+ ,61016
+ ,325930
+ ,0
+ ,520
+ ,0
+ ,58650
+ ,0
+ ,586500
+ ,318020
+ ,0
+ ,1766
+ ,0
+ ,60438
+ ,0
+ ,4411974
+ ,326389
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,458
+ ,0
+ ,58625
+ ,0
+ ,644875
+ ,302925
+ ,0
+ ,20
+ ,0
+ ,60938
+ ,0
+ ,182814
+ ,325540
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,98
+ ,0
+ ,61490
+ ,0
+ ,122980
+ ,326736
+ ,0
+ ,405
+ ,0
+ ,60845
+ ,0
+ ,425915
+ ,340580
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,483
+ ,0
+ ,60830
+ ,0
+ ,1642410
+ ,331828
+ ,0
+ ,454
+ ,0
+ ,63261
+ ,0
+ ,3226311
+ ,323299
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,757
+ ,0
+ ,45689
+ ,0
+ ,868091
+ ,387722
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,36
+ ,0
+ ,61564
+ ,0
+ ,246256
+ ,324598
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,203
+ ,0
+ ,61938
+ ,0
+ ,557442
+ ,328726
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,90
+ ,0
+ ,60951
+ ,0
+ ,487608
+ ,325043
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,71
+ ,71
+ ,60745
+ ,60745
+ ,60745
+ ,325806
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,972
+ ,0
+ ,71642
+ ,0
+ ,2435828
+ ,387732
+ ,1
+ ,531
+ ,531
+ ,71641
+ ,71641
+ ,716410
+ ,349729
+ ,0
+ ,604
+ ,0
+ ,55792
+ ,0
+ ,2120096
+ ,332202
+ ,1
+ ,283
+ ,283
+ ,71873
+ ,71873
+ ,718730
+ ,305442
+ ,1
+ ,23
+ ,23
+ ,62555
+ ,62555
+ ,312775
+ ,329537
+ ,1
+ ,638
+ ,638
+ ,60370
+ ,60370
+ ,845180
+ ,327055
+ ,1
+ ,699
+ ,699
+ ,64873
+ ,64873
+ ,1037968
+ ,356245
+ ,0
+ ,149
+ ,0
+ ,62041
+ ,0
+ ,310205
+ ,328451
+ ,0
+ ,226
+ ,0
+ ,65745
+ ,0
+ ,328725
+ ,307062
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,275
+ ,0
+ ,59500
+ ,0
+ ,238000
+ ,331345
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,141
+ ,0
+ ,61630
+ ,0
+ ,369780
+ ,331824
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,28
+ ,0
+ ,60890
+ ,0
+ ,121780
+ ,325685
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,2566
+ ,2566
+ ,113521
+ ,113521
+ ,9081680
+ ,404480
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,472
+ ,472
+ ,80045
+ ,80045
+ ,1600900
+ ,318314
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,203
+ ,203
+ ,50804
+ ,50804
+ ,1371708
+ ,311807
+ ,1
+ ,496
+ ,496
+ ,87390
+ ,87390
+ ,1485630
+ ,337724
+ ,1
+ ,10
+ ,10
+ ,61656
+ ,61656
+ ,123312
+ ,326431
+ ,1
+ ,63
+ ,63
+ ,65688
+ ,65688
+ ,262752
+ ,327556
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,1136
+ ,1136
+ ,48522
+ ,48522
+ ,1552704
+ ,356850
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,267
+ ,0
+ ,57640
+ ,0
+ ,1844480
+ ,322741
+ ,0
+ ,474
+ ,0
+ ,61977
+ ,0
+ ,1239540
+ ,310902
+ ,0
+ ,534
+ ,0
+ ,62620
+ ,0
+ ,438340
+ ,324295
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,15
+ ,0
+ ,60831
+ ,0
+ ,486648
+ ,326156
+ ,0
+ ,397
+ ,0
+ ,60646
+ ,0
+ ,1698088
+ ,326960
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,1061
+ ,0
+ ,56225
+ ,0
+ ,1124500
+ ,333411
+ ,0
+ ,288
+ ,0
+ ,60510
+ ,0
+ ,242040
+ ,297761
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,3
+ ,0
+ ,60698
+ ,0
+ ,121396
+ ,325536
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,20
+ ,0
+ ,60805
+ ,0
+ ,121610
+ ,325762
+ ,0
+ ,278
+ ,0
+ ,61404
+ ,0
+ ,1596504
+ ,327957
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,192
+ ,0
+ ,65276
+ ,0
+ ,261104
+ ,318521
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,317
+ ,0
+ ,63915
+ ,0
+ ,575235
+ ,319775
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,368
+ ,0
+ ,61686
+ ,0
+ ,1048662
+ ,332128
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,2
+ ,0
+ ,60743
+ ,0
+ ,60743
+ ,325486
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,0
+ ,60349
+ ,0
+ ,362094
+ ,325838
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,94
+ ,0
+ ,61360
+ ,0
+ ,184080
+ ,331767
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,24
+ ,0
+ ,59818
+ ,0
+ ,478544
+ ,324523
+ ,0
+ ,2332
+ ,0
+ ,72680
+ ,0
+ ,290720
+ ,339995
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,131
+ ,0
+ ,61808
+ ,0
+ ,679888
+ ,319582
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,206
+ ,0
+ ,53110
+ ,0
+ ,477990
+ ,307245
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,167
+ ,0
+ ,64245
+ ,0
+ ,128490
+ ,317967
+ ,0
+ ,622
+ ,0
+ ,73007
+ ,0
+ ,5329511
+ ,331488
+ ,0
+ ,885
+ ,0
+ ,82732
+ ,0
+ ,7032220
+ ,335452
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,365
+ ,0
+ ,54820
+ ,0
+ ,438560
+ ,334184
+ ,0
+ ,364
+ ,0
+ ,47705
+ ,0
+ ,1669675
+ ,313213
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,226
+ ,0
+ ,72835
+ ,0
+ ,874020
+ ,348678
+ ,0
+ ,307
+ ,0
+ ,58856
+ ,0
+ ,882840
+ ,328727
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,188
+ ,0
+ ,77655
+ ,0
+ ,854205
+ ,387978
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,138
+ ,0
+ ,69817
+ ,0
+ ,418902
+ ,336704
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,125
+ ,0
+ ,60798
+ ,0
+ ,729576
+ ,322076
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,282
+ ,0
+ ,62452
+ ,0
+ ,1873560
+ ,334272
+ ,0
+ ,335
+ ,0
+ ,64175
+ ,0
+ ,2117775
+ ,338197
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,813
+ ,813
+ ,67440
+ ,67440
+ ,5530080
+ ,321024
+ ,0
+ ,176
+ ,0
+ ,68136
+ ,0
+ ,1907808
+ ,322145
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,249
+ ,0
+ ,56726
+ ,0
+ ,4084272
+ ,323351
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,333
+ ,0
+ ,70811
+ ,0
+ ,920543
+ ,327748
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,30
+ ,0
+ ,62045
+ ,0
+ ,248180
+ ,328157
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,249
+ ,0
+ ,54323
+ ,0
+ ,3368026
+ ,311594
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,165
+ ,0
+ ,62841
+ ,0
+ ,1508184
+ ,335962
+ ,0
+ ,453
+ ,0
+ ,81125
+ ,0
+ ,1703625
+ ,372426
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,53
+ ,0
+ ,59506
+ ,0
+ ,833084
+ ,319844
+ ,1
+ ,382
+ ,382
+ ,59365
+ ,59365
+ ,1246665
+ ,355822
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,30
+ ,30
+ ,60798
+ ,60798
+ ,243192
+ ,324047
+ ,0
+ ,290
+ ,0
+ ,58790
+ ,0
+ ,117580
+ ,311464
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,366
+ ,0
+ ,61808
+ ,0
+ ,3275824
+ ,353417
+ ,0
+ ,2
+ ,0
+ ,60735
+ ,0
+ ,546615
+ ,325590
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,209
+ ,209
+ ,64016
+ ,64016
+ ,832208
+ ,328576
+ ,0
+ ,384
+ ,0
+ ,54683
+ ,0
+ ,1203026
+ ,326126
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,365
+ ,0
+ ,87192
+ ,0
+ ,7236936
+ ,369376
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,49
+ ,49
+ ,64107
+ ,64107
+ ,512856
+ ,332013
+ ,0
+ ,3
+ ,0
+ ,60761
+ ,0
+ ,243044
+ ,325871
+ ,0
+ ,133
+ ,0
+ ,65990
+ ,0
+ ,923860
+ ,342165
+ ,0
+ ,32
+ ,0
+ ,59988
+ ,0
+ ,59988
+ ,324967
+ ,0
+ ,368
+ ,0
+ ,61167
+ ,0
+ ,1039839
+ ,314832
+ ,0
+ ,1
+ ,0
+ ,60719
+ ,0
+ ,364314
+ ,325557
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,22
+ ,0
+ ,60722
+ ,0
+ ,121444
+ ,322649
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,60720
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,96
+ ,0
+ ,60379
+ ,0
+ ,301895
+ ,324598
+ ,0
+ ,1
+ ,0
+ ,60727
+ ,0
+ ,121454
+ ,325567
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,81
+ ,0
+ ,60925
+ ,0
+ ,426475
+ ,324005
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,26
+ ,0
+ ,60896
+ ,0
+ ,60896
+ ,325748
+ ,0
+ ,125
+ ,0
+ ,59734
+ ,0
+ ,776542
+ ,323385
+ ,0
+ ,304
+ ,0
+ ,62969
+ ,0
+ ,944535
+ ,315409
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,119
+ ,0
+ ,59118
+ ,0
+ ,354708
+ ,312275
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,312
+ ,0
+ ,58598
+ ,0
+ ,820372
+ ,320576
+ ,0
+ ,60
+ ,0
+ ,61124
+ ,0
+ ,611240
+ ,325246
+ ,0
+ ,587
+ ,0
+ ,59595
+ ,0
+ ,715140
+ ,332961
+ ,0
+ ,135
+ ,0
+ ,62065
+ ,0
+ ,124130
+ ,323010
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,514
+ ,0
+ ,78780
+ ,0
+ ,4096560
+ ,345253
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,1
+ ,0
+ ,60722
+ ,0
+ ,242888
+ ,325559
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,1
+ ,58
+ ,58
+ ,61600
+ ,61600
+ ,184800
+ ,319634
+ ,0
+ ,180
+ ,0
+ ,59635
+ ,0
+ ,655985
+ ,319951
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,448
+ ,0
+ ,59781
+ ,0
+ ,2391240
+ ,318519
+ ,0
+ ,227
+ ,0
+ ,76644
+ ,0
+ ,689796
+ ,343222
+ ,0
+ ,174
+ ,0
+ ,64820
+ ,0
+ ,64820
+ ,317234
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,121
+ ,0
+ ,56178
+ ,0
+ ,1348272
+ ,314025
+ ,0
+ ,607
+ ,0
+ ,60436
+ ,0
+ ,664796
+ ,320249
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,530
+ ,0
+ ,73433
+ ,0
+ ,4405980
+ ,349365
+ ,0
+ ,571
+ ,0
+ ,41477
+ ,0
+ ,3318160
+ ,289197
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,0
+ ,0
+ ,325560
+ ,0
+ ,78
+ ,0
+ ,62700
+ ,0
+ ,1003200
+ ,329245
+ ,0
+ ,2489
+ ,0
+ ,67804
+ ,0
+ ,2712160
+ ,240869
+ ,0
+ ,131
+ ,0
+ ,59661
+ ,0
+ ,357966
+ ,327182
+ ,0
+ ,923
+ ,0
+ ,58620
+ ,0
+ ,468960
+ ,322876
+ ,0
+ ,72
+ ,0
+ ,60398
+ ,0
+ ,181194
+ ,323117
+ ,0
+ ,572
+ ,0
+ ,58580
+ ,0
+ ,937280
+ ,306351
+ ,0
+ ,397
+ ,0
+ ,62710
+ ,0
+ ,627100
+ ,335137
+ ,0
+ ,450
+ ,0
+ ,59325
+ ,0
+ ,474600
+ ,308271
+ ,0
+ ,622
+ ,0
+ ,60950
+ ,0
+ ,426650
+ ,301731
+ ,0
+ ,694
+ ,0
+ ,68060
+ ,0
+ ,544480
+ ,382409
+ ,1
+ ,3425
+ ,3425
+ ,83620
+ ,83620
+ ,1003440
+ ,279230
+ ,0
+ ,562
+ ,0
+ ,58456
+ ,0
+ ,759928
+ ,298731
+ ,0
+ ,4917
+ ,0
+ ,52811
+ ,0
+ ,2218062
+ ,243650
+ ,1
+ ,1442
+ ,1442
+ ,121173
+ ,121173
+ ,14298414
+ ,532682
+ ,0
+ ,529
+ ,0
+ ,63870
+ ,0
+ ,574830
+ ,319771
+ ,1
+ ,2126
+ ,2126
+ ,21001
+ ,21001
+ ,2898138
+ ,171493
+ ,0
+ ,1061
+ ,0
+ ,70415
+ ,0
+ ,352075
+ ,347262
+ ,0
+ ,776
+ ,0
+ ,64230
+ ,0
+ ,578070
+ ,343945
+ ,0
+ ,611
+ ,0
+ ,59190
+ ,0
+ ,473520
+ ,311874
+ ,1
+ ,1526
+ ,1526
+ ,69351
+ ,69351
+ ,1733775
+ ,302211
+ ,0
+ ,592
+ ,0
+ ,64270
+ ,0
+ ,449890
+ ,316708
+ ,0
+ ,1182
+ ,0
+ ,70694
+ ,0
+ ,919022
+ ,333463
+ ,0
+ ,621
+ ,0
+ ,68005
+ ,0
+ ,1088080
+ ,344282
+ ,0
+ ,989
+ ,0
+ ,58930
+ ,0
+ ,648230
+ ,319635
+ ,0
+ ,438
+ ,0
+ ,58320
+ ,0
+ ,641520
+ ,301186
+ ,0
+ ,726
+ ,0
+ ,69980
+ ,0
+ ,209940
+ ,300381
+ ,0
+ ,1303
+ ,0
+ ,69863
+ ,0
+ ,4261643
+ ,318765
+ ,1
+ ,6341
+ ,6341
+ ,63255
+ ,63255
+ ,1518120
+ ,286146
+ ,1
+ ,1164
+ ,1164
+ ,57320
+ ,57320
+ ,974440
+ ,306844
+ ,1
+ ,3310
+ ,3310
+ ,75230
+ ,75230
+ ,2482590
+ ,307705
+ ,0
+ ,1366
+ ,0
+ ,79420
+ ,0
+ ,555940
+ ,312448
+ ,0
+ ,965
+ ,0
+ ,73490
+ ,0
+ ,220470
+ ,299715
+ ,0
+ ,3256
+ ,0
+ ,35250
+ ,0
+ ,2326500
+ ,373399
+ ,1
+ ,1135
+ ,1135
+ ,62285
+ ,62285
+ ,1058845
+ ,299446
+ ,0
+ ,1270
+ ,0
+ ,69206
+ ,0
+ ,1799356
+ ,325586
+ ,0
+ ,661
+ ,0
+ ,65920
+ ,0
+ ,197760
+ ,291221
+ ,0
+ ,1013
+ ,0
+ ,69770
+ ,0
+ ,139540
+ ,261173
+ ,0
+ ,2844
+ ,0
+ ,72683
+ ,0
+ ,4869761
+ ,255027
+ ,1
+ ,11528
+ ,11528
+ ,-14545
+ ,-14545
+ ,-1018150
+ ,-78375
+ ,0
+ ,6526
+ ,0
+ ,55830
+ ,0
+ ,1451580
+ ,-58143
+ ,0
+ ,2264
+ ,0
+ ,55174
+ ,0
+ ,1324176
+ ,227033
+ ,1
+ ,4461
+ ,4461
+ ,67038
+ ,67038
+ ,6301572
+ ,235098
+ ,0
+ ,3999
+ ,0
+ ,51252
+ ,0
+ ,1537560
+ ,21267
+ ,0
+ ,35624
+ ,0
+ ,157278
+ ,0
+ ,35072994
+ ,238675
+ ,0
+ ,9252
+ ,0
+ ,79510
+ ,0
+ ,3816480
+ ,197687
+ ,0
+ ,15236
+ ,0
+ ,77440
+ ,0
+ ,6969600
+ ,418341
+ ,0
+ ,18073
+ ,0
+ ,27284
+ ,0
+ ,4911120
+ ,-297706)
+ ,dim=c(7
+ ,431)
+ ,dimnames=list(c('Group'
+ ,'Costs'
+ ,'GrCosts'
+ ,'Dividends'
+ ,'GrDiv'
+ ,'TrDiv'
+ ,'Wealth
')
+ ,1:431))
> y <- array(NA,dim=c(7,431),dimnames=list(c('Group','Costs','GrCosts','Dividends','GrDiv','TrDiv','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 = '7'
> #'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\r Group Costs GrCosts Dividends GrDiv TrDiv
1 6282929 1 162556 162556 213118 213118 230380558
2 4324047 1 29790 29790 81767 81767 25266003
3 4108272 1 87550 87550 153198 153198 70164684
4 -1212617 0 84738 0 -26007 0 -15292116
5 1485329 1 54660 54660 126942 126942 37955658
6 1779876 1 42634 42634 157214 157214 24525384
7 1367203 0 40949 0 129352 0 62218312
8 2519076 1 42312 42312 234817 234817 75845891
9 912684 1 37704 37704 60448 60448 27322496
10 1443586 1 16275 16275 47818 47818 5212162
11 1220017 0 25830 0 245546 0 28237790
12 984885 0 12679 0 48020 0 5282200
13 1457425 1 18014 18014 -1710 -1710 -408690
14 -572920 0 43556 0 32648 0 8064056
15 929144 1 24524 24524 95350 95350 47388950
16 1151176 0 6532 0 151352 0 15589256
17 790090 0 7123 0 288170 0 31410530
18 774497 1 20813 20813 114337 114337 57397174
19 990576 1 37597 37597 37884 37884 9395232
20 454195 0 17821 0 122844 0 45820812
21 876607 1 12988 12988 82340 82340 9798460
22 711969 1 22330 22330 79801 79801 6703284
23 702380 0 13326 0 165548 0 16885896
24 264449 0 16189 0 116384 0 34333280
25 450033 0 7146 0 134028 0 14072940
26 541063 0 15824 0 63838 0 4085632
27 588864 1 26088 26088 74996 74996 20023932
28 -37216 0 11326 0 31080 0 4009320
29 783310 0 8568 0 32168 0 1190216
30 467359 0 14416 0 49857 0 17998377
31 688779 1 3369 3369 87161 87161 2440508
32 608419 1 11819 11819 106113 106113 9019605
33 696348 1 6620 6620 80570 80570 3545080
34 597793 1 4519 4519 102129 102129 5004321
35 821730 0 2220 0 301670 0 6636740
36 377934 0 18562 0 102313 0 15858515
37 651939 0 10327 0 88577 0 8060507
38 697458 1 5336 5336 112477 112477 9110637
39 700368 1 2365 2365 191778 191778 15150462
40 225986 0 4069 0 79804 0 11571580
41 348695 0 7710 0 128294 0 104687904
42 373683 0 13718 0 96448 0 5883328
43 501709 0 4525 0 93811 0 21201286
44 413743 0 6869 0 117520 0 12339600
45 379825 0 4628 0 69159 0 4287858
46 336260 1 3653 3653 101792 101792 2443008
47 636765 1 1265 1265 210568 210568 5474768
48 481231 1 7489 7489 136996 136996 44112712
49 469107 0 4901 0 121920 0 10241280
50 211928 0 2284 0 76403 0 2521299
51 563925 1 3160 3160 108094 108094 11674152
52 511939 1 4150 4150 134759 134759 20213850
53 521016 1 7285 7285 188873 188873 21720395
54 543856 1 1134 1134 146216 146216 23686992
55 329304 1 4658 4658 156608 156608 24744064
56 423262 0 2384 0 61348 0 5950756
57 509665 0 3748 0 50350 0 453150
58 455881 0 5371 0 87720 0 5789520
59 367772 0 1285 0 99489 0 10645323
60 406339 1 9327 9327 87419 87419 8829319
61 493408 1 5565 5565 94355 94355 4434685
62 232942 0 1528 0 60326 0 2292388
63 416002 1 3122 3122 94670 94670 3218780
64 337430 1 7317 7317 82425 82425 6923700
65 361517 0 2675 0 59017 0 4662343
66 360962 0 13253 0 90829 0 86015063
67 235561 0 880 0 80791 0 5978534
68 408247 1 2053 2053 100423 100423 5322419
69 450296 0 1424 0 131116 0 12324904
70 418799 1 4036 4036 100269 100269 6316947
71 247405 1 3045 3045 27330 27330 1585140
72 378519 0 5119 0 39039 0 1912911
73 326638 0 1431 0 106885 0 3634090
74 328233 0 554 0 79285 0 872135
75 386225 0 1975 0 118881 0 4160835
76 283662 1 1286 1286 77623 77623 1319591
77 370225 0 1012 0 114768 0 5394096
78 269236 0 810 0 74015 0 3182645
79 365732 0 1280 0 69465 0 8127405
80 420383 1 666 666 117869 117869 20155599
81 345811 0 1380 0 60982 0 1585532
82 431809 1 4608 4608 90131 90131 6579563
83 418876 0 876 0 138971 0 8199289
84 297476 0 814 0 39625 0 713250
85 416776 0 514 0 102725 0 1540875
86 357257 1 5692 5692 64239 64239 4625208
87 458343 0 3642 0 90262 0 7762532
88 388386 0 540 0 103960 0 1455440
89 358934 0 2099 0 106611 0 6823104
90 407560 0 567 0 103345 0 1136795
91 392558 0 2001 0 95551 0 4968652
92 373177 1 2949 2949 82903 82903 3399023
93 428370 0 2253 0 63593 0 6295707
94 369419 1 6533 6533 126910 126910 9518250
95 358649 0 1889 0 37527 0 1688715
96 376641 1 3055 3055 60247 60247 2590621
97 467427 0 272 0 112995 0 903960
98 364885 1 1414 1414 70184 70184 13896432
99 436230 0 2564 0 130140 0 2863080
100 329118 1 1383 1383 73221 73221 805431
101 317365 1 1261 1261 76114 76114 2511762
102 286849 0 975 0 90534 0 2082282
103 376685 0 3366 0 108479 0 8678320
104 407198 0 576 0 113761 0 2047698
105 377772 0 1306 0 68696 0 1923488
106 271483 0 746 0 71561 0 1645903
107 153661 1 3192 3192 59831 59831 3589860
108 513294 1 2045 2045 97890 97890 1957800
109 324881 0 5477 0 101481 0 5987379
110 264512 1 1932 1932 72954 72954 2626344
111 420968 0 936 0 67939 0 2038170
112 129302 1 3437 3437 48022 48022 2257034
113 191521 0 5131 0 86111 0 6113881
114 268673 1 2397 2397 74020 74020 1036280
115 353179 1 1389 1389 57530 57530 517770
116 354624 0 1503 0 56364 0 2198196
117 363713 0 402 0 84990 0 2209740
118 456657 0 2239 0 88590 0 1860390
119 211742 1 2234 2234 77200 77200 1235200
120 338381 0 837 0 61262 0 4227078
121 418530 0 10579 0 110309 0 10148428
122 351483 0 875 0 67000 0 938000
123 372928 0 1395 0 93099 0 9589197
124 485538 1 1659 1659 107577 107577 3119733
125 279268 1 2647 2647 62920 62920 2328040
126 219060 1 3294 3294 75832 75832 1744136
127 325560 0 0 0 60720 0 0
128 325314 0 94 0 60793 0 425551
129 322046 0 422 0 57935 0 1622180
130 325560 0 0 0 60720 0 0
131 325599 0 34 0 60630 0 485040
132 377028 0 1558 0 55637 0 3505131
133 325560 0 0 0 60720 0 0
134 323850 0 43 0 60887 0 182661
135 325560 0 0 0 60720 0 0
136 331514 0 316 0 60505 0 544545
137 325632 0 115 0 60945 0 792285
138 325560 0 0 0 60720 0 0
139 325560 0 0 0 60720 0 0
140 325560 0 0 0 60720 0 0
141 322265 0 389 0 58990 0 825860
142 325560 0 0 0 60720 0 0
143 325906 0 1002 0 56750 0 851250
144 325985 0 36 0 60894 0 182682
145 346145 0 460 0 63346 0 950190
146 325898 0 309 0 56535 0 621885
147 325560 0 0 0 60720 0 0
148 325356 0 9 0 60835 0 365010
149 325560 0 0 0 60720 0 0
150 325930 0 14 0 61016 0 61016
151 318020 0 520 0 58650 0 586500
152 326389 0 1766 0 60438 0 4411974
153 325560 0 0 0 60720 0 0
154 302925 0 458 0 58625 0 644875
155 325540 0 20 0 60938 0 182814
156 325560 0 0 0 60720 0 0
157 325560 0 0 0 60720 0 0
158 326736 0 98 0 61490 0 122980
159 340580 0 405 0 60845 0 425915
160 325560 0 0 0 60720 0 0
161 325560 0 0 0 60720 0 0
162 325560 0 0 0 60720 0 0
163 325560 0 0 0 60720 0 0
164 331828 0 483 0 60830 0 1642410
165 323299 0 454 0 63261 0 3226311
166 325560 0 0 0 60720 0 0
167 325560 0 0 0 60720 0 0
168 387722 0 757 0 45689 0 868091
169 325560 0 0 0 60720 0 0
170 325560 0 0 0 60720 0 0
171 325560 0 0 0 60720 0 0
172 324598 0 36 0 61564 0 246256
173 325560 0 0 0 60720 0 0
174 328726 0 203 0 61938 0 557442
175 325560 0 0 0 60720 0 0
176 325043 0 90 0 60951 0 487608
177 325560 0 0 0 60720 0 0
178 325806 1 71 71 60745 60745 60745
179 325560 0 0 0 60720 0 0
180 325560 0 0 0 60720 0 0
181 387732 0 972 0 71642 0 2435828
182 349729 1 531 531 71641 71641 716410
183 332202 0 604 0 55792 0 2120096
184 305442 1 283 283 71873 71873 718730
185 329537 1 23 23 62555 62555 312775
186 327055 1 638 638 60370 60370 845180
187 356245 1 699 699 64873 64873 1037968
188 328451 0 149 0 62041 0 310205
189 307062 0 226 0 65745 0 328725
190 325560 0 0 0 60720 0 0
191 331345 0 275 0 59500 0 238000
192 325560 0 0 0 60720 0 0
193 331824 0 141 0 61630 0 369780
194 325560 0 0 0 60720 0 0
195 325685 0 28 0 60890 0 121780
196 325560 1 0 0 60720 60720 0
197 404480 1 2566 2566 113521 113521 9081680
198 325560 1 0 0 60720 60720 0
199 325560 1 0 0 60720 60720 0
200 318314 1 472 472 80045 80045 1600900
201 325560 1 0 0 60720 60720 0
202 325560 1 0 0 60720 60720 0
203 325560 1 0 0 60720 60720 0
204 311807 1 203 203 50804 50804 1371708
205 337724 1 496 496 87390 87390 1485630
206 326431 1 10 10 61656 61656 123312
207 327556 1 63 63 65688 65688 262752
208 325560 1 0 0 60720 60720 0
209 356850 1 1136 1136 48522 48522 1552704
210 325560 1 0 0 60720 60720 0
211 325560 0 0 0 60720 0 0
212 325560 0 0 0 60720 0 0
213 322741 0 267 0 57640 0 1844480
214 310902 0 474 0 61977 0 1239540
215 324295 0 534 0 62620 0 438340
216 325560 0 0 0 60720 0 0
217 326156 0 15 0 60831 0 486648
218 326960 0 397 0 60646 0 1698088
219 325560 0 0 0 60720 0 0
220 333411 0 1061 0 56225 0 1124500
221 297761 0 288 0 60510 0 242040
222 325560 0 0 0 60720 0 0
223 325536 0 3 0 60698 0 121396
224 325560 0 0 0 60720 0 0
225 325762 0 20 0 60805 0 121610
226 327957 0 278 0 61404 0 1596504
227 325560 0 0 0 60720 0 0
228 325560 0 0 0 60720 0 0
229 318521 0 192 0 65276 0 261104
230 325560 0 0 0 60720 0 0
231 319775 0 317 0 63915 0 575235
232 325560 0 0 0 60720 0 0
233 325560 0 0 0 60720 0 0
234 332128 0 368 0 61686 0 1048662
235 325560 0 0 0 60720 0 0
236 325486 0 2 0 60743 0 60743
237 325560 0 0 0 60720 0 0
238 325838 0 53 0 60349 0 362094
239 325560 0 0 0 60720 0 0
240 325560 0 0 0 60720 0 0
241 325560 0 0 0 60720 0 0
242 331767 0 94 0 61360 0 184080
243 325560 0 0 0 60720 0 0
244 324523 0 24 0 59818 0 478544
245 339995 0 2332 0 72680 0 290720
246 325560 1 0 0 60720 60720 0
247 325560 1 0 0 60720 60720 0
248 319582 0 131 0 61808 0 679888
249 325560 1 0 0 60720 60720 0
250 325560 1 0 0 60720 60720 0
251 307245 0 206 0 53110 0 477990
252 325560 1 0 0 60720 60720 0
253 317967 0 167 0 64245 0 128490
254 331488 0 622 0 73007 0 5329511
255 335452 0 885 0 82732 0 7032220
256 325560 0 0 0 60720 0 0
257 334184 0 365 0 54820 0 438560
258 313213 0 364 0 47705 0 1669675
259 325560 1 0 0 60720 60720 0
260 325560 0 0 0 60720 0 0
261 325560 1 0 0 60720 60720 0
262 325560 1 0 0 60720 60720 0
263 348678 0 226 0 72835 0 874020
264 328727 0 307 0 58856 0 882840
265 325560 1 0 0 60720 60720 0
266 325560 0 0 0 60720 0 0
267 325560 1 0 0 60720 60720 0
268 387978 0 188 0 77655 0 854205
269 325560 1 0 0 60720 60720 0
270 336704 0 138 0 69817 0 418902
271 325560 1 0 0 60720 60720 0
272 325560 1 0 0 60720 60720 0
273 325560 1 0 0 60720 60720 0
274 322076 0 125 0 60798 0 729576
275 325560 0 0 0 60720 0 0
276 334272 0 282 0 62452 0 1873560
277 338197 0 335 0 64175 0 2117775
278 325560 0 0 0 60720 0 0
279 321024 1 813 813 67440 67440 5530080
280 322145 0 176 0 68136 0 1907808
281 325560 1 0 0 60720 60720 0
282 325560 0 0 0 60720 0 0
283 323351 0 249 0 56726 0 4084272
284 325560 0 0 0 60720 0 0
285 327748 0 333 0 70811 0 920543
286 325560 0 0 0 60720 0 0
287 325560 1 0 0 60720 60720 0
288 328157 0 30 0 62045 0 248180
289 325560 1 0 0 60720 60720 0
290 311594 0 249 0 54323 0 3368026
291 325560 1 0 0 60720 60720 0
292 335962 0 165 0 62841 0 1508184
293 372426 0 453 0 81125 0 1703625
294 325560 0 0 0 60720 0 0
295 319844 0 53 0 59506 0 833084
296 355822 1 382 382 59365 59365 1246665
297 325560 0 0 0 60720 0 0
298 325560 1 0 0 60720 60720 0
299 325560 0 0 0 60720 0 0
300 325560 0 0 0 60720 0 0
301 324047 1 30 30 60798 60798 243192
302 311464 0 290 0 58790 0 117580
303 325560 0 0 0 60720 0 0
304 325560 1 0 0 60720 60720 0
305 353417 0 366 0 61808 0 3275824
306 325590 0 2 0 60735 0 546615
307 325560 1 0 0 60720 60720 0
308 328576 1 209 209 64016 64016 832208
309 326126 0 384 0 54683 0 1203026
310 325560 1 0 0 60720 60720 0
311 325560 1 0 0 60720 60720 0
312 369376 0 365 0 87192 0 7236936
313 325560 1 0 0 60720 60720 0
314 332013 1 49 49 64107 64107 512856
315 325871 0 3 0 60761 0 243044
316 342165 0 133 0 65990 0 923860
317 324967 0 32 0 59988 0 59988
318 314832 0 368 0 61167 0 1039839
319 325557 0 1 0 60719 0 364314
320 325560 0 0 0 60720 0 0
321 325560 0 0 0 60720 0 0
322 325560 0 0 0 60720 0 0
323 325560 0 0 0 60720 0 0
324 325560 1 0 0 60720 60720 0
325 325560 0 0 0 60720 0 0
326 322649 0 22 0 60722 0 121444
327 325560 1 0 0 60720 60720 0
328 325560 1 0 0 60720 60720 0
329 325560 0 0 0 60720 0 0
330 325560 0 0 0 60720 0 0
331 325560 0 0 0 60720 0 0
332 325560 0 0 0 60720 0 0
333 325560 0 0 0 60720 0 0
334 324598 0 96 0 60379 0 301895
335 325567 0 1 0 60727 0 121454
336 325560 0 0 0 60720 0 0
337 324005 0 81 0 60925 0 426475
338 325560 0 0 0 60720 0 0
339 325748 0 26 0 60896 0 60896
340 323385 0 125 0 59734 0 776542
341 315409 0 304 0 62969 0 944535
342 325560 0 0 0 60720 0 0
343 325560 0 0 0 60720 0 0
344 325560 0 0 0 60720 0 0
345 325560 0 0 0 60720 0 0
346 325560 0 0 0 60720 0 0
347 312275 0 119 0 59118 0 354708
348 325560 0 0 0 60720 0 0
349 325560 0 0 0 60720 0 0
350 325560 0 0 0 60720 0 0
351 320576 0 312 0 58598 0 820372
352 325246 0 60 0 61124 0 611240
353 332961 0 587 0 59595 0 715140
354 323010 0 135 0 62065 0 124130
355 325560 0 0 0 60720 0 0
356 325560 0 0 0 60720 0 0
357 345253 0 514 0 78780 0 4096560
358 325560 0 0 0 60720 0 0
359 325560 0 0 0 60720 0 0
360 325560 0 0 0 60720 0 0
361 325559 0 1 0 60722 0 242888
362 325560 0 0 0 60720 0 0
363 325560 0 0 0 60720 0 0
364 319634 1 58 58 61600 61600 184800
365 319951 0 180 0 59635 0 655985
366 325560 0 0 0 60720 0 0
367 325560 0 0 0 60720 0 0
368 325560 0 0 0 60720 0 0
369 325560 0 0 0 60720 0 0
370 325560 0 0 0 60720 0 0
371 325560 0 0 0 60720 0 0
372 318519 0 448 0 59781 0 2391240
373 343222 0 227 0 76644 0 689796
374 317234 0 174 0 64820 0 64820
375 325560 0 0 0 60720 0 0
376 325560 0 0 0 60720 0 0
377 314025 0 121 0 56178 0 1348272
378 320249 0 607 0 60436 0 664796
379 325560 0 0 0 60720 0 0
380 325560 0 0 0 60720 0 0
381 325560 0 0 0 60720 0 0
382 349365 0 530 0 73433 0 4405980
383 289197 0 571 0 41477 0 3318160
384 325560 0 0 0 60720 0 0
385 329245 0 78 0 62700 0 1003200
386 240869 0 2489 0 67804 0 2712160
387 327182 0 131 0 59661 0 357966
388 322876 0 923 0 58620 0 468960
389 323117 0 72 0 60398 0 181194
390 306351 0 572 0 58580 0 937280
391 335137 0 397 0 62710 0 627100
392 308271 0 450 0 59325 0 474600
393 301731 0 622 0 60950 0 426650
394 382409 0 694 0 68060 0 544480
395 279230 1 3425 3425 83620 83620 1003440
396 298731 0 562 0 58456 0 759928
397 243650 0 4917 0 52811 0 2218062
398 532682 1 1442 1442 121173 121173 14298414
399 319771 0 529 0 63870 0 574830
400 171493 1 2126 2126 21001 21001 2898138
401 347262 0 1061 0 70415 0 352075
402 343945 0 776 0 64230 0 578070
403 311874 0 611 0 59190 0 473520
404 302211 1 1526 1526 69351 69351 1733775
405 316708 0 592 0 64270 0 449890
406 333463 0 1182 0 70694 0 919022
407 344282 0 621 0 68005 0 1088080
408 319635 0 989 0 58930 0 648230
409 301186 0 438 0 58320 0 641520
410 300381 0 726 0 69980 0 209940
411 318765 0 1303 0 69863 0 4261643
412 286146 1 6341 6341 63255 63255 1518120
413 306844 1 1164 1164 57320 57320 974440
414 307705 1 3310 3310 75230 75230 2482590
415 312448 0 1366 0 79420 0 555940
416 299715 0 965 0 73490 0 220470
417 373399 0 3256 0 35250 0 2326500
418 299446 1 1135 1135 62285 62285 1058845
419 325586 0 1270 0 69206 0 1799356
420 291221 0 661 0 65920 0 197760
421 261173 0 1013 0 69770 0 139540
422 255027 0 2844 0 72683 0 4869761
423 -78375 1 11528 11528 -14545 -14545 -1018150
424 -58143 0 6526 0 55830 0 1451580
425 227033 0 2264 0 55174 0 1324176
426 235098 1 4461 4461 67038 67038 6301572
427 21267 0 3999 0 51252 0 1537560
428 238675 0 35624 0 157278 0 35072994
429 197687 0 9252 0 79510 0 3816480
430 418341 0 15236 0 77440 0 6969600
431 -297706 0 18073 0 27284 0 4911120
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs GrCosts Dividends GrDiv
1.373e+05 5.217e+04 -8.643e+00 4.186e+01 3.063e+00 -1.623e+00
TrDiv
2.370e-03
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-792239 -17114 2279 20417 2967578
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.373e+05 3.260e+04 4.212 3.09e-05 ***
Group 5.217e+04 5.552e+04 0.940 0.3480
Costs -8.643e+00 1.790e+00 -4.830 1.91e-06 ***
GrCosts 4.186e+01 2.185e+00 19.155 < 2e-16 ***
Dividends 3.063e+00 4.550e-01 6.731 5.50e-11 ***
GrDiv -1.623e+00 6.978e-01 -2.326 0.0205 *
TrDiv 2.370e-03 1.218e-03 1.947 0.0522 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 209400 on 424 degrees of freedom
Multiple R-squared: 0.7841, Adjusted R-squared: 0.781
F-statistic: 256.6 on 6 and 424 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.452777e-88 1.726388e-88
[2,] 1.0000000 1.116770e-90 5.583848e-91
[3,] 1.0000000 1.013903e-97 5.069515e-98
[4,] 1.0000000 2.540066e-112 1.270033e-112
[5,] 1.0000000 3.368705e-121 1.684353e-121
[6,] 1.0000000 3.125641e-130 1.562821e-130
[7,] 1.0000000 1.723901e-140 8.619505e-141
[8,] 1.0000000 1.208109e-143 6.040546e-144
[9,] 1.0000000 8.019670e-149 4.009835e-149
[10,] 1.0000000 1.615477e-153 8.077383e-154
[11,] 1.0000000 5.550693e-154 2.775347e-154
[12,] 1.0000000 3.388307e-158 1.694153e-158
[13,] 1.0000000 2.823414e-160 1.411707e-160
[14,] 1.0000000 3.532197e-161 1.766099e-161
[15,] 1.0000000 9.437204e-162 4.718602e-162
[16,] 1.0000000 3.162630e-161 1.581315e-161
[17,] 1.0000000 4.853835e-163 2.426918e-163
[18,] 1.0000000 2.237874e-166 1.118937e-166
[19,] 1.0000000 1.978202e-168 9.891010e-169
[20,] 1.0000000 3.636812e-181 1.818406e-181
[21,] 1.0000000 3.021362e-183 1.510681e-183
[22,] 1.0000000 8.191577e-187 4.095788e-187
[23,] 1.0000000 2.583059e-188 1.291530e-188
[24,] 1.0000000 1.366914e-194 6.834571e-195
[25,] 1.0000000 5.616981e-196 2.808491e-196
[26,] 1.0000000 1.843496e-195 9.217478e-196
[27,] 1.0000000 3.605825e-195 1.802912e-195
[28,] 1.0000000 1.799987e-202 8.999936e-203
[29,] 1.0000000 1.953831e-207 9.769155e-208
[30,] 1.0000000 1.050539e-206 5.252694e-207
[31,] 1.0000000 3.132165e-207 1.566082e-207
[32,] 1.0000000 4.350185e-210 2.175092e-210
[33,] 1.0000000 3.917482e-210 1.958741e-210
[34,] 1.0000000 2.622678e-210 1.311339e-210
[35,] 1.0000000 1.706963e-209 8.534815e-210
[36,] 1.0000000 5.423556e-209 2.711778e-209
[37,] 1.0000000 2.498640e-208 1.249320e-208
[38,] 1.0000000 2.862180e-207 1.431090e-207
[39,] 1.0000000 9.773966e-207 4.886983e-207
[40,] 1.0000000 3.804399e-206 1.902199e-206
[41,] 1.0000000 1.378847e-206 6.894234e-207
[42,] 1.0000000 2.321784e-207 1.160892e-207
[43,] 1.0000000 1.687731e-206 8.438656e-207
[44,] 1.0000000 8.142705e-206 4.071352e-206
[45,] 1.0000000 9.629126e-205 4.814563e-205
[46,] 1.0000000 3.802545e-208 1.901272e-208
[47,] 1.0000000 3.461233e-208 1.730616e-208
[48,] 1.0000000 9.473927e-213 4.736963e-213
[49,] 1.0000000 2.054872e-213 1.027436e-213
[50,] 1.0000000 1.193569e-212 5.967844e-213
[51,] 1.0000000 3.916806e-213 1.958403e-213
[52,] 1.0000000 2.394822e-214 1.197411e-214
[53,] 1.0000000 5.187974e-214 2.593987e-214
[54,] 1.0000000 3.541166e-213 1.770583e-213
[55,] 1.0000000 8.071033e-213 4.035516e-213
[56,] 1.0000000 3.526063e-212 1.763031e-212
[57,] 1.0000000 8.156955e-212 4.078478e-212
[58,] 1.0000000 6.224956e-213 3.112478e-213
[59,] 1.0000000 7.098904e-212 3.549452e-212
[60,] 1.0000000 4.351694e-211 2.175847e-211
[61,] 1.0000000 2.333292e-210 1.166646e-210
[62,] 1.0000000 1.274030e-209 6.370152e-210
[63,] 1.0000000 2.408010e-211 1.204005e-211
[64,] 1.0000000 3.465855e-211 1.732927e-211
[65,] 1.0000000 2.859098e-210 1.429549e-210
[66,] 1.0000000 1.541681e-209 7.708406e-210
[67,] 1.0000000 7.022271e-209 3.511135e-209
[68,] 1.0000000 1.473284e-208 7.366421e-209
[69,] 1.0000000 2.260571e-208 1.130286e-208
[70,] 1.0000000 2.479037e-207 1.239519e-207
[71,] 1.0000000 4.314650e-208 2.157325e-208
[72,] 1.0000000 3.483099e-207 1.741549e-207
[73,] 1.0000000 3.733527e-207 1.866763e-207
[74,] 1.0000000 1.516982e-207 7.584909e-208
[75,] 1.0000000 1.314937e-206 6.574683e-207
[76,] 1.0000000 1.372789e-205 6.863944e-206
[77,] 1.0000000 2.620109e-206 1.310055e-206
[78,] 1.0000000 4.354709e-206 2.177355e-206
[79,] 1.0000000 3.434471e-205 1.717235e-205
[80,] 1.0000000 7.197697e-205 3.598849e-205
[81,] 1.0000000 7.193960e-204 3.596980e-204
[82,] 1.0000000 7.392572e-203 3.696286e-203
[83,] 1.0000000 4.185676e-202 2.092838e-202
[84,] 1.0000000 2.787888e-202 1.393944e-202
[85,] 1.0000000 1.773373e-201 8.866867e-202
[86,] 1.0000000 4.271945e-202 2.135972e-202
[87,] 1.0000000 1.468773e-202 7.343867e-203
[88,] 1.0000000 1.444069e-201 7.220344e-202
[89,] 1.0000000 1.300549e-200 6.502746e-201
[90,] 1.0000000 9.332337e-200 4.666168e-200
[91,] 1.0000000 9.652331e-199 4.826166e-199
[92,] 1.0000000 8.337135e-198 4.168568e-198
[93,] 1.0000000 2.551452e-198 1.275726e-198
[94,] 1.0000000 6.914171e-198 3.457086e-198
[95,] 1.0000000 1.499317e-197 7.496587e-198
[96,] 1.0000000 8.707725e-197 4.353863e-197
[97,] 1.0000000 1.305898e-196 6.529488e-197
[98,] 1.0000000 1.534625e-197 7.673125e-198
[99,] 1.0000000 2.001522e-199 1.000761e-199
[100,] 1.0000000 6.767041e-199 3.383520e-199
[101,] 1.0000000 2.626818e-198 1.313409e-198
[102,] 1.0000000 3.198219e-198 1.599109e-198
[103,] 1.0000000 5.192393e-199 2.596196e-199
[104,] 1.0000000 1.432547e-201 7.162734e-202
[105,] 1.0000000 9.318685e-201 4.659342e-201
[106,] 1.0000000 2.544089e-200 1.272044e-200
[107,] 1.0000000 7.448575e-200 3.724287e-200
[108,] 1.0000000 4.013331e-199 2.006666e-199
[109,] 1.0000000 6.254758e-199 3.127379e-199
[110,] 1.0000000 1.086817e-199 5.434086e-200
[111,] 1.0000000 1.157283e-198 5.786415e-199
[112,] 1.0000000 1.734923e-198 8.674614e-199
[113,] 1.0000000 1.577607e-197 7.888033e-198
[114,] 1.0000000 1.221069e-197 6.105345e-198
[115,] 1.0000000 5.784328e-198 2.892164e-198
[116,] 1.0000000 5.786748e-197 2.893374e-197
[117,] 1.0000000 7.318241e-197 3.659120e-197
[118,] 1.0000000 7.760798e-196 3.880399e-196
[119,] 1.0000000 8.217113e-195 4.108557e-195
[120,] 1.0000000 8.575385e-194 4.287693e-194
[121,] 1.0000000 8.937810e-193 4.468905e-193
[122,] 1.0000000 9.313576e-192 4.656788e-192
[123,] 1.0000000 1.031753e-191 5.158767e-192
[124,] 1.0000000 1.064540e-190 5.322701e-191
[125,] 1.0000000 1.098336e-189 5.491678e-190
[126,] 1.0000000 1.120567e-188 5.602836e-189
[127,] 1.0000000 1.072084e-187 5.360420e-188
[128,] 1.0000000 1.089305e-186 5.446527e-187
[129,] 1.0000000 1.092973e-185 5.464865e-186
[130,] 1.0000000 1.090585e-184 5.452927e-185
[131,] 1.0000000 1.082203e-183 5.411016e-184
[132,] 1.0000000 1.047805e-182 5.239025e-183
[133,] 1.0000000 1.028388e-181 5.141938e-182
[134,] 1.0000000 7.139632e-181 3.569816e-181
[135,] 1.0000000 6.951467e-180 3.475734e-180
[136,] 1.0000000 5.921156e-179 2.960578e-179
[137,] 1.0000000 4.950049e-178 2.475024e-178
[138,] 1.0000000 4.733987e-177 2.366994e-177
[139,] 1.0000000 4.534827e-176 2.267414e-176
[140,] 1.0000000 4.290700e-175 2.145350e-175
[141,] 1.0000000 4.047376e-174 2.023688e-174
[142,] 1.0000000 3.687696e-173 1.843848e-173
[143,] 1.0000000 3.319920e-172 1.659960e-172
[144,] 1.0000000 3.075520e-171 1.537760e-171
[145,] 1.0000000 2.810998e-170 1.405499e-170
[146,] 1.0000000 2.589188e-169 1.294594e-169
[147,] 1.0000000 2.361650e-168 1.180825e-168
[148,] 1.0000000 2.142856e-167 1.071428e-167
[149,] 1.0000000 1.935766e-166 9.678829e-167
[150,] 1.0000000 1.422724e-165 7.113620e-166
[151,] 1.0000000 1.271390e-164 6.356950e-165
[152,] 1.0000000 1.130281e-163 5.651406e-164
[153,] 1.0000000 9.996535e-163 4.998268e-163
[154,] 1.0000000 8.795756e-162 4.397878e-162
[155,] 1.0000000 7.540173e-161 3.770087e-161
[156,] 1.0000000 5.945882e-160 2.972941e-160
[157,] 1.0000000 5.155640e-159 2.577820e-159
[158,] 1.0000000 4.447652e-158 2.223826e-158
[159,] 1.0000000 1.550444e-159 7.752222e-160
[160,] 1.0000000 1.350113e-158 6.750564e-159
[161,] 1.0000000 1.169713e-157 5.848565e-158
[162,] 1.0000000 1.008297e-156 5.041485e-157
[163,] 1.0000000 8.748878e-156 4.374439e-156
[164,] 1.0000000 7.466298e-155 3.733149e-155
[165,] 1.0000000 6.355348e-154 3.177674e-154
[166,] 1.0000000 5.369411e-153 2.684705e-153
[167,] 1.0000000 4.552316e-152 2.276158e-152
[168,] 1.0000000 3.807847e-151 1.903924e-151
[169,] 1.0000000 3.121868e-150 1.560934e-150
[170,] 1.0000000 2.585894e-149 1.292947e-149
[171,] 1.0000000 2.131286e-148 1.065643e-148
[172,] 1.0000000 1.193891e-147 5.969454e-148
[173,] 1.0000000 9.282439e-147 4.641220e-147
[174,] 1.0000000 5.850345e-146 2.925172e-146
[175,] 1.0000000 2.937642e-145 1.468821e-145
[176,] 1.0000000 2.349733e-144 1.174867e-144
[177,] 1.0000000 1.819147e-143 9.095733e-144
[178,] 1.0000000 1.028350e-142 5.141752e-143
[179,] 1.0000000 8.223127e-142 4.111563e-142
[180,] 1.0000000 5.176239e-141 2.588119e-141
[181,] 1.0000000 4.094603e-140 2.047301e-140
[182,] 1.0000000 2.861336e-139 1.430668e-139
[183,] 1.0000000 2.241699e-138 1.120850e-138
[184,] 1.0000000 1.714773e-137 8.573863e-138
[185,] 1.0000000 1.330330e-136 6.651651e-137
[186,] 1.0000000 1.029794e-135 5.148971e-136
[187,] 1.0000000 7.866254e-135 3.933127e-135
[188,] 1.0000000 4.823948e-134 2.411974e-134
[189,] 1.0000000 3.656069e-133 1.828034e-133
[190,] 1.0000000 2.760555e-132 1.380277e-132
[191,] 1.0000000 1.084740e-131 5.423698e-132
[192,] 1.0000000 8.152679e-131 4.076340e-131
[193,] 1.0000000 6.103551e-130 3.051775e-130
[194,] 1.0000000 4.551389e-129 2.275695e-129
[195,] 1.0000000 3.377819e-128 1.688909e-128
[196,] 1.0000000 1.361859e-127 6.809295e-128
[197,] 1.0000000 1.006516e-126 5.032578e-127
[198,] 1.0000000 7.252788e-126 3.626394e-126
[199,] 1.0000000 5.314400e-125 2.657200e-125
[200,] 1.0000000 7.069468e-125 3.534734e-125
[201,] 1.0000000 5.170144e-124 2.585072e-124
[202,] 1.0000000 3.761664e-123 1.880832e-123
[203,] 1.0000000 2.723710e-122 1.361855e-122
[204,] 1.0000000 1.929688e-121 9.648439e-122
[205,] 1.0000000 1.356203e-120 6.781016e-121
[206,] 1.0000000 9.698371e-120 4.849186e-120
[207,] 1.0000000 6.892018e-119 3.446009e-119
[208,] 1.0000000 4.908200e-118 2.454100e-118
[209,] 1.0000000 3.453363e-117 1.726681e-117
[210,] 1.0000000 2.418952e-116 1.209476e-116
[211,] 1.0000000 9.568281e-116 4.784141e-116
[212,] 1.0000000 6.240020e-115 3.120010e-115
[213,] 1.0000000 4.320525e-114 2.160263e-114
[214,] 1.0000000 2.983022e-113 1.491511e-113
[215,] 1.0000000 2.045539e-112 1.022770e-112
[216,] 1.0000000 1.397904e-111 6.989519e-112
[217,] 1.0000000 9.589050e-111 4.794525e-111
[218,] 1.0000000 6.481020e-110 3.240510e-110
[219,] 1.0000000 4.359178e-109 2.179589e-109
[220,] 1.0000000 2.796340e-108 1.398170e-108
[221,] 1.0000000 1.864219e-107 9.321095e-108
[222,] 1.0000000 1.235201e-106 6.176003e-107
[223,] 1.0000000 8.158345e-106 4.079172e-106
[224,] 1.0000000 5.362319e-105 2.681160e-105
[225,] 1.0000000 3.426831e-104 1.713416e-104
[226,] 1.0000000 2.230604e-103 1.115302e-103
[227,] 1.0000000 1.447305e-102 7.236526e-103
[228,] 1.0000000 9.329293e-102 4.664647e-102
[229,] 1.0000000 5.966034e-101 2.983017e-101
[230,] 1.0000000 3.808229e-100 1.904115e-100
[231,] 1.0000000 2.418965e-99 1.209483e-99
[232,] 1.0000000 1.528982e-98 7.644909e-99
[233,] 1.0000000 9.382517e-98 4.691259e-98
[234,] 1.0000000 5.872363e-97 2.936182e-97
[235,] 1.0000000 3.654672e-96 1.827336e-96
[236,] 1.0000000 2.055958e-95 1.027979e-95
[237,] 1.0000000 1.276480e-94 6.382402e-95
[238,] 1.0000000 7.889419e-94 3.944710e-94
[239,] 1.0000000 4.887979e-93 2.443990e-93
[240,] 1.0000000 2.992800e-92 1.496400e-92
[241,] 1.0000000 1.823974e-91 9.119868e-92
[242,] 1.0000000 1.026842e-90 5.134212e-91
[243,] 1.0000000 6.202076e-90 3.101038e-90
[244,] 1.0000000 3.672052e-89 1.836026e-89
[245,] 1.0000000 1.263442e-88 6.317208e-89
[246,] 1.0000000 1.106290e-88 5.531451e-89
[247,] 1.0000000 6.616068e-88 3.308034e-88
[248,] 1.0000000 2.560101e-87 1.280050e-87
[249,] 1.0000000 1.027252e-86 5.136260e-87
[250,] 1.0000000 6.118778e-86 3.059389e-86
[251,] 1.0000000 3.598562e-85 1.799281e-85
[252,] 1.0000000 2.122783e-84 1.061391e-84
[253,] 1.0000000 1.246087e-83 6.230437e-84
[254,] 1.0000000 6.949625e-83 3.474813e-83
[255,] 1.0000000 3.739222e-82 1.869611e-82
[256,] 1.0000000 2.164561e-81 1.082281e-81
[257,] 1.0000000 1.236951e-80 6.184755e-81
[258,] 1.0000000 7.089634e-80 3.544817e-80
[259,] 1.0000000 4.052556e-79 2.026278e-79
[260,] 1.0000000 2.299356e-78 1.149678e-78
[261,] 1.0000000 1.246353e-77 6.231763e-78
[262,] 1.0000000 7.002271e-77 3.501136e-77
[263,] 1.0000000 3.913419e-76 1.956710e-76
[264,] 1.0000000 2.175509e-75 1.087755e-75
[265,] 1.0000000 1.212915e-74 6.064574e-75
[266,] 1.0000000 6.630226e-74 3.315113e-74
[267,] 1.0000000 3.648647e-73 1.824324e-73
[268,] 1.0000000 2.005628e-72 1.002814e-72
[269,] 1.0000000 1.079376e-71 5.396879e-72
[270,] 1.0000000 4.429022e-71 2.214511e-71
[271,] 1.0000000 1.928417e-70 9.642085e-71
[272,] 1.0000000 1.034509e-69 5.172546e-70
[273,] 1.0000000 5.474985e-69 2.737493e-69
[274,] 1.0000000 2.928042e-68 1.464021e-68
[275,] 1.0000000 1.533349e-67 7.666743e-68
[276,] 1.0000000 6.978559e-67 3.489279e-67
[277,] 1.0000000 3.623009e-66 1.811505e-66
[278,] 1.0000000 1.885446e-65 9.427232e-66
[279,] 1.0000000 9.779558e-65 4.889779e-65
[280,] 1.0000000 5.034169e-64 2.517084e-64
[281,] 1.0000000 2.596877e-63 1.298438e-63
[282,] 1.0000000 1.322051e-62 6.610256e-63
[283,] 1.0000000 6.724089e-62 3.362045e-62
[284,] 1.0000000 3.033616e-61 1.516808e-61
[285,] 1.0000000 1.511294e-60 7.556470e-61
[286,] 1.0000000 7.580025e-60 3.790013e-60
[287,] 1.0000000 2.842002e-59 1.421001e-59
[288,] 1.0000000 1.395798e-58 6.978992e-59
[289,] 1.0000000 6.826474e-58 3.413237e-58
[290,] 1.0000000 3.315645e-57 1.657823e-57
[291,] 1.0000000 1.601228e-56 8.006141e-57
[292,] 1.0000000 7.748859e-56 3.874430e-56
[293,] 1.0000000 3.718024e-55 1.859012e-55
[294,] 1.0000000 1.765136e-54 8.825679e-55
[295,] 1.0000000 8.333059e-54 4.166529e-54
[296,] 1.0000000 3.726724e-53 1.863362e-53
[297,] 1.0000000 1.755425e-52 8.777126e-53
[298,] 1.0000000 8.140036e-52 4.070018e-52
[299,] 1.0000000 3.802874e-51 1.901437e-51
[300,] 1.0000000 1.462117e-50 7.310584e-51
[301,] 1.0000000 6.658209e-50 3.329105e-50
[302,] 1.0000000 3.004474e-49 1.502237e-49
[303,] 1.0000000 3.016116e-49 1.508058e-49
[304,] 1.0000000 1.350223e-48 6.751113e-49
[305,] 1.0000000 6.083106e-48 3.041553e-48
[306,] 1.0000000 2.754141e-47 1.377070e-47
[307,] 1.0000000 1.250200e-46 6.251001e-47
[308,] 1.0000000 5.518911e-46 2.759455e-46
[309,] 1.0000000 2.467245e-45 1.233622e-45
[310,] 1.0000000 1.092645e-44 5.463227e-45
[311,] 1.0000000 4.775063e-44 2.387531e-44
[312,] 1.0000000 2.073515e-43 1.036757e-43
[313,] 1.0000000 8.946341e-43 4.473171e-43
[314,] 1.0000000 3.835090e-42 1.917545e-42
[315,] 1.0000000 1.584635e-41 7.923176e-42
[316,] 1.0000000 6.708925e-41 3.354462e-41
[317,] 1.0000000 2.844830e-40 1.422415e-40
[318,] 1.0000000 1.134893e-39 5.674466e-40
[319,] 1.0000000 4.363800e-39 2.181900e-39
[320,] 1.0000000 1.802409e-38 9.012043e-39
[321,] 1.0000000 7.394330e-38 3.697165e-38
[322,] 1.0000000 3.012865e-37 1.506432e-37
[323,] 1.0000000 1.219195e-36 6.095976e-37
[324,] 1.0000000 4.899553e-36 2.449777e-36
[325,] 1.0000000 1.955758e-35 9.778791e-36
[326,] 1.0000000 7.771476e-35 3.885738e-35
[327,] 1.0000000 3.057872e-34 1.528936e-34
[328,] 1.0000000 1.206241e-33 6.031205e-34
[329,] 1.0000000 4.678753e-33 2.339376e-33
[330,] 1.0000000 1.803952e-32 9.019758e-33
[331,] 1.0000000 6.933205e-32 3.466603e-32
[332,] 1.0000000 2.643434e-31 1.321717e-31
[333,] 1.0000000 9.963176e-31 4.981588e-31
[334,] 1.0000000 3.726956e-30 1.863478e-30
[335,] 1.0000000 1.383590e-29 6.917951e-30
[336,] 1.0000000 5.097136e-29 2.548568e-29
[337,] 1.0000000 1.863281e-28 9.316407e-29
[338,] 1.0000000 6.856102e-28 3.428051e-28
[339,] 1.0000000 2.467495e-27 1.233747e-27
[340,] 1.0000000 8.809834e-27 4.404917e-27
[341,] 1.0000000 3.120153e-26 1.560076e-26
[342,] 1.0000000 1.085361e-25 5.426805e-26
[343,] 1.0000000 3.829514e-25 1.914757e-25
[344,] 1.0000000 1.211460e-24 6.057298e-25
[345,] 1.0000000 4.210031e-24 2.105016e-24
[346,] 1.0000000 1.431276e-23 7.156381e-24
[347,] 1.0000000 4.824140e-23 2.412070e-23
[348,] 1.0000000 1.109348e-22 5.546740e-23
[349,] 1.0000000 3.709395e-22 1.854697e-22
[350,] 1.0000000 1.229326e-21 6.146631e-22
[351,] 1.0000000 4.037492e-21 2.018746e-21
[352,] 1.0000000 1.320852e-20 6.604261e-21
[353,] 1.0000000 4.259123e-20 2.129562e-20
[354,] 1.0000000 1.360543e-19 6.802715e-20
[355,] 1.0000000 4.179547e-19 2.089773e-19
[356,] 1.0000000 1.321716e-18 6.608581e-19
[357,] 1.0000000 4.107394e-18 2.053697e-18
[358,] 1.0000000 1.263860e-17 6.319299e-18
[359,] 1.0000000 3.850108e-17 1.925054e-17
[360,] 1.0000000 1.160975e-16 5.804877e-17
[361,] 1.0000000 3.464826e-16 1.732413e-16
[362,] 1.0000000 1.023240e-15 5.116200e-16
[363,] 1.0000000 3.058567e-15 1.529284e-15
[364,] 1.0000000 8.406479e-15 4.203240e-15
[365,] 1.0000000 2.456910e-14 1.228455e-14
[366,] 1.0000000 6.985504e-14 3.492752e-14
[367,] 1.0000000 1.963656e-13 9.818279e-14
[368,] 1.0000000 5.484102e-13 2.742051e-13
[369,] 1.0000000 1.494387e-12 7.471937e-13
[370,] 1.0000000 4.055302e-12 2.027651e-12
[371,] 1.0000000 1.087088e-11 5.435441e-12
[372,] 1.0000000 2.877936e-11 1.438968e-11
[373,] 1.0000000 7.099912e-11 3.549956e-11
[374,] 1.0000000 1.640420e-10 8.202102e-11
[375,] 1.0000000 4.205381e-10 2.102690e-10
[376,] 1.0000000 1.084734e-09 5.423671e-10
[377,] 1.0000000 2.168767e-09 1.084384e-09
[378,] 1.0000000 5.321133e-09 2.660566e-09
[379,] 1.0000000 1.211856e-08 6.059282e-09
[380,] 1.0000000 2.940525e-08 1.470263e-08
[381,] 1.0000000 7.087128e-08 3.543564e-08
[382,] 0.9999999 1.645997e-07 8.229986e-08
[383,] 0.9999998 3.847128e-07 1.923564e-07
[384,] 0.9999996 8.965229e-07 4.482614e-07
[385,] 0.9999991 1.745805e-06 8.729023e-07
[386,] 0.9999982 3.699572e-06 1.849786e-06
[387,] 0.9999959 8.199670e-06 4.099835e-06
[388,] 0.9999924 1.518332e-05 7.591661e-06
[389,] 0.9999835 3.295854e-05 1.647927e-05
[390,] 0.9999651 6.981304e-05 3.490652e-05
[391,] 0.9999277 1.446871e-04 7.234353e-05
[392,] 0.9998549 2.902640e-04 1.451320e-04
[393,] 0.9997280 5.439182e-04 2.719591e-04
[394,] 0.9994846 1.030736e-03 5.153678e-04
[395,] 0.9990009 1.998161e-03 9.990805e-04
[396,] 0.9981237 3.752500e-03 1.876250e-03
[397,] 0.9965443 6.911341e-03 3.455670e-03
[398,] 0.9938936 1.221289e-02 6.106446e-03
[399,] 0.9901006 1.979889e-02 9.899443e-03
[400,] 0.9834424 3.311517e-02 1.655758e-02
[401,] 0.9719489 5.610212e-02 2.805106e-02
[402,] 0.9539352 9.212952e-02 4.606476e-02
[403,] 0.9334672 1.330655e-01 6.653276e-02
[404,] 0.8964155 2.071691e-01 1.035845e-01
[405,] 0.8482379 3.035243e-01 1.517621e-01
[406,] 0.7799659 4.400682e-01 2.200341e-01
[407,] 0.6926359 6.147283e-01 3.073641e-01
[408,] 0.9079764 1.840472e-01 9.202358e-02
[409,] 0.8403355 3.193291e-01 1.596645e-01
[410,] 0.7553978 4.892044e-01 2.446022e-01
[411,] 0.6242974 7.514052e-01 3.757026e-01
[412,] 0.4614339 9.228678e-01 5.385661e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1avpv1291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/2kmog1291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/3kmog1291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/4kmog1291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/5vdnj1291330748.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
-1.582183e+05 2.967578e+06 6.242201e+05 -5.016268e+05 -7.922385e+05
6 7 8 9 10
-1.100347e+05 1.040185e+06 4.065020e+05 -6.808063e+05 6.323843e+05
11 12 13 14 15
4.869825e+05 7.975685e+05 6.730946e+05 -4.528743e+05 -3.244176e+05
16 17 18 19 20
5.698190e+05 -2.426991e+05 -4.068701e+05 -5.243865e+05 -1.392642e+04
21 22 23 24 25
1.140024e+05 -3.499120e+05 1.331922e+05 -1.707664e+05 -6.936383e+04
26 27 28 29 30
3.353181e+05 -6.224811e+05 -1.813264e+05 6.187109e+05 2.592913e+05
31 32 33 34 35
2.561397e+05 -1.477407e+05 1.626073e+05 9.933395e+04 -2.360687e+05
36 37 38 39 40
5.011134e+04 3.134921e+05 1.472344e+05 1.203332e+05 -1.480014e+05
41 42 43 44 45
-3.630277e+05 4.559789e+04 6.593958e+04 -5.337955e+04 6.053517e+04
46 47 48 49 50
-1.268817e+05 8.914354e+04 -2.587551e+05 -2.352834e+04 -1.456213e+05
51 52 53 54 55
8.620553e+04 -5.728963e+04 -2.338111e+05 5.006892e+04 -2.989893e+05
56 57 58 59 60
1.045591e+05 2.494654e+05 8.260581e+04 -8.837342e+04 -2.396914e+05
61 62 63 64 65
-2.724862e+04 -8.135870e+04 -2.108917e+04 -2.301377e+05 5.552223e+04
66 67 68 69 70
-1.438492e+05 -1.557565e+05 -6.608753e+03 -1.054951e+05 -6.405175e+04
71 72 73 74 75
-8.630760e+04 1.613525e+05 -1.342808e+05 -4.918722e+04 -1.079814e+05
76 77 78 79 80
-6.340797e+04 -1.226305e+05 -9.530678e+04 7.467767e+03 -8.680032e+03
81 82 83 84 85
2.989722e+04 -5.606581e+04 -1.559318e+05 4.414927e+04 -3.436637e+04
86 87 88 89 90
-1.247111e+05 5.766209e+04 -6.611169e+04 -1.029301e+05 -4.406551e+04
91 92 93 94 95
-3.188368e+04 -4.165470e+04 1.008414e+05 -2.423047e+05 1.187274e+05
96 97 98 99 100
-7.177381e+03 -1.575220e+04 -5.534011e+03 -8.429364e+04 -1.361739e+04
101 102 103 104 105
-2.952767e+04 -1.242537e+05 -8.434649e+04 -7.841049e+04 3.679140e+04
106 107 108 109 110
-8.245473e+04 -2.364768e+05 1.103250e+05 -9.009382e+04 -1.003881e+05
111 112 113 114 115
7.883615e+04 -2.488125e+05 -1.796693e+05 -1.094369e+05 3.351644e+04
116 117 118 119 120
5.246518e+04 -3.566432e+04 6.295913e+04 -1.660040e+05 1.065580e+04
121 122 123 124 125
1.075239e+04 1.430733e+04 -6.019245e+04 7.868870e+04 -9.422583e+04
126 127 128 129 130
-1.931276e+05 2.278928e+03 1.613222e+03 7.097509e+03 2.278928e+03
131 132 133 134 135
1.737873e+03 7.447368e+04 2.278928e+03 -3.816806e+00 2.278928e+03
136 137 138 139 140
1.033206e+04 7.780053e+02 2.278928e+03 2.278928e+03 2.278928e+03
141 142 143 144 145
5.687380e+03 2.278928e+03 2.142707e+04 2.049192e+03 1.654492e+04
146 147 148 149 150
1.663148e+04 2.278928e+03 9.354033e+02 2.278928e+03 1.718739e+03
151 152 153 154 155
4.183272e+03 8.778797e+03 2.278928e+03 -1.150939e+04 1.330827e+03
156 157 158 159 160
2.278928e+03 2.278928e+03 1.652149e+03 1.940711e+04 2.278928e+03
161 162 163 164 165
2.278928e+03 2.278928e+03 2.278928e+03 8.492057e+03 -1.148714e+04
166 167 168 169 170
2.278928e+03 2.278928e+03 1.149630e+05 2.278928e+03 2.278928e+03
171 172 173 174 175
2.278928e+03 -1.540542e+03 2.278928e+03 2.147856e+03 2.278928e+03
176 177 178 179 180
6.766517e+02 2.278928e+03 4.637155e+04 2.278928e+03 2.278928e+03
181 182 183 184 185
3.362737e+04 3.777604e+04 2.421000e+04 1.386134e+03 4.849355e+04
186 187 188 189 190
2.747001e+04 4.769418e+04 1.676626e+03 -3.043527e+04 2.278928e+03
191 192 193 194 195
1.361331e+04 2.278928e+03 6.098086e+03 2.278928e+03 1.836639e+03
196 197 198 199 200
4.866357e+04 -5.518042e+04 4.866357e+04 4.866357e+04 -5.874979e+03
201 202 203 204 205
4.866357e+04 4.866357e+04 4.866357e+04 3.919355e+04 2.436529e+03
206 207 208 209 210
4.756263e+04 4.079204e+04 4.866357e+04 5.610607e+04 4.866357e+04
211 212 213 214 215
2.278928e+03 2.278928e+03 6.829477e+03 -1.506991e+04 -1.228798e+03
216 217 218 219 220
2.278928e+03 1.511224e+03 3.312338e+03 2.278928e+03 3.040236e+04
221 222 223 224 225
-2.296131e+04 2.278928e+03 2.060524e+03 2.278928e+03 2.105233e+03
226 227 228 229 230
1.199979e+03 2.278928e+03 2.278928e+03 -1.767343e+04 2.278928e+03
231 232 233 234 235
-1.191510e+04 2.278928e+03 2.278928e+03 6.583576e+03 2.278928e+03
236 237 238 239 240
2.007807e+03 2.278928e+03 3.293122e+03 2.278928e+03 2.278928e+03
241 242 243 244 245
2.278928e+03 6.901927e+03 2.278928e+03 3.077812e+03 -4.500862e+02
246 247 248 249 250
4.866357e+04 4.866357e+04 -7.510493e+03 4.866357e+04 4.866357e+04
251 252 253 254 255
7.919280e+03 4.866357e+04 -1.497148e+04 -3.668061e+04 -6.426448e+04
256 257 258 259 260
2.278928e+03 3.108864e+04 2.898284e+04 4.866357e+04 2.278928e+03
261 262 263 264 265
4.866357e+04 4.866357e+04 -1.182673e+04 1.171601e+04 4.866357e+04
266 267 268 269 270
2.278928e+03 4.866357e+04 1.242921e+04 4.866357e+04 -1.423921e+04
271 272 273 274 275
4.866357e+04 4.866357e+04 4.866357e+04 -2.092711e+03 2.278928e+03
276 277 278 279 280
3.683124e+03 2.210244e+03 2.278928e+03 -5.655267e+03 -2.685004e+04
281 282 283 284 285
4.866357e+04 2.278928e+03 4.774845e+03 2.278928e+03 -2.574312e+04
286 287 288 289 290
2.278928e+03 4.866357e+04 4.888437e+02 4.866357e+04 2.075239e+03
291 292 293 294 295
4.866357e+04 4.036418e+03 -1.347336e+04 2.278928e+03 -1.235231e+03
296 297 298 299 300
6.523467e+04 2.278928e+03 4.866357e+04 2.278928e+03 2.278928e+03
301 302 303 304 305
4.546553e+04 -3.678073e+03 2.278928e+03 4.866357e+04 2.220313e+04
306 307 308 309 310
9.847648e+02 4.866357e+04 3.802061e+04 2.180264e+04 4.866357e+04
311 312 313 314 315
4.866357e+04 -4.898008e+04 4.866357e+04 4.739741e+04 1.914256e+03
316 317 318 319 320
1.703036e+03 4.062286e+03 -9.101933e+03 1.424190e+03 2.278928e+03
321 322 323 324 325
2.278928e+03 2.278928e+03 2.278928e+03 4.866357e+04 2.278928e+03
326 327 328 329 330
-7.358770e+02 4.866357e+04 4.866357e+04 2.278928e+03 2.278928e+03
331 332 333 334 335
2.278928e+03 2.278928e+03 2.278928e+03 2.475569e+03 1.985279e+03
336 337 338 339 340
2.278928e+03 -2.146158e+02 2.278928e+03 2.008275e+03 2.363771e+03
341 342 343 344 345
-1.437134e+04 2.278928e+03 2.278928e+03 2.278928e+03 2.278928e+03
346 347 348 349 350
2.278928e+03 -5.911648e+03 2.278928e+03 2.278928e+03 2.278928e+03
351 352 353 354 355
4.546473e+03 -2.025177e+02 1.650416e+04 -3.517876e+03 2.278928e+03
356 357 358 359 360
2.278928e+03 -3.860831e+04 2.278928e+03 2.278928e+03 2.278928e+03
361 362 363 364 365
1.704788e+03 2.278928e+03 2.278928e+03 3.910634e+04 -5.914782e+00
366 367 368 369 370
2.278928e+03 2.278928e+03 2.278928e+03 2.278928e+03 2.278928e+03
371 372 373 374 375
2.278928e+03 -3.681370e+03 -2.850358e+04 -1.725417e+04 2.278928e+03
376 377 378 379 380
2.278928e+03 2.505398e+03 1.508541e+03 2.278928e+03 2.278928e+03
381 382 383 384 385
2.278928e+03 -1.871470e+04 2.192394e+04 2.278928e+03 -1.803844e+03
386 387 388 389 390
-8.902396e+04 7.428261e+03 1.289292e+04 1.015009e+03 -7.653258e+03
391 392 393 394 395
7.706070e+03 -7.972913e+03 -1.788966e+04 4.135504e+04 -1.467653e+05
396 397 398 399 400
-1.455957e+04 -1.816617e+04 8.697146e+04 -9.947973e+03 -1.256977e+05
401 402 403 404 405
2.623251e+03 1.525061e+04 -2.562332e+03 -4.190231e+04 -1.339545e+04
406 407 408 409 410
-1.232814e+04 1.477185e+03 8.848026e+03 -1.247915e+04 -4.548403e+04
411 412 413 414 415
-3.135735e+04 -2.085961e+05 -6.125829e+03 -1.058975e+05 -5.761806e+04
416 417 418 419 420
-5.485962e+04 1.507550e+05 -1.990883e+04 -1.697357e+04 -4.274209e+04
421 422 423 424 425
-8.140141e+04 -9.185457e+04 -6.273686e+05 -3.134823e+05 -6.283221e+04
426 427 428 429 430
-2.139888e+05 -2.420959e+05 -1.555632e+05 -1.122225e+05 1.590189e+05
431
-3.740122e+05
> postscript(file="/var/www/html/freestat/rcomp/tmp/6vdnj1291330748.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 -1.582183e+05 NA
1 2.967578e+06 -1.582183e+05
2 6.242201e+05 2.967578e+06
3 -5.016268e+05 6.242201e+05
4 -7.922385e+05 -5.016268e+05
5 -1.100347e+05 -7.922385e+05
6 1.040185e+06 -1.100347e+05
7 4.065020e+05 1.040185e+06
8 -6.808063e+05 4.065020e+05
9 6.323843e+05 -6.808063e+05
10 4.869825e+05 6.323843e+05
11 7.975685e+05 4.869825e+05
12 6.730946e+05 7.975685e+05
13 -4.528743e+05 6.730946e+05
14 -3.244176e+05 -4.528743e+05
15 5.698190e+05 -3.244176e+05
16 -2.426991e+05 5.698190e+05
17 -4.068701e+05 -2.426991e+05
18 -5.243865e+05 -4.068701e+05
19 -1.392642e+04 -5.243865e+05
20 1.140024e+05 -1.392642e+04
21 -3.499120e+05 1.140024e+05
22 1.331922e+05 -3.499120e+05
23 -1.707664e+05 1.331922e+05
24 -6.936383e+04 -1.707664e+05
25 3.353181e+05 -6.936383e+04
26 -6.224811e+05 3.353181e+05
27 -1.813264e+05 -6.224811e+05
28 6.187109e+05 -1.813264e+05
29 2.592913e+05 6.187109e+05
30 2.561397e+05 2.592913e+05
31 -1.477407e+05 2.561397e+05
32 1.626073e+05 -1.477407e+05
33 9.933395e+04 1.626073e+05
34 -2.360687e+05 9.933395e+04
35 5.011134e+04 -2.360687e+05
36 3.134921e+05 5.011134e+04
37 1.472344e+05 3.134921e+05
38 1.203332e+05 1.472344e+05
39 -1.480014e+05 1.203332e+05
40 -3.630277e+05 -1.480014e+05
41 4.559789e+04 -3.630277e+05
42 6.593958e+04 4.559789e+04
43 -5.337955e+04 6.593958e+04
44 6.053517e+04 -5.337955e+04
45 -1.268817e+05 6.053517e+04
46 8.914354e+04 -1.268817e+05
47 -2.587551e+05 8.914354e+04
48 -2.352834e+04 -2.587551e+05
49 -1.456213e+05 -2.352834e+04
50 8.620553e+04 -1.456213e+05
51 -5.728963e+04 8.620553e+04
52 -2.338111e+05 -5.728963e+04
53 5.006892e+04 -2.338111e+05
54 -2.989893e+05 5.006892e+04
55 1.045591e+05 -2.989893e+05
56 2.494654e+05 1.045591e+05
57 8.260581e+04 2.494654e+05
58 -8.837342e+04 8.260581e+04
59 -2.396914e+05 -8.837342e+04
60 -2.724862e+04 -2.396914e+05
61 -8.135870e+04 -2.724862e+04
62 -2.108917e+04 -8.135870e+04
63 -2.301377e+05 -2.108917e+04
64 5.552223e+04 -2.301377e+05
65 -1.438492e+05 5.552223e+04
66 -1.557565e+05 -1.438492e+05
67 -6.608753e+03 -1.557565e+05
68 -1.054951e+05 -6.608753e+03
69 -6.405175e+04 -1.054951e+05
70 -8.630760e+04 -6.405175e+04
71 1.613525e+05 -8.630760e+04
72 -1.342808e+05 1.613525e+05
73 -4.918722e+04 -1.342808e+05
74 -1.079814e+05 -4.918722e+04
75 -6.340797e+04 -1.079814e+05
76 -1.226305e+05 -6.340797e+04
77 -9.530678e+04 -1.226305e+05
78 7.467767e+03 -9.530678e+04
79 -8.680032e+03 7.467767e+03
80 2.989722e+04 -8.680032e+03
81 -5.606581e+04 2.989722e+04
82 -1.559318e+05 -5.606581e+04
83 4.414927e+04 -1.559318e+05
84 -3.436637e+04 4.414927e+04
85 -1.247111e+05 -3.436637e+04
86 5.766209e+04 -1.247111e+05
87 -6.611169e+04 5.766209e+04
88 -1.029301e+05 -6.611169e+04
89 -4.406551e+04 -1.029301e+05
90 -3.188368e+04 -4.406551e+04
91 -4.165470e+04 -3.188368e+04
92 1.008414e+05 -4.165470e+04
93 -2.423047e+05 1.008414e+05
94 1.187274e+05 -2.423047e+05
95 -7.177381e+03 1.187274e+05
96 -1.575220e+04 -7.177381e+03
97 -5.534011e+03 -1.575220e+04
98 -8.429364e+04 -5.534011e+03
99 -1.361739e+04 -8.429364e+04
100 -2.952767e+04 -1.361739e+04
101 -1.242537e+05 -2.952767e+04
102 -8.434649e+04 -1.242537e+05
103 -7.841049e+04 -8.434649e+04
104 3.679140e+04 -7.841049e+04
105 -8.245473e+04 3.679140e+04
106 -2.364768e+05 -8.245473e+04
107 1.103250e+05 -2.364768e+05
108 -9.009382e+04 1.103250e+05
109 -1.003881e+05 -9.009382e+04
110 7.883615e+04 -1.003881e+05
111 -2.488125e+05 7.883615e+04
112 -1.796693e+05 -2.488125e+05
113 -1.094369e+05 -1.796693e+05
114 3.351644e+04 -1.094369e+05
115 5.246518e+04 3.351644e+04
116 -3.566432e+04 5.246518e+04
117 6.295913e+04 -3.566432e+04
118 -1.660040e+05 6.295913e+04
119 1.065580e+04 -1.660040e+05
120 1.075239e+04 1.065580e+04
121 1.430733e+04 1.075239e+04
122 -6.019245e+04 1.430733e+04
123 7.868870e+04 -6.019245e+04
124 -9.422583e+04 7.868870e+04
125 -1.931276e+05 -9.422583e+04
126 2.278928e+03 -1.931276e+05
127 1.613222e+03 2.278928e+03
128 7.097509e+03 1.613222e+03
129 2.278928e+03 7.097509e+03
130 1.737873e+03 2.278928e+03
131 7.447368e+04 1.737873e+03
132 2.278928e+03 7.447368e+04
133 -3.816806e+00 2.278928e+03
134 2.278928e+03 -3.816806e+00
135 1.033206e+04 2.278928e+03
136 7.780053e+02 1.033206e+04
137 2.278928e+03 7.780053e+02
138 2.278928e+03 2.278928e+03
139 2.278928e+03 2.278928e+03
140 5.687380e+03 2.278928e+03
141 2.278928e+03 5.687380e+03
142 2.142707e+04 2.278928e+03
143 2.049192e+03 2.142707e+04
144 1.654492e+04 2.049192e+03
145 1.663148e+04 1.654492e+04
146 2.278928e+03 1.663148e+04
147 9.354033e+02 2.278928e+03
148 2.278928e+03 9.354033e+02
149 1.718739e+03 2.278928e+03
150 4.183272e+03 1.718739e+03
151 8.778797e+03 4.183272e+03
152 2.278928e+03 8.778797e+03
153 -1.150939e+04 2.278928e+03
154 1.330827e+03 -1.150939e+04
155 2.278928e+03 1.330827e+03
156 2.278928e+03 2.278928e+03
157 1.652149e+03 2.278928e+03
158 1.940711e+04 1.652149e+03
159 2.278928e+03 1.940711e+04
160 2.278928e+03 2.278928e+03
161 2.278928e+03 2.278928e+03
162 2.278928e+03 2.278928e+03
163 8.492057e+03 2.278928e+03
164 -1.148714e+04 8.492057e+03
165 2.278928e+03 -1.148714e+04
166 2.278928e+03 2.278928e+03
167 1.149630e+05 2.278928e+03
168 2.278928e+03 1.149630e+05
169 2.278928e+03 2.278928e+03
170 2.278928e+03 2.278928e+03
171 -1.540542e+03 2.278928e+03
172 2.278928e+03 -1.540542e+03
173 2.147856e+03 2.278928e+03
174 2.278928e+03 2.147856e+03
175 6.766517e+02 2.278928e+03
176 2.278928e+03 6.766517e+02
177 4.637155e+04 2.278928e+03
178 2.278928e+03 4.637155e+04
179 2.278928e+03 2.278928e+03
180 3.362737e+04 2.278928e+03
181 3.777604e+04 3.362737e+04
182 2.421000e+04 3.777604e+04
183 1.386134e+03 2.421000e+04
184 4.849355e+04 1.386134e+03
185 2.747001e+04 4.849355e+04
186 4.769418e+04 2.747001e+04
187 1.676626e+03 4.769418e+04
188 -3.043527e+04 1.676626e+03
189 2.278928e+03 -3.043527e+04
190 1.361331e+04 2.278928e+03
191 2.278928e+03 1.361331e+04
192 6.098086e+03 2.278928e+03
193 2.278928e+03 6.098086e+03
194 1.836639e+03 2.278928e+03
195 4.866357e+04 1.836639e+03
196 -5.518042e+04 4.866357e+04
197 4.866357e+04 -5.518042e+04
198 4.866357e+04 4.866357e+04
199 -5.874979e+03 4.866357e+04
200 4.866357e+04 -5.874979e+03
201 4.866357e+04 4.866357e+04
202 4.866357e+04 4.866357e+04
203 3.919355e+04 4.866357e+04
204 2.436529e+03 3.919355e+04
205 4.756263e+04 2.436529e+03
206 4.079204e+04 4.756263e+04
207 4.866357e+04 4.079204e+04
208 5.610607e+04 4.866357e+04
209 4.866357e+04 5.610607e+04
210 2.278928e+03 4.866357e+04
211 2.278928e+03 2.278928e+03
212 6.829477e+03 2.278928e+03
213 -1.506991e+04 6.829477e+03
214 -1.228798e+03 -1.506991e+04
215 2.278928e+03 -1.228798e+03
216 1.511224e+03 2.278928e+03
217 3.312338e+03 1.511224e+03
218 2.278928e+03 3.312338e+03
219 3.040236e+04 2.278928e+03
220 -2.296131e+04 3.040236e+04
221 2.278928e+03 -2.296131e+04
222 2.060524e+03 2.278928e+03
223 2.278928e+03 2.060524e+03
224 2.105233e+03 2.278928e+03
225 1.199979e+03 2.105233e+03
226 2.278928e+03 1.199979e+03
227 2.278928e+03 2.278928e+03
228 -1.767343e+04 2.278928e+03
229 2.278928e+03 -1.767343e+04
230 -1.191510e+04 2.278928e+03
231 2.278928e+03 -1.191510e+04
232 2.278928e+03 2.278928e+03
233 6.583576e+03 2.278928e+03
234 2.278928e+03 6.583576e+03
235 2.007807e+03 2.278928e+03
236 2.278928e+03 2.007807e+03
237 3.293122e+03 2.278928e+03
238 2.278928e+03 3.293122e+03
239 2.278928e+03 2.278928e+03
240 2.278928e+03 2.278928e+03
241 6.901927e+03 2.278928e+03
242 2.278928e+03 6.901927e+03
243 3.077812e+03 2.278928e+03
244 -4.500862e+02 3.077812e+03
245 4.866357e+04 -4.500862e+02
246 4.866357e+04 4.866357e+04
247 -7.510493e+03 4.866357e+04
248 4.866357e+04 -7.510493e+03
249 4.866357e+04 4.866357e+04
250 7.919280e+03 4.866357e+04
251 4.866357e+04 7.919280e+03
252 -1.497148e+04 4.866357e+04
253 -3.668061e+04 -1.497148e+04
254 -6.426448e+04 -3.668061e+04
255 2.278928e+03 -6.426448e+04
256 3.108864e+04 2.278928e+03
257 2.898284e+04 3.108864e+04
258 4.866357e+04 2.898284e+04
259 2.278928e+03 4.866357e+04
260 4.866357e+04 2.278928e+03
261 4.866357e+04 4.866357e+04
262 -1.182673e+04 4.866357e+04
263 1.171601e+04 -1.182673e+04
264 4.866357e+04 1.171601e+04
265 2.278928e+03 4.866357e+04
266 4.866357e+04 2.278928e+03
267 1.242921e+04 4.866357e+04
268 4.866357e+04 1.242921e+04
269 -1.423921e+04 4.866357e+04
270 4.866357e+04 -1.423921e+04
271 4.866357e+04 4.866357e+04
272 4.866357e+04 4.866357e+04
273 -2.092711e+03 4.866357e+04
274 2.278928e+03 -2.092711e+03
275 3.683124e+03 2.278928e+03
276 2.210244e+03 3.683124e+03
277 2.278928e+03 2.210244e+03
278 -5.655267e+03 2.278928e+03
279 -2.685004e+04 -5.655267e+03
280 4.866357e+04 -2.685004e+04
281 2.278928e+03 4.866357e+04
282 4.774845e+03 2.278928e+03
283 2.278928e+03 4.774845e+03
284 -2.574312e+04 2.278928e+03
285 2.278928e+03 -2.574312e+04
286 4.866357e+04 2.278928e+03
287 4.888437e+02 4.866357e+04
288 4.866357e+04 4.888437e+02
289 2.075239e+03 4.866357e+04
290 4.866357e+04 2.075239e+03
291 4.036418e+03 4.866357e+04
292 -1.347336e+04 4.036418e+03
293 2.278928e+03 -1.347336e+04
294 -1.235231e+03 2.278928e+03
295 6.523467e+04 -1.235231e+03
296 2.278928e+03 6.523467e+04
297 4.866357e+04 2.278928e+03
298 2.278928e+03 4.866357e+04
299 2.278928e+03 2.278928e+03
300 4.546553e+04 2.278928e+03
301 -3.678073e+03 4.546553e+04
302 2.278928e+03 -3.678073e+03
303 4.866357e+04 2.278928e+03
304 2.220313e+04 4.866357e+04
305 9.847648e+02 2.220313e+04
306 4.866357e+04 9.847648e+02
307 3.802061e+04 4.866357e+04
308 2.180264e+04 3.802061e+04
309 4.866357e+04 2.180264e+04
310 4.866357e+04 4.866357e+04
311 -4.898008e+04 4.866357e+04
312 4.866357e+04 -4.898008e+04
313 4.739741e+04 4.866357e+04
314 1.914256e+03 4.739741e+04
315 1.703036e+03 1.914256e+03
316 4.062286e+03 1.703036e+03
317 -9.101933e+03 4.062286e+03
318 1.424190e+03 -9.101933e+03
319 2.278928e+03 1.424190e+03
320 2.278928e+03 2.278928e+03
321 2.278928e+03 2.278928e+03
322 2.278928e+03 2.278928e+03
323 4.866357e+04 2.278928e+03
324 2.278928e+03 4.866357e+04
325 -7.358770e+02 2.278928e+03
326 4.866357e+04 -7.358770e+02
327 4.866357e+04 4.866357e+04
328 2.278928e+03 4.866357e+04
329 2.278928e+03 2.278928e+03
330 2.278928e+03 2.278928e+03
331 2.278928e+03 2.278928e+03
332 2.278928e+03 2.278928e+03
333 2.475569e+03 2.278928e+03
334 1.985279e+03 2.475569e+03
335 2.278928e+03 1.985279e+03
336 -2.146158e+02 2.278928e+03
337 2.278928e+03 -2.146158e+02
338 2.008275e+03 2.278928e+03
339 2.363771e+03 2.008275e+03
340 -1.437134e+04 2.363771e+03
341 2.278928e+03 -1.437134e+04
342 2.278928e+03 2.278928e+03
343 2.278928e+03 2.278928e+03
344 2.278928e+03 2.278928e+03
345 2.278928e+03 2.278928e+03
346 -5.911648e+03 2.278928e+03
347 2.278928e+03 -5.911648e+03
348 2.278928e+03 2.278928e+03
349 2.278928e+03 2.278928e+03
350 4.546473e+03 2.278928e+03
351 -2.025177e+02 4.546473e+03
352 1.650416e+04 -2.025177e+02
353 -3.517876e+03 1.650416e+04
354 2.278928e+03 -3.517876e+03
355 2.278928e+03 2.278928e+03
356 -3.860831e+04 2.278928e+03
357 2.278928e+03 -3.860831e+04
358 2.278928e+03 2.278928e+03
359 2.278928e+03 2.278928e+03
360 1.704788e+03 2.278928e+03
361 2.278928e+03 1.704788e+03
362 2.278928e+03 2.278928e+03
363 3.910634e+04 2.278928e+03
364 -5.914782e+00 3.910634e+04
365 2.278928e+03 -5.914782e+00
366 2.278928e+03 2.278928e+03
367 2.278928e+03 2.278928e+03
368 2.278928e+03 2.278928e+03
369 2.278928e+03 2.278928e+03
370 2.278928e+03 2.278928e+03
371 -3.681370e+03 2.278928e+03
372 -2.850358e+04 -3.681370e+03
373 -1.725417e+04 -2.850358e+04
374 2.278928e+03 -1.725417e+04
375 2.278928e+03 2.278928e+03
376 2.505398e+03 2.278928e+03
377 1.508541e+03 2.505398e+03
378 2.278928e+03 1.508541e+03
379 2.278928e+03 2.278928e+03
380 2.278928e+03 2.278928e+03
381 -1.871470e+04 2.278928e+03
382 2.192394e+04 -1.871470e+04
383 2.278928e+03 2.192394e+04
384 -1.803844e+03 2.278928e+03
385 -8.902396e+04 -1.803844e+03
386 7.428261e+03 -8.902396e+04
387 1.289292e+04 7.428261e+03
388 1.015009e+03 1.289292e+04
389 -7.653258e+03 1.015009e+03
390 7.706070e+03 -7.653258e+03
391 -7.972913e+03 7.706070e+03
392 -1.788966e+04 -7.972913e+03
393 4.135504e+04 -1.788966e+04
394 -1.467653e+05 4.135504e+04
395 -1.455957e+04 -1.467653e+05
396 -1.816617e+04 -1.455957e+04
397 8.697146e+04 -1.816617e+04
398 -9.947973e+03 8.697146e+04
399 -1.256977e+05 -9.947973e+03
400 2.623251e+03 -1.256977e+05
401 1.525061e+04 2.623251e+03
402 -2.562332e+03 1.525061e+04
403 -4.190231e+04 -2.562332e+03
404 -1.339545e+04 -4.190231e+04
405 -1.232814e+04 -1.339545e+04
406 1.477185e+03 -1.232814e+04
407 8.848026e+03 1.477185e+03
408 -1.247915e+04 8.848026e+03
409 -4.548403e+04 -1.247915e+04
410 -3.135735e+04 -4.548403e+04
411 -2.085961e+05 -3.135735e+04
412 -6.125829e+03 -2.085961e+05
413 -1.058975e+05 -6.125829e+03
414 -5.761806e+04 -1.058975e+05
415 -5.485962e+04 -5.761806e+04
416 1.507550e+05 -5.485962e+04
417 -1.990883e+04 1.507550e+05
418 -1.697357e+04 -1.990883e+04
419 -4.274209e+04 -1.697357e+04
420 -8.140141e+04 -4.274209e+04
421 -9.185457e+04 -8.140141e+04
422 -6.273686e+05 -9.185457e+04
423 -3.134823e+05 -6.273686e+05
424 -6.283221e+04 -3.134823e+05
425 -2.139888e+05 -6.283221e+04
426 -2.420959e+05 -2.139888e+05
427 -1.555632e+05 -2.420959e+05
428 -1.122225e+05 -1.555632e+05
429 1.590189e+05 -1.122225e+05
430 -3.740122e+05 1.590189e+05
431 NA -3.740122e+05
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2.967578e+06 -1.582183e+05
[2,] 6.242201e+05 2.967578e+06
[3,] -5.016268e+05 6.242201e+05
[4,] -7.922385e+05 -5.016268e+05
[5,] -1.100347e+05 -7.922385e+05
[6,] 1.040185e+06 -1.100347e+05
[7,] 4.065020e+05 1.040185e+06
[8,] -6.808063e+05 4.065020e+05
[9,] 6.323843e+05 -6.808063e+05
[10,] 4.869825e+05 6.323843e+05
[11,] 7.975685e+05 4.869825e+05
[12,] 6.730946e+05 7.975685e+05
[13,] -4.528743e+05 6.730946e+05
[14,] -3.244176e+05 -4.528743e+05
[15,] 5.698190e+05 -3.244176e+05
[16,] -2.426991e+05 5.698190e+05
[17,] -4.068701e+05 -2.426991e+05
[18,] -5.243865e+05 -4.068701e+05
[19,] -1.392642e+04 -5.243865e+05
[20,] 1.140024e+05 -1.392642e+04
[21,] -3.499120e+05 1.140024e+05
[22,] 1.331922e+05 -3.499120e+05
[23,] -1.707664e+05 1.331922e+05
[24,] -6.936383e+04 -1.707664e+05
[25,] 3.353181e+05 -6.936383e+04
[26,] -6.224811e+05 3.353181e+05
[27,] -1.813264e+05 -6.224811e+05
[28,] 6.187109e+05 -1.813264e+05
[29,] 2.592913e+05 6.187109e+05
[30,] 2.561397e+05 2.592913e+05
[31,] -1.477407e+05 2.561397e+05
[32,] 1.626073e+05 -1.477407e+05
[33,] 9.933395e+04 1.626073e+05
[34,] -2.360687e+05 9.933395e+04
[35,] 5.011134e+04 -2.360687e+05
[36,] 3.134921e+05 5.011134e+04
[37,] 1.472344e+05 3.134921e+05
[38,] 1.203332e+05 1.472344e+05
[39,] -1.480014e+05 1.203332e+05
[40,] -3.630277e+05 -1.480014e+05
[41,] 4.559789e+04 -3.630277e+05
[42,] 6.593958e+04 4.559789e+04
[43,] -5.337955e+04 6.593958e+04
[44,] 6.053517e+04 -5.337955e+04
[45,] -1.268817e+05 6.053517e+04
[46,] 8.914354e+04 -1.268817e+05
[47,] -2.587551e+05 8.914354e+04
[48,] -2.352834e+04 -2.587551e+05
[49,] -1.456213e+05 -2.352834e+04
[50,] 8.620553e+04 -1.456213e+05
[51,] -5.728963e+04 8.620553e+04
[52,] -2.338111e+05 -5.728963e+04
[53,] 5.006892e+04 -2.338111e+05
[54,] -2.989893e+05 5.006892e+04
[55,] 1.045591e+05 -2.989893e+05
[56,] 2.494654e+05 1.045591e+05
[57,] 8.260581e+04 2.494654e+05
[58,] -8.837342e+04 8.260581e+04
[59,] -2.396914e+05 -8.837342e+04
[60,] -2.724862e+04 -2.396914e+05
[61,] -8.135870e+04 -2.724862e+04
[62,] -2.108917e+04 -8.135870e+04
[63,] -2.301377e+05 -2.108917e+04
[64,] 5.552223e+04 -2.301377e+05
[65,] -1.438492e+05 5.552223e+04
[66,] -1.557565e+05 -1.438492e+05
[67,] -6.608753e+03 -1.557565e+05
[68,] -1.054951e+05 -6.608753e+03
[69,] -6.405175e+04 -1.054951e+05
[70,] -8.630760e+04 -6.405175e+04
[71,] 1.613525e+05 -8.630760e+04
[72,] -1.342808e+05 1.613525e+05
[73,] -4.918722e+04 -1.342808e+05
[74,] -1.079814e+05 -4.918722e+04
[75,] -6.340797e+04 -1.079814e+05
[76,] -1.226305e+05 -6.340797e+04
[77,] -9.530678e+04 -1.226305e+05
[78,] 7.467767e+03 -9.530678e+04
[79,] -8.680032e+03 7.467767e+03
[80,] 2.989722e+04 -8.680032e+03
[81,] -5.606581e+04 2.989722e+04
[82,] -1.559318e+05 -5.606581e+04
[83,] 4.414927e+04 -1.559318e+05
[84,] -3.436637e+04 4.414927e+04
[85,] -1.247111e+05 -3.436637e+04
[86,] 5.766209e+04 -1.247111e+05
[87,] -6.611169e+04 5.766209e+04
[88,] -1.029301e+05 -6.611169e+04
[89,] -4.406551e+04 -1.029301e+05
[90,] -3.188368e+04 -4.406551e+04
[91,] -4.165470e+04 -3.188368e+04
[92,] 1.008414e+05 -4.165470e+04
[93,] -2.423047e+05 1.008414e+05
[94,] 1.187274e+05 -2.423047e+05
[95,] -7.177381e+03 1.187274e+05
[96,] -1.575220e+04 -7.177381e+03
[97,] -5.534011e+03 -1.575220e+04
[98,] -8.429364e+04 -5.534011e+03
[99,] -1.361739e+04 -8.429364e+04
[100,] -2.952767e+04 -1.361739e+04
[101,] -1.242537e+05 -2.952767e+04
[102,] -8.434649e+04 -1.242537e+05
[103,] -7.841049e+04 -8.434649e+04
[104,] 3.679140e+04 -7.841049e+04
[105,] -8.245473e+04 3.679140e+04
[106,] -2.364768e+05 -8.245473e+04
[107,] 1.103250e+05 -2.364768e+05
[108,] -9.009382e+04 1.103250e+05
[109,] -1.003881e+05 -9.009382e+04
[110,] 7.883615e+04 -1.003881e+05
[111,] -2.488125e+05 7.883615e+04
[112,] -1.796693e+05 -2.488125e+05
[113,] -1.094369e+05 -1.796693e+05
[114,] 3.351644e+04 -1.094369e+05
[115,] 5.246518e+04 3.351644e+04
[116,] -3.566432e+04 5.246518e+04
[117,] 6.295913e+04 -3.566432e+04
[118,] -1.660040e+05 6.295913e+04
[119,] 1.065580e+04 -1.660040e+05
[120,] 1.075239e+04 1.065580e+04
[121,] 1.430733e+04 1.075239e+04
[122,] -6.019245e+04 1.430733e+04
[123,] 7.868870e+04 -6.019245e+04
[124,] -9.422583e+04 7.868870e+04
[125,] -1.931276e+05 -9.422583e+04
[126,] 2.278928e+03 -1.931276e+05
[127,] 1.613222e+03 2.278928e+03
[128,] 7.097509e+03 1.613222e+03
[129,] 2.278928e+03 7.097509e+03
[130,] 1.737873e+03 2.278928e+03
[131,] 7.447368e+04 1.737873e+03
[132,] 2.278928e+03 7.447368e+04
[133,] -3.816806e+00 2.278928e+03
[134,] 2.278928e+03 -3.816806e+00
[135,] 1.033206e+04 2.278928e+03
[136,] 7.780053e+02 1.033206e+04
[137,] 2.278928e+03 7.780053e+02
[138,] 2.278928e+03 2.278928e+03
[139,] 2.278928e+03 2.278928e+03
[140,] 5.687380e+03 2.278928e+03
[141,] 2.278928e+03 5.687380e+03
[142,] 2.142707e+04 2.278928e+03
[143,] 2.049192e+03 2.142707e+04
[144,] 1.654492e+04 2.049192e+03
[145,] 1.663148e+04 1.654492e+04
[146,] 2.278928e+03 1.663148e+04
[147,] 9.354033e+02 2.278928e+03
[148,] 2.278928e+03 9.354033e+02
[149,] 1.718739e+03 2.278928e+03
[150,] 4.183272e+03 1.718739e+03
[151,] 8.778797e+03 4.183272e+03
[152,] 2.278928e+03 8.778797e+03
[153,] -1.150939e+04 2.278928e+03
[154,] 1.330827e+03 -1.150939e+04
[155,] 2.278928e+03 1.330827e+03
[156,] 2.278928e+03 2.278928e+03
[157,] 1.652149e+03 2.278928e+03
[158,] 1.940711e+04 1.652149e+03
[159,] 2.278928e+03 1.940711e+04
[160,] 2.278928e+03 2.278928e+03
[161,] 2.278928e+03 2.278928e+03
[162,] 2.278928e+03 2.278928e+03
[163,] 8.492057e+03 2.278928e+03
[164,] -1.148714e+04 8.492057e+03
[165,] 2.278928e+03 -1.148714e+04
[166,] 2.278928e+03 2.278928e+03
[167,] 1.149630e+05 2.278928e+03
[168,] 2.278928e+03 1.149630e+05
[169,] 2.278928e+03 2.278928e+03
[170,] 2.278928e+03 2.278928e+03
[171,] -1.540542e+03 2.278928e+03
[172,] 2.278928e+03 -1.540542e+03
[173,] 2.147856e+03 2.278928e+03
[174,] 2.278928e+03 2.147856e+03
[175,] 6.766517e+02 2.278928e+03
[176,] 2.278928e+03 6.766517e+02
[177,] 4.637155e+04 2.278928e+03
[178,] 2.278928e+03 4.637155e+04
[179,] 2.278928e+03 2.278928e+03
[180,] 3.362737e+04 2.278928e+03
[181,] 3.777604e+04 3.362737e+04
[182,] 2.421000e+04 3.777604e+04
[183,] 1.386134e+03 2.421000e+04
[184,] 4.849355e+04 1.386134e+03
[185,] 2.747001e+04 4.849355e+04
[186,] 4.769418e+04 2.747001e+04
[187,] 1.676626e+03 4.769418e+04
[188,] -3.043527e+04 1.676626e+03
[189,] 2.278928e+03 -3.043527e+04
[190,] 1.361331e+04 2.278928e+03
[191,] 2.278928e+03 1.361331e+04
[192,] 6.098086e+03 2.278928e+03
[193,] 2.278928e+03 6.098086e+03
[194,] 1.836639e+03 2.278928e+03
[195,] 4.866357e+04 1.836639e+03
[196,] -5.518042e+04 4.866357e+04
[197,] 4.866357e+04 -5.518042e+04
[198,] 4.866357e+04 4.866357e+04
[199,] -5.874979e+03 4.866357e+04
[200,] 4.866357e+04 -5.874979e+03
[201,] 4.866357e+04 4.866357e+04
[202,] 4.866357e+04 4.866357e+04
[203,] 3.919355e+04 4.866357e+04
[204,] 2.436529e+03 3.919355e+04
[205,] 4.756263e+04 2.436529e+03
[206,] 4.079204e+04 4.756263e+04
[207,] 4.866357e+04 4.079204e+04
[208,] 5.610607e+04 4.866357e+04
[209,] 4.866357e+04 5.610607e+04
[210,] 2.278928e+03 4.866357e+04
[211,] 2.278928e+03 2.278928e+03
[212,] 6.829477e+03 2.278928e+03
[213,] -1.506991e+04 6.829477e+03
[214,] -1.228798e+03 -1.506991e+04
[215,] 2.278928e+03 -1.228798e+03
[216,] 1.511224e+03 2.278928e+03
[217,] 3.312338e+03 1.511224e+03
[218,] 2.278928e+03 3.312338e+03
[219,] 3.040236e+04 2.278928e+03
[220,] -2.296131e+04 3.040236e+04
[221,] 2.278928e+03 -2.296131e+04
[222,] 2.060524e+03 2.278928e+03
[223,] 2.278928e+03 2.060524e+03
[224,] 2.105233e+03 2.278928e+03
[225,] 1.199979e+03 2.105233e+03
[226,] 2.278928e+03 1.199979e+03
[227,] 2.278928e+03 2.278928e+03
[228,] -1.767343e+04 2.278928e+03
[229,] 2.278928e+03 -1.767343e+04
[230,] -1.191510e+04 2.278928e+03
[231,] 2.278928e+03 -1.191510e+04
[232,] 2.278928e+03 2.278928e+03
[233,] 6.583576e+03 2.278928e+03
[234,] 2.278928e+03 6.583576e+03
[235,] 2.007807e+03 2.278928e+03
[236,] 2.278928e+03 2.007807e+03
[237,] 3.293122e+03 2.278928e+03
[238,] 2.278928e+03 3.293122e+03
[239,] 2.278928e+03 2.278928e+03
[240,] 2.278928e+03 2.278928e+03
[241,] 6.901927e+03 2.278928e+03
[242,] 2.278928e+03 6.901927e+03
[243,] 3.077812e+03 2.278928e+03
[244,] -4.500862e+02 3.077812e+03
[245,] 4.866357e+04 -4.500862e+02
[246,] 4.866357e+04 4.866357e+04
[247,] -7.510493e+03 4.866357e+04
[248,] 4.866357e+04 -7.510493e+03
[249,] 4.866357e+04 4.866357e+04
[250,] 7.919280e+03 4.866357e+04
[251,] 4.866357e+04 7.919280e+03
[252,] -1.497148e+04 4.866357e+04
[253,] -3.668061e+04 -1.497148e+04
[254,] -6.426448e+04 -3.668061e+04
[255,] 2.278928e+03 -6.426448e+04
[256,] 3.108864e+04 2.278928e+03
[257,] 2.898284e+04 3.108864e+04
[258,] 4.866357e+04 2.898284e+04
[259,] 2.278928e+03 4.866357e+04
[260,] 4.866357e+04 2.278928e+03
[261,] 4.866357e+04 4.866357e+04
[262,] -1.182673e+04 4.866357e+04
[263,] 1.171601e+04 -1.182673e+04
[264,] 4.866357e+04 1.171601e+04
[265,] 2.278928e+03 4.866357e+04
[266,] 4.866357e+04 2.278928e+03
[267,] 1.242921e+04 4.866357e+04
[268,] 4.866357e+04 1.242921e+04
[269,] -1.423921e+04 4.866357e+04
[270,] 4.866357e+04 -1.423921e+04
[271,] 4.866357e+04 4.866357e+04
[272,] 4.866357e+04 4.866357e+04
[273,] -2.092711e+03 4.866357e+04
[274,] 2.278928e+03 -2.092711e+03
[275,] 3.683124e+03 2.278928e+03
[276,] 2.210244e+03 3.683124e+03
[277,] 2.278928e+03 2.210244e+03
[278,] -5.655267e+03 2.278928e+03
[279,] -2.685004e+04 -5.655267e+03
[280,] 4.866357e+04 -2.685004e+04
[281,] 2.278928e+03 4.866357e+04
[282,] 4.774845e+03 2.278928e+03
[283,] 2.278928e+03 4.774845e+03
[284,] -2.574312e+04 2.278928e+03
[285,] 2.278928e+03 -2.574312e+04
[286,] 4.866357e+04 2.278928e+03
[287,] 4.888437e+02 4.866357e+04
[288,] 4.866357e+04 4.888437e+02
[289,] 2.075239e+03 4.866357e+04
[290,] 4.866357e+04 2.075239e+03
[291,] 4.036418e+03 4.866357e+04
[292,] -1.347336e+04 4.036418e+03
[293,] 2.278928e+03 -1.347336e+04
[294,] -1.235231e+03 2.278928e+03
[295,] 6.523467e+04 -1.235231e+03
[296,] 2.278928e+03 6.523467e+04
[297,] 4.866357e+04 2.278928e+03
[298,] 2.278928e+03 4.866357e+04
[299,] 2.278928e+03 2.278928e+03
[300,] 4.546553e+04 2.278928e+03
[301,] -3.678073e+03 4.546553e+04
[302,] 2.278928e+03 -3.678073e+03
[303,] 4.866357e+04 2.278928e+03
[304,] 2.220313e+04 4.866357e+04
[305,] 9.847648e+02 2.220313e+04
[306,] 4.866357e+04 9.847648e+02
[307,] 3.802061e+04 4.866357e+04
[308,] 2.180264e+04 3.802061e+04
[309,] 4.866357e+04 2.180264e+04
[310,] 4.866357e+04 4.866357e+04
[311,] -4.898008e+04 4.866357e+04
[312,] 4.866357e+04 -4.898008e+04
[313,] 4.739741e+04 4.866357e+04
[314,] 1.914256e+03 4.739741e+04
[315,] 1.703036e+03 1.914256e+03
[316,] 4.062286e+03 1.703036e+03
[317,] -9.101933e+03 4.062286e+03
[318,] 1.424190e+03 -9.101933e+03
[319,] 2.278928e+03 1.424190e+03
[320,] 2.278928e+03 2.278928e+03
[321,] 2.278928e+03 2.278928e+03
[322,] 2.278928e+03 2.278928e+03
[323,] 4.866357e+04 2.278928e+03
[324,] 2.278928e+03 4.866357e+04
[325,] -7.358770e+02 2.278928e+03
[326,] 4.866357e+04 -7.358770e+02
[327,] 4.866357e+04 4.866357e+04
[328,] 2.278928e+03 4.866357e+04
[329,] 2.278928e+03 2.278928e+03
[330,] 2.278928e+03 2.278928e+03
[331,] 2.278928e+03 2.278928e+03
[332,] 2.278928e+03 2.278928e+03
[333,] 2.475569e+03 2.278928e+03
[334,] 1.985279e+03 2.475569e+03
[335,] 2.278928e+03 1.985279e+03
[336,] -2.146158e+02 2.278928e+03
[337,] 2.278928e+03 -2.146158e+02
[338,] 2.008275e+03 2.278928e+03
[339,] 2.363771e+03 2.008275e+03
[340,] -1.437134e+04 2.363771e+03
[341,] 2.278928e+03 -1.437134e+04
[342,] 2.278928e+03 2.278928e+03
[343,] 2.278928e+03 2.278928e+03
[344,] 2.278928e+03 2.278928e+03
[345,] 2.278928e+03 2.278928e+03
[346,] -5.911648e+03 2.278928e+03
[347,] 2.278928e+03 -5.911648e+03
[348,] 2.278928e+03 2.278928e+03
[349,] 2.278928e+03 2.278928e+03
[350,] 4.546473e+03 2.278928e+03
[351,] -2.025177e+02 4.546473e+03
[352,] 1.650416e+04 -2.025177e+02
[353,] -3.517876e+03 1.650416e+04
[354,] 2.278928e+03 -3.517876e+03
[355,] 2.278928e+03 2.278928e+03
[356,] -3.860831e+04 2.278928e+03
[357,] 2.278928e+03 -3.860831e+04
[358,] 2.278928e+03 2.278928e+03
[359,] 2.278928e+03 2.278928e+03
[360,] 1.704788e+03 2.278928e+03
[361,] 2.278928e+03 1.704788e+03
[362,] 2.278928e+03 2.278928e+03
[363,] 3.910634e+04 2.278928e+03
[364,] -5.914782e+00 3.910634e+04
[365,] 2.278928e+03 -5.914782e+00
[366,] 2.278928e+03 2.278928e+03
[367,] 2.278928e+03 2.278928e+03
[368,] 2.278928e+03 2.278928e+03
[369,] 2.278928e+03 2.278928e+03
[370,] 2.278928e+03 2.278928e+03
[371,] -3.681370e+03 2.278928e+03
[372,] -2.850358e+04 -3.681370e+03
[373,] -1.725417e+04 -2.850358e+04
[374,] 2.278928e+03 -1.725417e+04
[375,] 2.278928e+03 2.278928e+03
[376,] 2.505398e+03 2.278928e+03
[377,] 1.508541e+03 2.505398e+03
[378,] 2.278928e+03 1.508541e+03
[379,] 2.278928e+03 2.278928e+03
[380,] 2.278928e+03 2.278928e+03
[381,] -1.871470e+04 2.278928e+03
[382,] 2.192394e+04 -1.871470e+04
[383,] 2.278928e+03 2.192394e+04
[384,] -1.803844e+03 2.278928e+03
[385,] -8.902396e+04 -1.803844e+03
[386,] 7.428261e+03 -8.902396e+04
[387,] 1.289292e+04 7.428261e+03
[388,] 1.015009e+03 1.289292e+04
[389,] -7.653258e+03 1.015009e+03
[390,] 7.706070e+03 -7.653258e+03
[391,] -7.972913e+03 7.706070e+03
[392,] -1.788966e+04 -7.972913e+03
[393,] 4.135504e+04 -1.788966e+04
[394,] -1.467653e+05 4.135504e+04
[395,] -1.455957e+04 -1.467653e+05
[396,] -1.816617e+04 -1.455957e+04
[397,] 8.697146e+04 -1.816617e+04
[398,] -9.947973e+03 8.697146e+04
[399,] -1.256977e+05 -9.947973e+03
[400,] 2.623251e+03 -1.256977e+05
[401,] 1.525061e+04 2.623251e+03
[402,] -2.562332e+03 1.525061e+04
[403,] -4.190231e+04 -2.562332e+03
[404,] -1.339545e+04 -4.190231e+04
[405,] -1.232814e+04 -1.339545e+04
[406,] 1.477185e+03 -1.232814e+04
[407,] 8.848026e+03 1.477185e+03
[408,] -1.247915e+04 8.848026e+03
[409,] -4.548403e+04 -1.247915e+04
[410,] -3.135735e+04 -4.548403e+04
[411,] -2.085961e+05 -3.135735e+04
[412,] -6.125829e+03 -2.085961e+05
[413,] -1.058975e+05 -6.125829e+03
[414,] -5.761806e+04 -1.058975e+05
[415,] -5.485962e+04 -5.761806e+04
[416,] 1.507550e+05 -5.485962e+04
[417,] -1.990883e+04 1.507550e+05
[418,] -1.697357e+04 -1.990883e+04
[419,] -4.274209e+04 -1.697357e+04
[420,] -8.140141e+04 -4.274209e+04
[421,] -9.185457e+04 -8.140141e+04
[422,] -6.273686e+05 -9.185457e+04
[423,] -3.134823e+05 -6.273686e+05
[424,] -6.283221e+04 -3.134823e+05
[425,] -2.139888e+05 -6.283221e+04
[426,] -2.420959e+05 -2.139888e+05
[427,] -1.555632e+05 -2.420959e+05
[428,] -1.122225e+05 -1.555632e+05
[429,] 1.590189e+05 -1.122225e+05
[430,] -3.740122e+05 1.590189e+05
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2.967578e+06 -1.582183e+05
2 6.242201e+05 2.967578e+06
3 -5.016268e+05 6.242201e+05
4 -7.922385e+05 -5.016268e+05
5 -1.100347e+05 -7.922385e+05
6 1.040185e+06 -1.100347e+05
7 4.065020e+05 1.040185e+06
8 -6.808063e+05 4.065020e+05
9 6.323843e+05 -6.808063e+05
10 4.869825e+05 6.323843e+05
11 7.975685e+05 4.869825e+05
12 6.730946e+05 7.975685e+05
13 -4.528743e+05 6.730946e+05
14 -3.244176e+05 -4.528743e+05
15 5.698190e+05 -3.244176e+05
16 -2.426991e+05 5.698190e+05
17 -4.068701e+05 -2.426991e+05
18 -5.243865e+05 -4.068701e+05
19 -1.392642e+04 -5.243865e+05
20 1.140024e+05 -1.392642e+04
21 -3.499120e+05 1.140024e+05
22 1.331922e+05 -3.499120e+05
23 -1.707664e+05 1.331922e+05
24 -6.936383e+04 -1.707664e+05
25 3.353181e+05 -6.936383e+04
26 -6.224811e+05 3.353181e+05
27 -1.813264e+05 -6.224811e+05
28 6.187109e+05 -1.813264e+05
29 2.592913e+05 6.187109e+05
30 2.561397e+05 2.592913e+05
31 -1.477407e+05 2.561397e+05
32 1.626073e+05 -1.477407e+05
33 9.933395e+04 1.626073e+05
34 -2.360687e+05 9.933395e+04
35 5.011134e+04 -2.360687e+05
36 3.134921e+05 5.011134e+04
37 1.472344e+05 3.134921e+05
38 1.203332e+05 1.472344e+05
39 -1.480014e+05 1.203332e+05
40 -3.630277e+05 -1.480014e+05
41 4.559789e+04 -3.630277e+05
42 6.593958e+04 4.559789e+04
43 -5.337955e+04 6.593958e+04
44 6.053517e+04 -5.337955e+04
45 -1.268817e+05 6.053517e+04
46 8.914354e+04 -1.268817e+05
47 -2.587551e+05 8.914354e+04
48 -2.352834e+04 -2.587551e+05
49 -1.456213e+05 -2.352834e+04
50 8.620553e+04 -1.456213e+05
51 -5.728963e+04 8.620553e+04
52 -2.338111e+05 -5.728963e+04
53 5.006892e+04 -2.338111e+05
54 -2.989893e+05 5.006892e+04
55 1.045591e+05 -2.989893e+05
56 2.494654e+05 1.045591e+05
57 8.260581e+04 2.494654e+05
58 -8.837342e+04 8.260581e+04
59 -2.396914e+05 -8.837342e+04
60 -2.724862e+04 -2.396914e+05
61 -8.135870e+04 -2.724862e+04
62 -2.108917e+04 -8.135870e+04
63 -2.301377e+05 -2.108917e+04
64 5.552223e+04 -2.301377e+05
65 -1.438492e+05 5.552223e+04
66 -1.557565e+05 -1.438492e+05
67 -6.608753e+03 -1.557565e+05
68 -1.054951e+05 -6.608753e+03
69 -6.405175e+04 -1.054951e+05
70 -8.630760e+04 -6.405175e+04
71 1.613525e+05 -8.630760e+04
72 -1.342808e+05 1.613525e+05
73 -4.918722e+04 -1.342808e+05
74 -1.079814e+05 -4.918722e+04
75 -6.340797e+04 -1.079814e+05
76 -1.226305e+05 -6.340797e+04
77 -9.530678e+04 -1.226305e+05
78 7.467767e+03 -9.530678e+04
79 -8.680032e+03 7.467767e+03
80 2.989722e+04 -8.680032e+03
81 -5.606581e+04 2.989722e+04
82 -1.559318e+05 -5.606581e+04
83 4.414927e+04 -1.559318e+05
84 -3.436637e+04 4.414927e+04
85 -1.247111e+05 -3.436637e+04
86 5.766209e+04 -1.247111e+05
87 -6.611169e+04 5.766209e+04
88 -1.029301e+05 -6.611169e+04
89 -4.406551e+04 -1.029301e+05
90 -3.188368e+04 -4.406551e+04
91 -4.165470e+04 -3.188368e+04
92 1.008414e+05 -4.165470e+04
93 -2.423047e+05 1.008414e+05
94 1.187274e+05 -2.423047e+05
95 -7.177381e+03 1.187274e+05
96 -1.575220e+04 -7.177381e+03
97 -5.534011e+03 -1.575220e+04
98 -8.429364e+04 -5.534011e+03
99 -1.361739e+04 -8.429364e+04
100 -2.952767e+04 -1.361739e+04
101 -1.242537e+05 -2.952767e+04
102 -8.434649e+04 -1.242537e+05
103 -7.841049e+04 -8.434649e+04
104 3.679140e+04 -7.841049e+04
105 -8.245473e+04 3.679140e+04
106 -2.364768e+05 -8.245473e+04
107 1.103250e+05 -2.364768e+05
108 -9.009382e+04 1.103250e+05
109 -1.003881e+05 -9.009382e+04
110 7.883615e+04 -1.003881e+05
111 -2.488125e+05 7.883615e+04
112 -1.796693e+05 -2.488125e+05
113 -1.094369e+05 -1.796693e+05
114 3.351644e+04 -1.094369e+05
115 5.246518e+04 3.351644e+04
116 -3.566432e+04 5.246518e+04
117 6.295913e+04 -3.566432e+04
118 -1.660040e+05 6.295913e+04
119 1.065580e+04 -1.660040e+05
120 1.075239e+04 1.065580e+04
121 1.430733e+04 1.075239e+04
122 -6.019245e+04 1.430733e+04
123 7.868870e+04 -6.019245e+04
124 -9.422583e+04 7.868870e+04
125 -1.931276e+05 -9.422583e+04
126 2.278928e+03 -1.931276e+05
127 1.613222e+03 2.278928e+03
128 7.097509e+03 1.613222e+03
129 2.278928e+03 7.097509e+03
130 1.737873e+03 2.278928e+03
131 7.447368e+04 1.737873e+03
132 2.278928e+03 7.447368e+04
133 -3.816806e+00 2.278928e+03
134 2.278928e+03 -3.816806e+00
135 1.033206e+04 2.278928e+03
136 7.780053e+02 1.033206e+04
137 2.278928e+03 7.780053e+02
138 2.278928e+03 2.278928e+03
139 2.278928e+03 2.278928e+03
140 5.687380e+03 2.278928e+03
141 2.278928e+03 5.687380e+03
142 2.142707e+04 2.278928e+03
143 2.049192e+03 2.142707e+04
144 1.654492e+04 2.049192e+03
145 1.663148e+04 1.654492e+04
146 2.278928e+03 1.663148e+04
147 9.354033e+02 2.278928e+03
148 2.278928e+03 9.354033e+02
149 1.718739e+03 2.278928e+03
150 4.183272e+03 1.718739e+03
151 8.778797e+03 4.183272e+03
152 2.278928e+03 8.778797e+03
153 -1.150939e+04 2.278928e+03
154 1.330827e+03 -1.150939e+04
155 2.278928e+03 1.330827e+03
156 2.278928e+03 2.278928e+03
157 1.652149e+03 2.278928e+03
158 1.940711e+04 1.652149e+03
159 2.278928e+03 1.940711e+04
160 2.278928e+03 2.278928e+03
161 2.278928e+03 2.278928e+03
162 2.278928e+03 2.278928e+03
163 8.492057e+03 2.278928e+03
164 -1.148714e+04 8.492057e+03
165 2.278928e+03 -1.148714e+04
166 2.278928e+03 2.278928e+03
167 1.149630e+05 2.278928e+03
168 2.278928e+03 1.149630e+05
169 2.278928e+03 2.278928e+03
170 2.278928e+03 2.278928e+03
171 -1.540542e+03 2.278928e+03
172 2.278928e+03 -1.540542e+03
173 2.147856e+03 2.278928e+03
174 2.278928e+03 2.147856e+03
175 6.766517e+02 2.278928e+03
176 2.278928e+03 6.766517e+02
177 4.637155e+04 2.278928e+03
178 2.278928e+03 4.637155e+04
179 2.278928e+03 2.278928e+03
180 3.362737e+04 2.278928e+03
181 3.777604e+04 3.362737e+04
182 2.421000e+04 3.777604e+04
183 1.386134e+03 2.421000e+04
184 4.849355e+04 1.386134e+03
185 2.747001e+04 4.849355e+04
186 4.769418e+04 2.747001e+04
187 1.676626e+03 4.769418e+04
188 -3.043527e+04 1.676626e+03
189 2.278928e+03 -3.043527e+04
190 1.361331e+04 2.278928e+03
191 2.278928e+03 1.361331e+04
192 6.098086e+03 2.278928e+03
193 2.278928e+03 6.098086e+03
194 1.836639e+03 2.278928e+03
195 4.866357e+04 1.836639e+03
196 -5.518042e+04 4.866357e+04
197 4.866357e+04 -5.518042e+04
198 4.866357e+04 4.866357e+04
199 -5.874979e+03 4.866357e+04
200 4.866357e+04 -5.874979e+03
201 4.866357e+04 4.866357e+04
202 4.866357e+04 4.866357e+04
203 3.919355e+04 4.866357e+04
204 2.436529e+03 3.919355e+04
205 4.756263e+04 2.436529e+03
206 4.079204e+04 4.756263e+04
207 4.866357e+04 4.079204e+04
208 5.610607e+04 4.866357e+04
209 4.866357e+04 5.610607e+04
210 2.278928e+03 4.866357e+04
211 2.278928e+03 2.278928e+03
212 6.829477e+03 2.278928e+03
213 -1.506991e+04 6.829477e+03
214 -1.228798e+03 -1.506991e+04
215 2.278928e+03 -1.228798e+03
216 1.511224e+03 2.278928e+03
217 3.312338e+03 1.511224e+03
218 2.278928e+03 3.312338e+03
219 3.040236e+04 2.278928e+03
220 -2.296131e+04 3.040236e+04
221 2.278928e+03 -2.296131e+04
222 2.060524e+03 2.278928e+03
223 2.278928e+03 2.060524e+03
224 2.105233e+03 2.278928e+03
225 1.199979e+03 2.105233e+03
226 2.278928e+03 1.199979e+03
227 2.278928e+03 2.278928e+03
228 -1.767343e+04 2.278928e+03
229 2.278928e+03 -1.767343e+04
230 -1.191510e+04 2.278928e+03
231 2.278928e+03 -1.191510e+04
232 2.278928e+03 2.278928e+03
233 6.583576e+03 2.278928e+03
234 2.278928e+03 6.583576e+03
235 2.007807e+03 2.278928e+03
236 2.278928e+03 2.007807e+03
237 3.293122e+03 2.278928e+03
238 2.278928e+03 3.293122e+03
239 2.278928e+03 2.278928e+03
240 2.278928e+03 2.278928e+03
241 6.901927e+03 2.278928e+03
242 2.278928e+03 6.901927e+03
243 3.077812e+03 2.278928e+03
244 -4.500862e+02 3.077812e+03
245 4.866357e+04 -4.500862e+02
246 4.866357e+04 4.866357e+04
247 -7.510493e+03 4.866357e+04
248 4.866357e+04 -7.510493e+03
249 4.866357e+04 4.866357e+04
250 7.919280e+03 4.866357e+04
251 4.866357e+04 7.919280e+03
252 -1.497148e+04 4.866357e+04
253 -3.668061e+04 -1.497148e+04
254 -6.426448e+04 -3.668061e+04
255 2.278928e+03 -6.426448e+04
256 3.108864e+04 2.278928e+03
257 2.898284e+04 3.108864e+04
258 4.866357e+04 2.898284e+04
259 2.278928e+03 4.866357e+04
260 4.866357e+04 2.278928e+03
261 4.866357e+04 4.866357e+04
262 -1.182673e+04 4.866357e+04
263 1.171601e+04 -1.182673e+04
264 4.866357e+04 1.171601e+04
265 2.278928e+03 4.866357e+04
266 4.866357e+04 2.278928e+03
267 1.242921e+04 4.866357e+04
268 4.866357e+04 1.242921e+04
269 -1.423921e+04 4.866357e+04
270 4.866357e+04 -1.423921e+04
271 4.866357e+04 4.866357e+04
272 4.866357e+04 4.866357e+04
273 -2.092711e+03 4.866357e+04
274 2.278928e+03 -2.092711e+03
275 3.683124e+03 2.278928e+03
276 2.210244e+03 3.683124e+03
277 2.278928e+03 2.210244e+03
278 -5.655267e+03 2.278928e+03
279 -2.685004e+04 -5.655267e+03
280 4.866357e+04 -2.685004e+04
281 2.278928e+03 4.866357e+04
282 4.774845e+03 2.278928e+03
283 2.278928e+03 4.774845e+03
284 -2.574312e+04 2.278928e+03
285 2.278928e+03 -2.574312e+04
286 4.866357e+04 2.278928e+03
287 4.888437e+02 4.866357e+04
288 4.866357e+04 4.888437e+02
289 2.075239e+03 4.866357e+04
290 4.866357e+04 2.075239e+03
291 4.036418e+03 4.866357e+04
292 -1.347336e+04 4.036418e+03
293 2.278928e+03 -1.347336e+04
294 -1.235231e+03 2.278928e+03
295 6.523467e+04 -1.235231e+03
296 2.278928e+03 6.523467e+04
297 4.866357e+04 2.278928e+03
298 2.278928e+03 4.866357e+04
299 2.278928e+03 2.278928e+03
300 4.546553e+04 2.278928e+03
301 -3.678073e+03 4.546553e+04
302 2.278928e+03 -3.678073e+03
303 4.866357e+04 2.278928e+03
304 2.220313e+04 4.866357e+04
305 9.847648e+02 2.220313e+04
306 4.866357e+04 9.847648e+02
307 3.802061e+04 4.866357e+04
308 2.180264e+04 3.802061e+04
309 4.866357e+04 2.180264e+04
310 4.866357e+04 4.866357e+04
311 -4.898008e+04 4.866357e+04
312 4.866357e+04 -4.898008e+04
313 4.739741e+04 4.866357e+04
314 1.914256e+03 4.739741e+04
315 1.703036e+03 1.914256e+03
316 4.062286e+03 1.703036e+03
317 -9.101933e+03 4.062286e+03
318 1.424190e+03 -9.101933e+03
319 2.278928e+03 1.424190e+03
320 2.278928e+03 2.278928e+03
321 2.278928e+03 2.278928e+03
322 2.278928e+03 2.278928e+03
323 4.866357e+04 2.278928e+03
324 2.278928e+03 4.866357e+04
325 -7.358770e+02 2.278928e+03
326 4.866357e+04 -7.358770e+02
327 4.866357e+04 4.866357e+04
328 2.278928e+03 4.866357e+04
329 2.278928e+03 2.278928e+03
330 2.278928e+03 2.278928e+03
331 2.278928e+03 2.278928e+03
332 2.278928e+03 2.278928e+03
333 2.475569e+03 2.278928e+03
334 1.985279e+03 2.475569e+03
335 2.278928e+03 1.985279e+03
336 -2.146158e+02 2.278928e+03
337 2.278928e+03 -2.146158e+02
338 2.008275e+03 2.278928e+03
339 2.363771e+03 2.008275e+03
340 -1.437134e+04 2.363771e+03
341 2.278928e+03 -1.437134e+04
342 2.278928e+03 2.278928e+03
343 2.278928e+03 2.278928e+03
344 2.278928e+03 2.278928e+03
345 2.278928e+03 2.278928e+03
346 -5.911648e+03 2.278928e+03
347 2.278928e+03 -5.911648e+03
348 2.278928e+03 2.278928e+03
349 2.278928e+03 2.278928e+03
350 4.546473e+03 2.278928e+03
351 -2.025177e+02 4.546473e+03
352 1.650416e+04 -2.025177e+02
353 -3.517876e+03 1.650416e+04
354 2.278928e+03 -3.517876e+03
355 2.278928e+03 2.278928e+03
356 -3.860831e+04 2.278928e+03
357 2.278928e+03 -3.860831e+04
358 2.278928e+03 2.278928e+03
359 2.278928e+03 2.278928e+03
360 1.704788e+03 2.278928e+03
361 2.278928e+03 1.704788e+03
362 2.278928e+03 2.278928e+03
363 3.910634e+04 2.278928e+03
364 -5.914782e+00 3.910634e+04
365 2.278928e+03 -5.914782e+00
366 2.278928e+03 2.278928e+03
367 2.278928e+03 2.278928e+03
368 2.278928e+03 2.278928e+03
369 2.278928e+03 2.278928e+03
370 2.278928e+03 2.278928e+03
371 -3.681370e+03 2.278928e+03
372 -2.850358e+04 -3.681370e+03
373 -1.725417e+04 -2.850358e+04
374 2.278928e+03 -1.725417e+04
375 2.278928e+03 2.278928e+03
376 2.505398e+03 2.278928e+03
377 1.508541e+03 2.505398e+03
378 2.278928e+03 1.508541e+03
379 2.278928e+03 2.278928e+03
380 2.278928e+03 2.278928e+03
381 -1.871470e+04 2.278928e+03
382 2.192394e+04 -1.871470e+04
383 2.278928e+03 2.192394e+04
384 -1.803844e+03 2.278928e+03
385 -8.902396e+04 -1.803844e+03
386 7.428261e+03 -8.902396e+04
387 1.289292e+04 7.428261e+03
388 1.015009e+03 1.289292e+04
389 -7.653258e+03 1.015009e+03
390 7.706070e+03 -7.653258e+03
391 -7.972913e+03 7.706070e+03
392 -1.788966e+04 -7.972913e+03
393 4.135504e+04 -1.788966e+04
394 -1.467653e+05 4.135504e+04
395 -1.455957e+04 -1.467653e+05
396 -1.816617e+04 -1.455957e+04
397 8.697146e+04 -1.816617e+04
398 -9.947973e+03 8.697146e+04
399 -1.256977e+05 -9.947973e+03
400 2.623251e+03 -1.256977e+05
401 1.525061e+04 2.623251e+03
402 -2.562332e+03 1.525061e+04
403 -4.190231e+04 -2.562332e+03
404 -1.339545e+04 -4.190231e+04
405 -1.232814e+04 -1.339545e+04
406 1.477185e+03 -1.232814e+04
407 8.848026e+03 1.477185e+03
408 -1.247915e+04 8.848026e+03
409 -4.548403e+04 -1.247915e+04
410 -3.135735e+04 -4.548403e+04
411 -2.085961e+05 -3.135735e+04
412 -6.125829e+03 -2.085961e+05
413 -1.058975e+05 -6.125829e+03
414 -5.761806e+04 -1.058975e+05
415 -5.485962e+04 -5.761806e+04
416 1.507550e+05 -5.485962e+04
417 -1.990883e+04 1.507550e+05
418 -1.697357e+04 -1.990883e+04
419 -4.274209e+04 -1.697357e+04
420 -8.140141e+04 -4.274209e+04
421 -9.185457e+04 -8.140141e+04
422 -6.273686e+05 -9.185457e+04
423 -3.134823e+05 -6.273686e+05
424 -6.283221e+04 -3.134823e+05
425 -2.139888e+05 -6.283221e+04
426 -2.420959e+05 -2.139888e+05
427 -1.555632e+05 -2.420959e+05
428 -1.122225e+05 -1.555632e+05
429 1.590189e+05 -1.122225e+05
430 -3.740122e+05 1.590189e+05
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/76n4m1291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/86n4m1291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/9zem71291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/freestat/rcomp/tmp/10zem71291330748.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/11kekv1291330748.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/12nx111291330748.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/13jpgr1291330748.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/145pfx1291330748.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/15q8dl1291330748.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/164itc1291330748.tab")
+ }
>
> try(system("convert tmp/1avpv1291330748.ps tmp/1avpv1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/2kmog1291330748.ps tmp/2kmog1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/3kmog1291330748.ps tmp/3kmog1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/4kmog1291330748.ps tmp/4kmog1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/5vdnj1291330748.ps tmp/5vdnj1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/6vdnj1291330748.ps tmp/6vdnj1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/76n4m1291330748.ps tmp/76n4m1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/86n4m1291330748.ps tmp/86n4m1291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/9zem71291330748.ps tmp/9zem71291330748.png",intern=TRUE))
character(0)
> try(system("convert tmp/10zem71291330748.ps tmp/10zem71291330748.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
15.704 3.040 16.122