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(6282929
+ ,1
+ ,162556
+ ,1081
+ ,213118
+ ,4324047
+ ,1
+ ,29790
+ ,309
+ ,81767
+ ,4108272
+ ,1
+ ,87550
+ ,458
+ ,153198
+ ,-1212617
+ ,0
+ ,84738
+ ,588
+ ,-26007
+ ,1485329
+ ,1
+ ,54660
+ ,299
+ ,126942
+ ,1779876
+ ,1
+ ,42634
+ ,156
+ ,157214
+ ,1367203
+ ,0
+ ,40949
+ ,481
+ ,129352
+ ,2519076
+ ,1
+ ,42312
+ ,323
+ ,234817
+ ,912684
+ ,1
+ ,37704
+ ,452
+ ,60448
+ ,1443586
+ ,1
+ ,16275
+ ,109
+ ,47818
+ ,1220017
+ ,0
+ ,25830
+ ,115
+ ,245546
+ ,984885
+ ,0
+ ,12679
+ ,110
+ ,48020
+ ,1457425
+ ,1
+ ,18014
+ ,239
+ ,-1710
+ ,-572920
+ ,0
+ ,43556
+ ,247
+ ,32648
+ ,929144
+ ,1
+ ,24524
+ ,497
+ ,95350
+ ,1151176
+ ,0
+ ,6532
+ ,103
+ ,151352
+ ,790090
+ ,0
+ ,7123
+ ,109
+ ,288170
+ ,774497
+ ,1
+ ,20813
+ ,502
+ ,114337
+ ,990576
+ ,1
+ ,37597
+ ,248
+ ,37884
+ ,454195
+ ,0
+ ,17821
+ ,373
+ ,122844
+ ,876607
+ ,1
+ ,12988
+ ,119
+ ,82340
+ ,711969
+ ,1
+ ,22330
+ ,84
+ ,79801
+ ,702380
+ ,0
+ ,13326
+ ,102
+ ,165548
+ ,264449
+ ,0
+ ,16189
+ ,295
+ ,116384
+ ,450033
+ ,0
+ ,7146
+ ,105
+ ,134028
+ ,541063
+ ,0
+ ,15824
+ ,64
+ ,63838
+ ,588864
+ ,1
+ ,26088
+ ,267
+ ,74996
+ ,-37216
+ ,0
+ ,11326
+ ,129
+ ,31080
+ ,783310
+ ,0
+ ,8568
+ ,37
+ ,32168
+ ,467359
+ ,0
+ ,14416
+ ,361
+ ,49857
+ ,688779
+ ,1
+ ,3369
+ ,28
+ ,87161
+ ,608419
+ ,1
+ ,11819
+ ,85
+ ,106113
+ ,696348
+ ,1
+ ,6620
+ ,44
+ ,80570
+ ,597793
+ ,1
+ ,4519
+ ,49
+ ,102129
+ ,821730
+ ,0
+ ,2220
+ ,22
+ ,301670
+ ,377934
+ ,0
+ ,18562
+ ,155
+ ,102313
+ ,651939
+ ,0
+ ,10327
+ ,91
+ ,88577
+ ,697458
+ ,1
+ ,5336
+ ,81
+ ,112477
+ ,700368
+ ,1
+ ,2365
+ ,79
+ ,191778
+ ,225986
+ ,0
+ ,4069
+ ,145
+ ,79804
+ ,348695
+ ,0
+ ,7710
+ ,816
+ ,128294
+ ,373683
+ ,0
+ ,13718
+ ,61
+ ,96448
+ ,501709
+ ,0
+ ,4525
+ ,226
+ ,93811
+ ,413743
+ ,0
+ ,6869
+ ,105
+ ,117520
+ ,379825
+ ,0
+ ,4628
+ ,62
+ ,69159
+ ,336260
+ ,1
+ ,3653
+ ,24
+ ,101792
+ ,636765
+ ,1
+ ,1265
+ ,26
+ ,210568
+ ,481231
+ ,1
+ ,7489
+ ,322
+ ,136996
+ ,469107
+ ,0
+ ,4901
+ ,84
+ ,121920
+ ,211928
+ ,0
+ ,2284
+ ,33
+ ,76403
+ ,563925
+ ,1
+ ,3160
+ ,108
+ ,108094
+ ,511939
+ ,1
+ ,4150
+ ,150
+ ,134759
+ ,521016
+ ,1
+ ,7285
+ ,115
+ ,188873
+ ,543856
+ ,1
+ ,1134
+ ,162
+ ,146216
+ ,329304
+ ,1
+ ,4658
+ ,158
+ ,156608
+ ,423262
+ ,0
+ ,2384
+ ,97
+ ,61348
+ ,509665
+ ,0
+ ,3748
+ ,9
+ ,50350
+ ,455881
+ ,0
+ ,5371
+ ,66
+ ,87720
+ ,367772
+ ,0
+ ,1285
+ ,107
+ ,99489
+ ,406339
+ ,1
+ ,9327
+ ,101
+ ,87419
+ ,493408
+ ,1
+ ,5565
+ ,47
+ ,94355
+ ,232942
+ ,0
+ ,1528
+ ,38
+ ,60326
+ ,416002
+ ,1
+ ,3122
+ ,34
+ ,94670
+ ,337430
+ ,1
+ ,7317
+ ,84
+ ,82425
+ ,361517
+ ,0
+ ,2675
+ ,79
+ ,59017
+ ,360962
+ ,0
+ ,13253
+ ,947
+ ,90829
+ ,235561
+ ,0
+ ,880
+ ,74
+ ,80791
+ ,408247
+ ,1
+ ,2053
+ ,53
+ ,100423
+ ,450296
+ ,0
+ ,1424
+ ,94
+ ,131116
+ ,418799
+ ,1
+ ,4036
+ ,63
+ ,100269
+ ,247405
+ ,1
+ ,3045
+ ,58
+ ,27330
+ ,378519
+ ,0
+ ,5119
+ ,49
+ ,39039
+ ,326638
+ ,0
+ ,1431
+ ,34
+ ,106885
+ ,328233
+ ,0
+ ,554
+ ,11
+ ,79285
+ ,386225
+ ,0
+ ,1975
+ ,35
+ ,118881
+ ,283662
+ ,1
+ ,1286
+ ,17
+ ,77623
+ ,370225
+ ,0
+ ,1012
+ ,47
+ ,114768
+ ,269236
+ ,0
+ ,810
+ ,43
+ ,74015
+ ,365732
+ ,0
+ ,1280
+ ,117
+ ,69465
+ ,420383
+ ,1
+ ,666
+ ,171
+ ,117869
+ ,345811
+ ,0
+ ,1380
+ ,26
+ ,60982
+ ,431809
+ ,1
+ ,4608
+ ,73
+ ,90131
+ ,418876
+ ,0
+ ,876
+ ,59
+ ,138971
+ ,297476
+ ,0
+ ,814
+ ,18
+ ,39625
+ ,416776
+ ,0
+ ,514
+ ,15
+ ,102725
+ ,357257
+ ,1
+ ,5692
+ ,72
+ ,64239
+ ,458343
+ ,0
+ ,3642
+ ,86
+ ,90262
+ ,388386
+ ,0
+ ,540
+ ,14
+ ,103960
+ ,358934
+ ,0
+ ,2099
+ ,64
+ ,106611
+ ,407560
+ ,0
+ ,567
+ ,11
+ ,103345
+ ,392558
+ ,0
+ ,2001
+ ,52
+ ,95551
+ ,373177
+ ,1
+ ,2949
+ ,41
+ ,82903
+ ,428370
+ ,0
+ ,2253
+ ,99
+ ,63593
+ ,369419
+ ,1
+ ,6533
+ ,75
+ ,126910
+ ,358649
+ ,0
+ ,1889
+ ,45
+ ,37527
+ ,376641
+ ,1
+ ,3055
+ ,43
+ ,60247
+ ,467427
+ ,0
+ ,272
+ ,8
+ ,112995
+ ,364885
+ ,1
+ ,1414
+ ,198
+ ,70184
+ ,436230
+ ,0
+ ,2564
+ ,22
+ ,130140
+ ,329118
+ ,1
+ ,1383
+ ,11
+ ,73221
+ ,317365
+ ,1
+ ,1261
+ ,33
+ ,76114
+ ,286849
+ ,0
+ ,975
+ ,23
+ ,90534
+ ,376685
+ ,0
+ ,3366
+ ,80
+ ,108479
+ ,407198
+ ,0
+ ,576
+ ,18
+ ,113761
+ ,377772
+ ,0
+ ,1306
+ ,28
+ ,68696
+ ,271483
+ ,0
+ ,746
+ ,23
+ ,71561
+ ,153661
+ ,1
+ ,3192
+ ,60
+ ,59831
+ ,513294
+ ,1
+ ,2045
+ ,20
+ ,97890
+ ,324881
+ ,0
+ ,5477
+ ,59
+ ,101481
+ ,264512
+ ,1
+ ,1932
+ ,36
+ ,72954
+ ,420968
+ ,0
+ ,936
+ ,30
+ ,67939
+ ,129302
+ ,1
+ ,3437
+ ,47
+ ,48022
+ ,191521
+ ,0
+ ,5131
+ ,71
+ ,86111
+ ,268673
+ ,1
+ ,2397
+ ,14
+ ,74020
+ ,353179
+ ,1
+ ,1389
+ ,9
+ ,57530
+ ,354624
+ ,0
+ ,1503
+ ,39
+ ,56364
+ ,363713
+ ,0
+ ,402
+ ,26
+ ,84990
+ ,456657
+ ,0
+ ,2239
+ ,21
+ ,88590
+ ,211742
+ ,1
+ ,2234
+ ,16
+ ,77200
+ ,338381
+ ,0
+ ,837
+ ,69
+ ,61262
+ ,418530
+ ,0
+ ,10579
+ ,92
+ ,110309
+ ,351483
+ ,0
+ ,875
+ ,14
+ ,67000
+ ,372928
+ ,0
+ ,1395
+ ,103
+ ,93099
+ ,485538
+ ,1
+ ,1659
+ ,29
+ ,107577
+ ,279268
+ ,1
+ ,2647
+ ,37
+ ,62920
+ ,219060
+ ,1
+ ,3294
+ ,23
+ ,75832
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325314
+ ,0
+ ,94
+ ,7
+ ,60793
+ ,322046
+ ,0
+ ,422
+ ,28
+ ,57935
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325599
+ ,0
+ ,34
+ ,8
+ ,60630
+ ,377028
+ ,0
+ ,1558
+ ,63
+ ,55637
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,323850
+ ,0
+ ,43
+ ,3
+ ,60887
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331514
+ ,0
+ ,316
+ ,9
+ ,60505
+ ,325632
+ ,0
+ ,115
+ ,13
+ ,60945
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,322265
+ ,0
+ ,389
+ ,14
+ ,58990
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325906
+ ,0
+ ,1002
+ ,15
+ ,56750
+ ,325985
+ ,0
+ ,36
+ ,3
+ ,60894
+ ,346145
+ ,0
+ ,460
+ ,15
+ ,63346
+ ,325898
+ ,0
+ ,309
+ ,11
+ ,56535
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325356
+ ,0
+ ,9
+ ,6
+ ,60835
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325930
+ ,0
+ ,14
+ ,1
+ ,61016
+ ,318020
+ ,0
+ ,520
+ ,10
+ ,58650
+ ,326389
+ ,0
+ ,1766
+ ,73
+ ,60438
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,302925
+ ,0
+ ,458
+ ,11
+ ,58625
+ ,325540
+ ,0
+ ,20
+ ,3
+ ,60938
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,326736
+ ,0
+ ,98
+ ,2
+ ,61490
+ ,340580
+ ,0
+ ,405
+ ,7
+ ,60845
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331828
+ ,0
+ ,483
+ ,27
+ ,60830
+ ,323299
+ ,0
+ ,454
+ ,51
+ ,63261
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,387722
+ ,0
+ ,757
+ ,19
+ ,45689
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324598
+ ,0
+ ,36
+ ,4
+ ,61564
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,328726
+ ,0
+ ,203
+ ,9
+ ,61938
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325043
+ ,0
+ ,90
+ ,8
+ ,60951
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325806
+ ,1
+ ,71
+ ,1
+ ,60745
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,387732
+ ,0
+ ,972
+ ,34
+ ,71642
+ ,349729
+ ,1
+ ,531
+ ,10
+ ,71641
+ ,332202
+ ,0
+ ,604
+ ,38
+ ,55792
+ ,305442
+ ,1
+ ,283
+ ,10
+ ,71873
+ ,329537
+ ,1
+ ,23
+ ,5
+ ,62555
+ ,327055
+ ,1
+ ,638
+ ,14
+ ,60370
+ ,356245
+ ,1
+ ,699
+ ,16
+ ,64873
+ ,328451
+ ,0
+ ,149
+ ,5
+ ,62041
+ ,307062
+ ,0
+ ,226
+ ,5
+ ,65745
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331345
+ ,0
+ ,275
+ ,4
+ ,59500
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331824
+ ,0
+ ,141
+ ,6
+ ,61630
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325685
+ ,0
+ ,28
+ ,2
+ ,60890
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,404480
+ ,1
+ ,2566
+ ,80
+ ,113521
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,318314
+ ,1
+ ,472
+ ,20
+ ,80045
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,311807
+ ,1
+ ,203
+ ,27
+ ,50804
+ ,337724
+ ,1
+ ,496
+ ,17
+ ,87390
+ ,326431
+ ,1
+ ,10
+ ,2
+ ,61656
+ ,327556
+ ,1
+ ,63
+ ,4
+ ,65688
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,356850
+ ,1
+ ,1136
+ ,32
+ ,48522
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,322741
+ ,0
+ ,267
+ ,32
+ ,57640
+ ,310902
+ ,0
+ ,474
+ ,20
+ ,61977
+ ,324295
+ ,0
+ ,534
+ ,7
+ ,62620
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,326156
+ ,0
+ ,15
+ ,8
+ ,60831
+ ,326960
+ ,0
+ ,397
+ ,28
+ ,60646
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,333411
+ ,0
+ ,1061
+ ,20
+ ,56225
+ ,297761
+ ,0
+ ,288
+ ,4
+ ,60510
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325536
+ ,0
+ ,3
+ ,2
+ ,60698
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325762
+ ,0
+ ,20
+ ,2
+ ,60805
+ ,327957
+ ,0
+ ,278
+ ,26
+ ,61404
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,318521
+ ,0
+ ,192
+ ,4
+ ,65276
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,319775
+ ,0
+ ,317
+ ,9
+ ,63915
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,332128
+ ,0
+ ,368
+ ,17
+ ,61686
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325486
+ ,0
+ ,2
+ ,1
+ ,60743
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325838
+ ,0
+ ,53
+ ,6
+ ,60349
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331767
+ ,0
+ ,94
+ ,3
+ ,61360
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324523
+ ,0
+ ,24
+ ,8
+ ,59818
+ ,339995
+ ,0
+ ,2332
+ ,4
+ ,72680
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,319582
+ ,0
+ ,131
+ ,11
+ ,61808
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,307245
+ ,0
+ ,206
+ ,9
+ ,53110
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,317967
+ ,0
+ ,167
+ ,2
+ ,64245
+ ,331488
+ ,0
+ ,622
+ ,73
+ ,73007
+ ,335452
+ ,0
+ ,885
+ ,85
+ ,82732
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,334184
+ ,0
+ ,365
+ ,8
+ ,54820
+ ,313213
+ ,0
+ ,364
+ ,35
+ ,47705
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,348678
+ ,0
+ ,226
+ ,12
+ ,72835
+ ,328727
+ ,0
+ ,307
+ ,15
+ ,58856
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,387978
+ ,0
+ ,188
+ ,11
+ ,77655
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,336704
+ ,0
+ ,138
+ ,6
+ ,69817
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,322076
+ ,0
+ ,125
+ ,12
+ ,60798
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,334272
+ ,0
+ ,282
+ ,30
+ ,62452
+ ,338197
+ ,0
+ ,335
+ ,33
+ ,64175
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,321024
+ ,1
+ ,813
+ ,82
+ ,67440
+ ,322145
+ ,0
+ ,176
+ ,28
+ ,68136
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,323351
+ ,0
+ ,249
+ ,72
+ ,56726
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,327748
+ ,0
+ ,333
+ ,13
+ ,70811
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,328157
+ ,0
+ ,30
+ ,4
+ ,62045
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,311594
+ ,0
+ ,249
+ ,62
+ ,54323
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,335962
+ ,0
+ ,165
+ ,24
+ ,62841
+ ,372426
+ ,0
+ ,453
+ ,21
+ ,81125
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,319844
+ ,0
+ ,53
+ ,14
+ ,59506
+ ,355822
+ ,1
+ ,382
+ ,21
+ ,59365
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324047
+ ,1
+ ,30
+ ,4
+ ,60798
+ ,311464
+ ,0
+ ,290
+ ,2
+ ,58790
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,353417
+ ,0
+ ,366
+ ,53
+ ,61808
+ ,325590
+ ,0
+ ,2
+ ,9
+ ,60735
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,328576
+ ,1
+ ,209
+ ,13
+ ,64016
+ ,326126
+ ,0
+ ,384
+ ,22
+ ,54683
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,369376
+ ,0
+ ,365
+ ,83
+ ,87192
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,332013
+ ,1
+ ,49
+ ,8
+ ,64107
+ ,325871
+ ,0
+ ,3
+ ,4
+ ,60761
+ ,342165
+ ,0
+ ,133
+ ,14
+ ,65990
+ ,324967
+ ,0
+ ,32
+ ,1
+ ,59988
+ ,314832
+ ,0
+ ,368
+ ,17
+ ,61167
+ ,325557
+ ,0
+ ,1
+ ,6
+ ,60719
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,322649
+ ,0
+ ,22
+ ,2
+ ,60722
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324598
+ ,0
+ ,96
+ ,5
+ ,60379
+ ,325567
+ ,0
+ ,1
+ ,2
+ ,60727
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324005
+ ,0
+ ,81
+ ,7
+ ,60925
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325748
+ ,0
+ ,26
+ ,1
+ ,60896
+ ,323385
+ ,0
+ ,125
+ ,13
+ ,59734
+ ,315409
+ ,0
+ ,304
+ ,15
+ ,62969
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,312275
+ ,0
+ ,119
+ ,6
+ ,59118
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,320576
+ ,0
+ ,312
+ ,14
+ ,58598
+ ,325246
+ ,0
+ ,60
+ ,10
+ ,61124
+ ,332961
+ ,0
+ ,587
+ ,12
+ ,59595
+ ,323010
+ ,0
+ ,135
+ ,2
+ ,62065
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,345253
+ ,0
+ ,514
+ ,52
+ ,78780
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325559
+ ,0
+ ,1
+ ,4
+ ,60722
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,319634
+ ,1
+ ,58
+ ,3
+ ,61600
+ ,319951
+ ,0
+ ,180
+ ,11
+ ,59635
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,318519
+ ,0
+ ,448
+ ,40
+ ,59781
+ ,343222
+ ,0
+ ,227
+ ,9
+ ,76644
+ ,317234
+ ,0
+ ,174
+ ,1
+ ,64820
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,314025
+ ,0
+ ,121
+ ,24
+ ,56178
+ ,320249
+ ,0
+ ,607
+ ,11
+ ,60436
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,349365
+ ,0
+ ,530
+ ,60
+ ,73433
+ ,289197
+ ,0
+ ,571
+ ,80
+ ,41477
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,329245
+ ,0
+ ,78
+ ,16
+ ,62700
+ ,240869
+ ,0
+ ,2489
+ ,40
+ ,67804
+ ,327182
+ ,0
+ ,131
+ ,6
+ ,59661
+ ,322876
+ ,0
+ ,923
+ ,8
+ ,58620
+ ,323117
+ ,0
+ ,72
+ ,3
+ ,60398
+ ,306351
+ ,0
+ ,572
+ ,16
+ ,58580
+ ,335137
+ ,0
+ ,397
+ ,10
+ ,62710
+ ,308271
+ ,0
+ ,450
+ ,8
+ ,59325
+ ,301731
+ ,0
+ ,622
+ ,7
+ ,60950
+ ,382409
+ ,0
+ ,694
+ ,8
+ ,68060
+ ,279230
+ ,1
+ ,3425
+ ,12
+ ,83620
+ ,298731
+ ,0
+ ,562
+ ,13
+ ,58456
+ ,243650
+ ,0
+ ,4917
+ ,42
+ ,52811
+ ,532682
+ ,1
+ ,1442
+ ,118
+ ,121173
+ ,319771
+ ,0
+ ,529
+ ,9
+ ,63870
+ ,171493
+ ,1
+ ,2126
+ ,138
+ ,21001
+ ,347262
+ ,0
+ ,1061
+ ,5
+ ,70415
+ ,343945
+ ,0
+ ,776
+ ,9
+ ,64230
+ ,311874
+ ,0
+ ,611
+ ,8
+ ,59190
+ ,302211
+ ,1
+ ,1526
+ ,25
+ ,69351
+ ,316708
+ ,0
+ ,592
+ ,7
+ ,64270
+ ,333463
+ ,0
+ ,1182
+ ,13
+ ,70694
+ ,344282
+ ,0
+ ,621
+ ,16
+ ,68005
+ ,319635
+ ,0
+ ,989
+ ,11
+ ,58930
+ ,301186
+ ,0
+ ,438
+ ,11
+ ,58320
+ ,300381
+ ,0
+ ,726
+ ,3
+ ,69980
+ ,318765
+ ,0
+ ,1303
+ ,61
+ ,69863
+ ,286146
+ ,1
+ ,6341
+ ,24
+ ,63255
+ ,306844
+ ,1
+ ,1164
+ ,17
+ ,57320
+ ,307705
+ ,1
+ ,3310
+ ,33
+ ,75230
+ ,312448
+ ,0
+ ,1366
+ ,7
+ ,79420
+ ,299715
+ ,0
+ ,965
+ ,3
+ ,73490
+ ,373399
+ ,0
+ ,3256
+ ,66
+ ,35250
+ ,299446
+ ,1
+ ,1135
+ ,17
+ ,62285
+ ,325586
+ ,0
+ ,1270
+ ,26
+ ,69206
+ ,291221
+ ,0
+ ,661
+ ,3
+ ,65920
+ ,261173
+ ,0
+ ,1013
+ ,2
+ ,69770
+ ,255027
+ ,0
+ ,2844
+ ,67
+ ,72683
+ ,-78375
+ ,1
+ ,11528
+ ,70
+ ,-14545
+ ,-58143
+ ,0
+ ,6526
+ ,26
+ ,55830
+ ,227033
+ ,0
+ ,2264
+ ,24
+ ,55174
+ ,235098
+ ,1
+ ,4461
+ ,94
+ ,67038
+ ,21267
+ ,0
+ ,3999
+ ,30
+ ,51252
+ ,238675
+ ,0
+ ,35624
+ ,223
+ ,157278
+ ,197687
+ ,0
+ ,9252
+ ,48
+ ,79510
+ ,418341
+ ,0
+ ,15236
+ ,90
+ ,77440
+ ,-297706
+ ,0
+ ,18073
+ ,180
+ ,27284)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('Wealth'
+ ,'Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Dividends')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('Wealth','Group','Costs','Trades','Dividends'),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 = '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Wealth Group Costs Trades Dividends
1 6282929 1 162556 1081 213118
2 4324047 1 29790 309 81767
3 4108272 1 87550 458 153198
4 -1212617 0 84738 588 -26007
5 1485329 1 54660 299 126942
6 1779876 1 42634 156 157214
7 1367203 0 40949 481 129352
8 2519076 1 42312 323 234817
9 912684 1 37704 452 60448
10 1443586 1 16275 109 47818
11 1220017 0 25830 115 245546
12 984885 0 12679 110 48020
13 1457425 1 18014 239 -1710
14 -572920 0 43556 247 32648
15 929144 1 24524 497 95350
16 1151176 0 6532 103 151352
17 790090 0 7123 109 288170
18 774497 1 20813 502 114337
19 990576 1 37597 248 37884
20 454195 0 17821 373 122844
21 876607 1 12988 119 82340
22 711969 1 22330 84 79801
23 702380 0 13326 102 165548
24 264449 0 16189 295 116384
25 450033 0 7146 105 134028
26 541063 0 15824 64 63838
27 588864 1 26088 267 74996
28 -37216 0 11326 129 31080
29 783310 0 8568 37 32168
30 467359 0 14416 361 49857
31 688779 1 3369 28 87161
32 608419 1 11819 85 106113
33 696348 1 6620 44 80570
34 597793 1 4519 49 102129
35 821730 0 2220 22 301670
36 377934 0 18562 155 102313
37 651939 0 10327 91 88577
38 697458 1 5336 81 112477
39 700368 1 2365 79 191778
40 225986 0 4069 145 79804
41 348695 0 7710 816 128294
42 373683 0 13718 61 96448
43 501709 0 4525 226 93811
44 413743 0 6869 105 117520
45 379825 0 4628 62 69159
46 336260 1 3653 24 101792
47 636765 1 1265 26 210568
48 481231 1 7489 322 136996
49 469107 0 4901 84 121920
50 211928 0 2284 33 76403
51 563925 1 3160 108 108094
52 511939 1 4150 150 134759
53 521016 1 7285 115 188873
54 543856 1 1134 162 146216
55 329304 1 4658 158 156608
56 423262 0 2384 97 61348
57 509665 0 3748 9 50350
58 455881 0 5371 66 87720
59 367772 0 1285 107 99489
60 406339 1 9327 101 87419
61 493408 1 5565 47 94355
62 232942 0 1528 38 60326
63 416002 1 3122 34 94670
64 337430 1 7317 84 82425
65 361517 0 2675 79 59017
66 360962 0 13253 947 90829
67 235561 0 880 74 80791
68 408247 1 2053 53 100423
69 450296 0 1424 94 131116
70 418799 1 4036 63 100269
71 247405 1 3045 58 27330
72 378519 0 5119 49 39039
73 326638 0 1431 34 106885
74 328233 0 554 11 79285
75 386225 0 1975 35 118881
76 283662 1 1286 17 77623
77 370225 0 1012 47 114768
78 269236 0 810 43 74015
79 365732 0 1280 117 69465
80 420383 1 666 171 117869
81 345811 0 1380 26 60982
82 431809 1 4608 73 90131
83 418876 0 876 59 138971
84 297476 0 814 18 39625
85 416776 0 514 15 102725
86 357257 1 5692 72 64239
87 458343 0 3642 86 90262
88 388386 0 540 14 103960
89 358934 0 2099 64 106611
90 407560 0 567 11 103345
91 392558 0 2001 52 95551
92 373177 1 2949 41 82903
93 428370 0 2253 99 63593
94 369419 1 6533 75 126910
95 358649 0 1889 45 37527
96 376641 1 3055 43 60247
97 467427 0 272 8 112995
98 364885 1 1414 198 70184
99 436230 0 2564 22 130140
100 329118 1 1383 11 73221
101 317365 1 1261 33 76114
102 286849 0 975 23 90534
103 376685 0 3366 80 108479
104 407198 0 576 18 113761
105 377772 0 1306 28 68696
106 271483 0 746 23 71561
107 153661 1 3192 60 59831
108 513294 1 2045 20 97890
109 324881 0 5477 59 101481
110 264512 1 1932 36 72954
111 420968 0 936 30 67939
112 129302 1 3437 47 48022
113 191521 0 5131 71 86111
114 268673 1 2397 14 74020
115 353179 1 1389 9 57530
116 354624 0 1503 39 56364
117 363713 0 402 26 84990
118 456657 0 2239 21 88590
119 211742 1 2234 16 77200
120 338381 0 837 69 61262
121 418530 0 10579 92 110309
122 351483 0 875 14 67000
123 372928 0 1395 103 93099
124 485538 1 1659 29 107577
125 279268 1 2647 37 62920
126 219060 1 3294 23 75832
127 325560 0 0 0 60720
128 325314 0 94 7 60793
129 322046 0 422 28 57935
130 325560 0 0 0 60720
131 325599 0 34 8 60630
132 377028 0 1558 63 55637
133 325560 0 0 0 60720
134 323850 0 43 3 60887
135 325560 0 0 0 60720
136 331514 0 316 9 60505
137 325632 0 115 13 60945
138 325560 0 0 0 60720
139 325560 0 0 0 60720
140 325560 0 0 0 60720
141 322265 0 389 14 58990
142 325560 0 0 0 60720
143 325906 0 1002 15 56750
144 325985 0 36 3 60894
145 346145 0 460 15 63346
146 325898 0 309 11 56535
147 325560 0 0 0 60720
148 325356 0 9 6 60835
149 325560 0 0 0 60720
150 325930 0 14 1 61016
151 318020 0 520 10 58650
152 326389 0 1766 73 60438
153 325560 0 0 0 60720
154 302925 0 458 11 58625
155 325540 0 20 3 60938
156 325560 0 0 0 60720
157 325560 0 0 0 60720
158 326736 0 98 2 61490
159 340580 0 405 7 60845
160 325560 0 0 0 60720
161 325560 0 0 0 60720
162 325560 0 0 0 60720
163 325560 0 0 0 60720
164 331828 0 483 27 60830
165 323299 0 454 51 63261
166 325560 0 0 0 60720
167 325560 0 0 0 60720
168 387722 0 757 19 45689
169 325560 0 0 0 60720
170 325560 0 0 0 60720
171 325560 0 0 0 60720
172 324598 0 36 4 61564
173 325560 0 0 0 60720
174 328726 0 203 9 61938
175 325560 0 0 0 60720
176 325043 0 90 8 60951
177 325560 0 0 0 60720
178 325806 1 71 1 60745
179 325560 0 0 0 60720
180 325560 0 0 0 60720
181 387732 0 972 34 71642
182 349729 1 531 10 71641
183 332202 0 604 38 55792
184 305442 1 283 10 71873
185 329537 1 23 5 62555
186 327055 1 638 14 60370
187 356245 1 699 16 64873
188 328451 0 149 5 62041
189 307062 0 226 5 65745
190 325560 0 0 0 60720
191 331345 0 275 4 59500
192 325560 0 0 0 60720
193 331824 0 141 6 61630
194 325560 0 0 0 60720
195 325685 0 28 2 60890
196 325560 1 0 0 60720
197 404480 1 2566 80 113521
198 325560 1 0 0 60720
199 325560 1 0 0 60720
200 318314 1 472 20 80045
201 325560 1 0 0 60720
202 325560 1 0 0 60720
203 325560 1 0 0 60720
204 311807 1 203 27 50804
205 337724 1 496 17 87390
206 326431 1 10 2 61656
207 327556 1 63 4 65688
208 325560 1 0 0 60720
209 356850 1 1136 32 48522
210 325560 1 0 0 60720
211 325560 0 0 0 60720
212 325560 0 0 0 60720
213 322741 0 267 32 57640
214 310902 0 474 20 61977
215 324295 0 534 7 62620
216 325560 0 0 0 60720
217 326156 0 15 8 60831
218 326960 0 397 28 60646
219 325560 0 0 0 60720
220 333411 0 1061 20 56225
221 297761 0 288 4 60510
222 325560 0 0 0 60720
223 325536 0 3 2 60698
224 325560 0 0 0 60720
225 325762 0 20 2 60805
226 327957 0 278 26 61404
227 325560 0 0 0 60720
228 325560 0 0 0 60720
229 318521 0 192 4 65276
230 325560 0 0 0 60720
231 319775 0 317 9 63915
232 325560 0 0 0 60720
233 325560 0 0 0 60720
234 332128 0 368 17 61686
235 325560 0 0 0 60720
236 325486 0 2 1 60743
237 325560 0 0 0 60720
238 325838 0 53 6 60349
239 325560 0 0 0 60720
240 325560 0 0 0 60720
241 325560 0 0 0 60720
242 331767 0 94 3 61360
243 325560 0 0 0 60720
244 324523 0 24 8 59818
245 339995 0 2332 4 72680
246 325560 1 0 0 60720
247 325560 1 0 0 60720
248 319582 0 131 11 61808
249 325560 1 0 0 60720
250 325560 1 0 0 60720
251 307245 0 206 9 53110
252 325560 1 0 0 60720
253 317967 0 167 2 64245
254 331488 0 622 73 73007
255 335452 0 885 85 82732
256 325560 0 0 0 60720
257 334184 0 365 8 54820
258 313213 0 364 35 47705
259 325560 1 0 0 60720
260 325560 0 0 0 60720
261 325560 1 0 0 60720
262 325560 1 0 0 60720
263 348678 0 226 12 72835
264 328727 0 307 15 58856
265 325560 1 0 0 60720
266 325560 0 0 0 60720
267 325560 1 0 0 60720
268 387978 0 188 11 77655
269 325560 1 0 0 60720
270 336704 0 138 6 69817
271 325560 1 0 0 60720
272 325560 1 0 0 60720
273 325560 1 0 0 60720
274 322076 0 125 12 60798
275 325560 0 0 0 60720
276 334272 0 282 30 62452
277 338197 0 335 33 64175
278 325560 0 0 0 60720
279 321024 1 813 82 67440
280 322145 0 176 28 68136
281 325560 1 0 0 60720
282 325560 0 0 0 60720
283 323351 0 249 72 56726
284 325560 0 0 0 60720
285 327748 0 333 13 70811
286 325560 0 0 0 60720
287 325560 1 0 0 60720
288 328157 0 30 4 62045
289 325560 1 0 0 60720
290 311594 0 249 62 54323
291 325560 1 0 0 60720
292 335962 0 165 24 62841
293 372426 0 453 21 81125
294 325560 0 0 0 60720
295 319844 0 53 14 59506
296 355822 1 382 21 59365
297 325560 0 0 0 60720
298 325560 1 0 0 60720
299 325560 0 0 0 60720
300 325560 0 0 0 60720
301 324047 1 30 4 60798
302 311464 0 290 2 58790
303 325560 0 0 0 60720
304 325560 1 0 0 60720
305 353417 0 366 53 61808
306 325590 0 2 9 60735
307 325560 1 0 0 60720
308 328576 1 209 13 64016
309 326126 0 384 22 54683
310 325560 1 0 0 60720
311 325560 1 0 0 60720
312 369376 0 365 83 87192
313 325560 1 0 0 60720
314 332013 1 49 8 64107
315 325871 0 3 4 60761
316 342165 0 133 14 65990
317 324967 0 32 1 59988
318 314832 0 368 17 61167
319 325557 0 1 6 60719
320 325560 0 0 0 60720
321 325560 0 0 0 60720
322 325560 0 0 0 60720
323 325560 0 0 0 60720
324 325560 1 0 0 60720
325 325560 0 0 0 60720
326 322649 0 22 2 60722
327 325560 1 0 0 60720
328 325560 1 0 0 60720
329 325560 0 0 0 60720
330 325560 0 0 0 60720
331 325560 0 0 0 60720
332 325560 0 0 0 60720
333 325560 0 0 0 60720
334 324598 0 96 5 60379
335 325567 0 1 2 60727
336 325560 0 0 0 60720
337 324005 0 81 7 60925
338 325560 0 0 0 60720
339 325748 0 26 1 60896
340 323385 0 125 13 59734
341 315409 0 304 15 62969
342 325560 0 0 0 60720
343 325560 0 0 0 60720
344 325560 0 0 0 60720
345 325560 0 0 0 60720
346 325560 0 0 0 60720
347 312275 0 119 6 59118
348 325560 0 0 0 60720
349 325560 0 0 0 60720
350 325560 0 0 0 60720
351 320576 0 312 14 58598
352 325246 0 60 10 61124
353 332961 0 587 12 59595
354 323010 0 135 2 62065
355 325560 0 0 0 60720
356 325560 0 0 0 60720
357 345253 0 514 52 78780
358 325560 0 0 0 60720
359 325560 0 0 0 60720
360 325560 0 0 0 60720
361 325559 0 1 4 60722
362 325560 0 0 0 60720
363 325560 0 0 0 60720
364 319634 1 58 3 61600
365 319951 0 180 11 59635
366 325560 0 0 0 60720
367 325560 0 0 0 60720
368 325560 0 0 0 60720
369 325560 0 0 0 60720
370 325560 0 0 0 60720
371 325560 0 0 0 60720
372 318519 0 448 40 59781
373 343222 0 227 9 76644
374 317234 0 174 1 64820
375 325560 0 0 0 60720
376 325560 0 0 0 60720
377 314025 0 121 24 56178
378 320249 0 607 11 60436
379 325560 0 0 0 60720
380 325560 0 0 0 60720
381 325560 0 0 0 60720
382 349365 0 530 60 73433
383 289197 0 571 80 41477
384 325560 0 0 0 60720
385 329245 0 78 16 62700
386 240869 0 2489 40 67804
387 327182 0 131 6 59661
388 322876 0 923 8 58620
389 323117 0 72 3 60398
390 306351 0 572 16 58580
391 335137 0 397 10 62710
392 308271 0 450 8 59325
393 301731 0 622 7 60950
394 382409 0 694 8 68060
395 279230 1 3425 12 83620
396 298731 0 562 13 58456
397 243650 0 4917 42 52811
398 532682 1 1442 118 121173
399 319771 0 529 9 63870
400 171493 1 2126 138 21001
401 347262 0 1061 5 70415
402 343945 0 776 9 64230
403 311874 0 611 8 59190
404 302211 1 1526 25 69351
405 316708 0 592 7 64270
406 333463 0 1182 13 70694
407 344282 0 621 16 68005
408 319635 0 989 11 58930
409 301186 0 438 11 58320
410 300381 0 726 3 69980
411 318765 0 1303 61 69863
412 286146 1 6341 24 63255
413 306844 1 1164 17 57320
414 307705 1 3310 33 75230
415 312448 0 1366 7 79420
416 299715 0 965 3 73490
417 373399 0 3256 66 35250
418 299446 1 1135 17 62285
419 325586 0 1270 26 69206
420 291221 0 661 3 65920
421 261173 0 1013 2 69770
422 255027 0 2844 67 72683
423 -78375 1 11528 70 -14545
424 -58143 0 6526 26 55830
425 227033 0 2264 24 55174
426 235098 1 4461 94 67038
427 21267 0 3999 30 51252
428 238675 0 35624 223 157278
429 197687 0 9252 48 79510
430 418341 0 15236 90 77440
431 -297706 0 18073 180 27284
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Dividends
56500.784 70846.109 26.482 -470.732 3.497
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-3145447 -20447 47198 56744 3267340
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 56500.7845 36533.7773 1.547 0.1227
Group 70846.1094 32881.3355 2.155 0.0318 *
Costs 26.4823 1.9134 13.840 < 2e-16 ***
Trades -470.7317 206.1013 -2.284 0.0229 *
Dividends 3.4966 0.4833 7.235 2.18e-12 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 299400 on 426 degrees of freedom
Multiple R-squared: 0.5565, Adjusted R-squared: 0.5523
F-statistic: 133.6 on 4 and 426 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 2.471451e-46 1.235726e-46
[2,] 1.0000000 4.814834e-79 2.407417e-79
[3,] 1.0000000 1.705622e-87 8.528108e-88
[4,] 1.0000000 2.115741e-89 1.057870e-89
[5,] 1.0000000 3.297289e-107 1.648645e-107
[6,] 1.0000000 3.034309e-124 1.517155e-124
[7,] 1.0000000 2.132831e-139 1.066416e-139
[8,] 1.0000000 5.694279e-151 2.847139e-151
[9,] 1.0000000 5.329587e-161 2.664794e-161
[10,] 1.0000000 2.283699e-162 1.141850e-162
[11,] 1.0000000 1.567413e-167 7.837067e-168
[12,] 1.0000000 2.995015e-177 1.497507e-177
[13,] 1.0000000 2.275749e-176 1.137875e-176
[14,] 1.0000000 1.479932e-183 7.399660e-184
[15,] 1.0000000 1.182463e-188 5.912316e-189
[16,] 1.0000000 9.793603e-189 4.896802e-189
[17,] 1.0000000 3.129126e-188 1.564563e-188
[18,] 1.0000000 2.304281e-187 1.152140e-187
[19,] 1.0000000 3.182174e-190 1.591087e-190
[20,] 1.0000000 2.314927e-195 1.157463e-195
[21,] 1.0000000 3.668875e-197 1.834437e-197
[22,] 1.0000000 2.909101e-214 1.454551e-214
[23,] 1.0000000 1.964745e-217 9.823726e-218
[24,] 1.0000000 1.014881e-221 5.074407e-222
[25,] 1.0000000 1.613932e-225 8.069660e-226
[26,] 1.0000000 1.852849e-232 9.264246e-233
[27,] 1.0000000 1.357399e-234 6.786993e-235
[28,] 1.0000000 1.918873e-234 9.594364e-235
[29,] 1.0000000 5.321350e-234 2.660675e-234
[30,] 1.0000000 4.075108e-242 2.037554e-242
[31,] 1.0000000 6.493496e-248 3.246748e-248
[32,] 1.0000000 4.203875e-249 2.101938e-249
[33,] 1.0000000 5.151715e-249 2.575858e-249
[34,] 1.0000000 6.188678e-249 3.094339e-249
[35,] 1.0000000 8.847765e-249 4.423883e-249
[36,] 1.0000000 2.189682e-249 1.094841e-249
[37,] 1.0000000 2.642796e-248 1.321398e-248
[38,] 1.0000000 3.267642e-248 1.633821e-248
[39,] 1.0000000 4.951756e-248 2.475878e-248
[40,] 1.0000000 1.665862e-248 8.329310e-249
[41,] 1.0000000 1.843172e-248 9.215858e-249
[42,] 1.0000000 1.218084e-247 6.090419e-248
[43,] 1.0000000 7.435821e-248 3.717911e-248
[44,] 1.0000000 6.502774e-249 3.251387e-249
[45,] 1.0000000 1.526112e-248 7.630560e-249
[46,] 1.0000000 6.527344e-249 3.263672e-249
[47,] 1.0000000 4.357798e-248 2.178899e-248
[48,] 1.0000000 1.944273e-250 9.721364e-251
[49,] 1.0000000 4.331229e-251 2.165614e-251
[50,] 1.0000000 1.638086e-256 8.190431e-257
[51,] 1.0000000 4.754339e-257 2.377169e-257
[52,] 1.0000000 3.496478e-256 1.748239e-256
[53,] 1.0000000 8.880228e-257 4.440114e-257
[54,] 1.0000000 4.096102e-258 2.048051e-258
[55,] 1.0000000 1.194414e-257 5.972071e-258
[56,] 1.0000000 7.649750e-257 3.824875e-257
[57,] 1.0000000 3.592465e-256 1.796233e-256
[58,] 1.0000000 8.856101e-256 4.428050e-256
[59,] 1.0000000 8.307795e-255 4.153898e-255
[60,] 1.0000000 2.358103e-255 1.179052e-255
[61,] 1.0000000 2.727739e-254 1.363870e-254
[62,] 1.0000000 3.077846e-253 1.538923e-253
[63,] 1.0000000 1.823895e-252 9.119477e-253
[64,] 1.0000000 1.887804e-251 9.439018e-252
[65,] 1.0000000 1.274720e-253 6.373598e-254
[66,] 1.0000000 2.184527e-253 1.092263e-253
[67,] 1.0000000 1.706563e-252 8.532816e-253
[68,] 1.0000000 1.049128e-251 5.245640e-252
[69,] 1.0000000 4.766407e-251 2.383203e-251
[70,] 1.0000000 1.301089e-250 6.505447e-251
[71,] 1.0000000 2.990510e-250 1.495255e-250
[72,] 1.0000000 2.050286e-249 1.025143e-249
[73,] 1.0000000 1.271508e-248 6.357538e-249
[74,] 1.0000000 7.826428e-248 3.913214e-248
[75,] 1.0000000 1.171875e-247 5.859377e-248
[76,] 1.0000000 9.946689e-248 4.973345e-248
[77,] 1.0000000 4.855968e-247 2.427984e-247
[78,] 1.0000000 5.363695e-246 2.681847e-246
[79,] 1.0000000 4.447679e-246 2.223839e-246
[80,] 1.0000000 3.891692e-246 1.945846e-246
[81,] 1.0000000 2.750131e-245 1.375066e-245
[82,] 1.0000000 1.183763e-244 5.918817e-245
[83,] 1.0000000 1.166529e-243 5.832646e-244
[84,] 1.0000000 1.450654e-242 7.253271e-243
[85,] 1.0000000 1.370590e-241 6.852949e-242
[86,] 1.0000000 1.833357e-242 9.166784e-243
[87,] 1.0000000 5.306553e-242 2.653277e-242
[88,] 1.0000000 4.338170e-243 2.169085e-243
[89,] 1.0000000 3.268875e-243 1.634438e-243
[90,] 1.0000000 3.915739e-242 1.957870e-242
[91,] 1.0000000 2.754674e-241 1.377337e-241
[92,] 1.0000000 1.327974e-240 6.639872e-241
[93,] 1.0000000 1.707697e-239 8.538486e-240
[94,] 1.0000000 1.763743e-238 8.818716e-239
[95,] 1.0000000 5.705230e-239 2.852615e-239
[96,] 1.0000000 4.298996e-238 2.149498e-238
[97,] 1.0000000 8.962346e-238 4.481173e-238
[98,] 1.0000000 4.929217e-237 2.464609e-237
[99,] 1.0000000 1.128768e-236 5.643840e-237
[100,] 1.0000000 2.713777e-237 1.356888e-237
[101,] 1.0000000 1.474139e-237 7.370694e-238
[102,] 1.0000000 8.989644e-237 4.494822e-237
[103,] 1.0000000 3.069808e-236 1.534904e-236
[104,] 1.0000000 2.801644e-236 1.400822e-236
[105,] 1.0000000 1.431704e-236 7.158518e-237
[106,] 1.0000000 3.081348e-238 1.540674e-238
[107,] 1.0000000 1.116047e-237 5.580237e-238
[108,] 1.0000000 4.254427e-237 2.127214e-237
[109,] 1.0000000 1.048865e-236 5.244324e-237
[110,] 1.0000000 8.603885e-236 4.301943e-236
[111,] 1.0000000 1.826082e-235 9.130412e-236
[112,] 1.0000000 5.788313e-237 2.894156e-237
[113,] 1.0000000 4.967795e-236 2.483897e-236
[114,] 1.0000000 6.263342e-236 3.131671e-236
[115,] 1.0000000 6.307007e-235 3.153503e-235
[116,] 1.0000000 5.573238e-234 2.786619e-234
[117,] 1.0000000 4.707670e-233 2.353835e-233
[118,] 1.0000000 5.924328e-232 2.962164e-232
[119,] 1.0000000 1.345509e-232 6.727543e-233
[120,] 1.0000000 1.487906e-231 7.439530e-232
[121,] 1.0000000 1.656035e-230 8.280173e-231
[122,] 1.0000000 1.720653e-229 8.603263e-230
[123,] 1.0000000 1.931005e-228 9.655027e-229
[124,] 1.0000000 2.164572e-227 1.082286e-227
[125,] 1.0000000 1.002074e-227 5.010368e-228
[126,] 1.0000000 1.149275e-226 5.746376e-227
[127,] 1.0000000 1.327898e-225 6.639489e-226
[128,] 1.0000000 1.529012e-224 7.645061e-225
[129,] 1.0000000 1.658907e-223 8.294537e-224
[130,] 1.0000000 1.908379e-222 9.541896e-223
[131,] 1.0000000 2.204916e-221 1.102458e-221
[132,] 1.0000000 2.548805e-220 1.274403e-220
[133,] 1.0000000 2.946973e-219 1.473487e-219
[134,] 1.0000000 3.299717e-218 1.649858e-218
[135,] 1.0000000 3.813818e-217 1.906909e-217
[136,] 1.0000000 3.287214e-216 1.643607e-216
[137,] 1.0000000 3.795535e-215 1.897768e-215
[138,] 1.0000000 3.752052e-214 1.876026e-214
[139,] 1.0000000 3.672779e-213 1.836390e-213
[140,] 1.0000000 4.233931e-212 2.116966e-212
[141,] 1.0000000 4.865316e-211 2.432658e-211
[142,] 1.0000000 5.591674e-210 2.795837e-210
[143,] 1.0000000 6.424036e-209 3.212018e-209
[144,] 1.0000000 7.251424e-208 3.625712e-208
[145,] 1.0000000 5.802227e-207 2.901113e-207
[146,] 1.0000000 6.621118e-206 3.310559e-206
[147,] 1.0000000 7.624161e-205 3.812080e-205
[148,] 1.0000000 8.662204e-204 4.331102e-204
[149,] 1.0000000 9.813242e-203 4.906621e-203
[150,] 1.0000000 1.108966e-201 5.544831e-202
[151,] 1.0000000 1.254566e-200 6.272828e-201
[152,] 1.0000000 1.185018e-199 5.925088e-200
[153,] 1.0000000 1.329196e-198 6.645981e-199
[154,] 1.0000000 1.486727e-197 7.433635e-198
[155,] 1.0000000 1.658135e-196 8.290677e-197
[156,] 1.0000000 1.843851e-195 9.219256e-196
[157,] 1.0000000 1.858139e-194 9.290694e-195
[158,] 1.0000000 2.058735e-193 1.029368e-193
[159,] 1.0000000 2.268535e-192 1.134267e-192
[160,] 1.0000000 2.491774e-191 1.245887e-191
[161,] 1.0000000 1.405656e-192 7.028278e-193
[162,] 1.0000000 1.566379e-191 7.831895e-192
[163,] 1.0000000 1.739600e-190 8.697999e-191
[164,] 1.0000000 1.925385e-189 9.626923e-190
[165,] 1.0000000 2.135505e-188 1.067753e-188
[166,] 1.0000000 2.347107e-187 1.173553e-187
[167,] 1.0000000 2.557155e-186 1.278578e-186
[168,] 1.0000000 2.790494e-185 1.395247e-185
[169,] 1.0000000 3.030016e-184 1.515008e-184
[170,] 1.0000000 3.282594e-183 1.641297e-183
[171,] 1.0000000 3.679499e-182 1.839749e-182
[172,] 1.0000000 3.952615e-181 1.976307e-181
[173,] 1.0000000 4.230357e-180 2.115178e-180
[174,] 1.0000000 2.357407e-179 1.178704e-179
[175,] 1.0000000 2.586125e-178 1.293063e-178
[176,] 1.0000000 1.685892e-177 8.429458e-178
[177,] 1.0000000 8.074023e-177 4.037011e-177
[178,] 1.0000000 8.867282e-176 4.433641e-176
[179,] 1.0000000 9.083079e-175 4.541540e-175
[180,] 1.0000000 7.528753e-174 3.764377e-174
[181,] 1.0000000 7.854207e-173 3.927104e-173
[182,] 1.0000000 6.544078e-172 3.272039e-172
[183,] 1.0000000 6.780985e-171 3.390493e-171
[184,] 1.0000000 6.387348e-170 3.193674e-170
[185,] 1.0000000 6.569784e-169 3.284892e-169
[186,] 1.0000000 6.576009e-168 3.288005e-168
[187,] 1.0000000 6.711888e-167 3.355944e-167
[188,] 1.0000000 6.822949e-166 3.411474e-166
[189,] 1.0000000 7.130121e-165 3.565060e-165
[190,] 1.0000000 2.617246e-164 1.308623e-164
[191,] 1.0000000 2.725197e-163 1.362598e-163
[192,] 1.0000000 2.824335e-162 1.412168e-162
[193,] 1.0000000 7.898688e-162 3.949344e-162
[194,] 1.0000000 8.164706e-161 4.082353e-161
[195,] 1.0000000 8.400210e-160 4.200105e-160
[196,] 1.0000000 8.602070e-159 4.301035e-159
[197,] 1.0000000 7.120288e-158 3.560144e-158
[198,] 1.0000000 1.373262e-157 6.866309e-158
[199,] 1.0000000 1.394830e-156 6.974148e-157
[200,] 1.0000000 1.312448e-155 6.562240e-156
[201,] 1.0000000 1.324345e-154 6.621726e-155
[202,] 1.0000000 9.818584e-155 4.909292e-155
[203,] 1.0000000 9.937227e-154 4.968614e-154
[204,] 1.0000000 9.680101e-153 4.840051e-153
[205,] 1.0000000 9.392578e-152 4.696289e-152
[206,] 1.0000000 8.083842e-151 4.041921e-151
[207,] 1.0000000 7.816793e-150 3.908396e-150
[208,] 1.0000000 7.576390e-149 3.788195e-149
[209,] 1.0000000 7.237398e-148 3.618699e-148
[210,] 1.0000000 6.848968e-147 3.424484e-147
[211,] 1.0000000 6.139933e-146 3.069966e-146
[212,] 1.0000000 5.795550e-145 2.897775e-145
[213,] 1.0000000 3.040793e-144 1.520396e-144
[214,] 1.0000000 2.647716e-143 1.323858e-143
[215,] 1.0000000 2.475822e-142 1.237911e-142
[216,] 1.0000000 2.302548e-141 1.151274e-141
[217,] 1.0000000 2.135280e-140 1.067640e-140
[218,] 1.0000000 1.969661e-139 9.848305e-140
[219,] 1.0000000 1.756347e-138 8.781733e-139
[220,] 1.0000000 1.608571e-137 8.042855e-138
[221,] 1.0000000 1.467017e-136 7.335087e-137
[222,] 1.0000000 1.261361e-135 6.306807e-136
[223,] 1.0000000 1.141241e-134 5.706205e-135
[224,] 1.0000000 1.026456e-133 5.132281e-134
[225,] 1.0000000 9.209166e-133 4.604583e-133
[226,] 1.0000000 8.226863e-132 4.113432e-132
[227,] 1.0000000 6.951779e-131 3.475890e-131
[228,] 1.0000000 6.157668e-130 3.078834e-130
[229,] 1.0000000 5.430229e-129 2.715114e-129
[230,] 1.0000000 4.768267e-128 2.384134e-128
[231,] 1.0000000 4.122384e-127 2.061192e-127
[232,] 1.0000000 3.588340e-126 1.794170e-126
[233,] 1.0000000 3.109754e-125 1.554877e-125
[234,] 1.0000000 2.683114e-124 1.341557e-124
[235,] 1.0000000 2.253351e-123 1.126675e-123
[236,] 1.0000000 1.927053e-122 9.635264e-123
[237,] 1.0000000 1.617826e-121 8.089131e-122
[238,] 1.0000000 1.376229e-120 6.881147e-121
[239,] 1.0000000 1.184151e-119 5.920753e-120
[240,] 1.0000000 1.014023e-118 5.070116e-119
[241,] 1.0000000 8.545465e-118 4.272732e-118
[242,] 1.0000000 7.248002e-117 3.624001e-117
[243,] 1.0000000 6.118058e-116 3.059029e-116
[244,] 1.0000000 4.654220e-115 2.327110e-115
[245,] 1.0000000 3.892974e-114 1.946487e-114
[246,] 1.0000000 3.098996e-113 1.549498e-113
[247,] 1.0000000 2.357610e-112 1.178805e-112
[248,] 1.0000000 1.130051e-111 5.650257e-112
[249,] 1.0000000 9.135237e-111 4.567619e-111
[250,] 1.0000000 4.865065e-110 2.432532e-110
[251,] 1.0000000 2.180001e-109 1.090001e-109
[252,] 1.0000000 1.780029e-108 8.900146e-109
[253,] 1.0000000 1.418969e-107 7.094844e-108
[254,] 1.0000000 1.147383e-106 5.736913e-107
[255,] 1.0000000 9.232096e-106 4.616048e-106
[256,] 1.0000000 7.093659e-105 3.546830e-105
[257,] 1.0000000 5.022811e-104 2.511405e-104
[258,] 1.0000000 3.985665e-103 1.992832e-103
[259,] 1.0000000 3.089676e-102 1.544838e-102
[260,] 1.0000000 2.427561e-101 1.213781e-101
[261,] 1.0000000 1.849520e-100 9.247598e-101
[262,] 1.0000000 1.438875e-99 7.194375e-100
[263,] 1.0000000 1.058158e-98 5.290790e-99
[264,] 1.0000000 8.152281e-98 4.076141e-98
[265,] 1.0000000 6.248897e-97 3.124449e-97
[266,] 1.0000000 4.765572e-96 2.382786e-96
[267,] 1.0000000 3.566662e-95 1.783331e-95
[268,] 1.0000000 2.647390e-94 1.323695e-94
[269,] 1.0000000 1.882924e-93 9.414622e-94
[270,] 1.0000000 1.345300e-92 6.726500e-93
[271,] 1.0000000 9.842506e-92 4.921253e-92
[272,] 1.0000000 7.072962e-91 3.536481e-91
[273,] 1.0000000 4.739310e-90 2.369655e-90
[274,] 1.0000000 3.475645e-89 1.737823e-89
[275,] 1.0000000 2.494409e-88 1.247205e-88
[276,] 1.0000000 1.523831e-87 7.619155e-88
[277,] 1.0000000 1.083975e-86 5.419873e-87
[278,] 1.0000000 6.945388e-86 3.472694e-86
[279,] 1.0000000 4.896810e-85 2.448405e-85
[280,] 1.0000000 3.488369e-84 1.744184e-84
[281,] 1.0000000 2.443689e-83 1.221844e-83
[282,] 1.0000000 1.722462e-82 8.612309e-83
[283,] 1.0000000 1.049381e-81 5.246904e-82
[284,] 1.0000000 7.326215e-81 3.663107e-81
[285,] 1.0000000 4.867411e-80 2.433705e-80
[286,] 1.0000000 3.201639e-79 1.600820e-79
[287,] 1.0000000 2.168561e-78 1.084281e-78
[288,] 1.0000000 1.456044e-77 7.280218e-78
[289,] 1.0000000 6.839108e-77 3.419554e-77
[290,] 1.0000000 4.571728e-76 2.285864e-76
[291,] 1.0000000 3.079540e-75 1.539770e-75
[292,] 1.0000000 2.036328e-74 1.018164e-74
[293,] 1.0000000 1.339189e-73 6.695944e-74
[294,] 1.0000000 8.879890e-73 4.439945e-73
[295,] 1.0000000 5.821210e-72 2.910605e-72
[296,] 1.0000000 3.765404e-71 1.882702e-71
[297,] 1.0000000 2.452459e-70 1.226230e-70
[298,] 1.0000000 1.162454e-69 5.812268e-70
[299,] 1.0000000 7.382815e-69 3.691408e-69
[300,] 1.0000000 4.741520e-68 2.370760e-68
[301,] 1.0000000 3.028420e-67 1.514210e-67
[302,] 1.0000000 1.475971e-66 7.379853e-67
[303,] 1.0000000 9.338041e-66 4.669021e-66
[304,] 1.0000000 5.872146e-65 2.936073e-65
[305,] 1.0000000 2.936485e-64 1.468243e-64
[306,] 1.0000000 1.825881e-63 9.129403e-64
[307,] 1.0000000 1.131785e-62 5.658926e-63
[308,] 1.0000000 6.863178e-62 3.431589e-62
[309,] 1.0000000 4.151007e-61 2.075504e-61
[310,] 1.0000000 2.470891e-60 1.235446e-60
[311,] 1.0000000 1.491235e-59 7.456173e-60
[312,] 1.0000000 8.828956e-59 4.414478e-59
[313,] 1.0000000 5.204148e-58 2.602074e-58
[314,] 1.0000000 3.048734e-57 1.524367e-57
[315,] 1.0000000 1.775007e-56 8.875034e-57
[316,] 1.0000000 1.027004e-55 5.135020e-56
[317,] 1.0000000 5.971449e-55 2.985725e-55
[318,] 1.0000000 3.411872e-54 1.705936e-54
[319,] 1.0000000 1.948093e-53 9.740467e-54
[320,] 1.0000000 1.110510e-52 5.552548e-53
[321,] 1.0000000 6.283688e-52 3.141844e-52
[322,] 1.0000000 3.499337e-51 1.749669e-51
[323,] 1.0000000 1.936023e-50 9.680116e-51
[324,] 1.0000000 1.064060e-49 5.320300e-50
[325,] 1.0000000 5.809375e-49 2.904688e-49
[326,] 1.0000000 3.150480e-48 1.575240e-48
[327,] 1.0000000 1.687921e-47 8.439605e-48
[328,] 1.0000000 9.027065e-47 4.513532e-47
[329,] 1.0000000 4.796323e-46 2.398162e-46
[330,] 1.0000000 2.536588e-45 1.268294e-45
[331,] 1.0000000 1.329130e-44 6.645648e-45
[332,] 1.0000000 6.917399e-44 3.458699e-44
[333,] 1.0000000 3.529195e-43 1.764597e-43
[334,] 1.0000000 1.834838e-42 9.174190e-43
[335,] 1.0000000 9.346627e-42 4.673314e-42
[336,] 1.0000000 4.726385e-41 2.363192e-41
[337,] 1.0000000 2.372407e-40 1.186204e-40
[338,] 1.0000000 1.181962e-39 5.909810e-40
[339,] 1.0000000 5.844391e-39 2.922195e-39
[340,] 1.0000000 2.899921e-38 1.449961e-38
[341,] 1.0000000 1.412231e-37 7.061154e-38
[342,] 1.0000000 6.824082e-37 3.412041e-37
[343,] 1.0000000 3.271637e-36 1.635819e-36
[344,] 1.0000000 1.513677e-35 7.568386e-36
[345,] 1.0000000 7.155764e-35 3.577882e-35
[346,] 1.0000000 3.037773e-34 1.518887e-34
[347,] 1.0000000 1.429665e-33 7.148323e-34
[348,] 1.0000000 6.586527e-33 3.293264e-33
[349,] 1.0000000 3.008965e-32 1.504482e-32
[350,] 1.0000000 1.274788e-31 6.373938e-32
[351,] 1.0000000 5.737952e-31 2.868976e-31
[352,] 1.0000000 2.560267e-30 1.280133e-30
[353,] 1.0000000 1.132330e-29 5.661650e-30
[354,] 1.0000000 4.963472e-29 2.481736e-29
[355,] 1.0000000 2.156024e-28 1.078012e-28
[356,] 1.0000000 9.279436e-28 4.639718e-28
[357,] 1.0000000 4.043815e-27 2.021908e-27
[358,] 1.0000000 1.704953e-26 8.524763e-27
[359,] 1.0000000 7.137058e-26 3.568529e-26
[360,] 1.0000000 2.958685e-25 1.479343e-25
[361,] 1.0000000 1.214464e-24 6.072322e-25
[362,] 1.0000000 4.935236e-24 2.467618e-24
[363,] 1.0000000 1.985162e-23 9.925808e-24
[364,] 1.0000000 7.902651e-23 3.951326e-23
[365,] 1.0000000 3.103054e-22 1.551527e-22
[366,] 1.0000000 1.190104e-21 5.950520e-22
[367,] 1.0000000 4.714785e-21 2.357393e-21
[368,] 1.0000000 1.806563e-20 9.032814e-21
[369,] 1.0000000 6.844780e-20 3.422390e-20
[370,] 1.0000000 2.506761e-19 1.253380e-19
[371,] 1.0000000 9.182415e-19 4.591207e-19
[372,] 1.0000000 3.360725e-18 1.680363e-18
[373,] 1.0000000 1.215114e-17 6.075571e-18
[374,] 1.0000000 4.338996e-17 2.169498e-17
[375,] 1.0000000 1.578238e-16 7.891192e-17
[376,] 1.0000000 4.187702e-16 2.093851e-16
[377,] 1.0000000 1.449459e-15 7.247296e-16
[378,] 1.0000000 4.981144e-15 2.490572e-15
[379,] 1.0000000 1.456507e-14 7.282534e-15
[380,] 1.0000000 4.766819e-14 2.383410e-14
[381,] 1.0000000 1.469654e-13 7.348271e-14
[382,] 1.0000000 4.790807e-13 2.395404e-13
[383,] 1.0000000 1.550189e-12 7.750946e-13
[384,] 1.0000000 4.786014e-12 2.393007e-12
[385,] 1.0000000 1.512044e-11 7.560220e-12
[386,] 1.0000000 4.800763e-11 2.400382e-11
[387,] 1.0000000 1.208213e-10 6.041065e-11
[388,] 1.0000000 2.583391e-10 1.291696e-10
[389,] 1.0000000 7.845136e-10 3.922568e-10
[390,] 1.0000000 2.118420e-09 1.059210e-09
[391,] 1.0000000 6.322213e-09 3.161106e-09
[392,] 1.0000000 1.832641e-08 9.163205e-09
[393,] 1.0000000 4.564867e-08 2.282434e-08
[394,] 0.9999999 1.274466e-07 6.372329e-08
[395,] 0.9999998 3.248774e-07 1.624387e-07
[396,] 0.9999996 8.473296e-07 4.236648e-07
[397,] 0.9999989 2.212691e-06 1.106345e-06
[398,] 0.9999971 5.793021e-06 2.896510e-06
[399,] 0.9999926 1.487337e-05 7.436683e-06
[400,] 0.9999821 3.578491e-05 1.789245e-05
[401,] 0.9999606 7.886807e-05 3.943403e-05
[402,] 0.9999093 1.814322e-04 9.071609e-05
[403,] 0.9997811 4.378506e-04 2.189253e-04
[404,] 0.9995127 9.746259e-04 4.873129e-04
[405,] 0.9989106 2.178836e-03 1.089418e-03
[406,] 0.9976088 4.782364e-03 2.391182e-03
[407,] 0.9949801 1.003978e-02 5.019890e-03
[408,] 0.9896587 2.068257e-02 1.034129e-02
[409,] 0.9794636 4.107281e-02 2.053640e-02
[410,] 0.9925859 1.482821e-02 7.414103e-03
[411,] 0.9843129 3.137426e-02 1.568713e-02
[412,] 0.9701958 5.960839e-02 2.980419e-02
[413,] 0.9392267 1.215467e-01 6.077333e-02
[414,] 0.8791527 2.416947e-01 1.208473e-01
[415,] 0.7984940 4.030121e-01 2.015060e-01
[416,] 0.6489841 7.020319e-01 3.510159e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/122zl1291400370.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/2vbzo1291400370.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/3vbzo1291400370.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/4vbzo1291400370.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/562gr1291400370.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
1614392.9083 3267339.8005 1342318.8812 -3145447.1126 -392660.3716
6 7 8 9 10
47198.0501 404.7394 602188.5762 -211744.7650 769347.7162
11 12 13 14 15
-324970.7246 476487.5596 971510.2238 -1780770.7218 52895.3919
16 17 18 19 20
440956.3216 -411358.2358 -67512.6769 -148150.4768 -328203.7754
21 22 23 24 25
173412.6685 -246220.5692 -237869.2592 -488859.5214 -214929.7252
26 27 28 29 30
-127584.6004 -365900.8069 -441606.0856 404846.4220 24692.1759
31 32 33 34 35
180623.9925 -162947.7947 132677.0574 16731.1667 -338033.7242
36 37 38 39 40
-454918.3637 55071.1870 73640.4234 -122998.4009 -149060.1825
41 42 43 44 45
23536.2045 -354630.1910 103738.8678 -186161.7635 -11873.8924
46 47 48 49 50
-232458.0978 -248121.3384 -171890.4711 -103951.1514 -156677.2080
51 52 53 54 55
25768.3763 -125901.9966 -405539.2242 -48526.5126 -394622.0142
56 57 58 59 60
134777.1577 182089.8545 -18512.2564 -20266.4483 -226136.2158
61 62 63 64 65
-89112.9730 -57073.6173 -109043.6851 -232356.0789 65003.2838
66 67 68 69 70
81678.9589 -91906.2811 -99661.3241 -58130.9141 -136377.9104
71 72 73 74 75
-28840.9276 73016.2831 -125491.3547 -14991.2288 -121785.5624
76 77 78 79 80
-141157.5865 -92252.6956 -47277.0402 87516.0975 -56250.2459
81 82 83 84 85
51772.1918 -98359.6297 -118980.2572 89337.8404 -5467.0157
86 87 88 89 90
-111554.4237 30263.8277 -39334.6250 -95805.4969 -20137.4147
91 92 93 94 95
-26562.2862 -102847.2821 136445.8672 -339389.1301 142088.0632
96 97 98 99 100
-22029.2939 12387.1874 47889.5549 -132866.7012 -85702.5943
101 102 103 104 105
-93984.4068 -101209.0829 -110607.5576 -53863.5304 59661.3416
106 107 108 109 110
-44169.0784 -239180.3303 -1079.6495 -203731.6259 -152145.4789
111 112 113 114 115
116244.2019 -234855.3016 -268536.7701 -174382.2503 -7876.3314
116 117 118 119 120
79594.8195 11626.7819 40981.2854 -237174.4564 77984.4846
121 122 123 124 125
-260529.3804 44126.2494 2437.0379 -48248.7533 -120768.4011
126 127 128 129 130
-249849.1739 56743.8463 57048.3791 64972.9222 56743.8463
131 132 133 134 135
59962.9990 114381.9051 56743.8463 54723.3658 56743.8463
136 137 138 139 140
59317.8032 59103.1536 56743.8463 56743.8463 56743.8463
141 142 143 144 145
55786.6489 56743.8463 51497.1868 57019.2654 63025.8180
146 147 148 149 150
68710.2635 56743.8463 58723.7837 56743.8463 56178.8235
151 152 153 154 155
47378.3959 46154.5855 56743.8463 34483.4455 56844.1304
156 157 158 159 160
56743.8463 56743.8463 53573.6402 63896.5616 56743.8463
161 162 163 164 165
56743.8463 56743.8463 56743.8463 62546.0270 57582.2679
166 167 168 169 170
56743.8463 56743.8463 160360.4989 56743.8463 56743.8463
171 172 173 174 175
56743.8463 53760.2551 56743.8463 54511.6315 56743.8463
176 177 178 179 180
56801.5725 56743.8463 -15353.1898 56743.8463 56743.8463
181 182 183 184 185
70989.7457 -37474.7381 82509.7404 -76005.3479 -14797.0131
186 187 188 189 190
-21688.9005 -8918.1814 53423.5954 17043.9417 56743.8463
191 192 193 194 195
61395.0315 56743.8463 58916.3004 56743.8463 56474.3785
196 197 198 199 200
-14102.2631 -150102.8416 -14102.2631 -14102.2631 -92005.6434
201 202 203 204 205
-14102.2631 -14102.2631 -14102.2631 14151.1709 -100326.1605
206 207 208 209 210
-15827.4682 -29262.9780 -14102.2631 44819.1609 -14102.2631
211 212 213 214 215
56743.8463 56743.8463 72687.1102 34552.6110 37988.8278
216 217 218 219 220
56743.8463 60320.3399 61069.6158 56743.8463 61629.1210
221 222 223 224 225
23935.1655 56743.8463 57658.7888 56743.8463 57060.4504
226 227 228 229 230
61626.0995 56743.8463 56743.8463 30572.5273 56743.8463
231 232 233 234 235
35628.8129 56743.8463 56743.8463 58191.0580 56743.8463
236 237 238 239 240
57007.1910 56743.8463 59739.9250 56743.8463 56743.8463
241 242 243 244 245
56743.8463 59635.8630 56743.8463 61991.0854 -30514.6252
246 247 248 249 250
-14102.2631 -14102.2631 48670.3819 -14102.2631 -14102.2631
251 252 253 254 255
63819.4336 -14102.2631 33344.1466 37600.1858 6243.3979
256 257 258 259 260
56743.8463 80097.7803 96741.5415 -14102.2631 56743.8463
261 262 263 264 265
-14102.2631 -14102.2631 37163.9577 65359.4770 -14102.2631
266 267 268 269 270
56743.8463 -14102.2631 60145.7968 -14102.2631 35248.8380
271 272 273 274 275
-14102.2631 -14102.2631 -14102.2631 55325.6036 56743.8463
276 277 278 279 280
66053.6292 63962.5695 56743.8463 -25065.7169 35917.4441
281 282 283 284 285
-14102.2631 56743.8463 95798.9808 56743.8463 20948.2632
286 287 288 289 290
56743.8463 -14102.2631 55796.2699 -14102.2631 87737.0651
291 292 293 294 295
-14102.2631 66657.4778 30150.0010 56743.8463 60459.4379
296 297 298 299 300
20666.8015 56743.8463 -14102.2631 56743.8463 56743.8463
301 302 303 304 305
-14799.5421 42657.9409 56743.8463 -14102.2631 96052.7763
306 307 308 309 310
60905.0179 -14102.2631 -22026.4417 78605.8994 -14102.2631
311 312 313 314 315
-14102.2631 37401.7567 -14102.2631 -17024.1271 58714.9646
316 317 318 319 320
57989.7061 58333.6778 42709.8089 59542.2511 56743.8463
321 322 323 324 325
56743.8463 56743.8463 56743.8463 -14102.2631 56743.8463
326 327 328 329 330
54184.7061 -14102.2631 -14102.2631 56743.8463 56743.8463
331 332 333 334 335
56743.8463 56743.8463 56743.8463 56785.5558 57641.3511
336 337 338 339 340
56743.8463 55622.0937 56743.8463 56098.6316 60825.7496
341 342 343 344 345
37739.2850 56743.8463 56743.8463 56743.8463 56743.8463
346 347 348 349 350
56743.8463 48733.4452 56743.8463 56743.8463 56743.8463
351 352 353 354 355
57507.4643 58135.5877 58182.2306 46857.2332 56743.8463
356 357 358 359 360
56743.8463 24153.8623 56743.8463 56743.8463 56743.8463
361 362 363 364 365
58592.2977 56743.8463 56743.8463 -23229.0752 55339.9264
366 367 368 369 370
56743.8463 56743.8463 56743.8463 56743.8463 56743.8463
371 372 373 374 375
56743.8463 59951.3846 16950.6168 29944.4766 56743.8463
376 377 378 379 380
56743.8463 69183.7439 41529.1872 56743.8463 56743.8463
381 382 383 384 385
56743.8463 50304.4797 110203.6465 56743.8463 58971.6081
386 387 388 389 390
-99802.4349 61423.9876 40725.4677 54932.2314 37401.4713
391 392 393 394 395
53556.4003 36181.4677 18933.7581 73314.7260 -225558.1571
396 397 398 399 400
29067.6811 -107954.0035 -1003.1517 30167.9153 -20626.9915
401 402 403 404 405
18801.9662 46542.0025 35992.8637 -96274.3599 23096.4155
406 407 408 409 410
4588.9031 41079.1019 36064.8764 34340.5635 -18627.8954
411 412 413 414 415
12187.3687 -219006.8707 -43752.6695 -154815.6019 -54634.8218
416 417 418 419 420
-37896.3342 138483.9631 -67743.4506 5703.9594 -11870.2290
421 422 423 424 425
-65172.7527 -99395.9490 -427200.0501 -470445.0446 -71049.1911
426 427 428 429 430
-200544.6914 -306223.7944 -1206200.7325 -359249.8705 -270057.1479
431
-843491.5827
> postscript(file="/var/www/html/freestat/rcomp/tmp/662gr1291400370.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 1614392.9083 NA
1 3267339.8005 1614392.9083
2 1342318.8812 3267339.8005
3 -3145447.1126 1342318.8812
4 -392660.3716 -3145447.1126
5 47198.0501 -392660.3716
6 404.7394 47198.0501
7 602188.5762 404.7394
8 -211744.7650 602188.5762
9 769347.7162 -211744.7650
10 -324970.7246 769347.7162
11 476487.5596 -324970.7246
12 971510.2238 476487.5596
13 -1780770.7218 971510.2238
14 52895.3919 -1780770.7218
15 440956.3216 52895.3919
16 -411358.2358 440956.3216
17 -67512.6769 -411358.2358
18 -148150.4768 -67512.6769
19 -328203.7754 -148150.4768
20 173412.6685 -328203.7754
21 -246220.5692 173412.6685
22 -237869.2592 -246220.5692
23 -488859.5214 -237869.2592
24 -214929.7252 -488859.5214
25 -127584.6004 -214929.7252
26 -365900.8069 -127584.6004
27 -441606.0856 -365900.8069
28 404846.4220 -441606.0856
29 24692.1759 404846.4220
30 180623.9925 24692.1759
31 -162947.7947 180623.9925
32 132677.0574 -162947.7947
33 16731.1667 132677.0574
34 -338033.7242 16731.1667
35 -454918.3637 -338033.7242
36 55071.1870 -454918.3637
37 73640.4234 55071.1870
38 -122998.4009 73640.4234
39 -149060.1825 -122998.4009
40 23536.2045 -149060.1825
41 -354630.1910 23536.2045
42 103738.8678 -354630.1910
43 -186161.7635 103738.8678
44 -11873.8924 -186161.7635
45 -232458.0978 -11873.8924
46 -248121.3384 -232458.0978
47 -171890.4711 -248121.3384
48 -103951.1514 -171890.4711
49 -156677.2080 -103951.1514
50 25768.3763 -156677.2080
51 -125901.9966 25768.3763
52 -405539.2242 -125901.9966
53 -48526.5126 -405539.2242
54 -394622.0142 -48526.5126
55 134777.1577 -394622.0142
56 182089.8545 134777.1577
57 -18512.2564 182089.8545
58 -20266.4483 -18512.2564
59 -226136.2158 -20266.4483
60 -89112.9730 -226136.2158
61 -57073.6173 -89112.9730
62 -109043.6851 -57073.6173
63 -232356.0789 -109043.6851
64 65003.2838 -232356.0789
65 81678.9589 65003.2838
66 -91906.2811 81678.9589
67 -99661.3241 -91906.2811
68 -58130.9141 -99661.3241
69 -136377.9104 -58130.9141
70 -28840.9276 -136377.9104
71 73016.2831 -28840.9276
72 -125491.3547 73016.2831
73 -14991.2288 -125491.3547
74 -121785.5624 -14991.2288
75 -141157.5865 -121785.5624
76 -92252.6956 -141157.5865
77 -47277.0402 -92252.6956
78 87516.0975 -47277.0402
79 -56250.2459 87516.0975
80 51772.1918 -56250.2459
81 -98359.6297 51772.1918
82 -118980.2572 -98359.6297
83 89337.8404 -118980.2572
84 -5467.0157 89337.8404
85 -111554.4237 -5467.0157
86 30263.8277 -111554.4237
87 -39334.6250 30263.8277
88 -95805.4969 -39334.6250
89 -20137.4147 -95805.4969
90 -26562.2862 -20137.4147
91 -102847.2821 -26562.2862
92 136445.8672 -102847.2821
93 -339389.1301 136445.8672
94 142088.0632 -339389.1301
95 -22029.2939 142088.0632
96 12387.1874 -22029.2939
97 47889.5549 12387.1874
98 -132866.7012 47889.5549
99 -85702.5943 -132866.7012
100 -93984.4068 -85702.5943
101 -101209.0829 -93984.4068
102 -110607.5576 -101209.0829
103 -53863.5304 -110607.5576
104 59661.3416 -53863.5304
105 -44169.0784 59661.3416
106 -239180.3303 -44169.0784
107 -1079.6495 -239180.3303
108 -203731.6259 -1079.6495
109 -152145.4789 -203731.6259
110 116244.2019 -152145.4789
111 -234855.3016 116244.2019
112 -268536.7701 -234855.3016
113 -174382.2503 -268536.7701
114 -7876.3314 -174382.2503
115 79594.8195 -7876.3314
116 11626.7819 79594.8195
117 40981.2854 11626.7819
118 -237174.4564 40981.2854
119 77984.4846 -237174.4564
120 -260529.3804 77984.4846
121 44126.2494 -260529.3804
122 2437.0379 44126.2494
123 -48248.7533 2437.0379
124 -120768.4011 -48248.7533
125 -249849.1739 -120768.4011
126 56743.8463 -249849.1739
127 57048.3791 56743.8463
128 64972.9222 57048.3791
129 56743.8463 64972.9222
130 59962.9990 56743.8463
131 114381.9051 59962.9990
132 56743.8463 114381.9051
133 54723.3658 56743.8463
134 56743.8463 54723.3658
135 59317.8032 56743.8463
136 59103.1536 59317.8032
137 56743.8463 59103.1536
138 56743.8463 56743.8463
139 56743.8463 56743.8463
140 55786.6489 56743.8463
141 56743.8463 55786.6489
142 51497.1868 56743.8463
143 57019.2654 51497.1868
144 63025.8180 57019.2654
145 68710.2635 63025.8180
146 56743.8463 68710.2635
147 58723.7837 56743.8463
148 56743.8463 58723.7837
149 56178.8235 56743.8463
150 47378.3959 56178.8235
151 46154.5855 47378.3959
152 56743.8463 46154.5855
153 34483.4455 56743.8463
154 56844.1304 34483.4455
155 56743.8463 56844.1304
156 56743.8463 56743.8463
157 53573.6402 56743.8463
158 63896.5616 53573.6402
159 56743.8463 63896.5616
160 56743.8463 56743.8463
161 56743.8463 56743.8463
162 56743.8463 56743.8463
163 62546.0270 56743.8463
164 57582.2679 62546.0270
165 56743.8463 57582.2679
166 56743.8463 56743.8463
167 160360.4989 56743.8463
168 56743.8463 160360.4989
169 56743.8463 56743.8463
170 56743.8463 56743.8463
171 53760.2551 56743.8463
172 56743.8463 53760.2551
173 54511.6315 56743.8463
174 56743.8463 54511.6315
175 56801.5725 56743.8463
176 56743.8463 56801.5725
177 -15353.1898 56743.8463
178 56743.8463 -15353.1898
179 56743.8463 56743.8463
180 70989.7457 56743.8463
181 -37474.7381 70989.7457
182 82509.7404 -37474.7381
183 -76005.3479 82509.7404
184 -14797.0131 -76005.3479
185 -21688.9005 -14797.0131
186 -8918.1814 -21688.9005
187 53423.5954 -8918.1814
188 17043.9417 53423.5954
189 56743.8463 17043.9417
190 61395.0315 56743.8463
191 56743.8463 61395.0315
192 58916.3004 56743.8463
193 56743.8463 58916.3004
194 56474.3785 56743.8463
195 -14102.2631 56474.3785
196 -150102.8416 -14102.2631
197 -14102.2631 -150102.8416
198 -14102.2631 -14102.2631
199 -92005.6434 -14102.2631
200 -14102.2631 -92005.6434
201 -14102.2631 -14102.2631
202 -14102.2631 -14102.2631
203 14151.1709 -14102.2631
204 -100326.1605 14151.1709
205 -15827.4682 -100326.1605
206 -29262.9780 -15827.4682
207 -14102.2631 -29262.9780
208 44819.1609 -14102.2631
209 -14102.2631 44819.1609
210 56743.8463 -14102.2631
211 56743.8463 56743.8463
212 72687.1102 56743.8463
213 34552.6110 72687.1102
214 37988.8278 34552.6110
215 56743.8463 37988.8278
216 60320.3399 56743.8463
217 61069.6158 60320.3399
218 56743.8463 61069.6158
219 61629.1210 56743.8463
220 23935.1655 61629.1210
221 56743.8463 23935.1655
222 57658.7888 56743.8463
223 56743.8463 57658.7888
224 57060.4504 56743.8463
225 61626.0995 57060.4504
226 56743.8463 61626.0995
227 56743.8463 56743.8463
228 30572.5273 56743.8463
229 56743.8463 30572.5273
230 35628.8129 56743.8463
231 56743.8463 35628.8129
232 56743.8463 56743.8463
233 58191.0580 56743.8463
234 56743.8463 58191.0580
235 57007.1910 56743.8463
236 56743.8463 57007.1910
237 59739.9250 56743.8463
238 56743.8463 59739.9250
239 56743.8463 56743.8463
240 56743.8463 56743.8463
241 59635.8630 56743.8463
242 56743.8463 59635.8630
243 61991.0854 56743.8463
244 -30514.6252 61991.0854
245 -14102.2631 -30514.6252
246 -14102.2631 -14102.2631
247 48670.3819 -14102.2631
248 -14102.2631 48670.3819
249 -14102.2631 -14102.2631
250 63819.4336 -14102.2631
251 -14102.2631 63819.4336
252 33344.1466 -14102.2631
253 37600.1858 33344.1466
254 6243.3979 37600.1858
255 56743.8463 6243.3979
256 80097.7803 56743.8463
257 96741.5415 80097.7803
258 -14102.2631 96741.5415
259 56743.8463 -14102.2631
260 -14102.2631 56743.8463
261 -14102.2631 -14102.2631
262 37163.9577 -14102.2631
263 65359.4770 37163.9577
264 -14102.2631 65359.4770
265 56743.8463 -14102.2631
266 -14102.2631 56743.8463
267 60145.7968 -14102.2631
268 -14102.2631 60145.7968
269 35248.8380 -14102.2631
270 -14102.2631 35248.8380
271 -14102.2631 -14102.2631
272 -14102.2631 -14102.2631
273 55325.6036 -14102.2631
274 56743.8463 55325.6036
275 66053.6292 56743.8463
276 63962.5695 66053.6292
277 56743.8463 63962.5695
278 -25065.7169 56743.8463
279 35917.4441 -25065.7169
280 -14102.2631 35917.4441
281 56743.8463 -14102.2631
282 95798.9808 56743.8463
283 56743.8463 95798.9808
284 20948.2632 56743.8463
285 56743.8463 20948.2632
286 -14102.2631 56743.8463
287 55796.2699 -14102.2631
288 -14102.2631 55796.2699
289 87737.0651 -14102.2631
290 -14102.2631 87737.0651
291 66657.4778 -14102.2631
292 30150.0010 66657.4778
293 56743.8463 30150.0010
294 60459.4379 56743.8463
295 20666.8015 60459.4379
296 56743.8463 20666.8015
297 -14102.2631 56743.8463
298 56743.8463 -14102.2631
299 56743.8463 56743.8463
300 -14799.5421 56743.8463
301 42657.9409 -14799.5421
302 56743.8463 42657.9409
303 -14102.2631 56743.8463
304 96052.7763 -14102.2631
305 60905.0179 96052.7763
306 -14102.2631 60905.0179
307 -22026.4417 -14102.2631
308 78605.8994 -22026.4417
309 -14102.2631 78605.8994
310 -14102.2631 -14102.2631
311 37401.7567 -14102.2631
312 -14102.2631 37401.7567
313 -17024.1271 -14102.2631
314 58714.9646 -17024.1271
315 57989.7061 58714.9646
316 58333.6778 57989.7061
317 42709.8089 58333.6778
318 59542.2511 42709.8089
319 56743.8463 59542.2511
320 56743.8463 56743.8463
321 56743.8463 56743.8463
322 56743.8463 56743.8463
323 -14102.2631 56743.8463
324 56743.8463 -14102.2631
325 54184.7061 56743.8463
326 -14102.2631 54184.7061
327 -14102.2631 -14102.2631
328 56743.8463 -14102.2631
329 56743.8463 56743.8463
330 56743.8463 56743.8463
331 56743.8463 56743.8463
332 56743.8463 56743.8463
333 56785.5558 56743.8463
334 57641.3511 56785.5558
335 56743.8463 57641.3511
336 55622.0937 56743.8463
337 56743.8463 55622.0937
338 56098.6316 56743.8463
339 60825.7496 56098.6316
340 37739.2850 60825.7496
341 56743.8463 37739.2850
342 56743.8463 56743.8463
343 56743.8463 56743.8463
344 56743.8463 56743.8463
345 56743.8463 56743.8463
346 48733.4452 56743.8463
347 56743.8463 48733.4452
348 56743.8463 56743.8463
349 56743.8463 56743.8463
350 57507.4643 56743.8463
351 58135.5877 57507.4643
352 58182.2306 58135.5877
353 46857.2332 58182.2306
354 56743.8463 46857.2332
355 56743.8463 56743.8463
356 24153.8623 56743.8463
357 56743.8463 24153.8623
358 56743.8463 56743.8463
359 56743.8463 56743.8463
360 58592.2977 56743.8463
361 56743.8463 58592.2977
362 56743.8463 56743.8463
363 -23229.0752 56743.8463
364 55339.9264 -23229.0752
365 56743.8463 55339.9264
366 56743.8463 56743.8463
367 56743.8463 56743.8463
368 56743.8463 56743.8463
369 56743.8463 56743.8463
370 56743.8463 56743.8463
371 59951.3846 56743.8463
372 16950.6168 59951.3846
373 29944.4766 16950.6168
374 56743.8463 29944.4766
375 56743.8463 56743.8463
376 69183.7439 56743.8463
377 41529.1872 69183.7439
378 56743.8463 41529.1872
379 56743.8463 56743.8463
380 56743.8463 56743.8463
381 50304.4797 56743.8463
382 110203.6465 50304.4797
383 56743.8463 110203.6465
384 58971.6081 56743.8463
385 -99802.4349 58971.6081
386 61423.9876 -99802.4349
387 40725.4677 61423.9876
388 54932.2314 40725.4677
389 37401.4713 54932.2314
390 53556.4003 37401.4713
391 36181.4677 53556.4003
392 18933.7581 36181.4677
393 73314.7260 18933.7581
394 -225558.1571 73314.7260
395 29067.6811 -225558.1571
396 -107954.0035 29067.6811
397 -1003.1517 -107954.0035
398 30167.9153 -1003.1517
399 -20626.9915 30167.9153
400 18801.9662 -20626.9915
401 46542.0025 18801.9662
402 35992.8637 46542.0025
403 -96274.3599 35992.8637
404 23096.4155 -96274.3599
405 4588.9031 23096.4155
406 41079.1019 4588.9031
407 36064.8764 41079.1019
408 34340.5635 36064.8764
409 -18627.8954 34340.5635
410 12187.3687 -18627.8954
411 -219006.8707 12187.3687
412 -43752.6695 -219006.8707
413 -154815.6019 -43752.6695
414 -54634.8218 -154815.6019
415 -37896.3342 -54634.8218
416 138483.9631 -37896.3342
417 -67743.4506 138483.9631
418 5703.9594 -67743.4506
419 -11870.2290 5703.9594
420 -65172.7527 -11870.2290
421 -99395.9490 -65172.7527
422 -427200.0501 -99395.9490
423 -470445.0446 -427200.0501
424 -71049.1911 -470445.0446
425 -200544.6914 -71049.1911
426 -306223.7944 -200544.6914
427 -1206200.7325 -306223.7944
428 -359249.8705 -1206200.7325
429 -270057.1479 -359249.8705
430 -843491.5827 -270057.1479
431 NA -843491.5827
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3267339.8005 1614392.9083
[2,] 1342318.8812 3267339.8005
[3,] -3145447.1126 1342318.8812
[4,] -392660.3716 -3145447.1126
[5,] 47198.0501 -392660.3716
[6,] 404.7394 47198.0501
[7,] 602188.5762 404.7394
[8,] -211744.7650 602188.5762
[9,] 769347.7162 -211744.7650
[10,] -324970.7246 769347.7162
[11,] 476487.5596 -324970.7246
[12,] 971510.2238 476487.5596
[13,] -1780770.7218 971510.2238
[14,] 52895.3919 -1780770.7218
[15,] 440956.3216 52895.3919
[16,] -411358.2358 440956.3216
[17,] -67512.6769 -411358.2358
[18,] -148150.4768 -67512.6769
[19,] -328203.7754 -148150.4768
[20,] 173412.6685 -328203.7754
[21,] -246220.5692 173412.6685
[22,] -237869.2592 -246220.5692
[23,] -488859.5214 -237869.2592
[24,] -214929.7252 -488859.5214
[25,] -127584.6004 -214929.7252
[26,] -365900.8069 -127584.6004
[27,] -441606.0856 -365900.8069
[28,] 404846.4220 -441606.0856
[29,] 24692.1759 404846.4220
[30,] 180623.9925 24692.1759
[31,] -162947.7947 180623.9925
[32,] 132677.0574 -162947.7947
[33,] 16731.1667 132677.0574
[34,] -338033.7242 16731.1667
[35,] -454918.3637 -338033.7242
[36,] 55071.1870 -454918.3637
[37,] 73640.4234 55071.1870
[38,] -122998.4009 73640.4234
[39,] -149060.1825 -122998.4009
[40,] 23536.2045 -149060.1825
[41,] -354630.1910 23536.2045
[42,] 103738.8678 -354630.1910
[43,] -186161.7635 103738.8678
[44,] -11873.8924 -186161.7635
[45,] -232458.0978 -11873.8924
[46,] -248121.3384 -232458.0978
[47,] -171890.4711 -248121.3384
[48,] -103951.1514 -171890.4711
[49,] -156677.2080 -103951.1514
[50,] 25768.3763 -156677.2080
[51,] -125901.9966 25768.3763
[52,] -405539.2242 -125901.9966
[53,] -48526.5126 -405539.2242
[54,] -394622.0142 -48526.5126
[55,] 134777.1577 -394622.0142
[56,] 182089.8545 134777.1577
[57,] -18512.2564 182089.8545
[58,] -20266.4483 -18512.2564
[59,] -226136.2158 -20266.4483
[60,] -89112.9730 -226136.2158
[61,] -57073.6173 -89112.9730
[62,] -109043.6851 -57073.6173
[63,] -232356.0789 -109043.6851
[64,] 65003.2838 -232356.0789
[65,] 81678.9589 65003.2838
[66,] -91906.2811 81678.9589
[67,] -99661.3241 -91906.2811
[68,] -58130.9141 -99661.3241
[69,] -136377.9104 -58130.9141
[70,] -28840.9276 -136377.9104
[71,] 73016.2831 -28840.9276
[72,] -125491.3547 73016.2831
[73,] -14991.2288 -125491.3547
[74,] -121785.5624 -14991.2288
[75,] -141157.5865 -121785.5624
[76,] -92252.6956 -141157.5865
[77,] -47277.0402 -92252.6956
[78,] 87516.0975 -47277.0402
[79,] -56250.2459 87516.0975
[80,] 51772.1918 -56250.2459
[81,] -98359.6297 51772.1918
[82,] -118980.2572 -98359.6297
[83,] 89337.8404 -118980.2572
[84,] -5467.0157 89337.8404
[85,] -111554.4237 -5467.0157
[86,] 30263.8277 -111554.4237
[87,] -39334.6250 30263.8277
[88,] -95805.4969 -39334.6250
[89,] -20137.4147 -95805.4969
[90,] -26562.2862 -20137.4147
[91,] -102847.2821 -26562.2862
[92,] 136445.8672 -102847.2821
[93,] -339389.1301 136445.8672
[94,] 142088.0632 -339389.1301
[95,] -22029.2939 142088.0632
[96,] 12387.1874 -22029.2939
[97,] 47889.5549 12387.1874
[98,] -132866.7012 47889.5549
[99,] -85702.5943 -132866.7012
[100,] -93984.4068 -85702.5943
[101,] -101209.0829 -93984.4068
[102,] -110607.5576 -101209.0829
[103,] -53863.5304 -110607.5576
[104,] 59661.3416 -53863.5304
[105,] -44169.0784 59661.3416
[106,] -239180.3303 -44169.0784
[107,] -1079.6495 -239180.3303
[108,] -203731.6259 -1079.6495
[109,] -152145.4789 -203731.6259
[110,] 116244.2019 -152145.4789
[111,] -234855.3016 116244.2019
[112,] -268536.7701 -234855.3016
[113,] -174382.2503 -268536.7701
[114,] -7876.3314 -174382.2503
[115,] 79594.8195 -7876.3314
[116,] 11626.7819 79594.8195
[117,] 40981.2854 11626.7819
[118,] -237174.4564 40981.2854
[119,] 77984.4846 -237174.4564
[120,] -260529.3804 77984.4846
[121,] 44126.2494 -260529.3804
[122,] 2437.0379 44126.2494
[123,] -48248.7533 2437.0379
[124,] -120768.4011 -48248.7533
[125,] -249849.1739 -120768.4011
[126,] 56743.8463 -249849.1739
[127,] 57048.3791 56743.8463
[128,] 64972.9222 57048.3791
[129,] 56743.8463 64972.9222
[130,] 59962.9990 56743.8463
[131,] 114381.9051 59962.9990
[132,] 56743.8463 114381.9051
[133,] 54723.3658 56743.8463
[134,] 56743.8463 54723.3658
[135,] 59317.8032 56743.8463
[136,] 59103.1536 59317.8032
[137,] 56743.8463 59103.1536
[138,] 56743.8463 56743.8463
[139,] 56743.8463 56743.8463
[140,] 55786.6489 56743.8463
[141,] 56743.8463 55786.6489
[142,] 51497.1868 56743.8463
[143,] 57019.2654 51497.1868
[144,] 63025.8180 57019.2654
[145,] 68710.2635 63025.8180
[146,] 56743.8463 68710.2635
[147,] 58723.7837 56743.8463
[148,] 56743.8463 58723.7837
[149,] 56178.8235 56743.8463
[150,] 47378.3959 56178.8235
[151,] 46154.5855 47378.3959
[152,] 56743.8463 46154.5855
[153,] 34483.4455 56743.8463
[154,] 56844.1304 34483.4455
[155,] 56743.8463 56844.1304
[156,] 56743.8463 56743.8463
[157,] 53573.6402 56743.8463
[158,] 63896.5616 53573.6402
[159,] 56743.8463 63896.5616
[160,] 56743.8463 56743.8463
[161,] 56743.8463 56743.8463
[162,] 56743.8463 56743.8463
[163,] 62546.0270 56743.8463
[164,] 57582.2679 62546.0270
[165,] 56743.8463 57582.2679
[166,] 56743.8463 56743.8463
[167,] 160360.4989 56743.8463
[168,] 56743.8463 160360.4989
[169,] 56743.8463 56743.8463
[170,] 56743.8463 56743.8463
[171,] 53760.2551 56743.8463
[172,] 56743.8463 53760.2551
[173,] 54511.6315 56743.8463
[174,] 56743.8463 54511.6315
[175,] 56801.5725 56743.8463
[176,] 56743.8463 56801.5725
[177,] -15353.1898 56743.8463
[178,] 56743.8463 -15353.1898
[179,] 56743.8463 56743.8463
[180,] 70989.7457 56743.8463
[181,] -37474.7381 70989.7457
[182,] 82509.7404 -37474.7381
[183,] -76005.3479 82509.7404
[184,] -14797.0131 -76005.3479
[185,] -21688.9005 -14797.0131
[186,] -8918.1814 -21688.9005
[187,] 53423.5954 -8918.1814
[188,] 17043.9417 53423.5954
[189,] 56743.8463 17043.9417
[190,] 61395.0315 56743.8463
[191,] 56743.8463 61395.0315
[192,] 58916.3004 56743.8463
[193,] 56743.8463 58916.3004
[194,] 56474.3785 56743.8463
[195,] -14102.2631 56474.3785
[196,] -150102.8416 -14102.2631
[197,] -14102.2631 -150102.8416
[198,] -14102.2631 -14102.2631
[199,] -92005.6434 -14102.2631
[200,] -14102.2631 -92005.6434
[201,] -14102.2631 -14102.2631
[202,] -14102.2631 -14102.2631
[203,] 14151.1709 -14102.2631
[204,] -100326.1605 14151.1709
[205,] -15827.4682 -100326.1605
[206,] -29262.9780 -15827.4682
[207,] -14102.2631 -29262.9780
[208,] 44819.1609 -14102.2631
[209,] -14102.2631 44819.1609
[210,] 56743.8463 -14102.2631
[211,] 56743.8463 56743.8463
[212,] 72687.1102 56743.8463
[213,] 34552.6110 72687.1102
[214,] 37988.8278 34552.6110
[215,] 56743.8463 37988.8278
[216,] 60320.3399 56743.8463
[217,] 61069.6158 60320.3399
[218,] 56743.8463 61069.6158
[219,] 61629.1210 56743.8463
[220,] 23935.1655 61629.1210
[221,] 56743.8463 23935.1655
[222,] 57658.7888 56743.8463
[223,] 56743.8463 57658.7888
[224,] 57060.4504 56743.8463
[225,] 61626.0995 57060.4504
[226,] 56743.8463 61626.0995
[227,] 56743.8463 56743.8463
[228,] 30572.5273 56743.8463
[229,] 56743.8463 30572.5273
[230,] 35628.8129 56743.8463
[231,] 56743.8463 35628.8129
[232,] 56743.8463 56743.8463
[233,] 58191.0580 56743.8463
[234,] 56743.8463 58191.0580
[235,] 57007.1910 56743.8463
[236,] 56743.8463 57007.1910
[237,] 59739.9250 56743.8463
[238,] 56743.8463 59739.9250
[239,] 56743.8463 56743.8463
[240,] 56743.8463 56743.8463
[241,] 59635.8630 56743.8463
[242,] 56743.8463 59635.8630
[243,] 61991.0854 56743.8463
[244,] -30514.6252 61991.0854
[245,] -14102.2631 -30514.6252
[246,] -14102.2631 -14102.2631
[247,] 48670.3819 -14102.2631
[248,] -14102.2631 48670.3819
[249,] -14102.2631 -14102.2631
[250,] 63819.4336 -14102.2631
[251,] -14102.2631 63819.4336
[252,] 33344.1466 -14102.2631
[253,] 37600.1858 33344.1466
[254,] 6243.3979 37600.1858
[255,] 56743.8463 6243.3979
[256,] 80097.7803 56743.8463
[257,] 96741.5415 80097.7803
[258,] -14102.2631 96741.5415
[259,] 56743.8463 -14102.2631
[260,] -14102.2631 56743.8463
[261,] -14102.2631 -14102.2631
[262,] 37163.9577 -14102.2631
[263,] 65359.4770 37163.9577
[264,] -14102.2631 65359.4770
[265,] 56743.8463 -14102.2631
[266,] -14102.2631 56743.8463
[267,] 60145.7968 -14102.2631
[268,] -14102.2631 60145.7968
[269,] 35248.8380 -14102.2631
[270,] -14102.2631 35248.8380
[271,] -14102.2631 -14102.2631
[272,] -14102.2631 -14102.2631
[273,] 55325.6036 -14102.2631
[274,] 56743.8463 55325.6036
[275,] 66053.6292 56743.8463
[276,] 63962.5695 66053.6292
[277,] 56743.8463 63962.5695
[278,] -25065.7169 56743.8463
[279,] 35917.4441 -25065.7169
[280,] -14102.2631 35917.4441
[281,] 56743.8463 -14102.2631
[282,] 95798.9808 56743.8463
[283,] 56743.8463 95798.9808
[284,] 20948.2632 56743.8463
[285,] 56743.8463 20948.2632
[286,] -14102.2631 56743.8463
[287,] 55796.2699 -14102.2631
[288,] -14102.2631 55796.2699
[289,] 87737.0651 -14102.2631
[290,] -14102.2631 87737.0651
[291,] 66657.4778 -14102.2631
[292,] 30150.0010 66657.4778
[293,] 56743.8463 30150.0010
[294,] 60459.4379 56743.8463
[295,] 20666.8015 60459.4379
[296,] 56743.8463 20666.8015
[297,] -14102.2631 56743.8463
[298,] 56743.8463 -14102.2631
[299,] 56743.8463 56743.8463
[300,] -14799.5421 56743.8463
[301,] 42657.9409 -14799.5421
[302,] 56743.8463 42657.9409
[303,] -14102.2631 56743.8463
[304,] 96052.7763 -14102.2631
[305,] 60905.0179 96052.7763
[306,] -14102.2631 60905.0179
[307,] -22026.4417 -14102.2631
[308,] 78605.8994 -22026.4417
[309,] -14102.2631 78605.8994
[310,] -14102.2631 -14102.2631
[311,] 37401.7567 -14102.2631
[312,] -14102.2631 37401.7567
[313,] -17024.1271 -14102.2631
[314,] 58714.9646 -17024.1271
[315,] 57989.7061 58714.9646
[316,] 58333.6778 57989.7061
[317,] 42709.8089 58333.6778
[318,] 59542.2511 42709.8089
[319,] 56743.8463 59542.2511
[320,] 56743.8463 56743.8463
[321,] 56743.8463 56743.8463
[322,] 56743.8463 56743.8463
[323,] -14102.2631 56743.8463
[324,] 56743.8463 -14102.2631
[325,] 54184.7061 56743.8463
[326,] -14102.2631 54184.7061
[327,] -14102.2631 -14102.2631
[328,] 56743.8463 -14102.2631
[329,] 56743.8463 56743.8463
[330,] 56743.8463 56743.8463
[331,] 56743.8463 56743.8463
[332,] 56743.8463 56743.8463
[333,] 56785.5558 56743.8463
[334,] 57641.3511 56785.5558
[335,] 56743.8463 57641.3511
[336,] 55622.0937 56743.8463
[337,] 56743.8463 55622.0937
[338,] 56098.6316 56743.8463
[339,] 60825.7496 56098.6316
[340,] 37739.2850 60825.7496
[341,] 56743.8463 37739.2850
[342,] 56743.8463 56743.8463
[343,] 56743.8463 56743.8463
[344,] 56743.8463 56743.8463
[345,] 56743.8463 56743.8463
[346,] 48733.4452 56743.8463
[347,] 56743.8463 48733.4452
[348,] 56743.8463 56743.8463
[349,] 56743.8463 56743.8463
[350,] 57507.4643 56743.8463
[351,] 58135.5877 57507.4643
[352,] 58182.2306 58135.5877
[353,] 46857.2332 58182.2306
[354,] 56743.8463 46857.2332
[355,] 56743.8463 56743.8463
[356,] 24153.8623 56743.8463
[357,] 56743.8463 24153.8623
[358,] 56743.8463 56743.8463
[359,] 56743.8463 56743.8463
[360,] 58592.2977 56743.8463
[361,] 56743.8463 58592.2977
[362,] 56743.8463 56743.8463
[363,] -23229.0752 56743.8463
[364,] 55339.9264 -23229.0752
[365,] 56743.8463 55339.9264
[366,] 56743.8463 56743.8463
[367,] 56743.8463 56743.8463
[368,] 56743.8463 56743.8463
[369,] 56743.8463 56743.8463
[370,] 56743.8463 56743.8463
[371,] 59951.3846 56743.8463
[372,] 16950.6168 59951.3846
[373,] 29944.4766 16950.6168
[374,] 56743.8463 29944.4766
[375,] 56743.8463 56743.8463
[376,] 69183.7439 56743.8463
[377,] 41529.1872 69183.7439
[378,] 56743.8463 41529.1872
[379,] 56743.8463 56743.8463
[380,] 56743.8463 56743.8463
[381,] 50304.4797 56743.8463
[382,] 110203.6465 50304.4797
[383,] 56743.8463 110203.6465
[384,] 58971.6081 56743.8463
[385,] -99802.4349 58971.6081
[386,] 61423.9876 -99802.4349
[387,] 40725.4677 61423.9876
[388,] 54932.2314 40725.4677
[389,] 37401.4713 54932.2314
[390,] 53556.4003 37401.4713
[391,] 36181.4677 53556.4003
[392,] 18933.7581 36181.4677
[393,] 73314.7260 18933.7581
[394,] -225558.1571 73314.7260
[395,] 29067.6811 -225558.1571
[396,] -107954.0035 29067.6811
[397,] -1003.1517 -107954.0035
[398,] 30167.9153 -1003.1517
[399,] -20626.9915 30167.9153
[400,] 18801.9662 -20626.9915
[401,] 46542.0025 18801.9662
[402,] 35992.8637 46542.0025
[403,] -96274.3599 35992.8637
[404,] 23096.4155 -96274.3599
[405,] 4588.9031 23096.4155
[406,] 41079.1019 4588.9031
[407,] 36064.8764 41079.1019
[408,] 34340.5635 36064.8764
[409,] -18627.8954 34340.5635
[410,] 12187.3687 -18627.8954
[411,] -219006.8707 12187.3687
[412,] -43752.6695 -219006.8707
[413,] -154815.6019 -43752.6695
[414,] -54634.8218 -154815.6019
[415,] -37896.3342 -54634.8218
[416,] 138483.9631 -37896.3342
[417,] -67743.4506 138483.9631
[418,] 5703.9594 -67743.4506
[419,] -11870.2290 5703.9594
[420,] -65172.7527 -11870.2290
[421,] -99395.9490 -65172.7527
[422,] -427200.0501 -99395.9490
[423,] -470445.0446 -427200.0501
[424,] -71049.1911 -470445.0446
[425,] -200544.6914 -71049.1911
[426,] -306223.7944 -200544.6914
[427,] -1206200.7325 -306223.7944
[428,] -359249.8705 -1206200.7325
[429,] -270057.1479 -359249.8705
[430,] -843491.5827 -270057.1479
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3267339.8005 1614392.9083
2 1342318.8812 3267339.8005
3 -3145447.1126 1342318.8812
4 -392660.3716 -3145447.1126
5 47198.0501 -392660.3716
6 404.7394 47198.0501
7 602188.5762 404.7394
8 -211744.7650 602188.5762
9 769347.7162 -211744.7650
10 -324970.7246 769347.7162
11 476487.5596 -324970.7246
12 971510.2238 476487.5596
13 -1780770.7218 971510.2238
14 52895.3919 -1780770.7218
15 440956.3216 52895.3919
16 -411358.2358 440956.3216
17 -67512.6769 -411358.2358
18 -148150.4768 -67512.6769
19 -328203.7754 -148150.4768
20 173412.6685 -328203.7754
21 -246220.5692 173412.6685
22 -237869.2592 -246220.5692
23 -488859.5214 -237869.2592
24 -214929.7252 -488859.5214
25 -127584.6004 -214929.7252
26 -365900.8069 -127584.6004
27 -441606.0856 -365900.8069
28 404846.4220 -441606.0856
29 24692.1759 404846.4220
30 180623.9925 24692.1759
31 -162947.7947 180623.9925
32 132677.0574 -162947.7947
33 16731.1667 132677.0574
34 -338033.7242 16731.1667
35 -454918.3637 -338033.7242
36 55071.1870 -454918.3637
37 73640.4234 55071.1870
38 -122998.4009 73640.4234
39 -149060.1825 -122998.4009
40 23536.2045 -149060.1825
41 -354630.1910 23536.2045
42 103738.8678 -354630.1910
43 -186161.7635 103738.8678
44 -11873.8924 -186161.7635
45 -232458.0978 -11873.8924
46 -248121.3384 -232458.0978
47 -171890.4711 -248121.3384
48 -103951.1514 -171890.4711
49 -156677.2080 -103951.1514
50 25768.3763 -156677.2080
51 -125901.9966 25768.3763
52 -405539.2242 -125901.9966
53 -48526.5126 -405539.2242
54 -394622.0142 -48526.5126
55 134777.1577 -394622.0142
56 182089.8545 134777.1577
57 -18512.2564 182089.8545
58 -20266.4483 -18512.2564
59 -226136.2158 -20266.4483
60 -89112.9730 -226136.2158
61 -57073.6173 -89112.9730
62 -109043.6851 -57073.6173
63 -232356.0789 -109043.6851
64 65003.2838 -232356.0789
65 81678.9589 65003.2838
66 -91906.2811 81678.9589
67 -99661.3241 -91906.2811
68 -58130.9141 -99661.3241
69 -136377.9104 -58130.9141
70 -28840.9276 -136377.9104
71 73016.2831 -28840.9276
72 -125491.3547 73016.2831
73 -14991.2288 -125491.3547
74 -121785.5624 -14991.2288
75 -141157.5865 -121785.5624
76 -92252.6956 -141157.5865
77 -47277.0402 -92252.6956
78 87516.0975 -47277.0402
79 -56250.2459 87516.0975
80 51772.1918 -56250.2459
81 -98359.6297 51772.1918
82 -118980.2572 -98359.6297
83 89337.8404 -118980.2572
84 -5467.0157 89337.8404
85 -111554.4237 -5467.0157
86 30263.8277 -111554.4237
87 -39334.6250 30263.8277
88 -95805.4969 -39334.6250
89 -20137.4147 -95805.4969
90 -26562.2862 -20137.4147
91 -102847.2821 -26562.2862
92 136445.8672 -102847.2821
93 -339389.1301 136445.8672
94 142088.0632 -339389.1301
95 -22029.2939 142088.0632
96 12387.1874 -22029.2939
97 47889.5549 12387.1874
98 -132866.7012 47889.5549
99 -85702.5943 -132866.7012
100 -93984.4068 -85702.5943
101 -101209.0829 -93984.4068
102 -110607.5576 -101209.0829
103 -53863.5304 -110607.5576
104 59661.3416 -53863.5304
105 -44169.0784 59661.3416
106 -239180.3303 -44169.0784
107 -1079.6495 -239180.3303
108 -203731.6259 -1079.6495
109 -152145.4789 -203731.6259
110 116244.2019 -152145.4789
111 -234855.3016 116244.2019
112 -268536.7701 -234855.3016
113 -174382.2503 -268536.7701
114 -7876.3314 -174382.2503
115 79594.8195 -7876.3314
116 11626.7819 79594.8195
117 40981.2854 11626.7819
118 -237174.4564 40981.2854
119 77984.4846 -237174.4564
120 -260529.3804 77984.4846
121 44126.2494 -260529.3804
122 2437.0379 44126.2494
123 -48248.7533 2437.0379
124 -120768.4011 -48248.7533
125 -249849.1739 -120768.4011
126 56743.8463 -249849.1739
127 57048.3791 56743.8463
128 64972.9222 57048.3791
129 56743.8463 64972.9222
130 59962.9990 56743.8463
131 114381.9051 59962.9990
132 56743.8463 114381.9051
133 54723.3658 56743.8463
134 56743.8463 54723.3658
135 59317.8032 56743.8463
136 59103.1536 59317.8032
137 56743.8463 59103.1536
138 56743.8463 56743.8463
139 56743.8463 56743.8463
140 55786.6489 56743.8463
141 56743.8463 55786.6489
142 51497.1868 56743.8463
143 57019.2654 51497.1868
144 63025.8180 57019.2654
145 68710.2635 63025.8180
146 56743.8463 68710.2635
147 58723.7837 56743.8463
148 56743.8463 58723.7837
149 56178.8235 56743.8463
150 47378.3959 56178.8235
151 46154.5855 47378.3959
152 56743.8463 46154.5855
153 34483.4455 56743.8463
154 56844.1304 34483.4455
155 56743.8463 56844.1304
156 56743.8463 56743.8463
157 53573.6402 56743.8463
158 63896.5616 53573.6402
159 56743.8463 63896.5616
160 56743.8463 56743.8463
161 56743.8463 56743.8463
162 56743.8463 56743.8463
163 62546.0270 56743.8463
164 57582.2679 62546.0270
165 56743.8463 57582.2679
166 56743.8463 56743.8463
167 160360.4989 56743.8463
168 56743.8463 160360.4989
169 56743.8463 56743.8463
170 56743.8463 56743.8463
171 53760.2551 56743.8463
172 56743.8463 53760.2551
173 54511.6315 56743.8463
174 56743.8463 54511.6315
175 56801.5725 56743.8463
176 56743.8463 56801.5725
177 -15353.1898 56743.8463
178 56743.8463 -15353.1898
179 56743.8463 56743.8463
180 70989.7457 56743.8463
181 -37474.7381 70989.7457
182 82509.7404 -37474.7381
183 -76005.3479 82509.7404
184 -14797.0131 -76005.3479
185 -21688.9005 -14797.0131
186 -8918.1814 -21688.9005
187 53423.5954 -8918.1814
188 17043.9417 53423.5954
189 56743.8463 17043.9417
190 61395.0315 56743.8463
191 56743.8463 61395.0315
192 58916.3004 56743.8463
193 56743.8463 58916.3004
194 56474.3785 56743.8463
195 -14102.2631 56474.3785
196 -150102.8416 -14102.2631
197 -14102.2631 -150102.8416
198 -14102.2631 -14102.2631
199 -92005.6434 -14102.2631
200 -14102.2631 -92005.6434
201 -14102.2631 -14102.2631
202 -14102.2631 -14102.2631
203 14151.1709 -14102.2631
204 -100326.1605 14151.1709
205 -15827.4682 -100326.1605
206 -29262.9780 -15827.4682
207 -14102.2631 -29262.9780
208 44819.1609 -14102.2631
209 -14102.2631 44819.1609
210 56743.8463 -14102.2631
211 56743.8463 56743.8463
212 72687.1102 56743.8463
213 34552.6110 72687.1102
214 37988.8278 34552.6110
215 56743.8463 37988.8278
216 60320.3399 56743.8463
217 61069.6158 60320.3399
218 56743.8463 61069.6158
219 61629.1210 56743.8463
220 23935.1655 61629.1210
221 56743.8463 23935.1655
222 57658.7888 56743.8463
223 56743.8463 57658.7888
224 57060.4504 56743.8463
225 61626.0995 57060.4504
226 56743.8463 61626.0995
227 56743.8463 56743.8463
228 30572.5273 56743.8463
229 56743.8463 30572.5273
230 35628.8129 56743.8463
231 56743.8463 35628.8129
232 56743.8463 56743.8463
233 58191.0580 56743.8463
234 56743.8463 58191.0580
235 57007.1910 56743.8463
236 56743.8463 57007.1910
237 59739.9250 56743.8463
238 56743.8463 59739.9250
239 56743.8463 56743.8463
240 56743.8463 56743.8463
241 59635.8630 56743.8463
242 56743.8463 59635.8630
243 61991.0854 56743.8463
244 -30514.6252 61991.0854
245 -14102.2631 -30514.6252
246 -14102.2631 -14102.2631
247 48670.3819 -14102.2631
248 -14102.2631 48670.3819
249 -14102.2631 -14102.2631
250 63819.4336 -14102.2631
251 -14102.2631 63819.4336
252 33344.1466 -14102.2631
253 37600.1858 33344.1466
254 6243.3979 37600.1858
255 56743.8463 6243.3979
256 80097.7803 56743.8463
257 96741.5415 80097.7803
258 -14102.2631 96741.5415
259 56743.8463 -14102.2631
260 -14102.2631 56743.8463
261 -14102.2631 -14102.2631
262 37163.9577 -14102.2631
263 65359.4770 37163.9577
264 -14102.2631 65359.4770
265 56743.8463 -14102.2631
266 -14102.2631 56743.8463
267 60145.7968 -14102.2631
268 -14102.2631 60145.7968
269 35248.8380 -14102.2631
270 -14102.2631 35248.8380
271 -14102.2631 -14102.2631
272 -14102.2631 -14102.2631
273 55325.6036 -14102.2631
274 56743.8463 55325.6036
275 66053.6292 56743.8463
276 63962.5695 66053.6292
277 56743.8463 63962.5695
278 -25065.7169 56743.8463
279 35917.4441 -25065.7169
280 -14102.2631 35917.4441
281 56743.8463 -14102.2631
282 95798.9808 56743.8463
283 56743.8463 95798.9808
284 20948.2632 56743.8463
285 56743.8463 20948.2632
286 -14102.2631 56743.8463
287 55796.2699 -14102.2631
288 -14102.2631 55796.2699
289 87737.0651 -14102.2631
290 -14102.2631 87737.0651
291 66657.4778 -14102.2631
292 30150.0010 66657.4778
293 56743.8463 30150.0010
294 60459.4379 56743.8463
295 20666.8015 60459.4379
296 56743.8463 20666.8015
297 -14102.2631 56743.8463
298 56743.8463 -14102.2631
299 56743.8463 56743.8463
300 -14799.5421 56743.8463
301 42657.9409 -14799.5421
302 56743.8463 42657.9409
303 -14102.2631 56743.8463
304 96052.7763 -14102.2631
305 60905.0179 96052.7763
306 -14102.2631 60905.0179
307 -22026.4417 -14102.2631
308 78605.8994 -22026.4417
309 -14102.2631 78605.8994
310 -14102.2631 -14102.2631
311 37401.7567 -14102.2631
312 -14102.2631 37401.7567
313 -17024.1271 -14102.2631
314 58714.9646 -17024.1271
315 57989.7061 58714.9646
316 58333.6778 57989.7061
317 42709.8089 58333.6778
318 59542.2511 42709.8089
319 56743.8463 59542.2511
320 56743.8463 56743.8463
321 56743.8463 56743.8463
322 56743.8463 56743.8463
323 -14102.2631 56743.8463
324 56743.8463 -14102.2631
325 54184.7061 56743.8463
326 -14102.2631 54184.7061
327 -14102.2631 -14102.2631
328 56743.8463 -14102.2631
329 56743.8463 56743.8463
330 56743.8463 56743.8463
331 56743.8463 56743.8463
332 56743.8463 56743.8463
333 56785.5558 56743.8463
334 57641.3511 56785.5558
335 56743.8463 57641.3511
336 55622.0937 56743.8463
337 56743.8463 55622.0937
338 56098.6316 56743.8463
339 60825.7496 56098.6316
340 37739.2850 60825.7496
341 56743.8463 37739.2850
342 56743.8463 56743.8463
343 56743.8463 56743.8463
344 56743.8463 56743.8463
345 56743.8463 56743.8463
346 48733.4452 56743.8463
347 56743.8463 48733.4452
348 56743.8463 56743.8463
349 56743.8463 56743.8463
350 57507.4643 56743.8463
351 58135.5877 57507.4643
352 58182.2306 58135.5877
353 46857.2332 58182.2306
354 56743.8463 46857.2332
355 56743.8463 56743.8463
356 24153.8623 56743.8463
357 56743.8463 24153.8623
358 56743.8463 56743.8463
359 56743.8463 56743.8463
360 58592.2977 56743.8463
361 56743.8463 58592.2977
362 56743.8463 56743.8463
363 -23229.0752 56743.8463
364 55339.9264 -23229.0752
365 56743.8463 55339.9264
366 56743.8463 56743.8463
367 56743.8463 56743.8463
368 56743.8463 56743.8463
369 56743.8463 56743.8463
370 56743.8463 56743.8463
371 59951.3846 56743.8463
372 16950.6168 59951.3846
373 29944.4766 16950.6168
374 56743.8463 29944.4766
375 56743.8463 56743.8463
376 69183.7439 56743.8463
377 41529.1872 69183.7439
378 56743.8463 41529.1872
379 56743.8463 56743.8463
380 56743.8463 56743.8463
381 50304.4797 56743.8463
382 110203.6465 50304.4797
383 56743.8463 110203.6465
384 58971.6081 56743.8463
385 -99802.4349 58971.6081
386 61423.9876 -99802.4349
387 40725.4677 61423.9876
388 54932.2314 40725.4677
389 37401.4713 54932.2314
390 53556.4003 37401.4713
391 36181.4677 53556.4003
392 18933.7581 36181.4677
393 73314.7260 18933.7581
394 -225558.1571 73314.7260
395 29067.6811 -225558.1571
396 -107954.0035 29067.6811
397 -1003.1517 -107954.0035
398 30167.9153 -1003.1517
399 -20626.9915 30167.9153
400 18801.9662 -20626.9915
401 46542.0025 18801.9662
402 35992.8637 46542.0025
403 -96274.3599 35992.8637
404 23096.4155 -96274.3599
405 4588.9031 23096.4155
406 41079.1019 4588.9031
407 36064.8764 41079.1019
408 34340.5635 36064.8764
409 -18627.8954 34340.5635
410 12187.3687 -18627.8954
411 -219006.8707 12187.3687
412 -43752.6695 -219006.8707
413 -154815.6019 -43752.6695
414 -54634.8218 -154815.6019
415 -37896.3342 -54634.8218
416 138483.9631 -37896.3342
417 -67743.4506 138483.9631
418 5703.9594 -67743.4506
419 -11870.2290 5703.9594
420 -65172.7527 -11870.2290
421 -99395.9490 -65172.7527
422 -427200.0501 -99395.9490
423 -470445.0446 -427200.0501
424 -71049.1911 -470445.0446
425 -200544.6914 -71049.1911
426 -306223.7944 -200544.6914
427 -1206200.7325 -306223.7944
428 -359249.8705 -1206200.7325
429 -270057.1479 -359249.8705
430 -843491.5827 -270057.1479
> 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/7gbfc1291400370.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/8gbfc1291400370.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/9gbfc1291400370.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/10r3ef1291400370.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/11dldk1291400370.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/12g4t81291400370.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/13ud9h1291400370.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/14fw8n1291400370.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/15jf6t1291400370.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/164fnz1291400370.tab")
+ }
>
> try(system("convert tmp/122zl1291400370.ps tmp/122zl1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/2vbzo1291400370.ps tmp/2vbzo1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/3vbzo1291400370.ps tmp/3vbzo1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/4vbzo1291400370.ps tmp/4vbzo1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/562gr1291400370.ps tmp/562gr1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/662gr1291400370.ps tmp/662gr1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/7gbfc1291400370.ps tmp/7gbfc1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/8gbfc1291400370.ps tmp/8gbfc1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/9gbfc1291400370.ps tmp/9gbfc1291400370.png",intern=TRUE))
character(0)
> try(system("convert tmp/10r3ef1291400370.ps tmp/10r3ef1291400370.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
14.453 3.229 15.209