R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(6282929
+ ,1
+ ,162556
+ ,1081
+ ,213118
+ ,4324047
+ ,1
+ ,29790
+ ,309
+ ,81767
+ ,4108272
+ ,1
+ ,87550
+ ,458
+ ,153198
+ ,-1212617
+ ,0
+ ,84738
+ ,588
+ ,-26007
+ ,1485329
+ ,1
+ ,54660
+ ,299
+ ,126942
+ ,1779876
+ ,1
+ ,42634
+ ,156
+ ,157214
+ ,1367203
+ ,0
+ ,40949
+ ,481
+ ,129352
+ ,2519076
+ ,1
+ ,42312
+ ,323
+ ,234817
+ ,912684
+ ,1
+ ,37704
+ ,452
+ ,60448
+ ,1443586
+ ,1
+ ,16275
+ ,109
+ ,47818
+ ,1220017
+ ,0
+ ,25830
+ ,115
+ ,245546
+ ,984885
+ ,0
+ ,12679
+ ,110
+ ,48020
+ ,1457425
+ ,1
+ ,18014
+ ,239
+ ,-1710
+ ,-572920
+ ,0
+ ,43556
+ ,247
+ ,32648
+ ,929144
+ ,1
+ ,24524
+ ,497
+ ,95350
+ ,1151176
+ ,0
+ ,6532
+ ,103
+ ,151352
+ ,790090
+ ,0
+ ,7123
+ ,109
+ ,288170
+ ,774497
+ ,1
+ ,20813
+ ,502
+ ,114337
+ ,990576
+ ,1
+ ,37597
+ ,248
+ ,37884
+ ,454195
+ ,0
+ ,17821
+ ,373
+ ,122844
+ ,876607
+ ,1
+ ,12988
+ ,119
+ ,82340
+ ,711969
+ ,1
+ ,22330
+ ,84
+ ,79801
+ ,702380
+ ,0
+ ,13326
+ ,102
+ ,165548
+ ,264449
+ ,0
+ ,16189
+ ,295
+ ,116384
+ ,450033
+ ,0
+ ,7146
+ ,105
+ ,134028
+ ,541063
+ ,0
+ ,15824
+ ,64
+ ,63838
+ ,588864
+ ,1
+ ,26088
+ ,267
+ ,74996
+ ,-37216
+ ,0
+ ,11326
+ ,129
+ ,31080
+ ,783310
+ ,0
+ ,8568
+ ,37
+ ,32168
+ ,467359
+ ,0
+ ,14416
+ ,361
+ ,49857
+ ,688779
+ ,1
+ ,3369
+ ,28
+ ,87161
+ ,608419
+ ,1
+ ,11819
+ ,85
+ ,106113
+ ,696348
+ ,1
+ ,6620
+ ,44
+ ,80570
+ ,597793
+ ,1
+ ,4519
+ ,49
+ ,102129
+ ,821730
+ ,0
+ ,2220
+ ,22
+ ,301670
+ ,377934
+ ,0
+ ,18562
+ ,155
+ ,102313
+ ,651939
+ ,0
+ ,10327
+ ,91
+ ,88577
+ ,697458
+ ,1
+ ,5336
+ ,81
+ ,112477
+ ,700368
+ ,1
+ ,2365
+ ,79
+ ,191778
+ ,225986
+ ,0
+ ,4069
+ ,145
+ ,79804
+ ,348695
+ ,0
+ ,7710
+ ,816
+ ,128294
+ ,373683
+ ,0
+ ,13718
+ ,61
+ ,96448
+ ,501709
+ ,0
+ ,4525
+ ,226
+ ,93811
+ ,413743
+ ,0
+ ,6869
+ ,105
+ ,117520
+ ,379825
+ ,0
+ ,4628
+ ,62
+ ,69159
+ ,336260
+ ,1
+ ,3653
+ ,24
+ ,101792
+ ,636765
+ ,1
+ ,1265
+ ,26
+ ,210568
+ ,481231
+ ,1
+ ,7489
+ ,322
+ ,136996
+ ,469107
+ ,0
+ ,4901
+ ,84
+ ,121920
+ ,211928
+ ,0
+ ,2284
+ ,33
+ ,76403
+ ,563925
+ ,1
+ ,3160
+ ,108
+ ,108094
+ ,511939
+ ,1
+ ,4150
+ ,150
+ ,134759
+ ,521016
+ ,1
+ ,7285
+ ,115
+ ,188873
+ ,543856
+ ,1
+ ,1134
+ ,162
+ ,146216
+ ,329304
+ ,1
+ ,4658
+ ,158
+ ,156608
+ ,423262
+ ,0
+ ,2384
+ ,97
+ ,61348
+ ,509665
+ ,0
+ ,3748
+ ,9
+ ,50350
+ ,455881
+ ,0
+ ,5371
+ ,66
+ ,87720
+ ,367772
+ ,0
+ ,1285
+ ,107
+ ,99489
+ ,406339
+ ,1
+ ,9327
+ ,101
+ ,87419
+ ,493408
+ ,1
+ ,5565
+ ,47
+ ,94355
+ ,232942
+ ,0
+ ,1528
+ ,38
+ ,60326
+ ,416002
+ ,1
+ ,3122
+ ,34
+ ,94670
+ ,337430
+ ,1
+ ,7317
+ ,84
+ ,82425
+ ,361517
+ ,0
+ ,2675
+ ,79
+ ,59017
+ ,360962
+ ,0
+ ,13253
+ ,947
+ ,90829
+ ,235561
+ ,0
+ ,880
+ ,74
+ ,80791
+ ,408247
+ ,1
+ ,2053
+ ,53
+ ,100423
+ ,450296
+ ,0
+ ,1424
+ ,94
+ ,131116
+ ,418799
+ ,1
+ ,4036
+ ,63
+ ,100269
+ ,247405
+ ,1
+ ,3045
+ ,58
+ ,27330
+ ,378519
+ ,0
+ ,5119
+ ,49
+ ,39039
+ ,326638
+ ,0
+ ,1431
+ ,34
+ ,106885
+ ,328233
+ ,0
+ ,554
+ ,11
+ ,79285
+ ,386225
+ ,0
+ ,1975
+ ,35
+ ,118881
+ ,283662
+ ,1
+ ,1286
+ ,17
+ ,77623
+ ,370225
+ ,0
+ ,1012
+ ,47
+ ,114768
+ ,269236
+ ,0
+ ,810
+ ,43
+ ,74015
+ ,365732
+ ,0
+ ,1280
+ ,117
+ ,69465
+ ,420383
+ ,1
+ ,666
+ ,171
+ ,117869
+ ,345811
+ ,0
+ ,1380
+ ,26
+ ,60982
+ ,431809
+ ,1
+ ,4608
+ ,73
+ ,90131
+ ,418876
+ ,0
+ ,876
+ ,59
+ ,138971
+ ,297476
+ ,0
+ ,814
+ ,18
+ ,39625
+ ,416776
+ ,0
+ ,514
+ ,15
+ ,102725
+ ,357257
+ ,1
+ ,5692
+ ,72
+ ,64239
+ ,458343
+ ,0
+ ,3642
+ ,86
+ ,90262
+ ,388386
+ ,0
+ ,540
+ ,14
+ ,103960
+ ,358934
+ ,0
+ ,2099
+ ,64
+ ,106611
+ ,407560
+ ,0
+ ,567
+ ,11
+ ,103345
+ ,392558
+ ,0
+ ,2001
+ ,52
+ ,95551
+ ,373177
+ ,1
+ ,2949
+ ,41
+ ,82903
+ ,428370
+ ,0
+ ,2253
+ ,99
+ ,63593
+ ,369419
+ ,1
+ ,6533
+ ,75
+ ,126910
+ ,358649
+ ,0
+ ,1889
+ ,45
+ ,37527
+ ,376641
+ ,1
+ ,3055
+ ,43
+ ,60247
+ ,467427
+ ,0
+ ,272
+ ,8
+ ,112995
+ ,364885
+ ,1
+ ,1414
+ ,198
+ ,70184
+ ,436230
+ ,0
+ ,2564
+ ,22
+ ,130140
+ ,329118
+ ,1
+ ,1383
+ ,11
+ ,73221
+ ,317365
+ ,1
+ ,1261
+ ,33
+ ,76114
+ ,286849
+ ,0
+ ,975
+ ,23
+ ,90534
+ ,376685
+ ,0
+ ,3366
+ ,80
+ ,108479
+ ,407198
+ ,0
+ ,576
+ ,18
+ ,113761
+ ,377772
+ ,0
+ ,1306
+ ,28
+ ,68696
+ ,271483
+ ,0
+ ,746
+ ,23
+ ,71561
+ ,153661
+ ,1
+ ,3192
+ ,60
+ ,59831
+ ,513294
+ ,1
+ ,2045
+ ,20
+ ,97890
+ ,324881
+ ,0
+ ,5477
+ ,59
+ ,101481
+ ,264512
+ ,1
+ ,1932
+ ,36
+ ,72954
+ ,420968
+ ,0
+ ,936
+ ,30
+ ,67939
+ ,129302
+ ,1
+ ,3437
+ ,47
+ ,48022
+ ,191521
+ ,0
+ ,5131
+ ,71
+ ,86111
+ ,268673
+ ,1
+ ,2397
+ ,14
+ ,74020
+ ,353179
+ ,1
+ ,1389
+ ,9
+ ,57530
+ ,354624
+ ,0
+ ,1503
+ ,39
+ ,56364
+ ,363713
+ ,0
+ ,402
+ ,26
+ ,84990
+ ,456657
+ ,0
+ ,2239
+ ,21
+ ,88590
+ ,211742
+ ,1
+ ,2234
+ ,16
+ ,77200
+ ,338381
+ ,0
+ ,837
+ ,69
+ ,61262
+ ,418530
+ ,0
+ ,10579
+ ,92
+ ,110309
+ ,351483
+ ,0
+ ,875
+ ,14
+ ,67000
+ ,372928
+ ,0
+ ,1395
+ ,103
+ ,93099
+ ,485538
+ ,1
+ ,1659
+ ,29
+ ,107577
+ ,279268
+ ,1
+ ,2647
+ ,37
+ ,62920
+ ,219060
+ ,1
+ ,3294
+ ,23
+ ,75832
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325314
+ ,0
+ ,94
+ ,7
+ ,60793
+ ,322046
+ ,0
+ ,422
+ ,28
+ ,57935
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325599
+ ,0
+ ,34
+ ,8
+ ,60630
+ ,377028
+ ,0
+ ,1558
+ ,63
+ ,55637
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,323850
+ ,0
+ ,43
+ ,3
+ ,60887
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331514
+ ,0
+ ,316
+ ,9
+ ,60505
+ ,325632
+ ,0
+ ,115
+ ,13
+ ,60945
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,322265
+ ,0
+ ,389
+ ,14
+ ,58990
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325906
+ ,0
+ ,1002
+ ,15
+ ,56750
+ ,325985
+ ,0
+ ,36
+ ,3
+ ,60894
+ ,346145
+ ,0
+ ,460
+ ,15
+ ,63346
+ ,325898
+ ,0
+ ,309
+ ,11
+ ,56535
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325356
+ ,0
+ ,9
+ ,6
+ ,60835
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325930
+ ,0
+ ,14
+ ,1
+ ,61016
+ ,318020
+ ,0
+ ,520
+ ,10
+ ,58650
+ ,326389
+ ,0
+ ,1766
+ ,73
+ ,60438
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,302925
+ ,0
+ ,458
+ ,11
+ ,58625
+ ,325540
+ ,0
+ ,20
+ ,3
+ ,60938
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,326736
+ ,0
+ ,98
+ ,2
+ ,61490
+ ,340580
+ ,0
+ ,405
+ ,7
+ ,60845
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331828
+ ,0
+ ,483
+ ,27
+ ,60830
+ ,323299
+ ,0
+ ,454
+ ,51
+ ,63261
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,387722
+ ,0
+ ,757
+ ,19
+ ,45689
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324598
+ ,0
+ ,36
+ ,4
+ ,61564
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,328726
+ ,0
+ ,203
+ ,9
+ ,61938
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325043
+ ,0
+ ,90
+ ,8
+ ,60951
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325806
+ ,1
+ ,71
+ ,1
+ ,60745
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,387732
+ ,0
+ ,972
+ ,34
+ ,71642
+ ,349729
+ ,1
+ ,531
+ ,10
+ ,71641
+ ,332202
+ ,0
+ ,604
+ ,38
+ ,55792
+ ,305442
+ ,1
+ ,283
+ ,10
+ ,71873
+ ,329537
+ ,1
+ ,23
+ ,5
+ ,62555
+ ,327055
+ ,1
+ ,638
+ ,14
+ ,60370
+ ,356245
+ ,1
+ ,699
+ ,16
+ ,64873
+ ,328451
+ ,0
+ ,149
+ ,5
+ ,62041
+ ,307062
+ ,0
+ ,226
+ ,5
+ ,65745
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331345
+ ,0
+ ,275
+ ,4
+ ,59500
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331824
+ ,0
+ ,141
+ ,6
+ ,61630
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325685
+ ,0
+ ,28
+ ,2
+ ,60890
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,404480
+ ,1
+ ,2566
+ ,80
+ ,113521
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,318314
+ ,1
+ ,472
+ ,20
+ ,80045
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,311807
+ ,1
+ ,203
+ ,27
+ ,50804
+ ,337724
+ ,1
+ ,496
+ ,17
+ ,87390
+ ,326431
+ ,1
+ ,10
+ ,2
+ ,61656
+ ,327556
+ ,1
+ ,63
+ ,4
+ ,65688
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,356850
+ ,1
+ ,1136
+ ,32
+ ,48522
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,322741
+ ,0
+ ,267
+ ,32
+ ,57640
+ ,310902
+ ,0
+ ,474
+ ,20
+ ,61977
+ ,324295
+ ,0
+ ,534
+ ,7
+ ,62620
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,326156
+ ,0
+ ,15
+ ,8
+ ,60831
+ ,326960
+ ,0
+ ,397
+ ,28
+ ,60646
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,333411
+ ,0
+ ,1061
+ ,20
+ ,56225
+ ,297761
+ ,0
+ ,288
+ ,4
+ ,60510
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325536
+ ,0
+ ,3
+ ,2
+ ,60698
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325762
+ ,0
+ ,20
+ ,2
+ ,60805
+ ,327957
+ ,0
+ ,278
+ ,26
+ ,61404
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,318521
+ ,0
+ ,192
+ ,4
+ ,65276
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,319775
+ ,0
+ ,317
+ ,9
+ ,63915
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,332128
+ ,0
+ ,368
+ ,17
+ ,61686
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325486
+ ,0
+ ,2
+ ,1
+ ,60743
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325838
+ ,0
+ ,53
+ ,6
+ ,60349
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,331767
+ ,0
+ ,94
+ ,3
+ ,61360
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324523
+ ,0
+ ,24
+ ,8
+ ,59818
+ ,339995
+ ,0
+ ,2332
+ ,4
+ ,72680
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,319582
+ ,0
+ ,131
+ ,11
+ ,61808
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,307245
+ ,0
+ ,206
+ ,9
+ ,53110
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,317967
+ ,0
+ ,167
+ ,2
+ ,64245
+ ,331488
+ ,0
+ ,622
+ ,73
+ ,73007
+ ,335452
+ ,0
+ ,885
+ ,85
+ ,82732
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,334184
+ ,0
+ ,365
+ ,8
+ ,54820
+ ,313213
+ ,0
+ ,364
+ ,35
+ ,47705
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,348678
+ ,0
+ ,226
+ ,12
+ ,72835
+ ,328727
+ ,0
+ ,307
+ ,15
+ ,58856
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,387978
+ ,0
+ ,188
+ ,11
+ ,77655
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,336704
+ ,0
+ ,138
+ ,6
+ ,69817
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,322076
+ ,0
+ ,125
+ ,12
+ ,60798
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,334272
+ ,0
+ ,282
+ ,30
+ ,62452
+ ,338197
+ ,0
+ ,335
+ ,33
+ ,64175
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,321024
+ ,1
+ ,813
+ ,82
+ ,67440
+ ,322145
+ ,0
+ ,176
+ ,28
+ ,68136
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,323351
+ ,0
+ ,249
+ ,72
+ ,56726
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,327748
+ ,0
+ ,333
+ ,13
+ ,70811
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,328157
+ ,0
+ ,30
+ ,4
+ ,62045
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,311594
+ ,0
+ ,249
+ ,62
+ ,54323
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,335962
+ ,0
+ ,165
+ ,24
+ ,62841
+ ,372426
+ ,0
+ ,453
+ ,21
+ ,81125
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,319844
+ ,0
+ ,53
+ ,14
+ ,59506
+ ,355822
+ ,1
+ ,382
+ ,21
+ ,59365
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324047
+ ,1
+ ,30
+ ,4
+ ,60798
+ ,311464
+ ,0
+ ,290
+ ,2
+ ,58790
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,353417
+ ,0
+ ,366
+ ,53
+ ,61808
+ ,325590
+ ,0
+ ,2
+ ,9
+ ,60735
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,328576
+ ,1
+ ,209
+ ,13
+ ,64016
+ ,326126
+ ,0
+ ,384
+ ,22
+ ,54683
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,369376
+ ,0
+ ,365
+ ,83
+ ,87192
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,332013
+ ,1
+ ,49
+ ,8
+ ,64107
+ ,325871
+ ,0
+ ,3
+ ,4
+ ,60761
+ ,342165
+ ,0
+ ,133
+ ,14
+ ,65990
+ ,324967
+ ,0
+ ,32
+ ,1
+ ,59988
+ ,314832
+ ,0
+ ,368
+ ,17
+ ,61167
+ ,325557
+ ,0
+ ,1
+ ,6
+ ,60719
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,322649
+ ,0
+ ,22
+ ,2
+ ,60722
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,1
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324598
+ ,0
+ ,96
+ ,5
+ ,60379
+ ,325567
+ ,0
+ ,1
+ ,2
+ ,60727
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,324005
+ ,0
+ ,81
+ ,7
+ ,60925
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325748
+ ,0
+ ,26
+ ,1
+ ,60896
+ ,323385
+ ,0
+ ,125
+ ,13
+ ,59734
+ ,315409
+ ,0
+ ,304
+ ,15
+ ,62969
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,312275
+ ,0
+ ,119
+ ,6
+ ,59118
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,320576
+ ,0
+ ,312
+ ,14
+ ,58598
+ ,325246
+ ,0
+ ,60
+ ,10
+ ,61124
+ ,332961
+ ,0
+ ,587
+ ,12
+ ,59595
+ ,323010
+ ,0
+ ,135
+ ,2
+ ,62065
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,345253
+ ,0
+ ,514
+ ,52
+ ,78780
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325559
+ ,0
+ ,1
+ ,4
+ ,60722
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,319634
+ ,1
+ ,58
+ ,3
+ ,61600
+ ,319951
+ ,0
+ ,180
+ ,11
+ ,59635
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,318519
+ ,0
+ ,448
+ ,40
+ ,59781
+ ,343222
+ ,0
+ ,227
+ ,9
+ ,76644
+ ,317234
+ ,0
+ ,174
+ ,1
+ ,64820
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,314025
+ ,0
+ ,121
+ ,24
+ ,56178
+ ,320249
+ ,0
+ ,607
+ ,11
+ ,60436
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,349365
+ ,0
+ ,530
+ ,60
+ ,73433
+ ,289197
+ ,0
+ ,571
+ ,80
+ ,41477
+ ,325560
+ ,0
+ ,0
+ ,0
+ ,60720
+ ,329245
+ ,0
+ ,78
+ ,16
+ ,62700
+ ,240869
+ ,0
+ ,2489
+ ,40
+ ,67804
+ ,327182
+ ,0
+ ,131
+ ,6
+ ,59661
+ ,322876
+ ,0
+ ,923
+ ,8
+ ,58620
+ ,323117
+ ,0
+ ,72
+ ,3
+ ,60398
+ ,306351
+ ,0
+ ,572
+ ,16
+ ,58580
+ ,335137
+ ,0
+ ,397
+ ,10
+ ,62710
+ ,308271
+ ,0
+ ,450
+ ,8
+ ,59325
+ ,301731
+ ,0
+ ,622
+ ,7
+ ,60950
+ ,382409
+ ,0
+ ,694
+ ,8
+ ,68060
+ ,279230
+ ,1
+ ,3425
+ ,12
+ ,83620
+ ,298731
+ ,0
+ ,562
+ ,13
+ ,58456
+ ,243650
+ ,0
+ ,4917
+ ,42
+ ,52811
+ ,532682
+ ,1
+ ,1442
+ ,118
+ ,121173
+ ,319771
+ ,0
+ ,529
+ ,9
+ ,63870
+ ,171493
+ ,1
+ ,2126
+ ,138
+ ,21001
+ ,347262
+ ,0
+ ,1061
+ ,5
+ ,70415
+ ,343945
+ ,0
+ ,776
+ ,9
+ ,64230
+ ,311874
+ ,0
+ ,611
+ ,8
+ ,59190
+ ,302211
+ ,1
+ ,1526
+ ,25
+ ,69351
+ ,316708
+ ,0
+ ,592
+ ,7
+ ,64270
+ ,333463
+ ,0
+ ,1182
+ ,13
+ ,70694
+ ,344282
+ ,0
+ ,621
+ ,16
+ ,68005
+ ,319635
+ ,0
+ ,989
+ ,11
+ ,58930
+ ,301186
+ ,0
+ ,438
+ ,11
+ ,58320
+ ,300381
+ ,0
+ ,726
+ ,3
+ ,69980
+ ,318765
+ ,0
+ ,1303
+ ,61
+ ,69863
+ ,286146
+ ,1
+ ,6341
+ ,24
+ ,63255
+ ,306844
+ ,1
+ ,1164
+ ,17
+ ,57320
+ ,307705
+ ,1
+ ,3310
+ ,33
+ ,75230
+ ,312448
+ ,0
+ ,1366
+ ,7
+ ,79420
+ ,299715
+ ,0
+ ,965
+ ,3
+ ,73490
+ ,373399
+ ,0
+ ,3256
+ ,66
+ ,35250
+ ,299446
+ ,1
+ ,1135
+ ,17
+ ,62285
+ ,325586
+ ,0
+ ,1270
+ ,26
+ ,69206
+ ,291221
+ ,0
+ ,661
+ ,3
+ ,65920
+ ,261173
+ ,0
+ ,1013
+ ,2
+ ,69770
+ ,255027
+ ,0
+ ,2844
+ ,67
+ ,72683
+ ,-78375
+ ,1
+ ,11528
+ ,70
+ ,-14545
+ ,-58143
+ ,0
+ ,6526
+ ,26
+ ,55830
+ ,227033
+ ,0
+ ,2264
+ ,24
+ ,55174
+ ,235098
+ ,1
+ ,4461
+ ,94
+ ,67038
+ ,21267
+ ,0
+ ,3999
+ ,30
+ ,51252
+ ,238675
+ ,0
+ ,35624
+ ,223
+ ,157278
+ ,197687
+ ,0
+ ,9252
+ ,48
+ ,79510
+ ,418341
+ ,0
+ ,15236
+ ,90
+ ,77440
+ ,-297706
+ ,0
+ ,18073
+ ,180
+ ,27284)
+ ,dim=c(5
+ ,431)
+ ,dimnames=list(c('Wealth'
+ ,'Group'
+ ,'Costs'
+ ,'Trades'
+ ,'Dividends')
+ ,1:431))
> y <- array(NA,dim=c(5,431),dimnames=list(c('Wealth','Group','Costs','Trades','Dividends'),1:431))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = '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 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
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Group Costs Trades Dividends M1
7.218e+04 7.870e+04 2.628e+01 -4.662e+02 3.513e+00 4.531e+04
M2 M3 M4 M5 M6 M7
3.026e+04 1.364e+04 -1.349e+05 -3.977e+04 -2.718e+04 -1.341e+04
M8 M9 M10 M11
-3.472e+04 1.513e+04 -1.311e+04 -6.342e+04
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-3011398 -45709 24947 67249 3216805
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 7.218e+04 5.942e+04 1.215 0.2252
Group 7.870e+04 3.328e+04 2.365 0.0185 *
Costs 2.628e+01 1.933e+00 13.600 < 2e-16 ***
Trades -4.662e+02 2.087e+02 -2.234 0.0260 *
Dividends 3.513e+00 4.921e-01 7.138 4.25e-12 ***
M1 4.531e+04 7.152e+04 0.634 0.5267
M2 3.026e+04 7.117e+04 0.425 0.6710
M3 1.364e+04 7.138e+04 0.191 0.8486
M4 -1.349e+05 7.144e+04 -1.889 0.0596 .
M5 -3.977e+04 7.143e+04 -0.557 0.5780
M6 -2.718e+04 7.137e+04 -0.381 0.7035
M7 -1.341e+04 7.119e+04 -0.188 0.8507
M8 -3.472e+04 7.142e+04 -0.486 0.6271
M9 1.513e+04 7.115e+04 0.213 0.8318
M10 -1.311e+04 7.157e+04 -0.183 0.8547
M11 -6.342e+04 7.166e+04 -0.885 0.3767
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 299700 on 415 degrees of freedom
Multiple R-squared: 0.567, Adjusted R-squared: 0.5514
F-statistic: 36.24 on 15 and 415 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 1.0000000 2.941805e-69 1.470903e-69
[2,] 1.0000000 2.801495e-106 1.400747e-106
[3,] 1.0000000 1.401116e-110 7.005578e-111
[4,] 1.0000000 2.975697e-120 1.487849e-120
[5,] 1.0000000 2.954009e-123 1.477004e-123
[6,] 1.0000000 1.568911e-135 7.844553e-136
[7,] 1.0000000 8.413709e-137 4.206855e-137
[8,] 1.0000000 1.099167e-137 5.495837e-138
[9,] 1.0000000 3.869193e-143 1.934596e-143
[10,] 1.0000000 2.344841e-147 1.172420e-147
[11,] 1.0000000 1.278893e-191 6.394467e-192
[12,] 1.0000000 1.577066e-202 7.885331e-203
[13,] 1.0000000 1.943662e-207 9.718312e-208
[14,] 1.0000000 2.375517e-210 1.187758e-210
[15,] 1.0000000 2.330796e-216 1.165398e-216
[16,] 1.0000000 2.256672e-218 1.128336e-218
[17,] 1.0000000 3.221510e-218 1.610755e-218
[18,] 1.0000000 4.898575e-218 2.449288e-218
[19,] 1.0000000 4.742914e-225 2.371457e-225
[20,] 1.0000000 1.948744e-236 9.743721e-237
[21,] 1.0000000 1.092052e-237 5.460258e-238
[22,] 1.0000000 2.275289e-238 1.137645e-238
[23,] 1.0000000 3.015279e-238 1.507639e-238
[24,] 1.0000000 1.051123e-238 5.255616e-239
[25,] 1.0000000 1.294814e-239 6.474068e-240
[26,] 1.0000000 6.679025e-239 3.339513e-239
[27,] 1.0000000 2.201441e-239 1.100721e-239
[28,] 1.0000000 4.107364e-239 2.053682e-239
[29,] 1.0000000 1.101589e-238 5.507944e-239
[30,] 1.0000000 2.793677e-239 1.396839e-239
[31,] 1.0000000 5.256679e-239 2.628339e-239
[32,] 1.0000000 4.475597e-239 2.237798e-239
[33,] 1.0000000 6.485537e-240 3.242769e-240
[34,] 1.0000000 3.921773e-239 1.960887e-239
[35,] 1.0000000 1.869029e-239 9.345144e-240
[36,] 1.0000000 1.153079e-238 5.765393e-239
[37,] 1.0000000 2.102602e-241 1.051301e-241
[38,] 1.0000000 4.328715e-242 2.164357e-242
[39,] 1.0000000 4.340002e-247 2.170001e-247
[40,] 1.0000000 2.925833e-247 1.462917e-247
[41,] 1.0000000 6.758724e-247 3.379362e-247
[42,] 1.0000000 8.012793e-247 4.006396e-247
[43,] 1.0000000 4.344185e-249 2.172093e-249
[44,] 1.0000000 2.003691e-248 1.001846e-248
[45,] 1.0000000 1.628414e-247 8.142070e-248
[46,] 1.0000000 9.002596e-247 4.501298e-247
[47,] 1.0000000 6.586270e-247 3.293135e-247
[48,] 1.0000000 5.501797e-246 2.750898e-246
[49,] 1.0000000 2.022656e-246 1.011328e-246
[50,] 1.0000000 2.494422e-245 1.247211e-245
[51,] 1.0000000 3.006567e-244 1.503284e-244
[52,] 1.0000000 2.157902e-243 1.078951e-243
[53,] 1.0000000 1.233786e-242 6.168932e-243
[54,] 1.0000000 1.934497e-244 9.672486e-245
[55,] 1.0000000 2.381179e-244 1.190590e-244
[56,] 1.0000000 2.654891e-243 1.327445e-243
[57,] 1.0000000 1.529057e-242 7.645283e-243
[58,] 1.0000000 7.378846e-242 3.689423e-242
[59,] 1.0000000 2.118425e-241 1.059213e-241
[60,] 1.0000000 5.504914e-241 2.752457e-241
[61,] 1.0000000 3.285034e-240 1.642517e-240
[62,] 1.0000000 1.828522e-239 9.142609e-240
[63,] 1.0000000 1.742933e-238 8.714663e-239
[64,] 1.0000000 4.386384e-238 2.193192e-238
[65,] 1.0000000 6.644184e-238 3.322092e-238
[66,] 1.0000000 3.532648e-237 1.766324e-237
[67,] 1.0000000 4.481949e-236 2.240974e-236
[68,] 1.0000000 1.343845e-236 6.719223e-237
[69,] 1.0000000 1.260994e-236 6.304970e-237
[70,] 1.0000000 2.049242e-236 1.024621e-236
[71,] 1.0000000 9.502894e-236 4.751447e-236
[72,] 1.0000000 1.004283e-234 5.021416e-235
[73,] 1.0000000 1.321405e-233 6.607025e-234
[74,] 1.0000000 1.487256e-232 7.436278e-233
[75,] 1.0000000 4.274861e-233 2.137430e-233
[76,] 1.0000000 9.845275e-233 4.922637e-233
[77,] 1.0000000 1.366558e-234 6.832788e-235
[78,] 1.0000000 2.551564e-234 1.275782e-234
[79,] 1.0000000 3.348561e-233 1.674280e-233
[80,] 1.0000000 1.629226e-232 8.146132e-233
[81,] 1.0000000 7.759503e-232 3.879751e-232
[82,] 1.0000000 7.580770e-231 3.790385e-231
[83,] 1.0000000 8.440279e-230 4.220139e-230
[84,] 1.0000000 4.101452e-230 2.050726e-230
[85,] 1.0000000 3.608508e-229 1.804254e-229
[86,] 1.0000000 6.447989e-229 3.223994e-229
[87,] 1.0000000 4.952683e-228 2.476341e-228
[88,] 1.0000000 5.543845e-228 2.771922e-228
[89,] 1.0000000 3.807794e-228 1.903897e-228
[90,] 1.0000000 4.145838e-228 2.072919e-228
[91,] 1.0000000 1.606167e-227 8.030835e-228
[92,] 1.0000000 3.877497e-227 1.938748e-227
[93,] 1.0000000 4.175775e-227 2.087887e-227
[94,] 1.0000000 3.528736e-227 1.764368e-227
[95,] 1.0000000 8.107497e-229 4.053748e-229
[96,] 1.0000000 3.521886e-228 1.760943e-228
[97,] 1.0000000 1.113767e-227 5.568834e-228
[98,] 1.0000000 3.047058e-227 1.523529e-227
[99,] 1.0000000 2.713747e-226 1.356874e-226
[100,] 1.0000000 9.420159e-226 4.710079e-226
[101,] 1.0000000 9.621586e-227 4.810793e-227
[102,] 1.0000000 9.092470e-226 4.546235e-226
[103,] 1.0000000 5.911602e-226 2.955801e-226
[104,] 1.0000000 6.783466e-225 3.391733e-225
[105,] 1.0000000 6.426902e-224 3.213451e-224
[106,] 1.0000000 3.500811e-223 1.750406e-223
[107,] 1.0000000 4.505707e-222 2.252854e-222
[108,] 1.0000000 1.212727e-222 6.063636e-223
[109,] 1.0000000 1.320725e-221 6.603624e-222
[110,] 1.0000000 1.450489e-220 7.252443e-221
[111,] 1.0000000 1.865639e-219 9.328197e-220
[112,] 1.0000000 1.933378e-218 9.666892e-219
[113,] 1.0000000 1.441552e-217 7.207758e-218
[114,] 1.0000000 1.102007e-217 5.510033e-218
[115,] 1.0000000 1.494501e-216 7.472504e-217
[116,] 1.0000000 1.991758e-215 9.958792e-216
[117,] 1.0000000 2.463684e-214 1.231842e-214
[118,] 1.0000000 7.972680e-214 3.986340e-214
[119,] 1.0000000 7.339706e-213 3.669853e-213
[120,] 1.0000000 7.351229e-212 3.675615e-212
[121,] 1.0000000 8.384224e-211 4.192112e-211
[122,] 1.0000000 9.745841e-210 4.872920e-210
[123,] 1.0000000 1.246152e-208 6.230759e-209
[124,] 1.0000000 1.359925e-207 6.799627e-208
[125,] 1.0000000 7.536415e-207 3.768207e-207
[126,] 1.0000000 8.592287e-206 4.296143e-206
[127,] 1.0000000 8.976044e-205 4.488022e-205
[128,] 1.0000000 9.182890e-204 4.591445e-204
[129,] 1.0000000 1.127261e-202 5.636306e-203
[130,] 1.0000000 5.467021e-202 2.733511e-202
[131,] 1.0000000 5.409526e-201 2.704763e-201
[132,] 1.0000000 5.700696e-200 2.850348e-200
[133,] 1.0000000 6.142554e-199 3.071277e-199
[134,] 1.0000000 5.338507e-198 2.669254e-198
[135,] 1.0000000 6.586177e-197 3.293088e-197
[136,] 1.0000000 6.702569e-196 3.351285e-196
[137,] 1.0000000 6.185776e-195 3.092888e-195
[138,] 1.0000000 6.927063e-194 3.463531e-194
[139,] 1.0000000 8.350915e-193 4.175457e-193
[140,] 1.0000000 9.971353e-192 4.985677e-192
[141,] 1.0000000 1.004645e-190 5.023224e-191
[142,] 1.0000000 5.770491e-190 2.885246e-190
[143,] 1.0000000 5.810473e-189 2.905237e-189
[144,] 1.0000000 6.123356e-188 3.061678e-188
[145,] 1.0000000 6.685978e-187 3.342989e-187
[146,] 1.0000000 6.976678e-186 3.488339e-186
[147,] 1.0000000 8.167824e-185 4.083912e-185
[148,] 1.0000000 8.585861e-184 4.292930e-184
[149,] 1.0000000 7.971005e-183 3.985502e-183
[150,] 1.0000000 7.126498e-184 3.563249e-184
[151,] 1.0000000 8.242743e-183 4.121372e-183
[152,] 1.0000000 9.466727e-182 4.733363e-182
[153,] 1.0000000 1.090390e-180 5.451948e-181
[154,] 1.0000000 7.231917e-180 3.615959e-180
[155,] 1.0000000 7.359436e-179 3.679718e-179
[156,] 1.0000000 7.753882e-178 3.876941e-178
[157,] 1.0000000 8.265849e-177 4.132925e-177
[158,] 1.0000000 9.096596e-176 4.548298e-176
[159,] 1.0000000 1.026884e-174 5.134419e-175
[160,] 1.0000000 1.104579e-173 5.522897e-174
[161,] 1.0000000 1.028828e-172 5.144139e-173
[162,] 1.0000000 1.105544e-171 5.527721e-172
[163,] 1.0000000 5.535408e-171 2.767704e-171
[164,] 1.0000000 5.543179e-170 2.771589e-170
[165,] 1.0000000 3.887074e-169 1.943537e-169
[166,] 1.0000000 2.199133e-168 1.099567e-168
[167,] 1.0000000 2.393133e-167 1.196567e-167
[168,] 1.0000000 2.506453e-166 1.253227e-166
[169,] 1.0000000 1.748126e-165 8.740628e-166
[170,] 1.0000000 1.847810e-164 9.239051e-165
[171,] 1.0000000 1.381292e-163 6.906462e-164
[172,] 1.0000000 1.354472e-162 6.772359e-163
[173,] 1.0000000 1.026893e-161 5.134467e-162
[174,] 1.0000000 1.062020e-160 5.310102e-161
[175,] 1.0000000 1.040761e-159 5.203804e-160
[176,] 1.0000000 1.077060e-158 5.385302e-159
[177,] 1.0000000 1.123061e-157 5.615303e-158
[178,] 1.0000000 1.020016e-156 5.100078e-157
[179,] 1.0000000 3.069498e-156 1.534749e-156
[180,] 1.0000000 3.201214e-155 1.600607e-155
[181,] 1.0000000 3.184626e-154 1.592313e-154
[182,] 1.0000000 7.871464e-154 3.935732e-154
[183,] 1.0000000 7.693796e-153 3.846898e-153
[184,] 1.0000000 7.319398e-152 3.659699e-152
[185,] 1.0000000 7.182237e-151 3.591118e-151
[186,] 1.0000000 6.130241e-150 3.065120e-150
[187,] 1.0000000 1.357700e-149 6.788502e-150
[188,] 1.0000000 1.307978e-148 6.539892e-149
[189,] 1.0000000 1.195505e-147 5.977524e-148
[190,] 1.0000000 1.063572e-146 5.317858e-147
[191,] 1.0000000 1.159517e-146 5.797584e-147
[192,] 1.0000000 1.171457e-145 5.857285e-146
[193,] 1.0000000 1.083696e-144 5.418481e-145
[194,] 1.0000000 1.059788e-143 5.298939e-144
[195,] 1.0000000 1.012633e-142 5.063164e-143
[196,] 1.0000000 8.286056e-142 4.143028e-142
[197,] 1.0000000 7.177807e-141 3.588904e-141
[198,] 1.0000000 6.869895e-140 3.434948e-140
[199,] 1.0000000 6.540345e-139 3.270172e-139
[200,] 1.0000000 5.612609e-138 2.806305e-138
[201,] 1.0000000 5.398515e-137 2.699258e-137
[202,] 1.0000000 1.834842e-136 9.174211e-137
[203,] 1.0000000 1.436863e-135 7.184317e-136
[204,] 1.0000000 1.311197e-134 6.555985e-135
[205,] 1.0000000 1.154118e-133 5.770588e-134
[206,] 1.0000000 1.080232e-132 5.401160e-133
[207,] 1.0000000 1.018488e-131 5.092438e-132
[208,] 1.0000000 9.175052e-131 4.587526e-131
[209,] 1.0000000 7.495193e-130 3.747597e-130
[210,] 1.0000000 6.852965e-129 3.426483e-129
[211,] 1.0000000 6.101010e-128 3.050505e-128
[212,] 1.0000000 5.397793e-127 2.698896e-127
[213,] 1.0000000 4.868744e-126 2.434372e-126
[214,] 1.0000000 3.366144e-125 1.683072e-125
[215,] 1.0000000 2.892967e-124 1.446483e-124
[216,] 1.0000000 2.436942e-123 1.218471e-123
[217,] 1.0000000 2.035962e-122 1.017981e-122
[218,] 1.0000000 1.809788e-121 9.048938e-122
[219,] 1.0000000 1.613978e-120 8.069888e-121
[220,] 1.0000000 1.360039e-119 6.800194e-120
[221,] 1.0000000 1.058236e-118 5.291179e-119
[222,] 1.0000000 9.182995e-118 4.591498e-118
[223,] 1.0000000 7.908310e-117 3.954155e-117
[224,] 1.0000000 6.284929e-116 3.142465e-116
[225,] 1.0000000 5.444925e-115 2.722463e-115
[226,] 1.0000000 3.584225e-114 1.792113e-114
[227,] 1.0000000 3.070324e-113 1.535162e-113
[228,] 1.0000000 2.636516e-112 1.318258e-112
[229,] 1.0000000 2.153023e-111 1.076511e-111
[230,] 1.0000000 1.790512e-110 8.952560e-111
[231,] 1.0000000 1.465542e-109 7.327712e-110
[232,] 1.0000000 1.124025e-108 5.620125e-109
[233,] 1.0000000 6.997735e-108 3.498868e-108
[234,] 1.0000000 5.817358e-107 2.908679e-107
[235,] 1.0000000 4.723700e-106 2.361850e-106
[236,] 1.0000000 3.790274e-105 1.895137e-105
[237,] 1.0000000 1.602432e-104 8.012161e-105
[238,] 1.0000000 1.059120e-103 5.295602e-104
[239,] 1.0000000 6.070422e-103 3.035211e-103
[240,] 1.0000000 2.846074e-102 1.423037e-102
[241,] 1.0000000 2.212079e-101 1.106040e-101
[242,] 1.0000000 1.780672e-100 8.903361e-101
[243,] 1.0000000 1.396310e-99 6.981548e-100
[244,] 1.0000000 1.015339e-98 5.076694e-99
[245,] 1.0000000 7.921415e-98 3.960707e-98
[246,] 1.0000000 5.658439e-97 2.829219e-97
[247,] 1.0000000 4.301282e-96 2.150641e-96
[248,] 1.0000000 3.194853e-95 1.597427e-95
[249,] 1.0000000 2.456480e-94 1.228240e-94
[250,] 1.0000000 1.541662e-93 7.708311e-94
[251,] 1.0000000 1.192863e-92 5.964313e-93
[252,] 1.0000000 8.703443e-92 4.351722e-92
[253,] 1.0000000 6.393698e-91 3.196849e-91
[254,] 1.0000000 4.846484e-90 2.423242e-90
[255,] 1.0000000 3.583137e-89 1.791569e-89
[256,] 1.0000000 2.476972e-88 1.238486e-88
[257,] 1.0000000 1.633560e-87 8.167800e-88
[258,] 1.0000000 1.165308e-86 5.826541e-87
[259,] 1.0000000 8.273244e-86 4.136622e-86
[260,] 1.0000000 5.761573e-85 2.880787e-85
[261,] 1.0000000 3.746063e-84 1.873032e-84
[262,] 1.0000000 2.462775e-83 1.231387e-83
[263,] 1.0000000 1.792387e-82 8.961933e-83
[264,] 1.0000000 1.271505e-81 6.357527e-82
[265,] 1.0000000 6.760946e-81 3.380473e-81
[266,] 1.0000000 4.846771e-80 2.423386e-80
[267,] 1.0000000 2.669705e-79 1.334853e-79
[268,] 1.0000000 1.765201e-78 8.826003e-79
[269,] 1.0000000 1.176541e-77 5.882705e-78
[270,] 1.0000000 8.210523e-77 4.105262e-77
[271,] 1.0000000 5.607142e-76 2.803571e-76
[272,] 1.0000000 3.082963e-75 1.541482e-75
[273,] 1.0000000 2.104632e-74 1.052316e-74
[274,] 1.0000000 1.186672e-73 5.933358e-74
[275,] 1.0000000 7.315189e-73 3.657594e-73
[276,] 1.0000000 4.916931e-72 2.458466e-72
[277,] 1.0000000 2.996756e-71 1.498378e-71
[278,] 1.0000000 1.546938e-70 7.734689e-71
[279,] 1.0000000 1.032840e-69 5.164201e-70
[280,] 1.0000000 6.108002e-69 3.054001e-69
[281,] 1.0000000 3.577420e-68 1.788710e-68
[282,] 1.0000000 2.342832e-67 1.171416e-67
[283,] 1.0000000 1.520375e-66 7.601876e-67
[284,] 1.0000000 9.336960e-66 4.668480e-66
[285,] 1.0000000 6.100591e-65 3.050296e-65
[286,] 1.0000000 3.755335e-64 1.877668e-64
[287,] 1.0000000 1.960133e-63 9.800663e-64
[288,] 1.0000000 1.236220e-62 6.181099e-63
[289,] 1.0000000 7.419573e-62 3.709786e-62
[290,] 1.0000000 4.670125e-61 2.335063e-61
[291,] 1.0000000 2.697774e-60 1.348887e-60
[292,] 1.0000000 1.464117e-59 7.320585e-60
[293,] 1.0000000 8.575283e-59 4.287641e-59
[294,] 1.0000000 4.007466e-58 2.003733e-58
[295,] 1.0000000 2.430232e-57 1.215116e-57
[296,] 1.0000000 1.427726e-56 7.138629e-57
[297,] 1.0000000 8.727023e-56 4.363512e-56
[298,] 1.0000000 4.723454e-55 2.361727e-55
[299,] 1.0000000 2.804461e-54 1.402231e-54
[300,] 1.0000000 1.677392e-53 8.386958e-54
[301,] 1.0000000 8.851561e-53 4.425781e-53
[302,] 1.0000000 5.226507e-52 2.613254e-52
[303,] 1.0000000 3.065840e-51 1.532920e-51
[304,] 1.0000000 1.625513e-50 8.127563e-51
[305,] 1.0000000 8.152298e-50 4.076149e-50
[306,] 1.0000000 4.719095e-49 2.359548e-49
[307,] 1.0000000 2.643374e-48 1.321687e-48
[308,] 1.0000000 1.384213e-47 6.921064e-48
[309,] 1.0000000 7.756228e-47 3.878114e-47
[310,] 1.0000000 4.090284e-46 2.045142e-46
[311,] 1.0000000 2.262971e-45 1.131486e-45
[312,] 1.0000000 1.239189e-44 6.195947e-45
[313,] 1.0000000 5.981734e-44 2.990867e-44
[314,] 1.0000000 3.265608e-43 1.632804e-43
[315,] 1.0000000 1.770469e-42 8.852347e-43
[316,] 1.0000000 8.724730e-42 4.362365e-42
[317,] 1.0000000 3.971450e-41 1.985725e-41
[318,] 1.0000000 2.110905e-40 1.055452e-40
[319,] 1.0000000 1.091248e-39 5.456238e-40
[320,] 1.0000000 5.067294e-39 2.533647e-39
[321,] 1.0000000 2.645497e-38 1.322749e-38
[322,] 1.0000000 1.182055e-37 5.910277e-38
[323,] 1.0000000 5.953598e-37 2.976799e-37
[324,] 1.0000000 2.994441e-36 1.497220e-36
[325,] 1.0000000 1.308535e-35 6.542674e-36
[326,] 1.0000000 6.558209e-35 3.279105e-35
[327,] 1.0000000 3.265558e-34 1.632779e-34
[328,] 1.0000000 1.452124e-33 7.260619e-34
[329,] 1.0000000 6.231999e-33 3.115999e-33
[330,] 1.0000000 3.033429e-32 1.516714e-32
[331,] 1.0000000 1.419982e-31 7.099910e-32
[332,] 1.0000000 5.855368e-31 2.927684e-31
[333,] 1.0000000 2.796118e-30 1.398059e-30
[334,] 1.0000000 1.162091e-29 5.810454e-30
[335,] 1.0000000 5.128022e-29 2.564011e-29
[336,] 1.0000000 2.374154e-28 1.187077e-28
[337,] 1.0000000 9.256266e-28 4.628133e-28
[338,] 1.0000000 4.209962e-27 2.104981e-27
[339,] 1.0000000 1.362964e-26 6.814819e-27
[340,] 1.0000000 5.524688e-26 2.762344e-26
[341,] 1.0000000 1.976074e-25 9.880368e-26
[342,] 1.0000000 8.764058e-25 4.382029e-25
[343,] 1.0000000 3.708098e-24 1.854049e-24
[344,] 1.0000000 1.287340e-23 6.436700e-24
[345,] 1.0000000 5.602286e-23 2.801143e-23
[346,] 1.0000000 2.283529e-22 1.141765e-22
[347,] 1.0000000 9.657493e-22 4.828746e-22
[348,] 1.0000000 3.961051e-21 1.980525e-21
[349,] 1.0000000 1.350552e-20 6.752761e-21
[350,] 1.0000000 5.507475e-20 2.753738e-20
[351,] 1.0000000 2.214132e-19 1.107066e-19
[352,] 1.0000000 8.000809e-19 4.000405e-19
[353,] 1.0000000 2.408746e-18 1.204373e-18
[354,] 1.0000000 9.590255e-18 4.795127e-18
[355,] 1.0000000 3.662111e-17 1.831056e-17
[356,] 1.0000000 1.250622e-16 6.253109e-17
[357,] 1.0000000 4.819349e-16 2.409674e-16
[358,] 1.0000000 1.519787e-15 7.598933e-16
[359,] 1.0000000 5.647197e-15 2.823599e-15
[360,] 1.0000000 1.985017e-14 9.925087e-15
[361,] 1.0000000 5.721411e-14 2.860706e-14
[362,] 1.0000000 2.018040e-13 1.009020e-13
[363,] 1.0000000 7.126554e-13 3.563277e-13
[364,] 1.0000000 1.556350e-12 7.781749e-13
[365,] 1.0000000 2.807860e-12 1.403930e-12
[366,] 1.0000000 9.898953e-12 4.949476e-12
[367,] 1.0000000 3.416875e-11 1.708438e-11
[368,] 1.0000000 1.157801e-10 5.789003e-11
[369,] 1.0000000 3.915537e-10 1.957769e-10
[370,] 1.0000000 8.296856e-10 4.148428e-10
[371,] 1.0000000 2.705631e-09 1.352816e-09
[372,] 1.0000000 8.215896e-09 4.107948e-09
[373,] 1.0000000 1.847582e-08 9.237910e-09
[374,] 1.0000000 5.532005e-08 2.766003e-08
[375,] 0.9999999 1.614701e-07 8.073506e-08
[376,] 0.9999998 4.908308e-07 2.454154e-07
[377,] 0.9999993 1.421330e-06 7.106648e-07
[378,] 0.9999979 4.172868e-06 2.086434e-06
[379,] 0.9999946 1.071998e-05 5.359990e-06
[380,] 0.9999852 2.953232e-05 1.476616e-05
[381,] 0.9999598 8.031072e-05 4.015536e-05
[382,] 0.9999400 1.199593e-04 5.997965e-05
[383,] 0.9998402 3.195598e-04 1.597799e-04
[384,] 0.9996486 7.027633e-04 3.513816e-04
[385,] 0.9994059 1.188147e-03 5.940734e-04
[386,] 0.9985228 2.954388e-03 1.477194e-03
[387,] 0.9967882 6.423508e-03 3.211754e-03
[388,] 0.9932633 1.347344e-02 6.736722e-03
[389,] 0.9864646 2.707088e-02 1.353544e-02
[390,] 0.9710111 5.797781e-02 2.898890e-02
[391,] 0.9422911 1.154179e-01 5.770893e-02
[392,] 0.8841626 2.316748e-01 1.158374e-01
[393,] 0.8113009 3.773983e-01 1.886991e-01
[394,] 0.8270985 3.458031e-01 1.729015e-01
> postscript(file="/var/www/html/freestat/rcomp/tmp/18y2n1291320718.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/28y2n1291320718.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/38y2n1291320718.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/4jp281291320718.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/5jp281291320718.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
1569858.1477 3216804.5650 1318169.8718 -3011398.4292 -368870.5277
6 7 8 9 10
56132.3433 2088.9969 616587.3876 -245850.0480 760920.0958
11 12 13 14 15
-276567.2123 462075.4986 905234.7776 -1819604.9030 16860.5273
16 17 18 19 20
558579.8501 -391037.7520 -63811.2844 -152457.7990 -309261.7393
21 22 23 24 25
135480.0583 -253842.5020 -190611.9735 -504512.5871 -277146.2879
26 27 28 29 30
-171676.1145 -400257.9907 -321172.5938 429960.1254 36650.9527
31 32 33 34 35
169625.5040 -151503.2213 93830.9434 5325.1372 -294870.5359
36 37 38 39 40
-469237.4376 -5701.6948 18715.6838 -163183.1173 -30951.2419
41 42 43 44 45
43415.1163 -342227.2146 99832.9212 -168130.2910 -43157.5823
46 47 48 49 50
-243920.4219 -211537.0486 -197604.9207 -166326.6387 -203553.7208
51 52 53 54 55
-13020.2667 -16553.3606 -391448.4984 -47772.7110 -407078.5585
56 57 58 59 60
152861.1425 151174.2325 -21731.3504 25629.2641 -249681.1403
61 62 63 64 65
-158595.9187 -103861.8336 -147288.9587 -121221.0244 88306.9426
66 67 68 69 70
90107.6762 -95650.0029 -89941.3646 -91213.1364 -147913.3273
71 72 73 74 75
10949.3791 57506.3400 -188093.9212 -62163.5410 -152804.1801
76 77 78 79 80
-30854.7275 -70049.8374 -37018.9676 83844.8336 -47622.7565
81 82 83 84 85
20131.2434 -109659.4849 -73596.1474 73092.8848 -68100.5949
86 87 88 89 90
-165579.6510 -181.4139 78258.6304 -73327.3995 -10263.5817
91 92 93 94 95
-30223.2831 -92607.5185 104610.2651 -350911.9747 189394.8250
96 97 98 99 100
-46129.8832 -50431.2966 -7656.7252 -163892.6059 24718.5248
101 102 103 104 105
-78894.3948 -91097.9397 -114331.1095 -36644.2181 27870.6053
106 107 108 109 110
-47860.2078 -199900.1793 -25894.5879 -265545.6115 -206906.6017
111 112 113 114 115
85871.1995 -123771.6851 -245146.4930 -171532.3565 -18703.3689
116 117 118 119 120
97843.3432 -20602.5726 37320.8522 -198173.0011 61162.1655
121 122 123 124 125
-321611.0043 -2794.4570 -28581.9495 61586.0824 -105202.5728
126 127 128 129 130
-246889.1051 53482.9591 75085.1371 33180.3946 53183.0694
131 132 133 134 135
106685.9316 97823.1152 -5239.7914 7784.7346 26435.2879
136 137 138 139 140
177598.0709 82157.2142 67249.1672 53482.9591 74794.3078
141 142 143 144 145
24033.0410 53183.0694 98446.4628 40338.0453 1024.4114
146 147 148 149 150
21860.2380 26435.2879 176950.4688 79836.7782 66677.6384
151 152 153 154 155
44210.8872 64236.9466 24946.6516 30999.5221 103581.6414
156 157 158 159 160
40071.6901 -5239.7914 6640.6904 33635.8992 174997.5049
161 162 163 164 165
79836.7782 67249.1672 53482.9591 80570.6160 25606.0709
166 167 168 169 170
53183.0694 103494.3525 144000.5954 -5239.7914 9812.7590
171 172 173 174 175
26435.2879 171989.4257 79836.7782 64997.4719 53482.9591
176 177 178 179 180
74830.4635 24946.6516 -26761.3976 103494.3525 40071.6901
181 182 183 184 185
8870.4918 -92379.1625 52232.5495 34221.3870 391.3591
186 187 188 189 190
-18969.3233 -20035.0601 71459.9865 -14812.3441 53183.0694
191 192 193 194 195
108202.7436 40071.6901 -3080.7867 9812.7590 26159.6990
196 197 198 199 200
96294.5998 -135572.3399 -11453.7379 -25219.9460 -82122.4322
201 202 203 204 205
-53756.2535 -25519.8357 24791.4474 -10296.1637 -170578.7136
206 207 208 209 210
-70637.5945 -67514.7389 96294.5998 60339.1038 -11453.7379
211 212 213 214 215
53482.9591 74794.3078 40850.3406 30976.8117 84784.1493
216 217 218 219 220
40071.6901 -1697.9686 14093.9080 26435.2879 180079.5502
221 222 223 224 225
47071.4215 67249.1672 54389.8956 74794.3078 25256.9161
226 227 228 229 230
57993.3973 103494.3525 40071.6901 -31464.8866 9812.7590
231 232 233 234 235
5291.3828 174997.5049 79836.7782 68678.2748 53482.9591
236 237 238 239 240
75053.1947 24946.6516 56168.9458 103494.3525 40071.6901
241 242 243 244 245
-5239.7914 12699.7489 26435.2879 180228.4215 -7166.6580
246 247 248 249 250
-11453.7379 -25219.9460 66680.0605 -53756.2535 -25519.8357
251 252 253 254 255
110695.2132 -38631.2150 -28672.4846 -9733.9543 -24628.0072
256 257 258 259 260
174997.5049 103324.4741 107375.5412 -25219.9460 74794.3078
261 262 263 264 265
-53756.2535 -25519.8357 83708.2143 48712.1463 -83942.6965
266 267 268 269 270
9812.7590 -52267.6172 178111.3539 1133.8731 45606.4224
271 272 273 274 275
-25219.9460 -3908.5973 -53756.2535 51734.8475 103494.3525
276 277 278 279 280
49275.3029 1841.8182 9812.7590 -63545.2223 153959.8020
281 282 283 284 285
1133.8731 67249.1672 92330.4772 74794.3078 -11005.1296
286 287 288 289 290
53183.0694 24791.4474 39090.5666 -83942.6965 40682.4378
291 292 293 294 295
-52267.6172 184802.0237 52906.5441 67249.1672 57166.2433
296 297 298 299 300
30865.1386 24946.6516 -25519.8357 103494.3525 40071.6901
301 302 303 304 305
-84653.1602 -4192.3673 26435.2879 96294.5998 118963.7940
306 307 308 309 310
71370.1413 -25219.9460 -11902.9183 46885.7458 -25519.8357
311 312 313 314 315
24791.4474 19998.4788 -83942.6965 -71893.3472 28388.4039
316 317 318 319 320
176121.2446 81440.5043 53205.5005 56254.6782 74794.3078
321 322 323 324 325
24946.6516 53183.0694 103494.3525 -38631.2150 -5239.7914
326 327 328 329 330
7249.0362 -52267.6172 96294.5998 79836.7782 67249.1672
331 332 333 334 335
53482.9591 74794.3078 24946.6516 53227.2049 104382.9762
336 337 338 339 340
40071.6901 -6380.0132 9812.7590 25787.9362 179062.3190
341 342 343 344 345
60789.2793 67249.1672 53482.9591 74794.3078 24946.6516
346 347 348 349 350
53183.0694 95507.1026 40071.6901 -5239.7914 9812.7590
351 352 353 354 355
27233.4321 176349.8594 81357.5962 57358.7320 53482.9591
356 357 358 359 360
74794.3078 -8068.1788 53183.0694 103494.3525 40071.6901
361 362 363 364 365
-3409.1070 9812.7590 26435.2879 87151.6133 78437.3973
366 367 368 369 370
67249.1672 53482.9591 74794.3078 24946.6516 53183.0694
371 372 373 374 375
103494.3525 43205.1677 -45287.8150 -17023.1039 26435.2879
376 377 378 379 380
174997.5049 92267.5274 52111.7090 53482.9591 74794.3078
381 382 383 384 385
24946.6516 46373.5097 157024.3228 40071.6901 -3100.4391
386 387 388 389 390
-146528.7581 31132.1234 159162.8844 78031.4273 47984.8525
391 392 393 394 395
50297.8980 54309.1973 -12773.6816 69737.5795 -186404.4461
396 397 398 399 400
12487.0534 -169009.4846 -57017.7005 -126.2104 90226.5753
401 402 403 404 405
41927.2186 57105.4751 32843.7778 -86027.3038 -8671.2602
406 407 408 409 410
1044.3173 87763.6077 19571.2178 -27565.2565 -65577.8594
411 412 413 414 415
-18282.3546 -107485.7438 -28303.5090 -151887.3283 -57958.1541
416 417 418 419 420
-19874.0729 107460.5380 -79034.8951 52454.3350 -28508.0185
421 422 423 424 425
-127109.7247 -146251.7195 -462135.4596 -350917.6996 -47518.7003
426 427 428 429 430
-197525.0426 -308661.5966 -1183573.5658 -389711.3115 -271235.1059
431
-793373.1387
> postscript(file="/var/www/html/freestat/rcomp/tmp/6cg1b1291320718.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 1569858.1477 NA
1 3216804.5650 1569858.1477
2 1318169.8718 3216804.5650
3 -3011398.4292 1318169.8718
4 -368870.5277 -3011398.4292
5 56132.3433 -368870.5277
6 2088.9969 56132.3433
7 616587.3876 2088.9969
8 -245850.0480 616587.3876
9 760920.0958 -245850.0480
10 -276567.2123 760920.0958
11 462075.4986 -276567.2123
12 905234.7776 462075.4986
13 -1819604.9030 905234.7776
14 16860.5273 -1819604.9030
15 558579.8501 16860.5273
16 -391037.7520 558579.8501
17 -63811.2844 -391037.7520
18 -152457.7990 -63811.2844
19 -309261.7393 -152457.7990
20 135480.0583 -309261.7393
21 -253842.5020 135480.0583
22 -190611.9735 -253842.5020
23 -504512.5871 -190611.9735
24 -277146.2879 -504512.5871
25 -171676.1145 -277146.2879
26 -400257.9907 -171676.1145
27 -321172.5938 -400257.9907
28 429960.1254 -321172.5938
29 36650.9527 429960.1254
30 169625.5040 36650.9527
31 -151503.2213 169625.5040
32 93830.9434 -151503.2213
33 5325.1372 93830.9434
34 -294870.5359 5325.1372
35 -469237.4376 -294870.5359
36 -5701.6948 -469237.4376
37 18715.6838 -5701.6948
38 -163183.1173 18715.6838
39 -30951.2419 -163183.1173
40 43415.1163 -30951.2419
41 -342227.2146 43415.1163
42 99832.9212 -342227.2146
43 -168130.2910 99832.9212
44 -43157.5823 -168130.2910
45 -243920.4219 -43157.5823
46 -211537.0486 -243920.4219
47 -197604.9207 -211537.0486
48 -166326.6387 -197604.9207
49 -203553.7208 -166326.6387
50 -13020.2667 -203553.7208
51 -16553.3606 -13020.2667
52 -391448.4984 -16553.3606
53 -47772.7110 -391448.4984
54 -407078.5585 -47772.7110
55 152861.1425 -407078.5585
56 151174.2325 152861.1425
57 -21731.3504 151174.2325
58 25629.2641 -21731.3504
59 -249681.1403 25629.2641
60 -158595.9187 -249681.1403
61 -103861.8336 -158595.9187
62 -147288.9587 -103861.8336
63 -121221.0244 -147288.9587
64 88306.9426 -121221.0244
65 90107.6762 88306.9426
66 -95650.0029 90107.6762
67 -89941.3646 -95650.0029
68 -91213.1364 -89941.3646
69 -147913.3273 -91213.1364
70 10949.3791 -147913.3273
71 57506.3400 10949.3791
72 -188093.9212 57506.3400
73 -62163.5410 -188093.9212
74 -152804.1801 -62163.5410
75 -30854.7275 -152804.1801
76 -70049.8374 -30854.7275
77 -37018.9676 -70049.8374
78 83844.8336 -37018.9676
79 -47622.7565 83844.8336
80 20131.2434 -47622.7565
81 -109659.4849 20131.2434
82 -73596.1474 -109659.4849
83 73092.8848 -73596.1474
84 -68100.5949 73092.8848
85 -165579.6510 -68100.5949
86 -181.4139 -165579.6510
87 78258.6304 -181.4139
88 -73327.3995 78258.6304
89 -10263.5817 -73327.3995
90 -30223.2831 -10263.5817
91 -92607.5185 -30223.2831
92 104610.2651 -92607.5185
93 -350911.9747 104610.2651
94 189394.8250 -350911.9747
95 -46129.8832 189394.8250
96 -50431.2966 -46129.8832
97 -7656.7252 -50431.2966
98 -163892.6059 -7656.7252
99 24718.5248 -163892.6059
100 -78894.3948 24718.5248
101 -91097.9397 -78894.3948
102 -114331.1095 -91097.9397
103 -36644.2181 -114331.1095
104 27870.6053 -36644.2181
105 -47860.2078 27870.6053
106 -199900.1793 -47860.2078
107 -25894.5879 -199900.1793
108 -265545.6115 -25894.5879
109 -206906.6017 -265545.6115
110 85871.1995 -206906.6017
111 -123771.6851 85871.1995
112 -245146.4930 -123771.6851
113 -171532.3565 -245146.4930
114 -18703.3689 -171532.3565
115 97843.3432 -18703.3689
116 -20602.5726 97843.3432
117 37320.8522 -20602.5726
118 -198173.0011 37320.8522
119 61162.1655 -198173.0011
120 -321611.0043 61162.1655
121 -2794.4570 -321611.0043
122 -28581.9495 -2794.4570
123 61586.0824 -28581.9495
124 -105202.5728 61586.0824
125 -246889.1051 -105202.5728
126 53482.9591 -246889.1051
127 75085.1371 53482.9591
128 33180.3946 75085.1371
129 53183.0694 33180.3946
130 106685.9316 53183.0694
131 97823.1152 106685.9316
132 -5239.7914 97823.1152
133 7784.7346 -5239.7914
134 26435.2879 7784.7346
135 177598.0709 26435.2879
136 82157.2142 177598.0709
137 67249.1672 82157.2142
138 53482.9591 67249.1672
139 74794.3078 53482.9591
140 24033.0410 74794.3078
141 53183.0694 24033.0410
142 98446.4628 53183.0694
143 40338.0453 98446.4628
144 1024.4114 40338.0453
145 21860.2380 1024.4114
146 26435.2879 21860.2380
147 176950.4688 26435.2879
148 79836.7782 176950.4688
149 66677.6384 79836.7782
150 44210.8872 66677.6384
151 64236.9466 44210.8872
152 24946.6516 64236.9466
153 30999.5221 24946.6516
154 103581.6414 30999.5221
155 40071.6901 103581.6414
156 -5239.7914 40071.6901
157 6640.6904 -5239.7914
158 33635.8992 6640.6904
159 174997.5049 33635.8992
160 79836.7782 174997.5049
161 67249.1672 79836.7782
162 53482.9591 67249.1672
163 80570.6160 53482.9591
164 25606.0709 80570.6160
165 53183.0694 25606.0709
166 103494.3525 53183.0694
167 144000.5954 103494.3525
168 -5239.7914 144000.5954
169 9812.7590 -5239.7914
170 26435.2879 9812.7590
171 171989.4257 26435.2879
172 79836.7782 171989.4257
173 64997.4719 79836.7782
174 53482.9591 64997.4719
175 74830.4635 53482.9591
176 24946.6516 74830.4635
177 -26761.3976 24946.6516
178 103494.3525 -26761.3976
179 40071.6901 103494.3525
180 8870.4918 40071.6901
181 -92379.1625 8870.4918
182 52232.5495 -92379.1625
183 34221.3870 52232.5495
184 391.3591 34221.3870
185 -18969.3233 391.3591
186 -20035.0601 -18969.3233
187 71459.9865 -20035.0601
188 -14812.3441 71459.9865
189 53183.0694 -14812.3441
190 108202.7436 53183.0694
191 40071.6901 108202.7436
192 -3080.7867 40071.6901
193 9812.7590 -3080.7867
194 26159.6990 9812.7590
195 96294.5998 26159.6990
196 -135572.3399 96294.5998
197 -11453.7379 -135572.3399
198 -25219.9460 -11453.7379
199 -82122.4322 -25219.9460
200 -53756.2535 -82122.4322
201 -25519.8357 -53756.2535
202 24791.4474 -25519.8357
203 -10296.1637 24791.4474
204 -170578.7136 -10296.1637
205 -70637.5945 -170578.7136
206 -67514.7389 -70637.5945
207 96294.5998 -67514.7389
208 60339.1038 96294.5998
209 -11453.7379 60339.1038
210 53482.9591 -11453.7379
211 74794.3078 53482.9591
212 40850.3406 74794.3078
213 30976.8117 40850.3406
214 84784.1493 30976.8117
215 40071.6901 84784.1493
216 -1697.9686 40071.6901
217 14093.9080 -1697.9686
218 26435.2879 14093.9080
219 180079.5502 26435.2879
220 47071.4215 180079.5502
221 67249.1672 47071.4215
222 54389.8956 67249.1672
223 74794.3078 54389.8956
224 25256.9161 74794.3078
225 57993.3973 25256.9161
226 103494.3525 57993.3973
227 40071.6901 103494.3525
228 -31464.8866 40071.6901
229 9812.7590 -31464.8866
230 5291.3828 9812.7590
231 174997.5049 5291.3828
232 79836.7782 174997.5049
233 68678.2748 79836.7782
234 53482.9591 68678.2748
235 75053.1947 53482.9591
236 24946.6516 75053.1947
237 56168.9458 24946.6516
238 103494.3525 56168.9458
239 40071.6901 103494.3525
240 -5239.7914 40071.6901
241 12699.7489 -5239.7914
242 26435.2879 12699.7489
243 180228.4215 26435.2879
244 -7166.6580 180228.4215
245 -11453.7379 -7166.6580
246 -25219.9460 -11453.7379
247 66680.0605 -25219.9460
248 -53756.2535 66680.0605
249 -25519.8357 -53756.2535
250 110695.2132 -25519.8357
251 -38631.2150 110695.2132
252 -28672.4846 -38631.2150
253 -9733.9543 -28672.4846
254 -24628.0072 -9733.9543
255 174997.5049 -24628.0072
256 103324.4741 174997.5049
257 107375.5412 103324.4741
258 -25219.9460 107375.5412
259 74794.3078 -25219.9460
260 -53756.2535 74794.3078
261 -25519.8357 -53756.2535
262 83708.2143 -25519.8357
263 48712.1463 83708.2143
264 -83942.6965 48712.1463
265 9812.7590 -83942.6965
266 -52267.6172 9812.7590
267 178111.3539 -52267.6172
268 1133.8731 178111.3539
269 45606.4224 1133.8731
270 -25219.9460 45606.4224
271 -3908.5973 -25219.9460
272 -53756.2535 -3908.5973
273 51734.8475 -53756.2535
274 103494.3525 51734.8475
275 49275.3029 103494.3525
276 1841.8182 49275.3029
277 9812.7590 1841.8182
278 -63545.2223 9812.7590
279 153959.8020 -63545.2223
280 1133.8731 153959.8020
281 67249.1672 1133.8731
282 92330.4772 67249.1672
283 74794.3078 92330.4772
284 -11005.1296 74794.3078
285 53183.0694 -11005.1296
286 24791.4474 53183.0694
287 39090.5666 24791.4474
288 -83942.6965 39090.5666
289 40682.4378 -83942.6965
290 -52267.6172 40682.4378
291 184802.0237 -52267.6172
292 52906.5441 184802.0237
293 67249.1672 52906.5441
294 57166.2433 67249.1672
295 30865.1386 57166.2433
296 24946.6516 30865.1386
297 -25519.8357 24946.6516
298 103494.3525 -25519.8357
299 40071.6901 103494.3525
300 -84653.1602 40071.6901
301 -4192.3673 -84653.1602
302 26435.2879 -4192.3673
303 96294.5998 26435.2879
304 118963.7940 96294.5998
305 71370.1413 118963.7940
306 -25219.9460 71370.1413
307 -11902.9183 -25219.9460
308 46885.7458 -11902.9183
309 -25519.8357 46885.7458
310 24791.4474 -25519.8357
311 19998.4788 24791.4474
312 -83942.6965 19998.4788
313 -71893.3472 -83942.6965
314 28388.4039 -71893.3472
315 176121.2446 28388.4039
316 81440.5043 176121.2446
317 53205.5005 81440.5043
318 56254.6782 53205.5005
319 74794.3078 56254.6782
320 24946.6516 74794.3078
321 53183.0694 24946.6516
322 103494.3525 53183.0694
323 -38631.2150 103494.3525
324 -5239.7914 -38631.2150
325 7249.0362 -5239.7914
326 -52267.6172 7249.0362
327 96294.5998 -52267.6172
328 79836.7782 96294.5998
329 67249.1672 79836.7782
330 53482.9591 67249.1672
331 74794.3078 53482.9591
332 24946.6516 74794.3078
333 53227.2049 24946.6516
334 104382.9762 53227.2049
335 40071.6901 104382.9762
336 -6380.0132 40071.6901
337 9812.7590 -6380.0132
338 25787.9362 9812.7590
339 179062.3190 25787.9362
340 60789.2793 179062.3190
341 67249.1672 60789.2793
342 53482.9591 67249.1672
343 74794.3078 53482.9591
344 24946.6516 74794.3078
345 53183.0694 24946.6516
346 95507.1026 53183.0694
347 40071.6901 95507.1026
348 -5239.7914 40071.6901
349 9812.7590 -5239.7914
350 27233.4321 9812.7590
351 176349.8594 27233.4321
352 81357.5962 176349.8594
353 57358.7320 81357.5962
354 53482.9591 57358.7320
355 74794.3078 53482.9591
356 -8068.1788 74794.3078
357 53183.0694 -8068.1788
358 103494.3525 53183.0694
359 40071.6901 103494.3525
360 -3409.1070 40071.6901
361 9812.7590 -3409.1070
362 26435.2879 9812.7590
363 87151.6133 26435.2879
364 78437.3973 87151.6133
365 67249.1672 78437.3973
366 53482.9591 67249.1672
367 74794.3078 53482.9591
368 24946.6516 74794.3078
369 53183.0694 24946.6516
370 103494.3525 53183.0694
371 43205.1677 103494.3525
372 -45287.8150 43205.1677
373 -17023.1039 -45287.8150
374 26435.2879 -17023.1039
375 174997.5049 26435.2879
376 92267.5274 174997.5049
377 52111.7090 92267.5274
378 53482.9591 52111.7090
379 74794.3078 53482.9591
380 24946.6516 74794.3078
381 46373.5097 24946.6516
382 157024.3228 46373.5097
383 40071.6901 157024.3228
384 -3100.4391 40071.6901
385 -146528.7581 -3100.4391
386 31132.1234 -146528.7581
387 159162.8844 31132.1234
388 78031.4273 159162.8844
389 47984.8525 78031.4273
390 50297.8980 47984.8525
391 54309.1973 50297.8980
392 -12773.6816 54309.1973
393 69737.5795 -12773.6816
394 -186404.4461 69737.5795
395 12487.0534 -186404.4461
396 -169009.4846 12487.0534
397 -57017.7005 -169009.4846
398 -126.2104 -57017.7005
399 90226.5753 -126.2104
400 41927.2186 90226.5753
401 57105.4751 41927.2186
402 32843.7778 57105.4751
403 -86027.3038 32843.7778
404 -8671.2602 -86027.3038
405 1044.3173 -8671.2602
406 87763.6077 1044.3173
407 19571.2178 87763.6077
408 -27565.2565 19571.2178
409 -65577.8594 -27565.2565
410 -18282.3546 -65577.8594
411 -107485.7438 -18282.3546
412 -28303.5090 -107485.7438
413 -151887.3283 -28303.5090
414 -57958.1541 -151887.3283
415 -19874.0729 -57958.1541
416 107460.5380 -19874.0729
417 -79034.8951 107460.5380
418 52454.3350 -79034.8951
419 -28508.0185 52454.3350
420 -127109.7247 -28508.0185
421 -146251.7195 -127109.7247
422 -462135.4596 -146251.7195
423 -350917.6996 -462135.4596
424 -47518.7003 -350917.6996
425 -197525.0426 -47518.7003
426 -308661.5966 -197525.0426
427 -1183573.5658 -308661.5966
428 -389711.3115 -1183573.5658
429 -271235.1059 -389711.3115
430 -793373.1387 -271235.1059
431 NA -793373.1387
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 3216804.5650 1569858.1477
[2,] 1318169.8718 3216804.5650
[3,] -3011398.4292 1318169.8718
[4,] -368870.5277 -3011398.4292
[5,] 56132.3433 -368870.5277
[6,] 2088.9969 56132.3433
[7,] 616587.3876 2088.9969
[8,] -245850.0480 616587.3876
[9,] 760920.0958 -245850.0480
[10,] -276567.2123 760920.0958
[11,] 462075.4986 -276567.2123
[12,] 905234.7776 462075.4986
[13,] -1819604.9030 905234.7776
[14,] 16860.5273 -1819604.9030
[15,] 558579.8501 16860.5273
[16,] -391037.7520 558579.8501
[17,] -63811.2844 -391037.7520
[18,] -152457.7990 -63811.2844
[19,] -309261.7393 -152457.7990
[20,] 135480.0583 -309261.7393
[21,] -253842.5020 135480.0583
[22,] -190611.9735 -253842.5020
[23,] -504512.5871 -190611.9735
[24,] -277146.2879 -504512.5871
[25,] -171676.1145 -277146.2879
[26,] -400257.9907 -171676.1145
[27,] -321172.5938 -400257.9907
[28,] 429960.1254 -321172.5938
[29,] 36650.9527 429960.1254
[30,] 169625.5040 36650.9527
[31,] -151503.2213 169625.5040
[32,] 93830.9434 -151503.2213
[33,] 5325.1372 93830.9434
[34,] -294870.5359 5325.1372
[35,] -469237.4376 -294870.5359
[36,] -5701.6948 -469237.4376
[37,] 18715.6838 -5701.6948
[38,] -163183.1173 18715.6838
[39,] -30951.2419 -163183.1173
[40,] 43415.1163 -30951.2419
[41,] -342227.2146 43415.1163
[42,] 99832.9212 -342227.2146
[43,] -168130.2910 99832.9212
[44,] -43157.5823 -168130.2910
[45,] -243920.4219 -43157.5823
[46,] -211537.0486 -243920.4219
[47,] -197604.9207 -211537.0486
[48,] -166326.6387 -197604.9207
[49,] -203553.7208 -166326.6387
[50,] -13020.2667 -203553.7208
[51,] -16553.3606 -13020.2667
[52,] -391448.4984 -16553.3606
[53,] -47772.7110 -391448.4984
[54,] -407078.5585 -47772.7110
[55,] 152861.1425 -407078.5585
[56,] 151174.2325 152861.1425
[57,] -21731.3504 151174.2325
[58,] 25629.2641 -21731.3504
[59,] -249681.1403 25629.2641
[60,] -158595.9187 -249681.1403
[61,] -103861.8336 -158595.9187
[62,] -147288.9587 -103861.8336
[63,] -121221.0244 -147288.9587
[64,] 88306.9426 -121221.0244
[65,] 90107.6762 88306.9426
[66,] -95650.0029 90107.6762
[67,] -89941.3646 -95650.0029
[68,] -91213.1364 -89941.3646
[69,] -147913.3273 -91213.1364
[70,] 10949.3791 -147913.3273
[71,] 57506.3400 10949.3791
[72,] -188093.9212 57506.3400
[73,] -62163.5410 -188093.9212
[74,] -152804.1801 -62163.5410
[75,] -30854.7275 -152804.1801
[76,] -70049.8374 -30854.7275
[77,] -37018.9676 -70049.8374
[78,] 83844.8336 -37018.9676
[79,] -47622.7565 83844.8336
[80,] 20131.2434 -47622.7565
[81,] -109659.4849 20131.2434
[82,] -73596.1474 -109659.4849
[83,] 73092.8848 -73596.1474
[84,] -68100.5949 73092.8848
[85,] -165579.6510 -68100.5949
[86,] -181.4139 -165579.6510
[87,] 78258.6304 -181.4139
[88,] -73327.3995 78258.6304
[89,] -10263.5817 -73327.3995
[90,] -30223.2831 -10263.5817
[91,] -92607.5185 -30223.2831
[92,] 104610.2651 -92607.5185
[93,] -350911.9747 104610.2651
[94,] 189394.8250 -350911.9747
[95,] -46129.8832 189394.8250
[96,] -50431.2966 -46129.8832
[97,] -7656.7252 -50431.2966
[98,] -163892.6059 -7656.7252
[99,] 24718.5248 -163892.6059
[100,] -78894.3948 24718.5248
[101,] -91097.9397 -78894.3948
[102,] -114331.1095 -91097.9397
[103,] -36644.2181 -114331.1095
[104,] 27870.6053 -36644.2181
[105,] -47860.2078 27870.6053
[106,] -199900.1793 -47860.2078
[107,] -25894.5879 -199900.1793
[108,] -265545.6115 -25894.5879
[109,] -206906.6017 -265545.6115
[110,] 85871.1995 -206906.6017
[111,] -123771.6851 85871.1995
[112,] -245146.4930 -123771.6851
[113,] -171532.3565 -245146.4930
[114,] -18703.3689 -171532.3565
[115,] 97843.3432 -18703.3689
[116,] -20602.5726 97843.3432
[117,] 37320.8522 -20602.5726
[118,] -198173.0011 37320.8522
[119,] 61162.1655 -198173.0011
[120,] -321611.0043 61162.1655
[121,] -2794.4570 -321611.0043
[122,] -28581.9495 -2794.4570
[123,] 61586.0824 -28581.9495
[124,] -105202.5728 61586.0824
[125,] -246889.1051 -105202.5728
[126,] 53482.9591 -246889.1051
[127,] 75085.1371 53482.9591
[128,] 33180.3946 75085.1371
[129,] 53183.0694 33180.3946
[130,] 106685.9316 53183.0694
[131,] 97823.1152 106685.9316
[132,] -5239.7914 97823.1152
[133,] 7784.7346 -5239.7914
[134,] 26435.2879 7784.7346
[135,] 177598.0709 26435.2879
[136,] 82157.2142 177598.0709
[137,] 67249.1672 82157.2142
[138,] 53482.9591 67249.1672
[139,] 74794.3078 53482.9591
[140,] 24033.0410 74794.3078
[141,] 53183.0694 24033.0410
[142,] 98446.4628 53183.0694
[143,] 40338.0453 98446.4628
[144,] 1024.4114 40338.0453
[145,] 21860.2380 1024.4114
[146,] 26435.2879 21860.2380
[147,] 176950.4688 26435.2879
[148,] 79836.7782 176950.4688
[149,] 66677.6384 79836.7782
[150,] 44210.8872 66677.6384
[151,] 64236.9466 44210.8872
[152,] 24946.6516 64236.9466
[153,] 30999.5221 24946.6516
[154,] 103581.6414 30999.5221
[155,] 40071.6901 103581.6414
[156,] -5239.7914 40071.6901
[157,] 6640.6904 -5239.7914
[158,] 33635.8992 6640.6904
[159,] 174997.5049 33635.8992
[160,] 79836.7782 174997.5049
[161,] 67249.1672 79836.7782
[162,] 53482.9591 67249.1672
[163,] 80570.6160 53482.9591
[164,] 25606.0709 80570.6160
[165,] 53183.0694 25606.0709
[166,] 103494.3525 53183.0694
[167,] 144000.5954 103494.3525
[168,] -5239.7914 144000.5954
[169,] 9812.7590 -5239.7914
[170,] 26435.2879 9812.7590
[171,] 171989.4257 26435.2879
[172,] 79836.7782 171989.4257
[173,] 64997.4719 79836.7782
[174,] 53482.9591 64997.4719
[175,] 74830.4635 53482.9591
[176,] 24946.6516 74830.4635
[177,] -26761.3976 24946.6516
[178,] 103494.3525 -26761.3976
[179,] 40071.6901 103494.3525
[180,] 8870.4918 40071.6901
[181,] -92379.1625 8870.4918
[182,] 52232.5495 -92379.1625
[183,] 34221.3870 52232.5495
[184,] 391.3591 34221.3870
[185,] -18969.3233 391.3591
[186,] -20035.0601 -18969.3233
[187,] 71459.9865 -20035.0601
[188,] -14812.3441 71459.9865
[189,] 53183.0694 -14812.3441
[190,] 108202.7436 53183.0694
[191,] 40071.6901 108202.7436
[192,] -3080.7867 40071.6901
[193,] 9812.7590 -3080.7867
[194,] 26159.6990 9812.7590
[195,] 96294.5998 26159.6990
[196,] -135572.3399 96294.5998
[197,] -11453.7379 -135572.3399
[198,] -25219.9460 -11453.7379
[199,] -82122.4322 -25219.9460
[200,] -53756.2535 -82122.4322
[201,] -25519.8357 -53756.2535
[202,] 24791.4474 -25519.8357
[203,] -10296.1637 24791.4474
[204,] -170578.7136 -10296.1637
[205,] -70637.5945 -170578.7136
[206,] -67514.7389 -70637.5945
[207,] 96294.5998 -67514.7389
[208,] 60339.1038 96294.5998
[209,] -11453.7379 60339.1038
[210,] 53482.9591 -11453.7379
[211,] 74794.3078 53482.9591
[212,] 40850.3406 74794.3078
[213,] 30976.8117 40850.3406
[214,] 84784.1493 30976.8117
[215,] 40071.6901 84784.1493
[216,] -1697.9686 40071.6901
[217,] 14093.9080 -1697.9686
[218,] 26435.2879 14093.9080
[219,] 180079.5502 26435.2879
[220,] 47071.4215 180079.5502
[221,] 67249.1672 47071.4215
[222,] 54389.8956 67249.1672
[223,] 74794.3078 54389.8956
[224,] 25256.9161 74794.3078
[225,] 57993.3973 25256.9161
[226,] 103494.3525 57993.3973
[227,] 40071.6901 103494.3525
[228,] -31464.8866 40071.6901
[229,] 9812.7590 -31464.8866
[230,] 5291.3828 9812.7590
[231,] 174997.5049 5291.3828
[232,] 79836.7782 174997.5049
[233,] 68678.2748 79836.7782
[234,] 53482.9591 68678.2748
[235,] 75053.1947 53482.9591
[236,] 24946.6516 75053.1947
[237,] 56168.9458 24946.6516
[238,] 103494.3525 56168.9458
[239,] 40071.6901 103494.3525
[240,] -5239.7914 40071.6901
[241,] 12699.7489 -5239.7914
[242,] 26435.2879 12699.7489
[243,] 180228.4215 26435.2879
[244,] -7166.6580 180228.4215
[245,] -11453.7379 -7166.6580
[246,] -25219.9460 -11453.7379
[247,] 66680.0605 -25219.9460
[248,] -53756.2535 66680.0605
[249,] -25519.8357 -53756.2535
[250,] 110695.2132 -25519.8357
[251,] -38631.2150 110695.2132
[252,] -28672.4846 -38631.2150
[253,] -9733.9543 -28672.4846
[254,] -24628.0072 -9733.9543
[255,] 174997.5049 -24628.0072
[256,] 103324.4741 174997.5049
[257,] 107375.5412 103324.4741
[258,] -25219.9460 107375.5412
[259,] 74794.3078 -25219.9460
[260,] -53756.2535 74794.3078
[261,] -25519.8357 -53756.2535
[262,] 83708.2143 -25519.8357
[263,] 48712.1463 83708.2143
[264,] -83942.6965 48712.1463
[265,] 9812.7590 -83942.6965
[266,] -52267.6172 9812.7590
[267,] 178111.3539 -52267.6172
[268,] 1133.8731 178111.3539
[269,] 45606.4224 1133.8731
[270,] -25219.9460 45606.4224
[271,] -3908.5973 -25219.9460
[272,] -53756.2535 -3908.5973
[273,] 51734.8475 -53756.2535
[274,] 103494.3525 51734.8475
[275,] 49275.3029 103494.3525
[276,] 1841.8182 49275.3029
[277,] 9812.7590 1841.8182
[278,] -63545.2223 9812.7590
[279,] 153959.8020 -63545.2223
[280,] 1133.8731 153959.8020
[281,] 67249.1672 1133.8731
[282,] 92330.4772 67249.1672
[283,] 74794.3078 92330.4772
[284,] -11005.1296 74794.3078
[285,] 53183.0694 -11005.1296
[286,] 24791.4474 53183.0694
[287,] 39090.5666 24791.4474
[288,] -83942.6965 39090.5666
[289,] 40682.4378 -83942.6965
[290,] -52267.6172 40682.4378
[291,] 184802.0237 -52267.6172
[292,] 52906.5441 184802.0237
[293,] 67249.1672 52906.5441
[294,] 57166.2433 67249.1672
[295,] 30865.1386 57166.2433
[296,] 24946.6516 30865.1386
[297,] -25519.8357 24946.6516
[298,] 103494.3525 -25519.8357
[299,] 40071.6901 103494.3525
[300,] -84653.1602 40071.6901
[301,] -4192.3673 -84653.1602
[302,] 26435.2879 -4192.3673
[303,] 96294.5998 26435.2879
[304,] 118963.7940 96294.5998
[305,] 71370.1413 118963.7940
[306,] -25219.9460 71370.1413
[307,] -11902.9183 -25219.9460
[308,] 46885.7458 -11902.9183
[309,] -25519.8357 46885.7458
[310,] 24791.4474 -25519.8357
[311,] 19998.4788 24791.4474
[312,] -83942.6965 19998.4788
[313,] -71893.3472 -83942.6965
[314,] 28388.4039 -71893.3472
[315,] 176121.2446 28388.4039
[316,] 81440.5043 176121.2446
[317,] 53205.5005 81440.5043
[318,] 56254.6782 53205.5005
[319,] 74794.3078 56254.6782
[320,] 24946.6516 74794.3078
[321,] 53183.0694 24946.6516
[322,] 103494.3525 53183.0694
[323,] -38631.2150 103494.3525
[324,] -5239.7914 -38631.2150
[325,] 7249.0362 -5239.7914
[326,] -52267.6172 7249.0362
[327,] 96294.5998 -52267.6172
[328,] 79836.7782 96294.5998
[329,] 67249.1672 79836.7782
[330,] 53482.9591 67249.1672
[331,] 74794.3078 53482.9591
[332,] 24946.6516 74794.3078
[333,] 53227.2049 24946.6516
[334,] 104382.9762 53227.2049
[335,] 40071.6901 104382.9762
[336,] -6380.0132 40071.6901
[337,] 9812.7590 -6380.0132
[338,] 25787.9362 9812.7590
[339,] 179062.3190 25787.9362
[340,] 60789.2793 179062.3190
[341,] 67249.1672 60789.2793
[342,] 53482.9591 67249.1672
[343,] 74794.3078 53482.9591
[344,] 24946.6516 74794.3078
[345,] 53183.0694 24946.6516
[346,] 95507.1026 53183.0694
[347,] 40071.6901 95507.1026
[348,] -5239.7914 40071.6901
[349,] 9812.7590 -5239.7914
[350,] 27233.4321 9812.7590
[351,] 176349.8594 27233.4321
[352,] 81357.5962 176349.8594
[353,] 57358.7320 81357.5962
[354,] 53482.9591 57358.7320
[355,] 74794.3078 53482.9591
[356,] -8068.1788 74794.3078
[357,] 53183.0694 -8068.1788
[358,] 103494.3525 53183.0694
[359,] 40071.6901 103494.3525
[360,] -3409.1070 40071.6901
[361,] 9812.7590 -3409.1070
[362,] 26435.2879 9812.7590
[363,] 87151.6133 26435.2879
[364,] 78437.3973 87151.6133
[365,] 67249.1672 78437.3973
[366,] 53482.9591 67249.1672
[367,] 74794.3078 53482.9591
[368,] 24946.6516 74794.3078
[369,] 53183.0694 24946.6516
[370,] 103494.3525 53183.0694
[371,] 43205.1677 103494.3525
[372,] -45287.8150 43205.1677
[373,] -17023.1039 -45287.8150
[374,] 26435.2879 -17023.1039
[375,] 174997.5049 26435.2879
[376,] 92267.5274 174997.5049
[377,] 52111.7090 92267.5274
[378,] 53482.9591 52111.7090
[379,] 74794.3078 53482.9591
[380,] 24946.6516 74794.3078
[381,] 46373.5097 24946.6516
[382,] 157024.3228 46373.5097
[383,] 40071.6901 157024.3228
[384,] -3100.4391 40071.6901
[385,] -146528.7581 -3100.4391
[386,] 31132.1234 -146528.7581
[387,] 159162.8844 31132.1234
[388,] 78031.4273 159162.8844
[389,] 47984.8525 78031.4273
[390,] 50297.8980 47984.8525
[391,] 54309.1973 50297.8980
[392,] -12773.6816 54309.1973
[393,] 69737.5795 -12773.6816
[394,] -186404.4461 69737.5795
[395,] 12487.0534 -186404.4461
[396,] -169009.4846 12487.0534
[397,] -57017.7005 -169009.4846
[398,] -126.2104 -57017.7005
[399,] 90226.5753 -126.2104
[400,] 41927.2186 90226.5753
[401,] 57105.4751 41927.2186
[402,] 32843.7778 57105.4751
[403,] -86027.3038 32843.7778
[404,] -8671.2602 -86027.3038
[405,] 1044.3173 -8671.2602
[406,] 87763.6077 1044.3173
[407,] 19571.2178 87763.6077
[408,] -27565.2565 19571.2178
[409,] -65577.8594 -27565.2565
[410,] -18282.3546 -65577.8594
[411,] -107485.7438 -18282.3546
[412,] -28303.5090 -107485.7438
[413,] -151887.3283 -28303.5090
[414,] -57958.1541 -151887.3283
[415,] -19874.0729 -57958.1541
[416,] 107460.5380 -19874.0729
[417,] -79034.8951 107460.5380
[418,] 52454.3350 -79034.8951
[419,] -28508.0185 52454.3350
[420,] -127109.7247 -28508.0185
[421,] -146251.7195 -127109.7247
[422,] -462135.4596 -146251.7195
[423,] -350917.6996 -462135.4596
[424,] -47518.7003 -350917.6996
[425,] -197525.0426 -47518.7003
[426,] -308661.5966 -197525.0426
[427,] -1183573.5658 -308661.5966
[428,] -389711.3115 -1183573.5658
[429,] -271235.1059 -389711.3115
[430,] -793373.1387 -271235.1059
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 3216804.5650 1569858.1477
2 1318169.8718 3216804.5650
3 -3011398.4292 1318169.8718
4 -368870.5277 -3011398.4292
5 56132.3433 -368870.5277
6 2088.9969 56132.3433
7 616587.3876 2088.9969
8 -245850.0480 616587.3876
9 760920.0958 -245850.0480
10 -276567.2123 760920.0958
11 462075.4986 -276567.2123
12 905234.7776 462075.4986
13 -1819604.9030 905234.7776
14 16860.5273 -1819604.9030
15 558579.8501 16860.5273
16 -391037.7520 558579.8501
17 -63811.2844 -391037.7520
18 -152457.7990 -63811.2844
19 -309261.7393 -152457.7990
20 135480.0583 -309261.7393
21 -253842.5020 135480.0583
22 -190611.9735 -253842.5020
23 -504512.5871 -190611.9735
24 -277146.2879 -504512.5871
25 -171676.1145 -277146.2879
26 -400257.9907 -171676.1145
27 -321172.5938 -400257.9907
28 429960.1254 -321172.5938
29 36650.9527 429960.1254
30 169625.5040 36650.9527
31 -151503.2213 169625.5040
32 93830.9434 -151503.2213
33 5325.1372 93830.9434
34 -294870.5359 5325.1372
35 -469237.4376 -294870.5359
36 -5701.6948 -469237.4376
37 18715.6838 -5701.6948
38 -163183.1173 18715.6838
39 -30951.2419 -163183.1173
40 43415.1163 -30951.2419
41 -342227.2146 43415.1163
42 99832.9212 -342227.2146
43 -168130.2910 99832.9212
44 -43157.5823 -168130.2910
45 -243920.4219 -43157.5823
46 -211537.0486 -243920.4219
47 -197604.9207 -211537.0486
48 -166326.6387 -197604.9207
49 -203553.7208 -166326.6387
50 -13020.2667 -203553.7208
51 -16553.3606 -13020.2667
52 -391448.4984 -16553.3606
53 -47772.7110 -391448.4984
54 -407078.5585 -47772.7110
55 152861.1425 -407078.5585
56 151174.2325 152861.1425
57 -21731.3504 151174.2325
58 25629.2641 -21731.3504
59 -249681.1403 25629.2641
60 -158595.9187 -249681.1403
61 -103861.8336 -158595.9187
62 -147288.9587 -103861.8336
63 -121221.0244 -147288.9587
64 88306.9426 -121221.0244
65 90107.6762 88306.9426
66 -95650.0029 90107.6762
67 -89941.3646 -95650.0029
68 -91213.1364 -89941.3646
69 -147913.3273 -91213.1364
70 10949.3791 -147913.3273
71 57506.3400 10949.3791
72 -188093.9212 57506.3400
73 -62163.5410 -188093.9212
74 -152804.1801 -62163.5410
75 -30854.7275 -152804.1801
76 -70049.8374 -30854.7275
77 -37018.9676 -70049.8374
78 83844.8336 -37018.9676
79 -47622.7565 83844.8336
80 20131.2434 -47622.7565
81 -109659.4849 20131.2434
82 -73596.1474 -109659.4849
83 73092.8848 -73596.1474
84 -68100.5949 73092.8848
85 -165579.6510 -68100.5949
86 -181.4139 -165579.6510
87 78258.6304 -181.4139
88 -73327.3995 78258.6304
89 -10263.5817 -73327.3995
90 -30223.2831 -10263.5817
91 -92607.5185 -30223.2831
92 104610.2651 -92607.5185
93 -350911.9747 104610.2651
94 189394.8250 -350911.9747
95 -46129.8832 189394.8250
96 -50431.2966 -46129.8832
97 -7656.7252 -50431.2966
98 -163892.6059 -7656.7252
99 24718.5248 -163892.6059
100 -78894.3948 24718.5248
101 -91097.9397 -78894.3948
102 -114331.1095 -91097.9397
103 -36644.2181 -114331.1095
104 27870.6053 -36644.2181
105 -47860.2078 27870.6053
106 -199900.1793 -47860.2078
107 -25894.5879 -199900.1793
108 -265545.6115 -25894.5879
109 -206906.6017 -265545.6115
110 85871.1995 -206906.6017
111 -123771.6851 85871.1995
112 -245146.4930 -123771.6851
113 -171532.3565 -245146.4930
114 -18703.3689 -171532.3565
115 97843.3432 -18703.3689
116 -20602.5726 97843.3432
117 37320.8522 -20602.5726
118 -198173.0011 37320.8522
119 61162.1655 -198173.0011
120 -321611.0043 61162.1655
121 -2794.4570 -321611.0043
122 -28581.9495 -2794.4570
123 61586.0824 -28581.9495
124 -105202.5728 61586.0824
125 -246889.1051 -105202.5728
126 53482.9591 -246889.1051
127 75085.1371 53482.9591
128 33180.3946 75085.1371
129 53183.0694 33180.3946
130 106685.9316 53183.0694
131 97823.1152 106685.9316
132 -5239.7914 97823.1152
133 7784.7346 -5239.7914
134 26435.2879 7784.7346
135 177598.0709 26435.2879
136 82157.2142 177598.0709
137 67249.1672 82157.2142
138 53482.9591 67249.1672
139 74794.3078 53482.9591
140 24033.0410 74794.3078
141 53183.0694 24033.0410
142 98446.4628 53183.0694
143 40338.0453 98446.4628
144 1024.4114 40338.0453
145 21860.2380 1024.4114
146 26435.2879 21860.2380
147 176950.4688 26435.2879
148 79836.7782 176950.4688
149 66677.6384 79836.7782
150 44210.8872 66677.6384
151 64236.9466 44210.8872
152 24946.6516 64236.9466
153 30999.5221 24946.6516
154 103581.6414 30999.5221
155 40071.6901 103581.6414
156 -5239.7914 40071.6901
157 6640.6904 -5239.7914
158 33635.8992 6640.6904
159 174997.5049 33635.8992
160 79836.7782 174997.5049
161 67249.1672 79836.7782
162 53482.9591 67249.1672
163 80570.6160 53482.9591
164 25606.0709 80570.6160
165 53183.0694 25606.0709
166 103494.3525 53183.0694
167 144000.5954 103494.3525
168 -5239.7914 144000.5954
169 9812.7590 -5239.7914
170 26435.2879 9812.7590
171 171989.4257 26435.2879
172 79836.7782 171989.4257
173 64997.4719 79836.7782
174 53482.9591 64997.4719
175 74830.4635 53482.9591
176 24946.6516 74830.4635
177 -26761.3976 24946.6516
178 103494.3525 -26761.3976
179 40071.6901 103494.3525
180 8870.4918 40071.6901
181 -92379.1625 8870.4918
182 52232.5495 -92379.1625
183 34221.3870 52232.5495
184 391.3591 34221.3870
185 -18969.3233 391.3591
186 -20035.0601 -18969.3233
187 71459.9865 -20035.0601
188 -14812.3441 71459.9865
189 53183.0694 -14812.3441
190 108202.7436 53183.0694
191 40071.6901 108202.7436
192 -3080.7867 40071.6901
193 9812.7590 -3080.7867
194 26159.6990 9812.7590
195 96294.5998 26159.6990
196 -135572.3399 96294.5998
197 -11453.7379 -135572.3399
198 -25219.9460 -11453.7379
199 -82122.4322 -25219.9460
200 -53756.2535 -82122.4322
201 -25519.8357 -53756.2535
202 24791.4474 -25519.8357
203 -10296.1637 24791.4474
204 -170578.7136 -10296.1637
205 -70637.5945 -170578.7136
206 -67514.7389 -70637.5945
207 96294.5998 -67514.7389
208 60339.1038 96294.5998
209 -11453.7379 60339.1038
210 53482.9591 -11453.7379
211 74794.3078 53482.9591
212 40850.3406 74794.3078
213 30976.8117 40850.3406
214 84784.1493 30976.8117
215 40071.6901 84784.1493
216 -1697.9686 40071.6901
217 14093.9080 -1697.9686
218 26435.2879 14093.9080
219 180079.5502 26435.2879
220 47071.4215 180079.5502
221 67249.1672 47071.4215
222 54389.8956 67249.1672
223 74794.3078 54389.8956
224 25256.9161 74794.3078
225 57993.3973 25256.9161
226 103494.3525 57993.3973
227 40071.6901 103494.3525
228 -31464.8866 40071.6901
229 9812.7590 -31464.8866
230 5291.3828 9812.7590
231 174997.5049 5291.3828
232 79836.7782 174997.5049
233 68678.2748 79836.7782
234 53482.9591 68678.2748
235 75053.1947 53482.9591
236 24946.6516 75053.1947
237 56168.9458 24946.6516
238 103494.3525 56168.9458
239 40071.6901 103494.3525
240 -5239.7914 40071.6901
241 12699.7489 -5239.7914
242 26435.2879 12699.7489
243 180228.4215 26435.2879
244 -7166.6580 180228.4215
245 -11453.7379 -7166.6580
246 -25219.9460 -11453.7379
247 66680.0605 -25219.9460
248 -53756.2535 66680.0605
249 -25519.8357 -53756.2535
250 110695.2132 -25519.8357
251 -38631.2150 110695.2132
252 -28672.4846 -38631.2150
253 -9733.9543 -28672.4846
254 -24628.0072 -9733.9543
255 174997.5049 -24628.0072
256 103324.4741 174997.5049
257 107375.5412 103324.4741
258 -25219.9460 107375.5412
259 74794.3078 -25219.9460
260 -53756.2535 74794.3078
261 -25519.8357 -53756.2535
262 83708.2143 -25519.8357
263 48712.1463 83708.2143
264 -83942.6965 48712.1463
265 9812.7590 -83942.6965
266 -52267.6172 9812.7590
267 178111.3539 -52267.6172
268 1133.8731 178111.3539
269 45606.4224 1133.8731
270 -25219.9460 45606.4224
271 -3908.5973 -25219.9460
272 -53756.2535 -3908.5973
273 51734.8475 -53756.2535
274 103494.3525 51734.8475
275 49275.3029 103494.3525
276 1841.8182 49275.3029
277 9812.7590 1841.8182
278 -63545.2223 9812.7590
279 153959.8020 -63545.2223
280 1133.8731 153959.8020
281 67249.1672 1133.8731
282 92330.4772 67249.1672
283 74794.3078 92330.4772
284 -11005.1296 74794.3078
285 53183.0694 -11005.1296
286 24791.4474 53183.0694
287 39090.5666 24791.4474
288 -83942.6965 39090.5666
289 40682.4378 -83942.6965
290 -52267.6172 40682.4378
291 184802.0237 -52267.6172
292 52906.5441 184802.0237
293 67249.1672 52906.5441
294 57166.2433 67249.1672
295 30865.1386 57166.2433
296 24946.6516 30865.1386
297 -25519.8357 24946.6516
298 103494.3525 -25519.8357
299 40071.6901 103494.3525
300 -84653.1602 40071.6901
301 -4192.3673 -84653.1602
302 26435.2879 -4192.3673
303 96294.5998 26435.2879
304 118963.7940 96294.5998
305 71370.1413 118963.7940
306 -25219.9460 71370.1413
307 -11902.9183 -25219.9460
308 46885.7458 -11902.9183
309 -25519.8357 46885.7458
310 24791.4474 -25519.8357
311 19998.4788 24791.4474
312 -83942.6965 19998.4788
313 -71893.3472 -83942.6965
314 28388.4039 -71893.3472
315 176121.2446 28388.4039
316 81440.5043 176121.2446
317 53205.5005 81440.5043
318 56254.6782 53205.5005
319 74794.3078 56254.6782
320 24946.6516 74794.3078
321 53183.0694 24946.6516
322 103494.3525 53183.0694
323 -38631.2150 103494.3525
324 -5239.7914 -38631.2150
325 7249.0362 -5239.7914
326 -52267.6172 7249.0362
327 96294.5998 -52267.6172
328 79836.7782 96294.5998
329 67249.1672 79836.7782
330 53482.9591 67249.1672
331 74794.3078 53482.9591
332 24946.6516 74794.3078
333 53227.2049 24946.6516
334 104382.9762 53227.2049
335 40071.6901 104382.9762
336 -6380.0132 40071.6901
337 9812.7590 -6380.0132
338 25787.9362 9812.7590
339 179062.3190 25787.9362
340 60789.2793 179062.3190
341 67249.1672 60789.2793
342 53482.9591 67249.1672
343 74794.3078 53482.9591
344 24946.6516 74794.3078
345 53183.0694 24946.6516
346 95507.1026 53183.0694
347 40071.6901 95507.1026
348 -5239.7914 40071.6901
349 9812.7590 -5239.7914
350 27233.4321 9812.7590
351 176349.8594 27233.4321
352 81357.5962 176349.8594
353 57358.7320 81357.5962
354 53482.9591 57358.7320
355 74794.3078 53482.9591
356 -8068.1788 74794.3078
357 53183.0694 -8068.1788
358 103494.3525 53183.0694
359 40071.6901 103494.3525
360 -3409.1070 40071.6901
361 9812.7590 -3409.1070
362 26435.2879 9812.7590
363 87151.6133 26435.2879
364 78437.3973 87151.6133
365 67249.1672 78437.3973
366 53482.9591 67249.1672
367 74794.3078 53482.9591
368 24946.6516 74794.3078
369 53183.0694 24946.6516
370 103494.3525 53183.0694
371 43205.1677 103494.3525
372 -45287.8150 43205.1677
373 -17023.1039 -45287.8150
374 26435.2879 -17023.1039
375 174997.5049 26435.2879
376 92267.5274 174997.5049
377 52111.7090 92267.5274
378 53482.9591 52111.7090
379 74794.3078 53482.9591
380 24946.6516 74794.3078
381 46373.5097 24946.6516
382 157024.3228 46373.5097
383 40071.6901 157024.3228
384 -3100.4391 40071.6901
385 -146528.7581 -3100.4391
386 31132.1234 -146528.7581
387 159162.8844 31132.1234
388 78031.4273 159162.8844
389 47984.8525 78031.4273
390 50297.8980 47984.8525
391 54309.1973 50297.8980
392 -12773.6816 54309.1973
393 69737.5795 -12773.6816
394 -186404.4461 69737.5795
395 12487.0534 -186404.4461
396 -169009.4846 12487.0534
397 -57017.7005 -169009.4846
398 -126.2104 -57017.7005
399 90226.5753 -126.2104
400 41927.2186 90226.5753
401 57105.4751 41927.2186
402 32843.7778 57105.4751
403 -86027.3038 32843.7778
404 -8671.2602 -86027.3038
405 1044.3173 -8671.2602
406 87763.6077 1044.3173
407 19571.2178 87763.6077
408 -27565.2565 19571.2178
409 -65577.8594 -27565.2565
410 -18282.3546 -65577.8594
411 -107485.7438 -18282.3546
412 -28303.5090 -107485.7438
413 -151887.3283 -28303.5090
414 -57958.1541 -151887.3283
415 -19874.0729 -57958.1541
416 107460.5380 -19874.0729
417 -79034.8951 107460.5380
418 52454.3350 -79034.8951
419 -28508.0185 52454.3350
420 -127109.7247 -28508.0185
421 -146251.7195 -127109.7247
422 -462135.4596 -146251.7195
423 -350917.6996 -462135.4596
424 -47518.7003 -350917.6996
425 -197525.0426 -47518.7003
426 -308661.5966 -197525.0426
427 -1183573.5658 -308661.5966
428 -389711.3115 -1183573.5658
429 -271235.1059 -389711.3115
430 -793373.1387 -271235.1059
> 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/74q0e1291320718.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/84q0e1291320718.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/94q0e1291320718.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/10xhzz1291320718.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/110hy51291320718.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/124iwb1291320718.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/13sjb41291320718.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/143at71291320718.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/15ptrv1291320718.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/163lpm1291320718.tab")
+ }
> try(system("convert tmp/18y2n1291320718.ps tmp/18y2n1291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/28y2n1291320718.ps tmp/28y2n1291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/38y2n1291320718.ps tmp/38y2n1291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/4jp281291320718.ps tmp/4jp281291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jp281291320718.ps tmp/5jp281291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/6cg1b1291320718.ps tmp/6cg1b1291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/74q0e1291320718.ps tmp/74q0e1291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/84q0e1291320718.ps tmp/84q0e1291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/94q0e1291320718.ps tmp/94q0e1291320718.png",intern=TRUE))
character(0)
> try(system("convert tmp/10xhzz1291320718.ps tmp/10xhzz1291320718.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
15.944 3.122 16.442