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 = 'Linear Trend'
> par2 = 'Include Monthly 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\r M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 M11
1 6282929 1 162556 1081 213118 1 0 0 0 0 0 0 0 0 0 0
2 4324047 1 29790 309 81767 0 1 0 0 0 0 0 0 0 0 0
3 4108272 1 87550 458 153198 0 0 1 0 0 0 0 0 0 0 0
4 -1212617 0 84738 588 -26007 0 0 0 1 0 0 0 0 0 0 0
5 1485329 1 54660 299 126942 0 0 0 0 1 0 0 0 0 0 0
6 1779876 1 42634 156 157214 0 0 0 0 0 1 0 0 0 0 0
7 1367203 0 40949 481 129352 0 0 0 0 0 0 1 0 0 0 0
8 2519076 1 42312 323 234817 0 0 0 0 0 0 0 1 0 0 0
9 912684 1 37704 452 60448 0 0 0 0 0 0 0 0 1 0 0
10 1443586 1 16275 109 47818 0 0 0 0 0 0 0 0 0 1 0
11 1220017 0 25830 115 245546 0 0 0 0 0 0 0 0 0 0 1
12 984885 0 12679 110 48020 0 0 0 0 0 0 0 0 0 0 0
13 1457425 1 18014 239 -1710 1 0 0 0 0 0 0 0 0 0 0
14 -572920 0 43556 247 32648 0 1 0 0 0 0 0 0 0 0 0
15 929144 1 24524 497 95350 0 0 1 0 0 0 0 0 0 0 0
16 1151176 0 6532 103 151352 0 0 0 1 0 0 0 0 0 0 0
17 790090 0 7123 109 288170 0 0 0 0 1 0 0 0 0 0 0
18 774497 1 20813 502 114337 0 0 0 0 0 1 0 0 0 0 0
19 990576 1 37597 248 37884 0 0 0 0 0 0 1 0 0 0 0
20 454195 0 17821 373 122844 0 0 0 0 0 0 0 1 0 0 0
21 876607 1 12988 119 82340 0 0 0 0 0 0 0 0 1 0 0
22 711969 1 22330 84 79801 0 0 0 0 0 0 0 0 0 1 0
23 702380 0 13326 102 165548 0 0 0 0 0 0 0 0 0 0 1
24 264449 0 16189 295 116384 0 0 0 0 0 0 0 0 0 0 0
25 450033 0 7146 105 134028 1 0 0 0 0 0 0 0 0 0 0
26 541063 0 15824 64 63838 0 1 0 0 0 0 0 0 0 0 0
27 588864 1 26088 267 74996 0 0 1 0 0 0 0 0 0 0 0
28 -37216 0 11326 129 31080 0 0 0 1 0 0 0 0 0 0 0
29 783310 0 8568 37 32168 0 0 0 0 1 0 0 0 0 0 0
30 467359 0 14416 361 49857 0 0 0 0 0 1 0 0 0 0 0
31 688779 1 3369 28 87161 0 0 0 0 0 0 1 0 0 0 0
32 608419 1 11819 85 106113 0 0 0 0 0 0 0 1 0 0 0
33 696348 1 6620 44 80570 0 0 0 0 0 0 0 0 1 0 0
34 597793 1 4519 49 102129 0 0 0 0 0 0 0 0 0 1 0
35 821730 0 2220 22 301670 0 0 0 0 0 0 0 0 0 0 1
36 377934 0 18562 155 102313 0 0 0 0 0 0 0 0 0 0 0
37 651939 0 10327 91 88577 1 0 0 0 0 0 0 0 0 0 0
38 697458 1 5336 81 112477 0 1 0 0 0 0 0 0 0 0 0
39 700368 1 2365 79 191778 0 0 1 0 0 0 0 0 0 0 0
40 225986 0 4069 145 79804 0 0 0 1 0 0 0 0 0 0 0
41 348695 0 7710 816 128294 0 0 0 0 1 0 0 0 0 0 0
42 373683 0 13718 61 96448 0 0 0 0 0 1 0 0 0 0 0
43 501709 0 4525 226 93811 0 0 0 0 0 0 1 0 0 0 0
44 413743 0 6869 105 117520 0 0 0 0 0 0 0 1 0 0 0
45 379825 0 4628 62 69159 0 0 0 0 0 0 0 0 1 0 0
46 336260 1 3653 24 101792 0 0 0 0 0 0 0 0 0 1 0
47 636765 1 1265 26 210568 0 0 0 0 0 0 0 0 0 0 1
48 481231 1 7489 322 136996 0 0 0 0 0 0 0 0 0 0 0
49 469107 0 4901 84 121920 1 0 0 0 0 0 0 0 0 0 0
50 211928 0 2284 33 76403 0 1 0 0 0 0 0 0 0 0 0
51 563925 1 3160 108 108094 0 0 1 0 0 0 0 0 0 0 0
52 511939 1 4150 150 134759 0 0 0 1 0 0 0 0 0 0 0
53 521016 1 7285 115 188873 0 0 0 0 1 0 0 0 0 0 0
54 543856 1 1134 162 146216 0 0 0 0 0 1 0 0 0 0 0
55 329304 1 4658 158 156608 0 0 0 0 0 0 1 0 0 0 0
56 423262 0 2384 97 61348 0 0 0 0 0 0 0 1 0 0 0
57 509665 0 3748 9 50350 0 0 0 0 0 0 0 0 1 0 0
58 455881 0 5371 66 87720 0 0 0 0 0 0 0 0 0 1 0
59 367772 0 1285 107 99489 0 0 0 0 0 0 0 0 0 0 1
60 406339 1 9327 101 87419 0 0 0 0 0 0 0 0 0 0 0
61 493408 1 5565 47 94355 1 0 0 0 0 0 0 0 0 0 0
62 232942 0 1528 38 60326 0 1 0 0 0 0 0 0 0 0 0
63 416002 1 3122 34 94670 0 0 1 0 0 0 0 0 0 0 0
64 337430 1 7317 84 82425 0 0 0 1 0 0 0 0 0 0 0
65 361517 0 2675 79 59017 0 0 0 0 1 0 0 0 0 0 0
66 360962 0 13253 947 90829 0 0 0 0 0 1 0 0 0 0 0
67 235561 0 880 74 80791 0 0 0 0 0 0 1 0 0 0 0
68 408247 1 2053 53 100423 0 0 0 0 0 0 0 1 0 0 0
69 450296 0 1424 94 131116 0 0 0 0 0 0 0 0 1 0 0
70 418799 1 4036 63 100269 0 0 0 0 0 0 0 0 0 1 0
71 247405 1 3045 58 27330 0 0 0 0 0 0 0 0 0 0 1
72 378519 0 5119 49 39039 0 0 0 0 0 0 0 0 0 0 0
73 326638 0 1431 34 106885 1 0 0 0 0 0 0 0 0 0 0
74 328233 0 554 11 79285 0 1 0 0 0 0 0 0 0 0 0
75 386225 0 1975 35 118881 0 0 1 0 0 0 0 0 0 0 0
76 283662 1 1286 17 77623 0 0 0 1 0 0 0 0 0 0 0
77 370225 0 1012 47 114768 0 0 0 0 1 0 0 0 0 0 0
78 269236 0 810 43 74015 0 0 0 0 0 1 0 0 0 0 0
79 365732 0 1280 117 69465 0 0 0 0 0 0 1 0 0 0 0
80 420383 1 666 171 117869 0 0 0 0 0 0 0 1 0 0 0
81 345811 0 1380 26 60982 0 0 0 0 0 0 0 0 1 0 0
82 431809 1 4608 73 90131 0 0 0 0 0 0 0 0 0 1 0
83 418876 0 876 59 138971 0 0 0 0 0 0 0 0 0 0 1
84 297476 0 814 18 39625 0 0 0 0 0 0 0 0 0 0 0
85 416776 0 514 15 102725 1 0 0 0 0 0 0 0 0 0 0
86 357257 1 5692 72 64239 0 1 0 0 0 0 0 0 0 0 0
87 458343 0 3642 86 90262 0 0 1 0 0 0 0 0 0 0 0
88 388386 0 540 14 103960 0 0 0 1 0 0 0 0 0 0 0
89 358934 0 2099 64 106611 0 0 0 0 1 0 0 0 0 0 0
90 407560 0 567 11 103345 0 0 0 0 0 1 0 0 0 0 0
91 392558 0 2001 52 95551 0 0 0 0 0 0 1 0 0 0 0
92 373177 1 2949 41 82903 0 0 0 0 0 0 0 1 0 0 0
93 428370 0 2253 99 63593 0 0 0 0 0 0 0 0 1 0 0
94 369419 1 6533 75 126910 0 0 0 0 0 0 0 0 0 1 0
95 358649 0 1889 45 37527 0 0 0 0 0 0 0 0 0 0 1
96 376641 1 3055 43 60247 0 0 0 0 0 0 0 0 0 0 0
97 467427 0 272 8 112995 1 0 0 0 0 0 0 0 0 0 0
98 364885 1 1414 198 70184 0 1 0 0 0 0 0 0 0 0 0
99 436230 0 2564 22 130140 0 0 1 0 0 0 0 0 0 0 0
100 329118 1 1383 11 73221 0 0 0 1 0 0 0 0 0 0 0
101 317365 1 1261 33 76114 0 0 0 0 1 0 0 0 0 0 0
102 286849 0 975 23 90534 0 0 0 0 0 1 0 0 0 0 0
103 376685 0 3366 80 108479 0 0 0 0 0 0 1 0 0 0 0
104 407198 0 576 18 113761 0 0 0 0 0 0 0 1 0 0 0
105 377772 0 1306 28 68696 0 0 0 0 0 0 0 0 1 0 0
106 271483 0 746 23 71561 0 0 0 0 0 0 0 0 0 1 0
107 153661 1 3192 60 59831 0 0 0 0 0 0 0 0 0 0 1
108 513294 1 2045 20 97890 0 0 0 0 0 0 0 0 0 0 0
109 324881 0 5477 59 101481 1 0 0 0 0 0 0 0 0 0 0
110 264512 1 1932 36 72954 0 1 0 0 0 0 0 0 0 0 0
111 420968 0 936 30 67939 0 0 1 0 0 0 0 0 0 0 0
112 129302 1 3437 47 48022 0 0 0 1 0 0 0 0 0 0 0
113 191521 0 5131 71 86111 0 0 0 0 1 0 0 0 0 0 0
114 268673 1 2397 14 74020 0 0 0 0 0 1 0 0 0 0 0
115 353179 1 1389 9 57530 0 0 0 0 0 0 1 0 0 0 0
116 354624 0 1503 39 56364 0 0 0 0 0 0 0 1 0 0 0
117 363713 0 402 26 84990 0 0 0 0 0 0 0 0 1 0 0
118 456657 0 2239 21 88590 0 0 0 0 0 0 0 0 0 1 0
119 211742 1 2234 16 77200 0 0 0 0 0 0 0 0 0 0 1
120 338381 0 837 69 61262 0 0 0 0 0 0 0 0 0 0 0
121 418530 0 10579 92 110309 1 0 0 0 0 0 0 0 0 0 0
122 351483 0 875 14 67000 0 1 0 0 0 0 0 0 0 0 0
123 372928 0 1395 103 93099 0 0 1 0 0 0 0 0 0 0 0
124 485538 1 1659 29 107577 0 0 0 1 0 0 0 0 0 0 0
125 279268 1 2647 37 62920 0 0 0 0 1 0 0 0 0 0 0
126 219060 1 3294 23 75832 0 0 0 0 0 1 0 0 0 0 0
127 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
128 325314 0 94 7 60793 0 0 0 0 0 0 0 1 0 0 0
129 322046 0 422 28 57935 0 0 0 0 0 0 0 0 1 0 0
130 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
131 325599 0 34 8 60630 0 0 0 0 0 0 0 0 0 0 1
132 377028 0 1558 63 55637 0 0 0 0 0 0 0 0 0 0 0
133 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
134 323850 0 43 3 60887 0 1 0 0 0 0 0 0 0 0 0
135 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
136 331514 0 316 9 60505 0 0 0 1 0 0 0 0 0 0 0
137 325632 0 115 13 60945 0 0 0 0 1 0 0 0 0 0 0
138 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
139 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
140 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
141 322265 0 389 14 58990 0 0 0 0 0 0 0 0 1 0 0
142 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
143 325906 0 1002 15 56750 0 0 0 0 0 0 0 0 0 0 1
144 325985 0 36 3 60894 0 0 0 0 0 0 0 0 0 0 0
145 346145 0 460 15 63346 1 0 0 0 0 0 0 0 0 0 0
146 325898 0 309 11 56535 0 1 0 0 0 0 0 0 0 0 0
147 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
148 325356 0 9 6 60835 0 0 0 1 0 0 0 0 0 0 0
149 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0
150 325930 0 14 1 61016 0 0 0 0 0 1 0 0 0 0 0
151 318020 0 520 10 58650 0 0 0 0 0 0 1 0 0 0 0
152 326389 0 1766 73 60438 0 0 0 0 0 0 0 1 0 0 0
153 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
154 302925 0 458 11 58625 0 0 0 0 0 0 0 0 0 1 0
155 325540 0 20 3 60938 0 0 0 0 0 0 0 0 0 0 1
156 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
157 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
158 326736 0 98 2 61490 0 1 0 0 0 0 0 0 0 0 0
159 340580 0 405 7 60845 0 0 1 0 0 0 0 0 0 0 0
160 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
161 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0
162 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
163 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
164 331828 0 483 27 60830 0 0 0 0 0 0 0 1 0 0 0
165 323299 0 454 51 63261 0 0 0 0 0 0 0 0 1 0 0
166 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
167 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
168 387722 0 757 19 45689 0 0 0 0 0 0 0 0 0 0 0
169 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
170 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
171 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
172 324598 0 36 4 61564 0 0 0 1 0 0 0 0 0 0 0
173 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0
174 328726 0 203 9 61938 0 0 0 0 0 1 0 0 0 0 0
175 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
176 325043 0 90 8 60951 0 0 0 0 0 0 0 1 0 0 0
177 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
178 325806 1 71 1 60745 0 0 0 0 0 0 0 0 0 1 0
179 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
180 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
181 387732 0 972 34 71642 1 0 0 0 0 0 0 0 0 0 0
182 349729 1 531 10 71641 0 1 0 0 0 0 0 0 0 0 0
183 332202 0 604 38 55792 0 0 1 0 0 0 0 0 0 0 0
184 305442 1 283 10 71873 0 0 0 1 0 0 0 0 0 0 0
185 329537 1 23 5 62555 0 0 0 0 1 0 0 0 0 0 0
186 327055 1 638 14 60370 0 0 0 0 0 1 0 0 0 0 0
187 356245 1 699 16 64873 0 0 0 0 0 0 1 0 0 0 0
188 328451 0 149 5 62041 0 0 0 0 0 0 0 1 0 0 0
189 307062 0 226 5 65745 0 0 0 0 0 0 0 0 1 0 0
190 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
191 331345 0 275 4 59500 0 0 0 0 0 0 0 0 0 0 1
192 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
193 331824 0 141 6 61630 1 0 0 0 0 0 0 0 0 0 0
194 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
195 325685 0 28 2 60890 0 0 1 0 0 0 0 0 0 0 0
196 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
197 404480 1 2566 80 113521 0 0 0 0 1 0 0 0 0 0 0
198 325560 1 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
199 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
200 318314 1 472 20 80045 0 0 0 0 0 0 0 1 0 0 0
201 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
202 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
203 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
204 311807 1 203 27 50804 0 0 0 0 0 0 0 0 0 0 0
205 337724 1 496 17 87390 1 0 0 0 0 0 0 0 0 0 0
206 326431 1 10 2 61656 0 1 0 0 0 0 0 0 0 0 0
207 327556 1 63 4 65688 0 0 1 0 0 0 0 0 0 0 0
208 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
209 356850 1 1136 32 48522 0 0 0 0 1 0 0 0 0 0 0
210 325560 1 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
211 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
212 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
213 322741 0 267 32 57640 0 0 0 0 0 0 0 0 1 0 0
214 310902 0 474 20 61977 0 0 0 0 0 0 0 0 0 1 0
215 324295 0 534 7 62620 0 0 0 0 0 0 0 0 0 0 1
216 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
217 326156 0 15 8 60831 1 0 0 0 0 0 0 0 0 0 0
218 326960 0 397 28 60646 0 1 0 0 0 0 0 0 0 0 0
219 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
220 333411 0 1061 20 56225 0 0 0 1 0 0 0 0 0 0 0
221 297761 0 288 4 60510 0 0 0 0 1 0 0 0 0 0 0
222 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
223 325536 0 3 2 60698 0 0 0 0 0 0 1 0 0 0 0
224 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
225 325762 0 20 2 60805 0 0 0 0 0 0 0 0 1 0 0
226 327957 0 278 26 61404 0 0 0 0 0 0 0 0 0 1 0
227 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
228 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
229 318521 0 192 4 65276 1 0 0 0 0 0 0 0 0 0 0
230 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
231 319775 0 317 9 63915 0 0 1 0 0 0 0 0 0 0 0
232 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
233 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0
234 332128 0 368 17 61686 0 0 0 0 0 1 0 0 0 0 0
235 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
236 325486 0 2 1 60743 0 0 0 0 0 0 0 1 0 0 0
237 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
238 325838 0 53 6 60349 0 0 0 0 0 0 0 0 0 1 0
239 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
240 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
241 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
242 331767 0 94 3 61360 0 1 0 0 0 0 0 0 0 0 0
243 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
244 324523 0 24 8 59818 0 0 0 1 0 0 0 0 0 0 0
245 339995 0 2332 4 72680 0 0 0 0 1 0 0 0 0 0 0
246 325560 1 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
247 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
248 319582 0 131 11 61808 0 0 0 0 0 0 0 1 0 0 0
249 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
250 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
251 307245 0 206 9 53110 0 0 0 0 0 0 0 0 0 0 1
252 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
253 317967 0 167 2 64245 1 0 0 0 0 0 0 0 0 0 0
254 331488 0 622 73 73007 0 1 0 0 0 0 0 0 0 0 0
255 335452 0 885 85 82732 0 0 1 0 0 0 0 0 0 0 0
256 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
257 334184 0 365 8 54820 0 0 0 0 1 0 0 0 0 0 0
258 313213 0 364 35 47705 0 0 0 0 0 1 0 0 0 0 0
259 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
260 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
261 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
262 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
263 348678 0 226 12 72835 0 0 0 0 0 0 0 0 0 0 1
264 328727 0 307 15 58856 0 0 0 0 0 0 0 0 0 0 0
265 325560 1 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
266 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
267 325560 1 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
268 387978 0 188 11 77655 0 0 0 1 0 0 0 0 0 0 0
269 325560 1 0 0 60720 0 0 0 0 1 0 0 0 0 0 0
270 336704 0 138 6 69817 0 0 0 0 0 1 0 0 0 0 0
271 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
272 325560 1 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
273 325560 1 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
274 322076 0 125 12 60798 0 0 0 0 0 0 0 0 0 1 0
275 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
276 334272 0 282 30 62452 0 0 0 0 0 0 0 0 0 0 0
277 338197 0 335 33 64175 1 0 0 0 0 0 0 0 0 0 0
278 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
279 321024 1 813 82 67440 0 0 1 0 0 0 0 0 0 0 0
280 322145 0 176 28 68136 0 0 0 1 0 0 0 0 0 0 0
281 325560 1 0 0 60720 0 0 0 0 1 0 0 0 0 0 0
282 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
283 323351 0 249 72 56726 0 0 0 0 0 0 1 0 0 0 0
284 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
285 327748 0 333 13 70811 0 0 0 0 0 0 0 0 1 0 0
286 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
287 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
288 328157 0 30 4 62045 0 0 0 0 0 0 0 0 0 0 0
289 325560 1 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
290 311594 0 249 62 54323 0 1 0 0 0 0 0 0 0 0 0
291 325560 1 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
292 335962 0 165 24 62841 0 0 0 1 0 0 0 0 0 0 0
293 372426 0 453 21 81125 0 0 0 0 1 0 0 0 0 0 0
294 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
295 319844 0 53 14 59506 0 0 0 0 0 0 1 0 0 0 0
296 355822 1 382 21 59365 0 0 0 0 0 0 0 1 0 0 0
297 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
298 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
299 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
300 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
301 324047 1 30 4 60798 1 0 0 0 0 0 0 0 0 0 0
302 311464 0 290 2 58790 0 1 0 0 0 0 0 0 0 0 0
303 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
304 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
305 353417 0 366 53 61808 0 0 0 0 1 0 0 0 0 0 0
306 325590 0 2 9 60735 0 0 0 0 0 1 0 0 0 0 0
307 325560 1 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
308 328576 1 209 13 64016 0 0 0 0 0 0 0 1 0 0 0
309 326126 0 384 22 54683 0 0 0 0 0 0 0 0 1 0 0
310 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
311 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
312 369376 0 365 83 87192 0 0 0 0 0 0 0 0 0 0 0
313 325560 1 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
314 332013 1 49 8 64107 0 1 0 0 0 0 0 0 0 0 0
315 325871 0 3 4 60761 0 0 1 0 0 0 0 0 0 0 0
316 342165 0 133 14 65990 0 0 0 1 0 0 0 0 0 0 0
317 324967 0 32 1 59988 0 0 0 0 1 0 0 0 0 0 0
318 314832 0 368 17 61167 0 0 0 0 0 1 0 0 0 0 0
319 325557 0 1 6 60719 0 0 0 0 0 0 1 0 0 0 0
320 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
321 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
322 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
323 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
324 325560 1 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
325 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
326 322649 0 22 2 60722 0 1 0 0 0 0 0 0 0 0 0
327 325560 1 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
328 325560 1 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
329 325560 0 0 0 60720 0 0 0 0 1 0 0 0 0 0 0
330 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
331 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
332 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
333 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
334 324598 0 96 5 60379 0 0 0 0 0 0 0 0 0 1 0
335 325567 0 1 2 60727 0 0 0 0 0 0 0 0 0 0 1
336 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
337 324005 0 81 7 60925 1 0 0 0 0 0 0 0 0 0 0
338 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
339 325748 0 26 1 60896 0 0 1 0 0 0 0 0 0 0 0
340 323385 0 125 13 59734 0 0 0 1 0 0 0 0 0 0 0
341 315409 0 304 15 62969 0 0 0 0 1 0 0 0 0 0 0
342 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
343 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
344 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
345 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
346 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
347 312275 0 119 6 59118 0 0 0 0 0 0 0 0 0 0 1
348 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
349 325560 0 0 0 60720 1 0 0 0 0 0 0 0 0 0 0
350 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
351 320576 0 312 14 58598 0 0 1 0 0 0 0 0 0 0 0
352 325246 0 60 10 61124 0 0 0 1 0 0 0 0 0 0 0
353 332961 0 587 12 59595 0 0 0 0 1 0 0 0 0 0 0
354 323010 0 135 2 62065 0 0 0 0 0 1 0 0 0 0 0
355 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
356 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
357 345253 0 514 52 78780 0 0 0 0 0 0 0 0 1 0 0
358 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
359 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
360 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
361 325559 0 1 4 60722 1 0 0 0 0 0 0 0 0 0 0
362 325560 0 0 0 60720 0 1 0 0 0 0 0 0 0 0 0
363 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
364 319634 1 58 3 61600 0 0 0 1 0 0 0 0 0 0 0
365 319951 0 180 11 59635 0 0 0 0 1 0 0 0 0 0 0
366 325560 0 0 0 60720 0 0 0 0 0 1 0 0 0 0 0
367 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
368 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
369 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
370 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 1 0
371 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 1
372 318519 0 448 40 59781 0 0 0 0 0 0 0 0 0 0 0
373 343222 0 227 9 76644 1 0 0 0 0 0 0 0 0 0 0
374 317234 0 174 1 64820 0 1 0 0 0 0 0 0 0 0 0
375 325560 0 0 0 60720 0 0 1 0 0 0 0 0 0 0 0
376 325560 0 0 0 60720 0 0 0 1 0 0 0 0 0 0 0
377 314025 0 121 24 56178 0 0 0 0 1 0 0 0 0 0 0
378 320249 0 607 11 60436 0 0 0 0 0 1 0 0 0 0 0
379 325560 0 0 0 60720 0 0 0 0 0 0 1 0 0 0 0
380 325560 0 0 0 60720 0 0 0 0 0 0 0 1 0 0 0
381 325560 0 0 0 60720 0 0 0 0 0 0 0 0 1 0 0
382 349365 0 530 60 73433 0 0 0 0 0 0 0 0 0 1 0
383 289197 0 571 80 41477 0 0 0 0 0 0 0 0 0 0 1
384 325560 0 0 0 60720 0 0 0 0 0 0 0 0 0 0 0
385 329245 0 78 16 62700 1 0 0 0 0 0 0 0 0 0 0
386 240869 0 2489 40 67804 0 1 0 0 0 0 0 0 0 0 0
387 327182 0 131 6 59661 0 0 1 0 0 0 0 0 0 0 0
388 322876 0 923 8 58620 0 0 0 1 0 0 0 0 0 0 0
389 323117 0 72 3 60398 0 0 0 0 1 0 0 0 0 0 0
390 306351 0 572 16 58580 0 0 0 0 0 1 0 0 0 0 0
391 335137 0 397 10 62710 0 0 0 0 0 0 1 0 0 0 0
392 308271 0 450 8 59325 0 0 0 0 0 0 0 1 0 0 0
393 301731 0 622 7 60950 0 0 0 0 0 0 0 0 1 0 0
394 382409 0 694 8 68060 0 0 0 0 0 0 0 0 0 1 0
395 279230 1 3425 12 83620 0 0 0 0 0 0 0 0 0 0 1
396 298731 0 562 13 58456 0 0 0 0 0 0 0 0 0 0 0
397 243650 0 4917 42 52811 1 0 0 0 0 0 0 0 0 0 0
398 532682 1 1442 118 121173 0 1 0 0 0 0 0 0 0 0 0
399 319771 0 529 9 63870 0 0 1 0 0 0 0 0 0 0 0
400 171493 1 2126 138 21001 0 0 0 1 0 0 0 0 0 0 0
401 347262 0 1061 5 70415 0 0 0 0 1 0 0 0 0 0 0
402 343945 0 776 9 64230 0 0 0 0 0 1 0 0 0 0 0
403 311874 0 611 8 59190 0 0 0 0 0 0 1 0 0 0 0
404 302211 1 1526 25 69351 0 0 0 0 0 0 0 1 0 0 0
405 316708 0 592 7 64270 0 0 0 0 0 0 0 0 1 0 0
406 333463 0 1182 13 70694 0 0 0 0 0 0 0 0 0 1 0
407 344282 0 621 16 68005 0 0 0 0 0 0 0 0 0 0 1
408 319635 0 989 11 58930 0 0 0 0 0 0 0 0 0 0 0
409 301186 0 438 11 58320 1 0 0 0 0 0 0 0 0 0 0
410 300381 0 726 3 69980 0 1 0 0 0 0 0 0 0 0 0
411 318765 0 1303 61 69863 0 0 1 0 0 0 0 0 0 0 0
412 286146 1 6341 24 63255 0 0 0 1 0 0 0 0 0 0 0
413 306844 1 1164 17 57320 0 0 0 0 1 0 0 0 0 0 0
414 307705 1 3310 33 75230 0 0 0 0 0 1 0 0 0 0 0
415 312448 0 1366 7 79420 0 0 0 0 0 0 1 0 0 0 0
416 299715 0 965 3 73490 0 0 0 0 0 0 0 1 0 0 0
417 373399 0 3256 66 35250 0 0 0 0 0 0 0 0 1 0 0
418 299446 1 1135 17 62285 0 0 0 0 0 0 0 0 0 1 0
419 325586 0 1270 26 69206 0 0 0 0 0 0 0 0 0 0 1
420 291221 0 661 3 65920 0 0 0 0 0 0 0 0 0 0 0
421 261173 0 1013 2 69770 1 0 0 0 0 0 0 0 0 0 0
422 255027 0 2844 67 72683 0 1 0 0 0 0 0 0 0 0 0
423 -78375 1 11528 70 -14545 0 0 1 0 0 0 0 0 0 0 0
424 -58143 0 6526 26 55830 0 0 0 1 0 0 0 0 0 0 0
425 227033 0 2264 24 55174 0 0 0 0 1 0 0 0 0 0 0
426 235098 1 4461 94 67038 0 0 0 0 0 1 0 0 0 0 0
427 21267 0 3999 30 51252 0 0 0 0 0 0 1 0 0 0 0
428 238675 0 35624 223 157278 0 0 0 0 0 0 0 1 0 0 0
429 197687 0 9252 48 79510 0 0 0 0 0 0 0 0 1 0 0
430 418341 0 15236 90 77440 0 0 0 0 0 0 0 0 0 1 0
431 -297706 0 18073 180 27284 0 0 0 0 0 0 0 0 0 0 1
t
1 1
2 2
3 3
4 4
5 5
6 6
7 7
8 8
9 9
10 10
11 11
12 12
13 13
14 14
15 15
16 16
17 17
18 18
19 19
20 20
21 21
22 22
23 23
24 24
25 25
26 26
27 27
28 28
29 29
30 30
31 31
32 32
33 33
34 34
35 35
36 36
37 37
38 38
39 39
40 40
41 41
42 42
43 43
44 44
45 45
46 46
47 47
48 48
49 49
50 50
51 51
52 52
53 53
54 54
55 55
56 56
57 57
58 58
59 59
60 60
61 61
62 62
63 63
64 64
65 65
66 66
67 67
68 68
69 69
70 70
71 71
72 72
73 73
74 74
75 75
76 76
77 77
78 78
79 79
80 80
81 81
82 82
83 83
84 84
85 85
86 86
87 87
88 88
89 89
90 90
91 91
92 92
93 93
94 94
95 95
96 96
97 97
98 98
99 99
100 100
101 101
102 102
103 103
104 104
105 105
106 106
107 107
108 108
109 109
110 110
111 111
112 112
113 113
114 114
115 115
116 116
117 117
118 118
119 119
120 120
121 121
122 122
123 123
124 124
125 125
126 126
127 127
128 128
129 129
130 130
131 131
132 132
133 133
134 134
135 135
136 136
137 137
138 138
139 139
140 140
141 141
142 142
143 143
144 144
145 145
146 146
147 147
148 148
149 149
150 150
151 151
152 152
153 153
154 154
155 155
156 156
157 157
158 158
159 159
160 160
161 161
162 162
163 163
164 164
165 165
166 166
167 167
168 168
169 169
170 170
171 171
172 172
173 173
174 174
175 175
176 176
177 177
178 178
179 179
180 180
181 181
182 182
183 183
184 184
185 185
186 186
187 187
188 188
189 189
190 190
191 191
192 192
193 193
194 194
195 195
196 196
197 197
198 198
199 199
200 200
201 201
202 202
203 203
204 204
205 205
206 206
207 207
208 208
209 209
210 210
211 211
212 212
213 213
214 214
215 215
216 216
217 217
218 218
219 219
220 220
221 221
222 222
223 223
224 224
225 225
226 226
227 227
228 228
229 229
230 230
231 231
232 232
233 233
234 234
235 235
236 236
237 237
238 238
239 239
240 240
241 241
242 242
243 243
244 244
245 245
246 246
247 247
248 248
249 249
250 250
251 251
252 252
253 253
254 254
255 255
256 256
257 257
258 258
259 259
260 260
261 261
262 262
263 263
264 264
265 265
266 266
267 267
268 268
269 269
270 270
271 271
272 272
273 273
274 274
275 275
276 276
277 277
278 278
279 279
280 280
281 281
282 282
283 283
284 284
285 285
286 286
287 287
288 288
289 289
290 290
291 291
292 292
293 293
294 294
295 295
296 296
297 297
298 298
299 299
300 300
301 301
302 302
303 303
304 304
305 305
306 306
307 307
308 308
309 309
310 310
311 311
312 312
313 313
314 314
315 315
316 316
317 317
318 318
319 319
320 320
321 321
322 322
323 323
324 324
325 325
326 326
327 327
328 328
329 329
330 330
331 331
332 332
333 333
334 334
335 335
336 336
337 337
338 338
339 339
340 340
341 341
342 342
343 343
344 344
345 345
346 346
347 347
348 348
349 349
350 350
351 351
352 352
353 353
354 354
355 355
356 356
357 357
358 358
359 359
360 360
361 361
362 362
363 363
364 364
365 365
366 366
367 367
368 368
369 369
370 370
371 371
372 372
373 373
374 374
375 375
376 376
377 377
378 378
379 379
380 380
381 381
382 382
383 383
384 384
385 385
386 386
387 387
388 388
389 389
390 390
391 391
392 392
393 393
394 394
395 395
396 396
397 397
398 398
399 399
400 400
401 401
402 402
403 403
404 404
405 405
406 406
407 407
408 408
409 409
410 410
411 411
412 412
413 413
414 414
415 415
416 416
417 417
418 418
419 419
420 420
421 421
422 422
423 423
424 424
425 425
426 426
427 427
428 428
429 429
430 430
431 431
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades `Dividends\r`
6.172e+04 7.975e+04 2.628e+01 -4.552e+02 3.552e+00
M1 M2 M3 M4 M5
4.485e+04 3.030e+04 1.306e+04 -1.351e+05 -4.050e+04
M6 M7 M8 M9 M10
-2.787e+04 -1.370e+04 -3.544e+04 1.512e+04 -1.340e+04
M11 t
-6.406e+04 3.328e+01
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-3006346 -43233 23134 67114 3219486
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.172e+04 7.326e+04 0.842 0.4000
Group 7.975e+04 3.359e+04 2.374 0.0180 *
Costs 2.628e+01 1.935e+00 13.584 < 2e-16 ***
Trades -4.552e+02 2.138e+02 -2.130 0.0338 *
`Dividends\r` 3.552e+00 5.183e-01 6.853 2.63e-11 ***
M1 4.485e+04 7.163e+04 0.626 0.5315
M2 3.030e+04 7.126e+04 0.425 0.6709
M3 1.306e+04 7.150e+04 0.183 0.8552
M4 -1.351e+05 7.152e+04 -1.889 0.0596 .
M5 -4.050e+04 7.157e+04 -0.566 0.5718
M6 -2.787e+04 7.150e+04 -0.390 0.6969
M7 -1.370e+04 7.129e+04 -0.192 0.8477
M8 -3.544e+04 7.156e+04 -0.495 0.6207
M9 1.512e+04 7.123e+04 0.212 0.8320
M10 -1.340e+04 7.166e+04 -0.187 0.8517
M11 -6.406e+04 7.179e+04 -0.892 0.3728
t 3.328e+01 1.360e+02 0.245 0.8068
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 3e+05 on 414 degrees of freedom
Multiple R-squared: 0.5671, Adjusted R-squared: 0.5504
F-statistic: 33.9 on 16 and 414 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 1.265917e-98 6.329587e-99
[2,] 1.0000000 6.716163e-105 3.358081e-105
[3,] 1.0000000 1.144361e-107 5.721805e-108
[4,] 1.0000000 1.863224e-118 9.316121e-119
[5,] 1.0000000 2.489475e-123 1.244738e-123
[6,] 1.0000000 1.190507e-123 5.952536e-124
[7,] 1.0000000 2.484028e-127 1.242014e-127
[8,] 1.0000000 2.545736e-129 1.272868e-129
[9,] 1.0000000 1.468990e-143 7.344950e-144
[10,] 1.0000000 1.841361e-195 9.206804e-196
[11,] 1.0000000 3.526058e-207 1.763029e-207
[12,] 1.0000000 4.472957e-211 2.236478e-211
[13,] 1.0000000 4.146188e-213 2.073094e-213
[14,] 1.0000000 6.641307e-219 3.320654e-219
[15,] 1.0000000 5.880558e-220 2.940279e-220
[16,] 1.0000000 6.401568e-219 3.200784e-219
[17,] 1.0000000 5.170429e-218 2.585215e-218
[18,] 1.0000000 2.280225e-224 1.140112e-224
[19,] 1.0000000 3.944189e-232 1.972095e-232
[20,] 1.0000000 2.341792e-232 1.170896e-232
[21,] 1.0000000 8.911331e-236 4.455665e-236
[22,] 1.0000000 8.690849e-236 4.345424e-236
[23,] 1.0000000 2.876262e-237 1.438131e-237
[24,] 1.0000000 7.101429e-239 3.550714e-239
[25,] 1.0000000 7.666977e-239 3.833489e-239
[26,] 1.0000000 7.503909e-240 3.751955e-240
[27,] 1.0000000 3.560756e-239 1.780378e-239
[28,] 1.0000000 4.664942e-238 2.332471e-238
[29,] 1.0000000 1.684288e-237 8.421442e-238
[30,] 1.0000000 1.074707e-236 5.373536e-237
[31,] 1.0000000 5.550422e-237 2.775211e-237
[32,] 1.0000000 1.048366e-237 5.241829e-238
[33,] 1.0000000 5.302879e-237 2.651439e-237
[34,] 1.0000000 3.360438e-236 1.680219e-236
[35,] 1.0000000 3.685753e-235 1.842876e-235
[36,] 1.0000000 8.287888e-237 4.143944e-237
[37,] 1.0000000 4.713826e-239 2.356913e-239
[38,] 1.0000000 1.069625e-244 5.348126e-245
[39,] 1.0000000 2.751407e-245 1.375704e-245
[40,] 1.0000000 5.928617e-246 2.964309e-246
[41,] 1.0000000 1.506898e-245 7.534488e-246
[42,] 1.0000000 4.158058e-247 2.079029e-247
[43,] 1.0000000 1.248983e-246 6.244916e-247
[44,] 1.0000000 1.312187e-245 6.560934e-246
[45,] 1.0000000 3.785175e-245 1.892588e-245
[46,] 1.0000000 7.019365e-246 3.509683e-246
[47,] 1.0000000 3.332265e-245 1.666133e-245
[48,] 1.0000000 4.510172e-246 2.255086e-246
[49,] 1.0000000 7.192212e-245 3.596106e-245
[50,] 1.0000000 7.868185e-244 3.934093e-244
[51,] 1.0000000 9.241163e-243 4.620581e-243
[52,] 1.0000000 4.299666e-242 2.149833e-242
[53,] 1.0000000 6.499896e-244 3.249948e-244
[54,] 1.0000000 8.375360e-244 4.187680e-244
[55,] 1.0000000 8.357762e-243 4.178881e-243
[56,] 1.0000000 3.425931e-242 1.712966e-242
[57,] 1.0000000 8.573595e-242 4.286797e-242
[58,] 1.0000000 1.545703e-241 7.728514e-242
[59,] 1.0000000 2.139412e-241 1.069706e-241
[60,] 1.0000000 1.080832e-240 5.404159e-241
[61,] 1.0000000 7.198123e-240 3.599062e-240
[62,] 1.0000000 6.576394e-239 3.288197e-239
[63,] 1.0000000 2.397219e-238 1.198609e-238
[64,] 1.0000000 3.027891e-238 1.513946e-238
[65,] 1.0000000 1.472809e-237 7.364046e-238
[66,] 1.0000000 1.885607e-236 9.428034e-237
[67,] 1.0000000 1.243230e-236 6.216152e-237
[68,] 1.0000000 1.318039e-236 6.590196e-237
[69,] 1.0000000 1.472417e-236 7.362084e-237
[70,] 1.0000000 5.630436e-236 2.815218e-236
[71,] 1.0000000 5.404747e-235 2.702374e-235
[72,] 1.0000000 6.955409e-234 3.477705e-234
[73,] 1.0000000 8.688381e-233 4.344191e-233
[74,] 1.0000000 3.588058e-233 1.794029e-233
[75,] 1.0000000 8.313263e-233 4.156631e-233
[76,] 1.0000000 1.865136e-234 9.325679e-235
[77,] 1.0000000 4.431713e-234 2.215856e-234
[78,] 1.0000000 5.904564e-233 2.952282e-233
[79,] 1.0000000 3.367609e-232 1.683804e-232
[80,] 1.0000000 1.525834e-231 7.629169e-232
[81,] 1.0000000 1.390739e-230 6.953696e-231
[82,] 1.0000000 1.479621e-229 7.398105e-230
[83,] 1.0000000 5.535371e-230 2.767686e-230
[84,] 1.0000000 4.415769e-229 2.207885e-229
[85,] 1.0000000 6.542160e-229 3.271080e-229
[86,] 1.0000000 5.752034e-228 2.876017e-228
[87,] 1.0000000 4.570108e-228 2.285054e-228
[88,] 1.0000000 1.963915e-228 9.819574e-229
[89,] 1.0000000 2.827842e-228 1.413921e-228
[90,] 1.0000000 8.968893e-228 4.484446e-228
[91,] 1.0000000 1.758811e-227 8.794056e-228
[92,] 1.0000000 2.626792e-227 1.313396e-227
[93,] 1.0000000 1.167800e-227 5.839002e-228
[94,] 1.0000000 8.747547e-230 4.373774e-230
[95,] 1.0000000 2.511665e-229 1.255833e-229
[96,] 1.0000000 1.247349e-228 6.236743e-229
[97,] 1.0000000 4.887978e-228 2.443989e-228
[98,] 1.0000000 3.416018e-227 1.708009e-227
[99,] 1.0000000 1.630125e-226 8.150627e-227
[100,] 1.0000000 6.630893e-228 3.315446e-228
[101,] 1.0000000 6.825004e-227 3.412502e-227
[102,] 1.0000000 9.660145e-227 4.830072e-227
[103,] 1.0000000 1.266026e-225 6.330131e-226
[104,] 1.0000000 9.071619e-225 4.535809e-225
[105,] 1.0000000 6.226392e-224 3.113196e-224
[106,] 1.0000000 6.947406e-223 3.473703e-223
[107,] 1.0000000 6.950678e-224 3.475339e-224
[108,] 1.0000000 7.778032e-223 3.889016e-223
[109,] 1.0000000 8.337801e-222 4.168901e-222
[110,] 1.0000000 1.128702e-220 5.643509e-221
[111,] 1.0000000 1.063518e-219 5.317590e-220
[112,] 1.0000000 8.727400e-219 4.363700e-219
[113,] 1.0000000 1.388885e-218 6.944425e-219
[114,] 1.0000000 1.917029e-217 9.585146e-218
[115,] 1.0000000 2.608724e-216 1.304362e-216
[116,] 1.0000000 3.253573e-215 1.626787e-215
[117,] 1.0000000 1.357398e-214 6.786992e-215
[118,] 1.0000000 1.304266e-213 6.521331e-214
[119,] 1.0000000 1.364725e-212 6.823624e-213
[120,] 1.0000000 1.701567e-211 8.507833e-212
[121,] 1.0000000 2.024847e-210 1.012423e-210
[122,] 1.0000000 2.684764e-209 1.342382e-209
[123,] 1.0000000 2.813268e-208 1.406634e-208
[124,] 1.0000000 2.236971e-207 1.118485e-207
[125,] 1.0000000 2.587885e-206 1.293943e-206
[126,] 1.0000000 3.126583e-205 1.563292e-205
[127,] 1.0000000 3.767932e-204 1.883966e-204
[128,] 1.0000000 4.697213e-203 2.348606e-203
[129,] 1.0000000 2.956765e-202 1.478383e-202
[130,] 1.0000000 3.093843e-201 1.546922e-201
[131,] 1.0000000 3.444551e-200 1.722276e-200
[132,] 1.0000000 4.272231e-199 2.136115e-199
[133,] 1.0000000 4.782359e-198 2.391179e-198
[134,] 1.0000000 5.713611e-197 2.856805e-197
[135,] 1.0000000 5.125367e-196 2.562683e-196
[136,] 1.0000000 5.585457e-195 2.792728e-195
[137,] 1.0000000 6.308486e-194 3.154243e-194
[138,] 1.0000000 7.534196e-193 3.767098e-193
[139,] 1.0000000 9.064075e-192 4.532037e-192
[140,] 1.0000000 1.056978e-190 5.284891e-191
[141,] 1.0000000 8.016863e-190 4.008432e-190
[142,] 1.0000000 8.487974e-189 4.243987e-189
[143,] 1.0000000 9.398941e-188 4.699470e-188
[144,] 1.0000000 1.118970e-186 5.594848e-187
[145,] 1.0000000 1.318741e-185 6.593706e-186
[146,] 1.0000000 1.430832e-184 7.154159e-185
[147,] 1.0000000 1.415082e-183 7.075411e-184
[148,] 1.0000000 1.544485e-182 7.722424e-183
[149,] 1.0000000 3.397231e-183 1.698616e-183
[150,] 1.0000000 3.828107e-182 1.914054e-182
[151,] 1.0000000 4.358312e-181 2.179156e-181
[152,] 1.0000000 5.003801e-180 2.501900e-180
[153,] 1.0000000 4.286715e-179 2.143357e-179
[154,] 1.0000000 4.563312e-178 2.281656e-178
[155,] 1.0000000 5.089968e-177 2.544984e-177
[156,] 1.0000000 5.833133e-176 2.916567e-176
[157,] 1.0000000 6.527925e-175 3.263962e-175
[158,] 1.0000000 6.901340e-174 3.450670e-174
[159,] 1.0000000 6.141913e-173 3.070957e-173
[160,] 1.0000000 6.614356e-172 3.307178e-172
[161,] 1.0000000 7.222604e-171 3.611302e-171
[162,] 1.0000000 4.558470e-170 2.279235e-170
[163,] 1.0000000 3.476067e-169 1.738033e-169
[164,] 1.0000000 3.019281e-168 1.509641e-168
[165,] 1.0000000 1.509561e-167 7.547803e-168
[166,] 1.0000000 1.524008e-166 7.620042e-167
[167,] 1.0000000 1.565598e-165 7.827992e-166
[168,] 1.0000000 1.156650e-164 5.783250e-165
[169,] 1.0000000 1.240400e-163 6.202000e-164
[170,] 1.0000000 7.068506e-163 3.534253e-163
[171,] 1.0000000 6.496466e-162 3.248233e-162
[172,] 1.0000000 6.084647e-161 3.042324e-161
[173,] 1.0000000 6.353754e-160 3.176877e-160
[174,] 1.0000000 6.378297e-159 3.189149e-159
[175,] 1.0000000 6.576200e-158 3.288100e-158
[176,] 1.0000000 6.814731e-157 3.407365e-157
[177,] 1.0000000 6.990730e-156 3.495365e-156
[178,] 1.0000000 1.233696e-155 6.168481e-156
[179,] 1.0000000 1.220846e-154 6.104230e-155
[180,] 1.0000000 1.171999e-153 5.859993e-154
[181,] 1.0000000 1.797175e-153 8.985876e-154
[182,] 1.0000000 1.494595e-152 7.472977e-153
[183,] 1.0000000 1.215747e-151 6.078733e-152
[184,] 1.0000000 1.234686e-150 6.173431e-151
[185,] 1.0000000 1.139060e-149 5.695302e-150
[186,] 1.0000000 1.413391e-149 7.066954e-150
[187,] 1.0000000 1.255418e-148 6.277092e-149
[188,] 1.0000000 9.727058e-148 4.863529e-148
[189,] 1.0000000 9.409626e-147 4.704813e-147
[190,] 1.0000000 1.774533e-146 8.872663e-147
[191,] 1.0000000 1.726426e-145 8.632129e-146
[192,] 1.0000000 1.738832e-144 8.694158e-145
[193,] 1.0000000 1.695318e-143 8.476590e-144
[194,] 1.0000000 1.692170e-142 8.460850e-143
[195,] 1.0000000 1.150242e-141 5.751209e-142
[196,] 1.0000000 1.118025e-140 5.590126e-141
[197,] 1.0000000 1.086300e-139 5.431498e-140
[198,] 1.0000000 1.048246e-138 5.241229e-139
[199,] 1.0000000 9.878389e-138 4.939194e-138
[200,] 1.0000000 9.474896e-137 4.737448e-137
[201,] 1.0000000 4.904683e-136 2.452342e-136
[202,] 1.0000000 3.334568e-135 1.667284e-135
[203,] 1.0000000 3.109045e-134 1.554522e-134
[204,] 1.0000000 2.973265e-133 1.486633e-133
[205,] 1.0000000 2.756539e-132 1.378269e-132
[206,] 1.0000000 2.457399e-131 1.228700e-131
[207,] 1.0000000 2.064627e-130 1.032314e-130
[208,] 1.0000000 1.895560e-129 9.477801e-130
[209,] 1.0000000 1.748845e-128 8.744225e-129
[210,] 1.0000000 1.348922e-127 6.744612e-128
[211,] 1.0000000 1.233142e-126 6.165708e-127
[212,] 1.0000000 1.006987e-125 5.034937e-126
[213,] 1.0000000 8.191082e-125 4.095541e-125
[214,] 1.0000000 7.009679e-124 3.504840e-124
[215,] 1.0000000 6.357440e-123 3.178720e-123
[216,] 1.0000000 5.755164e-122 2.877582e-122
[217,] 1.0000000 5.017235e-121 2.508617e-121
[218,] 1.0000000 4.177225e-120 2.088613e-120
[219,] 1.0000000 3.165921e-119 1.582961e-119
[220,] 1.0000000 2.770133e-118 1.385066e-118
[221,] 1.0000000 2.406088e-117 1.203044e-117
[222,] 1.0000000 2.061616e-116 1.030808e-116
[223,] 1.0000000 1.756536e-115 8.782681e-116
[224,] 1.0000000 1.494084e-114 7.470421e-115
[225,] 1.0000000 1.173694e-113 5.868469e-114
[226,] 1.0000000 9.785125e-113 4.892562e-113
[227,] 1.0000000 8.045456e-112 4.022728e-112
[228,] 1.0000000 6.543373e-111 3.271687e-111
[229,] 1.0000000 5.044873e-110 2.522436e-110
[230,] 1.0000000 3.650419e-109 1.825210e-109
[231,] 1.0000000 2.388237e-108 1.194118e-108
[232,] 1.0000000 1.823016e-107 9.115079e-108
[233,] 1.0000000 1.427722e-106 7.138609e-107
[234,] 1.0000000 1.020411e-105 5.102054e-106
[235,] 1.0000000 7.039147e-105 3.519574e-105
[236,] 1.0000000 1.796720e-104 8.983602e-105
[237,] 1.0000000 1.379959e-103 6.899797e-104
[238,] 1.0000000 1.018043e-102 5.090213e-103
[239,] 1.0000000 6.604199e-102 3.302100e-102
[240,] 1.0000000 5.190228e-101 2.595114e-101
[241,] 1.0000000 4.061107e-100 2.030553e-100
[242,] 1.0000000 2.863600e-99 1.431800e-99
[243,] 1.0000000 1.784347e-98 8.921737e-99
[244,] 1.0000000 1.394565e-97 6.972824e-98
[245,] 1.0000000 1.106766e-96 5.533829e-97
[246,] 1.0000000 7.982078e-96 3.991039e-96
[247,] 1.0000000 6.228839e-95 3.114420e-95
[248,] 1.0000000 4.444483e-94 2.222242e-94
[249,] 1.0000000 3.272412e-93 1.636206e-93
[250,] 1.0000000 2.426491e-92 1.213246e-92
[251,] 1.0000000 1.571147e-91 7.855735e-92
[252,] 1.0000000 1.181357e-90 5.906786e-91
[253,] 1.0000000 8.411522e-90 4.205761e-90
[254,] 1.0000000 5.583746e-89 2.791873e-89
[255,] 1.0000000 3.260437e-88 1.630218e-88
[256,] 1.0000000 2.434540e-87 1.217270e-87
[257,] 1.0000000 1.835846e-86 9.179229e-87
[258,] 1.0000000 1.359016e-85 6.795081e-86
[259,] 1.0000000 1.001033e-84 5.005167e-85
[260,] 1.0000000 5.257473e-84 2.628736e-84
[261,] 1.0000000 3.242805e-83 1.621403e-83
[262,] 1.0000000 2.243271e-82 1.121636e-82
[263,] 1.0000000 1.610920e-81 8.054599e-82
[264,] 1.0000000 1.097873e-80 5.489364e-81
[265,] 1.0000000 7.666484e-80 3.833242e-80
[266,] 1.0000000 3.096284e-79 1.548142e-79
[267,] 1.0000000 1.706625e-78 8.533123e-79
[268,] 1.0000000 1.210495e-77 6.052475e-78
[269,] 1.0000000 8.422003e-77 4.211002e-77
[270,] 1.0000000 5.563761e-76 2.781881e-76
[271,] 1.0000000 3.758791e-75 1.879396e-75
[272,] 1.0000000 2.353066e-74 1.176533e-74
[273,] 1.0000000 1.567976e-73 7.839879e-74
[274,] 1.0000000 7.898929e-73 3.949464e-73
[275,] 1.0000000 5.342107e-72 2.671053e-72
[276,] 1.0000000 3.671467e-71 1.835734e-71
[277,] 1.0000000 2.212815e-70 1.106407e-70
[278,] 1.0000000 1.340599e-69 6.702993e-70
[279,] 1.0000000 6.500872e-69 3.250436e-69
[280,] 1.0000000 4.327674e-68 2.163837e-68
[281,] 1.0000000 2.864225e-67 1.432113e-67
[282,] 1.0000000 1.783834e-66 8.919170e-67
[283,] 1.0000000 1.168735e-65 5.843675e-66
[284,] 1.0000000 7.214148e-65 3.607074e-65
[285,] 1.0000000 4.724043e-64 2.362022e-64
[286,] 1.0000000 2.994938e-63 1.497469e-63
[287,] 1.0000000 1.907197e-62 9.535986e-63
[288,] 1.0000000 1.208518e-61 6.042592e-62
[289,] 1.0000000 6.991633e-61 3.495816e-61
[290,] 1.0000000 4.443278e-60 2.221639e-60
[291,] 1.0000000 1.909895e-59 9.549476e-60
[292,] 1.0000000 1.200422e-58 6.002112e-59
[293,] 1.0000000 3.476208e-58 1.738104e-58
[294,] 1.0000000 2.036656e-57 1.018328e-57
[295,] 1.0000000 1.212292e-56 6.061460e-57
[296,] 1.0000000 6.963111e-56 3.481556e-56
[297,] 1.0000000 4.198624e-55 2.099312e-55
[298,] 1.0000000 2.467746e-54 1.233873e-54
[299,] 1.0000000 1.399564e-53 6.997822e-54
[300,] 1.0000000 8.456213e-53 4.228107e-53
[301,] 1.0000000 4.955160e-52 2.477580e-52
[302,] 1.0000000 2.635257e-51 1.317629e-51
[303,] 1.0000000 1.116421e-50 5.582104e-51
[304,] 1.0000000 6.489547e-50 3.244774e-50
[305,] 1.0000000 3.374513e-49 1.687257e-49
[306,] 1.0000000 1.939148e-48 9.695742e-49
[307,] 1.0000000 1.116631e-47 5.583155e-48
[308,] 1.0000000 5.218546e-47 2.609273e-47
[309,] 1.0000000 2.956542e-46 1.478271e-46
[310,] 1.0000000 1.580848e-45 7.904238e-46
[311,] 1.0000000 8.732961e-45 4.366480e-45
[312,] 1.0000000 4.897232e-44 2.448616e-44
[313,] 1.0000000 2.661885e-43 1.330942e-43
[314,] 1.0000000 1.283929e-42 6.419647e-43
[315,] 1.0000000 4.924936e-42 2.462468e-42
[316,] 1.0000000 2.654068e-41 1.327034e-41
[317,] 1.0000000 1.358982e-40 6.794910e-41
[318,] 1.0000000 7.160235e-40 3.580118e-40
[319,] 1.0000000 3.795607e-39 1.897804e-39
[320,] 1.0000000 1.842006e-38 9.210031e-39
[321,] 1.0000000 9.475968e-38 4.737984e-38
[322,] 1.0000000 4.071590e-37 2.035795e-37
[323,] 1.0000000 2.077340e-36 1.038670e-36
[324,] 1.0000000 1.072282e-35 5.361409e-36
[325,] 1.0000000 5.385084e-35 2.692542e-35
[326,] 1.0000000 2.368721e-34 1.184360e-34
[327,] 1.0000000 7.973038e-34 3.986519e-34
[328,] 1.0000000 4.000591e-33 2.000296e-33
[329,] 1.0000000 1.872919e-32 9.364596e-33
[330,] 1.0000000 9.125407e-32 4.562703e-32
[331,] 1.0000000 4.422509e-31 2.211255e-31
[332,] 1.0000000 2.039787e-30 1.019894e-30
[333,] 1.0000000 9.658328e-30 4.829164e-30
[334,] 1.0000000 4.612269e-29 2.306135e-29
[335,] 1.0000000 2.108958e-28 1.054479e-28
[336,] 1.0000000 9.841508e-28 4.920754e-28
[337,] 1.0000000 4.527169e-27 2.263585e-27
[338,] 1.0000000 8.678941e-27 4.339470e-27
[339,] 1.0000000 2.580537e-26 1.290269e-26
[340,] 1.0000000 1.155544e-25 5.777721e-26
[341,] 1.0000000 4.957932e-25 2.478966e-25
[342,] 1.0000000 2.218917e-24 1.109458e-24
[343,] 1.0000000 9.636545e-24 4.818273e-24
[344,] 1.0000000 3.948443e-23 1.974221e-23
[345,] 1.0000000 1.713758e-22 8.568792e-23
[346,] 1.0000000 6.974291e-22 3.487146e-22
[347,] 1.0000000 2.973028e-21 1.486514e-21
[348,] 1.0000000 1.259580e-20 6.297900e-21
[349,] 1.0000000 5.294536e-20 2.647268e-20
[350,] 1.0000000 1.762634e-19 8.813168e-20
[351,] 1.0000000 4.588891e-19 2.294446e-19
[352,] 1.0000000 1.805793e-18 9.028967e-19
[353,] 1.0000000 6.736422e-18 3.368211e-18
[354,] 1.0000000 2.083866e-17 1.041933e-17
[355,] 1.0000000 8.367980e-17 4.183990e-17
[356,] 1.0000000 3.054858e-16 1.527429e-16
[357,] 1.0000000 1.159483e-15 5.797417e-16
[358,] 1.0000000 4.277151e-15 2.138576e-15
[359,] 1.0000000 1.629658e-14 8.148288e-15
[360,] 1.0000000 6.035324e-14 3.017662e-14
[361,] 1.0000000 2.249465e-13 1.124733e-13
[362,] 1.0000000 6.604261e-13 3.302131e-13
[363,] 1.0000000 8.515367e-13 4.257684e-13
[364,] 1.0000000 2.423195e-12 1.211598e-12
[365,] 1.0000000 8.440083e-12 4.220041e-12
[366,] 1.0000000 2.950632e-11 1.475316e-11
[367,] 1.0000000 7.810577e-11 3.905288e-11
[368,] 1.0000000 2.638444e-10 1.319222e-10
[369,] 1.0000000 7.710004e-10 3.855002e-10
[370,] 1.0000000 2.523063e-09 1.261532e-09
[371,] 1.0000000 8.418425e-09 4.209213e-09
[372,] 1.0000000 2.595268e-08 1.297634e-08
[373,] 1.0000000 8.419120e-08 4.209560e-08
[374,] 0.9999999 1.867601e-07 9.338007e-08
[375,] 0.9999998 4.662505e-07 2.331252e-07
[376,] 0.9999995 1.072430e-06 5.362152e-07
[377,] 0.9999986 2.888089e-06 1.444044e-06
[378,] 0.9999958 8.306592e-06 4.153296e-06
[379,] 0.9999891 2.182516e-05 1.091258e-05
[380,] 0.9999714 5.729000e-05 2.864500e-05
[381,] 0.9999295 1.410838e-04 7.054190e-05
[382,] 0.9998171 3.658612e-04 1.829306e-04
[383,] 0.9995081 9.837792e-04 4.918896e-04
[384,] 0.9988479 2.304145e-03 1.152073e-03
[385,] 0.9972410 5.517943e-03 2.758971e-03
[386,] 0.9952724 9.455210e-03 4.727605e-03
[387,] 0.9946851 1.062987e-02 5.314934e-03
[388,] 0.9864874 2.702516e-02 1.351258e-02
[389,] 0.9689986 6.200287e-02 3.100143e-02
[390,] 0.9318243 1.363515e-01 6.817573e-02
[391,] 0.8780505 2.438991e-01 1.219495e-01
[392,] 0.8272068 3.455863e-01 1.727932e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/1ptne1291311351.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/2il4z1291311351.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/3il4z1291311351.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/4il4z1291311351.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/5su321291311351.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
1559382.1240 3219486.3610 1316976.5266 -3006346.2188 -367191.0825
6 7 8 9 10
58124.8679 2212.4809 613640.7565 -244097.0091 767206.3699
11 12 13 14 15
-276783.0698 469034.9237 912111.2466 -1813656.1773 17121.3969
16 17 18 19 20
561589.9066 -392968.6550 -64339.2419 -147609.0330 -307696.1383
21 22 23 24 25
139641.1243 -248941.0390 -187930.2372 -502685.5525 -273494.3746
26 27 28 29 30
-165338.3774 -397060.1855 -314108.1510 438513.2652 40866.8211
31 32 33 34 35
174558.3890 -147551.9334 98489.0883 9333.7365 -297071.3273
36 37 38 39 40
-465712.7664 -503.6238 21502.3355 -162914.9788 -26382.5030
41 42 43 44 45
39197.1954 -336941.6386 102970.6058 -164200.1884 -37599.3789
46 47 48 49 50
-240022.4042 -211638.3635 -198733.9311 -162764.7183 -198167.6920
51 52 53 54 55
-10172.6875 -15652.9742 -391777.7316 -47008.1210 -407114.0059
56 57 58 59 60
158694.2498 157658.4218 -17093.0133 29662.1339 -246820.5820
61 62 63 64 65
-154987.1448 -98296.5473 -143496.3063 -117930.2557 94141.3178
66 67 68 69 70
85054.1162 -91122.9404 -86610.8068 -89248.0610 -145183.7168
71 72 73 74 75
16918.9103 63495.4408 -184187.0440 -57447.3363 -149328.2290
76 77 78 79 80
-27035.5362 -66459.5557 -31848.2961 87945.1537 -46679.2297
81 82 83 84 85
25210.4788 -107039.8571 -71389.2716 79001.3082 -64219.7230
86 87 88 89 90
-162390.0116 3461.1806 81721.3766 -70002.3326 -6295.7068
91 92 93 94 95
-26834.3556 -89252.9649 108382.8758 -350163.4635 195354.9034
96 97 98 99 100
-42757.0653 -47277.4307 -6488.8067 -161515.9341 27978.6061
101 102 103 104 105
-75472.9421 -87156.7202 -112159.6626 -33603.8601 31825.0501
106 107 108 109 110
-43708.3939 -196431.7739 -24151.3414 -262899.4675 -204462.3356
111 112 113 114 115
90212.0722 -120314.5081 -241889.4258 -168288.3944 -15188.6548
116 117 118 119 120
102515.1909 -17667.6225 40424.0791 -195303.7037 64457.7881
121 122 123 124 125
-320076.0226 775.4535 -26436.4370 62494.9877 -102103.9437
126 127 128 129 130
-244215.1272 57619.6952 79535.7074 36760.3084 57216.9115
131 132 133 134 135
110945.1043 101007.1573 -1135.4498 11317.4683 30592.5650
136 137 138 139 140
181231.2433 86246.7689 71420.8796 57220.3376 78925.5282
141 142 143 144 145
27326.2681 56817.5539 102382.0307 43576.6836 4460.6000
146 147 148 149 150
25076.9885 30193.2074 180204.3464 83679.0830 70427.3078
151 152 153 154 155
47520.2898 67175.5227 27926.6053 34196.0012 107085.0498
156 157 158 159 160
42950.8845 -1934.1650 9361.9570 36912.3878 177922.6656
161 162 163 164 165
83279.7253 70622.1643 56421.6223 83601.2759 27524.5756
166 167 168 169 170
56018.8386 106640.0504 146863.4910 -2333.5227 12187.0581
171 172 173 174 175
29394.4921 174437.8893 82880.3677 67823.9325 56022.2646
176 177 178 179 180
77666.3580 27127.8901 -25384.8721 106240.6928 42152.1693
181 182 183 184 185
10572.2577 -91992.7523 54567.9171 34750.1699 1860.2953
186 187 188 189 190
-17583.4095 -19283.0027 73886.6124 -13283.6223 55220.1234
191 192 193 194 195
110553.7260 41752.8116 -1075.2154 11388.3428 28291.4501
196 197 198 199 200
96973.8085 -137338.0361 -10326.6928 -24527.2348 -82017.9469
201 202 203 204 205
-53421.6094 -24930.0185 25691.1933 -9968.9404 -171156.9114
206 207 208 209 210
-70568.1759 -67041.3728 96574.4509 61264.8860 -10726.0504
211 212 213 214 215
54824.1917 76529.3824 41602.3143 31945.2229 86180.3762
216 217 218 219 220
40954.0963 -481.6499 14565.1699 27797.0616 180964.6773
221 222 223 224 225
48481.7709 68625.3761 55310.6017 76130.0247 25815.3369
226 227 228 229 230
58518.9376 104643.2622 40554.7387 -30779.0600 10190.2699
231 232 233 234 235
6028.6918 175526.5198 80883.5795 69429.7276 54025.4765
236 237 238 239 240
75977.6294 25131.1019 56557.0814 104243.9046 40155.3811
241 242 243 244 245
-4729.6685 12619.6182 26998.3463 180305.4857 -7034.8025
246 247 248 249 250
-11924.1233 -26124.6654 67052.9749 -55019.0399 -26527.4490
251 252 253 254 255
111246.1796 -39994.7608 -28722.6998 -11443.7872 -26268.4696
256 257 258 259 260
174727.8045 103716.9431 107681.0092 -26524.0230 74931.9518
261 262 263 264 265
-55418.3976 -26926.8066 83049.3179 47905.3055 -85279.1680
266 267 268 269 270
8992.1970 -53551.1532 176653.6049 -65.2776 44960.5327
271 272 273 274 275
-26923.3806 -5218.1900 -55817.7552 50841.1017 103045.8317
276 277 278 279 280
47762.1019 654.0813 8592.8394 -66396.4955 152290.5768
281 282 283 284 285
-464.6352 66628.5880 90639.6794 74133.2365 -12959.0918
286 287 288 289 290
52025.2623 22895.6899 37480.5295 -86077.8832 38632.1684
291 292 293 294 295
-54349.8684 182986.2128 50920.9047 66229.2303 55605.5674
296 297 298 299 300
28578.8406 23134.3137 -28124.8796 102247.1164 38158.5929
301 302 303 304 305
-87234.8524 -6156.9831 25001.5581 93379.5898 116987.5958
306 307 308 309 310
69851.0939 -28121.4535 -14683.7381 44669.5805 -28524.2372
311 312 313 314 315
22096.9746 15728.1348 -86876.5985 -75580.8618 26509.6285
316 317 318 319 320
173492.7836 79709.6339 51181.9059 53935.6237 72935.1636
321 322 323 324 325
22335.5984 50827.1894 101448.4011 -42390.9066 -7525.1720
326 327 328 329 330
4409.5701 -55547.9414 92580.8745 77688.7184 65031.1574
331 332 333 334 335
50830.6154 72535.8060 21936.2408 50430.3128 101915.3553
336 337 338 339 340
36960.5200 -9149.9613 6596.0512 23138.1773 175892.7380
341 342 343 344 345
57987.9338 64631.7998 50431.2577 72136.4484 21536.8832
346 347 348 349 350
50028.4741 92659.4641 36561.1623 -8323.8872 6196.6935
351 352 353 354 355
24131.6419 172759.1919 78322.9168 54266.9516 50031.9001
356 357 358 359 360
71737.0907 -13162.0903 49629.1164 100250.3282 36161.8047
361 362 363 364 365
-6936.7112 5797.3359 23004.7699 82172.0736 75012.8225
366 367 368 369 370
63833.0845 49632.5425 71337.7331 20738.1679 49229.7588
371 372 373 374 375
99850.9706 38492.1934 -49897.4353 -21610.5108 22605.4123
376 377 378 379 380
170734.2282 88436.6202 48186.2337 49233.1848 70938.3755
381 382 383 384 385
20338.8102 40858.6605 152858.5696 35363.0894 -7636.9399
386 387 388 389 390
-152062.9544 26878.5157 154494.8296 73866.2132 43678.0849
391 392 393 394 395
45460.1314 50020.7238 -17867.1035 64608.0722 -192929.4043
396 397 398 399 400
7325.0672 -173842.2880 -66961.9915 -4978.1456 84161.6197
401 402 403 404 405
36945.7457 52253.7733 27767.4174 -92345.5438 -14294.8941
406 407 408 409 410
-4643.4801 82458.6994 14013.2055 -32672.7658 -71588.7820
411 412 413 414 415
-24342.8319 -114359.6437 -34348.3071 -158884.3866 -64220.2473
416 417 418 419 420
-25464.5040 101931.1580 -85882.5624 46592.5195 -34652.7350
421 422 423 424 425
-132968.7551 -153473.7950 -466415.9626 -356872.4710 -52907.5755
426 427 428 429 430
-205270.7298 -314466.3688 -1195291.5387 -397186.5168 -278836.6108
431
-799679.6048
> postscript(file="/var/www/html/freestat/rcomp/tmp/6su321291311351.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 1559382.1240 NA
1 3219486.3610 1559382.1240
2 1316976.5266 3219486.3610
3 -3006346.2188 1316976.5266
4 -367191.0825 -3006346.2188
5 58124.8679 -367191.0825
6 2212.4809 58124.8679
7 613640.7565 2212.4809
8 -244097.0091 613640.7565
9 767206.3699 -244097.0091
10 -276783.0698 767206.3699
11 469034.9237 -276783.0698
12 912111.2466 469034.9237
13 -1813656.1773 912111.2466
14 17121.3969 -1813656.1773
15 561589.9066 17121.3969
16 -392968.6550 561589.9066
17 -64339.2419 -392968.6550
18 -147609.0330 -64339.2419
19 -307696.1383 -147609.0330
20 139641.1243 -307696.1383
21 -248941.0390 139641.1243
22 -187930.2372 -248941.0390
23 -502685.5525 -187930.2372
24 -273494.3746 -502685.5525
25 -165338.3774 -273494.3746
26 -397060.1855 -165338.3774
27 -314108.1510 -397060.1855
28 438513.2652 -314108.1510
29 40866.8211 438513.2652
30 174558.3890 40866.8211
31 -147551.9334 174558.3890
32 98489.0883 -147551.9334
33 9333.7365 98489.0883
34 -297071.3273 9333.7365
35 -465712.7664 -297071.3273
36 -503.6238 -465712.7664
37 21502.3355 -503.6238
38 -162914.9788 21502.3355
39 -26382.5030 -162914.9788
40 39197.1954 -26382.5030
41 -336941.6386 39197.1954
42 102970.6058 -336941.6386
43 -164200.1884 102970.6058
44 -37599.3789 -164200.1884
45 -240022.4042 -37599.3789
46 -211638.3635 -240022.4042
47 -198733.9311 -211638.3635
48 -162764.7183 -198733.9311
49 -198167.6920 -162764.7183
50 -10172.6875 -198167.6920
51 -15652.9742 -10172.6875
52 -391777.7316 -15652.9742
53 -47008.1210 -391777.7316
54 -407114.0059 -47008.1210
55 158694.2498 -407114.0059
56 157658.4218 158694.2498
57 -17093.0133 157658.4218
58 29662.1339 -17093.0133
59 -246820.5820 29662.1339
60 -154987.1448 -246820.5820
61 -98296.5473 -154987.1448
62 -143496.3063 -98296.5473
63 -117930.2557 -143496.3063
64 94141.3178 -117930.2557
65 85054.1162 94141.3178
66 -91122.9404 85054.1162
67 -86610.8068 -91122.9404
68 -89248.0610 -86610.8068
69 -145183.7168 -89248.0610
70 16918.9103 -145183.7168
71 63495.4408 16918.9103
72 -184187.0440 63495.4408
73 -57447.3363 -184187.0440
74 -149328.2290 -57447.3363
75 -27035.5362 -149328.2290
76 -66459.5557 -27035.5362
77 -31848.2961 -66459.5557
78 87945.1537 -31848.2961
79 -46679.2297 87945.1537
80 25210.4788 -46679.2297
81 -107039.8571 25210.4788
82 -71389.2716 -107039.8571
83 79001.3082 -71389.2716
84 -64219.7230 79001.3082
85 -162390.0116 -64219.7230
86 3461.1806 -162390.0116
87 81721.3766 3461.1806
88 -70002.3326 81721.3766
89 -6295.7068 -70002.3326
90 -26834.3556 -6295.7068
91 -89252.9649 -26834.3556
92 108382.8758 -89252.9649
93 -350163.4635 108382.8758
94 195354.9034 -350163.4635
95 -42757.0653 195354.9034
96 -47277.4307 -42757.0653
97 -6488.8067 -47277.4307
98 -161515.9341 -6488.8067
99 27978.6061 -161515.9341
100 -75472.9421 27978.6061
101 -87156.7202 -75472.9421
102 -112159.6626 -87156.7202
103 -33603.8601 -112159.6626
104 31825.0501 -33603.8601
105 -43708.3939 31825.0501
106 -196431.7739 -43708.3939
107 -24151.3414 -196431.7739
108 -262899.4675 -24151.3414
109 -204462.3356 -262899.4675
110 90212.0722 -204462.3356
111 -120314.5081 90212.0722
112 -241889.4258 -120314.5081
113 -168288.3944 -241889.4258
114 -15188.6548 -168288.3944
115 102515.1909 -15188.6548
116 -17667.6225 102515.1909
117 40424.0791 -17667.6225
118 -195303.7037 40424.0791
119 64457.7881 -195303.7037
120 -320076.0226 64457.7881
121 775.4535 -320076.0226
122 -26436.4370 775.4535
123 62494.9877 -26436.4370
124 -102103.9437 62494.9877
125 -244215.1272 -102103.9437
126 57619.6952 -244215.1272
127 79535.7074 57619.6952
128 36760.3084 79535.7074
129 57216.9115 36760.3084
130 110945.1043 57216.9115
131 101007.1573 110945.1043
132 -1135.4498 101007.1573
133 11317.4683 -1135.4498
134 30592.5650 11317.4683
135 181231.2433 30592.5650
136 86246.7689 181231.2433
137 71420.8796 86246.7689
138 57220.3376 71420.8796
139 78925.5282 57220.3376
140 27326.2681 78925.5282
141 56817.5539 27326.2681
142 102382.0307 56817.5539
143 43576.6836 102382.0307
144 4460.6000 43576.6836
145 25076.9885 4460.6000
146 30193.2074 25076.9885
147 180204.3464 30193.2074
148 83679.0830 180204.3464
149 70427.3078 83679.0830
150 47520.2898 70427.3078
151 67175.5227 47520.2898
152 27926.6053 67175.5227
153 34196.0012 27926.6053
154 107085.0498 34196.0012
155 42950.8845 107085.0498
156 -1934.1650 42950.8845
157 9361.9570 -1934.1650
158 36912.3878 9361.9570
159 177922.6656 36912.3878
160 83279.7253 177922.6656
161 70622.1643 83279.7253
162 56421.6223 70622.1643
163 83601.2759 56421.6223
164 27524.5756 83601.2759
165 56018.8386 27524.5756
166 106640.0504 56018.8386
167 146863.4910 106640.0504
168 -2333.5227 146863.4910
169 12187.0581 -2333.5227
170 29394.4921 12187.0581
171 174437.8893 29394.4921
172 82880.3677 174437.8893
173 67823.9325 82880.3677
174 56022.2646 67823.9325
175 77666.3580 56022.2646
176 27127.8901 77666.3580
177 -25384.8721 27127.8901
178 106240.6928 -25384.8721
179 42152.1693 106240.6928
180 10572.2577 42152.1693
181 -91992.7523 10572.2577
182 54567.9171 -91992.7523
183 34750.1699 54567.9171
184 1860.2953 34750.1699
185 -17583.4095 1860.2953
186 -19283.0027 -17583.4095
187 73886.6124 -19283.0027
188 -13283.6223 73886.6124
189 55220.1234 -13283.6223
190 110553.7260 55220.1234
191 41752.8116 110553.7260
192 -1075.2154 41752.8116
193 11388.3428 -1075.2154
194 28291.4501 11388.3428
195 96973.8085 28291.4501
196 -137338.0361 96973.8085
197 -10326.6928 -137338.0361
198 -24527.2348 -10326.6928
199 -82017.9469 -24527.2348
200 -53421.6094 -82017.9469
201 -24930.0185 -53421.6094
202 25691.1933 -24930.0185
203 -9968.9404 25691.1933
204 -171156.9114 -9968.9404
205 -70568.1759 -171156.9114
206 -67041.3728 -70568.1759
207 96574.4509 -67041.3728
208 61264.8860 96574.4509
209 -10726.0504 61264.8860
210 54824.1917 -10726.0504
211 76529.3824 54824.1917
212 41602.3143 76529.3824
213 31945.2229 41602.3143
214 86180.3762 31945.2229
215 40954.0963 86180.3762
216 -481.6499 40954.0963
217 14565.1699 -481.6499
218 27797.0616 14565.1699
219 180964.6773 27797.0616
220 48481.7709 180964.6773
221 68625.3761 48481.7709
222 55310.6017 68625.3761
223 76130.0247 55310.6017
224 25815.3369 76130.0247
225 58518.9376 25815.3369
226 104643.2622 58518.9376
227 40554.7387 104643.2622
228 -30779.0600 40554.7387
229 10190.2699 -30779.0600
230 6028.6918 10190.2699
231 175526.5198 6028.6918
232 80883.5795 175526.5198
233 69429.7276 80883.5795
234 54025.4765 69429.7276
235 75977.6294 54025.4765
236 25131.1019 75977.6294
237 56557.0814 25131.1019
238 104243.9046 56557.0814
239 40155.3811 104243.9046
240 -4729.6685 40155.3811
241 12619.6182 -4729.6685
242 26998.3463 12619.6182
243 180305.4857 26998.3463
244 -7034.8025 180305.4857
245 -11924.1233 -7034.8025
246 -26124.6654 -11924.1233
247 67052.9749 -26124.6654
248 -55019.0399 67052.9749
249 -26527.4490 -55019.0399
250 111246.1796 -26527.4490
251 -39994.7608 111246.1796
252 -28722.6998 -39994.7608
253 -11443.7872 -28722.6998
254 -26268.4696 -11443.7872
255 174727.8045 -26268.4696
256 103716.9431 174727.8045
257 107681.0092 103716.9431
258 -26524.0230 107681.0092
259 74931.9518 -26524.0230
260 -55418.3976 74931.9518
261 -26926.8066 -55418.3976
262 83049.3179 -26926.8066
263 47905.3055 83049.3179
264 -85279.1680 47905.3055
265 8992.1970 -85279.1680
266 -53551.1532 8992.1970
267 176653.6049 -53551.1532
268 -65.2776 176653.6049
269 44960.5327 -65.2776
270 -26923.3806 44960.5327
271 -5218.1900 -26923.3806
272 -55817.7552 -5218.1900
273 50841.1017 -55817.7552
274 103045.8317 50841.1017
275 47762.1019 103045.8317
276 654.0813 47762.1019
277 8592.8394 654.0813
278 -66396.4955 8592.8394
279 152290.5768 -66396.4955
280 -464.6352 152290.5768
281 66628.5880 -464.6352
282 90639.6794 66628.5880
283 74133.2365 90639.6794
284 -12959.0918 74133.2365
285 52025.2623 -12959.0918
286 22895.6899 52025.2623
287 37480.5295 22895.6899
288 -86077.8832 37480.5295
289 38632.1684 -86077.8832
290 -54349.8684 38632.1684
291 182986.2128 -54349.8684
292 50920.9047 182986.2128
293 66229.2303 50920.9047
294 55605.5674 66229.2303
295 28578.8406 55605.5674
296 23134.3137 28578.8406
297 -28124.8796 23134.3137
298 102247.1164 -28124.8796
299 38158.5929 102247.1164
300 -87234.8524 38158.5929
301 -6156.9831 -87234.8524
302 25001.5581 -6156.9831
303 93379.5898 25001.5581
304 116987.5958 93379.5898
305 69851.0939 116987.5958
306 -28121.4535 69851.0939
307 -14683.7381 -28121.4535
308 44669.5805 -14683.7381
309 -28524.2372 44669.5805
310 22096.9746 -28524.2372
311 15728.1348 22096.9746
312 -86876.5985 15728.1348
313 -75580.8618 -86876.5985
314 26509.6285 -75580.8618
315 173492.7836 26509.6285
316 79709.6339 173492.7836
317 51181.9059 79709.6339
318 53935.6237 51181.9059
319 72935.1636 53935.6237
320 22335.5984 72935.1636
321 50827.1894 22335.5984
322 101448.4011 50827.1894
323 -42390.9066 101448.4011
324 -7525.1720 -42390.9066
325 4409.5701 -7525.1720
326 -55547.9414 4409.5701
327 92580.8745 -55547.9414
328 77688.7184 92580.8745
329 65031.1574 77688.7184
330 50830.6154 65031.1574
331 72535.8060 50830.6154
332 21936.2408 72535.8060
333 50430.3128 21936.2408
334 101915.3553 50430.3128
335 36960.5200 101915.3553
336 -9149.9613 36960.5200
337 6596.0512 -9149.9613
338 23138.1773 6596.0512
339 175892.7380 23138.1773
340 57987.9338 175892.7380
341 64631.7998 57987.9338
342 50431.2577 64631.7998
343 72136.4484 50431.2577
344 21536.8832 72136.4484
345 50028.4741 21536.8832
346 92659.4641 50028.4741
347 36561.1623 92659.4641
348 -8323.8872 36561.1623
349 6196.6935 -8323.8872
350 24131.6419 6196.6935
351 172759.1919 24131.6419
352 78322.9168 172759.1919
353 54266.9516 78322.9168
354 50031.9001 54266.9516
355 71737.0907 50031.9001
356 -13162.0903 71737.0907
357 49629.1164 -13162.0903
358 100250.3282 49629.1164
359 36161.8047 100250.3282
360 -6936.7112 36161.8047
361 5797.3359 -6936.7112
362 23004.7699 5797.3359
363 82172.0736 23004.7699
364 75012.8225 82172.0736
365 63833.0845 75012.8225
366 49632.5425 63833.0845
367 71337.7331 49632.5425
368 20738.1679 71337.7331
369 49229.7588 20738.1679
370 99850.9706 49229.7588
371 38492.1934 99850.9706
372 -49897.4353 38492.1934
373 -21610.5108 -49897.4353
374 22605.4123 -21610.5108
375 170734.2282 22605.4123
376 88436.6202 170734.2282
377 48186.2337 88436.6202
378 49233.1848 48186.2337
379 70938.3755 49233.1848
380 20338.8102 70938.3755
381 40858.6605 20338.8102
382 152858.5696 40858.6605
383 35363.0894 152858.5696
384 -7636.9399 35363.0894
385 -152062.9544 -7636.9399
386 26878.5157 -152062.9544
387 154494.8296 26878.5157
388 73866.2132 154494.8296
389 43678.0849 73866.2132
390 45460.1314 43678.0849
391 50020.7238 45460.1314
392 -17867.1035 50020.7238
393 64608.0722 -17867.1035
394 -192929.4043 64608.0722
395 7325.0672 -192929.4043
396 -173842.2880 7325.0672
397 -66961.9915 -173842.2880
398 -4978.1456 -66961.9915
399 84161.6197 -4978.1456
400 36945.7457 84161.6197
401 52253.7733 36945.7457
402 27767.4174 52253.7733
403 -92345.5438 27767.4174
404 -14294.8941 -92345.5438
405 -4643.4801 -14294.8941
406 82458.6994 -4643.4801
407 14013.2055 82458.6994
408 -32672.7658 14013.2055
409 -71588.7820 -32672.7658
410 -24342.8319 -71588.7820
411 -114359.6437 -24342.8319
412 -34348.3071 -114359.6437
413 -158884.3866 -34348.3071
414 -64220.2473 -158884.3866
415 -25464.5040 -64220.2473
416 101931.1580 -25464.5040
417 -85882.5624 101931.1580
418 46592.5195 -85882.5624
419 -34652.7350 46592.5195
420 -132968.7551 -34652.7350
421 -153473.7950 -132968.7551
422 -466415.9626 -153473.7950
423 -356872.4710 -466415.9626
424 -52907.5755 -356872.4710
425 -205270.7298 -52907.5755
426 -314466.3688 -205270.7298
427 -1195291.5387 -314466.3688
428 -397186.5168 -1195291.5387
429 -278836.6108 -397186.5168
430 -799679.6048 -278836.6108
431 NA -799679.6048
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3219486.3610 1559382.1240
[2,] 1316976.5266 3219486.3610
[3,] -3006346.2188 1316976.5266
[4,] -367191.0825 -3006346.2188
[5,] 58124.8679 -367191.0825
[6,] 2212.4809 58124.8679
[7,] 613640.7565 2212.4809
[8,] -244097.0091 613640.7565
[9,] 767206.3699 -244097.0091
[10,] -276783.0698 767206.3699
[11,] 469034.9237 -276783.0698
[12,] 912111.2466 469034.9237
[13,] -1813656.1773 912111.2466
[14,] 17121.3969 -1813656.1773
[15,] 561589.9066 17121.3969
[16,] -392968.6550 561589.9066
[17,] -64339.2419 -392968.6550
[18,] -147609.0330 -64339.2419
[19,] -307696.1383 -147609.0330
[20,] 139641.1243 -307696.1383
[21,] -248941.0390 139641.1243
[22,] -187930.2372 -248941.0390
[23,] -502685.5525 -187930.2372
[24,] -273494.3746 -502685.5525
[25,] -165338.3774 -273494.3746
[26,] -397060.1855 -165338.3774
[27,] -314108.1510 -397060.1855
[28,] 438513.2652 -314108.1510
[29,] 40866.8211 438513.2652
[30,] 174558.3890 40866.8211
[31,] -147551.9334 174558.3890
[32,] 98489.0883 -147551.9334
[33,] 9333.7365 98489.0883
[34,] -297071.3273 9333.7365
[35,] -465712.7664 -297071.3273
[36,] -503.6238 -465712.7664
[37,] 21502.3355 -503.6238
[38,] -162914.9788 21502.3355
[39,] -26382.5030 -162914.9788
[40,] 39197.1954 -26382.5030
[41,] -336941.6386 39197.1954
[42,] 102970.6058 -336941.6386
[43,] -164200.1884 102970.6058
[44,] -37599.3789 -164200.1884
[45,] -240022.4042 -37599.3789
[46,] -211638.3635 -240022.4042
[47,] -198733.9311 -211638.3635
[48,] -162764.7183 -198733.9311
[49,] -198167.6920 -162764.7183
[50,] -10172.6875 -198167.6920
[51,] -15652.9742 -10172.6875
[52,] -391777.7316 -15652.9742
[53,] -47008.1210 -391777.7316
[54,] -407114.0059 -47008.1210
[55,] 158694.2498 -407114.0059
[56,] 157658.4218 158694.2498
[57,] -17093.0133 157658.4218
[58,] 29662.1339 -17093.0133
[59,] -246820.5820 29662.1339
[60,] -154987.1448 -246820.5820
[61,] -98296.5473 -154987.1448
[62,] -143496.3063 -98296.5473
[63,] -117930.2557 -143496.3063
[64,] 94141.3178 -117930.2557
[65,] 85054.1162 94141.3178
[66,] -91122.9404 85054.1162
[67,] -86610.8068 -91122.9404
[68,] -89248.0610 -86610.8068
[69,] -145183.7168 -89248.0610
[70,] 16918.9103 -145183.7168
[71,] 63495.4408 16918.9103
[72,] -184187.0440 63495.4408
[73,] -57447.3363 -184187.0440
[74,] -149328.2290 -57447.3363
[75,] -27035.5362 -149328.2290
[76,] -66459.5557 -27035.5362
[77,] -31848.2961 -66459.5557
[78,] 87945.1537 -31848.2961
[79,] -46679.2297 87945.1537
[80,] 25210.4788 -46679.2297
[81,] -107039.8571 25210.4788
[82,] -71389.2716 -107039.8571
[83,] 79001.3082 -71389.2716
[84,] -64219.7230 79001.3082
[85,] -162390.0116 -64219.7230
[86,] 3461.1806 -162390.0116
[87,] 81721.3766 3461.1806
[88,] -70002.3326 81721.3766
[89,] -6295.7068 -70002.3326
[90,] -26834.3556 -6295.7068
[91,] -89252.9649 -26834.3556
[92,] 108382.8758 -89252.9649
[93,] -350163.4635 108382.8758
[94,] 195354.9034 -350163.4635
[95,] -42757.0653 195354.9034
[96,] -47277.4307 -42757.0653
[97,] -6488.8067 -47277.4307
[98,] -161515.9341 -6488.8067
[99,] 27978.6061 -161515.9341
[100,] -75472.9421 27978.6061
[101,] -87156.7202 -75472.9421
[102,] -112159.6626 -87156.7202
[103,] -33603.8601 -112159.6626
[104,] 31825.0501 -33603.8601
[105,] -43708.3939 31825.0501
[106,] -196431.7739 -43708.3939
[107,] -24151.3414 -196431.7739
[108,] -262899.4675 -24151.3414
[109,] -204462.3356 -262899.4675
[110,] 90212.0722 -204462.3356
[111,] -120314.5081 90212.0722
[112,] -241889.4258 -120314.5081
[113,] -168288.3944 -241889.4258
[114,] -15188.6548 -168288.3944
[115,] 102515.1909 -15188.6548
[116,] -17667.6225 102515.1909
[117,] 40424.0791 -17667.6225
[118,] -195303.7037 40424.0791
[119,] 64457.7881 -195303.7037
[120,] -320076.0226 64457.7881
[121,] 775.4535 -320076.0226
[122,] -26436.4370 775.4535
[123,] 62494.9877 -26436.4370
[124,] -102103.9437 62494.9877
[125,] -244215.1272 -102103.9437
[126,] 57619.6952 -244215.1272
[127,] 79535.7074 57619.6952
[128,] 36760.3084 79535.7074
[129,] 57216.9115 36760.3084
[130,] 110945.1043 57216.9115
[131,] 101007.1573 110945.1043
[132,] -1135.4498 101007.1573
[133,] 11317.4683 -1135.4498
[134,] 30592.5650 11317.4683
[135,] 181231.2433 30592.5650
[136,] 86246.7689 181231.2433
[137,] 71420.8796 86246.7689
[138,] 57220.3376 71420.8796
[139,] 78925.5282 57220.3376
[140,] 27326.2681 78925.5282
[141,] 56817.5539 27326.2681
[142,] 102382.0307 56817.5539
[143,] 43576.6836 102382.0307
[144,] 4460.6000 43576.6836
[145,] 25076.9885 4460.6000
[146,] 30193.2074 25076.9885
[147,] 180204.3464 30193.2074
[148,] 83679.0830 180204.3464
[149,] 70427.3078 83679.0830
[150,] 47520.2898 70427.3078
[151,] 67175.5227 47520.2898
[152,] 27926.6053 67175.5227
[153,] 34196.0012 27926.6053
[154,] 107085.0498 34196.0012
[155,] 42950.8845 107085.0498
[156,] -1934.1650 42950.8845
[157,] 9361.9570 -1934.1650
[158,] 36912.3878 9361.9570
[159,] 177922.6656 36912.3878
[160,] 83279.7253 177922.6656
[161,] 70622.1643 83279.7253
[162,] 56421.6223 70622.1643
[163,] 83601.2759 56421.6223
[164,] 27524.5756 83601.2759
[165,] 56018.8386 27524.5756
[166,] 106640.0504 56018.8386
[167,] 146863.4910 106640.0504
[168,] -2333.5227 146863.4910
[169,] 12187.0581 -2333.5227
[170,] 29394.4921 12187.0581
[171,] 174437.8893 29394.4921
[172,] 82880.3677 174437.8893
[173,] 67823.9325 82880.3677
[174,] 56022.2646 67823.9325
[175,] 77666.3580 56022.2646
[176,] 27127.8901 77666.3580
[177,] -25384.8721 27127.8901
[178,] 106240.6928 -25384.8721
[179,] 42152.1693 106240.6928
[180,] 10572.2577 42152.1693
[181,] -91992.7523 10572.2577
[182,] 54567.9171 -91992.7523
[183,] 34750.1699 54567.9171
[184,] 1860.2953 34750.1699
[185,] -17583.4095 1860.2953
[186,] -19283.0027 -17583.4095
[187,] 73886.6124 -19283.0027
[188,] -13283.6223 73886.6124
[189,] 55220.1234 -13283.6223
[190,] 110553.7260 55220.1234
[191,] 41752.8116 110553.7260
[192,] -1075.2154 41752.8116
[193,] 11388.3428 -1075.2154
[194,] 28291.4501 11388.3428
[195,] 96973.8085 28291.4501
[196,] -137338.0361 96973.8085
[197,] -10326.6928 -137338.0361
[198,] -24527.2348 -10326.6928
[199,] -82017.9469 -24527.2348
[200,] -53421.6094 -82017.9469
[201,] -24930.0185 -53421.6094
[202,] 25691.1933 -24930.0185
[203,] -9968.9404 25691.1933
[204,] -171156.9114 -9968.9404
[205,] -70568.1759 -171156.9114
[206,] -67041.3728 -70568.1759
[207,] 96574.4509 -67041.3728
[208,] 61264.8860 96574.4509
[209,] -10726.0504 61264.8860
[210,] 54824.1917 -10726.0504
[211,] 76529.3824 54824.1917
[212,] 41602.3143 76529.3824
[213,] 31945.2229 41602.3143
[214,] 86180.3762 31945.2229
[215,] 40954.0963 86180.3762
[216,] -481.6499 40954.0963
[217,] 14565.1699 -481.6499
[218,] 27797.0616 14565.1699
[219,] 180964.6773 27797.0616
[220,] 48481.7709 180964.6773
[221,] 68625.3761 48481.7709
[222,] 55310.6017 68625.3761
[223,] 76130.0247 55310.6017
[224,] 25815.3369 76130.0247
[225,] 58518.9376 25815.3369
[226,] 104643.2622 58518.9376
[227,] 40554.7387 104643.2622
[228,] -30779.0600 40554.7387
[229,] 10190.2699 -30779.0600
[230,] 6028.6918 10190.2699
[231,] 175526.5198 6028.6918
[232,] 80883.5795 175526.5198
[233,] 69429.7276 80883.5795
[234,] 54025.4765 69429.7276
[235,] 75977.6294 54025.4765
[236,] 25131.1019 75977.6294
[237,] 56557.0814 25131.1019
[238,] 104243.9046 56557.0814
[239,] 40155.3811 104243.9046
[240,] -4729.6685 40155.3811
[241,] 12619.6182 -4729.6685
[242,] 26998.3463 12619.6182
[243,] 180305.4857 26998.3463
[244,] -7034.8025 180305.4857
[245,] -11924.1233 -7034.8025
[246,] -26124.6654 -11924.1233
[247,] 67052.9749 -26124.6654
[248,] -55019.0399 67052.9749
[249,] -26527.4490 -55019.0399
[250,] 111246.1796 -26527.4490
[251,] -39994.7608 111246.1796
[252,] -28722.6998 -39994.7608
[253,] -11443.7872 -28722.6998
[254,] -26268.4696 -11443.7872
[255,] 174727.8045 -26268.4696
[256,] 103716.9431 174727.8045
[257,] 107681.0092 103716.9431
[258,] -26524.0230 107681.0092
[259,] 74931.9518 -26524.0230
[260,] -55418.3976 74931.9518
[261,] -26926.8066 -55418.3976
[262,] 83049.3179 -26926.8066
[263,] 47905.3055 83049.3179
[264,] -85279.1680 47905.3055
[265,] 8992.1970 -85279.1680
[266,] -53551.1532 8992.1970
[267,] 176653.6049 -53551.1532
[268,] -65.2776 176653.6049
[269,] 44960.5327 -65.2776
[270,] -26923.3806 44960.5327
[271,] -5218.1900 -26923.3806
[272,] -55817.7552 -5218.1900
[273,] 50841.1017 -55817.7552
[274,] 103045.8317 50841.1017
[275,] 47762.1019 103045.8317
[276,] 654.0813 47762.1019
[277,] 8592.8394 654.0813
[278,] -66396.4955 8592.8394
[279,] 152290.5768 -66396.4955
[280,] -464.6352 152290.5768
[281,] 66628.5880 -464.6352
[282,] 90639.6794 66628.5880
[283,] 74133.2365 90639.6794
[284,] -12959.0918 74133.2365
[285,] 52025.2623 -12959.0918
[286,] 22895.6899 52025.2623
[287,] 37480.5295 22895.6899
[288,] -86077.8832 37480.5295
[289,] 38632.1684 -86077.8832
[290,] -54349.8684 38632.1684
[291,] 182986.2128 -54349.8684
[292,] 50920.9047 182986.2128
[293,] 66229.2303 50920.9047
[294,] 55605.5674 66229.2303
[295,] 28578.8406 55605.5674
[296,] 23134.3137 28578.8406
[297,] -28124.8796 23134.3137
[298,] 102247.1164 -28124.8796
[299,] 38158.5929 102247.1164
[300,] -87234.8524 38158.5929
[301,] -6156.9831 -87234.8524
[302,] 25001.5581 -6156.9831
[303,] 93379.5898 25001.5581
[304,] 116987.5958 93379.5898
[305,] 69851.0939 116987.5958
[306,] -28121.4535 69851.0939
[307,] -14683.7381 -28121.4535
[308,] 44669.5805 -14683.7381
[309,] -28524.2372 44669.5805
[310,] 22096.9746 -28524.2372
[311,] 15728.1348 22096.9746
[312,] -86876.5985 15728.1348
[313,] -75580.8618 -86876.5985
[314,] 26509.6285 -75580.8618
[315,] 173492.7836 26509.6285
[316,] 79709.6339 173492.7836
[317,] 51181.9059 79709.6339
[318,] 53935.6237 51181.9059
[319,] 72935.1636 53935.6237
[320,] 22335.5984 72935.1636
[321,] 50827.1894 22335.5984
[322,] 101448.4011 50827.1894
[323,] -42390.9066 101448.4011
[324,] -7525.1720 -42390.9066
[325,] 4409.5701 -7525.1720
[326,] -55547.9414 4409.5701
[327,] 92580.8745 -55547.9414
[328,] 77688.7184 92580.8745
[329,] 65031.1574 77688.7184
[330,] 50830.6154 65031.1574
[331,] 72535.8060 50830.6154
[332,] 21936.2408 72535.8060
[333,] 50430.3128 21936.2408
[334,] 101915.3553 50430.3128
[335,] 36960.5200 101915.3553
[336,] -9149.9613 36960.5200
[337,] 6596.0512 -9149.9613
[338,] 23138.1773 6596.0512
[339,] 175892.7380 23138.1773
[340,] 57987.9338 175892.7380
[341,] 64631.7998 57987.9338
[342,] 50431.2577 64631.7998
[343,] 72136.4484 50431.2577
[344,] 21536.8832 72136.4484
[345,] 50028.4741 21536.8832
[346,] 92659.4641 50028.4741
[347,] 36561.1623 92659.4641
[348,] -8323.8872 36561.1623
[349,] 6196.6935 -8323.8872
[350,] 24131.6419 6196.6935
[351,] 172759.1919 24131.6419
[352,] 78322.9168 172759.1919
[353,] 54266.9516 78322.9168
[354,] 50031.9001 54266.9516
[355,] 71737.0907 50031.9001
[356,] -13162.0903 71737.0907
[357,] 49629.1164 -13162.0903
[358,] 100250.3282 49629.1164
[359,] 36161.8047 100250.3282
[360,] -6936.7112 36161.8047
[361,] 5797.3359 -6936.7112
[362,] 23004.7699 5797.3359
[363,] 82172.0736 23004.7699
[364,] 75012.8225 82172.0736
[365,] 63833.0845 75012.8225
[366,] 49632.5425 63833.0845
[367,] 71337.7331 49632.5425
[368,] 20738.1679 71337.7331
[369,] 49229.7588 20738.1679
[370,] 99850.9706 49229.7588
[371,] 38492.1934 99850.9706
[372,] -49897.4353 38492.1934
[373,] -21610.5108 -49897.4353
[374,] 22605.4123 -21610.5108
[375,] 170734.2282 22605.4123
[376,] 88436.6202 170734.2282
[377,] 48186.2337 88436.6202
[378,] 49233.1848 48186.2337
[379,] 70938.3755 49233.1848
[380,] 20338.8102 70938.3755
[381,] 40858.6605 20338.8102
[382,] 152858.5696 40858.6605
[383,] 35363.0894 152858.5696
[384,] -7636.9399 35363.0894
[385,] -152062.9544 -7636.9399
[386,] 26878.5157 -152062.9544
[387,] 154494.8296 26878.5157
[388,] 73866.2132 154494.8296
[389,] 43678.0849 73866.2132
[390,] 45460.1314 43678.0849
[391,] 50020.7238 45460.1314
[392,] -17867.1035 50020.7238
[393,] 64608.0722 -17867.1035
[394,] -192929.4043 64608.0722
[395,] 7325.0672 -192929.4043
[396,] -173842.2880 7325.0672
[397,] -66961.9915 -173842.2880
[398,] -4978.1456 -66961.9915
[399,] 84161.6197 -4978.1456
[400,] 36945.7457 84161.6197
[401,] 52253.7733 36945.7457
[402,] 27767.4174 52253.7733
[403,] -92345.5438 27767.4174
[404,] -14294.8941 -92345.5438
[405,] -4643.4801 -14294.8941
[406,] 82458.6994 -4643.4801
[407,] 14013.2055 82458.6994
[408,] -32672.7658 14013.2055
[409,] -71588.7820 -32672.7658
[410,] -24342.8319 -71588.7820
[411,] -114359.6437 -24342.8319
[412,] -34348.3071 -114359.6437
[413,] -158884.3866 -34348.3071
[414,] -64220.2473 -158884.3866
[415,] -25464.5040 -64220.2473
[416,] 101931.1580 -25464.5040
[417,] -85882.5624 101931.1580
[418,] 46592.5195 -85882.5624
[419,] -34652.7350 46592.5195
[420,] -132968.7551 -34652.7350
[421,] -153473.7950 -132968.7551
[422,] -466415.9626 -153473.7950
[423,] -356872.4710 -466415.9626
[424,] -52907.5755 -356872.4710
[425,] -205270.7298 -52907.5755
[426,] -314466.3688 -205270.7298
[427,] -1195291.5387 -314466.3688
[428,] -397186.5168 -1195291.5387
[429,] -278836.6108 -397186.5168
[430,] -799679.6048 -278836.6108
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3219486.3610 1559382.1240
2 1316976.5266 3219486.3610
3 -3006346.2188 1316976.5266
4 -367191.0825 -3006346.2188
5 58124.8679 -367191.0825
6 2212.4809 58124.8679
7 613640.7565 2212.4809
8 -244097.0091 613640.7565
9 767206.3699 -244097.0091
10 -276783.0698 767206.3699
11 469034.9237 -276783.0698
12 912111.2466 469034.9237
13 -1813656.1773 912111.2466
14 17121.3969 -1813656.1773
15 561589.9066 17121.3969
16 -392968.6550 561589.9066
17 -64339.2419 -392968.6550
18 -147609.0330 -64339.2419
19 -307696.1383 -147609.0330
20 139641.1243 -307696.1383
21 -248941.0390 139641.1243
22 -187930.2372 -248941.0390
23 -502685.5525 -187930.2372
24 -273494.3746 -502685.5525
25 -165338.3774 -273494.3746
26 -397060.1855 -165338.3774
27 -314108.1510 -397060.1855
28 438513.2652 -314108.1510
29 40866.8211 438513.2652
30 174558.3890 40866.8211
31 -147551.9334 174558.3890
32 98489.0883 -147551.9334
33 9333.7365 98489.0883
34 -297071.3273 9333.7365
35 -465712.7664 -297071.3273
36 -503.6238 -465712.7664
37 21502.3355 -503.6238
38 -162914.9788 21502.3355
39 -26382.5030 -162914.9788
40 39197.1954 -26382.5030
41 -336941.6386 39197.1954
42 102970.6058 -336941.6386
43 -164200.1884 102970.6058
44 -37599.3789 -164200.1884
45 -240022.4042 -37599.3789
46 -211638.3635 -240022.4042
47 -198733.9311 -211638.3635
48 -162764.7183 -198733.9311
49 -198167.6920 -162764.7183
50 -10172.6875 -198167.6920
51 -15652.9742 -10172.6875
52 -391777.7316 -15652.9742
53 -47008.1210 -391777.7316
54 -407114.0059 -47008.1210
55 158694.2498 -407114.0059
56 157658.4218 158694.2498
57 -17093.0133 157658.4218
58 29662.1339 -17093.0133
59 -246820.5820 29662.1339
60 -154987.1448 -246820.5820
61 -98296.5473 -154987.1448
62 -143496.3063 -98296.5473
63 -117930.2557 -143496.3063
64 94141.3178 -117930.2557
65 85054.1162 94141.3178
66 -91122.9404 85054.1162
67 -86610.8068 -91122.9404
68 -89248.0610 -86610.8068
69 -145183.7168 -89248.0610
70 16918.9103 -145183.7168
71 63495.4408 16918.9103
72 -184187.0440 63495.4408
73 -57447.3363 -184187.0440
74 -149328.2290 -57447.3363
75 -27035.5362 -149328.2290
76 -66459.5557 -27035.5362
77 -31848.2961 -66459.5557
78 87945.1537 -31848.2961
79 -46679.2297 87945.1537
80 25210.4788 -46679.2297
81 -107039.8571 25210.4788
82 -71389.2716 -107039.8571
83 79001.3082 -71389.2716
84 -64219.7230 79001.3082
85 -162390.0116 -64219.7230
86 3461.1806 -162390.0116
87 81721.3766 3461.1806
88 -70002.3326 81721.3766
89 -6295.7068 -70002.3326
90 -26834.3556 -6295.7068
91 -89252.9649 -26834.3556
92 108382.8758 -89252.9649
93 -350163.4635 108382.8758
94 195354.9034 -350163.4635
95 -42757.0653 195354.9034
96 -47277.4307 -42757.0653
97 -6488.8067 -47277.4307
98 -161515.9341 -6488.8067
99 27978.6061 -161515.9341
100 -75472.9421 27978.6061
101 -87156.7202 -75472.9421
102 -112159.6626 -87156.7202
103 -33603.8601 -112159.6626
104 31825.0501 -33603.8601
105 -43708.3939 31825.0501
106 -196431.7739 -43708.3939
107 -24151.3414 -196431.7739
108 -262899.4675 -24151.3414
109 -204462.3356 -262899.4675
110 90212.0722 -204462.3356
111 -120314.5081 90212.0722
112 -241889.4258 -120314.5081
113 -168288.3944 -241889.4258
114 -15188.6548 -168288.3944
115 102515.1909 -15188.6548
116 -17667.6225 102515.1909
117 40424.0791 -17667.6225
118 -195303.7037 40424.0791
119 64457.7881 -195303.7037
120 -320076.0226 64457.7881
121 775.4535 -320076.0226
122 -26436.4370 775.4535
123 62494.9877 -26436.4370
124 -102103.9437 62494.9877
125 -244215.1272 -102103.9437
126 57619.6952 -244215.1272
127 79535.7074 57619.6952
128 36760.3084 79535.7074
129 57216.9115 36760.3084
130 110945.1043 57216.9115
131 101007.1573 110945.1043
132 -1135.4498 101007.1573
133 11317.4683 -1135.4498
134 30592.5650 11317.4683
135 181231.2433 30592.5650
136 86246.7689 181231.2433
137 71420.8796 86246.7689
138 57220.3376 71420.8796
139 78925.5282 57220.3376
140 27326.2681 78925.5282
141 56817.5539 27326.2681
142 102382.0307 56817.5539
143 43576.6836 102382.0307
144 4460.6000 43576.6836
145 25076.9885 4460.6000
146 30193.2074 25076.9885
147 180204.3464 30193.2074
148 83679.0830 180204.3464
149 70427.3078 83679.0830
150 47520.2898 70427.3078
151 67175.5227 47520.2898
152 27926.6053 67175.5227
153 34196.0012 27926.6053
154 107085.0498 34196.0012
155 42950.8845 107085.0498
156 -1934.1650 42950.8845
157 9361.9570 -1934.1650
158 36912.3878 9361.9570
159 177922.6656 36912.3878
160 83279.7253 177922.6656
161 70622.1643 83279.7253
162 56421.6223 70622.1643
163 83601.2759 56421.6223
164 27524.5756 83601.2759
165 56018.8386 27524.5756
166 106640.0504 56018.8386
167 146863.4910 106640.0504
168 -2333.5227 146863.4910
169 12187.0581 -2333.5227
170 29394.4921 12187.0581
171 174437.8893 29394.4921
172 82880.3677 174437.8893
173 67823.9325 82880.3677
174 56022.2646 67823.9325
175 77666.3580 56022.2646
176 27127.8901 77666.3580
177 -25384.8721 27127.8901
178 106240.6928 -25384.8721
179 42152.1693 106240.6928
180 10572.2577 42152.1693
181 -91992.7523 10572.2577
182 54567.9171 -91992.7523
183 34750.1699 54567.9171
184 1860.2953 34750.1699
185 -17583.4095 1860.2953
186 -19283.0027 -17583.4095
187 73886.6124 -19283.0027
188 -13283.6223 73886.6124
189 55220.1234 -13283.6223
190 110553.7260 55220.1234
191 41752.8116 110553.7260
192 -1075.2154 41752.8116
193 11388.3428 -1075.2154
194 28291.4501 11388.3428
195 96973.8085 28291.4501
196 -137338.0361 96973.8085
197 -10326.6928 -137338.0361
198 -24527.2348 -10326.6928
199 -82017.9469 -24527.2348
200 -53421.6094 -82017.9469
201 -24930.0185 -53421.6094
202 25691.1933 -24930.0185
203 -9968.9404 25691.1933
204 -171156.9114 -9968.9404
205 -70568.1759 -171156.9114
206 -67041.3728 -70568.1759
207 96574.4509 -67041.3728
208 61264.8860 96574.4509
209 -10726.0504 61264.8860
210 54824.1917 -10726.0504
211 76529.3824 54824.1917
212 41602.3143 76529.3824
213 31945.2229 41602.3143
214 86180.3762 31945.2229
215 40954.0963 86180.3762
216 -481.6499 40954.0963
217 14565.1699 -481.6499
218 27797.0616 14565.1699
219 180964.6773 27797.0616
220 48481.7709 180964.6773
221 68625.3761 48481.7709
222 55310.6017 68625.3761
223 76130.0247 55310.6017
224 25815.3369 76130.0247
225 58518.9376 25815.3369
226 104643.2622 58518.9376
227 40554.7387 104643.2622
228 -30779.0600 40554.7387
229 10190.2699 -30779.0600
230 6028.6918 10190.2699
231 175526.5198 6028.6918
232 80883.5795 175526.5198
233 69429.7276 80883.5795
234 54025.4765 69429.7276
235 75977.6294 54025.4765
236 25131.1019 75977.6294
237 56557.0814 25131.1019
238 104243.9046 56557.0814
239 40155.3811 104243.9046
240 -4729.6685 40155.3811
241 12619.6182 -4729.6685
242 26998.3463 12619.6182
243 180305.4857 26998.3463
244 -7034.8025 180305.4857
245 -11924.1233 -7034.8025
246 -26124.6654 -11924.1233
247 67052.9749 -26124.6654
248 -55019.0399 67052.9749
249 -26527.4490 -55019.0399
250 111246.1796 -26527.4490
251 -39994.7608 111246.1796
252 -28722.6998 -39994.7608
253 -11443.7872 -28722.6998
254 -26268.4696 -11443.7872
255 174727.8045 -26268.4696
256 103716.9431 174727.8045
257 107681.0092 103716.9431
258 -26524.0230 107681.0092
259 74931.9518 -26524.0230
260 -55418.3976 74931.9518
261 -26926.8066 -55418.3976
262 83049.3179 -26926.8066
263 47905.3055 83049.3179
264 -85279.1680 47905.3055
265 8992.1970 -85279.1680
266 -53551.1532 8992.1970
267 176653.6049 -53551.1532
268 -65.2776 176653.6049
269 44960.5327 -65.2776
270 -26923.3806 44960.5327
271 -5218.1900 -26923.3806
272 -55817.7552 -5218.1900
273 50841.1017 -55817.7552
274 103045.8317 50841.1017
275 47762.1019 103045.8317
276 654.0813 47762.1019
277 8592.8394 654.0813
278 -66396.4955 8592.8394
279 152290.5768 -66396.4955
280 -464.6352 152290.5768
281 66628.5880 -464.6352
282 90639.6794 66628.5880
283 74133.2365 90639.6794
284 -12959.0918 74133.2365
285 52025.2623 -12959.0918
286 22895.6899 52025.2623
287 37480.5295 22895.6899
288 -86077.8832 37480.5295
289 38632.1684 -86077.8832
290 -54349.8684 38632.1684
291 182986.2128 -54349.8684
292 50920.9047 182986.2128
293 66229.2303 50920.9047
294 55605.5674 66229.2303
295 28578.8406 55605.5674
296 23134.3137 28578.8406
297 -28124.8796 23134.3137
298 102247.1164 -28124.8796
299 38158.5929 102247.1164
300 -87234.8524 38158.5929
301 -6156.9831 -87234.8524
302 25001.5581 -6156.9831
303 93379.5898 25001.5581
304 116987.5958 93379.5898
305 69851.0939 116987.5958
306 -28121.4535 69851.0939
307 -14683.7381 -28121.4535
308 44669.5805 -14683.7381
309 -28524.2372 44669.5805
310 22096.9746 -28524.2372
311 15728.1348 22096.9746
312 -86876.5985 15728.1348
313 -75580.8618 -86876.5985
314 26509.6285 -75580.8618
315 173492.7836 26509.6285
316 79709.6339 173492.7836
317 51181.9059 79709.6339
318 53935.6237 51181.9059
319 72935.1636 53935.6237
320 22335.5984 72935.1636
321 50827.1894 22335.5984
322 101448.4011 50827.1894
323 -42390.9066 101448.4011
324 -7525.1720 -42390.9066
325 4409.5701 -7525.1720
326 -55547.9414 4409.5701
327 92580.8745 -55547.9414
328 77688.7184 92580.8745
329 65031.1574 77688.7184
330 50830.6154 65031.1574
331 72535.8060 50830.6154
332 21936.2408 72535.8060
333 50430.3128 21936.2408
334 101915.3553 50430.3128
335 36960.5200 101915.3553
336 -9149.9613 36960.5200
337 6596.0512 -9149.9613
338 23138.1773 6596.0512
339 175892.7380 23138.1773
340 57987.9338 175892.7380
341 64631.7998 57987.9338
342 50431.2577 64631.7998
343 72136.4484 50431.2577
344 21536.8832 72136.4484
345 50028.4741 21536.8832
346 92659.4641 50028.4741
347 36561.1623 92659.4641
348 -8323.8872 36561.1623
349 6196.6935 -8323.8872
350 24131.6419 6196.6935
351 172759.1919 24131.6419
352 78322.9168 172759.1919
353 54266.9516 78322.9168
354 50031.9001 54266.9516
355 71737.0907 50031.9001
356 -13162.0903 71737.0907
357 49629.1164 -13162.0903
358 100250.3282 49629.1164
359 36161.8047 100250.3282
360 -6936.7112 36161.8047
361 5797.3359 -6936.7112
362 23004.7699 5797.3359
363 82172.0736 23004.7699
364 75012.8225 82172.0736
365 63833.0845 75012.8225
366 49632.5425 63833.0845
367 71337.7331 49632.5425
368 20738.1679 71337.7331
369 49229.7588 20738.1679
370 99850.9706 49229.7588
371 38492.1934 99850.9706
372 -49897.4353 38492.1934
373 -21610.5108 -49897.4353
374 22605.4123 -21610.5108
375 170734.2282 22605.4123
376 88436.6202 170734.2282
377 48186.2337 88436.6202
378 49233.1848 48186.2337
379 70938.3755 49233.1848
380 20338.8102 70938.3755
381 40858.6605 20338.8102
382 152858.5696 40858.6605
383 35363.0894 152858.5696
384 -7636.9399 35363.0894
385 -152062.9544 -7636.9399
386 26878.5157 -152062.9544
387 154494.8296 26878.5157
388 73866.2132 154494.8296
389 43678.0849 73866.2132
390 45460.1314 43678.0849
391 50020.7238 45460.1314
392 -17867.1035 50020.7238
393 64608.0722 -17867.1035
394 -192929.4043 64608.0722
395 7325.0672 -192929.4043
396 -173842.2880 7325.0672
397 -66961.9915 -173842.2880
398 -4978.1456 -66961.9915
399 84161.6197 -4978.1456
400 36945.7457 84161.6197
401 52253.7733 36945.7457
402 27767.4174 52253.7733
403 -92345.5438 27767.4174
404 -14294.8941 -92345.5438
405 -4643.4801 -14294.8941
406 82458.6994 -4643.4801
407 14013.2055 82458.6994
408 -32672.7658 14013.2055
409 -71588.7820 -32672.7658
410 -24342.8319 -71588.7820
411 -114359.6437 -24342.8319
412 -34348.3071 -114359.6437
413 -158884.3866 -34348.3071
414 -64220.2473 -158884.3866
415 -25464.5040 -64220.2473
416 101931.1580 -25464.5040
417 -85882.5624 101931.1580
418 46592.5195 -85882.5624
419 -34652.7350 46592.5195
420 -132968.7551 -34652.7350
421 -153473.7950 -132968.7551
422 -466415.9626 -153473.7950
423 -356872.4710 -466415.9626
424 -52907.5755 -356872.4710
425 -205270.7298 -52907.5755
426 -314466.3688 -205270.7298
427 -1195291.5387 -314466.3688
428 -397186.5168 -1195291.5387
429 -278836.6108 -397186.5168
430 -799679.6048 -278836.6108
> 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/73lkm1291311351.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/83lkm1291311351.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/9wuk71291311351.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/10wuk71291311351.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/11zdiv1291311351.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/123vzj1291311351.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/139eed1291311351.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/142ody1291311351.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/15not41291311351.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/161g9u1291311351.tab")
+ }
>
> try(system("convert tmp/1ptne1291311351.ps tmp/1ptne1291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/2il4z1291311351.ps tmp/2il4z1291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/3il4z1291311351.ps tmp/3il4z1291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/4il4z1291311351.ps tmp/4il4z1291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/5su321291311351.ps tmp/5su321291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/6su321291311351.ps tmp/6su321291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/73lkm1291311351.ps tmp/73lkm1291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/83lkm1291311351.ps tmp/83lkm1291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/9wuk71291311351.ps tmp/9wuk71291311351.png",intern=TRUE))
character(0)
> try(system("convert tmp/10wuk71291311351.ps tmp/10wuk71291311351.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
16.210 3.221 16.583